RESEARCH ARTICLE        

Genetic diversity and phylogeny of the genus Cinchona in Cutervo National Park, Peru

Diversidad genética y filogenia del género Cinchona en el Parque Nacional de Cutervo, Perú

Karla Andrea Peña Yacila1, * ; Erick Antonio Suarez-Peña2 ; Carlos Alberto Torres Pera3;

Luis Alberto Bermejo Requena1 ; Luis Xavier Llacsa Sánchez4 ; Iris Zárate Rodríguez5;

Oscar Junior Paredes-Vilca6 ; José Yasser Dávila García6

1  Universidad Nacional de Tumbes, Laboratorio Microbiología, Av. Universitaria S/N Pampa Grande – Tumbes. Peru.

2 Instituto Nacional de Innovación Agraria (INIA), Estación Experimental Agraria Los Cedros, km 12-Panamericana Norte-Caserío los Cedro, Tumbes. Peru.

3 Universidad Nacional Mayor de San marcos, Facultad de Ciencias Biológicas, Av. Carlos Germán Amézaga 375 - Edificio Jorge Basadre, Ciudad Universitaria, Lima. Peru.

4 Inca’Biotec SAC- Perú, departamento de Sanidad Vegetal, Jr. Filipinas 212, Tumbes, Peru.

5 Servicio Nacional de Áreas Naturales Protegidas por el Estado-SERNANP, Av. San Juan 724, Cutervo-Cajamarca. Peru.

6 Universidad Nacional Agraria La Molina-UNALM, Av. La Molina s/n, La Molina. Peru.

* Corresponding author: k.p.yacila@gmail.com (K. A. Peña Yacila).

Received: 13 July 2021. Accepted: 18 November 2021. Published: 15 December 2021.

Abstract

Species of the genus Cinchona are recognized for their effectiveness in treating malaria. For Peru, this genus is of greater importance regarding forestry and civic aspects, because it represents the richness of the plant resource. Nevertheless, nowadays Cinchona populations are strongly threatened by various anthropogenic activities. Cutervo National Park (PNC) is the Peruvian's first Protected Natural Area (ANP) and is considered as a conservation unit for the genus Cinchona. Objective: Determine the genetic diversity and phylogenetic relationships of the species of the genus Cinchona that inhabit the PNC. Methods: We tested 13 samples using the matK marker and a phylogenetic tree was generated using the maximum likelihood method and bootstrap measures of phylogenetic support. Results: 8 number of haplotypes (h) were recorded, with haplotype diversity (Hd) 0.859 and nucleotide diversity (π) 0.35%, also identifying: C. calisaya, C. pubescens and C. pitayensis for the PNC; the phylogenetic tree showed C. pitayensis as conserved specie, meanwhile C. pubsecens is present in many groups. Conclusions: The high genetic diversity of Cinchona and haplotype record of C. calisaya, species with a national conservation category of "Vulnerable", increases the importance for biodiversity conservation of the PNC.

Keywords: Conservation; Cinchona tree; haplotype diversity; nucleotide diversity; matK marker.

Resumen

Las especies del género Cinchona son reconocidas mundialmente por su eficacia en el tratamiento contra la malaria. Para el Perú, este género reviste mayor importancia a los aspectos forestal y cívico, por representar la riqueza del recurso vegetal. Sin embargo, en la actualidad sus poblaciones se encuentran fuertemente amenazadas por diversas actividades antropogénicas. El Parque Nacional de Cutervo (PNC) es la primer Área Natural Protegida (ANP) del Perú y es considerado como unidad de conservación para el género Cinchona. Objetivo: Determinar la diversidad genética y las relaciones filogenéticas de las especies del género Cinchona presentes en el PNC. Métodos: Nosotros evaluamos 13 muestras empleando el marcador matK y se generó un árbol filogenético empleando el método de máxima verosimilitud con soporte de filogenia bootstrap. Resultados: 8 haplotipos fueron registrados, con una diversidad haplotípica (Hd) 0.859 y nucleotídica (π) 0.35 %, identificándose además a: C. calisaya, C. pubescens y C. pitayensis para el PNC; el árbol filogenético mostró que la especie C. pitayensis se encuentra bastante conservada, mientras que C. pubsecens se encuentra presente en distintos grupos con similitud variable. Conclusiones: La elevada diversidad genética de Cinchona y el registro del haplotipo de C. calisaya, especie con categoría de conservación nacional “Vulnerable”, incrementa la importancia para la conservación de la biodiversidad del PNC.

Palabras clave: Conservación; Árbol de la quina; diversidad haplotípica; diversidad nucleotídica; marcador matK.

DOI: https://dx.doi.org/10.17268/sci.agropecu.2021.056

Cite this article:

Peña Yacila, K. A., Suarez-Peña, E. A., Torres Pera, C. A., Bermejo Requena, L. A., Llacsa Sánchez, L. X., Zárate Rodríguez, I., Paredes-Vilca, O. J., & Dávila García, J. Y. (2021). Diversidad genética y filogenia del género Cinchona en el Parque Nacional de Cutervo, Perú. Scientia Agropecuaria, 12(4), 517-524.

1. Introducción

La biodiversidad se comprende en diferentes niveles (la variedad de ecosistemas, especies y genes) y la compleji­dad de sus interacciones conforma el sistema de soporte de la vida (De Andrade, 2013; Rawat & Agarwal, 2015), sus muchos beneficios han sido confirmados en múltiples es­tudios y aún en la actualidad continúan revelándose, lle­gando incluso a proponerse a la biodiversidad como con­trolador en la manifestación de enfermedades emergen­tes por zoonosis (Ezenwa et al., 2006; Johnson & Thieltges, 2010; Keesing et al., 2010).

No obstante, a pesar de todos sus beneficios, se registra una alarmante pérdida de la biodiversidad, ligado a un marcado efecto antropogénico sobre todos los niveles de variación biológica (Hoffmann et al., 2010), reportándose tasas de extinción entre 100 y 1000 veces superiores a las tasas de los cincos episodios de extinción masiva anterio­res a la aparición del ser humano (Pimm et al., 1995). En este contexto, es una prioridad global acelerar la descrip­ción de la biodiversidad y gestionar su conservación (Marchese, 2015; Martins et al., 2020).

La diversidad de especies y ecosistemas deriva de proce­sos evolutivos a nivel genético, por ello, esta es la base de la biodiversidad y un requisito primordial para la evolución (Collada & Jiménez, 2008; Piñero, 2008). Si la diversidad genética de una población se reduce, indefectiblemente disminuye su capacidad de adaptarse ante potenciales cambios ambientales, por lo que, mantener la variabilidad genética es vital para responder exitosamente a las pre­siones ambientales de selección.

El estudio de la diversidad genética se basa en el análisis de potentes marcadores moleculares de naturaleza pro­teica y de ADN (nuclear y citoplasmático), en particular para el grupo de las plantas, los marcadores del cloro­plasto (ADN cp), destacan en la investigación de la historia de las poblaciones, flujo genético, hibridación e identifica­ción de áreas con elevada diversidad genética (Caruso et al., 2015; Collada & Jiménez, 2008; Salvador et al., 2000). El grupo de trabajo de plantas del “Consortium for the Barcode of Life” (CBOL), recomienda el uso de los genes cloroplastidícos matK + rbcL, para desarrollar códigos de barras (Little, 2014), también se ha reportado que el marcador matK tiene una tasa evolutiva rápida y provee una alta señal de filogenia para resolver la sistemática en plantas (Hilu & Liang, 1997).

El género Cinchona está representado por 24 especies botánicas, comúnmente conocidas como “quina” o “cas­carilla”, de importancia en el tratamiento contra la malaria por su composición de alcaloides de quinina en la corteza (Aymard, 2019; Cueva-Agila et al., 2019; Zevallos, 1989; Nair, 2010). Las especies de este género, usualmente, ha­bitan en los bosques de montaña de los sectores central y noreste de cordillera andina (Antonelli et al., 2009; Aymard, 2019).

El Perú registra 18 de las 24 especies reconocidas, siendo el país con mayor cantidad de especies para el género Cinchona (Albán, 2013; Aymard, 2019; Huamán, Albán, & Chilquillo, 2019). Este género reviste importancia para el país, más allá de los niveles forestal y cívico, por representar la riqueza del recurso vegetal (Huamán et al., 2019; Zevallos, 1989). A pesar de esto, actualmente sus poblaciones se encuentran fuertemente amenazadas por actividades, como la sobre explotación del recurso, deforestación y degradación del hábitat por crecimiento urbano (Huamán et al., 2019; López, 2016; Mendoza, Ramírez, & Jiménez, 2014).

A nivel internacional diversos estudios moleculares en el género Cinchona y taxones relacionados emplearon el marcador matK, obteniendo información representativa en sistemática y filogenia (Andersson & Antonelli, 2005; Chilquillo, 2016; Maldonado et al., 2017). En contraparte, para el Perú, si bien los datos biológicos existentes sobre el género Cinchona registran aspectos históricos, taxonó­micos, farmacológicos y ecológicos (Albán, 2013; Huamán et al., 2019; Nair, 2010; Zevallos, 1989), se carece de estu­dios moleculares para este género y del conocimiento de su historia evolutiva, datos fundamentales para el correcto establecimiento de actividades de restauración sin perjuicio de la diversidad local (Albán-Castillo et al., 2020).

El PNC, establecido el 20 de setiembre de 1961, es la pri­mer ANP de esta categoría en el Perú y se enfoca en la conservación del ave “guácharo” Steatornis caripensis y ecosistemas de tipo páramo, bosque seco y bosque hú­medo (SERNANP, 2014; SERNANP 2017). Adicionalmente, es considerado una de las cuatro unidades de conserva­ción para el género Cinchona en el Perú (Zevallos, 1989), pese a ello, hasta la actualidad se carece de investigacio­nes para este género en el PNC. Por consiguiente, el objetivo del presente estudio fue determinar la diversidad genética y las relaciones filogenéticas de las especies del género Cinchona que habitan en el PNC, con la finalidad de generar información que permita la conservación de este recurso.

2. Materiales y métodos

Área de estudio y recolección de la muestra: Corresponde al PNC, el cual se encuentra localizado en el área deno­minada Cordillera de Tarros (2200-3500 m.s.n.m.), se ubica en los andes del norte del Perú, protegiendo bos­ques montanos, páramos, y principalmente los bosques de neblina de las nacientes más altas de la cuenca del Marañón en la Región Cajamarca.

El material botánico fue colectado durante los meses de febrero a marzo del 2016, en siete centros poblados (Figura. 1), distribuidos en los distritos de: San Andrés de Cutervo, Santo Tomás, Santo Domingo de la Capilla, ubi­cados en el departamento de Cajamarca, Perú (Tabla 1).

La identificación preliminar en campo de los individuos del género Cinchona, se basó en caracteres morfológicos, si­guiendo las claves establecidas Zeballos (1989) y biblio­grafía científica especializada (Aymard, 2019; Gómez, Beraun, Gómez, & Llatas, 2016; Killeen, García & Beck, 1993; Mostacero, Mejía & Gamarra, 2009). Con la finalidad de evaluar hojas sanas, se colectaron dos muestras de ho­jas maduras por individuo, extraídas de la parte superior de la planta, a una altura aproximada de 2,5m, solo se procesó una hoja por individuo, la segunda hoja fue al­macenada como material de respaldo. Durante la colecta se registraron las coordenadas geográficas, el material biológico fue rotulado y transportado en bolsas herméticas estériles al laboratorio de la empresa Inca´Biotec SAC para su posterior procesamiento.

Tabla 1

Coordenadas de muestreo en el Parque Nacional Cutervo, Cajamarca-Perú

Datum WGS84, Zona 17S.

Figura 1. Localización de los puntos de muestreo.

Extracción de ADN cloroplastidial (ADN cp): Se realizó la extracción de ADN cp, siguiendo el protocolo establecido por Doyle & Doyle (1987), adaptado para el género Cinchona, mediante ensayos preliminares en el laboratorio. La superficie de las hojas fue desinfectada con alcohol al 70%, consecutivamente, se extrajo un corte de 2 cm2, evitando las nervaduras. Las muestras fueron trituradas en un mortero, con nitrógeno líquido, transfiriéndose 0,1g del triturado a un microtubo, en el cual se agregaron 590 µl del buffer de extracción Bromuro de Cetiltrimetilamonio (CTAB 2X) previamente calentado a 60 °C, y 10 µl de b-Mercaptoetanol. Seguidamente, la mezcla fue homogeneizada en Vortex (GB Lab®) durante 1 minuto, mantenida a una temperatura de 60°C por 60 minutos (uniformizando cada 15minutos), luego se adicionaron 300 µl de fenol concentrado y 300 µl de Cloroformo/Alcohol isoamílico (24:1).

La mezcla fue homogeneizada y centrifugada a 12 mil rpm por 10min a 4 °C, posteriormente se transfirió el sobrenadante a un nuevo microtubo, evitando tomar la interfaz, y se adicionó un volumen igual de Cloroformo /Alcohol Isoamílico (24:1), para ser centrifugada secuencialmente a 13 mil rpm durante 10 min a 4 °C y a 12 mil rpm por 10 minutos a 4 °C, recuperándose el sobrenadante a un nuevo microtubo. Rápidamente, se agregó Isopropanol helado (6/10) y se incubó a -20 °C por un lapso de 30-60 minutos. Se centrifugó secuencialmente la muestra a 10 mil rpm a 4 °C por 10 y 5 minutos, eliminando el sobrenadante y agregando 500 µl de etanol al 75% después del primer centrifugado, posteriormente se eliminó el sobrenadante.

Finalmente, el pellet se dejó secar a temperatura ambiente, se suspendió en 50µl de T.E (Tris HCL 10mM y EDTA 1Mm pH 8.0) y se agregó 1 µl de RNasa (10 mg/ml) e incubó a 37 °C por 30 minutos, el ADN se conservó a -20 °C. La cuantificación del ADN se realizó por espectrofotometría, con un volumen total de 60 μl de la muestra y una dilución de 5 μl ADN + 55 μl buffer TE pH: 8.

Amplificación del gen actiplant: Se amplificó la región específica del gen de actina del ADNcp mediante PCR, empleándose los siguientes primers: Forward (5’-CAC TAC TAC TGC TAA ACG GG AAA-3’) y Reverse (5’-ACA TCT GCT GGA AGG TGC TG-3’); y los volúmenes de reacción: 2,5 µL de Buffer Premix 10X, 0,5 µL de dNTP’s, 1 µL de MgCl₂ a 25 mM, 1 µL de Primer ACTplant F, 1 µL de Primer ACTplant R, 0,2µL de Taq DNA Polymerase (PlatinumR), 16,8µL de agua libre de nucleasas y 2 µL de ADNcp. Los ciclos térmicos para la amplificación consistieron en una desnaturalización inicial a 95 °C por 5 minutos, seguido por 35 ciclos compuestos por una desnaturalización a 95 °C por 5 segundos, hibridación a 54 °C por 45 segundos, elongación a 72 °C por 45 segundos, concluyendo con una extensión final a 72 °C por 6 minutos.

Amplificación del gen matK: La amplificación de la región del gen matK del ADNcp se realizó con los primers Forward (5’-CCT TCG ATA CTG GGT RAA AGA T-3’) Y Reverse (5’-CGT ACA GTA YTY TTG TGT TTA CGA-3’), para lo cual se emplearon los siguientes volúmenes de rección: 2,50µl Buffer Premix 10X, 0,50 µL de dNTPs (10mM), 0,50 µL de BSA (100X), Taq Polimerasa (0,5 U / μl) 0,10 µL, Primer Reverse matK3B (15 pMol / μl) 0,60 µL, Primer Forward matK2A (15 pMol / μl) 0,60 µL, agua ultra pura 18,20 µL, ADN 2,00 µL. Los ciclos térmicos para la amplificación consistieron en una desnaturalización inicial a 94 °C por 5 min, seguido por 34 ciclos compuestos por una desnaturalización 94 °C por 30 segundos, hibridación a 49 °C por 45 s; extensión 72 °C por 45seg, concluyendo con una extensión final a 72 °C durante 6min.

Validación de extracción de material genético: La validación de la extracción de ADN se llevó a cabo mediante la migración del ADN genómico en un gel de agarosa al 1%, agregándose 0,5 g de agarosa y 50 ml de TAE 1X; mientras que para los productos de los genes actiplant y matK se migraron geles de agarosa al 1,5%, agregándose 0,9 g de agarosa y 60 ml de TAE 1X. Ambos productos emplearon Blue/Orange Loading Dye 6X como tampón de depósito a una relación 5 : 1 μl por cada muestra. La migración del ADN se llevó a 90V por 15 min, colocando posteriormente el gel de agarosa en el transiluminador UV.

Secuenciamiento y análisis filogenético: Los productos de la PCR fueron preparados para posteriormente ser enviados a secuenciar a MACROGEN®. Las secuencias fueron editadas y alineadas con el algoritmo del software MAFFT (https://mafft.cbrc.jp/alignment/software/), en el software Mr Bayes 3,2 (Ronquist et al., 2012). Para la identificación de especies se consideró un límite máximo de divergencia genética de 1%.

Las secuencias del gen matK obtenidas por PCR fueron comparadas con las secuencias disponibles en la base de datos del GenBank usando el Basic Local Alignment Search Tool (BLAST) de la National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/BLAST)

El análisis filogenético se efectúo empleando el método de máxima verosimilitud (ML), la prueba de Filogenia Bootstrap se efectuó con 1000 repeticiones. La raíz del del filograma se infirió empleando el haplotipo generado de la secuencia del gen matK para la especie Remijia pedunculata como grupo de referencia.

El análisis diversidad genética se realizó mediante el polimorfismo de secuencias de ADN (DnaSP), evaluando como indicadores el número de haplotipos (h), la diversidad haplotípica (Hd) y la diversidad nucleotídica (π), los valores de los indicadores se calcularon con el software DnaSP (Librado & Rozas, 2009).

3. Resultados y discusión

Identificación de secuencias mediante el gen matK: Por medio del gen matK en la base de datos del GenBanK, se logró identificar 4 especies con una alta similitud: C. pitayensi, C. pubescens, C. calisaya y Remijia pedunculata (99%), con un índice de cobertura (Query Cover) promedio de 100%, corroborando una identificación precisa para todos los individuos (Tabla 2).

Tabla 2

Diversidad genética de especies del género Cinchona mediante el gen matK (En todos los casos se registró un índice de cobertura del 100% y similitud al 99%)

Descripción de centros poblados: Tarros (TR); La Flor (LF1); La Flor (LF2); La Granadilla (LG1); Zona de Amortiguamiento (Puente El Suro-límite entre el distrito de San Andrés de Cutervo y Santo Tomás) (ZA); La; Granadilla (LG2); Playa Grande (PG1); Playa Grande (PG2); Las Grutas (LGR1); Las Grutas (LGR2); Las Grutas (LGR3); El Suro (ES1) y El Suro (ES2).

Análisis filogenético: El árbol filogenético muestra que las muestras identificadas como C. pitayensis se agrupan en un solo clado con un elevado valor de soporte de bootstrap (<98%); por otra parte, destaca que las mues­tras identificadas como C. pubescens se agrupan en dis­tintas ramas, algunas muestras de C. pubescens se agru­pan en un solo clado con un valor elevado de soporte de Bootstrap (100% y 90%), sin embargo, son disimiles entre grupos. Finalmente, otro gran grupo conformado por C. pubescens, C. calisaya y R. pedunculata se encuentran agrupados (bootstrap 79), mientras que C. calisaya, se encuentra en una politomía sin resolver (Figura 2).

Distribución de sitios polimórficos: Se registraron 45 sitios polimórficos para el grupo de secuencias, descritas como: 56, 59, 114, 122, 123, 180, 183, 199, 205, 216, 44, 45, 47, 58, 112, 123, 180, 183 ,199, 205, 216, 55, 124, 158, 182, 46, 138, 140, 152, 154, 156, 159, 160, 179, 200, 201, 210, 211, 212, 213, 214, 215, 217, 218, 219, 110, 153, 204.

Análisis de diversidad genética: Entre los 13 organismos analizados se registró un total de 8 haplotipos (h), estando el primer haplotipo conformado por las muestras 1, 8, 11, 30 y 32; el segundo haplotipo, por las muestras 17 y 18; mientras que el resto de los haplotipos son correspondientes de manera individual a las muestras 12, 14, 15, 21, 23 y 25 (Tabla 3); finalmente se comprobó que los resultados muestran un alto nivel de diversidad haplotípica (Hd = 0,859) y nucleotídica (π = 0,35 %) para las poblaciones de las especies del género Cinchona dentro del PNC.

El centro poblado de Las Grutas (LG) registró la mayor cantidad de haplotipos por localidad, correspondiente a la especie C. pubescens, por otra parte, el centro poblado de Playa Grande registro el único haplotipo de la especie C. calisaya con dos muestras analizadas; C. pitayensis se identificó en tres muestras de un mismo haplotipo para los centros poblados de Tarros y La Flor; de manera resaltante, la especie C. pubescens se identificó en siete muestras para todos los centros poblados con excepción de Tarro, La Flor y Playa Grande, sin embargo, se agrupó en cinco haplotipos, siendo las secuencias de las muestras del centro poblado El Suro (30 y 32) compartidas con el haplotipo de C. pitayensis.

Para el PNC se logró identificar a la especie R. pedunculata y tres especies del género Cinchona: C. calisaya, C. pubescens y C. pitayensis; representando aproximadamente el 17% de la riqueza de especies reconocidas para el género Cinchona en el Perú, la cual asciende a 18 especies (Albán, 2013; Aymard, 2019; Huamán et al., 2019). De acuerdo con el Plan Maestro del PNC y publicaciones científicas referentes, se reconoce para esta ANP la presencia de C. humboldtiana, C. offinalis y C. pubescens (Sánchez & Calderón, 2010; SERNANP, 2014), por lo que nuestros resultados representan el primer registro de las especies C. pitayensis y C. calisaya para el PNC.

Nuestros resultados muestran conflicto con la idea ampliamente extendida de que R. pedunculata es la especie hermana del clado monofilético conformado por el género Cinchona (Andersson & Antonelli, 2005; Andersson, 1995; Aymard, 2019; Chilquillo, 2016). Cabe mencionar que la presencia de caracteres mixto entre Ladenbergia y Cinchona, determinó la reclasificación de Remijia pedunculata dentro del nuevo género Ciliosemina (Andersson & Antonelli, 2005; Aymard, 2019; Chilquillo, 2016; Manns & Bremer, 2010).

Tabla 3

Distribución de haplotipos identificadas con el marcador matK

Figura 2. Árbol filogenético obtenido por Máxima verosimilitud de individuos del género Cinchona del PNC (haplotipos citoplasmáticos del gen matK). Los números en los nodos representan el valor de soporte de Bootstrap. Grupo hermano (“outgroup”): Remijia pendunculata.

El árbol filogenético reporta que las muestras identificadas como C. pitayensis se agrupan en un solo clado con un elevado valor de soporte de bootstrap (<98%); por otra parte, destaca que las muestras identificadas como C. pubescens se agrupan en distintas ramas, algunas de estas muestras se agrupan en un solo clado con un valor elevado de soporte de bootstrap (100% y 90%), sin embargo, estas mismas son disimiles entre grupos. Finalmente, otro gran grupo conformado por C. pubescens, C. calisaya y R. pedunculata se encuentran vinculados (bootstrap 79%), mientras que C. calisaya se encuentra en una politomía sin resolver (Figura 2).

La robusta agrupación de C. pitayensis demuestra que esta especie se encuentra conservada en el PNC. Por otra parte, la distribución de C. pubescens en diversos clados independientes y en conjunto con C. calisaya y R. pedunculata, indican que esta especie tendría variedades hibridando con otras especies en determinadas zonas y conservadas en otras.

Patrones similares a los reportados en el presente estudio fueron reportados en el árbol filogenético de la tribu Cinchoneae, generado través del método de máxima parsimonia empleando cinco marcadores (región ITS, gen matK, gen rbcL, intrón rps16 y la región trnL-F) (Andersson & Antonelli, 2005); y para el árbol filogenético generado para Cinchoneae e Isertieae a través del método de inferencia bayesiana empleando siete marcadores (ITS1, ITS2, la región espaciadora atpB-rbcL, ndhF, rbcL, el intrón rps 16 y la región trn-T-L-F) (Manns & Bremer, 2010); esto a pesar de que el uso de diferentes marcadores moleculares y métodos de evaluación pueden generar diferentes valores de diversidad.

El registro de una mayor cantidad de haplotipos de C. pubescens sería un indicador de la presencia de variedades originadas por la presión ambiental, la cual ocasiona el desarrollo de caracteres adaptativos con tendencia a la especiación (García, 2012). Lo anterior mencionado sería correspondiente con lo manifestado por Albán-Castillo et al. (2020), quienes mencionan que la alta plasticidad morfológica presente en varias especies de Cinchona puede ser atribuida a gradientes ambientales, así como a eventos de hibridación o duplicación de sus genomas.

La diversidad haplotípica (Hd = 0859) y nucleotídica (π =0,35%) del género Cinchona en el PNC es elevada respecto del área del PNC (8,214.23 hectáreas) (SERNANP, 2017). Se conoce que la diversidad genética está afectada por eventos demográficos históricos y contemporáneos que pueden ser atribuibles a un contacto secundario entre linajes alopátricos previamente diferenciados, o bien a una gran población estable con una larga historia evolutiva (Grant & Bowen, 1998). A nivel local, la distribución de la diversidad genética registró al centro poblado de “Las Grutas” con la mayor cantidad de haplotipos (C. pubescens) y a “Playa Grande” con un único haplotipo (C. calisaya); mientras que la especie C. pitayensis con un único haplotipo se registró en los centros poblados de “Tarros” y “La Flor” y C. pubescens del centro poblado “El Suro” (30 y 32), esto se relacionaría con el comportamiento diferencial de las poblaciones, debido a factores actuales como la agricultura y ganadería local; y factores históricos como la formación de la geografía presente.

Los resultados aquí expuestos respaldan la hipótesis de que la distribución biogeográfica de C. calisaya, al sur del Perú y en Bolivia, con un reciente registro puntual de ampliación hacia el norte, representaría un caso de vicarianza respecto de la especie C. pubescens (bootstrap 79%), de amplia distribución, tal como lo propusieron Andersson & Antonelli (2005). No obstante, para respaldar esta hipótesis se recomienda complementar los estudios con el uso de más marcadores moleculares.

Con respecto a C. calisaya, la cual se encuentra en categoría de conservación “Vulnerable” para la legislación nacional (MINAGRI, 2006) y representa endemismo para el Perú y Bolivia, esta es reconocida como la especie con mayor cantidad de alcaloides (Nair, 2010), mientras que, recientes investigaciones indican que su composición de alcaloides anti malaria a nivel local estaría siendo dirigida por su filogenia (Maldonado et al., 2017). En consecuencia, el haplotipo de C. calisaya descrito para el PNC podría ser una variante particular de su diversidad genética, por lo que continuar con estudios en esta especie es preponderante.

4. Conclusiones

Mediante este estudio, se incrementó el número de especies del género Cinchona en el PNC con la identificación molecular de C. calisaya, C. pubescens y C. pitayensis, asimismo se confirmó la presencia de R. pedunculata. La resolución de la topología del árbol filogenético confirma que C. pitayensis se encuentra conservada dentro del PNC, mientras que la identificación de 7 haplotipos distribuidas en tres grupos de Cinchona demuestra la elevada diversidad genética presente en el PNC. La especie C. pubsecens se encuentra presente en distintos grupos con similitud variable. Los valores de la diversidad haplotípica (Hd = 0,859) y nucleotídica (π = 0,35%) del género Cinchona y el registro del haplotipo de la especie C. calisaya, endémica para el sur andino y con categoría de conservación nacional de “vulnerable”, registrada para el PNC permitirá diseñar y detectar nuevas especies y comprender su distribución, además de concientizar a la población sobre la importancia en la conservación y expansión de esta especie amenazada.

Agradecimientos

Nuestro especial agradecimiento al equipo técnico de la empresa de investigación y servicios Inca’Biotec SAC- Perú y al personal del Parque Nacional de Cutervo por las facilidades brindadas en la ejecución del presente estudio.

ORCID

K. A. Peña Yacila  https://orcid.org/0000-0002-3804-3837

E. A. Suarez-Peña  https://orcid.org/0000-0003-0137-8251

C. A. Torres Pera  https://orcid.org/0000-0003-1608-6622

C. A. Bermejo Requena  https://orcid.org/0000-0001-5294-7034

L. X. Llacsa Sánchez  https://orcid.org/0000-0002-7463-4394

O. J. Paredes-Vilca  https://orcid.org/0000-0003-2301-3675

J. Y. Dávila García  https://orcid.org/0000-0002-1310-1690

Referencias bibliográficas

Albán-Castillo, J., Chilquillo, E., Melchor-Castro, B., Arakaki, M., León, B., & Suni, M. (2020). Cinchona L.:" Arbol de la Quina": repoblamiento y reforestacion en el Peru. Revista peruana de biología, 27(3), 423-427.

Albán, J. (2013). Etnobotánica de rubiáceas peruanas (Doctoral dissertation, Tesis para optar al grado de Doctor. Escuela de Posgrado, Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima-Perú. 617 p.

Andersson, L., & Antonelli, A. (2005). Phylogeny of the Tribe Cinchoneae (Rubiaceae), Its Position in Cinchonoideae, and Description of a New Genus, Ciliosemina. Taxon, 54(1), 17-28.

Andersson, L. (1995). Tribes and genera of the Cinchoneae complex (Rubiaceae). Annals of the Missouri Botanical Garden, 409-427.

Antonelli, A., Nylander, J. A. A., Persson, C., & Sanmartín, I. (2009). Tracing the impact of the Andean uplift on Neotropical plant evolution. Proceedings of the National Academy of Sciences of the United States of America, 106(24), 9749-9754.

Aymard, G. A. (2019). Breve reseña de los aspectos taxonómicos y nomenclaturales actuales del género Cinchona (Rubiaceae-Cinchoneae). Revista de la Academia Colombiana de Ciencias Exactas, Físicas y Naturales, 234-241.

Caruso, G., Broglia, V., & Pocovi, M. (2015). Diversidad genética. Importancia y aplicaciones en el mejoramiento vegetal. Publicación del Instituto de Ecología y Ambiente Humano, 4(1), 45-50.

Chilquillo, E. A. (2016). Filogenia do gênero Ladenbergia Klotzsch (Rubiaceae: Cinchoneae) e taxonomia atualizada das espécies que ocorrem nos Andes tropicais. Mestre em Biologia Vegeta, Instituto de Biologia da Universidade Estadual de Campinas. Brasil.

Collada, C., & Jiménez, P. (2008). Técnicas para la evaluación de la diversidad genética y su uso en los programas de conservación. Forest Systems, 9(4), 237-248. DOI: 10.5424/706.

Cueva-Agila, A., Vélez-Mora, D., Arias, D., Curto, M., Meimberg, H., & Brinegar, C. (2019). Genetic characterization of fragmented populations of Cinchona officinalis L. (Rubiaceae), a threatened tree of the northern Andean cloud forests. Tree Genetics and Genomes, 15(6), 81.

De Andrade, J. L. (2013). The concept of biodiversity and the history of conservation biology: from wilderness preservation to biodiversity conservation. História, 32(2), 21-47.

Doyle, J., & Doyle, L. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin, 11-15.

Ezenwa, V. O., Godsey, M. S., King, R. J., & Guptill, S. C. (2006). Avian diversity and West Nile virus: Testing associations between biodiversity and infectious disease risk. Proceedings of the Royal Society B: Biological Sciences, 273(1582), 109-117.

García, E. C. (2012). Mecanismos de especiación ecológica en plantas y animales. Biológicas Revista de la DES Ciencias Biológico Agropecuarias, 14(2), 7-13.

Gómez, A., Beraun, L. A., Gómez, O. J., & Llatas, E. (2016). Procesos de regeneración natural de la quina o cascarilla (Cinchona spp.) en los bosques de neblina del distrito de Kañaris, region Lambayeque. INIA Instituto Nacional de Innovación Agraria. Ministerio de Agricultura. MINAGRI. Lima-Perú

Grant, W., & Bowen, B. . (1998). Shallow Population Histories in Deep Evolutionary Lineages of Marine Fishes: Insights From Sardines and Anchovies and Lessons for Conservation. The Journal of Heredity, 89(5), 415-426.

Hilu, K. W., & Liang, H. (1997). The matK gene sequence variation and application in plant systematics. American Journal of Botany, 84(6), 830-839.

Hoffmann, M., Hilton-Taylor, C., Angulo, A., Böhm, M., Brooks, T. M., et al. (2010). The impact of conservation on the status of the world’s vertebrates. Science, 330(6010), 1503-1509.

Huamán, L., Albán, J., & Chilquillo, E. (2019). Aspectos Taxonómicos y avances en el conocimiento del estado actual del árbol de la Quina (Cinchona officinalis L.) en el Norte de Perú. Ecología Aplicada, 18(2), 145.

Johnson, P. T. J., & Thieltges, D. W. (2010). Diversity, decoys and the dilution effect: How ecological communities affect disease risk. Journal of Experimental Biology, 213(6), 961-970.

Keesing, F., Belden, L. K., Daszak, P., Dobson, A., Harvell, C. D., et al. (2010). Impacts of biodiversity on the emergence and transmission of infectious diseases. Nature, 468(7324), 647-652.

Killeen, T., García, E., & Beck, S. G. (1993). Guía de Árboles de Bolivia. Publicado por Herbario Nacional de Bolivia Missouri Botanical Garden. Impresores Quipus S.R.L., La Paz, Bolivia. Depósito Legal, 693-93. 958.

Little, D. P. (2014). A DNA mini-barcode for land plants. Molecular Ecology Resources, 14(3), 437-446.

López, N. (2016). Evaluación del paisaje y recursos escénicos después de 350 años de explotación de la “cascarilla” o “quina” Cinchona officinalis L. (Rubiaceae) en el sector Cajanuma-Rumishitana, Ecuador. Arnaldoa, 23(2), 461-474.

Maldonado, C., Barnes, C. J., Cornett, C., Holmfred, E., Hansen, S. H., et al. (2017). Phylogeny predicts the quantity of anti-malarial alkaloids within the iconic yellow Cinchona bark (Rubiaceae: Cinchona calisaya). Frontiers in Plant Science, 8, 391.

Manns, U., & Bremer, B. (2010). Towards a better understanding of intertribal relationships and stable tribal delimitations within Cinchonoideae s.s. (Rubiaceae). Molecular Phylogenetics and Evolution, 56(1), 21-39.

Marchese, C. (2015). Biodiversity hotspots: A shortcut for a more complicated concept. Global Ecology and Conservation, 3, 297-309.

Martins, I. S., Navarro, L. M., Pereira, H. M., & Rosa, I. M. D. (2020). Alternative pathways to a sustainable future lead to contrasting biodiversity responses. Global Ecology and Conservation, 22, e01028.

Mendoza, H., Ramírez, B., & Jiménez, L. C. (2004). Rubiaceae de Colombia. Guía ilustrada de géneros. Bogotá, Colombia: Instituto de Investigación de Recursos Biológicos Alexander von Humboldt. ISBN: 8151-31-7.

MINAGRI (2006). Aprueban Categorización de Especies Amenazadas de Flora Silvestre. Decreto Supremo N° 043-2006-AG. Ministerio de Agricultura - MINAGRI. Limas-Perú.

Mostacero, L., Mejía , C., & Gamarra, T. (2009). Fanerógamas del Perú - Taxonomía, utilidad y ecogeografía, Universidad Nacional de Trujillo. CONCYTEC, 1331. ISBN: 978-9972-50-099-2

Nair, K. P. (2010). The Agronomy and Economy of Important Tree Crops of the Developing. Elsevier, 111-129.

Pimm, S. L., Russell, G. J., Gittleman, J. L., & Brooks, T. M. (1995). The future of biodiversity. Science, 269(5222), 347-350.

Piñero, D., Caballero-Mellado, J., Cabrera-Toledo, D., Canteros, C. E., Casas, A., & Castañeda-Sortibrán, A. (2008). La diversidad genética como instrumento para la conservación y el aprovechamiento de la biodiversidad: estudios en especies mexicanas. Capital natural de México, 1, 437-494.

Rawat, U. S., & Agarwal, N. K. (2015). Biodiversity: Concept, threats and conservation. Environment Conservation Journal, 16(3), 19-28.

Salvador, L., Alía, R., Agúndez, D., & Gil, L. (2000). Genetic variation and migration pathways of maritime pine (Pinus pinaster Ait) in the Iberian peninsula. Theoretical and Applied Genetics, 100(1), 89-95.

Sánchez, M., & Calderón, A. (2010). Evaluación preliminar de orquídeas en el Parque Nacional Cutervo, Cajamarca-Perú. Ecología Aplicada, 9(1-2), 1.

SERNANP. (2014). Diagnóstico del Proceso de Elaboración del Plan Maestro 2009-2014. Servicio Nacional de Areas Naturales Protegidas por el Estado (SERNANP). Lima, Perú.

SERNANP. (2017). Actualización de Plan Maestro del Parque Nacional de Cutervo 2017 - 2021. Aprobado mediante Resolución Presidencial RP N° 031-2017-SERNANP. Servicio Nacional de Areas Naturales Protegidas por el Estado (SERNANP). Lima - Perú.

Zevallos, P. A. (1989). Taxonomía, distribución geográfica y status del género Cinchona en el Perú. Centro de datos para la conservación Facultad de Ciencias Forestales. Edit. Virginia Isayama Okamoto.