A library of human gut bacterial isolates paired with longitudinal multiomics data enables mechanistic microbiome research.

Our understanding of how the gut microbiome interacts with its human host has been restrained by limited access to longitudinal datasets to examine stability and dynamics, and by having only a few isolates to test mechanistic hypotheses. Here, we present the Broad Institute-OpenBiome Microbiome Library (BIO-ML), a comprehensive collection of 7,758 gut bacterial isolates paired with 3,632 genome sequences and longitudinal multi-omics data. We show that microbial species maintain stable population sizes within and across humans and that commonly used 'omics' survey methods are more reliable when using averages over multiple days of sampling. Variation of gut metabolites within people over time is associated with amino acid levels, and differences across people are associated with differences in bile acids. Finally, we show that genomic diversification can be used to infer eco-evolutionary dynamics and in vivo selection pressures for strains within individuals. The BIO-ML is a unique resource designed to enable hypothesis-driven microbiome research.

Mobile genes in the human microbiome are structured from global to individual scales.

Recent work has underscored the importance of the microbiome in human health, and has largely attributed differences in phenotype to differences in the species present among individuals. However, mobile genes can confer profoundly different phenotypes on different strains of the same species. Little is known about the function and distribution of mobile genes in the human microbiome, and in particular whether the gene pool is globally homogenous or constrained by human population structure. Here, we investigate this question by comparing the mobile genes found in the microbiomes of 81 metropolitan North Americans with those of 172 agrarian Fiji islanders using a combination of single-cell genomics and metagenomics. We find large differences in mobile gene content between the Fijian and North American microbiomes, with functional variation that mirrors known dietary differences such as the excess of plant-based starch degradation genes found in Fijian individuals. Notably, we also observed differences between the mobile gene pools of neighbouring Fijian villages, even though microbiome composition across villages is similar. Finally, we observe high rates of recombination leading to individual-specific mobile elements, suggesting that the abundance of some genes may reflect environmental selection rather than dispersal limitation. Together, these data support the hypothesis that human activities and behaviours provide selective pressures that shape mobile gene pools, and that acquisition of mobile genes is important for colonizing specific human populations.

Clostridium difficile associated risk of death score (CARDS): a novel severity score to predict mortality among hospitalised patients with C. difficile infection.

BACKGROUND: Clostridium difficile infection (CDI) is a public health threat and associated with significant mortality. However, there is a paucity of objectively derived CDI severity scoring systems to predict mortality. AIM: To develop a novel CDI risk score to predict mortality entitled: Clostridium difficile associated risk of death score (CARDS). METHODS: We obtained data from the United States 2011 Nationwide Inpatient Sample (NIS) database. All CDI-associated hospitalisations were identified using discharge codes (ICD-9-CM, 008.45). Multivariate logistic regression was utilised to identify independent predictors of mortality. Clostridium difficile associated risk of death score was calculated by assigning a numeric weight to each parameter based on their odds ratio in the final logistic model. Predictive properties of model discrimination were assessed using the c-statistic and validated in an independent sample using the 2010 NIS database. RESULTS: We identified 77 776 hospitalisations, yielding an estimate of 374 747 cases with an associated diagnosis of CDI in the US, 8% of whom died in the hospital. The eight severity score predictors were identified on multivariate analysis: age, cardiopulmonary disease, malignancy, diabetes, inflammatory bowel disease, acute renal failure, liver disease and ICU admission, with weights ranging from -1 (for diabetes) to 5 (for ICU admission). The overall risk score in the cohort ranged from 0 to 18. Mortality increased significantly as CARDS increased. CDI-associated mortality was 1.2% with a CARDS of 0 compared to 100% with CARDS of 18. The model performed equally well in our validation cohort. CONCLUSION: Clostridium difficile associated risk of death score is a promising simple severity score to predict mortality among those hospitalised with C. difficile infection.

Ancestral Reconstruction of a Pre-LUCA Aminoacyl-tRNA Synthetase Ancestor Supports the Late Addition of Trp to the Genetic Code.

The genetic code was likely complete in its current form by the time of the last universal common ancestor (LUCA). Several scenarios have been proposed for explaining the code's pre-LUCA emergence and expansion, and the relative order of the appearance of amino acids used in translation. One co-evolutionary model of genetic code expansion proposes that at least some amino acids were added to the code by the ancient divergence of aminoacyl-tRNA synthetase (aaRS) families. Of all the amino acids used within the genetic code, Trp is most frequently claimed as a relatively recent addition. We observe that, since TrpRS and TyrRS are paralogous protein families retaining significant sequence similarity, the inferred sequence composition of their ancestor can be used to evaluate this co-evolutionary model of genetic code expansion. We show that ancestral sequence reconstructions of the pre-LUCA paralog ancestor of TyrRS and TrpRS have several sites containing Tyr, yet a complete absence of sites containing Trp. This is consistent with the paralog ancestor being specific for the utilization of Tyr, with Trp being a subsequent addition to the genetic code facilitated by a process of aaRS divergence and neofunctionalization. Only after this divergence could Trp be specifically encoded and incorporated into proteins, including the TyrRS and TrpRS descendant lineages themselves. This early absence of Trp is observed under both homogeneous and non-homogeneous models of ancestral sequence reconstruction. Simulations support that this observed absence of Trp is unlikely to be due to chance or model bias. These results support that the final stages of genetic code evolution occurred well within the "protein world," and that the presence-absence of Trp within conserved sites of ancient protein domains is a likely measure of their relative antiquity, permitting the relative timing of extremely early events within protein evolution before LUCA.

Temporal transcriptomic analysis as Desulfovibrio vulgaris Hildenborough transitions into stationary phase during electron donor depletion.

Desulfovibrio vulgaris was cultivated in a defined medium, and biomass was sampled for approximately 70 h to characterize the shifts in gene expression as cells transitioned from the exponential to the stationary phase during electron donor depletion. In addition to temporal transcriptomics, total protein, carbohydrate, lactate, acetate, and sulfate levels were measured. The microarray data were examined for statistically significant expression changes, hierarchical cluster analysis, and promoter element prediction and were validated by quantitative PCR. As the cells transitioned from the exponential phase to the stationary phase, a majority of the down-expressed genes were involved in translation and transcription, and this trend continued at the remaining times. There were general increases in relative expression for intracellular trafficking and secretion, ion transport, and coenzyme metabolism as the cells entered the stationary phase. As expected, the DNA replication machinery was down-expressed, and the expression of genes involved in DNA repair increased during the stationary phase. Genes involved in amino acid acquisition, carbohydrate metabolism, energy production, and cell envelope biogenesis did not exhibit uniform transcriptional responses. Interestingly, most phage-related genes were up-expressed at the onset of the stationary phase. This result suggested that nutrient depletion may affect community dynamics and DNA transfer mechanisms of sulfate-reducing bacteria via the phage cycle. The putative feoAB system (in addition to other presumptive iron metabolism genes) was significantly up-expressed, and this suggested the possible importance of Fe2+ acquisition under metal-reducing conditions. The expression of a large subset of carbohydrate-related genes was altered, and the total cellular carbohydrate levels declined during the growth phase transition. Interestingly, the D. vulgaris genome does not contain a putative rpoS gene, a common attribute of the delta-Proteobacteria genomes sequenced to date, and the transcription profiles of other putative rpo genes were not significantly altered. Our results indicated that in addition to expected changes (e.g., energy conversion, protein turnover, translation, transcription, and DNA replication and repair), genes related to phage, stress response, carbohydrate flux, the outer envelope, and iron homeostasis played important roles as D. vulgaris cells experienced electron donor depletion.

Global analysis of heat shock response in Desulfovibrio vulgaris Hildenborough.

Desulfovibrio vulgaris Hildenborough belongs to a class of sulfate-reducing bacteria (SRB) and is found ubiquitously in nature. Given the importance of SRB-mediated reduction for bioremediation of metal ion contaminants, ongoing research on D. vulgaris has been in the direction of elucidating regulatory mechanisms for this organism under a variety of stress conditions. This work presents a global view of this organism's response to elevated growth temperature using whole-cell transcriptomics and proteomics tools. Transcriptional response (1.7-fold change or greater; Z >/= 1.5) ranged from 1,135 genes at 15 min to 1,463 genes at 120 min for a temperature up-shift of 13 degrees C from a growth temperature of 37 degrees C for this organism and suggested both direct and indirect modes of heat sensing. Clusters of orthologous group categories that were significantly affected included posttranslational modifications; protein turnover and chaperones (up-regulated); energy production and conversion (down-regulated), nucleotide transport, metabolism (down-regulated), and translation; ribosomal structure; and biogenesis (down-regulated). Analysis of the genome sequence revealed the presence of features of both negative and positive regulation which included the CIRCE element and promoter sequences corresponding to the alternate sigma factors sigma(32) and sigma(54). While mechanisms of heat shock control for some genes appeared to coincide with those established for Escherichia coli and Bacillus subtilis, the presence of unique control schemes for several other genes was also evident. Analysis of protein expression levels using differential in-gel electrophoresis suggested good agreement with transcriptional profiles of several heat shock proteins, including DnaK (DVU0811), HtpG (DVU2643), HtrA (DVU1468), and AhpC (DVU2247). The proteomics study also suggested the possibility of posttranslational modifications in the chaperones DnaK, AhpC, GroES (DVU1977), and GroEL (DVU1976) and also several periplasmic ABC transporters.

Mechanisms of protein folding.

The strong correlation between protein folding rates and the contact order suggests that folding rates are largely determined by the topology of the native structure. However, for a given topology, there may be several possible low free energy paths to the native state and the path that is chosen (the lowest free energy path) may depend on differences in interaction energies and local free energies of ordering in different parts of the structure. For larger proteins whose folding is assisted by chaperones, such as the Escherichia coli chaperonin GroEL, advances have been made in understanding both the aspects of an unfolded protein that GroEL recognizes and the mode of binding to the chaperonin. The possibility that GroEL can remove non-native proteins from kinetic traps by unfolding them either during polypeptide binding to the chaperonin or during the subsequent ATP-dependent formation of folding-active complexes with the co-chaperonin GroES has also been explored.

NMR characterization of residual structure in the denatured state of protein L.

Triple-resonance NMR experiments were used to assign the (13)C(alpha), (13)C(beta), (15)N and NH resonances for all the residues in the denatured state of a destabilized protein L variant in 2 M guanidine. The chemical shifts of most resonances were very close to their random coil values. Significant deviations were observed for G22, L38 and K39; increasing the denaturant concentration shifted the chemical shifts of these residues towards theory random coil values. Medium-range nuclear Overhauser enhancements were detected in segments corresponding to the turn between the first two strands, the end of the second strand through the turn between the second strand and the helix, and the turn between the helix and the third strand in 3D H(1), N(15)-HSQC-NOESY-HSQC experiments on perdeuterated samples. Longer-range interactions were probed by measuring the paramagnetic relaxation enhancement produced by nitroxide spin labels introduced via cysteine residues at five sites around the molecule. Damped oscillations in the magnitude of the paramagnetic relaxation enhancement as a function of distance along the sequence suggested native-like chain reversals in the same three turn regions. The more extensive interactions within the region corresponding to the first beta-turn than in the region corresponding to the second beta-turn suggests that the asymmetry in the folding reaction evident in previous studies of the protein L folding transition state is already established in the denatured state.