Caracterización molecular de líneas de Bombyx mori (Lepidoptera: Bombycidae) mediante AFLP

Molecular characterization of Bombyx mori (Lepidoptera: Bombycidae) strains by AFLP

Resumen. Veintitrés líneas de gusano de seda Bombyx mori L, (12 líneas de origen japonés y 11 líneas de origen chino), pertenecientes al banco de germoplasma de gusano de seda del CDTS (Centro de desarrollo tecnológico de la Sericultura) se analizaron utilizando la técnica de Polimorfismo en la Longitud de los Fragmentos Amplificados (AFLP). Se emplearon cuatro combinaciones de cebadores. Las cuatro combinaciones ensayadas, arrojaron 247 bandas con un polimorfismo promedio de 54,6%, las cuales después de ser corregidas según el criterio de Lynch y Milligan resultaron en 168 bandas con un polimorfismo de 90,4%. El dendrograma construido utilizando el algori^Tmo UPGMA para los datos de disimilitud, reveló dos grupos muy bien definidos que corresponden a las razas geográficas japonesa y china. El promedio de heterocigocidad fue mayor en las líneas de la raza japonesa, observándose una dispersión mayor en el análisis de UPGMA y de los componentes principales. El estudio de la estructura poblacional, determina que aunque las dos razas no son muy diferentes entre sí, existe un componente de varianza mayor entre los grupos establecidos que dentro de éstos, lo que muestra un grado significativo de diferenciación genética.

Palabras clave: Sericultura. Estructura poblacional. Marcadores moleculares.

Abstract. Twenty-three silkworm Bombyx mori L. strains (twelve Japanese and eleven Chinese), belonging to the CDTS`s (Technological Development Center for Sericulture) silkworm germplasm bank, were analyzed using the Amplified Fragment Length Polymorphism (AFLP). Four combinations of primers were used. The assayed combinations threw 247 bands with an average polymorphism of 54.6%, these bands were corrected using Lynch and Milligan´s method, resulting in 168 bands with a polymorphism of 90,4%. The dendrogram was built using the UPGMA algorithm for the dissimilarity data, revealing two very well defined groups that correspond to the Japanese and Chinese geographical races. The heterocigocity average was higher in the Japanese race lines, showing a bigger dispersion in the UPGMA and the main components analysis. The population structure study, determines that although the two races are not very different to each other, it exists a component of greater variance, between the established groups than inside them; which shows a significant grade of genetic differentiation.

Keywords: Sericulture. Population structure. Molecular markers.

Introducción

Dentro de las diversas variedades de polillas que secretan fibras de seda, se incluyen Bombyx mori de la familia Bombycidae y las especies silvestres de la familia Saturniidae: Antherae mylitta, Antherae pernyi, Antherae assama, Antherae yamamai. Con una historia de más de 5000 años de domesticación, el gusano de seda, B mori, es la especie más importante a nivel agronómico como productor de seda; es el lepidóptero mejor estudiado debido a su rico repertorio de mutaciones bien caracterizadas que afectan virtualmente cada caracter de la morfología del organismo como también de su comportamiento y desarrollo. Su genoma haploide, con un tamaño de 560 Mb, es 3,8 veces más grande que el de Drosophila. melanogaster (Diptera: Drosophilidae) y se distribuye en 28 cromosomas holocéntricos de tamaño muy reducido (Nagaraju 2002).

El “stock” genético del gusano de seda domesticado lo componen un gran número de ecotipos y líneas autocruzadas artificialmente que se distribuyen en regiones templadas (líneas que presentan diapausa) y tropicales (líneas sin diapausa). Estas variedades difieren, entre otros, en los caracteres cualitativos y cuantitativos que afectan la calidad de la seda. Las líneas sin diapausa son pobres productoras de seda, presentan crecimiento rápido, peso corporal bajo y fibras de seda gruesas y cortas, son resistentes y se sabe que se pueden desarrollar aun en condiciones desfavorables como las encontradas en regiones tropicales, temperaturas altas y humedad alta propicias para condiciones de sanidad baja y el desarrollo de pestes y enfermedades. Las líneas tropicales de gusano de seda, sin embargo, presentan un gran número de crías al año; es decir, son polivoltinas, lo que de alguna manera compensa su productividad baja. Las líneas con diapausa, producen mayor cantidad de seda y de mejor calidad, presentan una vida larval larga, peso corporal alto y fibras de seda largas y delgadas aunque sólo se pueden obtener una o máximo dos crías al año (uni o bivoltinas). En estas ultimas, a pesar de ser mejores productoras de seda, su cría en regiones tropicales no es viable dada su alta susceptibilidad a pestes y malas condiciones de cría (Sharma et al. 1990).

En el gusano de seda, se han empleado muchas técnicas para identificar el polimorfismo a nivel molecular, como el uso de marcadores enzimáticos, estudiando amilasas de hemolinfa en líneas de gusanos tropicales y de regiones templadas (Nagaraju 2000). Los análisis permitieron identificar polimorfismo electroforético de las amilasas anódicas como cuatro bandas en las líneas sin diapausa y cinco bandas en las líneas con diapausa, la banda extra en las líneas de regiones templadas corresponde a una banda de migración lenta. Cuando la separación es catódica, se identificó una sola banda en las líneas sin diapausa y ninguna en las líneas con diapausa. En este estudio también se pudo determinar que las líneas de diapausa excretan cinco veces más almidón que las líneas sin diapausa (Nagaraju 2000). Los estudios moleculares de ADN basados tanto en marcadores que usan PCR como en aquellos que no, analizando trece líneas de gusano de seda seis con diapausa y siete sin diapausa, permitieron la separación de los individuos claramente en dos ramas de un dendrograma, aquella de los individuos con diapausa y aquellos que no presentaban diapausa, al igual que permitió la identificación de marcadores exclusivos para cada uno de los dos grupos (Nagaraju et al. 1995). El análisis de hibridización con sonda Bkm-2 (GATA₆₆TA), en las mismas trece líneas de gusano de seda, separó claramente los dos grupos y la identificación de bandas especificas para cada uno, además de la identificación de bandas específicas de sexo (Nagaraju et al. 1995). Damodar et al. (1999) informaron sobre el análisis de 28 loci microsatélites (SSR) que permitió la identificación de un marcador, sat 211, con alelos específicos para líneas con y sin diapausa, al igual que la identificación de un marcador muy polimórfico, sat 2763, con un total de 17 alelos y con una heterocigocidad elevada (90%) que permitiría fácilmente identificar las relaciones de parentesco y geográficas, muy útil, por ejemplo, en el análisis de variedades silvestres o de especies de gusano en peligro de extinción. Análisis similares sobre las mismas líneas, usando otras técnicas, han arrojado resultados parecidos (Reddy et al. 1999). Un estudio comparativo de técnicas de multilocus arrojó resultados importantes a la hora de escoger un sistema de marcadores moleculares para un trabajo en gusano de seda y aunque los agrupamientos fueron similares y se explicaban por la presencia de diapausa, origen geográfico y relaciones de parentesco de los individuos, estos agrupamientos diferían levemente en el tipo y grado de polimorfismo detectado (Nagaraju et al. 2001). El presente estudio tuvo como objetivos analizar la características poblacionales y grado de diversificación presente entre las líneas de gusano de seda B. mori mantenidas en el CDTS.

Materiales y Métodos

El estudio se desarrolló sobre 12 líneas de origen Japonés (KOI, K02, KO5, K10, K20, K30, K40, K522, SG3, SG2, NG, KNA) y 11 líneas de origen Chino (CA, CBS, CC, CJ, CGS, CHS, CLS, CTS, SC1, SC2, SC3), todas ellas pertenecientes al Centro de Desarrollo Tecnológico de la Sericultura (CDTS) mantenidas en la finca “El Pilamo”, Risaralda.

Extracción de ADN

El ADN genómico fue extraído de la porción posterior de la glándula serigena (Suzuki et al.1972). Se tomaron 0,3 g de glándula de seda de cada una de las líneas a analizar, manteniéndolas refrigeradas mientras se procesaban. El tejido se maceró en un tubo de microcentrífuga con amortiguador de extracción (NaCl 50mM, EDTA 50mM, SDS 1%) a 50°C. A cada muestra se agregaron 10 µl de proteinasa K (100 µg/ml). Las muestras se maceraron periódicamente hasta no observar fragmentos, durante la incubación en el baño. Se llevó a cabo extracción con un volumen de fenol-cloroformo-alcohol isoamilico (25-24-1) durante 5 min; posteriormente se centrifugó a 14.000 rpm durante 5 minutos y se recuperó la fase acuosa en otro tubo al que se le agregó un volumen de cloroformo-alcohol isoamilico (24-1), extrayéndose durante 5 min. Después de 5 min de centrifugación se tomó la fase acuosa y se le agregaron dos volúmenes de etanol absoluto y 1/10 de acetato de sodio 3M para precipitar los ácidos nucleicos. Posteriormente se resuspendió en amortiguador TE (10mM de Tris-HCl, 1mM de EDTA, pH 8,0) y se incubó a 37°C durante 1 h después de la adición de ARNasa (100 µg/ml). El ADN genómico fue resuspendido y cuantificado por espectrofotometría y observado por electroforesis en gel de agarosa al 0,8%.

AFLP y selección de iniciadores

Las reacciones de AFLP (Vos et al. 1995) se llevaron a cabo usando el estuche comercial ¨AFLP Analysis system I¨ No catálogo 10544-013 de la casa comercial INVITROGENE, siguiendo las recomendaciones del fabricante. La amplificación selectiva por la reacción en cadena de la polimerasa (PCR) se realizó usando una extensión de tres nucleótidos para ambos iniciadores. Las cuatro combinaciones de iniciadores empleadas, se escogieron por producir un número alto de polimorfismos no ambiguos, a partir del ensayo de 10 combinaciones sobre cuatro líneas de gusano. Los productos de PCR se mezclaron con 10 µl de solución de carga para secuenciación (98% formamida, 10mM EDTA, 0,025% xylen cyanol, 0,025% azul de bromofenol), se calentaron por 3 min a 90°C y rápidamente se enfriaron en hielo. 5 µl de la mezcla se agregaron en el gel de secuencia de poliacrilamida al 6% y se sometieron a 90 W durante 90 min. Después de terminada la corrida, el resultado se reveló mediante tinción con plata del estuche comercial “Silver sequence^Tm DNA sequencing system” de la casa comercial Promega No catálogo Q-4130, siguiendo las recomendaciones del fabricante. El registro permanente de estos resultados se realizó mediante fotografía digital transfiriendo las imágenes a formato TIFF con la cámara Nikon D´100.

Análisis estadístico

Similaridad y análisis de agrupamiento

La identificación de las relaciones para las 23 líneas de gusano analizadas se basó en las matrices construidas con los índices de similaridad (Nei y Li 1979). Los valores de disimilaridad se analizaron con el método de los pares de grupos no ponderados con media ari^Tmética (UPGMA), y se construyó un dendrograma. Un análisis de componentes principales (Gower 1966) también se llevó a cabo para mostrar la distribución de las relaciones en tres dimensiones. Todos estos análisis se hicieron usando NTSYS-pc ver.2.01 (Rohlf 1998).

Estadística descriptiva o diversidad. Antes de correr cada uno de los análisis el supuesto de homoscedasticidad (homogeneidad de varianzas) fue revisado mediante una prueba de Bartlett (Sokal y Rohlf 1981). Los programas POPGENE, ver 1.32 (Yeh et al. 1997) y TFPGA (Tools for Populaton Genetic Analisys), ver.1.3 (Miller 1999) se usaron para calcular, bajo el supuesto de equilibrio de Hardy-Weinberg, los siguientes parámetros: (i) frecuencia del alelo recesivo (0) y consecuentemente del dominante, con el método de la expansión de Taylor (Lynch y Milligan 1994), (ii) porcentaje de loci polimórficos, (iii) promedio de heterocigocidad no sesgada de todos los loci (Nei 1978), prueba de neutralidad, prueba exacta de diferenciación poblacional (Raymond y Rousset 1995). Estos análisis se realizaron tanto con las muestras crudas, como con las cuales se les había aplicado una restricción al análisis para evitar el sesgo en bandas que presentasen una frecuencia mayor a 1-(3/N) (Lynch y Milligan 1994) datos no mostrados

Estructura poblacional. El programa AMOVA-PREP (Miller 1998) se usó para transformar el conjunto de datos dominantes para el programa de Análisis Molecular de Varianza (AMOVA, ver.1.55.; Excofier et al. 1992). Los datos se analizaron como dos poblaciones y además como dos sub-poblaciones dentro de cada una de las poblaciones, para determinar el componente aportado por cada una de las subdivisiones a los componentes de varianza total. Los análisis se realizaron tomando tanto los datos crudos como los datos corregidos con el factor 1-(3/N).

Distancia genética. La distancia genética se evaluó sobre las poblaciones y sub-poblaciones con el programa TFPGA (Tools for Population Genetic Analisys), ver.1.3 (Miller 1999), empleando el cálculo descrito por Nei (1978).

Resultados

Similaridad y análisis de agrupamiento

La organización de los datos se llevó a cabo definiendo subpoblaciones dentro de cada una de las razas determinado a través de la relación en el dendrograma y analizando además solamente dos poblaciones correspondientes a cada una de las razas geográficas. La significación en la escogencia de estos grupos se comprobó estadísticamente con la prueba exacta de diferenciación de poblaciones (Raymond y Rousset 1995), encontrándose que el nivel de estructuración correspondiente a subpoblaciones no era significativo, pero el nivel de población si resultó serlo. Los análisis individuales de cada una de las combinaciones usadas arrojaron resultados similares y presentaron índices de correlación entre los datos de las diferentes matrices superiores al 90%, lo que significa que el número de combinaciones empleadas es suficiente para explorar la diversidad presente en las líneas de gusano de seda. El nivel de polimorfismo para las combinaciones estuvo en el rango 65,6-75,8%, aunque se generó un gran número de fragmentos de restricción (247 bandas), a partir de las muestras individuales con la técnica de AFLP (Tabla 1). Solamente aquellos fragmentos que fueron polimórficos y se registraron sin ambigüedad se usaron para el análisis, esta corrección se hizo según el criterio establecido por Lynch y Milligan en 1994 analizándose finalmente 168 bandas de los datos totales con un promedio de polimorfismo de 90,4% (Tabla 2). El resto de los análisis se llevaron a cabo solamente para los datos corregidos.

]]>