El huanglongbing (HLB) es considerada, actualmente, la enfermedad de mayor importancia a nivel mundial en la producción citrícola, ya que causa la muerte de los árboles algunos años posterior a la infección, además de que, a la fecha, no existe ningún método de control efectivo y económicamente viable, así como tampoco variedades resistentes al alcance del productor. Se relaciona a la enfermedad con la infección de bacterias del género Candidatus Liberibacter (Ca. L.), así, se ha observado que el HLB es causado por Ca. L. africanus en Sudáfrica (Roberts y Pietersen, 2017), por Ca. L. americanus en Brasil (Wulff et al., 2014) y más extensamente por Ca. L. asiaticus, por ejemplo en Estados Unidos de América (Manjunath et al., 2008) y México (Flores-Sánchez et al., 2015). Estas bacterias son patógenos obligados en el floema de los cítricos y endosimbiontes en diferentes sistemas de su vector, Diaphorina citri (Pelz-Stelinski y Killiny, 2016; Killiny et al., 2018), por lo cual la obtención de cultivos axénicos es muy compleja y dificulta los estudios de patología (Ha et al., 2019; Merfa et al., 2019). En aras de comprender los mecanismos moleculares que permiten el desarrollo del HLB, se ha secuenciado el genoma de varias especies de Ca. Liberibacter spp. (Duan et al., 2009; Lin et al., 2013; Fagen et al., 2014; Zheng et al., 2014, 2015; Cai et al., 2018). Esto ha permitido descubrir genes con alto potencial de patogenicidad en la interacción con la planta. Por ejemplo, Vahling et al. (2010) detectaron la presencia de un gen codificante para una enzima ATP/ADP translocasa funcional y con alta afinidad al ATP y ADP, lo cual resulta importante ya que se ha documentado que una de las respuestas de las plantas a la infección por HLB es la acumulación excesiva de ATP (Pitino et al., 2017). Esta enzima ha sido ampliamente relacionada con endosimbiontes intracelulares obligados, la cual tiene la capacidad de obtener ATP desde la célula del hospedante e importarlo dentro de su propia célula. Este mecanismo se conoce como parasitismo energético (Schmitz-Esser et al., 2004).Lo anterior resulta especialmente de interés debido a que previamente no se conocía este tipo de habilidades parasitarias en fitopatógenos, sino únicamente en patógenos humanos y de animales (Trentmann et al., 2007). La historia evolutiva de esta característica patogénica sugiere que el origen de este gen proviene de un endosimbionte ancestro al actual grupo de Ricketssiales que, tentativamente, se mueve de forma horizontal entre otros grupos de patógenos obligados (Emelyanov, 2007).
El presente trabajo tiene como objetivo analizar las relaciones evolutivas entre genes ATP/ADP translocasa de diferentes patógenos y endosimbiontes humanos y animales y los encontrados en el género Ca. Liberibacter, incluyendo los patógenos causantes del HLB de los cítricos. Para esto, utilizando como base la secuencia de aminoácidos reportada por Vahling et al. (2010), se realizó un análisis BLASTp para identificar las secuencias proteicas con mayor similitud de acuerdo a un e-value <0.001 y un porcentaje de identidad >90% en la base de datos del Genbank pertenecientes a Ca. Liberibacter y a grupos como Ricketssia spp. y Chlamydia spp.
Una vez obtenidas las secuencias de aminoácidos se realizó un alineamiento múltiple utilizando el algoritmo ClustalW (gap open=15; gap extend=3) mediante el software Bioedit (Hall, 1999). Posteriormente, una búsqueda del árbol filogenético óptimo mediante el criterio de la Máxima Parsimonia, para lo cual se utilizó el software T.N.T. (Goloboff et al., 2008), para esto se utilizó una búsqueda mediante algoritmos compuestos (ratchet + sectorial search + drift) y se construyó un árbol consenso estricto. Se realizó un remuestreo mediante bootstrap (1000 iteraciones) utilizando las opciones ofrecidas por el software. En seguida, se elaboró un análisis de conservación entre las secuencias de aminoácidos de la enzima en estudio. Para esto, el alineamiento fue analizado para determinar la conservación por sitios específicos mediante el software Jalview (Waterhouse et al., 2009), adicionalmente se realizó un análisis de entropía posicional mediante la función:
Dónde f(b;l) es la frecuencia con que cada aminoácido aparece en cada columna, siendo mínima cuando la posición está totalmente conservada y máxima cuando todos los aminoácidos están igualmente representados. El análisis fue realizado en el software Bioedit.
Se utilizaron en total 39 secuencias de la enzima ATP/ADP translocasa, como grupo externo para enraizar el árbol filogenético, se incorporó la secuencia de la enzima ATP/ADP translocasa de plastidios de Arabidopsis thaliana (Genbank CAA89201). El análisis filogenético permitió recuperar ocho árboles igualmente parsimoniosos (L=2502), con los cuales se construyó un árbol de consenso estricto (Figura 1). Los resultados del análisis muestran que todas las secuencias provenientes de Ca. Liberibacter spp. se aglomeran en un solo clado (remarcado en azul), este clado es grupo hermana con la enzima proveniente del endosimbionte Candidatus Midichloria mitochondrii, asociado a las garrapatas causantes de diversas enfermedades alrededor del mundo, la cual tiene una característica interesante al ser una bacteria patógena intramitocondrial (Sassera et al., 2006).
Es destacable que, dentro del clado antes mencionado, se ubican subclados bien definidos y con soporte bootstrap alto de acuerdo a cada especie del género. El ancestro hipotético basal del clado se comparte entre Candidatus Midichloria mitochondrii y Ca. Liberibacter aeuropaeus, esto resulta de interés ya que ambos organismos son considerados endosimbiontes no patogénicos en garrapatas (Najm et al., 2012) y en perales (Raddadi et al., 2011; Camerota et al., 2012), respectivamente. Los subclados de Ca. Liberibacter americanus y Ca. Liberibacter africanus se encuentran igualmente bien definidos, los cuales son agentes patógenos causantes de HLB en Brasil (do Carmo Teixeira et al., 2005; Teixeira et al., 2008) y en África (Lin et al., 2006). Igualmente, la presencia de Ca. L. solanacearum y Ca. L. asiaticus como subclados bien definidos y soportados permite plantear dos hipótesis: a) la capacidad patogénica de algunas especies de Ca. Liberibacter pudo haber sido desarrollada en un ancestro endófito no patogénico, esto debido a que Ca. L. aeuropaeus es muy basal dentro del clado de los Ca. Liberibacter y b) los subclados bien definidos por especie sugieren una coevolución de la patogenicidad con los hospederos.
Las secuencias del clado Ca. Liberibacter muestran un alto grado de conservación en algunas regiones, y el motivo relacionado con el transporte de ATP se encuentra intacto en todas ellas (datos no mostrados). Sin embargo, algunas regiones son altamente variables, siguiendo un patrón de acuerdo a la especie a que pertenece la secuencia de la enzima, lo cual podría generar una cambio en la dinámica estructural y en la función de la proteína (Liu y Bahar, 2012). En la Figura 2A se muestran los sitios con mayor nivel de conservación y los que presentan una mayor variabilidad de acuerdo a la frecuencia en el alineamiento. El análisis de entropía (Figura 2B) demuestra que varias posiciones de la secuencia de la enzima ATP/ADP translocasa contiene un nivel entrópico que sugiere una alta variabilidad en otras posibles secuencias, por tanto, de encontrarse nuevas especies dentro de Ca. Liberibacter con hospederos no reportados, la variación en las secuencias podría ser específicas a cada hospedero. Ahondar en el estudio de la función y la evolución de este gen permitirá comprender mejor los mecanismos por los cuales se genera el daño en la planta durante el desarrollo del HLB y coadyuvar al desarrollo de posibles estrategias de control de la enfermedad. Los resultados presentados en este estudio sugieren que la capacidad de parasitar ATP que tiene Ca. Liberibacter spp. es una característica compartida con endosimbiontes y endopatógenos; sin embargo, la secuencia de la enzima responsable por este parasitismo varía de acuerdo a cada especie dentro del género lo que sugiere eventos de coevolución del parasitismo con el hospedante.