Loss of OsEAF6, a Subunit of the Histone Acetyltransferase Complex, Causes Hybrid Breakdown in Intersubspecific Rice Crosses.

Gene duplication plays an important role in genetic diversification, adaptive evolution, and speciation. Understanding the mechanisms and effects of postzygotic isolation genes is important for further studies of speciation and crop breeding. The duplicate recessive genes hwe1 and hwe2 cause hybrid breakdown, characterized by poor vegetative growth and reproductive dysgenesis in intersubspecific crosses between Oryza sativa ssp. indica and japonica. Using a map-based cloning strategy, we found that HWE1 and HWE2 encode the Esa1-associated factor 6 (EAF6) protein, a component of histone acetyltransferase complexes. The indica hwe1 and japonica hwe2 alleles lacked functional EAF6, demonstrating that the double recessive homozygote causes hybrid breakdown. Morphological and physiological observations showed that weak plants with double recessive homozygotes had serious morphological defects with a wide range of effects on development and organs, leading to leaves with reduced chlorophyll content, flower and pistil malformation, and anomalies of gametogenesis. These findings suggest that EAF6 plays a pivotal role in the transcriptional regulation of essential genes during the vegetative and reproductive development of rice.

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.

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

The Oryza officinalis complex is the largest species group in Oryza, with more than nine species from four continents, and is a tertiary gene pool that can be exploited in breeding programs for the improvement of cultivated rice. Most diploid and tetraploid members of this group have a C genome. Using a new reference C genome for the diploid species O. officinalis, and draft genomes for two other C genome diploid species Oryza eichingeri and Oryza rhizomatis, we examine the influence of transposable elements on genome structure and provide a detailed phylogeny and evolutionary history of the Oryza C genomes. The O. officinalis genome is 1.6 times larger than the A genome of cultivated Oryza sativa, mostly due to proliferation of Gypsy type long-terminal repeat transposable elements, but overall syntenic relationships are maintained with other Oryza genomes (A, B, and F). Draft genome assemblies of the two other C genome diploid species, Oryza eichingeri and Oryza rhizomatis, and short-read resequencing of a series of other C genome species and accessions reveal that after the divergence of the C genome progenitor, there was still a substantial degree of variation within the C genome species through proliferation and loss of both DNA and long-terminal repeat transposable elements. We provide a detailed phylogeny and evolutionary history of the Oryza C genomes and a genomic resource for the exploitation of the Oryza tertiary gene pool.

Potential of Oryza officinalis to augment the cold tolerance genetic mechanisms of Oryza sativa by network complementation.

Oryza officinalis is an accessible alien donor for genetic improvement of rice. Comparison across a representative panel of Oryza species showed that the wild O. officinalis and cultivated O. sativa ssp. japonica have similar cold tolerance potentials. The possibility that either distinct or similar genetic mechanisms are involved in the low temperature responses of each species was addressed by comparing their transcriptional networks. General similarities were supported by shared transcriptomic signatures indicative of equivalent metabolic, hormonal, and defense status. However, O. officinalis has maintained an elaborate cold-responsive brassinosteroid-regulated BES1-network that appeared to have been fragmented in O. sativa. BES1-network is potentially important for integrating growth-related responses with physiological adjustments and defenses through the protection of photosynthetic machinery and maintenance of stomatal aperture, oxidative defenses, and osmotic adjustment. Equivalent physiological processes are functional in O. sativa but their genetic mechanisms are under the direct control of ABA-dependent, DREB-dependent and/or oxidative-mediated networks uncoupled to BES1. While O. officinalis and O. sativa represent long periods of speciation and domestication, their comparable cold tolerance potentials involve equivalent physiological processes but distinct genetic networks. BES1-network represents a novel attribute of O. officinalis with potential applications in diversifying or complementing other mechanisms in the cultivated germplasm.

Assembling the genome of the African wild rice Oryza longistaminata by exploiting synteny in closely related Oryza species.

The African wild rice species Oryza longistaminata has several beneficial traits compared to cultivated rice species, such as resistance to biotic stresses, clonal propagation via rhizomes, and increased biomass production. To facilitate breeding efforts and functional genomics studies, we de-novo assembled a high-quality, haploid-phased genome. Here, we present our assembly, with a total length of 351 Mb, of which 92.2% was anchored onto 12 chromosomes. We detected 34,389 genes and 38.1% of the genome consisted of repetitive content. We validated our assembly by a comparative linkage analysis and by examining well-characterized gene families. This genome assembly will be a useful resource to exploit beneficial alleles found in O. longistaminata. Our results also show that it is possible to generate a high-quality, functionally complete rice genome assembly from moderate SMRT read coverage by exploiting synteny in a closely related Oryza species.

Characterization of diterpene synthase genes in the wild rice species Oryza brachyatha provides evolutionary insight into rice phytoalexin biosynthesis.

Cultivated rice (Oryza sativa; Os) produces a variety of labdane-related diterpenoids; not only phytohormone gibberellins (GAs) but also phytoalexins for defense including phytocassanes, momilactones and oryzalexins. Their carbon skeleton diterpenes are constructed from geranylgeranyl diphosphate via ent-copalyl diphosphate (ent-CDP) or its diastereomer syn-CDP. These two-step reactions are successively catalyzed by homologs of the two diterpene synthases, ent-CDP synthase (ent-CPS) and ent-kaurene synthase (KS) that are responsible for the biosynthesis of GAs; e.g. OsCPS4 and OsKSL8 that are involved in the biosynthesis of oryzalexin S, a rice phytoalexin. Oryza brachyantha (Ob) is the most distant wild rice species from Os among the Oryza genus. We previously reported that the Ob genome contains ObCPS_11g, ObKSL8-a, ObKSL8-b and ObKSL8-c for specialized metabolism at a locus similar to the OsKSL8 locus on chromosome 11. These Ob genes are closely related to OsCPS4 and OsKSL8, respectively. We herein characterize the diterpene synthase genes in Ob, using functional analyses and expression analysis. Recombinant OsKSL8 and ObKSL8-a showed the same in vitro function when syn-CDP or normal-CDP were used as substrates. Nonetheless, our results suggest that Ob produces normal-CDP-related diterpenoid phytoalexins, presumably via ObKSL8-a, while Os produces a syn-CDP-related phytoalexin, oryzalexin S, via OsKSL8. This difference must be due to the kinds of CPS that are present in each species; Os has OsCPS4 encoding syn-CPS, while Ob has ObCPS_11g encoding normal-CPS. Thus, we propose the evolutionary history underlying oryzalexin S biosynthesis: the gain of a syn-CPS was a critical event allowing the biosynthesis of oryzalexin S.

Time-Course Transcriptomics Analysis Reveals Key Responses of Submerged Deepwater Rice to Flooding.

Water submergence is an environmental factor that limits plant growth and survival. Deepwater rice (Oryza sativa) adapts to submergence by rapidly elongating its internodes and thereby maintaining its leaves above the water surface. We performed a comparative RNA sequencing transcriptome analysis of the shoot base region, including basal nodes, internodes, and shoot apices of seedlings at two developmental stages from two varieties with contrasting deepwater growth responses. A transcriptomic comparison between deepwater rice cv C9285 and nondeepwater rice cv Taichung 65 revealed both similar and differential expression patterns between the two genotypes during submergence. The expression of genes related to gibberellin biosynthesis, trehalose biosynthesis, anaerobic fermentation, cell wall modification, and transcription factors that include ethylene-responsive factors was significantly different between the varieties. Interestingly, in both varieties, the jasmonic acid content at the shoot base decreased during submergence, while exogenous jasmonic acid inhibited submergence-induced internode elongation in cv C9285, suggesting that jasmonic acid plays a role in the submergence response of rice. Furthermore, a targeted de novo transcript assembly revealed transcripts that were specific to cv C9285, including submergence-induced biotic stress-related genes. Our multifaceted transcriptome approach using the rice shoot base region illustrates a differential response to submergence between deepwater and nondeepwater rice. Jasmonic acid metabolism appears to participate in the submergence-mediated internode elongation response of deepwater rice.

Cytokinin-induced expression of OSH1 in a shoot-regenerating rice callus.

The expression of a KNOX class 1 gene OSH1 is induced by cytokinin during regeneration of shoots from callus in Oryza sativa L. (rice). This cytokinin-induced expression was enhanced by overexpression of homologues of cytokinin-signalling phosphorelay genes such as a histidine kinase gene OHK3, a phosphotransmitter gene OHP2 and a response regulator gene ORR1 in cultured cells. Regionally overlapped expression of these genes and OSH1 was observed in shoot apex. These results suggest that these cytokinin-signalling genes are positive regulators of the expression of OSH1, and mediate the OSH expression upon shoot regeneration from callus in rice.

Overcoming the species hybridization barrier by ploidy manipulation in the genus Oryza.

In most eudicot and monocot species, interspecific and interploidy crosses generally display abnormalities in the endosperm that are the major cause of a post-zygotic hybridization barrier. In some eudicot species, however, this type of hybridization barrier can be overcome by the manipulation of ploidy levels of one parental species, suggesting that the molecular mechanisms underlying the species hybridization barrier can be circumvented by genome dosage. We previously demonstrated that endosperm barriers in interspecific and interploidy crosses in the genus Oryza involve overlapping but different mechanisms. This result contrasts with those in the genus Arabidopsis, which shows similar outcomes in both interploidy and interspecific crosses. Therefore, we postulated that an exploration of pathways for overcoming the species hybridization barrier in Oryza endosperm, by manipulating the ploidy levels in one parental species, might provide novel insights into molecular mechanisms. We showed that fertile hybrid seeds could be produced by an interspecific cross of female tetraploid Oryza sativa and male diploid Oryza longistaminata. Although the rate of nuclear divisions did not return to normal levels in the hybrid endosperm, the timing of cellularization, nucellus degeneration and the accumulation of storage products were close to normal levels. In addition, the expression patterns of the imprinted gene MADS87 and YUCCA11 were changed when the species barrier was overcome. These results suggest that the regulatory machinery for developmental transitions and imprinted gene expression are likely to play a central role in overcoming species hybridization barriers by genome dosage in the genus Oryza.

Genetic characterization and fine mapping of S25, a hybrid male sterility gene, on rice chromosome 12.

Hybrid male sterility genes are important factors in creating postzygotic reproductive isolation barriers in plants. One such gene, S25, is known to cause severe transmission ratio distortion in inter-subspecific progeny of cultivated rice Oryza sativa ssp. indica and japonica. To further characterize the S25 gene, we fine-mapped and genetically characterized the S25 gene using near-isogenic lines with reciprocal genetic backgrounds. We mapped the S25 locus within the 0.67-1.02 Mb region on rice chromosome 12. Further genetic analyses revealed that S25 substantially reduced male fertility in the japonica background, but not in the indica background. In first-generation hybrid progeny, S25 had a milder effect than it had in the japonica background. These results suggest that the expression of S25 is epistatically regulated by at least one partially dominant gene present in the indica genome. This finding supports our previous studies showing that hybrid male sterility due to pollen killer genes results from epistatic interaction with other genes that are hidden in the genetic background.

DNApod: DNA polymorphism annotation database from next-generation sequence read archives.

With the rapid advances in next-generation sequencing (NGS), datasets for DNA polymorphisms among various species and strains have been produced, stored, and distributed. However, reliability varies among these datasets because the experimental and analytical conditions used differ among assays. Furthermore, such datasets have been frequently distributed from the websites of individual sequencing projects. It is desirable to integrate DNA polymorphism data into one database featuring uniform quality control that is distributed from a single platform at a single place. DNA polymorphism annotation database (DNApod; http://tga.nig.ac.jp/dnapod/) is an integrated database that stores genome-wide DNA polymorphism datasets acquired under uniform analytical conditions, and this includes uniformity in the quality of the raw data, the reference genome version, and evaluation algorithms. DNApod genotypic data are re-analyzed whole-genome shotgun datasets extracted from sequence read archives, and DNApod distributes genome-wide DNA polymorphism datasets and known-gene annotations for each DNA polymorphism. This new database was developed for storing genome-wide DNA polymorphism datasets of plants, with crops being the first priority. Here, we describe our analyzed data for 679, 404, and 66 strains of rice, maize, and sorghum, respectively. The analytical methods are available as a DNApod workflow in an NGS annotation system of the DNA Data Bank of Japan and a virtual machine image. Furthermore, DNApod provides tables of links of identifiers between DNApod genotypic data and public phenotypic data. To advance the sharing of organism knowledge, DNApod offers basic and ubiquitous functions for multiple alignment and phylogenetic tree construction by using orthologous gene information.

Characterization and evolutionary analysis of ent-kaurene synthase like genes from the wild rice species Oryza rufipogon.

Cultivated rice (Oryza sativa) possesses various labdane-related diterpene synthase genes, homologs of ent-copalyl diphosphate synthase (CPS) and ent-kaurene synthase (KS) that are responsible for the biosynthesis of phytohormone gibberellins. The CPS homologs and KS like (KSL) homologs successively converted geranylgeranyl diphosphate to cyclic diterpene hydrocarbons via ent-copalyl diphosphate or syn-copalyl diphosphate in O. sativa. Consequently, a variety of labdane-related diterpenoids, including phytoalexin phytocassanes, momilactones and oryzalexins, have been identified from cultivated rice. Our previous report indicated that the biosynthesis of phytocassanes and momilactones is conserved in Oryza rufipogon, the progenitor of Asian cultivated rice. Moreover, their biosynthetic gene clusters, containing OsCPS2 and OsKSL7 for phytocassane biosynthesis and OsCPS4 and OsKSL4 for momilactone biosynthesis, are also present in the O. rufipogon genome. We herein characterized O. rufipogon homologs of OsKSL5, OsKSL6, OsKSL8 responsible for oryzalexin S biosynthesis, and OsKSL10 responsible for oryzalexins A-F biosynthesis, to obtain more evolutionary insight into diterpenoid biosynthesis in O. sativa. Our phytoalexin analyses showed that no accumulation of oryzalexins was detected in extracts from O. rufipogon leaf blades. In vitro functional analyses indicated that unlike OsKSL10, O. rufipogon KSL10 functions as an ent-miltiradiene synthase, which explains the lack of accumulation of oryzalexins A-F in O. rufipogon. The different functions of KSL5 and KSL8 in O. sativa japonica to those in indica are conserved in each type of O. rufipogon, while KSL6 functions (ent-isokaurene synthases) are well conserved. Our study suggests that O. sativa japonica has evolved distinct specialized diterpenoid metabolism, including the biosynthesis of oryzalexins.

Effect of Wild and Cultivated Rice Genotypes on Rhizosphere Bacterial Community Composition.

BACKGROUND: Deposition and secretion from roots influences the composition of the microbial communities surrounding them in the rhizosphere, and microbial activities influence the growth and health of the plant. Different host plant genotypes result in differences in those microbial communities. Crop genomes may have a narrow genetic base because of bottlenecks that occurred when domesticated crops were derived from small populations within the progenitor species. Desirable traits influencing root-associated microbial communities might therefore have been lost in the transition from wild species to modern cultivars. To investigate the diversity of bacterial communities associated with wild and cultivated rice, we surveyed a series of plant species and cultivars spanning the Oryza genus, growing them in the same nutrient-poor soil and assessing the bacterial composition of their rhizospheres and the surrounding soil using 16S rDNA sequencing. RESULTS: Root-associated bacterial communities showed small but significant differences dependent on the plant genotype. We found that differences between bacteria associated with differing plant genotypes were only weakly correlated with the phylogenetic distance between the Oryza wild species and cultivars. In ordination plots, domesticated and wild samples could be separated on the basis of their associated bacterial communities. Taxa of the Anaerolineae were overrepresented in wild samples compared to domesticated ones. Certain methanotrophs were overrepresented in the earliest diverged part of the Oryza genus. CONCLUSIONS: Bacterial populations associated with the rhizosphere of wild rice species displayed differences with those associated with cultivars, suggesting that root traits selected in domestication could have significant influence on the rhizosphere microbiota composition. Variation within the genus seems to influence the representation of methanotrophs. This suggests that greenhouse emissions from paddy fields could be altered by manipulating plant genotypes through the introgression of wild rice genetic material.

Evolutionary trajectory of phytoalexin biosynthetic gene clusters in rice.

Plants frequently possess operon-like gene clusters for specialized metabolism. Cultivated rice, Oryza sativa, produces antimicrobial diterpene phytoalexins represented by phytocassanes and momilactones, and the majority of their biosynthetic genes are clustered on chromosomes 2 and 4, respectively. These labdane-related diterpene phytoalexins are biosynthesized from geranylgeranyl diphosphate via ent-copalyl diphosphate or syn-copalyl diphosphate. The two gene clusters consist of genes encoding diterpene synthases and chemical-modification enzymes including P450s. In contrast, genes for the biosynthesis of gibberellins, which are labdane-related phytohormones, are scattered throughout the rice genome similar to other plant genomes. The mechanism of operon-like gene cluster formation remains undefined despite previous studies in other plant species. Here we show an evolutionary insight into the rice gene clusters by a comparison with wild Oryza species. Comparative genomics and biochemical studies using wild rice species from the AA genome lineage, including Oryza barthii, Oryza glumaepatula, Oryza meridionalis and the progenitor of Asian cultivated rice Oryza rufipogon indicate that gene clustering for biosynthesis of momilactones and phytocassanes had already been accomplished before the domestication of rice. Similar studies using the species Oryza punctata from the BB genome lineage, the distant FF genome lineage species Oryza brachyantha and an outgroup species Leersia perrieri suggest that the phytocassane biosynthetic gene cluster was present in the common ancestor of the Oryza species despite the different locations, directions and numbers of their member genes. However, the momilactone biosynthetic gene cluster evolved within Oryza before the divergence of the BB genome via assembly of ancestral genes.

Pollen Killer Gene S35 Function Requires Interaction with an Activator That Maps Close to S24, Another Pollen Killer Gene in Rice.

Pollen killer genes disable noncarrier pollens, and are responsible for male sterility and segregation distortion in hybrid populations of distantly related plant species. The genetic networks and the molecular mechanisms underlying the pollen killer system remain largely unknown. Two pollen killer genes, S24 and S35, have been found in an intersubspecific cross of Oryza sativa ssp. indica and japonica The effect of S24 is counteracted by an unlinked locus EFS Additionally, S35 has been proposed to interact with S24 to induce pollen sterility. These genetic interactions are suggestive of a single S24-centric genetic pathway (EFS-S24-S35) for the pollen killer system. To examine this hypothetical genetic pathway, the S35 and the S24 regions were further characterized and genetically dissected in this study. Our results indicated that S35 causes pollen sterility independently of both the EFS and S24 genes, but is dependent on a novel gene close to the S24 locus, named incentive for killing pollen (INK). We confirmed the phenotypic effect of the INK gene separately from the S24 gene, and identified the INK locus within an interval of less than 0.6 Mb on rice chromosome 5. This study characterized the genetic effect of the two independent genetic pathways of INK-S35 and EFS-S24 in indica-japonica hybrid progeny. Our results provide clear evidence that hybrid male sterility in rice is caused by several pollen killer networks with multiple factors positively and negatively regulating pollen killer genes.

OryzaGenome: Genome Diversity Database of Wild Oryza Species.

The species in the genus Oryza, encompassing nine genome types and 23 species, are a rich genetic resource and may have applications in deeper genomic analyses aiming to understand the evolution of plant genomes. With the advancement of next-generation sequencing (NGS) technology, a flood of Oryza species reference genomes and genomic variation information has become available in recent years. This genomic information, combined with the comprehensive phenotypic information that we are accumulating in our Oryzabase, can serve as an excellent genotype-phenotype association resource for analyzing rice functional and structural evolution, and the associated diversity of the Oryza genus. Here we integrate our previous and future phenotypic/habitat information and newly determined genotype information into a united repository, named OryzaGenome, providing the variant information with hyperlinks to Oryzabase. The current version of OryzaGenome includes genotype information of 446 O. rufipogon accessions derived by imputation and of 17 accessions derived by imputation-free deep sequencing. Two variant viewers are implemented: SNP Viewer as a conventional genome browser interface and Variant Table as a text-based browser for precise inspection of each variant one by one. Portable VCF (variant call format) file or tab-delimited file download is also available. Following these SNP (single nucleotide polymorphism) data, reference pseudomolecules/scaffolds/contigs and genome-wide variation information for almost all of the closely and distantly related wild Oryza species from the NIG Wild Rice Collection will be available in future releases. All of the resources can be accessed through http://viewer.shigen.info/oryzagenome/.

Two Tightly Linked Genes at the hsa1 Locus Cause Both F1 and F2 Hybrid Sterility in Rice.

Molecular mechanisms of hybrid breakdown associated with sterility (F2 sterility) are poorly understood as compared with those of F1 hybrid sterility. Previously, we characterized three unlinked epistatic loci, hybrid sterility-a1 (hsa1), hsa2, and hsa3, responsible for the F2 sterility in a cross between Oryza sativa ssp. indica and japonica. In this study, we identified that the hsa1 locus contains two interacting genes, HSA1a and HSA1b, within a 30-kb region. HSA1a-j (japonica allele) encodes a highly conserved plant-specific domain of unknown function protein (DUF1618), whereas the indica allele (HSA1a-i(s)) has two deletion mutations that cause disruption of domain structure. The second gene, HSA1b-i(s), encodes an uncharacterized protein with some similarity to a nucleotide-binding protein. Homozygous introgression of indica HSA1a-i(s)-HSA1b-i(s) alleles into japonica showed female gamete abortion at an early mitotic stage. The fact that the recombinant haplotype HSA1a-j-HSA1b-i(s) caused semi-sterility in the heterozygous state with the HSA1a-i(s)-HSA1b-i(s) haplotype suggests that variation in the hsa1 locus is a possible cause of the wide-spectrum sterility barriers seen in F1 hybrids and successive generations in rice. We propose a simple genetic model to explain how a single causal mechanism can drive both F1 and F2 hybrid sterility.

Global expression differences and tissue specific expression differences in rice evolution result in two contrasting types of differentially expressed genes.

BACKGROUND: Since the development of transcriptome analysis systems, many expression evolution studies characterized evolutionary forces acting on gene expression, without explicit discrimination between global expression differences and tissue specific expression differences. However, different types of gene expression alteration should have different effects on an organism, the evolutionary forces that act on them might be different, and different types of genes might show different types of differential expression between species. To confirm this, we studied differentially expressed (DE) genes among closely related groups that have extensive gene expression atlases, and clarified characteristics of different types of DE genes including the identification of regulating loci for differential expression using expression quantitative loci (eQTL) analysis data. RESULTS: We detected differentially expressed (DE) genes between rice subspecies in five homologous tissues that were verified using japonica and indica transcriptome atlases in public databases. Using the transcriptome atlases, we classified DE genes into two types, global DE genes and changed-tissues DE genes. Global type DE genes were not expressed in any tissues in the atlas of one subspecies, however changed-tissues type DE genes were expressed in both subspecies with different tissue specificity. For the five tissues in the two japonica-indica combinations, 4.6 ± 0.8 and 5.9 ± 1.5 % of highly expressed genes were global and changed-tissues DE genes, respectively. Changed-tissues DE genes varied in number between tissues, increasing linearly with the abundance of tissue specifically expressed genes in the tissue. Molecular evolution of global DE genes was rapid, unlike that of changed-tissues DE genes. Based on gene ontology, global and changed-tissues DE genes were different, having no common GO terms. Expression differences of most global DE genes were regulated by cis-eQTLs. Expression evolution of changed-tissues DE genes was rapid in tissue specifically expressed genes and those rapidly evolved changed-tissues DE genes were regulated not by cis-eQTLs, but by complicated trans-eQTLs. CONCLUSIONS: Global DE genes and changed-tissues DE genes had contrasting characteristics. The two contrasting types of DE genes provide possible explanations for the previous controversial conclusions about the relationships between molecular evolution and expression evolution of genes in different species, and the relationship between expression breadth and expression conservation in evolution.

PCR-based INDEL markers co-dominant between Oryza sativa, japonica cultivars and closely-related wild Oryza species.

Wild relatives genetically close to cultivars are precious genetic resources for plant breeding. Oryza rufipogon, O. barthii, O. glumaepatula, O. meridionalis and O. longistaminata are such wild species, and are also categorized as AA genome species based on their structural similarities. Chromosome segment substitution lines (CSSLs) are a powerful resource in breeding and genetics, and numerous rice CSSLs have been produced. This study aimed to develop DNA markers for evaluation of CSSLs directly by PCR and subsequent gel electrophoresis. We confirmed that up to 155 of 188 markers developed for detection of japonica-indica INDELs could also detect INDELs between rice cultivars and wild AA-species accessions. Percentages of applicable markers were higher in O. rufipogon accessions (61.7 to 85.6%), and lower in accessions of other four AA species (39.8 to 51.4%). These markers were distributed throughout the rice chromosomes, and will be useful for genotyping of CSSLs and other genetic resources derived from crosses between rice cultivars and closely related wild species.

RiTE database: a resource database for genus-wide rice genomics and evolutionary biology.

BACKGROUND: Comparative evolutionary analysis of whole genomes requires not only accurate annotation of gene space, but also proper annotation of the repetitive fraction which is often the largest component of most if not all genomes larger than 50 kb in size. RESULTS: Here we present the Rice TE database (RiTE-db)--a genus-wide collection of transposable elements and repeated sequences across 11 diploid species of the genus Oryza and the closely-related out-group Leersia perrieri. The database consists of more than 170,000 entries divided into three main types: (i) a classified and curated set of publicly-available repeated sequences, (ii) a set of consensus assemblies of highly-repetitive sequences obtained from genome sequencing surveys of 12 species; and (iii) a set of full-length TEs, identified and extracted from 12 whole genome assemblies. CONCLUSIONS: This is the first report of a repeat dataset that spans the majority of repeat variability within an entire genus, and one that includes complete elements as well as unassembled repeats. The database allows sequence browsing, downloading, and similarity searches. Because of the strategy adopted, the RiTE-db opens a new path to unprecedented direct comparative studies that span the entire nuclear repeat content of 15 million years of Oryza diversity.