Genome–environment associations (GEAs) are a powerful strategy for the study of adaptive traits in wild plant populations, yet they still lack behind in the use of modern statistical methods as the ones suggested for genome-wide association studies (GWASs). In order to bridge this gap, we couple GEA with last-generation GWAS algorithms in common bean to identify novel sources of heat tolerance across naturally heterogeneous ecosystems. Common bean (Phaseolus vulgaris L.) is the most important legume for human consumption, and breeding it for resistance to heat stress is key because annual increases in atmospheric temperature are causing decreases in yield of up to 9% for every 1°C. A total of 78 geo-referenced wild accessions, spanning the two gene pools of common bean, were genotyped by sequencing (GBS), leading to the discovery of 23,373 single-nucleotide polymorphism (SNP) markers. Three indices of heat stress were developed for each accession and inputted in last-generation algorithms (i.e. SUPER, FarmCPU, and BLINK) to identify putative associated loci with the environmental heterogeneity in heat stress. Best-fit models revealed 120 significantly associated alleles distributed in all 11 common bean chromosomes. Flanking candidate genes were identified using 1-kb genomic windows centered in each associated SNP marker. Some of these genes were directly linked to heat-responsive pathways, such as the activation of heat shock proteins (MED23, MED25, HSFB1, HSP40, and HSP20). We also found protein domains related to thermostability in plants such as S1 and Zinc finger A20 and AN1. Other genes were related to biological processes that may correlate with plant tolerance to high temperature, such as time to flowering (MED25, MBD9, and PAP), germination and seedling development (Pkinase_Tyr, Ankyrin-B, and Family Glicosil-hydrolase), cell wall stability (GAE6), and signaling pathway of abiotic stress via abscisic acid (histone-like transcription factors NFYB and phospholipase C) and auxin (Auxin response factor and AUX_IAA). This work offers putative associated loci for marker-assisted and genomic selection for heat tolerance in common bean. It also demonstrates that it is feasible to identify genome-wide environmental associations with modest sample sizes by using a combination of various carefully chosen environmental indices and last-generation GWAS algorithms.

Warming and drought are reducing global crop production with a potential to substantially worsen global malnutrition. As with the green revolution in the last century, plant genetics may offer concrete opportunities to increase yield and crop adaptability. However, the rate at which the threat is happening requires powering new strategies in order to meet the global food demand. In this review, we highlight major recent ‘big data’ developments from both empirical and theoretical genomics that may speed up the identification, conservation, and breeding of exotic and elite crop varieties with the potential to feed humans. We first emphasize the major bottlenecks to capture and utilize novel sources of variation in abiotic stress (i.e., heat and drought) tolerance. We argue that adaptation of crop wild relatives to dry environments could be informative on how plant phenotypes may react to a drier climate because natural selection has already tested more options than humans ever will. Because isolated pockets of cryptic diversity may still persist in remote semi-arid regions, we encourage new habitat-based population-guided collections for genebanks. We continue discussing how to systematically study abiotic stress tolerance in these crop collections of wild and landraces using geo-referencing and extensive environmental data. By uncovering the genes that underlie the tolerance adaptive trait, natural variation has the potential to be introgressed into elite cultivars. However, unlocking adaptive genetic variation hidden in related wild species and early landraces remains a major challenge for complex traits that, as abiotic stress tolerance, are polygenic (i.e., regulated by many low-effect genes). Therefore, we finish prospecting modern analytical approaches that will serve to overcome this issue. Concretely, genomic prediction, machine learning, and multi-trait gene editing, all offer innovative alternatives to speed up more accurate pre- and breeding efforts toward the increase in crop adaptability and yield, while matching future global food demands in the face of increased heat and drought. In order for these ‘big data’ approaches to succeed, we advocate for a trans-disciplinary approach with open-source data and long-term funding. The recent developments and perspectives discussed throughout this review ultimately aim to contribute to increased crop adaptability and yield in the face of heat waves and drought events.

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations’ responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder’s equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.

Some of the major impacts of climate change are expected in regions where drought stress is already an issue. Grain legumes are generally drought susceptible. However, tepary bean and its wild relatives within Phaseolus acutifolius or P. parvifolius are from arid areas between Mexico and the United States. Therefore, we hypothesize that these bean accessions have diversity signals indicative of adaptation to drought at key candidate genes such as: Asr2, Dreb2B, and ERECTA. By sequencing alleles of these genes and comparing to estimates of drought tolerance indices from climate data for the collection site of geo-referenced, tepary bean accessions, we determined the genotype x environmental association (GEA) of each gene. Diversity analysis found that cultivated and wild P. acutifolius were intermingled with var. tenuifolius and P. parvifolius, signifying that allele diversity was ample in the wild and cultivated clade over a broad sense (sensu lato) evaluation. Genes Dreb2B and ERECTA harbored signatures of directional selection, represented by six SNPs correlated with the environmental drought indices. This suggests that wild tepary bean is a reservoir of novel alleles at genes for drought tolerance, as expected for a species that originated in arid environments. Our study corroborated that candidate gene approach was effective for marker validation across a broad genetic base of wild tepary accessions.

Heat and drought are major stresses that significantly reduce seed yield of the common bean (Phaseolus vulgaris L.). In turn, this affects the profitability of the crop in climatic-vulnerable tropical arid regions, which happen to be the poorest and in most need of legume proteins. Therefore, it is imperative to broaden the sources of heat and drought resistance in the common bean by examining closely related species from warmer and drier environments (i.e., Tepary bean, P. acutifolius A. Gray), while harnessing such variation, typically polygenic, throughout advanced interspecific crossing schemes. As part of this study, interspecific congruity backcrosses for high temperature and drought tolerance conditions were characterized across four localities in coastal Colombia. Genotypes with high values of CO2 assimilation (>24 µmol CO2 m−2 s−1), promising yield scores (>19 g/plant), and high seed mineral content (Fe > 100 mg/kg) were identified at the warmest locality, Motilonia. At the driest locality, Caribia, one intercrossed genotype (i.e., 85) and the P. acutifolius G40001 control exhibited sufficient yield for commercial production (17.76 g/plant and 12.76 g/plant, respectively). Meanwhile, at southernmost Turipaná and Carmen de Bolívar localities, two clusters of genotypes exhibited high mean yield scores with 33.31 g/plant and 17.89 g/plant, respectively, and one genotype had an increased Fe content (109.7 mg/kg). Overall, a multi-environment AMMI analysis revealed that genotypes 13, 27, 82, and 84 were environmentally stable with higher yield scores compared to the Tepary control G40001. Ultimately, this study allows us to conclude that advanced common bean × Tepary bean interspecific congruity backcrosses are capable of pyramiding sufficient polygenic tolerance responses for the extreme weather conditions of coastal Colombia, which are likely to worsen due to climate change. Furthermore, some particular recombination events (i.e., genotype 68) show that there may be potential to couple breeding for heat and drought tolerance with Fe mineral biofortification, despite a prevalent trade-off, as a way to fight malnutrition of marginalized communities in tropical regions.

Leveraging innovative tools to speed up prebreeding and discovery of genotypic sources of adaptation from landraces, crop wild relatives, and orphan crops is a key prerequisite to accelerate genetic gain of abiotic stress tolerance in annual crops such as legumes and cereals, many of which are still orphan species despite advances in major row crops. Here, we review a novel, interdisciplinary approach to combine ecological climate data with evolutionary genomics under the paradigm of a new field of study: genome–environment associations (GEAs). We first exemplify how GEA utilizes in situ georeferencing from genotypically characterized, gene bank accessions to pinpoint genomic signatures of natural selection. We later discuss the necessity to update the current GEA models to predict both regional- and local- or micro-habitat–based adaptation with mechanistic ecophysiological climate indices and cutting-edge GWAS-type genetic association models. Furthermore, to account for polygenic evolutionary adaptation, we encourage the community to start gathering genomic estimated adaptive values (GEAVs) for genomic prediction (GP) and multi-dimensional machine learning (ML) models. The latter two should ideally be weighted by de novo GWAS-based GEA estimates and optimized for a scalable marker subset. We end the review by envisioning avenues to make adaptation inferences more robust through the merging of high-resolution data sources, such as environmental remote sensing and summary statistics of the genomic site frequency spectrum, with the epigenetic molecular functionality responsible for plastic inheritance in the wild. Ultimately, we believe that coupling evolutionary adaptive predictions with innovations in ecological genomics such as GEA will help capture hidden genetic adaptations to abiotic stresses based on crop germplasm resources to assist responses to climate change.“I shall endeavor to find out how nature’s forces act upon one another, and in what manner the geographic environment exerts its influence on animals and plants. In short, I must find out about the harmony in nature” Alexander von Humboldt—Letter to Karl Freiesleben, June 1799.

Mint (Mentha L., Lamiaceae) is a strongly scented herb of the family Lamiaceae that is grown mostly by clonal propagation, making it a valuable species for the study of somaclonal variation and its phenotypic consequences. The recent introduction of a few species of mint in South America, followed by a presumably rampant propagation, make this region particularly ideal for studying the extent of somaclonal genetic diversity. Hence, the objective of this work was to offer a preliminary characterization of somaclonal genetically coding diversity of the mint in the northern Andes in order to address the question of whether somaclonal variants may have emerged despite relatively recent introductions in a region where mint is not native. A total of 29 clonally propagated specimens, collected in mint export farms in the province of Antioquia, a major region for mint production in the northwest Andes of Colombia, were genotyped using RNA sequencing (RNA-Seq). SNP calling was carried out from the leaves’ transcriptome profiles of each plant by combining the GATK4 and TRINITY protocols, obtaining a total of 2033 loci across 912 transcripts with a minimum read depth of 20X and 4% of missing data. Unsupervised machine learning algorithms considered the K-means, AGNES and UPGMA approaches, all of which suggested three genetic clusters for M. spicata and a unique cluster for M. × piperita. The results indicate that at least two different origins of M. spicata reached the eastern region of the Antioquia province, clonally propagated in the locality ever since for local consumption and export. One of these ancestries had more population structure, possibly due to environmental or anthropological pressures that intervened in the fragmentation of this genetic group or to a higher somaclonal mutation rate. This work offers a first step into the study of the accumulation and transmission of presumably quasi-neutral somatic mutations at coding regions in an herbaceous clonally propagated scented species such as mint, likely favored by an expected population expansion after its Andean introduction. These ad hoc hypotheses warrant further study as part of future research.

Grafting induces precocity and maintains clonal integrity in fruit tree crops. However, the complex rootstock × scion interaction often precludes understanding how the tree phenotype is shaped, limiting the potential to select optimum rootstocks. Therefore, it is necessary to assess (1) how seedling progenies inherit trait variation from elite ‘plus trees’, and (2) whether such family superiority may be transferred after grafting to the clonal scion. To bridge this gap, we quantified additive genetic parameters (i.e., narrow sense heritability—h2, and genetic-estimated breeding values—GEBVs) across landraces, “criollo”, “plus trees” of the super-food fruit tree crop avocado (Persea americanaMill.), and their open-pollinated (OP) half-sib seedling families. Specifically, we used a genomic best linear unbiased prediction (G-BLUP) model to merge phenotypic characterization of 17 morpho-agronomic traits with genetic screening of 13 highly polymorphic SSR markers in a diverse panel of 104 avocado “criollo” “plus trees.” Estimated additive genetic parameters were validated at a 5-year-old common garden trial (i.e., provenance test), in which 22 OP half-sib seedlings from 82 elite “plus trees” served as rootstocks for the cv. Hass clone. Heritability (h2) scores in the “criollo” “plus trees” ranged from 0.28 to 0.51. The highesth2values were observed for ribbed petiole and adaxial veins with 0.47 (CI 95%0.2–0.8) and 0.51 (CI 0.2–0.8), respectively. Theh2scores for the agronomic traits ranged from 0.34 (CI 0.2–0.6) to 0.39 (CI 0.2–0.6) for seed weight, fruit weight, and total volume, respectively. When inspecting yield variation across 5-year-old grafted avocado cv. Hass trees with elite OP half-sib seedling rootstocks, the traits total number of fruits and fruits’ weight, respectively, exhibitedh2scores of 0.36 (± 0.23) and 0.11 (± 0.09). Our results indicate that elite “criollo” “plus trees” may serve as promissory donors of seedling rootstocks for avocado cv. Hass orchards due to the inheritance of their outstanding trait values. This reinforces the feasibility to leverage natural variation from “plus trees”viaOP half-sib seedling rootstock families. By jointly estimating half-sib family effects and rootstock-mediated heritability, this study promises boosting seedling rootstock breeding programs, while better discerning the consequences of grafting in fruit tree crops.