Seed is consumed as food and animal feed, providing more than 70% of caloric intake around the world, additionally it is also a fundamental component of the plant life cycle, as they store the genetic information necessary for the next generation of plants to disperse, establish, develop and eventually reproduce to maintain the species  . Seed vigor is an important and complex agronomic trait, determined by several factors including genetic and physical purity, mechanical damage and physiological condition, characterized by maintaining a high seed vigor and stable content after storage    , and required to ensure the rapid and uniform emergence of plants in the field under different environmental conditions. High vigor seeds make a great advantage for growth and production potential, which can enhance germination rates, resistance to environmental stresses, and crop yields   . Therefore, farmers and growers are constantly looking for high quality seeds able to ensure uniform germination and growth in field and to increase production.
Seed vigor essentially depends on the ability to withstand prolonged storage and the deleterious effects of aging. Seed vigor during storage can be defined as the maximum time period that pure seeds retain germination viability when stored under ideal environmental conditions and therefore represents an important trait for the conservation of seed resources. It varies among the different species due to natural variability and is usually regarded to be related with seed longevity or seed storability traits   . A reliable assay is essential to accurately phenotype the response to seed storability. However, studies of seed longevity under conventional or optimal storage conditions would take years to complete and therefore so-called accelerated aging or controlled deterioration tests (CDT) have been developed to assess the vigor of seed lots and to predict their relative longevity by aging seeds rapidly at elevated temperature and relative humidity (RH) as an alternative to analyze this property more efficiently    .
Although the environment during seed formation, harvest, and especially storage is important for seed vigor, genetic factors also largely affect seed vigor     . Genetics provides a powerful approach such as linkage analysis and, more recently, association mapping for genetic dissection of physiological and molecular bases of phenotypic traits such as seed longevity  . The former relies on trait segregation in a population derived from a bi-parental cross, and has been used to identify QTL for seed vigor under conventional storage conditions or CDT in rice, barley, wheat, oilseed rape and model plant Arabidopsis thaliana         . The latter is a population-based method that the mapping population consists of a set of unrelated accessions by the detection of linkage disequilibrium between a trait and a genetic marker  . However, seed vigor was not reported so far in all species. In addition, proteome analysis of seed vigor in Arabidopsis thaliana, maize revealed common features the CDT or conventionally aged seeds    .
In the present study, 208 and 212 F10 RILs derived from the two crosses between Yu82 and Shen137, Yu537A and Shen137 were used to detect QTL for four traits of seed vigor under control and three aging treatment conditions. The first aim of this research was to identify the QTL traits of seed vigor. The second aim was to integrate QTLs detected across two RIL populations to identify true QTLs, and furthermore was to integrate candidate gene analyses with related traits of seed vigor QTL mapping across two populations to test the effects of numerous candidate genes for the traits known from other species on the natural variations for the traits in maize.
2. Materials and Methods
2.1. Plant Materials and Artificial Aging Treatments
First, confirm that you have the correct template for your paper size. This template has been tailored for output on the custom paper size (21 cm * 28.5 cm). The two connected populations used in the study consisted of 208 and 212 F10 RILs derived by single-seed descent from two crosses of Yu82 × Shen137 and Yu537A × Shen137, which were referred as Population 1 (Pop 1) and Population 2 (Pop 2) and used to identify QTLs for related trait of seed vigor, respectively.
The artificial aging treatment was used the same method described by Zeng et al.  . The seeds of two populations and three parents were reproduced in the winter in Hainan Province in 2015. After harvest, the seeds were fully dried under natural conditions. The seeds of each genotype were divided into four portions (60 seeds choosing to ensure sowing quality of every portion) for artificial aging treatments. All the seeds were placed in Nylon mesh belt firstly, then were treated at 45˚C ± 1˚C and 90% relative humidity for 0, 2, 4, and 6 days (0d, 2d, 4d and 6d) by using a thermostatic moisture regulator, respectively. Every treatment followed a randomized complete block design with three replications. Among treatments, 0d treatment acted as control.
2.2. Germination Experiment and Related Trait of Seed Vigor Evaluation
The template is used to format your paper and style the text. All margins, column widths, line spaces, and text fonts are prescribed; please do not alter them. You may note peculiarities. For example, the head margin in this template measures proportionately more than is customary. This measurement and others are deliberate, using specifications that anticipate your paper as one part of the entire journals, and not as an independent document. Please do not revise any of the current designations.
The germination experiment conducted at 25˚C in artificial climate chamber in 2016. The method of germination experiment was as follows: the first, selecting diameter of 0.05 - 0.2 mm of fine sand as sprout bed and the sand was treated by high-handed sterilization pan at 120˚C for two hours; the second, using a germination container of 16 × 8 holes that the diameter of each hole was 40 mm; the third, each hole was filled with 3.5 cm thick sand and put 2 seeds in it, and then used 1.5 cm thick sand to cover them; the last, the germination containers sowing seeds were left in artificial climate chamber for a temperature 25˚C, a relative humidity 65% and illumination conditions 4000 lx, the photoperiod was 14/10 (day/night). The number of germinated seeds was counted daily. The data of related traits for 8 days after sowing were used for QTL analysis when obvious differences between the parents were observed. After daily statistics finished, 5 plants of each RIL were selected randomly to measure the seedling length, respectively. The germination percentage (GP) was calculated as GP = n/N × 100%, where n is the total number of germination seeds, N is the total number of seeds. The germination index (GI) was calculated as: GI = ΣGt/Dt, where the Dt is the germination time, Gt is the number of germinated seeds on the time. The vigor index (VI) was calculated as: VI = GI × SL, where the SL is the seedling length on day 8. The simple vigor index (SVI) was calculated as: SVI = GP × SL. The mean germination time (MGT) was calculated as: MGT = ΣGt × Dt/GP, where the sense of Gt and Dt as above.
2.3. Statistical Analysis of Phenotypic Data
The trait values for each RIL were reported as the average from five plants in each replication. The overall performance was the average over the three replications from each artificial aging treatment. Analysis of variance (ANOVA) was carried out to estimate genetic variation for all the measured traits among the RILs using the general linear model procedure of the statistical software SPSS 17.0. Descriptive statistics and simple correlation coefficients (r) between the traits were calculated using the above statistical software.
2.4. Construction of Genetic Linkage Map
A total of 3072 pairs of single nucleotide polymorphism (SNP) markers were selected from the more than 800,000 SNPs to genotype the 420 RILs and three parents. We analyzed polymorphisms of 3072 SNP markers between two pairs of parents, Yu82/Shen137 and Yu537A/Shen137. Ultimately, 1397 and 1371 SNP markers had polyphisms between the two parents, respectively. Chi-square values were generated for 2768 SNP markers, 225 and 232 SNP markers showed serious segregation distortion and failed to be assigned to any linkage in the two populations. The linkage analysis was done with JoinMap version 4.0. Two genetic linkage maps were constructed with 1172 and 1139 SNP markers using Joinmap version 4.0  , and the total length 1629.61 cM with an average interval of 1.39 cM for Pop.1 and 1681.75 cM with an average interval of 1.48 cM for Pop.2  .
2.5. QTL Analysis
QTL analysis was conducted using composite interval mapping (CIM) with WinQTLcart 2.5 software  . For CIM, Model 6 of the Zmapqtl dodule was employed for detecting QTL and their effects, specifying the five markers identified by stepwise regression that explained most of the variation for a given trait as forward and backward parameters and a window size of 10 cM on either side of the markers flanking the test site  . To identify an accurate significance threshold for each trait, an empirical threshold was determined by performing 1000 random permutations  . QTL position was assigned to relevant region at the point of the maximum likelihood odds ratio (LOD). QTL confidence interval was calculated by subtracting one LOD unit on each side from the maximum LOD position  .
For the additive effects of QTL, positive and negative values indicated that alleles from the normal maize inbred lines Yu82/Yu537A and the maize inbred line Shen137 increased the trait scores, respectively. QTL were named according to ‘‘q’’ +‘‘artificial aging treatment days” + ‘‘trait abbreviation’’ + ‘‘population code’’ + ‘‘−’’ + ‘‘chromosome number’’ + ‘‘QTL number’’.
2.6. Meta-QTL Analysis
To integrate QTLs information for the measured traits located in the two connected RIL populations, the genetic linkage maps were integrated and consensus QTLs were identified by meta-analysis   . The QTLs mapped in the two connected RIL populations were projected on the integrated map using their positions and confidence intervals shared by two linkage maps. Some controversial markers between two linkage maps were deleted, which could effectively improve the accuracy of projection.
Meta-analysis was performed by using BioMercator2.1 software  . The Akaike Information Criterion (AIC) was used to select the QTL model on each chromosome  . According to this, the QTL model with the lowest AIC value is considered a significant model indicating the number of meta-QTL. The number of mQTLs that best fitted the results on a given linkage group was determined based on a modified Akaike criterion  . Meta-QTL were named according to ‘‘q’’ +‘‘artificial aging treatment days” + ‘‘trait abbreviation’’ + ‘‘population code’’ + ‘‘−’’ + ‘‘chromosome number’’ + ‘‘QTL number’’.
3. Results & Discussion
3.1. Phenotypic Performance of Traits Associated with Seed Vigor in Three Parents and Connected Two RILs
The values of GI, VI and SVI were obviously decreased and the values of MGT were markedly increased after three treatment conditions compared with control in parents and two populations. For three parents, the values of GI, VI and SVI were higher for Yu82 and Yu537A than Shen137 under four aging treatments, while the reverse was true for MGT. trait differences were also found among three parents under each treatment. For RILs, the values presented a large range of variability with transgressive segregation exceeding values of high values parent. All traits showed normal distribution in the two RIL populations and differed substantially under various treatment conditions (Table 1).
Significant positive correlations were consistently observed for GI, VI and SVI from two RIL populations under control and after various aging treatments except for between VI and SVI from Pop. 1 under 4d aging treatment, while MGT and GI, VI, SVI showed significant negative correlations except for between MGT and SVI from Pop. 2 under 0 and 2d aging treatments (Table 2).
3.2. QTL Detection for Each Trait in Two Connected Populations
A total of 74 QTLs for GI, VI, SVI and MGT were detected in two connected populations under control and after three aging treatment conditions, with 40 QTLs in Pop. 1 and 34 QTL in Pop. 2 (Table 3). These QTLs were located on all chromosomes. The contributions to phenotypic variations for a single QTL ranged from 5.33% to 13.74%, with 10 QTLs over 10% and 1 QTL over 13%.
Nine QTLs in Pop.1 and ten QTLs in Pop.2 were identified and located on all chromosomes except for chromosomes 2 and 10 under four aging treatment conditions. The contribution rates of these QTLs ranged from 5.56% to 13.74% of total phenotypic variance (Table 3). The positive alleles of q2GI1-3, q2GI1-5-1, q2GI1-5-2, q4GI2-4 and q6GI2-9 were derived from Shen137 to contribute towards an increase in values of GI. There were qGI1-6 from Pop.1 consistently mapped in the same marker interval SYN31854-PZE-106102131 after 2d and 4d aging treatments, qGI2-8-2 from Pop.2 in the interval PZB00865.2-PZE-108073195 after 0d and 2d aging treatments, and qGI2-8-1 from Pop.2 in the interval PZE-105077135-PZE-105082252 after 0d, 2d and 4d aging treatments. Among these QTLs, QTL qGI2-8-1 was responsible for 10.25, 11.73 and 8.22% of phenotypic variance, and qGI2-8-2 responsible for 10.82 and 7.07% of phenotypic variance, respectively.
Eighteen QTLs were mapped for VI under control and after three aging treatments in the two populations, nine in Pop.1 and nine in Pop.2. They were distributed across the whole genome, except for chromosomes 9 and 10 with contribution to phenotypic variation for a single QTL from 5.39 to 8.89% (Table 3). The positive alleles of qnVI1-1-1, q2VI1-1, q4VI1-1, q4VI21-6 and q6VI1-7 in Pop.1 and of qnVI2-8 in Pop.2 were contributed by Yu82/Yu537A. However, there was no QTL identified at same marker intervals under different aging treatment conditions in the two populations.
Eleven QTLs for SVI were detected, with four in Pop.1 and seven in Pop.2 under all aging treatment environments. They were distributed across chromosomes 2, 3, 4 5 and 6, and explained 5.72 to 12.11% of the phenotypic variation (Table 3). Among these QTLs, two in Pop.1 and one in Pop.2 were derived from Yu82/Yu537A to increase in the trait values. The positive alleles of qSVI2-2
Table 1. Phenotypic performance of seed vigor related traits of three parents and two ROLs under favorable and three artificial aging conditions.
Table 2. Phenotypic correlations among seed vigor related traits for the two RILs based on average under favorable and three artificial aging conditions.
The table above diagonal line mean the phenotypic correlations of Yu82 × Shen137; the table below diagonal line mean the phenotypic correlations of Yu537 × Shen137. a. artificial aging time; b. GI germination index, VI vigor index, SVI simple vigor index, MGT mean germination time.*Significant at P = 0.05, **Significant at P = 0.01.
Table 3. QTL detected for nine traits in the two RIL populations under favorable and three artificial aging conditions.
a. GI germination index, VI vigor index, SVI simple vigor index, MGT mean germination time; b. A is represent for additive effect.
derived from Shen137 in Pop.2 were consistently identified at the marker interval SYN8399-PZE-102122951 after 2d and 6d aging treatment, respectively. And QTL qSVI2-2 accounted for 12.11% and 10.67% of phenotypic variance, respectively.
Under control and after three aging treatments, nineteen and nine QTLs for MGT were detected on all chromosomes except for chromosome 10 in Pop.1 and Pop.2 respectively, which explaining from 5.33% to 11.33% of the phenotypic variation (Table 3). The all alleles except for eleven QTL in Pop1 and of six in Pop.2 were derived from Shen137 to increase in the trait values. qMGT2-2 and q MGT2-5 in Pop.2 were consistently detected at the same marker interval PZE-102131962-SYN15147 and PZE-102131962-PZE-102136708 under control and after 2d aging treatment, respectively.
Component heads identify the different components of your paper and are not topically subordinate to each other. Examples include Acknowledgements and References and, for these, the correct style to use is “Heading 5”. Use “figure caption” for your Figure captions, and “table head” for your table title. Run-in heads, such as “Abstract”, will require you to apply a style (in this case, non-italic) in addition to the style provided by the drop down menu to differentiate the head from the text.
Text heads organize the topics on a relational, hierarchical basis. For example, the paper title is the primary text head because all subsequent material relates and elaborates on this one topic. If there are two or more sub-topics, the next level head should be used and, conversely, if there are not at least two sub-topics, then no subheads should be introduced. Styles named “Heading 1”, “Heading 2”, “Heading 3”, and “Heading 4” are prescribed.
3.3. Genetic Map Integration and mQTL Analysis from the Experimental Results
The consensus genetic map included 1712 SNP makers and was 1712.6 cM long with an average of 1.00 cM between markers on the basis of two connected populations. Meta-analysis identified mQTL that was associated with the variation of multiple traits measured. So 20 mQTLs were detected from the seventy-four initial mapped QTLs in the two RIL populations for four traits measured by using meta-analysis (Table 4). Seventy-two initial QTLs (97.3%) were integrated in these regions. The 20 mQTLs were located on all chromosomes except for 10, four on chromosome 4, three on chromosomes 1, 3 and 5, two on chromosome 2, and one on chromosomes 6, 7, 8 and 9. On average, one mQTL included 3.68 initial QTLs with a range from two to nine for 1-3 traits. It is worth noting that all initial QTLs with R2 > 10% were integrated in 7 mQTLs: mQTL1-1, mQTL2-1, mQTL2-2, mQTL5-1, mQTL5-3, mQTL5-4 and mQTL8.
Initial QTLs included in mQTL8 were all detected for GI under control and after three aging treatments, that in mQTL2-2 for three traits, that mQTL3-2 for two traits and those in mQTL1-3, mQTL3-1, mQTL4-4, mQTL5-1, mQTL5-3, mQTL5-4 for one traits under two conditions, respectively, and those in the rest of the mQTLs for two traits under one condition.
Seed vigor depends on their physiological and genetic conservation potential and on conditions encountered during storage   . Seed vigor strongly influences on plant stand establishment, which challenge crop breeders to produce high quality seeds for stabilizing crop yield. A key to achieving the challenge associated with seed vigor is elucidation of the molecular mechanisms underlying the traits. However, there are few papers about the research of the molecular mechanisms associated with the traits in maize. In this study, two sets of connected RIL populations were evaluated for four traits of seed vigor under control and three aging treatment conditions. Seventy-four QTLs were located in total, with 40 QTLs in Pop.1 and 34 in Pop.2. Individual QTL explained from 5.33 to 13.74% of the traits associated with seed vigor, with 11 QTLs over 10%. Seventy-two initial QTLs were integrated in 20 meta-QTLs (mQTL) using meta-analysis and all initial QTLs with contributions (R2) > 10% were integrated in mQTLs. Twenty candidate genes for seed vigor were consistently mapped in 13 corresponding mQTLs regions and mainly involved in glycolytic pathway and protein metabolism. These QTLs could provide useful information for marker-assisted selection in improving performance of seed vigor. At the same time, the results had important reference for the fine mapping of the major QTLs and validation of the potential candidate genes, and analyzing the molecular mechanism of seed vigor in maize.
Table 4. mQTLs for nine traits in the two connected RIL populations under favorable and three artificial aging conditions.
4.1. Comparison of QTL Located and Synthesis of Initial QTL in the Two Connected RIL Populations
Comparing the results in both populations, the research found that two common QTLs were located in the same chromosome regions, one for GI at bin 1.09 and the other for MGT at bin 10.03; two common QTLs were located near chromosome regions, one for VI at bin 4.08 - 4.09 and the other for SVI at bin 4.08 - 4.09, respectively. These QTLs showed great consistency across both populations, which might deserve further study in molecular marker-assisted selection (MAS). Although the two populations exhibited certain similarity because they share one of the parental lines, population-specific QTLs were also found. Three QTLs on chromosomes 9 and 10 were detected in Pop1 under 1 - 2 treatments. These QTLs were contributed by the unique parental line Yu82. One factor of inconsistency of QTL was the use of different bi-parent populations. Difference in the genetic background between the two populations exists because they share only one common parental line.
Twenty mQTLs were detected from the total of 74 initial QTLs for the traits measured using mQTL analysis proposed by Goffinet and Gerber  in this study, containing all QTLs of R2 over 10%. Genomic regions for the traits were mainly focused on chromosomes 1, 3, 4 and 5, even though the mQTLs were mapped on all chromosomes except for chromosome 10. Among the mQTLs, mQTL2-2, mQTL5-4 and mQTL8 with a major effect (R2 of initial QTL > 10%) included in 6 - 9 initial QTLs for 1 - 3 traits under 2 - 4 aging treatment conditions in one population or both populations. So the genomic regions might be hot spots of the important QTLs for the traits. Fine mapping of the mQTLs and validation of the potential candidate genes were a reliable and feasible strategy for QTL cloning. Therefore, Near-isogenic lines for three mQTLs are now in construction using marker-assisted selection, and the mQTLs might be used to improve maize seed vigor in the near future.
4.2. Associations between QTL and Candidate Genes in Maize
Maize is grown widely throughout the world in a range of agro-ecological environments. Companies sell many different registered hybrids in the same ecological region because of the high profits from producing hybrid maize. The amount of hybrid seeds produced by a few large companies could be excess to the requirements of market, so that the extra seeds stored to sell the following year. In the process of seed storage, the aspects of adverse factors associated with high temperature, high moisture content, and high oxygen gas pressure would induce free radical-mediated lipid peroxidation, enzyme inactivation, protein degradation, disruption of cellular membranes, and damage to genetic (nucleic acids) integrity    . Therefore, they could accelerate seed deterioration, and decrease seed quality and vigor. To understand further the genetic basis of seed vigor variation during seed storage, the association between QTLs and genes known to be involved in seed vigor in Arabidopsis  and maize  were investigated through a bioinformatics approach in maize. Twenty-two candidate genes for seed vigor were mapped in 15 mQTL intervals (Table 4). The candidate genes were involved in response to stress, molecular chaperones, hydrolase, energy, cell growth/division, protein destination and storage, signal transduction, translation and protein metabolism, amino acid metabolism and other processes.
4.3. Protein Metabolism Is Major Components of Seed Vigor during Seed Accelerated Aging
Seed germinations had requirement for some protein metabolisms regrouped several functions during seed accelerated aging, containing protein translocation, folding, thermotolerance, oligomeric assembly, and switching between active and inactive protein conformations  . The simultaneous impairment of these functions was closely linked with the loss of seed vigor  . Ten candidate genes identified within corresponding mQTLs detected under aging treatments were related to the protein metabolism in this study. 226494943 (Glutathione S-transferase) for defence was associated with the mQTL1-3 region affecting GI and VI; 146325682 (thioredoxin peroxidase) corresponded to the mQTL7 region for VI and MET. These gene families are major antioxidant enzymes in plant cells, as they were significantly accumulated in germinating seeds   , suggesting the antioxidant enzymes play an important role in seed viability maintenance. Plenty of works also revealed that the activities of anti-oxidases (like SOD, CAT, APX) and the content of antioxidant (like L-ascorbic acid and glutathione) changed upon aging treatment in seeds of soybean, rice, sunﬂower, maize etc. AT1G57720 (elongation factor 1-g2) for protein metabolism was mapped to the mQTL4-3 interval affecting SVI and MGT. AT1G09640 (Elongation factor 1B-g) for translation and protein metabolism corresponded to the same mQTL6-1 interval associated with GI, VI and MGT. At5g19550 (Asp aminotransferase) for protein metabolism was mapped to the same mQTL3-2 interval for VI, SVI and MGT. 242056533 (Hypothetical ACD_ScHsp26_like), 195605946 (HSP16.9) and 162459222 (HSP17.2) for molecular chaperone were located in the same mQTL3-1 region affecting GI and MGT; At5g51440 (hsp20/alpha crystallin family protein) for molecular chaperone was mapped to the mQTL3-3 interval associated with VI and SVI. The candidate genes (Hsp26, hsp20, HSP16.9 and HSP17.2) with chaperone activities in aged maize seeds might disturb signal transduction such as in responses to stresses like heat shock   and also favored targets for oxidation, presumably because they act as shields protecting other proteins against ROS damage  . 298201206 (Stress-related protein) was located in the mQTL8 interval for GI, VI and MGT. It has specific role in the detoxication of a wide range of exogenous and endogenous toxicants in artificially aged maize seeds  .
In addition, three candidate genes identified within mQTL detected under aging treatments were associated with the storage protein and protease in this study. 22284 (Vicilin-like embryo storage protein) were mapped to the same mQTL3-3 region associated with VI, SVI and MGT; 226508796 (CAAX prenyl protease 1) for hydrolase were associated with the mQTL4-2 interval affecting VI and MGT; At5g67360 (Cucumisin-like Ser protease) for protease corresponded with the mQTL5-3 region for MGT. The results showed artificial aging would increase proteases and breake down stored proteins, impaired metabolism and energy supply, and ultimately resulted in seed deterioration   .
4.4. The Glycolytic Pathway Is Affected during Seed Accelerated Aging
In this study, six candidate genes identified within mQTL detected under aging treatments were associated with the glycolytic pathway. 226532762 (glyoxalase family protein) in responses to stresses were mapped in the mQTL1-1 region affecting VI and MGT; 45238345 (Aldehyde dehydrogenase) involved in response to stress and At1g70730 (Phosphoglucomutase) for energy were located in the mQTL6 interval connected with GI, VI and MGT; 326509331 (V-type (H+)-ATPase domain) for energy was mapped in the same mQTL4-4 interval associated with SVI and MGT; 305671643 (ATP synthase beta subunit) and 197132370 (ATP synthase CF1 beta subunit) associated with energy were located in the mQTL2-1 and mQTL5-1 regions, respectively. The key enzymes/proteins participated in glycolysis, TCA cycle, electron transport chain and oxidative phosphorylation glycolysis by seed aging, which played a major role in the maintenance of the intracellular redoxstate and the maintenance of seed vigor  . The results showed that seeds experienced an oxidative stress during aging treatments and mounted a protective response through modification of the glycolytic pathway    . We believe that the candidate genes participate in metabolism and energy supply play important roles in seed aging and seed vigor.
4.5. Other Pathways Were Involved in Seed Vigor during Seed Accelerated Aging
In this way, other pathways were also involved in seed vigor after accelerated aging except for the above two pathways. For example, the embryo cell undergoes active division and expansion during seed germination. The events might be affected by accelerated aging treatment, as seeds germinated at a much lower speed after aging  . In our study, 224061823 (Predicted cyclin-dependent kinase A, CDK) for cell growth/division was mapped in the same mQTL1-3 interval associated with GI and VI; 162459414 (MEK homolog1) for cell growth/division were mapped to the same mQTL3-2 interval associated with VI, MGT and SVI, In addition, 195658029 (Lipoprotein) for lipid metabolism was mapped to the same mQTL4-2 interval associated with VI and MGT. CDK and MEK homolog 1 played a pivotal role in the regulation of the eukaryotic cell cycle.
In conclusion, the consistency of the QTL and the candidate genes identified in this study provided valuable information for further finding quantitative trait genes. The alleles for seed vigor could be useful targets for marker-assisted selection to produce germplasm of aging resistance.
This work was supported by grants from the National Natural Science Foundation of China (No.U1504315), Corn Industry Technology System in Henan Province (No.S2010-02) and the Scientific Research Foundation of Henan University of Science and Technology (No.13480067).
 Hilhorst, H.W., Finch-Savage, W.E., Buitink, J., et al. (2010) Dormancy in Plant Seeds: Dormancy and Resistance in Harsh Environments. Springer, Berlin, 43-67.
 Gupta, M.L., George, D.L. and Parwata, I.G. (2005) Effect of Harvest Time and Drying on Super Sweet Corn Seed Quality. Seed Science & Technology, 33, 167-176.
 Bentsink, L., Alonso-Blanco, C., Vreugdenhil, D., et al. (2000) Genetic Analysis of Seed-Soluble Oligosaccharides in Relation to Seed Storability of Arabidopsis. Plant Physiology, 124, 1595-1604.
 Clerkx, E.J., Blankestijn-De Vries, H., Ruys, G.J., et al. (2004a) Genetic Differences in Seed Longevity of Various Arabidopsis mutants. Physiology Plant, 121, 448-461.
 Miura, K., Lyn, S.Y., Yano, M. and Nagamine, T. (2002) Mapping Quantitative Trait Loci Controlling Seed Longevity in Rice (Oryza sativa L.). Theoretical and Applied Genetic, 104, 981-986.
 Xue, Y., Zhang, S.Q., Yao, Q.H., et al. (2008) Identification of Quantitative Trait Loci for Seed Storability in Rice (Oryza sativa L.). Euphytica, 164, 739-744.
 Nagel, M., Rosenhauer, M., Willner, E., et al. (2011) Seed Longevity in Oilseed Rape (Brassica napus L.) Genetic Variation and QTL Mapping. Plant Genetic Research, 9, 260-263.
 Rehman Arif, M.A., Nagel, M., Neumann, K., et al. (2012) Genetic Studies of Seed Longevity in Hexaploid Wheat using Segregation and Association Mapping Approaches. Euphytica, 186, 1-13.
 Rajjou, L., Lovigny, Y., Groot, S.P., et al. (2008) Proteome-Wide Characterization of Seed Aging in Arabidopsis: A Comparison between Artificial and Natural Aging Protocols. Plant Physiology, 148, 620-641.
 Xin, X., Lin, X.H., Zhou, Y.C., et al. (2011) Proteome Analysis of Maize Seeds: The Effect of Artificial Aging. Physiologia Plantarum, 143, 126-138.
 Wu, X.L., Liu, H.Y., Wang, W., et al. (2011) Proteomic Analysis of Seed Viability in Maize. Acta Physiologiae Plantarum, 33, 181-191.
 Lander, E.S., Green, P., Abrahamson, J., et al. (1987) MAPMAKER: An Interactive Computer Package for Constructing Primary Genetic Linkage Maps of Experimental and Natural Populations. Genomics, 1, 174-181.
 Han, Z., Ku, L., Zhang, Z., et al. (2014) QTLs for Seed Vigor-Related Traits Identified in Maize Seeds Germinated under Artificial Aging Conditions. PLoS ONE, 9, e92535.
 Trachsel, S., Messmer, R., Stamp, P., et al. (2009) Mapping of QTLs for Lateral and Axile Root Growth of Tropical Maize. Theoretical and Applied Genetic, 119, 1413-1424.
 Arcade, A., Labourdette, A., Falque, M., et al. (2004) BioMercator: Integrating Genetic Maps and QTL towards Discovery of Candidate Genes. Bioinformatics, 20, 2324-2326.
 Chardon, F.B., Virlon, L., Moreau, M., et al. (2004) Genetic Architecture of Flowering Time in Maize as Inferred from Quantitative Trait Loci Meta-Analysis and Synteny Conservation with Rice Genome. Genetics, 168, 2169-2185.
 Decourcelle, M., Perez-Fons, L., Baulande, S., et al. (2015) Combined Transcript, Proteome, and Metabolite Analysis of Transgenic Maize Seeds Engineered for Enhanced Carotenoid Synthesis Reveals Pleotropic Effects in Core Metabolism. Journal of Experimental Botany, 66, 3141-3150.
 Hyndman, D., Bauman, D.R. and Heredia, V.V. (2003) The aldo/keto Reductase Superfamily Homepage. Chemico-Biological Interactions, 143-144, 621-631.
 Kim, S.T., Wang, Y., Kang, S.Y., et al. (2009) Developing Rice Embryo Proteomics Reveals Essential Role for Embryonic Proteins in Regulation of Seed Germination. Journal of Proteome Research, 8, 3598-3605.
 Stone, J.M. and Walker, J.C. (1995) Plant Protein Kinase Families and Signal Transduction. Plant Physiology, 106, 451-457.
 Cabiscol, E., Piulats, E., Echave, P., et al. (2000) Oxidative Stress Promotes Specific Protein Damage in Saccharomyces cerevisiae. Journal of Biology Chemistry, 275, 27393-27398.
 Ching, T.M. and Schoolcraft, I. (1968) Physiological and Chemical Differences in Aged Seeds. Crop Science, 8, 407-409.
 Kocsy, G., Von Ballmoos, P., Ruegsegger, A., et al. (2001) Increasing the Glutathione Content in a Chilling-Sensitive Maize Genotype using Safeners Increased Protection against Chilling-Induced Injury. Plant Physiology, 127, 1147-1156.
 Sattler, S.E., Gilliland, L.U., Magallanes-Lundback, M., et al. (2004) Vitamin E Is Essential for Seed Longevity and for Preventing Lipid Peroxidation during Germination. Plant Cell, 16, 1419-1432.
 Job, C., Rajjou, L., Lovigny, Y., et al. (2005) Patterns of Protein Oxidation in Arabidopsis Seeds and during Germination. Plant Physiology, 138, 790-802.
 Kibinza, S., Vinel, D., Come, D., et al. (2006) Sunflower Seed Deterioration as Related to Moisture Content during Aging, Energy Metabolism and Active Oxygen Species Scavenging. Physiology Plant, 128, 496-506.
 Bingham, I.J. and Merritt, G.J. (1999) Effect of Seed Aging on Early Post-Germination Root Extension in Maize: A Spatial and Histological Analysis of the Growth-Zone. Seed Science & Technology, 27, 151-162.