Abstract

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.

Introduction

Crop breeding is crucial for guaranteeing food security and sustainable human population growth. Cultivated rice (Oryza sativa L.) is one of the most important crops worldwide and a model species for functional genomics of monocots. During the period of domestication, cultivated rice has lost many genes controlling important agronomic traits, such as resistance to abiotic and biotic stresses, which are potentially very useful for modern rice breeding (1,2). Many studies have revealed that genes regulating these traits are maintained in two closely related groups of the cultivated rice (3–5). One group includes species such as O. rufipogon that is phylogenetically closed to the cultivated rice (6). The other group consists of paddy weeds that have a similar ecological niche as rice and are highly competitive and readily adapt to the agroecosystem (5,7). Species of both groups are potential gene resources for modern rice molecular breeding programs aiming for improvement of agronomic traits.

With the availability of sequencing technologies, genome-based molecular approaches can increase the efficiency of rice breeding (8–12). The first rice genome was sequenced in 2002, making it the first crop genome to be deciphered (13,14). To date, dozens of de novo sequenced genomes of rice and its relatives have been published (7,15–20). According to the statistics from the TimeTree database (21), 13 of the 46 species in the Oryzoideae subfamily have already been sequenced and/or reconstructed to the pseudo-chromosome level, including seven rice relative species sharing the AA genome (16,20,22,23). Apart from the species phylogenetically close to O. sativa, the genome sequence of Echinochloa crus-galli, one of the most pervasive paddy weeds, has also been investigated (7). It demonstrated that E. crus-galli interacts with O. sativa reciprocally through allelopathy, which gives us an example of crop–weed interaction on the genomic level (7). Genomic data of rice and its relatives have provided significant and important insights into the molecular mechanisms for stress resistance, regulation of agronomic traits and interactions between cultivated rice and other plants (e.g. allelopathy) in its immediate environment (5).

To facilitate the accessibility of information from such enormous genomic data, several Oryza genome databases have been created to accommodate the genome data and various other types of data. The current online genomic resources of rice can be roughly divided into three categories depending on the main resources included. One is the de novo genome data (e.g. RAP-DB, MSU-RGAP, RIGW, RIS and RPAN) (24–27), another is rice genomic diversity data (e.g. SNP-Seek, RiceVarMap and OryzaGenome) (28–30) and the third is integrated databases (e.g. IC4R, Oryzabase and Gramene) (31–33). However, numerous genomic information of the species related to cultivated rice is still waiting to be incorporated into databases.

In this work, we constructed a new and user-friendly database termed RiceRelativesGD (http://ibi.zju.edu.cn/ricerelativesgd/), with the aim to serve it as a comprehensive genomic resource of rice relatives useful for rice breeding. We firstly integrated publicly available genomic resources from 2 cultivated rice and 11 rice relatives and identified 208 321 specific genes. We also identified valuable genomic resource for rice research, such as stress-related genes and photosynthesis genes, which is key to breed rice varieties with stress resistant and high efficiency of photosynthesis. Finally, practical bioinformatics online tools were provided to allow researchers to analyze genes with potential value for rice community.

Materials and methods

Data collection, classification and annotation

The current database included genomic datasets from 13 rice relatives that are publically available (Table 1). From the raw protein files of 13 genomes, the longest protein of each orthologous gene termed as ‘primary protein’ was extracted for gene family clustering analysis. Based on the Markov Cluster algorithm, 34 570 gene families were identified using Orthofinder v2.2.7 with sequence search program ‘diamond’ (34). And according to the gene family clustering results, genes present in multiple species were defined as ‘multi-species family genes’, genes present only in a single species were defined as ‘species-specific family genes’, while genes that could not be clustered to any gene family were defined as ‘orphan genes’. Species-specific family genes, orphan genes and multiple-species family genes that could not be clustered with O. sativa genes in the gene family analysis were further defined as ‘specific genes’ or, in other words, specific genes in RiceRelativesGD refer to genes without paralogs in O. sativa (japonica group) or O. sativa (indica group).

Table 1

List of sequenced genomes of rice and its relatives

SpeciesGenome typesGenome sizePredicted gene numberAccessionReference
(Estimated/assembly)
O. sativa (japonica group)AA (2n = 24)420 Mb/390 Mb35 825Nipponbare(15)
O. sativa (indica group)AA(2n = 24)−−/396 Mb38 729R498(19)
O. rufipogonAA(2n = 24)−−/338 Mb37 071W1943(20)
O. nivaraAA(2n = 24)−−/338 Mb36 313IRGC100897(20)
O. glaberrimaAA(2n = 24)−−/316 Mb33 164CG14(17)
O. barthiiAA(2n = 24)−−/308 Mb34 575IRGC105608(20)
O. glumaepatulaAA(2n = 24)−−/373 Mb38 149GEN1233_2(20)
O. meridionalisAA(2n = 24)−−/336 Mb34 897W2112(20)
O. puntacaBB(2n = 24)−−/394 Mb31 762IRGC105690(20)
O. brachyanthaFF(2n = 24)297 Mb/261 Mb32 037IRGC101232(16)
Leersia perrieri−−/(2n = 24)−−/267 Mb29 078IRGC105164(20)
Zizania latifolia−−/(2n = 34)594 Mb/604 Mb43 703HSD2(18)
Echinochloa crus-galli−−/(2n = 6x = 54)1.40Gb/1.27Gb108 771STB08(7)
SpeciesGenome typesGenome sizePredicted gene numberAccessionReference
(Estimated/assembly)
O. sativa (japonica group)AA (2n = 24)420 Mb/390 Mb35 825Nipponbare(15)
O. sativa (indica group)AA(2n = 24)−−/396 Mb38 729R498(19)
O. rufipogonAA(2n = 24)−−/338 Mb37 071W1943(20)
O. nivaraAA(2n = 24)−−/338 Mb36 313IRGC100897(20)
O. glaberrimaAA(2n = 24)−−/316 Mb33 164CG14(17)
O. barthiiAA(2n = 24)−−/308 Mb34 575IRGC105608(20)
O. glumaepatulaAA(2n = 24)−−/373 Mb38 149GEN1233_2(20)
O. meridionalisAA(2n = 24)−−/336 Mb34 897W2112(20)
O. puntacaBB(2n = 24)−−/394 Mb31 762IRGC105690(20)
O. brachyanthaFF(2n = 24)297 Mb/261 Mb32 037IRGC101232(16)
Leersia perrieri−−/(2n = 24)−−/267 Mb29 078IRGC105164(20)
Zizania latifolia−−/(2n = 34)594 Mb/604 Mb43 703HSD2(18)
Echinochloa crus-galli−−/(2n = 6x = 54)1.40Gb/1.27Gb108 771STB08(7)
Table 1

List of sequenced genomes of rice and its relatives

SpeciesGenome typesGenome sizePredicted gene numberAccessionReference
(Estimated/assembly)
O. sativa (japonica group)AA (2n = 24)420 Mb/390 Mb35 825Nipponbare(15)
O. sativa (indica group)AA(2n = 24)−−/396 Mb38 729R498(19)
O. rufipogonAA(2n = 24)−−/338 Mb37 071W1943(20)
O. nivaraAA(2n = 24)−−/338 Mb36 313IRGC100897(20)
O. glaberrimaAA(2n = 24)−−/316 Mb33 164CG14(17)
O. barthiiAA(2n = 24)−−/308 Mb34 575IRGC105608(20)
O. glumaepatulaAA(2n = 24)−−/373 Mb38 149GEN1233_2(20)
O. meridionalisAA(2n = 24)−−/336 Mb34 897W2112(20)
O. puntacaBB(2n = 24)−−/394 Mb31 762IRGC105690(20)
O. brachyanthaFF(2n = 24)297 Mb/261 Mb32 037IRGC101232(16)
Leersia perrieri−−/(2n = 24)−−/267 Mb29 078IRGC105164(20)
Zizania latifolia−−/(2n = 34)594 Mb/604 Mb43 703HSD2(18)
Echinochloa crus-galli−−/(2n = 6x = 54)1.40Gb/1.27Gb108 771STB08(7)
SpeciesGenome typesGenome sizePredicted gene numberAccessionReference
(Estimated/assembly)
O. sativa (japonica group)AA (2n = 24)420 Mb/390 Mb35 825Nipponbare(15)
O. sativa (indica group)AA(2n = 24)−−/396 Mb38 729R498(19)
O. rufipogonAA(2n = 24)−−/338 Mb37 071W1943(20)
O. nivaraAA(2n = 24)−−/338 Mb36 313IRGC100897(20)
O. glaberrimaAA(2n = 24)−−/316 Mb33 164CG14(17)
O. barthiiAA(2n = 24)−−/308 Mb34 575IRGC105608(20)
O. glumaepatulaAA(2n = 24)−−/373 Mb38 149GEN1233_2(20)
O. meridionalisAA(2n = 24)−−/336 Mb34 897W2112(20)
O. puntacaBB(2n = 24)−−/394 Mb31 762IRGC105690(20)
O. brachyanthaFF(2n = 24)297 Mb/261 Mb32 037IRGC101232(16)
Leersia perrieri−−/(2n = 24)−−/267 Mb29 078IRGC105164(20)
Zizania latifolia−−/(2n = 34)594 Mb/604 Mb43 703HSD2(18)
Echinochloa crus-galli−−/(2n = 6x = 54)1.40Gb/1.27Gb108 771STB08(7)
The flow diagram to show design and construction of RiceRelativeGD. (A) Data collection: genomic data from the cultivated rice and its relatives. (B) Data analysis: a flow chart to show the analyses used in the curation of the data included in the database, including the identification of specific genes and genes related to photosynthesis and stress tolerance. (C) Database construction: tools used in construction of the database and building functional modules of the database.
Figure 1

The flow diagram to show design and construction of RiceRelativeGD. (A) Data collection: genomic data from the cultivated rice and its relatives. (B) Data analysis: a flow chart to show the analyses used in the curation of the data included in the database, including the identification of specific genes and genes related to photosynthesis and stress tolerance. (C) Database construction: tools used in construction of the database and building functional modules of the database.

Detailed information on each gene, including organisms of origin, genomic location, family ID, gene structure, function descriptions and sequences, was displayed on the website. We used two approaches to annotate the exact function of each gene and its family. First, the primary protein sequences from 13 species were annotated with Pfam domain and GO terms using InterProscan v5.24–63.0 (35) and Mercator v4.0 with the default settings (36). Particularly, to better annotate gene functions of the genes included in RiceRelativesGD, we collected and used GO annotations and functional description of O. sativa (japonica group) from Rap-DB and MSU-RGAP. Second, we aligned the primary protein sequences to Swiss-prot dataset from Uniport (37) using ncbi-blast+ 2.6.0 with ‘-evalue 1e-5’ (38) and the best hit was extracted. The annotation results from both methods were integrated into our database. For each group of specific genes from different rice relatives, GOATOOLS (39) was applied to perform GO enrichment analysis with the GO annotation from InterProscan.

Phylogenetic tree and genomic presence and absence variations

Based on the paralogs and orthologs identified by Orthofinder v2.2.7 (34), single-copy genes were identified and their primary protein sequences were aligned using MAFFT v7.310 (40). The tree was built using FastTree 2.1.7 with the default setting (41) and was drawn with ITOL (42).

To display genomic PAV in the database, genomes of the 11 rice relatives were aligned to japonica or indica rice genome, using the Nucmer program in the MUMer 3.23 package with the default setting (43); the alignment results were further filtered using the delta-filter program with parameter ‘-q’ and transformed into human readable format using the show-coords program with parameter ‘-lrTH’ in the MUMer 3.23 package (43). Finally, gff3 files were generated from the obtained results and visualized in the genome browser.

Implementation of RiceRelativesGD

This database was built as a web-based system consisting of two major parts: one is data storage and management, and the other is high-level web interfaces displaying visualization functions (Figure 1). The backend of this system was constructed with MySQL 15.1, Perl 5.16 and PHP 5.6. HTML5, and JavaScript operating on the Apache web server was used to construct web interfaces. For the implementation of online tools in the database, several bioinformatics tools and external databases were integrated in this system. In the ‘search’ module, BLAST+ 2.2.29 (38) was used to perform sequence search. In the ‘Genome browse’ module, JBrowse 1.13.1 (44) was used to visualize genomic PAV (presence and absence variations) between cultivated rice and its relatives and specific genes. In the ‘build phyogenetic tree’ module, Mafft v7.310 (40), FastTree 2.1.7 (41) and PhyD3 (45) were used to perform sequence alignment, phylogenetic tree building and visualization.

Statistics of the data included in RiceRelativeGD. (A) Stacked bar diagram to show the proportion of the three types of genes in each species. (B) Bar diagram to show the number of stress-related genes and photosynthesis genes in each species. (C) Venn diagram to show the number of genes specific to japonica or indica rice.
Figure 2

Statistics of the data included in RiceRelativeGD. (A) Stacked bar diagram to show the proportion of the three types of genes in each species. (B) Bar diagram to show the number of stress-related genes and photosynthesis genes in each species. (C) Venn diagram to show the number of genes specific to japonica or indica rice.

A case study for the application of RiceRelativesGD (showing the related information of the EC_V6.g004014 gene). (A) General information about the gene. (B) Functional annotation of the gene. (C) Description of the gene family to which EC_V6.g004014 belongs. (D) Phylogenetic tree of the genes in the gene family to which EC_V6.g004014 belongs.
Figure 3

A case study for the application of RiceRelativesGD (showing the related information of the EC_V6.g004014 gene). (A) General information about the gene. (B) Functional annotation of the gene. (C) Description of the gene family to which EC_V6.g004014 belongs. (D) Phylogenetic tree of the genes in the gene family to which EC_V6.g004014 belongs.

Results

The overall design of RiceRelativesGD

The overall design of RiceRelativesGD is shown in Figure 1. To better help users browse the data in the RiceRelativesGD, we set up four sections in the ‘Browse module’, including ‘Specific genes’, ‘Gene families’, ‘Stress-related genes’ and ‘Photosynthesis genes’. In the ‘Specific genes’ section, users can browse genes specific to indica rice (R498) or japonica rice (Nipponbare) by species (Figure 1). In the interface of each species, detailed information about the gene, such as gene family entry, statistics results and GO enrichment results, are shown. In the ‘Stress-related genes’ and ‘Photosynthesis genes’ sections, genes are further classified by their sub-classifications. For example, in the ‘Herbicide’ of ‘Stress-related genes’ sections, genes are further listed in four categories, including Cytochrome P450 monooxygenases, Glutathione S-transferases, ABC transporters and Glycosyltransferases. All these datasets can be downloaded in the ‘Download’ module along with the statistics of the gene of interest or in each species.

RiceRelativesGD provides three online tools for the users: (1) the ‘Search’ tool performing sequence and keyword (gene name, GO term, Pfam and Swissprot term) search; (2) the ‘Genome browse’ tool providing a genome browser for each species. In the genome browser interfaces of ‘indica rice’ and ‘japonica rice’, genomic PAV of each rice relative and genes specific to indica or japonica rice can be observed, (3) the ‘Build phylogenetic tree’ tool providing an online tool for phylogenetic analysis. Users can input the gene names from the RiceRelativesGD database or paste sequences from external sources to construct a phylogenetic tree (Figure 1).

The overall analysis results

RiceRelativesGD incorporates publicly available genomic information from 12 species of the Oryzoideae subfamily (13,16–18,20,46) and the paddy weed E. crus-galli (7). A total of 34 570 gene families were identified, including 34 147 gene families with members from multiple species (multi-species families) and 423 gene families with members from only one species (species-specific families). Additionally, 91 085 orphan genes were found (Figure 2A). From the 11 rice relatives, a total of 208 321 genes were found to be specific to indica or japonica rice genome (Figure 2C).

Besides, the database also provides potentially valuable genes for rice breeding. We identified a total of 13 643 genes related to external stimulus responses (light, gravity, temperature, drought, salinity and abiotic stress) or photosynthesis (photophosphorylation, Calvin cycle, photorespiration and C4 photosynthesis) from the rice relatives (Figure 2B). The stress-related genes from rice relatives provide potential breeding targets for enhancing stress resistance of rice. Genes identified to be involved in the C4 photosynthesis pathways could be useful targets for improving the efficiency of rice production with less input of chemical fertilizers.

A case study for the application of RiceRelativesGD

The database provides detailed information for each gene, including organism of origin, genomic location, family ID, gene structure, function annotations and sequences. As an example EC_v6.g004014 is shown in Figure 3, which was defined as a gene encoding Cytochrome P450 based on Pfam domain, Swissprot and Mercator analyses. EC_v6.g004014 belongs to the OG0001327 gene family, which has 34 genes from 7 species. In the interface of the gene family OG0001327, users can find function annotations, overview of the gene family and the phylogenetic tree of all genes of the family as well as download sequences. The phylogenetic tree of the gene family includes not only family members but also the Pfam domain information of each gene allowing researchers to easily observe the differences and similarities among the family members. The phylogenetic tree can be downloaded as image or tree file by clicking ‘Download Tree’ (Figure 3).

Discussion

Rice relatives have become increasingly important for future improvement of rice varieties as they retain many competitive agronomic traits that have lost in rice during domestication and breeding with intensive artificial selection. Re-introducing these genetic elements back into the rice genetic background would not only enhance the performance of rice but also alleviate the increased genetic load caused by domestication and breeding (47). Genomic data of rice relatives are essential and crucial sources for uncovering the genes lost in the cultivated rice. Currently, most rice genomic databases do not provide information on rice relatives. Even though a few databases (e.g. Gramene, Ensemble Plants or PLAZA) integrate genomic resources of some rice relatives, they mainly focus on providing general information on the genomes, such as their orthologs and paralogs, gene gain/loss tree and genomic alignments, and pay no attention to a particular group of rice relatives or analysis of specific genes in rice relatives that could be valuable for modern rice breeding programs. RiceRelativesGD fills the gap by providing not only more comprehensive genomic information of rice relatives for the rice community but also specific genes from rice relatives, including stress-related genes, photosynthesis genes and so on. RiceRelativesGD collected and organized published genomic data of rice and its relatives from relevant literatures. Currently, a total of 208 321 specific genes from rice relatives are deposited in RiceRelativesGD.

In RiceRelativesGD, disagreement can be found between phylogenetic relationship (Fig. 1A) and specific gene number. Different genome assembly and gene prediction categories applied in each genome could affect the specific gene numbers. However, the difference of gene numbers or specific gene numbers may not be a convincing evidence to measure the phylogenetic distance among species. The protein-coding gene numbers were dynamically changed when the duplication, divergence or recombination events happened. For example, approximately 50 de novo genes generated and retained in recent divergence of Oryza per million years was validated (48). With the rapid development of sequencing technology and the declining sequencing cost, an increasing number of more complete genomes of rice relatives will be available, which could eliminate this disharmony. To keep pace with this trend, RiceRelativesGD should and will be updated to include all published genomic data of rice relatives. RiceRelativesGD is the first database focusing on providing comprehensive information on specific genes of diverse rice relatives. It will facilitate researchers mining genes with potential value for modern rice breeding. We will not only regularly update RiceRelativesGD with newly published data but also enhance the functionality with the aim of serving it as a foundation for future studies on rice relatives.

Authors’ contributions

L.F., L.M. and X.C. conceptualized the study. L.M. and M.C. conducted the study. L.M. and X.C. drafted the manuscript. L.F., Q.H.Z., X.C. and C.Y. supervised the study and revised the manuscript. Q.C., L.J., X.K. and D.W. were involved in data analysis and database construction. All authors read and approved the final manuscript.

Acknowledgements

National Natural Science Foundation of China (9143511/31401453), Natural Science Foundation of Zhejiang Province (LZ17C130001), Fundamental Research Funds for the Central Universities of China (2017QNA6013), 111 Project of China (B17039), China Postdoctoral Science Foundation (2017 M612002) and Jiangsu Collaborative Innovation Center for Modern Crop Production (JCIC-MCP).

Conflict of interest. None declared.

Database URL:http://ibi.zju.edu.cn/ricerelativesgd/

Reference

1.

Zhang
,
Q.-J.
,
Zhu
,
T.
,
Xia
,
E.-H.
et al.  (
2014
)
Rapid diversification of five Oryza AA genomes associated with rice adaptation
.
Proc. Natl. Acad. Sci.
,
111
,
E4954
E4962
.

2.

Sakai
,
H.
and
Itoh
,
T.
(
2010
)
Massive gene losses in Asian cultivated rice unveiled by comparative genome analysis
.
BMC Genomics
,
11
,
121
.

3.

Giuliani
,
R.
,
Koteyeva
,
N.
,
Voznesenskaya
,
E.
et al.  (
2013
)
Coordination of leaf photosynthesis, transpiration, and structural traits in rice and wild relatives (genus Oryza)
.
Plant Physiol.
,
162
,
1632
1651
.

4.

Jena
,
K.K.
(
2010
)
The species of the genus Oryza and transfer of useful genes from wild species into cultivated rice, O. sativa
.
Breed. Sci.
,
60
,
518
523
.

5.

Guo
,
L.
,
Qiu
,
J.
,
Li
,
L.F.
et al.  (
2018
)
Genomic clues for crop–weed interactions and evolution
.
Trends Plant Sci.
,
23
,
1102
1115
.

6.

Kellogg
,
E.A.
(
2009
)
The evolutionary history of Ehrhartoideae, Oryzeae, and Oryza
.
Rice
,
2
,
1
14
.

7.

Guo
,
L.
,
Qiu
,
J.
,
Ye
,
C.
et al.  (
2017
)
Echinochloa crus-galli genome analysis provides insight into its adaptation and invasiveness as a weed
.
Nat. Commun.
,
8
,
1031
.

8.

Huang
,
X.
,
Wei
,
X.
,
Sang
,
T.
et al.  (
2010
)
Genome-wide asociation studies of 14 agronomic traits in rice landraces
.
Nat. Genet.
,
42
,
961
967
.

9.

Huang
,
X.
,
Zhao
,
Y.
,
Wei
,
X.
et al.  (
2012
)
Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm
.
Nat. Genet.
,
44
,
32
39
.

10.

Huang
,
X.
,
Yang
,
S.
,
Gong
,
J.
et al.  (
2016
)
Genomic architecture of heterosis for yield traits in rice
.
Nature
,
537
,
629
633
.

11.

Yano
,
K.
,
Yamamoto
,
E.
,
Aya
,
K.
et al.  (
2016
)
Genome-wide association study using whole-genome sequencing rapidly identifies new genes influencing agronomic traits in rice
.
Nat. Genet.
,
48
,
927
934
.

12.

Zhao
,
Q.
,
Feng
,
Q.
,
Lu
,
H.
et al.  (
2018
)
Pan-genome analysis highlihgts the extent of genomic variation in cultivated and wild rice
.
Nat. Genet.
,
50
,
278
284
.

13.

Yu
,
J.
,
Hu
,
S.
,
Wang
,
J.
et al.  (
2002
)
A draft sequence of the rice genome (Oryza sativa L. ssp. indica)
.
Science
,
296
,
79
92
.

14.

Goff
,
S.A.
,
Ricke
,
D.
,
Lan
,
T.-H.
et al.  (
2002
)
A draft sequence of the rice genome (Oryza sativa L. ssp. japonica)
.
Science
,
296
,
92
100
.

15.

Kawahara
,
Y.
,
de la Bastide
,
M.
,
Hamilton
,
J.P.
et al.  (
2013
)
Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data
.
Rice
,
6
,
1
10
.

16.

Chen
,
J.
,
Huang
,
Q.
,
Gao
,
D.
et al.  (
2013
)
Whole-genome sequencing of Oryza brachyantha reveals mechanisms underlying Oryza genome evolution
.
Nat. Commun.
,
4
,
1595
.

17.

Wang
,
M.
,
Yu
,
Y.
,
Haberer
,
G.
et al.  (
2014
)
The genome sequence of African rice (Oryza glaberrima) and evidence for independent domestication
.
Nat. Genet.
,
46
,
982
988
.

18.

Guo
,
L.
,
Qiu
,
J.
,
Han
,
Z.
et al.  (
2015
)
A host plant genome (Zizania latifolia) after a century-long endophyte infection
.
Plant J.
,
83
,
600
609
.

19.

Du
,
H.
,
Yu
,
Y.
,
Ma
,
Y.
et al.  (
2017
)
Sequencing and de novo assembly of a near complete indica rice genome
.
Nat. Commun.
,
8
,
15324
.

20.

Stein
,
J.C.
,
Yu
,
Y.
,
Copetti
,
D.
et al.  (
2018
)
Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza
.
Nat. Genet.
,
50
,
285
296
.

21.

Kumar
,
S.
,
Stecher
,
G.
,
Suleski
,
M.
et al.  (
2017
)
TimeTree: a resource for timelines, Timetrees, and divergence times
.
Mol. Biol. Evol.
,
34
,
1812
1819
.

22.

Kawahara
,
Y.
,
de la Bastide
,
M.
,
Hamilton
,
J.P.
et al.  (
2013
)
Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data
.
Rice
,
6
,
4
.

23.

Cubry
,
P.
,
Tranchant-Dubreuil
,
C.
,
Thuillet
,
A.C.
et al.  (
2018
)
The rise and fall of African rice cultivation revealed by analysis of 246 new genomes
.
Curr. Biol.
,
28
,
2274
2282.e6
.

24.

Sakai
,
H.
,
Lee
,
S.S.
,
Tanaka
,
T.
et al.  (
2013
)
Rice annotation project database (RAP-DB): an integrative and interactive database for rice genomics
.
Plant Cell Physiol.
,
54
,
e6
.

25.

Song
,
J.M.
,
Lei
,
Y.
,
Shu
,
C.C.
et al.  (
2018
)
Rice information GateWay: a comprehensive bioinformatics platform for Indica rice genomes
.
Mol. Plant
,
11
,
505
507
.

26.

Zhao
,
W.
(
2004
)
BGI-RIS: an integrated information resource and comparative analysis workbench for rice genomics
.
Nucleic Acids Res.
,
32
,
D377
D382
.

27.

Sun
,
C.
,
Hu
,
Z.
,
Zheng
,
T.
et al.  (
2017
)
RPAN: rice pan-genome browser for ∼3000 rice genomes
.
Nucleic Acids Res.
,
45
,
597
605
.

28.

Mansueto
,
L.
,
Fuentes
,
R.R.
,
Chebotarov
,
D.
et al.  (
2016
)
SNP-seek II: a resource for allele mining and analysis of big genomic data in Oryza sativa
.
Curr. Plant Biol.
,
7–8
,
16
25
.

29.

Zhao
,
H.
,
Yao
,
W.
,
Ouyang
,
Y.
et al.  (
2015
)
RiceVarMap: a comprehensive database of rice genomic variations
.
Nucleic Acids Res.
,
43
,
D1018
D1022
.

30.

Ohyanagi
,
H.
,
Ebata
,
T.
,
Huang
,
X.
et al.  (
2016
)
OryzaGenome: genome diversity database of wild Oryza species
.
Plant Cell Physiol.
,
57
,
e1
.

31.

Yamazaki
,
Y.
,
Sakaniwa
,
S.
,
Tsuchiya
,
R.
et al.  (
2010
)
Oryzabase: an integrated information resource for rice science
.
Breed. Sci.
,
60
,
544
548
.

32.

Zhang
,
Z.
,
Hu
,
S.
,
He
,
H.
et al.  (
2016
)
Information commons for Rice (IC4R)
.
Nucleic Acids Res.
,
44
,
D1172
D1180
.

33.

Tello-Ruiz
,
M.K.
,
Naithani
,
S.
,
Stein
,
J.C.
et al.  (
2018
)
Gramene 2018: unifying comparative genomics and pathway resources for plant research
.
Nucleic Acids Res.
,
46
,
D1181
D1189
.

34.

Emms
,
D.M.
and
Kelly
,
S.
(
2015
)
OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy
.
Genome Biol.
,
16
,
157
.

35.

Jones
,
P.
,
Binns
,
D.
,
Chang
,
H.Y.
et al.  (
2014
)
InterProScan 5: genome-scale protein function classification
.
Bioinformatics
,
30
,
1236
1240
.

36.

Lohse
,
M.
,
Nagel
,
A.
,
Herter
,
T.
et al.  (
2014
)
Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data
.
Plant Cell Environ.
,
37
,
1250
1258
.

37.

Bateman
,
A.
,
Martin
,
M.J.
,
O’Donovan
,
C.
et al.  (
2015
)
UniProt: a hub for protein information
.
Nucleic Acids Res.
,
43
,
D204
D212
.

38.

Camacho
,
C.
,
Coulouris
,
G.
,
Avagyan
,
V.
et al.  (
2009
)
BLAST+: architecture and applications
.
BMC Bioinformatics
,
10
,
421
.

39.

Klopfenstein
,
D.V.
,
Zhang
,
L.
,
Pedersen
,
B.S.
et al.  (
2018
)
GOATOOLS: a python library for gene ontology analyses
.
Sci. Rep.
,
8
,
10872
.

40.

Katoh
,
K.
and
Standley
,
D.M.
(
2013
)
MAFFT multiple sequence alignment software version 7: improvements in performance and usability
.
Mol. Biol. Evol.
,
30
,
772
780
.

41.

Price
,
M.N.
,
Dehal
,
P.S.
and
Arkin
,
A.P.
(
2010
)
FastTree 2 - approximately maximum-likelihood trees for large alignments
.
PLoS One
,
5
,
e9490
.

42.

Letunic
,
I.
and
Bork
,
P.
(
2007
)
Interactive tree of life (iTOL): an online tool for phylogenetic tree display and annotation
.
Bioinformatics
,
23
,
127
128
.

43.

Kurtz
,
S.
,
Phillippy
,
A.
,
Delcher
,
A.L.
et al.  (
2004
)
Versatile and open software for comparing large genomes
.
Genome Biol.
,
5
,
R12
.

44.

Skinner
,
M.E.
,
Uzilov
,
A.V.
,
Stein
,
L.D.
et al.  (
2009
)
JBrowse: a next-generation genome browser
.
Genome Res.
,
19
,
1630
1638
.

45.

Kreft
,
L.
,
Botzki
,
A.
,
Coppens
,
F.
et al.  (
2017
)
PhyD3: a phylogenetic tree viewer with extended phyloXML support for functional genomics data visualization
.
Bioinformatics
,
33
,
2946
2947
.

46.

Matsumoto
,
T.
,
Wu
,
J.
,
Kanamori
,
H.
et al.  (
2005
)
The map-based sequence of the rice genome
.
Nature
,
436
,
793
800
.

47.

Liu
,
Q.
,
Zhou
,
Y.
,
Morrell
,
P.L.
et al.  (
2017
)
Deleterious variants in Asian Rice and the potential cost of domestication
.
Mol. Biol. Evol.
,
34
,
908
924
.

48.

Zhang
,
L.
,
Ren
,
Y.
,
Yang
,
T.
et al.  (
2019
)
Rapid evolution of protein diversity by de novo origination in Oryza
.
Nat. Ecol. Evol.
,
3
,
679
690
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.