Improved Whole Genome Sequence Resource of Phytophthora Melonis CJ26, Using Long-Read Nanopore Sequencing.
Shengdan Lu,Jiandong Bao,Zhixin Wang,Zhiting Li,Wanzhen Feng,Chengdong Yang,Qinghe Chen
DOI: https://doi.org/10.1094/mpmi-04-22-0096-a
IF: 3.4218
2022-01-01
Molecular Plant-Microbe Interactions
Abstract:HomeMolecular Plant-Microbe Interactions®Vol. 35, No. 10Improved Whole Genome Sequence Resource of Phytophthora melonis CJ26, Using Long-Read Nanopore Sequencing Previous RESOURCE ANNOUNCEMENT OPENOpen Access licenseImproved Whole Genome Sequence Resource of Phytophthora melonis CJ26, Using Long-Read Nanopore SequencingShengdan Lu, Jiandong Bao, Zhixin Wang, Zhiting Li, Wanzhen Feng, Chengdong Yang, and Qinghe ChenShengdan LuSanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, ChinaSearch for more papers by this author, Jiandong Baohttp://orcid.org/0000-0003-3423-5118The Institute of Plant Protection and Microbiology, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, ChinaSearch for more papers by this author, Zhixin WangSanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, ChinaSearch for more papers by this author, Zhiting LiSanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, ChinaSearch for more papers by this author, Wanzhen FengSanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, ChinaKey Laboratory of Green Prevention and Control of Tropical Plant Diseases and Pests, Ministry of Education, Hainan University, Haikou 570228, ChinaSearch for more papers by this author, Chengdong YangSanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, ChinaKey Laboratory of Green Prevention and Control of Tropical Plant Diseases and Pests, Ministry of Education, Hainan University, Haikou 570228, ChinaSearch for more papers by this author, and Qinghe Chen†Corresponding author: Q. Chen; E-mail Address: qhchen@hainanu.edu.cnhttp://orcid.org/0000-0002-1797-7614Sanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, ChinaKey Laboratory of Green Prevention and Control of Tropical Plant Diseases and Pests, Ministry of Education, Hainan University, Haikou 570228, ChinaSearch for more papers by this authorAffiliationsAuthors and Affiliations Shengdan Lu1 Jiandong Bao2 Zhixin Wang1 Zhiting Li1 Wanzhen Feng1 3 Chengdong Yang1 3 Qinghe Chen1 3 † 1Sanya Nanfan Research Institute, College of Plant Protection, Hainan University, Hainan Yazhou Bay Seed Laboratory, Sanya 572025, China 2The Institute of Plant Protection and Microbiology, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, China 3Key Laboratory of Green Prevention and Control of Tropical Plant Diseases and Pests, Ministry of Education, Hainan University, Haikou 570228, China Published Online:27 Sep 2022https://doi.org/10.1094/MPMI-04-22-0096-AAboutSectionsPDF ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat Genome AnnouncementThe oomycete Phytophthora melonis, which was first described by Katsura (1976), is one of the most destructive pathogens of cucumber and other cucurbits, causing damping-off and significantly crop yield reductions worldwide (Erwin and Ribeiro,1996). Symptoms of P. melonis infection include root collar rot, stem lesions, foliar blight and fruit rot, which finally resulting in plant death (Hashemi et al. 2019). The pathogen causes blight, dieback, and crown rot disease (Mirabolfathy et al. 2001) on several economically important crops, such as watermelon (Citrullus lanatus), muskmelon (Cucumis melo L.), wax gourd (Benincasa hispida), zucchini (Cucurbita pepo L.), gourd (Trichosanthes dioica Roxb.), and pistachio (Pistacia vera L.). Despite intensive research efforts made in the past years, there is still no effective cure currently available for this disease and little is known about the molecular and genetic mechanisms governing disease development.Until recently, only one genome assembly of P. melonis CPHST-BL 23 (BL23) (host: Cucumis sativus) has been reported (Srivastava et al. 2020). The genome of strain BL23 was sequenced, using the Oxford Nanopore Technologies (ONT) MinION platform, and was assembled into a 112-Mb genome assembly consisting of 2,091 contigs with N50 of 107 kb (Table 1). To better understand the pathogenicity of P. melonis, more high-quality genome resources of P. melonis were urgently needed.Table 1. Genomic features of Phytophthora melonisFeaturesCJ26CPHST BL 23aOxford Nanopore Technologies (ONT) long reads (Gb)6.815.13N50 of ONT reads13,220–Coverage of ONT reads65×50×Illumina short reads (bp)3,945,876,600 (about 38×)–Mapping rate of Illumina short reads96.16%95.29%Assembly size (bp)105,001,437112,010,074Contig number942,091Contig N50 (bp)2,135,025107,148Contig L5016276Maximum contig length (bp)5,998,312720,402GC content (%)54.9254.73Repeat sequence (%)49.0246.46Protein-coding genes29,80134,563Pfamb14,593–GOb7,571–KEGGb5,635–KOGb17,028–Carbohydrate-active enzymesc444–Cytochrome P450 enzymes27–Candidate RXLR effectorsd687–Pathogen-host interaction genes2,010–aPrevious reported genome assembly (National Center for Biotechnology Information accession GCA_008553785.1). A ‘–’ indicates data unavailable.bPfam and GO were annotated with InterProScan V5.52-86.0, and KEGG and KOG were conducted with EggNogMapper V2.cCarbohydrate-active enzymes were annotated by dbCAN2.dRXLR effectors were identified using the open reading frame–based method following our previous study (Wang et al. 2021).Table 1. Genomic features of Phytophthora melonisView as image HTML Here, we isolated a strain, P. melonis CJ26, from infected stems of muskmelon at Changjiang County, Hainan Province, China. This strain was confirmed as P. melonis by combining morphological characteristics and multi-loci phylogenetic analyses based on the sequences of internal transcribed spacer (EU088256.2, 100%), β-tubulin (EF372616.1, 99.79%), and transcription elongation factor (EF372614.1, 99.78%) genes. Colony morphology showed slightly radiate to stellate patterns, with cottony aerial mycelium, on V8 agar plates (200 ml of clarified V8 juice, 800 ml of distilled water, and 20 g of agar). For whole-genome sequencing, high-quality genomic DNA was extracted from the mycelia grown on V8 media. We obtained 6.81 Gb (about 65×; N50 13,220 bp) of long reads, using the ONT PromethION sequencing platform, and 3.95 Gb (about 38×) of 150-bp paired-end short reads (fragment size about 350 bp), using the Illumina HiSeq 3000 sequencing platform (Table 1).The de novo draft genome assembly of strain CJ26 was generated by NextDenovo V2.5.0, using ONT long reads, and was further polished by NextPolish v1.3.1 (Hu et al. 2020), using both ONT long reads and Illumina short reads. Finally, we obtained a 105-Mb polished genome assembly (GC content 54.92%) consisting of 94 contigs (Table 1). The N50 length of contigs was 2.14 Mb (L50 16), while the maximum contig length of this assembly was 6.00 Mb (Table 1). We found that CJ26 has a similar assembly size as BL23 (105 vs. 112 Mb), while with an approximately 95% reduced number of contigs (94 vs. 2,091), approximately 20-fold increased contig N50 (2.14 Mb vs. 107 kb), and approximately eightfold increased maximum contig length (6.00 Mb vs. 720 kb) (Table 1), suggesting significant improvement of the assembly quality.Genome completeness evaluation with BUSCO v5.3.0 (Manni et al. 2021) showed that the genome assembly of strain CJ26 contained 153 (89.47%), 229 (89.80%), and 93 (93.00%) complete single-copy orthologs at alveolata (n = 171), eukaryota (n = 255), and stramenopiles (n = 100) level, respectively (Fig. 1A). The BUSCO genome completeness of strain CJ26 is slightly higher than that of the strain BL23 (Fig. 1A). Genome completeness was also tested by the mapping rate of Illumina short reads, using BWA v0.7.17 (Li et al. 2009). A total of 96.16% of strain CJ26 Illumina short reads were mapped to the CJ26 genome assembly (Table 1). To evaluate genome similarities between CJ26 and BL23, we mapped CJ26 Illumina short reads to BL23 genome assembly, and the mapping rate was up to 95.29%. These results indicate high completeness of the CJ26 genome assembly and high genome similarity with strain BL23.Fig. 1. Comparative genomic analyses between strains CJ26 and CPHST-BL 23 (BL23). A, Genome completeness evaluated by BUSCO v5.3.0 with alveolata, eukaryota, and stramenopiles lineages. B, Whole-genome orthologous gene clusters analysis conducted by OrthoVenn2. The red number in brackets indicates the single-copy orthologous gene clusters. C, Whole-genome collinearity between CJ26 and BL23 using MCScanX and TBtools (shows only the large contigs >1 Mb in CJ26).Download as PowerPointA total of 51.47 Mb (49.02%) repeat sequences was identified in the genome assembly of CJ26 using RepeatMasker v4.1.2 with a de novo repeat library generated by RepeatModeler v2.02 (Table 1). The repeats were mainly enriched with interspersed repeats (49.87 Mb, 96.89%), including 31.97 Mb of long terminal repeats (62.11%), 7.21 Mb of DNA elements (14.01%), and 9.67 Mb of unclassified interspersed repeats (18.79%). The repeat-masked genome assembly was used for gene annotation, and a total of 29,801 protein-coding genes were identified by BRAKER v2.1.6 (Bruna et al. 2021), which integrated gene structure evidence from Augustus v3.4.0 (Stanke et al. 2008) and GeneMark-ET (Bruna et al. 2020) with the available Phytophthora homologous proteins as a training dataset (Table 1). The gene number of strain CJ26 is approximately 16% (29,801 vs. 34,563) lower than strain BL23 annotated using the same gene prediction pipeline.To figure out the genome difference between strains from different host plants, whole-genome orthologs analysis was conducted with OrthoVenn2 (Xu et al. 2019), and the result showed that strains CJ26 and BL23 shared 18,712 (approximately 94%) orthologous gene clusters, among which, 15,810 clusters (approximately 80%) were single-copy orthologs (Fig. 1B). Gene-based genome collinearity analysis using MCScanX (Wang et al. 2012) and TBtools (Chen et al. 2020) showed that strains CJ26 and BL23 not only have high genetic similarity but also high gene synteny (Fig. 1C).Gene functional annotation was conducted with eggNOG-mapper v2.1.6 (Letunic et al. 2021) and InterProScan v5.52-86.0 (Jones et al. 2014) with a set of popular bioinformatic databases including Pfam v33.1, Gene Ontology (GO) release 2021, Kyoto Encyclopedia of Genes and Genomes (KEGG) release 100.0, and Eukaryotic Orthologous Groups (KOG), which annotated functions of 14,593, 7,571, 5,635, and 17,028 genes, respectively (Table 1). We also identified a set of pathogenicity-related genes, including 2,010 genes from Pathogen-Host Interaction base (PHI-base) V4.12, 444 carbohydrate-active enzymes detected by dbCAN2 (Lombard et al. 2014), 687 putative RXLR effectors, following our previous open reading frame (ORF)-based pipeline (Wang et al. 2021), and 27 cytochrome P450 enzymes from Pfam annotations (Table 1).Data AvailabilityThe whole-genome assembly and genes reported in this paper have been deposited in the Genome Warehouse (accession number GWHBGBS00000000) at the National Genomics Data Center, China National Center for Bioinformation (CNCB-NGDC Members and Partners 2021). The genome assembly of CJ26 has also been submitted to National Center for Biotechnology Information assembly database (accession number JAMXKT000000000 under BioProject PRJNA847706). The raw sequence reads also are publicly accessible at the CNCB-NGDC Genome Sequence Archive (https://ngdc.cncb.ac.cn/gsa/) under accession number CRA005555 (BioProject PRJCA007515).Author-Recommended Internet ResourcesBUSCO v5.3.0 data: https://busco-data.ezlab.org/v5/dataGenome Sequence Archive: https://ngdc.cncb.ac.cn/gsaGenome Warehouse: https://ngdc.cncb.ac.cn/gwhGO database: http://geneontology.orgNextDenovo V2.5.0: https://github.com/NextomicsKEGG database: https://www.genome.jp/keggKOG database: https://ftp.ncbi.nih.gov/pub/COG/KOGPfam database: https://pfam.xfam.orgPHI-base: http://www.phi-base.orgRepeatMasker v4.1.2: http://www.repeatmasker.org/RepeatModeler v2.02: http://www.repeatmasker.org/RepeatModelerThe author(s) declare no conflict of interest.Literature CitedBruna, T., Hoff, K. J., Lomsadze, A., Stanke, M., and Borodovsky, M. 2021. BRAKER2: Automatic eukaryotic genome annotation with GeneMark-EP+ and AUGUSTUS supported by a protein database. NAR Genomics Bioinforma. 3:lqaa108. https://doi.org/10.1093/nargab/lqaa108 Crossref, Medline, Google ScholarBruna, T., Lomsadze, A., and Borodovsky, M. 2020. GeneMark-EP+: Eukaryotic gene prediction with self-training in the space of genes and proteins. NAR Genomics Bioinforma. 2:lqaa026. https://doi.org/10.1093/nargab/lqaa026 Crossref, Medline, Google ScholarChen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., and Xia, R. 2020. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 13:1194-1202. https://doi.org/10.1016/j.molp.2020.06.009 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. https://doi.org/10.1093/nar/gkaa1022 Crossref, Medline, ISI, Google ScholarErwin, D. C., and Ribeiro, O. K. 1996. Phytophthora Diseases Worldwide. APS Press, St. Paul, MN, U.S.A. Google ScholarHashemi, L., Golparvar, A. Reza., Esfahani, M. N., and Golabadi, M. 2019. Correlation between cucumber genotype and resistance to damping-off disease caused by Phytophthora melonis. Biotechnol. Biotec. Eq. 33:1494-1504. https://doi.org/10.1080/13102818.2019.1675535 Crossref, ISI, Google ScholarHu, J., Fan, J., Sun, Z., and Liu, S. 2020. NextPolish: A fast and efficient genome polishing tool for long-read assembly. Bioinformatics 36:2253-2255. https://doi.org/10.1093/bioinformatics/btz891 Crossref, Medline, ISI, Google ScholarJones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., McWilliam, H., Maslen, J., Mitchell, A., Nuka, G., Pesseat, S., Quinn, A. F., Sangrador-Vegas, A., Scheremetjew, M., Yong, S. Y., Lopez, R., and Hunter, S. 2014. InterProScan 5: Genome-scale protein function classification. Bioinformatics 30:1236-1240. https://doi.org/10.1093/bioinformatics/btu031 Crossref, Medline, ISI, Google ScholarKatsura, K. 1976. Two new species of Phytophthora causing damping-off of cucumber and trunk rot of chestnut. Trans. Mycol. Soc. Jap. 17:238-242. Google ScholarLetunic, I., Bork, P., and Huerta-Cepas, J. 2021. eggNOG-mapper v2: Functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38:5825-5829. https://doi.org/10.1093/molbev/msab293 Crossref, Medline, ISI, Google ScholarLi, H., and Durbin, R. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25:1754-1760. https://doi.org/10.1093/bioinformatics/btp324 Crossref, Medline, ISI, Google ScholarLombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., and Henrissat, B. 2014. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42:D490-D495. https://doi.org/10.1093/nar/gkt1178 Crossref, Medline, ISI, Google ScholarManni, M., Berkeley, M. R., Seppey, M., Simão, F. A., and Zdobnov, E. M. 2021. BUSCO update: Novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol. Biol. Evol. 38:4647-4654. https://doi.org/10.1093/molbev/msab199 Crossref, Medline, ISI, Google ScholarMirabolfathy, M., Cooke, D. E. L., Duncan, J. M., Williams, N. A., Ershad, D., and Alizadeh, A. 2001. Phytophthora pistaciae sp. nov. and P. melonis: The principal causes of pistachio gummosis in Iran. Mycol. Res.105:1166-1175. https://doi.org/10.1016/S0953-7562(08)61987-5 Crossref, ISI, Google ScholarSrivastava, S. K., Abad, Z. G., Knight, L. M., Zeller, K., Mavrodieva, V., and Nakhla, M. 2020. Draft genome resource for the ex-types of Phytophthora ramorum, P. kernoviae, and P. melonis, species of regulatory concern, using ultra-long read MinION nanopore sequencing. Mol. Plant-Microbe In. 33:794-797. https://doi.org/10.1094/MPMI-12-19-0342-A Link, ISI, Google ScholarStanke, M., Diekhans, M., Baertsch, R., and Haussler D. 2008. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24:637-644. https://doi.org/10.1093/bioinformatics/btn013 Crossref, Medline, ISI, Google ScholarWang, Y., Tang, H., Debarry, J. D., Tan, X., Li, J., Wang, X., Lee, T. H., Jin, H., Marler, B., Guo, H., Kissinger, J. C., and Paterson, A. H. 2012. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40:e49. https://doi.org/10.1093/nar/gkr1293 Crossref, Medline, ISI, Google ScholarWang, Z., Bao, J., Lv, L., Lin, L., Li, Z., Shi, M., Huang, Y., Wang, R., Li, B., Liu, P., and Chen, Q. 2021. Genome sequence resource of Phytophthora colocasiae from China using Nanopore sequencing technology. Plant Dis. 105:4141-4145. https://doi.org/10.1094/PDIS-11-20-2327-A Link, ISI, Google ScholarXu, L., Dong, Z., Fang, L., Luo, Y., Wei, Z., Guo, H., Zhang, G., Gu, Y. Q., Coleman-Derr, D., Xia, Q., and Wang, Y. 2019. OrthoVenn2: A web server for whole-genome comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 47:W52-W58. https://doi.org/10.1093/nar/gkz333 Crossref, Medline, ISI, Google ScholarS. Lu and J. Bao contributed equally to this work.Funding: This work was supported by grants from National Natural Science Foundation of China (32160614), Hainan Provincial Natural Science Foundation of China (321CXTD437), the National Key R and D Program (2021YFC2600602), the Major Science and Technology Plan of Hainan Province (ZDKJ2021016), and the Scientific Research Foundation of Hainan University (KYQD(ZR)-20080).The author(s) declare no conflict of interest. The author(s) have dedicated the work to the public domain under the Creative Commons CC0 “No Rights Reserved” license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2022.DetailsFiguresLiterature CitedRelated Vol. 35, No. 10 October 2022ISSN:0894-0282e-ISSN:1943-7706 Download Metrics Article History Issue Date: 8 Oct 2022Published: 27 Sep 2022Accepted: 14 Jul 2022 Pages: 959-962 InformationThe author(s) have dedicated the work to the public domain under the Creative Commons CC0 “No Rights Reserved” license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2022.Funding National Natural Science Foundation of ChinaGrant/Award Number: 32160614 Hainan Provincial Natural Science Foundation of ChinaGrant/Award Number: 321CXTD437 National Key R and D ProgramGrant/Award Number: 2021YFC2600602 Major Science and Technology Plan of Hainan ProvinceGrant/Award Number: ZDKJ2021016 Scientific Research Foundation of Hainan UniversityGrant/Award Number: KYQD(ZR)-20080 Keywordsgenome assemblylong-read sequencingPhytophthora melonisThe author(s) declare no conflict of interest.PDF download