Evaluation of an Adapted Semi-Automated DNA Extraction for Human Salivary Shotgun Metagenomics.

Recent attention has highlighted the importance of oral microbiota in human health and disease, e.g., in Parkinson's disease, notably using shotgun metagenomics. One key aspect for efficient shotgun metagenomic analysis relies on optimal microbial sampling and DNA extraction, generally implementing commercial solutions developed to improve sample collection and preservation, and provide high DNA quality and quantity for downstream analysis. As metagenomic studies are today performed on a large number of samples, the next evolution to increase study throughput is with DNA extraction automation. In this study, we proposed a semi-automated DNA extraction protocol for human salivary samples collected with a commercial kit, and compared the outcomes with the DNA extraction recommended by the manufacturer. While similar DNA yields were observed between the protocols, our semi-automated DNA protocol generated significantly higher DNA fragment sizes. Moreover, we showed that the oral microbiome composition was equivalent between DNA extraction methods, even at the species level. This study demonstrates that our semi-automated protocol is suitable for shotgun metagenomic analysis, while allowing for improved sample treatment logistics with reduced technical variability and without compromising the structure of the oral microbiome.

Cross talk between Paneth and tuft cells drives dysbiosis and inflammation in the gut mucosa.

Gut microbiota imbalance (dysbiosis) is increasingly associated with pathological conditions, both within and outside the gastrointestinal tract. Intestinal Paneth cells are considered to be guardians of the gut microbiota, but the events linking Paneth cell dysfunction with dysbiosis remain unclear. We report a three-step mechanism for dysbiosis initiation. Initial alterations in Paneth cells, as frequently observed in obese and inflammatorybowel diseases patients, cause a mild remodeling of microbiota, with amplification of succinate-producing species. SucnR1-dependent activation of epithelial tuft cells triggers a type 2 immune response that, in turn, aggravates the Paneth cell defaults, promoting dysbiosis and chronic inflammation. We thus reveal a function of tuft cells in promoting dysbiosis following Paneth cell deficiency and an unappreciated essential role of Paneth cells in maintaining a balanced microbiota to prevent inappropriate activation of tuft cells and deleterious dysbiosis. This succinate-tuft cell inflammation circuit may also contribute to the chronic dysbiosis observed in patients.

In Vitro Modelling of Oral Microbial Invasion in the Human Colon.

Recent advances in the human microbiome characterization have revealed significant oral microbial detection in stools of dysbiotic patients. However, little is known about the potential interactions of these invasive oral microorganisms with commensal intestinal microbiota and the host. In this proof-of-concept study, we proposed a new model of oral-to-gut invasion by the combined use of an in vitro model simulating both the physicochemical and microbial (lumen- and mucus-associated microbes) parameters of the human colon (M-ARCOL), a salivary enrichment protocol, and whole-metagenome shotgun sequencing. Oral invasion of the intestinal microbiota was simulated by injection of enriched saliva in the in vitro colon model inoculated with a fecal sample from the same healthy adult donor. The mucosal compartment of M-ARCOL was able to retain the highest species richness levels over time, while species richness levels decreased in the luminal compartment. This study also showed that oral microorganisms preferably colonized the mucosal microenvironment, suggesting potential oral-to-intestinal mucosal competitions. This new model of oral-to-gut invasion can provide useful mechanistic insights into the role of oral microbiome in various disease processes. IMPORTANCE Here, we propose a new model of oral-to-gut invasion by the combined use of an in vitro model simulating both the physicochemical and microbial (lumen- and mucus-associated microbes) parameters of the human colon (M-ARCOL), a salivary enrichment protocol, and whole-metagenome shotgun sequencing. Our study revealed the importance of integrating the mucus compartment, which retained higher microbial richness during fermentation, showed the preference of oral microbial invaders for the mucosal resources, and indicated potential oral-to-intestinal mucosal competitions. It also underlined promising opportunities to further understand mechanisms of oral invasion into the human gut microbiome, define microbe-microbe and mucus-microbe interactions in a compartmentalized fashion, and help to better characterize the potential of oral microbial invasion and their persistence in the gut.

Benchmarking second and third-generation sequencing platforms for microbial metagenomics.

Shotgun metagenomic sequencing is a common approach for studying the taxonomic diversity and metabolic potential of complex microbial communities. Current methods primarily use second generation short read sequencing, yet advances in third generation long read technologies provide opportunities to overcome some of the limitations of short read sequencing. Here, we compared seven platforms, encompassing second generation sequencers (Illumina HiSeq 300, MGI DNBSEQ-G400 and DNBSEQ-T7, ThermoFisher Ion GeneStudio S5 and Ion Proton P1) and third generation sequencers (Oxford Nanopore Technologies MinION R9 and Pacific Biosciences Sequel II). We constructed three uneven synthetic microbial communities composed of up to 87 genomic microbial strains DNAs per mock, spanning 29 bacterial and archaeal phyla, and representing the most complex and diverse synthetic communities used for sequencing technology comparisons. Our results demonstrate that third generation sequencing have advantages over second generation platforms in analyzing complex microbial communities, but require careful sequencing library preparation for optimal quantitative metagenomic analysis. Our sequencing data also provides a valuable resource for testing and benchmarking bioinformatics software for metagenomics.

Recovery of Metagenome-Assembled Genomes from a Human Fecal Sample with Pacific Biosciences High-Fidelity Sequencing.

Here, we report the recovery of 89 metagenome-assembled genomes (MAGs) derived from a human fecal sample subjected to Pacific Biosciences (PacBio) high-fidelity (HiFi) sequencing. A total of 9 MAGs consisted of complete circular contigs, and 45 MAGs were high-quality draft genomes according to the minimum information about a metagenome-assembled genome (MIMAG) standards.

StrainFLAIR: strain-level profiling of metagenomic samples using variation graphs.

Current studies are shifting from the use of single linear references to representation of multiple genomes organised in pangenome graphs or variation graphs. Meanwhile, in metagenomic samples, resolving strain-level abundances is a major step in microbiome studies, as associations between strain variants and phenotype are of great interest for diagnostic and therapeutic purposes. We developed StrainFLAIR with the aim of showing the feasibility of using variation graphs for indexing highly similar genomic sequences up to the strain level, and for characterizing a set of unknown sequenced genomes by querying this graph. On simulated data composed of mixtures of strains from the same bacterial species Escherichia coli, results show that StrainFLAIR was able to distinguish and estimate the abundances of close strains, as well as to highlight the presence of a new strain close to a referenced one and to estimate its abundance. On a real dataset composed of a mix of several bacterial species and several strains for the same species, results show that in a more complex configuration StrainFLAIR correctly estimates the abundance of each strain. Hence, results demonstrated how graph representation of multiple close genomes can be used as a reference to characterize a sample at the strain level.

Bacterial Consumption of T4 Phages.

The bacterial consumption of viruses not been reported on as of yet even though bacteria feed on almost anything. Viruses are widely distributed but have no acknowledged active biocontrol. Viral biomass undoubtedly reintegrates trophic cycles; however, the mechanisms of this phase still remain unknown. 13C-labelled T4 phages monitor the increase of the density of the bacterial DNA concomitant with the decrease of plaque forming units. We used 12C T4 phages as a control. T4 phage disappearance in wastewater sludge was found to occur mainly through predation by Aeromonadacea. Phage consumption also favours significant in situ bacterial growth. Furthermore, an isolated strain of Aeromonas was observed to grow on T4 phages as sole the source of carbon, nitrogen, and phosphorus. Bacterial species are capable of consuming bacteriophages in situ, which is likely a widespread and underestimated type of biocontrol. This assay is anticipated as a starting point for harnessing the bacterial potential in limiting the diffusion of harmful viruses within environments such as in the gut or in water.

Three transmission events of Vibrio cholerae O1 into Lusaka, Zambia.

BACKGROUND: Cholera has been present and recurring in Zambia since 1977. However, there is a paucity of data on genetic relatedness and diversity of the Vibrio cholerae isolates responsible for these outbreaks. Understanding whether the outbreaks are seeded from existing local isolates or if the outbreaks represent separate transmission events can inform public health decisions. RESULTS: Seventy-two V. cholerae isolates from outbreaks in 2009/2010, 2016, and 2017/2018 in Zambia were characterized using multilocus variable number tandem repeat analysis (MLVA) and whole genome sequencing (WGS). The isolates had eight distinct MLVA genotypes that clustered into three MLVA clonal complexes (CCs). Each CC contained isolates from only one outbreak. The results from WGS revealed both clustered and dispersed single nucleotide variants. The genetic relatedness of isolates based on WGS was consistent with the MLVA, each CC was a distinct genetic lineage and had nearest neighbors from other East African countries. In Lusaka, isolates from the same outbreak were more closely related to themselves and isolates from other countries than to isolates from other outbreaks in other years. CONCLUSIONS: Our observations are consistent with i) the presence of random mutation and alternative mechanisms of nucleotide variation, and ii) three separate transmission events of V. cholerae into Lusaka, Zambia. We suggest that locally, case-area targeted invention strategies and regionally, well-coordinated plans be in place to effectively control future cholera outbreaks.

Comparative methods for fecal sample storage to preserve gut microbial structure and function in an in vitro model of the human colon.

In vitro gut models, such as the mucosal artificial colon (M-ARCOL), provide timely and cost-efficient alternatives to in vivo assays allowing mechanistic studies to better understand the role of human microbiome in health and disease. Using such models inoculated with human fecal samples may require a critical step of stool storage. The effects of preservation methods on microbial structure and function in in vitro gut models have been poorly investigated. This study aimed to assess the impact of three commonly used preserving methods, compared with fresh fecal samples used as a control, on the kinetics of lumen and mucus-associated microbiota colonization in the M-ARCOL model. Feces from two healthy donors were frozen 48 h at - 80 °C with or without cryoprotectant (10% glycerol) or lyophilized with maltodextrin and trehalose prior to inoculation of four parallel bioreactors (e.g., fresh stool, raw stool stored at - 80 °C, stool stored at - 80 °C with glycerol and lyophilized stool). Microbiota composition and diversity (qPCR and 16S metabarcoding) as well as metabolic activity (gases and short chain fatty acids) were monitored throughout the fermentation process (9 days). All the preservative treatments allowed the maintaining inside the M-ARCOL of a complex and functional microbiota, but considering stabilization time of microbial profiles and activities (and not technical constraints associated with the supply of frozen material), our results highlighted 48 h freezing at - 80 °C without cryoprotectant as the most efficient method. These results will help scientists to determine the most accurate method for fecal storage prior to inoculation of in vitro gut microbiome models. KEY POINTS: • In vitro ARCOL model reproduces luminal and mucosal human microbiome. • Short-term storage of fecal sample influences microbial stabilization and activity. • 48 h freezing at - 80°C: most efficient method to preserve microbial ecosystem. • Scientific and technical requirements: influencers of preservation method.

Paenibacillus infection with frequent viral coinfection contributes to postinfectious hydrocephalus in Ugandan infants.

Postinfectious hydrocephalus (PIH), which often follows neonatal sepsis, is the most common cause of pediatric hydrocephalus worldwide, yet the microbial pathogens underlying this disease remain to be elucidated. Characterization of the microbial agents causing PIH would enable a shift from surgical palliation of cerebrospinal fluid (CSF) accumulation to prevention of the disease. Here, we examined blood and CSF samples collected from 100 consecutive infant cases of PIH and control cases comprising infants with non-postinfectious hydrocephalus in Uganda. Genomic sequencing of samples was undertaken to test for bacterial, fungal, and parasitic DNA; DNA and RNA sequencing was used to identify viruses; and bacterial culture recovery was used to identify potential causative organisms. We found that infection with the bacterium Paenibacillus, together with frequent cytomegalovirus (CMV) coinfection, was associated with PIH in our infant cohort. Assembly of the genome of a facultative anaerobic bacterial isolate recovered from cultures of CSF samples from PIH cases identified a strain of Paenibacillus thiaminolyticus This strain, designated Mbale, was lethal when injected into mice in contrast to the benign reference Paenibacillus strain. These findings show that an unbiased pan-microbial approach enabled characterization of Paenibacillus in CSF samples from PIH cases, and point toward a pathway of more optimal treatment and prevention for PIH and other proximate neonatal infections.

A Retrospective Case-Control Study of the Relationship between the Gut Microbiota, Enteropathy, and Child Growth.

The microbial communities residing in the child gut are thought to play an important role in child growth, although the relationship is not well understood. We examined a cohort of young children from Mirzapur, Bangladesh, prospectively over 18 months. Four fecal markers of environmental enteropathy (EE) (high levels of alpha-1-antitrypsin, calprotectin, myeloperoxidase, and neopterin) were examined and anthropometric measures obtained from a cohort of 68 children. The 16S rRNA gene of bacterial DNA was sequenced from stool samples and used to estimate amplicon sequence variants (ASVs). We age-matched children with poor growth to children with normal growth within 1 month and compared the change in abundance and diversity of ASVs over time. Elevated EE markers and poor linear growth in children were associated with changes in microbial communities in the gut. There were increased amounts of Escherichia/Shigella and Proteobacteria and decreased amounts of Prevotella associated with poorly growing children consistent with the mounting evidence supporting the relationship between intestinal inflammation, child growth, and changes in gut microbiota composition. Future research is needed to investigate this association among young children in low- and middle-income countries.

Complete Genome Sequences of the Human Pathogen Paenibacillus thiaminolyticus Mbale and Type Strain P. thiaminolyticus NRRL B-4156.

We have isolated a likely bacterial pathogen from cerebrospinal fluid from a Ugandan infant suffering from hydrocephalus. Whole-genome sequencing and assembly of the genome of the clinical isolate, as well as that of a previously deposited reference strain, identified the isolate as Paenibacillus thiaminolyticus, which has not been associated with widespread human infections.

MSPminer: abundance-based reconstitution of microbial pan-genomes from shotgun metagenomic data.

MOTIVATION: Analysis toolkits for shotgun metagenomic data achieve strain-level characterization of complex microbial communities by capturing intra-species gene content variation. Yet, these tools are hampered by the extent of reference genomes that are far from covering all microbial variability, as many species are still not sequenced or have only few strains available. Binning co-abundant genes obtained from de novo assembly is a powerful reference-free technique to discover and reconstitute gene repertoire of microbial species. While current methods accurately identify species core parts, they miss many accessory genes or split them into small gene groups that remain unassociated to core clusters. RESULTS: We introduce MSPminer, a computationally efficient software tool that reconstitutes Metagenomic Species Pan-genomes (MSPs) by binning co-abundant genes across metagenomic samples. MSPminer relies on a new robust measure of proportionality coupled with an empirical classifier to group and distinguish not only species core genes but accessory genes also. Applied to a large scale metagenomic dataset, MSPminer successfully delineates in a few hours the gene repertoires of 1661 microbial species with similar specificity and higher sensitivity than existing tools. The taxonomic annotation of MSPs reveals microorganisms hitherto unknown and brings coherence in the nomenclature of the species of the human gut microbiota. The provided MSPs can be readily used for taxonomic profiling and biomarkers discovery in human gut metagenomic samples. In addition, MSPminer can be applied on gene count tables from other ecosystems to perform similar analyses. AVAILABILITY AND IMPLEMENTATION: The binary is freely available for non-commercial users at www.enterome.com/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Molecular characterization of Vibrio cholerae responsible for cholera epidemics in Uganda by PCR, MLVA and WGS.

BACKGROUND: For almost 50 years sub-Saharan Africa, including Uganda, has experienced several outbreaks due to Vibrio cholerae. Our aim was to determine the genetic relatedness and spread of strains responsible for cholera outbreaks in Uganda. METHODOLOGY/PRINCIPAL FINDINGS: Sixty-three V. cholerae isolates collected from outbreaks in Uganda between 2014 and 2016 were tested using multiplex polymerase chain reaction (PCR), multi-locus variable number of tandem repeat analysis (MLVA) and whole genome sequencing (WGS). Three closely related MLVA clonal complexes (CC) were identified: CC1, 32% (20/63); CC2, 40% (25/63) and CC3, 28% (18/63). Each CC contained isolates from a different WGS clade. These clades were contained in the third wave of the 7th cholera pandemic strain, two clades were contained in the transmission event (T)10 lineage and other in T13. Analysing the dates and genetic relatedness revealed that V. cholerae genetic lineages spread between districts within Uganda and across national borders. CONCLUSION: The V. cholerae strains showed local and regional transmission within Uganda and the East African region. To prevent, control and eliminate cholera, these countries should implement strong cross-border collaboration and regional coordination of preventive activities.

Microbial ecology perturbation in human IgA deficiency.

Paradoxically, loss of immunoglobulin A (IgA), one of the most abundant antibodies, does not irrevocably lead to severe infections in humans but rather is associated with relatively mild respiratory infections, atopy, and autoimmunity. IgA might therefore also play covert roles, not uniquely associated with control of pathogens. We show that human IgA deficiency is not associated with massive quantitative perturbations of gut microbial ecology. Metagenomic analysis highlights an expected pathobiont expansion but a less expected depletion in some typically beneficial symbionts. Gut colonization by species usually present in the oropharynx is also reminiscent of spatial microbiota disorganization. IgM only partially rescues IgA deficiency because not all typical IgA targets are efficiently bound by IgM in the intestinal lumen. Together, IgA appears to play a nonredundant role at the forefront of the immune/microbial interface, away from the intestinal barrier, ranging from pathobiont control and regulation of systemic inflammation to preservation of commensal diversity and community networks.

Direct Detection of Shigella in Stool Specimens by Use of a Metagenomic Approach.

The underestimation of Shigella species as a cause of childhood diarrhea disease has become increasingly apparent with quantitative PCR (qPCR)-based diagnostic methods versus culture. We sought to confirm qPCR-based detection of Shigella via a metagenomics approach. Three groups of samples were selected from diarrheal cases from the Global Enteric Multicenter Study: nine Shigella culture-positive and qPCR-positive (culture+ qPCR+) samples, nine culture-negative but qPCR-positive (culture- qPCR+) samples, and nine culture-negative and qPCR-negative (culture- qPCR-) samples. Fecal DNA was sequenced using paired-end Illumina HiSeq, whereby 3.26 × 108 ± 5.6 × 107 high-quality reads were generated for each sample. We used Kraken software to compare the read counts specific to "Shigella" among the three groups. The proportions of Shigella-specific nonhuman sequence reads between culture+ qPCR+ (0.65 ± 0.42%) and culture- qPCR+ (0.55 ± 0.31%) samples were similar (Mann-Whitney U test, P = 0.627) and distinct from the culture- qPCR- group (0.17 ± 0.15%, P < 0.05). The read counts of sequences previously targeted by Shigella/enteroinvasive Escherichia coli (EIEC) qPCR assays, namely, ipaH, virA, virG, ial, ShET2, and ipaH3, were also similar between the culture+ qPCR+ and culture- qPCR+ groups and distinct from the culture- qPCR- groups (P < 0.001). Kraken performed well versus other methods: its precision and recall of Shigella were excellent at the genus level but variable at the species level. In summary, metagenomic sequencing indicates that Shigella/EIEC qPCR-positive samples are similar to those of Shigella culture-positive samples in Shigella sequence composition, thus supporting qPCR as an accurate method for detecting Shigella.

Genetic relatedness of Vibrio cholerae isolates within and between households during outbreaks in Dhaka, Bangladesh.

BACKGROUND: Household contacts of cholera patients have a 100 times higher risk of developing a cholera infection than the general population. To compare the genetic relatedness of clinical and water source Vibrio cholerae isolates from cholera patients' households across three outbreaks, we analyzed these isolates using whole-genome-sequencing (WGS) and multilocus variable-number tandem-repeat analysis (MLVA). RESULTS: The WGS analyses revealed that 80% of households had source water isolates that were more closely related to clinical isolates from the same household than to any other isolates. While in another 20% of households an isolate from a person was more closely related to clinical isolates from another household than to source water isolates from their own household. The mean pairwise differences in single nucleotide-variant (SNV) counts for isolates from the same household were significantly lower than those for different households (2.4 vs. 7.7 p < 0.0001), and isolates from the same outbreak had significantly fewer mean pairwise differences compared to isolates from different outbreaks (mean: 6.2 vs. 8.0, p < 0.0001). Based on MLVA in outbreak 1, we observed that the majority of households had clinical isolates with MLVA genotypes related to other clinical isolates and unrelated to water source isolates from the same household. While in outbreak 3, there were different MLVA genotypes between households, however within the majority of households, the clinical and water source isolates had the same MLVA genotypes. The beginning of outbreak 2 resembled outbreak 1 and the latter part resembled outbreak 3. We validated our use of MLVA by comparing it to WGS. Isolates with the identical MLVA genotype had significantly fewer mean pairwise SNV differences than those isolates with different MLVA genotypes (mean: 4.8 vs. 7.7, p < 0.0001). Furthermore, consistent with WGS results, the number of pairwise differences in the five MLVA loci for isolates within the same household was significantly lower than isolates from different households (mean: 1.6 vs. 3.0, p < 0.0001). CONCLUSION: These results suggest that transmission patterns for cholera are a combination of person-to-person and water-to-person cholera transmission with the proportions of the two modes varying within and between outbreaks.

Minimal genetic change in Vibrio cholerae in Mozambique over time: Multilocus variable number tandem repeat analysis and whole genome sequencing.

Although cholera is a major public health concern in Mozambique, its transmission patterns remain unknown. We surveyed the genetic relatedness of 75 Vibrio cholerae isolates from patients at Manhiça District Hospital between 2002-2012 and 3 isolates from river using multilocus variable-number tandem-repeat analysis (MLVA) and whole genome sequencing (WGS). MLVA revealed 22 genotypes in two clonal complexes and four unrelated genotypes. WGS revealed i) the presence of recombination, ii) 67 isolates descended monophyletically from a single source connected to Wave 3 of the Seventh Pandemic, and iii) four clinical isolates lacking the cholera toxin gene. This Wave 3 strain persisted for at least eight years in either an environmental reservoir or circulating within the human population. Our data raises important questions related to where these isolates persist and how identical isolates can be collected years apart despite our understanding of high change rate of MLVA loci and the V. cholerae molecular clock.

Comparison of inferred relatedness based on multilocus variable-number tandem-repeat analysis and whole genome sequencing of Vibrio cholerae O1.

Vibrio cholerae causes cholera, a severe diarrheal disease. Understanding the local genetic diversity and transmission of V. cholerae will improve our ability to control cholera. Vibrio cholerae isolates clustered in genetically related groups (clonal complexes, CC) by multilocus variable tandem-repeat analysis (MLVA) were compared by whole genome sequencing (WGS). Isolates in CC1 had been isolated from two geographical locations. Isolates in a second genetically distinct group, CC2, were isolated only at one location. Using WGS, CC1 isolates from both locations revealed, on average, 43.8 nucleotide differences, while those strains comprising CC2 averaged 19.7 differences. Strains from both MLVA-CCs had an average difference of 106.6. Thus, isolates comprising CC1 were more closely related (P < 10(-6)) to each other than to isolates in CC2. Within a MLVA-CC, after removing all paralogs, alternative alleles were found in all possible combinations on separate chromosomes indicative of recombination within the core genome. Including recombination did not affect the distinctiveness of the MLVA-CCs when measured by WGS. We found that WGS generally reflected the same genetic relatedness of isolates as MLVA, indicating that isolates from the same MLVA-CC shared a more recent common ancestor than isolates from the same location that clustered in a distinct MLVA-CC.

Capturing the most wanted taxa through cross-sample correlations.

The Human Microbiome Project (HMP) identified the 16S rRNA gene sequences of 'most wanted' taxa-prevalent in the healthy human microbiota but distant from previously known sequences. Since 2012, few of the corresponding genomes have been isolated and sequenced, and only through advanced isolation techniques. We demonstrate that the genomes of the most wanted taxa can be identified computationally through their correlation in abundance across multiple public metagenomic data sets. We link over 200 most wanted sequences with nearly complete genome sequences, including half of the taxa identified as high-priority targets by the HMP. The genomes we identify have strong similarity to genomes reconstructed through expensive isolation techniques, and provide a more complete functional characterization of these organisms than can be extrapolated from their 16S rRNA gene. We also provide insights into the function of organisms for which 16S rRNA gene signatures were recently reported to be associated with health and host genetic factors.