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.

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.

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.

Dynamical Expression of MicroRNA-127-3p in Proliferating and Differentiating C2C12 Cells.

MicroRNAs (miRNAs) are highly conserved, short non-coding RNAs that regulate gene expression at the posttranscriptional level. Although many miRNAs are identified in muscles and muscle cells, their individual roles are still not fully understood. In the present study, we investigated a muscle highly-expressed miRNA, miR-127-3p, in C2C12 myoblasts and tissues of goats with different muscle phenotypes (Boer vs Wushan black goats). Our results demonstrated that i) miR-127-3p was extensively expressed in tissues of goats; ii) miR-127-3p was higher expressed in muscle, spleen, heart, and skin in the muscular goats (Boer goats) than the control (Wushan black goats). Then we further characterized the dynamical expression of miR-127-3p, MyoD, MyoG, Myf5, Mef2c, and Myosin in the proliferating and differentiating C2C12 myoblasts at day of 0, 1, 3, 5, and 7 in culture mediums. Especially, we found that miR-127-3p was significantly higher expressed in the proliferating than differentiating cells. Our findings suggest that miR-127-3p probably plays roles in the proliferation and differentiation of myoblasts, which further underlies regulation of muscle phenotype in goats.

Identification of copy number variations in Qinchuan cattle using BovineHD Genotyping Beadchip array.

In recent years, copy number variations (CNVs), which associate with complex traits such as disease and quantitative phenotypes, are increasingly recognized as an important and abundant source of genetic variation and phenotypic diversity. CNVs have been studied in several breeds of cattle with the goal of improving selection methods for commercial use; however, little is known about the extent to which CNVs contribute to genetic variation in Qinchuan cattle. The BovineHD Genotyping BeadChip array was used for analyzing the whole genomic CNVs of Qinchuan cattle breed; we discovered 367 unique CNV events from 6 Qinchuan cattle. Accounting for overlapping regions, a total of 365 autosomal copy number variation regions (CNVRs) (131 losses and 234 gains) were identified with an average number of 60.8 CNV events per individual, which covered 13.13 Mb of the cattle genomic sequence corresponding to 0.4 % of the whole cattle genome. The average and median sizes of CNVRs were 35.07 and 18.56 kb, respectively. The CNVRs map of Qinchuan cattle was first constructed based on the BovineHD Genotyping Beadchip array. Functional analysis indicated that most genes in CNVRs that were significantly enriched are involved in environmental stress. Comparison of CNVRs in ten published studies and the 365 CNVRs identified in our study overlapped 0.7-42.7 %. These findings are the first report of CNVs mapping in Qinchuan cattle and contribute to the greater understanding of CNV genetics in commercial cattle phenotypes.

Polymorphisms of the Ovine BMPR-IB, BMP-15 and FSHR and Their Associations with Litter Size in Two Chinese Indigenous Sheep Breeds.

The Small Tailed Han sheep and Hu sheep are two prolific local sheep in China. In this study, the polymorphisms of BMPR-IB (Bone morphogenetic protein receptor IB), BMP-15 (Bone morphogenetic protein 15) and FSHR (follicle stimulating hormone receptor) were investigated to check whether they are associated with litter size in Small Tailed Han sheep and Hu sheep. Consequently, three polymorphisms, FecB mutation in BMPR-IB (c.746A>G), FecG mutation in BMP-15 (c.718C>T) and the mutation (g. 47C>T) in FSHR were found in the above two sheep breeds with a total number of 1630 individuals. The single marker association analysis showed that the three mutations were significantly associated with litter size. The ewes with genotype FecBB/FecBB and FecBB/FecB+ had 0.78 and 0.58 more lambs (p < 0.01) than those with genotype FecB+/FecB+, respectively. The heterozygous Han and Hu ewes with FecXG/FecX+ genotype showed 0.30 (p = 0.05) more lambs than those with the FecX+/FecX+ genotype. For FSHR gene, the ewes with genotype CC had 0.52 (p < 0.01) and 0.75 (p < 0.01) more lambs than those with genotypes TC and TT, respectively. Combined effect analyses indicated an extremely significant interaction (p < 0.01) between the random combinations of BMPR-IB, BMP-15 and FSHR genes on litter size. In addition, the Han and Hu ewes with BB/G+/CC genotype harbor the highest litter size among ewes analyzed in current study. In conclusion, BMPR-IB, BMP-15 and FSHR polymorphisms could be used as genetic markers in multi-gene pyramiding for improving litter size in sheep husbandry.

Expression of Ghrelin in gastrointestinal tract and the effect of early weaning on Ghrelin expression in lambs.

Ghrelin, an endogenous ligand for growth hormone secretagogue receptor, plays an important role in stimulating hormone secretion, development of gastrointestinal tract, food intake and regulating energy balance of animals. In this study we isolated the cDNA sequence of ovine Ghrelin from the abomasums of 7-day-aged lambs. Real-time PCR was used to determine the abundance of Ghrelin mRNA in lamb gastrointestinal tract, and analyze the development changes of abomasums Ghrelin mRNA expression in 0-56 days lambs, as well as find the effects of 42-day weaning on Ghrelin mRNA expression in lamb abomasums. The results showed that: (1) Ghrelin mRNA was expressed widely in gastrointestinal tract and was significantly higher in the abomasums than in other tissues (rumen, reticulum, omasum, duodenum, jejunum, ileum) (P < 0.01); (2) The expression of abomasums Ghrelin mRNA in lamb increased with the growth of age, it reached a plateau at the age of 49 days, however, got a slightly decrease at the age of 56 days; (3) The expression of abomasums Ghrelin mRNA of the 42 days-weaned groups were significantly lower than the no-weaned groups (P < 0.05), and the Ghrelin mRNA expression of the two treatments reached a maximum at the age of 49 days; (4) Correlation analysis indicated that the linear correlativity between abomasums Ghrelin mRNA expression and abomasums weight was very prominent (R(2) = 0.647, P = 0.009). Our results suggested that ovine Ghrelin gene may play an important role in the development of lamb abomasums and 42-day weaning could down regulate the expression of abomasum Ghrelin mRNA, but the mechanism of these needs further research.

Molecular characterization, expression profiles of the ovine FSHR gene and its association with litter size.

Follicle-stimulating hormone receptor (FSHR) is a glycoprotein family member of G proteins coupled receptor super-family, and plays important roles in animal follicular development. In this study, we cloned the 5' flanking region of ovine FSHR gene and analyzed its genomic structure. One special domain named the LRRNT (Leucine rich repeat N-terminal domain) was found in the prediction amino acid sequence of ovine FSHR. RT-qPCR showed that ovine FSHR was expressed widely in detected tissues and was significantly higher in sexual gland (testicle and ovary) than in other tissues (heart, liver, spleen, lung, kidney, rumen, duodenum, muscle, fat, hypothalamus and pituitary) (P < 0.05). In addition, synonymous mutation g.-47C > T of the FSHR gene was detected and confirmed to be significantly associated with the litter size (P < 0.01), the genotype of CC had 0.42 (P < 0.01) and 0.53 (P < 0.01) lambs more than TC and TT genotype, respectively. Our results indicated the association of ovine FSHR with the litter size in sheep and FSHR could be used as a candidate gene for improving reproductive traits in industrial sheep breeding.

Molecular cloning, bioinformatics analysis and expression of insulin-like growth factor 2 from Tianzhu white yak, Bos grunniens.

The IGF family is essential for normal embryonic and postnatal development and plays important roles in the immune system, myogenesis, bone metabolism and other physiological functions, which makes the study of its structure and biological characteristics important. Tianzhu white yak (Bos grunniens) domesticated under alpine hypoxia environments, is well adapted to survive and grow against severe hypoxia and cold temperatures for extended periods. In this study, a full coding sequence of the IGF2 gene of Tianzhu white yak was amplified by reverse transcription PCR and rapid-amplification of cDNA ends (RACE) for the first time. The cDNA sequence revealed an open reading frame of 450 nucleotides, encoding a protein with 179 amino acids. Its expression in different tissues was also studied by Real time PCR. Phylogenetic tree analysis indicated that yak IGF2 was similar to Bos taurus, and 3D structure showed high similarity with the human IGF2. The putative full CDS of yak IGF2 was amplified by PCR in five tissues, and cDNA sequence analysis showed high homology to bovine IGF2. Moreover the super secondary structure prediction showed a similar 3D structure with human IGF2. Its conservation in sequence and structure has facilitated research on IGF2 and its physiological function in yak.

Effects of adding Allium mongolicum Regel powder and yeast cultures to diet on rumen microbial flora of Tibetan sheep (Ovis aries).

The purpose of this experiment was to study the effect of Allium mongolicum Regel powder (AMR) and yeast cultures (YC) on rumen microbial diversity in Tibetan sheep in different Ecological niches. A total of 40 male Tibetan lambs with an initial weight of 18.56 ± 1.49 kg (6 months old) were selected and divided into four groups (10 sheep/pen; n = 10). In the Control Group, each animal was grazed for 8 h per day, in Group I, each animal was supplemented with 200 g of concentrate per day, in Group II, each animal was supplemented with 200 g of concentrate and 10 g of AMR per day, in Group III, each animal was supplemented with 200 g of concentrate and 20 g of YC per day. The experiment lasted 82 days and consisted of a 7-day per-feeding period and a 75-day formal period. The results indicated that at the phylum level, the abundance of Bacteroidota and Verrucomimicrobiota in L-Group II and L-Group III was increased, while the abundance of Proteobacteria was decreased in the LA (Liquid-Associated) groups. The proportion of F/B in S-Group II and S-Group III was increased compared to S-Group I and S-CON in the SA (Soild-Associated) group. At the genus level, the abundance of uncultured_rumen_bacterium and Eubacterium_ruminantium_group in L-Group II and L-Group III was increased. Furthermore, while the abundance of Rikenellaceae_RC9_gut_group was decreased in the LA, the abundance of Prevotella and Eubacterium_ruminantium_group was increased in the S-Group II and S-Group III compared to S-Group I and S-CON. The abundance of probable_genus_10 was the highest in S-Group II in the SA group. After the addition of YC and AMR, there was an increase in rumen microbial abundance, which was found to be beneficial for the stability of rumen flora and had a positive impact on rumen health.

Rumen Development of Tianhua Mutton Sheep Was Better than That of Gansu Alpine Fine Wool Sheep under Grazing Conditions.

The purpose of this experiment was to investigate the differences in rumen tissue morphology, volatile fatty acid content, and rumen microflora between Tianhua mutton sheep and Gansu alpine fine wool sheep under the same grazing conditions. Twelve 30-day-old lambs were randomly selected from two different flocks in Duolong Village and grazed together for a period of 150 days. The rumen tissue was fixed with 4% paraformaldehyde and brought back to the laboratory for H&E staining, the volatile fatty acid content of the rumen contents was detected by gas chromatography, and the rumen flora structure was sequenced by full-length sequencing of the bacterial 16S rRNA gene using the PacBio sequencing platform. The acetic acid and total acid contents of the rumen contents of Tianhua mutton sheep were significantly higher than those of Gansu alpine fine wool sheep (p < 0.05). The rumen papillae height of Tianhua mutton sheep was significantly higher than that of Gansu alpine fine wool sheep (p < 0.05). The diversity and richness of the rumen flora of Tianhua mutton sheep were higher than those of Gansu alpine fine wool sheep, and Beta analysis showed that the microflora structure of the two fine wool sheep was significantly different. At the phylum level, Firmicutes and Bacteroidetes dominated the rumen flora of Tianhua mutton sheep and Gansu alpine fine wool sheep. At the genus level, the dominant strains were Christensenellaceae_R_7_group and Rikenellaceae_RC9_gut_group. LEfSe analysis showed that Prevotella was a highly abundant differential species in Tianhua mutton sheep and lachnospiraccac was a highly abundant differential species in Gansu alpine fine wool sheep. Finally, both the KEGG and COG databases showed that the enrichment of biometabolic pathways, such as replication and repair and translation, were significantly higher in Tianhua mutton sheep than in Gansu alpine fine wool sheep (p < 0.05). In general, there were some similarities between Tianhua mutton sheep and Gansu alpine fine wool sheep in the rumen tissue morphology, rumen fermentation ability, and rumen flora structure. However, Tianhua mutton sheep had a better performance in the rumen acetic acid content, rumen papillae height, and beneficial bacteria content. These differences may be one of the reasons why Tianhua mutton sheep are more suitable for growing in alpine pastoral areas than Gansu alpine fine wool sheep.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations.

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

Prickly ash seeds can promote healthy production of sheep by regulating the rumen microbial community.

This study aimed to investigate the effect of prickly ash seeds (PAS) on the microbial community found in rumen microbes of Hu sheep by adding different percentages of prickly ash seeds and to carry out research on the relation between rumen flora and production performance. Twenty-seven male lambs of Hu sheep were classified into three groups based on the content of prickly ash seeds (PAS) fed for 90 days, i.e., 0%, 3%, and 6%. At the end of the feeding trial, rumen fluid samples were collected from six sheep in each group for 16S amplicon sequencing. The results showed that the addition of prickly ash seeds significantly increased both Chao1 and ACE indices (P < 0.05), and the differences between groups were greater than those within groups. The relative content of Bacteriodota decreased, and the relative content of Fusobacteriota, Proteobacteria, Acidobacteriota, and Euryarchaeota increased. The relative content of Papillibacter and Saccharofermentans was increased at the genus level, and the relative content of Bacteroides and Ruminococcus was decreased. The test group given 3% of prickly ash seeds was superior to the test group given 6% of prickly ash seeds. In addition, the addition of 3% of prickly ash seeds improved the metabolism or immunity of sheep. Fusobacteriota and Acidobacteriota were positively correlated with total weight, dressing percentage, and average daily gain (ADG) and negatively correlated with average daily feed intake (ADFI), feed-to-gain ratio (F/G), and lightness (L*). Methanobrevibacter and Saccharofermentans were positively correlated with ADG and negatively correlated with ADFI and L*. In conclusion, under the present experimental conditions, the addition of prickly ash seeds increased the abundance and diversity of rumen microorganisms in Hu sheep and changed the relative abundance of some genera. However, the addition of 6% prickly ash seeds may negatively affect the digestive and immune functions in sheep rumen.