Original Research ARTICLE
Ecogeography and utility to plant breeding of the crop wild relatives of sunflower (Helianthus annuus L.)
- 1Biodiversity Research Centre and Department of Botany, University of British Columbia, Vancouver, BC, Canada
- 2Department of Agronomy and Plant Genetics, University of Minnesota, St. Paul, MN, USA
- 3International Center for Tropical Agriculture, Cali, Colombia
- 4Centre for Crop Systems Analysis, Wageningen University, Wageningen, Netherlands
- 5School of Biosciences, University of Birmingham, Birmingham, UK
- 6Department of Horticulture and Crop Science, The Ohio State University, Columbus, OH, USA
- 7Department of Ecology and Evolutionary Biology, University of Colorado at Boulder, Boulder, CO, USA
- 8Agronomy Department, North Central Regional Plant Introduction Station, Iowa State University and United States Department of Agriculture Agricultural Research Service, Ames, IA, USA
- 9Northern Crop Science Laboratory, United States Department of Agriculture Agricultural Research Service, Fargo, ND, USA
- 10Department of Biology, Indiana University, Bloomington, IN, USA
Crop wild relatives (CWR) are a rich source of genetic diversity for crop improvement. Combining ecogeographic and phylogenetic techniques can inform both conservation and breeding. Geographic occurrence, bioclimatic, and biophysical data were used to predict species distributions, range overlap and niche occupancy in 36 taxa closely related to sunflower (Helianthus annuus L.). Taxa lacking comprehensive ex situ conservation were identified. The predicted distributions for 36 Helianthus taxa identified substantial range overlap, range asymmetry and niche conservatism. Specific taxa (e.g., Helianthus deblis Nutt., Helianthus anomalus Blake, and Helianthus divaricatus L.) were identified as targets for traits of interest, particularly for abiotic stress tolerance, and adaptation to extreme soil properties. The combination of techniques demonstrates the potential for publicly available ecogeographic and phylogenetic data to facilitate the identification of possible sources of abiotic stress traits for plant breeding programs. Much of the primary genepool (wild H. annuus) occurs in extreme environments indicating that introgression of targeted traits may be relatively straightforward. Sister taxa in Helianthus have greater range overlap than more distantly related taxa within the genus. This adds to a growing body of literature suggesting that in plants (unlike some animal groups), geographic isolation may not be necessary for speciation.
Plant genetic resources represent the biological foundation for maintaining and improving crop productivity having played a central role in crop development from antiquity (Porter et al., 2014). Crop wild relatives (CWR) are an important source of useful traits for plant breeding (Hoisington et al., 1999; Hajjar and Hodgkin, 2007). With the world's population projected to increase the need to produce more food while using fewer natural resource inputs under increasingly stochastic climatic conditions is a major challenge (Butler and Huybers, 2013; Challinor et al., 2014). CWR conservation and utilization focusing on the use of improving technologies (high throughput phenotyping, genotyping, and geographical information systems), has been proposed as a way to acquire a greater knowledge of conservation needs and lead to more targeted use of CWR germplasm (Khoury et al., 2010; Cabrera-Bosquet et al., 2012; McCouch et al., 2013). Targeted collecting for ex situ conservation has become a priority as rapid changes in both climate and land use patterns increasingly threaten CWR in their natural habitats (Jarvis et al., 2008; McCouch et al., 2013).
CWR have traditionally been categorized based on crossing relationships with domesticates; the primary germplasm contains no crossing barriers, the secondary contains some meiotic abnormalities, and the tertiary requires special techniques such as embryo rescue (Harlan and de Wet, 1971; Harlan, 1976). Such classifications may be supplemented by molecular, bioclimatic, and biophysical data to aid in the identification of candidate taxa for breeding, although such efforts have been constrained by challenges in comprehensively generating and integrating these data (Ricklefs and Jenkins, 2011).
The genus Helianthus L. contains 52 species comprising 67 taxa (Schilling, 2006; Stebbins et al., 2013). Native to North America, the taxa occupy a variety of habitats ranging from open plains to salt marshes (Seiler and Marek, 2011; Kane et al., 2013). Sunflower (Helianthus annuus L.) is the most economically important species from the genus, with ~26 million hectares in production worldwide and a substantial private sector breeding effort, particularly for oil production (FAOSTAT, 2013). Domesticated approximately 4000 years ago in east central North America, sunflower has a typical domestication syndrome; i.e., it does not branch, does not have seed dormancy, has a predictable flowering time, and does not shatter (Harlan et al., 1973; Harter et al., 2004; Blackman et al., 2011). The crop has undergone both selection and genetic drift during domestication and improvement, which has reduced genetic diversity (Tang and Knapp, 2003; Liu and Burke, 2006), with modern cultivars retaining 50–67% of the diversity present in wild H. annuus populations (Kolkman et al., 2007; Mandel et al., 2011).
Sunflower has often utilized CWR in breeding efforts, with many of the taxa hybridizing well with the crop (Table S1; Table 1) (Long, 1960; Chandler et al., 1986). Despite the historical use, CWR of sunflower are considered to be relatively untapped, particularly in regard to adaptation to abiotic stresses. To contribute to an enhanced understanding of the CWR of sunflower, this studies' objectives were to (1) create geographical distribution models for 36 CWR taxa, and (2) explore niche habitation through comparisons of ecogeographic and phylogenetic data, to identify taxa occurring in extreme environments of potential interest to sunflower breeding.
Table 1. Taxa examined in this study, recommendation, position in germplasm, environmental cluster, life history, and potential extreme characteristics.
Materials and Methods
Species Distribution Modeling
A modified gap analysis (Ramírez-Villegas et al., 2010) was used to determine the conservation status of 36 taxa within Helianthus selected based upon their potential to provide useful traits for sunflower breeding. Briefly, (1) target taxa were identified, and geographic occurrence data were gathered and verified, (2) the overall representation of CWR in germplasm collections was estimated, (3) potential distribution models were produced for taxa with sufficient samples with coordinates, (4) the geographic and ecological representation of germplasm collections were assessed for each taxon by comparing potential distribution models to existing germplasm collection locations, (5) taxa were prioritized for further collecting based upon the average of their overall, geographic, and ecological coverage results, and (6) gap analysis results were correlated with the subjective assessments of collection priorities from crop experts.
The selection of taxa for analysis was based on membership within the primary or secondary genepools of sunflower (Vincent et al., 2013) with the addition of all taxa from the tertiary genepool indicated in publications to be confirmed or potential trait donors (Table S1). A total of 12,737 occurrence records for the 36 taxa, sourced from 31 herbaria and five genebanks, were used for distribution models and conservation analysis (Table S2), including 4705 records with geographic coordinates. The overall representation of taxa in genebank collections was estimated using the “Sampling Representativeness Score” (SRS), calculated as the number of germplasm samples (GS) divided by the total number of samples (GS plus reference records). After eliminating duplicate records, potential distributions were calculated using Maxent (Phillips et al., 2006), with a k-5 cross-validation option and 10,000 background points for model training over North America (Phillips, 2008; VanDerWal et al., 2009). We included 19 bioclimatic variables derived from the WorldClim database (Nix, 1986; Hijmans et al., 2005a,b), seven biophysical variables from the ISRIC—World Soil Information database (http://soilgrids1km.isric.org) at a resolution of 2.5 arc-min, and the occurrence information (coordinates) for each taxon as inputs (Table S3). For edaphic data we calculated a weighted mean from five depths (0–5, 5–15, 15–30, 30–60, 60–100 cm) to generate a single value for the first meter of soil for each layer, and then resampled the data from 1 to 2.5 arc min resolution to match the WorldClim dataset, using the raster package in R and ArcGIS Desktop 10.1 (Hengl et al., 2014). Distributions were further restricted by applying a taxon independent threshold, based on the Receiver Operating Characteristic (ROC) curve (Liu et al., 2005). GRIN distribution data was used to ensure that taxa distributions were not overinflated beyond known native boundaries (USDA, 2007). Soil cover data from GlobCover 2009 (Global Land Cover Map) (http://due.esrin.esa.int/page_globcover.php) further refined the maxent outputs and collecting maps by excluding urban areas, water bodies, bare areas, and permanent snow and ice regions.
Potential distribution models were considered accurate if they complied with the following conditions: (i) five-fold average area under the test ROC curve (ATAUC) is greater than 0.7, (ii) the standard deviation of ATAUC (STAUC) is less than 0.15, and (iii) At least 10% of grids for each model has standard deviation less than 0.15 (ASD15). For taxa whose Maxent model did not comply, potential distributions were estimated by forming a circular buffer of 50 km around each occurrence point for each species.
Geographic representativeness of taxa in genebank collections was calculated using the “Geographic Representativeness Score” (GRS), comparing the spatial overlap of a circular buffer surrounding each accession record (50 Km radius as described in Hijmans et al., 2001) against the potential distribution of the taxon. Ecological gaps in genebank collections were calculated using the “Ecological Representativeness Score” (ERS), calculated by comparing records to the full environmental range of the modeled taxon across ecosystem types (Olson et al., 2001). The overall priority for further collecting for ex situ conservation for each taxon was determined by averaging the SRS, GRS, and ERS with equal weight to obtain a final prioritization score (FPS), classified according to the following ranges: 1., high priority (FPS between 0 and 3); 2., medium priority (FPS between 3.01 and 5); 3., low priority (FPS between 5.01 and 7.5); and 4., and well conserved taxa (FPS between 7.51 and 10).
Expert Evaluation of Conservation Assessment Results
Predicted taxon distributions based on genebank and herbarium records were compared to the knowledge of four crop experts with experience with Helianthus distributions, systematics, conservation, and diversity. Helianthus experts were asked to evaluate of the adequacy of germplasm collections per species based on their knowledge of total accessions conserved, geographic, and environmental gaps. This assessment was given an expert priority score (EPS), analogous to the FPS score. A second score was generated, the contextual EPS, which based on additional knowledge such as in situ threats and utility to crop breeding. After initial evaluation the experts were asked to review the quantitative results, occurrence data, potential distribution models, and maps of collecting priorities. Following expert input, occurrence data were refined through elimination of incorrect points and adjustment native areas. Potential distribution modeling and gap analyses were then conducted using refined datasets to create more accurate species distribution maps. Potential zones for collecting were identified for each high priority taxon, and then combined to create maps depicting areas where multiple taxa of high priority for conservation could be collected (Figure 1).
Figure 1. Synthesis of gap analysis results and expert assessments for each of the 36 Helianthus CWR taxa surveyed. Taxa are listed by descending priority for further collecting by category: HPS, high priority taxa; MPS, medium priority taxa; LPS, low priority taxa; NFCR, no further collecting recommended. The final priority scores (FPS, black circle) is the mean of the sampling representativeness score (SRS, blue circle), geographic representativeness score (GRS, red circle), and ecological representativeness score (ERS, green circle).
Ecogeographic Niche Overlap and Phylogenetic Analyses
Potential distribution probability outputs were used when Maxent models performed well and CA50 sample buffers when Maxent models did not pass the validation criteria, to calculate niche overlap based on Schoener's D and Hellinger's I as outlined in Warren et al. (2008), and implemented in the R package Phyloclim (Heibl, 2011). Both indices utilize probability distributions in geographic space, with statistics ranging from 0 (no niche overlap) to 1 (complete niche overlap). First pairwise niche overlap was examined, then niche overlap between allopatric/sympatric taxa separately, annual/perennial taxa separately, and lastly allopatric/sympatric sister taxa. Geographic range overlap for all pairwise combinations (630 comparisons) was calculated in two ways, with respect to the larger range [(2*number of shared grid cells)/(number of grid cells in taxa A + number of grid cells in taxa B)] and with respect to the smaller range [(2*number of shared grid cells)/(Total number of grid cells in taxa A + Total number of grid cells in taxa B)]/(Total potential number of shared grid cells) [(2*total number of grid cells in species with the smaller range)/(Total number of species A + Total number of species B)].
Principal component analyses (PCA) were used to assess the importance of ecogeographic variables (Table S3) to variation in occurrence data of distribution models per taxon. A hierarchical cluster of principal components (HCPC) identified climatic clusters using R package FactoMineR (Husson et al., 2014). Boxplots for each bioclimatic and biophysical layer were created based on occurrence data points (Figure S1). Ecogeographic variables for cultivated sunflower were extracted from the area of species distribution maps (Monfreda et al., 2008) at a resolution of 5 arc-min, with a random sample of 1000 points weighted by harvested area taken from major production regions.
We downloaded the publically available 18S-26S Ribosomal DNA sequence from the external transcribed spacer (ETS) from GenBank (NCBI-http://www.ncbi.nlm.nih.gov/) for 28 of the 36 Helianthus taxa, aligned the sequences using ClustalW, and constructed a maximum likelihood phylogeny with 1000 bootstrap replications, using MEGA6 with a Jukes-Cantor nucleotide substitution model (Tamura et al., 2013). We performed a Mantel test in R utilizing the ade4 package to explore the relationship between geography and genetics (Dray and Dufour, 2007). We estimated phylogenetic signal of individual ecogeographic traits utilizing Blomberg's K (Blomberg et al., 2003), using the multiphylosignal command with 1000 permutations in Picante (Kembel et al., 2010).
Geographic Distributions of Sunflower Crop Wild Relatives
Predicted distribution maps were produced for 36 Helianthus taxa, along with taxon richness and collecting hotspot maps (Figure 2; Figure S2). Thirty of the 36 taxa (83%) produced valid maxent models with utilization of soil pH and percent sand greatly improving the accuracy of distribution models, as assessed by expert opinion (Figure 3). Five hotspots (areas of high taxon-level diversity) were identified in the USA, including the southeastern gulf coast, the south-central, the midwest, the north central, and the central east coast (Figure 2A). Our results suggest that half of the 36 taxa are in urgent need of further collecting (high priority species—HPS), along with 28% in moderate need (medium priority species—MPS), 6% of low priority (LPS), and 17% that are well represented in existing germplasm collections and thus do not require urgent additional collecting (Table 1). While the primary genepool taxa has been well collected, only 10% of the taxa in the secondary genepool are well represented across their geographic, climatic, and edaphic ranges. Likewise, only 7% of taxa in the tertiary genepool were assessed as well-conserved (Figure 1; Table 1). These results contrasted with those of expert reviewers, who classified more species as LPS. The discrepancy between the results and expert opinion was due in part to overly optimistic distribution models regarding likelihood of occurrence, in comparison to expert realities of existence of populations in these regions. Additionally, experts assessed some taxa, such as Helianthus debilis ssp. cucumerifolius, at lower priority because distributions have expanded recently as weedy populations invade new areas, and such regions were not considered by the experts as of particular priority.
Figure 2. Map of North America showing (A) taxon richness of sunflower and (B) hotspots for further collecting of high priority taxa.
Figure 3. (A) Geographic niche overlap based on ecogeographic variables, Schoener's D (above diagonal) and modified Hellinger's I (below diagonal). Taxa are grouped by phylogenetic relationship. Values range from 0 (no overlap; purple) to 1 (complete overlap; orange); (B) Occurrence points for assessed taxa grouped based on the first three principle components of biophysical and bioclimatic variables. Clusters share homogeneous bioclimatic and biophysical conditions.
Ecological Niches of Sunflower Crop Wild Relatives
Three ecogeographic clusters differentiate the taxa, with the first three PCs accounted for 74.3% of the variation (Figure 3B; Table S4). Clusters broadly corresponded to plain, desert, and woodland ecosystems (Table 1). Cluster one was mostly composed of the secondary germplasm and differentiated by temperature, while cluster two was mostly the tertiary germplasm and differentiated by precipitation. Cluster three was differentiated by soil and was evenly split between the secondary and tertiary germplasm (Table S3). It is important to note that PCA can increase type one error, so ecological niches must be carefully examined and validated (Revell, 2009; Uyeda et al., 2015). Schoener's D and Hellinger's I identified substantial niche overlap with few taxa showing niche divergence (Figure 3; Table 1).
Potential geographic distributions of crop wild relative taxa were examined for overlap with wild H. annuus (Figure S1); most (81%) taxa exhibited some geographic range overlap with H. annuus (Table 1). Among CWR taxa, 39% of pairwise comparisons had overlapping geographic distributions (sympatry), while 61% were allopatric (Table S5; Figure S3). Eight of the 12 sister taxa pairs among the CWR showed some level of sympatry (Table S6). There was considerable range asymmetry between taxa (Figure S1), with the amount of overlap depending on the direction of the comparison, where the smaller range showed 26% more overlap on average than the larger range (Table S5).
There was general niche conservatism even for sister-taxa (Figure 3; Table 2). While ecogeographic niches were fairly similar for many variables, occasionally there was substantial divergence (Figure 4; Figures S1, S4). Phylogenetic niche conservatism was found in ~54% of variables (Figure 5). Divergence was found in several soil variables suggesting an important role of soil in Helianthus diversification. A Mantel's test using Mahalanobis distance (r = 0.1423, p = 0.01), indicated that taxa that are geographically close are generally more closely related genetically. Notable exceptions to this were H. maximilliani, H. grosseserratus, and H. giganteus, which are sympatric with H. annuus, but are distantly related.
Figure 4. Climatic niches for (A) mean diurnal range and annual precipitation, (B) Soil pH and mean annual precipitation, (C) mean diurnal range and annual precipitation. Niches per taxa represent the middle 90% of occurrence points, i.e., 10% outliers are not included. Red boxes show the niche of wild H. annuus and black boxes show the niche of cultivated H. annuus in North America.
Figure 5. Test of phylogenetic signal utilizing the K for 25 of 36 taxa analyzed with complete genetic and environmental information (Blomberg et al., 2003). K measures phylogenetic signal in traits, where K-values below 1 indicates low dependence of traits on evolutionary history (not conserved between taxa) and K-values above 1 indicates trait conservation over evolutionary history (traits conserved over evolutionary time). *indicates K significantly greater than 1 (p < 0.05).
There has been increased effort to digitize data related to plant species in general and CWR in particular. The public databases (GBIF, ISRIC, WorldClim, National Germplasm repositories, DivSeek) that archive these data are an increasingly important tool to conservationists, evolutionary biologists and plant breeders. Utilizing public data can reduce the research costs in terms of people hours and consumables to achieve desired environmental and food production goals. Exploring public databases can provide a targeted way to identify accessions for introgression that can then be used to validate predicted extreme variation. This may be a way to more quickly utilize germplasm collections and provide a link to international initiatives aimed at facilitating more use of plant genetic resources (www.DivSeek.org). Here we have used geographic occurrence, bioclimatic, and biophysical data to predict species distributions, range overlap, and niche occupancy in 36 Helianthus taxa that are cross-compatible with cultivated sunflower and thus likely to be useful in crop breeding. As discussed briefly below, our results not only have implications for conservation genetics and breeding in Helianthus, but they also impact our understanding of the role of geography in the origin of species in this group.
Implications for Conservation and Plant Breeding
Our approach is both new and complementary to previous work on Helianthus species distributions and CWR in the literature (Thompson et al., 1981; Rogers et al., 1982). The method of constraining ranges to known native distributions may have limited our identification of some the extreme variation. Despite this, many taxa that diverge ecologically from cultivated sunflower were identified (Figure 4; Table 1). It was also possible to identify extreme populations within taxa that showed potential adaptation to different ecological niches.
Taxa with larger ranges tend to have greater resilience to changes in environmental conditions than taxa with more limited distributions (Sexton et al., 2014; Sheth and Angert, 2014). Thus, the latter may be considered a primary priority for conservation. Several taxa have expanded far beyond their historical ranges, including H. annuus, H. petiolaris Nutt., H. argophyllus Torrey and Gray, H. giganteus L. and H. tuberosus L. While taxa from the non-native parts of their ranges have not been prioritized, existing accessions from such ranges are acknowledged, and may be worthwhile for exploration for traits useful in crop breeding.
Clustering of CWR by environmental variables has great utility by allowing genetic resources to be exploited in a more targeted manner. For example, with respect to soil pH the taxa H. atrorubens, H. resinosus, and H. deserticola occupy different ecological space from cultivated H. annuus (Figure 4). These taxa represent potential candidates for tolerance to acid or alkaline soils, particularly to improve the ability of the crop to accumulate heavy metals for phytoremediation (Fassler et al., 2010). Surprisingly, when examining the properties of the primary, secondary, and tertiary germplasm, often extreme profiles are found in the primary germplasm. This is fortuitous since introgression from primary germplasm is more likely to be successful (Figure 4; Figure S1; Table S7). Approximately 650 wild H. annuus accessions are conserved in genebanks which occur outside the ecological parameters of the cultivar (Table S7). The general reduction of environmental diversity occupied by the cultivated sunflower relative to wild H. annuus may indicate the reduction in genetic diversity occurring through domestication.
Recent advances in plant and animal breeding (e.g., marker assisted selection, genomic selection) have been facilitated by low cost molecular marker technologies resulting in new tools that can be used to broaden the genetic base in crops (Tester and Langridge, 2010). These methods can shorten breeding cycles, increasing genetic gain per unit time, and allow for wider crosses to be utilized by minimizing linkage drag (Bernardo, 2008). The recent development of genome wide marker sets (Bowers et al., 2012; Renaut et al., 2013) and release of the H. annuus genome (Kane et al., 2011; http://www.sunflowergenome.org) facilitate the use of marker assisted selection (Iftekharuddaula et al., 2011) by decreasing costs and increasing data resolution. Further, if germplasm collections are genotyped, these data can be used to associate particular allelic variants with environmental adaptation (Fang et al., 2014).
Range Overlap of Wild Relatives of Sunflower
Sister species in Helianthus often have overlapping ranges, an observation that is consistent with sympatric and “budding” speciation (parapatric or peripheral range speciation). Substantial range asymmetry among some (but not all) sister species is also consistent with a budding speciation scenario (Table S6). The amount of range overlap between sister taxa in Helianthus is similar to recent reports from other plant genera, but different from many animal groups, where allopatry tends to be the rule in speciation (Mayr, 1954; Soltis et al., 2004; Quenouille et al., 2011; Anacker and Strauss, 2014). This may suggest that geographic isolation is less critical to plant than animal speciation, perhaps because of the low vagility of many plant species.
Unlike sympatric congeners in other plant groups (Anacker and Strauss, 2014; Grossenbacher et al., 2014), Helianthus sister taxa typically lack strong ecological divergence. This observation is inconsistent with most models of speciation involving gene flow, which assume divergent ecological selection (Via, 2009). Possibly, our analyses lacked sufficient resolution or focus on key ecological attributes to detect real differences between the ecological niches of these species. For example, it is possible that there has been pollinator and phenological divergence between sister species that was not included in our analyses. Alternatively, local niche differences between sympatric populations may have been masked by substantial ecological heterogeneity among populations of the more widely ranging species. Additionally, the approach used was designed to analyze potential habitat in the historical, native range, rather than recent range expansions, which in many cases may be recent introductions facilitated by humans, perhaps accounting for observations of limited ecological divergence.
Our analyses imply that many Helianthus taxa have similar ecological niches and exhibit niche conservatism. Under niche conservatism, greater allopatric and parapatric speciation is predicted, as habitat fragmentation is expected to contribute to reproductive isolation (Loera et al., 2012). While such a speciation strategy would be surprising given the overlap in geographic range of sister species within Helianthus, this trend has been observed in North American Ephedra (Loera et al., 2012). That larger amount of niche conservatism observed here than in other systems may be due to properties of the K-statistic, which can have inflated values in polyphyletic phylogenies and in the presence of incomplete lineage sorting, both of which occur in Helianthus (Rosenthal et al., 2002; Gross and Rieseberg, 2005; Horandl and Stuessy, 2010; Davies et al., 2012).
Using a combination of gap analysis, environmental niche modeling, and phylogenetic approaches 36 CWR of sunflower were examined. Taxa that are under-represented in germplasm collections as well as species and populations inhabiting environmental niches with extreme phenotypes that may possess traits of value to crop improvement were identified. In Helianthus, sister taxa appear to occur more frequently in sympatry than allopatry, possibly suggesting that speciation may occur in the presence of gene flow. Finally, much of the primary genepool occurs in extreme environments indicating that utilization of wild H. annuus for the breeding of abiotic stress tolerance may produce quick gains with minimal effort.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was undertaken as part of the initiative “Adapting Agriculture to Climate Change: Conserving, Protecting and Preparing Crop Wild Relatives” which is supported by the Government of Norway. The project is managed by the Global Crop Diversity Trust with the Millennium Seed Bank of the Royal Botanic Gardens, Kew UK and implemented in partnership with national and international genebanks and plant breeding institutes around the world. For further information, go to the project website: http://www.cwrdiversity.org/. Funding was provided by the aforementioned initiative, The National Sunflower Association, The U.S. National Science Foundation, The Natural Sciences and Engineering Research Council of Canada, Genome BC, and Genome Canada.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.00841
Anacker, B. L., and Strauss, S. Y. (2014). The geography and ecology of plant speciation: range overlap and niche divergence in sister species. Proc. R. Soc. B 281:20132980. doi: 10.1098/rspb.2013.2980
Blackman, B. K., Scascitelli, M., Kane, N. C., Luton, H. H., Rasmussen, D. A., Bye, R. A., et al. (2011). Sunflower domestication alleles support single domestication center in eastern North America. Proc. Natl. Acad. Sci. U.S.A. 108, 14360–14365. doi: 10.1073/pnas.1104853108
Blomberg, S. P., Garland, T., and Ives, A. R. (2003). Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution 57, 717–745. doi: 10.1111/j.0014-3820.2003.tb00285.x
Bowers, J. E., Nambeesan, S., Corbi, J., Barker, M. S., Rieseberg, L. H., Knapp, S. J., et al. (2012). Development of an ultra-dense genetic map of the sunflower genome based on single-feature polymorphisms. PLoS ONE 7:e51360. doi: 10.1371/journal.pone.0051360
Cabrera-Bosquet, L., Crossa, J., von Zitzewitz, J., Serret, M. D., and Araus, J. L. (2012). High−throughput phenotyping and genomic selection: the frontiers of crop breeding converge. J. Integr. Plant Biol. 54, 312–320. doi: 10.1111/j.1744-7909.2012.01116.x
Challinor, A. J., Watson, J., Lobell, D. B., Howden, S. M., Smith, D. R., and Chhetri, N. (2014). A meta-analysis of crop yield under climate change and adaptation. Nat. Clim. Change 4, 287–291. doi: 10.1038/nclimate2153
Davies, T. J., Kraft, N. J. B., Salamin, N., and Wolkovich, E. M. (2012). Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism. Ecology 93, 242–247. doi: 10.1890/11-1360.1
Fang, Z., Gonzales, A. M., Clegg, M. T., Smith, K. P., Muehlbauer, G. J., Steffenson, B. J., et al. (2014). Two genomic regions contribute disproportionately to geographic differentiation in wild barley. G3 4, 1193–1203. doi: 10.1534/g3.114.010561
FAOSTAT. (2013). Final Data 2013. Available online at: http://faostat.fao.org (Accessed May, 2015).
Fassler, E., Robinson, B. H., Stauffer, W., Gupta, S. K., Papritz, A., and Schulin, R. (2010). Phytomanagement of metal contaminated agricultural land using sunflower, maize and tobacco. Agric. Ecosyst. Environ. 136, 49–58. doi: 10.1016/j.agee.2009.11.007
Grossenbacher, D. L., Veloz, S. D., and Sexton, J. P. (2014). Niche and range size patterns suggest that speciation begins in small, ecologically diverged populations in North American monkeyflowers (Mimulus spp.). Evolution 68, 1270–1280. doi: 10.1111/evo.12355
Harter, A. V., Gardner, K. A., Falush, D., Lentz, D. L., Bye, R. A., and Rieseberg, L. H. (2004). Origin of extant domesticated sunflowers in eastern North America. Nature 430, 201–205. doi: 10.1038/nature02710
Heibl, C. (2011). Webcite Phyloclim: Integrating Phylogenetics and Climatic Niche Modelling. Available online at: http://cran.r-project.org/web/packages/phyloclim/index.html
Hengl, T., de Jesus, J. M., MacMillan, R. A., Batjes, N. H., Heuvelink, G. B. M., Ribeiro, E., et al. (2014). SoilGrids1km—global soil information based on automated mapping. PLoS ONE 9:e105992. doi: 10.1371/journal.pone.0105992
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., and Jarvis, A. (2005b). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978. doi: 10.1002/joc.1276
Hijmans, R. J., Guarino, L., Jarvis, A., O'Brien, R., Mathur, P., Bussink, C., et al. (2005a). DIVA-GIS Version 5.2 Manual. Available online at: http://www.diva-gis.org/
Hoisington, D., Khairallah, M., Reeves, T., Ribault, J. M., Skovmand, B., Taba, S., et al. (1999). Plant genetic resources: what can they contribute toward increased crop productivity? Proc. Natl. Acad. Sci. U.S.A. 96, 5937–5943.
Husson, F., Josse, J., Le, S., and Mazet, J. (2014). FactoMineR: Multivariate Exploratory Data Analysis and Data Mining with R. R package version 1.27. Available online at: http://CRAN.R-project.org/package=FactoMineR
Iftekharuddaula, K. M., Newaz, M. A., Salam, M. A., Ahmed, H. U., Mahbub, M. A. A., Septiningsih, E. M., et al. (2011). Rapid and high-precision marker assisted backcrossing to introgress the SUB1 QTL into BR11, the rainfed lowland rice mega variety of Bangladesh. Euphytica 178, 83–97. doi: 10.1007/s10681-010-0272-2
Kembel, S. W., Cowan, P. D., Helmus, M. R., Cornwell, W. K., Morlon, H., Ackerly, D. D., et al. (2010). Picante: R tools for integrating phylogenies and ecology. Bioinformatics 26, 1463–1464. doi: 10.1093/bioinformatics/btq166
Khoury, C., Laliberté, B., and Guarino, L. (2010). Trends in ex situ conservation of plant genetic resources: a review of global crop and regional conservation strategies. Genet. Resour. Crop Evol. 57, 625–639. doi: 10.1007/s10722-010-9534-z
Kolkman, J. M., Berry, S. T., Leon, A. J., Slabaugh, M. B., Tang, S., Gao, W., et al. (2007). Single nucleotide polymorphisms and linkage disequilibrium in sunflower. Genetics 177, 457–468. doi: 10.1534/genetics.107.074054
Liu, C., Berry, P. M., Dawson, T. P., and Pearson, R. G. (2005). Selecting thresholds of occurrence in the prediction of species distributions. Ecography 28, 385–393. doi: 10.1111/j.0906-7590.2005.03957.x
Loera, I., Sosa, V., and Ickert-Bond, S. M. (2012). Diversification in North American arid lands: Niche conservatism, divergence and expansion of habitat explain speciation in the genus Ephedra. Mol. Phylogenet. Evol. 65, 437–450. doi: 10.1016/j.ympev.2012.06.025
Mandel, J. R., Dechaine, J. M., Marek, L. F., and Burke, J. M. (2011). Genetic diversity and population structure in cultivated sunflower and comparison to its wild progenitor Helianthus annuus L. Theor. Appl. Genet. 123, 693–704. doi: 10.1007/s00122-011-1619-3
Monfreda, C., Ramankutty, N., and Foley, J. A. (2008). “Farming the planet: 2. Geographic distribution of crop areas, yields, physiological types, and net primary production in the year 2000,” in Global Biogeochemical Cycles 22: GB1022. Available online at: http://www.geog.mcgill.ca/landuse/pub/Data/175crops2000/
Olson, D. M., Dinerstein, E., Wikramanayake, E. D., Burgess, N. D., Powell, G. V. N., Underwood, E. C., et al. (2001). Terrestrial ecoregions of the world: a new map of life on earth. BioScience 51, 933–938. doi: 10.1641/0006-3568(2001)051[0933:TEOTWA]2.0.CO;2
Phillips, S. J. (2008). Transferability, sample selection bias and background data in presence-only modeling: a response to Peterson et al. (2007). Ecography 31, 272–278. doi: 10.1111/j.0906-7590.2008.5378.x
Porter, J. R., Xie, L., Challinor, A. J., Cochrane, K., Howden, S. M., Iqbal, M. M., et al. (2014). “Food security and food production systems,” in Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds C. B. Field, V. R. Barros, D. J. Dokken, K. J. Mach, M. D. Mastrandrea, T. E. Bilir, M. Chatterjee, K. L. Ebi, Y. O. Estrada, R. C. Genova, B. Girma, E. S. Kissel, A. N. Levy, S. MacCracken, and P. R. Mastrandrea (Cambridge; New York, NY: Cambridge University Press), 485–533.
Quenouille, B., Hubert, N., Bermingham, E., and Planes, S. (2011). Speciation in tropical seas: allopatry followed by range change. Mol. Phylogenet. Evol. 58, 546–552. doi: 10.1016/j.ympev.2010.12.009
Ramírez-Villegas, J., Khoury, C., Jarvis, A., Debouck, D. G., and Guarino, L. (2010). A gap analysis methodology for collecting crop genepools: a case study with Phaseolus beans. PLoS ONE 5:e13497. doi: 10.1371/journal.pone.0013497
Renaut, S., Grassa, C. J., Yeaman, S., Moyers, B. T., Lai, Z., Kane, N. C., et al. (2013). Genomic islands of divergence are not affected by geography of speciation in sunflowers. Nat. Commun. 4:1827. doi: 10.1038/ncomms2833
Rosenthal, D. M., Schwarzbach, A. E., Donovan, L. A., Raymond, O., and Rieseberg, L. H. (2002). Phenotypic differentiation between three ancient hybrid taxa and their parental species. Int. J. Plant Sci. 163, 387–398. doi: 10.1086/339237
Sheth, S. N., and Angert, A. L. (2014). The evolution of environmental tolerance and range size: a comparison of geographically restricted and widespread Mimulus. Evolution 68, 2917–2931. doi: 10.1111/evo.12494
Stebbins, J. C., Winchell, C. J., and Constable, J. V. H. (2013). Helianthus winteri (Asteraceae), a new perennial species from the southern Sierra Nevada foothills, California. Aliso 31, 19–24. doi: 10.5642/aliso.20133101.04
Tang, S., and Knapp, S. J. (2003). Microsatellites uncover extraordinary diversity in native American land races and wild populations of cultivated sunflower. Theor. App. Genet. 106, 990–1003. doi: 10.1007/s00122-002-1127-6
USDA (2012). USDA, ARS, National Genetic Resources Program, Germplasm Resources Information Network−(GRIN). [Online Database]. Beltsville, MD: National Germplasm Resources Laboratory. Available online at: http://www.ars-grin.gov/cgi-bin/npgs/acc/display.pl?1346916
VanDerWal, J., Shoo, L. P., Graham, C. H., and Williams, S. E. (2009). Selecting pseudo-absence data for presence-only distribution modeling: how far should you stray from what you know? Ecol. Model. 220, 589–594. doi: 10.1016/j.ecolmodel.2008.11.010
Vincent, H., Wiersema, J., Kell, S., Fielder, H., Dobbie, S., Castañeda-Álvarez, N. P., et al. (2013). A prioritized crop wild relative inventory to help underpin global food security. Biol. Conserv. 167, 265–275. doi: 10.1016/j.biocon.2013.08.011
Keywords: conservation, climate change, crop wild relatives, ecological niche modeling, plant breeding, plant genetic resources, publicly available data sources
Citation: Kantar MB, Sosa CC, Khoury CK, Castañeda-Álvarez NP, Achicanoy HA, Bernau V, Kane NC, Marek L, Seiler G and Rieseberg LH (2015) Ecogeography and utility to plant breeding of the crop wild relatives of sunflower (Helianthus annuus L.). Front. Plant Sci. 6:841. doi: 10.3389/fpls.2015.00841
Received: 01 August 2015; Accepted: 24 September 2015;
Published: 08 October 2015.
Edited by:Jaime Prohens, Universitat Politècnica de València, Spain
Reviewed by:Inger Martinussen, Norwegian Institute of Bioeconomy Research (NIBIO), Norway
Mohammad Ehsan Dulloo, Bioversity International, Italy
Copyright © 2015 Kantar, Sosa, Khoury, Castañeda-Álvarez, Achicanoy, Bernau, Kane, Marek, Seiler and Rieseberg. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael B. Kantar, email@example.com