Dataset Information

Comprehensive transcriptome and methylome analysis delineates the biological basis of hair follicle development and wool-related traits in Merino sheep.

ABSTRACT:

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, and PTPN1) 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.

SUBMITTER: Zhao B

PROVIDER: S-EPMC8427949 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Wool is the critical textile raw material which is produced by the hair follicle of sheep. Therefore, it has important implications to investigate the molecular mechanism governing hair follicle development. Due to high cellular heterogeneity as well as the insufficient cellular, molecular, and spatial characterization of hair follicles on sheep, the molecular mechanisms involved in hair follicle development and wool curvature of sheep remains largely unknown. Single-cell RNA sequencing (scRNA-seq) technologies have made it possible to comprehensively dissect the cellular composition of complex skin tissues and unveil the differentiation and spatial signatures of epidermal and hair follicle development. However, such studies are lacking so far in sheep. Here, single-cell suspensions from the curly wool and straight wool lambskins were prepared for unbiased scRNA-seq. Based on UAMP dimension reduction analysis, we identified 19 distinct cell populations from 15,830 single-cell transcriptomes and characterized their cellular identity according to specific gene expression profiles. Furthermore, novel marker gene was applied in identifying dermal papilla cells isolated in vitro. By using pseudotime ordering analysis, we constructed the matrix cell lineage differentiation trajectory and revealed the dynamic gene expression profiles of matrix progenitors' commitment to the hair shaft and inner root sheath (IRS) cells. Meanwhile, intercellular communication between mesenchymal and epithelial cells was inferred based on CellChat and the prior knowledge of ligand-receptor pairs. As a result, strong intercellular communication and associated signaling pathways were revealed. Besides, to clarify the molecular mechanism of wool curvature, differentially expressed genes in specific cells between straight wool and curly wool were identified and analyzed. Our findings here provided an unbiased and systematic view of the molecular anatomy of sheep hair follicle comprising 19 clusters; revealed the differentiation, spatial signatures, and intercellular communication underlying sheep hair follicle development; and at the same time revealed the potential molecular mechanism of wool curvature, which will give important new insights into the biology of the sheep hair follicle and has implications for sheep breeding.

Project description:BackgroundWith the aim of identifying selection signals in three Merino sheep lines that are highly specialized for fine wool production (Australian Industry Merino, Australian Merino and Australian Poll Merino) and considering that these lines have been subjected to selection not only for wool traits but also for growth and carcass traits and parasite resistance, we contrasted the OvineSNP50 BeadChip (50 K-chip) pooled genotypes of these Merino lines with the genotypes of a coarse-wool breed, phylogenetically related breed, Spanish Churra dairy sheep. Genome re-sequencing datasets of the two breeds were analyzed to further explore the genetic variation of the regions initially identified as putative selection signals.ResultsBased on the 50 K-chip genotypes, we used the overlapping selection signals (SS) identified by four selection sweep mapping analyses (that detect genetic differentiation, reduced heterozygosity and patterns of haplotype diversity) to define 18 convergence candidate regions (CCR), five associated with positive selection in Australian Merino and the remainder indicating positive selection in Churra. Subsequent analysis of whole-genome sequences from 15 Churra and 13 Merino samples identified 142,400 genetic variants (139,745 bi-allelic SNPs and 2655 indels) within the 18 defined CCR. Annotation of 1291 variants that were significantly associated with breed identity between Churra and Merino samples identified 257 intragenic variants that caused 296 functional annotation variants, 275 of which were located across 31 coding genes. Among these, four synonymous and four missense variants (NPR2_His847Arg, NCAPG_Ser585Phe, LCORL_Asp1214Glu and LCORL_Ile1441Leu) were included.ConclusionsHere, we report the mapping and genetic variation of 18 selection signatures that were identified between Australian Merino and Spanish Churra sheep breeds, which were validated by an additional contrast between Spanish Merino and Churra genotypes. Analysis of whole-genome sequencing datasets allowed us to identify divergent variants that may be viewed as candidates involved in the phenotypic differences for wool, growth and meat production/quality traits between the breeds analyzed. The four missense variants located in the NPR2, NCAPG and LCORL genes may be related to selection sweep regions previously identified and various QTL reported in sheep in relation to growth traits and carcass composition.

Project description:Hormonal assessment tools are important for determining the reproductive success of production animals. This study used non-invasive wool assessment to quantify changes in progesterone and cortisol levels in reproducing female merino sheep. Wool samples were collected from a group of n = 46 maiden merino ewes (22-25 months old), naturally joined under natural light conditions in southern New South Wales (NSW), Australia. Three shearing opportunities were conducted as part of standard on-farm management practices. The wool samples were collected at three different dates during 2017, January (prior to rams being put out with the mob and to provide a baseline level since previous shearing in May 2016), September (during very late stages of gestation-approximately 2 weeks prior to parturition) and December (ewes had given birth and ~2-month-old lambs were at foot). Analysis of cortisol and progesterone was conducted concurrently from the same sample of wool. The hormones in wool samples quantified using commercially available cortisol and progesterone enzyme-immunoassay kits. Wool cortisol concentrations increased significantly (p = 3.04E-14) from pre-joining in January (1.33±0.12 ng/g) to late gestation in September (3.59±0.12 ng/g). Concentration of wool cortisol post-lambing in December (3.27±0.14 ng/g) did not decline significantly (p = 0.124) after gestation however remained significantly higher (p = 3.82E-10) than pre-joining levels. Wool progesterone (PG) concentrations increased significantly (p = 1.83E-33) from pre-joining (0.04±0.005 ng/g) in January to late gestation in September (5.53±0.13 ng/g) with a significant (p = 5.44E-59) decline observed in December (0.05±0.003 ng/g) to post- pregnancy concentrations. No significant difference was shown between pre-joining and post lambing PG concentrations (p = 0.057). Our results showed that non-invasive assessment of hormones in Merino sheep wool reflected significant increase in both cortisol and progesterone guided by pregnancy.

Project description:An individual merino sheep's output of wool production is influenced by synergistic interactions of sheep genetics, climate, farm management, and nutrition available to the whole flock. The price paid to the producer for this wool commodity is determined via numerous tested parameters and /or subjective appraisal of the raw greasy wool. This research investigated the level of variation in wool cortisol (a physiological stress biomarker) and wool micron (MIC) in Merino ewes (Ovis aries), pre-partum and post-lambing (lactation/lambs at foot), using maiden ewe (n = 38) managed in an outdoor paddock in a commercial farm. The key findings of this study are; (1) wool quality indicators showed a significant variation between pre- and post- parturition including significant reduction in MIC and (2) there was a negative correlation between wool cortisol levels and wool micron pre-parturition (rs = - 0.179, p < 0.05). This relationship between wool cortisol and wool micron was positive (rs = + 0.29, p < 0.05) during post-parturition suggesting that ewes with lambs at foot ended up with finer wool (reduction in fibre diameter) but they also maintained high levels of wool cortisol. Furthermore, the comfort factor, curvature, standard deviation and spin fineness of the wool were also significantly reduced post-parturition. The results of this study show that metabolic resources partitioning in ewe associated with pregnancy and lambing can result in a reduction in wool quality indices. The activity of the HPA-axis is attenuated during late gestation and parturition as a maternal adaptation; however, the results of our study show that wool cortisol remained similar between pre- and post- lambing. This result indicates that environmental stressors that may have been operating on farm (e.g., cold winter period) could influence on maternal physiological stress response however the exact level of influence of environment conditions on ewe stress levels and productivity traits (e.g., lambing success and wool quality) warrants further investigation. In conclusion, the use of top-knot wool sampling in combination with wool cortisol analysis provides researchers with a convenient method to quantify wool quality and physiological stress simultaneously under commercial sheep production.