Evolution and Diversity of the Wild Rice Oryza officinalis Complex, across Continents, Genome Types, and Ploidy Levels

Oryza genome evolution through a tetraploid lens

Oryza is a remarkable genus comprising 27 species and 11 genome types, with ~3.4-fold genome size variation, that possesses a virtually untapped reservoir of genes that can be used for crop improvement and neodomestication. Here we present 11 chromosome-level assemblies (nine tetraploid, two diploid) in the context of ~15 million years of evolution and show that the core Oryza (sub)genome is only ~200 Mb and largely syntenic, whereas the remaining nuclear fractions (~80-600 Mb) are intermingled, plastic and rapidly evolving. For the halophyte Oryza coarctata, we found that despite detection of gene fractionation in the subgenomes, homoeologous genes were expressed at higher levels in one subgenome over the other in a mosaic form, demonstrating subgenome equivalence. The integration of these 11 new reference genomes with previously published genome datasets provides a nearly complete view of the consequences of evolution for genome diversification across the genus.

A pangenome reference of wild and cultivated rice

Oryza rufipogon, the wild progenitor of Asian cultivated rice Oryza sativa, is an important resource for rice breeding1. Here we present a wild-cultivated rice pangenome based on 145 chromosome-level assemblies, comprising 129 genetically diverse O. rufipogon accessions and 16 diverse varieties of O. sativa. This pangenome contains 3.87 Gb of sequences that are absent from the O. sativa ssp. japonica cv. Nipponbare reference genome. We captured alternate assemblies that include heterozygous information missing in the primary assemblies, and identified a total of 69,531 pan-genes, with 28,907 core genes and 13,728 wild-rice-specific genes. We observed a higher abundance and a significantly greater diversity of resistance-gene analogues in wild rice than in cultivars. Our analysis indicates that two cultivated subpopulations, intro-indica and basmati, were generated through gene flows among cultivars in South Asia. We also provide strong evidence to support the theory that the initial domestication of all Asian cultivated rice occurred only once. Furthermore, we captured 855,122 differentiated single-nucleotide polymorphisms and 13,853 differentiated presence-absence variations between indica and japonica, which could be traced to the divergence of their respective ancestors and the existence of a larger genetic bottleneck in japonica. This study provides reference resources for enhancing rice breeding, and enriches our understanding of the origins and domestication process of rice.

Revealing Genomic Traits and Evolutionary Insights of Oryza officinalis from Southern China Through Genome Assembly and Transcriptome Analysis

Wild rice, as the ancestor of cultivated rice, has accumulated a wide range of beneficial traits through prolonged natural selection and evolution. Oryza officinalis, belonging to the CC genome, differs significantly from the AA genome. In this study, we utilized second- and third-generation sequencing, along with Hi-C technology, to assemble the genome of MT10 (O. officinalis). The assembled genome is 552.58 Mb, with contigs and scaffold N50 values of 40.04 and 44.48 Mb, respectively, and 96.73% of the sequences anchored to 12 chromosomes. A total of 33,813 genes were annotated, and repetitive sequences account for 54.24% of the MT10 genome. The number of unique genes in MT10 exceeds that in the O. officinalis genome from Thailand, and their divergence time is estimated at 1.6 million years ago. The MT10 genome exhibits fewer expanded gene families compared to contracted ones, with the expanded families predominantly associated with disease and pest resistance. Comparative genomic analysis of MT10 and Nipponbare reveals sequence variations in biotic and abiotic resistance-related genes. In particular, the presence of R genes and cystatin gene families in MT10 may contribute to its unique insect resistance. Transcriptome analyses indicate that flavonoid biosynthesis and MAPK-related genes are expressed in response to brown planthopper infestation. This study represents the first chromosome-level genome assembly of MT10, providing a reference sequence for the efficient cloning of beneficial genes from O. officinalis, which holds significant potential for the genetic improvement of cultivated rice.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

D, Sk H, Solanki M, Gokulan CG, Das A, Miriyala A, Gonuguntala R, Elumalai P, M. B. V. N K, S. K. M, Kumboju C, Arra Y, G. S. L, Chirravuri NN, Patel HK, Ghazi IA, Kim SR, Jena KK, Hanumanth SR, Oliva R, Mangrauthia SK, Sundaram RM. Fine mapping and sequence analysis reveal a promising candidate gene encoding a novel NB-ARC domain derived from wild rice (Oryza officinalis) that confers bacterial blight resistance

Bacterial blight disease of rice caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most serious constraints in rice production. The most sustainable strategy to combat the disease is the deployment of host plant resistance. Earlier, we identified an introgression line, IR 75084-15-3-B-B, derived from Oryza officinalis possessing broad-spectrum resistance against Xoo. In order to understand the inheritance of resistance in the O. officinalis accession and identify genomic region(s) associated with resistance, a recombinant inbred line (RIL) mapping population was developed from the cross Samba Mahsuri (susceptible to bacterial blight) × IR 75084-15-3-B-B (resistant to bacterial blight). The F2 population derived from the cross segregated in a phenotypic ratio of 3: 1 (resistant susceptible) implying that resistance in IR 75084-15-3-B-B is controlled by a single dominant gene/quantitative trait locus (QTL). In the F7 generation, a set of 47 homozygous resistant lines and 47 homozygous susceptible lines was used to study the association between phenotypic data obtained through screening with Xoo and genotypic data obtained through analysis of 7K rice single-nucleotide polymorphism (SNP) chip. Through composite interval mapping, a major locus was detected in the midst of two flanking SNP markers, viz., Chr11.27817978 and Chr11.27994133, on chromosome 11L with a logarithm of the odds (LOD) score of 10.21 and 35.93% of phenotypic variation, and the locus has been named Xa48t. In silico search in the genomic region between the two markers flanking Xa48t identified 10 putatively expressed genes located in the region of interest. The quantitative expression and DNA sequence analysis of these genes from contrasting parents identified the Os11g0687900 encoding an NB-ARC domain-containing protein as the most promising gene associated with resistance. Interestingly, a 16-bp insertion was noticed in the untranslated region (UTR) of the gene in the resistant parent, IR 75084-15-3-B-B, which was absent in Samba Mahsuri. The association of Os11g0687900 with resistance phenotype was further established by sequence-based DNA marker analysis in the RIL population. A co-segregating PCR-based INDEL marker, Marker_Xa48, has been developed for use in the marker-assisted breeding of Xa48t.

The International Oryza Map Alignment Project (IOMAP): the Americas-past achievements and future directions

The wild relatives of rice hold unexplored genetic diversity that can be employed to feed an estimated population of 10 billion by 2050. The Oryza Map Alignment Project (OMAP) initiated in 2003 has provided comprehensive genomic resources for comparative, evolutionary, and functional characterization of the wild relatives of rice, facilitating the cloning of >600 rice genes, including those for grain width (GW5) and submergence tolerance (SUB1A). Following in the footsteps of the original project, the goal of 'IOMAP: the Americas' is to investigate the present and historic genetic diversity of wild Oryza species endemic to the Americas through the sequencing of herbaria and in situ specimens. The generation of a large diversity panel describing past and current genetic status and potential erosion of genetic variation in the populations will provide useful knowledge for the conservation of the biodiversity in these species. The wild relatives of rice in the Americas present a wide range of resistance traits useful for crop improvement and neodomestication approaches. In the race against time for a sustainable food future, the neodomestication of the first cereal species recently accomplished in O. alta opens the door to the potential neodomestication of the other wild Oryza species in Americas.

Chromosome-Level Genome Assembly of a Fragrant Japonica Rice Cultivar 'Changxianggeng 1813' Provides Insights into Genomic Variations between Fragrant and Non-Fragrant Japonica Rice

East Asia has an abundant resource of fragrant japonica rice that is gaining increasing interest among both consumers and producers. However, genomic resources and in particular complete genome sequences currently available for the breeding of fragrant japonica rice are still scarce. Here, integrating Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods, we presented a high-quality chromosome-level genome assembly (~378.78 Mb) for a new fragrant japonica cultivar ‘Changxianggeng 1813’, with 31,671 predicated protein-coding genes. Based on the annotated genome sequence, we demonstrated that it was the badh2-E2 type of deletion (a 7-bp deletion in the second exon) that caused fragrance in ‘Changxianggeng 1813’. Comparative genomic analyses revealed that multiple gene families involved in the abiotic stress response were expanded in the ‘Changxianggeng 1813’ genome, which further supported the previous finding that no generalized loss of abiotic stress tolerance associated with the fragrance phenotype. Although the ‘Changxianggeng 1813’ genome showed high genomic synteny with the genome of the non-fragrant japonica rice cultivar Nipponbare, a total of 289,970 single nucleotide polymorphisms (SNPs), 96,093 small insertion-deletion polymorphisms (InDels), and 8690 large structure variants (SVs, >1000 bp) were identified between them. Together, these genomic resources will be valuable for elucidating the mechanisms underlying economically important traits and have wide-ranging implications for genomics-assisted breeding in fragrant japonica rice.

The Original Form of C4-Photosynthetic Phosphoenolpyruvate Carboxylase Is Retained in Pooids but Lost in Rice

Poaceae is the most prominent monocot family that contains the primary cereal crops wheat, rice, and maize. These cereal species exhibit physiological diversity, such as different photosynthetic systems and environmental stress tolerance. Phosphoenolpyruvate carboxylase (PEPC) in Poaceae is encoded by a small multigene family and plays a central role in C4-photosynthesis and dicarboxylic acid metabolism. Here, to better understand the molecular basis of the cereal species diversity, we analyzed the PEPC gene family in wheat together with other grass species. We could designate seven plant-type and one bacterial-type grass PEPC groups, ppc1a, ppc1b, ppc2a, ppc2b, ppc3, ppc4, ppcC4, and ppc-b, respectively, among which ppc1b is an uncharacterized type of PEPC. Evolutionary inference revealed that these PEPCs were derived from five types of ancient PEPCs (ppc1, ppc2, ppc3, ppc4, and ppc-b) in three chromosomal blocks of the ancestral Poaceae genome. C4-photosynthetic PEPC (ppcC4 ) had evolved from ppc1b, which seemed to be arisen by a chromosomal duplication event. We observed that ppc1b was lost in many Oryza species but preserved in Pooideae after natural selection. In silico analysis of cereal RNA-Seq data highlighted the preferential expression of ppc1b in upper ground organs, selective up-regulation of ppc1b under osmotic stress conditions, and nitrogen response of ppc1b. Characterization of wheat ppc1b showed high levels of gene expression in young leaves, transcriptional responses under nitrogen and abiotic stress, and the presence of a Dof1 binding site, similar to ppcC4 in maize. Our results indicate the evolving status of Poaceae PEPCs and suggest the functional association of ppc1-derivatives with adaptation to environmental changes.

The first long-read nuclear genome assembly of Oryza australiensis, a wild rice from northern Australia

Oryza australiensis is a wild rice native to monsoonal northern Australia. The International Oryza Map Alignment Project emphasises its significance as the sole representative of the EE genome clade. Assembly of the O. australiensis genome has previously been challenging due to its high Long Terminal Repeat (LTR) retrotransposon (RT) content. Oxford Nanopore long reads were combined with Illumina short reads to generate a high-quality ~ 858 Mbp genome assembly within 850 contigs with 46× long read coverage. Reference-guided scaffolding increased genome contiguity, placing 88.2% of contigs into 12 pseudomolecules. After alignment to the Oryza sativa cv. Nipponbare genome, we observed several structural variations. PacBio Iso-Seq data were generated for five distinct tissues to improve the functional annotation of 34,587 protein-coding genes and 42,329 transcripts. We also report SNV numbers for three additional O. australiensis genotypes based on Illumina re-sequencing. Although genetic similarity reflected geographical separation, the density of SNVs also correlated with our previous report on variations in salinity tolerance. This genome re-confirms the genetic remoteness of the O. australiensis lineage within the O. officinalis genome complex. Assembly of a high-quality genome for O. australiensis provides an important resource for the discovery of critical genes involved in development and stress tolerance.

Unveiling a Novel Source of Resistance to Bacterial Blight in Medicinal Wild Rice, Oryza officinalis

Bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae (Xoo) is among the oldest known bacterial diseases found for rice in Asia. It is the most serious bacterial disease in many rice growing regions of the world. A total of 47 resistance (R) genes (Xa1 to Xa47) have been identified. Nonetheless, these R genes could possibly be defeated to lose their qualitative nature and express intermediate phenotypes. The identification of sources of novel genetic loci regulating host plant resistance is crucial to develop an efficient control strategy. Wild ancestors of cultivated rice are a natural genetic resource contain a large number of excellent genes. Medicinal wild rice (Oryza officinalis) belongs to the CC genome and is a well-known wild rice in south China. In this study, O. officinalis was crossed with cultivated rice HY-8 and their hybrids were screened for BB resistance genes deployed through natural selection in wild rice germplasm. The molecular markers linked to R genes for BB were used to screen the genomic regions in wild parents and their recombinants. The gene coding and promoter regions of major R genes were inconsistently found in O. officinalis and its progenies. Oryza officinalis showed resistance to all thirty inoculated Xoo strains with non-availability of various known R genes. The results indicated the presence of novel genomic regions for BB resistance in O. officinalis. The present study not only provides a reference to investigate medicinal rice for R gene(s) identification against BB but also identified it as a new breeding material for BB resistance.

Water stress resilient cereal crops: Lessons from wild relatives

Cereal crops are significant contributors to global diets. As climate change disrupts weather patterns and wreaks havoc on crops, the need for generating stress-resilient, high-yielding varieties is more urgent than ever. One extremely promising avenue in this regard is to exploit the tremendous genetic diversity expressed by the wild ancestors of current day crop species. These crop wild relatives thrive in a range of environments and accordingly often harbor an array of traits that allow them to do so. The identification and introgression of these traits into our staple cereal crops can lessen yield losses in stressful environments. In the last decades, a surge in extreme drought and flooding events have severely impacted cereal crop production. Climate models predict a persistence of this trend, thus reinforcing the need for research on water stress resilience. Here we review: (i) how water stress (drought and flooding) impacts crop performance; and (ii) how identification of tolerance traits and mechanisms from wild relatives of the main cereal crops, that is, rice, maize, wheat, and barley, can lead to improved survival and sustained yields in these crops under water stress conditions.

Proto Kranz-like leaf traits and cellular ionic regulation are associated with salinity tolerance in a halophytic wild rice

Species of wild rice (Oryza spp.) possess a wide range of stress tolerance traits that can be potentially utilized in breeding climate-resilient cultivated rice cultivars (Oryza sativa) thereby aiding global food security. In this study, we conducted a greenhouse trial to evaluate the salinity tolerance of six wild rice species, one cultivated rice cultivar (IR64) and one landrace (Pokkali) using a range of electrophysiological, imaging, and whole-plant physiological techniques. Three wild species (O. latifolia, O. officinalis and O. coarctata) were found to possess superior salinity stress tolerance. The underlying mechanisms, however, were strikingly different. Na+ accumulation in leaves of O. latifolia, O. officinalis and O. coarctata were significantly higher than the tolerant landrace, Pokkali. Na+ accumulation in mesophyll cells was only observed in O. coarctata, suggesting that O. officinalis and O. latifolia avoid Na+ accumulation in mesophyll by allocating Na+ to other parts of the leaf. The finding also suggests that O. coarctata might be able to employ Na+ as osmolyte without affecting its growth. Further study of Na+ allocation in leaves will be helpful to understand the mechanisms of Na+ accumulation in these species. In addition, O. coarctata showed Proto Kranz-like leaf anatomy (enlarged bundle sheath cells and lower numbers of mesophyll cells), and higher expression of C4-related genes (e.g., NADPME, PPDK) and was a clear outlier with respect to salinity tolerance among the studied wild and cultivated Oryza species. The unique phylogenetic relationship of O. coarctata with C4 grasses suggests the potential of this species for breeding rice with high photosynthetic rate under salinity stress in the future.

Genome Size Variation and Evolution Driven by Transposable Elements in the Genus Oryza

Genome size variation and evolutionary forces behind have been long pursued in flowering plants. The genus Oryza, consisting of approximately 25 wild species and two cultivated rice, harbors eleven extant genome types, six of which are diploid (AA, BB, CC, EE, FF, and GG) and five of which are tetraploid (BBCC, CCDD, HHJJ, HHKK, and KKLL). To obtain the most comprehensive knowledge of genome size variation in the genus Oryza, we performed flow cytometry experiments and estimated genome sizes of 166 accessions belonging to 16 non-AA genome Oryza species. k-mer analyses were followed to verify the experimental results of the two accessions for each species. Our results showed that genome sizes largely varied fourfold in the genus Oryza, ranging from 279 Mb in Oryza brachyantha (FF) to 1,203 Mb in Oryza ridleyi (HHJJ). There was a 2-fold variation (ranging from 570 to 1,203 Mb) in genome size among the tetraploid species, while the diploid species had 3-fold variation, ranging from 279 Mb in Oryza brachyantha (FF) to 905 Mb in Oryza australiensis (EE). The genome sizes of the tetraploid species were not always two times larger than those of the diploid species, and some diploid species even had larger genome sizes than those of tetraploids. Nevertheless, we found that genome sizes of newly formed allotetraploids (BBCC-) were almost equal to totaling genome sizes of their parental progenitors. Our results showed that the species belonging to the same genome types had similar genome sizes, while genome sizes exhibited a gradually decreased trend during the evolutionary process in the clade with AA, BB, CC, and EE genome types. Comparative genomic analyses further showed that the species with different rice genome types may had experienced dissimilar amplification histories of retrotransposons, resulting in remarkably different genome sizes. On the other hand, the closely related rice species may have experienced similar amplification history. We observed that the contents of transposable elements, long terminal repeats (LTR) retrotransposons, and particularly LTR/Gypsy retrotransposons varied largely but were significantly correlated with genome sizes. Therefore, this study demonstrated that LTR retrotransposons act as an active driver of genome size variation in the genus Oryza.

Collection, preservation and distribution of Oryza genetic resources by the National Bioresource Project RICE (NBRP-RICE)

Biological resources are the basic infrastructure of bioscience research. Rice (Oryza sativa L.) is a good experimental model for research in cereal crops and monocots and includes important genetic materials used in breeding. The availability of genetic materials, including mutants, is important for rice research. In addition, Oryza species are attractive to researchers for both finding useful genes for breeding and for understanding the mechanism of genome evolution that enables wild plants to adapt to their own habitats. NBRP-RICE contributes to rice research by promoting the usage of genetic materials, especially wild Oryza accessions and mutant lines. Our activity includes collection, preservation and distribution of those materials and the provision of basic information on them, such as morphological and physiological traits and genomic information. In this review paper, we introduce the activities of NBRP-RICE and our database, Oryzabase, which facilitates the access to NBRP-RICE resources and their genomic sequences as well as the current situation of wild Oryza genome sequencing efforts by NBRP-RICE and other institutes.

OryzaGenome2.1: Database of Diverse Genotypes in Wild Oryza Species

BACKGROUND

OryzaGenome ( http://viewer.shigen.info/oryzagenome21detail/index.xhtml ), a feature within Oryzabase ( https://shigen.nig.ac.jp/rice/oryzabase/ ), is a genomic database for wild Oryza species that provides comparative and evolutionary genomics approaches for the rice research community.

RESULTS

Here we release OryzaGenome2.1, the first major update of OryzaGenome. The main feature in this version is the inclusion of newly sequenced genotypes and their meta-information, giving a total of 217 accessions of 19 wild Oryza species (O. rufipogon, O. barthii, O. longistaminata, O. meridionalis, O. glumaepatula, O. punctata, O. minuta, O. officinalis, O. rhizomatis, O. eichingeri, O. latifolia, O. alta, O. grandiglumis, O. australiensis, O. brachyantha, O. granulata, O. meyeriana, O. ridleyi, and O. longiglumis). These 19 wild species belong to 9 genome types (AA, BB, CC, BBCC, CCDD, EE, FF, GG, and HHJJ), representing wide genomic diversity in the genus. Using the genotype information, we analyzed the genome diversity of Oryza species. Other features of OryzaGenome facilitate the use of information on single nucleotide polymorphisms (SNPs) between O. sativa and its wild progenitor O. rufipogon in rice research, including breeding as well as basic science. For example, we provide Variant Call Format (VCF) files for genome-wide SNPs of 33 O. rufipogon accessions against the O. sativa reference genome, IRGSP1.0. In addition, we provide a new SNP Effect Table function, allowing users to identify SNPs or small insertion/deletion polymorphisms in the 33 O. rufipogon accessions and to search for the effect of these polymorphisms on protein function if they reside in the coding region (e.g., are missense or nonsense mutations). Furthermore, the SNP Viewer for 446 O. rufipogon accessions was updated by implementing new tracks for possible selective sweep regions and highly mutated regions that were potentially exposed to selective pressures during the process of domestication.

CONCLUSION

OryzaGenome2.1 focuses on comparative genomic analysis of diverse wild Oryza accessions collected around the world and on the development of resources to speed up the identification of critical trait-related genes, especially from O. rufipogon. It aims to promote the use of genotype information from wild accessions in rice breeding and potential future crop improvements. Diverse genotypes will be a key resource for evolutionary studies in Oryza, including polyploid biology.

Transcriptome analysis of diploid and triploid Populus tomentosa

Triploid Chinese white poplar (Populus tomentosa Carr., Salicaceae) has stronger advantages in growth and better stress resistance and wood quality than diploid P. tomentosa. Using transcriptome sequencing technology to identify candidate transcriptome-based markers for growth vigor in young tree tissue is of great significance for the breeding of P. tomentosa varieties in the future. In this study, the cuttings of diploid and triploid P. tomentosa were used as plant materials, transcriptome sequencing was carried out, and their tissue culture materials were used for RT-qPCR verification of the expression of genes. The results showed that 12,240 differentially expressed genes in diploid and triploid P. tomentosa transcripts were annotated and enriched into 135 metabolic pathways. The top six pathways that enriched the most significantly different genes were plant-pathogen interaction, phenylpropanoid biosynthesis, MAPK signalling pathway-plant, ascorbate and aldarate metabolism, diterpenoid biosynthesis, and the betalain biosynthesis pathway. Ten growth-related genes were selected from pathways of plant hormone signal transduction and carbon fixation in photosynthetic organisms for RT-qPCR verification. The expression levels of MDH and CYCD3 in tissue-cultured and greenhouse planted triploid P. tomentosa were higher than those in tissue-cultured diploid P. tomentosa, which was consist ent with the TMM values calculated by transcriptome.

Agrobacterium-Mediated Genetic Transformation of Wild Oryza Species Using Immature Embryos

Genetic transformation is one of the most important technologies for revealing or modulating gene function. It is used widely in both functional genomics and molecular breeding of rice. Demands on its use in wild Oryza species is increasing because of their high genetic diversity. Given the difficulties in genetic crosses between distantly related species, genetic transformation offers a way to alter or transfer genetic traits in wild rice accessions. However, transformation of wild Oryza accessions by conventional methods using calli induced from scutellum tissue of embryos in mature seeds often fails. Here, we report methods using immature embryos for the genetic transformation of a broad range of Oryza species. First, we investigated the ability of callus induction and regeneration from immature embryos of 192 accessions in 20 species under several culture conditions. We regenerated plants from immature embryos of 90 accessions in 16 species. Next, we optimized the conditions of Agrobacterium infection using a vector carrying the GFP gene driven by the maize ubiquitin promoter. GFP signals were observed in 51 accessions in 11 species. We analyzed the growth and seed set of transgenic plants of O. barthii, O. glumaepatula, O. rufipogon, and O. brachyantha. The plants grew to maturity and set seeds normally. Southern blot analyses using DNA from T0 plants showed that all GFP plants were derived from independent transformation events. We confirmed that the T-DNAs were transmitted to the next generation through the segregation of GFP signals in the T1 generation. These results show that many Oryza species can be transformed by using modified immature-embryo methods. This will accelerate the use of wild Oryza accessions in molecular genetic analyses and molecular breeding.