Metagenomic Methods for Addressing NASA's Planetary Protection Policy Requirements on Future Missions: A Workshop Report

Molecular biology methods and technologies have advanced substantially over the past decade. These new molecular methods should be incorporated among the standard tools of planetary protection (PP) and could be validated for incorporation by 2026. To address the feasibility of applying modern molecular techniques to such an application, NASA conducted a technology workshop with private industry partners, academics, and government agency stakeholders, along with NASA staff and contractors. The technical discussions and presentations of the Multi-Mission Metagenomics Technology Development Workshop focused on modernizing and supplementing the current PP assays. The goals of the workshop were to assess the state of metagenomics and other advanced molecular techniques in the context of providing a validated framework to supplement the bacterial endospore-based NASA Standard Assay and to identify knowledge and technology gaps. In particular, workshop participants were tasked with discussing metagenomics as a stand-alone technology to provide rapid and comprehensive analysis of total nucleic acids and viable microorganisms on spacecraft surfaces, thereby allowing for the development of tailored and cost-effective microbial reduction plans for each hardware item on a spacecraft. Workshop participants recommended metagenomics approaches as the only data source that can adequately feed into quantitative microbial risk assessment models for evaluating the risk of forward (exploring extraterrestrial planet) and back (Earth harmful biological) contamination. Participants were unanimous that a metagenomics workflow, in tandem with rapid targeted quantitative (digital) PCR, represents a revolutionary advance over existing methods for the assessment of microbial bioburden on spacecraft surfaces. The workshop highlighted low biomass sampling, reagent contamination, and inconsistent bioinformatics data analysis as key areas for technology development. Finally, it was concluded that implementing metagenomics as an additional workflow for addressing concerns of NASA's robotic mission will represent a dramatic improvement in technology advancement for PP and will benefit future missions where mission success is affected by backward and forward contamination.

Dataset of 143 metagenome-assembled genomes from the Arctic and Atlantic Oceans, including 21 for eukaryotic organisms

This article presents metagenome-assembled genomes (MAGs) for both eukaryotic and prokaryotic organisms originating from the Arctic and Atlantic oceans, along with gene prediction and functional annotation for MAGs from both domains. Eleven samples from the chlorophyll-a maximum layer of the surface ocean were collected during two cruises in 2012; six from the Arctic in June-July on ARK-XXVII/1 (PS80), and five from the Atlantic in November on ANT-XXIX/1 (PS81). Sequencing and assembly was carried out by the Joint Genome Institute (JGI), who provide annotation of the assembled sequences, and 122 MAGs for prokaryotic organisms. A subsequent binning process identified 21 MAGs for eukaryotic organisms, mostly identified as Mamiellophyceae or Bacillariophyceae. The data for each MAG includes sequences in FASTA format, and tables of functional annotation of genes. For eukaryotic MAGs, transcript and protein sequences for predicted genes are available. A spreadsheet is provided summarising quality measures and taxonomic classifications for each MAG. These data provide draft genomes for uncultured marine microbes, including some of the first MAGs for polar eukaryotes, and can provide reference genetic data for these environments, or used in genomics-based comparison between environments.

Substrate availability and not thermal acclimation controls microbial temperature sensitivity response to long-term warming

Microbes are responsible for cycling carbon (C) through soils, and predicted changes in soil C stocks under climate change are highly sensitive to shifts in the mechanisms assumed to control the microbial physiological response to warming. Two mechanisms have been suggested to explain the long-term warming impact on microbial physiology: microbial thermal acclimation and changes in the quantity and quality of substrates available for microbial metabolism. Yet studies disentangling these two mechanisms are lacking. To resolve the drivers of changes in microbial physiology in response to long-term warming, we sampled soils from 13- and 28-year-old soil warming experiments in different seasons. We performed short-term laboratory incubations across a range of temperatures to measure the relationships between temperature sensitivity of physiology (growth, respiration, carbon use efficiency, and extracellular enzyme activity) and the chemical composition of soil organic matter. We observed apparent thermal acclimation of microbial respiration, but only in summer, when warming had exacerbated the seasonally-induced, already small dissolved organic matter pools. Irrespective of warming, greater quantity and quality of soil carbon increased the extracellular enzymatic pool and its temperature sensitivity. We propose that fresh litter input into the system seasonally cancels apparent thermal acclimation of C-cycling processes to decadal warming. Our findings reveal that long-term warming has indirectly affected microbial physiology via reduced C availability in this system, implying that earth system models including these negative feedbacks may be best suited to describe long-term warming effects on these soils.

Metagenome-assembled genomes of phytoplankton microbiomes from the Arctic and Atlantic Oceans

BACKGROUND

Phytoplankton communities significantly contribute to global biogeochemical cycles of elements and underpin marine food webs. Although their uncultured genomic diversity has been estimated by planetary-scale metagenome sequencing and subsequent reconstruction of metagenome-assembled genomes (MAGs), this approach has yet to be applied for complex phytoplankton microbiomes from polar and non-polar oceans consisting of microbial eukaryotes and their associated prokaryotes.

RESULTS

Here, we have assembled MAGs from chlorophyll a maximum layers in the surface of the Arctic and Atlantic Oceans enriched for species associations (microbiomes) with a focus on pico- and nanophytoplankton and their associated heterotrophic prokaryotes. From 679 Gbp and estimated 50 million genes in total, we recovered 143 MAGs of medium to high quality. Although there was a strict demarcation between Arctic and Atlantic MAGs, adjacent sampling stations in each ocean had 51-88% MAGs in common with most species associations between Prasinophytes and Proteobacteria. Phylogenetic placement revealed eukaryotic MAGs to be more diverse in the Arctic whereas prokaryotic MAGs were more diverse in the Atlantic Ocean. Approximately 70% of protein families were shared between Arctic and Atlantic MAGs for both prokaryotes and eukaryotes. However, eukaryotic MAGs had more protein families unique to the Arctic whereas prokaryotic MAGs had more families unique to the Atlantic.

CONCLUSION

Our study provides a genomic context to complex phytoplankton microbiomes to reveal that their community structure was likely driven by significant differences in environmental conditions between the polar Arctic and warm surface waters of the tropical and subtropical Atlantic Ocean. Video Abstract.

The biogeographic differentiation of algal microbiomes in the upper ocean from pole to pole

Eukaryotic phytoplankton are responsible for at least 20% of annual global carbon fixation. Their diversity and activity are shaped by interactions with prokaryotes as part of complex microbiomes. Although differences in their local species diversity have been estimated, we still have a limited understanding of environmental conditions responsible for compositional differences between local species communities on a large scale from pole to pole. Here, we show, based on pole-to-pole phytoplankton metatranscriptomes and microbial rDNA sequencing, that environmental differences between polar and non-polar upper oceans most strongly impact the large-scale spatial pattern of biodiversity and gene activity in algal microbiomes. The geographic differentiation of co-occurring microbes in algal microbiomes can be well explained by the latitudinal temperature gradient and associated break points in their beta diversity, with an average breakpoint at 14 °C ± 4.3, separating cold and warm upper oceans. As global warming impacts upper ocean temperatures, we project that break points of beta diversity move markedly pole-wards. Hence, abrupt regime shifts in algal microbiomes could be caused by anthropogenic climate change.

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan

Thermoflexus hugenholtzii JAD2^T, the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2^T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2^T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing “Candidatus Thermoflexus japonica,” “Candidatus Thermoflexus tengchongensis,” and “Candidatus Thermoflexus sinensis.” Genomics was integrated with targeted exometabolomics and ¹³C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2^T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2^T using ¹³C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific ¹³C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

Insights into the dynamics between viruses and their hosts in a hot spring microbial mat

Our current knowledge of host-virus interactions in biofilms is limited to computational predictions based on laboratory experiments with a small number of cultured bacteria. However, natural biofilms are diverse and chiefly composed of uncultured bacteria and archaea with no viral infection patterns and lifestyle predictions described to date. Herein, we predict the first DNA sequence-based host-virus interactions in a natural biofilm. Using single-cell genomics and metagenomics applied to a hot spring mat of the Cone Pool in Mono County, California, we provide insights into virus-host range, lifestyle and distribution across different mat layers. Thirty-four out of 130 single cells contained at least one viral contig (26%), which, together with the metagenome-assembled genomes, resulted in detection of 59 viruses linked to 34 host species. Analysis of single-cell amplification kinetics revealed a lack of active viral replication on the single-cell level. These findings were further supported by mapping metagenomic reads from different mat layers to the obtained host-virus pairs, which indicated a low copy number of viral genomes compared to their hosts. Lastly, the metagenomic data revealed high layer specificity of viruses, suggesting limited diffusion to other mat layers. Taken together, these observations indicate that in low mobility environments with high microbial abundance, lysogeny is the predominant viral lifestyle, in line with the previously proposed "Piggyback-the-Winner" theory.

Terabase-scale metagenome coassembly with MetaHipMer

Metagenome sequence datasets can contain terabytes of reads, too many to be coassembled together on a single shared-memory computer; consequently, they have only been assembled sample by sample (multiassembly) and combining the results is challenging. We can now perform coassembly of the largest datasets using MetaHipMer, a metagenome assembler designed to run on supercomputers and large clusters of compute nodes. We have reported on the implementation of MetaHipMer previously; in this paper we focus on analyzing the impact of very large coassembly. In particular, we show that coassembly recovers a larger genome fraction than multiassembly and enables the discovery of more complete genomes, with lower error rates, whereas multiassembly recovers more dominant strain variation. Being able to coassemble a large dataset does not preclude one from multiassembly; rather, having a fast, scalable metagenome assembler enables a user to more easily perform coassembly and multiassembly, and assemble both abundant, high strain variation genomes, and low-abundance, rare genomes. We present several assemblies of terabyte datasets that could never be coassembled before, demonstrating MetaHipMer’s scaling power. MetaHipMer is available for public use under an open source license and all datasets used in the paper are available for public download.

Functional Signatures of the Epiphytic Prokaryotic Microbiome of Agaves and Cacti

Microbial symbionts account for survival, development, fitness and evolution of eukaryotic hosts. These microorganisms together with their host form a biological unit known as holobiont. Recent studies have revealed that the holobiont of agaves and cacti comprises a diverse and structured microbiome, which might be important for its adaptation to drylands. Here, we investigated the functional signatures of the prokaryotic communities of the soil and the episphere, that includes the rhizosphere and phyllosphere, associated with the cultivated Agave tequilana and the native and sympatric Agave salmiana, Opuntia robusta and Myrtillocactus geometrizans by mining shotgun metagenomic data. Consistent with previous phylogenetic profiling, we found that Proteobacteria, Actinobacteria and Firmicutes were the main represented phyla in the episphere of agaves and cacti, and that clustering of metagenomes correlated with the plant compartment. In native plants, genes related to aerobic anoxygenic phototrophy and photosynthesis were enriched in the phyllosphere and soil, while genes coding for biofilm formation and quorum sensing were enriched in both epiphytic communities. In the episphere of cultivated A. tequilana fewer genes were identified, but they belonged to similar pathways than those found in native plants. A. tequilana showed a depletion in several genes belonging to carbon metabolism, secondary metabolite biosynthesis and xenobiotic degradation suggesting that its lower microbial diversity might be linked to functional losses. However, this species also showed an enrichment in biofilm and quorum sensing in the epiphytic compartments, and evidence for nitrogen fixation in the rhizosphere. Aerobic anoxygenic phototrophic markers were represented by Rhizobiales (Methylobacterium) and Rhodospirillales (Belnapia) in the phyllosphere, while photosystem genes were widespread in Bacillales and Cyanobacteria. Nitrogen fixation and biofilm formation genes were mostly related to Proteobacteria. These analyses support the idea of niche differentiation in the rhizosphere and phyllosphere of agaves and cacti and shed light on the potential mechanisms by which epiphytic microbial communities survive and colonize plants of arid and semiarid ecosystems. This study establishes a guideline for testing the relevance of the identified functional traits on the microbial community and the plant fitness.

Triplicate PCR reactions for 16S rRNA gene amplicon sequencing are unnecessary

Conventional wisdom holds that PCR amplification for sequencing should employ pooled replicate reactions to reduce bias due to jackpot effects and chimera formation. However, modern amplicon data analysis employs methods that may be less sensitive to such artifacts. Here we directly compare results from single versus triplicate reactions for 16S amplicon sequencing and find no significant impact of adopting a less labor-intensive single-reaction protocol.

Effect of the macroalgae Asparagopsis taxiformis on methane production and rumen microbiome assemblage

BACKGROUND

Recent studies using batch-fermentation suggest that the red macroalgae Asparagopsis taxiformis has the potential to reduce methane (CH4) production from beef cattle by up to ~ 99% when added to Rhodes grass hay; a common feed in the Australian beef industry. These experiments have shown significant reductions in CH4 without compromising other fermentation parameters (i.e. volatile fatty acid production) with A. taxiformis organic matter (OM) inclusion rates of up to 5%. In the study presented here, A. taxiformis was evaluated for its ability to reduce methane production from dairy cattle fed a mixed ration widely utilized in California, the largest milk producing state in the US.

RESULTS

Fermentation in a semi-continuous in-vitro rumen system suggests that A. taxiformis can reduce methane production from enteric fermentation in dairy cattle by 95% when added at a 5% OM inclusion rate without any obvious negative impacts on volatile fatty acid production. High-throughput 16S ribosomal RNA (rRNA) gene amplicon sequencing showed that seaweed amendment effects rumen microbiome consistent with the Anna Karenina hypothesis, with increased β-diversity, over time scales of approximately 3 days. The relative abundance of methanogens in the fermentation vessels amended with A. taxiformis decreased significantly compared to control vessels, but this reduction in methanogen abundance was only significant when averaged over the course of the experiment. Alternatively, significant reductions of CH4 in the A. taxiformis amended vessels was measured in the early stages of the experiment. This suggests that A. taxiformis has an immediate effect on the metabolic functionality of rumen methanogens whereas its impact on microbiome assemblage, specifically methanogen abundance, is delayed.

CONCLUSIONS

The methane reducing effect of A. taxiformis during rumen fermentation makes this macroalgae a promising candidate as a biotic methane mitigation strategy for dairy cattle. But its effect in-vivo (i.e. in dairy cattle) remains to be investigated in animal trials. Furthermore, to obtain a holistic understanding of the biochemistry responsible for the significant reduction of methane, gene expression profiles of the rumen microbiome and the host animal are warranted.

Aerosol Mycobacterium tuberculosis infection causes rapid loss of diversity in gut microbiota

Mycobacterium tuberculosis is an important human pathogen, and yet diagnosis remains challenging. Little research has focused on the impact of M. tuberculosis on the gut microbiota, despite the significant immunological and homeostatic functions of the gastrointestinal tract. To determine the effect of M. tuberculosis infection on the gut microbiota, we followed mice from M. tuberculosis aerosol infection until death, using 16S rRNA sequencing. We saw a rapid change in the gut microbiota in response to infection, with all mice showing a loss and then recovery of microbial community diversity, and found that pre-infection samples clustered separately from post-infection samples, using ecological beta-diversity measures. The effect on the fecal microbiota was observed as rapidly as six days following lung infection. Analysis of additional mice infected by a different M. tuberculosis strain corroborated these results, together demonstrating that the mouse gut microbiota significantly changes with M. tuberculosis infection.