Multisubstrate DNA stable isotope probing reveals guild structure of bacteria that mediate soil carbon cycling.

Soil microorganisms determine the fate of soil organic matter (SOM), and their activities compose a major component of the global carbon (C) cycle. We employed a multisubstrate, DNA-stable isotope probing experiment to track bacterial assimilation of C derived from distinct sources that varied in bioavailability. This approach allowed us to measure microbial contributions to SOM processing by measuring the C assimilation dynamics of diverse microorganisms as they interacted within soil. We identified and tracked 1,286 bacterial taxa that assimilated 13C in an agricultural soil over a period of 48 d. Overall 13C-assimilation dynamics of bacterial taxa, defined by the source and timing of the 13C they assimilated, exhibited low phylogenetic conservation. We identified bacterial guilds composed of taxa that had similar 13C assimilation dynamics. We show that C-source bioavailability explained significant variation in both C mineralization dynamics and guild structure, and that the growth dynamics of bacterial guilds differed significantly in response to C addition. We also demonstrate that the guild structure explains significant variation in the biogeographical distribution of bacteria at continental and global scales. These results suggest that an understanding of in situ growth dynamics is essential for understanding microbial contributions to soil C cycling. We interpret these findings in the context of bacterial life history strategies and their relationship to terrestrial C cycling.

Land use alters bacterial growth dynamics in soil.

Microbial growth and mortality are major determinants of soil carbon cycling. We measured in situ growth dynamics of individual bacterial taxa in cropped and successional soils in response to a resource pulse. We hypothesized that land use imposes selection pressures on growth characteristics. We estimated growth and death for 453 and 73 taxa, respectively. The average generation time was 5.04 ± 6.28 (SD; range 0.7-63.5) days. Lag times were shorter in cultivated than successional soils and resource amendment decreased lag times. Taxa exhibiting the greatest growth response also exhibited the greatest mortality, indicative of boom-and-bust dynamics. We observed a bimodal growth rate distribution, representing fast- and slow-growing clusters. Both clusters grew more rapidly in successional soils, which had more organic matter, than cultivated soils. Resource amendment increased the growth rate of the slower growing but not the faster-growing cluster via a mixture of increased growth rates and species turnover, indicating that competitive dynamics constrain growth rates in situ. These two clusters show that copiotrophic bacteria in soils may be subdivided into different life history groups and that these subgroups respond independently to land use and resource availability.

Metagenomic stable isotope probing reveals bacteriophage participation in soil carbon cycling.

Soil viruses are important components of the carbon (C) cycle, yet we still know little about viral ecology in soils. We added diverse 13 C-labelled carbon sources to soil and we used metagenomic-SIP to detect 13 C assimilation by viruses and their putative bacterial hosts. These data allowed us to link a 13 C-labelled bacteriophage to its 13 C-labelled Streptomyces putative host, and we used qPCR to track the dynamics of the putative host and phage in response to C inputs. Following C addition, putative host numbers increased rapidly for 3 days, and then more gradually, reaching maximal abundance on Day 6. Viral abundance and virus:host ratio increased dramatically over 6 days, and remained high thereafter (8.42 ± 2.94). From Days 6 to 30, virus:host ratio remained high, while putative host numbers declined more than 50%. Putative host populations were 13 C-labelled on Days 3-30, while 13 C-labelling of phage was detected on Days 14 and 30. This dynamic suggests rapid growth and 13 C-labelling of the host fueled by new C inputs, followed by extensive host mortality driven by phage lysis. These findings indicate that the viral shunt promotes microbial turnover in soil following new C inputs, thereby altering microbial community dynamics, and facilitating soil organic matter production.

The RNA virome of echinoderms.

Echinoderms are a phylum of marine invertebrates that include model organisms, keystone species, and animals commercially harvested for seafood. Despite their scientific, ecological, and economic importance, there is little known about the diversity of RNA viruses that infect echinoderms compared to other invertebrates. We screened over 900 transcriptomes and viral metagenomes to characterize the RNA virome of 38 echinoderm species from all five classes (Crinoidea, Holothuroidea, Asteroidea, Ophiuroidea and Echinoidea). We identified 347 viral genome fragments that were classified to genera and families within nine viral orders - Picornavirales, Durnavirales, Martellivirales, Nodamuvirales, Reovirales, Amarillovirales, Ghabrivirales, Mononegavirales, and Hepelivirales. We compared the relative viral representation across three life stages (embryo, larvae, adult) and characterized the gene content of contigs which encoded complete or near-complete genomes. The proportion of viral reads in a given transcriptome was not found to significantly differ between life stages though the majority of viral contigs were discovered from transcriptomes of adult tissue. This study illuminates the biodiversity of RNA viruses from echinoderms, revealing the occurrence of viral groups in natural populations.

Conditional filamentation as an adaptive trait of bacteria and its ecological significance in soils.

Bacteria can regulate cell morphology in response to environmental conditions, altering their physiological and metabolic characteristics to improve survival. Conditional filamentation, in which cells suspend division while continuing lateral growth, is a strategy with a range of adaptive benefits. Here, we review the causes and consequences of conditional filamentation with respect to bacterial physiology, ecology and evolution. We describe four major benefits from conditional filamentation: stress tolerance, surface colonization, gradient spanning and the facilitation of biotic interactions. Adopting a filamentous growth habit involves fitness trade-offs which are also examined. We focus on the role of conditional filamentation in soil habitats, where filamentous morphotypes are highly prevalent and where environmental heterogeneity can benefit a conditional response. To illustrate the use of information presented in our review, we tested the conditions regulating filamentation by the forest soil isolate Paraburkholderia elongata 5NT . Filamentation by P. elongata was induced at elevated phosphate concentrations, and was associated with the accumulation of intracellular polyphosphate, highlighting the role of filamentation in a phosphate-solubilizing bacterium. Conditional filamentation enables bacteria to optimize their growth and metabolism in environments that are highly variable, a trait that can impact succession, symbioses, and biogeochemistry in soil environments.

Bacterial community dynamics explain carbon mineralization and assimilation in soils of different land-use history.

Soil dwelling microorganisms are key players in the terrestrial carbon cycle, driving both the degradation and stabilization of soil organic matter. Bacterial community structure and function vary with respect to land use; yet the ecological drivers of this variation remain poorly described and difficult to predict. We conducted a multi-substrate DNA-stable isotope probing experiment across cropland, old-field, and forest habitats to link carbon mineralization dynamics with the dynamics of bacterial growth and carbon assimilation. We tracked the movement of 13 C derived from five distinct carbon sources as it was assimilated into bacterial DNA over time. We show that carbon mineralization, community composition, and carbon assimilation dynamics all differed with respect to land use. We also show that microbial community dynamics affect carbon assimilation dynamics and are associated with soil DNA content. Soil DNA yield is easy to measure and may be useful in predicting microbial community dynamics linked to soil carbon cycling. Soil dwelling microorganisms are key players in the terrestrial carbon cycle, driving both the degradation and stabilization of soil organic matter. Microbial communities vary with respect to land use, but we still have an incomplete understanding of how variation in community structure links to variation in community function. DNA stable isotope probing (DNA-SIP) is a high-resolution method that can identify specific microbial taxa that assimilate carbon in situ. We conducted a large-scale multi-substrate DNA-SIP experiment to explore differences in bacterial activity across land-use regimes. We show that microbial community dynamics vary with land use, that these dynamics are linked to soil carbon cycling, and that they are associated with easily measured soil properties.

Watershed-scale liming reveals the short- and long-term effects of pH on the forest soil microbiome and carbon cycling.

Soil microbial community composition routinely correlates with pH, reflecting both direct pH effects on microbial physiology and long-term biogeochemical feedbacks. We used two watershed-scale liming experiments to identify short- (2 years) and long-term (25 years) changes in the structure and function of bacterial and fungal communities in organic horizons (Oe and Oa ) of acid forest soils. Liming increased soil pH, extractable calcium, and soil carbon stocks, reduced biomass-specific respiration, and caused major changes in the soil microbiome in the short and long term. More taxa responded to liming in the short term (70%) than in the long term (30%), with most showing consistent directional responses at both sites. The ratio of change in relative abundance between limed and reference sites was twofold higher at the long than the short-term site, indicating that the effects of liming grew over time. Liming impacts were most pronounced in fungi, as steep declines of dominant ectomycorrhizal fungi (Cenococcum and Russula) occurred at both sites. Liming favoured neutrophilic bacteria over acidophilic populations according to estimated environmental pH optima. Collectively, these results demonstrate that a liming-induced change of one pH unit has an immediate and persistent effect on the structure and function of microbial communities in acid forest soils. The corresponding suppression of respiration indicates that anthropogenic alterations of soil pH, as driven by acid deposition or liming, can affect forest floor C stocks due to pH-driven shifts in community structure.

Tracing Carbon Metabolism with Stable Isotope Metabolomics Reveals the Legacy of Diverse Carbon Sources in Soil.

Tracking the metabolic activity of whole soil communities can improve our understanding of the transformation and fate of carbon in soils. We used stable isotope metabolomics to trace 13C from nine labeled carbon sources into the water-soluble metabolite pool of an agricultural soil over time. Soil was amended with a mixture of all nine sources, with one source isotopically labeled in each treatment. We compared changes in the 13C enrichment of metabolites with respect to carbon source and time over a 48-day incubation and contrasted differences between soluble sources (glucose, xylose, amino acids, etc.) and insoluble sources (cellulose and palmitic acid). Whole soil metabolite profiles varied singularly by time, while the composition of 13C-labeled metabolites differed primarily by carbon source (R2 = 0.68) rather than time (R2 = 0.07), with source-specific differences persisting throughout incubations. The 13C labeling of metabolites from insoluble carbon sources occurred slower than that from soluble sources but yielded a higher average atom percent (atom%) 13C in metabolite markers of biomass (amino acids and nucleic acids). The 13C enrichment of metabolite markers of biomass stabilized between 5 and 15 atom% 13C by the end of incubations. Temporal patterns in the 13C enrichment of tricarboxylic acid cycle intermediates, nucleobases (uracil and thymine), and by-products of DNA salvage (allantoin) closely tracked microbial activity. Our results demonstrate that metabolite production in soils is driven by the carbon source supplied to the community and that the fate of carbon in metabolites do not generally converge over time as a result of ongoing microbial processing and recycling. IMPORTANCE Carbon metabolism in soil remains poorly described due to the inherent difficulty of obtaining information on the microbial metabolites produced by complex soil communities. Our study demonstrates the use of stable isotope probing (SIP) to study carbon metabolism in soil by tracking 13C from supplied carbon sources into metabolite pools and biomass. We show that differences in the metabolism of sources influence the fate of carbon in soils. Heterogeneity in 13C-labeled metabolite profiles corresponded with compositional differences in the metabolically active populations, providing a basis for how microbial community composition correlates with the quality of soil carbon. Our study demonstrates the application of SIP-metabolomics in studying soils and identifies several metabolite markers of growth, activity, and other aspects of microbial function.

Microbial community shifts correspond with suppression of decomposition 25 years after liming of acidic forest soils.

Microbial community structure and function regularly covary with soil pH, yet effects of these interactions on soil carbon are rarely tested experimentally within natural ecosystems. We investigated the enduring (25 year) impacts of liming on microbial community structure and decomposition at an acidic northern hardwood forest, where experimental liming increased pH one unit and surprisingly doubled the organic carbon stocks of the forest floor. We show that this increase in carbon storage corresponded with restructuring of the bacterial and fungal communities that drive decomposition. In the Oe horizon, liming reduced the activities of five extracellular enzymes that mediate decomposition, while the Oa horizon showed an especially large (64%) reduction in the activity of a sixth, peroxidase, which is an oxidative enzyme central to lignocellulose degradation. Decreased enzyme activities corresponded with loss of microbial taxa important for lignocellulose decay, including large reductions in the dominant ectomycorrhizal genera Russula and Cenococcum, saprotrophic and wood decaying fungi, and Actinobacteria (Thermomonosporaceae). These results demonstrate the importance of pH as a dominant regulator of microbial community structure and illustrate how changes to this structure can produce large, otherwise unexpected increases in carbon storage in forest soils.

Streptomyces apricus sp. nov., isolated from soil.

A novel Streptomyces strain, SUN51T, was isolated from soils sampled in Wisconsin, USA, as part of a Streptomyces biogeography survey. Genome sequencing revealed that this strain had less than 90 % average nucleotide identity (ANI) to type species of Streptomyces: SUN51T was most closely related to Streptomyces dioscori A217T (99.5 % 16S rRNA gene identity, 89.4 % ANI). Genome size was estimated at 8.81 Mb, and the genome DNA G+C content was 72 mol%. The strain possessed the cellular fatty acids anteiso-C15 : 0, iso-C16 : 0, 16 : 1 ω7c, anteiso-C17 : 0, iso-C14 : 0 and C16 : 0. The predominant menaquinones were MK-9 H4, MK-9 H6 and MK-9 H8. Strain SUN51T contained the polar lipids phosphatidic acid, phosphatidyl ethanolamine, phosphatidyl glycerol and diphosphatidyl glycerol. The cell wall contained ll-diaminopimelic acid. The strain could grow on a broad range of carbon sources and tolerate temperatures of up to 40 °C. The results of the polyphasic study confirmed that this isolate represents a novel species of the genus Streptomyces, for which the name Streptomyces apricus sp. nov. is proposed. The type strain of this species is SUN51T (=NRRL B-65543T=JCM 33736T).

The biogeography of Streptomyces in New Zealand enabled by high-throughput sequencing of genus-specific rpoB amplicons.

We evaluated Streptomyces biogeography in soils along a 1200 km latitudinal transect across New Zealand (NZ). Streptomyces diversity was examined using high-throughput sequencing of rpoB amplicons generated with a Streptomyces specific primer set. We detected 1287 Streptomyces rpoB operational taxonomic units (OTUs) with 159 ± 92 (average ± SD) rpoB OTUs per site. Only 12% (n = 149) of these OTUs matched rpoB sequences from cultured specimens (99% nucleotide identity cutoff). Streptomyces phylogenetic diversity (Faith's PD) was correlated with soil pH, mean annual temperature and plant community richness (Spearman's r: 0.77, 0.64 and -0.79, respectively; P < 0.05), but not with latitude. In addition, soil pH and plant community richness both explained significant variation in Streptomyces beta diversity. Streptomyces communities exhibited both high dissimilarity and strong dominance of one or a few species at each site. Taken together, these results suggest that dispersal limitation due to competitive interactions limits the colonization success of spores that relocate to new sites. Cultivated Streptomyces isolates represent a major source of clinically useful antibiotics, but only a small fraction of extant diversity within the genus have been identified and most species of Streptomyces have yet to be described.

Simulating metagenomic stable isotope probing datasets with MetaSIPSim.

BACKGROUND: DNA-stable isotope probing (DNA-SIP) links microorganisms to their in-situ function in diverse environmental samples. Combining DNA-SIP and metagenomics (metagenomic-SIP) allows us to link genomes from complex communities to their specific functions and improves the assembly and binning of these targeted genomes. However, empirical development of metagenomic-SIP methods is hindered by the complexity and cost of these studies. We developed a toolkit, 'MetaSIPSim,' to simulate sequencing read libraries for metagenomic-SIP experiments. MetaSIPSim is intended to generate datasets for method development and testing. To this end, we used MetaSIPSim generated data to demonstrate the advantages of metagenomic-SIP over a conventional shotgun metagenomic sequencing experiment. RESULTS: Through simulation we show that metagenomic-SIP improves the assembly and binning of isotopically labeled genomes relative to a conventional metagenomic approach. Improvements were dependent on experimental parameters and on sequencing depth. Community level G + C content impacted the assembly of labeled genomes and subsequent binning, where high community G + C generally reduced the benefits of metagenomic-SIP. Furthermore, when a high proportion of the community is isotopically labeled, the benefits of metagenomic-SIP decline. Finally, the choice of gradient fractions to sequence greatly influences method performance. CONCLUSIONS: Metagenomic-SIP is a valuable method for recovering isotopically labeled genomes from complex communities. We show that metagenomic-SIP performance depends on optimization of experimental parameters. MetaSIPSim allows for simulation of metagenomic-SIP datasets which facilitates the optimization and development of metagenomic-SIP experiments and analytical approaches for dealing with these data.

Plant growth rate and nitrogen uptake shape rhizosphere bacterial community composition and activity in an agricultural field.

Plant-microbial interactions in the rhizosphere are an essential link in soil nitrogen (N) cycling and plant N supply. Plant phenotype and genotype interact with the soil environment to determine rhizosphere community structure and activity. However, the relative contributions of plant identity, phenology and soil resource availability in shaping rhizosphere effects are not well understood. Four summer annuals and a collection of maize hybrids were grown in a common garden experiment conducted at two levels of organic nutrient availability. Plant biomass, N accumulation, rhizosphere bacterial community composition, and rhizosphere potential extracellular enzyme activity were assessed at vegetative, flowering and grain-filling stages of maize. Plant N uptake was strongly coupled with protease activity in the rhizosphere. Temporal trends in rhizosphere community composition varied between plant species. Changes in rhizosphere community composition could be explained by variation in plant growth dynamics. These findings indicate that species-level variation in plant growth dynamics and resource acquisition drive variation in rhizosphere bacterial community composition and activity linked to plant N uptake.

Paraburkholderia madseniana sp. nov., a phenolic acid-degrading bacterium isolated from acidic forest soil.

RP11T was isolated from forest soil following enrichment with 4-hydroxybenzoic acid. Cells of RP11T are aerobic, non-sporulating, exhibit swimming motility, and are rods (0.8 µm by 1.4 µm) that often occur as diplobacillus or in short chains (3-4 cells). Optimal growth on minimal media containing 4-hydroxybenzoic acid (µ=0.216 hr-1) occurred at 30 °C, pH 6.5 or 7.0 and 0% salinity. Comparative chemotaxonomic, genomic and phylogenetic analyses revealed the isolate was distinct from its closest relative type strains identified as Paraburkholderia aspalathi LMG 27731T, Paraburkholderia fungorum LMG 16225T and Paraburkholderia caffeinilytica CF1T. Strain RP11T is genetically distinct from P. aspalathi, its closest relative, in terms of 16S rRNA gene sequence similarity (98.7%), genomic average nucleotide identity (94%) and in silico DNA-DNA hybridization (56.7 %±2.8). The composition of fatty acids and substrate utilization pattern differentiated strain RP11T from its closest relatives, including growth on phthalic acid. Strain RP11T encoded the greatest number of aromatic degradation genes of all eleven closely related type strains and uniquely encoded a phthalic acid dioxygenase and paralog of the 3-hydroxybenzoate 4-monooxygenase. The only ubiquinone detected in strain RP11T was Q-8, and the major cellular fatty acids were C16 : 0, 3OH-C16 : 0, C17 : 0 cyclo, C19 : 0 cyclo ω8c, and summed feature 8 (C18 : 1 ω7c/ω6c). On the basis of this polyphasic approach, it was determined that strain RP11T represents a novel species from the genus Paraburkholderia for which the name Paraburkholderia madseniana sp. nov. is proposed. The type strain is RP11T (=DSM 110123T=LMG 31517T).

Paraburkholderia solitsugae sp. nov. and Paraburkholderia elongata sp. nov., phenolic acid-degrading bacteria isolated from forest soil and emended description of Paraburkholderia madseniana.

Two bacterial strains, 1NT and 5NT, were isolated from hemlock forest soil using a soluble organic matter enrichment. Cells of 1NT (0.65×1.85 µm) and 5NT (0.6×1.85 µm) are Gram-stain-negative, aerobic, motile, non-sporulating and exist as single rods, diplobacilli or in chains of varying length. During growth in dilute media (≤0.1× tryptic soy broth; TSB), cells are primarily motile with flagella. At higher concentrations (≥0.3× TSB), cells of both strains increasingly form non-motile chains, and cells of 5NT elongate (0.57×~7 µm) and form especially long filaments. Optimum growth of 1NT and 5NT occurred at 25-30 °C, pH 6.5-7.0 and <0.5% salinity. Results of comparative chemotaxonomic, genomic and phylogenetic analyses revealed that 1NT and 5NT were distinct from one another and their closest related type strains: Paraburkholderia madseniana RP11T, Paraburkholderia aspalathi LMG 27731T and Paraburkholderia caffeinilytica CF1T. The genomes of 1NT and 5NT had an average nucleotide identity (91.6 and 91.3%) and in silico DNA-DNA hybridization values (45.8%±2.6 and 45.5%±2.5) and differed in functional gene content from their closest related type strains. The composition of fatty acids and patterns of substrate use, including the catabolism of phenolic acids, also differentiated strains 1NT and 5NT from each other and their closest relatives. The only ubiquinone present in strains 1NT and 5NT was Q-8. The major cellular fatty acids were C16 : 0, 3OH-C16 : 0, C17 : 0 cyclo, C19 : 0 cyclo ω8c and summed features 2 (3OH-C14 : 0 / C16 : 1 iso I), 3 (C16 : 1 ω6c/ω7c) and 8 (C18 : 1 ω7c/ω6c). A third bacterium, strain RL16-012-BIC-B, was isolated from soil associated with shallow roots and was determined to be a strain of P. madseniana (ANI, 98.8%; 16S rRNA gene similarity, 100%). Characterizations of strain RL16-012-BIC-B (DSM 110723=LMG 31706) led to proposed emendments to the species description of P. madseniana. Our polyphasic approach demonstrated that strains 1NT and 5NT represent novel species from the genus Paraburkholderia for which the names Paraburkholderia solitsugae sp. nov. (type strain 1NT=DSM 110721T=LMG 31704T) and Paraburkholderia elongata sp. nov. (type strain 5NT=DSM 110722T=LMG 31705T) are proposed.

The Use of Degenerate Primers in qPCR Analysis of Functional Genes Can Cause Dramatic Quantification Bias as Revealed by Investigation of nifH Primer Performance.

The measurement of functional gene abundance in diverse microbial communities often employs quantitative PCR (qPCR) with highly degenerate oligonucleotide primers. While degenerate PCR primers have been demonstrated to cause template-specific bias in PCR applications, the effect of such bias on qPCR has been less well explored. We used a set of diverse, full-length nifH gene standards to test the performance of several universal nifH primer sets in qPCR. We found significant template-specific bias in all but the PolF/PolR primer set. Template-specific bias caused more than 1000-fold mis-estimation of nifH gene copy number for three of the primer sets and one primer set resulted in more than 10,000-fold mis-estimation. Furthermore, such template-specific bias will cause qPCR estimates to vary in response to beta-diversity, thereby causing mis-estimation of changes in gene copy number. A reduction in bias was achieved by increasing the primer concentration. We conclude that degenerate primers should be evaluated across a range of templates, annealing temperatures, and primer concentrations to evaluate the potential for template-specific bias prior to their use in qPCR.

Latitude delineates patterns of biogeography in terrestrial Streptomyces.

The biogeography of Streptomyces was examined at regional spatial scales to identify factors that govern patterns of microbial diversity. Streptomyces are spore forming filamentous bacteria which are widespread in soil. Streptomyces strains were isolated from perennial grass habitats sampled across a spatial scale of more than 6000 km. Previous analysis of this geographically explicit culture collection provided evidence for a latitudinal diversity gradient in Streptomyces species. Here the hypothesis that this latitudinal diversity gradient is a result of evolutionary dynamics associated with historical demographic processes was evaluated. Historical demographic phenomena have genetic consequences that can be evaluated through analysis of population genetics. Population genetic approaches were applied to analyze population structure in six of the most numerically abundant and geographically widespread Streptomyces phylogroups from our culture collection. Streptomyces population structure varied at regional spatial scales, and allelic diversity correlated with geographic distance. In addition, allelic diversity and gene flow are partitioned by latitude. Finally, it was found that nucleotide diversity within phylogroups was negatively correlated with latitude. These results indicate that phylogroup diversification is constrained by dispersal limitation at regional spatial scales, and they are consistent with the hypothesis that historical demographic processes have influenced the contemporary biogeography of Streptomyces.

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin.

BACKGROUND: Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized. RESULTS: We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system. CONCLUSIONS: S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.

Microscale sulfur cycling in the phototrophic pink berry consortia of the Sippewissett Salt Marsh.

Microbial metabolism is the engine that drives global biogeochemical cycles, yet many key transformations are carried out by microbial consortia over short spatiotemporal scales that elude detection by traditional analytical approaches. We investigate syntrophic sulfur cycling in the 'pink berry' consortia of the Sippewissett Salt Marsh through an integrative study at the microbial scale. The pink berries are macroscopic, photosynthetic microbial aggregates composed primarily of two closely associated species: sulfide-oxidizing purple sulfur bacteria (PB-PSB1) and sulfate-reducing bacteria (PB-SRB1). Using metagenomic sequencing and (34) S-enriched sulfate stable isotope probing coupled with nanoSIMS, we demonstrate interspecies transfer of reduced sulfur metabolites from PB-SRB1 to PB-PSB1. The pink berries catalyse net sulfide oxidation and maintain internal sulfide concentrations of 0-500 µm. Sulfide within the berries, captured on silver wires and analysed using secondary ion mass spectrometer, increased in abundance towards the berry interior, while δ(34) S-sulfide decreased from 6‰ to -31‰ from the exterior to interior of the berry. These values correspond to sulfate-sulfide isotopic fractionations (15-53‰) consistent with either sulfate reduction or a mixture of reductive and oxidative metabolisms. Together this combined metagenomic and high-resolution isotopic analysis demonstrates active sulfur cycling at the microscale within well-structured macroscopic consortia consisting of sulfide-oxidizing anoxygenic phototrophs and sulfate-reducing bacteria.