Original Research ARTICLE
Favorable Marker Alleles for Panicle Exsertion Length in Rice (Oryza sativa L.) Mined by Association Mapping and the RSTEP-LRT Method
- State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, China
The panicle exsertion length (PEL) in rice (Oryza sativa L.) is an important trait for hybrid seed production. We investigated the PEL in a chromosome segment substitution line (CSSL) population consisting of 66 lines and a natural population composed of 540 varieties. In the CSSL population, a total of seven QTLs for PEL were detected across two environments. The percentage of phenotypic variance explained (PVE) ranged from 10.22 to 50.18%, and the additive effect ranged from −1.77 to 6.47 cm. Among the seven QTLs, qPEL10.2 had the largest PVE, 44.05 and 50.18%, with an additive effect of 5.91 and 6.47 cm in 2015 and in 2016, respectively. In the natural population, 13 SSR marker loci were detected that were associated with PEL in all four environments, with the PVE ranging from 1.20 to 6.26%. Among the 13 loci, 7 were novel. The RM5746-170 bp allele had the largest phenotypic effect (5.11 cm), and the typical carrier variety was Qiaobinghuang. An RM5620-RM6100 region harboring the EUI2 locus on chromosome 10 was detected in both populations. The sequencing results showed that the accessions with a shorter PEL contained the A base, while the accessions with a longer PEL contained the G base at the 1,475 bp location of the EUI2 gene.
Rice (Oryza sativa L.), one of the most important food crops, feeds more than half of the world's population (Sasaki and Burr, 2000). With rapid global population growth and the decreasing availability of arable land, food security has encountered great challenges, and the demand for a high-yielding rice cultivar is a high priority for breeders. Utilization of heterosis is considered one of the most effective strategies for increasing rice yield. To date, hybrid rice technology has been adopted in 27 rice-growing countries in the world (Xie, 2009). In China, hybrid rice has been planted in up to 50% of the total area each year since 1985 (Yuan and Virmani, 1988). The rice hybrid seed production area is ~150,000 hectares annually (Lu and Hong, 1999; Cheng et al., 2004).
In the process of hybrid rice seed production (Figure 1), cytoplasmic male-sterile (CMS) lines often have a defect of panicle enclosure (Figure 2), which greatly reduces the number of spikelets that can be pollinated by male parents (Shen et al., 1987; Yang et al., 2002). Panicle enclosure in CMS lines is caused by a reduction of the indigenous hormone GA3 supply to the panicle, resulting in restricted elongation of the upper internodes. Spraying exogenous gibberellin (GA3) on the CMS plants at the initial heading stage in the F1 seed production field can eliminate the panicle enclosure. However, spraying a large amount of GA3 increases the seed production cost, environmental pollution, and the occurrence of rice kernel smut (Leung et al., 2003; Tsuda et al., 2006; Brooks et al., 2009; Chen et al., 2016). Rutger and Carnahan (1981) found a recessive rice mutant (No. 76:4512) with an upper internode elongation trait derived from the progenies of a japonica rice hybrid and designated it elongated uppermost internode (eui) mutant. The mutant shows a dramatically elongated upper internode and is an important germplasm for developing unsheathed-panicle male-sterile lines in hybrid rice breeding. Although the sheathed-panicle in the CMS line is an undesirable phenomenon induced by the sterile cytoplasm, if the panicle exsertion length (PEL) of the maintainer line were increased, panicle enclosure of the CMS line derived from the maintainer line would be alleviated. Therefore, it is important to identify and evaluate more genetic resources and favorable marker alleles of the PEL in different germplasms for hybrid rice seed production.
Figure 2. Phenotypes of the maintainer line Zhen 257B, showing no panicle enclosure, and the male-sterile line Zhen 257A, showing panicle enclosure. (A) Plants of Zhen 257A and Zhen 257B at the filling stage on the plot. Individual plants of Zhen 257A (B) and Zhen 257B (C) at the filling stage. Scale bar, 20 cm.
Following the report by Rutger and Carnahan (1981), several long panicle exsertion germplasms, such as mh-1, Grlc, 02428h, and Xieqingzao eB2, were identified (Wu and Zhang, 1983; Liao et al., 1988; Li et al., 1992; Yang, 1998). Genetic and functional analyses revealed that EUI1 was located on the long arm of chromosome 5 and encodes a putative cytochrome P450 monooxygenase (Ma et al., 2004, 2006; Luo et al., 2006; Zhu et al., 2006; Magome et al., 2013). In addition to EUI1, several other genes related to upmost internode elongation have been cloned and their functions elucidated. EUI2 encodes a homologous protein of epoxide hydrolase (Zhu, 2003). Four SUI (shortened upmost internode) family genes—SUI1, SUI2, SUI3, and SUI4—encode base-exchange types of phosphatidylserine synthases (Zhu et al., 2011; Yin et al., 2013; Ji et al., 2014). Further, 33 QTLs associated with PEL have been identified using bi-parental segregating populations. These QTL were distributed throughout all 12 chromosomes (Yamamoto et al., 2001; Hittalmani et al., 2003; Qiao et al., 2007; Xiao et al., 2008; Yang et al., 2009; Guan et al., 2011; Ji et al., 2014; Zhao et al., 2015).
QTLs mapped using bi-parental segregating populations reflect the differences of alleles between parents at these loci for the target trait. Association mapping (AM) based on linkage disequilibrium (LD) has been frequently used to search for more than two alleles at a gene locus or additional gene loci and explore the natural genomic diversity at a population level (Zhu et al., 2008; Brachi et al., 2011; Weigel, 2012). AM has been widely applied to exploit favorable marker alleles in rice for many traits, including yield 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; McCouch et al., 2016), quality traits (Borba et al., 2010; Huang et al., 2010; Jin et al., 2010; Feng et al., 2016), resistance traits (Famoso et al., 2011; Jia et al., 2012, 2016; Cui et al., 2013, seed vigor traits Cui et al., 2013; Dang et al., 2014; Rebolledo et al., 2015; Phung et al., 2016), and outcrossing traits (Yan et al., 2009; Huang et al., 2010; Dang et al., 2016). Although several studies on associations between markers and the PEL trait in rice have been reported (Begum et al., 2015; Crowell et al., 2016), no significant association has been detected between PEL and molecular markers. It is necessary to enlarge the geographical range of accessions for identifying favorable marker alleles.
The objectives of the present study were (1) to detect QTLs for PEL in a chromosome segment substitution line (CSSL) population using a likelihood ratio test based on stepwise regression (RSTEP-LRT) method; (2) to identify favorable marker alleles of PEL in a natural population composed of 540 accessions collected from 18°N to 54°N by AM; and (3) to compare the relationship between common chromosome regions controlling PEL detected by the two categories of population and the known gene locus.
Materials and Methods
A CSSL population and a natural population were used. The CSSL population consisted of 66 lines with different A7444 chromosomal segments in the II-32B genetic background (called the IIA-CSSL population hereafter). The total chromosome length of the donor parent A7444 segments in the IIA-CSSL population was ~989.93 Mb (non-idealized CSS lines carrying several substituted segments from the donor parents), covering 94.98% of the complete A7444 genome (She et al., 2017). The natural population was composed of 540 varieties and was the same as that used by Dang et al. (2014). In self-pollinated crop, a natural population is defined as a population that consists of various pure line varieties bred in various places and times within the limits of species. Among the 540 rice accessions, 419 were from China, and 121 were from Vietnam.
Field Planting and Trait Measurement
For the CSSL population, field experiments were conducted in two environments, i.e., Nanjing (32°07′N, 118°64′E) 2015 (E1) and Nanjing 2016 (E2). For the natural population, field experiments were conducted in four environments, i.e., E1, E2, Yuanyang (35°05′N, 113°97′E) 2015 (E3), and Yuanyang 2016 (E4). The seeds were sown on May 12, and the seedlings were transplanted on June 12 in the E1 and E2 environments. The corresponding dates were May 5 and June 10 in the E3 and E4 environments. The field trials followed a complete randomized block design with two replicates in each environment. Each plot contained 40 plants, with five rows and eight plants in each row. The density was 17 cm between plants, and 20 cm between rows. Standard agronomic management practices were followed.
Twenty to thirty days after the heading dates, five plants in the middle row of each plot were selected to measure the PEL of the main stem panicle. We defined the PEL as follows. If the panicle base node was below the collar of the flag leaf, the PEL was recorded as a negative value; this case is shown in Figure 2B. If the panicle base node was above the collar of the flag leaf, the PEL was recorded as a positive value; this case is shown in Figure 2C. The average value of the five measurements for each accession was used for further analysis.
DNA Extract and PCR Amplification
At the tillering stage, a fresh leaf of each accession was collected for DNA extraction. The method of DNA extraction was the same as that described by Gross et al. (2009). The protocols used for SSR amplification, separation of PCR products, and allele recording were the same as those described by Dang et al. (2014).
All basic statistical analyses were performed in SAS package (SAS Institute Inc., Cary, NC, USA). The calculation method of broad-sense heritability () was the same as that described by Dang et al. (2014).
QTL Detection in the IIA-CSSL Population
Since the IIA-CSSL population is not an idealized one, in that each line has a single segment from the donor parent, the standard t-test is not suitable for detecting QTLs. Stepwise regression was used to select the most important segments for the trait of interest, and the likelihood ratio test was used to calculate the LOD score of each chromosome segment. This method is statistically equivalent to the standard t-test with idealized CSS lines and is designated the RSTEP-LRT (likelihood ratio test based on stepwise regression) mapping method (Wang et al., 2006). The LOD score threshold was arbitrarily set to 3.0, and the QTL detection was completed using the RSTEP-LRT method in the CSL program of QTL IciMapping 4.1 software (Wang et al., 2006). Further, the additive effect and the percentage of phenotypic variance explained (PVE) of each QTL were analyzed. The QTL nomenclature followed the principle reported by McCouch (2008).
Association Analysis in Natural Population
Rare alleles, i.e., alleles with frequencies of 5% or less, were removed from the dataset before the associations of traits and markers were analyzed. The association analysis of traits and markers was performed with a mixed linear model (MLM) using TASSEL 3.0 (Bradbury et al., 2007). An MLM method can significantly reduce spurious marker-trait associations (type I errors, or false positives) resulting from the population structure because the Q and K matrices are used as covariates. The Q matrix has been reported by Dang et al. (2014), and the K matrix is shown in Table S1. A false discovery rate (FDR) of 0.05 was used as a threshold for significant associations according to the method described by Benjamini and Hochberg (1995). The “null allele” (non-amplified allele, Figure S1) was used to determine the phenotypic effects of other alleles using the identified association locus (Breseghello and Sorrells, 2006).
Sequencing analysis of the allele variation of the EUI2 gene (Gene ID: LOC_Os10g35490) was performed using 27 accessions, and the Nipponbare sequence was used as the reference sequence. Gene-specific PCR primers (forward primer: TCTTTTGCCCCATGCGTTGA; reverse primer: TTCAAACCAAAGCTTCAGCAC) were designed using Primer-BLAST from NCBI (http://www.ncbi.nlm.nih.gov/tools/primer-blast) to amplify the target DNA fragment.
Each 50 μl PCR reaction consisted of 20 mM Tris-HCl (pH 8.0), 100 mM NaCl, 1 mM ethylene diamine tetraacetic acid, 0.1% sodium dodecyl sulfate, 2 nM dNTPs, 0.14 pM forward primers, 0.14 pM reverse primers, 1 U of Super-Fidelity DNA polymerase, and 20 ng of genomic DNA. DNA amplification was performed using a PTC-100TM Peltier Thermal Cycler (MJ Research™ Incorporated, USA) under the following conditions: (1) pre-denaturation at 95°C for 3 min; (2) 35 cycles of denaturation at 95°C for 15 s, annealing at 50°C for 15 s, extension at 72°C for 2 min; and (3) final extension at 72°C for 5 min. The PCR product was gel-purified and sequenced by TsingKe Biological Technology Ltd., Nanjing, China. All sequences were performed using the Sanger double chain termination method by the 3730XL Sequencer of Applied Biosystems Company, USA. Multiple gene and amino acid sequence alignment was performed using the software DNAman 5.2.9.
Characteristics of the PEL Trait in Parents and IIA-CSSL Populations
Table 1 shows that the PEL of A7444 (3.2 cm) was significantly longer than that of II-32B (2.0 cm) in both environments. The mean value of the IIA-CSSL population is 2.9 cm, ranging from −3.1 to 12.2 cm. Bidirectional transgressive segregation of the PEL trait was observed in the IIA-CSSL population (Figure S2). The PEL of CSSL64 was −3.3 cm, which was shorter than that of II-32B, and the PEL of CSSL57 was 10.4 cm, which was longer than that of A7444. The results of variance analysis indicated highly significant differences in the PEL trait among the 66 CSSL lines.
QTL Detection in the IIA-CSSL Population
Seven QTLs for PEL were detected on chromosomes 1, 8, 10, 11, and 12 in the two environments (Table 2; Figure S3). The PVE ranged from 10.22 to 50.18%, and the additive effect ranged from −1.77 to 6.47 cm. The favorable marker alleles at qPEL1.1, qPEL8, and qPEL10.1 loci were from II-32B, and the favorable marker alleles at qPEL1.2, qPEL10.2, qPEL11, and qPEL12 loci were from A7444. The QTLs qPEL1.2 and qPEL10.2 were identified in both E1 and E2. The QTL qPEL10.2 had the largest PVE, and the PVE was 44.05 and 50.18%, respectively, in E1 and E2. In addition, the additive effect was 5.91 and 6.47 cm in E1 and E2, respectively. The average additive effect over two environments was 6.19 cm, suggesting that the PEL of II-32B would increase 6.19 cm if the two alleles at the qPEL10.2 locus in II-32B were replaced by the two alleles in A7444.
Table 2. QTLs detected for the PEL trait and their genetic parameters in IIA-CSSL population across two environments.
Phenotypic Variations in the Natural Population
The average values over the 540 rice varieties for the PEL were 6.4, 6.1, 6.1, and 6.1 cm for E1, E2, E3, and E4, respectively (Table 3). The rice panicles exhibiting the different values of PEL are shown in Figure 3. The shortest PEL (−5.8 cm) in four environments was found in the variety Hongyin 1009 (Figure 3A). The PELs of Nannongjing 62401 (Figure 3B), Wujing 68 (Figure 3C), and Xudao 2hao (Figure 3D) were 1.4, 4.8, and 10.5 cm, respectively. The largest average value (16.6 cm) of the PEL in four environments was found in the variety Shenlenuo (Figure 3E). A normal distribution for the PEL trait was observed based on the skewness and kurtosis statistics (Table 3, Figure S4). The broad-sense heritability value averaged over four environments was 92.4% for the PEL trait.
Table 3. Phenotypic characteristics of the PEL (cm) based on 540 rice accessions in natural population across four environments.
Figure 3. Morphology of the panicle showing the PEL in rice. (A) Panicle of accession Hongyin 1009, showing a minimum value of the PEL. (B) Panicle of accession Nannongjing 62401. (C) Panicle of accession Wujing 68. (D) Panicle of accession Xudao 2hao. (E) Panicle of accession Shenlenuo, showing a maximum value of the PEL. Scale bar, 3 cm.
SSR Loci Associated with PEL in the Natural Population
To reduce the spurious marker-trait associations, we deleted the accessions with Q-value larger than 0.9. Then, the 483 accessions were used to detect the association between the marker and the trait by MLM, in which the Q and K matrices were used as covariates. Thirteen significant association loci were detected in all four environments by MLM under an FDR at the 0.05 level (Table 4). The 13 marker loci were located on chromosomes 1, 2, 3, 5, 6, 8, 9, 10, and 12. The PVE ranged from 1.20 to 6.26%. RM6811 on chromosome 6 had the highest PVE for PEL, explaining phenotypic variance of 5.36% in E1, 4.83% in E2, 6.26% in E3, and 5.69% in E4.
Table 4. Marker associated with PEL at P-values < 0.05, their equivalent false discovery rate probability (FDR), proportion of PVE, and marker position on chromosome derived from 262 markers and 483 rice accessions.
Favorable Marker Alleles for PEL and Their Typical Carriers
In this study, the marker alleles with positive effects were considered favorable marker alleles for the PEL. The larger the positive effect value, the better the marker allele. A summary of the favorable marker alleles and their typical carrier varieties is shown in Table S2. A total of 26 favorable marker alleles for the PEL were detected across the entire population. The observed allele frequencies for RM16-170 bp and RM6100-145 bp were 82.8 and 83.0%, respectively. The phenotypic effect values of RM524-195 bp, RM159-240 bp, RM269-165 bp, RM152-135 bp, and RM5746-180 bp were more than 3.0 cm on the PEL. The allele RM524-195 bp and RM152-135 bp showed 3.3 and 4.5 cm of the phenotypic effect on the PEL, respectively. In addition, the typical carrier varieties were both Zhongshuyangzhongdao. The phenotypic effect value of RM159-240 bp and RM269-165 bp was 3.6 and 4.2 cm on the PEL, respectively. The typical carrier variety for both of these was Shelenuo. The allele RM5746-180 bp showed the largest phenotypic effect (5.1 cm) on the PEL, and the typical carrier variety was Juhuahuang.
Differences in DNA Sequences at the EUI2 Gene Locus between Parents of CSSL and among Accessions with Contrasting PEL
In the IIA-CSSL population, we detected a QTL qPEL10.2 located on the region between RM5620 (17.93 Mb) and RM6100 (19.35 Mb). In the natural population, we detected a region between RM269 (17.40 Mb) and RM6100 (19.35 Mb) on chromosome 10, associated with the PEL trait. The RM5620-RM6100 region is known to harbor a gene EUI2 locus on chromosome 10. The total length of this gene is 2,114 bp, including five exons. The gene encodes a homologous protein of epoxide hydrolase with 311 amino acids (Zhu et al., 2006). In the IIA-CSSL population, we sequenced DNA of EUI2 in II-32B and A7444 and aligned it with the reference genome Nipponbare. No difference between A7444 and Nipponbare was found (Figure 4A). Seven variants were found between II-32B and Nipponbare (Figure 4A). The seven variants included six single nucleotide substitutions and one deletion of two bases (Figure 4A). Base A of the fourth exon 1,475 bp in II-32B is replaced by G in A7444, resulting in the transformation of amino acids from asparagine acid to aspartic acid (Figure 4B). The remaining mutations were synonymous. In the natural population, 16 accessions with high values of PEL carrying the RM6100-145 bp allele and nine accessions with low values of PEL carrying the RM6100-135 bp allele were sequenced for the EUI2 gene. The PEL-values of the 25 accessions are listed in Table S3. After aligning the DNA sequences of EUI2, we found that the gene sequences of the accessions with high PEL-values are the same as those of Nipponbare and that the gene sequences of the accessions with low PEL-values show seven variants, in contrast to Nipponbare, which includes six single nucleotide substitutions and one deletion of two bases (Figure S5A). The one deletion of two bases was located in a non-coding region. Among the six single nucleotide substitutions, four were located in a non-coding region and two were located in a coding region. Within the two variants in the coding region, the base A of the accessions with low PEL-values is replaced by G in the accessions with high PEL-values at the location of 1,475 bp, resulting in the transformation of amino acid from asparagine acid to aspartic acid, and the single nucleotide substitutions were synonymous at the location of 1,771 bp (Figure S5B).
Figure 4. Gene sequence and deduced peptide sequence alignment of EUI2 between two parents (A7444 and II-32B) and Nipponbare. (A) Gene sequence. We show only the alignment results of the variance sequence; the alignment results of the no-difference sequence are not shown. (B) Deduced peptide sequence.
In the IIA-CSSL population, we detected seven QTLs for PEL in two environments (Table 2, Figure S3). Among the seven QTLs detected, four QTLs, namely, qPEL1.2, qPEL10.1, qPEL10.2, and qPEL12, were located on sites near the chromosome location of QTL reported, and three QTLs, qPEL1.1, qPEL8, and qPEL11, were found to be novel after reviewing previous reports. For the four QTLs reported previously, qPEL1.2 was located on a site near the MOC2 gene cloned by Koumoto et al. (2013). The moc2 mutant showed a low level of cytosolic fructose-1,6-bisphosphatase 1. Defective fructose-1,6-bisphosphatase 1 activity leads to a shortage of sucrose supply, which probably causes the inhibition of tissue elongation. The QTL qPEL10.1 was located on a site near the TAW1 gene cloned by Yoshida et al. (2013). TAW1 encodes a nuclear protein of unknown function and shows high levels of expression in the shoot apical meristem, the inflorescence meristem and the branch meristem. The QTL qPEL10.2, which was the main QTL (44.05% PVE in E1 and 50.18% in E2), was located on a site near the EUI2 gene cloned by Zhu (2003). The QTL qPEL12, closely linked with the marker RM1246, was in the same region as that reported by Zhao et al. (2015). For the three new QTLs detected in the present study, the additive effect of QTLs qPEL1.1 and qPEL8 was a negative value. The synergistic allelic variations were from II-32B. The PVE of qPEL1.1 and qPEL8 was 18.92 and 13.94%, respectively. The QTL qPEL11 was closely linked with marker RM224 and the favorable marker allele for PEL was RM224-120 bp from parent A7444 (japonica). This favorable marker allele can be used directly for improving the PEL of the japonica rice variety by marker-assisted breeding.
Among the 13 association loci of marker-traits in the natural population identified for the PEL, three marker loci (RM7288 on chromosome 2, RM6811 on chromosome 6, and RM5746 on chromosome 12) exhibited a PVE of more than 5% over four environments (Table 4). Marker locus RM6811 had the largest PVE, with three favorable marker alleles. Allele RM6811-150 bp showed the largest phenotypic effect value (1.7 cm). Fifty of the 483 accessions carried this favorable marker allele, and Qiaobinghuang was the typical carrier variety (possessing the largest phenotypic value among the 50 accessions). Of the two favorable marker alleles found at marker locus RM7288, marker allele RM7288-170 bp had the larger phenotypic effect value (2.14 cm). Forty-seven of the accessions carried this favorable marker allele, and Xiaobaidao was the typical carrier variety. RM5746-170 bp and RM5746-180 bp were the two favorable marker alleles at the RM5746 locus. The former showed the larger phenotypic effect value (5.11 cm), with the typical carrier variety being Qiaobinghuang, and this allele was carried by 36 accessions.
Six of the 13 associations were located on the chromosomal regions where the QTLs (qPE2, qPE3-2, qPE9, qUIL3b, qEUI10, and qUIL10b) for the PEL had been identified by Yang et al. (2001), Ma et al. (2006), Qiao et al. (2007), Xiao et al. (2008) and Yang et al. (2009) (http://www.gramene.org/). These QTLs are listed in Table S4. The remaining seven association loci (RM283 on chromosome 1, RM159 on chromosome 5, RM276 on chromosome 6, RM6811 on chromosome 6, RM152 on chromosome 8, RM524 on chromosome 9, and RM5746 on chromosome 12) were newly identified in this study. For the seven new loci, 18 favorable marker alleles were identified. Among the 18 alleles, allele RM5746-170 bp showed the largest phenotypic effect value (5.11 cm). Thirty-six of the 483 accessions carried this allele, and Qiaobinghuang was the typical carrier variety (possessing the largest phenotypic value among the 36 accessions).
Marker locus RM6100, which is associated with the PEL, was detected in both the natural population and the IIA-CSSL population. In the IIA-CSSL population, RM6100-145 bp was the favorable marker allele derived from A7444. In the natural population, RM6100-145 bp was also the favorable marker allele. This marker allele was carried by 401 accessions (occupied 83%), and the typical carrier was Shenlenuo. Although the detection principle and method of the two populations for mapping QTLs are different, RM6100-145 bp can be detected consistently, indicating that the allele RM6100-145 bp is conservative.
The EUI gene encodes a previously uncharacterized cytochrome P450 monooxygenase, CYP714D1. The monooxygenase catalyzes 16α, 17-epoxidation of non-13-hydroxylated GAs. The 16α, 17-epoxidation reduces the biological activity of GA4 in rice. EUI functions as a GA-deactivating enzyme (Zhu et al., 2006). Our results of sequencing and alignment of EUI gene showed that a base A of the fourth exon 1475 bp was replaced by G, resulting in the transformation of amino acids from asparagine acid to aspartic acid. We hypothesize that the transformation of amino acids causes deactivation of cytochrome P450 monooxygenase, which hinders 16α,17-epoxidation to non-13-hydroxylated GAs, resulting in accumulation of GA4 in upmost internode cells. The accumulation of GA4 leads to elongation of the PEL.
Regarding the new loci detected in this study, it is necessary to further narrow the chromosome region and clone the favorable alleles. The mechanism of the PEL genes needs further clarification. It is also necessary to pyramid the favorable alleles into one variety by convergent crossing to breed longer PEL parents for hybrid seed production.
According to the favorable marker allele information, we predict parental combinations to improve the PEL by cross-breeding. The names of the parental combinations and the numbers of favorable alleles that we expect by pyramiding are listed in Table S5. The favorable marker alleles carried by the parents in favorable crosses and the corresponding phenotypic effects are listed in Table S6. Figure S6 shows the morphologies of the panicles of accessions Qiaobinghuang, Zhongshuyangzhongdao, Wanzhongqiu, and Yanglingdao, which are considered superior parents for improving the PEL trait.
DH: Planned and designed the research; XD, BF, DL, and OS: Performed the field experiment; XD, EL, XC, ZD, DS, GW, and YL: Conducted the molecular experiment; XD: Analyzed the data and 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.
This work was supported by a grant from the National Natural Science Foundation of China to XD (31571743) and to DH (31671658).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.02112/full#supplementary-material
Table S1. The K-matrix obtained in this study.
Table S2. Positive favorable marker alleles, phenotypic effect values, and typical varieties for the PEL.
Table S3. The list of the phenotypic value of PEL, of the 25 accessions sequenced.
Table S4. The list of common QTL identified in the natural population of this study and previous studies.
Table S5. Parental combinations and numbers of favorable marker alleles after combinations predicted from AM of the PEL.
Table S6. Favorable marker alleles carried by the superior parents for the PEL and the corresponding phenotypic effect.
Figure S1. Profile amplified by RM6215 using total DNA of partial varieties (Variety A385 to A432) M: Marker, 1-48: rice varieties: 1, Yuedao 114; 2, Yuedao 115; 3, Yuedao 116; 4, Yuedao 117; 5, Yuedao 118; 6, Yuedao 119; 7, Yuedao 120; 8, Yuedao 121; 9, Jia 45; 10, Nannongjing 3786; 11, 24248; 12, Nannongjing 4004; 13, Nannongjing 4016; 14, Zijianwujing; 15, Ningjing 2hao; 16, Wuxiang 99-8; 17, Wuyunjing 8hao; 18, Nannongjing 002; 19, Nannongjing 004; 20, Huaidao 5hao; 21, Zhongzuo 93; 22, Yandao 9hao; 23, Lianjing 4hao; 24, Jindao 1007; 25, Huajing 5hao; 26, Huajing 6hao; 27, Yangfujing 7hao; 28, Yangfujing 8hao; 29, Zhendao 99; 30, Nanjing 42; 31, Lianjing 9823; 32, Huifeng 1hao; 33, Huifeng 2hao; 34, Yandao 8hao; 35, Wuyunjing 21hao; 36, Shashani; 37, Muzhan 4hao; 38, Mudanjiang 29; 39, Mudanjiang 28; 40, Mudanjiang 27; 41, Kenzhan 2hao; 42, Heijing 8hao; 43, Hejing 1hao; 44, Beidao 4hao; 45, Beidao 3hao; 46, Suijing 12hao; 47, Songjing 12; 48, Songjing 11.
Figure S2. Frequency distribution of the PEL in the CSSL population under two growing environments. The parents of II-32B and A7444 were noted above the histograms.
Figure S3. Chromosome locations of the QTLs for the PEL trait detected in the IIA-CSSL population and other reported genes of related traits. The black symbols stand for QTLs detected in E1; the white symbols stand for QTLs detected in E2.
Figure S4. Frequency distribution of the PEL in the natural population under four growing environments.
Figure S5. Gene sequence and deduced peptide sequence alignment of EUI2 between 25 accessions and Nipponbare. (A) Gene sequence. We show only the alignment results of the variance sequence; the alignment results of the no-difference sequence are not shown. (B) Deduced peptide sequence.
Figure S6. Morphology of the panicle for four superior parents. Scale bar, 3 cm.
Begum, H., Spindel, J. E., Lalusin, A., Borromeo, T., Gregorio, G., Hemandez, J., et al. (2015). Genome-wide association mapping for yield and other agronomic traits in an elite breeding population of tropical rice (Oryza sativa L.). PLoS ONE 10:e0119873. doi: 10.1371/journal.pone.0119873
Borba, T., Brondani, R., Breseghello, F., Coelho, A., Mendonca, J., Rangel, P., et al. (2010). Association mapping for yield and grain quality traits in rice (Oryza sativa L.). Genet. Mol. Biol. 33, 515–524. doi: 10.1590/S1415-47572010005000065
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
Brooks, S. A., Anders, M. M., and Yeater, K. M. (2009). Effect of cultural management practices on the severity of false smut and kernel smut of rice. Plant Dis. 93, 1202–1208. doi: 10.1094/PDIS-93-11-1202
Chen, Y., Yang, X., Yao, J., Kyaw, E. P., Zhang, A. F., Li, Y. F., et al. (2016). Simple and rapid detection of Tilletia horrida causing rice kernel smut in rice seeds. Sci. Rep. 6:33258. doi: 10.1038/srep33258
Cheng, S. H., Cao, L. Y., Yang, S. H., and Zhai, H. Q. (2004). Forty years' development of hybrid rice: China's experience. Rice Sci. 11, 225–230. Available online at: http://www.ricescience.org
Crowell, S., Korniliev, P., Falcao, A., Ismail, A., Gregorio, G., Mezey, J., et al. (2016). Genome-wide association and high-resolution phenotyping link Oryza sativa panicle traits to numerous trait-specific QTL clusters. Nat. Commun. 7:10527. doi: 10.1038/ncomms10527
Cui, D., Xu, C. Y., Tang, C. F., Yang, C. G., Yu, T. Q., 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., Liu, E. B., Liang, Y. F., Liu, Q. M., Breria, C. M., and Hong, D. L. (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
Dang, X. J., Tran, T. G., Dong, G. S., Wang, H., Edzesi, 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., Tran, T. G., Edzesi, 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
Famoso, A. N., Zhao, K. Y., Clark, R. T., Tung, C. W., Wright, M. H., Bustamante, C., et al. (2011). Genetic architecture of aluminum tolerance in rice (Oryza sativa) determined through genome-wide association analysis ans QTL mapping. PLoS Genet. 7:e1002221. doi: 10.1371/journal.pgen.1002221
Feng, Y., Lu, Q., Zhai, R. R., Zhang, M. C., Xu, Q., Yang, Y. L., et al. (2016). Genome wide association mapping for grain shape traits in indica rice. Planta 244, 819–830. doi: 10.1007/s00425-016-2548-9
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
Guan, H. Z., Duan, Y. L., Liu, H. Q., Chen, Z. W., Zhuo, M., Zhuang, L. J., et al. (2011). Genetic ananlysis and fine mapping of an enclosed panicle mutant esp2 in rice (Oryza sativa L.). Chin. Agr. Sci. Bull 56, 1476–1480. doi: 10.1007/s11434-011-4552-9
Hittalmani, S., Huang, N., Courtois, B., Venuprasad, R., Shashidhar, H. E., Zhuang, J. Y., et al. (2003). Identification of QTL for growth- and grain yield-related traits in rice across nine locations of Asia. Theor. Appl. Genet. 107, 679–690. doi: 10.1007/s00122-003-1269-1
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
Ji, H., Kim, H., Yun, D. W., Yoon, U. H., Kim, T. H., Eun, M. Y., et al. (2014). Characterization and fine mapping of a shortened uppermost internode mutant in rice. Plant Biotechnol. Rep. 8, 125–134. doi: 10.1007/s11816-013-0280-5
Jia, L. M., Yan, W. G., Zhu, C. S., Agrama, H. A., Jackson, A., Yeater, K., et al. (2012). Allelic analysis of sheath blight resistance with association mapping in rice. PLoS ONE 7:e32703. doi: 10.1371/journal.pone.0032703
Jia, Y. L., Zhou, E. X., Lee, S., and Bianco, T. (2016). Coevolutionary dynamics of rice blast resistance gene Pi-ta and magnaporthe oryzae avirulence gen AVR-Pita 1. Phytopathology 106, 676–683. doi: 10.1094/PHYTO-02-16-0057-RVW
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
Koumoto, T., Shimada, H., Kusano, H., She, K., Iwamoto, M., and Takano, M. (2013). Rice monoculm mutation moc2, which inhibits outgrowth of the second tillers, is ascribed to lack of a fructose-1, 6-bisphosphatase. Plant Biotechnol. 30, 47–56. doi: 10.5511/plantbiotechnology.12.1210a
Leung, H., Fan, J. X., Zhu, Y. Y., Chen, H. R., Revilla-Molina, I., Pangga, I., et al. (2003). Using genetic diversity to achieve sustainable rice disease management. Plant Dis. 87, 1156–1169. doi: 10.1094/PDIS.2003.87.10.1156
Li, X. B., Yan, W. G., Agrama, H., 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
Liao, C. L., Ni, K. Y., Liu, Y. K., and Huang, Z. H. (1988). Studies on heredity and utilization of Grlc recessive tall type in rice I. Characters of plant and culm type of Grlc and its testcross F1. Southwest China J. Agr. Sci. 1, 43–46.
Lu, Z. M., and Hong, D. L. (1999). “Advances in hybrid rice seed production techniques,” in Heterosis and Hybrid Seed Producton in Agronomic Crops, ed A. S. Basra (New York, NY: Food Products Press, an imprint of the Haworth Press, Inc), 65–79.
Luo, A. D., Qian, Q., Yin, H. F., Liu, X. Q., Yin, C. X., Lan, Y., et al. (2006). EUI1, encoding a putative cytochrome P450 monooxygenase, regulates internode elongation by modulating gibberellin respomses in rice. Plant Cell Physiol. 47, 181–191. doi: 10.1093/pcp/pci233
Ma, H. L., Zhang, S. B., Ji, L., Zhu, H. B., Yang, S. L., Fang, X. J., et al. (2006). Fine mapping and in silico isolation of the EUI1 gene controlling upper internode elongation in rice. Plant Mol. Biol. 60, 87–94. doi: 10.1007/s11103-005-2762-5
Ma, H. L., Zhang, S. B., Lu, Q., Fang, X. J., Yang, S. L., and Yang, R. C. (2004). Genetic analysis and molecular mapping of eui1 gene conferring elongated uppermost internode in the india rice mutant Xieqingzao eB1. Chin. J. Agr. Biotechnol. 1, 103–107. doi: 10.1079/CJB200418
Magome, H., Nomura, T., Hanada, A., Takeda-Kamiya, N., Ohnishi, T., Shinma, Y., et al. (2013). CYP714B1 and CYP714B2 encode gibberellin 13-oxidases that reduce gibberellin activity in rice. Proc. Natl. Acad. Sci. U.S.A. 110, 1947–1952. doi: 10.1073/pnas.1215788110
McCouch, S. R., Wright, M. H., Tung, C. W., Maron, L. G., McNally, K. L., Fitzgerald, M., et al. (2016). Open access resources for genome-wide association mapping in rice. Nat. Commun. 7:10532. doi: 10.1038/ncomms10532
Phung, N. T., Mai, C. D., Thi, H. G., Truong, H., Lavarenne, J., Gonin, M., et al. (2016). Genome-wide association mapping for root traits in a panel of rice accessions from Vietnam. BMC Plant Biol. 16:64. doi: 10.1186/s12870-016-0747-y
Qiao, B. J., Wang, Y. Y., Zhu, X. B., and Hong, D. L. (2007). QTL analysis of the uppermost internode length in rice under different growing environments. Hereditas 29, 1001–1007. doi: 10.1360/yc-007-1001
Rebolledo, M. C., Dingkuhn, M., Courtois, B., Gibon, Y., Clement-Vidal, A., Cruz, D., 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
Sasaki, T., and Burr, B. (2000). International rice genome sequencing project: the effort to completely sequence the rice genome. Curr. Opin. Plant Biol. 3, 138–141. doi: 10.1016/S1369-5266(99)00047-3
She, D., Liu, Q. M., Li, D. L., Liang, Y. F., Liu, E. B., Dang, X. J., et al. (2017). QTL mapping of seven panicle traits in rice (Oryza sativa L.) using chromosome segment substitution lines derived from II-32B/A7444. Acta Agron. Sin. 43, 658–668. doi: 10.3724/sp.J.1006.2017.00658
Tsuda, M., Sasahara, M., Ohara, T., and Kato, S. (2006). Optimal application timing of simeconazole granules for control of rice kernel smut and false smut. J. Gen. Plant Pathol. 72, 301–304. doi: 10.1007/s10327-006-0288-6
Wang, J. K., Wan, X. Y., Crossa, J., Crouch, J., Weng, J. F., Zhai, H. Q., et al. (2006). QTL mapping of grain length in rice (Oryza sativa L.) using chromosome segment substitution lines. Genet. Res. 88, 93–104. doi: 10.1017/S0016672306008408
Xiao, K., Zuo, H. L., Gong, Y. J., Zhang, J. Z., Zhang, Y. J., and Dong, Y. J. (2008). Location of QTLs controlling panicle exsertion and plant height in rice. Chin. Agr. Sci. Bull 24, 95–99. doi: 10.1016/j.gene.2015.06.065
Yamamoto, T., Taguchi-Shiobara, F., Ukai, Y., Sasaki, T., and Yano, M. (2001). Mapping quantitative trait loci for days-to heading, and culm, panicle and internode lengths in a BC1F3 population using an elite rice variety, Koshihikari, as the recurrent parent. Breeding Sci. 51, 63–71. doi: 10.1270/jsbbs.51.63
Yan, W. G., Li, Y., Agrama, H., Luo, D. G., Gao, F. Y., Lu, X. J., et al. (2009). Association mapping of stigma and spikelet characteristics in rice (Oryza sativa L.). Mol. Breeding 24, 277–292. doi: 10.1007/s11032-009-9290-y
Yang, D. W., Zhu, Z., Zhang, Y. D., Lin, J., Chen, T., Zhao, L., et al. (2009). Substitution mapping of QTL for panicle exertion using CSSL in rice (Oryza sativa L.). Hereditas 7, 741–747. doi: 10.3724/SP.J.1005.2009.00741
Yang, S. L., Yang, R. C., Qu, X. P., Zhang, Q. Q., Huang, R. H., and Wang, B. (2001). Genetic and microsatellite analyses of a new elongated uppermost internode gene eui2 of rice. Acta Botanica Sin. 43, 67–71.
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
Yin, H., Gao, P., Liu, C. W., Yang, J., Liu, Z. C., and Luo, D. (2013). SUI-family genes encode phosphatidylserine synthases and regulate stem development in rice. Planta 237, 15–27. doi: 10.1007/s00425-012-1736-5
Yoshida, A., Sasao, M., Yasuno, N., Takagi, K., Daimon, Y., Chen, R., et al. (2013). TAWAWA1, a regulator of rice inflorescence architecture, function through the suppression of meristem phase transition. Proc. Natl. Acad. Sci. U.S.A. 110, 767–772. doi: 10.1073/pnas.1216151110
Zhao, C. F., Qiang, X. T., Zhang, Y. D., Zhu, Z., Chen, T., Zhao, Q. Y., et al. (2015). Genetic analysis of panicle neck length for three advanced backcross substitution lines in rice. Chin. Bull. Bot. 50, 32–39. doi: 10.3724/SP.J.1259.2015.00032
Zhao, K. Y., Tung, C. W., Eizenga, G., 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, L., Hu, J., Zhu, K. M., Fang, Y. X., Gao, Z. Y., He, Y. H., et al. (2011). Identification and characterization of SHORTENED UPPERMOST INTERNODE 1, a gene negatively regulating uppermost internode elongation in rice. Plant Mol. Biol. 77, 475–487. doi: 10.1007/s11103-011-9825-6
Zhu, Y. Y., Nomura, T., Xu, Y. H., Zhang, Y. Y., Peng, Y., Mao, B. Z., et al. (2006). ELONGATED UPPERMOST INTERNODE encodes a cytochrome P450 monooxygenase that epoxidizes gibberellins in a novel deactivation reaction in rice. Plant Cell 18, 442–456. doi: 10.1105/tpc.105.038455
Keywords: chromosome segment substitution line population, favorable marker allele, hybrid seed production, natural population, Oryza sativa, panicle exsertion length
Citation: Dang X, Fang B, Chen X, Li D, Sowadan O, Dong Z, Liu E, She D, Wu G, Liang Y and Hong D (2017) Favorable Marker Alleles for Panicle Exsertion Length in Rice (Oryza sativa L.) Mined by Association Mapping and the RSTEP-LRT Method. Front. Plant Sci. 8:2112. doi: 10.3389/fpls.2017.02112
Received: 11 July 2017; Accepted: 27 November 2017;
Published: 12 December 2017.
Edited by:Giampiero Valè, Cereal Research Centre, CRA-CER, Italy
Reviewed by:Enrico Francia, University of Modena and Reggio Emilia, Italy
Tian-Qing Zheng, Institute of Crop Sciences (CAAS), China
Copyright © 2017 Dang, Fang, Chen, Li, Sowadan, Dong, Liu, She, Wu, Liang 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, email@example.com