- Open Access
A comprehensive study of the genomic differentiation between temperate Dent and Flint maize
© The Author(s). 2016
- Received: 16 December 2015
- Accepted: 15 June 2016
- Published: 8 July 2016
Dent and Flint represent two major germplasm pools exploited in maize breeding. Several traits differentiate the two pools, like cold tolerance, early vigor, and flowering time. A comparative investigation of their genomic architecture relevant for quantitative trait expression has not been reported so far. Understanding the genomic differences between germplasm pools may contribute to a better understanding of the complementarity in heterotic patterns exploited in hybrid breeding and of mechanisms involved in adaptation to different environments.
We perform whole-genome screens for signatures of selection specific to temperate Dent and Flint maize by comparing high-density genotyping data of 70 American and European Dent and 66 European Flint inbred lines. We find 2.2 % and 1.4 % of the genes are under selective pressure, respectively, and identify candidate genes associated with agronomic traits known to differ between the two pools. Taking flowering time as an example for the differentiation between Dent and Flint, we investigate candidate genes involved in the flowering network by phenotypic analyses in a Dent–Flint introgression library and find that the Flint haplotypes of the candidates promote earlier flowering. Within the flowering network, the majority of Flint candidates are associated with endogenous pathways in contrast to Dent candidate genes, which are mainly involved in response to environmental factors like light and photoperiod. The diversity patterns of the candidates in a unique panel of more than 900 individuals from 38 European landraces indicate a major contribution of landraces from France, Germany, and Spain to the candidate gene diversity of the Flint elite lines.
In this study, we report the investigation of pool-specific differences between temperate Dent and Flint on a genome-wide scale. The identified candidate genes represent a promising source for the functional investigation of pool-specific haplotypes in different genetic backgrounds and for the evaluation of their potential for future crop improvement like the adaptation to specific environments.
- Population genetics
- Genome-wide screen
Selection creates specific patterns of diversity in the genome  and these signatures can be used for the detection of regions under selection. When a favorable, new (derived compared to the ancestral) allele rises in frequency within a population, selective sweeps are generated, which are characterized by a local reduction in nucleotide diversity and high derived allele frequencies [9–11]. In addition, strong and recent sweeps will display large blocks with high linkage disequilibrium surrounding the derived mutation as the dispersal of the new allele will be faster than recombination is able to break down linkage disequilibrium [12, 13]. The identification of selection signatures through genome-wide screens provides an efficient way to detect selection candidates and methods for their detection are often combined to reduce the number of false-positives [14–16]. In maize, genome-wide screens for selection signatures were successfully applied to identify genes involved in domestication and improvement and allowed insights into evolutionary processes shaping the genome diversity of maize [17–20]. Taking advantage of the characteristics of selective sweeps and using high-density genotyping data from a maize 600 k single nucleotide polymorphism (SNP) array , we screened a panel of 136 temperate Dent and Flint elite lines for extreme allele frequencies over extended linked sites to identify genomic regions under selective pressure and to gain insights into the genomic variation underlying the differentiation of Dent and Flint. We included outgroup information from Sorghum bicolor to further support the identified candidate genes based on derived allele frequencies. We furthermore investigated the candidate genes based on whole-genome sequence data of 40 Dent and Flint lines [21, 22] and examined if genic and upstream regions contributed equally to the differentiation between temperate Dent and Flint.
The elite line panel under study comprised frequently used and important founder lines exploited in breeding programs for temperate climates. The Dent lines in our panel represent U.S. Corn Belt and European material, whereas most of the Flint lines originated from European breeding programs. Based on the selection screens, we examined pool-specific enrichment of candidate genes for metabolic pathways and investigated candidates associated with traits that are known to differentiate Dent and Flint like cold tolerance and flowering time [23, 24]. Flowering time is essential for local adaptation and represents a major determinant for other agronomic traits, such as grain filling and yield. The complex genetic architecture of flowering time has been studied in maize in a large number of studies mapping QTL with a meta-QTL analysis revealing 62 flowering time consensus QTL . Phenotypic differences in maize flowering time are mainly caused by the accumulation of many small-effect QTL  and only a few large-effect genes have been characterized so far [27–30]. Hundreds of homologs to A. thaliana flowering time genes have been found in the maize genome , but in most cases their functional roles in the maize flowering network remain to be elucidated [25, 26, 29, 32–37]. In this study, we identified candidate genes from the flowering network with haplotypes near fixation or fixed in either of the two elite pools. We used this set of genes as an example to characterize genomic differentiation between Dent and Flint in more detail. We evaluated the effect of these genes on flowering time in a Dent–Flint introgression library and investigated their assignment to different pathways within the flowering network. To assess the congruency of the allelic composition of the candidate genes between elite lines and landraces, we expanded our candidate gene analysis to a large dataset of 38 European landraces that comprises more than 900 individuals. By exploring this unique resource, we gained insights into the genetic variation of the selection candidates between landraces and elite lines and investigated, which landraces likely contributed to the observed candidate gene diversity in the elite lines and if haplotypes not yet exploited in breeding could be detected. Taken together, our study allowed insights into patterns of differentiation between temperate Dent and Flint germplasm and provided candidates for follow-up studies to characterize their biological and molecular functions, to investigate their impact on phenotypes, and to assess their potential use for further crop improvement.
Characterization of the Dent and Flint panels
We genotyped a diverse panel of 136 temperate inbred lines (Additional file 1: Table S1) at high density with the Axiom® Maize Genotyping Array . The array comprises more than 600 k SNP markers, which were identified based on mid- to high-coverage whole-genome sequence data of 30 representative temperate Dent and Flint maize lines . Markers were filtered according to quality scores and stable performance on the array, thus representing high-confidence sequence variants, and their final distribution followed the average recombination rate along the chromosomes . After stringent quality filtering of the 616,201 markers included on the array, 547,412 high-quality SNPs (88.8 %) remained for analysis. These SNPs tagged 19,759 genes (49.8 % of the annotated gene set of maize) with, on average, two SNPs in their coding region (52.6 % synonymous and 47.4 % non-synonymous). Slightly more SNPs were polymorphic in the Flint compared to the Dent panel (95.4 % versus 93.1 %), but the majority of SNPs segregated in both germplasm pools (88.6 %).
The panel of 136 temperate Dent and Flint inbred lines comprised frequently used and important founder lines exploited in breeding programs in Europe and the U.S., including lines which were used as parents for the U.S. and European NAM panels [6, 38, 39]. The 70 Dent lines were selected according to available pedigree information and their frequency of use and citation [40, 41] to assemble a representative set of lines. Besides 16 European Dent lines, the lines represent U.S. Corn Belt Dent and include lines from the Maize Association Population  and the list of inbred lines with expired U.S. plant variety protection . The 66 Flint lines investigated in this study comprised important founder lines of European breeding programs like F2 and F7 originating from the French landrace Lacaune, EP1 from the Spanish landrace Lizargarate, and derivatives of the German landrace Gelber Badischer Landmais . The Flints comprised in total 34 lines from France, 20 from Germany, four from Spain, three from Italy, three from North America, as well as one from Switzerland and Austria. Between the elite lines of the two germplasm pools, we observed a clear separation of pools (Fig. 1b) and a high genome-wide level of differentiation (F ST = 0.14), which is consistent with the long-term genetic differentiation between Dent-type and Flint-type maize [2, 3].
Genome-wide screens for selection signals
Taking advantage of the characteristics of selective sweeps, we screened the genome for extreme allele frequencies over extended linked sites to detect regions under differential selective pressure between Dent and Flint. Signatures of selection in only one of the two pools, Dent or Flint, were detected based on low levels of nucleotide diversity (π)  and Tajima’s D (TD)  in the respective pool. In addition, a signature had to be supported by a high value of the composite likelihood ratio (CLR) test  within the respective pool, which indicates a deviation of the allelic composition of a genetic region compared to a neutrally evolving sequence determined by the genomic background. To ensure that the selection signature was specific for one of the two pools, it had to be associated with a high level of differentiation between Dent and Flint measured by the fixation index F ST . Except for the CLR statistic, which was calculated for non-overlapping grids of 150 kb, we applied a sliding window approach averaging data over windows of 40 SNPs (sliding by 10 %) and filtered for regions below the 10 % quantile for π and TD and above the 90 % quantile for F ST and CLR (Additional file 1: Table S2). Following the approach reported by , adjacent windows passing the threshold for all four statistics were grouped together for candidate gene analysis, as the observed changes in allele frequency were likely caused by the same selective sweep event. This resulted in a filtered set of 265 windows for Dent and 158 windows for Flint, with an average length of 331.40 kb and 267.80 kb, respectively, and thus comparable to the length of domestication windows found in a previous study . An example of a signature of differential selection in Dent and Flint determined by all four metrics (π, TD, CLR, and F ST ) is shown in Fig. 1c for a region on chromosome 8 harboring two candidate genes. The underlying genetic region was composed of four major haplotypes. The first three haplotypes occurred at intermediate frequencies in Dent, whereas the fourth haplotype was almost exclusive for Flint.
Genome-wide patterns of diversity and the resulting distribution of selection signatures in the Dent and Flint panels are given in Additional file 2: Figure S1. Within the filtered set of windows, which covered 4.3 % of the total length of the maize genome for Dent and 2.1 % for Flint, we identified 876 genes as candidates under differential selective pressure in Dent and 545 genes for Flint with 14 genes common to both candidate genes sets (Additional file 3: Table S3). This corresponded to 2.2 % and 1.4 % of the filtered gene set of maize, respectively, and is in the same order of magnitude as the estimated number of genes under selective pressure during maize domestication and improvement . When comparing the candidate gene sets with the 571 improvement candidates reported by , 26 genes overlapped with the list of Dent candidates but only one gene with the Flint candidate gene set. Considering that the genetic material studied in  comprised mainly U.S. Dent and (sub-) tropical lines and that pool-specific sequence variation in temperate Dent and Flint has been reported here and, for example, by , these results emphasize the relevance of a representative panel of lines belonging to divergent germplasm pools to obtain a comprehensive picture of the genomic diversity in maize.
In genome-wide screens for signatures of positive selection, also other forces than selection, such as heterogeneous mutation and recombination rates along the genome, past demographic history and background selection shape the genomic diversity and can give rise to false-positive signals. It is beyond the scope of this paper to infer a full demographic history of maize for the elite lines and landraces as the breeding history of maize is complex and violates several assumptions of the classic population genetics models (e.g. discussed in ), as, for example, the assumption of panmictic populations and applicability of the coalescent at short time scales. We therefore applied the CLR test , which detects selective sweeps based on the comparison of the site-frequency spectrum within a specific genomic region to the average site-frequency spectrum over the genome, a method which has been successfully used in human and other species to detect selective sweeps [11, 47, 48]. To further decrease the rate of false-positives, the CLR test was combined with three additional metrics (π, TD, and F ST ) and we identified signatures of positive selection based on this conservative approach with an overlap of genome-wide extreme values per metric. The high level of linkage disequilibrium in temperate Dent and Flint elite lines  facilitates the detection of selective sweep signals over sufficiently large genomic regions by the CLR test. On the other hand, the extent of linkage disequilibrium may decrease the power to discriminate between signals caused by genetic hitchhiking due to positive selection and negative background selection in regions with reduced levels of recombination [49, 50]. To assess the number of false-positives due to this effect, we explored the recombination landscape in the Dent and Flint panels by estimating lower bounds of historical recombination events . The proportion of candidate genes located in regions with strongly reduced recombination rates and high linkage disequilibrium like (peri-) centromeric regions was then estimated. We found that 74.8 % of the Dent and 80.9 % of the Flint candidates were not located in regions with low levels of recombination (10 % quantile per chromosome; Additional file 1: Figure S2) indicating that the majority of candidates represent targets of selection rather than false-positive signals. Furthermore, in a classic selective sweep scenario (in contrast to background selection) targets of selection are to be enriched for derived alleles. As an additional test of our candidate regions, we included information from Sorghum bicolor to distinguish between ancestral and derived alleles. The Dent and Flint candidate gene sets revealed significantly higher derived allele frequencies compared to the remaining genes as measured by Fay and Wu’s normalized H  (p < 2.2e-16; Additional file 1: Table S4), which also supported positive selection as the driving force of the observed allele frequency changes.
Gene ontology and pathway analyses of candidate gene sets
Considering genetic differentiation and distinct phenotypic characteristics of Dent and Flint, we tested whether the candidate gene sets were enriched for specific biological processes or pathways. Gene ontology (GO) terms associated with the identified genes were available for around 40 % of the candidates (333 for Dent and 214 for Flint). No significant GO term enrichment of biological processes, cellular components, and molecular functions could be detected for either of the two sets (Additional file 1: Figure S3). To investigate if candidate genes revealed a pool-specific enrichment for metabolic pathways, we performed pathway analyses using MapMan . Based on information available for 58 Dent and 40 Flint candidate genes, we observed a grouping of genes associated with tetrapyrroles (chlorophyll and heme precursors) for Dent and for terpenoid metabolism for Flint (Additional file 1: Figure S4). The latter included the two genes ZmPPS7.3 (GRMZM2G014508) and ZmPPS8.2 (GRMZM2G483889), which encode a large and a small subunit of the geranyl diphosphate synthase complex in maize, respectively . Like their homologues in A. thaliana , they are assumed to be involved in the biosynthesis of precursors of hormones from the isoprenoid pathway (e.g. gibberellins, brassinosteroids, and abscisic acid). The ability to produce other downstream products of this enzyme, namely β-caryophyllenes, has been shown to differ between European Flint and U.S. Dent lines and suggested that this defense response signal against herbivores was largely lost in temperate U.S. Dent [56, 57]. The analysis of candidates associated with other traits that are known to differentiate Dent and Flint revealed six Flint candidates that, according to GO terms, are related to cold tolerance, a trait that is characteristic for temperate Flint . For two of the candidates, differential expression upon exposure to chilling temperature has been reported in maize (GRMZM2G035584  and GRMZM2G095562 ) as well as for the homologous gene of GRMZM2G139680 in rice . The molecular and functional characterization of the identified candidate genes in maize and the investigation of differences between Dent and Flint in the regulation of phytohormone pathways or secondary metabolism may provide further insights in the adaptation of maize to different environments. Up to now, comprehensive RNA expression data across various developmental stages and tissues are mainly available for U.S. Dent lines like B73, which underlines the need for a better structural and functional genomic characterization of the Flint germplasm pool and its unique properties.
Assessing the phenotypic effects of candidate genes on flowering time in a Dent–Flint introgression library
In the genome-wide selection screens, we identified 18 candidates for Dent and 12 candidates for Flint, which could be assigned to the flowering pathway based on previous reports in maize, GO terms, and/or sequence homology to flowering genes characterized in other species [30, 32–34, 61, 62]. We focused exemplarily on candidate genes associated with the flowering network in maize as flowering time is an important agronomic trait that differentiates temperate Dent and Flint. However, functional studies of these genes in maize were available for only 30 % of the candidates (Additional file 4: Table S5). Here, we investigated the effect of the flowering time candidate genes in more detail using a maize introgression library.
Of the 97 lines, 22 carried a Flint introgression harboring one or several of the flowering time candidates identified in the selection screens (Additional file 5: Table S6). Fourteen of the 30 candidates were represented in these 22 introgression lines. Seven lines carried a segment with one or more of seven flowering time candidates identified in Flint and nine lines carried one or more of six flowering time candidates identified in Dent. Six lines carried a segment with a combination of Dent and Flint candidates. Although 75 lines did not carry one of the flowering time candidates identified in our selection screens, they may carry other flowering time genes with alleles differing between the Dent and the Flint parent of the introgression library. Lines carrying the Flint haplotype of a Flint candidate differed significantly from the 75 lines which did not carry one of the flowering time candidates from the selection screens (93.1 versus 96.1 days, p value = 0.011; Fig. 2a). For the lines which carried the Flint haplotype of a Dent selection candidate, this difference was not significant. The results indicate that in the genetic material under study, the Flint haplotypes of Flint candidates promoted flowering time more than the Flint haplotypes of Dent candidates.
Of the six lines with significant difference in male flowering compared to the Dent parent after correcting for multiple testing (α = 0.05/97), two carried Flint haplotypes of Flint candidates and one a Flint haplotype of a Dent candidate. One of the lines included the well-characterized large-effect region comprising the ethylene-responsive transcription factor Rap2 (related to APETALA2 7, ZmRap2.7, Rap2, GRMZM2G700665) and its regulatory upstream locus Vgt1 , a major QTL for flowering time in maize [29, 30]. The other line contained Zcn1 (one of several members of the ZEA CENTRORADIALIS or TERMINAL FLOWER1 (TFL1)-like gene family ; also Phosphatidylethanolamine-binding protein1, Pebp1, GRMZM2G092008), for which so far only a moderate effect on flowering time was reported in maize . This gene is related to TFL1 in A. thaliana , which is an antagonist of the FLOWERING LOCUS T (FT) [64, 65] and required for the maintenance of an indeterminate inflorescence meristem identity and the regulation of flowering time in A. thaliana and maize [63–66]. The line with the Dent candidate carried Zmm22 (MADS-transcription factor 69, Mads69, GRMZM2G171650), which was recently reported to be associated with variation in flowering time in maize  and is considered a candidate for maize domestication and/or improvement [68, 69].
Overall, our findings in the introgression library support the relevance of the investigated genomic regions and their associated candidates for promoting flowering time and confirm the quantitative nature of flowering time in maize, determined by many genes with small effects  and only few genes with larger effects. Here, the effects of Zcn1 and Zmm22 were stronger than reported previously, which may be attributed to a stronger substitution effect when replacing a Dent haplotype with a Flint haplotype. We will target potential expression differences of flowering time candidates in the Dent–Flint introgression library in future studies to characterize possible differences in the regulation of the flowering network between germplasm pools adapted to different environments.
Differential selection on components of the flowering network within temperate maize
Diversity of flowering time candidates in elite lines and European landraces
Consistent with the hypothesis that Flint elite lines and Flint-type landraces have a common history, significantly higher levels of differentiation were observed for Dent candidate genes compared to all remaining genes (F ST = 0.111 versus 0.095, p value = 0.017; Fig. 4b, A versus C). Together, these findings indicated that the reduced diversity observed for Flint candidate genes in Flint elite lines was already present in a broad panel of European landraces and that the candidate gene diversity of the Flint elite lines originate from a limited number of Flint-type landraces used for elite line development in some historically important breeding centers .
The remaining seven landraces displayed at least partially Dent-type kernels. These landraces revealed high levels of diversity for Dent and Flint flowering time candidates (π = 0.225 and 0.260, respectively; Fig. 3b, Additional file 6: Table S8) and showed a high level of differentiation with Dent elite lines for the Dent flowering time candidates (F ST = 0.170; Fig. 4c, FT). The same pattern was found in the analysis of the entire Dent candidate gene set, which revealed significantly higher levels of differentiation compared to all remaining genes (F ST = 0.164 versus 0.111, p value = 0.026; Fig. 4c, A versus C), but no significant difference for Flint candidates compared to all remaining genes (F ST = 0.138 versus 0.112, p value = 0.209; Fig. 4d, A versus C; Additional file 7: Table S9). These results indicated that the European Dent-type landraces exhibit a different allelic composition in the Dent candidates compared to the Dent elite lines and did most likely not contribute to the Dent elite material under study.
Selection on upstream and genic regions of the candidates
To examine how specific elements of the genic regions contributed to the differentiation between Dent and Flint, we compared levels of differentiation for 5 kb and 500 bp upstream regions, genic regions, and exons between the candidate gene sets and all remaining genes. To increase the resolution of our analyses, we investigated the candidate gene sets based on whole-genome sequence data of 40 temperate elite lines (21 Dent and 19 Flint) [21, 22], which were part of the panel of 136 elite lines genotyped with the 600 k array with the exception of three lines (Additional file 1: Table S1). Based on 13,246,294 bi-allelic SNPs, we observed a significant reduction of mean π and TD in 727 Dent and 403 Flint candidate genes tagged by at least five SNPs (of in total 876 and 545 candidates, respectively) compared to 31,163 remaining genes (p value < 2.2e-16; Additional file 1: Figure S5 and Additional file 1: Table S4). F ST values calculated between Dent and Flint were significantly higher for candidate gene sets compared to all remaining genes for 5 kb and 500 bp upstream as well as genic and exonic regions. Together, these findings supported the results obtained from the selection screens in the panel of 136 temperate inbred lines genotyped with the 600 k array.
Previous studies in maize suggested an important role of the divergence of regulatory elements in the context of domestication [73–75]. In our study, distributions of F ST values were comparable for 5 kb and 500 bp upstream as well as genic and exonic regions in each of the two candidate gene sets (Additional file 1: Figure S6 and Additional file 1: Table S4). However, the power to resolve whether selection acted differentially in upstream and genic regions was probably limited by the high level of linkage disequilibrium observed in temperate Dent and Flint lines . The outcome of ongoing large-scale whole genome and transcriptome sequencing will allow the investigation of the impact of selection on the regulation of gene activity in the two pools and its consequence for the genomic differentiation between Dent and Flint.
In this study, we report genomic differentiation between two major temperate maize germplasm pools, Dent and Flint. By comparing a representative panel of Dent and Flint elite lines, we identified candidate genes under differential selective pressure in Dent and Flint. The significant enrichment in derived allele frequencies for these genes provided strong indication that the candidate regions represented selective sweeps. Candidate genes associated with agronomic traits known to differ between Dent and Flint could be identified. Most of the detected flowering time candidates have not yet been functionally characterized in maize. Investigating the effect of the flowering time candidates in a Dent–Flint introgression library, we found that Flint haplotypes of these candidates promoted earlier flowering. Within the flowering network of maize, a Flint-specific enrichment of genes associated with endogenous signaling, hormone-dependent pathways, and developmental processes was discovered in contrast to Dent, where selection seemed to act predominantly on genes involved in the response to environmental factors. Low levels of differentiation of Flint flowering time candidate genes between European Flint elite lines and European landraces indicated that the allelic composition of the elite lines was comparable to those of the Flint-type landraces and suggested a major contribution of landraces from France, Germany, and Spain to the candidate gene diversity in the Flint elite lines. Our findings highlight the role of genomic regions that have undergone intense selection and contributed to the differentiation of temperate Dent and Flint with likely effects on different agronomic traits. The identification of pool-specific selection signatures enabled insights into different patterns of diversity between temperate Dent and Flint and provides new targets for future functional analyses and crop improvement.
Plant material and genotyping of elite lines and landraces
The 136 elite inbred lines (Additional file 1: Table S1) were selected to represent the genetic diversity of European and American temperate maize and were genotyped with the 600 k Affymetrix® Axiom® Maize Array described in . Landraces were selected to reflect the genetic and phenotypic diversity within Central and Western Europe (Additional file 1: Table S7) and were represented by 22 to 24 plants each. All 906 landrace individuals were genotyped using the 600 k array and the genotype cluster model file (http://www.affymetrix.com/catalog/prod820010/AFFY/Axiom%26%23174%3B-Maize-Genotyping-Array).
If not denoted otherwise, analyses were performed using R version 3.0.1 . SNP positions were assigned to the reference sequence B73 v2  for all datasets. Analyses of the elite line panel were based on 566,961 best quality SNPs  with heterozygous calls masked as missing. Indels, unmapped markers, SNPs with ≥ 10 % missing values in the elite line panel , markers designed to specifically differentiate between two Dent lines , and monomorphic SNPs were excluded, resulting in 547,412 SNPs for analyses of temperate elite lines. Analyses including the landraces were based on a subset of 486,208 SNPs.
Population structure analysis and estimation of historical recombination rates
For population structure analyses of the genotyped panel of 136 elite lines, missing data were imputed using Beagle  version 3.3.1 via the R package “synbreed”  version 0.10-3. A linkage disequilibrium pruning step was performed applying an r 2 threshold of 0.8 followed by the estimation of ancestry using ADMIXTURE  version 1.23. To obtain an estimate of the historical recombination rates in Dent and Flint based on the genotype data of the elite line panel, the four-gamete test  was calculated. This test gives a conservative estimate (i.e. the minimum number) of recombination events in the history of a sample. Values for the pairwise tests of neighboring SNPs were averaged over 1000 sites and the mean of each 1000 site bin was plotted. Reported were recombination events per Mb and regions of low recombination rates were defined as regions exhibiting rates within the 10 % quantile per chromosome.
Screens for selection signatures
Nucleotide diversity π  and Tajima’s D (TD)  were calculated for each panel of inbred lines using a customized script. The fixation index F ST  was calculated across the two panels using PLINK  version v1.90b2m. Metrics were calculated per SNP and averaged over windows of 40 SNPs (sliding by 10 % and corresponding to an average physical distance of 10 kb), using the R package “zoo” . The grid-based composite likelihood ratio (CLR) test was calculated for each panel as implemented in SweepFinder . For CLR, the size of the non-overlapping grids was 150 kb, which is the same magnitude as the maximal distance between two SNPs with r 2 > 0.2 in the elite line panel . Windows exhibiting values within the 10 % quantile (π, TD) and the 90 % quantile (F ST , CLR) were submitted for candidate gene analysis based on the B73 v2  annotation, version 5b60 (ftp://ftp.gramene.org/pub/gramene/maizesequence.org/release-5b/filtered-set/), containing 39,656 gene models. GO terms were tested for enrichment using the resources developed by  based on maize gene IDs by applying a hypergeometric test with a Benjamini–Yekutieli correction (; FDR = 0.05) to account for multiple tests. Pathway analysis was performed using MapMan version 3.5.1  based on the mapping of the first transcript of each gene to the file Zm_B73_5b_FGS_cds_2012.m02 downloaded from the MapMan webpage (http://mapman.gabipd.org/web/guest/mapmanstore). For visualization, arbitrary values of –4.5 and 4.5 were assigned to Dent and Flint candidate gene transcripts, respectively.
Whole-genome sequence datasets
For Fay and Wu’s normalized H , the maize-sorghum genome alignment (http://pipeline.lbl.gov/downloads.shtml) was parsed with a custom Perl script to obtain the nucleotide in sorghum representing the ancestral maize allele for 298,388 SNPs of the genotyped panel of elite lines. Whole-genome sequence data for 30 elite lines were used to call SNPs and small indels by employing an integrative analysis pipeline  and filtering for high mapping (MQ ≥ 30) and genotyping quality (GQ > 5), major allele frequency ≥ 90 %, a minimal distance of 3 bp between adjacent SNPs, and a minimal and maximal coverage by three and 60 MQ30 reads for lines sequenced to medium coverage, respectively. The latter criterion was adjusted for four deep sequenced lines requiring ten and 300 MQ30 reads, respectively. We combined the obtained marker set with whole-genome sequence data from ten temperate inbred Dent and Flint lines from the maize HapMap2 project [17, 22] resulting in a VCF file including 13,246,294 bi-allelic SNPs. We filtered for SNPs with ≤ 50 % missing values across the 40 lines (F ST ) and within germplasm pools (π and TD). The combined VCF file was converted to hapmap format with a customized Perl script, and π and TD per gene were obtained with Variscan version 2  with runmode “12”.
For gene-wise calculations based on the genotyped panels of 136 elite lines and 38 European landraces (normalized H, F ST ) or on whole-genome sequence data of 40 elite lines (π and TD), the genic region of the longest protein-coding transcript, including 5 kb upstream, was used and metrics were calculated if at least 5 SNPs were available for analysis. For a separate analysis of exonic, genic, 500 bp, and 5 kb upstream regions, F ST was determined in case of at least 5 SNPs per region based on whole-genome sequence data using vcftools v0.1.11 . Two-sided Wilcoxon rank sum tests  were performed on gene-wise metrics to test for differences between candidate genes and remaining genes within pools.
To represent the Flint genome as introgression segments in a Dent genetic background, a European Dent inbred line was crossed with a European Flint inbred line, followed by backcrossing, marker-assisted selection, and several rounds of selfing . The Dent parental line originated from south-eastern Europe and exhibits a high general combining ability for kernel yield. The Flint parental line was selected for high general combining ability for biomass yield and good performance in Central European climates. The introgression library used in the present study comprised 535 lines, with 97 lines carrying a single segment of the Flint parent. To estimate the length of individual donor genome fragments, the distance between markers on the respective donor genome fragment plus half the distance to the adjacent marker flanking the donor genome fragment on either side of the fragment was calculated. For the introgression lines, the two parental lines, and a check, male and female flowering times were obtained from two field experiments conducted in 2014 at the German trial locations Roggenstein (N 48°11′13.24″, E 11°19′50.86″, 517 m AMSL, average temperature in 2014: 9.8 °C) and Freising (N 48°24′11.62″, E 11°43′21.99″, 480 m AMSL, average temperature in 2014: 9.7 °C). Each experiment was laid out as an α–lattice design with two replications, except for parental lines and the check that were repeated three and five times, respectively. Male and female flowering time was recorded as days after sowing until 50 % of plants per plot exhibited emerged anthers and silks, respectively. Adjusted means for flowering time were calculated using Plabstat . The difference in adjusted means between the 97 single-segment introgression lines and the Dent parental line was tested at a significance level of α = 0.05 and α = 0.05/97 to correct for multiple testing. Adjusted means of flowering times for seven lines carrying Flint haplotypes of Flint flowering time candidates, nine lines carrying Flint haplotypes of Dent flowering time candidates, and the 75 lines not carrying a flowering time candidate gene were tested for significant differences by calculating Student’s t-test.
We are grateful to Angel Alvarez (CSIC, Estación Experimental de Aula Dei, Zaragoza, Spain), Alain Charcosset (INRA, UMR de Génétique Quantitative et Evolution, Gif-Sur-Yvette, France), Barbara Eder and Joachim Eder (Bavarian State Research Center, Institute for Crop Science and Plant Breeding, Freising, Germany), Albrecht E. Melchinger (University of Hohenheim, Institute of Plant Breeding, Seed Science and Population Genetics, Stuttgart, Germany), Mark J. Millard (USDA-ARS, Ames, USA), Jesús Moreno González (CIAM, A Coruña, Spain), Rita Redaelli (CRA-MAC Maize Research Unit, Bergamo, Italy), Pedro Revilla (CSIC, Pontevedra, Spain), and Anne Zanetto (INRA, UE de Melgueil DiaScope, Mauguio, France) for providing elite line and/or landrace seeds. The original introgression library was kindly provided by Milena Ouzunova and Thomas Presterl (KWS SAAT SE, Einbeck, Germany). We thank Manfred Schönleben for his help with field trials and their analysis and Stefan Schwertfirm for his technical assistance, as well as Ruedi Fries (Technical University of Munich, Animal Breeding, Freising, Germany) for providing genotyping facilities. Furthermore, we thank our colleagues from the University of Hohenheim for kindly providing pictures of maize cobs and Joachim Hermisson (University of Vienna, Mathematics and BioSciences Group, Vienna, Austria) for fruitful discussions and helpful remarks during a research stay of SU in Vienna. We thank three anonymous reviewers for their valuable comments and suggestions.
This study was funded by the Federal Ministry of Education and Research (BMBF, Germany) within the AgroClustEr Synbreed - Synergistic plant and animal breeding (grant 0315528) and by the Spanish Ministry of Economy and Competitiveness (project RF2011-00022-C02-01 “Regeneration and rationalization of the maize landraces from the Iberian Peninsula”). Bernardo Ordas acknowledges a grant from the program “Ramón y Cajal” of the Spanish Ministry of Economy and Competitiveness.
Availability of data and materials
Genotype calls for elite lines and landraces are available at NCBI dbSNP upon Build B149 under accession number ss1998449529 to ss1999065412 and ss1999065413 to ss1999681296, respectively. Whole-genome sequence data for 30 elite lines are available under accession number PRJNA260788 at the NCBI BioProject webpage and the combined and filtered VCF file of whole-genome sequence data for 30 elite lines and the ten lines from the maize HapMap2 project [17, 22] are available at NCBI dbSNP under accession number ss1999681297 to ss2012927636 (upon Build B149). All three datasets are available at https://0-dx-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.3427040.v1.
CCS, EB, AT, and SU conceived the study and discussed the results; SU investigated genotypic and outgroup data; SU and MM analyzed genotypic data of landraces; MAS, GH, and KFXM carried out variant calling; SP extracted ancestral allele information; RP, SP, and SU performed sequence analyses; PW conducted the introgression library experiments; SU and EB investigated the introgression library data; BO contributed part of the Spanish landrace data; HP performed genotyping of elite lines and landraces; SU drafted the manuscript; EB, CCS, and AT edited the manuscript; all authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Open AccessThis article is 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 you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lobell DB, Tebaldi C. Getting caught with our plants down: the risks of a global crop yield slowdown from climate trends in the next two decades. Environ Res Lett. 2014;9:074003.View ArticleGoogle Scholar
- Smith CW, Betrán J, Runge ECA. Corn: origin, history, technology, and production. 1st ed. Hoboken, NJ: John Wiley & Sons Inc.; 2004.Google Scholar
- Brown WL, Anderson E. The northern flint corn. Ann Mo Bot Gard. 1947;34:1–28.View ArticleGoogle Scholar
- Rebourg C, Chastanet M, Gouesnard B, Welcker C, Dubreuil P, Charcosset A. Maize introduction into Europe: the history reviewed in the light of molecular data. Theor Appl Genet. 2003;106:895–903.PubMedGoogle Scholar
- Dubreuil P, Dufour P, Krejci E, Causse M, deVienne D, Gallais A, et al. Organization of RFLP diversity among inbred lines of maize representing the most significant heterotic groups. Crop Sci. 1996;36:790–9.View ArticleGoogle Scholar
- Bauer E, Falque M, Walter H, Bauland C, Camisan C, Campo L, et al. Intraspecific variation of recombination rate in maize. Genome Biol. 2013;14:R103.View ArticlePubMedPubMed CentralGoogle Scholar
- Giraud H, Lehermeier C, Bauer E, Falque M, Segura V, Bauland C, et al. Linkage disequilibrium with linkage analysis of multiline crosses reveals different multiallelic QTL for hybrid performance in the flint and dent heterotic groups of maize. Genetics. 2014;198:1717–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Nielsen R. Molecular signatures of natural selection. Annu Rev Genet. 2005;39:197–218.View ArticlePubMedGoogle Scholar
- Tajima F. Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983;105:437–60.PubMedPubMed CentralGoogle Scholar
- Tajima F. Statistical-method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMedPubMed CentralGoogle Scholar
- Nielsen R, Williamson S, Kim Y, Hubisz MJ, Clark AG, Bustamante C. Genomic scans for selective sweeps using SNP data. Genome Res. 2005;15:1566–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.View ArticlePubMedGoogle Scholar
- Voight BF, Kudaravalli S, Wen XQ, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:446–58.View ArticleGoogle Scholar
- Long Q, Rabanal FA, Meng D, Huber CD, Farlow A, Platzer A, et al. Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden. Nat Genet. 2013;45:884–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, et al. Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009;19:826–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Qanbari S, Gianola D, Hayes B, Schenkel F, Miller S, Moore S, et al. Application of site and haplotype-frequency based approaches for detecting selection signatures in cattle. BMC Genomics. 2011;12:318.View ArticlePubMedPubMed CentralGoogle Scholar
- Hufford MB, Xu X, van Heerwaarden J, Pyhäjärvi T, Chia JM, Cartwright RA, et al. Comparative population genomics of maize domestication and improvement. Nat Genet. 2012;44:808–U118.View ArticlePubMedGoogle Scholar
- Jiao YP, Zhao HN, Ren LH, Song WB, Zeng B, Guo JJ, et al. Genome-wide genetic changes during modern breeding of maize. Nat Genet. 2012;44:812–U124.View ArticlePubMedGoogle Scholar
- Vigouroux Y, McMullen M, Hittinger CT, Houchins K, Schulz L, Kresovich S, et al. Identifying genes of agronomic importance in maize by screening microsatellites for evidence of selection during domestication. Proc Natl Acad Sci U S A. 2002;99:9650–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamasaki M, Tenaillon MI, Bi IV, Schroeder SG, Sanchez-Villeda H, Doebley JF, et al. A large-scale screen for artificial selection in maize identifies candidate agronomic loci for domestication and crop improvement. Plant Cell. 2005;17:2859–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Unterseer S, Bauer E, Haberer G, Seidel M, Knaak C, Ouzunova M, et al. A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array. BMC Genomics. 2014;15:823.View ArticlePubMedPubMed CentralGoogle Scholar
- Chia JM, Song C, Bradbury PJ, Costich D, de Leon N, Doebley J, et al. Maize HapMap2 identifies extant variation from a genome in flux. Nat Genet. 2012;44:803–7.View ArticlePubMedGoogle Scholar
- Camus-Kulandaivelu L, Veyrieras JB, Madur D, Combes V, Fourmann M, Barraud S, et al. Maize adaptation to temperate climate: relationship between population structure and polymorphism in the Dwarf8 gene. Genetics. 2006;172:2449–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Revilla P, Rodríguez VM, Ordás A, Rincent R, Charcosset A, Giauffret C, et al. Cold tolerance in two large maize inbred panels adapted to European climates. Crop Sci. 2014;54:1981–91.View ArticleGoogle Scholar
- Chardon F, Virlon B, Moreau L, Falque M, Joets J, Decousset L, et al. Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome. Genetics. 2004;168:2169–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Buckler ES, Holland JB, Bradbury PJ, Acharya CB, Brown PJ, Browne C, et al. The genetic architecture of maize flowering time. Science. 2009;325:714–8.View ArticlePubMedGoogle Scholar
- Colasanti J, Yuan Z, Sundaresan V. The indeterminate gene encodes a zinc finger protein and regulates a leaf-generated signal required for the transition to flowering in maize. Cell. 1998;93:593–603.View ArticlePubMedGoogle Scholar
- Muszynski MG, Dam T, Li B, Shirbroun DM, Hou Z, Bruggemann E, et al. Delayed flowering1 encodes a basic leucine zipper protein that mediates floral inductive signals at the shoot apex in maize. Plant Physiol. 2006;142:1523–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Salvi S, Tuberosa R, Chiapparino E, Maccaferri M, Veillet S, van Beuningen L, et al. Toward positional cloning of Vgt1, a QTL controlling the transition from the vegetative to the reproductive phase in maize. Plant Mol Biol. 2002;48:601–13.View ArticlePubMedGoogle Scholar
- Vladutu C, McLaughlin J, Phillips RL. Fine mapping and characterization of linked quantitative trait loci involved in the transition of the maize apical meristem from vegetative to generative structures. Genetics. 1999;153:993–1007.PubMedPubMed CentralGoogle Scholar
- Chen C, DeClerck G, Tian F, Spooner W, McCouch S, Buckler E. PICARA, an analytical pipeline providing probabilistic inference about a priori candidates genes underlying genome-wide association QTL in plants. PLoS One. 2012;7:e46596.View ArticlePubMedPubMed CentralGoogle Scholar
- Danilevskaya ON, Meng X, Hou Z, Ananiev EV, Simmons CR. A genomic and expression compendium of the expanded PEBP gene family from maize. Plant Physiol. 2008;146:250–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Dong Z, Danilevskaya O, Abadie T, Messina C, Coles N, Cooper M. A gene regulatory network model for floral transition of the shoot apex in maize and its dynamic modeling. PLoS One. 2012;7:e43450.View ArticlePubMedPubMed CentralGoogle Scholar
- Ducrocq S, Madur D, Veyrieras JB, Camus-Kulandaivelu L, Kloiber-Maitz M, Presterl T, et al. Key impact of Vgt1 on flowering time adaptation in maize: evidence from association mapping and ecogeographical information. Genetics. 2008;178:2433–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu J, Liu Y, Liu J, Cao M, Wang J, Lan H, et al. The genetic architecture of flowering time and photoperiod sensitivity in maize as revealed by QTL review and meta analysis. J Integr Plant Biol. 2012;54:358–373.View ArticlePubMedGoogle Scholar
- Romay MC, Millard MJ, Glaubitz JC, Peiffer JA, Swarts KL, Casstevens TM, et al. Comprehensive genotyping of the USA national maize inbred seed bank. Genome Biol. 2013;14:R55.View ArticlePubMedPubMed CentralGoogle Scholar
- Castelletti S, Tuberosa R, Pindo M, Salvi S. A MITE transposon insertion is associated with differential methylation at the maize flowering time QTL Vgt1. G3 (Bethesda). 2014;4:805–12.View ArticleGoogle Scholar
- McMullen MD, Kresovich S, Villeda HS, Bradbury P, Li H, Sun Q, et al. Genetic properties of the maize nested association mapping population. Science. 2009;325:737–40.View ArticlePubMedGoogle Scholar
- Lehermeier C, Krämer N, Bauer E, Bauland C, Camisan C, Campo L, et al. Usefulness of multiparental populations of maize (Zea mays L.) for genome-based prediction. Genetics. 2014;198:3–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Mikel MA. Availability and analysis of proprietary dent corn inbred lines with expired US plant variety protection. Crop Sci. 2006;46:2555–60.View ArticleGoogle Scholar
- Mikel MA, Dudley JW. Evolution of North American dent corn from public to proprietary germplasm. Crop Sci. 2006;46:1193–205.View ArticleGoogle Scholar
- Flint-Garcia SA, Thuillet AC, Yu J, Pressoir G, Romero SM, Mitchell SE, et al. Maize association population: a high-resolution platform for quantitative trait locus dissection. Plant J. 2005;44:1054–64.View ArticlePubMedGoogle Scholar
- Nelson PT, Coles ND, Holland JB, Bubeck DM, Smith S, Goodman MM. Molecular characterization of maize inbreds with expired US plant variety protection. Crop Sci. 2008;48:1673–85.View ArticleGoogle Scholar
- Barrière Y, Alber D, Dolstra O, Lapierre C, Motto M, Ordas A, et al. Past and prospects of forage maize breeding in Europe: II. History, germplasm evolution and correlative agronomic changes. Maydica. 2006;51:435–49.Google Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.View ArticleGoogle Scholar
- Parat F, Schwertfirm G, Rudolph U, Miedaner T, Korzun V, Bauer E, et al. Geography and end use drive the diversification of worldwide winter rye populations. Mol Ecol. 2016;25:500–14.View ArticlePubMedGoogle Scholar
- Qanbari S, Pausch H, Jansen S, Somel M, Strom TM, Fries R, et al. Classic selective sweeps revealed by massive sequencing in cattle. PLoS Genet. 2014;10:e1004148.View ArticlePubMedPubMed CentralGoogle Scholar
- Horton MW, Hancock AM, Huang YS, Toomajian C, Atwell S, Auton A, et al. Genome-wide patterns of genetic variation in worldwide Arabidopsis thaliana accessions from the RegMap panel. Nat Genet. 2012;44:212–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Stephan W. Genetic hitchhiking versus background selection: the controversy and its implications. Philos Trans R Soc Lond B Biol Sci. 2010;365:1245–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Charlesworth B, Morgan MT, Charlesworth D. The effect of deleterious mutations on neutral molecular variation. Genetics. 1993;134:1289–303.PubMedPubMed CentralGoogle Scholar
- Hudson RR, Kaplan NL. Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics. 1985;111:147–64.PubMedPubMed CentralGoogle Scholar
- Zeng K, Fu YX, Shi SH, Wu CI. Statistical tests for detecting positive selection by utilizing high-frequency variants. Genetics. 2006;174:1431–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37:914–39.View ArticlePubMedGoogle Scholar
- Zhou M, Zhang Q, Wang C, Chen L, Sun Z, Zhu X, et al. Characterization of genes involved in isoprenoid diphosphate biosynthesis in maize. J Plant Growth Regul. 2015;34:294–308.View ArticleGoogle Scholar
- van Schie CC, Ament K, Schmidt A, Lange T, Haring MA, Schuurink RC. Geranyl diphosphate synthase is required for biosynthesis of gibberellins. Plant J. 2007;52:752–62.View ArticlePubMedGoogle Scholar
- Degen T, Dillmann C, Marion-Poll F, Turlings TC. High genetic variability of herbivore-induced volatile emission within a broad range of maize inbred lines. Plant Physiol. 2004;135:1928–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Kollner TG, Held M, Lenk C, Hiltpold I, Turlings TC, Gershenzon J, et al. A maize (E)-beta-caryophyllene synthase implicated in indirect defense responses against herbivores is not expressed in most American maize varieties. Plant Cell. 2008;20:482–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Trzcinska-Danielewicz J, Bilska A, Fronk J, Zielenkiewicz P, Jarochowska E, Roszczyk M, et al. Global analysis of gene expression in maize leaves treated with low temperature I. Moderate chilling (14 °C). Plant Sci. 2009;177:648–58.View ArticleGoogle Scholar
- Sobkowiak A, Jonczyk M, Jarochowska E, Biecek P, Trzcinska-Danielewicz J, Leipner J, et al. Genome-wide transcriptomic analysis of response to low temperature reveals candidate genes determining divergent cold-sensitivity of maize inbred lines. Plant Mol Biol. 2014;85:317–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Yan SP, Zhang QY, Tang ZC, Su WA, Sun WN. Comparative proteomic analysis provides new insights into chilling stress responses in rice. Mol Cell Proteomics. 2006;5:484–96.View ArticlePubMedGoogle Scholar
- Hung HY, Shannon LM, Tian F, Bradbury PJ, Chen C, Flint-Garcia SA, et al. ZmCCT and the genetic basis of day-length adaptation underlying the postdomestication spread of maize. Proc Natl Acad Sci U S A. 2012;109:E1913–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Salvi S, Sponza G, Morgante M, Tomes D, Niu X, Fengler KA, et al. Conserved noncoding genomic sequences associated with a flowering-time quantitative trait locus in maize. Proc Natl Acad Sci U S A. 2007;104:11376–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Danilevskaya ON, Meng X, Ananiev EV. Concerted modification of flowering time and inflorescence architecture by ectopic expression of TFL1-like genes in maize. Plant Physiol. 2010;153:238–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Kobayashi Y, Kaya H, Goto K, Iwabuchi M, Araki T. A pair of related genes with antagonistic roles in mediating flowering signals. Science. 1999;286:1960–2.View ArticlePubMedGoogle Scholar
- Koornneef M, Hanhart CJ, van der Veen JH. A genetic and physiological analysis of late flowering mutants in Arabidopsis thaliana. Mol Gen Genet. 1991;229:57–66.View ArticlePubMedGoogle Scholar
- Hanano S, Goto K. Arabidopsis TERMINAL FLOWER1 is involved in the regulation of flowering time and inflorescence development through transcriptional repression. Plant Cell. 2011;23:3172–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Hirsch CN, Foerster JM, Johnson JM, Sekhon RS, Muttoni G, Vaillancourt B, et al. Insights into the maize pan-genome and pan-transcriptome. Plant Cell. 2014;26:121–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Briggs WH, McMullen MD, Gaut BS, Doebley J. Linkage mapping of domestication loci in a large maize teosinte backcross resource. Genetics. 2007;177:1915–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao Q, Weber AL, McMullen MD, Guill K, Doebley J. MADS-box genes of maize: frequent targets of selection during domestication. Genet Res (Camb). 2011;93:65–75.View ArticleGoogle Scholar
- Sasaki A, Ashikari M, Ueguchi-Tanaka M, Itoh H, Nishimura A, Swapan D, et al. Green revolution: A mutant gibberellin-synthesis gene in rice - New insight into the rice variant that helped to avert famine over thirty years ago. Nature. 2002;416:701–2.View ArticlePubMedGoogle Scholar
- Bolduc N, Hake S. The maize transcription factor KNOTTED1 directly regulates the gibberellin catabolism gene ga2ox1. Plant Cell. 2009;21:1647–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller TA, Muslin EH, Dorweiler JE. A maize CONSTANS-like gene, conz1, exhibits distinct diurnal expression patterns in varied photoperiods. Planta. 2008;227:1377–88.View ArticlePubMedGoogle Scholar
- Wang RL, Stec A, Hey J, Lukens L, Doebley J. The limits of selection during maize domestication. Nature. 1999;398:236–9.View ArticlePubMedGoogle Scholar
- Wang H, Nussbaum-Wagler T, Li B, Zhao Q, Vigouroux Y, Faller M, et al. The origin of the naked grains of maize. Nature. 2005;436:714–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lemmon ZH, Bukowski R, Sun Q, Doebley JF. The role of cis regulatory evolution in maize domestication. PLoS Genet. 2014;10:e1004745.View ArticlePubMedPubMed CentralGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation; 2013.Google Scholar
- Ganal MW, Durstewitz G, Polley A, Berard A, Buckler ES, Charcosset A, et al. A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome. PLoS One. 2011;6:e28334.View ArticlePubMedPubMed CentralGoogle Scholar
- Browning BL, Browning SR. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009;84:210–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Wimmer V, Albrecht T, Auinger HJ, Schön C-C. Synbreed: a framework for the analysis of genomic prediction data using R. Bioinformatics. 2012;28:2086–7.View ArticlePubMedGoogle Scholar
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.View ArticlePubMedPubMed CentralGoogle Scholar
- Zeileis A, Grothendieck G. zoo: S3 infrastructure for regular and irregular time series. J Stat Softw. 2005;14:1–27.View ArticleGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:W64–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Statist. 2001;29:1165–88.View ArticleGoogle Scholar
- Hutter S, Vilella AJ, Rozas J. Genome-wide DNA polymorphism analyses using VariScan. BMC Bioinformatics. 2006;7:409.View ArticlePubMedPubMed CentralGoogle Scholar
- Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27:2156–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilcoxon F. Individual comparisons by ranking methods. Biometrics Bull. 1945;1:80–3.View ArticleGoogle Scholar
- Gresset S, Westermeier P, Rademacher S, Ouzunova M, Presterl T, Westhoff P, et al. Stable carbon isotope discrimination is under genetic control in the C4 species maize with several genomic regions influencing trait expression. Plant Physiol. 2014;164:131–43.View ArticlePubMedGoogle Scholar
- Utz HF. PLABSTAT - A computer program for statistical analysis of Plant Bred experiments. Version 3A. Stuttgart: Universität Hohenheim; 2011.Google Scholar