GENOME SEQUENCES Draft Genome Sequence of Fusarium oxysporum f. sp. cubense Tropical Race 4 from Peru, Obtained by Nanopore and Illumina Hybrid Assembly Ana M. Leiva,a Mathieu Rouard,b Diana Lopez-Alvarez,a,c Alberto Cenci,b Catherine Breton,b Rosalyn Acuña,d Juan Carlos Rojas,e Miguel Dita,f Wilmer J. Cuellara aCrops for Nutrition and Health, International Center for Tropical Agriculture (CIAT), Palmira, Colombia bBioversity International, Montpellier, France cFacultad de Ciencias Agropecuarias, Universidad Nacional de Colombia, Palmira, Colombia dServicio Nacional de Sanidad Agraria (SENASA), Lima, Peru eInstituto Nacional de Innovación Agraria (INIA), Lima, Peru fBioversity International, Cali, Colombia Ana M. Leiva and Mathieu Rouard contributed equally to this article. ABSTRACT Fusarium oxysporum f. sp. cubense tropical race 4 (Foc TR4) is the causal agent of Fusarium wilt, a major threat to the banana industry worldwide. Here, we report the genome of a Foc TR4 strain from Peru, sequenced using a combination of Illumina and Oxford Nanopore Technologies. F usarium wilt of banana, caused by Fusarium oxysporum f. sp. cubense, is a devastatingfungal disease affecting bananas worldwide. The pathogen population is divided into four races, but tropical race 4 (Foc TR4) is by far the most aggressive, as it attacks several banana types, including Cavendish, which dominates the global banana export economy (1). In the last 6 years, Foc TR4 has spread from Asia into the Middle East and Africa (2), and in 2019, Foc TR4 reached Latin America, in the north of Colombia (3). In April 2021, banana plants (Musa acuminata group AAA, subgroup Cavendish) showing symptoms of Fusarium wilt were observed on a farm in Querecotillo, Peru (4°43954.840S, 80°33945.000W). Diagnostic analyses confirmed the identity of the pathogen as Foc TR4 (4). Pseudostem strands from symptomatic plants were transferred to potato dextrose TABLE 1 Representative subset of SIX gene homologs detected in isolates of Fusarium oxysporum f. sp. cubensea Presence of SIX gene: Race VCGb BRIP accession codec 1 2 4 6 7 8 9 13 R1 0123 62895 xd,f xb xb xa xa R2 01214 25609 xf xa,c xa Editor Antonis Rokas, Vanderbilt University R4 0122 62892 xc,i xa3 xa xc Copyright © 2022 Leiva et al. This is an open- STR4 0120 44012 xg xd xa xa xa3,b xa access article distributed under the terms of TR4 01213 40340 xa,h,i xa xc xa xa1,a2 xa a,e the Creative Commons Attribution 4.0 x International license. TR4d xa,h,i xa xc xa xa1 xa xa,e Address correspondence to Miguel Dita, a Isolates shown were reported in reference 12, with the addition of the TR4 Peruvian allelic variants. Sequences m.dita@cgiar.org, or Wilmer J. Cuellar, were searched on the assembly using BLAST matching with high similarity. X denotes the presence of a gene, w.cuellar@cgiar.org. while the superscript letters correspond to allelic variants of the gene (TR4 SIX1a, h, i: GenBank accession The authors declare no conflict of interest. numbers KX434991, KX434998, KX434999, respectively; SIX2a: KX435000; SIX4c: KX435006; SIX6a: SIX8a1: KX435011, KX435012; SIX9a: KX435015; SIX13: KX435019, KX435023). SIX, secreted in xylem. Received 11 April 2022 b VCG, vegetative compatibility group. Accepted 12 July 2022 c BRIP, Queensland Plant Pathology Herbarium. Published 8 August 2022 d Peruvian samples. September 2022 Volume 11 Issue 9 10.1128/mra.00347-22 1 Downloaded from https://journals.asm.org/journal/mra on 20 March 2023 by 2001:1388:541:2adc:5442:aa2f:fdbc:5905. Announcement Microbiology Resource Announcements FIG 1 Neighbor-joining phylogenetic tree of 17 Fusarium oxysporum f. sp. cubense TR4 isolates with 6,462 SNP variants. Support bootstrap values are indicated as a percentage of the 1,000 replicates. The sequences included were from Colombia (SRA accession number SRR10747097, SRR10125423, SRR10103605), India (SRR13311628), Indonesia (SRR10054446), Israel (SRR10054450), Jordan (SRR10054448), Laos (SRR7226878), Lebanon (SRR7226880), Pakistan (SRR7226883), Peru (SRR15514269 to SRR15514272), and the Philippines (SRR10054447) and were mapped onto strain UK0001 (GenBank accession number GCA_007994515.1). agar (PDA) medium and incubated at 25°C. Single-spore isolates from fungal colonies iden- tified as Fusarium oxysporum species complex were further purified and used for DNA extraction (5). DNA from 4 samples (PerS1 to PerS4) was extracted using the Illumina DNA prep kit and sequenced using the MiSeq platform (2  151 bp). The same DNA sample from PerS4 was further used for sequencing with Oxford Nanopore Technology (FLOW- MIN111, R10.3 chemistry, LSK109 kit) (6). A total of 46,707,802 Illumina and 379,956 Nanopore reads (average length, 2,783 bp) were obtained. The filtered reads (Illumina, 93.60% . Q30) were combined using Unicycler v0.4.8 (7) to make a hybrid genome assembly with a total length of 46,361,425 bp distributed in 115 contigs (G1C content, 47.59%; N50, 1.63 Mbp). The September 2022 Volume 11 Issue 9 10.1128/mra.00347-22 2 Downloaded from https://journals.asm.org/journal/mra on 20 March 2023 by 2001:1388:541:2adc:5442:aa2f:fdbc:5905. Announcement Microbiology Resource Announcements consensus sequence quality, checked using Qualimap v2.2.1 (8), resulted in an average depth of 19.25 with 6,718,534 reads mapped. The sequence assembly, mapped using QUAST v5.0.2 (9), showed high contiguity and a total aligned length of 45.9 Mb (94.7% genome fraction) with the highest-quality genome sequence available, strain UK0001 (10). Gene space assessment was performed using BUSCO v5.2.2 (11) (hypocreales odb10), which reported 97.7% completeness (single copy, 97.2%; duplicate, 0.5%; fragmented, 0.5%; missing, 1.8%; n = 4,494). The hybrid assembly was used to identify a family of secreted in xylem (SIX) genes. The presence or absence of the SIX homo- logs was checked using BLASTN v2.2.26 (Table 1) to identify the F. oxysporum f. sp. cubense sequences and matched the expected allelic variants identified in Foc TR4 for SIX1, SIX6, and SIX8 (12). Moreover, the PerS4 reads, combined with previously reported Foc TR4 strains (2, 3, 13–15), were mapped on UK0001 (10) using BWA v0.7.15 (16), and single nucleotide polymor- phism (SNP) calling was conducted using GATK v4.1.6 (17). A dissimilarity matrix (simple matching index) and a neighbor-joining phylogenetic tree were subsequently computed using Darwin v6 (18) (Fig. 1). Peruvian samples were clustered together, separated from other strains, including those from Colombia, suggesting independent incursions of Foc TR4 in the Americas (Fig. 1). Default parameters were used for all software unless otherwise specified. The availability of complete genome sequences and their comparative analysis will contribute to a better understanding of the Foc TR4 population biology, disease epide- miology, and management. Data availability. The sequence reads have been deposited at GenBank under the BioProject accession number PRJNA755905, and the assembly sequence has been de- posited under the GenBank accession number GCA_021237285.1. ACKNOWLEDGMENTS We acknowledge GenLab del Perú S.A.C. for their support with the Illumina sequencing. We thank Viviana Dominguez from the Cassava Crop Protection Team for technical support. Special thanks to the Plant Protection team of SENASA Peru. The genomic analysis was partially realized with the support of MESO@LR and the South Green Platform. This work was supported in part by grants from STC-CGIAR, Ministerio de Desarrollo Agrario y Riego (MIDAGRI) from Peru, the CGIAR Research Program on Roots, Tubers and Bananas (RTB) and CGIAR's Plant Health (PHI) initiative. REFERENCES 1. Ploetz RC. 2015. Management of Fusarium wilt of banana: a review with genome assemblies from short and long sequencing reads. PLoS Comput special reference to tropical race 4. Crop Prot 73:7–15. https://doi.org/10.1016/ Biol 13:e1005595. https://doi.org/10.1371/journal.pcbi.1005595. j.cropro.2015.01.007. 8. Okonechnikov K, Conesa A, García-Alcalde F. 2016. Qualimap 2: advanced 2. Zheng S-J, García-Bastidas FA, Li X, Zeng L, Bai T, Xu S, Yin K, Li H, Fu G, Yu Y, multi-sample quality control for high-throughput sequencing data. Bioin- Yang L, Nguyen HC, Douangboupha B, Khaing AA, Drenth A, Seidl MF, Meijer formatics 32:292–294. https://doi.org/10.1093/bioinformatics/btv566. HJG, Kema GHJ. 2018. New geographical insights of the latest expansion of Fu- 9. Gurevich A, Saveliev V, Vyahhi N, Tesler G. 2013. QUAST: quality assess- sarium oxysporum f. sp. cubense tropical race 4 into the Greater Mekong subre- ment tool for genome assemblies. Bioinformatics 29:1072–1075. https:// gion. Front Plant Sci 9:457. https://doi.org/10.3389/fpls.2018.00457. doi.org/10.1093/bioinformatics/btt086. 3. García-Bastidas FA, Quintero-Vargas JC, Ayala-Vasquez M, Schermer T, 10. Warmington RJ, Kay W, Jeffries A, O’Neill P, Farbos A, Moore K, Bebber DP, Seidl MF, Santos-Paiva M, Noguera AM, Aguilera-Galvez C, Wittenberg A, Studholme DJ. 2019. High-quality draft genome sequence of the causal Hofstede R, Sørensen A, Kema GHJ. 2020. First report of Fusarium wilt tropical agent of the current Panama disease epidemic. Microbiol Resour Announc 8: race 4 in Cavendish bananas caused by Fusarium odoratissimum in Colombia. e00904-19. https://doi.org/10.1128/MRA.00904-19. Plant Dis 104:994–994. https://doi.org/10.1094/PDIS-09-19-1922-PDN. 11. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. 4. Acuña R, Rouard M, Leiva AM, Marques C, Olortegui A, Ureta C, Cabrera- 2015. BUSCO: assessing genome assembly and annotation completeness Pintado RM, Rojas JC, Lopez-Alvarez D, Cenci A, Cuellar WJ, Dita M. 29 June with single-copy orthologs. Bioinformatics 31:3210–3212. https://doi.org/ 2022. First report of Fusarium oxysporum f. sp. cubense tropical race 4 causing 10.1093/bioinformatics/btv351. Fusarium wilt in Cavendish bananas in Peru. Plant Dis https://doi.org/10.1094/ 12. Czislowski E, Fraser-Smith S, Zander M, O'Neill WT, Meldrum RA, Tran- PDIS-09-21-1951-PDN. Nguyen LTT, Batley J, Aitken EAB. 2018. Investigation of the diversity of effector 5. Dita MA, Waalwijk C, Buddenhagen IW, Souza MT, Jr, Kema GHJ. 2010. A mo- genes in the banana pathogen, Fusarium oxysporum f. sp. cubense, reveals evi- lecular diagnostic for tropical race 4 of the banana fusarium wilt pathogen. dence of horizontal gene transfer. Mol Plant Pathol 19:1155–1171. https://doi Plant Pathol 59:348–357. https://doi.org/10.1111/j.1365-3059.2009.02221.x. .org/10.1111/mpp.12594. 6. López-Alvarez D, Leiva AM, Barrantes I, Pardo JM, Dominguez V, Cuellar 13. Thangavelu R, Edwin Raj E, Pushpakanth P, Loganathan M, Uma S. 2021. WJ. 2020. Complete genome sequence of the plant pathogen Ralstonia Draft genome of Fusarium oxysporum f. sp. cubense strain tropical race-4 solanacearum strain CIAT-078, isolated in Colombia, obtained using Oxford infecting Cavendish (AAA) group of banana in India. Plant Dis 105: Nanopore Technology. Microbiol Resour Announc 9:e00448-20. https://doi 481–483. https://doi.org/10.1094/PDIS-06-20-1170-A. .org/10.1128/MRA.00448-20. 14. Maryani N, Lombard L, Poerba YS, Subandiyah S, Crous PW, Kema GHJ. 2019. 7. Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycler: resolving bacterial Phylogeny and genetic diversity of the banana Fusarium wilt pathogen September 2022 Volume 11 Issue 9 10.1128/mra.00347-22 3 Downloaded from https://journals.asm.org/journal/mra on 20 March 2023 by 2001:1388:541:2adc:5442:aa2f:fdbc:5905. Announcement Microbiology Resource Announcements Fusarium oxysporum f. sp. cubense in the Indonesian centre of origin. Stud Wheeler transform. Bioinformatics 26:589–595. https://doi.org/10.1093/ Mycol 92:155–194. https://doi.org/10.1016/j.simyco.2018.06.003. bioinformatics/btp698. 15. Maymon M, Sela N, Shpatz U, Galpaz N, Freeman S. 2020. The origin and 17. McKenna A, HannaM, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella current situation of Fusarium oxysporum f. sp. cubense tropical race 4 in K, Altshuler D, Gabriel S, Daly M, DePristo MA. 2010. The Genome Analysis Tool- Israel and the Middle East. Sci Rep 10:1590. https://doi.org/10.1038/s41598 kit: a MapReduce framework for analyzing next-generation DNA sequencing -020-58378-9. data. Genome Res 20:1297–1303. https://doi.org/10.1101/gr.107524.110. 16. Li H, Durbin R. 2010. Fast and accurate long-read alignment with Burrows- 18. Perrier X, Jacquemoud-Collet JP. 2006. DARwin. http://darwin.cirad.fr/. September 2022 Volume 11 Issue 9 10.1128/mra.00347-22 4 Downloaded from https://journals.asm.org/journal/mra on 20 March 2023 by 2001:1388:541:2adc:5442:aa2f:fdbc:5905.