PROGRAMA FONDECYT

INFORME FINAL

ETAPA 2016

COMISIÓN NACIONAL DE INVESTIGACION CIENTÍFICA Y TECNOLÓGICA

VERSION OFICIAL Nº 3

FECHA: 06/06/2017

Nº PROYECTO : 3140610 DURACIÓN : 3 años AÑO ETAPA : 2016 TÍTULO PROYECTO : FILOGEOGRAFIA COMPARADA DE MOLUSCOS LITORALES DE LA COSTA CHILENA: DETERMINADO LA RELEVANCIA DE LAS ESTRATEGIAS DE DESARROLLO LARVAL EN SU HISTORIA EVOLUTIVA

DISCIPLINA PRINCIPAL : GENETICA Y EVOLUCION GRUPO DE ESTUDIO : BIOLOGIA 1 INVESTIGADOR(A) RESPONSABLE : MARIA CECILIA PARDO GANDARILLAS DIRECCIÓN : COMUNA : CIUDAD : SANTIAGO REGIÓN : METROPOLITANA

FONDO NACIONAL DE DESARROLLO CIENTIFICO Y TECNOLOGICO (FONDECYT) Moneda 1375, Santiago de Chile - casilla 297-V, Santiago 21 Telefono: 2435 4350 FAX 2365 4435 Email: [email protected] INFORME FINAL PROYECTO FONDECYT POSTDOCTORADO

OBJETIVOS

Cumplimiento de los Objetivos planteados en la etapa final, o pendientes de cumplir. Recuerde que en esta sección debe referirse a objetivos desarrollados, NO listar actividades desarrolladas. Nº OBJETIVOS CUMPLIMIENTO FUNDAMENTO 1 1. Objetivo: Estimar la diversidad y TOTAL El cumplimiento de este objetivo fue total estructuración genética de cinco especies de considerando las etapas y actividades en el plan moluscos a lo largo de su rango de distribución, a de trabajo en la propuesta del proyecto para través de distintos marcadores moleculares. ejecutar durante el año 2016. En el curso de la etapa 2016 se aumentó el numero de localidades abarcando todo el rango de distribución de cuatro especies propuestas en el proyecto. En esta última etapa se incorporó una cuarta especie, obteniendo resultados de la diversidad y estructuración genética, con el gen mitocondrial citocromo oxidasa I (COI) y del gen nuclear espaciador transcrito interno (ITS1) de esta última especie a lo largo de todo Chile. En este contexto se tuvo resultados de cuatro especies con el gen COI y con ITS1, de las cinco propuestas en este proyecto (la última especie fue imposible obtener buenos resultados de amplificación de DNA). De esta manera se pudo mostrar resultados comparativos de sus diversidad y estructura genética desde distintas localidades a lo largo de la costa del Pacifico sureste. 2 2.- Objetivo: Inferir la historia demográfica de las TOTAL El cumplimiento de este objetivo fue total según poblaciones de cada especie, a partir de los las etapas y actividades para el año en ejecución patrones de diversidad genética obtenidos con 2016 del proyecto. En el curso de la etapa 2016, distintos marcadores moleculares. se avanzó aumentando el número de localidades desde todo Chile y Perú de dos especie, y se obtuvieron nuevos datos de dos especies más para evaluar la historia demográfica, para finalmente tener resultados de cuatro especies (usando el gen mitocondrial citocromo oxidasa I, COI) (gen ITS1 fue poco diverso e imposible hacer análisis demográficos). Así se obtuvo la inferencia demográfica de cuatro de las cinco especies propuestas en el proyecto (la última especie no se pudo tener datos de amplificación de DNA), a partir de datos de diversidad genética. 3 3.- Objetivo: Realizar paleo-reconstrucciones del TOTAL El cumplimiento de este objetivo fue total, rango de distribución de las especies en el pasado. considerando en esto las etapas y actividades en el plan de trabajo para el año en ejecución 2016 en la propuesta del proyecto. Se tiene mapeado los puntos geo-referenciados de 30 puntos a lo largo de todo Chile y Perú como dato contemporáneo del rango de distribución de cuatro de las cinco especies propuesta en el proyecto. Para la paleo reconstrucción del rango de distribución en el pasado se generó modelos de nicho ecológico, usando datos de variables bioclimáticas oceánicas que fueron seleccionadas and bajadas desde website: Marine Spatial Ecology (MARSPEC) y desde Paleoclimatic Marine Spatial Ecology (Paleo-MARSPEC) dataset. De este modos se obtuvo la paleo-reconstrucción del rango de distribución de una especie. 4 4.- Objetivo: Evaluar si las poblaciones de las TOTAL El cumplimiento de este objetivo fue total, cinco especies de moluscos (con larva considerando en esto las etapas y actividades en lecitotrófica y planctotrófica) que habitan en la el plan de trabajo para el año 2016. Se tiene zona norte-centro y sur de Chile, muestran un resultados de cuatro de la cinco especies patrón distinto de diversidad y estructuración propuestas (la última no se pudo obtener genética. amplificación de DNA). Tegula atra y Chiton magnificus y Fissurella picta con larva lecitotrófica, y la especie lessonii con larva planctotrófica desde varias localidades a lo largo de Chile u Perú. Tegula atra (con larva lecitotrófica) fue la que mostró niveles más bajo de diversidad genética (con el gen COI e ITS1) seguido por Chiton magnificus y luego Fissurella picta (ambas con larva lecitotrófica) y finalmente con una muy alta diversidad la lapa Siphonaria lessonii (con larva planctotrófica). Las tres especies T. atra y C. magnificus y S. ñessonii mostraron una estructuración genética de dos unidades poblacionales, mientras que la especie F. picta (con larva lecitotrófica) no mostró significativa estructuración genética. 5 5.- Objetivo: Evaluar si las poblaciones de las TOTAL El cumplimiento de este objetivo fue total, cinco especies de moluscos que habitan en la considerando en esto las etapas y actividades en zona norte-centro y sur de Chile, muestran un el plan de trabajo del año 2016. Se tuvo patrón demográfico y de paleo-distribución resultados de la inferencia demográfica distintas distinto. de todas las poblaciones de cuatro de la cinco especies propuestas, desde la zona norte y centro sur de Chile y Perú. Las especies Tegula atra y Chiton magnificus, con larva lecitotrófica, mostraron un patrón demográfico más reciente o joven en as poblaciones del sur, mientras Siphonaria lessonii (con larva planctotrófica) y Fissurella picta (con larva lecitotrófica) mostraron una expansión mucho más antigua. Los modelos de paleo-distribución también mostraron un patrón distinto en las poblaciones del extremo norte y sur del rango de distribución, como es el caso de S. lessonii.

Otro(s) aspecto(s) que Ud. considere importante(s) en la evaluación del cumplimiento de objetivos planteados en la propuesta original o en las modificaciones autorizadas por los Consejos. RESULTADOS OBTENIDOS: Para cada uno de los objetivos específicos, describa o resuma los resultados. Relacione las publicaciones y /o manuscritos enviados a publicación con los objetivos específicos. En la sección Anexos incluya información adicional que considere pertinente para efectos de la evaluación. La extensión máxima de esta sección es de 5 páginas (letra tamaño 10, Arial o Verdana).

1. Objetivo: Estimar la diversidad y estructuración genética de cinco especies de moluscos a lo largo de su rango de distribución, a través de distintos marcadores moleculares. 1.1 Diversidad genética: 1.1. a Método utilizado: Se calcularon los índices de diversidad genética estándares como: número de haplotipos (K), número de sitios polimórficos (S), diversidad de haplotipos (Hd), número promedio de diferencias entre secuencias (), diversidad nucleotídica (), en cada localidad usando el programa DnaSP versión 5.0 (Librado & Rozas 2009).

Especie: Tegula atra, Resultados: Se encontraron 42 haplotipos en fragmentos de 663 pares de bases (pb), del gen mitocondrial gen citocromo oxidasa I (COI) de un total de 447 individuos Tegula atra de 16 localidades entre los 8° S y 39° S de latitud por el Pacifico sureste. Cada localidad mostró una diversidad haplotípica y Nucleotídica intermedia a baja, y pocos sitios polimórficos (Tabla 1 a). Las localidades hacia la zona central y sur mostraron los valores más altos de diversidad genética (Tabla 1 a). Con el gen nuclear ITS1 se encontró entre baja a intermedia diversidad genética en las nueve localidades entre los 8° S y 43° S. Cuatro haplotipos se encontraron en fragmentos de 263 pb de 118 individuos T. atra. Todos los índices de diversidad mostraron un incremento desde localidades del extremo norte hacia el sur (Tabla 1 b). Especie: Siphonaria lessonii, Resultados: Se encontraron 241 haplotipos con fragmentos de 676 pb del gen mitocondrial COI en los 561 individuos S. lessonii. En promedio se evidenció una muy alta diversidad haplotípica y nucleotídica, un número alto de sitios polimórficos y número de diferencias entre pares de secuencias en las 16 localidades muestreadas desde los -11° hasta -55°Sur por el Pacífico sureste. La diversidad genética mostró una cierta reducción latitudinal hacia el sur del rango de distribución (tabla 2 a). Mientras, el fragmento de 640 pb del gen nuclear denominado Espaciador transcrito interno (ITS 1) reveló una muy baja diversidad genética, constituida por 5 haplotipos únicamente. Además, solo dos localidades se encontró variabilidad muy reducida en diversidad haplotípica, nucleotídica, en el número sitios polimórficos y número de diferencias entre pares de secuencias en las 14 localidades entre los 18° a los 53°S de latitud (Tabla 2b). Se está finalizando un manuscrito para enviar pronto a publicar sobre: Molecular footprint and niche modelling reveal the effects of the ice ages on the pulmonate Siphonaria lessonii (: ) evolutionary lineages in Southeastern Pacific Ocean.

Especie: Fissurella picta, Resultados: Un Total de sesenta y cuatro haplotipos fueron encontrados en fragmentos de 584 pb del gen COI de un total de 106 individuos Fissurella picta desde 11 localidades desde los 27° S hasta los 54° S por el Pacifico Sureste. En cada localidad se evidenció una muy alta diversidad haplotípica y nucleotídica, un número alto de sitios polimórficos y número de diferencias entre pares de secuencias (Tabla 3 a). Todos los índices mostraron una leve reducción latitudinal hacia el sur de la diversidad genética. Con el gen nuclear ITS1 se encontraron 3 haplotipos en fragmentos de 340 pb de 88 individuos de F. picta desde siete localidades desde los 32° a los 54° S de latitud. Todos los índices mostraron una muy baja diversidad genética en la zona central y casi nula o cero hacia las localidades del norte y sur del rango de distribución (Tabla 3 b). Especie: Chiton magnificus, Resultados: Un Total de treinta y ocho haplotipos fueron encontrados en los fragmentos de 680 pb del gen COI en un total de 330 individuos Chiton magnificus desde 15 localidades desde los 28° a los 53° S por el Pacifico sureste. En cada localidad se evidenció una diversidad haplotípica y nucleotídica intermedia, intermedio número de sitios polimórficos y relativamente pocas diferencias entre pares de secuencias (Tabla 4 a). Todos los índices mostraron una tendencia latitudinal de significativa reducción de la diversidad genética hacia al sur del rango de distribución (Tabla 4 a). Con el gen nuclear ITS1 se encontraron 5 haplotipos en fragmentos de 334 pb de 150 individuos de C. magnificus de 13 localidades desde 28°S hasta 43°S de latitud. Todos los índices mostraron una muy baja o nula diversidad genética en casi todas las localidades del norte y sur del rango de distribución (Tabla 4 b). ______1.1. b Método utilizado: La reconstrucción filogenéticas con datos intra-específicos fue representado a través de una red de haplotipos (datos poblacionales sin recombinación), construida usando el algoritmo median-joining (MJ) para datos binarios implementado en el software Network versión 4.6.1.2 (Bandelt et al. 1999).

Especie: Tegula atra, Resultados: Las relaciones genealógicas con el gen COI, realizado mediante Median-Joining Network mostraron un patrón caracterizado por dos haplotipos dominantes separados por 1 paso mutacional, cada uno rodeado de haplotipos de baja frecuencia y únicos (Fig. 1 a). El primer haplotipo dominante es el enlace central de 25 haplotipos únicos, de 7 haplotipos de baja frecuencia, todos separados por uno o dos pasos mutacionales de. El segundo haplotipo más frecuente se conecta con siete haplotipos únicos, separados por uno y dos mutaciones, y de un haplotipo de baja frecuencia. El primer haplotipo dominante se encontró en 249 (59.4 %) de los 471 caracoles analizados de todas las localidades muestreadas desde 8° S hasta los 39° S de latitud. El segundo haplotipo dominante se encontró en 120 (28.6 %) de los 471 caracoles desde los 29° S hasta los 39°S (Fig. 1a). La genealogía con el gen ITS 1 mostró 4 haplotipos relativamente frecuentes en casi las 9 localidades, siendo el más frecuente encontrado en la zona norte y el centro del rango de distribución (Fig. 1 b). Especie Siphonaria lessonii, Resultados: Las genealogía de los 249 haplotipos con el gen COI mostró un patrón de cuatro haplotipos dominantes H49 seguido por H7, H24 y H161 (Fig. 2 a). El haplotipo dominante H49 se encontró rodeado por más de 70 haplotipos, entre ellos, uno de los que les sigue en dominancia H161, además de algunos haplotipos baja frecuencia y varios únicos unidos por uno, dos hasta cuatro pasos mutacionales. Luego lo siguiente haplotipos dominantes H7 y H24 se encontraron conectados por un paso mutacional, y de 8 y 9 pasos mutacionales del haplotipo más frecuente H49. Ambos haplotipos H7 y H24 se encontraron constituyente la base de la formando de una red genealógica de 179 haplotipos únicos y de baja frecuencia (Fig. 2 a). El haplotipo dominante H95 fue encontrado en 85 (15.6%) de las 546 lapas analizadas, mientras que los siguiente haplotipos fueron encontrados en menos de 25 (4.6%) de las lapas analizadas. La genealogía de los 5 haplotipos con el gen nuclear ITS 1 mostró dos haplotipos frecuentes H1 y H3 unidos por tres haplotipos únicos H2, H3 y H5, más dos pasos mutacionles (Fig. 2 b). De los 170 individuos analizados, 103 (60.6%) individuos tenía el haplotipo dominante H1, y 64 (37.6%) tenían el haplotipo H3. Especie Fissurella picta, Resultados: Las relaciones genealógicas de los 64 haplotipos de F. picta mostraron un patrón de tres haplotipos dominantes H29 el más frecuente seguido por el segundo haplotipo frecuente H19 y un tercero H13, todos unidos por un paso mutacional (Fig. 3 a). El haplotipo más frecuente H29 se encontró rodeado por 3 haplotipos únicos, el H19 rodeado por más de 10 haplotipos únicos y de baja frecuencia conectados por uno o dos paaso mutacionales. El haplotipos H13 se encontró rodeado por los restantes haplotipos que fueron más de veinte (Fig. 3 a). El haplotipo H29 se presentó en 17 (16.04 %) de los individuos del total, y el segundo haplotipo H19 en 15 individuos (14.15 %) analizados, y el tercer haplotipo H13 se encontró en 5 individuos (4.71%). Con el gen ITS 1, las relaciones genealógicas de los 3 haplotipos mostraron que el haplotipo H1 más frecuente o dominantes se conectó con un paso mutacional a los otros dos haplotipos H2 y H3 (Fig. 3 b). El haplotipo H1 fue encontrados en todas las localidades, mientras que H2 y H3 en dos y tres localidades. Especie Chiton magnificus, Resultados: Las relaciones genealógicas de los 38 haplotipos de C. magnificus mostraron un patrón de dos haplotipos dominantes H14 el más frecuente seguido por el segundo haplotipo frecuente H1, ambos unidos por un haplotipo H23 y dos pasos mutacionales (Fig. 4 a). El haplotipo más frecuente H14 se encontró rodeado por 14 haplotipos únicos (uno de ellos conectado mediante un median vector), y 3 haplotipos de baja frecuencia H19, H24 y H35. El segundo haplotipo frecuente H1 estuvo conectado por uno o más pasos mutacionales de 13 haplotipos únicos y 7 haplotipos de menor frecuencia (Fig. 4 a). El haplotipo H14 se presentó en 151 (45.8 %) de los individuos del total, y el segundo haplotipo H1 en 106 (32.1 %) chitones analizados. La genealogía de los haplotipos con el gen ITS 1 mostró al haplotipo dominante H4 conectado por un paso mutacional a los haplotipos H2, H3 y H5. El Haplotipo H2 conectado por un paso mutacional con el haplotipo H1. El haplotipo H4 fue encontrado en 74 individuos (49.3%), H4 en 25 (16.7%), H2 en 24 (16%) y H1 en 16 (10.7%) de los 150 chitones totales.

1.2.- Estructura genética: 1.2. a Método utilizado: la estructuración genética se calculó con el índice st (calculado con el n° de diferencia entre pares de secuencias) o Fst (calculado con el n° de diferencia en frecuencia haplotípica entre pares de secuencias desde distintas localidades), implementado en el programa Arlequin 3.1 (Excoffier et al. 2005). Todas las hipótesis comparativas fueron evaluadas con corrección de Bonferroni.

Especie: Tegula atra, Resultados: Los valores de st (con el gen COI) entre pares de localidades, mostraron una diferenciación genética significativa entre las localidades del norte con las localidades desde aproximadamente los 30° Sur hacia el sur (Tabla 5a). Con el gen ITS 1 no se encontraron diferencias significativas entre pares de localidades, excepto uno que no tuvo un valor alto de st (Tabla 5b). Especie: Siphonaria lessonii, Resultados: La matriz de st entre pares de localidades mostró que las diferencias genéticas con el gen COI aumentaron significativamente entre localidades más distantes geográficamente. Particularmente la clina genética comenzó desde Metri y aumentaron las diferencias entre localidades del extremo norte con las del extremo sur de Chile (Tabla 6 a). El gen nuclear ITS 1, los resultados de st mostraron también diferencia significativa mayormente entre las localidades hacia los extremo del rango de distribución (Tabla 6 b). Se está finalizando un manuscrito para enviar pronto a publicar sobre: Molecular footprint and niche modelling reveal the effects of the ice ages on the pulmonate limpet Siphonaria lessonii (Gastropoda: Siphonariidae) evolutionary lineages in Southeastern Pacific Ocean.

Especie: Fissurella picta, Resultados: Los valores de st (con el gen COI) entre pares de localidades no mostraron una diferenciación genética significativa (Tabla 7 a). Con el gen ITS 1 tampoco se encontraron diferencias genéticas significativas entre pares de localidades (Tabla 7 b). Especie: Chiton magnificus, Resultados: Los valores de st (con el gen COI) mostraron que las diferencias genéticas significativas entre pares de localidades siendo mayor entre localidades más distantes geográficamente. La matriz de st muestra una clina genética (Tabla 8 a) que comienza desde el centro y aumenta hacia el sur de Chile. Los valores de st (con el gen ITS1) también mostró diferencia genéticas significativas que aumentó entre localidades más distantes (Tabla 8 b) ______1.2. b Método utilizado: Mediante inferencia bayesiana se estimó la probabilidad de detectar estructuración genética poblacional con su respectiva organización espacial a través del uso de datos los datos de Dna mitocondrial y un modelo espacial geo-referenciado, implementado en el programa Geneland 4.0.2 (Guillot et al. 2008).

Especie: Tegula atra, Resultados: El análisis de estructuración genética (con el gen COI) con inferencia bayesiana (con el modelo de frecuencias no-correlacionadas, elegido mediante Bayes Factor= 0.079), evidenció una máxima estimación a posterior de dos cluster o poblaciones georreferenciadas de T. atra a lo largo del Pacifico sureste (Fig. 5). Especie: Siphonaria lessonii, Resultados: El análisis de estructuración genética (con el gen COI) con inferencia bayesiana (con el modelo de frecuencias no-correlacionadas, elegido mediante Bayes Factor= 2.78), evidenció una máxima probabilidad a posterior de 2 cluster poblaciones georreferenciadas de S. lessonii a lo largo del Pacifico sureste (Fig. 6). Especie: Fissurella picta, Resultados: El análisis de estructuración genética (con el gen COI) con inferencia bayesiana (con el modelo de frecuencias no-correlacionadas, elegido mediante Bayes Factor= 2.61) evidenció 4 cluster o poblaciones georreferenciadas de F. picta pero con valores muy bajos y no significativos de probabilidad a posterior a lo largo del Pacifico sureste (Fig. 7). Especie: Chiton magnificus, Resultados: Los análisis bayesianos (con el modelo de frecuencias correlacionadas, elegido mediante Bayes Factor= 27.934), evidenciaron una máxima probabilidad a posterior de ocurrencia de dos cluster o poblaciones georreferenciadas (Fig. 8). La posición espacial de la discontinuidad genética sugiere una población del norte con un límite entre Pichilemu-Duao, y otra población del sur. La asignación espacial de cada secuencia de pertenencia a cada población georreferenciada muestra valores de probabilidad a posterior > a 0,9) (Fig. 8). ______1.2. c Método utilizado: Se realizaron análisis de correlación de Mantel para evaluar aislamiento genético mediado por distancia geográfica con el gen COI solamente. Se realizó un test de mantel (Mantel 1967), se utilizó el modelo unidimensional, para grandes distancias (hábitats extendidos), que es la relación lineal entre Fst/1-fst con la distancia geográfica (km) (Rousset 1997), implementado en Arlequin 3.11; y el nivel de significancia del coeficiente de correlación fue calculado mediante 10000 permutaciones paramétricas dentro de las matrices de distancia.

Especie: Tegula atra, Resultados: El test de Mantel con un 66 % de correlación entre la matriz de distancia genética y de distancia geográfica, que fue significativa (P= 0.006, Fig. 9). Por lo tanto, los individuos T atra de localidades cercanas mostraron mayor asociación genética. Así, las localidades de los extremos del rango de distribución por el Pacifico sureste mostraron mayor distancia genética (Fig. 9). Especie: Sifonaria lessonii, Resultados: El test de Mantel mostro una asociación significativa entre la matriz de distancia genética y la de distancia geográfica (Fig. 10). Esto significa que existe un aislamiento por distancia, donde las sifonarias de localidades vecina mostraron menor diferencia genética entre ellas, y entre sifonarias de localidades más distantes, de los extremos del rango de distribución en el Pacifico sureste, mostraron las mayores diferencias genéticas (Fig. 10). Especie: Fissurella picta, Resultados: El test de Mantel con un -6 % de correlación entre la matriz de distancia genética y de distancia geográfica, y no significativo (P= 0.566, Fig. 11). Así, los individuos F. picta de las localidades hacia los extremos del rango de distribución por el Pacifico sureste no mostraron mayor distancia genética que entre localidades cercanas (Fig. 11). Especie: Chiton magnificus, Resultados: El test de Mantel dio un 59 % de correlación significativa (P= 0.0008) entre la matriz de distancia genética y de distancia geográfica. Esto significa que los chitones de localidades cercanas tenían una mayor relación genética, y a medida que aumentaba la distancia entre localidades fue mayor la diferencia genéticas entre ellas (Fig. 12). ______2.- Objetivo: Inferir la historia demográfica de las poblaciones de cada especie, a partir de los patrones de diversidad genética obtenidos con distintos marcadores moleculares.

2.1 Inferencia demográfica

2.1. a Método utilizado: Se realizaron inferencia demográfica de las poblaciones usando el métodos Bayesian Skyride (BSR, Minin et al. 2008) en el software BEAST version 1.7.2 (Drummond & Rambaut, 2007). Este análisis van muestreando, en el campo aleatorio Markov Gausiano en el Skyride, para generar una distribución a posteriori del parámetro Tamaño Efectivo Poblacional (Ne) a través del tiempo, y estimar el tiempo del ancestro común más reciente (tMRCA) en las genealogías de coalescencia. Especie: Tegula atra, Resultados: El análisis Bayesian Skyride (elegido con Bayes Factor, BF = 3.541), realizado con reloj molecular Log Normal (BF = 3,186), y con el modelo HKY de sustitución nucleotídica (-lnL = 757.7646, BIC = 4002.5659, programa jModelTest) y sin un modelo paramétrico pre-especificado, fue el utilizado para hacer la inferencia demográfica de las poblaciónes de Tegula atra. Este análisis reveló que la población I (Fig. 13) habría comenzado un proceso de expansión hace aproximadamente 20000 años atrás, incrementando en aproximadamente dos órdenes de magnitud en crecimiento del tamaño efectivo. Mientras, la población II habría comenzado el proceso de expansión hace 50000 años atrás, también incrementando en dos órdenes de magnitud en tamaño efectivo (Fig. 13). El tiempo estimado del ancestro común más reciente (MRCAt) de la población 1 sería de  22000 años, mientras que la población2 tendría un ancestro común más reciente de 50000 años, señalando una población más antigua, que la anterior (Fig. 13). Especie: Siphonaria lessonii, Resultados: El análisis Bayesian Skyride (elegido con Bayes Factor, BF = 3.413), realizado con reloj molecular Log Normal (BF = 3,657), y con el modelo HKY de sustitución nucleotídica (-lnL = 657.9346, BIC = 3195.5687, programa jModelTest) y sin un modelo paramétrico pre- especificado, fue el utilizado para hacer la inferencia demográfica de las poblaciones de Siphonaria lessonii. Este análisis reveló que la población 1 hacia el norte habría comenzado un proceso de expansión hace aproximadamente 200000 años atrás, incrementando en aproximadamente un orden de magnitud en tamaño efectivo (Fig. 14). La población 2 del sur también mostraría un comienzo de expansión de hace aproximadamente 200.000 años (Fig. 14). El tiempo estimado del ancestro común más reciente (TMRCA) señalaría que la población 1 sería más ancestral que la 2, con un tiempo estimado mucho mayor del ancestro común más reciente de aproximadamente 410.000 años y la población 2 del sur de aproximadamente 230.000 años (Fig. 14). Se está finalizando un manuscrito para enviar pronto a publicar sobre: Molecular footprint and niche modelling reveal the effects of the ice ages on the pulmonate limpet Siphonaria lessonii (Gastropoda: Siphonariidae) evolutionary lineages in Southeastern Pacific Ocean.

Especie: Fissurella picta, Resultados: El análisis Bayesian Skyride (elegido con Bayes Factor, BF = 2.43), realizado con reloj molecular Log Normal (BF = 3,51), y con el modelo JC de sustitución nucleotídica (-lnL = 663.345, BIC = 4953.687, programa jModelTest) y sin un modelo paramétrico pre-especificado, fue el utilizado para hacer la inferencia demográfica de la única población Fissurella picta por el sureste del Pacifico. El análisis reveló que la población habría comenzado un proceso de expansión hace aproximadamente 250000 años atrás, incrementando en casi dos órdenes de magnitud en tamaño efectivo (Fig. 15). El tiempo estimado del ancestro común más reciente (TMRCA) además señala una población muy ancestral, con un tiempo estimado de aproximadamente 310.000 años (Fig. 15). Especie: Chiton magnificus, Resultados: El análisis Bayesian Skyride (elegido con Bayes Factor, BF = 3.643), realizado con reloj molecular Log Normal (BF = 3,034), y con el modelo HKY de sustitución nucleotídica (-lnL = 677.7346, BIC = 4105.5557, programa jModelTest) y sin un modelo paramétrico pre- especificado, fue el utilizado para hacer la inferencia demográfica de las poblaciones de Chiton magnificus. Este análisis reveló que la población 1 habría comenzado un proceso de expansión hace aproximadamente 70000 años atrás, incrementando en un orden de magnitud en tamaño efectivo (Fig. 16). Mientras, la población 2 mostraría un comienzo de expansión en su tamaño efectivo hace aproximadamente 20000 años (Fig. 16). El tiempo estimado del ancestro común más reciente (TMRCA) también señala que la población 1 del norte sería una población más ancestral con una edad estimada de  77.653 años, mientras la población 2 del sur sería relativamente reciente con un ancestro de  21.322 años (Fig. 16). ______3.- Objetivo: Realizar paleo-reconstrucciones del rango de distribución de las especies en el pasado.

3.1 Reconstrucción del actual nicho del rango distribucional 3. 1. a Resultados de las cuatro especies: Se realizó un mapeo de puntos geo-referenciados de alrededor de 20 localidades a los largo del Pacífico sureste, donde se han realizado los muestreos de aproximadamente 40 individuos por cada especies. De las cincos especies, cuatro se realizaron para hacer el estudio de este objetivo que fueron: Tegula atra, Chiton magnificus, Sifonaria lessoni y Fissurella picta, durante los años 2014, 2015 y 2016 (Fig. 17). Asociado a la recolección de muestra del proyecto por la zona norte de Perú, se hizo un catastro exhaustivo sobre la diversidad de poliplacóforos (chitones) en la zona sur de la Provincia Panameña (Ecuador-Perú). Se recolectaron quitones de cinco localidades y se estimó su diversidad y similitud de los ensambles. De esta recolección se hizo publicación ISI titulada: Diversidad de poliplacóforos tropicales del sur de la Provincia Panameña. Lat. Am. J. Aquat. Res., 44(4): 807-814, 2016.DOI: 10.3856/vol44-issue4-fulltext-16. Asociado con los datos morfológicos realizados a las especies trabajadas en este proyecto. Se ha enviado una publicación asociado a los datos morfológicos de las especies Chitones y Fissurellas, y que es titulado: Mollusc Freaks: New teratological data on marine molluscs from the South Pacific Ocean, enviado a la Revista Marine Biodiversity: Manuscript Number: MARB-D-16-00329 (ver en anexos).

3.1 Paleoreconstrucción del nicho distribucional

3.1. b Método utilizado Datos paleo climáticos del Last Glacial Maximum (LGM; 21 ka), 10 variables océano bioclimáticas mencionadas anteriormente en un 5’ de resolución espacial que fueron generados desde Paleoclimatic Marine Spatial Ecology (Paleo-MARSPEC) dataset (Sbrocco et al 2014)

Especie: Siphonaria lessonii, Resultados: Después de la selección de variables no correlacionadas modernas para la calibración del modelo actual, de cuatro variables bioclimáticas oceanicas, finalmente se han seleccionado para el proceso de modelamiento que represente el máximo mensual de salinidad superficial del mar (SSS), la variación anual de SSS, la temperatura media anual de la temperatura superficial del mar (TSM) y la varianza anual de TSM. Un valor AICc apropiado fue seleccionado y el modelo final y una proyección hacia escenario pasado fueron realizados. Los modelos de nicho mostraron una moderada y homogénea en los valores ambientales de alta idoneidad en S. lessonii para casi todas las zonas costeras de Chile y Perú en el escenario actual (Fig. 18 a), similar a la distribución actual de esta especie descrita en este estudio. Mientras, la distribución pasada en el potencial LGM predicho, usando capas paleo climáticas exhibió una transición notable desde valores ambientales inadecuados a los adecuados dentro de las regiones chilenas del centro-sur (El Maule y regiones Bio-Bio, Fig. 18 a). Este escenario fragmentado se parecía a la ruptura genética detectada dentro de las poblaciones de S. lessonii en el sur de Chile en la zona de la Patagonia (Fig 18 b), donde hay evidencia de la formación de una capa de hielo durante el LGM se (Hulton et al., 2002).

Especie: Tegula atra, Fisurella picta, Chiton magnificus Resultados: Están en proceso de análisis ______4.- Objetivo: Evaluar si las poblaciones de las cinco especies de moluscos (con larva lecitotrófica y planctotrófica) que habitan en la zona norte-centro y sur de Chile, muestran un patrón distinto de diversidad y estructuración genética.

Resultados: Entre las cuatro especies que se realizaron los análisis se pudo determinar que: comparativamente el caracol Tegula atra (con larva lecitotrófica) fue el que mostró valores más bajo de diversidad genética (tabla 1a y Fig. 1a) con el gen mitocondrial COI, seguido después por el Chiton magnificus (con larva lecitotrófica) (tabla 4 y Fig. 4, a y b), luego Fissurella picta (Tabla 3 a y Fig. 3a) y finalmente con una muy alta diversidad la lapa Siphonaria lessonii (con larva planctotrófica) (tabla 2a y Fig. 2a). Con el gen ITS1 fue distinto, S. lessoni mostró menor diversidad genética, seguido por F. picta, luego C. magnificus y T. atra ( Tabla 1b, 2b, 3b y 4b). Las cuatro especies mostraron variación latitudinal en su diversidad genética (aunque algunas nos significativas), donde T. atra aumentó hacia el sur, S. lessonii y C. magnificus redujeron hacia el sur y F. picta no mostró variación. Tres especies, T. atra. S. lessonii y C. magnificus mostraron diferencia genética significativa entre pares de localidades y particularmente entre las de zona norte y con el sur (Tablas 5a, 6a y 8a), mientras F. picta no mostró diferencias significativas. Sin embargo, S. lessonii (con larva planctotrófica) mostró valores de st relativamente bajos (Tabla 3) comparado con las otras especies, sugiriendo una diferencia muy baja entre localidades. Por otra parte F. picta (con larva lecitotrófica) no mostró diferencia genéticas entre pares de localidades, lo que no se espera para una especie con menor potencial de dispersión que una especie larva planctotrófica. Las especies T atra, C. magnificus (ambas con larva lecitotrófica) y S. lessonii (con larva planctotrófica) mostraron un patrón de estructuración genética constituido por dos unidades poblacionales con valores significativos de probabilidad a posterior (Fig. 5, 6 y 8), mientras que F. picta mostró una estructuración genética de 4 unidades poblacionales que no fueron significativos, por los bajos valores de probabilidad a posterior (Fig. 7). Además, las especies T. atra, S. lessonii y C. magnificus mostraron un patrón de aislamiento genético por distancia geográfica (Fig. 9, 10 y 12), contrario a la ausencia de este patrón en F. picta (Fig. 11), lo que significó una homogeneidad genética en todo su rango de distribución (Fig. 3 a y b). ______5.- Objetivo: Evaluar si las poblaciones de las cinco especies de moluscos que habitan en la zona norte-centro y sur de Chile, muestran un patrón demográfico y de paleo-distribución distinto.

Resultados: Las inferencias demográficas con el gen COI realizadas en las especies Tegula atra y Chiton magnificus, mostraron un patrón una expansión demográfica semejante en sus poblaciones del norte (Fig. 13 y Fig. 16) y del sur (Fig. 13 y Fig. 16), diferenciando en el tiempo en que comenzó cada expansión. Ambas especies revelaron una expansión demografica en sus poblaciones; siendo las más recientes expansiones en las dos poblaciones de T. atra y la población del sur de C magnificus (entre 20000 y 70000 años, Fig. 13 y 16), siguiendo con una expansión más antigua en la población del norte ( 70000 años atrás Fig. 16) de C. magnificus. Fissurella picta y Siphonaria lessonii mostraron los patrones de expansión poblacional más antiguo (entre 200000 a 250000 años atrás con gen COI, Fig. 14 y 15). En este contexto, estos resultados sugieren que las tres especies, con desarrollo larval diferente, muestran evidencia demográfica y genética de haber afrontado de manera diferente el último período glacial del Pleistoceno, aunque no por consecuencia de su desarrollo biológico diferencial. Resultados: Los resultados de paleodistribución muestran en S. lessonii un modelo del nicho pasado distinto al contemporáneo (Fig. 18 a y b), afectando de manera distinta la distribución geográfica de las poblaciones del norte y sur durante los eventos glaciales del Pleistoceno. Esto implicó una reducción significativa del nicho geográfico de la población del sur, hacia latitudes más cálidas, lo que tiene correspondencia con la menor y distinta diversidad genética en la población del sur, y con el proceso de expansión poblacional en el sur de su rango de distribución. En este contexto, ambos análisis muestran un patrón de inferencia demográfica y de paleo distribución diferentes entre la zona norte y sur del rango de distribución. Se está finalizando un manuscrito para enviar pronto a publicar sobre: Molecular footprint and niche modelling reveal the effects of the ice ages on the pulmonate limpet Siphonaria lessonii (Gastropoda: Siphonariidae) evolutionary lineages in Southeastern Pacific Ocean. DESTAQUE OTROS LOGROS DEL PROYECTO TALES COMO: • Estadías de investigación. • Actividades de difusión y/o extensión en la temática del proyecto. • Cualquier otro logro no contemplado en los ítem anteriores y que Ud. quiera destacar. La extensión máxima de esta sección es de 1 página (letra tamaño 10, Arial o Verdana).

< Formación Recursos Humanos

En el proyecto Fondecyt 3140610, en la etapa en curso correspondiente al 2016 se ha desarrollado la formación científica de recurso humano, teniendo a cargo a un estudiante de pregrado de la carrera de Biología Marina de la U. Andrés Bello: Felipe Torres quien realizó su práctica profesional trabajando en el proyecto el año pasado 2015. Su práctica profesional consistió en el trabajo de laboratorio donde el alumno aprendió protocolo de trabajo de identificación de especies (usando claves) y toma de mediciones y análisis morfológicos de todas la especies del proyecto. El realizó mediciones morfológicas a todos los ejemplares por localidad por cada especie. Adicionalmente todos los datos morfológicos recopilados fueron anotados en un cuaderno adicionando los datos de cada localidad, latitud, longitud. Luego generaron una base de datos en Excel incorporando todos los datos por cada ejemplar y por especie. Posterior a eso, les realizaron corte de tejido, el cual fue guardado para posterior análisis genético. El alumno debió aprender usar claves taxonómicas para identificar las especies. Además debió aprender a usar algunos programas estadísticos como el Past.exe para aprender a analizar los datos morfológicos (en anexos va una constancia de la práctica profesional realizado bajo el proyecto y mi tutela). Este Este años 2016, el estudiante Felipe Torres se encuentra realizando su tesis de magister bajo mi tutela, asociado a unos de los objeticos del proyecto, para pronto defender su tesis de tesis de magister en Biología Marina, que deberá presentar en diciembre de este año 2016, en la Universidad Andrés Bello (en anexos va una constancia del desarrollo de su tesis de magister asociado al proyecto Fondecyt 3140610 y mi tutela).

Atentamente Dra. María Cecilia Pardo Gandarillas.

INFORME DE EVALUACION DEL (DE LA) INVESTIGADOR(A) PATROCINANTE

NOMBRE: MARCO A. MÉNDEZ

Por medio de la presente informo del estado de avance del tercer y último año del proyecto FONDECYT Nº 3140610 titulado “Filogeografía comparada de moluscos litorales de la costa chilena: determinando la relevancia de las estrategias de desarrollo larval en su historia evolutiva”, cuyo investigador responsable es la Dra. M. Cecilia Pardo Gandarillas. Este proyecto fue adjudicado por un periodo de tres años (2014-2016). A mi juicio, en relación con las actividades desarrolladas en esta etapa 2016, se avanzó satisfactoriamente, teniendo resultados que pueden responder a los objetivos específicos del proyecto. Lo anterior se fundamenta en que las actividades de terreno, la obtención de muestras y sus datos de georreferenciación, y las actividades relacionadas al trabajo en laboratorio, con el posterior análisis de datos, se han desarrollado de acuerdo a lo programado. Cada objetivo del proyecto se ha logrado de manera satisfactoria, teniendo además una publicación ISI titulada: Diversidad de poliplacóforos tropicales del sur de la Provincia Panameña; y una publicación enviada : Mollusc Freaks: New teratological data on marine molluscs from the South Pacific Ocean, y un manuscrito a enviar pronto titulado: Molecular footprint and niche modelling reveal the effects of the ice ages on the pulmonate limpet Siphonaria lessonii (Gastropoda: Siphonariidae) evolutionary lineages in Southeastern Pacific Ocean. Adicionalmente destaco que el proyecto en su ejecución avanzó más allá de lo originalmente planteado en el plan de trabajo (carta Gantt), que inicialmente fue propuesto.

La Dra. Pardo ha realizado un esfuerzo de muestreo del material biológico de las cinco especies propuestas en el proyecto, que incluyen un total de aproximadamente 25 localidades toda la zona norte, centro y de hasta el extremo sur de Chile y Perú también. En el laboratorio se ha logrado obtener secuencias del gen mitocondrial citocromo oxidasa subunidad I (COI) de 471 individuos de la especie Tegula atra, 330 de Chiton magnificus, 546 de Siphonaria lessonii y 106 Fissurella picta (cambiada por F. nigra, por problema de obtención de muestras). Con respecto a Austrolittorina araucana aunque se obtuvieron muestras de 15 localidades, fue imposible la amplificación de los genes COI e ITS1. Las otras cuatro especies antes mencionadas se logró obtener también secuencias del gen nuclear ITS1 (internal transcripted spacer 1) de 170 individuos de la especie Siphonaria lessonii a lo largo de todo Chile y Perú, 89 individuos de Chiton magnificus, 94 individuos de Tegula atra y 88 de Fissurella picta. Adicionalmente, se realizaron mediciones morfométricas de tres de las cinco especies propuestas en el proyecto, a saber: un total de 1164 individuos de Siphonaria lessoni, 602 de Chiton magnificus y todas las especies de Fissurellas, incluyendo 152 de Fissurella picta, sumando 1196 lapas medidas de todos los muestreos realizados durante estos tres años de proyecto. También se avanzó en los análisis de modelamiento de nicho ambiental para la especie

1 Siphonaria lessoni, y se está desarrollando lo mismo para las otras especies. En relación con la participación a congresos nacionales e internacionales durante la etapa 2016. La Dra. Pardo ha participado en el 19th International Congress of UNITAS MALACOLOGICA, The World Congress of Malacology (WCM) 2016, realizado en Julio del 2016, en Penang, Malaysia. Además, presentó junto con su estudiante (Felipe Torres) otros dos trabajos realizados con datos del proyecto en la LIX Reunión Anual de la Sociedad de Biología de Chile, en conjunto con la X Reunión Anual de la Sociedad Chilena de Evolución y XXVII Reunión de Anual de la Sociedad de Botánica, realizado desde el 8 al 10 de noviembre, en Concón, Valparaíso, Chile.

Finalmente, es importante destacar que la Dra. Pardo, ha impulsado la formación científica de recursos humanos a través del desarrollo de su proyecto. Durante el curso de la etapa 2015 ha tenido a su cargo a un estudiante de la carrera de Biología Marina de la U. Andrés Bello desarrollando su práctica profesional realizando mediciones morfométricas y análisis moleculares; y después en el 2016 el mismo estudiante se encuentra realizando una tesis de Magister de Biología Marina de la U. Andrés Bello, titulada: Modelamiento de nicho, diversidad genética e historia demográfica de Fissurella picta (Gmelin, 1791) (, Vetigastropoda) en el Sureste del Pacífico.

De acuerdo a los antecedentes presentados en el informe de avance Nº 3 correspondiente al desarrollo de la etapa 2016 del Proyecto FONDECYT 3140610. Es mi convicción que la Dra. Pardo ha avanzado y finalizado satisfactoriamente en su proyecto de Postdoctorado.

Les saluda cordialmente, Dr. Marco Méndez Torres.

Firma Investigador(a) Patrocinante ______Fecha: 31 de Octubre 2016____

2 PRODUCTOS

ARTÍCULOS Para trabajos en Prensa/ Aceptados/Enviados adjunte copia de carta de aceptación o de recepción.

Nº : 1 Autor (a)(es/as) : Ibáñez, C; Sellanes J; Pardo-Gandarillas MC Nombre Completo de la Revista : Latin American Journal of aquatic Research Título (Idioma original) : Diversidad de poliplacóforos tropicales del sur de la Provincia Panameña Indexación : ISI ISSN : 0718-560X Año : 2016 Vol. : 44 Nº : 4 Páginas : 8 Estado de la publicación a la fecha : Publicada Otras Fuentes de financiamiento, si las hay :

Envía documento en papel : no Archivo(s) Asociado(s) al artículo : Ibanez_et_al_2016_Diversidad_de_Poliplacoforos.pdf

Nº : 2 Autor (a)(es/as) : TORRES, FI; IBÁÑEZ, CM; SANHUEZA, VE; PARDO-GANDARILLAS MC Nombre Completo de la Revista : Molluscan Research Título (Idioma original) : Mollusc Freaks: New teratological cases on marine molluscs from the South Pacific Ocean Indexación : ISI ISSN : 13235818 Año : Vol. : Nº : Páginas : Estado de la publicación a la fecha : Enviada Otras Fuentes de financiamiento, si las hay :

Envía documento en papel : no Archivo(s) Asociado(s) al artículo : TMOS_2017_0028_Teratology.pdf OTRAS PUBLICACIONES / PRODUCTOS

Sin información ingresada.

CONGRESOS

Nº : 1 Autor (a)(es/as) : Pardo-Gandarillas, C; Torres, F; Fabres, A; Ibáñez C. Título (Idioma original) : MOLECULAR AND MORPHOLOGICAL GEOGRAPHIC VARIATION OF SIPHONARIA LESSONI ALONG THE SOUTHEASTERN PACIFIC OCEAN Nombre del Congreso : 48th Annual Western Society of Malacologist Meeting País : ESTADOS UNIDOS DE AMERICA Ciudad : Fullerton, Los Angeles, California Fecha Inicio : 25/06/2015 Fecha Término : 28/06/2015 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : Resumen_USA.pdf

Nº : 2 Autor (a)(es/as) : Fabres, A; Sanhueza, V; Ibáñez, CM.; Méndez, MA.; Pardo-Gandarillas, MC. Título (Idioma original) : Influencia del desarrollo larval en la historia evolutiva reciente de Chiton magnificus Deshayes, 1827 Nombre del Congreso : IX Reunión Anual de la Sociedad Chilena de Evolución País : CHILE Ciudad : Talca Fecha Inicio : 13/10/2015 Fecha Término : 15/10/2015 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : Resumen_Fabre_et_al_Talca_Chile.pdf Nº : 3 Autor (a)(es/as) : Torres, FI.; Hernández, MF.; Varas, DG.; Ibáñez, CM.; Pardo-Gandarillas, MC. Título (Idioma original) : Alometría y morfología del género Fissurella (Mollusca: Vetigastropoda) del Pacífico Sureste Nombre del Congreso : IX Reunión Anual de la Sociedad Chilena de Evolución País : CHILE Ciudad : Talca Fecha Inicio : 13/10/2015 Fecha Término : 15/10/2015 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : Resumen_Torre_et_al_Talca_Chile.pdf

Nº : 4 Autor (a)(es/as) : Pardo-Gandarillas C; Campos C; Hernández F; Salazar J; Torres F; Varas D; Sepúlveda R & Ibáñez C. Título (Idioma original) : Variación latitudinal del tamaño corporal y morfología de Siphonaria lessoni a lo largo de la costa sureste del Pacifico Nombre del Congreso : LVII Reunión Anual de la Sociedad de Biología y de las Sociedades de Genética, Evolución, Ecología y Botánica País : CHILE Ciudad : Puerto Varas Fecha Inicio : 25/11/2014 Fecha Término : 28/11/2014 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado :

Nº : 5 Autor (a)(es/as) : Pardo-Gandarillas MC; Sanhueza, V; Torres, F; Méndez, MA; Ibáñez CM Título (Idioma original) : COMPARATIVE PHYLOGEOGRAPHY IN LITTORAL MOLLUSKS FROM SOUTHEASTERN PACIFIC OCEAN: THE IMPORTANCE OF LARVAL DEVELOPMENT STRATEGIES IN THEIR EVOLUTIONARY HISTORY Nombre del Congreso : World Congress of Malacology País : MALASIA Ciudad : Penang Fecha Inicio : 18/07/2016 Fecha Término : 24/07/2016 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : pardo.pdf

Nº : 6 Autor (a)(es/as) : Ibáñez, C; Eernisse, DJ; Méndez, MA; Fabres, A; Sanhueza, V; Sellanes, J; Pardo-Gandarillas MC Título (Idioma original) : Phylogeny, divergent time and delimitation of Tonicia (Mollusca: Polyplacophora) from the Eastern Pacific Ocean Nombre del Congreso : World Congress of Malacology País : MALASIA Ciudad : Penang Fecha Inicio : 18/07/2016 Fecha Término : 24/07/2016 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : ibanez.pdf

Nº : 7 Autor (a)(es/as) : Pardo-Gandarillas, MC; Sanhueza, V; Torres, F; Méndez, MA; Ibáñez Christian CM Título (Idioma original) : Filogeografía comparativa en moluscos costeros del Pacifico Sureste: relevancia del desarrollo larval en su historia evolutiva Nombre del Congreso : LIX REUNIÓN ANUAL 2016, X Reunión Anual Sociedad Chilena de Evolución País : CHILE Ciudad : Concón Fecha Inicio : 08/11/2016 Fecha Término : 10/11/2016 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : Pardo2016.pdf

Nº : 8 Autor (a)(es/as) : Torres, F I1; Ibáñez, CM; Pardo-Gandarillas, MC Título (Idioma original) : DIVERSIDAD GENÉTICA E HISTORIA DEMOGRAFICA Nombre del Congreso : LIX REUNIÓN ANUAL 2016, X Reunión Anual Sociedad Chilena de Evolución País : CHILE Ciudad : Concón Fecha Inicio : 08/11/2016 Fecha Término : 10/11/2016 Nombre Publicación : Año : Vol. : Nº : Páginas : Envía documento en papel : no Archivo Asociado : Torres2016.pdf

ANEXOS

Nº : 1 Archivo Asociado : InformeSeguimientoBioetico__130717_001_FCS_UCHs.pdf

Nº : 2 Archivo Asociado : Anexo_relevantes__Tablas__figurasInforme__2016.pdf

Nº : 3 Archivo Asociado : PractProfesioF.Torre_asociP3140610_.pdf

Nº : 4 Archivo Asociado : TesisMagisterAsociaP3140610.pdf

Nº : 5 Archivo Asociado : CERTIF._COMPT_._HORRARIA_MAR_A_CECILIA_PARDO_.pdf

Nº : 6 Archivo Asociado : Permiso_pesca_Pardo_Gandarillas.pdf

A continuación se detallan los anexos físicos/papel que no se incluyen en el informe en formato PDF.