Whole-Genome Analysis Revealed the Positively Selected Genes during the Differentiation of indica and Temperate japonica Rice

To investigate the selective pressures acting on the protein-coding genes during the differentiation of indica and japonica, all of the possible orthologous genes between the Nipponbare and 93–11 genomes were identified and compared with each other. Among these genes, 8,530 pairs had identical sequences, and 27,384 pairs shared more than 90% sequence identity. Only 2,678 pairs of genes displaying a Ka/Ks ratio significantly greater than one were revealed, and most of these genes contained only nonsynonymous sites. The genes without synonymous site were further analyzed with the SNP data of 1529 O. sativa and O. rufipogon accessions, and 1068 genes were identified to be under positive selection during the differentiation of indica and temperate japonica. The positively selected genes (PSGs) are unevenly distributed on 12 chromosomes, and the proteins encoded by the PSGs are dominant with binding, transferase and hydrolase activities, and especially enriched in the plant responses to stimuli, biological regulations, and transport processes. Meanwhile, the most PSGs of the known function and/or expression were involved in the regulation of biotic/abiotic stresses. The evidence of pervasive positive selection suggested that many factors drove the differentiation of indica and japonica, which has already started in wild rice but is much lower than in cultivated rice. Lower differentiation and less PSGs revealed between the Or-It and Or-IIIt wild rice groups implied that artificial selection provides greater contribution on the differentiation than natural selection. In addition, the phylogenetic tree constructed with positively selected sites showed that the japonica varieties exhibited more diversity than indica on differentiation, and Or-III of O. rufipogon exhibited more than Or-I.     

Essential amino acids of starch synthase IIa differentiate amylopectin structure and starch quality between <Emphasis Type="Italic">japonica</Emphasis> and <Emphasis Type="Italic">indica</Emphasis> rice varieties

Four amino acids were variable between the ‘active’ indica-type and ‘inactive’ japonica-type soluble starch synthase IIa (SSIIa) of rice plants; Glu-88 and Gly-604 in SSIIa of indica-cultivars IR36 and Kasalath were replaced by Asp-88 and Ser-604, respectively, in both japonica cultivars Nipponbare and Kinmaze SSIIa, whereas Val-737 and Leu-781 in indica SSIIa were replaced by Met-737 in cv. Nipponbare and Phe-781 in cv. Kinmaze SSIIa, respectively. The SSIIa gene fragments shuffling experiments revealed that Val-737 and Leu-781 are essential not only for the optimal SSIIa activity, but also for the capacity to synthesize indica-type amylopectin. Surprisingly, however, a combination of Phe-781 and Gly-604 could restore about 44% of the SSIIa activity provided that Val-737 was conserved. The introduction of the ‘active’ indica-type SSIIa gene enabled the japonica-type cv. Kinmaze to synthesize indica-type amylopectin. The starch in the transformed japonica rice plants exhibited gelatinization-resistant properties that are characteristic of indica-rice starch. Transformed lines expressing different levels of the IR36 SSIIa protein produced a variety of starches with amylopectin chain-length distribution patterns that correlated well with their onset temperatures of gelatinization. The present study confirmed that the SSIIa activity determines the type of amylopectin structure of rice starch to be either the typical indica-type or japonica-type, by playing a specific role in the synthesis of the long B₁ chains by elongating short A and B₁ chains, notwithstanding the presence of functional two additional SSII genes, a single SSI gene, two SSIII genes, and two SSIV genes in rice plants.     

Genome-wide selection and introgression of Chinese rice varieties during breeding

China is the largest rice-producing country, but the genomic landscape of rice diversity has not yet been clarified. In this study, we re-sequence 1070 rice varieties collected from China (400) and other regions in Asia (670). Among the six major rice groups (aus, indica-I, indica-II, aromatic, temperate japonica, and tropical japonica), almost all Chinese varieties belong to the indica-II or temperate japonica group. Most Chinese indica varieties belong to indica-II, which consists of two subgroups developed during different phases of rice breeding. The genomic segments underlying the differences between these subgroups span 36.32 Mb. The Chinese japonica rice varieties fall into the temperate japonica group, consisting of two subgroups based on their geographical distribution. The genomic segments underlying the differences between these subgroups span 27.69 Mb. These differentiated segments in the Chinese indica varieties span 45 genes with nonsynonymous mutations that are closely related to variations in plant height and grain width. Fifty-four genes with nonsynonymous mutations are associated with the differences in heading date between the two Chinese japonica subgroups. These findings provide new insights into rice diversity in China that will facilitate the molecular breeding.     

Quantitative Trait Locus Mapping and Candidate Gene Analysis for Plant Architecture Traits Using Whole Genome Re-Sequencing in Rice

Plant breeders have focused on improving plant architecture as an effective means to increase crop yield. Here, we identify the main-effect quantitative trait loci (QTLs) for plant shape-related traits in rice (Oryza sativa) and find candidate genes by applying whole genome re-sequencing of two parental cultivars using next-generation sequencing. To identify QTLs influencing plant shape, we analyzed six traits: plant height, tiller number, panicle diameter, panicle length, flag leaf length, and flag leaf width. We performed QTL analysis with 178 F₇ recombinant in-bred lines (RILs) from a cross of japonica rice line ‘SNUSG1’ and indica rice line ‘Milyang23’. Using 131 molecular markers, including 28 insertion/deletion markers, we identified 11 main- and 16 minor-effect QTLs for the six traits with a threshold LOD value > 2.8. Our sequence analysis identified fifty-four candidate genes for the main-effect QTLs. By further comparison of coding sequences and meta-expression profiles between japonica and indica rice varieties, we finally chose 15 strong candidate genes for the 11 main-effect QTLs. Our study shows that the whole-genome sequence data substantially enhanced the efficiency of polymorphic marker development for QTL fine-mapping and the identification of possible candidate genes. This yields useful genetic resources for breeding high-yielding rice cultivars with improved plant architecture.     

Gene Expression Profiles Deciphering Rice Phenotypic Variation between Nipponbare (Japonica) and 93-11 (Indica) during Oxidative Stress

Rice is a very important food staple that feeds more than half the world''s population. Two major Asian cultivated rice (Oryza sativa L.) subspecies, japonica and indica, show significant phenotypic variation in their stress responses. However, the molecular mechanisms underlying this phenotypic variation are still largely unknown. A common link among different stresses is that they produce an oxidative burst and result in an increase of reactive oxygen species (ROS). In this study, methyl viologen (MV) as a ROS agent was applied to investigate the rice oxidative stress response. We observed that 93-11 (indica) seedlings exhibited leaf senescence with severe lesions under MV treatment compared to Nipponbare (japonica). Whole-genome microarray experiments were conducted, and 1,062 probe sets were identified with gene expression level polymorphisms between the two rice cultivars in addition to differential expression under MV treatment, which were assigned as Core Intersectional Probesets (CIPs). These CIPs were analyzed by gene ontology (GO) and highlighted with enrichment GO terms related to toxin and oxidative stress responses as well as other responses. These GO term-enriched genes of the CIPs include glutathine S-transferases (GSTs), P450, plant defense genes, and secondary metabolism related genes such as chalcone synthase (CHS). Further insertion/deletion (InDel) and regulatory element analyses for these identified CIPs suggested that there may be some eQTL hotspots related to oxidative stress in the rice genome, such as GST genes encoded on chromosome 10. In addition, we identified a group of marker genes individuating the japonica and indica subspecies. In summary, we developed a new strategy combining biological experiments and data mining to study the possible molecular mechanism of phenotypic variation during oxidative stress between Nipponbare and 93-11. This study will aid in the analysis of the molecular basis of quantitative traits.     

Development and characterization of <Emphasis Type="Italic">japonica</Emphasis> rice lines carrying the brown planthopper-resistance genes <Emphasis Type="Italic">BPH12</Emphasis> and <Emphasis Type="Italic">BPH6</Emphasis>

The brown planthopper (Nilaparvata lugens Stål; BPH) has become a severe constraint on rice production. Identification and pyramiding BPH-resistance genes is an economical and effective solution to increase the resistance level of rice varieties. All the BPH-resistance genes identified to date have been from indica rice or wild species. The BPH12 gene in the indica rice accession B14 is derived from the wild species Oryza latifolia. Using an F₂ population from a cross between the indica cultivar 93-11 and B14, we mapped the BPH12 gene to a 1.9-cM region on chromosome 4, flanked by the markers RM16459 and RM1305. In this population, BPH12 appeared to be partially dominant and explained 73.8% of the phenotypic variance in BPH resistance. A near-isogenic line (NIL) containing the BPH12 locus in the background of the susceptible japonica variety Nipponbare was developed and crossed with a NIL carrying BPH6 to generate a pyramid line (PYL) with both genes. BPH insects showed significant differences in non-preference in comparisons between the lines harboring resistance genes (NILs and PYL) and Nipponbare. BPH growth and development were inhibited and survival rates were lower on the NIL-BPH12 and NIL-BPH6 plants compared to the recurrent parent Nipponbare. PYL-BPH6 + BPH12 exhibited 46.4, 26.8 and 72.1% reductions in population growth rates (PGR) compared to NIL-BPH12, NIL-BPH6 and Nipponbare, respectively. Furthermore, insect survival rates were the lowest on the PYL-BPH6 + BPH12 plants. These results demonstrated that pyramiding different BPH-resistance genes resulted in stronger antixenotic and antibiotic effects on the BPH insects. This gene pyramiding strategy should be of great benefit for the breeding of BPH-resistant japonica rice varieties.     

Recombinational and mutational analyses of subspecific differentiation in Oryza sativa L.

Recombinational and mutational analyses were performed in an attempt to elucidate the nature of the genic differences and the ancestral relationship between the indica and japonica subspecies of Oryza sativa L. by taking representative varieties from each group. The absence of cytological aberrations at meiosis in indica×japonica hybrids suggest that genetic differences are predominantly genic rather than chromosomal. The pattern of variation induced in the key characters, and the frequency and spectra of mutations obtained in each group, failed to support the hypothesis that indica has given rise to japonica, and suggested that the two subspecies have originated from a common ancestor.     

Genetic diversity and classification of Oryza sativa with emphasis on Chinese rice germplasm

Despite extensive studies on cultivated rice, the genetic structure and subdivision of this crop remain unclear at both global and local scales. Using 84 nuclear simple sequence repeat markers, we genotyped a panel of 153 global rice cultivars covering all previously recognized groups and 826 cultivars representing the diversity of Chinese rice germplasm. On the basis of model-based grouping, neighbour-joining tree and principal coordinate analysis, we confirmed the widely accepted five major groups of rice cultivars (indica, aus, aromatic, temperate japonica and tropical japonica), and demonstrated that rayada rice was unique in genealogy and should be treated as a new (the sixth) major group of rice germplasm. With reference to the global classification of rice cultivars, we identified three major groups (indica, temperate japonica and tropical japonica) in Chinese rice germplasm and showed that Chinese temperate japonica contained higher diversity than that of global samples, whereas Chinese indica and tropical japonica maintained slightly lower diversity than that present in the global samples. Particularly, we observed that all seasonal, drought-tolerant and endosperm types occurred within each of three major groups of Chinese cultivars, which does not support previous claims that seasonal differentiation exists in Indica and drought-tolerant differentiation is present in Japonica. It is most likely that differentiation of cultivar types arose multiple times stemming from artificial selection for adaptation to local environments.     

Construction of Pseudomolecule Sequences of the aus Rice Cultivar Kasalath for Comparative Genomics of Asian Cultivated Rice

Having a deep genetic structure evolved during its domestication and adaptation, the Asian cultivated rice (Oryza sativa) displays considerable physiological and morphological variations. Here, we describe deep whole-genome sequencing of the aus rice cultivar Kasalath by using the advanced next-generation sequencing (NGS) technologies to gain a better understanding of the sequence and structural changes among highly differentiated cultivars. The de novo assembled Kasalath sequences represented 91.1% (330.55 Mb) of the genome and contained 35 139 expressed loci annotated by RNA-Seq analysis. We detected 2 787 250 single-nucleotide polymorphisms (SNPs) and 7393 large insertion/deletion (indel) sites (>100 bp) between Kasalath and Nipponbare, and 2 216 251 SNPs and 3780 large indels between Kasalath and 93-11. Extensive comparison of the gene contents among these cultivars revealed similar rates of gene gain and loss. We detected at least 7.39 Mb of inserted sequences and 40.75 Mb of unmapped sequences in the Kasalath genome in comparison with the Nipponbare reference genome. Mapping of the publicly available NGS short reads from 50 rice accessions proved the necessity and the value of using the Kasalath whole-genome sequence as an additional reference to capture the sequence polymorphisms that cannot be discovered by using the Nipponbare sequence alone.     

Genomic regions involved in yield potential detected by genome-wide association analysis in Japanese high-yielding rice cultivars

Background

High-yielding cultivars of rice (Oryza sativa L.) have been developed in Japan from crosses between overseas indica and domestic japonica cultivars. Recently, next-generation sequencing technology and high-throughput genotyping systems have shown many single-nucleotide polymorphisms (SNPs) that are proving useful for detailed analysis of genome composition. These SNPs can be used in genome-wide association studies to detect candidate genome regions associated with economically important traits. In this study, we used a custom SNP set to identify introgressed chromosomal regions in a set of high-yielding Japanese rice cultivars, and we performed an association study to identify genome regions associated with yield.

Results

An informative set of 1152 SNPs was established by screening 14 high-yielding or primary ancestral cultivars for 5760 validated SNPs. Analysis of the population structure of high-yielding cultivars showed three genome types: japonica-type, indica-type and a mixture of the two. SNP allele frequencies showed several regions derived predominantly from one of the two parental genome types. Distinct regions skewed for the presence of parental alleles were observed on chromosomes 1, 2, 7, 8, 11 and 12 (indica) and on chromosomes 1, 2 and 6 (japonica). A possible relationship between these introgressed regions and six yield traits (blast susceptibility, heading date, length of unhusked seeds, number of panicles, surface area of unhusked seeds and 1000-grain weight) was detected in eight genome regions dominated by alleles of one parental origin. Two of these regions were near Ghd7, a heading date locus, and Pi-ta, a blast resistance locus. The allele types (i.e., japonica or indica) of significant SNPs coincided with those previously reported for candidate genes Ghd7 and Pi-ta.

Conclusions

Introgression breeding is an established strategy for the accumulation of QTLs and genes controlling high yield. Our custom SNP set is an effective tool for the identification of introgressed genome regions from a particular genetic background. This study demonstrates that changes in genome structure occurred during artificial selection for high yield, and provides information on several genomic regions associated with yield performance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-346) contains supplementary material, which is available to authorized users.     

Construction of a Microsatellite Linkage Map with Two Sequenced Rice Varieties

Based on the successful development of new microsatellite markers from the data of two whole-sequenced rice varieties, japonica variety Nipponbare and indica variety 9311, an F₂ population of 90 lines, which was derived from a single cross between Nipponbare and 9311, was applied to construct a genetic linkage framework map. The map covered 2 455.7 cM of total genomic length, and consisted of 152 simple sequence repeats (SSRs) loci including 46 pairs of new SSR primers developed by our research institute. The average genetic distance between two markers was 16.16 cM. In addition, markers RM345 and RM494, which have not been mapped on the Temnykh's map et al. (2001) were anchored on the sixth chromosome of this map. We compared this research with maps of Temnykh et al.(2001) and LAN et al. (2003) regarding the aspects of type and size of population, type and quantity of markers, and the marker arrangement order on chromosome, etc. Results indicated that the similarity of marker linear alignment was 93.81% between this map and T-map, Finally, the important significance of using sequenced rice varieties to construct linkage map was also discussed.     

Cloning and Mapping of Telomere-Associated Sequences from Rice

We have isolated three telomere-associated sequences from riceusing cassette-ligation-mediated polymerase chain reaction (PCR).Each of the obtained clones hybridized to the terminal of oneor several rice chromosome arms. The telomeres recognized bythe clones displayed a high level of polymorphism between tworice varieties, Nipponbare (a japonica variety) and Kasalath(an indica variety). Variability in the chromosome termini wasalso detected among individual F₂ progeny plants, which werederived from a cross between the two rice varieties. One clonecontaining telomere-associated sequences was located to oneend of chromosome 5, and another clone to one end of chromosome11. For another clone, non-allelic segregation of polymorphichybridization bands was observed between japonica and indicarice; this clone was mapped to one end of chromosome 12 in japonicaand to one end of chromosome 11 in indica rice. This indicatesan exchange of termini between nonhomologous chromosomes.     

Non-specific resistance in brown planthopper resistant indica rice varieties against Plodia interpunctella (Lepidoptera: Phycitidae)

Mechanisms of host plant resistance against insect pests can be manifold. Resistance screenings generally use single target insect pests, but the resistance thus screened may not always be specific to the target insect species. We conducted a test for non‐specific resistance in indica rice varieties with resistance genes against brown planthopper (BPH), by using the Indian meal moth, Plodia interpunctella. The test system was very simple, and only required the non‐pest moth to be reared on rice flour. We compared the survival rate, developmental period and adult weight of the moth on three rice varieties: ‘Nipponbare’, a BPH‐susceptible japonica variety, and ‘Thai Collection 11’ and ‘Pokkali’, two resistant indica varieties. Our results were straightforward and demonstrate that resistance in the two resistant rice varieties is not BPH specific, because development of the moth was retarded and adult body weight was reduced.     

Genome-Wide Distribution,Organisation and Functional Characterization of Disease Resistance and Defence Response Genes across Rice Species

The resistance (R) genes and defense response (DR) genes have become very important resources for the development of disease resistant cultivars. In the present investigation, genome-wide identification, expression, phylogenetic and synteny analysis was done for R and DR-genes across three species of rice viz: Oryza sativa ssp indica cv 93-11, Oryza sativa ssp japonica and wild rice species, Oryza brachyantha. We used the in silico approach to identify and map 786 R -genes and 167 DR-genes, 672 R-genes and 142 DR-genes, 251 R-genes and 86 DR-genes in the japonica, indica and O. brachyanth a genomes, respectively. Our analysis showed that 60.5% and 55.6% of the R-genes are tandemly repeated within clusters and distributed over all the rice chromosomes in indica and japonica genomes, respectively. The phylogenetic analysis along with motif distribution shows high degree of conservation of R- and DR-genes in clusters. In silico expression analysis of R-genes and DR-genes showed more than 85% were expressed genes showing corresponding EST matches in the databases. This study gave special emphasis on mechanisms of gene evolution and duplication for R and DR genes across species. Analysis of paralogs across rice species indicated 17% and 4.38% R-genes, 29% and 11.63% DR-genes duplication in indica and Oryza brachyantha, as compared to 20% and 26% duplication of R-genes and DR-genes in japonica respectively. We found that during the course of duplication only 9.5% of R- and DR-genes changed their function and rest of the genes have maintained their identity. Syntenic relationship across three genomes inferred that more orthology is shared between indica and japonica genomes as compared to brachyantha genome. Genome wide identification of R-genes and DR-genes in the rice genome will help in allele mining and functional validation of these genes, and to understand molecular mechanism of disease resistance and their evolution in rice and related species.     

Pyramiding transgenes for multiple resistance in rice against bacterial blight,yellow stem borer and sheath blight

The present investigation revealed that the alk and gel(t) genes, which cause the differences between a japonica rice variety Nipponbare and an indica rice variety Kasalath in terms of the disintegration of endosperm starch granules in alkali solution and their gelatinisation in a 4 M urea solution, respectively, cosegregated in backcross inbred lines derived from a cross between the two varieties. The segregation pattern of the profile for amylopectin chain-length, which was distinguished by enrichment in short chains of DP≦11 and depletion in intermediate-size chains of 12≦DP≦24 in japonica as compared with indica, was exactly the same as those of the above physico-chemical properties of starch granules, and the gene was designated as acl(t). Gene-mapping analysis showed that the starch synthase IIa (SSIIa) gene is located at the alk locus on chromosome 6 in the rice genome. These results lead us to the possibility that different alleles of the SSIIa gene are responsible for differences in amylopectin structure between the two varieties, in that SSIIa plays a distinct role in the elongation of short chains within clusters (A+B₁ chains) of amylopectin. It is proposed that the activity of SSIIa in japonica rice is reduced in amount or functional capacity relative to the activity of this enzyme in indica rice. This, in turn, would explain why starch from japonica rice has a lower gelatinisation temperature than starch from indica rice and is more susceptible to disintegration in alkali or urea. The evidence for this hypothesis is that the alk(t), gel(t), acl(t) and SSIIa genes all map to the same locus. Received: 29 January 2001 / Accepted: 12 April 2001     

杂草稻叶绿体籼粳分化的多重PCR分析

通过分析籼稻93-11和粳稻培矮64S的叶绿体全基因组,优化和构建了籼粳分化的叶绿体分子标记ORF100和ORF29-TrnC^GCA的多重PCR。应用这个多重PCR对200余份世界各地杂草稻和其它水稻材料进行分析。结果表明:杂草稻中有明显的叶绿体籼粳分化,表现出明显的地域性,且与传统的中国栽培稻的南籼北粳能较好的对应。推测粳型杂草稻可能是栽培稻突变或粳型水稻(作母本)与其它类型水稻材料杂交而形成的。     

Genome-wide Investigation of Intron Length Polymorphisms and Their Potential as Molecular Markers in Rice (Oryza sativa L.)

Intron length polymorphisms (ILPs) have been used as geneticmarkers in some studies. However, a systematic investigationand large-scale exploitation of ILP markers has not been reported.In this study, we performed a genome-wide search of ILPs betweentwo subspecies (indica and japonica) in rice using the draftgenomic sequences of cultivars 93-11 (indica) and Nipponbare(japonica) and 32 127 full-length cDNA sequences of Nipponbareobtained from public databases. We identified 13 308 putativeILPs. Based on these putative ILPs, we developed 5811 candidateILP markers via electronic-PCR with primers designed in flankingexons. We further conducted experiment to verify the candidateILP markers. Out of 215 candidate ILP markers tested on 93-11,Nipponbare and their hybrid, we successfully exploited 173 codominantILP markers. Further analyses on 10 rice accessions showed thatthese ILP markers were widely applicable and most (71.1%) exhibitedsubspecies specificity. This feature suggests that ILPs wouldbe useful for the studies of genome evolution and inter-subspeciesheterosis and for cross-subspecies marker-assisted selectionin rice. In addition, by testing 51 pairs of the ILP primerson five Gramineae plants and three dicot plants, we found anotherdesirable characteristic of rice ILP markers that they havehigh transferability to other plants.