Genome wide sequence analysis grants unbiased definition of species boundaries in "Candidatus Phytoplasma"

Multiple factors driving the acquisition efficiency of apple proliferation phytoplasma in Cacopsylla melanoneura

Phytoplasmas are bacterial pathogens located in the plant's phloem that are responsible for several plant diseases and are mainly transmitted by phloem-sucking insects. Apple proliferation (AP) is an economically important disease associated with the presence of 'Candidatus Phytoplasma mali' which is transmitted by two psyllid species. While Cacopsylla picta is a vector in different regions, the vector efficiency of C. melanoneura varies between different populations. This species is considered the main AP vector in Northwestern Italy but plays a minor role in Northeastern Italy and other European regions. To investigate whether the psyllid and/or the phytoplasma subtype drive the phytoplasma acquisition in C. melanoneura, a phytoplasma acquisition experiment was set up using single mating couples of overwintered individuals from different psyllid populations and phytoplasma subtypes. All analyzed insect populations acquired phytoplasma, but with different efficiencies and concentrations. The main factors driving the acquisition were the phytoplasma subtype and its concentration in the leaves of the infected trees together with the psyllid lineage. The phytoplasma concentration in the psyllids was again influenced by the phytoplasma subtype, the psyllid lineage and the region of origin, whereas the phytoplasma concentration in the leaves and the psyllid haplotype defined with the cytochrome oxidase I gene had only a minor impact on the phytoplasma concentration. This is the first study evaluating the roles of both the psyllid haplotype and the phytoplasma subtype on the acquisition process and highlights the importance of C. melanoneura as an additional AP vector.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10340-023-01699-1.

From sequences to species: Charting the phytoplasma classification and taxonomy in the era of taxogenomics

Phytoplasma taxonomy has been a topic of discussion for the last two and half decades. Since the Japanese scientists discovered the phytoplasma bodies in 1967, the phytoplasma taxonomy was limited to disease symptomology for a long time. The advances in DNA-based markers and sequencing improved phytoplasma classification. In 2004, the International Research Programme on Comparative Mycoplasmology (IRPCM)- Phytoplasma/Spiroplasma Working Team – Phytoplasma taxonomy group provided the description of the provisional genus ‘Candidatus Phytoplasma’ with guidelines to describe the new provisional phytoplasma species. The unintentional consequences of these guidelines led to the description of many phytoplasma species where species characterization was restricted to a partial sequence of the 16S rRNA gene alone. Additionally, the lack of a complete set of housekeeping gene sequences or genome sequences, as well as the heterogeneity among closely related phytoplasmas limited the development of a comprehensive Multi-Locus Sequence Typing (MLST) system. To address these issues, researchers tried deducing the definition of phytoplasma species using phytoplasmas genome sequences and the average nucleotide identity (ANI). In another attempts, a new phytoplasma species were described based on the Overall Genome relatedness Values (OGRI) values fetched from the genome sequences. These studies align with the attempts to standardize the classification and nomenclature of ‘Candidatus’ bacteria. With a brief historical account of phytoplasma taxonomy and recent developments, this review highlights the current issues and provides recommendations for a comprehensive system for phytoplasma taxonomy until phytoplasma retains ‘Candidatus’ status.

Multilocus sequence typing of diverse phytoplasmas using hybridization probe-based sequence capture provides high resolution strain differentiation

Phytoplasmas are insect-vectored, difficult-to-culture bacterial pathogens that infect a wide variety of crop and non-crop plants, and are associated with diseases that can lead to significant yield losses in agricultural production worldwide. Phytoplasmas are currently grouped in the provisional genus ‘Candidatus Phytoplasma’, which includes 49 ‘Candidatus’ species. Further differentiation of phytoplasmas into ribosomal groups is based on the restriction fragment length polymorphism (RFLP) pattern of the 16S rRNA-encoding operon, with more than 36 ribosomal groups (16Sr) and over 100 subgroups reported. Since disease symptoms on plants are not associated with phytoplasma identity, accurate diagnostics is of critical importance to manage disease associated with these microorganisms. Phytoplasmas are typically detected from plant and insect tissue using PCR-based methods targeting universal taxonomic markers. Although these methods are relatively sensitive, specific and are widely used, they have limitations, since they provide limited resolution of phytoplasma strains, thus necessitating further assessment of biological properties and delaying implementation of mitigation measures. Moreover, the design of PCR primers that can target multiple loci from phytoplasmas that differ at the sequence level can be a significant challenge. To overcome these limitations, a PCR-independent, multilocus sequence typing (MLST) assay to characterize an array of phytoplasmas was developed. Hybridization probe s targeting cpn60, tuf, secA, secY, and nusA genes, as well as 16S and rp operons, were designed and used to enrich DNA extracts from phytoplasma-infected samples for DNA fragments corresponding to these markers prior to Illumina sequencing. This method was tested using different phytoplasmas including ‘Ca. P. asteris’ (16SrI-B), ‘Ca. P. pruni’ (16SrIII-A),‘Ca. P. prunorum’ (16SrX-B), ‘Ca. P. pyri’ (16SrX-C), ‘Ca. P. mali’ (16SrX-A), and ‘Ca. P. solani’ (16SrXII-A). Thousands of reads were obtained for each gene with multiple overlapping fragments, which were assembled to generate full-length (typically >2 kb), high-quality sequences. Phytoplasma groups and subgroups were accurately determined based on 16S ribosomal RNA and cpn60 gene sequences. Hybridization-based MLST facilitates the enrichment of target genes of phytoplasmas and allows the simultaneous determination of sequences corresponding to seven different markers. In this proof-of-concept study, hybridization-based MLST was demonstrated to be an efficient way to generate data regarding ‘Ca. Phytoplasma’ species/strain differentiation.

Iodixanol density gradients as an effective phytoplasma enrichment approach to improve genome sequencing

Obtaining complete phytoplasma genomes is difficult due to the lack of a culture system for these bacteria. To improve genome assembly, a non-ionic, low- and iso-osmotic iodixanol (Optiprep™) density gradient centrifugation method was developed to enrich for phytoplasma cells and deplete plant host tissues prior to deoxyribonucleic acid (DNA) extraction and high-throughput sequencing (HTS). After density gradient enrichment, potato infected with a ‘Candidatus Phytoplasma australasia’-related strain showed a ∼14-fold increase in phytoplasma HTS reads, with a ∼1.7-fold decrease in host genomic reads compared to the DNA extracted from the same sample without density gradient centrifugation enrichment. Additionally, phytoplasma genome assemblies from libraries equalized to 5 million reads were, on average, ∼15,000 bp larger and more contiguous (N50 ∼14,800 bp larger) than assemblies from the DNA extracted from the infected potato without enrichment. The method was repeated on capsicum infected with Sweet Potato Little Leaf phytoplasma (‘Ca. Phytoplasma australasia’-related strain) with a lower phytoplasma titer than the potato. In capsicum, ∼threefold more phytoplasma reads and ∼twofold less host genomic reads were obtained, with the genome assembly size and N50 values from libraries equalized to 3.4 million reads ∼137,000 and ∼4,000 bp larger, respectively, compared to the DNA extracted from infected capsicum without enrichment. Phytoplasmas from potato and capsicum were both enriched at a density of 1.049–1.058 g/ml. Finally, we present two highly contiguous ‘Ca. Phytoplasma australasia’ phytoplasma reference genomes sequenced from naturally infected Solanaceae hosts in Australia. Obtaining high-quality phytoplasma genomes from naturally infected hosts will improve insights into phytoplasma taxonomy, which will improve their detection and disease management.

The agent associated with blue dwarf disease in wheat represents a new phytoplasma taxon, 'Candidatus Phytoplasma tritici'

Wheat blue dwarf (WBD) is one of the most economically damaging cereal crop diseases in northwestern PR China. The agent associated with the WBD disease is a phytoplasma affiliated with the aster yellows (AY) group, subgroup C (16SrI-C). Since phytoplasma strains within the AY group are ecologically and genetically diverse, it has been conceived that the AY phytoplasma group may consist of more than one species. This communication presents evidence to demonstrate that, while each of the two 16 rRNA genes of the WBD phytoplasma shares >97.5 % sequence similarity with that of the 'Candidatus Phytoplasma asteris' reference strain, the WBD phytoplasma clearly represents an ecologically separated lineage: the WBD phytoplasma not only has its unique transmitting vector (Psammotettix striatus) but also elicits a distinctive symptom in its predominant plant host (wheat). In addition, the WBD phytoplasma possesses molecular characteristics that further manifest its significant divergence from 'Ca. P. asteris'. Such molecular characteristics include lineage-specific antigenic membrane proteins and a lower than 95 % genome-wide average nucleotide identity score with 'Ca. P. asteris'. These ecological, molecular and genomic evidences justify the recognition of the WBD phytoplasma as a novel taxon, 'Candidatus Phytoplasma tritici'.

Competitive Exclusion of Flavescence dorée Phytoplasma Strains in Catharanthus roseus Plants

Flavescence dorée phytoplasmas (FDp, 16SrV-C and -D) are plant pathogenic non-cultivable bacteria associated with a severe grapevine disease. The incidence of the two reference strains on cultivated grapevines is unbalanced, and mixed infections are rare. To investigate the interaction between the two strains, Catharanthus roseus plants were graft-infected with both strains, either simultaneously or sequentially. Different combinations of lateral and apical grafting were applied to avoid possible benefits due to graft position. The infection was monitored for four months through a new diagnostic protocol developed for differentiation and relative quantification of the two strains. Regardless of the temporal or spatial advantage at grafting, FD-C generally outcompeted FD-D. The prevalence of FD-C increased over time and, at the end of the experiment, FD-C was the unique strain detected in the aerial part and the roots of 74% and 90% of grafted plants, respectively. These data indicate that the interaction between the two strains results in competitive exclusion. Understanding the bases of the competition between FD-C and FD-D may contribute to explain the biology of the coexistence of different FDp strains under field conditions, aiming at identifying potential suppressor strains, which can provide alternative and environmentally sustainable solutions for FD control.

Species Boundaries and Molecular Markers for the Classification of 16SrI Phytoplasmas Inferred by Genome Analysis

Phytoplasmas are plant-pathogenic bacteria that impact agriculture worldwide. The commonly adopted classification system for phytoplasmas is based on the restriction fragment length polymorphism (RFLP) analysis of their 16S rRNA genes. With the increased availability of phytoplasma genome sequences, the classification system can now be refined. This work examined 11 strains in the 16SrI group within the genus ‘Candidatus Phytoplasma’ and investigated the possible species boundaries. We confirmed that the RFLP classification method is problematic due to intragenomic variation of the 16S rRNA genes and uneven weighing of different nucleotide positions. Importantly, our results based on the molecular phylogeny, differentiations in chromosomal segments and gene content, and divergence in homologous sequences, all supported that these strains may be classified into multiple operational taxonomic units (OTUs) equivalent to species. Strains assigned to the same OTU share >97% genome-wide average nucleotide identity (ANI) and >78% of their protein-coding genes. In comparison, strains assigned to different OTUs share < 94% ANI and < 75% of their genes. Reduction in homologous recombination between OTUs is one possible explanation for the discontinuity in genome similarities, and these findings supported the proposal that 95% ANI could serve as a cutoff for distinguishing species in bacteria. Additionally, critical examination of these results and the raw sequencing reads led to the identification of one genome that was presumably mis-assembled by combining two sequencing libraries built from phytoplasmas belonging to different OTUs. This finding provided a cautionary tale for working on uncultivated bacteria. Based on the new understanding of phytoplasma divergence and the current genome availability, we developed five molecular markers that could be used for multilocus sequence analysis (MLSA). By selecting markers that are short yet highly informative, and are distributed evenly across the chromosome, these markers provided a cost-effective system that is robust against recombination. Finally, examination of the effector gene distribution further confirmed the rapid gains and losses of these genes, as well as the involvement of potential mobile units (PMUs) in their molecular evolution. Future improvements on the taxon sampling of phytoplasma genomes will allow further expansions of similar analysis, and thus contribute to phytoplasma taxonomy and diagnostics.

Flavescence Dorée Phytoplasma Has Multiple ftsH Genes that Are Differentially Expressed in Plants and Insects

Flavescence dorée (FD) is a severe epidemic disease of grapevines caused by FD phytoplasma (FDP) transmitted by the leafhopper vector Scaphoideus titanus. The recent sequencing of the 647-kbp FDP genome highlighted an unusual number of genes encoding ATP-dependent zinc proteases FtsH, which have been linked to variations in the virulence of "Candidatus Phytoplasma mali" strains. The aims of the present study were to predict the FtsH repertoire of FDP, to predict the functional domains and topologies of the encoded proteins in the phytoplasma membrane and to measure the expression profiles in different hosts. Eight complete ftsH genes have been identified in the FDP genome. In addition to ftsH6, which appeared to be the original bacterial ortholog, the other seven gene copies were clustered on a common distinct phylogenetic branch, suggesting intra-genome duplication of ftsH. The expression of these proteins, quantified in plants and insect vectors in natural and experimental pathosystems, appeared to be modulated in a host-dependent manner. Two of the eight FtsH C-tails were predicted by Phobius software to be extracellular and, therefore, in direct contact with the host cellular content. As phytoplasmas cannot synthesize amino acids, our data raised questions regarding the involvement of FtsH in the adaptation to hosts via potentially enhanced recycling of phytoplasma cellular proteins and host protein degradation.

CReasPy-Cloning: A Method for Simultaneous Cloning and Engineering of Megabase-Sized Genomes in Yeast Using the CRISPR-Cas9 System

Over the past decade, a new strategy was developed to bypass the difficulties to genetically engineer some microbial species by transferring (or "cloning") their genome into another organism that is amenable to efficient genetic modifications and therefore acts as a living workbench. As such, the yeast Saccharomyces cerevisiae has been used to clone and engineer genomes from viruses, bacteria, and algae. The cloning step requires the insertion of yeast genetic elements in the genome of interest, in order to drive its replication and maintenance as an artificial chromosome in the host cell. Current methods used to introduce these genetic elements are still unsatisfactory, due either to their random nature (transposon) or the requirement for unique restriction sites at specific positions (TAR cloning). Here we describe the CReasPy-cloning, a new method that combines both the ability of Cas9 to cleave DNA at a user-specified locus and the yeast's highly efficient homologous recombination to simultaneously clone and engineer a bacterial chromosome in yeast. Using the 0.816 Mbp genome of Mycoplasma pneumoniae as a proof of concept, we demonstrate that our method can be used to introduce the yeast genetic element at any location in the bacterial chromosome while simultaneously deleting various genes or group of genes. We also show that CReasPy-cloning can be used to edit up to three independent genomic loci at the same time with an efficiency high enough to warrant the screening of a small (<50) number of clones, allowing for significantly shortened genome engineering cycle times.

Genetic Diversity of Flavescence Dorée Phytoplasmas at the Vineyard Scale

To study the role of wild areas around the vineyards in the epidemiology of flavescence dorée (FD) and track the origin of new foci, two phytoplasma genetic markers, dnaK and malG, were developed for FD phytoplasma (FDp) characterization. The two genes and the vmpA locus were used to genetically characterize FDp populations at seven agroecosystems of a wine-growing Italian region. Vitis vinifera, "gone-wild" V. vinifera and rootstocks, Clematis spp., and Scaphoideus titanus adults were sampled within and outside the vineyards. A range of genotypes infecting the different hosts of the FDp epidemiological cycle was found. Type FD-C isolates were fairly homogeneous compared to type FD-D ones. Most of the FD-D variability was correlated with the malG sequence, and a duplication of this locus was demonstrated for this strain. Coinfection with FD-C and FD-D strains was rare, suggesting possible competition between the two. Similar levels of FDp genetic variation recorded for grapevines or leafhoppers of cultivated and wild areas and co-occurrence of many FDp genotypes inside and outside the vineyards supported the idea of the importance of wild or abandoned Vitis plants and associated S. titanus insects in the epidemiology of the disease. Genetic profiles of FDp found in Clematis were never found in the other hosts, indicating that this species does not take part in the disease cycle in the area. Due to the robustness of analyses using dnaK for discriminating between FD-C and FD-D strains and the high variability of malG sequences, these are efficient markers to study FDp populations and epidemiology at a small geographical scale.IMPORTANCE Flavescence dorée, a threatening disease of grapevine caused by FD phytoplasma (FDp), is distributed within the most important wine-producing areas of Europe and has severe effects on both vineyard productivity and landscape management. FDp is a quarantine pest in Europe, and despite the efforts to contain the pathogen, the disease is still spreading. In this work, new genetic markers for the fine genetic characterization of FDp at local scale are presented. Our findings improve the knowledge of FDp epidemiological cycle and offer the possibility of tracking the route of the FDp infection. In particular, due to its high genetic variability, one of the newly developed markers could be sufficient to track the origin of new infection foci, either from the wild areas or from nurseries.

Multilocus Sequence Analysis Reveals Three Distinct Populations of "Candidatus Phytoplasma palmicola" with a Specific Geographical Distribution on the African Continent

Coconut is an important crop for both industry and small stakeholders in many intertropical countries. Phytoplasma-associated lethal yellowing-like diseases have become one of the major pests that limit coconut cultivation as they have emerged in different parts of the world. We developed a multilocus sequence typing scheme (MLST) for tracking epidemics of “Ca. Phytoplasma palmicola,” which is responsible for coconut lethal yellowing disease (CLYD) on the African continent. MLST analysis applied to diseased coconut samples collected in western and eastern African countries also showed the existence of three distinct populations of “Ca. Phytoplasma palmicola” with low intrapopulation diversity. The reasons for the observed strong geographic patterns remain to be established but could result from the lethality of CLYD and the dominance of short-distance insect-mediated transmission.

To sustain epidemiological studies on coconut lethal yellowing disease (CLYD), a devastating disease in Africa caused by a phytoplasma, we developed a multilocus sequence typing (MLST) scheme for “Candidatus Phytoplasma palmicola” based on eight housekeeping genes. At the continental level, eight different sequence types were identified among 132 “Candidatus Phytoplasma palmicola”-infected coconuts collected in Ghana, Nigeria, and Mozambique, where CLYD epidemics are still very active. “Candidatus Phytoplasma palmicola” appeared to be a bacterium that is subject to strong bottlenecks, reducing the fixation of positively selected beneficial mutations into the bacterial population. This phenomenon, as well as a limited plant host range, might explain the observed country-specific distribution of the eight haplotypes. As an alternative means to increase fitness, bacteria can also undergo genetic exchange; however, no evidence for such recombination events was found for “Candidatus Phytoplasma palmicola.” The implications for CLYD epidemiology and prophylactic control are discussed. The usefulness of seven housekeeping genes to investigate the genetic diversity in the genus “Candidatus Phytoplasma” is underlined.

IMPORTANCE Coconut is an important crop for both industry and small stakeholders in many intertropical countries. Phytoplasma-associated lethal yellowing-like diseases have become one of the major pests that limit coconut cultivation as they have emerged in different parts of the world. We developed a multilocus sequence typing scheme (MLST) for tracking epidemics of “Ca. Phytoplasma palmicola,” which is responsible for coconut lethal yellowing disease (CLYD) on the African continent. MLST analysis applied to diseased coconut samples collected in western and eastern African countries also showed the existence of three distinct populations of “Ca. Phytoplasma palmicola” with low intrapopulation diversity. The reasons for the observed strong geographic patterns remain to be established but could result from the lethality of CLYD and the dominance of short-distance insect-mediated transmission.

An Effective Pipeline Based on Relative Coverage for the Genome Assembly of Phytoplasmas and Other Fastidious Prokaryotes

BACKGROUND

For the plant pathogenic phytoplasmas, as well as for several fastidious prokaryotes, axenic cultivation is extremely difficult or not possible yet; therefore, even with second generation sequencing methods, obtaining the sequence of their genomes is challenging due to host sequence contamination.

OBJECTIVE

With the Phytoassembly pipeline here presented, we aim to provide a method to obtain high quality genome drafts for the phytoplasmas and other uncultivable plant pathogens, by exploiting the coverage differential in the ILLUMINA sequences from the pathogen and the host, and using the sequencing of a healthy, isogenic plant as a filter.

VALIDATION

The pipeline has been benchmarked using simulated and real ILLUMINA runs from phytoplasmas whose genome is known, and it was then used to obtain high quality drafts for three new phytoplasma genomes.

CONCLUSION

For phytoplasma infected samples containing >2-4% of pathogen DNA and an isogenic reference healthy sample, the resulting assemblies can be next to complete. The Phytoassembly source code is available on GitHub at https://github.com/cpolano/phytoassembly.

Molecular diagnostic tools for detection and differentiation of phytoplasmas based on chaperonin-60 reveal differences in host plant infection patterns

Phytoplasmas (‘Candidatus Phytoplasma’ spp.) are insect-vectored bacteria that infect a wide variety of plants, including many agriculturally important species. The infections can cause devastating yield losses by inducing morphological changes that dramatically alter inflorescence development. Detection of phytoplasma infection typically utilizes sequences located within the 16S–23S rRNA-encoding locus, and these sequences are necessary for strain identification by currently accepted standards for phytoplasma classification. However, these methods can generate PCR products >1400 bp that are less divergent in sequence than protein-encoding genes, limiting strain resolution in certain cases. We describe a method for accessing the chaperonin-60 (cpn60) gene sequence from a diverse array of ‘Ca.Phytoplasma’ spp. Two degenerate primer sets were designed based on the known sequence diversity of cpn60 from ‘Ca.Phytoplasma’ spp. and used to amplify cpn60 gene fragments from various reference samples and infected plant tissues. Forty three cpn60 sequences were thereby determined. The cpn60 PCR-gel electrophoresis method was highly sensitive compared to 16S-23S-targeted PCR-gel electrophoresis. The topology of a phylogenetic tree generated using cpn60 sequences was congruent with that reported for 16S rRNA-encoding genes. The cpn60 sequences were used to design a hybridization array using oligonucleotide-coupled fluorescent microspheres, providing rapid diagnosis and typing of phytoplasma infections. The oligonucleotide-coupled fluorescent microsphere assay revealed samples that were infected simultaneously with two subtypes of phytoplasma. These tools were applied to show that two host plants, Brassica napus and Camelina sativa, displayed different phytoplasma infection patterns.