Population genomic analysis of strain variation in Leptospirillum group II bacteria involved in acid mine drainage formation

Modern microbiology: Embracing complexity through integration across scales

Microbes were the only form of life on Earth for most of its history, and they still account for the vast majority of life's diversity. They convert rocks to soil, produce much of the oxygen we breathe, remediate our sewage, and sustain agriculture. Microbes are vital to planetary health as they maintain biogeochemical cycles that produce and consume major greenhouse gases and support large food webs. Modern microbiologists analyze nucleic acids, proteins, and metabolites; leverage sophisticated genetic tools, software, and bioinformatic algorithms; and process and integrate complex and heterogeneous datasets so that microbial systems may be harnessed to address contemporary challenges in health, the environment, and basic science. Here, we consider an inevitably incomplete list of emergent themes in our discipline and highlight those that we recognize as the archetypes of its modern era that aim to address the most pressing problems of the 21st century.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Metapangenomic investigation provides insight into niche differentiation of methanogenic populations from the subsurface serpentinizing environment, Samail Ophiolite, Oman

Serpentinization reactions produce highly reduced waters that have hyperalkaline pH and that can have high concentrations of H₂ and CH₄. Putatively autotrophic methanogenic archaea have been identified in the subsurface waters of the Samail Ophiolite, Sultanate of Oman, though the strategies to overcome hyperalkaline pH and dissolved inorganic carbon limitation remain to be fully understood. Here, we recovered metagenome assembled genomes (MAGs) and applied a metapangenomic approach to three different Methanobacterium populations to assess habitat-specific functional gene distribution. A Type I population was identified in the fluids with neutral pH, while a Type II and "Mixed" population were identified in the most hyperalkaline fluids (pH 11.63). The core genome of all Methanobacterium populations highlighted potential DNA scavenging techniques to overcome phosphate or nitrogen limitation induced by environmental conditions. With particular emphasis on the Mixed and Type II population found in the most hyperalkaline fluids, the accessory genomes unique to each population reflected adaptation mechanisms suggesting lifestyles that minimize niche overlap. In addition to previously reported metabolic capability to utilize formate as an electron donor and generate intracellular CO₂, the Type II population possessed genes relevant to defense against antimicrobials and assimilating potential osmoprotectants to provide cellular stability. The accessory genome of the Mixed population was enriched in genes for multiple glycosyltransferases suggesting reduced energetic costs by adhering to mineral surfaces or to other microorganisms, and fostering a non-motile lifestyle. These results highlight the niche differentiation of distinct Methanobacterium populations to circumvent the challenges of serpentinization impacted fluids through coexistence strategies, supporting our ability to understand controls on methanogenic lifestyles and adaptations within the serpentinizing subsurface fluids of the Samail Ophiolite.

Sulfur cycling connects microbiomes and biogeochemistry in deep-sea hydrothermal plumes

In globally distributed deep-sea hydrothermal vent plumes, microbiomes are shaped by the redox energy landscapes created by reduced hydrothermal vent fluids mixing with oxidized seawater. Plumes can disperse over thousands of kilometers and their characteristics are determined by geochemical sources from vents, e.g., hydrothermal inputs, nutrients, and trace metals. However, the impacts of plume biogeochemistry on the oceans are poorly constrained due to a lack of integrated understanding of microbiomes, population genetics, and geochemistry. Here, we use microbial genomes to understand links between biogeography, evolution, and metabolic connectivity, and elucidate their impacts on biogeochemical cycling in the deep sea. Using data from 36 diverse plume samples from seven ocean basins, we show that sulfur metabolism defines the core microbiome of plumes and drives metabolic connectivity in the microbial community. Sulfur-dominated geochemistry influences energy landscapes and promotes microbial growth, while other energy sources influence local energy landscapes. We further demonstrated the consistency of links among geochemistry, function, and taxonomy. Amongst all microbial metabolisms, sulfur transformations had the highest MW-score, a measure of metabolic connectivity in microbial communities. Additionally, plume microbial populations have low diversity, short migration history, and gene-specific sweep patterns after migrating from background seawater. Selected functions include nutrient uptake, aerobic oxidation, sulfur oxidation for higher energy yields, and stress responses for adaptation. Our findings provide the ecological and evolutionary bases of change in sulfur-driven microbial communities and their population genetics in adaptation to changing geochemical gradients in the oceans.

Viruses Regulate Microbial Community Assembly Together with Environmental Factors in Acid Mine Drainage

Viruses are widespread in various ecosystems, and they play important roles in regulating the microbial community via host-virus interactions. Recently, metagenomic studies showed that there are extremely diverse viruses in different environments from the ocean to the human gut, but the influences of viral communities on microbial communities are poorly understood, especially in extreme environments. Here, we used metagenomics to characterize microbial communities and viral communities in acid mine drainage (AMD) and evaluated how viruses shape microbial community constrained by the harsh environments. Our results showed that AMD viral communities are significantly associated with the microbial communities, and viral diversity has positive correlations with microbial diversity. Viral community explained more variations of microbial community composition than environmental factors in AMD of a polymetallic mine. Moreover, we found that viruses harboring adaptive genes regulate a relative abundance of hosts under the modulation of environmental factors, such as pH. We also observed that viral diversity has significant correlations with the global properties of microbial cooccurrence networks, such as modularity. In addition, the results of null modeling analyses revealed that viruses significantly affect microbial community phylogeny and play important roles in regulating ecological processes of community assembly, such as dispersal limitation and homogenous dispersal. Together, these results revealed that AMD viruses are critical forces driving microbial network and community assembly via host-virus interactions. IMPORTANCE Viruses as mobile genetic elements play critical roles in the adaptive evolution of their hosts in extreme environments. However, how viruses further influence microbial community structure and assembly is still unclear. A recent metagenomic study observed diverse viruses unexplored in acid mine drainage, revealing the associations between the viral community and environmental factors. Here, we showed that viruses together with environmental factors can constrain the relative abundance of host and microbial community assembly in AMD of copper mines and polymetallic mines. Our results highlight the importance of viruses in shaping the microbial community from the individual host level to the community level.

Structure-informed microbial population genetics elucidate selective pressures that shape protein evolution

Comprehensive sampling of natural genetic diversity with metagenomics enables highly resolved insights into the interplay between ecology and evolution. However, resolving adaptive, neutral, or purifying processes of evolution from intrapopulation genomic variation remains a challenge, partly due to the sole reliance on gene sequences to interpret variants. Here, we describe an approach to analyze genetic variation in the context of predicted protein structures and apply it to a marine microbial population within the SAR11 subclade 1a.3.V, which dominates low-latitude surface oceans. Our analyses reveal a tight association between genetic variation and protein structure. In a central gene in nitrogen metabolism, we observe decreased occurrence of nonsynonymous variants from ligand-binding sites as a function of nitrate concentrations, revealing genetic targets of distinct evolutionary pressures maintained by nutrient availability. Our work yields insights into the governing principles of evolution and enables structure-aware investigations of microbial population genetics.

Panguiarchaeum symbiosum, a potential hyperthermophilic symbiont in the TACK superphylum

The biology of Korarchaeia remains elusive due to the lack of genome representatives. Here, we reconstruct 10 closely related metagenome-assembled genomes from hot spring habitats and place them into a single species, proposed herein as Panguiarchaeum symbiosum. Functional investigation suggests that Panguiarchaeum symbiosum is strictly anaerobic and grows exclusively in thermal habitats by fermenting peptides coupled with sulfide and hydrogen production to dispose of electrons. Due to its inability to biosynthesize archaeal membranes, amino acids, and purines, this species likely exists in a symbiotic lifestyle similar to DPANN archaea. Population metagenomics and metatranscriptomic analyses demonstrated that genes associated with amino acid/peptide uptake and cell attachment exhibited positive selection and were highly expressed, supporting the proposed proteolytic catabolism and symbiotic lifestyle. Our study sheds light on the metabolism, evolution, and potential symbiotic lifestyle of Panguiarchaeum symbiosum, which may be a unique host-dependent archaeon within the TACK superphylum.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

Virioplankton assemblages from challenger deep, the deepest place in the oceans

Hadal ocean biosphere, that is, the deepest part of the world's oceans, harbors a unique microbial community, suggesting a potential uncovered co-occurring virioplankton assemblage. Herein, we reveal the unique virioplankton assemblages of the Challenger Deep, comprising 95,813 non-redundant viral contigs from the surface to the hadal zone. Almost all of the dominant viral contigs in the hadal zone were unclassified, potentially related to Alteromonadales and Oceanospirillales. 2,586 viral auxiliary metabolic genes from 132 different KEGG orthologous groups were mainly related to the carbon, nitrogen, sulfur, and arsenic metabolism. Lysogenic viral production and integrase genes were augmented in the hadal zone, suggesting the prevalence of viral lysogenic life strategy. Abundant rve genes in the hadal zone, which function as transposase in the caudoviruses, further suggest the prevalence of viral-mediated horizontal gene transfer. This study provides fundamental insights into the virioplankton assemblages of the hadal zone, reinforcing the necessity of incorporating virioplankton into the hadal biogeochemical cycles.

Marine Bacteria and Omic Approaches: A Novel and Potential Repository for Bioremediation Assessment

Marine environments accommodating diverse assortments of life constitute a great pool of differentiated natural resources. The cumulative need to remedy unpropitious effects of anthropogenic activities on estuaries, and coastal marine ecosystems has propelled the development of effective bioremediation strategies. Marine bacteria producing biosurfactants are promising agents for bio-remediating oil pollution in marine environments, making them prospective candidates for enhancing oil recovery. Molecular omics technologies are considered an emerging field of research in ecological and diversity assessment owing to their utility in environmental surveillance and bioremediation of polluted sites. A thorough literature review was undertaken to understand the applicability of different omic techniques employed for bioremediation assessment using marine bacteria. This review further establishes that for bioremediation of environmental pollutants (i.e., heavy metals, hydrocarbons, xenobiotic and numerous recalcitrant compounds), organisms isolated from marine environments can be better utilized for their removal. The literature survey shows that omics approaches can provide exemplary knowledge about microbial communities and their role in the bioremediation of environmental pollutants. This review centres on applications of marine bacteria in enhanced bioremediation, utilizing the omics approaches that can be a vital biological contrivance in environmental monitoring to tackle environmental degradation. The paper aims to identify the gaps in investigations involving marine bacteria to help researchers, ecologists, and decision-makers to develop a holistic understanding regarding their utility in bioremediation assessment.

Microbial diversity in extreme environments

A wide array of microorganisms, including many novel, phylogenetically deeply rooted taxa, survive and thrive in extreme environments. These unique and reduced-complexity ecosystems offer a tremendous opportunity for studying the structure, function and evolution of natural microbial communities. Marker gene surveys have resolved patterns and ecological drivers of these extremophile assemblages, revealing a vast uncultured microbial diversity and the often predominance of archaea in the most extreme conditions. New omics studies have uncovered linkages between community function and environmental variables, and have enabled discovery and genomic characterization of major new lineages that substantially expand microbial diversity and change the structure of the tree of life. These efforts have significantly advanced our understanding of the diversity, ecology and evolution of microorganisms populating Earth's extreme environments, and have facilitated the exploration of microbiota and processes in more complex ecosystems.

Eco-evolutionary significance of domesticated retroelements in microbial genomes

Since the first discovery of reverse transcriptase in bacteria, and later in archaea, bacterial and archaeal retroelements have been defined by their common enzyme that coordinates diverse functions. Yet, evolutionary refinement has produced distinct retroelements across the tree of microbial life that are perhaps best described in terms of their programmed RNA-a compact sequence that preserves core information for a sophisticated mechanism. From this perspective, reverse transcriptase has been selected as the modular tool for carrying out nature's instructions in various RNA templates. Beneficial retroelements-those that can provide a fitness advantage to their host-evolved to their extant forms in a wide array of microorganisms and their viruses, spanning nearly all habitats. Within each specialized retroelement class, several universal features seem to be shared across diverse taxa, while specific functional and mechanistic insights are based on only a few model retroelement systems from clinical isolates. Currently, little is known about the diversity of cellular functions and ecological significance of retroelements across different biomes. With increasing availability of isolate, metagenome-assembled, and single-amplified genomes, the taxonomic and functional breadth of prokaryotic retroelements is coming into clearer view. This review explores the recently characterized classes of beneficial, yet accessory retroelements of bacteria and archaea. We describe how these specialized mechanisms exploit a form of fixed mobility, whereby the retroelements do not appear to proliferate selfishly throughout the genome. Moreover, we discuss computational approaches for systematic identification of retroelements from vast sequence repositories and highlight recent discoveries in terms of their apparent distribution and ecological significance in nature. Lastly, we present a new perspective on the eco-evolutionary significance of these genetic elements in marine bacteria and demonstrate approaches that enable the characterization of their environmental diversity through metagenomics.

Integrative Genomics Sheds Light on Evolutionary Forces Shaping the Acidithiobacillia Class Acidophilic Lifestyle

Extreme acidophiles thrive in environments rich in protons (pH values <3) and often high levels of dissolved heavy metals. They are distributed across the three domains of the Tree of Life including members of the Proteobacteria. The Acidithiobacillia class is formed by the neutrophilic genus Thermithiobacillus along with the extremely acidophilic genera Fervidacidithiobacillus, Igneacidithiobacillus, Ambacidithiobacillus, and Acidithiobacillus. Phylogenomic reconstruction revealed a division in the Acidithiobacillia class correlating with the different pH optima that suggested that the acidophilic genera evolved from an ancestral neutrophile within the Acidithiobacillia. Genes and mechanisms denominated as "first line of defense" were key to explaining the Acidithiobacillia acidophilic lifestyle including preventing proton influx that allows the cell to maintain a near-neutral cytoplasmic pH and differ from the neutrophilic Acidithiobacillia ancestors that lacked these systems. Additional differences between the neutrophilic and acidophilic Acidithiobacillia included the higher number of gene copies in the acidophilic genera coding for "second line of defense" systems that neutralize and/or expel protons from cell. Gain of genes such as hopanoid biosynthesis involved in membrane stabilization at low pH and the functional redundancy for generating an internal positive membrane potential revealed the transition from neutrophilic properties to a new acidophilic lifestyle by shaping the Acidithiobacillaceae genomic structure. The presence of a pool of accessory genes with functional redundancy provides the opportunity to "hedge bet" in rapidly changing acidic environments. Although a core of mechanisms for acid resistance was inherited vertically from an inferred neutrophilic ancestor, the majority of mechanisms, especially those potentially involved in resistance to extremely low pH, were obtained from other extreme acidophiles by horizontal gene transfer (HGT) events.

inStrain profiles population microdiversity from metagenomic data and sensitively detects shared microbial strains

Coexisting microbial cells of the same species often exhibit genetic variation that can affect phenotypes ranging from nutrient preference to pathogenicity. Here we present inStrain, a program that uses metagenomic paired reads to profile intra-population genetic diversity (microdiversity) across whole genomes and compares microbial populations in a microdiversity-aware manner, greatly increasing the accuracy of genomic comparisons when benchmarked against existing methods. We use inStrain to profile >1,000 fecal metagenomes from newborn premature infants and find that siblings share significantly more strains than unrelated infants, although identical twins share no more strains than fraternal siblings. Infants born by cesarean section harbor Klebsiella with significantly higher nucleotide diversity than infants delivered vaginally, potentially reflecting acquisition from hospital rather than maternal microbiomes. Genomic loci that show diversity in individual infants include variants found between other infants, possibly reflecting inoculation from diverse hospital-associated sources. inStrain can be applied to any metagenomic dataset for microdiversity analysis and rigorous strain comparison.

Metagenomic profiling of host-associated bacteria from 8 datasets of the red alga Porphyra purpurea with MetaPhlAn3

Microbial communities play fundamental roles in association with marine algae; in fact, they are recognized to be actively involved in growth and morphogenesis of the algae. Porphyra purpurea is a red alga commonly found in the intertidal zone with a high economic value, however little is known about the bacterial species associated with this genus. Here we report the bacterial-associated diversity of P. purpurea in four different localities (Ireland, Italy United Kingdom and the USA) from analyzing eight publicly available metagenomic datasets. These were analyzed with Methaplan3 to identify the putative bacterial taxonomies and their relative abundances. Furthermore, we compared these results to the 16S rRNA metagenomic analysis pipeline of the MGnify database to evaluate both methods. Kraken2 was used to verify and support the results, as a complementary classification method to Metaphlan3. This approach highlighted the different taxonomic resolution of a 16S rRNA OTU-based method compared to the pan-genome approach deployed by Metaphlan3 and complemented by Kraken2. The results presented here provide valuable preliminary data on the putative host-associated bacterial species of P. purpurea.

Metapangenomics of the oral microbiome provides insights into habitat adaptation and cultivar diversity

BACKGROUND

The increasing availability of microbial genomes and environmental shotgun metagenomes provides unprecedented access to the genomic differences within related bacteria. The human oral microbiome with its diverse habitats and abundant, relatively well-characterized microbial inhabitants presents an opportunity to investigate bacterial population structures at an ecosystem scale.

RESULTS

Here, we employ a metapangenomic approach that combines public genomes with Human Microbiome Project (HMP) metagenomes to study the diversity of microbial residents of three oral habitats: tongue dorsum, buccal mucosa, and supragingival plaque. For two exemplar taxa, Haemophilus parainfluenzae and the genus Rothia, metapangenomes reveal distinct genomic groups based on shared genome content. H. parainfluenzae genomes separate into three distinct subgroups with differential abundance between oral habitats. Functional enrichment analyses identify an operon encoding oxaloacetate decarboxylase as diagnostic for the tongue-abundant subgroup. For the genus Rothia, grouping by shared genome content recapitulates species-level taxonomy and habitat preferences. However, while most R. mucilaginosa are restricted to the tongue as expected, two genomes represent a cryptic population of R. mucilaginosa in many buccal mucosa samples. For both H. parainfluenzae and the genus Rothia, we identify not only limitations in the ability of cultivated organisms to represent populations in their native environment, but also specifically which cultivar gene sequences are absent or ubiquitous.

CONCLUSIONS

Our findings provide insights into population structure and biogeography in the mouth and form specific hypotheses about habitat adaptation. These results illustrate the power of combining metagenomes and pangenomes to investigate the ecology and evolution of bacteria across analytical scales.

Evolution of Predicted Acid Resistance Mechanisms in the Extremely Acidophilic Leptospirillum Genus

Organisms that thrive in extremely acidic environments (≤pH 3.5) are of widespread importance in industrial applications, environmental issues, and evolutionary studies. Leptospirillum spp. constitute the only extremely acidophilic microbes in the phylogenetically deep-rooted bacterial phylum Nitrospirae. Leptospirilli are Gram-negative, obligatory chemolithoautotrophic, aerobic, ferrous iron oxidizers. This paper predicts genes that Leptospirilli use to survive at low pH and infers their evolutionary trajectory. Phylogenetic and other bioinformatic approaches suggest that these genes can be classified into (i) "first line of defense", involved in the prevention of the entry of protons into the cell, and (ii) neutralization or expulsion of protons that enter the cell. The first line of defense includes potassium transporters, predicted to form an inside positive membrane potential, spermidines, hopanoids, and Slps (starvation-inducible outer membrane proteins). The "second line of defense" includes proton pumps and enzymes that consume protons. Maximum parsimony, clustering methods, and gene alignments are used to infer the evolutionary trajectory that potentially enabled the ancestral Leptospirillum to transition from a postulated circum-neutral pH environment to an extremely acidic one. The hypothesized trajectory includes gene gains/loss events driven extensively by horizontal gene transfer, gene duplications, gene mutations, and genomic rearrangements.

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Single-amino acid variants reveal evolutionary processes that shape the biogeography of a global SAR11 subclade

Members of the SAR11 order Pelagibacterales dominate the surface oceans. Their extensive diversity challenges emerging operational boundaries defined for microbial 'species' and complicates efforts of population genetics to study their evolution. Here, we employed single-amino acid variants (SAAVs) to investigate ecological and evolutionary forces that maintain the genomic heterogeneity within ubiquitous SAR11 populations we accessed through metagenomic read recruitment using a single isolate genome. Integrating amino acid and protein biochemistry with metagenomics revealed that systematic purifying selection against deleterious variants governs non-synonymous variation among very closely related populations of SAR11. SAAVs partitioned metagenomes into two main groups matching large-scale oceanic current temperatures, and six finer proteotypes that connect distant oceanic regions. These findings suggest that environmentally-mediated selection plays a critical role in the journey of cosmopolitan surface ocean microbial populations, and the idea 'everything is everywhere but the environment selects' has credence even at the finest resolutions.

Beyond cells - The virome in the human holobiont

Viromics, or viral metagenomics, is a relatively new and burgeoning field of research that studies the complete collection of viruses forming part of the microbiota in any given niche. It has strong foundations rooted in over a century of discoveries in the field of virology and recent advances in molecular biology and sequencing technologies. Historically, most studies have deconstructed the concept of viruses into a simplified perception of viral agents as mere pathogens, which demerits the scope of large-scale viromic analyses. Viruses are, in fact, much more than regular parasites. They are by far the most dynamic and abundant entity and the greatest killers on the planet, as well as the most effective geo-transforming genetic engineers and resource recyclers, acting on all life strata in any habitat. Yet, most of this uncanny viral world remains vastly unexplored to date, greatly hindered by the bewildering complexity inherent to such studies and the methodological and conceptual limitations. Viromic studies are just starting to address some of these issues but they still lag behind microbial metagenomics. In recent years, however, higher-throughput analysis and resequencing have rekindled interest in a field that is just starting to show its true potential. In this review, we take a look at the scientific and technological developments that led to the advent of viral and bacterial metagenomics with a particular, but not exclusive, focus on human viromics from an ecological perspective. We also address some of the most relevant challenges that current viral studies face and ponder on the future directions of the field.

MetaSMC: a coalescent-based shotgun sequence simulator for evolving microbial populations

MOTIVATION

High-throughput sequencing technology has revolutionized the study of metagenomics and cancer evolution. In a relatively simple environment, a metagenomics sequencing data is dominated by a few species. By analyzing the alignment of reads from microbial species, single nucleotide polymorphisms can be discovered and the evolutionary history of the populations can be reconstructed. The ever-increasing read length will allow more detailed analysis about the evolutionary history of microbial or tumor cell population. A simulator of shotgun sequences from such populations will be helpful in the development or evaluation of analysis algorithms.

RESULTS

Here, we described an efficient algorithm, MetaSMC, which simulates reads from evolving microbial populations. Based on the coalescent theory, our simulator supports all evolutionary scenarios supported by other coalescent simulators. In addition, the simulator supports various substitution models, including Jukes-Cantor, HKY85 and generalized time-reversible models. The simulator also supports mutator phenotypes by allowing different mutation rates and substitution models in different subpopulations. Our algorithm ignores unnecessary chromosomal segments and thus is more efficient than standard coalescent when recombination is frequent. We showed that the process behind our algorithm is equivalent to Sequentially Markov Coalescent with an incomplete sample. The accuracy of our algorithm was evaluated by summary statistics and likelihood curves derived from Monte Carlo integration over large number of random genealogies.

AVAILABILITY AND IMPLEMENTATION

MetaSMC is written in C. The source code is available at https://github.com/tarjxvf/metasmc.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparative genomics of 16 Microbacterium spp. that tolerate multiple heavy metals and antibiotics

A total of 16 different strains of Microbacterium spp. were isolated from contaminated soil and enriched on the carcinogen, hexavalent chromium [Cr(VI)]. The majority of the isolates (11 of the 16) were able to tolerate concentrations (0.1 mM) of cobalt, cadmium, and nickel, in addition to Cr(VI) (0.5–20 mM). Interestingly, these bacteria were also able to tolerate three different antibiotics (ranges: ampicillin 0–16 μg ml⁻¹, chloramphenicol 0–24 μg ml⁻¹, and vancomycin 0–24 μg ml⁻¹). To gain genetic insight into these tolerance pathways, the genomes of these isolates were assembled and annotated. The genomes of these isolates not only have some shared genes (core genome) but also have a large amount of variability. The genomes also contained an annotated Cr(VI) reductase (chrR) that could be related to Cr(VI) reduction. Further, various heavy metal tolerance (e.g., Co/Zn/Cd efflux system) and antibiotic resistance genes were identified, which provide insight into the isolates’ ability to tolerate metals and antibiotics. Overall, these isolates showed a wide range of tolerances to heavy metals and antibiotics and genetic diversity, which was likely required of this population to thrive in a contaminated environment.

Review of Constructed Wetlands for Acid Mine Drainage Treatment

The mining industry is the major producer of acid mine drainage (AMD). The problem of AMD concerns at active and abandoned mine sites. Acid mine drainage needs to be treated since it can contaminate surface water. Constructed wetlands (CW), a passive treatment technology, combines naturally-occurring biogeochemical, geochemical, and physical processes. This technology can be used for the long-term remediation of AMD. The challenge is to overcome some factors, for instance, chemical characteristics of AMD such a high acidity and toxic metals concentrations, to achieve efficient CW systems. Design criteria, conformational arrangements, and careful selection of each component must be considered to achieve the treatment. The main objective of this review is to summarize the current advances, applications, and the prevalent difficulties and opportunities to apply the CW technology for AMD treatment. According to the cited literature, sub-surface CW (SS-CW) systems are suggested for an efficient AMD treatment. The synergistic interactions between CW components determine heavy metal removal from water solution. The microorganism-plant interaction is considered the most important since it implies symbiosis mechanisms for heavy metal removal and tolerance. In addition, formation of litter and biofilm layers contributes to heavy metal removal by adsorption mechanisms. The addition of organic amendments to the substrate material and AMD bacterial consortium inoculation are some of the strategies to improve heavy metal removal. Adequate experimental design from laboratory to full scale systems need to be used to optimize equilibria between CW components selection and construction and operational costs. The principal limitations for CW treating AMD is the toxicity effect that heavy metals produce on CW plants and microorganisms. However, these aspects can be solved partially by choosing carefully constructed wetlands components suitable for the AMD characteristics. From the economic point of view, a variety of factors affects the cost of constructed wetlands, such as: detention time, treatment goals, media type, pretreatment type, number of cells, source, and availability of gravel media, and land requirements, among others.

Biological species in the viral world

Due to their dependence on cellular organisms for metabolism and replication, viruses are typically named and assigned to species according to their genome structure and the original host that they infect. But because viruses often infect multiple hosts and the numbers of distinct lineages within a host can be vast, their delineation into species is often dictated by arbitrary sequence thresholds, which are highly inconsistent across lineages. Here we apply an approach to determine the boundaries of viral species based on the detection of gene flow within populations, thereby defining viral species according to the biological species concept (BSC). Despite the potential for gene transfer between highly divergent genomes, viruses, like the cellular organisms they infect, assort into reproductively isolated groups and can be organized into biological species. This approach revealed that BSC-defined viral species are often congruent with the taxonomic partitioning based on shared gene contents and host tropism, and that bacteriophages can similarly be classified in biological species. These results open the possibility to use a single, universal definition of species that is applicable across cellular and acellular lifeforms.

Microdiversity of an Abundant Terrestrial Bacterium Encompasses Extensive Variation in Ecologically Relevant Traits

Much genetic diversity within a bacterial community is likely obscured by microdiversity within operational taxonomic units (OTUs) defined by 16S rRNA gene sequences. However, it is unclear how variation within this microdiversity influences ecologically relevant traits. Here, we employ a multifaceted approach to investigate microdiversity within the dominant leaf litter bacterium, Curtobacterium, which comprises 7.8% of the bacterial community at a grassland site undergoing global change manipulations. We use cultured bacterial isolates to interpret metagenomic data, collected in situ over 2 years, together with lab-based physiological assays to determine the extent of trait variation within this abundant OTU. The response of Curtobacterium to seasonal variability and the global change manipulations, specifically an increase in relative abundance under decreased water availability, appeared to be conserved across six Curtobacterium lineages identified at this site. Genomic and physiological analyses in the lab revealed that degradation of abundant polymeric carbohydrates within leaf litter, cellulose and xylan, is nearly universal across the genus, which may contribute to its high abundance in grassland leaf litter. However, the degree of carbohydrate utilization and temperature preference for this degradation varied greatly among clades. Overall, we find that traits within Curtobacterium are conserved at different phylogenetic depths. We speculate that similar to bacteria in marine systems, diverse microbes within this taxon may be structured in distinct ecotypes that are key to understanding Curtobacterium abundance and distribution in the environment.

Despite the plummeting costs of sequencing, characterizing the fine-scale genetic diversity of a microbial community—and interpreting its functional importance—remains a challenge. Indeed, most studies, particularly studies of soil, assess community composition at a broad genetic level by classifying diversity into taxa (OTUs) defined by 16S rRNA sequence similarity. However, these classifications potentially obscure variation in traits that result in fine-scale ecological differentiation among closely related strains. Here, we investigated “microdiversity” in a highly diverse and poorly characterized soil system (leaf litter in a southern Californian grassland). We focused on the most abundant bacterium, Curtobacterium, which by standard methods is grouped into only one OTU. We find that the degree of carbohydrate usage and temperature preference vary within the OTU, whereas its responses to changes in precipitation are relatively uniform. These results suggest that microdiversity may be key to understanding how soil bacterial diversity is linked to ecosystem functioning.

Geobiological feedbacks and the evolution of thermoacidophiles

Oxygen-dependent microbial oxidation of sulfur compounds leads to the acidification of natural waters. How acidophiles and their acidic habitats evolved, however, is largely unknown. Using 16S rRNA gene abundance and composition data from 72 hot springs in Yellowstone National Park, Wyoming, we show that hyperacidic (pH<3.0) hydrothermal ecosystems are dominated by a limited number of archaeal lineages with an inferred ability to respire O₂. Phylogenomic analyses of 584 existing archaeal genomes revealed that hyperacidophiles evolved independently multiple times within the Archaea, each coincident with the emergence of the ability to respire O₂, and that these events likely occurred in the recent evolutionary past. Comparative genomic analyses indicated that archaeal thermoacidophiles from independent lineages are enriched in similar protein-coding genes, consistent with convergent evolution aided by horizontal gene transfer. Because the generation of acidic environments and their successful habitation characteristically require O₂, these results suggest that thermoacidophilic Archaea and the acidity of their habitats co-evolved after the evolution of oxygenic photosynthesis. Moreover, it is likely that dissolved O₂ concentrations in thermal waters likely did not reach levels capable of sustaining aerobic thermoacidophiles and their acidifying activity until ~0.8 Ga, when present day atmospheric levels were reached, a time period that is supported by our estimation of divergence times for archaeal thermoacidophilic clades.

Microbial copper resistance: importance in biohydrometallurgy

Industrial biomining has been extensively used for many years to recover valuable metals such as copper, gold, uranium and others. Furthermore, microorganisms involved in these processes can also be used to bioremediate places contaminated with acid and metals. These uses are possible due to the great metal resistance that these extreme acidophilic microorganisms possess. In this review, the most recent findings related to copper resistance mechanisms of bacteria and archaea related to biohydrometallurgy are described. The recent search for novel metal resistance determinants is not only of scientific interest but also of industrial importance, as reflected by the genomic sequencing of microorganisms present in mining operations and the search of those bacteria with extreme metal resistance to improve the extraction processes used by the biomining companies.

Genome-wide selective sweeps and gene-specific sweeps in natural bacterial populations

Multiple models describe the formation and evolution of distinct microbial phylogenetic groups. These evolutionary models make different predictions regarding how adaptive alleles spread through populations and how genetic diversity is maintained. Processes predicted by competing evolutionary models, for example, genome-wide selective sweeps vs gene-specific sweeps, could be captured in natural populations using time-series metagenomics if the approach were applied over a sufficiently long time frame. Direct observations of either process would help resolve how distinct microbial groups evolve. Here, from a 9-year metagenomic study of a freshwater lake (2005-2013), we explore changes in single-nucleotide polymorphism (SNP) frequencies and patterns of gene gain and loss in 30 bacterial populations. SNP analyses revealed substantial genetic heterogeneity within these populations, although the degree of heterogeneity varied by >1000-fold among populations. SNP allele frequencies also changed dramatically over time within some populations. Interestingly, nearly all SNP variants were slowly purged over several years from one population of green sulfur bacteria, while at the same time multiple genes either swept through or were lost from this population. These patterns were consistent with a genome-wide selective sweep in progress, a process predicted by the 'ecotype model' of speciation but not previously observed in nature. In contrast, other populations contained large, SNP-free genomic regions that appear to have swept independently through the populations prior to the study without purging diversity elsewhere in the genome. Evidence for both genome-wide and gene-specific sweeps suggests that different models of bacterial speciation may apply to different populations coexisting in the same environment.

Microbial Ecology and Evolution in the Acid Mine Drainage Model System

Acid mine drainage (AMD) is a unique ecological niche for acid- and toxic-metals-adapted microorganisms. These low-complexity systems offer a special opportunity for the ecological and evolutionary analyses of natural microbial assemblages. The last decade has witnessed an unprecedented interest in the study of AMD communities using 16S rRNA high-throughput sequencing and community genomic and postgenomic methodologies, significantly advancing our understanding of microbial diversity, community function, and evolution in acidic environments. This review describes new data on AMD microbial ecology and evolution, especially dynamics of microbial diversity, community functions, and population genomes, and further identifies gaps in our current knowledge that future research, with integrated applications of meta-omics technologies, will fill.

Draft genome sequence of chloride-tolerant Leptospirillum ferriphilum Sp-Cl from industrial bioleaching operations in northern Chile

Leptospirillum ferriphilum Sp-Cl is a Gram negative, thermotolerant, curved, rod-shaped bacterium, isolated from an industrial bioleaching operation in northern Chile, where chalcocite is the major copper mineral and copper hydroxychloride atacamite is present in variable proportions in the ore. This strain has unique features as compared to the other members of the species, namely resistance to elevated concentrations of chloride, sulfate and metals. Basic microbiological features and genomic properties of this biotechnologically relevant strain are described in this work. The 2,475,669 bp draft genome is arranged into 74 scaffolds of 74 contigs. A total of 48 RNA genes and 2,834 protein coding genes were predicted from its annotation; 55 % of these were assigned a putative function. Release of the genome sequence of this strain will provide further understanding of the mechanisms used by acidophilic bacteria to endure high osmotic stress and high chloride levels and of the role of chloride-tolerant iron-oxidizers in industrial bioleaching operations.

Illuminating structural proteins in viral "dark matter" with metaproteomics

Viruses are ecologically important, yet environmental virology is limited by dominance of unannotated genomic sequences representing taxonomic and functional "viral dark matter." Although recent analytical advances are rapidly improving taxonomic annotations, identifying functional dark matter remains problematic. Here, we apply paired metaproteomics and dsDNA-targeted metagenomics to identify 1,875 virion-associated proteins from the ocean. Over one-half of these proteins were newly functionally annotated and represent abundant and widespread viral metagenome-derived protein clusters (PCs). One primarily unannotated PC dominated the dataset, but structural modeling and genomic context identified this PC as a previously unidentified capsid protein from multiple uncultivated tailed virus families. Furthermore, four of the five most abundant PCs in the metaproteome represent capsid proteins containing the HK97-like protein fold previously found in many viruses that infect all three domains of life. The dominance of these proteins within our dataset, as well as their global distribution throughout the world's oceans and seas, supports prior hypotheses that this HK97-like protein fold is the most abundant biological structure on Earth. Together, these culture-independent analyses improve virion-associated protein annotations, facilitate the investigation of proteins within natural viral communities, and offer a high-throughput means of illuminating functional viral dark matter.

Migration and horizontal gene transfer divide microbial genomes into multiple niches

Horizontal gene transfer is central to microbial evolution, because it enables genetic regions to spread horizontally through diverse communities. However, how gene transfer exerts such a strong effect is not understood. Here we develop an eco-evolutionary model and show how genetic transfer, even when rare, can transform the evolution and ecology of microbes. We recapitulate existing models, which suggest that asexual reproduction will overpower horizontal transfer and greatly limit its effects. We then show that allowing immigration completely changes these predictions. With migration, the rates and impacts of horizontal transfer are greatly increased, and transfer is most frequent for loci under positive natural selection. Our analysis explains how ecologically important loci can sweep through competing strains and species. In this way, microbial genomes can evolve to become ecologically diverse where different genomic regions encode for partially overlapping, but distinct, ecologies. Under these conditions ecological species do not exist, because genes, not species, inhabit niches.

Horizontal gene transfer is central to microbial evolution. Here, the authors develop an eco-evolutionary model and show that migration can greatly promote horizontal gene transfer, which explains how ecologically-important loci can sweep through the species in a microbial community.

Metabolic and genomic analysis elucidates strain-level variation in Microbacterium spp. isolated from chromate contaminated sediment

Hexavalent chromium [Cr(VI)] is a soluble carcinogen that has caused widespread contamination of soil and water in many industrial nations. Bacteria have the potential to aid remediation as certain strains can catalyze the reduction of Cr(VI) to insoluble and less toxic Cr(III). Here, we examine Cr(VI) reducing Microbacterium spp. (Cr-K1W, Cr-K20, Cr-K29, and Cr-K32) isolated from contaminated sediment (Seymore, Indiana) and show varying chromate responses despite the isolates' phylogenetic similarity (i.e., identical 16S rRNA gene sequences). Detailed analysis identified differences based on genomic metabolic potential, growth and general metabolic capabilities, and capacity to resist and reduce Cr(VI). Taken together, the discrepancies between the isolates demonstrate the complexity inter-strain variation can have on microbial physiology and related biogeochemical processes.

Anvi'o: an advanced analysis and visualization platform for 'omics data

Advances in high-throughput sequencing and ‘omics technologies are revolutionizing studies of naturally occurring microbial communities. Comprehensive investigations of microbial lifestyles require the ability to interactively organize and visualize genetic information and to incorporate subtle differences that enable greater resolution of complex data. Here we introduce anvi’o, an advanced analysis and visualization platform that offers automated and human-guided characterization of microbial genomes in metagenomic assemblies, with interactive interfaces that can link ‘omics data from multiple sources into a single, intuitive display. Its extensible visualization approach distills multiple dimensions of information about each contig, offering a dynamic and unified work environment for data exploration, manipulation, and reporting. Using anvi’o, we re-analyzed publicly available datasets and explored temporal genomic changes within naturally occurring microbial populations through de novo characterization of single nucleotide variations, and linked cultivar and single-cell genomes with metagenomic and metatranscriptomic data. Anvi’o is an open-source platform that empowers researchers without extensive bioinformatics skills to perform and communicate in-depth analyses on large ‘omics datasets.

Metagenomic reconstructions of bacterial CRISPR loci constrain population histories

Bacterial CRISPR-Cas systems provide insight into recent population history because they rapidly incorporate, in a unidirectional manner, short fragments (spacers) from coexisting infective virus populations into host chromosomes. Immunity is achieved by sequence identity between transcripts of spacers and their targets. Here, we used metagenomics to study the stability and dynamics of the type I-E CRISPR-Cas locus of Leptospirillum group II bacteria in biofilms sampled over 5 years from an acid mine drainage (AMD) system. Despite recovery of 452,686 spacers from CRISPR amplicons and metagenomic data, rarefaction curves of spacers show no saturation. The vast repertoire of spacers is attributed to phage/plasmid population diversity and retention of old spacers, despite rapid evolution of the targeted phage/plasmid genome regions (proto-spacers). The oldest spacers (spacers found at the trailer end) are conserved for at least 5 years, and 12% of these retain perfect or near-perfect matches to proto-spacer targets. The majority of proto-spacer regions contain an AAG proto-spacer adjacent motif (PAM). Spacers throughout the locus target the same phage population (AMDV1), but there are blocks of consecutive spacers without AMDV1 target sequences. Results suggest long-term coexistence of Leptospirillum with AMDV1 and periods when AMDV1 was less dominant. Metagenomics can be applied to millions of cells in a single sample to provide an extremely large spacer inventory, allow identification of phage/plasmids and enable analysis of previous phage/plasmid exposure. Thus, this approach can provide insights into prior bacterial environment and genetic interplay between hosts and their viruses.

Microbial Speciation

What are species? How do they arise? These questions are not easy to answer and have been particularly controversial in microbiology. Yet, for those microbiologists studying environmental questions or dealing with clinical issues, the ability to name and recognize species, widely considered the fundamental units of ecology, can be practically useful. On a more fundamental level, the speciation problem, the focus here, is more mechanistic and conceptual. What is the origin of microbial species, and what evolutionary and ecological mechanisms keep them separate once they begin to diverge? To what extent are these mechanisms universal across diverse types of microbes, and more broadly across the entire the tree of life? Here, we propose that microbial speciation must be viewed in light of gene flow, which defines units of genetic similarity, and of natural selection, which defines units of phenotype and ecological function. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and population genomics studies. These studies suggest a continuous "speciation spectrum," which microbial populations traverse in different ways depending on their balance of gene flow and natural selection.

Proteomic tools to decipher microbial community structure and functioning

Recent advances in microbial ecology allow studying microorganisms in their environment, without laboratory cultivation, in order to get access to the large uncultivable microbial community. With this aim, environmental proteomics has emerged as an appropriate complementary approach to metagenomics providing information on key players that carry out main metabolic functions and addressing the adaptation capacities of living organisms in situ. In this review, a wide range of proteomic approaches applied to investigate the structure and functioning of microbial communities as well as recent examples of such studies are presented.

Toxic metal resistance in biofilms: diversity of microbial responses and their evolution

Since biofilms are an important issue in the fields of medicine and health, several recent microbiological studies have focused on their formation and their contribution to toxic compound resistance mechanisms. In this review, we describe how metals impact biofilm formation and resistance, and how biofilms can help cells resist toxic metals. First, the organic matrix acts as a barrier isolating the cells from many environmental stresses. Secondly, the metabolism of the cells changes, and a slowly-growing or non-growing sub-population of cells known as persisters emerges. Thirdly, in the case of multispecies biofilms, metabolic interactions are developed, allowing cells to be more persistent or to have greater capacity to survive than a single species biofilm. Finally, we discuss how the high density of the cells may promote horizontal gene transfer processes, resulting in the acquisition of new features. All these crucial mechanisms enable microorganisms to survive and colonize toxic environments, and probably accelerate ongoing evolutionary processes.

Mapping the inner workings of the microbiome: genomic- and metagenomic-based study of metabolism and metabolic interactions in the human microbiome

The human gut microbiome is a major contributor to human metabolism and health, yet the metabolic processes that are carried out by various community members, the way these members interact with each other and with the host, and the impact of such interactions on the overall metabolic machinery of the microbiome have not yet been mapped. Here, we discuss recent efforts to study the metabolic inner workings of this complex ecosystem. We will specifically highlight two interrelated lines of work, the first aiming to deconvolve the microbiome and to characterize the metabolic capacity of various microbiome species and the second aiming to utilize computational modeling to infer and study metabolic interactions between these species.

Comparative genomics of planktonic Flavobacteriaceae from the Gulf of Maine using metagenomic data

BACKGROUND

The Gulf of Maine is an important biological province of the Northwest Atlantic with high productivity year round. From an environmental Sanger-based metagenome, sampled in summer and winter, we were able to assemble and explore the partial environmental genomes of uncultured members of the class Flavobacteria. Each of the environmental genomes represents organisms that compose less than 1% of the total microbial metagenome.

RESULTS

Four partial environmental genomes were assembled with varying degrees of estimated completeness (37%-84% complete) and were analyzed from a perspective of gathering information regarding niche partitioning between co-occurring organisms. Comparative genomics revealed potentially important niche partitioning genomic variations, including iron transporters and genes associated with cell attachment and polymer degradation. Analysis of large syntenic regions helped reveal potentially ecologically relevant variations for Flavobacteriaceae in the Gulf of Maine, such as arginine biosynthesis, and identify a putative genomic island incorporating novel exogenous genes from the environment.

CONCLUSIONS

Biogeographic analysis revealed flavobacteria species with distinct abundance patterns suggesting the presence of local blooms relative to the other species, as well as seasonally selected organisms. The analysis of genomic content for the Gulf of Maine Flavobacteria supports the hypothesis of a particle-associated lifestyle and specifically highlights a number of putative coding sequences that may play a role in the remineralization of particulate organic matter. And lastly, analysis of the underlying sequences for each assembled genome revealed seasonal and nonseasonal variants of specific genes implicating a dynamic interaction between individuals within the species.

Functionally relevant diversity of closely related Nitrospira in activated sludge

Nitrospira are chemolithoautotrophic nitrite-oxidizing bacteria that catalyze the second step of nitrification in most oxic habitats and are important for excess nitrogen removal from sewage in wastewater treatment plants (WWTPs). To date, little is known about their diversity and ecological niche partitioning within complex communities. In this study, the fine-scale community structure and function of Nitrospira was analyzed in two full-scale WWTPs as model ecosystems. In Nitrospira-specific 16S rRNA clone libraries retrieved from each plant, closely related phylogenetic clusters (16S rRNA identities between clusters ranged from 95.8% to 99.6%) within Nitrospira lineages I and II were found. Newly designed probes for fluorescence in situ hybridization (FISH) allowed the specific detection of several of these clusters, whose coexistence in the WWTPs was shown for prolonged periods of several years. In situ ecophysiological analyses based on FISH, relative abundance and spatial arrangement quantification, as well as microautoradiography revealed functional differences of these Nitrospira clusters regarding the preferred nitrite concentration, the utilization of formate as substrate and the spatial coaggregation with ammonia-oxidizing bacteria as symbiotic partners. Amplicon pyrosequencing of the nxrB gene, which encodes subunit beta of nitrite oxidoreductase of Nitrospira, revealed in one of the WWTPs as many as 121 species-level nxrB operational taxonomic units with highly uneven relative abundances in the amplicon library. These results show a previously unrecognized high diversity of Nitrospira in engineered systems, which is at least partially linked to niche differentiation and may have important implications for process stability.

Selective whole genome amplification for resequencing target microbial species from complex natural samples

Population genomic analyses have demonstrated power to address major questions in evolutionary and molecular microbiology. Collecting populations of genomes is hindered in many microbial species by the absence of a cost effective and practical method to collect ample quantities of sufficiently pure genomic DNA for next-generation sequencing. Here we present a simple method to amplify genomes of a target microbial species present in a complex, natural sample. The selective whole genome amplification (SWGA) technique amplifies target genomes using nucleotide sequence motifs that are common in the target microbe genome, but rare in the background genomes, to prime the highly processive phi29 polymerase. SWGA thus selectively amplifies the target genome from samples in which it originally represented a minor fraction of the total DNA. The post-SWGA samples are enriched in target genomic DNA, which are ideal for population resequencing. We demonstrate the efficacy of SWGA using both laboratory-prepared mixtures of cultured microbes as well as a natural host-microbe association. Targeted amplification of Borrelia burgdorferi mixed with Escherichia coli at genome ratios of 1:2000 resulted in >10(5)-fold amplification of the target genomes with <6.7-fold amplification of the background. SWGA-treated genomic extracts from Wolbachia pipientis-infected Drosophila melanogaster resulted in up to 70% of high-throughput resequencing reads mapping to the W. pipientis genome. By contrast, 2-9% of sequencing reads were derived from W. pipientis without prior amplification. The SWGA technique results in high sequencing coverage at a fraction of the sequencing effort, thus allowing population genomic studies at affordable costs.

Inter-species interconnections in acid mine drainage microbial communities

Metagenomic studies are revolutionizing our understanding of microbes in the biosphere. They have uncovered numerous proteins of unknown function in tens of essentially unstudied lineages that lack cultivated representatives. Notably, few of these microorganisms have been visualized, and even fewer have been described ultra-structurally in their essentially intact, physiologically relevant states. Here, we present cryogenic transmission electron microscope (cryo-TEM) 2D images and 3D tomographic datasets for archaeal species from natural acid mine drainage (AMD) microbial communities. Ultrastructural findings indicate the importance of microbial interconnectedness via a range of mechanisms, including direct cytoplasmic bridges and pervasive pili. The data also suggest a variety of biological structures associated with cell-cell interfaces that lack explanation. Some may play roles in inter-species interactions. Interdependences amongst the archaea may have confounded prior isolation efforts. Overall, the findings underline knowledge gaps related to archaeal cell components and highlight the likely importance of co-evolution in shaping microbial lineages.

Elevated temperature alters proteomic responses of individual organisms within a biofilm community

Microbial communities that underpin global biogeochemical cycles will likely be influenced by elevated temperature associated with environmental change. Here, we test an approach to measure how elevated temperature impacts the physiology of individual microbial groups in a community context, using a model microbial-based ecosystem. The study is the first application of tandem mass tag (TMT)-based proteomics to a microbial community. We accurately, precisely and reproducibly quantified thousands of proteins in biofilms growing at 40, 43 and 46 °C. Elevated temperature led to upregulation of proteins involved in amino-acid metabolism at the level of individual organisms and the entire community. Proteins from related organisms differed in their relative abundance and functional responses to temperature. Elevated temperature repressed carbon fixation proteins from two Leptospirillum genotypes, whereas carbon fixation proteins were significantly upregulated at higher temperature by a third member of this genus. Leptospirillum group III bacteria may have been subject to viral stress at elevated temperature, which could lead to greater carbon turnover in the microbial food web through the release of viral lysate. Overall, these findings highlight the utility of proteomics-enabled community-based physiology studies, and provide a methodological framework for possible extension to additional mixed culture and environmental sample analyses.

Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus

Extensive genomic diversity within coexisting members of a microbial species has been revealed through selected cultured isolates and metagenomic assemblies. Yet, the cell-by-cell genomic composition of wild uncultured populations of co-occurring cells is largely unknown. In this work, we applied large-scale single-cell genomics to study populations of the globally abundant marine cyanobacterium Prochlorococcus. We show that they are composed of hundreds of subpopulations with distinct "genomic backbones," each backbone consisting of a different set of core gene alleles linked to a small distinctive set of flexible genes. These subpopulations are estimated to have diverged at least a few million years ago, suggesting ancient, stable niche partitioning. Such a large set of coexisting subpopulations may be a general feature of free-living bacterial species with huge populations in highly mixed habitats.

Ordering microbial diversity into ecologically and genetically cohesive units

We propose that microbial diversity must be viewed in light of gene flow and selection, which define units of genetic similarity, and of phenotype and ecological function, respectively. We discuss to what extent ecological and genetic units overlap to form cohesive populations in the wild, based on recent evolutionary modeling and on evidence from some of the first microbial populations studied with genomics. These show that if recombination is frequent and selection moderate, ecologically adaptive mutations or genes can spread within populations independently of their original genomic background (gene-specific sweeps). Alternatively, if the effect of recombination is smaller than selection, genome-wide selective sweeps should occur. In both cases, however, distinct units of overlapping ecological and genotypic similarity will form if microgeographic separation, likely involving ecological tradeoffs, induces barriers to gene flow. These predictions are supported by (meta)genomic data, which suggest that a 'reverse ecology' approach, in which genomic and gene flow information is used to make predictions about the nature of ecological units, is a powerful approach to ordering microbial diversity.