Juan José Arranz

Revealing the History of Sheep Domestication Using Retrovirus Integrations

Sheep retroviruses can be used to map the selective preferences of early farmers and trace livestock movements across Europe. Not Just Dinner on Legs Several thousand years ago, human beings realized the virtues of domesticating wild animals as easy meat. Soon other possibilities became apparent, and as revealed in a series of papers in this issue, early pastoralists became selective about breeding for wool, leather, milk, and muscle power. In two papers, Gibbs et al. report on the bovine genome sequence (p. 522; see the cover, the Perspective by Lewin, and the Policy Forum by Roberts) and trace the diversity and genetic history of cattle (p. 528), while Chessa et al. (p. 532) survey the occurrence of endogenous retroviruses in sheep and map their distribution to historical waves of human selection and dispersal across Europe. Finally, Ludwig et al. (p. 485) note the origins of variation in the coat-color of horses and suggest that it is most likely to have been selected for by humans in need of good-looking transport. The domestication of livestock represented a crucial step in human history. By using endogenous retroviruses as genetic markers, we found that sheep differentiated on the basis of their “retrotype” and morphological traits dispersed across Eurasia and Africa via separate migratory episodes. Relicts of the first migrations include the Mouflon, as well as breeds previously recognized as “primitive” on the basis of their morphology, such as the Orkney, Soay, and the Nordic short-tailed sheep now confined to the periphery of northwest Europe. A later migratory episode, involving sheep with improved production traits, shaped the great majority of present-day breeds. The ability to differentiate genetically primitive sheep from more modern breeds provides valuable insights into the history of sheep domestication.

Differences in the expression of the ASIP gene are involved in the recessive black coat colour pattern in sheep: evidence from the rare Xalda sheep breed

Here we have tested the hypothesis of association between different levels of agouti signalling peptide (ASIP) mRNA and the recessive black coat colour in the rare Xalda breed of sheep. To deal with this task, we first tested the possible action of both the dominant black extension allele (E(D)) and a 5-bp deletion (X99692:c.100_104del; A(del)) in the ovine ASIP coding sequence on the black coat colour pattern in 188 Xalda individuals. The E(D) allele was not present in the sample and only 11 individuals were homozygous for the A(del)ASIP allele. All Xalda individuals carrying the A(del)/A(del) genotype were phenotypically black. However, most black-coated individuals (109 out of 120) were not homozygous for the 5-bp deletion, thus rejecting the A(del)/A(del) genotype as the sole cause of recessive black coat colour in sheep. Differences in expression of ASIP mRNA were assessed via RT-PCR in 14 black-coated and 10 white-coated Xalda individuals showing different ASIP genotypes (A(wt)/A(wt), A(wt)/A(del) and A(del)/A(del)). Levels of expression in black animals were significantly (P < 0.0001) lower than those assessed for white-coated individuals. However, the ASIP genotype did not influence the ASIP mRNA level of expression. The consistency of these findings with those recently reported in humans is discussed, and the need to isolate the promoter region of ovine ASIP to obtain further evidence for a role of ASIP in recessive black ovine pigmentation is pointed out.

An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds

This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as “divergence signals,” which may be useful in characterizing and protecting livestock genetic diversity.

Application of selection mapping to identify genomic regions associated with dairy production in sheep

In Europe, especially in Mediterranean areas, the sheep has been traditionally exploited as a dual purpose species, with income from both meat and milk. Modernization of husbandry methods and the establishment of breeding schemes focused on milk production have led to the development of “dairy breeds.” This study investigated selective sweeps specifically related to dairy production in sheep by searching for regions commonly identified in different European dairy breeds. With this aim, genotypes from 44,545 SNP markers covering the sheep autosomes were analysed in both European dairy and non-dairy sheep breeds using two approaches: (i) identification of genomic regions showing extreme genetic differentiation between each dairy breed and a closely related non-dairy breed, and (ii) identification of regions with reduced variation (heterozygosity) in the dairy breeds using two methods. Regions detected in at least two breeds (breed pairs) by the two approaches (genetic differentiation and at least one of the heterozygosity-based analyses) were labeled as core candidate convergence regions and further investigated for candidate genes. Following this approach six regions were detected. For some of them, strong candidate genes have been proposed (e.g. ABCG2, SPP1), whereas some other genes designated as candidates based on their association with sheep and cattle dairy traits (e.g. LALBA, DGAT1A) were not associated with a detectable sweep signal. Few of the identified regions were coincident with QTL previously reported in sheep, although many of them corresponded to orthologous regions in cattle where QTL for dairy traits have been identified. Due to the limited number of QTL studies reported in sheep compared with cattle, the results illustrate the potential value of selection mapping to identify genomic regions associated with dairy traits in sheep.

Characterization and Comparative Analysis of the Milk Transcriptome in Two Dairy Sheep Breeds using RNA Sequencing.

This study presents a dynamic characterization of the sheep milk transcriptome aiming at achieving a better understanding of the sheep lactating mammary gland. Transcriptome sequencing (RNA-seq) was performed on total RNA extracted from milk somatic cells from ewes on days 10, 50, 120 and 150 after lambing. The experiment was performed in Spanish Churra and Assaf breeds, which differ in their milk production traits. Nearly 67% of the annotated genes in the reference genome (Oar_v3.1) were expressed in ovine milk somatic cells. For the two breeds and across the four lactation stages studied, the most highly expressed genes encoded caseins and whey proteins. We detected 573 differentially expressed genes (DEGs) across lactation points, with the largest differences being found, between day 10 and day 150. Upregulated GO terms at late lactation stages were linked mainly to developmental processes linked to extracellular matrix remodeling. A total of 256 annotated DEGs were detected in the Assaf and Churra comparison. Some genes selectively upregulated in the Churra breed grouped under the endopeptidase and channel activity GO terms. These genes could be related to the higher cheese yield of this breed. Overall, this study provides the first integrated overview on sheep milk gene expression.

Using Regulatory and Epistatic Networks to Extend the Findings of a Genome Scan: Identifying the Gene Drivers of Pigmentation in Merino Sheep

Extending genome wide association analysis by the inclusion of gene expression data may assist in the dissection of complex traits. We examined piebald, a pigmentation phenotype in both human and Merino sheep, by analysing multiple data types using a systems approach. First, a case control analysis of 49,034 ovine SNP was performed which confirmed a multigenic basis for the condition. We combined these results with gene expression data from five tissue types analysed with a skin-specific microarray. Promoter sequence analysis of differentially expressed genes allowed us to reverse-engineer a regulatory network. Likewise, by testing two-loci models derived from all pair-wise comparisons across piebald-associated SNP, we generated an epistatic network. At the intersection of both networks, we identified thirteen genes with insulin-like growth factor binding protein 7 (IGFBP7), platelet-derived growth factor alpha (PDGFRA) and the tetraspanin platelet activator CD9 at the kernel of the intersection. Further, we report a number of differentially expressed genes in regions containing highly associated SNP including ATRN, DOCK7, FGFR1OP, GLI3, SILV and TBX15. The application of network theory facilitated co-analysis of genetic variation with gene expression, recapitulated aspects of the known molecular biology of skin pigmentation and provided insights into the transcription regulation and epistatic interactions involved in piebald Merino sheep.

Identification of a 31-bp deletion in the RELN gene causing lissencephaly with cerebellar hypoplasia in sheep

Lissencephaly is an inherited developmental disorder in which neuronal migration is impaired. A type of lissencephaly associated with cerebellar hypoplasia (LCH) was diagnosed in a commercial flock of Spanish Churra sheep. The genotyping of 7 affected animals and 33 controls with the OvineSNP50 BeadChip enabled the localization of the causative mutation for ovine LCH to a 4.8-Mb interval on sheep chromosome 4 using genome-wide association and homozygosity mapping. The RELN gene, which is located within this interval, was considered a strong positional and functional candidate because it plays critical roles in neuronal migration and layer formation. By performing a sequencing analysis of this gene’s specific mRNA in a control lamb, we obtained the complete CDS of the ovine RELN gene. The cDNA sequence from an LCH-affected lamb revealed a deletion of 31 bp (c.5410_5440del) in predicted exon 36 of RELN, resulting in a premature termination codon. A functional analysis of this mutation revealed decreased levels of RELN mRNA and a lack of reelin protein in the brain cortex and blood of affected lambs. This mutation showed a complete concordance with the Mendelian recessive pattern of inheritance observed for the disease. The identification of the causal mutation of LCH in Churra sheep will facilitate the implementation of gene-assisted selection to detect heterozygous mutants, which will help breeders avoid at-risk matings in their flocks. Moreover, the identification of this naturally occurring RELN mutation provides an opportunity to use Churra sheep as a genetically characterized large animal model for the study of reelin functions in the developing and mature brain.

GENETIC STRUCTURE AND DIVERSITY OF FOUR CHILEAN SHEEP BREEDS

Chilean sheep breeds have a significant founder effect since they descend from very few parents, in some cases population inbreeding has increased and in others there have been significant differentiation processes as compared to their original population. The objective of this study was to estimate the current status of genetic diversity in sheep herds, which comprise the main sheep breeds in Chile, using molecular methods. Animals from four breeding herds were sampled and analyzed. The breeds studied were Corriedale, Suffolk Down, Romney Marsh and Chilota, these are the predominant numerical pure breeds in Chile. From each breed population 40 animals were sampled. Genetic characterization was done using nine microsatellite molecular markers (SSR) recommended by FAO-ISAG. Results showed that in the Chilean sheep herds there are low genetic complementarities among breeds and a high risk of losing genetic diversity due to inbreeding in Chilota and Romney Marsh breeds.

Elucidating fish oil-induced milk fat depression in dairy sheep: Milk somatic cell transcriptome analysis

In this study, RNA sequencing was used to obtain a comprehensive profile of the transcriptomic changes occurring in the mammary gland of lactating sheep suffering from fish oil-induced milk fat depression (FO-MFD). The milk somatic cell transcriptome analysis of four control and four FO-MFD ewes generated an average of 42 million paired-end reads per sample. In both conditions, less than 220 genes constitute approximately 89% of the total counts. These genes, which are considered as core genes, were mainly involved in cytoplasmic ribosomal proteins and electron transport chain pathways. In total, 117 genes were upregulated, and 96 genes were downregulated in FO-MFD samples. Functional analysis of the latter indicated a downregulation of genes involved in the SREBP signaling pathway (e.g., ACACA, ACSL, and ACSS) and Gene Ontology terms related to lipid metabolism and lipid biosynthetic processes. Integrated interpretation of upregulated genes indicated enrichment in genes encoding plasma membrane proteins and proteins regulating protein kinase activity. Overall, our results indicate that FO-MFD is associated with the downregulation of key genes involved in the mammary lipogenesis process. In addition, the results also suggest that this syndrome may be related to upregulation of other genes implicated in signal transduction and codification of transcription factors.

The β-amyloid peptide compromises Reelin signaling in Alzheimer’s disease

Reelin is a signaling protein that plays a crucial role in synaptic function, which expression is influenced by β-amyloid (Aβ). We show that Reelin and Aβ oligomers co-immunoprecipitated in human brain extracts and were present in the same size-exclusion chromatography fractions. Aβ treatment of cells led to increase expression of Reelin, but secreted Reelin results trapped together with Aβ aggregates. In frontal cortex extracts an increase in Reelin mRNA, and in soluble and insoluble (guanidine-extractable) Reelin protein, was associated with late Braak stages of Alzheimer’s disease (AD), while expression of its receptor, ApoER2, did not change. However, Reelin-dependent induction of Dab1 phosphorylation appeared reduced in AD. In cells, Aβ reduced the capacity of Reelin to induce internalization of biotinylated ApoER2 and ApoER2 processing. Soluble proteolytic fragments of ApoER2 generated after Reelin binding can be detected in cerebrospinal fluid (CSF). Quantification of these soluble fragments in CSF could be a tool to evaluate the efficiency of Reelin signaling in the brain. These CSF-ApoER2 fragments correlated with Reelin levels only in control subjects, not in AD, where these fragments diminished. We conclude that while Reelin expression is enhanced in the Alzheimer’s brain, the interaction of Reelin with Aβ hinders its biological activity.