Diversidad genética de Lilaea scilloides (Juncaginaceae) en el centro de México

Genetic diversity of Lilaea scilloides (Juncaginaceae) in central Mexico

Fabiola Magallán¹, Mahinda Martínez¹*, Luis Hernández-Sandoval¹, Antonio González-Rodríguez² y Ken Oyama^2,3

¹ Universidad Autónoma de Querétaro. Av. de las Ciencias s/n, 76230, Juriquilla, Querétaro, México. *mahinda@uaq.mx

³ Escuela Nacional de Estudios Superiores. Unidad Morelia, Universidad Nacional Autónoma de México. Antigua carretera a Pátzcuaro Núm. 8701 Col. Ex-Hacienda de San José La Huerta, 58190 Morelia, Michoacán, México.

Recibido: 17 agosto 2010
Aceptado: 03 octubre 2012

Resumen

Lilaea scilloides (Juncaginaceae) es una especie acuática de humedales temporales. Es una herbácea anual emergente que va desde Canadá hasta Argentina. Tiene reproducción sexual y asexual y es polinizada por viento. Su heteromorfismo floral es notorio con 5 diferentes tipos de flores. Se planteó la hipótesis de que, debido a su polinización, amplio intervalo de distribución y abundancia local, tendría alta diversidad genética, baja endogamia, baja diferenciación genética y alto flujo génico. Los objetivos de este estudio fueron determinar la diversidad y estructura genética de L. scilloides y compararlas con especies con características ecológicas y de historia de vida similares. Se analizaron 40 individuos en cada una de las 12 poblaciones con la técnica de electroforesis de isoenzimas en geles de almidón. Los resultados muestran baja diversidad genética dentro de las poblaciones (A= 1.77, Ae= 1.15, P= 56.6, Ho= 0.064 y He= 0.103), altos índices de endogamia (F_IS= 0.288), baja diferenciación genética entre las poblaciones (F_ST= 0.077) y alto flujo génico (Nm= 4.42). Se encontró evidencia de altos niveles de propagación clonal. Los resultados indican que la reproducción asexual es extremadamente importante, así como la posibilidad de que L. scilloides presente hidrofilia, un tipo de polinización no registrado para la especie.

Palabras clave: anemofilia, humedales temporales, Lilaea, plantas acuáticas.

Lilaea scilloides (Juncaginaceae) is an aquatic species inhabiting temporary wetlands. It is an annual herbaceous emergent plant distributed from Canada to Argentina. This species reproduces both sexually and asexually, and is wind-pollinated. Flowers are highly heteromorphic with 5 different types of flowers. We hypothesized that, because of its pollination syndrome, wide distribution and local abundance, the species would have high genetic diversity, low endogamy, low genetic differentiation, and high gene flow. The objectives of this study were to determine the genetic diversity and structure of L. scilloides and compare with species with similar ecological and life history traits. We analyzed 40 individuals from 12 populations using isozyme electrophoresis in starch gels. The results showed a low within population genetic diversity (A= 1.77, Ae= 1.15, P= 56.6, Ho= 0.064 and He= 0.103), a high endogamy index (F_IS= 0.288), low genetic differentiation among populations (F_ST= 0.077), and high gene flow (Nm= 4.42). We also found evidence of a high degree of clonal propagation in the species. These results indicate that asexual reproduction is extremely important for the species, as well as the possibility that Lilaea scilloides has hydrophilous pollination, which has not been previously reported.

Key words: aneomophily, temporary wetlands, Lilaea, aquatic plants.

Introducción

La diversidad y estructura genética están correlacionadas con las características ecológicas y de historia de vida de las plantas (e.g., Hamrick et al., 1979; Loveless y Hamrick, 1984; Hamrick y Godt, 1989; Hamrick y Godt, 1996). Sin embargo, las generalizaciones propuestas por estos autores se basan principalmente en estudios de plantas terrestres. Algunos estudios de plantas acuáticas muestran que en la categoría de especie generalmente se presenta mayor diversidad genética que en las plantas terrestres; por el contrario, a nivel de población presentan menor diversidad genética (Barrett et al., 1993; Laushman, 1993; Elam, 1998). La alta diversidad genética de las plantas acuáticas a nivel de especie se debe a los altos niveles de diferenciación genética entre poblaciones (Edwards y Sharitz, 2000; Chen et al., 2008; Wang et al., 2010), lo cual está relacionado con la distribución de sus hábitats, que generalmente se encuentran aislados dentro de paisajes terrestres (Elam, 1998). Esta distribución disminuye la probabilidad de flujo génico entre las poblaciones y aumenta la diferenciación a través de la deriva génica (Barrett et al., 1993; Elam, 1998; Chen et al., 2008). Otra característica de las plantas acuáticas es la diversidad clonal, que tiende a ser más alta en las poblaciones que presentan múltiples morfos florales (Shibayama y Kadono, 2007). Sin embargo, aún no se observan patrones claros para las especies de plantas acuáticas y es probable que su estructura y diversidad genética tenga mayor relación con su historia de vida particular, que con el hecho de vivir en ambientes acuáticos (Magallán et al., 2009; Honnay et al., 2010).

En el centro de México se ha descrito un tipo de humedales temporales con una alta diversidad florística (Martínez y García, 2001). Estos humedales, se distribuyen a manera de parches, lo cual es un aspecto de este tipo de ecosistemas que los hace atractivos para realizar estudios de genética de poblaciones y biología de la conservación (Elam, 1998). Estos humedales presentan una estacionalidad marcada, con un período de inundación de 6 a 8 meses y uno de sequía de 4 a 6 meses; se forman en depresiones poco profundas de suelos impermeables que se llenan únicamente por agua de lluvia, sin drenaje de ríos o arroyos, por lo que se pueden considerar como ecosistemas aislados. El período de inundación coincide con condiciones de temperatura ideal para el crecimiento de las plantas y el nivel de agua disminuye hasta secarse completamente entre diciembre y marzo dependiendo de la precipitación anual.

Lilaea scilloides (Poir.) Hauman (Juncaginaceae) es una monocotiledónea acuática que habita este tipo de humedales. Es una planta enraizada emergente anual con distribución desde Canadá hasta Argentina (Lot et al., 1999; Novelo y Lot, 2001). Esta especie ocupa una posición clave en la evolución floral de las monocotiledóneas (Posluszny et al., 1986). Presenta un notable heteromorfismo floral pues tiene 5 diferentes tipos de flores, 4 en la inflorescencia y una más en la base del escapo. La inflorescencia es una espiga con: i) las flores masculinas unisexuales en la punta de la espiga, ii) unas flores bisexuales con estilo corto en la parte media de la inflorescencia, iii) unas flores femeninas unisexuales con perianto y estilo corto en la base de la espiga y iv) unas flores femeninas unisexuales sin perianto con estilo corto. La quinta flor está en la base del escapo y es femenina con un estilo filamentoso muy largo (Posluszny et al., 1986). Lilaea scilloides es una especie con el síndrome de polinización anemófila debido a su morfología floral (Cook, 1988; Lot et al., 1999). Lilaea scilloides presenta 2 tipos de aquenios, los de la inflorescencia tienen 1 ala dorsal y 2 laterales, mientras que los de la base de la inflorescencia tienen estructuras parecidas a cuernos en su parte superior y no tienen alas (Lot et al., 1999; Novelo y Lot, 2001); la dispersión es hidrócora y ornitócora (Cook, 1988; Cook, 1996; Lot et al., 1999).

Los objetivos de este estudio fueron determinar la diversidad y estructura genética de las poblaciones de L. scilloides y compararlas con los patrones reportados para especies con características ecológicas y de historia de vida similares. Debido a su polinización anemófila, su amplia distribución y su abundancia en las poblaciones locales se hipotetizó que L. scilloides tendría alta diversidad genética, bajos índices de endogamia, baja diferenciación genética y alto flujo génico interpoblacional.

Se recolectaron 480 individuos de L. scilloides en 12 poblaciones; de éstas, 11 se distribuyen en Querétaro y 1 en Guanajuato, en el centro de México (Fig. 1). Las plantas se colectaron de forma aleatoria separadas entre sí al menos 1 m, se transportaron a 4° C y en el laboratorio se almacenaron en bolsas de plástico a -70° C.

Se usó la técnica de electroforesis horizontal de isoenzimas en geles de almidón (Cheliak y Pitel, 1984; Werth, 1985; Wendel y Weeden, 1989). Para la extracción de enzimas se usó el buffer "Microbuffer" (Werth, 1985), el procedimiento consistió en macerar un fragmento de hoja con 2 gotas del buffer de extracción y una pizca de arena esterilizada, con el extracto obtenido, se impregnaron los "wicks" hechos con papel de cromatografía Whatman 3MM. Los "wicks" se almacenaron a -70° C para ser utilizados al siguiente día. Los geles se prepararon a concentraciones del 12% de almidón (Wendel y Weeden, 1989). Se usaron 8 enzimas en 2 sistemas: C (APX, GOT, ME y PGI) y D (GDH, G-6PD, PGD y PGM). El sistema C se corrió a 65 mA durante 8 horas y el sistema D se corrió a 30 mA durante 9 horas. Los geles se tiñeron a 37° C. La nomenclatura de los loci y de los alelos se estableció de la siguiente manera: a los loci que estuvieran más cerca del ánodo se les asignó el número 1, los que seguían el número 2 y así sucesivamente conforme se acercaran al cátodo, se usó el mismo criterio para los alelos, asignando letras. Al momento de introducir los "wicks" en el gel se incluyó un "wick" testigo con genotipo conocido para asignar correctamente la nomenclatura de los loci y alelos. Se usaron "wicks" impregnados con amido black para determinar el frente de migración. La información de la estructura cuaternaria de las enzimas, que sirvió para llevar a cabo la adecuada interpretación de las bandas, se obtuvo de Wendel y Weeden (1989). Se llevó a cabo una primera lectura de los patrones de bandeo con los geles en fresco; éstos se fijaron en alcohol al 50% o en una solución 5:5:1 de metanol, agua y ácido acético.

Para los análisis estadísticos, se obtuvieron los genotipos individuales para cada población y se conformaron matrices de datos. Se utilizó el programa FSTAT (Goudet, 1995) para estimar en cada una de las poblaciones muestreadas el número promedio de alelos por locus (A), el número efectivo de alelos por locus (Ae), el porcentaje de loci polimórficos (P), la heterocigosidad observada (Ho), la heterocigosidad esperada (He) y el índice de fijación (F). La significancia de las desviaciones con respecto al equilibrio Hardy-Weinberg para cada locus en cada población se determinó mediante una prueba exacta, y la significancia de los valores globales mediante la prueba U, implementadas en el programa GENEPOP (Raymond y Rousset, 1995).

La base de datos se analizó para identificar individuos con genotipos multilocus idénticos (posibles clones), mediante el programa GENALEX 6.2 (Peakall y Smouse, 2006). Para cada genotipo multilocus repetido se calculó el valor de Psex, que es la probabilidad de observar n repeticiones de un genotipo dado en una muestra de tamaño N mediante reproducción sexual en una población con apareamientos aleatorios.

La diferenciación genética entre las poblaciones se estimó utilizando el método de Weir y Cockerham (1984) para calcular los estadísticos F de Wright usando el programa GENETIX 4 (Belkhir et al., 2004). El flujo génico entre poblaciones se calculó a partir del valor promedio de Fst mediante la fórmula Nm= (1 - Fst) / 4 Fst (Wright, 1951).

Adicionalmente, se realizó un análisis de varianza molecular (AMOVA) mediante el programa ARLEQUIN 3.11 (Excoffier et al., 2005), el cual subdivide la varianza molecular total en sus componentes inter-poblacional e intra-poblacional, mediante el cálculo de los estadísticos Φ, que son análogos a los estadísticos F de Wright. La significancia de cada uno de los componentes de la varianza se estimó a partir de distribuciones generadas mediante 10 000 permutaciones al azar.

La estructura genética también se evaluó mediante el análisis de asignación Bayesiano implementado en el programa STRUCTURE ver. 2.3.1 (Falush et al., 2003). En este análisis los individuos son asignados probabilística-mente a uno de K grupos genéticos, donde K puede tomar sucesivos valores desde K= 1 hasta un número definido por el usuario. Al final del análisis se determina el valor de K que tiene mayor probabilidad dados los datos (Falush et al., 2003). En este caso, se probaron valores desde K= 1 hasta K= 12 y se realizaron 10 corridas para cada valor de K. Se utilizó el modelo que permite mezcla entre los grupos genéticos y frecuencias alélicas correlacionadas entre poblaciones. No se proporcionó al programa información previa sobre la población de origen de los individuos. Cada corrida tuvo un "burn-in" de 500 000 pasos, seguidos de 10⁶ iteraciones. El número de grupos genéticos más probable se encontró determinando el valor máximo del estadístico ΔK, según el método de Evanno et al. (2005).

Finalmente, las relaciones de similitud entre las poblaciones se analizaron mediante la identidad genética (I) de Nei (1978) entre todos los pares posibles de poblaciones.

Se registraron un total de 45 alelos en los 14 loci definidos. Todos los loci fueron polimórficos excepto PGD-1. El número promedio de alelos por locus (A), incluyendo los 14 loci, fue de 1.77 con un intervalo de 1.29 a 2.07, mientras que el promedio del número efectivo de alelos por locus (Ae) fue de 1.15 con un intervalo de 1.09 a 1.25. En todas las poblaciones Ae fue menor que A, lo cual indica presencia de alelos raros. El porcentaje de loci polimórficos (P) fue de 56.6% con un intervalo de 28.6 a 71.4%. El promedio de heterocigosidad observada (Ho) fue de 0.064 con un intervalo de 0.009 a 0.115 y el promedio de heterocigosidad esperada (He) fue de 0.103 con un intervalo de 0.057 a 0.156. Estos valores se muestran en el cuadro 1.

Se obtuvieron 95 valores de índice de fijación o endogamia (F) en los 13 loci polimórficos de las 12 poblaciones estudiadas, de los cuales 30 fueron significativamente diferentes de cero y positivos, indicando exceso de individuos homocigotos, y 2 fueron significativamente diferentes de cero y negativos, indicando exceso de individuos hetero-cigotos. Las poblaciones Pedro Escobedo 1, Huimilpan 6, Huimilpan 7, San Juan 9, Amealco 10 y Amealco 11 no presentan valores significativamente diferentes de cero, indicando que estas poblaciones se encuentran en equilibrio de Hardy-Weinberg, mientras que las poblaciones Pedro Escobedo 2, Huimilpan 3, Huimilpan 4, Huimilpan 5, Huimilpan 8 y J. Rosas 12 presentan valores positivos y significativamente diferentes de cero, indicando endogamia en estas poblaciones. Los valores promedio de F por población se muestran en el cuadro 1.

El análisis de identificación de clones mostró un total de 182 genotipos multilocus únicos y 52 genotipos repetidos, cada uno entre 2 y 78 veces. Las población que mostró mayor presencia de clonalidad es Pedro Escobedo 1 (33 individuos tienen al menos un clon), mientras que las que tuvieron menor clonalidad fueron las poblaciones Humilpan 4 y J. Rosas 12 (con 8 individuos con algún clon) (Cuadro 1). Es interesante destacar que los genotipos multilocus con mayor número de repeticiones estuvieron presentes en varias poblaciones y no en una sola. Por ejemplo, el genotipo repetido más abundante se encontró en todas las poblaciones, excepto en la población 12. En todos los casos la probabilidad de repetición de los genotipos multilocus a través de reproducción sexual (Psex) fue baja (> 0.05).

Los estadísticos F de Wright: F_IS F_IT y F_ST calculados mediante el método de Weir y Cockerham (1984) se muestran en el cuadro 2. Las estimaciones de F_IS fueron positivas y significativamente diferentes de cero en 8 loci, es decir, en 61.53% de los loci polimórficos existe exceso de homocigotos. El promedio de F_IS fue de 0.288 y significativamente diferente de cero, indicando que existe endogamia al nivel de los individuos dentro de las poblaciones locales. Las estimaciones de F_IT fueron positivas y significativamente diferentes de cero en 8 loci, indicando que 61.53% de los loci polimórficos presentan exceso de homocigotos. El promedio de F_IT presentó un valor de 0.342 y fue significativamente diferente de cero. El promedio de F_ST tuvo un valor de 0.077 y fue significativamente diferente de cero, lo cual indica diferenciación genética moderada entre las poblaciones (Hartl y Clark, 1997).

El flujo génico entre poblaciones (Nm) presentó un valor de 4.421 (calculado a partir del promedio de F_ST). Como el valor de Nm es mayor a uno se considera que existe alto flujo génico entre las 12 poblaciones locales estudiadas de L. scilloides (Cuadro 3).

El análisis de varianza molecular (AMOVA) indicó que el 15% de la varianza genética total se encuentra entre poblaciones (Φ_st= 0.154; p= 0.01), mientras que el 85% de la variación genética se encuentra dentro de las poblaciones.

Los resultados de STRUCTURE, analizados mediante el estadístico ΔK mostraron que K= 2 es el número de grupos genéticos más probable dados los datos (Figs. 2, 3). Ambos grupos genéticos se encuentran presentes en todas las poblaciones, aunque destaca que la población 12 (J. Rosas, Guanajuato) tiene la mayor proporción del grupo genético representado en rojo, mientras que la población 9 (San Juan, Querétaro) muestra una preponderancia del grupo genético representado en verde. La proporción de individuos que muestran evidencia de ancestría mezclada entre los 2 grupos genéticos es relativamente baja, lo cual puede deberse a una mayor frecuencia de reproducción clonal en comparación con la reproducción sexual, como indicó el análisis de los genotipos multilocus.

Las identidades genéticas de Nei (1978) tienen un intervalo de 0.973 a 1.000, con un promedio de I= 0.990, indicando que las 12 poblaciones locales estudiadas son muy parecidas genéticamente (Cuadro 4).

Discusión

Al comparar la diversidad genética a nivel de población de L. scilloides (A= 1.77, P= 56.60, He= 0.103) con el promedio registrado por Hamrick y Godt (1989) para especies terrestres con características similares, tales como especies anuales (A= 1.48, P= 30.20, He= 0.105), con distribución amplia (A= 1.72, P= 43.00, He= 0.159) y con reproducción sexual y asexual (A= 1.47, P= 29.40, He= 0.103), se observa que L. scilloides presenta niveles similares de diversidad genética. En cambio, si se compara con especies que presentan polinización por viento (A= 1.79, P= 49.70, He= 0.148) se observa que sus niveles de diversidad genética son más bajos, lo que parece indicar que existen otros factores, además de la anemofilia, que están determinando la diversidad genética de la especie. La diversidad genética dentro de las poblaciones de L. scilloides (A= 1.77, P= 56.60, He= 0.103) fue comparada con los estudios en especies acuáticas anemófilas, como Carex basiantha (A= 1.46, P= 40.21, He= 0.138) y Carex willdenowii (A= 1.50, P= 40.00, He= 0.148) (Ford et al., 1998) y Typha minima (A= 1.50, He= 0.146) (Galeuchet et al., 2002), observándose que L. scilloides presenta valores muy similares en el número de alelos por locus y porcentaje de loci polimórficos y más bajos en la heterocigosidad esperada, que no corresponde con lo esperado para especies anemófilas.

Lilaea scilloides se considera anemófila debido a las características florales de las espigas (Cook, 1988), sin embargo no se ha considerado el gran heteromorfismo floral de la especie, que permite la posibilidad de autogamia, anemofilia e hidrofilia. Esta especie presenta un tipo de flores femeninas que crecen en la base de la inflorescencia y que nunca emergen del agua, en algunas de las cuales observamos la presencia de tubos polínicos y frutos maduros (Magallán, 2004), por lo que es posible suponer que L. scilloides presenta también hidrofilia. Las especies hidrófilas generalmente exhiben modificaciones estructurales en su polen, entre las que destacan el polen de forma globular que forma cadenas filiformes o granos muy alargados, probablemente para maximizar la eficiencia del transporte (Les, 1988). El polen de L. scilloides es elíptico, tripantocolpado, de 40 jim en su eje ecuatorial y de 50 jim en su eje polar y ornamentación reticulada (Magallán, 2004) el cual difiere del polen globular típico de las especies hidrófilas por lo que se requieren estudios ecológicos de la polinización más detallados en L. scilloides para explicar la discordancia entre el síndrome de polinización y el tipo de polen.

Los resultados muestran que existe alta endogamia dentro de las poblaciones de L. scilloides (F_IS= 0.288), lo cual es raro en especies anemófilas, por lo que es probable que exista autogamia en las flores hermafroditas de la espiga, dado que se observó que el androceo y el gineceo de dichas flores alcanzan su madurez al mismo tiempo. También es posible que el polen de las flores de las espigas llegue a través del agua a las flores femeninas sumergidas. En resumen, los diferentes modos de polinización que presenta L. scilloides junto con su modo de propagación vegetativa están influyendo de manera importante en la estructura y diversidad genética de la especie a nivel intra e interpoblacional.

El valor de F_ST (0.077) señala que existe poca diferenciación genética entre poblaciones lo cual contrasta con los patrones sugeridos para plantas acuáticas (Edwards y Sharitz, 2000; Chen et al., 2008; Wang et al., 2010). Asimismo, los valores de identidad genética entre las poblaciones locales de L. scilloides (I= 0.990) señalan que las poblaciones locales son muy parecidas genéticamente. La alta similitud genética puede deberse al alto flujo génico (Nm= 4.42) o a cuestiones históricas (polimorfismos ancestrales). Todos los humedales en los que L. scilloides fue colectada son ecosistemas aislados, algunos muy cercanos, como es el caso de la poblaciones 1 a 8, que como máxima distancia de separación tienen 6 km. Por el contrario, otros humedales se ubican hasta 40 km de distancia, como la población 12 que se ubica en el estado de Guanajuato (Fig. 1). La poca diferenciación genética entre poblaciones podría explicarse, por un lado, por el amplio alcance que tiene la dispersión de las diásporas, a través del agua y/o incluso por aves (Novelo y Lot, 2001); así mismo, la dispersión de los frutos que presentan estructuras aladas y están situados en las espigas, pueden ser llevados por el viento a largas distancias. Por otro lado, es posible que las poblaciones locales fueran fundadas por un pequeño número de individuos provenientes de una sola fuente, que después tuvieron una extensiva reproducción clonal y que históricamente no han sufrido procesos de colonización - extinción importantes. La reproducción clonal es muy común en plantas acuáticas llegándose a casos extremos de monoclonalidad en ciertas poblaciones estables en contraste con poblaciones policlonales que han estado sujetas a disturbios por actividades humanas (Lambertini et al., 2008).

Dado que L. scilloides habita ecosistemas aislados, se esperaba que la deriva génica fuera un proceso determinante en su estructura genética, sin embargo, esta fuerza evolutiva tiene como consecuencia alta diferenciación genética entre poblaciones (Hartl y Clark, 1997), que no se observa en las poblaciones estudiadas, por lo que es posible que no esté afectando la estructura poblacional de la especie. Es probable que la presencia de un banco de semillas esté eliminando los efectos de la deriva génica, debido a que los humedales temporales de estudio pasan por períodos de sequía - inundación, es necesario que las plantas tengan mecanismos que les permiten sobrevivir de un periodo a otro (banco de semillas). La persistencia del banco de semillas amortigua los "cuellos de botella" y permite que la diversidad genética alta se mantenga en las poblaciones (Edwards y Sharitz, 2000)

Los resultados obtenidos en este estudio no corroboran las hipótesis de alta diversidad y baja endogamia basadas en su distribución y polinización anemófila. Por el contrario, L. scilloides presenta baja diversidad genética y alta endogamia, debido a que no se habían tomado en cuenta variables como la polinización por autogamia, la posible hidrofilia y la extensa propagación clonal de la especie. Otro aspecto contrastante es que la mayoría de los estudios en plantas acuáticas muestran alta diferenciación poblacional (Barrett et al., 1993; Laushman, 1993; Edwards y Sharitz, 2000; Chen et al., 2008; Wang et al., 2010), mientras que L. scilloides presenta poblaciones muy parecidas genéticamente. Las características florales poco comunes de L. scilloides la hacen un objeto de estudio muy interesante desde el punto de vista evolutivo, por lo que es recomendable llevar a cabo estudios de su biología floral, mecanismos de dispersión y estudios de diversidad genética incluyendo más poblaciones, ya que esto contribuiría ampliamente a su entendimiento. Estudios en plantas acuáticas con heterostilia muestran que se puede mantener una alta diversidad y estructura genética en las poblaciones con distintas proporciones de morfos reproductivos, si se presentan mecanismos que garanticen la reproducción cruzada en poblaciones adyacentes (Shibayama y Kadono, 2007).

Finalmente, estudios con marcadores moleculares de ADN en plantas acuáticas pueden dilucidar las relaciones genealógicas e históricas entre los haplotipos distribuidos en poblaciones a diferentes escalas espaciales (Koga et al., 2008; Tan et al., 2008; Fehrmann et al., 2012), distinguir procesos de hibridación e introgresión (Ito et al., 2010; Yoshino et al., 2011) y robustecer los estudios taxonómicos y filogenéticos (Dkhar et al., 2010; Na et al., 2010; Telford et al., 2011). Son muy pocos los estudios de esta naturaleza realizados en plantas acuáticas mexicanas.

Agradecimientos

Los autores agradecen a D. Cabrera y a un revisor anónimo por sus valiosos comentarios y sugerencias; a P. Balderas por su ayuda en campo y a A. Beltrán por su ayuda en el laboratorio. A R. Luna-Reyes y N. Pérez-Naser por sus enseñanzas en cuanto a la técnica de electroforesis. C. Rojas elaboró el mapa de distribución. Esta investigación fue financiada por el Consejo Nacional de Ciencia y Tecnología (CONACYT 29104N) y por el Consejo de Ciencia y Tecnología de Querétaro (CONCYTEQ).

Barrett, S., C. Eckert y B. Husband 1993. Evolutionary processes in aquatic plant populations. Aquatic Botany 44:105-145.         [ Links ]

Belkhir, K., P. Borsa, L. Chikhi, N. Raufaste y F. Bonhomme. 1996-2004. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Montpellier: Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II.         [ Links ]

Cheliak, W. y J. Pitel. 1984. Techniques for starch gel electrophoresis of enzymes from forest tree species. Information Report PI-X-42. Petawa National Forestry Institute, Berkeley, California.         [ Links ]

Chen, Y. Y., X. L. Li, L. Y. Yin y W. Li. 2008. Genetic diversity of the threatened aquatic plant Ottelia alismoides in the Yangtze river. Aquatic Botany 88:10-16.         [ Links ]

Cook, C. 1988. Wind pollination in aquatic angiosperms. Annals of the Missouri Botanical Garden 75:768-777.         [ Links ]

Cook, C. 1996. Aquatic Plant Book. SPB Academic Publ. 2a. Edición, Amsterdam. 220 p.         [ Links ]

Dkhar, J., S. Kumaria, S. R. Rao y P. Tandon. 2010. Molecular phylogenetics and taxonomic reassessment of four Indian representatives of the genus Nymphaea. Aquatic Botany 93:135-139.         [ Links ]

Edwards, A. y R. Sharitz. 2000. Population genetics of two rare perennials in isolated wetlands: Sagittaria isoetiformis and S. teres (Alismataceae). American Journal of Botany 87:1147-1158.         [ Links ]

Elam, D. 1998. Population genetics of vernal pool plants: theory, data and conservation implications. In Ecology, conservation and management of vernal pool ecosystems, C. Witham, E. Bauder, D. Belk, W. Ferren y R. Ornduff (eds.). California Native Plants Society. Sacramento. p. 180-189.         [ Links ]

Evanno, G., S. Regnaut y J. Goudet. 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Molecular Ecology 14:2611-2620.         [ Links ]

Excoffier, L., G. Laval y S. Schneider. 2005. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evolutionary Bioinformatics 1:47-50.         [ Links ]

Falush, D., M. Stephens y J. K. Pritchard. 2003. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164:15671587.         [ Links ]

Fehrmann, S., C. T. Philbrick y R. Halliburton. 2012. Intraspecific variation in Podostemum ceratophylum (Podostomaceae): evidence of refugia and colonization since the last glacial maximum. American Journal of Botany 99:145-151.         [ Links ]

Ford, B., D. McQueen, R. Naczi y A. Reznicek. 1998. Allozyme variation and genetic relationships among species in the Carex willdenowii complex (Cyperaceae). American Journal of Botany 85:546-552.         [ Links ]

Galeuchet, D., R. Holderegger, R. Rutishauser y J. Schneller. 2002. Isozyme diversity and reproduction of Thypa minima populations on the upper River Rhine. Aquatic Botany 74:19-32.         [ Links ]

Goudet, J. 1995. Fstat version 1.2: a computer program to calculate F statistics. Journal of Heredity 86:485-486.         [ Links ]

Hamrick, J., Y. Linhart y J. Mitton. 1979. Relationships between life history characteristics and electrophoretically detectable genetic variation in plants. Annual Review of Ecology and Systematics 10:173-200.         [ Links ]

Hamrick, J. y M. Godt. 1989. Allozyme diversity in plant species. In Plant population genetics, breeding, and genetic resources, A. D. Brown, M. T. Clegg, A. L. Kahler y B.S. Weir (eds.). Sinauer Associates, Sunderland, Massachusetts. p. 43-63.         [ Links ]

Hamrick, J. y M. Godt. 1996. Effects of history traits on genetic diversity in plant species. Philosophical Transactions of the Royal Society of London 351:1291-1298.         [ Links ]

Hartl, D. y A. Clark 1997. Principles of Population Genetics. Sinauer Associates, Sunderland, Massachusetts. 542 p.         [ Links ]

Hofstra, D., J. Clayton, J. Green y K. Adam. 2000. RAPD profiling and isozyme analysis of New Zeland Hydrilla verticillata. Aquatic Botany 66:153-166.         [ Links ]

Honnay, O., H. Jacquemyn, K. Nackaerts, P. Breyne y K. Van Looy. 2010. Patterns of population genetic diversity in riparian and acquatic plant species along rivers. Journal of Biogeography 37:1730-1739.         [ Links ]

Ito, Y., T. Ohi-Toma, J. Murata y N. Tanaka. 2010. Hybridization and poliploidy of an aquatic plant, Ruppia (Ruppiaceae), inferred from plastid and nuclear DNA phylogenies. American Journal of Botany 97:1157-1167.         [ Links ]

Koga, K., Y. Kadono y H. Setoguchi. 2008. Phylogeography of the Japanese water crowfoot based on chloroplast DNA haplotypes. Aquatic Botany 89:1-8.         [ Links ]

Lambertini, C., M. H. G. Gustaffson, J. Frydenberg, M. Speranza y H. Brix. 2008. Genetic diversity patterns in Phragmites australis at the population, regional and continental scales. Aquatic Botany 88:160-170.         [ Links ]

Laushman, R. 1993. Population genetics of hydrophilous angiosperms. Aquatic Botany 44:147-158.         [ Links ]

Les, D. 1988. Breeding systems, populations structure and evolution in hydrophylous angiosperms. Annals of the Missouri Botanical Garden 75:819-835.         [ Links ]

Les, D. 1991. Genetic diversity in the monoecious hydrophile Ceratophyllum (Ceratophyllaceae). American Journal of Botany 78:1070-1082.         [ Links ]

Lot, A., A. Novelo, M. Olvera y P. Ramírez-García. 1999. Catálogo de angiospermas acuáticas de México. Hidrófitas estrictas emergentes, sumergidas y flotantes. Cuaderno No. 33. Instituto de Biología, Universidad Nacional Autónoma de México. México. 161 p.         [ Links ]

Loveless, M. y J. Hamrick. 1984. Ecological determinants of genetic structure in plant populations. Annual Review of Ecology and Systematics 15:65-95.         [ Links ]

Magallán, F. 2004. Diversidad genética de plantas acuáticas en los charcos temporales de Amealco y Huimilpan, Querétaro. Tesis doctorado, Universidad Autónoma de Querétaro. Querétaro. 121 p.         [ Links ]

Magallán, F., M. Martínez, L. Hernández y K. Oyama. 2009. Estructura genética de poblaciones de Eriocaulon bilobatum (Eriocaulaceae): una especie amenazada de humedales temporales. Boletín de la Sociedad Botánica de México 85:81-88.         [ Links ]

Martínez, M. y A. García. 2001. Flora y vegetación acuáticas de localidades selectas del estado de Querétaro. Acta Botanica Mexicana 54:1-23.         [ Links ]

Na, H. R., C. Kim y H. K. Choi. 2010. Genetic relationships and genetic diversity among Typha taxa from East Asia based on AFLP markers. Aquatic Botany 92:207-213.         [ Links ]

Nei, M. 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583-590.         [ Links ]

Novelo, A. y A. Lot. 2001. Lilaeaceae. In Flora fanerogámica del Valle de México, G. Calderón y J. Rzedowski (eds.). Instituto de Ecología A.C., México, D. F. p. 990-993.         [ Links ]

Peakall, R. y P. E. Smouse. 2006. GENALEX 6, genetic analysis in Excel: population genetics software for teaching and research. Molecular Ecology Notes 6:288-295.         [ Links ]

Posluszny, U., W. Charlton y D. Jain. 1986. Morphology and development of the reproductive shoots of Lilaea scilloides (Poir.) Hauman (Alismatidae). Botanical Journal of the Linnean Society 92:323-342.         [ Links ]

Raymond, M. y F. Rousset. 1995. GENEPOP (Version 1.2): population genetics software for exact tests and ecumenicism. Journal of Heredity 86:248-249.         [ Links ]

Shibayama, Y. y Y. Kadono. 2007. Reproductive success and genetic structure of populations of the heterostylous aquatic plant Nymphoides indica (L.) Kuntze (Menyanthaceae). Aquatic Botany 86:1-8.         [ Links ]

Tan, B., K. Liu, X. L. Yue, F. Liu, J. M. Chen y Q. F. Wang. 2008. Chloroplast DNA variation and phylogeographic patterns in the Chinese endemic marsh herb Sagittaria potamogetifolia. Aquatic Botany 89:372-378.         [ Links ]

Telford, A., M. T. O'Hare, S. Cavers y N. Holmes. 2011. Can genetic bar-coding be used to identify Ranunculus L. subgenus Batrachium (DC) A. Gray? A test using some species from the British Isles. Aquatic Botany 95:65-70.         [ Links ]

Triest, L. 1989. Electrophoretic polymorphism and divergence in Najas marina L. (Najadaceae). Molecular markers for individuals, hybrids, cytodemes, lower taxa, ecodemes and conservation of genetic diversity. Aquatic Botany 33:301-380.         [ Links ]

Wang, Y. H., J. M. Chen, C. Xu, X. Li, Q. F. Wang y T. J. Motley. 2010. Population genetic structure of an aquatic herb Batrachium bungei (Ranunculaceae) in the Hengduan Mountains of China. Aquatic Botany 92:221-225.         [ Links ]

Weir, B. S. y C. C. Cockerham. 1984. Estimating F-Statistics for the analysis of population structure. Evolution 38:1358-1370.         [ Links ]

Wendel, J. y N. Weeden. 1989. Visualization and interpretation of plant isozymes. In Isozymes in Plant Biology, D. Soltis y P. Soltis (eds.). Dioscorides, Portland, Oregon. p. 5-45.         [ Links ]

Werth, C. 1985. Implementing an isozyme laboratory at a field station. Virginia Journal of Science 36:53-76.         [ Links ]

Wright, S. 1951. The genetical structure of populations. Annals of Eugenics 15:322-354.         [ Links ]

Yoshino, N., G. X. Wang, A. Uchino y T. Tominaga. 2011. Hybridizations and genetic relationships among Lindernia species (Scrophulariaceae): L. procumbens and two subspecies of L. dubia. Aquatic Botany 94:165-171.         [ Links ]