Combined analysis of microbial metagenomic and metatranscriptomic sequencing data to assess in situ physiological conditions in the premature infant gut

Microbes alter their transcriptomic profiles in response to the environment. The physiological conditions experienced by a microbial community can thus be inferred using meta-transcriptomic sequencing by comparing transcription levels of specifically chosen genes. However, this analysis requires accurate reference genomes to identify the specific genes from which RNA reads originate. In addition, such an analysis should avoid biases in transcript counts related to differences in organism abundance. In this study we describe an approach to address these difficulties. Sample-specific meta-genomic assembled genomes (MAGs) were used as reference genomes to accurately identify the origin of RNA reads, and transcript ratios of genes with opposite transcription responses were compared to eliminate biases related to differences in organismal abundance, an approach hereafter named the “diametric ratio” method. We used this approach to probe the environmental conditions experienced by Escherichia spp. in the gut of 4 premature infants, 2 of whom developed necrotizing enterocolitis (NEC), a severe inflammatory intestinal disease. We analyzed twenty fecal samples taken from four premature infants (4–6 time points from each infant), and found significantly higher diametric ratios of genes associated with low oxygen levels in samples of infants later diagnosed with NEC than in samples without NEC. We also show this method can be used for examining other physiological conditions, such as exposure to nitric oxide and osmotic pressure. These study results should be treated with caution, due to the presence of confounding factors that might also distinguish between NEC and control infants. Nevertheless, together with benchmarking analyses, we show here that the diametric ratio approach can be applied for evaluating the physiological conditions experienced by microbes in situ. Results from similar studies can be further applied for designing diagnostic methods to detect NEC in its early developmental stages.

[1]  H. Vlamakis,et al.  Microbial genes and pathways in inflammatory bowel disease , 2019, Nature Reviews Microbiology.

[2]  Brian C. Thomas,et al.  Mediterranean grassland soil C–N compound turnover is dependent on rainfall and depth, and is mediated by genomically divergent microorganisms , 2019, Nature Microbiology.

[3]  R. Gibbs,et al.  Temporal development of the gut microbiome in early childhood from the TEDDY study , 2018, Nature.

[4]  D. Hackam,et al.  Necrotizing enterocolitis: Pathophysiology from a historical context. , 2018, Seminars in pediatric surgery.

[5]  David A. Drew,et al.  Metatranscriptome of human fecal microbial communities in a cohort of adult men , 2018, Nature Microbiology.

[6]  C. Huttenhower,et al.  Dynamics of metatranscription in the inflammatory bowel disease gut microbiome , 2018, Nature Microbiology.

[7]  Brian C. Thomas,et al.  Hospitalized Premature Infants Are Colonized by Related Bacterial Strains with Distinct Proteomic Profiles , 2017, mBio.

[8]  J. Banfield,et al.  Strain-resolved analysis of hospital rooms and infants reveals overlap between the human and room microbiome , 2017, Nature Communications.

[9]  A. Bäumler,et al.  Dysbiotic Proteobacteria expansion: a microbial signature of epithelial dysfunction. , 2017, Current opinion in microbiology.

[10]  J. Banfield,et al.  dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication , 2017, The ISME Journal.

[11]  D. Rasko,et al.  Comparative genomics and transcriptomics of Escherichia coli isolates carrying virulence factors of both enteropathogenic and enterotoxigenic E. coli , 2017, Scientific Reports.

[12]  Daniel Patrick Smith,et al.  Cesarean or Vaginal Birth Does Not Impact the Longitudinal Development of the Gut Microbiome in a Cohort of Exclusively Preterm Infants , 2017, Front. Microbiol..

[13]  Michael Wagner,et al.  A peripheral epigenetic signature of immune system genes is linked to neocortical thickness and memory , 2017, Nature Communications.

[14]  Fabian Rivera-Chávez,et al.  Oxygen as a driver of gut dysbiosis. , 2017, Free radical biology & medicine.

[15]  J. Neu,et al.  Intestinal dysbiosis in preterm infants preceding necrotizing enterocolitis: a systematic review and meta-analysis , 2017, Microbiome.

[16]  A. Bhardwaj,et al.  In situ click chemistry generation of cyclooxygenase-2 inhibitors , 2017, Nature Communications.

[17]  Alexander J Probst,et al.  Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy , 2017, Nature Microbiology.

[18]  F. Turroni,et al.  Ancient bacteria of the Ötzi’s microbiome: a genomic tale from the Copper Age , 2017, Microbiome.

[19]  Z. Zeng,et al.  Development of a Novel Quantum Dots and Graphene Oxide Based FRET Assay for Rapid Detection of invA Gene of Salmonella , 2017, Front. Microbiol..

[20]  Sudipto Roy,et al.  The role of ischemia in necrotizing enterocolitis. , 2016, Journal of pediatric surgery.

[21]  Martin J. Blaser,et al.  Antibiotics, birth mode, and diet shape microbiome maturation during early life , 2016, Science Translational Medicine.

[22]  Brian C. Thomas,et al.  Evidence for persistent and shared bacterial strains against a background of largely unique gut colonization in hospitalized premature infants , 2016, The ISME Journal.

[23]  C. Bogdan.  Nitric oxide synthase in innate and adaptive immunity: an update. , 2015, Trends in immunology.

[24]  G. Braus,et al.  One Juliet and four Romeos: VeA and its methyltransferases , 2015, Front. Microbiol..

[25]  G. Weinstock,et al.  Erratum: Patterned progression of bacterial populations in the premature infant gut (Proceedings of the National Academy of Sciences of the United States of America (2014) 111 (12522-12527) DOI: 10.1073/pnas.1409497111) , 2014 .

[26]  F. Bushman,et al.  Correlation between intraluminal oxygen gradient and radial partitioning of intestinal microbiota. , 2014, Gastroenterology.

[27]  William D. Shannon,et al.  Patterned progression of bacterial populations in the premature infant gut , 2014, Proceedings of the National Academy of Sciences.

[28]  J. Clemente,et al.  The Long-Term Stability of the Human Gut Microbiota , 2013, Science.

[29]  Charlotte E. Egan,et al.  Endothelial TLR4 activation impairs intestinal microcirculatory perfusion in necrotizing enterocolitis via eNOS–NO–nitrite signaling , 2013, Proceedings of the National Academy of Sciences.

[30]  N. Ambalavanan,et al.  Feeding practices and necrotizing enterocolitis. , 2013, Clinics in perinatology.

[31]  Brian C. Thomas,et al.  Time series community genomics analysis reveals rapid shifts in bacterial species, strains, and phage during infant gut colonization , 2013, Genome research.

[32]  C. Brocker,et al.  The role of hyperosmotic stress in inflammation and disease , 2012, Biomolecular concepts.

[33]  Siu-Ming Yiu,et al.  IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth , 2012, Bioinform..

[34]  M. Teixeira,et al.  Oxidative Stress Modulates the Nitric Oxide Defense Promoted by Escherichia coli Flavorubredoxin , 2012, Journal of bacteriology.

[35]  James E Baumgardner,et al.  Oxygen-dependent regulation of nitric oxide production by inducible nitric oxide synthase. , 2011, Free radical biology & medicine.

[36]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[37]  J. Neu,et al.  Ischemia-reperfusion and neonatal intestinal injury. , 2011, The Journal of pediatrics.

[38]  Vincent J. Denef,et al.  Strain-resolved community genomic analysis of gut microbial colonization in a premature infant , 2010, Proceedings of the National Academy of Sciences.

[39]  M. Caplan,et al.  Redefining the Role of Intestinal Microbes in the Pathogenesis of Necrotizing Enterocolitis , 2010, Pediatrics.

[40]  Miriam L. Land,et al.  Trace: Tennessee Research and Creative Exchange Prodigal: Prokaryotic Gene Recognition and Translation Initiation Site Identification Recommended Citation Prodigal: Prokaryotic Gene Recognition and Translation Initiation Site Identification , 2022 .

[41]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[42]  James C Liao,et al.  Integrated network analysis identifies nitric oxide response networks and dihydroxyacid dehydratase as a crucial target in Escherichia coli , 2007, Proceedings of the National Academy of Sciences.

[43]  B. Sjöberg,et al.  NrdR Controls Differential Expression of the Escherichia coli Ribonucleotide Reductase Genes , 2007, Journal of bacteriology.

[44]  Guido Sanguinetti,et al.  Transition of Escherichia coli from Aerobic to Micro-aerobic Conditions Involves Fast and Slow Reacting Regulatory Components* , 2007, Journal of Biological Chemistry.

[45]  Nathalie Rolhion,et al.  OmpC and the σE regulatory pathway are involved in adhesion and invasion of the Crohn's disease‐associated Escherichia coli strain LF82 , 2007, Molecular microbiology.

[46]  R. Thomson,,et al.  The Roles of Bacteria and TLR4 in Rat and Murine Models of Necrotizing Enterocolitis1 , 2006, The Journal of Immunology.

[47]  M. Inouye,et al.  Transcription Regulation of ompF and ompC by a Single Transcription Factor, OmpR* , 2006, Journal of Biological Chemistry.

[48]  J. Neu.  The ‘Myth’ of Asphyxia and Hypoxia-Ischemia as Primary Causes of Necrotizing Enterocolitis , 2005, Neonatology.

[49]  E. Kolker,et al.  Transcriptome analysis of Escherichia coli using high-density oligonucleotide probe arrays. , 2002, Nucleic acids research.

[50]  G. Unden,et al.  Control of FNR function of Escherichia coli by O2 and reducing conditions. , 2002, Journal of molecular microbiology and biotechnology.

[51]  C. Nankervis,et al.  Age‐Dependent Changes in the Postnatal Intestinal Microcirculation , 2001, Microcirculation.

[52]  R. Poole,et al.  Roles of respiratory oxidases in protecting Escherichia coli K12 from oxidative stress , 2000, Antonie van Leeuwenhoek.

[53]  R. Poole,et al.  The cytochrome bd quinol oxidase in Escherichia coli has an extremely high oxygen affinity and two oxygen-binding haems: implications for regulation of activity in vivo by oxygen inhibition. , 1996, Microbiology.

[54]  R. Gunsalus,et al.  Effect of microaerophilic cell growth conditions on expression of the aerobic (cyoABCDE and cydAB) and anaerobic (narGHJI, frdABCD, and dmsABC) respiratory pathway genes in Escherichia coli , 1996, Journal of bacteriology.

[55]  R. Poole,et al.  The oxygen affinity of cytochrome bo' in Escherichia coli determined by the deoxygenation of oxyleghemoglobin and oxymyoglobin: Km values for oxygen are in the submicromolar range , 1995, Journal of bacteriology.

[56]  C. Nankervis,et al.  The role of the circulation in the pathogenesis of necrotizing enterocolitis. , 1994, Clinics in perinatology.

[57]  K. Crissinger,et al.  Regulation of hemodynamics and oxygenation in developing intestine: insight into the pathogenesis of necrotizing enterocolitis , 1994, Acta paediatrica (Oslo, Norway : 1992). Supplement.

[58]  D. Touati,et al.  Anaerobic activation of arcA transcription in Escherichia coli: roles of Fnr and ArcA , 1994, Molecular microbiology.

[59]  R. Kliegman,et al.  Pathology of neonatal necrotizing enterocolitis: a ten-year experience. , 1990, The Journal of pediatrics.

[60]  J. Guest,et al.  Regulation and over-expression of the fnr gene of Escherichia coli. , 1987, Journal of general microbiology.

[61]  J. Liao,et al.  Systems Approaches to Unraveling Nitric Oxide Response Networks in Prokaryotes , 2010 .

[62]  M. Henry,et al.  Necrotizing Enterocolitis , 1977, Pediatric Surgery.

[63]  Brian C. Thomas,et al.  Supplemental Figures Supplemental Figure 1 -esom Showing Best Blast Hit of Scaffolds , 2022 .