Complete mitochondrial genome and phylogeny of the causal agent of Bayoud disease on date palm, Fusarium oxysporum f. sp. albedinis

Abstract The complete mitogenome of Fusarium oxysporum f. sp. albedinis (FOA), the causal agent of the destructive fusarium wilt in date palm, is sequenced and assembled. The circular mitogenome of isolate Foa44 is 51,601 bp in length and contains 26 transfer RNA (tRNA) genes, one ribosomal RNA (rRNA), and 28 protein-coding genes. A mitogenome-based phylogenetic analysis of Fusarium revealed that FOA is congruent with previous nuclear-gene phylogenetic results.

Introduction Fusarium oxysporum f. sp. albedinis (Killian & Maire) Malenc¸on is a fungal pathogen causing Bayoud disease, also known as Fusarium wilt, on date palm (Phoenix dactylifera). This disease represents a major threat to the date palm industry in North African countries. In fact, Bayoud disease has decimated more than 10 million trees in Morocco over the last century (Djerbi 1982;Sedra M 2003;Sedra MH 2003a, 2003bSedra 2007). As a soil-borne pathogen, FOA spores and mycelium colonize date palm roots, spread internally through the vascular system, and cause external symptoms characterized by both external hemiplegia character and dried palm leaves having the appearance of wet feathers and ultimately resulting in date palm death (Bouhssini and Faleiro 2018). The nuclear genome of F. oxysporum f. sp. albedinis has been recently sequenced and assembled (Khayi et al. 2020 Total genomic DNA of Foa44 strain was extracted from freeze-dried mycelium using the cetyltrimethylammonium bromide (CTAB) method (M€ oller et al. 1992). A paired-end library was prepared using a Nextera DNA Flex library kit from total genomic DNA (0.5 lg), following the manufacturer's protocol. The library was sequenced (2 Â 150 bp) on a NovaSeq 6000 platform (Illumina, San Diego, CA). Adapters and low-quality reads were removed using CLC Genomics Workbench V12.
De novo assembly was performed using MaSurca V3.4 with default parameters (Zimin et al. 2013). The obtained assembly was subject to BLASTN searches against the reference mitochondrial genome sequence for Fusarium oxysporum f. sp. lycopersici (GenBank accession number CM010346). The BLASTN output identified a single homologous contig of 67,535 bp in length. This contained an $14 kbp flanking inverted repeat region, as determined by dot plot analysis in CLC genomics Workbench. The mitogenome sequence was circularized by fragmenting the contig in half and then assembling the two fragments using CLC Genomics Workbench. The raw reads were remapped on the assembled sequence to correct any conflicts created during the assembly step. The final length of FOA mitogenome is 51,601 bp with 31% in G þ C content. The mitogenome of Foa44 was deposited at GenBank under accession number MW493386. The annotation was performed using GeSeq (Tillich et al. 2017) and MFannot pipelines (https://github.com/BFL-lab/ Mfannot). The resulting annotations were manually inspected and curated in comparison to published F. oxysporum mitogenomes. In total, 28 protein-coding genes, one ribosomal RNA (rRNA), and 26 transfer RNA (tRNA) genes were predicted. Among 28 predicted genes, 14 genes are commonly found in MT genomes including those implicated in ATP production (atp6, atp8, and atp9), oxidative phosphorylation (nad1-6 and nad4L), apocytochrome b (cob), and cytochrome C oxidase subunits (cox1-3). The 14 remaining genes are of unknown functions except for four genes (orf529, orf304, orf292, and orf348) that are coding for Homing Endonucleases (GIY-YIG and LAGLIDADG) that are commonly found in fungal MT genomes (Megarioti and Kouvelis 2020).
To highlight the phylogenetic position of F. oxysporum f. sp. albedinis Foa44 within Fusarium oxysporum complex, a whole MT genome alignment of 16 Fusarium mitogenome sequences retrieved from Genbank, was performed with MAFTT software (Katoh and Standley 2013). The Maximumlikelihood phylogenetic tree was constructed using MEGA7, under a Kimura 2-parameter model (Kumar et al. 2016). The resulting phylogenetic tree shows that our mitogenome is clearly clustered with two other strains F. oxysporum f. sp. lycopersici 4287 (CM010346) and F. oxysporum F11 (NC017930) (Figure 1). Variant calling showed that there are 369 and 500 SNP/Indels differentiating the FOA mitogenome from these related mitogenomes, respectively.

Disclosure statement
No potential conflict of interest was reported by the author(s).

Data availability statement
The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at https://www.ncbi.nlm.nih.gov/ under the accession no. MW493386. The associated ÃÃ BioProject ÃÃ , ÃÃ BioSample ÃÃ , and ÃÃ SRA ÃÃ numbers are PRJNA658960, SAMN15893572, and SRP313423, respectively.

Funding
The author(s) reported there is no funding associated with the work featured in this article. , F. oxysporum f. sp. radices-cucumerinum (LT906343), and F. graminearum (NC_009493) as an outgroup). The tree with the highest log likelihood (À64,062.86) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach and then selecting the topology with a superior log-likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. All positions containing gaps and missing data were eliminated. There were a total of 27,505 positions in the final dataset. Evolutionary analyses were conducted in MEGA7 (Kumar et al. 2016).