Long-term genetic studies of wild populations are very scarce, but are essential for connecting ecological and population genetics models, and for understanding the dynamics of biodiversity. We present a study of a wild wheat population sampled over a 36-year period at high spatial resolution. We genotyped 832 individuals from regular sampling along transects during the course of the experiment. Genotypes were clustered into ecological microhabitats over scales of tens of metres, and this clustering was remarkably stable over the 36 generations of the study. Simulations show that it is difficult to determine whether this spatial and temporal stability reflects extremely limited dispersal or fine-scale local adaptation to ecological parameters. Using a common-garden experiment, we showed that the genotypes found in distinct microhabitats differ phenotypically. Our results provide a rare insight into the population genetics of a natural population over a long monitoring period. 

Plant roots integrate environmental signals with development using exquisite spatiotemporal control. This is apparent in the deposition of suberin, an apoplastic diffusion barrier, which regulates flow of water, solutes and gases, and is environmentally plastic. Suberin is considered a hallmark of endodermal differentiation but is absent in the tomato endodermis. Instead, suberin is present in the exodermis, a cell type that is absent in the model organism Arabidopsis thaliana. Here we demonstrate that the suberin regulatory network has the same parts driving suberin production in the tomato exodermis and the Arabidopsis endodermis. Despite this co-option of network components, the network has undergone rewiring to drive distinct spatial expression and with distinct contributions of specific genes. Functional genetic analyses of the tomato MYB92 transcription factor and ASFT enzyme demonstrate the importance of exodermal suberin for a plant water-deficit response and that the exodermal barrier serves an equivalent function to that of the endodermis and can act in its place. 

Genetic dissection of a QTL from wild emmer wheat, QGpc.huj.uh-5B.2, introgressed into bread wheat, identified candidate genes associated with tolerance to nitrogen deficiency, and potentially useful for improving nitrogen-use efficiency. Nitrogen (N) is an important macronutrient critical to wheat growth and development; its deficiency is one of the main factors causing reductions in grain yield and quality. N availability is significantly affected by drought or flooding, that are dependent on additional factors including soil type or duration and severity of stress. In a previous study, we identified a high grain protein content QTL (QGpc.huj.uh-5B.2) derived from the 5B chromosome of wild emmer wheat, that showed a higher proportion of explained variation under water-stress conditions. We hypothesized that this QTL is associated with tolerance to N deficiency as a possible mechanism underlying the higher effect under stress. To validate this hypothesis, we introgressed the QTL into the elite bread wheat var. Ruta, and showed that under N-deficient field conditions the introgression IL99 had a 33% increase in GPC (p < 0.05) compared to the recipient parent. Furthermore, evaluation of IL99 response to severe N deficiency (10% N) for 14 days, applied using a semi-hydroponic system under controlled conditions, confirmed its tolerance to N deficiency. Fine-mapping of the QTL resulted in 26 homozygous near-isogenic lines (BC4F5) segregating to N-deficiency tolerance. The QTL was delimited from - 28.28 to - 1.29 Mb and included 13 candidate genes, most associated with N-stress response, N transport, and abiotic stress responses. These genes may improve N-use efficiency under severely N-deficient environments. Our study demonstrates the importance of WEW as a source of novel candidate genes for sustainable improvement in tolerance to N deficiency in wheat.

Hypothesizing that CpG codon dyads, formed by consecutive codons containing a cytosine-guanine pair (NNC-GNN), may play a crucial role in gene function, we conducted an extensive analysis to investigate their distribution and conservation within mammalian genes. Our findings reveal that genes characterized by a high density of CpG codon dyads are notably associated with homeobox domains and RNA polymerase II transcription factors. Conversely, genes exhibiting low CpG codon dyad density have links to DNA damage repair and mitosis. Importantly, our study identifies a re-markable increase in expressed genes that harbor CpG during embryonic development, suggesting their potential involvement in gene regulation at these developmental stages. These results under-score the functional significance of CpG codon dyads in DNA methylation and gene expression, fur-ther demonstrating the coevolution of consecutive codons and their contribution to codon usage bias. 

Zokors, an Asiatic group of subterranean rodents, originated in lowlands and colonized high-elevational zones following the uplift of the Qinghai–Tibet plateau about 3.6 million years ago. Zokors live at high elevation in subterranean burrows and experience hypobaric hypoxia, including both hypoxia (low oxygen concentration) and hypercapnia (elevated partial pressure of CO2). Here we report a genomic analysis of six zokor species (genus Eospalax) with different elevational ranges to identify structural variants (deletions and inversions) that may have contributed to high-elevation adaptation. Based on an assembly of a chromosome-level genome of the high-elevation species, Eospalax baileyi, we identified 18 large inversions that distinguished this species from congeners native to lower elevations. Small-scale structural variants in the introns of EGLN1, HIF1A, HSF1 and SFTPD of E. baileyi were associated with the upregulated expression of those genes. A rearrangement on chromosome 1 was associated with altered chromatin accessibility, leading to modified gene expression profiles of key genes involved in the physiological response to hypoxia. Multigene families that underwent copy-number expansions in E. baileyi were enriched for autophagy, HIF1 signalling and immune response. E. baileyi show a significantly larger lung mass than those of other Eospalax species. These findings highlight the key role of structural variants underlying hypoxia adaptation of high-elevation species in Eospalax. 

The diversity of Lycoperdaceae in Israel was studied. Molecular phylogenetic relationships within the family, and genus Lycoperdon in particular, were inferred using original ITS rDNA sequences of 58 samples belonging to 25 species from Israel and six other countries, together with 66 sequences stored in the GenBank database. The current molecular phylogenetic study recovered the family Lycoperdaceae as a monophyletic group, which was supported in both maximum likelihood and Bayesian analyses. The clades corresponding to the genera Apioperdon, Bovista, Calvatia, Disciseda, and Lycoperdon were revealed. The taxonomic structure of the named genera was partially resolved. Within the genus Lycoperdon, some species received significant statistical support; however, their relationships, as well as the problem of the genus monophyly, mostly remained questionable. As a result of a thorough literature survey, extensive sample collection, and studies of the material stored in the Herbarium of the Institute of Evolution, University of Haifa (HAI, Israel), fifteen species representing five genera were found in the territory of Israel. Six species, namely Apioperdon pyriforme, Bovista aestivalis, Calvatia candida, Lycoperdon decipiens, L. niveum, and L. perlatum, are new additions to the diversity of Lycoperdaceae in Israel. Detailed macro- and micromorphological descriptions, ecology, geography, and critical notes, together with light microscopy photos and SEM micrographs, are provided. In-depth discussion on some taxonomically challenging species is presented. 

Two studies describe kinase fusion proteins (KFPs) that regulate the perception and deception of wheat pathogens. These highlight the emergence of KFPs as plant immune regulators and emphasize the importance of crop wild relatives as a reservoir for resistance breeding and global food security. 

Wheat, an essential crop for global food security, is well adapted to a wide variety of soils. However, the gene networks shaping different root architectures remain poorly understood. We report here that dosage differences in a cluster of monocot-specific 12-OXOPHYTODIENOATE REDUCTASE genes from subfamily III (OPRIII) modulate key differences in wheat root architecture, which are associated with grain yield under water-limited conditions. Wheat plants with loss-of-function mutations in OPRIII show longer seminal roots, whereas increased OPRIII dosage or transgenic over-expression result in reduced seminal root growth, precocious development of lateral roots and increased jasmonic acid (JA and JA-Ile). Pharmacological inhibition of JA-biosynthesis abolishes root length differences, consistent with a JA-mediated mechanism. Transcriptome analyses of transgenic and wild-type lines show significant enriched JA-biosynthetic and reactive oxygen species (ROS) pathways, which parallel changes in ROS distribution. OPRIII genes provide a useful entry point to engineer root architecture in wheat and other cereals 

The aim of the present study was to identify the structure of active compounds in Cyathus stratus that previously demonstrated anti-pancreatic cancer activity. The active compounds were purified from a crude extract by a series of RP-18 preparative chromatography using homemade octadecyl silica gel column. HPLC injection of the crude extract revealed a chromatogram with three main peaks with retention times (RT) 15.6, 18.2, and 22.5 min. Each fraction that exhibited promising activity in vitro was further separated using various available chromatographic techniques. The purified compound with the ultimate anti-cancer activity appeared at RT of 15.8 in the HPLC chromatogram with more than 90% purity. The main peak at the mass spectra appeared at m/z = 446.2304 with the calculated molecular formula of C25H34O7. One- and two-dimensional NMR analyses indicated that the structure of the active molecule (peak 15.8 min in HPLC) was identified as striatal C. Exposure of human pancreatic cancer cells to purified striatal C resulted in induction of apoptosis. Further studies are needed in order to develop a method for the synthesis of striatal in order to use it in clinical studies for treatment of cancer. 

Species that hibernate generally live longer than would be expected based solely on their body size. Hibernation is characterized by long periods of metabolic suppression (torpor) interspersed by short periods of increased metabolism (arousal). The torpor–arousal cycles occur multiple times during hibernation, and it has been suggested that processes controlling the transition between torpor and arousal states cause ageing suppression. Metabolic rate is also a known correlate of longevity; we thus proposed the ‘hibernation–ageing hypothesis’ whereby ageing is suspended during hibernation. We tested this hypothesis in a well-studied population of yellow-bellied marmots (Marmota flaviventer), which spend 7–8 months per year hibernating. We used two approaches to estimate epigenetic age: the epigenetic clock and the epigenetic pacemaker. Variation in epigenetic age of 149 samples collected throughout the life of 73 females was modelled using generalized additive mixed models (GAMM), where season (cyclic cubic spline) and chronological age (cubic spline) were fixed effects. As expected, the GAMM using epigenetic ages calculated from the epigenetic pacemaker was better able to detect nonlinear patterns in epigenetic ageing over time. We observed a logarithmic curve of epigenetic age with time, where the epigenetic age increased at a higher rate until females reached sexual maturity (two years old). With respect to circannual patterns, the epigenetic age increased during the active season and essentially stalled during the hibernation period. Taken together, our results are consistent with the hibernation–ageing hypothesis and may explain the enhanced longevity in hibernators. 

Stripe rust is a devastating disease in wheat that causes substantial yield loss around the world. The most effective strategy for mitigating yield loss is to develop resistant cultivars. The wild relatives of wheat are good sources of resistance to fungal pathogens. Here, we used a genome-wide association study (GWAS) to identify loci associated with stripe rust (causal agent: Puccinia striiformis f. sp. tritici) resistance in wild emmer (Triticum dicoccoides) at the seedling stage, in the greenhouse, and at the adult plant stage, in the field. We found that the two major loci contributing to resistance in our wild emmer panel were the previously cloned seedling-stage resistance gene, Yr15, and the adult-plant-stage resistance gene, Yr36. Nevertheless, we detected 12 additional minor QTLs that additionally contribute to adult plant resistance and mapped a locus on chromosome 3AS that tentatively harbors a novel seedling resistance gene. The genotype and phenotype data generated for the wild emmer panel, together with the detected SNPs associated with resistance to stripe rust, provide a valuable resource for disease-resistance breeding in durum and bread wheat 

While it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the 'HbS mutation' when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible 

The boletoid genera Butyriboletus and Exsudoporus have recently been suggested by some researchers to constitute a single genus, and Exsudoporus was merged into Butyriboletus as a later synonym. However, no convincing arguments have yet provided significant evidence for this congeneric placement. In this study, we analyze material from Exsudoporus species and closely related taxa to assess taxonomic and phylogenetic boundaries between these genera and to clarify species delimitation within Exsudoporus. Outcomes from a multilocus phylogenetic analysis (ITS, nrLSU, tef1-α and rpb2) clearly resolve Exsudoporus as a monophyletic, homogenous and independent genus that is sister to Butyriboletus. An accurate morphological description, comprehensive sampling, type studies, line drawings and a historical overview on the nomenclatural issues of the type species E. permagnificus are provided. Furthermore, this species is documented for the first time from Israel in association with Quercus calliprinos. The previously described North American species Exsudoporus frostii and E. floridanus are molecularly confirmed as representatives of Exsudoporus, and E. floridanus is epitypified. The eastern Asian species Leccinum rubrum is assigned here to Exsudoporus based on molecular evidence, and a new combination is proposed. Sequence data from the original material of the Japanese Boletus kermesinus were generated, and its conspecificity with L. rubrum is inferred as formerly presumed based on morphology. Four additional cryptic species from North and Central America previously misdetermined as either B. frostii or B. floridanus are phylogenetically placed but remain undescribed due to the paucity of available material. Boletus weberi (syn. B. pseudofrostii) and Xerocomus cf. mcrobbii cluster outside of Exsudoporus and are herein assigned to the recently described genus Amoenoboletus. Biogeographic distribution patterns are elucidated, and a dichotomous key to all known species of Exsudoporus worldwide is presented. 

This paper presents data compilation for indicator species of organic pollution and trophic state of continental waters. Information was collected from research papers, monographs, electronic resources, and our research. Altogether 316 species of higher plants, plants, mosses, charophytes, protozoan, and bacteria from 11 taxonomical phyla are represented with ecological preferences for saprobity indicators with saprobity index (S) and indicators of trophic state. This comprehensive data can be used for the purpose of aquatic ecosystem assessment and monitoring of water quality based on bioindication methods. 

The river basin of Siversky Donets is of great scientific interest since this river runs through a territory with heavy industry (in particular, coal mining, chemical processing and metal industries). Within the basin, rivers of different sizes were explored (small, medium, large and extra-large) that flow through siliceous and calcareous rocks on the same elevation (lowland—below 200 m a.s.l.). Phytoplankton, as one of the Biological Quality Element, was used to perform the assessment of ecological status of the water bodies within the Siversky Donets river basin in 2019. The state monitoring program based on the updated approaches has been implemented in the river basin for the first time. The composition of phytoplankton species in the basin comprised 167 species (168 intraspecies taxa), mainly Bacillariophyta (63%) and Chlorophyta (22%) with the presence of other species (Cyanobacteria, Charophyta, Chrysophyta, Dinophyta and Euglenophyta). High species diversity and divisions amount are a distinctive property of the smaller rivers, while the bigger rivers show lower number of divisions. The “bloom” events, which are important ecological factors, were not detected in the Siversky Donets river basin. Algal species composition in plankton samples of the basin was identified and series of ecological parameters, such as habitat preferences, temperature, pH, salinity, oxygenation and organic water pollution according to Watanabe and Sládeček’s index of saprobity (S) trophic state and nitrogen uptake metabolism were analyzed. The ecological conclusions were also verified by a canonical correspondence analysis (CCA). The significance of the Canonical Correspondence Analysis (CCA) results was estimated of by a Monte-Carlo permutation test. The high concentrations of inorganic phosphorus compounds (permanganate index (CODMn)) and nitrite ions favored the diversity of Chlorophyta and Cyanobacteria diversity correlated with the levels of bicarbonate and CODMn. High diversity of diatoms was facilitated by the total amount of dissolved solids and chemical oxygen demand (COD). It was found that low water quality could be associated with conditions leading to predominant growth of the mentioned groups of algae. According to the analysis, the highest water quality was characterized by balanced phytoplankton composition and optimal values of the environmental variables. The sites with reference conditions are proposed for future monitoring. 

The destructive wheat powdery mildew disease is caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt). PmG3M, derived from wild emmer wheat Triticum dicoccoides accession G305-3M, is a major gene providing a wide-spectrum resistance against Bgt. PmG3M was previously mapped to wheat chromosome 6B using an F6 recombinant inbred line (RIL) mapping population generated by crossing G305-3M with the susceptible T. durum wheat cultivar Langdon (LDN). In the current study, we aimed to explore the defense mechanisms conferred by PmG3M against Bgt. Histopathology of fungal development was characterized in artificially inoculated leaves of G305-3M, LDN, and homozygous RILs using fluorescence and light microscopy. G305-3M exhibited H2O2 accumulation typical of a hypersensitive response, which resulted in programmed cell death (PCD) in Bgt-penetrated epidermal cells, while LDN showed well-developed colonies without PCD. In addition, we observed a post-haustorial resistance mechanism that arrested the development of fungal feeding structures and pathogen growth in both G305-3M and resistant RIL, while LDN and a susceptible RIL displayed fully developed digitated haustoria and massive accumulation of fungal biomass. In contrast, both G305-3M and LDN exhibited callose deposition in attempt to prevent fungal invasion, supporting this as a mechanism of a basal defense response not associated with PmG3M resistance mechanism per se. The presented results shed light on the resistance mechanisms conferred by PmG3M against wheat powdery mildew. 

Net blotch (NB) is a major disease of barley caused by the fungus Pyrenophora teres f. teres (Ptt), and P. teres f. maculata (Ptm). Ptt and Ptm infect the cultivated crop (Hordeum vulgare) and its wild relatives (H. vulgare ssp. spontaneum and H. murinum ssp. glaucum). The main goal of this research was to study the NB-causing pathogen in the crop center of origin. To address this, we have constructed a Ptt (n = 15) and Ptm (n = 12) collection isolated from three barley species across Israel. Isolates were characterized genetically and phenotypically. Aggressiveness of the isolates was determined based on necrotrophic growth rate on detached leaves of barley. In addition, isolates were genetically characterized by the mating type, followed by phylogenetic analysis, clustering them into seven groups. The analysis showed no significant differentiation of isolates based on either geographic origin, host of origin or form (Ptt vs. Ptm). Nevertheless, there was a significant difference in aggressiveness among the isolates regardless of host species, geographic location or sampling site. Moreover, it was apparent that the isolates derived from wild hosts were more variable in their necrotrophic growth rate, compared to isolates sampled from cultivated hosts, thereby suggesting that NB plays a major role in epidemiology at the center of barley origin where most of the diversity lies. Ptm has significantly higher necrotrophic and saprotrophic growth rates than Ptt, and for both a significant negative correlation was found between light intensity exposure and growth rates 

We studied the effect of microclimatic gradient on the spatial (8 habitats) and seasonal (autumn, winter, and summer) variations of culturable soil microfungal communities at the north-facing slope (NFS) and the south-facing slope (SFS) of Nahal Boker, the Central Negev Desert, Israel. Using the soil dilution plate method, a total of 122 species from 56 genera was isolated. The communities' composition was subjected to pronounced interslope variations in each season. While xerotolerant melanin-containing species (mainly with the protective multicellular spore morphology) predominated in the soil of NFS, peaking up in the summer, thermotolerant and thermophilic Aspergillus spp. were especially abundant at the SFS. Aspergilli also dominated the thermotolerant communities isolated at 37 °C from both slopes. The more severe microclimatic regime at the SFS resulted in the least pronounced spatial and seasonal intraslope variability of the community composition as compared to the more benign and heterogeneous NFS. The quantitative parameter—the density of microfungal isolates, exhibited significant and positive relationship with organic matter and moisture content. In seasonal dynamics, the severe summer conditions significantly reduced the differences in isolate density at the NFS, which were markedly expressed in the autumn and winter between habitats receiving additional water by runoff and the adjacent habitats that did not benefit from runoff water.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail. 

Wild emmer wheat (WEW), the tetraploid progenitor of durum and bread wheat, is a valuable genetic resource for resistance to powdery mildew fungal disease caused by Blumeria graminis f. sp. tritici (Bgt). PmG16 gene, derived from WEW, confers high resistance to most tested Bgt isolates. We mapped PmG16 to a 1.4-cM interval between the flanking markers uhw386 and uhw390 on Chromosome 7AL. Based on gene annotation of WEW reference genome Zavitan_V1, 34 predicted genes were identified within the ~ 3.48-Mb target region. Six genes were annotated as associated with disease resistance, of which TRIDC7AG077150.1 was found to be highly similar to Pm60, previously cloned from Triticum urartu, and resides in the same syntenic region. The functional molecular marker (FMM) for Pm60 (M-Pm60-S1) co-segregated with PmG16, suggesting the Pm60 ortholog from WEW (designated here as TdPm60) as a strong candidate for PmG16. Sequence alignment identified only eight SNPs that differentiate between TdPm60 and TuPm60. Furthermore, TdPm60 was found to be present also in the WEW donor lines of the powdery mildew resistance genes MlIW172 and MlIW72, mapped to the same region of Chromosome 7AL as PmG16, suggesting that TdPm60 constitutes a candidate also for these genes. Furthermore, screening of additional 230 WEW accessions with Pm60 specific markers revealed 58 resistant accessions from the Southern Levant that harbored TdPm60, while none of the susceptible accessions showed the presence of this gene. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta conferred resistance against several local Bgt isolates.

In spite of the obvious climate changes effects on the Carpathian Basin hydrographic nets fish fauna, studies on their potential refuge habitats in drought periods are scarce. Multiannual (2016–2021) research of fish in some streams located in the Saxon Villages area during hydrological drought periods identified, mapped, and revealed the refuge aquatic habitats presence, management needs, and importance for fish diversity and abundance for small rivers. The impact of increasing global temperature and other human activities induced hydrologic net and habitats alteration, decreased the refuge habitats needed by freshwater fish, diminished the fish abundance, and influenced the spatial and temporal variation in fish assemblage structure in the studied area. The sites more than one meter in depth in the studied lotic system were inventoried and all 500 m of these lotic systems were also checked to see what species and how many individuals were present, and if there is was difference in their abundance between refuge and non-refuge 500 m sectors. The scarce number of these refuges due to relatively high soil erosion and clogging in those basins and the cumulative effects of other human types of impact induced a high degree of pressure on the fish fauna. Overall, it reduced the role of these lotic systems as a refuge and for reproduction for the fish of downstream Târnava Mare River, into which all of them flow. Management elements were proposed to maintain and improve these refuges’ ecological support capacity.

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. 

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.