Background: Characterization of the molecular mechanisms underlying hair follicle development is of paramount importance in the genetic improvement of wool-related traits in sheep and skin-related traits in humans. The Merino is the most important breed of fine-wooled sheep in the world. In this study, we systematically investigated the complexity of sheep hair follicle development by integrating transcriptome and methylome datasets from Merino sheep skin.Results: We analysed 72 sequence datasets, including DNA methylome and the whole transcriptome of four gene types, i.e. protein-coding genes (PCGs), lncRNAs, circRNAs, and miRNAs, across four embryonic days (E65, E85, E105, and E135) and two postnatal days (P7 and P30) from the skin tissue of 18 Merino sheep. We revealed distinct expression profiles of these four gene types across six hair follicle developmental stages, and demonstrated their complex interactions with DNA methylation. PCGs with stage-specific expression or regulated by stage-specific lncRNAs, circRNAs, and miRNAs were significantly enriched in epithelial differentiation and hair follicle morphogenesis. Regulatory network and gene co-expression analyses identified key transcripts controlling hair follicle development. We further predicted transcriptional factors (e.g. KLF4, LEF1, HOXC13, RBPJ, VDR, RARA, and STAT3) with stage-specific involvement in hair follicle morphogenesis. Through integrating these stage-specific genomic features with results from genome-wide association studies (GWAS) of five wool-related traits in 7135 Merino sheep, we detected developmental stages and genes that were relevant with wool-related traits in sheep. For instance, genes that were specifically upregulated at E105 were significantly associated with most of wool-related traits. A phenome-wide association study (PheWAS) demonstrated that candidate genes of wool-related traits (e.g.SPHK1,GHR,PPP1R27,CSRP2,EEF1A2, andPTPN1) in sheep were also significantly associated with dermatological, metabolic, and immune traits in humans.Conclusions: Our study provides novel insights into the molecular basis of hair follicle morphogenesis and will serve as a foundation to improve breeding for wool traits in sheep. It also indicates the importance of studying gene expression in the normal development of organs in understanding the genetic architecture of economically important traits in livestock. The datasets generated here are useful resources for functionally annotating the sheep genome, and for elucidating early skin development in mammals, including humans.

Background: Tibetan cashmere goats are served as a valuable model for high altitude adaptation and hypoxia complications related studies, while the cashmere produced by these goats is an important source of income for the herders. The aim of this study was to investigate the differences in protein abundance underlying the fine (average 12.20 ± 0.03 μm of mean fiber diameter) and coarse cashmere (average 14.67 ± 0.05 μm of mean fiber diameter) producing by Tibetan cashmere goats. We systematically investigated the genetic determinants of fiber diameter by integrated analysis with proteomic and transcriptomic datasets from skin tissues of Tibetan cashmere goats. Results: We identified 1980 proteins using a label-free proteomics approach. They were annotated to three different databases, while 1730 proteins were mapped to the original protein coding genes (PCGs) of the transcriptomic study. Comparative analyses of cashmere with extremely fine vs. coarse phenotypes yielded 29 differentially expressed proteins (DEPs), for instance, APOH, GANAB, AEBP1, CP, CPB2, GPR142, VTN, IMPA1, CTSZ, GLB1, and HMCN1. Functional enrichment analysis of these DEPs revealed their involvement in oxidation-reduction process, cell redox homeostasis, metabolic, PI3K-Akt, MAPK, and Wnt signaling pathways. Transcription factors enrichment analysis revealed the proteins mainly belong to NF-YB family, HMG family, CSD family. We further validated the protein abundance of four DEPs (GC, VTN, AEBP1, and GPR142) through western blot, and considered they were the most potential candidate genes for cashmere traits in Tibetan cashmere goats. Conclusions: These analyses indicated that the major biological variations underlying the difference of cashmere fiber diameter in Tibetan cashmere goats were attributed to the inherent adaptations related to metabolic, hypoxic, and stress response differences. This study provided novel insights into the breeding strategies for cashmere traits and enhance the understanding of the biological and genetic mechanisms of cashmere traits in Tibetan cashmere goats.

Chicken is a valuable model for understanding fundamental biology, vertebrate evolution and diseases, as well as a major source of nutrient-dense and lean-protein-enriched food globally. Although it is the first non-mammalian amniote genome to be sequenced, the chicken genome still lacks a systematic characterization of functional impacts of genetic variants. Here, through integrating 7,015 RNA-Seq and 2,869 whole-genome sequence data, the Chicken Genotype- Tissue Expression (ChickenGTEx) project presents the pilot reference of regulatory variants in 28 chicken tissue transcriptomes, including millions of regulatory effects on primary expression (including protein-coding genes, lncRNA and exon) and post-transcriptional modifications (alternative splicing and 3’ untranslated region alternative polyadenylation). We explored the tissue-sharing and context-specificity of these regulatory variants, their underlying molecular mechanisms of action, and their utility in interpreting adaptation and genome-wide associations of 108 chicken complex traits. Finally, we illustrated shared and lineage-specific features of gene regulation between chickens and mammals, and demonstrated how the ChickenGTEx resource can further assist with translating genetic findings across species. One-Sentence Summary: The ChickenGTEx provides a multi-tissue reference of regulatory variants for chicken genetics and genomics, functional genomics, precision breeding, veterinary medicine, vertebrate evolution and even human biomedicine.

Genomic evaluations are a method for improving the accuracy of breeding value estimation. This study aimed to compare estimates of genetic parameters and the accuracy of breeding values for wool traits in Merino sheep between pedigree-based best linear unbiased prediction (PBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP) using Bayesian inference. Data were collected from 28,391 yearlings of Chinese Merino sheep (classified in 1992–2018) at the Xinjiang Gonaisi Fine Wool Sheep-Breeding Farm, China. Subjectively-assessed wool traits, namely, spinning count (SC), crimp definition (CRIM), oil (OIL), and body size (BS), and objectively-measured traits, namely, fleece length (FL), greasy fleece weight (GFW), mean fiber diameter (MFD), crimp number (CN), and body weight pre-shearing (BWPS), were analyzed. The estimates of heritability for wool traits were low to moderate. The largest h² values were observed for FL (0.277) and MFD (0.290) with ssGBLUP. The heritabilities estimated for wool traits with ssGBLUP were slightly higher than those obtained with PBLUP. The accuracies of breeding values were low to moderate, ranging from 0.362 to 0.573 for the whole population and from 0.318 to 0.676 for the genotyped subpopulation. The correlation between the estimated breeding values (EBVs) and genomic EBVs (GEBVs) ranged from 0.717 to 0.862 for the whole population, and the relative increase in accuracy when comparing EBVs with GEBVs ranged from 0.372% to 7.486% for these traits. However, in the genotyped population, the rank correlation between the estimates obtained with PBLUP and ssGBLUP was reduced to 0.525 to 0.769, with increases in average accuracy of 3.016% to 11.736% for the GEBVs in relation to the EBVs. Thus, genomic information could allow us to more accurately estimate the relationships between animals and improve estimates of heritability and the accuracy of breeding values by ssGBLUP.

Pituitary gonadotropins perform essential functions in mammalian reproduction by stimulating gametogenesis and steroidogenesis in the ovaries and testicles. EZH2 is a histone methyltransferase that inhibits proliferation and aggravates apoptosis in stem cells subjected to pathological stimuli. However, the expression and molecular mechanisms of EZH2 in pituitary cells in vitro have not been extensively studied. In this study, the relative abundances of EZH2 mRNA (p < 0.01) and protein (p < 0.05) expression were larger in the pituitary cells of Hu sheep with relatively greater fecundity (GF) compared to those with lesser fecundity (LF). Loss-of-function examinations demonstrated that EZH2 gene knockdown led to an earlier induction of apoptosis in sheep pituitary cells (PCs). The relative abundance of CASP3, CASP9, and BAX was increased (p < 0.01), while BCL2’s abundance was less decreased (p < 0.01) in PCs where there was EZH2 gene knockdown. Additionally, cell proliferation (p < 0.01) and viability (p < 0.01) were decreased in EZH2-knockdown sheep PCs, and the cell cycle was blocked compared to a negative control (NC). Notably, EZH2 gene knockdown led to reduced abundances of gonadotropin subunit gene transcripts (FSHβ, p < 0.05) and reduced FSH release (p < 0.01) from PCs. EZH2 gene knockdown led to reduced phosphorylation of AKT, ERK, and mTOR (p < 0.01). The results suggest that EZH2 regulates pituitary cell proliferation, apoptosis, and FSH secretion through modulation of the AKT/ERK signaling pathway, providing a foundation for further study of pituitary cell functions.

Understanding the genetic characteristics of indigenous goat breeds is vital for their conservation and breeding. Haimen goats, native to China’s Yangtze River Delta, possess distinctive traits such as white hair, moderate growth rate, high-quality meat, and small body size. However, knowledge regarding the genetic structure and germplasm characteristics of Haimen goats remains limited. In this study, we performed 20× whole-genome resequencing of 90 goats (60 Haimen goats and 30 Boer goats) to identify single-nucleotide polymorphisms (SNPs) and insertions/deletions (Indels) associated with growth traits. Here, we analyzed population genetic structure and genome-wide selection signatures between the Haimen and Boer goats based on whole-genome resequencing data. The principal component analysis (PCA) and neighbor-joining (N-J) tree results demonstrated significant genetic differentiation between the Haimen and Boer goats. The nucleotide diversity (Pi) and linkage disequilibrium (LD) decay results indicated higher genomic diversity in the Haimen goat population. Furthermore, selective sweep analysis identified candidate genes associated with growth traits. These genes exhibited strong selection signatures and were related to body size (DONSON, BMPR1B, and EPHA5), muscle development (GART, VGLL3, MYH15), and fat metabolism (ADAMTS5, LRP6, XDH, CPT1A, and GPD1). We also identified growth-related candidate genes (NCOR1, DPP6, NOTCH2, and FGGY) specific to Haimen goats. Among these genes, pancreatic lipase-related protein 1 (PNLIPRP1) emerged as the primary candidate gene influencing growth phenotypes. Further analysis revealed that a 26 bp Indel in PNLIPRP1 increased its gene expression, suggesting that this Indel could serve as a molecular marker for early marker-assisted selection, potentially enhancing early growth in goats. These findings provide valuable molecular markers and candidate genes for improving growth traits in Haimen goat breeding.