Differential gene expression analysis between anagen and telogen of Capra hircus skin based on the de novo assembled transcriptome sequence

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.

Identification and Characterization of Circular RNAs (circRNAs) Using RNA-Seq in Two Breeds of Cashmere Goats

Circular RNA (circRNA) is a type of non-coding RNA generated from back-splicing the reactions of linear RNA. It plays an important role in various cellular and biological processes. However, there are few studies about the regulatory effect of circRNAs on cashmere fiber traits in cashmere goats. In this study, the expression profiles of circRNAs in skin tissue were compared between Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats, with a significant difference in cashmere fiber yield, cashmere fiber diameter, and cashmere fiber color, using RNA-seq. A total of 11,613 circRNAs were expressed in the caprine skin tissue, and their type, chromosomal distribution, and length distribution were characterized. A total of 115 up-regulated circRNAs and 146 down-regulated circRNAs in LC goats were screened compared to ZB goats. The authenticity of 10 differentially expressed circRNAs was validated by detecting their expression levels and the head-to-tail splice junction using RT-PCR and DNA sequencing, respectively. The parent genes of differentially expressed circRNA were mainly enriched in some Gene Ontology (GO) terms and pathways related to cashmere fiber traits, such as the canonical Wnt signaling pathway, which is involved in the regulation of cell promotion, stem cell proliferation, Wnt signaling pathway regulation, epithelial morphogenesis, MAPK signaling pathway, and cell adhesion molecules pathway. Eight differentially expressed circRNAs were further selected to construct a circRNA-miRNA network, and some miRNAs that were previously reported as related to fiber traits were found in the network. This study provides a deep understanding of the roles of circRNAs in the regulation of cashmere fiber traits in cashmere goats and the involvement of differential splicing in phenotypic expression according to breed and region.

Effect of Sox18 on the Induction Ability of Dermal Papilla Cells in Hu Sheep

Sox18 is a developmental gene that encodes transcription factors. It has been indicated as be a key gene affecting the growth and development of hair follicles, in which dermal papilla cells (DPCs) have been demonstrated to play an important role through their ability to induce the formation of hair follicles. Pre-laboratory studies have found that Sox18 is differentially expressed in the dermal papilla cells of different pattern types of Hu sheep. We speculated that Sox18 plays an important role in the dermal papilla cells of Hu sheep. In our study, we analyzed the effect of Sox18 on the induction ability of DPCs in order to elucidate the function and molecular mechanism of Sox18 in the DPCs of Hu sheep. We first identified the expression of Sox18 in the DPCs of Hu sheep by immunofluorescence staining. We then used alkaline phosphatase staining, cell morphology observations and RT-PCR to detect the effect of Sox18 on the induction of DPCs after overexpression of or interference with Sox18. We also used RT-PCR, WB and immunofluorescence staining to detect the effect of Sox18 on the Wnt/β-catenin signal pathway in DPCs. We found that Sox18 was specifically expressed in the DPCs of Hu sheep, and that Sox18 could enhance the alkaline phosphatase activity in the DPCs of Hu sheep and accelerate cell agglutination. The results of RT-PCR revealed that Sox18 promoted the mRNA expression of Versican, HHIP and FGFRI, and inhibited the mRNA expression of BMP4 and WIF1. Further studies showed that Sox18 promoted the expression of β-catenin and activated the Wnt/β-catenin signal pathway in DPCs. When the Wnt/β-catenin signal pathway of DPCs was activated, the induction ability of DPCs was enhanced. Overall, we believe that Sox18 could enhance the induction ability of DPCs in Hu sheep and regulate the induction ability of DPCs through the Wnt/β-catenin signal pathway.

Time-course RNA-seq analysis reveals stage-specific and melatonin-triggered gene expression patterns during the hair follicle growth cycle in Capra hircus

Background

Cashmere goat is famous for its high-quality fibers. The growth of cashmere in secondary hair follicles exhibits a seasonal pattern arising from circannual changes in the natural photoperiod. Although several studies have compared and analyzed the differences in gene expression between different hair follicle growth stages, the selection of samples in these studies relies on research experience or morphological evidence. Distinguishing hair follicle growth cycle according to gene expression patterns may help to explore the regulation mechanisms related to cashmere growth and the effect of melatonin from a molecular level more accurately.

Results

In this study, we applied RNA-sequencing to the hair follicles of three normal and three melatonin-treated Inner Mongolian cashmere goats sampled every month during a whole hair follicle growth cycle. A total of 3559 and 988 genes were subjected as seasonal changing genes (SCGs) in the control and treated groups, respectively. The SCGs in the normal group were divided into three clusters, and their specific expression patterns help to group the hair follicle growth cycle into anagen, catagen and telogen stages. Some canonical pathways such as Wnt, TGF-beta and Hippo signaling pathways were detected as promoting the hair follicle growth, while Cell adhesion molecules (CAMs), Cytokine-cytokine receptor interaction, Jak-STAT, Fc epsilon RI, NOD-like receptor, Rap1, PI3K-Akt, cAMP, NF-kappa B and many immune-related pathways were detected in the catagen and telogen stages. The PI3K-Akt signaling, ECM-receptor interaction and Focal adhesion were found in the transition stage between telogen to anagen, which may serve as candidate biomarkers for telogen-anagen regeneration. A total of 16 signaling pathways, 145 pathway mRNAs, and 93 lncRNAs were enrolled to construct the pathway-mRNA-lncRNA network, which indicated the function of lncRNAs through interacting with their co-expressed mRNAs. Pairwise comparisons between the control and melatonin-treated groups also indicated 941 monthly differentially expressed genes (monthly DEGs). These monthly DEGs were mainly distributed from April and September, which revealed a potential signal pathway map regulating the anagen stage triggered by melatonin. Enrichment analysis showed that Wnt, Hedgehog, ECM, Chemokines and NF-kappa B signaling pathways may be involved in the regulation of non-quiescence and secondary shedding under the influence of melatonin.

Conclusions

Our study decoded the key regulators of the whole hair follicle growth cycle, laying the foundation for the control of hair follicle growth and improvement of cashmere yield.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08331-z.

Genome-wide transcriptome profiling uncovers differential miRNAs and lncRNAs in ovaries of Hu sheep at different developmental stages

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

Comparative proteomics reveals genetic mechanisms underlying secondary hair follicle development in fine wool sheep during the fetal stage

The aim of this study was to investigate the genetic mechanisms underlying wool production by characterizing the skin protein profile and determining the proteomic changes that occur as a consequence of development in wool-producing sheep using a label-free proteomics approach. Samples were collected at four stages during gestation (87, 96, 102, and 138 days), and every two consecutive stages were statistically compared (87 versus 96, 96 versus 102, and 102 versus 138 days). We identified 227 specific proteins in the sheep proteome that were present in all four stages, and 123 differentially abundant proteins (DAPs). We also observed that the microstructure of the secondary follicles changed significantly during the development of the fetal skin hair follicle. The screened DAPs were strictly related to metabolic and skin development pathways, and were associated with pathways such as the glycolysis/gluconeogenesis. These analyses indicated that the wool production of fine wool sheep is regulated via a variety of pathways. These findings provide an important resource that can be used in future studies of the genetic mechanisms underlying wool traits in fine wool sheep, and the identified DAPs should be further investigated as candidate markers for predicting wool traits in sheep. SIGNIFICANCE: Wool quality (fiber diameter, length, etc.) is an important economic trait of fine wool sheep that is determined by secondary follicle differentiation and re-differentiation. Secondary follicles of fine wool sheep developed from a bud (87 days), and underwent differentiation (96 days) and rapid growth (102 days) until maturity (138 days) during gestation. Comparative analysis based on differential proteomics of these four periods could provide a better understanding of the wool growth mechanism of fine wool sheep and offer novel strategies for improving fine wool quality by breeding.

Exploring differentially expressed genes associated with coat color in goat skin using RNA-seq

Fur color in domestic goats is an important, genetically determined characteristic that is associated with economic value. This study was designed to perform a comprehensive expression profiling of genes expressed in the skin tissues from Laiwu Black goat and Lubei White goat. Comparisons of black and white goat skin transcriptomes revealed 102 differentially expressed genes (DEGs), of which 38 were upregulated and 64 downregulated in black skin compared with white skin. Among the DEGs, we identified six genes involved in pigmentation, including agouti signaling protein (ASIP), CAMP responsive element binding protein 3-like 1 (CREB3L1), dopachrome tautomerase (DCT), premelanosome protein (PMEL), transient receptor potential cation channel subfamily M member 1 (TRPM1), and tyrosinase-related protein 1 (TYRP1). Notably, there were no significant differences in the expression of melanocortin 1 receptor, microphthalmia-associated transcription factor, tyrosinase, and KIT proto-oncogene receptor tyrosine kinase between the black and white skin samples, whereas ASIP expression was detected only in white skin. PMEL, TRPM1, TYRP1, and DCT showed higher expression in black goat skin, but ASIP and CREB3L1 had higher expression in white goat skin. Quantitative polymerase chain reaction results for PMEL, TRPM1, DCT, TYRP1, and CREB3L1 expression were consistent with those for RNA-seq. These results will expand our understanding of the complex molecular mechanisms of skin physiology and melanogenesis in goats, and provide a foundation for future studies.

An Integrated Analysis of Cashmere Fineness lncRNAs in Cashmere Goats

Animal growth and development are regulated by long non-coding RNAs (lncRNAs). However, the functions of lncRNAs in regulating cashmere fineness are poorly understood. To identify the key lncRNAs that are related to cashmere fineness in skin, we have collected skin samples of Liaoning cashmere goats (LCG) and Inner Mongolia cashmere goats (MCG) in the anagen phase, and have performed RNA sequencing (RNA-seq) approach on these samples. The high-throughput sequencing and bioinformatics analyses identified 437 novel lncRNAs, including 93 differentially expressed lncRNAs. We also identified 3084 differentially expressed messenger RNAs (mRNAs) out of 27,947 mRNAs. Gene ontology (GO) analyses of lncRNAs and target genes in cis show a predominant enrichment of targets that are related to intermediate filament and intermediate filament cytoskeleton. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, sphingolipid metabolism is a significant pathway for lncRNA targets. In addition, this is the first report to reveal the possible lncRNA–mRNA regulatory network for cashmere fineness in cashmere goats. We also found that lncRNA XLOC_008679 and its target gene, KRT35, may be related to cashmere fineness in the anagen phase. The characterization and expression analyses of lncRNAs will facilitate future studies on the potential value of fiber development in LCG.

Transcriptomic analysis reveals critical genes for the hair follicle of Inner Mongolia cashmere goat from catagen to telogen

There are two main types of hair follicle in Inner Mongolia Cashmere goats, the primary hair follicle (PHF) producing hair fibers and the secondary hair follicle (SHF) producing cashmere fibers. Of both fibers from cashmere-bearing goats in Aerbasi, Inner Mongolia, the timing of cyclical phases for the cashmere have been well clarified but hair fibers have been less noticeable. Herein, we evaluated transcriptome of PHF and SHF from the same three goats in Aerbasi at the catagen- and telogen phase of cashmere growth. We totally found 1977 DEGs between PHFs at the telogen and catagen phases of SHF, 1199 DEGs between telogen- and catagen SHF, 2629 DEGs between PHF at the catagen phase of SHF and catagen SHF, and 755 DEGs between PHF at the telogen phase of SHF and telogen SHF. By analyzing gene functions based on GO and KEGG database, we found that the DEGs have functions in muscle contraction and muscle filament sliding between catagen- and telogen SHF, indicating that arrector pilli muscles might play a role on the transition from catagen to telogen. Moreover, considering that the enriched GO and KEGG categories of the DEGs between PHF and SHF, we suggested that part of PHF might rest in their own anagen phase when SHF are at catagen, but PHF might enter into the telogen phase at SHF’s telogen. Finally, we highly recommended the several potential genes acting as the regulators of the transition between growth phases including IL17RB and eight members of ZNF. These results provide insight into molecular mechanisms on the transition of SHF from catagen to telogen together with PHF’s growth situation at SHF’s catagen and telogen in Inner Mongolia Cashmere goats.

Whole-genome bisulfite sequencing of goat skins identifies signatures associated with hair cycling

Background

Hair follicles (HFs), upon development, undergo repetitive cycles of growth (anagen), regression (catagen), and rest (telogen). The transition between the stages is determined by multiple molecular signals, including DNA methylation, which plays important roles in mammalian cellular identity and is essential for the development of HFs. Secondary hair follicles (SHFs) in cashmere goat exhibit classic cyclic hair development, and little has been done on a genome-wide scale to examine potentially methylated genes involved in the hair cyclic transition.

Results

Genome-wide DNA methylation profiles between skin tissues sampled during the anagen and telogen stages in cashmere goats were investigated using whole-genome bisulfite sequencing (WGBS). The methylation status was observed to be higher in the skin samples with HFs in the telogen than those in the anagen stage. A total of 1311 differentially methylated regions (DMRs) were identified between the two groups, which contained 493 fully annotated DMR-related genes (DMGs) (269 Hyper- DMGs and 224 Hypo-DMGs). Furthermore, a significant over-representation of the functional categories for DMGs related to immune response and intercellular crosstalk during hair cycling was observed. By integrating DNA methylation and mRNA expression data, we revealed that four genes (FMN1, PCOLCE, SPTLC3, and COL5A1) are crucial factors for elucidating epigenetic mechanisms contributing to the telogen-to-anagen transition.

Conclusion

Our study provided systematic methylome maps pertaining to the hair cycling stages (anagen vs telogen) at a single-base resolution, and revealed stage-specific methylation loci during cashmere growth or quiescence. Furthermore, we identified epigenetically regulated genes that are potentially involved in HF development and growth in cashmere goats, and likely in other mammal species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5002-5) contains supplementary material, which is available to authorized users.

Dynamic transcriptome profiling towards understanding the morphogenesis and development of diverse feather in domestic duck

Background

Feathers with complex and fine structure are hallmark avian integument appendages, which have contributed significantly to the survival and breeding for birds. Here, we aimed to explore the differentiation, morphogenesis and development of diverse feathers in the domestic duck.

Results

Transcriptome profiles of skin owing feather follicle from two body parts at three physiological stages were constructed to understand the molecular network and excavate the candidate genes associated with the development of plumulaceous and flight feather structures. The venn analysis of differentially expressed genes (DEGs) between abdomen and wing skin tissues at three developmental stages showed that 38 genes owing identical differentially expression pattern. Together, our data suggest that feather morphological and structural diversity can be possibly related to the homeobox proteins. The key series-clusters, many candidate biological processes and genes were identified for the morphogenesis, growth and development of two feather types. Through comparing the results of developmental transcriptomes from plumulaceous and flight feather, we found that DEGs belonging to the family of WNT, FGF and BMP have certain differences; even the consistent DEGs of skin and feather follicle transcriptomes from abdomen and wing have the different expression patterns.

Conclusions

Overall, this study detected many functional genes and showed differences in the molecular mechanisms of diverse feather developments. The findings in WNT, FGF and BMP, which were consistent with biological experiments, showed more possible complex modulations. A correlative role of HOX genes was also suggested but future biological verification experiments are required. This work provided valuable information for subsequent research on the morphogenesis of feathers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4778-7) contains supplementary material, which is available to authorized users.

Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin

Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-β, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.

A comparison of transcriptomic patterns measured in the skin of Chinese fine and coarse wool sheep breeds

We characterised wool traits, and skin gene expression profiles of fine wool Super Merino (SM) and coarse wool Small Tail Han (STH) sheep. SM sheep had a significantly higher total density of wool follicles, heavier fleeces, finer fibre diameter, and increased crimp frequency, staple length and wool grease (lanolin) production. We found 435 genes were expressed at significantly different levels in the skin of the two breeds (127 genes more highly in SM and 308 genes more highly in STH sheep). Classification of the genes more highly expressed in SM sheep revealed numerous lipid metabolic genes as well as genes encoding keratins, keratin-associated proteins, and wool follicle stem cell markers. In contrast, mammalian epidermal development complex genes and other genes associated with skin cornification and muscle function were more highly expressed in STH sheep. Genes identified in this study may be further evaluated for inclusion in breeding programs, or as targets for therapeutic or genetic interventions, aimed at altering wool quality or yield. Expression of the lipid metabolic genes in the skin of sheep may be used as a novel trait with the potential to alter the content or properties of lanolin or the fleece.

Gene expression profiling analysis reveals fur development in rex rabbits (Oryctolagus cuniculus)

Fur is an important economic trait in rabbits. The identification of genes that influence fur development and knowledge regarding the actions of these genes provides useful tools for improving fur quality. However, the mechanism of fur development is unclear. To obtain candidate genes related to fur development, the transcriptomes of tissues from backs and bellies of Chinchilla rex rabbits were compared. Of the genes analyzed, 336 showed altered expression in the two groups (285 upregulated and 51 downregulated, P ≤ 0.05, fold-change ≥2 or ≤0.5). Using GO and KEGG to obtain gene classes that were differentially enriched, we found several genes to be involved in many important biological processes. In addition, we identified several signaling pathways involved in fur development, including the Wnt and MAPK signaling pathways, revealing mechanisms of skin and hair follicle development, and epidermal cell and keratinocytes differentiation. The obtained rabbit transcriptome and differentially expressed gene profiling data provided comprehensive gene expression information for SFRP2, FRZB, CACNG1, SLC25A4, and SLC16A3. To validate the RNA-seq data, the expression levels of eight differentially expressed genes involved in fur development were confirmed by qRT-PCR. The results of rabbit transcriptomic profiling provide a basis for understanding the molecular mechanisms of fur development.

Genome-wide analysis of DNA methylation in hypothalamus and ovary of Capra hircus

BACKGROUND

DNA methylation is a frequently studied epigenetic modification due to its role in regulating gene expression and hence in biological processes and in determining phenotypic plasticity in organisms. Rudimentary DNA methylation patterns for some livestock species are publically available: among these, goat methylome deserves to be further explored.

RESULTS

Genome-wide DNA methylation maps of the hypothalamus and ovary from Saanen goats were generated using Methyl-CpG binding domain protein sequencing (MBD-seq). Analysis of DNA methylation patterns indicate that the majority of methylation peaks found within genes are located gene body regions, for both organs. Analysis of the distribution of methylated sites per chromosome showed that chromosome X had the lowest number of methylation peaks. The X chromosome has one of the highest percentages of methylated CpG islands in both organs, and approximately 50% of the CpG islands in the goat epigenome are methylated in hypothalamus and ovary. Organ-specific Differentially Methylated Genes (DMGs) were correlated with the expression levels.

CONCLUSIONS

The comparison between transcriptome and methylome in hypothalamus and ovary showed that a higher level of methylation is not accompanied by a higher gene suppression. The genome-wide DNA methylation map for two goat organs produced here is a valuable starting point for studying the involvement of epigenetic modifications in regulating goat reproduction performance.

Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats

Cashmere goat skin contains two types of hair follicles (HF): primary hair follicles (PHF) and secondary hair follicles (SHF). Although multiple genetic determinants associated with HF formation have been identified, the molecules that determine the independent morphogenesis of HF in cashmere goats remain elusive. The growth and development of SHF directly influence the quantity and quality of cashmere production. Here, we report the transcriptome profiling analysis of nine skin samples from cashmere goats using 60- and 120-day-old embryos (E60 and E120, respectively), as well as newborns (NB), through RNA-sequencing (RNA-seq). HF morphological changes indicated that PHF were initiated at E60, with maturation from E120, while differentiation of SHF was identified at E120 until formation of cashmere occurred after birth (NB). The RNA-sequencing analysis generated over 20.6 million clean reads from each mRNA library. The number of differentially expressed genes (DEGs) in E60 vs. E120, E120 vs. NB, and E60 vs. NB were 1,024, 0 and 1,801, respectively, indicating that no significant differences were found at transcriptomic levels between E120 and NB. Key genes including B4GALT4, TNC, a-integrin, and FGFR1, were up-regulated and expressed in HF initiation from E60 to E120, while regulatory genes such as GPRC5D, PAD3, HOXC13, PRR9, VSIG8, LRRC15, LHX2, MSX-2, and FOXN1 were up-regulated and expressed in HF keratinisation and hair shaft differentiation from E120 and NB to E60. Several genes belonging to the KRT and KRTAP gene families were detected throughout the three HF developmental stages. The transcriptional trajectory analyses of all DEGs indicated that immune privilege, glycosaminoglycan biosynthesis, extracellular matrix receptor interaction, and growth factor receptors all played dominant roles in the epithelial-mesenchymal interface and HF formation. We found that the Wnt, transforming growth factor-beta/bone morphogenetic protein, and Notch family members played vital roles in HF differentiation and maturation. The DEGs we found could be attributed to the generation and development of HF, and thus will be critically important for improving the quantity and quality of fleece production in animals for fibres.

A Microarray-Based Analysis Reveals that a Short Photoperiod Promotes Hair Growth in the Arbas Cashmere Goat

Many animals exhibit different behaviors in different seasons. The photoperiod can have effects on migration, breeding, fur growth, and other processes. The cyclic growth of the fur and feathers of some species of mammals and birds, respectively, is stimulated by the photoperiod as a result of hormone-dependent regulation of the nervous system. To further examine this phenomenon, we evaluated the Arbas Cashmere goat (Capra hircus), a species that is often used in this type of research. The goats were exposed to an experimentally controlled short photoperiod to study the regulation of cyclic cashmere growth. Exposure to a short photoperiod extended the anagen phase of the Cashmere goat hair follicle to increase cashmere production. Assessments of tissue sections indicated that the short photoperiod significantly induced cashmere growth. This conclusion was supported by a comparison of the differences in gene expression between the short photoperiod and natural conditions using gene chip technology. Using the gene chip data, we identified genes that showed altered expression under the short photoperiod compared to natural conditions, and these genes were found to be involved in the biological processes of hair follicle growth, structural composition of the hair follicle, and the morphogenesis of the surrounding skin appendages. Knowledge about differences in the expression of these genes as well as their functions and periodic regulation patterns increases our understanding of Cashmere goat hair follicle growth. This study also provides preliminary data that may be useful for the development of an artificial method to improve cashmere production by controlling the light cycle, which has practical significance for livestock breeding.

Expression profiling reveals genes involved in the regulation of wool follicle bulb regression and regeneration in sheep

Wool is an important material in textile manufacturing. In order to investigate the intrinsic factors that regulate wool follicle cycling and wool fiber properties, Illumina sequencing was performed on wool follicle bulb samples from the middle anagen, catagen and late telogen/early anagen phases. In total, 13,898 genes were identified. KRTs and KRTAPs are the most highly expressed gene families in wool follicle bulb. In addition, 438 and 203 genes were identified to be differentially expressed in wool follicle bulb samples from the middle anagen phase compared to the catagen phase and the samples from the catagen phase compared to the late telogen/early anagen phase, respectively. Finally, our data revealed that two groups of genes presenting distinct expression patterns during the phase transformation may have important roles for wool follicle bulb regression and regeneration. In conclusion, our results demonstrated the gene expression patterns in the wool follicle bulb and add new data towards an understanding of the mechanisms involved in wool fiber growth in sheep.

Transcriptome analysis of the Capra hircus ovary

BACKGROUND

Capra hircus is an important economic livestock animal, and therefore, it is necessary to discover transcriptome information about their reproductive performance. In this study, we performed de novo transcriptome sequencing to produce the first transcriptome dataset for the goat ovary using high-throughput sequencing technologies. The result will contribute to research on goat reproductive performance.

METHOD AND RESULTS

RNA-seq analysis generated more than 38.8 million clean paired end (PE) reads, which were assembled into 80,069 unigenes (mean size = 619 bp). Based on sequence similarity searches, 64,824 (60.6%) genes were identified, among which 29,444 and 11,271 unigenes were assigned to Gene Ontology (GO) categories and Clusters of Orthologous Groups (COG), respectively. Searches in the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) showed that 27,766 (63.4%) unigenes were mapped to 258 KEGG pathways. Furthermore, we investigated the transcriptome differences of goat ovaries at two different ages using a tag-based digital gene expression system. We obtained a sequencing depth of over 5.6 million and 5.8 million tags for the two ages and identified a large number of genes associated with reproductive hormones, ovulatory cycle and follicle. Moreover, many antisense transcripts and novel transcripts were found; clusters with similar differential expression patterns, enriched GO terms and metabolic pathways were revealed for the first time with regard to the differentially expressed genes.

CONCLUSIONS

The transcriptome provides invaluable new data for a functional genomic resource and future biological research in Capra hircus, and it is essential for the in-depth study of candidate genes in breeding programs.

Differential gene expression in ovaries of Qira black sheep and Hetian sheep using RNA-Seq technique

The Qira black sheep and the Hetian sheep are two local breeds in the Northwest of China, which are characterized by high-fecundity and low-fecundity breed respectively. The elucidation of mRNA expression profiles in the ovaries among different sheep breeds representing fecundity extremes will helpful for identification and utilization of major prolificacy genes in sheep. In the present study, we performed RNA-seq technology to compare the difference in ovarian mRNA expression profiles between Qira black sheep and Hetian sheep. From the Qira black sheep and the Hetian sheep libraries, we obtained a total of 11,747,582 and 11,879,968 sequencing reads, respectively. After aligning to the reference sequences, the two libraries included 16,763 and 16,814 genes respectively. A total of 1,252 genes were significantly differentially expressed at Hetian sheep compared with Qira black sheep. Eight differentially expressed genes were randomly selected for validation by real-time RT-PCR. This study provides a basic data for future research of the sheep reproduction.

Cyclic expression of Lhx2 is involved in secondary hair follicle development in cashmere goat

Lhx2, a member of LIM homeobox transcription factors, plays a key role in normal tissue development. However, the molecular mechanism of Lhx2 gene in the regulation of the secondary hair follicle cycling in cashmere goat remains largely unknown. In the present study, the Lhx2 gene was cloned and characterized in cashmere goat. The cloned cDNA of Lhx2 was 1233 bp in length, encoding for proteins of 406 amino acids which contained all functionally important domains conserved among vertebrate Lhx2 gene. Tissue distribution analysis showed that Lhx2 mRNA was highly expressed in the skin and low expressed in all other tissues. Immunohistochemical localization revealed that Lhx2 was expressed in secondary hair follicles. Analysis of expression profiles of Lhx2 mRNA during different development stages in secondary hair follicles showed that the highest expression was observed at the anagen stage, while the lowest expression was detected at the telogen stage. The expression tendency during the development stages was that it increased from telogen to anagen, decreased from anagen to catagen, and decreased from catagen to telogen. The expression pattern of Lhx2 protein and mRNA was similar. The mRNA and protein expression of Lhx2 were consistent throughout the development cycle in secondary hair follicles. These findings provided a better understanding of the function of Lhx2 and suggested that the cyclic expression of Lhx2 might play important roles during secondary hair follicle development in cashmere goat.

Expression analysis of KAP9.2 and Hoxc13 genes during different cashmere growth stages by qRT-PCR method

Keratin-associated protein 9.2 (KAP9.2) and Homeobox C13 (Hoxc13) genes were chosen to study because of their biological functions involving hair formation. KAP9.2 gene belongs to the ultra high sulfur KAPs, which is important for hair formation and may have association with cashmere. Hoxc13 takes part in the formation of cashmere keratin and maintaining the normal structure of follicle. It has been reported that Hoxc13 gene exists binding site of KP and KAP genes at its promoter regions in mouse. So the expression of KAP9.2 and Hoxc13 genes was detected at anagen stage vs telogen stage by qRT-PCR. The data showed that KAP9.2 and Hoxc13 gene had similar expression trend at different stages, which indicated that there was interaction between them. KAP9.2 and Hoxc13 gene had lower expression level in anagen than that of in telogen of cashmere growth. In anagen, KAP9.2 and Hoxc13 expressed lower in high cashmere yield individuals than that of in low cashmere yield ones. In telogen, the result was reverse. The study would provide the evidence of involvement of KAP9.2 and Hoxc13 in hair periodic growth.

Transcriptome sequencing reveals differences between anagen and telogen secondary hair follicle-derived dermal papilla cells of the Cashmere goat (Capra hircus)

Dermal papilla is considered the control center of hair follicle growth and hair cycle. The secondary hair follicle (producing cashmere) growth cycle of the Cashmere goat (Capra hircus) is circannual, and each growth phase can be easily distinguished by its long duration. To identify gene expression patterns and differences of the dermal papilla cell (DPC) between the anagen and telogen phases, we established two DPC lines: ana-DPCs (DPCs derived from the anagen secondary hair follicle) and tel-DPCs (DPCs derived from the telogen secondary hair follicle). Compared with the ana-DPCs, the tel-DPCs lost the capacity to form cell aggregates and showed lower cell proliferation rate. Transcriptome sequencing revealed that 825 genes were differentially expressed by at least threefold between the two DPC lines. These genes were significantly enriched in cell cycle control, cell division, and chromosome partitioning from the Eukaryotic Orthologous Groups of proteins (KOG) database and in cell cycle, cell adhesion molecules, cytokine-cytokine receptor interaction, and p53 signaling pathway from the Kyoto Encyclopedia of Gene and Genomes (KEGG) database. Enrichment analyses revealed that in the middle of the telogen the DPCs of secondary hair follicles (SHFs) seemed on the one hand to promote the degeneration of SHFs and cessation of cashmere growth, while on the other hand to resist self-apoptosis and prepare for the regeneration or revivification of fully functional dermal papillae. These findings provide a better understanding of hair follicle growth and will be useful for identification of novel molecules associated with the control of hair growth cycle.

Characterization of Liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis

Background

Liaoning cashmere goat is a famous goat breed for cashmere wool. In order to increase the transcriptome data and accelerate genetic improvement for this breed, we performed denovo transcriptome sequencing to generate the first expressed sequence tag dataset for the Liaoning cashmere goat, using next-generation sequencing technology.

Results

Transcriptome sequencing of Liaoning cashmere goat on a Roche 454 platform yielded 804,601 high-quality reads. Clustering and assembly of these reads produced a non-redundant set of 117,854 unigenes, comprising 13,194 isotigs and 104,660 singletons. Based on similarity searches with known proteins, 17,356 unigenes were assigned to 6,700 GO categories, and the terms were summarized into three main GO categories and 59 sub-categories. 3,548 and 46,778 unigenes had significant similarity to existing sequences in the KEGG and COG databases, respectively. Comparative analysis revealed that 42,254 unigenes were aligned to 17,532 different sequences in NCBI non-redundant nucleotide databases. 97,236 (82.51%) unigenes were mapped to the 30 goat chromosomes. 35,551 (30.17%) unigenes were matched to 11,438 reported goat protein-coding genes. The remaining non-matched unigenes were further compared with cattle and human reference genes, 67 putative new goat genes were discovered. Additionally, 2,781 potential simple sequence repeats were initially identified from all unigenes.

Conclusion

The transcriptome of Liaoning cashmere goat was deep sequenced, denovo assembled, and annotated, providing abundant data to better understand the Liaoning cashmere goat transcriptome. The potential simple sequence repeats provide a material basis for future genetic linkage and quantitative trait loci analyses.

Transcriptome sequencing reveals differences between primary and secondary hair follicle-derived dermal papilla cells of the Cashmere goat (Capra hircus)

The dermal papilla is thought to establish the character and control the size of hair follicles. Inner Mongolia Cashmere goats (Capra hircus) have a double coat comprising the primary and secondary hair follicles, which have dramatically different sizes and textures. The Cashmere goat is rapidly becoming a potent model for hair follicle morphogenesis research. In this study, we established two dermal papilla cell lines during the anagen phase of the hair growth cycle from the primary and secondary hair follicles and clarified the similarities and differences in their morphology and growth characteristics. High-throughput transcriptome sequencing was used to identify gene expression differences between the two dermal papilla cell lines. Many of the differentially expressed genes are involved in vascularization, ECM-receptor interaction and Wnt/β-catenin/Lef1 signaling pathways, which intimately associated with hair follicle morphogenesis. These findings provide valuable information for research on postnatal morphogenesis of hair follicles.