Ecogenomics of virophages and their giant virus hosts assessed through time series metagenomics.

Virophages are small viruses that co-infect eukaryotic cells alongside giant viruses (Mimiviridae) and hijack their machinery to replicate. While two types of virophages have been isolated, their genomic diversity and ecology remain largely unknown. Here we use time series metagenomics to identify and study the dynamics of 25 uncultivated virophage populations, 17 of which represented by complete or near-complete genomes, in two North American freshwater lakes. Taxonomic analysis suggests that these freshwater virophages represent at least three new candidate genera. Ecologically, virophage populations are repeatedly detected over years and evolutionary stable, yet their distinct abundance profiles and gene content suggest that virophage genera occupy different ecological niches. Co-occurrence analyses reveal 11 virophages strongly associated with uncultivated Mimiviridae, and three associated with eukaryotes among the Dinophyceae, Rhizaria, Alveolata, and Cryptophyceae groups. Together, these findings significantly augment virophage databases, help refine virophage taxonomy, and establish baseline ecological hypotheses and tools to study virophages in nature.Virophages are recently-identified small viruses that infect larger viruses, yet their diversity and ecological roles are poorly understood. Here, Roux and colleagues present time series metagenomics data revealing new virophage genera and their putative ecological interactions in two freshwater lakes.

Ecophysiology of Freshwater Verrucomicrobia Inferred from Metagenome-Assembled Genomes.

Microbes are critical in carbon and nutrient cycling in freshwater ecosystems. Members of the Verrucomicrobia are ubiquitous in such systems, and yet their roles and ecophysiology are not well understood. In this study, we recovered 19 Verrucomicrobia draft genomes by sequencing 184 time-series metagenomes from a eutrophic lake and a humic bog that differ in carbon source and nutrient availabilities. These genomes span four of the seven previously defined Verrucomicrobia subdivisions and greatly expand knowledge of the genomic diversity of freshwater Verrucomicrobia. Genome analysis revealed their potential role as (poly)saccharide degraders in freshwater, uncovered interesting genomic features for this lifestyle, and suggested their adaptation to nutrient availabilities in their environments. Verrucomicrobia populations differ significantly between the two lakes in glycoside hydrolase gene abundance and functional profiles, reflecting the autochthonous and terrestrially derived allochthonous carbon sources of the two ecosystems, respectively. Interestingly, a number of genomes recovered from the bog contained gene clusters that potentially encode a novel porin-multiheme cytochrome c complex and might be involved in extracellular electron transfer in the anoxic humus-rich environment. Notably, most epilimnion genomes have large numbers of so-called "Planctomycete-specific" cytochrome c-encoding genes, which exhibited distribution patterns nearly opposite to those seen with glycoside hydrolase genes, probably associated with the different levels of environmental oxygen availability and carbohydrate complexity between lakes/layers. Overall, the recovered genomes represent a major step toward understanding the role, ecophysiology, and distribution of Verrucomicrobia in freshwater. IMPORTANCE Freshwater Verrucomicrobia spp. are cosmopolitan in lakes and rivers, and yet their roles and ecophysiology are not well understood, as cultured freshwater Verrucomicrobia spp. are restricted to one subdivision of this phylum. Here, we greatly expanded the known genomic diversity of this freshwater lineage by recovering 19 Verrucomicrobia draft genomes from 184 metagenomes collected from a eutrophic lake and a humic bog across multiple years. Most of these genomes represent the first freshwater representatives of several Verrucomicrobia subdivisions. Genomic analysis revealed Verrucomicrobia to be potential (poly)saccharide degraders and suggested their adaptation to carbon sources of different origins in the two contrasting ecosystems. We identified putative extracellular electron transfer genes and so-called "Planctomycete-specific" cytochrome c-encoding genes and identified their distinct distribution patterns between the lakes/layers. Overall, our analysis greatly advances the understanding of the function, ecophysiology, and distribution of freshwater Verrucomicrobia, while highlighting their potential role in freshwater carbon cycling.

Nanogrid single-nucleus RNA sequencing reveals phenotypic diversity in breast cancer.

Single cell RNA sequencing has emerged as a powerful tool for resolving transcriptional diversity in tumors, but is limited by throughput, cost and the ability to process archival frozen tissue samples. Here we develop a high-throughput 3' single-nucleus RNA sequencing approach that combines nanogrid technology, automated imaging, and cell selection to sequence up to ~1800 single nuclei in parallel. We compare the transcriptomes of 485 single nuclei to 424 single cells in a breast cancer cell line, which shows a high concordance (93.34%) in gene levels and abundance. We also analyze 416 nuclei from a frozen breast tumor sample and 380 nuclei from normal breast tissue. These data reveal heterogeneity in cancer cell phenotypes, including angiogenesis, proliferation, and stemness, and a minor subpopulation (19%) with many overexpressed cancer genes. Our studies demonstrate the utility of nanogrid single-nucleus RNA sequencing for studying the transcriptional programs of tumor nuclei in frozen archival tissue samples.Single cell RNA sequencing is a powerful tool for understanding cellular diversity but is limited by cost, throughput and sample preparation. Here the authors use nanogrid technology with integrated imaging to sequence thousands of cancer nuclei in parallel from fresh or frozen tissue.

Experimental Identification of Small Non-Coding RNAs in the Model Marine Bacterium Ruegeria pomeroyi DSS-3.

In oligotrophic ocean waters where bacteria are often subjected to chronic nutrient limitation, community transcriptome sequencing has pointed to the presence of highly abundant small RNAs (sRNAs). The role of sRNAs in regulating response to nutrient stress was investigated in a model heterotrophic marine bacterium Ruegeria pomeroyi grown in continuous culture under carbon (C) and nitrogen (N) limitation. RNAseq analysis identified 99 putative sRNAs. Sixty-nine were cis-encoded and located antisense to a presumed target gene. Thirty were trans-encoded and initial target prediction was performed computationally. The most prevalent functional roles of genes anti-sense to the cis-sRNAs were transport, cell-cell interactions, signal transduction, and transcriptional regulation. Most sRNAs were transcribed equally under both C and N limitation, and may be involved in a general stress response. However, 14 were regulated differentially between the C and N treatments and may respond to specific nutrient limitations. A network analysis of the predicted target genes of the R. pomeroyi cis-sRNAs indicated that they average fewer connections than typical protein-encoding genes, and appear to be more important in peripheral or niche-defining functions encoded in the pan genome.

Genome-wide selective sweeps and gene-specific sweeps in natural bacterial populations.

Multiple models describe the formation and evolution of distinct microbial phylogenetic groups. These evolutionary models make different predictions regarding how adaptive alleles spread through populations and how genetic diversity is maintained. Processes predicted by competing evolutionary models, for example, genome-wide selective sweeps vs gene-specific sweeps, could be captured in natural populations using time-series metagenomics if the approach were applied over a sufficiently long time frame. Direct observations of either process would help resolve how distinct microbial groups evolve. Here, from a 9-year metagenomic study of a freshwater lake (2005-2013), we explore changes in single-nucleotide polymorphism (SNP) frequencies and patterns of gene gain and loss in 30 bacterial populations. SNP analyses revealed substantial genetic heterogeneity within these populations, although the degree of heterogeneity varied by >1000-fold among populations. SNP allele frequencies also changed dramatically over time within some populations. Interestingly, nearly all SNP variants were slowly purged over several years from one population of green sulfur bacteria, while at the same time multiple genes either swept through or were lost from this population. These patterns were consistent with a genome-wide selective sweep in progress, a process predicted by the 'ecotype model' of speciation but not previously observed in nature. In contrast, other populations contained large, SNP-free genomic regions that appear to have swept independently through the populations prior to the study without purging diversity elsewhere in the genome. Evidence for both genome-wide and gene-specific sweeps suggests that different models of bacterial speciation may apply to different populations coexisting in the same environment.

Bromodeoxyuridine labelling and fluorescence-activated cell sorting of polyamine-transforming bacterioplankton in coastal seawater.

Polyamines (PAs) are a group of nitrogen-rich dissolved organic nitrogen (DON) compounds that are ubiquitously distributed in marine environments. To identify bacteria that are involved in PA transformations, coastal bacterioplankton microcosms were amended with a single PA model compound, i.e. putrescine (PUT) or spermidine (SPD), or with no addition as controls (CTRs). Bromodeoxyuridine (BrdU) was added to all the microcosms to label newly synthesized DNAs. Fluorescence-activated cell sorting (FACS) analysis indicated significant increases in numbers of total cells and cells with both high and low levels of BrdU incorporation in the PUT and SPD microcosms, but not in the CTRs. 16S rDNA pyrotag sequencing of FACS-sorted cells indicated that PUT- and SPD-transforming bacteria were composed similarly of a diverse group of taxa affiliated with Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria (especially Roseobacter of its alpha lineage). Broad taxonomic distribution of PA-transforming bacteria was also indicated by the abundance and distribution of PA transporter gene homologues in a survey of sequenced marine bacterial genomes. Our results suggest that PAs may be common DON substrates for marine bacterioplankton, in line with the hypothesis that bacterially mediated PA transformation accounts for an important proportion of marine DON flux.

Comparative single-cell genomics reveals potential ecological niches for the freshwater acI Actinobacteria lineage.

Members of the acI lineage of Actinobacteria are the most abundant microorganisms in most freshwater lakes; however, our understanding of the keys to their success and their role in carbon and nutrient cycling in freshwater systems has been hampered by the lack of pure cultures and genomes. We obtained draft genome assemblies from 11 single cells representing three acI tribes (acI-A1, acI-A7, acI-B1) from four temperate lakes in the United States and Europe. Comparative analysis of acI SAGs and other available freshwater bacterial genomes showed that acI has more gene content directed toward carbohydrate acquisition as compared to Polynucleobacter and LD12 Alphaproteobacteria, which seem to specialize more on carboxylic acids. The acI genomes contain actinorhodopsin as well as some genes involved in anaplerotic carbon fixation indicating the capacity to supplement their known heterotrophic lifestyle. Genome-level differences between the acI-A and acI-B clades suggest specialization at the clade level for carbon substrate acquisition. Overall, the acI genomes appear to be highly streamlined versions of Actinobacteria that include some genes allowing it to take advantage of sunlight and N-rich organic compounds such as polyamines, di- and oligopeptides, branched-chain amino acids and cyanophycin. This work significantly expands the known metabolic potential of the cosmopolitan freshwater acI lineage and its ecological and genetic traits.

Sizing up metatranscriptomics.

A typical marine bacterial cell in coastal seawater contains only ∼200 molecules of mRNA, each of which lasts only a few minutes before being degraded. Such a surprisingly small and dynamic cellular mRNA reservoir has important implications for understanding the bacterium's responses to environmental signals, as well as for our ability to measure those responses. In this perspective, we review the available data on transcript dynamics in environmental bacteria, and then consider the consequences of a small and transient mRNA inventory for functional metagenomic studies of microbial communities.

Transcriptional changes underlying elemental stoichiometry shifts in a marine heterotrophic bacterium.

Marine bacteria drive the biogeochemical processing of oceanic dissolved organic carbon (DOC), a 750-Tg C reservoir that is a critical component of the global C cycle. Catabolism of DOC is thought to be regulated by the biomass composition of heterotrophic bacteria, as cells maintain a C:N:P ratio of ∼50:10:1 during DOC processing. Yet a complicating factor in stoichiometry-based analyses is that bacteria can change the C:N:P ratio of their biomass in response to resource composition. We investigated the physiological mechanisms of resource-driven shifts in biomass stoichiometry in continuous cultures of the marine heterotrophic bacterium Ruegeria pomeroyi (a member of the Roseobacter clade) under four element limitation regimes (C, N, P, and S). Microarray analysis indicated that the bacterium scavenged for alternate sources of the scarce element when cells were C-, N-, or P-limited; reworked the ratios of biomolecules when C- and P- limited; and exerted tighter control over import/export and cytoplasmic pools when N-limited. Under S limitation, a scenario not existing naturally for surface ocean microbes, stress responses dominated transcriptional changes. Resource-driven changes in C:N ratios of up to 2.5-fold and in C:P ratios of up to sixfold were measured in R. pomeroyi biomass. These changes were best explained if the C and P content of the cells was flexible in the face of shifting resources but N content was not, achieved through the net balance of different transcriptional strategies. The cellular-level metabolic trade-offs that govern biomass stoichiometry in R. pomeroyi may have implications for global carbon cycling if extendable to other heterotrophic bacteria. Strong homeostatic responses to N limitation by marine bacteria would intensify competition with autotrophs. Modification of cellular inventories in C- and P-limited heterotrophs would vary the elemental ratio of particulate organic matter sequestered in the deep ocean.

Beyond the genome: functional studies of phototrophic sulfur oxidation.

The increasing availability of complete genomic sequences for cultured phototrophic bacteria and assembled metagenomes from environments dominated by phototrophs has reinforced the need for a "post-genomic" analytical effort to test models of cellular structure and function proposed from genomic data. Comparative genomics has produced a testable model for pathways of sulfur compound oxidation in the phototrophic bacteria. In the case of sulfide, two enzymes are predicted to oxidize sulfide: sulfide:quinone oxidoreductase and flavocytochrome c sulfide dehydrogenase. However, these models do not predict which enzyme is important under what conditions. In Chlorobaculum tepidum, a model green sulfur bacterium, a combination of genetics and physiological analysis of mutant strains has led to the realization that this organism contains at least two active sulfide:quinone oxidoreductases and that there is significant interaction between sulfide oxidation and light harvesting. In the case of elemental sulfur, an organothiol intermediate of unknown structure has been proposed to activate elemental sulfur for transport into the cytoplasm where it can be oxidized or assimilated, and recent approaches using classical metabolite analysis have begun to shed light on this issue both in C. tepidum and the purple sulfur bacterium Allochromatium vinosum.

Chlorobaculum tepidum regulates chlorosome structure and function in response to temperature and electron donor availability.

Green sulfur bacteria (GSB) rely on the chlorosome, a light-harvesting apparatus comprised almost entirely of self-organizing arrays of bacteriochlorophyll (BChl) molecules, to harvest light energy and pass it to the reaction center. In Chlorobaculum tepidum, over 97% of the total BChl is made up of a mixture of four BChl c homologs in the chlorosome that differ in the number and identity of alkyl side chains attached to the chlorin ring. C. tepidum has been reported to vary the distribution of BChl c homologs with growth light intensity, with the highest degree of BChl c alkylation observed under low-light conditions. Here, we provide evidence that this functional response at the level of the chlorosome can be induced not only by light intensity, but also by temperature and a mutation that prevents phototrophic thiosulfate oxidation. Furthermore, we show that in conjunction with these functional adjustments, the fraction of cellular volume occupied by chlorosomes was altered in response to environmental conditions that perturb the balance between energy absorbed by the light-harvesting apparatus and energy utilized by downstream metabolic reactions.

Functional analysis of three sulfide:quinone oxidoreductase homologs in Chlorobaculum tepidum.

Sulfide:quinone oxidoreductase (SQR) catalyzes sulfide oxidation during sulfide-dependent chemo- and phototrophic growth in bacteria. The green sulfur bacterium Chlorobaculum tepidum (formerly Chlorobium tepidum) can grow on sulfide as the sole electron donor and sulfur source. C. tepidum contains genes encoding three SQR homologs: CT0117, CT0876, and CT1087. This study examined which, if any, of the SQR homologs possess sulfide-dependent ubiquinone reduction activity and are required for growth on sulfide. In contrast to CT0117 and CT0876, transcripts of CT1087 were detected only when cells actively oxidized sulfide. Mutation of CT0117 or CT1087 in C. tepidum decreased SQR activity in membrane fractions, and the CT1087 mutant could not grow with >or=6 mM sulfide. Mutation of both CT0117 and CT1087 in C. tepidum completely abolished SQR activity, and the double mutant failed to grow with >or=4 mM sulfide. A C-terminal His(6)-tagged CT1087 protein was membrane localized, as was SQR activity. Epitope-tagged CT1087 was detected only when sulfide was actively consumed by cells. Recombinantly produced CT1087 and CT0117 proteins had SQR activity, while CT0876 did not. In summary, we conclude that, under the conditions tested, both CT0117 and CT1087 function as SQR proteins in C. tepidum. CT0876 may support the growth of C. tepidum at low sulfide concentrations, but no evidence was found for SQR activity associated with this protein.

A genomic region required for phototrophic thiosulfate oxidation in the green sulfur bacterium Chlorobium tepidum (syn. Chlorobaculum tepidum).

The specific enzymes employed by Chlorobium tepidum for the anaerobic oxidation of thiosulfate, sulfide and elemental sulfur during anoxygenic photosynthesis are not well defined. In particular, it is unclear how C. tepidum completely oxidizes thiosulfate. A C. tepidum genomic region, encoding a putative quinone-interacting membrane-bound oxidoreductase (Qmo) complex (CT0866-0868), hypothetical proteins (CT0869-0875) and a sulfide : quinone oxidoreductase (SQR) homologue (CT0876), was analysed for its role in anaerobic sulfur oxidation. Transcripts of genes encoding the Qmo complex, which is similar to archaeal heterodisulfide reductases, were detected by RT-PCR only while sulfide or elemental sulfur were being oxidized, whereas the SQR homologue and CT0872 were expressed during thiosulfate oxidation and into early stationary phase. A mutant of C. tepidum was obtained in which the region between CT0868 and CT0876 was replaced by a transposon insertion resulting in the truncation or deletion of nine genes. This strain, C5, was completely defective for growth on thiosulfate as the sole electron donor in C. tepidum, but only slightly defective for growth on sulfide or thiosulfate plus sulfide. Strain C5 did not oxidize thiosulfate and also displayed a defect in acetate assimilation under all growth conditions. A gene of unknown function, CT0872, deleted in strain C5 that is conserved in chemolithotrophic sulfur-oxidizing bacteria and archaea is the most likely candidate for the thiosulfate oxidation phenotype observed in this strain. The defect in acetate assimilation may be explained by deletion of CT0874, which encodes a homologue of 3-oxoacyl acyl carrier protein synthase.