Introducción

El cambio climático (CC) acelerado está ocasionando la extinción rápida de ectotermos a nivel mundial (^{Carpenter et al., 2008}; ^{Thomas et al., 2004}) y uno de los grupos con mayor riesgo de extinción es el de las lagartijas (^{Huey et al., 2010}; ^{Sinervo et al., 2010}). Algunas especies de lagartijas son heliotermas y regulan su temperatura corporal (T _b ) exponiéndose u ocultándose de la radiación solar (^{Hertz et al., 1993}; ^{Huey y Slatkin, 1976}). Otras especies son termoconformistas y su T _b varía de acuerdo con la temperatura ambiental (^{Huey et al., 2009}). En ambos casos, cuando la temperatura ambiental es alta, los individuos se refugian para evitar que su T _b exceda su límite superior de tolerancia térmica (temperatura crítica máxima; Ct _max ) y mueran por sobrecalentamiento (^{Cowles y Bogert, 1944}; ^{Sinervo et al., 2017}). Sin embargo, si los individuos permanecen refugiados (horas de restricción térmica; h _r ) durante su horario de actividad (h _a ) el tiempo para que realicen sus actividades biológicas básicas se acorta (e.g., tiempo para forrajeo), lo cual es especialmente importante durante la temporada de reproducción, ya que las tasas de reproducción descienden, el reclutamiento poblacional disminuye o no existe y en consecuencia las poblaciones decrecen hasta su extinción (^{Sinervo et al., 2011}).

^{Sinervo et al. (2010}) desarrollaron un modelo "genérico" que incorpora aproximaciones mecanicistas y correlativas para proyectar el riesgo de extinción de lagartijas heliotermas o termoconformistas por medio del cálculo de las h _r y encontraron que el 20% de las especies de lagartijas a nivel mundial se podrían extinguir para el año 2080, debido al aumento de la temperatura ambiental. De acuerdo con las proyecciones, algunas especies de lagartijas podrían tener un riesgo de extinción de hasta el 91% para el año 2080 (e.g., Helodermatidae), mientras que otras especies tendrían un riesgo de extinción nulo (e.g., Gekkonidae). Las diferencias en la susceptibilidad al CC entre especies se deben a 2 aspectos principales. El primero es la sensibilidad térmica, las especies con T _b < 30 °C son más susceptibles que las especies con T _b > 30 °C (^{Sinervo et al., 2010}). El segundo es la distribución geográfica, a menor altitud o latitud, la temperatura ambiental suele ser más alta y las especies suelen presentar un mayor riesgo de extinción (^{Deutsch et al., 2008}; ^{Sinervo et al., 2018}).

A pesar de que las proyecciones generales de riesgo de extinción a nivel de familias de lagartijas son alarmantes, éstas podrían subestimar la probabilidad de extinción a nivel de especie. Por ejemplo, ^{Sinervo et al. (2010}) proyectaron que el riesgo de extinción para los gecos de la familia Phyllodactylidae será nulo para el año 2080, pero ^{Kubisch et al. (2016}) proyectaron que el riesgo de extinción del geco Homonota darwinii (Phyllodactylidae) será del 20% para el mismo año. Por lo tanto, consideramos que es necesario proyectar y validar el modelo de riesgo de extinción a nivel específico, principalmente para las especies que, por su sensibilidad térmica y su distribución geográfica, podrían ser altamente susceptibles al CC.

En este trabajo, seleccionamos como modelo de estudio al geco diurno semiarborícola Gonatodes concinnatus (Sphaerodactylidae; Fig. 1), el cual habita en la región occidental de la selva amazónica de Ecuador y Perú (^{Sturaro y Ávila-Pires, 2011}). Gonatodes concinnatus alcanza una longitud hocico-cloaca (LHC) de entre 47.6 mm (♀) y 48.3 mm (♂), tiene una T _b promedio < 30 °C (^{Fitch, 1968}) y habita en un intervalo altitudinal de 100-600 m y latitudinal de 0-5° (^{Sturaro y Avila-Pires, 2011}). De acuerdo con las proyecciones de ^{Sinervo et al. (2010}), el riesgo de extinción para las especies de la familia Sphaerodactyilidae será nulo para 2050 y 2080. Sin embargo, por su baja T _b y por su estrecha distribución altitudinal y latitudinal planteamos la hipótesis de que G. concinnatus tendrá un riesgo de extinción alto ante el CC. Por lo tanto, los objetivos del presente trabajo fueron evaluar la estrategia de termorregulación de G. concinnatus, proyectar su riesgo de extinción para 2050 y 2070 bajo un escenario de cambio climático y definir si el riesgo de extinción está correlacionado con la distribución altitudinal de las poblaciones.

Figura 1 Hembra (a) y macho (b) de Gonatodes concinnatus fotografíados en la Estación Científica Amazónica Juri-Juri Kawsay. 

Materiales y métodos

Realizamos trabajo de campo durante julio de 2017 en los alrededores de la Estación Científica Amazónica Juri-Juri Kawsay de la Universidad Central del Ecuador, ubicada en el Bosque Protector del Oglán Alto, Cantón Arajuno, provincia de Pastaza, Ecuador (01°19'26.95" S, 77°41'18.9" O; 604 m snm). El tipo de vegetación corresponde a un bosque pluvial premontano caracterizado por árboles con copas emergentes, entre los que destacan Ceiba pentandra, Pachira insignis, Ficus perisiana y Otoba parviflora, donde hay abundantes musgos y hepáticas en las hojas y ramas de la vegetación arbórea y arbustiva (^{Albuja et al., 2012}; ^{Cañadas-Cruz, 1983}; ^{Cerón et al., 2007}). La temperatura promedio anual en el sitio de estudio varía entre 18 y 24 °C y la precipitación oscila entre los 4,000 y 8,000 mm anuales (^{Cerón et al., 2007}).

Para estimar los requerimientos térmicos y la estrategia de termorregulación de G. concinnatus (Fig. 1) recolectamos de forma directa individuos adultos (~ 48 mm LHC; ^{Carvajal-Campos, 2017}) entre las 0700-1900 h, horario establecido con base en la información sobre actividad diurna en otras especies del género Gonatodes (^{Gamble et al., 2015}; ^{Vitt et al., 2000}) y confirmado para este taxón durante el trabajo de campo. Inmediatamente después de capturar a cada individuo (en un lapso no mayor a 10 s) registramos su T _b (cloacal) mediante un termopar tipo K conectado a un termómetro digital de lectura rápida (Fluke 51-II; ± 0.1 °C). Posteriormente, trasladamos a todas las lagartijas recolectadas a las instalaciones de la estación científica Juri-Juri Kawsay y las colocamos en un gradiente térmico dentro del intervalo de tiempo de 0900-1600 h para registrar la temperatura seleccionada (T_sel; ^{Hertz et al., 1993}) de cada individuo. El gradiente térmico consistió en una caja de policarbonato (100 cm de ancho χ 100 cm de largo χ 50 cm de alto) dividida en 10 carriles, en el cual ofrecimos a las lagartijas un intervalo de temperatura de ~ 15 °C a ~ 40 °C. Para obtener el extremo frío colocamos bolsas de hielo debajo de un extremo del gradiente y para obtener el extremo caliente colocamos 10 focos incandescentes de 100 W suspendidos a 30 cm de altura en el otro extremo del gradiente (un foco por carril). Asignamos un carril a cada individuo para evitar interacción conductual (^{Lara-Reséndiz et al., 2015}) y registramos su T _sel cada hora con el mismo termómetro digital de lectura rápida usado en el trabajo de campo. Al terminar los registros de T _sel liberamos a los ejemplares en el sitio exacto donde los recolectamos.

Durante 6 meses consecutivos (20 de julio al 21 de diciembre de 2017), registramos las temperaturas operativas (T _e ;^{Bakken, 1992}) en los microhábitats usados por G. concinnatus en el sitio de estudio. La T _e se define como el producto de las interacciones entre factores biofísicos y morfológicos que influye en la temperatura corporal de un ectotermo y representa la temperatura de equilibrio que tendrían los individuos en ausencia de termorregulación conductual y fisiológica (^{Hertz et al., 1993}; ^{Lara-Reséndiz et al., 2015}). Para registrar la T _e , utilizamos 2 modelos biofísicos previamente calibrados que consistieron en tubos de policloruro de vinilo (PVC) de 100 mm de largo y 20 mm de diámetro, pintados de color gris y que emularon las tasas de ganancia y pérdida de calor (r = 0.89, p < 0.05) de los individuos adultos de la especie en estudio (ver detalles de calibración en ^{Arenas-Moreno et al. [2018}] y ^{Lara-Reséndiz et al. [2015]}). Dentro de cada modelo biofísico colocamos un sensor de temperatura ambiental modelo Hobo® U23-003 (HOBO Temp, ONSET Computer Corporation, Massachusetts) programado para registrar la temperatura cada 10 min durante el horario de actividad de los individuos. Con el objetivo de representar la variación de la temperatura microambiental colocamos los modelos de PVC en 2 sitios diferentes. El primero en un tronco bajo sombra de vegetación natural y el segundo en un tronco con exposición entre sol y sombra a través del día. Ambos microhábitats son lugares de actividad de G. concinnatus de acuerdo con nuestras observaciones durante el trabajo de campo. Es importante mencionar que, en los sitios de latitud o altitud baja, la temperatura ambiental es constante y poco variable a lo largo del año (^{Huey et al., 2009}; ^{Janzen, 1967}) y la T _e es relativamente homogénea en tiempo y espacio (^{Hertz, 1992}; ^{Medina et al., 2016}). Por lo anterior, consideramos que los datos obtenidos durante 6 meses con 2 modelos nulos representan adecuadamente la variación anual de la T _e .

Adicionalmente, registramos la temperatura ambiental (temperatura del aire a 2 m del suelo) en el sitio de estudio durante el mismo periodo e intervalos en que registramos la T _e , mediante un registrador de temperatura ambiental modelo Hobo® U23-002. El registro de la temperatura ambiental lo utilizamos para validar el uso del método senoidal en el modelo de riesgo de extinción descrito a continuación.

Para proyectar el riesgo de extinción de G. concinnatus utilizamos el modelo mecanicista y correlativo de riesgo de extinción (MRE) propuesto por ^{Sinervo et al. (2010}), el cual está basado tanto en registros de presencia como en atributos ecofisiológicos. Para ello, empleamos la paquetería Mapinguari Versión 0.4.1 (^{Caetano et al., 2017}) del lenguaje de programación R (^{R Development Core Team, 2018}). El método de Sinervo et al. (2010) proyecta la probabilidad de presencia (distribución potencial) al presente y la probabilidad de presencia al 2050 y 2070 de lagartijas heliotermas o termoconformistas con base en el cálculo de horas de actividad (h _a ) y de restricción (h _r ) de los individuos durante su temporada de reproducción, ambas estimadas a partir de información microclimática, variables climáticas, distribución geográfica y datos fisiológicos de los organismos (^{Pontes-da Silva et al., 2018}). Para más información del método ver detalles en ^{Sinervo et al. (2010}, ²⁰¹¹, ²⁰¹⁷) y ^{Caetano et al. (2017)}. Para construir el MRE empleamos el método senoidal que estima las h _a y h _r bajo el supuesto de que los organismos son termoconformistas y mantienen una T _b igual o cercana a la temperatura del aire (T_air). Este supuesto fue validado por medio de la estimación de la eficiencia en la termorregulación de los individuos y mediante un análisis de correlación de Spearman entre la temperatura ambiental del sitio de estudio y la T _e del modelo biofísico 1 (r = 0.98) y la temperatura ambiental y la T _e del modelo biofísico 2 (r = 0.96).

Los datos de presencia para estimar el MRE fueron obtenidos de los registros de distribución de G. concinnatus considerando la localidad de estudio y otras 39 localidades alejadas entre sí por al menos 2 km (criterio de ^{Sinervo et al., 2010}), de acuerdo con las localidades mencionadas para la especie de estudio en la revisión taxonómica de ^{Sturaro y Ávila-Pirez (2011}) y de las localidades registradas en el portal FaunaWebEcuador-Reptiles (^{Torres-Carvajal, et al., 2018}; disponible en https://bioweb.bio/faunaweb/reptiliaweb/). Después obtuvimos las capas bioclimáticas de altitud y de temperatura mensual promedio mínima y máxima (T _min y T _max ) y precipitación anual para los periodos presente y futuro (2050 y 2070) a una resolución espacial de 2.5 minutos de arco. La información climática fue descargada del portal WorldClim versión 1.4 (^{Hijmans et al., 2005}; http://www.worldclim.org/version1) para el presente (mediante interpolaciones de datos de 1960 a 1990) y para el futuro (2050 y 2070) a partir de capas generadas por el Instituto de Meteorología MaxPlanck (MPI-ESM-LR) correspondientes a la trayectoria de concentración representativa 4.5 (Representative Concentration Pathway, RCP 4.5), misma que implica un escenario de mitigación de gases de efecto invernadero, con un pico de emisión al año 2100, seguido de una estabilización en el forzamiento radiativo y una anomalía térmica media 2.4 °C por encima de los niveles preindustriales (^{Thomson et al., 2011}). Posteriormente, indicamos en el MRE el límite inferior (T _lower ) y superior (T _upper ) de tolerancia térmica de G. concinnatus, que de acuerdo con el criterio de ^{Pontes-da Silva et al. (2018}) corresponde con la T _b mínima y máxima de individuos activos en campo (20.4 °C y 29.8 °C, respectivamente). Debido a que la especie de estudio presenta reproducción continua (^{Carvajal-Campos, 2017}) creamos capas ecofisiológicas de las h _a y las h _r de G. concinnatus, considerando los 12 meses del año. Estimamos la variación diaria de la T _air asumiendo una curva senoidal entre la T _min . y la T _max y definimos las h _a como el número de horas en que T _air > 20.4 °C (T _lower ) y las h _r como el número de horas en que T _air > 29.8 °C (T _upper ) (^{Caetano et al., 2017}; ^{Sinervo et al., 2017}).

Estimamos la probabilidad de presencia para cada pixel dentro del área de estudio para los periodos presente y futuro por medio de un modelo lineal generalizado (GLM, por sus siglas en inglés) y un modelo aditivo generalizado (GAM, por sus siglas en inglés) utilizando el paquete biomod2 (^{Thuiller et al., 2016}), para lo cual usamos como insumos las capas ecofisiológicas de h _a y h _r , las capas bioclimáticas y los puntos de presencia de la especie con 1000 pseudoausencias (recomendadas por ^{Barbet-Massin et al. (2012}) para modelos GLM y GAM) y validados por medio del criterio de la curva de "característica operativa del receptor" (ROC, por sus siglas en inglés), con lo cual seleccionamos únicamente los modelos con valores mayores a 0.75 (^{Pontes-da Silva et al., 2018}; ^{Vicenzi et al., 2017}). Para estimar el porcentaje de riesgo de extinción comparamos la probabilidad de presencia al presente (de los 40 registros de distribución G. concinnatus utilizados en este estudio) con la probabilidad de presencia en 2050 y 2070 (Pontes-da Silva et al., 2018).

Calculamos el promedio de la calidad térmica del hábitat (de¯), el promedio de la precisión térmica (db¯) y la eficiencia en la termorregulación de los individuos (E) mediante el protocolo de ^{Hertz et al. (1993}; E = 1-( db¯ / de¯) y Blouin-Demmers y Weatherhead (2001; E = de¯-db¯). Para calcular los índices propuestos por Hertz et al. (1993) y ^{Blouin-Demers y Weatherhead (2001}) utilizamos los registros de T _b de los individuos, los registros de T _e y el intervalo set-point (T_set; cuartiles del 25 % y 75 % de los datos centrales) de las T _sel (T _set25 y T _set75 , respectivamente).

En primer lugar, de representa el promedio de la diferencia absoluta de cada T _e con respecto a la T _set , donde los valores cercanos a cero indican calidad térmica del hábitat alta, es decir, que la temperatura del hábitat está dentro o cerca del intervalo de T _set de los individuos. Si T _e < T _set25 , entonces de = T _set25 -T _e ; si T _e está dentro de los cuartiles de T _set , entonces de = 0; y si T _e > T _set75 , entonces de = T _e -T _set75 . En segundo lugar, db¯ representa el promedio de la diferencia absoluta de la T _b de cada individuo con respecto a la T _set , donde los valores cercanos a cero indican una precisión térmica alta, es decir, los individuos mantienen su T _b dentro o cerca del intervalo de T _sef Si la T _b < T _set25 , entonces db = T _set25 -T _b ; si T _b está dentro de los cuartiles de T _set , entonces db = 0; y si T _b > T _set75 , entonces db = T _b -T _set75 Finalmente, E representa qué tan eficientes son los individuos para termorregular. Utilizamos 2 protocolos para estimar E, ya que tienen enfoques ligeramente distintos (^{Lara-Reséndiz et al., 2015}). Por una parte, ^{Hertz et al. (1993}) establecieron que los valores de E cercanos a 0 indican termoconformismo y los valores cercanos a 1 indican termorregulación activa. Por otra parte, ^{Blouin-Demers y Weatherhead (2001}) establecieron que E representa un valor absoluto (°C) del esfuerzo de termorregulación de los individuos. En ambos casos, los valores negativos sugieren una baja explotación del ambiente térmico en relación con la oferta térmica disponible.

Mediante análisis descriptivos calculamos el promedio y la desviación estándar (±) de T_b, T_e, temperatura ambiental del sitio de estudio y T _se¡ . Finalmente, realizamos un análisis de correlación de Spearman entre el valor de riesgo de extinción para 2050 y 2070 de 40 pixeles (correspondientes con los registros de distribución de la especie de estudio) y la altitud, para poner a prueba si el riesgo de extinción cambia a través de la elevación.

Resultados

La T _b promedio de los gecos (n = 18) fue de 26.8 ± 2.5 °C (intervalo de 20.4-29.8 °C; Fig. 2), la T _e (n = 7690) fue de 24.6 ± 2.9 °C (intervalo de 17.9-35.9 °C; Fig. 2), la T _sel (n = 18) fue de 24.6 ± 1.6 °C con cuartil T _set25 de 23.4 °C y cuartil T _set75 de 25.7 °C (Fig. 2) y la temperatura ambiental (n = 3845) fue de 22.95 ± 2.6 °C. La calidad térmica del hábitat (de¯) en el sitio de estudio fue de 1.42 °C y la precisión térmica (db¯) de los organismos fue de 2.01 °C. De acuerdo con el protocolo de ^{Hertz et al. (1993}), el índice de eficiencia en la termorregulación (E) fue de -0.41 y de acuerdo con el protocolo de ^{Blouin-Demers y Weatherhead (2001}) fue de -0.59.

Figura 2 Intervalo de temperaturas corporales (T_b; sombreado diagonal), cuartiles de temperatura seleccionada (T _se225 y T _se775 ; líneas punteadas verticales) y temperaturas operativas (T_e; barras verticales) de Gonatodes concinnatus en la Estación Científica Amazónica Juri-Juru Kawsay. 

La probabilidad de presencia al presente (Fig. 3A) es mayor que la proyectada a 2050 (Fig. 3B) y a 2070 (Fig. 3C). El porcentaje de riesgo de extinción para las poblaciones analizadas de G. concinnatus tuvo un intervalo de 12-97% (promedio de 84.2%) para 2050 y de 26-97% (promedio de 86.4%) para 2070 (Tabla 1). La población del Oglán Alto fue la que presentó un menor porcentaje de riesgo de extinción ante el CC, con 12 y 26%, a 2050 y 2070, respectivamente (símbolo de estrella en la Figura 3). Encontramos una correlación negativa entre la altitud de las poblaciones y su riesgo de extinción a 2050 y 2070 (r = -0.93, p < 0.05 y r = -0.96, p < 0.05, respectivamente), es decir, a menor altitud, mayor riesgo de extinción (Fig. 4).

Figura 3 Proyecciones de probabilidad de presencia de Gonatodes concinnatus. A) Presente; B) año 2050 (RCP 4.5); C) año 2070 (RCP 4.5); D) mapa de referencia del área proyectada (rectángulo púrpura). La escala de color representa la probabilidad de presencia en porcentaje (tonalidades frías: alta probabilidad; tonalidades cálidas: baja probabilidad). Los puntos indican registros de la especie y la estrella indica el sitio de estudio. 

Figura 4 Análisis de correlación de Spearman entre el porcentaje de riesgo de extinción y la altitud de las poblaciones de Gonatodes concinnatus. Los círculos negros representan el riesgo de extinción para el año 2050 (RCP 4.5) y los blancos para el año 2070 (RCP 4.5). 

Discusión

En el presente trabajo determinamos la estrategia de termorregulación de una población del geco G. concinnatus y proyectamos su riesgo de extinción para los años 2050 y 2070 (RCP 4.5) a través de su distribución geográfica en la región occidental de la selva amazónica. Junto con el trabajo de ^{Pontes-da-Silva et al. (2018}), el presente estudio ofrece una de las primeras proyecciones de riesgo de extinción por CC a nivel de especie en lagartijas de la selva amazónica bajo un enfoque correlativo, mecanicista y biofísico. Además, es el primero en mostrar a nivel intraespecífico una correlación negativa entre el riesgo de extinción y la altitud. Consideramos que nuestros resultados podrían servir como base para nuevos estudios de vulnerabilidad ecofisiológica al CC en organismos ectotérmicos.

La conducta de termorregulación es un continuo de estrategias conductuales con 2 extremos: la termorregulación activa y el termoconformismo (^{Hertz et al., 1993}). Por una parte, la termorregulación activa ocurre cuando los individuos mantienen su T _b dentro o cerca de su intervalo de T _set (valores de db¯ cercanos a 0) en ambientes donde la calidad térmica del hábitat es baja. En el caso de termorreguladores activos el valor de db siempre es menor que de¯ y se ve reflejado en valores de E cercanos a 1 (^{Artacho et al., 2017}; Hertz et al., 1993) o en diferencias absolutas (positivas) altas entre de¯ y db¯ (^{Blouin-Demers y Weatherhead, 2001}). Por otra parte, el termoconformiso ocurre cuando la T _b de los individuos es similar a la temperatura ambiental. En el caso de termoconformistas el valor de db¯ siempre es similar al de de¯ y se ve reflejado con valores de E (o en diferencias absolutas entre de¯ y db¯) cerca de 0 (Blouin-Demers y Weatherhead, 2001; Hertz et al., 1993). El termoconformismo puede ocurrir debido a 2 razones principales. La primera es que el costo de la termorregulación activa supera los beneficios (^{Huey y Slatkin, 1976}). La segunda es que los individuos habitan ambientes con calidad térmica alta y no necesitan termorregular. En el caso de G. concinnatus el ambiente ofrece una calidad térmica alta (de¯ = 1.42 °C) y por ello los individuos no necesitan termorregular de forma activa, ya que el ambiente ofrece temperaturas cercanas a su T _set . Sin embargo, encontramos que la precisión térmica de los individuos (db¯ = 2.01 °C) fue ligeramente menor a la calidad térmica del hábitat y como resultado el valor de eficiencia en la termorregulación resultó negativo (E = -0.41; Hertz et al., 1993) con una diferencia absoluta entre de¯ y db¯ de -0.58 °C (Blouin-Demers y Weatherhead, 2001). De acuerdo con Hertz et al. (1993) y con Blouin-Demers y Weatherhead (2001) los resultados negativos de E se pueden deber a que los individuos evaden temperaturas microambientales idóneas, como una posible forma de evasión de depredadores o por competencia intraespecífica. Consideramos que en el caso de G. concinnatus los valores negativos de E también se podrían deber a que los individuos eligen microhabitats húmedos, pero térmicamente poco favorables, ya que de forma general las lagartijas de la familia Sphaerodactylidae presentan tasas altas de deshidratación (^{Snyder, 1979}; ^{Nava, 2004}). Es importante considerar que el tamaño de muestra y el periodo de estudio podrían limitar nuestras observaciones, sin embargo, nuestros resultados son apoyados con trabajos previos de especies tropicales (Hertz et al., 1993; ^{Tewksbury et al., 2008}). Consideramos que G. concinnatus es un termoconformista, pues el valor de E es cercano a cero (termoconformismo) con una diferencia absoluta entre de¯ y db¯ que no supera 1 °C. Además, nuestra consideración es consistente con resultados registrados en otras especies de gecos de la familia Sphaerodactylidae, también consideradas termoconformistas y con patrones generales de termorregulación en lagartijas que señalan que, en sitios tropicales, las lagartijas tienen tendencia termoconformista (^{Huey et al., 2009}; ^{Medina et al., 2016}).

En un estudio previo (^{Sinervo et al., 2010}), se proyectó que la probabilidad de extinción ante el CC para poblaciones de gecos de la familia Sphaerodactylidae es nula para el presente y futuro; sin embargo, nuestros resultados indican que la mayoría de las poblaciones de G. concinnatus tienen una probabilidad de extinción superior a 80% para 2050 y 2070 (Tabla 1). La amplia diferencia entre los resultados de las proyecciones de riesgo de extinción a nivel de familia y a nivel de especies particulares denota la importancia de evaluar la vulnerabilidad de lagartijas a niveles taxonómicos específicos (e.g., ^{Kubisch et al., 2016}; ^{Lara-Reséndiz et al., 2015}; ^{Pontes-da-Silva et al., 2018}). La diferencia entre las proyecciones de ^{Sinervo et al. (2010)} para la familia Sphaerodactylidae y nuestras proyecciones para G. concinnatus se deben principalmente a los valores distintos de T _lower y T _upper utilizados entre ambas proyecciones. Para proyectar el riesgo de extinción de la familia Sphaerodactylidae, ^{Sinervo et al. (2010)} utilizaron como T _lower y T _upper los valores de T _b más bajos y altos registrados para la familia Sphaerodactylidae; 25.3 °C para Teratoscincus przewalskii (^{Dimitri y Borkin, 1992}) y 38.6 °C para Pristurus carteri (^{Arnold, 1993}); mientras que nosotros, al proyectar el riesgo de extinción de G. concinnatus consideramos apropiado utilizar como T _lower y T _upper los valores de T _b mínima (20.4 °C) y máxima (29.8 °C) de individuos activos en campo (^{Kubisch et al., 2016}; ^{Pontes-da-Silva et al., 2018}). En otras palabras, las diferencias entre las proyecciones de riesgo de extinción de ^{Sinervo et al. (2010)} y las nuestras se deben a que los límites de tolerancia térmica son mayores entre especies que entre individuos de la misma especie. Al utilizar límites de tolerancia térmica amplios, las proyecciones pueden ser precisas a nivel de familia, pero se puede subestimar el riesgo de extinción intraespecífico (^{Kubisch et al., 2016}; ^{Pontes-da-Silva et al., 2018}). Es importante señalar que delimitamos T _lower y T _upper para G. concinnatus con base en datos de individuos de una sola población, por ello es posible que estemos omitiendo variación interpoblacional (^{Angert et al., 2011}), lo cual podría tener efecto en nuestras proyecciones. Sin embargo, la T _b documentada en otras poblaciones de G. concinnatus (^{Fitch, 1968}) está dentro del intervalo que observamos en la población de estudio e inclusive es similar a la de otras especies del género Gonatodes (^{Brusch et al., 2016}; ^{Vitt et al., 2000}), por lo cual consideramos que el intervalo de tolerancias térmicas que utilizamos es adecuado.

El riesgo de extinción elevado de G. concinnatus concuerda con el patrón típico de susceptibilidad ante el CC de ectotermos de baja latitud y de ambientes tropicales (^{Deutsch et al., 2008}). Por ejemplo, ^{Huey et al. (2009}) realizaron un análisis macrofisiológico de especies de lagartijas que habitan en diferentes latitudes y encontraron que las especies de ambientes tropicales son las más vulnerables al calentamiento global. La alta susceptibilidad de G. concinnatus al CC se debe a que en los ambientes tropicales existe poca heterogeneidad térmica (^{Hertz, 1992}; ^{Janzen, 1967}; ^{Medina et al., 2016}) y las lagartijas tienen pocas opciones de regular su T _b de forma conductual, de manera que presentan tendencia al termoconformismo (^{Huey et al., 2009}). Además, los ectotermos de ambientes tropicales suelen presentar una T _b cercana a su Ct _max (margen de seguridad térmica bajo; ^{Sunday et al., 2014}), son especialistas térmicos con capacidades de aclimatación limitadas (^{Tewksbury et al., 2008}) y un aumento de la temperatura ambiental podría ocasionar que los individuos mueran por sobrecalentamiento o aumenten sus horas de restricción (^{Sinervo et al., 2010}).

Al correlacionar el riesgo de extinción de las poblaciones de G. concinnatus con la altitud, encontramos que a menor elevación mayor riesgo de extinción. Nuestros resultados a nivel intraespecífico son similares a los resultados encontrados entre especies. Por ejemplo, ^{Sinervo et al. (2018}), proyectaron que el riesgo de extinción de las especies de lagartijas del género Phrynocephalus disminuye de menor a mayor altitud. Por lo anterior, podríamos considerar que los sitios de mayor elevación serían sitios térmicamente adecuados para la especie de estudio y probablemente para otros ectotermos con características fisiológicas similares. Del total de localidades para las que proyectamos el riesgo de extinción de G. concinnatus, únicamente la población de mayor elevación tendría un riesgo de extinción bajo ante el CC, incluso para 2070; mientras que el resto de las poblaciones tendrían un riesgo de extinción alto. En este punto es importante señalar que es posible que aún existan poblaciones de G. concinnatus no documentadas en la selva amazónica, especialmente en la región norte de Perú. Sumado al riesgo de extinción ante el CC, existe una elevada tasa de deforestación de la selva amazónica (^{Coe et al., 2013}), la cual podría acelerar la extinción de diversas especies.

A pesar del escenario desalentador de G. concinnatus ante el CC, se debe considerar que podrían existir al menos 3 mecanismos de respuesta por parte de los organismos. El primero es mediante cambios en la distribución de las poblaciones a ambientes térmicamente favorables (^{Bickford et al., 2010}). En este escenario, es posible que algunas poblaciones de G. concinnatus puedan desplazarse y colonizar sitios de mayor altitud, principalmente hacia las zonas cercanas a la cordillera de los Andes (Fig. 3). El segundo es mediante ajustes conductuales, que involucraría la modificación del horario de actividad de los individuos (optar por una estrategia catemeral o nocturna) según la calidad térmica del hábitat (^{Arenas-Moreno et al., 2018}; ^{Diele-Viegas et al., 2018}). El tercero es mediante plasticidad fisiológica, en caso de que los individuos de G. concinnatus puedan aumentar su tolerancia al calor como respuesta al aumento de la temperatura ambiental, como puede ocurrir en otros ectotermos (^{Palumbi et al., 2014}; ^{Pörtner et al., 2009}). De continuar con la tendencia actual de CC (^{IPCC, 2018}), la pérdida y modificación del hábitat en la selva amazónica (^{Coe et al., 2013}), y de no existir mecanismos de respuesta adecuados, G. concinnatus se podría sumar a las especies extintas en el presente siglo debido al cambio climático acelerado. Cabe mencionar que de acuerdo con la Lista Roja de Reptiles del Ecuador (^{Carrillo et al., 2005}), G. concinnatus es una especie considerada como de preocupación menor (LC); por ello sugerimos reevaluar su estado de conservación, al igual que el de otros ectotermos tropicales.

En conclusión, la estrategia de termorregulación de los individuos de G. concinnatus tiende hacia el termoconformismo, es decir, depende de la temperatura y la calidad térmica del hábitat. De acuerdo con nuestras proyecciones, G. concinnatus será una especie con alto riesgo de extinción para los años 2050 y 2070. El riesgo de extinción de G. concinnatus será mayor en las poblaciones que habitan a baja altitud, y para el año 2050 es posible que la única población que persista sea la del bosque protector del Oglán Alto. Ante la susceptibilidad de los ectotermos al CC, es urgente mejorar las estrategias de mitigación y conservación para hacer frente a la pérdida de biodiversidad, especialmente en zonas megadiversas.