Genome-wide DNA polymorphisms in seven rice cultivars of temperate and tropical japonica groups

A combination of upstream alleles involved in rice heading hastens natural long-day responses

BACKGROUND

The female parental line Jinbuol (JBO, early heading) and two recombinant isogenic lines, JSRIL1 and JSRIL2, have been shown to flower 44, 34 and 16 days earlier, respectively, than the male parental line Samgwang (SG, late heading) in paddy fields.

OBJECTIVE

To explore how photoperiodicity-related genes are involved in differential heading among these lines.

METHODS

Deep sequencing was conducted for these lines, photoperiodicity-related genes (71) were categorized, and qRT-PCR was performed for some key genes.

RESULTS

Deep sequencing revealed a nearly even contribution of parental groups, with 48.5% and 45% of the chromosomes in JSRIL1 and JSRIL2, respectively, inherited from the female parent JBO; however, Chr6 contained the most biased parental contribution, with 99.4% inherited from the female parent. The variation in single-nucleotide polymorphisms (SNPs) among many known flower-inducing genes, including rice GIGANTEA (OsGI); grain number, plant height and heading date 7 (Ghd7); and EARLY HEADING DATE 1 (Ehd1), was minimal. In the JSRILs, HEADING DATE 1 (Hd1) and VERNALIZATION INSENSITIVE 3-LIKE 1 (OsVIL2) originated from JBO, whereas FLAVIN-BINDING, KELCH REPEAT, F BOX 1 (OsFKF1) originated from SG. Interestingly, HEN1 suppressor 1 (OsHESO1) originated from SG in JSRIL1 and JBO in JSRIL2. RNA sequencing and qRT‒PCR analyses of plants at the floral meristem stage revealed that transcriptional regulation through chromosomal restructuring and posttranscriptional regulation might control minute gene regulation, resulting in delayed heading in JSRILs.

CONCLUSION

Our gene expression and SNP analyses of elite recombinant isogenic lines could be helpful in understanding how photoperiodicity-related genes in rice are modulated.

Hybrid-derived weedy rice maintains adaptive combinations of alleles associated with seed dormancy

Plant hybridization is a pathway for the evolution of adaptive traits. However, hybridization between adapted and nonadapted populations may affect the persistence of combinations of adaptive alleles evolved through natural selection. Seed dormancy is an adaptive trait for weedy rice because it regulates the timing of seed germination and the persistence of the soil seed bank. Hybridization between weedy and cultivated rice has been confirmed with an adaptive introgression of deep seed dormancy alleles from cultivated rice. Here, we explored the influence of hybridization on the conservation of adaptive allele combinations by evaluating natural variation and genetic structure in seed dormancy-associated genomic regions. Based on sequence variation in the genomic regions associated with seed dormancy, hybrid-derived weedy rice strains maintained most of the adaptive combinations for this trait observed in the parental weedy rice, despite equal representation of the parental weedy and cultivated rice in the whole genome sequence. Moreover, hybrid-derived weedy rice strains were more dormant than their parental weedy rice strains, and this trait was strongly influenced by the environment. This study suggests that hybridization between weedy rice (adaptive allelic combinations for seed dormancy) and cultivated rice (nonadaptive combinations) generates weedy rice strains expressing deep seed dormancy caused by genome stabilization through the removal of alleles derived from cultivated rice, in addition to the adaptive introgression of deep seed dormancy alleles derived from cultivated rice. Thus, hybridization between adapted and nonadapted populations appears to be reinforcing the trajectory towards the evolution of adaptive traits.

Genomic Variation in Korean japonica Rice Varieties

Next-generation sequencing technologies have enabled the discovery of numerous sequence variations among closely related crop varieties. We analyzed genome resequencing data from 24 Korean temperate japonica rice varieties and discovered 954,233 sequence variations, including 791,121 single nucleotide polymorphisms (SNPs) and 163,112 insertions/deletions (InDels). On average, there was one variant per 391 base-pairs (bp), a variant density of 2.6 per 1 kbp. Of the InDels, 10,860 were longer than 20 bp, which enabled conversion to markers resolvable on an agarose gel. The effect of each variant on gene function was predicted using the SnpEff program. The variants were categorized into four groups according to their impact: high, moderate, low, and modifier. These groups contained 3524 (0.4%), 27,656 (2.9%), 24,875 (2.6%), and 898,178 (94.1%) variants, respectively. To test the accuracy of these data, eight InDels from a pre-harvest sprouting resistance QTL (qPHS11) target region, four highly polymorphic InDels, and four functional sequence variations in known agronomically important genes were selected and successfully developed into markers. These results will be useful to develop markers for marker-assisted selection, to select candidate genes in map-based cloning, and to produce efficient high-throughput genome-wide genotyping systems for Korean temperate japonica rice varieties.

Genomic divergence during feralization reveals both conserved and distinct mechanisms of parallel weediness evolution

Agricultural weeds are the most important biotic constraints to global crop production, and chief among these is weedy rice. Despite increasing yield losses from weedy rice in recent years worldwide, the genetic basis of weediness evolution remains unclear. Using whole-genome sequence analyses, we examined the origins and adaptation of Japanese weedy rice. We find evidence for a weed origin from tropical japonica crop ancestry, which has not previously been documented in surveys of weedy rice worldwide. We further show that adaptation occurs largely through different genetic mechanisms between independently-evolved temperate japonica- and tropical japonica-derived strains; most genomic signatures of positive selection are unique within weed types. In addition, some weedy rice strains have evolved through hybridization between weedy and cultivated rice with adaptive introgression from the crop. Surprisingly, introgression from cultivated rice confers not only crop-like adaptive traits (such as shorter plant height, facilitating crop mimicry) but also weedy-like traits (such as seed dormancy). These findings reveal how hybridization with cultivated rice can promote persistence and proliferation of weedy rice.

Translation of continuous artificial selection on phenotype into genotype during rice breeding programs

Understanding genetic diversity among local populations is a primary goal of modern crop breeding programs. Here, we demonstrated the genetic relationships of rice varieties in Hokkaido, Japan, one of the northern limits of rice cultivation around the world. Furthermore, artificial selection during rice breeding programs has been characterized using genome sequences. We utilized 8,565 single nucleotide polymorphisms and insertion/deletion markers distributed across the genome in genotype-by-sequencing for genetic diversity analyses. Phylogenetics, genetic population structure, and principal component analysis showed that a total of 110 varieties were classified into four distinct clusters according to different populations geographically and historically. Furthermore, the genome sequences of 19 rice varieties along with historic representations in Hokkaido, nucleotide diversity and F_ST values in each cluster revealed that artificial selection of elite phenotypes focused on chromosomal regions. These results clearly demonstrated the history of the selections on agronomic traits as genome sequences among current rice varieties from Hokkaido.

Landraces of temperate japonica rice have superior alleles for improving culm strength associated with lodging resistance

Lodging can reduce grain yield and quality in cereal crops including rice (Oryza sativa L.). To achieve both high biomass production and lodging resistance, the breeding of new cultivars with strong culms is a promising strategy. However, little is known about the diversity of culm strength in temperate japonica rice and underlying genetic factors. Here, we report a wide variation of culm strength among 135 temperate japonica cultivars, and some landraces having the strongest culms among these cultivars. The genome-wide association study (GWAS) identified 55 quantitative trait loci for culm strength and morphological traits, and revealed several candidate genes. The superior allele of candidate gene for culm thickness, OsRLCK191, was found in many landraces but had not inherited to the modern improved cultivars. Our results suggest that landraces of temperate japonica rice have unutilized superior alleles for contributing future improvements of culm strength and lodging resistance.

Discovery of Functional SNPs via Genome-Wide Exploration of Malaysian Pigmented Rice Varieties

Recently, rice breeding program has shown increased interests on the pigmented rice varieties due to their benefits to human health. However, the genetic variation of pigmented rice varieties is still scarce and remains unexplored. Hence, we performed genome-wide SNP analysis from the genome resequencing of four Malaysian pigmented rice varieties, representing two black and two red rice varieties. The genome of four pigmented varieties was mapped against Nipponbare reference genome sequences, and 1.9 million SNPs were discovered. Of these, 622 SNPs with polymorphic sites were identified in 258 protein-coding genes related to metabolism, stress response, and transporter. Comparative analysis of 622 SNPs with polymorphic sites against six rice SNP datasets from the Ensembl Plants variation database was performed, and 70 SNPs were identified as novel SNPs. Analysis of SNPs in the flavonoid biosynthetic genes revealed 40 nonsynonymous SNPs, which has potential as molecular markers for rice seed colour identification. The highlighted SNPs in this study show effort in producing valuable genomic resources for application in the rice breeding program, towards the genetic improvement of new and improved pigmented rice varieties.

Genome-wide analysis of DNA polymorphisms, the methylome and transcriptome revealed that multiple factors are associated with low pollen fertility in autotetraploid rice

Autotetraploid rice is a useful germplasm with high biomass production; however, low fertility is the main barrier in commercial utilization. In our previous study, differential expression of meiosis-related miRNAs was found to be involved in the pollen sterility of autotetraploid rice. However, genome-wide DNA variations and methylomes associated with low fertility of autotetraploid rice are still poorly understood. Here, we measured both global DNA variations and the methylome and compared them with the transcriptome during pollen development in autotetraploid rice by high-throughput sequencing. A total of 34416 SNPs, 6993 InDels, 1003 SVs and 25 CNVs were detected, and 11367 and 41117 differentially methylated regions showed hypermethylation and hypomethylation in 02428-4x. In total, 1110 genes displayed differentially expression in 02428-4x during meiosis, of these six harbored CNVs, including four upregulated genes with gain CNVs, such as LOC_Os11g38620. We identified 122 genes by comparing with the previous data that might be associated with low fertility during pollen development in 02428-4x. Of the 122 gens, 98 were displayed methylation and differential expression, including OsMADS98, CYP703A3 and OsABCG26. The downregulation of these three genes were confirmed by qPCR during meiosis of 02428-4x, which played pivotal roles in pollen fertility. These results indicate that the low fertility of autotetraploid rice is not only caused by the differential expression of genes involved in pollen development, but also by sequence variation and differential methylation, suggesting that the reason for pollen sterility in autotetraploid rice is complex and might be affected by multiple factors.

Tracking the genome-wide outcomes of a transposable element burst over decades of amplification

Rice (Oryza sativa) has a unique combination of attributes that made it an ideal host to track the natural behavior of very active transposable elements (TEs) over generations. In this study, we have exploited its small genome and propagation by self or sibling pollination to identify and characterize two strain pairs, EG4/HEG4 and A119/A123, undergoing bursts of the nonautonomous miniature inverted repeat transposable element mPing. Comparative sequence analyses of these strains have advanced our understanding of (i) factors that contribute to sustaining a TE burst for decades, (ii) features that distinguish a natural TE burst from bursts in cell culture or mutant backgrounds, and (iii) the extent to which TEs can rapidly diversify the genome of an inbred organism.

To understand the success strategies of transposable elements (TEs) that attain high copy numbers, we analyzed two pairs of rice (Oryza sativa) strains, EG4/HEG4 and A119/A123, undergoing decades of rapid amplification (bursts) of the class 2 autonomous Ping element and the nonautonomous miniature inverted repeat transposable element (MITE) mPing. Comparative analyses of whole-genome sequences of the two strain pairs validated that each pair has been maintained for decades as inbreds since divergence from their respective last common ancestor. Strains EG4 and HEG4 differ by fewer than 160 SNPs and a total of 264 new mPing insertions. Similarly, strains A119 and A123 exhibited about half as many SNPs (277) as new mPing insertions (518). Examination of all other potentially active TEs in these genomes revealed only a single new insertion out of ∼40,000 loci surveyed. The virtual absence of any new TE insertions in these strains outside the mPing bursts demonstrates that the Ping/mPing family gradually attains high copy numbers by maintaining activity and evading host detection for dozens of generations. Evasion is possible because host recognition of mPing sequences appears to have no impact on initiation or maintenance of the burst. Ping is actively transcribed, and both Ping and mPing can transpose despite methylation of terminal sequences. This finding suggests that an important feature of MITE success is that host recognition does not lead to the silencing of the source of transposase.

Analysis of Genes with Alternatively Spliced Transcripts in the Leaf, Root, Panicle and Seed of Rice Using a Long Oligomer Microarray and RNA-Seq

Pre-mRNA splicing further increases protein diversity acquired through evolution. The underlying driving forces for this phenomenon are unknown, especially in terms of gene expression. A rice alternatively spliced transcript detection microarray (ASDM) and RNA sequencing (RNA-Seq) were applied to differentiate the transcriptome of 4 representative organs of Oryza sativa L. cv. Ilmi: leaves, roots, 1-cm-stage panicles and young seeds at 21 days after pollination. Comparison of data obtained by microarray and RNA-Seq showed a bell-shaped distribution and a co-lineation for highly expressed genes. Transcripts were classified according to the degree of organ enrichment using a coefficient value (CV, the ratio of the standard deviation to the mean values): highly variable (CVI), variable (CVII), and constitutive (CVIII) groups. A higher index of the portion of loci with alternatively spliced transcripts in a group (IAST) value was observed for the constitutive group. Genes of the highly variable group showed the characteristics of the examined organs, and alternatively spliced transcripts tended to exhibit the same organ specificity or less organ preferences, with avoidance of 'organ distinctness'. In addition, within a locus, a tendency of higher expression was found for transcripts with a longer coding sequence (CDS), and a spliced intron was the most commonly found type of alternative splicing for an extended CDS. Thus, pre-mRNA splicing might have evolved to retain maximum functionality in terms of organ preference and multiplicity.

Genome wide re-sequencing of newly developed Rice Lines from common wild rice (Oryza rufipogon Griff.) for the identification of NBS-LRR genes

Common wild rice (Oryza rufipogon Griff.) is an important germplasm for rice breeding, which contains many resistance genes. Re-sequencing provides an unprecedented opportunity to explore the abundant useful genes at whole genome level. Here, we identified the nucleotide-binding site leucine-rich repeat (NBS-LRR) encoding genes by re-sequencing of two wild rice lines (i.e. Huaye 1 and Huaye 2) that were developed from common wild rice. We obtained 128 to 147 million reads with approximately 32.5-fold coverage depth, and uniquely covered more than 89.6% (> = 1 fold) of reference genomes. Two wild rice lines showed high SNP (single-nucleotide polymorphisms) variation rate in 12 chromosomes against the reference genomes of Nipponbare (japonica cultivar) and 93-11 (indica cultivar). InDels (insertion/deletion polymorphisms) count-length distribution exhibited normal distribution in the two lines, and most of the InDels were ranged from -5 to 5 bp. With reference to the Nipponbare genome sequence, we detected a total of 1,209,308 SNPs, 161,117 InDels and 4,192 SVs (structural variations) in Huaye 1, and 1,387,959 SNPs, 180,226 InDels and 5,305 SVs in Huaye 2. A total of 44.9% and 46.9% genes exhibited sequence variations in two wild rice lines compared to the Nipponbare and 93-11 reference genomes, respectively. Analysis of NBS-LRR mutant candidate genes showed that they were mainly distributed on chromosome 11, and NBS domain was more conserved than LRR domain in both wild rice lines. NBS genes depicted higher levels of genetic diversity in Huaye 1 than that found in Huaye 2. Furthermore, protein-protein interaction analysis showed that NBS genes mostly interacted with the cytochrome C protein (Os05g0420600, Os01g0885000 and BGIOSGA038922), while some NBS genes interacted with heat shock protein, DNA-binding activity, Phosphoinositide 3-kinase and a coiled coil region. We explored abundant NBS-LRR encoding genes in two common wild rice lines through genome wide re-sequencing, which proved to be a useful tool to exploit elite NBS-LRR genes in wild rice. The data here provide a foundation for future work aimed at dissecting the genetic basis of disease resistance in rice, and the two wild rice lines will be useful germplasm for the molecular improvement of cultivated rice.

Whole-Genome Characteristics and Polymorphic Analysis of Vietnamese Rice Landraces as a Comprehensive Information Resource for Marker-Assisted Selection

Next generation sequencing technologies have provided numerous opportunities for application in the study of whole plant genomes. In this study, we present the sequencing and bioinformatic analyses of five typical rice landraces including three indica and two japonica with potential blast resistance. A total of 688.4 million 100 bp paired-end reads have yielded approximately 30-fold coverage to compare with the Nipponbare reference genome. Among them, a small number of reads were mapped to both chromosomes and organellar genomes. Over two million and eight hundred thousand single nucleotide polymorphisms (SNPs) and insertions and deletions (InDels) in indica and japonica lines have been determined, which potentially have significant impacts on multiple transcripts of genes. SNP deserts, contiguous SNP-low regions, were found on chromosomes 1, 4, and 5 of all genomes of rice examined. Based on the distribution of SNPs per 100 kilobase pairs, the phylogenetic relationships among the landraces have been constructed. This is the first step towards revealing several salient features of rice genomes in Vietnam and providing significant information resources to further marker-assisted selection (MAS) in rice breeding programs.

Development of genome-wide PCR-based markers from insertion, deletion and single nucleotide polymorphisms for closely related Japanese rice cultivars and identification of QTLs for the appearance of cooked rice and polished rice

Appearance of rice grain is an important property, affecting its acceptance by consumers. Moreover, appearance is a complex characteristic involving many components, including glossiness and whiteness. The genetic bases for the glossiness of cooked rice and the whiteness of polished rice (WPR) were determined using 133 recombinant inbred lines (RILs) derived from a cross between two closely related cultivars from Hokkaido, Joiku462, with high glossiness and whiteness, and Yukihikari, an ancestor of Joiku462 with low glossiness and whiteness. Analyses identified 167 genome-wide InDel markers, five cleaved amplified polymorphic sequences (CAPS) and eight derived CAPS markers differentiating the parental lines. The glossiness area (GLA) and glossiness strength (GLS) of cooked rice and WPR were determined for RILs in two locations, Pippu and Sapporo, Hokkaido. Four QTLs were detected. qGLA10 and qGLS9 were detected on chromosomes 10 and 9, respectively, with both being significant at both geographic locations. qWPR1 on chromosome 1 was significant at Pippu, and qWPR4 on chromosome 4 was significant at Sapporo. The Joiku462 alleles at all QTLs increased each trait. The PCR-based markers flanking these four QTLs may be useful for improvement of GLA, GLS and WPR.

Nucleotide diversity analysis highlights functionally important genomic regions

We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3′ UTRs, the area rich with regulatory regions.

First-generation HapMap in Cajanus spp. reveals untapped variations in parental lines of mapping populations

Whole genome re-sequencing (WGRS) was conducted on a panel of 20 Cajanus spp. accessions (crossing parentals of recombinant inbred lines, introgression lines, multiparent advanced generation intercross and nested association mapping population) comprising of two wild species and 18 cultivated species accessions. A total of 791.77 million paired-end reads were generated with an effective mapping depth of ~12X per accession. Analysis of WGRS data provided 5 465 676 genome-wide variations including 4 686 422 SNPs and 779 254 InDels across the accessions. Large structural variations in the form of copy number variations (2598) and presence and absence variations (970) were also identified. Additionally, 2 630 904 accession-specific variations comprising of 2 278 571 SNPs (86.6%), 166 243 deletions (6.3%) and 186 090 insertions (7.1%) were also reported. Identified polymorphic sites in this study provide the first-generation HapMap in Cajanus spp. which will be useful in mapping the genomic regions responsible for important traits.

Genetic shift in local rice populations during rice breeding programs in the northern limit of rice cultivation in the world

The rapid accumulation of pre-existing mutations may play major roles in the establishment and shaping of adaptability for local regions in current rice breeding programs. The cultivated rice, Oryza sativa L., which originated from tropical regions, is now grown worldwide due to the concerted efforts of breeding programs. However, the process of establishing local populations and their origins remain unclear. In the present study, we characterized DNA polymorphisms in the rice variety KITAAKE from Hokkaido, one of the northern limits of rice cultivation in the world. Indel polymorphisms were attributed to transposable element-like insertions, tandem duplications, and non-TE deletions as the original mutation events in the NIPPONBARE and KITAAKE genomes. The allele frequencies of the KITAAKE alleles markedly shifted to the current variety types among the local population from Hokkaido in the last two decades. The KITAAKE alleles widely distributed throughout wild rice and cultivated rice over the world. These have accumulated in the local population from Hokkaido via Japanese landraces as the ancestral population of Hokkaido. These results strongly suggested that combinations of pre-existing mutations played a role in the establishment of adaptability. This approach using the re-sequencing of local varieties in unique environmental conditions will be useful as a genetic resource in plant breeding programs in local regions.

Whole genome sequencing and analysis of Swarna, a widely cultivated indica rice variety with low glycemic index

Swarna is a popular cultivated indica rice variety with low glycemic index (GI) but its genetic basis is not known. The whole genome of Swarna was sequenced using Illumina’s paired-end technology, and the reads were mapped to the Nipponbare reference genome. Overall, 65,984 non-synonymous SNPs were identified in 20,350 genes, and in silico analysis predicted that 4,847 of them in 2,214 genes may have deleterious effect on protein functions. Polymorphisms were found in all the starch biosynthesis genes, except the gene for branching enzyme IIa. It was found that T/G SNP at position 246, ‘A’ at position 2,386, and ‘C’ at position 3,378 in the granule bound starch synthase I gene, and C/T SNP at position 1,188 in the glucose-6-phosphate translocator gene may contribute to the low GI phenotype in Swarna. All these variants were also found in the genome of another low GI indica rice variety from Columbia, Fedearroz 50. The whole genome analysis of Swarna helped to understand the genetic basis of GI in rice, which is a complex trait involving multiple factors.

Genome-wide indel markers shared by diverse Asian rice cultivars compared to Japanese rice cultivar 'Koshihikari'

Insertion-deletion (indel) polymorphisms, such as simple sequence repeats, have been widely used as DNA markers to identify QTLs and genes and to facilitate rice breeding. Recently, next-generation sequencing has produced deep sequences that allow genome-wide detection of indels. These polymorphisms can potentially be used to develop high-accuracy polymerase chain reaction (PCR)-based markers. Here, re-sequencing of 5 indica, 2 aus, and 3 tropical japonica cultivars and Japanese elite cultivar 'Koshihikari' was performed to extract regions containing large indels (10-51 bp) shared by diverse cultivars. To design indel markers for the discrimination of genomic regions between 'Koshihikari' and other diverse cultivars, we subtracted the indel regions detected in 'Koshihikari' from those shared in other cultivars. Two sets of indel markers, KNJ8-indel (shared in eight or more cultivars, including 'Khao Nam Jen' as a representative tropical japonica cultivar) and C5-indel (shared in five to eight cultivars), were established, with 915 and 9,899 indel regions, respectively. Validation of the two marker sets by using 23 diverse cultivars showed a high PCR success rate (≥95%) for 83.3% of the KNJ8-indel markers and 73.9% of the C5-indel markers. The marker sets will therefore be useful for the effective breeding of Japanese rice cultivars.

Whole genome sequencing of elite rice cultivars as a comprehensive information resource for marker assisted selection

Current advances in sequencing technologies and bioinformatics revealed the genomic background of rice, a staple food for the poor people, and provided the basis to develop large genomic variation databases for thousands of cultivars. Proper analysis of this massive resource is expected to give novel insights into the structure, function, and evolution of the rice genome, and to aid the development of rice varieties through marker assisted selection or genomic selection. In this work we present sequencing and bioinformatics analyses of 104 rice varieties belonging to the major subspecies of Oryza sativa. We identified repetitive elements and recurrent copy number variation covering about 200 Mbp of the rice genome. Genotyping of over 18 million polymorphic locations within O. sativa allowed us to reconstruct the individual haplotype patterns shaping the genomic background of elite varieties used by farmers throughout the Americas. Based on a reconstruction of the alleles for the gene GBSSI, we could identify novel genetic markers for selection of varieties with high amylose content. We expect that both the analysis methods and the genomic information described here would be of great use for the rice research community and for other groups carrying on similar sequencing efforts in other crops.