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.

Transcriptomes Divergence of Ricotia lunaria Between the Two Micro-Climatic Divergent Slopes at "Evolution Canyon" I, Israel

As one of the hotspot regions for sympatric speciation studies, Evolution Canyon (EC) became an ideal place for its high level of microclimatic divergence interslopes. In this study, to highlight the genetic mechanisms of sympatric speciation, phenotypic variation on flowering time and transcriptomic divergence were investigated between two ecotypes of Ricotia lunaria, which inhabit the opposite temperate and tropical slopes of EC I (Lower Nahal Oren, Mount Carmel, Israel) separated by 100 m at the bottom of the slopes. Growth chamber results showed that flowering time of the ecotype from south-facing slope population # 3 (SFS 3) was significantly 3 months ahead of the north-facing slope population # 5 (NFS 5). At the same floral development stage, transcriptome analysis showed that 1,064 unigenes were differentially expressed between the two ecotypes, which enriched in the four main pathways involved in abiotic and/or biotic stresses responses, including flavonoid biosynthesis, α-linolenic acid metabolism, plant-pathogen interaction and linoleic acid metabolism. Furthermore, based on Ka/Ks analysis, nine genes were suggested to be involved in the ecological divergence between the two ecotypes, whose homologs functioned in RNA editing, ABA signaling, photoprotective response, chloroplasts protein-conducting channel, and carbohydrate metabolism in Arabidopsis thaliana. Among them, four genes, namely, SPDS1, FCLY, Tic21 and BGLU25, also showed adaptive divergence between R. lunaria and A. thaliana, suggesting that these genes could play an important role in plant speciation, at least in Brassicaceae. Based on results of both the phenotype of flowering time and comparative transcriptome, we hypothesize that, after long-time local adaptations to their interslope microclimatic environments, the molecular functions of these nine genes could have been diverged between the two ecotypes. They might differentially regulate the expression of the downstream genes and pathways that are involved in the interslope abiotic stresses, which could further diverge the flowering time between the two ecotypes, and finally induce the reproductive isolation establishment by natural selection overruling interslope gene flow, promoting sympatric speciation.

Evolution in action: adaptation and incipient sympatric speciation with gene flow across life at “Evolution Canyon”, Israel

Various major evolutionary problems are still open, controversial or unsettled. These include even the basic evolutionary processes of adaptation and speciation. The “Evolution Canyon” model is a microscale natural laboratory that can highlight some of the basic problems requiring clarification (Nevo list of “Evolution Canyon” publications at http://evolution.haifa.ac.il). This is especially true if an interdisciplinary approach is practiced including ecological functional genomics, transcriptomics, proteomics, metabolomics and phenomics. Here I overview and reanalyze the incipient sympatric adaptive ecological speciation of five model organisms at “Evolution Canyon”, across life: the soil bacterium, Bacillus simplex; wild barley, the progenitor of cultivated barley, Hordeum spontaneum; the tiny beetle Oryzaephilus surinamensis; the cosmopolitan fruit-fly, Drosophila melanogaster, and the Africa-originated spiny mouse, Acomys cahirinus. All five models of organisms display evolution in action of microclimatic adaptation and incipient sympatric adaptive ecological speciation on the tropical and temperate abutting slopes, separated on average by only 250 meters. Some distant species converge in their micro-climatic adaptations to the hot and dry “African”, south-facing slope (SFS or AS) and to the cool and humid “European”, north-facing slope (NSF or ES). Natural selection overrules ongoing inter-slope gene-flow between the free interbreeding populations within and between slopes, and leads to adaptive incipient sympatric ecological speciation on the dramatically opposite abutting xeric savannoid and mesic forested slopes.

Information-theoretic uncertainty of SCFG-modeled folding space of the non-coding RNA

RNA secondary structure ensembles define probability distributions for alternative equilibrium secondary structures of an RNA sequence. Shannon's entropy is a measure for the amount of diversity present in any ensemble. In this work, Shannon's entropy of the SCFG ensemble on an RNA sequence is derived and implemented in polynomial time for both structurally ambiguous and unambiguous grammars. Micro RNA sequences generally have low folding entropy, as previously discovered. Surprisingly, signs of significantly high folding entropy were observed in certain ncRNA families. More effective models coupled with targeted randomization tests can lead to a better insight into folding features of these families.

AVAILABILITY

URL http://www.plantbio.uga.edu/~russell/index.php?s=1&n=5&r=0.

"Evolution Canyon," a potential microscale monitor of global warming across life

Climatic change and stress is a major driving force of evolution. The effects of climate change on living organisms have been shown primarily on regional and global scales. Here I propose the "Evolution Canyon" (EC) microscale model as a potential life monitor of global warming in Israel and the rest of the world. The EC model reveals evolution in action at a microscale involving biodiversity divergence, adaptation, and incipient sympatric speciation across life from viruses and bacteria through fungi, plants, and animals. The EC consists of two abutting slopes separated, on average, by 200 m. The tropical, xeric, savannoid, "African" south-facing slope (AS = SFS) abuts the forested "European" north-facing slope (ES = NFS). The AS receives 200-800% higher solar radiation than the ES. The ES represents the south European forested maquis. The AS and ES exhibit drought and shade stress, respectively. Major adaptations on the AS are because of solar radiation, heat, and drought, whereas those on the ES relate to light stress and photosynthesis. Preliminary evidence suggests the extinction of some European species on the ES and AS. In Drosophila, a 10-fold higher migration was recorded in 2003 from the AS to ES. I advance some predictions that could be followed in diverse species in EC. The EC microclimatic model is optimal to track global warming at a microscale across life from viruses and bacteria to mammals in Israel, and in additional ECs across the planet.