Identification of Proteome-Based Immune Subtypes of Early Hepatocellular Carcinoma and Analysis of Potential Metabolic Drivers.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related mortality worldwide, ranking fourth in frequency. The relationship between metabolic reprogramming and immune infiltration has been identified as having a crucial impact on HCC progression. However, a deeper understanding of the interplay between the immune system and metabolism in the HCC microenvironment is required. In this study, we used a proteomic dataset to identify three immune subtypes (IM1-IM3) in HCC, each of which has distinctive clinical, immune, and metabolic characteristics. Among these subtypes, IM3 was found to have the poorest prognosis, with the highest levels of immune infiltration and T-cell exhaustion. Furthermore, IM3 showed elevated glycolysis and reduced bile acid metabolism, which was strongly correlated with CD8 T cell exhaustion and regulatory T cell accumulation. Our study presents the proteomic immune stratification of HCC, revealing the possible link between immune cells and reprogramming of HCC glycolysis and bile acid metabolism, which may be a viable therapeutic strategy to improve HCC immunotherapy.

AnimalMetaOmics: a multi-omics data resources for exploring animal microbial genomes and microbiomes.

The Animal Meta-omics landscape database (AnimalMetaOmics, https://yanglab.hzau.edu.cn/animalmetaomics#/) is a comprehensive and freely available resource that includes metagenomic, metatranscriptomic, and metaproteomic data from various non-human animal species and provides abundant information on animal microbiomes, including cluster analysis of microbial cognate genes, functional gene annotations, active microbiota composition, gene expression abundance, and microbial protein identification. In this work, 55 898 microbial genomes were annotated from 581 animal species, including 42 924 bacterial genomes, 12 336 virus genomes, 496 archaea genomes and 142 fungi genomes. Moreover, 321 metatranscriptomic datasets were analyzed from 31 animal species and 326 metaproteomic datasets from four animal species, as well as the pan-genomic dynamics and compositional characteristics of 679 bacterial species and 13 archaea species from animal hosts. Researchers can efficiently access and acquire the information of cross-host microbiota through a user-friendly interface, such as species, genomes, activity levels, expressed protein sequences and functions, and pan-genome composition. These valuable resources provide an important reference for better exploring the classification, functional diversity, biological process diversity and functional genes of animal microbiota.

CRAMdb: a comprehensive database for composition and roles of microbiome in animals.

CRAMdb (a database for composition and roles of animal microbiome) is a comprehensive resource of curated and consistently annotated metagenomes for non-human animals. It focuses on the composition and roles of the microbiome in various animal species. The main goal of the CRAMdb is to facilitate the reuse of animal metagenomic data, and enable cross-host and cross-phenotype comparisons. To this end, we consistently annotated microbiomes (including 16S, 18S, ITS and metagenomics sequencing data) of 516 animals from 475 projects spanning 43 phenotype pairs to construct the database that is equipped with 9430 bacteria, 278 archaea, 2216 fungi and 458 viruses. CRAMdb provides two main contents: microbiome composition data, illustrating the landscape of the microbiota (bacteria, archaea, fungi, and viruses) in various animal species, and microbiome association data, revealing the relationships between the microbiota and various phenotypes across different animal species. More importantly, users can quickly compare the composition of the microbiota of interest cross-host or body site and the associated taxa that differ between phenotype pairs cross-host or cross-phenotype. CRAMdb is freely available at (http://www.ehbio.com/CRAMdb).

Identification of (E)-1-((1H-indol-3-yl)methylene)-4-substitute-thiosemicarbazones as potential anti-hepatic fibrosis agents.

Liver fibrosis remains a global health challenge due to its rapidly rising prevalence and limited treatment options. The orphan nuclear receptor Nur77 has been implicated in regulation of autophagy and liver fibrosis. Targeting Nur77-mediated autophagic flux may thus be a new promising strategy against hepatic fibrosis. In this study, we synthesized four types of Nur77-based thiourea derivatives to determine their anti-hepatic fibrosis activity. Among the synthesized thiourea derivatives, 9e was the most potent inhibitor of hepatic stellate cells (HSCs) proliferation and activation. This compound could directly bind to Nur77 and inhibit TGF-ß1-induced α-SMA and COLA1 expression in a Nur77-dependent manner. In vivo, 9e significantly reduced CCl4-mediated hepatic inflammation response and extracellular matrix (ECM) production, revealing that 9e is capable of blocking the progression of hepatic fibrosis. Mechanistically, 9e induced Nur77 expression and enhanced autophagic flux by inhibiting the mTORC1 signaling pathway in vitro and in vivo. Thus, the Nur77-targeted lead 9e may serve as a promising candidate for treatment of chronic liver fibrosis.

Efficacy of H2O2 inactivated bovine virus diarrhoea virus (BVDV) type 1 vaccine in mice.

BACKGROUND: Bovine viral diarrhea (BVD) is an acute febrile infectious disease caused by the bovine viral diarrhea virus (BVDV), which has brought huge economic losses to the world's cattle industry. At present, commercial inactivated BVDV vaccines may cause some adverse reactions during use. This study aims to develop a safer and more efficient inactivated BVDV vaccine. METHODS: Here, we described the generation and preclinical efficacy of a hydrogen peroxide (H2O2) inactivated BVDV type 1 vaccine in mice, and administered it separately with commercial vaccine (formaldehyde inactivated) in mice to study its efficacy. RESULTS: The BVDV type 1 IgG, IFN- γ, IL-4 and neutralizing antibody in the serum of the H2O2 inactivated vaccine group can be maintained in mice for 70 days. The IgG level reached its maximum value of 0.67 on the 42nd day, significantly higher than the commercial formaldehyde inactivated BVDV type 1 vaccine. IFN- γ and IL-4 reached their maximum values on the 28th day after immunization, at 123.16 pg/ml and 143.80 pg/ml, respectively, slightly higher than commercial vaccines, but the effect was not significant. At the same time the BVDV-1 neutralizing antibody titer reached a maximum of 12 Nu on the 42nd day post vaccination. CONCLUSIONS: The H2O2 inactivated BVDV vaccine has good safety and immunogenicity, which provides a potential solution for the further development of an efficient and safe BVDV vaccine.

Analysis of the eighth intron polymorphism of NR6A1 gene in sheep and its correlation with lumbar spine number.

For revealing molecular markers related to the traits of multiple lumbar vertebrae in sheep, we analyze the relationship between NR6A1 gene polymorphism and lumbar vertebrae number traits in Xinjiang Kazakh sheep. Lumbar muscle tissues were collected from 6-lumbar spine (L6) Kazak sheep and 7-lumbar spine (L7) Kazak sheep and the intron-8 of NR6A1 gene was amplified by PCR. The SNP locus was detected by the PCR-SSCP method. One-Way ANOVA and an Independent Chi-square Test is adopted to analyze the genotype association with lumbar spine number variation. There were two SNP loci in the intron-8 of the NR6A1 gene: IVS8-188 and IVS8-281. One-Way ANOVA and Independent Chi-square Test indicated a significant association between IVS8-281 and lumbar spine number. The SNP locus of NR6A1 gene intron 8 (IVS8-281G > A) could play a certain role in the variation of lumbar spine number in Xinjiang Kazakh sheep and demonstrates potential to accelerate the sheep breeding of selection process.

Sulforaphane activates anti-inflammatory microglia, modulating stress resilience associated with BDNF transcription.

Sulforaphane (SFN) is an organic isothiocyanate and an NF-E2-related factor-2 (Nrf2) inducer that exerts prophylactic effects on depression-like behavior in mice. However, the underlying mechanisms remain poorly understood. Brain-derived neurotrophic factor (BDNF), a neurotrophin, is widely accepted for its antidepressant effects and role in stress resilience. Here, we show that SFN confers stress resilience via BDNF upregulation and changes in abnormal dendritic spine morphology in stressed mice, which is accompanied by rectifying the irregular levels of inflammatory cytokines. Mechanistic studies demonstrated that SFN activated Nrf2 to promote BDNF transcription by binding to the exon I promoter, which is associated with increased Nrf2, and decreased methyl-CpG binding protein-2 (MeCP2), a transcriptional suppressor of BDNF, in BV2 microglial cells. Furthermore, SFN inhibited the pro-inflammatory phenotype and activated the anti-inflammatory phenotype of microglia, which was associated with increased Nrf2 and decreased MeCP2 expression in microglia of stressed mice. Hence, our findings support that Nrf2 induces BDNF transcription via upregulation of Nrf2 and downregulation of MeCP2 in microglia, which is associated with changes in the morphology of damaged dendritic spines in stressed mice. Meanwhile, the data presented here provide evidence for the application of SFN as a candidate for the prevention and intervention of depression.

Nematode Genome Announcement: Draft Genome of Meloidogyne chitwoodi, an Economically Important Pest of Potato in the Pacific Northwest.

Meloidogyne chitwoodi is one of the most devastating pests of potato in the U.S. Pacific Northwest (PNW). Nematode-infected tubers develop external as well as internal defects, making the potato tubers unmarketable, and resulting in economic losses. Draft genome assemblies of three M. chitwoodi genotypes-race 1, race 2 and race 1 pathotype Roza-were generated using Illumina and PacBio Sequel RS II sequencing. The final assemblies consist of 30, 39, and 38 polished contigs for race 1, race 2 and race 1 pathotype Roza, respectively, with average N50 of 2.37 Mb and average assembled genome size of approximately 47.41 Mb. On average, 10,508 genes were annotated for each genome. Benchmarking universal single-copy ortholog (BUSCO) analysis indicated that 69.80% of the BUSCOs were complete whereas 68.80, 0.93, and 12.67% were single copy, duplicated, and fragmented, respectively. These highly contiguous genomes will enrich resources to study potato-nematode interactions and enhance breeding efforts to develop nematode-resistant potato varieties for the PNW.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

ITRAQ-based quantitative proteomics reveals the proteome profiles of MDBK cells infected with bovine viral diarrhea virus.

BACKGROUND: Bovine viral diarrhea (BVD) which is caused by Bovine viral diarrhea virus (BVDV), is an acute, contagious disease. In spite of the use of vaccines and elimination projects, BVDV still causes severe economic losses to the cattle industry for the past few years. The current study presents a preliminary analysis of the pathogenic mechanisms from the perspective of protein expression levels in infected host cells at different points in time to elucidate the infection process associated with BVDV. METHODS: We used the isobaric tags for relative and absolute quantitation (iTRAQ) technology coupled with liquid chromatography-tandem mass spectrometric (LC-MS/MS) approach for a quantitative proteomics comparison of BVDV NADL-infected MDBK cells and non-infected cells. The functions of the proteins were deduced by functional annotation and their involvement in metabolic processes explored by KEGG pathway analysis to identify their interactions. RESULTS: There were 357 (47.6% downregulated, 52.4% upregulated infected vs. control), 101 (52.5% downregulated, 47.5% upregulated infected vs. control), and 66 (21.2% downregulated, 78.8% upregulated infected vs. control) proteins were differentially expressed (fold change > 1.5 or < 0.67) in the BVDV NADL-infected MDBK cells at 12, 24, and 48 h after infection. GO analysis showed that the differentially expressed proteins (DEPs) are mainly involved in metabolic processes, biological regulation and localization. KEGG enrichment analysis showed that some signaling pathways that involved in the regulation of BVDV NADL-infection and host resistance are significantly (P < 0.05) enriched at different stages of the BVDV NADL-infection, such as Endocytosis signaling pathway, FoxO signaling pathway, Homologous recombination signaling pathway and Lysosome pathway. CONCLUSIONS: These results revealed that the DEPs in BVDV NADL-infected MDBK cells have a wide range of regulatory effects; in addition, they provide a lot of resources for the study of host cell proteomics after BVDV infection.

MicroRNA-197-3p mediates damage to human coronary artery endothelial cells via targeting TIMP3 in Kawasaki disease.

Kawasaki disease (KD) causes cardiovascular system injury in children. However, the pathogenic mechanisms of KD have not been well defined. Recently, strong correlation between aberrant microRNAs and KD nosogenesis has been revealed. A role of microRNA-197-3p (miR-197-3p) in the pathogenesis of KD is identified in the present study. Cell proliferation assay showed human coronary artery endothelial cells (HCAECs) were suppressed by serum from KD patients, which was correlated with high levels of miR-197-3p in both KD serum and HCAECs cultured with KD serum. The inhibition of HCAECs by miR-197-3p was confirmed by cells expressing miR-197-3p mimic and miR-197-3p inhibitor. Comparative proteomics analysis and Ingenuity Pathway Analysis (IPA) revealed TIMP3 as a potential target of miR-197-3p, which was demonstrated by western blot and dual-luciferase reporter assays. Subsequently, by detecting the endothelium damage markers THBS1, VWF, and HSPG2, the role of miR-197-3p/TIMP3 in KD-induced damage to HCAECs was confirmed, which was further validated by a KD mouse model in vivo. The expressions of miR-197-3p and its target, TIMP3, are dramatically variational in KD serum and HCAECs cultured with KD serum. Increased miR-197-3p induces HCAECs abnormal by restraining TIMP3 expression directly. Hence, dysregulation of miR-197-3p/TIMP3 expression in HCAECs may be an important mechanism in cardiovascular endothelium injury in KD patients, which offers a feasible therapeutic target for KD treatment.

Screening and evaluating of long non-coding RNAs in prenatal and postnatal pituitary gland of sheep.

Long non-coding RNAs are transcribed into RNA molecules that are >200 nucleotides in length. However, the expression and function analysis of lncRNAs in the sheep pituitary gland are still lacking. In this study, we identified 1755 lncRNAs (545 annotated lncRNAs and 1210 novel lncRNAs) from RNA-seq data in the pituitary gland of embryonic and adult sheep. A total of 235 lncRNAs were differentially expressed between embryonic and adult group. We verified the presence of some lncRNAs using RT-PCR and DNA sequencing, and identified some differentially expressed lncRNAs using qPCR. We also investigated the role of cis-acting lncRNAs on target genes. GO and KEGG enrichment analysis revealed that the target genes of lncRNAs were involved in the regulation of hormones secretion and some signaling pathways in the sheep pituitary gland. Our study provides comprehensive expression profiles of lncRNAs and valuable resource for understanding their function in the pituitary gland.

Analysis of pituitary transcriptomics indicates that lncRNAs are involved in the regulation of sheep estrus.

Seasonal estrus is a key factor limiting animal fertility, and understanding the molecular mechanisms that regulate animal estrus is important for improving animal fertility. The pituitary gland, which is the most important endocrine gland in mammals, plays an important role in regulating the physiological processes such as growth, development, and reproduction of animals. Here, we used RNA-seq technology to study the expression profile of lncRNAs in the anterior pituitary of sheep during estrus and anestrus. In this study, we identified a total of 995 lncRNAs, of which 335 lncRNAs were differentially expressed in two states (including 38 up-regulated and 297 down-regulated lncRNAs). RT-qPCR verified the expression levels of several lncRNAs. Target predictive analysis revealed that these lncRNAs can act in cis or trans and regulate the expression of genes involved in the regulation of sheep estrus. Target gene enrichment analysis of differentially expressed lncRNAs indicates that these lncRNAs can regulate sheep estrus by regulating hormone metabolism and energy metabolism. Through our research, we provide the expression profile of lncRNAs in the pituitary of sheep, which provides a valuable resource for further understanding of the genetic regulation of seasonal estrus in sheep from the perspective of lncRNAs.

Identification and characterization of long non-coding RNA in prenatal and postnatal skeletal muscle of sheep.

lncRNAs are a class of transcriptional RNA molecules of >200 nucleotides in length. However, the overall expression pattern and function of lncRNAs in sheep muscle is not clear. Here, we identified 1566 lncRNAs and 404 differentially expressed lncRNAs in sheep muscle from prenatal (110 days of fetus) and postnatal (2 to 3 years old of adult sheep) developmental stages by using RNA-seq technology. Several lncRNAs were identified by using RT-PCR and DNA sequencing. The expression levels of several lncRNAs were confirmed by qRT-PCR. We analyzed the effect of lncRNAs that act cis to the target genes. lncRNA targeting genes were involved in signaling pathways associated with growth and development of muscle by GO and KEGG enrichment analysis. Through our study, we provide a comprehensive expression profile of muscle lncRNAs in sheep, which provides valuable resources for further understanding genetic regulation of muscle growth and development from the perspective of lncRNA.

Expression profiles of microRNAs in skeletal muscle of sheep by deep sequencing.

OBJECTIVE: MicroRNAs are a class of endogenous small regulatory RNAs that regulate cell proliferation, differentiation and apoptosis. Recent studies on miRNAs are mainly focused on mice, human and pig. However, the studies on miRNAs in skeletal muscle of sheep are not comprehensive. METHODS: RNA-seq technology was used to perform genomic analysis of miRNAs in prenatal and postnatal skeletal muscle of sheep. Targeted genes were predicted using miRanda software and miRNA-mRNA interactions were verified by quantitative real-time polymerase chain reaction. To further investigate the function of miRNAs, candidate targeted genes were enriched for analysis using gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment. RESULTS: The results showed total of 1,086 known miRNAs and 40 new candidate miRNAs were detected in prenatal and postnatal skeletal muscle of sheep. In addition, 345 miRNAs (151 up-regulated, 94 down-regulated) were differentially expressed. Moreover, miRanda software was performed to predict targeted genes of miRNAs, resulting in a total of 2,833 predicted targets, especially miR-381 which targeted multiple muscle-related mRNAs. Furthermore, GO and KEGG pathway analysis confirmed that targeted genes of miRNAs were involved in development of skeletal muscles. CONCLUSION: This study supplements the miRNA database of sheep, which provides valuable information for further study of the biological function of miRNAs in sheep skeletal muscle.

Expression profiles of circular RNAs in sheep skeletal muscle.

OBJECTIVE: Circular RNAs (circRNAs) are a newfound class of non-coding RNA in animals and plants. Recent studies have revealed that circRNAs play important roles in cell proliferation, differentiation, autophagy and apoptosis during development. However, there are few reports about muscle development-related circRNAs in livestock. METHODS: RNA sequencing analysis was employed to identify and annotate circRNAs from longissimus dorsi of sheep. Reverse transcription followed by real-time quantitative (q) polymerase chain reaction (PCR) analysis verified the presence of these circRNAs. Targetscan7.0 and miRanda were used to analyse the interaction of circRNA-microRNA (miRNA). To investigate the function of circRNAs, an experiment was conducted to perform enrichment analysis hosting genes of circRNAs using gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathways. RESULTS: About 75.5 million sequences were obtained from RNA libraries of sheep skeletal muscle. These sequences were mapped to 729 genes in the sheep reference genome. We identified 886 circRNAs, including numerous circular intronic RNAs and exonic circRNAs. Reverse transcription PCR (RT-PCR) and DNA sequencing analysis confirmed the presence of several circRNAs. Real-Time RT-PCR analysis exhibited resistance of sheep circRNAs to RNase R digestion. We found that many circRNAs interacted with muscle-specific miRNAs involved in growth and development of muscle, especially circ776. The GO and KEGG enrichment analysis showed that hosting genes of circRNAs was involved in muscle cell development and signaling pathway. CONCLUSION: The study provides comprehensive expression profiles of circRNAs in sheep skeletal muscle. Our study offers a large number of circRNAs to facilitate a better understanding of their roles in muscle growth. Meanwhile, we suggested that circ776 could be analyzed in future study.

Targeted Editing of Myostatin Gene in Sheep by Transcription Activator-like Effector Nucleases.

Myostatin (MSTN) is a secreted growth factor expressed in skeletal muscle and adipose tissue that negatively regulates skeletal muscle mass. Gene knockout of MSTN can result in increasing muscle mass in sheep. The objectives were to investigate whether myostatin gene can be edited in sheep by transcription activator-like effector nucleases (TALENs) in tandem with single-stranded DNA oligonucleotides (ssODNs). We designed a pair of TALENs to target a highly conserved sequence in the coding region of the sheep MSTN gene. The activity of the TALENs was verified by using luciferase single-strand annealing reporter assay in HEK 293T cell line. Co-transfection of TALENs and ssODNs oligonucleotides induced precise gene editing of myostatin gene in sheep primary fibroblasts. MSTN gene-edited cells were successfully used as nuclear donors for generating cloned embryos. TALENs combined with ssDNA oligonucleotides provide a useful approach for precise gene modification in livestock animals.

Characterization of a novel multifunctional glycoside hydrolase family in the metagenome-assembled genomes of horse gut.

The gut microbiota is a treasure trove of carbohydrate-active enzymes (CAZymes). To explore novel and efficient CAZymes, we analyzed the 4,142 metagenome-assembled genomes (MAGs) of the horse gut microbiota and found the MAG117.bin13 genome (Bacteroides fragilis) contains the highest number of polysaccharide utilisation loci sites (PULs), indicating its high capability for carbohydrate degradation. Bioinformatics analysis indicate that the PULs region of the MAG117.bin13 genome encodes many hypothetical proteins, which are important sources for exploring novel CAZymes. Interestingly, we discovered a hypothetical protein (595 amino acids). This protein exhibits potential CAZymes activity and has a lower similarity to CAZymes, we named it BfLac2275. We purified the protein using prokaryotic expression technology and studied its enzymatic function. The hydrolysis experiment of the polysaccharide substrate showed that the BfLac2275 protein has the ability to degrade α-lactose (156.94 U/mg), maltose (92.59 U/mg), raffinose (86.81 U/mg), and hyaluronic acid (5.71 U/mg). The enzyme activity is optimal at pH 5.0 and 30 ℃, indicating that the hypothetical protein BfLac2275 is a novel and multifunctional CAZymes in the glycoside hydrolases (GHs). These properties indicate that BfLac2275 has broad application prospects in many fields such as plant polysaccharide decomposition, food industry, animal feed additives and enzyme preparations. This study not only serves as a reference for exploring novel CAZymes encoded by gut microbiota but also provides an example for further studying the functional annotation of hypothetical genes in metagenomic assembly genomes.

Transcriptomic analysis provides an insight into the function of CmGH9B3, a key gene of ß-1, 4-glucanase, during the graft union healing of oriental melon scion grafted onto squash rootstock.

The melon (Cucumis melo L.) is a globally cherished and economically significant crop. The grafting technique has been widely used in the vegetative propagation of melon to promote environmental tolerance and disease resistance. However, mechanisms governing graft healing and potential incompatibilities in melons following the grafting process remain unknown. To uncover the molecular mechanism of healing of grafted melon seedlings, melon wild type (Control) and TRV-CmGH9B3 lines were obtained and grafted onto the squash rootstocks (C. moschata). Anatomical differences indicated that the healing process of the TRV-CmGH9B3 plants was slower than that of the control. A total of 335 significantly differentially expressed genes (DEGs) were detected between two transcriptomes. Most of these DEGs were down-regulated in TRV-CmGH9B3 grafted seedlings. GO and KEGG analysis showed that many metabolic, physiological, and hormonal responses were involved in graft healing, including metabolic processes, plant hormone signaling, plant MAPK pathway, and sucrose starch pathway. During the healing process of TRV-CmGH9B3 grafted seedlings, gene synthesis related to hormone signal transduction (auxin, cytokinin, gibberellin, brassinolide) was delayed. At the same time, it was found that most of the DEGs related to the sucrose pathway were down-regulated in TRV-CmGH9B3 grafted seedlings. The results showed that sugar was also involved in the healing process of melon grafted onto squash. These results deepened our understanding of the molecular mechanism of GH9B3, a key gene of ß-1, 4-glucanase. It also provided a reference for elucidating the gene mechanism and function analysis of CmGH9B3 in the process of graft union healing.

Isolation and functional analysis of acid-producing bacteria from bovine rumen.

Ruminants such as cattle rely mainly on microbes in the rumen to digest cellulose and hemicellulose from forage, and the digestion products are mainly absorbed and utilized by the host in the form of short chain fatty acids (SCFAs). This study aimed to isolate acid-producing strains from the cattle rumen and investigate their functions. A total of 980 strains of acid-producing bacteria were isolated from cattle rumen contents using a medium supplemented with bromocresol green. Combined with the test of acid production ability and 16S rRNA amplicon sequencing technology, five strains were selected based on their ability to produce relatively high levels of acid, including Bacillus pumillus, Enterococcus hirae, Enterococcus faecium, and Bacillus subtilis. Sheep were treated by gavage with a mixed bacterial suspension. The results showed that mixed bacteria significantly increased the body weight gain and feed conversion rate of sheep. To investigate the function of acid-producing bacteria in sheep, we used 16S rDNA sequencing technology to analyze the rumen microbes of sheep. We found that mixed bacteria changed the composition and abundance of sheep rumen bacteria. Among them, the abundance of Bacteroidota, Actinobacteriota, Acidobacteriota, and Proteobacteria was significantly increased, and the abundance of Firmicutes was significantly decreased, indicating that the changes in gut microbiota changed the function of the sheep rumen. The acid-producing bacteria isolated in this study can effectively promote the growth of ruminants, such as cattle and sheep, and can be used as additives to improve breeding efficiency, which lays a foundation for subsequent research on probiotics.

Expanded catalogue of metagenome-assembled genomes reveals resistome characteristics and athletic performance-associated microbes in horse.

BACKGROUND: As a domesticated species vital to humans, horses are raised worldwide as a source of mechanical energy for sports, leisure, food production, and transportation. The gut microbiota plays an important role in the health, diseases, athletic performance, and behaviour of horses. RESULTS: Here, using approximately 2.2 Tb of metagenomic sequencing data from gut samples from 242 horses, including 110 samples from the caecum and 132 samples from the rectum (faeces), we assembled 4142 microbial metagenome-assembled genomes (MAG), 4015 (96.93%) of which appear to correspond to new species. From long-read data, we successfully assembled 13 circular whole-chromosome bacterial genomes representing novel species. The MAG contained over 313,568 predicted carbohydrate-active enzymes (CAZy), over 59.77% of which had low similarity match in CAZy public databases. High abundance and diversity of antibiotic resistance genes (ARG) were identified in the MAG, likely showing the wide use of antibiotics in the management of horse. The abundances of at least 36 MAG (e.g. MAG belonging to Lachnospiraceae, Oscillospiraceae, and Ruminococcus) were higher in racehorses than in nonracehorses. These MAG enriched in racehorses contained every gene in a major pathway for producing acetate and butyrate by fibre fermentation, presenting potential for greater amount of short-chain fatty acids available to fuel athletic performance. CONCLUSION: Overall, we assembled 4142 MAG from short- and long-read sequence data in the horse gut. Our dataset represents an exhaustive microbial genome catalogue for the horse gut microbiome and provides a valuable resource for discovery of performance-enhancing microbes and studies of horse gut microbiome. Video Abstract.