Genome sequence of the model rice variety KitaakeX

doi:10.21203/rs.2.10528/v2

Download PDF

Research article

Genome sequence of the model rice variety KitaakeX

https://doi.org/10.21203/rs.2.10528/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 27 Nov, 2019

Read the published version in BMC Genomics →

You are reading this older preprint version

Read the latest preprint version →

Background: The availability of thousands of complete rice genome sequences from diverse varieties and accessions has laid the foundation for in-depth exploration of the rice genome. One drawback to these collections is that most of these rice varieties have long life cycles, and/or low transformation efficiencies, which limits their usefulness as model organisms for functional genomics studies. In contrast, the rice variety Kitaake has a rapid life cycle (9 weeks seed to seed) and is easy to propagate. For these reasons, Kitaake has emerged as a model for studies of diverse monocotyledonous species. Results: Here, we report the de novo genome sequencing and analysis of Oryza sativa ssp. japonica variety KitaakeX, a Kitaake plant carrying the rice XA21 immune receptor. Our KitaakeX sequence assembly contains 377.6 Mb, consisting of 33 scaffolds (476 contigs) with a contig N50 of 1.4 Mb. Complementing the assembly are detailed gene annotations of 35,594 protein coding genes. We identified 331,335 genomic variations between KitaakeX and Nipponbare (ssp. japonica), and 2,785,991 variations between KitaakeX and Zhenshan97 (ssp. indica). We also compared Kitaake resequencing reads to the KitaakeX assembly and identified 219 small variations. The high-quality genome of the model rice plant KitaakeX will accelerate rice functional genomics. Conclusions: The high quality, de novo assembly of the KitaakeX genome will serve as a useful reference genome for rice and will accelerate functional genomics studies of rice and other species.

Plant Molecular Biology and Genetics

Epigenetics & Genomics

Rice

Kitaake

KitaakeX

XA21 immune receptor

Whole genome sequence

De novo genome assembly

Nipponbare

Zhenshan97

Rice (Oryza sativa) provides food for more than half of the world’s population [1] and also serves as a model for studies of other monocotyledonous species. Cultivated rice contains two major types of O. sativa, the O. sativa indica/Xian group and the O. sativa japonica/Geng group. Using genomic markers, two additional minor types have been recognized, the circum-Aus group and the circum-Basmati group [2]. More than 3,000 rice varieties and species have been sequenced, including Nipponbare [3], 93-11 [4], DJ 123, IR64 [5], Zhenshan97, Minghui 63 [6], Shuhui498 [7], Oryza glaberrima [8] [2]. The availability of these genomes has laid a strong foundation for basic rice research and breeding [2]. However, the use of these sequenced varieties for functional genomics analyses is limited by their long life cycles or low transformation efficiencies. For example, it takes up to 6 months for Nipponbare to produce seeds under winter conditions. The Indica varieties typically have relatively low transformation efficiencies [9].

The Kitaake cultivar (ssp. japonica), which originated at the northern limit of rice cultivation in Hokkaido, Japan [10], has emerged as a model for rice research [11]. Kitaake is insensitive to day length, easy to propagate, relatively cold tolerant, short in stature and completes its life cycle in about 9 weeks [12] [13]. These properties make it easy to cultivate under typical greenhouse conditions. Kitaake is also highly amenable to transformation [14]. Several hundred genes have been overexpressed or silenced in KitaakeX [14]. The transformation efficiency of Kitaake is comparable to that of that Dongjin, a cultivar that historically transforms well [12]. Kitaake has been used to establish multiple mutant populations, including an RNAi mutant collection [15], T-DNA insertion collections [11], [16], and a whole-genome sequenced mutant population of KitaakeX, a Kitaake variety carrying the Xa21 immune receptor gene (formerly called X.Kitaake) [17, 18]. Kitaake has been used to explore diverse aspects of rice biology, including flowering time [19], disease resistance [20], [21], [22], small RNA biology [23], and the CRISPR-Cas9 and TALEN technologies [24], [25].

The unavailability of the Kitaake genome sequence has posed an obstacle to the use of Kitaake in rice research. For example, analysis of a fast-neutron (FN) induced mutant population in KitaakeX, a Kitaake plant carrying the rice XA21 gene [17], required the use of Nipponbare (ssp. japonica) as the reference genome. Additionally, CRISPR/Cas9 guide RNAs cannot be accurately designed for Kitaake without a complete sequence. To address these issues, we assembled a high-quality genome sequence of KitaakeX, compared its genome to the genomes of rice varieties Nipponbare and Zhenshan97 (ssp. indica), and identified genomic variations. The XA21 gene confers resistance to the bacterial pathogen, Xanthomonas oryzae pv. oryzae, making KItaakeX a model for studies of infectious disease [18].

KitaakeX flowers significantly earlier than other sequenced rice varieties

Kitaake has long been recognized as a rapid life-cycle variety [26], but it has yet to be systematically compared to other rice varieties. We compared the flowering time of KitaakeX with other sequenced rice varieties under long-day conditions (14 h light/10 h dark). Consistent with other studies, we found that KitaakeX flowers much earlier than other varieties (Fig. 1a, 1b), heading at 54 days after germination. Other rice varieties Nipponbare, 93-11 (ssp. indica), IR64 (ssp. indica), Zhenshan 97, Minghui 63 (ssp. indica), and Kasalath (aus rice cultivar) start heading at 134, 99, 107, 79, 125, and 84 days after germination, respectively (Fig. 1b).

We next assessed how KitaakeX is related to other rice varieties using a phylogenetic approach based on the rice population structure and diversity published for 3,010 varieties [2]. The 3,010 sequenced accessions were classified into nine subpopulations, most of which could be connected to geographical origins. The phylogenetic tree reveals that KitaakeX and Nipponbare are closely related within the same subpopulation (Fig. 1c).

Genome sequencing and assembly

To obtain a high-quality, de novo genome assembly, we sequenced the KitaakeX genome using a strategy that combines short-read and long-read sequencing. Sequencing reads were collected using Illumina, 10x Genomics, PACBIO, and Sanger platforms at the Joint Genome Institute (JGI) and the HudsonAlpha Institute. The current release is version 3.0, which is a combination of a MECAT (Mapping, Error Correction and de novo Assembly Tools) PACBIO based assembly and an Illumina sequenced 10x genomics SuperNova assembly. The assembled sequence contains 377.6 Mb, consisting of 33 scaffolds (476 contigs) with a contig N50 of 1.4 Mb, covering a total of 99.67% of assembled bases in chromosomes (Table 1).

We assessed the quality of the KitaakeX assembly for sequence completeness and accuracy. Completeness of the assembly was assessed by aligning the 34,651 annotated genes from the v7.0 Nipponbare to the KitaakeX assembly using BLAT [27]. The alignments indicate that 98.94% (34,285 of genes) genes completely aligned to the KitaakeX assembly, 0.75% (259 genes) partially aligned, and 0.31% (107 genes) were not detected. A bacterial artificial chromosome (BAC) library was constructed and a set of 346 BAC clones (9.2x clone coverage) was sequenced using PACBIO sequencing. A range of variants was detected by comparing the BAC clones to the assembly. Alignments were of high quality (<0.1% of error) in 271 clones (Additional file 1: Figure S13). Sixty BACs indicate a higher error rate (0.45% of error) due mainly to their placement in repetitive regions (Additional file 1: figure S14). Fifteen BAC clones indicate a rearrangement (10 clones) or a putative overlap on adjacent contigs (5 clones) (Additional file 1: figure S15). The overall error rate in the BAC clones is 0.09%, indicating the high quality of this assembly (for detailed information, see Additional file 1).

Genome annotation

We predicted 35,594 protein-coding genes in the KitaakeX genome (Table 1, Additional file 2: Table S12), representing 31.5% genic space of the assembled genome size (Table 1). There is some transcriptome support for 89.5% (31,854/35,594) of the KitaakeX genes, and 81.6% (29,039/35,594) genes are fully supported by the transcriptome (Additional file 2: Table S11). The predicted protein-coding genes are distributed unevenly across each chromosome; gene density tends to be higher toward chromosome ends (Fig. 2f). The average GC content of the genome is 43.7% (Fig. 2e, Table 1).

To assess the quality of the annotation of KitaakeX genes, we compared the KitaakeX annotation to those of other completed rice genomes using the BUSCO v2 method, which is based on a set of 1,440 conserved plant genes. The results confirm 99.0% completeness of the KitaakeX genome annotation (Table1, Additional file 2: Table S7). To further evaluate the quality of the annotation, we studied the extent of conservation of functional genes in KitaakeX. We selected 291 genes (Additional file 3: Table S13) from three pathways associated with stress resistance, flowering time and response to light [8], and then searched for orthologous genes in the KitaakeX genome. We found that 275 of 291 (94.5%) of the selected KitaakeX genes show greater than 90% identity with the corresponding Nipponbare genes at the protein level. Twenty-three out of the 291 show 100% identity on the genome level but not on the protein level. Of these 23 genes, the KitaakeX gene model for 16 genes has better transcriptomic evidence than the Nipponbare gene model. One of the 291 KitaakeX genes is slightly shorter than its Nipponbare ortholog due to an alternative transcript (Additional file 3: Table S13). These results indicate the high quality of the annotation, and conservation between the KitaakeX and Nipponbare japonica rice varieties.

Using SynMap, we identified 2,469 pairs of colinear genes (88 blocks) in the KitaakeX genome (Fig. 2g). These results correlate with already published findings [28]. We used RepeatMaker and Blaster to identify transposable elements (TEs) in the KitaakeX genome, and identified 122.2 Mb of sequence corresponding to TEs (32.0% of the genome). DNA transposons account for ~33 Mb; retrotransposons account for ~90 Mb. The TEs belong mostly to the Gypsy and Copia retroelement families, and account for 23% of the genome (Additional file 2: Table S8), as is true in the Nipponbare and Zhenshan97 genomes [6].

Table 1 Summary of the KitaakeX genome assembly and annotation

	Estimated genome size	409.5 Mb
	Assembled contigs	377.6 Mb
	Contig N50	1.4 Mb
Genome assembly	Longest contig	8.6 Mb
	Assembled scaffolds	381.6 Mb
	Scaffold N50	30.3 Mb
	Longest scaffold	44.3 Mb
	GC content	43.7%
	Retrotransposons	89.6 Mb
Transposable elements	DNA transposons	32.6 Mb
	Total	122.2 Mb
	Protein-coding genes	35,594
	Complete BUSCOs	99.0%
Genome annotation	Average transcript length	1,874 bp
	Average coding sequence length	1,222 bp
	Functionally annotated	33,226

Genomic variations between KitaakeX and other rice varieties

We compared the genome of KitaakeX to the Nipponbare and Zhenshan97 genomes to detect genomic variations, including single nucleotide polymorphisms (SNPs), insertions and deletions under 30 bp (InDels), presence/absence variations (PAVs), and inversions using MUMmer [29]. We found 331,335 variations between KitaakeX and Nipponbare (Additional file 4), and nearly 10 times as many (2,785,991) variations between KitaakeX and Zhenshan97 (Additional file 5). There are 253,295 SNPs and 75,183 InDels between KitaakeX and Nipponbare, and 2,328,319 SNPs and 442,962 InDels between KitaakeX and Zhenshan97 (Additional files 6 and Additional file 2: Table S3). With respect to SNPs in both intersubspecies (japonica vs. indica) as well as intrasubspecies (japonica vs. japonica) comparisons, transitions (Tss) (G ->A and C ->T) are about twice as abundant as transversions (Tvs) (G ->C and C ->G) (Additional file 2: Table S10). Genomic variations between KitaakeX and Nipponbare are highly concentrated in some genomic regions (Fig. 2b), but variations between KitaakeX and Zhenshan97 are spread evenly through the genome (Fig. 2c). Intersubspecies genomic variations, then, are much more extensive than intrasubspecies variations. We also detected multiple genomic inversions using comparative genomics (Additional files 4 and 5).

For variations occurring in the genic regions, we found that single-base and 3 bp (without frame shift) InDels are much more abundant than others (Additional file 7: Figure S16a), suggesting that these genetic variations have been functionally selected. We carried out detailed analysis of gene structure alterations that exist as a consequence of SNPs and InDels between KitaakeX and Nipponbare and Kitaake and Zhenshan97. Between KitaakeX and Nipponbare, we identified 2,092 frameshifts, 78 changes affecting splice-site acceptors, 71 changes affecting splice-site donors, 19 lost start codons, 161 gained stop codons, and 15 lost stop codons. In the comparison of KitaakeX to Zhenshan97, 6,809 unique genes in KitaakeX are affected by 8,640 frameshifts (Additional file 7: Figure S16b), 531 changes affecting splice-site acceptors, 530 changes affecting splice-site donors, 185 lost start codons, 902 gained stop codons and 269 lost stop codons (Additional file 7: Figure S16b).

Based on PAV analysis, we identified 456 loci that are specific to KitaakeX (Additional file 4) compared with Nipponbare. Pfam analysis of KitaakeX-specific regions revealed 275 proteins. Out of these 275 genes, 148 genes are from 19 different gene families with more than 2 genes in those regions. These gene families include protein kinases, leucine-rich repeat proteins, NB-ARC domain-containing proteins, F-box domain containing proteins, protein tyrosine kinases, Myb/SANt-like DNA binding domain proteins, transferase family proteins, xylanase inhibitor C-terminal protein, and plant proteins of unknown function (Additional file 7: Figure S16c). We identified 4589 loci specific to KitaakeX compared with Zhenshan97 (Additional file 5).

We also compared our de novo assembly of KitaakeX genome with Kitaake resequencing reads using an established pipeline [17] . This analysis revealed 219 small variations (200 SNPs and 19 INDELs) between the two genomes (Additional file 8). These variations affect 9 genes in KitaakeX besides the Ubi-Xa21 transgene, including the selectable marker encoding a hygromycin B phosphotransferase on chromosome 6 (Additional file 8, Additional file 9: Figure S17).

In 2005 the Nipponbare genome was sequenced and annotated to a high-quality level (International Rice Genome Sequencing and Sasaki 2005). Since that time, it has served as a reference genome for many rice genomic studies [30]. Despite its use, the long life cycle of Nipponbare makes it time-consuming for most genetic analyses.

Here we report the de novo assembly and annotation of KitaakeX, an early-flowering rice variety with a rapid life cycle that is easy to propagate under greenhouse conditions. We predict that KitaakeX contains 35,594 protein-coding genes, comparable to the published genomes (39,045 for Nipponbare and 34,610 for Zhenshan97) (Additional file 4 and Additional file 5). The availability of a high-quality genome and annotation for KitaakeX will be useful for associating traits of interest with genetic variations, and for identifying the genes controlling those traits.

We identified 219 SNPs and InDels between the KitaakeX and Kitaake genomes. These variations may have resulted from somatic mutations that arose during tissue culture and regeneration, or they may be spontaneous mutations [31]. For rice, 150 mutations are typically induced during tissue culture and 41 mutations occur spontaneously per three generations [31]. These numbers are consistent with the independent propagation of KitaakeX and Kitaake over approximately 10 generations in the greenhouse.

The KitaakeX genome will be useful for variety of studies. For example, we recently published the whole genome sequences of 1,504 FN-mutated KitaakeX rice lines. Mutations were identified by aligning reads of the KitaakeX mutants to the Nipponbare reference genome [17]. On average, 97% of the Nipponbare genome is covered by the KitaakeX reads. However, in some regions, the KitaakeX genome diverges from Nipponbare to such an extent that no variants can be confidently identified. These appear either as gaps in coverage or as regions containing a concentration of natural variations between KitaakeX and Nipponbare. We can now use the KitaakeX sequence as the direct reference genome and detect mutations in highly variable regions. This approach will simplify analysis and increase confidence in the identification of FN-induced mutations. Because there are only 219 small variations between KitaakeX and Kitaake (Additional file 8), the KitaakeX genome can also be used as the reference genome for Kitaake.

Conclusions:

The de novo assembly of the KitaakeX genome serves as a useful reference genome for the model rice variety Kitaake and will facilitate investigations into the genetic basis of diverse traits critical for rice biology and genetic improvement.

Plant materials and growth conditions

Dr. Thomas W. Okita from Washington State University provided the Kitaake seeds, which were originally obtained from Dr. Hiroyuki Ito, Akita National College of Technology, Japan. Dr. Jan E. Leach at Colorado State University provided seeds for Zhenshan 97, Minghui 63, IR64 and 93-11. Seeds of Kasalath were provided by the USDA Dale Bumpers National Rice Research Center, Stuttgart, Arkansas. Seeds were germinated on 1/2x MS (Murashige and Skoog) medium. Seedlings were transferred to a greenhouse and planted 3 plants/pot during the springtime (Mar. 2, 2017) in Davis, California. The light intensity was set at approximately 250 μmol m^− 2 s^− 1. The day/night period was set to 14/10 h, and the temperature was set between 28 and 30 °C [32]. Rice plants were grown in sandy soil supplemented with nutrient water. The day when the first panicle of the plant emerged was recorded as the heading date for that plant. Kasalath seeds were received later, and the heading date was recorded in the same way. The experiment was repeated in winter.

Construction of a phylogenetic tree

We obtained 178,496 evenly distributed SNPs by dividing the genome into 3.8 kb bins and selecting one or two SNPs per bin randomly according to the SNP density of the bin. Genotypes of all the rice accessions, including 3,010 accessions of the 3K Rice Genomes Project and additional noted accessions, were fetched from the SNP database RiceVarMap v2.0 [33] and related genomic data [34] and used to calculate an IBS distance matrix which was then applied to construct a phylogenetic tree by the unweighted neighbor-joining method, implemented in the R package APE [35]. Branches of the phylogenetic tree were colored according to the classification of the 3,010 rice accessions [2].

Genome Sequencing and Assembly

High molecular weight DNA from young leaves of KitaakeX was isolated and used in sequencing.

See (Additional file 1) for further details.

Annotation of Protein-Coding Genes

To obtain high-quality annotations, we performed high throughput RNA-seq analysis of libraries from diverse rice tissues (leaf, stem, panicle, and root). Approximately 683 million pairs of 2x151 paired-end RNA-seq reads were obtained and assembled using a comprehensive pipeline PERTRAN (unpublished). Gene models were predicted by combining ab initio gene prediction, protein-based homology searches, experimentally cloned cDNAs/expressed-sequence tags (ESTs) and assembled transcripts from the RNA-seq data. Gene functions were further annotated according to the best-matched proteins from the SwissProt and TrEMBL databases [36] using BLASTP (E value < 10^-5) (Additional file 11). Genes without hits in these databases were annotated as “hypothetical proteins”. Gene Ontology (GO) [37] term assignments and protein domains and motifs were extracted with InterPro [38]. Pathway analysis was derived from the best-match eukaryotic protein in the Kyoto encyclopedia of genes and genomes (KEGG) database [39] using BLASTP (E value<1.0e^-10).

Genome Synteny

We used SynMap (CoGe, www.genomevolution.org) to identify collinearity blocks using homologous CDS pairs with parameters according to Daccord et al [40] and visualized collinearity blocks using Circos [41].

Repeat Annotation

The fraction of transposable elements and repeated sequences in the assembly was obtained merging the output of RepeatMasker (http://www.repeatmasker.org/, v. 3.3.0) and Blaster (a component of the REPET package) [42]. The two programs were run using nucleotide libraries (PReDa and RepeatExplorer) from RiTE-db [43] and an in-house curated collection of transposable element (TE) proteins, respectively. Reconciliation of masked repeats was carried out using custom Perl scripts and formatted in gff3 files. Infernal [44] was adopted to identify non-coding RNAs (ncRNAs) using the Rfam library Rfam.cm.12.2 [45]. Results with scores lower than the family-specific gathering threshold were removed; when loci on both strands were predicted, only the hit with the highest score was kept. Transfer RNAs were also predicted using tRNAscan-SE [46] at default parameters. Repeat density was calculated from the file that contains the reconciled annotation (Additional file 10).

Analysis of Genomic Variations

Analysis of SNPs and InDels: We used MUMmer (version 3.23) [29] to align the Nipponbare and Zhenshan97 genomes to the KitaakeX genome using parameters -maxmatch -c 90 -l 40. To filter the alignment results, we used the delta -filter -1 parameter with the one-to-one alignment block option. To identify SNPs and InDels we used show-snp option with parameter (-Clr TH). We used snpEff [47] to annotate the effects of SNPs and InDels. Distribution of SNPs and InDels along the KitaakeX genome was visualized using Circos [41].

Analysis of PAVs and Inversions: We used the show-coords option of MUMmer (version 3.23) with parameters -TrHcl to identify gap regions and PAVs above 86 bp in size from the alignment blocks. We used the inverted alignment blocks with ≥98% identity from the show-coords output file to identify inversions.

To identify genomic variations between Kitaake and KitaakeX we sequenced and compared the sequences using the established pipeline [48].

BAC library construction

Arrayed BAC libraries were constructed using established protocols [49]. Please see Additional file 1 for further details.

Genome size estimation

We used the following methodology to estimate KitaakeX genome size:

(1) Using the Illumina fragment library, we created a histogram of 24mer frequencies. This was performed by first counting the frequency of all 24mers. The number of kmers at each frequency was tallied, and a histogram was created. (2) The kmer histogram generally indicates a peak value at a particular frequency corresponding to the average coverage of 24mers on the genome. (3) We then took the peak value representing the coverage on the genome, and computed the total bases in the Illumina library. Further dividing the total bases by the coverage, provided an estimate of the genome size. This value is generally accurate to +/- 10%.

BAC: Bacterial Artificial Chromosome

FN: Fast Neutron

BLAST: Basic Local Alignment Search Tool

BLAT: BLAST-like alignment tool

MECAT: Mapping, Error Correction and de novo Assembly Tools

BUSCO: Benchmarking Universal Single-Copy Orthologs

TEs: Transposable Elements

SNP: Single Nucleotide Polymorphisms

PAVs: Presence/Absence Variations

GO: Gene Ontology

KEGG: Kyoto encyclopedia of genes and genomes

NB-ARC: Nucleotide-Binding Adaptor shared by APAF-1, R proteins, and CED-4

EST: Expressed-Sequence Tags

MS: Murashige and Skoog

Ethics approval and consent to participate:

Not applicable

Consent to Publication:

Not applicable

Availability of data and material:

The genome sequencing reads and assembly have been deposited under GenBank under accession number PRJNA234782 and PRJNA448171 respectively. The assembly and annotation of the Kitaake genome are available at Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html). The RNA-Seq reads of KitaakeX leaf, panicle, stem and root have been deposited under GenBank accession numbers SRP182736, SRP182738, SRP182741, and SRP182737 respectively. Genome sequencing reads for Kitaake have been deposited under GenBank under accession number SRP193308.

Competing interests:

The authors declare no conflict of interest.

Funding:

This work was supported by NIH (GM59962) NIH (GM122968) and NSF (IOS-1237975) grants to PCR. It was also supported in part by the U. S. Department of Energy, Office of Science, Office of Biological and Environmental Research, through contract DE-AC02-05CH11231 between Lawrence Berkeley National Laboratory and the U. S. Department of Energy. The work conducted by the US Department of Energy Joint Genome Institute (JGI) was supported by the Office of Science of the US Department of Energy under Contract no. DE-AC02-05CH11231.

Author’s Contributions

R.J., G.L., M.C. and P.C.R. conceived and initiated the study. R.J. and G.L. carried out sequencing, assembly and annotation in collaboration with J.J., S.S., D.A.K., J.T., D.D., D.B., J.G., D.C., K.W.B., R.A.W., N.T.P., and J.S. R.J. and G.L. conducted comparative genomics analyses. F.X. and W.X. contributed to the phylogenetic tree. R.J., P.C.R., G.L and C.R.N. wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Rick A. Rios, Maria E. Hernandez, and Natasha Brown for assistance in genomic DNA isolation and submission and seed organization.

Gross BL, Zhao Z: Archaeological and genetic insights into the origins of domesticated rice. 2014, 111(17):6190-6197.
Wang W, Mauleon R, Hu Z, Chebotarov D, Tai S, Wu Z, Li M, Zheng T, Fuentes RR, Zhang FJN: Genomic variation in 3,010 diverse accessions of Asian cultivated rice. 2018, 557(7703):43.
Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, Schwartz DC, Tanaka T, Wu J, Zhou S et al: Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice (N Y) 2013, 6(1):4.
Yu J, Hu S, Wang J, Wong GK-S, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang XJs: A draft sequence of the rice genome (Oryza sativa L. ssp. indica). 2002, 296(5565):79-92.
Schatz MC, Maron LG, Stein JC, Wences AH, Gurtowski J, Biggers E, Lee H, Kramer M, Antoniou E, Ghiban EJGb: Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica. 2014, 15(11):506.
Zhang J, Chen L-L, Xing F, Kudrna DA, Yao W, Copetti D, Mu T, Li W, Song J-M, Xie W et al: Extensive sequence divergence between the reference genomes of two elite <em>indica</em> rice varieties Zhenshan 97 and Minghui 63. 2016, 113(35):E5163-E5171.
Du H, Yu Y, Ma Y, Gao Q, Cao Y, Chen Z, Ma B, Qi M, Li Y, Zhao XJNc: Sequencing and de novo assembly of a near complete indica rice genome. 2017, 8:15324.
Wang M, Yu Y, Haberer G, Marri PR, Fan C, Goicoechea JL, Zuccolo A, Song X, Kudrna D, Ammiraju JS et al: The genome sequence of African rice (Oryza glaberrima) and evidence for independent domestication. Nat Genet 2014, 46(9):982-988.
Tie W, Zhou F, Wang L, Xie W, Chen H, Li X, Lin Y: Reasons for lower transformation efficiency in indica rice using Agrobacterium tumefaciens-mediated transformation: lessons from transformation assays and genome-wide expression profiling. Plant molecular biology 2012, 78(1-2):1-18.
lchitani K, Okumoto Y, Tanisaka T: Photoperiod sensitivity gene of se-1 locus found in photoperiod insensitive rice cultivars of the northern limit region of rice cultivation. Breeding science 1997, 47:8.
Kim SL, Choi M, Jung KH, An G: Analysis of the early-flowering mechanisms and generation of T-DNA tagging lines in Kitaake, a model rice cultivar. J Exp Bot 2013, 64(14):4169-4182.
Kim SL, Choi M, Jung K-H, An GJJoeb: Analysis of the early-flowering mechanisms and generation of T-DNA tagging lines in Kitaake, a model rice cultivar. 2013, 64(14):4169-4182.
Kunihiro Y, Ebe Y, Wada S, Shinbashi N, Honma A, Sasaki T, Sasaki K, Numao Y, Morimura K, Tan No H: The new rice variety Kita-ake. Bulletin of Hokkaido prefectural agricultural experiment stations 1989, 59:4.
Jung KH, An G, Ronald PC: Towards a better bowl of rice: assigning function to tens of thousands of rice genes. Nat Rev Genet 2008, 9(2):91-101.
Wang L, Zheng J, Luo Y, Xu T, Zhang Q, Zhang L, Xu M, Wan J, Wang MB, Zhang CJPbj: Construction of a genomewide RNA i mutant library in rice. 2013, 11(8):997-1005.
Gao H, Zheng XM, Fei G, Chen J, Jin M, Ren Y, Wu W, Zhou K, Sheng P, Zhou F et al: Ehd4 encodes a novel and Oryza-genus-specific regulator of photoperiodic flowering in rice. PLoS Genet 2013, 9(2):e1003281.
Li G, Jain R, Chern M, Pham NT, Martin JA, Wei T, Schackwitz WS, Lipzen AM, Duong PQ, Jones KC et al: The Sequences of 1504 Mutants in the Model Rice Variety Kitaake Facilitate Rapid Functional Genomic Studies. Plant Cell 2017, 29(6):1218-1231.
Song WY, Wang GL, Chen LL, Kim HS, Pi LY, Holsten T, Gardner J, Wang B, Zhai WX, Zhu LH et al: A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science 1995, 270(5243):1804-1806.
Gao H, Jin M, Zheng X-M, Chen J, Yuan D, Xin Y, Wang M, Huang D, Zhang Z, Zhou KJPotNAoS: Days to heading 7, a major quantitative locus determining photoperiod sensitivity and regional adaptation in rice. 2014, 111(46):16337-16342.
Ronald PC, Beutler B: Plant and animal sensors of conserved microbial signatures. Science 2010, 330(6007):1061-1064.
Liu Y, Wu H, Chen H, Liu Y, He J, Kang H, Sun Z, Pan G, Wang Q, Hu JJNb: A gene cluster encoding lectin receptor kinases confers broad-spectrum and durable insect resistance in rice. 2015, 33(3):301.
Zhou X, Liao H, Chern M, Yin J, Chen Y, Wang J, Zhu X, Chen Z, Yuan C, Zhao W et al: Loss of function of a rice TPR-domain RNA-binding protein confers broad-spectrum disease resistance. 2018, 115(12):3174-3179.
Rodrigues JA, Ruan R, Nishimura T, Sharma MK, Sharma R, Ronald PC, Fischer RL, Zilberman D: Imprinted expression of genes and small RNA is associated with localized hypomethylation of the maternal genome in rice endosperm. Proc Natl Acad Sci U S A 2013, 110(19):7934-7939.
Li T, Liu B, Spalding MH, Weeks DP, Yang B: High-efficiency TALEN-based gene editing produces disease-resistant rice. Nat Biotechnol 2012, 30(5):390-392.
Xie K, Minkenberg B, Yang Y: Boosting CRISPR/Cas9 multiplex editing capability with the endogenous tRNA-processing system. Proc Natl Acad Sci U S A 2015, 112(11):3570-3575.
Jung K-H, An G, Ronald PC: Towards a better bowl of rice: assigning function to tens of thousands of rice genes. Nature Reviews Genetics 2008, 9:91.
Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res 2002, 12(4):656-664.
Guyot R, Keller BJG: Ancestral genome duplication in rice. 2004, 47(3):610-614.
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SLJGb: Versatile and open software for comparing large genomes. 2004, 5(2):R12.
Matsumoto T, Wu J, Itoh T, Numa H, Antonio B, Sasaki T: The Nipponbare genome and the next-generation of rice genomics research in Japan. Rice (New York, NY) 2016, 9(1):33-33.
Tang X, Liu G, Zhou J, Ren Q, You Q, Tian L, Xin X, Zhong Z, Liu B, Zheng XJGb: A large-scale whole-genome sequencing analysis reveals highly specific genome editing by both Cas9 and Cpf1 (Cas12a) nucleases in rice. 2018, 19(1):84.
Schwessinger B, Bahar O, Thomas N, Holton N, Nekrasov V, Ruan D, Canlas PE, Daudi A, Petzold CJ, Singan VR et al: Transgenic expression of the dicotyledonous pattern recognition receptor EFR in rice leads to ligand-dependent activation of defense responses. PLoS Pathog 2015, 11(3):e1004809.
Zhao H, Yao W, Ouyang Y, Yang W, Wang G, Lian X, Xing Y, Chen L, Xie W: RiceVarMap: a comprehensive database of rice genomic variations. Nucleic Acids Res 2015, 43(Database issue):D1018-1022.
Li G, Chern M, Jain R, Martin JA, Schackwitz WS, Jiang L, Vega-Sanchez ME, Lipzen AM, Barry KW, Schmutz J et al: Genome-Wide Sequencing of 41 Rice (Oryza sativa L.) Mutated Lines Reveals Diverse Mutations Induced by Fast-Neutron Irradiation. Mol Plant 2016, 9(7):1078-1081.
Paradis E, Claude J, Strimmer K: APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics 2004, 20(2):289-290.
Bairoch A, Apweiler R: The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res 2000, 28(1):45-48.
Dutkowski J, Kramer M, Surma MA, Balakrishnan R, Cherry JM, Krogan NJ, Ideker T: A gene ontology inferred from molecular networks. Nat Biotechnol 2013, 31(1):38-45.
Finn RD, Attwood TK, Babbitt PC, Bateman A, Bork P, Bridge AJ, Chang HY, Dosztanyi Z, El-Gebali S, Fraser M et al: InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res 2017, 45(D1):D190-D199.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K: KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res 2017, 45(D1):D353-D361.
Daccord N, Celton JM, Linsmith G, Becker C, Choisne N, Schijlen E, van de Geest H, Bianco L, Micheletti D, Velasco R et al: High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat Genet 2017, 49(7):1099-1106.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res 2009, 19(9):1639-1645.
Flutre T, Duprat E, Feuillet C, Quesneville H: Considering transposable element diversification in de novo annotation approaches. PLoS One 2011, 6(1):e16526.
Copetti D, Zhang J, El Baidouri M, Gao D, Wang J, Barghini E, Cossu RM, Angelova A, Maldonado LC, Roffler S et al: RiTE database: a resource database for genus-wide rice genomics and evolutionary biology. BMC Genomics 2015, 16:538.
Nawrocki EP, Eddy SR: Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 2013, 29(22):2933-2935.
Nawrocki EP, Burge SW, Bateman A, Daub J, Eberhardt RY, Eddy SR, Floden EW, Gardner PP, Jones TA, Tate J et al: Rfam 12.0: updates to the RNA families database. Nucleic Acids Res 2015, 43(Database issue):D130-137.
Schattner P, Brooks AN, Lowe TM: The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res 2005, 33(Web Server issue):W686-689.
Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 2012, 6(2):80-92.
Li G, Jain R, Chern M, Pham NT, Martin JA, Wei T, Schackwitz WS, Lipzen AM, Duong PQ, Jones KCJTPC: The sequences of 1,504 mutants in the model rice variety Kitaake facilitate rapid functional genomic studies. 2017:tpc. 00154.02017.
Luo M, Wing RA: An Improved Method for Plant BAC Library Construction. In: Plant Functional Genomics. Edited by Grotewold E. Totowa, NJ: Humana Press; 2003: 3-19.

Additional file 1: Supplementary Methods: Sequencing, genome assembly and construction of pseudomolecule chromosomes, and BAC library construction. Table S1: Genomic libraries included in the Oryza sativa (KitaakeX) genome assembly and their respective assembled sequence coverage levels in the final release. Table S2: PACBIO library statistics for single pass yield of the 42 chips included in the Oryza sativa (KitaakeX) genome assembly and their respective assembled sequence coverage levels. Table S3: Summary statistics of the output of the SuperNova whole genome shotgun assembly prior to integration with the PACBIO assembly. Table S4: Summary statistics of the raw output of the MECAT whole genome shotgun assembly. Figure S1. Syntenic Japonica sequence placements on the Oryza sativa (var. KitaakeX) chromosome 1. Figure S2. Syntenic Japonica sequence placements on the Oryza sativa (var. KitaakeX) chromosome 2. Figure S3: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 3. Figure S4: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 4. Figure S5: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 5. Figure S6: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 6. Figure S7: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 7. Figure S8: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 8. Figure S9: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 9. Figure S10: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 10. Figure S11: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 11. Figure S12: Syntenic Japonica sequence placements on the Oryza sativa (KitaakeX) chromosome 12. Table S5: Final summary assembly statistics for chromosome scale assembly. Figure S13: Dot plot of BAC clone 119492 on a region of Chr_02. Figure S14: Dot plot of BAC clone 120743 on a region of Chr_12. Figure S15: Dot plot of BAC clone 119503 in a region of Chr_06. Table S6: Rice, KitaakeX, BAC libraries.

Additional file 2: Table S7: BUSCO analysis of KitaakeX and comparison with other rice genomes. Table S8: Summary of transposable elements in KitaakeX, Nipponbare, and Zhenshan97. Table S9: Comparison of SNPs and INDELs between three rice genomes. Table S10: Comparison of single base substitutions between three rice genomes. Table S11: Oryza sativa KitaakeX annotation v3.1 on assembly v3.0. Table S12: Sequence length of pseudomolecules, number of genes and gene models for each of the 12 rice chromosomes

Additional file 3: Table S13: Genes used in annotation quality control

Additional file 4: Comparative genomic analysis between KitaakeX and Nipponbare. SNPs, InDels, PAVs, Inversions, and genes affected by SNPs, IndDels, PAVs and Inversions are listed in this file.

Additional file 5: Comparative genomic analysis between KitaakeX and Zhenshan97. SNPs, InDels, PAVs, Inversions, and genes affected by SNPs, IndDels, PAVs and Inversions are listed in this file.

Additional file 6: SNPs between KitaakeX and Zhenshan97

Additional file 7: Figure S16: Genomic variation showing gene variations between KitaakeX and Nipponbare and ZS97. a. Length distribution of InDels in protein-coding regions. b. SNPs and InDels that cause high-impact gene variations between KitaakeX and Nipponbare and ZS97. c. Gene enrichment in KitaakeX unique present regions compared with Nipponbare.

Additional file 8. Genomic variations between KitaakeX and Kitaake

Additional file 9: Figure S17: Integrative genomics viewer (IGV) snapshot showing presence of XA21 transgene and selectable marker encoding a hygromycin B phosphotransferase on chromosome 6 of KitaakeX.

Additional file 10: Repeat annotation of KitaakeX genome

Additional file 11: Functional annotation of KitaakeX genome

Download PDF

Journal Publication

published 27 Nov, 2019

Read the published version in BMC Genomics →

Editorial decision: Minor revision
19 Oct, 2019
Review #4 received at journal
16 Oct, 2019
Reviewer #4 agreed at journal
02 Oct, 2019
Reviewers invited by journal
29 Sep, 2019
Reviewer #1 agreed at journal
29 Sep, 2019
Reviewer #2 agreed at journal
29 Sep, 2019
Reviewer #3 agreed at journal
29 Sep, 2019
Review #1 received at journal
29 Sep, 2019
Review #2 received at journal
29 Sep, 2019
Review #3 received at journal
29 Sep, 2019
Editor assigned by journal
23 Sep, 2019
Submission checks completed at journal
22 Sep, 2019
Editor invited by journal
22 Sep, 2019

You are reading this older preprint version

Read the latest preprint version →

Genome sequence of the model rice variety KitaakeX

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Results

Discussion

Methods

Abbreviations

Declarations

References

Additional file legend

Supplementary Files

Status:

Journal Publication

Version 2