Publications

@article{M2021,

title = {Genomic Variation Influences Methanothermococcus Fitness in Marine Hydrothermal Systems},

author = {Hoffert, M. and Anderson, R.E. and Reveillaud, J. and Murphy, L.G. and Stepanauskas, R. and Huber, J.A.},

url = {https://www.frontiersin.org/articles/10.3389/fmicb.2021.714920/full},

doi = {https://doi.org/10.3389/fmicb.2021.714920},

year  = {2021},

date = {2021-08-20},

urldate = {2021-08-20},

journal = {Frontiers in Microbiology},

volume = {12},

pages = {714920},

abstract = {Hydrogenotrophic methanogens are ubiquitous chemoautotrophic archaea inhabiting globally distributed deep-sea hydrothermal vent ecosystems and associated subseafloor niches within the rocky subseafloor, yet little is known about how they adapt and diversify in these habitats. To determine genomic variation and selection pressure within methanogenic populations at vents, we examined five Methanothermococcus single cell amplified genomes (SAGs) in conjunction with 15 metagenomes and 10 metatranscriptomes from venting fluids at two geochemically distinct hydrothermal vent fields on the Mid-Cayman Rise in the Caribbean Sea. We observed that some Methanothermococcus lineages and their transcripts were more abundant than others in individual vent sites, indicating differential fitness among lineages. The relative abundances of lineages represented by SAGs in each of the samples matched phylogenetic relationships based on single-copy universal genes, and genes related to nitrogen fixation and the CRISPR/Cas immune system were among those differentiating the clades. Lineages possessing these genes were less abundant than those missing that genomic region. Overall, patterns in nucleotide variation indicated that the population dynamics of Methanothermococcus were not governed by clonal expansions or selective sweeps, at least in the habitats and sampling times included in this study. Together, our results show that although specific lineages of Methanothermococcus co-exist in these habitats, some outcompete others, and possession of accessory metabolic functions does not necessarily provide a fitness advantage in these habitats in all conditions. This work highlights the power of combining single-cell, metagenomic, and metatranscriptomic datasets to determine how evolution shapes microbial abundance and diversity in hydrothermal vent ecosystems.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Goordial2021,

title = {Microbial diversity and function in shallow subsurface sediment and oceanic lithosphere of the atlantis massif},

author = {Goordial, J. and D’angelo, T. and Labont\'{e}, J.M. and Poulton, N.J. and Brown, J.M. and Stepanauskas, R. and Fr\"{u}h-Green, G.L. and Orcutt, B.N.},

url = {https://journals.asm.org/doi/10.1128/mBio.00490-21},

doi = {https://doi.org/10.1128/mBio.00490-21 },

year  = {2021},

date = {2021-08-03},

journal = {mBio},

volume = {12},

number = {4},

pages = {e00490-21},

abstract = {The marine lithospheric subsurface is one of the largest biospheres on Earth; however, little is known about the identity and ecological function of microorganisms found in low abundance in this habitat, though these organisms impact global-scale biogeochemical cycling. Here, we describe the diversity and metabolic potential of sediment and endolithic (within rock) microbial communities found in ultrasmall amounts (101 to 104 cells cm−3) in the subsurface of the Atlantis Massif, an oceanic core complex on the Mid-Atlantic Ridge that was sampled on International Ocean Discovery Program (IODP) Expedition 357. This study used fluorescence-activated cell sorting (FACS) to enable the first amplicon, metagenomic, and single-cell genomic study of the shallow (<20 m below seafloor) subsurface of an actively serpentinizing marine system. The shallow subsurface biosphere of the Atlantis Massif was found to be distinct from communities observed in the nearby Lost City alkaline hydrothermal fluids and chimneys, yet similar to other low-temperature, aerobic subsurface settings. Genes associated with autotrophy were rare, although heterotrophy and aerobic carbon monoxide and formate cycling metabolisms were identified. Overall, this study reveals that the shallow subsurface of an oceanic core complex hosts a biosphere that is not fueled by active serpentinization reactions and by-products.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{\v{S}ul\v{c}ius2021,

title = {Exploring viral diversity in a gypsum karst lake ecosystem using targeted single-cell genomics},

author = {\v{S}ul\v{c}ius, S. and Alzbutas, G. and Juknevi\v{c}i\={u}t\.{e}, V. and \v{S}imoli\={u}nas, E. and Venckus, P. and \v{S}imoli\={u}nien\.{e}, M. and Pa\v{s}kauskas, R.},

url = {https://www.mdpi.com/2073-4425/12/6/886},

doi = {https://doi.org/10.3390/genes12060886},

year  = {2021},

date = {2021-06-08},

journal = {Genes},

volume = {12},

number = {6},

pages = {886},

abstract = {Little is known about the diversity and distribution of viruses infecting green sulfur bacteria (GSB) thriving in euxinic (sulfuric and anoxic) habitats, including gypsum karst lake ecosystems. In this study, we used targeted cell sorting combined with single-cell sequencing to gain insights into the gene content and genomic potential of viruses infecting sulfur-oxidizing bacteria Chlorobium clathratiforme, obtained from water samples collected during summer stratification in gypsum karst Lake Kirkilai (Lithuania). In total, 82 viral contigs were bioinformatically identified in 62 single amplified genomes (SAGs) of C. clathratiforme. The majority of viral gene and protein sequences showed little to no similarity with phage sequences in public databases, uncovering the vast diversity of previously undescribed GSB viruses. We observed a high level of lysogenization in the C. clathratiforme population, as 87% SAGs contained intact prophages. Among the thirty identified auxiliary metabolic genes (AMGs), two, thiosulfate sulfurtransferase (TST) and thioredoxin-dependent phosphoadenosine phosphosulfate (PAPS) reductase (cysH), were found to be involved in the oxidation of inorganic sulfur compounds, suggesting that viruses can influence the metabolism and cycling of this essential element. Finally, the analysis of CRISPR spacers retrieved from the consensus C. clathratiforme genome imply persistent and active virus\textendashhost interactions for several putative phages prevalent among C. clathratiforme SAGs. Overall, this study provides a glimpse into the diversity of phages associated with naturally occurring and highly abundant sulfur-oxidizing bacteria.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Murphy2021,

title = {Genomic analysis of the yet-uncultured binatota reveals broad methylotrophic, alkane-degradation, and pigment production capacities},

author = {Murphy, C.L. and Sheremet, A. and Dunfield, P.F. and Spear, J.R. and Stepanauskas, R. and Woyke, T. and Elshahed, M.S. and Youssef, N.H.},

url = {https://pubmed.ncbi.nlm.nih.gov/34006650/},

doi = {10.1128/mBio.00985-21},

year  = {2021},

date = {2021-05-18},

journal = {mBio},

volume = {12},

number = {3},

pages = {e00985-21},

abstract = {The recent leveraging of genome-resolved metagenomics has generated an enormous number of genomes from novel uncultured microbial lineages yet left many clades undescribed. Here, we present a global analysis of genomes belonging to Binatota (UBP10), a globally distributed, yet-uncharacterized bacterial phylum. All orders in Binatota encoded the capacity for aerobic methylotrophy using methanol, methylamine, sulfomethanes, and chloromethanes as the substrates. Methylotrophy in Binatota was characterized by order-specific substrate degradation preferences, as well as extensive metabolic versatility, i.e., the utilization of diverse sets of genes, pathways, and combinations to achieve a specific metabolic goal. The genomes also encoded multiple alkane hydroxylases and monooxygenases, potentially enabling growth on a wide range of alkanes and fatty acids. Pigmentation is inferred from a complete pathway for carotenoids (lycopene, β- and γ-carotenes, xanthins, chlorobactenes, and spheroidenes) production. Further, the majority of genes involved in bacteriochlorophyll a, c, and d biosynthesis were identified, although absence of key genes and failure to identify a photosynthetic reaction center preclude proposing phototrophic capacities. Analysis of 16S rRNA databases showed the preferences of Binatota to terrestrial and freshwater ecosystems, hydrocarbon-rich habitats, and sponges, supporting their potential role in mitigating methanol and methane emissions, breakdown of alkanes, and their association with sponges. Our results expand the lists of methylotrophic, aerobic alkane-degrading, and pigment-producing lineages. We also highlight the consistent encountering of incomplete biosynthetic pathways in microbial genomes, a phenomenon necessitating careful assessment when assigning putative functions based on a set-threshold of pathway completion.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{WE2021,

title = {Synthase-Selective Exploration of a Tunicate Microbiome by Activity-Guided Single-Cell Genomics},

author = {Kim WE and Charov K and Dzunkova M and Becraft ED and Brown J and Schulz F and Woyke T and La Clair JJ and Stepanauskas R and Burkart MD },

url = {https://pubs.acs.org/doi/full/10.1021/acschembio.1c00157},

year  = {2021},

date = {2021-05-06},

organization = {ACS Chem Biol},

abstract = {While thousands of environmental metagenomes have been mined for the presence of novel biosynthetic gene clusters, such computational predictions do not provide evidence of their in vivo biosynthetic functionality. Using fluorescent in situ enzyme assay targeting carrier proteins common to polyketide (PKS) and nonribosomal peptide synthetases (NRPS), we applied fluorescence-activated cell sorting to tunicate microbiome to enrich for microbes with active secondary metabolic capabilities. Single-cell genomics uncovered the genetic basis for a wide biosynthetic diversity in the enzyme-active cells and revealed a member of marine Oceanospirillales harboring a novel NRPS gene cluster with high similarity to phylogenetically distant marine and terrestrial bacteria. Interestingly, this synthase belongs to a larger class of siderophore biosynthetic gene clusters commonly associated with pestilence and disease. This demonstrates activity-guided single-cell genomics as a tool to guide novel biosynthetic discovery.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{S2021,

title = {A Genomic Catalogue of Earth’s Microbiomes},

author = {Nayfach S and Roux A and Seshadri R and Udwary D and Varghese N and Schulz F and Wu D and Paez-Espino D and Chen I-M and Huntemann M and Palaniappan K and Ladau J and Mukherjee S and Reddy TBK and Nielsen T and Kirton E and Faria EP and Edirisinghe JN and Henry CS and Jungbluth SP and Chivian D and Dehal P and Wood-Charlson EM and Arkin AP and Tringe S and Visel A and IMG/M Data Consortium and Woyke T and Mouncey NJ and Ivanova NN and Kyrpides NC and Eloe-Fadrosh EA},

doi = {https://doi.org/10.1038/s41587-020-0718-6},

year  = {2021},

date = {2021-04-12},

journal = {Nature Biotechnology},

volume = {39},

pages = {499-509},

abstract = {The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this approach to >10,000 metagenomes collected from diverse habitats covering all of Earth’s continents and oceans, including metagenomes from human and animal hosts, engineered environments, and natural and agricultural soils, to capture extant microbial, metabolic and functional potential. This comprehensive catalog includes 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units spanning 135 phyla. The catalog expands the known phylogenetic diversity of bacteria and archaea by 44% and is broadly available for streamlined comparative analyses, interactive exploration, metabolic modeling and bulk download. We demonstrate the utility of this collection for understanding secondary-metabolite biosynthetic potential and for resolving thousands of new host linkages to uncultivated viruses. This resource underscores the value of genome-centric approaches for revealing genomic properties of uncultivated microorganisms that affect ecosystem processes.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ED2021,

title = {Evolutionary stasis of a deep subsurface microbial lineage},

author = {Becraft ED and Lau MCY and Bezuidt OKI and Brown JM and Labont\'{e} JM and Kauneckaite-Griguole K and Salkauskaite R and Alzbutas G and Sackett JD and Kruger BR and Kadnikov V and Van Heerden E and Moser D and Ravin N and Onstott T and Stepanauskas R },

url = {https://www.nature.com/articles/s41396-021-00965-3#citeas},

doi = {https://doi.org/10.1038/s41396-021-00965-3},

year  = {2021},

date = {2021-04-06},

journal = {The ISME Journal},

abstract = {Sulfate-reducing bacteria Candidatus Desulforudis audaxviator (CDA) were originally discovered in deep fracture fluids accessed via South African gold mines and have since been found in geographically widespread deep subsurface locations. In order to constrain models for subsurface microbial evolution, we compared CDA genomes from Africa, North America and Eurasia using single cell genomics. Unexpectedly, 126 partial single amplified genomes from the three continents, a complete genome from of an isolate from Eurasia, and metagenome-assembled genomes from Africa and Eurasia shared >99.2% average nucleotide identity, low frequency of SNP’s, and near-perfectly conserved prophages and CRISPRs. Our analyses reject sample cross-contamination, recent natural dispersal, and unusually strong purifying selection as likely explanations for these unexpected results. We therefore conclude that the analyzed CDA populations underwent only minimal evolution since their physical separation, potentially as far back as the breakup of Pangea between 165 and 55 Ma ago. High-fidelity DNA replication and repair mechanisms are the most plausible explanation for the highly conserved genome of CDA. CDA presents a stark contrast to the current model organisms in microbial evolutionary studies, which often develop adaptive traits over far shorter periods of time.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{EM2021,

title = {Diversification of methanogens into hyperalkaline serpentinizing environments through adaptations to minimize oxidant limitation},

author = {Fones EM and Colman DR and Kraus EA and Stepanauskas R and Templeton AS and Spear JR and Boyd ES },

url = {https://www.nature.com/articles/s41396-020-00838-1#citeas},

doi = {https://doi.org/10.1038/s41396-020-00838-1},

year  = {2021},

date = {2021-04-01},

journal = {The ISME Journal},

volume = {15},

pages = {1121-1135},

abstract = {Metagenome assembled genomes (MAGs) and single amplified genomes (SAGs) affiliated with two distinct Methanobacterium lineages were recovered from subsurface fracture waters of the Samail Ophiolite, Sultanate of Oman. Lineage Type I was abundant in waters with circumneutral pH, whereas lineage Type II was abundant in hydrogen rich, hyperalkaline waters. Type I encoded proteins to couple hydrogen oxidation to CO2 reduction, typical of hydrogenotrophic methanogens. Surprisingly, Type II, which branched from the Type I lineage, lacked homologs of two key oxidative [NiFe]-hydrogenases. These functions were presumably replaced by formate dehydrogenases that oxidize formate to yield reductant and cytoplasmic CO2 via a pathway that was unique among characterized Methanobacteria, allowing cells to overcome CO2/oxidant limitation in high pH waters. This prediction was supported by microcosm-based radiotracer experiments that showed significant biological methane generation from formate, but not bicarbonate, in waters where the Type II lineage was detected in highest relative abundance. Phylogenetic analyses and variability in gene content suggested that recent and ongoing diversification of the Type II lineage was enabled by gene transfer, loss, and transposition. These data indicate that selection imposed by CO2/oxidant availability drove recent methanogen diversification into hyperalkaline waters that are heavily impacted by serpentinization.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{O2021,

title = {The cyanobacterium Prochlorococcus has divergent light-harvesting antennae and may have evolved in a low-oxygen ocean},

author = {Ulloa O and Henr\'{i}quez-Castillo C and Ram\'{i}rez-Flandes S and Plominsky AM and Murillo AA and Morgan-Lang C and Hallam SJ and Stepanauskas R },

url = {https://www.pnas.org/content/118/11/e2025638118},

doi = {https://doi.org/10.1073/pnas.2025638118},

issn = {1091-6490},

year  = {2021},

date = {2021-03-16},

journal = {PNAS},

volume = {118},

number = {11},

abstract = {Marine picocyanobacteria of the genus Prochlorococcus are the most abundant photosynthetic organisms in the modern ocean, where they exert a profound influence on elemental cycling and energy flow. The use of transmembrane chlorophyll complexes instead of phycobilisomes as light-harvesting antennae is considered a defining attribute of Prochlorococcus. Its ecology and evolution are understood in terms of light, temperature, and nutrients. Here, we report single-cell genomic information on previously uncharacterized phylogenetic lineages of this genus from nutrient-rich anoxic waters of the eastern tropical North and South Pacific Ocean. The most basal lineages exhibit optical and genotypic properties of phycobilisome-containing cyanobacteria, indicating that the characteristic light-harvesting antenna of the group is not an ancestral attribute. Additionally, we found that all the indigenous lineages analyzed encode genes for pigment biosynthesis under oxygen-limited conditions, a trait shared with other freshwater and coastal marine cyanobacteria. Our findings thus suggest that Prochlorococcus diverged from other cyanobacteria under low-oxygen conditions before transitioning from phycobilisomes to transmembrane chlorophyll complexes and may have contributed to the oxidation of the ancient ocean.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JM2020,

title = {Single Cell Genomics Reveals Viruses Consumed by Marine Protists},

author = {Brown JM and Labont\'{e} JM and Brown J and Record NR and Poulton NJ and Sieracki ME and Logares R and Stepanauskas R},

url = {https://www.frontiersin.org/articles/10.3389/fmicb.2020.524828/full},

doi = {https://doi.org/10.3389/fmicb.2020.524828},

year  = {2020},

date = {2020-09-24},

journal = {Frontiers in Microbiology},

volume = {11},

number = {524828},

abstract = {The predominant model of the role of viruses in the marine trophic web is that of the “viral shunt,” where viral infection funnels a substantial fraction of the microbial primary and secondary production back to the pool of dissolved organic matter. Here, we analyzed the composition of non-eukaryotic DNA associated with individual cells of small, planktonic protists in the Gulf of Maine (GoM) and the Mediterranean Sea. We found viral DNA associated with a substantial fraction cells from the GoM (51%) and the Mediterranean Sea (35%). While Mediterranean SAGs contained a larger proportion of cells containing bacterial sequences (49%), a smaller fraction of cells contained bacterial sequences in the GoM (19%). In GoM cells, nearly identical bacteriophage and ssDNA virus sequences where found across diverse lineages of protists, suggesting many of these viruses are non-infective. The fraction of cells containing viral DNA varied among protistan lineages and reached 100% in Picozoa and Choanozoa. These two groups also contained significantly higher numbers of viral sequences than other identified taxa. We consider mechanisms that may explain the presence of viral DNA in protistan cells and conclude that protistan predation on free viral particles contributed to the observed patterns. These findings confirm prior experiments with protistan isolates and indicate that the viral shunt is complemented by a viral link in the marine microbial food web. This link may constitute a sink of viral particles in the ocean and has implications for the flow of carbon through the microbial food web.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JP2020,

title = {Ancestral absence of electron transport chains in Patescibacteria and DPANN},

author = {Beam JP and Becraft ED and Brown JM and Schulz F and Jarett JK and Bezuidt O and Poulton NJ and Clark K and Dunfield PF and Ravin NV and Spear JR and Hedlund BP and Kormas KA and Sievert SM and Elshahed MS and Barton HA and Stott MB and Eisen JA and Moser DP and Onstott TC and Woyke T and Stepanauskas R},

url = {https://www.frontiersin.org/articles/10.3389/fmicb.2020.01848/full},

doi = {https://doi.org/10.3389/fmicb.2020.01848},

year  = {2020},

date = {2020-08-17},

journal = {Frontiers in Microbiology},

volume = {11},

number = {1848},

abstract = {Recent discoveries suggest that the candidate superphyla Patescibacteria and DPANN constitute a large fraction of the phylogenetic diversity of Bacteria and Archaea. Their small genomes and limited coding potential have been hypothesized to be ancestral adaptations to obligate symbiotic lifestyles. To test this hypothesis, we performed cell\textendashcell association, genomic, and phylogenetic analyses on 4,829 individual cells of Bacteria and Archaea from 46 globally distributed surface and subsurface field samples. This confirmed the ubiquity and abundance of Patescibacteria and DPANN in subsurface environments, the small size of their genomes and cells, and the divergence of their gene content from other Bacteria and Archaea. Our analyses suggest that most Patescibacteria and DPANN in the studied subsurface environments do not form specific physical associations with other microorganisms. These data also suggest that their unusual genomic features and prevalent auxotrophies may be a result of ancestral, minimal cellular energy transduction mechanisms that lack respiration, thus relying solely on fermentation for energy conservation.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JK2020,

title = {Insights into the dynamics between viruses and their hosts in a hot spring microbial mat},

author = {Jarett JK and D\v{z}unkov\'{a} M and Schulz F and Roux S and Paez-Espino D and Eloe-Fadrosh E and Jungbluth SP and Ivanova N and Spear JR and Carr SA and Trivedi CB and Corsetti FA and Johnson HA and Becraft E and Kyrpides N and Stepanauskas R and Woyke T},

url = {https://www.nature.com/articles/s41396-020-0705-4#citeas},

doi = {https://doi.org/10.1038/s41396-020-0705-4},

year  = {2020},

date = {2020-07-07},

journal = {The ISME Journal},

volume = {14},

pages = {2527-2541},

abstract = {Our current knowledge of host\textendashvirus 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\textendashvirus 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\textendashhost 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\textendashvirus 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.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ML2020,

title = {Hiding in Plain Sight: The Globally Distributed Bacterial Candidate Phylum PAUC34f},

author = {Chen ML and Becraft ED and Pachiadaki M and Brown JM and Jarett JK and Gasol JM and Ravin NV and Moser DP and Nunoura T and Herndl GJ and Woyke T and Stepanauskas R },

url = {https://www.frontiersin.org/articles/10.3389/fmicb.2020.00376/full},

doi = {https://doi.org/10.3389/fmicb.2020.00376},

year  = {2020},

date = {2020-03-12},

journal = {Frontiers in Microbiology},

volume = {11},

number = {376},

abstract = {Bacterial candidate phylum PAUC34f was originally discovered in marine sponges and is widely considered to be composed of sponge symbionts. Here, we report 21 single amplified genomes (SAGs) of PAUC34f from a variety of environments, including the dark ocean, lake sediments, and a terrestrial aquifer. The diverse origins of the SAGs and the results of metagenome fragment recruitment suggest that some PAUC34f lineages represent relatively abundant, free-living cells in environments other than sponge microbiomes, including the deep ocean. Both phylogenetic and biogeographic patterns, as well as genome content analyses suggest that PAUC34f associations with hosts evolved independently multiple times, while free-living lineages of PAUC34f are distinct and relatively abundant in a wide range of environments.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JHW2020,

title = {Pangenomics Analysis Reveals Diversification of Enzyme Families and Niche Specialization in Globally Abundant SAR202 Bacteria},

author = {Saw JHW and Nunoura T and Hirai M and Takaki Y and Parsons R and Michelsen M and Longnecker K and Kujawinski EB and Stepanauskas R and Landry Z and Carlson CA and Giovannoni SJ },

url = {https://mbio.asm.org/content/11/1/e02975-19},

doi = {10.1128/mBio.02975-19},

year  = {2020},

date = {2020-01-07},

journal = {mBio},

volume = {11},

number = {1},

pages = {e02975-19},

abstract = {It has been hypothesized that the abundant heterotrophic ocean bacterioplankton in the SAR202 clade of the phylum Chloroflexi evolved specialized metabolisms for the oxidation of organic compounds that are resistant to microbial degradation via common metabolic pathways. Expansions of paralogous enzymes were reported and implicated in hypothetical metabolism involving monooxygenase and dioxygenase enzymes. In the proposed metabolic schemes, the paralogs serve the purpose of diversifying the range of organic molecules that cells can utilize. To further explore SAR202 evolution and metabolism, we reconstructed single amplified genomes and metagenome-assembled genomes from locations around the world that included the deepest ocean trenches. In an analysis of 122 SAR202 genomes that included seven subclades spanning SAR202 diversity, we observed additional evidence of paralog expansions that correlated with evolutionary history, as well as further evidence of metabolic specialization. Consistent with previous reports, families of flavin-dependent monooxygenases were observed mainly in the group III SAR202 genomes, and expansions of dioxygenase enzymes were prevalent in those of group VII. We found that group I SAR202 genomes encode expansions of racemases in the enolase superfamily, which we propose evolved for the degradation of compounds that resist biological oxidation because of chiral complexity. Supporting the conclusion that the paralog expansions indicate metabolic specialization, fragment recruitment and fluorescent in situ hybridization (FISH) with phylogenetic probes showed that SAR202 subclades are indigenous to different ocean depths and geographical regions. Surprisingly, some of the subclades were abundant in surface waters and contained rhodopsin genes, altering our understanding of the ecological role of SAR202 species in stratified water columns.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MG2019,

title = {Charting the complexity of the marine microbiome through single cell genomics},

author = {Pachiadaki MG and Brown JM and Brown J and Bezuidt O and Berube PM and Biller SJ and Poulton NJ and Burkart MD and La Clair JJ and Chisholm SW and Stepanauskas R},

url = {https://reader.elsevier.com/reader/sd/pii/S0092867419312735?token=5CCA0F9F7B2364A01F02747FA24E677BBD690529309015082F4BD56865547E22999EEFA8E287D658003474FFC76566E3},

doi = {https://doi.org/10.1016/j.cell.2019.11.017},

year  = {2019},

date = {2019-12-12},

journal = {Cell},

volume = {179},

number = {7},

pages = {1623-1635},

abstract = {Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JM2019,

title = {Single cell genomics-based analysis of gene content and expression of prophages in a diffuse-flow deep-sea hydrothermal system},

author = {Labont\'{e} JM and Pachiadaki M and Fergusson E and McNichol J and Grosche A and Gulmann LK and Vetriani C and Sievert SM and Stepanauskas R},

url = {https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6581674/},

doi = {10.3389/fmicb.2019.01262},

year  = {2019},

date = {2019-06-12},

journal = {Frontiers in Microbiology},

volume = {10},

pages = {1262},

abstract = {Phage\textendashhost interactions likely play a major role in the composition and functioning of many microbiomes, yet remain poorly understood. Here, we employed single cell genomics to investigate phage\textendashhost interactions in a diffuse-flow, low-temperature hydrothermal vent that may be reflective of a broadly distributed biosphere in the subseafloor. We identified putative prophages in 13 of 126 sequenced single amplified genomes (SAGs), with no evidence for lytic infections, which is in stark contrast to findings in the surface ocean. Most were distantly related to known prophages, while their hosts included bacterial phyla Campylobacterota, Bacteroidetes, Chlorobi, Proteobacteria, Lentisphaerae, Spirochaetes, and Thermotogae. Our results suggest the predominance of lysogeny over lytic interaction in diffuse-flow, deep-sea hydrothermal vents, despite the high activity of the dominant Campylobacteria that would favor lytic infections. We show that some of the identified lysogens have co-evolved with their host over geological time scales and that their genes are transcribed in the environment. Functional annotations of lysogeny-related genes suggest involvement in horizontal gene transfer enabling host’s protection against toxic metals and antibacterial compounds.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{LR2019,

title = {Red Sea SAR11 and Prochlorococcus Single-Cell Genomes Reflect Globally Distributed Pangenomes},

author = {Thompson LR and Haroon MF and Shibl AA and Cahill MJ and Ngugi DK and Williams GJ and Morton JT and Knight R and Goodwin KD and Stingl U},

url = {https://aem.asm.org/content/85/13/e00369-19},

doi = {10.1128/AEM.00369-19},

year  = {2019},

date = {2019-04-26},

journal = {Applied and Environmental Microbiology},

volume = {85},

abstract = {Evidence suggests many marine bacteria are cosmopolitan, with widespread but sparse strains poised to seed abundant populations under conducive growth conditions. However, studies supporting this “microbial seed bank” hypothesis have analyzed taxonomic marker genes rather than whole genomes/metagenomes, leaving open the possibility that disparate ocean regions harbor endemic gene content. The Red Sea is isolated geographically from the rest of the ocean and has a combination of high irradiance, high temperature, and high salinity that is unique among the oceans; we therefore asked whether it harbors endemic gene content. We sequenced and assembled single-cell genomes of 21 SAR11 (subclades Ia, Ib, Id, and II) and 5 Prochlorococcus (ecotype HLII) samples from the Red Sea and combined them with globally sourced reference genomes to cluster genes into ortholog groups (OGs). Ordination of OG composition could distinguish clades, including phylogenetically cryptic Prochlorococcus ecotypes LLII and LLIII. Compared with reference genomes, 1% of Prochlorococcus and 17% of SAR11 OGs were unique to the Red Sea genomes (RS-OGs). Most (83%) RS-OGs had no annotated function, but 65% of RS-OGs were expressed in diel Red Sea metatranscriptomes, suggesting they are functional. Searching Tara Oceans metagenomes, RS-OGs were as likely to be found as non-RS-OGs; nevertheless, Red Sea and other warm samples could be distinguished from cooler samples using the relative abundances of OGs. The results suggest that the prevalence of OGs in these surface ocean bacteria is largely cosmopolitan, with differences in population metagenomes manifested by differences in relative abundance rather than complete presence/absence of OGs.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{M2019b,

title = {Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution},

author = {Ren M and Feng X and Huang Y and Wang H and Hu Z and Clingenpeel S and Swan BK and Fonseca MM and Posada D and Stepanauskas R and Hollibaugh JT and Foster PG and Woyke T and Luo H},

url = {https://www.nature.com/articles/s41396-019-0418-8},

doi = {https://doi.org/10.1038/s41396-019-0418-8},

year  = {2019},

date = {2019-04-25},

journal = {ISME},

volume = {13},

pages = {2150-2161},

abstract = {Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in marine and terrestrial habitats, playing a major role in the global nitrogen cycle. However, their evolutionary history remains unexplored, which limits our understanding of their adaptation mechanisms. Here, our comprehensive phylogenomic tree of Thaumarchaeota supports three sequential events: origin of AOA from terrestrial non-AOA ancestors, colonization of the shallow ocean, and expansion to the deep ocean. Careful molecular dating suggests that these events coincided with the Great Oxygenation Event around 2300 million years ago (Mya), and oxygenation of the shallow and deep ocean around 800 and 635\textendash560 Mya, respectively. The first transition was likely enabled by the gain of an aerobic pathway for energy production by ammonia oxidation and biosynthetic pathways for cobalamin and biotin that act as cofactors in aerobic metabolism. The first transition was also accompanied by the loss of dissimilatory nitrate and sulfate reduction, loss of oxygen-sensitive pyruvate oxidoreductase, which reduces pyruvate to acetyl-CoA, and loss of the Wood\textendashLjungdahl pathway for anaerobic carbon fixation. The second transition involved gain of a K+ transporter and of the biosynthetic pathway for ectoine, which may function as an osmoprotectant. The third transition was accompanied by the loss of the uvr system for repairing ultraviolet light-induced DNA lesions. We conclude that oxygen availability drove the terrestrial origin of AOA and their expansion to the photic and dark oceans, and that the stressors encountered during these events were partially overcome by gene acquisitions from Euryarchaeota and Bacteria, among other sources.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ME2019,

title = {Single cell genomics yields a wide diversity of small planktonic protists across major ocean ecosystems},

author = {Sieracki ME and Poulton NJ and Jaillon O and Wincker P and de Vargas C and Rubinat-Ripoll L and Stepanauskas R and Logares R and Massana R},

url = {https://www.nature.com/articles/s41598-019-42487-1},

doi = {https://doi.org/10.1038/s41598-019-42487-1},

year  = {2019},

date = {2019-04-15},

journal = {Scientific Reports},

volume = {9},

pages = {6025},

abstract = {Marine planktonic protists are critical components of ocean ecosystems and are highly diverse. Molecular sequencing methods are being used to describe this diversity and reveal new associations and metabolisms that are important to how these ecosystems function. We describe here the use of the single cell genomics approach to sample and interrogate the diversity of the smaller (pico- and nano-sized) protists from a range of oceanic samples. We created over 900 single amplified genomes (SAGs) from 8 Tara Ocean samples across the Indian Ocean and the Mediterranean Sea. We show that flow cytometric sorting of single cells effectively distinguishes plastidic and aplastidic cell types that agree with our understanding of protist phylogeny. Yields of genomic DNA with PCR-identifiable 18S rRNA gene sequence from single cells was low (15% of aplastidic cell sorts, and 7% of plastidic sorts) and tests with alternate primers and comparisons to metabarcoding did not reveal phylogenetic bias in the major protist groups. There was little evidence of significant bias against or in favor of any phylogenetic group expected or known to be present. The four open ocean stations in the Indian Ocean had similar communities, despite ranging from 14°N to 20°S latitude, and they differed from the Mediterranean station. Single cell genomics of protists suggests that the taxonomic diversity of the dominant taxa found in only several hundreds of microliters of surface seawater is similar to that found in molecular surveys where liters of sample are filtered.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{NH2019,

title = {Genomic characterization of candidate division LCP-89 reveals an atypical cell wall structure, microcompartment production, and dual respiratory and fermentative capacities},

author = {Youssef NH and Farag IF and Hahn CR and Jarett J and Becraft E and Eloe-Fadrosh E and Lightfoot J and Bourgeois A and Cole T and Ferrante S and Truelock M and Marsh W and Jamaleddine M and Ricketts S and Simpson R and McFadden A and Hoff W and Ravin NV and Sievert S and Stepanauskas R and Woyke T and Elshahed M},

url = {https://aem.asm.org/content/85/10/e00110-19},

doi = {10.1128/AEM.00110-19},

year  = {2019},

date = {2019-03-22},

journal = {Applied and Environmental Microbiology},

volume = {85},

pages = {e00110-19},

abstract = {Recent experimental and bioinformatic advances enable the recovery of genomes belonging to yet-uncultured microbial lineages directly from environmental samples. Here, we report on the recovery and characterization of single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs) representing candidate phylum LCP-89, previously defined based on 16S rRNA gene sequences. Analysis of LCP-89 genomes recovered from Zodletone Spring, an anoxic spring in Oklahoma, predicts slow-growing, rod-shaped organisms. LCP-89 genomes contain genes for cell wall lipopolysaccharide (LPS) production but lack the entire machinery for peptidoglycan biosynthesis, suggesting an atypical cell wall structure. The genomes, however, encode S-layer homology domain-containing proteins, as well as machinery for the biosynthesis of CMP-legionaminate, inferring the possession of an S-layer glycoprotein. A nearly complete chemotaxis machinery coupled to the absence of flagellar synthesis and assembly genes argues for the utilization of alternative types of motility. A strict anaerobic lifestyle is predicted, with dual respiratory (nitrite ammonification) and fermentative capacities. Predicted substrates include a wide range of sugars and sugar alcohols and a few amino acids. The capability of rhamnose metabolism is confirmed by the identification of bacterial microcompartment genes to sequester the toxic intermediates generated. Comparative genomic analysis identified differences in oxygen sensitivities, respiratory capabilities, substrate utilization preferences, and fermentation end products between LCP-89 genomes and those belonging to its four sister phyla (Calditrichota, SM32-31, AABM5-125-24, and KSB1) within the broader FCB (Fibrobacteres-Chlorobi-Bacteroidetes) superphylum. Our results provide a detailed characterization of members of the candidate division LCP-89 and highlight the importance of reconciling 16S rRNA-based and genome-based phylogenies.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JD2019,

title = {Four draft single-cell genome sequences of novel, nearly identical Kiritimatiellaeota strains isolated from the continental deep subsurface},

author = {Sackett JD and Kruger BR and Becraft ED and Jarett JK and Stepanauskas R and Woyke T and Moser DP},

url = {https://mra.asm.org/content/8/11/e01249-18},

doi = {10.1128/MRA.01249-18},

year  = {2019},

date = {2019-03-14},

journal = {Microbiology Resource Announcements},

volume = {8},

pages = {e01249-18},

abstract = {The recently proposed bacterial phylum Kiritimatiellaeota represents a globally distributed monophyletic clade distinct from other members of the Planctomycetes, Verrucomicrobia, and Chlamydiae (PVC) superphylum. Here, we present four phylogenetically distinct single-cell genome sequences from within the Kiritimatiellaeota lineage sampled from deep continental subsurface aquifer fluids of the Death Valley Regional Flow System in the United States.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{SA2019,

title = {Carboxydotrophy potential of uncultivated Hydrothermarchaeota from the subseafloor crustal biosphere},

author = {Carr SA and Jungbluth SP and Eloe-Fadrosh EA and Stepanauskas R and Woyke T and Rapp\'{e} MS and Orcutt BN},

url = {https://www.nature.com/articles/s41396-019-0352-9},

doi = {https://doi.org/10.1038/s41396-019-0352-9},

issn = {1751-7370},

year  = {2019},

date = {2019-02-07},

journal = {The ISME Journal},

volume = {13},

pages = {1457-1468},

abstract = {The exploration of Earth’s terrestrial subsurface biosphere has led to the discovery of several new archaeal lineages of evolutionary significance. Similarly, the deep subseafloor crustal biosphere also harbors many unique, uncultured archaeal taxa, including those belonging to Candidatus Hydrothermarchaeota, formerly known as Marine Benthic Group-E. Recently, Hydrothermarchaeota was identified as an abundant lineage of Juan de Fuca Ridge flank crustal fluids, suggesting its adaptation to this extreme environment. Through the investigation of single-cell and metagenome-assembled genomes, we provide insight into the lineage’s evolutionary history and metabolic potential. Phylogenomic analysis reveals the Hydrothermarchaeota to be an early-branching archaeal phylum, branching between the superphylum DPANN, Euryarchaeota, and Asgard lineages. Hydrothermarchaeota genomes suggest a potential for dissimilative and assimilative carbon monoxide oxidation (carboxydotrophy), as well as sulfate and nitrate reduction. There is also a prevalence of chemotaxis and motility genes, indicating adaptive strategies for this nutrient-limited fluid-rock environment. These findings provide the first genomic interpretations of the Hydrothermarchaeota phylum and highlight the anoxic, hot, deep marine crustal biosphere as an important habitat for understanding the evolution of early life.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{PM2019,

title = {Emergence of trait variability through the lens of nitrogen assimilation in Prochlorococcus},

author = {Berube PM and Rasmussen A and Braakman R and Stepanauskas R and Chisholm SW},

url = {https://elifesciences.org/articles/41043},

doi = {10.7554/eLife.41043},

year  = {2019},

date = {2019-02-01},

journal = {eLife},

volume = {8},

abstract = {Intraspecific trait variability has important consequences for the function and stability of marine ecosystems. Here we examine variation in the ability to use nitrate across hundreds of Prochlorococcus genomes to better understand the modes of evolution influencing intraspecific allocation of ecologically important functions. Nitrate assimilation genes are absent in basal lineages but occur at an intermediate frequency that is randomly distributed within recently emerged clades. The distribution of nitrate assimilation genes within clades appears largely governed by vertical inheritance, gene loss, and homologous recombination. By mapping this process onto a model of Prochlorococcus’ macroevolution, we propose that niche-constructing adaptive radiations and subsequent niche partitioning set the stage for loss of nitrate assimilation genes from basal lineages as they specialized to lower light levels. Retention of these genes in recently emerged lineages has likely been facilitated by selection as they sequentially partitioned into niches where nitrate assimilation conferred a fitness benefit.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Y2019,

title = {Genomics insights into ecotype formation of ammonia-oxidizing archaea in the deep ocean},

author = {Wang Y and Huang JM and Cui GJ and Nunoura T and Takaki Y and Li WL and Li J and Gao ZM and Takai K and Zhang AQ and Stepanauskas R},

url = {https://onlinelibrary.wiley.com/doi/abs/10.1111/1462-2920.14518},

doi = {https://doi.org/10.1111/1462-2920.14518},

year  = {2019},

date = {2019-01-31},

journal = {Environmental Microbiology },

volume = {21},

number = {2},

pages = {716-729},

abstract = {Various lineages of ammonia‐oxidizing archaea (AOA) are present in deep waters, but the mechanisms that determine ecotype formation are obscure. We studied 18 high‐quality genomes of the marine group I AOA lineages (alpha, gamma and delta) from the Mariana and Ogasawara trenches. The genomes of alpha AOA resembled each other, while those of gamma and delta lineages were more divergent and had even undergone insertion of some phage genes. The instability of the gamma and delta AOA genomes could be partially due to the loss of DNA polymerase B (polB) and methyladenine DNA glycosylase (tag) genes responsible for the repair of point mutations. The alpha AOA genomes harbour genes encoding a thrombospondin‐like outer membrane structure that probably serves as a barrier to gene flow. Moreover, the gamma and alpha AOA lineages rely on vitamin B12‐independent MetE and B12‐dependent MetH, respectively, for methionine synthesis. The delta AOA genome contains genes involved in uptake of sugar and peptide perhaps for heterotrophic lifestyle. Our study provides insights into co‐occurrence of cladogenesis and anagenesis in the formation of AOA ecotypes that perform differently in nitrogen and carbon cycling in dark oceans.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{PB2019,

title = {Hydrogen-based metabolism as an ancestral trait in lineages sibling to the Cyanobacteria},

author = {Matheus Carnevali PB and Schulz F and Castelle CJ and Kantor RS and Shih PM and Sharon I and Santini JM and Olm MR and Amano Y and Thomas BC and Anantharaman K and Burstein D and Becraft ED and Stepanauskas R and Woyke T and Banfield JF},

url = {https://www.nature.com/articles/s41467-018-08246-y},

doi = {https://doi.org/10.1038/s41467-018-08246-y},

issn = {2041-1723},

year  = {2019},

date = {2019-01-28},

journal = {Nature Communications},

volume = {10},

number = {463},

abstract = {The evolution of aerobic respiration was likely linked to the origins of oxygenic Cyanobacteria. Close phylogenetic neighbors to Cyanobacteria, such as Margulisbacteria (RBX-1 and ZB3), Saganbacteria (WOR-1), Melainabacteria and Sericytochromatia, may constrain the metabolic platform in which aerobic respiration arose. Here, we analyze genomic sequences and predict that sediment-associated Margulisbacteria have a fermentation-based metabolism featuring a variety of hydrogenases, a streamlined nitrogenase, and electron bifurcating complexes involved in cycling of reducing equivalents. The genomes of ocean-associated Margulisbacteria encode an electron transport chain that may support aerobic growth. Some Saganbacteria genomes encode various hydrogenases, and others may be able to use O2 under certain conditions via a putative novel type of heme copper O2 reductase. Similarly, Melainabacteria have diverse energy metabolisms and are capable of fermentation and aerobic or anaerobic respiration. The ancestor of all these groups may have been an anaerobe in which fermentation and H2 metabolism were central metabolic features. The ability to use O2 as a terminal electron acceptor must have been subsequently acquired by these lineages.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{YM2018,

title = {Genomic insight into the predominance of candidate phylum Atribacteria JS1 lineage in marine sediments},

author = {Lee YM and Hwang K and Lee JI and Kim M and Hwang CY and Noh HJ and Choi H and Lee HK and Chun J and Hong SG and Shin SC },

url = {https://www.frontiersin.org/articles/10.3389/fmicb.2018.02909/full},

doi = { 10.3389/fmicb.2018.02909},

year  = {2018},

date = {2018-11-29},

journal = {Frontiers in Microbiology},

volume = {9},

abstract = {Candidate phylum Atribacteria JS1 lineage is one of the predominant bacterial groups in anoxic subseafloor sediments, especially in organic-rich or gas hydrate-containing sediments. However, due to the lack of axenic culture representatives, metabolic potential and biogeochemical roles of this phylum have remained elusive. Here, we examined the microbial communities of marine sediments of the Ross Sea, Antarctica, and found candidate phylum Atribacteria JS1 lineage was the most abundant candidate phylum accounting for 9.8\textendash40.8% of the bacterial communities with a single dominant operational taxonomic unit (OTU). To elucidate the metabolic potential and ecological function of this species, we applied a single-cell genomic approach and obtained 18 single-cell amplified genomes presumably from a single species that was consistent with the dominant OTU throughout the sediments. The composite genome constructed by co-assembly showed the highest genome completeness among available Atribacteria JS1 genomes. Metabolic reconstruction suggested fermentative potential using various substrates and syntrophic acetate oxidation coupled with hydrogen or formate scavenging methanogens. This metabolic potential supports the predominance of Atribacteria JS1 in anoxic environments expanding our knowledge of the ecological function of this uncultivated group.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JK2018,

title = {Single-cell genomics of co-sorted Nanoarchaeota suggests novel putative host associations and diversification of proteins involved in symbiosis},

author = {Jarett JK and Nayfach S and Podar M and Inskeep W and Ivanova NN and Munson-Mcgee J and Schulz F and Young M and Jay ZJ and Beam JP and Kyrpides NC and Malmstrom RR and Stepanauskas R and Woyke T},

url = {https://microbiomejournal.biomedcentral.com/articles/10.1186/s40168-018-0539-8},

doi = {https://doi.org/10.1186/s40168-018-0539-8},

issn = {2049-2618},

year  = {2018},

date = {2018-09-17},

journal = {Microbiome},

volume = {6},

number = {161},

abstract = {Nanoarchaeota are obligate symbionts of other Archaea first discovered 16 years ago, yet little is known about this largely uncultivated taxon. While Nanoarchaeota diversity has been detected in a variety of habitats using 16S rRNA gene surveys, genome sequences have been available for only three Nanoarchaeota and their hosts. The host range and adaptation of Nanoarchaeota to a wide range of environmental conditions has thus largely remained elusive. Single-cell genomics is an ideal approach to address these questions as Nanoarchaeota can be isolated while still attached to putative hosts, enabling the exploration of cell-cell interactions and fine-scale genomic diversity.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{PM2018,

title = {Single cell genomes of Prochlorococcus, Synechococcus, and sympatric microbes from diverse marine environments},

author = {Berube PM and Biller SJ and Hackl T and Hogle SL and Satinsky BM and Becker JW and Braakman R and Collins SB and Kelly L and Berta-Thompson J and Coe A and Bergauer K and Bouman HA and Browning TJ and De Corte D and Hassler C and Hulata Y and Jacquot JE and Maas EW and Reinthaler T and Sintes E and Yokokawa T and Lindell D and Stepanauskas R and Chisholm SW},

url = {https://www.nature.com/articles/sdata2018154},

doi = {https://doi.org/10.1038/sdata.2018.154},

issn = {2052-4463},

year  = {2018},

date = {2018-09-04},

journal = {Scientific Data},

volume = {5},

abstract = {Prochlorococcus and Synechococcus are the dominant primary producers in marine ecosystems and perform a significant fraction of ocean carbon fixation. These cyanobacteria interact with a diverse microbial community that coexists with them. Comparative genomics of cultivated isolates has helped address questions regarding patterns of evolution and diversity among microbes, but the fraction that can be cultivated is miniscule compared to the diversity in the wild. To further probe the diversity of these groups and extend the utility of reference sequence databases, we report a data set of single cell genomes for 489 Prochlorococcus, 50 Synechococcus, 9 extracellular virus particles, and 190 additional microorganisms from a diverse range of bacterial, archaeal, and viral groups. Many of these uncultivated single cell genomes are derived from samples obtained on GEOTRACES cruises and at well-studied oceanographic stations, each with extensive suites of physical, chemical, and biological measurements. The genomic data reported here greatly increases the number of available Prochlorococcus genomes and will facilitate studies on evolutionary biology, microbial ecology, and biological oceanography.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{DK2018,

title = {High-quality draft single-cell genome sequence of the NS5 marine group from the coastal Red Sea},

author = {Ngugi DK and Stingl U },

url = {https://www.ncbi.nlm.nih.gov/pubmed/29930069},

doi = {10.1128/genomeA.00565-18},

year  = {2018},

date = {2018-06-21},

journal = {Genome Announcements},

volume = {6},

number = {25},

abstract = {The uncultured NS5 marine group represents one of the most ubiquitous flavobacterial bacterioplankton associated with marine blooms in the pelagic ocean. Here, we present a single-cell genome sampled from coastal waters in the Red Sea that represents the first high-quality draft genome sequence within the NS5 lineage.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{M2018,

title = {Metabolic potential and in situ transcriptomic profiles of previously uncharacterized key microbial groups involved in coupled carbon, nitrogen and sulfur cycling in anoxic marine zones},

author = {Plominsky AM and Trefault N and Podell S and Blanton JM and De la Iglesia R and Allen EE and von Dassow P and Ulloa O},

url = {https://www.ncbi.nlm.nih.gov/pubmed/29575531},

doi = {10.1111/1462-2920.14109},

year  = {2018},

date = {2018-03-24},

journal = {Enviornmental Microbiology},

abstract = {Anoxic marine zones (AMZs) impact biogeochemical cycles at the global scale, particularly the nitrogen cycle. Key microbial players from AMZs have been identified, but the majority remains unrecognized or uncharacterized. Thirty-one single-cell amplified genomes (SAGs) from the eastern tropical North and South Pacific AMZs were sequenced to gain insight into the distribution, metabolic potential and contribution to the community transcriptional profile of these uncharacterized bacterial and archaeal groups. Detailed analyses focused on SAG-bins assigned to three of these groups that presented 79%-100% estimated genome completeness: the putative sulphur-oxidizing Gamaproteobacteria EOSA II clade, a Marinimicrobia member of the recently recognized PN262000N21 clade found to be abundant in AMZ anoxic cores, and a representative of the Marine Benthic Group A Thaumarchaeota. Community-based analyses revealed that these three groups are significantly more abundant and transcriptionally more active in the AMZ microbial communities than previously described phylogenetically related microbial groups. Collectively, these groups have the potential to link biogeochemically relevant processes by coupling the carbon, nitrogen and sulfur cycles. Together, these results increase our understanding of key microbial components inhabiting AMZs and other oxygen-deficient marine environments, enhancing our capacity to predict the impact of the expansion of these ecosystems due to climate change.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{SL2017,

title = {Contrasting patterns of genome-level diversity across distinct co-occurring bacterial populations},

author = {Garcia SL and Stevens SLR and Crary B and Martinez-Garcia M and Stepanauskas R and Woyke T and Tringe SG and Andersson SGE and Bertilsson S and Malmstrom RR and McMahon KD},

url = {https://www.nature.com/articles/s41396-017-0001-0},

doi = {10.1038/s41396-017-0001-0},

year  = {2018},

date = {2018-03-01},

journal = {The ISME Journal},

volume = {12},

number = {3},

pages = {742-755},

abstract = {To understand the forces driving differentiation and diversification in wild bacterial populations, we must be able to delineate and track ecologically relevant units through space and time. Mapping metagenomic sequences to reference genomes derived from the same environment can reveal genetic heterogeneity within populations, and in some cases, be used to identify boundaries between genetically similar, but ecologically distinct, populations. Here we examine population-level heterogeneity within abundant and ubiquitous freshwater bacterial groups such as the acI Actinobacteria and LD12 Alphaproteobacteria (the freshwater sister clade to the marine SAR11) using 33 single-cell genomes and a 5-year metagenomic time series. The single-cell genomes grouped into 15 monophyletic clusters (termed "tribes") that share at least 97.9% 16S rRNA identity. Distinct populations were identified within most tribes based on the patterns of metagenomic read recruitments to single-cell genomes representing these tribes. Genetically distinct populations within tribes of the acI Actinobacterial lineage living in the same lake had different seasonal abundance patterns, suggesting these populations were also ecologically distinct. In contrast, sympatric LD12 populations were less genetically differentiated. This suggests that within one lake, some freshwater lineages harbor genetically discrete (but still closely related) and ecologically distinct populations, while other lineages are composed of less differentiated populations with overlapping niches. Our results point at an interplay of evolutionary and ecological forces acting on these communities that can be observed in real time.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JH2018,

title = {A virus or more in (nearly) every cell: ubiquitous networks of virus\textendashhost interactions in extreme environments},

author = {Munson-McGee JH and Peng S and Dewerff S and Stepanauskas R and Whitaker RJ and Weitz JS and Young MJ},

url = {https://www.nature.com/articles/s41396-018-0071-7},

doi = {10.1038/s41396-018-0071-7},

year  = {2018},

date = {2018-02-21},

journal = {The ISME Journal},

volume = {12},

number = {7},

pages = {1706-1714},

abstract = {The application of viral and cellular metagenomics to natural environments has expanded our understanding of the structure, functioning, and diversity of microbial and viral communities. The high diversity of many communities, e.g., soils, surface ocean waters, and animal-associated microbiomes, make it difficult to establish virus-host associations at the single cell (rather than population) level, assign cellular hosts, or determine the extent of viral host range from metagenomics studies alone. Here, we combine single-cell sequencing with environmental metagenomics to characterize the structure of virus\textendashhost associations in a Yellowstone National Park (YNP) hot spring microbial community. Leveraging the relatively low diversity of the YNP environment, we are able to overlay evidence at the single-cell level with contextualized viral and cellular community structure. Combining evidence from hexanucelotide analysis, single cell read mapping, network-based analytics, and CRISPR-based inference, we conservatively estimate that >60% of cells contain at least one virus type and a majority of these cells contain two or more virus types. Of the detected virus types, nearly 50% were found in more than 2 cellular clades, indicative of a broad host range. The new lens provided by the combination of metaviromics and single-cell genomics reveals a network of virus\textendashhost interactions in extreme environments, provides evidence that extensive virus\textendashhost associations are common, and further expands the unseen impact of viruses on cellular life.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Y2018,

title = {Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans},

author = {Seeleuthner Y and Mondy S and Lombard V and Carradec Q and Pelletier E and Wessner M and Leconte J and Mangot JF and Poulain J and Labadie K and Logares R and Sunagawa S and De Berardinis V and Salanoubat M and Dimier C and Kandels-Lewis S and Picheral M and Searson S and Acinas SG and Boss E and Follows M and Gorsky G and Grimsley N and Karp-Boss L and Krzic U and Not F and Ogata H and Raes J and Reynaud EG and Sardet C and Speich S and Stemmann L and Velayoudon D and Weissenbach J and Pesant S and Poulton N and Stepanauskas R and Bork P and Bowler C and Hingamp P and Sullivan MB and Iudicone D and Massana R and Aury JM and Henrissat B and Karsenti E and Jaillon O and Sieracki M and De Vargas C and Wincker P},

url = {https://www.nature.com/articles/s41467-017-02235-3},

doi = {10.1038/s41467-017-02235-3},

year  = {2018},

date = {2018-01-22},

journal = {Nature Communications},

volume = {9},

abstract = {Single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans. While their roles as primary producers and grazers are well appreciated, other aspects of their life histories remain obscure due to challenges in culturing and sequencing their natural diversity. Here, we exploit single-cell genomics and metagenomics data from the circumglobal Tara Oceans expedition to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles. Based on the available data, each sequenced genome or genotype appears to have a specific oceanic distribution, principally correlated with water temperature and depth. The genome content provides hypotheses for specialization in terms of cell motility, food spectra, and trophic stages, including the potential impact on their lifestyles of horizontal gene transfer from prokaryotes. Our results support the idea that prominent heterotrophic marine protists perform diverse functions in ocean ecology.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{K2017b,

title = {Organic matter processing by microbial communities throughout the Atlantic water column as revealed by metaproteomics},

author = {Bergauer K and Fernandez-Guerra A and Garcia JAL and Sprenger RR and Stepanauskas R and Pachiadaki MG and Jensen ON and Herndl GJ},

url = {http://www.pnas.org/content/early/2017/12/13/1708779115},

doi = {https://doi.org/10.1073/pnas.1708779115},

year  = {2018},

date = {2018-01-16},

journal = {Proceedings of the National Academy of Sciences of the United States of America},

volume = {115},

number = {3},

pages = {E400-E408},

abstract = {The phylogenetic composition of the heterotrophic microbial community is depth stratified in the oceanic water column down to abyssopelagic layers. In the layers below the euphotic zone, it has been suggested that heterotrophic microbes rely largely on solubilized particulate organic matter as a carbon and energy source rather than on dissolved organic matter. To decipher whether changes in the phylogenetic composition with depth are reflected in changes in the bacterial and archaeal transporter proteins, we generated an extensive metaproteomic and metagenomic dataset of microbial communities collected from 100- to 5,000-m depth in the Atlantic Ocean. By identifying which compounds of the organic matter pool are absorbed, transported, and incorporated into microbial cells, intriguing insights into organic matter transformation in the deep ocean emerged. On average, solute transporters accounted for 23% of identified protein sequences in the lower euphotic and ∼39% in the bathypelagic layer, indicating the central role of heterotrophy in the dark ocean. In the bathypelagic layer, substrate affinities of expressed transporters suggest that, in addition to amino acids, peptides and carbohydrates, carboxylic acids and compatible solutes may be essential substrates for the microbial community. Key players with highest expression of solute transporters were Alphaproteobacteria, Gammaproteobacteria, and Deltaproteobacteria, accounting for 40%, 11%, and 10%, respectively, of relative protein abundances. The in situ expression of solute transporters indicates that the heterotrophic prokaryotic community is geared toward the utilization of similar organic compounds throughout the water column, with yet higher abundances of transporters targeting aromatic compounds in the bathypelagic realm.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{HL2017,

title = {Homoacetogenesis in Deep-Sea Chloroflexi, as Inferred by Single-Cell Genomics, Provides a Link to Reductive Dehalogenation in Terrestrial Dehalococcoidetes},

author = {Sewell HL and Kaster AK and Spormann AM},

url = {http://mbio.asm.org/content/8/6/e02022-17.abstract},

doi = {10.1128/mBio.02022-17},

year  = {2017},

date = {2017-12-19},

journal = {mBio},

volume = {8},

number = {6},

abstract = {The deep marine subsurface is one of the largest unexplored biospheres on Earth and is widely inhabited by members of the phylum Chloroflexi. In this report, we investigated genomes of single cells obtained from deep-sea sediments of the Peruvian Margin, which are enriched in such Chloroflexi. 16S rRNA gene sequence analysis placed two of these single-cell-derived genomes (DscP3 and Dsc4) in a clade of subphylum I Chloroflexi which were previously recovered from deep-sea sediment in the Okinawa Trough and a third (DscP2-2) as a member of the previously reported DscP2 population from Peruvian Margin site 1230. The presence of genes encoding enzymes of a complete Wood-Ljungdahl pathway, glycolysis/gluconeogenesis, a Rhodobacter nitrogen fixation (Rnf) complex, glyosyltransferases, and formate dehydrogenases in the single-cell genomes of DscP3 and Dsc4 and the presence of an NADH-dependent reduced ferredoxin:NADP oxidoreductase (Nfn) and Rnf in the genome of DscP2-2 imply a homoacetogenic lifestyle of these abundant marine Chloroflexi. We also report here the first complete pathway for anaerobic benzoate oxidation to acetyl coenzyme A (CoA) in the phylum Chloroflexi (DscP3 and Dsc4), including a class I benzoyl-CoA reductase. Of remarkable evolutionary significance, we discovered a gene encoding a formate dehydrogenase (FdnI) with reciprocal closest identity to the formate dehydrogenase-like protein (complex iron-sulfur molybdoenzyme [CISM], DET0187) of terrestrial Dehalococcoides/Dehalogenimonas spp. This formate dehydrogenase-like protein has been shown to lack formate dehydrogenase activity in Dehalococcoides/Dehalogenimonas spp. and is instead hypothesized to couple HupL hydrogenase to a reductive dehalogenase in the catabolic reductive dehalogenation pathway. This finding of a close functional homologue provides an important missing link for understanding the origin and the metabolic core of terrestrial Dehalococcoides/Dehalogenimonas spp. and of reductive dehalogenation, as well as the biology of abundant deep-sea Chloroflexi.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ED2017,

title = {Rokubacteria: Genomic giants among the uncultured bacterial phyla},

author = {Becraft ED, Woyke T, Jarett J, Ivanova N, Godoy-Vitorino F, Poulton N, Brown JM, Brown J, Lau M, Onstott T, Eisen JA, Moser D, Stepanauskas R},

url = {https://www.frontiersin.org/articles/10.3389/fmicb.2017.02264/full},

doi = {10.3389/fmicb.2017.02264},

year  = {2017},

date = {2017-11-28},

journal = {Frontiers in Microbiology },

volume = {8},

abstract = {Recent advances in single-cell genomic and metagenomic techniques have facilitated the discovery of numerous previously unknown, deep branches of the tree of life that lack cultured representatives. Many of these candidate phyla are composed of microorganisms with minimalistic, streamlined genomes lacking some core metabolic pathways, which may contribute to their resistance to growth in pure culture. Here we analyzed single-cell genomes and metagenome bins to show that the “Candidate phylum Rokubacteria,” formerly known as SPAM, represents an interesting exception, by having large genomes (6\textendash8 Mbps), high GC content (66\textendash71%), and the potential for a versatile, mixotrophic metabolism. We also observed an unusually high genomic heterogeneity among individual Rokubacteria cells in the studied samples. These features may have contributed to the limited recovery of sequences of this candidate phylum in prior cultivation and metagenomic studies. Our analyses suggest that Rokubacteria are distributed globally in diverse terrestrial ecosystems, including soils, the rhizosphere, volcanic mud, oil wells, aquifers, and the deep subsurface, with no reports from marine environments to date.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{MG2017,

title = {Major role of nitrite-oxidizing bacteria in dark ocean carbon fixation},

author = {Pachiadaki MG and Sintes E and Bergauer K and Brown JM and Record NR and Swan BK and Mathyer ME and Hallam SJ and Lopez-Garcia P and Takaki Y and Nunoura T and Woyke T and Herndl GJ and Stepanauskas R},

url = {http://science.sciencemag.org/content/358/6366/1046.long},

doi = {10.1126/science.aan8260},

year  = {2017},

date = {2017-11-24},

journal = {Science},

volume = {358},

pages = {1046-1051},

abstract = {Carbon fixation by chemoautotrophic microorganisms in the dark ocean has a major impact on global carbon cycling and ecological relationships in the ocean’s interior, but the relevant taxa and energy sources remain enigmatic. We show evidence that nitrite-oxidizing bacteria affiliated with the Nitrospinae phylum are important in dark ocean chemoautotrophy. Single-cell genomics and community metagenomics revealed that Nitrospinae are the most abundant and globally distributed nitrite-oxidizing bacteria in the ocean. Metaproteomics and metatranscriptomics analyses suggest that nitrite oxidation is the main pathway of energy production in Nitrospinae. Microautoradiography, linked with catalyzed reporter deposition fluorescence in situ hybridization, indicated that Nitrospinae fix 15 to 45% of inorganic carbon in the mesopelagic western North Atlantic. Nitrite oxidation may have a greater impact on the carbon cycle than previously assumed.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{D2017,

title = {Characterization of two cryptic plasmids isolated in Haiti from clinical Vibrio cholerae non-O1/non-O139},

author = {Ceccarelli D and Garriss G and Choi SY and Hasan NA and Stepanauskas R and Pop M and Huq A and Colwell RR},

url = {https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5703827/},

doi = {10.3389/fmicb.2017.02283},

year  = {2017},

date = {2017-11-23},

journal = {Frontiers in Microbiology},

volume = {8},

abstract = {We report the complete sequence of two novel plasmids, pSDH-1 and pSDH-2, isolated from clinical Vibrio cholerae non-O1/non-O139 during the early phase of the 2010 Haitian cholera epidemic. Plasmids were revealed by employing single-cell genomics and their genome content suggests self-mobilization and, for pSDH-2, a toxin-antitoxin (TA) system for plasmid stabilization was identified. The putative origin of replication of pSDH-2 was mapped suggesting it replicates following the ColE1 model of plasmid replication. pSDH-1 and pSDH-2 were widespread among environmental V. cholerae non-O1/non-O139 with variable prevalence in four Haitian Departments. pSDH-2 was the most common element, either alone or with pSDH-1. The two plasmids detection adds to the composite scenario of mobile genetic elements (MGEs) observed in V. cholerae in Haiti. The role these small cryptic plasmids circulating in Vibrio spp. play in bacterial fitness or pathogenicity merits further investigation.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{AK2017,

title = {Diverse Marinimicrobia bacteria may mediate coupled biogeochemical cycles along eco-thermodynamic gradients},

author = {Hawley AK and Nobu MK and Wright JJ and Durno WE and Morgan-Lang C and Sage B and Schwientek P and Swan BK and Rinke C and Torres-Beltr\'{a}n M and Mewis K and Liu WT and Stepanauskas R and Woyke T and Hallam SJ},

url = {https://www.nature.com/articles/s41467-017-01376-9},

doi = {10.1038/s41467-017-01376-9},

year  = {2017},

date = {2017-11-15},

journal = {Nature Communications},

volume = {8},

abstract = {Microbial communities drive biogeochemical cycles through networks of metabolite exchange that are structured along energetic gradients. As energy yields become limiting, these networks favor co-metabolic interactions to maximize energy disequilibria. Here we apply single-cell genomics, metagenomics, and metatranscriptomics to study bacterial populations of the abundant “microbial dark matter” phylum Marinimicrobia along defined energy gradients. We show that evolutionary diversification of major Marinimicrobia clades appears to be closely related to energy yields, with increased co-metabolic interactions in more deeply branching clades. Several of these clades appear to participate in the biogeochemical cycling of sulfur and nitrogen, filling previously unassigned niches in the ocean. Notably, two Marinimicrobia clades, occupying different energetic niches, express nitrous oxide reductase, potentially acting as a global sink for the greenhouse gas nitrous oxide.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{ED2017b,

title = {Genomic comparison of two family-level groups of the uncultivated NAG1 archaeal lineage from chemically and geographically disparate hot springs},

author = {Becraft ED and Dodsworth JA and Murugapiran SK and Thomas SC and Ohlsson JI and Stepanauskas R and Hedlund BP and Swingley WD},

url = {https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5671600/},

doi = {10.3389/fmicb.2017.02082},

year  = {2017},

date = {2017-10-31},

journal = {Frontiers in Microbiology },

volume = {8},

abstract = {Recent progress based on single-cell genomics and metagenomic investigations of archaea in a variety of extreme environments has led to significant advances in our understanding of the diversity, evolution, and metabolic potential of archaea, yet the vast majority of archaeal diversity remains undersampled. In this work, we coordinated single-cell genomics with metagenomics in order to construct a near-complete genome from a deeply branching uncultivated archaeal lineage sampled from Great Boiling Spring (GBS) in the U.S. Great Basin, Nevada. This taxon is distantly related (distinct families) to an archaeal genome, designated “Novel Archaeal Group 1” (NAG1), which was extracted from a metagenome recovered from an acidic iron spring in Yellowstone National Park (YNP). We compared the metabolic predictions of the NAG1 lineage to better understand how these archaea could inhabit such chemically distinct environments. Similar to the NAG1 population previously studied in YNP, the NAG1 population from GBS is predicted to utilize proteins as a primary carbon source, ferment simple carbon sources, and use oxygen as a terminal electron acceptor under oxic conditions. However, GBS NAG1 populations contained distinct genes involved in central carbon metabolism and electron transfer, including nitrite reductase, which could confer the ability to reduce nitrite under anaerobic conditions. Despite inhabiting chemically distinct environments with large variations in pH, GBS NAG1 populations shared many core genomic and metabolic features with the archaeon identified from YNP, yet were able to carve out a distinct niche at GBS.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{J.J2017,

title = {Metabolic Network Analysis and Metatranscriptomics Reveal Auxotrophies and Nutrient Sources of the Cosmopolitan Freshwater Microbial Lineage acI},

author = {Hamilton J.J and Garcia S.L and Brown B.S and Oyserman B.O and Moya-Flores F and Bertilsson S and Malmstrom R.R and Forest K and McMahon K },

url = {https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5574706/pdf/mSystems.00091-17.pdf},

doi = {10.1128/mSystems.00091-17},

year  = {2017},

date = {2017-08-29},

journal = {American Society for Microbiology},

volume = {2},

number = {4},

pages = {e00091-17},

abstract = {An explosion in the number of available genome sequences obtained through metagenomics and single-cell genomics has enabled a new view of the diversity of microbial life, yet we know surprisingly little about how microbes interact with each other or their environment. In fact, the majority of microbial species remain uncultivated, while our perception of their ecological niches is based on reconstruction of their metabolic potential. In this work, we demonstrate how the "seed set framework," which computes the set of compounds that an organism must acquire from its environment (E. Borenstein, M. Kupiec, M. W. Feldman, and E. Ruppin, Proc Natl Acad Sci U S A 105:14482-14487, 2008, https://doi.org/10.1073/pnas.0806162105), enables computational analysis of metabolic reconstructions while providing new insights into a microbe's metabolic capabilities, such as nutrient use and auxotrophies. We apply this framework to members of the ubiquitous freshwater actinobacterial lineage acI, confirming and extending previous experimental and genomic observations implying that acI bacteria are heterotrophs reliant on peptides and saccharides. We also present the first metatranscriptomic study of the acI lineage, revealing high expression of transport proteins and the light-harvesting protein actinorhodopsin. Putative transport proteins complement predictions of nutrients and essential metabolites while providing additional support of the hypothesis that members of the acI are photoheterotrophs. IMPORTANCE The metabolic activity of uncultivated microorganisms contributes to numerous ecosystem processes, ranging from nutrient cycling in the environment to influencing human health and disease. Advances in sequencing technology have enabled the assembly of genomes for these microorganisms, but our ability to generate reference genomes far outstrips our ability to analyze them. Common approaches to analyzing microbial metabolism require reconstructing the entirety of an organism's metabolic pathways or performing targeted searches for genes involved in a specific process. This paper presents a third approach, in which draft metabolic reconstructions are used to identify compounds through which an organism may interact with its environment. These compounds can then guide more-intensive metabolic reconstruction efforts and can also provide new hypotheses about the specific contributions that microbes make to ecosystem-scale metabolic processes.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{JF2017b,

title = {Physiological and ecological implications of an iron- or hydrogen-oxidizing member of the Zetaproteobacteria, Ghiorsea bivora, gen. nov., sp. nov},

author = {Mori JF and Scott JJ and Hager KW and Moyer CL and Kusel K and Emerson D},

url = {https://www.nature.com/articles/ismej2017132},

doi = {10.1038/ismej.2017.132},

year  = {2017},

date = {2017-08-18},

journal = {The ISME Journal},

volume = {11},

pages = {2624-2636},

abstract = {Chemosynthetic Fe-oxidizing communities are common at diffuse-flow hydrothermal vents throughout the world's oceans. The foundational members of these communities are the Zetaproteobacteria, a class of Proteobacteria that is primarily associated with ecosystems fueled by ferrous iron, Fe(II). We report here the discovery of two new isolates of Zetaproteobacteria isolated from the Mid-Atlantic Ridge (TAG-1), and the Mariana back-arc (SV-108), that are unique in that they can utilize either Fe(II) or molecular hydrogen (H2) as sole electron donor and oxygen as terminal electron acceptor for growth. Both strains precipitated Fe-oxyhydroxides as amorphous particulates. The cell doubling time on H2 vs Fe(II) for TAG-1 was 14.1 vs 21.8 h, and for SV-108 it was 16.3 vs 20 h, and it appeared both strains could use either H2 or Fe(II) simultaneously. The strains were close relatives, based on genomic analysis, and both possessed genes for the uptake NiFe-hydrogenase required for growth on H2. These two strains belong to Zetaproteobacteria operational taxonomic unit 9 (ZetaOTU9). A meta-analysis of public databases found ZetaOTU9 was only associated with Fe(II)-rich habitats, and not in other environments where known H2-oxidizers exist. These results expand the metabolic repertoire of the Zetaproteobacteria, yet confirm that Fe(II) metabolism is the primary driver of their physiology and ecology.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{RM2017b,

title = {Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea},

author = {Bowers RM and Kyrpides NC and Stepanauskas R and Harmon-Smith M and Doud D and Reddy TBK and Schulz F and Jarett J and Rivers AR and Eloe-Fadrosh EA and Tringe SG and Ivanova NN and Copeland A and Clum A and Becraft ED and Malmstrom RR and Birren B and Podar M and Bork P and Weinstock GM and Garrity GM and Dodsworth JA and Yooseph S and Sutton G and Gl\"{o}ckner FO and Gilbert JA and Nelson WC and Hallam SJ and Jungbluth SP and Ettema TJG and Tighe S and Konstantinidis KT and Liu WT and Baker BJ and Rattei T and Eisen JA and Hedlund B and McMahon KD and Fierer N and Knight R and Finn R and Cochrane G and Karsch-Mizrachi I and Tyson GW and Rinke C and Lapidus A and Meyer F and Yilmaz P and Parks DH and Eren AM and Schriml L and Banfield JF and Hugenholtz P and Woyke T},

url = {https://www.nature.com/articles/nbt.3893},

doi = {10.1038/nbt.3893},

year  = {2017},

date = {2017-08-08},

journal = {Nature Biotechnology},

volume = {35},

pages = {725-731},

abstract = {We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{A2017b,

title = {Viral to metazoan marine plankton nucleotide sequences from the Tara Oceans expedition},

author = {Alberti A and Poulain J and Engelen S and Labadie K and Romac S and Ferrera I and Albini G and Aury JM and Belser C and Bertrand A and Cruaud C and Da Silva C and Dossat C and Gavory F and Gas S and Guy J and Haquelle M and Jacoby E and Jaillon O and Lemainque A and Pelletier E and Samson G and Wessner M and Acinas SG and Royo-Llonch M and Cornejo-Castillo FM and Logares R and Fern\'{a}ndez-G\'{o}mez B and Bowler C and Cochrane G and Amid C and Hoopen PT and De Vargas C and Grimsley N and Desgranges E and Kandels-Lewis S and Ogata H and Poulton N and Sieracki ME and Stepanauskas R and Sullivan MB and Brum JR and Duhaime MB and Poulos BT and Hurwitz BL and Pesant S and Karsenti E and Wincker P and Bazire P and Beluche O and Bertrand L and Besnard-Gonnet M and Bordelais I and Boutard M and Dubois M and Dumont C and Ettedgui E and Fernandez P and Garcia E and Aiach NG and Guerin T and Hamon C and Brun E and Lebled S and Lenoble P and Louesse C and Mahieu E and Mairey B and Martins N and Megret C and Milani C and Muanga J and Orvain C and Payen E and Perroud P and Petit E and Robert D and Ronsin M and Vacherie B and Bork P and Boss E and Follows M and Gorsky G and Hingamp P and Iudicone D and Karp-Boss L and Not F and Raes J and Sardet C and Speich S and Stemmann L and Sunagawa S},

url = {https://www.nature.com/articles/sdata201793},

doi = {10.1038/sdata.2017.93},

year  = {2017},

date = {2017-08-01},

journal = {Scientific Data},

volume = {4},

abstract = {A unique collection of oceanic samples was gathered by the Tara Oceans expeditions (2009\textendash2013), targeting plankton organisms ranging from viruses to metazoans, and providing rich environmental context measurements. Thanks to recent advances in the field of genomics, extensive sequencing has been performed for a deep genomic analysis of this huge collection of samples. A strategy based on different approaches, such as metabarcoding, metagenomics, single-cell genomics and metatranscriptomics, has been chosen for analysis of size-fractionated plankton communities. Here, we provide detailed procedures applied for genomic data generation, from nucleic acids extraction to sequence production, and we describe registries of genomics datasets available at the European Nucleotide Archive (ENA, www.ebi.ac.uk/ena). The association of these metadata to the experimental procedures applied for their generation will help the scientific community to access these data and facilitate their analysis. This paper complements other efforts to provide a full description of experiments and open science resources generated from the Tara Oceans project, further extending their value for the study of the world’s planktonic ecosystems.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{M2017b,

title = {Exploring microdiversity in novel Kordia sp. (Bacteroidetes) with proteorhodopsin from the tropical Indian Ocean via Single Amplified Genomes},

author = {Royo-Llonch M and Ferrera I and Cornejo-Castillo FM and S\'{a}nchez P and Salazar G and Stepanauskas R and Gonz\'{a}lez JM and Sieracki ME and Speich S and Stemmann L and Pedr\'{o}s-Ali\'{o} C and Acinas SG},

url = {https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5525439/},

doi = {10.3389/fmicb.2017.01317},

year  = {2017},

date = {2017-07-25},

journal = {Frontiers in Microbiology},

volume = {8},

abstract = {Marine Bacteroidetes constitute a very abundant bacterioplankton group in the oceans that plays a key role in recycling particulate organic matter and includes several photoheterotrophic members containing proteorhodopsin. Relatively few marine Bacteroidetes species have been described and, moreover, they correspond to cultured isolates, which in most cases do not represent the actual abundant or ecologically relevant microorganisms in the natural environment. In this study, we explored the microdiversity of 98 Single Amplified Genomes (SAGs) retrieved from the surface waters of the underexplored North Indian Ocean, whose most closely related isolate is Kordia algicida OT-1. Using Multi Locus Sequencing Analysis (MLSA) we found no microdiversity in the tested conserved phylogenetic markers (16S rRNA and 23S rRNA genes), the fast-evolving Internal Transcribed Spacer and the functional markers proteorhodopsin and the beta-subunit of RNA polymerase. Furthermore, we carried out a Fragment Recruitment Analysis (FRA) with marine metagenomes to learn about the distribution and dynamics of this microorganism in different locations, depths and size fractions. This analysis indicated that this taxon belongs to the rare biosphere, showing its highest abundance after upwelling-induced phytoplankton blooms and sinking to the deep ocean with large organic matter particles. This uncultured Kordia lineage likely represents a novel Kordia species (Kordia sp. CFSAG39SUR) that contains the proteorhodopsin gene and has a widespread spatial and vertical distribution. The combination of SAGs and MLSA makes a valuable approach to infer putative ecological roles of uncultured abundant microorganisms.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{R2017,

title = {Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles},

author = {Stepanauskas R and Fergusson EA and Brown J and Poulton NJ and Tupper B and Labont\'{e} JM and Becraft ED and Brown JM and Pachiadaki MJ and Povilaitis T and Thompson BP and Mascena CJ and Bellows WK and Lubys A},

url = {https://www.nature.com/articles/s41467-017-00128-z},

doi = {10.1038/s41467-017-00128-z},

year  = {2017},

date = {2017-07-20},

journal = {Nature Communications},

volume = {8},

number = {84},

abstract = {Microbial single-cell genomics can be used to provide insights into the metabolic potential, interactions, and evolution of uncultured microorganisms. Here we present WGA-X, a method based on multiple displacement amplification of DNA that utilizes a thermostable mutant of the phi29 polymerase. WGA-X enhances genome recovery from individual microbial cells and viral particles while maintaining ease of use and scalability. The greatest improvements are observed when amplifying high G+C content templates, such as those belonging to the predominant bacteria in agricultural soils. By integrating WGA-X with calibrated index-cell sorting and high-throughput genomic sequencing, we are able to analyze genomic sequences and cell sizes of hundreds of individual, uncultured bacteria, archaea, protists, and viral particles, obtained directly from marine and soil samples, in a single experiment. This approach may find diverse applications in microbiology and in biomedical and forensic studies of humans and other multicellular organisms.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{A2017,

title = {Unexpected genomic features in widespread intracellular bacteria: evidence for motility of marine chlamydiae},

author = {Collingro A and Köstlbacher S and Mussmann M and Stepanauskas R and Hallam S and Horn M},

url = {https://www.nature.com/ismej/journal/vaop/ncurrent/full/ismej201795a.html},

doi = {10.1038/ismej.2017.95},

year  = {2017},

date = {2017-06-23},

journal = {ISME},

abstract = {Chlamydiae are obligate intracellular bacteria comprising important human pathogens and symbionts of protists. Molecular evidence indicates a tremendous diversity of chlamydiae particularly in marine environments, yet our current knowledge is based mainly on terrestrial representatives. Here we provide first insights into the biology of marine chlamydiae representing three divergent clades. Our analysis of single-cell amplified genomes revealed hallmarks of the chlamydial lifestyle, supporting the ancient origin of their characteristic developmental cycle and major virulence mechanisms. Surprisingly, these chlamydial genomes encode a complete flagellar apparatus, a previously unreported feature. We show that flagella are an ancient trait that was subject to differential gene loss among extant chlamydiae. Together with a chemotaxis system, these marine chlamydiae are likely motile, with flagella potentially playing a role during host cell infection. This study broadens our view on chlamydial biology and indicates a largely underestimated potential to adapt to different hosts and environments.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{Luo2017,

title = {Excess of non-conservative amino acid changes in marine bacterioplankton lineages with reduced genomes},

author = {Haiwei Luo and Yongjie Huang and Ramunas Stepanauskas and Jijun Tang},

url = {https://www.nature.com/articles/nmicrobiol201791},

doi = {10.1007/s11096-011-9488-z},

year  = {2017},

date = {2017-06-12},

journal = {Nature Microbiology},

volume = {2},

number = {17091},

abstract = {Surface ocean waters are dominated by planktonic bacterial lineages with highly reduced genomes. The best examples are the cyanobacterial genus Prochlorococcus, the alphaproteobacterial clade SAR11 and the gammaproteobacterial clade SAR86, which together represent over 50% of the cells in surface oceans. Several studies have identified signatures of selection on these lineages in today's ocean and have postulated selection as the primary force throughout their evolutionary history. However, massive loss of genomic DNA in these lineages often occurred in the distant past, and the selective pressures underlying these ancient events have not been assessed. Here, we probe ancient selective pressures by computing %GC-corrected rates of conservative and radical nonsynonymous nucleotide substitutions. Surprisingly, we found an excess of radical changes in several of these lineages in comparison to their relatives with larger genomes. Furthermore, analyses of allelic genome sequences of several populations within these lineages consistently supported that radical replacements are more likely to be deleterious than conservative changes. Our results suggest coincidence of massive genomic DNA losses and increased power of genetic drift, but we also suggest that additional evidence independent of the nucleotide substitution analyses is needed to support a primary role of genetic drift driving ancient genome reduction of marine bacterioplankton lineages.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{N2017,

title = {Fundamental differences in diversity and genomic population structure between Atlantic and Pacific Prochlorococcus},

author = {Kashtan N and Roggensack SE and Berta-Thompson JW and Grinberg M and Stepanauskas R and Chisholm SW},

url = {http://www.nature.com/ismej/journal/vaop/ncurrent/full/ismej201764a.html},

doi = {10.1038/ismej.2017.64},

year  = {2017},

date = {2017-05-19},

journal = {ISME},

abstract = {The Atlantic and Pacific Oceans represent different biogeochemical regimes in which the abundant marine cyanobacterium Prochlorococcus thrives. We have shown that Prochlorococcus populations in the Atlantic are composed of hundreds of genomically, and likely ecologically, distinct coexisting subpopulations with distinct genomic backbones. Here we ask if differences in the ecology and selection pressures between the Atlantic and Pacific are reflected in the diversity and genomic composition of their indigenous Prochlorococcus populations. We applied large-scale single-cell genomics and compared the cell-by-cell genomic composition of wild populations of co-occurring cells from samples from Station ALOHA off Hawaii, and from Bermuda Atlantic Time Series Station off Bermuda. We reveal fundamental differences in diversity and genomic structure of populations between the sites. The Pacific populations are more diverse than those in the Atlantic, composed of significantly more coexisting subpopulations and lacking dominant subpopulations. Prochlorococcus from the two sites seem to be composed of mostly non-overlapping distinct sets of subpopulations with different genomic backbones\textemdashlikely reflecting different sets of ocean-specific micro-niches. Furthermore, phylogenetically closely related strains carry ocean-associated nutrient acquisition genes likely reflecting differences in major selection pressures between the oceans. This differential selection, along with geographic separation, clearly has a significant role in shaping these populations.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}

@article{WH2017,

title = {Genomic exploration of individual giant ocean viruses},

author = {Wilson WH and Gilg IC and Moniruzzaman M and Field EK and Koren S and LeCleir GR and Martínez Martínez J and Poulton NJ and Swan BK and Stepanauskas R and Wilhelm SW},

url = {https://www.nature.com/ismej/journal/vaop/ncurrent/full/ismej201761a.html},

doi = {10.1038/ismej.2017.61},

year  = {2017},

date = {2017-05-12},

journal = {ISME},

abstract = {Viruses are major pathogens in all biological systems. Virus propagation and downstream analysis remains a challenge, particularly in the ocean where the majority of their microbial hosts remain recalcitrant to current culturing techniques. We used a cultivation-independent approach to isolate and sequence individual viruses. The protocol uses high-speed fluorescence-activated virus sorting flow cytometry, multiple displacement amplification (MDA), and downstream genomic sequencing. We focused on ‘giant viruses’ that are readily distinguishable by flow cytometry. From a single-milliliter sample of seawater collected from off the dock at Boothbay Harbor, ME, USA, we sorted almost 700 single virus particles, and subsequently focused on a detailed genome analysis of 12. A wide diversity of viruses was identified that included Iridoviridae, extended Mimiviridae and even a taxonomically novel (unresolved) giant virus. We discovered a viral metacaspase homolog in one of our sorted virus particles and discussed its implications in rewiring host metabolism to enhance infection. In addition, we demonstrated that viral metacaspases are widespread in the ocean. We also discovered a virus that contains both a reverse transcriptase and a transposase; although highly speculative, we suggest such a genetic complement would potentially allow this virus to exploit a latency propagation mechanism. Application of single virus genomics provides a powerful opportunity to circumvent cultivation of viruses, moving directly to genomic investigation of naturally occurring viruses, with the assurance that the sequence data is virus-specific, non-chimeric and contains no cellular contamination.},

keywords = {},

pubstate = {published},

tppubtype = {article}

}