evolución-ago-2025

Aldo Gómez-Benitez ^{a, b,} *, Erika Adriana Reyes-Velázquez ^{b, c}, Justin Rheubert ^d, Oswaldo Hernández-Gallegos ^c

^a Universidad Autónoma Metropolitana, Unidad Lerma, División de Ciencias Biológicas y de la Salud, Departamento de Ciencias Ambientales, Av. de las Garzas No. 10, El Panteón, 52005 Lerma, Estado de México, Mexico

^b Red de Investigación y Divulgación de Anfibios y Reptiles MX, Guadalupe Victoria No. 33, 56800 Ozumba de Alzate, Estado de México, México

^c Universidad Autónoma del Estado de México, Facultad de Ciencias, Laboratorio de Herpetología, Instituto Literario No. 100, Colonia Centro, 50000 Toluca, Estado de México, Mexico

^d University of Findlay, Department of Biology, 1000 N. Main St., Findlay, Ohio, 45840 USA 

*Corresponding author: gobeal940814@gmail.com (A. Gómez-Benitez)

Abstract

The evolutionary role of ontogenetic allometry is a topic of debate, as it can act as an evolutionary constraint or as a mechanism for functional optimization. We assessed the allometric changes in head shape of Crotalus triseriatus (Wagler, 1830) to understand how ontogenetic shifts in head morphology adapts to ecological demands as it grows. The sampling was conducted in Toluca de Lerdo, Estado de México, where dorsal-view photographs of the head were taken for the 22 collected specimens. We digitized 22 landmarks that represent the C. triseriatus head shape. We calculated the common allometric component and the residual shape components (RSC). Results showed statistically significant positive allometry, considering shape as a function of size, in the head shape of C. triseriatus. The main changes include a contraction along the sagittal plane, primarily in the snout region, and a marked expansion in the posterior part of the head. Size being the main source of variation appears to be a common pattern in snakes, though the functional source of the variation may relate to the natural history aspects of each species. In C. triseriatus, these changes seem to be associated with predatory efficiency.

Keywords: Crotalus triseriatus; Common allometric component; Geometric morphometrics; Natural history; Predatory efficiency; RSC

Alometría ontogenética y variación no dependiente del tamaño en la forma de la cabeza de la Serpiente de Cascabel Transvolcánica

Resumen 

El papel evolutivo de la alometría ontogenética es un tema debatido, pues puede actuar como una restricción evolutiva o como un mecanismo de optimización funcional. Evaluamos los cambios alométricos en la forma de la cabeza de Crotalus triseriatus (Wagler, 1830) para entender cómo los cambios ontogenéticos se ajustan a las demandas ecológicas a lo largo del crecimiento. El muestreo se llevó a cabo en Toluca de Lerdo, Estado de México, se tomaron fotografías dorsales de la cabeza de los 22 especímenes recolectados. Digitalizamos 22 puntos de referencia que representan la forma de la cabeza de C. triseriatus. Calculamos el componente alométrico común y los componentes residuales de la forma (CRF). Los resultados mostraron una alometría positiva estadísticamente significativa, considerando la forma como función del tamaño, en la cabeza de C. triseriatus. Los principales cambios incluyen una contracción a lo largo del plano sagital, en la región del hocico y una expansión marcada en la parte posterior de la cabeza. La talla como principal fuente de variación es un patrón común en serpientes, aunque su origen funcional puede depender de la historia natural de cada especie. En C. triseriatus, estos cambios están relacionados con la eficiencia depredadora.

Palabras clave: Crotalus triseriatus; Componente alométrico común; Morfometría geométrica; Historia natural; Eficiencia depredadora; CRF

Introduction

Allometry refers to the study of size and its biological consequences in morphology and function, typically focusing on how the proportions of body structures relate to each other at specific determined time or throughout development (Gould, 1966; Pretzsch, 2009). Modern studies of allometry are based on Huxley’s (1932) research, where he introduced the allometric equation, the concept of relative growth —different body parts grow at different rates— and hypothesized its evolutionary significance. Allometry is classified into 3 types: 1) static allometry which describes the relationship between body parts within a population at a given developmental stage, 2) evolutionary allometry which examines these relationships across populations or species, and 3) ontogenetic allometry which investigates changes in body proportions during growth, providing key insights into developmental processes (Pélabon et al., 2014). The evolutionary significance of allometry is debated, as it can act both as an evolutionary constraint, by which differential growth relationships restrict the range of phenotypic variation available for selection, and as a mechanism for functional optimization (Le Verger et al., 2023; Voje et al., 2022). While the allometric relationships observed during organisms’ development can act as evolutionary constraints, guiding species along fixed developmental trajectories (Voje et al., 2022), these same relationships are also shaped by natural selection to optimize functional traits, highlighting an adaptive dimension of allometry (Le Verger et al., 2023; Strelin & Diggle, 2022). 

In reptiles, allometric patterns vary widely. In the fossorial species Cynisca leucura (Amphisbaenia), allometry acts as an evolutionary constraint, limiting cranial morphology to maintain proportions suitable for head-first digging (Hipsley et al., 2016). Sexual selection also influences allometry, as seen in Anolis species, where divergent dewlap growth trajectories are linked to reproductive pressures, though all exhibit positive allometry (Stroud et al., 2023). Similarly, in snakes of the genus Malacophagus, sexual dimorphism is reflected in distinct allometric trajectories between males and females (Dos Santos et al., 2022). In crocodilians, dietary demands drive allometry, enabling the consumption of larger prey as the animal grows (Balaguera-Reina et al., 2024; Dodson, 1975). These examples demonstrate that reptile allometry exhibits diverse patterns, making this group ideal for studying allometry. A robust method to assess the relationship between shape and size is the calculation of the Common Allometric Component (CAC), as described by Mitteroecker et al. (2004). The CAC captures allometric variation by summarizing how shape changes as a function of size across all specimens. This metric has been previously applied in reptile studies to evaluate evolutionary shifts in squamate cranial morphology during development and its relationship with whole-head integration (Hipsley & Müller, 2017; Ollonen et al., 2024), proving to be an excellent tool for studying allometry. 

We selected the Mexican Dusky Rattlesnake (Crotalus triseriatus) as a study model because, as a fully terrestrial predator with clearly defined foraging and dietary habits, it offers a well‑characterized system to link ontogenetic shape changes with ecological and biotic interactions. Crotalus triseriatus (Wagler, 1830), is a medium-sized venomous snake endemic to Mexico (Fig. 1A), found in the Trans-Mexican Volcanic Belt across the states of Veracruz, Puebla, Tlaxcala, Morelos, Estado de México, and Michoacán, at elevations between 1,905 and 4,572 m asl (Campbell & Lamar, 2004; Castillo-Juárez et al., 2021). Their primary habitat is coniferous forests and oak forest (primarily pine-oak forests), grasslands, montane cloud forests, and agroecosystems (Campbell & Lamar, 2004; Castillo-Juárez et al., 2021; Gómez-Benitez, 2023). This species is diurnal to crepuscular, preying on lizards, rodents, and salamanders; additionally, C. triseriatus exhibits a distinct ontogenetic dietary shift, transitioning from ectothermic prey to mammals as the snakes grow (Mociño-Deloya et al., 2014; Ramírez-Bautista et al., 2009).

In the present study, we analyzed the allometric changes in head shape of C. triseriatus using geometric morphometrics, a set of multivariate statistical methods that quantify organismal shape by analyzing the Cartesian coordinates of homologous anatomical landmarks, to characterize the ontogenetic variation in head morphology of a species known for exhibiting ontogenetic shifts in its ecological interactions. Understanding these morphological changes provides valuable insights into how this species adapts to ecological demands as it grows and contributes to a greater understanding of allometry in an evolutionary sense.

Materials and methods

We conducted fieldwork in an agroecosystem located in El Cerrillo, Piedras Blancas, within the municipality of Toluca de Lerdo, Estado de México, Mexico. The study area is characterized by an arboreal stratum primarily composed of Salix babilonica, Cupressus sp. and Casuarina equisetifolia, the herbaceous stratum contains over 118 species of angiosperms (Álvarez-Lopeztello et al., 2016). The study site is disturbed due to agricultural practices (many of them unregulated), livestock, constant human presence and intentional killing of snakes. We searched for rattlesnakes in various microhabitats, including under rocks, inside and alongside a cement irrigation canal, and under polypropylene plastic sheets previously used as waterproofing membranes in agricultural activities at the site. Field researchers wore protective snake gaiters and bait gloves, and rattlesnakes were handled using herpetological hooks. The research team has a detailed snakebite emergency protocol, which includes always carrying sufficient vials of antivenom to ensure support dose administration if required. We took dorsal photographs of the collected individuals heads using a Nikon D3500 camera (settings: 1/200 shutter speed, f/22 aperture, ISO 100-400) equipped with a Nikkor 70-300 mm lens at its maximum focal length. To avoid the limitation of the lens minimum focusing distance, we used macro extension tubes. The camera was mounted on a tripod and positioned 60 cm above the specimens. Additionally, the photographs were illuminated using the camera built-in flash to ensure consistent lighting. Only one dorsal photograph of the head was taken per individual for geometric morphometric analysis. Snakes were handled following standard safety protocols: individuals were manipulated using herpetological hooks and gently positioned over a millimetric sheet, ensuring that the head lay flat and straight. When necessary, restraining tubes were used to guide the head into position. All photographs were taken in the field; specimens were not transported to the laboratory. Instead, to guarantee photograph quality, individuals were briefly placed in cloth bags and moved to a flat surface present at the study site. All individuals were sexed through hemipenes eversion, measured in SVL and released where they were captured.

Using TPSDig 2.31, we digitized 22 anatomical landmarks on the dorsal view of the head of C. triseriatus (Fig. 1B, Table 1), in duplicate. Landmarks were placed at intersections of external cranial scales to avoid invasive procedures, as all specimens were alive. These points were selected based on their anatomical relevance and their ability to capture overall head shape, with particular emphasis on regions of functional and evolutionary interest such as the snout, the orbital area, and the venom delivery system. In addition, the chosen landmarks enable the analysis of head proportions, including width and length. Because our aim was to study ontogenetic allometry, individuals were analyzed jointly regardless of age class, in order to capture variation in shape throughout development. A Procrustes superimposition was then performed in MorphoJ (Klingenberg, 2011) to remove differences in size, position, and orientation according with the methodology proposed by Dryden and Mardia (1998). The data were projected onto the tangent space of the shape space, enabling linear analysis of shape variation (Dryden & Mardia, 1998). This method retains centroid size information for subsequent analyses, such as allometry studies.

Table 1

Position of the 22 landmarks selected to represent the head shape of Crotalus triseriatus (see Fig. 1B).

Landmarks	Description
1	Between the internasal scales, posterior to rostral
2, 3	Intersection between the first nasal and internasal scales
4, 7	Intersection between the second nasal, intercanthal and prefrontal scales
5, 6	Intersection between internasal, intercanthal and prefrontal scales
8, 9	Intersection between intercanthal, prefrontal and the first row of frontal scales
10, 13	Distal intersection between prefrontal and supraocular scales
11, 12	Intersection between prefrontal, supraocular and second row of frontal scales
14, 15	Anterior intersection of supraocular and intersupraocular scales
16, 17	Intersection between the supraocular scales and the first row of the scales in the crown of the head
18, 19	Posterior part of the most distal fourth row of the scales in the crown of the head
20, 21	Posterior part of the most distal eighth row of the scales in the crown of the head
22	Intersection between both intercanthal and first row of frontal scales

After Procrustes superimposition, we conducted a multiple linear regression in MorphoJ with the Cartesian coordinates of the 22 landmarks (44 dependent variables: X and Y for each landmark) as the response and the log‑transformed centroid size as the predictor. This approach allowed us to quantify how head shape (landmark coordinates) changes in relation to overall head size (log‑transformed centroid size). To assess the statistical power of our regression analysis, we conducted a power analysis using the pwr package (Champely, 2020) in R (R Core Team, 2021). Results were considered significative at an α = 0.05.

To assess the changes related to size in C. triseriatus, and the non-size-dependent variation, we employed the CAC and the residual shape components (RSC). The CAC was calculated conducting a linear regression of shape as a function of size (Mitteroecker et al., 2004). The RSC represents the variation in shape that is not attributable to size. The residuals obtained during the linear regression were analyzed using Principal Component Analysis (PCA) to extract RSC (Mitteroecker et al., 2004). To correctly display the changes detected with the CAC and RSC, we calculated the Jacobian expansion factor using PAST4 version 4.17 (Hammer et al., 2001). This calculation involves deriving a Jacobian matrix from a thin-plate spline deformation (or warp), which quantifies local shape changes by relating the derivatives of the deformed configuration to a consensus configuration, or from smaller to bigger individuals in the case of CAC. The expansion factor, determined as the determinant of the Jacobian matrix, reflects the local area distortion resulting from the deformation, providing insight into the degree of expansion or contraction across regions of the shape.

Finally, we plotted the CAC against log-transformed centroid size and against the RSC. We also calculated the standard deviation to understand the patterns of variation related to allometry and non-size-dependent factors. These exploratory analyses allow us to visualize the relationships between size and shape variations, and how non-allometric factors contribute to morphological changes in C. triseriatus.

Results

We collected a total of 22 individuals, including 13 juveniles and 9 adults, based on SVL, including 7 females and 15 males. The multiple linear regression results show a significant relationship between head shape and log centroid size (R² = 0.4992, p = 0.0122), indicating a positive allometry (Fig. 2A). Our statistical power analysis revealed that the regression model had a power of 0.993, indicating a high resolution of our model with a 99.3% probability of detecting a significant relationship between head shape and log centroid size in C. triseriatus. In other words, there is a very low probability (0.7%) of committing a Type II error (Whitlock & Schluter, 2008).

The CAC analysis (Fig. 2A) revealed the following changes in the head shape of C. triseriatus: a contraction between the internasal scales and the rostral, resulting in a more square-shaped snout in adults than in juveniles (landmarks 1-3); a contraction of the prefrontal scales over the intercanthals, reducing the proportion of space they occupy (landmarks 8, 9, and 22); a slight expansion on both sides of the intra supraocular scales (landmarks 8, 11, and 14 on the right side, and landmarks 9, 12, and 15 on the left side); another significant contraction in the scales at the crown of the head, right in the center at the level of the back of the supraocular scales (landmarks 16 and 17); and finally, the most notable expansion in the anterior sector of the head (landmarks 18-21) (Fig. 2A).

To interpret the variation not related to size in the head shape of C. triseriatus, we used the first 3 RSCs, as these explained 66.30% of the variation (RSC1 = 37.44%, RSC2 = 15.73%, RSC3 = 13.13%). The principal sources of variation in RSC1 are an expansion of the prefrontal scales to the external edges (landmarks 10-13), a significant contraction between the first and fourth row of scales on the crown of the head (landmarks 16-19), and the most important source of variation in this component was an expansion in the posterior sector of the head (landmarks 20 and 21) (Fig. 3), similar to the pattern observed in the CAC (Fig. 2B). In RSC2, 3 important expansions occur throughout the head of C. triseriatus: one in the intercanthal scales (landmarks 5, 6, 8, 9, and 22), a second in the supraocular scales (landmarks 10, 11, 14, and 16 on the right side and landmarks 12, 13, 15, and 17 on the left side of the head), and the last one in the middle of the head in the scales of the crown (Fig. 3). Lastly, in RSC3, an important expansion occurs at the center of the snout (landmarks 1-9), delimited on both sides by a contraction in the prefrontal scales (landmarks 4, 11, and 10 on the right side and landmarks 7, 12, and 13 on the left side); another important expansion occurs in the first 4 rows of the scales on the crown of the head (landmarks 16-19), followed posteriorly by a contraction (landmarks 20 and 21) (Fig. 3).

The standard deviations of the CAC (σ = 0.0288) and the first 3 RSCs (RSC1 σ = 0.0356, RSC2 σ = 0.0231, RSC3 σ = 0.0211) indicate that there are sources of variation in the head shape of C. triseriatus that are significant just as the size-related variation, particularly those contained in RSC1 (Fig. 3).

Discussion 

According to our results, Crotalus triseriatus exhibits positive allometry in head shape, with important ontogenetic changes occurring throughout the entire head. Our analysis revealed that the CAC explains approximately 49.92% of the variation in head shape of C. triseriatus, indicating that size is a major source of variation. This appears to be a common pattern in snakes (Abegg et al., 2020; Lucchini et al., 2020; Murta-Fonseca & Fernandes, 2016), though the specific ways in which the head shape varies can differ within species. In C. triseriatus, the primary expansions occur in the posterior region of the head, a pattern consistent with other viperids (Lucchini et al., 2020), while contractions are observed mainly along the sagittal plane. In contrast, species from other families display more prominent snout expansions and modifications in the mid anterior-posterior section of the head (Abegg et al., 2020; Murta-Fonseca & Fernandes, 2016). These disparities may be attributed to the natural history of different species, as allometric trajectories are shaped by ecological factors such as diet or microhabitat preferences (Kaliontzopoulou et al., 2012; Simonsen et al., 2017). In snakes for which allometry in head shape has been analyzed through geometric morphometrics, it has been shown that semi-fossorial and aquatic habits primarily drive changes in the mid-section of the head (Abegg et al., 2020; Murta-Fonseca & Fernandes, 2016). In contrast C. triseriatus, a completely terrestrial species (Campbell & Lamar, 2004), shows more pronounced changes in the posterior part of the head. This suggests that environmental factors may impact head shape modifications in species adapted to different habitats, concretely in C. triseriatus, recent work on multiple populations has confirmed a significant environment–head shape association within this species (Caballero‑Viñas et al., 2025).

In C. triseriatus, the allometric changes observed in the snout may also be related to ontogenetic dietary shifts. Neonates primarily consume ectothermic prey, such as amphibians and lizards, while adults predominantly prey on mammals (Mociño-Deloya et al., 2014). Ontogenetic cranial changes associated with diet are well-documented in snakes, as these modifications are necessary to develop a wider gape and enhanced jaw musculature for prey capture (Patterson et al., 2021; Vincent et al., 2007). Such changes highlight the importance of head shape adaptations throughout an organism life to optimize survival. The greater expansion in the posterior part of the head likely reflects the development of venom glands and the muscularis compressor glandulae, which in adults may occupy a proportionally larger area of the head. Further studies conducting a direct anatomical assessment of allometric changes in the venom delivery system though dissection, are required to confirm that relationship. This anatomical adaptation likely enhances predatory efficiency, allowing adults to capture and subdue larger prey (such as mammals) as they grow. The hypothesis that ontogenetic dietary shifts or broader predatory requirements are the primary drivers of allometric variation in C. triseriatus is further supported by evidence that in Crotalus, allometry in head shape is linked to the kinematics of predatory strikes, specifically, maximum strike velocity (positive allometry) and the proportion of the body involved in the strike, as well as timing variables (negative allometry), change as rattlesnakes grow (LaDuc, 2003).

The residual shape components in C. triseriatus head shape, which account for the remaining 50.07% variation that is not explained by size, highlights significant shape differences unrelated to size. This suggests that while size is a key factor, other important sources of variation influence head morphology. In amphibians and reptiles, head shape variation has been linked to factors such as phylogenetic history, habitat use (Openshaw & Keogh, 2014), sexual dimorphism (Alarcón-Ríos et al., 2017), and dietary or ecological interactions (Alarcón-Ríos et al., 2017; Segall et al., 2020; Urošević et al., 2014). This raises a fundamental question: if allometry explains a significant portion (half) of head shape variation in C. triseriatus and this is driven by predatory requirements, is it not ultimately ecology that explains this variation more than size by itself? In other words, do ecological interactions adapt to changes in size, or does size adapt to ecological niches available? The answer lies in a bidirectional relationship. Ecological interactions can shape allometric trajectories (Maestri et al., 2015); that is, certain traits may adapt to available resources, predation pressures, or competition as an organism grows, thus growth itself is limited by ecological constraints. However, allometry often acts as an evolutionary constraint, limiting morphological variation by imposing functional requirements (Huxley, 1932). In C. triseriatus, the strong correlation between size and head shape likely reflects functional limitations, particularly regarding prey capture and venom delivery. These evolutionary constraints imply that although head shape evolves in response to ecological pressures, these changes are not entirely free to vary but are shaped by underlying allometric relationships that ensure functionality. Further research into the “what came first?” question is necessary to better understand the role of allometry as an evolutionary constraint and the extent to which ecology constrains allometry.

The allometric changes in the head shape of C. triseriatus demonstrated to be a major source of variation in this species, highlighting the importance of allometry as an evolutionary constraint. In a predator such as rattlesnakes, predatory efficiency is closely linked to morphological adaptations, particularly in the head and jaw regions. Our findings suggest that ontogenetic changes in head shape are not merely a consequence of growth but instead reflect functional adjustments necessary for successful predation. As individuals grow, their ability to capture and process larger prey increases, which likely enhanced survival, just as has been stated in crocodilians (Barrios-Quiroz et al., 2012; Falcón-Espitia & Jerez, 2021). Moreover, the observed posterior expansion of the head in C. triseriatus, coupled with the contractions along the sagittal plane, could underscore the specialized adaptations linked to venom delivery and jaw musculature development. These morphological adjustments not only facilitate the consumption of larger mammalian prey but may also influence predator-prey dynamics and human-snake interactions by modulating strike performance and venom efficiency. Understanding these allometric trajectories is therefore crucial for both ecologically informed management of rattlesnake populations and for mitigating human-wildlife conflict.

Acknowledgements

This study was conducted under the Semarnat permit SGPA/DGVS/10352/21. We are also grateful to Conahcyt for the doctoral scholarship awarded to AGB. Finally, we thank the students of the Herpetology Laboratory for their support during fieldwork.

References

Abegg, A., Passos, P., Mario-da-Rosa, C., Azevedo, W., Malta-Borges, L., & Bubadué, J. (2020). Sexual dimorphism, ontogeny and static allometry of a semi-fossorial snake (genus Atractus). Zoologischer Anzeiger, 287, 95–104. https://doi.org/10.1016/j.jcz.2020.05.008 

Alarcón-Ríos, L., Velo-Antón, G., & Kaliontzopoulou, A. (2017). A non-invasive geometric morphometrics method for exploring variation in dorsal head shape in urodeles: Sexual dimorphism and geographic variation in Salamandra salamandra. Journal of Morphology, 278, 475–485. https://doi.org/10.1002/jmor.20643 

Álvarez-Lopeztello, J., Rivas-Manzano, I. V., Aguilera-Gómez, L. I., & González-Ledesma, M. (2016). Diversity and structure of a grassland at El Cerrillo, Piedras Blancas, Estado de México, Mexico. Revista Mexicana de Biodiversidad, 87, 980–989. https://doi.org/10.1016/j.rmb.2016.06.006 

Balaguera-Reina, S. A., Mason, B. M., Brandt, L. A., Hernández, N. D., Daykin, B. L., McCaffrey, K. L. et al. (2024). Ecological implications of allometric relationships in American alligators (Alligator mississippiensis). Scientific Reports, 14, 6140. https://doi.org/10.1038/s41598-024-56798-5 

Barrios-Quiroz, G., Casas-Andreu, G., & Escobedo-Galván, A. H. (2012). Sexual Size dimorphism and allometric growth of Morelet’s crocodiles in captivity. Zoological Science, 29, 198–203. https://doi.org/10.2108/zsj.29.198

Caballero-Viñas, C., Arenas, S., Jaramillo-Alba, J. L., Pérez-Mendoza, H. A., Manjarrez, J., Domínguez-Vega, H. et al. (2025). Mexican dusky rattlesnakes (Crotalus triseriatus) on Mexican sky islands: morphometric variation and operative temperature relationships with local environmental variables. Biologia, 80, 1389–1406. https://doi.org/10.1007/s11756-025-01905-8 

Campbell, J. A., & Lamar, W. W. (2004). The venomous reptiles of the western hemisphere. Ithaca, New York: Comstock Publishing Associates.

Castillo-Juárez, J. L., Vásquez-Cruz, V., Taval-Velázquez, L. P., Avalos-Vela, R., & Lara-Hernández, F. A. (2021). Notas de distribución e historia natural de Crotalus triseriatus (Viperidae) en la región de las altas montañas, Veracruz y en el estado de Puebla. Revista Latinoamericana de Herpetología, 4, 133–139. https://doi.org/10.22201/fc.25942158e.2021.02.213 

Champely, S. (2020). pwr: Basic Functions for Power Analysis. Université de Lyon, Lyon. Retrieved 13 may, 2025 from: https://cran.r-project.org/web/packages/pwr/index.html 

Dodson, P. (1975). Functional and ecological significance of relative growth in Alligator. Journal of Zoology, 175, 315–355. https://doi.org/10.1111/j.1469-7998.1975.tb01405.x 

Dos Santos, M. M., Klaczko J., & da Costa-Prudente, A. L. (2022). Sexual dimorphism and allometry in Malacophagus snakes (Dipsadidae: Dipsadinae). Zoology (Jena), 153, 126026. https://doi.org/10.1016/j.zool.2022.126026 

Dryden, I. L., & Mardia, K. V. (1998). Statistical shape analysis. Chichester, England: Wiley.

Falcón-Espitia, N., & Jerez, A. (2021). The skull of Caiman crocodilus fuscus: Allometric andontogenetic shifts. Acta Biológica Colombiano, 27, 458–463. https://doi.org/10.15446/abc.v27n3.90810 

Gómez-Benitez, A. (2023). Inestabilidad en el desarrollo y canalización de una comunidad de reptiles en un hábitat perturbado (Ph.D. Thesis). Facultad de Ciencias, Universidad Autónoma del Estado de México. Estado de México, Mexico. 

Gould, S. J. (1966). Allometry and size in ontogeny and phylogeny. Biological Review, 41, 587–640. https://doi.org/10.1111/j.1469-185X.1966.tb01624.x 

Hammer, Ø., Harper, D. A. T., & Paul, D. R. (2001). Past: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontologia Electronica, 4, 1–9.

Hipsley, C. A., & Müller, J. (2017). Developmental dynamics of ecomorphological convergence in a transcontinental lizard radiation. Evolution, 71, 936–948. https://doi.org/10.1111/evo.13186 

Hipsley, C. A., Rentinck, M. N., Rödel, M. O., & Müller, J. (2016). Ontogenetic allometry constrains cranial shape of the head-first burrowing worm lizard Cynisca leucura (Squamata: Amphisbaenidae). Journal of Morphology, 277, 1159–1167. https://doi.org/10.1002/jmor.20564

Huxley, J. S. S. (1932). Problems of relative growth. Baltimore, Maryland: Johns Hopkins University Press.

Kaliontzopoulou, A., Adams, D. C., van der Meijden, A., Perera, A., & Carretero, M. A. (2012). Relationships between head morphology, bite performance and ecology in two species of Podarcis wall lizards. Evolutionary Ecology, 26, 825–845. https://doi.org/10.1007/s10682-011-9538-y 

Klingenberg, C. P. (2011). MorphoJ: an integrated software package for geometric morphometrics. Molecular Ecology Resources, 11, 353–357. https://doi.org/10.1111/j.1755-0998.2010.02924.x 

LaDuc, T. J. (2003). Allometry and size evolution in the rattlesnake, with emphasis on predatory strike performance (Ph.D. Thesis). The University of Texas at Austin, USA. 

Le Verger, K., Hautier, L., Gerber, S., Bardin, J. P., Delsuc, F., González Ruiz, L. R. et al. (2023). Pervasive cranial allometry at different anatomical scales and variational levels in extant armadillos. Evolution, 78, 423–441. https://doi.org/10.1093/evolut/qpad214 

Lucchini, N., Kaliontzopoulou, A., Val, G., & Martínez-Freiría, F. (2020). Sources of intraspecific morphological variation in Vipera seoanei: Allometry, sex, and colour phenotype. Amphibia-Reptilia, 42, 1–16. https://doi.org/10.
1163/15685381-bja10024 

Maestri, R., Fornel, R., Freitas, T. R. O., & Marinho, J. R. (2015). Ontogenetic allometry in the foot size of Oligoryzomys flavescens (Waterhouse, 1837) (Rodentia, Sigmodontinae). Brazilian Journal of Biology, 75, 435–441. https://doi.org/10.1590/1519-6984.16613

Mitteroecker, P., Gunz, P., Bernhard, M., Schaefer, K., & Bookstein, F. L. (2004). Comparison of cranial ontogenetic trajectories among great apes and humans. Journal of Human Evolution, 46, 679–698. https://doi.org/10.1016/j.jhevol.2004.03.006 

Mociño-Deloya, E., Setser, K., & Pérez-Ramos, E. (2014). Observations on the diet of Crotalus triseriatus (Mexican Dusky Rattlesnake). Revista Mexicana de Biodiversidad, 85, 1289–1291. https://doi.org/10.7550/rmb.43908 

Murta-Fonseca, R., & Fernandes, D. (2016). The skull of Hydrodynastes gigas (Duméril, Bibron & Duméril, 1854) (Serpentes: Dipsadidae) as a model of snake ontogenetic allometry inferred by geometric morphometrics. Zoomorphology, 135, 233–241. https://doi.org/10.1007/s00435-015-0297-0 

Ollonen, J., Khannoon, E. R., Macrì, S., Vergilov, V., Kuurne, J., Saarikivi, J. et al. (2024). Dynamic evolutionary interplay between ontogenetic skull patterning and whole-head integration. Nature Ecology & Evolution, 8, 536–551. https://doi.org/10.1038/s41559-023-02295-3 

Openshaw, G., & Keogh, J. (2014). Head shape evolution in monitor lizards (Varanus): Interactions between extreme size disparity, phylogeny and ecology. Journal of Evolutionary Biology, 27, 363–373. https://doi.org/10.1111/jeb.12299 

Patterson, M. B., Wolfe, A. K., Fleming, P., Bateman, P., Martin, M. L., Sherratt, E. et al. (2021). Ontogenetic shift in diet of a large elapid snake is facilitated by allometric change in skull morphology. Evolutionary Ecology, 36, 489–509. https://doi.org/10.1007/s10682-022-10164-x 

Pélabon, C., Firmat, C., Bolstad, G., Voje, K., Houle, D., Cassara, J. A. et al. (2014). Evolution of morphological allometry. Annals of the New York Academy of Sciences, 1320, 58–75. https://doi.org/10.1007/10.1111/nyas.12470 

Pretzsch, H. (2009). Re-evaluation of allometry: State-of-the-art and perspective regarding individuals and stands of woody plants. In U. Lüttge, W. Beyschlag, B. Büdel, & D. Francis (Eds.), Progress in botany (pp. 339–369). Heidelberg, Germany: Springer.

R Core Team (2021). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. recovered 13 may, 2025 from: https://www.R-project.org/ 

Ramírez-Bautista, A., Hernández-Salinas, U., García-Vázquez, U. O., Leyte-Manrique, A., & Canseco-Márquez, L. (2009). Herpetofauna del valle de México: diversidad y conservación. Pachuca, Hidalgo: Universidad Autónoma del Estado de Hidalgo/ Comisión Nacional para el Conocimiento y Uso de la Biodiversidad.

Segall, M., Cornette, R., Godoy-Diana, R., & Herrel, A. (2020). Exploring the functional meaning of head shape disparity in aquatic snakes. Ecology and Evolution, 10, 6993–7005. https://doi.org/10.1002/ece3.6380 

Simonsen, M. K., Siwertsson, A., Adams, C. E., Amundsen, P. A., Præbel, K., & Knudsen, R. (2017). Allometric trajectories of body and head morphology in three sympatric Arctic charr (Salvelinus alpinus (L.)) morphs. Ecology and Evolution, 7, 7277–7289. https://doi.org/10.1002/ece3.3224 

Strelin, M. M., & Diggle, P. K. (2022). Within-individual leaf allometry and the evolution of leaf morphology:
A multilevel analysis of leaf allometry in temperate Viburnum (Adoxaceae) species. Evolution & Development, 24, 145–157. https://doi.org/10.1111/ede.12414 

Stroud, J. T., Petherick, A., Krasnoff, B., Walker, K., Suh, J. J., & Losos, J. B. (2023). Signal size allometry in Anolis lizard dewlaps. Biology Letters, 19, 20230160. https://doi.org/10.1098/rsbl.2023.0160 

Urošević, A., Ljubisavljević, K., & Ivanović, A. (2014). Variation in skull size and shape of the Common wall lizard (Podarcis muralis): Allometric and non-allometric shape changes. Contributions to Zoology, 83, 67–77. https://doi.org/10.1163/18759866-08301003 

Vincent, S., Moon, B. R., Herrel, A., & Kley, N. (2007). Are ontogenetic shifts in diet linked to shifts in feeding mechanics? Scaling of the feeding apparatus in the banded watersnake Nerodia fasciata. Journal of Experimental Biology, 210, 2057–2069. https://doi.org/10.1242/jeb.02779 

Voje, K. L., Bell, M. A., & Stuart, Y. E. (2022). Evolution of static allometry and constraint on evolutionary allometry in a fossil stickleback. Journal of Evolutionary Biology, 35, 423–438. https://doi.org/10.1111/jeb.13984 

Whitlock, M. C., & Schluter, D. (2008). The analysis of biological data. New York: Roberts and Company Publishers.

Marina Simota-Ruiz ^a, Oscar Mikery ^b,Juan Cisneros ^a, Oscar Carmona-Castro ^a, Alfredo Castillo-Vera ^a, *, Alma E. García-Morales ^c

^a El Colegio de la Frontera Sur, Ecología de Artrópodos y Manejo de Plagas, Carretera Antiguo Aeropuerto Km 2.5, 30700 Tapachula, Chiapas, Mexico

^b Ortomolecular Mexicana BC, Grupo Molecular, Investigación, Calle Fray Justo Sierra No. 9216, Colonia Marrón, 22015 Tijuana, Baja California, Mexico

^c El Colegio de la Frontera Sur, Departamento de Sistemática y Ecología Acuática, Avenida Centenario Km. 5.5, 77014 Chetumal, Quintana Roo, Mexico

*Corresponding author: acastill@ecosur.mx (A. Castillo-Vera)

Received: 28 February 2024; accepted: 14 February 2025

Abstract

The nematode Metaparasitylenchus hypothenemi (Poinar) (Tylenchida: Allantonematidae) is a parasite of the coffee berry borer (Hypothenemus hampei), one of the most significant pests affecting coffee crops. This study aimed to analyze the genetic variation of M. hypothenemi using the mitochondrial cytochrome oxidase subunit I (COI) gene. Reproductive females of M. hypothenemi were extracted from parasitized wild coffee berry borers for DNA extraction, amplification, and sequencing. Phylogenetic analyses revealed 2 well-differentiated lineages and 6 haplotypes distributed across 18 populations along a 100 km transect. Global genetic diversity was moderate (Hd = 0.52 ± 0.06), with a dominant haplotype present in 16 populations and others being representative or geographically isolated. Populations showed high genetic differentiation and restricted gene flow. These findings confirm the utility of the COI gene for analyzing the genetic variation of M. hypothenemi and provide a molecular basis for future studies on its biology and management.

Keywords: Obligate endoparasite; Insect pest; Genetic variation; COI Gene; Phylogenetic tree; Genetic structure

Variabilidad genética del nemátodo Metaparasitylenchus hypothenemi (Tylenchida:Allantonematidae), parásito natural de la broca del café (Coleoptera: Curculionidae: Scolytinae) usando secuencias de COI

Resumen

El nemátodo Metaparasitylenchus hypothenemi (Poinar) (Tylenchida: Allantonematidae) es un parásito de la broca del café (Hypothenemus hampei), una de las principales plagas del cultivo de café. Este estudio tuvo como objetivo analizar la variación genética de M. hypothenemi mediante el gen mitocondrial citocromo oxidasa subunidad I (COI). Se obtuvieron hembras reproductoras de M. hypothenemi de brocas silvestres parasitadas para la extracción, amplificación y secuenciación de DNA. Los análisis filogenéticos revelaron 2 linajes bien diferenciados y 6 haplotipos distribuidos en 18 poblaciones a lo largo de un transecto de 100 km. La diversidad genética global fue moderada (Hd = 0.52 ± 0.06), con un haplotipo dominante en 16 poblaciones y otros representativos o aislados geográficamente. Las poblaciones mostraron una alta diferenciación genética y flujo genético restringido. Estos resultados confirman la utilidad del gen COI para analizar la variación genética de M. hypothenemi y proporcionan bases moleculares para futuros estudios sobre su biología y manejo.

Palabras clave: Endoparásito obligado; Plaga de insectos; Variación genética; Gen COI; Árbol filogenético; Estructura genética

Introduction

The nematode, Metaparasitylenchus hypothenemi (Poinar) (Tylenchida: Allantonematidae), is an obligate endoparasite of the coffee berry borer (CBB), Hypo-
thenemus hampei (Ferrari) (Coleoptera: Curculionidae, Scolytinae), the most important insect pest of coffee worldwide (Le Pelley, 1968). Metaparasitylenchus hypothenemi was discovered in a commercial coffee plantation in southeastern Mexico attacking H. hampei adults (Castillo et al., 2002). M. hypothenemi reduce the fecundity and longevity of the females H. hampei (Castillo et al., 2019). This parasite-host relationship highlights the importance of studying its genetic variability and understanding its ecological role. The study of genetic variability in parasitic nematodes like M. hypothenemi requires molecular tools that enable precise population characterization. Among these, the mitochondrial gene cytochrome oxidase subunit I (COI) has proven to be an efficient marker for identifying specimens, clarifying evolutionary relationships among species, and analyzing population divergence (Hebert et al., 2003; Marsjan & Oldenbroek, 2007). This marker has been successfully used to analyze population structure and infer phylogenetic relationships in related nematodes, including parasites and entomopathogens such as Heterorhabditis marelatus and Deladenus proximus (Blouin et al., 1999; Fitza et al., 2019; Hartshorn et al., 2017; Saeb & David, 2014). However, there are few nematode sequences in both the BOLD Systems database and GenBank, and those corresponding to insect-parasitic nematodes are even scarcer, with M. hypothenemi being a species without previous genetic studies. The objective of this study was to analyze the genetic variability, in terms of the phylogeny, frequency of haplotypes and genetic differentiation between popula-
tions of the nematode Metaparasitylenchus hypothenemi collected from Soconusco, the only region in the world where this nematode has been systematically recorded.

Materials and methods

Five reproductive females of M. hypothenemi were collected in each of 18 localities from the Soconusco region, Chiapas, Mexico (Table 1, Fig. 1). Reproductive females of M. hypothenemi were extracted from the abdominal cavity of H. hampei, which were obtained from 100 coffee berries collected from coffee plants per locality. Reproductive females were selected for sampling due to their larger size compared to males, which facilitated DNA extraction and ensured the viability of the samples. Moreover, their genetic material plays a key role in transmitting genetic diversity to subsequent generations (Stewart & Larsson, 2014). The nematodes were stored in 2.5 ml microtubes containing 96% alcohol at -20 °C, until molecular analysis. The samples obtained in each sampled locality were considered as a population.

DNA extraction, amplification, and sequencing. A total of 90 nematodes of M. hypothenemi were used for DNA extraction. Total DNA from each specimen was extracted from the body complete of the nematode using a standard glass fiber method (Ivanova et al., 2006). A barcode region of approximately 658 base pairs (bp) of the mitochondrial COI gene was amplified using the primers ZplankF1_t1 (5’-TGTAAAACGACGGCCAGTTCTASWAATCATAARGATATTGG-3’) and ZplankR1_t1 (5’-CAGGAAACAGCTATGACTTCAGGRTGRCCRAARAATCA-3’) (Prosser et al., 2013). The final volume of the PCR mix was 12.5 µl, containing 0.12 µl of each primer (0.01 µM), 0.6 µl of MgCl₂ (50 mM), 2 µl of ultrapure water, 0.06 µl of each dNTP (0.05 mM), 1.25 µl of 10X PCR buffer, 6.25 µl of trehalose 10%, 0.06 µl of Taq DNA polymerase (5U/µl) (Platinum^® Taq, Invitrogen), and 3.5 µl of DNA template (Hajibabaei et al., 2005). The PCR amplification was carried out using the following conditions: 1 min at 94°C, followed by 5 cycles of 94˚C for 40 sec, 45˚C for 40 sec and 72˚C for 1 min, followed by 35 cycles of 94˚C for 40 sec, 51˚C for 40 sec and 72˚C for 1 min, with a final extension of 72˚C for 5 minutes. PCR products were checked for quality and length using electrophoresis on ethidium bromide stained, 2% agarose gels (E-Gel 96 Invitrogen, Carlsbad, CA) and were bidirectionally sequenced by Eurofins Genomics (USA). The sequences were edited using Codon Code Aligner v. 8.0.1 (Codon Code Corporation) and uploaded to both the Barcode of Life Data System (BOLD, www.boldsystems.org) and GenBank. In BOLD, they were labeled as Parasitic Nematode of coffee berry borer (PNBC), and in GenBank, they were assigned the accession numbers MT520707 to MT520790. 

Genetic diversity and population structure. The sequences were aligned using the ClustalW algorithm in the MEGA X v. 7.0.26 (Kumar et al., 2018) software, optimizing global similarity to ensure accuracy in phylogenetic inferences. Genetic diversity was estimated by calculating the number of segregation sites (S), number of haplotypes (h), haplotype diversity (Hd), nucleotide diversity (Π) and average number of nucleotide differences (K) using the software DnaSP v. 6.12.03 (Rozas & Librado, 2009). To test the hypothesis of neutral evolution and to understand potential demographic events, such as expansion trends, Tajima’s D (Tajima, 1989) and Fu and Li’s D and F indices were estimated (Fu & Li, 1993). Genetic differentiation among populations was analyzed using the pairwise fixation index (F_ST). An analysis of molecular variance (AMOVA) was also performed to assess the statistical significance of F_ST values (Excoffier & Lischer, 2010). Two separate AMOVA analyses were conducted to determine whether genetic differentiation between these lineages was statistically significant: i) considering the populations as a panmictic group and ii) the populations were grouped based on the 2 phylogenetic lineages identified in the maximum likelihood tree. Furthermore, the number of migrants per generation (Nm), an indirect measure of gene flow between populations, was estimated. All analyses were performed using Arlequin v. 3.5 (Excoffier & Lischer, 2010). The relationship between genetic differences (F_ST) and geographic distances among populations was analyzed using the Mantel test (Mantel, 1967). Two matrices were compared: genetic differences between population pairs and geographic distances (in kilometers). Geographic distances were calculated using QGIS v3.10.3 software (https://qgis.org). These distances were used to perform a Mantel test to assess the potential relationship between genetic differentiation and geographic separation in the context of an isolation-by-distance pattern. The statistical analysis was conducted in R (R Core Team, 2020) using the Vegan package (Oksanen et al., 2010), with 1,000 permutations to assess significance.

Table 1

Geographical locations where CBB samples infected with M. hypothenemi were collected in Mexico and Guatemala.

Municipality	Locality	Country	Code	Geographic coordinates
				Latitude	Longitude
Tapachula	Santa Lucía	Mexico	SL	15°04’42.5’’	92°13’42.5
	Finca Brasil	Mexico	BR	15°05’50.3’’	92°18’43.1’’
Cacahoatán	Salvador Urbina	Mexico	SU	15°02’26.9’’	92°12’03.9
	La Alianza	Mexico	LA	15°02’40.8’’	92°11’03.2’’
	El Zapote	Mexico	EZ	15°02’08.4’’	92°10’18.7’’
	La Unidad	Mexico	LU	15°00’31.4’’	92°09’53.2
	Dos de Mayo	Mexico	DM	15°02’55.0’’	92°09’14.5’’
	Faja de Oro	Mexico	FO	15°02’56.4’’	92°09’14.4’’
	Rosario Ixtal	Mexico	RI	15°00’37.2’’	92°09’28.0’’
	San Antonio	Mexico	SA	15°00’24.2’’	92°09’04.7’’
Unión Juárez	Santo Domingo	Mexico	SD	15°01’37.5’’	92°06’06.9’’
	San Rafael	Mexico	SR	15°02’29.8’’	92°07’10.3’’
	San Jerónimo	Mexico	SJ	15°02’26.9’’	92°08’08.9’’
	Monte Perla	Mexico	MP	15°02’42.1’’	92°05’17.8’’
	Río Suchiate	Mexico	RS	15°03’26.1’’	92°04’13.2’’
	Once de Abril	Mexico	OA	15°03’14.8’	92°08’30.9’’
Acacoyagua	Los Cacaos	Mexico	LC	15°23’23’’	92°39’13.0’’
San Pablo	BuenaVista	Guatemala	BV	14°57’53.1	91°59’48.7’’

The phylogeny of M. hypothenemi was inferred using RAxML v. 2.0. (Stamatakis, 2006). The General Time Reversible model with a proportion of invariant sites (GTR+I) was selected as the best approximation to the (TIM3+I) (-ln L 1,588.857; Delta-AIC = 0.0; AIC = 3,573.631) model obtained from the Akaike information criterion (AIC) (Akaike, 1974) analysis performed in JModelTest v. 0.1.1 (Posada, 2008). The analysis included 10,000 bootstrap replicates to ensure robust nodal support. Sequences of the nematode Howardula aoronymphium (Welch) (Tylenchida: Allantonematidae, GenBank: AY589466) was used as external group. The relationship between the identified haplotypes was analysed using a haplotype network generated in the PopART software v. 1.7 (Leigh & Bryant, 2015), applying the Median-Joining criterion.

Results

A total of 76 sequences of 560 bp from M. hypothenemi were analysed revealing 93 polymorphic sites (16.6%), whose A-T content (67.6%) was higher than the G-C content (32.4%). Global genetic diversity was moderate (Hd = 0.52 ± 0.06) (± SD) (Table 2). Ten populations presented zero diversity, while the SU and LC populations presented the highest genetic diversity from 18 populations (Table 2). The global neutrality tests were statistically significant according to Tajima values (D = -2.062, p < 0.05) and Fu and Li (D = 2.167, p < 0.02) (Table 2). Tajima and Fu and Li values were not significant inside of each population. A total of 6 haplotypes were identified, including one dominant haplotype (H1 = 50) and 4 unique ones (H3 = 2, H4 = 1, H5 = 5, and H6 = 1). The dominant haplotype was present in 16 from 18 populations sampled, while the second most common haplotype (H2 = 17) was recorded in 8 populations (Table 3). 

High genetic differentiation was observed in M. hypothenemi, with a global value F_ST of 0.66 (p < 0.05) and a low number of migrants per generation (Nm = 0.50). The EZ, BR, and LC populations showed greater genetic differentiation than the rest of the populations, with an average distance between these and the other populations of 1.7-66.5 km, 9.5-50 km, and 50.5-87.5 km, respectively (Table 4). The smallest geographical distance was between the FO and DM populations (0.04 km), with no genetic differentiation between them. The AMOVA analysis conducted under panmictic indicate that genetic variation was greater between populations (66.61%) than within populations (33.39%). A second AMOVA, grouping populations according to the 2 phylogenetic lineages identified in the maximum likelihood tree, revealed a higher level of genetic differentiation between the lineages, with an F_ST of 0.99 (p < 0.01), indicating that the differentiation between the 2 lineages was statistically significant. In this second AMOVA, genetic variation was mainly distributed between lineages (99.42%) rather than within them (0.58%). The Mantel test revealed a significant correlation between the genetic pairwise distances and the geographic distances (km) among M. hypothenemi populations (r = 0.4728, p < 0.0001), consistent with a pattern of isolation by distance. 

The tree topology revealed 2 independent clades with high support (> 90%). Clade I formed by 3 monophyletic groups, which corresponds to 96% of the specimens distributed in the municipalities of Tapachula, Cacahoatán, Unión Juárez, and San Pablo (Guatemala). Inside clade I, 2 subgroups were observed. The first subgroup had a bootstrap support value of 82%, while the second subgroup showed a higher support, with a bootstrap value of 99%. Clade II formed by 3 specimens belonging to the LC (Acacoyagua) population. A phylogenetic separation (16.5%) was observed between clade I and clade II specimens (Fig. 2).

The most frequent haplotype was H1, followed by H2. Both haplotypes were found in most of the studied sites, being particularly frequent in the municipalities of Tapachula, Cacahoatán, Unión Juárez, Acacoyagua, and San Pablo. On the other hand, H6 was exclusively detected in SU, in the municipality of Cacahoatán. H3 and H4 were restricted to LC, in Acacoyagua, while H5 was only found in BR, in Tapachula. The distribution of haplotypes by municipality and their frequencies are shown in Table 3. The inferred haplotype network identified H1 as the central, while the other haplotypes appeared as its derivatives. H4 diverged from H1 by 90 mutational steps, whereas H2, H5, and H6 diverged from H1 by only 1 mutational step (Fig. 3C).

Discussion

This study presents the first genetic description of M. hypothenemi populations, a nematode endoparasite of H. hampei, whose presence has only been naturally recorded in several coffee plantations in Mexico (Pérez et al., 2014). The genetic variability of M. hypothenemi suggests significant differentiation among its populations. Although not all populations exhibit the same degree of differentiation, the results support the hypothesis that this parasite is endemic to southeastern Mexico (Poinar et al., 2004). The genetic complexity of the populations of this parasite is similar to that observed in populations of endemic species (Blouin et al., 1999; Stock 2009). 

This study also includes the first genetic sequences of M. hypothenemi in the BOLD System database. The COI gene has not yet been widely adopted as a tool for identifying parasitic nematode species, a high genetic structure in M. hypothenemi, an obligate parasite of H. hampei. 

Table 2

Estimation of the genetic diversity in 18 populations of the nematode M. hypothenemi collected in the Soconusco region, Chiapas, Mexico.

Population	N	S	H	Su	Hd (± SD)	Π (± SD)	K	Tajima	Fu and Li
								D	D	F
SL	5	1	2	0	0.600 (0.175)	0.0010 (0.0003)	0.600	1.224	1.224	1.157
BR	5	0	1	0	0	0	0	0	0	0
SU	4	2	3	2	0.833 (0.222)	0.0017 (0.0006)	1.000	-0.709	-0.709	-0.604
LA	2	1	2	1	1 (0.500)	0.0017 (0.0008)	1	–	–	–
EZ	4	0	1	0	0	0	0	0	0	0
LU	4	1	2	1	0.500 (0.265)	0.0008 (0.0004)	0.500	-0.612	-0.612	-0.478
DM	5	1	2	1	0.400 (0.237)	0.0007 (0.0004)	0.400	-0.816	-0.816	-0.771
FO	3	1	2	1	0.666 (0.314)	0.0011 (0.0000)	0.666	–	–	–
RI	5	0	1	0	0	0	0	0	0	0
SA	5	1	2	1	0.400 (0.237)	0.0007 (0.0000)	0.400	-0.816	-0.816	-0.771
SD	2	0	1	0	0	0	0	–	–	–
SR	5	0	1	0	0	0	0	0	0	0
SJ	4	0	1	0	0	0	0	0	0	0
MP	5	0	1	0	0	0	0	0	0	0
RS	5	0	1	0	0	0	0	0	0	0
OA	4	0	1	0	0	0	0	0	0	0
LC	4	92	3	91	0.833 (0.222)	0.0824 (0.0431)	46.166	-0.836	-0.836	-0.895
BV	5	0	1	0	0	0	0	0	0	0
Total	76	93	6	0	0.519 (0.056)	0.0134 (0.0067)	7.6291	-2.0622*	2.1674**	0.5539

N, Number of sequences; S, number of segregating sites; H, number of haplotypes; Su, number of unique sites; Hd, haplotype diversity; Π, nucleotide diversity; K, average number of nucleotide differences. ±SD: Standard deviation. The D Tajima index (Tajima, 1989), D and F of Fu and Li (Fu & Li, 1993). The hyphen represents an estimate not performed due to a limited number of sequences used. * p < 0.05 ** p < 0.02.

Identifying and detecting M. hypothenemi using classical taxonomy remains challenging due to its small size and the subtlety of its morphological characteristics, requiring a fast and reliable method such as DNA barcoding (Gonçalves et al., 2021). For example, 11 new marine nematode species from the Mexican Caribbean have been identified using this technique (Martínez et al., 2020). However, the COI gene has not yet been widely used to identify parasitic nematode species. To date, only 2 species have been reported naturally parasitizing adult coffee berry borers: Panagrolaimus sp. in India (Varaprasad, 1994) and M. hypothenemi in Mexico (Castillo et al., 2002). Nevertheless, the presence of M. hypothenemi has only been detected in this region of the world, and its identification through DNA barcoding could expand our understanding of this species’ geographic range.

Molecular analyses using the COI gene identified 6 haplotypes distributed among the sampled populations, with moderate global genetic diversity (Hd = 0.52) (Wang, 2020). This level of genetic diversity, although lower than that reported for other nematodes such as Heterorhabditis bacteriophora (Hd = 0.8; Saeb & David, 2014), is notable considering the limited geographic distribution of M. hypothenemi and the environmental pressures it faces in Soconusco coffee plantations (Simota et al., 2024). However, this genetic diversity was higher compared to other studies that used 40-120 specimens per population (Klimpel et al., 2007; Powers et al., 2018), possibly due to the number of samples used per population or the high evolutionary rate of mitochondrial DNA (Avise et al., 1987; Brown et al., 1979). Therefore, our sample size may have influenced the estimation of population parameters, and a larger number of specimens per population is recommended in future studies. The observed genetic diversity reflects the adaptive capacity of populations to environmental changes or selective pressure (Clarke, 1979), key aspects for the survival of M. hypothenemi in a fragmented environment. In this sense, understanding the genetic variation of this nematode is essential for designing management and conservation strategies aimed at preserving its role as a potential biological regulator of H. hampei.

Table 3

Frequency of female haplotypes for 18 populations of the nematode Metaparasitylenchus hypothenemi collected in the Soconusco region, Chiapas, Mexico and San Pablo, Guatemala.

Municipality	Population	*N	Haplotipe frequency
Tapachula	SL	5	H1(2), H2(3)
	BR	5	H5(5)
Cacahoatán	SU	4	H1(1), H2(2), H6(1)
	LA	2	H1(1), H2(1)
	EZ	4	H2(4)
	LU	4	H1(3), H2(1)
	DM	5	H1(4), H2(1)
	FO	3	H1(2), H2(1)
	RI	5	H1(5)
	SA	5	H1(1), H2(4)
Unión Juárez	SD	2	H1(2)
	SR	5	H1(5)
	SJ	4	H1(4)
	MP	5	H1(5)
	RS	5	H1(5)
	OA	4	H1(4)
Acacoyagua	LC	4	H1(1), H3(2), H4(1)
San Pablo	BV	5	H1(5)

* Number of sequences used in the analysis; the number in parenthesis is the female specimens observed for each haplotype.

Selective neutrality analyses showed significant negative values at a global level, which commonly occur in populations undergoing demographic expansion (Schmidt, 2002), a process that likely began when this pest invaded Mexico in 1978 (Baker, 1984). Selective neutrality analyses revealed negative values for the DM, SA, and LC populations, while the SL population presented a positive value. Neutrality values within populations were not significant, indicating that the studied populations are not under the effects of natural selection (Schmidt, 2002). Therefore, it is feasible that the parasitism of H. hampei by M. hypothenemi represents a novel ecological interaction, as often occurs when a parasite is introduced into a new environment (Bush et al., 2001).

The determination of the origin of the parasite and the age of a parasite-host interaction is complex, due to the variety of mechanisms involved in the historical evolution of parasitism, including the possibility of an initial accidental association (Rico, 2011). The diversity of scolytines associated with coffee plantations in the Soconusco region is very high (Equihua, 1992), and H. hampei is endemic to Africa (Le Pelley, 1968), where populations can reach up to 11 million adults per hectare in mixed plantations of Robusta and Arabica coffee (Baker & Barrera, 1993). Under these conditions, it is likely that an endemic parasitic nematode from this region formed one or more new host-parasite associations with H. hampei (Bickford et al., 2007). Indeed, Poinar et al. (2004) hypothesized that M. hypothenemi first infected the insect in the New World.

Genetic differentiation analyses (F_ST) revealed high genetic differentiation among M. hypothenemi populations (F_ST = 0.66, p < 0.05), with gene flow that, in general, is very limited (Nm = 0.50), possibly due to the nematode’s dependence on host movement and the physical barriers imposed by a fragmented landscape (Ali et al., 2016; Fonseca & Netto, 2006; McGaughran et al., 2014). AMOVA confirmed that most genetic variation is found between populations (66.61%), indicating significant genetic structure and rejecting the hypothesis of panmixia. The observed level of genetic differentiation (FST = 0.66) suggests significant population structure, making panmixia unlikely. Additionally, the grouping of populations according to the 2 identified phylogenetic lineages revealed even greater genetic differentiation (F_ST = 0.99, p < 0.01). This high level of differentiation suggests strong genetic structuring, further reinforcing the rejection of panmixia. These findings highlight the need for further biological studies on M. hypothenemi to explore the existence of 2 lineages or even
distinct species.

Table 4

Genetic (F_ST) and geographic distances (km, above the diagonal in bold) calculated for 18 populations of the nematode M. hypothenemi collected in the Soconusco region, Chiapas, Mexico and San Pablo, Guatemala. 

	LA	SD	EZ	LU	DM	FO	SR	SJ	MP	RI	RS	OA	SL	*BV	LC	BR	SU	SA
LA	0	9.36	1.70	4.54	3.38	3.39	7.20	5.40	10.68	4.82	12.75	4.82	6.19	22.67	64.85	15.38	1.92	5.59
SD	0	0	7.84	7.29	6.27	6.29	2.54	4.07	2.50	6.49	4.86	5.37	15.20	13.59	73.49	24.65	11.14	5.94
EZ	0.38	1	0	3.10	2.45	2.48	5.86	4.05	9.36	3.22	11.55	3.91	7.90	21.01	66.55	17.03	3.30	3.95
LU	0	0	0.66	0	4.59	4.64	6.2	4.81	9.42	0.79	11.82	5.65	10.51	19.32	68.99	19.12	5.39	1.51
DM	0	0	0.72	0	0	0.04	3.91	2.20	7.33	4.28	9.36	1.48	8.92	19.83	67.34	18.39	5.30	4.67
FO	0	0	0.57	0	0	0	3.92	2.22	7.32	4.32	9.35	1.46	8.91	19.84	67.31	18.38	5.32	4.71
SR	0.47	0	1	0.06	0	0.18	0	1.81	3.49	5.50	5.74	2.85	12.80	16.11	70.96	22.30	9.07	5.25
SJ	0.38	0	1	0	0	0.11	0	0	5.31	4.18	7.51	1.62	11.13	17.63	69.50	20.59	7.26	4.16
MP	0.47	0	1	0.06	0	0.18	0	0	0	8.64	2.41	6.05	16.04	13.54	73.70	25.57	12.56	8.21
RI	0.47	0	1	0.06	0	0.18	0	0	0	0	11.04	5.18	10.93	18.61	69.50	19.70	5.89	0.82
RS	0.47	0	1	0.06	0	0.18	0	0	0	0	0	7.97	17.76	13.15	74.73	27.26	14.66	10.62
OA	0.38	0	1	0	0	0.11	0	0	0	0	0	0	10.00	18.96	68.12	19.53	6.75	5.37
SL	0	0.28	0.19	0	0.10	0	0.5	0.45	0.5	0.5	0.5	0.45	0	28.72	58.65	9.52	5.18	11.72
*BV	0.47	0	1	0.06	0	0.18	0	0	0	0	0	0	0.5	0	87.08	38.05	24.25	17.81
LC	0.52	0.52	0.66	0.66	0.70	0.60	0.70	0.66	0.70	0.70	0.70	0.66	0.70	0.70	0	50.05	63.61	70.31
BR	0.85	1	1	0.82	0.83	0.82	1	1	1	1	1	1	0.81	1	0.70	0	13.85	20.53
SU	0	0.31	0	0.14	0.24	0	0.55	0.5	0.55	0.55	0.55	0.5	0	0.55	0.66	0.78	0	6.71
SA	0	0.62	0	0.31	0.41	0.15	0.75	0.72	0.75	0.75	0.75	0.72	0	0.75	0.70	0.88	0	0

* Site located in San Pablo, Guatemala.

Pairwise F_ST values showed that some populations, such as EZ, BR, and LC, exhibit high genetic differentiation compared to other populations, despite moderate geographic distances (average of 1.7-87.5 km). This suggests that these geographic distances might be sufficient to induce isolation by distance, although the presence of specific environmental characteristics or geographic barriers could also play a role. However, the presence of one common haplotype alongside 3 highly divergent ones could indicate the existence of cryptic lineages, migration from other populations, or the persistence of ancestral genetic diversity (Blouin et al., 1995; Chávez-González et al., 2022; Nieberding et al., 2005). Further studies with a larger sample size and broader geographic coverage are needed to clarify these patterns. On the other hand, the FO and DM populations, separated by only 0.04 km, showed no genetic differentiation, suggesting a high degree of gene flow between them, likely facilitated by their geographic proximity and a more homogeneous landscape.

The high genetic differentiation detected in this study could have important implications for the conservation of M. hypothenemi. On one hand, it suggests that populations are undergoing local adaptation processes, which could result in genetically distinct lineages that might be lost if their specific habitats are not protected. The significant correlation between genetic and geographic distances (r = 0.4728, p < 0.0001) obtained through the Mantel test supports the hypothesis of a pattern of isolation by distance. However, the moderate correlation value suggests that, in addition to geographic distance, other factors such as local environmental conditions, coffee plantation management practices, or landscape fragmentation may also contribute to the observed genetic differentiation (Diniz-Filho et al., 2013).

The presence of 2 well-differentiated clades with high phylogenetic support (> 90%) suggests complex evolutionary processes. Clade I, which includes 96% of specimens distributed in the municipalities of Tapachula, Cacahoatán, Unión Juárez, and San Pablo, also showed internal substructures, indicating genetic differentiation within these populations. The formation of these subgroups could be interpreted as a possible response to geographic isolation or ecological barriers that have favored genetic divergence among populations (Wright, 1943). On the other hand, Clade II, composed only of 3 specimens from the LC population, showed a 16.4% phylogenetic separation from Clade I. This differentiation is considerably high and exceeds the intraspecific threshold generally accepted for genetic barcoding (1-2%; Lanteri, 2007). However, studies on nematodes have reported higher genetic divergences between species, as observed in Rhabdochona spp., where genetic differentiation range from 11% to 15% (Caspeta-Mandujano et al., 2021; Santacruz et al., 2020). This suggests that further comparisons with closely related taxa are needed to determine whether this divergence corresponds to an ongoing speciation process or represents intraspecific variation. The existence of 2 well-differentiated clades raises interesting questions about the evolutionary and ecological processes shaping M. hypothenemi populations. One possibility is that these differences reflect a conserved ancestral lineage in certain geographically isolated populations, such as LC. Another possibility is that the exclusive haplotypes are the result of recent local adaptation events or genetic drift in response to landscape fragmentation (Cheptou et al., 2017; Safran & Nosil, 2012). It is also plausible that these differences reflect the presence of a different species or a cryptic species with conserved morphology but genetic divergence. To confirm the occurrence of this genetic divergence process, further genetic studies with greater representation of the isolated population, as well as additional morphological studies, would be necessary. The use of nuclear markers would be particularly useful to determine whether the observed genetic differentiation is restricted to mitochondrial DNA or also occurs at the genome-wide level.

The haplotype analysis supports the hypothesis of genetic differentiation among populations. The high frequency and wide distribution of haplotype H1 found in most of the studied sites, suggest that it may represent the ancestral haplotype and the origin of other derived haplotypes. Nevertheless, the high divergence observed between haplogroups, particularly the separation of H4 by 90 mutational steps, indicates significant genetic structuring and highlights the complexity of inferring ancestral states (Avise, 2000). Therefore, it is not possible to confirm H1 as the ancestral haplotype without additional phylogenetic or coalescent analyses. In contrast, the exclusivity of H3 and H4 in the LC population, as well as H6 in SU and H5 in BR, reflects possible geographic isolation, local adaptation, and genetic drift processes (Forster, 2004).

These patterns are consistent with the theory that M. hypothenemi populations are influenced by geographic barriers that restrict gene flow and promote differentiation among localities. The search for this parasitism in the entire study region is required to determine the possibility of a greater number of M. hypothenemi haplotypes. The greater genetic difference of M. hypothenemi was observed in the population with greater geographical isolation (100 km), suggesting that more haplotypes could be located as more sites with the presence of the parasite are detected. Consequently, a complete genetic description of the M. hypothenemi populations require of the location of new sites with presence of this parasitism, whether in the study region or other coffee-growing regions of Chiapas and Central America. 

This study added 6 new geographic records for M. hypothenemi, in addition to those previously known (Pérez et al., 2015). However, we believe that the identification of new sites might be influenced by the result of random interactions caused by the transport of infested fruits with the pest and the adaptation of the parasite to local climatic conditions.

This study provides the first evidence of the American origin of this nematode. Our study also opens possibilities for future research on the genetic and geographic diversity of this species, its relationship with environmental variables. These aspects are essential for designing conservation strategies for this species and its management for the control of H. hampei, an important global pest.

Acknowledgements

Special thanks to Guadalupe Eugenia Zarza Franco for her assistance in interpreting the results, to Fernando E. Vega for the comments on the manuscript, and to Jose Higinio López Urbina for making the maps.

References

Akaike, H. (1974). A new look at the statistical model identification. IEEE Transactions on Automatic Control, 19, 716–723. https://doi.org/10.1109/TAC.1974.1100705

Ali, N., Tavoillot, J., Chapuis, E., & Mateille, T. (2016). Trend to explain the distribution of root–knot nematodes Meloidogyne spp. associated with olive trees in Morocco. Agriculture, Ecosystems and Environment, 225, 22–32. https://doi.org/10.1016/j.agee.2016.03.042 

Avise, J. C. (2000). Phylogeography: the history and formation of species (Vol. 447). Cambridge, Massachusetts, London England: Harvard University Press. 

Avise, J. C., Arnold, J., Ball, R. M., Bermingham, E., Lamb, T., Neige, J. E., Reeb, C. A., & Saunders, N. C. (1987). Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics. Annual Review of Ecology and Systematics, 18, 489–522. https://doi.org/10.1146/annurev.es.18.110187.002421 

Baker, P. S. (1984). Some aspects of the behaviour of the coffee berry borer in relation to its control in southern Mexico (Coleoptera: Scolytidae). Folia Entomológica Mexicana, 61, 9–24. 

Baker, P. S., & Barrera, J. F. (1993). A field study of a population of coffee berry borer, Hypothenemus hampei (Coleoptera; Scolytidae), in Chiapas, Mexico. Tropical Agriculture, 70, 351–355.

Bickford, D., Lohman, D. J., Sodhi, N. S., Ng, P. K., Meier, R., Winker, K., Ingram, K. K., & Das, I. (2007). Cryptic species as a window on diversity and conservation. Trends in Ecology and Evolution, 22, 148–155. https://doi.org/10.1016/j.tree.2006.11.004 

Blouin, M. S., Liu, J., & Berry, R. E. (1999). Life cycle variation and the genetic structure of nematode populations. Heredity, 83, 253–259. https://doi.org/10.1038/sj.hdy.6885420 

Blouin, M. S., Yowell, C. A., Courtney, C. H., & Dame, J. B. (1995). Host movement and the genetic structure of populations of parasitic nematodes. Genetics, 141, 1007–1014. https://doi.org/10.1093/genetics/141.3.1007 

Brown, W. M., George, M., & Wilson, A.C. (1979). Rapid evolution of animal mitochondrial DNA. Proceedings of the National Academy of Sciences, 76, 1967–1971. https://doi.org/10.1073/pnas.76.4.1967 

Bush, A. O., Fernández, J. C., Esch, G. W., Seed, J. R., & Ndez, J. C. F. (2001). Parasitism: the diversity and ecology of animal parasites. Cambridge, UK: Cambridge University Press. 

Castillo, A., Infante, F., Barrera, J. F., Carta, L., & Vega, F. E. (2002). First field report of a nematode (Tylenchida: Sphaerularioidea) attacking the coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Scolytidae) in the Americas. Journal of Invertebrate Pathology, 79, 99–202. https://doi.org/10.1016/S00222011(02)00024-1 

Castillo, A., Martínez, F., Gómez, J., Cisneros, J., & Vega, F. E. (2019). Sterility of the coffee berry borer, Hypothenemus hampei (Coleoptera: Curculionidae), caused by the nematode Metaparasitylenchus hypothenemi (Tylenchidae: Allantonematidae). Biocontrol Science and Technology, 29, 786–795. https://doi.org/10.1080/09583157.2019.1597336 

Caspeta-Mandujano, J. M., Salinas-Ocampo, J. C., Suárez-Rodríguez, R., Martínez-Ramírez, C., & Matamoros, W. A. (2021). Morphological and molecular evidence for a new rhabdochonid species, Rhabdochona (Nematoda: Rhabdochonidae), parasitizing Eugerres mexicanus (Perciformes: Gerreidae), from the Lacantún River in the Biosphere Reserve of Montes Azules, Chiapas, Mexico. Revista Mexicana de Biodiversidad, 92, 01–11. https://doi.org/10.22201/ib.20078706e.2021.92.3266 

Cháves-González, L. E., Morales-Calvo, F., Mora, J., Solano-Barquero, A., Verocai, G. G., & Rojas, A. (2022). What lies behind the curtain: Cryptic diversity in helminth parasites of human and veterinary importance. Current Research in Parasitology & Vector-Borne Diseases, 2, 100094. https://doi.org/10.1016/j.crpvbd.2022.100094 

Cheptou, P. O., Hargreaves, A. L., Bonte, D., & Jacquemyn, H. (2017). Adaptation to fragmentation: evolutionary dynamics driven by human influences. Philosophical Transactions of the Royal Society B: Biological Sciences, 372, 20160037. https://doi.org/10.1098/rstb.2016.0037 

Clarke, B. C. (1979). The evolution of genetic diversity. Proceedings of the Royal Society of London. Series B. Biological Sciences, 205, 453–474. https://doi.org/10.1098/rspb.1979.0079 

Diniz-Filho, J. A. F., Soares, T. N., Lima, J. S., Dobrovolski, R., Landeiro, V. L., Telles, M. P. D. C., & Bini, L. M. (2013). Mantel test in population genetics. Genetics and Molecular Biology, 36, 475–485. https://doi.org/10.1590/S141547572013000400002 

Equihua, M. A. (1992). Coleópteros Scolytidae atraídos a trampas NTP–80 en el Soconusco, Chiapas, México. Folia Entomológica Mexicana, 84, 55–66.

Excoffier, L., & Lischer, H. E. L. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 10, 564–567. https://doi.org/10.1111/j.
1755-0998.2010.02847.x 

Fitza, K. N., Garnas, J. R., Lombardero, M. J., Ayres, M. P., Krivak–Tetley, F. E., Ahumada, R., & Slippers, B. (2019). The global diversity of Deladenus siricidicola in native and non–native populations. Biological Control, 132, 57–65. https://doi.org/10.1016/j.biocontrol.2019.01.014 

Fonseca, G., & Netto, S. A. (2006). Shallow sublittoral benthic communities of the Laguna estuarine system, Brazilian Journal of Oceanography, 54, 41–54. 

Forster, P. (2004). Ice ages and the mitochondrial DNA chronology of human dispersals: a review. Philosophical Transactions of the Royal Society, 359, 255–264. https://doi.org/10.1098/rstb.2003.1394 

Fu, Y. X., & Li, W. H. (1993). Statistical tests of neutrality of mutations. Genetics, 133, 693–709. https://doi.org/10.1093/genetics/133.3.693 

Gonçalves, L. T., Bianchi, F. M., Deprá, M., & Calegaro-Marques, C. (2021). Barcoding a can of worms: testing cox1 performance as a DNA barcode of Nematoda. Genome, 64, 705–717. https://doi.org/10.1139/gen-2020-0140 

Hajibabaei, M., DeWaard, J. R., Ivanova, N. V., Ratnasingham, S., Dooh, R. T., Kirk, S. L. et al. (2005). Critical factors for assembling a high volume of DNA barcodes. Philosophical Transactions of the Royal Society B: Biological Sciences, 360, 1959–1967. https://doi.org/10.1098/rstb.2005.1727 

Hartshorn, J. A., Fisher, J. R., Riggins, J. J., & Stephen, F. M. (2017). Molecular identification of Deladenus proximus Bedding, 1974 (Tylenchida: Neotylenchidae), a parasite of Sirex nigricornis (Hymenoptera: Siricidae). Nematology, 19, 15–20. https://doi.org/10.1163/15685411-00003029 

Hebert, P. D., Cywinska, A., Ball, S. L., & Dewaard, J. R. (2003). Biological identifications through DNA barcodes. Proceedings of the Royal Society of London. Series B: Biological Sciences, 270, 313–321. https://doi.org/10.1098/rspb.2002.2218 

Ivanova, N. V., Dewaard, J. R., & Hebert, P. D. (2006). An inexpensive, automation-friendly protocol for recovering high-quality DNA. Molecular Ecology Notes, 6, 998–1002. https://doi.org/10.1111/j.1471-8286.2006.01428.x 

Klimpel, S., Kleinertz, S., Hanel, R., & Rückert, S. (2007). Genetic variability in Hysterothylacium aduncum, a raphidascarid nematode isolated from sprat (Sprattus sprattus) of different geographical areas of the northeastern Atlantic. Jounal Parasitology Research, 101, 1425–1430. https://doi.org/10.1007/s00436-007-06620 

Kumar S., Stecher G., Li M., Knyaz C., & Tamura K. (2018). MEGA X: molecular Evolutionary Genetics Analysis across computing platforms. Molecular Biology Evolution, 35, 1547–1549. https://doi.org/10.1093/molbev/msy096 

Lanteri, A. A. (2007). Código de barras del ADN y sus posibles aplicaciones en el campo de la entomología. Revista de la Sociedad Entomológica Argentina, 66, 15–25.

Le Pelley, R. H. (1968). Pests of coffee. London: Longmans, Green and Co.

Leigh, J. W., Bryant, D., & Nakagawa, S. (2015). POPART: full–feature software for haplotype network construction. Methods in Ecology & Evolution, 6, 1110–1116. https://doi.org/10.1111/2041-210X.12410 

Mantel, N. (1967). The detection of disease clustering and a generalized regression approach. Cancer Research, 27, 209–220.

Marsjan, P. A., & Oldenbroek, J. K. (2007). Molecular markers, a tool for exploring genetic diversity. In B. Rischkowsky, & D.Pilling (Eds.), The state of the world’s animal genetic resources for food and agriculture. Rome: Food and Agriculture Organization.

Martínez, A. A., Jesús, N. D., & Leasi, F. (2020). DNA Barcoding for delimitation of putative Mexican marine nematodes species. Diversity, 12, 107. https://doi.org/10.3390/d12030107

McGaughran, A., Morgan, K., & Sommer, R. J. (2014). Environmental variables explain genetic structure in a beetle–associated nematode. Plos One, 9, 87317. https://doi.org/10.1371/journal.pone.0087317 

Nieberding, C., Libois, R., Douady, C. J., Morand, S., & Michaux, J. R. (2005). Phylogeography of a nematode (Heligmosomoides polygyrus) in the western Palearctic region: persistence of northern cryptic populations during ice ages? Molecular Ecology, 14, 765–779. https://doi.org/10.1111/j.1365294X.2005.02440.x 

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., O’hara, R. B., Simpson, G. L., & Wagner, H. (2010). Vegan: community ecology package. R package version 1.17–4. http://CRAN.R–project.org/package=vegan 

Pérez, J., Infante, F., Poinar, G. Jr., Castillo, A., & Vega, F. E. (2015). Natural parasitism of Metaparasitylenchus hypothenemi (Tylenchida: Allantonematidae) on the coffee berry borer in Chiapas, Mexico. Biocontrol Science and Technology, 25, 608–612. https://doi.org/10.1080/09583157.2014.991910 

Poinar, G. Jr., Vega, F. E., Castillo, A., Chavéz, I. E., & Infante, F. (2004). Metaparasitylenchus hypothenemi Sp. (Nematoda: Allantonematidae), a parasite of the co-ee berry borer, Hypothenemus hampei (Curculionidae: Scolytinae). Journal of Parasitology, 90, 1106–1110. https://doi.org/10.1645/GE-3369 

Posada, D. (2008). Model Test: phylogenetic model averaging. Molecular Biology and Evolution, 25, 1253–1256. https://doi.org/10.1093/molbev/msn083 

Powers, T., Harris, T., Higgins, R., Mullin, P., & Powers, K. (2018). Discovery and identification of Meloidogyne species using COI DNA barcoding. Journal of Nematology, 50, 399–412. https://doi.org/10.21307/jofnem-2018-029 

Prosser, S., Martínez, A. A., & Elías, G. M. (2013). A new set of primers for COI amplification from freshwater microcrustaceans. Molecular Ecology Resources, 13, 151–1155. https://doi.org/10.1111/1755-0998.12132 

Rico, G. (2011). The evolution of host-parasite interactions: coevolution, sexual selection and other suggested theories. Revista UDCA Actualidad and Divulgación Científica, 14, 119–130.

R Core Team (2020). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R–project.org/ 

Rozas, J., & Librado, P. (2009). DnaSP v. 5. A software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25, 1451–1452. https://doi.org/10.1093/bioinformatics/btp187 

Saeb, A. T., & David, S. K. (2014). Comparative estimation of genetic diversity in population studies using molecular sampling and traditional sampling methods. Biomedical Informatics, 10, 347–352. https://doi.org/10.6026/9732063
0010347 

Safran, R. J., & Nosil, P. (2012). Speciation: the origin of new species. Nature Education Knowledge, 3, 17. 

Santacruz, A., Ornelas-García, C. P., & De León, G. P. P. (2020). Diversity of Rhabdochona mexicana (Nematoda: Rhabdochonidae), a parasite of Astyanax spp. (Characidae) in Mexico and Guatemala, using mitochondrial and nuclear genes, with the description of a new species. Journal of Helminthology, 94, 1–11. https://doi.org/10.1017/S0022149X19000014

Simota-Ruiz, M., Castillo, A., Cisneros-Hernández, J., & Carmona-Castro, O. (2024). Predicted distribution of Metaparasitylenchus hypothenemi (Tylenchida: Allantonematidae), parasite of the coffee berry borer. Journal of Nematology, 56, 1–10. https://doi.org/10.2478/jofnem-2024-0031 

Schmidt, D., & Pool, J. (2002). The effect of population history on the distribution of the Tajima’s D statistic. Population English Edition, 1, 1–8. 

Stamatakis A. (2006). RAxML-VI-HPC: maximum likelihood–based phylogenetic analyses with thousands of taxa and mixed models, Bioinformatics, 22, 2688–2690. https://doi.org/10.1093/bioinformatics/btl446 

Stewart, J. B., & Larsson, N. G. (2014). Keeping mtDNA in shape between generations. Plos Genetics, 10, 1–7. https://doi.org/10.1371/journal.pgen.1004670 

Stock, S. P. (2009). Molecular approaches and the taxonomy of insect–parasitic and pathogenic nematodes. In S. P. Stock, J. Vandenberg, I. Glazer, & N. Boemare (Eds.), Insect pathogens: molecular approaches and techniques (pp. 71–100). Wallingford UK: CABI. https://doi.org/10.1079/9781845934781.0071

Tajima, F. (1989). Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 123, 585–595. https://doi.org/10.1093/genetics/123.3.585 

Varaprasad, K. S., Balasubramanian, S., Diwakar, B. J., & Ramarao, C. V. (1994). First report of an entomogenous nematode, Panagrolaimus sp. from coffee berry borer, Hypothenemus hampei (Ferrari) from Karnataka, India. Plant Protection Bulletin (Faridabad), 46, 34–35.

Wang, S. Q. (2020). Genetic diversity and population structure of the endangered species Paeonia decomposita endemic to China and implications for its conservation. BMC Plant Biology, 20, 1–14. https://doi.org/10.1186/s12870-020-02682-z 

Wright, S. (1943). Isolation by distance. Genetics, 28, 114–138. https://doi.org/10.1093/genetics/28.2.114