Massive gene losses in Asian cultivated rice unveiled by comparative genome analysis

Integrated transcriptome and metabolome analysis unveil the response mechanism in wild rice (Zizania latifolia griseb.) against sheath rot infection

Sheath rot disease (SRD) is one of the most devastating diseases of Manchurian wild rice (MWR) (Zizania latifolia Griseb). Pilot experiments in our laboratory have shown that an MWR cultivar "Zhejiao NO.7"exhibits signs of SRD tolerance. To explore the responses of Zhejiao No. 7 to SRD infection, we used a combined transcriptome and metabolome analysis approach. A total of 136 differentially accumulated metabolites (DAMs, 114 up- and 22 down-accumulated in FA compared to CK) were detected. These up-accumulated metabolites were enriched in tryptophan metabolism, amino acid biosynthesis, flavonoids, and phytohormone signaling. Transcriptome sequencing results showed the differential expression of 11,280 genes (DEGs, 5,933 up-, and 5,347 downregulated in FA compared to CK). The genes expressed in tryptophan metabolism, amino acid biosynthesis, phytohormone biosynthesis and signaling, and reactive oxygen species homeostasis confirmed the metabolite results. In addition, genes related to the cell wall, carbohydrate metabolism, and plant-pathogen interaction (especially hypersensitive response) showed changes in expression in response to SRD infection. These results provide a basis for understanding the response mechanisms in MWR to FA attack that can be used for breeding SRD-tolerant MWR.

Genome-Wide Analysis of microRNAs and Their Target Genes in Dongxiang Wild Rice (Oryza rufipogon Griff.) Responding to Salt Stress

Rice (Oryza sativa) is a staple food for more than half of the world's population, and its production is critical for global food security. Moreover, rice yield decreases when exposed to abiotic stresses, such as salinity, which is one of the most detrimental factors for rice production. According to recent trends, as global temperatures continue to rise due to climate change, more rice fields may become saltier. Dongxiang wild rice (Oryza rufipogon Griff., DXWR) is a progenitor of cultivated rice and has a high tolerance to salt stress, making it useful for studying the regulatory mechanisms of salt stress tolerance. However, the regulatory mechanism of miRNA-mediated salt stress response in DXWR remains unclear. In this study, miRNA sequencing was performed to identify miRNAs and their putative target genes in response to salt stress in order to better understand the roles of miRNAs in DXWR salt stress tolerance. A total of 874 known and 476 novel miRNAs were identified, and the expression levels of 164 miRNAs were found to be significantly altered under salt stress. The stem-loop quantitative real-time PCR (qRT-PCR) expression levels of randomly selected miRNAs were largely consistent with the miRNA sequencing results, suggesting that the sequencing results were reliable. The gene ontology (GO) analysis indicated that the predicted target genes of salt-responsive miRNAs were involved in diverse biological pathways of stress tolerance. This study contributes to our understanding of DXWR salt tolerance mechanisms regulated by miRNAs and may ultimately improve salt tolerance in cultivated rice breeding using genetic methods in the future.

Potentially Useful Dwarfing or Semi-dwarfing Genes in Rice Breeding in Addition to the sd1 Gene

The "Green revolution" gene sd1 has been used widely in the breeding of modern rice varieties for over half a century. The application of this gene has increased rice yields and thereby supported a significant proportion of the global population. The use of a single gene, however, has raised concerns in the scientific community regarding its durability, especially given the bottleneck in genetic background and the need for large input of fertilizer. New dwarfing or semi-dwarfing genes are needed to alleviate our dependence on the sole "Green revolution" gene. In the past few years, several new dwarfing and semi-dwarfing genes as well as their mutants have been reported. Here, we provide an extensive review of the recent discoveries concerning newly identified genes that are potentially useful in rice breeding, including methods employed to create and effectively screen new rice mutants, the phenotypic characteristics of the new dwarfing and semi-dwarfing mutants, potential values of the new dwarfing and semi-dwarfing genes in rice breeding, and potential molecular mechanisms associated with the newly identified genes.

Expanding Gene-Editing Potential in Crop Improvement with Pangenomes

Pangenomes aim to represent the complete repertoire of the genome diversity present within a species or cohort of species, capturing the genomic structural variance between individuals. This genomic information coupled with phenotypic data can be applied to identify genes and alleles involved with abiotic stress tolerance, disease resistance, and other desirable traits. The characterisation of novel structural variants from pangenomes can support genome editing approaches such as Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR associated protein Cas (CRISPR-Cas), providing functional information on gene sequences and new target sites in variant-specific genes with increased efficiency. This review discusses the application of pangenomes in genome editing and crop improvement, focusing on the potential of pangenomes to accurately identify target genes for CRISPR-Cas editing of plant genomes while avoiding adverse off-target effects. We consider the limitations of applying CRISPR-Cas editing with pangenome references and potential solutions to overcome these limitations.

Evolutionary Significance of NHX Family and NHX1 in Salinity Stress Adaptation in the Genus Oryza

Rice (Oryza sativa), a staple crop for a substantial part of the world’s population, is highly sensitive to soil salinity; however, some wild Oryza relatives can survive in highly saline environments. Sodium/hydrogen antiporter (NHX) family members contribute to Na⁺ homeostasis in plants and play a major role in conferring salinity tolerance. In this study, we analyzed the evolution of NHX family members using phylogeny, conserved domains, tertiary structures, expression patterns, and physiology of cultivated and wild Oryza species to decipher the role of NHXs in salt tolerance in Oryza. Phylogenetic analysis showed that the NHX family can be classified into three subfamilies directly related to their subcellular localization: endomembrane, plasma membrane, and tonoplast (vacuolar subfamily, vNHX1). Phylogenetic and structural analysis showed that vNHX1s have evolved from streptophyte algae (e.g., Klebsormidium nitens) and are abundant and highly conserved in all major land plant lineages, including Oryza. Moreover, we showed that tissue tolerance is a crucial trait conferring tolerance to salinity in wild rice species. Higher Na⁺ accumulation and reduced Na⁺ effluxes in leaf mesophyll were observed in the salt-tolerant wild rice species O. alta, O. latifolia, and O. coarctata. Among the key genes affecting tissue tolerance, expression of NHX1 and SOS1/NHX7 exhibited significant correlation with salt tolerance among the rice species and cultivars. This study provides insights into the evolutionary origin of plant NHXs and their role in tissue tolerance of Oryza species and facilitates the inclusion of this trait during the development of salinity-tolerant rice cultivars.

Novel disease resistance gene paralogs created by CRISPR/Cas9 in soy

KEY MESSAGE

Novel disease resistance gene paralogues are generated by targeted chromosome cleavage of tandem duplicated NBS-LRR gene complexes and subsequent DNA repair in soybean. This study demonstrates accelerated diversification of innate immunity of plants using CRISPR. Nucleotide-binding-site-leucine-rich-repeat (NBS-LRR) gene families are key components of effector-triggered immunity. They are often arranged in tandem duplicated arrays in the genome, a configuration that is conducive to recombinations that will lead to new, chimeric genes. These rearrangements have been recognized as major sources of novel disease resistance phenotypes. Targeted chromosome cleavage by CRISPR/Cas9 can conceivably induce rearrangements and thus emergence of new resistance gene paralogues. Two NBS-LRR families of soy have been selected to demonstrate this concept: a four-copy family in the Rpp1 region (Rpp1L) and a large, complex locus, Rps1 with 22 copies. Copy-number variations suggesting large-scale, CRISPR/Cas9-mediated chromosome rearrangements in the Rpp1L and Rps1 complexes were detected in up to 58.8% of progenies of primary transformants using droplet-digital PCR. Sequencing confirmed development of novel, chimeric paralogs with intact open reading frames. These novel paralogs may confer new disease resistance specificities. This method to diversify innate immunity of plants by genome editing is readily applicable to other disease resistance genes or other repetitive loci.

Validation of a QTL for Grain Size and Weight Using an Introgression Line from a Cross between Oryza sativa and Oryza minuta

BACKGROUND

Grain size and weight are important target traits determining grain yield and quality in rice. Wild rice species possess substantial elite genes that can be served as an important resource for genetic improvement of rice. In this study, we identify and validate a novel QTL on chromosome 7 affecting the grain size and weight using introgression lines from cross of Oryza sativa and Oryza minuta.

RESULTS

An introgression line 'IL188' has been achieved from a wild species Oryza minuta (2n = 48, BBCC, W303) into O. sativa japonica Nipponbare. The F₂ and F_2:3 populations derived from a cross between IL188 and Nipponbare were used to map QTLs for five grain size traits, including grain length (GL), grain width (GW), grain length to width ratio (LWR), grain thickness (GT) and thousand grain weight (TGW). A total of 12 QTLs for the five grain traits were identified on chromosomes 1, 2, 3, 6, 7, and 8. The QTL-qGL7 controlling GL on chromosome 7 was detected stably in the F₂ and F_2:3 populations, and explained 15.09-16.30% of the phenotypic variance. To validate the effect of qGL7, eight residual heterozygous line (RHL) populations were developed through selfing four F_2:3 and four F_2:4 plants with different heterozygous segments for the target region. By further developing SSR and Indel markers in the target interval, qGL7 was delimited to a ~ 261 kb region between Indel marker Y7-12 and SSR marker Y7-38, which also showed significant effects on grain width and thousand grain weight. Comparing with the reference genome of Nipponbare, stop or frameshift mutations in the exon of the three putative genes LOC_Os07g36830, LOC_Os07g36900 and LOC_Os07g36910 encoding F-box domain-containing proteins may be the candidate genes for qGL7. Scanning electron microscopy analysis of the glume's epidermal cells showed that the cell length and width of NIL-qGL7^IL188 was higher than NIL-qGL7^Nip, indicating that qGL7 increases grain size and weight by regulating cell expansion.

CONCLUSIONS

In this study, we detected 12 QTLs regulating grain size and weight using an introgression line from a cross between Oryza sativa and Oryza minuta. Of these loci, we confirmed and delimited the qGL7 to a ~ 261 kb region. Three putative genes, LOC_Os07g36830, LOC_Os07g36900 and LOC_Os07g36910 encoding F-box domain-containing proteins may be the candidate genes for qGL7. These results provide a basis for map-based cloning of the qGL7 gene and useful information for marker assisted selection in rice grain quality improvement.

Genetics and Genomics of African Rice (Oryza glaberrima Steud) Domestication

African rice (Oryza glaberrima Steud) is one of the two independently domesticated rice species, the other one being Asian rice (Oryza sativa L.). Despite major progress being made in understanding the evolutionary and domestication history of African rice, key outstanding issues remain controversial. There appears to be an underlying difficulty in identifying the domestication centre and number of times the crop has been domesticated. Advances in genomics have provided unprecedented opportunities for understanding the genetic architecture of domestication related traits. For most of the domestication traits, the underlying genes and mutations have been identified. Comparative analysis of domestication genes between Asian and African rice has revealed that the two species went through an independent but convergent evolution process. The genetic and developmental basis of some of the domestic traits are conserved not only between Asian and African rice but also with other domesticated crop species. Analysis of genome data and its interpretation is emerging as a major challenge facing studies of domestication in African rice as key studies continue giving contradictory findings and conclusions. Insights obtained on the domestication of this species are vital for guiding crop improvement efforts.

RiceRelativesGD: a genomic database of rice relatives for rice research

Rice (Oryza sativa L.) is one of the most important crops worldwide. Its relatives, including phylogenetically related species of rice and paddy weeds with a similar ecological niche, can provide crucial genetic resources (such as resistance to biotic and abiotic stresses and high photosynthetic efficiency) for rice research. Although many rice genomic databases have been constructed, a database providing large-scale curated genomic data from rice relatives and offering specific gene resources is still lacking. Here, we present RiceRelativesGD, a user-friendly genomic database of rice relatives. RiceRelativesGD integrates large-scale genomic resources from 2 cultivated rice and 11 rice relatives, including 208 321 specific genes and 13 643 genes related to photosynthesis and responsive to external stimuli. Diverse bioinformatics tools are embedded in the database, which allow users to search, visualize and download the information of interest. To our knowledge, this is the first genomic database providing a centralized genetic resource of rice relatives. RiceRelativesGD will serve as a significant and comprehensive knowledgebase for the rice community.

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.

Genome-wide association of yield traits in a nested association mapping population of barley reveals new gene diversity for future breeding

To explore wild barley as a source of useful alleles for yield improvement in breeding, we have carried out a genome-wide association scan using the nested association mapping population HEB-25, which contains 25 diverse exotic barley genomes superimposed on an ~70% genetic background of cultivated barley. A total of 1420 HEB-25 lines were trialled for nine yield-related grain traits for 2 years in Germany and Scotland, with varying N fertilizer application. The phenotypic data were related to genotype scores for 5398 gene-based single nucleotide polymorphism (SNP) markers. A total of 96 quantitative trait locus (QTL) regions were identified across all measured traits, the majority of which co-localize with known major genes controlling flowering time (Ppd-H2, HvCEN, HvGI, VRN-H1, and VRN-H3) and spike morphology (VRS3, VRS1, VRS4, and INT-C) in barley. Fourteen QTL hotspots, with at least three traits coinciding, were also identified, several of which co-localize with barley orthologues of genes controlling grain dimensions in rice. Most of the allele effects are specific to geographical location and/or exotic parental genotype. This study shows the existence of beneficial alleles for yield-related traits in exotic barley germplasm and provides candidate alleles for future improvement of these traits by the breeder.

Effects of different treatments of fly ash and mining soil on growth and antioxidant protection of Indian wild rice

The aim of the present study was investigation of the effects of fly ash and mining soil on growth and antioxidant protection of two cultivars of Indian wild rice (Oryza nivara and Oryza rufipogon) for possible phytoremediation and restoration of metal-contaminated site. In this study, Indian wild rice showed significant changes in germination, growth, and biochemical parameters after exposure to different ratio of fly ash and mining soil with garden soil. There was significant reduction of germination, fresh weight, dry weight, leaf chlorophyll content, leaf area, Special Analysis Device Chlorophyll (SPAD) Index, proteins, and activities of antioxidant enzymes in both cultivars of the wild rice grown in 100% fly ash and mining soil compared to the plants grown in 100% garden soil. Results from this study showed that in both cultivars of wild rice, all growth and antioxidant parameters increased when grown in 50% fly ash and mining soil. Taken together, Indian wild rice has the capacity to tolerate 50% of fly ash and mining soil, and can be considered as a good candidate for possible phytoremediation of contaminated soils.

A walk on the wild side: Oryza species as source for rice abiotic stress tolerance

Oryza sativa, the common cultivated rice, is one of the most important crops for human consumption, but production is increasingly threatened by abiotic stresses. Although many efforts have resulted in breeding rice cultivars that are relatively tolerant to their local environments, climate changes and population increase are expected to soon call for new, fast generation of stress tolerant rice germplasm, and current within-species rice diversity might not be enough to overcome such needs. The Oryza genus contains other 23 wild species, with only Oryza glaberrima being also domesticated. Rice domestication was performed with a narrow genetic diversity, and the other Oryza species are a virtually untapped genetic resource for rice stress tolerance improvement. Here we review the origin of domesticated Oryza sativa from wild progenitors, the ecological and genomic diversity of the Oryza genus, and the stress tolerance variation observed for wild Oryza species, including the genetic basis underlying the tolerance mechanisms found. The summary provided here is important to indicate how we should move forward to unlock the full potential of these germplasms for rice improvement.

Novel Sources of Witchweed (Striga) Resistance from Wild Sorghum Accessions

Sorghum is a major food staple in sub-Saharan Africa (SSA), but its production is constrained by the parasitic plant Striga that attaches to the roots of many cereals crops and causes severe stunting and loss of yield. Away from cultivated farmland, wild sorghum accessions grow as weedy plants and have shown remarkable immunity to Striga. We sought to determine the extent of the resistance to Striga in wild sorghum plants. Our screening strategy involved controlled laboratory assays of rhizotrons, where we artificially infected sorghum with Striga, as well as field experiments at three sites, where we grew sorghum with a natural Striga infestation. We tested the resistance response of seven accessions of wild sorghum of the aethiopicum, drummondii, and arundinaceum races against N13, which is a cultivated Striga resistant landrace. The susceptible control was farmer-preferred variety, Ochuti. From the laboratory experiments, we found three wild sorghum accessions (WSA-1, WSE-1, and WSA-2) that had significantly higher resistance than N13. These accessions had the lowest Striga biomass and the fewest and smallest Striga attached to them. Further microscopic and histological analysis of attached Striga haustorium showed that wild sorghum accessions hindered the ingression of Striga haustorium into the host endodermis. In one of the resistant accessions (WSE-1), host and parasite interaction led to the accumulation of large amounts of secondary metabolites that formed a dark coloration at the interphase. Field experiments confirmed the laboratory screening experiments in that these same accessions were found to have resistance against Striga. In the field, wild sorghum had low Area under the Striga Number Progressive curve (AUSNPC), which measures emergence of Striga from a host over time. We concluded that wild sorghum accessions are an important reservoir for Striga resistance that could be used to expand the genetic basis of cultivated sorghum for resistance to the parasite.

Genome-wide identification of conserved microRNA and their response to drought stress in Dongxiang wild rice (Oryza rufipogon Griff.)

OBJECTIVES

To identify drought stress-responsive conserved microRNA (miRNA) from Dongxiang wild rice (Oryza rufipogon Griff., DXWR) on a genome-wide scale, high-throughput sequencing technology was used to sequence libraries of DXWR samples, treated with and without drought stress.

RESULTS

505 conserved miRNAs corresponding to 215 families were identified. 17 were significantly down-regulated and 16 were up-regulated under drought stress. Stem-loop qRT-PCR revealed the same expression patterns as high-throughput sequencing, suggesting the accuracy of the sequencing result was high. Potential target genes of the drought-responsive miRNA were predicted to be involved in diverse biological processes. Furthermore, 16 miRNA families were first identified to be involved in drought stress response from plants.

CONCLUSION

These results present a comprehensive view of the conserved miRNA and their expression patterns under drought stress for DXWR, which will provide valuable information and sequence resources for future basis studies.

To grow or not to grow: a stressful decision for plants

Progress in improving abiotic stress tolerance of crop plants using classic breeding and selection approaches has been slow. This has generally been blamed on the lack of reliable traits and phenotyping methods for stress tolerance. In crops, abiotic stress tolerance is most often measured in terms of yield-capacity under adverse weather conditions. "Yield" is a complex trait and is determined by growth and developmental processes which are controlled by environmental signals throughout the life cycle of the plant. The use of model systems has allowed us to gradually unravel how plants grow and develop, but our understanding of the flexibility and opportunistic nature of plant development and its capacity to adapt growth to environmental cues is still evolving. There is genetic variability for the capacity to maintain yield and productivity under abiotic stress conditions in crop plants such as cereals. Technological progress in various domains has made it increasingly possible to mine that genetic variability and develop a better understanding about the basic mechanism of plant growth and abiotic stress tolerance. The aim of this paper is not to give a detailed account of all current research progress, but instead to highlight some of the current research trends that may ultimately lead to strategies for stress-proofing crop species. The focus will be on abiotic stresses that are most often associated with climate change (drought, heat and cold) and those crops that are most important for human nutrition, the cereals.

Rapid diversification of five Oryza AA genomes associated with rice adaptation

Comparative genomic analyses among closely related species can greatly enhance our understanding of plant gene and genome evolution. We report de novo-assembled AA-genome sequences for Oryza nivara, Oryza glaberrima, Oryza barthii, Oryza glumaepatula, and Oryza meridionalis. Our analyses reveal massive levels of genomic structural variation, including segmental duplication and rapid gene family turnover, with particularly high instability in defense-related genes. We show, on a genomic scale, how lineage-specific expansion or contraction of gene families has led to their morphological and reproductive diversification, thus enlightening the evolutionary process of speciation and adaptation. Despite strong purifying selective pressures on most Oryza genes, we documented a large number of positively selected genes, especially those genes involved in flower development, reproduction, and resistance-related processes. These diversifying genes are expected to have played key roles in adaptations to their ecological niches in Asia, South America, Africa and Australia. Extensive variation in noncoding RNA gene numbers, function enrichment, and rates of sequence divergence might also help account for the different genetic adaptations of these rice species. Collectively, these resources provide new opportunities for evolutionary genomics, numerous insights into recent speciation, a valuable database of functional variation for crop improvement, and tools for efficient conservation of wild rice germplasm.

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.

Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data

Background

Rice research has been enabled by access to the high quality reference genome sequence generated in 2005 by the International Rice Genome Sequencing Project (IRGSP). To further facilitate genomic-enabled research, we have updated and validated the genome assembly and sequence for the Nipponbare cultivar of Oryza sativa (japonica group).

Results

The Nipponbare genome assembly was updated by revising and validating the minimal tiling path of clones with the optical map for rice. Sequencing errors in the revised genome assembly were identified by re-sequencing the genome of two different Nipponbare individuals using the Illumina Genome Analyzer II/IIx platform. A total of 4,886 sequencing errors were identified in 321 Mb of the assembled genome indicating an error rate in the original IRGSP assembly of only 0.15 per 10,000 nucleotides. A small number (five) of insertions/deletions were identified using longer reads generated using the Roche 454 pyrosequencing platform. As the re-sequencing data were generated from two different individuals, we were able to identify a number of allelic differences between the original individual used in the IRGSP effort and the two individuals used in the re-sequencing effort. The revised assembly, termed Os-Nipponbare-Reference-IRGSP-1.0, is now being used in updated releases of the Rice Annotation Project and the Michigan State University Rice Genome Annotation Project, thereby providing a unified set of pseudomolecules for the rice community.

Conclusions

A revised, error-corrected, and validated assembly of the Nipponbare cultivar of rice was generated using optical map data, re-sequencing data, and manual curation that will facilitate on-going and future research in rice. Detection of polymorphisms between three different Nipponbare individuals highlights that allelic differences between individuals should be considered in diversity studies.

Electronic supplementary material

The online version of this article (doi:10.1186/1939-8433-6-4) contains supplementary material, which is available to authorized users.

Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data

BACKGROUND

Rice research has been enabled by access to the high quality reference genome sequence generated in 2005 by the International Rice Genome Sequencing Project (IRGSP). To further facilitate genomic-enabled research, we have updated and validated the genome assembly and sequence for the Nipponbare cultivar of Oryza sativa (japonica group).

RESULTS

The Nipponbare genome assembly was updated by revising and validating the minimal tiling path of clones with the optical map for rice. Sequencing errors in the revised genome assembly were identified by re-sequencing the genome of two different Nipponbare individuals using the Illumina Genome Analyzer II/IIx platform. A total of 4,886 sequencing errors were identified in 321 Mb of the assembled genome indicating an error rate in the original IRGSP assembly of only 0.15 per 10,000 nucleotides. A small number (five) of insertions/deletions were identified using longer reads generated using the Roche 454 pyrosequencing platform. As the re-sequencing data were generated from two different individuals, we were able to identify a number of allelic differences between the original individual used in the IRGSP effort and the two individuals used in the re-sequencing effort. The revised assembly, termed Os-Nipponbare-Reference-IRGSP-1.0, is now being used in updated releases of the Rice Annotation Project and the Michigan State University Rice Genome Annotation Project, thereby providing a unified set of pseudomolecules for the rice community.

CONCLUSIONS

A revised, error-corrected, and validated assembly of the Nipponbare cultivar of rice was generated using optical map data, re-sequencing data, and manual curation that will facilitate on-going and future research in rice. Detection of polymorphisms between three different Nipponbare individuals highlights that allelic differences between individuals should be considered in diversity studies.

Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data

BACKGROUND

Rice research has been enabled by access to the high quality reference genome sequence generated in 2005 by the International Rice Genome Sequencing Project (IRGSP). To further facilitate genomic-enabled research, we have updated and validated the genome assembly and sequence for the Nipponbare cultivar of Oryza sativa (japonica group).

RESULTS

The Nipponbare genome assembly was updated by revising and validating the minimal tiling path of clones with the optical map for rice. Sequencing errors in the revised genome assembly were identified by re-sequencing the genome of two different Nipponbare individuals using the Illumina Genome Analyzer II/IIx platform. A total of 4,886 sequencing errors were identified in 321 Mb of the assembled genome indicating an error rate in the original IRGSP assembly of only 0.15 per 10,000 nucleotides. A small number (five) of insertions/deletions were identified using longer reads generated using the Roche 454 pyrosequencing platform. As the re-sequencing data were generated from two different individuals, we were able to identify a number of allelic differences between the original individual used in the IRGSP effort and the two individuals used in the re-sequencing effort. The revised assembly, termed Os-Nipponbare-Reference-IRGSP-1.0, is now being used in updated releases of the Rice Annotation Project and the Michigan State University Rice Genome Annotation Project, thereby providing a unified set of pseudomolecules for the rice community.

CONCLUSIONS

A revised, error-corrected, and validated assembly of the Nipponbare cultivar of rice was generated using optical map data, re-sequencing data, and manual curation that will facilitate on-going and future research in rice. Detection of polymorphisms between three different Nipponbare individuals highlights that allelic differences between individuals should be considered in diversity studies.

Rice Annotation Project Database (RAP-DB): an integrative and interactive database for rice genomics

The Rice Annotation Project Database (RAP-DB, http://rapdb.dna.affrc.go.jp/) has been providing a comprehensive set of gene annotations for the genome sequence of rice, Oryza sativa (japonica group) cv. Nipponbare. Since the first release in 2005, RAP-DB has been updated several times along with the genome assembly updates. Here, we present our newest RAP-DB based on the latest genome assembly, Os-Nipponbare-Reference-IRGSP-1.0 (IRGSP-1.0), which was released in 2011. We detected 37,869 loci by mapping transcript and protein sequences of 150 monocot species. To provide plant researchers with highly reliable and up to date rice gene annotations, we have been incorporating literature-based manually curated data, and 1,626 loci currently incorporate literature-based annotation data, including commonly used gene names or gene symbols. Transcriptional activities are shown at the nucleotide level by mapping RNA-Seq reads derived from 27 samples. We also mapped the Illumina reads of a Japanese leading japonica cultivar, Koshihikari, and a Chinese indica cultivar, Guangluai-4, to the genome and show alignments together with the single nucleotide polymorphisms (SNPs) and gene functional annotations through a newly developed browser, Short-Read Assembly Browser (S-RAB). We have developed two satellite databases, Plant Gene Family Database (PGFD) and Integrative Database of Cereal Gene Phylogeny (IDCGP), which display gene family and homologous gene relationships among diverse plant species. RAP-DB and the satellite databases offer simple and user-friendly web interfaces, enabling plant and genome researchers to access the data easily and facilitating a broad range of plant research topics.

Cis-regulatory signatures of orthologous stress-associated bZIP transcription factors from rice, sorghum and Arabidopsis based on phylogenetic footprints

Background

The potential contribution of upstream sequence variation to the unique features of orthologous genes is just beginning to be unraveled. A core subset of stress-associated bZIP transcription factors from rice (Oryza sativa) formed ten clusters of orthologous groups (COG) with genes from the monocot sorghum (Sorghum bicolor) and dicot Arabidopsis (Arabidopsis thaliana). The total cis-regulatory information content of each stress-associated COG was examined by phylogenetic footprinting to reveal ortholog-specific, lineage-specific and species-specific conservation patterns.

Results

The most apparent pattern observed was the occurrence of spatially conserved ‘core modules’ among the COGs but not among paralogs. These core modules are comprised of various combinations of two to four putative transcription factor binding site (TFBS) classes associated with either developmental or stress-related functions. Outside the core modules are specific stress (ABA, oxidative, abiotic, biotic) or organ-associated signals, which may be functioning as ‘regulatory fine-tuners’ and further define lineage-specific and species-specific cis-regulatory signatures. Orthologous monocot and dicot promoters have distinct TFBS classes involved in disease and oxidative-regulated expression, while the orthologous rice and sorghum promoters have distinct combinations of root-specific signals, a pattern that is not particularly conserved in Arabidopsis.

Conclusions

Patterns of cis-regulatory conservation imply that each ortholog has distinct signatures, further suggesting that they are potentially unique in a regulatory context despite the presumed conservation of broad biological function during speciation. Based on the observed patterns of conservation, we postulate that core modules are likely primary determinants of basal developmental programming, which may be integrated with and further elaborated by additional intrinsic or extrinsic signals in conjunction with lineage-specific or species-specific regulatory fine-tuners. This synergy may be critical for finer-scale spatio-temporal regulation, hence unique expression profiles of homologous transcription factors from different species with distinct zones of ecological adaptation such as rice, sorghum and Arabidopsis. The patterns revealed from these comparisons set the stage for further empirical validation by functional genomics.

Independent domestication of Asian rice followed by gene flow from japonica to indica

Results from studies on the domestication process of Asian rice Oryza sativa have been controversial because of its complicated evolutionary history. Previous studies have yielded two alternative hypotheses about the origin(s) of the two major groups of O. sativa: japonica and indica. One study proposes a single common wild ancestor, whereas the other suggests that there were multiple domestication events of different types of wild rice. Here, we provide clear evidence of the independent domestication of japonica and indica obtained via high-throughput sequencing and a large-scale comparative analysis of two wild rice accessions (W1943 and W0106) and two cultivars (a japonica cultivar called "Nipponbare" and an indica cultivar called "Guangluai-4"). The different domestication processes of the two cultivar groups appear to have led to distinct patterns of molecular evolution in protein-coding regions. The intensity of purifying selection was relaxed only in the japonica group, possibly because of a bottleneck effect. Moreover, a genome-wide comparison between Nipponbare, Guangluai-4, and another indica cultivar (93-11) suggests multiple hybridization events between japonica and indica, both before and after the divergence of the indica cultivars. We found that a large amount of genomic DNA, including domestication-related genes, was transferred from japonica to indica, which might have been important in the development of modern rice. Our study provides an overview of the dynamic process of Asian rice domestication, including independent domestication events and subsequent gene flow.

Cell death mediated by the N-terminal domains of a unique and highly conserved class of NB-LRR protein

Plant genomes encode large numbers of nucleotide-binding, leucine-rich repeat (NB-LRR) proteins, many of which are active in pathogen detection and defense response induction. NB-LRR proteins fall into two broad classes: those with a Toll and interleukin-1 receptor (TIR) domain at their N-terminus and those with a coiled-coil (CC) domain at the N-terminus. Within CC-NB-LRR-encoding genes, one basal clade is distinguished by having CC domains resembling the Arabidopsis thaliana RPW8 protein, which we refer to as CCR domains. Here, we show that CCR-NB-LRR-encoding genes are present in the genomes of all higher plants surveyed, and that they comprise two distinct subgroups: one typified by the Nicotiana benthamiana N-required gene 1 (NRG1) protein and the other typified by the Arabidopsis activated disease resistance gene 1 (ADR1) protein. We further report that, in contrast to CC-NB-LRR proteins, the CCR domains of both NRG1- and ADR1-like proteins are sufficient for the induction of defense responses, and that this activity appears to be SGT1-independent. Additionally, we report the apparent absence of both NRG1 homologs and TIR-NB-LRR-encoding genes from the dicot Aquilegia coerulea and the dicotyledonous order Lamiales as well as from monocotyledonous species. This strong correlation in occurrence is suggestive of a functional relationship between these two classes of NB-LRR proteins.

Distinct evolutionary patterns of Oryza glaberrima deciphered by genome sequencing and comparative analysis

Here we present the genomic sequence of the African cultivated rice, Oryza glaberrima, and compare these data with the genome sequence of Asian cultivated rice, Oryza sativa. We obtained gene-enriched sequences of O. glaberrima that correspond to about 25% of the gene regions of the O. sativa (japonica) genome by methylation filtration and subtractive hybridization of repetitive sequences. While patterns of amino acid changes did not differ between the two species in terms of the biochemical properties, genes of O. glaberrima generally showed a larger synonymous-nonsynonymous substitution ratio, suggesting that O. glaberrima has undergone a genome-wide relaxation of purifying selection. We further investigated nucleotide substitutions around splice sites and found that eight genes of O. sativa experienced changes at splice sites after the divergence from O. glaberrima. These changes produced novel introns that partially truncated functional domains, suggesting that these newly emerged introns affect gene function. We also identified 2451 simple sequence repeats (SSRs) from the genomes of O. glaberrima and O. sativa. Although tri-nucleotide repeats were most common among the SSRs and were overrepresented in the protein-coding sequences, we found that selection against indels of tri-nucleotide repeats was relatively weak in both African and Asian rice. Our genome-wide sequencing of O. glaberrima and in-depth analyses provide rice researchers not only with useful genomic resources for future breeding but also with new insights into the genomic evolution of the African and Asian rice species.

Gene tree discordance of wild and cultivated Asian rice deciphered by genome-wide sequence comparison

Although a large number of genes are expected to correctly solve a phylogenetic relationship, inconsistent gene tree topologies have been observed. This conflicting evidence in gene tree topologies, known as gene tree discordance, becomes increasingly important as advanced sequencing technologies produce an enormous amount of sequence information for phylogenomic studies among closely related species. Here, we aim to characterize the gene tree discordance of the Asian cultivated rice Oryza sativa and its progenitor, O. rufipogon, which will be an ideal case study of gene tree discordance. Using genome and cDNA sequences of O. sativa and O. rufipogon, we have conducted the first in-depth analyses of gene tree discordance in Asian rice. Our comparison of full-length cDNA sequences of O. rufipogon with the genome sequences of the japonica and indica cultivars of O. sativa revealed that 60% of the gene trees showed a topology consistent with the expected one, whereas the remaining genes supported significantly different topologies. Moreover, the proportions of the topologies deviated significantly from expectation, suggesting at least one hybridization event between the two subgroups of O. sativa, japonica and indica. In fact, a genome-wide alignment between japonica and indica indicated that significant portions of the indica genome are derived from japonica. In addition, literature concerning the pedigree of the indica cultivar strongly supported the hybridization hypothesis. Our molecular evolutionary analyses deciphered complicated evolutionary processes in closely related species. They also demonstrated the importance of gene tree discordance in the era of high-speed DNA sequencing.

Numbers of genes in the NBS and RLK families vary by more than four-fold within a plant species and are regulated by multiple factors

Many genes exist in the form of families; however, little is known about their size variation, evolution and biology. Here, we present the size variation and evolution of the nucleotide-binding site (NBS)-encoding gene family and receptor-like kinase (RLK) gene family in Oryza, Glycine and Gossypium. The sizes of both families vary by numeral fold, not only among species, surprisingly, also within a species. The size variations of the gene families are shown to correlate with each other, indicating their interactions, and driven by natural selection, artificial selection and genome size variation, but likely not by polyploidization. The numbers of genes in the families in a polyploid species are similar to those of one of its diploid donors, suggesting that polyploidization plays little roles in the expansion of the gene families and that organisms tend not to maintain their ‘surplus’ genes in the course of evolution. Furthermore, it is found that the size variations of both gene families are associated with organisms’ phylogeny, suggesting their roles in speciation and evolution. Since both selection and speciation act on organism’s morphological, physiological and biological variation, our results indicate that the variation of gene family size provides a source of genetic variation and evolution.