Comparative Genome Analysis of ‘Candidatus Phytoplasma luffae’ Reveals the Influential Roles of Potential Mobile Units in Phytoplasma Evolution

Virulence effectors encoded in the rice yellow dwarf phytoplasma genome participate in pathogenesis

Bacteria-like phytoplasmas alternate between plant and insect hosts, secreting proteins that disrupt host development. In this study, we sequenced the complete genome of "Candidatus Phytoplasma oryzae" strain HN2022, associated with rice yellow dwarf (RYD) disease, using PacBio HiFi technology. The strain was classified within the 16Sr XI-B subgroup. Through SignalP v5.0 for prediction and subsequent expression analysis of secreted proteins in Nicotiana benthamiana and rice (Oryza sativa L.), we identified the key virulence effector proteins RY348 and RY378. RY348, a homolog of Secreted Aster Yellows Phytoplasma Effector 54 (SAP54), targets and degrades the MADS-box transcription factors MADS1 and MADS15, causing pollen sterility. Meanwhile, RY378 impacts the strigolactone and auxin signaling pathways, substantially increasing tillering. These findings offer insights into the interactions between plants and phytoplasmas.

Phytoplasma DNA Enrichment from Sugarcane White Leaves for Shotgun Sequencing Improvement

Sugarcane white leaf (SCWL) disease, caused by Candidatus Phytoplasma sacchari, poses a significant threat to sugarcane cultivation. An obligate parasite, phytoplasma is difficult to culture in laboratory conditions, making the isolation of its DNA from the massive amount of plant host DNA extremely challenging. Yet, the appropriate amount and quality of plant microbiome-derived DNA are key for high-quality DNA sequencing data. Here, a simple, cost-effective, alternative method for DNA isolation was applied using a guanidine-HCl-hydroxylated silica (GuHCl-Silica)-based method and microbiome DNA enrichment based on size-selective low-molecular-weight (LMW) DNA by PEG/NaCl precipitation. qPCR analysis revealed a significant enrichment of phytoplasma DNA in the LMW fraction. Additionally, the NEBNext Microbiome DNA enrichment kit was utilized to further enrich microbial DNA, demonstrating a remarkable increase in the relative abundance of phytoplasma DNA to host DNA. Shotgun sequencing of the isolated DNA gave high-quality data on the metagenome assembly genome (MAG) of Ca. Phytoplasma sacchari SCWL with completeness at 95.85 and 215× coverage. The results indicate that this combined approach of PEG/NaCl size selection and microbiome enrichment is effective for obtaining high-quality genomic data from phytoplasma, surpassing previous methods in efficiency and resource utilization. This low-cost method not only enhances the recovery of microbiome DNA from plant hosts but also provides a robust framework for studying plant pathogens in complex plant models.

Phytoplasma: A plant pathogen that cannot be ignored in agricultural production-Research progress and outlook

Phytoplasmas are phloem-restricted plant-pathogenic bacteria transmitted by insects. They cause diseases in a wide range of host plants, resulting in significant economic and ecological losses worldwide. Research on phytoplasmas has a long history, with significant progress being made in the past 30 years. Notably, with the rapid development of phytoplasma research, scientists have identified the primary agents involved in phytoplasma transmission, established classification and detection systems for phytoplasmas, and 243 genomes have been sequenced and assembled completely or to draft quality. Multiple possible phytoplasma effectors have been investigated, elucidating the molecular mechanisms by which phytoplasmas manipulate their hosts. This review summarizes recent advances in phytoplasma research, including identification techniques, host range studies, whole- or draft-genome sequencing, effector pathogenesis and disease control methods. Additionally, future research directions in the field of phytoplasma research are discussed.

Improving the Comprehension of Pathogenicity and Phylogeny in 'Candidatus Phytoplasma meliae' through Genome Characterization

'Candidatus Phytoplasma meliae' is a pathogen associated with chinaberry yellowing disease, which has become a major phytosanitary problem for chinaberry forestry production in Argentina. Despite its economic impact, no genome information of this phytoplasma has been published, which has hindered its characterization at the genomic level. In this study, we used a metagenomics approach to analyze the draft genome of the 'Ca. P. meliae' strain ChTYXIII. The draft assembly consisted of twenty-one contigs with a total length of 751.949 bp, and annotation revealed 669 CDSs, 34 tRNAs, and 1 set of rRNA operons. The metabolic pathways analysis showed that ChTYXIII contains the complete core genes for glycolysis and a functional Sec system for protein translocation. Our phylogenomic analysis based on 133 single-copy genes and genome-to-genome metrics supports the classification as unique 'Ca. P. species' within the MPV clade. We also identified 31 putative effectors, including a homolog to SAP11 and others that have only been described in this pathogen. Our ortholog analysis revealed 37 PMU core genes in the genome of 'Ca. P. meliae' ChTYXIII, leading to the identification of 2 intact PMUs. Our work provides important genomic information for 'Ca. P. meliae' and others phytoplasmas for the 16SrXIII (MPV) group.

PhyEffector, the First Algorithm That Identifies Classical and Non-Classical Effectors in Phytoplasmas

Phytoplasmas are the causal agents of more than 100 plant diseases in economically important crops. Eleven genomes have been fully sequenced and have allowed us to gain a better understanding of the biology and evolution of phytoplasmas. Effectors are key players in pathogenicity and virulence, and their identification and description are becoming an essential practice in the description of phytoplasma genomes. This is of particular importance because effectors are possible candidates for the development of new strategies for the control of plant diseases. To date, the prediction of effectors in phytoplasmas has been a great challenge; the reliable comparison of effectoromes has been hindered because research teams have used the combination of different programs in their predictions. This is not trivial since significant differences in the results can arise, depending on the predictive pipeline used. Here, we tested different predictive pipelines to create the PhyEffector algorithm; the average value of the F1 score for PhyEffector was 0.9761 when applied to different databases or genomes, demonstrating its robustness as a predictive tool. PhyEffector can recover both classical and non-classical phytoplasma effectors, making it an invaluable tool to accelerate effectoromics in phytoplasmas.

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.

Complete genome sequence of "Candidatus Phytoplasma sacchari" obtained using a filter-based DNA enrichment method and Nanopore sequencing

Phytoplasmas are phloem-limited plant pathogens, such as sugarcane white leaf (SCWL) phytoplasma, which are responsible for heavy economic losses to the sugarcane industry. Characterization of phytoplasmas has been limited because they cannot be cultured in vitro. However, with the advent of genome sequencing, different aspects of phytoplasmas are being investigated. In this study, we developed a DNA enrichment method for sugarcane white leaf (SCWL) phytoplasma, evaluated the effect of DNA enrichment via Illumina sequencing technologies, and utilized Illumina and Nanopore sequencing technologies to obtain the complete genome sequence of the "Candidatus Phytoplasma sacchari" isolate SCWL1 that is associated with sugarcane white leaf in China. Illumina sequencing analysis elucidated that only 1.21% of the sequencing reads from total leaf DNA were mapped to the SCWL1 genome, whereas 40.97% of the sequencing reads from the enriched DNA were mapped to the SCWL1 genome. The genome of isolate SCWL1 consists of a 538,951 bp and 2976 bp long circular chromosome and plasmid, respectively. We identified 459 protein-encoding genes, 2 complete 5S-23S-16S rRNA gene operons, 27 tRNA genes, and an incomplete potential mobile unit (PMU) in the circular chromosome. Phylogenetic analyses and average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values based on the sequenced genome revealed that SCWL phytoplasma and sugarcane grassy shoot (SCGS) phytoplasma belonged to the same phytoplasma species. This study provides a genomic DNA enrichment method for phytoplasma sequencing. Moreover, we report the first complete genome of a "Ca. Phytoplasma sacchari" isolate, thus contributing to future studies on the evolutionary relationships and pathogenic mechanisms of "Ca. Phytoplasma sacchari" isolates.

Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

'Candidatus Phytoplasma ziziphi' Changes the Metabolite Composition of Jujube Tree Leaves and Affects the Feeding Behavior of Its Insect Vector Hishimonus hamatus Kuoh

Hishimonus hamatus Kuoh is a leafhopper species native to China that feeds on Chinese jujube leaves. This leafhopper species has been verified to transmit jujube witches' broom (JWB) disease, caused by phytoplasma, a fatal plant pathogen, which belongs to the phytoplasma subgroup 16SrV-B. The transmission of JWB phytoplasma largely relies on the feeding behavior of piercing-sucking leafhoppers. However, the specific mechanisms behind how and why the infection of JWB influences the feeding behavior of these leafhoppers are not fully understood. To address this, a study was conducted to compare the feeding patterns of H. hamatus when feeding JWB-infested jujube leaves to healthy leaves using the electrical penetration graph (EPG) technique. Then, a widely targeted metabolome analysis was performed to identify differences in the metabolite composition of JWB-infected jujube leaves and that of healthy jujube leaves. The results of EPG analyses revealed that when feeding on JWB-infected jujube leaves, H. hamatus exhibited an increased frequency of phloem ingestion and spent longer in the phloem feeding phase compared to when feeding on healthy leaves. In addition, the results of metabolomic analyses showed that JWB-infected leaves accumulated higher levels of small-molecular carbohydrates, free amino acids, and free fatty acids, as well as lower levels of lignans, coumarins and triterpenoids compared to healthy leaves. The above results indicated that the H. hamatus preferentially fed on the phloem of infected leaves, which seems to be linked to the transmission of the JWB phytoplasma. The results of metabolomic analyses partially imply that the chemical compounds might play a role in making the infected leaves more attractive to H. hamatus for feeding.

The observation of taxonomic boundaries for the 16SrII and 16SrXXV phytoplasmas using genome-based delimitation

Within the 16SrII phytoplasma group, subgroups A-X have been classified based on restriction fragment length polymorphism of their 16S rRNA gene, and two species have been described, namely 'Candidatus Phytoplasma aurantifolia' and 'Ca. Phytoplasma australasia'. Strains of 16SrII phytoplasmas are detected across a broad geographic range within Africa, Asia, Australia, Europe and North and South America. Historically, all members of the 16SrII group share ≥97.5 % nucleotide sequence identity of their 16S rRNA gene. In this study, we used whole genome sequences to identify the species boundaries within the 16SrII group. Whole genome analyses were done using 42 phytoplasma strains classified into seven 16SrII subgroups, five 16SrII taxa without official 16Sr subgroup classifications, and one 16SrXXV-A phytoplasma strain used as an outgroup taxon. Based on phylogenomic analyses as well as whole genome average nucleotide and average amino acid identity (ANI and AAI), eight distinct 16SrII taxa equivalent to species were identified, six of which are novel descriptions. Strains within the same species had ANI and AAI values of >97 %, and shared ≥80 % of their genomic segments based on the ANI analysis. Species also had distinct biological and/or ecological features. A 16SrII subgroup often represented a distinct species, e.g., the 16SrII-B subgroup members. Members classified within the 16SrII-A, 16SrII-D, and 16SrII-V subgroups as well as strains classified as sweet potato little leaf phytoplasmas fulfilled criteria to be included as members of a single species, but with subspecies-level relationships with each other. The 16SrXXV-A taxon was also described as a novel phytoplasma species and, based on criteria used for other bacterial families, provided evidence that it could be classified as a distinct genus from the 16SrII phytoplasmas. As more phytoplasma genome sequences become available, the classification system of these bacteria can be further refined at the genus, species, and subspecies taxonomic ranks.

Potential mobile units drive the horizontal transfer of phytoplasma effector phyllogen genes

Phytoplasmas are obligate intracellular plant pathogenic bacteria that can induce phyllody, which is a type of abnormal floral organ development. Phytoplasmas possess phyllogens, which are effector proteins that cause phyllody in plants. Phylogenetic comparisons of phyllogen and 16S rRNA genes have suggested that phyllogen genes undergo horizontal transfer between phytoplasma species and strains. However, the mechanisms and evolutionary implications of this horizontal gene transfer are unclear. Here, we analyzed synteny in phyllogen flanking genomic regions from 17 phytoplasma strains that were related to six 'Candidatus' species, including three strains newly sequenced in this study. Many of the phyllogens were flanked by multicopy genes within potential mobile units (PMUs), which are putative transposable elements found in phytoplasmas. The multicopy genes exhibited two distinct patterns of synteny that correlated with the linked phyllogens. The low level of sequence identities and partial truncations found among these phyllogen flanking genes indicate that the PMU sequences are deteriorating, whereas the highly conserved sequences and functions (e.g., inducing phyllody) of the phyllogens suggest that the latter are important for phytoplasma fitness. Furthermore, although their phyllogens were similar, PMUs in strains related to 'Ca. P. asteris' were often located in different regions of the genome. These findings strongly indicate that PMUs drive the horizontal transfer of phyllogens among phytoplasma species and strains. These insights improve our understanding of how symptom-determinant genes have been shared among phytoplasmas.

Detection, Identification and Molecular Characterization of the 16SrII-V Subgroup Phytoplasma Strain Associated with Pisum sativum and Parthenium hysterophorus L

Two unrelated plant species, green pea and parthenium weed, harboring typical phytoplasma symptoms, were discovered in Yunlin, Taiwan. Green pea (Pisum sativum.) and parthenium weed (Parthenium hysterophorus L.) are both herbaceous annual plants belonging to the Fabaceae and Asteraceae families, respectively. Displayed symptoms were witches' broom, phyllody and virescence, which are typical indications of phytoplasma infection. Pleomorphic phytoplasma-like bodies were observed under the transmission electron microscope in the sieve elements of symptomatic green pea and parthenium weed. The iPhyClassifier-based virtual RFLP study demonstrated that the phytoplasma associated with the diseased plants belongs to the 16SrII-V subgroup. The disease symptoms of both plants can be explained by the identification of PHYL1 and SAP11 effectors, identical to those of peanut witches' broom phytoplasma. The phytoplasma strains identified in this study present a very close phylogenetic relationship with other 16SrII-V subgroup phytoplasma strains discovered in Taiwan. These results not only convey the local status of the 16SrII-V subgroup phytoplasma strains but also encourage attention to be given to preventing the spread of this threat before it becomes pervasive.

Genome-Informed Design of a LAMP Assay for the Specific Detection of the Strain of 'Candidatus Phytoplasma asteris' Phytoplasma Occurring in Grapevines in South Africa

Grapevine yellows is one of the most damaging phytoplasma-associated diseases worldwide. It is linked to several phytoplasma species, which can vary regionally due to phytoplasma and insect-vector diversity. Specific, rapid, and reliable detection of the grapevine yellows pathogen has an important role in phytoplasma control. The purpose of this study was to develop and validate a specific loop-mediated isothermal amplification (LAMP) assay for detection of a distinct strain of grapevine 'Candidatus Phytoplasma asteris' that is present in South Africa, through implementation of a genome-informed test design approach. Several freely available, user-friendly, web-based tools were coupled to design the specific LAMP assays. The criteria for selection of the assays were set for each step of the process, which resulted in four experimentally operative LAMP assays that targeted the ftsH/hflB gene region, specific to the aster yellows phytoplasma strain from South Africa. A real-time PCR was developed, targeting the same genetic region, to provide extensive validation of the LAMP assay. The validated molecular assays are highly specific to the targeted aster yellows phytoplasma strain from South Africa, with good sensitivity and reproducibility. We show a genome-informed molecular test design and an efficient validation approach for molecular tests if reference and sample materials are sparse and hard to obtain.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Phytoplasma Taxonomy: Nomenclature, Classification, and Identification

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 Complete Genome of the “Flavescence Dorée” Phytoplasma Reveals Characteristics of Low Genome Plasticity

Members of the genus ‘Candidatus Phytoplasma’ are obligate intracellular bacteria restricted to phloem sieve elements and are able to colonize several tissues and the hemolymph in their insect vectors. The current unfeasibility of axenic culture and the low complexity of genomic sequences are obstacles in assembling complete chromosomes. Here, a method combining pathogen DNA enrichment from infected insects and dual deep-sequencing technologies was used to obtain the complete genome of a phytoplasma causing Grapevine Flavescence dorée. The de novo assembly generated a circular chromosome of 654,223 bp containing 506 protein-coding genes. Quality assessment of the draft showed a high degree of completeness. Comparative analysis with other phytoplasmas revealed the absence of potential mobile units and a reduced amount of putative phage-derived segments, suggesting a low genome plasticity. Phylogenetic analyses identified Candidatus Phytoplasma ziziphi as the closest fully sequenced relative. The “Flavescence dorée” phytoplasma strain CH genome also encoded for several putative effector proteins potentially playing a role in pathogen virulence. The availability of this genome provides the basis for the study of the pathogenicity mechanisms and evolution of the Flavescence dorée phytoplasma.