# Genome assembly ofMusa beccariishows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

Zheng-Feng Wang, Mathieu Rouard, Gaetan Droc, Xue-Jun Ge

► To cite this version:

Zheng-Feng Wang, Mathieu Rouard, Gaetan Droc, Xue-Jun Ge. Genome assembly ofMusa beccariishows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes. GigaScience, 2023, 12, 10.1093/gigascience/giad005. http://hal-inrae.fr/hal-04199202

HAL Id: hal-04199202

https://hal.inrae.fr/hal-04199202

Submitted on 7 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genome assembly of *Musa beccarii* shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

Zheng-Feng Wang<sup>1,2,3</sup>, Mathieu Rouard<sup>1,4</sup>, Gaetan Droc<sup>1,5,6</sup>, Pat (J. S.) Heslop-Harrison<sup>1,7,8</sup> and Xue-Jun Ge<sup>1,7,*</sup>

## Abstract

**Background:** *Musa beccarii* (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of *Musa* species is *x* = 7, 10, or 11; however, *M. beccarii* has a basic chromosome number of *x* = 9 (2*n* = 2*x* = 18), which is the same basic chromosome number of species in the sister genera *Ensete* and *Musella*. *Musa beccarii* is in the section *Callimusa*, which is sister to the section *Musa*. We generated a high-quality chromosome-scale genome assembly of *M. beccarii* to better understand the evolution and diversity of genomes within the family Musaceae.

**Findings:** The *M. beccarii* genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long iso- seq and short RNA-seq reads. The size of *M. beccarii* was the largest among all known Musaceae assemblies (~570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between *M. beccarii* and the other *Musa* and *Ensete* species, far beyond those expected from differ- ences in chromosome number. Within Musaceae, *M. beccarii* showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in *M. beccarii*. BGCs were not conserved in Musaceae genomes.

**Conclusions:** The genome assembly of *M. beccarii* is the first chromosome-scale genome assembly in the *Callimusa* section in *Musa*, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

**Keywords:** ancestral genome reconstruction, biosynthetic gene cluster, comparative genome, gene family, Musaceae, transcription factors, whole-genome duplication

## Introduction

Bananas are one of the most well-known and highly consumed fruits in the world. Phylogenetic studies of the genus *Musa* (family Musaceae) have shown that the genus comprises 2 sections: sect. *Musa* and sect. *Callimusa* [1–3]. The basic number of chromosomes in sect. *Musa* (c. 33–50 species) is *x* = 11 (wild accessions are 2*n* = 2*x* = 22) while the basic number of chromosomes in members of sect. *Callimusa* can be *x* = 7, *x* = 9, and *x* = 10 [1, 4]. The approx- imately 38 species in sect. *Callimusa* are 2*n* = 2*x* = 20, and lower numbers are found in *Musa ingens* (2*n* = 2*x* = 14) and *Musa beccarii* (2*n* = 2*x* = 18) [5]. Although the basic number of chromosomes of *M. beccarii* is unique among members within the genus *Musa*, *x* = 9 is shared among species in the 2 sister genera in the family: *Ensete* and *Musella* [5]. *Musa beccarii* is closely related to *Musa maclayi* and *Musa peekelii*, and these 3 taxa form a subclade sister to another

subclade including *Musa gracilis*, while *M. ingens* is sister to these 2 subclades.

*Musa beccarii* (NCBI:txid574481) is endemic in Borneo [5, 6]. Its leaves are long, narrow, bright green, and pest-free, and the verti- cal inflorescence has large, bright red bracts (Fig. 1). *Musa beccarii* begins flowering after 6 to 8 months, and the height of the plants ranges from 1 to 3 m. It is more compact than most other mem- bers of the family Musaceae, and it can be grown indoors as an ornamental plant [7]. The long-lasting bright red and attractive flowers [8] can be used as cut flowers. The conservation status of *M. beccarii* is currently “least concern” [9]; however, some have con- sidered this species to be endangered because of habitat loss [5, 6] and its small, isolated populations in the wild. It can be prop- agated by suckers; tissue culture has also been used in *M. beccarii* [7, 8], and this has aided its conservation.

<footnote n="1">¹Guangdong Provincial Key Laboratory of Applied Botany, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China</footnote>

<footnote n="2">²Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 511458, China</footnote>

<footnote n="3">³Key Laboratory of Vegetation Restoration and Management of Degraded Ecosystems, Key Laboratory of Carbon Sequestration in Terrestrial Ecosystem, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China</footnote>

<footnote n="4">⁴Biodiversity International, Parc Scientifique Agropolis II, 34397 Montpellier, France</footnote>

<footnote n="5">⁵CIRAD, UMR AGAP Institut, F-34398 Montpellier, France</footnote>

<footnote n="6">⁶UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, France</footnote>

<footnote n="7">⁷Key Laboratory of Plant Resources Conservation and Sustainable Utilization, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China</footnote>

<footnote n="8">⁸Department of Genetics and Genome Biology, University of Leicester, Leicester LE1 7RH, UK</footnote>

<footnote>*Correspondence address. Xue-Jun Ge, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou 510650, China. E-mail: xjge@scbg.ac.cn</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1:** Picture of a *Musa beccarii* flower.</figcaption></figure>

There are currently 12 fully assembled and annotated Musaceae genomes, including *Ensete glaucum*, *Musa acuminata*, *Musa balbisiana*, *Musa itinerans*, and *Musa schizocarpa*, according to "The Banana Genome Hub" [10]. Only half of these genomes have been assembled at the chromosome scale. The genome of *M. acuminata* was the first to be assembled with its "DH Pa-hang" genome sequence [11], and its genome was updated in 2021 [12]. The assembly revealed that 3 rounds of ancient whole-genome duplications (WGDs) have occurred in *Musa*. Following WGD, many genes involved in transcription regulation, signal transduction, and translational elongation were retained. Comparisons of published genomes have revealed that genes associated with transcription factors (TFs), defense-related proteins, enzymes involved in cell wall biosynthesis, and enzymes involved in secondary metabolism are *Musa* lineage specific. Subsequently, the genome of *M. balbisiana* was assembled [13], and it was updated using a double haploid [14]. Most edible banana cultivars are triploids derived from *M. acuminata* and *M. balbisiana* ancestors. Compared to *M. acuminata*, *M. balbisiana* shows more genome fractionation (gene loss) but contains more biotic and abiotic stress resistance properties [14]. *Musa itinerans* was the third species with its genome assembled [15]. *Musa itinerans* is a wild banana native to southeast Asia and one of the cold- and disease-resistant *Musa* species [15]. *Musa schizocarpa* was the fourth species in the

genus *Musa* to have its genome assembled [16]. It is native to Papua New Guinea, and a small proportion of its genome has introgressed into many cultivated edible bananas [17]. However, this genome is a draft assembly, and no in-depth comparative genomics has been conducted to data. Recently, a draft genome assembly of *Musa textilis*, an important fiber plant, was published [18]. However, this genome assembly is not appropriate for accurate comparative genomic analyses because it is fragmented and incomplete (e.g., only 78.2% of complete BUSCO genes were retrieved).

All previously assembled *Musa* genomes were from members of the section *Musa*; no genome assemblies have been generated from members of the sect. *Callimusa*. In the *Ensete* sister group, the first chromosome-scale genome assembly of *E. glaucum* (x = 9) was recently published by Wang et al. [19]. This genome assembly provided insights into the chromosome rearrangements and fusions that occurred between sister genera. Given that it has the same number of chromosomes as members of the genus *Ensete*, *M. beccarii* might have the most conserved genome structure with respect to the common ancestor between *Musa* and *Ensete* [5]; thus, a genome assembly of *M. beccarii* would provide an excellent resource for studies of genome evolution in the family Musaceae and enhance our knowledge of the pangenome and structural variants in the *Callimusa* section.
# Materials and Methods

## Sample collection and sequencing

One *M. beccarii* N. W. Simmonds individual, planted in a greenhouse in the South China Botanical Garden, Guangdong Province, China, was used for genome sequencing. The orientation of the greenhouse was from north to south. A fan was installed on the southern wall to lower the temperature in summer. The indoor temperature was maintained between 10°C and 35°C. The individual was a seedling approximately 50 cm in height and was cultivated in a plastic pot (diameter: 37 cm; height: 30 cm). The soil in the pot was a 1:1 mixture of Jiffy's TPS fine peat substrate (made from Estonian peat moss, pH 5.8) and sands; no fertilization was applied. The plant was automatically irrigated 2 times daily at 09:00 and 15:00. A sunshade net was placed on the roof of the greenhouse during the growing period. No other special treatments were applied. The plant was collected for genome sequencing between 14:00 and 15:00 on 21 September 2020.

Briefly, genomic DNA was extracted from fresh leaves using the cetyl trimethylammonium bromide method. Quality control was carried out using a NanoDrop 2000 microspectrophotometer (Thermo Fisher Scientific, Carlsbad, CA, USA), Qubit fluorometers (Thermo Fisher Scientific), and gel electrophoresis. High-quality DNA was used to build 1 short-read (Illumina, San Diego, CA, USA) and 2 long-read (Nanopore Oxford, UK and PacBio HiFi, San Diego, CA, USA) whole-genome sequencing (WGS) libraries. To perform Hi-C scaffolding, the genomic DNA was cross-linked with formaldehyde and extracted for Hi-C library preparation. Additionally, total RNA from *M. beccarii* leaves of the same individual was extracted and reverse-transcribed to complementary DNA (cDNA) for the construction of PacBio full-length cDNA sequencing (Iso-seq) and short-read cDNA fragment sequencing libraries; both libraries were used for genome annotation. HiFi WGS library construction and sequencing were conducted by Annoroad Gene Technology (AGT, Beijing, China), and the rest of the sequencing was conducted by GrandOmics Biosciences (GB, Wuhan, China).

The DNA sample was used to prepare a whole-genome shotgun paired-end (2 × 150 bp) Illumina library using the Truseq Nano DNA HT Sample Preparation Kit (Illumina). The library was sequenced on the Illumina HiSeq X Ten platform (RRID:SCR_01638 5). For Nanopore WGS library construction and sequencing, the DNA fragments were size-selected using the BluePippin system (Sage Science, Beverly, MA, USA). A sequencing library was prepared using size-selected fragments with the SQK-LSK109 Ligation Sequencing Kit (Oxford Nanopore Technologies, Oxford, UK), and sequencing was conducted using the Nanopore PromethION sequencer. The Hi-C library was prepared following the procedure described in a previous study [20] with some modifications. Briefly, fresh *M. beccarii* leaves were cut into 2-cm pieces and immersed in a nuclei isolation buffer with 2% formaldehyde for fixation. Vacuum infiltration was conducted for 20 minutes in this step. Glycine was added to stop fixation, and vacuum infiltration was conducted for another 15 minutes. Fixed tissue was rinsed in chilled water, dried on paper, and frozen in liquid nitrogen until nuclei were isolated. The isolated nuclei were digested with 100 units of DpnII (New England Biolabs, Ipswich, MA, USA), and then biotin was marked with biotin-14-dCTP. Extra Biotin-14-dCTP was removed, and ligation was conducted using T4 DNA polymerase (New England Biolabs). The ligated DNA was sheared into 300- to 600-bp fragments, blunt-end repaired, A-tailed, and purified using biotin-streptavidin-mediated pull-down. Finally, the Hi-C libraries were paired-end sequenced (2 × 150 bp) using the Illumina HiSeq X Ten or MGI DNBSEQ-T7 (RRID:SCR_017981) (MGI

Tech, Shenzhen, China) sequencing platforms. For HiFi WGS library construction and sequencing, a total of 50 µg extracted genomic DNA was sheared to approximately 10 kb using Covaris g-Tubes (Covaris, Woburn, MA, USA). The sheared DNA was purified and concentrated using AMPure PB magnetic beads (Cultek, Madrid, Spain). HiFi sequencing libraries were then prepared using Pacific Biosciences SMRTbell Template Prep Kit 1.0 (Pacific Biosciences, San Diego, CA, USA). The constructed library was further size-selected electrophoretically using SageELF systems from Sage Science. Primer annealing was then performed using the constructed library, and SMRTbell templates were bound to polymerases using the Sequel Binding Kit (Pacific Biosciences, San Diego, CA, USA). Finally, the Pacific Bioscience Sequel II platform (RRID:SCR_017990) was used for sequencing.

For RNA sequencing (RNA-seq) library construction and sequencing, total RNA was extracted from *M. beccarii* leaves using the TRNzol Universal RNA Extraction Kit (Tiangen, Beijing, China). A NanoDrop One UV-Vis spectrophotometer (Thermo Fisher Scientific), Qubit Fluorometer (Thermo Fisher Scientific), and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) were used to evaluate the quality and integrity of RNA. A TruSeq RNA Library Preparation Kit (Illumina) was then used to generate sequencing libraries following the manufacturer's instructions. The Illumina HiSeq X Ten platform was used with paired-end sequence (2 × 150 bp) cDNA libraries. For PacBio Iso-seq library construction and sequencing, total RNA was extracted, and the quality of the RNA was assessed using the method described above. RNA was then reverse transcribed into cDNA using the NEBNext Single Cell/Low Input cDNA Synthesis & Amplification Module (New England Biolabs) and the Iso-Seq Express Oligo Kit (Pacific Biosciences). The cDNA was then purified using the ProNex Beads; the SMRTbell Express Template Prep Kit 2.0 (Pacific Biosciences) was used to prepare the cDNA library. The Sequel Binding Kit was used to anneal sequencing primers to the SMRTbell templates and promote binding to polymerases. Sequencing was performed on the Pacific Biosciences Sequel II platform.

The libraries, sequencing, platforms, and data are summarized in Supplementary Table S1.

## Data preprocessing

After sequencing, Sickle v1.33 (RRID:SCR_006800) [21] was used to quality trim both short WGS and Hi-C reads by removing reads with base quality values less than 30 and lengths shorter than 80 bp. RECKONER v1.1 [22] was used to further error-correct short WGS reads. The CCS algorithm v6.0.0 (RRID:SCR_021174) [23] was used to process PacBio HiFi reads and obtain consensus reads. The error-corrected short WGS reads and/or HiFi reads were used to estimate the genome size of *M. beccarii* via KmerGenie v1.7044 [24], GenomeScope 2.0 (RRID:SCR_017014) [25], findGSE [26], GCE v1.0.2 [27], MGSE [28], and Gnodes [29]. Both MGSE and Gnodes provide mapping-based genome size estimations, while the others are k-mer based. PacBio Iso-seq reads were processed using IsoSeq v3.0 [30] to obtain full-length transcripts. Basecalling of raw Nanopore sequencing data (FAST5 format) was conducted using Guppy 3.2.10 (Oxford Nanopore Technologies) with default parameters to convert them to the FASTQ format. Porchop v0.2.4 [31] and HiFiAdapterFilt v1.0.0 [32] were used to remove the adapters in Nanopore and PacBio consensus long reads.

## Genome assembly

The initial assembly was optimized using different assemblers. NextDenovo v2.3.1 [33], Flye 2.9.1 (RRID:SCR_017016) [34], and
orthologous groups) that had potentially undergone expansions or contractions. When running CAFE, families that were not at the phylogenetic root were filtered.

For the above gene sets (family/genus/species specific and expanded/contracted), GO and KEGG enrichment analyses were performed using TBtools v1.098669 [78]. In the enrichment analysis, all the predicted genes with their GO/KEGG annotations were used as the background/reference gene set; the query gene set comprised the genes obtained from the above analyses (e.g., expanded/contracted genes). P values were obtained from hypergeometric tests and corrected using the Benjamini–Hochberg (BH) method. Significantly enriched (P <0.05 following BH adjustment) GO terms (the name of the term, i.e., biological process, among others) were further grouped and visualized using a treemap generated in REVIGO (RRID:SCR_005825) [79].

## Whole-genome duplication

Ancient WGD events in M. beccarii and the other 5 species in Musaceae were detected using wgd v1.2 [80]. Ksrates v1.1.1 [81] was used to characterize the timing of WGD events with respect to speciation events between M. beccarii and the other Musa species. Ksrates is based on the wgd package, but it rescales the syntenous nucleotide substitution (Ks) estimation by considering different Ks rates among lineages in a given phylogenetic tree, which permits more accurate inferences of speciation events. A simplified phylogeny obtained from the gene family analysis above was used in the ksrates analysis; only species in the family Musaceae were considered, and E. glaucum was used as an out-group species.

The DupGen_finder pipeline [82] was used to determine the number of duplications derived from WGD and other types of duplication events. Through searches of homologous gene pairs, DupGen_finder also identified possible tandem duplications (TDs), proximal duplications (PDs), transposed duplications (TRDs), and dispersed duplications (DSDs). TDs are defined as duplications in which the duplicated sequence is next to the original sequence (separated by 5 or fewer genes), PDs are defined as duplications in which the duplicated sequence is 10 or fewer genes away from the original sequence, TRDs are defined as transposable element mediated duplications, and DSDs are random nonneighboring duplications. GO and KEGG enrichment analyses were conducted using TBtools v1.098669 for genes derived from each type of duplication event. Significantly enriched GO terms were further grouped and visualized using a treemap generated in REVIGO if needed.

## Whole-genome alignment and synteny analysis

Syntenic blocks within the M. beccarii assembly and between the Musaceae assemblies were analyzed using MCScan (RRID:SCR_017650, Python version) in the jcvi package and visualized using the jcvi v1.1.19 [83] and Shinycircos [84] packages. The default parameters for the synteny analysis in MCScan were used, with the exception that the parameter "minimum number of anchors" was set to 10. MCScanX (RRID:SCR_022067, match score 3 and match size 10) [85] was also used, and results were visualized using SynVisio [86]. Dot plot alignments between Musaceae genome assemblies were generated and visualized using D-GENIES v1.2.0 (RRID:SCR_018967) [87].

## Biosynthetic gene clusters

Biosynthetic gene clusters (BGCs) in Musaceae species were identified by plantiSMASH v1.0 (Plant Secondary Metabolite Analysis Shell) [88]. The libraries used in PhytoClust [89] were also used

when running the plantiSMASH tool to enhance BGC identifications.

## Nucleotide-binding site-leucine-rich repeat gene identification

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes are the major plant resistance genes serving as an active defense against pathogens [90]. There are generally 3 main types of NBS-LRR genes [91]: Toll/interleukin-1 receptor NBS-LRR (TNL), N-terminal coiled-coil motif NBS-LRR (CNL), and resistance to powdery mildew 8 NBS-LRR (RNL) genes; TNL genes are absent in monocots [91, 92]. According to the results of InterPro/Pfam annotation in the tested Musaceae species, NBS-LRR genes were identified using the following protein domains: IPR03800, PF00931/IPR002182, PF13855/PF00560/IPR032675, and PF05659/IPR008808.

NBS-LRR genes were also detected using NLR-Annotator [93] with default settings. Rather than using annotated proteins predicted by gene models and transcriptomic data, NLR-Annotator directly uses genomic sequences to identify possible NLR genes, which were confirmed to be most efficient in NLR gene identification. After NLR genes were detected, NLR-Annotator categorized NLR genes as "complete," "complete (pseudogene)," "partial," or "partial (pseudogene)" according to the properties of each gene.

## Ancestral genome reconstruction

AnChro [94] was used to reconstruct ancestral genomes of Musaceae members with ginger (Zingiber officinale, GenBank accession number of GCA_018446,385.1) as an out-group. SynChro [95] was used to identify conserved syntenic blocks between different pairs of genomes; the blocks in the two genomes (with the shortest path connecting them in the phylogenetic tree) were then used to infer ancestral gene order by comparing them with the reference genomes. During syntenic block inferences, the stringency parameter, which determines the number of reciprocal best hits within a syntenic block, was set to 3.

## Results

The size of the M. beccarii genome was inferred to be 554,284,138 bp by KmerGenie under the best-selected k-mer size of 87 following comparisons of different k-mer spectra. Using different programs, the estimated genome size of M. beccarii ranged from 547,121,747 bp to 746,096,492 bp (Supplementary Table S5); the size of the genome was predicted to be between 565,341,681 bp and 661,920,459 bp according to the 2 mapping-based genome estimation programs, MGSE and Gnodes. The level of heterozygosity in the genome estimated by GenomeScope ranged from 0.287% to 0.815%.

The assembly sizes using different assemblers ranged from 607,623,222 bp (NextDenovo) to 816,255,026 bp (Canu) with Nanopore reads and 636,694,734 bp (Hifiasm) to 1,247,860,321 bp (HiCanu) with HiFi reads (Supplementary Table S6). The NextDenovo and Hifiasm assemblers displayed superior contig numbers, average and minimum lengths, and N50 values for Nanopore and HiFi reads, respectively; these 2 assemblies were used in subsequent steps of the genome assembly process. Details of the NextDenovo and Hifiasm assembly results are shown in Table 1.

After Hi-C read scaffolding, the assembly size was 569,617,942 bp with 449 scaffolds, the N50 value was 67,088,101 bp, and 551,683,906 bp (96.85%) of the sequences were assembled into 9 chromosomes (Table 1, Fig. 2A). The largest
<table><caption>**Table 1:** Statistics of the genome assembly of *Musa beccarii*</caption><thead><tr><th colspan="2">Contig statistics of initial assembly using Nanopore reads</th><th colspan="2">Contig statistics of initial assembly using PacBio HiFi reads</th><th colspan="2">Scaffold statistics after Hi-C scaffolding</th><th rowspan="2">Chromosome</th><th rowspan="2">Length</th></tr><tr><th>Length of the sequence (bp)</th><th>Order of the sequence length</th><th>Length of the sequence (bp)</th><th>Order of the sequence length</th><th>Length of the sequence (bp)</th><th>Order of the sequence length</th></tr></thead><tbody><tr><td>N10 = 48,080,317</td><td>L10 = 2</td><td>N10 = 8,700,004</td><td>L10 = 7</td><td>N10 = 79,885,826</td><td>L10 = 1</td><td>chr1</td><td>79,367,759</td></tr><tr><td>N20 = 39,700,570</td><td>L20 = 3</td><td>N20 = 5,180,184</td><td>L20 = 17</td><td>N20 = 79,367,759</td><td>L20 = 2</td><td>chr2</td><td>79,885,826</td></tr><tr><td>N30 = 27,992,656</td><td>L30 = 5</td><td>N30 = 3,933,706</td><td>L30 = 31</td><td>N30 = 73,517,995</td><td>L30 = 3</td><td>chr3</td><td>67,088,101</td></tr><tr><td>N40 = 21,895,089</td><td>L40 = 7</td><td>N40 = 3,192,498</td><td>L40 = 48</td><td>N40 = 73,517,995</td><td>L40 = 3</td><td>chr4</td><td>57,442,642</td></tr><tr><td>N50 = 18,949,966</td><td>L50 = 11</td><td>N50 = 2,546,178</td><td>L50 = 70</td><td>N50 = 67,088,101</td><td>L50 = 4</td><td>chr5</td><td>73,517,995</td></tr><tr><td>N60 = 15,652,116</td><td>L60 = 14</td><td>N60 = 2,007,927</td><td>L60 = 99</td><td>N60 = 60,040,564</td><td>L60 = 5</td><td>chr6</td><td>60,040,564</td></tr><tr><td>N70 = 12,145,816</td><td>L70 = 18</td><td>N70 = 1,507,786</td><td>L70 = 136</td><td>N70 = 57,442,642</td><td>L70 = 6</td><td>chr7</td><td>53,040,366</td></tr><tr><td>N80 = 8,091,256</td><td>L80 = 25</td><td>N80 = 1,059,843</td><td>L80 = 186</td><td>N80 = 53,040,366</td><td>L80 = 7</td><td>chr8</td><td>42,891,246</td></tr><tr><td>N90 = 1,849,914</td><td>L90 = 40</td><td>N90 = 527,812</td><td>L90 = 271</td><td>N90 = 42,891,246</td><td>L90 = 8</td><td>chr9</td><td>38,409,407</td></tr><tr><td>N100 = 21,817</td><td>L100 = 306</td><td>N100 = 12,368</td><td>L100 = 811</td><td>N100 = 1,000</td><td>L100 = 449</td><td></td><td></td></tr><tr><td>Total length, bp</td><td>607,623,222</td><td>636,694,734</td><td></td><td>569,617,942</td><td></td><td></td><td></td></tr><tr><td>Average length, bp</td><td>1,985,696.80</td><td>785,073.65</td><td></td><td>1,268,636.84</td><td></td><td></td><td></td></tr><tr><td>Largest length, bp</td><td>52,524,701</td><td>11,573,678</td><td></td><td>79,885,826</td><td></td><td></td><td></td></tr><tr><td>Minimum length, bp</td><td>21,817</td><td>12,368</td><td></td><td>1,000</td><td></td><td></td><td></td></tr></tbody></table>

chromosome (chr2) was 79,885,826 bp, and the shortest chromosome (chr9) was approximately 2 times smaller (38,409,407 bp) (Table 1).

BUSCO assessment of the final genome sequence assembly yielded a completeness score of 98.4% for the Embryophyta (1,614 core genes) datasets, including 1,510 (93.6%) complete and single-copy genes and 78 (4.8%) complete and duplicated genes. Another 11 genes (0.7%) were reported as fragmented, and 15 (0.9%) were reported as missing. The integrity of the genome sequence assembly was evaluated by mapping the Illumina WGS reads using BWA: 99.83% of the reads were mapped reads, and 95.26% of the reads were correctly paired.

The mapping quality metrics generated by SQUAT revealed 81.5% uniquely mapped reads (62.1% perfectly matched, 14.5% with substitution errors, 2.3% containing clips, and 2.6% with other errors), 18.3% multiply mapped reads, and 0.2% unmapped reads in the BWA-MEM mode; 79% uniquely mapped reads (62.5% perfectly matched, 14.9% with substitution errors, and 1.6% with other errors), 16.1% multiply mapped reads, and 4.9% unmapped reads were identified in the BWA-backtrack mode. The overall percentage of poorly mapped reads in the 2 modes was 2.5%.

## Annotation of repeat sequences

Repetitive regions comprised 51.79% (295,005,341 bp), and 51.45% (293,068,842 bp) of the genome assembly was identified according to EDTA and RED software. EDTA indicated that the most abundant repetitive sequences were long terminal repeat (LTR) retrotransposons, which accounted for 43.47% (247,628,340 bp) of the assembly, followed by terminal inverted repeats (TIRs), which accounted for 5.53% (31,478,316 bp) of the assembly (Supplementary Table S7). The largest proportion of LTR elements were Copia-like (144,383,969 bp, 25.35%) and Gypsy-like (51,691,527 bp, 9.07%) sequences.

A total of 318,946,703 bp (55.99%) of the genome sequence was annotated and masked as repetitive components when the results of EDTA and RED software were combined. The density of repeat sequences in the assembly is shown in Fig. 2B. Comparative analysis indicated that *M. beccarii* contained the highest number of repetitive sequences and the longest repetitive sequences (Fig. 3, Supplementary Table S7), which were mainly LTR retrotransposons and a small number of non-TIR helitrons.

The “seed” Nanica repetitive sequence from *M. acuminata* is 5,291 bp long. A BLAST search revealed 822 Nanica-like sequences in *M. beccarii* with lengths ranging from 55 to 3,891 bp, and 668 of them were longer than 1,000 bp. However, these sequences were not only concentrated in the centromere (Supplementary Fig. S2A). The “seed” sequence of Egcn repeats from *E. glaucum* was 134 bp; no similar sequences were detected in *M. beccarii* via BLAST searches or short/long WGS read-mapping of the assembly.

Three 45S rDNA repeats (18S, 5.8S, and 26S ribosomal RNA genes as well as intergenic spacers) were detected on chr5 (around bp 21,300,000), chr9 (around bp 38,200,000 near the telomere), and chr7 (around bp 44,400,000). The consensus monomer sequence was 10,402 bp long with a GC content of 60%. The consensus sequence included 17 copies of a tandem repeat (MuTR; GenBank AM905874 to AM905898), although the number of MuTR repeats varied between rDNA monomers in the Nanopore long-molecule reads. Excluding MuTR repeats, present at multiple genome sequence sites, the 45S rDNA repeat was present in 5.0% of the examined short sequence reads. The 5S rDNA monomer was 432 bp long with a GC content of 55.8%, and it was present in 0.11% of the reads; the major sites were located on chr8 at bp 16,100,000 and chr7 at bp 37,424,000. Peaks associated with the higher GC content of rDNA sequences (assembly average 38.7% GC) are shown in the GC content plot (Fig. 2B).

## Gene prediction and annotation

A total of 39,112 genes encoding 45,461 proteins were predicted in *M. beccarii*. Of these genes, 38,756 (85.25%) were functionally annotated (Supplementary Table S8) with a completeness score of 94.8% in the embryophyta_odb10 database, according to the BUSCO assessment.

Alternative splicing events were detected in 4,602 genes. Retained introns were the most frequent events (2,847), followed by alternative 3' splice sites events (830). Other detected events include skipping exon events (313), alternative 5' splice sites (424), mutually exclusive exons (8), alternative first exons (109), and alternative last exons (103).

A total of 3,168 genes encoding TF in *M. beccarii* were identified, and this was similar to the number of TF genes identified in the genomes of *M. acuminata* and *E. glaucum* (Fig. 2C; Supplementary Table S9). Among these genes, MYB genes were the most
<figure><img src="image_3.png" /><figcaption>**Figure 2:** (A) Hi-C interaction heatmap (bin length 10,000 bp) for the *Musa beccarii* genome. (B) Genome features across *M. beccarii* chromosomes. (C) Inferred phylogenetic tree and contracted (-) and expanded (+) gene families in *M. beccarii* and other species in Liliopsida. Gene families within the most recent common ancestor are denoted at the root. Numbers following each species indicate the statistics of different genes</figcaption></figure>

abundant in *M. beccarii* and in the other Musaceae species (Supplementary Table S10). Using the MYB_annotator, a total of 292 MYB genes were identified in *M. acuminata*; this is similar to 294 previously reported [96]; 268 MYB genes were identified in the *M. beccarii* genome (Supplementary Tables S10 and S11). MYB genes associated with "axillary meristem, root growth"; "cell wall, lignin, seed oil, axillary meristem"; "defense, stress response"; "repres- sor phenylpropanoid, sinapate, lignin"; and "stress response, hor- mone signaling" were more abundant in the Musaceae genomes than in the genomes of other species (Supplementary Table S12). The 3 anthocyanin genes, MB_008808-T1, MB_018229-T1, and MB_003891-T1, in *M. beccarii* were orthologous to *MusaMYB-α*, -β, and -γ in *M. acuminata*, and this is associated with the transcrip- tional activation of anthocyanin biosynthesis in banana [97].

## Orthogroup identification and gene enrichment

A total of 32,123 orthogroups were identified from a set of 495,640 genes from selected monocots species. In the *M. beccarii* genome,

83.90% (32,815/39,112) of the genes were assigned to 50.03% (16,070/32,123) of the gene families, and 248 gene families com- posed of 671 genes were specific to *M. beccarii* (Supplementary Ta- ble S13).

A total of 7,810 gene families, which included 3,531 genes in *M. beccarii*, were specific to Musaceae. GO and KEGG enrich- ment analysis revealed that these Musaceae-specific genes in *M. beccarii* were involved in the regulation of protein modifica- tion, transcription, and cell wall in the GO Biological Process (BP) category (Supplementary Table S14, Supplementary Fig. S3) and flavonoid biosynthesis, tryptophan metabolism, and phenyl- propanoid biosynthesis in KEGG (Supplementary Table S15).

The 5 *Musa* species shared 22,000 gene families; 11,136 of these genes were shared among all 5 *Musa* species, and these genes were considered core families (Supplementary Fig. S4). There were 25,754 *M. beccarii* genes in the core families. The GO analysis re- vealed that these genes were mainly involved in cellular and metabolic processes in the BP category and binding and catalytic
<figure><img src="image_4.png" /><figcaption>**Figure 3:** (A) Number and (B) length of transposable elements in Musaceae species. LTR: long terminal repeat; TIR: terminal inverted repeat; nonTIR: non-terminal inverted repeats.</figcaption></figure>

activity in the Molecular Function (MF) category (Supplementary Fig. S5). Transcription regulator and transporter activities were the other 2 main functions in the MF category. A total of 1,062 gene families comprising 1,617 genes were specific to the *M. beccarii* genome (Supplementary Fig. S4).

A phylogenetic tree (Fig. 2C) revealed that *M. beccarii* had an estimated divergence time from the other *Musa* species of approximately 25.26 (95% CI: 8.25–54.84) million years ago. A total of 12,211 gene families were used in the gene family expansion and contraction analysis. In the *M. beccarii* genome, 1,518 gene families have expanded, and 885 gene families have contracted; 84 of these expansions were significant (P < 0.05), and 50 of these contractions were significant. Significantly expanded gene families

were mainly involved in transcription, carbohydrate metabolism, and membrane transport (Supplementary Tables S16 and S17, Supplementary Fig. S6). Significantly contracted gene families were mainly involved in defense response according to the GO analysis and (mono)terpenoid biosynthesis and translation factors according to the KEGG analysis (Supplementary Tables S18 and S19).

Genes with alternative splicing were mainly involved in messenger RNA 3'-end processing, amino acid catabolic processes, phosphorus metabolic processes, response to stress (such as DNA repair), and taurine and hypotaurine metabolism according to the GO and KEGG analyses (Supplementary Tables S20 and S21, Supplementary Fig. S7)
## Gene duplications

All Musaceae species have undergone the same 3 ancient WGD events (Fig. 4A), and the 5 *Musa* species diverged following the WGD events (Fig. 4B).

Gene duplication analyses in *M. beccarii* revealed 11,244 gene pairs that possibly derived from WGDs, 531 pairs derived from tandem duplications, 646 pairs derived from proximal duplications, 2,313 pairs derived from transposed duplications, and 7,690 pairs derived from dispersed duplications. The enrichment analysis revealed that duplicated genes derived from WGDs were mainly involved in transcription, signaling, defense, environment adaptation, and root development (Supplementary Tables S22 and S23, Supplementary Figs. S8 and S9A). Genes derived from tandem duplications were mainly involved in various metabolic processes related to stress responses (e.g., glutathione and phenylpropanoid metabolism) and defense (e.g., cell wall formation and membrane transport) (Supplementary Tables S24 and S25, Supplementary Figs. S10 and S9B). Genes derived from proximal duplications were mainly involved in benzoxazinoid, terpenoid, and flavonoid biosynthesis; membrane transport; and cell wall formation (Supplementary Tables S26 and S27, Supplementary Fig. S9C). Genes derived from transposed duplications were mainly involved in ion transport (Supplementary Table S28). Genes derived from dispersed duplications were mainly involved in DNA repair, monosaccharide metabolic process, and prokaryotic defense system (Supplementary Tables S29 and S30, Supplementary Figs. S11 and S9D).

## Whole-genome alignment and syntenic analysis

Overall, syntenic relationships of *M. beccarii* with the other Musaceae assemblies were observed, including several major syntenic blocks of genes with extensive rearrangements, including fusions, fissions, and translocations (Fig. 4C and Fig. 5). Chr4 of *M. beccarii* showed the highest conserved relationship with chr4 of other *Musa* species (Fig. 5). Chr5, which was the only one conserved between *E. glaucum* and *M. acuminata*, was divided into chr3 and chr5 in *M. beccarii* (Supplementary Fig. S12); this indicates that this fission is specific to *M. beccarii*. The other chromosomes of *M. beccarii* have undergone various rearrangements; for example, chr3 is syntenic with 5 chromosomes of both *E. glaucum* and *M. acuminata* (Supplementary Fig. S12).

Syntenic analysis of the ancient WGDs within *M. beccarii* revealed 233 syntenic blocks containing 9,594 genes and 5,512 gene pairs. The size of the longest syntenic block was 8,645,195 bp, and it contained 44 gene pairs between chr5 and chr6; the smallest syntenic block was 161,642 bp, and it contained 15 gene pairs between chr2 and chr3 (Supplementary Table S31). The syntenic relationships are shown in the CIRCOS plot (Fig. 2B).

A total of 196, 111, 155, and 141 syntenic blocks in *M. beccarii* were shared with *E. glaucum*, *M. balbisiana*, *M. acuminata*, and *M. schizocarpa*, respectively. The largest blocks occurred in the chr4s of *M. beccarii* and the other *Musa* species (Fig. 4C and Fig. 5). These largest blocks contained 1,776, 2,495, and 1,675 gene pairs between *M. beccarii* and *M. balbisiana*, *M. acuminata*, and *M. schizocarpa*, respectively, including a total of 2,602 genes in *M. beccarii* chr4.

## Biosynthetic gene clusters

*M. beccarii* contained 66 possible BGCs (Table 2 and Supplementary Table S32), the second largest number of BGCs in sequenced genomes of Musaceae species after *M. acuminata* (72 clusters). The most abundant BGCs in all Musaceae species were similar to type III polyketide synthase (T3PKS) and tomatine clusters. The BGCs

in chr4 of *M. beccarii* and their syntenic genes in the chr4s in the other *Musa* species are shown in Fig. 4D; these findings indicate that BGCs are not conserved in *Musa*. Substantial gains and losses of BGC and genes were observed.

## Identification of NBS-LRR genes

Using the annotated protein sequences in Musaceae species, the highest number of CNL genes was detected in *M. itinerans* (59 genes), and the lowest number of CNL genes was detected in *M. balbisiana* (14). A total of 31 CNL genes were detected in *M. beccarii* (Fig. 2C). Only 1 RNL gene was detected in each Musaceae species.

NLR-Annotator was used to identify the most complete and highest overall number of NBS-LRR genes in *M. beccarii* (Supplementary Fig. S2B, Supplementary Table S33); the most NBS-LRR genes were observed in chr6, and no NBS-LRR genes were identified in chr5 (Supplementary Fig. S2B). Only 74 NBS-LRR genes were identified in *M. beccarii* when the genes were detected using predicted genes (Supplementary Table S34), much less than the 179 genes obtained from the genome sequence. The number of NBS-LRR genes detected using predicted genes was also low in *M. balbisiana* compared with the number of genes identified from its genome sequence (43/96). In *M. acuminata*, which has the most complete genome sequence, the numbers of NBS-LRR genes identified from predicted genes and the genome sequence were 111 and 128. In *M. itinerans*, the numbers were 149 and 138, most likely stems from the fragmented draft assembly.

## Ancestral genome reconstruction

The ancestral genome reconstruction revealed 86 contigs in the genome of the last common ancestor (LCA) and between 19 and 40 contigs in the genome of the intermediate ancestors in Musaceae (Fig. 6). Although these ancestral genomes are fragmented, the complex chromosomal rearrangements that occurred between *Ensete*, *Musa* sect. *Musa*, and *Musa* sect. *Callimusa* are evident in these genomes. The number of macro- and micro-rearrangements in these Musaceae species was high; the number of contig rearrangements ranged from 18 for *M. acuminata* to 91 for *M. beccarii*, which is consistent with the phylogeny.

## Discussion

### Genome size

The chromosome-scale assembly of *M. beccarii* identified 9 pseudomolecules between 38 and 79 Mbp long (Fig. 2, Table 1), which is more variable than those in other *Musa* species (e.g., between 35 and 51 Mbp in *M. acuminata*) [12]. *Musa beccarii* has the largest genome (assembly size ≈ 570 Mbp, consistent with various estimates from k-mers) among species within the family Musaceae [5, 98]. The genome size of *M. beccarii* was estimated to be between 764 and 804 Mb according to DNA flow cytometry [5, 99], which might be an overestimation compared with the k-mer/mapping-based genome size estimation. Similar to other *Musa* species, the size of the genome estimated by flow cytometry [98–100] was larger than the genome assembly (534–578 vs. ~457 Mbp in *M. balbisiana*, 591–646 vs. ~469 Mbp in *M. acuminata*, and 704 vs. 515 Mbp in *M. schizocarpa*). Variation in the size of the genome estimated among methods might be affected by many factors, including the accuracy of flow cytometry, reference values used based on chemical measurements [101–103], variation among samples [100, 104], reference genome staining, and errors in genome assembly. Regardless of the method used to estimate the genome size, the data indicate that *M. beccarii* has the largest genome size among
<figure><img src="image_5.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_8.png" /></figure>

**Figure 4:** (A) Density distribution of synonymous nucleotide substitutions (Ks) in whole-genome duplication analysis. (B) Speciation event (red line) detection using the rate-adjusted Ks distribution for *Musa beccarii* with the ksrates package. The background was the whole-paranome Ks distribution (light gray histogram and KDE curve) and anchor-pair Ks distribution (dark gray histogram and KDE curve) for *M. beccarii*. The shared number in the red circle indicates the same speciation event between *M. beccarii* and the other *Musa* species. The numbers and arrows in the parentheses of the 4 *Musa* species in the panel legend indicate Ks value shifts after substitution rate adjustments by ksrates. (C) Syntenic blocks between Musaceae species. The largest blocks in *Musa* are highlighted in orange. (D) Biosynthetic gene clusters (BGCs) in chr4 in *M. beccarii* and gene synteny with the other 3 *Musa* species in their chr4s. The regulatory genes in the BGCs are not shown.

*Musa* species [5, 98]. The relationship between genome sizes in Musaceae species and their biological functions requires further study. Furthermore, variation in the structure of the genome between *M. beccarii* and *M. acuminata* indicates that the large genome of *M. beccarii* stemmed from changes in several chromosomes (Supplementary Table S35).

Repetitive sequences, especially transposable elements (TEs), are important elements driving genome expansion [105–107], at least in species with genomes smaller than 5 Gbp [108]. We detected a clear increase in the number and length of TEs in *M. beccarii*, which is approximately 30 to 106 Mbp larger than those of *E. glaucum* and the other 3 *Musa* species (Supplementary Table S7; Fig. 3). Previous studies using low-coverage sequencing have shown that *M. beccarii* contained the highest repetitive sequences among the 5 tested *Musa* species, including *M. balbisiana* and *M. acuminata* [98]; this suggests that increases in repetitive sequences might contribute to the larger genome size of *M. beccarii*. Nevertheless, for each TE, only LTR-unknown and helitron were consistently more abundant in the *M. beccarii* genome (51,552,844 bp and 15,898,685 bp, respectively) than in all genomes examined (Supplementary Table S7). Using LTR markers, Häkkinen et al. [6] identified rich and distinct LTRs in *M. beccarii*, suggesting diversification of LTRs, resulting in unknown LTRs in *M. beccarii*.

The assembly revealed the presence of 3 pairs of 45S rDNA loci, compared with only 1 in the other *Musa* and *Ensete* species assemblies, which is consistent with the *in situ* hybridization results of Bartoš et al. [99]. The 45S rDNA on 3 chromosomes accounted for 5.0% of the Illumina sequence reads in *M. beccarii* compared with 1.2% in *E. glaucum* on 1 chromosome [19], so increased rDNA copy number is responsible for some of the increase in genome sequence size.

## Gene family evolution

Gene family expansion due to duplications in *Musa*, including *M. beccarii*, was mainly caused by ancient WGD events. TFs are abundant in Musaceae species and higher than all the other monocots used in our dataset, with the exception of *Zingiber officinale* (Fig. 2C, Supplementary Table S9), which is a tetraploid species. The enrichment analysis of *M. beccarii* revealed that genes encoding TFs are some of the main genes retained following WGD events (Supplementary Table S23; Supplementary Figs. S8 and S9A). The fact that Musaceae species have experienced the same WGD events indicates that the retained genes encoding TFs play key roles in mediating adaptation to stress by TFs in Musaceae, which is consistent with results of previous studies of the genomes of *M. acuminata* and *M. itinerans* [11, 15]. To identify possible
<figure><img src="image_9.png" /><figcaption>**Figure 5:** Dot plots of *Musa beccarii* and the other 4 species in Musaceae made using D-GENIES. The dot colors correspond to the similarity values, which are binned into 4 groups. Highly conserved chromosomes between *M. beccarii* and the other *Musa* genomes are highlighted in red boxes.</figcaption></figure>

duplications caused by polyploidies or aneuploidies, the mapping results of Illumina WGS reads in the *M. beccarii* assembly were examined, including the mapping coverage distribution along the chromosomes and the allele frequencies of single-nucleotide variants according to Busche et al. [109]. The mapping coverages were constant along all chromosomes (Supplementary Fig. S13), and a peak allele frequency of 0.5 was observed in all the chromosomes (Supplementary Fig. S14), suggesting that no large segmental duplications have occurred in diploid *M. beccarii*.

We detected a contraction of gene families containing genes involved in defense response, monoterpenoid biosynthesis, and terpenoid backbone biosynthesis (Fig. 2C, and *M. beccarii* had fewer genes than other Musaceae species; Supplementary Tables S18 and S19). Terpenoids are important natural products [110–112]. They comprise diverse components and have various applications, particularly for defense [111], acting as toxic com-

pounds against biological stress agents in plants. Aside from terpenoid backbone and monoterpenoid biosynthesis, according to KEGG, there are other terpenoid biosynthesis-related pathways, such as steroid biosynthesis (ko00100), ubiquinone and other terpenoid-quinone biosynthesis (ko00130), limonene and pinene degradation (ko00903), diterpenoid biosynthesis (ko00904), brassinosteroid biosynthesis (ko00905), carotenoid biosynthesis (ko00906), zeatin biosynthesis (ko00908), and sesquiterpenoid and triterpenoid biosynthesis (ko00909). No gene families in these pathways have undergone contraction. Because typical terpenoid synthase genes are characterized by 2 conserved domains with Pfam ID PF01397 and PF03936 [112], the comparison of the genes with these domains in Musaceae revealed that the number of primary terpenoid synthase genes in *M. beccarii* is low (Supplementary Table S9, Fig. 2C). Gene families that have undergone significant expansion in *M. beccarii*, which are mostly
<figure><img src="image_10.png" /><figcaption>**Figure 6:** The chromosomal history of the Musaceae genomes shows changes in the genome structure from the last common ancestor (LCA, the topmost one above the genome in the picture) to 3 intermediate ancestors and 5 studied species. The genes in the LCA are represented with lines and the same colors if they are in the same contigs, with the exception of genes in the contigs containing fewer than 100 genes, which are stacked into 1 super contig and all colored in black. The genes in the intermediate ancestors and in the studied species are colored based on orthologous genes in the LCA (determined by reciprocal best hits in SynChro analysis); otherwise, they are white to indicate the lack of homology with genes in the LCA.</figcaption></figure>

involved in transcription, carbohydrate metabolism, and membrane transport (Supplementary Tables S16, S17, Supplementary Fig. S6), mediate various processes related to plant growth, development, and defenses [113–116]. Therefore, genes in expanded families might help balance the growth, development, and defense in *M. beccarii*. We also examined NBS-LRR genes in *M. beccarii* and other species in the phylogeny (Fig. 2C). These genes in *M. beccarii* were more abundant (179 genes according to sensitive NLR-Annotator analysis [93] in the assembly) than in other Musaceae (67–138 genes), although they were less abundant than in most species out of Musaceae (Fig. 2C, Supplementary Table S9). Therefore, combing Pfam annotation indicated that NBS-LRR genes might not be a priority in disease defense in the family Musaceae.

## The cell wall as a defense (lipid metabolism and adenosine triphosphate-binding cassette transporters)

Although they do not play a direct role in defense, genes encoding proteins that mediate the synthesis of substances such as cutin, suberin, and wax play a critical role in “physical defense” in plants [117–119]. They are lipids [117, 120] formed by fatty acids and glycerol. They comprise the extracellular hydrophobic layer of cell walls in plants, provide mechanical support, and protect plants from desiccation, extreme temperatures, UV light, and attack by pathogens and pests [117, 119]. According to the KEGG PATHWAY database [121], lipid metabolism comprises 16 pathways, and genes involved in 14 of these pathways were identified in our study species (Supplementary Table S9). Comparative analyses revealed that Musaceae species do not contain more

genes involved in these pathways than other species, and the numbers of genes involved in fatty acid elongation (ko00062); cutin, suberin, and wax biosynthesis (ko00073); glycerophospholipid metabolism (ko00564); and ether lipid metabolism (ko00565) pathways were only slightly higher in the genome of *M. beccarii* than in the genomes of other Musaceae species (Fig. 2C).

Cell wall lipids need to be exported to the plant surface to be synthesized in epidermal cells. Adenosine triphosphate (ATP)-binding cassette (ABC) transporters are essentially required [117, 122]. G family ABC transporters are responsible for the secretion of lipids [123]. We identified a large number of genes in the G family of ABC transporters in all Musaceae species (Supplementary Table S36, both InterPro and eggNOG annotations). The number of ABC transporter genes was highest in Musaceae species according to the eggNOG annotation; however, this was not observed when other annotation pipelines were used. *M. beccarii* had the highest number of ABC transporter genes according to the InterPro annotation. ABC transporters are one of the largest protein families in nature [124]. They bind and hydrolyze ATP and mediate cellular transport processes [124, 125]. They transport molecules such as ions, amino acids, sugars, lipids, peptides, proteins, and antibiotics. The expansion of ABC transporters has been shown to be associated with increases in defenses against abiotic and biotic stress in plants, which promotes adaptation [124, 126, 127].

## Flavonoid biosynthesis as a defense

Flavonoids are a ubiquitous group of polyphenolic compounds in plants. They are important secondary metabolites that have been studied extensively from their biosynthesis to their biological activities [128–131], also in Musaceae species [132, 133].
However, no comparative genomic studies of flavonoids have been conducted among Musaceae species. We noticed significant gene enrichment in flavonoid biosynthesis in M. beccarii (Supplementary Table S15). Genes involved in flavonoid biosynthesis are high in Musaceae species, and the genome of M. beccarii contained the second highest number of flavonoid biosynthesis genes (Supplementary Table S9, Fig. 2C). The high number of flavonoid biosynthesis genes in Musaceae species is majorly derived from the genes encoding naringenin 7-O-methyltransferase (NOMT; KEGG Orthology term K22440). NOMT can catalyze the methylation of naringenin to produce sakuranetin (Supplementary Fig. S15), a phytoalexin with strong antifungal activity [134]. Therefore, the accumulation of NOMT genes involved in flavonoid biosynthesis suggests that flavonoids are functionally important in the resistance to disease in Musaceae and are potentially valuable components of harvested Musa crop [135]. ABC transporters are key mediators of flavonoid transport in plants [136, 137]. Both flavonoid biosynthesis and ABC transporter-related genes were enriched in tandem and proximal duplications in M. beccarii (Supplementary Tables S25 and S27). These 2 duplicates promote the evolution of self-defense in plants [82]. Future studies are needed to evaluate the presence/absence of flavonoid biosynthesis-related genes among Musaceae and other species, as well as the expression patterns and the biological/phenotypic effects of these genes, as such studies will enhance our understanding of defense mechanisms and other metabolic processes in Musaceae.

Musa beccarii has bright red flowers (Fig. 1). Anthocyanins are important substances that affect flower colors [128, 138], and their biosynthesis is closely associated with the flavonoid biosynthesis pathway. Anthocyanins have been shown to affect the red peel of Musa AAA Cavendish cv. Baxi [139] and the purple peel of M. itinerans [140] fruits. We found that Musaceae species show little variation in the number of genes involved in anthocyanin biosynthesis compared with the other species (Supplementary Table S9). Given that flavonoid/anthocyanin biosynthesis is mostly regulated at the transcriptional level [132, 141, 142], transcriptome comparisons among different tissues and species are needed to clarify the formation of the red color in the flowers of M. beccarii.

## Biosynthetic gene clusters

We identified diverse BGCs in Musaceae. BGCs are nonrandomly ordered genes along chromosomes that may optimize the synthesis of natural products in living organisms [143, 144]. The most developed cluster in Musaceae BGCs comprises T3PKSs (Table 2, Supplementary Table S32). T3PKSs are homodimer ketosynthases widely distributed in plants, fungi, and bacteria [145]. They take part in various important biosynthesis of secondary metabolites related to polyketides, facilitate the production of various natural products [146, 147], and play a role in defense responses and development [147–149]. In Musa, T3PKSs can initiate phenylphenalenone biosynthesis; phenylphenalenones are major phytoalexins involved in the defense against multiple pathogens in Musa [149]. The T3PKS BGCs in Musa merit further investigation because of their potential for enhancing defense systems.

Tomatine alike clusters were also abundant BGCs in Musa species. Tomatine is a steroidal glycoalkaloid saponin in tomatoes and other Solanum species [150, 151]. It has antipathogen and antiherbivore properties and serves as a natural defense in plants [151–153]. The primary genes involved in tomatine biosynthesis are glycosyltransferases (GTs) [154, 155]. GTs mediate the glycosylation of tomatidine, which is phytotoxic and a steroidal alkaloid (SA), to promote tomatine formation, and they reduce the toxicity

of SA metabolites to the plant cell [155–157]. Although a tomatine alike steroidal saponin that has been shown to promote resistance to black Sigatoka has been reported in M. acuminata [157], the BGC containing the genes that mediate its synthesis remains poorly resolved [158]. The BGCs characterized in the current study provide valuable resources for future studies of the biosynthesis of this saponin in Musa.

GTs are ubiquitous enzymes that are involved in the synthesis of various secondary metabolites in plants [159]. They generally function by glycosylating substrates with sugar moieties attached to aglycones, which then form glycosidic bonds. Their acceptor substrates can be sugars, lipids, proteins, nucleic acids, antibiotics, or small molecules [160, 161]. Glycosylation is rendered highly diverse by using various sugar moieties, and it can play diverse roles in plant growth, development, and defense responses [162]. In M. beccarii, we detected highly abundant GT-related genes (Supplementary Table S37, Fig. 2C), which might underlie its ability to adapt to environmental conditions. Furthermore, GT family 61 (GT61) genes encode proteins involved in xylan biosynthesis for the cell wall [163, 164]. Xylans are hemicelluloses that can affect cell wall recalcitrance and aid defense against herbivores and pathogens [165]. The number of GT61 genes was highest in M. beccarii among Musaceae species (Supplementary Table S37, Fig. 2C). DupGen_finder revealed that 32 (66.7%) of these genes were derived from ancient WGDs; this finding is consistent with the results of a previous study showing that the major duplication of GT families stemmed from WGD events [164, 166]. Experimental studies are needed to fully characterize BGCs.

## Chromosomal rearrangements

We detected substantial numbers of chromosomal reorganization events involving chromosomal fusion and fission in M. beccarii and the other 3 Musa species (Figs. 4C and 6), and only 1 chromosome remained largely intact. The extensive fusion/fission events between M. beccarii in the Musa section *Callimusa* and the other Musa section *Musa* species are similar in number to those detected between *Musa* and *Ensete* (Figs. 4C, 5, and 6 and Supplementary Fig. S12) and are not only a consequence of the reduced number of chromosomes (x = 9 vs. x = 11). This finding strongly supports the division of M. beccarii and the other studied *Musa* species into 2 different sections. Chromosomal fusion and fission are important mechanisms of speciation [167–169]. However, our current ancestral reconstructions did not permit the numbers of ancestral chromosomes between *Musa* sections or between *Ensete* species to be inferred.

The high abundance and expansion of TEs might facilitate evolutionary genome rearrangements in M. beccarii and be responsible for the large structural differences between M. beccarii and other studied Musa species. Structural rearrangements mediated by various families of TE elements have been reported in other plants [170, 171]; chromosome-scale assemblies anchored by long-molecule sequencing will enable further study of the association of TEs to chromosomal rearrangements.

Our findings indicated that Egcen centromeric tandemly repeated sequences of *Ensete* are absent in *Musa* [19]. Because these repeats were also detected in *Musella*, which is the third genus in the family Musaceae, one possible reason for the absence of Egcen repeats in *Musa* might be the loss of segments stemming from ancestral centromere breakage; this is plausible given that centromeres are hotspots of chromosome rearrangements [172, 173].
<table><caption>**Table 2:** Possible biosynthetic gene clusters identified in Musaceae species. Tomatine 1 and 2 are tomatine clusters that locate in different chromosomes when previously identified.</caption><thead><tr><th>Cluster</th><th>M. beccarii</th><th>E. glaucum</th><th>M. balbisiana</th><th>M. acuminata</th><th>M. schizocarpa</th></tr></thead><tbody><tr><td>Saccharide</td><td>2</td><td>3</td><td>3</td><td>2</td><td>3</td></tr><tr><td>Solanum l tomatine 1</td><td>11</td><td>7</td><td>10</td><td>13</td><td>9</td></tr><tr><td>Solanum l tomatine 1-Saccharide</td><td>0</td><td>1</td><td>0</td><td>1</td><td>1</td></tr><tr><td>Solanum l tomatine 1-Tomatine 2</td><td>2</td><td>1</td><td>2</td><td>2</td><td>3</td></tr><tr><td>T3PKS</td><td>21</td><td>17</td><td>18</td><td>22</td><td>14</td></tr><tr><td>T3PKS-Saccharide</td><td>2</td><td>1</td><td>0</td><td>1</td><td>1</td></tr><tr><td>T3PKS-Solanum l tomatine 1</td><td>6</td><td>8</td><td>9</td><td>9</td><td>11</td></tr><tr><td>T3PKS-Solanum l tomatine 1-Tomatine 2</td><td>3</td><td>1</td><td>0</td><td>2</td><td>3</td></tr><tr><td>T3PKS-Terpene</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td></tr><tr><td>T3PKS-Terpene-Solanum l tomatine 1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td></tr><tr><td>T3PKS-Tomatine 2</td><td>2</td><td>2</td><td>3</td><td>1</td><td>0</td></tr><tr><td>Terpene</td><td>3</td><td>6</td><td>8</td><td>8</td><td>7</td></tr><tr><td>Terpene-Solanum l tomatine 1</td><td>3</td><td>5</td><td>2</td><td>2</td><td>2</td></tr><tr><td>Terpene-Solanum l tomatine 1-Tomatine 2</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td></tr><tr><td>Tomatine 2</td><td>8</td><td>4</td><td>6</td><td>6</td><td>6</td></tr><tr><td>Other</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0</td></tr><tr><td>Other-Solanum l tomatine 1-Tomatine 2</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td></tr><tr><td>Total</td><td>66</td><td>59</td><td>64</td><td>72</td><td>62</td></tr></tbody></table>

## AS related to DNA repair systems

We detected AS in at least 11.7% of all the genes in *M. beccarii*. However, this is underestimated because we did not examine the transcriptomes of different tissues and different developmental stages. The enrichment analysis of these genes with AS revealed that they were involved in important cellular responses and DNA repair systems, including DNA repair, nucleotide excision repair, and replication and repair (Supplementary Tables S20 and S21). A wide variety of stress conditions can induce DNA damage. DNA repair systems are therefore important for maintaining the stability of chromosomes in eukaryotic cells [174–176]. AS is a post-transcriptional mechanism that produces many functional proteins from a limited number of genes. AS in *M. beccarii* plays a key role in DNA repair pathways and other processes that mediate evolutionary adaptation.

## Conclusion

The assembly of a genome from a member of the section *Calimusa* in the genus *Musa* is important for the development of a pangenome model of Musaceae. The new data reveal extensive rearrangements and expansions, and they provide new insights into the range of structural chromosome variation within the family Musaceae. The genes and TFs identified and our structural analysis of the genome are important for conserving biodiversity within the genus *Musa*. Our findings also have implications for breeding novel variants and addressing some of the major challenges faced by banana crop production.

## Additional Files

**Supplementary Table S1.** Libraries for genome assembly and annotation for *Musa beccarii*.

**Supplementary Table S2.** Protein sequences of 3 species used for gene prediction.

**Supplementary Table S3.** Species used for comparative genomics.

**Supplementary Table S4.** Species pairs and their estimated divergence times used for time calibration points to infer the time-calibrated phylogeny of *Musa beccarii*.

**Supplementary Table S5.** Genome size (bp) estimation using different programs.

**Supplementary Table S6.** Assembly results using different assemblers.

**Supplementary Table S7.** Repeat content of assemblies of Musaceae.

**Supplementary Table S8.** Summary of gene functional annotations of the *Musa beccarii* genome using different databases.

**Supplementary Table S9.** Summary of genes in Musaceae and the other compared species.

**Supplementary Table S10.** Summary of MYB transcription factors identified by MYB_annotator and their comparison with Or thermoFinder results.

**Supplementary Table S11.** MYB genes identified by MYB_annotator in *Musa beccarii*.

**Supplementary Table S12.** Summer of MYB gene functions.

**Supplementary Table S13.** Statistics of gene families in different species.

**Supplementary Table S14.** GO enrichment results for Musaceae-specific gene families in *Musa beccarii*.

**Supplementary Table S15.** KEGG enrichment results for Musaceae-specific gene families in *Musa beccarii*.

**Supplementary Table S16.** GO enrichment results for significantly expanded gene families.

**Supplementary Table S17.** KEGG enrichment results for significantly expanded gene families.

**Supplementary Table S18.** GO enrichment results for significantly contracted gene families.

**Supplementary Table S19.** KEGG enrichment results for significantly contracted gene families.

**Supplementary Table S20.** GO enrichment results for alternative splicing genes.

**Supplementary Table S21.** KEGG enrichment results for alternative splicing genes.

**Supplementary Table S22.** GO enrichment results for WGD-type genes.

**Supplementary Table S23.** KEGG enrichment results for WGD-type genes.

**Supplementary Table S24.** GO enrichment results for TD-type genes.
**Supplementary Table S25.** KEGG enrichment results for TD-type genes.

**Supplementary Table S26.** GO enrichment results for PD-type genes.

**Supplementary Table S27.** KEGG enrichment results for PD-type genes.

**Supplementary Table S28.** GO enrichment results for TRD-type genes.

**Supplementary Table S29.** GO enrichment results for DSD-type genes.

**Supplementary Table S30.** KEGG enrichment results for DSD-type genes.

**Supplementary Table S31.** Results of the syntenic block analysis.

**Supplementary Table S32.** Biosynthetic gene clusters in Musaceae.

**Supplementary Table S33.** NLR-Annotator results in *Musa beccarii*.

**Supplementary Table S34.** NLR-Annotator results in Musaceae using their predicted genes and comparison of the results using the genomes.

**Supplementary Table S35.** Structure variations (SVs) detected by SVanalyzer between *Musa beccarii* and *Musa acuminata*.

**Supplementary Table S36.** Summary of ABC transporter-related genes.

**Supplementary Table S37.** Summary of glycosyltransferase genes.

**Supplementary Fig. S1.** A scheme showing *Musa beccarii* genome assembly.

**Supplementary Fig. S2.** (A) Nanica repetitive sequences and (B) NBS-LRR genes identified by NLR-Annotator across the *Musa beccarii* genome.

**Supplementary Fig. S3.** Treemap showing hierarchy for GO terms enriched for Musaceae-specific gene families of *Musa beccarii* in biological process.

**Supplementary Fig. S4.** Upset plot showing the intersection of the gene family in *Musa*. *Musa* species are presented in the rows, and the bar beside each species shows its total gene number. Black circles and vertical lines between the rows represent the intersection of gene families between species. Barplot indicates the total gene family count in each intersection.

**Supplementary Fig. S5.** Summary of GO annotations of Musaceae shared genes.

**Supplementary Fig. S6.** Treemap showing hierarchy for GO terms enriched with significantly expanded gene families in biological process.

**Supplementary Fig. S7.** Treemap showing hierarchy for GO terms enriched with alternative splicing genes in biological process.

**Supplementary Fig. S8.** Treemap showing hierarchy for GO terms enriched with genes related to WGD in biological process.

**Supplementary Fig. S9.** KEGG enrichment for genes related to (A) whole-genome duplication, (B) tandem duplications, (C) proximal duplications, and (D) dispersed duplications tested by DupGen_Finder.

**Supplementary Fig. S10.** Treemap showing the hierarchy for GO terms enriched with genes related to tandem duplications in biological process.

**Supplementary Fig. S11.** Treemap showing the hierarchy for GO terms enriched with genes related to dispersed duplications in biological process.

**Supplementary Fig. S12.** Synteny plot (Synvisio) between *Ensete glaucum* (egxx) and *Musa acuminata* (mpxx), *Musa beccarii* (bexx) genomes. Syntenic blocks of high homology are indicated by uniformly colored areas in the graphs.

**Supplementary Fig. S13.** Coverage distributions of Illumina WGS reads mapping against *Musa beccarii* assembly.

**Supplementary Fig. S14.** Histograms showing the allele frequencies of SNVs (single-nucleotide variants) obtained from the Illumina WGS reads mapping against the *Musa beccarii* assembly.

**Supplementary Fig. S15.** KEGG pathway of flavonoid biosynthesis. In the pathway, the genes identified in *Musa beccarii* are boxed in green or pink colors.

## Abbreviations

ABC: adenosine triphosphate-binding cassette; AS: alternative splicing; ATP: adenosine triphosphate; BGC: biosynthetic gene cluster; BH: Benjamini-Hochberg; BLAST: Basic Local Alignment Search Tool; bp: base pair; BUSCO: Benchmarking Universal Single-Copy Orthologs; cDNA: complementary DNA; CNL: N-terminal coiled-coil motif NBS-LRR; DSD: dispersed duplication; Gbp: gigabase pair; GT: glycosyltransferase; KEGG: Kyoto Encyclopedia of Genes and Genomes; LCA: last common ancestor; LTR: long terminal repeat; Mbp: megabase pair; NBS-LRR: nucleotide-binding site-leucine-rich repeat (NBS-LRR); NOMT: naringenin 7-O-methyltransferase; PD: proximal duplication; rDNA: ribosomal DNA; RNL: resistance to powdery mildew 8 NBS-LRR; SA: steroidal alkaloid; TD: tandem duplication; TF: transcription factor; TIR: terminal inverted repeat; TNL: Toll/interleukin-1 receptor NBS-LRR; TRD: transposed duplication; WGD: whole-genome duplication; WGS: whole-genome sequencing.

## Data Availability

We deposited the sequenced reads to NCBI Sequence Read Archive under the accession number SRR16526886 for the Nanopore reads, SRR16526885 for PacBio HiFi reads, SRR16526887 for the Illumina WGS reads, SRR16588090 and SRR16588091 for the Illumina Hi-C reads, SRR16351760 for the Illumina RNA-seq reads, and SRR16351759 for the PacBio Iso-seq reads. The high-quality genome sequence was submitted to GenBank under the accession number JAIWVJ000000000. Genome Assembly, gene annotation data, and transcriptomic data are also available on the Banaana Genome Hub [177] for download or exploration via a dedicated genome browser (Jbrowse) and syntenic browser (SynVisio). All supporting data are available in Figshare [38] and the GigaScience GigaDB database [178].

## Ethics Approval and Consent to Participate

No ethical approval/permission is required to obtain the materials and perform the research in this study. Plant material was collected before 1993 and was received at the South China Botanical Garden under an agreement from a European botanical garden.

## Competing Interests

The authors declare that they have no competing interests.

## Funding

This work was financially supported by the National Natural Science Foundation of China (No. 32070237, 31261140366), and the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDB3100000).
## Authors' Contributions

X.-J.G. and Z.-F.W. designed this experiment. Z.-F.W., M.R., G.D., and P.H.H. conducted genetic work and data analyses. Z.-F.W., X.-J.G., M.R., and P.H.H. drafted the manuscript. Z.-F.W., M.R., G.D., P.H.H., and X.-J.G. revised the manuscript. All authors gave final approval of the paper.

## References

1. Li, L-F, Häkkinen, M, Yuan, Y-M, et al. Molecular phylogeny and systematics of the banana family (Musaceae) inferred from multiple nuclear and chloroplast DNA fragments, with a special reference to the genus *Musa*. *Mol Phylogenet Evol* 2010;**57**(1):1–10.

2. Häkkinen, M. Reappraisal of sectional taxonomy in *Musa* (Musaceae). *Taxon* 2013;**62**:809–13.

3. Fu, N, Ji, M, Rouard, M, et al. Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships. *BMC Genomics* 2022;**23**(1):223.

4. Christelová, P, Valárik, M, Hřibová, E, et al. A multi gene sequence-based phylogeny of the Musaceae (banana) family. *BMC Evol Biol* 2011;**11**:103.

5. Häkkinen, M, Suchuánková, P, Doleželová, M, et al. Karyological observation in *Musa beccarii* var. *Hottana* (Musaceae). *Acta Phytotax Geobot* 2007;**58**(2/3):112–8.

6. Häkkinen, M, Teo, CH, Othman, YR. Genome constitution for *Musa beccarii* (Musaceae) varieties. *Acta Phytotaxon Sin* 2007;**45**(1):69–74.

7. Natarajan, N, Sundararajan, S, Ramalingam, S, et al. Efficient and rapid in-vitro plantlet regeneration via somatic embryogenesis in ornamental bananas (*Musa* spp.). *Biologia (Bratisl)* 2020;**75**:317–26.

8. Rashid, K, Nezhadahmadi, A, Othman, RY, et al. Micropropagation of ornamental plant *Musa beccarii* through tissue culture technique using suckers and male buds as explants. *Life Sci J* 2012;**9**(4):2046–53.

9. Allen, R. *Musa beccarii*. The IUCN Red List of Threatened Species 2019: e.T121033043A121033225. 2019. https://www.iucnredlist.org/species/121033043/121033225. Accessed 13 July 2022.

10. Droc, G, Martin, G, Guignon, V, et al. The banana genome hub: a community database for genomics in the Musaceae. *Hortic Res* 2022;**9**:uhac221.

11. D'Hont, A, Denoeud, F, Aury, JM, et al. The banana (*Musa acuminata*) genome and the evolution of monocotyledonous plants. *Nature* 2012;**488**(7410):213–7.

12. Belser, C, Baurens, FC, Noel, B, et al. Telomere-to-telomere gapless chromosomes of banana using nanopore sequencing. *Commun Biol* 2021;**4**(1):1047.

13. Davey, MW, Gudimella, R, Harikrishna, JA, et al. A draft *Musa balbisiana* genome sequence for molecular genetics in polyploid, inter- and intra-specific *Musa* hybrids. *BMC Genomics* 2013; **14**:683.

14. Wang, Z, Miao, H, Liu, J, et al. *Musa balbisiana* genome reveals subgenome evolution and functional divergence. *Nat Plants* 2019;**5**:810–21.

15. Wu, W, Yang, YL, He, WM, et al. Whole genome sequencing of a banana wild relative *Musa itinerans* provides insights into lineage-specific diversification of the *Musa* genus. *Sci Rep* 2016;**6**:31586.

16. Belser, C, Istace, B, Denis, E, et al. Chromosome-scale assemblies of plant genomes using Nanopore long reads and optical maps. *Nat Plants* 2018;**4**(11):879–87.

17. Eyland, D, Breton, C, Sardos, J, et al. Filling the gaps in gene banks: collecting, characterizing, and phenotyping wild banana relatives of Papua New Guinea. *Crop Sci* 2021;**61**: 137–49.

18. Galvez, LC, Koh, RBL, Barbosa, CFC, et al. Sequencing and de novo assembly of Abaca (*Musa textilis* Nee) var. Abuab genome. *Genes* 2021;**12**(8):1202.

19. Wang, Z, Rouard, M, Biswas, MK, et al. A chromosome-level reference genome of *Ensete glaucum* gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae. *Gigascience* 2022;**11**:giac027.

20. Belton, JM, McCord, RP, Gibcus, JH, et al. Hi-C: a comprehensive technique to capture the conformation of genomes. *Methods* 2012;**58**(3):268–76.

21. Joshi, NA, Fass, JN. Sickle: a sliding-window, adaptive, quality-based trimming tool for FastQ files (version 1.33). 2011. https://github.com/najoshi/sickle. Accessed 3 September 2021.

22. Długosz, M, Deorowicz, S. RECKONER: read error corrector based on KMC. *Bioinformatics* 2017;**33**:1086–9.

23. CCS (Version 6.0.0). https://github.com/PacificBiosciences/ccs/releases/tag/v6.0.0. Accessed 12 February 2021.

24. Chikhi, R, Medvedev, P. Informed and automated k-mer size selection for genome assembly. *Bioinformatics* 2014;**30**:31–7.

25. Vulture, GW, Sedlazeck, FJ, Nattestad, M, et al. GenomeScope: fast reference-free genome profiling from short reads. *Bioinformatics* 2017;**33**:2202–4.

26. Sun, H, Ding, J, Piednoël, M, et al. findGSE: estimating genome size variation within human and Arabidopsis using k-mer frequencies. *Bioinformatics* 2018;**34**:550–7.

27. Liu, B, Shi, Y, Yuan, J, et al. Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome project. *arXiv* 2013:1308.2012. https://arxiv.org/abs/1308.2012. Accessed 8 November 2022.

28. Pucker, B. Mapping-based genome size estimation. *Biorxiv* 2019. https://doi.org/10.1101/607390. Accessed 8 November 2022.

29. Gilbert, DG. Genes ruler for genomes, gnodes, measures assembly accuracy in animals and plants. *Biorxiv* 2022. https://doi.org/10.1101/2022.05.13.491861.

30. IsoSeq (Version 3.0). https://github.com/ben-lerch/IsoSeq-3.0. Accessed 12 October 2019.

31. Porechop (Version 0.2.4). https://github.com/rrwick/Porechop/rleleases/tag/v0.2.4. Accessed 8 January 2019.

32. Sim, SB, Corpuz, RL, Simmonds, TJ, et al. HiFiAdapterFilt, a memory efficient read processing pipeline, prevents occurrence of adapter sequence in PacBio HiFi reads and their negative impacts on genome assembly. *BMC Genom* 2022;**23**(1): 157.

33. NextDenovo (Version 2.3.1). https://github.com/Nextomics/NextDenovo/releases/tag/v2.3.1. Accessed 22 September 2020.

34. Kolmogorov, M, Yuan, J, Lin, Y, et al. Assembly of long, error-prone reads using repeat graphs. *Nat Biotechnol* 2019;**37**:540–6.

35. Koren, S, Walenz, BP, Berlin, K, et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. *Genome Res* 2017;**27**(5):722–36.

36. Cheng, H, Concepcion, GT, Feng, X, et al. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. *Nat Methods* 2021;**18**(2):170–5.

37. Nurk, S, Walenz, BP, Rhiea, A, et al. HiCanu: accurate assembly of segmental duplications, satellites, and allelic variants from high-fidelity long reads. *Genome Res* 2020;**30**(9):1291–305.

38. Wang, Z-F. Genome assembly of *Musa beccarii*. Figshare. Dataset. 2022. https://doi.org/10.6084/m9.figshare.19165280.v13.
39. Vaser, R, Sović, I, Nagarajan, N, et al. Fast and accurate de novo genome assembly from long uncorrected reads. *Genome Res* 2017;**27**(5):737-46.

40. Aury, JM, Istace, B. Hapo-G, haplotype-aware polishing of genome assemblies with accurate reads. *NAR Genom Bioinform* 2021;**3**(2):lqab034.

41. Pseudohaploid. https://github.com/schatzlab/pseudohaploid. Accessed 28 August 2020.

42. Guan, DF, McCarthy, SA, Wood, J, et al. Identifying and removing haplotypic duplication in primary genome assemblies. *Bioinformatics* 2020;**36**:2896-8.

43. Chen, Y, Zhang, Y, Wang, AY, et al. Accurate long-read de novo assembly evaluation with Inspector. *Genome Biol* 2021;**22**:312.

44. Alonge, M, Lebeigle, L, Kirsche, M, et al. Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing. *Genome Biol* 2022;**23**:258.

45. Scaffhic (Version 1.1). https://github.com/wtsi-hpag/scaffHiC. Accessed 28 August 2020.

46. Durand, NC, Shamim, MS, Machol, I, et al. Juicer provides a one-click system for analyzing loop-resolution hi-C experiments. *Cell Syst* 2016;**3**(1):95-8.

47. Dudchenko, O, Batra, SS, Omer, AD, et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. *Science* 2017;**356**(6333):92-5.

48. Xu, M, Guo, L, Gu, S, et al. TGS-GapCloser: a fast and accurate gap closer for large genomes with low coverage of error-prone long reads. *Gigascience* 2020;**9**(9):giaa094.

49. Seppey, M, Manni, M, Zdobnov, EM. BUSCO: assessing genome assembly and annotation completeness. *Methods Mol Biol* 2019;**1962**:227-45.

50. Li, H, Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler Transform. *Bioinformatics* 2009;**25**:1754-60.

51. Danecek, P, Bonfield, JK, Liddle, J, et al. Twelve years of SAMtools and BCFtools. *Gigascience* 2021;**10**(2):giab008.

52. Yang, LA, Chang, YJ, Chen, SH, et al. SQUAT: a sequencing quality assessment tool for data quality assessments of genome assemblies, BMC Genom* 2019;**19**:238.

53. Ou, S, Su, W, Liao, Y, et al. Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline. *Genome Biol* 2019;**20**:275.

54. Girgis, HZ. Red: an intelligent, rapid, accurate tool for detecting repeats de-novo on the genomic scale. BMC Bioinf 2015;**16**(1):227.

55. Quinlan, AR, Hall, IM. BEDTools: a flexible suite of utilities for comparing genomic features. *Bioinformatics* 2010;**26**(6):841-2.

56. Camacho, C, Coulouris, G, Avagyan, V, et al. BLAST+: architecture and applications. BMC Bioinf* 2009;**10**:421.

57. Nanica sequences. *Banana Genome Hub*. 2022. https://bana.genome-hub.southgreen.fr/node/50/353. Accessed 15 April 2022.

58. Cook, DE, Valle-Inclan, JE, Pajoro, A, et al. Long-read annotation: automated eukaryotic genome annotation based on long-read cDNA sequencing. *Plant Physiol* 2019;**179**(1):38-54.

59. Funannotate (Version 1.8.7). https://github.com/nextgenusfs/funannotate/releases/tag/v1.8.7. Accessed 28 April 2021.

60. Zhang, H, Tanner, Y, Huang, L, et al. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. *Nucleic Acids Res* 2018;**46**:W95-W101.

61. Huerta-Cepas, J, Forslund, K, Coelho, LP, et al. Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. *Mol Biol Evol* 2017;**34**:2115-22.

62. The Gene Ontology Consortium. The gene ontology resource: 20 years and still going strong. *Nucleic Acids Res* 2019;**47**(D1):D330-8.

63. Ashburner, M, Ball, CA, Blake, JA, et al. Gene ontology: tool for the unification of biology. *Nat Genet* 2000;**25**:25-9.

64. Kanehisa, M, Soto, Y, Kawashima, M, et al. KEGG as a reference resource for gene and protein annotation. *Nucleic Acids Res* 2016;**44**(D1):D457-62.

65. Mitchell, AL, Attwood, TK, Babbitt, PC, et al. InterPro in 2019: improving coverage, classification and access to protein sequence annotations. *Nucleic Acids Res* 2019;**47**(D1):D351-60.

66. Rawlings, ND, Barrett, AJ, Thomas, PD, et al. The merops database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. *Nucleic Acids Res* 2018;**46**(D1):D624-32.

67. El-Gebali, S, Mistry, J, Bateman, A, et al. The Pfam protein families database in 2019. *Nucleic Acids Res* 2019;**47**(D1):D427-32.

68. The UniProt Consortium. UniProt: a worldwide hub of protein knowledge. *Nucleic Acids Res* 2019;**47**(D1):D506-15.

69. Trincado, JL, Entizne, JC, Hysenaj, G, et al. SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions. *Genome Biol* 2018;**19**:40.

70. Zheng, Y, Jiao, C, Sun, H, et al. iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. *Mol Plant* 2016;**9**:1667-70.

71. Pucker, B. Automatic identification and annotation of MYB gene family members in plants. *BMC Genom* 2022;**23**:220

72. Emms, DM, Kelly, S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. *Genome Biol* 2015;**16**:157.

73. Emms, DM, Kelly, S. OrthoFinder: phylogenetic orthology inference for comparative genomics. *Genome Biol* 2019;**20**:238.

74. Kozlov, AM, Darriba, D, Flouri, T, et al. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. *Bioinformatics* 2019;**35**(21):4453-5.

75. Darriba, D, Posada, D, Kozlov, AM, et al. ModelTest-NG: a new and scalable tool for the selection of DNA and protein evolutionary models. *Mol Biol Evol* 2020;**37**(1):291-4.

76. dos Reis, M, Zhu, T, Yang, Z. The impact of the rate prior on Bayesian estimation of divergence times with multiple loci. *Syst Biol* 2014;**63**:555-65

77. Han, MV, Thomas, GWC, Jose, LM, et al. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using cafe 3. *Mol Biol Evol* 2013;**30**(8):1987-97.

78. Chen, CJ, Chen, H, Zhang, Y, et al. TBtools- an integrative toolkit developed for interactive analyses of big biological data. *Mol Plant* 2020;**13**(8):1194-202.

79. Supek, F, Bošnjak, M, Škunca, N, et al. REVIGO summarizes and visualizes long lists of gene ontology terms. *PLoS One* 2011;**6**(7):e21800.

80. Zwaenepoel, A, de Peer, YV. WGD-simple command line tools for the analysis of ancient whole-genome duplications. *Bioinformatics* 2019;**35**:2153-5.

81. Sensalari, C, Maere, S, Lohaus, R. Ksrates: positioning whole-genome duplications relative to speciation events in *K. S* distributions. *Bioinformatics* 2022;**38**(2):530-2.

82. Qiao, X, Li, QH, Yin, H, et al. Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants. *Genome Biol* 2019;**20**:38.

83. JCVI (Version 1.1.19). https://github.com/tanghaibao/jcvi/releases/tag/v1.1.19. Accessed 8 November 2021.
84. Yu, Y, Ouyang, Y, Yao, W. shinyCircos: an R/Shiny application for interactive creation of Circos plot. *Bioinformatics* 2018; **34**(7):1229–31.

85. Wang, Y, Tang, H, Debarry, JD, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. *Nucleic Acids Res* 2012; **40**(7):e49

86. Bandi, V, Gutwin, C. Interactive exploration of genomic conservation. In *Proceedings of the 46th Graphics Interface Conference on Proceedings of Graphics Interface 2020 (GI'20)*. Waterloo: Canadian Human-Computer Communications Society; 2020.

87. Cabanettes, F, Köpp, C. D-GENIES: dot plot large genomes in an interactive, efficient and simple way. *PeerJ* 2018; **6**:e4958.

88. Kautsar, SA, Duran, HGS, Blin, K, et al. plantiSMASH: automated identification, annotation and expression analysis of plant biosynthetic gene clusters. *Nucleic Acids Res* 2017; **45**(W1):W55–63.

89. Töpfer, N, Fuchs, LM, Aharoni, A. The PhytoClust tool for metabolic gene clusters discovery in plant genomes. *Nucleic Acids Res* 2017; **45**(12):7049–63

90. McHale, L, Tan, X, Koehl, P, et al. Plant NBS-LRR proteins: adaptable guards. *Genome Biol* 2006; **7**(4):212.

91. Shao, Z-Q, Xue, J-Y, Wu, P, et al. Large-scale analyses of angiosperm nucleotide-binding site-leucine-rich repeat genes reveal three anciently diverged classes with distinct evolutionary patterns. *Plant Physiol* 2016; **170**(4):2095–109.

92. Guo, X, Fang, D, Sahu, SK, et al. Chloranthus genome provides insights into the early diversification of angiosperms. *Nat Commun* 2021; **12**(1):6930.

93. Steuernagel, B, Witek, K, Krattinger, SG, et al. The NLR-annotator tool enables annotation of the intracellular immune receptor repertoire. *Plant Physiol* 2020; **183**:468–82.

94. Vakirlis, N, Sarilar, V, Drillon, G, et al. Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus. *Genome Res* 2016; **26**(7):918–32.

95. Drillon, G, Carbone, A, Fischer, G. SynChro: a fast and easy tool to reconstruct and visualize synteny blocks along eukaryotic chromosomes. *PLoS One* 2014; **9**(3):e92621.

96. Pucker, B, Pandey, A, Weisshaar, B, et al. The R2R3-MYB gene family in banana (*Musa acuminata*): genome-wide identification, classification and expression patterns. *PLoS One* 2020; **15**(10):e0239275.

97. Busche, M, Pucker, B, Weisshaar, B, et al. Three R2R3-MYB transcription factors from banana (*Musa spp.*) activate structural anthocyanin biosynthesis genes as part of an MBW complex. *Biorxiv* 2022. https://doi.org/10.1101/2022.08.15.503939. Accessed 25 October 2022.

98. Novák, P, Hřibová, E, Neumann, P, et al. Genome-wide analysis of repeat diversity across the family Musaceae. *PLoS One* 2014; **9**(6):e98918.

99. Bartoš, J, Alkhimova, O, Doleželová, M, et al. Nuclear genome size and genomic distribution of ribosomal DNA in *Musa* and *Ensete* (Musaceae): taxonomic implications. *Cytogenet Genome Res* 2005; **109**(1–3):50–57.

100. Lysák, MA, Doleželová, M, Horry, JP, et al. Flow cytometric analysis of nuclear DNA content in *Musa*. *Theor Appl Genet* 1999; **98**:1344–50.

101. Doležel, J, Bartoš, J. Plant DNA flow cytometry and estimation of nuclear genome size. *Ann Bot* 2005; **95**(1):99–110.

102. Pellicer, J, Powell, RF, Leitch, IJ. The application of flow cytometry for estimating genome size, ploidy level endopolyploidy, and reproductive modes in plants. *Methods Mol Biol* 2021; **2222**:325–61.

103. Van't Ho, J, Sparrow, AH. A relationship between DNA content, nuclear volume, and minimum mitotic cycle time. *Proc Natl Acad Sci USA* 1963; **49**:897–902.

104. Šmarda, P, Horová, L, Bureš, P, et al. Stabilizing selection on genome size in a population of *Festuca pallens* under conditions of intensive intraspecific competition. *New Phytol* 2010; **187**(4):1195–204.

105. Piegu, B, Guyot, R, Picault, N, et al. Doubling genome size without polyploidization: dynamics of retrotransposition-driven genomic expansions in *Oryza australiensis*, a wild relative of rice. *Genome Res* 2006; **16**(10):1262–9.

106. Macas, J, Novák, P, Pellicer, J, et al. In depth characterization of repetitive DNA in 23 plant genomes reveals sources of genome size variation in the Legume Tribe Fabeae. *PLoS One* 2015; **10**(11):e0143424.

107. Wicker, T, Gundlach, H, Spannagl, M, et al. Impact of transposable elements on genome structure and evolution in bread wheat. *Genome Biol* 2018; **19**(1):103.

108. Novák, P, Guignard, MS, Neumann, P, et al. Repeat-sequence turnover shifts fundamentally in species with large genomes. *Nat Plants* 2020; **6**(11):1325–9.

109. Busche, M, Pucker, B, Viehöver, P, et al. Genome sequencing of *Musa acuminata* Dwarf Cavendish reveals a duplication of a large segment of chromosome 2. *G3 (Bethesda)* 2020; **10**(1):37–42.

110. Chen, F, Tholl, D, Bohlmann, J, et al. The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. *Plant J* 2011; **66**(1):212–9.

111. Pichersky, E, Raguso, RA. Why do plants produce so many terpenoid compounds? *New Phytol* 2016; **220**(3):655–8.

112. Jiang, SY, Jin, JJ, Sarojam, R, et al. A comprehensive survey on the terpene synthase gene family provides new insight into its evolutionary patterns. *Genome Biol Evol* 2019; **11**(8):2078–98.

113. Rojas, CM, Senthil-Kumar, M, Tzin, V, et al. Regulation of primary plant metabolism during plant-pathogen interactions and its contribution to plant defense. *Front Plant Sci* 2014; **5**:17.

114. Meshi, T, Iwabuchi, M. Plant transcription factors. *Plant Cell Physiol* 1995; **36**(8):1405–20.

115. Amorim, LLB, da Fonseca, DSR, Neto, JPB, et al. Transcription factors involved in plant resistance to pathogens. *Curr Protein Pept Sci* 2017; **18**(4):335–51.

116. Gani, U, Vishwakarma, RA, Misra, P. Membrane transporters: the key drivers of transport of secondary metabolites in plants. *Plant Cell Rep* 2021; **40**(1):1–18.

117. Pollard, M, Beisson, F, Li, Y, et al. Building lipid barriers: biosynthesis of cutin and suberin. *Trends Plant Sci* 2008; **13**(5):236–46.

118. Wang, A, Zha, Z, Yin, D, et al. Comparative transcriptome analysis of *Tilletia horrida* infection in resistant and susceptible rice (*Oryza sativa* L.) male sterile lines reveals potential candidate genes and resistance mechanisms. *Genomics* 2020; **112**(6):5214–26.

119. Ziv, C, Zhao, Z, Gao, YG, et al. Multifunctional roles of plant cuticle during plant-pathogen interactions. *Front Plant Sci* 2018; **9**:1088.

120. Baales, J, Zeisler-Diehl, VV, Schreiber, L. Analysis of extracellular cell wall lipids: wax, cutin, and suberin in leaves, roots, fruits, and seeds. *Methods Mol Biol* 2021; **2295**:275–93.

121. KEGG PATHWAY Database. https://www.genome.jp/kegg/patway.html. Accessed 27 January 2022.

122. Pighin, JA, Zheng, H, Balakshin, LJ, et al. Plant cuticular lipid export requires an ABC transporter. *Science* 2004; **306**(5696):702–4.
123. Elejalde-Palmett, C, Segundo, IMS, Garroum, I, et al. ABCG transporters export cutin precursors for the formation of the plant cuticle. *Curr Biol* 2021;31(10):2111-23.e9.

124. Kang, J, Park, J, Choi, H, et al. Plant ABC transporters. *Arabidopsis Book* 2011;9:e0153.

125. Bailly, A. Structure-function of plant ABC-transporters. In: M Geisler, editor. *Plant ABC Transporters*. Cham, Switzerland: Springer; 2014; p. 219-40.

126. Do, THT, Martinoia, E, Lee, Y. Functions of ABC transporters in plant growth and development. *Curr Opin Plant Biol* 2018;41:32-8.

127. Banasiak, J, Jasiński, M. ATP-binding cassette transporters in nonmodel plants. *New Phytol* 2022;233:1597-612.

128. Winkel-Shirley, B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. *Plant Physiol* 2001;126(2):485-93.

129. Pucker, B, Reiher, F, Schilbert, HM. Automatic identification of players in the flavonoid biosynthesis with application on the biomedicinal plant *Croton tiglum*. *Plants* (Basel) 2020;9(9):1103.

130. Liu, W, Feng, Y, Yu, S, et al. The flavonoid biosynthesis network in plants. *Int J Mol Sci* 2021;22:12824.

131. Shen, N, Wang, T, Gan, Q, et al. Plant flavonoids: classification, distribution, biosynthesis, and antioxidant activity. *Food Chem* 2022;383:132531.

132. Pandey, A, Alok, A, Lakhwani, D, et al. Genome-wide expression analysis and metabolite profiling elucidate transcriptional regulation of flavonoid biosynthesis and modulation under abiotic stresses in Banana. *Sci Rep* 2016;6:31361.

133. Sun, X, Gao, P, Zhang, J, et al. Characteristics of flavonoids biosynthesis and the differential expression analysis of the key enzyme genes in *Musa AAA* Group cv Brazilian fruit pulp. *Mol Plant Breeding* 2018;16(7):2116-23.

134. Murata, K, Kitano, T, Yoshimoto, R, et al. Natural variation in the expression and catalytic activity of a naringenin 7-O-methyltransferase influences antifungal defenses in diverse rice cultivars. *Plant J* 2019;101(5):1103-17.

135. Busche, M, Acatay, C, Martens, S, et al. Functional characterisation of banana (*Musa* spp.) 2-oxoglutarate-dependent dioxygenases involved in flavonoid biosynthesis. *Front Plant Sci* 2021;12:701780.

136. Petrussa, E, Braidot, E, Zancani, M, et al. Plant flavonoids—biosynthesis, transport and involvement in stress responses. *Int J Mol Sci* 2013;14:14950-73.

137. Pucker, B, Selmar, D. Biochemistry and molecular basis of intracellular flavonoid transport in plants. *Plants* 2022;11:963.

138. Pazmiño-Durán, EA, Giusti, MM, Wrolstad, RE, et al. Anthocyanins from banana bracts (*Musa X paradisiaca*) as potential food colorants. *Food Chem* 2001;73(3):327-32.

139. Fu, X, Cheng, S, Liao, Y, et al. Comparative analysis of pigments in red and yellow banana fruit. *Food Chem* 2018;239:1009-18.

140. Deng, S, Cheng, C, Liu, Z, et al. Comparative transcriptome analysis reveals a role for anthocyanin biosynthesis genes in the formation of purple peel in Minhou wild banana (*Musa itinerans* Cheesman). *J Hortic Sci Biotechnol* 2019;94(2):184-200.

141. Jiao, F, Zhao, L, Wu, X, et al. Metabolome and transcriptome analyses of the molecular mechanisms of flower color mutation in tobacco. *BMC Genom* 2020;21:611.

142. Zhang, X, Lin, S, Peng, D, et al. Integrated multi-omic data and analyses reveal the pathways underlying key ornamental traits in carnation flowers. *Plant Biotechnol J* 2022;20(6):1182-96.

143. Nützmann, H-W, Huang, A, Osbourn, A. Plant metabolic clusters—from genetics to genomics. *New Phytol* 2016;211(3):771-89.

144. Polturak, G, Osbourn, A. The emerging role of biosynthetic gene clusters in plant defense and plant interactions. *PLoS Pathog* 2021;17(7):e1009698.

145. Yu, D, Xu, F, Zeng, J, et al. Type III polyketide synthases in natural product biosynthesis. *IUBMB Life* 2012;64(4):285-95.

146. Flores-Sanchez, IJ, Verpoorte, R. Plant polyketide synthases: a fascinating group of enzymes. *Plant Physiol Biochem* 2009;47(3):167-74.

147. Mhlanga, M. Plant polyketides. *Nat Biotechnol* 1999;17:9.

148. Rajesh, T, Tiwari, MK, Thiagarajan, S, et al. Type III polyketide synthases: current state and perspectives. In: P Arora, editor. *Microbial Technology for the Welfare of Society*. Singapore: Springer; 2019; p. 183-200.

149. Pothiraj, R, Ravikumar, MJ, Suthanthiram, B, et al. Genome-scale analyses of polyketide synthases in banana: phylogenetics and expression profiling forecast their candidacy in specialized metabolism. *Gene* 2021;778:145472.

150. Gröger, D. Terpenoid and steroid alkaloids. In: F Constabel, IK Vasileditors. *Phytochemicals in Plant Cell Cultures*. New York: Academic Press; 1988; p. 435-48.

151. Piasecka, A, Jedrzejczak-Rey, N, Bednarek, P. Secondary metabolites in plant innate immunity: conserved function of divergent chemicals. *New Phytol* 2015;206(3):948-64.

152. Hoagland, RE. Toxicity of tomatine and tomatidine on weeds, crops and phytopathogens fungi. *Allelopathy J* 2009;23(2):425-35.

153. Nakayasu, M, Akiyama, R, Kobayashi, M, et al. Identification of α-tomatine 23-hydroxylase involved in the detoxification of a bitter glycoalkaloid. *Plant Cell Physiol* 2020;61(1):21-8.

154. Itkin, M, Rogachev, I, Rogachev, I, et al. GLYCOALKALOID METABOLISM1 is required for steroidal alkaloid glycosylation and prevention of phytotoxicity in tomato. *Plant Cell* 2011;23(12):4507-25.

155. Itkin, M, Heinig, U, Tzfadia, O, et al. Biosynthesis of antinutritional alkaloids in solanaceous crops is mediated by clustered genes. *Science* 2013;341(6142):175-9.

156. You, Y, van Kan, JAL. Bitter and sweet make tomato hard to (b)eat. *New Phytol* 2020;230(1):90-100.

157. Cruz-Cruz, CA, Ramírez-Tec, G, García-Sosa, K, et al. Phytoanticipins from banana (*Musa acuminata* cv. Grande Naine) plants, with antifungal activity against *mycosphaerella fijiensis*, the causal agent of black Sigatoka. *Eur J Plant Pathol* 2010;126(4):459-63.

158. Soares, JMS, Rocha, AJ, Nascimento, FS, et al. Genetic improvement for resistance to black Sigatoka in bananas: a systematic review. *Front Plant Sci* 2021;12:657916.

159. Gachon, CM, Langlois-Meurinne, M, Saindrenan, P. Plant secondary metabolism glycosyltransferases: the emerging functional analysis. *Trends Plant Sci* 2005;10(11):542-9.

160. Lairson, LL, Henrissat, B, Davies, GJ, et al. Glycosyltransferases: structures, functions, and mechanisms. *Annu Rev Biochem* 2008;77:521-55.

161. He, B, Bai, X, Tan, Y, et al. Glycosyltransferases: mining, engineering and applications in biosynthesis of glycosylated plant natural products. *Synth Syst Biotechnol* 2022;7:602-20.

162. Wang, J, Hou, B-K. Glycosyltransferases: key players involved in the modification of plant secondary metabolites. *Front Biol China* 2009;4(1):39-46.

163. Phan, JL, Tucker, MR, Khor, SF, et al. Differences in glycosyltransferase family 61 accompany variation in seed coat mucilage composition in *Plantago* spp. *J Exp Bot* 2016;67(22):6481-95.
164. Cenci, A, Chantret, N, Rouard, M. Glycosyltransferase family 61 in Liliopsida (Monocot): the story of a gene family expansion. *Front Plant Sci* 2018;9:1843.

165. Rennie, EA, Scheller, HV. Xylan biosynthesis. *Curr Opin Biotechnol* 2014;**26**:100-7.

166. Yu, J, Hu, F, Dossa, K, et al. Genome-wide analysis of UDP-glycosyltransferase super family in *Brassica rapa* and *Brassica oleracea* reveals its evolutionary history and functional characterization. *BMC Genom* 2017;**18**(1):474.

167. Hou, J, Ye, N, Dong, ZY, et al. Major chromosomal rearrangements distinguish willow and poplar after the ancestral "Salicoid" genome duplication. *Genome Biol Evol* 2016;8:1868-75.

168. Susek, K, Bielski, WK, Hasterok, R, et al. A first glimpse of wild lupin karyotype variation as revealed by comparative cytogenetic mapping. *Front Plant Sci* 2016;7:1152.

169. Ma, X, Vaistij, FE, Li, Y, et al. A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop. *Plant J* 2021;107:613-28.

170. Bennetzen, JL. Transposable elements, gene creation and genome rearrangement in flowering plants. *Curr Opin Genet Dev* 2005;15(6):621-7.

171. Kalendar, R, Sabot, F, Rodriguez, F, et al. Editorial: mobile elements and plant genome evolution, comparative

analyzes and computational tools. *Front Plant Sci* 2021;12:735134.

172. Lysák, MA, Schubert, I. Mechanisms of chromosome rearrangements. In: J Greilhuber, J Dolezel, J Wendel, editors. *Plant Genome Diversity. Volume 2*. Vienna: Springer; 2013; p. 137-47.

173. Barra, V, Fachinetti, D. The dark side of centromeres: types, causes and consequences of structural abnormalities implicating centromeric DNA. *Nat Commun* 2018;9(1):4340.

174. Manova, V, Gruszka, D. DNA damage and repair in plants—from models to crops. *Front Plant Sci* 2015;6:885.

175. Wood, RD. DNA repaired in eukaryotes. *Annu Rev Biochem* 1996;65:135-67.

176. Nisa, M-U, Huang, Y, Benhamed, M, et al. The plant DNA damage response: signaling pathways leading to growth inhibition and putative role in response to stress conditions. *Front Plant Sci* 2019;10:653.

177. Banana Genome Hub. *Musa beccarii* genome. 2022. https://bana.genome-hub.southgreen.fr/node/50/1440401. Accessed 20 November 2022.

178. Wang, Z, Rouard, M, Droc, G, et al. Supporting data for "Genome assembly of *Musa beccarii* shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes." *GigaScience Database*. 2023. http://dx.doi.org/10.5524/102347.