Ancient balancing selection on heterocyst function in a cosmopolitan cyanobacterium

Genomic Legacies of Ancient Adaptation Illuminate GC-Content Evolution in Bacteria

Bacterial evolution is characterized by strong purifying selection as well as rapid adaptive evolution in changing environments. In this context, the genomic GC content (genomic GC) varies greatly but presents some level of phylogenetic stability, making it challenging to explain based on current hypotheses. To illuminate the evolutionary mechanisms of the genomic GC, we analyzed the base composition and functional inventory of 11,083 representative genomes. A phylogenetically constrained bimodal distribution of the genomic GC, which mainly originated from parallel divergences in the early evolution, was demonstrated. Such variation of the genomic GC can be well explained by DNA replication and repair (DRR), in which multiple pathways correlate with the genomic GC. Furthermore, the biased conservation of various stress-related genes, especially the DRR-related ones, implies distinct adaptive processes in the ancestral lineages of high- or low-GC clades which are likely induced by major environmental changes. Our findings support that the mutational biases resulting from these legacies of ancient adaptation have changed the course of adaptive evolution and generated great variation in the genomic GC. This highlights the importance of indirect effects of natural selection, which indicates a new model for bacterial evolution. IMPORTANCE GC content has been shown to be an important factor in microbial ecology and evolution, and the genomic GC of bacteria can be characterized by great intergenomic heterogeneity, high intragenomic homogeneity, and strong phylogenetic inertia, as well as being associated with the environment. Current hypotheses concerning direct selection or mutational biases cannot well explain these features simultaneously. Our findings of the genomic GC showing that ancient adaptations have transformed the DRR system and that the resulting mutational biases further contributed to a bimodal distribution of it offer a more reasonable scenario for the mechanism. This would imply that, when thinking about the evolution of life, diverse processes of adaptation exist, and combined effects of natural selection should be considered.

Compensatory Transcriptional Response of Fischerella thermalis to Thermal Damage of the Photosynthetic Electron Transfer Chain

Key organisms in the environment, such as oxygenic photosynthetic primary producers (photosynthetic eukaryotes and cyanobacteria), are responsible for fixing most of the carbon globally. However, they are affected by environmental conditions, such as temperature, which in turn affect their distribution. Globally, the cyanobacterium Fischerella thermalis is one of the main primary producers in terrestrial hot springs with thermal gradients up to 60 °C, but the mechanisms by which F. thermalis maintains its photosynthetic activity at these high temperatures are not known. In this study, we used molecular approaches and bioinformatics, in addition to photophysiological analyses, to determine the genetic activity associated with the energy metabolism of F. thermalis both in situ and in high-temperature (40 °C to 65 °C) cultures. Our results show that photosynthesis of F. thermalis decays with temperature, while increased transcriptional activity of genes encoding photosystem II reaction center proteins, such as PsbA (D1), could help overcome thermal damage at up to 60 °C. We observed that F. thermalis tends to lose copies of the standard G4 D1 isoform while maintaining the recently described D1INT isoform, suggesting a preference for photoresistant isoforms in response to the thermal gradient. The transcriptional activity and metabolic characteristics of F. thermalis, as measured by metatranscriptomics, further suggest that carbon metabolism occurs in parallel with photosynthesis, thereby assisting in energy acquisition under high temperatures at which other photosynthetic organisms cannot survive. This study reveals that, to cope with the harsh conditions of hot springs, F. thermalis has several compensatory adaptations, and provides emerging evidence for mixotrophic metabolism as being potentially relevant to the thermotolerance of this species. Ultimately, this work increases our knowledge about thermal adaptation strategies of cyanobacteria.

Genomic attributes of thermophilic and hyperthermophilic bacteria and archaea

Thermophiles and hyperthermophiles are immensely useful in understanding the evolution of life, besides their utility in environmental and industrial biotechnology. Advancements in sequencing technologies have revolutionized the field of microbial genomics. The massive generation of data enhances the sequencing coverage multi-fold and allows to analyse the entire genomic features of microbes efficiently and accurately. The mandate of a pure isolate can also be bypassed where whole metagenome-assembled genomes and single cell-based sequencing have fulfilled the majority of the criteria to decode various attributes of microbial genomes. A boom has, therefore, been seen in analysing the extremophilic bacteria and archaea using sequence-based approaches. Due to extensive sequence analysis, it becomes easier to understand the gene flow and their evolution among the members of bacteria and archaea. For instance, sequencing unveiled that Thermotoga maritima shares around 24% of genes of archaeal origin. Comparative and functional genomics provide an analytical view to understanding the microbial diversity of thermophilic bacteria and archaea, their interactions with other microbes, their adaptations, gene flow, and evolution over time. In this review, the genomic features of thermophilic bacteria and archaea are dealt with comprehensively.

BetaScan2: Standardized Statistics to Detect Balancing Selection Utilizing Substitution Data

Long-term balancing selection results in a build-up of alleles at similar frequencies and a deficit of substitutions when compared with an outgroup at a locus. The previously published β(1) statistics detect balancing selection using only polymorphism data. We now propose the β(2) statistic which detects balancing selection using both polymorphism and substitution data. In addition, we derive the variance of all β statistics, allowing for their standardization and thereby reducing the influence of parameters which can confound other selection tests. The standardized β statistics outperform existing summary statistics in simulations, indicating β is a well-powered and widely applicable approach for detecting balancing selection. We apply the β(2) statistic to 1000 Genomes data and report two missense mutations with high β scores in the ACSBG2 gene. An implementation of all β statistics and their standardization are available in the BetaScan2 software package at https://github.com/ksiewert/BetaScan.

Taxonomic Novelty and Distinctive Genomic Features of Hot Spring Cyanobacteria

Several cyanobacterial species are dominant primary producers in hot spring microbial mats. To date, hot spring cyanobacterial taxonomy, as well as the evolution of their genomic adaptations to high temperatures, are poorly understood, with genomic information currently available for only a few dominant genera, including Fischerella and Synechococcus. To address this knowledge gap, the present study expands the genomic landscape of hot spring cyanobacteria and traces the phylum-wide genomic consequences of evolution in high temperature environments. From 21 globally distributed hot spring metagenomes, with temperatures between 32 and 75°C, 57 medium- and high-quality cyanobacterial metagenome-assembled genomes were recovered, representing taxonomic novelty for 1 order, 3 families, 15 genera and 36 species. Comparative genomics of 93 hot spring genomes (including the 57 metagenome-assembled genomes) and 66 non-thermal genomes, showed that the former have smaller genomes and a higher GC content, as well as shorter proteins that are more hydrophilic and basic, when compared to the non-thermal genomes. Additionally, the core accessory orthogroups from the hot spring genomes of some genera had a greater abundance of functional categories, such as inorganic ion metabolism, translation and post-translational modifications. Moreover, hot spring genomes showed increased abundances of inorganic ion transport and amino acid metabolism, as well as less replication and transcription functions in the protein coding sequences. Furthermore, they showed a higher dependence on the CRISPR-Cas defense system against exogenous nucleic acids, and a reduction in secondary metabolism biosynthetic gene clusters. This suggests differences in the cyanobacterial response to environment-specific microbial communities. This phylum-wide study provides new insights into cyanobacterial genomic adaptations to a specific niche where they are dominant, which could be essential to trace bacterial evolution pathways in a warmer world, such as the current global warming scenario.

Cellular Innovation of the Cyanobacterial Heterocyst by the Adaptive Loss of Plasticity

Cellular innovation is central to biological diversification, yet its underlying mechanisms remain poorly understood [1]. One potential source of new cellular traits is environmentally induced phenotypic variation, or phenotypic plasticity. The plasticity-first hypothesis [2-4] proposes that natural selection can improve upon an ancestrally plastic phenotype to produce a locally adaptive trait, but the role of plasticity for adaptive evolution is still unclear [5-10]. Here, we show that a structurally novel form of the heterocyst, the specialized nitrogen-fixing cell of the multicellular cyanobacterium Fischerella thermalis, has evolved multiple times from ancestrally plastic developmental variation during adaptation to high temperature. Heterocyst glycolipids (HGs) provide an extracellular gas diffusion barrier that protects oxygen-sensitive nitrogenase [11, 12], and cyanobacteria typically exhibit temperature-induced plasticity in HG composition that modulates heterocyst permeability [13, 14]. By contrast, high-temperature specialists of F. thermalis constitutively overproduce glycolipid isomers associated with high temperature to levels unattained by plastic strains. This results in a less-permeable heterocyst, which is advantageous at high temperature but deleterious at low temperature for both nitrogen fixation activity and fitness. Our study illustrates how the origin of a novel cellular phenotype by the genetic assimilation and adaptive refinement of a plastic trait can be a source of biological diversity and contribute to ecological specialization.

Fischerella thermalis: a model organism to study thermophilic diazotrophy, photosynthesis and multicellularity in cyanobacteria

The true-branching cyanobacterium Fischerella thermalis (also known as Mastigocladus laminosus) is widely distributed in hot springs around the world. Morphologically, it has been described as early as 1837. However, its taxonomic placement remains controversial. F. thermalis belongs to the same genus as mesophilic Fischerella species but forms a monophyletic clade of thermophilic Fischerella strains and sequences from hot springs. Their recent divergence from freshwater or soil true-branching species and the ongoing process of specialization inside the thermal gradient make them an interesting evolutionary model to study. F. thermalis is one of the most complex prokaryotes. It forms a cellular network in which the main trichome and branches exchange metabolites and regulators via septal junctions. This species can adapt to a variety of environmental conditions, with its photosynthetic apparatus remaining active in a temperature range from 15 to 58 °C. Together with its nitrogen-fixing ability, this allows it to dominate in hot spring microbial mats and contribute significantly to the de novo carbon and nitrogen input. Here, we review the current knowledge on the taxonomy and distribution of F. thermalis, its morphological complexity, and its physiological adaptations to an extreme environment.

A genome-wide scan for genes under balancing selection in the plant pathogen Ralstonia solanacearum

Background

Plant pathogens are under significant selective pressure by the plant host. Consequently, they are expected to have adapted to this condition or contribute to evading plant defenses. In order to acquire long-term fitness, plant bacterial pathogens are usually forced to maintain advantageous genetic diversity in populations. This strategy ensures that different alleles in the pathogen’s gene pool are maintained in a population at frequencies larger than expected under neutral evolution. This selective process, known as balancing selection, is the subject of this work in the context of a common bacterial phytopathogen. We performed a genome-wide scan of Ralstonia solanacearum species complex, an aggressive plant bacterial pathogen that shows broad host range and causes a devastating disease called ‘bacterial wilt’.

Results

Using a sliding window approach, we analyzed 57 genomes from three phylotypes of the R. solanacearum species complex to detect signatures of balancing selection. A total of 161 windows showed extreme values in three summary statistics of population genetics: Tajima’s D, θ_w and Fu & Li’s D*. We discarded any confounding effects due to demographic events by means of coalescent simulations of genetic data. The prospective windows correspond to 78 genes with known function that map in any of the two main replicons (1.7% of total number of genes). The candidate genes under balancing selection are related to primary metabolism and other basal activities (51.3%) or directly associated to virulence (48.7%), the latter being involved in key functions targeted to dismantle plant defenses or to participate in critical stages in the pathogenic process.

Conclusions

We identified various genes under balancing selection that play a significant role in basic metabolism as well as in virulence of the R. solanacearum species complex. These genes are useful to understand and monitor the evolution of bacterial pathogen populations and emerge as potential candidates for future treatments to induce specific plant immune responses.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1456-6) contains supplementary material, which is available to authorized users.

Active Crossfire Between Cyanobacteria and Cyanophages in Phototrophic Mat Communities Within Hot Springs

Cyanophages are viruses with a wide distribution in aquatic ecosystems, that specifically infect Cyanobacteria. These viruses can be readily isolated from marine and fresh waters environments; however, their presence in cosmopolitan thermophilic phototrophic mats remains largely unknown. This study investigates the morphological diversity (TEM), taxonomic composition (metagenomics), and active infectivity (metatranscriptomics) of viral communities over a thermal gradient in hot spring phototrophic mats from Northern Patagonia (Chile). The mats were dominated (up to 53%) by cosmopolitan thermophilic filamentous true-branching cyanobacteria from the genus Mastigocladus, the associated viral community was predominantly composed of Caudovirales (70%), with most of the active infections driven by cyanophages (up to 90% of Caudovirales transcripts). Metagenomic assembly lead to the first full genome description of a T7-like Thermophilic Cyanophage recovered from a hot spring (Porcelana Hot Spring, Chile), with a temperature of 58°C (TC-CHP58). This could potentially represent a world-wide thermophilic lineage of podoviruses that infect cyanobacteria. In the hot spring, TC-CHP58 was active over a temperature gradient from 48 to 66°C, showing a high population variability represented by 1979 single nucleotide variants (SNVs). TC-CHP58 was associated to the Mastigocladus spp. by CRISPR spacers. Marked differences in metagenomic CRISPR loci number and spacers diversity, as well as SNVs, in the TC-CHP58 proto-spacers at different temperatures, reinforce the theory of co-evolution between natural virus populations and cyanobacterial hosts. Considering the importance of cyanobacteria in hot spring biogeochemical cycles, the description of this new cyanopodovirus lineage may have global implications for the functioning of these extreme ecosystems.

Temperature modulates Fischerella thermalis ecotypes in Porcelana Hot Spring

In the Porcelana Hot Spring (Northern Patagonia), true-branching cyanobacteria are the dominant primary producers in microbial mats, and they are mainly responsible for carbon and nitrogen fixation. However, little is known about their metabolic and genomic adaptations at high temperatures. Therefore, in this study, a total of 81 Fischerella thermalis strains (also known as Mastigocladus laminosus) were isolated from mat samples in a thermal gradient between 61-46°C. The complementary use of proteomic comparisons from these strains, and comparative genomics of F. thermalis pangenomes, suggested that at least two different ecotypes were present within these populations. MALDI-TOF MS analysis separated the strains into three clusters; two with strains obtained from mats within the upper temperature range (61 and 54°C), and a third obtained from mats within the lower temperature range (51 and 46°C). Both groups possessed different but synonymous nifH alleles. The main proteomic differences were associated with the abundance of photosynthesis-related proteins. Three F. thermalis metagenome assembled genomes (MAGs) were described from 66, 58 and 48°C metagenomes. These pangenomes indicated a divergence of orthologous genes and a high abundance of exclusive genes at 66°C. These results improved the current understanding of thermal adaptation of F. thermalis and the evolution of these thermophilic cyanobacterial species.