Community-wide analysis of microbial genome sequence signatures

Comparative genomic analysis of the 'pseudofungus' Hyphochytrium catenoides

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete ‘pseudofungus’; Hyphochytrium catenoides. Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.

Virus Genomes from Deep Sea Sediments Expand the Ocean Megavirome and Support Independent Origins of Viral Gigantism

Genomics and evolution of giant viruses are two of the most vigorously developing areas of virus research. Lately, metagenomics has become the main source of new virus genomes. Here we describe a metagenomic analysis of the genomes of large and giant viruses from deep sea sediments. The assembled new virus genomes substantially expand the known diversity of the nucleocytoplasmic large DNA viruses of eukaryotes. The results support the concept of independent evolution of giant viruses from smaller ancestors in different virus branches.

The nucleocytoplasmic large DNA viruses (NCLDV) of eukaryotes (proposed order, “Megavirales”) include the families Poxviridae, Asfarviridae, Iridoviridae, Ascoviridae, Phycodnaviridae, Marseilleviridae, and Mimiviridae, as well as still unclassified pithoviruses, pandoraviruses, molliviruses, and faustoviruses. Several of these virus groups include giant viruses, with genome and particle sizes exceeding those of many bacterial and archaeal cells. We explored the diversity of the NCLDV in deep sea sediments from the Loki’s Castle hydrothermal vent area. Using metagenomics, we reconstructed 23 high-quality genomic bins of novel NCLDV, 15 of which are related to pithoviruses, 5 to marseilleviruses, 1 to iridoviruses, and 2 to klosneuviruses. Some of the identified pithovirus-like and marseillevirus-like genomes belong to deep branches in the phylogenetic tree of core NCLDV genes, substantially expanding the diversity and phylogenetic depth of the respective groups. The discovered viruses, including putative giant members of the family Marseilleviridae, have a broad range of apparent genome sizes, in agreement with the multiple, independent origins of gigantism in different branches of the NCLDV. Phylogenomic analysis reaffirms the monophyly of the pithovirus-iridovirus-marseillevirus branch of the NCLDV. Similarly to other giant viruses, the pithovirus-like viruses from Loki’s Castle encode translation systems components. Phylogenetic analysis of these genes indicates a greater bacterial contribution than had been detected previously. Genome comparison suggests extensive gene exchange between members of the pithovirus-like viruses and Mimiviridae. Further exploration of the genomic diversity of Megavirales in additional sediment samples is expected to yield new insights into the evolution of giant viruses and the composition of the ocean megavirome.

Co-occurring genomic capacity for anaerobic methane and dissimilatory sulfur metabolisms discovered in the Korarchaeota

Phylogenetic and geological evidence supports the hypothesis that life on Earth originated in thermal environments and conserved energy through methanogenesis or sulfur reduction. Here we describe two populations of the deeply rooted archaeal phylum Korarchaeota, which were retrieved from the metagenome of a circumneutral, suboxic hot spring that contains high levels of sulfate, sulfide, methane, hydrogen and carbon dioxide. One population is closely related to 'Candidatus Korarchaeum cryptofilum OPF8', while the more abundant korarchaeote, 'Candidatus Methanodesulfokores washburnensis', contains genes that are necessary for anaerobic methane and dissimilatory sulfur metabolisms. Phylogenetic and ancestral reconstruction analyses suggest that methane metabolism originated in the Korarchaeota, whereas genes for dissimilatory sulfite reduction were horizontally transferred to the Korarchaeota from the Firmicutes. Interactions among enzymes involved in both metabolisms could facilitate exergonic, sulfite-dependent, anaerobic oxidation of methane to methanol; alternatively, 'Ca. M. washburnensis' could conduct methanogenesis and sulfur reduction independently. Metabolic reconstruction suggests that 'Ca. M. washburnensis' is a mixotroph, capable of amino acid uptake, assimilation of methane-derived carbon and/or CO₂ fixation by archaeal type III-b RuBisCO for scavenging ribose carbon. Our findings link anaerobic methane metabolism and dissimilatory sulfur reduction within a single deeply rooted archaeal population and have implications for the evolution of these traits throughout the Archaea.

Expanding anaerobic alkane metabolism in the domain of Archaea

Methanogenesis and anaerobic methane oxidation through methyl-coenzyme M reductase (MCR) as a key enzyme have been suggested to be basal pathways of archaea¹. How widespread MCR-based alkane metabolism is among archaea, where it occurs and how it evolved remain elusive. Here, we performed a global survey of MCR-encoding genomes based on metagenomic data from various environments. Eleven high-quality mcr-containing metagenomic-assembled genomes were obtained belonging to the Archaeoglobi in the Euryarchaeota, Hadesarchaeota and different TACK superphylum archaea, including the Nezhaarchaeota, Korarchaeota and Verstraetearchaeota. Archaeoglobi WYZ-LMO1 and WYZ-LMO3 and Korarchaeota WYZ-LMO9 encode both the (reverse) methanogenesis and the dissimilatory sulfate reduction pathway, suggesting that they have the genomic potential to couple both pathways in individual organisms. The Hadesarchaeota WYZ-LMO4-6 and Archaeoglobi JdFR-42 encode highly divergent MCRs, enzymes that may enable them to thrive on non-methane alkanes. The occurrence of mcr genes in different archaeal phyla indicates that MCR-based alkane metabolism is common in the domain of Archaea.

Biogeochemical Regimes in Shallow Aquifers Reflect the Metabolic Coupling of the Elements Nitrogen, Sulfur, and Carbon

Near-surface groundwaters are prone to receive (in)organic matter input from their recharge areas and are known to harbor autotrophic microbial communities linked to nitrogen and sulfur metabolism. Here, we use multi-omic profiling to gain holistic insights into the turnover of inorganic nitrogen compounds, carbon fixation processes, and organic matter processing in groundwater. We sampled microbial biomass from two superimposed aquifers via monitoring wells that follow groundwater flow from its recharge area through differences in hydrogeochemical settings and land use. Functional profiling revealed that groundwater microbiomes are mainly driven by nitrogen (nitrification, denitrification, and ammonium oxidation [anammox]) and to a lesser extent sulfur cycling (sulfur oxidation and sulfate reduction), depending on local hydrochemical differences. Surprisingly, the differentiation potential of the groundwater microbiome surpasses that of hydrochemistry for individual monitoring wells. Being dominated by a few phyla (Bacteroidetes, Proteobacteria, Planctomycetes, and Thaumarchaeota), the taxonomic profiling of groundwater metagenomes and metatranscriptomes revealed pronounced differences between merely present microbiome members and those actively participating in community gene expression and biogeochemical cycling. Unexpectedly, we observed a constitutive expression of carbohydrate-active enzymes encoded by different microbiome members, along with the groundwater flow path. The turnover of organic carbon apparently complements for lithoautotrophic carbon assimilation pathways mainly used by the groundwater microbiome depending on the availability of oxygen and inorganic electron donors, like ammonium.IMPORTANCE Groundwater is a key resource for drinking water production and irrigation. The interplay between geological setting, hydrochemistry, carbon storage, and groundwater microbiome ecosystem functioning is crucial for our understanding of these important ecosystem services. We targeted the encoded and expressed metabolic potential of groundwater microbiomes along an aquifer transect that diversifies in terms of hydrochemistry and land use. Our results showed that the groundwater microbiome has a higher spatial differentiation potential than does hydrochemistry.

Viral population analysis of the taiga tick, Ixodes persulcatus, by using Batch Learning Self-Organizing Maps and BLAST search

Ticks transmit a wide range of viral, bacterial, and protozoal pathogens, which are often zoonotic. Several novel tick-borne viral pathogens have been reported during the past few years. The aim of this study was to investigate a diversity of tick viral populations, which may contain as-yet unidentified viruses, using a combination of high throughput pyrosequencing and Batch Learning Self-Organizing Map (BLSOM) program, which enables phylogenetic estimation based on the similarity of oligonucleotide frequencies. DNA/cDNA prepared from virus-enriched fractions obtained from Ixodes persulcatus ticks was pyrosequenced. After de novo assembly, contigs were cataloged by the BLSOM program. In total 41 different viral families and order including those previously associated with human and animal diseases such as Bunyavirales, Flaviviridae, and Reoviridae, were detected. Therefore, our strategy is applicable for viral population analysis of other arthropods of medical and veterinary importance, such as mosquitos and lice. The results lead to the contribution to the prediction of emerging tick-borne viral diseases. A sufficient understanding of tick viral populations will also empower to analyze and understand tick biology including vector competency and interactions with other pathogens.

Reconstruction of the Genomes of Drug-Resistant Pathogens for Outbreak Investigation through Metagenomic Sequencing

The study results reported here perfectly demonstrate the power and promise of clinical metagenomics to recover genome sequences of important drug-resistant bacteria and to rapidly provide rich data that inform outbreak investigations and treatment decisions, independently of the need to culture the organisms.

Culture-independent methods that target genome fragments have shown promise in identifying certain pathogens, but the holy grail of comprehensive pathogen genome detection from microbiologically complex samples for subsequent forensic analyses remains a challenge. In the context of an investigation of a nosocomial outbreak, we used shotgun metagenomic sequencing of a human fecal sample and a neural network algorithm based on tetranucleotide frequency profiling to reconstruct microbial genomes and tested the same approach using rectal swabs from a second patient. The approach rapidly and readily detected the genome of Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae in the patient fecal specimen and in the rectal swab sample, achieving a level of strain resolution that was sufficient for confident transmission inference during a highly clonal outbreak. The analysis also detected previously unrecognized colonization of the patient by vancomycin-resistant Enterococcus faecium, another multidrug-resistant bacterium.

IMPORTANCE The study results reported here perfectly demonstrate the power and promise of clinical metagenomics to recover genome sequences of important drug-resistant bacteria and to rapidly provide rich data that inform outbreak investigations and treatment decisions, independently of the need to culture the organisms.

Genome-scale phylogenetics reveals a monophyletic Zoopagales (Zoopagomycota, Fungi)

Previous genome-scale phylogenetic analyses of Fungi have under sampled taxa from Zoopagales; this order contains many predacious or parasitic genera, and most have never been grown in pure culture. We sequenced the genomes of 4 zoopagalean taxa that are predators of amoebae, nematodes, or rotifers and the genome of one taxon that is a parasite of amoebae using single cell sequencing methods with whole genome amplification. Each genome was a metagenome, which was assembled and binned using multiple techniques to identify the target genomes. We inferred phylogenies with both super matrix and coalescent approaches using 192 conserved proteins mined from the target genomes and performed ancestral state reconstructions to determine the ancestral trophic lifestyle of the clade. Our results indicate that Zoopagales is monophyletic. Ancestral state reconstructions provide moderate support for mycoparasitism being the ancestral state of the clade.

"Candidatus Thermonerobacter thiotrophicus," A Non-phototrophic Member of the Bacteroidetes/Chlorobi With Dissimilatory Sulfur Metabolism in Hot Spring Mat Communities

In this study we present evidence for a novel, thermophilic bacterium with dissimilatory sulfur metabolism, tentatively named “Candidatus Thermonerobacter thiotrophicus,” which is affiliated with the Bacteroides/Ignavibacteria/Chlorobi and which we predict to be a sulfate reducer. Dissimilatory sulfate reduction (DSR) is an important and ancient metabolic process for energy conservation with global importance for geochemical sulfur and carbon cycling. Characterized sulfate-reducing microorganisms (SRM) are found in a limited number of bacterial and archaeal phyla. However, based on highly diverse environmental dsrAB sequences, a variety of uncultivated and unidentified SRM must exist. The recent development of high-throughput sequencing methods allows the phylogenetic identification of some of these uncultured SRM. In this study, we identified a novel putative SRM inhabiting hot spring microbial mats that is a member of the OPB56 clade (“Ca. Kapabacteria”) within the Bacteroidetes/Chlorobi superphylum. Partial genomes for this new organism were retrieved from metagenomes from three different hot springs in Yellowstone National Park, United States, and Japan. Supporting the prediction of a sulfate-reducing metabolism for this organism during period of anoxia, diel metatranscriptomic analyses indicate highest relative transcript levels in situ for all DSR-related genes at night. The presence of terminal oxidases, which are transcribed during the day, further suggests that these organisms might also perform aerobic respiration. The relative phylogenetic proximity to the sulfur-oxidizing, chlorophototrophic Chlorobi further raises new questions about the evolution of dissimilatory sulfur metabolism.

Multiple Data Analyses and Statistical Approaches for Analyzing Data from Metagenomic Studies and Clinical Trials

Metagenomics, also known as environmental genomics, is the study of the genomic content of a sample of organisms (microbes) obtained from a common habitat. Metagenomics and other "omics" disciplines have captured the attention of researchers for several decades. The effect of microbes in our body is a relevant concern for health studies. There are plenty of studies using metagenomics which examine microorganisms that inhabit niches in the human body, sometimes causing disease, and are often correlated with multiple treatment conditions. No matter from which environment it comes, the analyses are often aimed at determining either the presence or absence of specific species of interest in a given metagenome or comparing the biological diversity and the functional activity of a wider range of microorganisms within their communities. The importance increases for comparison within different environments such as multiple patients with different conditions, multiple drugs, and multiple time points of same treatment or same patient. Thus, no matter how many hypotheses we have, we need a good understanding of genomics, bioinformatics, and statistics to work together to analyze and interpret these datasets in a meaningful way. This chapter provides an overview of different data analyses and statistical approaches (with example scenarios) to analyze metagenomics samples from different medical projects or clinical trials.

Elucidating functional microorganisms and metabolic mechanisms in a novel engineered ecosystem integrating C, N, P and S biotransformation by metagenomics

Denitrifying sulfur conversion-associated enhanced biological phosphorous removal (DS-EBPR) system is not only a novel wastewater treatment process, but also an ideal model for microbial ecology in a community context. However, it exists the knowledge gap on the roles and interactions of functional microorganisms in the DS-EBPR system for carbon (C), nitrogen (N), phosphorus (P) and sulfur (S) bioconversions. We use genome-resolved metagenomics to build up an ecological model of microbial communities in a lab-scale DS-EBPR system with stable operation for more than 400 days. Our results yield 11 near-complete draft genomes that represent a substantial portion of the microbial community (39.4%). Sulfate-reducing bacteria (SRB) and sulfide-oxidizing bacteria (SOB) promote complex metabolic processes and interactions for C, N, P and S conversions. Bins 1-4 and 10 are considered as new potential polyphosphate-accumulating organisms (PAOs), in which Bins 1-4 can be considered as S-related PAOs (S-PAOs) with no previously cultivated or reported members. Our findings give an insight into a new ecological system with C, N, P and S simultaneous bioconversions and improve the understanding of interactions among SRB, SOB, denitrifiers and PAOs within a community context.

Expansive microbial metabolic versatility and biodiversity in dynamic Guaymas Basin hydrothermal sediments

Microbes in Guaymas Basin (Gulf of California) hydrothermal sediments thrive on hydrocarbons and sulfur and experience steep, fluctuating temperature and chemical gradients. The functional capacities of communities inhabiting this dynamic habitat are largely unknown. Here, we reconstructed 551 genomes from hydrothermally influenced, and nearby cold sediments belonging to 56 phyla (40 uncultured). These genomes comprise 22 unique lineages, including five new candidate phyla. In contrast to findings from cold hydrocarbon seeps, hydrothermal-associated communities are more diverse and archaea dominate over bacteria. Genome-based metabolic inferences provide first insights into the ecological niches of these uncultured microbes, including methane cycling in new Crenarchaeota and alkane utilization in ANME-1. These communities are shaped by a high biodiversity, partitioning among nitrogen and sulfur pathways and redundancy in core carbon-processing pathways. The dynamic sediments select for distinctive microbial communities that stand out by expansive biodiversity, and open up new physiological perspectives into hydrothermal ecosystem function.

The diversity and function of microbes inhabiting hydrothermal areas is a topic of active interest in marine microbiology. Here, the authors assemble genomes from Guaymas Basin hydrothermal sediments and describe the metabolic roles of the bacterial community, which includes five new bacterial candidate phyla

Members of the Genus Methylobacter Are Inferred To Account for the Majority of Aerobic Methane Oxidation in Oxic Soils from a Freshwater Wetland

Microbial carbon degradation and methanogenesis in wetland soils generate a large proportion of atmospheric methane, a highly potent greenhouse gas. Despite their potential to mitigate greenhouse gas emissions, knowledge about methane-consuming methanotrophs is often limited to lower-resolution single-gene surveys that fail to capture the taxonomic and metabolic diversity of these microorganisms in soils. Here our objective was to use genome-enabled approaches to investigate methanotroph membership, distribution, and in situ activity across spatial and seasonal gradients in a freshwater wetland near Lake Erie. 16S rRNA gene analyses demonstrated that members of the methanotrophic Methylococcales were dominant, with the dominance largely driven by the relative abundance of four taxa, and enriched in oxic surface soils. Three methanotroph genomes from assembled soil metagenomes were assigned to the genus Methylobacter and represented the most abundant methanotrophs across the wetland. Paired metatranscriptomes confirmed that these Old Woman Creek (OWC) Methylobacter members accounted for nearly all the aerobic methanotrophic activity across two seasons. In addition to having the capacity to couple methane oxidation to aerobic respiration, these new genomes encoded denitrification potential that may sustain energy generation in soils with lower dissolved oxygen concentrations. We further show that Methylobacter members that were closely related to the OWC members were present in many other high-methane-emitting freshwater and soil sites, suggesting that this lineage could participate in methane consumption in analogous ecosystems. This work contributes to the growing body of research suggesting that Methylobacter may represent critical mediators of methane fluxes in freshwater saturated sediments and soils worldwide.IMPORTANCE Here we used soil metagenomics and metatranscriptomics to uncover novel members within the genus Methylobacter We denote these closely related genomes as members of the lineage OWC Methylobacter Despite the incredibly high microbial diversity in soils, here we present findings that unexpectedly showed that methane cycling was primarily mediated by a single genus for both methane production ("Candidatus Methanothrix paradoxum") and methane consumption (OWC Methylobacter). Metatranscriptomic analyses revealed that decreased methanotrophic activity rather than increased methanogenic activity possibly contributed to the greater methane emissions that we had previously observed in summer months, findings important for biogeochemical methane models. Although members of this Methylococcales order have been cultivated for decades, multi-omic approaches continue to illuminate the methanotroph phylogenetic and metabolic diversity harbored in terrestrial and marine ecosystems.

Discovery of several novel, widespread, and ecologically distinct marine Thaumarchaeota viruses that encode amoC nitrification genes

Much of the diversity of prokaryotic viruses has yet to be described. In particular, there are no viral isolates that infect abundant, globally significant marine archaea including the phylum Thaumarchaeota. This phylum oxidizes ammonia, fixes inorganic carbon, and thus contributes to globally significant nitrogen and carbon cycles in the oceans. Metagenomics provides an alternative to culture-dependent means for identifying and characterizing viral diversity. Some viruses carry auxiliary metabolic genes (AMGs) that are acquired via horizontal gene transfer from their host(s), allowing inference of what host a virus infects. Here we present the discovery of 15 new genomically and ecologically distinct Thaumarchaeota virus populations, identified as contigs that encode viral capsid and thaumarchaeal ammonia monooxygenase genes (amoC). These viruses exhibit depth and latitude partitioning and are distributed globally in various marine habitats including pelagic waters, estuarine habitats, and hydrothermal plume water and sediments. We found evidence of viral amoC expression and that viral amoC AMGs sometimes comprise up to half of total amoC DNA copies in cellular fraction metagenomes, highlighting the potential impact of these viruses on N cycling in the oceans. Phylogenetics suggest they are potentially tailed viruses and share a common ancestor with related marine Euryarchaeota viruses. This work significantly expands our view of viruses of globally important marine Thaumarchaeota.

Genomes from uncultivated prokaryotes: a comparison of metagenome-assembled and single-amplified genomes

BACKGROUND

Prokaryotes dominate the biosphere and regulate biogeochemical processes essential to all life. Yet, our knowledge about their biology is for the most part limited to the minority that has been successfully cultured. Molecular techniques now allow for obtaining genome sequences of uncultivated prokaryotic taxa, facilitating in-depth analyses that may ultimately improve our understanding of these key organisms.

RESULTS

We compared results from two culture-independent strategies for recovering bacterial genomes: single-amplified genomes and metagenome-assembled genomes. Single-amplified genomes were obtained from samples collected at an offshore station in the Baltic Sea Proper and compared to previously obtained metagenome-assembled genomes from a time series at the same station. Among 16 single-amplified genomes analyzed, seven were found to match metagenome-assembled genomes, affiliated with a diverse set of taxa. Notably, genome pairs between the two approaches were nearly identical (average 99.51% sequence identity; range 98.77-99.84%) across overlapping regions (30-80% of each genome). Within matching pairs, the single-amplified genomes were consistently smaller and less complete, whereas the genetic functional profiles were maintained. For the metagenome-assembled genomes, only on average 3.6% of the bases were estimated to be missing from the genomes due to wrongly binned contigs.

CONCLUSIONS

The strong agreement between the single-amplified and metagenome-assembled genomes emphasizes that both methods generate accurate genome information from uncultivated bacteria. Importantly, this implies that the research questions and the available resources are allowed to determine the selection of genomics approach for microbiome studies.

Metagenomic analysis with strain-level resolution reveals fine-scale variation in the human pregnancy microbiome

Recent studies suggest that the microbiome has an impact on gestational health and outcome. However, characterization of the pregnancy-associated microbiome has largely relied on 16S rRNA gene amplicon-based surveys. Here, we describe an assembly-driven, metagenomics-based, longitudinal study of the vaginal, gut, and oral microbiomes in 292 samples from 10 subjects sampled every three weeks throughout pregnancy. Nonhuman sequences in the amount of 1.53 Gb were assembled into scaffolds, and functional genes were predicted for gene- and pathway-based analyses. Vaginal assemblies were binned into 97 draft quality genomes. Redundancy analysis (RDA) of microbial community composition at all three body sites revealed gestational age to be a significant source of variation in patterns of gene abundance. In addition, health complications were associated with variation in community functional gene composition in the mouth and gut. The diversity of Lactobacillus iners-dominated communities in the vagina, unlike most other vaginal community types, significantly increased with gestational age. The genomes of co-occurring Gardnerella vaginalis strains with predicted distinct functions were recovered in samples from two subjects. In seven subjects, gut samples contained strains of the same Lactobacillus species that dominated the vaginal community of that same subject and not other Lactobacillus species; however, these within-host strains were divergent. CRISPR spacer analysis suggested shared phage and plasmid populations across body sites and individuals. This work underscores the dynamic behavior of the microbiome during pregnancy and suggests the potential importance of understanding the sources of this behavior for fetal development and gestational outcome.

Metabolic marker gene mining provides insight in global mcrA diversity and, coupled with targeted genome reconstruction, sheds further light on metabolic potential of the Methanomassiliicoccales

Over the past years, metagenomics has revolutionized our view of microbial diversity. Moreover, extracting near-complete genomes from metagenomes has led to the discovery of known metabolic traits in unsuspected lineages. Genome-resolved metagenomics relies on assembly of the sequencing reads and subsequent binning of assembled contigs, which might be hampered by strain heterogeneity or low abundance of a target organism. Here we present a complementary approach, metagenome marker gene mining, and use it to assess the global diversity of archaeal methane metabolism through the mcrA gene. To this end, we have screened 18,465 metagenomes for the presence of reads matching a database representative of all known mcrA proteins and reconstructed gene sequences from the matching reads. We use our mcrA dataset to assess the environmental distribution of the Methanomassiliicoccales and reconstruct and analyze a draft genome belonging to the ‘Lake Pavin cluster’, an uncultivated environmental clade of the Methanomassiliicoccales. Analysis of the ‘Lake Pavin cluster’ draft genome suggests that this organism has a more restricted capacity for hydrogenotrophic methylotrophic methanogenesis than previously studied Methanomassiliicoccales, with only genes for growth on methanol present. However, the presence of the soluble subunits of methyltetrahydromethanopterin:CoM methyltransferase (mtrAH) provide hypothetical pathways for methanol fermentation, and aceticlastic methanogenesis that await experimental verification. Thus, we show that marker gene mining can enhance the discovery power of metagenomics, by identifying novel lineages and aiding selection of targets for in-depth analyses. Marker gene mining is less sensitive to strain heterogeneity and has a lower abundance threshold than genome-resolved metagenomics, as it only requires short contigs and there is no binning step. Additionally, it is computationally cheaper than genome resolved metagenomics, since only a small subset of reads needs to be assembled. It is therefore a suitable approach to extract knowledge from the many publicly available sequencing projects.

Metagenomics reveals niche partitioning within the phototrophic zone of a microbial mat

Hypersaline photosynthetic microbial mats are stratified microbial communities known for their taxonomic and metabolic diversity and strong light-driven day-night environmental gradients. In this study of the upper photosynthetic zone of hypersaline microbial mats of Elkhorn Slough, California (USA), we show how metagenome sequencing can be used to meaningfully assess microbial ecology and genetic partitioning in these complex microbial systems. Mapping of metagenome reads to the dominant Cyanobacteria observed in the system, Coleofasciculus (Microcoleus) chthonoplastes, was used to examine strain variants within these metagenomes. Highly conserved gene subsystems indicated a core genome for the species, and a number of variant genes and subsystems suggested strain level differentiation, especially for nutrient utilization and stress response. Metagenome sequence coverage binning was used to assess ecosystem partitioning of remaining microbes to both reconstruct the model organisms in silico and identify their ecosystem functions as well as to identify novel clades and propose their role in the biogeochemical cycling of mats. Functional gene annotation of these bins (primarily of Proteobacteria, Bacteroidetes, and Cyanobacteria) recapitulated the known biogeochemical functions in microbial mats using a genetic basis, and revealed significant diversity in the Bacteroidetes, presumably in heterotrophic cycling. This analysis also revealed evidence of putative phototrophs within the Gemmatimonadetes and Gammaproteobacteria residing in microbial mats. This study shows that metagenomic analysis can produce insights into the systems biology of microbial ecosystems from a genetic perspective and to suggest further studies of novel microbes.

Complex Evolutionary History of Translation Elongation Factor 2 and Diphthamide Biosynthesis in Archaea and Parabasalids

Diphthamide is a modified histidine residue which is uniquely present in archaeal and eukaryotic elongation factor 2 (EF-2), an essential GTPase responsible for catalyzing the coordinated translocation of tRNA and mRNA through the ribosome. In part due to the role of diphthamide in maintaining translational fidelity, it was previously assumed that diphthamide biosynthesis genes (dph) are conserved across all eukaryotes and archaea. Here, comparative analysis of new and existing genomes reveals that some archaea (i.e., members of the Asgard superphylum, Geoarchaea, and Korarchaeota) and eukaryotes (i.e., parabasalids) lack dph. In addition, while EF-2 was thought to exist as a single copy in archaea, many of these dph-lacking archaeal genomes encode a second EF-2 paralog missing key residues required for diphthamide modification and for normal translocase function, perhaps suggesting functional divergence linked to loss of diphthamide biosynthesis. Interestingly, some Heimdallarchaeota previously suggested to be most closely related to the eukaryotic ancestor maintain dph genes and a single gene encoding canonical EF-2. Our findings reveal that the ability to produce diphthamide, once thought to be a universal feature in archaea and eukaryotes, has been lost multiple times during evolution, and suggest that anticipated compensatory mechanisms evolved independently.

Occurrence, integrity and functionality of AcaML1-like viruses infecting extreme acidophiles of the Acidithiobacillus species complex

General knowledge on the diversity and biology of microbial viruses infecting bacterial hosts from extreme acidic environments lags behind most other econiches. In this study, we analyse the AcaML1 virus occurrence in the taxon, its genetic composition and infective behaviour under standard acidic and SOS-inducing conditions to assess its integrity and functionality. Occurrence analysis in sequenced acidithiobacilli showed that AcaML1-like proviruses are confined to the mesothermophiles Acidithiobacillus caldus and Thermithiobacillus tepidarius. Among A. caldus strains and isolates this provirus had a modest prevalence (30%). Comparative genomic analysis revealed a significant conservation with the T. tepidarius AcaML1-like provirus, excepting the tail genes, and a high conservation of the virus across strains of the A. caldus species. Such conservation extends from the modules architecture to the gene level, suggesting that organization and composition of these viruses are preserved for functional reasons. Accordingly, the AcaML1 proviruses were demonstrated to excise from their host genomes under DNA-damaging conditions triggering the SOS-response and to produce DNA-containing VLPs. Despite this fact, under the conditions evaluated (acidic) the VLPs obtained from A. caldus ATCC 51756 could not produce productive infections of a candidate sensitive strain (#6) nor trigger it lysis.

Biphenyl-Metabolizing Microbial Community and a Functional Operon Revealed in E-Waste-Contaminated Soil

Primitive electronic waste (e-waste) recycling activities release massive amounts of persistent organic pollutants (POPs) and heavy metals into surrounding soils, posing a major threat to the ecosystem and human health. Microbes capable of metabolizing POPs play important roles in POPs remediation in soils, but their phylotypes and functions remain unclear. Polychlorinated biphenyls (PCBs), one of the main pollutants in e-waste contaminated soils, have drawn increasing attention due to their high persistence, toxicity, and bioaccumulation. In the present study, we employed the culture-independent method of DNA stable-isotope probing to identify active biphenyl and PCB degraders in e-waste-contaminated soil. A total of 19 rare operational taxonomic units and three dominant bacterial genera ( Ralstonia, Cupriavidus, and uncultured bacterium DA101) were enriched in the ¹³C heavy DNA fraction, confirming their functions in PCBs metabolism. Additionally, a 13.8 kb bph operon was amplified, containing a bphA gene labeled by ¹³C that was concentrated in the heavy DNA fraction. The tetranucleotide signature characteristics of the bph operon suggest that it originated from Ralstonia. The bph operon may be shared by horizontal gene transfer because it contains a transposon gene and is found in various bacterial species. This study gives us a deeper understanding of PCB-degrading mechanisms and provides a potential resource for the bioremediation of PCBs-contaminated soils.

Investigation of mitochondrial-derived plastome sequences in the Paspalum lineage (Panicoideae; Poaceae)

BACKGROUND

The grass family (Poaceae), ca. 12,075 species, is a focal point of many recent studies that aim to use complete plastomes to reveal and strengthen relationships within the family. The use of Next Generation Sequencing technology has revealed intricate details in many Poaceae plastomes; specifically the trnI - trnL intergenic spacer region. This study investigates this region and the putative mitochondrial inserts within it in complete plastomes of Paspalum and other Poaceae.

RESULTS

Nine newly sequenced plastomes, seven of which contain an insert within the trnI - trnL intergenic spacer, were combined into plastome phylogenomic and divergence date analyses with 52 other species. A robust Paspalum topology was recovered, originating at 10.6 Ma, with the insert arising at 8.7 Ma. The alignment of the insert across Paspalum reveals 21 subregions with pairwise homology in 19. In an analysis of emergent self-organizing maps of tetranucleotide frequencies, the Paspalum insert grouped with mitochondrial DNA.

CONCLUSIONS

A hypothetical ancestral insert, 17,685 bp in size, was found in the trnI - trnL intergenic spacer for the Paspalum lineage. A different insert, 2808 bp, was found in the same region for Paraneurachne muelleri. Seven different intrastrand deletion events were found within the Paspalum lineage, suggesting selective pressures to remove large portions of noncoding DNA. Finally, a tetranucleotide frequency analysis was used to determine that the origin of the insert in the Paspalum lineage is mitochondrial DNA.

Genomic inference of the metabolism and evolution of the archaeal phylum Aigarchaeota

Microbes of the phylum Aigarchaeota are widely distributed in geothermal environments, but their physiological and ecological roles are poorly understood. Here we analyze six Aigarchaeota metagenomic bins from two circumneutral hot springs in Tengchong, China, to reveal that they are either strict or facultative anaerobes, and most are chemolithotrophs that can perform sulfide oxidation. Applying comparative genomics to the Thaumarchaeota and Aigarchaeota, we find that they both originated from thermal habitats, sharing 1154 genes with their common ancestor. Horizontal gene transfer played a crucial role in shaping genetic diversity of Aigarchaeota and led to functional partitioning and ecological divergence among sympatric microbes, as several key functional innovations were endowed by Bacteria, including dissimilatory sulfite reduction and possibly carbon monoxide oxidation. Our study expands our knowledge of the possible ecological roles of the Aigarchaeota and clarifies their evolutionary relationship to their sister lineage Thaumarchaeota.

Host-linked soil viral ecology along a permafrost thaw gradient

Climate change threatens to release abundant carbon that is sequestered at high latitudes, but the constraints on microbial metabolisms that mediate the release of methane and carbon dioxide are poorly understood^1–7. The role of viruses, which are known to affect microbial dynamics, metabolism and biogeochemistry in the oceans^8–10, remains largely unexplored in soil. Here, we aimed to investigate how viruses influence microbial ecology and carbon metabolism in peatland soils along a permafrost thaw gradient in Sweden. We recovered 1,907 viral populations (genomes and large genome fragments) from 197 bulk soil and size-fractionated metagenomes, 58% of which were detected in metatranscriptomes and presumed to be active. In silico predictions linked 35% of the viruses to microbial host populations, highlighting likely viral predators of key carbon-cycling microorganisms, including methanogens and methanotrophs. Lineage-specific virus/host ratios varied, suggesting that viral infection dynamics may differentially impact microbial responses to a changing climate. Virus-encoded glycoside hydrolases, including an endomannanase with confirmed functional activity, indicated that viruses influence complex carbon degradation and that viral abundances were significant predictors of methane dynamics. These findings suggest that viruses may impact ecosystem function in climate-critical, terrestrial habitats and identify multiple potential viral contributions to soil carbon cycling.

The recovery of viral populations from peatland soils across a permafrost thaw gradient provides insights into soil viral diversity, their hosts and the potential impacts on carbon cycling in this environment.

Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy

Microbial communities are critical to ecosystem function. A key objective of metagenomic studies is to analyse organism-specific metabolic pathways and reconstruct community interaction networks. This requires accurate assignment of assembled genome fragments to genomes. Existing binning methods often fail to reconstruct a reasonable number of genomes and report many bins of low quality and completeness. Furthermore, the performance of existing algorithms varies between samples and biotopes. Here, we present a dereplication, aggregation and scoring strategy, DAS Tool, that combines the strengths of a flexible set of established binning algorithms. DAS Tool applied to a constructed community generated more accurate bins than any automated method. Indeed, when applied to environmental and host-associated samples of different complexity, DAS Tool recovered substantially more near-complete genomes, including previously unreported lineages, than any single binning method alone. The ability to reconstruct many near-complete genomes from metagenomics data will greatly advance genome-centric analyses of ecosystems.

Alignment-Free Sequence Analysis and Applications

Genome and metagenome comparisons based on large amounts of next generation sequencing (NGS) data pose significant challenges for alignment-based approaches due to the huge data size and the relatively short length of the reads. Alignment-free approaches based on the counts of word patterns in NGS data do not depend on the complete genome and are generally computationally efficient. Thus, they contribute significantly to genome and metagenome comparison. Recently, novel statistical approaches have been developed for the comparison of both long and shotgun sequences. These approaches have been applied to many problems including the comparison of gene regulatory regions, genome sequences, metagenomes, binning contigs in metagenomic data, identification of virus-host interactions, and detection of horizontal gene transfers. We provide an updated review of these applications and other related developments of word-count based approaches for alignment-free sequence analysis.

A closely-related clade of globally distributed bloom-forming cyanobacteria within the Nostocales

In order to better understand the relationships among current Nostocales cyanobacterial blooms, eight genomes were sequenced from cultured isolates or from environmental metagenomes of recent planktonic Nostocales blooms. Phylogenomic analysis of publicly available sequences placed the new genomes among a group of 15 genomes from four continents in a distinct ADA clade (Anabaena/Dolichospermum/Aphanizomenon) within the Nostocales. This clade contains four species-level groups, two of which include members with both Anabaena-like and Aphanizomenon flos-aquae-like morphology. The genomes contain many repetitive genetic elements and a sizable pangenome, in which ABC-type transporters are highly represented. Alongside common core genes for photosynthesis, the differentiation of N₂-fixing heterocysts, and the uptake and incorporation of the major nutrients P, N and S, we identified several gene pathways in the pangenome that may contribute to niche partitioning. Genes for problematic secondary metabolites-cyanotoxins and taste-and-odor compounds-were sporadically present, as were other polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) gene clusters. By contrast, genes predicted to encode the ribosomally generated bacteriocin peptides were found in all genomes.

Hydrocarbon degradation and response of seafloor sediment bacterial community in the northern Gulf of Mexico to light Louisiana sweet crude oil

The Deepwater Horizon (DWH) blowout resulted in the deposition to the seafloor of up to 4.9% of 200 million gallons of oil released into the Gulf of Mexico. The petroleum hydrocarbon concentrations near the wellhead were high immediately after the spill, but returned to background levels a few years after the spill. Microbial communities in the seafloor are thought to be responsible for the degradation of hydrocarbons, however, our knowledge is primarily based upon gene diversity surveys and hydrocarbon concentration in field sediment samples. Here, we investigated the oil degradation potential and changes in bacterial community by amending seafloor sediment collected near the DWH site with crude oil and both oil and Corexit dispersant. Polycyclic aromatic hydrocarbons were rapidly degraded during the first 30 days of incubation, while alkanes were degraded more slowly. With the degradation of hydrocarbons, the relative abundances of Colwelliaceae, Alteromonadaceae, Methylococales, Alcanivorax, Bacteriovorax, and Phaeobacter increased remarkably. However, the abundances of oil-degrading bacteria changed with oil chemistry. Colwelliaceae decreased with increasing oil degradation, whereas Alcanivorax and Methylococcales increased considerably. We assembled seven genomes from the metagenome, including ones belonging to Colwellia, Alteromonadaceae, Rhodobacteraceae, the newly reported genus Woeseia, and candidate phylum NC10, all of which possess a repertoire of genes for hydrocarbon degradation. Moreover, genes related to hydrocarbon degradation were highly enriched in the oiled treatment, suggesting that the hydrocarbons were biodegraded, and that the indigenous microflora have a remarkable potential for the natural attenuation of spilled oil in the deep-sea surface sediment.

Insights from the Genomes of Microbes Thriving in Uranium-Enriched Sediments

Elevated uranium dose (4 g kg^-1) causes a shift in billabong sediment communities that result in the enrichment of five bacterial species. These taxa include Geobacter, Geothrix and Dyella species, as well as a novel-potentially predatory-Bacteroidetes species, and a new member of class Anaerolineae (Chloroflexi). Additionally, a population of methanogenic Methanocella species was also identified. Genomic reconstruction and metabolic examination of these taxa reveal a host of divergent life strategies and putative niche partitioning. Resistance-nodulation-division heavy metal efflux (RND-HME) transporters are implicated as potential uranium tolerance strategies among the bacterial taxa. Potential interactions, uranium tolerance and ecologically relevant catabolism are presented in a conceptual model of life in this environment.

Draft Genome Assembly of the Sheep Scab Mite, Psoroptes ovis

Sheep scab, caused by infestation with Psoroptes ovis, is highly contagious, results in intense pruritus, and represents a major welfare and economic concern. Here, we report the first draft genome assembly and gene prediction of P. ovis based on PacBio de novo sequencing. The ∼63.2-Mb genome encodes 12,041 protein-coding genes.

Hospitalized Premature Infants Are Colonized by Related Bacterial Strains with Distinct Proteomic Profiles

During the first weeks of life, microbial colonization of the gut impacts human immune system maturation and other developmental processes. In premature infants, aberrant colonization has been implicated in the onset of necrotizing enterocolitis (NEC), a life-threatening intestinal disease. To study the premature infant gut colonization process, genome-resolved metagenomics was conducted on 343 fecal samples collected during the first 3 months of life from 35 premature infants housed in a neonatal intensive care unit, 14 of whom developed NEC, and metaproteomic measurements were made on 87 samples. Microbial community composition and proteomic profiles remained relatively stable on the time scale of a week, but the proteome was more variable. Although genetically similar organisms colonized many infants, most infants were colonized by distinct strains with metabolic profiles that could be distinguished using metaproteomics. Microbiome composition correlated with infant, antibiotics administration, and NEC diagnosis. Communities were found to cluster into seven primary types, and community type switched within infants, sometimes multiple times. Interestingly, some communities sampled from the same infant at subsequent time points clustered with those of other infants. In some cases, switches preceded onset of NEC; however, no species or community type could account for NEC across the majority of infants. In addition to a correlation of protein abundances with organism replication rates, we found that organism proteomes correlated with overall community composition. Thus, this genome-resolved proteomics study demonstrated that the contributions of individual organisms to microbiome development depend on microbial community context.IMPORTANCE Humans are colonized by microbes at birth, a process that is important to health and development. However, much remains to be known about the fine-scale microbial dynamics that occur during the colonization period. We conducted a genome-resolved study of microbial community composition, replication rates, and proteomes during the first 3 months of life of both healthy and sick premature infants. Infants were found to be colonized by similar microbes, but each underwent a distinct colonization trajectory. Interestingly, related microbes colonizing different infants were found to have distinct proteomes, indicating that microbiome function is not only driven by which organisms are present, but also largely depends on microbial responses to the unique set of physiological conditions in the infant gut.

Bacteria That Make a Meal of Sulfonamide Antibiotics: Blind Spots and Emerging Opportunities

The release of sulfonamide antibiotics into the environment is alarming because the existence of these antibiotics in the environment may promote resistance in clinically relevant microorganisms, which is a potential threat to the effectiveness of antibiotic therapies. Controllable biodegradation processes are of particular significance for the inexpensive yet effective restoration of sulfonamide-contaminated environments. Cultivation-based techniques have already made great strides in successfully isolating bacteria with promising sulfonamide degradation abilities, but the implementation of these isolates in bioremediation has been limited by unknown microbial diversity, vast population responsiveness, and the impact of perturbations from open and complex environments. Advances in DNA sequencing technologies and metagenomic analyses are being used to complement the information derived from cultivation-based procedures. In this Review, we provide an overview of the growing understanding of isolated sulfonamide degraders and identify shortcomings of the prevalent literature in this field. In addition, we propose a technical paradigm that integrates experimental testing with metagenomic analysis to pave the way for improved understanding and exploitation of these ecologically important isolates. Overall, this Review aims to outline how metagenomic studies of isolated sulfonamide degraders are being applied for the advancement of bioremediation strategies for sulfonamide contamination.

Cultivation and genomics of the first freshwater SAR11 (LD12) isolate

Evolutionary transitions between fresh and salt water happen infrequently among bacterioplankton. Within the ubiquitous and highly abundant heterotrophic Alphaproteobacteria order Pelagibacterales (SAR11), most members live in marine habitats, but the LD12 subclade has evolved as a unique freshwater lineage. LD12 cells occur as some of the most dominant freshwater bacterioplankton, yet this group has remained elusive to cultivation, hampering a more thorough understanding of its biology. Here, we report the first successful isolation of an LD12 representative, strain LSUCC0530, using high-throughput dilution-to-extinction cultivation methods, and its complete genome sequence. Growth experiments corroborate ecological data suggesting active populations of LD12 in brackish water up to salinities of ~5. LSUCC0530 has the smallest closed genome thus far reported for a SAR11 strain (1.16 Mbp). The genome affirms many previous metabolic predictions from cultivation-independent analyses, like a complete Embden–Meyerhof–Parnas glycolysis pathway, but also provides novel insights, such as the first isocitrate dehydrogenase in LD12, a likely homologous recombination of malate synthase from outside of the SAR11 clade, and analogous substitutions of ion transporters with others that occur throughout the rest of the SAR11 clade. Growth data support metagenomic recruitment results suggesting temperature-based ecotype diversification within LD12. Key gene losses for osmolyte uptake provide a succinct hypothesis for the evolutionary transition of LD12 from salt to freshwater. For strain LSUCC0530, we propose the provisional nomenclature Candidatus fonsibacter ubiquis.

Phylogenetic Diversity and Single-Cell Genome Analysis of "Melainabacteria", a Non-Photosynthetic Cyanobacterial Group, in the Termite Gut

Termite guts harbor diverse yet-uncultured bacteria, including a non-photosynthetic cyanobacterial group, the class "Melainabacteria". We herein reported the phylogenetic diversity of "Melainabacteria" in the guts of diverse termites and conducted a single-cell genome analysis of a melainabacterium obtained from the gut of the termite Termes propinquus. We performed amplicon sequencing of 16S rRNA genes from the guts of 60 termite and eight cockroach species, and detected melainabacterial sequences in 48 out of the 68 insect species, albeit with low abundances (0.02-1.90%). Most of the melainabacterial sequences obtained were assigned to the order "Gastranaerophilales" and appeared to form clusters unique to termites and cockroaches. A single-cell genome of a melainabacterium, designated phylotype Tpq-Mel-01, was obtained using a fluorescence-activated cell sorter and whole genome amplification. The genome shared basic features with other melainabacterial genomes previously reconstructed from the metagenomes of human and koala feces. The bacterium had a small genome (~1.6 Mb) and possessed fermentative pathways possibly using sugars and chitobiose as carbon and energy sources, while the pathways for photosynthesis and carbon fixation were not found. The genome contained genes for flagellar components and chemotaxis; therefore, the bacterium is likely motile. A fluorescence in situ hybridization analysis showed that the cells of Tpq-Mel-01 and/or its close relatives are short rods with the dimensions of 1.1±0.2 μm by 0.5±0.1 μm; for these bacteria, we propose the novel species, "Candidatus Gastranaerophilus termiticola". Our results provide fundamental information on "Melainabacteria" in the termite gut and expand our knowledge on this underrepresented, non-photosynthetic cyanobacterial group.

Metabolic versatility of small archaea Micrarchaeota and Parvarchaeota

Small acidophilic archaea belonging to Micrarchaeota and Parvarchaeota phyla are known to physically interact with some Thermoplasmatales members in nature. However, due to a lack of cultivation and limited genomes on hand, their biodiversity, metabolisms, and physiologies remain largely unresolved. Here, we obtained 39 genomes from acid mine drainage (AMD) and hot spring environments around the world. 16S rRNA gene based analyses revealed that Parvarchaeota were only detected in AMD and hot spring habitats, while Micrarchaeota were also detected in others including soil, peat, hypersaline mat, and freshwater, suggesting a considerable higher diversity and broader than expected habitat distribution for this phylum. Despite their small genomes (0.64-1.08 Mb), these archaea may contribute to carbon and nitrogen cycling by degrading multiple saccharides and proteins, and produce ATP via aerobic respiration and fermentation. Additionally, we identified several syntenic genes with homology to those involved in iron oxidation in six Parvarchaeota genomes, suggesting their potential role in iron cycling. However, both phyla lack biosynthetic pathways for amino acids and nucleotides, suggesting that they likely scavenge these biomolecules from the environment and/or other community members. Moreover, low-oxygen enrichments in laboratory confirmed our speculation that both phyla are microaerobic/anaerobic, based on several specific genes identified in them. Furthermore, phylogenetic analyses provide insights into the close evolutionary history of energy related functionalities between both phyla with Thermoplasmatales. These results expand our understanding of these elusive archaea by revealing their involvement in carbon, nitrogen, and iron cycling, and suggest their potential interactions with Thermoplasmatales on genomic scale.

Draft Genome Sequence of a Novel Thermofilum sp. Strain from a New Zealand Hot Spring Enrichment Culture

A draft genome of a new Thermofilum sp. strain was obtained from an enrichment culture metagenome. Like its relatives, Thermofilum sp. strain NZ13 is adapted to organic-rich thermal environments and has to depend on other organisms and the environment for some key amino acids, purines, and cofactors.

Comparative genomic inference suggests mixotrophic lifestyle for Thorarchaeota

Thorarchaeota are a new archaeal phylum within the Asgard superphylum, whose ancestors have been proposed to play possible ecological roles in cellular evolution. However, little is known about the lifestyles of these uncultured archaea. To provide a better resolution of the ecological roles and metabolic capacity of Thorarchaeota, we obtained Thorarchaeota genomes reconstructed from metagenomes of different depth layers in mangrove and mudflat sediments. These genomes from deep anoxic layers suggest the presence of Thorarchaeota with the potential to degrade organic matter, fix inorganic carbon, reduce sulfur/sulfate and produce acetate. In particular, Thorarchaeota may be involved in ethanol production, nitrogen fixation, nitrite reduction, and arsenic detoxification. Interestingly, these Thorarchaeotal genomes are inferred to contain the tetrahydromethanopterin and tetrahydrofolate Wood-Ljungdahl (WL) pathways for CO₂ reduction, and the latter WL pathway appears to have originated from bacteria. These archaea are predicted to be able to use various inorganic and organic carbon sources, possessing genes inferred to encode ribulose bisphosphate carboxylase-like proteins (normally without RuBisCO activity) and a near-complete Calvin-Benson-Bassham cycle. The existence of eukaryotic selenocysteine insertion sequences and many genes for proteins previously considered eukaryote-specific in Thorarchaeota genomes provide new insights into their evolutionary roles in the origin of eukaryotic cellular complexity. Resolving the metabolic capacities of these enigmatic archaea and their origins will enhance our understanding of the origins of eukaryotes and their roles in ecosystems.

Microbiomes in extremely acidic environments: functionalities and interactions that allow survival and growth of prokaryotes at low pH

Extremely acidic environments have global distribution and can have natural or, increasingly, anthropogenic origins. Extreme acidophiles grow optimally at pH 3 or less, have multiple strategies for tolerating stresses that accompany high levels of acidity and are scattered in all three domains of the tree of life. Metagenomic studies have expanded knowledge of the diversity of extreme acidophile communities, their ecological networks and their metabolic potentials, both confirmed and inferred. High resolution compositional and functional profiling of these microbiomes have begun to reveal spatial diversity patterns at global, regional, local, zonal and micro-scales. Future integration of genomic and other meta-omic data will offer new opportunities to utilize acidic microbiomes and to engineer beneficial interactions within them in biotechnologies.

A novel Chromatiales bacterium is a potential sulfide oxidizer in multiple orders of marine sponges

Sponges are benthic filter feeders that play pivotal roles in coupling benthic-pelagic processes in the oceans that involve transformation of dissolved and particulate organic carbon and nitrogen into biomass. While the contribution of sponge holobionts to the nitrogen cycle has been recognized in past years, their importance in the sulfur cycle, both oceanic and physiological, has only recently gained attention. Sponges in general, and Theonella swinhoei in particular, harbour a multitude of associated microorganisms that could affect sulfur cycling within the holobiont. We reconstructed the genome of a Chromatiales (class Gammaproteobacteria) bacterium from a metagenomic sequence dataset of a T. swinhoei-associated microbial community. This relatively abundant bacterium has the metabolic capability to oxidize sulfide yet displays reduced metabolic potential suggestive of its lifestyle as an obligatory symbiont. This bacterium was detected in multiple sponge orders, according to similarities in key genes such as 16S rRNA and polyketide synthase genes. Due to its sulfide oxidation metabolism and occurrence in many members of the Porifera phylum, we suggest naming the newly described taxon Candidatus Porisulfidus.

Giant viruses with an expanded complement of translation system components

The discovery of giant viruses blurred the sharp division between viruses and cellular life. Giant virus genomes encode proteins considered as signatures of cellular organisms, particularly translation system components, prompting hypotheses that these viruses derived from a fourth domain of cellular life. Here we report the discovery of a group of giant viruses (Klosneuviruses) in metagenomic data. Compared with other giant viruses, the Klosneuviruses encode an expanded translation machinery, including aminoacyl transfer RNA synthetases with specificities for all 20 amino acids. Notwithstanding the prevalence of translation system components, comprehensive phylogenomic analysis of these genes indicates that Klosneuviruses did not evolve from a cellular ancestor but rather are derived from a much smaller virus through extensive gain of host genes.

Decontaminating eukaryotic genome assemblies with machine learning

BACKGROUND

High-throughput sequencing has made it theoretically possible to obtain high-quality de novo assembled genome sequences but in practice DNA extracts are often contaminated with sequences from other organisms. Currently, there are few existing methods for rigorously decontaminating eukaryotic assemblies. Those that do exist filter sequences based on nucleotide similarity to contaminants and risk eliminating sequences from the target organism.

RESULTS

We introduce a novel application of an established machine learning method, a decision tree, that can rigorously classify sequences. The major strength of the decision tree is that it can take any measured feature as input and does not require a priori identification of significant descriptors. We use the decision tree to classify de novo assembled sequences and compare the method to published protocols.

CONCLUSIONS

A decision tree performs better than existing methods when classifying sequences in eukaryotic de novo assemblies. It is efficient, readily implemented, and accurately identifies target and contaminant sequences. Importantly, a decision tree can be used to classify sequences according to measured descriptors and has potentially many uses in distilling biological datasets.

The genome of an intranuclear parasite, Paramicrosporidium saccamoebae, reveals alternative adaptations to obligate intracellular parasitism

Intracellular parasitism often results in gene loss, genome reduction, and dependence upon the host for cellular functioning. Rozellomycota is a clade comprising many such parasites and is related to the diverse, highly reduced, animal parasites, Microsporidia. We sequenced the nuclear and mitochondrial genomes of Paramicrosporidium saccamoebae [Rozellomycota], an intranuclear parasite of amoebae. A canonical fungal mitochondrial genome was recovered from P. saccamoebae that encodes genes necessary for the complete oxidative phosphorylation pathway including Complex I, differentiating it from most endoparasites including its sequenced relatives in Rozellomycota and Microsporidia. Comparative analysis revealed that P. saccamoebae shares more gene content with distantly related Fungi than with its closest relatives, suggesting that genome evolution in Rozellomycota and Microsporidia has been affected by repeated and independent gene losses, possibly as a result of variation in parasitic strategies (e.g. host and subcellular localization) or due to multiple transitions to parasitism.

Novel Microbial Diversity and Functional Potential in the Marine Mammal Oral Microbiome

The vast majority of bacterial diversity lies within phylum-level lineages called "candidate phyla," which lack isolated representatives and are poorly understood. These bacteria are surprisingly abundant in the oral cavity of marine mammals. We employed a genome-resolved metagenomic approach to recover and characterize genomes and functional potential from microbes in the oral gingival sulcus of two bottlenose dolphins (Tursiops truncatus). We detected organisms from 24 known bacterial phyla and one archaeal phylum. We also recovered genomes from two deep-branching, previously uncharacterized phylum-level lineages (here named "Candidatus Delphibacteria" and "Candidatus Fertabacteria"). The Delphibacteria lineage is found in both managed and wild dolphins; its metabolic profile suggests a capacity for denitrification and a possible role in dolphin health. We uncovered a rich diversity of predicted Cas9 proteins, including the two longest predicted Cas9 proteins to date. Notably, we identified the first type II CRISPR-Cas systems encoded by members of the Candidate Phyla Radiation. Using their spacer sequences, we subsequently identified and assembled a complete Saccharibacteria phage genome. These findings underscore the immense microbial diversity and functional potential that await discovery in previously unexplored environments.

The Threat of Antimicrobial Resistance on the Human Microbiome

Ubiquitous in nature, antimicrobial resistance (AMR) has existed long before the golden age of antimicrobials. While antimicrobial agents are beneficial to combat infection, their widespread use contributes to the increase in and emergence of novel resistant microbes in virtually all environmental niches. The human microbiome is an important reservoir of AMR with initial exposure occurring in early life. Once seeded with AMR, commensal organisms may be key contributors to the dissemination of resistance due to the interconnectedness of microbial communities. When acquired by pathogens however, AMR becomes a serious public health threat worldwide. Our ability to combat the threat of emerging resistance relies on accurate AMR detection methods and the development of therapeutics that function despite the presence of antimicrobial resistance.

Genome sequences of lower Great Lakes Microcystis sp. reveal strain-specific genes that are present and expressed in western Lake Erie blooms

Blooms of the potentially toxic cyanobacterium Microcystis are increasing worldwide. In the Laurentian Great Lakes they pose major socioeconomic, ecological, and human health threats, particularly in western Lake Erie. However, the interpretation of "omics" data is constrained by the highly variable genome of Microcystis and the small number of reference genome sequences from strains isolated from the Great Lakes. To address this, we sequenced two Microcystis isolates from Lake Erie (Microcystis aeruginosa LE3 and M. wesenbergii LE013-01) and one from upstream Lake St. Clair (M. cf aeruginosa LSC13-02), and compared these data to the genomes of seventeen Microcystis spp. from across the globe as well as one metagenome and seven metatranscriptomes from a 2014 Lake Erie Microcystis bloom. For the publically available strains analyzed, the core genome is ~1900 genes, representing ~11% of total genes in the pan-genome and ~45% of each strain's genome. The flexible genome content was related to Microcystis subclades defined by phylogenetic analysis of both housekeeping genes and total core genes. To our knowledge this is the first evidence that the flexible genome is linked to the core genome of the Microcystis species complex. The majority of strain-specific genes were present and expressed in bloom communities in Lake Erie. Roughly 8% of these genes from the lower Great Lakes are involved in genome plasticity (rapid gain, loss, or rearrangement of genes) and resistance to foreign genetic elements (such as CRISPR-Cas systems). Intriguingly, strain-specific genes from Microcystis cultured from around the world were also present and expressed in the Lake Erie blooms, suggesting that the Microcystis pangenome is truly global. The presence and expression of flexible genes, including strain-specific genes, suggests that strain-level genomic diversity may be important in maintaining Microcystis abundance during bloom events.