The first complete T2T Assemblies of Cattle and Sheep Y-Chromosomes uncover remarkable divergence in structure and gene content.

Reference genomes of cattle and sheep have lacked contiguous assemblies of the sex-determining Y chromosome. We assembled complete and gapless telomere to telomere (T2T) Y chromosomes for these species. The pseudo-autosomal regions were similar in length, but the total chromosome size was substantially different, with the cattle Y more than twice the length of the sheep Y. The length disparity was accounted for by expanded ampliconic region in cattle. The genic amplification in cattle contrasts with pseudogenization in sheep suggesting opposite evolutionary mechanisms since their divergence 18MYA. The centromeres also differed dramatically despite the close relationship between these species at the overall genome sequence level. These Y chromosome have been added to the current reference assemblies in GenBank opening new opportunities for the study of evolution and variation while supporting efforts to improve sustainability in these important livestock species that generally use sire-driven genetic improvement strategies.

Genome-wide association with footrot in hair and wool sheep.

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum. Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH. We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1. These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete.

Analysis of Polycerate Mutants Reveals the Evolutionary Co-option of HOXD1 for Horn Patterning in Bovidae.

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as "headgear," which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

Genome-wide association reveals the locus responsible for four-horned ruminant.

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo-Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome-wide association mapping was performed using 125 animals from the two breeds that contain two- and four-horned individuals. A case-control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9-132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two-horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four-horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.

Ectopic Expression of Retrotransposon-Derived PEG11/RTL1 Contributes to the Callipyge Muscular Hypertrophy.

The callipyge phenotype is an ovine muscular hypertrophy characterized by polar overdominance: only heterozygous +Mat/CLPGPat animals receiving the CLPG mutation from their father express the phenotype. +Mat/CLPGPat animals are characterized by postnatal, ectopic expression of Delta-like 1 homologue (DLK1) and Paternally expressed gene 11/Retrotransposon-like 1 (PEG11/RTL1) proteins in skeletal muscle. We showed previously in transgenic mice that ectopic expression of DLK1 alone induces a muscular hypertrophy, hence demonstrating a role for DLK1 in determining the callipyge hypertrophy. We herein describe newly generated transgenic mice that ectopically express PEG11 in skeletal muscle, and show that they also exhibit a muscular hypertrophy phenotype. Our data suggest that both DLK1 and PEG11 act together in causing the muscular hypertrophy of callipyge sheep.

Experimental evaluation does not reveal a direct effect of microRNA from the callipyge locus on DLK1 expression.

BACKGROUND: Polar overdominance at the ovine callipyge (CLPG) locus involves the post-transcriptional trans-inhibition of DLK1 in skeletal muscle of CLPG/CLPG sheep. The abundant maternally expressed microRNAs (miRNAs) mapping to the imprinted DLK1-GTL2 domain are prime candidate mediators of this trans-effect. RESULTS: We have tested the affinity of 121 miRNAs processed from this locus for DLK1 by co-transfecting COS1 cells with a vector expressing the full-length ovine DLK1 with corresponding mimic miRNAs. None of the tested miRNAs was able to down regulate DLK1 to the extent observed in vivo. CONCLUSIONS: This suggests that other factors, with or without these miRNAs, are involved in mediating the observed trans-effect.

Impacts of the Callipyge mutation on ovine plasma metabolites and muscle fibre type.

The ovine Callipyge mutation causes postnatal muscle hypertrophy localized to the pelvic limbs and torso, as well as body leanness. The mechanism underpinning enhanced muscle mass is unclear, as is the systemic impact of the mutation. Using muscle fibre typing immunohistochemistry, we confirmed muscle specific effects and demonstrated that affected muscles had greater prevalence and hypertrophy of type 2X fast twitch glycolytic fibres and decreased representation of types 1, 2C, 2A and/or 2AX fibres. To investigate potential systemic effects of the mutation, proton NMR spectra of plasma taken from lambs at 8 and 12 weeks of age were measured. Multivariate statistical analysis of plasma metabolite profiles demonstrated effects of development and genotype but not gender. Plasma from Callipyge lambs at 12 weeks of age, but not 8 weeks, was characterized by a metabolic profile consistent with contributions from the affected hypertrophic fast twitch glycolytic muscle fibres. Microarray analysis of the perirenal adipose tissue depot did not reveal a transcriptional effect of the mutation in this tissue. We conclude that there is an indirect systemic effect of the Callipyge mutation in skeletal muscle in the form of changes of blood metabolites, which may contribute to secondary phenotypes such as body leanness.

Tracking the emergence of a new breed using 49,034 SNP in sheep.

Domestic animals are unique in that they have been organised into managed populations called breeds. The strength of genetic divergence between breeds may vary dependent on the age of the breed, the scenario under which it emerged and the strength of reproductive isolation it has from other breeds. In this study, we investigated the Gulf Coast Native breed of sheep to determine if it contains lines of animals that are sufficiently divergent to be considered separate breeds. Allele sharing and principal component analysis (PCA) using nearly 50,000 SNP loci revealed a clear genetic division that corresponded with membership of either the Florida or Louisiana Native lines. Subsequent analysis aimed to determine if the strength of the divergence exceeded that found between recognised breed pairs. Genotypes from 14 breeds sampled from Europe and Asia were used to obtain estimates of pair-wise population divergence measured as F(ST). The divergence separating the Florida and Louisiana Native (F(ST) = 6.2%) was approximately 50% higher than the average divergence separating breeds developed within the same region of Europe (F(ST) = 4.2%). This strongly indicated that the two Gulf Coast Native lines are sufficiently different to be considered separate breeds. PCA using small SNP sets successfully distinguished between the Florida and Louisiana Native animals, suggesting that allele frequency differences have accumulated across the genome. This is consistent with a population history involving geographic separation and genetic drift. Suggestive evidence was detected for divergence at the poll locus on sheep chromosome 10; however drift at neutral markers has been the largest contributor to the genetic separation observed. These results document the emergence of populations that can be considered separate breeds, with practical consequences for bio-conservation priorities, animal registration and the establishment of separate breed societies.

Assessing the effect of the CLPG mutation on the microRNA catalog of skeletal muscle using high-throughput sequencing.

The callipyge phenotype is a monogenic muscular hypertrophy that is only expressed in heterozygous sheep receiving the CLPG mutation from their sire. The wild-type phenotype of CLPG/CLPG animals is thought to result from translational inhibition of paternally expressed DLK1 transcripts by maternally expressed miRNAs. To identify the miRNA responsible for this trans effect, we used high-throughput sequencing to exhaustively catalog miRNAs expressed in skeletal muscle of sheep of the four CLPG genotypes. We have identified 747 miRNA species of which 110 map to the DLK1-GTL2 or callipyge domain. We demonstrate that the latter are imprinted and preferentially expressed from the maternal allele. We show that the CLPG mutation affects their level of expression in cis (∼3.2-fold increase) as well as in trans (∼1.8-fold increase). In CLPG/CLPG animals, miRNAs from the DLK1-GTL2 domain account for ∼20% of miRNAs in skeletal muscle. We show that the CLPG genotype affects the levels of A-to-I editing of at least five pri-miRNAs of the DLK1-GTL2 domain, but that levels of editing of mature miRNAs are always minor. We present suggestive evidence that the miRNAs from the domain target the ORF of DLK1, thereby causing the trans inhibition underlying polar overdominance. We highlight the limitations of high-throughput sequencing for digital gene expression profiling as a result of biased and inconsistent amplification of specific miRNAs.

A gene network switch enhances the oxidative capacity of ovine skeletal muscle during late fetal development.

BACKGROUND: The developmental transition between the late fetus and a newborn animal is associated with profound changes in skeletal muscle function as it adapts to the new physiological demands of locomotion and postural support against gravity. The mechanisms underpinning this adaption process are unclear but are likely to be initiated by changes in hormone levels. We tested the hypothesis that this developmental transition is associated with large coordinated changes in the transcription of skeletal muscle genes. RESULTS: Using an ovine model, transcriptional profiling was performed on Longissimus dorsi skeletal muscle taken at three fetal developmental time points (80, 100 and 120 d of fetal development) and two postnatal time points, one approximately 3 days postpartum and a second at 3 months of age. The developmental time course was dominated by large changes in expression of 2,471 genes during the interval between late fetal development (120 d fetal development) and 1-3 days postpartum. Analysis of the functions of genes that were uniquely up-regulated in this interval showed strong enrichment for oxidative metabolism and the tricarboxylic acid cycle indicating enhanced mitochondrial activity. Histological examination of tissues from these developmental time points directly confirmed a marked increase in mitochondrial activity between the late fetal and early postnatal samples. The promoters of genes that were up-regulated during this fetal to neonatal transition were enriched for estrogen receptor 1 and estrogen related receptor alpha cis-regulatory motifs. The genes down-regulated during this interval highlighted de-emphasis of an array of functions including Wnt signaling, cell adhesion and differentiation. There were also changes in gene expression prior to this late fetal--postnatal transition and between the two postnatal time points. The former genes were enriched for functions involving the extracellular matrix and immune response while the latter principally involved functions associated with transcriptional regulation of metabolic processes. CONCLUSIONS: It is concluded that during late skeletal muscle development there are substantial and coordinated changes in the transcription of a large number of genes many of which are probably triggered by increased estrogen levels. These changes probably underpin the adaption of muscle to new physiological demands in the postnatal environment.

Cytogenetic anchoring of radiation hybrid and virtual maps of sheep chromosome X and comparison of X chromosomes in sheep, cattle, and human.

A comprehensive physical map was generated for Ovis aries chromosome X (OARX) based on a cytogenomics approach. DNA probes were prepared from bacterial artificial chromosome (BAC) clones from the CHORI-243 sheep library and were assigned to G-banded metaphase spreads via fluorescence in-situ hybridization (FISH). A total of 22 BACs gave a single hybridization signal to the X chromosome and were assigned out of 32 tested. The positioned BACs contained 16 genes and a microsatellite marker which represent new cytogenetically mapped loci in the sheep genome. The gene and microsatellite loci serve to anchor between the existing radiation hybrid (RH) and virtual sheep genome (VSG) maps to the cytogenetic OARX map, whilst the BACs themselves also serve as anchors between the VSG and the cytogenetic maps. An additional 17 links between the RH and cytogenetic maps are provided by BAC end sequence (BES) derived markers that have also been positioned on the RH map. Comparison of the map orders for the cytogenetic, RH, and virtual maps reveals that the orders for the cytogenetic and RH maps are most similar, with only one locus, represented by BAC CH243-330E18, mapping to relatively different positions. Several discrepancies, including an inverted segment are found when comparing both the cytogenetic and RH maps with the virtual map. These discrepancies highlight the value of using physical mapping methods to inform the process of future in silico map construction. A detailed comparative analysis of sheep, human, and cattle mapping data allowed the construction of a comparative map that confirms and expands the knowledge about evolutionary conservation and break points between the X chromosomes of the three mammalian species.

The callipyge mutation and other genes that affect muscle hypertrophy in sheep.

Genetic strategies to improve the profitability of sheep operations have generally focused on traits for reproduction. However, natural mutations exist in sheep that affect muscle growth and development, and the exploitation of these mutations in breeding strategies has the potential to significantly improve lamb-meat quality. The best-documented mutation for muscle development in sheep is callipyge (CLPG), which causes a postnatal muscle hypertrophy that is localized to the pelvic limbs and loin. Enhanced skeletal muscle growth is also observed in animals with the Carwell (or rib-eye muscling) mutation, and a double-muscling phenotype has been documented for animals of the Texel sheep breed. However, the actual mutations responsible for these muscular hypertrophy phenotypes in sheep have yet to be identified, and further characterization of the genetic basis for these phenotypes will provide insight into the biological control of muscle growth and body composition.

Ectopic expression of DLK1 protein in skeletal muscle of padumnal heterozygotes causes the callipyge phenotype.

The callipyge (CLPG) phenotype is an inherited skeletal muscle hypertrophy described in sheep. It is characterized by an unusual mode of inheritance ("polar overdominance") in which only heterozygous individuals having received the CLPG mutation from their father (+(MAT)/CLPG(PAT)) express the phenotype . +(MAT)/CLPG(PAT) individuals are born normal and develop the muscular hypertrophy at approximately 1 month of age. The CLPG mutation was identified as an A to G transition in a highly conserved dodecamer motif located between the imprinted DLK1 and GTL2 genes . This motif is thought to be part of a long-range control element (LRCE) because the CLPG mutation was shown, in postnatal skeletal muscle, to enhance the transcript levels of the DLK1, PEG11, GTL2, and MEG8 genes in cis without altering their imprinting status . As a result, the +(MAT)/CLPG(PAT) individuals have a unique expression profile thought to underlie the callipyge phenotype: an overexpression of the paternally expressed protein encoding DLK1 (Figure 1A) and PEG11 transcripts in the absence of an overexpression of the maternally expressed noncoding GTL2 and MEG8 transcripts . However, the way in which this distinct expression profile causes the callipyge muscular hypertrophy has remained unclear. Herein, we demonstrate that the callipyge phenotype is perfectly correlated with ectopic expression of DLK1 protein in hypertrophied muscle of +(MAT)/CLPG(PAT) sheep. We demonstrate the causality of this association by inducing a generalized muscular hypertrophy in transgenic mice that express DLK1 in skeletal muscle. The absence of DLK1 protein in skeletal muscle of CLPG/CLPG animals, despite the presence of DLK1 mRNA, supports a trans inhibition mediated by noncoding RNAs expressed from the maternal allele.

Expression of PEG11 and PEG11AS transcripts in normal and callipyge sheep.

BACKGROUND: The callipyge mutation is located within an imprinted gene cluster on ovine chromosome 18. The callipyge trait exhibits polar overdominant inheritance due to the fact that only heterozygotes inheriting a mutant paternal allele (paternal heterozygotes) have a phenotype of muscle hypertrophy, reduced fat and a more compact skeleton. The mutation is a single A to G transition in an intergenic region that results in the increased expression of several genes within the imprinted cluster without changing their parent-of-origin allele-specific expression. RESULTS: There was a significant effect of genotype (p < 0.0001) on the transcript abundance of DLK1, PEG11, and MEG8 in the muscles of lambs with the callipyge allele. DLK1 and PEG11 transcript levels were elevated in the hypertrophied muscles of paternal heterozygous animals relative to animals of the other three genotypes. The PEG11 locus produces a single 6.5 kb transcript and two smaller antisense strand transcripts, referred to as PEG11AS, in skeletal muscle. PEG11AS transcripts were detectable over a 5.5 kb region beginning 1.2 kb upstream of the PEG11 start codon and spanning the entire open reading frame. Analysis of PEG11 expression by quantitative PCR shows a 200-fold induction in the hypertrophied muscles of paternal heterozygous animals and a 13-fold induction in homozygous callipyge animals. PEG11 transcripts were 14-fold more abundant than PEG11AS transcripts in the gluteus medius of paternal heterozygous animals. PEG11AS transcripts were expressed at higher levels than PEG11 transcripts in the gluteus medius of animals of the other three genotypes. CONCLUSIONS: The effect of the callipyge mutation has been to alter the expression of DLK1, GTL2, PEG11 and MEG8 in the hypertrophied skeletal muscles. Transcript abundance of DLK1 and PEG11 was highest in paternal heterozygous animals and exhibited polar overdominant gene expression patterns; therefore, both genes are candidates for causing skeletal muscle hypertrophy. There was unique relationship of PEG11 and PEG11AS transcript abundance in the paternal heterozygous animals that suggests a RNA interference mechanism may have a role in PEG11 gene regulation and polar overdominance in callipyge sheep.