Expansion of the global RNA virome reveals diverse clades of bacteriophages.

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Reference-free structural variant detection in microbiomes via long-read co-assembly graphs.

MOTIVATION: The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining. RESULTS: We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux. AVAILABILITY AND IMPLEMENTATION: rhea is open source and available at: https://github.com/treangenlab/rhea.

Retrotransposition facilitated the establishment of a primary plastid in the thecate amoeba Paulinella.

The evolution of eukaryotic life was predicated on the development of organelles such as mitochondria and plastids. During this complex process of organellogenesis, the host cell and the engulfed prokaryote became genetically codependent, with the integration of genes from the endosymbiont into the host nuclear genome and subsequent gene loss from the endosymbiont. This process required that horizontally transferred genes become active and properly regulated despite inherent differences in genetic features between donor (endosymbiont) and recipient (host). Although this genetic reorganization is considered critical for early stages of organellogenesis, we have little knowledge about the mechanisms governing this process. The photosynthetic amoeba Paulinella micropora offers a unique opportunity to study early evolutionary events associated with organellogenesis and primary endosymbiosis. This amoeba harbors a "chromatophore," a nascent photosynthetic organelle derived from a relatively recent cyanobacterial association (∼120 million years ago) that is independent of the evolution of primary plastids in plants (initiated ∼1.5 billion years ago). Analysis of the genome and transcriptome of Paulinella revealed that retrotransposition of endosymbiont-derived nuclear genes was critical for their domestication in the host. These retrocopied genes involved in photoprotection in cyanobacteria became expanded gene families and were "rewired," acquiring light-responsive regulatory elements that function in the host. The establishment of host control of endosymbiont-derived genes likely enabled the cell to withstand photo-oxidative stress generated by oxygenic photosynthesis in the nascent organelle. These results provide insights into the genetic mechanisms and evolutionary pressures that facilitated the metabolic integration of the host­endosymbiont association and sustained the evolution of a photosynthetic organelle.

Differential Phototactic Behavior of Closely Related Cyanobacterial Isolates from Yellowstone Hot Spring Biofilms.

Phototrophic biofilms in most environments experience major changes in light levels throughout a diel cycle. Phototaxis can be a useful strategy for optimizing light exposure under these conditions, but little is known about its role in cyanobacteria from thermal springs. We examined two closely related Synechococcus isolates (Synechococcus OS-A dominates at 60 to 65°C and OS-B' at 50 to 55°C) from outflows of Octopus Spring in Yellowstone National Park. Both isolates exhibited phototaxis and photokinesis in white light, but with differences in speed and motility bias. OS-B' exhibited phototaxis toward UVA, blue, green, and red wavelengths, while OS-A primarily exhibited phototaxis toward red and green. OS-A also exhibited negative phototaxis under certain conditions. The repertoires of photoreceptors and signal transduction elements in both isolates were quite different from those characterized in other unicellular cyanobacteria. These differences in the photoresponses between OS-A and OS-B' in conjunction with in situ observations indicate that phototactic strategies may be quite versatile and finely tuned to the light and local environment. IMPORTANCE Optimizing light absorption is of paramount importance to photosynthetic organisms. Some photosynthetic microbes have evolved a sophisticated process called phototaxis to move toward or away from a light source. In many hot springs in Yellowstone National Park, cyanobacteria thrive in thick, laminated biofilms or microbial mats, where small movements can result in large changes in light exposure. We quantified the light-dependent motility behaviors in isolates representing two of the most abundant and closely related cyanobacterial species from these springs. We found that they exhibited unexpected differences in their speed, directionality, and responses to different intensities or qualities of light. An examination of their genomes revealed several variations from well-studied phototaxis-related genes. Studying these recently isolated cyanobacteria reveals that diverse phototactic strategies can exist even among close relatives in the same environment. It also provides insights into the importance of phototaxis for growth and survival in microbial biofilm communities.

Deciphering proteolysis pathways for the error-prone DNA polymerase in cyanobacteria.

Protein quality control pathways require AAA+ proteases, such as Clp and Lon. Lon protease maintains UmuD, an important component of the error-prone DNA repair polymerase (Pol V), at very low levels in E. coli. Most members of the phylum Cyanobacteria lack Lon (including the model cyanobacterium, Synechocystis sp. PCC6803), so maintenance of UmuD at low levels must employ different proteases. We demonstrate that the first 19 residues from the N-terminus of UmuD (Sug1-19 ) fused to a reporter protein are adequate to trigger complete proteolysis and that mutation of a single leucine residue (L6) to aspartic acid inhibits proteolysis. This process appears to follow the N-end rule and is mediated by ClpA/P protease and the ClpS adaptor. Additionally, mutations of arginine residues in the Sug1-19 tag suggest that the ClpX/P pathway also plays a role in proteolysis. We propose that there is a dual degron at the N-terminus of the UmuD protein in Synechocystis sp. PCC6803, which is distinct from the degron required for degradation of UmuD in E. coli. The use of two proteolysis pathways to tune levels of UmuD might reflect how a photosynthetic organism responds to multiple environmental stressors.

The role of three-tandem Pho Boxes in the control of the C-P lyase operon in a thermophilic cyanobacterium.

Cyanobacteria have an inherited advantage in phosphonate phytoremediation. However, studies on phosphonate metabolism in cyanobacteria are rare and mostly focus on physiology and ecology. Here, C-P lyase gene cluster regulation in an undomesticated thermophilic Synechococcus OS-B' was examined in Synechocystis sp. PCC6803, a unicellular cyanobacterial model. Phylogenetic and cluster synteny analysis of C-P lyase genes revealed a closer relationship between Syn OS-B' and Thermus thermophilus, than with other cyanobacteria. Pho boxes were identified in the 5'-end-flanking region of the C-P lyase gene cluster, through which the downstream gene expression was regulated in a phosphate concentration-dependent manner. Unexpectedly, the phosphate concentration that thoroughly inhibited Pho boxes was almost two orders of magnitude higher than that of any natural or anthropogenic wastewater reported so far. The Pho boxes mediated regulation was achieved through the Pho regulon two-component system, and the absence of either SphS or SphR ablated the cell's ability to sense ambient phosphate changes. The three tandems of Pho boxes maintained inequivalent roles, of which the third tandem was not essential; however, it played a role in adjusting Pho boxes response in both positive and negative manner under phosphorus limitation.

CRISPR-Cas systems in bacteria and archaea: versatile small RNAs for adaptive defense and regulation.

Bacteria and archaea have evolved defense and regulatory mechanisms to cope with various environmental stressors, including virus attack. This arsenal has been expanded by the recent discovery of the versatile CRISPR-Cas system, which has two novel features. First, the host can specifically incorporate short sequences from invading genetic elements (virus or plasmid) into a region of its genome that is distinguished by clustered regularly interspaced short palindromic repeats (CRISPRs). Second, when these sequences are transcribed and precisely processed into small RNAs, they guide a multifunctional protein complex (Cas proteins) to recognize and cleave incoming foreign genetic material. This adaptive immunity system, which uses a library of small noncoding RNAs as a potent weapon against fast-evolving viruses, is also used as a regulatory system by the host. Exciting breakthroughs in understanding the mechanisms of the CRISPR-Cas system and its potential for biotechnological applications and understanding evolutionary dynamics are discussed.

Long-term microfluidic tracking of coccoid cyanobacterial cells reveals robust control of division timing.

BACKGROUND: Cyanobacteria are important agents in global carbon and nitrogen cycling and hold great promise for biotechnological applications. Model organisms such as Synechocystis sp. and Synechococcus sp. have advanced our understanding of photosynthetic capacity and circadian behavior, mostly using population-level measurements in which the behavior of individuals cannot be monitored. Synechocystis sp. cells are small and divide slowly, requiring long-term experiments to track single cells. Thus, the cumulative effects of drift over long periods can cause difficulties in monitoring and quantifying cell growth and division dynamics. RESULTS: To overcome this challenge, we enhanced a microfluidic cell-culture device and developed an image analysis pipeline for robust lineage reconstruction. This allowed simultaneous tracking of many cells over multiple generations, and revealed that cells expand exponentially throughout their cell cycle. Generation times were highly correlated for sister cells, but not between mother and daughter cells. Relationships between birth size, division size, and generation time indicated that cell-size control was inconsistent with the "sizer" rule, where division timing is based on cell size, or the "timer" rule, where division occurs after a fixed time interval. Instead, single cell growth statistics were most consistent with the "adder" rule, in which division occurs after a constant increment in cell volume. Cells exposed to light-dark cycles exhibited growth and division only during the light period; dark phases pause but do not disrupt cell-cycle control. CONCLUSIONS: Our analyses revealed that the "adder" model can explain both the growth-related statistics of single Synechocystis cells and the correlation between sister cell generation times. We also observed rapid phenotypic response to light-dark transitions at the single cell level, highlighting the critical role of light in cyanobacterial cell-cycle control. Our findings suggest that by monitoring the growth kinetics of individual cells we can build testable models of circadian control of the cell cycle in cyanobacteria.

Gordon research conference on the dynamics and regulation of photosynthesis: from the origin of bio-catalysis to innovative solar conversion.

We provide here a news report on the 2015 Gordon Research Conference "Dynamics and regulation of photosynthesis: from the origin of biocatalysis to innovative solar conversion.'' It was held at Bentley University, Waltham, MA, USA, June 28-July 3, 2015, and offered a mix of traditional and emerging areas that highlighted new directions and methods of analyses. A major innovation was short (1 min) poster highlights that added an exciting dynamic to the interactions. Following the end of the formal sessions, three young scientists (Andrian Gutu, of Harvard University, USA; Alizée Malnoë, of University of California, Berkeley, USA; and Yuval Mazor of Tel Aviv University, Israel) were recognized for their research; they also each received a recent volume of "Advances in photosynthesis and respiration including bioenergy and related processes" from Govindjee. We also provide at the end a brief report on the Gordon Research Seminar that preceded the conference.

Maintenance of motility bias during cyanobacterial phototaxis.

Signal transduction in bacteria is complex, ranging across scales from molecular signal detectors and effectors to cellular and community responses to stimuli. The unicellular, photosynthetic cyanobacterium Synechocystis sp. PCC6803 transduces a light stimulus into directional movement known as phototaxis. This response occurs via a biased random walk toward or away from a directional light source, which is sensed by intracellular photoreceptors and mediated by Type IV pili. It is unknown how quickly cells can respond to changes in the presence or directionality of light, or how photoreceptors affect single-cell motility behavior. In this study, we use time-lapse microscopy coupled with quantitative single-cell tracking to investigate the timescale of the cellular response to various light conditions and to characterize the contribution of the photoreceptor TaxD1 (PixJ1) to phototaxis. We first demonstrate that a community of cells exhibits both spatial and population heterogeneity in its phototactic response. We then show that individual cells respond within minutes to changes in light conditions, and that movement directionality is conferred only by the current light directionality, rather than by a long-term memory of previous conditions. Our measurements indicate that motility bias likely results from the polarization of pilus activity, yielding variable levels of movement in different directions. Experiments with a photoreceptor (taxD1) mutant suggest a supplementary role of TaxD1 in enhancing movement directionality, in addition to its previously identified role in promoting positive phototaxis. Motivated by the behavior of the taxD1 mutant, we demonstrate using a reaction-diffusion model that diffusion anisotropy is sufficient to produce the observed changes in the pattern of collective motility. Taken together, our results establish that single-cell tracking can be used to determine the factors that affect motility bias, which can then be coupled with biophysical simulations to connect changes in motility behaviors at the cellular scale with group dynamics.

Challenges of metagenomics and single-cell genomics approaches for exploring cyanobacterial diversity.

Cyanobacteria have played a crucial role in the history of early earth and continue to be instrumental in shaping our planet, yet applications of cutting edge technology have not yet been widely used to explore cyanobacterial diversity. To provide adequate background, we briefly review current sequencing technologies and their innovative uses in genomics and metagenomics. Next, we focus on current cell capture technologies and the challenges of using them with cyanobacteria. We illustrate the utility in coupling breakthroughs in DNA amplification with cell capture platforms, with an example of microfluidic isolation and subsequent targeted amplicon sequencing from individual terrestrial thermophilic cyanobacteria. Single cells of thermophilic, unicellular Synechococcus sp. JA-2-3-B'a(2-13) (Syn OS-B') were sorted in a microfluidic device, lysed, and subjected to whole genome amplification by multiple displacement amplification. We amplified regions from specific CRISPR spacer arrays, which are known to be highly diverse, contain semi-palindromic repeats which form secondary structure, and can be difficult to amplify. Cell capture, lysis, and genome amplification on a microfluidic device have been optimized, setting a stage for further investigations of individual cyanobacterial cells isolated directly from natural populations.

Motility enhancement through surface modification is sufficient for cyanobacterial community organization during phototaxis.

The emergent behaviors of communities of genotypically identical cells cannot be easily predicted from the behaviors of individual cells. In many cases, it is thought that direct cell-cell communication plays a critical role in the transition from individual to community behaviors. In the unicellular photosynthetic cyanobacterium Synechocystis sp. PCC 6803, individual cells exhibit light-directed motility ("phototaxis") over surfaces, resulting in the emergence of dynamic spatial organization of multicellular communities. To probe this striking community behavior, we carried out time-lapse video microscopy coupled with quantitative analysis of single-cell dynamics under varying light conditions. These analyses suggest that cells secrete an extracellular substance that modifies the physical properties of the substrate, leading to enhanced motility and the ability for groups of cells to passively guide one another. We developed a biophysical model that demonstrates that this form of indirect, surface-based communication is sufficient to create distinct motile groups whose shape, velocity, and dynamics qualitatively match our experimental observations, even in the absence of direct cellular interactions or changes in single-cell behavior. Our computational analysis of the predicted community behavior, across a matrix of cellular concentrations and light biases, demonstrates that spatial patterning follows robust scaling laws and provides a useful resource for the generation of testable hypotheses regarding phototactic behavior. In addition, we predict that degradation of the surface modification may account for the secondary patterns occasionally observed after the initial formation of a community structure. Taken together, our modeling and experiments provide a framework to show that the emergent spatial organization of phototactic communities requires modification of the substrate, and this form of surface-based communication could provide insight into the behavior of a wide array of biological communities.

Draft genome of Chloroflexus sp. MS-CIW-1, of the Chloroflexus sp. MS-G group from Mushroom Spring, Yellowstone National Park.

Chloroflexus sp. MS-CIW-1 was isolated from a phototrophic mat in Mushroom Spring, an alkaline hot spring in Yellowstone National Park, WY, USA. We report the draft genome of 4.8 Mb consisting of 6 contigs with 3755 protein-coding genes and a GC content of 54.45%.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs.

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

Role of polyphosphate in thermophilic Synechococcus sp. from microbial mats.

Synechococcus OS-B', a thermophilic unicellular cyanobacterium, recently isolated from the microbial mats in Octopus Spring (Yellowstone National Park), induces a suite of genes, including phosphatases and transporters, in response to phosphorus (P) starvation. Here we describe two different approaches to examine the ability of Synechococcus OS-B' to synthesize and break down polyphosphate (poly P), a key storage compound in many prokaryotes. First, we developed a transformation protocol to create mutants in the polyphosphate kinase (ppk), the major enzyme responsible for the synthesis of poly P. The ppk mutant exhibited a pleiotropic phenotype with defects in poly P accumulation, aberrant levels of Pho regulon transcripts, growth defects, and changes in cell size and exopolysaccharide levels, among others. Second, we measured transcripts of ppk and ppx (encoding the polyphosphatase) directly from mat samples and found that the levels varied dramatically over a diel cycle. We also used Western blot analysis to quantify levels of PPK and PPX and found that these enzymes differentially accumulated during the diel cycle. Levels of polyphosphate kinase peaked at night, while polyphosphatase levels were highest during the early morning hours. We hypothesize that the opposing activities of these two enzymes allow cells to store and utilize poly P to optimize growth over a diel cycle.

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population.

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

Novel miniature transposable elements in thermophilic Synechococcus strains and their impact on an environmental population.

The genomes of the two closely related freshwater thermophilic cyanobacteria Synechococcus sp. strain JA-3-3Ab and Synechococcus sp. strain JA-2-3B'a(2-13) each host several families of insertion sequences (ISSoc families) at various copy numbers, resulting in an overall high abundance of insertion sequences in the genomes. In addition to full-length copies, a large number of internal deletion variants have been identified. ISSoc2 has two variants (ISSoc2∂-1 and ISSoc2∂-2) that are observed to have multiple near-exact copies. Comparison of environmental metagenomic sequences to the Synechococcus genomes reveals novel placement of copies of ISSoc2, ISSoc2∂-1, and ISSoc2∂-2. Thus, ISSoc2∂-1 and ISSoc2∂-2 appear to be active nonautonomous mobile elements derived by internal deletion from ISSoc2. Insertion sites interrupting genes that are likely critical for cell viability were detected; however, most insertions either were intergenic or were within genes of unknown function. Most novel insertions detected in the metagenome were rare, suggesting a stringent selective environment. Evidence for mobility of internal deletion variants of other insertion sequences in these isolates suggests that this is a general mechanism for the formation of miniature insertion sequences.

Modeling local interactions during the motion of cyanobacteria.

Synechocystis sp., a common unicellular freshwater cyanobacterium, has been used as a model organism to study phototaxis, an ability to move in the direction of a light source. This microorganism displays a number of additional characteristics such as delayed motion, surface dependence, and a quasi-random motion, where cells move in a seemingly disordered fashion instead of in the direction of the light source, a global force on the system. These unexplained motions are thought to be modulated by local interactions between cells such as intercellular communication. In this paper, we consider only local interactions of these phototactic cells in order to mathematically model this quasi-random motion. We analyze an experimental data set to illustrate the presence of quasi-random motion and then derive a stochastic dynamic particle system modeling interacting phototactic cells. The simulations of our model are consistent with experimentally observed phototactic motion.

Analysis of insertion sequences in thermophilic cyanobacteria: exploring the mechanisms of establishing, maintaining, and withstanding high insertion sequence abundance.

Insertion sequences (ISs) are simple mobile genetic elements capable of relocating within a genome. Through this transposition activity, they are known to create mutations which are mostly deleterious to the cell, although occasionally they are beneficial. Two closely related isolates of thermophilic Synechococcus species from hot spring microbial mats are known to harbor a large number of diverse ISs. To explore the mechanism of IS acquisition within natural populations and survival in the face of high IS abundance, we examined IS content and location in natural populations of Synechococcus by comparing metagenomic data to the genomes of fully sequenced cultured isolates. The observed IS distribution in the metagenome was equivalent to the distribution in the isolates, indicating that the cultured isolates are appropriate models for the environmental population. High sequence conservation between IS families shared between the two isolates suggests that ISs are able to move between individuals within populations and between species via lateral gene transfer, consistent with models for IS family accumulation. Most IS families show evidence of recent activity, and interruption of critical genes in some individuals was observed, demonstrating that transposition is an ongoing mutational force in the populations.

Whole gene amplification and protein separation from a few cells.

Despite the growing interest to explore untapped microbial gene and protein diversity, no single platform has been able to acquire both gene and protein information from just a few cells. We present a microfluidic system that simultaneously performs on-chip capillary electrophoresis for protein analysis and whole genome amplification (WGA), and we demonstrate this by doing both for the same cohort of cyanobacterial cells. This technology opens avenues for studying protein profiles of precious environmental microbial samples and simultaneously accessing genomic information based on WGA.