Next Article in Journal
Sodium Polyoxotungstate Inhibits the Replication of Influenza Virus by Blocking the Nuclear Import of vRNP
Next Article in Special Issue
Phytoplasma-Associated Diseases in South America: Thirty Years of Research
Previous Article in Journal
Exploring Gut Microbiota Alterations with Trimethoprim-Sulfamethoxazole and Dexamethasone in a Humanized Microbiome Mouse Model
Previous Article in Special Issue
Paulownia Witches’ Broom Disease: A Comprehensive Review
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 12
Issue 5
10.3390/microorganisms12051016
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Divergence within the Taxon ‘Candidatus Phytoplasma asteris’ Confirmed by Comparative Genome Analysis of Carrot Strains

by
Rafael Toth
1,
Anna-Marie Ilic
1,
Bruno Huettel
2,
Bojan Duduk
3 and
Michael Kube
1,*
1
Department of Integrative Infection Biology Crops-Livestock, University of Hohenheim, 70599 Stuttgart, Germany
2
Max Planck Genome Center, 50829 Cologne, Germany
3
Institute of Pesticides and Environmental Protection, 11080 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Microorganisms 2024, 12(5), 1016; https://doi.org/10.3390/microorganisms12051016
Submission received: 19 March 2024 / Revised: 30 April 2024 / Accepted: 14 May 2024 / Published: 17 May 2024
(This article belongs to the Special Issue Phloem Localized Insect Transmitted Bacteria Associated with Plant Diseases)
Browse Figures
Versions Notes

Abstract

:
Phytoplasmas are linked to diseases in hundreds of economically important crops, including carrots. In carrots, phytoplasmosis is associated with leaf chlorosis and necrosis, coupled with inhibited root system development, ultimately leading to significant economic losses. During a field study conducted in Baden-Württemberg (Germany), two strains of the provisional taxon ‘Candidatus Phytoplasma asteris’ were identified within a carrot plot. For further analysis, strains M8 and M33 underwent shotgun sequencing, utilising single-molecule-real-time (SMRT) long-read sequencing and sequencing-by-synthesis (SBS) paired-end short-read sequencing techniques. Hybrid assemblies resulted in complete de novo assemblies of two genomes harboring circular chromosomes and two plasmids. Analyses, including average nucleotide identity and sequence comparisons of established marker genes, confirmed the phylogenetic divergence of ‘Ca. P. asteris’ and a different assignment of strains to the 16S rRNA subgroup I-A for M33 and I-B for M8. These groups exhibited unique features, encompassing virulence factors and genes, associated with the mobilome. In contrast, pan-genome analysis revealed a highly conserved gene set related to metabolism across these strains. This analysis of the Aster Yellows (AY) group reaffirms the perception of phytoplasmas as bacteria that have undergone extensive genome reduction during their co-evolution with the host and an increase of genome size by mobilome.

1. Introduction

Bacteria of the provisional taxon ‘Candidatus Phytoplasma’ are phloem-limited parasites [1] that infect a wide range of plant species, including many important crops as well as ornamentals and wild hosts. These parasites of the Mollicutes class are transmitted by phloem-sucking insect vectors in general [2,3], by vegetative propagation and in some rare cases by seed transmission [4,5]. Phytoplasmosis often manifests as leaf yellowing, hindered organogenesis, growth abnormalities, and general plant decline. As a consequence, this can lead to drastic yield losses in agriculture [6]. This is also true for phytoplasmosis of carrots (Daucus carota ssp. sativus L.), which has been reported in various cropping areas worldwide. Molecular screenings using PCR-based methods and phylogenetic analysis show that infections in European carrot production areas are mainly associated with ‘Ca. P. asteris’ strains that are also known as aster yellows (AY) phytoplasmas [7,8,9,10,11,12].
Infected carrot plants exhibit symptoms such as yellowing, reddening, necrosis of leaves, proliferation, and reduced taproot size. Such symptoms have been linked to the secretion of several phytoplasma effector proteins [13,14,15,16]. ‘Ca. P. asteris’- caused diseases have been reported in more than 300 species in 38 families of plants [3], and this taxon is the first and most extensively analysed phytoplasma species in terms of genome and pathogen–host interaction. Among the analysed genomes are ‘Ca. P. asteris’ strains infecting various hosts, including onion, lettuce [17,18], maize [19], grapevine (acc. no. CP035949), rapeseed [20], paulownia [21], and mulberry [22], but despite its importance, no carrot-associated strain has been characterised to date.
Herein, we provide new insights into the AY group by analysing the genomes of the two ‘Ca. P. asteris’ strains M8 and M33, were obtained during a field study in a carrot cultivation area in Baden-Württemberg, Germany. The complete genomes of the strains were reconstructed and compared with respect to their phylogenetic assignment and the common and distinct genetic repertoires of the previously described complete genomes of the AY group. Further, our results provide complex insights into the pathogen–host interaction and dependency of the phytoplasmas causing aster yellows disease.

2. Materials and Methods

2.1. Plant Material and DNA Extraction

Phytoplasmosis of carrots associated with ‘Ca. P. asteris’ strains are frequently ob-served in the south of Germany [23]. In August 2019, symptomatic carrots (Daucus carota subsp. sativus) were collected from a single carrot plot with 0.5% symptomatic plants in Langenau (Baden-Württemberg, Germany). Symptomatic carrots were collected and individually tested for ‘Ca. P. asteris’ by PCR and confirmed by sequencing the amplicons [24]. Two samples were selected for genome sequencing (Figure 1). Total DNA was purified using the cetyltrimethylammonium bromide (CTAB) plant DNA extraction protocol [25]. The DNA concentrations of the extracted DNA samples were quantified using a Qubit fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). Infection was detected through endpoint PCR utilising the universal primers P1 and P7 for partial amplification of the phytoplasma rRNA operon [26,27].

2.2. Genome Sequencing

Shotgun sequencing of extracted DNA from infected carrot samples M8 and M33 (Figure 1) was performed separately for each sample using two sequencing technologies. To generate short reads with paired ends, Illumina sequencing [28] was carried out on the HiSeq 2500 platform (Illumina, San Diego, CA, USA). For long-read sequencing, single-molecule real-time sequencing (SMRT) sequencing [29] was performed on the Sequel IIe platform (Pacific Biosciences, Menlo Park, CA, USA). First, DNA was enriched for longer fragments through a 0.45% (v/v) PB AMPure bead purification step (Pacific Biosciences). Barcoded libraries were then prepared according to the manufacturer’s protocol “Preparing HiFi Libraries from Low DNA Input Using SMRTbell Express Template Prep Kit 2.0” (Pacific Biosciences). Libraries were equimolarly pooled and sequenced on a Sequel II device (Pacific Biosciences) using a Sequel II binding kit 2.0, Sequel II sequencing chemistry 2.0, and an 8M ZMW SMRT cell for 30 h (Pacific Biosciences). Sequencing data were demultiplexed and high-fidelity data were generated using the SMRTlink Suite v.9.0 (Pacific Biosciences) with default settings. The NGS approaches were conducted by the Max Planck Genome Centre Cologne (Cologne, Germany).

2.3. Hybrid Genome Assembly and Quality Assessment

The following analyses were also carried out individually for each strain. The short reads derived from each Illumina sequencing were mapped to the genome of ‘Ca. P. asteris’ strain RP166, which served as a reference to assign the reads obtained from the metagenomic DNA templates. Mapping and extraction of the reads were performed using the short-read to reference genome mapping tool as part of the CLC Genomic Workbench v. 22. (QIAGEN, Aarhus, Denmark) with default mapping parameters.
To reconstruct the chromosomes and plasmids of ‘Ca. P. asteris’ strains M8 and M33, the selection of mapped Illumina read pairs along with all generated SMRT reads were corrected, trimmed, and assembled to contigs using Canu assembler v. 1.9 [30]. The number of incorporated reads and the calculation of the sequencing coverage were taken from the Canu report files. Default correction and trimming parameters were used, while hybrid assembly was performed using the parameters ‘haplotype’ for Illumina reads and ‘pacbio-corrected’ for SMRT reads, setting an estimated genome size of 0.8 Mb. In order to identify phytoplasma replicons, all contigs were sorted by length using seqtk v. 1.3 (https://github.com/lh3/seqtk (accessed on 27 June 2022)) and split into two separate data sets for further evaluation. Contigs >5 kb kilobases were compared against a custom-made subset of protein databases containing all proteins assigned to D. carota subsp. sativus and to the phylum Mycoplasmatota using BLASTX [31]. Contigs <5 kb were compared with the UniProt Reference Cluster 100 (UniREF100) [32] using the DIAMOND high-throughput aligner v. 2.0.15 [33], applying the ‘fast’ parameter to increase the throughput of the larger dataset of contigs <5 kb. BLASTX and DIAMOND outputs were parsed using the Metagenome Analyzer (MEGAN) v. 6 [34] for taxonomic binning with default parameters enabling the identification of phytoplasma chromosomes and plasmids. Sequence coverage was extracted from the assembly output of Canu to assess the reliability and integrity of the assemblies [30]. Overlaps of circular constructs were confirmed using BLASTN [31], with default settings, and manually removed in the Artemis Genome Browser [35]. The chromosome start was set to the gene dnaA based on the cumulative G+C skew minimum, which was calculated in the Artemis Genome Browser [35]. For plasmids, the gene rep was set to position one.
Genome completeness was also estimated by comparative analysis of the protein content using the Benchmarking Universal Single-Copy Orthologs (BUSCO) v. 5.4.6 software [36]. The analysis was performed using the nucleotide sequences of the chromosomes of ‘Ca. P. asteris’ M8 and M33 along with all available complete chromosome sequences of ‘Ca. P. asteris’ and compared with a dataset of 151 conserved orthologues of the class Mollicutes on the Galaxy Europe server (https://usegalaxy.eu (accessed 12 August 2022)). In the absence of any further specifications, the default parameters were used.

2.4. Functional Annotation

The chromosomes and plasmids were functionally annotated using the Rapid Annotation using Subsystem Technology (RAST) v.2.0 pipeline [37]. Annotation parameters were set to translation table eleven and the NCBI taxonomy ID: 85,620 of ‘Ca. P. asteris’. Annotations were manually curated for missing, overlapping, or frameshifted ORFs. This analysis was conducted using BLASTP [31] comparison of the deduced protein sequence against the NRProt (from GenBank Release 251, www.ncbi.nlm.nih.gov (accessed on 19 August 2022)) and UniRef100 [32] databases. Adjustments were made in Artemis [35]. The deduced proteins were functionally characterised and used for metabolic reconstruction by applying BlastKOALA [38], the MetaCyc database v. 23.1 [39] and InterPro v. 86.0 [40]. Structural non-coding RNAs, such as tRNAs and rRNAs, were predicted using the software tRNAscan-SE v2.0.10 [41] and RNAmmer [42]. Membrane and secreted proteins were predicted with Phobius v. 1.01 [43] and InterPro [40]. Finally, the annotations of ‘Ca. P. asteris’ M8 and M33 were evaluated again using the Microbial Genome Submission check service from NCBI (https://www.ncbi.nlm.nih.gov/genomes/frameshifts/frameshifts.cgi (accessed on 6 June 2023)). Default parameters were used unless otherwise noted.

2.5. Phylogenetic and Functional Comparison

The M8 and M33 chromosomes were compared with all available complete chromosome sequences belonging to the taxon ‘Ca. P. asteris’, retrieved from NCBI (taxonomy ID: 85620) in terms of phylogeny and potential functions. For a phylogenetic analysis at the chromosome level, the average nucleotide identity (ANI) [44] was calculated based on a whole chromosome alignment and the neighbour-joining method [45] within the CLC Genomic Workbench v. 22 (QIAGEN Aarhus, Denmark) with default parameters. To verify the ANI results a whole chromosome sequence synteny analysis was conducted with the multiple genome aligner Mauve v. 20150226 [46] and the Artemis Comparison Tool (ACT) [35]. ACT analysis was calculated based on BLASTN [31] comparison M8-format outputs. Moreover, sequence identity analysis at the single gene level was performed with recently published marker genes (namely, 16S rRNA, tufB, groEL, secY, and secA) and the chromosomal region comprising the genes rplV to rpsC, according to current recommendations [47]. For each marker gene, an identity matrix based on multiple sequence alignments was calculated using BioEdit v. 7.2.5 [48] on default settings. The 16S rDNA phylogeny was analysed using Molecular Evolutionary Genetic Analysis (MEGA) v. 10.2 [49] with the maximum likelihood [50] and neighbour-joining method [45] based on a multiple sequence alignment generated in MEGA. Bootstrapping was performed with 1000 iterations. The threshold for a significant score to be included in the phylogenetic tree was set to 70%. OrthoFinder v. 2.5.5 [51] was used for pan-genome analysis, which included predicting orthologs, paralogs, and unique coding sequences (CDS) from assigned orthogroups that included all coding sequences (without pseudogenes) of the AY group members studied. Unassigned CDS were considered as unique CDS. Default parameters have been used unless otherwise stated.

3. Results

3.1. Sequencing and Hybrid Assembly

Illumina sequencing generated 2,116,564 paired-end short reads for strain M8 and 2,168,250 paired-end short reads for strain M33. A total of 223,496 reads for M8 were positively mapped on the genome of the strain RP166, whereas for M33 480,718 reads were selected for genome assembly input (Table 1). SMRT sequencing produced 191,988 long reads for strain M8 and 105,166 long reads for M33. For chromosome assembly, 10,116 reads for strain M8 and 11,776 reads for strain M33 were used. Final phytoplasma chromosome contigs showed a total length of 772,691 bp for M8 and 657,324 bp for M33. Sequencing coverage of the chromosome contigs was 116-fold for M8 and 176-fold for M33.

3.2. Quality Assessment

The final quality assessment of the hybrid system using BUSCO analysis revealed a total of 151 identified orthologues in each of the analysed genome sequences that match orthologues in the BUSCO database (Figure 2). This indicates high quality according to completeness concerning the evaluation of conserved orthologs within Mollicutes and supports membership of the reconstructed chromosomes in the AY group.

3.3. Genomic Benchmarks of the Taxon ‘Ca. P. asteris’

The circular chromosomes of M8 and M33 differ by ~115 kb bp in size (Table 2, Figure 3). M8, with a chromosome size of ~773 kb, is among the larger chromosomes of the AY group (an average of 734.4 kb). For instance, ‘Ca. P. asteris’ Zhengzhou possesses the longest chromosome sequence, measuring roughly 892 kb, while M33, with a chromosome size of 657,324 bp, represents strains with smaller chromosome types. Among these, M3 has the smallest chromosome, with a length of approximately 576 kb. The total number, i.e., 741 CDS, in the M8 genome also differed compared to ‘Ca. P. asteris’ strain M33, which had 595 CDS (Table 2, Figure 3). However, the coding density of 0.958 kb/gene for strain M8 is above average (0.909 kb/gene). Moreover, strain M33 showed the lowest G+C content among the AY group, with a value of 26.8%, in contrast to other AY group members whose G+C content ranged from 27% to 29% (an average of 27.79%). Structural RNAs are encoded including the typical two rRNA operons and 32 tRNAs of phytoplasmas, but also rnpB, ffs, and ssrA [52]. M8 and M33, harbour a plasmid, a feature that is not described for all complete genomes within the AY group. In addition to strains M8 and M33, single plasmids were identified in strains OY-M, De Villa, and Zhengzhou, whereas AYWB had four reported plasmids. The QS2020 strain exhibited the highest coding density within the AY group, with 0.981 CDS/kb (Table 2). In summary, the completely reconstructed genomes of M8 and M33 display typical genomic features associated with the AY group.

3.4. Phylogenetic Analysis

3.4.1. Average Nucleotide Identity

The chromosomes of strains M8, M33, and the other analysed complete chromosomes of the AY group underwent ANI analysis, resulting in the formation of two distinct clusters. Despite originating from the same field, M8 and M33 clustered into different groups and this analysis did not reveal any obvious correlation in the geographical distribution of ‘Ca. P. asteris’ strains (Figure 4).
The strains M33 and AYWB displayed identity values ranging from 92% to 94% when compared to the other chromosomes within the AY group (Table 3), which is less than the 95% threshold recommended for taxon affiliation in recent requirements for the taxonomic revision of the genus ‘Ca. Phytoplasma’ [47].
Analyses highlight the assignment of these two strains to different clusters which may represent separate taxons, which is in accordance with the results of former studies using ANI analysis to investigate genomic divergence among 16SrI phytoplasmas [53]. The division ‘Ca. P. asteris’ strains into two groups via ANI analysis is further supported by blocks of conserved sequence synteny observed in Mauve (Figure 5).

3.4.2. Single Gene Analysis Supports Cluster Formation

ANI analysis revealed identities lower than 95% for the ribosomal cluster 16SrI-A comprising the strains M33 and AYWB compared to the other AY group members, which indicates a questionable affiliation of this cluster to the provisional taxon ‘Ca. P. asteris’. The phylogenetic split is also supported by the analysis of the 16S rRNA gene with both maximum likelihood (Figure 6) and neighbour-joining method (Figure S1). The 16S rRNA phylogeny assigned the M33 and M8 strains to I-A and I-B ribosomal subgroups, respectively [53].
The AY phytoplasma divergence was presented by several previously performed sequence analyses [53,54,55]. To verify if these results are also reflected at the single-gene level; a sequence identity was calculated based on the nucleotide sequences of the marker genes 16S rRNA, tufB, groEL, secA, and secY, and a chromosomal region containing the genes rplV and rpsC (Table 4). The pairwise values of the sequence identity matrix for the 16S rRNA gene showed no significant deviation from the threshold value of 98.65% of phytoplasmas [47] for the 16SrI-A group containing the strains M33 and AYWB. Pairwise comparison of the chromosomal region spanning from the gene rplV to rpsC resulted only in the AY strains DeVilla, M3, and Zhengzhou, with identity values that crossed the threshold of 97.50% and did not support the ANI assignment. Sequence identity analysis with secY appeared for strain M33 and AYWB identity values crossed the threshold of 95% and therefore supported the ANI cluster assignment. Similar support was achieved with the analysis of the genes secA, groEL, and tufB. Even though the threshold identity values were higher than the suggested thresholds, the values exceeded the threshold values, albeit to a minimal extent only. To confirm a significant overlap in species membership ‘Ca. P. rubi’ strain RS [56] was used as an outgroup in the phylogenetic analysis and exhibited significant divergence from the AY group clusters for all elicited markers compared to the cluster formed by strains M33 and AYWB (Table 4).

3.5. Pan-Genome Analyses

A comparison of the shared and unique features of the AY group was conducted via pan-genome analysis. In total 6717 CDS were used as queries for the ortholog prediction. Out of this total, 6505 CDS (96.8%) were assigned to 726 orthogroups with at least two members. The pairwise comparison of shared orthologs for each strain illustrates that ‘Ca. P. asteris’ M33 and AYWB shared the highest number of orthogroups (Table 5), whereas the carrot strains shared 401 orthogroups comprising a set of 503 proteins for M33 and 535 for M8, representing 84.5% and 72.2%, respectively. However, 326 single-copy orthologs were predicted to be shared by the asteris strains.
Therefore, the pan-genome analysis illustrated again the impact of paralog-associated information and supported the phylogenetic cluster assignment obtained from ANI, whole chromosome synteny, and single-gene analysis, and supports previous analysis [53]. It is notable that despite the number of shared orthologs the number of strain-specific features is low, with 212 CDS for the ten asteris chromosomes in total.
The strains M33 and M8 encode for 47 and 41 CDS not encoded by the other strains (Table S1, Figure S2). The deduced hypothetical proteins are characterised by a transmembrane domain (M33: 21/47; M8, 10/41). Six and five hypothetical proteins carry a signal peptide alone and represent putative secreted proteins (M33: M33023_01430, M33023_01470, M33023_03890, M33023_04050, M33023_04100, M33023_0677; M8: QN326_01250, QN326_04910, QN326_06420, QN326_06560, QN326_07200). According to these findings, the pan-genome analysis indicated that phytoplasmas of the AY group mainly differ in terms of features associated with pathogen–host interactions.
Multi-copy genes were predicted to investigate whether the genome size of phytoplasmas is associated with the occurrence of multi-copy genes (Table 6). The chromosome of strain M3 represents the smallest chromosome of the investigated AY group members with ~576 kb encoding 36 multicopy genes, whereas the Zhengzhou strain, with the largest chromosome of ~892 kb, possesses 385. These gene numbers represent 3.37% and 24.6% of the chromosome length, highlighting the impact. Strain M8 codes for 188 and M33 for 112 multi-copy genes, which puts them in the midfield with 14.1% and 9% of their chromosome length, respectively.
The impressive number of multi-copy genes is also characterised by a decreased G+C content for many asteris chromosomes, with the exception of strain MDGZ-01. It has the highest G+C content of the chromosomes examined and is separated from the other multi-copy genes in the other asteris chromosomes by an increased G+C content, which suggests horizontal gene transfer event(s) (Table 6). If considering the encoded function of the multi-copy genes, it was found that they are predominantly associated with the mobilome; for instance, genes known to be present on transposable elements are also described as potential mobile units (PMUs) or with genes originating from phage insertions as previously described for asteris [17] and phytoplasma genomes [57,58]. Instability is also supported by the detection of replicative forms [59] and recently for the origin and the coding of phyllogen genes in phytoplasmas [60].

3.6. Functional Reconstruction and Comparison

3.6.1. Key Metabolism and Membrane Transport

The analysed chromosomes encode the core metabolic pathways of phytoplasmas, i.e., the core module of the glycolysis (starting from glucose 6-phosphate), pyruvate oxidation to acetyl coenzyme A (CoA) and acetate, and phosphatidylethanolamine biosynthesis as part of glycerophospholipid metabolism (Figure 7). In addition, both genomes encode the malate-acetate pathway conserved in phytoplasmas [57]. However, the asteris strains lack the utilisation of lactate as reported for 16SrV phytoplasmas [61]. In the chromosomes of M8 and M33, 53 and 39 CDS were identified that code for subunits for ABC transporters for sugars, amino acids, peptides, polyamides, vitamins, and bivalent cations. In total, seven putative ABC transporters were identified in the examined genomes. The ABC-type multiple sugar transport system is the only sugar uptake system identified in phytoplasmas, while a phosphoenolpyruvate-dependent sugar phosphotransferase system is missing and leaves the problem of substrate phosphorylation unsolved [57]. Furthermore, three P-type ATPases were encoded. All P-type ATPases were assigned to the function of exporting cations such as sodium, potassium, and calcium, or other non-selective bivalent cations. Symporters involved in regulation, like the magnesium-cobalt exporter CorC and the malate-sodium symporter MaeN, were identified. Conserved genes coding for antiporters were represented by the MATE family efflux transporter and the large conductance mechanosensitive channel MscL.

3.6.2. Secretome and Characteristic Effector Proteins

M8 and M33 share with the other asteris strains the Sec-dependent secretion pathway encoded by the genes secA, secE, secY, yidC, ffh, and ftsY and representing the major functional secretion system identified within the investigated chromosomes of the AY group. Furthermore, components of the signal recognition particle (SRP) pathway were also identified. The SRP is involved in targeting ribosomes while translating and guiding them to the SecYEG pore complex of the Sec-dependent secretion pathway and targeting integral membrane proteins for co-translational integration into the membrane. The major SRP complex is formed by the SRP protein, a 4.5S RNA encoded by the gene ffs and the protein encoded by the gene ffh [62]. For a functional translocation system, the proteins YidC and FtsY are also needed, which were encoded in all AY genomes [63]. In total, 252 proteins for M8 and 231 proteins for M33 were predicted to be secreted by or integrated into the membrane. Predictions of the potential secretome with Phobius [43] comprised 43 proteins for the M8 strain and 33 for M33 (Table 7), which contained only a signal peptide domain (SP), 51 of which (67.11%) were shared by M8 and M33. Within the M8 chromosome, 201 proteins contained only transmembrane (TM) domain(s), whereas, for M33, 191 proteins were found, 297 (75.77%) of which were shared by strains M8 and M33. Eight proteins that had both (SP+TM) were predicted in M8 and seven were predicted in M33 (Figure 3). M8 and M33 shared eight out of a total of 15 proteins containing the SP and TM domains. Protein sequences that possess a signal peptide only were further investigated to identify the hitherto well-described and experimentally approved effector proteins TENGU, SAP05, SAP11, and SAP54. For all analysed chromosomes of the AY group, a TENGU gene coding for the effector was found, which indicates that tengu is conserved encoded in the AY group, which is in accordance with previous findings [53]. However, SAP05, SAP11, and SAP54 were not consistently encoded by the analysed complete ‘Ca. P. asteris’ chromosomes (Table S2). Notably, the genes coding for the effectors SAP11 and SAP54 were not identified for strain M33 whereas M8 encodes both genes. The lack of SAP11 stands in accordance with earlier findings from field studies supporting their absence in the 16SrI-B group [64]. Similar to the multi-copy gene analysis, these results also indicate that the numbers of genes coding for secreted or putative secreted proteins are higher within the larger chromosomes than in the smaller ones. This also supports the claim that interaction with the mobilome influences chromosome size, because secreted proteins such as the effectors from the SAP group are often present on transposable elements [59].

3.6.3. Immunodominant Membrane Proteins

The immunodominant membrane protein (Imp) and the antigenic membrane protein (Amp) are involved in physical interactions with the respective host organism of phytoplasmas. Imp is suggested to interact with the host plant actin filament, whereas Amp interacts with the actin filaments and the beta subunit of the ATPase within the insect environment [65,66,67,68]. Our analysis showed that the genes amp and imp were encoded in all AY group members. Konnerth and colleagues showed the genomic context for imp, which is directly adjacent to the gene pyrG coding for CTP synthase and the dnaG that codes for a protein involved in DNA replication [68]. However, in this work, we identified a different genomic context for imp within the ANI cluster of ‘Ca. P. asteris’ AYWB and M33, whereby only dnaG was adjacent. In contrast to the other ANI cluster comprising ‘Ca. P. asteris’ M8, OY-M, M3, RP166, De Villa, MDGZ-01, and Zhengzhou, this cluster has a pseudogene assigned to glucose-1-phosphatase flanked to imp. For amp, the same genomic context was found as reported by Konnerth and associates [68]. This region possesses the bordering genes nadE and groEL, as well as amp.

3.6.4. Adhesine P38

Another important group of cell surface factors comprises adhesins, which have hitherto been poorly studied in phytoplasmas. The adhesine P38 is suggested to interact with the insect host, as well as weakly with the plant host. Adhesin P38 was first described within the genome of ‘Ca. P. asteris’ OY-M [69]. A gene encoding adhesine P38 was found in all chromosomes of the analysed strains and showed a conserved genomic context. Flanking genes pyk and pepV coded for pyruvate kinase and the dipeptidase PepV.

3.6.5. Bax-Inhibitor 1

Bax-inhibitor 1 (BI-1) is a protein that has been characterised to reduce programmed cell death (PCD) [70] and it has also been suggested as being present in phytoplasma genomes [71]. However, how this PCD suppressor works in phytoplasmas is still not clear. Bax-inhibitor 1 was found in all analysed AY genomes. Within the chromosome of strain M33, two identical copies of the BI-1 gene were encoded, whereas the other strains showed only single genes. All genes coding for BI-1 revealed a conserved genomic context and were flanked by the genes tufB, coding for the elongation factor Tu, and rsmG encoding a 16S rRNA (guanine(527)-N(7))-methyltransferase. Annotations for rsmG of ‘Ca. P. asteris’ OY-M, AYWB, and M3 are biased compared to other AY group members and are deposited as gidB.

3.6.6. Superoxide Dismutase

Another notable feature is a gene-encoding superoxide dismutase (SOD), an important enzyme thought to be involved in protecting phytoplasmas against the plant defence response in the form of an oxidative burst created by enhanced reactive oxygen species (ROS) production [72,73]. Within the analysed AY genomes, a manganese/iron-dependent SOD was identified with a conserved genomic context, flanked by the genes pdhA coding for the pyruvate dehydrogenase E1 alpha subunit and nusB, which encodes a protein involved in transcription termination and antitermination.

3.7. Extrachromosomal Elements

Two plasmids, named pM8-6959 and pM33-16, were identified (Figure 8). The plasmids showed similar sequence lengths, with 5617 bp for pM8-6959 and 5045 bp for pM33-16. Eight ORFs were predicted for pM8-6959, whereas for pM33-16 six ORFs were identified. The exact positions of the OriC were not identified. The encoded rep genes were set at position 1 of the plasmids. In addition to, the replication-associated genes rep and ssb, a gene for a putative transmembrane domain-containing protein (QN326_00080) were found on pM8-6959, whereas for pM33-16, no encoded secreted proteins were identified. Moreover, genes that were assigned as putative copy number control proteins were found on both plasmids (QN326_00050; M33023_00050; M33023_00060). All remaining ORFs were assigned to hypothetical proteins that are not characterised by their function. Both plasmids were compared with each other and with all AY genomes at nucleotide and amino acid levels to identify potential interactions between the genomes in the form of horizontal gene transfer. The analysis showed no significant match at either the nucleotide or amino acid level. Results from the comparisons of the plasmids pM8-6959 and pM33-16 with the plasmid sequences from the other analysed complete genomes of ‘Ca. P. asteris’ showed that the rep and ssb gene sequences of pM8-6959 had the highest similarities with genes encoded on the plasmids of the onion-yellows (OY) group, whereas the sequences of the genes rep and ssb of pM33-16 shared the highest consensus with the plasmids pAWYB-I to IV of ‘Ca. P. asteris’ AYWB. This supports the assignment of M33 to 16SrIA and M8 to the 16Sr-IB group also at the extrachromosomal level.

4. Discussion

We added two complete genomes to the provisional taxon ‘Ca. P. asteris’. The strains ‘Ca. P. asteris’ M33 and M8—originating from the same field and host variety—were assigned to different phylogenetic clusters. The clusters reflect the phytoplasma phylogeny with the clearly diverged 16SrI-A and I-B subclades [74]. In view of this, an examination of the taxon status was proposed [53]. ANI cluster formation was also supported by single-gene phylogeny and whole-chromosome synteny analysis. An ongoing debate on the taxonomic revision of the genus ‘Candidatus Phytoplasma’ has recently started to consider ANI identity values, which again underlines the importance of this method for the taxonomic classification of newly discovered phytoplasma species. Moreover, the ANI method has been suggested as shifting the gold standard for the classification of prokaryotic species, especially in the era of big data, with an extreme increase in complete prokaryotic genome sequences [47,75,76,77,78]. In contrast to the three strains MDGZ-01, which infect mulberry (Morus alba), Zhengzhou [21], which is associated with paulownia (Paulownia fortunei), and QS2022 [18], which has been reported in lettuce (Lactuca sativa) in China, the M8 and M33 strains originated from the same field and host plant. These results are consistent with the occurrence of different 16SrI subgroups on individual fields during the season. Furthermore, our research confirms the disparities in the coded effectors of the SAP group within the 16SrI groups, which were illustrated by Clements and colleagues [64]. The phylogenetic analysis did not provide information on a shared geographical origin regarding the introduction into Germany (Figure 4). M8 shows a close relationship with ‘Ca. P. asteris’ RP166, which causes rapeseed phyllody in Poland; this may hint at the geographical location. The genetic distance of M8 and M33 might be linked to different vector populations or species. Cicadas of the genus Macrosteles are potential vectors. Another study from Germany showed that Macrosteles sexnotaus, M. laevis and M. cristatus caught in carrot fields infected with ‘Ca. Phytoplasma asteris’ also carried the pathogen, indicating that these cicadas are possible vectors [23].
Phytoplasmas have a narrow host range as far as their insect hosts are concerned, and since transovarial transmission is not a common case [79], the presence of an insect host is a crucial factor in their survival in the natural environment [80,81,82,83,84,85,86,87,88]. It cannot be excluded that this is also determined by the plasmids. Ishii and colleagues, for instance, provided evidence that the plasmid pOYNIM of ‘Ca. P. asteris’ OY-NIM had lost its orf3 and could therefore not be transmitted by its insect vector. The plasmid sequences of pM8-6959 and pM33-16 encoded no similar sequence to orf3. An exact functional characterisation of this ORF was not provided, but it was established that orf3 encodes a transmembrane protein representative on the cell surface [81,82,83]. Our results indicate that only the plasmid pM8-6959 possesses a gene coding for a protein with a length of 160 amino acids, which, according to Phobius prediction, is a transmembrane protein (QN326_00080), thus indicating that this gene could be involved in pathogen–host interaction (Figure 8). However, the exact role of this protein is still not clear, and whether the two strains ‘Ca. P. asteris’ M8 and M33 share one insect vector remains unclear. It has also been reported that phytoplasma transmission is possible through seed material from infected plants, which has been shown for carrots, and other important crop plants such as corn, (Zea mays), tomatoes (Solanum lycopersicum), winter oilseed rape (Brassica napus), and limes (Citrus aurantifolia) [4,89,90,91]. Mixing of different seed lots in the field of origin of M8 and M33 cannot be excluded but is unlikely to be the reason for the occurrence of the two different subgroups in one field since seed transmission is a rare case. Seed exports may represent a crucial issue considering the geographical distribution of these phytoplasmas, due to the fact that phytoplasmas in propagating material are not considered a risk either in the quarantine protocols of plant protection or by seed producers [79]. It is notable that the sampling field for M8 and M33 had poplar trees in their surroundings. Poplar trees have been reported as natural host plants for phytoplasmas in Europe, including black poplar (Populus nigra ‘Italica’ and P. canadensis) in Bulgaria, the Netherlands, Croatia, and Serbia; grey poplar (P. canescens) and white poplar (P. alba) from France and trembling poplar (P. tremula) in Germany. For France, Germany, and Serbia, Populus witches’ broom (PopWB) disease has been assigned to phytoplasmas of the aster yellows group [24,92,93,94,95,96,97,98]. Moreover, this was also demonstrated for some weed plants such as wild carrot (Daucus carota subsp. carota), hemlock (Conium maculatum), coast tarweed (Madia sativa), field bindweed (Convolvulus arvensis), and field madder (Sherardia arvensis). Therefore, screening for ‘Ca. P. asteris’ M8 and M33 and potential vectors associated with the poplar trees and weed plants in the tested cropping area should be considered.

5. Conclusions

In this study, we added two complete genome sequences to the provisional taxon ‘Ca. P. asteris’ sharing the host carrot. We confirmed the phylogenetic differentiation of the 16Sr I-A and I-B subclades in this taxon at the whole-genome level. Moreover, it was confirmed that the basic repertoire of genes coding for proteins with metabolic functions is highly conserved. Genomic plasticity regarding chromosome size is therefore not associated with extended metabolic functions but rather with duplication events and mobilome interactions. Pan-genome analysis of the AY group established that the unique features contributed to phytoplasma effector variability.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms12051016/s1. Table S1: CDS only encoded in ‘Ca. P. asteris’ strain M8 and M33, respectively, Table S2: Coding of the experimentally approved effector proteins within the asteris group, Figure S1. Phylogenetic tree of the ‘Ca. P. asteris’ strains were constructed using the neighbour-joining method, using 16S rDNA sequences of the employing ‘Ca. P. rubi’ strain RS as the outgroup. Numbers on the branches are bootstrap values obtained for 1000 replicates (only values above 70% are shown), Figure S2. Number of shared and unique CDS within the asteris group.

Author Contributions

Conceptualisation M.K.; methodology R.T., B.H., A.-M.I. and M.K.; software, R.T., A.-M.I., B.D. and M.K.; validation, A.-M.I., R.T. and M.K.; formal analysis, R.T., A.-M.I., B.H. and B.D.; investigation, R.T., A.-M.I., B.H. and B.D.; data curation, R.T. and B.D.; writing—original draft preparation, R.T. and M.K.; writing—review and editing, R.T., A.-M.I., B.H., B.D. and M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Annotation for ‘Ca. P. asteris’ M8 and M33 were deposited in the GenBank database under accession numbers CP128397.1 and CP128414.1.

Acknowledgments

We want to thank Daniela Rivatelli for technical support and Mathias Reindl for the organisation of the provided sample material from Langenau near the Schlammenbach (Baden-Württemberg, Germany).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Doi, Y.; Teranaka, M.; Yora, K.; Asuyama, H. Mycoplasma- or PLT Group-like Microorganisms Found in the Phloem Elements of Plants Infected with Mulberry Dwarf, Potato Witches’ Broom, Aster Yellows, or Paulownia Witches’ Broom. Jpn. J. Phytopathol. 1967, 33, 259–266. [Google Scholar] [CrossRef]
  2. Lee, I.M.; Davis, R.E.; Gundersen-Rindal, D.E. Phytoplasma: Phytopathogenic mollicutes. Annl. Rev. Microbiol. 2000, 54, 221–255. [Google Scholar] [CrossRef] [PubMed]
  3. Davis, M.R.; Raid, R.N. Compendium of Umbelliferous Crop Diseases, 1st ed.; American Phytopathological Society: St. Paul, MN, USA, 2002; pp. 58–59. [Google Scholar]
  4. Carminati, G.; Satta, E.; Paltrinieri, S.; Bertaccini, A. Simultaneous evaluation of ‘Candidatus Phytoplasma’ and ‘Candidatus Liberibacter solanacearum’ seed transmission in carrot. Phytopathog. Mollicutes 2019, 9, 141–142. [Google Scholar] [CrossRef]
  5. Randa Zelyüt, F.; Ertunç, F.; Şenal, D. The association of 16SrVI and 16SrI phytoplasma groups with carrot seeds and weeds in Ankara and Konya provinces in Turkey. Plant Prot. Bull. 2022, 62, 24–33. [Google Scholar] [CrossRef]
  6. Strauss, E. Phytoplasma research begins to bloom. Science 2009, 325, 388–390. [Google Scholar] [CrossRef] [PubMed]
  7. Font, I.; Abad, P.; Albinana, M.; Espino, A.I.; Dally, E.L.; Davis, R.E.; Jorda, C. Amarilleos y enrojecimientos en zanahoria: Una enfermedad a diagnóstico. Boletín De Sanid. Veg. Plagas 1999, 25, 405–415. [Google Scholar]
  8. Duduk, B.; Perić, P.; Marčić, D.; Drobnjaković, T.; Picciau, L.; Alma, A.; Bertaccini, A. Phytoplasmas in carrots: Disease and potential vectors in Serbia. Bull. Insectol. 2008, 61, 327–331. [Google Scholar]
  9. Munyaneza, J.E.; Lemmetty, A.; Nissinen, A.I.; Sengoda, V.G.; Fisher, T.W. Molecular detection of aster yellows phytoplasma and ‘Candidatus Liberibacter solanacearum’ in carrots affected by the psyllid Trioza apicalis (Hemiptera: Triozidae) in Finland. J. Plant Pathol. 2011, 93, 697–700. [Google Scholar]
  10. Nisbet, C.; Ross, S.; Monger, W.A.; Highet, F.; Jeffries, C. First report of ‘Candidatus Phytoplasma asteris’ in commercial carrots in the United Kingdom. New Dis. Rep. 2014, 30, 16. [Google Scholar] [CrossRef]
  11. Jonghe, K.D.; Roo, I.D.; Goedefroit, T. A survey in carrot reveals a widespread aster yellows infection, and a first case of ‘Candidatus Liberibacter solanacearum’ in Belgium. Phytopathog. Mollicutes 2019, 9, 139–140. [Google Scholar] [CrossRef]
  12. Zwolińska, A.; Borodynko-Filas, N. Intra and extragenomic variation between 16S rRNA genes found in 16SrI-B- related phytopathogenic phytoplasma strains. Ann. Appl. Biol. 2021, 179, 368–381. [Google Scholar] [CrossRef]
  13. Oshima, K.; Kakizawa, S.; Nishigawa, H.; Jung, H.-Y.; Wei, W.; Suzuki, S.; Arashida, R.; Nakata, D.; Miyata, S.; Ugaki, M.; et al. Reductive evolution suggested from the complete genome sequence of a plant-pathogenic phytoplasma. Nat. Genet. 2004, 36, 27–29. [Google Scholar] [CrossRef]
  14. Bai, X.; Correa, V.R.; Toruño, T.Y.; Ammar, E.-D.; Kamoun, S.; Hogenhout, S.A. AY-WB phytoplasma secretes a protein that targets plant cell nuclei. Mol. Plant Microbe Interact. 2009, 22, 18–30. [Google Scholar] [CrossRef]
  15. Wang, N.; Yang, H.; Yin, Z.; Liu, W.; Sun, L.; Wu, Y. Phytoplasma effector SWP1 induces witches’ broom symptom by destabilizing the TCP transcription factor BRANCHED1. Mol. Plant Pathol. 2018, 19, 2623–2634. [Google Scholar] [CrossRef]
  16. Huang, W.; MacLean, A.M.; Sugio, A.; Maqbool, A.; Busscher, M.; Cho, S.-T.; Kamoun, S.; Kuo, C.-H.; Immink, R.G.H.; Hogenhout, S.A. Parasitic modulation of host development by ubiquitin-independent protein degradation. Cell 2021, 184, 5201–5214.e12. [Google Scholar] [CrossRef]
  17. Bai, X.; Zhang, J.; Ewing, A.; Miller, S.A.; Jancso Radek, A.; Shevchenko, D.V.; Tsukerman, K.; Walunas, T.; Lapidus, A.; Campbell, J.W.; et al. Living with genome instability: The adaptation of phytoplasmas to diverse environments of their insect and plant hosts. J. Bacteriol. 2006, 188, 3682–3696. [Google Scholar] [CrossRef]
  18. Yan, X.-H.; Lin, J.; Liu, Y.; Huang, P.; Liu, J.; Hu, Q.; Li, Y.; Pei, S.-C.; Huang, W.; Kuo, C.-H. Complete Genome Sequence of ‘Candidatus Phytoplasma asteris’ QS2022, a Plant Pathogen Associated with Lettuce Chlorotic Leaf Rot Disease in China. Microbiol. Resour. Announc. 2023, 12, e0030623. [Google Scholar] [CrossRef]
  19. Orlovskis, Z.; Canale, M.C.; Haryono, M.; Lopes, J.R.S.; Kuo, C.-H.; Hogenhout, S.A. A few sequence polymorphisms among isolates of Maize bushy stunt phytoplasma associate with organ proliferation symptoms of infected maize plants. Ann. Bot. 2017, 119, 869–884. [Google Scholar] [CrossRef]
  20. Cho, S.-T.; Zwolińska, A.; Huang, W.; Wouters, R.H.M.; Mugford, S.T.; Hogenhout, S.A.; Kuo, C.-H. Complete Genome Sequence of ‘Candidatus Phytoplasma asteris’ RP166, a Plant Pathogen Associated with Rapeseed Phyllody Disease in Poland. Microbiol. Resour. Announc. 2020, 9, e00760-20. [Google Scholar] [CrossRef]
  21. Cao, Y.; Sun, G.; Zhai, X.; Xu, P.; Ma, L.; Deng, M.; Zhao, Z.; Yang, H.; Dong, Y.; Shang, Z.; et al. Genomic insights into the fast growth of paulownias and the formation of Paulownia witches’ broom. Mol. Plant 2021, 14, 1668–1682. [Google Scholar] [CrossRef]
  22. Luo, L.; Zhang, X.; Meng, F.; Wang, Y.; Zhou, Y.; Liu, J. Draft Genome Sequences Resources of Mulberry Dwarf Phytoplasma Strain MDGZ-01 Associated with Mulberry Yellow Dwarf (MYD) Diseases. Plant Dis. 2022, 106, 2239–2242. [Google Scholar] [CrossRef]
  23. Schneller, H.; Rißler, D.; Zgraja, G.; Zunker, M.; Zimmermann, O.; Kost, W.; Lasch, E.; Schrameyer, K. Erster Nachweis von ‚Aster-Yellows-Disease’ an Möhren (Phytoplasmen bedingte Möhrenröte) und in der Kleinzikade Macrosteles sexnotatus (Fallén 1806) in Deutschland –Monitoring und Diagnose. J. Kulturpflanz. 2016, 68, 281–294. [Google Scholar]
  24. Mitrović, J.; Paltrinieri, S.; Contaldo, N.; Bertaccini, A.; Duduk, B. Occurrence of two ‘Candidatus Phytoplasma asteris’—Related phytoplasmas in poplar trees in Serbia. Bull. Insectol. 2011, 64, 57–58. [Google Scholar]
  25. Murray, M.G.; Thompson, W.F. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980, 8, 4321–4325. [Google Scholar] [CrossRef]
  26. Deng, S.; Hiruki, C. Amplification of 16S rRNA genes from culturable and nonculturable Mollicutes. J. Microbiol. Meth. 1991, 14, 53–61. [Google Scholar] [CrossRef]
  27. Schneider, B.; Seemueller, E.; Smart, C.D.; Kirkpatrick, B.C. Phylogenetic Classification of Plant Pathogenic Mycoplasma-like Organisms or Phytoplasmas, 1st ed.; Academic Press: London, UK, 1995; pp. 369–380. [Google Scholar]
  28. Bentley, D.R.; Balasubramanian, S.; Swerdlow, H.P.; Smith, G.P.; Milton, J.; Brown, C.G.; Hall, K.P.; Evers, D.J.; Barnes, C.L.; Bignell, H.R.; et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature 2008, 456, 53–59. [Google Scholar] [CrossRef] [PubMed]
  29. Eid, J.; Fehr, A.; Gray, J.; Luong, K.; Lyle, J.; Otto, G.; Peluso, P.; Rank, D.; Baybayan, P.; Bettman, B.; et al. Real-time DNA sequencing from single polymerase molecules. Science 2009, 323, 133–138. [Google Scholar] [CrossRef]
  30. Koren, S.; Walenz, B.P.; Berlin, K.; Miller, J.R.; Bergman, N.H.; Phillippy, A.M. Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017, 27, 722–736. [Google Scholar] [CrossRef] [PubMed]
  31. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  32. Suzek, B.E.; Huang, H.; McGarvey, P.; Mazumder, R.; Wu, C.H. UniRef: Comprehensive and non-redundant UniProt reference clusters. Bioinformatics 2007, 23, 1282–1288. [Google Scholar] [CrossRef]
  33. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  34. Huson, D.H.; Auch, A.F.; Qi, J.; Schuster, S.C. MEGAN analysis of metagenomic data. Genome Res. 2007, 17, 377–386. [Google Scholar] [CrossRef]
  35. Carver, T.; Harris, S.R.; Berriman, M.; Parkhill, J.; McQuillan, J.A. Artemis: An integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics 2012, 28, 464–469. [Google Scholar] [CrossRef]
  36. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef]
  37. Aziz, R.K.; Bartels, D.; Best, A.A.; DeJongh, M.; Disz, T.; Edwards, R.A.; Formsma, K.; Gerdes, S.; Glass, E.M.; Kubal, M.; et al. The RAST Server: Rapid annotations using subsystems technology. BMC Genom. 2008, 9, 75. [Google Scholar] [CrossRef]
  38. Kanehisa, M.; Sato, Y.; Morishima, K. BlastKOALA and GhostKOALA: KEGG Tools for Functional Characterization of Genome and Metagenome Sequences. J. Mol. Biol. 2016, 428, 726–731. [Google Scholar] [CrossRef]
  39. Caspi, R.; Altman, T.; Billington, R.; Dreher, K.; Foerster, H.; Fulcher, C.A.; Holland, T.A.; Keseler, I.M.; Kothari, A.; Kubo, A.; et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res. 2014, 42, 459–471. [Google Scholar] [CrossRef]
  40. Paysan-Lafosse, T.; Blum, M.; Chuguransky, S.; Grego, T.; Pinto, B.L.; Salazar, G.A.; Bileschi, M.L.; Bork, P.; Bridge, A.; Colwell, L.; et al. InterPro in 2022. Nucleic Acids Res. 2023, 51, 418–427. [Google Scholar] [CrossRef]
  41. Lowe, T.M.; Chan, P.P. tRNAscan-SE On-line: Integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016, 44, 54–57. [Google Scholar] [CrossRef]
  42. Lagesen, K.; Hallin, P.; Rødland, E.A.; Staerfeldt, H.-H.; Rognes, T.; Ussery, D.W. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35, 3100–3108. [Google Scholar] [CrossRef]
  43. Käll, L.; Krogh, A.; Sonnhammer, E.L.L. Advantages of combined transmembrane topology and signal peptide prediction-the Phobius web server. Nucleic Acids Res. 2007, 35, 429–432. [Google Scholar] [CrossRef]
  44. Goris, J.; Konstantinidis, K.T.; Klappenbach, J.A.; Coenye, T.; Vandamme, P.; Tiedje, J.M. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int. J. Syst. Evol. Microbiol. 2007, 57, 81–91. [Google Scholar] [CrossRef]
  45. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar]
  46. Darling, A.C.E.; Mau, B.; Blattner, F.R.; Perna, N.T. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14, 1394–1403. [Google Scholar] [CrossRef]
  47. Bertaccini, A.; Arocha-Rosete, Y.; Contaldo, N.; Duduk, B.; Fiore, N.; Montano, H.G.; Kube, M.; Kuo, C.-H.; Martini, M.; Oshima, K.; et al. Revision of the ‘Candidatus Phytoplasma’ species description guidelines. Int. J. Syst. Evol. Microbiol. 2022, 72, 5353. [Google Scholar] [CrossRef]
  48. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  49. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  50. Fischer, R. On the mathematical foundations of theoretical statistics. Philos. Trans. R. Soc. A 1922, 222, 309–368. [Google Scholar]
  51. Emms, D.M.; Kelly, S. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biol. 2019, 20, 238. [Google Scholar] [CrossRef] [PubMed]
  52. Zübert, C.; Ilic, A.-M.; Duduk, B.; Kube, M. The Genome Reduction Excludes the Ribosomal Rescue System in Acholeplasmataceae. Microb. Physiol. 2022, 32, 45–56. [Google Scholar] [CrossRef]
  53. Cho, S.-T.; Kung, H.-J.; Huang, W.; Hogenhout, S.A.; Kuo, C.-H. Species Boundaries and Molecular Markers for the Classification of 16SrI Phytoplasmas Inferred by Genome Analysis. Front. Microbiol. 2020, 11, 1531. [Google Scholar] [CrossRef]
  54. Duduk, B.; Calari, A.; Paltrinieri, S.; Duduk, N.; Bertaccini, A. Multigene analysis for differentiation of aster yellows phytoplasmas infecting carrots in Serbia. Ann. Appl. Biol. 2009, 154, 219–229. [Google Scholar] [CrossRef]
  55. Mitrović, J.; Kakizawa, S.; Duduk, B.; Oshima, K.; Namba, S.; Bertaccini, A. The groEL gene as an additional marker for finer differentiation of ‘Candidatus Phytoplasma asteris’-related strains. Ann. Appl. Biol. 2011, 159, 41–48. [Google Scholar] [CrossRef]
  56. Duckeck, D.; Zübert, C.; Böhm, J.W.; Carminati, G.; Schneider, B.; Kube, M. Complete Genome of ‘Candidatus Phytoplasma rubi’ RS, a Phytopathogenic Bacterium Associated with Rubus Stunt Disease. Microbiol. Resour. Announc. 2023, 12, e0130322. [Google Scholar] [CrossRef]
  57. Kube, M.; Mitrovic, J.; Duduk, B.; Rabus, R.; Seemüller, E. Current view on phytoplasma genomes and encoded metabolism. Sci. World J. 2012, 2012, 185942. [Google Scholar] [CrossRef]
  58. Ku, C.; Lo, W.-S.; Kuo, C.-H. Horizontal transfer of potential mobile units in phytoplasmas. Mob. Genet. Elements 2013, 3, e26145. [Google Scholar] [CrossRef]
  59. Toruño, T.Y.; Musić, M.S.; Simi, S.; Nicolaisen, M.; Hogenhout, S.A. Phytoplasma PMU1 exists as linear chromosomal and circular extrachromosomal elements and has enhanced expression in insect vectors compared with plant hosts. Mol. Microbiol. 2010, 77, 1406–1415. [Google Scholar] [CrossRef]
  60. Tokuda, R.; Iwabuchi, N.; Kitazawa, Y.; Nijo, T.; Suzuki, M.; Maejima, K.; Oshima, K.; Namba, S.; Yamaji, Y. Potential mobile units drive the horizontal transfer of phytoplasma effector phyllogen genes. Front. Genet. 2023, 14, 1132432. [Google Scholar] [CrossRef]
  61. Böhm, J.W.; Duckeck, D.; Duduk, B.; Schneider, B.; Kube, M. Genome Comparison of ‘Candidatus Phytoplasma rubi’ with Genomes of Other 16SrV Phytoplasmas Highlights Special Group Features. Appl. Microbiol. 2023, 3, 1083–1100. [Google Scholar] [CrossRef]
  62. Koch, H.G.; Hengelage, T.; Neumann-Haefelin, C.; MacFarlane, J.; Hoffschulte, H.K.; Schimz, K.L.; Mechler, B.; Müller, M. In vitro studies with purified components reveal signal recognition particle (SRP) and SecA/SecB as constituents of two independent protein-targeting pathways of Escherichia coli. Mol. Biol. Cell 1999, 10, 2163–2173. [Google Scholar] [CrossRef]
  63. Koch, H.-G.; Moser, M.; Müller, M. Signal recognition particle-depencent protein targeting, universal to all kingdoms of life. In Reviews of Physiology, Biochemistry and Pharmacology, 1st ed.; Amara, S.G., Bamberg, E., Blaustein, M.P., Grunicke, H., Jahn, R., Lederer, W.J., Miyajima, A., Murer, H., Offermanns, S., Pfanner, N., et al., Eds.; Springer: Berlin/Heidelberg, Germany, 2003; pp. 55–94. [Google Scholar]
  64. Clements, J.; Bradford, B.Z.; Garcia, M.; Piper, S.; Huang, W.; Zwolinska, A.; Lamour, K.; Hogenhout, S.; Groves, R.L. ‘Candidatus Phytoplasma asteris’ subgroups display distinct disease progression dynamics during the carrot growing season. PLoS ONE 2021, 16, e0239956. [Google Scholar] [CrossRef]
  65. Barbara, D.J.; Morton, A.; Clark, M.F.; Davies, D.L. Immunodominant membrane proteins from two phytoplasmas in the aster yellows clade (chlorante aster yellows and clover phyllody) are highly divergent in the major hydrophilic region. Microbiology 2002, 148, 157–167. [Google Scholar] [CrossRef]
  66. Kakizawa, S.; Oshima, K.; Nishigawa, H.; Jung, H.-Y.; Wei, W.; Suzuki, S.; Tanaka, M.; Miyata, S.; Ugaki, M.; Namba, S. Secretion of immunodominant membrane protein from onion yellows phytoplasma through the Sec protein-translocation system in Escherichia coli. Microbiology 2004, 150, 135–142. [Google Scholar] [CrossRef]
  67. Kakizawa, S.; Oshima, K.; Ishii, Y.; Hoshi, A.; Maejima, K.; Jung, H.-Y.; Yamaji, Y.; Namba, S. Cloning of immunodominant membrane protein genes of phytoplasmas and their in planta expression. FEMS Microbiol. Lett. 2009, 293, 92–101. [Google Scholar] [CrossRef] [PubMed]
  68. Konnerth, A.; Krczal, G.; Boonrod, K. Immunodominant membrane proteins of phytoplasmas. Microbiology 2016, 162, 1267–1273. [Google Scholar] [CrossRef]
  69. Neriya, Y.; Maejima, K.; Nijo, T.; Tomomitsu, T.; Yusa, A.; Himeno, M.; Netsu, O.; Hamamoto, H.; Oshima, K.; Namba, S. Onion yellow phytoplasma P38 protein plays a role in adhesion to the hosts. FEMS Microbiol. Lett. 2014, 361, 115–122. [Google Scholar] [CrossRef]
  70. Hückelhoven, R. BAX Inhibitor-1, an ancient cell death suppressor in animals and plants with prokaryotic relatives. Apoptosis 2004, 9, 299–307. [Google Scholar] [CrossRef]
  71. Quaglino, F.; Kube, M.; Jawhari, M.; Abou-Jawdah, Y.; Siewert, C.; Choueiri, E.; Sobh, H.; Casati, P.; Tedeschi, R.; Lova, M.M.; et al. ‘Candidatus Phytoplasma phoenicium’ associated with almond witches’-broom disease: From draft genome to genetic diversity among strain populations. BMC Microbiol. 2015, 15, 148. [Google Scholar] [CrossRef]
  72. Wojtaszek, P. Oxidative burst: An early plant response to pathogen infection. Biochem. J. 1997, 322 Pt 3, 681–692. [Google Scholar] [CrossRef]
  73. Miura, C.; Sugawara, K.; Neriya, Y.; Minato, N.; Keima, T.; Himeno, M.; Maejima, K.; Komatsu, K.; Yamaji, Y.; Oshima, K.; et al. Functional characterization and gene expression profiling of superoxide dismutase from plant pathogenic phytoplasma. Gene 2012, 510, 107–112. [Google Scholar] [CrossRef]
  74. Lee, I.-M.; Gundersen-Rindal, D.E.; Davis, R.E.; Bottner, K.D.; Marcone, C.; Seemüller, E. ‘Candidatus Phytoplasma asteris’, a novel phytoplasma taxon associated with aster yellows and related diseases. Int. J. Syst. Evol. Microbiol. 2004, 54, 1037–1048. [Google Scholar] [CrossRef]
  75. Richter, M.; Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 2009, 106, 19126–19131. [Google Scholar] [CrossRef]
  76. Jain, C.; Rodriguez-R, L.M.; Phillippy, A.M.; Konstantinidis, K.T.; Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 2018, 9, 5114. [Google Scholar] [CrossRef]
  77. Hugenholtz, P.; Chuvochina, M.; Oren, A.; Parks, D.H.; Soo, R.M. Prokaryotic taxonomy and nomenclature in the age of big sequence data. ISME J. 2021, 15, 1879–1892. [Google Scholar] [CrossRef]
  78. Wei, W.; Zhao, Y. Phytoplasma Taxonomy: Nomenclature, Classification, and Identification. Biology 2022, 11, 1119. [Google Scholar] [CrossRef]
  79. Bertaccini, A.; Weintraub, P.G.; Rao, G.P.; Mori, N.; Bertaccini, A. Phytoplasmas: Plant Pathogenic Bacteria—II: Transmission and Management of Phytoplasma—Associated Diseases, 1st ed.; Springer Nature: Singapore, 2019; pp. 5–9. [Google Scholar]
  80. Nishimura, N.; Nakajima, S.; Sawayanagi, T.; Namba, S.; Shiomi, T.; Matsuda, I.; Tsuchizaki, T. Transmission of Cryptotaenia japonica Witches’ Broom and Onion Yellows Phytoplasmas by Hishimonus sellatus Uhler. Jpn. J. Phytopathol. 1998, 64, 474–477. [Google Scholar] [CrossRef]
  81. Nishigawa, H.; Oshima, K.; Kakizawa, S.; Jung, H.-Y.; Kuboyama, T.; Miyata, S.-I.; Ugaki, M.; Namba, S. A plasmid from a non-insect-transmissible line of a phytoplasma lacks two open reading frames that exist in the plasmid from the wild-type line. Gene 2002, 298, 195–201. [Google Scholar] [CrossRef]
  82. Nishigawa, H.; Oshima, K.; Kakizawa, S.; Jung, H.-Y.; Kuboyama, T.; Miyata, S.-I.; Ugaki, M.; Namba, S. Evidence of intermolecular recombination between extrachromosomal DNAs in phytoplasma: A trigger for the biological diversity of phytoplasma? Microbiology 2002, 148, 1389–1396. [Google Scholar] [CrossRef]
  83. Nishigawa, H.; Oshima, K.; Miyata, S.-I.; Ugaki, M.; Namba, S. Complete set of extrachromosomal DNAs from three pathogenic lines of onion yellows phytoplasma and use of PCR to differentiate each line. J. Gen. Plant Pathol. 2003, 69, 194–198. [Google Scholar]
  84. Nishimura, N.; Nakajima, S.; Kawakita, H.; Sato, M.; Namba, S.; Fujisawa, I.; Tsuchizaki, T. Transmission of Cryptotaenia japonica Witches’ Broom and Onion Yellows by Hishimonoides sellatiformis. Jpn. J. Phytopathol. 2004, 70, 22–25. [Google Scholar] [CrossRef]
  85. Wei, W.; Kakizawa, S.; Suzuki, S.; Jung, H.-Y.; Nishigawa, H.; Miyata, S.; Oshima, K.; Ugaki, M.; Hibi, T.; Namba, S. In planta dynamic analysis of onion yellows phytoplasma using localized inoculation by insect transmission. Phytopathology 2004, 94, 244–250. [Google Scholar] [CrossRef] [PubMed]
  86. Weintraub, P.G.; Beanland, L. Insect vectors of phytoplasmas. Annu. Rev. Entomol. 2006, 51, 91–111. [Google Scholar] [CrossRef]
  87. Ishii, Y.; Oshima, K.; Kakizawa, S.; Hoshi, A.; Maejima, K.; Kagiwada, S.; Yamaji, Y.; Namba, S. Process of reductive evolution during 10 years in plasmids of a non-insect-transmissible phytoplasma. Gene 2009, 446, 51–57. [Google Scholar] [CrossRef]
  88. Bertaccini, A.; Oshima, K.; Kakizawa, S.; Duduk, B.; Namba, S. Dissecting the Multifaceted Mechanisms That Drive Leafhopper Host–Phytoplasma Specificity. In Vector-Mediated Transmission of Plant Pathogens, 1st ed.; Judith, K.B., Ed.; American Phytopathological Society: St. Paul, MN, USA, 2016; pp. 21–28. [Google Scholar]
  89. Calari, A.; Paltrinieri, S.; Contaldo, N.; Sakalieva, D.; Mori, N.; Duduk, B.; Bertaccini, A. Molecular evidence of phytoplasmas in winter oilseed rape, tomato and corn seedlings. Bull. Insectol. 2011, 64, 151–158. [Google Scholar]
  90. Satta, E.; Ramirez, A.S.; Paltrinieri, S.; Contaldo, N.; Benito, P.; Poveda, J.B.; Bertaccini, A. Simultaneous detection of mixed ‘Candidatus Phytoplasma asteris’ and ‘Ca. Liberibacter solanacearum’ infection in carrot. Phytopathol. Mediterr. 2016, 55, 401–409. [Google Scholar]
  91. Satta, E.; Carminati, G.; Bertaccini, A. Phytoplasma presence in carrot seedlings. Australas. Plant Dis. Notes 2020, 15, 11. [Google Scholar] [CrossRef]
  92. Atanasoff, D. Stammhexenbesen bei Ulmen und anderen Bäumen. Arch. Pflanzenschutz 1973, 9, 241–243. [Google Scholar] [CrossRef]
  93. Van der Meer, F.A. Witches’ brooms in poplars? Populier 1980, 17, 42–43. [Google Scholar]
  94. Sharma, A.K.; Cousin, M.-T. Mycoplasmalike Organisms (MLOs) Associated With the Witches’ Broom Disease of Poplar. J. Phytopathol. 1986, 117, 349–356. [Google Scholar] [CrossRef]
  95. Seemüller, E.; Lederer, W. MLO-Associated Decline of Alnus glutinosa, Populus tremula and Crataegus monogyna. J. Phytopathol. 1988, 121, 33–39. [Google Scholar] [CrossRef]
  96. Berges, R.; Cousin, M.-T.; Roux, J.; Mäurer, R.; Seemüller, E. Detection of phytoplasma infections in declining Populus nigra ‘Italica’ trees and molecular differentiation of the aster yellows phytoplasmas identified in various Populus species. Eur. J. For. Pathol. 1997, 27, 33–43. [Google Scholar] [CrossRef]
  97. Marcone, C.; Ragozzino, A.; Seemuller, E. Dodder transmission of alder yellows phytoplasma to the experimental host Catharanthus roseus (periwinkle). Eur. J. For. Pathol. 1997, 27, 347–350. [Google Scholar] [CrossRef]
  98. Šeruga, M.; Škorić, D.; Botti, S.; Paltrinieri, S.; Juretić, N.; Bertaccini, A. Molecular identification of a phytoplasma naturally infecting Populus nigra L. cv. Italica trees in Croatia. Plant Prot. Sci. 2002, 38, 28–30. [Google Scholar] [CrossRef]
Microorganisms 12 01016 g001
Figure 1. Carrots infected with ‘Ca. P. asteris’ strains M8 and M33 with leaves showing chlorosis, partial reddening, and terminal necrosis. In addition, infection with M33 shows a thin root (“rat tail” symptom).
Figure 1. Carrots infected with ‘Ca. P. asteris’ strains M8 and M33 with leaves showing chlorosis, partial reddening, and terminal necrosis. In addition, infection with M33 shows a thin root (“rat tail” symptom).
Microorganisms 12 01016 g001
Microorganisms 12 01016 g002
Figure 2. BUSCO results for the ‘Ca. P. asteris’ strains. Results are separated into (C) complete single copy orthologues comprising (S) single and (D) duplicated copy orthologues (blue), (F) fragmented orthologues (yellow), and (M) missing orthologues (red).
Figure 2. BUSCO results for the ‘Ca. P. asteris’ strains. Results are separated into (C) complete single copy orthologues comprising (S) single and (D) duplicated copy orthologues (blue), (F) fragmented orthologues (yellow), and (M) missing orthologues (red).
Microorganisms 12 01016 g002
Microorganisms 12 01016 g003
Figure 3. Genomic organisation of the circular chromosomes of ‘Ca. P. asteris’ strains M8 (A) and M33 (B). Circular patterns (from outside): 1 (outer ring), scale in base pairs of the chromosome; 2 (cyan), predicted protein-coding sequences; 2 (grey), hypothetical proteins; 3 (multi-coloured), predicted orthologues shared with all investigated genomes of the asteris group (dark blue), unique protein-coding sequences(brown), predicted paralogs (red); 4 (multi-coloured), predicted membrane proteins: with only a signal peptide (orange), only transmembrane domains (yellow), and both signal peptide and transmembrane domains (mid red); 5 (multi-coloured), RNAs: Signal recognition particle RNA and RNase P RNA component class B (magenta), transfer RNAs (green), and rRNA operons (black),transfer-messenger RNAs (light blue); 6 G+C skew (olive and pink).
Figure 3. Genomic organisation of the circular chromosomes of ‘Ca. P. asteris’ strains M8 (A) and M33 (B). Circular patterns (from outside): 1 (outer ring), scale in base pairs of the chromosome; 2 (cyan), predicted protein-coding sequences; 2 (grey), hypothetical proteins; 3 (multi-coloured), predicted orthologues shared with all investigated genomes of the asteris group (dark blue), unique protein-coding sequences(brown), predicted paralogs (red); 4 (multi-coloured), predicted membrane proteins: with only a signal peptide (orange), only transmembrane domains (yellow), and both signal peptide and transmembrane domains (mid red); 5 (multi-coloured), RNAs: Signal recognition particle RNA and RNase P RNA component class B (magenta), transfer RNAs (green), and rRNA operons (black),transfer-messenger RNAs (light blue); 6 G+C skew (olive and pink).
Microorganisms 12 01016 g003
Microorganisms 12 01016 g004
Figure 4. Phylogenetic tree based on ANI comparison from the analysed chromosomes of ‘Ca. P. asteris’ strains using the neighbour-joining method. Branch lengths are measured via the number of substitutions per digit.
Figure 4. Phylogenetic tree based on ANI comparison from the analysed chromosomes of ‘Ca. P. asteris’ strains using the neighbour-joining method. Branch lengths are measured via the number of substitutions per digit.
Microorganisms 12 01016 g004
Microorganisms 12 01016 g005
Figure 5. Sequence synteny analysis of the ‘Ca. P. asteris’ strains using Mauve. The outer blue and red box contain selected representatives of the respective ANI clusters of ‘Ca. P. asteris’ M8 and M33 as well as their associated 16SrI subgroups. Inner blocks with identical colours show sequence synteny.
Figure 5. Sequence synteny analysis of the ‘Ca. P. asteris’ strains using Mauve. The outer blue and red box contain selected representatives of the respective ANI clusters of ‘Ca. P. asteris’ M8 and M33 as well as their associated 16SrI subgroups. Inner blocks with identical colours show sequence synteny.
Microorganisms 12 01016 g005
Microorganisms 12 01016 g006
Figure 6. Phylogenetic tree of the ‘Ca. P. asteris’ strains constructed using the maximum likelihood method, using 16S rDNA sequences of the employing ‘Ca. P. rubi’ strain RS as the outgroup. Numbers on the branches are bootstrap values obtained for 1000 replicates (only values above 70% are shown). Strains from this study and their corresponding 16SrI subgroup are highlighted in bold.
Figure 6. Phylogenetic tree of the ‘Ca. P. asteris’ strains constructed using the maximum likelihood method, using 16S rDNA sequences of the employing ‘Ca. P. rubi’ strain RS as the outgroup. Numbers on the branches are bootstrap values obtained for 1000 replicates (only values above 70% are shown). Strains from this study and their corresponding 16SrI subgroup are highlighted in bold.
Microorganisms 12 01016 g006
Microorganisms 12 01016 g007
Figure 7. Schematic overview of the complete metabolic pathways suggesting membrane transport, and membrane proteins involved in pathogen–host interaction with ‘Ca. P. asteris’. Curved arrows indicate ATP hydrolysis. The unclear mechanism of phosphorylation and substrate supply for glycolysis have not been clarified in detail and have been labelled with a question mark in consequence.
Figure 7. Schematic overview of the complete metabolic pathways suggesting membrane transport, and membrane proteins involved in pathogen–host interaction with ‘Ca. P. asteris’. Curved arrows indicate ATP hydrolysis. The unclear mechanism of phosphorylation and substrate supply for glycolysis have not been clarified in detail and have been labelled with a question mark in consequence.
Microorganisms 12 01016 g007
Microorganisms 12 01016 g008
Figure 8. Plasmid map of pM8-6959 (A) and pM33-16 (B). Colour code: rep coding for the replication protein (red); ssb coding for single-strand binding protein (blue); putative secreted proteins (green); genes coding for putative copy number controlling proteins (pink); (grey) CDS assigned as hypothetical proteins.
Figure 8. Plasmid map of pM8-6959 (A) and pM33-16 (B). Colour code: rep coding for the replication protein (red); ssb coding for single-strand binding protein (blue); putative secreted proteins (green); genes coding for putative copy number controlling proteins (pink); (grey) CDS assigned as hypothetical proteins.
Microorganisms 12 01016 g008
Table 1. Canu hybrid assembly statistics of the complete genomes of ‘Ca. P. asteris’ M8 and M33.
Table 1. Canu hybrid assembly statistics of the complete genomes of ‘Ca. P. asteris’ M8 and M33.
M8M33
Assembly input
Selected/mapped Illumina reads (%)223,496 (10.6)480,718 (22.2)
Average read length (nt)250250
SMRT reads total191,988105,166
Average read length (nt)86949275
Incorporated reads
No. of Illumina & SMRT reads used for chromosome reconstruction (sequencing coverage)10,11611,776
(115.86-fold) (176.15-fold)
No. of reads used for plasmid reconstruction (sequencing coverage)30745
(95.09-fold)(11.97-fold)
Genome organization (complete)One circular chromosome, one plasmid One circular chromosome, one plasmid
Table 2. Genomic benchmarks of the analysed complete chromosomes within the taxon ‘Ca. P. asteris’ according to annotation.
Table 2. Genomic benchmarks of the analysed complete chromosomes within the taxon ‘Ca. P. asteris’ according to annotation.
16Sr-GroupStrainCountry of OriginPlant HostLength (bp)G+C
Content (%)		CDS
(+ Pseudo Genes)		Coding Density (per kb)Acc. no. 1
16SrI-AM33GermanyCarrot657,32426.79595 (51)0.905CP128397.1
AYWBUSALettuce706,56926.89586 (62)0.829NC_007716.1
16SrI-BM8GermanyCarrot772,69127.87741 (77)0.958CP128414.1
OY-MJapanOnion853,09227.76752 (-)0.881AP006628.2
M3BrazilMaize576,11828.46485 (-)0.841CP015149.1
DeVillaSouth AfricaPeriwinkle600,11628.44545 (10)0.908CP035949.1
RP166PolandRapeseed829,54627.70753 (-)0.907CP055264.1
ZhengzhouChinaPaulownia891,64127.35906 (68)1.016NZ_CP066882.1
MDGZ-01ChinaMulberry622,35829.09535 (10)0.859NZ_CP085837.1
QS2022ChinaLettuce834,30327.57819 (-)0.981CP120448.1
1 Data from RefSeq were used with regard to continuous curation (indicated by underline) if available. Data from this study is highlighted bold.
Table 3. ANI matrix of the compared ‘Ca. P. asteris’ chromosome sequences.
Table 3. ANI matrix of the compared ‘Ca. P. asteris’ chromosome sequences.
StrainM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022
M33 99.4493.1992.9792.8892.8593.2992.9592.9392.91
AYWB99.44 93.2093.0092.8092.8093.1992.9192.9392.81
M893.1993.20 99.0298.5498.6499.9398.6698.9198.44
OY-M92.9793.0099.02 98.4498.4999.0198.4998.7099.12
M392.8892.8098.5498.44 99.1598.5698.2498.3798.23
DeVilla92.8592.8098.6498.4999.15 98.6398.2998.4098.31
RP16693.2993.1999.9399.0198.5698.63 98.6498.9198.45
Zhengzhou92.9592.9198.6698.4998.2498.2998.64 98.3598.12
MDGZ-0192.9392.9398.9198.7098.3798.4098.9198.35 98.51
QS202292.9192.8198.4499.1298.2398.3198.4598.1298.51
Identity values in percentages. Bold values show identities crossing the species affiliation threshold [47].
Table 4. Sequence identity matrices for the genes 16S rRNA, tufB, groEL, secA, secY, and the genomic region spanning between rplV to rpsC.
Table 4. Sequence identity matrices for the genes 16S rRNA, tufB, groEL, secA, secY, and the genomic region spanning between rplV to rpsC.
16SrRNAM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022RS
M33-99.499.398.799.298.89998.998.899.189.1
AYWB99.4-98.899.298.699.298.599.499.498.689.4
M899.398.8-99.199.49999.799.299.199.488.9
OY-M98.799.299.1-98.899.698.899.699.698.889.3
M399.298.699.498.8-9999.7999999.888.7
DeVilla98.899.29999.699-98.899.699.598.889.5
RP1669998.599.798.899.798.8-98.998.899.688.6
Zhengzhou98.999.499.299.69999.698.9-99.89989.3
MDGZ-0198.899.499.199.69999.598.899.8-98.989.2
QS202299.198.699.498.899.898.899.69998.9-88.6
RS89.189.488.989.388.789.588.689.389.288.6-
tufBM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022RS
M33-99.896.796.796.796.796.796.496.796.868
AYWB99.8-96.896.896.896.796.896.596.896.968
M896.796.8-10099.699.410099.610099.767.8
OY-M96.796.8100-99.699.410099.610099.767.8
M396.796.899.699.6-99.799.699.399.699.567.8
DeVilla96.796.799.499.499.7-99.49999.499.367.9
RP16696.796.810010099.699.4-99.610099.767.8
Zhengzhou96.496.599.699.699.39999.6-99.699.467.7
MDGZ-0196.796.810010099.699.410099.6-99.767.8
QS202296.896.999.799.799.599.399.799.499.7-67.8
RS686867.867.867.867.967.867.767.867.8-
groELM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022RS
M33-99.89797.196.896.99796.896.896.869.5
AYWB99.8-97.297.2979797.296.996.99769.5
M89797.2-99.999.899.710099.799.799.669.4
OY-M97.197.299.9-99.799.899.999.699.699.569.4
M396.89799.899.7-99.899.899.599.599.469.5
DeVilla96.99799.799.899.8-99.799.599.599.369.4
RP1669797.210099.999.899.7-99.799.799.669.4
Zhengzhou96.896.999.799.699.599.599.7-99.599.369.3
MDGZ-0196.896.999.799.699.599.599.799.5-99.369.6
QS202296.89799.699.599.499.399.699.399.3-69.4
RS69.569.569.469.469.569.469.469.369.669.4-
secAM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022RS
M33-99.995.295.295.195.195.295.195.395.267.7
AYWB99.9-95.295.395.195.195.295.295.395.367.8
M895.295.2-99.699.299.210099.499.699.367.7
OY-M95.295.399.6-99.399.399.699.599.699.467.8
M395.195.199.299.3-99.499.299.299.499.268.1
DeVilla95.195.199.299.399.4-99.299.299.499.267.9
RP16695.295.210099.699.299.2-99.499.699.367.7
Zhengzhou95.195.299.499.599.299.299.4-99.599.367.8
MDGZ-0195.395.399.699.699.499.499.699.5-99.467.8
QS202295.295.399.399.499.299.299.399.399.4-67.9
RS67.767.867.767.868.167.967.767.867.867.9-
secYM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022RS
M33-99.894.794.294.594.494.794.293.994.255
AYWB99.8-94.994.494.694.694.994.494.194.355.2
M894.794.9-99.499.399.310099.19999.355.5
OY-M94.294.499.4-99.199.199.498.998.799.755.4
M394.594.699.399.1-99.599.398.898.79955.5
DeVilla94.494.699.399.199.5-99.39998.79955.6
RP16694.794.910099.499.399.3-99.19999.355.5
Zhengzhou94.294.499.198.998.89999.1-98.598.855.5
MDGZ-0193.994.19998.798.798.79998.5-98.755.5
QS202294.294.399.399.7999999.398.898.7-55.3
RS5555.255.555.455.555.655.555.555.555.3-
rplV-rpsCM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022RS
M33-10097.797.597.197.397.79797.597.563.6
AYWB100-97.797.597.197.397.79797.597.563.6
M897.797.7-99.899.499.410099.399.899.863.2
OY-M97.597.599.8-99.299.299.899.299.699.663.2
M397.197.199.499.2-99.699.498.899.299.263.1
DeVilla97.397.399.499.299.6-99.498.899.299.263.2
RP16697.797.710099.899.499.4-99.399.899.863.2
Zhengzhou979799.399.298.898.899.3-99.299.263
MDGZ-0197.597.599.899.699.299.299.899.2-99.663.1
QS202297.597.599.899.699.299.299.899.299.6-63.2
RS63.663.663.263.263.163.263.26363.163.2-
Identity values in percentages. Bold values show identities crossing the species affiliation threshold [47].
Table 5. Pairwise comparison of shared ‘Ca. P. asteris’ orthogroups.
Table 5. Pairwise comparison of shared ‘Ca. P. asteris’ orthogroups.
StrainM33AYWBM8OY-MM3DeVillaRP166ZhengzhouMDGZ-01QS2022
M33-421401382370370398381365403
AYWB421-410389377373411388376404
M8401410-437418432597448420470
OY-M382389437-399411437441407498
M3370377418399-429414413400414
DeVilla370373432411429-421443413426
RP166398411597437414421-443411463
Zhengzhou381388448441413443443-420456
MDGZ-01365376420407400413411420-406
QS2022403404470498414426463456406-
Data from this study is highlighted in bold.
Table 6. Impact of multi-copy genes in complete asteris chromosomes.
Table 6. Impact of multi-copy genes in complete asteris chromosomes.
16Sr-GroupStrainAll GenesMulti-Copy Genes
No.No.% of All Σ Length in bpChromosome (%)
16SrI-AM3359511218.8259,4519.04
AYWB58615526.4585,75512.14
16SrI-BM874118825.37108,65114.06
OY-M75225133.38185,03721.69
M3485367.4219,4253.37
DeVilla5457814.3148,1748.03
RP16675326935.72177,86121.44
Zhengzhou90638542.49219,27024.59
MDGZ-015356612.3446,8007.52
QS202281929135.53172,76720.71
Data from this study is highlighted in bold.
Table 7. Differentiation of proteins with respect to coding or absence of a signal and/or one or more membrane domain(s).
Table 7. Differentiation of proteins with respect to coding or absence of a signal and/or one or more membrane domain(s).
StrainsCytosolic Signal PeptideSignal Peptide and Transmembrane DomaineTransmembrane Domaine Only
M33369337186
AYWB399319147
M8498438192
OY-M5173712186
M3333207125
DeVilla387197132
RP166493478205
Zhengzhou630309237
MDGZ-013612611137
QS20225194411245
Data from this study is highlighted in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Toth, R.; Ilic, A.-M.; Huettel, B.; Duduk, B.; Kube, M. Divergence within the Taxon ‘Candidatus Phytoplasma asteris’ Confirmed by Comparative Genome Analysis of Carrot Strains. Microorganisms 2024, 12, 1016. https://doi.org/10.3390/microorganisms12051016

AMA Style

Toth R, Ilic A-M, Huettel B, Duduk B, Kube M. Divergence within the Taxon ‘Candidatus Phytoplasma asteris’ Confirmed by Comparative Genome Analysis of Carrot Strains. Microorganisms. 2024; 12(5):1016. https://doi.org/10.3390/microorganisms12051016

Chicago/Turabian Style

Toth, Rafael, Anna-Marie Ilic, Bruno Huettel, Bojan Duduk, and Michael Kube. 2024. "Divergence within the Taxon ‘Candidatus Phytoplasma asteris’ Confirmed by Comparative Genome Analysis of Carrot Strains" Microorganisms 12, no. 5: 1016. https://doi.org/10.3390/microorganisms12051016

APA Style

Toth, R., Ilic, A. -M., Huettel, B., Duduk, B., & Kube, M. (2024). Divergence within the Taxon ‘Candidatus Phytoplasma asteris’ Confirmed by Comparative Genome Analysis of Carrot Strains. Microorganisms, 12(5), 1016. https://doi.org/10.3390/microorganisms12051016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop