Original Research ARTICLE
Association Analysis of Grain-setting Rates in Apical and Basal Spikelets in Bread Wheat (Triticum aestivum L.)
- 1Key Laboratory of Wheat Biology and Genetic Improvement for Low and Middle Yangtze Valley (Ministry of Agriculture), Lixiahe Agricultural Institute of Jiangsu Province, Yangzhou, China
- 2Key Laboratory of Crop Gene Resources and Germplasm Enhancement, Ministry of Agriculture/The National Key Facility for Crop Gene Resources and Genetic Improvement/Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing, China
- 3Institute for Chemical Ecology, Shanxi Agricultural University, Taigu, China
- 4Agricultural Science, School of Veterinary and Life Sciences, Murdoch University, Murdoch, WA, Australia
- 5Hubei Collaborative Innovation Centre for Grain Industry/College of Agriculture, Yangtze University, Jingzhou, China
The rates of grain-setting in apical and basal spikelets in wheat directly affect the kernel number per spike (KNPS). In this study, 220 wheat lines from 18 Chinese provinces and five foreign countries were used as a natural population. Phenotypic analysis showed differences in grain-setting rates between apical and basal spikelets. The broad-sense heritability of grain-setting rate in apical spikelets (18.7–21.0%) was higher than that for basal spikelets (9.4–16.4%). Significant correlations were found between KNPS and grain numbers in apical (R2 = 0.40–0.45, P < 0.01) and basal (R2 = 0.41–0.56, P < 0.01) spikelets. Seventy two of 106 SSR markers were associated with grain setting, 32 for apical spikelets, and 34 for basal spikelets. The SSR loci were located on 17 chromosomes, except 3A, 3D, 4A, and 7D, and explained 3.7–22.9% of the phenotypic variance. Four markers, Xcfa2153-1A202, Xgwm186-5A118, Xgwm156-3B319, and Xgwm537-7B210, showed the largest effects on grain numbers in apical and basal spikelets. High grain numbers in apical and basal spikelets were associated with elite alleles. Ningmai 9, Ning 0569, and Yangmai 18 with high grain-setting rates carried larger numbers of favorable alleles. Comparison of grain numbers in basal and apical spikelets of 35 Yangmai and Ningmai lines indicated that the Ningmai lines had better grain-setting rates (mean 21.4) than the Yangmai lines (16.5).
Wheat is one of the most important cereal crops in the world. As the third-largest crop in China, the average grain yield of wheat has reached 4986 kg/ha, and the maximum was more than 12,000 kg/ha (National Bureau of Statistics of China, http://data.stats.gov.cn). Improved grain yield is still the most important objective in wheat breeding. Spike number per square meter, kernel number per spike (KNPS) and thousand grain weight (TGW) are important grain yield components. At a fixed spike number per square meter and fixed TGW, improvements in grain yield will depend largely on increasing KNPS (Slafer and Andrade, 1989; Fischer, 1993, 2008, 2011; Magrin et al., 1993). The numbers of grain produced in apical and basal spikelets become key factors for achieving high kernel numbers per spike (Arisnabarreta and Miralles, 2006; Acreche et al., 2008).
The physiological processes of grain formation in apical and basal splikelets include floret differentiation, grain development, and floret degeneration. The number of fertile florets is the key factor determining KNPS (Stockman et al., 1983; Sibony and Pinthus, 1988; González et al., 2003; Shitsukawa et al., 2009). Flower development in wheat spikes is not synchronized. The central spikelets develop first, followed by the distal spikelets. Florets that developed first can continue to develop and form seeds whereas the later florets are more subject to non-viability and loss of fertility (Kirby, 1988; Bancal, 2009; Ferrante et al., 2013). Moreover, floret fertility varies with position on the spike (González-Navarro et al., 2015). Previous studies elucidated causes of low fertility in apical and basal spikelets from the perspective of physiology. Floret differentiation and development are also influenced by environmental factors, and grain-setting rates are affected by various factors such as cultivation practices (Fischer, 1993; Demotes-Mainard et al., 1999; Ferrante et al., 2010), temperature (Dreccer et al., 2014), and light (González et al., 2003).
Low floret fertility in the apical and basal regions of the ear is a very common biological phenomenon in cereal crops, as shown by the “bald” phenomenon in maize (Satoh-Nagasawa et al., 2006; Meng et al., 2007; Gallavotti et al., 2011). More detailed studies have been conducted in rice (Yamagishi et al., 2004; Li et al., 2009; Cheng et al., 2011; Tan et al., 2011; Akter et al., 2014). Gene SP1 was shown to be a PTR transporter following map-based cloning and mutation caused severe degradation of basal spikelets (Li et al., 2009). A recessive gene for apical floret degradation in rice was fine mapped (Akter et al., 2014). There are almost no reports on isolation of genes involved in apical spikelet degradation. In wheat, there are few reports on grain-setting rates of apical and basal spikelets. However, there are many studies on spikelet infertility (Slafer and Andrade, 1989; Li et al., 2007; Ma et al., 2007; Wu et al., 2012; Xu et al., 2014). Those reports showed that additive effects play a major role in floret fertility through the study of heritability of yield-related traits and the mode of gene function using two crosses of common wheat (Salman and Heyne, 1987). QTLs associated with floret sterility were detected in 13 chromosomal regions using a recombinant inbreb line (RIL) population derived from the cross Nanda 2419 × Wangshuibai and an immortalized F2(IF2) population (Ma et al., 2007). QTLs QSspn.nau-4A and QSspn.nau-5A were detected in both populations. One QTL for floret infertility was detected in the interval Xwmc112-Xbarc168 in chromosome 2D using a Xiaoyan 54 × Jing 411 RIL population grown under different levels of nitrogen and phosphorus (Xu et al., 2014). Reduced height (Rht) genes also indirectly affect grain-setting rates (Miralles et al., 1998; Álvaro et al., 2008). One study showed that semi-dwarf varieties with Rht genes had significantly increased kernel numbers per spike, especially kernel numbers in the apical and basal spikelets (Miralles et al., 1998). Therefore, genetic information on control of grain-setting rates is of significant theoretical and applied value in understanding the underlying mechanisms.
In this study, we used 220 wheat varieties from 18 Chinese provinces and five foreign countries as a natural population to analyse the heritability of grain-setting rates in apical and basal spikelets in different environments. SSR loci associated with grain-setting rates in apical and basal spikelets were studied using association analysis. To some extent the study revealed possible reasons for significant differences in grain-setting rate between two popular series of wheat varieties, namely Yangmai and Ningmai cultivars, which have been widely cultivated in the lower Yangtze River region. The study provides evidence that increased KNPS can be achieved by improvement of grain-setting rates in the apical and basal spikelets. Our association results indicate that marker-assisted selection can be used to increase levels of grain setting in wheat.
Materials and Methods
A natural population formed by 220 wheat varieties was used for association analysis. There were 204 varieties from 18 provinces in China and seven foreign varieties. The origins of nine varieties were unknown. The seven foreign varieties from Italy (3), the U.S.A. (1), Mexico (1), Chile (1), and Japan (1) are historic donors that made important contributions to wheat breeding in China (Zhuang, 2003). The 18 provincial origins for the Chinese varieties were Jiangsu (65), Henan (25), Shandong (19), Shaanxi (18), Sichuan (16), Anhui (13) Hunan (10), Hubei (6), Beijing (7), Hebei (7), Gansu (4), Zhejiang (3), Fujian (3), Shanxi (3), Heilongjiang (2), Jiangxi (1), Guizhou (1), and Yunnan (1). All grow normally in the Yangtze River region. Variety details are provided in Table S1.
Phenotyping and Genotyping
The 220 wheat varieties were planted in Jingzhou (JZ), Hubei province, and Yangzhou (YZ), Jiangsu province in 2013 and 2014. The notations 13JZ, 13YZ, 14JZ, and 14YZ represent the four experiments conducted over 2 years. Field trial entries were randomized and three replicates were used for each variety per experiment. Each plot was 1.33 m long and 0.75 m wide. Seeds were evenly sown in three rows at approximately 40 seeds per row. To achieve a plant density of 30 plants per row some seedlings were removed. At the maturity stage, plants in the middle row of each plot were selected for the measurement of KNPS. Kernel number for each spikelet from the apex downwards was recorded as GNAS1, GNAS2, and GNAS3, and kernel number for each spikelet from the base upwards was recorded as GNBS1, GNBS2, and GNBS3 (Figure S1).
Genomic DNA was extracted using the CTAB method (Sharp et al., 1989) and 106 SSR markers distributed on all 21 chromosomes were used. Among the 106 SSR markers, 21 markers were previously reported to be associated with the KNPS (Zhang et al., 2012) and 17 SSR markers were chosen based on linkage with KNPS and spikelet infertility (Guo et al., 2015).
DNA fragment sizes were detected using an ABI 3730 DNA sequencer (Applied Biosystems, Foster City, USA) and software Genemapper V 3.7 (Applied Biosystems; http://www.appliedbiosystems.com.cn/) was used for data recording. SPSS 21.0 (Genetic Diversity Evaluation http://www.brothersoft.com/ibm-spss-statistics-469577.html) was used in statistical analysis of phenotypic variation. PowerMarker V 3.25 software was used for evaluation of genetic diversity (Liu and Muse, 2005) and Structure V 2.3.2 was used for evaluation of genetic structure of the natural population using the 106 SSR markers (Pritchard et al., 2000). The number of subsets was determined with the ΔK model (Evanno et al., 2005). SPAGeDi software was used for calculation of genetic covariance between pairs of individuals (Loiselle et al., 1995; Yu et al., 2005). Based on the model of Q + K (Zhang et al., 2010b), association analysis of grain-setting rate and SSR markers was performed using TASSEL 2.1 (Bradbury et al., 2007) (http://www.maizegenetics.net/). The genetic effects of favorable alleles were calculated using ANOVA program. Based on the data of phenotype and kinship between individuals, broad heritabilities in different environments were calculated according to the formula of h2 /() using TASSEL 2.1, where was genetic variance and represented the residual variance.
Related phenotypes for grain-setting rates in apical and basal spikelets (GNAS1, GNAS2, GNAS3, GNBS1, GNBS2, GNBS3, and KNPS) were investigated in four environments (13JZ, 13YZ, 14JZ, and 14YZ). The average coefficient of variation for each trait within the population was between 10.89 and 87.85%, indicating wide variation in grain-setting rate in the natural population, especially in regard to GNBS1. The average grain numbers per spike for GNAS1, GNAS2, and GNAS3 in the four environments were 1.39, 1.54, and 1.78, whereas those for GNBS1, GNBS2, GNBS3, were 0.46, 1.44, and 2.60, respectively (Table 1). The results indicated differences in grain-setting rates between apical and basal spikelets. Similar numbers of grains were produced in all apical spikelets. However, a significant factor was that the grain-setting rate was obviously low for GNBS1.
Broad-sense heritabilities for apical spikelets, namely GNAS1, GNAS2, GNAS3 were 18.68, 20.49, and 21.00%, respectively, whereas for basal spikelets the corresponding values were 9.44, 12.40, and 16.37% for GNBS1, GNBS2, and GNBS3, respectively (Table 1). Broad-sense heritabilities for GNAS1 and GNBS1 were lower than for other spikelets, indicating that the terminal spikelets were subject to greater environmental influence. In general, the heritabilities for apical spikelets were higher than for basal spikelets and increased from the apex downwards and from the base upwards.
Correlation Between Grain-setting Rates in Apical and Basal Spikelets, and Kernel Number per Spike
Analyses of regression between KNPS and total kernel number of the three spikelets from apical or basal spikes were performed in all four environments. Strong positive correlation (R2) of 0.40–0.45 (P < 0.01) and 0.41–0.56 (P < 0.01) in apical and basal spikelets, respectively, occurred in all four environments.
Correlation analysis of kernel number was also performed between the apical and basal spikelets, as well as with KNPS. There were significant correlations (P < 0.01) in all comparisons (Figure 1, Table S2). On average, correlations between KNPS and grain number of apical and basal spikelets were in the range of 0.31–0.48. The average correlation of grain number within apical spikelets was in the range of 0.87–0.94 compared to 0.69–0.89 within basal spikelets.
Figure 1. Analysis of regression between grain-setting in the whole spike, and apical and basal spikelets, respectively. Blue dots represent the phenotypic value and the red dotted lines represent the predicted regression line. The correlation between grain-setting in apical spikelets and kernel number per spike in 13JZ (A), 13YZ (B), 14JZ (C), and 14YZ (D); Correlations between the grain-setting in basal spikelets and kernel number per spike in 13JZ (E), 13YZ (F), 14JZ (G), and 14YZ (H).
Allelic Diversity and Genetic Structure Analysis
Eight hundred and seventy nine allelic variations were detected in the population using 106 SSR markers and the range in major allele frequency (MAF) was from 0.164 to 0.986, with an average of 0.541. The range in allele number per locus was 2–24, with an average of 8.3. The range of polymorphism information content (PIC) was from 0.027 to 0.908, with an average of 0.555 (Table S3) indicating that the natural population has a high level of genetic diversity.
To eliminate spurious associations, genetic structure (Q-value) and the individual kinship coefficient (K-value) were evaluated. Two subgroups were evident (Figure 2A). The maximum value of ΔK at the point of K = 2 further demonstrated two subgroups (Figure 2B). For the whole population, 76.5% of individual kinship coefficients were in the range of 0–0.05 (Figure S3). Low coefficients between individuals indicate no or little relationship between them.
Figure 2. Population structure of 220 wheat cultivars based on 106 genome-wide SSR markers. (A), genetic structure produced by Structure V2.3.2; (B), number of sub-populations estimated by ΔK at a range of K-values.
Association Analysis Between Grain Setting-related Phenotypes and SSR Markers
In the association analysis for grain-setting rates of apical and basal spikelets and KNPS using 106 SSR markers (including 879 alleles), 72 significant association signals were detected at 36 SSR loci, including 32 signals for apical grain-setting rates at 18 SSR loci, 34 signals for the basal grain-setting rate at 20 SSR loci, and six signals for KNPS at five SSR markers. The range in phenotypic explanation rates was 3.71–22.93%. The 36 SSR loci were distributed on chromosomes except 3A, 3D, 4A, and 7D (Table 2, Figure S2). Significant association signals detected in more than two environments included Xcfa2153-1A (GNAS2 and GNAS3), Xgwm148-2B (GNAS1), Xgwm102-2D (GNAS2), Xgwm108-3B (GNBS3), Xbarc56-5A (GNBS3), Xgwm126-5A (GNBS1), Xgwm186-5A (GNBS1), Xgwm415-5A (GNBS3), Xgwm55-6D (GNBS3), Xwmc168-7A (GNBS1), and Xgwm537-7B (KNPS, GNBS2, and GNBS3). Thirteen SSR markers were significantly associated with more than two phenotypes including Xcfa2153-1A (GNAS1, GNAS2, and GNAS3), Xgwm403-1B (GNAS3 and GNBS3), Xwmc147-1D (GNAS1, GNBS1, and GNBS2), Xgwm312-2A (GNAS2 and GNAS3), Xgwm148-2B (GNAS1, GNAS2, and GNAS3), Xgwm429-2B (GNAS1 and GNAS3), Xgwm102-2D (GNAS1, GNAS2, and GNAS3), Xgwm156-3B (GNAS2 and GNBS2), Xgwm149-4B (GNAS3 and GNBS1), Xcfd52-5D (GNAS1, GNAS2, and GNAS3), Xgwm190-5D (GNAS3 and GNBS3), Xgwm292-5D (KNPS and GNAS3), and Xgwm537-7B (KNPS, GNBS1, GNBS2, and GNBS3). Xgwm102-2D, Xgwm108-3B, Xgwm186-5A, Xgwm415-5A, Xcfd52-5D, and Xgwm55-6D were associated with grain-setting rate QTLs in previous studies (Peng et al., 2003; Zhang et al., 2010a, 2012; Guo et al., 2015; Table 2).
Table 2. Seventy-two significant association signals (P < 0.01) involving 36 SSR loci and seven phenotypic traits.
Favorable Alleles and their Genetic Effects
Favorable alleles and their genetic effects for each locus were further analyzed following association analysis (Table 2, Table S4). The largest genetic effect for KNPS was associated with Xgwm537-7B210 (2.18, 14YZ). Xcfa2153-1A202 had the largest genetic effect on apical spikelets, the values of GNAS1, GNAS2, GNAS3 were 0.31 (13YZ), 0.20/0.29 (13YZ and 14YZ), and 0.20/0.34 (13YZ and 14YZ), respectively, while basal spikelets genetic effects for GNBS1, GNBS2, and GNBS3 were 0.23 (13YZ) at Xgwm186-5A118, 0.87 (13JZ) at Xgwm156-3B319, and 0.31/0.33 (13JZ and 13YZ) at Xgwm537-7B210, respectively. Favorable alleles with frequencies greater than 50% were Xgwm403-1B136 (GNAS3 and GNBS3), Xgwm102-2D144 (GNAS1, GNAS2, and GNAS3), Xbarc56-5A131 (GNBS3), Xgwm499-5B124 (GNBS3), Xcfd52-5D282 (GNAS1, GNAS2, and GNAS3), Xgwm190-5D210 (GNBS3), Xgwm212-5D96 (GNBS1), and Xgwm292-5D208 (KNPS and GNAS3). The high frequencies suggest that varieties possessing these alleles had been subjected to strong selection pressure in modern breeding.
A yield QTL was reported previously at Xgwm186 on chromosome 5AL (Peng et al., 2003; Figure 3A). The same locus was detected in this study for grain-setting rate in GNBS1 in environments 13JZ and 13YZ (Figure 3B). The superior allele had 118 bp and its frequency in the population was 9.35% (Figure 3C). It conferred the highest genetic effect among those affecting the first basal spikelet (Table S4) and increased grain number by 0.14 and 0.23 in environments 13JZ and 13YZ, respectively (Figure 3D).
Figure 3. Favorable allele of Xgwm186 and genetic effect on grain-setting in the basal spikelets. (A) QTL locus at Xgwm186-5AL (Peng et al., 2003); (B) Association signals in the natural population using a mixed linear model (P < 0.01); green dots represents the associated signal of Xgwm186 in different environments; (C) Frequency distribution of Xgwm186-5A alleles in the population; green bar represents favorable allele of 118 bp; (D) Genetic effects of Xgwm186-5A118 in environments 13JZ and 13YZ.
The relationship between favorable allele number and grain-setting rates in both apical and basal spikelets in the entire population is portrayed in Figure 4. In general, higher numbers of favorable alleles led to higher grain-setting rates in apical and basal spikelets; for example, Ningmai 9, Ning 0569, and Yangmai 18 (an offspring of Ningmai 9) carried relatively high numbers of favorable alleles at 34, 33, and 32, respectively (Table S5).
Figure 4. Relationships between grain numbers in apical and basal spikelets and the number of favorable alleles. (A) The relatedness between the grain numbers in apical spikelets and number of favorable alleles. (B) The relatedness between the grain numbers in basal spikelets and number of favorable alleles.
For rate of grain-setting in apical spikelets, there were 15, 16, and 16 favorable alleles carried by Yangmai 18, Ning 0569, and Ning 9, respectively. Apical spikelet grain-setting rates in those varieties were better than others in the Yangmai and Ningmai series (Figure 5A). Similarly, high grain-setting rates also occurred in basal spikelets of Yangmai 18, Ningmai 9, and Ning 0569 with 17, 17, and 18 favorable alleles, respectively (Figure 5B).
Figure 5. Distribution of favorable alleles associated with grain numbers in apical and basal spikelets in Yangmai and Ningmai varieties. (A) Distribution of favorable alleles associated with grain-setting in apical spikelets in Yangmai and Ningmai varieties. (B) Distribution of favorable alleles associated with grain-setting in basal spikelets in Yangmai and Ningmai varieties.
Dissection of Differences in Grain-setting Rates in Apical and Basal Spikelets Between Yangmai and Ningmai Varieties
In the Yangtze River region Yangmai varieties predominate over Ningmai varieties although both groups occupy large areas. Grain-setting rates in the Ningmai group were higher than the Yangmai series (Table S5). Apical and basal spikelet grain-setting rates in the four environments were compared between 13 Ningmai and 22 Yangmai varieties. Apart from GNBS in environment 14JZ, the grain-setting rates for the Ningmai series were significantly higher than those for Yangmai series (Figure 6, Figure S4). The numbers of favorable alleles associated with grain-setting in apical and basal spikelets, and the sums of both in the Ningmai series (10.57, 10.79, 21.36) were higher than for the Yangmai series (8.14, 8.33, 16.48), respectively (Figure 6).
Figure 6. Relationships between average grain-setting rates in apical and basal spikelets, and numbers of favorable alleles.
Grain-setting Rates in Apical and Basal Spikelets in Wheat
Grain yield of wheat will be enhanced if the KNPS is increased by 1 or 2 while other yield factors remain constant (Arisnabarreta and Miralles, 2006; Acreche et al., 2008). There are several ways to increase KNPS; for example, by reducing plant height (Rebetzke and Richards, 2000), increasing spikelet number per spike (Rawson, 1970), and decreasing tiller number (Guo and Schnurbusch, 2015). The approach of the present study was to increase kernel numbers by increasing the numbers of grains in apical and basal spikelets. Significant correlations were shown between grain numbers in the apical and basal spikelets and total numbers per spike (P < 0.01) (Figure 1, Table S2). The results showed different grain-setting rates between apical and basal spikelets. Similar numbers of grains were produced in all apical spikelets. However, a significant factor was that the grain-setting rate was especially low for GNBS1 (Table 1). Our results thus confirmed that grain-setting rate in basal spikelets has high potential for increasing the KNPS.
There are relatively few studies on heritability of grain-setting rates. Maich and Zumelzú (1998) reported that narrow-sense heritability of spikelet fertility was 23% using an F2:3 population of hexaploid triticale. According to Martino et al. (2015) spikelet fertility is a moderately heritable trait. In the present study broad-sense heritability of grain-setting in apical and basal spikelets was 18.68–21.00% and 9.44–16.37%, respectively. Owing to the low heritability it would be difficult to obtain rapid progress by phenotypic selection at the early stages of a breeding program.
Analysis of the Grain-setting Rates in Apical and Basal Spikelets Partitions Variation in the Kernel Number per Spike
Past QTL analyses of spikelet fertility focused on overall kernel numbers per spike (Slafer and Andrade, 1989; Li et al., 2007; Ma et al., 2007; Wu et al., 2012; Xu et al., 2014). In this study we addressed grain numbers in apical and basal spikelets and analyzed the genetic effects using 106 SSR markers. Thirty two markers were associated with grain-setting rate in apical spikelets and 34 were associated with grain-setting rate in basal spikelets (Table 2, Figure S2). Zhang et al. (2012) reported associations of KNPS with Xwmc24-1A, Xcfd52-5D, and Xgdm127-6D. In our study these loci were also related to grain numbers in apical spikelets (Table 2). Guo et al. (2015) found associations of Xgwm337-1D, Xgwm108-3B, Xgwm186-5A, Xgwm415-5A, Xgwm234-5B, and Xgwm55-6D with spikelet fertility, and again these loci were associated with grain numbers in basal spikelets in the present work. Thus, high spikelet fertility might be the consequence of higher grain-setting rate in basal spikelets. Instead, a gene near the Xgwm102-2D locus may contribute to grain number in the apical spikelets there by reducing spikelet infertility (Guo et al., 2015). Therefore, some SSR loci reported for KNPS may play a role through regulation of grain-setting rates in apical and basal spikelets. Thus, a significant part of the genetic basis of increased grain-setting per spike may be related to the fertility of apical and basal spikelets.
Prospects for Molecular Breeding for Kernel Number per Spike in Wheat
The germplasm in this study included several Yangmai and Ningmai varieties. Despite the greater use of Yangmai vaieties by farmers in the Yangtse River region, grain-setting rates in apical and basal spikelets were higher in Ningmai varieties than in Yangmai varieties (Figure 4). Thus, improvement of KNPS in Yangmai varieties should improve yields. Overall, the relatedness of the phenotype and the number of favorable alleles carried by 220 varieties showed that the grain-setting rates in apical and basal spikelets tended to increase with increasing numbers of favorable alleles (Figure 4). Heritability studies of grain-setting rates provided evidence of simple additive effects, consistent with previous results (Salman and Heyne, 1987; Mirabella et al., 2015). Of all varieties Ning 0569 and Ningmai 9 carried the highest numbers of favorable alleles (32 and 34, respectively) associated with grain-setting rates in apical and basal spikelets (Table S5). Therefore, Ning 0569 and Ningmai 9 could be suitable donors of genes for high grain-setting in wheat breeding.
So far, relatively few genetic studies have been conducted on grain-setting rates in wheat, and almost no effective method exists for solving issues of grain-setting rates in the apical and basal spikelets. In this study, 36 SSR loci were associated with grain numbers in apical and basal spikelets (Table 2). To some extent, the results also revealed the basis of higher KNPS in Ningmai varieties compared to the Yangmai group (Figure 6). Xgwm102-2D, Xgwm108-3B, Xgwm186-5A, Xgwm415-5A, Xcfd52-5D, and Xgwm55-6D were associated with multiple traits in multiple environments. These loci were reported previously to be associated with grain-setting QTLs on a total spike basis (Peng et al., 2003; Zhang et al., 2010a, 2012, 2014; Guo et al., 2015). Among favorable alleles associated with grain-setting rates in apical and basal spikelets the largest effects were conferred by Xcfa2153-1A202, Xgwm186-5A118, Xgwm156-3B319, and Xgwm537-7B210 and their frequencies within the population were 6.51, 9.35, 6.13, and 22.12%, respectively. Low frequencies indicate that particular favorable alleles have not been strongly selected in past wheat breeding. Therefore, to achieve high grain-setting rates more favorable alleles can be introduced to breeding populations and could be tracked by specific molecular markers.
The results from this study were based on a larger mass of phenotypic data obtained through association analysis between grain-setting rates in apical and basal spikelets and SSR markers. In genotypic analysis results may be constrained by limited marker types and marker number. The high density 90K iSelect SNP chip (Wang et al., 2014) could be used to identify further loci associated with grain-setting rates in apical and basal spikelets. Meanwhile, Ning 0569 and Ningmai 9 can be used as donors for accumulating QTLs in breeding lines (Lu et al., 2010; Korir et al., 2013).
Conceived and designed the experiments: JG, CH, and SC. Performed the experiments: JG, WS, YX. Analyzed the data: JG and CH. Contributed reagents/materials/analysis tools: YZ, BZ, XC, KC, and XZ. Wrote the manuscript: JG, JZ, and CH.
This work was supported by grants from the Chinese Ministry of Science and Technology (2011CB100100), the China Agricultural Research System (CARS-3-1-1), the National High-Tech R & D Program of China (2012AA101105), the Six Talent Peaks Project in Jiangsu Province (2013NY-049), and the Key R & D Program in Jiangsu Province (BE2015352-3).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We gratefully acknowledge help from Professor Robert A. McIntosh, University of Sydney, with English editing.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.01029
Figure S1. A schematic representation of the grain set on apical and basal spikelets.
Figure S2. Associations of seven phenotypic traits with 106 SSR markers illustrated as dot plots of compressed MLM at P < 0.01. The red dotted line indicates the threshold value of significant association.
Figure S3. Distribution of pairwise kinship coefficients among 220 bread wheat cultivars based on 106 whole genome SSR markers.
Figure S4. The grain number of the apical and basal spikelets among different environments in series of cultivars named Yangmai and Ningmai.
Table S1. Information of 220 wheat accessions used in this study.
Table S2. Pearson's correlation coefficients between phenotypic traits in the different environments.
Table S3. Allele number, MAF, and PIC of 106 polymorphic SSR markers detected in the association panel.
Table S4. Favorable alleles and genetic effects of 56 SSR loci significantly (P < 0.01) associated with seven phenotypic traits.
Table S5. The number of favorable alleles carried by 220 wheat accessions used in this study.
KNPS, Kernel number per spike; GNAS1, Grain number in the first apical spikelet; GNAS2, Grain number in the second apical spikelet; GNAS3, Grain number in the third apical spikelet; GNBS1, Grain number in the first basal spikelet; GNBS2, Grain number in the second basal spikelet; GNBS3, Grain number in the third basal spikelet.
Acreche, M. M., Briceño-Félix, G., Sánchez, J. A. M., and Slafer, G. A. (2008). Physiological bases of genetic gains in Mediterranean bread wheat yield in Spain. Eur. J. Agron. 28, 162–170. doi: 10.1016/j.eja.2007.07.001
Akter, M. B., Piao, R., Kim, B., Lee, Y., Koh, E., and Koh, H. J. (2014). Fine mapping and candidate gene analysis of a new mutant gene for panicle apical abortion in rice. Euphytica 197, 387–398. doi: 10.1007/s10681-014-1074-8
Álvaro, F., Isidro, J., Villegas, D., García del Moral, L. F., and Royo, C. (2008). Old and modern durum wheat varieties from Italy and Spain differ in main spike components. Field Crops Res. 106, 86–93. doi: 10.1016/j.fcr.2007.11.003
Arisnabarreta, S., and Miralles, D. J. (2006). Floret development and grain setting in near isogenic two- and six-rowed barley lines (Hordeum vulgare L.). Field Crops Res. 96, 466–476. doi: 10.1016/j.fcr.2005.09.004
Bradbury, P. J., Zhang, Z. W., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Cheng, Z. J., Mao, B. G., Gao, S. W., Zhang, L., Wang, J. L., Lei, C. L., et al. (2011). Fine mapping of qPAA8, a gene controlling panicle apical development in rice. J. Integr. Plant Biol. 53, 710–718. doi: 10.1111/j.1744-7909.2011.01055.x
Demotes-Mainard, S., Jeuffroy, M. H., and Robin, S. (1999). Spike dry matter and nitrogen accumulation before anthesis in wheat as affected by nitrogen fertilizer: relationship to kernels per spike. Field Crops Res. 64, 249–259. doi: 10.1016/S0378-4290(99)00046-5
Dreccer, M. F., Wockner, K. B., Palta, J. A., McIntyre, C. L., Borgognone, M. G., Bourgault, M., et al. (2014). More fertile florets and grains per spike can be achieved at higher temperature in wheat lines with high spike biomass and sugar content at booting. Funct. Plant. Biol. 41, 482. doi: 10.1071/FP13232
Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Gallavotti, A., Malcomber, S., Gaines, C., Stanfield, S., Whipple, C., Kellogg, E., et al. (2011). BARREN STALK FASTIGIATE1 is an AT-hook protein required for the formation of maize ears. Plant Cell 23, 1756–1771. doi: 10.1105/tpc.111.084590
González, F. G., Slafer, G. A., and Miralles, D. J. (2003). Floret development and spike growth as affected by photoperiod during stem elongation in wheat. Field Crops Res. 81, 29–38. doi: 10.1016/S0378-4290(02)00196-X
González-Navarro, O. E., Griffiths, S., Molero, G., Reynolds, M. P., and Slafer, G. A. (2015). Dynamics of floret development determining differences in spike fertility in an elite population of wheat. Field Crops Res. 172, 21–31. doi: 10.1016/j.fcr.2014.12.001
Guo, J., Hao, C. Y., Zhang, Y., Zhang, B. Q., Cheng, X. M., Qin, L., et al. (2015). Association and validation of yield-favored alleles in Chinese cultivars of common wheat (Triticum aestivum L.). PLoS ONE 10:e130029. doi: 10.1371/journal.pone.0130029
Korir, P. C., Zhang, J., Wu, K. J., Zhao, T. J., and Gai, J. Y. (2013). Association mapping combined with linkage analysis for aluminum tolerance among soybean cultivars released in the Yellow and Changjiang River Valleys in China. Theor. Appl. Genet. 126, 1659–1675. doi: 10.1007/s00122-013-2082-0
Li, S. B., Qian, Q., Fu, Z. M., Zeng, D. L., Meng, X. B., Kyozuka, J., et al. (2009). Short panicle1 encodes a putative PTR family transporter and determines rice panicle size. Plant J. 58, 592–605. doi: 10.1111/j.1365-313X.2009.03799.x
Li, S. S., Jia, J. Z., Wei, X. Y., Zhang, X. C., Li, L. Z., Chen, H. M., et al. (2007). A intervarietal genetic map and QTL analysis for yield traits in wheat. Mol. Breed. 20, 167–178. doi: 10.1007/s11032-007-9080-3
Lu, Y. L., Zhang, S. H., Shah, T., Xie, C. X., Hao, Z. F., Li, X. H., et al. (2010). Joint linkage-linkage disequilibrium mapping is a powerful approach to detecting quantitative trait loci underlying drought tolerance in maize. Proc. Natl. Acad. Sci. U.S.A. 107, 19585–19590. doi: 10.1073/pnas.1006105107
Ma, Z. Q., Zhao, D. M., Zhang, C. Q., Zhang, Z. Z., Xue, S. L., Lin, F., et al. (2007). Molecular genetic analysis of five spike-related traits in wheat using RIL and immortalized F2populations. Mol. Genet. Genomics 277, 31–42. doi: 10.1007/s00438-006-0166-0
Magrin, G. O., Hall, A. J., Baldy, C., and Grondona, M. O. (1993). Spatial and interannual variations in the photothermal quotient: implications for the potential kernel number of wheat crops in Argentina. Agr. Forest Meteorol. 67, 29–41. doi: 10.1016/0168-1923(93)90048-M
Maich, R., and Zumelzú, D. M. D. (1998). Spikelet fertility and meiotic index in hexaploid triticale. Realized heritabilities and correlated responses for agronomical traits. Cereal Res. Commun. 26, 433–437.
Martino, D. L., Abbate, P. E., Cendoya, M. G., Gutheim, F., Mirabella, N. E., and Pontaroli, A. C. (2015). Wheat spike fertility: inheritance and relationship with spike yield components in early generations. Plant Breed. 134, 264–270. doi: 10.1111/pbr.12262
Meng, Z. D., Zhang, F. J., Ding, Z. H., Sun, Q., Wang, L. M., Guo, Q. F., et al. (2007). Inheritance of ear tip-barrenness trait in maize. Agr. Sci. China 6, 628–633. doi: 10.1016/S1671-2927(07)60092-3
Mirabella, N. E., Abbate, P. E., Ramirez, I. A., and Pontaroli, A. C. (2015). Genetic variation for wheat spike fertility in cultivars and early breeding materials. J. Agr. Sci. doi: 10.1017/S0021859614001245. [Epub ahead of print].
Miralles, D. J., Katz, S. D., Colloca, A., and Slafer, G. A. (1998). Floret development in near isogenic wheat lines differing in plant height. Field Crops Res. 59, 21–30. doi: 10.1016/S0378-4290(98)00103-8
Peng, J. H., Ronin, Y., Fahima, T., Roder, M. S., Li, Y. C., Nevo, E., et al. (2003). Domestication quantitative trait loci in Triticum dicoccoides, the progenitor of wheat. Proc. Natl. Acad. Sci. U.S.A. 100, 2489–2494. doi: 10.1073/pnas.252763199
Rebetzke, G. J., and Richards, R. A. (2000). Gibberellic acid-sensitive dwarfing genes reduce plant height to increase kernel number and grain yield of wheat. Crop Pasture Sci. 51, 235–246. doi: 10.1071/AR99043
Satoh-Nagasawa, N., Nagasawa, N., Malcomber, S., Sakai, H., and Jackson, D. (2006). A trehalose metabolic enzyme controls inflorescence architecture in maize. Nature 441, 227–230. doi: 10.1038/nature04725
Sharp, P. J., Chao, S., Desai, S., and Gale, M. D. (1989). The isolation, characterization and application in the Triticeae of a set of wheat RFLP probes identifying each homoeologous chromosome arm. Theor. Appl. Genet. 78, 342–348. doi: 10.1007/BF00265294
Shitsukawa, N., Kinjo, H., Takumi, S., and Murai, K. (2009). Heterochronic development of the floret meristem determines grain number per spikelet in diploid, tetraploid and hexaploid wheats. Ann. Bot. 104, 243–251. doi: 10.1093/aob/mcp129
Stockman, Y. M., Fischer, R. A., and Brittain, E. G. (1983). Assimilate supply and floret development within the spike of wheat (Triticum aestivum L.). Funct. Plant Biol. 10, 585–594. doi: 10.1071/PP9830585
Tan, C. J., Sun, Y. J., Xu, H. S., and Yu, S. B. (2011). Identification of quantitative trait locus and epistatic interaction for degenerated spikelets on the top of panicle in rice. Plant Breed. 130, 177–184. doi: 10.1111/j.1439-0523.2010.01770.x
Wang, S. C., Wong, D., Forrest, K., Allen, A., Chao, S., Huang, B. E., et al. (2014). Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. Plant Biotechnol. J. 12, 787–796. doi: 10.1111/pbi.12183
Xu, Y. F., Wang, R. F., Tong, Y. P., Zhao, H. T., Xie, Q. G., Liu, D. C., et al. (2014). Mapping QTLs for yield and nitrogen-related traits in wheat: influence of nitrogen and phosphorus fertilization on QTL expression. Theor. Appl. Genet. 127, 59–72. doi: 10.1007/s00122-013-2201-y
Yamagishi, J., Miyamoto, N., Hirotsu, S., Laza, R. C., and Nemoto, K. (2004). QTLs for branching, floret formation, and pre-flowering floret abortion of rice panicle in a temperate japonica × tropical japonica cross. Theor. Appl. Genet. 109, 1555–1561. doi: 10.1007/s00122-004-1795-5
Yu, J. M., Pressoir, G., Briggs, W. H., Bi, I. V., Yamasaki, M., Doebley, J. F., et al. (2005). A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 38, 203–208. doi: 10.1038/ng1702
Zhang, D. L., Hao, C. Y., Wang, L. F., and Zhang, X. Y. (2012). Identifying loci influencing grain number by microsatellite screening in bread wheat (Triticum aestivum L.). Planta 236, 1507–1517. doi: 10.1007/s00425-012-1708-9
Zhang, J., Dell, B., Biddulph, B., Khan, N., Xu, Y. J., Luo, H., et al. (2014). Vernalization gene combination to maximize grain yield in bread wheat (Triticum aestivum L.) in diverse environments. Euphytica 198, 439–454. doi: 10.1007/s10681-014-1120-6
Zhang, L. Y., Liu, D. C., Guo, X. L., Yang, W. L., Sun, J. Z., Wang, D. W., et al. (2010a). Genomic distribution of quantitative trait loci for yield and yield-related traits in common wheat. J. Integr. Plant Biol. 52, 996–1007. doi: 10.1111/j.1744-7909.2010.00967.x
Zhang, Z. W., Ersoz, E., Lai, C. Q., Todhunter, R. J., Tiwari, H. K., Gore, M. A., et al. (2010b). Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 42, 355–360. doi: 10.1038/ng.546
Keywords: bread wheat, released cultivars, grain numbers in apical spikelets, grain numbers in basal spikelets, association analysis
Citation: Guo J, Zhang Y, Shi W, Zhang B, Zhang J, Xu Y, Cheng X, Cheng K, Zhang X, Hao C and Cheng S (2015) Association Analysis of Grain-setting Rates in Apical and Basal Spikelets in Bread Wheat (Triticum aestivum L.). Front. Plant Sci. 6:1029. doi: 10.3389/fpls.2015.01029
Received: 11 September 2015; Accepted: 06 November 2015;
Published: 20 November 2015.
Edited by:Tiegang Lu, Chinese Academy of Agricultural Sciences, China
Reviewed by:Meixue Zhou, University of Tasmania, Australia
Xiao Han, Chinese Academy of Agricultural Sciences, China
Copyright © 2015 Guo, Zhang, Shi, Zhang, Zhang, Xu, Cheng, Cheng, Zhang, Hao and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work.