Comparative study on seasonal hair follicle cycling by analysis of the transcriptomes from cashmere and milk goats

Comparative ultrastructural and transcriptomic profile analysis of skin tissues from indigenous, improved meat, and dairy goat breeds

BACKGROUND

High-quality goatskins are valuable byproducts usually produced by indigenous goat breeds with poorer production performance in Asia and Africa. However, the genetic and molecular mechanisms underpinning goatskin's biomechanical properties (e.g., tensile strength) remain elusive. Mechanistic exploration of these traits could greatly aid the genetic improvement and genetic resource conservation of native breeds in these regions. To fulfill this purpose, we collected skin tissues from three goat breeds: Huai goat (HG), a Chinese native variety producing high-quality goatskins; Yudong meat goat (YDMG), a crossbreed of HG and Boer goat; Henan dairy goat (HNDG), a dairy goat breed.

RESULTS

Scanning electronic microscopy analysis of skin tissues found that the collagen fiber diameters, collagen fibril diameters, and crimps significantly differed among the three goat breeds; however, collagen fibril diameters are similar in HG and HNDG. A sum of 230, 775, and 86 differentially expressed genes (DEGs) were identified from YDMG versus HNDG, HG versus HNDG, and YDMG versus HG, respectively. Functional enrichment analysis suggested that signaling pathways involved in fatty acid, retinol, steroid metabolisms, and GO items related to the physical properties of the skin (e.g., collagen-containing extracellular matrix) are significantly overrepresented in DEGs identified from meat versus dairy goats. Furthermore, 106 DEGs (e.g., COL1A1, COL1A2, and SPARC) showed specific expression patterns in HG and YDMG versus HNDG. Items about biophysical features of skin (e.g., extracellular matrix organization and ECM proteoglycans) are markedly enriched. Protein-protein interaction analysis suggested that two growth factors (IGF1 and PDGFD) are latent collagen and other ECM protein expression modulators.

CONCLUSION

Ultrastructural analysis of goat skin tissues suggested that collagen fibril diameter is not a major factor affecting goatskin quality. Transcriptomic profiles unveiled core genes and associated biological processes potentially involved in regulating goatskin quality. These discoveries shined new light on deeper understanding the mechanisms of hide-related traits in goat and other livestock.

Effects of lncRNA MTC on protein expression in skin fibroblasts of Liaoning Cashmere goat based on iTRAQ technique

Existing experiments have found a new intergenic lncRNA activated by melatonin, which is called lncRNA MTC. However, the regulatory mechanism of lncRNA MTC in Liaoning Cashmere goat skin fibroblasts has not been clarified. Specific knockdown of lncRNA MTC inhibits cell proliferation and increases apoptosis. iTRAQ reagent was used for relative and absolute quantification of proteins in lncRNA MTC-KD and NC groups to evaluate changes in protein expression during dermal fibroblast development following lncRNA MTC deletion. A total of 5931 proteins were found in Liaoning Cashmere goat skin fibroblasts, of which 123 were differentially expressed, including 32 up-regulated proteins and 91 down-regulated proteins. Of the 91 down-regulated proteins, 32 act mainly through related pathways (e.g., cell cycle, mitochondrial function, ribosomal structure, vesicular transport, cytoskeletal components and skin morphogenesis). LncRNA MTC facilitates the proliferation of Liaoning Cashmere goat skin fibroblasts by regulating ITGB5, TlN2, CTSS, POLG, RAP1B, CHAF1A, CDCA8 and other proteins involved in cell proliferation. The results of this study provide some candidate proteins for the in-depth investigation of the molecular mechanism of lncRNA MTC, which facilitates hair growth in cashmere goats and provides more insights into their regulatory networks and biochemical pathways.

Association between the cashmere production performance, milk production performance, and body size traits and polymorphism of COL6A5 and LOC102181374 genes in Liaoning cashmere goats

The purpose of this study was to analyze the relationship between COL6A5 (collagen type VI alpha 5 chain) and LOC102181374 (alcohol dehydrogenase 1) genes and the production performance of Liaoning cashmere goats by single nucleotide polymorphism (SNP). We have searched for SNP loci of COL6A5 and LOC102181374 genes through sequence alignment and PCR experiments, and have used SPSS and SHEsis software to analyze production data. We obtained five SNP loci in total, including three SNP loci (G50985A, G51140T, G51175A) in COL6A5 gene and two SNP loci (A10067G, T10108C) in LOC102181374 gene. The genotypes G50985A (AG), G51140T (GT), G51175A (AA), A10067G (AA), and T10108C (CC) of these loci have certain advantages in improving the production performance of Liaoning cashmere goats. The haplotype combinations that can improve production performance in COL6A5 gene were H1H5:AGGGAG, H4H4:GGGGAA, and H4H4:GGGGAA. H3H3:GGCC and H2H4:AGTT were the dominant combinations in LOC102181374 gene. At G51175A and A10067G loci, we found that H1H2:AAAG and H1H3:AGAA have dominant effects. These results may provide some support for the molecular breeding of production traits in Liaoning cashmere goats.

Basic principles of hair follicle structure, morphogenesis, and regeneration

Diseases affecting the hair follicle are common in domestic animals, but despite the importance of an intact skin barrier and a fully functional hair coat, knowledge about the detailed morphological features and the diversity of these complex mini-organs are often limited, although mandatory to evaluate skin biopsies with a history of alopecia. The factors that regulate the innate hair follicle formation and the postnatal hair cycle are still not completely understood in rodents, only rudimentarily known in humans, and are poorly understood in our companion animals. This review aims to summarize the current knowledge about hair follicle and hair shaft anatomy, the arrangement of hair follicles, hair follicle morphogenesis in the embryo, and the lifelong regeneration during the postnatal hair cycle in domestic animals. The role of follicular stem cells and the need for a multitude of interacting signaling events during hair follicle morphogenesis and regeneration is unquestioned. Because of the lack of state of the art methods that can be applied in rodents but are not feasible in companion animals, most of the information in this review is based on rodent studies. However, the few data from domestic animals that are available will be discussed, and it can be assumed that at least the principal molecular mechanisms are similar in rodents and other species.

Hedgehog signaling is required for endometrial remodeling and myometrial homeostasis in the cycling mouse uterus

Decades of work demonstrate that the mammalian estrous cycle is controlled by cycling steroid hormones. However, the signaling mechanisms that act downstream, linking hormonal action to the physical remodeling of the cycling uterus, remain unclear. To address this issue, we analyzed gene expression at all stages of the mouse estrous cycle. Strikingly, we found that several genetic programs well-known to control tissue morphogenesis in developing embryos displayed cyclical patterns of expression. We find that most of the genetic architectures of Hedgehog signaling (ligands, receptors, effectors, and transcription factors) are transcribed cyclically in the uterus, and that conditional disruption of the Hedgehog receptor smoothened not only elicits a failure of normal cyclical thickening of the endometrial lining but also induces aberrant deformation of the uterine smooth muscle. Together, our data shed light on the mechanisms underlying normal uterine remodeling specifically and cyclical gene expression generally.

Combination of Transcriptomics and Proteomics Reveals Differentially Expressed Genes and Proteins in the Skin of EDAR Gene-Targeted and Wildtype Cashmere Goats

Cashmere goats play a pivotal role in the animal hair industry and are economically valuable. Cashmere is produced through the periodic growth of secondary hair follicles. To improve their yield of cashmere, the regulatory mechanisms of cashmere follicle growth and development need to be analysed. Therefore, in this study, EDAR gene-targeted cashmere goats were used as an animal model to observe the phenotypic characteristics of abnormal hair growth and development at the top of the head. Transcriptomic and proteomic techniques were used to screen for differentially expressed genes and proteins. In total, 732 differentially expressed genes were identified, including 395 upregulated and 337 downregulated genes. In addition, 140 differentially expressed proteins were identified, including 69 upregulated and 71 downregulated proteins. These results provide a research target for elucidating the mechanism through which EDAR regulates hair follicle growth in cashmere goats. It also enriches the available data on the regulatory network involved in hair follicle growth.

Melatonin promotes the development of the secondary hair follicles by regulating circMPP5

BACKGROUND

The quality and yield of cashmere fibre are closely related to the differentiation and development of secondary hair follicles in the skin of cashmere goats. The higher the density of secondary hair follicles, the higher the quality and yield of cashmere from the fleece. Development of secondary hair follicles commences in the embryonic stage of life and is completed 6 months after birth. Preliminary experimental results from our laboratory showed that melatonin (MT) treatment of goat kids after their birth could increase the density of secondary hair follicles and, thus, improve the subsequent yield and quality of cashmere. These changes in the secondary hair follicles resulted from increases in levels of antioxidant and expression of anti-apoptotic protein, and from a reduction in apoptosis. The present study was conducted to explore the molecular mechanism of MT-induced secondary hair follicle differentiation and development by using whole-genome analysis.

RESULTS

MT had no adverse effect on the growth performance of cashmere kids but significantly improved the character of the secondary hair follicles and the quality of cashmere, and this dominant effect continued to the second year. Melatonin promotes the proliferation of secondary hair follicle cells at an early age. The formation of secondary hair follicles in the MT group was earlier than that in the control group in the second year. The genome-wide data results involved KEGG analysis of 1044 DEmRNAs, 91 DElncRNAs, 1054 DEcircRNAs, and 61 DEmiRNAs which revealed that the mitogen-activated protein kinase (MAPK) signaling pathway is involved in the development of secondary hair follicles, with key genes (FGF2, FGF21, FGFR3, MAPK3 (ERK1)) being up-regulated and expressed. We also found that the circMPP5 could sponged miR-211 and regulate the expression of MAPK3.

CONCLUSIONS

We conclude that MT achieves its effects by regulating the MAPK pathway through the circMPP5 sponged the miR-211, regulating the expression of MAPK3, to induce the differentiation and proliferation of secondary hair follicle cells. In addition there is up-regulation of expression of the anti-apoptotic protein causing reduced apoptosis of hair follicle cells. Collectively, these events increase the numbers of secondary hair follicles, thus improving the production of cashmere from these goats.

Characterization of Growth Secondary Hair in Min Pig Activated by Follicle Stem Cell Stimulated by Wnt and BMP Signaling Pathway

In China, the national-level protected pig, the Min pig, is characterized by the development of secondary hairs and hair follicles in winter. Factors that dominate the genotype in the growth of secondary hairs are not clear through the concrete cell signaling pathways. This study compared hair phenotypes based on morphological structure, transcriptomics, and potential targeting molecules in the breeds of Min, Berkshire, and Yorkshire pigs. The results indicated that Min pigs have specific characteristics for the growth of secondary hairs compared with the Berkshire and Yorkshire pigs. The transcriptome analyses and quantitative reverse transcription-polymerase chain reaction results revealed that secondary hair growth was activated by follicle stem cells. The specific inhibitors of Wnt and BMP were studied using respective signals. The density of follicles, activity of follicle stem cells, and relative gene expression results have shown that Wnt and BMP stimulate the activity of follicle stem cells, and the Wnt signaling molecule has a significantly better effect than the BMP signaling molecule on stem cells. Wnt and BMP can promote the growth of local secondary hair and gene expression. Therefore, this study was conducted to verify the development mechanisms of secondary hairs, which have potential applications in laboratory animals and comparative medicine.

Identification of key genes and signaling pathways related to Hetian sheep wool density by RNA-seq technology

Hetian sheep is a breed of sheep unique to the Hetian area of Xinjiang whose wool is used for producing blankets. Individual differences and hair follicle density are the key factors affecting wool production. Therefore, this study aimed to assess the Hetian sheep having different wool densities to statistically analyze the wool traits and hair follicle parameters. Furthermore, the transcriptome sequencing analysis was performed on the skins with different wool densities. The results showed that wool quantity and total hair follicle density of the high wool density sheep was significantly higher than low wool density sheep. The sheepskin with high wool density was found to grow more and finer wool than sheepskin with low wool density. A total of 1,452 differentially expressed genes were screened from the two sets of samples, including 754 upregulated and 698 downregulated genes. The differentially expressed genes were involved in the TGF-β/BMP and MAPK signaling pathways related to hair growth. Eleven differentially expressed genes belonging to the KAPs and KIFs might affect the fineness of the wool. The key genes, like the TNF, MAP2K2, INHBA, FST, PTPN11, MAP3K7, KIT, and BMPR1A, were found to probably affect the growth and density of the wool. The qPCR verified eight genes related to the MAPK pathway whose gene expression trends were consistent with the transcriptome sequencing results. This study furnishes valuable resources for enhancing the quality and production of wool in the Hetian sheep.

Integrated analysis of lncRNAs and mRNAs by RNA-Seq in secondary hair follicle development and cycling (anagen, catagen and telogen) of Jiangnan cashmere goat (Capra hircus)

BACKGROUND

Among the world's finest natural fiber composites is derived from the secondary hair follicles (SHFs) of cashmere goats yield one of the world's best natural fibres. Their development and cycling are characterized by photoperiodism with diverse, well-orchestrated stimulatory and inhibitory signals. Long non-coding RNA (lncRNAs) and mRNAs play important roles in hair follicle (HF) development. However, not many studies have explored their specific functions in cashmere development and cycling. This study detected mRNAs and lncRNAs with their candidate genes and related pathways in SHF development and cycling of cashmere goat. We utilized RNA sequencing (RNA-Seq) and bioinformatics analysis on lncRNA and mRNA expressions in goat hair follicles to discover candidate genes and metabolic pathways that could affect development and cycling (anagen, catagen, and telogen).

RESULTS

We identified 228 differentially expressed (DE) mRNAs and 256 DE lncRNA. For mRNAs, catagen and anagen had 16 upregulated and 35 downregulated DEGs, catagen and telogen had 18 upregulated and 9 downregulated DEGs and telogen and anagen had 52 upregulated and 98 downregulated DEGs. LncRNA witnessed 22 upregulated and 39 downregulated DEGs for catagen and anagen, 36 upregulated and 29 downregulated DEGs for catagen and telogen as well as 66 upregulated and 97 downregulated DEGs for telogen and anagen. Several key genes, including MSTRG.5451.2, MSTRG.45465.3, MSTRG.11609.2, CHST1, SH3BP4, CDKN1A, GAREM1, GSK-3β, DEFB103A KRTAP9-2, YAP1, S100A7A, FA2H, LOC102190037, LOC102179090, LOC102173866, KRT2, KRT39, FAM167A, FAT4 and EGFL6 were shown to be potentially important in hair follicle development and cycling. They were related to, WNT/β-catenin, mTORC1, ERK/MAPK, Hedgehog, TGFβ, NFkB/p38MAPK, caspase-1, and interleukin (IL)-1a signaling pathways.

CONCLUSION

This work adds to existing understanding of the regulation of HF development and cycling in cashmere goats via lncRNAs and mRNAs. It also serves as theoretical foundation for future SHF research in cashmere goats.

Skin transcriptome analysis identifies the key genes underlying fur development in Chinese Tan sheep in the birth and Er-mao periods

Hair follicle development in Tan sheep differs significantly between the birth and Er-mao periods, but the underlying molecular mechanism is still unclear. We profiled the skin transcriptomes of Tan sheep in the birth and Er-mao periods via RNA-seq technology. The Tan sheep examined consisted of three sheep in the birth period and three sheep in the Er-mao period. A total of 364 differentially expressed genes (DEGs) in the skin of Tan sheep between the birth period and the Er-mao period were identified, among which 168 were upregulated and 196 were downregulated. Interestingly, the FOS proto-oncogene (FOS) (fold change = 22.67, P value = 2.15*10^-44) was the most significantly differentially expressed gene. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the FOS gene was significantly enriched in the signaling pathway related to hair follicle development. Immunohistochemical analysis showed that the FOS gene was expressed in the skin of Chinese Tan sheep at the birth and Er-mao periods, with significantly higher expression in the Er-mao period. Our findings suggest that the FOS gene promotes hair follicle development in Tan sheep.

The fiber diameter traits of Tibetan cashmere goats are governed by the inherent differences in stress, hypoxic, and metabolic adaptations: an integrative study of proteome and transcriptome

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08422-x.

Melatonin Regulates the Periodic Growth of Cashmere by Upregulating the Expression of Wnt10b and β-catenin in Inner Mongolia Cashmere Goats

Secondary hair follicle growth in cashmere goats has seasonal cycle changes, and melatonin (MT) has a regulatory effect on the cashmere growth cycle. In this study, the growth length of cashmere was measured by implanting MT in live cashmere goats. The results indicated that the continuous implantation of MT promoted cashmere to enter the anagen 2 months earlier and induce secondary hair follicle development. HE staining of skin tissues showed that the number of secondary hair follicles in the MT-implanted goats was significantly higher than that in the control goats (P < 0.05). Transcriptome sequencing of the skin tissue of cashmere goats was used to identify differentially expressed genes: 532 in February, 641 in October, and 305 in December. Fluorescence quantitative PCR and Western blotting results showed that MT had a significant effect on the expression of Wnt10b, β-catenin, and proteins in the skin tissue of Inner Mongolia cashmere goats. This finding suggested that MT alters the cycle of secondary hair follicle development by changing the expression of related genes. This research lays the foundation for further study on the mechanism by which MT regulates cashmere growth.

Integrated analysis of lncRNA and mRNA reveals novel insights into cashmere fineness in Tibetan cashmere goats

Tibetan cashmere goats are famous for producing the finest, softest and lightest cashmere fiber in China. The growth and development of skin are closely related to fineness and are the key factors affecting the quality of cashmere. To investigate the specific role of long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) in regulating cashmere fineness of Tibetan Cashmere goats in the anagen phase, we conducted high-throughput RNA sequencing of fine-type and coarse-type skin tissues. We identified 2,059 lncRNA candidates (1,589 lncRNAs annotated, 470 lncRNAs novel), and 80 differentially expressed (DE) lncRNAs and their potential targets were predicted. We also identified 384 DE messenger RNAs (mRNAs) out of 29,119 mRNAs. Several key genes in KRT26, KRT28, KRT39, IFT88, JAK3, NOTCH2 and NOTCH3 and a series of lncRNAs, including ENSCHIT00000009853, MSTRG.16794.17, MSTRG.17532.2, were shown to be potentially important for regulating cashmere fineness. GO and KEGG enrichment analyses of DE mRNAs and DE lncRNAs targets significantly enriched in positive regulation of the canonical Wnt signaling pathway, regulation of protein processing and metabolism processes. The mRNA-mRNA and lncRNA-mRNA regulatory networks further revealed potential transcripts involved in cashmere fineness. We further validated the expression patterns of DE mRNAs and DE lncRNAs by quantitative real-time PCR (qRT-PCR), and the results were consistent with the sequencing data. This study will shed new light on selective cashmere goat breeding, and these lncRNAs and mRNAs that were found to be enriched in Capra hircus RNA database.

Transcriptome-based selection and validation of optimal house-keeping genes for skin research in goats (Capra hircus)

BACKGROUND

In quantitative real-time polymerase chain reaction (qRT-PCR) experiments, accurate and reliable target gene expression results are dependent on optimal amplification of house-keeping genes (HKGs). RNA-seq technology offers a novel approach to detect new HKGs with improved stability. Goat (Capra hircus) is an economically important livestock species and plays an indispensable role in the world animal fiber and meat industry. Unfortunately, uniform and reliable HKGs for skin research have not been identified in goat. Therefore, this study seeks to identify a set of stable HKGs for the skin tissue of C. hircus using high-throughput sequencing technology.

RESULTS

Based on the transcriptome dataset of 39 goat skin tissue samples, 8 genes (SRP68, NCBP3, RRAGA, EIF4H, CTBP2, PTPRA, CNBP, and EEF2) with relatively stable expression levels were identified and selected as new candidate HKGs. Commonly used HKGs including SDHA and YWHAZ from a previous study, and 2 conventional genes (ACTB and GAPDH) were also examined. Four different experimental variables: (1) different development stages, (2) hair follicle cycle stages, (3) breeds, and (4) sampling sites were used for determination and validation. Four algorithms (geNorm, NormFinder, BestKeeper, and ΔCt method) and a comprehensive algorithm (ComprFinder, developed in-house) were used to assess the stability of each HKG. It was shown that NCBP3 + SDHA + PTPRA were more stably expressed than previously used genes in all conditions analysis, and that this combination was effective at normalizing target gene expression. Moreover, a new algorithm for comprehensive analysis, ComprFinder, was developed and released.

CONCLUSION

This study presents the first list of candidate HKGs for C. hircus skin tissues based on an RNA-seq dataset. We propose that the NCBP3 + SDHA + PTPRA combination could be regarded as a triplet set of HKGs in skin molecular biology experiments in C. hircus and other closely related species. In addition, we also encourage researchers who perform candidate HKG evaluations and who require comprehensive analysis to adopt our new algorithm, ComprFinder.

Transcriptomic Analysis Reveals the Involvement of lncRNA-miRNA-mRNA Networks in Hair Follicle Induction in Aohan Fine Wool Sheep Skin

Long non-coding RNAs (lncRNA) and microRNAs (miRNA) are new found classes of non-coding RNAs (ncRNAs) that are not translated into proteins but regulate various cellular and biological processes. In this study, we conducted a transcriptomic analysis of ncRNA and mRNA expression in Aohan fine wool sheep (AFWS) at different growth stages (embryonic day 90, embryonic day 120, and the day of birth), and explored their relationship with wool follicle growth. In total, 461 lncRNAs, 106 miRNAs, and 1,009 mRNAs were found to be differentially expressed during the three stages of wool follicle development. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to clarify the roles of the differentially expressed lncRNA, miRNA, and mRNA in the different stages of wool follicle development. Quantitative real-time PCR (qRT-PCR) was used to validate the results of RNA-seq analysis. lncRNA (MSTRG.223165) was found to act as a competing endogenous RNA (ceRNA) and may participate in wool follicle development by acting as an miR-21 sponge. Network prediction implicated the MSTRG.223165-miR-21-SOX6 axis in the wool follicle development. The targeting relationships of miR-21 with SOX6 and MSTRG.223165 were validated in dual-luciferase assays. This is the first report indicating the association of the lncRNA-miRNA-mRNA network with wool follicle development in AFWS. This study provides new insights into the regulation of the wool follicle growth and represents a solid foundation for wool sheep breeding programs.

Screening the key genes of hair follicle growth cycle in Inner Mongolian Cashmere goat based on RNA sequencing

Inner Mongolian Cashmere goat is an excellent local breed selected for the dual-purpose of cashmere and meat. There are three lines of Inner Mongolian Cashmere goat: Erlangshan, Alashan and Aerbasi. Cashmere is a kind of precious textile raw material with a high price. Cashmere is derived from secondary hair follicle (SHF), while hair is derived from primary hair follicle (PHF). The growth cycle of SHF of cashmere goat is 1 year, and it can be divided into three different stages: anagen, catagen and telogen. In this study, we tried to find some important influence factors of SHF growth cycle in skin tissue from Inner Mongolian Cashmere goats by RNA sequencing (RNA-Seq). Three female Aerbasi Inner Mongolian Cashmere goats (2 years old) were used as experimental samples in this study. Skin samples were collected in September (anagen), December (catagen) and March (telogen) at dorsal side from cashmere goats. Results showed that over 511 396 044 raw reads and 487 729 890 clean reads were obtained from sequence data. In total, 51 different expression genes (DEGs) including 29 downregulated genes and 22 upregulated genes were enriched in anagen–catagen comparing group. The 443 DEGs contained 117 downregulated genes and 326 upregulated genes that were enriched in catagen–telogen comparing group. In telogen–anagen comparing group, 779 DEGs were enriched including 582 downregulated genes and 197 upregulated genes. The result of gene ontology (GO) annotation showed that DEGs are in different growth cycle periods, and enriched GO items are mostly related to the transformation of cell and protein. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment result indicated that metabolic process has a great impact on SHF growth cycle. Based on the results of a comprehensive analysis of differentially expressed genes, GO enrichment and KEGG enrichment, we found that FGF5, FGFR1 and RRAS had an effect on the hair follicle growth cycle. The results of this study may provide a theoretical basis for further research on the growth and development of SHF in Inner Mongolian Cashmere goats.

Preliminary study on microR-148a and microR-10a in dermal papilla cells of Hu sheep

Background

Hu sheep, a unique Chinese breed with high reproductive performance, are also well known for their rare white lambskin in China. The quality of lambskin is affected by hair follicles, and dermal papilla cells are an important component of hair follicles that plays a key role in hair follicle growth and development. This study helps elucidate the effect of miR-148a and miR-10a on hair follicle growth and development.

Results

Based on the results of gene chip and high-throughput sequencing, bone morphogenetic protein 7 (BMP7) was used as a research object. Bioinformatics analysis and the dual-luciferase reporter system indicated that, along with Western blot and quantitative real-time polymerase chain reaction (qRT-PCR) that miR-148a and miR-10a target relationships with BMP7. BMP7 was the target gene both for miR-148a and miR-10a by the dual-luciferase reporter system and Western blot. Hu sheep dermal papilla cells were successfully isolated and purified, and after transfecting miR-148a/miR-10a mimics and inhibitors into dermal papilla cells, a Cell Counting Kit-8 (CCK-8) was used to determine that miR-148a/miR-10a inhibited the proliferation of Hu sheep dermal papilla cells. In addition, after the overexpression of miR-148a, the expression levels of Smad3 (P < 0.05), Smad6 (P < 0.05), Smad4 (P < 0.01), and Smad5 (P < 0.01) were significantly higher than those of the control groups. After the inhibition of miR-148a, the expression levels of Smad3 (P < 0.05), Smad4 (P < 0.05), and TGF-β (P < 0.01) were significantly lower than those of the control groups. After the overexpression of miR-10a, the expression levels of Smad1 (P < 0.01), Smad2 (P < 0.05), Smad4 (P < 0.01), Smad5 (P < 0.01), and TGF-β (P < 0.05) were significantly lower than those of the control groups. After the inhibition of miR-10a, the expression levels of Smad1 (P < 0.01) and Smad2 (P < 0.05) were significantly lower than those of the control groups.

Conclusions

These results revealed the target relationship between miR-148a, miR-10a and BMP7, and the effect of miR-148a and miR-10a on the proliferation of dermal papilla cells. They will provide the basis for a follow-up study on how miR-148a, and miR-10a mediate BMP7 regulation of hair follicle growth and development.