Introducción

El psílido asiático de los cítricos (PAC) Diaphorina citri Kuwayama (Hemiptera: Liviidae), es el vector de las bacterias gram-negativas ‘Candidatus Liberibacter asiaticus’, ‘Ca. L. americanus’, ‘Ca. L. africanus’ y ‘Ca. L. caribbeanus’, esta última recientemente encontrada en Colombia (^{Keremane et al., 2015}). Estas bacterias son agentes causales de la enfermedad conocida como Huanglongbing (HLB) de los cítricos que tiene un efecto devastador a escala mundial (^{Bové, 2006}; ^{Hall, Richardson et al., 2012}). Considerando el riesgo de la enfermedad y su vector, se recomienda como estrategias controlar las poblaciones del PAC con insecticidas sintéticos, eliminar árboles enfermos y realizar siembra de plantas libres de enfermedades (^{Brlansky et al., 2006}). El control químico constituye una importante herramienta para reducir la población de D. citri; sin embargo, las aplicaciones foliares dirigidas al control de adultos y ninfas, proporcionan un periodo corto de protección y tienen un efecto nocivo sobre la entomofauna benéfica (^{Qureshi et al., 2009}). Para limitar la dispersión del HLB y disminuir la densidad poblacional de D. citri se requiere de un esquema de manejo integrado, ya que el uso continuo de insecticidas químicos puede originar la selección de psílidos resistentes (^{Hoy, 1998}). Lo anterior plantea buscar y desarrollar alternativas como el control biológico, esta estrategia aspira, a la conservación de enemigos naturales de las plagas en los sistemas agrícolas. Los agentes de control biológico de D. citri reportados son varios artrópodos depredadores y parasitoides (Hall, Richardson et al., 2012), además de la infección natural de diferentes especies de hongos entomopatógenos (HE) como Hirsutella citriformis Speare (^{Étienne et al., 2001}; ^{Rivero-Aragón y Grillo-Ravelo, 2000}; ^{Subandiyah et al., 2000}); Isaria fumosorosea Wize (≡ Paecilomyces fumosoroseus) (^{Samson, 1974}; ^{Subandiyah et al., 2000}); Isaria javanica (Friedrichs y Bally) Samson y Hywel-Jones (≡ Paecilomyces javanicus) (^{Yang et al., 2006}); Lecanicillium lecanii (Zimmerman) Zare y W. Gams (≡ Verticillium lecanii) (^{Rivero-Aragón y Grillo-Ravelo, 2000}; ^{Xie et al., 1988}; ^{Yang et al., 2006}) y Beauveria bassiana (Balsamo) Vuillemin (^{Rivero-Aragón y Grillo-Ravelo, 2000}; ^{Yang et al., 2006}).

Las alternativas a los plaguicidas sintéticos como el control biológico son relevantes y pueden ser aprovechadas para llevar a cabo un manejo integrado de plagas (MIP) o cualquier tipo de acción fitosanitaria. En México, el Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria (Senasica) a través de la Dirección General de Sanidad Vegetal (DGSV), desarrolló un programa de control biológico para el manejo de D. citri a través del uso de HE y del parasitoide Tamarixia radiata (Waterston) (^{Sánchez-González et al., 2015}). Con el propósito de contar con una mayor cantidad de enemigos naturales de esta plaga, el presente trabajo tuvo como objetivo realizar exploraciones para conocer las especies de HE asociados a D. citri en huertos comerciales de cítricos y en plantas de Murraya paniculata (L) Jack en el estado de Colima, México.

Materiales y métodos

Las exploraciones para la recolección de los insectos se realizaron en huertos comerciales de limón mexicano (Citrus aurantifolia Christm. Swingle), limón persa (Citrus latifolia Tanaka) y plantas de mirto o limonaria (M. paniculata) del área urbana en los municipios de Tecomán, Armería y Manzanillo, en Colima, México, desde el año 2012 hasta 2014. El muestreo fue dirigido a los diferentes estadios del insecto. Para recolectar huevos y ninfas, se tomaron brotes tiernos de la planta, los cuales fueron colocados en bolsas de papel. En insectos adultos se buscaron signos y síntomas de infección por hongos, incluyendo alteraciones en la conducta, cambios de coloración, momificación o presencia de micelio y conidios en la superficie de su cuerpo. Los insectos con estas características se recolectaron en forma individual tomando la hoja o parte del sustrato en que se encontraban y se depositaron individualmente en contenedores secos, evitando las partículas de tierra y otros contaminantes. También se incluyó la recolecta de insectos muertos sin indicios aparentes de infección; el material fue colocado en viales entomológicos debidamente etiquetados y depositados en una hielera de poliestireno con geles refrigerantes y llevados al laboratorio para su análisis.

Los insectos recolectados en campo se examinaron en un microscopio estereoscópico (Stemi 2000-C, Carl Zeiss), para visualizar los signos externos característicos de las especies de HE. Los insectos muertos carentes de signos de infección fueron colocados en cámara húmeda (K₂SO₄, al 11%) a 25 °C para favorecer el desarrollo de posibles HE (^{Goettel e Inglis, 1997}); una vez que el hongo se desarrolló, se procedió a revisar las estructuras morfológicas. Los especímenes con signos de micosis se examinaron mediante la preparación de laminillas con lactofenol azul de algodón (10:1) o aceto orceína según el HE. Las preparaciones fueron observadas en un microscopio óptico (AXIO SCOPE A1, Carl Zeiss) para analizar las características morfológicas. Se consideraron las estructuras de reproducción asexual (conidios y conidióforos) para su identificación taxonómica, tomando como referencia las claves reportadas por ^{Humber (2012)}. Para realizar la identificación específica se procedió a evaluar su morfometría; para ello se realizaron preparaciones de insectos micosados y microcultivos de acuerdo con ^{Ayala-Zermeño et al. (2015)}. Se tomaron fotografías con el objetivo 100× (AXIO SCOPE A1) y una cámara digital AxioCam ICc 1 (Carl Zeiss) utilizando el software AxionVision 4.9.1 (03-2013). Por cada uno de los aislados se midieron 30 conidios y 30 fiálides considerando largo y ancho de cada estructura.

Para purificar los aislados a partir del insecto, se seleccionaron especímenes con micosis utilizando la transferencia directa de esporas y el método de descarga de conidios acorde a la especie de hongo causante de la infección (^{Hayek et al., 2012}). Se utilizaron los medios de cultivo Sabouraud dextrosa agar (SDA) suplementado con extracto de levadura [g L^-1: dextrosa 40, peptona de carne 5, peptona de caseína 5, extracto de levadura 10, agar 15] y SDA suplementado con yema de huevo (120 mL) y leche (80 mL); la incubación se realizó a 25 ± 2 °C hasta su esporulación. Los hongos aislados fueron depositados en la Colección de Hongos Entomopatógenos (CHE) del Centro Nacional de Referencia de Control Biológico (CNRCB-DGSV), con el acrónimo CHE-CNRCB correspondiente.

El ADN genómico de la mayoría de los aislados fue extraído a partir de 10 mg de micelio liofilizado y macerado con un balín de tungsteno (3 mm) o acero (5 mm) en el equipo TissueLyser II™ (Qiagen Valencia, USA) y usando el kit DNeasy Plant Mini (Qiagen, Valencia, USA) de acuerdo a las instrucciones del fabricante. La extracción de ADN de los aislados CHE-CNRCB 373 y 375 se llevó a cabo con 10 mg de micelio liofilizado macerado con nitrógeno líquido en un mortero y el procedimiento convencional descrito por ^{Amer et al. (2011)}. Los parámetros de calidad del ADN genómico extraído (i.e., pureza, concentración e integridad) se determinaron según ^{Gallou et al. (2016)}. Acorde con el género, las regiones utilizadas para la identificación fueron las siguientes: la región intergénica Bloc (5.3F/B5.1R; ^{Rehner et al., 2006}), el factor de elongación de traducción 1α (TEF) (983F/2218R; ^{Carbone y Kohn, 1999}; ^{Rehner y Buckley, 2005}) y los espaciadores internos transcritos (ITS) (ITS5/ITS4; ^{White et al., 1990}) para los aislados del género Beauveria. Las regiones TEF e ITS para los aislados del género Isaria. La subunidad pequeña ribosomal (SSU) (NS1/NS4; ^{White et al., 1990}), la subunidad 1 de la ribopolimerasa II (RPB1) (CRPB1A/RPB1Cr; ^{Castlebury et al., 2004}), el TEF y la región ITS para los aislados del género Hirsutella, y la región ITS para el género Simplicillium. La mezcla de reacción de PCR para cada región a amplificar consistió de Polimerasa Phire Hot Start II (Thermo Fisher Scientific, Waltham, USA), 1X del buffer Phire Reaction (ThermoFisher Scientific), 0.2 mM del mix de dNTPs, 0.4 µM de cada primer y 1-20 ng de ADN genómico, ajustados a un volumen final de 25 µL. Los programas de amplificación por PCR se efectuaron como sigue: para Bloc, 1 ciclo de desnaturalización inicial a 98 °C por 30 s, seguido de 35 ciclos de desnaturalización a 98 °C por 5 s, alineación a 56 °C por 5 s y extensión a 72 °C por 22 s, con 1 ciclo de extensión final a 72 °C por 1 min. TEF, 1 ciclo de desnaturalización inicial a 98 °C por 30 s, seguido de 10 ciclos de: 98 °C por 5 s, 66 °C (disminuyendo 1 °C por ciclo) por 5 s y 72 °C por 15 s, en seguida 30 ciclos de: 98 °C por 5 s, 56 °C por 5 s y 72 °C por 15 s y por último, un ciclo de extensión final a 72 °C por 1 min. ITS, 1 ciclo de desnaturalización inicial a 98 °C por 30 s, seguido de 36 ciclos a 98 °C por 5 s, 55 °C por 5 s, 72 °C por 9 s y 1 ciclo de extensión final a 72 °C por 1 min. SSU, 1 ciclo de desnaturalización inicial a 98 °C por 30 s, seguido de 4 ciclos de: 98 °C por 5 s, 55 °C (disminuyendo 1 °C por ciclo) por 5 s y 72 °C por 15 s, continuando con 42 ciclos de: 98 °C por 5 s, 50.5 °C por 5 s y 72 °C por 15 s, con 1 ciclo de extensión final a 72 °C por 1 min. RPB1, 1 ciclo de desnaturalización inicial a 98 °C por 30 s, seguido de 4 ciclos de: 98 °C por 5 s, 47 °C por 5 s y 72 °C por 10 s, continuando con 40 ciclos de: 98 °C por 5 s, 47 °C por 5 s y 72 °C por 10 s, con 1 ciclo de extensión final a 72 °C por 1 min. La purificación de los productos de PCR se realizó con el kit Wizard^® SV Gel y PCR Clean-Up System Kit (Promega, Madison, EUA). Para la secuenciación de los amplicones se utilizaron los primers de PCR y en algunos casos, primers adicionales específicos a las regiones Bloc Brint, Bfn6, B822Ldg, B22Udg, Bfn2 y BRn2 (^{Rehner et al., 2011}) y TEF 1577F y 1567R (^{Rehner y Buckley, 2005}) para los géneros Beauveria e Isaria (Macrogen, Seúl, Corea del Sur).

La edición y el ensamble de las secuencias se realizó en el programa BioEdit 7.2.0 (^{Hall, 1999}). Las secuencias correspondientes a cada región fueron depositadas en el GenBank (Tabla 1). Los alineamientos múltiples se efectuaron con MUSCLE (^{Edgar, 2004}) para el género Simplicillium y MAFFT (^{Katoh et al., 2002}) para los géneros Beauveria, Isaria e Hirsutella. El análisis filogenético para el género Beauveria se realizó por inferencia bayesiana (IB) derivada de la concatenación de las regiones Bloc, TEF e ITS. Se utilizaron 2 modelos de evolución de acuerdo con el programa jModelTest 2.1.10 (^{Darriba et al., 2012}): Hasegawa-Kishino-Yano (HKY) con una distribución gama para las regiones Bloc e ITS, mientras que para el TEF, se usó el modelo “general time reversible” (GTR) con sitios invariantes y una distribución gama. El análisis fue llevado a cabo con 1,000,000 de generaciones y un muestreo cada 100 árboles en el programa MrBayes 3.2.5 (^{Ronquist y Huelsenbeck, 2003}). Las secuencias de las cepas de referencia tipo se obtuvieron de ^{Rehner et al. (2011)}, mientras que las secuencias de Isaria cicadae ARSEF 7260, I. tenuipes ARSEF 4096, I. farinosa ARSEF 4029 y Cordyceps militaris ARSEF 5050 se utilizaron como grupo externo. El análisis filogenético para el género Isaria se realizó por IB con las regiones TEF e ITS de forma individual. Se utilizaron 2 modelos de evolución acorde al programa jModelTest 2.1.10 (^{Darriba et al., 2012}): Kimura de 2 parámetros (K80) con una distribución gama para el TEF, mientras que para el ITS, se usó el modelo HKY con una distribución gama. El análisis fue llevado a cabo con 1,400,000 generaciones y un muestreo cada 100 árboles para el TEF y con 1,000,000 de generaciones y un muestreo cada 100 árboles para el ITS en el programa MrBayes 3.2.5 (^{Ronquist y Huelsenbeck, 2003}). Las cepas de referencia tipo se obtuvieron de ^{Cabanillas et al. (2013)}, mientras que las cepas Metarhizium marquandii CBS 182.27, M. robertsii ARSEF 727 e Hirsutella sp. OSC 128575 se utilizaron como grupo externo. El análisis filogenético para el género Hirsutella se realizó por IB derivada de la concatenación de las regiones SSU, TEF y RPB1. Se utilizó el modelo GTR con una distribución gama y sitios invariantes para las 3 regiones con base al jModelTest 2.1.10 (^{Darriba et al., 2012}). El análisis bayesiano se realizó con 1,000,000 de generaciones y un muestreo cada 100 árboles en el programa MrBayes 3.2.5 (^{Ronquist y Huelsenbeck, 2003}). Las cepas de referencia tipo se obtuvieron de ^{Quandt et al. (2014)}, mientras que las cepas Metarhizium guizhouense ARSEF 5714, M. majus ARSEF 3145, Pochonia chlamydosporia CBS 504.66 y P. chlamydosporia CBS 101244 se utilizaron como grupo externo. Para el IB de los géneros Beauveria, Isaria e Hirsutella, los árboles resultantes se visualizaron en FigTree 1.4.2 (^{Rambaut, 2014}) y se calculó un árbol de consenso del 50%. El análisis filogenético para el género Simplicillium se realizó con el método de máxima verosimilitud (MV) de la región ITS. Se utilizó el modelo de evolución de Tamura de 3 parámetros con base al programa MEGA 6 (^{Tamura et al., 2013}). Se aplicó una distribución gamma para el análisis con 1,000 réplicas de bootstrap. Las cepas de referencia tipo se obtuvieron de ^{Nonaka et al. (2013}), mientras que las cepas Torrubiella luteorostrata CBS 398.86, Verticillium chlamydosporium CBS 103.65 y Pochonia chlamydosporia ARSEF 2218 se utilizaron como grupo externo.

Resultados

Las exploraciones realizadas en huertos comerciales de limón y plantas de M. paniculata permitieron recolectar insectos infectados de diferentes estadios de desarrollo con signos de micosis. Los HE asociados a ninfas y adultos de D. citri corresponden a 6 géneros de HE, de los cuales 4 pertenecen al orden Hypocreales: Hirsutella (19), Isaria (5), Beauveria (5) y Simplicillium (1), obteniendo un total de 30 aislados (Tabla 2). Además se identificaron 2 géneros del orden Entomophthorales: Batkoa y Entomophthora, de los cuales no fue posible obtener los aislados.

Hirsutella citriformis fue el hongo con mayor incidencia y distribución, se encontró en todos los huertos muestreados infectando adultos, aunque en el verano de 2014 se encontraron 3 ninfas infectadas por H. citriformis en la comunidad de Coalatilla en Armería, siendo el primer reporte de infección de esta especie sobre ninfas en México. Los PAC infectados con este hongo se encontraban adheridos a hojas y ramas en limón y limonaria, cubiertos por micelio blanco a grisáceo, su característica particular fue la presencia de sinemas largos y cilíndricos de tamaño variable que emergían del tórax del insecto, en algunos casos ramificados, regularmente con protuberancias cónicas en la base del sinema (Fig. 1a, b). La estructura conidiógena mostró forma de botella con la base ensanchada y la parte distal con un cuello delgado y alargado, distribuidas a lo largo del sinema (Fig. 1c, d) y los conidios mostraron una forma fusiforme (Fig. 1e).

Figura 1 Diaphorina citri infectada con hongos entomopatógenos. Hirsutella citriformis infección sobre adulto (a) y ninfa (b). Estructuras morfológicas de H. citriformis: sinema con fiálides (c), fialide esporulando (d), conidios sueltos (e). Isaria javanica en adulto (f) y ninfa (g), estructuras morfológicas: fiálides (h) y conidios (i). Beauveria bassiana sobre adulto (j), morfología de las fiálides y conidios (k), conidios en raquiz (l). Simplicillium lanosoniveum en ninfa (m), fiálide con esporas (n), conidios en cabezas (o). Barras de figuras (c) = 20 μm, (d) = 10 μm, (e, h, i, k, l, n y o) = 5 μm. 

Del material recolectado se obtuvieron 19 aislados que se registraron con los acrónimos CHE-CNRCB 328, 329, 330, 331,332, 333, 335, 337, 338, 339, 340, 341, 342, 345, 346, 347, 350 y 360, todos ellos de adultos y el aislado 375 obtenido a partir de ninfas del psílido. Considerando la morfometría de 5 aislados, presentaron variabilidad en sus dimensiones (Tabla 3). La longitud de las fiálides mostraron un intervalo de 24.61 - 29.32 µm de cuello y la base de las fiálides mostró valores entre 5.82 - 6.73 µm. Las dimensiones de los conidios fueron de 6.95 - 7.68 × 2.09 - 2.54 µm (largo × ancho). El análisis filogenético por IB derivada de la concatenación de las regiones SSU, TEF y RPB1 de las muestras de Hirsutella spp. mostró que los aislados CHE-CNRCB 335, 338, 339 y 375 forman parte del clado de H. citriformis y en el que se incluyen los 2 aislados H. citriformis ARSEF 1035 y ARSEF 1446 (Fig. 2). La rama tiene un apoyo fuerte y significativo de 100% de probabilidades posteriores por IB (PP IB).

Figura 2 Inferencia bayesiana derivada de la concatenación de las regiones SSU (subunidad pequeña ribosomal), TEF (factor de elongación de traducción 1α) y RPB1 (subunidad 1 de la RNA polimerasa II) (~ 2652 pb) de las cepas CHE-CNRCB 335, 338, 339 y 375. Las probabilidades posteriores (PP) > 80 se señalan en los nodos de cada rama del árbol. Las cepas de referencia proceden de ^{Quandt et al. (2014)}. 

Isaria javanica. Durante los muestreos de campo se encontraron ninfas y adultos del PAC infectadas con I. javanica en 2 huertos de limón del municipio de Tecomán, los insectos momificados tenían apariencia polvosa grisácea. Los adultos momificados se observaron adheridos al follaje en posición de alimentación con el abdomen cubierto por el hongo (Fig. 1f). En ninfas se observó micosis en la parte ventral sobresaliendo por sus costados (Fig. 1g). Se obtuvieron 5 aislados CHE-CNRCB 356, 357, 358, 359 y 363. La morfología de las células conidiógenas fue en forma de botella con la porción basal ensanchada con un ahusamiento que da origen a un cuello característico (Fig. 1h), las dimensiones de las fiálides fueron de 5.5 × 2.3 µm (4.6 - 6.3 ± 0.45 × 1.7 - 2.9 ± 0.29 µm, n = 90) y los conidios en cadenas de forma cilíndrica a fusiforme (Fig. 1i) de 4.5 ×1.6 µm (4.0 - 5.4 ± 0.41 × 1.3 - 2.0 ± 0.15, n = 90), características acordes con la especie I. javanica (^{Gallou et al., 2016}). Los análisis filogenéticos individuales por IB de las regiones ITS y TEF de las cepas de Isaria spp. mostraron que los aislados CHE-CNRCB 357, 358 y 363 forman parte del clado de I. javanica y en el que se incluye el aislado de referencia tipo I. javanica CBS 134.22 (Fig. 3). Para las 2 regiones, las ramas tienen un apoyo significativo de 100% de PP IB.

Figura 3 Inferencia bayesiana de la región ITS (espaciadores internos transcritos) (~ 600 pb) de las cepas CHE-CNRCB 357, 358 y 363. Las probabilidades posteriores (PP) > 80 se señalan en los nodos de cada rama del árbol. La cepa de referencia tipo de Isaria javanica (^{Cabanillas et al., 2013}) se señala con un (•). 

Beauveria bassiana. Se recolectaron adultos de D. citri infectados con B. bassiana en plantas de M. paniculata y en huertos de limón en los municipios de Armería y Tecomán, los insectos mostraron crecimiento micelial en el abdomen que se fue tornando polvoso, de coloración blanca a amarillo cremoso (Fig. 1j). Se observaron estructuras de reproducción asexual típicas del género, micelio septado con un diámetro de 1.8 µm (1.3 - 2.3 ± 0.29, n = 40), células conidiógenas de 2.42 µm (1.8 - 3.3 ± 0.33, n = 90) con una base globosa en forma de botella adelgazadas en la parte superior, formando un raquis curveado en zig-zag que se extiende apicalmente con un conidio por dentículo (Fig. 1k); conidios globosos a sub-globosos (Fig. 1l) de 2.5 × 2.2 µm (2.1 - 3 ± 0.23 × 1.9 - 2.6 ± 0.15, n = 90). Cinco aislados fueron obtenidos (i.e., CHE-CNRCB 167, 168, 169, 170 y 365). El análisis filogenético por IB derivada de la concatenación de las regiones Bloc, TEF e ITS de las muestras de Beauveria spp. mostró que los aislados CHE-CNRCB 167, 168, 169, 170 y 365 forman parte del clado de B. bassiana y en el que se incluye el aislado de referencia tipo B. bassiana ARSEF 1564. La rama tiene un apoyo fuerte y significativo de 100% de PP IB. Sin embargo, dentro de la especie B. bassiana se observaron 2 clados distintos con un apoyo del 100% y del 99% de PP IB, separando los aislados CHE-CNRCB 167, 168 y 365 de los aislados CHE-CNRCB 169, 170 y ARSEF 1564. Esta diferencia genética fue soportada únicamente por la región intergénica nuclear Bloc, la cual fue desarrollada inicialmente para realizar investigaciones de especiación críptica de B. bassiana (^{Rehner et al., 2006}) y posteriormente utilizada en un análisis filogenético multilocus para evaluar la diversidad del género Beauveria y el estatus taxonómico de sus especies (Rehner et al., 2011). Por ejemplo, la secuencia Bloc del aislado CHE-CNRCB 365 presentó 96% de similitud (1,452/1,510 pb) respecto a la secuencia Bloc del aislado tipo ARSEF 1564 (EF222317), mientras que las secuencias ITS y TEF del mismo aislado presentaron más de 99% de similitud (515/519 pb y 984/986 pb, respectivamente) con las secuencias ITS (NR_111594) y TEF (HQ880974) del aislado tipo ARSEF 1564. Adicionalmente, el análisis realizado con la herramienta “basic logic alignment tool” (BLAST) en la base de datos del National Center for Biotechnology Information (NCBI) mostró que la secuencia Bloc del aislado CHE-CNRCB 365 presentó 100% de similitud con la secuencia Bloc del aislado B. bassiana ARSEF 1966 (DQ384384), confirmando su identificación como B. bassiana. Debido a lo anterior, el contenido de linajes hermanos o como se observa en la Figura 4, la presencia de clados dentro de la especie B. bassiana sugiere una gran variación intraespecífica que es determinada principalmente por la región intergénica Bloc.

Figura 4 Inferencia bayesiana derivada de la concatenación de las regiones Bloc, TEF (exón del factor de elongación 1-α) e ITS (espaciadores internos transcritos) (~ 3024 pb) de las cepas CHE-CNRCB 167, 168, 169, 170 y 365. Las probabilidades posteriores (PP) > 80 se señalan en los nodos de cada rama del árbol. Las cepas de referencia tipo (^{Rehner et al., 2011}) se señalan con un (•). 

Simplicillium lanosoniveum. La especie se localizó en un huerto de Armería sobre ninfas cubiertas por micelio blanco y de aspecto granuloso (Fig. 1m), sin sinemas; al aislado se le asignó el acrónimo CHE-CNRCB 373. Esta especie se caracteriza por presentar fiálides largas y delgadas generalmente solitarias y raramente verticiladas (Fig. 1n); las dimensiones de sus fiálides fueron 29.3 × 1.2 µm (20 - 44 ± 0.57 × 0.8 - 1.7 ± 0.21, n = 30), los conidios cilíndricos se mostraron agregados en cabezas en los vértices de cada fiálide (Fig. 1n, o), constituidos por una sola célula de 3.1 × 1.3 µm (2.4 - 4.3 ± 0.43 × 0.9 - 1.6 ± 0.17, n = 30). El análisis filogenético por IB de la región ITS de la muestra de Simplicillium spp. mostró que el aislado CHE-CNRCB 373 forma parte del clado de S. lanosoniveum y en el que se incluyen los 2 aislados S. lanosoniveum CBS 704.86 y CBS 962.72 (Fig. 5). La rama tiene un apoyo significativo de 100% de PP IB.

Figura 5 Filograma derivado del método de máxima verosimilitud de la región ITS (~ 600 pb) de la cepa CHE-CNRCB 373. El soporte de las ramas del árbol procedente del análisis de bootstrap con 1,000 réplicas se representa como valor numérico en cada nodo. Las cepas de referencia tipo (^{Nonaka et al., 2013}) se señalan con un (•). 

El género Batkoa se encontró sobre ninfas de D. citri, el crecimiento del hongo exhibió una textura granulosa dando la apariencia que se encontraba cubierta por azúcar (Fig. 6a), con conidios primarios subesféricos con una papila característica de la especie (Fig. 6b), con dimensiones de 29.5 ± 2.5 µm de largo × 24.14 ± 2 µm (n = 15) ancho, conteniendo un promedio de 8 núcleos por conidio con un diámetro de 4.25 ± 0.67µm (n = 60). No se logró aislar este género por lo cual no se realizó el análisis molecular correspondiente.

Figura 6 Diaphorina citri infectada con Entomophthorales. Batkoa sp. sobre ninfa infectada (a) y conidio primario (b). Entomophthora sp. en adulto (c) y conidio primario (d). Núcleos teñidos con aceto-orceina. Barra figura (b) = 10 μm y figura (d) = 5 μm. 

El género Entomophthora fue identificado sobre adultos; las características de los insectos infectados con este hongo fueron abdomen hinchado, las alas extendidas dorso-lateralmente y elevadas por encima del abdomen, con una masa de hifas y conidióforos aglutinados emergiendo a través de las membranas intersegmentales del abdomen y el tórax (Fig. 6c). Se observaron conidios primarios con una base plana, rodeada con un halo de mucílago, con punto apical prominente y papila basal amplia (Fig. 6d), típicos del género Entomophthora (^{Keller y Dhoi, 2007}), con 5 núcleos por conidio con un diámetro de 2.83 ± 0.32 µm (n = 40). El promedio de las dimensiones de los conidios primarios fue de 14.59 µm ± 0.73 × 11.9 µm ± 0.61 (n = 20) de largo por ancho respectivamente. No se obtuvieron aislados de este género.

Discusión

En México, se ha reportado a H. citriformis infectando a D. citri en la mayoría de las zonas citrícolas que incluyen Chiapas, Quintana Roo, San Luis Potosí, Tabasco, Veracruz, Yucatán (^{Pérez-González et al., 2015}), Nuevo León y Tamaulipas (^{Casique-Valdes et al., 2011}). Hirsutella citriformis se encontró en todos los huertos muestreados, por lo que fue el hongo con mayor incidencia y distribución en Colima. La presencia de insectos infectados se dio durante los 3 años de recolecta, aunque se observó una mayor abundancia de septiembre a febrero. Al respecto, ^{Hall, Hentz et al. (2012)} reportan mayor abundancia de psílidos micosados con H. citriformis en otoño e invierno y afirman que existe una correlación positiva entre la humedad relativa (HR) y la proporción de adultos infectados, cuando esta condición se mantiene estable por 7 días. La infección de H. citriformis sobre adultos es más común que en ninfas, durante la exploración se encontraron ninfas micosadas por el hongo en huertos de limón, en la localidad de Coalatilla en Armería. El desarrollo del hongo sobre las ninfas del psílido fue favorecido por un periodo de lluvias intensas que generaron condiciones de alta HR. Estas observaciones coinciden con lo reportado por ^{Étienne et al. (2001)}, quienes encontraron a H. citriformis atacando ninfas y adultos en Guadeloupe, Francia durante periodos con lluvias intensas provocando una mortalidad de más de 80% de psílidos. Bajo condiciones similares, ^{Álvarez et al. (2003)}, indican la infección del hongo sobre ninfas de PAC en Jovellanos, Matanzas, Cuba.

La morfometría de los conidios de H. citriformis es similar a la descrita por ^{Mains (1951)}, ^{Subandiyah et al. (2000}), ^{Álvarez et al. (2003}), ^{Meyer et al. (2007}) y ^{Toledo et al. (2013)} (Tabla 3); sin embargo, mostró mayores dimensiones conidiales que con otras cepas aisladas en México reportadas por ^{Pérez-González et al. (2015)}. Estos autores afirman que los aislados de H. citriformis recolectados en México presentan variabilidad en las dimensiones de sus estructuras.

La presencia de I. javanica fue menos abundante, se ubicó únicamente en 2 huertos de limón en Tecomán, en el primero se encontró infectando ninfas, en el segundo sobre adultos. Isaria fumosorosea Wize ha sido asociada a D. citri en México (Casique-Valdés et al., 2015; ^{Lezama-Gutiérrez et al., 2011}) y en otros Países (^{Subandiyah et al., 2000}). Las cepas de I. javanica obtenidas en nuestro estudio presentaron estructuras de mayor tamaño a las reportadas para la cepa DCP03 de I. fumosorosea cuyas dimensiones para las fiálides son de 4.7 × 1.9 µm y conidios de 2.9 × 1.4 µm (Subandiyah et al., 2000), pero similares a las cepas de I. javanica caracterizadas por ^{Gallou et al. (2016)}.

Aunque B. bassiana es una especie que ocasiona infección sobre muchas especies de insectos (^{Zimmermann, 2007}), su presencia fue baja. ^{Lezama-Gutiérrez et al. (2011)} encontraron ninfas y adultos de D. citri infectados con B. bassiana en Colima, mientras que en el presente estudio, solo se hallaron adultos infectados en un huerto de limón y en plantas de limonaria en los municipios de Armería y Tecomán. Igualmente se ha detectado la presencia de este hongo en el estado de Nayarit (^{Rodríguez-Palomera et al., 2012}) en limón persa con bajos porcentajes de infección.

El género Simplicillium tiene parentesco con el género Lecanicillium, ambos géneros producen fiálides delgadas en forma de espina, conidios que se adhieren en cabezas viscosas mucilaginosas o formando cadenas. Simplicillium y Lecanicillium se pueden distinguir por la disposición de fiálides, el primero generalmente posee una disposición de fiálides solitarias con poca diferencia entre las hifas, mientras que este último tiene fiálides verticiladas. Los caracteres morfológicos del aislamiento CHE-CNRCB 373 coinciden con la descripción de S. lanosoniveum (^{Zare y Gams, 2001}), en este género se incluyen aislados de diversos artrópodos, otros hongos y suelo (^{Bałazy et al., 2008}; ^{Nonaka et al., 2013}; ^{Zare y Gams, 2001}). Se reporta por primera vez la presencia de S. lanosoniveum sobre D. citri.

Los Entomophthorales Batkoa sp. y Entomophthora sp. fueron detectados sobre ninfas y adultos de D. citri respectivamente, este es el primer reporte de Batkoa sp. asociado al PAC en Colima, México. Estos géneros se presentaron en los meses de diciembre a febrero cuando las condiciones ambientales son frescas por la mañana con presencia de rocío. Los 2 géneros fueron identificados solo morfológicamente ya que no fue posible su aislamiento. Las características morfométricas de Batkoa sp. se confirmaron considerando las dimensiones reportadas por ^{Keller y Dhoi (2007}) y para Entomophthora sp. la forma de campana del conidio primario, número de núcleos y las dimensiones concuerdan con la descripción de ^{Keller (2007)}. Los conidios primarios (CP) y sus núcleos presentaron dimensiones similares a la espora de Entomophthora philippinensis Villacarlos y Wilding (CP: 11.5 - 20.5 × 9 - 18.5 μm; con 2 - 6 núcleos de 2 - 4.5 × 2 - 3.5 μm) reportado infectando al psílido de la leucaena, Heteropsylla cubana Crawford, en las Filipinas (Villacarlos y Wilding 1994), y más pequeños que los reportados sobre D. citri (CP: 16.6 × 22.7 μm) en Veracruz, México (^{Guízar-Guzmán y Sánchez-Peña, 2013}). Los insectos infectados con Batkoa sp. se encontraron en un solo huerto, mientras que Entomophthora sp. estaba presente en huertos de los 3 municipios explorados encontrando pocos insectos infectados por huerto. En algunos de los casos los insectos que inicialmente fueron infectados por Entomophthora presentaron estructuras típicas del género Cladosporium, lo anterior se evidenció al encontrar insectos con estructuras de los 2 hongos, pero con los síntomas característicos de infección por Entomophthorales, al respecto ^{Guízar-Guzmán y Sánchez-Peña (2013)} afirman que el porcentaje de infección causado por Entomophthora es posiblemente subestimado por el rápido crecimiento de hongos secundarios. Los hongos del orden Entomophthorales asociados a D. citri reportados previamente en México son Zoophthora radicans (Brefeld) Batko para Colima (^{Lezama-Gutiérrez et al., 2011}), el cual no fue encontrado en esta búsqueda y Entomophthora sp. en el estado de Veracruz (^{Guízar-Guzmán y Sánchez-Peña, 2013}).

Considerando los resultados obtenidos y otros reportes de hongos asociados al PAC en México se cuenta con una gran diversidad de HE. La obtención de aislados de hongos patógenos presentes en las poblaciones de D. citri bajo determinadas condiciones ambientales, deben ser consideradas para generar tecnología que en un futuro permita incluir acciones adicionales para el manejo de la plaga como parte del programa de control biológico del PAC implementado por la DGSV. En condiciones naturales se observó a H. citriformis con un alto potencial infectivo sobre D. citri, lo que sugiere el desarrollo de técnicas para su reproducción masiva; además se requiere desarrollar estrategias para el uso de Entomophthorales con el fin de explotar su potencial como agentes de control de la plaga.