High-throughput sequencing of hair follicle development-related micrornas in cashmere goat at various fetal periods

MiR-21 regulated hair follicle cycle development in Cashmere goats by targeting FGF18 and SMAD7

Increasing Cashmere production can add value because it is the primary product of Cashmere goats. Recent years, peoples find miRNAs are crucial in regulating the development of hair follicle. Following Solexa sequencing, many miRNAs were distinguishingly expressed in telogen skin samples of goats and sheep in earlier study. But the method through which miR-21 controls the growth of hair follicles is still ambiguous. Bioinformatics analysis was used to predict the target genes of miR-21. The mRNA level of miR-21 in telogen Cashmere goat skins was higher than in anagen, according to the results of qRT-PCR, and the target genes expressed similarly with miR-21. Western blot showed similar trend, the protein expression of FGF18 and SMAD7 were lower in anagen samples. The Dual-Luciferase reporter assay confirmed miRNA-21's relationship with its target gene, and the consequences indicated found FGF18 and SMAD7 have positive correlations with miR-21. Western blot and qRT-PCR distinguished the expression of protein and mRNA in miR-21 and its target genes. According to the consequence, we found that target genes expression was increased by miR-21 in HaCaT cells. This study identified that miR-21 might take part in the development of Cashmere goat's hair follicles by targeting FGF18 and SMAD7.

Regulation of cashmere fineness traits by noncoding RNA in Jiangnan cashmere goats

BACKGROUND

Cashmere has long been used as the raw material for wool textiles. The diameter of the cashmere fibre determines its quality and economic value. However, the regulatory role of noncoding RNAs (ncRNAs) in cashmere fineness remains unclear, especially regarding the interaction between ncRNAs and coding RNAs.

RESULTS

Transcriptome sequencing was used to identify the expression profiles of long noncoding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in the skin tissues of Jiangnan cashmere goats with different cashmere fineness levels. Integration analysis of ncRNA and coding RNA was performed in combination with previous research results. The results showed that 16,437 lncRNAs, 2234 circRNAs, and 1322 miRNAs were identified in 8 skin samples of cashmere goats. A total of 403 differentially expressed (DE) lncRNAs, 62 DE circRNAs and 30 DE miRNAs were identified in the skin tissues of the fine groups (Fe) and coarse groups (Ce). We predicted the target gene of DE lncRNA, the target gene of DE miRNA and the host gene of DE circRNA. Based on functional annotation and enrichment analysis of target genes, we found that DE lncRNAs could be involved in regulating the fineness traits of cashmere. The most potential lncRNAs were MSTRG.42054.1, MSTRG.18602.3, and MSTRG.2199.13.

CONCLUSIONS

The data from this study enriched the cashmere goat noncoding RNA database and helped to supplement the annotation of the goat genome. The results provided a new direction for the breeding of cashmere characters.

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.

The Identification and Characteristics of miRNAs Related to Cashmere Fiber Traits in Skin Tissue of Cashmere Goats

microRNAs (miRNAs) are involved in the regulation of biological phenomena by down-regulating the expression of mRNAs. In this study, Liaoning cashmere (LC) goats (n = 6) and Ziwuling black (ZB) goats (n = 6) with different cashmere fiber production performances were selected. We supposed that miRNAs are responsible for the cashmere fiber trait differences. To test the hypothesis, the expression profiles of miRNAs from the skin tissue of the two caprine breeds were compared using small RNA sequencing (RNA-seq). A total of 1293 miRNAs were expressed in the caprine skin samples, including 399 known caprine miRNAs, 691 known species-conserved miRNAs, and 203 novel miRNAs. Compared with ZB goats, 112 up-regulated miRNAs, and 32 down-regulated miRNAs were found in LC goats. The target genes of the differentially expressed miRNAs were remarkably concentrated on some terms and pathways associated with cashmere fiber performance, including binding, cell, cellular protein modification process, and Wnt, Notch, and MAPK signaling pathways. The miRNA-mRNA interaction network found that 14 miRNAs selected may contribute to cashmere fiber traits regulation by targeting functional genes associated with hair follicle activities. The results have reinforced others leading to a solid foundation for further investigation of the influences of individual miRNAs on cashmere fiber traits in cashmere goats.

Integrated analysis of lncRNA and gene expression in longissimus dorsi muscle at two developmental stages of Hainan black goats

It is deemed that meat quality of kids’ is better than that of adults’ for Hainan black goat. Generally, meat quality is affected by many indicators, such as intramuscular fat (IMF) content, muscle fiber diameter and shear force. It is indicated that long non-coding RNAs (lncRNAs) play essential roles in meat quality of goats. However, it is unclear whether and how lncRNAs and genes play their roles in meat quality of Hainan Black goats. Here, we firstly compared the meat quality between two-month-old kids (kids) and adult goats (adults). Then, the lncRNA-seq and RNA-seq data were integrated and analyzed to explore the potential functions of lncRNAs and genes. The results showed that adults’ IMF content and muscle fiber diameter were extremely significantly higher than that of kids (P<0.01). For the sequenced data, average 84,970,398, and 83,691,250 clean reads were obtained respectively for Kids and adults, among which ~96% were mapped to the reference genome of goats. Through analyzing, 18,242 goat annotated genes, 1,429 goat annotated lncRNAs and 2,967 novel lncRNAs were obtained. Analysis of differential expression genes (DEGs) and lncRNAs (DELs) showed that 328 DEGs and 98 DELs existed between kids and adults. Furthermore, functional enrichment analysis revealed that a number of DEGs and DELs were mainly associated with IMF. Primarily, DGAT2 expressed higher in adults than that in kids and CPT1A expressed higher in kids than that in adults. Both of them were overlapped by DEGs and targets of DELs, suggesting the two DEGs and the DELs targeted by the two DEGs might be the potential regulators of goat IMF deposition. Taken together, our results provide basic support for further understanding the function and mechanism of lncRNAs and genes in meat quality of Hainan black goats.

Gene network analysis reveals candidate genes related with the hair follicle development in sheep

Background

Merino sheep are the most famous fine wool sheep in the world. They have high wool production and excellent wool quality and have attracted worldwide attention. The fleece of the Merino sheep is composed predominantly of wool fibers grown from secondary wool follicles. Therefore, it is necessary to study the development of hair follicles to understand the mechanism of wool production. The hair follicle is a complex biological system involved in a dynamic process governed by gene regulation. The hair follicle development process is very complex and poorly understood. The purpose of our research is to identify candidate genes related to hair follicle development, provide a theoretical molecular breeding basis for the cultivation of fine wool sheep, and provide a reference for the problems of hair loss and alopecia areata that affect human beings.

Results

We analyzed mRNAs data in skin tissues of 18 Merino sheep at four embryonic days (E65, E85, E105 and E135) and two postnatal days (P7 and P30). G1 to G6 represent hair follicles developmental at six stages (i.e. E65 to P30). We identified 7879 differentially expressed genes (DEGs) and 12623 novel DEGs, revealed different expression patterns of these DEGs at six stages of hair follicle development, and demonstrated their complex interactions. DEGs with stage-specific expression were significantly enriched in epidermal differentiation and development, hair follicle development and hair follicle morphogenesis and were enriched in many pathways related to hair follicle development. The key genes (LAMA5, WNT10A, KRT25, SOSTDC1, ZDHHC21, FZD1, BMP7, LRP4, TGFβ2, TMEM79, SOX10, ITGB4, KRT14, ITGA6, and GLI2) affecting hair follicle morphogenesis were identified by network analysis.

Conclusion

This study provides a new reference for the molecular basis of hair follicle development and lays a foundation for further improving sheep hair follicle breeding. Candidate genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine wool sheep. These results are a valuable resource for biological investigations of fleece evolution in animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08552-2.

Revisiting the miR-200 Family: A Clan of Five Siblings with Essential Roles in Development and Disease

Over two decades of studies on small noncoding RNA molecules illustrate the significance of microRNAs (miRNAs/miRs) in controlling multiple physiological and pathological functions through post-transcriptional and spatiotemporal gene expression. Among the plethora of miRs that are essential during animal embryonic development, in this review, we elaborate the indispensable role of the miR-200 family (comprising miR-200a, -200b, 200c, -141, and -429) in governing the cellular functions associated with epithelial homeostasis, such as epithelial differentiation and neurogenesis. Additionally, in pathological contexts, miR-200 family members are primarily involved in tumor-suppressive roles, including the reversal of the cancer-associated epithelial–mesenchymal transition dedifferentiation process, and are dysregulated during organ fibrosis. Moreover, recent eminent studies have elucidated the crucial roles of miR-200s in the pathophysiology of multiple neurodegenerative diseases and tissue fibrosis. Lastly, we summarize the key studies that have recognized the potential use of miR-200 members as biomarkers for the diagnosis and prognosis of cancers, elaborating the application of these small biomolecules in aiding early cancer detection and intervention.

Integrated Hair Follicle Profiles of microRNAs and mRNAs to Reveal the Pattern Formation of Hu Sheep Lambskin

Hair follicle development is closely associated with wool curvature. Current studies reveal the crucial role of microRNAs (miRNAs) in hair follicle growth and development. However, few studies are known regarding their role in wool curvature. To reveal the potential roles of miRNAs in Hu sheep lambskin with different patterns, a total of 37 differentially expressed (DE) miRNAs were identified in hair follicles between small waves (SM) and straight wool (ST) groups using RNA-seq. Through functional enrichment and miRNA-mRNA co-expression analysis, some key miRNAs (oar-miR-143, oar-miR-200b, oar-miR-10a, oar-miR-181a, oar-miR-10b, oar-miR-125b, etc.) and miRNA-mRNA pairs (miR-125b target CD34, miR-181a target FGF12, LMO3, miR-200b target ZNF536, etc.) were identified. Though direct or indirect ways affecting hair follicle development, these miRNAs and mRNAs may have possible effects on wool curvature, and this study thus provides valuable insight on potential pattern formation.

Genetic Signatures of Selection for Cashmere Traits in Chinese Goats

Simple Summary

Cashmere goats are a unique husbandry resource in China. These goats are well known for producing the highest cashmere yield and best fiber quality in the world. Although cashmere is highly valued and also known as “fiber gem” and “soft gold”, few studies have examined the genetic basis of cashmere traits in cashmere goats. Here, we identified selection signals by comparing Fst and XP-EHH (the cross population extend haplotype homozygosity test) of a non-cashmere breed (Huanghuai goat) with those of two cashmere breeds (Inner Mongolia and Liaoning cashmere goats). Two genes (WNT10A and CSN3) were potentially associated with cashmere traits. This information may be valuable for studying the genetic uniqueness of cashmere goats and elucidating the mechanisms underlying cashmere traits in cashmere goats.

Abstract

Inner Mongolia and Liaoning cashmere goats in China are well-known for their cashmere quality and yield. Thus, they are great models for identifying genomic regions associated with cashmere traits. Herein, 53 Inner Mongolia cashmere goats, Liaoning cashmere goats and Huanghuai goats were genotyped, and 53,347 single-nucleotide polymorphisms (SNPs) were produced using the Illumina Caprine 50K SNP chip. Additionally, we identified some positively selected SNPs by analyzing Fst and XP-EHH. The top 5% of SNPs had selection signatures. After gene annotation, 222 and 173 candidate genes were identified in Inner Mongolia and Liaoning cashmere goats, respectively. Several genes were related to hair follicle development, such as TRPS1, WDR74, LRRC14, SPTLC3, IGF1R, PADI2, FOXP1, WNT10A and CSN3. Gene enrichment analysis of these cashmere trait-associated genes related 67 enriched signaling pathways that mainly participate in hair follicle development and stem cell pluripotency regulation. Furthermore, we identified 20 overlapping genes that were selected in both cashmere goat breeds. Among these overlapping genes, WNT10A and CSN3, which are associated with hair follicle development, are potentially involved in cashmere production. These findings may improve molecular breeding of cashmere goats in the future.

Effect of Nutritional Restriction on the Hair Follicles Development and Skin Transcriptome of Chinese Merino Sheep

The high concentration of secondary branched wool follicles is a distinctive feature of the Merino sheep. At present, the molecular control of the development and branching of secondary wool follicles (SF) remains elusive. To reveal the potential genes associated with the development of hair follicles, we investigated the characteristics of prenatal and postnatal development of wool follicles, and the transcriptional expression profile in fetuses/lambs from dams under either maternal maintenance or sub-maintenance (75% maintenance) nutrition. The density of SF and the ratio of SF to primary wool follicles (PF) were reduced (p < 0.05) in fetuses from day 105 to 135 of gestation under sub-maintenance nutrition. Differentially expressed genes were enriched in the binding, single-organism process, cellular process, cell and cell part Gene Ontology (GO) functional categories and metabolism, apoptosis, and ribosome pathways. Four candidate genes, SFRP4, PITX1, BAMBI, and KRT16, which were involved in secondary wool follicles branching and development, were identified. Our results indicate that nutritional intervention imposed on pregnant ewes by short-term sub-maintenance nutrition could provide a strategy for the study of wool follicle development. Overall insight into the global gene expression associated with SF development can be used to investigate the underlying mechanisms of SF branching in Merino sheep.

Conservation Assessment of the State Goat Farms by Using SNP Genotyping Data

Conservation of genetic resources is of great concern globally to maintain genetic diversity for sustainable food security. Comprehensive identification of the breed composition, estimation of inbreeding and effective population size are essential for the effective management of farm animal genetic resources and to prevent the animals from genetic erosion. The Zhongwei male (ZWM), Arbas Cashmere male (ACM) and Jining Grey male (JGM) goats are conserved in three different state goat farms in China but their family information, level of inbreeding and effective population size are unknown. We investigated the genomic relationship, inbreeding coefficient and effective population size in these three breeds from three state goat farms using the Illumina goat SNP50 BeadChip. Genomic relationships and phylogenetic analysis revealed that the breeds are clearly separated and formed separate clusters based on their genetic relationship. We obtained a high proportion of informative SNPs, ranging from 91.8% in the Arbas Cashmere male to 96.2% in the Jining Grey male goat breeds with an average mean of 96.8%. Inbreeding, as measured by F_ROH, ranged from 1.79% in ZWM to 8.62% in ACM goat populations. High F_ROH values, elevated genomic coverage of very long ROH (>30 Mb) and severe decline in effective population size were recorded in ACM goat farm. The existence of a high correlation between F_HOM and F_ROH indicates that F_ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records. The N_e estimates 13 generations ago were 166, 69 and 79 for ZWM, ACM and JGM goat farm, respectively indicating that these goat breeds were strongly affected by selection pressure or genetic drift. This study provides insight into the genomic relationship, levels of inbreeding and effective population size in the studied goat populations conserved in the state goat farms which will be valuable in prioritizing populations for conservation and for developing suitable management practices for further genetic improvement of these Chinese male goats.

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.

Genome-Wide Linkage Disequilibrium and the Extent of Effective Population Sizes in Six Chinese Goat Populations Using a 50K Single Nucleotide Polymorphism Panel

Simple Summary

Information on linkage disequilibrium (LD) and the extent of effective population size (Ne) has important implications for exploring the degree of biological diversity, for predicting underlying selection pressure, and for designing animal breeding programs. In this study, we assessed LD, Ne, and the distribution of minor allele frequency in six goat populations. Accordingly, the results of LD and Ne using a single nucleotide polymorphism (SNP) panel (Caprine SNP 50K BeadChip, Lincoln, NE, USA) are helpful for the sustainable conservation, proper management, and utilization of Chinese goat populations.

Abstract

Genome-wide linkage disequilibrium is a useful parameter to study quantitative trait locus (QTL) mapping and genetic selection. In many genomic methodologies, effective population size is an important genetic parameter because of its relationship to the loss of genetic variation, increases in inbreeding, the accumulation of mutations, and the effectiveness of selection. In this study, a total of 193 individuals were genotyped to assess the extent of LD and Ne in six Chinese goat populations using the SNP 50K BeadChip. Across the determined autosomal chromosomes, we found an average of 0.02 and 0.23 for r² and D’ values, respectively. The average r² between all the populations varied little and ranged from 0.055 r² for the Jining Grey to 0.128 r² for the Guangfeng, with an overall mean of 0.083. Across the 29 autosomal chromosomes, minor allele frequency (MAF) was highest on chromosome 1 (0.321) and lowest on chromosome 25 (0.309), with an average MAF of 0.317, and showing the lowest (25.5% for Louping) and highest (28.8% for Qingeda) SNP proportions at MAF values > 0.3. The inbreeding coefficient ranged from 0.064 to 0.085, with a mean of 0.075 for all the autosomes. The Jining Grey and Qingeda populations showed higher Ne estimates, highlighting that these animals could have been influenced by artificial selection. Furthermore, a declining recent Ne was distinguished for the Arbas Cashmere and Guangfeng populations, and their estimated values were closer to 64 and 95, respectively, 13 generations ago, which indicates that these breeds were exposed to strong selection. This study provides an insight into valuable genetic information and will open up the opportunity for further genomic selection analysis of Chinese goat populations.

Genetic diversity and population structure in multiple Chinese goat populations using a SNP panel

Information about genetic diversity and population structure among goat breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of goat breeds. Here, we measured genetic diversity and population structure in multiple Chinese goat populations, namely, Nanjiang, Qinggeda, Arbas Cashmere, Jining Grey, Luoping Yellow and Guangfeng goats. A total of 193 individuals were genotyped for about 47 401 autosomal single nucleotide polymorphisms (SNPs). We found a high proportion of informative SNPs, ranging from 69.5% in the Luoping Yellow to 93.9% in the Jining Grey goat breeds with an average mean of 84.7%. Diversity, as measured by expected heterozygosity, ranged from 0.371 in Luoping Yellow to 0.405 in Jining Grey goat populations. The average estimated pair-wise genetic differentiation (F_ST ) among the populations was 8.6%, ranging from 0.2% to 16% and indicating low to moderate genetic differentiation. Principal component analysis, genetic structure and phylogenetic tree analysis revealed a clustering of six Chinese goat populations according to geographic distribution. The results from this study can contribute valuable genetic information and can properly assist with within-breed diversity, which provides a good opportunity for sustainable utilization of and maintenance of genetic resource improvements in the Chinese goat populations.