Abstract
Banana tree bacterial wilt is caused by the Ralstonia solanacearum Moko ecotype. These strains vary in their symptom progression in banana, and are classified as typical Moko variants (phylotype IIA and IIB strains from across Central and South America), Bugtok variant (Philippines), and Sergipe facies (the states of Sergipe and Alagoas, Brazil). This study used comparative genomic and phylogenomic approaches to identify a correlation between the symptom progression of the Moko ecotypes based on the analysis of 23 available genomes. Average nucleotide identity and in silico DNA-DNA hybridization revealed a high correlation (>96% and >78%, respectively) between the genomes of Moko variants. Pan-genome analysis identified 21.3% of inheritable regions between representatives of the typical Moko and Sergipe facies variants, which could be traced to an abundance of exclusive homolog clusters. Moko ecotype genomes shared 1,951 orthologous genes, but representatives with typical symptoms did not display unique orthologues. Moreover, Bugtok disease and Sergipe facies genomes did not share any unique genes, suggesting convergent evolution to a shared symptom progression. Overall, genomic and phylogenomic analyses were insufficient to differentiate the Moko variants based on symptom progression.
Keywords
Genomics; Musa spp; Sergipe facies; Bugtok; symptomatology
Introduction
The Ralstonia solanacearum Moko ecotype was first described by Schomburgk during his travels to British Guiana in 1840-1844, and gained economic importance after it devastated banana plantations of the ‘Moko’ cultivar (Musa ABB, Bluggoe subgroup) on the island of Trinidad in the early Twentieth century (Sequeira, 1998Sequeira L (1998) Bacterial wilt: The missing element in international banana improvement programs. In: Prior P, Allen C and Elphinstone J (eds) Bacterial Wilt Disease: Molecular and Ecological Aspects. INRA Editions, Paris, pp 6-14.). ‘Ecotype’ describes a group highly adapted to its host, and the Moko ecotype is a polyphyletic group intimately associated with the Musa spp. (Fegan and Prior, 2006Fegan M and Prior PB (2006) Diverse members of the Ralstonia solanacearum species complex cause bacterial wilts of banana. Australas Plant Pathol 35:93-101.; Ailloud et al., 2015Ailloud F, Lowe T, Cellier G, Roche D, Allen C and Prior P (2015) Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity. BMC Genomics 16:270.).
The pathological outcome caused by the Moko ecotype depends on which of the three variants is involved: Moko, Bugtok, or Sergipe facies (Albuquerque et al., 2014Albuquerque GMR, Santos LA, Felix KCS, Rollemberg CL, Silva AMF, Souza EB, Cellier G, Prior P and Mariano RLR (2014) Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil extend known diversity in paraphyletic phylotype II. Phytopathology 104:1175-1182.; Blomme et al., 2017Blomme G, Dita M, Jacobsen KS, Vicente LP, Molina A, Walter O, Poussier S and Prior P (2017) Bacterial diseases of bananas and enset: Current state of knowledge and integrated approaches toward sustainable management. Front Plant Sci 8:1290.). The typical symptoms of Moko, which is the most widespread variant, begin in the rhizomes and then move towards the pseudostem, vascular darkening is observed in the central region, the leaves turn yellow and wilt, the fruits become deformed (Albuquerque et al., 2014Albuquerque GMR, Santos LA, Felix KCS, Rollemberg CL, Silva AMF, Souza EB, Cellier G, Prior P and Mariano RLR (2014) Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil extend known diversity in paraphyletic phylotype II. Phytopathology 104:1175-1182.), and finally the whole plant wilts (Sequeira, 1998Sequeira L (1998) Bacterial wilt: The missing element in international banana improvement programs. In: Prior P, Allen C and Elphinstone J (eds) Bacterial Wilt Disease: Molecular and Ecological Aspects. INRA Editions, Paris, pp 6-14.). In the Philippines, the sequevar IIB-3 of R. solanacearum attacks the cultivars ‘Saba’ (Musa BBB) and ‘Cardaba’ (Musa ABB), causing a disease called Bugtok or Tapurok (Roperos, 1965Roperos NI (1965) Note on the occurrence of a new disease of cooking banana in the Phillipines. Coffee Cacau J 8:135-136.; Soguilon et al., 1995Soguilon CE, Magnaye LV, Natural MP (1995) Bugtok disease of banana. INIBAP 33:4. ), with infection starting in the inflorescence due to inoculation by vector insects (Blomme et al., 2017). Bugtok symptoms are restricted to the stalk, rachis, and fruits, together with an occasional reddish-brown discoloration of the vascular tissue of the pseudostem, which rarely extends to the rhizome (Blomme et al., 2017). Sergipe facies occurs only in northeastern Brazil, in the states of Sergipe and Alagoas, where it is associated with sequevar IIA-53 of R. solanacearum. The symptomatological picture is similar to that of Bugtok, but the fruits are uneven, ripen prematurely (Albuquerque et al., 2014Albuquerque GMR, Santos LA, Felix KCS, Rollemberg CL, Silva AMF, Souza EB, Cellier G, Prior P and Mariano RLR (2014) Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil extend known diversity in paraphyletic phylotype II. Phytopathology 104:1175-1182.), and display external necrosis (Albuquerque et al., 2021aAlbuquerque GMR, Silva AMF, Silva JRD, Melo EA, Mariano RLR, Lemos MC, Ferraz E and Souza EB (2021a) Sequevar distribution of Ralstonia spp. in Solanaceae in the semiarid climate of the Pernambuco state, Brazil. Eur J Plant Pathol159:13-25.,bAlbuquerque GMR, Silva AMF, Pais AKL, Santos LVSD, Gama MASD and Souza EB (2021b) Avanços na identificação, variabilidade e genômica comparativa de Ralstonia solanacearum ecotipo Moko no Brasil. RAPP 27:118-137.), although both are a result of inoculation by insects vectors.
In spite of Moko ecotype identification (Albuquerque et al., 2014Albuquerque GMR, Santos LA, Felix KCS, Rollemberg CL, Silva AMF, Souza EB, Cellier G, Prior P and Mariano RLR (2014) Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil extend known diversity in paraphyletic phylotype II. Phytopathology 104:1175-1182., 2021aAlbuquerque GMR, Silva AMF, Silva JRD, Melo EA, Mariano RLR, Lemos MC, Ferraz E and Souza EB (2021a) Sequevar distribution of Ralstonia spp. in Solanaceae in the semiarid climate of the Pernambuco state, Brazil. Eur J Plant Pathol159:13-25.) and genomic characterization of representative strains by Silva et al. (2020Silva JR, Pais AKL, Albuquerque GMR, Silva AMF, Silva Junior WJ, Balbino VQ, Fonseca MEN, Gama MAS, Souza EB and Mariano RLR (2020) Genomic sequencing of two isolates of Ralstonia solanacearum causing Sergipe facies and comparative analysis with Bugtok disease isolates. Genet Mol Biol 43:e20200155. ) and Pais et al. (2021Pais AKL, Silva JR, Santos LVS, Albuquerque GMR, Farias ARG, Silva Junior WJ, Balbino VQ, Silva AMF, Gama MAS and Souza EB (2021) Genomic sequencing of different sequevars of Ralstonia solanacearum belonging to the Moko ecotype. Genet Mol Biol 44:e20200172.), no genomic studies have looked the different symptomatology caused by the R. solanacearum Moko ecotype. Therefore, the objective of this study was to investigate the representatives of the Moko ecotype that induce different symptoms in banana trees using comparative genomics and phylogenomic approaches.
Material and Methods
Genomic sequences
The genomes of R. solanacearum available from the National Center for Biotechnology Information (NCBI) at January 2021 were filtered to obtain a final dataset containing only representatives of the three symptomatological variants of the Moko ecotype. These included the ‘M’ group causing typical Moko symptoms (19 genomes of sequevars IIA-6, IIA-24, IIB-3, IIB-4, and IIB-25), the ‘B’ group causing Bugtok symptoms (two genomes of sequevar IIB-3), and the ‘SE’ group causing Sergipe facies (two genomes of sequevar IIA-53). The origin and characteristics of the genomes used in this study are listed in Table 1. Those defined as Bugtok genomes are all representatives of the phylotype/sequevar IIB-3 that originated in the Philippines.
Comparative genomics of the R. solanacearum Moko ecotype
The average nucleotide identity (ANI) and tetranucleotide signature (TETRA) correlation indices were calculated using the pyani 0.2.7 Python3 module (Pritchard et al., 2016Pritchard L, Glover RH, Humphris S, Elphinstone JG and Tothc IK (2016) Genomics and taxonomy in diagnostics for food security: Soft-rotting enterobacterial plant pathogens. Anal Methods 8:12-24.). The ANI was calculated by global alignment of the MUMmer algorithm (ANIm; Kurtz et al., 2004Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C and Salzberg SL (2004) Versatile and open software for comparing large genomes. Genome Biol 5:R12. ). In silico DNA-DNA hybridization (isDDH) values were calculated using the Genome-to-Genome Distance CalculatorGenome-to-Genome Distance Calculator, Genome-to-Genome Distance Calculator, http://ggdc.dsmz.de/ggdc.php (accessed 9 December 2021)
http://ggdc.dsmz.de/ggdc.php ...
platform 2.1 (Meier-Kolthoff et al., 2013Meier-Kolthoff JP, Auch AF, Klenk HP and Goker M (2013) Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics 14:60.) by applying formula 2 for incomplete genomes. Accordingly, isDDH estimates were based on identities/high-scoring pair length. The similarity matrices obtained by ANIm and isDDH were converted into a heatmap using the Morpheus platform.Morpheus platform, Morpheus platform, https://software.broadinstitute.org/morpheus (accessed 13 December 2021)
https://software.broadinstitute.org/morp...
Pan-genome and phylogenetic analysis of the R. solanacearum Moko ecotype
Pan-genome analysis of strains was performed in Roary v. 3.13.0 (Page et al., 2015Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, Fookes M, Falush D, Keane JA and Parkhill J (2015) Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 31:3691-3693.) using genome sequences obtained from RefSeq/NCBI. The resulting orthologous genes were classified as core (genes common to all genomes), softcore (genes contained in 95% of genomes), shell (moderately conserved genes present in various genomes), and clouds (rare genes present in only a few genomes) according to the default settings of the software. The functions and descriptions of gene clusters were acquired from the UniProt platformUniProt, UniProt, https://www.uniprot.org/ (accessed 16 December 2021)
https://www.uniprot.org/ ...
(The UniProt Consortium, 2021The UniProt Consortium (2021) UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res 49:D480-D489.). The core gene pool was automatically aligned using MAFFT v. 7.3102 (Katoh and Standley, 2013Katoh K and Standley DM (2013) MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol 30:772-780.) and implemented in Roary using the -mafft flag. This created a multi-FASTA nucleotide sequence alignment of all core genes. The phylogenomic tree of the core genes was constructed using multi-FASTA alignment with the maximum likelihood method in IQ-TREE v. 2.0.4 (Nguyen et al., 2015Nguyen LT, Schmidt HA, Von Haeseler A and Minh BQ (2015) IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268-274.). ModelFinder (Kalyaanamoorthy et al., 2017Kalyaanamoorthy S, Minh BQ, Wong TKF, Haeseler AV and Jermiin LS (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat Methods 14:587-589.) was employed to select the best evolutionary model. Node support was determined by ultrafast bootstrap (Minh et al., 2013Minh BQ, Nguyen MAT and Von Haeseler A (2013) Ultrafast approximation for phylogenetic bootstrap. Mol Biol Evol 30:1188-1195. ) with 100,000 repetitions. The maximum likelihood tree was viewed in Figtree v.1.4.4 (Rambaut, 2009Rambaut A (2009) FigTree v1.3.1, 1, http://tree.bio.ed.ac.uk (accessed 9 December 2021)
http://tree.bio.ed.ac.uk ...
). IQ-TREE v. 2.0.4 was used to construct the phylogeny matrix summarizing the presence/absence of a gene with the same bootstrap configuration and following the same process described earlier for visualization.
Results
Genomic sequences
The analyses were performed with the 23 genomes available in this group, however some limitations were found regarding the level of assembly and representation of Bugtok disease and Sergipe facies, since both contained only two representatives (Table 1).
Comparative genomics of the R. solanacearum Moko ecotype
ANIm analysis revealed an average similarity of 97.6% among the 23 genomes of the R. solanacearum Moko ecotype (Table 2). A strong similarity was observed in the Moko genomes of the three symptomatic variants (Figure S1 Figure S1 - Heatmap generated from in silico DNA-DNA hybridization and average nucleotide identity using the MUMmer algorithm of Ralstonia solanacearum Moko ecotype genomes and its symptomatological variants. *The upper triangle refers to the values of isDDH and the lower triangle to the values of ANIm. ). Specifically, ‘M’ genomes presented 97.3% similarity with ‘B’ and 97.4% with ‘SE’ genomes; whereas ‘B’ genomes presented 96.3% similarity with their ‘SE’ counterparts. TETRA values confirmed the elevated similarity (99.9%-100%) between the genomes (Table 2).
Based on isDDH results, the genomes of the Moko ecotype presented an average similarity of 78.8%, with a variation of 66.1%-100% (Table 2). The mean values for phylotypes IIA and IIB were 94% and 87.8%, respectively, with a similarity of 67.1% among them. Group ‘M’ presented 73.2% similarity, whereas groups ‘B’ and ‘SE’ exhibited more than 99% similarity. Intergroup comparison revealed 77.2% and 67.3% similarity between ‘M’ and ‘B’ or ‘SE’ groups, respectively, as well as 67.7% similarity between ‘B’ and ‘SE’ groups (Table 2).
A heatmap was constructed with the means calculated by the ANIm and isDDH of the Moko ecotype (Figure S1 Figure S1 - Heatmap generated from in silico DNA-DNA hybridization and average nucleotide identity using the MUMmer algorithm of Ralstonia solanacearum Moko ecotype genomes and its symptomatological variants. *The upper triangle refers to the values of isDDH and the lower triangle to the values of ANIm. ). Accordingly, the strains were grouped into two large phylotypes (IIA and IIB), which were further divided into four subgroups: IIA(α), IIA(β), IIB(α), and IIB(β). ‘M’ genomes were present mainly in subgroups IIA(α) and IIB(β), although some genomes of this group were detected also in other subgroups. ‘B’ genomes displayed significant similarity and constituted a group together with some ‘M’ genomes. Similar results were observed for ‘SE’ genomes.
Pan-genome and multilocus sequence analysis of the R. solanacearum Moko ecotype
Pan-genome analyses revealed the presence of 9,164 clusters of information, of which 1,951 were identified as core genes (Figure 1A and Table S1 Table S1 - Core gene clusters obtained by pan-genome analysis of the Ralstonia solanacearum Moko ecotype and its symptomatological variants. ). A total of 3,308 clusters were found among representatives with typical Moko symptoms, but none were common to all representatives (Figure 1B and Table S2 Table S2 - Clusters obtained by pan-genome analysis of the Ralstonia solanacearum Moko ecotype and its symptomatological variants. ). ‘B’ genomes contained 135 clusters, of which two insertions and two unknown sequences were detected in this group. (Tables S2 and S3). The identified gene clusters included the transposase families IS3 and IS5, as well as two proteins of unknown function (one containing the domain DUF4158 plus a hypothetical protein). ‘SE’ genomes revealed 113 clusters, of which 60 were present in all representatives of this group and were related to biological processes, molecular functions, relationship with binding molecules, the type three secretion system (T3SS), insertion sequences, and CRISPR (Tables S2 and S3).
Pan-genome representation of the Ralstonia solanacearum Moko ecotype generated by Roary software. (A) Gene categories (core, softcore, shell, and cloud) present in genomes were identified with 90% percent identity. (B) Venn diagram showing clusters present in the genome of R. solanacearum Moko ecotype strains and its symptomatological variants. Strains of Moko ecotype causing Bugtok symptoms (CIP417 and Molk2), Sergipe facies (IBSBF2570 and SFC), and typical symptoms of Moko (all other strains) are shown.
‘M’ and ‘B’ genomes shared 15 unique clusters, while ‘M’ and ‘SE’ genomes shared only one cluster (Figure 1B and Tables S2 and S3). The information shared between the genomes associated with typical Moko and Bugtok symptoms related to biological processes, molecular functions, and cellular components. The genome associated with Sergipe facies symptoms shared information related to molecular function and/or biological processes. No clusters were shared exclusively between the Bugtok and Sergipe facies genomes (Figure 1B and Table S3 Table S3 - Exclusive and shared gene clusters obtained by pan-genome analysis of the Ralstonia solanacearum Moko ecotype and its symptomatological variants. ).
The phylogenomic trees showed strong bootstrap support in all branches, indicating a robust phylogeny for R. solanacearum. As indicated by ANIm and isDDH, the phylogenomic tree based on core genes distinguished clearly the two phylotypes, IIA and IIB, and the four subgroups IIA(α), IIA(β), IIB(α), and IIB(β) (Figure 2A). However, inference based on the presence and absence of gene clusters could not group the genomes of the Moko ecotype in the same clusters as done previously in this work, because the genomes IBSBF2570 and SFC (phylotype IIA) did not group with other representatives of subgroup IIA(β) (Figure 2B).
(A) Maximum likelihood phylogenetic tree of core gene sequences annotated for the Ralstonia solanacearum Moko ecotype genomes and its symptomatological variants. (B) Maximum likelihood phylogenetic tree based on the presence or absence of orthologous genes of R. solanacearum Moko ecotype genomes and its symptomatological variants. * Subgroup IIA(α) is shown in green, IIA(β) in black, IIB(α) in blue, and IIB(β) in red. Strains of Sergipe facies (IBSBF2570 and SFC, IIA- β) are shown yellow. Strains of Moko ecotype causing Bugtok symptoms (CIP417 and Molk2), Sergipe facies (IBSBF2570 and SFC), and typical symptoms of Moko (all other strains) are shown.
Discussion
Based on these results, the ANIm and isDDH analyses were not sufficient to discriminate between strains with different Moko ecotype symptoms (Figure S1 Figure S1 - Heatmap generated from in silico DNA-DNA hybridization and average nucleotide identity using the MUMmer algorithm of Ralstonia solanacearum Moko ecotype genomes and its symptomatological variants. *The upper triangle refers to the values of isDDH and the lower triangle to the values of ANIm. ). However, the genomes of Bugtok disease and Sergipe facies variants presented high similarity between their representatives, proving a high genomic homogeneity within these groups (Table 2). Thus, the ANIm and isDDH analyses proved that the similarity of symptoms did not correspond with greater genomic proximity between organisms.
However, the pan-genome of the strain of typical Moko symptoms, Bugtok disease, and Sergipe facies variants shared 1,951 homologous genes, or 21.3% of inheritable regions for the entire Moko ecotype. Based on this result and the hypotheses proposed by Ailloud et al. (2015Ailloud F, Lowe T, Cellier G, Roche D, Allen C and Prior P (2015) Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity. BMC Genomics 16:270.), we conclude that the Moko ecotype may have inherited pathogenic traits from a recent common ancestor by sharing some homologous genes. Ailloud et al. (2015Ailloud F, Lowe T, Cellier G, Roche D, Allen C and Prior P (2015) Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity. BMC Genomics 16:270.) evaluated groups of pathogens highly adapted to hosts, which included representatives of the R. solanacearum Moko strains. These strains were characterized by the absence of exclusive homologous regions, suggesting that this ecotype might have arisen from the convergent evolution of several strains, which led to the ability to infect banana trees.
For R. solanacearum representatives with typical Moko symptoms, no clusters of homologous genes were found. Instead, Sergipe facies representatives contained the largest number (113) of exclusive clusters, of which 51.1% were observed in both strains, indicating greater diversity within the group. This observation may be related to the high rate of mutation necessary to ensure the prevalence of this trait in the environment, considering that it is the most recently reported symptomatological group (Albuquerque et al., 2014Albuquerque GMR, Santos LA, Felix KCS, Rollemberg CL, Silva AMF, Souza EB, Cellier G, Prior P and Mariano RLR (2014) Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil extend known diversity in paraphyletic phylotype II. Phytopathology 104:1175-1182.; Silva et al., 2020Silva JR, Pais AKL, Albuquerque GMR, Silva AMF, Silva Junior WJ, Balbino VQ, Fonseca MEN, Gama MAS, Souza EB and Mariano RLR (2020) Genomic sequencing of two isolates of Ralstonia solanacearum causing Sergipe facies and comparative analysis with Bugtok disease isolates. Genet Mol Biol 43:e20200155. ). Within the exclusive clusters, 13.3% were associated with T3SS and 10% with insertion sequences, which could favor various genetic rearrangements. R. solanacearum uses a T3SS to deliver effector proteins, which manipulate the host physiology to increase pathogen success. Insertion sequences are mobile genetic elements that are commonly present in bacteria. They can lead to gene activation or repression, as well as DNA rearrangements, resulting in deletions, inversions, and amplification of genes (Chandler and Mahillon, 2002Chandler M and Mahillon J (2002) Insertion Sequences Revisited. In: Craig N, Craigie R, Gellert M and Lambowitz A (eds), Mobile DNA II. ASM Press, Washington, pp 305-366.). These insertion sequences demonstrate the process of horizontal transfer of genetic information, known to be an important mechanism for the evolution of the bacterial genome, as evidenced in Blood Disease Bacterium (BDB; Remenant et al., 2011Remenant B, Cambiaire J-C, Cellier G, Jacobs JM, Mangenot S, Barbe V, Lajus A, Vallenet D, Medigue C, Fegan M et al. (2011) Ralstonia syzygii, the Blood Disease Bacterium and some Asian R. solanacearum strains form a single genomic species despite divergent lifestyles. PLoS One 6:e24356.), corroborating this process in the genomes of Sergipe facies.
Strains with typical Moko and Bugtok symptoms not were differentiated by phylogenetic analyses of the endoglucanase (egl) gene (Fegan and Prior, 2006Fegan M and Prior PB (2006) Diverse members of the Ralstonia solanacearum species complex cause bacterial wilts of banana. Australas Plant Pathol 35:93-101.). The sequence of the egl genes is used to determine a phylogenetic relationship among isolates of R. solanacearum, differentiating them by sequevar. However, the results obtained in the current study identified four unique ortholog clusters in strains representative of Bugtok symptoms, which can be used to distinguish these two symptomatological conditions. The identified gene clusters included two transposase families, as well as two proteins of unknown function. Various insertion sequence families have been identified among R. solanacearum, most of them are scattered throughout the single strains. In addition, closely related strains tend to have similar insertion sequence patterns (Gonçalves et al., 2020Gonçalves OS, Campos KF, Assis JCS, Fernandes AS, Souza TS, Rodrigues LGC, Queiroz MV and Santana MF (2020) Transposable elements contribute to the genome plasticity of Ralstonia solanacearum species complex. Microb Genom 6:e000374.).
Phylogenetic analysis of the genus Ralstonia successfully distinguished strains from phylotypes IIA and IIB (Zhang and Qiu, 2016Zhang Y and Qiu S (2016) Phylogenomic analysis of the genus Ralstonia based on 686 single-copy genes. Antonie Van Leeuwenhoek 109:71-82.). Phylotype IIA has been reported to be highly recombinogenic and diverse, with ongoing species expansion. In contrast, multilocus sequence analysis of nine loci has suggested the almost clonal character of phylotype IIB (Wicker et al., 2012Wicker E, Lefeuvre P, Cambiaire JC, Lemaire C, Poussier S and Prior P (2012) Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA. ISME J 6:961-974. ). The description of phylotype IIA as recombinogenic and diverse may hint at the behavior of the genomes of Sergipe facies strains (IIA-53) observed in both phylogenomic analyses (Figure 2B) performed in the present study. In contrast, Bugtok strains (IIB-3) showed strong genetic similarity, which seems to confirm previous descriptions of phylotype IIB by Wicker et al. (2012Wicker E, Lefeuvre P, Cambiaire JC, Lemaire C, Poussier S and Prior P (2012) Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA. ISME J 6:961-974. ). Finally, it is important to highlight that the analyses carried out in this study may have limitations due to the low quality of some genomes and the discrepancy in the number of representatives of the symptomatological group.
Conclusions
The genome of the R. solanacearum Moko ecotype are strains that present high genomic similarity, chiefly among variants expressing Sergipe facies and Bugtok symptoms. Here, pan-genome analysis identified 21.3% of inheritable regions among representatives of the Moko ecotype, and the symptomatological Sergipe facies variant stood out for presenting the largest number of clusters of exclusive homologues. Accordingly, the similarity among symptomatic cases of Bugtok disease and Sergipe facies does not correspond with genomic or phylogenomic properties. Other approaches, possibly focusing on pathogenicity, virulence, and ecological factors, should be employed to determine a common denominator of different Moko ecotype symptoms. For example, little is known about the interactions between bacterial strains and insects responsible for the transmission of Bugtok and Sergipe facies. Knowing that both diseases occur by infection via inflorescence, this point may help understand the peculiarities of Moko pathogenesis and its symptoms.
Acknowledgements
We thank the Conselho Nacional de Ciência e Tecnologia, Brazil and Fundação de Amparo a Ciência e Tecnologia do Estado de Pernambuco, Brazil. We also thank the Laboratory of Bioinformatics and Evolutionary Biology of Universidade Federal de Pernambuco for making their servers available for bioinformatics analysis.
References
- Ailloud F, Lowe T, Cellier G, Roche D, Allen C and Prior P (2015) Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity. BMC Genomics 16:270.
- Albuquerque GMR, Santos LA, Felix KCS, Rollemberg CL, Silva AMF, Souza EB, Cellier G, Prior P and Mariano RLR (2014) Moko Disease-Causing Strains of Ralstonia solanacearum from Brazil extend known diversity in paraphyletic phylotype II. Phytopathology 104:1175-1182.
- Albuquerque GMR, Silva AMF, Silva JRD, Melo EA, Mariano RLR, Lemos MC, Ferraz E and Souza EB (2021a) Sequevar distribution of Ralstonia spp. in Solanaceae in the semiarid climate of the Pernambuco state, Brazil. Eur J Plant Pathol159:13-25.
- Albuquerque GMR, Silva AMF, Pais AKL, Santos LVSD, Gama MASD and Souza EB (2021b) Avanços na identificação, variabilidade e genômica comparativa de Ralstonia solanacearum ecotipo Moko no Brasil. RAPP 27:118-137.
- Blomme G, Dita M, Jacobsen KS, Vicente LP, Molina A, Walter O, Poussier S and Prior P (2017) Bacterial diseases of bananas and enset: Current state of knowledge and integrated approaches toward sustainable management. Front Plant Sci 8:1290.
- Chandler M and Mahillon J (2002) Insertion Sequences Revisited. In: Craig N, Craigie R, Gellert M and Lambowitz A (eds), Mobile DNA II. ASM Press, Washington, pp 305-366.
- Fegan M and Prior PB (2006) Diverse members of the Ralstonia solanacearum species complex cause bacterial wilts of banana. Australas Plant Pathol 35:93-101.
- Gonçalves OS, Campos KF, Assis JCS, Fernandes AS, Souza TS, Rodrigues LGC, Queiroz MV and Santana MF (2020) Transposable elements contribute to the genome plasticity of Ralstonia solanacearum species complex. Microb Genom 6:e000374.
- Kalyaanamoorthy S, Minh BQ, Wong TKF, Haeseler AV and Jermiin LS (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat Methods 14:587-589.
- Katoh K and Standley DM (2013) MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol 30:772-780.
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C and Salzberg SL (2004) Versatile and open software for comparing large genomes. Genome Biol 5:R12.
- Meier-Kolthoff JP, Auch AF, Klenk HP and Goker M (2013) Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics 14:60.
- Minh BQ, Nguyen MAT and Von Haeseler A (2013) Ultrafast approximation for phylogenetic bootstrap. Mol Biol Evol 30:1188-1195.
- Nguyen LT, Schmidt HA, Von Haeseler A and Minh BQ (2015) IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268-274.
- Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, Fookes M, Falush D, Keane JA and Parkhill J (2015) Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 31:3691-3693.
- Pais AKL, Silva JR, Santos LVS, Albuquerque GMR, Farias ARG, Silva Junior WJ, Balbino VQ, Silva AMF, Gama MAS and Souza EB (2021) Genomic sequencing of different sequevars of Ralstonia solanacearum belonging to the Moko ecotype. Genet Mol Biol 44:e20200172.
- Pritchard L, Glover RH, Humphris S, Elphinstone JG and Tothc IK (2016) Genomics and taxonomy in diagnostics for food security: Soft-rotting enterobacterial plant pathogens. Anal Methods 8:12-24.
- Remenant B, Cambiaire J-C, Cellier G, Jacobs JM, Mangenot S, Barbe V, Lajus A, Vallenet D, Medigue C, Fegan M et al (2011) Ralstonia syzygii, the Blood Disease Bacterium and some Asian R. solanacearum strains form a single genomic species despite divergent lifestyles. PLoS One 6:e24356.
- Roperos NI (1965) Note on the occurrence of a new disease of cooking banana in the Phillipines. Coffee Cacau J 8:135-136.
- Sequeira L (1998) Bacterial wilt: The missing element in international banana improvement programs. In: Prior P, Allen C and Elphinstone J (eds) Bacterial Wilt Disease: Molecular and Ecological Aspects. INRA Editions, Paris, pp 6-14.
- Silva JR, Pais AKL, Albuquerque GMR, Silva AMF, Silva Junior WJ, Balbino VQ, Fonseca MEN, Gama MAS, Souza EB and Mariano RLR (2020) Genomic sequencing of two isolates of Ralstonia solanacearum causing Sergipe facies and comparative analysis with Bugtok disease isolates. Genet Mol Biol 43:e20200155.
- Soguilon CE, Magnaye LV, Natural MP (1995) Bugtok disease of banana. INIBAP 33:4.
- The UniProt Consortium (2021) UniProt: The universal protein knowledgebase in 2021. Nucleic Acids Res 49:D480-D489.
- Wicker E, Lefeuvre P, Cambiaire JC, Lemaire C, Poussier S and Prior P (2012) Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA. ISME J 6:961-974.
- Zhang Y and Qiu S (2016) Phylogenomic analysis of the genus Ralstonia based on 686 single-copy genes. Antonie Van Leeuwenhoek 109:71-82.
Internet Resources
- Genome-to-Genome Distance Calculator, Genome-to-Genome Distance Calculator, http://ggdc.dsmz.de/ggdc.php (accessed 9 December 2021)
» http://ggdc.dsmz.de/ggdc.php - Morpheus platform, Morpheus platform, https://software.broadinstitute.org/morpheus (accessed 13 December 2021)
» https://software.broadinstitute.org/morpheus - Rambaut A (2009) FigTree v1.3.1, 1, http://tree.bio.ed.ac.uk (accessed 9 December 2021)
» http://tree.bio.ed.ac.uk - UniProt, UniProt, https://www.uniprot.org/ (accessed 16 December 2021)
» https://www.uniprot.org/
Supplementary material
The following online material is available for this article:
Edited by
Associate Editor
Publication Dates
-
Publication in this collection
02 Dec 2022 -
Date of issue
2022
History
-
Received
21 Feb 2022 -
Accepted
14 Oct 2022