Reference genomes are essential for metagenomic analyses and functional characterization of the human gut microbiota. We present the Culturable Genome Reference (CGR), a collection of 1,520 nonredundant, high-quality draft genomes generated from >6,000 bacteria cultivated from fecal samples of healthy humans. Of the 1,520 genomes, which were chosen to cover all major bacterial phyla and genera in the human gut, 264 are not represented in existing reference genome catalogs. We show that this increase in the number of reference bacterial genomes improves the rate of mapping metagenomic sequencing reads from 50% to >70%, enabling higher-resolution descriptions of the human gut microbiome. We use the CGR genomes to annotate functions of 338 bacterial species, showing the utility of this resource for functional studies. We also carry out a pan-genome analysis of 38 important human gut species, which reveals the diversity and specificity of functional enrichment between their core and dispensable genomes.

Human gut microbiome is a promising target for managing type 2 diabetes (T2D). Measures altering gut microbiota like oral intake of probiotics or berberine (BBR), a bacteriostatic agent, merit metabolic homoeostasis. We hence conducted a randomized, double-blind, placebo-controlled trial with newly diagnosed T2D patients from 20 centres in China. Four-hundred-nine eligible participants were enroled, randomly assigned (1:1:1:1) and completed a 12-week treatment of either BBR-alone, probiotics+BBR, probiotics-alone, or placebo, after a one-week run-in of gentamycin pretreatment. The changes in glycated haemoglobin, as the primary outcome, in the probiotics+BBR (least-squares mean [95% CI], −1.04[−1.19, −0.89]%) and BBR-alone group (−0.99[−1.16, −0.83]%) were significantly greater than that in the placebo and probiotics-alone groups (−0.59[−0.75, −0.44]%, −0.53[−0.68, −0.37]%, P < 0.001). BBR treatment induced more gastrointestinal side effects. Further metagenomics and metabolomic studies found that the hypoglycaemic effect of BBR is mediated by the inhibition of DCA biotransformation by Ruminococcus bromii. Therefore, our study reports a human microbial related mechanism underlying the antidiabetic effect of BBR on T2D. (Clinicaltrial.gov Identifier: NCT02861261).

The gut microbiome has been implicated in a variety of physiological states, but controversy over causality remains unresolved. Here, we performed bidirectional Mendelian randomization (MR) analyses on 3,432 Chinese individuals with whole genome, whole metagenome, anthropometric, and blood metabolic trait data. We identified 58 causal relationships between the gut microbiome and blood metabolites, and replicated 43 of them. Increased relative abundances of fecal Oscillibacter and Alistipes were causally linked to decreased triglyceride concentration. Conversely, blood metabolites such as glutamic acid appeared to decrease fecal Oxalobacter, and members of Proteobacteria were influenced by metabolites such as 5-methyltetrahydrofolic acid, alanine, glutamate, and selenium. Two-sample MR with data from Biobank Japan partly corroborated results with triglyceride and with uric acid, and also provided causal support for published fecal bacterial markers for cancer and cardiovascular diseases. This study illustrates the value of human genetic information to help prioritize gut microbial features for mechanistic and clinical studies.Metagenome-wide association studies (MWAS) using human stool samples, as well as animal models, especially germ-free mice, have pointed to a potential role of the gut microbiome in diseases such as cardiometabolic, autoimmune, neuropsychiatric disorders and cancer, with mechanistic investigations for diseases such as obesity, colorectal cancer and schizophrenia [1][2][3][4] . Twin-based heritability estimation and more recent metagenome-genome-wide association studies (M-GWAS) have questioned the traditional view of the gut microbiota as a purely environmental factor 5-9 , although the extent of the genetic influence remains controversial 7,10 . Yet, all these published cohorts, except for human sequences in the metagenomic data of HMP (Human Microbiome Project), utilized array data for human genetics, and most of them had 16S rRNA gene amplicon sequencing for the fecal microbiota [5][6][7][8][9] .As the gut microbiome is considered to be highly dynamic, causality has been an unresolved issue in the field. Mendelian randomization (MR) 11 offers an opportunity to distinguish between causal and non-causal effects from cross-sectional data, without animal studies or randomized controlled trials. An early study used MR to look at the gut microbiota and ischemic heart disease 12 . Recently, a study used MR to confirm that increased relative abundance of bacteria producing the fecal volatile short-chain fatty acid (SCFA) butyrate was causally linked to improved insulin response to oral glucose challenge; in contrast, another fecal SCFA, propionate, was causally related to an increased risk of T2D 13 . However, both studies used genotype data, and it was not clear to what extent the genetic factors explained the microbial feature of interest.In this study, we present a large-scale M-GWAS using whole genome and fecal microbiome, followed by bidirectional MR for the fecal microbiome and anthropometric...

Dickeya zeae is the causal agent of rice foot rot and maize stalk rot diseases, which could cause severe economic losses. The pathogen is known to produce two phytotoxins known as zeamine and zeamine II which are also potent antibiotics against both gram-positive and gram-negative bacteria pathogens. Zeamine II is a long-chain aminated polyketide and zeamine shares the same polyketide structure as zeamine II, with an extra valine derivative moiety conjugated to the primary amino group of zeamine II. In this study, we have identified a gene designated as zmsK encoding a putative nonribosomal peptide synthase (NRPS) by screening of the transposon mutants defective in zeamine production. Different from most known NRPS enzymes, which are commonly multidomain proteins, ZmsK contains only a condensation domain. High-performance liquid chromatography and mass spectrometry analyses showed that the ZmsK deletion mutant produced only zeamine II but not zeamine, suggesting that ZmsK catalyzes the amide bond formation by using zeamine II as a substrate to generate zeamine. We also present evidence that a partially conserved catalytic motif within the condensation domain is critical for zeamine production. Furthermore, we show that deletion of zmsK substantially decreased the total antimicrobial activity and virulence of D. zeae. Our findings provide a new insight into the biosynthesis pathway of zeamines and the virulence mechanisms of the bacterial pathogen D. zeae.

Butyrate-producing bacteria can biosynthesize butyrate and alleviate inflammatory diseases. However, few studies have reported that the genus Collinsella has the ability to produce butyric acid. Here, our study depicts a Collinsella strain, which is a rod-shaped obligate anaerobe that is able to produce butyric acid. This microorganism was isolated from a human gut, and the optimal growth conditions were found to be 37 °C on PYG medium with pH 6.5. The 16S rRNA gene sequence demonstrated that this microorganism shared 99.93% similarity with C. aerofaciens ATCC 25986T, which was higher than the threshold (98.65%) for differentiating two species. Digital DNA–DNA hybridization and average nucleotide identity values also supported that this microorganism belonged to the species C. aerofaciens. Distinct phenotypic characteristics between TF06-26 and the type strain of C. aerofaciens, such as the fermentation of D-lactose, D-fructose and D-maltose, positive growth under pH 5 and 0.2% (w/v) cholate, suggested this strain was a novel subspecies. Comparative genome analysis revealed that butyric acid kinase and phosphate butyryltransferase enzymes were coded exclusively by this strain, indicating a specific butyric acid-producing function of this C. aerofaciens subspecies within the genus Collinsella. Thus, Collinsella aerofaciens subsp. shenzhenensis subsp. nov. was proposed, with set strain TF06-26T (=CGMCC 1.5216T = DSM 105138T) as the type strain.

The oral microbiota contains billions of microbial cells, which could contribute to diseases in many body sites. Challenged by eating, drinking, and dental hygiene on a daily basis, the oral microbiota is regarded as highly dynamic. Here, we report significant human genomic associations with the oral metagenome from more than 1915 individuals, for both the tongue dorsum (n = 2017) and saliva (n = 1915). We identified five genetic loci associated with oral microbiota at study-wide significance (p < 3.16 × 10−11). Four of the five associations were well replicated in an independent cohort of 1439 individuals: rs1196764 at APPL2 with Prevotella jejuni, Oribacterium uSGB 3339 and Solobacterium uSGB 315; rs3775944 at the serum uric acid transporter SLC2A9 with Oribacterium uSGB 1215, Oribacterium uSGB 489 and Lachnoanaerobaculum umeaense; rs4911713 near OR11H1 with species F0422 uSGB 392; and rs36186689 at LOC105371703 with Eggerthia. Further analyses confirmed 84% (386/455 for tongue dorsum) and 85% (391/466 for saliva) of host genome-microbiome associations including six genome-wide significant associations mutually validated between the two niches. As many of the oral microbiome-associated genetic variants lie near miRNA genes, we tentatively validated the potential of host miRNAs to modulate the growth of specific oral bacteria. Human genetics accounted for at least 10% of oral microbiome compositions between individuals. Machine learning models showed that polygenetic risk scores dominated over oral microbiome in predicting risk of dental diseases such as dental calculus and gingival bleeding. These findings indicate that human genetic differences are one explanation for a stable or recurrent oral microbiome in each individual.

The life-threatening coronaviruses MERS-CoV, SARS-CoV-1 and SARS-CoV-2 (SARS-CoV-1/2) have caused and will continue to cause enormous morbidity and mortality to humans. Virus-encoded noncoding RNAs are poorly understood in coronaviruses. Data mining of viral-infection-related RNA-sequencing data has resulted in the identification of 28 754, 720 and 3437 circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2, respectively. MERS-CoV exhibits much more prominent ability to encode circRNAs in all genomic regions than those of SARS-CoV-1/2. Viral circRNAs typically exhibit low expression levels. Moreover, majority of the viral circRNAs exhibit expressions only in the late stage of viral infection. Analysis of the competitive interactions of viral circRNAs, human miRNAs and mRNAs in MERS-CoV infections reveals that viral circRNAs up-regulated genes related to mRNA splicing and processing in the early stage of viral infection, and regulated genes involved in diverse functions including cancer, metabolism, autophagy, viral infection in the late stage of viral infection. Similar analysis in SARS-CoV-2 infections reveals that its viral circRNAs down-regulated genes associated with metabolic processes of cholesterol, alcohol, fatty acid and up-regulated genes associated with cellular responses to oxidative stress in the late stage of viral infection. A few genes regulated by viral circRNAs from both MERS-CoV and SARS-CoV-2 were enriched in several biological processes such as response to reactive oxygen and centrosome localization. This study provides the first glimpse into viral circRNAs in three deadly coronaviruses and would serve as a valuable resource for further studies of circRNAs in coronaviruses.

No abstract