Genomic evidence for climate-linked diversity loss and increased vulnerability of wild barley spanning 28 years of climate warming

The information on how plant populations respond genetically to climate warming is scarce. Here, landscape genomic and machine learning approaches were integrated to assess genetic response of 10 wild barley (Hordeum vulgare ssp. spontaneum; WB) populations in the past and future, using whole genomic sequencing (WGS) data. The WB populations were sampled in 1980 and again in 2008. Phylogeny of accessions was roughly in conformity with sampling sites, which accompanied by admixture/introgressions. The 28-y climate warming resulted in decreased genetic diversity, increased selection pressure, and an increase in deleterious single nucleotide polymorphism (dSNP) numbers, heterozygous deleterious and total deleterious burdens for WB. Genome-environment associations identified some candidate genes belonging to peroxidase family (HORVU2Hr1G057450, HORVU4Hr1G052060 and HORVU4Hr1G057210) and heat shock protein 70 family (HORVU2Hr1G112630). The gene HORVU2Hr1G120170 identified by selective sweep analysis was under strong selection during the climate warming of the 28-y, and its derived haplotypes were fixed by WB when faced with the 28-y increasingly severe environment. Temperature variables were found to be more important than precipitation variables in influencing genomic variation, with an eco-physiological index gdd5 (growing degree-days at the baseline threshold temperature of 5 °C) being the most important determinant. Gradient forest modelling revealed higher predicted genomic vulnerability in Sede Boqer under future climate scenarios at 2041-2070 and 2071-2100. Additionally, estimates of effective population size (Ne) tracing back to 250 years indicated a forward decline in all populations over time. Our assessment about past genetic response and future vulnerability of WB under climate warming is crucial for informing conservation efforts for wild cereals and rational use strategies.