Identification of Genes Related to Hair Follicle Cycle Development in Inner Mongolia Cashmere Goat by WGCNA

Stress adaptation in Tibetan cashmere goats is governed by inherent metabolic differences and manifested through variable cashmere phenotypes

Tibetan cashmere goats are not only served as a valuable model for studying adaptation to hypoxia and high-altitude conditions but also playing a pivotal role in bolstering local economies through the provision of premium quality cashmere yarn. In this study, we performed an integration and network analysis of metabolomic, transcriptomic and proteomic to elucidate the role of differentially expressed genes, important metabolites, and relevant cellular and metabolic pathways between the fine (average 12.04 ± 0.03 μm of mean fiber diameter) and coarse cashmere (average 14.88 ± 0.05 μm of mean fber diameter) producing by Tibetan cashmere goats. We identified a distinction of 56 and 71 differential metabolites (DMs) between the F and C cashmere groups under positive and negative ion modes, respectively. The KEGG pathway enrichment analysis of these DMs highlighted numerous pathways predominantly involved in amino acid and protein metabolism, as indicated by the finding that the most impactful pathway was the mammalian target of rapamycin (mTOR) signalling pathway. In the F group, we identified a distinctive metabolic profile where amino acid metabolites including serine, histidine, asparagine, glutamic acid, arginine, valine, aspartic acid, tyrosine, and methionine were upregulated, while lysine, isoleucine, glutamine, tryptophan, and threonine were downregulated. The regulatory network and gene co-expression network revealed crucial genes, metabolites, and metabolic pathways. The integrative omics analysis revealed a high enrichment of several pathways, notably encompassing protein digestion and absorption, sphingolipid signalling, and the synaptic vesicle cycle. Within the sphere of our integrative analysis, DNMT3B was identified as a paramount gene, intricately associated with significant proteins such as HMCN1, CPB2, GNG12, and LRP1. Our present study delineated the molecular underpinnings governing the variations in cashmere characteristics by conducting comprehensive analyses across metabolomic, transcriptomic, and proteomic dimensions. This research provided newly insights into the mechanisms regulating cashmere traits and facilitated the advancement of selective breeding programs aimed at cultivating high-quality superfine Tibetan cashmere goats.

The demographic history and adaptation of Canarian goat breeds to environmental conditions through the use of genome-wide SNP data

BACKGROUND

The presence of goats in the Canary Islands dates back to the late 1st millennium BC, which coincides with the colonization by the Amazigh settlers. However, the exact geographic origin of Canarian goats is uncertain since the Amazigh peoples were distributed over a wide spatial range. Nowadays, three Canarian breeds (Palmera, Majorera and Tinerfeña) are officially recognized, along with two distinct South and North Tinerfeña ecotypes, with the South Tinerfeña and Majorera goats thriving in arid and dry semi-desertic environments and the Palmera and North Tinerfeña goats are adapted to humid and temperate areas that are influenced by trade winds. Genotypes for 224 Canarian goats were generated using the Illumina Goat single nucleotide polymorphism (SNP)50 BeadChip. By merging these data with the genotypes from 1007 individuals of African and Southern European ancestry, our aim was to ascertain the geographic origin of the Canarian goats and identify genes associated with adaptation to diverse environmental conditions.

RESULTS

The diversity indices of the Canarian breeds align with most of those of the analyzed local breeds from Africa and Europe, except for the Palmera goats that showed lower levels of genetic variation. The Canarian breeds demonstrate a significant genetic differentiation compared to other populations, which indicates a history of prolonged geographic isolation. Moreover, the phylogenetic reconstruction indicated that the ancestry of the Canarian goats is fundamentally North African rather than West African. The ADMIXTURE and the TreeMix analyses showed no evidence of gene flow between Canarian goats and other continental breeds. The analysis of runs of homozygosity (ROH) identified 13 ROH islands while the window-based FST method detected 25 genomic regions under selection. Major signals of selection were found on Capra hircus (CHI) chromosomes 6, 7, and 10 using various comparisons and methods.

CONCLUSIONS

This genome-wide analysis sheds new light on the evolutionary history of the four breeds that inhabit the Canary Islands. Our findings suggest a North African origin of the Canarian goats. In addition, within the genomic regions highlighted by the ROH and FST approaches, several genes related to body size and heat tolerance were identified.

Developmental Mapping of Hair Follicles in the Embryonic Stages of Cashmere Goats Using Proteomic and Metabolomic Construction

The hair follicle (HF) is the fundamental unit for fleece and cashmere production in cashmere goats and is crucial in determining cashmere yield and quality. The mechanisms regulating HF development in cashmere goats during the embryonic period remain unclear. Growing evidence suggests that HF development involves complex developmental stages and critical events, and identifying the underlying factors can improve our understanding of HF development. In this study, samples were collected from embryonic day 75 (E75) to E125, the major HF developmental stages. The embryonic HFs of cashmere goats were subjected to proteomic and metabolomic analyses, which revealed dynamic changes in the key factors and signalling pathways controlling HF development at the protein and metabolic levels. Gene ontology and the Kyoto Encyclopaedia of Genes and Genomes were used to functionally annotate 1784 significantly differentially expressed proteins and 454 significantly differentially expressed metabolites enriched in different HF developmental stages. A joint analysis revealed that the oxytocin signalling pathway plays a sustained role in embryonic HF development by activating the MAPK and Ca2+ signalling pathways, and a related regulatory network map was constructed. This study provides a global perspective on the mechanism of HF development in cashmere goats and enriches our understanding of embryonic HF development.

Effect of Fibroblast Growth Factor 10 and an Interacting Non-Coding RNA on Secondary Hair Follicle Dermal Papilla Cells in Cashmere Goats' Follicle Development Assessed by Whole-Transcriptome Sequencing Technology

Cashmere, a keratinised product of secondary hair follicles (SHFs) in cashmere goats, holds an important place in international high-end textiles. However, research on the complex molecular and signal regulation during the development and growth of hair follicles (HFs), which is essential for the development of the cashmere industry, is limited. Moreover, increasing evidence indicates that non-coding RNAs (ncRNAs) participate in HF development. Herein, we systematically investigated a competing endogenous RNA (ceRNA) regulatory network mediated by circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in skin samples of cashmere goat embryos, using whole-transcriptome sequencing technology. We obtained 6468, 394, and 239 significantly differentially expressed mRNAs, circRNAs, and miRNAs, respectively. These identified RNAs were further used to construct a ceRNA regulatory network, mediated by circRNAs, for cashmere goats at a late stage of HF development. Among the molecular species identified, miR-184 and fibroblast growth factor (FGF) 10 exhibited competitive targeted interactions. In secondary HF dermal papilla cells (SHF-DPCs), miR-184 promotes proliferation, inhibits apoptosis, and alters the cell cycle via the competitive release of FGF10. This study reports that FGF10 and its interaction with ncRNAs significantly affect SHF-DPCs, providing a reference for research on the biology of HFs in cashmere goats and other mammals.

Melatonin Promotes the Development of Secondary Hair Follicles in Adult Cashmere Goats by Activating the Keap1-Nrf2 Signaling Pathway and Inhibiting the Inflammatory Transcription Factors NFκB and AP-1

Exogenous melatonin (MT) has been used to promote the growth of secondary hair follicles and improve cashmere fiber quality, but the specific cellular-level mechanisms involved are unclear. This study was carried out to investigate the effect of MT on the development of secondary hair follicles and on cashmere fiber quality in cashmere goats. The results showed that MT improved secondary follicle numbers and function as well as enhanced cashmere fiber quality and yield. The MT-treated goat groups had high secondary-to-primary ratios (S:P) for hair follicles, greater in the elderly group (p < 0.05). Antioxidant capacities of secondary hair follicles improved fiber quality and yield in comparison with control groups (p < 0.05/0.01). Levels of reactive oxygen and nitrogen species (ROS, RNS) and malondialdehyde (MDA) were lowered (p < 0.05/0.01) by MT. There was significant upregulation of antioxidant genes (for SOD-3; GPX-1; NFE2L2) and the protein of nuclear factor (Nrf2), and downregulation of the Keap1 protein. There were significant differences in the expression of genes for secretory senescence-associated phenotype (SASP) cytokines (IL-1β, IL-6, MMP-9, MMP-27, CCL-21, CXCL-12, CXCL-14, TIMP-1,2,3) plus their protein of key transcription factors, nuclear factor kappa B (NFκB) and activator protein-1 (AP-1), in comparison with the controls. We concluded that MT could enhance antioxidant capacity and reduce ROS and RNS levels of secondary hair follicles through the Keap1-Nrf2 pathway in adult cashmere goats. Furthermore, MT reduced the expression of the SASP cytokines genes by inhibiting the protein of NFκB and AP-1 in the secondary hair follicles in older cashmere goats, thus delaying skin aging, improving follicle survival, and increasing the number of secondary hair follicles. Collectively, these effects of exogenous MT enhanced the quality and yield of cashmere fibers, especially at 5-7 years old.

Transcriptome analysis of pika heart tissue reveals mechanisms underlying the adaptation of a keystone species on the roof of the world

High-altitude environments impose intense stresses on living organisms and drive striking phenotypic and genetic adaptations, such as hypoxia resistance, cold tolerance, and increases in metabolic capacity and body mass. As one of the most successful and dominant mammals on the Qinghai-Tibetan Plateau (QHTP), the plateau pika (Ochotona curzoniae) has adapted to the extreme environments of the highest altitudes of this region and exhibits tolerance to cold and hypoxia, in contrast to closely related species that inhabit the peripheral alpine bush or forests. To explore the potential genetic mechanisms underlying the adaptation of O. curzoniae to a high-altitude environment, we sequenced the heart tissue transcriptomes of adult plateau pikas (comparing specimens from sites at two different altitudes) and Gansu pikas (O. cansus). Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were used to identify differentially expressed genes (DEGs) and their primary functions. Key genes and pathways related to high-altitude adaptation were identified. In addition to the biological processes of signal transduction, energy metabolism and material transport, the identified plateau pika genes were mainly enriched in biological pathways such as the negative regulation of smooth muscle cell proliferation, the apoptosis signalling pathway, the cellular response to DNA damage stimulus, and ossification involved in bone maturation and heart development. Our results showed that the plateau pika has adapted to the extreme environments of the QHTP via protection against cardiomyopathy, tissue structure alterations and improvements in the blood circulation system and energy metabolism. These adaptations shed light on how pikas thrive on the roof of the world.

Lef1 and Dlx3 May Facilitate the Maturation of Secondary Hair Follicles in the Skin of Gansu Alpine Merino

Lymphatic enhancer factor 1 (Lef1) and distal-less homeobox 3 (Dlx3) are the transcription factors involved in regulating hair follicle development in mice, goats, and other animals. Their deletion can lead to hair follicle deficiency. In this study, hematoxylin−eosin staining (HE), real-time quantitative PCR (RT-qPCR), immunohistochemistry, and immunofluorescence were used to analyze the expression, location, and biological functions of Lef1 and Dlx3 in the lateral skin of Gansu Alpine Merino aged 1, 30, 60, and 90 days. The results revealed that the number of hair follicles decreased with age and was significantly higher at 1 day than in the other three age groups (p < 0.05). The mRNA levels of Lef1 and Dlx3 in the skin of 30-day old Gansu Alpine Merino were significantly higher than those in the other three age groups (p < 0.05). Protein expression of Lef1 and Dlx3 was lowest at 1 day (p < 0.05) and peaked at 60 days. Lef1 and Dlx3 exhibited a high density and strong positive expression in the dermal papillae; additionally, Dlx3 exhibited a high density and strong positive expression in the inner and outer root sheaths. Collectively, Lef1 and Dlx3 may facilitate the maturation of secondary hair follicles, which is mainly achieved through the dermal papillae and inner and outer root sheaths.