Search

Additional file 6: Fig. S3. Go and KEGG network diagram during the hair follicle morphogenesis. (a-f) represents the comparison group of G1-G6. Circle is represented gene, square is represented KEGG; triangle is represented GO term.

Additional file 4: Figure S3. Enrichment analysis of miRNA target gene KEGG pathway at different stages of hair follicle development.

Additional file 3: Figure S2. Analysis of miRNA target genes go enrichment at different stages of hair follicle development. (a) GO enrichment analysis of E65 vs. E85 target genes. (b) GO enrichment analysis of E85 vs. E105 target genes. (c) GO enrichment analysis of E105 vs. E135 target genes. (d) GO enrichment analysis of E135 vs. D7 target genes. (e) GO enrichment analysis of D7 vs. D30 target genes. (f) GO enrichment analysis of E65 vs. D30 target genes. where red is BP (biological processes), green is CC (Cellular components) and blue is MF (Molecular functions).

Additional file 7: Fig. S4. K-means clustering analysis of differentially expressed mRNAs among the six comparison groups.

Additional file 5: Figure S4. (a) Western blot level of ACVR1B (target strip in red box dimension). (b) Western blot level of GAPDH (target strip in red box dimension). (c) Western blot level of WNT10A (target strip in red box dimension). (d) Western blot level of ACVR1B (target strip in red box dimension). All the images of western blots are cut prior to hybridization with antibodies. Because the markers are marked on both sides of the glue, the notch in the upper right corner (upper left corner) is to mark the front and back.

Additional file 2: Figure S1. The distribution of the read from the sequencing data. (a) Read redundancy statistics of different samples. (b) Read length distribution of different samples.

Additional file 2: Fig. S1. Heatmap of DEGs during the hair follicle morphogenesis. The x-axis represents the sample number, where S1 to S3 is G1; S4 to S6 is G2; S7 to S9 is G3; S10 to S12 is G4; S13 to S15 is G5; S16 to S18 is G6.

Hair follicle growth and development are a complex and long-term physiological process, which is regulated by a variety of physical factors and signal pathways. Increasing the understanding of the epigenetic regulation and function of candidate genes related to hair follicle development will help to better understand the molecular regulatory mechanisms of hair follicle development. In this study, the methylated DNA immunoprecipitation sequencing (MeDIP-seq) was used to obtain the genome-wide methylation map of the hair follicular development of Super Merino sheep in six stages (fetal skin tissue at 65d, 85d, 105d, 135d, 7d, and 30d after birth). Combined with the results of previous RNA-sequencing, 65 genes were screened out that were both differential methylation and differential expression, including EDN1, LAMC2, NR1D1, RORB, MyOZ3, and WNT2 gene. Differential methylation genes were enriched in Wnt, TNF, TGF-beta, and other signaling pathways related to hair follicle development. The bisulfite sequencing PCR results and MeDIP-seq were basically consistent, indicating that the sequencing results were accurate. As a key gene in the Wnt signaling pathway, both differential methylation and expression gene identified by MeDIP-seq and RNA-seq, further exploration of the function of WNT2 gene revealed that the DNA methylation of exon 5 (CpG11 site) promoted the expression of WNT2 gene. The overexpression vector of lentivirus pLEX-MCS-WNT2 was constructed, and WNT2 gene effectively promoted the proliferation of sheep skin fibroblasts. The results showed that WNT2 gene could promote the growth and development of skin and hair follicles. The results of this study will provide a theoretical basis for further research on sheep hair follicle development and gene regulation mechanisms.

BackgroundCharacterization 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.ResultsWe 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.ConclusionsOur 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 Genetic improvement of wool and growth traits is a major goal in the sheep industry, but their underlying genetic architecture remains elusive. To improve our understanding of these mechanisms, we conducted a weighted single-step genome-wide association study (WssGWAS) and then integrated the results with large-scale transcriptome data for five wool traits and one growth trait in Merino sheep: mean fibre diameter (MFD), coefficient of variation of the fibre diameter (CVFD), crimp number (CN), mean staple length (MSL), greasy fleece weight (GFW), and live weight (LW). Results Our dataset comprised 7135 individuals with phenotype data, among which 1217 had high-density (HD) genotype data (n = 372,534). The genotypes of 707 of these animals were imputed from the Illumina Ovine single nucleotide polymorphism (SNP) 54 BeadChip to the HD Array. The heritability of these traits ranged from 0.05 (CVFD) to 0.36 (MFD), and between-trait genetic correlations ranged from − 0.44 (CN vs. LW) to 0.77 (GFW vs. LW). By integrating the GWAS signals with RNA-seq data from 500 samples (representing 87 tissue types from 16 animals), we detected tissues that were relevant to each of the six traits, e.g. liver, muscle and the gastrointestinal (GI) tract were the most relevant tissues for LW, and leukocytes and macrophages were the most relevant cells for CN. For the six traits, 54 quantitative trait loci (QTL) were identified covering 81 candidate genes on 21 ovine autosomes. Multiple candidate genes showed strong tissue-specific expression, e.g. BNC1 (associated with MFD) and CHRNB1 (LW) were specifically expressed in skin and muscle, respectively. By conducting phenome-wide association studies (PheWAS) in humans, we found that orthologues of several of these candidate genes were significantly (FDR BNC1 was significantly associated with MFD in sheep and with hair colour in humans, and CHRNB1 was significantly associated with LW in sheep and with body mass index in humans. Conclusions Our findings provide novel insights into the biological and genetic mechanisms underlying wool and growth traits, and thus will contribute to the genetic improvement and gene mapping of complex traits in sheep.

As a member of the Frizzled family, Frizzled3 (FZD3) is a receptor of the canonical Wnt signaling pathway and plays a vital role in mammalian hair follicle developmental processes. However, its effects on wool traits are not clear. The objectives of this study were to identify the single nucleotide polymorphisms (SNPs) and the expression patterns of FZD3 gene, and then to determine whether it affected wool traits of Chinese Merino sheep (Xinjiang Type) or not. PCR-single stranded conformational polymorphism (PCR-SSCP) and sequencing were used to identify mutation loci, and general linear model (GLM) with SAS 9.1 was used for the association analysis between wool traits and SNPs. Quantitative real-time PCR (qRT-PCR) was used to investigate FZD3 gene expression levels. The results showed that six exons of FZD3 gene were amplified and two mutation loci were identified in exon 1 (NC_019459.2: g.101771685 T>C (SNP1)) and exon 3 (NC_019459.2: g.101810848, A>C (SNP2)), respectively. Association analysis showed that SNP1 was significantly associated with mean fiber diameter (MFD) (P=0.04) and live weight (LW) (P=0.0004), SNP2 was significantly associated with greasy fleece weight (GFW) (P=0.04). The expression level of FZD3 gene in skin tissues of the superfine wool (SF) group was significantly lower (P

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&ndash, 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 h2 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.

Objective The purpose of this study was to investigate the genetic effects of six keratin (KRT) genes on the wool traits of 418 Chinese Merino (Xinjiang type) (CMXT) individuals. Methods To explore the effects and association of six KRT genes on sheep wool traits, The polymerase chain reaction-based single-strand conformation polymorphism (PCR-SSCP), DNA sequencing, and the gene pyramiding effect methods were used. Results We report 20 mutation sites (single-nucleotide polymorphisms) within the six KRT genes, in which twelve induced silent mutations; five induced missense mutations and resulted in Ile→Thr, Glu→Asp, Gly→Ala, Ala→Ser, Se→His; two were nonsense mutations and one was a same-sense mutation. Association analysis showed that two genotypes of the KRT31 gene were significantly associated with fiber diameter (p

This study was conducted to estimate variance components and genetic parameters for 305 days mature equivalent milk yields (MEM) and calving interval (CAI) of registered Xinjiang Brown cattle. The total records were 3940 including 2579 for 305-day MEM and 1970 for CAI, which were collected from Xinjiang Urumqi Cattle Breeding farm in China with calving records from 1990 to 2008. Genetic parameters were obtained by multiple trait derivative-free restricted maximum likelihood (MTDFREML) using animal model. The model included year, season, parity and calving interval of calving for 305 days MEM, and year, season and parity of birth for CAI as fixed effects. Heritability for 305 days MEM was moderate (0.39) and in the same range of parameters estimated in management systems with medium production levels. Heritability of calving interval was small (0.02) as fertility traits for Xinjiang Brown cattle. Estimates of genetic and environmental correlations between 305 days MEM and CAI were 0.47 and 0.37, respectively. Estimates of genetic variation and heritability indicated that selection would result in genetic improvement of production traits. Estimates of both heritability and genetic variation for CAI were small, which indicates that genetic improvement would be difficult.

Long non-coding RNAs (lncRNAs), >200 nt in length, are transcribed from mammalian genomes. They play important regulatory roles in various biological processes; However, the function and expression profile of lncRNAs involved in the development of hair follicles in the fetus, have been relatively under-explored area. To investigate the specific role of lncRNAs and mRNAs that regulate hair follicle development, we herein performed a comprehensive study on the lncRNA and mRNA expression profiles of sheep at multiple embryonic days (E65, E85, E105, and E135) and six lambs aged one week (D7) and one month (D30) using RNA-seq technology. The number of genes (471 lncRNAs and 12,812 mRNAs) differentially expressed and potential targets of differentially expressed lncRNAs were predicted. Differentially expressed lncRNAs were grouped into 10 clusters based on their expression pattern by K-means clustering. Moreover, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that some differentially expressed mRNAs, such as DKK1, DSG4, FOXE1, Hoxc13, SFRP1, SFRP2, and Wnt10A overlapped with lncRNAs targets, and enriched in important hair follicle developmental pathways, including Wnt, TNF, and MAPK signaling pathways. In addition, 9 differentially expressed lncRNAs and 4 differentially expressed mRNAs were validated using quantitative real-time PCR (qRT-PCR). This study helps enrich the Ovis lncRNA databases and provides a comprehensive lncRNA transcriptome profile of fetal and postnatal skin of sheep. Additionally, it provides a foundation for further experiments on the role of lncRNAs in the regulation of hair growth in sheep.