The Tibetan cashmere goat (Capra hircus), one of the most ancient breeds in China, has historically been a critical source of meat and cashmere production for local farmers. To adapt to the high-altitude area, extremely harsh climate, and hypoxic environment that the Tibetan cashmere goat lives in, this goat has developed distinct phenotypic traits compared to lowland breeds. However, the genetic components underlying this phenotypic adaptation remain largely unknown. We obtained 118,700 autosomal SNPs through exome sequencing of 330 cashmere goats located at a wide geographic range, including the Tibetan Plateau and low-altitude regions in China. The great majority of SNPs showed low genetic differentiation among populations; however, approximately 2-3 % of the loci showed more genetic differentiation than expected under a selectively neutral model. Together with a combined analysis of high- and low-altitude breeds, we revealed 339 genes potentially under high-altitude selection. Genes associated with cardiovascular system development were significantly enriched in our study. Among these genes, the most evident one was endothelial PAS domain protein 1 (EPAS1), which has been previously reported to be involved in complex oxygen sensing and significantly associated with high-altitude adaptation of human, dog, and grey wolf. The missense mutation Q579L that we identified in EPAS1, which occurs next to the Hypoxia-Inducible Factor-1 (HIF-1) domain, was exclusively enriched in the high-altitude populations. Our study provides insights concerning the population variation in six different cashmere goat populations in China. The variants in cardiovascular system-related genes may explain the observed phenotypic adaptation of the Tibetan cashmere goat.

The growth of cashmere exhibits a seasonal pattern arising from photoperiod change. However, the underlying molecular mechanism remains unclear. We profiled the skin transcriptome of six goats at seven time points during hair follicle cycling via RNA-seq. The six goats comprised three goats exposed to a natural photoperiod and three exposed to a shortened photoperiod. During hair cycle transition, 1713 genes showed differential expression, and 332 genes showed a pattern of periodic expression. Moreover, a short photoperiod induced the hair follicle to enter anagen early, and 246 genes overlapped with the periodic genes. Among these key genes, cold-shock domain containing C2 (CSDC2) was highly expressed in the epidermis and dermis of Cashmere goat skin, although its function in hair-follicle development remains unknown. CSDC2 silencing in mouse fibroblasts resulted in the decreased mRNA expression of two key hair-follicle factors, leading to reduced cell numbers and a lower cell density. Cashmere growth or molting might be controlled by a set of periodic regulatory genes. The appropriate management of short light exposure can induce hair follicles to enter full anagen early through the activation of these regulators. The CSDC2 gene is a potentially important transcription factor in the hair growth cycle.

Long intergenic noncoding RNAs (lincRNAs) have been recognized in recent years as key regulators of biological processes. However, lincRNAs in goat remain poorly characterized both across various tissues and during different developmental stages in goat (Capra hircus). We performed the genome-wide discovery of the lincRNAs in goat by combining the RNA-seq dataset that were generated from 28 cashmere goat skin samples and the 12 datasets of goat tissues downloaded from the NCBI database. We identified a total of 5546 potential lincRNA transcripts that overlapped 3641 lincRNA genes. These lincRNAs exhibited a tissue-specific pattern. Specifically, there are 584 lincRNAs expressed exclusively in only one tissue, and 91 were highly expressed in hair follicle (HF). In addition, 2350 protein-coding genes and 492 lincRNAs were differentially expressed in the skin of goat. The majority exhibited the remarkable differential expression during the transition of the goat skin from the May–June to August–October time point, which covered the different seasons. Fundamental biological processes, such as skin development, were significantly enriched in these genes. Furthermore, we identified several lincRNAs highly expressed in the HF, which exhibited not only the co-expression pattern with the key factors to the HF development but also the activated expression in the August to October time point. Intriguingly, one of spatiotemporal lincRNAs, linc-chig1598 could be a potential regulator of distal-less homeobox 3 expression during the secondary hair follicle growth. This study will facilitate future studies aimed at unravelling the function of lincRNAs in hair follicle development.

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs’ expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17β4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

The Yili goose is the only indigenous goose breed that originates from Anser anser in China, known for its adaptability, strong flying ability, and tender meat with a low body lipid content. The liver plays a crucial role in lipid and glucose metabolism, including the intake, secretion, transportation, and storage of fatty acids (FAs). In this study, RNA-sequencing (RNA-seq) technology was performed to analyze the liver differentially expressed genes of Yili geese and their hybrid geese to investigate differences in liver lipid and glucose metabolism. A total of 452 differentially expressed genes (Q-value < 0.05) were identified. Notably, in KEGG enrichment analysis, four pathways (Q-value < 0.05) were enriched to be associated with lipid and glucose metabolism, including the metabolic pathway, PI3K–Akt signaling pathway, glycolysis/gluconeogenesis, and steroid biosynthesis. This study provides insights into potential candidate genes and metabolic pathways that affect the liver lipid metabolism of Yili goose. These findings provide a better understanding of animal liver lipid deposition and metabolism.