Integrative ATAC-seq and RNA-seq analysis of myogenic differentiation of ovine skeletal muscle satellite cell.

Skeletal muscle satellite cells (SMSCs) play an important role in regulating muscle growth and regeneration. Chromatin accessibility allows physical interactions that synergistically regulate gene expression through enhancers, promoters, insulators, and chromatin binding factors. However, the chromatin accessibility altas and its regulatory role in ovine myoblast differentiation is still unclear. Therefore, ATAC-seq and RNA-seq analysis were performed on ovine SMSCs at the proliferation stage (SCG) and differentiation stage (SCD). 17,460 DARs (differential accessibility regions) and 3732 DEGs (differentially expressed genes) were identified. Based on joint analysis of ATAC-seq and RNA-seq, we revealed that PI3K-Akt, TGF-ß and other signaling pathways regulated SMSCs differentiation. We identified two novel candidate genes, FZD5 and MAP2K6, which may affect the proliferation and differentiation of SMSCs. Our data identify potential cis regulatory elements of ovine SMSCs. This study can provide a reference for exploring the mechanisms of the differentiation and regeneration of SMSCs in the future.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds.

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

Genome-Wide Analysis of Circular RNAs Reveals circCHRNG Regulates Sheep Myoblast Proliferation via miR-133/SRF and MEF2A Axis.

As relatively new members of the non-coding RNA family, circRNAs play important roles in a variety of biological processes. However, the temporal expression pattern and the function of circRNAs during sheep skeletal muscle development remains unclear. This study aimed to identify circRNAs related to sheep skeletal muscle development and explore their roles in myoblast proliferation. The circRNA expression profiles of longissimus dorsi of sheep from F90, L30, and A3Y were obtained by the RNA-seq method. The function and mechanisms of the novel circCHRNG in muscle satellite cell proliferation were explored using CCK-8 assay, Western blot, qPCR, and dual-luciferase reporter assay. We identified 12,375 circRNAs, including 476, 133, and 233 DEcircRNAs found among three comparative groups. KEGG results showed that DEcircRNAs were enriched in muscle contraction, the regulation of cell proliferation, and the AMPK, insulin, and PI3K-Akt signaling pathways. Notably, a novel circRNA, termed circRNA CHRNG, acts as a miR-133 sponge to promote skeletal muscle satellite cell proliferation. Our study provides a systematic description of circRNAs of ovine skeletal muscle across fetal, lamb, and adult stages. GO and KEGG analyses showed that DEcircRNAs were enriched in multiple pathways associated with muscle development, such as the PI3K-Akt and AMPK signaling pathways. In addition, we propose that circCHRNG acts as a miR-133 sponge to upregulate the expression levels of SRF and MEF2A, thereby promoting myoblast proliferation.

EPAS1 Gain-of-Function Mutation Contributes to High-Altitude Adaptation in Tibetan Horses.

High altitude represents some of the most extreme environments worldwide. The genetic changes underlying adaptation to such environments have been recently identified in multiple animals but remain unknown in horses. Here, we sequence the complete genome of 138 domestic horses encompassing a whole altitudinal range across China to uncover the genetic basis for adaptation to high-altitude hypoxia. Our genome data set includes 65 lowland animals across ten Chinese native breeds, 61 horses living at least 3,300 m above sea level across seven locations along Qinghai-Tibetan Plateau, as well as 7 Thoroughbred and 5 Przewalski's horses added for comparison. We find that Tibetan horses do not descend from Przewalski's horses but were most likely introduced from a distinct horse lineage, following the emergence of pastoral nomadism in Northwestern China ∼3,700 years ago. We identify that the endothelial PAS domain protein 1 gene (EPAS1, also HIF2A) shows the strongest signature for positive selection in the Tibetan horse genome. Two missense mutations at this locus appear strongly associated with blood physiological parameters facilitating blood circulation as well as oxygen transportation and consumption in hypoxic conditions. Functional validation through protein mutagenesis shows that these mutations increase EPAS1 stability and its hetero dimerization affinity to ARNT (HIF1B). Our study demonstrates that missense mutations in the EPAS1 gene provided key evolutionary molecular adaptation to Tibetan horses living in high-altitude hypoxic environments. It reveals possible targets for genomic selection programs aimed at increasing hypoxia tolerance in livestock and provides a textbook example of evolutionary convergence across independent mammal lineages.

Expression Profiling and Functional Characterization of miR-26a and miR-130a in Regulating Zhongwei Goat Hair Development via the TGF-ß/SMAD Pathway.

The Zhongwei goat is an important and unique goat breed indigenous to China. It has a natural hair curling phenotype at birth, but the degree of curling gradually decreases with growth. The molecular mechanism underlying the dynamic changes in the wool curvature in Zhongwei goats is poorly understood. MicroRNAs (miRNAs) play important roles in many biological processes, including hair growth and development. In this study, we selected skins from Zhongwei goats at different ages (45 and 108 days) that exhibited different levels of hair curvature and performed miRNA sequencing to explore the molecular mechanism of hair bending. In total, 28 significantly differentially expressed miRNAs (DE miRNAs) were identified in the three groups of samples between the two developmental stages. An analysis of the target genes of the above-mentioned DE miRNAs by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that the DE miRNAs were involved in signal pathways which were previously associated with hair bending and hair follicle development, such as the TGF-ß/SMAD, PI3K-Akt, JAK-STAT, and MAPK pathways. A comprehensive analysis of the correlations between the miRNA-seq results and issued transcriptional findings indicated that SMAD1 was a target gene of miR-26a and SMAD5 was a target gene of miR-130a. Furthermore, goat dermal papilla cells were successfully isolated and purified to determine the role of miRNAs in follicle development in vitro. The study results demonstrated that miR-130a and miR-26a had significant effects on the proliferation of dermal papilla cells. In addition, the detection results of mRNA and protein levels indicate that the overexpression of miR-26a can promote the expression of related genes in the TGF-ß/SMAD pathway, while miR-130a has the opposite substitution effect. The dual luciferase report test showed that miR-26a targeted the SMAD1 gene and reduced the expression of the SMAD1 protein in hair papillary cells. Our results identified DE microRNAs which perhaps change at the time of hair straightening in Zhongwei goats and explore the role of miR-26a and miR-130a in dermal papilla cells proliferation. The present study provided a theoretical basis to explore the mechanisms underlying the Zhongwei hair growth and curly phenotype.

Local origin or external input: modern horse origin in East Asia.

BACKGROUND: Despite decades of research, the horse domestication scenario in East Asia remains poorly understood. RESULTS: The study identified 16 haplogroups with fine-scale phylogenetic resolution using mitochondrial genomes of 317 horse samples. The time to the most recent common ancestor of the 16 haplogroups ranges from [0.8-3.1] thousand years ago (KYA) to [7.9-27.1] KYA. With combined analyses of the mitochondrial control region for 35 extant Przewalski's horses, 3544 modern and 203 ancient horses across the world, researchers provide evidence for that East Asian prevalent haplogroups Q and R were indigenously domesticated or they were involved in numerous distinct genetic components from wild horses in the southern part of East Asia. These events of haplotypes Q and R occurred during 4.7 to 16.3 KYA and 2.1 to 11.5 KYA, respectively. The diffusion of preponderant European haplogroups L from west to East Asia is consistent with the external gene input. Furthermore, genetic differences were detected between northern East Asia and southern East Asia cohorts by Principal Component Analysis, Analysis of Molecular Variance test, the χ2 test and phylogeographic analyses. CONCLUSIONS: All results suggest a complex picture of horse domestication, as well as geographic pattern in East Asia. Both local origin and external input occurred in East Asia horse populations. And besides, there are at least two different domestication or hybridization centers in East Asia.

Insights into the genetic diversity of indigenous goats and their conservation priorities.

Objective: An experiment was conducted to evaluate genetic diversity of 26 Chinese indigenous goats by 30 microsatellite markers, and then to define conservation priorities to set up the protection programs according to the weight given to within- and between-breed genetic diversity. Methods: Twenty-six representative populations of Chinese indigenous goats, 1351 total, were sampled from different geographic regions of China. Within-breed genetic diversity and marker polymorphism were estimated calculating the mean number of alleles (MNA), observed heterozygosities (HO), expected heterozygosities (He), fixation index (FIS), effective number of alleles (NE) and allelic richness (Rt). Conservation priorities were analyzed by statistical methods. Results: A relatively high level of genetic diversity was found in twenty-four populations, the exceptions were in the Daiyun and Fuqing goat populations. Within-breed kinship coefficients matrix identified seven highly inbred breeds which should be concerned. Of these, six breeds will get a negative contribution to heterozygosity when the method was based on proportional contribution to heterozygosity. Based on Weitzman or Piyasatian and Kinghorn methods, the breeds distant from others i.e. Inner Mongolia Cashmere goat, Chengdu Brown goat and Leizhou goat will obtain a high ranking. Evidence from Caballero and Toro and Fabuel et al. method prioritized Jining Gray goat, Liaoning Cashmere goat and Inner Mongolia Cashmere goat, that in agreement with results from Kinship-based methods. Conclusion: Conservation priorities had planned according to multiple methods. Our results suggest Inner Mongolia Cashmere goat (most methods), Jining Gray goat and Liaoning Cashmere goat (high contribution to heterozygosity and total diversity) should be prioritized based on most methods. Furthermore, Daiyun goat and Shannan White goat also should be prioritized based on consideration of effective population size. However, if one breed could survive in changing conditions all the time, the straightforward approach is to increase its utilization and attraction for production via mining breed germplasm characteristic.

Genome-wide discovery of lincRNAs with spatiotemporal expression patterns in the skin of goat during the cashmere growth cycle.

BACKGROUND: 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). RESULTS: 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. CONCLUSIONS: This study will facilitate future studies aimed at unravelling the function of lincRNAs in hair follicle development.

Exome sequencing reveals genetic differentiation due to high-altitude adaptation in the Tibetan cashmere goat (Capra hircus).

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Mitogenomic Meta-Analysis Identifies Two Phases of Migration in the History of Eastern Eurasian Sheep.

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637-0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8-11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20-60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5-6.8 ka (lineages A and B: ∼6.4-6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West-East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a "transportation hub" in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.

Runs of Homozygosity Detection and Selection Signature Analysis for Local Goat Breeds in Yunnan, China.

Runs of Homozygosity (ROH) are continuous homozygous DNA segments in diploid genomes, which have been used to estimate the genetic diversity, inbreeding levels, and genes associated with specific traits in livestock. In this study, we analyzed the resequencing data from 10 local goat breeds in Yunnan province of China and five additional goat populations obtained from a public database. The ROH analysis revealed 21,029 ROH segments across the 15 populations, with an average length of 1.27 Mb, a pattern of ROH, and the assessment of the inbreeding coefficient indicating genetic diversity and varying levels of inbreeding. iHS (integrated haplotype score) was used to analyze high-frequency Single-Nucleotide Polymorphisms (SNPs) in ROH regions, specific genes related to economic traits such as coat color and weight variation. These candidate genes include OCA2 (OCA2 melanosomal transmembrane protein) and MLPH (melanophilin) associated with coat color, EPHA6 (EPH receptor A6) involved in litter size, CDKAL1 (CDK5 regulatory subunit associated protein 1 like 1) and POMC (proopiomelanocortin) linked to weight variation and some putative genes associated with high-altitude adaptability and immune. This study uncovers genetic diversity and inbreeding levels within local goat breeds in Yunnan province, China. The identification of specific genes associated with economic traits and adaptability provides actionable insights for utilization and conservation efforts.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat.

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Correction: Luan et al. Identification of Critical Genes for Ovine Horn Development Based on Transcriptome during the Embryonic Period. Biology 2023, 12, 591.

In the original publication [1], there were mistakes in the order of the references, which were as follows: [...].

Identification of Critical Genes for Ovine Horn Development Based on Transcriptome during the Embryonic Period.

Horns, also known as headgear, are a unique structure of ruminants. As ruminants are globally distributed, the study of horn formation is critical not only for increasing our understanding of natural and sexual selection but also for the breeding of polled sheep breeds to facilitate modern sheep farming. Despite this, a significant number of the underlying genetic pathways in sheep horn remain unclear. In this study, to clarify the gene expression profile of horn buds and investigate the key genes in horn bud formation, RNA-sequencing (RNA-seq) technology was utilized to investigate differential gene expression in the horn buds and adjacent forehead skin of Altay sheep fetuses. There were only 68 differentially expressed genes (DEGs) identified, consisting of 58 up-regulated genes and 10 down-regulated genes. RXFP2 was differentially up-regulated in the horn buds and had the highest significance (p-value = 7.42 × 10-14). In addition, 32 DEGs were horn-related genes identified in previous studies, such as RXFP2, FOXL2, SFRP4, SFRP2, KRT1, KRT10, WNT7B, and WNT3. Further, Gene Ontology (GO) analysis showed that the DEGs were mainly enriched with regard to growth, development, and cell differentiation. Pathway analysis revealed that the Wnt signaling pathway may be responsible for horn development. Further, through combining the protein-protein interaction networks of the DEGs, it was found that the top five hub genes, namely, ACAN, SFRP2, SFRP4, WNT3, and WNT7B, were also associated with horn development. Our results suggest that only a few key genes, including RXFP2, are involved in bud formation. This study not only validates the expression of candidate genes identified at the transcriptome level in previous studies but also provides new possible marker genes for horn development, which may promote our understanding of the genetic mechanisms of horn formation.

Transcriptomic Analysis Reveals mRNA and Alternative Splicing Events in Ovine Skeletal Muscle Satellite Cells during Proliferation and Differentiation.

Skeletal muscle satellite cells (SMSCs), which are highly multifunctional muscle-derived stem cells, play an essential role in myogenesis and regeneration. Here, the transcriptional profile of SMSCs during proliferation and differentiation were constructed using the RNA-Seq method. A total of 1954 differentially expressed genes (DEGs) and 1092 differentially alternative splicing genes (DAGs) were identified including 1288 upregulated genes as well as 666 downregulated genes. GO and KEGG analyses showed that the DEGs and DAGs were enriched in the MAPK (mitogen-activated protein kinase) signaling pathway, the PI3K-Akt (phosphatidylinositol-tris-phosphate kinase 3/protein kinase B) signaling pathway, the Wnt signaling pathway, and the Ras signaling pathway. In total, 1479 alternative splice events (AS) were also identified during SMSC proliferation and differentiation. Among them, a unique AS event was the major per-mRNA splicing type, and SE was the predominant splicing pattern. Furthermore, transcription factors with AS were scanned during SMSC differentiation such as myocyte enhancer factor-2C (MEF2C) and the nuclear receptor subfamily 4 group A member 2 (NR4A2). Our results imply that MEF2C and NR4A2 can interact, and we speculate that NR4A2 and MEF2C might regulate the myogenesis of ovine SMSCs through interaction. Together, our study provides useful information on the transcriptional regulation of SMSCs during proliferation and differentiation at the transcriptional level, and provides a valuable resource for understanding the molecular mechanism of myogenesis and muscle development.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array.

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

Integration Analysis of Transcriptome and Proteome Reveal the Mechanisms of Goat Wool Bending.

Zhongwei goat is a unique Chinese native goat breed for excellent lamb fur. The pattern of flower spikes of the lamb fur was significantly reduced due to the reduction of the bending of the hair strands with growth. In order to explore the molecular mechanism underlying hair bending with growth, we performed the comprehensive analysis of transcriptome and proteome of skins from 45-days, 108-days and 365-days goat based on TMT-based quantitative proteomics and RNA-seq methods. In the three comparison groups, 356, 592 and 282 differentially expressed proteins (DEPs) were screened, respectively. KEGG pathway analysis indicated that DEPs were significantly enriched in a set of signaling pathways related to wool growth and bending, such as ECM-receptor interaction, PI3K-Akt signaling pathway, PPAR signaling pathway, protein digestion and absorption, and metabolic pathways. In addition, 20 DEPs abundance of goat skin at three development stages were examined by PRM method, which validated the reliability of proteomic data. Among them, KRT and collagen alpha family may play an important role in the development of goat hair follicle and wool bending. COL6A1, COL6A2, CRNN, TNC and LOC102178129 were identified as candidate genes based on combined analysis of transcriptome and proteome data and PRM quantification. Our results identify the differential expressed proteins as well as pathways related to the wool bending of Zhongwei goats and provide a theoretical basis for further revealing the molecular mechanism underlying wool bending of goats.

Early Weaning and Milk Substitutes Affect the Gut Microbiome, Metabolomics, and Antibody Profile in Goat Kids Suffering From Diarrhea.

Early weaning and milk substitutes increase the incidence of diarrhea in young ruminants, which may modify their gut microbiota, metabolism, immunity, and health. The aim of the study was to determine if early weaning and milk substitutes affect the gut microbiota, metabolism, and immunological status of goat kids suffering from diarrhea. The 16S rRNA gene and metagenomic sequencing in feces and serum metabolomics of early-weaned and artificially reared goat kids suffering from diarrhea (DK group) and healthy goat kids reared by their mothers (HK group) were analyzed. The serum biochemistry and immunoglobulin concentration were also determined. Several probiotics, such as Streptococcus and Lactobacillus, were higher in the feces of the DK group than in feces of the HK group. Ruminococcus sp. was elevated in the feces of HKs, likely being a biomarker for goat health. Taking all the carbohydrate-active enzyme (CAZyme) families into consideration, 20 CAZyme families were different between the groups. Compared with the DK group, the relative quantity of glycoside hydrolases (GH) and glycosyltransferase (GT) families in the HK group decreased. GT70 was only identified in HK kids participating in the activity of ß-glucuronosyltransferase during the carbohydrate metabolism. Overall, 24 metabolites were different between the groups, which were mainly involved in protein digestion and absorption, cyanoamino acid metabolism, and cholesterol metabolism. The concentrations of immunoglobulins G and M were significantly lower in the DK than in the HK group. In conclusion, our study characterized the fecal microbiota, metabolism, and immunological status of early-weaned and artificially reared goat kids suffering from diarrhea.

Brief communication: Y-chromosome haplogroup analysis indicates that Chinese Tuvans share distinctive affinity with Siberian Tuvans.

Tuvans are mainly distributed in Siberia (the Republic of Tuva), Mongolia, and China. The genetic origin of Chinese Tuvans remains controversial. The Tuvans in China were classified as Mongolians in the early 1950s by the National Ethnic Affairs Commission of China, but they defined themselves as a separate group. To resolve this dispute and determine their genetic relationships with the peoples in Central Asia, we randomly selected 150 male subjects from the Tuvans in the Altai region of Xinjiang Uygur Autonomous Region in China. Fourteen Y chromosomal markers were genotyped using the RFLP method or direct sequencing. These haplogroup data were combined with public data for 15 populations in South Siberia and Central Asia. Tuvans in both China and the Republic of Tuva had the highest frequencies of haplogroups K-M9 and Q-M242. Principal component analysis demonstrated that the Tuvans in China were of a distinct cluster, separated from their neighbors, the Mongolians and Kazakhs, which finding was consistent with the Analysis of Molecular Variances. Further population tree analysis revealed that Tuvans were on a far-separated cluster from their neighbors. Based on these results, we propose that the Tuvans (in both China and the Republic of Tuva) constitute a group distinct from Mongolians and from other Central Asia populations. However, the genetic results might be the consequence of some evolutionary forces like genetic drift and founder effect, and do not necessarily reflect their ultimate origin.