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

Lactobacillus rhamnosus GG-supplemented formula expands butyrate-producing bacterial strains in food allergic infants.

Dietary intervention with extensively hydrolyzed casein formula supplemented with Lactobacillus rhamnosus GG (EHCF+LGG) accelerates tolerance acquisition in infants with cow's milk allergy (CMA). We examined whether this effect is attributable, at least in part, to an influence on the gut microbiota. Fecal samples from healthy controls (n=20) and from CMA infants (n=19) before and after treatment with EHCF with (n=12) and without (n=7) supplementation with LGG were compared by 16S rRNA-based operational taxonomic unit clustering and oligotyping. Differential feature selection and generalized linear model fitting revealed that the CMA infants have a diverse gut microbial community structure dominated by Lachnospiraceae (20.5±9.7%) and Ruminococcaceae (16.2±9.1%). Blautia, Roseburia and Coprococcus were significantly enriched following treatment with EHCF and LGG, but only one genus, Oscillospira, was significantly different between infants that became tolerant and those that remained allergic. However, most tolerant infants showed a significant increase in fecal butyrate levels, and those taxa that were significantly enriched in these samples, Blautia and Roseburia, exhibited specific strain-level demarcations between tolerant and allergic infants. Our data suggest that EHCF+LGG promotes tolerance in infants with CMA, in part, by influencing the strain-level bacterial community structure of the infant gut

Anti-GD2 CAR-NKT cells in patients with relapsed or refractory neuroblastoma: an interim analysis

Vα24-invariant natural killer T (NKT) cells have shown potent anti-tumor properties in murine tumor models and have been linked to favorable outcomes in patients with cancer. However, low numbers of these cells in humans have hindered their clinical applications. Here we report interim results from all three patients enrolled on dose level 1 in a phase 1 dose-escalation trial of autologous NKT cells engineered to co-express a GD2-specific chimeric antigen receptor (CAR) with interleukin-15 in children with relapsed or resistant neuroblastoma (NCT03294954). Primary and secondary objectives were to assess safety and anti-tumor responses, respectively, with immune response evaluation as an additional objective. We ex vivo expanded highly pure NKT cells (mean ± s.d., 94.7 ± 3.8%) and treated patients with 3 × 106 CAR-NKT cells per square meter of body surface area after lymphodepleting conditioning with cyclophosphamide/fludarabine (Cy/Flu). Cy/Flu conditioning was the probable cause for grade 3–4 hematologic adverse events, as they occurred before CAR-NKT cell infusion, and no dose-limiting toxicities were observed. CAR-NKT cells expanded in vivo, localized to tumors and, in one patient, induced an objective response with regression of bone metastatic lesions. These initial results suggest that CAR-NKT cells can be expanded to clinical scale and safely applied to treat patients with cancer

Temporal bacterial and metabolic development of the preterm gut reveals specific signatures in health and disease

Background - The preterm microbiome is crucial to gut health and may contribute to necrotising enterocolitis (NEC), which represents the most significant pathology affecting preterm infants. From a cohort of 318 infants, <32 weeks gestation, we selected 7 infants who developed NEC (defined rigorously) and 28 matched controls. We performed detailed temporal bacterial (n = 641) and metabolomic (n = 75) profiling of the gut microbiome throughout the disease. Results - A core community of Klebsiella, Escherichia, Staphyloccocus, and Enterococcus was present in all samples. Gut microbiota profiles grouped into six distinct clusters, termed preterm gut community types (PGCTs). Each PGCT reflected dominance by the core operational taxonomic units (OTUs), except of PGCT 6, which had high diversity and was dominant in bifidobacteria. While PGCTs 1–5 were present in infants prior to NEC diagnosis, PGCT 6 was comprised exclusively of healthy samples. NEC infants had significantly more PGCT transitions prior to diagnosis. Metabolomic profiling identified significant pathways associated with NEC onset, with metabolites involved in linoleate metabolism significantly associated with NEC diagnosis. Notably, metabolites associated with NEC were the lowest in PGCT 6. Conclusions - This is the first study to integrate sequence and metabolomic stool analysis in preterm neonates, demonstrating that NEC does not have a uniform microbial signature. However, a diverse gut microbiome with a high abundance of bifidobacteria may protect preterm infants from disease. These results may inform biomarker development and improve understanding of gut-mediated mechanisms of NEC

The Influence of Transcription Factor Competition on the Relationship between Occupancy and Affinity

Transcription factors (TFs) are proteins that bind to specific sites on the DNA and regulate gene activity. Identifying where TF molecules bind and how much time they spend on their target sites is key to understanding transcriptional regulation. It is usually assumed that the free energy of binding of a TF to the DNA (the affinity of the site) is highly correlated to the amount of time the TF remains bound (the occupancy of the site). However, knowing the binding energy is not sufficient to infer actual binding site occupancy. This mismatch between the occupancy predicted by the affinity and the observed occupancy may be caused by various factors, such as TF abundance, competition between TFs or the arrangement of the sites on the DNA. We investigated the relationship between the affinity of a TF for a set of binding sites and their occupancy. In particular, we considered the case of the transcription factor lac repressor (lacI) in E.coli, and performed stochastic simulations of the TF dynamics on the DNA for various combinations of lacI abundance and competing TFs that contribute to macromolecular crowding. We also investigated the relationship of site occupancy and the information content of position weight matrices (PWMs) used to represent binding sites. Our results showed that for medium and high affinity sites, TF competition does not play a significant role for genomic occupancy except in cases when the abundance of the TF is significantly increased, or when the PWM displays relatively low information content. Nevertheless, for medium and low affinity sites, an increase in TF abundance (for both cognate and non-cognate molecules) leads to an increase in occupancy at several sites. © 2013 Zabet et al

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

Modeling Reveals Bistability and Low-Pass Filtering in the Network Module Determining Blood Stem Cell Fate

Combinatorial regulation of gene expression is ubiquitous in eukaryotes with multiple inputs converging on regulatory control elements. The dynamic properties of these elements determine the functionality of genetic networks regulating differentiation and development. Here we propose a method to quantitatively characterize the regulatory output of distant enhancers with a biophysical approach that recursively determines free energies of protein-protein and protein-DNA interactions from experimental analysis of transcriptional reporter libraries. We apply this method to model the Scl-Gata2-Fli1 triad—a network module important for cell fate specification of hematopoietic stem cells. We show that this triad module is inherently bistable with irreversible transitions in response to physiologically relevant signals such as Notch, Bmp4 and Gata1 and we use the model to predict the sensitivity of the network to mutations. We also show that the triad acts as a low-pass filter by switching between steady states only in response to signals that persist for longer than a minimum duration threshold. We have found that the auto-regulation loops connecting the slow-degrading Scl to Gata2 and Fli1 are crucial for this low-pass filtering property. Taken together our analysis not only reveals new insights into hematopoietic stem cell regulatory network functionality but also provides a novel and widely applicable strategy to incorporate experimental measurements into dynamical network models

Quantitative Models of the Mechanisms That Control Genome-Wide Patterns of Transcription Factor Binding during Early Drosophila Development

Transcription factors that drive complex patterns of gene expression during animal development bind to thousands of genomic regions, with quantitative differences in binding across bound regions mediating their activity. While we now have tools to characterize the DNA affinities of these proteins and to precisely measure their genome-wide distribution in vivo, our understanding of the forces that determine where, when, and to what extent they bind remains primitive. Here we use a thermodynamic model of transcription factor binding to evaluate the contribution of different biophysical forces to the binding of five regulators of early embryonic anterior-posterior patterning in Drosophila melanogaster. Predictions based on DNA sequence and in vitro protein-DNA affinities alone achieve a correlation of ∼0.4 with experimental measurements of in vivo binding. Incorporating cooperativity and competition among the five factors, and accounting for spatial patterning by modeling binding in every nucleus independently, had little effect on prediction accuracy. A major source of error was the prediction of binding events that do not occur in vivo, which we hypothesized reflected reduced accessibility of chromatin. To test this, we incorporated experimental measurements of genome-wide DNA accessibility into our model, effectively restricting predicted binding to regions of open chromatin. This dramatically improved our predictions to a correlation of 0.6–0.9 for various factors across known target genes. Finally, we used our model to quantify the roles of DNA sequence, accessibility, and binding competition and cooperativity. Our results show that, in regions of open chromatin, binding can be predicted almost exclusively by the sequence specificity of individual factors, with a minimal role for protein interactions. We suggest that a combination of experimentally determined chromatin accessibility data and simple computational models of transcription factor binding may be used to predict the binding landscape of any animal transcription factor with significant precision

Current understanding of the human microbiome

Author Posting. © The Author(s), 2018. This is the author's version of the work. It is posted here by permission of Nature Publishing Group for personal use, not for redistribution. The definitive version was published in Nature Medicine 24 (2018): 392–400, doi:10.1038/nm.4517.Our understanding of the link between the human microbiome and disease, including obesity, inflammatory bowel disease, arthritis and autism, is rapidly expanding. Improvements in the throughput and accuracy of DNA sequencing of the genomes of microbial communities associated with human samples, complemented by analysis of transcriptomes, proteomes, metabolomes and immunomes, and mechanistic experiments in model systems, have vastly improved our ability to understand the structure and function of the microbiome in both diseased and healthy states. However, many challenges remain. In this Review, we focus on studies in humans to describe these challenges, and propose strategies that leverage existing knowledge to move rapidly from correlation to causation, and ultimately to translation.Many of the studies described here in our laboratories were supported by the NIH, NSF, DOE, and the Alfred P. Sloan Foundation.2018-10-1