Comprehensive analysis of the expression profiles of mRNA, lncRNA, circRNA, and miRNA in primary hair follicles of coarse sheep fetal skin.

BACKGROUND: The Qinghai Tibetan sheep, a local breed renowned for its long hair, has experienced significant deterioration in wool characteristics due to the absence of systematic breeding practices. Therefore, it is imperative to investigate the molecular mechanisms underlying follicle development in order to genetically enhance wool-related traits and safeguard the sustainable utilization of valuable germplasm resources. However, our understanding of the regulatory roles played by coding and non-coding RNAs in hair follicle development remains largely elusive. RESULTS: A total of 20,874 mRNAs, 25,831 circRNAs, 4087 lncRNAs, and 794 miRNAs were annotated. Among them, we identified 58 DE lncRNAs, 325 DE circRNAs, 924 DE mRNAs, and 228 DE miRNAs during the development of medullary primary hair follicle development. GO and KEGG functional enrichment analyses revealed that the JAK-STAT, TGF-ß, Hedgehog, PPAR, cGMP-PKG signaling pathway play crucial roles in regulating fibroblast and epithelial development during skin and hair follicle induction. Furthermore, the interactive network analysis additionally identified several crucial mRNA, circRNA, and lncRNA molecules associated with the process of primary hair follicle development. Ultimately, by investigating DEmir's role in the ceRNA regulatory network mechanism, we identified 113 circRNA-miRNA pairs and 14 miRNA-mRNA pairs, including IGF2BP1-miR-23-x-novel-circ-01998-MSTRG.7111.3, DPT-miR-370-y-novel-circ-005802-MSTRG.14857.1 and TSPEAR-oar-miR-370-3p-novel-circ-005802- MSTRG.10527.1. CONCLUSIONS: Our study offers novel insights into the distinct expression patterns of various transcription types during hair follicle morphogenesis, establishing a solid foundation for unraveling the molecular mechanisms that drive hair development and providing a scientific basis for selectively breeding desirable wool-related traits in this specific breed.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep.

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Genome-wide identification of candidate copy number polymorphism genes associated with complex traits of Tibetan-sheep.

Copy number variation (CNV) is a genetic structural polymorphism important for phenotypic diversity and important economic traits of livestock breeds, and it plays an important role in the desired genetic variation. This study used whole genome sequencing to detect the CNV variation in the genome of 6 local Tibetan sheep groups. We detected 69,166 CNV events and 7230 copy number variable regions (CNVRs) after merging the overlapping CNVs, accounting for 2.72% of the reference genome. The CNVR length detected ranged from 1.1 to 1693.5 Kb, with a total length of 118.69 Mb and an average length of 16.42 Kb per CNVR. Functional GO cluster analysis showed that the CNVR genes were mainly involved in sensory perception systems, response to stimulus, and signal transduction. Through CNVR-based Vst analysis, we found that the CACNA2D3 and CTBP1 genes related to hypoxia adaptation, the HTR1A gene related to coat color, and the TRNAS-GGA and PIK3C3 genes related to body weight were all strongly selected. The findings of our study will contribute novel insights into the genetic structural variation underlying hypoxia adaptation and economically important traits in Tibetan sheep.

Genetic diversity and selection of Tibetan sheep breeds revealed by whole-genome resequencing.

OBJECTIVE: This study aimed to elucidate the underlying gene regions responsible for productive, phenotypic or adaptive traits in different ecological types of Tibetan sheep and the discovery of important genes encoding valuable traits. METHODS: We used whole-genome resequencing to explore the genetic relationships, phylogenetic tree, and population genetic structure analysis. In addition, we identified 28 representative Tibetan sheep single-nucleotide polymorphisms (SNPs) and genomic selective sweep regions with different traits in Tibetan sheep by fixation index (Fst) and the nucleotide diversity (θπ) ratio. RESULTS: The genetic relationships analysis showed that each breed partitioned into its own clades and had close genetic relationships. We also identified many potential breed-specific selective sweep regions, including genes associated with hypoxic adaptability (MTOR, TRHDE, PDK1, PTPN9, TMTC2, SOX9, EPAS1, PDGFD, SOCS3, TGFBR3), coat color (MITF, MC1R, ERCC2, TCF25, ITCH, TYR, RALY, KIT), wool traits (COL4A2, ERC2, NOTCH2, ROCK1, FGF5, SOX9), and horn phenotypes (RXFP2). In particular, a horn-related gene, RXFP2, showed the four most significantly associated SNP loci (g. 29481646 A>G, g. 29469024 T>C, g. 29462010 C>T, g. 29461968 C>T) and haplotypes. CONCLUSION: This finding demonstrates the potential for genetic markers in future molecular breeding programs to improve selection for horn phenotypes. The results will facilitate the understanding of the genetic basis of production and adaptive unique traits in Chinese indigenous Tibetan sheep taxa and offer a reference for the molecular breeding of Tibetan sheep.

Identification of SSTR5 Gene Polymorphisms and Their Association With Growth Traits in Hulun Buir Sheep.

The aim of this study was to locate SSTR5 polymorphisms and evaluate their association with growth traits in Hulun Buir sheep. The study followed up 884 Hulun Buir sheep from birth to 16 months of age, which were born in the same pasture and the same year, and a consistent grazing management strategy was maintained. The birth weight (BRW) was recorded at birth, and body weight (BW), body height (BH), body length (BL), chest circumference (ChC), chest depth (ChD), chest width (ChW), hip width (HW), and cannon circumference (CaC) were measured at 4 and 9 months of age. BW, BH, BL, ChD, HW, and CaC were also recorded at 16 months of age. Based on the growth traits, 233 sheep were selected as experimental animals. Sanger sequencing was performed, and seven single-nucleotide polymorphisms (SNPs) were identified. Association analyses of the SNPs and the growth traits were then conducted. Seven SNPs of the SSTR5 exhibited moderate polymorphism (0.25

Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep.

The identification of candidate genes and genetic variations associated with growth traits is important for sheep breeding. Insulin like growth factor 1 (IGF1) and insulin like growth factor 1 receptor (IGF1R) are well-accepted candidate genes that affect animal growth and development. The current study attempted to assess the association between IGF1 and IGF1R genetic polymorphisms and growth traits in Hulun Buir sheep. To achieve this goal, we first identified three and ten single nucleotide polymorphisms (SNPs) in exons of IGF1 and IGF1R in Hulun Buir sheep and then constructed six haplotypes of IGF1R based on linkage disequilibrium, respectively. Association studies were performed between SNPs and haplotypes of IGF1 and IGF1R with twelve growth traits in a population encompassing 229 Hulun Buir sheep using a general linear model. Our result indicated three SNPs in IGF1 were significantly associated with four growth traits (p < 0.05). In IGF1R, three SNPs and two haplotype blocks were significantly associated with twelve growth traits (p < 0.05). The combined haplotype H5H5 and H5H6 in IGF1R showed the strong association with 12 superior growth traits in Hulun Buir sheep (p < 0.05). In conclusion, we identified SNPs and haplotype combinations associated with the growth traits, which provided genetic resources for marker-assisted selection (MAS) in Hulun Buir sheep breeding.

Identification of Somatostatin Receptor Subtype 1 (SSTR1) Gene Polymorphism and Their Association with Growth Traits in Hulun Buir Sheep.

This study was conducted to evaluate SSTR1 gene polymorphisms and their association with growth traits in Hulun Buir sheep. We followed 233 Hulun Buir sheep from birth to 16 months of age, born in the same pasture and on the same year under a consistent grazing conditions. The body weight (BW), body height (BH), body length (BL), chest circumference (ChC), chest depth (ChD), chest width (ChW), hip width (HW), and cannon circumference (CaC) were measured and recorded at birth, 4 months, 9 months, and 16 months of age. The polymorphisms of the SSTR1 gene in Hulun Buir sheep were excavated using exon sequencing, and association analyses of between SNPs and growth traits at each growth stage were conducted. The results showed that there were four SNPs in Exon 2 of the SSTR1 gene, SNP1, SNP2, and SNP3 were low mutation sites, and SNP4 was a moderate mutation site. Four SNPs were consistent with Hardy-Weinberg equilibrium, and all of them were synonymous mutations. The association analyses found that the genotypes of SNP2 were significantly associated with WW and BH at 4 months of age, BW, BL, ChC, and HW at 9 months of age (p < 0.05), and extremely significantly associated with ChD at 4 and 9 months of age (p < 0.01). There were significant associations between SNP3 and BH at 9 months of age, between SNP4 and ChD, ChW, and CaC at 9 months of age, and BW and ChC at 16 months of age (p < 0.05). There were no detectable associations with growth traits among the seven haplotypes between the SNP1, 3, and 4 of a strong linkage disequilibrium (p > 0.05). These results indicated that SNP2, SNP3, and SNP4 may be used as molecular markers for growth traits of Hulun Buir sheep.

Comparative transcriptome of reproductive axis in Chinese indigenous sheep with different FecB genotypes and prolificacies.

In the present study, there were 128, 263, and 139 differentially expressed genes (DEGs) from the ovary, pituitary, and hypothalamus, respectively, with the FecB B+/FecB ++ genotype. Results indicated that there was a synchronized, and differential increase in mRNA transcript abundances for genes encoding for ovarian steroidogenesis-signaling proteins in the ovary, TGF-ß-signaling in hypothalamus, cAMP-signaling-transduction, and dopaminergic-synapse pathways in the pituitary. The values for these variables were associated with mean prolificacy in indigenous sheep and results provided indications for candidate regulator proteins that coordinate processes related to prolificacy in ewes. Furthermore, the results when there was evaluation of the FecB-genotype-specific co-expression modules implied that CYP17 was the hub gene connecting the pathways and for modulation of folliculogenesis and ovulation in the ewe with the FecB B+ genotype.