Warming is expected to lead to drier environments worldwide, especially in the tropics, and it is unclear how crops will react. Drought tolerance often varies at small spatial scales in natural ecosystems, where many of the wild relatives and landraces of the main crops have been collected. Through a series of examples, we will show that collections of wild relatives and landraces, many of those deposited at germplasm banks, may represent this desired source of variation, as they are genetically diverse and phenotypically variable. For instance, using a spectrum of genotyping and phenotyping approaches, we have studied the extent of genetic and phenotypic diversity for drought tolerance in wild and landraces of common bean (Phaseolus vulgaris L.) and compared it with the one available at cultivated varieties. Not surprisingly, most of the naturally available variation to cope with drought in the natural environments was lost through domestication and recent plant breeding. It is therefore imperative to exploit the reservoir of wild relatives and landraces to make crops more tolerant. Yet, it remains to be seen if the rate at which this naturally available variation can be incorporated into the cultivated varieties may keep pace with the rate of climate change.
- drought tolerance
- environmental adaptation
- genomic signatures of selection
- agro-ecological models
- divergent selection
1. Introduction: Common bean – a model to explore the usefulness of wild relatives and landraces as a resource for the future
In this chapter, we review the utility of genome-environment association approaches to infer the potential of wild accessions and landraces to make tropical crops more resistant to climate change, using the food crop common bean (
Common bean is a source of nutrients and protein for over 500 million people in Latin America and Africa, and more than 4.5 out of 23 million hectares are grown in zones where drought is severe, such as in northeastern Brazil, coastal Peru, the central and northern highlands of Mexico, and in Eastern and Southern Africa [9, 10]. This situation may worsen as increased drought due to climate change will reduce global crop production in greater than 10% by 2050 . Increasing drought tolerance in common bean varieties is therefore needed. Characterizing geo-referenced landraces and wild accessions of common bean at the genetic level (e.g., Figure 1) and quantifying single-nucleotide polymorphism (SNP) allelic associations with a bioclimatic-based drought index offer an efficient path to identify adaptive variations suitable to breed new drought-tolerant varieties.
In the following two sections, we fist explain the theoretical bases behind genome-environment associations, as well as its caveats (Section 2), and later we exemplify it with concrete cases that used geo-referenced landraces and wild accessions of common bean to infer naturally available adaptive variations (Section 3).
2. Strategies to infer adaptability of wild relatives and landraces to their natural habitats
Understanding the genomic signatures associated with environmental variation provides insights into how species adapt to their environment [13, 14, 15]. Recent genomic studies in wild populations have demonstrated that genome-environment associations, which are associations between SNP alleles and accessions’ environment of origin, can indeed be used to identify adaptive loci and predict phenotypic variation. For instance, Turner
This genome-environment association approach has also been explored in some crop accessions as a prospection strategy of germplasm, alternative to traditional phenotyping. For example, Yoder
Nonetheless, since genomic signatures associated with habitat heterogeneity can result from causes other than adaptation and selection [24, 25], for example, random genetic drift (Figure 2), and are also influenced by differences in ancestral variation and recombination in the genome [27, 28, 29], some further approaches need to be undertaken to clarify the truthful nature of the divergent regions. For instance, the origin of habitat-associated variants from novel or standing genetic variation leads to distinctively different patterns of genomic divergence [30, 31, 32]. One approach that can help to distinguish these underlying causes of divergence is comparing summary statistics (i.e., Tajima’s D) from different genomic sections because demographic processes usually leave genome-wide signatures while selection tends to imprint more localized regions . Specifically, habitat-mediated purifying selection is associated with localized low values of nucleotide diversity (π)  and Tajima’s D  and high scores of the Watterson’s theta (θ) estimator  because only low-frequency polymorphisms can avoid being eliminated by widespread directional selection. Although recent population bottlenecks tend to achieve the same reduction in nucleotide variation, this pattern is expected at a more genome-wide level. Similarly, local adaptation tends to homogenize haplotypes within the same niche, fix polymorphisms in different populations, and eliminate low-frequency polymorphism. Consequently, few haplotypes with high frequency are retained, corresponding to high values of nucleotide diversity (π) and Tajima’s D and low scores of the Watterson’s theta (θ) estimator . While independent domestication events, extensive population structure, and population expansions after bottlenecks can produce the same patterns, these demographic processes also imprint genomes at a more genome-wide level.
In the following two sub-sections, we explain how to implement genome-environment associations in order to infer the adaptability of wild relatives and landraces to their natural habitats (Section 2.1) and discuss ways to account for causes, other than adaptation and selection, that may be shaping the genomic landscape of signatures associated with habitat heterogeneity (Section 2.2).
2.1. Using genome-environment association scans to identify loci associated with bioclimatic-based indexes
First of all, in order to account for possible demographic effects, subpopulation structure must be determined in geo-referenced landraces and wild accessions using principal coordinates analysis (PCoA) implemented in the software Trait Analysis by aSSociation, Evolution and Linkage, TASSEL v.5 . The same dataset and software can be used to perform association analyses between the SNP markers and bioclimatic-based indexes (e.g. [12, 38, 39]).
As a rule of thumb, a total of ten generalized (GLM) and mixed (MLM) linear models should be compared . Within each model family, five models are usually built as follows: (1) models with the gene pool identity and the first two PCoA axes scores as covariates; (2) models with the within-gene pool subpopulation identity (e.g. ) and the first two PCoA axes scores as covariates; (3) models with the first two PCoA axes scores as covariates; (4) models with the within-gene pool subpopulation identity (e.g. ) as covariates; and (5) models with the gene pool identity as covariates. All five MLMs usually use a centered IBS kinship matrix as a random effect to control for genomic background implementing the EMMA and P3D algorithms to reduce computing time . QQ-plots of the P-values should be inspected to assess whether excessive numbers of false positives are generated and choose in this way the optimum model. Significant associations are determined using strict Bonferroni corrections of P-values at alpha = 0.001, leading, for example, to a significance threshold of 4.4 × 10 in a usual dataset of ca. 23,000 SNP markers (0.001 divided by the number of markers) or −log10(4.4 × 10 ) = 7.36. The construction of customized PCoA and Manhattan diagrams can be carried out with the software R v.3.3.1 (R Core Team).
Finally, candidate genes for habitat adaptation can be identified within the 1000 bp sections, flanking each SNP marker that is associated with a bioclimatic-based index by using the corresponding reference genome (e.g. ) and the PhytoMine and BioMart tools in Phytozome v.12 (phytozome.jgi.doe.gov).
2.2. Accounting for genomic constrains by inspecting genome-wide patterns of variation
In order to identify causes other than adaptation and selection that may be shaping the genomic landscape of signatures associated with habitat heterogeneity (i.e., genomic constrains and genetic drift), sliding window approaches (e.g., window size = 1 × 106 bps, step size = 200 kb) can be implemented to describe patterns of variation and overall divergence across the genome. For instance, SNP density, nucleotide diversity as measured by π , Watterson’s theta (θ) estimator , and Tajima’s D  can be computed using the software TASSEL v.5  and customized R scripts. Results of all windowed analyses are usually plotted against window midpoints in millions of base pairs (Mb) in the software R v.3.3.1 (R Core Team). The centromeres can be marked to visualize the extent of the centromeric repeats and its correlation with overall patterns of diversity and divergence.
It is advisable to calculate bootstrap-based means and 95% confidence intervals around the mean for some summary statistics (i.e., SNP density, π, θ, and Tajima’s D) when computed in sliding windows that contained or did not contain at least one marker that was associated with a bioclimatic-based index. For this, each summary statistic of windows containing and not containing associated SNPs should be randomly resampled with replacement (bootstrapping) across windows within grouping factors (associated vs. not associated). The overall mean is then stored for each grouping factor. This step should be iterated at least 1000 times using customized R scripts. Bootstrapping must be performed independently for each summary statistic in order to eliminate correlations among these.
3. The adaptive potential of wild relatives and landraces in common bean
In common bean, ecological gradients related with drought stress are associated with divergent selection at the genetic level, after accounting for gene pool and subpopulation structure. This divergent selective pressure might be a consequence of local-level rainfall patterns. Specifically, in tropical environments near the equator with bimodal rainfall, a mid-season dry period occurs that can last 2–4 weeks. In contrast, in the sub-tropics, a dry period of 3 or more months can occur. In response to this mid-cycle drought of the sub-tropics,
In the next two sub-sections, we summarize the concrete evidence supporting these statements (Section 3.1) and explain how we can discard other fortuitous causes that may also explain the same pattern (Section 3.2), based on the approaches that we introduced in the previous section (Section 2).
3.1. The signatures of adaptation in common bean are widespread throughout the genome
SNP markers are good at recovering the well-described Andean and Mesoamerican gene pool structure and the five within-gene pool subpopulations observed in wild common bean . Because of this, in a previous research by us with more than 22,000 SNP markers, QQ-plots, from the association analyses between those SNP markers and a bioclimatic-based drought index , indicated that GLM analyses likely had excessive rates of false positives whereas MLM models controlling for population structure and using a kinship matrix reduced more effectively the false positive rate.
In that particular case, the MLM model with the first two PCoA axes scores used as covariates was the best at controlling for false positives. This model yielded a total of 115 SNP markers associated with the bioclimatic-based drought index at a Bonferroni-corrected significance threshold of 7.36 –log10(P-value). These markers explained on average 51.3% ± 0.4 of the variation in the bioclimatic-based drought index. The 115 SNPs were clustered in 90 different regions, defined as overlapping 1000 bp sections that flanked associated markers (Figure 3). Associated SNPs and regions were widespread in all 11 common bean chromosomes.
Following the previous example, chromosomes Pv3 and Pv8 had the highest number of associated SNPs with 21 and 32 SNPs clustered in 16 and 21 different regions, respectively. Chromosomes Pv1, Pv2, Pv4, Pv5, Pv6, and Pv9 contained an intermediate number of associated SNPs with 11, 6, 11, 7, 12, and 9 SNPs clustered in 11, 6, 8, 6, 8, and 9 different regions, respectively. Chromosomes Pv7, Pv10, and Pv11 had the fewest number of associated SNPs with 3, 2, and 1 SNPs clustered in 3, 1, and 1 different regions, respectively. Chromosome Pv8 had more regions with at least 2 associated SNPs than any other chromosome, and these regions had more associated SNPs than in any other chromosome for a total of 5 regions with an average number of associated SNPs of 3.2. The single region that contained more associated SNPs was also situated in chromosome Pv8 with 6 SNPs explaining on average 51.1% ± 0.3 of the variation in the bioclimatic-based drought index. After chromosome Pv8, Pv3 was also outstanding, having 4 regions (with at least 2 associated SNPs) with an average number of associated SNPs of 2.5. Therefore, a total of 75 regions, comprising 99 SNP markers associated with the bioclimatic-based drought index, contained at least 1 gene for a total of 77 genes. Most genes were in chromosomes Pv1, Pv3, and Pv8 with 11, 14, and 16 genes. Only two regions, at chromosomes Pv1 and Pv8 and containing a total of seven different SNPs, spanned two or more genes. The one in Pv8 was the region with more associated SNPs (six in total). One of the two genes in this region encoded an Ankyrin repeat-containing protein, which was associated with osmotic regulation via the assembly of cation channels in the membranes . Among other identified candidate genes, there was a phototropic-responsive NPH3 gene  in Pv3.
3.2. Rampant divergent selection: interpreting genomic signatures of adaptation in common bean beyond genomic constrains
As a follow-up of the previous example, associated genomic windows were enriched for SNP density and positive Tajima’s D scores. This conclusion was achieved after implementing a sliding window analysis to explore the patterns of genome-wide diversity (Figure 3). Marker density decayed drastically toward the centromeres. This decay in diversity proportional to the decay in the rate of recombination was first described in
The average marker density was 44 SNPs per million base pairs (95% CI, 4–143). The average nucleotide diversity as measured by π was 0.3 per million base pairs (95% CI, 0.2–0.4). The average Watterson’s theta (θ) was 0.20 per million base pairs (95% CI, 0.19–0.21). The average Tajima’s D was 0.68 per million base pairs (95% CI, 0.05–1.22). These very same statistics were compared between 1 Mb sliding windows that contained (associated) or did not contain (no associated) at least one marker that was associated with the bioclimatic-based drought index. Genomic windows containing at least one associated SNP had an overall higher SNP density (79 ± 6 vs. 39 ± 2), lower values for Watterson’s theta (θ) scores (0.2016 ± 0.0001 vs. 0.2026 ± 0001), and more positive Tajima’s D scores (0.71 ± 0.02 vs. 0.678 ± 0.009) than windows without associated markers. Nucleotide diversity, as measured by π, was slightly elevated in associated windows when compared with no associated windows (0.322 ± 0.006 vs. 0.317 ± 0.003).
Selective process, such as purifying selection and local adaptation (divergent selection), differentially imprint regions within the same genome, causing a heterogeneous departure of genetic variation from the neutral expectations and from the background trend . Divergent selection tends to homogenize haplotypes within the same niche, fix polymorphisms in different populations, and eliminate low-frequency polymorphism. Consequently, few haplotypes with high frequency are retained, corresponding to high values of nucleotide diversity and Tajima’s D and low scores of the Watterson’s theta (θ) estimator . We have identified these signatures in the various genomic regions associated with a bioclimatic-based drought index. Therefore, it is unlikely that independent domestication events, extensive population structure, and population expansions after bottlenecks are responsible for these patterns because the mixed linear model that we used to identify the genome-environment associations accounted for population structure, while demographic processes would leave genome-wide signatures in both, associated and no-associated windows.
Wild accessions and landraces of common bean occupy more geographical regions with extreme ecologies  and extensive drought stress  than cultivated accessions. Those regions include the arid areas of Peru, Bolivia and Argentina, and the valleys of northwest Mexico. Hence, a broad habitat distribution for wild common bean has exposed these genotypes to both dry and wetter conditions, while cultivated common bean has a narrower distribution and is traditionally considered susceptible to drought. These differences in the ecologies of wild and cultivated common bean have been associated with higher genetic diversity in the former group when surveying candidate genes for drought tolerance such as the ASR , DREB , and ERECTA  gene families, once the population structure  and the background distribution of genetic diversity have been accounted for.
Also, as identified through the genome-environment association approach that was illustrated in this chapter, there are notorious differences between the adaptations of wild accessions and landraces found in arid and more humid environments, in congruence with natural divergent selection acting for thousands of years. Several of these differences might be valuable for plant breeding. Therefore, we reinforce, as was envisioned by Acosta , that wild accessions and landraces of common bean be taken into account to exploit naturally available divergent variations for drought tolerance. We envision that this lesson from common bean will inspire the exploitation of wild relatives and landraces of other crops to face the threats imposed by current climate change.
This chapter ultimately illustrates that genomic signatures of environmental adaptation (e.g. ) are useful for germplasm characterization, potentially enhancing future marker-assisted selection and crop improvement. We envision that genome-environment association studies coupled with estimates of genome-wide diversity will become more common in the upcoming years. These types of studies will likely go beyond estimates of drought tolerance, as exemplified here, to also include estimates regarding frost stress (i.e. [52, 53, 54]), nutrient limitation [55, 56], as well as other threats imposed by climate change [57, 58] in different types of ecosystems (e.g. ) and screened by a variety wide range of genotyping techniques [60, 61, 62, 63]. Genomic selection models  could also incorporate at some point environmental variables in order to improve the prediction of phenotypic variation and the estimation of the genotype-by-environment interactions as well as phenotypic characterizations through novel high-throughput methods such as remote sensing and image analysis , and novel models such as genome-environment associations  in the light of linkage disequilibrium (LD)  and various stochastic approaches [68, 69].
We are grateful with Daniel G. Debouck and the Genetic Resources Unit at the International Center for Tropical Agriculture for donating the seeds and sharing the geographic coordinates of the plant material that was referred to in this study. This research was supported by the Lundell and Tullberg grants to AJC. The Geneco mobility fund from Lund University is acknowledged for making possible the synergistic meeting between AJC and MWB in the Spring of 2015 at Nashville, Tennessee. Writing time was possible thanks to the grants CTS14:55 from the Carl Trygger Foundation to S. Berlin and 4.1-2016-00418 from Vetenskapsrådet (VR) to AJC.
Bitocchi E, et al. Mesoamerican origin of the common bean ( Phaseolus vulgarisL.) is revealed by sequence data. Proceedings of the National Academy of Sciences. 2012; 109(14):E788-E796
Rodriguez M, et al. Landscape genetics, adaptive diversity and population structure in Phaseolus vulgaris. The New Phytologist. 2016; 209(4):1781-1794
Gepts P, Debouck D. Origin, domestication and evolution of the common bean ( Phaseolus vulgarisL.). In: Van Shoonhoven A, Voysest O, editors. Common Beans: Research for Crop Improvement. Commonwealth Agricultural Bureau: Wallingford; 1991. p. 7-53
Kwak M, Gepts P. Structure of genetic diversity in the two major gene pools of common bean ( Phaseolus vulgarisL., Fabaceae). Theoretical and Applied Genetics. 2009; 118(5):979-992
Schmutz J, et al. A reference genome for common bean and genome-wide analysis of dual domestications. Nature Genetics. 2014; 46(7):707-713
Bitocchi E, et al. Molecular analysis of the parallel domestication of the common bean ( Phaseolus vulgaris) in Mesoamerica and the Andes. New Phytologist. 2013; 197(1):300-313
Cortés AJ. On the origin of the common bean ( Phaseolus vulgarisL.). American Journal of Plant Sciences. 2013; 4(10):1998-2000
Asfaw A, Blair MW, Almekinders C. Genetic diversity and population structure of common bean ( Phaseolus vulgarisL.) landraces from the East African highlands. Theoretical and Applied Genetics. 2009; 120(1-12)
Broughton WJ, et al. Beans ( Phaseolusspp.) – Model food legumes. Plant and Soil. 2003; 252(1):55-128
Darkwa K, et al. Evaluation of common bean ( Phaseolus vulgarisL.) genotypes for drought stress adaptation in Ethiopia. The Crop Journal. 2016; 4(5):367-376
Tai APK, Martin MV, Heald CL. Threat to future global food security from climate change and ozone air pollution. Nature Climate Change. 2014; 4:817-821
Cortés AJ, et al. Drought tolerance in wild plant populations: The case of common beans ( Phaseolus vulgarisL.). PLoS One. 2013; 8(5):e62898
Hoffmann AA, Sgro CM. Climate change and evolutionary adaptation. Nature. 2011; 470(7335):479-485
Stinchcombe JR, Hoekstra HE. Combining population genomics and quantitative genetics: Finding the genes underlying ecologically important traits. Heredity. 2008; 100:158-170
Savolainen O, Lascoux M, Merilä J. Ecological genomics of local adaptation. Nature Reviews Genetics. 2013; 14(11):807-820
Turner TL, et al. Population resequencing reveals local adaptation of Arabidopsis lyratato serpentine soils. Nature Genetics. 2010; 42(3):260-263
Hancock AM, et al. Adaptation to climate across the Arabidopsis thalianagenome. Science. 2011; 334(6052):83-86
Fischer MC, et al. Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis hallerifrom the Alps. Molecular Ecology. 2013; 22(22):5594-5607
Pluess AR, et al. Genome-environment association study suggests local adaptation to climate at the regional scale in Fagus sylvatica. The New Phytologist. 2016; 210(2):589-601
Yeaman S, et al. Convergent local adaptation to climate in distantly related conifers. Science. 2016; 353(6306):1431-1433
Yoder JB, et al. Genomic signature of adaptation to climate in Medicago truncatula. Genetics. 2014; 196(4):1263-1275
Lasky JR, et al. Genome-environment associations in sorghum landraces predict adaptive traits. Science Advances. 2015; 1(6):e1400218
Berthouly-Salazar C, et al. Genome scan reveals selection acting on genes linked to stress response in wild pearl millet. Molecular Ecology. 2016; 25(21):5500-5512
Nei M. The neutral theory of molecular evolution in the genomic era. Annual Review of Genomics and Human Genetics. 2010; 11:265-289
Cortés AJ, et al. Genome-wide patterns of microhabitat-driven divergence in the alpine dwarf shrub Salix herbaceaL. In: On The Big Challenges of a Small Shrub: Ecological Genetics of Salix herbaceaL. Uppsala: Acta Universitatis Upsaliensis; 2015
Ravinet M, et al. Interpreting the genomic landscape of speciation: A road map for finding barriers to gene flow. Journal of Evolutionary Biology. 2017; 30(8):1450-1477
Strasburg JL, et al. What can patterns of differentiation across plant genomes tell us about adaptation and speciation? Philosophical Transactions of the Royal Society B: Biological Sciences. 2011; 367(1587):364-373
Ellegren H, Galtier N. Determinants of genetic diversity. Nature Reviews. Genetics. 2016; 17(7):422-433
Wolf JB, Ellegren H. Making sense of genomic islands of differentiation in light of speciation. Nature Reviews Genetics. 2017; 18:87-100
Barrett RD, Schluter D. Adaptation from standing genetic variation. Trends in Ecology & Evolution. 2008; 23(1):38-44
Pritchard JK, Pickrell JK, Coop G. The genetics of human adaptation: Hard sweeps, soft sweeps, and polygenic adaptation. Current Biology. 2010; 20(4):R208-R215
Hermisson J, Pennings PS. Soft sweeps: Molecular population genetics of adaptation from standing genetic variation. Genetics. 2005; 169(4):2335-2352
Wakeley J. Coalescent Theory: An Introduction. Cambridge: Harvard University; 2008
Nei M. Molecular Evolutionary Genetics. New York: Columbia University Press; 1987
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989; 123(3):585-595
Watterson GA. Number of segregating sites in genetic models without recombination. Theoretical Population Biology. 1975; 7:256-276
Bradbury PJ, et al. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23(19):2633-2635
Cortés AJ. Environmental Heterogeneity at a Fine Scale: Ecological and Genetic Implications in a Changing World. Vol. 60. Saarbrücken, Germany: LAP Lambert Academic Publishing; 2016
Cortés AJ. Local scale genetic diversity and its role in coping with changing climate. In: Bitz L, editor. Genetic Diversity. Rijeka: InTech; 2017. 140 p
Cortés AJ, et al. Maintenance of female-bias in a polygenic sex determination system is consistent with genomic conflict. In: On The Big Challenges of a Small Shrub: Ecological Genetics of Salix herbaceaL. Uppsala: Acta Universitatis Upsaliensis; 2015
Blair MW, Soler A, Cortés AJ. Diversification and population structure in common beans ( Phaseolus vulgarisL.). PLoS One. 2012; 7(11):e49488
Zhang Z, et al. Mixed linear model approach adapted for genome-wide association studies. Nature Genetics. 2010; 42(4):355-360
Beebe S, et al. Selection for drought resistance in common bean also improves yield in phosphorus limited and Favorable environments. Crop Science. 2008; 48:582-592
Cortés AJ, et al. Small-scale patterns in snowmelt timing affect gene flow and the distribution of genetic diversity in the alpine dwarf shrub Salix herbacea. Heredity. 2014; 113:233-239
Voronin DA, Kiseleva EV. Functional role of proteins containing ankyrin repeats. Cell and Tissue Biology. 2008; 2(1):1-12
Pedmale UV, Liscum E. Regulation of phototropic signaling in Arabidopsis via phosphorylation state changes in the phototropin 1-interacting protein NPH3. Journal of Biological Chemestry. 2007; 282:19992-20001
Cortés AJ, et al. Molecular ecology and selection in the drought-related Asr gene polymorphisms in wild and cultivated common bean ( Phaseolus vulgarisL.). BMC Genetics. 2012; 13:58
Cortés AJ, et al. Nucleotide diversity patterns at the drought-related DREB2 encoding genes in wild and cultivated common bean ( Phaseolus vulgarisL.). Theoretical and Applied Genetics. 2012; 125(5):1069-1085
Blair MW, Cortés AJ, This D. Identification of an ERECTA gene and its drought adaptation associations with wild and cultivated common bean. Plant Science. 2016; 242:250-259
Acosta JA, Kelly JD, Gepts P. Prebreeding in common bean and use of genetic diversity from wild germplasm. Crop Science. 2007; 47(S3):S44-S59
Cortés AJ. On the Big Challenges of a Small Shrub: Ecological Genetics of Salix herbaceaL.Uppsala: Acta Universitatis Upsaliensis; 2015
Wheeler JA, et al. The snow and the willows: accelerated spring snowmelt reduces performance in the low-lying alpine shrub Salix herbacea. Journal of Ecology. 2016; 104(4):1041-1050
Wheeler JA, et al. Increased spring freezing vulnerability for alpine shrubs under early snowmelt. Oecologia. 2014; 175(1):219-229
Wheeler JA, et al. With a little help from my friends: Community facilitation increases performance in the dwarf shrub Salix herbacea. Basic and Applied Ecology. 2015; 16:202-209
Little CJ, et al. Small-scale drivers: The importance of nutrient availability and snowmelt timing on performance of the alpine shrub Salix herbacea. Oecologia. 2016; 180(4):1015-1024
Sedlacek J, et al. What role do plant-soil interactions play in the habitat suitability and potential range expansion of the alpine dwarf shrub Salix herbacea? Basic and Applied Ecology. 2014; 15(4):305-315
Sedlacek J, et al. Evolutionary potential in the alpine: Trait heritabilities and performance variation of the dwarf willow Salix herbaceafrom different elevations and microhabitats. Ecology and Evolution. 2016; 6(12):3940-3952
Sedlacek J, et al. The response of the alpine dwarf shrub Salix herbaceato altered snowmelt timing: Lessons from a multi-site transplant experiment. PLoS One. 2015; 10(4):e0122395
Madriñán S, Cortés AJ, Richardson JE. Páramo is the world’s fastest evolving and coolest biodiversity hotspot. Frontiers in Genetics. 2013; 4(192)
Blair MW, et al. A high-throughput SNP marker system for parental polymorphism screening, and diversity analysis in common bean ( Phaseolus vulgarisL.). Theoretical and Applied Genetics. 2013; 126(2):535-548
Kelleher CT, et al. SNP discovery, gene diversity, and linkage disequilibrium in wild populations of Populus tremuloides. Tree Genetics & Genomes. 2012; 8(4):821-829
Galeano CH, et al. Gene-based single nucleotide polymorphism markers for genetic and association mapping in common bean. BMC Genetics. 2012; 13(1):48
Cortés AJ, Chavarro MC, Blair MW. SNP marker diversity in common bean ( Phaseolus vulgarisL.). Theoretical and Applied Genetics. 2011; 123(5):827-845
Desta ZA, Ortiz R. Genomic selection: Genome-wide prediction in plant improvement. Trends in Plant Science. 2014; 19(9):592-601
Lo Bianco M, Grillo O, Sarigu M, Cremonini R, Venora G. Characterisation of Italian bean landraces ( Phaseolus vulgarisL.) using seed image analysis and texture descriptors. Australian Journal of Crop Science. 2015; 9(11):1022-1034
Cortés AJ, Blair MW. Genotyping by sequencing and genome – environment associations in wild common bean predict widespread divergent adaptation to drought. Frontiers in Plant Science. 2018; 9(128)
Blair MW, Cortés AJ, Farmer A, Assefa T, Penmetsa RV, Cook D, et al. Uneven recombination rate and linkage disequilibrium across a reference SNP map for common bean. PLOS ONE. 2018; DOI: pone.0189597
Cortés AJ. Prevalence in MSM is enhanced by role versatility. In: Mazari A, editor. Big Data Analytics in HIV/AIDS Research. Hershey: IGI Global; 2018
Cortés AJ. On how role versatility boosts an STI. Journal of Theoretical Biology. 2017; 440:66-69