Widespread horse-based mobility arose around 2,200 BCE in Eurasia.

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and widespread integration as a means of transportation remains contentious2-4. Here we assemble a large collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged ~2,200 BCE (Before Common Era), through close kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than ~2,700 BCE, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly-held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe ~3,000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai ~3,500 BCE, a settlement from Central Asia associated with corrals and a subsistence economy centered on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

Unraveling the cross-talk between N6-methyladenosine modification and non-coding RNAs in breast cancer: Mechanisms and clinical implications.

In recent years, breast cancer (BC) has surpassed lung cancer as the most common malignant tumor worldwide and remains the leading cause of cancer death in women. The etiology of BC usually involves dysregulation of epigenetic mechanisms and aberrant expression of certain non-coding RNAs (ncRNAs). N6-methyladenosine (m6A), the most prevalent RNA modification in eukaryotes, widely exists in ncRNAs to affect its biosynthesis and function, and is an important regulator of tumor-related signaling pathways. Interestingly, ncRNAs can also regulate or target m6A modification, playing a key role in cancer progression. However, the m6A-ncRNAs regulatory network in BC has not been fully elucidated, especially the regulation of m6A modification by ncRNAs. Therefore, in this review, we comprehensively summarize the interaction mechanisms and biological significance of m6A modifications and ncRNAs in BC. Meanwhile, we also focused on the clinical application value of m6A modification in BC diagnosis and prognosis, intending to explore new biomarkers and potential therapeutic targets.

Porcine ZBED6 regulates growth of skeletal muscle and internal organs via multiple targets.

ZBED6 (zinc finger BED domain containing protein 6) is a transcription factor unique to placental mammals and its interaction with the IGF2 (insulin-like growth factor 2) locus plays a prominent role in the regulation of postnatal skeletal muscle growth. Here, we generated lean Bama miniature pigs by generating ZBED6-knockout (ZBED6-/-) and investigated the mechanism underlying ZBED6 in growth of muscle and internal organs of placental mammals. ZBED6-/- pigs show markedly higher lean mass, lean mass rate, larger muscle fiber area and heavier internal organs (heart and liver) than wild-type (WT) pigs. The striking phenotypic changes of ZBED6-/- pigs coincided with remarkable upregulation of IGF2 mRNA and protein expression across three tissues (gastrocnemius muscle, longissimus dorsi, heart). Despite a significant increase in liver weight, ZBED6-/- pigs show comparable levels of IGF2 expression to those of WT controls. A mechanistic study revealed that elevated methylation in the liver abrogates ZBED6 binding at the IGF2 locus, explaining the unaltered hepatic IGF2 expression in ZBED6-/- pigs. These results indicate that a ZBED6-IGF2-independent regulatory pathway exists in the liver. Transcriptome analysis and ChIP-PCR revealed new ZBED6 target genes other than IGF2, including cyclin dependent kinase inhibitor 1A (CDKN1A) and tsukushi, small leucine rich proteoglycan (TSKU), that regulates growth of muscle and liver, respectively.

Single-cell transcriptomics provide insight into metastasis-related subsets of breast cancer.

Breast cancer metastasis is a complex, multi-step process, with high cellular heterogeneity between primary and metastatic breast cancer, and more complex interactions between metastatic cancer cells and other cells in the tumor microenvironment. High-resolution single-cell transcriptome sequencing technology can visualize the heterogeneity of malignant and non-malignant cells in the tumor microenvironment in real time, especially combined with spatial transcriptome analysis, which can directly compare changes between different stages of metastatic samples. Therefore, this study takes single-cell analysis as the first perspective to deeply explore special or rare cell subpopulations related to breast cancer metastasis, systematically summarizes their functions, molecular features, and corresponding treatment strategies, which will contribute to accurately identify, understand, and target tumor metastasis-related driving events, provide a research basis for the mechanistic study of breast cancer metastasis, and provide new clues for its personalized precision treatment.

mapDATAge: a ShinyR package to chart ancient DNA data through space and time.

SUMMARY: Ancient DNA datasets are increasingly difficult to visualize for users lacking computational experience. Here, we describe mapDATAge, which aims to provide user-friendly automated modules for the interactive mapping of allele, haplogroup and/or ancestry distributions through space and time. mapDATAge enhances collaborative data sharing while assisting the assessment and reporting of spatiotemporal patterns of genetic changes. AVAILABILITY AND IMPLEMENTATION: mapDATAge is a Shiny R application designed for exploring spatiotemporal patterns in ancient DNA data through a graphical user interface. It is freely available under GNU Public License in Github: https://github.com/xuefenfei712/mapDATAge. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A study on the pesticides-loading capacity of dendritic fibrous nano silica synthesized from 1-pentanol-water microemulsion with a low oil-water ratio.

Dendritic fibrous nano silica (DFNS) represents an optimal carrier material for pesticide constituents, due to its radial accessibility channels and high specific surface area. A low-energy methodology for synthesizing DFNS at a low volume ratio of oil to water is provided by employing 1-pentanol as the oil solvent in the microemulsion synthesis system, renowned for its remarkable stability and exceptional solubility. The DFNS@KM nano-pesticide was fabricated using a diffusion supported loading (DiSupLo) method and kresoxim-methyl (KM) as the template drug. Findings from Fourier-transform infrared spectroscopy, XRD, thermogravimetric, differential thermal analysis, and Brunauer-Emmet-Teller analyzes revealed the physical adsorption of KM onto the synthesized DFNS without any chemical bonding, with KM mainly existing in an amorphous state within the channels. High-performance liquid chromatography measurements demonstrated that only the loading amount of DFNS@KM was primarily dependent on the KM to DFNS ratio, with minimal effects observed from loading temperature and time. The loading amount and encapsulation efficiency of DFNS@KM were found to be 63.09% and 84.12%, respectively. Furthermore, DFNS effectively prolonged the release of KM, with a cumulative release rate of 85.43% over 180 h. The successful loading of pesticide components into DFNS synthesized with a low oil-to-water ratio provides theoretical support for the industrialization of nano-pesticides, with significant implications for enhancing pesticide utilization, reducing pesticide dosage, augmenting agricultural efficiency, and promoting sustainable agricultural development.

Whole-genome identification of transposable elements reveals the equine repetitive element insertion polymorphism in Chinese horses.

Transposable elements (TEs) are diverse, abundant, and complicated in genomes. They not only can drive the genome evolution process but can also act as special resources for adaptation. However, little is known about the evolutionary processes that shaped horses. In this work, 126 horse assemblages involved in most horse breeds in China were used to investigate the patterns of TE variation for the first time. By using RepeatMasker and melt software, we found that the horse-specific short interspersed repetitive elements family, equine repetitive elements (ERE1), exhibited polymorphisms in horse genomes. Phylogenetic analysis based on these ERE1 loci (minor allele frequency ≥0.05) revealed three major horse groups, namely, those in northern China, southern China, and Qinghai-Tibetan, which mirrors the result determined by SNPs to some extent. The present ERE1 family emerged ~0.26 to 1.77 Mya ago, with an activity peak at ~0.49 Mya, which matches the early stage of the horse lineage and decreases after the divergence of Equus caballus and Equus ferus przewalskii. To detect the functional ERE1(s) associated with adaptation, locus-specific branch length, genome-wide association study, and absolute allele frequency difference analyses were conducted and resulted in two common protein-coding genes annotated by candidate ERE1s. They were clustered into the vascular smooth muscle contraction (p = 0.01, EDNRA) and apelin signalling pathways (p = 0.02, NRF1). Notably, ERE1 insertion into the EDNRA gene showed a higher association with adaptation among southern China horses and other horses in 15 populations and 451 individuals (p = 4.55 e-8). Our results provide a comprehensive understanding of TE variations to analyse the phylogenetic relationships and traits relevant to adaptive evolution in horses.

Tunable synthesis of dendritic fibrous nano silica using 1-pentanol-water microemulsion at low oil to water ratio.

Dendritic fibrous nanosilica (DFNS) is a suitable nano-carrier for loading pesticides with radially oriented pores and a large surface area. The microemulsion method is standard method to prepare DFNS, and 1-pentanol is taken to replace cyclohexane as an oil solvent due to its high stability and nontoxic property. The results showed that the volume ratio of 1-pentanol (oil) to water (O/W) and the molar ratio of hexadecyltrimethylammonium bromide (CTAB) to tetraethylorthosilicate (TEOS) had effected on morphology and adsorption properties of DFNS in the water-CTAB-1-pentanol-ethanol-trimethylbenzene (TMB) microemulsion system. DFNS with bicontinuous concentric lamellar morphologies can be synthesized in this microemulsion at the meager O/W volume ratio (0.025-0.045). It features a tight mesoporous structure with a thin dendritic fibrous in 0.03 to 0.04 O/W volume ratio. The particle sizes, surface areas, and porosity of DFNS were positively correlated with the addition of the silica precursor TEOS. The size of DFNS increased from 123 to about 220 nm with the CTAB/TEOS molar ratio decreasing from 0.119 to 0.050. When the molar ratio of CTAB to TEOS  = 0.119, DFNS has a smaller particle size (123 nm) with a larger surface area and abundant honeycomb mesopores; the low O/W volume ratio strategy provides theoretical support for the industrialization development of DFNS and nano-pesticides, which plays a profound role in promoting the sustainable development of pesticide reduction, efficiency and green agriculture.

A deletion variant within the FGF5 gene in goats is associated with gene expression levels and cashmere growth.

The FGF5 gene has been associated with the regulation of fibre length in mammals, including cashmere goats. A deletion variant at ~14 kb downstream of the FGF5 gene showed significant divergence between cashmere and non-cashmere goats in previous studies. In this study, we designed specific primers to genotype the deletion variant. The results of gel electrophoresis and Sanger sequencing revealed that a 507-bp deletion mutation is located at 95 454 685-95 455 191 of chromosome 6 in goats. Genotyping data from a large panel of 288 goats showed that the deletion at the FGF5 gene locus appeared to be associated with cashmere length. The deletion variant was close to fixation (frequency 0.97) in cashmere goats. Furthermore, electrophoretic mobility shift assays for evaluating DNA-protein interaction and mRNA expression levels of FGF5 suggested that the deletion variant may serve as a cis-acting element by specifically binding transcription factors to mediate quantitative changes in FGF5 mRNA expression. Our study illustrates how a structural mutation of the FGF5 gene has contributed to the cashmere growth phenotype in domestic goats. The deletion mutation within the FGF5 gene could potentially serve as a molecular marker of cashmere growth in cashmere goat breeding.

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.

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.

Identification of copy number variations in three Chinese horse breeds using 70K single nucleotide polymorphism BeadChip array.

Copy number variation (CNV), an essential form of genetic variation, has been increasingly recognized as one promising genetic marker in the analysis of animal genomes. Here, we used the Equine 70K single nucleotide polymorphism genotyping array for the genome-wide detection of CNVs in 96 horses from three diverse Chinese breeds: Debao pony (DB), Mongolian horse (MG) and Yili horse (YL). A total of 287 CNVs were determined and merged into 122 CNV regions (CNVRs) ranging from 199 bp to 2344 kb in size and distributed in a heterogeneous manner on chromosomes. These CNVRs were integrated with seven existing reports to generate a composite genome-wide dataset of 1558 equine CNVRs, revealing 69 (56.6%) novel CNVRs. The majority (69.7%) of the 122 CNVRs overlapped with 438 genes, whereas 30.3% were located in intergenic regions. Most of these genes were associated with common CNVRs, which were shared by divergent horse breeds. As many as 60, 42 and 91 genes overlapping with the breed-specific ss were identified in DB, MG and YL respectively. Among these genes, FGF11, SPEM1, PPARG, CIDEB, HIVEP1 and GALR may have potential relevance to breed-specific traits. These findings provide valuable information for understanding the equine genome and facilitating association studies of economically important traits with equine CNVRs in the future.

CERKL interacts with mitochondrial TRX2 and protects retinal cells from oxidative stress-induced apoptosis.

Mutations in the ceramide kinase-like gene (CERKL) are associated with severe retinal degeneration. However, the exact function of the encoded protein (CERKL) remains unknown. Here we show that CERKL interacts with mitochondrial thioredoxin 2 (TRX2) and maintains TRX2 in the reduced redox state. Overexpression of CERKL protects cells from apoptosis under oxidative stress, whereas suppressing CERKL renders cells more sensitive to oxidative stress. In zebrafish, CERKL protein prominently locates in the outer segment and inner segment of the photoreceptor of the retina. Knockdown of CERKL in the zebrafish leads to an increase of retinal cell death, including cone and rod photoreceptor degeneration. Signs of oxidative damage to macromolecules were also detected in CERKL deficient zebrafish retina. Our results show that CERKL interacts with TRX2 and plays a novel key role in the regulation of the TRX2 antioxidant pathway and, for the first time, provides an explanation of how mutations in CERKL may lead to retinal cell death.

The involvement of circulating miR-146a and miR-27a in patients with atherosclerotic cardiovascular disease after SARS-CoV-2 infection.

BACKGROUND: Atherosclerotic cardiovascular disease (ASCVD) is a group of clinical diseases based on pathology of atherosclerosis that is the leading cause of mortality worldwide. There is a bidirectional interaction between ASCVD and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Alterations in circulating miRNAs levels are involved in the development of ASCVD in patients infected with SARS-CoV-2, however, the correlation between ASCVD co-infection with SARS-CoV-2 and alterations of cardiac-specific miRNAs is not well understood. HYPOTHESIS: The circulating miR-146a and miR-27a are involved in bidirectional interactions between ASCVD and SARS-CoV-2 infections. METHODS: Circulating miR-146a and miR-27a levels were measured in serum and PBMCs deriving from ASCVD patients and controls after SARS-CoV-2 infection by qRT-PCR analysis. The levels of neutralizing antibodies-resistant SARS-CoV-2 in human serum was determined by competitive magnetic particle chemiluminescence method. Interleukin (IL)-6 levels were detected by automatic biochemical analyzer using electrochemiluminescence. RESULTS: Significant downregulation of circulating miR-146a and upregulation of miR-27a in ASCVD patients after infection with SARS-CoV-2 compared with controls were observed, among which the alterations were more evident in ASCVD patients comorbid with hyperlipidemia and diabetes mellitus. Consistently, correlation analysis revealed that serum miR-146a and miR-27a levels were associated with the levels of lipids and glucose, inflammatory response, and immune function in ASCVD patients. Remarkably, SARS-CoV-2 S protein RBD stimulation of PBMCs derived from both ASCVD and controls significantly downregulated miR-146a, upregulated miR-27a expression levels, and promoted IL-6 release in vitro. CONCLUSIONS: The circulating miR-146a and miR-27a are involved in metabolism, inflammation, and immune levels in patients with ASCVD after SARS-CoV-2 infection, laying the foundation for the development of strategies to prevent the risk of SARS-CoV-2 infection in ASCVD patients.

Carrier frequencies, trends, and geographical distribution of hearing loss variants in China: The pooled analysis of 2,161,984 newborns.

The aim of this study is to comprehensively investigate the prevalence and distribution patterns of three common genetic variants associated with hearing loss (HL) in Chinese neonatal population. Methods: Prior to June 30, 2023, an extensive search and screening process was conducted across multiple literature databases. R software was utilized for conducting meta-analyses, cartography, and correlation analyses. Results: Firstly, our study identified a total of 99 studies meeting the inclusion criteria. Notably, provinces such as Qinghai, Tibet, Jilin, and Heilongjiang lack large-scale genetic screening data for neonatal deafness. Secondly, in Chinese newborns, the carrier frequencies of GJB2 variants (c.235delC, c.299_300delAT) were 1.63 % (95 %CI 1.52 %-1.76 %) and 0.33 % (95 %CI 0.30 %-0.37 %); While SLC26A4 variants (c.919-2A > G, c.2168A > G) exhibited carrier rates of 0.95 % (95 %CI 0.86 %-1.04 %) and 0.17 % (95 %CI 0.15 %-0.19 %); Additionally, Mt 12S rRNA m.1555 A > G variant was found at a rate of 0.24 % (95 % CI 0.22 %-0.26 %). Thirdly, the mutation rate of GJB2 c.235delC was higher in the east of the Heihe-Tengchong line, whereas the mutation rate of Mt 12S rRNA m.1555 A > G variant exhibited the opposite pattern. Forthly, no significant correlation exhibited the opposite pattern of GJB2 variants, but there was a notable correlation among SLC26A4 variants. Lastly, strong regional distribution correlations were evident between mutation sites from different genes, particularly between SLC26A4 (c.919-2A > G and c.2168A > G) and GJB c.299_300delAT. Conclusions: The most prevalent deafness genes among Chinese neonates were GJB2 c.235delC variant, followed by SLC26A4 c.919-2A > G variant. These gene mutation rates exhibit significant regional distribution characteristics. Consequently, it is imperative to enhance genetic screening efforts to reduce the incidence of deafness in high-risk areas.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Comparison of muscle metabolomics between two Chinese horse breeds.

With their enormous muscle mass and athletic ability, horses are well-positioned as model organisms for understanding muscle metabolism. There are two different types of horse breeds-Guanzhong (GZ) horses, an athletic breed with a larger body height (~148.7 cm), and the Ningqiang pony (NQ) horses, a lower height breed generally used for ornamental purposes-both inhabited in the same region of China with obvious differences in muscle content. The main objective of this study was to evaluate the breed-specific mechanisms controlling muscle metabolism. In this study, we observed muscle glycogen, enzyme activities, and LC-MS/MS untargeted metabolomics in the gluteus medius muscle of six, each of GZ and NQ horses, to explore differentiated metabolites that are related to the development of two muscles. As expected, the glycogen content, citrate synthase, and hexokinase activity of muscle were significantly higher in GZ horses. To alleviate the false positive rate, we used both MS1 and MS2 ions for metabolite classification and differential analysis. As a result, a total of 51,535 MS1 and 541 MS2 metabolites were identified, and these metabolites can separate these two groups from each other. Notably, 40% of these metabolites were clustered into lipids and lipid-like molecules. Furthermore, 13 significant metabolites were differentially detected between GZ and NQ horses (fold change [FC] value ≥ 2, variable important in projection value ≥1, and Q value ≤ 0.05). They are primarily clustered into glutathione metabolism (GSH, p = 0.01), taurine, and hypotaurine metabolism (p < 0.05) pathways. Seven of the 13 metabolites were also found in thoroughbred racing horses, suggesting that metabolites related to antioxidants, amino acids, and lipids played a key role in the development of skeleton muscle in horses. Those metabolites related to muscle development shed a light on racing horses' routine maintenance and improvement of athletic performance.

Machine learning: a powerful tool for identifying key microbial agents associated with specific cancer types.

Machine learning (ML) includes a broad class of computer programs that improve with experience and shows unique strengths in performing tasks such as clustering, classification and regression. Over the past decade, microbial communities have been implicated in influencing the onset, progression, metastasis, and therapeutic response of multiple cancers. Host-microbe interaction may be a physiological pathway contributing to cancer development. With the accumulation of a large number of high-throughput data, ML has been successfully applied to the study of human cancer microbiomics in an attempt to reveal the complex mechanism behind cancer. In this review, we begin with a brief overview of the data sources included in cancer microbiomics studies. Then, the characteristics of the ML algorithm are briefly introduced. Secondly, the application progress of ML in cancer microbiomics is also reviewed. Finally, we highlight the challenges and future prospects facing ML in cancer microbiomics. On this basis, we conclude that the development of cancer microbiomics can not be achieved without ML, and that ML can be used to develop tumor-targeting microbial therapies, ultimately contributing to personalized and precision medicine.

Diversity of the fecal microbiota in Chinese ponies.

Introduction: The gut microbiomes of equine are plentiful and intricate, which plays an important part in the growth. However, there is a relative lack of information on the microbial diversity in the pony's gut. Methods: In this article, 118 fecal samples from DeBa pony, NiQi pony and GuZh horse were studied by 16S rRNA amplicon sequencing. Results: Diversity analysis was used to determine the difference of gut microbiota composition among different breeds. Alpha diversity analysis showed that the gut microbiota of NiQi ponies were abundant and various. Beta diversity analysis showed that the microorganisms constitution of DeBa ponies was more similar to that of NiQi ponies. LDA Effect Size (LEfSe) analysis result that the microorganism biomarkers for NiQi pony at the genus level were Phascolarctobacterium, Paludibacter, and Fibrobacter; the bacterial biomarker for DeBa pony was Streptococcus and Prevotella; and the bacterial biomarkers for GuZh horses was Treponema, Treponema Mogibacterium, Adlercreutzia, and Blautia. The correlation analysis between genera with >1% abundance and horse height found that Streptococcus (P < 0.01), Treponema (P < 0.01), Coprococcus (P < 0.01), Prevotella (P < 0.01), Phascolarctobacterium (P < 0.01), and Mogibacterium (P < 0.01) were significantly associated with horses' height. The functional prediction results indicated that DeBa pony have a microbiota functional more similar to NiQi pony. Discussion: For the first time, our results announce the species composition and structure of the gut microbiota in Chinese ponies. At the same time, our results can provide theoretical reference for further understanding the healthy breeding, feeding management and disease prevention of horses.

Fecal Microbiota Was Reshaped in UCP1 Knock-In Pigs via the Adipose-Liver-Gut Axis and Contributed to Less Fat Deposition.

The relationship between the host gut microbiota and obesity has been well documented in humans and mice; however, few studies reported the association between the gut microbiota and fat deposition in pigs. In a previous study, we generated uncoupling protein 1 (UCP1) knock-in pigs (UCP1 pigs), which exhibited a lower fat deposition phenotype. Whether the gut microbiota was reshaped in these pigs and whether the reshaped gut microbiota contributes to the lower fat content remain unknown. Here, we revealed that the fecal microbiota composition and metabolites were significantly altered under both chow diet (CD) and high-fat/high-cholesterol (HFHC) diet conditions in UCP1 pigs compared to those in wild-type (WT) pigs. The abundance of Oscillospira and Coprococcus and the level of metabolite hyodeoxycholic acid (HDCA) from feces were observed to be significantly increased in UCP1 pigs. An association analysis revealed that Oscillospira and Coprococcus were significantly negatively related to backfat thickness. In addition, after fecal microbiota transplantation (FMT), the mice that were orally gavaged with feces from UCP1 pigs exhibited less fat deposition under both CD and high-fat diet (HFD) conditions, suggesting that the fecal microbes of UCP1 pigs participate in regulating host lipid metabolism. Consistently, HDCA-treated mice also exhibited reduced fat content. Mechanistically, we found that UCP1 expression in white adipose tissue alters the gut microbiota via the adipose-liver-gut axis in pigs. Our study provides new data concerning the cross talk between host genetic variations and the gut microbiota and paves the way for the potential application of microbes or their metabolites in the regulation of fat deposition in pigs. IMPORTANCE This article investigated the effect of the ectopic expression of UCP1 on the regulation of fecal microbiota composition and metabolites and which alters the fat deposition phenotype. Bacteria, including Oscillospira and Coprococcus, and the metabolite HDCA were found to be significantly increased in feces of UCP1 pigs and had a negative relationship with backfat thickness. Mice with fecal microbiota transplantation phenocopied the UCP1 pigs under both CD and HFD conditions, suggesting that the fecal microbes of UCP1 pigs participate in regulating host lipid metabolism. Our study provides new data regarding the cross talk between host genetic variations and the gut microbiota and paves the way for the potential application of microbes or their metabolic production in the regulation of fat deposition in pigs.