Genetic dissection and validation of candidate genes for flag leaf size in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Key message

Two major loci with functional candidate genes were identified and validated affecting flag leaf size, which offer desirable genes to improve leaf architecture and photosynthetic capacity in rice.

Abstract

Leaf size is a major determinant of plant architecture and yield potential in crops. However, the genetic and molecular mechanisms regulating leaf size remain largely elusive. In this study, quantitative trait loci (QTLs) for flag leaf length and flag leaf width in rice were detected with high-density single nucleotide polymorphism genotyping of a chromosomal segment substitution line (CSSL) population, in which each line carries one or a few chromosomal segments from the japonica cultivar Nipponbare in a common background of the indica variety Zhenshan 97. In total, 14 QTLs for flag leaf length and nine QTLs for flag leaf width were identified in the CSSL population. Among them, qFW4-2 for flag leaf width was mapped to a 37-kb interval, with the most likely candidate gene being the previously characterized NAL1. Another major QTL for both flag leaf width and length was delimited by substitution mapping to a small region of 13.5 kb that contains a single gene, Ghd7.1. Mutants of Ghd7.1 generated using CRISPR/CAS9 approach showed reduced leaf size. Allelic variation analyses also validated Ghd7.1 as a functional candidate gene for leaf size, photosynthetic capacity and other yield-related traits. These results provide useful genetic information for the improvement of leaf size and yield in rice breeding programs.

     

Population structure analysis and association mapping of bacterial blight resistance in <Emphasis Type="Italic">indica</Emphasis> rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) accessions

Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is a serious disease in rice production worldwide. To understand the genetic diversity of bacterial blight resistance a population consisting of 175 indica accessions from nine countries was collected and detected their association between SSR (Simple Sequence Repeat) markers and resistance to six bacterial races. The resistance phenotypes of various rice accessions were evaluated through artificial inoculation under controlled conditions in 2013 and 2014. Association analysis showed that 17 SSR markers were significantly associated with resistance to four bacterial races and the phenotypic variations explained (PVE) ranged from 7.43 to 15.05%. Among the 17 associated SSR markers, two SSR markers located in previously reported genes regions, and 15 SSR markers were newly identified in this study. These results validated a new approach to map resistance genes of rice to bacterial blight. These markers could be used for marker-assisted selection (MAS) in rice bacterial blight resistance breeding programs.     

Genetic mapping and salt tolerance of a novel <Emphasis Type="Italic">D1</Emphasis>-allelic mutant of rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Plant height and grain shape are important traits that may affect yield in rice, and they therefore have enormous importance in breeding. A dwarf small-grain mutant (S525) was identified among progeny of the Indica rice restorer line ‘Xida 1B’ (wild type) raised from seeds treated with ethyl methanesulfonate. The dwarf and small-grain phenotypes were stably inherited after multi-generation selfing. Field-grown mutant plants showed the phenotypes of dwarfism, broad leaves, and small round grains. Genetic mapping and sequencing confirmed that S525 was a novel d1-allelic mutant. A single-base transition (G to A) in the functional dwarfism gene D1 at the conjunction site of the 11th intron caused excision or duplication of the 11th exon in the mRNA and resulted in translation of a defective G_α protein. The S525 showed enhanced salt tolerance compared with the wild type (WT), and the expression of genes associated with salt tolerance quantitatively increased in response to treatment with 200 mM NaCl. The S525 may be useful for future investigation of G_α functions and in the breeding of new dwarf rice cultivars.     

Fine genetic mapping and physical delimitation of the lesion mimic gene <Emphasis Type="Italic">spotted leaf 5</Emphasis> (<Emphasis Type="Italic">spl5</Emphasis>) in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Spotted leaf 5 (spl5), a lesion mimic mutant, was first identified in rice (Oryza sativa L.) japonica cv. Norin8 in 1978. This mutant exhibits spontaneous disease-like lesions in the absence of any pathogens and resistance to rice blast and bacterial blight; however, the target gene has not yet been isolated. In the present study, we employed a map-based cloning strategy to finely map the spl5 gene. In an initial mapping with 100 F₂ individuals (spl5/spl5) derived from a cross between the spl5 mutant and indica cv. 93-11, the spl5 gene was located in a 3.3-cM region on chromosome 7 using six simple sequence repeat (SSR) markers. In a high-resolution genetic mapping, two F₂ populations with 3,149 individuals (spl5/spl5) were derived from two crosses between spl5 mutant and two indica cvs. 93-11 and Zhefu802 and six sequence-tagged site (STS) markers were newly developed. Finally, the spl5 gene was mapped to a region of 0.048 cM between two markers SSR7 and RM7121. One BAC/PAC contig map covering these markers’ loci and the spl5 gene was constructed through Pairwise BLAST analysis. Our bioinformatics analysis shows that the spl5 gene is located in the 80-kb region between two markers SSR7 and RM7121 with a high average ratio of physical to genetic distance (1.67 Mb/cM) and eighteen candidate genes. The analysis of these candidate genes indicates that the spl5 gene represents a novel class of regulators controlling cell death and resistance response in plants.     

Characterization and fine mapping of <Emphasis Type="Italic">osh15</Emphasis>(<Emphasis Type="Italic">t</Emphasis>), a novel dwarf mutant gene in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Plant height is one of the most important agronomic traits of plant architecture, and also affects grain yield in rice. In this study, we obtained a novel dwarf rice mutant of japonica variety Shennong9816, designated Shennong9816d. Compared with wild-type, the Shennong9816d plant height was significantly reduced, and the tiller number significantly increased. Additionally, the mutant yield component, and the number of large and small vascular bundles were significantly decreased compared with wild-type. Genetic analysis indicated that the Shennong9816d dwarf phenotype was controlled by a recessive nuclear gene, while the plant was shown to be sensitive to gibberellic acid. Using a large F₂ population derived from a cross between Shennong9816d and the indica rice variety Habataki, the osh15(t) gene was fine mapped between RM20891 and RM20898, within a physical distance of 73.78 kb. Sequencing analysis showed that Shennong9816d carries a 1 bp mutation and a 30 bp insertion in the OSH15 region. These results suggest that osh15(t) is a novel allelic mutant originally derived from japonica variety Shennong9816, which may be useful for introducing the semi-dwarf phenotype to improve plant architecture in rice breeding practice.     

Transposition behavior of nonautonomous a <Emphasis Type="Italic">hAT</Emphasis> superfamily transposon <Emphasis Type="Italic">nDart</Emphasis> in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Transposable elements (TEs) have a significant impact on the evolution of gene function and genome structures. An endogenous nonautonomous transposable element nDart was discovered in an albino mutant that had an insertion in the Mg-protoporphyrin IX methyltransferase gene in rice. In this study, we elucidated the transposition behavior of nDart, the frequency of nDart transposition and characterized the footprint of nDart. Novel independent nDart insertions in backcrossed progenies were detected by DNA blotting analysis. In addition, germinal excision of nDart occurred at very low frequency compared with that of somatic excision, 0–13.3%, in the nDart1-4(3-2) and nDart1-A loci by a locus-specific PCR strategy. A total of 253 clones from somatic excision at five nDart loci in 10 varieties were determined. nDart rarely caused deletions beyond target site duplication (TSD). The footprint of nDart contained few transversions of nucleotides flanking to both sides of the TSD. The predominant footprint of nDart was an 8-bp addition. Precise excision of nDart was detected at a rate of only 2.2%, which occurred at two loci among the five loci examined. Furthermore, the results in this study revealed that a highly conserved mechanism of transposition is involved between maize Ac/Ds and rice Dart/nDart, which are two-component transposon systems of the hAT superfamily transposons in plant species.     

Genetic differentiation for nuclear,mitochondrial and chloroplast genomes in common wild rice (<Emphasis Type="Italic">Oryza rufipogon</Emphasis> Griff.) and cultivated rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

The genetic differentiation of nuclear, mitochondrial (mt) and chloroplast (cp) genomes was investigated by Southern and PCR analysis using 75 varieties of cultivated rice ( Oryza sativa L.) and 118 strains of common wild rice (CWR, Oryza rufipogon Griff.) from ten countries of Asia. The distinguishing differences between the Indica and Japonica cultivars were detected both in the nuclear genome and the cytoplasmic genome, confirming that the Indica-Japonica differentiation is of major importance for the three different classes of genome in cultivated rice. This differentiation was also detected in common wild rice with some differences among the genome compartments and the various regions. For nuclear DNA variation, both Indica-like and Japonica-like types were observed in the Chinese CWR, with the latter more-frequent than the former. No Japonica-like type was found in South Asia, and only two strains of the Japonica-like type were detected in Southeast Asia, thus the Indica-like type is the major type among South and Southeast Asian CWR. For mtDNA, only a few strains of the Japonica-like type were detected in CWR. For cpDNA, the Japonica type was predominant among the CWR strains from China, Bangladesh and Burma, while the Indica type was predominant among the CWR strains from Thailand, Malaysia, Cambodia and Sri Lanka, and both types were found in similar frequencies among the Indian CWR. Altogether, however, the degree of Indica-Japonica differentiation in common wild rice was much-less important than that in cultivated rice. Cluster analyses for nuclear and mitochondrial DNA variation revealed that some CWR strains showed large genetic distances from cultivated rice and formed clusters distinct from cultivated rice. Coincidence in the genetic differentiation between the three different classes of genome was much higher in cultivated rice than in CWR. Among the 75 cultivars, about 3/4 entries were "homoeotype" showing congruent results for nuclear, mt and cpDNA regarding the Indica-Japonica differentiation. In CWR, the proportions of homoeotypes were 5.7%, 15% and 48.8% in China, South Asia and Southeast Asia, respectively. Based on the average genetic distance among all the strains of CWR and cultivated rice for nuclear and mitochondrial genomes, the variability of the nuclear genome was found to be higher than that of the mitochondrial genome. The global pattern based on all genomes shows much-more diversification in CWR than that in cultivated rice.     

Functional markers for <Emphasis Type="Italic">xa5</Emphasis>-mediated resistance in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis>, L.)

The recent cloning of several agronomically important genes has facilitated the development of functional markers. These markers reside within the target genes themselves and can be used with great reliability and efficiency to identify favorable alleles in a breeding program. Bacterial blight (BB) is a severe rice disease throughout the world that is controlled primarily through use of resistant cultivars. xa5 is a race-specific, recessive gene mediating resistance to BB. It is widely used in rice breeding programs throughout the tropics. Due to its recessive nature, phenotypic selection for xa5-mediated resistance is both slow and costly. Previously, marker assisted selection (MAS) for this resistance gene was not efficient because it involved markers that were only indirectly linked to xa5 and ran the risk of being separated from the trait by recombination. Recently, the cloning of the gene underlying this trait made it possible to develop functional markers. Here we present a set of CAPS markers for easy, quick and direct identification of cultivars or progeny carrying xa5-mediated resistance and provide evidence that these markers are 100% predictive of the presence of the xa5 allele. These markers are expected to enhance the reliability and cost-effectiveness of MAS for xa5-mediated resistance.     

Fine mapping of the <Emphasis Type="Italic">qCTS4</Emphasis> locus associated with seedling cold tolerance in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Rice seedlings are sensitive to low temperatures (≤15°C) and under prolonged or repeated exposure, yellowing and stunting are commonly observed. Damage to seedlings results in poor stand establishment and delayed maturation, which can cause significant reductions in yield. In general, japonica rice varieties exhibit more cold tolerance than indica varieties. Earlier genetic analysis of the California rice variety M202 revealed several quantitative trait loci (QTL) that contribute to its tolerance to low temperatures in comparison to the indica rice variety IR50. Among these QTL, qCTS4 is associated with tolerance to yellowing and stunting of rice seedlings and accounts for 40% of the phenotypic variation. Here we report on the fine mapping of qCTS4 to a 128 kb region of chromosome 4, which is highly suppressed for recombination in our mapping populations. Our results provide the necessary foundation for identifying the gene(s) underlying qCTS4 and the markers developed here may be used to introgress this region into indica varieties to improve seedling tolerance to low temperatures. The mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.     

High-resolution genetic mapping at the <Emphasis Type="Italic">Bph15</Emphasis> locus for brown planthopper resistance in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Resistance to the brown planthopper (BPH), Nilaparvata lugens Stål, a devastating sucking insect pest of rice, is an important breeding objective in rice improvement programs. Bph15, one of the 17 major BPH resistance genes so far identified in both cultivated and wild rice, has been identified in an introgression line, B5, and mapped on chromosome 4 flanked by restriction fragment length polymorphism markers C820 and S11182. In order to pave the way for positional cloning of this gene, we have developed a high-resolution genetic map of Bph15 by positioning 21 DNA markers in the target chromosomal region. Mapping was based on a PCR-based screening of 9,472 F₂ individuals derived from a cross between RI93, a selected recombinant inbred line of B5 bearing the resistance gene Bph15, and a susceptible variety, Taichung Native 1, in order to identify recombinant plants within the Bph15 region. Recombinant F₂ individuals with the Bph15 genotype were determined by phenotype evaluation. Analysis of recombination events in the Bph15 region delimited the gene locus to an interval between markers RG1 and RG2 that co-segregated with the M1 marker. A genomic library of B5 was screened using these markers, and bacterial artificial chromosome clones spanning the Bph15 chromosome region were obtained. An assay of the recombinants using the sub-clones of these clones in combination with sequence analysis delimited the Bph15 gene to a genomic segment of approximately 47 kb. This result should serve as the basis for eventual isolation of the Bph15 resistance gene.     

Genetic analysis of starch paste viscosity parameters in glutinous rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Starch paste viscosity plays an important role in estimating the cooking, eating, and processing quality of rice. The inheritance of starch paste viscosity in glutinous rice remains undefined. In the present study, 118 glutinous rice accessions were collected, and the genotypes of 17 starch synthesis-related genes (SSRG) were analyzed by using 43 gene-specific molecular markers. Association analysis indicated that 10 of 17 SSRGs were involved in controlling the rapid visco analyzer (RVA) profile parameters. Among these, the PUL gene was identified to play an important role in control of peak viscosity (PKV), hot paste viscosity (HPV), cool paste viscosity (CPV), breakdown viscosity (BDV), peak time (PeT), and paste temperature (PaT) in glutinous rice. Other SSRGs involved only a few RVA profile parameters. Furthermore, interactions between SSRGs were found being responsible for PeT, PaT, and BDV. Some of the RVA parameters, including PKV, HPV, CPV, CSV, and PaT, were mainly governed by single SSRG, whereas other parameters, such as BDV, SBV, and PeT, were controlled by a few SSRGs, functioning cooperatively. Further, three near-isogenic lines (NIL) of a japonica glutinous cv. Suyunuo as genetic background, with PUL, SSIII-1, and SSIII-2 alleles replaced with those of indica cv. Guichao 2, were employed to verify the genetic effects of the various genes, and the results were consistent with those obtained from the association analysis. These findings indicated that starch paste viscosity in glutinous rice had a complex genetic system, and the PUL gene played an important role in determining the RVA profile parameters in glutinous rice. These results provide important information for potentially improving the quality of glutinous rice.     

Genetic diversity in basmati rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) germplasm as revealed by microsatellite (SSR) markers

Genetic diversity among rice genotypes, including 15 indica basmati advance lines and 5 basmati improved varieties were investigated by 28 SSR markers including one indel marker. The SSRs covered all the 12 chromosomes that distributed across the rice genomes. The mean number of alleles per locus was 3.60, showing average number of polymorphism information content was 0.48. A total of 101 alleles were also identified from the microsatellite marker loci. A number of SSR markers were also identified that could be utilized to differentiate between rice genotypes. Pair wise Nei’s genetic distance between rice genotypes ranged from 0.07 to 0.95. The dendrogram based on cluster analysis by using SSR polymorphism that grouped the 20 genotypes of rice in to five clusters based on their genetic similarity. The result could be useful for the identification and selection of the diverse genotypes for the future cross breeding program and development of new rice varieties.     

Fine mapping of <Emphasis Type="Italic">S31</Emphasis>, a gene responsible for hybrid embryo-sac abortion in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Partial abortion of female gametes and the resulting semi-sterility of indica × japonica inter-subspecific rice hybrids have been ascribed to an allelic interaction, which can be avoided by the use of wide compatibility varieties. To further understand the genetic mechanism of hybrid sterility, we have constructed two sets of hybrids, using as male parent either the typical japonica variety Asominori, or the wide compatibility variety 02428; and as female, a set of 66 chromosome segment substitution lines in which various chromosomal segments from the indica variety IR24 have been introduced into a common genetic background of Asominori. Spikelet semi-sterility was observed in hybrid between CSSL34 and Asominori, which is known to carry the sterility gene S31 (Zhao et al. in Euphytica 151:331–337, 2006). Cytological analysis revealed that the semi-sterility of the CSSL34 × Asominori hybrid was caused primarily by partial abortion of the embryo sac at the stage of the mitosis of the functional megaspore. A population of 1,630 progeny of the three-way cross (CSSL34 × 02428) × Asominori was developed to map S31. Based on the physical location of linked molecular markers, S31 was thereby delimited to a 54-kb region on rice chromsome 5. This fragment contains eight predicted open reading frames, four of which encode known proteins and four putative proteins. These results are relevant to the map-based cloning of S31, and the development of marker-assisted transfer of non-sterility allele inducing alleles to breeding germplasm, to allow for a more efficient exploitation of heterosis in hybrid rice.     

SNP haplotypes of the <Emphasis Type="Italic">BADH1</Emphasis> gene and their association with aroma in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Betaine aldehyde dehydrogenase (BADH) is a key enzyme involved in the synthesis of glycinebetaine—a powerful osmoprotectant against salt and drought stress in a large number of species. Rice is not known to accumulate glycinebetaine but it has two functional genes coding for the BADH enzyme. A non-functional allele of the BADH2 gene located on chromosome 8 is a major factor associated with rice aroma. However, similar information is not available regarding the BADH1 gene located on chromosome 4 despite the similar biochemical function of the two genes. Here we report on the discovery and validation of SNPs in the BADH1 gene by re-sequencing of diverse rice varieties differing in aroma and salt tolerance. There were 17 SNPs in introns with an average density of one per 171 bp, but only three SNPs in exons at a density of one per 505 bp. Each of the three exonic SNPs led to changes in amino acids with functional significance. Multiplex SNP assays were used for genotyping of 127 diverse rice varieties and landraces. In total 15 SNP haplotypes were identified but only four of these, corresponding to two protein haplotypes, were common, representing more than 85% of the cultivars. Determination of population structure using 54 random SNPs classified the varieties into two groups broadly corresponding to indica and japonica cultivar groups, aromatic varieties clustering with the japonica group. There was no association between salt tolerance and the common BADH1 haplotypes, but aromatic varieties showed specific association with a BADH1 protein haplotype (PH2) having lysine₁₄₄ to asparagine₁₄₄ and lysine₃₄₅ to glutamine₃₄₅ substitutions. Protein modeling and ligand docking studies show that these two substitutions lead to reduction in the substrate binding capacity of the BADH1 enzyme towards gamma-aminobutyraldehyde (GABald), which is a precursor of the major aroma compound 2-acetyl-1-pyrroline (2-AP). This association requires further validation in segregating populations for potential utilization in the rice breeding programs.     

Quantitative proteomics analysis of proteins involved in leaf senescence of rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Delaying leaf senescence and prolonging the available time for photosynthesis is one of the important approaches to increase grain yield of rice (Oryza sativa L.). Here, iTRAQ-based quantitative proteomics approach was used to comparative analyze the expression profiles of proteins in rice leaves in response to senescence. Totally 5067 proteins were identified. Compared with the proteins in the flag leaves at early stage of grain-filling in rice Liang-You-Pei 9 (LYP9), 240 and 188 proteins were up-regulated and down-regulated in the flag leaves at middle stage of grain-filling, and 387 and 202 proteins were up-regulated and down-regulated in the flag leaves at late stage of grain-filling, respectively. In addition, 39 and 18 identified proteins were constantly up-regulated and down-regulated in the leaves from early to middle and late stages of grain-filling, respectively. Among them, chloroplast chaperonin 10, geranylgeranyl diphosphate reductase, Mg chelatase subunit ChLD, porphobilinogen deaminase, protochlorophyllide reductase B and thioredoxin-like protein CITRX might have involved in the senescence of leaves. This study provided important information for understanding the age-sensitive mechanism of LYP9, and offered a foundation for future studying and improving it.     

Assessment of genetic diversity of <Emphasis Type="Italic">Saltol</Emphasis> QTL among the rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) genotypes

Eight Saltol quantitative trait locus (QTL) linked simple sequence repeat (SSR) markers of rice (Oryza sativa L.) were used to study the polymorphism of this QTL in 142 diverse rice genotypes that comprised salt tolerant as well as sensitive genotypes. The SSR profiles of the eight markers generated 99 alleles including 20rare alleles and 16 null alleles. RM8094 showed the highest number (13) of alleles followed by RM3412 (12), RM562 (11), RM493 (9) and RM1287 (8) while as, RM10764 and RM10745 showed the lowest number (6) of alleles. Based on the highest number of alleles and PIC value (0.991), we identified RM8094 as suitable marker for discerning salt tolerant genotypes from the sensitive ones. Based upon the haplotype analysis using FL478 as a reference (salt tolerant genotypes containing Saltol QTL), we short listed 68 rice genotypes that may have at least one allele of FL478 haplotype. Further study may confirm that some of these genotypes might have Saltol QTL and can be used as alternative donors in salt tolerant rice breeding programmes.     

Identification and fine mapping of <Emphasis Type="Italic">qWBPH11</Emphasis> conferring resistance to whitebacked planthopper (<Emphasis Type="Italic">Sogatella furcifera</Emphasis> Horvath) in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

The whitebacked planthopper (WBPH), Sogatella furcifera Horvath, is one of the most destructive pests in rice (Oryza sativa L.) production. Host-plant resistance has been considered as an efficient and eco-friendly strategy to reduce yield losses caused by WBPH. In this study, we found that an indica rice cultivar IR54751-2-44-15-24-2 (IR54751) displayed high resistance to WBPH at both seedling and tillering stages. The resistance of IR54751 was mainly contributed by antixenosis and tolerance rather than antibiosis. An F₂ population derived from a cross between IR54751 and a susceptible japonica cultivar 02428 was constructed to detect the quantitative trait loci (QTLs) conferring the resistance to WBPH. In total, four QTLs including qWBPH3.1, qWBPH3.2, qWBPH11, and qWBPH12 were identified and distributed on three different chromosomes. The four QTLs had LOD scores of 3.8, 8.2, 5.8, and 3.9, accounting for 8.2, 21.5, 13.9, and 10.4% of the phenotypic variation, respectively. Except for qWBPH3.1, the resistance alleles of the other three QTLs were all from IR54751. Further, a secondary population harboring only single qWBPH11 locus was developed from the F₂ population by marker-assisted selection. Finally, qWBPH11 was delimited in a 450-kb region between markers DJ53973 and SNP56. The identification of WBPH resistance QTLs and the fine mapping of qWBPH11 will be helpful for cloning resistance genes and breeding resistant rice cultivars.     

Map-based cloning and functional analysis of the chromogen gene <Emphasis Type="Italic">C</Emphasis> in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

The chromogen gene C is critical for anthocyanin regulation in rice, and apiculus color is an important agronomic trait in selective breeding and variety purification. Mapbased cloning and in-depth functional analysis of the C gene will be useful for understanding the molecular mechanism of anthocyanin biosynthesis and for rice breeding. Japonica landrace Lijiangxintuanheigu (LTH) has red apiculi and purple stigmas. Genetic analysis showed that red apiculus and purple stigma in LTH co-segregated indicating control by a single dominant gene, or by two completely linked genes. Using 1,851 recessive individuals from two F₂ populations, the target gene OsC was delimited to a 70.8 kb interval on chromosome 6 that contains the rice homologue of the maize anthocyanin regulatory gene C1. When the entire OsC gene and its full-length cDNA cloned from LTH were transformed into japonica cultivar Kitaake with colorless apiculi and stigmas all positive transformants had red apiculi but non-colored stigmas, validating that OsC alone was responsible for the apiculus color and represented the functional C gene. OsC was constitutively expressed in all tissues examined, with strongest expression in leaf blades. These results set a foundation to clarify the regulatory mechanisms of OsC in the anthocyanin biosynthetic pathway.