Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection.

BACKGROUND: Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding. RESULTS: CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences. CONCLUSIONS: Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep.

BACKGROUND: The detection of selective traits in different populations can not only reveal current mechanisms of artificial selection for breeding, but also provide new insights into phenotypic variation in new varieties and the search for genes associated with important traits. Panou sheep is a cultivated breed of Tibetan sheep in China with stable genetic performance, consistent appearance and fast growth and development after decades of artificial selection and cultivation. Due to long-term adaptation to the high altitude, cold and hypoxic environment in the plateau area, they may have formed a unique gene pool that is different from other Tibetan sheep breeds. To explore the genetic resources of Panou sheep, we used next-generation sequencing technology for the first time to investigate the genome-wide population structure, genetic diversity, and candidate signatures of positive selection in Panou sheep. RESULTS: Comparative genomic analysis with the closely related species Oula sheep (a native breed of Tibetan sheep in China) was used to screen the population selection signal of Panou sheep. Principal component analysis and neighbor joining tree showed that Panou sheep and Oula sheep had differences in population differentiation. Furthermore, analyses of population structure, they came from the same ancestor, and when K = 2, the two populations could be distinguished. Panou sheep exhibit genetic diversity comparable to Oula sheep, as shown by observed heterozygosity, expected heterozygosity and runs of homozygosity. Genome-wide scanning using the Fst and π ratio methods revealed a list of potentially selected related genes in Panou sheep compared to Oula sheep, including histone deacetylase 9 (HDAC9), protein tyrosine kinase 2 (PTK2), microphthalmia-related transcription factor (MITF), vesicular amine transporter 1 (VAT1), trichohyalin-like 1 (TCHHL1), amine oxidase, copper containing 3 (AOC3), interferon-inducible protein 35 (IFI35). CONCLUSIONS: The results suggest that traits related to growth and development and plateau adaptation may be selection targets for the domestication and breeding improvement of Tibetan sheep. This study provides the fundamental footprints for Panou sheep breeding and management.

Differential tissue expression of sex steroid-synthesizing enzyme CYP11A1 in male Tibetan sheep (Ovis aries).

Steroid metabolism is a fundament to testicular development and function. The cytochrome P450, family 11, subfamily A, polypeptide 1 (CYP11A1) is a key rate-limiting enzyme for catalyzing the conversion of cholesterol to pregnenolone. However, despite its importance, what expression and roles of CYP11A1 possesses and how it regulates the testicular development and spermatogenesis in Tibetan sheep remains largely unknown. Based on this, we evaluated the expression and localization patterns of CYP11A1 in testes and epididymides of Tibetan sheep at three developmental stages (three-month-old, pre-puberty; one-year-old, sexual maturity and three-year-old, adult) by quantitative real-time PCR (qPCR), western blot and immunofluorescence. The results showed that CYP11A1 mRNA and protein were expressed in testes and epididymides throughout the development stages and obviously more intense in one- and three-year-old groups than three-month-old group (except for the caput epididymidis). Immunofluorescence assay showed that the CYP11A1 protein was mainly located in Leydig cells and epididymal epithelial cells. In addition, positive signals of CYP11A1 protein were observed in germ cells, epididymal connective tissue and sperms stored in the epididymal lumen. Collectively, these results suggested that the CYP11A1 gene might be mainly involved in regulating spermatogenesis and androgen synthesis in developmental Tibetan sheep testis and epididymis.

Truncation of the Murine Neonatal Fc Receptor Cytoplasmic Tail Does Not Alter IgG Metabolism or Transport In Vivo.

The neonatal Fc receptor (FcRn) is involved in IgG metabolism and transport in placental mammals. However, whether FcRn is responsible for IgG transfer from maternal serum to colostrum/milk is controversial. Interestingly, large domestic animals, such as cows, pigs, sheep, and horses, in which passive IgG transfer is exclusively completed via colostrum/milk, all express an FcRn α-chain that is shorter in the cytoplasmic tail (CYT) than its counterparts in humans and rodents. To address whether the length variation has any functional significance, we performed in vitro experiments using the Transwell system with the MDCK cell line stably transfected with various FcRn constructs; these clearly suggested that truncation of the CYT tail caused a polar change in IgG transfer. However, we observed no evidence supporting functional changes in IgG in vivo using mice in which the FcRn CYT was precisely truncated. These data suggest that the length variation in FcRn is not functionally associated with passive IgG transfer routes in mammals.

Cloning, Molecular Characterization and Expression Patterns of DMRTC2 Implicated in Germ Cell Development of Male Tibetan Sheep.

The double sex and mab-3-related transcription factors like family C2 (DMRTC2) gene is indispensable for mammalian testicular function and spermatogenesis. Despite its importance, what expression and roles of DMRTC2 possesses and how it regulates the testicular development and spermatogenesis in sheep, especially in Tibetan sheep, remains largely unknown. In this study, DMRTC2 cDNA from testes of Tibetan sheep was firstly cloned by the RT-PCR method, and its molecular characterization was identified. Subsequently, the expression and localization patterns of DMRTC2 were evaluated by quantitative real-time PCR (qPCR), Western blot, and immunofluorescence. The cloning and sequence analysis showed that the Tibetan sheep DMRTC2 cDNA fragment contained 1113 bp open reading frame (ORF) capable of encoding 370 amino acids, and displayed high identities with some other mammals, which shared an identical DM domain sequence of 47 amino acids ranged from residues 38 to 84. qPCR and Western blot results showed that DMRTC2 was expressed in testes throughout the development stages while not in epididymides (caput, corpus, and cauda), with higher mRNA and protein abundance in Tibetan sheep testes of one- and three-year-old (post-puberty) compared with that of three-month-old (pre-puberty). Immunofluorescence results revealed that immune staining for DMRTC2 protein was observed in spermatids and spermatogonia from post-puberty Tibetan sheep testes, and gonocytes from pre-puberty Tibetan sheep testes. Together, these results demonstrated, for the first time, in sheep, that DMRTC2, as a highly conserved gene in mammals, is essential for sheep spermatogenesis by regulating the proliferation or differentiation of gonocytes and development of spermatids in ram testes at different stages of maturity.

Identification and verification of potential piRNAs from domesticated yak testis.

PIWI-interacting RNAs (piRNA) are small non-coding RNA molecules expressed in animal germ cells that interact with PIWI family proteins to form RNA-protein complexes involved in epigenetic and post-transcriptional gene silencing of retrotransposons and other genetic elements in germ line cells, including reproductive stem cell self-sustainment, differentiation, meiosis and spermatogenesis. In the present study, we performed high-throughput sequencing of piRNAs in testis samples from yaks in different stages of sexual maturity. Deep sequencing of the small RNAs (18-40 nt in length) yielded 4,900,538 unique reads from a total of 53,035,635 reads. We identified yak small RNAs (18-30 nt) and performed functional characterization. Yak small RNAs showed a bimodal length distribution, with two peaks at 22 nt and >28 nt. More than 80% of the 3,106,033 putative piRNAs were mapped to 4637 piRNA-producing genomic clusters using RPKM. 6388 candidate piRNAs were identified from clean reads and the annotations were compared with the yak reference genome repeat region. Integrated network analysis suggested that some differentially expressed genes were involved in spermatogenesis through ECM-receptor interaction and PI3K-Akt signaling pathways. Our data provide novel insights into the molecular expression and regulation similarities and diversities in spermatogenesis and testicular development in yaks at different stages of sexual maturity.

Expression of melatonin and its related synthase and membrane receptors in the oestrous corpus luteum and corpus luteum verum of sheep.

Melatonin is an important factor involved in regulating reproduction; it is synthesized enzymatically by the sequential action of melatonin-synthesizing enzymes, arylalkylamine N-acetyltransferase (AANAT) and hydroxyindole-O-methyltransferase (HIOMT), and exerts its biological functions mainly through receptor-mediated action. To evaluate the expression of melatonin, two melatonin-synthesizing enzymes (HIOMT and AANAT), and membrane receptors (MT1 and MT2) in oestrous corpus luteum (CL) and CL verum of sheep (Ovis aries), we performed ELISA, qRT-PCR, western blotting and immunohistochemistry. The quantitative results showed that melatonin, HIOMT and AANAT levels in the CL verum were significantly higher than those in oestrous CL (p < 0.05), whereas MT1 and MT2 exhibited no change between the oestrous CL and CL verum (p > 0.05); moreover, the localization results showed that HIOMT, AANAT, MT1 and MT2 were mainly expressed in large luteal cells (LLCs). In summary, the above results suggested that sheep CL has potential for the synthesis of melatonin; meanwhile, they also suggested that CL is one of the targets of melatonin. These results provide not only a basis for whether sheep CL can synthesize melatonin but also provide a reference for further study on the mechanism of melatonin in the CL.

Comprehensive Analysis of MicroRNA⁻Messenger RNA from White Yak Testis Reveals the Differentially Expressed Molecules Involved in Development and Reproduction.

Testis development is a vital and tightly regulated process in mammals. Understanding the biological mechanisms underlying testis development will benefit the animal reproduction industry. Expression changes in microRNA and messenger RNA in response to dynamic regulation effects have been associated with this process. However, very little is known about the roles of these molecules in yak development. Using whole-genome small RNA and messenger RNA sequencing, we performed a comprehensive analysis of the microRNA⁻messenger RNA interaction network expression in the testicles of Tianzhu white yaks during three developmental stages. Using Short Time-series Expression Miner analysis we identified 589 differentially expressed microRNAs (DERs) and 3383 differentially expressed messenger RNAs (DEGs) in the three age groups. A total of 93 unique DEGs are primarily involved in reproduction and testis development. Subsequently, four integration networks were constructed according to the DEGs and DERs in three biological processes. Nineteen DEGs were potentially regulated by 60 DERs, of which miR-574 and target gene AURKA played a crucial role in yak testis development and reproduction. The results of this study provide a basis for further exploration of the microRNA⁻messenger RNA interactions in testis development and reproduction and aid in uncovering the molecular mechanisms of spermatogenesis in male mammals.

Differential proteome association study of freeze-thaw damage in ram sperm.

In this study proteomics analysis was performed to investigate damage caused to ram sperm by the freeze-thaw process. Sperm motility, viability, reactive oxygen species (ROS) and adenosine triphosphate (ATP) content were measured to evaluate sperm quality. Compared with fresh groups, motility, viability and ATP content were all lower in freeze-thawed sperm (P < 0.001), and ROS content was higher (P < 0.001). Moreover, 25 differential protein spots were detected in two-dimensional gels using PDQuest 8.0 software and the corresponding proteins were identified using matrix-assisted laser desorption/ionization tandem time-of-flight mass spectrometry (MALDI-TOF-TOF MS) coupled with searching of the NCBI protein sequence database. Among these proteins, hexokinase1 (HXK1), the enzyme that catalyzes the first step of glycolysis in the sperm glycolytic pathway, is known to be associated with sperm motility. Casein kinase II subunit alpha (CSNK2A2), a serine/threonine-selective protein kinase, is associated with sperm apoptosis. We used immunoblotting and immunofluorescence to analyze the expression and localization of these two proteins. HXK1 and CSNK2A2 expression levels in fresh sperm were significantly higher than that in freeze-thawed sperm (P < 0.001). HXK1 and CSNK2A2 were detected in the main part of the sperm flagellum, and the immunofluorescence signal from these proteins was weakened in the freeze-thawed group. Decreased expression of HXK1 and CSNK2A2 may be associated with decreased sperm motility and viability following freeze-thawing.

Analysis of copy number variations by SNP50 BeadChip array in Chinese sheep.

The sheep (Ovis aries) plays a major socio-economic role in the world. Copy number variations (CNVs) are increasingly recognized as a key and potent source of genetic variation and phenotypic diversity, but little is known about the extent to which CNVs contribute to genetic variation in Chinese sheep breeds. Analyses of CNVs in the genomes of eight sheep breeds were performed using the sheep SNP50 BeadChip genotyping array. A total of 111 CNV regions (CNVRs) were obtained from 160 Chinese sheep breeds. These CNVRs covered 13.75Mb of the sheep genome sequence. A total of 22 Go terms and 17 candidate genes were obtained from the functional analysis. Ten CNVRs were selected for validation, of which 7 CNVRs were further experimentally confirmed by quantitative PCR. Four candidate genes were selected to confirm the results of the functional analysis. These results provide a resource for furthering understanding of ruminant biology, and for further improving the genetic quality of sheep breeds.