The Transcriptional Response of Soil Bacteria to Long-Term Warming and Short-Term Seasonal Fluctuations in a Terrestrial Forest.

Terrestrial ecosystems are an important carbon store, and this carbon is vulnerable to microbial degradation with climate warming. After 30 years of experimental warming, carbon stocks in a temperate mixed deciduous forest were observed to be reduced by 30% in the heated plots relative to the controls. In addition, soil respiration was seasonal, as was the warming treatment effect. We therefore hypothesized that long-term warming will have higher expressions of genes related to carbohydrate and lipid metabolism due to increased utilization of recalcitrant carbon pools compared to controls. Because of the seasonal effect of soil respiration and the warming treatment, we further hypothesized that these patterns will be seasonal. We used RNA sequencing to show how the microbial community responds to long-term warming (~30 years) in Harvard Forest, MA. Total RNA was extracted from mineral and organic soil types from two treatment plots (+5°C heated and ambient control), at two time points (June and October) and sequenced using Illumina NextSeq technology. Treatment had a larger effect size on KEGG annotated transcripts than on CAZymes, while soil types more strongly affected CAZymes than KEGG annotated transcripts, though effect sizes overall were small. Although, warming showed a small effect on overall CAZymes expression, several carbohydrate-associated enzymes showed increased expression in heated soils (~68% of all differentially expressed transcripts). Further, exploratory analysis using an unconstrained method showed increased abundances of enzymes related to polysaccharide and lipid metabolism and decomposition in heated soils. Compared to long-term warming, we detected a relatively small effect of seasonal variation on community gene expression. Together, these results indicate that the higher carbohydrate degrading potential of bacteria in heated plots can possibly accelerate a self-reinforcing carbon cycle-temperature feedback in a warming climate.

Reclassification of the Clostridium clostridioforme and Clostridium sphenoides clades as Enterocloster gen. nov. and Lacrimispora gen. nov., including reclassification of 15 taxa.

The genus Clostridium belongs to the family Clostridiaceae. However, many species with the genus name Clostridium are found in different families and even crossing into a different phylum. Motivated by recently completed genome sequences, we propose the reclassification of two separate clades that include misclassified Clostridium species which phylogenetically lie within the family Lachnospiraceae, known for being benign members of gut microbiomes and for their plant-degrading capabilities. We use several phylogenetic and phylogenomic perspectives as well as phenotypic comparisons to gain insight into the evolutionary history of these taxa. One clade, which includes Clostridium clostridioforme, Clostridium aldenense, Clostridium asparagiforme, Clostridium bolteae, Clostridium citroniae and Clostridium lavalense, we propose to reclassify as Enterocloster gen. nov., and reclassify the species as Enterocloster clostridioformis comb. nov., Enterocloster aldensis comb. nov., Enterocloster asparagiformis comb. nov., Enterocloster bolteae comb. nov., Enterocloster citroniae comb. nov. and Enterocloster lavalensis comb. nov. The other clade comprises Clostridium sphenoides, Clostridium aerotolerans, Clostridium algidixylanolyticum, Clostridium amygdalinum, Clostridium celerecrescens, Clostridium indolis, Clostridium saccharolyticum, Clostridium xylanolyticum and Desulfotomaculum guttoideum, and we propose to reclassify it as Lacrimispora gen. nov., including reclassification of the members as Lacrimispora sphenoides comb. nov., Lacrimispora aerotolerans comb. nov., Lacrimispora algidixylanolytica comb. nov., Lacrimispora amygdalina comb. nov., Lacrimispora celerecrescens comb. nov., Lacrimispora indolis comb. nov., Lacrimispora saccharolytica comb. nov. and Lacrimispora xylanolytica comb. nov. We emend the description of D. guttoideum to reflect that it is a later heterotypic synonym of Clostridiums phenoides, which we have reclassified as Lacrimispora sphenoides.

Kineothrix alysoides, gen. nov., sp. nov., a saccharolytic butyrate-producer within the family Lachnospiraceae.

An anaerobic, saccharolytic, spore-forming, butyrate-producing bacterium, strain KNHs209T, was isolated from a switchgrass microcosm seeded with forest soil. Cells were highly motile rods, often forming long filamentous chains which were easily observed moving under the microscope. Its closest phylogenetic relative was Eisenbergiella tayi (16S rRNA gene sequence identity 94.2 %), although it was easily distinguishable based on its morphology and physiology. Whole-genome sequencing enabled development of a minimal medium, and also suggested that the organism is capable of fixing nitrogen. Its wide variety of growth substrates was mirrored by a high number of encoded chemotaxis receptors (45, the highest in the family Lachnospiraceae). Strain KNHs209T utilized a wide variety of carbohydrates, but not cellulose or xylan. Fermentation products included formate, acetate and butyrate; sulfur compounds and nitrate were not reduced. Strain KNHs209T grew optimally at 35-40 °C and pH 7. The genomic DNA G+C content was 42.74 mol%; the major membrane fatty acids were C14 : 0 and C16 : 0. Based on phenotypic, genomic, phylogenetic and chemotaxonomic analyses, this organism represents a novel genus and species within the family Lachnospiraceae for which the name Kineothrix alysoides, gen. nov., sp. nov. is proposed. The type strain is KNHs209T (=ATCC TSD-26T=DSM 100556T).

Tempo and mode of genome evolution in a 50,000-generation experiment.

Adaptation by natural selection depends on the rates, effects and interactions of many mutations, making it difficult to determine what proportion of mutations in an evolving lineage are beneficial. Here we analysed 264 complete genomes from 12 Escherichia coli populations to characterize their dynamics over 50,000 generations. The populations that retained the ancestral mutation rate support a model in which most fixed mutations are beneficial, the fraction of beneficial mutations declines as fitness rises, and neutral mutations accumulate at a constant rate. We also compared these populations to mutation-accumulation lines evolved under a bottlenecking regime that minimizes selection. Nonsynonymous mutations, intergenic mutations, insertions and deletions are overrepresented in the long-term populations, further supporting the inference that most mutations that reached high frequency were favoured by selection. These results illuminate the shifting balance of forces that govern genome evolution in populations adapting to a new environment.

Genome and Transcriptome of Clostridium phytofermentans, Catalyst for the Direct Conversion of Plant Feedstocks to Fuels.

Clostridium phytofermentans was isolated from forest soil and is distinguished by its capacity to directly ferment plant cell wall polysaccharides into ethanol as the primary product, suggesting that it possesses unusual catabolic pathways. The objective of the present study was to understand the molecular mechanisms of biomass conversion to ethanol in a single organism, Clostridium phytofermentans, by analyzing its complete genome and transcriptome during growth on plant carbohydrates. The saccharolytic versatility of C. phytofermentans is reflected in a diversity of genes encoding ATP-binding cassette sugar transporters and glycoside hydrolases, many of which may have been acquired through horizontal gene transfer. These genes are frequently organized as operons that may be controlled individually by the many transcriptional regulators identified in the genome. Preferential ethanol production may be due to high levels of expression of multiple ethanol dehydrogenases and additional pathways maximizing ethanol yield. The genome also encodes three different proteinaceous bacterial microcompartments with the capacity to compartmentalize pathways that divert fermentation intermediates to various products. These characteristics make C. phytofermentans an attractive resource for improving the efficiency and speed of biomass conversion to biofuels.

Long-term forest soil warming alters microbial communities in temperate forest soils.

Soil microbes are major drivers of soil carbon cycling, yet we lack an understanding of how climate warming will affect microbial communities. Three ongoing field studies at the Harvard Forest Long-term Ecological Research (LTER) site (Petersham, MA) have warmed soils 5°C above ambient temperatures for 5, 8, and 20 years. We used this chronosequence to test the hypothesis that soil microbial communities have changed in response to chronic warming. Bacterial community composition was studied using Illumina sequencing of the 16S ribosomal RNA gene, and bacterial and fungal abundance were assessed using quantitative PCR. Only the 20-year warmed site exhibited significant change in bacterial community structure in the organic soil horizon, with no significant changes in the mineral soil. The dominant taxa, abundant at 0.1% or greater, represented 0.3% of the richness but nearly 50% of the observations (sequences). Individual members of the Actinobacteria, Alphaproteobacteria and Acidobacteria showed strong warming responses, with one Actinomycete decreasing from 4.5 to 1% relative abundance with warming. Ribosomal RNA copy number can obfuscate community profiles, but is also correlated with maximum growth rate or trophic strategy among bacteria. Ribosomal RNA copy number correction did not affect community profiles, but rRNA copy number was significantly decreased in warming plots compared to controls. Increased bacterial evenness, shifting beta diversity, decreased fungal abundance and increased abundance of bacteria with low rRNA operon copy number, including Alphaproteobacteria and Acidobacteria, together suggest that more or alternative niche space is being created over the course of long-term warming.

The complete genome sequence of Clostridium indolis DSM 755(T.).

Clostridium indolis DSM 755(T) is a bacterium commonly found in soils and the feces of birds and mammals. Despite its prevalence, little is known about the ecology or physiology of this species. However, close relatives, C. saccharolyticum and C. hathewayi, have demonstrated interesting metabolic potentials related to plant degradation and human health. The genome of C. indolis DSM 755(T) reveals an abundance of genes in functional groups associated with the transport and utilization of carbohydrates, as well as citrate, lactate, and aromatics. Ecologically relevant gene clusters related to nitrogen fixation and a unique type of bacterial microcompartment, the CoAT BMC, are also detected. Our genome analysis suggests hypotheses to be tested in future culture based work to better understand the physiology of this poorly described species.

Strong genome-wide selection early in the evolution of Prochlorococcus resulted in a reduced genome through the loss of a large number of small effect genes.

The smallest genomes of any photosynthetic organisms are found in a group of free-living marine cyanobacteria, Prochlorococcus. To determine the underlying evolutionary mechanisms, we developed a new method to reconstruct the steps leading to the Prochlorococcus genome reduction using 12 Prochlorococcus and 6 marine Synechococcus genomes. Our results reveal that small genome sizes within Prochlorococcus were largely determined shortly after the split of Prochlorococcus and Synechococcus (an early big shrink) and thus for the first time decouple the genome reduction from Prochlorococcus diversification. A maximum likelihood approach was then used to estimate changes of nucleotide substitution rate and selection strength along Prochlorococcus evolution in a phylogenetic framework. Strong genome wide purifying selection was associated with the loss of many genes in the early evolutionary stage. The deleted genes were distributed around the genome, participated in many different functional categories and in general had been under relaxed selection pressure. We propose that shortly after Prochlorococcus diverged from its common ancestor with marine Synechococcus, its population size increased quickly thus increasing efficacy of selection. Due to limited nutrients and a relatively constant environment, selection favored a streamlined genome for maximum economy. Strong genome wide selection subsequently caused the loss of genes with small functional effect including the loss of some DNA repair genes. In summary, genome reduction in Prochlorococcus resulted in genome features that are similar to symbiotic bacteria and pathogens, however, the small genome sizes resulted from an increase in genome wide selection rather than a consequence of a reduced ecological niche or relaxed selection due to genetic drift.

Population level analysis of evolved mutations underlying improvements in plant hemicellulose and cellulose fermentation by Clostridium phytofermentans.

BACKGROUND: The complexity of plant cell walls creates many challenges for microbial decomposition. Clostridium phytofermentans, an anaerobic bacterium isolated from forest soil, directly breaks down and utilizes many plant cell wall carbohydrates. The objective of this research is to understand constraints on rates of plant decomposition by Clostridium phytofermentans and identify molecular mechanisms that may overcome these limitations. RESULTS: Experimental evolution via repeated serial transfers during exponential growth was used to select for C. phytofermentans genotypes that grow more rapidly on cellobiose, cellulose and xylan. To identify the underlying mutations an average of 13,600,000 paired-end reads were generated per population resulting in ∼300 fold coverage of each site in the genome. Mutations with allele frequencies of 5% or greater could be identified with statistical confidence. Many mutations are in carbohydrate-related genes including the promoter regions of glycoside hydrolases and amino acid substitutions in ABC transport proteins involved in carbohydrate uptake, signal transduction sensors that detect specific carbohydrates, proteins that affect the export of extracellular enzymes, and regulators of unknown specificity. Structural modeling of the ABC transporter complex proteins suggests that mutations in these genes may alter the recognition of carbohydrates by substrate-binding proteins and communication between the intercellular face of the transmembrane and the ATPase binding proteins. CONCLUSIONS: Experimental evolution was effective in identifying molecular constraints on the rate of hemicellulose and cellulose fermentation and selected for putative gain of function mutations that do not typically appear in traditional molecular genetic screens. The results reveal new strategies for evolving and engineering microorganisms for faster growth on plant carbohydrates.

An in vitro model of the horse gut microbiome enables identification of lactate-utilizing bacteria that differentially respond to starch induction.

Laminitis is a chronic, crippling disease triggered by the sudden influx of dietary starch. Starch reaches the hindgut resulting in enrichment of lactic acid bacteria, lactate accumulation, and acidification of the gut contents. Bacterial products enter the bloodstream and precipitate systemic inflammation. Hindgut lactate levels are normally low because specific bacterial groups convert lactate to short chain fatty acids. Why this mechanism fails when lactate levels rapidly rise, and why some hindgut communities can recover is unknown. Fecal samples from three adult horses eating identical diets provided bacterial communities for this in vitro study. Triplicate microcosms of fecal slurries were enriched with lactate and/or starch. Metabolic products (short chain fatty acids, headspace gases, and hydrogen sulfide) were measured and microbial community compositions determined using Illumina 16S rRNA sequencing over 12-hour intervals. We report that patterns of change in short chain fatty acid levels and pH in our in vitro system are similar to those seen in in vivo laminitis induction models. Community differences between microcosms with disparate abilities to clear excess lactate suggest profiles conferring resistance of starch-induction conditions. Where lactate levels recover following starch induction conditions, propionate and acetate levels rise correspondingly and taxa related to Megasphaeraelsdenii reach levels exceeding 70% relative abundance. In lactate and control cultures, taxa related to Veillonellamontpellierensis are enriched as lactate levels fall. Understanding these community differences and factors promoting the growth of specific lactate utilizing taxa may be useful to prevent acidosis under starch-induction conditions.

Involvement of a bacterial microcompartment in the metabolism of fucose and rhamnose by Clostridium phytofermentans.

BACKGROUND: Clostridium phytofermentans, an anaerobic soil bacterium, can directly convert plant biomass into biofuels. The genome of C. phytofermentans contains three loci with genes encoding shell proteins of bacterial microcompartments (BMC), organelles composed entirely of proteins. METHODOLOGY AND PRINCIPAL FINDINGS: One of the BMC loci has homology to a BMC-encoding locus implicated in the conversion of fucose to propanol and propionate in a human gut commensal, Roseburia inulinivorans. We hypothesized that it had a similar role in C. phytofermentans. When C. phytofermentans was grown on fucose, the major products identified were ethanol, propanol and propionate. Transmission electron microscopy of fucose- and rhamnose-grown cultures revealed polyhedral structures, presumably BMCs. Microarray analysis indicated that during growth on fucose, operons coding for the BMC locus, fucose dissimilatory enzymes, and an ATP-binding cassette transporter became the dominant transcripts. These data are consistent with fucose fermentation producing a 1,2-propanediol intermediate that is further metabolized in the microcompartment encoded in the BMC locus. Growth on another deoxyhexose sugar, rhamnose, resulted in the expression of the same BMC locus and similar fermentation products. However, a different set of dissimilatory enzymes and transport system genes were induced. Quite surprisingly, growth on fucose or rhamnose also led to the expression of a diverse array of complex plant polysaccharide-degrading enzymes. CONCLUSIONS/SIGNIFICANCE: Based on physiological, genomic, and microarray analyses, we propose a model for the fermentation of fucose and rhamnose in C. phytofermentans that includes enzymes encoded in the same BMC locus. Comparative genomic analysis suggests that this BMC may be present in other clostridial species.

Environmental sequence data from the Sargasso Sea reveal that the characteristics of genome reduction in Prochlorococcus are not a harbinger for an escalation in genetic drift.

The marine cyanobacterium Prochlorococcus MED4 has the smallest sequenced genome of any photosynthetic organism. Prochlorococcus MED4 shares many genomic characteristics with chloroplasts and bacterial endosymbionts, including a reduced coding capacity, missing DNA repair genes, a minimal transcriptional regulatory network, a marked AT% bias, and an accelerated rate of amino acid changes. In chloroplasts and endosymbionts, these molecular phenotypes appear to be symptomatic of a relative increase in genetic drift due to restrictions on effective population size in the host environment. As a free-living bacterium, Prochlorococcus MED4 is not known to be subject to similar ecological constraints. To test whether the high-light-adapted Prochlorococcus MED4 is experiencing a reduction in selection efficiency resulting from genetic drift, we examine two data sets, namely, the environmental genome shotgun sequencing data from the Sargasso Sea and a set of cyanobacterial genome sequences. After integrating these data sets, we compare the evolutionary profile of a high-light Prochlorococcus group to that of a group of Synechococcus (a closely related group of marine cyanobacteria) that does not exhibit a similar small-genome syndrome. The average pairwise dN/dS ratios in the high-light-adapted Prochlorococcus group are significantly lower than those in the Synechococcus group, leading us to reject the hypothesis that the Prochlorococcus group is currently experiencing higher levels of genetic drift.

Transcriptional responses to estrogen and progesterone in mammary gland identify networks regulating p53 activity.

Estrogen and progestins are essential for mammary growth and differentiation but also enhance the activity of the p53 tumor suppressor protein in the mammary epithelium. However, the pathways by which these hormones regulate p53 activity are unknown. Microarrays were used to profile the transcriptional changes within the mammary gland after administration of either vehicle, 17beta-estradiol (E), or progesterone (P) individually and combined (EP). Treatment with EP yielded 1182 unique genes that were differentially expressed compared to the vehicle-treated group. Although 30% of genes were responsive to either E or P individually, combined treatment with both EP had a synergistic effect accounting for 60% of the differentially regulated genes. Analysis of protein-protein interactions identified p53, RelA, Snw1, and Igfals as common targets of genes regulated by EP. RelA and p53 form hubs within a network connected by genes that are regulated by EP and that may coordinate the competing functions of RelA and p53 in proliferation and survival of cells. Induction of early growth response 1 (Egr1) and Stratifin (Sfn) (also known as 14-3-3sigma) by EP was confirmed by reverse transcription-quantitative PCR and shown to be p53 independent. In luciferase reporter assays, Egr1 was shown to enhance transcriptional activation by p53 and inhibit nuclear factor kappaB activity. These results identify a gene expression network that provides redundant activation of RelA to support proliferation as well as sensitize p53 to ensure proper surveillance and integration of their competing functions through factors such as Egr1, which both enhance p53 and inhibit RelA.

Rapid changes in gene expression dynamics in response to superoxide reveal SoxRS-dependent and independent transcriptional networks.

BACKGROUND: SoxR and SoxS constitute an intracellular signal response system that rapidly detects changes in superoxide levels and modulates gene expression in E. coli. A time series microarray design was used to identify co-regulated SoxRS-dependent and independent genes modulated by superoxide minutes after exposure to stress. METHODOLOGY/PRINCIPAL FINDINGS: soxS mRNA levels surged to near maximal levels within the first few minutes of exposure to paraquat, a superoxide-producing compound, followed by a rise in mRNA levels of known SoxS-regulated genes. Based on a new method for determining the biological significance of clustering results, a total of 138 genic regions, including several transcription factors and putative sRNAs were identified as being regulated through the SoxRS signaling pathway within 10 minutes of paraquat treatment. A statistically significant two-block SoxS motif was identified through analysis of the SoxS-regulated genes. The SoxRS-independent response included members of the OxyR, CysB, IscR, BirA and Fur regulons. Finally, the relative sensitivity to superoxide was measured in 94 strains carrying deletions in individual, superoxide-regulated genes. CONCLUSIONS/SIGNIFICANCE: By integrating our microarray time series results with other microarray data, E. coli databases and the primary literature, we propose a model of the primary transcriptional response containing 226 protein-coding and sRNA sequences. From the SoxS dependent network the first statistically significant SoxS-related motif was identified.

Evolutionary origins of the eukaryotic shikimate pathway: gene fusions, horizontal gene transfer, and endosymbiotic replacements.

Currently the shikimate pathway is reported as a metabolic feature of prokaryotes, ascomycete fungi, apicomplexans, and plants. The plant shikimate pathway enzymes have similarities to prokaryote homologues and are largely active in chloroplasts, suggesting ancestry from the plastid progenitor genome. Toxoplasma gondii, which also possesses an alga-derived plastid organelle, encodes a shikimate pathway with similarities to ascomycete genes, including a five-enzyme pentafunctional arom. These data suggests that the shikimate pathway and the pentafunctional arom either had an ancient origin in the eukaryotes or was conveyed by eukaryote-to-eukaryote horizontal gene transfer (HGT). We expand sampling and analyses of the shikimate pathway genes to include the oomycetes, ciliates, diatoms, basidiomycetes, zygomycetes, and the green and red algae. Sequencing of cDNA from Tetrahymena thermophila confirmed the presence of a pentafused arom, as in fungi and T. gondii. Phylogenies and taxon distribution suggest that the arom gene fusion event may be an ancient eukaryotic innovation. Conversely, the Plantae lineage (represented here by both Viridaeplantae and the red algae) acquired different prokaryotic genes for all seven steps of the shikimate pathway. Two of the phylogenies suggest a derivation of the Plantae genes from the cyanobacterial plastid progenitor genome, but if the full Plantae pathway was originally of cyanobacterial origin, then the five other shikimate pathway genes were obtained from a minimum of two other eubacterial genomes. Thus, the phylogenies demonstrate both separate HGTs and shared derived HGTs within the Plantae clade either by primary HGT transfer or secondarily via the plastid progenitor genome. The shared derived characters support the holophyly of the Plantae lineage and a single ancestral primary plastid endosymbiosis. Our analyses also pinpoints a minimum of 50 gene/domain loss events, demonstrating that loss and replacement events have been an important process in eukaryote genome evolution.

Perspective: Evolution and detection of genetic robustness.

Robustness is the invariance of phenotypes in the face of perturbation. The robustness of phenotypes appears at various levels of biological organization, including gene expression, protein folding, metabolic flux, physiological homeostasis, development, and even organismal fitness. The mechanisms underlying robustness are diverse, ranging from thermodynamic stability at the RNA and protein level to behavior at the organismal level. Phenotypes can be robust either against heritable perturbations (e.g., mutations) or nonheritable perturbations (e.g., the weather). Here we primarily focus on the first kind of robustness--genetic robustness--and survey three growing avenues of research: (1) measuring genetic robustness in nature and in the laboratory; (2) understanding the evolution of genetic robustness: and (3) exploring the implications of genetic robustness for future evolution.

Evidence that plant-like genes in Chlamydia species reflect an ancestral relationship between Chlamydiaceae, cyanobacteria, and the chloroplast.

An unusually high proportion of proteins encoded in Chlamydia genomes are most similar to plant proteins, leading to proposals that a Chlamydia ancestor obtained genes from a plant or plant-like host organism by horizontal gene transfer. However, during an analysis of bacterial-eukaryotic protein similarities, we found that the vast majority of plant-like sequences in Chlamydia are most similar to plant proteins that are targeted to the chloroplast, an organelle derived from a cyanobacterium. We present further evidence suggesting that plant-like genes in Chlamydia, and other Chlamydiaceae, are likely a reflection of an unappreciated evolutionary relationship between the Chlamydiaceae and the cyanobacteria-chloroplast lineage. Further analyses of bacterial and eukaryotic genomes indicates the importance of evaluating organellar ancestry of eukaryotic proteins when identifying bacteria-eukaryote homologs or horizontal gene transfer and supports the proposal that Chlamydiaceae, which are obligate intracellular bacterial pathogens of animals, are not likely exchanging DNA with their hosts.