Genomic and Metabolic Diversity of Marine Group I Thaumarchaeota in the Mesopelagic of Two Subtropical Gyres

Marine Group I (MGI) Thaumarchaeota are one of the most abundant and cosmopolitan chemoautotrophs within the global dark ocean. To date, no representatives of this archaeal group retrieved from the dark ocean have been successfully cultured. We used single cell genomics to investigate the genomic and metabolic diversity of thaumarchaea within the mesopelagic of the subtropical North Pacific and South Atlantic Ocean. Phylogenetic and metagenomic recruitment analysis revealed that MGI single amplified genomes (SAGs) are genetically and biogeographically distinct from existing thaumarchaea cultures obtained from surface waters. Confirming prior studies, we found genes encoding proteins for aerobic ammonia oxidation and the hydrolysis of urea, which may be used for energy production, as well as genes involved in 3-hydroxypropionate/4-hydroxybutyrate and oxidative tricarboxylic acid pathways. A large proportion of protein sequences identified in MGI SAGs were absent in the marine cultures Cenarchaeum symbiosum and Nitrosopumilus maritimus, thus expanding the predicted protein space for this archaeal group. Identifiable genes located on genomic islands with low metagenome recruitment capacity were enriched in cellular defense functions, likely in response to viral infections or grazing. We show that MGI Thaumarchaeota in the dark ocean may have more flexibility in potential energy sources and adaptations to biotic interactions than the existing, surface-ocean cultures.     

The genome of the ammonia‐oxidizing Candidatus Nitrososphaera gargensis: insights into metabolic versatility and environmental adaptations

The cohort of the ammonia‐oxidizing archaea (AOA) of the phylum Thaumarchaeota is a diverse, widespread and functionally important group of microorganisms in many ecosystems. However, our understanding of their biology is still very rudimentary in part because all available genome sequences of this phylum are from members of the Nitrosopumilus cluster. Here we report on the complete genome sequence of Candidatus Nitrososphaera gargensis obtained from an enrichment culture, representing a different evolutionary lineage of AOA frequently found in high numbers in many terrestrial environments. With its 2.83 Mb the genome is much larger than that of other AOA. The presence of a high number of (active) IS elements/transposases, genomic islands, gene duplications and a complete CRISPR/Cas defence system testifies to its dynamic evolution consistent with low degree of synteny with other thaumarchaeal genomes. As expected, the repertoire of conserved enzymes proposed to be required for archaeal ammonia oxidation is encoded by N. gargensis, but it can also use urea and possibly cyanate as alternative ammonia sources. Furthermore, its carbon metabolism is more flexible at the central pyruvate switch point, encompasses the ability to take up small organic compounds and might even include an oxidative pentose phosphate pathway. Furthermore, we show that thaumarchaeota produce cofactor F420 as well as polyhydroxyalkanoates. Lateral gene transfer from bacteria and euryarchaeota has contributed to the metabolic versatility of N. gargensis. This organisms is well adapted to its niche in a heavy metal‐containing thermal spring by encoding a multitude of heavy metal resistance genes, chaperones and mannosylglycerate as compatible solute and has the genetic ability to respond to environmental changes by signal transduction via a large number of two‐component systems, by chemotaxis and flagella‐mediated motility and possibly even by gas vacuole formation. These findings extend our understanding of thaumarchaeal evolution and physiology and offer many testable hypotheses for future experimental research on these nitrifiers.     

Southern Appalachian Peatlands Support High Archaeal Diversity

Mid-latitude peatlands with a temperate climate are sparsely studied and as such represent a gap in the current knowledge base regarding archaeal populations present and their roles in these environments. Phylogenetic analysis of the archaeal populations among three peatlands in the Southern Appalachians reveal not only methanogenic species but also significant populations of thaumarchaeal and crenarchaeal-related organisms of the uncultured miscellaneous crenarchaeotal group (MCG) and the terrestrial group 1.1c, as well as deep-branching Euryarchaeota primarily within the Lake Dagow sediment and rice cluster V lineages. The Thaum/Crenarchaea and deep-branching Euryarchaea represented approximately 24–83 % and 2–18 %, respectively, of the total SSU rRNA clones retrieved in each library, and methanogens represented approximately 14–72 % of the clones retrieved. Several taxa that are either rare or novel to acidic peatlands were detected including the euryarchaeal SM1K20 cluster and thaumarchaeal/crenarchaeal-related clusters 1.1a, C3, SAGMCG-1, pSL12, and AK59. All three major groups (methanogens, Thaumarchaea/Crenarchaea, and deep-branching Euryarchaea) were detected in the RNA library, suggesting at least a minimum level of maintenance activity. Compared to their northern counterparts, Southern Appalachian peatlands appear to harbor a relatively high diversity of Archaea and exhibit a high level of intra-site heterogeneity.     

Assembly-Driven Community Genomics of a Hypersaline Microbial Ecosystem

Microbial populations inhabiting a natural hypersaline lake ecosystem in Lake Tyrrell, Victoria, Australia, have been characterized using deep metagenomic sampling, iterative de novo assembly, and multidimensional phylogenetic binning. Composite genomes representing habitat-specific microbial populations were reconstructed for eleven different archaea and one bacterium, comprising between 0.6 and 14.1% of the planktonic community. Eight of the eleven archaeal genomes were from microbial species without previously cultured representatives. These new genomes provide habitat-specific reference sequences enabling detailed, lineage-specific compartmentalization of predicted functional capabilities and cellular properties associated with both dominant and less abundant community members, including organisms previously known only by their 16S rRNA sequences. Together, these data provide a comprehensive, culture-independent genomic blueprint for ecosystem-wide analysis of protein functions, population structure, and lifestyles of co-existing, co-evolving microbial groups within the same natural habitat. The “assembly-driven” community genomic approach demonstrated in this study advances our ability to push beyond single gene investigations, and promotes genome-scale reconstructions as a tangible goal in the quest to define the metabolic, ecological, and evolutionary dynamics that underpin environmental microbial diversity.     

Genomics, environmental genomics and the issue of microbial species

A microbial species concept is crucial for interpreting the variation detected by genomics and environmental genomics among cultivated microorganisms and within natural microbial populations. Comparative genomic analyses of prokaryotic species as they are presently described and named have led to the provocative idea that prokaryotes may not form species as we think about them for plants and animals. There are good reasons to doubt whether presently recognized prokaryotic species are truly species. To achieve a better understanding of microbial species, we believe it is necessary to (i) re-evaluate traditional approaches in light of evolutionary and ecological theory, (ii) consider that different microbial species may have evolved in different ways and (iii) integrate genomic, metagenomic and genome-wide expression approaches with ecological and evolutionary theory. Here, we outline how we are using genomic methods to (i) identify ecologically distinct populations (ecotypes) predicted by theory to be species-like fundamental units of microbial communities, and (ii) test their species-like character through in situ distribution and gene expression studies. By comparing metagenomic sequences obtained from well-studied hot spring cyanobacterial mats with genomic sequences of two cultivated cyanobacterial ecotypes, closely related to predominant native populations, we can conduct in situ population genetics studies that identify putative ecotypes and functional genes that determine the ecotypes' ecological distinctness. If individuals within microbial communities are found to be grouped into ecologically distinct, species-like populations, knowing about such populations should guide us to a better understanding of how genomic variation is linked to community function.     

The Thaumarchaeota: an emerging view of their phylogeny and ecophysiology

Thaumarchaeota range among the most abundant archaea on Earth. Initially classified as 'mesophilic Crenarchaeota', comparative genomics has recently revealed that they form a separate and deep-branching phylum within the Archaea. This novel phylum comprises in 16S rRNA gene trees not only all known archaeal ammonia oxidizers but also several clusters of environmental sequences representing microorganisms with unknown energy metabolism. Ecophysiological studies of ammonia-oxidizing Thaumarchaeota suggest adaptation to low ammonia concentrations and an autotrophic or possibly mixotrophic lifestyle. Extrapolating from the wide substrate range of copper-containing membrane-bound monooxygenases, to which the thaumarchaeal ammonia monooxygenases belong, the use of substrates other than ammonia for generating energy by some members of the Thaumarchaeota seems likely.     

Genome repository of oil systems: An interactive and searchable database that expands the catalogued diversity of crude oil-associated microbes

Microbial communities ultimately control the fate of petroleum hydrocarbons (PHCs) that enter the natural environment, but the interactions of microbes with PHCs and the environment are highly complex and poorly understood. Genome-resolved metagenomics can help unravel these complex interactions. However, the lack of a comprehensive database that integrates existing genomic/metagenomic data from oil environments with physicochemical parameters known to regulate the fate of PHCs currently limits data analysis and interpretations. Here, we curated a comprehensive, searchable database that documents microbial populations in natural oil ecosystems and oil spills, along with available underlying physicochemical data, geocoded via geographic information system to reveal their geographic distribution patterns. Analysis of the ~2000 metagenome-assembled genomes (MAGs) available in the database revealed strong ecological niche specialization within habitats. Over 95% of the recovered MAGs represented novel taxa underscoring the limited representation of cultured organisms from oil-contaminated and oil reservoir ecosystems. The majority of MAGs linked to oil-contaminated ecosystems were detectable in non-oiled samples from the Gulf of Mexico but not in comparable samples from elsewhere, indicating that the Gulf is primed for oil biodegradation. The repository should facilitate future work toward a predictive understanding of the microbial taxa and their activities that control the fate of oil spills.     

Comparative metagenomics of toxic freshwater cyanobacteria bloom communities on two continents

Toxic cyanobacterial blooms have persisted in freshwater systems around the world for centuries and appear to be globally increasing in frequency and severity. Toxins produced by bloom-associated cyanobacteria can have drastic impacts on the ecosystem and surrounding communities, and bloom biomass can disrupt aquatic food webs and act as a driver for hypoxia. Little is currently known regarding the genomic content of the Microcystis strains that form blooms or the companion heterotrophic community associated with bloom events. To address these issues, we examined the bloom-associated microbial communities in single samples from Lake Erie (North America), Lake Tai (Taihu, China), and Grand Lakes St. Marys (OH, USA) using comparative metagenomics. Together the Cyanobacteria and Proteobacteria comprised >90% of each bloom bacterial community sample, although the dominant phylum varied between systems. Relative to the existing Microcystis aeruginosa NIES 843 genome, sequences from Lake Erie and Taihu revealed a number of metagenomic islands that were absent in the environmental samples. Moreover, despite variation in the phylogenetic assignments of bloom-associated organisms, the functional potential of bloom members remained relatively constant between systems. This pattern was particularly noticeable in the genomic contribution of nitrogen assimilation genes. In Taihu, the genetic elements associated with the assimilation and metabolism of nitrogen were predominantly associated with Proteobacteria, while these functions in the North American lakes were primarily contributed to by the Cyanobacteria. Our observations build on an emerging body of metagenomic surveys describing the functional potential of microbial communities as more highly conserved than that of their phylogenetic makeup within natural systems.     

Stimulation of thaumarchaeal ammonia oxidation by ammonia derived from organic nitrogen but not added inorganic nitrogen

Ammonia oxidation, the first step in nitrification, is performed by autotrophic bacteria and thaumarchaea, whose relative contributions vary in different soils. Distinctive environmental niches for the two groups have not been identified, but evidence from previous studies suggests that activity of thaumarchaea, unlike that of bacterial ammonia oxidizers, is unaffected by addition of inorganic N fertilizer and that they preferentially utilize ammonia generated from the mineralization of organic N. This hypothesis was tested by determining the influence of both inorganic and organic N sources on nitrification rate and ammonia oxidizer growth and community structure in microcosms containing acidic, forest soil in which ammonia oxidation was dominated by thaumarchaea. Nitrification rate was unaffected by the incubation of soil with inorganic ammonium but was significantly stimulated by the addition of organic N. Oxidation of ammonia generated from native soil organic matter or added organic N, but not added inorganic N, was accompanied by increases in abundance of the thaumarchaeal amoA gene, a functional gene for ammonia oxidation, but changes in community structure were not observed. Bacterial amoA genes could not be detected. Ammonia oxidation was completely inhibited by 0.01% acetylene in all treatments, indicating ammonia monooxygenase-dependent activity. The findings have implications for current models of soil nitrification and for nitrification control strategies to minimize fertilizer loss and nitrous oxide production.     

A novel bioinformatics strategy for searching industrially useful genome resources from metagenomic sequence libraries

Although remarkable progress in metagenomic sequencing of various environmental samples has been made, large numbers of fragment sequences have been registered in the international DNA databanks, primarily without information on gene function and phylotype, and thus with limited usefulness. Industrial useful biological activity is often carried out by a set of genes, such as those constituting an operon. In this connection, metagenomic approaches have a weakness because sets of the genes are usually split up, since the sequences obtained by metagenome analyses are fragmented into 1-kb or much shorter segments. Therefore, even when a set of genes responsible for an industrially useful function is found in one metagenome library, it is usually difficult to know whether a single genome harbors the entire gene set or whether different genomes have individual genes. By modifying Self-Organizing Map (SOM), we previously developed BLSOM for oligonucleotide composition, which allowed classification (self-organization) of sequence fragments according to genomes. Because BLSOM could reassociate genomic fragments according to genomes, BLSOM may ameliorate the abovementioned weakness of metagenome analyses. Here, we have developed a strategy for clustering of metagenomic sequences according to phylotypes and genomes, by testing a gene set contributing to environment preservation.     

Ammonia-oxidizing archaea have more important role than ammonia-oxidizing bacteria in ammonia oxidation of strongly acidic soils

Increasing evidence demonstrated the involvement of ammonia-oxidizing archaea (AOA) in the global nitrogen cycle, but the relative contributions of AOA and ammonia-oxidizing bacteria (AOB) to ammonia oxidation are still in debate. Previous studies suggest that AOA would be more adapted to ammonia-limited oligotrophic conditions, which seems to be favored by protonation of ammonia, turning into ammonium in low-pH environments. Here, we investigated the autotrophic nitrification activity of AOA and AOB in five strongly acidic soils (pH<4.50) during microcosm incubation for 30 days. Significantly positive correlations between nitrate concentration and amoA gene abundance of AOA, but not of AOB, were observed during the active nitrification. ¹³CO₂-DNA-stable isotope probing results showed significant assimilation of ¹³C-labeled carbon source into the amoA gene of AOA, but not of AOB, in one of the selected soil samples. High levels of thaumarchaeal amoA gene abundance were observed during the active nitrification, coupled with increasing intensity of two denaturing gradient gel electrophoresis bands for specific thaumarchaeal community. Addition of the nitrification inhibitor dicyandiamide (DCD) completely inhibited the nitrification activity and CO₂ fixation by AOA, accompanied by decreasing thaumarchaeal amoA gene abundance. Bacterial amoA gene abundance decreased in all microcosms irrespective of DCD addition, and mostly showed no correlation with nitrate concentrations. Phylogenetic analysis of thaumarchaeal amoA gene and 16S rRNA gene revealed active ¹³CO₂-labeled AOA belonged to groups 1.1a-associated and 1.1b. Taken together, these results provided strong evidence that AOA have a more important role than AOB in autotrophic ammonia oxidation in strongly acidic soils.     

The Sorcerer II Global Ocean Sampling Expedition: Northwest Atlantic through Eastern Tropical Pacific

The world's oceans contain a complex mixture of micro-organisms that are for the most part, uncharacterized both genetically and biochemically. We report here a metagenomic study of the marine planktonic microbiota in which surface (mostly marine) water samples were analyzed as part of the Sorcerer II Global Ocean Sampling expedition. These samples, collected across a several-thousand km transect from the North Atlantic through the Panama Canal and ending in the South Pacific yielded an extensive dataset consisting of 7.7 million sequencing reads (6.3 billion bp). Though a few major microbial clades dominate the planktonic marine niche, the dataset contains great diversity with 85% of the assembled sequence and 57% of the unassembled data being unique at a 98% sequence identity cutoff. Using the metadata associated with each sample and sequencing library, we developed new comparative genomic and assembly methods. One comparative genomic method, termed “fragment recruitment,” addressed questions of genome structure, evolution, and taxonomic or phylogenetic diversity, as well as the biochemical diversity of genes and gene families. A second method, termed “extreme assembly,” made possible the assembly and reconstruction of large segments of abundant but clearly nonclonal organisms. Within all abundant populations analyzed, we found extensive intra-ribotype diversity in several forms: (1) extensive sequence variation within orthologous regions throughout a given genome; despite coverage of individual ribotypes approaching 500-fold, most individual sequencing reads are unique; (2) numerous changes in gene content some with direct adaptive implications; and (3) hypervariable genomic islands that are too variable to assemble. The intra-ribotype diversity is organized into genetically isolated populations that have overlapping but independent distributions, implying distinct environmental preference. We present novel methods for measuring the genomic similarity between metagenomic samples and show how they may be grouped into several community types. Specific functional adaptations can be identified both within individual ribotypes and across the entire community, including proteorhodopsin spectral tuning and the presence or absence of the phosphate-binding gene PstS.     

The Sorcerer II Global Ocean Sampling Expedition: Northwest Atlantic through Eastern Tropical Pacific

The world's oceans contain a complex mixture of micro-organisms that are for the most part, uncharacterized both genetically and biochemically. We report here a metagenomic study of the marine planktonic microbiota in which surface (mostly marine) water samples were analyzed as part of the Sorcerer II Global Ocean Sampling expedition. These samples, collected across a several-thousand km transect from the North Atlantic through the Panama Canal and ending in the South Pacific yielded an extensive dataset consisting of 7.7 million sequencing reads (6.3 billion bp). Though a few major microbial clades dominate the planktonic marine niche, the dataset contains great diversity with 85% of the assembled sequence and 57% of the unassembled data being unique at a 98% sequence identity cutoff. Using the metadata associated with each sample and sequencing library, we developed new comparative genomic and assembly methods. One comparative genomic method, termed “fragment recruitment,” addressed questions of genome structure, evolution, and taxonomic or phylogenetic diversity, as well as the biochemical diversity of genes and gene families. A second method, termed “extreme assembly,” made possible the assembly and reconstruction of large segments of abundant but clearly nonclonal organisms. Within all abundant populations analyzed, we found extensive intra-ribotype diversity in several forms: (1) extensive sequence variation within orthologous regions throughout a given genome; despite coverage of individual ribotypes approaching 500-fold, most individual sequencing reads are unique; (2) numerous changes in gene content some with direct adaptive implications; and (3) hypervariable genomic islands that are too variable to assemble. The intra-ribotype diversity is organized into genetically isolated populations that have overlapping but independent distributions, implying distinct environmental preference. We present novel methods for measuring the genomic similarity between metagenomic samples and show how they may be grouped into several community types. Specific functional adaptations can be identified both within individual ribotypes and across the entire community, including proteorhodopsin spectral tuning and the presence or absence of the phosphate-binding gene PstS.     

Phylogeography and cryptic introduction of the ragworm Hediste diversicolor (Annelida,Nereididae) in the Northwest Atlantic

Many benthic marine invertebrates show striking range disjunctions across broad spatial scales. Without direct evidence for endemism or introduction, these species remain cryptogenic. The common ragworm Hediste diversicolor plays a pivotal role in sedimentary littoral ecosystems of the North Atlantic as an abundant prey item and ecosystem engineer, but exhibits a restricted dispersal capacity that may limit connectivity at both evolutionary and ecological time scales. In Europe, H. diversicolor is subdivided into cryptic taxa and genetic lineages whose distributions have been modified by recent invasions. Its origin in the northwest Atlantic has not been adequately addressed. To trace the age and origin of North American ragworm populations, we analyzed mtDNA sequence data (COI) from the Gulf of Maine and Bay of Fundy (n=73 individuals) and compared our findings with published data from the northeast Atlantic. Our results together with previous data indicate that two species of the H. diversicolor complex have independently colonized the northwest Atlantic at least three different times, resulting in two distinct conspecific assemblages in the Bay of Fundy and Gulf of Maine (respectively) that are different from the species found in the Gulf of St. Lawrence. North American populations had significantly lower genetic diversity compared with populations in the northeast Atlantic, and based on patterns of shared identity, populations in the Bay of Fundy originated from the Baltic Sea and North Sea. Populations from the Gulf of Maine were phylogenetically distinct and most likely originated from unsampled European populations. Analyses of the North American populations revealed patterns of post‐colonization gene flow among populations within the Gulf of Maine and Bay of Fundy. However, we failed to detect shared haplotypes between the two regions, and this pattern of complete isolation corroborates a strong phylogeographic break observed in other species.     

Single-cell genomic sequencing using Multiple Displacement Amplification

Single microbial cells can now be sequenced using DNA amplified by the Multiple Displacement Amplification (MDA) reaction. The few femtograms of DNA in a bacterium are amplified into micrograms of high molecular weight DNA suitable for DNA library construction and Sanger sequencing. The MDA-generated DNA also performs well when used directly as template for pyrosequencing by the 454 Life Sciences method. While MDA from single cells loses some of the genomic sequence, this approach will greatly accelerate the pace of sequencing from uncultured microbes. The genetically linked sequences from single cells are also a powerful tool to be used in guiding genomic assembly of shotgun sequences of multiple organisms from environmental DNA extracts (metagenomic sequences).     

Genome biology of Actinobacillus pleuropneumoniae JL03, an isolate of serotype 3 prevalent in China

Actinobacillus pleuropneumoniae is the etiologic agent of porcine contagious pleuropneumonia, a cause of considerable world wide economic losses in the swine industry. We sequenced the complete genome of A. pleuropneumoniae, JL03, an isolate of serotype 3 prevalent in China. Its genome is a single chromosome of 2,242,062 base pairs containing 2,097 predicted protein-coding sequences, six ribosomal rRNA operons, and 63 tRNA genes. Preliminary analysis of the genomic sequence and the functions of the encoded proteins not only confirmed the present physiological and pathological knowledge but also offered new insights into the metabolic and virulence characteristics of this important pathogen. We identified a full spectrum of genes related to its characteristic chemoheterotrophic catabolism of fermentation and respiration with an incomplete TCA system for anabolism. In addition to confirming the lack of ApxI toxin, identification of a nonsense mutation in apxIVA and a 5'-proximal truncation of the flp operon deleting both its promoter and the flp1flp2tadV genes have provided convincing scenarios for the low virulence property of JL03. Comparative genomic analysis using the available sequences of other serotypes, probable strain (serotype)-specific genomic islands related to capsular polysaccharides and lipopolysaccharide O-antigen biosyntheses were identified in JL03, which provides a foundation for future research into the mechanisms of serotypic diversity of A. pleuropneumoniae.     

Comparative genomics of microbial drug efflux systems

The complete genome sequences of 36 microorganisms have now been published and this wealth of genome data has enabled the development of comparative genomic and functional genomic approaches to investigate the biology of these organisms. Comparative genomic analyses of membrane transport systems have revealed that transporter substrate specificities correlate with an organism's lifestyle. The types and numbers of predicted drug efflux systems vary dramatically amongst sequenced organisms. Microarray and gene knockout studies to date have suggested that predicted drug efflux genes often appear to be a) non-essential and b) expressed at detectable levels under standard laboratory growth conditions.