Abstract:HomeMolecular Plant-Microbe Interactions®Vol. 35, No. 8Genome Sequence of Magnaporthe oryzae EA18 Virulent to Multiple Widely Used Rice Varieties PreviousNext RESOURCE ANNOUNCEMENT OPENOpen Access licenseGenome Sequence of Magnaporthe oryzae EA18 Virulent to Multiple Widely Used Rice VarietiesYin Wang, Lei Yang, Chang Ma, Yaru Zhou, Mengfei Zhao, Ruiqing Bi, Xiaofei Liang, You-Liang Peng, Jun Yang, Zhensheng Kang, and Guotian LiYin WangState Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China, Lei YangState Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China, Chang MaMinistry of Agriculture and Rural Affairs Key Laboratory of Pest Monitoring and Green Management, Department of Plant Biosecurity, College of Plant Protection, China Agricultural University, Beijing 100193, China, Yaru ZhouState Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China, Mengfei ZhaoState Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China, Ruiqing BiState Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China, Xiaofei Lianghttp://orcid.org/0000-0002-3279-5894State Key Laboratory of Crop Stress Biology for Arid Areas and College of Plant Protection, Northwest A&F University, Yangling 712100, Shaanxi China, You-Liang PengThe MOA Key Laboratory of Molecular Plant Pathology, Department of Plant Pathology, China Agricultural University, Beijing 100094, China, Jun YangMinistry of Agriculture and Rural Affairs Key Laboratory of Pest Monitoring and Green Management, Department of Plant Biosecurity, College of Plant Protection, China Agricultural University, Beijing 100193, China, Zhensheng Kanghttp://orcid.org/0000-0001-5575-0122State Key Laboratory of Crop Stress Biology for Arid Areas and College of Plant Protection, Northwest A&F University, Yangling 712100, Shaanxi China, and Guotian Li†Corresponding author: G. Li; E-mail Address: li4@mail.hzau.edu.cnhttp://orcid.org/0000-0001-6780-7085State Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, ChinaAffiliationsAuthors and Affiliations Yin Wang1 Lei Yang1 Chang Ma2 Yaru Zhou1 Mengfei Zhao1 Ruiqing Bi1 Xiaofei Liang3 You-Liang Peng4 Jun Yang2 Zhensheng Kang3 Guotian Li1 † 1State Key Laboratory of Agricultural Microbiology, Hubei Hongshan Laboratory, the Provincial Key Laboratory of Plant Pathology of Hubei Province, College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China 2Ministry of Agriculture and Rural Affairs Key Laboratory of Pest Monitoring and Green Management, Department of Plant Biosecurity, College of Plant Protection, China Agricultural University, Beijing 100193, China 3State Key Laboratory of Crop Stress Biology for Arid Areas and College of Plant Protection, Northwest A&F University, Yangling 712100, Shaanxi China 4The MOA Key Laboratory of Molecular Plant Pathology, Department of Plant Pathology, China Agricultural University, Beijing 100094, China Published Online:6 Jun 2022https://doi.org/10.1094/MPMI-01-22-0030-AAboutSectionsView articlePDFSupplemental ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat View articleGenome AnnouncementRice blast caused by the filamentous ascomycete fungus Magnaporthe oryzae poses a serious threat to global food security, with worldwide, annual yield losses sufficient to feed more than 60 million people (Cook et al. 2011). In addition to its economic importance, the M. oryzae–rice pathosystem ranks the first in the top 10 of model systems for studies of fungal phytopathogen-host interactions (Dean et al. 2012; Yang et al. 2022). Additionally, M. oryzae is amenable to genetic manipulations, and genomes for both the host and the pathogen are available (Dean et al. 2005; Kawahara et al. 2013). Recently, we sequenced the genome of rice variety Kitaake, which is an ideal candidate as a model system because it has a very short life cycle of 9 weeks (Jain et al. 2019; Li et al. 2017). Lacking was the genome information for M. oryzae strains compatible with Kitaake. Therefore, we sequenced and assembled a high-quality genome for M. oryzae EA18, which is virulent to Kitaake and other commonly used rice varieties.We isolated M. oryzae strains from Kitaake plants that were grown in a local rice blast nursery. Typical blast lesions with abundant conidia on the leaf were sampled, and 55 single-spore strains were isolated, out of which 53 strains developed typical lesions when artificially inoculated on Kitaake in the greenhouse. We then screened for strains that were also virulent to widely used rice varieties Nipponbare and Zhonghua11 (Fig. 1). Finally, strain EA18 with normal conidiation and vegetative growth was selected as the best candidate for genome sequencing (Supplementary Fig. S1; Supplementary Table S1).Fig. 1. The EA18 genome and comparisons with different Magnaporthe oryzae strains. A, Seven EA18 chromosomes on a megabase-pair scale; B, gene density; C, repeat density; and D, GC content. Single-nucleotide polymorphisms and insertion and deletions of EA18 compared with E, 70-15; F, Guy11; and G, P131. The innermost cycle indicates punch inoculation assays on widely used rice varieties Kitaake (Kit), Nipponbare (Nip), and Zhonghua11 (ZH11). Data in circles are displayed in nonoverlapping 100-kb intervals.Download as PowerPointWe sequenced the EA18 genome using Oxford Nanopore technology and the next-generation sequencing (NGS) technique (Wang et al. 2021). In total, 9.7-Gb Nanopore long reads (205× in coverage) and 3.5-Gb NGS short reads (72× in coverage) were used for genome assembly with NextDenovo. The draft genome was polished using both long reads and NGS short reads. The final genome assembly is 45.7 Mb, consisting of 14 contigs (N50 = 5.5 Mb) (Table 1), which covers 96% of the genome estimated using k-mer analysis (Supplementary Fig. S2). Using the RAGOO scaffolding approach (Alonge et al. 2019), these contigs were anchored into seven scaffolds, from which four contigs are gap-free, including chromosomes 2, 4, 5, and 6. Among the seven scaffolds, telomere sequences were identified for both ends of chromosomes 3, 5, and 6 and one end for chromosomes 1 and 4 (Supplementary Table S2). In total, 7.6 Mb of the genome was identified and masked as repetitive sequences using Repeatmasker (v4.0.7) and RepeatModeler (v2.0.1) (Table 1; Supplementary Table S3). BUSCO (v4.0.1) was used to evaluate the completeness of the assembly at the gene space (Seppey et al. 2019), and the results show that the assembly covers 98.3% of the genes in the standard database (Table 1). We mapped the NGS short reads to the genome assembly and filtered the low-quality reads using the parameter '- q 20' (Danecek et al. 2021), and achieved a mapping rate of 100%, indicating the high quality of the EA18 genome.Table 1. Genomic features of EA18VariablesStatisticsSequence coverage205×Genome size (Mb)45.7Contigs14Contig N50 size (bp)5,510,206Maximum contig size (bp)7,860,214Average contig size (bp)3,257,893Number of predicted genes10,988Secreted proteins1,379BUSCO completeness (%)98.3Table 1. Genomic features of EA18View as image HTML To fully annotate genes in EA18, 3.9-Gb RNA-sequencing data were generated from conidia, germ tubes, appressoria, mycelia, and infected tissues at various developmental stages (Zhao et al. 2022). We used an integrative pipeline for gene annotations, in which results from ab initio prediction, homology searches, and transcript-based prediction were integrated using Maker (Holt and Yandell 2011). In total, 10,988 genes were predicted (Table 1), 9,642 of them (87.8%) were annotated using eggNOG (Huerta-Cepas et al. 2019) and 10,436 of them (95.0%) had candidates in M. oryzae 70-15. Genome-wide distribution of genes shows that the majority of genes are located toward the chromosomal ends (Fig. 1B). Using an established pipeline (Dong et al. 2015), we annotated 1,379 secreted proteins in EA18, and 537 of them are candidate effectors, which is similar to the rice blast strain MZ5-1-6 (Luciano et al. 2019).We performed a genome-wide alignment analysis with two high-quality M. oryzae genomes, MZ5-1-6 (Luciano et al. 2019) and 70-15, and the resultant dot plots show that the linear relationships between the genomes are highly established but not for some short blocks, indicating a few rearrangement events present within the M. oryzae genomes (Supplementary Figs. S3 to S5). We compared the genome of EA18 to genomes of three commonly used M. oryzae strains, P131, Guy 11, and 70-15 (Bao et al. 2017; Dean et al. 2005; Xue et al. 2012). Using comparative genomics tools, we identified 33,592, 33,290, and 30,403 small genetic variations, including single-nucleotide polymorphisms (SNPs) and small insertions and deletions (InDels), between EA18 and P131, Guy 11, and 70-15, respectively. Most SNPs and InDels occurred in the upstream and downstream sequences of the gene and the gene spacer (Supplementary Fig. S6). We used Lumpy-SV (Layer et al. 2014) to call structural variants (SVs), which include large insertions, deletions, translocations, inversions, and duplications, between EA18 and other M. oryzae genomes. Comparing with EA18, we called 14,607, 1,989, and 2,148 SVs from P131, Guy 11, and 70-15, respectively, and the P131 strain had more SVs than other strains. In comparison with 70-15, most SVs affect genes and transcript regions (Supplementary Fig. S7). Carbohydrate-active enzymes (CAZymes) are important for fungal infection. They may degrade the plant cell wall and modify the fungal cell wall. We used the dbCAN pipeline (Zhang et al. 2018) to annotate putative CAZymes in the EA18 strain. In total, 771 CAZymes were annotated, including 135 auxiliary activity family proteins, 314 glycoside hydrolases, 103 glycosyltransferases, 6 polysaccharide lyases, 68 carbohydrate esterases, and 57 noncatalytic carbohydrate-binding proteins. In addition, 26 genes may have dual functions and were annotated into two categories of CAZymes.In the current study, we present the high-quality genome of M. oryzae EA18, virulent to multiple widely used rice varieties, including Kitaake, Nipponbare, and Zhonghua11. Particularly with the fast-cycling rice variety Kitaake, the EA18 strain, including its genome, will facilitate studies of the M. oryzae–rice pathosystem and development of innovative disease management strategies for better control of rice blast, a constant threat to rice production worldwide.Data AvailabilityThe EA18 genome file was available at NCBI under the project accession number PRJNA799857, and the annotation file, the variation file, the sequence file of candidate effectors, and the file of comparative analysis to the 70-15 genome have been deposited to the Figshare database. All other sequencing data are available at the National Genomics Data Center Genome Sequence Archive (Chen et al. 2021; CNCB-NGDC Members and Partners 2021) under the accession number PRJCA007759.Author-Recommended Internet ResourcesEnsemblGenomes: http://ensemblgenomes.orgFigshare database: https://doi.org/10.6084/m9.figshare.17704970.v4Kitaake Genome at Phytozome: https://phytozome-next.jgi.doe.govNextDenovo: https://github.com/Nextomics/NextDenovoPathogen Host Interactions database: http://www.phi-base.orgRepeatmasker: http://www.repeatmasker.orgRepeatModeler: https://www.repeatmasker.org/RepeatModeler/Rice Variation Map at RiceVarMap v2.0: http://ricevarmap.ncpgr.cnWhole-genome Sequenced Mutants in Kitaake at KitBase: https://kitbase.ucdavis.eduThe author(s) declare no conflict of interest.Literature CitedAlonge, M., Soyk, S., Ramakrishnan, S., Wang, X., Goodwin, S., Sedlazeck, F. J., Lippman, Z. B., and Schatz, M. C. 2019. RaGOO: Fast and accurate reference-guided scaffolding of draft genomes. Genome Biol. 20:224. https://doi.org/10.1186/s13059-019-1829-6 Crossref, Medline, ISI, Google ScholarBao, J. D., Chen, M. L., Zhong, Z. H., Tang, W., Lin, L. Y., Zhang, X. T., Jiang, H. L., Zhang, D. Y., Miao, C. Y., Tang, H. B., Zhang, J. S., Lu, G. D., Ming, R., Norvienyeku, J., Wang, B. H., and Wang, Z. H. 2017. PacBio sequencing reveals transposable elements as a key contributor to genomic plasticity and virulence variation in Magnaporthe oryzae. Mol. Plant 10:1465-1468. https://doi.org/10.1016/j.molp.2017.08.008 Crossref, Medline, ISI, Google ScholarChen, T., Chen, X., Zhang, S., Zhu, J., Tang, B., Wang, A., Dong, L., Zhang, Z., Yu, C., Sun, Y., Chi, L., Chen, H., Zhai, S., Sun, Y., Lan, L., Zhang, X., Xiao, J., Bao, Y., Wang, Y., Zhang, Z., and Zhao, W. 2021. The genome sequence archive family: Toward explosive data growth and diverse data types. Genomics Proteomics Bioinf. 19:578-583. https://doi.org/10.1016/j.gpb.2021.08.001 Crossref, Medline, ISI, Google ScholarCNCB-NGDC Members and Partners. 2021. Database resources of the National Genomics Data Center, China National Center for Bioinformation in 2021. Nucleic Acids Res. 49:D18-D28. Crossref, Medline, ISI, Google ScholarCook, D. C., Fraser, R. W., Paini, D. R., Warden, A. C., Lonsdale, W. M., and De Barro, P. J. 2011. Biosecurity and yield improvement technologies are strategic complements in the fight against food insecurity. PLoS One 6:e26084. https://doi.org/10.1371/journal.pone.0026084 Crossref, Medline, ISI, Google ScholarDanecek, P., Bonfield, J. K., Liddle, J., Marshall, J., Ohan, V., Pollard, M. O., Whitwham, A., Keane, T., McCarthy, S. A., Davies, R. M., and Li, H. 2021. Twelve years of SAMtools and BCFtools. Gigascience 10:giab008. https://doi.org/10.1093/gigascience/giab008 Crossref, Medline, ISI, Google ScholarDean, R., Van Kan, J. A., Pretorius, Z. A., Hammond-Kosack, K. E., Di Pietro, A., Spanu, P. D., Rudd, J. J., Dickman, M., Kahmann, R., Ellis, J., and Foster, G. D. 2012. The top 10 fungal pathogens in molecular plant pathology. Mol. Plant Pathol. 13:414-430. https://doi.org/10.1111/j.1364-3703.2011.00783.x Crossref, Medline, ISI, Google ScholarDean, R. A., Talbot, N. J., Ebbole, D. J., Farman, M. L., Mitchell, T. K., Orbach, M. J., Thon, M., Kulkarni, R., Xu, J. R., Pan, H., Read, N. D., Lee, Y. H., Carbone, I., Brown, D., Oh, Y. Y., Donofrio, N., Jeong, J. S., Soanes, D. M., Djonovic, S., Kolomiets, E., Rehmeyer, C., Li, W., Harding, M., Kim, S., Lebrun, M. H., Bohnert, H., Coughlan, S., Butler, J., Calvo, S., Ma, L. J., Nicol, R., Purcell, S., Nusbaum, C., Galagan, J. E., and Birren, B. W. 2005. The genome sequence of the rice blast fungus Magnaporthe grisea. Nature 434:980-986. https://doi.org/10.1038/nature03449 Crossref, Medline, ISI, Google ScholarDong, Y., Li, Y., Zhao, M., Jing, M., Liu, X., Liu, M., Guo, X., Zhang, X., Chen, Y., Liu, Y., Liu, Y., Ye, W., Zhang, H., Wang, Y., Zheng, X., Wang, P., and Zhang, Z. 2015. Global genome and transcriptome analyses of Magnaporthe oryzae epidemic isolate 98-06 uncover novel effectors and pathogenicity-related genes, revealing gene gain and lose dynamics in genome evolution. PLoS Pathog. 11:e1004801. https://doi.org/10.1371/journal.ppat.1004801 Crossref, Medline, ISI, Google ScholarHolt, C., and Yandell, M. 2011. MAKER2: An annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinf. 12:491. https://doi.org/10.1186/1471-2105-12-491 Crossref, Medline, ISI, Google ScholarHuerta-Cepas, J., Szklarczyk, D., Heller, D., Hernandez-Plaza, A., Forslund, S. K., Cook, H., Mende, D. R., Letunic, I., Rattei, T., Jensen, L. J., von Mering, C., and Bork, P. 2019. eggNOG 5.0: A hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47:D309-D314. https://doi.org/10.1093/nar/gky1085 Crossref, Medline, ISI, Google ScholarJain, R., Jenkins, J., Shu, S., Chern, M., Martin, J. A., Copetti, D., Duong, P. Q., Pham, N. T., Kudrna, D. A., Talag, J., Schackwitz, W. S., Lipzen, A. M., Dilworth, D., Bauer, D., Grimwood, J., Nelson, C. R., Xing, F., Xie, W., Barry, K. W., Wing, R. A., Schmutz, J., Li, G., and Ronald, P. C. 2019. Genome sequence of the model rice variety KitaakeX. BMC Genomics 20:905. https://doi.org/10.1186/s12864-019-6262-4 Crossref, Medline, ISI, Google ScholarKawahara, Y., de la Bastide, M., Hamilton, J. P., Kanamori, H., McCombie, W. R., Ouyang, S., Schwartz, D. C., Tanaka, T., Wu, J. Z., Zhou, S. G., Childs, K. L., Davidson, R. M., Lin, H. N., Quesada-Ocampo, L., Vaillancourt, B., Sakai, H., Lee, S. S., Kim, J., Numa, H., Itoh, T., Buell, C. R., and Matsumoto, T. 2013. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice (N Y) 6:4. https://doi.org/10.1186/1939-8433-6-4 Google ScholarLayer, R. M., Chiang, C., Quinlan, A. R., and Hall, I. M. 2014. LUMPY: A probabilistic framework for structural variant discovery. Genome Biol. 15:R84. https://doi.org/10.1186/gb-2014-15-6-r84 Crossref, Medline, ISI, Google ScholarLi, G., Jain, R., Chern, M., Pham, N. T., Martin, J. A., Wei, T., Schackwitz, W. S., Lipzen, A. M., Duong, P. Q., Jones, K. C., Jiang, L., Ruan, D., Bauer, D., Peng, Y., Barry, K. W., Schmutz, J., and Ronald, P. C. 2017. The sequences of 1504 mutants in the model rice variety Kitaake facilitate rapid functional genomic studies. Plant Cell 29:1218-1231. https://doi.org/10.1105/tpc.17.00154 Crossref, Medline, ISI, Google ScholarLuciano, L. B. G., Tsai, I. J., Chuma, I., Tosa, Y., Chen, Y. H., Li, J. Y., Li, M. Y., Lu, M. Y. J., Nakayashiki, H., and Li, W. H. 2019. Blast fungal genomes show frequent chromosomal changes, gene gains and losses, and effector gene turnover. Mol. Biol. Evol. 36:1148-1161. https://doi.org/10.1093/molbev/msz045 Crossref, Medline, ISI, Google ScholarSeppey, M., Manni, M., and Zdobnov, E. M. 2019. BUSCO: Assessing genome assembly and annotation completeness. Methods Mol. Biol. 1962:227-245. https://doi.org/10.1007/978-1-4939-9173-0_14 Crossref, Medline, Google ScholarWang, H., Dong, Y., Liao, W., Zhang, X., Wang, Q., Li, G., Xu, J. R., and Liu, H. 2021. High-quality genome resource of Clonostachys rosea Strain CanS41 by oxford nanopore long-read sequencing. Plant Dis. 105:2231-2234. https://doi.org/10.1094/PDIS-12-20-2615-A Link, ISI, Google ScholarXue, M., Yang, J., Li, Z., Hu, S., Yao, N., Dean, R. A., Zhao, W., Shen, M., Zhang, H., Li, C., Liu, L., Cao, L., Xu, X., Xing, Y., Hsiang, T., Zhang, Z., Xu, J. R., and Peng, Y. L. 2012. Comparative analysis of the genomes of two field isolates of the rice blast fungus Magnaporthe oryzae. PLoS Genet. 8:e1002869. https://doi.org/10.1371/journal.pgen.1002869 Crossref, Medline, ISI, Google ScholarYang, L., Zhao, M., Sha, G., Sun, Q., Gong, Q., Yang, Q., Xie, K., Yuan, M., Mortimer, J. C., Xie, W., Wei, T., Kang, Z., and Li, G. 2022. The genome of the rice variety LTH provides insight into its universal susceptibility mechanism to worldwide rice blast fungal strains. Comput. Struct. Biotechnol. J. 20:1012-1026. https://doi.org/10.1016/j.csbj.2022.01.030 Crossref, Medline, ISI, Google ScholarZhang, H., Yohe, T., Huang, L., Entwistle, S., Wu, P., Yang, Z., Busk, P. K., Xu, Y., and Yin, Y. 2018. dbCAN2: A meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 46:W95-W101. https://doi.org/10.1093/nar/gky418 Crossref, Medline, ISI, Google ScholarZhao, J., Sun, P., Sun, Q., Li, R., Qin, Z., Sha, G., et al. 2022. The MoPah1 phosphatidate phosphatase is involved in lipid metabolism, development, and pathogenesis in Magnaporthe oryzae. Mol. Plant Pathol. 23:720-732. https://doi.org/10.1111/mpp.13193 Crossref, Medline, ISI, Google ScholarY. Wang and L. Yang contributed equally to this work.Funding: This work was supported by Fundamental Research Funds for the Central Universities (2662020ZKPY006) and the National Natural Science Foundation of China (32172373 and 31801723) to G. Li, and by Hubei Hongshan Laboratory.The author(s) declare no conflict of interest. Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.DetailsFiguresLiterature CitedRelated Vol. 35, No. 8 August 2022ISSN:0894-0282e-ISSN:1943-7706 Download Metrics Article History Issue Date: 31 Aug 2022Published: 6 Jun 2022Accepted: 3 Apr 2022 Pages: 727-730 InformationCopyright © 2022 The Author(s).This is an open access article distributed under the CC BY-NC-ND 4.0 International license.Funding Fundamental Research Funds for the Central UniversitiesGrant/Award Number: 2662020ZKPY006 National Natural Science Foundation of ChinaGrant/Award Number: 32172373Grant/Award Number: 31801723 Hubei Hongshan Laboratory Keywordsgenome assemblyKitaakeNanoporeOryza sativaPyricularia oryzaerice blastThe author(s) declare no conflict of interest.PDF downloadCited byA High-Quality Genome of Rhizoctonia solani, a Devastating Fungal Pathogen with a Wide Host RangeQun Yang, Lei Yang, Yin Wang, Yilyu Chen, Keming Hu, Wei Yang, Shimin Zuo, Jiandi Xu, Zhensheng Kang, Xueqiong Xiao, and Guotian Li29 September 2022 | Molecular Plant-Microbe Interactions, Vol. 35, No. 10Editorial: Molecular interactions between crops and phytopathogens, volume II: Rice29 July 2022 | Frontiers in Plant Science, Vol. 13

Abstract:As one of the most important staple food crops worldwide, rice is among the first plant species whose genomes were sequenced. The reference genome for Oryza sativa ssp. japonica variety Nipponbare was completed by the International Rice Genome Sequencing Project (IRGSP) using a bacterial artificial chromosome (BAC)-based cloning strategy (Matsumoto et al., 2005Matsumoto T. Wu J.Z. Kanamori H. Katayose Y. Fujisawa M. Namiki N. Mizuno H. Yamamoto K. Antonio B.A. Baba T. et al.The map-based sequence of the rice genome.Nature. 2005; 436: 793-800Crossref PubMed Scopus (2968) Google Scholar). IRGSP1.0 was the latest release (Kawahara et al., 2013Kawahara Y. de la Bastide M. Hamilton J.P. Kanamori H. McCombie W.R. Ouyang S. Schwartz D.C. Tanaka T. Wu J.Z. Zhou S.G. et al.Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data.Rice. 2013; 6: 4Crossref Scopus (1023) Google Scholar) as the result of iterative updates. Rice Oryza sativa ssp. indica is estimated to account for over 70% of the world's rice production. A representative variety of indica, 93-11, was sequenced by Beijing Institute of Genome using a whole-genome shotgun approach and also underwent updates (Yu et al., 2002Yu J. Hu S. Wang J. Wong G.K. Li S. Liu B. Deng Y. Dai L. Zhou Y. Zhang X. et al.A draft sequence of the rice genome (Oryza sativa L. ssp. indica).Science. 2002; 296: 79-92Crossref PubMed Scopus (2557) Google Scholar, Gao et al., 2013Gao Z.Y. Zhao S.C. He W.M. Guo L.B. Peng Y.L. Wang J.J. Guo X.S. Zhang X.M. Rao Y.C. Zhang C. et al.Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences.Proc. Natl. Acad. Sci. USA. 2013; 110: 14492-14497Crossref PubMed Scopus (130) Google Scholar). Assembly errors and DNA structural defects are major issues dogging sequenced genomes today, more so for indica 93-11 (Yu et al., 2006Yu J. Ni P.X. Wong G.K.S. Comparing the whole-genome-shotgun and map-based sequences of the rice genome.Trends Plant Sci. 2006; 11: 387-391Abstract Full Text Full Text PDF PubMed Scopus (14) Google Scholar), as issues at sequence structural level, i.e., insertions, deletions, translocations, and copy number variations are difficult to resolve with short-read sequencing technology. Unfilled gaps that contain complex sequences (often referred to as dark matter) remained with both rice genomes. Array-based comparative genomic hybridization (aCGH) (Bruce et al., 2009Bruce M. Hess A. Bai J.F. Mauleon R. Diaz M.G. Sugiyama N. Bordeos A. Wang G.L. Leung H. Leach J.E. Detection of genomic deletions in rice using oligonucleotide microarrays.BMC Genomics. 2009; 10: 129Crossref PubMed Scopus (37) Google Scholar) and physical mapping with BAC libraries (Oryza Map Alignment Project) (Jacquemin et al., 2013Jacquemin J. Bhatia D. Singh K. Wing R.A. The International Oryza Map Alignment Project: development of a genus-wide comparative genomics platform to help solve the 9 billion-people question.Curr. Opin. Plant Biol. 2013; 16: 147-156Crossref PubMed Scopus (92) Google Scholar) were among the first technologies used to detect and catalog sequence structural features for O. sativa varieties and relatives in genus Oryza. Another approach, optical mapping, represented by the latest BioNano technology (Lam et al., 2012Lam E.T. Hastie A. Lin C. Ehrlich D. Das S.K. Austin M.D. Deshpande P. Cao H. Nagarajan N. Xiao M. et al.Genome mapping on nanochannel arrays for structural variation analysis and sequence assembly.Nat. Biotechnol. 2012; 30: 771-776Crossref PubMed Scopus (436) Google Scholar), offered significant advantages by revealing enzyme recognition site patterns across long DNA molecules. BioNano enabled high resolution and high throughput using nanochannel chips and fluorescent labeling, and rapidly expanded its applications to the study of large complex genomes (Cao et al., 2014Cao H.Z. Hastie A.R. Cao D.D. Lam E.T. Sun Y.H. Huang H.D. Liu X. Lin L.Y. Andrews W. Chan S. et al.Rapid detection of structural variation in a human genome using nanochannel-based genome mapping technology.Gigascience. 2014; 3: 34Crossref PubMed Scopus (105) Google Scholar, Stankova et al., 2016Stankova H. Hastie A.R. Chan S. Vrana J. Tulpova Z. Kubalakova M. Visendi P. Hayashi S. Luo M.C. Batley J. et al.BioNano genome mapping of individual chromosomes supports physical mapping and sequence assembly in complex plant genomes.Plant Biotechnol. J. 2016; 14: 1523-1531Crossref PubMed Scopus (68) Google Scholar). We also found that the BAC-based approach was sensitive to detect small-size structural changes, which was complementary to BioNano (Supplemental Notes). Using BioNano technology, we first built high-resolution optical maps for the genomes of Nipponbare and 93-11. We generated 71.4 and 64.7 Gb BioNano single-molecule data for Nipponbare and 93-11, respectively (Supplemental Table 1). To ensure the quality of BioNano data for subsequent analyses, the single-molecule data were filtered with a quality filter before they were used to assemble genome maps (termed BioNano maps) with stringent parameters (Supplemental Methods). For Nipponbare, a consensus genome map of 377.4 Mb was generated with 482 contigs and an N50 of 1.1 Mb (Supplemental Table 1 and Supplemental Dataset 1). It matched IRGSP1.0 at ∼96.0% (Supplemental Table 2). Similarly, a consensus genome map for 93-11 of 393.5 Mb was assembled with 394 contigs and an N50 of 1.4 Mb (Supplemental Table 1 and Supplemental Dataset 2). However, it had a much lower mapping rate (85.2%) (Supplemental Table 2) than the 93-11 reference, 93-11v2 (Gao et al., 2013Gao Z.Y. Zhao S.C. He W.M. Guo L.B. Peng Y.L. Wang J.J. Guo X.S. Zhang X.M. Rao Y.C. Zhang C. et al.Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences.Proc. Natl. Acad. Sci. USA. 2013; 110: 14492-14497Crossref PubMed Scopus (130) Google Scholar). Given the differential mapping rates, we further compared the strings of labeling motifs in BioNano maps against corresponding genome regions and investigated discrepancies from diverging signals that represent likely assembly defects (Supplemental Methods). In IRGSP1.0, a catalog of 262 discrepancies with a total size of ∼8.78 Mb were found, which were categorized into four groups: insertions (I), deletions (II), inversions or translocations (III), and other complex discordances (IV) marked by restriction patterns distinctly different from BioNano maps (Figure 1A left and Supplemental Dataset 3). Note that chromosomes 1 and 4 had the most discrepancies (38 and 37, respectively), whereas chromosomes 3 and 6 had the fewest. To trace how Nipponbare reference evolved, the old version of Nipponbare reference, IRGSPv5 (Matsumoto et al., 2005Matsumoto T. Wu J.Z. Kanamori H. Katayose Y. Fujisawa M. Namiki N. Mizuno H. Yamamoto K. Antonio B.A. Baba T. et al.The map-based sequence of the rice genome.Nature. 2005; 436: 793-800Crossref PubMed Scopus (2968) Google Scholar), was also analyzed. A total of 269 discrepancies with a total size of ∼11.04 Mb were identified (Supplemental Dataset 3). Among them, five were corrected in IRGSP1.0. For example, BioNano contig29 aligned perfectly with IRGSP1.0 in region chr10:7961923-8042411, but diverged from IRGSPv5 in region chr10:7894177-7974664 (Figure 1B, upper left). Unexpectedly, we found newly introduced assembly errors in IRGSP1.0 that did not exist in IRGSPv5 (Figure 1B, lower left). Compared with IRGSP1.0, 93-11v2 had a greater number of discrepancies with its BioNano map; 3419 were identified (∼105.17 Mb total region size) and categorized into the same four groups as defined for IRGSP1.0 (Figure 1A, right, and Supplemental Dataset 4). Chromosomes 1 and 3 had the most discrepancies (428 and 357, respectively), whereas chromosomes 9 and 12 had the fewest. Further analysis of the old version of the 93-11 reference, 93-11v1 (Yu et al., 2006Yu J. Ni P.X. Wong G.K.S. Comparing the whole-genome-shotgun and map-based sequences of the rice genome.Trends Plant Sci. 2006; 11: 387-391Abstract Full Text Full Text PDF PubMed Scopus (14) Google Scholar) with the 93-11 BioNano map, identified 2566 discrepancies with a total size of ∼106.77 Mb (Supplemental Dataset 4). Among them, 333 (∼13.51 Mb) were corrected in 93-11v2. For example, the BioNano contig4 aligned well with 93-11v2 in region chr02:31626328-31669232, but diverged from 93-11v1 in region chr02:33054748-33127908 because of an erroneous insertion (Figure 1B, upper right). Surprisingly, 463 discrepancies (∼13.16 Mb) were newly introduced into 93-11v2 that did not exist in 93-11v1. One example pointed to a translocation in 93-11v2 (Figure 1B, lower right). With the newly available Zhenshan 97 and Minghui 63 genomes (Zhang et al., 2016Zhang J. Chen L.-L. Xing F. Kudrna D.A. Yao W. Copetti D. Mu T. Li W. Song J.-M. Xie W. Extensive sequence divergence between the reference genomes of two elite indica rice varieties Zhenshan 97 and Minghui 63.Proc. Natl. Acad. Sci. USA. 2016; 113: E5163-E5171Crossref PubMed Scopus (154) Google Scholar), we found 1728 and 1710 of the 93-11v2 discrepancy regions had good matches between the 93-11 BioNano map and the genome of either Zhenshan 97 or Minghui 63, respectively (Supplemental Notes), pointing to common genomic elements among indica rice and affirming the accuracy of the BioNano map. In retrospect of both Nipponbare and 93-11 genomes, our results indicated that although their quality at nucleotide sequence level greatly improved with the advance of sequencing technology, progress was limited on issues associated with structural complexity in rice genomes. We next focused on three specific types of sequence structural issues in rice genomes, namely complex sequence structure in gaps, highly repetitive sequences, and BAC assembly errors. The rice reference genomes, particularly 93-11v2, have large gaps that often contain complex sequences. There were 39 and 9399 large gaps in Nipponbare and 93-11 genomes that were analyzed with BioNano maps (Supplemental Dataset 5). In IRGSP1.0, a gap in chr05:27174270-27175270 was found to have a 40-kb-long inverted repeat, unraveled by Nipponbare BioNano map (Figure 1C top). In 93-11v2, an example was found with chr09, where a 544-kb large insertion contained several 35-kb or 47-kb repeats, arranged alternately (Figure 1C, bottom). Note these complex sequences in gaps were supported by high coverage of single-molecule data. Highly repetitive sequences are almost intractable to short-read sequencing technologies. As a result, assemblies are heavily biased against repeats and duplications. To map out the tandem repeats in rice genomes, we systematically scanned the BioNano genome maps (Supplemental Methods) and found that both rice genomes, IRGSP1.0 and 93-11v2, underestimated frequencies of tandem repeats (Figure 1D and Supplemental Dataset 6). The peaks of 8.2-kb observed in Nipponbare and of 9.1-kb in 93-11 were likely repeats of different types of rDNA genes, concentrating on chromosome 9 in Nipponbare, and chromosomes 9 and 10 in 93-11 (Supplemental Figure 1). Based on BioNano data, we estimated the frequency of total repeats for rDNA genes in Nipponbare was ∼141, and in 93-11 was ∼215. Notably BioNano maps also helped refine different types of repeats/duplications in rice genomes (Supplemental Figure 2). Historically, genome assemblies relied on the BAC scaffolding approach. And BAC scaffolding errors could lead to sequence structural defects in rice genomes. One case was found with IRGSP1.0 in region Chr11:28547360-28644282, where BAC_3335 (AC120507) was missing (Figure 1E). Based on BioNano map contig640, BAC_3335 that was constructed from 23 contigs to form four repeats units locates between BAC_3334 and BAC_3336. Lastly, we developed a bioinformatics pipeline, BASCGF (BioNano Assisted Sequence Correction and Gap Filling), to reconstruct the 93-11 reference genome, taking advantage of the robust BioNano map resource and the homologous sequences shared between 93-11 and Nipponbare genomes. We reasoned it was possible to correct the defects in the 93-11 genome by referencing the high-quality Nipponbare sequences via the BioNano genome maps. We analyzed the 3419 discrepancies (Supplemental Dataset 4) identified in 93-11v2 by comparing their BioNano map fragments with the corresponding regions in the Nipponbare genome. Remarkably, 873 of 3419 were found to match between the 93-11 BioNano map and the Nipponbare genome, IRGSP1.0. An example is shown with the 93-11 BioNano map contig788, which aligned perfectly with IRGSP1.0 in chr12:1174639-1197895, but did not align with 93-11v2 (Figure 1F). Based on these results, BASCGF was designed to correct errors and fill in gaps in the 93-11 genome with the homologous Nipponbare sequences that were readily validated and polished with Illumina sequencing data generated from 93-11 DNA (Supplemental Methods and Supplemental Figure 3). The reconstructed 93-11 genome, 93-11v2Plus, had a total size of 372.48 Mb (Supplemental Dataset 8). It had 12.93-Mb sequences corrected, 5703 gaps closed, and 455 gaps partially filled (Supplemental Table 3), representing a significant improvement over 93-11v2. Changes were made to 864 genes, while 28 genes were newly added in filled gaps, and 836 were removed from existing loci in 93-11v2 and placed in the undetermined category in 93-11v2Plus. The above example (Figure 1F) was fixed in 93-11v2Plus by closing six gaps and removing the 23.3-kb erroneously anchored sequences. In summary, using BioNano technology, we built high-resolution genome map resources for rice Nipponbare and 93-11. We demonstrated the utility of BioNano maps in elucidating complex sequence structure in gaps, resolving highly repetitive sequences, and correcting BAC assembly errors. An improved 93-11 reference, 93-11v2Plus, was reconstructed with a bioinformatics pipeline, BASCGF. These results validate the value of the BioNano map resource and our strategy for BioNano map-based analysis. Considering the strengths of different sequencing technologies, we would suggest that the best-quality genome can be achieved by integrating BioNano mapping and PacBio and Illumina sequencing, i.e., building longer contigs with PacBio long reads, better scaffolding/assembly with BioNano map, and clearing sequence ambiguity/variation at single-nucleotide resolution with Illumina short reads (Supplemental Notes). Our study also pinpoints that dissecting genome structural complexity should be considered as a new priority in future genome studies. This work is supported in part by grants from the National Key Basic Research Program in China (2013CB127005), the National Natural Science Foundation of China (nos. 31401128, 31571310), and the Ministry of Agriculture of China (2016ZX08010-002), and by Special Fund for Strategic Pilot Technology Chinese Academy of Sciences (XDA08020104).

PacBio Sequencing Reveals Transposable Elements As a Key Contributor to Genomic Plasticity and Virulence Variation in Magnaporthe Oryzae

Transposable elements impact the population divergence of rice blast fungus Magnaporthe oryzae

Single-Molecule Sequencing Assists Genome Assembly Improvement and Structural Variation Inference

Genome Sequence of Magnaporthe Oryzae EA18 Virulent to Multiple Widely Used Rice Varieties.

Comparative Analysis of the Genomes of Two Field Isolates of the Rice Blast Fungus Magnaporthe Oryzae

Transposable elements impact the population divergence of rice blast fungus

Comparative Analysis of the Genome of the Field Isolate V86010 of the Rice Blast Fungus Magnaporthe Oryzae from Philippines

Dynamic Nucleotide-Binding Site and Leucine-Rich Repeat-Encoding Genes in the Grass Family

Genetic Variation Bias Toward Noncoding Regions and Secreted Proteins in the Rice Blast Fungus Magnaporthe Oryzae

Insertion of Transposable Elements in AVR-Pib of Magnaporthe oryzae Leading to LOSS of the Avirulent Function

Transposons are important contributors to gene expression variability under selection in rice populations

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

Global genome and transcriptome analyses of Magnaporthe oryzae epidemic isolate 98-06 uncover novel effectors and pathogenicity-related genes, revealing gene gain and lose dynamics in genome evolution.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Natural genomic variation in rice blast genomes is associated with specific heterochromatin modifications

BioNano Genome Map Resource for Oryza Sativa Ssp. Japonica and Indica and Its Application in Rice Genome Sequence Correction and Gap Filling

Variable genome evolution in fungi after transposon-mediated amplification of a housekeeping gene

Whole genome regulatory effect of MoISW2 and consequences for the evolution of the rice plant pathogenic fungus Magnaporthe oryzae

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria

Precise Editing Enables Crop Broad-Spectrum Resistance

Distinct genomic contexts predict gene presence-absence variation in different pathotypes of Magnaporthe oryzae