Should 'Candidatus Phytoplasma' be retained within the order Acholeplasmatales?

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.

Guilt by Association: DNA Barcoding-Based Identification of Potential Plant Hosts of Phytoplasmas from Their Insect Carriers

Phytoplasmas are small phloem-restricted and insect-transmissible bacteria that infect many plant species, including important crops and ornamental plants, causing severe economic losses. Our previous studies screened phytoplasmas in hundreds of leafhoppers collected from natural habitats worldwide and identified multiple genetically different phytoplasmas in seven leafhopper species (potential insect vectors). As an initial step toward determining the impact of these phytoplasmas on the ecosystem, ribulose 1,5-biphosphate carboxylase large subunit (rbcL), a commonly used plant DNA barcoding marker, was employed to identify the plant species that the phytoplasma-harboring leafhoppers feed on. The DNA of 17 individual leafhoppers was PCR amplified using universal rbcL primers. PCR products were cloned, and five clones per amplicon were randomly chosen for Sanger sequencing. Moreover, Illumina high-throughput sequencing on selected PCR products was conducted and confirmed no missing targets in Sanger sequencing. The nucleotide BLAST results revealed 14 plant species, including six well-known plant hosts of phytoplasmas such as tomato, alfalfa, and maize. The remaining species have not been documented as phytoplasma hosts, expanding our knowledge of potential plant hosts. Notably, the DNA of tomato and maize (apparently cultivated in well-managed croplands) was detected in some phytoplasma-harboring leafhopper species sampled in non-crop lands, suggesting the spillover/spillback risk of phytoplasma strains between crop and non-crop areas. Furthermore, our results indicate that barcoding (or metabarcoding) is a valuable tool to study the three-way interactions among phytoplasmas, plant hosts, and vectors. The findings contribute to a better understanding of phytoplasma host range, host shift, and disease epidemiology.

Structure and variation of root-associated bacterial communities of Cyperus rotundus L. in the contaminated soils around Pb/Zn mine sites

Soil contamination due to mining activities is a great concern in China. Although the effects of mining pollution resulting in changes of soil characteristics and the microbiome have been documented, studies on the responses of plant root-associated microbial assemblages remain scarce. In this work, we collected bulk soil, rhizosphere soil, and root endosphere samples of Cyperus rotundus L (Cyp) plants from two Pb/Zn mines, of which, one was abandoned (SL) and the other was active (GD), to investigate the bacterial community responses across different site contamination levels and Cyp plant compartments. For comparison, one unpolluted site (SD) was included. Results revealed that soils from the SL and GD sites were seriously contaminated by metal(loid)s, including Pb, Zn, As, and Sb. Bacterial richness and diversity depended on the sampling site and plant compartment. All sample types from the SL site had the lowest bacterial diversities and their bacterial communities also exhibited distinct patterns compared to GD and SD samples. As for the specific sampling site, bacterial communities from the root endosphere exhibited different patterns from those in bulk and rhizosphere soil. Compared to the GD and SD sites, the root endosphere and the rhizosphere soil from the SL site shared core microbes, including Halomonas, Pelagibacterium, and Chelativorans, suggesting that they play key roles in Cyp plant survival in such harsh environments.

Screening potential insect vectors in a museum biorepository reveals undiscovered diversity of plant pathogens in natural areas

Phytoplasmas (Mollicutes, Acholeplasmataceae), vector-borne obligate bacterial plant parasites, infect nearly 1,000 plant species and unknown numbers of insects, mainly leafhoppers (Hemiptera, Deltocephalinae), which play a key role in transmission and epidemiology. Although the plant-phytoplasma-insect association has been evolving for >300 million years, nearly all known phytoplasmas have been discovered as a result of the damage inflicted by phytoplasma diseases on crops. Few efforts have been made to study phytoplasmas occurring in noneconomically important plants in natural habitats. In this study, a subsample of leafhopper specimens preserved in a large museum biorepository was analyzed to unveil potential new associations. PCR screening for phytoplasmas performed on 227 phloem-feeding leafhoppers collected worldwide from natural habitats revealed the presence of 6 different previously unknown phytoplasma strains. This indicates that museum collections of herbivorous insects represent a rich and largely untapped resource for discovery of new plant pathogens, that natural areas worldwide harbor a diverse but largely undiscovered diversity of phytoplasmas and potential insect vectors, and that independent epidemiological cycles occur in such habitats, posing a potential threat of disease spillover into agricultural systems. Larger-scale future investigations will contribute to a better understanding of phytoplasma genetic diversity, insect host range, and insect-borne phytoplasma transmission and provide an early warning for the emergence of new phytoplasma diseases across global agroecosystems.

Rubidium chloride modulated the fecal microbiota community in mice

Background

The microbiota plays an important role in host health. Although rubidium (Rb) has been used to study its effects on depression and cancers, the interaction between microbial commensals and Rb is still unexplored. To gain the knowledge of the relationship between Rb and microbes, 51 mice receiving RbCl-based treatment and 13 untreated mice were evaluated for their characteristics and bacterial microbiome changes.

Results

The 16S ribosomal RNA gene sequencing of fecal microbiota showed that RbCl generally maintained fecal microbial community diversity, while the shifts in fecal microbial composition were apparent after RbCl exposure. RbCl significantly enhanced the abundances of Rikenellaceae, Alistipes, Clostridium XlVa and sulfate-reducing bacteria including Deltaproteobacteria, Desulfovibrionales, Desulfovibrionaceae and Desulfovibrio, but significantly inhibited the abundances of Tenericutes, Mollicutes, Anaeroplasmatales, Anaeroplasmataceae and Anaeroplasma lineages. With regarding to the archaea, we only observed two less richness archaea Sulfolobus and Acidiplasma at the genus level.

Conclusions

Changes of fecal microbes may in part contribute to the anticancer or anti-depressant effects of RbCl. These findings further validate that the microbiome could be a target for therapeutic intervention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02095-4.

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'.

Draft Genome Sequence of "Candidatus Phytoplasma pruni" (X-Disease Group, Subgroup 16SrIII-B) Strain ChTDIII from Argentina

Herein, we report the draft genome sequence of "Candidatus Phytoplasma pruni" strain ChTDIII (subgroup 16SrIII-B). The final assembly consists of 790,517 nucleotides organized in 67 contigs (minimal size, 1 kb), with a G+C content of 29.4% and encoding 672 proteins.

Phytoplasmas diversity and identification of new aster yellows subgroup (16SrI) associated with weed species in Argentina

Symptoms of phytoplasma infection were observed in different weed species, Bidens subalternans, Conyza bonariensis, Heterosperma ovatifolium and Conium maculatum, collected from diverse geographical regions in Argentina. To confirm the association of phytoplasma infection with symptomatic plants, PCR, RFLP and phylogenetic analyses based on 16S rRNA-encoding sequences were performed. In this work, we report the presence of phytoplasmas from group 16SrVII (subgroup 16VII-B) infecting C. bonariensis and B. subalternans and from group 16SrIII (subgroup 16SrIII-X) B. subalternans, H. ovatifolium, and C. maculatum. Phytoplasmas from the aster yellows group were detected infecting C. bonariensis and B. subalternans. Analysis of 16S rRNA-encoding genes revealed the presence of two distinct operons, rrnB (16SrI-B) and newly described rrnA, which is different from the reference RFLP patterns of all previously established 16SrI-subgroups. A single rp operon sequence analysis reveals the presence of simple infection and confirms a description of a novel subgroup. On the basis of these results we propose a designation of new subgroup 16SrI-(B/AJ) AJ (rp-AJ). To our knowledge, this is the first report of phytoplasmas infecting Bidens subalternans¸ Heterosperma ovatifolium and Conium maculatum.

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.

Genotyping Points to Divergent Evolution of 'Candidatus Phytoplasma asteris' Strains Causing North American Grapevine Yellows and Strains Causing Aster Yellows

Grapevine yellows diseases occur in cultivated grapevine (Vitis vinifera L.) on several continents, where the diseases are known by different names depending upon the identities of the causal phytoplasmas. In this study, phytoplasma strains associated with grapevine yellows disease (North American grapevine yellows [NAGY]) in vineyards of Pennsylvania were characterized as belonging to 16S ribosomal RNA (rRNA) gene restriction fragment length polymorphism group 16SrI (aster yellows phytoplasma group), subgroup 16SrI-B (I-B), and variant subgroup I-B*. The strains (NAGYI strains) were subjected to genotyping based on analyses of 16S rRNA and secY genes, and to in silico three-dimensional modeling of the SecY protein. Although the NAGYI strains are closely related to aster yellows (AY) phytoplasma strains and are classified like AY strains in subgroup I-B or in variant subgroup I-B*, the results from genotyping and protein modeling may signal ongoing evolutionary divergence of NAGYI strains from related strains in subgroup 16SrI-B.

Detection and Typing of "Candidatus Phytoplasma " spp. in Host DNA Extracts Using Oligonucleotide-Coupled Fluorescent Microspheres

The use of oligonucleotide-coupled fluorescent microspheres is a rapid, sequencing-independent, and reliable way to diagnose bacterial diseases. Previously described applications of oligonucleotide-coupled fluorescent microspheres for the detection and identification of bacteria in human clinical samples have been successfully adapted to detect and differentiate "Ca. Phytoplasma" species using as a target the chaperonin 60-encoding gene. In this chapter, we describe in detail the design and validation of oligonucleotide capture probes, and their application in the assay aiming to differentiate phytoplasma strains infecting Brassica napus and Camelina sativa plants grown in the same geographic location at the same time.

'Candidatus Phytoplasma wodyetiae', a new taxon associated with yellow decline disease of foxtail palm (Wodyetia bifurcata) in Malaysia

Landscape-grown foxtail palm (Wodyetia bifurcata A. K. Irvine) trees displaying symptoms of severe foliar chlorosis, stunting, general decline and mortality reminiscent of coconut yellow decline disease were observed in Bangi, Malaysia, during 2012. DNA samples from foliage tissues of 15 symptomatic palms were analysed by employing a nested PCR assay primed by phytoplasma universal ribosomal RNA operon primer pairs, P1/P7 followed by R16F2n/R2. The assay yielded amplicons of a single band of 1.25 kb from DNA samples of 11 symptomatic palms. Results from cloning and sequence analysis of the PCR-amplified 16S rRNA gene segments revealed that, in three palms, three mutually distinct phytoplasmas comprising strains related to 'Candidatus Phytoplasma asteris' and 'Candidatus Phytoplasma cynodontis', as well as a novel phytoplasma, were present as triple infections. The 16S rRNA gene sequence derived from the novel phytoplasma shared less than 96 % nucleotide sequence identity with that of each previously describedspecies of the provisional genus 'Ca. Phytoplasma', justifying its recognition as the reference strain of a new taxon, 'Candidatus Phytoplasma wodyetiae'. Virtual RFLP profiles of the R16F2n/R2 portion of the 16S rRNA gene and the pattern similarity coefficient value (0.74) supported the delineation of 'Ca. Phytoplasma wodyetiae' as the sole representative subgroup A member of a new phytoplasma ribosomal group, 16SrXXXVI.

Endophytic root bacteria associated with the natural vegetation growing at the hydrocarbon-contaminated Bitumount Provincial Historic site

The Bitumount Provincial Historic site is the location of 2 of the world’s first oil-extracting and -refining operations. Despite hydrocarbon levels ranging from 330 to 24 700 mg·(kg soil)−1, plants have been able to recolonize the site through means of natural revegetation. This study was designed to achieve a better understanding of the plant-root-associated bacterial partnerships occurring within naturally revegetated hydrocarbon-contaminated soils. Root endophytic bacterial communities were characterized from representative plant species throughout the site by both high-throughput sequencing and culturing techniques. Population abundance of rhizosphere and root endosphere bacteria was significantly influenced (p < 0.05) by plant species and sampling location. In general, members of the Actinomycetales, Rhizobiales, Pseudomonadales, Burkholderiales, and Sphingomonadales orders were the most commonly identified orders. Community structure of root-associated bacteria was influenced by both plant species and sampling location. Quantitative real-time polymerase chain reaction was used to determine the potential functional diversity of the root endophytic bacteria. The gene copy numbers of 16S rRNA and 2 hydrocarbon-degrading genes (CYP153 and alkB) were significantly affected (p < 0.05) by the interaction of plant species and sampling location. Our findings suggest that some of the bacterial communities detected are known to exhibit plant growth promotion characteristics.

Identification of putative effector genes and their transcripts in three strains related to 'Candidatus Phytoplasma aurantifolia'

Molecular mechanisms underlying phytoplasma interactions with host plants are largely unknown. In this study attempts were made to identify effectors of three phytoplasma strains related to 'Ca. P. aurantifolia', crotalaria phyllody (CrP), faba bean phyllody (FBP), and witches' broom disease of lime (WBDL), using information from draft genome of peanut witches' broom phytoplasma. Seven putative effectors were identified in WBDL genome (SAP11, SAP21, Eff64, Eff115, Eff197, Eff211 and EffSAP67), five (SAP11, SAP21, Eff64, Eff99 and Eff197) in CrP and two (SAP11, Eff64) in FBP. No homologs to Eff64, Eff197 and Eff211 in phytoplasmas of other phylogenetic groups were found. SAP11 and Eff64 homologs of 'Ca. P. aurantifolia' strains shared at least 95.9% identity and were detected in the three phytoplasmas, supporting their role within the group. Five of the putative effectors (SAP11, SAP21, Eff64, Eff115, and Eff99) were transcribed from total RNA extracts of periwinkle plants infected with these phytoplasmas. Transcription profiles of selected putative effectors of CrP, FBP and WBDL indicated that SAP11 transcripts were the most abundant in the three phytoplasmas. SAP21 transcript levels were comparable to those of SAP11 for CrP and not measurable for the other phytoplasmas. Eff64 had the lowest transcription level irrespective of sampling date and phytoplasma isolate. Eff115 transcript levels were the highest in WBDL infected plants. This work reports the first sequence information for 14 putative effectors in three strains related to 'Ca. P. aurantifolia', and offers novel insight into the transcription profile of five of them during infection of periwinkle.

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.

Metabolic Consequences of Infection of Grapevine (Vitis vinifera L.) cv. "Modra frankinja" with Flavescence Dorée Phytoplasma

Flavescence dorée, caused by the quarantine phytoplasma FDp, represents the most devastating of the grapevine yellows diseases in Europe. In an integrated study we have explored the FDp-grapevine interaction in infected grapevines of cv. "Modra frankinja" under natural conditions in the vineyard. In FDp-infected leaf vein-enriched tissues, the seasonal transcriptional profiles of 14 genes selected from various metabolic pathways showed an FDp-specific plant response compared to other grapevine yellows and uncovered a new association of the SWEET17a vacuolar transporter of fructose with pathogens. Non-targeted metabolome analysis from leaf vein-enriched tissues identified 22 significantly changed compounds with increased levels during infection. Several metabolites corroborated the gene expression study. Detailed investigation of the dynamics of carbohydrate metabolism revealed significant accumulation of sucrose and starch in the mesophyll of FDp-infected leaves, as well as significant up-regulation of genes involved in their biosynthesis. In addition, infected leaves had high activities of ADP-glucose pyrophosphorylase and, more significantly, sucrose synthase. The data support the conclusion that FDp infection inhibits phloem transport, resulting in accumulation of carbohydrates and secondary metabolites that provoke a source-sink transition and defense response status.

Role of the major antigenic membrane protein in phytoplasma transmission by two insect vector species

Background

Phytoplasmas are bacterial plant pathogens (class Mollicutes), transmitted by phloem feeding leafhoppers, planthoppers and psyllids in a persistent/propagative manner. Transmission of phytoplasmas is under the control of behavioral, environmental and geographical factors, but molecular interactions between membrane proteins of phytoplasma and vectors may also be involved. The aim of the work was to provide experimental evidence that in vivo interaction between phytoplasma antigenic membrane protein (Amp) and vector proteins has a role in the transmission process. In doing so, we also investigated the topology of the interaction at the gut epithelium and at the salivary glands, the two barriers encountered by the phytoplasma during vector colonization.

Methods

Experiments were performed on the ‘Candidatus Phytoplasma asteris’ chrysanthemum yellows strain (CYP), and the two leafhopper vectors Macrosteles quadripunctulatus Kirschbaum and Euscelidius variegatus Kirschbaum.

To specifically address the interaction of CYP Amp at the gut epithelium barrier, insects were artificially fed with media containing either the recombinant phytoplasma protein Amp, or the antibody (A416) or both, and transmission, acquisition and inoculation efficiencies were measured.

An abdominal microinjection protocol was employed to specifically address the interaction of CYP Amp at the salivary gland barrier. Phytoplasma suspension was added with Amp or A416 or both, injected into healthy E. variegatus adults and then infection and inoculation efficiencies were measured.

An internalization assay was developed, consisting of dissected salivary glands from healthy E. variegatus exposed to phytoplasma suspension alone or together with A416 antibody. The organs were then either observed in confocal microscopy or subjected to DNA extraction and phytoplasma quantification by qPCR, to visualize and quantify possible differences among treatments in localization/presence/number of CYP cells.

Results

Artificial feeding and abdominal microinjection protocols were developed to address the two barriers separately. The in vivo interactions between Amp of ‘Candidatus Phytoplasma asteris’ Chrysanthemum yellows strain (CYP) and vector proteins were studied by evaluating their effects on phytoplasma transmission by Euscelidius variegatus and Macrosteles quadripunctulatus leafhoppers. An internalization assay was developed, consisting of dissected salivary glands from healthy E. variegatus exposed to phytoplasma suspension alone or together with anti-Amp antibody. To visualize possible differences among treatments in localization/presence of CYP cells, the organs were observed in confocal microscopy. Pre-feeding of E. variegatus and M. quadripunctulatus on anti-Amp antibody resulted in a significant decrease of acquisition efficiencies in both species. Inoculation efficiency of microinjected E. variegatus with CYP suspension and anti-Amp antibody was significantly reduced compared to that of the control with phytoplasma suspension only. The possibility that this was due to reduced infection efficiency or antibody-mediated inhibition of phytoplasma multiplication was ruled out. These results provided the first indirect proof of the role of Amp in the transmission process.

Conclusion

Protocols were developed to assess the in vivo role of the phytoplasma native major antigenic membrane protein in two phases of the vector transmission process: movement through the midgut epithelium and colonization of the salivary glands. These methods will be useful also to characterize other phytoplasma-vector combinations. Results indicated for the first time that native CYP Amp is involved in vivo in specific crossing of the gut epithelium and salivary gland colonization during early phases of vector infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0522-5) contains supplementary material, which is available to authorized users.

Transcriptomics-based analysis using RNA-Seq of the coconut (Cocos nucifera) leaf in response to yellow decline phytoplasma infection

Invasive phytoplasmas wreak havoc on coconut palms worldwide, leading to high loss of income, food insecurity and extreme poverty of farmers in producing countries. Phytoplasmas as strictly biotrophic insect-transmitted bacterial pathogens instigate distinct changes in developmental processes and defence responses of the infected plants and manipulate plants to their own advantage; however, little is known about the cellular and molecular mechanisms underlying host-phytoplasma interactions. Further, phytoplasma-mediated transcriptional alterations in coconut palm genes have not yet been identified. This study evaluated the whole transcriptome profiles of naturally infected leaves of Cocos nucifera ecotype Malayan Red Dwarf in response to yellow decline phytoplasma from group 16SrXIV, using RNA-Seq technique. Transcriptomics-based analysis reported here identified genes involved in coconut innate immunity. The number of down-regulated genes in response to phytoplasma infection exceeded the number of genes up-regulated. Of the 39,873 differentially expressed unigenes, 21,860 unigenes were suppressed and 18,013 were induced following infection. Comparative analysis revealed that genes associated with defence signalling against biotic stimuli were significantly overexpressed in phytoplasma-infected leaves versus healthy coconut leaves. Genes involving cell rescue and defence, cellular transport, oxidative stress, hormone stimulus and metabolism, photosynthesis reduction, transcription and biosynthesis of secondary metabolites were differentially represented. Our transcriptome analysis unveiled a core set of genes associated with defence of coconut in response to phytoplasma attack, although several novel defence response candidate genes with unknown function have also been identified. This study constitutes valuable sequence resource for uncovering the resistance genes and/or susceptibility genes which can be used as genetic tools in disease resistance breeding.