Favorable Marker Alleles for Panicle Exsertion Length in Rice (Oryza sativa L.) Mined by Association Mapping and the RSTEP-LRT Method

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.

Introduction

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 GA₃ supply to the panicle, resulting in restricted elongation of the upper internodes. Spraying exogenous gibberellin (GA₃) on the CMS plants at the initial heading stage in the F₁ seed production field can eliminate the panicle enclosure. However, spraying a large amount of GA₃ 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 1

Figure 1. Scene of F₁ hybrid seed production in the filling stage.

FIGURE 2

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

Plant Materials

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).

Data Analysis

All basic statistical analyses were performed in SAS package (SAS Institute Inc., Cary, NC, USA). The calculation method of broad-sense heritability ($H_B^2$) 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).

Gene Sequencing

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-100^TM 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.

Results

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.

TABLE 1

Table 1. Phenotypic values of the PEL of parents and IIA-CSSLs across two environments.

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

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

Table 3. Phenotypic characteristics of the PEL (cm) based on 540 rice accessions in natural population across four environments.

FIGURE 3

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

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

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.

Discussion

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 GA₄ 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 GA₄ in upmost internode cells. The accumulation of GA₄ 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.

Author Contributions

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.

Acknowledgments

This work was supported by a grant from the National Natural Science Foundation of China to XD (31571743) and to DH (31671658).

Supplementary Material

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.

References

Agrama, H. A., Eizenga, G. C., and Yan, W. (2007). Association mapping of yield and its components in rice cultivars. Mol. Breeding 19, 341–356. doi: 10.1007/s11032-006-9066-6

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300.

Google Scholar

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

CrossRef Full Text | Google Scholar

Brachi, B., Morris, G. P., and Borevitz, J. O. (2011). Genome-wide association studies in plants: the missing heritability is in the field. Genome Biol. 12:232 doi: 10.1186/gb-2011-12-10-232

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Garris, A. J., Tai, T. H., Coburn, J., Kresovich, S., and McCouch, S. (2005). Genetic structure and diversity in Oryza sativa L. Genetics 169, 1631–1638. doi: 10.1534/genetics.104.035642

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X. H., Wei, X. H., Sang, T., Zhao, Q., Feng, Q., Zhao, Y., et al. (2010). Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 42, 961–969. doi: 10.1038/ng.695

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Li, H. B., Sun, L. H., Zou, J. S., Li, C. G., and Chen, Z. M. (1992). Identification and study of recessive tall stalk and wide compatibility rice germplasm 02428h. Jiangsu J. Agr. Sci. 8, 48–50.

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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 F₁. 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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

McCouch, S. R. (2008). Gene nomenclature system for rice. Rice 1, 72–84. doi: 10.1007/s12284-008-9004-9

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Rutger, J. N., and Carnahan, H. L. (1981). A fourth element to facilitate hybrid cereal production – a recessive tall in rice. Crop Sci. 21, 373–376. doi: 10.2135/cropsci1981.0011183X002100030005x

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text

Shen, Z. T., Yang, C. D., and He, Z. H. (1987). Studies on eliminating panicle enclosure in WA-type MS line of rice (Oryza sativa sibsp. indica). Chin. J. Rice Sci. 1, 95–99.

Google Scholar

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

CrossRef Full Text | Google Scholar

Vanniarajan, C., Vinod, K. K., and Pereira, A. (2012). Molecular evaluation of genetic diversity and association studies in rice (Oryza sativa L.). J. Genet. 91, 1–11. doi: 10.1007/s12041-012-0146-6

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Weigel, D. (2012). Natural variation in Arabidopsis: from molecular genetics to ecological genomics. Plant Physiol. 158, 2–22. doi: 10.1104/pp.111.189845

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, S. B., and Zhang, Q. H. (1983). Study on induced mutation leading to obtaining a recessive tall-culm gene in rice. Quanzhou Institute Agr. Sci. 1, 41–45.

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

CrossRef Full Text | Google Scholar

Xie, F. (2009). “Priorities of IRRI hybrid rice breeding,” in Accelerating Hybrid Rice Development, eds F. Xie and B. Hardy (Los Baños: International Rice Research Institute), 49–62.

Google Scholar

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 BC₁F₃ population using an elite rice variety, Koshihikari, as the recurrent parent. Breeding Sci. 51, 63–71. doi: 10.1270/jsbbs.51.63

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Yang, R. C. (1998). The breeding method of male-sterile line with no panicle enclosure in rice. Bull Sci. Technol. Ach. 14, 2–3.

Yang, R. C., Zhang, S. B., Huang, R. H., Yang, S. L., and Zhang, Q. Q. (2002). Breeding technology of eui-hybrids of rice. Agr. Sci. China 1, 359–363.

Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Yuan, L., and Virmani, S. (1988). “Status of hybrid rice research and development,” in Hybrid Rice, ed W. H. Smith (Manila: International Rice Research Institute), 7–24.

Google Scholar

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

CrossRef Full Text

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, C. S., Gore, M., Buckler, E. S., and Yu, J. M. (2008). Status and prospects of association mapping in plants. Plant Genome 1, 5–20. doi: 10.3835/plantgenome2008.02.0089

CrossRef Full Text | Google Scholar

Zhu, H. B. (2003). Fine Mapping and Cloning of the Rice Elongated Uppermost Internode Gene EUI2 (t). Dissertation, Fujian Agriculture and Forestry University.

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar