Skip to main content
  • Facebook
  • Twitter
  • YouTube
  • LinkedIn
  • Google Plus
  • Other GSA Resources
    • Genetics Society of America
    • Genetics
    • Genes to Genomes: The GSA Blog
    • GSA Conferences
    • GeneticsCareers.org
  • Log in
G3: Genes | Genomes | Genetics

Main menu

  • HOME
  • ISSUES
    • Current Issue
    • Early Online
    • Archive
  • ABOUT
    • About the journal
    • Why publish with us?
    • Editorial board
    • Contact us
  • SERIES
    • Genetics of Immunity
    • Genetics of Sex
    • Genomic Prediction
    • Multiparental Populations
  • ARTICLE TYPES
    • About Article Types
    • Genome Reports
    • Meeting Reports
    • Mutant Screen Reports
    • Software and Data Resources
  • PUBLISH & REVIEW
    • Scope & publication policies
    • Submission & review process
    • Article types
    • Prepare your manuscript
    • Submit your manuscript
    • After acceptance
    • Guidelines for reviewers
  • SUBSCRIBE
    • Email alerts
    • RSS feeds
  • Other GSA Resources
    • Genetics Society of America
    • Genetics
    • Genes to Genomes: The GSA Blog
    • GSA Conferences
    • GeneticsCareers.org

User menu

Search

  • Advanced search
G3: Genes | Genomes | Genetics

Advanced Search

  • HOME
  • ISSUES
    • Current Issue
    • Early Online
    • Archive
  • ABOUT
    • About the journal
    • Why publish with us?
    • Editorial board
    • Contact us
  • SERIES
    • Genetics of Immunity
    • Genetics of Sex
    • Genomic Prediction
    • Multiparental Populations
  • ARTICLE TYPES
    • About Article Types
    • Genome Reports
    • Meeting Reports
    • Mutant Screen Reports
    • Software and Data Resources
  • PUBLISH & REVIEW
    • Scope & publication policies
    • Submission & review process
    • Article types
    • Prepare your manuscript
    • Submit your manuscript
    • After acceptance
    • Guidelines for reviewers
  • SUBSCRIBE
    • Email alerts
    • RSS feeds
Previous ArticleNext Article

Genomic Signatures of Adaptation to a Precipitation Gradient in Nigerian Sorghum

Marcus O. Olatoye, Zhenbin Hu, Fanna Maina and View ORCID ProfileGeoffrey P. Morris
G3: Genes, Genomes, Genetics October 1, 2018 vol. 8 no. 10 3269-3281; https://doi.org/10.1534/g3.118.200551
Marcus O. Olatoye
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhenbin Hu
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fanna Maina
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Geoffrey P. Morris
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Geoffrey P. Morris
  • For correspondence: gpmorris@ksu.edu
  • Article
  • Figures & Data
  • Supplemental
  • Info & Metrics
Loading

Abstract

Evolution of plants under climatic gradients may lead to clinal adaptation. Understanding the genomic basis of clinal adaptation in crops species could facilitate breeding for climate resilience. We investigated signatures of clinal adaptation in the cereal crop sorghum (Sorghum bicolor L. [Moench]) to the precipitation gradient in West Africa using a panel (n = 607) of sorghum accessions from diverse agroclimatic zones of Nigeria. Significant correlations were observed between common-garden phenotypes of three putative climate-adaptive traits (flowering time, plant height, and panicle length) and climatic variables. The panel was characterized at >400,000 single nucleotide polymorphisms (SNPs) using genotyping-by-sequencing (GBS). Redundancy analysis indicated that a small proportion of SNP variation can be explained by climate (1%), space (1%), and climate collinear with space (3%). Discriminant analysis of principal components identified three genetic groups that are distributed differently along the precipitation gradient. Genome-wide association studies were conducted with phenotypes and three climatic variables (annual mean precipitation, precipitation in the driest quarter, and annual mean temperature). There was no overall enrichment of associations near a priori candidate genes implicated in flowering time, height, and inflorescence architecture in cereals, but several significant associations were found near a priori candidates including photoperiodic flowering regulators SbCN12 and Ma6. Together, the findings suggest that a small (3%) but significant proportion of nucleotide variation in Nigerian sorghum landraces reflects clinal adaptation along the West African precipitation gradient.

  • Crop evolution
  • climate adaptation
  • flowering time
  • West Africa
  • genome-wide association studies

Adaptation to environmental gradients can lead to evolution of clines, a continuous form of local adaptation (Slatkin 1973; Savolainen et al. 2013; Yoder et al. 2014). Climatic gradients that may shape trait adaptation include ultraviolet radiation, photoperiodicity, temperature, and precipitation (Hancock et al. 2011; Haussmann et al. 2012; Bastide et al. 2016). In the model plant Arabidopsis, clinal variation across latitude has been observed for seed dormancy, cold tolerance, height, and flowering time (Zhen and Ungerer 2008; Samis et al. 2012; Kronholm et al. 2012; Debieu et al. 2013). Likewise, in crop species, diffusion from tropical to temperate zones has led to clinal adaptation in flowering time (Ducrocq et al. 2008; Buckler et al. 2009; Wu et al. 2013; Kloosterman et al. 2013) and cold tolerance (Comadran et al. 2012; Ma et al. 2015). Local adaptation of traditional varieties has played a key role in smallholder crop production under adverse climatic conditions and low agricultural inputs (Vasconcelos et al. 2013). Locally-adapted crop landraces possess alleles that can be beneficial for the development of improved varieties to ensure food security under stressful climates (Zeven 1998; Soler et al. 2013; Lasky et al. 2015).

Understanding genetic diversity, population structure, and genotype-phenotype associations in crop landraces can guide germplasm conservation and breeding (Djè et al. 2000; Soler et al. 2013; Dwivedi et al. 2016). Recent advances in genotyping technology have facilitated studies of genomic diversity in crops, including studies of local and clinal adaptation using population and quantitative trait genomics (Myles et al. 2009; Morrell et al. 2012; Savolainen et al. 2013). Population genomics methods based on genome-wide patterns of nucleotide variation can identify loci with signatures of selection (Siol et al. 2010). These have been used to identify genomic targets of adaptation in many crops including maize (Gore et al. 2009; Hufford et al. 2012), rice (Meyer et al. 2016; Li et al. 2017), and sorghum (Morris et al. 2013; Mace et al. 2013). In sorghum, quantitative trait genomics approaches using mixed linear model have identified genomic regions associated with adaptive traits and climatic variables (Morris et al. 2013; Zhang et al. 2015; Lasky et al. 2015).

Sorghum (Sorghum bicolor L. [Moench]) is an essential staple cereal crop in dryland regions of the world (National Research Council 1996). It has adapted to a wide variety of climatic gradients and has abundant phenotypic variation for flowering time, plant morphology, and inflorescence morphology (Morris et al. 2013; Zhang et al. 2015; Lasky et al. 2015). Globally, the morphological types (botanical races) of sorghum are distributed according to precipitation zones, with open-panicle guinea types predominant in humid regions, semi-compact caudatum types predominant in semi-arid regions, and compact-panicle durra types predominant in arid regions (Doggett 1988; Morris et al. 2013). In West Africa, sorghum is found across a steep north-south precipitation gradient, ranging from semiarid grasslands bordering the Sahara Desert in the north (Sahelian zone), through subhumid savannah (Sudanian zone), to humid forest zones in the south (Guinean zone). These regions have been subject to major droughts for several millennia (Shanahan et al. 2009) and increased drought under climate change is expected to reduce sorghum yields in this region (Lobell et al. 2008).

The West African country of Nigeria is Africa’s most populous nation and its largest sorghum producer, with 5-10 million Mg of grain production per year (Nzeka and Akhidenor 2018). Sorghum is the major cereal in the northern Sudano-Sahelian region of Nigeria, which is characterized by prolonged dry seasons and short rainy seasons (National Research Council 1996). Sorghum, as a non-centric crop, has multiple centers of diversity and two of these overlap with the boundaries of Nigeria (Harlan 1971, 1992). The genetic diversity of Nigerian sorghum is poorly characterized compared to other African sorghum germplasm (Rao et al. 1985; Deu et al. 2008; Barro-Kondombo et al. 2010; Leiser et al. 2014). Identifying genomic regions underlying adaptation in Nigerian sorghum germplasm could facilitate the identification of adaptive traits and genetic diversity relevant to crop improvement.

Given that sorghum is distributed across the precipitation gradient in Nigeria, we hypothesized that Nigerian sorghum germplasm has been shaped by clinal adaptation. Under this hypothesis, we expect precipitation variables to be associated with both phenotype (putative climate-adaptive traits) and genotype (population structure and SNPs). Further, we expect that trait-associated and climate-associated genome regions will colocalize with genes involved in putative climate-adaptive traits. We investigated these predictions in a large panel of georeferenced Nigerian genebank accessions, which were previously phenotyped and which we genotyped at high-density using GBS. We characterized patterns of association among climatic, phenotypic, and genotypic variables and tested colocalization of associated genomic regions. Overall, the patterns are consistent with a small contribution of clinal adaptation shaping genomic variation in Nigerian sorghum.

Materials and Methods

Plant materials

Seeds for 553 Nigerian accessions were obtained from the USDA National Plant Germplasm System (NPGS) (https://www.ars-grin.gov/). Seedlings were raised in a greenhouse for two weeks and 50 mg of fresh leaf tissue was collected from each accession into 96-well plates. A control well was left empty on each plate. Leaf tissue was lyophilized (Labconco Freeze Dryer, Kansas City, MO, USA) for two days and then ground using 96-well plate plant tissue grinder (Retsch Mixer Mill, Haan, Germany). Genomic DNA was extracted using BioSprint 96 DNA Plant Kit (QIAGEN, Valencia CA, USA), quantified using Quant-iTTM PicoGreen dsDNA Assay Kit (ThermoFisher Scientific, Waltham MA, USA), then normalized to 10 ng/μl.

Genotyping-by-sequencing

GBS was conducted on 553 Nigerian accessions using methods previously described (Elshire et al. 2011; Morris et al. 2013). Briefly, individual DNA samples were digested using ApeKI restriction enzyme (NEB R0643L) followed by ligation of barcode and common adapters ligation using T4 DNA ligase (NEB M0202L). Ligated libraries were pooled (96-plex libraries) then amplified by polymerase chain reaction (PCR). Purification of libraries was performed using QIAquick PCR purification kit (QIAGEN, Valencia CA, USA). Library size distribution was obtained using a Bioanalyzer (Agilent Technologies 2100, Santa Clara CA, USA). Four 96-plex libraries were pooled to generate 384-plex sequencing libraries. Libraries were sequenced using single end 100-cycle sequencing using Illumina HiSeq2500 (Illumina, San Diego CA, USA) at the University of Kansas Medical Center, Kansas City MO, USA.

Sequence reads for Nigerian germplasm were combined with published sequence reads obtained for 1943 accessions (Lasky et al. 2015). The published sequence data were composed of globally diverse sorghum landraces with major representation by accessions from Africa and Asia. They were obtained from the United States NPGS-GRIN and the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) gene banks. From these 1943 georeferenced global accessions, sequence information from 158 Nigerian accessions was obtained and combined with the Nigerian NPGS set. Duplicated accessions and accessions with sorghum conversion (SC) numbers (i.e., with introgressions for early maturity and semi-dwarf genes) in the NPGS database were removed from the Nigerian germplasm. Thus, 607 Nigerian accessions (of which 443 were georeferenced; Figure 1A) and 1785 georeferenced global accessions were used for downstream analysis (Files S1 and S2). Reads were aligned to the sorghum reference genome v3.0 (McCormick et al. 2018) using Burrow Wheeler Alignment algorithm (Li and Durbin 2009). SNP calling was performed using TASSEL 5.0 GBS pipeline (Glaubitz et al. 2014). The SNPs were filtered for < 20% missingness, then missing data were imputed using BEAGLE 4.0 (Browning and Browning 2013).

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Geographic distribution of Nigerian accessions across precipitation zones. (A) Geographic distribution of georeferenced Nigerian accessions with known botanical race. (B) Panicle length variation across precipitation zones (size of circles commensurate to length of panicle in cm). (C) plant height variation across precipitation zones (size of circles commensurate to height of plant in cm). (D) flowering time variation across precipitation zones (size of circles commensurate to days until anthesis, 1 = earliest flowering and 6 = photoperiod sensitive). Background map is color-coded according to average annual precipitation (in mm).

Climate and phenotype data

Climate data (average from 1960 to 1990) were obtained from WorldClim 1.4 using the Raster package in R (Hijmans et al. 2016) based on the coordinate (latitude and longitude) for each of the 443 georeferenced Nigerian accessions (File S1) and 1785 global accessions (File S2). As proxies for precipitation gradients that are hypothesized to affect Nigerian sorghum we investigated “annual precipitation” and “precipitation in the driest quarter”. Common garden passport data for flowering time, plant height, and panicle length for the Nigerian accessions were obtained from the USDA-NPGS Germplasm Resource Information Network database (https://www.ars-grin.gov/). The passport data were based on evaluations in one or more common garden experiments in tropical latitudes (Puerto Rico and St. Croix, 17-18°N), so best linear unbiased predictors (BLUPs) were estimated for each trait for each accession (File S1). A term for common garden was fit as random in the BLUP estimation model using lmer function in LME4 package in R (Bates et al. 2014) as follows:Embedded Imagewhere Embedded Image is the vector of phenotypic observation of the ith accession, Embedded Image is the jth common garden where ith accession was evaluated, and Embedded Image is the residual or error term. Pearson correlations were calculated between BLUPs of three adaptive traits (flowering time, panicle length, and plant height), and environmental factors (latitude, temperature and precipitation). To reduce the influence of outlier sites with exceptional climatic variables, precipitation values of three geographical locations where sorghum is not commonly cultivated were removed. Analysis of variance (ANOVA) and Tukey HSD test in R were performed to identify precipitation differences among sites of origin for different genetic groups or botanical races.

Redundancy analysis

Redundancy analysis (RDA) was performed separately for global and Nigerian germplasm sets using the varpart function in the R vegan package (Oksanen et al. 2017). A multivariate model was fit using the genomic data (431,698 SNPs for the global accession and 279,689 SNPs for the Nigerian accessions, filtered for monomorphic and singleton markers) as response variable. Ten WorldClim 1.4 climatic variables (annual mean temperature, mean temperature wettest quarter, mean temperature driest quarter, mean temperature warmest quarter, mean temperature coldest quarter, annual precipitation, precipitation wettest quarter, precipitation driest quarter, precipitation in the warmest quarter, and precipitation in the coldest quarter) and geographical variables (latitude and longitude, which we refer to as “space”) were fitted as predictor terms. The “space” term is included to account for isolation-by-distance (Lasky et al. 2012). To test the significance of the proportion of variation explained by climate collinear with space in the Nigerian germplasm, the proportion of variation explained was compared to the distribution from 1000 permuted data sets. In each stage of the permutation, individuals (genotypes) were randomized and RDA regression fitted and repeated 1000 times.

Population structure and linkage disequilibrium analyses

Discriminant analysis of principal components (DAPC) was conducted with the find clusters function in Adegenet package in R (Jombart 2008; Jombart et al. 2010). Population differentiation (FST) between DAPC groups was estimated using –weir-fst-pop parameter (Weir and Cockerham’s FST) in VCFtools (Danecek et al. 2011). While nucleotide diversity within DAPC groups was estimated using –window-pi (1kb) in VCFtools. LD decay analysis for each DAPC group was performed by PopLDdecay (BGI-shenzhen 2017). For comparison with Nigerian germplasm, West African accessions were identified from the published global GBS data (Lasky et al. 2015). The published GBS data were composed of global accessions from 55 countries, predominantly representing landraces from sub-Saharan Africa and Asia. In the text, “global” refers to all accessions (including Nigerian and other West African accessions), unless otherwise noted that Nigerian or West African accessions have been removed.

Linkage disequilibrium decay for the genomic data for Nigerian, West African, and global germplasm was estimated by PopLDdecay (BGI-shenzhen 2017), with minor allele frequency parameter set at 0.05 and smoothing by the spline function in R. Principal component analysis (PCA) was performed using SNPrelate package in R (Zheng et al. 2012) with LD pruning threshold parameter set to 0.5 and minor allele frequency parameter set to 0.05. Neighbor-joining analysis was performed using TASSEL 5.0 and visualized in APE (Analyses of Phylogenetics and Evolution) package in R (Paradis et al. 2004). Population differentiation (FST) between Nigerian, West African and global germplasm was evaluated using –weir-fst-pop parameter (Weir and Cockerham’s FST) in VCFtools (Danecek et al. 2011). While nucleotide diversity within each germplasm, their inbreeding coefficients, and observed heterozygosity were estimated using –window-pi (1kb window), –het, and –hardy respectively, in VCFtools.

Genome-wide association studies

Genome-wide association studies (GWAS) were performed using BLUPs of traits (panicle length, n = 330; plant height, n = 332; and flowering time, n = 412). After filtering the Nigerian data for a minimum minor allele frequency (MAF) of 0.03, a total of 149,342 SNPs were used in the GWAS analysis. First, a multi-locus mixed linear model (MLMM) (Segura et al. 2012) with a fixed population term (Q) and a random polygenic term (K) was used to perform GWAS for the phenotypic traits. PCA components (first three PCs) used for Q term were estimated using TASSEL 5.0 (Bradbury et al. 2007) and kinship matrix used for the polygenic term was derived from GAPIT (Lipka et al. 2012). Bonferroni correction of 2.6e-07 (α/number of markers, where α = 0.05) was used to determine the cut-off threshold for the phenotypic associations. A set of a priori candidate genes was compiled from Phytozome including known sorghum genes, and sorghum homologs of rice and maize genes known to be involved in inflorescence morphology, maturity, and plant height (n = 169; File S3).

Genome scans

Three environmental variables (annual precipitation, precipitation in the driest quarter, and annual mean temperature) were used as proxies for the precipitation gradient (n = 443). A GLM, which does not include population structure and kinship terms, was used to perform an association scan for climatic variables to reduce false negatives (Bergelson and Roux 2010; Lasky et al. 2015). The top 1% outliers of the environmental associations were selected for enrichment analysis. Enrichment of a priori genes near association peak were performed using a chi-square test. Windows of 100 kb were used as conservative regions for colocalization between SNPs and a priori genes since LD decayed to background levels at > 100 kb. The genome wide Tajima’s D across 100 kb windows was tested using VCF tools in global germplasm (DGlobal), West African germplasm (DWestAfrica), and Nigerian germplasm (DNigeria). Enrichment analysis for a priori genes was performed by testing whether the DNigeria 100 kb windows were significantly enriched for our a priori candidate genes relative to a set of random genes derived from the sorghum genome version 3 gff3 gene file from Phytozome (Goodstein et al. 2011; McCormick et al. 2018) for 1000 whole genome permutations.

Data availability

Raw sequencing data are available from the NCBI Sequence Read Archive under project accession SRP132525 SNP genotype, phenotype, and geographic data are available at Dryad (doi:10.5061/dryad.g0141g7). All data are publicly available. File S1 contains detailed descriptions of Nigerian accessions, their passport data, georeference information, the BLUPs of phenotypes, climatic data, and DAPC groups. File S2 contains detailed descriptions of global accessions, their georeference information, and climatic data. File S3 contains a priori candidate genes list and literature sources. File S4 contains ANOVA and Tukey test results for race by precipitation analysis. File S5 contains detailed descriptions of a priori candidate genes associated with significant SNPs for MLMM and GLM GWAS results for the phenotypes. File S6 contains detailed descriptions of a priori candidate genes associated with outlier SNPs for GLM of environmental variables. File S7 contains detailed descriptions of a priori candidate genes associated with Nigerian germplasm Tajima’s D (DNigeria) windows. Supplemental material available at Figshare: https://doi.org/10.25387/g3.6942986.

Results

Trait and environment correlations

The georeferenced sorghum accessions from Nigeria originated across a wide precipitation gradient (Figure 1A). Annual precipitation at the locations of origin ranges from < 500 mm/year for the northern Nigerian accessions (Sahelian zone) to > 3500 mm/year for the southern Nigerian accessions (Guinean zone). Analysis of variance indicates a significant difference in annual precipitation at sites of origin among the botanical races (P-value < 0.001). Significant differences were found between guinea and caudatum types (P-value < 0.05) and between guinea and durra-caudatum types (P-value < 0.01) (File S4) in precipitation differences among sites of origin.

The correlation of annual precipitation with common garden phenotypes of georeferenced Nigerian germplasm was investigated for three traits (panicle length, plant height, and flowering time) (Figure 1B-D). For comparison, we also considered correlations with annual mean temperature and latitude. Significant positive relationships with annual mean precipitation were observed for panicle length (Figure 2A; r = 0.21, P-value < 0.001) and plant height (Figure 2B; r = 0.22, P-value < 0.001) but not flowering time (Figure 2C). By contrast, annual mean temperature had no correlation with panicle length (Figure 2D) or plant height (Figure 2F) but a significant negative relationship with flowering time (Figure 2F; r = - 0.19; P-value < 0.001). Latitude had negative relationships with panicle length (Figure 2G; r = - 0.19, P-value < 0.001) and plant height (Figure 2H; r = - 0.26, P-value < 0.001), but no relationship with flowering time (Figure 2I). Among the traits, flowering time had significant positive relationships with both panicle length and plant height (Fig. S1; r = 0.32 and 0.41, respectively; P-values < 0.001). Among the environmental variables, annual mean precipitation had a significant negative relationship with latitude (Fig. S1; r = -0.86, P-value < 0.001) and a weaker but significant negative correlation with annual mean temperature (Fig. S1; r = -0.23, P-value < 0.001).

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2

Relationships between traits and environmental variables in georeferenced germplasm. Pearson correlations between three environmental variables, annual mean precipitation (in mm) (A-C), annual temperature (in °C) (D-F), and latitude (G-I), vs. three phenotypic variables, panicle length (A, D, G), plant height (B, E, H), and flowering time (C, F, I) (n = 443). Phenotypes are BLUPs from common garden experiments. Significant correlation values are noted at 0.05, 0.01 and 0.001 (*, **, and ***).

Genome-wide nucleotide variation

To investigate genomic variation we developed a data set consisting of 431,698 SNPs genotyped across 2392 accessions (Nigeria, West Africa, and global). Most of the SNP variation in the three panels was rare; about 51% of the Nigerian genomic data were composed of SNPs with minor allele frequencies (MAF) < 0.01, 46% of the West African SNPs have MAF < 0.01, 36% of the global reference SNPs have MAF < 0.01, and 37% of global SNPs have MAF < 0.01 (Fig. S2A). The mean observed heterozygosity across loci in each of the germplasm is 0.02 (2%) (Fig. S2B). SNP density was higher in sub-telomeric regions and lower in sub-centromeric regions (Fig. S3). In the Nigerian germplasm (n = 607 accessions; Figure 1A), 279,689 SNPs were retained after removing monomorphic markers, singletons, and doubletons from the initial 431,698 SNPs. This corresponds to an average of 1 SNP per 2.7 kb. In the West African germplasm (n = 325 accessions), 311,786 SNPs were obtained after removing monomorphic markers, singletons, and doubletons.

The nucleotide diversity of Nigerian germplasm was similar to the West African and global germplasm. The average nucleotide diversity for global germplasm (πglobal) and West African germplasm (πWestAfrica) was 4.5 × 10−4. The nucleotide diversity in Nigerian germplasm (πNigeria) was somewhat lower at 4.0 × 10−4. The average inbreeding coefficients were 0.83, 0.82, and 0.80 for global germplasm, West African germplasm, and Nigerian germplasm, respectively (all significantly different from each other at P-value < 0.01) (Fig. S2C). The FST between Nigerian and West African germplasm was 0.007 (Fig. S4A) and FST between Nigeria and global germplasm was 0.07 (Fig. S4B). The genome-wide average rate of linkage disequilibrium decay differed among the panels. Of the three germplasm sets, the global germplasm had the fastest LD decay rate by reducing to r2 = 0.1 at 20 kb (Fig. S5). LD decayed to half its initial value at 12 kb and r2 = 0.1 at 50 kb in the Nigerian germplasm. The West African germplasm had a slowest LD decay rate among the three sets, with r2 = 0.1 at 100 kb.

Redundancy analysis

To estimate the proportion of SNP variation that has been shaped by climate vs. geographic distance (space) we carried out a redundancy analysis, a form of multivariate regression. The proportion of SNP variation explained by climate and space in global germplasm was substantially greater than in Nigerian germplasm. In global germplasm, the proportion of variation explained by the 10 climate variables, space (latitude and longitude), and their combination together are 4%, 8%, and 5%, respectively (Figure 3A). By contrast, in Nigerian germplasm, climate and geographic variables explained a smaller proportion of the total SNP variation (Figure 3B); climate and space alone each explained 1% of the SNP variation, while climate collinear with space explained 3% of the SNP variation. The proportion of variation explained by climate collinear with space was significantly greater (P < 0.001) than the null distribution from geographically permuted data (Figure 3C).

Figure 3
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3

Redundancy analysis of SNP variation explained by climatic and spatial variables. Multivariate redundancy analysis showing the proportion of genotypic variation explained by climate variables (precipitation and temperature) and space (latitude and longitude) in (A) global germplasm and (B) Nigerian germplasm. (C) Density plot showing the distribution of proportion of variation explained by climate collinear with space obtained for 1000 permutations. The green line indicates the position of the proportion of variation explained by climate collinear with space and the red line indicates the 95th percentile of the null distribution obtained by permutations.

Population structure analysis of the Nigerian germplasm

To characterize the genetic structure of Nigerian germplasm in relation to global sorghum diversity, we conducted PCA and DAPC. In the PCA, the first PC explained about 6% of the variation while the second PC explained about 4% of the variation. The Nigerian accessions formed mostly separate clusters relative to the global accessions. The West African and Nigerian germplasm clustered together in most cases (Fig. S6A). Neighbor joining analysis also showed that Nigerian accessions and West African accessions cluster together, separately from the rest of the global germplasm (Fig. S7A). Clustering by botanical race was also observed in the Nigerian germplasm (Fig. S6B and Fig. S7B).

DAPC analysis identified three genetic groups (Figure 4A–C). The DAPC groups were genetically differentiated from each other as follows: Group 1 vs. Group 2 (FST of 0.21), Group 1 vs. Group 3 (FST of 0.18), and Group 2 vs. Group 3 (FST of 0.22). Accessions in Group 2 originate from locations with higher precipitation than the accessions in Group 1 (P-value < 0.001) and Group 3 (P-value < 0.001) (Figure 4D). The average nucleotide diversity in 1kb windows for groups 1, 2, and 3 are 4.4 × 10−4, 3.3 × 10−4, and 3.8 × 10−4 respectively. Linkage disequilibrium (r2) level decayed to 0.1 at 30 kb in in Group 1, 80 kb in Group 2, and 90 kb in Group 3 (Fig. S8). Caudatum types are more predominant in Group 1, Guinea types are more predominant in Group 2, and Durra types are more predominant in Group 3 (Table 1).

Figure 4
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4

Genetic grouping of Nigerian germplasm in relation to the precipitation gradient. (A) Three genetic groups identified by the discriminant analysis of principal components (DAPC), (B) neighbor joining tree showing genetic relatedness among Nigerian accessions belonging to the three DAPC groups, (C) distribution of the georeferenced accessions based on their DAPC genetic groups across precipitation zones in Nigeria, and (D) density plots of precipitation distribution of georeferenced accessions based on their DAPC genetic groups. Group 2 accessions originate from locations with higher precipitation than Group 1 (P-value < 0.001) and Group 3 (P-value < 0.001).

Botanical RaceGroup 1Group 2Group 3
Bicolor631
Caudatum3976
Durra413
Guinea121790
Bicolor Caudatum953
Bicolor Durra710
Bicolor Guinea2132
Caudatum Durra281157
Guinea Caudatum554
Guinea Durra111
Table 1 Distribution of sorghum botanical races among DAPC groups

In order to characterize variation among genetic groups identified by DAPC, ANOVA, and Tukey test were performed for putative adaptive traits and environmental variables. Significant difference in the distribution of precipitation and temperature gradient were found between DAPC Groups 1 and 3 (P-value < 0.01) and Groups 2 and 3 (P-value < 0.01). Also, there were significant differences between Groups 2 and 1 (P-value < 0.001) and Groups 2 and 3 (P-value < 0.001). We also found significant differences for putative adaptive traits between DAPC groups. For panicle length, all DAPC groups comparisons were different at P-value < 0.01. Significant differences were found for flowering time distribution between Groups 1 and 2 (P-value < 0.001). For plant height, significant differences were found between Groups 1 and 2 (P-value < 0.001) and Groups 2 and 3 (P-value < 0.001).

Genome-wide association studies of phenotypes

To identify genomic regions associated with phenotypic variation, we conducted MLMM GWAS for panicle length, plant height, and flowering time using BLUPs of common-garden phenotypes. Several genomic regions associated with the traits were identified and two a priori candidate genes fell within 100 kb (Figure 5, File S5). For panicle length no associations were significant at the Bonferroni threshold (Figure 5A). For plant height, a single significant association was observed on chromosome 3 (Figure 5B). The single significant association for plant height (S3_62675143, MAF = 0.29) colocalized with photoperiodic flowering gene SbCN12 (Sobic.003G295300, 73 kb away) (Yang et al. 2014). For flowering time, nine significant associations were found on chromosomes 3, 6, 7, 9, and 10 (Figure 5C). The most significant flowering time association (S6_799609, MAF = 0.09) colocalized with the known sorghum flowering time and photoperiod sensitivity gene Maturity6 (Ma6/Ghd7, Sobic.006G004400, 99 kb from gene) (Murphy et al. 2014). With a naive model (GLM GWAS), nominally significant associations were found on all chromosomes for panicle length, plant height, and flowering time, and a large number of these colocalized with a priori candidate genes (Fig. S9).

Figure 5
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5

Genome wide association studies of phenotypes for putative adaptation traits. Genome-wide SNP associations for multi-locus mixed model with a fixed population term and a random polygenic term for (A) panicle length, (B) plant height, and (C) flowering time in Nigerian germplasm (n = 329, 331, and 411 respectively). Vertical dashed lines represent position of a priori candidate genes (Ma6 and SbCN12) and found within 150 kb region of the significant association. Horizontal broken red lines signify the Bonferroni correction threshold (P-value < 0.05).

Genome scans for adaptation

Associations with environmental variables were used to investigate possible genomic signature of climate adaptation. GLM outliers (top 1% of associations) were identified on all chromosomes for “annual precipitation”, “precipitation in the driest quarter,” and “annual temperature” (Figure 6 A-D and File S6). Genome-wide, 15% and 6% of a priori candidate genes were localized within 100 kb of a 1% outlier SNP for “annual precipitation” and “precipitation in the driest quarter”, respectively (vs. 17% and 16% of all genes). Therefore, there is no enrichment of a priori candidate genes near environmental association outliers. A few of the a priori candidate genes that localize near 1% outlier SNPs are as follows. The SNP S9_54870238 (MAF = 0.36; 99th percentile) associated with annual precipitation was 90 kb away from the regulator of photoperiodic flowering SbCN8 (Centroradialis8, Sobic.009G199800). S9_8022437 (MAF = 0.49; 99th percentile) associated with annual precipitation is 94 kb from the sorghum ortholog (Sobic.009G069700) of maize barren inflorescence4 (bif4) (Galli et al. 2015). S3_4891237 and S3_4750963 (MAF = 0.31, 98th percentile and MAF = 0.04, 99th percentile) associated with annual precipitation and precipitation in the driest quarter were 100 kb from Sbra2 (Sobic.003G052900), the sorghum ortholog of maize inflorescence gene ramosa2 (ra2) (Brown et al. 2006).

Figure 6
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6

Genome wide association studies of climatic variables. Genome-wide SNP associations for general linear model of (A) annual mean precipitation, (B) precipitation in the driest quarter, (C) and annual mean temperature. Horizontal broken lines indicate the cutoff point of top 1% outliers. Vertical broken lines indicate the position of a priori candidate genes related to inflorescence morphology (orange), plant height (blue), and flowering time (purple), which are within 100 kb from associated SNPs.

Genome-wide pattern of Tajima’s D in the Nigerian germplasm (DNigeria) across 1 kb windows ranged between -2.0 to 4.0 (Fig. S10 and Fig. S11). The average Tajima’s D value in the Nigerian germplasm (DNigeria) was -0.2 while the average genome wide Tajima’s D (across 1 kb windows) in the global germplasm (DGlobal) and West African germplasm (DWestAfrica) were 0.1 and 0.2, respectively. Positive DNigeria windows were significantly enriched for a priori candidate genes compared to the expectation under a null distribution (Fig. S12 and File S7). Some of the DNigeria windows (Fig. S10) contain genes that control for flowering time and inflorescence development (File S7). For instance, the sorghum flowering time gene Maturity6 (Ma6, Sobic.006G004400) colocalized with the genome wide DNigeria scan window at 0.697 Mb (Tajima’s D = 1.48, 89th percentile) on chromosome 6. Maturity1 (Ma1, Sobic.006G057866) colocalized with the genome wide DNigeria scan window (Tajima’s D = 2.9; in the 98th percentile) around 40 Mb on chromosome 6. The sorghum ortholog of branched silkless1 (Sobic.002G411000), a maize spikelet meristem identity gene (Chuck et al. 2002), colocalized with the genome wide DNigeria scan window (Tajima’s D = 2.2 in the 94th percentile) around 75.9 Mb on chromosome 2.

Discussion

Evidence for clinal adaptation to the precipitation gradient

Genome-wide studies of nucleotide variation can provide insights into patterns of genetic variation in crop landraces and the role of clinal adaptation in shaping this variation (Meyer et al. 2016). Overall, we found several lines of evidence that sorghum phenotypic variation across Nigeria has been shaped by clinal adaptation to precipitation. For panicle length, common garden variation in the Nigerian sorghum germplasm was correlated with annual precipitation (Figure 2, Fig. S1). Sorghum accessions originating from lower latitudes that have high precipitation had longer panicles than accessions originating from higher latitudes that have less precipitation (Figure 2A, 2G). The long panicle morphology is associated with open and lax primary branches, which is a key feature of guinea race which are predominant in humid to sub-humid regions of Nigeria (Figure 1A-B) and West Africa more generally (Deu et al. 2008; Barro-Kondombo et al. 2010; Lasky et al. 2015). This open panicle morphology is thought to allow airflow, reducing mold infection under high humidity (Doggett 1988), though this model has not been formally tested in diverse germplasm.

For plant height, sorghum accessions from lower latitudes associated with high precipitation were taller than sorghum accessions from higher latitudes associated with less precipitation (Figure 2B, 2H). This pattern is consistent with cross-species ecological studies of plant height, which identified precipitation as the best environmental predictor of within-species latitudinal variation of height (Moles et al. 2009). Given that higher latitudes in Nigeria have lower rainfall, reduced plant height and panicle length in dry regions may be an adaptation to increase yield stability under reduced water availability, as has been observed in West African pearl millet (Vigouroux et al. 2011). In our study, common garden variation in flowering time is associated with temperature but not precipitation at the location of origin. The negative relationship between flowering time and annual mean temperature (Figure 2F) suggests that sorghum in hot climates may flower early as an escape from high temperature and resulting water limitation (Tuinstra et al. 1997).

The significant proportion of nucleotide variation explained by climate collinear with space (3%; Figure 3) is consistent with clinal adaptation of sorghum in Nigeria. Redundancy analysis indicated that climate collinear with space explained more SNP variation in the Nigerian germplasm than either of climate and space (isolation-by-distance) alone. The finding that climate collinear with space explained more SNP variation than either climate or space alone is consistent with findings in global sorghum germplasm (Lasky et al. 2015) and regional germplasm in wild soybean (Glycine soja) (Leamy et al. 2016) and barley (Hordeum vulgare) landraces (Abebe et al. 2015). However, the proportion of SNP variation explained by climate collinear with space we observed in this study was much lower than what was observed in wild soybean and barley landraces (6–34% and 29–61%, respectively). Methodological differences that may contribute to the lower proportion of SNP variation explained in this study are the greater number of environmental variables used in the soybean and barley studies, and the use of ascertained SNPs in the soybean study.

Genome-wide nucleotide variation in Nigeria is also structured according to precipitation zones. The DAPC analysis identified genetic groups within the sorghum botanical races. These genetic groups showed differences in precipitation distributions (Figure 4C-D, Table 1). Group 1 was associated with the lowest annual mean precipitation, and composed predominantly of caudatum, caudatum intermediates and bicolor types and prevalent at higher latitudes in northeastern Nigeria characterized with lower annual precipitation. The northeastern part of Nigeria was classified as part of the center of diversity of caudatum and caudatum-durra (Harlan 1992). Group 2 was associated with higher annual precipitation distribution and more prevalent at lower latitudes. Most of the accessions in this group belong to the guinea and guinea intermediate racial types. Group 3 was predominantly made up of durra and caudatum-durra intermediates. Notably, there was a complete absence of the guinea race from this group. Consistent with the model that botanical races in sorghum are differentially adapted to precipitation-based agroclimatic zones (Harlan 1992; Morris et al. 2013), we found differences in botanical race distribution in precipitation zones (File S7). Precipitation distribution of guinea accessions was significantly different from precipitation distribution of caudatum (P-value < 0.05) and caudatum-durra races (P-value < 0.01).

Despite some evidence of clinal adaptation, we should note that the redundancy analysis (Figure 3) and theoretical considerations (Pavlidis et al. 2012; Meirmans 2015) suggest that bulk of the nucleotide variation observed in the Nigerian germplasm is neutral. In addition, the small proportion of phenotypic and genomic variation explained by the precipitation gradient (Figure 3B-C; Figure 2) suggest a modest role of clinal adaptation in shaping phenotypic diversity at the geographic scale of modern Nigeria. However, more detailed studies of panicle and vegetative morphology may reveal traits with stronger climate associations. Other factors that may have shaped the observed diversity patterns which we did not investigate include seed sharing based on ethnolinguistic grouping (Soler et al. 2013; Labeyrie et al. 2014) and historical processes of domestication and diffusion (Kimber et al. 2013; Morris et al. 2013). Consistent with a major role of cultural and historical factors shaping diversity at a Nigeria-wide scale, in the cases where multiple accessions were collected from single locations, multiple botanical races and genetic group were observed (Figure 1A, File S1).

Identifying putative loci underlying adaptation

The combination of phenotypic association, environmental association, and selection scans can provide multiple lines of evidence for the involvement of particular loci in adaptation (Meyer et al. 2016). In the Nigerian sorghum germplasm, we observed some cases where a priori candidate genes related to adaptive traits colocalized with GWAS signals (Figure 5-6; File S5-S6). The colocalization of two sorghum photoperiod sensitivity genes (SbCN12 and Ma6) with plant height and flowering time associations suggests that photoperiod sensitivity contributes to plant height and flowering time adaptation in Nigerian sorghum germplasm. This is consistent with previous studies that identified SbCN12 and Ma6/Ghd7 as major genes underlying natural variation in photoperiodic flowering in sorghum (Murphy et al. 2014; Yang et al. 2014; Bouchet et al. 2017). In sub-tropical latitudes, like the common gardens used by the USDA genebank (∼17° N), photoperiod sensitive sorghums are expected to have longer vegetative growth and attain greater heights than in tropical latitudes due to longer days at the higher latitudes (Murphy et al. 2014).

Given that strong association peaks for both environmental variables (annual precipitation and precipitation in the driest quarter; Figure 6) were found near SbCN8, this gene may be a candidate for clinal adaptation to the precipitation gradient in West Africa. The colocalization of inflorescence genes (ra2 and bif4) with associations for annual precipitation (Figure 6A; File S5-S6) is consistent with an adaptive role of inflorescence morphology across the precipitation gradient (Harlan 1992; Morris et al. 2013). Accounting for population structure in GWAS models when mapping phenotypes that are correlated with population structure can lead to false negatives (Bergelson and Roux 2010; Lasky et al. 2015). Notably there were no significant associations for panicle length after accounting for population structure (MLM; Figure 5A). This finding is consistent with the expectation that climate-associated traits will be confounded with population structure and there will be little power to map the genetic basis of these traits using mixed model association (Brachi et al. 2011). Genetic dissection of panicle morphology and other putative clinal adaptation traits should be more effective with a multi-parent mapping strategy that breaks up confounding population structure, such as nested association mapping (Buckler et al. 2009; Bouchet et al. 2017).

Long-term stable clinal adaptation is expected to be reflected in genomic signals of balancing selection (Novembre and Rienzo 2009; Yoder et al. 2014). Consistent with the phenotypic evidence of clinal trait adaptation (Figure 2A, 2B, 2F), there was evidence of balancing selection from genome-wide enrichment of DNigeria windows at a priori candidate genes for flowering time, plant height, and inflorescence morphology (Fig. S11 and Fig. S12). Further functional studies of candidate genes will be needed to establish if these candidate genes have a role in climate adaptation (Kesari et al. 2012; Romero Navarro et al. 2017). Resequencing of georeferenced germplasm should facilitate the identification of putative functional variants underlying clinal adaptation (Gore et al. 2009; Meyer et al. 2016). Demographic effects due to historical processes of domestication and differentiation affect the pattern of diversity under neutrality (Molina et al. 2011; Meyer and Purugganan 2013). While it is suspected that sorghum has a complex domestication history based on morphology (Harlan 1992), population structure (Deu et al. 2006), and shattering gene haplotypes (Shattering1) (Olsen 2012), the demographic scenarios have not been formally described or evaluated. Therefore, further study on demographic history in sorghum will be valuable to identify which genomic outlier loci are most likely to represent signatures of selection.

If variation at genes involved in putative adaptive traits (flowering time, plant height, and panicle length) underlie clinal variation, then we expect significant enrichment of the a priori genes with environment associations (Rellstab et al. 2015). However, there was no enrichment of a priori candidate genes near associations with “precipitation in the driest quarter” or “annual mean precipitation”. An alternative to the single-trait GWAS and colocalization approach used in this study would be a multi-trait GWAS approach (Korte et al. 2012; Zhou and Stephens 2014) simultaneously considering traits and environmental variables. However, biological interpretation of the synthetic or “eigen” traits may remain a challenge until more corroborating functional genetic data are available (Banerjee et al. 2008).

Resources for genomic-enabled breeding of clinal adaptation

The Nigerian germplasm harbors abundant nucleotide polymorphism (90% of the global nucleotide polymorphism based on the πNigeria vs. πglobal), consistent with West Africa as a center of diversity for sorghum (Harlan 1992). The high diversity could be a result of ancestral diversity in the Nigerian sorghum, and/or gene flow and diffusion from other regions of Africa (Westengen et al. 2014). Population structure and neighbor-joining tree analysis showed that majority of the Nigerian accessions and West African accessions clustered together (Fig. S6 and Fig. S7). The Nigerian germplasm was 10 times less differentiated from West African germplasm compared to the rest of global sorghum germplasm (FST Nigeria-West Africa = 0.007 and FST Nigeria-Global = 0.07). Low level of differentiation among regional germplasm in sorghum has been attributed to gene flow from human migration and agricultural trade (Menkir et al. 1997).

Characterization of LD patterns is critical for interpretation of genome scans since the local extent of LD decay determines resolution of mapping and long-range LD creates spurious associations (Myles et al. 2009). LD decay rate in Nigerian germplasm (half of maximum value at 12 kb) was slower than LD decay estimated in the global sorghum association panel (half of maximum value at 1 kb) (Morris et al. 2013). The global sorghum association panel may capture more historical recombination than the Nigerian germplasm because of its greater geographic diversity. The reduced long-range LD decay in the Nigerian germplasm compared to the West African germplasm (Fig. S5) may be due to the smaller geographic scale, which should reduce long range LD due to isolation-by-distance (Brachi et al. 2011). Given the observed LD decay rates, the mapping resolution for genome-wide scans in the Nigerian germplasm is expected to be less than in global sorghum panels but greater than in West African regional panels. Overall, the modest LD decay rate and high genetic diversity in the Nigerian germplasm make it suitable for genome-wide association studies.

The application of genomics for crop improvement and plant genetic resource management is still in the early stage for most national agricultural research systems in sub-Saharan Africa, including Nigeria and neighboring countries (Ezeaku and Gupta 2004; Leiser et al. 2014; Yohannes et al. 2015). The genomic resources developed in this study represent a step toward genomics-enabled breeding and germplasm management for sorghum landraces in Nigeria. The resources developed include a genome-wide catalog of SNP variation, a description of geographic population structure, and estimates of genetic properties including nucleotide diversity and LD decay. The genomic signatures of clinal adaptation identified in this study, if validated in managed stress and multi-environment mapping studies (Cooper et al. 2014; Lasky et al. 2015), could be another resource to facilitate genomics-enabled breeding for climate-resilience in West Africa.

Acknowledgements

This study is made possible by the support of the American People provided to the Feed the Future Innovation Lab for Collaborative Research on Sorghum and Millet through the United States Agency for International Development (USAID). The contents are the sole responsibility of the authors and do not necessarily reflect the views of USAID or the United States Government. Program activities are funded by the USAID under Cooperative Agreement No. AID-OAA-A-13-00047. We thank the editor and two anonymous reviewers for suggestions that improved the manuscript. This study made use of resources at the Beocat high-performance computing facility and Integrated Genomics Facility at Kansas State University. This study is contribution 18-633-J from the Kansas Agricultural Experiment Station.

Footnotes

  • Supplemental material available at Figshare: https://doi.org/10.25387/g3.6942986.

  • Communicating editor: P. Morrell

  • Received June 29, 2018.
  • Accepted August 4, 2018.
  • Copyright © 2018 Olatoye et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Literature Cited

  1. ↵
    1. Abebe T. D.,
    2. Naz A. A.,
    3. Léon J.
    , 2015 Landscape genomics reveal signatures of local adaptation in barley (Hordeum vulgare L.). Front. Plant Sci. 6: 813. doi:10.3389/fpls.2015.00813
    OpenUrlCrossRef
  2. ↵
    1. Banerjee S.,
    2. Yandell B. S.,
    3. Yi N.
    , 2008 Bayesian Quantitative Trait Loci Mapping for Multiple Traits. Genetics 179: 2275–2289. doi:10.1534/genetics.108.088427
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Barro-Kondombo C.,
    2. Sagnard F.,
    3. Chantereau J.,
    4. Deu M.,
    5. vom Brocke K.,
    6. et al.
    , 2010 Genetic structure among sorghum landraces as revealed by morphological variation and microsatellite markers in three agroclimatic regions of Burkina Faso. Theor. Appl. Genet. 120: 1511–1523. doi:10.1007/s00122-010-1272-2
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bastide H.,
    2. Lange J. D.,
    3. Lack J. B.,
    4. Yassin A.,
    5. Pool J. E.
    , 2016 A variable genetic architecture of melanic evolution in Drosophila melanogaster. Genetics genetics.116.192492. doi:10.1534/genetics.116.192492
    OpenUrlAbstract/FREE Full Text
  5. ↵
    Bates, D., M. Maechler, B. Bolker, S. Walker, R. H. B. Christensen et al., 2014 lme4: Linear mixed-effects models using Eigen and S4.
  6. ↵
    1. Bergelson J.,
    2. Roux F.
    , 2010 Towards identifying genes underlying ecologically relevant traits in Arabidopsis thaliana. Nat. Rev. Genet. 11: 867–879. doi:10.1038/nrg2896
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bouchet S.,
    2. Olatoye M. O.,
    3. Marla S. R.,
    4. Perumal R.,
    5. Tesso T.,
    6. et al.
    , 2017 Increased power to dissect adaptive traits in global sorghum diversity using a nested association mapping population. Genetics 206: 573–585. doi:10.1534/genetics.116.198499
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Brachi B.,
    2. Morris G. P.,
    3. Borevitz J. O.
    , 2011 Genome-wide association studies in plants: the missing heritability is in the field. Genome Biol. 12: 232. doi:10.1186/gb-2011-12-10-232
    OpenUrlCrossRefPubMed
  9. ↵
    1. Bradbury P. J.,
    2. Zhang Z.,
    3. Kroon D. E.,
    4. Casstevens T. M.,
    5. Ramdoss Y.,
    6. et al.
    , 2007 TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23: 2633–2635. doi:10.1093/bioinformatics/btm308
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    1. Brown P.,
    2. Klein P.,
    3. Bortiri E.,
    4. Acharya C.,
    5. Rooney W.,
    6. et al.
    , 2006 Inheritance of inflorescence architecture in sorghum. Theor. Appl. Genet. 113: 931–942. doi:10.1007/s00122-006-0352-9
    OpenUrlCrossRefPubMedWeb of Science
  11. ↵
    1. Browning B. L.,
    2. Browning S. R.
    , 2013 Improving the accuracy and efficiency of identity-by-descent detection in population data. Genetics 194: 459–471. doi:10.1534/genetics.113.150029
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Buckler E. S.,
    2. Holland J. B.,
    3. Bradbury P. J.,
    4. Acharya C. B.,
    5. Brown P. J.,
    6. et al.
    , 2009 The genetic architecture of maize flowering time. Science 325: 714–718. doi:10.1126/science.1174276
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Chuck G.,
    2. Muszynski M.,
    3. Kellogg E.,
    4. Hake S.,
    5. Schmidt R. J.
    , 2002 The control of spikelet meristem identity by the branched silkless1 gene in maize. Science 298: 1238–1241. doi:10.1126/science.1076920
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Comadran J.,
    2. Kilian B.,
    3. Russell J.,
    4. Ramsay L.,
    5. Stein N.,
    6. et al.
    , 2012 Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat. Genet. 44: 1388–1392. doi:10.1038/ng.2447
    OpenUrlCrossRefPubMed
  15. ↵
    1. Cooper M.,
    2. Messina C. D.,
    3. Podlich D.,
    4. Totir L. R.,
    5. Baumgarten A.,
    6. et al.
    , 2014 Predicting the future of plant breeding: complementing empirical evaluation with genetic prediction. Crop Pasture Sci. 65: 311–336. doi:10.1071/CP14007
    OpenUrlCrossRefWeb of Science
  16. ↵
    1. Danecek P.,
    2. Auton A.,
    3. Abecasis G.,
    4. Albers C. A.,
    5. Banks E.,
    6. et al.
    , 2011 The variant call format and VCFtools. Bioinformatics 27: 2156–2158. doi:10.1093/bioinformatics/btr330
    OpenUrlCrossRefPubMedWeb of Science
  17. ↵
    1. Debieu M.,
    2. Tang C.,
    3. Stich B.,
    4. Sikosek T.,
    5. Effgen S.,
    6. et al.
    , 2013 Co-variation between seed dormancy, growth rate and flowering time changes with latitude in Arabidopsis thaliana. PLoS One 8: e61075. doi:10.1371/journal.pone.0061075
    OpenUrlCrossRefPubMed
  18. ↵
    1. Deu M.,
    2. Rattunde F.,
    3. Chantereau J.
    , 2006 A global view of genetic diversity in cultivated sorghums using a core collection. Genome 49: 168–180. doi:10.1139/g05-092
    OpenUrlCrossRefPubMed
  19. ↵
    1. Deu M.,
    2. Sagnard F.,
    3. Chantereau J.,
    4. Calatayud C.,
    5. Hérault D.,
    6. et al.
    , 2008 Niger-wide assessment of in situ sorghum genetic diversity with microsatellite markers. Theor. Appl. Genet. 116: 903–913. doi:10.1007/s00122-008-0721-7
    OpenUrlCrossRefPubMed
  20. ↵
    1. Djè Y.,
    2. Heuertz M.,
    3. Lefèbvre C.,
    4. Vekemans X.
    , 2000 Assessment of genetic diversity within and among germplasm accessions in cultivated sorghum using microsatellite markers. Theor. Appl. Genet. 100: 918–925. doi:10.1007/s001220051371
    OpenUrlCrossRefWeb of Science
  21. ↵
    1. Doggett H.
    , 1988 Sorghum, 2d. ed. Tropical agricultural series. Longman Scientific, Essex, UK.
  22. ↵
    1. Ducrocq S.,
    2. Madur D.,
    3. Veyrieras J.-B.,
    4. Camus-Kulandaivelu L.,
    5. Kloiber-Maitz M.,
    6. et al.
    , 2008 Key impact of Vgt1 on flowering time adaptation in maize: Evidence from association mapping and ecogeographical information. Genetics 178: 2433–2437. doi:10.1534/genetics.107.084830
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Dwivedi S. L.,
    2. Ceccarelli S.,
    3. Blair M. W.,
    4. Upadhyaya H. D.,
    5. Are A. K.,
    6. et al.
    , 2016 Landrace germplasm for improving yield and abiotic stress adaptation. Trends Plant Sci. 21: 31–42. doi:10.1016/j.tplants.2015.10.012
    OpenUrlCrossRef
  24. ↵
    1. Elshire R. J.,
    2. Glaubitz J. C.,
    3. Sun Q.,
    4. Poland J. A.,
    5. Kawamoto K.,
    6. et al.
    , 2011 A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One 6: e19379. doi:10.1371/journal.pone.0019379
    OpenUrlCrossRefPubMed
  25. ↵
    1. Ezeaku I. E.,
    2. Gupta S. C.
    , 2004 Development of sorghum populations for resistance to Striga hermonthica in the Nigerian Sudan Savanna. Afr. J. Biotechnol. 3: 324–329. doi:10.5897/AJB2004.000-2059
    OpenUrlCrossRef
  26. ↵
    1. Galli M.,
    2. Liu Q.,
    3. Moss B. L.,
    4. Malcomber S.,
    5. Li W.,
    6. et al.
    , 2015 Auxin signaling modules regulate maize inflorescence architecture. Proc. Natl. Acad. Sci. USA 112: 13372–13377. doi:10.1073/pnas.1516473112
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Glaubitz J. C.,
    2. Casstevens T. M.,
    3. Lu F.,
    4. Harriman J.,
    5. Elshire R. J.,
    6. et al.
    , 2014 TASSEL-GBS: A high capacity genotyping by sequencing analysis pipeline. PLoS One 9: e90346. doi:10.1371/journal.pone.0090346
    OpenUrlCrossRefPubMed
  28. ↵
    1. Goodstein D. M.,
    2. Shu S.,
    3. Howson R.,
    4. Neupane R.,
    5. Hayes R. D.,
    6. et al.
    , 2011 Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 40: D1178–D1186. doi:10.1093/nar/gkr944
    OpenUrlCrossRefPubMedWeb of Science
  29. ↵
    1. Gore M. A.,
    2. Chia J.-M.,
    3. Elshire R. J.,
    4. Sun Q.,
    5. Ersoz E. S.,
    6. et al.
    , 2009 A first-generation haplotype map of maize. Science 326: 1115–1117. doi:10.1126/science.1177837
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Hancock A. M.,
    2. Witonsky D. B.,
    3. Alkorta-Aranburu G.,
    4. Beall C. M.,
    5. Gebremedhin A.,
    6. et al.
    , 2011 Adaptations to climate-mediated selective pressures in humans. PLoS Genet. 7: e1001375. doi:10.1371/journal.pgen.1001375
    OpenUrlCrossRefPubMed
  31. ↵
    1. Harlan J. R.
    , 1971 Agricultural origins: centers and noncenters. Science 174: 468–474. doi:10.1126/science.174.4008.468
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Harlan J. R.
    , 1992 Crops and Man, American Society of Agronomy, Madison, Wisconsin.
  33. ↵
    1. Haussmann B. I. G.,
    2. Fred Rattunde H.,
    3. Weltzien-Rattunde E.,
    4. Traoré P. S. C.,
    5. vom Brocke K.,
    6. et al.
    , 2012 Breeding strategies for adaptation of pearl millet and sorghum to climate variability and change in West Africa. J. Agron. Crop Sci. 198: 327–339. doi:10.1111/j.1439-037X.2012.00526.x
    OpenUrlCrossRef
  34. ↵
    Hijmans, R. J., J. van Etten, J. Cheng, M. Mattiuzzi, M. Sumner et al., 2016 raster: Geographic Data Analysis and Modeling.
  35. ↵
    1. Hufford M. B.,
    2. Xu X.,
    3. van Heerwaarden J.,
    4. Pyhäjärvi T.,
    5. Chia J.-M.,
    6. et al.
    , 2012 Comparative population genomics of maize domestication and improvement. Nat. Genet. 44: 808–811. doi:10.1038/ng.2309
    OpenUrlCrossRefPubMed
  36. ↵
    1. Jombart T.
    , 2008 Adegenet: A R Package for the Multivariate Analysis of Genetic Markers. Bioinformatics 24: 1403–1405. doi:10.1093/bioinformatics/btn129
    OpenUrlCrossRefPubMedWeb of Science
  37. ↵
    1. Jombart T.,
    2. Devillard S.,
    3. Balloux F.
    , 2010 Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11: 94. doi:10.1186/1471-2156-11-94
    OpenUrlCrossRefPubMed
  38. ↵
    1. Kesari R.,
    2. Lasky J. R.,
    3. Villamor J. G.,
    4. Marais D. L. D.,
    5. Chen Y.-J. C.,
    6. et al.
    , 2012 Intron-mediated alternative splicing of Arabidopsis P5CS1 and its association with natural variation in proline and climate adaptation. Proc. Natl. Acad. Sci. USA 109: 9197–9202. doi:10.1073/pnas.1203433109
    OpenUrlAbstract/FREE Full Text
  39. ↵
    Kimber, C. T., J. A. Dahlberg, and S. Kresovich, 2013 The Gene Pool of Sorghum bicolor and Its Improvement, pp. 23–41 in Genomics of the Saccharinae, edited by A. H. Paterson. Plant Genetics and Genomics: Crops and Models 11, Springer, New York.
  40. ↵
    1. Kloosterman B.,
    2. Abelenda J. A.,
    3. Gomez Mdel M.,
    4. Oortwijn M.,
    5. de Boer J. M.,
    6. et al.
    , 2013 Naturally occurring allele diversity allows potato cultivation in northern latitudes. Nature 495: 246–250. doi:10.1038/nature11912
    OpenUrlCrossRefPubMedWeb of Science
  41. ↵
    1. Korte A.,
    2. Vilhjálmsson B. J.,
    3. Segura V.,
    4. Platt A.,
    5. Long Q.,
    6. et al.
    , 2012 A mixed-model approach for genome-wide association studies of correlated traits in structured populations. Nat. Genet. 44: 1066–1071. doi:10.1038/ng.2376
    OpenUrlCrossRefPubMed
  42. ↵
    1. Kronholm I.,
    2. Picó F. X.,
    3. Alonso-Blanco C.,
    4. Goudet J.,
    5. de Meaux J.
    , 2012 Genetic basis of adaptation in Arabidopsis thaliana: Local adaptation at the seed dormancy QTL DOG1. Evolution 66: 2287–2302. doi:10.1111/j.1558-5646.2012.01590.x
    OpenUrlCrossRefPubMedWeb of Science
  43. ↵
    1. Labeyrie V.,
    2. Deu M.,
    3. Barnaud A.,
    4. Calatayud C.,
    5. Buiron M.,
    6. et al.
    , 2014 Influence of ethnolinguistic diversity on the sorghum genetic patterns in subsistence farming systems in eastern Kenya. PLoS One 9: e92178. doi:10.1371/journal.pone.0092178
    OpenUrlCrossRefPubMed
  44. ↵
    1. Lasky J. R.,
    2. Des Marais D. L.,
    3. McKay J. K.,
    4. Richards J. H.,
    5. Juenger T. E.,
    6. et al.
    , 2012 Characterizing genomic variation of Arabidopsis thaliana: the roles of geography and climate. Mol. Ecol. 21: 5512–5529. doi:10.1111/j.1365-294X.2012.05709.x
    OpenUrlCrossRefPubMedWeb of Science
  45. ↵
    1. Lasky J. R.,
    2. Upadhyaya H. D.,
    3. Ramu P.,
    4. Deshpande S.,
    5. Hash C. T.,
    6. et al.
    , 2015 Genome-environment associations in sorghum landraces predict adaptive traits. Sci. Adv. 1: e1400218. doi:10.1126/sciadv.1400218
    OpenUrlFREE Full Text
  46. ↵
    1. Leamy L. J.,
    2. Lee C.-R.,
    3. Song Q.,
    4. Mujacic I.,
    5. Luo Y.,
    6. et al.
    , 2016 Environmental vs. geographical effects on genomic variation in wild soybean (Glycine soja) across its native range in northeast Asia. Ecol. Evol. 6: 6332–6344. doi:10.1002/ece3.2351
    OpenUrlCrossRef
  47. ↵
    1. Leiser W. L.,
    2. Rattunde H. F.,
    3. Weltzien E.,
    4. Cisse N.,
    5. Abdou M.,
    6. et al.
    , 2014 Two in one sweep: aluminum tolerance and grain yield in P-limited soils are associated to the same genomic region in West African Sorghum. BMC Plant Biol. 14: 206. doi:10.1186/s12870-014-0206-6
    OpenUrlCrossRefPubMed
  48. ↵
    1. Li H.,
    2. Durbin R.
    , 2009 Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25: 1754–1760. doi:10.1093/bioinformatics/btp324
    OpenUrlCrossRefPubMedWeb of Science
  49. ↵
    1. Li L.-F.,
    2. Li Y.-L.,
    3. Jia Y.,
    4. Caicedo A. L.,
    5. Olsen K. M.
    , 2017 Signatures of adaptation in the weedy rice genome. Nat. Genet. 49: 811–814. doi:10.1038/ng.3825
    OpenUrlCrossRef
  50. ↵
    1. Lipka A. E.,
    2. Tian F.,
    3. Wang Q.,
    4. Peiffer J.,
    5. Li M.,
    6. et al.
    , 2012 GAPIT: genome association and prediction integrated tool. Bioinformatics 28: 2397–2399. doi:10.1093/bioinformatics/bts444
    OpenUrlCrossRefPubMedWeb of Science
  51. ↵
    1. Lobell D. B.,
    2. Burke M. B.,
    3. Tebaldi C.,
    4. Mastrandrea M. D.,
    5. Falcon W. P.,
    6. et al.
    , 2008 Prioritizing climate change adaptation needs for food security in 2030. Science 319: 607–610. doi:10.1126/science.1152339
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Ma Y.,
    2. Dai X.,
    3. Xu Y.,
    4. Luo W.,
    5. Zheng X.,
    6. et al.
    , 2015 COLD1 confers chilling tolerance in rice. Cell 160: 1209–1221. doi:10.1016/j.cell.2015.01.046
    OpenUrlCrossRefPubMed
  53. ↵
    1. Mace E. S.,
    2. Li Y.,
    3. Prentis P. J.,
    4. Bian L.,
    5. Campbell B. C.,
    6. et al.
    , 2013 Whole-genome sequencing reveals untapped genetic potential in Africa’s indigenous cereal crop sorghum. Nat. Commun. 4: 2320. doi:10.1038/ncomms3320
    OpenUrlCrossRefPubMed
  54. ↵
    1. McCormick R. F.,
    2. Truong S. K.,
    3. Sreedasyam A.,
    4. Jenkins J.,
    5. Shu S.,
    6. et al.
    , 2018 The Sorghum bicolor reference genome: improved assembly, gene annotations, a transcriptome atlas, and signatures of genome organization. Plant J. 93: 338–354. doi:10.1111/tpj.13781
    OpenUrlCrossRef
  55. ↵
    1. Meirmans P. G.
    , 2015 Seven common mistakes in population genetics and how to avoid them. Mol. Ecol. 24: 3223–3231. doi:10.1111/mec.13243
    OpenUrlCrossRef
  56. ↵
    1. Menkir A.,
    2. Goldsbrough P.,
    3. Ejeta G.
    , 1997 RAPD based assessment of genetic diversity in cultivated races of sorghum. Crop Sci. 37: 564–569. doi:10.2135/cropsci1997.0011183X003700020042x
    OpenUrlCrossRef
  57. ↵
    1. Meyer R. S.,
    2. Choi J. Y.,
    3. Sanches M.,
    4. Plessis A.,
    5. Flowers J. M.,
    6. et al.
    , 2016 Domestication history and geographical adaptation inferred from a SNP map of African rice. Nat. Genet. 48: 1083–1088. doi:10.1038/ng.3633
    OpenUrlCrossRef
  58. ↵
    1. Meyer R. S.,
    2. Purugganan M. D.
    , 2013 Evolution of crop species: genetics of domestication and diversification. Nat. Rev. Genet. 14: 840–852. doi:10.1038/nrg3605
    OpenUrlCrossRefPubMed
  59. ↵
    1. Moles A. T.,
    2. Warton D. I.,
    3. Warman L.,
    4. Swenson N. G.,
    5. Laffan S. W.,
    6. et al.
    , 2009 Global patterns in plant height. J. Ecol. 97: 923–932. doi:10.1111/j.1365-2745.2009.01526.x
    OpenUrlCrossRef
  60. ↵
    1. Molina J.,
    2. Sikora M.,
    3. Garud N.,
    4. Flowers J. M.,
    5. Rubinstein S.,
    6. et al.
    , 2011 Molecular evidence for a single evolutionary origin of domesticated rice. Proc. Natl. Acad. Sci. USA 108: 8351–8356. doi:10.1073/pnas.1104686108
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Morrell P. L.,
    2. Buckler E. S.,
    3. Ross-Ibarra J.
    , 2012 Crop genomics: advances and applications. Nat. Rev. Genet. 13: 85–96. doi:10.1038/nrg3097
    OpenUrlCrossRef
  62. ↵
    1. Morris G. P.,
    2. Ramu P.,
    3. Deshpande S. P.,
    4. Hash C. T.,
    5. Shah T.,
    6. et al.
    , 2013 Population genomic and genome-wide association studies of agroclimatic traits in sorghum. Proc. Natl. Acad. Sci. USA 110: 453–458. doi:10.1073/pnas.1215985110
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Murphy R. L.,
    2. Morishige D. T.,
    3. Brady J. A.,
    4. Rooney W. L.,
    5. Yang S.,
    6. et al.
    , 2014 Ghd7 (Ma6) represses sorghum flowering in long days: Ghd7 alleles enhance biomass accumulation and grain production. Plant Genome 7: 1–10. doi:10.3835/plantgenome2013.11.0040
    OpenUrlCrossRef
  64. ↵
    1. Myles S.,
    2. Peiffer J.,
    3. Brown P. J.,
    4. Ersoz E. S.,
    5. Zhang Z.,
    6. et al.
    , 2009 Association mapping: critical considerations shift from genotyping to experimental design. Plant Cell 21: 2194–2202. doi:10.1105/tpc.109.068437
    OpenUrlAbstract/FREE Full Text
  65. ↵
    1. National Research Council
    , 1996 Lost Crops of Africa: Volume I: Grains. National Academy Press, Washington, D.C.
  66. ↵
    1. Novembre J.,
    2. Rienzo A. D.
    , 2009 Spatial patterns of variation due to natural selection in humans. Nat. Rev. Genet. 10: 745–755. doi:10.1038/nrg2632
    OpenUrlCrossRefPubMedWeb of Science
  67. ↵
    1. Nzeka U.,
    2. Akhidenor J.
    , 2018 Nigeria: Grain and Feed Annual. GAIN Reports, USDA-FAS. https://gain.fas.usda.gov/Recent GAIN Publications/Grain and Feed Annual_Lagos_Nigeria_4-12-2018.pdf
  68. ↵
    Oksanen, J., F. G. Blanchet, M. Friendly, R. Kindt, P. Legendre et al., 2017 vegan: Community Ecology Package.
  69. ↵
    1. Olsen K. M.
    , 2012 One gene’s shattering effects. Nat. Genet. 44: 616–617. doi:10.1038/ng.2289
    OpenUrlCrossRefPubMed
  70. ↵
    1. Paradis E.,
    2. Claude J.,
    3. Strimmer K.
    , 2004 APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20: 289–290. doi:10.1093/bioinformatics/btg412
    OpenUrlCrossRefPubMedWeb of Science
  71. ↵
    1. Pavlidis P.,
    2. Jensen J. D.,
    3. Stephan W.,
    4. Stamatakis A.
    , 2012 A critical assessment of storytelling: Gene ontology categories and the importance of validating genomic scans. Mol. Biol. Evol. 29: 3237–3248. doi:10.1093/molbev/mss136
    OpenUrlCrossRefPubMedWeb of Science
  72. ↵
    1. Rao K. E. P.,
    2. Obilana A. T.,
    3. Mengesha M. H.
    , 1985 Collection of kaura, fara-fara and guineense sorghums in northern Nigeria. J. Agric. Tradit. Bot. Appl. 32: 73–81.
    OpenUrl
  73. ↵
    1. Rellstab C.,
    2. Gugerli F.,
    3. Eckert A. J.,
    4. Hancock A. M.,
    5. Holderegger R.
    , 2015 A practical guide to environmental association analysis in landscape genomics. Mol. Ecol. 24: 4348–4370. doi:10.1111/mec.13322
    OpenUrlCrossRef
  74. ↵
    1. Romero Navarro J. A.,
    2. Willcox M.,
    3. Burgueño J.,
    4. Romay C.,
    5. Swarts K.,
    6. et al.
    , 2017 A study of allelic diversity underlying flowering-time adaptation in maize landraces. Nat. Genet. 49: 476–480. doi:10.1038/ng.3784
    OpenUrlCrossRefPubMed
  75. ↵
    1. Samis K. E.,
    2. Murren C. J.,
    3. Bossdorf O.,
    4. Donohue K.,
    5. Fenster C. B.,
    6. et al.
    , 2012 Longitudinal trends in climate drive flowering time clines in North American Arabidopsis thaliana. Ecol. Evol. 2: 1162–1180. doi:10.1002/ece3.262
    OpenUrlCrossRefPubMed
  76. ↵
    1. Savolainen O.,
    2. Lascoux M.,
    3. Merilä J.
    , 2013 Ecological genomics of local adaptation. Nat. Rev. Genet. 14: 807–820. doi:10.1038/nrg3522
    OpenUrlCrossRefPubMed
  77. ↵
    1. Segura V.,
    2. Vilhjálmsson B. J.,
    3. Platt A.,
    4. Korte A.,
    5. Seren Ü.,
    6. et al.
    , 2012 An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations. Nat. Genet. 44: 825–830. doi:10.1038/ng.2314
    OpenUrlCrossRefPubMed
  78. ↵
    1. Shanahan T. M.,
    2. Overpeck J. T.,
    3. Anchukaitis K. J.,
    4. Beck J. W.,
    5. Cole J. E.,
    6. et al.
    , 2009 Atlantic forcing of persistent drought in West Africa. Science 324: 377–380. doi:10.1126/science.1166352
    OpenUrlAbstract/FREE Full Text
  79. ↵
    1. Siol M.,
    2. Wright S. I.,
    3. Barrett S. C. H.
    , 2010 The population genomics of plant adaptation. New Phytol. 188: 313–332. doi:10.1111/j.1469-8137.2010.03401.x
    OpenUrlCrossRefPubMedWeb of Science
  80. ↵
    1. Slatkin M.
    , 1973 Gene flow and selection in a cline. Genetics 75: 733–756.
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Soler C.,
    2. Saidou A.-A.,
    3. Hamadou T. V. C.,
    4. Pautasso M.,
    5. Wencelius J.,
    6. et al.
    , 2013 Correspondence between genetic structure and farmers’ taxonomy – a case study from dry-season sorghum landraces in northern Cameroon. Plant Genet. Resour. 11: 36–49. doi:10.1017/S1479262112000342
    OpenUrlCrossRef
  82. ↵
    1. Tuinstra M. R.,
    2. Grote E. M.,
    3. Goldsbrough P. B.,
    4. Ejeta G.
    , 1997 Genetic analysis of post-flowering drought tolerance and components of grain development in Sorghum bicolor (L.). Moench. Mol. Breed. 3: 439–448. doi:10.1023/A:1009673126345
    OpenUrlCrossRef
  83. ↵
    1. Vasconcelos A. C. F.,
    2. Bonatti M.,
    3. Schlindwein S. L.,
    4. D’Agostini L. R.,
    5. Homem L. R.,
    6. et al.
    , 2013 Landraces as an adaptation strategy to climate change for smallholders in Santa Catarina, Southern Brazil. Land Use Policy 34: 250–254. doi:10.1016/j.landusepol.2013.03.017
    OpenUrlCrossRef
  84. ↵
    1. Vigouroux Y.,
    2. Mariac C.,
    3. Mita S. D.,
    4. Pham J.-L.,
    5. Gérard B.,
    6. et al.
    , 2011 Selection for earlier flowering crop associated with climatic variations in the Sahel. PLoS One 6: e19563. doi:10.1371/journal.pone.0019563
    OpenUrlCrossRefPubMed
  85. ↵
    1. Westengen O. T.,
    2. Okongo M. A.,
    3. Onek L.,
    4. Berg T.,
    5. Upadhyaya H.,
    6. et al.
    , 2014 Ethnolinguistic structuring of sorghum genetic diversity in Africa and the role of local seed systems. Proc. Natl. Acad. Sci. USA 111: 14100–14105. doi:10.1105/tpc.109.068437
    OpenUrlAbstract/FREE Full Text
  86. ↵
    1. Wu W.,
    2. Zheng X.-M.,
    3. Lu G.,
    4. Zhong Z.,
    5. Gao H.,
    6. et al.
    , 2013 Association of functional nucleotide polymorphisms at DTH2 with the northward expansion of rice cultivation in Asia. Proc. Natl. Acad. Sci. USA 110: 2775–2780. doi:10.1073/pnas.1213962110
    OpenUrlAbstract/FREE Full Text
  87. ↵
    1. Yang S.,
    2. Murphy R. L.,
    3. Morishige D. T.,
    4. Klein P. E.,
    5. Rooney W. L.,
    6. et al.
    , 2014 Sorghum phytochrome B inhibits flowering in long days by activating expression of SbPRR37 and SbGHD7, repressors of SbEHD1, SbCN8 and SbCN12. PLoS One 9: e105352. doi:10.1371/journal.pone.0105352
    OpenUrlCrossRefPubMed
  88. ↵
    1. Yoder J. B.,
    2. Stanton-Geddes J.,
    3. Zhou P.,
    4. Briskine R.,
    5. Young N. D.,
    6. et al.
    , 2014 Genomic signature of adaptation to climate in Medicago truncatula. Genetics 196: 1263–1275. doi:10.1534/genetics.113.159319
    OpenUrlAbstract/FREE Full Text
  89. ↵
    1. Yohannes T.,
    2. Abraha T.,
    3. Kiambi D.,
    4. Folkertsma R.,
    5. Hash C. T.,
    6. et al.
    , 2015 Marker-assisted introgression improves Striga resistance in an Eritrean farmer-preferred sorghum variety. Field Crops Res. 173: 22–29. doi:10.1016/j.fcr.2014.12.008
    OpenUrlCrossRef
  90. ↵
    1. Zeven A. C.
    , 1998 Landraces: A review of definitions and classifications. Euphytica 104: 127–139. doi:10.1023/A:1018683119237
    OpenUrlCrossRefWeb of Science
  91. ↵
    1. Zhang D.,
    2. Kong W.,
    3. Robertson J.,
    4. Goff V. H.,
    5. Epps E.,
    6. et al.
    , 2015 Genetic analysis of inflorescence and plant height components in sorghum (Panicoidae) and comparative genetics with rice (Oryzoidae). BMC Plant Biol. 15: 107. doi:10.1186/s12870-015-0477-6
    OpenUrlCrossRef
  92. ↵
    1. Zhen Y.,
    2. Ungerer M. C.
    , 2008 Clinal variation in freezing tolerance among natural accessions of Arabidopsis thaliana. New Phytol. 177: 419–427.
    OpenUrlPubMedWeb of Science
  93. ↵
    1. Zheng X.,
    2. Levine D.,
    3. Shen J.,
    4. Gogarten S. M.,
    5. Laurie C.,
    6. et al.
    , 2012 A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics 28: 3326–3328. doi:10.1093/bioinformatics/bts606
    OpenUrlCrossRefPubMedWeb of Science
  94. ↵
    1. Zhou X.,
    2. Stephens M.
    , 2014 Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat. Methods 11: 407–409. doi:10.1038/nmeth.2848
    OpenUrlCrossRefPubMed
View Abstract
Previous ArticleNext Article
Back to top

PUBLICATION INFORMATION

Volume 8 Issue 10, October 2018

G3: Genes|Genomes|Genetics: 8 (10)

ARTICLE CLASSIFICATION

Investigations
View this article with LENS
Email

Thank you for sharing this G3: Genes | Genomes | Genetics article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Genomic Signatures of Adaptation to a Precipitation Gradient in Nigerian Sorghum
(Your Name) has forwarded a page to you from G3: Genes | Genomes | Genetics
(Your Name) thought you would be interested in this article in G3: Genes | Genomes | Genetics.
Print
Alerts
Enter your email below to set up alert notifications for new article, or to manage your existing alerts.
SIGN UP OR SIGN IN WITH YOUR EMAIL
View PDF
Share

Genomic Signatures of Adaptation to a Precipitation Gradient in Nigerian Sorghum

Marcus O. Olatoye, Zhenbin Hu, Fanna Maina and View ORCID ProfileGeoffrey P. Morris
G3: Genes, Genomes, Genetics October 1, 2018 vol. 8 no. 10 3269-3281; https://doi.org/10.1534/g3.118.200551
Marcus O. Olatoye
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhenbin Hu
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fanna Maina
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Geoffrey P. Morris
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Geoffrey P. Morris
  • For correspondence: gpmorris@ksu.edu
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Citation

Genomic Signatures of Adaptation to a Precipitation Gradient in Nigerian Sorghum

Marcus O. Olatoye, Zhenbin Hu, Fanna Maina and View ORCID ProfileGeoffrey P. Morris
G3: Genes, Genomes, Genetics October 1, 2018 vol. 8 no. 10 3269-3281; https://doi.org/10.1534/g3.118.200551
Marcus O. Olatoye
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhenbin Hu
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fanna Maina
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Geoffrey P. Morris
Department of Agronomy, Kansas State University, Manhattan KS 66506
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Geoffrey P. Morris
  • For correspondence: gpmorris@ksu.edu

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Related Articles

Cited By

More in this TOC Section

  • A Statistical Procedure for Genome-Wide Detection of QTL Hotspots Using Public Databases with Application to Rice
  • GAL4 Drivers Specific for Type Ib and Type Is Motor Neurons in Drosophila
  • A Very Oil Yellow1 Modifier of the Oil Yellow1-N1989 Allele Uncovers a Cryptic Phenotypic Impact of Cis-regulatory Variation in Maize
Show more Investigations
  • Top
  • Article
    • Abstract
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgements
    • Footnotes
    • Literature Cited
  • Figures & Data
  • Supplemental
  • Info & Metrics

GSA

The Genetics Society of America (GSA), founded in 1931, is the professional membership organization for scientific researchers and educators in the field of genetics. Our members work to advance knowledge in the basic mechanisms of inheritance, from the molecular to the population level.

Online ISSN: 2160-1836

  • For Authors
  • For Reviewers
  • For Advertisers
  • Submit a Manuscript
  • Editorial Board
  • Press Releases

SPPA Logo

GET CONNECTED

RSS  Subscribe with RSS.

email  Subscribe via email. Sign up to receive alert notifications of new articles.

  • Facebook
  • Twitter
  • YouTube
  • LinkedIn
  • Google Plus

Copyright © 2019 by the Genetics Society of America

  • About G3
  • Terms of use
  • Permissions
  • Contact us
  • International access