Development of 12 sets of chromosome segment substitution lines that enhance allele mining in Asian cultivated rice.

Many agronomic traits that are important in rice breeding are controlled by multiple genes. The extensive time and effort devoted so far to identifying and selecting such genes are still not enough to target multiple agronomic traits in practical breeding in Japan because of a lack of suitable plant materials in which to efficiently detect and validate beneficial alleles from diverse genetic resources. To facilitate the comprehensive analysis of genetic variation in agronomic traits among Asian cultivated rice, we developed 12 sets of chromosome segment substitution lines (CSSLs) with the japonica background, 11 of them in the same genetic background, using donors representing the genetic diversity of Asian cultivated rice. Using these materials, we overviewed the chromosomal locations of 1079 putative QTLs for seven agronomic traits and their allelic distribution in Asian cultivated rice through multiple linear regression analysis. The CSSLs will allow the effects of putative QTLs in the highly homogeneous japonica background to be validated.

Prediction of heading date, culm length, and biomass from canopy-height-related parameters derived from time-series UAV observations of rice.

Unmanned aerial vehicles (UAVs) are powerful tools for monitoring crops for high-throughput phenotyping. Time-series aerial photography of fields can record the whole process of crop growth. Canopy height (CH), which is vertical plant growth, has been used as an indicator for the evaluation of lodging tolerance and the prediction of biomass and yield. However, there have been few attempts to use UAV-derived time-series CH data for field testing of crop lines. Here we provide a novel framework for trait prediction using CH data in rice. We generated UAV-based digital surface models of crops to extract CH data of 30 Japanese rice cultivars in 2019, 2020, and 2021. CH-related parameters were calculated in a non-linear time-series model as an S-shaped plant growth curve. The maximum saturation CH value was the most important predictor for culm length. The time point at the maximum CH contributed to the prediction of days to heading, and was able to predict stem and leaf weight and aboveground weight, possibly reflecting the association of biomass with duration of vegetative growth. These results indicate that the CH-related parameters acquired by UAV can be useful as predictors of traits typically measured by hand.

Remote-Sensing-Combined Haplotype Analysis Using Multi-Parental Advanced Generation Inter-Cross Lines Reveals Phenology QTLs for Canopy Height in Rice.

High-throughput phenotyping systems with unmanned aerial vehicles (UAVs) enable observation of crop lines in the field. In this study, we show the ability of time-course monitoring of canopy height (CH) to identify quantitative trait loci (QTLs) and to characterise their pleiotropic effect on various traits. We generated a digital surface model from low-altitude UAV-captured colour digital images and investigated CH data of rice multi-parental advanced generation inter-cross (MAGIC) lines from tillering and heading to maturation. Genome-wide association studies (GWASs) using the CH data and haplotype information of the MAGIC lines revealed 11 QTLs for CH. Each QTL showed haplotype effects on different features of CH such as stage-specificity and constancy. Haplotype analysis revealed relationships at the QTL level between CH and, vegetation fraction and leaf colour [derived from UAV red-green-blue (RGB) data], and CH and yield-related traits. Noticeably, haplotypes with canopy lowering effects at qCH1-4, qCH2, and qCH10-2 increased the ratio of panicle weight to leaf and stem weight, suggesting biomass allocation to grain yield or others through growth regulation of CH. Allele mining using gene information with eight founders of the MAGIC lines revealed the possibility that qCH1-4 contains multiple alleles of semi-dwarf 1 (sd1), the IR-8 allele of which significantly contributed to the "green revolution" in rice. This use of remote-sensing-derived phenotyping data into genetics using the MAGIC lines gives insight into how rice plants grow, develop, and produce grains in phenology and provides information on effective haplotypes for breeding with ideal plant architecture and grain yield.

Genetic Elucidation for Response of Flowering Time to Ambient Temperatures in Asian Rice Cultivars.

Climate resilience of crops is critical for global food security. Understanding the genetic basis of plant responses to ambient environmental changes is key to developing resilient crops. To detect genetic factors that set flowering time according to seasonal temperature conditions, we evaluated differences of flowering time over years by using chromosome segment substitution lines (CSSLs) derived from japonica rice cultivars "Koshihikari" × "Khao Nam Jen", each with different robustness of flowering time to environmental fluctuations. The difference of flowering times in 9 years' field tests was large in "Khao Nam Jen" (36.7 days) but small in "Koshihikari" (9.9 days). Part of this difference was explained by two QTLs. A CSSL with a "Khao Nam Jen" segment on chromosome 11 showed 28.0 days' difference; this QTL would encode a novel flowering-time gene. Another CSSL with a segment from "Khao Nam Jen" in the region around Hd16 on chromosome 3 showed 23.4 days" difference. A near-isogenic line (NIL) for Hd16 showed 21.6 days' difference, suggesting Hd16 as a candidate for this QTL. RNA-seq analysis showed differential expression of several flowering-time genes between early and late flowering seasons. Low-temperature treatment at panicle initiation stage significantly delayed flowering in the CSSL and NIL compared with "Koshihikari". Our results unravel the molecular control of flowering time under ambient temperature fluctuations.

Genomic Regions Involved in Differences in Eating and Cooking Quality Other than Wx and Alk Genes between indica and japonica Rice Cultivars.

BACKGROUND: In temperate rice cultivation regions, japonica rice cultivars are grown preferentially because consumers deem them to have good eating quality, whereas indica rice cultivars have high grain yields and strong heat tolerance but are considered to have poor eating quality. To mitigate the effects of global warming on rice production, it is important to develop novel rice cultivars with both desirable eating quality and resilience to high temperatures. Eating quality and agronomic traits were evaluated in a reciprocal set of chromosome segment substitution lines derived from crosses between a japonica rice cultivar 'Koshihikari' and an indica rice cultivar 'Takanari'. RESULTS: We detected 112 QTLs for amylose and protein contents, whiteness, stickiness, hardness and eating quality of cooked rice grains. Almost of 'Koshihikari' chromosome segments consistently improved eating quality. Among detected QTLs, six QTLs on chromosomes 1-5 and 11 were detected that increased whiteness and stickiness of cooked grains or decreased their hardness for 3 years. The QTLs on chromosomes 2-4 were not associated with differences in amylose or protein contents. QTLs on chromosomes 1-5 did not coincide with QTLs for agronomic traits such as heading date, culm length, panicle length, spikelet fertility and grain yield. Genetic effects of the detected QTLs were confirmed in substitution lines carrying chromosome segments from five other indica cultivars in the 'Koshihikari' genetic background. CONCLUSION: The detected QTLs were associated with differences in eating quality between indica and japonica rice cultivars. These QTLs appear to be widely distributed among indica cultivars and to be novel genetic factors for eating quality traits because their chromosome regions differed from those of the GBSSI (Wx) and SSIIa (Alk) genes. The detected QTLs would be very useful for improvement of eating quality of indica rice cultivars in breeding programs.

Haplotype analysis from unmanned aerial vehicle imagery of rice MAGIC population for the trait dissection of biomass and plant architecture.

Unmanned aerial vehicles (UAVs) are popular tools for high-throughput phenotyping of crops in the field. However, their use for evaluation of individual lines is limited in crop breeding because research on what the UAV image data represent is still developing. Here, we investigated the connection between shoot biomass of rice plants and the vegetation fraction (VF) estimated from high-resolution orthomosaic images taken by a UAV 10 m above a field during the vegetative stage. Haplotype-based genome-wide association studies of multi-parental advanced generation inter-cross (MAGIC) lines revealed four quantitative trait loci (QTLs) for VF. VF was correlated with shoot biomass, but the haplotype effect on VF was better correlated with that on shoot biomass at these QTLs. Further genetic characterization revealed the relationships between these QTLs and plant spreading habit, final shoot biomass and panicle weight. Thus, genetic analysis using high-throughput phenotyping data derived from low-altitude, high-resolution UAV images during early stages of rice growing in the field provides insights into plant growth, architecture, final biomass, and yield.

Surveillance of panicle positions by unmanned aerial vehicle to reveal morphological features of rice.

Rice plant architecture affects biomass and grain yield. Thus, it is important to select rice genotypes with ideal plant architecture. High-throughput phenotyping by use of an unmanned aerial vehicle (UAV) allows all lines in a field to be observed in less time than with traditional procedures. However, discrimination of plants in dense plantings is difficult, especially during the reproductive stage, because leaves and panicles overlap. Here, we developed an original method that relies on using UAV to identify panicle positions for dissecting plant architecture and to distinguish rice lines by detecting red flags attached to panicle bases. The plant architecture of recombinant inbred lines derived from Japanese cultivars 'Hokuriku 193' and 'Mizuhochikara', which differ in plant architecture, was assessed using a commercial camera-UAV system. Orthomosaics were made from UAV digital images. The center of plants was plotted on the image during the vegetative stage. The horizontal distance from the center to the red flag during the reproductive stage was used as the panicle position (PP). The red flags enabled us to recognize the positions of the panicles at a rate of 92%. The PP phenotype was related to but was not identical with the phenotypes of the panicle base angle, leaf sheath angle, and score of spreading habit. These results indicate that PP on orthomosaics could be used as an index of plant architecture under field conditions.

Detection of heading date QTLs in advanced-backcross populations of an elite indica rice cultivar, IR64.

IR64 is one of the world's most popular rice cultivars. To collect genetic factors involved in controlling its heading date, we developed 70 reciprocal advanced-backcross populations with a total of 6284 individuals at the BC4F2 generation from crosses between Koshihikari and IR64. We detected 29 QTLs associated with heading date on chromosomes 3, 5-8, 10, and 12. Twenty QTLs were located in the same chromosome regions as previously isolated heading date genes (Hd1, Hd6, Hd16, Ghd7, DTH8, Hd17, and Hd18). The rest were located in other chromosome regions. We found more number of QTLs than previous studies using mapping populations of IR64. Fine mapping in additional advanced-backcross populations clearly revealed that QTLs on the long arm of chromosome 7 are overlapping and seem to be a novel genetic factor for heading date because of their different locations from OsPRR37. Our results suggest that the difference in heading date between IR64 and Koshihikari is genetically controlled by many factors, and that a non-functional allele of Hd1 contributes to early heading of IR64 in the genetic background of functional alleles of other heading date QTLs and genes such as Hd6 and Hd16.

Switching genetic effects of the flowering time gene Hd1 in LD conditions by Ghd7 and OsPRR37 in rice.

Flowering time control in plants is a major limiting factor on the range of species. Day length, perceived via the photoperiodic pathway, is a critical factor for the induction of flowering. The module of GIGANTEA (GI)-CONSTANS (CO)-FLOWERING LOCUS T in the long day (LD) plant Arabidopsis is conserved in diverse plant species including the short day (SD) plant rice, where this module comprises OsGI-Heading date 1 (Hd1)-Heading date 3a. Hd1, the rice ortholog of Arabidopsis CO, has dual functions in the regulation of flowering time, promoting flowering in SD conditions and delaying it in LD conditions. We herein show genetic interactions among three LD repressor genes: Hd1, Grain number, plant height and heading date 7 (Ghd7), and Oryza sativa Pseudo-Response Regulator37 (OsPRR37). Genetic analyses, including segregation analyses, evaluations of near isogenic lines, and transformation for flowering time demonstrated that Hd1 promoted flowering time in inductive SD and non-inductive LD conditions in genetic condition of loss-of-function Ghd7 and OsPRR37 (ghd7osprr37) in rice. Functional Ghd7 or OsPRR37 may switch the genetic effects of Hd1 from the promotion to the delay of flowering times in LD conditions.

Discovery of QTL Alleles for Grain Shape in the Japan-MAGIC Rice Population Using Haplotype Information.

A majority of traits are determined by multiple quantitative trait loci (QTL) that can have pleiotropic effects. A multi-parent advanced generation inter-cross (MAGIC) population is well suited for genetically analyzing the effects of multiple QTL on traits of interest because it contains a higher number of QTL alleles than a biparental population. We previously produced the JAPAN-MAGIC (JAM) population, derived from eight rice (Oryza sativa L.) cultivars with high yield and biomass in Japan, and developed the method of genome-wide association study (GWAS) using haplotype information on the JAM lines. This method was effective for identifying major genes such as Waxy for eating quality and Sd1 for culm length. Here, we show that haplotype-based GWAS is also effective for the evaluation of multiple QTL with small effects on rice grain shape in the JAM lines. Although both the haplotype- and SNP-based GWAS identified multiple QTL for grain length and width, the sum of the estimated trait values of each allele for the QTL detected by haplotype-based GWAS had higher correlation with observed values than those detected by SNP-based GWAS, indicating high-accuracy QTL detection in the haplotype-based GWAS. Furthermore, the study revealed pleiotropic effects of some QTL regions in regulation of grain shape, suggesting that the haplotype-based GWAS using the JAM lines is an effective means to evaluate the main and side effects of haplotypes at each QTL. Information on the pleiotropic effects of haplotypes on various traits will be useful for designing ideal lines in a breeding program.

Genetic dissection of pre-harvest sprouting resistance in an upland rice cultivar.

Seed dormancy is important in rice breeding because it confers resistance to pre-harvest sprouting (PHS). To detect quantitative trait loci (QTLs) for pre-harvest sprouting resistance, we used chromosome segment substitution lines (CSSLs) derived from a cross between the Japanese upland rice cultivar 'Owarihatamochi' and the lowland rice cultivar 'Koshihikari'. In the CSSLs, several chromosomal regions were associated with PHS resistance. Among these, the chromosome 9 segment from 'Owarihatamochi' had the greatest association with increased PHS resistance. Further QTL analysis using an advanced backcross population (BC4F2) derived from a 'Koshihikari' × 'Owarihatamochi' cross revealed two putative QTLs, here designated qSDR9.1 (Seed dormancy 9.1) and qSDR9.2, on chromosome 9. The 'Owarihatamochi' alleles of the two QTLs reduced germination. Further fine mapping revealed that qSDR9.1 and qSDR9.2 were located within 4.1-Mb and 2.3-Mb intervals (based on the 'Nipponbare' reference genome sequence) defined by the simple sequence repeat marker loci RM24039 and RM24260 and Indel_2 and RM24540, respectively. We thus identified two QTLs for PHS resistance in 'Owarihatamochi', even though resistance levels are relatively low in this cultivar. This unexpected finding suggests the advantages of using CSSLs for QTL detection.

Haplotype-based allele mining in the Japan-MAGIC rice population.

Multi-parent advanced generation inter-cross (MAGIC) lines have broader genetic variation than bi-parental recombinant inbred lines. Genome-wide association study (GWAS) using high number of DNA polymorphisms such as single-nucleotide polymorphisms (SNPs) is a popular tool for allele mining in MAGIC populations, in which the associations of phenotypes with SNPs are investigated; however, the effects of haplotypes from multiple founders on phenotypes are not considered. Here, we describe an improved method of allele mining using the newly developed Japan-MAGIC (JAM) population, which is derived from eight high-yielding rice cultivars in Japan. To obtain information on the haplotypes in the JAM lines, we predicted the haplotype blocks in the whole chromosomes using 16,345 SNPs identified via genotyping-by-sequencing analysis. Using haplotype-based GWAS, we clearly detected the loci controlling the glutinous endosperm and culm length traits. Information on the alleles of the eight founders, which was based on the effects of mutations revealed by the analysis of next-generation sequencing data, was used to narrow down the candidate genes and reveal the associations between alleles and phenotypes. The haplotype-based allele mining (HAM) proposed in this study is a promising approach to the detection of allelic variation in genes controlling agronomic traits in MAGIC populations.

Detection of novel QTLs qDTH4.5 and qDTH6.3, which confer late heading under short-day conditions, by SSR marker-based and QTL-seq analysis.

Heading date is one of the most important traits in rice breeding. It is governed by multiple genes, including known quantitative trait loci (QTLs). In general, almost all japonica cultivars, including Nipponbare, head early under short-day (SD) conditions, but some indica cultivars, including Kasalath, head late. To explain this difference, we identified QTLs controlling heading date under SD conditions. We used NILs, CSSLs, and BILs from a cross between Nipponbare and Kasalath, and evaluated days to heading (DTH) under SD conditions. No NILs or CSSLs showed late heading, but two BILs (BIL-55 and BIL-78) had almost the same DTH as Kasalath. We developed an F2 population from a cross between BIL-55 and Nipponbare and performed QTL analysis using SSR markers. The late-heading phenotype was controlled by two known genes and at least two novel QTLs on chromosomes 4 and 6, named qDTH4.5 and qDTH6.3. These QTLs were confirmed by QTL-seq. The QTLs and polymorphisms detected here will provide useful information for further genetic studies and breeding under SD conditions at lower latitudes.

Hd1,a CONSTANS ortholog in rice, functions as an Ehd1 repressor through interaction with monocot-specific CCT-domain protein Ghd7.

Flowering time is an important agronomic trait that affects crop yields. In cereals, several CCT-domain proteins unique to monocots, such as Grain number, plant height, and heading date 7 (Ghd7) gene, have been identified as key floral repressors, although the corresponding molecular mechanisms have been unknown until now. In rice, a short-day plant, Heading date 1 (Hd1) gene, a rice ortholog of Arabidopsis floral activator CONSTANS (CO), represses flowering under non-inductive long-day (LD) conditions and induces it under inductive short-day (SD) conditions. Here, we report biological interactions between Ghd7 and Hd1, which together repress Early heading date 1 (Ehd1), a key floral inducer under non-inductive LD conditions. In addition to this genetic interaction between them, Co-IP experiments further demonstrated that a Ghd7-Hd1 protein formed a complex in vivo and ChIP and luciferase reporter analyses suggested that this complex specifically binds to a cis-regulatory region in Ehd1 and represses its expression. These findings imply that Hd1, an evolutionally conserved transcriptional activator, can function as a strong transcriptional repressor within a monocot-specific flowering-time pathway through with Ghd7.

Advanced backcross QTL analysis reveals complicated genetic control of rice grain shape in a japonica × indica cross.

Grain shape is an important trait for improving rice yield. A number of quantitative trait loci (QTLs) for this trait have been identified by using primary F2 mapping populations and recombinant inbred lines, in which QTLs with a small effect are harder to detect than they would be in advanced generations. In this study, we developed two advanced mapping populations (chromosome segment substitution lines [CSSLs] and BC4F2 lines consisting of more than 2000 individuals) in the genetic backgrounds of two improved cultivars: a japonica cultivar (Koshihikari) with short, round grains, and an indica cultivar (IR64) with long, slender grains. We compared the ability of these materials to reveal QTLs for grain shape with that of an F2 population. Only 8 QTLs for grain length or grain width were detected in the F2 population, versus 47 in the CSSL population and 65 in the BC4F2 population. These results strongly suggest that advanced mapping populations can reveal QTLs for agronomic traits under complicated genetic control, and that DNA markers linked with the QTLs are useful for choosing superior allelic combinations to enhance grain shape in the Koshihikari and IR64 genetic backgrounds.

Natural Variation in the Flag Leaf Morphology of Rice Due to a Mutation of the NARROW LEAF 1 Gene in Oryza sativa L.

We investigated the natural variations in the flag leaf morphology of rice. We conducted a principal component analysis based on nine flag leaf morphology traits using 103 accessions from the National Institute of Agrobiological Sciences Core Collection. The first component explained 39% of total variance, and the variable with highest loading was the width of the flag leaf (WFL). A genome-wide association analysis of 102 diverse Japanese accessions revealed that marker RM6992 on chromosome 4 was highly associated with WFL. In analyses of progenies derived from a cross between Takanari and Akenohoshi, the most significant quantitative trait locus (QTL) for WFL was in a 10.3-kb region containing the NARROW LEAF 1 (NAL1) gene, located 0.4 Mb downstream of RM6992. Analyses of chromosomal segment substitution lines indicated that a mutation (G1509A single-nucleotide mutation, causing an R233H amino acid substitution in NAL1) was present at the QTL. This explained 13 and 20% of total variability in WFL and the distance between small vascular bundles, respectively. The mutation apparently occurred during rice domestication and spread into japonica, tropical japonica, and indica subgroups. Notably, one accession, Phulba, had a NAL1 allele encoding only the N-terminal, or one-fourth, of the wild-type peptide. Given that the Phulba allele and the histidine-type allele showed essentially the same phenotype, the histidine-type allele was regarded as malfunctional. The phenotypes of transgenic plants varied depending on the ratio of histidine-type alleles to arginine-type alleles, raising the possibility that H(233)-type products function differently from and compete with R(233)-type products.

Genetic architecture of variation in heading date among Asian rice accessions.

BACKGROUND: Heading date, a crucial factor determining regional and seasonal adaptation in rice (Oryza sativa L.), has been a major selection target in breeding programs. Although considerable progress has been made in our understanding of the molecular regulation of heading date in rice during last two decades, the previously isolated genes and identified quantitative trait loci (QTLs) cannot fully explain the natural variation for heading date in diverse rice accessions. RESULTS: To genetically dissect naturally occurring variation in rice heading date, we collected QTLs in advanced-backcross populations derived from multiple crosses of the japonica rice accession Koshihikari (as a common parental line) with 11 diverse rice accessions (5 indica, 3 aus, and 3 japonica) that originate from various regions of Asia. QTL analyses of over 14,000 backcrossed individuals revealed 255 QTLs distributed widely across the rice genome. Among the detected QTLs, 128 QTLs corresponded to genomic positions of heading date genes identified by previous studies, such as Hd1, Hd6, Hd3a, Ghd7, DTH8, and RFT1. The other 127 QTLs were detected in different chromosomal regions than heading date genes. CONCLUSIONS: Our results indicate that advanced-backcross progeny allowed us to detect and confirm QTLs with relatively small additive effects, and the natural variation in rice heading date could result from combinations of large- and small-effect QTLs. We also found differences in the genetic architecture of heading date (flowering time) among maize, Arabidopsis, and rice.

Genetic mechanisms underlying yield potential in the rice high-yielding cultivar Takanari, based on reciprocal chromosome segment substitution lines.

BACKGROUND: Increasing rice yield potential is a major objective in rice breeding programs, given the need for meeting the demands of population growth, especially in Asia. Genetic analysis using genomic information and high-yielding cultivars can facilitate understanding of the genetic mechanisms underlying rice yield potential. Chromosome segment substitution lines (CSSLs) are a powerful tool for the detection and precise mapping of quantitative trait loci (QTLs) that have both large and small effects. In addition, reciprocal CSSLs developed in both parental cultivar backgrounds may be appropriate for evaluating gene activity, as a single factor or in epistatic interactions. RESULTS: We developed reciprocal CSSLs derived from a cross between Takanari (one of the most productive indica cultivars) and a leading japonica cultivar, Koshihikari; both the cultivars were developed in Japan. Forty-one CSSLs covered most of the Takanari genome in the Koshihikari background and 39 CSSLs covered the Koshihikari genome in the Takanari background. Using the reciprocal CSSLs, we conducted yield trials under canopy conditions in paddy fields. While no CSSLs significantly exceeded the recurrent parent cultivar in yield, genetic analysis detected 48 and 47 QTLs for yield and its components in the Koshihikari and Takanari backgrounds, respectively. A number of QTLs showed a trade-off, in which the allele with increased sink-size traits (spikelet number per panicle or per square meter) was associated with decreased ripening percentage or 1000-grain weight. These results indicate that increased sink size is not sufficient to increase rice yield in both backgrounds. In addition, most QTLs were detected in either one of the two genetic backgrounds, suggesting that these loci may be under epistatic control with other gene(s). CONCLUSIONS: We demonstrated that the reciprocal CSSLs are a useful tool for understanding the genetic mechanisms underlying yield potential in the high-yielding rice cultivar Takanari. Our results suggest that sink-size QTLs in combination with QTLs for source strength or translocation capacity, as well as careful attention to epistatic interactions, are necessary for increasing rice yield. Thus, our findings provide a foundation for developing rice cultivars with higher yield potential in future breeding programs.

Detection of QTLs to reduce cadmium content in rice grains using LAC23/Koshihikari chromosome segment substitution lines.

To advance the identification of quantitative trait loci (QTLs) to reduce Cd content in rice (Oryza sativa L.) grains and breed low-Cd cultivars, we developed a novel population consisting of 46 chromosome segment substitution lines (CSSLs) in which donor segments of LAC23, a cultivar reported to have a low grain Cd content, were substituted into the Koshihikari genetic background. The parental cultivars and 32 CSSLs (the minimum set required for whole-genome coverage) were grown in two fields with different natural levels of soil Cd. QTL mapping by single-marker analysis using ANOVA indicated that eight chromosomal regions were associated with grain Cd content and detected a major QTL (qlGCd3) with a high F-test value in both fields (F = 9.19 and 5.60) on the long arm of chromosome 3. The LAC23 allele at qlGCd3 was associated with reduced grain Cd levels and appeared to reduce Cd transport from the shoots to the grains. Fine substitution mapping delimited qlGCd3 to a 3.5-Mbp region. Our results suggest that the low-Cd trait of LAC23 is controlled by multiple QTLs, and qlGCd3 is a promising candidate QTL to reduce the Cd level of rice grain.

Natural variation of the RICE FLOWERING LOCUS T 1 contributes to flowering time divergence in rice.

In rice (Oryza sativa L.), there is a diversity in flowering time that is strictly genetically regulated. Some indica cultivars show extremely late flowering under long-day conditions, but little is known about the gene(s) involved. Here, we demonstrate that functional defects in the florigen gene RFT1 are the main cause of late flowering in an indica cultivar, Nona Bokra. Mapping and complementation studies revealed that sequence polymorphisms in the RFT1 regulatory and coding regions are likely to cause late flowering under long-day conditions. We detected polymorphisms in the promoter region that lead to reduced expression levels of RFT1. We also identified an amino acid substitution (E105K) that leads to a functional defect in Nona Bokra RFT1. Sequencing of the RFT1 region in rice accessions from a global collection showed that the E105K mutation is found only in indica, and indicated a strong association between the RFT1 haplotype and extremely late flowering in a functional Hd1 background. Furthermore, SNPs in the regulatory region of RFT1 and the E105K substitution in 1,397 accessions show strong linkage disequilibrium with a flowering time-associated SNP. Although the defective E105K allele of RFT1 (but not of another florigen gene, Hd3a) is found in many cultivars, relative rate tests revealed no evidence for differential rate of evolution of these genes. The ratios of nonsynonymous to synonymous substitutions suggest that the E105K mutation resulting in the defect in RFT1 occurred relatively recently. These findings indicate that natural mutations in RFT1 provide flowering time divergence under long-day conditions.