Introducción
El tomate (Solanum lycopersicum) es una de las principales hortalizas en el mundo por su demanda en el mercado y las grandes ganancias económicas que genera su comercialización (FAOSTAT, 2016). Esta especie es una planta arbustiva que, en su forma silvestre, puede tener un ciclo de vida superior a un año; sin embargo, con la ayuda de la ingeniería genética y procesos de selección, se han desarrollado cultivares que, con una infraestructura adecuada para controlar las condiciones ambientales de cultivo - así como el uso de insumos como fertilizantes, insecticidas, herbicidas- permiten obtener mejores rendimientos en poco tiempo.
A pesar de que la especie es ampliamente cultivada en el mundo, su diversidad genética se considera restringida (Peralta, Knapp & Spooner 2005; Govindaraj, Vetriventhan, & Srinivasan, 2015; Rao & Hodgkin, 2002). Como consecuencia, son más susceptibles a plagas, enfermedades, estrés abiótico, así como a la pérdida de sus atributos sensoriales (Bolger, Lohse, & Usadel, 2014). Este rasgo hace que las poblaciones silvestres de tomate representen una fuente de genes que pueden conferir resistencia o tolerancia a plagas y enfermedades en las variedades actualmente cultivadas (Álvarez-Hernández, Cortez-Madrigal, & García-Ruiz, 2009; Bonilla-Barrientos et al., 2014).
Existe una forma silvestre S. l. var ceraciforme que es de particular importancia, ya que se considera el ancestro de los cultivares que actualmente existen en el mundo (Sánchez-Peña et al., 2006; Peralta & Spooner, 2007). Esta se distribuye desde América del Sur en países - como Perú, Ecuador, Chile e Islas Galápagos - (Peralta & Spooner, 2007; Nakazato, Warren, & Moyle, 2010) y continúa por toda América Central hasta México (Ranc, Muños, Santoni, & Causse, 2008; Flores-Hernández et al., 2017). Esta zona se ha considerado el centro de domesticación, particularmente en la vertiente del Golfo de México en los estados de Veracruz y Puebla (Long, 1995). En este proceso de dispersión desde el sur de América hasta México (Nakazato et al., 2010; Jenkins, 1948), la especie pudo establecerse en áreas con condiciones climáticas similares a las de su origen, pero también pudo adaptarse a condiciones climáticas distintas (Croteau, 2010), lo que implica que la especie tuvo que modificar sus requerimientos fisiológicos y características morfológicas en respuesta a las nuevas condiciones ambientales (Nakazato et al., 2010; Khimoun et al., 2013).
Para comparar dos ambientes en donde se ha desarrollado una especie, como lo es el centro de origen y el de domesticación de S. l. var ceraciforme, Broennimann et al. (2012) propusieron un análisis multivariado en donde se determina si las variables ambientales de dos regiones en particular son diferentes o similares entre sí, con base en un índice de similitud de Schoener (1968). Mediante la propiedad intrínseca de los individuos de una especie de adaptarse a diferentes condiciones ambientales, la cual está ligada a la variabilidad geográfica y ambiental en la que se desarrolla (Hoyt, 1992; Rao & Hodgkin, 2002; Zuriaga et al., 2009), es posible inferir si existe una diversificación de los patrones de distribución ambiental en la región de domesticación con respecto a su centro de origen del tomate.
Existen poblaciones de tomate en el área de domesticación cuyas condiciones ambientales podrían ser distintas a las de su centro de origen. El centro de Veracruz, por ejemplo, es de particular importancia, ya que allí se encuentran poblaciones en una gran variedad de ambientes y ecosistemas, lo que sugiere que varias se han adaptado y potencialmente representarían material fundamental como banco de germoplasma. En el presente trabajo, se caracterizan y comparan las condiciones ambientales de ambas regiones (origen y domesticación) para identificar los patrones de diversificación ambiental, así como determinar las áreas de distribución geográfica potencial de la especie en donde se podrían establecer medidas de conservación y rescate del tomate S. l. var ceraciforme.
Materiales y Métodos
Registros de presencia de S. l. var ceraciforme: Para identificar las áreas de distribución geográfica potencial del centro de origen y de domesticación de S. l. var ceraciforme, se generaron modelos de nicho ecológico (MNE). Para generar estos modelos se requiere de dos fuentes de información: los registros de la especie y coberturas o variables ambientales. Por ello, se recopilaron registros de presencia de S. l. var. cerasiforme en la bibliografía especializada (Arce, 2015; Álvarez-Hernández et al., 2009), la base de datos de la Infraestructura Global de Información sobre Biodiversidad (GBIF, 2018) y el trabajo de campo (como se describe más adelante).
Posteriormente, se realizó una base de datos con todos los registros disponibles: se eliminaron los registros duplicados que carecieran de precisión en la localidad de registro o que cayeran dentro de una misma celda o pixel (~ 1 km2) equivalente al tamaño (resolución) de las variables utilizadas para generar los modelos (como se detalla posteriormente). Esto se realizó con la ayuda de los Sistemas de Información Geográfica (SIG). Mediante este proceso, se obtuvo un total de 98 registros de S. l. var. cerasiforme para el centro de origen y 67 para el área de domesticación. En ambos casos se utilizó el 80 % de los puntos para generar los modelos; el otro 20 % se utilizó para validarlo con base en la Roc Parcial como se describe posteriormente.
Muestreo de áreas seleccionadas: El trabajo de campo se realizó en el norte, centro y sur del estado de Veracruz, México de septiembre 2016 a noviembre 2017 (Tabla 1, Fig. 1A, Fig. 1B). Considerando que S. l. var. cerasiforme se distribuye desde el nivel del mar hasta los 3 300 msnm (Álvarez-Hernández et al., 2009), los recorridos en campo se realizaron entre dicho rango altitudinal. En cada localidad, se registró la ubicación geográfica y topográfica con un GPS marca Garmin y se identificó el tipo de vegetación dominante en el área con la finalidad de reportar la vegetación asociada al tomate. Además, se obtuvo el perfil climático de temperatura y precipitación (con base en las variables ambientales que se describen a continuación) para poder caracterizar ambientalmente y de manera puntual los registros obtenidos en campo.
ID | Latitud | Longitud | Localidad | Ind | Región | Fecha | Descripción del sitio |
1 | 18.810362 | -97.277332 | Maltrata | 1 | Norte | 05/06/17 | Manglar, sabana y pastizal con algunas especies Crescentia cujete, Curatella americana, Byrsonima crassifolia, Sabal mexicana, Quercus oleoides dominantes; actividades agropecuarias: cultivo de naranja y pastos para ganadería. Altitud media (m.s.n.m.) = 18 (Mín: 14, Máx: 24) Temperatura promedio anual (°C) = 24 (Mín: 15, Máx. 33) Precipitación promedio anual (mm) = 1 166 |
2 | 20.871069 | -97.344062 | Países Bajos | 1 | Norte | 14/05/17 | - |
3 | 20.972661 | -97.527579 | Ojite | 4 | Norte | 14/05/17 | - |
4 | 20.951682 | -97.460643 | Tuxpan | 2 | Norte | 15/05/17 | - |
5 | 20.942754 | -97.453605 | Retoño | 1 | Norte | 15/05/17 | - |
6 | 19.255000 | -97.091944 | Tlanalapan | 3 | Centro | 01/11/17 | Bosques de encinos, coníferas y bosques húmedos de altas montañas, especies como: Enterolobium cyclocarpum, Tabebuia pentaphylla, Holiocereus sp. Cultivos de café, plátano, limón, naranja, chayote. Altitud media (m.s.n.m.) = 1 346 (Mín: 553, Máx: 2 100) Temperatura promedio anual (°C) = 19 (Mín: 5, Máx: 32) Precipitación promedio anual (mm): 1 791 |
7 | 19.076197 | -96.894200 | Tepatlaxco | 12 | Centro | 03/11/17 | - |
8 | 19.044236 | -96.908142 | Ocotitlán | 1 | Centro | 03/11/17 | - |
9 | 19.047061 | -96.972935 | Ixhuatlán del Café | 8 | Centro | 03/11/17 | - |
10 | 19.052289 | -96.935747 | Ocotitlán | 6 | Centro | 03/11/17 | - |
11 | 18.416775 | -95.102441 | Catemaco | 1 | Centro | 26/12/16 | - |
12 | 18.863122 | -96.901091 | Peñuela | 2 | Centro | 30/12/16 | - |
13 | 18.888717 | -96.936399 | Córdoba | 2 | Centro | 30/12/16 | - |
14 | 19.413464 | -96.915741 | Mahuixtlán | 3 | Centro | 02/04/17 | - |
15 | 19.410983 | -96.919018 | Enríquez | 2 | Centro | 03/04/17 | - |
16 | 19.042595 | -97.069957 | Coscomatepec | 3 | Centro | 08/11/16 | - |
17 | 19.566523 | -96.983029 | Tlalnelhuacoyan | 5 | Centro | 05/04/17 | - |
18 | 19.451676 | -96.969638 | Coatepec | 3 | Centro | 01/04/17 | - |
19 | 18.816502 | -97.225010 | Balastrera | 2 | Centro | 05/06/17 | - |
20 | 19.164444 | -97.081666 | Huatusco | 3 | Centro | 24/09/16 | - |
21 | 19.156944 | -97.081117 | Chavaxtla | 1 | Centro | 24/09/16 | - |
22 | 18.889200 | -96.806700 | Providencia | 3 | Centro | 15/01/17 | - |
23 | 18.744086 | -97.090591 | Tequila | 2 | Centro | 28/01/17 | - |
24 | 18.808572 | -97.092217 | Tlilapan | 2 | Sur | 10/11/17 | Selvas altas perennifolias con especies como Bursera simaruba, Cupania dentata, Ficus sp., Dentropanax arboreus. Cultivos de maíz, frijol, piña y pastos para ganadería. Altitud media (m.s.n.m.) = 392 (Mín: 138, Máx: 541) Temperatura promedio anual (°C) = 23 (Mín: 15, Máx: 34) Precipitación promedio anual (mm) = 2 260 |
25 | 18.384573 | -95.204378 | Salto de Eyipantla | 5 | Sur | 29/12/17 | - |
26 | 18.469031 | -95.215091 | Tonalapan | 1 | Sur | 27/12/17 | - |
27 | 18.468964 | -95.217206 | Tonalapan (1) | 1 | Sur | 27/12/17 | - |
Las variables ambientales se obtuvieron del proyecto Worldclim, el cual proporciona valores promedio para 1950-2000 (Hijmans, Cameron, Parra, Jones, & Jarvis, 2005) con una resolución espacial de 1 pixel = 0.0083 º (~ 1 km2). Para elegir las variables que se utilizaron para correr los modelos, primero se evaluó la multicolinealidad entre las 19 variables con base en un análisis de coeficiente de correlación (r) (Pearson & Dawson, 2003). Posteriormente, se determinó como umbral de corte un valor de correlación de r ≥ 0.80. Finalmente, se eliminaron las variables que estaban altamente correlacionadas (variables con un valor de r ≥ 0.80; Yang, Kushwaha, Saran, Xu, & Roy, 2013).
Este proceso se realizó con la finalidad de evitar la autocorrelación espacial entre las variables (Escobar, Qiao, Phelps, Wagner, & Larkin, 2016), lo cual puede provocar un sobreajuste en la construcción del modelo (Lozier, Aniello, & Hickerson, 2009). Las ocho variables que se utilizaron para correr los modelos fueron: bio1 (temperatura media anual, °C), bio2 (rango medio diurno de temperatura, °C), bio3 (isotermalidad, °C), bio4 (índice de variación estacional de temperatura), bio12 (precipitación media anual, mm), bio14 (precipitación del mes más seco, mm), bio15 (índice de variación estacional de la precipitación) y bio18 (precipitación del trimestre más seco, mm).
Construcción de los modelos: Los modelos se generaron con el algoritmo Maxent versión 3.4.1 que se basa en el principio de máxima entropía para determinar la idoneidad de las condiciones ambientales para la especie (Phillips, Anderson, & Schapire, 2006). Este algoritmo ha sido ampliamente utilizado y ha demostrado un buen desempeño (Elith et al., 2006). Maxent se configuró por defecto (Phillips & Dudik, 2008), se desactivaron los módulos de “extrapolate” y “do Clamping” y se utilizó un formato de salida logístico para el modelo. Los parámetros establecidos fueron los siguientes: un valor máximo de 1 000 iteraciones, un límite de convergencia de 0.00001 y el valor del factor de regularización en 1.
Seguidamente, se generaron dos MNE y sus correspondientes proyecciones, es decir, se corrió el MNE en México que representa el área de domesticación y se proyectó en Suramérica (que representa el centro de origen). Posteriormente, se generó un MNE en Suramérica (origen) y se proyectó a México (domesticación) con la finalidad de identificar la superposición geográfica de las condiciones ambientales en ambas regiones. Para correr los modelos, se delimitó el área de accesibilidad o “M” (sensuSoberón & Peterson, 2005) con base en las provincias biogeográficas definidas por Morrone (2014) y Löwenberg-Neto (2014). Específicamente, en el centro de origen se incluyó la región Neotropical formada por los países de Panamá, Venezuela, Colombia, Ecuador, Perú, Bolivia, Chile, Argentina Paraguay y Brasil. Para el caso del área de domesticación, la “M” se delimitó con base en las siguientes provincias biogeográficas: Xerófila mexicana, Pacífica, Mesoamericana de Montaña, Caribe y Guajira, desde México hasta Nicaragua.
Para evaluar los modelos, se utilizó el área bajo la curva (Area Under Curve o AUC) de la característica operativa del receptor (Receiver Operating Characteristic o ROC), implementada directamente en Maxent. Asimismo, para complementar la evaluación de los modelos, se usó el método de ROC parcial (Peterson, Papes, & Soberón, 2008) en el software propuesto por Barve (2008). El resultado probabilístico (0 a 1) de cada MNE obtenido se convirtió en un mapa binario de presencia-ausencia, cuyo valor de corte fue del 5 % de presencias de entrenamiento. Por último, se realizó un análisis de MOP (Mobility Operating Parity; Owens et al., 2013) con la finalidad de calcular la analogía climática que podría existir entre la “M” que se delimitó para el MNE del centro de origen y el espacio geográfico donde se proyectó, correspondiente al área de domesticación.
Comparación ambiental entre el centro de origen y el área de domesticación: Para comparar las condiciones ambientales entre el centro de origen y el de domesticación, se utilizaron las ocho variables con que se corrieron los MNE del tamaño (extensión) del área de accesibilidad o “M” previamente descritas. Posteriormente, se realizó la prueba propuesta por Broennimann et al. (2012), modificada por Di Cola et al. (2017). Esta propuesta se basa en un análisis multivariado (componentes principales) e incluye una primera prueba que es de equivalencia de nicho. En este caso particular, se evaluó si las condiciones ambientales en el centro de origen eran equivalentes (hipótesis nula) a las del área de domesticación o si eran menos equivalentes que lo esperado por azar (hipótesis alternativa, P < 0.05).
La segunda prueba realizada fue la de similitud de nicho, donde la hipótesis nula indica que los nichos son similares, mientras que la hipótesis alternativa indica que los nichos son menos similares que lo esperado por azar (P < 0.05). En ambas pruebas, el azar o modelos nulos están representados por 100 repeticiones. Por ello, utilizamos la opción “lower” para determinar la hipótesis alternativa (sensuBroennimann et al., 2012, Di Cola et al., 2017). Finalmente, las variables climáticas que se utilizaron para correr los MNE del centro de origen y de domesticación fueron sometidas a un análisis de varianza y una prueba de comparación de medias (Tukey; P ≤ 0.05) para detectar si existían diferencias significativas entre ambos conjuntos de variables.
Resultados
El MNE en el centro de origen tiene una extensión geográfica de 5 217 041 km2 y se distribuye en la vertiente del Pacífico en Suramérica desde Panamá hasta el norte de Chile, Argentina, Paraguay y el oeste de Brasil. Los valores de idoneidad más altos se ubicaron principalmente en Perú (Anexo 1). En el caso del área de domesticación, esta tuvo una extensión de 1 084 894 km2 que se distribuye desde la vertiente del Pacífico y el Golfo de México hasta Honduras. Los valores de idoneidad más altos se ubicaron principalmente en la vertiente del Golfo de México (Anexo 2).
Por su parte, la proyección del área de origen al de domesticación tuvo una extensión de 1 002 983 km2. De este total, las condiciones ambientales en 511 618 km2 fueron parecidas a las del centro de origen, mientras que en el resto (491 365 km2) las condiciones ambientales fueron diferentes. Por último, la proyección del área de domesticación a la de origen tuvo una extensión de 1 481 987 km2. De este total, las condiciones ambientales fueron parecidas en 392 529 km2 a las del área de domesticación; en el resto (1 089 458 km2), las condiciones ambientales resultaron diferentes (Fig. 1A, Fig. 1B). Con base en las 27 localidades con presencia de S. l. var cerasiforme en el estado de Veracruz, se obtuvieron ocho registros en el norte, 60 en el centro y ocho en el sur, en diferentes tipos de vegetación y un amplio gradiente altitudinal y de variables ambientales (Tabla 1).
Los valores de AUC del MNE de origen y de domesticación fueron 0.93 y 0.80 respectivamente para los datos de entrenamiento. Los valores de ROC-parcial fueron 1.86 y 1.71 respectivamente y, en ambos casos, fueron significativos (P = 0.0001). Esto implica que el modelo al azar fue un buen modelo predictivo. Por otro lado, el resultado del análisis de MOP reveló que el área geográfica con climas no análogos entre las condiciones climáticas entre el centro de origen y área de domesticación fue mínima (Anexo 3). Por lo tanto, las proyecciones del MNE de origen al MNE de domesticación fueron adecuadas.
Con respecto a la comparación ambiental entre el centro de origen y de domesticación, el análisis de componentes principales mostró que el primer componente explicó el 57.42 % de la variación total y las variables de temperatura (bio1, bio2, bio3, bio4) y de precipitación (bio14) fueron las más importantes. El segundo componente explicó solamente el 17.36 % de la variación y el índice de variación estacional de la precipitación (bio15) fue lo más importante (Fig. 2A). También observamos que las condiciones ambientales del centro de origen presentan una mayor amplitud con respecto al área de domesticación (Fig. 2B, Fig. 2C, Fig. 2D). La hipótesis nula de equivalencia fue rechazada (P = 0.003, Fig. 2E), es decir, los modelos de nicho ecológico entre ambas regiones fueron menos equivalentes que lo esperado por azar. Por otro lado, en el caso de la prueba de similitud, no se rechazó la hipótesis nula (P = 0.683; Fig. 2F). Esto nos indica que la similitud que existe entre ambas regiones también se puede obtener por azar.
Al comparar las variables del centro de origen contra el área de domesticación, se encontraron diferencias significativas (P ≤ 0.05) en tres variables de temperatura (bio2, bio3, bio4) y tres variables de precipitación (bio12, bio14, bio15), pero no se detectaron diferencias significativas al comparar temperatura media anual (bio1) y la precipitación del trimestre más seco (bio18; Tabla 2). Cabe señalar que en el centro de origen llueve tres veces más (bio14) que en su área de domesticación, mientras que el índice de variación de temperatura (bio4) y de precipitación (bio15) en esta zona es menor en comparación con el área de domesticación. Esto indica que el clima en el centro de origen es más uniforme a lo largo del año, mientras que en el área de domesticación está sujeto a cierta variación ambiental.
- | bio1 | bio2 | bio3 | bio4 | bio12 | bio14 | bio15 | bio18 |
Origen | 21.28a | 10.873b | 79.58a | 754.21b | 1768.73a | 74.01a | 49.35b | 445.31a |
Domesticación | 22.24a | 13.326a | 61.73b | 2644.05a | 1275.28b | 20.75b | 87.85a | 485.22a |
Las medias con letras diferentes indican diferencias significativas (P < 0.05). bio1: temperatura media anual (°C), bio2: rango medio diurno de temperatura (°C), bio3: isotermalidad, bio4: índice de variación estacional de temperatura, bio12: precipitación media anual (mm), bio14: precipitación del mes más seco (mm), bio15: índice de variación estacional de la precipitación, bio18: precipitación del trimestre más seco (mm). Means with different letters indicate significant differences (P < 0.05). bio1: annual mean temperature °C, bio2: mean diurnal range temperature °C, bio3: isothermality, bio4: temperature seasonality, bio12: annual precipitation, bio14: precipitation of driest month (mm), bio15: precipitation seasonality, and bio18: precipitation of warmest quarter (mm).
Discusión
Los MNE y sus correspondientes transferencias entre el centro de origen y de domesticación de S. l. var cerasiforme mostraron áreas geográficas donde se superponen ambas regiones, por ejemplo, el centro de México. Esto quiere decir que existen poblaciones de tomate cuyos requerimientos ambientales son muy parecidos tanto en Suramérica como en México. Sin embargo, también se observaron áreas geográficas en México (domesticación) cuyas condiciones ambientales no fueron equivalentes a las de su centro de origen, específicamente el noroeste de México y la vertiente del Golfo de México, como se demostró en la prueba de comparación de nichos (Broennimann et al., 2012; Di Cola et al., 2017).
Cabe señalar que se tiene registro histórico de presencia de S. l. var cerasiforme en ambos casos, incluso en la vertiente del Golfo de México se corroboró su presencia actual durante el trabajo de campo, lo que sugiere su viabilidad. El hecho de encontrar poblaciones de tomate en el área de domesticación con distintas condiciones ambientales a las de su origen podría sugerir que existe un proceso de adaptación y variación genética para sobrevivir en entornos diferentes a los de su origen, particularmente si consideramos que el ambiente es uno de los factores que provocan cambios que se pueden manifestar a nivel morfológico y genético de una especie (Schluter, 2000). Esta relación entre las variables climáticas y diferencias genotípicas en poblaciones de otras especies de tomates ya ha sido reportada por otros autores (Zuriaga et al., 2009; Sifres. Blanca, & Nuez, 2011).
Con respecto a las variables ambientales de temperatura, Vargas et al. (2005) registraron plantas de tomate a una temperatura media anual entre 18 y 22 ºC en el área de domesticación; mientras que Nakazato, Bogonovich, & Moyle (2008) reportaron una temperatura media anual de 20.8 ºC en Suramérica. Los valores de temperatura media anual obtenidos en este trabajo (21.28 ºC en su centro de origen y 22.24 ºC en el área de domesticación) no presentaron diferencias significativas entre ambas regiones y coinciden con los resultados reportados previamente por otros autores. Esto sugiere que en el área de domesticación existen las condiciones de temperatura media anual para la distribución y desarrollo del tomate.
Por otro lado, en las variables de temperatura donde sí se encontraron diferencias significativas, el índice de variación estacional fue mayor para el área de domesticación. Esto implica una mayor fluctuación de temperatura a largo del año en comparación con el centro de origen. Aunado a lo anterior, el rango medio diurno de temperatura, que hace referencia a una fluctuación de temperatura entre el día y la noche, mostró que para el área de domesticación (13 °C) fue 2 °C superior a la del centro de origen (10.8 °C). Esta variación de temperatura a diferentes escalas (estacional y durante el día y la noche) es importantes ya que esto afecta la floración, fructificación y crecimiento de la planta (Chaves-Barrantes & Gutiérrez-Soto, 2017). Asimismo, la temperatura controla la germinación de la semilla del tomate (Labouriau & Osborn, 1984), activando las reacciones enzimáticas que toman lugar durante el proceso de germinación de semillas de la planta (Liu et al., 2013).
Durante el trabajo de campo, pudimos observar algunas características como plantas con abundante follaje de color verde intenso, decoloración en los frutos- por ejemplo, frutos de color rojo con poca intensidad -, así como plantas con crecimiento restringido particularmente del tallo. Estas características podrían responder a la variación de la temperatura como ya lo han sugerido otros autores (Yang, Li, Li, Zhang & Thomas, 2016; San Martín-Hernández, Ordaz-Chaparro, Sánchez-García, Beryl Colinas-Leon, & Borges-Gómez, 2012). Sin embargo, estas características corresponden a observaciones generales de los ejemplares en campo, por lo que es necesario realizar un análisis formal y evaluar las características morfológicas tanto del centro de domesticación como del área de domesticación.
En el caso de las variables de precipitación, algunos reportes de precipitación media anual para el área de domesticación, como el de Vargas et al. (2005), indican que oscila entre 800 y 1 000 mm. Mientras tanto, Álvarez-Hernández et al. (2009) muestra una oscilación similar entre 751 y 1 014 mm. No obstante, en el presente trabajo la precipitación media anual obtenida fue de 1 275 mm. Todos estos valores fueron menores a los reportados para el centro de origen: Nakazato et al. (2010) reporta 1 800 mm, por ejemplo. Es importante aclarar que estos datos contabilizaron la precipitación del mes más seco, ya que en el área de domesticación llueve menos (20.75 mm) con respecto al centro de origen (74.01 mm).
Así también, el índice de variación estacional es significativamente menor en el centro de origen en comparación con el área de domesticación. Esto quiere decir que, además de existir una mayor precipitación en el en centro de origen, la lluvia es constante a lo largo del año. Estas diferencias son muy importantes ya que la lluvia interviene en el transporte de nutrientes, fotosíntesis, mantenimiento de la turgencia e hidratación de las células de la planta (Basurto-Sotelo, Núñez, Pérez, & Hernández, 2012), lo cual podría ser un indicador de diferenciación y adaptación de las poblaciones de S. l. var cerasiforme a nuevas condiciones de precipitación.
Por otro lado, los registros de S. l. var. Cerasiforme, específicamente en el estado de Veracruz, muestran cómo las poblaciones se encuentran en una variedad importe de tipos de vegetación desde zonas de manglar a baja altitud (14 m.s.n.m.) hasta bosques de encinos y coníferas (2 100 m.s.n.m.); la población pasa por áreas con cierto grado de perturbación como zonas agrícolas y ganaderas. También se obtuvieron registros en sitios con una variabilidad climática importante, por ejemplo, entre 1 166 y 2 260 mm de precipitación promedio anual y en intervalos de temperatura entre 5 y 34 °C. Este último dato es de particular importancia porque, bajo condiciones controladas, se sabe que una temperatura inferior a 12 °C limita el crecimiento y desarrollo de la planta (Hu, Zhou, Du, Xia, & Yu, 2006). La presencia de poblaciones de S. l. var. cerasiforme en este amplio gradiente de condiciones, incluso bajo su límite de tolerancia como es el caso de la temperatura, pone en evidencia su plasticidad y amplia capacidad de adaptación como ya ha sido reportado en otros estudios (Comlekcioglu, Simsek, Boncuk, & Aka-Kacar, 2010; Kamel, Soliman, Mandour, & Ahmed, 2010; Flores-Hernández et al., 2017).
Es importante aclarar que es posible que, de los ejemplares colectados durante el trabajo de campo en Veracruz, haya plantas domesticadas o semi-domesticadas; sin embargo; no nos es posible determinar su grado de domesticación en el presente estudio. No obstante, consideramos que la comparación entre el centro de origen en Suramérica y el de domesticación en México es adecuado ya que se comparan las condiciones ambientales de ambas regiones y los lugares donde actualmente hay poblaciones de S. l. var. cerasiforme.
Finalmente, el centro de origen y de domesticación de S. l. var. cerasiforme tienen características ambientales en común a pesar de la distancia geográfica; no obstante, existen áreas en la región de domesticación con condiciones ambientales diferentes a las de su centro de origen, particularmente en la vertiente del Golfo de México en Veracruz. Los registros obtenidos en esta región se encuentran en un gradiente ambiental muy amplio que sugiere una plasticidad y amplia capacidad de adaptación del tomate a estos ambientes. Consideramos que las localidades de registro en esta región requieren de medidas de manejo y conservación, ya que podrían representar un material fundamental para los bancos de germoplasma.