'Candidatus Phytoplasma phoenicium' associated with almond witches'-broom disease: from draft genome to genetic diversity among strain populations

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

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.

The genome of Candidatus phytoplasma ziziphi provides insights into their biological characteristics

Phytoplasmas are obligate cell wall-less prokaryotic bacteria that primarily multiply in plant phloem tissue. Jujube witches' broom (JWB) associated with phytoplasma is a destructive disease of jujube (Ziziphus jujuba Mill.). Here we report the complete 'Candidatus Phytoplasma ziziphi' chromosome of strain Hebei-2018, which is a circular genome of 764,108-base pairs with 735 predicted CDS. Notably, extra 19,825 bp (from 621,995 to 641,819 bp) compared to the previously reported one complements the genes involved in glycolysis, such as pdhA, pdhB, pdhC, pdhD, ackA, pduL and LDH. The synonymous codon usage bias (CUB) patterns by using comparative genomics analysis among the 9 phytoplasmas were similar for most codons. The ENc-GC3s analysis among the 9 phytoplasmas showed a greater effect under the selection on the CUBs of phytoplasmas genes than mutation and other factors. The genome exhibited a strongly reduced ability in metabolic synthesis, while the genes encoding transporter systems were well developed. The genes involved in sec-dependent protein translocation system were also identified.The expressions of nine FtsHs encoding membrane associated ATP-dependent Zn proteases and Mn-SodA with redox capacity in the Ca. P. ziziphi was positively correlated with the phytoplasma concentration. Taken together, the genome will not only expand the number of phytoplasma species and provide some new information about Ca. P. ziziphi, but also contribute to exploring its pathogenic mechanism.

Phytoplasma Taxonomy: Nomenclature, Classification, and Identification

Simple Summary

Phytoplasmas are vector-borne and graft-transmissible bacteria that cause various plant diseases, leading to severe economic losses. Since phytoplasmas cannot be cultured in cell-free media, their identification and taxonomy rely on molecular techniques and gene sequences. In this article, we summarize the recent advances in phytoplasma taxonomy from three different aspects, including (i) nomenclature (naming Candidatus Phytoplasma species); (ii) classification (group and subgroup assignment based on 16S rRNA gene sequences); and (iii) identification (fine differentiation of phytoplasma strains). In addition, some important issues, especially those related to recognizing new ‘Candidatus Phytoplasma’ species, are discussed. This information will be helpful for rapid diagnosis of phytoplasma diseases and accurate taxonomic identification of both emerging and known phytoplasma strains.

Abstract

Phytoplasmas are pleomorphic, wall-less intracellular bacteria that can cause devastating diseases in a wide variety of plant species. Rapid diagnosis and precise identification of phytoplasmas responsible for emerging plant diseases are crucial to preventing further spread of the diseases and reducing economic losses. Phytoplasma taxonomy (identification, nomenclature, and classification) has lagged in comparison to culturable bacteria, largely due to lack of axenic phytoplasma culture and consequent inaccessibility of phenotypic characteristics. However, the rapid expansion of molecular techniques and the advent of high throughput genome sequencing have tremendously enhanced the nucleotide sequence-based phytoplasma taxonomy. In this article, the key events and milestones that shaped the current phytoplasma taxonomy are highlighted. In addition, the distinctions and relatedness of two parallel systems of ‘Candidatus phytoplasma’ species/nomenclature system and group/subgroup classification system are clarified. Both systems are indispensable as they serve different purposes. Furthermore, some hot button issues in phytoplasma nomenclature are also discussed, especially those pertinent to the implementation of newly revised guidelines for ‘Candidatus Phytoplasma’ species description. To conclude, the challenges and future perspectives of phytoplasma taxonomy are briefly outlined.

The Diversity, Distribution and Status of Phytoplasma Diseases in China

Phytoplasmas are important prokaryotic pathogenic bacteria without cell walls, which were formerly known as mycoplasma-like organisms, and belong to the Mollicutes class, Candidatus Phytoplasma genus. They are widely distributed in plants and insects, and can cause serious diseases in important food crops, vegetables, fruit trees, ornamental plants and trees, resulting in huge economic losses. To date, more than 100 phytoplasma diseases have been reported in China, which are distributed throughout the country. Jujube witches'-broom, paulownia witches'-broom, wheat blue dwarf, banana bunchy top, sugarcane white leaf, rice orange leaf and mulberry dwarf represent the phytoplasma diseases causing the most serious damage in China. New phytoplasma diseases and their strains are being reported continuously, indicating that phytoplasmas are more diverse than previously thought. Phytoplasmas are mainly transmitted by insect vectors, such as leafhopper and planthopper, and can also be spread by grafting or Cuscuta australis (known as dodder). Mixed infections of phytoplasmas and viruses, bacteria, and spiroplasmas have also become a serious problem in several crops and are responsible for more synergistic losses. With the continuous development and improvement of technology, molecular biological detection has become the main technique for phytoplasma detection and identification. Currently, research on phytoplasma diseases in China mainly focuses on pathogen identification and classification, and insect vector and host diversity; however, there is less focus on pathogenicity, comparative genomics, and effect factors. More research attention has been paid to wheat blue dwarf phytoplasma, paulownia witches'-broom phytoplasma, jujube witches'-broom phytoplasma, and sugarcane white leaf phytoplasma. Other phytoplasma diseases have been reported; however, there have been no in-depth studies. In this paper, the history and present situation of phytoplasma research, and the status, distribution, and diversity of phytoplasma diseases are summarized, and some possible research directions of phytoplasma in the future in China are proposed.

Comparison of Current Methods for Signal Peptide Prediction in Phytoplasmas

Although phytoplasma studies are still hampered by the lack of axenic cultivation methods, the availability of genome sequences allowed dramatic advances in the characterization of the virulence mechanisms deployed by phytoplasmas, and highlighted the detection of signal peptides as a crucial step to identify effectors secreted by phytoplasmas. However, various signal peptide prediction methods have been used to mine phytoplasma genomes, and no general evaluation of these methods is available so far for phytoplasma sequences. In this work, we compared the prediction performance of SignalP versions 3.0, 4.0, 4.1, 5.0 and Phobius on several sequence datasets originating from all deposited phytoplasma sequences. SignalP 4.1 with specific parameters showed the most exhaustive and consistent prediction ability. However, the configuration of SignalP 4.1 for increased sensitivity induced a much higher rate of false positives on transmembrane domains located at N-terminus. Moreover, sensitive signal peptide predictions could similarly be achieved by the transmembrane domain prediction ability of TMHMM and Phobius, due to the relatedness between signal peptides and transmembrane regions. Beyond the results presented herein, the datasets assembled in this study form a valuable benchmark to compare and evaluate signal peptide predictors in a field where experimental evidence of secretion is scarce. Additionally, this study illustrates the utility of comparative genomics to strengthen confidence in bioinformatic predictions.

Functional variation in phyllogen, a phyllody-inducing phytoplasma effector family, attributable to a single amino acid polymorphism

Flower malformation represented by phyllody is a common symptom of phytoplasma infection induced by a novel family of phytoplasma effectors called phyllogens. Despite the accumulation of functional and structural phyllogen information, the molecular mechanisms of phyllody have not yet been integrated with their evolutionary aspects due to the limited data on their homologs across diverse phytoplasma lineages. Here, we developed a novel universal PCR-based approach to identify 25 phytoplasma phyllogens related to nine "Candidatus Phytoplasma" species, including four species whose phyllogens have not yet been identified. Phylogenetic analyses showed that the phyllogen family consists of four groups (phyl-A, -B, -C, and -D) and that the evolutionary relationships of phyllogens were significantly distinct from those of phytoplasmas, suggesting that phyllogens were transferred horizontally among phytoplasma strains and species. Although phyllogens belonging to the phyl-A, -C, and -D groups induced phyllody, the phyl-B group lacked the ability to induce phyllody. Comparative functional analyses of phyllogens revealed that a single amino acid polymorphism in phyl-B group phyllogens prevented interactions between phyllogens and A- and E-class MADS domain transcription factors (MTFs), resulting in the inability to degrade several MTFs and induce phyllody. Our finding of natural variation in the function of phytoplasma effectors provides new insights into molecular mechanisms underlying the aetiology of phytoplasma diseases.

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.

Assembly of Phytoplasma Genome Drafts from Illumina Reads Using Phytoassembly

Genome drafts for the phytoplasmas may be rapidly and efficiently assembled from NGS sequence data alone exploiting the proper bioinformatic tools and starting from properly collected samples. Here, we describe the use of the Phytoassembly pipeline ( https://github.com/cpolano/phytoassembly ), a fully automated tool that accepts as input row Illumina data from two samples (a phytoplasma infected sample and a healthy reference sample) to produce a phytoplasma genome draft, using the healthy plant host genome as a filter and profiting from the difference in reads coverage between the genome of the pathogen and that of the host. For phytoplasma infected samples containing >2% of pathogen DNA and an isogenic healthy reference sequence the resulting assemblies span the almost entire genomes.

The genome of 'Candidatus Phytoplasma solani' strain SA-1 is highly dynamic and prone to adopting foreign sequences

Bacteria of the genus 'Candidatus Phytoplasma' are uncultivated intracellular plant pathogens transmitted by phloem-feeding insects. They have small genomes lacking genes for essential metabolites, which they acquire from either plant or insect hosts. Nonetheless, some phytoplasmas, such as 'Ca. P. solani', have broad plant host range and are transmitted by several polyphagous insect species. To understand better how these obligate symbionts can colonize such a wide range of hosts, the genome of 'Ca. P. solani' strain SA-1 was sequenced from infected periwinkle via a metagenomics approach. The de novo assembly generated a draft genome with 19 contigs totalling 821,322bp, which corresponded to more than 80% of the estimated genome size. Further completion of the genome was challenging due to the high occurrence of repetitive sequences. The majority of repeats consisted of gene arrangements characteristic of phytoplasma potential mobile units (PMUs). These regions showed variation in gene orders intermixed with genes of unknown functions and lack of similarity to other phytoplasma genes, suggesting that they were prone to rearrangements and acquisition of new sequences via recombination. The availability of this high-quality draft genome also provided a foundation for genome-scale genotypic analysis (e.g., average nucleotide identity and average amino acid identity) and molecular phylogenetic analysis. Phylogenetic analyses provided evidence of horizontal transfer for PMU-like elements from various phytoplasmas, including distantly related ones. The 'Ca. P. solani' SA-1 genome also contained putative secreted protein/effector genes, including a homologue of SAP11, found in many other phytoplasma species.

Comparative genome analysis of jujube witches'-broom Phytoplasma, an obligate pathogen that causes jujube witches'-broom disease

BACKGROUND

JWB phytoplasma is a kind of insect-transmitted and uncultivable bacterial plant pathogen causeing a destructive Jujube disease. To date, no genome information about JWB phytoplasma has been published, which hindered its characterization at genomic level. To understand its pathogenicity and ecology, the genome of a JWB phytoplasma isolate jwb-nky was sequenced and compared with other phytoplasmas enabled us to explore the mechanisms of genomic rearrangement.

RESULTS

The complete genome sequence of JWB phytoplasma (jwb-nky) was determined, which consisting of one circular chromosome of 750,803 bp with a GC content of 23.3%. 694 protein-encoding genes, 2 operons for rRNA genes and 31 tRNA genes as well as 4 potential mobile units (PMUs) containing clusters of DNA repeats were identified. Based on PHIbaes analysis, a large number of genes were genome-specific and approximately 13% of JWB phytoplasma genes were predicted to be associated with virulence. Although transporters for maltose, dipeptides/oligopeptides, spermidine/putrescine, cobalt, Mn/Zn and methionine were identified, KEGG pathway analysis revealed the reduced metabolic capabilities of JWB phytoplasma. Comparative genome analyses between JWB phytoplasma and other phytoplasmas shows the occurrence of large-scale gene rearrangements. The low synteny with other phytoplasmas indicated that the expansion of multiple gene families/duplication probably occurred separately after differentiation.

CONCLUSIONS

In this study, the complete genome sequence of a JWB phytoplasma isolate jwb-nky that causing JWB disease was reported for the first time and a number of species-specific genes were identified in the genome. The study enhanced our understandings about genomic basis and the pathogenicity mechanism of this pathogen, which will aid in the development of improved strategies for efficient management of JWB diseases.

Draft genome sequence of the New Jersey aster yellows strain of 'Candidatus Phytoplasma asteris'

The NJAY (New Jersey aster yellows) strain of ‘Candidatus Phytoplasma asteris’ is a significant plant pathogen responsible for causing severe lettuce yellows in the U.S. state of New Jersey. A draft genome sequence was prepared for this organism. A total of 177,847 reads were assembled into 75 contigs > 518 bp with a total base value of 652,092 and an overall [G+C] content of 27.1%. A total of 733 protein coding genes were identified. This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the accession MAPF00000000. This draft genome was used for genome- and gene-based comparative phylogenetic analyses with other phytoplasmas, including the closely related ‘Ca. Phytoplasma asteris’ strain, aster yellows witches’- broom (AY-WB). NJAY and AY-WB exhibit approximately 0.5% dissimilarity at the nucleotide level among their shared genomic segments. Evidence indicated that NJAY harbors four plasmids homologous to those known to encode pathogenicity determinants in AY-WB, as well as a chromosome-encoded mobile unit. Apparent NJAY orthologs to the important AY-WB virulence factors, SAP11 and SAP54, were identified. A number of secreted proteins, both membrane-bound and soluble, were encoded, with many bearing similarity to known AY-WB effector molecules and others representing possible secreted proteins that may be novel to the NJAY lineage.

Molecular identification and characterization of the new 16SrIX-J and cpn60 UT IX-J phytoplasma subgroup associated with chicory bushy stunt disease in Saudi Arabia

Chicory (Cichorium intybus) is a perennial plant (Asteraceae) that grows wild in pasture fields in Saudi Arabia. Chicory plants displaying symptoms typically induced by phytoplasmas, such as bushy phenotype and stunt, were observed in the Mulayda region, Qassim governorate, Saudi Arabia. In this study we examined samples taken from three symptomatic chicory plants and confirmed the presence of phytoplasma DNA. Analysis of the 16S rRNA-encoding sequences showed that the plants were infected with a phytoplasma from the pigeon pea witches'-broom group (16SrIX). Sequencing of the 16S rRNA-encoding gene and the partial cpn60 sequence, computer-simulated RFLP analysis, and phylogenetic analysis of both markers revealed that the phytoplasma identified was representative of a new 16SrIX-J and cpn60 UT IX-IJ subgroup. The present study identified chicory plants as a novel host for phytoplasma strains within the pigeon pea witches'-broom phytoplasma group, and expanded the known diversity of this group.

A few sequence polymorphisms among isolates of Maize bushy stunt phytoplasma associate with organ proliferation symptoms of infected maize plants

Background and Aims

Maize bushy stunt phytoplasma (MBSP) is a bacterial pathogen of maize ( Zea mays L.) across Latin America. MBSP belongs to the 16SrI-B sub-group within the genus ' Candidatus Phytoplasma'. MBSP and its insect vector Dalbulus maidis (Hemiptera: Cicadellidae) are restricted to maize; both are thought to have coevolved with maize during its domestication from a teosinte-like ancestor. MBSP-infected maize plants show a diversity of symptoms. and it is likely that MBSP is under strong selection for increased virulence and insect transmission on maize hybrids that are widely grown in Brazil. In this study it was investigated whether the differences in genome sequences of MBSP isolates from two maize-growing regions in South-east Brazil explain variations in symptom severity of the MBSP isolates on various maize genotypes.

Methods

MBSP isolates were collected from maize production fields in Guaíra and Piracicaba in South-east Brazil for infection assays. One representative isolate was chosen for de novo whole-genome assembly and for the alignment of sequence reads from the genomes of other phytoplasma isolates to detect polymorphisms. Statistical methods were applied to investigate the correlation between variations in disease symptoms of infected maize plants and MBSP sequence polymorphisms.

Key Results

MBSP isolates contributed consistently to organ proliferation symptoms and maize genotype to leaf necrosis, reddening and yellowing of infected maize plants. The symptom differences are associated with polymorphisms in a phase-variable lipoprotein, which is a candidate effector, and an ATP-dependent lipoprotein ABC export protein, whereas no polymorphisms were observed in other candidate effector genes. Lipoproteins and ABC export proteins activate host defence responses, regulate pathogen attachment to host cells and activate effector secretion systems in other pathogens.

Conclusions

Polymorphisms in two putative virulence genes among MBSP isolates from maize-growing regions in South-east Brazil are associated with variations in organ proliferation symptoms of MBSP-infected maize plants.

Phytoplasma classification and phylogeny based on in silico and in vitro RFLP analysis of cpn60 universal target sequences

Phytoplasmas are unculturable, phytopathogenic bacteria that cause economic losses worldwide. As unculturable micro-organisms, phytoplasma taxonomy has been based on the use of the 16S rRNA-encoding gene to establish 16Sr groups and subgroups based on the restriction fragment length polymorphism (RFLP) pattern resulting from the digestion of amplicon (in vitro) or sequence (in silico) with seventeen restriction enzymes. Problems such as heterogeneity of the ribosomal operon and the inability to differentiate closely related phytoplasma strains has motivated the search for additional markers capable of providing finer differentiation of phytoplasma strains. In this study we developed and validated a scheme to classify phytoplasmas based on the use of cpn60 universal target (cpn60 UT) sequences. Ninety-six cpn60 UT sequences from strains belonging to 19 16Sr subgroups were subjected to in silico RFLP using pDRAW32 software, resulting in 25 distinctive RFLP profiles. Based on these results we delineated cpn60 UT groups and subgroups, and established a threshold similarity coefficient for groups and subgroups classifying all the strains analysed in this study. The nucleotide identity among the reference strains, the correspondence between in vitro and in silico RFLP, and the phylogenetic relationships of phytoplasma strains based on cpn60 UT sequences are also discussed.