PMS1T, producing phased small-interfering RNAs, regulates photoperiod-sensitive male sterility in rice.

Phased small-interfering RNAs (phasiRNAs) are a special class of small RNAs, which are generated in 21- or 24-nt intervals from transcripts of precursor RNAs. Although phasiRNAs have been found in a range of organisms, their biological functions in plants have yet to be uncovered. Here we show that phasiRNAs generated by the photopheriod-sensetive genic male sterility 1 (Pms1) locus were associated with photoperiod-sensitive male sterility (PSMS) in rice, a germplasm that started the two-line hybrid rice breeding. The Pms1 locus encodes a long-noncoding RNA PMS1T that was preferentially expressed in young panicles. PMS1T was targeted by miR2118 to produce 21-nt phasiRNAs that preferentially accumulated in the PSMS line under long-day conditions. A single nucleotide polymorphism in PMS1T nearby the miR2118 recognition site was critical for fertility change, likely leading to differential accumulation of the phasiRNAs. This result suggested possible roles of phasiRNAs in reproductive development of rice, demonstrating the potential importance of this RNA class as regulators in biological processes.

Coordinated regulation of vegetative and reproductive branching in rice.

Grasses produce tiller and panicle branching at vegetative and reproductive stages; the branching patterns largely define the diversity of grasses and constitute a major determinant for grain yield of many cereals. Here we show that a spatiotemporally coordinated gene network consisting of the MicroRNA 156 (miR156/)miR529/SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) and miR172/APETALA2 (AP2) pathways regulates tiller and panicle branching in rice. SPL genes negatively control tillering, but positively regulate inflorescence meristem and spikelet transition. Underproduction or overproduction of SPLs reduces panicle branching, but by distinct mechanisms: miR156 and miR529 fine-tune the SPL levels for optimal panicle size. miR172 regulates spikelet transition by targeting AP2-like genes, which does not affect tillering, and the AP2-like proteins play the roles by interacting with TOPLESS-related proteins (TPRs). SPLs modulate panicle branching by directly regulating the miR172/AP2 and PANICLE PHYTOMER2 (PAP2)/Rice TFL1/CEN homolog 1 (RCN1) pathways and also by integrating other regulators, most of which are not involved in tillering regulation. These findings may also have significant implications for understanding branching regulation of other grasses and for application in rice genetic improvement.

Breeding signatures of rice improvement revealed by a genomic variation map from a large germplasm collection.

Intensive rice breeding over the past 50 y has dramatically increased productivity especially in the indica subspecies, but our knowledge of the genomic changes associated with such improvement has been limited. In this study, we analyzed low-coverage sequencing data of 1,479 rice accessions from 73 countries, including landraces and modern cultivars. We identified two major subpopulations, indica I (IndI) and indica II (IndII), in the indica subspecies, which corresponded to the two putative heterotic groups resulting from independent breeding efforts. We detected 200 regions spanning 7.8% of the rice genome that had been differentially selected between IndI and IndII, and thus referred to as breeding signatures. These regions included large numbers of known functional genes and loci associated with important agronomic traits revealed by genome-wide association studies. Grain yield was positively correlated with the number of breeding signatures in a variety, suggesting that the number of breeding signatures in a line may be useful for predicting agronomic potential and the selected loci may provide targets for rice improvement.

Genetic analysis of the metabolome exemplified using a rice population.

Plant metabolites are crucial for both plant life and human nutrition. Despite recent advance in metabolomics, the genetic control of plant metabolome remains largely unknown. Here, we performed a genetic analysis of the rice metabolome that provided over 2,800 highly resolved metabolic quantitative trait loci for 900 metabolites. Distinct and overlapping accumulation patterns of metabolites were observed and complex genetic regulation of metabolism was revealed in two different tissues. We associated 24 candidate genes to various metabolic quantitative trait loci by data mining, including ones regulating important morphological traits and biological processes. The corresponding pathways were reconstructed by updating in vivo functions of previously identified and newly assigned genes. This study demonstrated a powerful tool and provided a vast amount of high-quality data for understanding the plasticity of plant metabolome, which may help bridge the gap between the genome and phenome.

Grain number, plant height, and heading date7 is a central regulator of growth, development, and stress response.

Grain number, plant height, and heading date7 (Ghd7) has been regarded as an important regulator of heading date and yield potential in rice (Oryza sativa). In this study, we investigated functions of Ghd7 in rice growth, development, and environmental response. As a long-day dependent negative regulator of heading date, the degree of phenotypic effect of Ghd7 on heading date and yield traits is quantitatively related to the transcript level and is also influenced by both environmental conditions and genetic backgrounds. Ghd7 regulates yield traits through modulating panicle branching independent of heading date. Ghd7 also regulates plasticity of tiller branching by mediating the PHYTOCHROME B-TEOSINTE BRANCHED1 pathway. Drought, abscisic acid, jasmonic acid, and high-temperature stress strongly repressed Ghd7 expression, whereas low temperature enhanced Ghd7 expression. Overexpression of Ghd7 increased drought sensitivity, whereas knock-down of Ghd7 enhanced drought tolerance. Gene chip analysis of expression profiles revealed that Ghd7 was involved in the regulation of multiple processes, including flowering time, hormone metabolism, and biotic and abiotic stresses. This study suggests that Ghd7 functions to integrate the dynamic environmental inputs with phase transition, architecture regulation, and stress response to maximize the reproductive success of the rice plant.

Differential expression of GS5 regulates grain size in rice.

Grain weight is a major determinant of grain yield. GS5 is a positive regulator of grain size such that grain width, filling, and weight are correlated with its expression level. Previous work suggested that polymorphisms of GS5 in the promoter region might be responsible for the variation in grain size. In this study, two single nucleotide polymorphisms (SNPs) between the wide-grain allele GS5-1 and the narrow-grain allele GS5-2 in the upstream region of the gene that were responsible for the differential expression in developing young panicles were identified. These two polymorphs altered the responses of the GS5 alleles to abscisic acid (ABA) treatments, resulting in higher expression of GS5-1 than of GS5-2 in developing young panicles. It was also shown that SNPs in light-responsive elements of the promoter altered the response to light induction, leading to higher expression of GS5-2 than GS5-1 in leaves. Enhanced expression of GS5 competitively inhibits the interaction between OsBAK1-7 and OsMSBP1 by occupying the extracellular leucine-rich repeat (LRR) domain of OsBAK1-7, thus preventing OsBAK1-7 from endocytosis caused by interacting with OsMSBP1, providing an explanation for the positive association between grain size and GS5 expression. These results advanced our understanding of the molecular mechanism by which GS5 controls grain size.

A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice.

Hybrid rice has greatly contributed to the global increase of rice productivity. A major component that facilitated the development of hybrids was a mutant showing photoperiod-sensitive male sterility (PSMS) with its fertility regulated by day length. Transcriptome studies have shown that large portions of the eukaryotic genomic sequences are transcribed to long noncoding RNAs (lncRNAs). However, the potential roles for only a few lncRNAs have been brought to light at present. Thus, great efforts have to be invested to understand the biological functions of lncRNAs. Here we show that a lncRNA of 1,236 bases in length, referred to as long-day-specific male-fertility-associated RNA (LDMAR), regulates PSMS in rice. We found that sufficient amount of the LDMAR transcript is required for normal pollen development of plants grown under long-day conditions. A spontaneous mutation causing a single nucleotide polymorphism (SNP) between the wild-type and mutant altered the secondary structure of LDMAR. This change brought about increased methylation in the putative promoter region of LDMAR, which reduced the transcription of LDMAR specifically under long-day conditions, resulting in premature programmed cell death (PCD) in developing anthers, thus causing PSMS. Thus, a lncRNA could directly exert a major effect on a trait like a structure gene, and a SNP could alter the function of a lncRNA similar to amino acid substitution in structural genes. Molecular elucidating of PSMS has important implications for understanding molecular mechanisms of photoperiod regulation of many biological processes and also for developing male sterile germplasms for hybrid crop breeding.

Rice WRKY13 regulates cross talk between abiotic and biotic stress signaling pathways by selective binding to different cis-elements.

Plants use a complex signal transduction network to regulate their adaptation to the ever-changing environment. Rice (Oryza sativa) WRKY13 plays a vital role in the cross talk between abiotic and biotic stress signaling pathways by suppressing abiotic stress resistance and activating disease resistance. However, it is not clear how WRKY13 directly regulates this cross talk. Here, we show that WRKY13 is a transcriptional repressor. During the rice responses to drought stress and bacterial infection, WRKY13 selectively bound to certain site- and sequence-specific cis-elements on the promoters of SNAC1 (for STRESS RESPONSIVE NO APICAL MERISTEM, ARABIDOPSIS TRANSCRIPTION ACTIVATION FACTOR1/2, CUP-SHAPED COTYLEDON), the overexpression of which increases drought resistance, and WRKY45-1, the knockout of which increases both bacterial disease and drought resistance. WRKY13 also bound to two cis-elements of its native promoter to autoregulate the balance of its gene expression in different physiological activities. WRKY13 was induced in leaf vascular tissue, where bacteria proliferate, during infection, and in guard cells, where the transcriptional factor SNAC1 enhances drought resistance, during both bacterial infection and drought stress. These results suggest that WRKY13 regulates the antagonistic cross talk between drought and disease resistance pathways by directly suppressing SNAC1 and WRKY45-1 and autoregulating its own expression via site- and sequence-specific cis-elements on the promoters of these genes in vascular tissue where bacteria proliferate and guard cells where the transcriptional factor SNAC1 mediates drought resistance by promoting stomatal closure.

An expression quantitative trait loci-guided co-expression analysis for constructing regulatory network using a rice recombinant inbred line population.

The ability to reveal the regulatory architecture of genes at the whole-genome level by constructing a regulatory network is critical for understanding the biological processes and developmental programmes of organisms. Here, we conducted an eQTL-guided function-related co-expression analysis to identify the putative regulators and construct gene regulatory network. We performed an eQTL analysis of 210 recombinant inbred lines (RILs) derived from a cross between two indica rice lines, Zhenshan 97 and Minghui 63, the parents of an elite hybrid, using data obtained by hybridizing RNA samples of flag leaves at the heading stage with Affymetrix whole-genome arrays. Making use of an ultrahigh-density single-nucleotide polymorphism bin map constructed by population sequencing, 13 647 eQTLs for 10 725 e-traits were detected, comprising 5079 cis-eQTLs (37.2%) and 8568 trans-eQTLs (62.8%). The analysis revealed 138 trans-eQTLs hotspots, each of which apparently regulates the expression variations of many genes. Co-expression analysis of functionally related genes within the framework of regulator-target relationships outlined by the eQTLs led to the identification of putative regulators in the system. The usefulness of the strategy was demonstrated with the genes known to be involved in flowering. We also applied this strategy to the analysis of QTLs for yield traits, which also suggested likely candidate genes. eQTL-guided co-expression analysis may provide a promising solution for outlining a framework for the complex regulatory network of an organism.

XIAO is involved in the control of organ size by contributing to the regulation of signaling and homeostasis of brassinosteroids and cell cycling in rice.

Organ size is determined by cell number and size, and involves two fundamental processes: cell proliferation and cell expansion. Although several plant hormones are known to play critical roles in shaping organ size by regulating the cell cycle, it is not known whether brassinosteroids (BRs) are also involved in regulating cell division. Here we identified a rice T-DNA insertion mutant for organ size, referred to as xiao, that displays dwarfism and erect leaves, typical BR-related phenotypes, together with reduced seed setting. XIAO is predicted to encode an LRR kinase. The small stature of the xiao mutant resulted from reduced organ sizes due to decreased cell numbers resulting from reduced cell division rate, as supported by the observed co-expression of XIAO with a number of genes involved in cell cycling. The xiao mutant displayed a tissue-specific enhanced BR response and greatly reduced BR contents at the whole-plant level. These results indicated that XIAO is a regulator of BR signaling and cell division. Thus, XIAO may provide a possible connection between BRs and cell-cycle regulation in controlling organ growth.

OsAP65, a rice aspartic protease, is essential for male fertility and plays a role in pollen germination and pollen tube growth.

Aspartic proteases (APs) comprise a large proteolytic enzyme family widely distributed in animals, microbes, viruses, and plants. The rice genome encodes 96 APs, of which only a few have been functionally characterized. Here, the identification and characterization of a novel AP gene, OsAP65, which plays an indispensable role in pollen tube growth in rice, is reported. The T-DNA insertion line of OsAP65 caused severe segregation distortion. In the progeny derived from an individual heterozygous for the T-DNA insertion, the wild type and T-DNA-carrying heterozygote segregated at a ratio close to 1:1, while homozygotes of disrupted OsAP65 (OsAP65-/-) were not recovered. Reciprocal crosses between heterozygotes and wild-type plants demonstrated that the mutant alleles could not be transmitted through the male gamete. Examination of the anthers from heterozygous plants revealed that the mutant pollen matured normally, but did not germinate or elongate. OsAP65 was expressed in various tissues and the transcript level in heterozygous plants was about half of the amount measured in the wild-type plants. The subcellular localization showed that OsAP65 is a pre-vacuolar compartment (PVC) protein. These results indicated that OsAP65 was essential for rice pollen germination and tube growth.

The bacterial pathogen Xanthomonas oryzae overcomes rice defenses by regulating host copper redistribution.

Pathogen effectors are virulence factors causing plant diseases. How the host targets of these effectors facilitate pathogen infection is largely unknown. An effector of Xanthomonas oryzae pv oryzae (Xoo) transcriptionally activates rice (Oryza sativa) susceptibility gene Xa13 to cause bacterial blight disease. Xa13 encodes an indispensable plasma membrane protein of the MtN3/saliva family, which is prevalent in eukaryotes with unknown biochemical function. We show that the XA13 protein cooperates with two other proteins, COPT1 and COPT5, to promote removal of copper from xylem vessels, where Xoo multiplies and spreads to cause disease. Copper, an essential micronutrient of plants and an important element for a number of pesticides in agriculture, suppresses Xoo growth. Xoo strain PXO99 is more sensitive to copper than other strains; its infection of rice is associated with activation of XA13, COPT1, and COPT5, which modulate copper redistribution in rice. The involvement of XA13 in copper redistribution has led us to propose a mechanism of bacterial virulence.

Linking differential domain functions of the GS3 protein to natural variation of grain size in rice.

Grain yield in many cereal crops is largely determined by grain size. Here we report the genetic and molecular characterization of GS3, a major quantitative trait locus for grain size. It functions as a negative regulator of grain size and organ size. The wild-type isoform is composed of four putative domains: a plant-specific organ size regulation (OSR) domain in the N terminus, a transmembrane domain, a tumor necrosis factor receptor/nerve growth factor receptor (TNFR/NGFR) family cysteine-rich domain, and a von Willebrand factor type C (VWFC) in the C terminus. These domains function differentially in grain size regulation. The OSR domain is both necessary and sufficient for functioning as a negative regulator. The wild-type allele corresponds to medium grain. Loss of function of OSR results in long grain. The C-terminal TNFR/NGFR and VWFC domains show an inhibitory effect on the OSR function; loss-of-function mutations of these domains produced very short grain. This study linked the functional domains of the GS3 protein to natural variation of grain size in rice.

Opposite functions of a rice mitogen-activated protein kinase during the process of resistance against Xanthomonas oryzae.

The pathogen-induced plant defense signaling network consists of multiple components, although only some of them are characterized. Most of the known components function either as activators or repressors in host-pathogen interactions. Here we report that a mitogen-activated protein kinase, OsMPK6, functions both as an activator and a repressor in rice resistance against Xanthomonas oryzae pv. oryzae (Xoo), the causal organism of bacterial blight disease. Activation of OsMPK6 resulted in the formation of lesion mimics and local resistance to Xoo, accompanied by the accumulation of salicylic acid (SA) and jasmonic acid (JA), and the induced expression of SA- and JA-signaling genes. Nuclear localization of OsMPK6 was essential for local resistance, suggesting that modulating the expression of defense-responsive genes through transcription regulators may be the primary mechanism of OsMPK6-mediated local resistance. The knock-out of OsMPK6 resulted in enhanced Xoo resistance, increased accumulation of SA and enhanced resistance to X. oryzae pv. oryzicola, the causal organism of bacterial streak disease, in systemic tissues. Xoo infection induced the expression of PR1a, the marker gene of systemic acquired resistance (SAR), in systemic health tissues of OsMPK6-knock-out plants. These results suggest that OsMPK6 negatively regulates SAR. Thus OsMPK6 is a two-faced player in the rice-Xoo interaction.

A global analysis of QTLs for expression variations in rice shoots at the early seedling stage.

Analyses of quantitative trait loci (QTLs) for expression levels (eQTLs) of genes reveal a genetic relationship between expression variation and the regulator, thus unlocking information for identifying the regulatory network. Oligo-nucleotide expression microarrays hybridized with RNA can simultaneously provide data for molecular markers and transcript abundance. In this study, we used an Affymetrix GeneChip Rice Genome Array to analyze eQTLs in rice shoots at 72 h after germination from 110 recombinant inbred lines (RILs) derived from a cross between Zhenshan 97 and Minghui 63. In total, 1632 single-feature polymorphisms (SFPs) plus 23 PCR markers were identified and placed into 601 recombinant bins, spanning 1459 cM in length, which were used as markers to genotype the RILs. We obtained 16,372 expression traits (e-traits) each with at least one eQTL, resulting in 26,051 eQTLs in total, including both cis- and trans-eQTLs. We also identified 171 eQTL hot spots in the rice genome, each of which controls transcript variations of many e-traits. Gene ontology analysis revealed an enrichment of certain functional categories of genes in some of the eQTL hot spots. In particular, eQTLs for e-traits involving the DNA metabolic process was significantly enriched in several eQTL hot spots on chromosomes 3, 5 and 10. Several e-traits co-localizing with cis-eQTLs showed significant correlations with hundreds of e-traits, indicating possible co-regulation. We also detected correlations between QTLs for shoot dry weight and eQTLs, revealing possible candidate genes for the trait. These results provided clues for the identification and characterization of the regulatory network in the whole genome at the transcriptional level.

A dynamic gene expression atlas covering the entire life cycle of rice.

Growth and development of a plant are controlled by programmed expression of suits of genes at the appropriate time, tissue and abundance. Although genomic resources have been developed rapidly in recent years in rice, a model plant for cereal genome research, data of gene expression profiling are still insufficient to relate the developmental processes to transcriptomes, leaving a large gap between the genome sequence and phenotype. In this study, we generated genome-wide expression data by hybridizing 190 Affymetrix GeneChip Rice Genome Arrays with RNA from 39 tissues collected throughout the life cycle of the rice plant from two varieties, Zhenshan 97 and Minghui 63. Analyses of the global transcriptomes revealed many interesting features of dynamic patterns of gene expression across the tissues and stages. In total, 38 793 probe sets were detected as expressed and 69% of the expressed transcripts showed significantly variable expression levels among tissues/organs. We found that similarity of transcriptomes among organs corresponded well to their developmental relatedness. About 5.2% of the expressed transcripts showed tissue-specific expression in one or both varieties and 22.7% of the transcripts exhibited constitutive expression including 19 genes with high and stable expression in all the tissues. This dataset provided a versatile resource for plant genomic research, which can be used for associating the transcriptomes to the developmental processes, understanding the regulatory network of these processes, tracing the expression profile of individual genes and identifying reference genes for quantitative expression analyses.

OsEDR1 negatively regulates rice bacterial resistance via activation of ethylene biosynthesis.

Rice OsEDR1 is a sequence ortholog of Arabidopsis EDR1. However, its molecular function is unknown. We show here that OsEDR1-suppressing/knockout (KO) plants, which developed spontaneous lesions on the leaves, have enhanced resistance to Xanthomonas oryzae pv. oryzae (Xoo) causing bacterial blight disease. This resistance was associated with increased accumulation of salicylic acid (SA) and jasmonic acid (JA), induced expression of SA- and JA-related genes and suppressed accumulation of 1-aminocyclopropane-1-carboxylic acid (ACC), the direct precursor of ethylene, and expression of ethylene-related genes. OsEDR1-KO plants also showed suppressed production of ethylene. Knockout of OsEDR1 suppressed the ACC synthase (ACS) gene family, which encodes the rate-limiting enzymes of ethylene biosynthesis by catalysing the formation of ACC. The lesion phenotype and enhanced bacterial resistance of the OsEDR1-KO plants was partly complemented by the treatment with ACC. ACC treatment was associated with decreased SA and JA biosynthesis in OsEDR1-KO plants. In contrast, aminoethoxyvinylglycine, the inhibitor of ethylene biosynthesis, promoted expression of SA and JA synthesis-related genes in OsEDR1-KO plants. These results suggest that ethylene is a negative signalling molecule in rice bacterial resistance. In the rice-Xoo interaction, OsEDR1 transcriptionally promotes the synthesis of ethylene that, in turn, suppresses SA- and JA-associated defence signalling.

A triallelic system of S5 is a major regulator of the reproductive barrier and compatibility of indica-japonica hybrids in rice.

Hybrid sterility is a major form of postzygotic reproductive isolation. Although reproductive isolation has been a key issue in evolutionary biology for many decades in a wide range of organisms, only very recently a few genes for reproductive isolation were identified. The Asian cultivated rice (Oryza sativa L.) is divided into two subspecies, indica and japonica. Hybrids between indica and japonica varieties are usually highly sterile. A special group of rice germplasm, referred to as wide-compatibility varieties, is able to produce highly fertile hybrids when crossed to both indica and japonica. In this study, we cloned S5, a major locus for indica-japonica hybrid sterility and wide compatibility, using a map-based cloning approach. We show that S5 encodes an aspartic protease conditioning embryo-sac fertility. The indica (S5-i) and japonica (S5-j) alleles differ by two nucleotides. The wide compatibility gene (S5-n) has a large deletion in the N terminus of the predicted S5 protein, causing subcellular mislocalization of the protein, and thus is presumably nonfunctional. This triallelic system has a profound implication in the evolution and artificial breeding of cultivated rice. Genetic differentiation between indica and japonica would have been enforced because of the reproductive barrier caused by S5-i and S5-j, and species coherence would have been maintained by gene flow enabled by the wide compatibility gene.

Pathogen-induced expressional loss of function is the key factor in race-specific bacterial resistance conferred by a recessive R gene xa13 in rice.

The fully recessive disease resistance (R) gene xa13, which mediates race-specific resistance to Xanthomonas oryzae pv. oryzae (Xoo), encodes a plasma membrane protein that differs by one amino acid from that encoded by its dominant (susceptible) allele Xa13. The molecular mechanism of xa13-mediated resistance is largely unknown. Here we show that, compared with its dominant allele, expressional non-reaction of xa13 to Xoo infection, not its protein composition, is the key factor for xa13-mediated resistance. We used the promoter (P(Xa13)) of the dominant Xa13, which was induced by only the incompatible Xoo strain for xa13, to regulate xa13 and xa13(Leu49) (a natural recessive allele of xa13) in the rice line IRBB13 carrying xa13. The transgenic plants showed the same level of susceptibility and bacterial growth rate as those of the rice line carrying dominant Xa13, accompanied by the induced accumulation of xa13 or xa13(Leu49) proteins. Constitutive expression of dominant XA13 or different xa13 proteins (xa13, xa13(Leu49), xa13(Ala85) or xa13(Val184)) in IRBB13 had no effect on Xoo infection in the transgenic plants. These results suggest that race-specific pathogen-induced Xa13 expression is critical for infection. Thus, xa13 stands out from other R genes in that its functions in disease resistance are due to only the loss of pathogen-induced transcriptional motivation caused by natural selection.

Four rice QTL controlling number of spikelets per panicle expressed the characteristics of single Mendelian gene in near isogenic backgrounds.

Development of quantitative trait loci (QTL) near isogenic lines is a crucial step to QTL isolation using the strategy of map-based cloning. In this study, a recombinant inbred line (RIL) population derived from two indica rice varieties, Zhenshan 97 and HR5, was employed to map QTL for spikelets per panicle (SPP). One major QTL (qSPP7) and three minor QTL (qSPP1, qSPP2 and qSPP3) were identified on chromosomes 7, 1, 2 and 3, respectively. Four sets of near isogenic lines (NILs) BC(4)F(2) targeted for the four QTL were developed by following a standard procedure of consecutive backcross, respectively. These QTL were not only validated in corresponding NILs, but also explained amounts of phenotypic variation with much larger LOD scores compared with those identified in RILs. SPP in the four QTL-NILs expressed bimodal or discontinuous distributions and followed the expected segregation ratio of single Mendelian factor by progeny test. Finally, qSPP1, qSPP2, qSPP3 and qSPP7 were respectively mapped to a locus, 0.5 cM from MRG2746, 0.6 cM from MRG2762, 0.8 cM from RM49 and 0.7 cM from MRG4436, as co-dominant markers on the basis of progeny tests. These results indicate no matter how small effect minor QTL is, QTL may still express the characteristics of single Mendelian factor in NILs and isolation of minor QTL will be possible using high quality NILs. Pyramiding these QTL into a variety will largely enhance rice grain yield.