Comparative genomics and phylogenomic analyses of lysine riboswitch distributions in bacteria

Incipient Sympatric Speciation and Evolution of Soil Bacteria Revealed by Metagenomic and Structured Non-Coding RNAs Analysis

Simple Summary

The microevolutionary dynamics of soil bacteria under microclimatic differences are largely unexplored in contrast to our improving knowledge of their vast diversity. In this study, we performed a comparative metagenomic analysis of two sharply divergent rocks and soil types at the Evolution Plateau (EP) in eastern Upper Galilee, Israel. We have identified the significant differences in bacterial taxonomic diversity, functions, and patterns of RNA-based gene regulation between the bacteria from two different soil types. Furthermore, we have identified several species with a significant genetic divergence of the same species between the two soil types, highlighting the soil bacteria’s incipient sympatric speciation.

Abstract

Soil bacteria respond rapidly to changes in new environmental conditions. For adaptation to the new environment, they could mutate their genome, which impacts the alternation of the functional and regulatory landscape. Sometimes, these genetic and ecological changes may drive the bacterial evolution and sympatric speciation. Although sympatric speciation has been controversial since Darwin suggested it in 1859, there are several strong theoretical or empirical evidences to support it. Sympatric speciation associated with soil bacteria remains largely unexplored. Here, we provide potential evidence of sympatric speciation of soil bacteria by comparison of metagenomics from two sharply contrasting abutting divergence rock and soil types (Senonian chalk and its rendzina soil, and abutting Pleistocene basalt rock and basalt soil). We identified several bacterial species with significant genetic differences in the same species between the two soil types and ecologies. We show that the bacterial community composition has significantly diverged between the two soils; correspondingly, their functions were differentiated in order to adapt to the local ecological stresses. The ecologies, such as water availability and pH value, shaped the adaptation and speciation of soil bacteria revealed by the clear-cut genetic divergence. Furthermore, by a novel analysis scheme of riboswitches, we highlight significant differences in structured non-coding RNAs between the soil bacteria from two divergence soil types, which could be an important driver for functional adaptation. Our study provides new insight into the evolutionary divergence and incipient sympatric speciation of soil bacteria under microclimatic ecological differences.

The application of riboswitch sequencing for human gut bacterial classification and identification

Bacteria, especially gut bacteria play important roles in human health and diseases. The classification of many bacterial genera by the 16S ribosomal RNA (rRNA) has failed due to its low inter-species resolution. Given the wide distribution of riboswitches in bacteria, they may help 16S rRNA differentiate closely related species. We found that among 28 groups of species that could not be distinguished by 16S rRNA, eight of them could be separated by the TPP riboswitch and other riboswitches. Moreover, the species in the 16S rRNA database and these riboswitch databases overlap, therefore, using riboswitch databases can help 16S rRNA better identify species. In addition, we used Klenow DNA polymerase and a pair of short primers to facilitate the library construction of TPP riboswitches for sequencing. The sequencing results showed that the TPP riboswitch could detect the major phyla similar to those detected by 16S rRNA. Therefore, the TPP riboswitch and other riboswitch classes could potentially be applied to gut bacteria classification.

Computational analysis of sense-antisense chimeric transcripts reveals their potential regulatory features and the landscape of expression in human cells

Many human genes are transcribed from both strands and produce sense-antisense gene pairs. Sense-antisense (SAS) chimeric transcripts are produced upon the coalescing of exons/introns from both sense and antisense transcripts of the same gene. SAS chimera was first reported in prostate cancer cells. Subsequently, numerous SAS chimeras have been reported in the ChiTaRS-2.1 database. However, the landscape of their expression in human cells and functional aspects are still unknown. We found that longer palindromic sequences are a unique feature of SAS chimeras. Structural analysis indicates that a long hairpin-like structure formed by many consecutive Watson-Crick base pairs appears because of these long palindromic sequences, which possibly play a similar role as double-stranded RNA (dsRNA), interfering with gene expression. RNA-RNA interaction analysis suggested that SAS chimeras could significantly interact with their parental mRNAs, indicating their potential regulatory features. Here, 267 SAS chimeras were mapped in RNA-seq data from 16 healthy human tissues, revealing their expression in normal cells. Evolutionary analysis suggested the positive selection favoring sense-antisense fusions that significantly impacted the evolution of their function and structure. Overall, our study provides detailed insight into the expression landscape of SAS chimeras in human cells and identifies potential regulatory features.

Ancestral state reconstruction of metabolic pathways across pangenome ensembles

As genome sequencing efforts are unveiling the genetic diversity of the biosphere with an unprecedented speed, there is a need to accurately describe the structural and functional properties of groups of extant species whose genomes have been sequenced, as well as their inferred ancestors, at any given taxonomic level of their phylogeny. Elaborate approaches for the reconstruction of ancestral states at the sequence level have been developed, subsequently augmented by methods based on gene content. While these approaches of sequence or gene-content reconstruction have been successfully deployed, there has been less progress on the explicit inference of functional properties of ancestral genomes, in terms of metabolic pathways and other cellular processes. Herein, we describe PathTrace, an efficient algorithm for parsimony-based reconstructions of the evolutionary history of individual metabolic pathways, pivotal representations of key functional modules of cellular function. The algorithm is implemented as a five-step process through which pathways are represented as fuzzy vectors, where each enzyme is associated with a taxonomic conservation value derived from the phylogenetic profile of its protein sequence. The method is evaluated with a selected benchmark set of pathways against collections of genome sequences from key data resources. By deploying a pangenome-driven approach for pathway sets, we demonstrate that the inferred patterns are largely insensitive to noise, as opposed to gene-content reconstruction methods. In addition, the resulting reconstructions are closely correlated with the evolutionary distance of the taxa under study, suggesting that a diligent selection of target pangenomes is essential for maintaining cohesiveness of the method and consistency of the inference, serving as an internal control for an arbitrary selection of queries. The PathTrace method is a first step towards the large-scale analysis of metabolic pathway evolution and our deeper understanding of functional relationships reflected in emerging pangenome collections.

The Development of Bacteriophage Resistance in Vibrio alginolyticus Depends on a Complex Metabolic Adaptation Strategy

Lytic bacteriophages have been well documented to play a pivotal role in microbial ecology due to their complex interactions with bacterial species, especially in aquatic habitats. Although the use of phages as antimicrobial agents, known as phage therapy, in the aquatic environment has been increasing, recent research has revealed drawbacks due to the development of phage-resistant strains among Gram-negative species. Acquired phage resistance in marine Vibrios has been proven to be a very complicated process utilizing biochemical, metabolic, and molecular adaptation strategies. The results of our multi-omics approach, incorporating transcriptome and metabolome analyses of Vibrio alginolyticus phage-resistant strains, corroborate this prospect. Our results provide insights into phage-tolerant strains diminishing the expression of phage receptors ompF, lamB, and btuB. The same pattern was observed for genes encoding natural nutrient channels, such as rbsA, ptsG, tryP, livH, lysE, and hisp, meaning that the cell needs to readjust its biochemistry to achieve phage resistance. The results showed reprogramming of bacterial metabolism by transcript regulations in key-metabolic pathways, such as the tricarboxylic acid cycle (TCA) and lysine biosynthesis, as well as the content of intracellular metabolites belonging to processes that could also significantly affect the cell physiology. Finally, SNP analysis in resistant strains revealed no evidence of amino acid alterations in the studied putative bacterial phage receptors, but several SNPs were detected in genes involved in transcriptional regulation. This phenomenon appears to be a phage-specific, fine-tuned metabolic engineering, imposed by the different phage genera the bacteria have interacted with, updating the role of lytic phages in microbial marine ecology.

Unprecedented tunability of riboswitch structure and regulatory function by sub-millimolar variations in physiological Mg2

Riboswitches are cis-acting regulatory RNA biosensors that rival the efficiency of those found in proteins. At the heart of their regulatory function is the formation of a highly specific aptamer–ligand complex. Understanding how these RNAs recognize the ligand to regulate gene expression at physiological concentrations of Mg²⁺ ions and ligand is critical given their broad impact on bacterial gene expression and their potential as antibiotic targets. In this work, we used single-molecule FRET and biochemical techniques to demonstrate that Mg²⁺ ions act as fine-tuning elements of the amino acid-sensing lysC aptamer's ligand-free structure in the mesophile Bacillus subtilis. Mg²⁺ interactions with the aptamer produce encounter complexes with strikingly different sensitivities to the ligand in different, yet equally accessible, physiological ionic conditions. Our results demonstrate that the aptamer adapts its structure and folding landscape on a Mg²⁺-tunable scale to efficiently respond to changes in intracellular lysine of more than two orders of magnitude. The remarkable tunability of the lysC aptamer by sub-millimolar variations in the physiological concentration of Mg²⁺ ions suggests that some single-aptamer riboswitches have exploited the coupling of cellular levels of ligand and divalent metal ions to tightly control gene expression.

From Root to Tips: Sporulation Evolution and Specialization in Bacillus subtilis and the Intestinal Pathogen Clostridioides difficile

Bacteria of the Firmicutes phylum are able to enter a developmental pathway that culminates with the formation of highly resistant, dormant endospores. Endospores allow environmental persistence, dissemination and for pathogens, are also infection vehicles. In both the model Bacillus subtilis, an aerobic organism, and in the intestinal pathogen Clostridioides difficile, an obligate anaerobe, sporulation mobilizes hundreds of genes. Their expression is coordinated between the forespore and the mother cell, the two cells that participate in the process, and is kept in close register with the course of morphogenesis. The evolutionary mechanisms by which sporulation emerged and evolved in these two species, and more broadly across Firmicutes, remain largely unknown. Here, we trace the origin and evolution of sporulation using the genes known to be involved in the process in B. subtilis and C. difficile, and estimating their gain-loss dynamics in a comprehensive bacterial macroevolutionary framework. We show that sporulation evolution was driven by two major gene gain events, the first at the base of the Firmicutes and the second at the base of the B. subtilis group and within the Peptostreptococcaceae family, which includes C. difficile. We also show that early and late sporulation regulons have been coevolving and that sporulation genes entail greater innovation in B. subtilis with many Bacilli lineage-restricted genes. In contrast, C. difficile more often recruits new sporulation genes by horizontal gene transfer, which reflects both its highly mobile genome, the complexity of the gut microbiota, and an adjustment of sporulation to the gut ecosystem.

RiboD: a comprehensive database for prokaryotic riboswitches

Summary

Riboswitches are cis-regulatory non-coding genomic segments that control the expression of downstream genes by undergoing conformational change upon ligand binding. We present a comprehensive database of prokaryotic riboswitches that allows the user to search for riboswitches using multiple criteria, extract information about riboswitch location and gene/operon it regulates. RiboD provides a very useful resource that can be utilized for the better understanding of riboswitch-based gene regulation in bacteria and archaea.

Availability and implementation

RiboD can be freely accessed on the web at http://ribod.iiserkol.ac.in/.

Recent advances and future trends of riboswitches: attractive regulatory tools

Bacterial genomes contain a huge amount of different genes. These genes are spatiotemporally expressed to accomplish some required functions within the organism. Inside the cell, any step of gene expression may be modulated at four possible places such as transcription initiation, translation regulation, mRNA stability and protein stability. To achieve this, there is a necessity of strong regulators either natural or synthetic which can fine-tune gene expression regarding the required function. In recent years, riboswitches as metabolite responsive control elements residing in the untranslated regions of certain messenger RNAs, have been known to control gene expression at transcription or translation level. Importantly, these control elements do not prescribe the involvement of protein factors for metabolite binding. However, they own their particular properties to sense intramolecular metabolites (ligands). Herein, we highlighted current important bacterial riboswitches, their applications to support genetic control, ligand-binding domain mechanisms and current progress in synthetic riboswitches.

Phylogenomic and comparative analysis of the distribution and regulatory patterns of TPP riboswitches in fungi

Riboswitches are metabolite or ion sensing cis-regulatory elements that regulate the expression of the associated genes involved in biosynthesis or transport of the corresponding metabolite. Among the nearly 40 different classes of riboswitches discovered in bacteria so far, only the TPP riboswitch has also been found in algae, plants, and in fungi where their presence has been experimentally validated in a few instances. We analyzed all the available complete fungal and related genomes and identified TPP riboswitch-based regulation systems in 138 fungi and 15 oomycetes. We find that TPP riboswitches are most abundant in Ascomycota and Basidiomycota where they regulate TPP biosynthesis and/or transporter genes. Many of these transporter genes were found to contain conserved domains consistent with nucleoside, urea and amino acid transporter gene families. The genomic location of TPP riboswitches when correlated with the intron structure of the regulated genes enabled prediction of the precise regulation mechanism employed by each riboswitch. Our comprehensive analysis of TPP riboswitches in fungi provides insights about the phylogenomic distribution, regulatory patterns and functioning mechanisms of TPP riboswitches across diverse fungal species and provides a useful resource that will enhance the understanding of RNA-based gene regulation in eukaryotes.