Transposon-triggered innate immune response confers cancer resistance to the blind mole rat.

Blind mole rats (BMRs) are small rodents, characterized by an exceptionally long lifespan (>21 years) and resistance to both spontaneous and induced tumorigenesis. Here we report that cancer resistance in the BMR is mediated by retrotransposable elements (RTEs). Cells and tissues of BMRs express very low levels of DNA methyltransferase 1. Following cell hyperplasia, the BMR genome DNA loses methylation, resulting in the activation of RTEs. Upregulated RTEs form cytoplasmic RNA-DNA hybrids, which activate the cGAS-STING pathway to induce cell death. Although this mechanism is enhanced in the BMR, we show that it functions in mice and humans. We propose that RTEs were co-opted to serve as tumor suppressors that monitor cell proliferation and are activated in premalignant cells to trigger cell death via activation of the innate immune response. Activation of RTEs is a double-edged sword, serving as a tumor suppressor but contributing to aging in late life via the induction of sterile inflammation.

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei.

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

Sympatric speciation of the spiny mouse from Evolution Canyon in Israel substantiated genomically and methylomically.

SignificanceWhether sympatric speciation (SS) is rare or common is still debated. Two populations of the spiny mouse, Acomys cahirinus, from Evolution Canyon I (EC I) in Israel have been depicted earlier as speciating sympatrically by molecular markers and transcriptome. Here, we investigated SS both genomically and methylomically, demonstrating that the opposite populations of spiny mice are sister taxa and split from the common ancestor around 20,000 years ago without an allopatric history. Mate choice, olfactory receptors, and speciation genes contributed to prezygotic/postzygotic reproductive isolation. The two populations showed different methylation patterns, facilitating adaptation to their local environment. They cope with abiotic and biotic stresses, due to high solar interslope radiation differences. We conclude that our new genomic and methylomic data substantiated SS.

Genomic insights into zokors' phylogeny and speciation in China.

The phylogeny and speciation of subterranean zokors in China are unclear, as previous studies on morphology and limited molecular markers have generated conflicting results. This study unraveled the complex evolutionary history of eight zokor species in China based on de novo assembly at chromosome level and whole-genome sequencing of 23 populations. We found extensive phylogenetic discordances between nuclear and mitochondrial phylogenies, and different coalescent phylogenies, which could be explained by introgression and incomplete lineage sorting (ILS). The recent Qinghai-Tibet Plateau uplift (∼3.60 million y ago; Mya) drove Eospalax to speciate into clade A and clade B (∼3.22 Mya), and discordant phylogenies in this node were mainly attributed to introgression rather than ILS. Clade A rapidly diverged into three lineages due to geographical isolation and glaciation, while glaciation and C4 plant expansion contributed to the speciation of clade B. ILS contributed to the discordances of two rapidly radiated nodes rather than introgression. The effective population sizes (Ne's) of all the species of Eospalax were affected by three glaciations. Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of all the species pairs. Positively selected genes putatively related to specific inhabitation adaptations were identified, such as heart development, neurogenesis, DNA repair, and immune response. Climate, geological tectonism, and C4 vegetation shaped the adaptation and speciation of zokors in China.

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.

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation.

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

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.

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 evolution of blind subterranean mole rats: Adaptive peripatric versus sympatric speciation.

Speciation mechanisms remain controversial. Two speciation models occur in Israeli subterranean mole rats, genus Spalax: a regional speciation cline southward of four peripatric climatic chromosomal species and a local, geologic-edaphic, genic, and sympatric speciation. Here we highlight their genome evolution. The five species were separated into five genetic clusters by single nucleotide polymorphisms, copy number variations (CNVs), repeatome, and methylome in sympatry. The regional interspecific divergence correspond to Pleistocene climatic cycles. Climate warmings caused chromosomal speciation. Triple effective population size, Ne , declines match glacial cold cycles. Adaptive genes evolved under positive selection to underground stresses and to divergent climates, involving interspecies reproductive isolation. Genomic islands evolved mainly due to adaptive evolution involving ancient polymorphisms. Repeatome, including both CNV and LINE1 repetitive elements, separated the five species. Methylation in sympatry identified geologically chalk-basalt species that differentially affect thermoregulation, hypoxia, DNA repair, P53, and other pathways. Genome adaptive evolution highlights climatic and geologic-edaphic stress evolution and the two speciation models, peripatric and sympatric.

Sympatric speciation of wild emmer wheat driven by ecology and chromosomal rearrangements.

In plants, the mechanism for ecological sympatric speciation (SS) is little known. Here, after ruling out the possibility of secondary contact, we show that wild emmer wheat, at the microclimatically divergent microsite of "Evolution Canyon" (EC), Mt. Carmel, Israel, underwent triple SS. Initially, it split following a bottleneck of an ancestral population, and further diversified to three isolated populations driven by disruptive ecological selection. Remarkably, two postzygotically isolated populations (SFS1 and SFS2) sympatrically branched within an area less than 30 m at the tropical hot and dry savannoid south-facing slope (SFS). A series of homozygous chromosomal rearrangements in the SFS1 population caused hybrid sterility with the SFS2 population. We demonstrate that these two populations developed divergent adaptive mechanisms against severe abiotic stresses on the tropical SFS. The SFS2 population evolved very early flowering, while the SFS1 population alternatively evolved a direct tolerance to irradiance by improved ROS scavenging activity that potentially accounts for its evolutionary fate with unstable chromosome status. Moreover, a third prezygotically isolated sympatric population adapted on the abutting temperate, humid, cool, and forested north-facing slope (NFS), separated by 250 m from the SFS wild emmer wheat populations. The NFS population evolved multiple resistant loci to fungal diseases, including powdery mildew and stripe rust. Our study illustrates how plants sympatrically adapt and speciate under disruptive ecological selection of abiotic and biotic stresses.

Thermal biology in the Upper Galili Mountain blind mole rat (Nannospalax galili) and an overview of spalacine energetics.

Several hundred mammalian species thrive in complex burrow systems, which protect them from climatic extremes and predation. At the same time, it is also a stressful environment due to low food supply, high humidity, and, in some cases, a hypoxic and hypercapnic atmosphere. To face such conditions, subterranean rodents have convergently evolved low basal metabolic rate, high minimal thermal conductance and low body temperature. Although these parameters have been intensively studied in the last decades, such information is far from being well-known in one of the most studied groups of subterranean rodents, the blind mole rats of the genus Nannospalax. The lack of information is particularly noticeable for parameters such as the upper critical temperature and the width of the thermoneutral zone. In our study, we analysed the energetics of the Upper Galilee Mountain blind mole rat Nannospalax galili and found its basal metabolic rate of 0.84 ± 0.10 mL O2×g-1 × h-1, thermoneutral zone between 28 and 35 °C, mean Tb within the zone of 36.3 ± 0.6 °C, and minimal thermal conductance equal to 0.082 mL O2×g-1 × h-1 × C-1. Nannospalax galili is a truly homeothermic rodent well adapted to face lower ambient temperatures, because its Tb was stable down to the lowest temperature measured (10 °C). At the same time, a relatively high basal metabolic rate and relatively low minimal thermal conductance for a subterranean rodent of such body mass, and the difficulty of surviving ambient temperatures slightly above upper critical temperature, indicates problems with sufficient heat dissipation at higher temperatures. This can easily lead to overheating, that is relevant mainly during the hot-dry season. These findings suggest that N. galili can be threatened by ongoing global climate change.

The Influence of Edaphic Factors on DNA Damage and Repair in Wild Wheat Triticum dicoccoides Körn. (Poaceae, Triticeae).

A complex DNA repair network maintains genome integrity and genetic stability. In this study, the influence of edaphic factors on DNA damage and repair in wild wheat Triticum dicoccoides was addressed. Plants inhabiting two abutting microsites with dry terra rossa and humid basalt soils were studied. The relative expression level of seven genes involved in DNA repair pathways-RAD51, BRCA1, LigIV, KU70, MLH1, MSH2, and MRE11-was assessed using quantitative real-time PCR (qPCR). Immunolocalization of RAD51, LigIV, γH2AX, RNA Polymerase II, and DNA-RNA hybrid [S9.6] (R-loops) in somatic interphase nuclei and metaphase chromosomes was carried out in parallel. The results showed a lower expression level of genes involved in DNA repair and a higher number of DNA double-strand breaks (DSBs) in interphase nuclei in plants growing in terra rossa soil compared with plants in basalt soil. Further, the number of DSBs and R-loops in metaphase chromosomes was also greater in plants growing on terra rossa soil. Finally, RAD51 and LigIV foci on chromosomes indicate ongoing DSB repair during the M-phase via the Homologous Recombination and Non-Homologous End Joining pathways. Together, these results show the impact of edaphic factors on DNA damage and repair in the wheat genome adapted to contrasting environments.

Local Adaptation of Bitter Taste and Ecological Speciation in a Wild Mammal.

Sensory systems are attractive evolutionary models to address how organisms adapt to local environments that can cause ecological speciation. However, tests of these evolutionary models have focused on visual, auditory, and olfactory senses. Here, we show local adaptation of bitter taste receptor genes in two neighboring populations of a wild mammal-the blind mole rat Spalax galili-that show ecological speciation in divergent soil environments. We found that basalt-type bitter receptors showed higher response intensity and sensitivity compared with chalk-type ones using both genetic and cell-based functional analyses. Such functional changes could help animals adapted to basalt soil select plants with less bitterness from diverse local foods, whereas a weaker reception to bitter taste may allow consumption of a greater range of plants for animals inhabiting chalk soil with a scarcity of food supply. Our study shows divergent selection on food resources through local adaptation of bitter receptors, and suggests that taste plays an important yet underappreciated role in speciation.

The CC-NB-LRR protein BSR1 from Brachypodium confers resistance to Barley stripe mosaic virus in gramineous plants by recognising TGB1 movement protein.

Although some nucleotide binding, leucine-rich repeat immune receptor (NLR) proteins conferring resistance to specific viruses have been identified in dicot plants, NLR proteins involved in viral resistance have not been described in monocots. We have used map-based cloning to isolate the CC-NB-LRR (CNL) Barley stripe mosaic virus (BSMV) resistance gene barley stripe resistance 1 (BSR1) from Brachypodium distachyon Bd3-1 inbred line. Stable BSR1 transgenic Brachypodium line Bd21-3, barley (Golden Promise) and wheat (Kenong 199) plants developed resistance against BSMV ND18 strain. Allelic variation analyses indicated that BSR1 is present in several Brachypodium accessions collected from countries in the Middle East. Protein domain swaps revealed that the intact LRR domain and the C-terminus of BSR1 are required for resistance. BSR1 interacts with the BSMV ND18 TGB1 protein in planta and shows temperature-sensitive antiviral resistance. The R390 and T392 residues of TGB1ND (ND18 strain) and the G196 and K197 residues within the BSR1 P-loop motif are key amino acids required for immune activation. BSR1 is the first cloned virus resistance gene encoding a typical CNL protein in monocots, highlighting the utility of the Brachypodium model for isolation and analysis of agronomically important genes for crop improvement.

Elevated mutation and selection in wild emmer wheat in response to 28 years of global warming.

Global warming has been documented to threaten wild plants with strong selection pressures, but how plant populations respond genetically to the threats remains poorly understood. We characterized the genetic responses of 10 wild emmer wheat (Triticum dicoccoides Koern.; WEW) populations in Israel, sampling them in 1980 and again in 2008, through an exome capture analysis. It was found that these WEW populations were under elevated selection, displayed reduced diversity and temporal divergence, and carried increased mutational burdens forward. However, some populations still showed the ability to acquire beneficial alleles via selection or de novo mutation for future adaptation. Grouping populations with mean annual rainfall and temperature revealed significant differences in most of the 14 genetic estimates in either sampling year or over the 28 y. The patterns of genetic response to rainfall and temperature varied and were complex. In general, temperature groups displayed more temporal differences in genetic response than rainfall groups. The highest temperature group had more deleterious single nucleotide polymorphisms (dSNPs), higher nucleotide diversity, fewer selective sweeps, lower differentiation, and lower mutational burden. The least rainfall group had more dSNPs, higher nucleotide diversity, lower differentiation and higher mutational burden. These characterized genetic responses are significant, allowing not only for better understanding of evolutionary changes in the threatened populations, but also for realistic modeling of plant population adaptability and vulnerability to global warming.

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.

Evolution of chloroplast retrograde signaling facilitates green plant adaptation to land.

Chloroplast retrograde signaling networks are vital for chloroplast biogenesis, operation, and signaling, including excess light and drought stress signaling. To date, retrograde signaling has been considered in the context of land plant adaptation, but not regarding the origin and evolution of signaling cascades linking chloroplast function to stomatal regulation. We show that key elements of the chloroplast retrograde signaling process, the nucleotide phosphatase (SAL1) and 3'-phosphoadenosine-5'-phosphate (PAP) metabolism, evolved in streptophyte algae-the algal ancestors of land plants. We discover an early evolution of SAL1-PAP chloroplast retrograde signaling in stomatal regulation based on conserved gene and protein structure, function, and enzyme activity and transit peptides of SAL1s in species including flowering plants, the fern Ceratopteris richardii, and the moss Physcomitrella patens Moreover, we demonstrate that PAP regulates stomatal closure via secondary messengers and ion transport in guard cells of these diverse lineages. The origin of stomata facilitated gas exchange in the earliest land plants. Our findings suggest that the conquest of land by plants was enabled by rapid response to drought stress through the deployment of an ancestral SAL1-PAP signaling pathway, intersecting with the core abscisic acid signaling in stomatal guard cells.

Evolution of rapid blue-light response linked to explosive diversification of ferns in angiosperm forests.

Ferns appear in the fossil record some 200 Myr before angiosperms. However, as angiosperm-dominated forest canopies emerged in the Cretaceous period there was an explosive diversification of modern (leptosporangiate) ferns, which thrived in low, blue-enhanced light beneath angiosperm canopies. A mechanistic explanation for this transformative event in the diversification of ferns has remained elusive. We used physiological assays, transcriptome analysis and evolutionary bioinformatics to investigate a potential connection between the evolution of enhanced stomatal sensitivity to blue light in modern ferns and the rise of angiosperm-dominated forests in the geological record. We demonstrate that members of the largest subclade of leptosporangiate ferns, Polypodiales, have significantly faster stomatal response to blue light than more ancient fern lineages and a representative angiosperm. We link this higher sensitivity to levels of differentially expressed genes in blue-light signaling, particularly in the cryptochrome (CRY) signaling pathway. Moreover, CRYs of the Polypodiales examined show gene duplication events between 212.9-196.9 and 164.4-151.8 Ma, when angiosperms were emerging, which are lacking in other major clades of extant land plants. These findings suggest that evolution of stomatal blue-light sensitivity helped modern ferns exploit the shady habitat beneath angiosperm forest canopies, fueling their Cretaceous hyperdiversification.

Phenotypic and genomic adaptations to the extremely high elevation in plateau zokor (Myospalax baileyi).

The evolutionary outcomes of high elevation adaptation have been extensively described. However, whether widely distributed high elevation endemic animals adopt uniform mechanisms during adaptation to different elevational environments remains unknown, especially with respect to extreme high elevation environments. To explore this, we analysed the phenotypic and genomic data of seven populations of plateau zokor (Myospalax baileyi) along elevations ranging from 2,700 to 4,300 m. Based on whole-genome sequencing data and demographic reconstruction of the evolutionary history, we show that two populations of plateau zokor living at elevations exceeding 3,700 m diverged from other populations nearly 10,000 years ago. Further, phenotypic comparisons reveal stress-dependent adaptation, as two populations living at elevations exceeding 3,700 m have elevated ratios of heart mass to body mass relative to other populations, and the highest population (4,300 m) displays alterations in erythrocytes. Correspondingly, genomic analysis of selective sweeps indicates that positive selection might contribute to the observed phenotypic alterations in these two extremely high elevation populations, with the adaptive cardiovascular phenotypes of both populations possibly evolving under the functional constrains of their common ancestral population. Taken together, phenotypic and genomic evidence demonstrates that heterogeneous stressors impact adaptations to extreme elevations and reveals stress-dependent and genetically constrained adaptation to hypoxia, collectively providing new insights into the high elevation adaptation.