Genome‐wide analysis of Italian sheep diversity reveals a strong geographic pattern and cryptic relationships between breeds

Italy counts several sheep breeds, arisen over centuries as a consequence of ancient and recent genetic and demographic events. To finely reconstruct genetic structure and relationships between Italian sheep, 496 subjects from 19 breeds were typed at 50K single nucleotide polymorphism loci. A subset of foreign breeds from the Sheep HapMap dataset was also included in the analyses. Genetic distances (as visualized either in a network or in a multidimensional scaling analysis of identical by state distances) closely reflected geographic proximity between breeds, with a clear north–south gradient, likely because of high levels of past gene flow and admixture all along the peninsula. Sardinian breeds diverged more from other breeds, a probable consequence of the combined effect of ancient sporadic introgression of feral mouflon and long‐lasting genetic isolation from continental sheep populations. The study allowed the detection of previously undocumented episodes of recent introgression (Delle Langhe into the endangered Altamurana breed) as well as signatures of known, or claimed, historical introgression (Merino into Sopravissana and Gentile di Puglia; Bergamasca into Fabrianese, Appenninica and, to a lesser extent, Leccese). Arguments that would question, from a genomic point of view, the current breed classification of Bergamasca and Biellese into two separate breeds are presented. Finally, a role for traditional transhumance practices in shaping the genetic makeup of Alpine sheep breeds is proposed. The study represents the first exhaustive analysis of Italian sheep diversity in an European context, and it bridges the gap in the previous HapMap panel between Western Mediterranean and Swiss breeds.     

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow – a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co‐occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment‐length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host‐associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.     

Splenic gene expression profiling in White Leghorn layer inoculated with the Salmonella enterica serovar Enteritidis

Salmonella enterica serovar Enteritidis (SE) is a foodborne pathogen that can threaten human health through contaminated poultry products. Live poultry, chicken eggs and meat are primary sources of human salmonellosis. To understand the genetic resistance of egg‐type chickens in response to SE inoculation, global gene expression in the spleen of 20‐week‐old White Leghorn was measured using the Agilent 4 × 44 K chicken microarray at 7 and 14 days following SE inoculation (dpi). Results showed that there were 1363 genes significantly differentially expressed between inoculated and non‐inoculated groups at 7 dpi (I7/N7), of which 682 were up‐regulated and 681 were down‐regulated genes. By contrast, 688 differentially expressed genes were observed at 14 dpi (I14/N14), of which 371 were up‐regulated genes and 317 were down‐regulated genes. There were 33 and 28 immune‐related genes significantly differentially expressed in the comparisons of I7/N7 and I14/N14 respectively. Functional annotation revealed that several Gene Ontology (GO) terms related to immunity were significantly enriched between the inoculated and non‐inoculated groups at 14 dpi but not at 7 dpi, despite a similar number of immune‐related genes identified between I7/N7 and I14/N14. The immune response to SE inoculation changes with different time points following SE inoculation. The complicated interaction between the immune system and metabolism contributes to the immune responses to SE inoculation of egg‐type chickens at 14 dpi at the onset of lay. GC, TNFSF8, CD86, CD274, BLB1 and BLB2 play important roles in response to SE inoculation. The results from this study will deepen the current understanding of the genetic response of the egg‐type chicken to SE inoculation at the onset of egg laying.     

DNA binding and bending by a chloroplast-encoded HU-like protein overexpressed in Escherichia coli

The Guillardia theta chloroplast hlpA gene encodes a protein resembling bacterial histone-like protein HU. This gene was cloned and overexpressed in Escherichia coli cells, and the resulting protein product, HlpA, was purified and characterized in vitro. In addition to exhibiting a general DNA-binding activity, the chloroplast HlpA protein also strongly facilitated cyclization of a short DNA fragment in the presence of T4 DNA ligase, indicating its ability to mediate very tight DNA curvatures.     

Comparative mitochondrial and chloroplast genomics of a genetically distinct form of Sargassum contributing to recent “Golden Tides” in the Western Atlantic

Over the past 5 years, massive accumulations of holopelagic species of the brown macroalga Sargassum in coastal areas of the Caribbean have created “golden tides” that threaten local biodiversity and trigger economic losses associated with beach deterioration and impact on fisheries and tourism. In 2015, the first report identifying the cause of these extreme events implicated a rare form of the holopelagic species Sargassum natans (form VIII). However, since the first mention of S. natans VIII in the 1930s, based solely on morphological characters, no molecular data have confirmed this identification. We generated full‐length mitogenomes and partial chloroplast genomes of all representative holopelagic Sargassum species, S. fluitans III and S. natans I alongside the putatively rare S. natans VIII, to demonstrate small but consistent differences between S. natans I and VIII (7 bp differences out of the 34,727). Our comparative analyses also revealed that both S. natans I and S. natans VIII share a very close phylogenetic relationship with S. fluitans III (94‐ and 96‐bp differences of 34,727). We designed novel primers that amplified regions of the cox2 and cox3 marker genes with consistent polymorphic sites that enabled differentiation between the two S. natans forms (I and VIII) from each other and both from S. fluitans III in over 150 Sargassum samples including those from the 2014 golden tide event. Despite remarkable gene synteny and sequence conservation, the three Sargassum forms differ in morphology, ecology, and distribution patterns, warranting more extensive interrogation of holopelagic Sargassum genomes as a whole.     

Direct characterization of bovine microsatellites from cosmids: polymorphism and synteny mapping

A (TG)₈ oligonucleotide probe was used to screen 186 cosmids from a commercial bovine cosmid library. Of the 56 positive discovered, 7 were sequenced in the region of the microsatellite and analysed for polymorphism. These microsatellites, IDVGA-2, -3, -7, -8, -9, -10, -11 showed the following number of alleles and polymorphism information content (PIC) values (7/0.616, 8/0.693, 6/0.641, 5/0.643, 2/0.239, 10/0.844, 6/0.720). The microsatellites were also assigned to synteny groups as follows: IDVGA-2/U17, IDVGA-3/U16, IDVGA-7/U7, IDVGA-8/U29, IDVGA-9/U3, IDVGA-10/U19, IDVGA-11/A (probably U18).