Genetic association with high-resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula

Adaptation and the Geographic Spread of Crop Species

Crops are plant species that were domesticated starting about 11,000 years ago from several centers of origin, most prominently the Fertile Crescent, East Asia, and Mesoamerica. From their domestication centers, these crops spread across the globe and had to adapt to differing environments as a result of this dispersal. We discuss broad patterns of crop spread, including the early diffusion of crops associated with the rise and spread of agriculture, the later movement via ancient trading networks, and the exchange between the Old and New Worlds over the last ∼550 years after the European colonization of the Americas. We also examine the various genetic mechanisms associated with the evolutionary adaptation of crops to their new environments after dispersal, most prominently seasonal adaptation associated with movement across latitudes, as well as altitudinal, temperature, and other environmental factors.

Ensemble species distribution modeling and multilocus phylogeography provide insight into the spatial genetic patterns and distribution dynamics of a keystone forest species, Quercus glauca

BACKGROUND

Forests are essential for maintaining species diversity, stabilizing local and global climate, and providing ecosystem services. Exploring the impact of paleogeographic events and climate change on the genetic structure and distribution dynamics of forest keystone species could help predict responses to future climate change. In this study, we combined an ensemble species distribution model (eSDM) and multilocus phylogeography to investigate the spatial genetic patterns and distribution change of Quercus glauca Thunb, a keystone of East Asian subtropical evergreen broad-leaved forest.

RESULTS

A total of 781 samples were collected from 77 populations, largely covering the natural distribution of Q. glauca. The eSDM showed that the suitable habitat experienced a significant expansion after the last glacial maximum (LGM) but will recede in the future under a general climate warming scenario. The distribution centroid will migrate toward the northeast as the climate warms. Using nuclear SSR data, two distinct lineages split between east and west were detected. Within-group genetic differentiation was higher in the West than in the East. Based on the identified 58 haplotypes, no clear phylogeographic structure was found. Populations in the Nanling Mountains, Wuyi Mountains, and the southwest region were found to have high genetic diversity.

CONCLUSIONS

A significant negative correlation between habitat stability and heterozygosity might be explained by the mixing of different lineages in the expansion region after LGM and/or hybridization between Q. glauca and closely related species. The Nanling Mountains may be important for organisms as a dispersal corridor in the west-east direction and as a refugium during the glacial period. This study provided new insights into spatial genetic patterns and distribution dynamics of Q. glauca.

Landscape genomics reveals adaptive genetic differentiation driven by multiple environmental variables in naked barley on the Qinghai-Tibetan Plateau

Understanding the local adaptation of crops has long been a concern of evolutionary biologists and molecular ecologists. Identifying the adaptive genetic variability in the genome is crucial not only to provide insights into the genetic mechanism of local adaptation but also to explore the adaptation potential of crops. This study aimed to identify the climatic drivers of naked barley landraces and putative adaptive loci driving local adaptation on the Qinghai-Tibetan Plateau (QTP). To this end, a total of 157 diverse naked barley accessions were genotyped using the genotyping-by-sequencing approach, which yielded 3123 high-quality SNPs for population structure analysis and partial redundancy analysis, and 37,636 SNPs for outlier analysis. The population structure analysis indicated that naked barley landraces could be divided into four groups. We found that the genomic diversity of naked barley landraces could be partly traced back to the geographical and environmental diversity of the landscape. In total, 136 signatures associated with temperature, precipitation, and ultraviolet radiation were identified, of which 13 had pleiotropic effects. We mapped 447 genes, including a known gene HvSs1. Some genes involved in cold stress and regulation of flowering time were detected near eight signatures. Taken together, these results highlight the existence of putative adaptive loci in naked barley on QTP and thus improve our current understanding of the genetic basis of local adaptation.

Major chromosome 5H haplotype switch structures the European two-rowed spring barley germplasm of the past 190 years

Selection over 70 years has led to almost complete fixation of a haplotype spanning ~ 250 Mbp of chomosome 5H in European two-rowed spring barleys, possibly originating from North Africa. Plant breeding and selection have shaped the genetic composition of modern crops over the past decades and centuries and have led to great improvements in agronomic and quality traits. Knowledge of the genetic composition of breeding germplasm is essential to make informed decisions in breeding programs. In this study, we characterized the structure and composition of 209 barley cultivars representative of the European two-rowed spring barley germplasm of the past 190 years. Utilizing high-density SNP marker data, we identified a distinct centromeric haplotype spanning a ~ 250 Mbp large region on chromosome 5H which likely was first introduced into the European breeding germplasm in the early to mid-twentieth century and has been non-recombining and under strong positive selection over the past 70 years. Almost all cultivars in our panel that were released after 2000 carry this new haplotype, suggesting that this region carries one or several genes conferring highly beneficial traits. Using the global barley collection of the German Federal ex situ gene bank at IPK Gatersleben, we found the new haplotype at high frequencies in six-rowed spring-type landraces from Northern Africa, from which it may have been introduced into modern European barley germplasm via southern European landraces. The presence of a 250 Mbp genomic region characterized by lack of recombination and high levels of fixation in modern barley germplasm has substantial implications for the genetic diversity of the modern barley germplasm and for barley breeding.

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus)

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future "Geminate Evolutionary Units" (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad (Anaxyrus canorus), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection (R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

Genetic diversity and population structure of barley landraces from Southern Ethiopia's Gumer district: Utilization for breeding and conservation

Barley is the fifth most important food crop in Ethiopia. The genetic relationship and population structure studies of barley are limited to gene bank collections. Therefore, this study fills a gap by investigating the selection, consumption, economic value, genetic diversity, and population structure of farm-collected barley from the Gumer district of the Gurage Zone, which has received little attention. The information on the use of barley in the study area was collected using semi-structured interviews and questionnaires. 124 households of 11 kebeles, the smallest community unit, were interviewed. Barley landraces collected were compared with those collected from Japan, the United States (USA), and other Ethiopian locations. Illumina iSelect (50K genotyping platform) was used to identify single nucleotide polymorphisms (SNP) (20,367). Thirty landraces were found in Gumer. Burdaenadenber had the highest on-farm Shannon index estimate (2.0), followed by Aselecha (1.97) and Enjefo (1.95). Aselecha and Fetazer had the highest (44%) and the lowest (29%) richness values, respectively. High and low Simpson index values were found in Aselecha (84%) and Wulbaragenateretero (79%), respectively. The neighbor-joining tree revealed that Gumer landraces formed a separate subcluster with a common ancestral node; a sister subcluster contained barley landraces from Japan. According to the population structure analysis, barley landraces from Gumer differed from Japan and the United States. The principal component analysis revealed that US barley was the most distant group from Gumer barley. The markers' allele frequencies ranged from 0.10 to 0.50, with an average value of 0.28. The mean values of Nei's gene diversity (0.38) and the polymorphic information content (0.30) indicated the presence of high genetic diversity in the samples. The clustering of accessions was not based on geographic origin. Significant genetic diversity calls for additional research and analysis of local barley diversity because the selection and use of barley in Ethiopia would have been affected by the preference of ethnic groups.

Climatic similarity and genomic background shape the extent of parallel adaptation in Timema stick insects

Evolution can repeat itself, resulting in parallel adaptations in independent lineages occupying similar environments. Moreover, parallel evolution sometimes, but not always, uses the same genes. Two main hypotheses have been put forth to explain the probability and extent of parallel evolution. First, parallel evolution is more likely when shared ecologies result in similar patterns of natural selection in different taxa. Second, parallelism is more likely when genomes are similar because of shared standing variation and similar mutational effects in closely related genomes. Here we combine ecological, genomic, experimental and phenotypic data with Bayesian modelling and randomization tests to quantify the degree of parallelism and its relationship with ecology and genetics. Our results show that the extent to which genomic regions associated with climate are parallel among species of Timema stick insects is shaped collectively by shared ecology and genomic background. Specifically, the extent of genomic parallelism decays with divergence in climatic conditions (that is, habitat or ecological similarity) and genomic similarity. Moreover, we find that climate-associated loci are likely subject to selection in a field experiment, overlap with genetic regions associated with cuticular hydrocarbon traits and are not strongly shaped by introgression between species. Our findings shed light on when evolution is most expected to repeat itself.

Worldwide Selection Footprints for Drought and Heat in Bread Wheat (Triticum aestivum L.)

Genome–environment Associations (GEA) or Environmental Genome-Wide Association scans (EnvGWAS) have been poorly applied for studying the genomics of adaptive traits in bread wheat landraces (Triticum aestivum L.). We analyzed 990 landraces and seven climatic variables (mean temperature, maximum temperature, precipitation, precipitation seasonality, heat index of mean temperature, heat index of maximum temperature, and drought index) in GEA using the FarmCPU approach with GAPIT. Historical temperature and precipitation values were obtained as monthly averages from 1970 to 2000. Based on 26,064 high-quality SNP loci, landraces were classified into ten subpopulations exhibiting high genetic differentiation. The GEA identified 59 SNPs and nearly 89 protein-encoding genes involved in the response processes to abiotic stress. Genes related to biosynthesis and signaling are mainly mediated by auxins, abscisic acid (ABA), ethylene (ET), salicylic acid (SA), and jasmonates (JA), which are known to operate together in modulation responses to heat stress and drought in plants. In addition, we identified some proteins associated with the response and tolerance to stress by high temperatures, water deficit, and cell wall functions. The results provide candidate regions for selection aimed to improve drought and heat tolerance in bread wheat and provide insights into the genetic mechanisms involved in adaptation to extreme environments.

Hybrids Provide More Options for Fine-Tuning Flowering Time Responses of Winter Barley

Crop adaptation requires matching resource availability to plant development. Tight coordination of the plant cycle with prevailing environmental conditions is crucial to maximizing yield. It is expected that winters in temperate areas will become warmer, so the vernalization requirements of current cultivars can be desynchronized with the environment's vernalizing potential. Therefore, current phenological ideotypes may not be optimum for future climatic conditions. Major genes conferring vernalization sensitivity and phenological responses in barley (Hordeum vulgare L.) are known, but some allelic combinations remain insufficiently evaluated. Furthermore, there is a lack of knowledge about flowering time in a hybrid context. To honor the promise of increased yield potentials, hybrid barley phenology must be studied, and the knowledge deployed in new cultivars. A set of three male and two female barley lines, as well as their six F₁ hybrids, were studied in growth chambers, subjected to three vernalization treatments: complete (8 weeks), moderate (4 weeks), and low (2 weeks). Development was recorded up to flowering, and expression of major genes was assayed at key stages. We observed a gradation in responses to vernalization, mostly additive, concentrated in the phase until the initiation of stem elongation, and proportional to the allele constitution and dosage present in VRN-H1. These responses were further modulated by the presence of PPD-H2. The duration of the late reproductive phase presented more dominance toward earliness and was affected by the rich variety of alleles at VRN-H3. Our results provide further opportunities for fine-tuning total and phasal growth duration in hybrid barley, beyond what is currently feasible in inbred cultivars.

Physical geography, isolation by distance and environmental variables shape genomic variation of wild barley (Hordeum vulgare L. ssp. spontaneum) in the Southern Levant

Determining the extent of genetic variation that reflects local adaptation in crop-wild relatives is of interest for the purpose of identifying useful genetic diversity for plant breeding. We investigated the association of genomic variation with geographical and environmental factors in wild barley (Hordeum vulgare L. ssp. spontaneum) populations of the Southern Levant using genotyping by sequencing (GBS) of 244 accessions in the Barley 1K+ collection. The inference of population structure resulted in four genetic clusters that corresponded to eco-geographical habitats and a significant association between lower gene flow rates and geographical barriers, e.g. the Judaean Mountains and the Sea of Galilee. Redundancy analysis (RDA) revealed that spatial autocorrelation explained 45% and environmental variables explained 15% of total genomic variation. Only 4.5% of genomic variation was solely attributed to environmental variation if the component confounded with spatial autocorrelation was excluded. A synthetic environmental variable combining latitude, solar radiation, and accumulated precipitation explained the highest proportion of genomic variation (3.9%). When conditioned on population structure, soil water capacity was the most important environmental variable explaining 1.18% of genomic variation. Genome scans with outlier analysis and genome-environment association studies were conducted to identify adaptation signatures. RDA and outlier methods jointly detected selection signatures in the pericentromeric regions, which have reduced recombination, of the chromosomes 3H, 4H, and 5H. However, selection signatures mostly disappeared after correction for population structure. In conclusion, adaptation to the highly diverse environments of the Southern Levant over short geographical ranges had a limited effect on the genomic diversity of wild barley. This highlighted the importance of nonselective forces in genetic differentiation.

Landscape genomics of the streamside salamander: Implications for species management in the face of environmental change

Understanding spatial patterns of genetic differentiation and local adaptation is critical in a period of rapid environmental change. Climate change and anthropogenic development have led to population declines and shifting geographic distributions in numerous species. The streamside salamander, Ambystoma barbouri, is an endemic amphibian with a small geographic range that predominantly inhabits small, ephemeral streams. As A. barbouri is listed as near‐threatened by the IUCN, we describe range‐wide patterns of genetic differentiation and adaptation to assess the species’ potential to respond to environmental change. We use outlier scans and genetic‐environment association analyses to identify genomic variation putatively underlying local adaptation across the species’ geographic range. We find evidence for adaptation with a polygenic architecture and a set of candidate SNPs that identify genes putatively contributing to local adaptation. Our results build on earlier work that suggests that some A. barbouri populations are locally adapted despite evidence for asymmetric gene flow between the range core and periphery. Taken together, the body of work describing the evolutionary genetics of range limits in A. barbouri suggests that the species may be unlikely to respond naturally to environmental challenges through a range shift or in situ adaptation. We suggest that management efforts such as assisted migration may be necessary in future.

Characterization of the genetic basis of local adaptation of wheat landraces from Iran and Pakistan using genome-wide association study

Characterization of genomic regions underlying adaptation of landraces can reveal a quantitative genetics framework for local wheat (Triticum aestivum L.) adaptability. A collection of 512 wheat landraces from the eastern edge of the Fertile Crescent in Iran and Pakistan were genotyped using genome-wide single nucleotide polymorphism markers generated by genotyping-by-sequencing. The minor allele frequency (MAF) and the heterozygosity (H) of Pakistani wheat landraces (MAF = 0.19, H = 0.008) were slightly higher than the Iranian wheat landraces (MAF = 0.17, H = 0.005), indicating that Pakistani landraces were slightly more genetically diverse. Population structure analysis clearly separated the Pakistani landraces from Iranian landraces, which indicates two separate adaptability trajectories. The large-scale agro-climatic data of seven variables were quite dissimilar between Iran and Pakistan as revealed by the correlation coefficients. Genome-wide association study identified 91 and 58 loci using agroclimatic data, which likely underpin local adaptability of the wheat landraces from Iran and Pakistan, respectively. Selective sweep analysis identified significant hits on chromosomes 4A, 4B, 6B, 7B, 2D, and 6D, which were colocalized with the loci associated with local adaptability and with some known genes related to flowering time and grain size. This study provides insight into the genetic diversity with emphasis on the genetic architecture of loci involved in adaptation to local environments, which has breeding implications.

Major flowering time genes of barley: allelic diversity, effects, and comparison with wheat

This review summarizes the allelic series, effects, interactions between genes and with the environment, for the major flowering time genes that drive phenological adaptation of barley. The optimization of phenology is a major goal of plant breeding addressing the production of high-yielding varieties adapted to changing climatic conditions. Flowering time in cereals is regulated by genetic networks that respond predominately to day length and temperature. Allelic diversity at these genes is at the basis of barley wide adaptation. Detailed knowledge of their effects, and genetic and environmental interactions will facilitate plant breeders manipulating flowering time in cereal germplasm enhancement, by exploiting appropriate gene combinations. This review describes a catalogue of alleles found in QTL studies by barley geneticists, corresponding to the genetic diversity at major flowering time genes, the main drivers of barley phenological adaptation: VRN-H1 (HvBM5A), VRN-H2 (HvZCCTa-c), VRN-H3 (HvFT1), PPD-H1 (HvPRR37), PPD-H2 (HvFT3), and eam6/eps2 (HvCEN). For each gene, allelic series, size and direction of QTL effects, interactions between genes and with the environment are presented. Pleiotropic effects on agronomically important traits such as grain yield are also discussed. The review includes brief comments on additional genes with large effects on phenology that became relevant in modern barley breeding. The parallelisms between flowering time allelic variation between the two most cultivated Triticeae species (barley and wheat) are also outlined. This work is mostly based on previously published data, although we added some new data and hypothesis supported by a number of studies. This review shows the wide variety of allelic effects that provide enormous plasticity in barley flowering behavior, which opens new avenues to breeders for fine-tuning phenology of the barley crop.

Localized environmental heterogeneity drives the population differentiation of two endangered and endemic Opisthopappus Shih species

BACKGROUND

Climate heterogeneity not only indirectly shapes the genetic structures of plant populations, but also drives adaptive divergence by impacting demographic dynamics. The variable localized climates and topographic complexity of the Taihang Mountains make them a major natural boundary in Northern China that influences the divergence of organisms distributed across this region. Opisthopappus is an endemic genus of the Taihang Mountains that includes only two spatially partitioned species Opisthopappus longilobus and Opisthopappus taihangensis. For this study, the mechanisms behind the genetic variations in Opisthopappus populations were investigated.

RESULTS

Using SNP and InDel data coupled with geographic and climatic information, significant genetic differentiation was found to exist either between Opisthopappus populations or two species. All studied populations were divided into two genetic groups with the differentiation of haplotypes between the groups. At approximately 17.44 Ma of the early Miocene, O. taihangensis differentiated from O. longilobus under differing precipitation regimes due to the intensification of the Asian monsoon. Subsequently, intraspecific divergence might be induced by the dramatic climatic transformation from the mid- to late Miocene. During the Pleistocene period, the rapid uplift of the Taihang Mountains coupled with violent climatic oscillations would further promote the diversity of the two species. Following the development of the Taihang Mountains, its complex topography created geographical and ecological heterogeneity, which could lead to spatiotemporal isolation between the Opisthopappus populations. Thus the adaptive divergence might occur within these intraspecific populations in the localized heterogeneous environment of the Taihang Mountains.

CONCLUSIONS

The localized environmental events through the integration of small-scale spatial effects impacted the demographic history and differentiation mechanism of Opisthopappus species in the Taihang Mountains. The results provide useful information for us to understand the ecology and evolution of organisms in the mountainous environment from population and species perspective.

Canonical correlations reveal adaptive loci and phenotypic responses to climate in perennial ryegrass

Germplasm from perennial ryegrass (Lolium perenne L.) natural populations is useful for breeding because of its adaptation to a wide range of climates. Climate-adaptive genes can be detected from associations between genotype, phenotype and climate but an integrated framework for the analysis of these three sources of information is lacking. We used two approaches to identify adaptive loci in perennial ryegrass and their effect on phenotypic traits. First, we combined Genome-Environment Association (GEA) and GWAS analyses. Then, we implemented a new test based on a Canonical Correlation Analysis (CANCOR) to detect adaptive loci. Furthermore, we improved the previous perennial ryegrass gene set by de novo gene prediction and functional annotation of 39,967 genes. GEA-GWAS revealed eight outlier loci associated with both environmental variables and phenotypic traits. CANCOR retrieved 633 outlier loci associated with two climatic gradients, characterized by cold-dry winter versus mild-wet winter and long rainy season versus long summer, and pointed out traits putatively conferring adaptation at the extremes of these gradients. Our CANCOR test also revealed the presence of both polygenic and oligogenic climatic adaptations. Our gene annotation revealed that 374 of the CANCOR outlier loci were positioned within or close to a gene. Co-association networks of outlier loci revealed a potential utility of CANCOR for investigating the interaction of genes involved in polygenic adaptations. The CANCOR test provides an integrated framework to analyse adaptive genomic diversity and phenotypic responses to environmental selection pressures that could be used to facilitate the adaptation of plant species to climate change.

Environmental Association Identifies Candidates for Tolerance to Low Temperature and Drought

Barley (Hordeum vulgare ssp. vulgare) is cultivated from the equator to the Arctic Circle. The wild progenitor species, Hordeum vulgare ssp. spontaneum, occupies a relatively narrow latitudinal range (∼30 - 40° N) primarily at low elevation (< 1,500 m). Adaptation to the range of cultivation has occurred over ∼8,000 years. The genetic basis of adaptation is amenable to study through environmental association. An advantage of environmental association in a well-characterized crop is that many loci that contribute to climatic adaptation and abiotic stress tolerance have already been identified. This provides the opportunity to determine if environmental association approaches effectively identify these loci of large effect. Using published genotyping from 7,864 SNPs in 803 barley landraces, we examined allele frequency differentiation across multiple partitions of the data and mixed model associations relative to bioclimatic variables. Using newly generated resequencing data from a subset of these landraces, we tested for linkage disequilibrium (LD) between SNPs queried in genotyping and SNPs in neighboring loci. Six loci previously reported to contribute to adaptive differences in flowering time and abiotic stress in barley and six loci previously identified in other plant species were identified in our analyses. In many cases, patterns of LD are consistent with the causative variant occurring in the immediate vicinity of the queried SNP. The identification of barley orthologs to well-characterized genes may provide a new understanding of the nature of adaptive variation and could permit a more targeted use of potentially adaptive variants in barley breeding and germplasm improvement.

Genetic association with high-resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula

Landraces are local populations of crop plants adapted to a particular environment. Extant landraces are surviving genetic archives, keeping signatures of the selection processes experienced by them until settling in their current niches. This study intends to establish relationships between genetic diversity of barley (Hordeum vulgare L.) landraces collected in Spain and the climate of their collection sites. A high-resolution climatic data set (5 × 5 km spatial, 1-day temporal grid) was computed from over 2,000 temperature and 7,000 precipitation stations across peninsular Spain. This data set, spanning the period 1981-2010, was used to derive agroclimatic variables meaningful for cereal production at the collection sites of 135 barley landraces. Variables summarize temperature, precipitation, evapotranspiration, potential vernalization and frost probability at different times of the year and time scales (season and month). SNP genotyping of the landraces was carried out combining Illumina Infinium assays and genotyping-by-sequencing, yielding 9,920 biallelic markers (7,479 with position on the barley reference genome). The association of these SNPs with agroclimatic variables was analysed at two levels of genetic diversity, with and without taking into account population structure. The whole data sets and analysis pipelines are documented and available at https://eead-csic-compbio.github.io/barley-agroclimatic-association. We found differential adaptation of the germplasm groups identified to be dominated by reactions to cold temperature and late-season frost occurrence, as well as to water availability. Several significant associations pointing at specific adaptations to agroclimatic features related to temperature and water availability were observed, and candidate genes underlying some of the main regions are proposed.