¹ Universidad del Valle, Facultad de Ciencias Exactas y Naturales, Programa de Biología. Biólogo. Cali, Colombia.

² Universidad del Cauca, Facultad de Ciencias Agropecuarias, Departamento de Ciencias Agropecuarias. Ph.D. Popayán, Colombia.

³ Universidad Nacional de Colombia, Facultad de Ciencias Agropecuarias. Ph.D. Palmira, Colombia.

⁴ Universidad del Valle, Facultad de Ciencias Exactas y Naturales, Sección Entomología. Ph.D. Cali, Colombia.

⁵ Corporación Agropecuaria de Colombia, Entomología. Ph.D. Palmira, Colombia.

Correspondencia: miguelangelmoncayodonoso@gmail.com

Recibido para evaluación: 23-02-2014. Aprobado para publicación: 07-04-2014.

RESUMEN

El psílido asiático Diaphorina citri es el principal vector de la bacteria Candidatus Liberibacter que causa la enfermedad Huanglongbing, considerada devastadora en la citricultura mundial y aún no reportada en Colombia. La variabilidad genética poblacional de D. citri fue estudiadautilizando la técnica de secuenciación del gen COI mitocondrial como marcador molecular. Se realizaron colectas de adultos en zonas productoras de cítricos del Valle del Cauca y el Quindío. La amplificación se realizó con dos parejas de cebadores específicos para hemípteros. Los productos de la PCR fueron secuenciados en Macrogen-Korea obteniéndose un total de 124 secuencias. En el análisis bioinformático se utilizaron los programas Vector NTI 11.5, Arlequín V 3.5, MEGA 5 y MAFFT 6. Los índices de diversidad molecular fueron similares entre poblaciones, mostrando un mismo origen, separación reciente de poblaciones y que no existe diferenciación genética significativa que pueda albergar variantes de la bacteria; sin embargo, el índice de diversidad haplotípica fue más alto que el de diversidad nucleotídica y este último índice presentó bajo número de sitios polimórficos, indicando que las poblaciones de D. citri están en expansión. El estudio de la variabilidad genética del vector representa una herramienta para predecir escenarios probables de propagación de la enfermedad.

PALABRAS CLAVE: ADN Mitocondrial, Reverdecimiento de los Cítricos, Insectos Vectores, Huanglongbing, Psylido Asiático de los Cítricos.

ABSTRACT

The Asiatic psyllid Diaphorina citri (Hemiptera: Psyllidae) is the main vector of Candidatus liberibacter, which causes the Huanglongbing HLB disease, known for devastating citrus in the world but not yet reported in Colombia. The genetic variability of the D. citri population was studied through sequencing the COI mitochondrial gene as molecular marker. Adults were collected in citrus producing zones of the Colombian Valle del Cauca and Quindío. Amplification was performed with two pairs of specific primers for Hemiptera. The PCR products were sequenced at Macrogen-Korea, obtaining a total of 124 sequences. For the bioinformatic analysis, the Vector NTI 11.5, Harlequin V 3.5, MEGA 5 and MAFFT 6 programs were used. The molecular diversity indices between populations were similar, revealing a common origin and a recent split of the populations excluding a significant genetic differentiation associated to variations of the bacterium, however the haplotype diversity index was higher than the nucleotide diversity index. The latter one showed a low number of polymorphic sites, indicating that the D. citri populations are expanding. The study of the vector’s genetic variability is a tool for the prediction of likely scenarios for the spread of diseases.

KEYWORDS: Mitochondrial DNA, Citrus Regreening, Vector Insects, Huanglongbing, Asiatic Citrus Psyllid.

RESUMO

O psílido asiático Diaphorina citri, é o principal vetor da bactéria Candidatus Liberibacter que causa a doença Huanglongbing HLB, considerada devastadora na citricultura mundial e que ainda não foi registrada na Colômbia. A variabilidade genética populacional de D. citri foi estudadautilizando-se a técnica de seqüenciamento do gene COI mitocondrial como marcador molecular. Foram feitas coletas de adultos nas áreas produtoras de citros do Valle del Cauca e Quindío. A amplificação foi realizada com dois pares de iniciadores específicos para hemípteros. Os produtos de PCR foram seqüenciados em Macrogen-Korea obtendo-se um total de 124 seqüencias. Para as análises de bioinformática foram utilizads os programas Vector NTI 11.5, Arlequín V 3.5, MEGA 5 e MAFFT 6. Os índices de diversidade molecular foram similares entre populações, apresentando uma mesma origem, separação recente de populações, e que não existe uma diferenciação genética importante que possa abrigar variantes da bactéria; no entanto, o índice da diversidade haplotípica foi mais alto do que o da diversidade nucleotídica e este último índice apresentou um baixo número de sítios polimórficos, indicando que as populações de D. citri estão em expansão. O estudo da variabilidade genética do vetor representa uma ferramenta para predizer cenários prováveis de propagação da doença.

PALAVRAS CHAVE: DNA Mitocondrial, Citros Greening, Insetos Vetores, Huanglongbing, Psilídio Asiático dos Citros.

INTRODUCCIÓN

El psílido asiático de los cítricos Diaphorina citri Kuwayama es importante por ser el vector de la bacteria Candidatus Liberibacter, agente causal de la Huanglongbing HLB o enverdecimiento de los cítricos “greening”, enfermedad devastadora de la citricultura mundial [1]. El insecto causa daños directos, al alimentarse de la planta e indirectos cuando transmite la bacteria causante de la enfermedad; afecta zonas de crecimiento de la planta ocasionando la caída y deformación de hojas, clorosis, necrosis, enanismo y disminución de la calidad de la fruta; las ninfas y los adultos ocasionan daños al extraer gran cantidad de floema de las hojas y los pecíolos; además, poblaciones altas del insecto causan la caída de flores y frutos pequeños [2]. La bacteria, se caracteriza como proteobacteria del genero Candidatus Liberibacter para la enfermedad Huanglongbing, enfermedad generada por Candidatus Liberibacter asiaticus, Candidatus Liberibacter africanus y Candidatus Liberibacter americus. La enfermedad tiene alta incidencia en Asia y África donde se señala como factor limitante para la producción de cítricos [3].

Existen diversos métodos y estrategias de control para enfrentar, tanto al vector como a la enfermedad, que incluyen la destrucción y eliminación de las fuentes de inóculo, la renovación de las plantaciones utilizando posturas sanas, el control o la erradicación del insecto vector mediante químicos y/o un programa de control integrado. En el caso de establecerse el vector y no detectarse la enfermedad, entonces se plantea la implementación de programas de control biológico. Sin embargo, hasta el momento no se conoce ningún método efectivo, aunque el control biológico clásico puede contribuir a la supresión de las poblaciones del psílido [6].

El insecto vector D. citri está ampliamente distribuido en las regiones tropicales y subtropicales de Asia [4]. A principios de 1998, fue descubierto en la isla de Guadalupe en el Caribe y luego en Argentina, Paraguay, Uruguay, Brasil, Venezuela, Colombia y México [5]. También, se encontró en la Florida y con alta probabilidad de dispersarse hacia los cítricos de Louisiana, Texas, Arizona y California [6]. El desarrollo del insecto es exclusivamente en plantas de la familia Rutaceae, particularmente en los géneros Citrus spp. y Murraya y en las especies, Swinglea glutinosa M., Limonia acidissima L. y Citropsis schweinfurthii S. y K. [1].

Colombia presenta una topografía y biomodalidad climática que favorece la diversidad de ambientes, además, las condiciones diferentes en que se cultivan los cítricos en el Valle del Cauca y en el Quindío, pueden generar diferentes ambientes que favorezcan la diferenciación génica en D. citri. Según [7], un análisis comparativo de las similitudes genotípicas de los haplotipos de poblaciones D. citri de Colombia con las del resto del mundo sugiere que las colombianas provienen del sudeste asiático, considerado el centro con mayor probabilidad de haber originado este psílido.

Es necesario conocer la diversidad genética de D. citri a nivel molecular para establecer estrategias de manejo. El ADN mitocondrial es una herramienta útil para estudios de filogenia, diversidad genética y estructura poblacional, dada su transmisión exclusiva por vía materna, sin recombinación y tasa de mutación relativamente rápida [9]. Esta última característica es conveniente para análisis genéticos intra e interespecíficos. Los genes de la citocromo oxidasa, subunidades I y II (COI y COII) son los más estudiados del ADN mitocondrial y las unidades terminales del análisis filogenético son los haplotipos mitocondriales o variantes del ADN mitocondrial. El número y estructura de los haplotipos presentes en las poblaciones de una especie aportan medidas de diversidad y estructura genética poblacional y permiten realizar inferencias sobre su historia evolutiva [3].

El presente estudio analizó la diversidad genética poblacional de D. citri mediante el análisis de las secuencias de los genes mitocondriales COI, en dos regiones productoras de cítricos en Colombia. Los resultados de estos análisis permitirán diseñar medidas de control de D. citri, identificar poblaciones de insectos con riesgo futuro de transmisión de la enfermedad HLB y establecer el potencial de dispersión y flujo genético de D. citri. Cabe destacar que este estudio es el primer trabajo de caracterización molecular de poblaciones de D. citri en Colombia.

MÉTODO

Material biológico

Los adultos de D. citri (48 ejemplares) fueron colectados manualmente en árboles de vivero y en producción, con pincel o con un aspirador, almacenados en viales con alcohol (96%) y rotulados con la información de localización, hospedero y fecha de colecta. En total ocho sitios de muestreo (cuadro 1).

Cuadro 1

Extracción de ADN

El protocolo de extracción de ADN efectivo fue el basado en sales modificado del descrito por [10], se basa en un tampón de homogenización estéril (0,4 M de NaCl, 10 mM de Tris pH 8,0 y 2 mM EDTA pH 8,0), SDS al 20%, proteinasa K 20 mg/mL y NaCl 6M. La cuantificación se realizó visualizando el ADN en geles de agarosa al 0,8% en tampón TBE 0,5X (tris-borato 0,045M; EDTA 0,001M) a 100V durante una hora. Los geles se tiñeron con bromuro de etidio a una concentración final de 0,5 ug/mL. La cantidad de ADN extraído de buena calidad, se cuantificó comparándolo con ADN del bacteriófago Lambda a una concentración de 30 ng/μL La comparación se realizó usando 1 μL, 2 μL y 3 μL de Lambda para obtener una concentración de 30, 60 y 90 ng/μL, respectivamente.

Amplificación del ADN por PCR

Los amplificados del gen COI se obtuvieron por PCR a partir del ADN genómico, se obtuvo un fragmento de 375 pb aproximadamente, determinado por la migración en el gel de agarosa con respecto a un estándar de tamaños moleculares. La concentración de ADN varió entre 30 y 40 ng/mL. El coctel de reacción se llevó a un volumen final de 50 μL, 1,5 mM MgCl2, 50-200 μM dNTP, 0,06 μM de cebadores C1-J-1763/C1-N-2096d, 1 ng de ADN, y 2 unidades de Taq polimerasa; en tampón de reacción 1X (50 mM KC1, 10 mM Tris, 0,1% Tritón X-100, pH 9,0 a 25°C). Se siguieron las condiciones de amplificación descritas por [11]: denaturación inicial a 94ºC por cuatro min, 25 ciclos a 95ºC por 40 seg; un minuto de alineamiento a 40-50ºC para el par de cebadores para la amplificación de ADN mitocondrial de D. citri y por último, un minuto de extensión a 72ºC. Los cebadores utilizados para la caracterización molecular (cuadro 2), amplifican para hemípteros, cubren las partes más variables de la COI y pueden emplearse para estudios genéticos a nivel intraespecífico [11]. La visualización de los productos de la PCR se realizó en un gel de agarosa al 1,2% y Teñido con bromuro de etidio en un transiluminador ultravioleta.

Cuadro 2

Secuenciación, edición y alineamiento de secuencias

Los productos de la PCR fueron secuenciados por la casa comercial Macrogen® (Seul, Korea) en un secuenciador automático 3730xl DNA analyzer™, en las dos direcciones por cada ejemplar, obteniéndose un total de 124 secuencias.

Las secuencias fueron editadas con el programa Vector 11.5, realizando una secuencia consenso para cada muestra a partir de las secuencias “reverse” y “forward”; y se corroboró su identidad biológica mediante un BLAST y en relación con las secuencias reportadas en las bases de datos para variantes génicas de la misma especie. Lo anterior se realizó con MAFFT vers. 6.

Análisis genético poblacional

El análisis estadístico de la información genética del locus COI, entre poblaciones de D. citri, fue inferido mediante la comparación de las secuencias de haplotipos, aplicando el programa Arlequín V. 3.5. Los índices de diversidad molecular aplicados fueron:

Índice de diversidad haplotípica (H). Estima la variabilidad genética en términos del número de haplotipos presentes en la muestra y sus respectivas frecuencias [12]. Es definido como la probabilidad de que dos haplotipos de la población, elegidos al azar, sean diferentes.

Índice de diversidad nucleotídica (pn). Estima la variabilidad genética en términos del número promedio de nucleótidos diferentes por sitio entre pares de secuencias de haplotipos. Es definido como la probabilidad de que dos nucleótidos con homología de posición escogidos al azar sean diferentes [13, 14].

Número promedio de diferencias de bases entre las secuencias (p). Estima la variabilidad genética de nucleótidos por sitio entre dos haplotipos, considera que no todas las sustituciones tienen la misma probabilidad de ocurrir y corrige por las diferencias entre las tasas de sustitución de transición y las de transversión.

RESULTADOS

Distribución de Haplotipos

Los 24 individuos de D. citri que constituyeron la muestra del Quindío presentaron mayor cantidad de haplotipos diferentes en comparación con la muestra de igual número de individuos del Valle del Cauca (17 y 14 haplotipos, respectivamente) (cuadro 3). Estos primeros resultados indican que cada región geográfica registra una historia de 17 y 14 linajes matrilineales, respectivamente, sugiriendo que son muestras representativas de puestas o eventos mutacionales diferentes, por lo tanto, los resultados son válidos para realizar análisis de variabilidad genética intraespecífica en las dos regiones.

El Valle del Cauca presenta mayor cantidad de individuos que comparten haplotipos en comparación con el Quindío: los 15 y 11 individuos, respectivamente, comparten 5 haplotipos en 4 poblaciones (cuadro 3); las poblaciones de Jamaica y La Rivera y las de Brasil y Corpoica comparten 4 haplotipos y 1 haplotipo, respectivamente, mientras que en el Quindío, las poblaciones de San Juan y Calarcá comparten 3 haplotipos y las de Caracolines y la vía a Calarcá comparten 2 haplotipos.

Los haplotipos compartidos con frecuencia alta entre individuos de diferentes regiones podrían indicar: primero, que son haplotipos característicos de la especie D. citri para estas regiones y mientras más antiguos sean estos haplotipos, mayor será su distribución geográfica [15]); segundo, que existe un mismo origen compartido; tercero, según [16] estos haplotipos han sido conservados por selección natural debido a un mayor éxito reproductivo. Finalmente, el hecho de que se encuentren varios individuos asociados a un mismo haplotipo (en general existe un haplotipo mitocondrial preponderante por cada individuo) es un indicio de homoplasmia mitocondrial.

El Valle del Cauca y Quindío presentan gran cantidad de haplotipos únicos (específicos o exclusivos), 9 y 12, respectivamente, y ningún haplotipo predominante. Según [17] los haplotipos únicos son resultado de mutaciones puntuales, y tienden a desaparecer de las poblaciones. [18] señala que la aparición de un número alto de haplotipos únicos es resultado de la ocurrencia de más de un caso de introducción dentro de las poblaciones en las cuales los grupos introducidos ya presentaban una variabilidad alta. [19] señalan que los plaguicidas ejercen una fuerte presión de selección y su aplicación periódica puede estar coincidiendo con la eliminación selectiva de haplotipos, lo que puede favorecer la expansión de los haplotipos únicos producidos por mutaciones puntuales sobre todo sí estos se asocian con los genes implicados en la resistencia a insecticidas.

Variabilidad Poblacional

Las dos regiones geográficas muestran alta variabilidad genética intrapoblacional que se refleja en sus altos índices de diversidad haplotica (H), siendo los índices H del Quindío más altos que los del Valle del Cauca (cuadro 4), resultado que concuerda con la mayor cantidad total de haplotipos presentes en la región del Quindío, y mayor cantidad de haplotipos únicos por región y por localidad. Es importante destacar que la similitud de los altos índices de diversidad haplotica de La Rivera en el Valle del Cauca con los índices de las localidades San Juan, Caracolines y Calarcá en el Quindío es debida probablemente a su relativa cercanía geográfica.

En la cuadro 4 se observa que a pesar de la alta diversidad haplotípica (H) de las dos regiones, la probabilidad de hallar un cambio nucleotídico en una posición concreta de la secuencia es reducida (baja tasa de sustitución nucleotídica); la región de Quindío presentó el mayor índice de diversidad nucleotídica en comparación con el Valle del Cauca, oscilando entre 0,0045±0,0039 y 0,0603±0,0346 y entre 0,0049±0,0036 y 0,0107±0,0081, respectivamente. Las poblaciones de San Juan, Calarcá y Caracolines (Quindío) presentaron relativamente los más altos índices de diversidad nucleotídica (0,0603±0,035, 0,024±0,014 y 0,0225±0,014, respectivamente); y los mayores promedios de diferencias de bases por sitio entre pares de secuencias (26,38±13,21, 8,86±4,65 y 8,33±4,50, respectivamente).

Los índices de diversidad nucleotídica obtenidos para D. citri son normales, según [20] las diferencias nucleotídicas dentro de las especies varía de 0 a 1,8%. Sin embargo, [18] indica que alta diversidad haplotípica acompañada de una baja diversidad nucleotídica significa que existen muchos haplotipos genéticamente similares, es decir, hay pocos cambios para proporcionar haplotipos diferentes, lo cual es característico de efectos de expansión de las poblaciones en un tiempo reciente. La expansión reciente de las poblaciones de D. citri también es apoyada por valores bajos de diversidad nucleotídica en relación al bajo número de sitios polimórficos, indicativo del crecimiento de las poblaciones o la presencia de un equilibrio en la selección [21].

Se observa una tendencia ascendente en relación con el número de sitios polimórficos en el corredor geográfico del Valle del Cauca al Eje Cafetero alcanzando 82 sitios en la población de San Juan (Quindío) (cuadro 4), debida probablemente a las diferentes condiciones climáticas/ ecológicas/geográficas y antrópicas de las dos regiones geográficas. El corredor se encuentra en gran parte de los Andes occidentales, al centro-occidente de Colombia entre las cordilleras Occidental y Central; presenta características de relieve andino donde predominan numerosas ciénagas y embalses y ríos de primer orden en el País como el Cauca y el Magdalena [22]. Además, el corredor tiene una fuerte influencia humana, el 24% de los cultivos de cítricos del País se encuentran entre el Valle del Cauca, Quindío y Antioquia, y existe una conexión en el transporte de cítricos, tanto de plántulas como de frutos [23]. Las condiciones mencionadas juegan un papel fundamental en la distribución y flujo génico de las poblaciones de una especie y favorecen la diversidad genética, aumentando el número de sitios polimórficos para las poblaciones de insectos ubicadas en éste corredor geográfico.

En las regiones estudiadas solo se observaron mutaciones puntuales siendo las más frecuentes, las transversiones que variaron entre 1 (Vía Calarcá) y 7 (Caracolines) en el Quindío mientras que las localidades de Jamaica y La Rivera en el Valle del Cauca no presentaron ninguna. Las transiciones fueron poco frecuentes: de 1 a 2 en el Valle del Cauca y de 2 a 3 en el Quindío (cuadro 4).

Los cambios en los niveles de diversidad genética de las especies o poblaciones son consecuencia de un balance entre procesos opuestos de ganancia (mutaciones y flujo génico) y perdida (deriva génica y selección natural) [3]. En insectos los cambios en la variabilidad genética, probablemente ocurren más por perdida que por ganancia, debido a diferentes factores de disturbio en el ambiente, como las presiones de selección que ejercen los insecticidas, la disponibilidad de alimento, la resistencia que pueden generar las propias plantas hospederas y factores biológicos que intervienen en la dinámica poblacional: patrones de apareamiento, tipo de reproducción, tiempo de duración de las generaciones, tasas de mortalidad, migraciones, etc. [24]. En el caso de D. citri el corto tiempo de las generaciones y el número alto de ovoposiciones siempre aseguran el tamaño de la población. Estas características entre otras le permiten a este insecto adaptarse a las condiciones ambientales variables de la región donde se ubique [4].

CONCLUSIONES

Los índices de diversidad genética de las poblaciones de D. citri fueron relativamente similares indicando una separación reciente de las poblaciones y aportan evidencias respecto a que las poblaciones de D. citri en el Valle del Cauca y el Quindío se encuentran en expansión, sin embargo no existe una diferenciación genética importante que pueda generar nuevos biotipos de esta especie capaces de adaptarse a nuevos ambientes y con la capacidad de albergar otras variantes de la bacteria Candidatus Liberibacter.

La variabilidad genética reducida disminuye la capacidad de una población para adaptarse a cambios ambientales generando la disminución de la adaptación de los individuos. Por lo tanto, el manejo de estas plagas requiere disminuir o mantener baja la variabilidad genética.

Es necesario hacer muestreos y monitoreo por todo el país, en general, y en particular en las localidades cítricas de Colombia, para poder tener así un panorama completo de la genética poblacional de la especie.

AGRADECIMIENTOS

Los autores agradecen el apoyo científico del Laboratorio de Biología Molecular de la Universidad Nacional de Colombia, sede Palmira, la asesoría y apoyo estadístico del Departamento de Ciencias Agropecuarias de la Universidad del Cauca y la financiación del proyecto e intercambio de conocimientos al grupo de Investigación de Diversidad Biológica de la Universidad Nacional de Colombia, sede Palmira.

REFERENCIAS

[1] CERMELI, M., MORALES, P., y GODOY, F. Presencia del psílido asiático de los cítricos Diaphorina citri Kuwayama (Hemiptera: Psyllidae) en Venezuela. Boletín Entomológico Venezuela, 15(2), 2000, p. 235-243.         [ Links ]

[2] REYES, R. Chicharrita (Diaphorina citri) amenaza para la citricultura en el Salvador. San Salvador (El Salvador): Programa Nacional de Frutas del Salvador, 2006, 4 p.         [ Links ]

[3] LANTERI, A., LOIÁCONO, M.S. y MARGARÍA, C. Aportes de La Biología Molecular a la Conservación de los Insectos. Red Iberoamericana de Biogeografía y Entomología Sistemática, 2, 2002, p. 207-220.         [ Links ]

[4] GARCÍA, C.Diaphorina citri Kuwayama (Hemíptera: Psyllidae), vector de la bacteria que causa el Huanglongbing (HLB-Greening). Buenos Aires (Argentina): Ministerio de la Producción, Servicio Nacional de Sanidad y Calidad Agroalimentaria, 2009, 18 p.         [ Links ]

[5] HALBERT, S.E. and MANJUNATH, K.L. Asian citrus Psyllids (Sternorrhyncha: Psyllidae) and greening disease of citrus: a literature review and assessment of risk in Florida. Florida Entomologist, 87(3), 2004, p. 330-353.         [ Links ]

[6] ALEMÁN, J., HEYKER, B. y RAVELO, J. Diaphorina citri y la enfermedad Huanglongbing: una combinación destructiva para la producción citrícola. Revista Protección Vegetal, 2(3), 2007, p. 154-165.         [ Links ]

[7] HALL, D.G. Biology, history and world status of Diaphorina citri. I Taller internacional sobre Huanglongbing de los citricos (Candidatus Liberibacter spp) y el psílido asiático de los cítricos (Diaphorina citri), Hermosillo (México): 2008, 11 p.         [ Links ]

[8] INSTITUTO COLOMBIANO AGROPECUARIO (ICA). Situación actual de HLB (huanglonbing) y su vector el psílido asiático de los cítricos (Diaphorina citri kuwayama) en Colombia. Boletín Epidemiológico, Bogotá (Colombia): 2010, 10 p.         [ Links ]

[9] TABERLET, P. The Use of Mitochondrial DNA Control Region Sequencing in Conservation Genetics. En: Molecular Genetic Approaches in Conservation. New York (USA): Smith, T. B. and Wayne, R. K. (Eds), Oxford University Press, 1996, p. 125-148.         [ Links ]

[10] ALJANABI, M. and MARTÍNEZ, I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Research, 25(22), 1997, p. 4692-4693.         [ Links ]
[11] ZHANG, D. and HEWITT, G.M. Assessment of the universality and utility of a set of conserved mitochondrial COI primers in insects. Insect Molecular Biology, 6(2), 1997, p. 143-150.         [ Links ]

[12] NEI, M. Molecular evolutionary genetics. New York, (USA): Columbia University Press, 1987, 512 p.         [ Links ]

[13] NEI, M. and LI, W.H. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences USA, 76, 1979, p. 5269-5273.         [ Links ]

[14] NEI, M., and TAJIMA, F. DNA polymorphism detectable by restriction endonucleases. Genetics, 97, 1981, p. 145-163.         [ Links ]

[15] NEIGEL, J.E., and AVISE, J.C. Application of a random walk model to geographic distributions of animal mitochondrial DNA variation. Genetics, 135, 1993, p. 1209-1220.         [ Links ]

[16] SAÉNZ, C.V. Genética de poblaciones y prevalencia del hongo patógeno Batrachochytrium dendrobatidis en anfibios de Papallacta y Guamaní, Ecuador [Tesis Ciencias Biológicas]. Quito (Ecuador): Pontificia Universidad Católica del Ecuador, Escuela de Ciencias Biológicas , 2011, 98 p.         [ Links ]

[17] BROWN, B.L., EPIFANIO, J.M., SMOUSE, P.E. and KOBAK, C.J. Temporal stability of mtDNA haplotype frequencies in American shad stocks: to pool or not to pool across years. Canadian Journal of Fisheries and Aquatic Sciences, 53(10), 1996, p. 2274-2283.         [ Links ]

[18] SARTORI, G.A. Diversidae genética de Diaphorina citri Kuwayama, (Hemiptera: Psyllidae) e caracterização molecular das linhagens de Wolbachia associadas [Tese de doutorado Entomología]. Sao Paulo (Brasil). Universidade de São Paulo, Escola Superior de Agricultura "Luiz de Queiroz", 2011,127 p.         [ Links ]

[19] BASS, C. and FIELD, L.M. Gene amplification and insecticide resistance. Pest Management Science, 67(8), 2011, p. 886-890.         [ Links ]

[20] MURRAY, D. and PASHLEY, D.P. Molecular phylogenetics and evolutionary history of the neotropical Satyrine Subtribe Euptyciina (Nymphalidae: Satyrinae). Molecular Phylogenetics and Evolution, 34, 2005, p. 67-80.         [ Links ]

[21] HAMILTON, M.B. Populations genetics. 1 ed. Chichester (England): Wiley-Blackwell, Oxford, 2009, 424 p.         [ Links ]

[22] BARÓN, J.D. Geografía económica de los Andes occidentales de Colombia. Bogotá (Colombia): Banco de la Republica, 2009, 95 p.         [ Links ]

[23] PEÑALOZA, M.C. y DÍAS, R.C. Así se maneja y controla el picudo de los cítricos Compsus sp. Bogotá (Colombia): Instituto Colombiano Agropecuario (ICA), Produmedios, 2004, 30 p.         [ Links ]

[24] SARMIENTO, L. Variabilidad genética de poblaciones mexicanas de la cochinilla rosada del hibisco, Maconellicoccus hirsutus (Green) (Hemiptera: Pseudococcidae) [Tesis Maestro em Ciencias en Biotecnología Genómica]. Tamaulipas (México): Instituto Politécnico Nacional, 2008, 61 p.         [ Links ]