MPMI PhytoFrontiers Phytobiomes all journals
RESEARCHOpen Access icon OPENOpen Access license

Detection and Identification of Phytoplasmas in the 16SrIV-A, -B, and -D Subgroups in Palms in Tabasco, Mexico

    Affiliations
    Authors and Affiliations
    • Eder Ramos Hernández1
    • Julia M. Lesher Gordillo2
    • Carlos Oropeza Salín3
    • Carlos F. Ortiz García4
    • Miguel A. Magaña Alejandro2
    • Saúl Sánchez Soto4
    • Yesenia García Estrada4
    1. 1Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias, Campo Experimental Huimanguillo, 86400, Huimanguillo, Tabasco, México
    2. 2Universidad Juárez Autónoma de Tabasco, División Académica de Ciencias Biológicas, 86150, Villahermosa, Tabasco, México
    3. 3Centro de Investigación Científica de Yucatán A.C., 97205, Chuburná de Hidalgo, Mérida, Yucatán, México
    4. 4Colegio de Postgraduados, 86500, Cárdenas, Tabasco, México

    Published Online:https://doi.org/10.1094/PDIS-09-18-1488-RE

    Abstract

    The 16SrIV-A phytoplasmas are associated with the devastating disease lethal yellowing (LY) of palms. In Tabasco (Mexico), the death of Cocos nucifera, Adonidia merrillii, and Attalea butyracea palms have been suspected to be associated with LY based on symptomatology. Samples from the trunk of both symptomatic and nonsymptomatic palms were collected in three different environments: two species of palms within a rural zone and the other within an urban zone. DNA was extracted to perform a nested PCR with phytoplasma primers P1/P7-LY16SF/R16R2. A 1,345-bp fragment was amplified from the DNA extracted from each of the 29 LY-symptomatic palms sampled. Phytoplasma identification was achieved by amplicon sequencing and virtual restriction fragment length polymorphism analyses. Three 16SrIV phytoplasma subgroups were detected: 16SrIV-A in C. nucifera, 16SrIV-B in A. merrillii, and 16SrIV-D in C. nucifera, A. merrillii, and A. butyracea. Phylogenetic analysis showed also that the three phytoplasma strains found in the palm species clustered with phytoplasmas reported in the literature in the three subgroups identified. This is the first report of phytoplasmas associated with these palm species in Tabasco.

    Phytoplasma diseases affect approximately 1,000 plant species belonging to 98 families (Gasparich 2010) and are distributed worldwide in 85 countries (Lee et al. 2000), although some of them may be limited to one continent or one specific region (Lee et al. 1998a). Phytoplasma diseases cause important economic losses to many crops (Bertaccini et al. 2014). These bacteria lack cell walls, colonize the plant phloem tissues, and are transmitted by insect vectors, mainly leaf-hoppers, plant-hoppers, and psyllids (Bertaccini 2007).

    The use of the polymerase chain reaction (PCR) and sequencing of the corresponding amplicons generated from the ribosomal RNA (rRNA) gene have allowed the identification and comparison of phytoplasma strains with techniques such as restriction fragment length polymorphism (RFLP) and phylogenetic analysis (Gundersen and Lee 1996; Lee et al. 1998b; Santos-Cervantes et al. 2008; Schneider et al. 1993; Wei et al. 2007). The interactive iPhyClassifier tool (https://plantpathology.ba.ars.usda.gov/cgi-bin/resource/iphyclassifier.cgi) performs a similarity analysis based on the sequence of the R16F2n/R2 fragment from the phytoplasma 16S rRNA gene (Zhao et al. 2013). This tool is used to do the preliminary classification of phytoplasmas at the group/subgroup level and/or to allocate them in ‘Candidatus Phytoplasma’ species (Nejat et al. 2013; Zhao et al. 2009a, 2013).

    In the Americas, the devastating lethal yellowing (LY) disease in coconut (Cocos nucifera L.) is mainly associated with the presence of phytoplasmas belonging to the 16SrIV-A group (Harrison et al. 2002) and is transmitted by the American palm cixiid, Haplaxius crudus (Van Duzee) (Howard et al. 1983), which is the only described vector of LY in the western hemisphere. The phytoplasmas associated with LY and LY-type diseases in palms in the Americas have been classified into different subgroups in the 16SrIV group (Fig. 1) (Ntushelo et al. 2013; Vázquez-Euán et al. 2011). Subgroups include (i) 16SrIV-A, detected in the majority of palms that have been killed by LY; (ii) 16SrIV-B (Tymon et al. 1998); (iii) 16SrIV-D (Bahder et al. 2018a, 2018b; Cordova et al. 2017; Harrison et al. 2008; Rodrigues et al. 2010) (the phytoplasma disease associated with 16SrIV-D, currently termed lethal bronzing disease, is synonymous with Texas phoenix palm decline [TPPD]) (Bahder et al. 2018b); (iv) 16SrIV-E (Martinez et al. 2008); and (v) 16SrIV-F (Harrison et al. 2008). 16SrIV-C has been detected only in Africa and is known as Tanzanian coconut lethal decline.

    Fig. 1.

    Fig. 1. Distribution map of phytoplasmas in the 16SrIV group associated with lethal yellowing diseases in the palms of North and Central America and the Caribbean Islands.

    Download as PowerPoint

    More recently, phytoplasmas of other ribosomal groups have been identified in Arecaceae in other parts of the world: ‘Candidatus Phytoplasma asteris’ (16SrI), a phytoplasma that is also reported in oil palm in Colombia (Alvarez et al. 2014); ‘Ca. P. oryzae’ (16SrXI) in India (Manimekalai et al. 2014); ‘Ca. P. cynodontis’ (16SrXIV) (Nejat et al. 2009), ‘Ca. P. malaysianum’ (16SrXXXII) (Nejat et al. 2013), and ‘Ca. P. wodyetiae’ in Malaysia (Naderali et al. 2017); ‘Ca. P. palmicola’ (16SrXXII) in Africa (Harrison et al. 2014); and ‘Ca. P. noviguineense’ in Papua New Guinea (Miyazaki et al. 2018).

    In an urban landscape in Yucatán, Mexico, where C. nucifera, Sabal mexicana Martius, Thrinax radiata Lodd. ex Schult. & Schult. f., and Pseudophoenix sargentii H. Wendl. ex Sarg. coexist, two phytoplasma subgroups,16SrIV-A and 16SrIV-D, have been found (Vázquez-Euán et al. 2011). In addition, these authors found that in a single palm species (S. mexicana or T. radiata), both phytoplasma subgroups can be found, although never within a single individual palm.

    In Tabasco, another region of Mexico, plants of different palm species, including C. nucifera, Adonidia merrillii (Becc.) Becc., and Attalea butyracea (Mutis ex L.f.) Wess. Boer, were found showing LY-type symptoms. This study reports the analyses carried out in these palms to determine if the observed syndromes are associated with phytoplasmas and classify the detected phytoplasmas.

    Materials and Methods

    Palm sampling.

    Samples were collected from both symptomatic and asymptomatic C. nucifera, A. merrillii, and A. butyracea. Symptoms included fruit drop, necrotic inflorescence, and yellowing of low and intermediate leaves. Samples were collected from palms in different environments used for various purposes (large-scale commercial agriculture, family farming, ornamental), in Tabasco State (Table 1). C. nucifera samples were collected from rural communities in Pailebot, Cárdenas, a site with monoculture, and family gardens in Miahuatlán, Cunduacán. Additionally, A. merrillii palms from urban environments were sampled from the gardens of the Academic Division of Biological Sciences, Juarez Autonomous University of Tabasco, and the Cárdenas center housing unit. A. butyracea samples were collected in grassland areas of Miahuatlán, Cunduacán.

    Table 1. Ribosomal subgroups of 16SrIV phytoplasmas identified in symptomatic plants with lethal yellowing listed by agricultural activities, palm species, environment, and sampling sites within Tabasco State (Mexico)

    DNA extraction.

    A total of 42 trunk-boring samples (3 to 6 g) were collected in 2017 from palms using a portable drill. Of each sample, 1 g was ground to a fine powder by placing it in liquid nitrogen and using a sterile mortar and pestle. Nucleic acid was extracted from the powdered tissue according to the method of Doyle and Doyle (1990). Aliquots from each final preparation were suspended in Tris-EDTA (TE) buffer (10 mM Tris, 1 mM EDTA, pH 8) and subsequently visualized in 1% agarose gels in TAE buffer (40 mM Tris-acetate, 1 mM EDTA).

    Nested PCR.

    The employed protocol was according to Harrison et al. (2002a). Each DNA sample was diluted 1:10 with sterile deionized water, and subsequently 5 µl of each dilution was used as the template for PCR amplification, including universal rRNA gene primers P1 (Deng and Hiruki 1991) and P7 (Schneider et al. 1995) for initial amplification (1,800 bp). The product of the reaction was diluted 1:20 or 1:40 with sterile deionized water, and 5 μl of each dilution was then reamplified (1,345 bp) using primers LY16SF (Harrison et al. 2002a) and R16R2 (Lee et al. 1993). PCR amplifications were performed with a C1000 Thermal Cycler (Bio-Rad, Hercules, CA) using 1 U of the Mango Taq DNA polymerase (Bioline, London, U.K.) in a 25-µl volume that contained 2 µl of dNTPs (1.25 mM), 5 µl of 10× PCR buffer, 5 U of Taq polymerase, and 2 µl of each primer (25 ng/µl). The reaction was performed with an initial 2-min cycle at 94°C, followed by 35 cycles of 1 min at 94°C, annealing for 50 s at 54°C (60°C for LY16SF/R16R2), and extension for 80 s at 72°C. Reactions were terminated by a 10-min extension step at 72°C and cooled to 4°C. The PCR products were analyzed by electrophoresis with a 1% agarose gel stained with ethidium bromide and visualized by using a UV transilluminator Molecular Image Gel Doc XR System (Bio-Rad). Negative controls were DNA extracted from nonsymptomatic plants and sterile water. Coconut palm DNA (GenBank no. GU473590) previously identified as infected by LY was used as a positive control (Vázquez-Euán et al. 2011).

    PCR product purification and sequencing.

    PCR amplicons from the LY16SF/R16R2 primers were purified in columns of the QIAquick Gel Extraction kit (cat. no. 28706; Qiagen, Hilden, Germany) according to the manufacturer’s instructions. All purified products were quantified in a nanodrop spectrophotometer (Jenway-Genova Nano, Stone, U.K.). The sequencing was performed with the 3500xl genetic analyzer (Applied Biosystems, Foster City, CA) at the Biotechnology Institute of the National Autonomous University of Mexico in both directions with primers used for the amplification plus the 503F primer (Harrison et al. 1999).

    Sequence analysis.

    Sequence editing and assembly were carried out using the BioEdit 7.2.5 software (Hall 1999). The sequences obtained in this study were compared with those from known phytoplasmas using the Basic Local Alignment Search Tool (BLAST) system of the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/). The sequences obtained were aligned using the ClustalX 2.1 software (Thompson et al. 1997). A phylogenetic tree was constructed with the phytoplasma sequences using Molecular Evolutionary Genetics Analysis (MEGA) 7.0.18 software (Kumar et al. 2016) and the neighbor-joining method, and it was evaluated with 2,000 replicates for the bootstrap analysis to estimate the stability of the inferred subclades. The 16S rRNA gene sequence from Bacillus subtilis (GenBank no. AJ544538) was used as an outgroup for the tree root.

    In silico RFLP.

    Partial sequences from the 16S rRNA gene corresponding to the R16F2n/R2 fragment (Lee et al. 1995) were analyzed by RFLP in silico using the iPhyClassifier tool.

    Results

    Phytoplasma detection in palms.

    All the samples from the 29 palms showing LY within rural and urban environments in Tabasco produced positive results for phytoplasma detection by nested PCR, whereas no detection was observed for the 13 asymptomatic palms analyzed (Table 1).

    Symptomatology.

    The symptoms manifested by the PCR-positive palms varied depending on the species. In the rural environment with large-scale commercial activity in Pailebot, seven C. nucifera palms were found with typical LY symptoms, showing fruit drop, necrotic inflorescences, and yellowing of basal leaves. In a rural environment in the Miahuatlán community in Cunduacán used for family farming activities and as ornamentals, different palm species were found with LY-type symptoms. These included six C. nucifera palms showing LY symptoms as described above; four symptomatic A. merrillii palms showing necrotic inflorescences, brown leaves with collapsed but persistent petioles (whereas younger leaves were green and vertical), and an A. butyracea palm with leaf yellowing symptoms and complete fruit drop. Finally, eight A. merrillii palms in a city garden in the urban environment and another three A. merrillii palms in a Cárdenas housing unit showed the same symptoms as described above for A. merrillii palms in Cunduacán.

    All C. nucifera and A. merrillii palms that were symptomatic and positive for phytoplasma presence in nested-PCR detection died during the study, whereas the positive A. butyracea palm showed only a slow decline over 49 months and has not died yet (Fig. 2).

    Fig. 2.

    Fig. 2. Attalea butyracea showed only a slow decline over 49 months: a, yellowing of older leaves, July 2014; b, continuous yellowing in intermediate leaflets, December 2014; c, yellowing on all leaves, June 2015; and d, foliage recovery and death of inflorescence, July 2018.

    Download as PowerPoint

    In silico RFLP analysis and phytoplasma identification.

    The ribosomal DNA was sequenced from 20 out of the 29 samples from positive palms, and 16 sequences could be obtained (GenBank accession nos. MG881770 to MG881785). These sequences were subjected to virtual RFLP analysis of the R16F2n/R2 region in the phytoplasma ribosomal gene, and the results showed that phytoplasmas enclosed in three different subgroups were present in the samples: 16SrIV-A, -B, and -D, with profiles identical to phytoplasmas of LY (GenBank no. U18747), Yucatan coconut lethal decline (GenBank no. U18753), and Sabal palmetto decline (GenBank no. HQ613895), respectively (Fig. 3).

    Fig. 3.

    Fig. 3. Distinctive restriction fragment length polymorphism patterns obtained with iPhyClassifier from in silico digestion of the 16S rRNA gene (Zhao et al. 2009b) to determine the subgroup in the 16SrIV group. GenBank accession numbers of reference sequences used were as follows: 16SrIV-A (AF498307), 16SrIV-B (AF500334), 16SrIV-C (X80117), and 16SrIV-D (AF237615). The separation of the 16SrIV-B and -D subgroups is achieved by comparing the AluI profiles. The separation of 16SrIV-A and -C is obtained by comparing the MseI band patterns. Restriction fragment polymorphisms are shown on a virtual 3% agarose gel. Lanes labeled MW represent HaeIII-digested phage ΦX174 DNA.

    Download as PowerPoint

    The virtual RFLP patterns obtained with AluI digestion clearly differentiated 16SrIV-B and 16SrIV-D subgroups between themselves, whereas those of 16SrIV-A and 16SrIV-C had the same pattern with AluI (Fig. 3) but were different with MseI from the patterns of 16SrIV-A and 16SrIV-C (Fig. 3).

    Phylogenetic analysis.

    The phylogenetic analysis on the 16SrRNA gene (Fig. 4) confirmed that all phytoplasmas from this study clustered within phytoplasmas in the 16SrIV clade and were separated into three subclades, corresponding to three ribosomal subgroups. The 16SrIV-A sequences were in the same subclade as the sequences of the S. mexicana (GenBank no. GU473589, Yucatán) and C. nucifera (GenBank no. AF498307, Jamaica) phytoplasmas. This phytoplasma subgroup was detected in diseased monocultures of C. nucifera palms within the agricultural area and family gardens located in the farming area. The 16SrIV-D sequences from A. merrillii, A. butyracea, and C. nucifera were in the same subclade as Texas Phoenix decline (GenBank no. AF434989). These plants were used for farming purposes and as ornamentals in Miahuatlán, Cunduacán, and in the case of A. merrillii they were located within an urban area in the Cárdenas housing unit. The strains of the third phytoplasma subgroup, 16SrIV-B in A. merrillii palms, clustered in the same subclade as Yucatan coconut lethal decline (GenBank nos. U18753 and AF500334). These plants were found in the central gardens at the Academic Division of Biological Sciences, Juarez Autonomous University of Tabasco.

    Fig. 4.

    Fig. 4. Phylogenetic tree derived from selected phytoplasma 16S rRNA gene sequences. The tree was constructed using the neighbor-joining method using the Tajima–Nei substitution model. Bacillus subtilis (GenBank no. AJ544538) was used as an outgroup to root the tree. The phytoplasma strain sequences obtained in this study are bolded and are grouped as indicated by circles, squares, and triangles. The branch length represents the distance between the sequences. The bootstrap values for 2,000 replicates are shown on the branches. The GenBank accession number for each sequence is indicated in parentheses followed by the phytoplasma subgroup.

    Download as PowerPoint

    Discussion

    Phytoplasmas in the 16SrIV group are associated with lethal diseases of the coconut palm with high economic importance, as well as other ornamental palm species of commercial value, in North and Central America and the Caribbean Islands. The phytoplasmas in the 16SrIV subgroups have been associated with LY and LY-type diseases with symptoms similar to those of LY, including early fruit drop, inflorescence necrosis, leaf discoloration, and decline (Myrie et al. 2014). In the present study, three 16SrIV phytoplasma subgroups are reported: 16SrIV-A detected in C. nucifera, 16SrIV-B in A. merrillii, and 16SrIV-D in C. nucifera, A. merrillii, and A. butyracea. All these cases were detected in Tabasco State in Mexico.

    The subgroup 16SrIV-A of phytoplasmas is associated with LY in C. nucifera and is the prevalent phytoplasma subgroup detected so far in several countries in the Caribbean region, including Mexico, with several reports confirming the presence of 16SrIV-A phytoplasmas in the states of Quintana Roo, Campeche, and Yucatán (Harrison et al. 2002b; Oropeza et al. 2011; Vázquez-Euán et al. 2011). The occurrence of 16SrIV-A phytoplasmas has also been described in other palms, mostly used for ornamental purposes: for example, A. merrillii, Phoenix canariensis H. Wildpret, and P. dactylifera L. in Florida (Tymon et al. 1998) and Acrocomia mexicana Karw. ex Mart., Coccothrinax readii H.J. Quero, S. mexicana, T. radiata, and Roystonea regia (Kunth) O.F. Cook in Yucatán (Cordova et al. 2017; Narvaez et al. 2006; Vázquez-Euán et al. 2011). The current finding of this subgroup in C. nucifera in Tabasco extends the knowledge of its distribution in Mexico, suggesting that this subgroup is certainly the most widely distributed in Mexico and the Americas. This can be attributed to the ubiquitous distribution of the highly susceptible Atlantic Tall variety.

    The 16SrIV-B phytoplasmas detected in A. merrillii palms were originally found in Tabasco in a C. nucifera palm showing atypical LY symptoms (Tymon et al. 1998). Subsequently, Harrison et al. (2002b) reported their presence in C. nucifera on the Pacific Coast in Mexico in Guerrero and Oaxaca. Later, a phytoplasma of the same subgroup was detected in Honduras in declining Acrocomia aculeata (Jacq.) Lodd. et al. ex Mart. palms (Roca et al. 2006). Other phytoplasma subgroups have been reported in A. merrillii. Recently, Cordova et al. (2017) reported 16SrIV-A and 16SrIV-D phytoplasmas occurring separately in different individuals of A. merrillii with LY symptoms in urban landscapes in Mérida in Yucatán city. Therefore, this is the first report of the occurrence of a 16SrIV-B subgroup in A. merrillii. This phytoplasma subgroup seems, therefore, to have a wide geographical distribution, infecting several palm species growing in various ecological conditions.

    Another phytoplasma that was also detected in the present study is 16SrIV-D. It was found in C. nucifera, A. merrillii, and A. butyracea palms that were planted for different purposes, in rural areas in Miahuatlán in Tabasco. This phytoplasma subgroup had already been reported in Mexico, in C. nucifera palms in the states of Guerrero and Oaxaca on the Pacific Coast of Mexico (Harrison et al. 2002b) and in A. merrillii, P. sargentii, S. mexicana, T. radiata, and Pritchardia pacifica Seem. & H. Wendl. in urban landscapes of Yucatán (Cordova et al. 2017; Narváez et al. 2017; Vázquez-Euán et al. 2011), and also in a nonpalm species, Carludovica palmata Ruiz & Pav., in Campeche State (Cordova et al. 2000).

    16SrIV-D has also been detected in different states and territories of the United States. In Florida, it has been found in Butiagrus nabonnandii (Prosch.) Vorster, P. canariensis, P. roebelenii O’Brien, P. dactylifera, P. sylvestris (L.) Roxb., P. reclinata Jacq., S. palmetto (Walter) Lodd. ex Schult. & Schult. f., Syagrus romanzoffiana (Cham.) Glassman, P. sylvestris × P. dactylifera Chab. hybrid, and P. sylvestris × P. canariensis hybrid Chab. (Bahder et al. 2018a, 2018b; Harrison et al. 2008); in Louisiana, it was found in Trachycarpus fortunei (Hook.) H. Wendl. (Singh and Ferguson 2017); in Texas, it was found in different P. canariensis and P. dactylifera species affected with TPPD (Harrison et al. 2002a; Ntushelo et al. 2013); in Puerto Rico, this subgroup was reported to affect the palms of the landscape and tourism industry, where the affected species were Roystonea sp. O.F. Cook, Carpentaria acuminata (H. Wendl. & Drude) Becc., and Caryota mitis Lour. (Rodrigues et al. 2010).

    The three species of host palms reported in this study coexist in Miahuatlán with at least two subgroups of phytoplasmas that could have favored the occurrence of mixed infections as reported by Harrison et al. (2008), in P. dactylifera with 16SrIV-A and 16SrIV-F; however, in this study, no mixed infection was detected.

    In conclusion, the present report is the first detecting the presence of phytoplasmas 16SrIV-A, 16SrIV-B, 16SrIV-D in palm plants in Tabasco, and also the first report of 16SrIV-B in A. merrillii and 16SrIV-D in A. butyracea. These results extend the knowledge of the host diversity of these phytoplasmas and widen their geographic distribution. There is a need for continuing research efforts on the diversity of phytoplasmas in different host plant species, and to study their geographic distribution together with the identification of other insect vectors.

    The author(s) declare no conflict of interest.

    Literature Cited

    The author(s) declare no conflict of interest.

    Funding: The authors thank the Colegio de Postgraduados (COLPOS) for providing the infrastructure needed to carry out a part of this research. The European Union provided funds under the Horizon 2020 project: Insect borne prokaryote-associated diseases in tropical and subtropical perennial crops (TROPICSAFE: Grant Agreement No. 727459).