Genome-wide association analysis of fleece traits in Northwest Xizang white cashmere goat.

Northwest Xizang White Cashmere Goat (NXWCG) is the first new breed of cashmere goat in the Xizang Autonomous Region. It has significant characteristics of extremely high fineness, gloss, and softness. Genome-wide association analysis is an effective biological method used to measure the consistency and correlation of genotype changes between two molecular markers in the genome. In addition, it can screen out the key genes affecting the complex traits of biological individuals. The aim of this study was to analyze the genetic mechanism of cashmere trait variation in NXWCG and to discover SNP locus and key genes closely related to traits such as superfine cashmere. Additionally, the key genes near the obtained significant SNPs were analyzed by gene function annotation and biological function mining. In this study, the phenotype data of the four traits (cashmere length, fiber length, cashmere diameter, and cashmere production) were collected. GGP_Goat_70K SNP chip was used for genotyping the ear tissue DNA of the experimental group. Subsequently, the association of phenotype data and genotype data was performed using Gemma-0.98.1 software. A linear mixed model was used for the association study. The results showed that four fleece traits were associated with 18 significant SNPs at the genome level and 232 SNPs at the chromosome level, through gene annotated from Capra hircus genome using assembly ARS1. A total of 107 candidate genes related to fleece traits were obtained. Combined with Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis, we can find that CLNS1A, CCSER1, RPS6KC1, PRLR, KCNRG, KCNK9, and CLYBL can be used as important candidate genes for fleece traits of NXWCG. We used Sanger sequencing and suitability chi-square test to further verify the significant loci and candidate genes screened by GWAS, and the results show that the base mutations loci on the five candidate genes, CCSER1 (snp12579, 34,449,796, A → G), RPS6KC1 (snp41503, 69,173,527, A → G), KCNRG (snp41082, 67,134,820, G → A), KCNK9 (14:78472665, 78,472,665, G → A), and CLYBL (12: 9705753, 9,705,753, C → T), significantly affect the fleece traits of NXWCG. The results provide a valuable basis for future research and contribute to a better understanding of the genetic structure variation of the goat.

Compendium of 5810 genomes of sheep and goat gut microbiomes provides new insights into the glycan and mucin utilization.

BACKGROUND: Ruminant gut microbiota are critical in ecological adaptation, evolution, and nutrition utilization because it regulates energy metabolism, promotes nutrient absorption, and improves immune function. To study the functional roles of key gut microbiota in sheep and goats, it is essential to construct reference microbial gene catalogs and high-quality microbial genomes database. RESULTS: A total of 320 fecal samples were collected from 21 different sheep and goat breeds, originating from 32 distinct farms. Metagenomic deep sequencing and binning assembly were utilized to construct a comprehensive microbial genome information database for the gut microbiota. We successfully generated the largest reference gene catalogs for gut microbiota in sheep and goats, containing over 162 million and 82 million nonredundant predicted genes, respectively, with 49 million shared nonredundant predicted genes and 1138 shared species. We found that the rearing environment has a greater impact on microbial composition and function than the host's species effect. Through subsequent assembly, we obtained 5810 medium- and high-quality metagenome-assembled genomes (MAGs), out of which 2661 were yet unidentified species. Among these MAGs, we identified 91 bacterial taxa that specifically colonize the sheep gut, which encode polysaccharide utilization loci for glycan and mucin degradation. CONCLUSIONS: By shedding light on the co-symbiotic microbial communities in the gut of small ruminants, our study significantly enhances the understanding of their nutrient degradation and disease susceptibility. Our findings emphasize the vast potential of untapped resources in functional bacterial species within ruminants, further expanding our knowledge of how the ruminant gut microbiota recognizes and processes glycan and mucins. Video Abstract.

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.

Impacts of Population Size and Domestication Process on Genetic Diversity and Genetic Load in Genus Ovis.

In theoretical biology, a prevailing hypothesis posits a profound interconnection between effective population size (Ne), genetic diversity, inbreeding, and genetic load. The domestication and improvement processes are believed to be pivotal in diminishing genetic diversity while elevating levels of inbreeding and increasing genetic load. In this study, we performed a whole genome analysis to quantity genetic diversity, inbreeding, and genetic load across seven wild Ovis species and five domesticated sheep breeds. Our research demonstrates that the genetic load and diversity of species in the genus Ovis have no discernible impact on recent Ne, and three species within the subgenus Pachyceros tend to carry a higher genetic load and lower genetic diversity patterns. The results coincide with these species' dramatic decline in population sizes within the subgenus Pachyceros ~80-250 thousand years ago. European mouflon presented with the lowest Ne, lower genetic diversity, and higher individual inbreeding coefficient but a lower genetic load (missense and LoF). This suggests that the small Ne of European mouflon could reduce harmful mutations compared to other species within the genus Ovis. We showed lower genetic diversity in domesticated sheep than in Asiatic mouflon, but counterintuitive patterns of genetic load, i.e., lower weak genetic load (missense mutation) and no significant difference in strong genetic load (LoF mutation) between domestic sheep and Asiatic mouflon. These findings reveal that the "cost of domestication" during domestication and improvement processes reduced genetic diversity and purified weak genetic load more efficiently than wild species.

Drivers of plateau adaptability in cashmere goats revealed by genomic and transcriptomic analyses.

BACKGROUND: The adaptive evolution of plateau indigenous animals is a current research focus. However, phenotypic adaptation is complex and may involve the interactions between multiple genes or pathways, many of which remain unclear. As a kind of livestock with important economic value, cashmere goat has a high ability of plateau adaptation, which provides us with good materials for studying the molecular regulation mechanism of animal plateau adaptation. RESULTS: In this study, 32 Jiangnan (J) and 32 Tibetan (T) cashmere goats were sequenced at an average of 10. Phylogenetic, population structure, and linkage disequilibrium analyses showed that natural selection or domestication has resulted in obvious differences in genome structure between the two breeds. Subsequently, 553 J vs. T and 608 T vs. J potential selected genes (PSGs) were screened. These PSGs showed potential relationships with various phenotypes, including myocardial development and activity (LOC106502520, ATP2A2, LOC102181869, LOC106502520, MYL2, ISL1, and LOC102181869 genes), pigmentation (MITF and KITLG genes), hair follicles/hair growth (YAP1, POGLUT1, AAK1, HES1, WNT1, PRKAA1, TNKS, WNT5A, VAX2, RSPO4, CSNK1G1, PHLPP2, CHRM2, PDGFRB, PRKAA1, MAP2K1, IRS1, LPAR1, PTEN, PRLR, IBSP, CCNE2, CHAD, ITGB7, TEK, JAK2, and FGF21 genes), and carcinogenesis (UBE2R2, PIGU, DIABLO, NOL4L, STK3, MAP4, ADGRG1, CDC25A, DSG3, LEPR, PRKAA1, IKBKB, and ABCG2 genes). Phenotypic analysis showed that Tibetan cashmere goats has finer cashmere than Jiangnan cashmere goats, which may allow cashmere goats to better adapt to the cold environment in the Tibetan plateau. Meanwhile, KRTs and KAPs expression in Jiangnan cashmere goat skin was significantly lower than in Tibetan cashmere goat. CONCLUSIONS: The mutations in these PSGs maybe closely related to the plateau adaptation ability of cashmere goats. In addition, the expression differences of KRTs and KAPs may directly determine phenotypic differences in cashmere fineness between the two breeds. In conclusion, this study provide a reference for further studying plateau adaptive mechanism in animals and goat breeding.

Gut Microbiota and Acylcarnitine Metabolites Connect the Beneficial Association between Estrogen and Lipid Metabolism Disorders in Ovariectomized Mice.

Decreased estrogen level is one of the main causes of lipid metabolism disorders and coronary heart disease in women after menopause. Exogenous estradiol benzoate is effective to some extent in alleviating lipid metabolism disorders caused by estrogen deficiency. However, the role of gut microbes in the regulation process is not yet appreciated. The objective of this study was to investigate the effects of estradiol benzoate supplementation on lipid metabolism, gut microbiota, and metabolites in ovariectomized (OVX) mice and to reveal the importance of gut microbes and metabolites in the regulation of lipid metabolism disorders. This study found that high doses of estradiol benzoate supplementation effectively attenuated fat accumulation in OVX mice. There was a significant increase in the expression of genes enriched in hepatic cholesterol metabolism and a concomitant decrease in the expression of genes enriched in unsaturated fatty acid metabolism pathways. Further screening of the gut for characteristic metabolites associated with improved lipid metabolism revealed that estradiol benzoate supplementation influenced major subsets of acylcarnitine metabolites. Ovariectomy significantly increased the abundance of characteristic microbes that are significantly negatively associated with acylcarnitine synthesis, such as Lactobacillus and Eubacterium ruminantium group bacteria, while estradiol benzoate supplementation significantly increased the abundance of characteristic microbes that are significantly positively associated with acylcarnitine synthesis, such as Ileibacterium and Bifidobacterium spp. The use of pseudosterile mice with gut microbial deficiency greatly facilitated the synthesis of acylcarnitine due to estradiol benzoate supplementation and also alleviated lipid metabolism disorders to a greater extent in OVX mice. IMPORTANCE Our findings establish a role for gut microbes in the progression of estrogen deficiency-induced lipid metabolism disorders and reveal key target bacteria that may have the potential to regulate acylcarnitine synthesis. These findings suggest a possible route for the use of microbes or acylcarnitine to regulate disorders of lipid metabolism induced by estrogen deficiency.

Polymorphisms and mRNA Expression Levels of IGF-1, FGF5, and KAP 1.4 in Tibetan Cashmere Goats.

The Tibetan cashmere goat is a precious breed in China and its cashmere is widely used in clothing and textiles. The genes IGF-1, FGF5, and KAP 1.4 have been shown to be crucial regulators of cashmere growth. In this study, we examined mRNA expression levels of these three genes and detected IGF-1, FGF5, and KAP 1.4 SNP loci in the Tibetan cashmere goat. After amplification and sequence alignment of the genes IGF-1, FGF5, and KAP 1.4 among 206 Tibetan cashmere goats, two new SNP loci were detected in gene KAP 1.4, while no SNP loci were found in amplified fragments of genes IGF-1 and FGF5. The expression levels of gene IGF-1 in Baingoin and Nyima counties were significantly higher than in other counties (p < 0.05). Moreover, the expression level of gene FGF5 in Gêrzê was significantly higher than in Rutog. The expression levels of mRNA in KAP 1.4 showed significant variation among seven counties. There were no significant differences in mRNA expression levels of IGF-1, FGF5, and KAP 1.4 in Tibetan cashmere goats when analysed by sex. The gene IGF-1 was slightly up-regulated in one to five-year-old cashmere goats, except in those that were 4 years old. The mRNA expression levels of FGF5 in one and two-year-old cashmere goats was lower compared with those in three to five-year-old cashmere goats. KAP 1.4 was up-regulated across one to five-year-old cashmere goats. In this study, SNP detection and mRNA expression analysis of IGF-1, FGF5, and KAP 1.4 genes was able to add data to genetic evolutionary analysis. Further studies should be carried out in SNPs to detect other fragments in genes IGF-1 and FGF5, as well as signal pathways and gene functions in protein levels of genes IGF-1, FGF5, and KAP 1.4 in the Tibetan cashmere goat.

Gut Microbial Succession Patterns and Metabolic Profiling during Pregnancy and Lactation in a Goat Model.

The maternal gut microbiome affects the duration of pregnancy, delivery, and lactation. It also coordinates the stability of maternal metabolism by regulating and modulating inflammatory cytokines and reproductive hormones. This has been shown in several species; however, the situation in ruminants remains a black box. Here, we aimed to elucidate the relationship between the hindgut microbiota, metabolism, and reproductive hormones in domestic goats (Capra hircus) during nonpregnancy, pregnancy, and lactation stages. The hindgut microbiota was altered during these three stages, with a drastic decrease in the abundance of Family_XIII_AD3011_group in the second and third trimesters of pregnancy. Additionally, a decline in the abundance of Christensenellaceae_R-7_group and Turicibacter was observed from the nonpregnancy stage to late gestation. Family_XIII_AD3011_group and Paeniclostridium were strongly correlated with decreased fecal estradiol and progesterone. Furthermore, we generated a metabolome atlas of the gut and serum from nonpregnancy to lactation to reveal the specific metabolic fingerprints of each physiological stage. Several specific gut metabolites, including carnitine C8:1, γ-aminobutyric acid, and indole-3-carboxylic acid, were negatively correlated with the fecal and serum estradiol concentrations. In contrast, 2'-deoxyinosine, deoxyadenosine, and 5'-deoxyadenosine were positively correlated with the fecal and serum estradiol concentrations. The levels of 2'-deoxyinosine, deoxyadenosine, and 5'-deoxyadenosine in fecal samples were positively correlated with Family_XIII_AD3011_group. Other serum metabolites, such as (±)12-HEPE (hydroxy eicosapentaenoic acid), (±)15-HEPE, (±)18-HEPE, cytidine, uracil, and 5-hydroxyindole-3-acetic acid, were negatively correlated with the serum concentrations of estradiol and progesterone. Finally, Corynebacterium and Clostridium_sensu_stricto_1 in the fecal samples were positively correlated with the abundance of 11,12-EET (epoxy-eicosatrienoic acid), (±)18-HEPE, (±)15-HEPE, and (±)12-HEPE in the serum. IMPORTANCE Our findings revealed that the activity of Family_XIII_AD3011_group and Corynebacterium is strongly correlated with the beneficial regulation of physiological hormones and metabolic changes during pregnancy and lactation. These findings are key for guiding targeted microbial therapeutic approaches to modulate microbiomes in gestating and lactating mammals.

Markhor-derived Introgression of a Genomic Region Encompassing PAPSS2 Confers High-altitude Adaptability in Tibetan Goats.

Understanding the genetic mechanism of how animals adapt to extreme conditions is fundamental to determine the relationship between molecular evolution and changing environments. Goat is one of the first domesticated species and has evolved rapidly to adapt to diverse environments, including harsh high-altitude conditions with low temperature and poor oxygen supply but strong ultraviolet radiation. Here, we analyzed 331 genomes of domestic goats and wild caprid species living at varying altitudes (high > 3000 m above sea level and low < 1200 m), along with a reference-guided chromosome-scale assembly (contig-N50: 90.4 Mb) of a female Tibetan goat genome based on PacBio HiFi long reads, to dissect the genetic determinants underlying their adaptation to harsh conditions on the Qinghai-Tibetan Plateau (QTP). Population genomic analyses combined with genome-wide association studies (GWAS) revealed a genomic region harboring the 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2) gene showing strong association with high-altitude adaptability (PGWAS = 3.62 × 10-25) in Tibetan goats. Transcriptomic data from 13 tissues revealed that PAPSS2 was implicated in hypoxia-related pathways in Tibetan goats. We further verified potential functional role of PAPSS2 in response to hypoxia in PAPSS2-deficient cells. Introgression analyses suggested that the PAPSS2 haplotype conferring the high-altitude adaptability in Tibetan goats originated from a recent hybridization between goats and a wild caprid species, the markhor (Capra falconeri). In conclusion, our results uncover a hitherto unknown contribution of PAPSS2 to high-altitude adaptability in Tibetan goats on QTP, following interspecific introgression and natural selection.

Gut microbiota-derived metabolites contribute negatively to hindgut barrier function development at the early weaning goat model.

Early weaning induces intestinal injury, leading to a series of long-term symptoms such as inflammation, malabsorption and diarrhea. In this study, we hypothesized that microbes and their metabolites modulate the host's inflammatory response to early weaning stress in a goat model. A total of 18 female Tibetan goat kids (n = 9) were weaned from their mothers at 28 d (D28) and 60 d (D60) postpartum. D60 and D28 groups were fed the same solid diet ad libitum from weaning to 75 d of age. The colonic epithelium was subject to RNA-sequencing, the caecal digesta metabolomics were assessed by liquid chromatography-tandem mass spectrometry (LC-MS/MS), and the caecal microbiota composition was analysed by 16S ribosomal RNA gene sequencing. We found that early weaning substantially increased the colonic pro-apoptotic gene expression of B-cell lymphoma associated X (Bax), caspase-9, and caspase-3, and decreased the expression of zonula occludens-1 (ZO-1) and claudin-1 (P < 0.01). In addition, a significant Bacteroides acidifaciens enrichment was observed in the hindgut of early-weaned goats (P < 0.01), which negatively correlated with lysophosphatidylcholine products. Similarly, the chemokine signaling, IL-17 signaling, and peroxisome proliferators-activated receptor (PPAR) signaling pathways were upregulated in the colonic mucosa of the early-weaned goats. By applying caecal microbiota transplantation from goats to defaunated C57/6J mice, we confirmed that caecal microbiota of D28 goat kids increased the relative abundance of B. acidifaciens and significantly up-regulated the genes of Bax, G protein-coupled receptor (GPR) 109A, GPR 43, fatty acid binding protein 6, nuclear receptor subfamily 1 group H member 3, angiotensin converting enzyme 2, and IL-6 expression (P < 0.05), and decreased ZO-1, and claudin-1 protein expression in the mice jejunum and colon (P < 0.001). These results proposed that the hindgut microbiota and metabolites mediate the barrier function weakening during early weaning, and the relative abundance of B. acidifaciens was negatively correlated with the hindgut barrier gene expression. This study demonstrates how weaning stress can affect key host-microbe interaction regulators in the hindgut, in a lysophosphatidylcholine dependent and independent manner. Furthermore, based on our mice data, these results are transferable to other mammal species.

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.

Taxonomic and functional adaption of the gastrointestinal microbiome of goats kept at high altitude (4800 m) under intensive or extensive rearing conditions.

The gut microbiota composition is influenced by the diet as well as the environment in both wild and domestic animals. We studied the effects of two feeding systems on the rumen and hindgut microbiome of semi-feral Tibetan goats kept at high altitude (∼4800 m) using 16S rRNA gene and metagenomic sequencing. Intensive drylot feeding resulted in significantly higher zootechnical performance, narrower ruminal acetate: propionate ratios and a drop in the average rumen pH at slaughter to ∼5.04. Hindgut microbial adaption appeared to be more diverse in the drylot group suggesting a higher influx of undegraded complex non-starch polysaccharides from the rumen. Despite their higher fiber levels in the diet, grazing goats exhibited lower counts of Methanobrevibacter and genes associated with the hydrogenotrophic methanogenesis pathway, presumably reflecting the scarce dietary conditions (low energy density) when rearing goats on pasture from extreme alpine environments. These conditions appeared to promote a relevant abundance of bacitracin genes. In parallel, we recognized a significant increase in the abundance of antibiotic resistance genes in the digestive tracts of drylot animals. In summary, this study provides a deeper insight into the metataxonomic and functional adaption of the gastrointestinal microbiome of goats subject to intensive drylot and extensive pasture rearing conditions at high altitude.

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.