Introducción

Las musarañas (Mammalia, Soricidae) son los mamíferos de menor tamaño en el mundo. Debido a sus demandas energéticas, son depredadores voraces de artrópodos y anélidos, pero también son presas comunes de aves rapaces, serpientes y mamíferos carnívoros (^{Churchfield, 1990}). Con cerca de 440 especies agrupadas en 26 géneros, representan a una de las familias de mamíferos con mayor distribución geográfica en el planeta, habitando todos los ecosistemas terrestres (^{Burgin et al., 2018}). Su capacidad de dispersión relativamente limitada, tiempo de generación corto y asociaciones estrechas con hábitats particulares, han contribuido a que las musarañas reciban una atención especial para evaluar el papel que han jugado el cambio climático del pasado y la orografía en la conformación de los patrones de variabilidad genética actual y los procesos de divergencia de linajes. Lo anterior ha sido particularmente importante en regiones templadas, boreales y paleotropicales (^{Giarla y Esselstyn, 2015}; ^{Hope et al., 2020}; ^{Jacquet et al., 2013}; ^{Vega et al., 2020}; ^{Vuilleumier y Fontanillas, 2007}). En el Neotrópico, la situación es diferente. El acervo relativamente limitado tanto de ejemplares en colecciones biológicas como de muestras para análisis genéticos, así como el número reducido de taxónomos, han impedido un avance sustancial en el conocimiento sobre la evolución y biogeografía del grupo (^{Lorenzo et al., 2019}; ^{Noguera-Urbano et al., 2019}).

Por su diversidad y distribución, las musarañas de orejas cortas (género Cryptotis) son el representante principal de musarañas neotropicales y el único género que habita en Norte y Sudamérica. Con pocas excepciones, las especies están restringidas a ecosistemas de media y alta montaña, en donde se combinan condiciones de humedad y precipitación alta con temperatura baja (^{Guevara et al., 2015}; ^{Noguera-Urbano et al., 2019}; ^{Woodman et al., 2012}). Dentro del género, el grupo de especies Cryptotis mexicanus comprende a musarañas exclusivas del norte del Neotrópico y estrechamente asociadas al bosque húmedo de montaña (o simplemente bosque de niebla), entre los 500 y los 3,000 m snm (^{Guevara et al., 2015}). Esta agrupación de especies se caracteriza por poseer extremidades anteriores grandes que han sido interpretadas como adaptaciones para explotar hábitos semifosoriales (^{Guevara, 2017}; ^{Woodman, 2011}).

En el grupo C. mexicanus, las relaciones filogenéticas muestran a un clado distribuido en el extremo boreal compuesto por 3 especies: C. mexicanus (Coues, 1877) y C. obscurus (Merriam, 1895) cercanamente relacionadas y hermanas de C. nelsoni (Merriam, 1895) (^{Guevara y Cervantes, 2014}; ^{He et al., 2015}). Cryptotis obscurus se distribuye en las provincias biogeográficas de Sierra Madre Oriental y el norte de la Veracruzana; C. mexicanus al este de la Faja Volcánica Transmexicana, noreste de la Sierra Madre del Sur y en las Tierras Altas de Chiapas; mientras que C. nelsoni es endémica de la Sierra de Los Tuxtlas, al este de la provincia Veracruzana (fig. 1). Dada su afinidad con el bosque de niebla de estas provincias, se ha especulado que la evolución de estas musarañas está sincronizada con la intrincada historia de expansión, conectividad y fragmentación de dicho ecosistema a través de los últimos ciclos glaciales-interglaciales del Cuaternario (^{Choate, 1970}; ^{Guevara y Sánchez-Cordero, 2018a}).

Figura 1 Registros conocidos del clado compuesto por C. mexicanus (círculos), C. obscurus (triángulos) y C. nelsoni (cuadrados) en México. Los símbolos rojos indican las muestras genéticas analizadas en el presente estudio. Las colores más oscuros representan mayor elevación. Se señalan las características fisiográficas principales que son mencionadas en el texto. 

Las glaciaciones, por un lado, habrían ocasionado el descenso en elevación y la expansión de los bosques de niebla (^{Caballero et al., 2010}), por lo que las tierras bajas servirían como puentes para permitir la conectividad de poblaciones; mientras que durante los interglaciales (como en el presente), estas zonas bajas fungen como barreras que aíslan a las poblaciones al no poseer condiciones idóneas para el desarrollo del bosque de niebla (^{Ramírez-Barahona y Eguiarte, 2013}). En específico, los valles y cañones que conducen a los ríos Tampaón, Moctezuma, Amajac al norte de la Sierra Madre Oriental, a los ríos Blanco y Santo Domingo al norte de la Sierra Madre del Sur, la planicie hacia la costa del Golfo de México y el Istmo de Tehuantepec, podrían actuar como barreras y corredores que aíslan y conectan a las poblaciones de estas especies durante las oscilaciones climáticas del Cuaternario (^{Guevara y Sánchez-Cordero, 2018a}).

Una revisión taxonómica reciente no encontró una diferenciación morfológica ni ecológica significativa entre C. obscurus y C. mexicanus, lo cual podría ser consecuencia de un proceso de divergencia reciente o una conectividad frecuente entre las poblaciones durante los ciclos glaciales-interglaciales. Esto ha dejado dudas sobre el estado taxonómico de ambas especies (^{Guevara y Sánchez-Cordero, 2018a}). Por el otro lado, la diferenciación morfológica notable de C. nelsoni podría deberse a su distribución restringida en el sistema montañoso de la Sierra de Los Tuxtlas, donde la especie se encuentra aislada y probablemente sin flujo con el resto de las especies por un tiempo largo (^{Guevara y Sánchez-Cordero, 2018b}; ^{Woodman y Timm, 1999}), incluso durante la expansión del bosque de niebla en los eventos glaciales más recientes (^{Guevara, 2020}).

Para poner a prueba si los valles y cañones mencionados anteriormente podrían estar promoviendo la divergencia dentro del clado C. mexicanus, C. obscurus y C. nelsoni, sería crucial contar con información genética. Afortunadamente, en los últimos años, el trabajo de campo dirigido a la recolecta de musarañas ha permitido incrementar la representación geográfica y el almacenamiento de muestras para estudios genéticos (^{Guevara y Sánchez-Cordero, 2018a}, ^b; ^{Mayén-Zaragoza et al., 2019}). Por lo anterior, en este estudio usamos secuencias del gen mitocondrial Citocromo b para evaluar la estructura genética y la diferenciación dentro del clado conformado por estas 3 especies. Esperamos que este primer estudio filogeográfico de musarañas neotropicales no solo contribuya a comprender la historia biogeográfica de este clado y otras especies del bosque de niebla en la norte del Neotrópico, sino que repercuta en el conocimiento taxonómico del grupo y su conservación.

Materiales y métodos

Los ejemplares de referencia se encuentran depositados en las siguientes colecciones científicas de México (apéndice): Colección Nacional de Mamíferos (CNMA), Ciudad de México; Museo de Zoología ‘Alfonso L. Herrera’ (MZFC), Ciudad de México; y Colección Mastozoológica de El Colegio de la Frontera Sur (ECO- SC-M), San Cristóbal de las Casas. Se extrajo DNA desde tejido con un kit comercial DNeasy® Tissue Qiagen, y la amplificación del gen mitocondrial Citocromo b (Citb) por medio de la Reacción en Cadena de la Polimerasa (PCR, por su siglas en inglés) usando el par de oligonucleótidos MVZ05 (5′-CGAAGCTTGATATGAAAAACCATCGTTG; ^{Smith y Patton, 1999}) y H15915 (5′-AACTGCAGTCATCTCCGGTTTACAAGAC; ^{Irwin et al., 1991}). El volumen final de reacción fue de 15 µl utilizando el buffer de la enzima en concentración 1X, 0.16 µM del cebador MVZ05, 0.16 µM del cebador H15915, 0.4 mM de cada dNTP, 0.6 U de Taq Polimerasa MyTaq^TM Bioline y una concentración entre 50-265 ng/µl de DNA.

El programa de la PCR incluyó 3 min de desnaturalización inicial a 95 °C; seguidos de 35 ciclos que incluyeron 30 segundos de desnaturalización a 95 °C, un min de alineamiento a 49 °C, y 2 min de extensión a 72 °C. Se visualizó el producto de PCR utilizando 3 µl mediante electroforesis en geles de agarosa al 1.5%. Se utilizó el kit de purificación Montage PCR de Millipore y el kit de secuenciación ABI PRISM Big Dye Terminator en un secuenciador ABI 3730XL (Applied Biosystems, Foster City, CA, USA). Se obtuvieron 46 secuencias genéticas del clado conformado por C. mexicanus (29), C. obscurus (7) y C. nelsoni (10). Se realizó la edición y limpieza de las secuencias por medio de una inspección visual con el programa MEGA X (^{Kumar et al., 2018}). Además, se obtuvieron 24 secuencias de la base de datos GenBank (C. mexicanus [9], C. obscurus [2] y C. nelsoni [5] y 8 utilizadas como grupos externos (apéndice).

Se identificó la distribución espacial de la variabilidad genética mediante el paquete Geneland ® v. 4.0.8 en R® v3.4.4, el cual localiza discontinuidades genéticas entre las muestras (^{Guillot et al., 2005}). Para ello, se utilizó el modelo haploide, un modelo de frecuencias alélicas no correlacionadas y un modelo espacial verdadero. También se usaron las cadenas de Markov Monte-Carlo con 6 corridas independientes de 5,000,000 iteraciones cada una, muestreando datos cada 1,000. Aplicamos un descarte ("burn-in") del 25%, comparamos las corridas independientes y seleccionamos el análisis con mayor probabilidad posterior.

Se construyó una red de haplotipos para establecer la relación entre todos los haplotipos mediante una red de unión de medianas ("median-joining") usando el programa Network v10.0 (^{Fluxus Technology Ltd., 2015}). Calculamos la diversidad haplotídica, nucleotídica, así como el número de sitios polimórficos en el programa DNAsp v6.12.02 (^{Librado y Rozas, 2009}).

Para medir la diferencia genética en el Citb entre los grupos identificados en el análisis de agrupación espacial bayesiana, se calcularon las distancias genéticas utilizando el programa MEGA X (^{Kumar et al., 2018}). Para ello, se agruparon las secuencias en los grupos identificados y se ejecutó un análisis de distancia neta entre el promedio de los grupos. Los parámetros definidos fueron 10,000 repeticiones de bootstrap mediante el modelo p-distance (^{Starcová et al., 2016}). Asimismo, se analizó la distribución de la variabilidad genética entre los grupos definidos previamente. Para ello, se realizó un análisis de varianza molecular (Amova) implementando 16,000 permutaciones para el cómputo de una matriz de distancias, bajo el modelo de proporción de diferencias, usando el programa Arlequin (^{Excoffier y Lischer, 2010}).

Se llevó a cabo un análisis de inferencia bayesiana para determinar las relaciones filogenéticas entre los grupos identificados en los análisis previos. Para estimar el mejor modelo evolutivo de sustitución nucleotídica utilizado para la inferencia bayesiana, utilizamos el software Smart model selection (SMS) en el entorno de PhyML (^{Anisimova y Gascuel, 2006}; ^{Lefort et al., 2017}) AIC. Según el criterio de información bayesiano (BIC), el mejor modelo de sustitución nucleotídica para los datos fue el modelo HKY85+G. La inferencia bayesiana se realizó en el programa Mr.Bayes v3.2.7 (^{Ronquist et al., 2012}). Se realizaron 2 corridas con 15 millones de generaciones e implementando 8 Cadenas de Markov Monte Carlo (MCMC; 7 cadenas calientes y una fría). Se muestrearon cada 500 generaciones, descartando el primer 25% de la muestra (burn-in). Finalmente, se generó el árbol de mayoría con los valores de probabilidad posterior para cada clado. Para analizar y visualizar los resultados, se usaron los programas Tracer v1.7.1 y FigTree v1.4.4 (^{Rambaut et al., 2018}).

Resultados

El análisis de agrupación espacial reveló 4 grupos genéticos dentro del clado conformado por C. mexicanus, C. obscurus y C. nelsoni (fig. 2). La probabilidad posterior de que existieran más de 4 grupos según el Citb fue extremadamente baja. El primer grupo (en adelante denominado NORTE) incluye a todas las muestras pertenecientes a C. obscurus y a los haplotipos previamente identificados como C. mexicanus del centro del país y de la Sierra de Otontepec, Veracruz, al este de la Sierra Madre Oriental (fig. 2A). El segundo grupo (CENTRO; fig. 2B) está formado por haplotipos de C. mexicanus desde los alrededores del volcán Pico de Orizaba, Veracruz hasta el norte del río Santo Domingo, Oaxaca. El tercer grupo (SUR) está compuesto por todas las muestras de C. mexicanus al sur del río Santo Domingo, incluyendo a las 2 secuencias al este del istmo de Tehuantepec, Oaxaca (fig. 2C). Finalmente, el cuarto grupo (LOS TUXTLAS) está formado exclusivamente por los haplotipos correspondientes a C. nelsoni, en el sistema montañoso aislado de Los Tuxtlas (fig. 2D). Los valores máximos de probabilidad posterior para la pertenencia poblacional de los haplotipos fueron cercanos a 50%, mientras que la probabilidad posterior de que los haplotipos de un grupo pertenecieran a cualquiera de las otras 3 agrupaciones fue considerablemente menor.

Figura 2 Estructura poblacional del clado Cryptotis mexicanus, C. obscurus y C. nelsoni: grupos NORTE (A), CENTRO (B), SUR (C) y LOS TUXTLAS (D). El gradiente de colores representa en los tonos más claros una mayor probabilidad posterior de pertenencia poblacional. 

Se obtuvieron 42 haplotipos en la muestra total (Hd: 0.9672; Pi: 0.0469). El número de sitios polimórficos fue de 211 en las 1,050 pb (20.11%). La red mostró haplogrupos similares a los resultados obtenidos por el análisis de agrupación espacial bayesiana, donde el grupo NORTE se formó con 19 secuencias distribuidas en 15 haplotipos. A su vez, el grupo CENTRO estuvo compuesto por 19 secuencias con 12 haplotipos y el SUR por igual número de secuencias y 8 haplotipos. Ambos están formados únicamente por haplotipos de C. mexicanus. En el caso de LOS TUXTLAS, se obtuvieron 7 haplotipos de C. nelsoni en las 16 secuencias genéticas examinadas.

En el grupo NORTE, se observó la fragmentación entre diferentes regiones geográficas incluidas en el estudio (fig. 3). Este grupo requirió al menos 48 pasos mutacionales respecto al grupo SUR y 62 pasos mutacionales respecto al grupo de LOS TUXTLAS. Para el caso del grupo SUR, los 2 haplotipos más distantes al resto del mismo grupo corresponden a aquellos localizados al este del istmo de Tehuantepec, requiriendo al menos 33 pasos mutacionales al nodo más próximo. Por otro lado, al menos 42 pasos mutacionales separan al grupo CENTRO del grupo SUR. En el caso del grupo de LOS TUXTLAS (C. nelsoni), los pocos haplotipos se mostraron cercanos entre sí y más relacionados con el grupo NORTE.

Figura 3 Red de haplotipos del gen Citb para el clado conformado por Cryptotis mexicanus, C. obscurus y C. nelsoni. Los colores representan a los 4 grupos encontrados en el análisis de agrupación espacial: NORTE (verde), CENTRO (amarillo), SUR (rojo) y LOS TUXTLAS (azul). El tamaño de los nodos es proporcional a las frecuencias de los haplotipos en la muestra, y los nodos intermedios turquesa representan nodos hipotéticos. Los segmentos en cada rama representan pasos mutacionales entre los haplotipos. 

El grupo genéticamente más distante del resto fue el de LOS TUXTLAS (> 5.1%). Dentro del par de especies C. mexicanus y C. obscurus, la mayor distancia resultó entre la población NORTE y CENTRO (> 4.4%), mientras que los grupos más cercanos genéticamente son el grupo CENTRO y el grupo SUR (2.4 %). El análisis de varianza molecular mostró que más de 70% de la variación genética está entre las poblaciones y no dentro de ellas. Asimismo, el análisis reveló que los valores de diferenciación entre los 4 grupos resultaron significativos (tabla 1).

Las frecuencias de cada base nitrogenada en los 1,050 pb examinados se distribuyeron bajo los siguientes porcentajes: 29.3% (T), 29.3% (C), 27.6% (A) y 13.8%

(G). El árbol consenso de la inferencia bayesiana mostró resultados altamente congruentes con los resultados de agrupación espacial bayesiana y la red de haplotipos (fig. 4). El grupo de LOS TUXTLAS es un clado bien diferenciado y apoyado (1.0), formado exclusivamente por los haplotipos de C. nelsoni. A su vez, este es grupo hermano del clado que incluye a todas las muestras previamente identificadas como C. mexicanus y C. obscurus. El grupo NORTE se conformó con todos los haplotipos de C. obscurus, además de los de C. mexicanus del centro de México y de la sierra de Otontepec, con alto valor de apoyo (1.0). Este grupo es hermano del clado que contiene a los grupos CENTRO y SUR con valor de apoyo alto (1.0); sin embargo, la monofilia del grupo SUR está débilmente apoyada (0.65).

Figura 4 Árbol de consenso bayesiano con 70 secuencias del gen Citb. Las primeras 8 secuencias corresponden al grupo externo (morados). Los colores en las ramas del grupo interno corresponden a los resultados de Geneland, el grupo NORTE se presenta en verde, el grupo SUR se presenta en rojo, el grupo CENTRO en amarillo y al grupo de LOS TUXTLAS en azul. Con fines comparativos, a la derecha se muestra una barra en colores que indica los 7 grupos mencionados por ^{Guevara y Sánchez Cordero (2018a)}. 

Discusión

La distribución actual en el espacio geográfico y la divergencia genética entre los 4 grupos dentro del clado integrado por Cryptotis mexicanus, C. obscurus y C. nelsoni sugieren una serie de eventos de expansión, contracción y fragmentación geográfica como consecuencia, probablemente, de las fluctuaciones climáticas a lo largo del Pleistoceno tardío (^{Ramírez-Barahona y Eguiarte, 2013}). Las oscilaciones climáticas y la orografía intrincada del norte del Neotrópico habrían favorecido periodos alternados de conectividad y aislamiento poblacional, dando como resultado el patrón de variabilidad genética observado en estas musarañas (^{Barber y Klicka, 2010}). Si bien ésto no es nuevo y está de acuerdo con lo encontrado en diversos estudios previos (e.g., ^{Guevara, 2020}; ^{Ornelas et al., 2010}, ²⁰¹³), el presente estudio destaca el impacto diferencial de las barreras geográficas a lo largo del área de estudio.

El grupo NORTE se observa como una agrupación cohesiva, pero fragmentada geográficamente y con estructura genética interna. La población en el extremo boreal, aislada en la Reserva de la Biosfera El Cielo, Tamaulipas, muestra una estrecha relación con las muestras al sur. Algunas zonas bajas como las subcuencas de los ríos Tampaón, Moctezuma y Amajac en los estados de Hidalgo y Querétaro, afluentes del Pánuco, parecen restringir de manera intermitente el flujo genético entre ambas regiones, evidenciando el papel de barreras ribereñas en musarañas (^{Jacquet et al., 2013}; tabla 1, figs. 2, 3). Al noreste, las tierras bajas húmedas de la provincia Veracruzana, a lo largo de la costa del golfo de México, han mantenido en aislamiento a la Sierra de Otontepec, un grupo de montañas aisladas de hasta 1,300 m de altitud (^{Mayen-Zaragoza et al., 2019}). A pesar de este notable aislamiento actual, las muestras de Otontepec también están estrechamente relacionadas con el centro y norte de la Sierra Madre Oriental.

En general, la agrupación NORTE es biogeográficamente interesante, ya que apoya la predicción del desplazamiento del bosque de niebla hacia la vertiente oriental de la Sierra Madre Oriental durante el Último Máximo Glacial y se acentúa con la deglaciación en el evento 1 de Heinrich (hace entre 17,000 y 15,000 años; ^{Guevara, 2019}). Parece prudente sugerir que las poblaciones dentro del grupo NORTE estuvieron en contacto durante las fases glaciales, siguiendo la expansión del bosque de niebla, pero fragmentándose durante las fases interglaciales.

Con respecto al límite meridional del grupo NORTE, la evidencia genética indica que la diferenciación entre los grupos NORTE y CENTRO es alta y significativa (tabla 1, fig. 4). Con base en el análisis de discontinuidades genéticas, es posible inferir que las zonas más elevadas y frías de la Faja Volcánica Transmexicana estén actuando como una barrera para estas musarañas (e.g., la cara norte del volcán Pico de Orizaba), y no así el cañón del río Blanco como había sido sugerido por ^{Guevara y Sánchez- Cordero (2018a)}. De ser así, el límite entre ambos grupos apoya lo sugerido originalmente por ^{Choate (1970)}, quien mencionó que alrededor de los estados de Hidalgo, Puebla y Veracruz podría estar una zona de contacto entre C. obscurus y C. mexicanus.

El cañón del río Santo Domingo (ca. 450 m snm), por su parte, estaría actuando como una barrera intermitente o reciente entre los grupos CENTRO y SUR. La divergencia genética relativamente baja entre ambos grupos, la monofilia débilmente apoyada del grupo SUR y la alta similitud morfológica (^{Guevara y Sánchez-Cordero, 2018a}), sugieren un efecto tenue del cañón del río Santo Domingo en la divergencia de ambos grupos. El río Santo Domingo parece haber tenido un impacto en la diversificación de algunos pequeños mamíferos, tal como roedores de los géneros Habromys (^{Carleton et al., 2002}) y Megadontomys (^{Vallejo y González-Cózatl, 2012}), pero no tan fuerte para el género Handleyomys (^{Almendra et al., 2014}). Los efectos diferenciales entre especies podrían deberse a los atributos de cada taxón, tal como el tamaño corporal, hábitos y vagilidad.

Al sureste, el istmo de Tehuantepec es una franja de tierras bajas de 250 km de ancho en su punto más angosto, con algunos lomeríos de baja altitud (ca. 250 m snm) (^{Barrier et al., 1998}; ^{Ornelas et al., 2013}). Históricamente, se ha sugerido que el istmo es una barrera geográfica importante para la divergencia de diversos grupos taxonómicos, incluyendo las musarañas (^{Choate, 1970}; ^{Flores-Villela y Martínez-Salazar, 2009}; ^{Peterson et al., 1999}). Por lo tanto, el agrupamiento de los haplotipos desde el sur del río Santo Domingo hasta las Tierras Altas de Chiapas fue inesperado (i.e., grupo SUR). Sin embargo, las relaciones halladas en la red de haplotipos y en la inferencia bayesiana mostraron a los haplotipos, al este del istmo, como un grupo notablemente más distanciado del resto del grupo SUR, lo que evidencia que esos haplotipos son únicos.

Este efecto de barrera intermitente del istmo de Tehuantepec ha sido observado en otros vertebrados (^{Barber y Klicka, 2010}; ^{Hardy et al., 2013}; ^{Sullivan et al., 1997}, ²⁰⁰⁰) y en plantas (^{Ornelas et al., 2010}). Los pasos mutacionales hallados en los análisis sugieren que los haplotipos al este del istmo de Tehuantepec divergieron del grupo SUR en un tiempo relativamente similar a la divergencia entre los grupos CENTRO y SUR. Hasta la fecha, se sabe poco sobre la posible conectividad de los bosques de niebla a través del istmo de Tehuantepec en los ciclos glaciales-interglaciales ocurridos a través de los últimos 400,000 años (^{Campbell, 1984}). El registro de un espécimen de musaraña identificado como C. mexicanus (^{Choate, 1970}) en una zona montañosa en medio del istmo de Tehuantepec (California Academy of Sciences, Mammalogy Collection, 14636), sugiere una distribución continua de las musarañas a través de un puente de bosque en el pasado. Un mayor tamaño de muestra de C. mexicanus en esta región del istmo, así como de las Tierras Altas de Chiapas, podría ayudar a esclarecer la posible conectividad de especies del bosque de niebla a través del Cuaternario (^{Mendelson III y Campbell, 1999}; ^{Hardy et al., 2013}).

Finalmente, las tierras bajas de la provincia Veracruzana separan a la Sierra de Los Tuxtlas de la Sierra Madre Oriental y de la Sierra Madre del Sur. Con base en la divergencia genética observada en C. nelsoni (> 5 %, LOS TUXTLAS), esta especie se ha mantenido en aislamiento geográfico, probablemente durante al menos el Cuaternario tardío (^{Guevara, 2019}; ^{Woodman y Timm, 1999}). Algo similar ha ocurrido con otros taxones endémicos de Los Tuxtlas (e.g., anfibios, reptiles y aves; ^{Flores-Villela y Martínez-Salazar, 2009}; ^{Ornelas et al., 2013}). Lo anterior destaca el papel de este sistema montañoso como un promotor de especiación por aislamiento geográfico (^{Chávez-Pesqueira y Núñez-Farfán, 2016}; ^{Ruiz-Sánchez y Ornelas, 2014}), debido a la extensa llanura de tierras bajas a lo largo de la provincia Veracruzana, que potencialmente ha dificultado la dispersión de especies asociadas a ecosistemas de montaña.

Hasta la fecha, las evaluaciones taxonómicas de las especies del clado Cryptotis mexicanus, C. obscurus y C. nelsoni se han basado en caracteres morfológicos (^{Woodman y Timm, 1999}), lo cual ha brindado un avance notable en su clasificación taxonómica. Sin embargo, debido a la relativa baja variación morfológica entre las especies, aún quedan dudas acerca del límite entre especies y la posible existencia de linajes crípticos (^{Guevara y Sánchez-Cordero, 2018a}). Los patrones filogeográficos del gen Citb revelados en el presente estudio ofrecen una visión complementaria para esclarecer la distribución espacial de los linajes dentro de este grupo y asistir en la toma de decisiones taxonómicas.

Primero, el grupo de LOS TUXTLAS es un linaje independiente de los otros grupos cercanamente relacionados. La monofilia y los altos valores de apoyo, los pasos mutacionales y distancia genética con respecto al grupo más cercano, y el aislamiento geográfico, aportan evidencia concluyente para el reconocimiento específico de C. nelsoni (^{Guevara y Cervantes, 2014}; ^{He et al., 2015}; ^{Woodman y Timm, 1999}), un taxón que previamente había sido considerado subespecie de C. mexicanus (^{Choate, 1970}). Dentro del par de especies C. mexicanus y C. obscurus, el panorama es más complicado. Considerando los resultados de estudios previos y la evidencia que el gen Citb aporta en este estudio, consideramos necesaria una revisión taxonómica, en la que se podría evaluar la existencia de al menos 2 linajes evolutivos distribuidos desde el norte de la Sierra Madre Oriental hasta las Tierras Altas de Chiapas (^{Guevara y Sánchez-Cordero, 2018a}; ^{Ramírez-Pulido et al., 2004}).

El primer linaje comprende al grupo NORTE, distribuyéndose desde la Reserva de la Biosfera El Cielo en la Sierra Madre Oriental, la Sierra de Otontepec al norte de la provincia Veracruzana y hasta el extremo oriental de la Faja Volcánica Transmexicana. Dado que las localidades tipo de C. mexicanus y C. obscurus están dentro del grupo NORTE y atendiendo al principio de prioridad, C. mexicanus sería el nombre válido para este linaje, y como consecuencia, C. obscurus sería un sinónimo junior de C. mexicanus. El segundo linaje comprendería a los grupos CENTRO y SUR, desde el sur del Pico de Orizaba hasta las Tierras Altas de Chiapas. Sin embargo, no se debe descartar la posibilidad de que los grupos CENTRO y SUR también representen unidades independientes, aunque el nivel de diferenciación no es de la misma magnitud que el observado entre los grupos LOS TUXTLAS, NORTE y SUR. Dentro de los grupos CENTRO y SUR, el río Santo Domingo marcaría el límite entre sus distribuciones. En cuanto al caso de la población de las Tierras Altas de Chiapas, este estudio genético no apoya el reconocimiento de su estado específico, tal como había sido sugerido anteriormente con base en evidencia morfológica (^{Guevara y Sánchez-Cordero, 2018a}). En cualquier caso, es importante mencionar que no existen nombres disponibles para los linajes que abarcarían los grupos CENTRO y SUR. Un paso siguiente para resolver las incertidumbres taxonómicas deberá ser la inclusión de evidencia genética adicional, tal como genes nucleares, así como un mayor número de muestras provenientes del centro y sur de la Sierra Madre Oriental y al este del istmo de Tehuantepec. Además, es recomendable la evaluación de caracteres morfológicos del esqueleto poscraneal, los cuáles han sido útiles para la delimitación de especies del género.