Integrative Transkingdom Analysis of the Gut Microbiome in Antibiotic Perturbation and Critical Illness

cited By 2Bacterial microbiota play a critical role in mediating local and systemic immunity, and shifts in these microbial communities have been linked to impaired outcomes in critical illness. Emerging data indicate that other intestinal organisms, including bacteriophages, viruses of eukaryotes, fungi, and protozoa, are closely interlinked with the bacterial microbiota and their host, yet their collective role during antibiotic perturbation and critical illness remains to be elucidated. We employed multi-omics factor analysis (MOFA) to systematically integrate the bacterial (16S rRNA), fungal (intergenic transcribed spacer 1 rRNA), and viral (virus discovery next generation sequencing) components of the intestinal microbiota of 33 critically ill patients with and without sepsis and 13 healthy volunteers. In addition, we quantified the absolute abundances of bacteria and fungi using 16S and 18S rRNA PCRs and characterized the short-chain fatty acids (SCFAs) butyrate, acetate, and propionate using nuclear magnetic resonance spectroscopy. We observe that a loss of the anaerobic intestinal environment is directly correlated with an overgrowth of aerobic pathobionts and their corresponding bacteriophages as well as an absolute enrichment of opportunistic yeasts capable of causing invasive disease. We also observed a strong depletion of SCFAs in both disease states, which was associated with an increased absolute abundance of fungi with respect to bacteria. Therefore, these findings illustrate the complexity of transkingdom changes following disruption of the intestinal bacterial microbiome. IMPORTANCE While numerous studies have characterized antibiotic-induced disruptions of the bacterial microbiome, few studies describe how these disruptions impact the composition of other kingdoms such as viruses, fungi, and protozoa. To address this knowledge gap, we employed MOFA to systematically integrate viral, fungal, and bacterial sequence data from critically ill patients (with and without sepsis) and healthy volunteers, both prior to and following exposure to broad-spectrum antibiotics. In doing so, we show that modulation of the bacterial component of the microbiome has implications extending beyond this kingdom alone, enabling the overgrowth of potentially invasive fungi and viruses. While numerous preclinical studies have described similar findings in vitro, we confirm these observations in humans using an integrative analytic approach. These findings underscore the potential value of multi-omics data integration tools in interrogating how different components of the microbiota contribute to disease states. In addition, our findings suggest that there is value in further studying potential adjunctive therapies using anaerobic bacteria or SCFAs to reduce fungal expansion after antibiotic exposure, which could ultimately lead to improved outcomes in the intensive care unit (ICU).Peer reviewe

Entamoeba and Giardia parasites implicated as hosts of CRESS viruses.

Metagenomic techniques have enabled genome sequencing of unknown viruses without isolation in cell culture, but information on the virus host is often lacking, preventing viral characterisation. High-throughput methods capable of identifying virus hosts based on genomic data alone would aid evaluation of their medical or biological relevance. Here, we address this by linking metagenomic discovery of three virus families in human stool samples with determination of probable hosts. Recombination between viruses provides evidence of a shared host, in which genetic exchange occurs. We utilise networks of viral recombination to delimit virus-host clusters, which are then anchored to specific hosts using (1) statistical association to a host organism in clinical samples, (2) endogenous viral elements in host genomes, and (3) evidence of host small RNA responses to these elements. This analysis suggests two CRESS virus families (Naryaviridae and Nenyaviridae) infect Entamoeba parasites, while a third (Vilyaviridae) infects Giardia duodenalis. The trio supplements five CRESS virus families already known to infect eukaryotes, extending the CRESS virus host range to protozoa. Phylogenetic analysis implies CRESS viruses infecting multicellular life have evolved independently on at least three occasions

A Sensitive Assay for Virus Discovery in Respiratory Clinical Samples

In 5–40% of respiratory infections in children, the diagnostics remain negative, suggesting that the patients might be infected with a yet unknown pathogen. Virus discovery cDNA-AFLP (VIDISCA) is a virus discovery method based on recognition of restriction enzyme cleavage sites, ligation of adaptors and subsequent amplification by PCR. However, direct discovery of unknown pathogens in nasopharyngeal swabs is difficult due to the high concentration of ribosomal RNA (rRNA) that acts as competitor. In the current study we optimized VIDISCA by adjusting the reverse transcription enzymes and decreasing rRNA amplification in the reverse transcription, using hexamer oligonucleotides that do not anneal to rRNA. Residual cDNA synthesis on rRNA templates was further reduced with oligonucleotides that anneal to rRNA but can not be extended due to 3′-dideoxy-C6-modification. With these modifications >90% reduction of rRNA amplification was established. Further improvement of the VIDISCA sensitivity was obtained by high throughput sequencing (VIDISCA-454). Eighteen nasopharyngeal swabs were analysed, all containing known respiratory viruses. We could identify the proper virus in the majority of samples tested (11/18). The median load in the VIDISCA-454 positive samples was 7.2 E5 viral genome copies/ml (ranging from 1.4 E3–7.7 E6). Our results show that optimization of VIDISCA and subsequent high-throughput-sequencing enhances sensitivity drastically and provides the opportunity to perform virus discovery directly in patient material

Two Novel Parvoviruses in Frugivorous New and Old World Bats

Bats, a globally distributed group of mammals with high ecological importance, are increasingly recognized as natural reservoir hosts for viral agents of significance to human and animal health. In the present study, we evaluated pools of blood samples obtained from two phylogenetically distant bat families, in particular from flying foxes (Pteropodidae), Eidolon helvum in West Africa, and from two species of New World leaf-nosed fruit bats (Phyllostomidae), Artibeus jamaicensis and Artibeus lituratus in Central America. A sequence-independent virus discovery technique (VIDISCA) was used in combination with high throughput sequencing to detect two novel parvoviruses: a PARV4-like virus named Eh-BtPV-1 in Eidolon helvum from Ghana and the first member of a putative new genus in Artibeus jamaicensis from Panama (Aj-BtPV-1). Those viruses were circulating in the corresponding bat colony at rates of 7–8%. Aj-BtPV-1 was also found in Artibeus lituratus (5.5%). Both viruses were detected in the blood of infected animals at high concentrations: up to 10E8 and to 10E10 copies/ml for Aj-BtPV-1 and Eh-BtPV-1 respectively. Eh-BtPV-1 was additionally detected in all organs collected from bats (brain, lungs, liver, spleen, kidneys and intestine) and spleen and kidneys were identified as the most likely sites where viral replication takes place. Our study shows that bat parvoviruses share common ancestors with known parvoviruses of humans and livestock. We also provide evidence that a variety of Parvovirinae are able to cause active infection in bats and that they are widely distributed in these animals with different geographic origin, ecologies and climatic ranges

Parasitic, bacterial, viral, immune-mediated, metabolic, and nutritional factors associated with Nodding syndrome

Nodding syndrome is a neglected, disabling and potentially fatal epileptic disorder of unknown aetiology affecting thousands of individuals mostly confined to Eastern sub-Saharan Africa. Previous studies have identified multiple associations – including O. volvulus, antileiomodin-1 antibodies, vitamin B6 deficiency, and measles virus infection – yet none is proven causal. We conducted a case-control study of children with early-stage Nodding syndrome (symptom onset <1 year). Cases and controls were identified through a household survey in the Greater Mundri area in South Sudan. A wide range of parasitic, bacterial, viral, immune-mediated, metabolic, and nutritional risk factors was investigated using conventional and state-of-the-art untargeted assays. Associations were examined by multiple logistic regression analysis and a hypothetical causal model was constructed using structural equation modelling. From 607 children with Nodding syndrome, 72 with early-stage disease were included as cases and matched to 65 household- and 44 community controls. Mansonella perstans infection (odds ratio [OR] 7.04, 95% confidence interval [CI] 2.28-21.7), Necator americanus infection (OR 2.33, 95% CI 1.02-5.3), higher antimalarial seroreactivity (OR 1.75, 95% CI 1.20-2.57), higher vitamin E concentration (OR 1.53 per standard deviation [SD] increase, 95% CI 1.07-2.19) and lower vitamin B12 concentration (OR 0.56 per SD increase, 95% CI 0.36-0.87) were associated with higher odds of NS. In a structural equation model, we hypothesized that M. perstans infection, higher vitamin E concentration and fewer viral exposures increased the risk of Nodding syndrome while lower vitamin B12 concentration, N. americanus and malaria infections resulted from having Nodding syndrome. We found no evidence that O. volvulus, antileiomodin-1 antibodies, vitamin B6 and other factors were associated with Nodding syndrome. Our results argue against several previous causal hypotheses including O. volvulus. Instead, Nodding syndrome may be caused by a complex interplay between multiple pathogens and nutrient levels. Further studies need to confirm these associations and determine the direction of effect

Complete genome characterization of two wild-type measles viruses from Vietnamese infants during the 2014 outbreak

A large measles virus outbreak occurred across Vietnam in 2014. We identified and obtained complete measles virus genomes in stool samples collected from two diarrheal pediatric patients in Dong Thap Province. These are the first complete genome sequences of circulating measles viruses in Vietnam during the 2014 measles outbreak

Genome sequences of a novel Vietnamese bat bunyavirus

To document the viral zoonotic risks in Vietnam, fecal samples were systematically collected from a number of mammals in southern Vietnam and subjected to agnostic deep sequencing. We describe here novel Vietnamese bunyavirus sequences detected in bat feces. The complete L and S segments from 14 viruses were determined

Enhanced bioinformatic profiling of VIDISCA libraries for virus detection and discovery

VIDISCA is a next-generation sequencing (NGS) library preparation method designed to enrich viral nucleic acids from samples before highly-multiplexed low depth sequencing. Reliable detection of known viruses and discovery of novel divergent viruses from NGS data require dedicated analysis tools that are both sensitive and accurate. Existing software was utilised to design a new bioinformatic workflow for high-throughput detection and discovery of viruses from VIDISCA data. The workflow leverages the VIDISCA library preparation molecular biology, specifically the use of Mse1 restriction enzyme which produces biological replicate library inserts from identical genomes. The workflow performs total metagenomic analysis for classification of non-viral sequence including parasites and host, and separately carries out virus specific analyses. Ribosomal RNA sequence is removed to increase downstream analysis speed and remaining reads are clustered at 100% identity. Known and novel viruses are sensitively detected via alignment to a virus-only protein database, and false positives are removed. A new cluster-profiling analysis takes advantage of the viral biological replicates produced by Mse1 digestion, using read clustering to flag the presence of short genomes at very high copy number. Importantly, this analysis ensures that highly repeated sequences are identified even if no homology is detected, as is shown here with the detection of a novel gokushovirus genome from human faecal matter. The workflow was validated using read data derived from serum and faeces samples taken from HIV-1 positive adults, and serum samples from pigs that were infected with atypical porcine pestivirus

Full genome virus detection in fecal samples using sensitive nucleic acid preparation, deep sequencing, and a novel iterative sequence classification algorithm

We have developed a full genome virus detection process that combines sensitive nucleic acid preparation optimised for virus identification in fecal material with Illumina MiSeq sequencing and a novel post-sequencing virus identification algorithm. Enriched viral nucleic acid was converted to double-stranded DNA and subjected to Illumina MiSeq sequencing. The resulting short reads were processed with a novel iterative Python algorithm SLIM for the identification of sequences with homology to known viruses. De novo assembly was then used to generate full viral genomes. The sensitivity of this process was demonstrated with a set of fecal samples from HIV-1 infected patients. A quantitative assessment of the mammalian, plant, and bacterial virus content of this compartment was generated and the deep sequencing data were sufficient to assembly 12 complete viral genomes from 6 virus families. The method detected high levels of enteropathic viruses that are normally controlled in healthy adults, but may be involved in the pathogenesis of HIV-1 infection and will provide a powerful tool for virus detection and for analyzing changes in the fecal virome associated with HIV-1 progression and pathogenesis