Horizontal gene transfer of the Mer operon is associated with large effects on the transcriptome and increased tolerance to mercury in nitrogen-fixing bacteria.

BACKGROUND: Mercury (Hg) is highly toxic and has the potential to cause severe health problems for humans and foraging animals when transported into edible plant parts. Soil rhizobia that form symbiosis with legumes may possess mechanisms to prevent heavy metal translocation from roots to shoots in plants by exporting metals from nodules or compartmentalizing metal ions inside nodules. Horizontal gene transfer has potential to confer immediate de novo adaptations to stress. We used comparative genomics of high quality de novo assemblies to identify structural differences in the genomes of nitrogen-fixing rhizobia that were isolated from a mercury (Hg) mine site that show high variation in their tolerance to Hg. RESULTS: Our analyses identified multiple structurally conserved merA homologs in the genomes of Sinorhizobium medicae and Rhizobium leguminosarum but only the strains that possessed a Mer operon exhibited 10-fold increased tolerance to Hg. RNAseq analysis revealed nearly all genes in the Mer operon were significantly up-regulated in response to Hg stress in free-living conditions and in nodules. In both free-living and nodule environments, we found the Hg-tolerant strains with a Mer operon exhibited the fewest number of differentially expressed genes (DEGs) in the genome, indicating a rapid and efficient detoxification of Hg from the cells that reduced general stress responses to the Hg-treatment. Expression changes in S. medicae while in bacteroids showed that both rhizobia strain and host-plant tolerance affected the number of DEGs. Aside from Mer operon genes, nif genes which are involved in nitrogenase activity in S. medicae showed significant up-regulation in the most Hg-tolerant strain while inside the most Hg-accumulating host-plant. Transfer of a plasmid containing the Mer operon from the most tolerant strain to low-tolerant strains resulted in an immediate increase in Hg tolerance, indicating that the Mer operon is able to confer hyper tolerance to Hg. CONCLUSIONS: Mer operons have not been previously reported in nitrogen-fixing rhizobia. This study demonstrates a pivotal role of the Mer operon in effective mercury detoxification and hypertolerance in nitrogen-fixing rhizobia. This finding has major implications not only for soil bioremediation, but also host plants growing in mercury contaminated soils.

Ecological and genomic analyses of candidate phylum WPS-2 bacteria in an unvegetated soil.

Members of the bacterial candidate phylum WPS-2 (or Eremiobacterota) are abundant in several dry, bare soil environments. In a bare soil deposited by an extinct iron-sulfur spring, we found that WPS-2 comprised up to 24% of the bacterial community and up to 108 cells per g of soil based on 16S rRNA gene sequencing and quantification. A single genus-level cluster (Ca. Rubrimentiphilum) predominated in bare soils but was less abundant in adjacent forest. Nearly complete genomes of Ca. Rubrimentiphilum were recovered as single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs). Surprisingly, given the abundance of WPS-2 in bare soils, the genomes did not indicate any capacity for autotrophy, phototrophy, or trace gas metabolism. Instead, they suggest a predominantly aerobic organoheterotrophic lifestyle, perhaps based on scavenging amino acids, nucleotides, and complex oligopeptides, along with lithotrophic capacity on thiosulfate. Network analyses of the entire community showed that some species of Chloroflexi, Actinobacteria, and candidate phylum AD3 (or Dormibacterota) co-occurred with Ca. Rubrimentiphilum and may represent ecological or metabolic partners. We propose that Ca. Rubrimentiphilum act as efficient heterotrophic scavengers. Combined with previous studies, these data suggest that the phylum WPS-2 includes bacteria with diverse metabolic capabilities.

Impact of library preparation protocols and template quantity on the metagenomic reconstruction of a mock microbial community.

BACKGROUND: The rapid development of sequencing technologies has provided access to environments that were either once thought inhospitable to life altogether or that contain too few cells to be analyzed using genomics approaches. While 16S rRNA gene microbial community sequencing has revolutionized our understanding of community composition and diversity over time and space, it only provides a crude estimate of microbial functional and metabolic potential. Alternatively, shotgun metagenomics allows comprehensive sampling of all genetic material in an environment, without any underlying primer biases. Until recently, one of the major bottlenecks of shotgun metagenomics has been the requirement for large initial DNA template quantities during library preparation. RESULTS: Here, we investigate the effects of varying template concentrations across three low biomass library preparation protocols on their ability to accurately reconstruct a mock microbial community of known composition. We analyze the effects of input DNA quantity and library preparation method on library insert size, GC content, community composition, assembly quality and metagenomic binning. We found that library preparation method and the amount of starting material had significant impacts on the mock community metagenomes. In particular, GC content shifted towards more GC rich sequences at the lower input quantities regardless of library prep method, the number of low quality reads that could not be mapped to the reference genomes increased with decreasing input quantities, and the different library preparation methods had an impact on overall metagenomic community composition. CONCLUSIONS: This benchmark study provides recommendations for library creation of representative and minimally biased metagenome shotgun sequencing, enabling insights into functional attributes of low biomass ecosystem microbial communities.

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake.

BACKGROUND: Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria. RESULTS: Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota. CONCLUSIONS: We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

scMicrobe PTA: near complete genomes from single bacterial cells.

Microbial genomes produced by standard single-cell amplification methods are largely incomplete. Here, we show that primary template-directed amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard multiple displacement amplification-based approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

scMicrobe PTA: Near Complete Genomes from Single Bacterial Cells.

Microbial genomes produced by single-cell amplification are largely incomplete. Here, we show that primary template amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard amplification approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

Genetic and Functional Diversity Help Explain Pathogenic, Weakly Pathogenic, and Commensal Lifestyles in the Genus Xanthomonas.

The genus Xanthomonas has been primarily studied for pathogenic interactions with plants. However, besides host and tissue-specific pathogenic strains, this genus also comprises nonpathogenic strains isolated from a broad range of hosts, sometimes in association with pathogenic strains, and other environments, including rainwater. Based on their incapacity or limited capacity to cause symptoms on the host of isolation, nonpathogenic xanthomonads can be further characterized as commensal and weakly pathogenic. This study aimed to understand the diversity and evolution of nonpathogenic xanthomonads compared to their pathogenic counterparts based on their cooccurrence and phylogenetic relationship and to identify genomic traits that form the basis of a life history framework that groups xanthomonads by ecological strategies. We sequenced genomes of 83 strains spanning the genus phylogeny and identified eight novel species, indicating unexplored diversity. While some nonpathogenic species have experienced a recent loss of a type III secretion system, specifically the hrp2 cluster, we observed an apparent lack of association of the hrp2 cluster with lifestyles of diverse species. We performed association analysis on a large data set of 337 Xanthomonas strains to explain how xanthomonads may have established association with the plants across the continuum of lifestyles from commensals to weak pathogens to pathogens. Presence of distinct transcriptional regulators, distinct nutrient utilization and assimilation genes, transcriptional regulators, and chemotaxis genes may explain lifestyle-specific adaptations of xanthomonads.

Seasonal variability in bacterial and fungal diversity of the near-surface atmosphere.

Bacteria and fungi are ubiquitous throughout the Earth's lower atmosphere where they often represent an important component of atmospheric aerosols with the potential to impact human health and atmospheric dynamics. However, the diversity, composition, and spatiotemporal dynamics of these airborne microbes remain poorly understood. We performed a comprehensive analysis of airborne microbes across two aerosol size fractions at urban and rural sites in the Colorado Front Range over a 14-month period. Coarse (PM10-2.5) and fine (PM2.5) particulate matter samples were collected at weekly intervals with both bacterial and fungal diversity assessed via high-throughput sequencing. The diversity and composition of the airborne communities varied across the sites, between the two size fractions, and over time. Bacteria were the dominant type of bioaerosol in the collected air samples, while fungi and plants (pollen) made up the remainder, with the relative abundances of fungi peaking during the spring and summer months. As bacteria made up the majority of bioaerosol particles, we analyzed the bacterial communities in greater detail using a bacterial-specific 16S rRNA gene sequencing approach. Overall, bacterial taxonomic richness and the relative abundances of specific bacterial taxa exhibited significant patterns of seasonality. Likewise, airborne bacterial communities varied significantly between sites and across aerosol size fractions. Source-tracking analyses indicate that soils and leaves represented important sources of bacteria to the near-surface atmosphere across all locations with cow fecal bacteria also representing an important source of bioaerosols at the more rural sites during early fall and early spring. Together, these data suggest that a complex set of environmental factors, including changes in atmospheric conditions and shifts in the relative importance of available microbial sources, act to control the composition of microbial bioaerosols in rural and urban environments.

Terabase-Scale Coassembly of a Tropical Soil Microbiome.

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.

Synthase-selected sorting approach identifies a beta-lactone synthase in a nudibranch symbiotic bacterium.

BACKGROUND: Nudibranchs comprise a group of > 6000 marine soft-bodied mollusk species known to use secondary metabolites (natural products) for chemical defense. The full diversity of these metabolites and whether symbiotic microbes are responsible for their synthesis remains unexplored. Another issue in searching for undiscovered natural products is that computational analysis of genomes of uncultured microbes can result in detection of novel biosynthetic gene clusters; however, their in vivo functionality is not guaranteed which limits further exploration of their pharmaceutical or industrial potential. To overcome these challenges, we used a fluorescent pantetheine probe, which produces a fluorescent CoA-analog employed in biosynthesis of secondary metabolites, to label and capture bacterial symbionts actively producing these compounds in the mantle of the nudibranch Doriopsilla fulva. RESULTS: We recovered the genome of Candidatus Doriopsillibacter californiensis from the Ca. Tethybacterales order, an uncultured lineage of sponge symbionts not found in nudibranchs previously. It forms part of the core skin microbiome of D. fulva and is nearly absent in its internal organs. We showed that crude extracts of D. fulva contained secondary metabolites that were consistent with the presence of a beta-lactone encoded in Ca. D. californiensis genome. Beta-lactones represent an underexplored group of secondary metabolites with pharmaceutical potential that have not been reported in nudibranchs previously. CONCLUSIONS: Altogether, this study shows how probe-based, targeted sorting approaches can capture bacterial symbionts producing secondary metabolites in vivo. Video Abstract.

Metagenomes and Metagenome-Assembled Genomes from Substrate-Amended Hot Spring Sediment Incubations from Yellowstone National Park.

Here, we report on eight sediment metagenomes obtained from an alkaline hot spring, with their corresponding metagenome-assembled genomes. Samples had been incubated for 48 h with various substrate amendments in conjunction with the amino acid analog l-homopropargylglycine in a study targeted at identifying anabolicly active uncultured thermophilic archaea and bacteria.

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution.

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Sources of bacteria in outdoor air across cities in the midwestern United States.

Bacteria are abundant in the atmosphere, where they often represent a major portion of the organic aerosols. Potential pathogens of plants and livestock are commonly dispersed through the atmosphere, and airborne bacteria can have important effects on human health as pathogens or triggers of allergic asthma and seasonal allergies. Despite their importance, the diversity and biogeography of airborne microorganisms remain poorly understood. We used high-throughput pyrosequencing to analyze bacterial communities present in the aerosol fraction containing fine particulate matter of ≤2.5 µm from 96 near-surface atmospheric samples collected from cities throughout the midwestern United States and found that the communities are surprisingly diverse and strongly affected by the season. We also directly compared the airborne communities to those found in hundreds of samples representing potential source environments. We show that, in addition to the more predictable sources (soils and leaf surfaces), fecal material, most likely dog feces, often represents an unexpected source of bacteria in the atmosphere at more urbanized locations during the winter. Airborne bacteria are clearly an important, but understudied, component of air quality that needs to be better integrated into efforts to measure and model pollutants in the atmosphere.

Effect of template on generating a standard curve for absolute quantification of an RNA virus by real-time reverse transcriptase-polymerase chain reaction.

The effect of different templates on generating standard curves that are needed for the absolute quantification of an RNA virus by real-time reverse transcriptase-PCR (RT-PCR) was evaluated. We used infectious pancreatic necrosis virus (IPNV), a major viral pathogen of wild and cultured salmon, as a RNA virus example for the study. A dilution series of four different templates representing the IPNV protease gene (two in vitro transcribed RNAs of 100 bases and 500 bases in length, a plasmid DNA and a DNA oligo) were used as template to produce standard curves to quantify IPNV load in rainbow trout. The slope, the goodness of fit (r(2)), and the efficiency (e) of PCR were statistically equivalent irrespective of the nature of template used in the PCR. Using a factorial ANOVA, no significant difference in IPNV copy number was observed using the four different standard curves for absolute quantification of IPNV in experimentally-challenged rainbow trout. However, when IPNV transcript abundance was less than 100 copies per reaction and when the template size was bigger than the amplicon size amplification was more variable. The data suggests that the size of the template used to generate standard curve should be very similar to the size of the amplicon. A synthetic DNA oligo template would be optimal for this purpose as it can be custom made and only requires the sequence information for its synthesis. However, if the standard curve is generated with template copy number in excess of 100 copies per reaction, the nature of the template has no effect on the standard curve, and therefore, the cheaper template would be the preferred choice of template over the other more expensive options.

High Potential for Biomass-Degrading Enzymes Revealed by Hot Spring Metagenomics.

Enzyme stability and activity at elevated temperatures are important aspects in biotechnological industries, such as the conversion of plant biomass into biofuels. In order to reduce the costs and increase the efficiency of biomass conversion, better enzymatic processing must be developed. Hot springs represent a treasure trove of underexplored microbiological and protein chemistry diversity. Herein, we conduct an exploratory study into the diversity of hot spring biomass-degrading potential. We describe the taxonomic diversity and carbohydrate active enzyme (CAZyme) coding potential in 71 publicly available metagenomic datasets from 58 globally distributed terrestrial geothermal features. Through taxonomic profiling, we detected a wide diversity of microbes unique to varying temperature and pH ranges. Biomass-degrading enzyme potential included all five classes of CAZymes and we described the presence or absence of genes encoding 19 glycosyl hydrolases hypothesized to be involved with cellulose, hemicellulose, and oligosaccharide degradation. Our results highlight hot springs as a promising system for the further discovery and development of thermo-stable biomass-degrading enzymes that can be applied toward generation of renewable biofuels. This study lays a foundation for future research to further investigate the functional diversity of hot spring biomass-degrading enzymes and their potential utility in biotechnological processing.

Zymomonas diversity and potential for biofuel production.

BACKGROUND: Zymomonas mobilis is an aerotolerant α-proteobacterium, which has been genetically engineered for industrial purposes for decades. However, a comprehensive comparison of existing strains on the genomic level in conjunction with phenotype analysis has yet to be carried out. We here performed whole-genome comparison of 17 strains including nine that were sequenced in this study. We then compared 15 available Zymomonas strains for their natural abilities to perform under conditions relevant to biofuel synthesis. We tested their growth in anaerobic rich media, as well as growth, ethanol production and xylose utilization in lignocellulosic hydrolysate. We additionally compared their tolerance to isobutanol, flocculation characteristics, and ability to uptake foreign DNA by electroporation and conjugation. RESULTS: Using clustering based on 99% average nucleotide identity (ANI), we classified 12 strains into four clusters based on sequence similarity, while five strains did not cluster with any other strain. Strains belonging to the same 99% ANI cluster showed similar performance while significant variation was observed between the clusters. Overall, conjugation and electroporation efficiencies were poor across all strains, which was consistent with our finding of coding potential for several DNA defense mechanisms, such as CRISPR and restriction-modification systems, across all genomes. We found that strain ATCC31821 (ZM4) had a more diverse plasmid profile than other strains, possibly leading to the unique phenotypes observed for this strain. ZM4 also showed the highest growth of any strain in both laboratory media and lignocellulosic hydrolysate and was among the top 3 strains for isobutanol tolerance and electroporation and conjugation efficiency. CONCLUSIONS: Our findings suggest that strain ZM4 has a unique combination of genetic and phenotypic traits that are beneficial for biofuel production and propose investing future efforts in further engineering of ZM4 for industrial purposes rather than exploring new Zymomonas isolates.

GAL08, an Uncultivated Group of Acidobacteria, Is a Dominant Bacterial Clade in a Neutral Hot Spring.

GAL08 are bacteria belonging to an uncultivated phylogenetic cluster within the phylum Acidobacteria. We detected a natural population of the GAL08 clade in sediment from a pH-neutral hot spring located in British Columbia, Canada. To shed light on the abundance and genomic potential of this clade, we collected and analyzed hot spring sediment samples over a temperature range of 24.2-79.8°C. Illumina sequencing of 16S rRNA gene amplicons and qPCR using a primer set developed specifically to detect the GAL08 16S rRNA gene revealed that absolute and relative abundances of GAL08 peaked at 65°C along three temperature gradients. Analysis of sediment collected over multiple years and locations revealed that the GAL08 group was consistently a dominant clade, comprising up to 29.2% of the microbial community based on relative read abundance and up to 4.7 × 105 16S rRNA gene copy numbers per gram of sediment based on qPCR. Using a medium quality threshold, 25 single amplified genomes (SAGs) representing these bacteria were generated from samples taken at 65 and 77°C, and seven metagenome-assembled genomes (MAGs) were reconstructed from samples collected at 45-77°C. Based on average nucleotide identity (ANI), these SAGs and MAGs represented three separate species, with an estimated average genome size of 3.17 Mb and GC content of 62.8%. Phylogenetic trees constructed from 16S rRNA gene sequences and a set of 56 concatenated phylogenetic marker genes both placed the three GAL08 bacteria as a distinct subgroup of the phylum Acidobacteria, representing a candidate order (Ca. Frugalibacteriales) within the class Blastocatellia. Metabolic reconstructions from genome data predicted a heterotrophic metabolism, with potential capability for aerobic respiration, as well as incomplete denitrification and fermentation. In laboratory cultivation efforts, GAL08 counts based on qPCR declined rapidly under atmospheric levels of oxygen but increased slightly at 1% (v/v) O2, suggesting a microaerophilic lifestyle.

The ecology of the phyllosphere: geographic and phylogenetic variability in the distribution of bacteria on tree leaves.

Large populations of bacteria live on leaf surfaces and these phyllosphere bacteria can have important effects on plant health. However, we currently have a limited understanding of bacterial diversity on tree leaves and the inter- and intra-specific variability in phyllosphere community structure. We used a barcoded pyrosequencing technique to characterize the bacterial communities from leaves of 56 tree species in Boulder, Colorado, USA, quantifying the intra- and inter-individual variability in the bacterial communities from 10 of these species. We also examined the geographic variability in phyllosphere communities on Pinus ponderosa from several locations across the globe. Individual tree species harboured high levels of bacterial diversity and there was considerable variability in community composition between trees. The bacterial communities were organized in patterns predictable from the relatedness of the trees as there was significant correspondence between tree phylogeny and bacterial community phylogeny. Inter-specific variability in bacterial community composition exceeded intra-specific variability, a pattern that held even across continents where we observed minimal geographic differentiation in the bacterial communities on P. ponderosa needles.

A pipeline for targeted metagenomics of environmental bacteria.

BACKGROUND: Metagenomics and single cell genomics provide a window into the genetic repertoire of yet uncultivated microorganisms, but both methods are usually taxonomically untargeted. The combination of fluorescence in situ hybridization (FISH) and fluorescence activated cell sorting (FACS) has the potential to enrich taxonomically well-defined clades for genomic analyses. METHODS: Cells hybridized with a taxon-specific FISH probe are enriched based on their fluorescence signal via flow cytometric cell sorting. A recently developed FISH procedure, the hybridization chain reaction (HCR)-FISH, provides the high signal intensities required for flow cytometric sorting while maintaining the integrity of the cellular DNA for subsequent genome sequencing. Sorted cells are subjected to shotgun sequencing, resulting in targeted metagenomes of low diversity. RESULTS: Pure cultures of different taxonomic groups were used to (1) adapt and optimize the HCR-FISH protocol and (2) assess the effects of various cell fixation methods on both the signal intensity for cell sorting and the quality of subsequent genome amplification and sequencing. Best results were obtained for ethanol-fixed cells in terms of both HCR-FISH signal intensity and genome assembly quality. Our newly developed pipeline was successfully applied to a marine plankton sample from the North Sea yielding good quality metagenome assembled genomes from a yet uncultivated flavobacterial clade. CONCLUSIONS: With the developed pipeline, targeted metagenomes at various taxonomic levels can be efficiently retrieved from environmental samples. The resulting metagenome assembled genomes allow for the description of yet uncharacterized microbial clades. Video abstract.

Function-driven single-cell genomics uncovers cellulose-degrading bacteria from the rare biosphere.

Assigning a functional role to a microorganism has historically relied on cultivation of isolates or detection of environmental genome-based biomarkers using a posteriori knowledge of function. However, the emerging field of function-driven single-cell genomics aims to expand this paradigm by identifying and capturing individual microbes based on their in situ functions or traits. To identify and characterize yet uncultivated microbial taxa involved in cellulose degradation, we developed and benchmarked a function-driven single-cell screen, which we applied to a microbial community inhabiting the Great Boiling Spring (GBS) Geothermal Field, northwest Nevada. Our approach involved recruiting microbes to fluorescently labeled cellulose particles, and then isolating single microbe-bound particles via fluorescence-activated cell sorting. The microbial community profiles prior to sorting were determined via bulk sample 16S rRNA gene amplicon sequencing. The flow-sorted cellulose-bound microbes were subjected to whole genome amplification and shotgun sequencing, followed by phylogenetic placement. Next, putative cellulase genes were identified, expressed and tested for activity against derivatives of cellulose and xylose. Alongside typical cellulose degraders, including members of the Actinobacteria, Bacteroidetes, and Chloroflexi, we found divergent cellulases encoded in the genome of a recently described candidate phylum from the rare biosphere, Goldbacteria, and validated their cellulase activity. As this genome represents a species-level organism with novel and phylogenetically distinct cellulolytic activity, we propose the name Candidatus 'Cellulosimonas argentiregionis'. We expect that this function-driven single-cell approach can be extended to a broad range of substrates, linking microbial taxonomy directly to in situ function.