Técnicas moleculares para la detección de cianobacterias en los embalses Riogrande II y La Fe, ColombiaResúmenes
In lentic water bodies as reservoirs occur eutrophication processes, originated mainly from human activities (i.e. agriculture, animal exploitation). This influx of nutrients in aquatic ecosystems could promote blooms of potentially toxic cyanobacteria. The purpose of this work is to detect the presence of cyanobacteria strains in water samples, using molecular techniques to help in preventive management of reservoirs dedicated to water purification. We used two molecular techniques to detect genes implied with the synthesis of hepatotoxic microcystins from potentially toxic cyanobacteria strains, and to evaluate the molecular diversity of cyanobacteria in water samples from two high-mountain reservoirs used for purification of drinking water for the metropolitan area of Medellin, Colombia. Between 2010-2011 collections of 12 water samples were taken and DNA extraction together with PCR and DGGE analyses where carried out. We amplified 22 sequences between 250-300bp of the genes <span name="style_italic">mcy</span>A and <span name="style_italic">mcy</span>E, and these sequences were related with several strains and cyanobacteria genera accessions from NCBI-GenBank databases. Moreover, sequence amplifications of the 16S small ribosomal RNA subunit - 16S rRNA- between 400-800bp were also performed in order to use them for the DGGE technique. The amplification products of DGGE were set in polyacrilamide gel with posterior denaturing electrophoresis, and the scanned images of the gel bands were analysed with the software <span name="style_italic">GelCompar II</span>. For Riogrande II and La Fe reservoirs we found 35 and 30 different DGGE bands, respectively, as a measurement of molecular diversity in these artificial ecosystems. Here, we demonstrated the utility of two molecular techniques for the detection of genes associated with toxicity and molecular diversity of cyanobacteria in reservoirs destined for drinking water in urban centers. We recommend strongly following with periodically molecular biology studies in these ecosystems combined with limnological and ecological data, as new tools for management of plants of water purification and for the prevention of potentially public health concerns. Rev. Biol. Trop. 62 (1): 403-419. Epub 2014 March 01.. For Riogrande II and La Fe reservoirs we found 35 and 30 different DGGE bands, respectively, as a measurement of molecular diversity in these artificial ecosystems. Here, we demonstrated the utility of two molecular techniques for the detection of genes associated with toxicity and molecular diversity of cyanobacteria in reservoirs destined for drinking water in urban centers. We recommend strongly following with periodically molecular biology studies in these ecosystems combined with limnological and ecological data, as new tools for management of plants of water purification and for the prevention of potentially public health concerns. Rev. Biol. Trop. 62 (1): 403-419. Epub 2014 March 01.
eutrofización; cianobacterias; toxicidad; microcistinas; diversidad; marcadores moleculares
En embalses, la eutrofización es consecuencia de procesos naturales y de actividades humanas, lo cual puede facilitar la aparición de afloramientos de cianobacterias potencialmente tóxicas. En este trabajo, se utilizaron dos técnicas moleculares en la detección de genes presentes en cepas potencialmente tóxicas de cianobacterias y evaluaciones de diversidad molecular de cianobacterias en los embalses para potabilización de agua, Riogrande II y La Fe, Colombia. Entre 2010-2011, se tomaron 12 muestras de agua en ambos embalses y se realizaron extracciones de ADN para un análisis de marcadores moleculares mediante PCR y DGGE. Se amplificaron secuencias entre 250-300pb de los genes <span name="style_italic">mcy</span>A y <span name="style_italic">mcy</span>E, implicados en la toxicidad de cepas de cianobacterias. Asimismo, se amplificaron secuencias de la región 16S del ARN ribosomal (422pb), para la técnica DGGE. Se corrieron geles de poliacrilamida en gradientes de desnaturalización, se realizó agrupamiento genético (<span name="style_italic">UPGMA</span>), y se separaron por patrones de bandeo las muestras de cada embalse evaluado. Se demuestra la utilidad de las técnicas moleculares en estudios relacionados con la búsqueda de genes asociados con toxicidad y diversidad molecular de cianobacterias en muestras de agua provenientes de embalses de agua con fines de potabilización para centros urbanos.), y se separaron por patrones de bandeo las muestras de cada embalse evaluado. Se demuestra la utilidad de las técnicas moleculares en estudios relacionados con la búsqueda de genes asociados con toxicidad y diversidad molecular de cianobacterias en muestras de agua provenientes de embalses de agua con fines de potabilización para centros urbanos.
eutrofización; cianobacterias; toxicidad; microcistinas; diversidad; marcadores moleculares
Técnicas moleculares para la detección de cianobacterias en los embalses Riogrande II y La Fe, Colombia
Molecular techniques for cyanobacteria detection at Riogrande II and la Fe water reservoirs, Colombia
Julio César Hurtado-Alarcón1* & Jaime Polanía-Vorenberg2*
*Dirección para correspondencia:
Abstract
In lentic water bodies as reservoirs occur eutrophication processes, originated mainly from human activities (i.e. agriculture, animal exploitation). This influx of nutrients in aquatic ecosystems could promote blooms of potentially toxic cyanobacteria. The purpose of this work is to detect the presence of cyanobacteria strains in water samples, using molecular techniques to help in preventive management of reservoirs dedicated to water purification. We used two molecular techniques to detect genes implied with the synthesis of hepatotoxic microcystins from potentially toxic cyanobacteria strains, and to evaluate the molecular diversity of cyanobacteria in water samples from two high-mountain reservoirs used for purification of drinking water for the metropolitan area of Medellin, Colombia. Between 2010-2011 collections of 12 water samples were taken and DNA extraction together with PCR and DGGE analyses where carried out. We amplified 22 sequences between 250-300bp of the genes mcyA and mcyE, and these sequences were related with several strains and cyanobacteria genera accessions from NCBI-GenBank databases. Moreover, sequence amplifications of the 16S small ribosomal RNA subunit - 16S rRNA- between 400-800bp were also performed in order to use them for the DGGE technique. The amplification products of DGGE were set in polyacrilamide gel with posterior denaturing electrophoresis, and the scanned images of the gel bands were analysed with the software GelCompar II. For Riogrande II and La Fe reservoirs we found 35 and 30 different DGGE bands, respectively, as a measurement of molecular diversity in these artificial ecosystems. Here, we demonstrated the utility of two molecular techniques for the detection of genes associated with toxicity and molecular diversity of cyanobacteria in reservoirs destined for drinking water in urban centers. We recommend strongly following with periodically molecular biology studies in these ecosystems combined with limnological and ecological data, as new tools for management of plants of water purification and for the prevention of potentially public health concerns. Rev. Biol. Trop. 62 (1): 403-419. Epub 2014 March 01.
Key words: eutrophication, cyanobacteria, toxicity, microcystins, diversity, molecular markers.
Resumen
En embalses, la eutrofización es consecuencia de procesos naturales y de actividades humanas, lo cual puede facilitar la aparición de afloramientos de cianobacterias potencialmente tóxicas. En este trabajo, se utilizaron dos técnicas moleculares en la detección de genes presentes en cepas potencialmente tóxicas de cianobacterias y evaluaciones de diversidad molecular de cianobacterias en los embalses para potabilización de agua, Riogrande II y La Fe, Colombia. Entre 2010-2011, se tomaron 12 muestras de agua en ambos embalses y se realizaron extracciones de ADN para un análisis de marcadores moleculares mediante PCR y DGGE. Se amplificaron secuencias entre 250-300pb de los genes mcyA y mcyE, implicados en la toxicidad de cepas de cianobacterias. Asimismo, se amplificaron secuencias de la región 16S del ARN ribosomal (422pb), para la técnica DGGE. Se corrieron geles de poliacrilamida en gradientes de desnaturalización, se realizó agrupamiento genético (UPGMA), y se separaron por patrones de bandeo las muestras de cada embalse evaluado. Se demuestra la utilidad de las técnicas moleculares en estudios relacionados con la búsqueda de genes asociados con toxicidad y diversidad molecular de cianobacterias en muestras de agua provenientes de embalses de agua con fines de potabilización para centros urbanos.
Palabras clave: eutrofización, cianobacterias, toxicidad, microcistinas, diversidad, marcadores moleculares.
Medellín, en el Departamento de Antioquia, es la segunda ciudad más importante de Colombia. Su área metropolitana se encuentra distribuida a lo largo y ancho del Valle de Aburrá, y está conformada por los municipios de Medellín, Barbosa, Bello, Caldas, Copacabana, Envigado, Girardota, Itagüí, La Estrella y Sabaneta. Para el censo poblacional de 2005 en Colombia, se estimaron más de dos millones de habitantes en Medellín y cerca de un millón de habitantes en el resto de municipios del área metropolitana (DANE, 2005). Asimismo, algunos de estos municipios comparten la utilización de embalses para obtención de agua potable y energía eléctrica. En particular, parte del agua que abastece el acueducto de municipios del Valle de Aburrá como Medellín y Envigado, proviene principalmente de los embalses Riogrande II y La Fe, en los municipios de Don Matías y El Retiro, localizados al norte y al oriente de Medellín en el departamento de Antioquia, respectivamente.
Estudios anteriores han mostrado que embalses de características similares presentan procesos de eutrofización artificial, como consecuencia de actividades humanas (Ramírez, Gutiérrez & Vargas, 2005). Tales actividades, particularmente las agrícolas, han incrementado la concentración de nutrientes y estimulado la productividad primaria en los cuerpos de agua continentales (Havens, James, East & Smith, 2003). Este proceso de eutrofización artificial se manifiesta en floraciones visibles de plantas flotantes, agregaciones de macrófitas bénticas y, en algunos casos, en florecimientos de cianobacterias con potenciales efectos directos e indirectos sobre diferentes organismos, incluidos los seres humanos (Hitzfeld, Hoger & Dietrich, 2000). Además de los nutrientes, la temperatura, el pH o las altas intensidades lumínicas también pueden influir sobre los florecimientos de cianobacterias (Runnegar, Kong & Berndt, 1993; Dawson, 1998), con implicaciones directas e indirectas sobre los ecosistemas, puesto que estos microorganismos tienen la potencialidad de producir toxinas como estrategia de defensa. De esta forma, los peces, las aves y los seres humanos, entre otras formas de vida, pueden ser afectados por toxinas cianobacterianas -cianotoxinas- (Tencalla & Dietrich, 1997; Azevedo et al., 2002). Autores como Azevedo & Carmouze (1994), Codd, Steffensen, Burch, & Baker (1994), Rodger, Turnbull, Edwards & Codd (1994), entre otros, han registrado altas mortalidades de peces y animales domésticos expuestos a cianotoxinas. Por otra parte, en Caruarú, noreste del Brasil, fallecieron 76 pacientes de una clínica de hemodiálisis por contaminación del agua con la cianotoxina hepatotóxica microcistina-LR (Jochimsen et al., 1998; Domingos, Rubim, Molica, Azevedo & Carmichael, 1999; Carmichael et al., 2001).
Para la detección de cianobacterias, el uso de técnicas de microscopía basada en su morfología ha sido una herramienta de uso común. No obstante, cianobacterias potencialmente tóxicas no tienen características morfométricas únicas que ayuden en su identificación (Ouellete, Handy & Wilhelm, 2006). En algunos casos, las cepas tóxicas de géneros cianobacterianos se diferencian de las no tóxicas porque presentan genes que codifican para la síntesis de péptidos potencialmente tóxicos, los cuales pueden detectarse por métodos moleculares (Hotto, Satchwell & Boyer, 2007). Uno de ellos, la Reacción en Cadena de la Polimerasa (PCR, Polymerase Chain Reaction), permite detectar genes del cluster mcy, algunos de los cuales codifican para proteínas del péptido hepatotóxico microcistina sintetasa (Hisbergues, Christiansen, Rouhiainen, Sivonen & Börner, 2003). Dentro de este cluster, la detección de las regiones mcyA y mcyE ha sido de particular importancia, puesto que son importantes en la identificación de géneros de cianobacterias potencialmente tóxicas (Rantala et al., 2006; Hotto et al., 2007). Así mismo, se han llevado a cabo diferentes estudios para analizar la ocurrencia natural de cianobacterias en ecosistemas acuáticos mediante la amplificación de secuencias de la subunidad 16S del ADN ribosomal, lo cual es útil en el monitoreo de la distribución de especies (Matsunaga, Takeyama & Nakayama, 2001), en estudios de biodiversidad (Ouellete et al., 2006), y en estudios filogenéticos y evolutivos (Genuário, Silva-Stenico, Welker, Moraesc & Fiore, 2010).
Dado que la Organización Mundial de la Salud (WHO) permite un máximo de 1µg/L de toxinas como microcistina-LR en agua para potabilizar (WHO, 2004), es de vital importancia detectar la presencia de cianobacterias con potencial toxicidad en cuerpos de agua lénticos dedicados a la obtención de agua potable para millones de personas. Considerando que Correa (2008) identificó dos géneros diferentes de cianobacterias productoras de Microcistina-LR en el embalse Riogrande II, Colombia, mediante técnicas de microscopía tradicional, en este artículo se describe la detección molecular de cepas de cianobacterias con potencial toxicidad en los embalses Riogrande II y La Fe, Colombia, entre los años 2010-11, mediante la amplificación de los genes mcyA y mcyE. Así mismo, se amplificó un fragmento de la subunidad 16S del ARN ribosomal (rRNA) cianobacteriano, y se realizó Electroforesis en Gradientes de Desnaturalización en Geles de poliacrilamida (DGGE), en un intento para detectar diferencias o similitudes en los patrones de bandeo obtenidos por esta técnica molecular con relación a los ecosistemas evaluados, como un estudio complementario para el monitoreo de la presencia de géneros y cepas de cianobacterias potencialmente tóxicas en estos ecosistemas artificiales en Colombia.
Materiales y métodos
Área de estudio: Se obtuvieron muestras de agua de los embalses Riogrande II y La Fe, Colombia, en seis salidas de campo realizadas entre 2010-2011. El embalse Riogrande II se encuentra en la zona central del Departamento de Antioquia en jurisdicción de los municipios de San Pedro, Entrerríos, Belmira, Don Matías y Santa Rosa de Osos (Otaya, Vásquez & Bustamante, 2008). Se construyó entre 1986-1988, y está formado por la confluencia de los ríos Chico y Grande, y la quebrada Las Ánimas; se encuentra a una altura promedio de 2 270msnm, y presenta un área superficial de 12 000km2, un volumen de 220Mm3, una profundidad máxima de 42m, una profundidad promedio de 37.9m y una longitud máxima de 10km (Montoya & Ramírez, 2007). El agua que se capta es aprovechada por las Empresas Públicas de Medellín –EPM- para el abastecimiento del acueducto de la ciudad y algunos municipios del Valle de Aburrá, y también es utilizada para la generación de energía en las centrales hidroeléctricas Tasajera y Niquía, generando alrededor de 324MW (CORANTIOQUIA, 1997). Por su parte, el embalse La Fe fue llenado en 1 973, y cubre un área de 173km2. Se sitúa en el municipio de El Retiro, con alturas entre los 2 175-3 000msnm. Este embalse es utilizado principalmente para acueducto, pero una parte está destinada a recreación. De su volumen total (aprox. 15Mm3), 12Mm3 regulan el río Pantanillo, y 8m3/s están destinados al acueducto metropolitano a través de la planta de tratamiento de La Ayurá, municipio de Envigado, Antioquia (Abuchaibe, Agudelo & Sañudo 1988).
Recolecta de muestras de agua: Mediante redes de arrastre de fitoplancton de 35µm de diámetro de poro, se tomaron entre 100-300mL de agua en la zona fótica en los embalses Riogrande II y La Fe. Dentro de cada embalse se tomaron muestras de agua en diferentes estaciones de muestreo (cuatro estaciones en Riogrande II: 1, 2, 3, 8; y seis estaciones en La Fe: 1, 2, 3, 4, 6, 7), las cuales aparecen con sus nombres sintetizados y coordenadas geográficas en el cuadro 1.
Las muestras de agua se pasaron por filtros de celulosa de 0.5µm y se refrigeraron a 4°C hasta la noche o hasta el día siguiente en el laboratorio del grupo GAIA –Grupo de Investigación, Gestión y Modelación Ambiental–, en la Sede de Investigación Universitaria de la Universidad de Antioquia, SIU. Posteriormente, se transportaron hacia el Laboratorio de Biología Molecular y Celular de la Universidad Nacional de Colombia, Sede Medellín. Allí, las muestras se homogenizaron y de cada una de ellas se tomaron volúmenes de 50mL por duplicado, y se centrifugaron en volúmenes decrecientes de 50, 15, y 2mL entre 200-3 000g durante 10-15min. De esta forma, el precipitado depositado en los viales de menor volumen (2mL), contenía la materia prima para llevar a cabo las extracciones de ADN, las cuales se realizaron usando el estuche QIAamp DNA Mini Kit (Qiagen), de acuerdo con las recomendaciones del fabricante. La concentración y pureza del ADN extraído se evaluó mediante un espectrofotómetro Nanodrop (Thermoscientific, USA). La calidad de las extracciones de ADN se verificó en geles de agarosa al 1.0% teñidos con bromuro de etidio (EtBr) (10mg/mL), y fueron visualizadas en transiluminador de luz ultravioleta.
Amplificación por PCR (Polymerase Chain Reaction, Reacción en Cadena de la Polimerasa) de regiones de interés: Inicialmente, se evaluaron las secuencias cebadoras cpcB-cpcA (PCβ-F:
5›-GGCTGCTTGTTTACGCGACA-3›, y PCα-R: 5›-CCAGTACCACCAGCAACTAA- 3› (Neilan, Jacobs & Goodman, 1995), para amplificar el gen que codifica para ficocianina-c específica de cianobacterias (Cuadro 2), como control positivo de la presencia de cianobacterias en las muestras de ADN obtenidas de los embalses. Para la amplificación posterior de fragmentos del cluster de genes mcy y de la región 16S del ARN ribosomal, solo se utilizaron muestras de ADN positivas para ficocianina-c.
Evaluación de genes relacionados con toxicidad de cianobacterias: Para identificar regiones relacionadas con toxicidad de cianobacterias, se evaluaron los cebadores mcyA-Cd (1F5’- AAAATTAAAAGCCGTATCAAA-3’;1R-5’AAAAGTGTTTTATTAGCGGCTCAT) (Hisbergues et al., 2003), y para la posible identificación de los géneros Microcystis y Anabaena se evaluaron cebadores específicos (mcyE-F2: 5’-GAAATTTGTGTAGAAGGTGC-3’ para ambos géneros; mcyE-R8 5’- CAATGGGAGCATAACGAG-3’, Microcystis; mcyE-12R 5’-CAATCTCGGTATAGCGGC-3’, Anabaena) (Vaitomaa et al., 2003; Rantala et al., 2006; Hotto et al., 2007) (Cuadro 2).
El volumen final de cada mezcla de PCR fue de 25mL. Las amplificaciones por PCR de las regiones de interés se realizaron de acuerdo con los protocolos publicados por Hisbergues et al. (2003) (gen mcyA), Vaitomaa et al. (2003) (gen mcyE) y Dos Anjos et al. (2006) (operon ficocianina cianobacteriana, cpcB-cpcA), con respecto a concentraciones de reactivos y perfiles térmicos. Las reacciones de PCR se llevaron a cabo en un termociclador T3 000 (Biometra, Alemania) en microtubos de 200μL. Los productos de amplificación se visualizaron mediante electroforesis en geles de agarosa (1.5%) teñidos con EtBr (10mg/mL). Las condiciones de corrido fueron: 60-80min en solución 1X TBE (Tris-ácido bórico-EDTA, pH8.0), a un voltaje constante (80V) en cámaras de electroforesis horizontales (Bio-Rad, E.U.A.). Los geles de agarosa se observaron sobre una pantalla de luz ultravioleta y se fotografiaron con un equipo para captura de imágenes (Biometra, Alemania) para confirmar la efectividad de la amplificación. Para determinar los tamaños de los alelos en cada corrida se insertó un marcador molecular de secuencia conocida en un rango de 100-1 000 pares de bases (GeneRuler Low Range DNA Ladder, Fermentas, Canada) (Hisbergues et al., 2003, Dos Anjos et al., 2006).
Purificación de las muestras positivas y secuenciación: Las muestras que resultaron positivas para genes de toxicidad (mcyA, mcyE) se purificaron utilizando el estuche comercial QIAquick PCR Purification Kit (QIAGEN, Alemania), siguiendo las recomendaciones del fabricante, las amplificaciones se verificaron en geles de agarosa (1.5%) y finalmente se enviaron a secuenciar (Macrogen, Corea). Se realizaron comparaciones individuales de las secuencias obtenidas de los genes mcyA y mcyE en las bases de datos del NCBI (National Center for Biotechnology Information, Centro Nacional de Información en Biotecnología, E.U.A.), usando la herramienta para alineamientos locales BLAST (Basic Local Alignment Search Tool).
Análisis de secuencias para los genes mcyA y mcyE: Se utilizó el editor de alineamientos gratuito BioEdit (Biological Sequence Alignment Editor) versión 7.0.9 (Hall, T.A. 1999), el cual se puede descargar gratuitamente en el sitio URL http://www.mbio.ncsu.edu/bioedit/bioedit.html. Este programa permite obtener secuencias consenso a partir de los resultados de las secuenciaciones de las cadenas Forward y Reverse de cada una de las muestras positivas para los genes mcyA y mcyE. Con las secuencias consenso, se realizó un alineamiento múltiple con el algoritmo Clustal W incluido en el programa de uso libre BioEdit. Adicionalmente, las secuencias consenso obtenidas fueron alineadas y analizadas junto con secuencias publicadas en las bases de datos del NCBI de algunos géneros y cepas de cianobacterias (Microcystis sp., Nostoc sp., Anabaena sp., Radiocystis sp.) que portan los genes mcyA y mcyE, relacionados con la síntesis de péptidos con potencial hepatotoxicidad (Cuadro 3). La escogencia del modelo de sustitución nucleotídica y los análisis de agrupamiento de las secuencias se llevaron a cabo con el programa de uso libre MEGA5, descargado del sitio URL http://www.megasoftware.net/ (Tamura et al., 2011).
Amplificación de secuencias de la subunidad 16S del ARN ribosomal (rRNA) para la evaluación de diversidad de cepas de cianobacterias mediante la técnica molecular DGGE (Denaturing gradient gel electrophoresis - Electroforesis en gel con gradiente de desnaturalización): Se evaluaron las secuencias cebadoras descritas por Nübel, Garcia-Pichel & Muyzer (1997), para amplificar regiones del gen que codifica para la subunidad 16S rRNA de cianobacterias (Cuadro 3), de aproximadamente 422 pares de bases (bp) (cebadores 359F-781R), a partir del ADN extraído de las muestras de agua de los embalses Riogrande II y La Fe (Cuadro 4). Se utilizó el cebador Forward 359F modificado en el extremo 3´, mediante la adición de 40 desoxirribonucleótidos fosfato guanina-citosina (GC), para modificar el comportamiento desnaturalizante del fragmento de interés y mejorar la resolución de los fragmentos de tamaño similar pero con contenido diferente de GC (Muyzer, De Waal & Uitterlinden, 1993). La mezcla de PCR se preparó de acuerdo con lo publicado por Nübel et al. (1997). Se tomaron 25rm de cada primer, 25hmol dNTPs, 10mL buffer PCR, (100mM Tris-HCl [pH 9.0], 1.5mM MgCl, 500mM KCl), 0.5U Taq polimerasa (FERMENTAS UAB, Canadá) y entre 30-100ηg de ADN molde obtenido de las extracciones, en un volumen final de 30-50mL por reacción. El perfil térmico de la PCR incluyó un ciclo inicial de desnaturalización del ADN a 95°C, seguido de 35 ciclos de un minuto cada uno a 94°C, 60°C y 72°C. Por último, se incluyó un ciclo de extensión final a 72°C por 10 min.
Evaluación de la técnica DGGE: Las muestras que amplificaron con los cebadores 359FGC-781R se precipitaron con etanol de grado reactivo al 96%, se llevaron a una concentración entre 2 000-3 000ηg y se diluyeron en un volumen final de 10μL con agua ultrapura estéril. Con estos productos de PCR se llevaron a cabo electroforesis en geles de poliacrilamida, en el sistema D-Code Mutation Detection System (Bio-Rad, E.U.A.). Para ello, se prepararon geles de poliacrilamida al 8% y 1mm de grosor, en dos gradientes desnaturalizantes de urea y formamida (30-60%), de acuerdo con los protocolos descritos por Nübel et al. (1997). Los geles de poliacrilamida se sometieron a electroforesis en solución tampón TAE (40mM Tris-HCl [pH 8.3], 20mM ácido acético, 1mM EDTA), durante 3.5 horas a 200V, a una temperatura constante de 60°C. Los geles de poliacrilamida se tiñeron con nitrato de plata (AgNO3) y se sumergieron en solución reveladora (NaOH, formamida y agua destilada) para visualizar bandas de ADN, (Sanguinetti, Dias Neto & Simpson, 1994). Posteriormente, los geles se fijaron con solución Etanol-ácido acético 95:5, se escanearon y se digitalizaron para el análisis de imágenes usando el programa GelCompar II 6.0 (Applied Maths, Bélgica). Este programa permite analizar agrupaciones (clusters) evaluando algoritmos para determinación de similitud de bandas de amplificaciones y polimorfismos en evaluaciones de diversidad molecular (Rademaker & de Bruijin, 2004). Dicha metodología es útil para evaluar similitudes o diferencias en los patrones de bandeo de las amplificaciones de la región 16S del ARN ribosomal cianobacteriano, obtenidos en ecosistemas lénticos artificiales e independientes como lo son los embalses Riogrande II y La Fe.
Por medio del programa GelCompar II 6.0, se calculó la matriz de distancia basada en las bandas de DGGE usando el coeficiente de similitud de Dice (Nei & Li, 1979), el cual es comúnmente usado para comparar especies entre ecosistemas diferentes. Dos perfiles idénticos dan como resultado un valor igual al 100%, mientras que dos perfiles completamente diferentes dan un valor igual al 0%. Se realizó un análisis de agrupamiento mediante el método de promedio aritmético no ponderado (UPGMA -Unweighted Pair Group Method using Arithmetic averages), y los resultados fueron presentados finalmente como un dendrograma.
Resultados
Extracción y cuantificación de ADN a partir de las muestras de agua tomadas en los embalses Riogrande II y La Fe: El trabajo de extracción de ADN usando el estuche comercial permitió obtener ADN de buena calidad en las muestras de agua obtenidas de los embalses Riogrande II y La Fe, Colombia. Así mismo, la cuantificación del ADN usando espectrofotómetro arrojó lecturas entre 7.3- 188.3ηg/mL. A partir de las muestras cuantificadas, se usaron concentraciones de 20ηg de ADN por cada muestra para las amplificaciones por PCR de las regiones de interés.
Amplificaciones de genes del cluster mcy, relacionados con toxicidad de cianobacterias:
Amplificación de los genes mcyA y mcyE: Para los dos embalses evaluados, se obtuvieron amplificaciones en las reacciones de PCR de los genes mcyA y mcyE en las muestras de agua obtenidas. Las amplificaciones de parte del dominio de condensación (Cd) del gen mcyA, tuvieron un tamaño aproximado al publicado por Hisbergues et al. (2003) de 300bp, utilizando el par de cebadores mcyA-Cd (Hisbergues et al., 2003). Para el gen mcyE, solo se obtuvieron amplificaciones por PCR con los cebadores específicos para Microcystis sp. (mic-mcyE R8). El tamaño obtenido para esta región fue diferente a lo encontrado por otros autores (aprox. 300bp) (Vaitomaa et al., 2003). Para Anabaena sp. no se obtuvieron amplificaciones positivas al usar los cebadores Ana-mcyE 12R.
Comparación de las secuencias consenso obtenidas por secuenciación de las amplificaciones de los genes mcyA y mcyE, usando la herramienta de alineamiento local BLAST: Para el gen mcyA que codifica para microcistina sintetasa, se hallaron correspondencias en el GenBank con varios clones para los embalses evaluados. Para el embalse Riogrande II, las consensos analizados correspondieron con un clon de cianobacterias sin identificar (Uncultured cyanobacterium clone 03LO-OS4, Nº accesión GenBank EF424345.1, cobertura 99%) y con un clon no cultivado de Microcystis (Uncultured Microcystis sp. clone MB20, Nº accesión GenBank JN108773.1, cobertura 100%). Para el embalse La Fe, de cuatro consensos dos se relacionaron con clones no cultivados de cianobacterias sin identificar, y dos consensos se relacionaron con uno de los clones hallados en el embalse Riogrande II (Uncultured cyanobacterium clone 03LO-OS4, Nº accesión GenBank EF424345.1, cobertura 99%; Uncultured cyanobacterium clone 03LOOR9, Nº accesión GenBank EF424337.1, cobertura 100%).
Con relación al gen mcyE que codifica para microcistina sintetasa para ambos embalses, se hallaron correspondencias similares. Para los embalses Riogrande II y La Fe, todos los 13 consensos diferentes se relacionaron con las mismas cepas de Microcystis (Microcystis sp. CYN06, Nº accesión GenBank FJ393327.1, cobertura 99-100%), y en una medida similar con una cepa de Microcystis aeruginosa (CPCC 299, Nº accesión GenBank JN936964.1, cobertura 99-100%, identidad >98%). No obstante para el embalse La Fe, de dos secuencias consenso, una de ellas tuvo altos valores de identidad (98%) pero con un valor bajo de cobertura de 44%.
Análisis de agrupamiento para las secuencias consenso obtenidas: El programa BioEdit, permitió encontrar siete secuencias consenso para el gen mcyA y 13 secuencias consenso para el gen mcyE. El programa MEGA 5 permitió obtener matrices de agrupamiento y de distancias genéticas para estos genes. Los análisis de las secuencias consenso de los genes mcyA y mcyE se basaron en el método de Máxima Verosimilitud (Maximum-likelihood, ML) usando el modelo de sustitución nucleotídica Kimura 2P (K2), y se obtuvieron dendrogramas de agrupamiento (Fig. 1 y Fig. 2), con topologías soportadas por análisis de Bootstrapping (1 000 réplicas). Las secuencias que se obtuvieron de las bases de datos del NCBI de algunos géneros y cepas de cianobacterias (Microcystis sp., Nostoc sp., Anabaena sp., Radiocystis sp.), se mostraron como outgroups para las secuencias obtenidas en el presente estudio, debido a altos valores de remuestreo mediante bootstrapping hallados entre los clados obtenidos (99%), exceptuando una secuencia obtenida del embalse La Fe para el gen mcyA.
Para el gen mcyA se encontraron dos agrupaciones de secuencias, la primera que comprende ocho de las nueve secuencias de los embalses Riogrande II y La Fe, y la segunda agrupación la cual incluye seis secuencias del NCBI más una sola secuencia del embalse La Fe. Para este gen, los valores de remuestreo bootstrapping entre la mayoría de las secuencias de los embalses antioqueños fueron muy bajos (1-11% de las replicaciones del remuestreo). Los valores de bootstrapping para la segunda agrupación de secuencias (una secuencia del embalse La Fe y las secuencias obtenidas del NCBI), fueron mayores a los valores obtenidos para la primera agrupación (41-99% de las replicaciones del remuestreo), siendo el valor más bajo (41%) el obtenido al comparar cinco de las seis muestras del NCBI con la secuencia del embalse La Fe.
Para el gen mcyE, se hallaron varias agrupaciones dentro de las muestras de los embalses evaluados, separando claramente dichas muestras de las secuencias obtenidas del NCBI para el gen mcyE. Los valores obtenidos de remuestreo mediante bootstrapping fueron muy bajos para 11 de las 13 secuencias de los embalses Riogrande II y La Fe (0-4% de las replicaciones del remuestreo), y tres de las 13 muestras de los embalses antioqueños tuvieron valores altos de bootstrapping (58-99% de las replicaciones del remuestreo).
Evaluación de diversidad de cepas cianobacterianas en los embalses Riogrande II y La Fe:
Evaluación de la amplificación de secuencias de la subunidad 16S del ARN ribosomal (rRNA) y su evaluación en geles de poliacrilamida para la técnica DGGE: Para los embalses Riogrande II y La Fe, se obtuvieron amplificaciones de la región 359FGC-781F(a,b) del gen 16S del ARN ribosomal (rRNA), las cuales se corrieron en geles de poliacrilamida para la técnica DGGE. Las bandas obtenidas se pudieron visualizar mediante los protocolos de tinción con nitrato de plata (Fig. 3).
Análisis de los agrupamientos por patrones de bandeo de las reacciones obtenidas por amplificación y DGGE de la región 16S rRNA: Los geles de poliacrilamida digitalizados y analizados con el programa GelCompar II permitieron realizar agrupamientos de bandas comunes y no comunes en los embalses Riogrande II y La Fe. Se detectaron 35 bandas diferentes en las muestras del embalse Riogrande II, y 30 bandas diferentes para el embalse La Fe, con las cuales se realizó el dendrograma descripto de agrupaciones genéticas (Fig. 4). Se agruparon la mayoría de sitios de muestreo dentro de cada embalse a partir de los geles de poliacrilamida obtenidos. El dendrograma obtenido separa dos agrupaciones genéticas, diferenciando con claridad las muestras del embalse Riogrande II de las del embalse La Fe, exceptuando las muestras del sitio 7A. (Fig. 5).
Discusión
En este proyecto, la mayoría de cebadores evaluados para amplificar regiones de ADN de interés, permitieron obtener resultados positivos en las muestras de agua de los embalses Riogrande II y La Fe, para genes que codifican para ficocianina de cianobacterias y algunos péptidos asociados con toxicidad (mcyA, mcyE). Así mismo, se obtuvo éxito en la amplificación de regiones del gen que codifica para la subunidad 16S del ARN ribosomal (rRNA). Estos resultados son importantes debido a que son indicadores de la presencia de cepas de cianobacterias con potencial toxicidad en los embalses Riogrande II y La Fe, destinados para la potabilización de agua en centros urbanos de la ciudad de Medellín, Colombia.
La amplificación de fragmentos de la región cpc-IGS específica de cianobacterias en muestras de agua de embalses, puede dar una idea inicial de la presencia de cianobacterias. No obstante, su identificación en muestras de agua no permite determinar cuáles cepas o géneros de cianobacterias están presentes en las muestras evaluadas. Así mismo, con la amplificación de esta región no se logra diferenciar cepas tóxicas (Bittencourt-Oliveira et al., 2001). Sin embargo, los cebadores evaluados para los genes relacionados con toxicidad mcyA y mcyE, ayudaron a detectar regiones relacionadas con cepas y géneros cianobacterianos como Microcystis con potencial toxicidad.
Los resultados de las amplificaciones inespecíficas de regiones del gen mcyA para muestras de agua de los embalses Riogrande II y La Fe eran de esperarse. La amplificación del gen mcyA solo apunta a describir la presencia de cepas potencialmente productoras de microcistinas sin discriminación de cepas o especies, lo cual ha sido relacionado en estudios anteriores (Hisbergues et al., 2003; Hotto et al., 2007). Esto es diferente de los resultados que pueden arrojar amplificaciones del gen mcyE, el cual usa cebadores específicos de cepas y especies (Vaitomaa et al., 2003; Rantala et al., 2006). De esta forma, para la región mcyE, también se encontró el género cianobacteriano Microcystis sp. a partir de la información de secuenciación contrastada en las bases de datos geneticos del NCBI, con representantes de dos cepas diferentes de M. aeruginosa para los dos embalses evaluados.
De acuerdo con los dendrogramas obtenidos para ambos genes, se indica un alto grado de similitud entre las mismas secuencias y poca divergencia genética, para las secuencias de los dos genes evaluados en los embalses antioqueños (mcyA, mcyE). Estos datos indican a su vez alta similitud genética entre la mayoría de secuencias analizadas en este trabajo, en dos embalses antioqueños separados geográficamente por menos de 50km linealmente y por 90km de carretera. Por otra parte, la mayoría de las secuencias de cianobacterias de diferentes ecosistemas y latitudes descargadas de la página del NCBI, tuvieron altos valores de remuestreo por bootstrapping entre sí. Por estas razones, estas secuencias fueron de utilidad como outgroups para las amplificaciones de los dos genes de toxicidad evaluados. Para el gen mcyA, el dato más sorprendente se relaciona con la ubicación de una secuencia del embalse La Fe (LAFE4A), en el agrupamiento de secuencias del NCBI. No obstante y como ya se ha mencionado, la amplificación de esta región para la mayoría de estudios relacionados ha mostrado un alto grado de inespecificidad debido a que el gen que codifica para la subunidad A del cluster que codifica para el heptapéptido hepatotóxico microcistina-LR, está presente en diversas clases de cianobacterias productoras de microcistinas (Hisbergues et al., 2003). Es posible que la secuencia del embalse La Fe hallada dentro de la agrupación tomada como outgroup, difiera del resto de secuencias encontradas en los embalses antioqueños, por motivos relacionados con aspectos de la cuenca que hayan permitido que se mantengan cepas diferentes en cada sitio de muestreo dentro de cada embalse. Por otra parte, para el gen mcyE los análisis obtenidos de alinear secuencias de los embalses antioqueños con secuencias del NCBI de diferentes ecosistemas en una escala mundial, fueron más contundentes para separar las muestras secuenciadas de los embalses Riogrande II y La Fe. Aunque las cepas relacionadas en este estudio sean similares a cepas descriptas en otros países, cabe mencionar que los dos embalses evaluados en Antioquia son relativamente nuevos (aproximadamente 25 años para Riogrande II y 40 años para La Fe) y relativamente cercanos como ya se describió, lo cual podría explicar la alta similitud entre secuencias entre ambos ecosistemas.
Parte de los resultados del presente estudio confirman datos publicados por Correa en 2008 para el embalse Riogrande II. Este autor mostró que para el embalse Riogrande II se encontraban por taxonomía clásica tres géneros de cianobacterias con potencial producción de la toxina Microcistina-LR: Microcystis aeruginosa, M. wesenbergii, Anabaena sp. Y Radiocystis sp. en diversos eventos de florecimientos de cianobacterias (Correa, 2008). Pese a lo anterior, la presencia de Radiocystis debería confirmarse en evaluaciones futuras de toxicidad y diversidad cianobacterianas en el embalse Riogrande II, dado que su presencia se ha limitado principalmente a estudios llevados a cabo en el Brasil. Es importante enfatizar el hecho de encontrar dos cepas diferentes de Microcystis en muestras de agua evaluadas para el gen mcyE en el embalse La Fe, puesto que es un sitio que en los últimos años no ha tenido problemas de presencia de factores que desencadenen problemas de eutrofización (Palacio, com. pers.). De esta forma, en el presente trabajo pudo detectarse la presencia de genes de toxicidad de cianobacterias, sin evidencias visibles de eventos de florecimientos masivos de cianobacterias potencialmente tóxicas, para los embalses Riogrande II y La Fe. Sin embargo, aunque en este trabajo no se cuenta con datos de mediciones de parámetros físico-químicos de las estaciones de muestreo, es importante la detección inicial de secuencias de genes asociados con toxinas de cianobacterias en embalses de importancia para potabilización de agua.
Las evaluaciones microscópicas de cianobacterias son lábiles a errorres en la identificación de taxones, lo cual depende en gran medida de la experiencia del investigador. Por ello, la caracterización molecular de comunidades cinaobacterianas usando secuencias parciales de la subunidad 16S del ARN ribosomal, puede ser una herramienta complementaria para la identificación de microorganismos que prevalecen en un hábitat en particular (Dadheech, Krienitz, Kotut, Ballot & Casper, 2009). De acuerdo con esta idea, los resultados obtenidos por DGGE (Electroforesis en Geles con Gradiente Desnaturalizante) pueden proporcionar una separación visual de la mezcla de ADN en las muestras (fingerprinting, huellas genéticas) en función del polimorfismo en las secuencias, sin necesidad de usar enzimas de restricción (Romaní, Artigas, Camacho, Graça & Pascoal, 2009).
En el presente trabajo, aunque se amplificaron la mayoría de muestras para los dos fragmentos evaluados de la subunidad 16S rRNA, no fue posible reamplificar las eluciones de ADN de bandas obtenidas de geles de poliacrilamida. No obstante, se hallaron asociaciones de las muestras amplificadas de la region 359-781 de este gen en los embalses Riogrande II y La Fe. Dado que en varios estudios de cianobacterias de muestras ambientales se han hallado diferencias en las cepas encontradas con potencial toxicidad (Dos Anjos et al., 2006; Rantala et al., 2006; Magana-Arachchi, Wanigatunge & Jeyanandarajah, 2008), es factible que se encuentren diferencias entre las cepas encontradas en este estudio.
En otros estudios, se han realizado análisis similares a partir de la información obtenida por los patrones de bandeo de geles de poliacrilamida en gradientes desnaturalizantes de urea y formamida, para la evaluación de diversidad molecular de la región 359-781 del gen 16S del ARN ribosomal (rRNA) de cianobacterias en ecosistemas de aguas continentales. En 2009, Dadheech y colaboradores, evaluaron la diversidad de cepas de cianobacterias en sedimentos y agua de varios lagos africanos, mediante perfiles de DGGE. Para cada ecosistema evaluado, encontraron entre 5-11 bandas para cada perfil. El patrón de bandeo obtenido de DGGE fue similar entre réplicas de las muestras. Así mismo, el análisis de cluster del dendrograma mediante el algoritmo UPGMA mostró cuatro agrupaciones principales, de las cuales algunas agrupaciones eran compartidas por muestras de varios lagos (Dadheech et al., 2009). Por su parte, en 2008 Dorador y colaboradores evaluaron la diversidad de cianobacterias en humedales salinos en el norte de Chile. Estos investigadores encontraron patrones de DGGE de entre 5 y 11 bandas en muestras de agua, y entre 9-14 patrones de bandas de DGGE en muestras de sedimentos, lo cual pudo utilizarse como indicadores de riqueza de especies. Estos datos revelaron una habilidad mejorada para la detección de cianobacterias usando métodos moleculares, al contrastar esta técnica con técnicas de detección de la taxonomía tradicional (Dorador, Vila, Imhoff & Witzel, 2008).
Para el presente trabajo, es importante resaltar que fue posible la detección de un alto número de bandas para los ecosistemas evaluados, lo que lo destaca al compararse con los otros estudios mencionados (Riogrande II, 35 bandas diferentes; La Fe, 30 bandas diferentes). Con relación a la reamplificación de las bandas eluídas de DGGE, es claro que no en todos los casos se pueden obtener secuencias de las bandas escogidas para evaluaciones de diversidad, tal y como sucedió en este estudio. Sin embargo, en el caso de tener éxito en la reamplificación de muestras de bandas de DGGE, es preciso mencionar que se puede contar con información de secuencias idénticas de bandas diferentes dentro de la misma población analizada (Boutte, Grubisic, Balthasart & Wilmotte, 2006). Esto puede suceder debido a que amplicones dominantes podrían distribuirse en posiciones diferentes dentro del mismo patrón de bandeo, probablemente porque varios dominios tienen propiedades bioquímicas similares, siendo lábiles a efectos estocásticos que podrían causar que un dominio se desnaturalice antes que otros (Nikolausz, Sipos, Revesz, Szekely & Marialigeti, 2005; Boutte et al., 2006). Independientemente de la imposibilidad de re-amplificar secuencias de la subunidad 16S rRNA, en este trabajo fue posible determinar el grado de similitud entre los hallazgos obtenidos en cada embalse, de acuerdo con las agrupaciones obtenidas por el programa de análisis de bandas GelCompar II. Este análisis genético permitió determinar la utilidad de la técnica DGGE, para el hallazgo de perfiles genéticos únicos para cada ecosistema evaluado.
Aunque este trabajo es de carácter informativo de la presencia o ausencia de cepas de cianobacterias con potencial toxicidad, es recomendable incluir información de carácter físico-químico y ambiental de los sitios de muestreo para correlacionar dichas evaluaciones con los datos obtenidos de amplificación de regiones de toxicidad en diferentes épocas del año en los embalses Riogrande II y La Fe. Así mismo, puede recomendarse para futuros trabajos la detección del número de copias de genes y su expresión por qPCR (Quantitative real-time PCR, PCR cuantitativa en tiempo real) de los diferentes genes que hacen parte del cluster mcy, algunos de ellos que codifican para péptidos involucrados con la potencial toxicidad de cepas de cianobacterias productoras de microcistinas (Ngwa, Madramootoo & Jabaji, 2012).
Por último, se podrían llevar a cabo evaluaciones adicionales de diversidad de cepas y géneros de cianobacterias en los embalses Riogrande II y La Fe, mediante reamplificación de bandas obtenidas por la técnica DGGE, o mediante la secuenciación de regiones del gen 16S en una aproximación de metagenomas. De esta forma, los métodos moleculares independientes de cultivos cianobacterianos, relacionados con extracción de ADN directamente de muestras ambientales, se están haciendo más comunes como formas de estudiar microorganismos provenientes de ecosistemas complejos con problemas de eutrofización (Pope & Patel, 2008). Estos trabajos pueden ser de gran utilidad para estimar la posible relación entre la presencia de diversas cepas de cianobacterias con las condiciones ambientales, y para determinar posibles relaciones evolutivas entre cepas presentes en el mismo ecosistema. Por su condición de procariotes y por los eventos de eutrofización artificial que se pueden presentar en estos ecosistemas, puede ser interesante estudiar la capacidad de evolución de las cepas de cianobacterias en términos de diferencias en las tasas de mutación y/o sustitución nucleotídica de secuencia específicas, lo cual puede ser indicador a su vez de su capacidad de adaptación a estos ecosistemas artificiales. Estudios adicionales de la diversidad molecular de cianobacterias tóxicas podrían incluir el uso de estas herramientas de la biología molecular y de la bioinformática, siendo de utilidad en el apoyo de estudios limnológicos para la detección espacial o temporal de estos microorganismos, lo cual puede ser a su vez una herramienta nueva para la toma de decisiones en plantas de potabilización de agua.
Agradecimientos
A la DIME - Dirección de Investigaciones de la Universidad Nacional de Colombia Sede Medellín, por la financiación del presente estudio, derivado del proyecto “Desarrollo de bioensayos para la detección de toxinas de cianobacterias”, Código 20201007168. A Jaime Palacio y Néstor Aguirre del Grupo GAIA de la Universidad de Antioquia, por su colaboración en la obtención de las muestras de agua utilizadas en este estudio. Al Laboratorio de Biología Molecular y Celular de la Universidad Nacional de Colombia Sede Medellín, por el uso de sus instalaciones y equipos para el componente de biología molecular. A Mónica Higuita, Claudia Moreno y Sebastián Reynaldi, de la Universidad Nacional de Colombia Sede Medellín, por su colaboración en los análisis de datos moleculares, por la asesoría en la técnica DGGE, y por su ayuda inicial en la escritura del proyecto y la consecución de recursos para el desarrollo de este trabajo, respectivamente.
Referencias
Molecular techniques for cyanobacteria detection at Riogrande II and la Fe water reservoirs, Colombia
Julio César Hurtado-Alarcón1* & Jaime Polanía-Vorenberg2*
*Dirección para correspondencia:
Abstract
In lentic water bodies as reservoirs occur eutrophication processes, originated mainly from human activities (i.e. agriculture, animal exploitation). This influx of nutrients in aquatic ecosystems could promote blooms of potentially toxic cyanobacteria. The purpose of this work is to detect the presence of cyanobacteria strains in water samples, using molecular techniques to help in preventive management of reservoirs dedicated to water purification. We used two molecular techniques to detect genes implied with the synthesis of hepatotoxic microcystins from potentially toxic cyanobacteria strains, and to evaluate the molecular diversity of cyanobacteria in water samples from two high-mountain reservoirs used for purification of drinking water for the metropolitan area of Medellin, Colombia. Between 2010-2011 collections of 12 water samples were taken and DNA extraction together with PCR and DGGE analyses where carried out. We amplified 22 sequences between 250-300bp of the genes mcyA and mcyE, and these sequences were related with several strains and cyanobacteria genera accessions from NCBI-GenBank databases. Moreover, sequence amplifications of the 16S small ribosomal RNA subunit - 16S rRNA- between 400-800bp were also performed in order to use them for the DGGE technique. The amplification products of DGGE were set in polyacrilamide gel with posterior denaturing electrophoresis, and the scanned images of the gel bands were analysed with the software GelCompar II. For Riogrande II and La Fe reservoirs we found 35 and 30 different DGGE bands, respectively, as a measurement of molecular diversity in these artificial ecosystems. Here, we demonstrated the utility of two molecular techniques for the detection of genes associated with toxicity and molecular diversity of cyanobacteria in reservoirs destined for drinking water in urban centers. We recommend strongly following with periodically molecular biology studies in these ecosystems combined with limnological and ecological data, as new tools for management of plants of water purification and for the prevention of potentially public health concerns. Rev. Biol. Trop. 62 (1): 403-419. Epub 2014 March 01.
Key words: eutrophication, cyanobacteria, toxicity, microcystins, diversity, molecular markers.
Resumen
En embalses, la eutrofización es consecuencia de procesos naturales y de actividades humanas, lo cual puede facilitar la aparición de afloramientos de cianobacterias potencialmente tóxicas. En este trabajo, se utilizaron dos técnicas moleculares en la detección de genes presentes en cepas potencialmente tóxicas de cianobacterias y evaluaciones de diversidad molecular de cianobacterias en los embalses para potabilización de agua, Riogrande II y La Fe, Colombia. Entre 2010-2011, se tomaron 12 muestras de agua en ambos embalses y se realizaron extracciones de ADN para un análisis de marcadores moleculares mediante PCR y DGGE. Se amplificaron secuencias entre 250-300pb de los genes mcyA y mcyE, implicados en la toxicidad de cepas de cianobacterias. Asimismo, se amplificaron secuencias de la región 16S del ARN ribosomal (422pb), para la técnica DGGE. Se corrieron geles de poliacrilamida en gradientes de desnaturalización, se realizó agrupamiento genético (UPGMA), y se separaron por patrones de bandeo las muestras de cada embalse evaluado. Se demuestra la utilidad de las técnicas moleculares en estudios relacionados con la búsqueda de genes asociados con toxicidad y diversidad molecular de cianobacterias en muestras de agua provenientes de embalses de agua con fines de potabilización para centros urbanos.
Palabras clave: eutrofización, cianobacterias, toxicidad, microcistinas, diversidad, marcadores moleculares.
Medellín, en el Departamento de Antioquia, es la segunda ciudad más importante de Colombia. Su área metropolitana se encuentra distribuida a lo largo y ancho del Valle de Aburrá, y está conformada por los municipios de Medellín, Barbosa, Bello, Caldas, Copacabana, Envigado, Girardota, Itagüí, La Estrella y Sabaneta. Para el censo poblacional de 2005 en Colombia, se estimaron más de dos millones de habitantes en Medellín y cerca de un millón de habitantes en el resto de municipios del área metropolitana (DANE, 2005). Asimismo, algunos de estos municipios comparten la utilización de embalses para obtención de agua potable y energía eléctrica. En particular, parte del agua que abastece el acueducto de municipios del Valle de Aburrá como Medellín y Envigado, proviene principalmente de los embalses Riogrande II y La Fe, en los municipios de Don Matías y El Retiro, localizados al norte y al oriente de Medellín en el departamento de Antioquia, respectivamente.
Estudios anteriores han mostrado que embalses de características similares presentan procesos de eutrofización artificial, como consecuencia de actividades humanas (Ramírez, Gutiérrez & Vargas, 2005). Tales actividades, particularmente las agrícolas, han incrementado la concentración de nutrientes y estimulado la productividad primaria en los cuerpos de agua continentales (Havens, James, East & Smith, 2003). Este proceso de eutrofización artificial se manifiesta en floraciones visibles de plantas flotantes, agregaciones de macrófitas bénticas y, en algunos casos, en florecimientos de cianobacterias con potenciales efectos directos e indirectos sobre diferentes organismos, incluidos los seres humanos (Hitzfeld, Hoger & Dietrich, 2000). Además de los nutrientes, la temperatura, el pH o las altas intensidades lumínicas también pueden influir sobre los florecimientos de cianobacterias (Runnegar, Kong & Berndt, 1993; Dawson, 1998), con implicaciones directas e indirectas sobre los ecosistemas, puesto que estos microorganismos tienen la potencialidad de producir toxinas como estrategia de defensa. De esta forma, los peces, las aves y los seres humanos, entre otras formas de vida, pueden ser afectados por toxinas cianobacterianas -cianotoxinas- (Tencalla & Dietrich, 1997; Azevedo et al., 2002). Autores como Azevedo & Carmouze (1994), Codd, Steffensen, Burch, & Baker (1994), Rodger, Turnbull, Edwards & Codd (1994), entre otros, han registrado altas mortalidades de peces y animales domésticos expuestos a cianotoxinas. Por otra parte, en Caruarú, noreste del Brasil, fallecieron 76 pacientes de una clínica de hemodiálisis por contaminación del agua con la cianotoxina hepatotóxica microcistina-LR (Jochimsen et al., 1998; Domingos, Rubim, Molica, Azevedo & Carmichael, 1999; Carmichael et al., 2001).
Para la detección de cianobacterias, el uso de técnicas de microscopía basada en su morfología ha sido una herramienta de uso común. No obstante, cianobacterias potencialmente tóxicas no tienen características morfométricas únicas que ayuden en su identificación (Ouellete, Handy & Wilhelm, 2006). En algunos casos, las cepas tóxicas de géneros cianobacterianos se diferencian de las no tóxicas porque presentan genes que codifican para la síntesis de péptidos potencialmente tóxicos, los cuales pueden detectarse por métodos moleculares (Hotto, Satchwell & Boyer, 2007). Uno de ellos, la Reacción en Cadena de la Polimerasa (PCR, Polymerase Chain Reaction), permite detectar genes del cluster mcy, algunos de los cuales codifican para proteínas del péptido hepatotóxico microcistina sintetasa (Hisbergues, Christiansen, Rouhiainen, Sivonen & Börner, 2003). Dentro de este cluster, la detección de las regiones mcyA y mcyE ha sido de particular importancia, puesto que son importantes en la identificación de géneros de cianobacterias potencialmente tóxicas (Rantala et al., 2006; Hotto et al., 2007). Así mismo, se han llevado a cabo diferentes estudios para analizar la ocurrencia natural de cianobacterias en ecosistemas acuáticos mediante la amplificación de secuencias de la subunidad 16S del ADN ribosomal, lo cual es útil en el monitoreo de la distribución de especies (Matsunaga, Takeyama & Nakayama, 2001), en estudios de biodiversidad (Ouellete et al., 2006), y en estudios filogenéticos y evolutivos (Genuário, Silva-Stenico, Welker, Moraesc & Fiore, 2010).
Dado que la Organización Mundial de la Salud (WHO) permite un máximo de 1µg/L de toxinas como microcistina-LR en agua para potabilizar (WHO, 2004), es de vital importancia detectar la presencia de cianobacterias con potencial toxicidad en cuerpos de agua lénticos dedicados a la obtención de agua potable para millones de personas. Considerando que Correa (2008) identificó dos géneros diferentes de cianobacterias productoras de Microcistina-LR en el embalse Riogrande II, Colombia, mediante técnicas de microscopía tradicional, en este artículo se describe la detección molecular de cepas de cianobacterias con potencial toxicidad en los embalses Riogrande II y La Fe, Colombia, entre los años 2010-11, mediante la amplificación de los genes mcyA y mcyE. Así mismo, se amplificó un fragmento de la subunidad 16S del ARN ribosomal (rRNA) cianobacteriano, y se realizó Electroforesis en Gradientes de Desnaturalización en Geles de poliacrilamida (DGGE), en un intento para detectar diferencias o similitudes en los patrones de bandeo obtenidos por esta técnica molecular con relación a los ecosistemas evaluados, como un estudio complementario para el monitoreo de la presencia de géneros y cepas de cianobacterias potencialmente tóxicas en estos ecosistemas artificiales en Colombia.
Materiales y métodos
Área de estudio: Se obtuvieron muestras de agua de los embalses Riogrande II y La Fe, Colombia, en seis salidas de campo realizadas entre 2010-2011. El embalse Riogrande II se encuentra en la zona central del Departamento de Antioquia en jurisdicción de los municipios de San Pedro, Entrerríos, Belmira, Don Matías y Santa Rosa de Osos (Otaya, Vásquez & Bustamante, 2008). Se construyó entre 1986-1988, y está formado por la confluencia de los ríos Chico y Grande, y la quebrada Las Ánimas; se encuentra a una altura promedio de 2 270msnm, y presenta un área superficial de 12 000km2, un volumen de 220Mm3, una profundidad máxima de 42m, una profundidad promedio de 37.9m y una longitud máxima de 10km (Montoya & Ramírez, 2007). El agua que se capta es aprovechada por las Empresas Públicas de Medellín –EPM- para el abastecimiento del acueducto de la ciudad y algunos municipios del Valle de Aburrá, y también es utilizada para la generación de energía en las centrales hidroeléctricas Tasajera y Niquía, generando alrededor de 324MW (CORANTIOQUIA, 1997). Por su parte, el embalse La Fe fue llenado en 1 973, y cubre un área de 173km2. Se sitúa en el municipio de El Retiro, con alturas entre los 2 175-3 000msnm. Este embalse es utilizado principalmente para acueducto, pero una parte está destinada a recreación. De su volumen total (aprox. 15Mm3), 12Mm3 regulan el río Pantanillo, y 8m3/s están destinados al acueducto metropolitano a través de la planta de tratamiento de La Ayurá, municipio de Envigado, Antioquia (Abuchaibe, Agudelo & Sañudo 1988).
Recolecta de muestras de agua: Mediante redes de arrastre de fitoplancton de 35µm de diámetro de poro, se tomaron entre 100-300mL de agua en la zona fótica en los embalses Riogrande II y La Fe. Dentro de cada embalse se tomaron muestras de agua en diferentes estaciones de muestreo (cuatro estaciones en Riogrande II: 1, 2, 3, 8; y seis estaciones en La Fe: 1, 2, 3, 4, 6, 7), las cuales aparecen con sus nombres sintetizados y coordenadas geográficas en el cuadro 1.
Las muestras de agua se pasaron por filtros de celulosa de 0.5µm y se refrigeraron a 4°C hasta la noche o hasta el día siguiente en el laboratorio del grupo GAIA –Grupo de Investigación, Gestión y Modelación Ambiental–, en la Sede de Investigación Universitaria de la Universidad de Antioquia, SIU. Posteriormente, se transportaron hacia el Laboratorio de Biología Molecular y Celular de la Universidad Nacional de Colombia, Sede Medellín. Allí, las muestras se homogenizaron y de cada una de ellas se tomaron volúmenes de 50mL por duplicado, y se centrifugaron en volúmenes decrecientes de 50, 15, y 2mL entre 200-3 000g durante 10-15min. De esta forma, el precipitado depositado en los viales de menor volumen (2mL), contenía la materia prima para llevar a cabo las extracciones de ADN, las cuales se realizaron usando el estuche QIAamp DNA Mini Kit (Qiagen), de acuerdo con las recomendaciones del fabricante. La concentración y pureza del ADN extraído se evaluó mediante un espectrofotómetro Nanodrop (Thermoscientific, USA). La calidad de las extracciones de ADN se verificó en geles de agarosa al 1.0% teñidos con bromuro de etidio (EtBr) (10mg/mL), y fueron visualizadas en transiluminador de luz ultravioleta.
Amplificación por PCR (Polymerase Chain Reaction, Reacción en Cadena de la Polimerasa) de regiones de interés: Inicialmente, se evaluaron las secuencias cebadoras cpcB-cpcA (PCβ-F:
5›-GGCTGCTTGTTTACGCGACA-3›, y PCα-R: 5›-CCAGTACCACCAGCAACTAA- 3› (Neilan, Jacobs & Goodman, 1995), para amplificar el gen que codifica para ficocianina-c específica de cianobacterias (Cuadro 2), como control positivo de la presencia de cianobacterias en las muestras de ADN obtenidas de los embalses. Para la amplificación posterior de fragmentos del cluster de genes mcy y de la región 16S del ARN ribosomal, solo se utilizaron muestras de ADN positivas para ficocianina-c.
Evaluación de genes relacionados con toxicidad de cianobacterias: Para identificar regiones relacionadas con toxicidad de cianobacterias, se evaluaron los cebadores mcyA-Cd (1F5’- AAAATTAAAAGCCGTATCAAA-3’;1R-5’AAAAGTGTTTTATTAGCGGCTCAT) (Hisbergues et al., 2003), y para la posible identificación de los géneros Microcystis y Anabaena se evaluaron cebadores específicos (mcyE-F2: 5’-GAAATTTGTGTAGAAGGTGC-3’ para ambos géneros; mcyE-R8 5’- CAATGGGAGCATAACGAG-3’, Microcystis; mcyE-12R 5’-CAATCTCGGTATAGCGGC-3’, Anabaena) (Vaitomaa et al., 2003; Rantala et al., 2006; Hotto et al., 2007) (Cuadro 2).
El volumen final de cada mezcla de PCR fue de 25mL. Las amplificaciones por PCR de las regiones de interés se realizaron de acuerdo con los protocolos publicados por Hisbergues et al. (2003) (gen mcyA), Vaitomaa et al. (2003) (gen mcyE) y Dos Anjos et al. (2006) (operon ficocianina cianobacteriana, cpcB-cpcA), con respecto a concentraciones de reactivos y perfiles térmicos. Las reacciones de PCR se llevaron a cabo en un termociclador T3 000 (Biometra, Alemania) en microtubos de 200μL. Los productos de amplificación se visualizaron mediante electroforesis en geles de agarosa (1.5%) teñidos con EtBr (10mg/mL). Las condiciones de corrido fueron: 60-80min en solución 1X TBE (Tris-ácido bórico-EDTA, pH8.0), a un voltaje constante (80V) en cámaras de electroforesis horizontales (Bio-Rad, E.U.A.). Los geles de agarosa se observaron sobre una pantalla de luz ultravioleta y se fotografiaron con un equipo para captura de imágenes (Biometra, Alemania) para confirmar la efectividad de la amplificación. Para determinar los tamaños de los alelos en cada corrida se insertó un marcador molecular de secuencia conocida en un rango de 100-1 000 pares de bases (GeneRuler Low Range DNA Ladder, Fermentas, Canada) (Hisbergues et al., 2003, Dos Anjos et al., 2006).
Purificación de las muestras positivas y secuenciación: Las muestras que resultaron positivas para genes de toxicidad (mcyA, mcyE) se purificaron utilizando el estuche comercial QIAquick PCR Purification Kit (QIAGEN, Alemania), siguiendo las recomendaciones del fabricante, las amplificaciones se verificaron en geles de agarosa (1.5%) y finalmente se enviaron a secuenciar (Macrogen, Corea). Se realizaron comparaciones individuales de las secuencias obtenidas de los genes mcyA y mcyE en las bases de datos del NCBI (National Center for Biotechnology Information, Centro Nacional de Información en Biotecnología, E.U.A.), usando la herramienta para alineamientos locales BLAST (Basic Local Alignment Search Tool).
Análisis de secuencias para los genes mcyA y mcyE: Se utilizó el editor de alineamientos gratuito BioEdit (Biological Sequence Alignment Editor) versión 7.0.9 (Hall, T.A. 1999), el cual se puede descargar gratuitamente en el sitio URL http://www.mbio.ncsu.edu/bioedit/bioedit.html. Este programa permite obtener secuencias consenso a partir de los resultados de las secuenciaciones de las cadenas Forward y Reverse de cada una de las muestras positivas para los genes mcyA y mcyE. Con las secuencias consenso, se realizó un alineamiento múltiple con el algoritmo Clustal W incluido en el programa de uso libre BioEdit. Adicionalmente, las secuencias consenso obtenidas fueron alineadas y analizadas junto con secuencias publicadas en las bases de datos del NCBI de algunos géneros y cepas de cianobacterias (Microcystis sp., Nostoc sp., Anabaena sp., Radiocystis sp.) que portan los genes mcyA y mcyE, relacionados con la síntesis de péptidos con potencial hepatotoxicidad (Cuadro 3). La escogencia del modelo de sustitución nucleotídica y los análisis de agrupamiento de las secuencias se llevaron a cabo con el programa de uso libre MEGA5, descargado del sitio URL http://www.megasoftware.net/ (Tamura et al., 2011).
Amplificación de secuencias de la subunidad 16S del ARN ribosomal (rRNA) para la evaluación de diversidad de cepas de cianobacterias mediante la técnica molecular DGGE (Denaturing gradient gel electrophoresis - Electroforesis en gel con gradiente de desnaturalización): Se evaluaron las secuencias cebadoras descritas por Nübel, Garcia-Pichel & Muyzer (1997), para amplificar regiones del gen que codifica para la subunidad 16S rRNA de cianobacterias (Cuadro 3), de aproximadamente 422 pares de bases (bp) (cebadores 359F-781R), a partir del ADN extraído de las muestras de agua de los embalses Riogrande II y La Fe (Cuadro 4). Se utilizó el cebador Forward 359F modificado en el extremo 3´, mediante la adición de 40 desoxirribonucleótidos fosfato guanina-citosina (GC), para modificar el comportamiento desnaturalizante del fragmento de interés y mejorar la resolución de los fragmentos de tamaño similar pero con contenido diferente de GC (Muyzer, De Waal & Uitterlinden, 1993). La mezcla de PCR se preparó de acuerdo con lo publicado por Nübel et al. (1997). Se tomaron 25rm de cada primer, 25hmol dNTPs, 10mL buffer PCR, (100mM Tris-HCl [pH 9.0], 1.5mM MgCl, 500mM KCl), 0.5U Taq polimerasa (FERMENTAS UAB, Canadá) y entre 30-100ηg de ADN molde obtenido de las extracciones, en un volumen final de 30-50mL por reacción. El perfil térmico de la PCR incluyó un ciclo inicial de desnaturalización del ADN a 95°C, seguido de 35 ciclos de un minuto cada uno a 94°C, 60°C y 72°C. Por último, se incluyó un ciclo de extensión final a 72°C por 10 min.
Evaluación de la técnica DGGE: Las muestras que amplificaron con los cebadores 359FGC-781R se precipitaron con etanol de grado reactivo al 96%, se llevaron a una concentración entre 2 000-3 000ηg y se diluyeron en un volumen final de 10μL con agua ultrapura estéril. Con estos productos de PCR se llevaron a cabo electroforesis en geles de poliacrilamida, en el sistema D-Code Mutation Detection System (Bio-Rad, E.U.A.). Para ello, se prepararon geles de poliacrilamida al 8% y 1mm de grosor, en dos gradientes desnaturalizantes de urea y formamida (30-60%), de acuerdo con los protocolos descritos por Nübel et al. (1997). Los geles de poliacrilamida se sometieron a electroforesis en solución tampón TAE (40mM Tris-HCl [pH 8.3], 20mM ácido acético, 1mM EDTA), durante 3.5 horas a 200V, a una temperatura constante de 60°C. Los geles de poliacrilamida se tiñeron con nitrato de plata (AgNO3) y se sumergieron en solución reveladora (NaOH, formamida y agua destilada) para visualizar bandas de ADN, (Sanguinetti, Dias Neto & Simpson, 1994). Posteriormente, los geles se fijaron con solución Etanol-ácido acético 95:5, se escanearon y se digitalizaron para el análisis de imágenes usando el programa GelCompar II 6.0 (Applied Maths, Bélgica). Este programa permite analizar agrupaciones (clusters) evaluando algoritmos para determinación de similitud de bandas de amplificaciones y polimorfismos en evaluaciones de diversidad molecular (Rademaker & de Bruijin, 2004). Dicha metodología es útil para evaluar similitudes o diferencias en los patrones de bandeo de las amplificaciones de la región 16S del ARN ribosomal cianobacteriano, obtenidos en ecosistemas lénticos artificiales e independientes como lo son los embalses Riogrande II y La Fe.
Por medio del programa GelCompar II 6.0, se calculó la matriz de distancia basada en las bandas de DGGE usando el coeficiente de similitud de Dice (Nei & Li, 1979), el cual es comúnmente usado para comparar especies entre ecosistemas diferentes. Dos perfiles idénticos dan como resultado un valor igual al 100%, mientras que dos perfiles completamente diferentes dan un valor igual al 0%. Se realizó un análisis de agrupamiento mediante el método de promedio aritmético no ponderado (UPGMA -Unweighted Pair Group Method using Arithmetic averages), y los resultados fueron presentados finalmente como un dendrograma.
Resultados
Extracción y cuantificación de ADN a partir de las muestras de agua tomadas en los embalses Riogrande II y La Fe: El trabajo de extracción de ADN usando el estuche comercial permitió obtener ADN de buena calidad en las muestras de agua obtenidas de los embalses Riogrande II y La Fe, Colombia. Así mismo, la cuantificación del ADN usando espectrofotómetro arrojó lecturas entre 7.3- 188.3ηg/mL. A partir de las muestras cuantificadas, se usaron concentraciones de 20ηg de ADN por cada muestra para las amplificaciones por PCR de las regiones de interés.
Amplificaciones de genes del cluster mcy, relacionados con toxicidad de cianobacterias:
Amplificación de los genes mcyA y mcyE: Para los dos embalses evaluados, se obtuvieron amplificaciones en las reacciones de PCR de los genes mcyA y mcyE en las muestras de agua obtenidas. Las amplificaciones de parte del dominio de condensación (Cd) del gen mcyA, tuvieron un tamaño aproximado al publicado por Hisbergues et al. (2003) de 300bp, utilizando el par de cebadores mcyA-Cd (Hisbergues et al., 2003). Para el gen mcyE, solo se obtuvieron amplificaciones por PCR con los cebadores específicos para Microcystis sp. (mic-mcyE R8). El tamaño obtenido para esta región fue diferente a lo encontrado por otros autores (aprox. 300bp) (Vaitomaa et al., 2003). Para Anabaena sp. no se obtuvieron amplificaciones positivas al usar los cebadores Ana-mcyE 12R.
Comparación de las secuencias consenso obtenidas por secuenciación de las amplificaciones de los genes mcyA y mcyE, usando la herramienta de alineamiento local BLAST: Para el gen mcyA que codifica para microcistina sintetasa, se hallaron correspondencias en el GenBank con varios clones para los embalses evaluados. Para el embalse Riogrande II, las consensos analizados correspondieron con un clon de cianobacterias sin identificar (Uncultured cyanobacterium clone 03LO-OS4, Nº accesión GenBank EF424345.1, cobertura 99%) y con un clon no cultivado de Microcystis (Uncultured Microcystis sp. clone MB20, Nº accesión GenBank JN108773.1, cobertura 100%). Para el embalse La Fe, de cuatro consensos dos se relacionaron con clones no cultivados de cianobacterias sin identificar, y dos consensos se relacionaron con uno de los clones hallados en el embalse Riogrande II (Uncultured cyanobacterium clone 03LO-OS4, Nº accesión GenBank EF424345.1, cobertura 99%; Uncultured cyanobacterium clone 03LOOR9, Nº accesión GenBank EF424337.1, cobertura 100%).
Con relación al gen mcyE que codifica para microcistina sintetasa para ambos embalses, se hallaron correspondencias similares. Para los embalses Riogrande II y La Fe, todos los 13 consensos diferentes se relacionaron con las mismas cepas de Microcystis (Microcystis sp. CYN06, Nº accesión GenBank FJ393327.1, cobertura 99-100%), y en una medida similar con una cepa de Microcystis aeruginosa (CPCC 299, Nº accesión GenBank JN936964.1, cobertura 99-100%, identidad >98%). No obstante para el embalse La Fe, de dos secuencias consenso, una de ellas tuvo altos valores de identidad (98%) pero con un valor bajo de cobertura de 44%.
Análisis de agrupamiento para las secuencias consenso obtenidas: El programa BioEdit, permitió encontrar siete secuencias consenso para el gen mcyA y 13 secuencias consenso para el gen mcyE. El programa MEGA 5 permitió obtener matrices de agrupamiento y de distancias genéticas para estos genes. Los análisis de las secuencias consenso de los genes mcyA y mcyE se basaron en el método de Máxima Verosimilitud (Maximum-likelihood, ML) usando el modelo de sustitución nucleotídica Kimura 2P (K2), y se obtuvieron dendrogramas de agrupamiento (Fig. 1 y Fig. 2), con topologías soportadas por análisis de Bootstrapping (1 000 réplicas). Las secuencias que se obtuvieron de las bases de datos del NCBI de algunos géneros y cepas de cianobacterias (Microcystis sp., Nostoc sp., Anabaena sp., Radiocystis sp.), se mostraron como outgroups para las secuencias obtenidas en el presente estudio, debido a altos valores de remuestreo mediante bootstrapping hallados entre los clados obtenidos (99%), exceptuando una secuencia obtenida del embalse La Fe para el gen mcyA.
Para el gen mcyA se encontraron dos agrupaciones de secuencias, la primera que comprende ocho de las nueve secuencias de los embalses Riogrande II y La Fe, y la segunda agrupación la cual incluye seis secuencias del NCBI más una sola secuencia del embalse La Fe. Para este gen, los valores de remuestreo bootstrapping entre la mayoría de las secuencias de los embalses antioqueños fueron muy bajos (1-11% de las replicaciones del remuestreo). Los valores de bootstrapping para la segunda agrupación de secuencias (una secuencia del embalse La Fe y las secuencias obtenidas del NCBI), fueron mayores a los valores obtenidos para la primera agrupación (41-99% de las replicaciones del remuestreo), siendo el valor más bajo (41%) el obtenido al comparar cinco de las seis muestras del NCBI con la secuencia del embalse La Fe.
Para el gen mcyE, se hallaron varias agrupaciones dentro de las muestras de los embalses evaluados, separando claramente dichas muestras de las secuencias obtenidas del NCBI para el gen mcyE. Los valores obtenidos de remuestreo mediante bootstrapping fueron muy bajos para 11 de las 13 secuencias de los embalses Riogrande II y La Fe (0-4% de las replicaciones del remuestreo), y tres de las 13 muestras de los embalses antioqueños tuvieron valores altos de bootstrapping (58-99% de las replicaciones del remuestreo).
Evaluación de diversidad de cepas cianobacterianas en los embalses Riogrande II y La Fe:
Evaluación de la amplificación de secuencias de la subunidad 16S del ARN ribosomal (rRNA) y su evaluación en geles de poliacrilamida para la técnica DGGE: Para los embalses Riogrande II y La Fe, se obtuvieron amplificaciones de la región 359FGC-781F(a,b) del gen 16S del ARN ribosomal (rRNA), las cuales se corrieron en geles de poliacrilamida para la técnica DGGE. Las bandas obtenidas se pudieron visualizar mediante los protocolos de tinción con nitrato de plata (Fig. 3).
Análisis de los agrupamientos por patrones de bandeo de las reacciones obtenidas por amplificación y DGGE de la región 16S rRNA: Los geles de poliacrilamida digitalizados y analizados con el programa GelCompar II permitieron realizar agrupamientos de bandas comunes y no comunes en los embalses Riogrande II y La Fe. Se detectaron 35 bandas diferentes en las muestras del embalse Riogrande II, y 30 bandas diferentes para el embalse La Fe, con las cuales se realizó el dendrograma descripto de agrupaciones genéticas (Fig. 4). Se agruparon la mayoría de sitios de muestreo dentro de cada embalse a partir de los geles de poliacrilamida obtenidos. El dendrograma obtenido separa dos agrupaciones genéticas, diferenciando con claridad las muestras del embalse Riogrande II de las del embalse La Fe, exceptuando las muestras del sitio 7A. (Fig. 5).
Discusión
En este proyecto, la mayoría de cebadores evaluados para amplificar regiones de ADN de interés, permitieron obtener resultados positivos en las muestras de agua de los embalses Riogrande II y La Fe, para genes que codifican para ficocianina de cianobacterias y algunos péptidos asociados con toxicidad (mcyA, mcyE). Así mismo, se obtuvo éxito en la amplificación de regiones del gen que codifica para la subunidad 16S del ARN ribosomal (rRNA). Estos resultados son importantes debido a que son indicadores de la presencia de cepas de cianobacterias con potencial toxicidad en los embalses Riogrande II y La Fe, destinados para la potabilización de agua en centros urbanos de la ciudad de Medellín, Colombia.
La amplificación de fragmentos de la región cpc-IGS específica de cianobacterias en muestras de agua de embalses, puede dar una idea inicial de la presencia de cianobacterias. No obstante, su identificación en muestras de agua no permite determinar cuáles cepas o géneros de cianobacterias están presentes en las muestras evaluadas. Así mismo, con la amplificación de esta región no se logra diferenciar cepas tóxicas (Bittencourt-Oliveira et al., 2001). Sin embargo, los cebadores evaluados para los genes relacionados con toxicidad mcyA y mcyE, ayudaron a detectar regiones relacionadas con cepas y géneros cianobacterianos como Microcystis con potencial toxicidad.
Los resultados de las amplificaciones inespecíficas de regiones del gen mcyA para muestras de agua de los embalses Riogrande II y La Fe eran de esperarse. La amplificación del gen mcyA solo apunta a describir la presencia de cepas potencialmente productoras de microcistinas sin discriminación de cepas o especies, lo cual ha sido relacionado en estudios anteriores (Hisbergues et al., 2003; Hotto et al., 2007). Esto es diferente de los resultados que pueden arrojar amplificaciones del gen mcyE, el cual usa cebadores específicos de cepas y especies (Vaitomaa et al., 2003; Rantala et al., 2006). De esta forma, para la región mcyE, también se encontró el género cianobacteriano Microcystis sp. a partir de la información de secuenciación contrastada en las bases de datos geneticos del NCBI, con representantes de dos cepas diferentes de M. aeruginosa para los dos embalses evaluados.
De acuerdo con los dendrogramas obtenidos para ambos genes, se indica un alto grado de similitud entre las mismas secuencias y poca divergencia genética, para las secuencias de los dos genes evaluados en los embalses antioqueños (mcyA, mcyE). Estos datos indican a su vez alta similitud genética entre la mayoría de secuencias analizadas en este trabajo, en dos embalses antioqueños separados geográficamente por menos de 50km linealmente y por 90km de carretera. Por otra parte, la mayoría de las secuencias de cianobacterias de diferentes ecosistemas y latitudes descargadas de la página del NCBI, tuvieron altos valores de remuestreo por bootstrapping entre sí. Por estas razones, estas secuencias fueron de utilidad como outgroups para las amplificaciones de los dos genes de toxicidad evaluados. Para el gen mcyA, el dato más sorprendente se relaciona con la ubicación de una secuencia del embalse La Fe (LAFE4A), en el agrupamiento de secuencias del NCBI. No obstante y como ya se ha mencionado, la amplificación de esta región para la mayoría de estudios relacionados ha mostrado un alto grado de inespecificidad debido a que el gen que codifica para la subunidad A del cluster que codifica para el heptapéptido hepatotóxico microcistina-LR, está presente en diversas clases de cianobacterias productoras de microcistinas (Hisbergues et al., 2003). Es posible que la secuencia del embalse La Fe hallada dentro de la agrupación tomada como outgroup, difiera del resto de secuencias encontradas en los embalses antioqueños, por motivos relacionados con aspectos de la cuenca que hayan permitido que se mantengan cepas diferentes en cada sitio de muestreo dentro de cada embalse. Por otra parte, para el gen mcyE los análisis obtenidos de alinear secuencias de los embalses antioqueños con secuencias del NCBI de diferentes ecosistemas en una escala mundial, fueron más contundentes para separar las muestras secuenciadas de los embalses Riogrande II y La Fe. Aunque las cepas relacionadas en este estudio sean similares a cepas descriptas en otros países, cabe mencionar que los dos embalses evaluados en Antioquia son relativamente nuevos (aproximadamente 25 años para Riogrande II y 40 años para La Fe) y relativamente cercanos como ya se describió, lo cual podría explicar la alta similitud entre secuencias entre ambos ecosistemas.
Parte de los resultados del presente estudio confirman datos publicados por Correa en 2008 para el embalse Riogrande II. Este autor mostró que para el embalse Riogrande II se encontraban por taxonomía clásica tres géneros de cianobacterias con potencial producción de la toxina Microcistina-LR: Microcystis aeruginosa, M. wesenbergii, Anabaena sp. Y Radiocystis sp. en diversos eventos de florecimientos de cianobacterias (Correa, 2008). Pese a lo anterior, la presencia de Radiocystis debería confirmarse en evaluaciones futuras de toxicidad y diversidad cianobacterianas en el embalse Riogrande II, dado que su presencia se ha limitado principalmente a estudios llevados a cabo en el Brasil. Es importante enfatizar el hecho de encontrar dos cepas diferentes de Microcystis en muestras de agua evaluadas para el gen mcyE en el embalse La Fe, puesto que es un sitio que en los últimos años no ha tenido problemas de presencia de factores que desencadenen problemas de eutrofización (Palacio, com. pers.). De esta forma, en el presente trabajo pudo detectarse la presencia de genes de toxicidad de cianobacterias, sin evidencias visibles de eventos de florecimientos masivos de cianobacterias potencialmente tóxicas, para los embalses Riogrande II y La Fe. Sin embargo, aunque en este trabajo no se cuenta con datos de mediciones de parámetros físico-químicos de las estaciones de muestreo, es importante la detección inicial de secuencias de genes asociados con toxinas de cianobacterias en embalses de importancia para potabilización de agua.
Las evaluaciones microscópicas de cianobacterias son lábiles a errorres en la identificación de taxones, lo cual depende en gran medida de la experiencia del investigador. Por ello, la caracterización molecular de comunidades cinaobacterianas usando secuencias parciales de la subunidad 16S del ARN ribosomal, puede ser una herramienta complementaria para la identificación de microorganismos que prevalecen en un hábitat en particular (Dadheech, Krienitz, Kotut, Ballot & Casper, 2009). De acuerdo con esta idea, los resultados obtenidos por DGGE (Electroforesis en Geles con Gradiente Desnaturalizante) pueden proporcionar una separación visual de la mezcla de ADN en las muestras (fingerprinting, huellas genéticas) en función del polimorfismo en las secuencias, sin necesidad de usar enzimas de restricción (Romaní, Artigas, Camacho, Graça & Pascoal, 2009).
En el presente trabajo, aunque se amplificaron la mayoría de muestras para los dos fragmentos evaluados de la subunidad 16S rRNA, no fue posible reamplificar las eluciones de ADN de bandas obtenidas de geles de poliacrilamida. No obstante, se hallaron asociaciones de las muestras amplificadas de la region 359-781 de este gen en los embalses Riogrande II y La Fe. Dado que en varios estudios de cianobacterias de muestras ambientales se han hallado diferencias en las cepas encontradas con potencial toxicidad (Dos Anjos et al., 2006; Rantala et al., 2006; Magana-Arachchi, Wanigatunge & Jeyanandarajah, 2008), es factible que se encuentren diferencias entre las cepas encontradas en este estudio.
En otros estudios, se han realizado análisis similares a partir de la información obtenida por los patrones de bandeo de geles de poliacrilamida en gradientes desnaturalizantes de urea y formamida, para la evaluación de diversidad molecular de la región 359-781 del gen 16S del ARN ribosomal (rRNA) de cianobacterias en ecosistemas de aguas continentales. En 2009, Dadheech y colaboradores, evaluaron la diversidad de cepas de cianobacterias en sedimentos y agua de varios lagos africanos, mediante perfiles de DGGE. Para cada ecosistema evaluado, encontraron entre 5-11 bandas para cada perfil. El patrón de bandeo obtenido de DGGE fue similar entre réplicas de las muestras. Así mismo, el análisis de cluster del dendrograma mediante el algoritmo UPGMA mostró cuatro agrupaciones principales, de las cuales algunas agrupaciones eran compartidas por muestras de varios lagos (Dadheech et al., 2009). Por su parte, en 2008 Dorador y colaboradores evaluaron la diversidad de cianobacterias en humedales salinos en el norte de Chile. Estos investigadores encontraron patrones de DGGE de entre 5 y 11 bandas en muestras de agua, y entre 9-14 patrones de bandas de DGGE en muestras de sedimentos, lo cual pudo utilizarse como indicadores de riqueza de especies. Estos datos revelaron una habilidad mejorada para la detección de cianobacterias usando métodos moleculares, al contrastar esta técnica con técnicas de detección de la taxonomía tradicional (Dorador, Vila, Imhoff & Witzel, 2008).
Para el presente trabajo, es importante resaltar que fue posible la detección de un alto número de bandas para los ecosistemas evaluados, lo que lo destaca al compararse con los otros estudios mencionados (Riogrande II, 35 bandas diferentes; La Fe, 30 bandas diferentes). Con relación a la reamplificación de las bandas eluídas de DGGE, es claro que no en todos los casos se pueden obtener secuencias de las bandas escogidas para evaluaciones de diversidad, tal y como sucedió en este estudio. Sin embargo, en el caso de tener éxito en la reamplificación de muestras de bandas de DGGE, es preciso mencionar que se puede contar con información de secuencias idénticas de bandas diferentes dentro de la misma población analizada (Boutte, Grubisic, Balthasart & Wilmotte, 2006). Esto puede suceder debido a que amplicones dominantes podrían distribuirse en posiciones diferentes dentro del mismo patrón de bandeo, probablemente porque varios dominios tienen propiedades bioquímicas similares, siendo lábiles a efectos estocásticos que podrían causar que un dominio se desnaturalice antes que otros (Nikolausz, Sipos, Revesz, Szekely & Marialigeti, 2005; Boutte et al., 2006). Independientemente de la imposibilidad de re-amplificar secuencias de la subunidad 16S rRNA, en este trabajo fue posible determinar el grado de similitud entre los hallazgos obtenidos en cada embalse, de acuerdo con las agrupaciones obtenidas por el programa de análisis de bandas GelCompar II. Este análisis genético permitió determinar la utilidad de la técnica DGGE, para el hallazgo de perfiles genéticos únicos para cada ecosistema evaluado.
Aunque este trabajo es de carácter informativo de la presencia o ausencia de cepas de cianobacterias con potencial toxicidad, es recomendable incluir información de carácter físico-químico y ambiental de los sitios de muestreo para correlacionar dichas evaluaciones con los datos obtenidos de amplificación de regiones de toxicidad en diferentes épocas del año en los embalses Riogrande II y La Fe. Así mismo, puede recomendarse para futuros trabajos la detección del número de copias de genes y su expresión por qPCR (Quantitative real-time PCR, PCR cuantitativa en tiempo real) de los diferentes genes que hacen parte del cluster mcy, algunos de ellos que codifican para péptidos involucrados con la potencial toxicidad de cepas de cianobacterias productoras de microcistinas (Ngwa, Madramootoo & Jabaji, 2012).
Por último, se podrían llevar a cabo evaluaciones adicionales de diversidad de cepas y géneros de cianobacterias en los embalses Riogrande II y La Fe, mediante reamplificación de bandas obtenidas por la técnica DGGE, o mediante la secuenciación de regiones del gen 16S en una aproximación de metagenomas. De esta forma, los métodos moleculares independientes de cultivos cianobacterianos, relacionados con extracción de ADN directamente de muestras ambientales, se están haciendo más comunes como formas de estudiar microorganismos provenientes de ecosistemas complejos con problemas de eutrofización (Pope & Patel, 2008). Estos trabajos pueden ser de gran utilidad para estimar la posible relación entre la presencia de diversas cepas de cianobacterias con las condiciones ambientales, y para determinar posibles relaciones evolutivas entre cepas presentes en el mismo ecosistema. Por su condición de procariotes y por los eventos de eutrofización artificial que se pueden presentar en estos ecosistemas, puede ser interesante estudiar la capacidad de evolución de las cepas de cianobacterias en términos de diferencias en las tasas de mutación y/o sustitución nucleotídica de secuencia específicas, lo cual puede ser indicador a su vez de su capacidad de adaptación a estos ecosistemas artificiales. Estudios adicionales de la diversidad molecular de cianobacterias tóxicas podrían incluir el uso de estas herramientas de la biología molecular y de la bioinformática, siendo de utilidad en el apoyo de estudios limnológicos para la detección espacial o temporal de estos microorganismos, lo cual puede ser a su vez una herramienta nueva para la toma de decisiones en plantas de potabilización de agua.
Agradecimientos
A la DIME - Dirección de Investigaciones de la Universidad Nacional de Colombia Sede Medellín, por la financiación del presente estudio, derivado del proyecto “Desarrollo de bioensayos para la detección de toxinas de cianobacterias”, Código 20201007168. A Jaime Palacio y Néstor Aguirre del Grupo GAIA de la Universidad de Antioquia, por su colaboración en la obtención de las muestras de agua utilizadas en este estudio. Al Laboratorio de Biología Molecular y Celular de la Universidad Nacional de Colombia Sede Medellín, por el uso de sus instalaciones y equipos para el componente de biología molecular. A Mónica Higuita, Claudia Moreno y Sebastián Reynaldi, de la Universidad Nacional de Colombia Sede Medellín, por su colaboración en los análisis de datos moleculares, por la asesoría en la técnica DGGE, y por su ayuda inicial en la escritura del proyecto y la consecución de recursos para el desarrollo de este trabajo, respectivamente.
Referencias
Referencias bibliográficas
- Abuchaibe, H., Agudelo, J., & Sañudo, C. (1998). Descripción general de las instalaciones del acueducto metropolitano. Revista Empresas Públicas de Medellín, 10, 35-97.
-
Azevedo, S. M. F. O. & Carmouze, J. (1994). Une mortalité de poissons dans une lagune tropicale (Brésil) durant une période de dominance de cyanophyceae. Coincidence ou conséquence?. Revue d’Hydrobiologie Tropicale, 27, 265-272. -
Azevedo, S. M., Carmichael, W. W., Jochimsen, E. M., Rinehart, K. L., Lau, S., Shaw, G. R., & Eaglesham, G. K. (2002). Human intoxication by microcystins during renal dialysis treatment in Caruaru-Brazil. Toxicology, 181-182, 441-446. -
Bittencourt-Oliveira, M. C., Oliveira, M. C., & Bolch, C. J. S. (2001). Genetic variability of some Brazilian strains of Microcystis aeruginosa complex (Cyanophyceae/Cyanobacteria) using the nucleotide sequence analysis of the intergenic spacer and flanking regions from cpcBA-phycocyanin operon. Journal of Phycology, 37, 810-818. -
Boutte, C., Grubisic, S., Balthasart, P., & Wilmotte, A. (2006). Testing of primers for the study of cyanobacterial molecular diversity by DGGE. Journal of Microbiological Methods, 65, 542-550. -
Carmichael, W. W., Azevedo, S. M. F. O., An, J. S., Molica, R. J. R., Jochimsen, E. M., Lau, S., Rinehart, K. L., Shaw, G. R., & Eaglesham, G. K. (2001). Human fatalities from cyanobacteria: chemical and biological evidence for cyanotoxins. Environmental Health Perspectives, 109, 663-668. -
Codd, G. A., Steffensen, D. A., Burch, M. D., & Baker, P. D. (1994). Toxic blooms of cyanobacteria in Lake Alexandrina, South Australia - learning from history. Australian Journal of Marine and Freshwater Research, 45, 731-736. -
Corporación Autónoma Regional del Centro de Antioquia - CORANTIOQUIA. (1997). Conservación, ordenamiento y manejo del sistema de páramos y bosques altoandinos del noroccidente medio antioqueño.(Tomo I: Diagnóstico biofísico, socio-económico y socio-cultural). Corantioquia, Medellín, Colombia. -
Correa, I. (2008). Toxicidad de florecimientos de cianobacterias en el embalse Riogrande II.(Tesis de maestría). Universidad de Antioquia, Medellín, Colombia. -
Dadheech, P. K., Krienitz, L., Kotut, K., Ballot, A., & Casper, P. (2009). Molecular detection of uncultured cyanobacteria and aminotransferase domains for cyanotoxin production in sediments of different Kenyan lakes. FEMS Microbiology Ecology, 68, 340-350. -
Departamento Administrativo Nacional de Estadística – DANE. (2005). Censo General 2005, Nivel Nacional.Bogotá, Colombia. -
Dawson, R. M. (1998). The toxicology of microcystins. Toxicon, 36(7), 953-962. -
Domingos, P., Rubim, T. K., Molica, R. J. R., Azevedo, S. M. F. O., & Carmichael, W. W. (1999). First report of microcystin production by picoplanktonic cyanobacteria isolated from a northeast Brazilian drinking water supply. Environmental Toxicology, 14, 31-35. -
Dorador, C., Vila, I., Imhoff, J., & Witzel, K. P. (2008). Cyanobacterial diversity in Salar de Huasco, a high altitude saline wetland in northern Chile: an example of geographical dispersion? FEMS Microbiology Ecology, 64, 419-432. -
Dos Anjos, M., Bittencourt-Oliveira, M., Zajac, M., Hiller, S., Christian, B., Erler, K., Luckas, B., & Pinto, E. (2006). Detection of harmful cyanobacteria and their toxins by both PCR amplification and LC-MS during a bloom event. Toxicon, 48, 239-45. -
Fewer, D. P., Halinen, K., Sipari, H., Bernardová, K., Mänttäri, M., Eronen, E., & Sivonen, K. (2011). Nonautonomous transposable elements associated with inactivation of microcystin gene clusters in strains of the genus Anabaena isolated from the Baltic Sea. Environmental Microbiology Reports, 3(2), 189-194. -
Gaevsky, N. A., Kolmakov, V. I., Belykh, O. I., Tikhonova, I. V., Joung, Y., Ahn, T. S., Nabatova, V. A., & Gladkikh, A. S. (2011). Ecological development and genetic diversity of Microcystis aeruginosa from artificial reservoir in Russia. The Journal of Microbiology, 49(5), 714-720. -
Gehringer, M. M., Adler, L., Roberts, A. A., Moffitt, M. C., Mihali, T. K., Mills, T. J., Fieker, C., & Neilan, B. A. (2012). Nodularin, a cyanobacterial toxin, is synthesized in planta by symbiotic Nostoc sp. The ISME Journal, 6(10), 1834-1847. -
Genuário, D. B., Silva-Stenico, M. E., Welker, M., Moraesc, L. A., & Fiore, M. F. (2010). Characterization of a microcystin and detection of microcystin synthetase genes from a Brazilian isolate of Nostoc. Toxicon, 55, 846-854. -
Hall, T.A. (1999). Bioedit. A user friendly biological sequence aligment editor and analysis program for windows 95/98/NT. Nucleic Acids Symposium Series 41: 95-98. -
Havens, K. E., James, R. T., East, T. L., & Smith, V. H. (2003). N:P ratios, light limitation, and cyanobacterial dominance in a subtropical lake impacted by non-point source nutrient pollution. Environmental Pollution, 122(3), 379-90. -
Hisbergues, M., Christiansen, G., Rouhiainen, L., Sivonen, K., & Börner, T. (2003). PCR-based identification of microcystin-producing genotypes of different cyanobacterial genera. Archives of Microbiology, 180, 402-410. -
Hitzfeld, B. C., Hoger, S. J., & Dietrich, D. R. (2000). Cyanobacterial toxins: removal during drinking water treatment, and human risk assessment. Environmental Health Perspectives, 108(1), 113-122. -
Hotto, A. M., Satchwell, M. F., & Boyer, G. L. (2007). Molecular characterization of potential microcystinproducing cyanobacteria in Lake Ontario embayments and nearshore waters. Applied and Environmental Microbiology, 73(14), 4570-4578. -
Jochimsen, E. M., Carmichael, W. W., An, J., Cardo, D., Cookson, S. T., Holmes, C. E. M., & Jarvis, W. R. (1998). Liver failure and death following exposure to microcystin toxins at a hemodialysis center in Brazil. The New England Journal of Medicine, 36, 373-378. -
Kaasalainen, U., Fewer, D. P., Jokela, J., Wahlsten, M., Sivonen, K., & Rikkinen, J. (2012). Cyanobacteria produce a high variety of hepatotoxic peptides in lichen simbiosis. Proceedings of the National Academy of Sciences of the United States of America, 109(15), 5886-5891. -
Lorenzi, A. S., Andreote, F. D., Furtado, A. L., Agujaro, L. F., Machado, J. B., Figueira, A., & Fiore, M. F. (s.f.). Development and field testing of a LUX real-time PCR for quantification of Microcystis and potential microcystin producers. Recuperado de http://www.ncbi.nlm.nih.gov/nuccore/FJ807671.1
» http://www.ncbi.nlm.nih.gov/nuccore/FJ807671.1 -
Magana-Arachchi, D., Wanigatunge, R., & Jeyanandarajah, P. (2008). Setting up a polymerase chain reaction assay for the detection of toxic cyanobacteria. Journal of the National Science Foundation of Sri Lanka, 36(3), 229-233. -
Matsunaga, T., Takeyama, H. & Nakayama, H. (2001). 16S rRNA-Targeted identification of cyanobacterial genera using oligonucleotide-probes immobilized on bacterial magnetic particles. Journal of Applied Phycology, 13, 389-394. -
Moisander, P. H., Lehman, P. W., Ochiai, M., & Corum, S. (2009). Diversity of Microcystis aeruginosa in the Klamath River and San Francisco Bay delta, California, USA. Aquatic Microbial Ecology, 57, 19-31. -
Montoya, Y. & Ramírez, J. J. (2007). Flujos de mineralización en el embalse tropical Río Grande II (Antioquia, Colombia). Limnetica, 26(1), 39-51. -
Muyzer, G., De Waal, E. C., & Uitterlinden, A. G. (1993). Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction amplified genes coding for 16S rRNA. Applied and Environmental Microbiology, 59(3), 695-700. -
Nei, M. & Li, W. H. (1979). Mathematical models for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences of the United States of America, 76(10), 5269-5273. -
Neilan, B., Jacobs, D., & Goodman, A. (1995). Genetic diversity and phylogeny of toxic cyanobacteria determined by DNA polymorphisms within the phycocyanin locus. Applied and Environmental Microbiology, 61, 3875-3883. -
Ngwa, F., Madramootoo, C., & Jabaji, S. (2012). Monitoring toxigenic Microcystis strains in the Missisquoi bay, Quebec, by PCR targeting multiple toxic gene loci. Environmental Toxicology. doi: 10.1002/tox.21770 -
Nikolausz, M., Sipos, R., Revesz, S., Szekely, A., & Marialigeti, K. (2005). Observation of bias associated with re-amplification of DNA isolated from denaturing gels. FEMS Microbiology Letters, 244, 385-390. -
Nübel, U., Garcia-Pichel, F., & Muyzer, G. (1997). PCR primers to amplify 16S rRNA genes from cyanobacteria. Applied and Environmental Microbiology, 63(8), 3327-3332. -
Otaya, L., Vásquez, G., & Bustamante, G. (2008). Estimación de la oferta hídrica con información escasa en ecosistemas estratégicos. Revista Facultad Nacional de Agronomía Medellín, 61(1), 4366-4380. -
Ouellette, A., Handy, S., & Wilhelm, S. (2006). Toxic microcystis is widespread in Lake Erie: PCR detection of toxin genes and molecular characterization of associated cyanobacterial communities. Microbial Ecology, 51, 154-165. -
Pope, P. B. & Patel, B. K. C. (2008). Metagenomic analysis of a freshwater toxic cyanobacteria bloom. FEMS Microbiology Ecology, 64, 9-27. -
Rademaker, J. L. W. & de Bruijn, F. J. (2004). Computerassisted analysis of molecular fingerprint profiles and database construction. In G. A. Kowalchuk, F. J. De Bruijn, I. M. Head, A. D. L. Akkermans, & J. D. Van Elsas (Eds.), Molecular Microbial Ecology Manual(pp. 1397-1446). Dordrecht, Nederlands: Springer. -
Ramírez, J. J., Gutiérrez, F. L., & Vargas, A. (2005). Respuesta de la comunidad fitoplanctónica a experimentos de eutrofización artificial realizados en la represa La Fe, El Retiro, Antioquia, Colombia. Caldasia, 27(1), 103-115. -
Rantala, A., Rajaniemi-Wacklin, P., Lyra, C., Lepisto, L., Rintala, J., Mankiewicz-Boczek, J., & Sivonen, K. (2006). Detection of microcystin-producing cyanobacteria in Finnish Lakes with genus-specific microcystin synthetase gene E (mcyE) PCR and associations with environmental factors. Applied and Environmental Microbiology, 72(9), 6101-6110. -
Rodger, H. D., Turnbull, T., Edwards, C., & Codd, G. A. (1994). Cyanobacterial (blue-green-algal) bloom associated pathology in brown trout Salmo trutta L. in Loch Leven, Scotland. Journal of Fish Diseases, 17, 177-181. -
Romaní, A., Artigas, J., Camacho, A., Graça, M., & Pascoal, C. (2009). La biota de los ríos: los microorganismos heterotróficos. In Conceptos y técnicas en ecología fluvial.España: Fundación BBVA. -
Runnegar, M. T., Kong, S., & Berndt, N. (1993). Protein phosphatase inhibition and in vivo hepatotoxicity of microcystins. The American Journal of Physiology, 265(2), G224-30. -
Sanguinetti, C. J., Dias Neto, E., & Simpson, A. J. G. (1994). Rapid silver staining and recovery of PCR products separated on polyacrylamide gels. Biotechniques, 17, 915-919. -
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., & Kumar, S. (2011). MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution, 28, 2731-2739. -
Tencalla, F. & Dietrich, D. (1997). Biochemical characterization of microcystin toxicity in rainbow trout (Oncorhynchus mykiss). Toxicon, 35, 583-595. -
Vaitomaa, J., Rantala, A., Halinen, K., Rouhiainen, L., Tallberg, P., Mokelke, L., & Sivonen, K. (2003). Quantitative real-time PCR for determination of microcystin synthetase E copy numbers for Microcystis and Anabaena in lakes. Applied and Environmental Microbiology, 69, 7289-7297. -
World Health Organization - WHO. (2004). Guidelines for drinking-water quality.Geneva, Switzerland: World Health Organization.
*Correspondencia a:
Julio César Hurtado-Alarcón
Julio César Hurtado-Alarcón
. Grupo BIODES-Biotecnología del Desarrollo Sustentable – Laboratorio de Biología Molecular y Celular, Universidad Nacional de Colombia Sede Medellín, Calle 59A Nº 63-20 Bloque 19A-310 Medellín, Colombia; jchurtad@unal.edu.co
Jaime Polanía-Vorenberg
. Posgrados en Bosques y Conservación Ambiental, Departamento de Ciencias Forestales, Universidad Nacional de Colombia Sede Medellín; jhpolaniav@unal.edu.co
1. Grupo BIODES-Biotecnología del Desarrollo Sustentable – Laboratorio de Biología Molecular y Celular, Universidad Nacional de Colombia Sede Medellín, Calle 59A Nº 63-20 Bloque 19A-310 Medellín, Colombia; jchurtad@unal.edu.co
2. Posgrados en Bosques y Conservación Ambiental, Departamento de Ciencias Forestales, Universidad Nacional de Colombia Sede Medellín; jhpolaniav@unal.edu.co
Recibido 03-IX-2012. Corregido 10-X-2013. Aceptado 13-Xi-2013.
Fechas de Publicación
-
Publicación en esta colección
09 Nov 2015 -
Fecha del número
Mar 2014
Histórico
-
Recibido
03 Set 2012 -
Revisado
10 Oct 2013 -
Acepto
13 Nov 2013
