A communal catalogue reveals Earth's multiscale microbial diversity.

Our growing awareness of the microbial world's importance and diversity contrasts starkly with our limited understanding of its fundamental structure. Despite recent advances in DNA sequencing, a lack of standardized protocols and common analytical frameworks impedes comparisons among studies, hindering the development of global inferences about microbial life on Earth. Here we present a meta-analysis of microbial community samples collected by hundreds of researchers for the Earth Microbiome Project. Coordinated protocols and new analytical methods, particularly the use of exact sequences instead of clustered operational taxonomic units, enable bacterial and archaeal ribosomal RNA gene sequences to be followed across multiple studies and allow us to explore patterns of diversity at an unprecedented scale. The result is both a reference database giving global context to DNA sequence data and a framework for incorporating data from future studies, fostering increasingly complete characterization of Earth's microbial diversity.

Red Sea SAR11 and Prochlorococcus Single-Cell Genomes Reflect Globally Distributed Pangenomes.

Evidence suggests many marine bacteria are cosmopolitan, with widespread but sparse strains poised to seed abundant populations under conducive growth conditions. However, studies supporting this "microbial seed bank" hypothesis have analyzed taxonomic marker genes rather than whole genomes/metagenomes, leaving open the possibility that disparate ocean regions harbor endemic gene content. The Red Sea is isolated geographically from the rest of the ocean and has a combination of high irradiance, high temperature, and high salinity that is unique among the oceans; we therefore asked whether it harbors endemic gene content. We sequenced and assembled single-cell genomes of 21 SAR11 (subclades Ia, Ib, Id, and II) and 5 Prochlorococcus (ecotype HLII) samples from the Red Sea and combined them with globally sourced reference genomes to cluster genes into ortholog groups (OGs). Ordination of OG composition could distinguish clades, including phylogenetically cryptic Prochlorococcus ecotypes LLII and LLIII. Compared with reference genomes, 1% of Prochlorococcus and 17% of SAR11 OGs were unique to the Red Sea genomes (RS-OGs). Most (83%) RS-OGs had no annotated function, but 65% of RS-OGs were expressed in diel Red Sea metatranscriptomes, suggesting they are functional. Searching Tara Oceans metagenomes, RS-OGs were as likely to be found as non-RS-OGs; nevertheless, Red Sea and other warm samples could be distinguished from cooler samples using the relative abundances of OGs. The results suggest that the prevalence of OGs in these surface ocean bacteria is largely cosmopolitan, with differences in population metagenomes manifested by differences in relative abundance rather than complete presence/absence of OGs.IMPORTANCE Studies have shown that as we sequence seawater from a selected environment deeper and deeper, we approach finding every bacterial taxon known for the ocean as a whole. However, such studies have focused on taxonomic marker genes rather than on whole genomes, raising the possibility that the lack of endemism results from the method of investigation. We took a geographically isolated water body, the Red Sea, and sequenced single cells from it. We compared those single-cell genomes to available genomes from around the ocean and to ocean-spanning metagenomes. We showed that gene ortholog groups found in Red Sea genomes but not in other genomes are nevertheless common across global ocean metagenomes. These results suggest that Baas Becking's hypothesis "everything is everywhere, but the environment selects" also applies to gene ortholog groups. This widely dispersed functional diversity may give oceanic microbial communities the functional capacity to respond rapidly to changing conditions.

Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage.

Anaerobic oxidation of methane (AOM) is critical for controlling the flux of methane from anoxic environments. AOM coupled to iron, manganese and sulphate reduction have been demonstrated in consortia containing anaerobic methanotrophic (ANME) archaea. More recently it has been shown that the bacterium Candidatus 'Methylomirabilis oxyfera' can couple AOM to nitrite reduction through an intra-aerobic methane oxidation pathway. Bioreactors capable of AOM coupled to denitrification have resulted in the enrichment of 'M. oxyfera' and a novel ANME lineage, ANME-2d. However, as 'M. oxyfera' can independently couple AOM to denitrification, the role of ANME-2d in the process is unresolved. Here, a bioreactor fed with nitrate, ammonium and methane was dominated by a single ANME-2d population performing nitrate-driven AOM. Metagenomic, single-cell genomic and metatranscriptomic analyses combined with bioreactor performance and (13)C- and (15)N-labelling experiments show that ANME-2d is capable of independent AOM through reverse methanogenesis using nitrate as the terminal electron acceptor. Comparative analyses reveal that the genes for nitrate reduction were transferred laterally from a bacterial donor, suggesting selection for this novel process within ANME-2d. Nitrite produced by ANME-2d is reduced to dinitrogen gas through a syntrophic relationship with an anaerobic ammonium-oxidizing bacterium, effectively outcompeting 'M. oxyfera' in the system. We propose the name Candidatus 'Methanoperedens nitroreducens' for the ANME-2d population and the family Candidatus 'Methanoperedenaceae' for the ANME-2d lineage. We predict that 'M. nitroreducens' and other members of the 'Methanoperedenaceae' have an important role in linking the global carbon and nitrogen cycles in anoxic environments.

Single-cell genomics reveals pyrrolysine-encoding potential in members of uncultivated archaeal candidate division MSBL1.

Pyrrolysine (Pyl), the 22nd canonical amino acid, is only decoded and synthesized by a limited number of organisms in the domains Archaea and Bacteria. Pyl is encoded by the amber codon UAG, typically a stop codon. To date, all known Pyl-decoding archaea are able to carry out methylotrophic methanogenesis. The functionality of methylamine methyltransferases, an important component of corrinoid-dependent methyltransfer reactions, depends on the presence of Pyl. Here, we present a putative pyl gene cluster obtained from single-cell genomes of the archaeal Mediterranean Sea Brine Lakes group 1 (MSBL1) from the Red Sea. Functional annotation of the MSBL1 single cell amplified genomes (SAGs) also revealed a complete corrinoid-dependent methyl-transfer pathway suggesting that members of MSBL1 may possibly be capable of synthesizing Pyl and metabolizing methylated amines.

Metagenomic covariation along densely sampled environmental gradients in the Red Sea.

Oceanic microbial diversity covaries with physicochemical parameters. Temperature, for example, explains approximately half of global variation in surface taxonomic abundance. It is unknown, however, whether covariation patterns hold over narrower parameter gradients and spatial scales, and extending to mesopelagic depths. We collected and sequenced 45 epipelagic and mesopelagic microbial metagenomes on a meridional transect through the eastern Red Sea. We asked which environmental parameters explain the most variation in relative abundances of taxonomic groups, gene ortholog groups, and pathways-at a spatial scale of <2000 km, along narrow but well-defined latitudinal and depth-dependent gradients. We also asked how microbes are adapted to gradients and extremes in irradiance, temperature, salinity, and nutrients, examining the responses of individual gene ortholog groups to these parameters. Functional and taxonomic metrics were equally well explained (75-79%) by environmental parameters. However, only functional and not taxonomic covariation patterns were conserved when comparing with an intruding water mass with different physicochemical properties. Temperature explained the most variation in each metric, followed by nitrate, chlorophyll, phosphate, and salinity. That nitrate explained more variation than phosphate suggested nitrogen limitation, consistent with low surface N:P ratios. Covariation of gene ortholog groups with environmental parameters revealed patterns of functional adaptation to the challenging Red Sea environment: high irradiance, temperature, salinity, and low nutrients. Nutrient-acquisition gene ortholog groups were anti-correlated with concentrations of their respective nutrient species, recapturing trends previously observed across much larger distances and environmental gradients. This dataset of metagenomic covariation along densely sampled environmental gradients includes online data exploration supplements, serving as a community resource for marine microbial ecology.

Draft Genome Sequence of Uncultured SAR324 Bacterium lautmerah10, Binned from a Red Sea Metagenome.

A draft genome of SAR324 bacterium lautmerah10 was assembled from a metagenome of a surface water sample from the Red Sea, Saudi Arabia. The genome is more complete and has a higher G+C content than that of previously sequenced SAR324 representatives. Its genomic information shows a versatile metabolism that confers an advantage to SAR324, which is reflected in its distribution throughout different depths of the marine water column.

A catalogue of 136 microbial draft genomes from Red Sea metagenomes.

Earth is expected to continue warming and the Red Sea is a model environment for understanding the effects of global warming on ocean microbiomes due to its unusually high temperature, salinity and solar irradiance. However, most microbial diversity analyses of the Red Sea have been limited to cultured representatives and single marker gene analyses, hence neglecting the substantial uncultured majority. Here, we report 136 microbial genomes (completion minus contamination is ≥50%) assembled from 45 metagenomes from eight stations spanning the Red Sea and taken from multiple depths between 10 to 500 m. Phylogenomic analysis showed that most of the retrieved genomes belong to seven different phyla of known marine microbes, but more than half representing currently uncultured species. The open-access data presented here is the largest number of Red Sea representative microbial genomes reported in a single study and will help facilitate future studies in understanding the physiology of these microorganisms and how they have adapted to the relatively harsh conditions of the Red Sea.

Phylogenomic analysis of Candidatus 'Izimaplasma' species: free-living representatives from a Tenericutes clade found in methane seeps.

Tenericutes are a unique class of bacteria that lack a cell wall and are typically parasites or commensals of eukaryotic hosts. Environmental 16S rDNA surveys have identified a number of tenericute clades in diverse environments, introducing the possibility that these Tenericutes may represent non-host-associated, free-living microorganisms. Metagenomic sequencing of deep-sea methane seep sediments resulted in the assembly of two genomes from a Tenericutes-affiliated clade currently known as 'NB1-n' (SILVA taxonomy) or 'RF3' (Greengenes taxonomy). Metabolic reconstruction revealed that, like cultured members of the Mollicutes, these 'NB1-n' representatives lack a tricarboxylic acid cycle and instead use anaerobic fermentation of simple sugars for substrate level phosphorylation. Notably, the genomes also contained a number of unique metabolic features including hydrogenases and a simplified electron transport chain containing an RNF complex, cytochrome bd oxidase and complex I. On the basis of the metabolic potential predicted from the annotated genomes, we devised an anaerobic enrichment media that stimulated the growth of these Tenericutes at 10 °C, resulting in a mixed culture where these organisms represented ~60% of the total cells by targeted fluorescence in situ hybridization (FISH). Visual identification by FISH confirmed these organisms were not directly associated with Eukaryotes and electron cryomicroscopy of cells in the enrichment culture confirmed an ultrastructure consistent with the defining phenotypic property of Tenericutes, with a single membrane and no cell wall. On the basis of their unique gene content, phylogenetic placement and ultrastructure, we propose these organisms represent a novel class within the Tenericutes, and suggest the names Candidatus 'Izimaplasma sp. HR1' and Candidatus 'Izimaplasma sp. HR2' for the two genome representatives.

Transcriptomics in the tropics: Total RNA-based profiling of Costa Rican bromeliad-associated communities.

RNA-Seq was used to examine the microbial, eukaryotic, and viral communities in water catchments ('tanks') formed by tropical bromeliads from Costa Rica. In total, transcripts with taxonomic affiliation to a wide array of bacteria, archaea, and eukaryotes, were observed, as well as RNA-viruses that appeared related to the specific presence of eukaryotes. Bacteria from 25 phyla appeared to comprise the majority of transcripts in one tank (Wg24), compared to only 14 phyla in the other (Wg25). Conversely, eukaryotes from only 16 classes comprised the majority of transcripts in Wg24, compared to 24 classes in the Wg25, revealing a greater eukaryote diversity in the latter. Given that these bromeliads had tanks of similar size (i.e. vertical oxygen gradient), and were neighboring with presumed similar light regime and acquisition of leaf litter through-fall, it is possible that pH was the factor governing these differences in bacterial and eukaryotic communities (Wg24 had a tank pH of 3.6 and Wg25 had a tank pH of 6.2). Archaeal diversity was similar in both tanks, represented by 7 orders, with the exception of Methanocellales transcripts uniquely recovered from Wg25. Based on measures of FPKG (fragments mapped per kilobase of gene length), genes involved in methanogenesis, in addition to a spirochaete flagellin gene, were among those most highly expressed in Wg25. Conversely, aldehyde dehydrogenase and monosaccharide-binding protein were among genes most highly expressed in Wg24. The ability to observe specific presence of insect, plant, and fungi-associated RNA-viruses was unexpected. As with other techniques, there are inherent biases in the use of RNA-Seq, however, these data suggest the possibility of understanding the entire community, including ecological interactions, via simultaneous analysis of microbial, eukaryotic, and viral transcripts.

A laboratory investigation of interactions between denitrifying anaerobic methane oxidation (DAMO) and anammox processes in anoxic environments.

This study investigates interactions between recently identified denitrifying anaerobic methane oxidation (DAMO) and anaerobic ammonium oxidation (anammox) processes in controlled anoxic laboratory reactors. Two reactors were seeded with the same inocula containing DAMO organisms Candidatus Methanoperedens nitroreducens and Candidatus Methylomirabilis oxyfera, and anammox organism Candidatus Kuenenia stuttgartiensis. Both were fed with ammonium and methane, but one was also fed with nitrate and the other with nitrite, providing anoxic environments with different electron acceptors. After steady state reached in several months, the DAMO process became solely/primarily responsible for nitrate reduction while the anammox process became solely responsible for nitrite reduction in both reactors. 16S rRNA gene amplicon sequencing showed that the nitrate-driven DAMO organism M. nitroreducens dominated both the nitrate-fed (~70%) and the nitrite-fed (~26%) reactors, while the nitrite-driven DAMO organism M. oxyfera disappeared in both communities. The elimination of M. oxyfera from both reactors was likely the results of this organism being outcompeted by anammox bacteria for nitrite. K. stuttgartiensis was detected at relatively low levels (1-3%) in both reactors.

In-solution fluorescence in situ hybridization and fluorescence-activated cell sorting for single cell and population genome recovery.

Over the past decade, technological advances in whole genome amplification, microfluidics, flow sorting, and high-throughput sequencing have led to the development of single-cell genomics. Single-cell genomic approaches are typically applied to anonymous microbial cells with only morphology providing clues to their identity. However, targeted separation of microorganisms based on phylogenetic markers, such as the 16S rRNA gene, is beginning to emerge in the single-cell genomics field. Here, we describe an in-solution fluorescence in situ hybridization (FISH) protocol which can be combined with fluorescence-activated cell sorting (FACS) for separation of single cells or populations of interest from environmental samples. Sequencing of DNA obtained from sorted cells can be used for the recovery of draft quality genomes, and when performed in parallel with deep metagenomics, can be used to validate and further scaffold metagenomic assemblies. We illustrate in this chapter the feasibility of this FISH-FACS approach by describing the targeted recovery of a novel anaerobic methanotrophic archaeon.

Fixation-free fluorescence in situ hybridization for targeted enrichment of microbial populations.

We modified the standard ribosomal RNA-targeted fluorescence in situ hybridization (FISH) protocol by removing the fixation steps to allow recovery of unmodified nucleic acids. Using this method, hybridized cells could be visualized in two bioreactor sludges and termite hindgut samples by epifluorescence microscopy. We then targeted one bacterial and one archaeal population in the sludge samples with group-specific oligonucleotide probes using in-solution fixation-free FISH and sorted hybridized populations using fluorescence-activated cell sorting (FACS). We could show that sorted populations were highly enriched for the target organisms based on 16S rRNA gene sequencing, thus confirming probe specificity using the modified FISH protocol. This approach should facilitate subsequent genomic sequencing and analysis of targeted populations as DNA is not compromised by crosslinking during fixation.