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

Competition, Drought, Season Length? Disentangling Key Factors for Local Adaptation in Two Mediterranean Annuals across Combined Macroclimatic and Microclimatic Aridity Gradients

Competition in mesic sites and drought stress combined with short growing seasons in drier sites are key environmental factors along macroclimatic aridity gradients. They impose a triangular trade-off for local adaptation. However, as experiments have rarely disentangled their effects on plant fitness, uncertainty remained whether mesic populations are indeed better competitors and drier populations better adapted to drought stress and short season length. Aridity differs also at microclimatic scale between north (more mesic) and south (more arid) exposed hill-slopes. Little is known whether local adaptation occurs among exposures and whether south exposures harbor conspecifics better adapted to drier climates that could provide adaptive reservoirs under climate change. We sampled two Mediterranean annuals (Brachypodium hybridum, Hedypnois rhagadioloides) in 15 sites along a macroclimatic aridity gradient (89-926 mm rainfall) on corresponding north and south exposures. In a large greenhouse experiment, we measured their fitness under drought stress, competition, and short vs. long growing seasons. Along the macroclimatic gradient, mesic populations were better competitors under benign conditions. Drier populations performed no better under drought stress per se but coped better with the short growing seasons typical for drier macroclimates. At microclimatic scale, north exposure plants were slightly better competitors in H. rhagadioloides; in B. hybridum, south exposure plants coped better with drought under short season length. We demonstrate that local adaptation to drier macroclimates is trading-off with competitive ability under benign conditions and vice-versa. Drought escape via short life-cycles was the primary adaptation to drier macroclimates, suggesting that intensified drought stress within the growing season under climate change challenges arid and mesic populations alike. Moreover, the drier microclimates at south exposures exhibited some potential as nearby reservoirs of drier-adapted genotypes. This potential needs further investigation, yet may assist populations to persist under climate change and lessen the need for long-distance migration.

Fine-scale genetic structure and phenotypic divergence of a passerine bird population inhabiting a continuous Mediterranean woodland

Genetic differentiation between populations inhabiting ecologically different habitats might appear because of limited dispersal and gene flow, which may lead to patterns of phenotypic divergence and local adaptation. In this study, we use dispersal, genotypic (24 microsatellite loci) and phenotypic (body size and clutch size) data to analyse patterns of genetic structuring and phenotypic divergence in a blue tit (Cyanistes caeruleus) population inhabiting a continuous and heterogeneous woodland along a valley. The two slopes of the valley differ in their forest formations and environmental conditions. Findings showed that most blue tits reproduced within their natal slope. Accordingly, microsatellite analyses revealed that populations of blue tits established in the two slopes show subtle genetic differentiation. The two genetic populations diverged in clutch size, exceeding the level of differentiation expected based on genetic drift, hence suggesting divergent selection (or other processes promoting divergence) on this life-history trait. Our findings reveal that restricted dispersal and spatial heterogeneity may lead to genetic differentiation among bird populations at a surprisingly small scale. In this respect, it is worth highlighting that such differentiation occurs for an organism with high dispersal capacity and within a continuous woodland. Moreover, we show that small-scale ecological differences, together with limited gene flow, can result in selection favouring different phenotypes even within the same continuum population.

Plant root mechanisms and their effects on carbon and nutrient accumulation in desert ecosystems under changes in land use and climate

Deserts represent key carbon reservoirs, yet as these systems are threatened this has implications for biodiversity and climate change. This review focuses on how these changes affect desert ecosystems, particularly plant root systems and their impact on carbon and mineral nutrient stocks. Desert plants have diverse root architectures shaped by water acquisition strategies, affecting plant biomass and overall carbon and nutrient stocks. Climate change can disrupt desert plant communities, with droughts impacting both shallow and deep-rooted plants as groundwater levels fluctuate. Vegetation management practices, like grazing, significantly influence plant communities, soil composition, root microorganisms, biomass, and nutrient stocks. Shallow-rooted plants are particularly susceptible to climate change and human interference. To safeguard desert ecosystems, understanding root architecture and deep soil layers is crucial. Implementing strategic management practices such as reducing grazing pressure, maintaining moderate harvesting levels, and adopting moderate fertilization can help preserve plant-soil systems. Employing socio-ecological approaches for community restoration enhances carbon and nutrient retention, limits desert expansion, and reduces CO2 emissions. This review underscores the importance of investigating belowground plant processes and their role in shaping desert landscapes, emphasizing the urgent need for a comprehensive understanding of desert ecosystems.

The north-south divide? Macroalgal functional trait diversity and redundancy varies with intertidal aspect

BACKGROUND AND AIMS

Marine macroalgae ('seaweeds') are critical to coastal ecosystem structure and function, but also vulnerable to the many environmental changes associated with anthropogenic climate change (ACC). The local habitat conditions underpinning observed and predicted ACC-driven changes in intertidal macroalgal communities are complex and probably site-specific and operate in addition to more commonly reported regional factors such as sea surface temperatures.

METHODS

We examined how the composition and functional trait expression of macroalgal communities in SW England varied with aspect (i.e. north-south orientation) at four sites with opposing Equator- (EF) and Pole-facing (PF) surfaces. Previous work at these sites had established that average annual (low tide) temperatures vary by 1.6 °C and that EF-surfaces experience six-fold more frequent extremes (i.e. >30 °C).

KEY RESULTS

PF macroalgal communities were consistently more taxon rich; 11 taxa were unique to PF habitats, with only one restricted to EF. Likewise, functional richness and dispersion were greater on PF-surfaces (dominated by algae with traits linked to rapid resource capture and utilization, but low desiccation tolerance), although differences in both taxon and functional richness were probably driven by the fact that less diverse EF-surfaces were dominated by desiccation-tolerant fucoids.

CONCLUSIONS

Although we cannot disentangle the influence of temperature variation on algal ecophysiology from the indirect effects of aspect on species interactions (niche pre-emption, competition, grazing, etc.), our study system provides an excellent model for understanding how environmental variation at local scales affects community composition and functioning. By virtue of enhanced taxonomic diversity, PF-aspects supported higher functional diversity and, consequently, greater effective functional redundancy. These differences may imbue PF-aspects with resilience against environmental perturbation, but if predicted increases in global temperatures are realized, some PF-sites may shift to a depauperate, desiccation-tolerant seaweed community with a concomitant loss of functional diversity and redundancy.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

The potential for plant growth-promoting bacteria to impact crop productivity in future agricultural systems is linked to understanding the principles of microbial ecology

Global climate change poses challenges to land use worldwide, and we need to reconsider agricultural practices. While it is generally accepted that biodiversity can be used as a biomarker for healthy agroecosystems, we must specify what specifically composes a healthy microbiome. Therefore, understanding how holobionts function in native, harsh, and wild habitats and how rhizobacteria mediate plant and ecosystem biodiversity in the systems enables us to identify key factors for plant fitness. A systems approach to engineering microbial communities by connecting host phenotype adaptive traits would help us understand the increased fitness of holobionts supported by genetic diversity. Identification of genetic loci controlling the interaction of beneficial microbiomes will allow the integration of genomic design into crop breeding programs. Bacteria beneficial to plants have traditionally been conceived as "promoting and regulating plant growth". The future perspective for agroecosystems should be that microbiomes, via multiple cascades, define plant phenotypes and provide genetic variability for agroecosystems.

Genetic Biodiversity and Posttranslational Modifications of Protease Serine Endopeptidase in Different Strains of Sordaria fimicola

Genetic variations (mutation, crossing over, and recombination) act as a source for the gradual alternation in phenotype along a geographic transect where the environment changes. Posttranslational modifications (PTMs) predicted modifications successfully in different and the same species of living organisms. Protein diversity of living organisms is predicted by PTMs. Environmental stresses change nucleotides to produce alternations in protein structures, and these alternations have been examined through bioinformatics tools. The goal of the current study is to search the diversity of genes and posttranslational modifications of protease serine endopeptidase in various strains of Sordaria fimicola. The S. fimicola's genomic DNA was utilized to magnify the protease serine endopeptidase (SP2) gene; the size of the product was 700 and 1400 base pairs. Neurospora crassa was taken as the reference strain for studying the multiple sequence alignment of the nucleotide sequence. Six polymorphic sites of six strains of S. fimicola with respect to N. crassa were under observation. Different bioinformatics tools, i.e., NetPhos 3.1, NetNES 1.1 Server, YinOYang1.2, and Mod Pred, to search phosphorylation sites, acetylation, nuclear export signals, O-glycosylation, and methylation, respectively, were used to predict PTMs. The findings of the current study were 35 phosphorylation sites on the residues of serine for protease SP2 in SFS and NFS strains of S. fimicola and N. crassa. The current study supported us to get the reality of genes involved in protease production in experimental fungi. Our study examined the genetic biodiversity in six strains of S. fimicola which were caused by stressful environments, and these variations are a strong reason for evolution. In this manuscript, we predicted posttranslational modifications of protease serine endopeptidase in S. fimicola obtained from different sites, for the first time, to see the effect of environmental stress on nucleotides, amino acids, and proteases and to study PTMs by using various bioinformatics tools. This research confirmed the genetic biodiversity and PTMs in six strains of S. fimicola, and the designed primers also provided strong evidence for the presence of protease serine endopeptidase in each strain of S. fimicola.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

Molecular evidence for adaptive evolution of drought tolerance in wild cereals

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

How do ectotherms perform in cold environments? Physiological and life-history traits in an Andean viviparous lizard

Both the mean and the variation in environmental temperature are increasing globally. Indeed, the predicted increases in temperature range from 2 to 4°C in the next 50 years. Ectotherms control body temperature by means of behavior selecting microsites with different temperatures, which makes them more susceptible to changes in climate. Nevertheless, lizards living in high mountain environments have developed several mechanisms to inhabit and colonize variable environments with extreme temperatures. These mechanisms include a high metabolism to be active at lower temperatures and viviparity to improve embryonic development. Despite behavioral thermoregulation acting as a buffer to changes in environmental temperature, other traits such as life-history traits may be less flexible. Consequently, in an attempt to understand how lizards cope with harsh habitats, we evaluated some physiological traits and responses of females of Liolaemus bellii from two contrasting slope sites with differences in environmental temperature and humidity, but at the same altitude in the southern Andes range. We collected pregnant females from opposite slopes and maintained them until parturition in a common-garden experiment. Females from the south-facing slope (S-slope) had higher preferred body temperature (Tpref) values before and after parturition and exhibited higher daily energy expenditure before parturition. Nevertheless, no difference in Tpref was shown by their offspring, suggesting a developmental plastic response or adaptation to lower environmental temperature. For instance, the higher metabolism during pregnancy could be associated with a shorter activity period on the snowy S-slope. Additionally, females from the S-slope had larger kidneys and gave birth later than N-slope females, likely due to developmental plasticity or genetic differentiation. How fixed these traits are, in individuals from the contrasting slopes, will determine the response capacity of the L. bellii population to climate change.

The phenomenon of red and yellow autumn leaves: Hypotheses, agreements and disagreements

Yellow and red autumn leaves are typical of many temperate/boreal woody plants. Since the 19th century, it has been either considered the non-functional outcome of chlorophyll degradation that unmasks the pre-existing yellow and red pigments or that the de novo synthesis of red anthocyanins in autumn leaves indicated that it should have a physiological function, although it was commonly ignored. Defending free amino acids and various other resources released especially following the breakdown of the photosynthetic system, and mobilizing them for storage in other organs before leaf fall, is the cornerstone of both the physiological and anti-herbivory hypotheses about the functions of yellow and red autumn leaf colouration. The complicated phenomenon of conspicuous autumn leaf colouration has received significant attention since the year 2000, especially because ecologists started paying attention to its anti-herbivory potential. The obvious imperfection of the hypotheses put forth in several papers stimulated many other scientists. Hot debates among physiologists, among ecologists, and between physiologists and ecologists have been common since the year 2000, first because the various functions of yellow and red autumn leaf colouration are non-exclusive, and second because many scientists were trained to focus on a single subject. Here, I will review the debates, especially between the photoprotective and the anti-herbivory hypotheses, and describe both the progress in their understanding and the required progress.

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.

Topographic and vegetation drivers of thermal heterogeneity along the boreal-grassland transition zone in western Canada: Implications for climate change refugia

Climate change refugia are areas that are relatively buffered from contemporary climate change and may be important safe havens for wildlife and plants under anthropogenic climate change. Topographic variation is an important driver of thermal heterogeneity, but it is limited in relatively flat landscapes, such as the boreal plain and prairie regions of western Canada. Topographic variation within this region is mostly restricted to river valleys and hill systems, and their effects on local climates are not well documented. We sought to quantify thermal heterogeneity as a function of topography and vegetation cover within major valleys and hill systems across the boreal-grassland transition zone. Using iButton data loggers, we monitored local temperature at four hills and 12 river valley systems that comprised a wide range of habitats and ecosystems in Alberta, Canada (N = 240), between 2014 and 2020. We then modeled monthly temperature by season as a function of topography and different vegetation cover types using general linear mixed effect models. Summer maximum temperatures (T max) varied nearly 6°C across the elevation gradient sampled. Local summer mean (T mean) and maximum (T max) temperatures on steep, north-facing slopes (i.e., low levels of potential solar radiation) were up to 0.70°C and 2.90°C cooler than highly exposed areas, respectively. T max in incised valleys was between 0.26 and 0.28°C cooler than other landforms, whereas areas with greater terrain roughness experienced maximum temperatures that were up to 1.62°C cooler. We also found that forest cover buffered temperatures locally, with coniferous and mixedwood forests decreasing summer T mean from 0.23 to 0.72°C and increasing winter T min by up to 2°C, relative to non-forested areas. Spatial predictions of temperatures from iButton data loggers were similar to a gridded climate product (ClimateNA), but the difference between them increased with potential solar radiation, vegetation cover, and terrain roughness. Species that can track their climate niche may be able to compensate for regional climate warming through local migrations to cooler microsites. Topographic and vegetation characteristics that are related to cooler local climates should be considered in the evaluation of future climate change impacts and to identify potential refugia from climate change.

Microbiome: A Tool for Plant Stress Management in Future Production Systems

Climate change, due to the altered composition of the global atmosphere from the “greenhouse effect”, is one of the biggest challenges to agricultural production systems [...]

Conserved ecophysiology despite disparate microclimatic conditions in a gecko

Microscale differences in the habitats organisms occupy can influence selection regimes and promote intraspecific variation of traits. Temperature-dependent traits can be locally adapted to climatic conditions or be highly conserved and insensitive to directional selection under all but the most extreme regimes, and thus be similar across populations. The opposing slopes of Nahal Oren canyon in the Carmel Mountains, Israel, are strikingly different: the south-facing slope receives intensive solar radiation, is hot and supports mostly annual vegetation, whereas the north-facing slope is ~10°C cooler, more humid, and supports Mediterranean woodland. We examined whether these differences manifest in the thermal physiology of a common gecko species Ptyodactylus guttatus in controlled laboratory conditions. We predicted that geckos from the hotter south-facing slope would prefer higher temperatures, have faster gut passage times, lower metabolic and evaporative water loss rates, and start diel activity earlier compared with north-facing slope conspecifics. Contrary to these predictions, there were no differences between any of the ecophysiological traits in geckos from the opposing slopes. Nevertheless, our data showed that individuals from the north-facing slope were generally more active in earlier hours of the afternoon compared with south-facing individuals. We suggest that P. guttatus individuals disperse between the slopes and either gene-flow or behavioral plasticity deter local adaptation, resulting in similar physiological traits. Perhaps a stronger contrast in climatic conditions and a stronger barrier are needed to result in interpopulation divergence in temperature-dependent traits.

Facing up to climate change: Community composition varies with aspect and surface temperature in the rocky intertidal

Marine rocky intertidal organisms are amongst those most affected by climate change with regional distributional changes observed for many species. Although often ascribed to increased sea surface temperatures, precise assessment of the local habitat conditions underpinning observed and predicted changes in community assembly is lacking. Here we examine how aspect (i.e. north-south orientation) affects intertidal community composition and how rock surface temperatures and stress responses of two dominant grazer species (Patella spp.) elucidate emergent differences. We quantified year-round temperature variation and surveyed intertidal community composition on paired natural rock gullies with Equator- (EF) and Pole-facing (PF) surfaces. We also investigated variation in limpet (Patella spp.) reproductive phenology and osmotic stress. Average annual temperatures were 0.8 °C (1.6 °C at low tide) higher, with six-fold more frequent extremes (i.e. > 30 °C) on EF than PF surfaces. Intertidal community composition varied with aspect across trophic levels with greater overall species richness, abundance of primary producers and grazers on PF-surfaces, and greater barnacle abundance on EF-surfaces. Although species richness of organisms from different biogeographical origins ('Boreal' or 'Lusitanian') did not vary, the Lusitanian limpet Patella depressa exhibited earlier reproductive development on EF-surfaces and both limpet species exhibited greater thermal stress on EF-surfaces. We argue that our study system provides a good model for understanding how temperature variation at local scales can affect community composition, as well as ecophysiological and ecological responses to climate change and so better inform and predict regional range shifts over coming decades.

Changes in "natural antibiotic" metabolite composition during tetraploid wheat domestication

Gramineous plants protect their seeds from a variety of biotic stresses by producing toxic and deterrent secondary metabolites such as benzoxazinoids. It is unclear how the composition and abundance of these natural toxins has changed over the course of crop-plant domestication. To address this uncertainty, we characterized differences in metabolic levels of benzoxazinoids and their derivatives, between four lines of tetraploid wheat: wild emmer wheat (WEW), the direct progenitor of modern wheat; non-fragile domesticated emmer wheat (DEW), which was first domesticated about 11,000 years ago; the subsequently developed non-fragile and free-threshing durum landraces (LD); and modern durum (MD) varieties. Three-dimensional principal component analysis of mass spectrometry data of wheat metabolites showed with high resolution clear differences between metabolic profiles of WEW, DEW, and durum (LD + MD) and similarity in the metabolic profiles of the two durum lines (LD and MD) that is coherent with the phylogenetic relationship between the corresponding wheat lines. Moreover, our results indicated that some secondary metabolites involved in plant defense mechanisms became significantly more abundant during wheat domestication, while other defensive metabolites decreased or were lost. These metabolic changes reflect the beneficial or detrimental roles the corresponding metabolites might play during the domestication of three taxonomic subspecies of tetraploid wheat (Triticum turgidum).

Multi-Omics Analysis Reveals the Mechanism Underlying the Edaphic Adaptation in Wild Barley at Evolution Slope (Tabigha)

At the microsite "Evolution Slope", Tabigha, Israel, wild barley (Hordeum spontaneum) populations adapted to dry Terra Rossa soil, and its derivative abutting wild barley population adapted to moist and fungi-rich Basalt soil. However, the mechanisms underlying the edaphic adaptation remain elusive. Accordingly, whole genome bisulfite sequencing, RNA-sequencing, and metabolome analysis are performed on ten wild barley accessions inhabiting Terra Rossa and Basalt soil. A total of 121 433 differentially methylated regions (DMRs) and 10 478 DMR-genes are identified between the two wild barley populations. DMR-genes in CG context (CG-DMR-genes) are enriched in the pathways related with the fundamental processes, and DMR-genes in CHH context (CHH-DMR-genes) are mainly associated with defense response. Transcriptome and metabolome analysis reveal that the primary and secondary metabolisms are more active in Terra Rossa and Basalt wild barley populations, respectively. Multi-omics analysis indicate that sugar metabolism facilitates the adaptation of wild barley to dry Terra Rossa soil, whereas the enhancement of phenylpropanoid/phenolamide biosynthesis is beneficial for wild barley to inhabit moist and fungi pathogen-rich Basalt soil. The current results make a deep insight into edaphic adaptation of wild barley and provide elite genetic and epigenetic resources for developing barley with high abiotic stress tolerance.

Silica Particles Trigger the Exopolysaccharide Production of Harsh Environment Isolates of Growth-Promoting Rhizobacteria and Increase Their Ability to Enhance Wheat Biomass in Drought-Stressed Soils

In coming decades, drought is expected to expand globally owing to increased evaporation and reduced rainfall. Understanding, predicting, and controlling crop plants’ rhizosphere has the potential to manipulate its responses to environmental stress. Our plant growth-promoting rhizobacteria (PGPR) are isolated from a natural laboratory, ‘The Evolution Canyon’, Israel, (EC), from the wild progenitors of cereals, where they have been co-habituating with their hosts for long periods of time. The study revealed that commercial TM50 silica particles (SN) triggered the PGPR production of exopolysaccharides (EPS) containing D-glucuronate (D-GA). The increased EPS content increased the PGPR water-holding capacity (WHC) and osmotic pressure of the biofilm matrix, which led to enhanced plant biomass in drought-stressed growth environments. Light- and cryo-electron- microscopic studies showed that, in the presence of silica (SN) particles, bacterial morphology is changed, indicating that SNs are associated with significant reprogramming in bacteria. The findings encourage the development of large-scale methods for isolate formulation with natural silicas that ensure higher WHC and hyperosmolarity under field conditions. Osmotic pressure involvement of holobiont cohabitation is also discussed.

Effects of Slope Aspect and Rainfall on Belowground Deep Fine Root Traits and Aboveground Tree Height

The vertical root distribution and rooting depth are the main belowground plant functional traits used to indicate drought resistance in arid and semiarid regions. The effects of the slope aspect on the aboveground traits are visible but not the belowground deep root traits. We aimed to investigate the fine root traits of the locust tree (Robinia pseudoacacia L.) planted on southerly and northerly aspects, and the variations in the rooting depth in regions with different rainfall, as well as assessing how deep rooting, might affect the response to drought in a loess region. We selected three study sites with different rainfall amounts, with six sampling plots at each site (three each with southerly and northerly aspects). Soil core samples were collected down to the depth where no roots were present. The locust trees tended to develop deeper fine roots rather than greater heights. The tree height and diameter were greater for locust trees on northerly aspects, whereas trees on southerly aspects had significantly deeper rooting depths. Fine root traits (root length, root area, and root dry weight density) were higher in the southerly aspect for both Changwu and Ansai, but lower in Suide. The ratio of the root front depth tree height ranged from 1.04 to 3.17, which was higher on southerly than northerly aspects, and it increased as the rainfall decreased. Locust tree growth traits (belowground fine root and aboveground tree height) were positively correlated with the mean annual rainfall. The soil moisture content of the topsoil decreased as the rainfall decreased, but the pattern varied in the deep layer. Our results suggest that the variations in the belowground rooting depth under different slope aspects may be related to plant survival strategies. The vertical extension of the rooting depth and tree height may be key functional traits that determine plant growth in drought-prone regions.

Cross-genera amplification of Cajanus spp. specific SSR markers in Clitoria ternatea (L.) and their application in genetic diversity studies

Clitoria ternatea (L.) is a medicinal leguminous plant and is cultivated to cater the need of herbal industries and asthetic purposes. The unavailability of steady molecular marker impedes the genetic improvement of C. ternatea. In the present study, transferability of 98 pairs of Cajanus spp. specific SSR primers were assessed among 14 genotypes of C. ternatea, varied for their flower color, floral architecture and bio-metabolite (taraxerol and delphinidin) content, and out of them 43 had successfully amplified the fragments. Among them, 36 pairs of primers showed 100% transferability, whereas rest seven varied from 42.86 to 92.85% transferability. The transferable 43 pairs of SSR primers generated 196 alleles across the 14 genotypes and the AMOVA analysis showed moderate genetic variation (55.1%) among the genotypes of C. ternatea, which was also reinforced by Nei's genetic distance and gene identity estimates derived haplotype matrix. Similarly, both the principal coordinate analysis and dendrogram grouped these 14 genotypes of C. ternatea into two major clusters based on SSR allele distribution and frequency, and the clustering pattern is in accordance with petal color but in contrast to floral architecture. MCheza based outlier analysis revealed 16 alleles for balancing selection, which are putatively involved in the maintenance of genetic polymorphism in C. ternatea. Moreover, the estimates of molecular diversity and bio-metabolite content revealed the possible use of these genotypes in future breeding programme of this species.

Current Progress in Understanding and Recovering the Wheat Genes Lost in Evolution and Domestication

The modern cultivated wheat has passed a long evolution involving origin of wild emmer (WEM), development of cultivated emmer, formation of spelt wheat and finally establishment of modern bread wheat and durum wheat. During this evolutionary process, rapid alterations and sporadic changes in wheat genome took place, due to hybridization, polyploidization, domestication, and mutation. This has resulted in some modifications and a high level of gene loss. As a result, the modern cultivated wheat does not contain all genes of their progenitors. These lost genes are novel for modern wheat improvement. Exploring wild progenitor for genetic variation of important traits is directly beneficial for wheat breeding. WEM wheat (Triticum dicoccoides) is a great genetic resource with huge diversity for traits. Few genes and quantitative trait loci (QTL) for agronomic, quantitative, biotic and abiotic stress-related traits have already been mapped from WEM. This resource can be utilized for modern wheat improvement by integrating identified genes or QTLs through breeding.

Soil Microbes Trade-Off Biogeochemical Cycling for Stress Tolerance Traits in Response to Year-Round Climate Change

Winter air temperatures are rising faster than summer air temperatures in high-latitude forests, increasing the frequency of soil freeze/thaw events in winter. To determine how climate warming and soil freeze/thaw cycles affect soil microbial communities and the ecosystem processes they drive, we leveraged the Climate Change across Seasons Experiment (CCASE) at the Hubbard Brook Experimental Forest in the northeastern United States, where replicate field plots receive one of three climate treatments: warming (+5°C above ambient in the growing season), warming in the growing season + winter freeze/thaw cycles (+5°C above ambient +4 freeze/thaw cycles during winter), and no treatment. Soil samples were taken from plots at six time points throughout the growing season and subjected to amplicon (rDNA) and metagenome sequencing. We found that soil fungal and bacterial community composition were affected by changes in soil temperature, where the taxonomic composition of microbial communities shifted more with the combination of growing-season warming and increased frequency of soil freeze/thaw cycles in winter than with warming alone. Warming increased the relative abundance of brown rot fungi and plant pathogens but decreased that of arbuscular mycorrhizal fungi, all of which recovered under combined growing-season warming and soil freeze/thaw cycles in winter. The abundance of animal parasites increased significantly under combined warming and freeze/thaw cycles. We also found that warming and soil freeze/thaw cycles suppressed bacterial taxa with the genetic potential for carbon (i.e., cellulose) decomposition and soil nitrogen cycling, such as N fixation and the final steps of denitrification. These new soil communities had higher genetic capacity for stress tolerance and lower genetic capacity to grow or reproduce, relative to the communities exposed to warming in the growing season alone. Our observations suggest that initial suppression of biogeochemical cycling with year-round climate change may be linked to the emergence of taxa that trade-off growth for stress tolerance traits.

Long Tandem Arrays of Cassandra Retroelements and Their Role in Genome Dynamics in Plants

Retrotransposable elements are widely distributed and diverse in eukaryotes. Their copy number increases through reverse-transcription-mediated propagation, while they can be lost through recombinational processes, generating genomic rearrangements. We previously identified extensive structurally uniform retrotransposon groups in which no member contains the gag, pol, or env internal domains. Because of the lack of protein-coding capacity, these groups are non-autonomous in replication, even if transcriptionally active. The Cassandra element belongs to the non-autonomous group called terminal-repeat retrotransposons in miniature (TRIM). It carries 5S RNA sequences with conserved RNA polymerase (pol) III promoters and terminators in its long terminal repeats (LTRs). Here, we identified multiple extended tandem arrays of Cassandra retrotransposons within different plant species, including ferns. At least 12 copies of repeated LTRs (as the tandem unit) and internal domain (as a spacer), giving a pattern that resembles the cellular 5S rRNA genes, were identified. A cytogenetic analysis revealed the specific chromosomal pattern of the Cassandra retrotransposon with prominent clustering at and around 5S rDNA loci. The secondary structure of the Cassandra retroelement RNA is predicted to form super-loops, in which the two LTRs are complementary to each other and can initiate local recombination, leading to the tandem arrays of Cassandra elements. The array structures are conserved for Cassandra retroelements of different species. We speculate that recombination events similar to those of 5S rRNA genes may explain the wide variation in Cassandra copy number. Likewise, the organization of 5S rRNA gene sequences is very variable in flowering plants; part of what is taken for 5S gene copy variation may be variation in Cassandra number. The role of the Cassandra 5S sequences remains to be established.

Small spaces, big impacts: contributions of micro-environmental variation to population persistence under climate change

Individuals within natural populations can experience very different abiotic and biotic conditions across small spatial scales owing to microtopography and other micro-environmental gradients. Ecological and evolutionary studies often ignore the effects of micro-environment on plant population and community dynamics. Here, we explore the extent to which fine-grained variation in abiotic and biotic conditions contributes to within-population variation in trait expression and genetic diversity in natural plant populations. Furthermore, we consider whether benign microhabitats could buffer local populations of some plant species from abiotic stresses imposed by rapid anthropogenic climate change. If microrefugia sustain local populations and communities in the short term, other eco-evolutionary processes, such as gene flow and adaptation, could enhance population stability in the longer term. We caution, however, that local populations may still decline in size as they contract into rare microhabitats and microrefugia. We encourage future research that explicitly examines the role of the micro-environment in maintaining genetic variation within local populations, favouring the evolution of phenotypic plasticity at local scales and enhancing population persistence under global change.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

Aggression and courtship differences found in Drosophila melanogaster from two different microclimates at Evolution Canyon, Israel

Aggression and courtship behavior were examined of wild Drosophila melanogaster flies isolated from two contrasting microclimates found at Evolution Canyon in Mt. Carmel, Israel: an African-like dry tropical Slope (AS) and a European-like humid temperate Slope (ES), separated by 250 meters. Studies were carried out to ask whether behavioral differences existed between the two populations obtained from opposite slopes with divergent microclimates in Israel. First, we measured and compared intraslope aggression between same sex fly pairings collected from the same slope. Both male and female flies displayed similar fighting abilities from both slopes. ES males, however, from the humid biome, showed a tendency to lunge more per aggressive encounter, compared with AS males from the dry biome. Next, we tested interslope aggression by pairing flies from opposite slopes. ES males displayed higher numbers of lunges, and won more fights against their AS opponents. We also observed enhanced courtship performances in ES compared to AS males. The fighting and courtship superiority seen in ES males could reinforce fitness and pre-mating reproductive isolation mechanisms that underlie incipient sympatric speciation. This may support an evolutionary advantage of adaptively divergent fruit fly aggression phenotypes from different environments.

Paenibacillus polymyxa biofilm polysaccharides antagonise Fusarium graminearum

Fusarium Head Blight (FHB) caused by Fusarium graminearum pathogens constitutes a major threat to agricultural production because it frequently reduces the yield and quality of the crop. The disease severity is predicted to increase in various regions owing to climate change. Integrated management where biocontrol plays an important role has been suggested in order to fight FHB. P. polymyxa A26 is known to be an effective antagonist against F. graminearum. Deeper understanding of the mode of action of P. polymyxa A26 is needed to develop strategies for its application under natural settings in order to effectively overcome the pathogenic effects. This study aims to re-evaluate a former study and reveal whether compounds other than non-ribosomal antibiotic lipopeptides could be responsible for the antagonistic effect, despite what is often reported. Wheat seedlings were grown to maturity and the spikes infected with the pathogen under greenhouse conditions. The development of FHB infection, quantified via the disease incidence severity and 100-kernel weight, was strongly correlated (r > 0.78, p < 0.01) with the content of the polysaccharide component D-glucuronic acid in the biofilm. Furthermore, while increased inoculum density from 10⁶ to 10⁸ cells/ml did not affect wild type performance, a significant increase was observed with the P. polymyxa mutant deficient in nonribosomal lipopeptide synthesis. Our results show that P. polymyxa A26 biofilm extracellular polysaccharides are capable of antagonizing F. graminearum and that the uronate content of the polysaccharides is of critical importance in the antagonism.

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.

Transposable elements, placental development, and oocyte activation: Cellular stress and AMPK links jumping genes with the creation of human life

Transposable elements (TEs), also known as "jumping genes", are DNA sequences first described by Nobel laureate Barbara McClintock that comprise nearly half of the human genome and are able to transpose or move from one genomic location to another. As McClintock also noted that a genome "shock" or stress may induce TE activation and transposition, accumulating evidence suggests that cellular stress (e.g. mediated by increases in intracellular reactive oxygen species [ROS] and calcium [Ca2+], etc.) induces TE mobilization in several model organisms and L1s (a member of the retrotransposon class of TEs) are active and capable of retrotransposition in human oocytes, human sperm, and in human neural progenitor cells. Cellular stress also plays a critical role in human placental development, with cytotrophoblast (CTB) differentiation leading to the formation of the syncytiotrophoblast (STB), a cellular layer that facilitates nutrient and gas exchange between the mother and the fetus. Syncytin-1, a protein that promotes fusion of CTB cells and is necessary for STB formation, and its receptor is found in human sperm and human oocytes, respectively, and increases in ROS and Ca2+ promote trophoblast differentiation and syncytin-1 expression. Cellular stress is also essential in promoting human oocyte maturation and activation which, similar to TE mobilization, can be induced by compounds that increase intracellular Ca2+ and ROS levels. AMPK is a master metabolic regulator activated by increases in ROS, Ca2+, and/or an AMP(ADP)/ATP ratio increase, etc. as well as compounds that induce L1 mobilization in human cells. AMPK knockdown inhibits trophoblast differentiation and AMPK-activating compounds that promote L1 mobility also enhance trophoblast differentiation. Cellular stressors that induce TE mobilization (e.g. heat shock) also promote oocyte maturation in an AMPK-dependent manner and the antibiotic ionomycin activates AMPK, promotes TE activation, and induces human oocyte activation, producing normal, healthy children. Metformin promotes AMPK-dependent telomerase activation (critical for telomere maintenance) and induces activation of the endonuclease RAG1 (promotes DNA cleavage and transposition) via AMPK. Both RAG1 and telomerase are derived from TEs. It is our hypothesis that cellular stress and AMPK links TE activation and transposition with placental development and oocyte activation, facilitating both human genome evolution and the creation of all human life. We also propose the novel observation that various cellular stress-inducing compounds (e.g. metformin, resveratrol, etc.) may facilitate beneficial TE activation and transposition and enhance fertilization and embryological development through a common mechanism of AMPK activation.

Genomic adaptation to drought in wild barley is driven by edaphic natural selection at the Tabigha Evolution Slope

Ecological divergence at a microsite suggests adaptive evolution, and this study examined two abutting wild barley populations, each 100 m across, differentially adapted to drought tolerance on two contrasting soil types, Terra Rossa and basalt at the Tabigha Evolution Slope, Israel. We resequenced the genomes of seven and six wild barley genotypes inhabiting the Terra Rossa and basalt soils, respectively, and identified a total of 69,192,653 single-nucleotide variants (SNVs) and insertions/deletions in comparison with a reference barley genome. Comparative genomic analysis between these abutting wild barley populations involved 19,615,087 high-quality SNVs. The results revealed dramatically different selection sweep regions relevant to drought tolerance driven by edaphic natural selection within 2,577 selected genes in these regions, including key drought-responsive genes associated with ABA synthesis and degradation (such as Cytochrome P450 protein) and ABA receptor complex (such as PYL2, SNF1-related kinase). The genetic diversity of the wild barley population inhabiting Terra Rossa soil is much higher than that from the basalt soil. Additionally, we identified different sets of genes for drought adaptation in the wild barley populations from Terra Rossa soil and from wild barley populations from Evolution Canyon I at Mount Carmel. These genes are associated with abscisic acid signaling, signaling and metabolism of reactive oxygen species, detoxification and antioxidative systems, rapid osmotic adjustment, and deep root morphology. The unique mechanisms for drought adaptation of the wild barley from the Tabigha Evolution Slope may be useful for crop improvement, particularly for breeding of barley cultivars with high drought tolerance.

Adaptive evolution of benzoxazinoids in wild emmer wheat, Triticum dicoccoides, at "Evolution Canyon", Mount Carmel, Israel

Background

"Evolution Canyon" (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unraveling evolution-in-action, highlighting the evolutionary processes of biodiversity evolution, adaptation, and incipient sympatric speciation. A major model organism in ECI is the tetraploid wild emmer wheat, Triticum dicoccoides (TD), the progenitor of cultivated emmer and durum wheat. TD displays dramatic interslope adaptive evolutionary divergence on the tropical, savannoid-hot and dry south-facing, "African" slope (AS), and on the temperate, forested, cool and humid, north-facing, "European" slope (ES), separated on average by 250 m. From the perspective of chemical evolution and metabolomics, it is important to unravel interslope divergence in biologically relevant secondary metabolites between the abutting slope populations. Here, in TD we examined hydroxamic acid (Hx), which is a family of secondary cereal metabolites, and plays a major role in defending the plant against fungi, insects and weeds.

Results

Our examination revealed that higher concentrations of DIBOA and DIMBOA were found in seedlings growing in the same greenhouse from seeds collected from the cool and humid forested ES, whereas the seedlings of seeds collected from the savannoid AS (both in root and shoot tissues), showed no DIMBOA. Remarkably, only DIBOA appears in both shoots and roots of the AS seedlings. It rises to a peak and then decreases in both organs and in seedlings from both slopes. The DIMBOA, which appears only in the ES seedlings, rises to a peak and decreases in the shoot, but increased and remained in a plateau in the root, till the end of the experiment.

Conculsions/Significance

The results suggest stronger genetic resistance of defense compounds DIBOA and DIMBOA against biotic stresses (fungi and other pathogens) by ES seedlings. However, AS seedlings responded earlier but were to the same biotic stresses. The genetic difference found in AS seedlings was caused by the main adaptive selection in AS, which was against climatic, abiotic stresses, and was weaker, or not at all, against biotic stresses. The distinct genetic interslope differences appear important and is very significant and are elaborated in the discussion.

Regulation of gene expression and RNA editing in Drosophila adapting to divergent microclimates

Determining the mechanisms by which a species adapts to its environment is a key endeavor in the study of evolution. In particular, relatively little is known about how transcriptional processes are fine-tuned to adjust to different environmental conditions. Here we study Drosophila melanogaster from 'Evolution Canyon' in Israel, which consists of two opposing slopes with divergent microclimates. We identify several hundred differentially expressed genes and dozens of differentially edited sites between flies from each slope, correlate these changes with genetic differences, and use CRISPR mutagenesis to validate that an intronic SNP in prominin regulates its editing levels. We also demonstrate that while temperature affects editing levels at more sites than genetic differences, genetically regulated sites tend to be less affected by temperature. This work shows the extent to which gene expression and RNA editing differ between flies from different microclimates, and provides insights into the regulation responsible for these differences.

Lack of population differentiation patterns of previously identified putatively adaptive transposable element insertions at microgeographic scales

Background

Transposable elements (TEs) play an important role in genome function and evolution. It has been shown that TEs are a considerable source of adaptive changes in the genome of Drosophila melanogaster. Specifically, footprints of selection at the DNA level, the presence of population differentiation patterns across environmental gradients, and detailed mechanistic and fitness analyses of a few candidate adaptive TEs pointed to the role of TEs in environmental adaptation. However, whether the population differentiation patterns observed at large geographic scales can be replicated at a microgeographic scale has never been assessed before.

Results

In this work, we explored the population patterns of putatively adaptive TEs at a micro-spatial scale level. We compared the frequencies of TEs, previously identified as putatively adaptive and putatively neutral, in populations collected in opposite slopes of the Evolution Canyon at Mt. Carmel in Israel separated by 200 m on average. However, the differentiation patterns previously observed across large geographic distances (2000–2200 km) were not replicated at the microscale level of the Evolution Canyon populations.

Conclusion

TE insertions previously associated with D. melanogaster environmental adaptation at a macro scale level do not play such a role at the microscale level of the Evolution Canyon populations. However, these results do not exclude a role of TEs in microgeographic adaptation because the dataset analyzed in this work is restricted to TEs identified in a single North American strain and as such is highly biased and incomplete.

Reviewers

This article was reviewed by Eugene Koonin, Limsoon Wong and Fyodor Kondrashov.

Putative adaptive inter-slope divergence of transposon frequency in fruit flies (Drosophila melanogaster) at "Evolution Canyon", Mount Carmel, Israel

Background

The current analysis of transposon elements (TE) in Drosophila melanogaster at Evolution Canyon, (EC), Israel, is based on data and analysis done by our collaborators (Drs. J. Gonzalez, J. Martinez and W. Makalowski, this issue). They estimated the frequencies of 28 TEs (transposon elements) in fruit flies (D. melanogaster) from the ecologically tropic, hot, and dry south-facing slope (SFS) or “African” slope (AS) of EC and compared it with the TE frequencies on the temperate-cool and humid north-facing slope (NFS) or “European” slope (ES), separated, on average, by 250 m. The flies were sampled from two stations on each slope. We received their results, including the frequencies of each TE on each slope, and the probabilities of the statistical analyses (G-tests) of each TE separately. We continued the analysis of the inter-slope differences of the frequencies of the TEs, and based our different conclusions on that analysis and on the difference between micro (=EC) and macro (2000 km.) comparisons [Gonzalez et al. 2015 doi:10.1186/s13062-015-0075-4].

Results

Our collaborators based all their conclusions on the non-significant results of each of the individual tests of the 28 TEs. We analysed also the distribution of the TE differences between the slopes, based on their results. Thirteen TEs were more frequent on the SFS, 11 were more frequent on the NFS, and four had equal frequencies. Because of the equalizing effect of the ongoing migration, only small and temporary differences between the slopes (0 – 0.06) were regarded by us as random fluctuations (drift). Three TEs were intermediate (0.08-0.09) and await additional research. The 11 TEs with large frequency differences (0.12 – 0.22) were regarded by us as putative adaptive TEs, because the equalizing power of ongoing migration will eliminate random large differences. Five of them were higher on the SFS and six were higher on the NFS. Gaps in the distribution of the differences distinguished between the large and small differences. The large gap among the 11 TEs favored on the NFS was significant and supports our rejection of drift as the only explanation of the distribution of the slope differences. The gaps in the distribution of the differences separated the putative TEs with strong enough selection from those TEs that couldn't overrule the migration.

The results are compared and contrasted with the directional effect of the frequencies of the same TEs in the study of global climatic comparisons across thousands of kilometers. From the 11 putative adaptive TEs in the local “Evolution Canyon,” six differentiate in the same direction as in the continental comparisons and four in the opposite direction. One TE, FBti0019144, differentiated in EC in the same direction as in Australia and in the opposite direction to that of North America.

Conclusions

We presume that the major divergent evolutionary driving force at the local EC microsite is natural selection overruling gene flow. Therefore, after we rejected drift as an explanation of all the large slope differences, we regarded them as putatively adaptive. In order to substantiate the individual TE adaptation, we need to increase the sample sizes and reveal the significant adaptive TEs.

The comparison of local and global studies show only partial similarity in the adaptation of the TEs, because of the dryness of the ecologically tropical climate in EC, in contrast to the wet tropical climate in the global compared climates. Moreover, adaptation of a TE may be expressed only in part of the time and specific localities.

Reviewers

Reviewed by Eugene Koonin, Limsoon Wong and Fyodor Kondrashov. For the full reviews, please go to the Reviewers’ comments section.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0074-5) contains supplementary material, which is available to authorized users.

Microclimatic Divergence in a Mediterranean Canyon Affects Richness, Composition, and Body Size in Saproxylic Beetle Assemblages

Large valleys with opposing slopes may act as a model system with which the effects of strong climatic gradients on biodiversity can be evaluated. The advantage of such comparisons is that the impact of a change of climate can be studied on the same species pool without the need to consider regional differences. The aim of this study was to compare the assemblage of saproxylic beetles on such opposing slopes at Lower Nahal Oren, Mt. Carmel, Israel (also known as "Evolution Canyon") with a 200-800% higher solar radiation on the south-facing (SFS) compared to the north-facing slope (NFS). We tested specific hypotheses of species richness patterns, assemblage structure, and body size resulting from interslope differences in microclimate. Fifteen flight-interception traps per slope were distributed over three elevation levels ranging from 50 to 100 m a.s.l. Richness of saproxylic beetles was on average 34% higher on the SFS compared with the NFS, with no detected influence of elevation levels. Both assemblage structure and average body size were determined by slope aspect, with more small-bodied beetles found on the SFS. Both the increase in species richness and the higher prevalence of small species on the SFS reflect ecological rules present on larger spatial grain (species-energy hypothesis and community body size shift hypothesis), and both can be explained by the metabolic theory of ecology. This is encouraging for the complementary use of micro- and macroclimatic gradients to study impacts of climate warming on biodiversity.

Natural selection causes adaptive genetic resistance in wild emmer wheat against powdery mildew at "Evolution Canyon" microsite, Mt. Carmel, Israel

Background

“Evolution Canyon” (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unraveling evolution in action highlighting the basic evolutionary processes of adaptation and speciation. A major model organism in ECI is wild emmer, Triticum dicoccoides, the progenitor of cultivated wheat, which displays dramatic interslope adaptive and speciational divergence on the tropical-xeric “African” slope (AS) and the temperate-mesic “European” slope (ES), separated on average by 250 m.

Methods

We examined 278 single sequence repeats (SSRs) and the phenotype diversity of the resistance to powdery mildew between the opposite slopes. Furthermore, 18 phenotypes on the AS and 20 phenotypes on the ES, were inoculated by both Bgt E09 and a mixture of powdery mildew races.

Results

In the experiment of genetic diversity, very little polymorphism was identified intra-slope in the accessions from both the AS or ES. By contrast, 148 pairs of SSR primers (53.23%) amplified polymorphic products between the phenotypes of AS and ES. There are some differences between the two wild emmer wheat genomes and the inter-slope SSR polymorphic products between genome A and B. Interestingly, all wild emmer types growing on the south-facing slope (SFS=AS) were susceptible to a composite of Blumeria graminis, while the ones growing on the north-facing slope (NFS=ES) were highly resistant to Blumeria graminis at both seedling and adult stages.

Conclusion/Significance

Remarkable inter-slope evolutionary divergent processes occur in wild emmer wheat, T. dicoccoides at EC I, despite the shot average distance of 250 meters. The AS, a dry and hot slope, did not develop resistance to powdery mildew, whereas the ES, a cool and humid slope, did develop resistance since the disease stress was strong there. This is a remarkable demonstration in host-pathogen interaction on how resistance develops when stress causes an adaptive result at a micro-scale distance.

High environmental stress yields greater tocotrienol content while changing vitamin e profiles of wild emmer wheat seeds

Vitamin E is an essential human nutrient that was first isolated from wheat. Emmer wheat, the cereal of Old World agriculture and a precursor to durum wheat, grows wild in the Fertile Crescent. Evolution Canyon, Israel, provides a microsite that models effects of contrasting environments. The north-facing and south-facing slopes exhibit low and high stress environments, respectively. Wild emmer wheat seeds were collected from both slopes and seed tocochromanol contents measured to test the hypothesis that high stress alters emmer wheat seed tocol-omics. Seeds from high stress areas contained more total vitamin E (108±15 nmol/g) than seeds from low stress environments (80±17 nmol/g, P=.0004). Vitamin E profiles within samples from these different environments revealed significant differences in isoform concentrations. Within each region, β- plus γ-tocotrienols represented the highest concentration of wheat tocotrienols (high stress, P<.0001; low stress, P<.0001), while α-tocopherol represented the highest concentration of the tocopherols (high stress, P=.0002; low stress, P<.0001). Percentages of both δ-tocotrienol and δ-tocopherol increased in high stress conditions. Changes under higher stress apparently are due to increased pathway flux toward more tocotrienol production. The production of more δ-isoforms suggests increased flow through a divergent path controlled by the VTE1 gene. Hence, stress conditions alter plant responses such that vitamin E profiles are changed, likely an attempt to provide additional antioxidant activity to promote seed viability and longevity.

Exploring natural selection to guide breeding for agriculture

Climate change threatens reduced crop production and poses major challenges to food security. The breeding of climate-resilient crop varieties is increasingly urgent. Wild plant populations evolve to cope with changes in their environment due to the forces of natural selection. This adaptation may be followed over time in populations at the same site or explored by examining differences between populations growing in different environments or across an environmental gradient. Survival in the wild has important differences to the objective of agriculture to maximize crop yields. However, understanding the nature of adaptation in wild populations at the whole genome level may suggest strategies for crop breeding to deliver agricultural production with more resilience to climate variability.

Divergence of Drosophila melanogaster repeatomes in response to a sharp microclimate contrast in Evolution Canyon, Israel

Repeat sequences, especially mobile elements, make up large portions of most eukaryotic genomes and provide enormous, albeit commonly underappreciated, evolutionary potential. We analyzed repeatomes of Drosophila melanogaster that have been diverging in response to a microclimate contrast in Evolution Canyon (Mount Carmel, Israel), a natural evolutionary laboratory with two abutting slopes at an average distance of only 200 m, which pose a constant ecological challenge to their local biotas. Flies inhabiting the colder and more humid north-facing slope carried about 6% more transposable elements than those from the hot and dry south-facing slope, in parallel to a suite of other genetic and phenotypic differences between the two populations. Nearly 50% of all mobile element insertions were slope unique, with many of them disrupting coding sequences of genes critical for cognition, olfaction, and thermotolerance, consistent with the observed patterns of thermotolerance differences and assortative mating.

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.

Increases in both acute and chronic temperature potentiate tocotrienol concentrations in wild barley at 'Evolution Canyon'

Biosynthesis of tocols (vitamin E isoforms) is linked to response to temperature in plants. 'Evolution Canyon', an ecogeographical microcosm extending over an average of 200 meters (range 100-400) wide area in the Carmel Mountains of northern Israel, has been suggested as a model for studying global warming. Both domestic (Hordeum vulgare) and wild (Hordeum spontaneum) barley compared with wheat, oat, corn, rice, and rye show high tocotrienol/tocopherol ratios. Therefore, we hypothesized that tocol distribution might change in response to global warming. α-, β-, γ-, and δ-tocopherol, and α-, β-, γ-, and δ-tocotrienol concentrations were measured in wild barley (H. spontaneum) seeds harvested from the xeric (African) and mesic (European) slopes of Evolution Canyon over a six-year period from 2005-2011. Additionally, we examined seeds from areas contiguous to and distant from the part of the Canyon severely burned during the Carmel Fire of December 2010. Increased α-tocotrienol (p<0.01) was correlated with 1) temperature increases, 2) to the hotter 'African' slope in contrast to the cooler 'European' slope, and 3) to propinquity to the fire. The study illustrates the role of α-tocotrienol in both chronic and acute temperature adaptation in wild barley and suggests future research into thermoregulatory mechanisms in plants.

Sympatric incipient speciation of spiny mice Acomys at "Evolution Canyon," Israel

Does the paucity of empirical evidence of sympatric speciation in nature reflect reality, despite theoretical support? Or is it due to inappropriate searches in nature with overly restrictive assumptions and an incorrect null hypothesis? Spiny mice, Acomys, described here at Evolution Canyon (EC) incipiently and sympatrically speciate owing to microclimatic interslope divergence. The opposite slopes at EC vary dramatically, physically and biotically, representing the dry and hot south-facing slope savannoid-African continent ["African" slope (AS)], abutting with the north-facing slope forested south-European continent ["European" slope (ES)]. African-originated spiny mice, of the Acomys cahirinus complex, colonized Israel 30,000 y ago based on fossils. Genotypically, we showed significantly higher genetic diversity of mtDNA and amplified fragment length polymorphism of Acomys on the AS compared with the ES. This is also true regionally across Israel. In complete mtDNA, 25% of the haplotypes at EC were slope-biased. Phenotypically, the opposite slope's populations also showed adaptive morphology, physiology, and behavior divergence paralleling regional populations across Israel. Preliminary tests indicate slope-specific mate choices. Colonization of Acomys at the EC first occurred on the AS and then moved to the ES. Strong slope-specific natural selection (both positive and negative) overrules low interslope gene flow. Both habitat slope selection and mate choices suggest ongoing incipient sympatric speciation. We conclude that Acomys at the EC is ecologically and genetically adaptively, incipiently, sympatrically speciating on the ES owing to adaptive microclimatic natural selection.

Genome differentiation of Drosophila melanogaster from a microclimate contrast in Evolution Canyon, Israel

The opposite slopes of "Evolution Canyon" in Israel have served as a natural model system of adaptation to a microclimate contrast. Long-term studies of Drosophila melanogaster populations inhabiting the canyon have exhibited significant interslope divergence in thermal and drought stress resistance, candidate genes, mobile elements, habitat choice, mating discrimination, and wing-shape variation, all despite close physical proximity of the contrasting habitats, as well as substantial interslope migration. To examine patterns of genetic differentiation at the genome-wide level, we used high coverage sequencing of the flies' genomes. A total of 572 genes were significantly different in allele frequency between the slopes, 106 out of which were associated with 74 significantly overrepresented gene ontology (GO) terms, particularly so with response to stimulus and developmental and reproductive processes, thus corroborating previous observations of interslope divergence in stress response, life history, and mating functions. There were at least 37 chromosomal "islands" of interslope divergence and low sequence polymorphism, plausible signatures of selective sweeps, more abundant in flies derived from one (north-facing) of the slopes. Positive correlation between local recombination rate and the level of nucleotide polymorphism was also found.

SNP-revealed genetic diversity in wild emmer wheat correlates with ecological factors

Background

Patterns of genetic diversity between and within natural plant populations and their driving forces are of great interest in evolutionary biology. However, few studies have been performed on the genetic structure and population divergence in wild emmer wheat using a large number of EST-related single nucleotide polymorphism (SNP) markers.

Results

In the present study, twenty-five natural wild emmer wheat populations representing a wide range of ecological conditions in Israel and Turkey were used. Genetic diversity and genetic structure were investigated using over 1,000 SNP markers. A moderate level of genetic diversity was detected due to the biallelic property of SNP markers. Clustering based on Bayesian model showed that grouping pattern is related to the geographical distribution of the wild emmer wheat. However, genetic differentiation between populations was not necessarily dependent on the geographical distances. A total of 33 outlier loci under positive selection were identified using a F_ST-outlier method. Significant correlations between loci and ecogeographical factors were observed.

Conclusions

Natural selection appears to play a major role in generating adaptive structures in wild emmer wheat. SNP markers are appropriate for detecting selectively-channeled adaptive genetic diversity in natural populations of wild emmer wheat. This adaptive genetic diversity is significantly associated with ecological factors.