Original Research ARTICLE
QTL Detection and Elite Alleles Mining for Stigma Traits in Oryza sativa by Association Mapping
- 1State Key Laboratory of Crop Genetics and Germplasm Enhancement, College of Agriculture, Nanjing Agricultural University, Nanjing, China
- 2Rice Research Institute, Chongqing Academy of Agricultural Sciences, Chongqing, China
Stigma traits are very important for hybrid seed production in Oryza sativa, which is a self-pollinated crop; however, the genetic mechanism controlling the traits is poorly understood. In this study, we investigated the phenotypic data of 227 accessions across 2 years and assessed their genotypic variation with 249 simple sequence repeat (SSR) markers. By combining phenotypic and genotypic data, a genome-wide association (GWA) map was generated. Large phenotypic variations in stigma length (STL), stigma brush-shaped part length (SBPL) and stigma non-brush-shaped part length (SNBPL) were found. Significant positive correlations were identified among stigma traits. In total, 2072 alleles were detected among 227 accessions, with an average of 8.3 alleles per SSR locus. GWA mapping detected 6 quantitative trait loci (QTLs) for the STL, 2 QTLs for the SBPL and 7 QTLs for the SNBPL. Eleven, 5, and 12 elite alleles were found for the STL, SBPL, and SNBPL, respectively. Optimal cross designs were predicted for improving the target traits. The detected genetic variation in stigma traits and QTLs provides helpful information for cloning candidate STL genes and breeding rice cultivars with longer STLs in the future.
Cultivated Asian rice (Oryza sativa L.) is one of the oldest domesticated crop species in the world and feeds more than one half of the world's population. This widespread utility is attributed to the large geographical range of O. sativa, which extends from 43°S (Australia) to 54°N (Mo River, North China) and from 7°E (Italy) to 117°W (California, USA). Two groups of genetically divergent subspecies, Oryza sativa subspecies indica and Oryza sativa subspecies japonica (Caicedo et al., 2007; Kumagai et al., 2010), are planted in large areas. Despite the continued debate regarding the origin(s) of domesticated Oryza sativa (Kovach et al., 2007; Sang and Ge, 2007; Molina et al., 2011; Huang et al., 2012; Civáñ et al., 2015; Huang and Han, 2015), both indica and japonica undoubtedly exhibit unique ecological distributions (Sang and Ge, 2007).
In China, rice is grown on 30 million hectares annually, and only India devotes a larger area to rice production. The annual rice production in China is 200.78 million tons, making China the largest rice producer in the world (Xie and Hardy, 2014). Each year, hybrid rice is planted on up to 50% of the total area in China devoted to rice production. The yield of hybrid rice cultivars is approximately 20% (1 ton per hectare) higher than that of conventional rice cultivars (Lu and Hong, 1999; Cheng et al., 2004). Despite the continuous improvements in cultivation techniques for F1 seed production over the last 10 years, the yield of hybrid rice seed production stagnated at 2.5 tons per hectare (Xie, 2009). A low stigma exsertion percentage is the main factor limiting further increases in the yield of F1 seed production in rice because it is a typical self-pollinating crop (Figure 1). To increase the percentage of rice stigma exsertion in rice, several studies have used QTL mapping to examine the genetic variation of stigma traits (Virmani and Athwal, 1973, 1974; Kato and Namai, 1987a; Virmani, 1994; Yamamoto et al., 2003; Uga et al., 2003a,b). To date, 50 QTLs that control stigma exsertion percentage have been mapped to all 12 rice chromosomes (Yamamoto et al., 2003; Uga et al., 2003a,b; Yu et al., 2006; Miyata et al., 2007; Hu et al., 2009; Yan et al., 2009; Li et al., 2014). Forty of the 50 QTLs were detected in the bi-parent derived segregating populations, and the remaining 10 were identified in the natural population. However, many environmental conditions, including high humidity, wind velocity, physical interruption, and low temperatures during the flowering period, influence stigma exsertion percentage (Beachell et al., 1938; Kato and Namai, 1987b), which may lower the accuracy of QTL mapping and complicate further studies. Therefore, to improve the outcrossing rate of the maternal parent, the trait stigma length (STL) can be more reliably measured than stigma exsertion in the mining of favorable alleles because STL is less influenced by external conditions than stigma exsertion percentage.
Figure 1. Scene of F1 hybrid seed production in two different growth stages and exerted stigmas after the palea and lemma enclosed in a single plant in rice. (A) Scene of F1 hybrid seed production in the heading stage. (B) Scene of F1 hybrid seed production in the filling stage. (C) Exerted stigmas (arrowheads) after the palea and lemma enclosed.
Several studies have shown that the stigma exsertion percentage is highly positively correlated with stigma length (Virmani and Athwal, 1973, 1974; Li and Chen, 1985; Kato and Namai, 1987a,b; Miyata et al., 2007). To the best of our knowledge, 23 QTLs that control stigma length have been detected, and these are distributed across all chromosomes with the exception of chromosomes 8 and 11 (Uga et al., 2003a, 2010; Yan et al., 2009; Liu et al., 2015). Most of these QTLs were identified using bi-parent-derived segregating populations. However, despite the availability of these scientific resources, linkage mapping is limited by the fact that only two alleles can be studied at any given locus in bi-parental crosses of inbred lines.
Recently, association mapping based on linkage disequilibrium (LD) has emerged as a popular method for mapping the loci responsible for natural variations, for mining natural genomic diversity, and for locating valuable genes (Zhu et al., 2008; Brachi et al., 2011; Weigel, 2012). In rice, association mapping has been used to exploit excellent alleles for many traits, including agronomic traits (Garris et al., 2005; Agrama et al., 2007; Huang et al., 2010, 2011; Zhao et al., 2011; Li et al., 2012; Vanniarajan et al., 2012; Dang et al., 2015; Yang et al., 2015), seed vigor traits (Cui et al., 2013; Dang et al., 2014; Rebolledo et al., 2015), and outcrossing traits (Yan et al., 2009). To the best of our knowledge, only one report (Yan et al., 2009) describes an association mapping-based study of stigma length.
The objectives of this study were to (1) evaluate the phenotypic diversity of stigma traits and the genetic architecture of the collection, (2) exploit elite alleles underlying stigma traits, and (3) predict the optimal cross combinations for improving stigma traits.
Materials and Methods
A total of 227 O. sativa accessions were obtained from China (154), Vietnam (57), and Japan (16). The details of these accessions, including their origins and accession IDs, are summarized in Table S1. The seeds of all accessions were collected, stored and supplied by the State Key Laboratory of Crop Genetics, and Germplasm Enhancement at Nanjing Agricultural University.
Field Planting and Trait Measurement
Plants from the 227 accessions were grown in the paddy field of Jiangpu Experimental Station at the Nanjing Agricultural University in Nanjing (32°07′N, 118°64′E), Jiangsu Province, China, from May to October in 2013, and 2014. The field experiments conducted in each of the two consecutive years was treated as two independent environments. The field trials followed a completely randomized block design with two replicates per year. Each plot contained five rows, with eight plants in each row. The rows were spaced 20 cm apart, and each plant was spaced at a distance of 17 cm according to standard agronomic management practices.
After pistil maturation, 10 spikelets were randomly collected from five plants of each accession prior to glume opening (approximately 10:00 a.m.). For each spikelet, the STL, SBPL, and SNBPL (Figure 2A) of the fertile flower were measured with a micrometer using a stereomicroscope (10 ×, MC50, Guangdong, China), and the average values served as the measurements for the accessions.
Figure 2. Morphology of pistil and the stigma in rice. (A) Names of rice pistil parts defined in this study. (B) Stigma morphology of Lianjing 4hao, with the minimum value of STL among 227 accessions. (C) Stigma morphology of Nongxiang 21, with the maximum value of STL among 227 accessions. (D) Stigma morphology of Yuedao 23. (E) Stigma morphology of Yuedao 34. The STL of Yuedao 23 and Yuedao 34 are equal. For Yuedao 23, SBPL is shorter than SNBPL. For Yuedao 34, SBPL is longer than SNBPL. Scale bar, 1 mm.
To understand the relationship between stigma traits and grain length (GL), fully filled grains were measured. In each year, five normally developed plants from each accession were harvested from the middle of each of the plots at the grain ripening stage and dried under natural conditions for GL investigation. Ten randomly selected grains (after the awns were removed) from each accession were lined up length-wise along an electronic digital Vernier caliper (http://www.guanglu.com.cn) to measure the grain length, and the average values served as the measurements for the accessions.
SSR Marker Genotyping
Genomic DNA was extracted from young and healthy leaf blades of each accession approximately 3 months after germination using the methods described by Gross et al. (2009). According to the rice molecular map and microsatellite database published by Temnykh et al. (2000) and McCouch et al. (2002), 249 SSRs scattered on 12 chromosomes were selected. The SSR primers were synthesized by Shanghai Generay Biotech Co., Ltd. (Shanghai, China).
Each 10-μl PCR reaction contained 10 mM Tris-HCl (pH 9.0), 50 mM KCl, 0.1% Triton X-100, 1.5 mM MgCl2, 0.5 nM dNTPs, 0.14 pM forward primer, 0.14 pM reverse primer, 0.5 U of Taq polymerase, and 20 ng of genomic DNA. DNA amplification was performed using a PTC-100™ Peltier Thermal Cycler (MJ Research™ Incorporated, USA) under the following conditions: (1) denaturation at 94°C for 5 min; (2) 34 cycles of denaturation at 94°C for 0.5 min, annealing at 55–61°C for 1 min, and extension at 72°C for 1 min; and (3) a final extension at 72°C for 10 min. The PCR products were separated on an 8% polyacrylamide gel for 1 h at 150 V and visualized using silver staining. One pair of SSR markers was used to detect one locus, and each polymorphic band at the same marker locus in the population was recorded as one allele. After screening the PAGE products, the size of each band was determined using Quantity One software (Bio-Rad Company, USA).
Population Genetic Structure
The genetic clusters in the 227 accessions were identified using STRUCTURE version 2.2 (Falush et al., 2007). This analysis was performed five times for each number of clusters (K) (from 2 to 10) with random starting points. We set the length of the burn-in period equal to 50,000 iterations and defined a run of 100,000 Markov Chain Monte Carlo (MCMC) replicates after burn in. A mean log-likelihood value over five runs at each K was used. If the mean log-likelihood value was positively correlated with the model parameter K, a suitable value for K could not be determined. In this situation, the optimal K value was determined through an ad hoc statistic, ΔK, based on the rate of change in [LnP(D)] between successive K values (Evanno et al., 2005). Non-admixed individuals in each genetic group were determined using a Q-matrix assignment greater than 0.9. A principal component analysis (PCA) was performed using the pcaMethods package (Stacklies et al., 2007) of R.2.11.1 (R Development Core Team, 2011) to examine the population structure. The first three principal components (PCs) were evaluated. The genetic distance was calculated based on 249 molecular markers using Nei's distance (Nei et al., 1983), and phylogenetic reconstruction was performed based on the neighbor-joining method implemented in PowerMarker version 3.25 (Liu and Muse, 2005). The tree used to visualize the phylogenetic distribution of the accessions and ancestry groups was constructed using MEGA version 4.0 (Tamura et al., 2007). The results from the three different approaches were then summarized and compared.
All basic statistical analyses were performed using the SAS package (SAS Institute Inc., Cary, NC, USA). The broad-sense heritability () was calculated based on the natural population through an analysis of variance using the formula /(σ+σ/n), where σ is the genetic variance, σ is the error variance, and n is the number of replicates.
The polymorphic information content (PIC) was used to measure the probability that two randomly selected alleles from a population were distinct. The number of alleles per locus, the gene diversity, and the PIC values were determined using PowerMarker version 3.25 (Liu and Muse, 2005).
The average standardized individual allele size of the SSRs was determined by calculating the mean standardized size of 249 SSR loci following the method reported by Vigouroux et al. (2003). This standardization ensured that each locus contributed equally to the average individual size of the genome.
The D′ value (Farnir et al., 2000) was calculated with TASSEL 3.0 software using 1000 permutations and was used to measure the level of linkage disequilibrium (LD) between loci (Bradbury et al., 2007). Before the association analysis, rare alleles with allele frequencies of 5% or less were removed from the dataset. According to the level of LD and the genetic distance among markers with intrachromosomal combinations, the regression equation of the LD as a function of changes in the genetic distance was calculated through a regression analysis. A LD decay plot was drawn to observe the relationship between LD and genetic distance.
A Mixed Linear Model (MLM) analysis, which can significantly reduce spurious marker-trait associations (type I errors showing false positives) resulting from the population structure, was performed using TASSEL 3.0 to calculate the associations between traits and markers (Bradbury et al., 2007). The matrices Q and K were used as covariants in the MLM analysis. The Q matrix was adapted from the analysis results obtained from Structure 2.2. The K matrix (kinship matrix) was obtained from the results of the relatedness analysis using SPAGeDi (Hardy and Vekemans, 2002). According to the correction method published by Benjamini and Hochberg (1995), a false discovery rate (FDR) of 0.05 was used as a threshold for significant associations. Based on the identified association locus, the “null allele” (non-amplified allele) was used to determine the phenotypic effects of other alleles (Breseghello and Sorrells, 2006).
To detect a phylogeographic signal in O. sativa, we divided the distribution range into six regions: Vietnam (8°30′–23°22′N, 102°8′–109°24′E), Southern China (20°45′–30°08′N, 97°31′–124°34′E), Eastern China (27°12′–35°20′N, 116°18′–123°00′E), Northern China (31°23′–42°37′N, 110°21′–119°45′E), Northeastern China (38°43′–53°33′N, 118°53′–135°05′E) and Japan (30°38′–44°09′N, 130°49′–145°15′E). This delimitation was based on the latitude of each region. We performed a locus-by-locus analysis of molecular variance (AMOVA) (Weir and Cockerham, 1984) based on genetic groups delimited by the Bayesian clustering method in the program Arlequin 3.5 (Excoffier and Lischer, 2010) to statistically verify the geographical structure using SSR and standard multi-locus frequency data.
Phenotypic Variations in the Rice Germplasm and Correlations among Traits
The mean value, standard deviation, skewness, and kurtosis for stigma traits measured in 227 rice varieties were calculated (Table 1). A significantly positive correlation was found between stigma traits in 2013 and 2014 (Figure S1). The materials exhibiting the minimum and maximum values are shown in Figures 2B,C. The stigmas exhibited longer SNBPLs and shorter SBPLs in some accessions (Figure 2D), and these trends were reversed in other accessions (Figure 2E). The average STL, SBPL, and SNBPL values of 227 accessions over the 2 years were 1.8 mm, 1.1 mm, and 0.7 mm, respectively. The phenotypic data for the STL, the SBPL, and the SNBPL in the studied population followed normal distributions (Table 1). Table 1 also shows that the average SBPL was greater than the average SNBPL. The broad-sense heritability values for the STL, SBPL, and SNBPL traits averaged over 2 years were 93.8, 92.8, and 88.6%, respectively. Furthermore, the three stigma traits exhibited significant positive correlations. The correlation coefficients between the STL and SBPL, the STL and SNBPL, and the SBPL and SNBPL were 0.803, 0.787, and 0.264, respectively (Table 2).
Table 1. Phenotypic characteristics for stigma length and grain length traits in 227 rice accessions across 2 years.
The highest average GL value over 2 years was 12.8 mm, the lowest average GL value over 2 years was 6.5 mm, and a skewed distribution was observed based on the skewness and kurtosis statistics (Table 1). The broad-sense heritability value averaged over 2 years was 98.3%. The correlation coefficients between the GL and STL, the GL and SBPL, and the GL and SNBPL were 0.710, 0.477, and 0.654, respectively (Table 2).
Genetic Variation in the Rice Germplasm
A marker analysis of the 227 accessions using 249 SSR molecular markers resulted in the detection of 2072 alleles. The numbers of alleles ranged from two (at locus RM7163_Chr11) to 18 (RM162_Chr6), with an average of 8.32 alleles per locus (Table S2). The genetic diversity averaged 0.699 and ranged from 0.083 (RM7163_Chr11) to 0.894 (RM162_Chr6) (Table S2). The PIC had a mean value of 0.668, ranged from 0.050 (RM7163_Chr11) to 0.885 (RM162_Chr6), and was mainly distributed between 0.509 and 0.809 (Table S2). A total of 210 markers (84%) were highly informative (PIC > 0.5), 27 (11%) were moderately informative (0.5 > PIC > 0.25), and 12 (5%) were slightly informative (PIC < 0.25).
The highest delta K values from the STRUCTURE analysis were obtained for eight clusters (Figure S2). The population structure data based on the Q matrix for each accession are summarized in Table S1, and the 227 accessions could be divided into eight subpopulations, viz. Group 1 to Group 8 and an admixed group (Figure S3A). Group 1 contained 38 accessions, which were modern improved varieties obtained mainly from Eastern China. Group 2 contained 45 accessions, which were mainly from Vietnam. Group 3 contained 20 accessions, and these were mainly obtained from Japan. Group 4 included 32 accessions, which consisted of landraces from Eastern China. Group 5 contained 21 accessions, which were mainly from Northern China. Group 6 contained 25 accessions, and these were mainly obtained from Southern China. Group 7 included 30 accessions, which were mainly from Northeastern China. Group 8 only contained four accessions obtained from North and Northeastern China. The admixed group contained the remaining 12 accessions.
Although the non-admixed individual cluster contained some admixed individuals, the neighbor-joining tree indicated that each of the eight groups corresponded to a distinct branch on the tree (Figures S3B, S4). These admixed individuals tended to be genetically similar to the group corresponding to their branch but exhibited a small amount of introgression from other groups.
The PCA results were essentially consistent with the STRUCTURE findings. The eight genetic groups were distinct in the principal component space of the three major PCs, whereas some individuals appeared mixed (Figure S3C). The three major PCs identified from the PCA explained 31.8% (11.4, 10.3, and 10.1%) of the total variance.
The following analysis did not include the accessions in the admixed group or Group 8 because Group 8 contained only four accessions.
Geography-Based Directional Evolution in the Allele Length of Subpopulations
To test whether the microsatellite size exhibited directional evolution in rice, we compared the average allele sizes in the geographically derived groups (Groups 1–7). Figure 3 shows the average individual microsatellite allele size for these seven groups. No significant differences were detected between Groups 1 and 4 (P > 0.05, two-tailed t-test), Groups 2 and 6 (P > 0.05, two-tailed t-test), or Groups 3 and 5 (P > 0.05, two-tailed t-test), whereas the rest groups exhibited significant differences, with P-values less than 0.001. The average individual SSR allele size in the subgroups from low latitudes (Groups 2 and 6) was 0.14, which is greater than the value obtained for subgroups from high latitudes (-0.04; Figure 3). These data indicate a decrease in allele size from south to north in the geographically derived groups.
Figure 3. Average individual microsatellite allele size for the seven subgroups (Groups 1–7). The standard error is presented.
A total of 18,120 of 30,876 pairs, including both interchromosomal, and intrachromosomal combinations, showed significant LD (based on D′, P < 0.05). Of these 18,120 pairs, 1530 were intrachromosomal pairs of SSR loci, which corresponds to a rate of 8.4%. The extent of LD was assessed for each of the seven groups (Table 3). For the average D′, Group 1 showed the lowest LD (0.399), whereas Group 5 exhibited the highest LD (0.468).
We used the D′ value corresponding to the intrachromosomal SSR loci and the genetic distance in each subpopulation to draw the attenuation map. Figure 4 shows that the D′ values decayed as the genetic distance (cM) increased. A regression analysis between the D′ value and the genetic distance of syntenic marker pairs revealed that the seven subpopulation genomes fitted the equation y = b lnx + c. The minimum distances of LD decay for Groups 1, 2, 3, 4, 5, 6, and 7 were 19.5 cM, 23.2 cM, 24.1 cM, 12.9 cM, 22.2 cM, 22.8 cM, and 21.7 cM, respectively. Thus, among the seven subpopulations, Group 4 exhibited the highest decay velocity with the shortest decay distance, whereas Group 3 demonstrated the lowest decay velocity.
Figure 4. Relationship between D′ value and genetic distance of syntenic marker pairs in subpopulations.
Genetic Differentiation across Subpopulations
The average FST among the seven subpopulations was 0.672, with the FST for each locus ranging from 0.197 for RM333_Chr 10 to 0.984 for RM6361_Chr 2. A results of a pairwise comparison based on the values of FST can be interpreted as the standardized population distances between two subpopulations. The pairwise FST values obtained in the present study ranged from 0.162 (between Groups 2 and 6) to 0.832 (between Groups 2 and 7), with an average value of 0.679 (Table 4). The AMOVA results indicated that 67.88% of the total genetic variation occurred between the subpopulations, whereas 32.12% occurred within the subpopulations (Table S3). These results indicate the existence of a high degree of genetic differentiation across the seven subpopulations.
Table 4. Pairwise estimates of FST based on 249 SSR loci among the seven model-based subpopulations.
Significant Marker-Trait Association Loci Detected across the Entire Population
A marker-trait association analysis based on an MLM revealed that six markers located on chromosomes 1, 2, 4, and 6 were associated with the STL (Table 5). The phenotypic variation explained (PVE) ranged from 2.9 to 12.5%. RM5753_Chr 6, which resides at 124.4 cM, had the highest PVE values for the STL, namely 12.5% in 2013, and 8.9% in 2014 (Table 5). Two markers distributed on chromosomes 4 and 6 were associated with the SBPL (Table 5), and of these, RM136_Chr 6 had the highest average PVE of 6.5% over the 2 years. Seven markers distributed on chromosomes 1, 2, 4, 6, and 11 were associated with the SNBPL (Table 5), and their PVE values ranged from 3.1 to 12.9%.
Table 5. Marker-trait associations with P-value less than 0.05, their equivalent false discovery rate probability (FDR), proportion of phenotypic variance explained (PVE), marker position on chromosome derived from 249 markers and 211 rice accessions.
Thirteen markers distributed on chromosomes 1, 2, 3, 4, 6, 7, 8, 10, and 11 were associated with the GL (Table 5), and their PVE values ranged from 2.2 to 22.9%. RM6327_Chr 11, which resides at 1.7 cM, had the maximum PVE for the GL, specifically 21.8% in 2013, and 22.9% in 2014 (Table 5).
Using the entire set of accessions, we identified 28 markers associated with stigma traits and grain length, including six markers associated with the STL, two related to the SBPL, seven associated with the SNBPL, and 13 associated with the GL. Nineteen of the 28 associations were in regions where the QTL associated with the given trait had been identified by Redoña and Mackill (1998), Tan et al. (2000), Uga et al. (2003a), Aluko et al. (2004), Li et al. (2004), Agrama et al. (2007), Bai et al. (2010), Uga et al. (2010), Wang et al. (2011, 2012), Dang et al. (2015) (http://www.gramene.org/); these QTLs are listed in Table S4. Nine loci identified in this study, including one for the STL, one for the SBPL, three for the SNBPL, and four for the GL, are novel (Table 5).
The marker RM280 was co-associated with the STL, SBPL and SNBPL. Five markers, i.e., RM5389, RM450, RM7598, RM280, and RM5753, were co-associated with the STL and SNBPL. Four markers, namely RM5389, RM450, RM7598, and RM280, were co-associated with the GL and STL and with the GL and SNBPL. Two markers, RM280 and RM136, were co-associated with the GL and SBPL (Table 5).
Discovery of Elite Alleles
The alleles with positive effects identified in this study were considered elite alleles for all four traits measured. A summary of the elite alleles and their typical carrier materials is shown in Table S5. The total numbers of elite alleles detected across the entire population for the STL, SBPL, and SNBPL were 11, 5, and 12, respectively. The allele RM450-135 bp showed the greatest phenotypic effect (0.235 mm) on the STL, and its typical carrier accession was Yuedao 32. The allele RM136-200 bp showed the greatest phenotypic effect (0.111 mm) on the SBPL, and its typical carrier accession was Yuexiangzhan. The allele RM450-135 bp showed its greatest phenotypic effect (0.233 mm) on the SNBPL, and its typical carrier accession was Yuedao 32.
A total of 30 elite alleles were detected for the GL. The allele RM136-200 bp exerted the greatest phenotypic effect (1.612 mm) on the GL, and its typical carrier accession was Yuexiangzhan (Table S5).
Correlations between the measured traits were observed, and the STL was significantly positively correlated with the SBPL and SNBPL. Furthermore, the GL was also significantly positively correlated with the STL, SBPL, and SNBPL. We identified one SSR marker co-associated with the STL and the SBPL, and the allele RM280-175 bp in this SSR marker locus simultaneously increased the phenotypic effect values of the STL and SBPL (Table 5 and Table S5). We identified five SSR markers co-associated with the STL and SNBPL, and the alleles RM5389-120, RM450-135, RM450-155, RM7598-115, RM280-175, RM5753-200, and RM5753-205 bp at these marker loci simultaneously increased the phenotypic effect values of the STL and SNBP (Table 5 and Table S5). We also identified four SSR markers co-associated with the STL and GL, and the alleles RM5389-120, RM450-135, RM450-155, RM7598-115, and RM280-175 bp at these marker loci simultaneously increased the phenotypic effect values of the STL and GL (Table 5 and Table S5). We further identified one SSR marker co-associated with the STL, SBPL, SNBPL, and GL, and the allele RM280-175 bp at this marker locus simultaneously increased the phenotypic effect values of the STL, SBPL, SNBPL, and GL (Table 5 and Table S5). These co-associated alleles have the correct sign with respect to trait correlations, and these data illustrate the genetic basis of trait correlations.
Optimal Cross Designs for Improving Target Traits
Based on the number of elite alleles that could be substituted into an individual plant and the expected phenotypic effects of the elite alleles that could be pyramided, the top five cross combinations for improving the STL, SBPL, SNBPL, and GL are proposed (Table S6). The elite alleles carried by the parents in excellent crosses and the corresponding phenotypic effects are listed in Table S7. Figure S5 shows the excellent parents in the superior cross for stigma traits. Certain accessions were found repeatedly in these proposed parental combinations (e.g., Yuedao 32 emerged four times in the combinations for the STL), indicating that these accessions possess unique elite alleles.
All three analyses, STRUCTURE, PCA and the neighbor-joining method, uncovered the same pattern of eight very distinct genetic groups in the samples of O. sativa investigated (Figure S3). Most accessions of these groups clearly belonged in unique subgroups, and minor accessions belonged in admixed subgroup (Table S1). In addition, the population structure was observed to be tied to the geographical origin, e.g., the accessions from Vietnam essentially clustered into Group 2, and the accessions from Northeastern China primarily clustered into Group 7. The distinctive geographical origins corresponding to the differences in ecological environments may have been partially responsible for the genetic differentiation. These findings demonstrate that natural populations of O. sativa, which is a self-pollinating species, show clear population subdivisions and a high amount of genetic diversity based on SSR markers resulting from historical mixtures among the subgroups.
No significant differences were found between Group 1 and Group 4, between Group 2 and Group 6, or among Group 3, Group 5, and Group 7 in the average standardized individual SSR allele sizes (Figure 3). The average standardized individual SSR allele sizes in Group 1 and Group 4 (mainly from Eastern China) is smaller than that in Group 2 and Group 6 (mainly from Southern China), but larger than that in Group 3, Group 5, and Group 7 (mainly from Northern China). This result indicates a tendency for directional evolution in SSRs during improvement of Oryza sativa from low (warmer climate) to high latitudes (cooler climate).
The average number of alleles per locus was 8.32 among the 227 accessions, which were genotyped using 249 markers (Table S2). This allele number per locus was higher than the values reported by Cho et al. (2000), Jain et al. (2004), Garris et al. (2005), Agrama et al. (2007), Jin et al. (2010), and Vanniarajan et al. (2012) but lower than those reported by Thomson et al. (2007), Borba et al. (2009), Li et al. (2012), and Dang et al. (2014, 2015). The average polymorphic information content (PIC) value obtained in this study was 0.668, which is the same as that reported by Thomson et al. (2007) and higher than those obtained in all previously published studies of rice populations (Ordonez et al., 2010; Vanniarajan et al., 2012) with the exception of those reported by Li et al. (2012), Dang et al. (2014, 2015) and Borba et al. (2009) (0.71, 0.71, and 0.75, respectively). The wide range of genetic diversity observed in this study among the accessions with a broad geographical origin makes the set of accessions included in this study one of the best collections for mining valuable genes in rice.
The studies conducted by Olsen et al. (2006), Mather et al. (2007) and Rakshit et al. (2007) using DNA sequences indicated that the LD decays at 1 cM or less in rice. In contrast, other studies using SSR markers, such as those performed by Agrama et al. (2007), Vanniarajan et al. (2012) and Dang et al. (2014, 2015), indicated that LD decays at 20–30 cM, 20–30 cM, 10–80 cM, and 10–30 cM, respectively. The analysis of the LD levels of Groups 1 to 7 in this study revealed values that were similar to those reported by Vanniarajan et al. (2012) and Dang et al. (2015). These results suggest that the extent of LD varies by genomic region (Mather et al., 2007), rice accession (Agrama and Eizenga, 2008), and marker. Moreover, Group 4 exhibited the fastest decay velocity, followed by Group 1, and Groups 2, 3, 5, 6, and 7 exhibited the lowest decay velocities in this study (Figure 4). The fast LD decay of Group 4 is mainly attributable to the high outcrossing rate of landraces. The global high LD in cultivated rice is mainly attributable to selfing and strong bottleneck (Mather et al., 2007; Zhu et al., 2007).
The FST values of the non-admixed individuals from the seven subgroups are exceedingly high (0.672 by AMOVA) compared with those observed in large-scale samples of other species, such as Oryza rufipogon (0.1, Huang et al., 2012), Arabidopsis thaliana (0.2, Long et al., 2013) and Setaria viridis (0.49, Huang et al., 2014), indicating substantial differentiation among subgroups, and genomic differentiation among species. The high value of FST for the seven subgroups in our study identified large differences between the accessions (Table 4). The markers with higher FST values had the greatest resolution power and produced more consistent genetic distance estimates, which is consistent with the results obtained from an analysis of rice conducted by Agrama et al. (2007) and a study of a human population performed by Watkins et al. (2003). The significant FST values suggest a real difference among the subgroups, and these differences might be used to predict heterotic crosses for improving yield in hybrid breeding.
We also found that the marker RM7598_Chr2, which co-associated with the STL and SNBPL, is located in the region (20,044,821—35,072,135 bp) in which a QTL for stigma exsertion was identified by Hu et al. (2009). The marker RM280_Chr4, which co-associated with the STL, SBPL, and SNBPL, is in the region (17,689,612—26,137,600 bp) in which a QTL for stigma exsertion was identified by Yamamoto et al. (2003). These results confirm the close relationship between stigma traits (STL, SBPL, and SNBPL) and stigma exsertion.
In addition, certain loci were mapped close to gene resolution, e.g., RM282 was close to GL3 (Wang et al., 2012), indicating that association analyses of rice accessions can provide an effective approach for gene identification.
Using the entire population to mine elite alleles can maintain information integrity (Dang et al., 2015). In this study, 11 elite alleles for the STL were mined at the six identified loci. Among these alleles, 18.2% were carried by accessions collected from Northeastern China, 27.3% were carried by accessions from Central China, and 54.5% were carried by accessions from Vietnam. Similarly, certain unique elite alleles for the SBPL, SNBPL, and GL were identified in various accessions. These results suggested that during the process of rice evolution from southern to northern latitudes, the loss of certain SSR alleles is likely caused by domestication bottleneck, whereas others were retained or appear for the first time in modern cultivars. For example, RM7598-115 bp is common among Vietnam accessions but is not found in Northeastern China accessions, whereas the allele RM7598-95 bp is found only in Northeastern China accessions.
For the STL trait, the broad-sense heritability averaged over the 2 years was 93%. Among the six SSR-associated markers detected for the STL, RM5753_Chr 6 exhibited the largest PVE (12.5% in 2013 and 8.9% in 2014). Of the two elite alleles found at this marker locus, RM5753-200 bp presented the largest phenotypic effect value (0.162 mm). This elite allele was carried by 57 accessions, and Yuedao 32 was the typical carrier material. Thus, the crosses described in Table S6 could significantly improve the STL.
For the SBPL trait, the broad-sense heritability averaged over the 2 years was 93%. Of the two SSR-associated markers detected for the SBPL, RM136_Chr 6 had the largest PVE (6.8% in 2013 and 6.1% in 2014). Among the four elite alleles found at this marker locus, RM136-200 bp presented the largest phenotypic effect value (0.111 mm). This elite allele was carried by 16 accessions, and Yuexiangzhan was the typical carrier material. Thus, the crosses described in Table S6 could significantly improve the SBPL.
The broad-sense heritability for the SNBPL trait averaged over the 2 years was 89%, which is also high. Among the seven SSR markers associated with the SNBPL, RM6327_Chr 11 had the largest PVE (12.9% in 2013 and 8.1% in 2014), and two elite alleles—RM6327-180 bp and RM6327-200 bp—were found at this marker locus. Thus, the SNBPL might be improved by the crosses listed in Table S6.
For the GL, the broad-sense heritability averaged over the 2 years was 98%, which is a considerably high value. Thus, the expected improvements in the GL could be obtained by marker-assisted selection. Among the 13 SSR-associated markers detected for the GL, RM6327_Chr 11 presented the largest PVE (21.8% in 2013 and 22.9% in 2014). Of the two elite alleles found at this marker locus, RM6327-190 bp exhibited the largest phenotypic effect value (0.787 mm). This elite allele was carried by 19 accessions, and Yuexiangzhan was the typical carrier material. Thus, the crosses described in Table S6 could significantly improve the GL. Because the STL was positively correlated with the GL, increases in the GL will also improve the STL. In addition, the STL could also be improved independently from the GL because the determinant coefficient between the STL and the GL was approximately 0.5.
If the target trait must be improved further, the best elite alleles could be pyramided into one cultivar through multi-round crossing. For example, 11 elite alleles were detected for STL, and the six best elite alleles could be pyramided or substituted by the combination of the accessions Yuzhenxiang, Yuedao 32, Yuedao 100, Yuexiangzhan, and Nongxiang 18 (Table S6).
DH planned and designed the research; XD, QL, CB, and YL. Performed the field experiment; XD, EL, QL, and CB. Conducted the molecular experiment; XD and EL. Analyzed the data; XD. Wrote the manuscript; and DH. Revised the manuscript. All authors read and approved the manuscript.
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.
Funding support was provided by a grant from the China national “863” programme (2010AA101301), a grant from the Doctoral Fund of the Educational Ministry (B0201300662) and a grant from the National Natural Science Foundation of China (31571743).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01188
Figure S1. Correlation of stigma traits between 2013 and 2014. A total of 227 rice accessions were grown in the same field in 2013 and 2014 growing season. (A) Correlation of STL between 2013 and 2014. (B) Correlation of SBPL between 2013 and 2014. (C) Correlation of SNBPL between 2013 and 2014. **Significantly correlation at the P < 0.01 (two-tailed).
Figure S2. Changes of the mean LnP (K) (A) and ΔK (B). (A) A graph with mean LnP (K) on Y-axis and number of subgroups on X-axis; (B) A graph with ΔK on Y-axis and number of subgroups on X-axis.
Figure S3. Population structure of 227 Oryza Sativa. (A) STRUCTURE result. The panel is a result of K = 8 (highest likelihood result among replicates). Each accession is represented by a vertical bar. The colored subsections within each vertical bar indicate membership coefficient (Q) of the accession to different clusters. Identified subpopulations are Group1 (red color), Group2 (green color), Group3 (navy blue color), Group4 (yellow color), Group5 (purple color), Group6 (light blue color), Group7 (brown color) and Group8 (caramel color). Accessions marked by group name are considered non-admixed (more than 0.9 assignment to one group in STRUCTURE analysis with K = 8). (B) Result of principal component analysis. Each point corresponds to an individual, and the different color slices correspond to the group assignment matrix (Q matrix) of the STRUCTURE result for K = 8. Numbers in parentheses beside each axis denote the amount of variance explained by that axis. (C) Neighbor-joining tree. The different colored groups approximately correspond to Groups 1–8. Colored branches represent the non-admixed individuals within each corresponding group. Gray branches respent admixed individuals (< 0.9 assignment assigned to any group in STRUCTURE analysis with K = 8).
Figure S4. Neighbor-joining tree constructed from Nei's (1983) genetic distance of 249 SSRs. Each branch is corresponding to each accession ID.
Figure S5. Rice stigma morphology of the 6 elite parents included in the elite crosses predicted.
Supp. Table S1. Accessions used and corresponding values of the Q matrix in this study.
Supp. Table S2. Summary statistics for the 249 SSR markers used in this study.
Supp. Table S3. Locus-by-locus AMOVA for O. sativa based on Structure groups.
Supp. Table S4. Positive elite alleles, phenotypic effect value and typical materials for stigma traits and grain length.
Supp. Table S5. Parental combinations and numbers of elite alleles after combinations predicted from association mapping of stigma traits and grain length.
Supp. Table S6. Elite alleles carried by the superior parents for stigma traits and grain length and corresponding phenotypic effect.
Supp. Table S7. The list for QTLs identified from this study and shared in previous studies.
Agrama, H. A., and Eizenga, G. C. (2008). Molecular diversity and genome-wide linkage disequilibrium patterns in a worldwide collection of Oryza sativa and its wild relatives. Euphytica 160, 339–355. doi: 10.1007/s10681-007-9535-y
Aluko, G., Martinez, C., Tohme, J., Castano, C., Bergman, C., and Oard, H. (2004). QTL mapping of grain quality traits from the interspecific cross Oryza sativa × O. glaberrima. Theor. Appl. Genet. 109, 630–639. doi: 10.1007/s00122-004-1668-y
Bai, X. F., Luo, L., Yan, W., Kovi, M. R., Zhan, W., and Xing, Y. (2010). Genetic dissection of rice grain shape using a recombinant inbred line population derived from two contrasting parents and fine mapping a pleiotropic quantitative trait locus qGL7. BMC Genet. 11:16. doi: 10.1186/1471-2156-11-16. Available online at: http://www.biomedcentral.com/1471-2156/11/16
Borba, T. C. O., Brondani, R. P. V., Rangel, P. H. N., and Brondani, C. (2009). Microsatellite marker-mediated analysis of the EMBRPA rice core collection genetic diversity. Genetica 137, 293–304. doi: 10.1007/s10709-009-9380-0
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 2, 2633–2635. doi: 10.1093/bioinformatics/btm308
Breseghello, F., and Sorrells, M. E. (2006). Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 172, 1165–1177. doi: 10.1534/genetics.105.044586
Caicedo, A. L., Williamson, S. H., Hernandez, R. D., Boyko, A., Fledel-Alon, A., York, T. L., et al. (2007). Genome-wide patterns of nucleotide polymorphism in domesticated rice. PLoS Genet. 3, 1745–1756. doi: 10.1371/journal.pgen.0030163
Cho, Y. G., Ishii, T., Temnykh, S., Chen, X., Lipovich, L., McCouch, S. R., et al. (2000). Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L.). Theor. Appl. Genet. 100, 713–722. doi: 10.1007/s001220051343
Cui, D., Xu, C. Y., Tang, C. F., Yang, C. G., Yu, T. Q., A, X. X., et al. (2013). Genetic structure and association mapping of cold tolerance in improved japonica rice germplasm at the booting stage. Euphytica 193, 369–382. doi: 10.1007/s10681-013-0935-x
Dang, X. J., Thi, T. G., Dong, G. S., Wang, H., Edzes, W. M., and Hong, D. L. (2014). Genetic diversity and association mapping of seed vigor in rice (Oryza sativa L.). Planta 239, 1309–1319. doi: 10.1007/s00425-014-2060-z
Dang, X. J., Thi, T. G., Edzes, W. M., Liang, L. J., Liu, Q. M., Liu, E. B., et al. (2015). Population genetic structure of Oryza sativa in East and Southeast Asia and the discovery of elite alleles for grain traits. Sci. Rep. 5:11254. doi: 10.1038/srep11254
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
Excoffier, L., and Lischer, H. E. L. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x
Falush, D., Stephens, M., and Pritchard, J. K. (2007). Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol. Ecol. Notes 7, 574–578. doi: 10.1111/j.1471-8286.2007.01758.x
Farnir, F., Coppieters, W., Arranz, J. J., Berzi, P., Cambisano, N., Grisart, B., et al. (2000). Extensive genome-wide linkage disequilibrium in cattle. Genome Res. 10, 220–227. doi: 10.1101/gr.10.2.220
Gross, B. L., Skare, K. J., and Olsen, K. M. (2009). Novel Phr1 mutations and the evolution of phenol reaction variation in US weedy rice (Oryza sativa L.). New Phytol. 184, 842–850. doi: 10.1111/j.1469-8137.2009.02957.x
Hardy, O., and Vekemans, X. (2002). SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol. Ecol. Notes 2, 618–620. doi: 10.1046/j.1471-8286.2002.00305.x
Hu, S., Zhou, Y., Zhang, L., Zhu, X., Wang, Z., Li, L., et al. (2009). QTL analysis of floral traits of rice (Oryza sativa L.) under well-watered and drought stress conditions. Genes Genom. 31, 173–181. doi: 10.1007/BF03191150
Huang, P., Molina, J., Flowers, J. M., Rubinstein, S., Jackson, S. A., Purugganan, M. D., et al. (2012). Phylogeography of Asian wild rice, Oryza rufipogon: a genome-wide view. Mol. Ecol. 21, 4593–4604. doi: 10.1111/j.1365-294x.2012.05625.x
Huang, X. H., Zhao, Y., Wei, X. H., Li, C. Y., Wang, A. H., Zhao, Q., et al. (2011). Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat. Genet. 44, 32–39. doi: 10.1038/ng.1018
Jain, S., Jain, R. K., and McCouch, S. R. (2004). Genetic analysis of Indian aromatic and quality rice (Oryza sativa L.) germplasm using panels of fluorescently-labeled microsatellite markers. Theor. Appl. Genet. 109, 965–977. doi: 10.1007/s00122-004-1700-2
Jin, L., Lu, Y., Xiao, P., Sun, M., Corke, H., and Bao, J. S. (2010). Genetic diversity and population structure of a diverse set of rice germplasm for association mapping. Theor. Appl. Genet. 121, 475–487. doi: 10.1007/s00122-010-1324-7
Kato, H., and Namai, H. (1987a). Intervarietal variations of floral characteristics with special reference to F1 seed production in Japonica rice (Oryza sativa L.). Jpn. J. Breed. 37, 75–87. doi: 10.1270/jsbbs1951.37.75
Kato, H., and Namai, H. (1987b). Floral characteristics and environmental factors for increasing natural outcrossing rate for F1hybrid seed production of rice (Oryza sativa L.). Jpn. J. Breed. 37, 318–330. doi: 10.1270/jsbbs1951.37.318
Kumagai, M., Wang, L., and Ueda, S. (2010). Genetic diversity and evolutionary relationships in genus Oryza revealed by using highly variable regions of chloroplast DNA. Gene 462, 44–51. doi: 10.1016/j.gene.2010.04.013
Li, J., Xiao, J., Grandillo, S., Jiang, L., Wan, Y., Deng, Q., et al. (2004). QTL detection for rice grain quality traits using an interspecific backcross population derived from cultivated Asian (O. sativa L.) and African (O. glaberrima S.) rice. Genome 47, 697–704. doi: 10.1139/g04-029
Li, P. B., Feng, F. C., Zhang, Q. L., Chao, Y., Gao, G. J., and He, Y. Q. (2014). Genetic mapping and validation of quantitative trait loci for stigma exsertion rate in rice. Mol. Breed. 34, 2131–2138. doi: 10.1007/s11032-014-0168-2
Li, X. B., Yan, W. G., Agrama, H. S., Jia, L. M., Jackson, A., Moldenhauer, K., et al. (2012). Unraveling the complex trait of harvest index with association mapping in rice (Oryza sativa L.). PLoS ONE 7:e29350. doi: 10.1371/journal.pone.0029350
Liu, Q. M., Qin, J. C., Li, T. W., Liu, E. B., Fan, D. J., Edzesi, W. M., et al. (2015). Fine mapping and candidate gene analysis of qSTL3, a stigma length-conditioning locus in rice (Oryza sativa L.). PLoS ONE 10:e0127938. doi: 10.1371/journal.pone.0127938
Long, Q., Rabanal, F. A., Meng, D., Huber, C. D., Farlow, A., Platzer, A., et al. (2013). Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden. Nat. Genet. 45, 884–890. doi: 10.1038/ng.2678
Lu, Z. M., and Hong, D. L. (1999). “Advances in hybrid rice seed production techniques,” in Heterosis and Hybrid Seed Production in Agronomic Crops, ed A. S. Basra (New York, NY: Food Products Press, an imprint of the Haworth Press, Inc.), 65–79.
Mather, K. A., Caicedo, A. L., Polato, N. R., Olsen, K. M., McCouch, S., and Purugganan, M. D. (2007). The extent of linkage disequilibrium in rice (Oryza sativa L.). Genetics 177, 2223–2232. doi: 10.1534/genetics.107.079616
McCouch, S. R., Teytelman, L., Xu, Y. B., Lobos, K. B., Clare, K., Walton, M., et al. (2002). Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.). DNA Res. 9, 199–207. doi: 10.1093/dnares/9.6.199
Miyata, M., Yamamoto, T., Komori, T., and Nitta, N. (2007). Marker-assisted selection and evaluation of the QTL for stigma exsertion under japonica rice genetic background. Theor. Appl. Genet. 114, 539–548. doi: 10.1007/s00122-006-0454-4
Molina, J., Sikora, M., Garud, N., Flowers, J. M., Rubinstein, S., Reynolds, A., et al. (2011). Molecular evidence for a single evolutionary origin of domesticated rice. Pro. Natl. Acad. Sci. U.S.A. 108, 8351–8356. doi: 10.1073/pnas.1104686108
Olsen, K. M., Caicedo, A. L., Polato, N. R., McClung, A., McCouch, S., and Purugganan, M. D. (2006). Selection under domestication: evidence for a sweep in the rice Waxy genomic region. Genetics 173, 975–983. doi: 10.1534/genetics.106.056473
Ordonez, S. A. Jr. Silva, J., and Oard, J. H. (2010). Association mapping of grain quality and flowering time in elite japonica rice germplasm. J. Cereal Sci. 51, 337–343. doi: 10.1016/j.jcs.2010.02.001
R Development Core Team (2011). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org
Rakshit, S., Rakshit, A., Matsumura, H., Takahashi, Y., Hasegawa, Y., Ito, A., et al. (2007). Large-scale DNA polymorphism study of Oryza sativa and O. rufipogon reveals the origin and divergence of Asian rice. Theor. Appl. Genet. 114, 731–743. doi: 10.1007/s00122-006-0473-1
Rebolledo, M. C., Dingkuhn, M., Courtois, B., Gibon, Y., Clément-Vidal, A., Cruz, D. F., et al. (2015). Phenotypic and genetic dissection of component traits for early vigour in rice using plant growth modelling, sugar content analyses and association mapping. J. Exp. Bot. 18, 5555–5566. doi: 10.1093/jxb/erv258
Stacklies, W., Redestig, H., Scholz, M., Walther, D., and Selbig, J. (2007). pcaMethods a bioconductor package providing PCA methods for incomplete data. Bioinformatics 23, 1164–1167. doi: 10.1093/bioinformatics/btm069
Tan, Y. F., Xing, Y. Z., Li, J. X., Yu, S. B., Xu, C. G., and Zhang, Q. F. (2000). Genetic bases of appearance quality of rice grains in Shanyou 63, an elite rice hybrid. Theor. Appl. Genet. 101, 823–829. doi: 10.1007/s001220051549
Temnykh, S., Park, W. D., Ayres, N., Cartinhour, S., Hauck, N., Lipovich, L., et al. (2000). Mapping and genome organization of microsatellite sequence in rice (Oryza sativa L.). Theor. Appl. Genet. 100, 697–712. doi: 10.1007/s001220051342
Thomson, M. J., Septiingsih, E. M., Suwardjo, F., Santoso, T. J., Silitonga, T. S., and McCouch, S. R. (2007). Genetic diversityanalysis of traditional and improved Indonesian rice (Oryza sativa L.) germplasm using microsatellite markers. Theor. Appl. Genet. 114, 559–568. doi: 10.1007/s00122-006-0457-1
Uga, Y., Fukuta, Y., Cai, H. W., Iwata, H., Ohsawa, R., Morishima, H., et al. (2003a). Mapping QTLs influencing rice floral morphology using recombinant inbred lines derived from a cross beween Oryza sativa L. and Oryza rufipogon Griff. Theor. Appl. Genet. 107, 218–226. doi: 10.1007/s00122-003-1227-y
Uga, Y., Fukuta, Y., Ohsawa, R., and Fujimura, T. (2003b). Variations of floral traits in Asian cultivated rice (Oryza sativa L.) and its wild relatives (O rufipogon Griff). Breed. Sci. 53, 345–352. doi: 10.1270/jsbbs.53.345
Uga, Y., Siangliw, M., Nagamine, T., Ohsawa, R., Fujimura, T., and Fukuta, Y. (2010). Comparative mapping of QTLs determining glume, pistil and stamen sizes in cultivated rice (Oryza sativa L.). Plant Breed. 129, 657–669. doi: 10.1111/j.1439-0523.2009.01765.x
Virmani, S. S. (1994). “Heterosis and hybrid rice breeding,” in Monographs on Theoretical and Applied Genetics eds R. Frankel, B. M. Grossman, W. H. F. Linskens, N. P. Maliga, and P. R. Riley (Berlin, Heidelberg: Springer-Verlag), 189.
Wang, L., Wang, A., Huang, X., Zhao, Q., Dong, G., Qian, Q., et al. (2011). Mapping 49 quantitative trait loci at high resolution through sequencing-based genotyping of rice recombinant inbred lines. Theor. Appl. Genet. 122, 327–340. doi: 10.1007/s00122-010-1449-8
Watkins, W. S., Rogers, A. R., Ostler, C. T., Wooding, S., Bamshad, M. J., Brassington, A. E., et al. (2003). Genetic variation among world populations: inferences from 100 Alu insertion polymorphisms. Genome Res. 13, 1607–1618. doi: 10.1101/gr.894603
Xie, F., and Hardy, B. (2014). “Public-private partnership for hybrid rice,” in Progress in Rice Breeding and Production in China, Proceedings of the 6th International Hybrid Rice Symposium, eds L. Cao and P. Yu (Los Baños, CA: International Rice Research Institute), 9–17.
Yamamoto, T., Takemori, N., Sue, N., and Nitta, N. (2003). QTL analysis of stigma exsertion in rice. Rice Genet. Newsl. 20, 33–34. Available online at: http://archive.gramene.org/newsletters/rice_genetics/rgn20/b13.html
Yan, W. G., Li, Y., Agrama, H. A., Luo, D., Gao, F., Lu, X., et al. (2009). Association mapping of stigma and spikelet characteristics in rice (Oryza sativa L.). Mol. Breed. 24, 277–292. doi: 10.1007/s11032-009-9290-y
Yang, W. N., Guo, Z. L., Huang, C. L., Wang, K., Jiang, N., Feng, H., et al. (2015). Genome-wide association study of rice (Oryza sativa L.) leaf traits with a high-throughput leaf scorer. J. Exp. Bot. 18, 5605–5615. doi: 10.1093/jxb/erv100
Yu, X. Q., Mei, H. W., Luo, L. J., Liu, G. L., Zou, G. H., Hu, S. P., et al. (2006). Dissection of additive, epistatic and Q × E interaction of quantitative trait loci influencing stigma exsertion under water stress in rice. Acta Genet. Sin. 33, 542–550. doi: 10.1016/S0379-4172(06)60083-8
Zhao, K. Y., Tung, C. W., Eizenga, G. C., Wright, M. H., Ali, M. L., Price, A. H., et al. (2011). Genome-wide association mapping reveals a rich genetic arachitecture of complex traits in Oryza sativa. Nat. Commun. 2, 467. doi: 10.1038/ncomms1467
Zhu, Q., Zheng, X., Luo, J., Gaut, B. S., and Ge, S. (2007). Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice. Mol. Biol. Evol. 24, 875–888. doi: 10.1093/molbev/msm005
Keywords: QTL detection, elite allele, genome-wide association mapping, Oryza sativa, stigma traits
Citation: Dang X, Liu E, Liang Y, Liu Q, Breria CM and Hong D (2016) QTL Detection and Elite Alleles Mining for Stigma Traits in Oryza sativa by Association Mapping. Front. Plant Sci. 7:1188. doi: 10.3389/fpls.2016.01188
Received: 20 May 2016; Accepted: 22 July 2016;
Published: 09 August 2016.
Edited by:Tian Tang, Sun Yat-sen University, China
Reviewed by:Jinfeng Chen, University of California, Riverside, USA
Guiquan Zhang, South China Agricultural University, China
Copyright © 2016 Dang, Liu, Liang, Liu, Breria and Hong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Delin Hong, firstname.lastname@example.org