Tectonic and geological setting influence hot spring microbiology.

Hydrothermal systems form at divergent and convergent boundaries of lithospheric plates and within plates due to weakened crust and mantle plumes, playing host to diverse microbial ecosystems. Little is known of how differences in tectonic setting influence the geochemical and microbial compositions of these hydrothermal ecosystems. Here, coordinated geochemical and microbial community analyses were conducted on 87 high-temperature (>65°C) water and sediment samples from hot springs in Yellowstone National Park, Wyoming, USA (n = 41; mantle plume setting), Iceland (n = 41, divergent boundary), and Japan (n = 5; convergent boundary). Region-specific variation in geochemistry and sediment-associated 16S rRNA gene amplicon sequence variant (ASV) composition was observed, with 16S rRNA gene assemblages being nearly completely distinguished by region and pH being the most explanatory parameter within regions. Several low abundance ASVs exhibited cosmopolitan distributions across regions, while most high-abundance ASVs were only identified in specific regions. The presence of some cosmopolitan ASVs across regions argues against dispersal limitation primarily shaping the distribution of taxa among regions. Rather, the results point to local tectonic and geologic characteristics shaping the geochemistry of continental hydrothermal systems that then select for distinct microbial assemblages. These results provide new insights into the co-evolution of hydrothermal systems and their microbial communities.

Multi-omics analysis provides insight to the Ignicoccus hospitalis-Nanoarchaeum equitans association.

BACKGROUND: Studies of interspecies interactions are inherently difficult due to the complex mechanisms which enable these relationships. A model system for studying interspecies interactions is the marine hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans. Recent independently-conducted 'omics' analyses have generated insights into the molecular factors modulating this association. However, significant questions remain about the nature of the interactions between these archaea. METHODS: We jointly analyzed multiple levels of omics datasets obtained from published, independent transcriptomics, proteomics, and metabolomics analyses. DAVID identified functionally-related groups enriched when I. hospitalis is grown alone or in co-culture with N. equitans. Enriched molecular pathways were subsequently visualized using interaction maps generated using STRING. RESULTS: Key findings of our multi-level omics analysis indicated that I. hospitalis provides precursors to N. equitans for energy metabolism. Analysis indicated an overall reduction in diversity of metabolic precursors in the I. hospitalis-N. equitans co-culture, which has been connected to the differential use of ribosomal subunits and was previously unnoticed. We also identified differences in precursors linked to amino acid metabolism, NADH metabolism, and carbon fixation, providing new insights into the metabolic adaptions of I. hospitalis enabling the growth of N. equitans. CONCLUSIONS: This multi-omics analysis builds upon previously identified cellular patterns while offering new insights into mechanisms that enable the I. hospitalis-N. equitans association. GENERAL SIGNIFICANCE: Our study applies statistical and visualization techniques to a mixed-source omics dataset to yield a more global insight into a complex system, that was not readily discernable from separate omics studies.

Rescuing Those Left Behind: Recovering and Characterizing Underdigested Membrane and Hydrophobic Proteins To Enhance Proteome Measurement Depth.

The marine archaeon Nanoarchaeum equitans is dependent on direct physical contact with its host, the hyperthermophile Ignicoccus hospitalis. As this interaction is thought to be membrane-associated, involving a myriad of membrane-anchored proteins, proteomic efforts to better characterize this difficult to analyze interface are paramount to uncovering the mechanism of their association. By extending multienzyme digestion strategies that use sample filtration to recover underdigested proteins for reprocessing/consecutive proteolytic digestion, we applied chymotrypsin to redigest the proteinaceous material left over after initial proteolysis with trypsin of sodium dodecyl sulfate (SDS)-extracted I. hospitalis-N. equitans proteins. Using this method, we show that proteins with increased hydrophobic character, including membrane proteins with multiple transmembrane helices, are enriched and recovered in the underdigested fraction. Chymotryptic reprocessing provided significant sequence coverage gains in both soluble and hydrophobic proteins alike, with the latter benefiting more so in terms of membrane protein representation. These gains were despite a large proportion of high-quality peptide spectra remaining unassigned in the underdigested fraction suggesting high levels of protein modification on these often surface-exposed proteins. Importantly, these gains were achieved without applying extensive fractionation strategies usually required for thorough characterization of membrane-associated proteins and were facilitated by the generation of a distinct, complementary set of peptides that aid in both the identification and quantitation of this important, under-represented class of proteins.

Characterization of Fatty Acids in Crenarchaeota by GC-MS and NMR.

Lipids composed of condensed isoprenyl units connected to glycerol backbones by ether linkages are a distinguishing feature of Archaea. Data suggesting that fatty acids with linear hydrocarbon chains are present in some Archaea have been available for decades. However, lack of genomic and biochemical evidence for the metabolic machinery required to synthesize and degrade fatty acids has left the field unclear on this potentially significant biochemical aspect. Because lipids are energy currency and cell signaling molecules, their presence in Archaea is significant for understanding archaeal biology. A recent large-scale bioinformatics analysis reignited the debate as to the importance of fatty acids in Archaea by presenting genetic evidence for the presence of enzymes required for anabolic and catabolic fatty acid metabolism across the archaeal domain. Here, we present direct biochemical evidence from gas chromatography-mass spectrometry (GC-MS) and nuclear magnetic resonance (NMR) spectroscopy for the presence of fatty acids in two members of the Crenarchaeota, Sulfolobus solfataricus and Ignicoccus hospitalis. This is the first report providing biochemical data for the existence of fatty acids in these Crenarchaeota, opening new discussions on energy balance and the potential for the discovery of new thermostable enzymes for industry.

Expression of a xanthine permease and phosphate transporter in cultures and field populations of the harmful alga Aureococcus anophagefferens: tracking nutritional deficiency during brown tides.

Targeted gene expression using quantitative reverse transcription polymerase chain reaction (qRT-PCR) was employed to track patterns in the expression of genes indicative of nitrogen or phosphorus deficiency in the brown tide-forming alga Aureococcus anophagefferens. During culture experiments, a xanthine/uracil/vitamin C permease (XUV) was upregulated ∼20-fold under nitrogen-deficient conditions relative to a nitrogen-replete control and rapidly returned to nitrogen-replete levels after nitrogen-deficient cells were resupplied with nitrate or ammonium. It was not responsive to phosphorus deficiency. Expression of an inorganic phosphate transporter (PTA3) was enriched ∼10-fold under phosphorus-deficient conditions relative to a phosphorus-replete control, and this signal was rapidly lost upon phosphate resupply. PTA3 was not upregulated by nitrogen deficiency. Natural A. anophagefferens populations from a dense brown tide that occurred in Long Island, NY, in 2009 were assayed for XUV and PTA3 expression and compared with nutrient concentrations over the peak of a bloom. Patterns in XUV expression were consistent with nitrogen-replete growth, never reaching the values observed in N-deficient cultures. PTA3 expression was highest prior to peak bloom stages, reaching expression levels within the range of P-deficient cultures. These data highlight the value of molecular-level assessments of nutrient deficiency and suggest that phosphorus deficiency could play a role in the dynamics of destructive A. anophagefferens blooms.

Niche of harmful alga Aureococcus anophagefferens revealed through ecogenomics.

Harmful algal blooms (HABs) cause significant economic and ecological damage worldwide. Despite considerable efforts, a comprehensive understanding of the factors that promote these blooms has been lacking, because the biochemical pathways that facilitate their dominance relative to other phytoplankton within specific environments have not been identified. Here, biogeochemical measurements showed that the harmful alga Aureococcus anophagefferens outcompeted co-occurring phytoplankton in estuaries with elevated levels of dissolved organic matter and turbidity and low levels of dissolved inorganic nitrogen. We subsequently sequenced the genome of A. anophagefferens and compared its gene complement with those of six competing phytoplankton species identified through metaproteomics. Using an ecogenomic approach, we specifically focused on gene sets that may facilitate dominance within the environmental conditions present during blooms. A. anophagefferens possesses a larger genome (56 Mbp) and has more genes involved in light harvesting, organic carbon and nitrogen use, and encoding selenium- and metal-requiring enzymes than competing phytoplankton. Genes for the synthesis of microbial deterrents likely permit the proliferation of this species, with reduced mortality losses during blooms. Collectively, these findings suggest that anthropogenic activities resulting in elevated levels of turbidity, organic matter, and metals have opened a niche within coastal ecosystems that ideally suits the unique genetic capacity of A. anophagefferens and thus, has facilitated the proliferation of this and potentially other HABs.

Environmental DNA sequencing dataset from Lake Erie algal blooms using Oxford Nanopore MinION.

Here we describe a publicly available environmental DNA (eDNA) sequence dataset, consisting of samples collected from a National Oceanic and Atmospheric Administration (NOAA) Great Lakes Environmental Research Laboratory (GLERL) on Lake Erie. We sequenced samples drawn from before, during, and after a 2019 Microcystis harmful algal bloom (HAB) using 3rd generation sequencing with the Oxford Nanopore MinION device. We classified the eDNA reads taxonomically, and estimated the abundances of all taxa in each sample. While the taxonomic data showed evidence of significant human and E. coli contamination, we found abundant Mycrocystis, especially in the samples drawn from bloom environments. The raw sequence data are available in the Sequence Read Archive (SRA) under accession number PRJNA812770. HABs pose a significant and increasing risk, both to human health and to the Blue Economy, and genomic approaches to early detection promise to help mitigate these risks. As such, this dataset could be of interest to freshwater ecology research teams, or any stakeholders interested in the detection and mitigation of HABs.

Nutrient-regulated transcriptional responses in the brown tide-forming alga Aureococcus anophagefferens.

Long-SAGE (serial analysis of gene expression) was used to profile the transcriptome of the brown tide-forming alga, Aureococcus anophagefferens, under nutrient replete (control), and nitrogen (N) and phosphorus (P) deficiency to understand how this organism responds at the transcriptional level to varying nutrient conditions. This approach has aided A. anophagefferens genome annotation efforts and identified a suite of genes upregulated by N and P deficiency, some of which have known roles in nutrient metabolism. Genes upregulated under N deficiency include an ammonium transporter, an acetamidase/formamidase and two peptidases. This suggests an ability to utilize reduced N compounds and dissolved organic nitrogen, supporting the hypothesized importance of these N sources in A. anophagefferens bloom formation. There are also a broad suite of P-regulated genes, including an alkaline phosphatase, and two 5'-nucleotidases, suggesting A. anophagefferens may use dissolved organic phosphorus under low phosphate conditions. These N- and P-regulated genes may be important targets for exploring nutrient controls on bloom formation in field populations.

Genomic signatures of Lake Erie bacteria suggest interaction in the Microcystis phycosphere.

Microbial interactions in harmful algal bloom (HAB) communities have been examined in marine systems, but are poorly studied in fresh waters. To investigate HAB-microbe interactions, we isolated bacteria with close associations to bloom-forming cyanobacteria, Microcystis spp., during a 2017 bloom in the western basin of Lake Erie. The genomes of five isolates (Exiguobacterium sp. JMULE1, Enterobacter sp. JMULE2, Deinococcus sp. JMULE3, Paenibacillus sp. JMULE4, and Acidovorax sp. JMULE5.) were sequenced on a PacBio Sequel system. These genomes ranged in size from 3.1 Mbp (Exiguobacterium sp. JMULE1) to 5.7 Mbp (Enterobacter sp. JMULE2). The genomes were analyzed for genes relating to critical metabolic functions, including nitrogen reduction and carbon utilization. All five of the sequenced genomes contained genes that could be used in potential signaling and nutrient exchange between the bacteria and cyanobacteria such as Microcystis. Gene expression signatures of algal-derived carbon utilization for two isolates were identified in Microcystis blooms in Lake Erie and Lake Tai (Taihu) at low levels, suggesting these organisms are active and may have a functional role during Microcystis blooms in aggregates, but were largely missing from whole water samples. These findings build on the growing evidence that the bacterial microbiome associated with bloom-forming algae have the functional potential to contribute to nutrient exchange within bloom communities and interact with important bloom formers like Microcystis.

Transcriptional Shifts Highlight the Role of Nutrients in Harmful Brown Tide Dynamics.

Harmful algal blooms (HABs) threaten ecosystems and human health worldwide. Controlling nitrogen inputs to coastal waters is a common HAB management strategy, as nutrient concentrations often suggest coastal blooms are nitrogen-limited. However, defining best nutrient management practices is a long-standing challenge: in part, because of difficulties in directly tracking the nutritional physiology of harmful species in mixed communities. Using metatranscriptome sequencing and incubation experiments, we addressed this challenge by assaying the in situ physiological ecology of the ecosystem destructive alga, Aureococcus anophagefferens. Here we show that gene markers of phosphorus deficiency were expressed in situ, and modulated by the enrichment of phosphorus, which was consistent with the observed growth rate responses. These data demonstrate the importance of phosphorus in controlling brown-tide dynamics, suggesting that phosphorus, in addition to nitrogen, should be evaluated in the management and mitigation of these blooms. Given that nutrient concentrations alone were suggestive of a nitrogen-limited ecosystem, this study underscores the value of directly assaying harmful algae in situ for the development of management strategies.

Virus-host relationships of marine single-celled eukaryotes resolved from metatranscriptomics.

Establishing virus-host relationships has historically relied on culture-dependent approaches. Here we report on the use of marine metatranscriptomics to probe virus-host relationships. Statistical co-occurrence analyses of dsDNA, ssRNA and dsRNA viral markers of polyadenylation-selected RNA sequences from microbial communities dominated by Aureococcus anophagefferens (Quantuck Bay, NY), and diatoms (Narragansett Bay, RI) show active infections by diverse giant viruses (NCLDVs) associated with algal and nonalgal hosts. Ongoing infections of A. anophagefferens by a known Mimiviridae (AaV) occur during bloom peak and decline. Bloom decline is also accompanied by increased activity of viruses other than AaV, including (+) ssRNA viruses. In Narragansett Bay, increased temporal resolution reveals active NCLDVs with both 'boom-and-bust' and 'steady-state infection'-like ecologies that include known as well as novel virus-host interactions. Our approach offers a method for screening active viral infections and develops links between viruses and their potential hosts in situ. Our observations further demonstrate that previously unknown virus-host relationships in marine systems are abundant.

A Complex Endomembrane System in the Archaeon Ignicoccus hospitalis Tapped by Nanoarchaeum equitans.

Based on serial sectioning, focused ion beam scanning electron microscopy (FIB/SEM), and electron tomography, we depict in detail the highly unusual anatomy of the marine hyperthermophilic crenarchaeon, Ignicoccus hospitalis. Our data support a complex and dynamic endomembrane system consisting of cytoplasmic protrusions, and with secretory function. Moreover, we reveal that the cytoplasm of the putative archaeal ectoparasite Nanoarchaeum equitans can get in direct contact with this endomembrane system, complementing and explaining recent proteomic, transcriptomic and metabolomic data on this inter-archaeal relationship. In addition, we identified a matrix of filamentous structures and/or tethers in the voluminous inter-membrane compartment (IMC) of I. hospitalis, which might be responsible for membrane dynamics. Overall, this unusual cellular compartmentalization, ultrastructure and dynamics in an archaeon that belongs to the recently proposed TACK superphylum prompts speculation that the eukaryotic endomembrane system might originate from Archaea.

Genomics-informed isolation and characterization of a symbiotic Nanoarchaeota system from a terrestrial geothermal environment.

Biological features can be inferred, based on genomic data, for many microbial lineages that remain uncultured. However, cultivation is important for characterizing an organism's physiology and testing its genome-encoded potential. Here we use single-cell genomics to infer cultivation conditions for the isolation of an ectosymbiotic Nanoarchaeota ('Nanopusillus acidilobi') and its host (Acidilobus, a crenarchaeote) from a terrestrial geothermal environment. The cells of 'Nanopusillus' are among the smallest known cellular organisms (100-300 nm). They appear to have a complete genetic information processing machinery, but lack almost all primary biosynthetic functions as well as respiration and ATP synthesis. Genomic and proteomic comparison with its distant relative, the marine Nanoarchaeum equitans illustrate an ancient, common evolutionary history of adaptation of the Nanoarchaeota to ectosymbiosis, so far unique among the Archaea.

Untargeted metabolomics studies employing NMR and LC-MS reveal metabolic coupling between Nanoarcheum equitans and its archaeal host Ignicoccus hospitalis.

Interspecies interactions are the basis of microbial community formation and infectious diseases. Systems biology enables the construction of complex models describing such interactions, leading to a better understanding of disease states and communities. However, before interactions between complex organisms can be understood, metabolic and energetic implications of simpler real-world host-microbe systems must be worked out. To this effect, untargeted metabolomics experiments were conducted and integrated with proteomics data to characterize key molecular-level interactions between two hyperthermophilic microbial species, both of which have reduced genomes. Metabolic changes and transfer of metabolites between the archaea Ignicoccus hospitalis and Nanoarcheum equitans were investigated using integrated LC-MS and NMR metabolomics. The study of such a system is challenging, as no genetic tools are available, growth in the laboratory is challenging, and mechanisms by which they interact are unknown. Together with information about relative enzyme levels obtained from shotgun proteomics, the metabolomics data provided useful insights into metabolic pathways and cellular networks of I. hospitalis that are impacted by the presence of N. equitans, including arginine, isoleucine, and CTP biosynthesis. On the organismal level, the data indicate that N. equitans exploits metabolites generated by I. hospitalis to satisfy its own metabolic needs. This finding is based on N. equitans's consumption of a significant fraction of the metabolite pool in I. hospitalis that cannot solely be attributed to increased biomass production for N. equitans. Combining LC-MS and NMR metabolomics datasets improved coverage of the metabolome and enhanced the identification and quantitation of cellular metabolites.

Life on the edge: functional genomic response of Ignicoccus hospitalis to the presence of Nanoarchaeum equitans.

The marine hyperthermophilic crenarchaeon Ignicoccus hospitalis supports the propagation on its surface of Nanoarchaeum equitans, an evolutionarily enigmatic archaeon that resembles highly derived parasitic and symbiotic bacteria. The cellular and molecular mechanisms that enable this interarchaea relationship and the intimate physiologic consequences to I. hospitalis are unknown. Here, we used concerted proteomic and transcriptomic analyses to probe into the functional genomic response of I. hospitalis as N. equitans multiplies on its surface. The expression of over 97% of the genes was detected at mRNA level and over 80% of the predicted proteins were identified and their relative abundance measured by proteomics. These indicate that little, if any, of the host genomic information is silenced during growth in the laboratory. The primary response to N. equitans was at the membrane level, with increases in relative abundance of most protein complexes involved in energy generation as well as that of several transporters and proteins involved in cellular membrane stabilization. Similar upregulation was observed for genes and proteins involved in key metabolic steps controlling nitrogen and carbon metabolism, although the overall biosynthetic pathways were marginally impacted. Proliferation of N. equitans resulted, however, in selective downregulation of genes coding for transcription factors and replication and cell cycle control proteins as I. hospitalis shifted its physiology from its own cellular growth to that of its ectosymbiont/parasite. The combination of these multiomic approaches provided an unprecedented level of detail regarding the dynamics of this interspecies interaction, which is especially pertinent as these organisms are not genetically tractable.

Proteome changes driven by phosphorus deficiency and recovery in the brown tide-forming alga Aureococcus anophagefferens.

Shotgun mass spectrometry was used to detect proteins in the harmful alga, Aureococcus anophagefferens, and monitor their relative abundance across nutrient replete (control), phosphate-deficient (-P) and -P refed with phosphate (P-refed) conditions. Spectral counting techniques identified differentially abundant proteins and demonstrated that under phosphate deficiency, A. anophagefferens increases proteins involved in both inorganic and organic phosphorus (P) scavenging, including a phosphate transporter, 5'-nucleotidase, and alkaline phosphatase. Additionally, an increase in abundance of a sulfolipid biosynthesis protein was detected in -P and P-refed conditions. Analysis of the polar membrane lipids showed that cellular concentrations of the sulfolipid sulphoquinovosyldiacylglycerol (SQDG) were nearly two-fold greater in the -P condition versus the control condition, while cellular phospholipids were approximately 8-fold less. Transcript and protein abundances were more tightly coupled for gene products involved in P metabolism compared to those involved in a range of other metabolic functions. Comparison of protein abundances between the -P and P-refed conditions identified differences in the timing of protein degradation and turnover. This suggests that culture studies examining nutrient starvation responses will be valuable in interpreting protein abundance patterns for cellular nutritional status and history in metaproteomic datasets.

Long serial analysis of gene expression for gene discovery and transcriptome profiling in the widespread marine coccolithophore Emiliania huxleyi.

The abundant and widespread coccolithophore Emiliania huxleyi plays an important role in mediating CO2 exchange between the ocean and the atmosphere through its impact on marine photosynthesis and calcification. Here, we use long serial analysis of gene expression (SAGE) to identify E. huxleyi genes responsive to nitrogen (N) or phosphorus (P) starvation. Long SAGE is an elegant approach for examining quantitative and comprehensive gene expression patterns without a priori knowledge of gene sequences via the detection of 21-bp nucleotide sequence tags. E. huxleyi appears to have a robust transcriptional-level response to macronutrient deficiency, with 42 tags uniquely present or up-regulated twofold or greater in the N-starved library and 128 tags uniquely present or up-regulated twofold or greater in the P-starved library. The expression patterns of several tags were validated with reverse transcriptase PCR. Roughly 48% of these differentially expressed tags could be mapped to publicly available genomic or expressed sequence tag (EST) sequence data. For example, in the P-starved library a number of the tags mapped to genes with a role in P scavenging, including a putative phosphate-repressible permease and a putative polyphosphate synthetase. In short, the long SAGE analyses have (i) identified many new differentially regulated gene sequences, (ii) assigned regulation data to EST sequences with no database homology and unknown function, and (iii) highlighted previously uncharacterized aspects of E. huxleyi N and P physiology. To this end, our long SAGE libraries provide a new public resource for gene discovery and transcriptional analysis in this biogeochemically important marine organism.