SARS-CoV-2 diagnostic testing rates determine the sensitivity of genomic surveillance programs

The first step in SARS-CoV-2 genomic surveillance is testing to identify people who are infected. However, global testing rates are falling as we emerge from the acute health emergency and remain low in many low- and middle-income countries (mean = 27 tests per 100,000 people per day). We simulated COVID-19 epidemics in a prototypical low- and middle-income country to investigate how testing rates, sampling strategies and sequencing proportions jointly impact surveillance outcomes, and showed that low testing rates and spatiotemporal biases delay time to detection of new variants by weeks to months and can lead to unreliable estimates of variant prevalence, even when the proportion of samples sequenced is increased. Accordingly, investments in wider access to diagnostics to support testing rates of approximately 100 tests per 100,000 people per day could enable more timely detection of new variants and reliable estimates of variant prevalence. The performance of global SARS-CoV-2 genomic surveillance programs is fundamentally limited by access to diagnostic testing.

Higher-order interactions shape microbial interactions as microbial community complexity increases

Non-pairwise interactions, or higher-order interactions (HOIs), in microbial communities have been described as significant drivers of emergent features in microbiomes. Yet, the re-organization of microbial interactions between pairwise cultures and larger communities remains largely unexplored from a molecular perspective but is central to our understanding and further manipulation of microbial communities. Here, we used a bottom-up approach to investigate microbial interaction mechanisms from pairwise cultures up to 4-species communities from a simple microbiome (Hafnia alvei, Geotrichum candidum, Pencillium camemberti and Escherichia coli). Specifically, we characterized the interaction landscape for each species combination involving E. coli by identifying E. coli's interaction-associated mutants using an RB-TnSeq-based interaction assay. We observed a deep reorganization of the interaction-associated mutants, with very few 2-species interactions conserved all the way up to a 4-species community and the emergence of multiple HOIs. We further used a quantitative genetics strategy to decipher how 2-species interactions were quantitatively conserved in higher community compositions. Epistasis-based analysis revealed that, of the interactions that are conserved at all levels of complexity, 82% follow an additive pattern. Altogether, we demonstrate the complex architecture of microbial interactions even within a simple microbiome, and provide a mechanistic and molecular explanation of HOIs.

Biofilms preserve the transmissibility of a multi-drug resistance plasmid

Self-transmissible multidrug resistance (MDR) plasmids are a major health concern because they can spread antibiotic resistance to pathogens. Even though most pathogens form biofilms, little is known about how MDR plasmids persist and evolve in biofilms. We hypothesize that (i) biofilms act as refugia of MDR plasmids by retaining them in the absence of antibiotics longer than well-mixed planktonic populations and that (ii) the evolutionary trajectories that account for the improvement of plasmid persistence over time differ between biofilms and planktonic populations. In this study, we evolved Acinetobacter baumannii with an MDR plasmid in biofilm and planktonic populations with and without antibiotic selection. In the absence of selection, biofilm populations were better able to maintain the MDR plasmid than planktonic populations. In planktonic populations, plasmid persistence improved rapidly but was accompanied by a loss of genes required for the horizontal transfer of plasmids. In contrast, in biofilms, most plasmids retained their transfer genes, but on average, plasmid, persistence improved less over time. Our results showed that biofilms can act as refugia of MDR plasmids and favor the horizontal mode of plasmid transfer, which has important implications for the spread of MDR.

An RNA sponge controls quorum sensing dynamics and biofilm formation in Vibrio cholerae

Small regulatory RNAs (sRNAs) acting in concert with the RNA chaperone Hfq are prevalent in many bacteria and typically act by base-pairing with multiple target transcripts. In the human pathogen Vibrio cholerae, sRNAs play roles in various processes including antibiotic tolerance, competence, and quorum sensing (QS). Here, we use RIL-seq (RNA-interaction-by-ligation-and-sequencing) to identify Hfq-interacting sRNAs and their targets in V. cholerae. We find hundreds of sRNA-mRNA interactions, as well as RNA duplexes formed between two sRNA regulators. Further analysis of these duplexes identifies an RNA sponge, termed QrrX, that base-pairs with and inactivates the Qrr1-4 sRNAs, which are known to modulate the QS pathway. Transcription of qrrX is activated by QrrT, a previously uncharacterized LysR-type transcriptional regulator. Our results indicate that QrrX and QrrT are required for rapid conversion from individual to community behaviours in V. cholerae.

Genomic insights into the phylogeny and biomass-degrading enzymes of rumen ciliates

Understanding the biodiversity and genetics of gut microbiomes has important implications for host physiology and industrial enzymes, whereas most studies have been focused on bacteria and archaea, and to a lesser extent on fungi and viruses. One group, still underexplored and elusive, is ciliated protozoa, despite its importance in shaping microbiota populations. Integrating single-cell sequencing and an assembly-and-identification pipeline, we acquired 52 high-quality ciliate genomes of 22 rumen morphospecies from 11 abundant morphogenera. With these genomes, we resolved the taxonomic and phylogenetic framework that revised the 22 morphospecies into 19 species spanning 13 genera and reassigned the genus Dasytricha from Isotrichidae to a new family Dasytrichidae. Comparative genomic analyses revealed that extensive horizontal gene transfers and gene family expansion provided rumen ciliate species with a broad array of carbohydrate-active enzymes (CAZymes) to degrade all major kinds of plant and microbial carbohydrates. In particular, the genomes of Diplodiniinae and Ophryoscolecinae species encode as many CAZymes as gut fungi, and ~80% of their degradative CAZymes act on plant cell-wall. The activities of horizontally transferred cellulase and xylanase of ciliates were experimentally verified and were 2-9 folds higher than those of the inferred corresponding bacterial donors. Additionally, the new ciliate dataset greatly facilitated rumen metagenomic analyses by allowing ~12% of the metagenomic sequencing reads to be classified as ciliate sequences.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Genome sequence and characterization of Streptomyces phage Pablito, representing a new species within the genus Janusvirus

Streptomycetes are ubiquitous soil bacteria. Here we report the complete and annotated genome sequence and characterization of Streptomyces phage Pablito, isolated from a soil sample in Haarlem, the Netherlands using Streptomyces lividans as host. This phage was able to infect a diverse range of Streptomyces strains, but none belonging to the sister genus Kitasatospora. Phage Pablito has double-stranded DNA with a genome length of 49,581 base pairs encoding 76 putative proteins, of which 26 could be predicted. The presence of a serine integrase protein indicated the lysogenic nature of phage Pablito. The phage remained stable over a wide range of temperatures (25-45 °C) and at pH ≥ 7.0, but lost infectivity at temperatures above 55 °C or when the pH dropped below 6.0. This newly isolated phage is closely related to Streptomyces phage Janus and Hank144 and considered a new species classified in the genus Janusvirus, within the subfamily Arquattrovirinae.

Genome-centric analysis of short and long read metagenomes reveals uncharacterized microbiome diversity in Southeast Asians

Despite extensive efforts to address it, the vastness of uncharacterized 'dark matter' microbial genetic diversity can impact short-read sequencing based metagenomic studies. Population-specific biases in genomic reference databases can further compound this problem. Leveraging advances in hybrid assembly (using short and long reads) and Hi-C technologies in a cross-sectional survey, we deeply characterized 109 gut microbiomes from three ethnicities in Singapore to comprehensively reconstruct 4497 medium and high-quality metagenome assembled genomes, 1708 of which were missing in short-read only analysis and with >28× N50 improvement. Species-level clustering identified 70 (>10% of total) novel gut species out of 685, improved reference genomes for 363 species (53% of total), and discovered 3413 strains unique to these populations. Among the top 10 most abundant gut bacteria in our study, one of the species and >80% of strains were unrepresented in existing databases. Annotation of biosynthetic gene clusters (BGCs) uncovered more than 27,000 BGCs with a large fraction (36-88%) unrepresented in current databases, and with several unique clusters predicted to produce bacteriocins that could significantly alter microbiome community structure. These results reveal significant uncharacterized gut microbial diversity in Southeast Asian populations and highlight the utility of hybrid metagenomic references for bioprospecting and disease-focused studies.

Estimation of the mutation rate of Mycobacterium tuberculosis in cases with recurrent tuberculosis using whole genome sequencing

The study of tuberculosis latency is problematic due to the difficulty of isolating the bacteria in the dormancy state. Despite this, several in vivo approaches have been taken to mimic the latency process. Our group has studied the evolution of the bacteria in 18 cases of recurrent tuberculosis. We found that HIV positive patients develop recurrent tuberculosis earlier, generally in the first two years (p value = 0.041). The genome of the 36 Mycobacterium tuberculosis paired isolates (first and relapsed isolates) showed that none of the SNPs found within each pair was observed more than once, indicating that they were not directly related to the recurrence process. Moreover, some IS6110 movements were found in the paired isolates, indicating the presence of different clones within the patient. Finally, our results suggest that the mutation rate remains constant during all the period as no correlation was found between the number of SNPs and the time to relapse.

McaA and McaB control the dynamic positioning of a bacterial magnetic organelle

Magnetotactic bacteria are a diverse group of microorganisms that use intracellular chains of ferrimagnetic nanocrystals, produced within magnetosome organelles, to align and navigate along the geomagnetic field. Several conserved genes for magnetosome formation have been described, but the mechanisms leading to distinct species-specific magnetosome chain configurations remain unclear. Here, we show that the fragmented nature of magnetosome chains in Magnetospirillum magneticum AMB-1 is controlled by genes mcaA and mcaB. McaA recognizes the positive curvature of the inner cell membrane, while McaB localizes to magnetosomes. Along with the MamK actin-like cytoskeleton, McaA and McaB create space for addition of new magnetosomes in between pre-existing magnetosomes. Phylogenetic analyses suggest that McaA and McaB homologs are widespread among magnetotactic bacteria and may represent an ancient strategy for magnetosome positioning.

Structural and Functional Analysis of Toxin and Small RNA Gene Promoter Regions in Bacillus anthracis

It was previously demonstrated that anthrax toxin activator (AtxA) binds directly to the σ^A-like promoter region of pagA (encoding protective antigen, PA) immediately upstream of the RNA polymerase binding site. In this study, using electrophoretic mobility shift assays and in vivo analyses, we identified AtxA-binding sites in the promoter regions of the lef and cya genes (encoding lethal and edema factors, respectively) and of two Bacillus anthracis small RNAs (XrrA and XrrB). Activities of all four newly studied promoters were enhanced in the presence of CO₂/bicarbonate and AtxA, as previously seen for the pagA promoter. Notably, the cya promoter was less activated by AtxA and CO₂/bicarbonate conditions. The putative promoter of a recently described third small RNA, XrrC, showed a negligible response to AtxA and CO₂/bicarbonate. RNA polymerase binding sites of the newly studied promoters show no consensus and differ from the σ^A-like promoter region of pagA. In silico analysis of the probable AtxA binding sites in the studied promoters revealed several palindromes. All the analyzed palindromes showed very little overlap with the σ^A-like pagA promoter. It remains unclear as to how AtxA and DNA-dependent RNA-polymerase identify such diverse DNA-sequences and differentially regulate promoter activation of the studied genes. IMPORTANCE Anthrax toxin activator (AtxA) is the major virulence regulator of Bacillus anthracis, the causative agent of anthrax. Understanding AtxA's mechanism of regulation could facilitate the development of therapeutics for B. anthracis infection. We provide evidence that AtxA binds to the promoters of the cya, lef, xrrA, and xrrB genes. In vivo assays confirmed the activities of all four promoters were enhanced in the presence of AtxA and CO₂/bicarbonate, as previously seen for the pagA promoter. The cya and lef genes encode important toxin components. The xrrA and xrrB genes encode sRNAs with a suggested function as cell physiology regulators. Our data provides further evidence for the direct regulatory role of AtxA that was previously shown with the pagA promoter.

The healthy equine uterus harbors a distinct core microbiome plus a rich and diverse microbiome that varies with geographical location

The goal of this study was to understand the composition and existence of the resident uterine microbiome in healthy mares and to establish the presence of a core microbiome for the healthy equine uterus. We analyzed the microbiomes of 35 healthy mares that are long-time residents of three farms in Oklahoma, Louisiana, and Australia as well as that of 19 mares purchased from scattered owners in the Southern Mid-Western states of the United States. Over 6 million paired-end reads of the V4 region of the 16S rRNA gene were obtained resulting in 19,542 unique Amplicon Sequence Variants (ASVs). ASVs were assigned to 17 known phyla and 213 known genera. Most abundant genera across all animals were Pseudomonas (27%) followed by Lonsdalea (8%), Lactobacillus (7.5%), Escherichia/Shigella (4.5%), and Prevotella (3%). Oklahoma and Louisiana samples were dominated by Pseudomonas (75%). Lonsdalea (28%) was the most abundant genus in the Australian samples but was not found in any other region. Microbial diversity, richness, and evenness of the equine uterine microbiome is largely dependent on the geographical location of the animal. However, we observed a core uterine microbiome consisting of Lactobacillus, Escherichia/Shigella, Streptococcus, Blautia, Staphylococcus, Klebsiella, Acinetobacter, and Peptoanaerobacter.

ToxR is a c-di-GMP binding protein that modulates surface-associated behaviour in Pseudomonas aeruginosa

Pseudomonas aeruginosa uses multiple protein regulators that work in tandem to control the production of a wide range of virulence factors and facilitate rapid adaptation to diverse environmental conditions. In this opportunistic pathogen, ToxR was known to positively regulate the production of the major virulence factor exotoxin A and now, through analysis of genetic changes between two sublines of P. aeruginosa PAO1 and functional complementation of swarming, we have identified a previously unknown role of ToxR in surface-associated motility in P. aeruginosa. Further analysis revealed that ToxR had an impact on swarming motility by regulating the Rhl quorum sensing system and subsequent production of rhamnolipid surfactants. Additionally, ToxR was found to tightly bind cyclic diguanylate (c-di-GMP) and negatively affect traits controlled by this second messenger including reducing biofilm formation and the expression of Psl and Pel exopolysaccharides, necessary for attachment and sessile communities matrix scaffolding, in P. aeruginosa. Moreover, a link between the post-transcriptional regulator RsmA and toxR expression via the alternative sigma factor PvdS, induced under iron-limiting conditions, is established. This study reveals the importance of ToxR in a sophisticated regulation of free-living and biofilm-associated lifestyles, appropriate for establishing acute or chronic P. aeruginosa infections.

AhrC Negatively Regulates Streptococcus mutans Arginine Biosynthesis

Streptococcus mutans is a primary cariogenic pathogen in humans. Arginine metabolism is required for bacterial growth. In S. mutans, however, the involvement of transcription factors in regulating arginine metabolism is unclear. The purpose of this study was to investigate the function and mechanism of ArgR family transcription factors in S. mutans. Here, we identified an ArgR (arginine repressor) family transcription factor named AhrC, which negatively regulates arginine biosynthesis and biofilm formation in S. mutans. The ahrC in-frame deletion strain exhibited slow growth and significantly increased intracellular arginine content. The strain overexpressing ahrC showed reduced intracellular arginine content, decreased biofilm biomass, reduced production of water-insoluble exopolysaccharides (EPS), and different biofilm structures. Furthermore, global gene expression profiles revealed differential expression levels of 233 genes in the ahrC-deficient strain, among which genes related to arginine biosynthesis (argJ, argB, argC, argD, argF, argG, argH) were significantly upregulated. In the ahrC overexpression strain, there are 89 differentially expressed genes, mostly related to arginine biosynthesis. The conserved DNA patterns bound by AhrC were identified by electrophoretic mobility shift assay (EMSA) and DNase I footprinting. In addition, the analysis of β-galactosidase activity showed that AhrC acted as a negative regulator. Taken together, our findings suggest that AhrC is an important transcription factor that regulates arginine biosynthesis gene expression and biofilm formation in S. mutans. These findings add new aspects to the complexity of regulating the expression of genes involved in arginine biosynthesis and biofilm formation in S. mutans.

IMPORTANCE Arginine metabolism is essential for bacterial growth. The regulation of intracellular arginine metabolism in Streptococcus mutans, one of the major pathogens of dental caries, is unclear. In this study, we found that the transcription factor AhrC can directly and negatively regulate the expression of N-acetyl-gamma-glutamyl-phosphate reductase (argC), thus regulating arginine biosynthesis in S. mutans. In addition, the ahrC overexpression strain exhibited a significant decrease in biofilm and water-insoluble extracellular polysaccharides (EPS). This study adds new support to our understanding of the regulation of intracellular arginine metabolism in S. mutans.

Privacy Risks in Microbiome Research: Public Perspectives before and during a Global Pandemic

We assessed public perspectives of microbiome research privacy risks before and after a nationwide emergency was declared in the United States regarding the Covid-19 pandemic. From January to July of 2020, we conducted an online survey of perceived privacy risks of microbiome research among U.S. adults. Among 3,106 participants (the preemergency group), most expressed that the microbiome posed privacy risks similar to those associated with DNA (60.3%) or medical records (50.6%) and that they would prefer detailed explanations (70.2%) of risk in consent materials. Only 8.9% reported moderate to high familiarity with microbiome privacy risks. In adjusted analyses, individuals who participated in the study after the Covid-19 emergency was declared (the Covid-19 emergency group) were less likely to express that microbiome privacy risks were similar to those of DNA or medical records and more likely to report familiarity with the privacy risks of microbiomes. There was a trend toward increased concern after the Covid-19 emergency was declared (p = 0.053). Overall, the study revealed that many U.S. adults believe that microbiome privacy risks are similar to those associated with DNA or medical records, and they prefer detailed explanations in consent documents. Individuals who participated after the Covid-19 emergency was declared reported greater knowledge of microbiome privacy risks but had more concern.

Impact of salivary and pancreatic amylase gene copy numbers on diabetes, obesity, and functional profiles of microbiome in Northern Japanese population

Amylase genes reside in a structurally complex locus, and their copy numbers vary greatly, and several studies have reported their association with obesity. The mechanism of this effect was partially explained by changes in the oral and gut microbiome compositions; however, a detailed mechanism has been unclarified. In this study, we showed their association with diabetes in addition to obesity, and further discovered a plausible mechanism of this association based on the function of commensal bacteria. First, we confirmed that the amylase copy number in the population tends to be larger than that reported in other studies and that there is a positive association between obesity and diabetes (p = 1.89E-2 and 8.63E-3). Second, we identified that relative abundance of some genus level microbiome, Capnocytophaga, Dialister, and previously reported bacteria, were significantly associated with amylase copy numbers. Finally, through functional gene-set analysis using shotgun sequencing, we observed that the abundance of genes in the Acarbose pathway in the gut microbiome was significantly decreased with an increase in the amylase copy number (p-value = 5.80E-4). Our findings can partly explain the mechanism underlying obesity and diabetes in populations with high amylase copy numbers.

Cross-kingdom expression of synthetic genetic elements promotes discovery of metabolites in the human microbiome

Small molecules encoded by biosynthetic pathways mediate cross-species interactions and harbor untapped potential, which has provided valuable compounds for medicine and biotechnology. Since studying biosynthetic gene clusters in their native context is often difficult, alternative efforts rely on heterologous expression, which is limited by host-specific metabolic capacity and regulation. Here, we describe a computational-experimental technology to redesign genes and their regulatory regions with hybrid elements for cross-species expression in Gram-negative and -positive bacteria and eukaryotes, decoupling biosynthetic capacity from host-range constraints to activate silenced pathways. These synthetic genetic elements enabled the discovery of a class of microbiome-derived nucleotide metabolites-tyrocitabines-from Lactobacillus iners. Tyrocitabines feature a remarkable orthoester-phosphate, inhibit translational activity, and invoke unexpected biosynthetic machinery, including a class of "Amadori synthases" and "abortive" tRNA synthetases. Our approach establishes a general strategy for the redesign, expression, mobilization, and characterization of genetic elements in diverse organisms and communities.

Prevalence and characterisation of antimicrobial resistance genes and class 1 and 2 integrons in multiresistant Escherichia coli isolated from poultry production

A global increase in the populations of drug resistant bacteria exerts negative effects on animal production and human health. Our study has been focused on the assessment of resistance determinants in relation to phenotypic resistance of the 74 commensal E. coli isolates present in different ecological environments. The samples were collected from poultry litter, feces, and neck skin. Among the microorganisms isolated from the poultry litter (group A), the highest resistance was noted against AMP and DOX (100%). In the E. coli extracts from the cloacal swabs (group B), the highest resistance was observed against AMP (100%) and CIP (92%). The meat samples (group C) were characterized by resistance to AMP (100%) and STX (94.7%). Genes encoding resistance to β-lactams (bla_TEM, bla_CTX-M), fluoroquinolones (qnrA, qnrB, qnrS), aminoglycosides (strA-strB, aphA1, aac(3)-II), sulfonamides (sul1, sul2, sul3), trimethoprim (dfr1, dfr5, dfr7/17) and tetracyclines (tetA, tetB) were detected in the studied bacterial isolates. The presence of class 1 and 2 integrons was confirmed in 75% of the MDR E. coli isolates (plasmid DNA), of which 60% contained class 1 integrons, 15% contained class 2 integrons, and 11.7% carried integrons of both classes. Thus, it may be concluded that integrons are the common mediators of antimicrobial resistance among commensal multidrug resistant Escherichia coli at important stages of poultry production.

PPE51 mediates uptake of trehalose across the mycomembrane of Mycobacterium tuberculosis

The disaccharide trehalose is essential for viability of Mycobacterium tuberculosis, which synthesizes trehalose de novo but can also utilize exogenous trehalose. The mycobacterial cell wall encompasses two permeability barriers, the cytoplasmic membrane and the outer mycolic acid-containing mycomembrane. The ABC transporter LpqY-SugA-SugB-SugC has previously been demonstrated to mediate the specific uptake of trehalose across the cytoplasmic membrane. However, it is still unclear how the transport of trehalose molecules across the mycomembrane is mediated. In this study, we harnessed the antimycobacterial activity of the analogue 6-azido trehalose to select for spontaneous resistant M. tuberculosis mutants in a merodiploid strain harbouring two LpqY-SugA-SugB-SugC copies. Mutations mediating resistance to 6-azido trehalose mapped to the proline-proline-glutamate (PPE) family member PPE51 (Rv3136), which has recently been shown to be an integral mycomembrane protein involved in uptake of low-molecular weight compounds. A site-specific ppe51 gene deletion mutant of M. tuberculosis was unable to grow on trehalose as the sole carbon source. Furthermore, bioorthogonal labelling of the M. tuberculosis Δppe51 mutant incubated with 6-azido trehalose corroborated the impaired internalization. Taken together, the results indicate that the transport of trehalose and trehalose analogues across the mycomembrane of M. tuberculosis is exclusively mediated by PPE51.

Bacteroidales species in the human gut are a reservoir of antibiotic resistance genes regulated by invertible promoters

Antibiotic-resistance genes (ARGs) regulated by invertible promoters can mitigate the fitness cost of maintaining ARGs in the absence of antibiotics and could potentially prolong the persistence of ARGs in bacterial populations. However, the origin, prevalence, and distribution of these ARGs regulated by invertible promoters remains poorly understood. Here, we sought to assess the threat posed by ARGs regulated by invertible promoters by systematically searching for ARGs regulated by invertible promoters in the human gut microbiome and examining their origin, prevalence, and distribution. Through metagenomic assembly of 2227 human gut metagenomes and genomic analysis of the Unified Human Gastrointestinal Genome (UHGG) collection, we identified ARGs regulated by invertible promoters and categorized them into three classes based on the invertase-regulating phase variation. In the human gut microbiome, ARGs regulated by invertible promoters are exclusively found in Bacteroidales species. Through genomic analysis, we observed that ARGs regulated by invertible promoters have convergently originated from ARG insertions into glycan-synthesis loci that were regulated by invertible promoters at least three times. Moreover, all three classes of invertible promoters regulating ARGs are located within integrative conjugative elements (ICEs). Therefore, horizontal transfer via ICEs could explain the wide taxonomic distribution of ARGs regulated by invertible promoters. Overall, these findings reveal that glycan-synthesis loci regulated by invertible promoters in Bacteroidales species are an important hotspot for the emergence of clinically-relevant ARGs regulated by invertible promoters.

Oxidative ornithine metabolism supports non-inflammatory C. difficile colonization

The enteric pathogen Clostridioides difficile (Cd) is responsible for a toxin-mediated infection that causes more than 200,000 recorded hospitalizations and 13,000 deaths in the United States every year1. However, Cd can colonize the gut in the absence of disease symptoms. Prevalence of asymptomatic colonization by toxigenic Cd in healthy populations is high; asymptomatic carriers are at increased risk of infection compared to noncolonized individuals and may be a reservoir for transmission of Cd infection2,3. Elucidating the molecular mechanisms by which Cd persists in the absence of disease is necessary for understanding pathogenesis and developing refined therapeutic strategies. Here, we show with gut microbiome metatranscriptomic analysis that mice recalcitrant to Cd infection and inflammation exhibit increased community-wide expression of arginine and ornithine metabolic pathways. To query Cd metabolism specifically, we leverage RNA sequencing in gnotobiotic mice infected with two wild-type strains (630 and R20291) and isogenic toxin-deficient mutants of these strains to differentiate inflammation-dependent versus -independent transcriptional states. A single operon encoding oxidative ornithine degradation is consistently upregulated across non-toxigenic Cd strains. Combining untargeted and targeted metabolomics with bacterial and host genetics, we demonstrate that both diet- and host-derived sources of ornithine provide a competitive advantage to Cd, suggesting a mechanism for Cd persistence within a non-inflammatory, healthy gut.

Pneumococcal within-host diversity during colonization, transmission and treatment

Characterizing the genetic diversity of pathogens within the host promises to greatly improve surveillance and reconstruction of transmission chains. For bacteria, it also informs our understanding of inter-strain competition and how this shapes the distribution of resistant and sensitive bacteria. Here we study the genetic diversity of Streptococcus pneumoniae within 468 infants and 145 of their mothers by deep sequencing whole pneumococcal populations from 3,761 longitudinal nasopharyngeal samples. We demonstrate that deep sequencing has unsurpassed sensitivity for detecting multiple colonization, doubling the rate at which highly invasive serotype 1 bacteria were detected in carriage compared with gold-standard methods. The greater resolution identified an elevated rate of transmission from mothers to their children in the first year of the child's life. Comprehensive treatment data demonstrated that infants were at an elevated risk of both the acquisition and persistent colonization of a multidrug-resistant bacterium following antimicrobial treatment. Some alleles were enriched after antimicrobial treatment, suggesting that they aided persistence, but generally purifying selection dominated within-host evolution. Rates of co-colonization imply that in the absence of treatment, susceptible lineages outcompeted resistant lineages within the host. These results demonstrate the many benefits of deep sequencing for the genomic surveillance of bacterial pathogens.

A catalogue of 1,167 genomes from the human gut archaeome

The human gut microbiome plays an important role in health, but its archaeal diversity remains largely unexplored. In the present study, we report the analysis of 1,167 nonredundant archaeal genomes (608 high-quality genomes) recovered from human gastrointestinal tract, sampled across 24 countries and rural and urban populations. We identified previously undescribed taxa including 3 genera, 15 species and 52 strains. Based on distinct genomic features, we justify the split of the Methanobrevibacter smithii clade into two separate species, with one represented by the previously undescribed 'Candidatus Methanobrevibacter intestini'. Patterns derived from 28,581 protein clusters showed significant associations with sociodemographic characteristics such as age groups and lifestyle. We additionally show that archaea are characterized by specific genomic and functional adaptations to the host and carry a complex virome. Our work expands our current understanding of the human archaeome and provides a large genome catalogue for future analyses to decipher its impact on human physiology.

Modulation of AggR levels reveals features of virulence regulation in enteroaggregative E. coli

Enteroaggregative Escherichia coli (EAEC) strains are one of the diarrheagenic pathotypes. EAEC strains harbor a virulence plasmid (pAA2) that encodes, among other virulence determinants, the aggR gene. The expression of the AggR protein leads to the expression of several virulence determinants in both plasmids and chromosomes. In this work, we describe a novel mechanism that influences AggR expression. Because of the absence of a Rho-independent terminator in the 3'UTR, aggR transcripts extend far beyond the aggR ORF. These transcripts are prone to PNPase-mediated degradation. Structural alterations in the 3'UTR result in increased aggR transcript stability, leading to increased AggR levels. We therefore investigated the effect of increased AggR levels on EAEC virulence. Upon finding the previously described AggR-dependent virulence factors, we detected novel AggR-regulated genes that may play relevant roles in EAEC virulence. Mutants exhibiting high AggR levels because of structural alterations in the aggR 3'UTR show increased mobility and increased pAA2 conjugation frequency. Furthermore, among the genes exhibiting increased fold change values, we could identify those of metabolic pathways that promote increased degradation of arginine, fatty acids and gamma-aminobutyric acid (GABA), respectively. In this paper, we discuss how the AggR-dependent increase in specific metabolic pathways activity may contribute to EAEC virulence.

Using Cartesian Doubt To Build a Sequencing-Based View of Microbiology

The technological leap of DNA sequencing generated a tension between modern metagenomics and historical microbiology. We are forcibly harmonizing the output of a modern tool with centuries of experimental knowledge derived from culture-based microbiology. As a thought experiment, we borrow the notion of Cartesian doubt from philosopher Rene Descartes, who used doubt to build a philosophical framework from his incorrigible statement that "I think therefore I am." We aim to cast away preconceived notions and conceptualize microorganisms through the lens of metagenomic sequencing alone. Specifically, we propose funding and building analysis and engineering methods that neither search for nor rely on the assumption of independent genomes bound by lipid barriers containing discrete functional roles and taxonomies. We propose that a view of microbial communities based in sequencing will engender novel insights into metagenomic structure and may capture functional biology not reflected within the current paradigm.