106urn:lsid:arphahub.com:pub:dddd9632-2f62-529d-aa08-fcb37c695039Acta Ichthyologica et PiscatoriaAIeP0137-15921734-1515Pensoft Publishers10.3897/aiep.52.8469484694Research ArticleOsteichthyesTaxonomyAmericasIntraspecific morphological variation in shads, Dorosomaanale and D.petenense (Actinopterygii: Clupeiformes: Clupeidae), in the Mexican Grijalva and Usumacinta river basinsMacossay-CortezAlberto1ConceptualizationData curationFormal analysisWriting - review and editingRodiles-HernándezRocío2ConceptualizationFormal analysisFunding acquisitionProject administrationResourcesValidationWriting - original draftGonzález-DíazAlfonso A.agonzalez@ecosur.mx2ConceptualizationFormal analysisMethodologySupervisionWriting - review and editingOrnelas-GarcíaC. Patricia3MethodologySupervisionValidationWriting - original draftGonzález-AcostaAdrián F.4MethodologySupervisionValidationWriting - original draftPosgrado en Ciencias en Ecología y Desarrollo Sustentable, El Colegio de la Frontera Sur, San Cristóbal de Las Casas, Chiapas, MexicoPosgrado en Ciencias en Ecología y Desarrollo Sustentable, El Colegio de la Frontera SurSan Cristóbal de Las CasasMexicoColección de Peces, Departamento de Conservación de la Biodiversidad, El Colegio de la Frontera Sur, San Cristóbal de Las Casas, Chiapas, MexicoColección de Peces, Departamento de Conservación de la Biodiversidad, El Colegio de la Frontera SurSan Cristóbal de Las CasasMexicoColección Nacional de Peces, Departamento de Zoología, Universidad Nacional Autónoma de México, Ciudad de México, MexicoUniversidad Nacional Autónoma de MéxicoCiudad de MéxicoMexicoInstituto Politécnico Nacional, Centro Interdisciplinario de Ciencias Marinas, Instituto Politécnico Nacional, La Paz, Baja California Sur, MexicoInstituto Politécnico Nacional, Centro Interdisciplinario de Ciencias MarinasLa PazMexico
Corresponding author: Alfonso A. González-Díaz (agonzalez@ecosur.mx)
Felipe Ottoni
202207072022522149158DDE94F32-A69E-5BE5-8B32-69FF17884215696CBE49-2B29-424D-9C3A-2281CD477D933103202217062022Alberto Macossay-Cortez, Rocío Rodiles-Hernández, Alfonso A. González-Díaz, C. Patricia Ornelas-García, Adrián F. González-AcostaThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/696CBE49-2B29-424D-9C3A-2281CD477D93
Historical hydrological changes and the environmental characteristics of northern Middle America have promoted diversification and determined the distribution of fishes in the Grijalva and Usumacinta river basins of Mexico. In several taxa with wide distributions, cryptic diversity has been identified through molecular and morphological analyses. This study evaluated the intraspecific morphological variation of Dorosomaanale Meek, 1904 and Dorosomapetenense (Günther, 1867) along the Grijalva and Usumacinta river basins through geometric morphometric and linear biometric analyses. Little intraspecific differentiation was detected for either species. However, differences were identified between populations in the Grijalva basin and those from the upper Usumacinta River basins with respect to body height, head size, pelvic fin position, and anal fin size. The phenotypic expression of these attributes appears to be closely related to habitat type and geographic isolation. The morphological differences within D.petenense support the molecular hypothesis of two lineages existing in the Usumacinta River basin.
body shape variationgeometric morphometricsMiddle-American fishphenotypic differentiationFinancial support for this study was received from the project “Conectividad y diversidad funcional de la Cuenca del río Usumacinta” (Fondo de investigación Científica y Desarrollo Tecnológico de El Colegio de la Frontera Sur, FID-784), coordinated by RRH.Introduction
The highly diverse ichthyofauna of northern Middle America has a complex biogeographic history. Frequent geological, volcanic, and climatic events from the Late Cretaceous to the Miocene and Pleistocene determined the diversification and distribution of fishes in the region (Elías et al. 2021). Biogeographic patterns and a large amount of endemism among freshwater species represent the clearest evidence of the hydrological history of northern Middle America (Albert et al. 2020; Elías et al. 2021).
The Grijalva–Usumacinta hydrological system provides an excellent model for understanding the effects of geological and climatic events on the evolution of fish communities in northern Middle America. This hydrological system is characterized by its diversity of fish species and large amount of endemism, which is predominantly observed in the Cichlidae and Poeciliidae families (Gómez-González et al. 2015; Velázquez-Velázquez et al. 2016; Soria-Barreto et al. 2018; Sánchez et al. 2019; Álvarez-Pliego et al. 2021). Constant events involving vicariance resulting from the capture of rivers and changes in sea level favored the diversification of several lineages through geographic isolation (Albert et al. 2020; Elías et al. 2021). The construction of dams on the Grijalva River can be considered an example to test the recent effects of hydrological and environmental changes on the dynamics of fish populations, particularly the consequences of a reduction or interruption in genetic flow for morphological and functional differentiation (Sánchez et al. 2019).
The lower regions of the Grijalva and Usumacinta basins share many species of the same ichthyofauna (Macossay-Cortez et al. 2011; Soria-Barreto et al. 2018). Nevertheless, the structure of the community in the upper regions of these rivers notably differs (Gómez-González et al. 2015; Velázquez-Velázquez et al. 2016; Soria-Barreto et al. 2018). Studies of systematic, biogeography and population genetics have revealed that the upper regions of the two basins have independent evolutionary histories (Elías et al. 2021; Beltrán-López et al. 2021; Terán-Martínez et al. 2021). The presence of relic endemic species such as the catfish of Chiapas (Lacantuniaenigmatica Rodiles-Hernández, Hendrickson et Lundberg, 2005) and the livebearer (Xenodexiactenolepis Hubbs, 1950) of the upper Usumacinta are evidence of the region’s unique biogeographic history (Rodiles-Hernández et al. 2005; Hrbek et al. 2007; Elías et al. 2021).
Although cichlids and poeciliids are the most diverse and abundant fish families within the Grijalva–Usumacinta system, other families also reflect the effects of the region’s historic events. Such is the case of the genus Dorosoma, for which genetic evidence shows that cryptic diversity exists throughout the distribution of the species Dorosomapetenense (Günther, 1867) in Middle America, which consists of several lineages (Elías et al. 2021). Two of these lineages converge in the Grijalva–Usumacinta basins. The lineage with the broadest distribution inhabits the region from the Río Grande to the lower Grijalva and Usumacinta rivers in Mexico and Guatemala, while the other is restricted to the upper Usumacinta River and Petén Itzá Lake in Guatemala (Elías et al. 2021). Curiously, in other Dorosoma species with similar distribution patterns, no significant genetic differences have been found (Elías et al. 2021). In this sense, the migratory behavior and salinity tolerance of each species could help to explain differences in genetic segregation and morphological differentiation (Bloom and Egan 2018). Dorosomapetenense is the most widespread and tolerant of different salinity ranges (Elías et al. 2021). In contrast, Dorosomaanale Meek, 1904 is considered a freshwater fish with little tolerance for salinity changes (Castro-Aguirre et al. 1999; Elías et al. 2021).
Based on the biogeographic and molecular precedents of the ichthyofauna in northern Middle America, we proposed an analysis of the intraspecific morphological variation of the shads Dorosomaanale and D.petenense throughout the Grijalva and Usumacinta rivers in Mexico. This study used linear biometric and geometric morphometric methods. Notably, both analyses are complementary and have been widely used in ichthyology to identify intra- and interspecific morphological differences and describe patterns of variation (Kerschbaumer and Sturmbauer 2011; Tripathy 2020; Lishchenko and Jones 2021). Morphological differentiation was expected due to historical vicariance events that occurred during river separation and sea level changes in northern Middle America. Furthermore, recent biological, environmental, and hydrological changes caused by human activities were expected to affect the phenotypic expression of morphological differences.
Methods
A total of 262 adult specimens were analyzed, corresponding to the species Dorosomaanale (n = 136, 71 males and 65 females) and D.petenense (n = 126, 44 males, 82 females). The specimens were deposited in the Fish Collection of El Colegio de la Frontera Sur (ECOSC; Table 1). The specimens originated from 18 sites distributed throughout the Grijalva and Usumacinta basins in Mexico (Fig. 1 and Table 2). For the Grijalva basin, specimens from sites located in the middle region were used. For the Usumacinta basin, specimens were selected from representative sites of the upper, middle, and lower regions, following the criteria proposed by Soria-Barreto et al. (2018; see Table 2). Sexual maturity was determined by reviewing the gonads.
Sample sites location of Dorosomaanale and D.petenense in the Grijalva and Usumacinta basins: 1 = Lacantún River, 2 = Lacanjá River, 3 = San Leandro Lagoon, 4 = Miramar Lagoon, 5 = Canitzán Lagoon, 6 = Chacamax River, 7 = Nueva Esperanza Lagoon, 8 = San Pedro River, 9 = Usumacinta River, 10 = San Isidro Lagoon, 11 = Pom Lagoon, 12 = Palancares Lagoon, 13 = Vapor Lagoon, 14 = Boca Chica estuary, 15 = Malpaso dam, 16 = Chicoasén dam, 17 = Peñitas dam, 18 = Tzendales River. Black triangles indicate the dams.
https://binary.pensoft.net/fig/712062
Location of the sampling sites and samples sizes of Dorosomaanale and D.petenese in Grijalva and Usumacinta basins.
Site
Region
Coordinates
D.anale (n = 136)
D.petenense (n = 126)
1 Lacantun River
U
16°32′13′′N, 090°41′52′′W
10
1
2 Lacanja River
U
16°24′21′′N, 090°47′54′′W
10
3
3 San Leandro Lagoon
U
16°15′28′′N, 090°52′31′′W
17
10
4 Miramar Lagoon
U
16°23′40′′N, 091°15′44′′W
—
13
5 Canitzan Lagoon
M
17°35′34′′N, 091°23′46′′W
7
10
6 Chacamax River
M
17°41′08′′N, 091°41′11′′W
9
2
7 Nueva Esperanza Lagoon
M
17°47′16′′N, 091°48′30′′W
10
—
8 San Pedro River
M
16°18′48′′N, 090°53′23′′W
—
15
9 Usumacinta River
M
17°29′34′′N, 091°26′26′′W
8
—
10 San Isidro Lagoon
L
18°24′26′′N, 092°28′09′′W
10
10
11 Pom Lagoon
L
18°33′33′′N, 092°13′31′′W
10
10
12 Palancares Lagoon
L
18°33′58′′N, 092°04′33′′W
10
—
13 Vapor Lagoon
L
18°22′42′′N, 091°49′52′′W
—
10
14 Boca Chica estuary
L
18°26′46′′N, 091°47′33′′W
—
8
15 Malpaso dam
G
17°06′38′′N, 093°29′59′′W
32
12
16 Chicoasen dam
G
16°53′26′′N, 093°07′00′′W
3
9
17 Peñitas dam
G
17°27′02′′N, 093°26′03′′W
—
9
18 Tzendales River
U
16°17′20′′N, 090°54′23′′W
—
4
U = Upper, M = Middle, L = Lower, G = Grijalva.
Procrustes distances values (above diagonal) and P-values (below diagonal) to pairwise comparison test between all sections of the Grijalva–Usumacinta rivers basin to Dorosomaanale and D.petenense.
Procrustes distances
D.anale
Lower
Middle
Upper
Grijalva
Lower
—
0.0075
0.0109
0.0142
Middle
0.4219
—
0.0121
0.0136
Upper
0.064
0.0083
—
0.0164
Grijalva
0.0012
<0.0001
<0.0001
—
D.petenense
Lower
Middle
Upper
Grijalva
Lower
—
0.0111
0.0272
0.0137
Middle
0.0295
—
0.0228
0.0139
Upper
<0.0001
<0.0001
—
0.0324
Grijalva
0.0004
0.0008
<0.0001
—
Bold font indicates statistically significant P-values.
Specimens were photographed from the left side of the body with a Sony DSC-HX300 digital camera (10 megapixels) using a 10 mm reference scale. To characterize the body shape, we used a configuration of 18 fixed landmarks (Fig. 2 and Table 2) digitalized using tpsDig2 software, version 2.16 (Rohlf 2015).
Location of fixed landmarks in two species of the Dorosoma genus (image modified from Hubbs and Lagler 1947, Whitehead 1985, and Farré et al. 2016). 1 = Anterior end of the upper maxilla, 2 = End of the supraoccipital bone, 3 = Start of the dorsal fin, 4 = End of the dorsal fin, 5 = Upper boundary of the caudal fin, 6 = Center of the caudal fin, 7 = Base of the caudal fin, 8 = End of the anal fin, 9 = Origin of the anal fin, 10 = Origin of the pelvic fin, 11 = Cleitral fusion, 12 = Anterior end of the lower maxilla, 13 = Posterior end of the maxilla, 14 = Left extreme of the sphenotic orbit, 15 = Right extreme of the sphenotic orbit, 16 = Upper end of the operculum, 17 = Most posterior end at the operculum, 18 = Dorsal insertion of the pectoral fin. The dashed lines are the intraspecific linear discriminant measures.
https://binary.pensoft.net/fig/712063
Morphometric and statistical analyses. To analyze geographic intraspecific variation, specimens were classified into four groups according to the collection site. One group consisted of specimens from the Grijalva basin, while the other three were from the upper, middle, and lower regions of the Usumacinta basin (Fig. 1 and Table 2).
An analysis of geometric morphometrics was conducted using MorphoJ software, version 1.07a (Klingenberg 2011). Based on the coordinates obtained from digitization, a generalized procrustes analysis was performed (Rohlf and Slice 1990; Dryden and Mardia 1998) to eliminate variation among the configurations of the reference points due to the effects of position, orientation, and specimen size. Additionally, to eliminate the effect of allometry due to variation in specimen size, multivariate regression was conducted based on the procrustes coordinates (shape variables) with the values of the logarithm of the centroid (size variable).
With the residual values of the multivariate regression, a principal components analysis (PCA) was performed to evaluate intraspecific variation. The first two principal components were used to explore the distribution of the specimens in the morphospace and describe variation in body shape based on the deformation grids. Later, we conducted a canonical variate analysis (CVA) to determine whether significant differences in body shape exist among the four groups. Additionally, we carried out paired comparisons based on the procrustes distances. Finally, we applied a discriminant function analysis (DFA) to perform cross-validation to determine the percentage of classification of the specimens in each group based on the Mahalanobis distances. All tests subjected the data to 10,000 permutations, when appropriate.
Additionally, based on the deformation grids, we identified the body sections for which greater variation existed. We then took linear measurements to evaluate whether they are discriminant among the four groups. Measurements were obtained by using the CoordGen8 program (IMP; Sheets 2014) with the photographs. The configuration of landmarks used in the geometric morphometric analysis followed the protocols of Hubbs and Lagler (1947), Whitehead (1985), and Farré et al. (2016). Measurements were standardized according to the method of Elliott et al. (1995) to eliminate the effect of allometry caused by variation in specimen size. Later, we conducted statistical analyses (ANOVA and Kruskal–Wallis) and a posteriori tests (Tukey, Mann–Whitney) to determine whether significant differences existed among the groups. For the measurements showing significant differences, box plots were elaborated to visualize their variation, which was expressed in proportions with respect to the standard length. For statistical analyses and the elaboration of box plots, we used the PAST program, version 3.14 (Hammer et al. 2001).
Museum catalog information. Catalogue number of the specimens used in the morphometric analysis. Dorosomaanale: ECOSC 612, 658, 1286, 1737, 3492, 4426, 6708, 6714, 10713, 10714, 11748, 11752, 12549, 12555, 12665, 12790, 13521, 13533 al 13535, 13546, 13549, 13561, 13564, 13565, 13976, 14290, 14300; D.petenense: ECOSC 7339, 8698–8707, 9882–9891, 12616, 12619, 12657, 12669, 12716, 12722, 13702, 13723, 13738, 14547, 14679, 14680.
Results
Intraspecific variation in Dorosomaanale. In the PCA, the first two components explained 42.5% of the total variance. In the morphospace, no formation of groups was observed given the extensive overlap among specimens (Fig. 3A). Throughout PC1 (28.35%), greater variation was obtained for specimens from the upper region were the most broadly distributed throughout this component. Although much less variation occurred for PC2 (14.18%), it was not possible to distinguish groups (Fig. 3A). Based on the deformation grids for PC1, specimens in the negative axis were observed to have larger heads, shorter caudal peduncles, and anteriorly displaced dorsal and pelvic fins. For the positive axis, specimens had shorter heads, longer caudal peduncles, and dorsal and pelvic fins that were slightly more forward (Fig. 3A). For PC2, the most evident variations were in body height and the fact that specimens reached maximum body depth on the positive axis, while they tended to be less deep toward the negative axis (Fig. 3A).
(A) Morphospace formed by PC1 (36.20%) and PC2 (18.94%) for Dorosomaanale. Squares represent the upper region, triangles represent the middle region, dots represent the lower region and stars represent the Grijalva region. Deformation grids are associated to the most negative and positive values of the PC1 and PC2. (B) discriminatory linear measures expressed in percent for D.anale. U = Upper, M = Middle, L = Lower, G = Grijalva.
https://binary.pensoft.net/fig/712064
Meanwhile, the CVA and paired tests revealed significant differences among groups (P < 0.05). Nevertheless, the only differences among the four groups were between the Grijalva group and the three groups of the Usumacinta basin. For the Usumacinta groups, only the middle and upper regions were significantly different (P < 0.05; Table 2). Cross-validation based on the DFA indicated that the assignment percentages were highest in the Grijalva group, while the highest percentages for the Usumacinta basin were found in the group of the upper region (Table 3).
Percentage of Dorosomaanale and D.petenense correctly classified to their a priori groups based on the discriminant function analysis.
D.anale
n
Lower
Middle
Upper
Grijalva
Lower
30
72.2
56.7
76.7
83.3
Middle
34
64.7
72.6
73.5
79.4
Upper
37
70.3
70.3
77.48
91.9
Grijalva
35
82.9
85.7
88.6
85.7
Total
136
D.petenense
n
Lower
Middle
Upper
Grijalva
Lower
38
76.31
68.42
86.84
73.68
Middle
27
70.37
72.85
88.89
59.29
Upper
31
77.42
83.87
84.95
93.55
Grijalva
30
70
60
96.67
75.56
Total
126
n = number of specimens; Lower, Middle, Upper, and Grijalva are respective sections in the Grijalva–Usumacinta rivers basin.
Based on the variation in body shape observed in the deformation grids, we selected five linear measurements to evaluate their capacity to discriminate among the groups. The selected measurements were as follows: I) anterior margin of the upper mandible to the posterior margin of the operculum (landmarks 1–17); II) posterior margin of the supraoccipital crest to the anterior insertion of the anal fin (2–9); III) anterior insertion of the dorsal fin to the anterior insertion of the pelvic fin (3–10); IV) posterior insertion of the dorsal fin to the anterior insertion of the anal fin (4–9); V) anterior insertion of the anal fin to the posterior insertion of the anal fin (8–9).
In D.anale, the statistical analyses (ANOVA, Kruskal–Wallis) and respective a posteriori tests (Tukey, Mann–Whitney) revealed that only the following three measurements could discriminate at least one of the groups (P < 0.05): I) anterior margin of the upper mandible to the posterior margin of the operculum (1–17); III) anterior insertion of the dorsal fin to the anterior insertion of the pelvic fin (3–10); IV) posterior insertion of the dorsal fin to the anterior insertion of the anal fin (4–9). The box plot of these three measurements (expressed in proportions) allowed us to determine that the groups with the greatest variation were from the upper regions of the Usumacinta and Grijalva basin (Fig. 3B).
Intraspecific variation in Dorosomapetenense. In the PCA, the first two components explained 42.73% of the total variance. In PC1 (25.19%), we observed a substantial overlap of the four groups on the negative axis. However, on the positive axis, specimens from the upper Usumacinta appeared to diverge, especially from Site 4 (Miramar Lagoon, Fig. 4A). Although no clear separation was found among the groups in PC2 (17.53%), the majority of specimens from the upper Usumacinta were located toward the positive extreme (Fig. 4A). The negative end of PC1 represented specimens with relatively deep bodies and fins that were anteriorly displaced relative to specimens that possessed positive scores (Fig. 4). In the deformation grids of PC2, shallow-bodied specimens were in the negative axis, while deeper-bodied specimens were located toward the positive part (Fig. 4B).
(A) Morphospace formed by PC1 (25.1%) and PC2 (17.5%) for Dorosomapetenense. Squares represent the upper region, triangles represent the middle region, dots represent the lower region and stars represent the Grijalva region. Deformation grids are associated to the most negative and positive values of the PC1 and PC2. (B) discriminatory linear measures for D.petenense. U = Upper, M = Middle, L = Lower, G = Grijalva.
https://binary.pensoft.net/fig/712065
The CVA and paired tests showed significant differences among the four groups (P < 0.05; Table 3. The greatest procrustes distances were obtained for the upper Usumacinta group (Table 2). The DFA corroborated that the specimens of the upper Usumacinta were the most different, while the cross-validation revealed that they had the highest percentage of allocation (Table 3).
Multivariate analyses (ANOVA, Kruskal–Wallis) of the five linear measurements and the corresponding a posteriori tests showed that four measurements are discriminant (P < 0.05): II) posterior margin of the supraoccipital crest to the anterior insertion of the anal fin (2–9); III) anterior insertion of the dorsal fin to the anterior insertion of the pelvic fin (3–10); IV) posterior insertion of the dorsal fin to the anterior insertion of the anal fin (4–9); V) anterior insertion of the anal fin to the posterior insertion of the anal fin (8–9). Statistical differences were observed in the groups of the upper Usumacinta and the Grijalva. The diagrams of the four measurements (expressed in proportions) suggest that the groups that differed most were from the upper Usumacinta and the Grijalva (Fig. 4B).
Discussion
Little intraspecific morphological differentiation was observed within Dorosomaanale and D.petenense throughout their distribution in the Grijalva and Usumacinta basins. Nevertheless, patterns of variation and morphological differences were identified in some of the geographic groups, which allowed us to assume that some regional historic and/or ecological processes were involved in creating and maintaining the phenotypic differentiation in both species. In the Clupeiformes and other fish taxa, migratory behavior and tolerance to salinity could have important implications for morphological differentiation (Aguirre and Bell 2012; Bloom and Egan 2018; Jiménez-Prado and Aguirre 2021). While it was difficult to test this hypothesis based on our results, this possibility should not be discarded.
In the morphospace of D.anale, no separation among groups by geographic location was observed. However, among specimens of the upper Usumacinta basin, variation existed in terms of head size, body depth, and fin position. This could be related to the type of habitat since the separation was observed among specimens from the river (sites 1, 2, and 18) and lake habitats (sites 3 and 4). Additionally, comparisons among geographic groups based on the statistical tests showed that the specimens of the Grijalva basin were differentiated by having shallower bodies.
A similar pattern of variation was observed in the morphospace of D.petenense, while the overlap between the four geographic groups was also found. Nevertheless, specimens of the upper Usumacinta tended to be differentiated from the rest of the groups by having a deeper body. Within the upper Usumacinta group, the separation from specimens of the Miramar Lagoon was notable (Fig. 1 and Table 1). The latter had shorter anal fins and their pelvic fins were further back along the body. Statistical tests confirmed that the upper Usumacinta group was the most different; however, specimens from the Grijalva tended to exhibit slender bodies.
In both species, the greatest intraspecific morphological variation was principally observed in body shape, head size, and fin position. For many diverse species of fish, it has been demonstrated that variation in these anatomic attributes has functional importance and has been correlated with environmental factors such as water current speed, habitat structure, and the presence of predators (Langerhans 2008; Langerhans et al. 2003; Bravi et al. 2013; Araújo et al. 2014; Peris-Tamayo et al. 2020). This has also been documented in certain fishes of Middle America, such as characins (Santos and Araújo 2015; Garita-Alvarado et al. 2018, 2021), cichlids (Feilich 2016; Barrientos-Villalobos et al. 2018; Gómez-González et al. 2018; Aguilar-Contreras et al. 2021), and poeciliids (Araújo et al. 2014; Jourdan et al. 2016). In many fish taxa, phenotypic plasticity is an important precursor to morphological differentiation through adaptation to changing environmental conditions (Oufiero and Whitlow 2016).
Nevertheless, despite morphological evidence indicating that ecological-environmental factors may be promoting phenotypic differentiation in both species between the Grijalva and upper Usumacinta groups, the effect of geographic isolation and distance should also be considered—particularly for specimens of the Grijalva, which are the most geographically isolated (Sánchez et al. 2019). Furthermore, the specimens analyzed come from dams on the principal course of the Grijalva River. Notably, there is evidence of environmental changes and reductions in gene flow being caused by dam construction, which may induce rapid morphological modifications in fish (Michel et al. 2008; Franssen 2011; Perazzo et al. 2019; Gilbert et al. 2020). Additionally, while the Grijalva and Usumacinta basins share some biological and ecological characteristics, each basin has an independent biogeographical history (Elías et al. 2021). Likewise, they have been subjected to different modifications and anthropic activities (Lázaro-Vázquez et al. 2018; Sánchez et al. 2019). Also, it should be noted that hydrogeomorphological and habitat changes can affect the migratory behavior of these species—especially in the lower regions of both basins (Bloom and Egan 2018).
The morphological differentiation among D.petenense specimens in the upper Usumacinta appears to corroborate part of the hydrological history of the basin (Rosen 1967, 1970, 1979; Elías et al. 2021). Although this geographic group was the most highly differentiated within the study area, a notable separation among specimens from the Miramar Lagoon was evident in this group. This morphological evidence supports the hypothesis that diverse lineages exist within D.petenense in the Usumacinta basin (Elías et al. 2021). As presented in other biogeographic studies of the fish of northern Middle America, the territory that includes the upper Usumacinta has a particular independent evolutionary history characterized by geological events that promoted the geographic isolation of fish populations (Elías et al. 2021), such as the reversal of river current, collapse of stream beds, and underground connectivity among rivers (Rosen 1967, 1970, 1979).
Although morphological variation in the same anatomical attributes was found in both species, D.anale is less variable than D.petenense, mainly observed in the upper Usumacinta populations. Contrary to our expectations, despite being closely related and having similar ecological requirements, the magnitude and direction of morphological changes were distinct. This has also been found for other groups of fish for which the level of morphological variation among species is related to the level of dietary specialization (Ornelas-García et al. 2018; Kentao and Jearranaiprepame 2021). Regarding cichlids, generalist species were found to have greater variation in body shape than specialists (Kassam et al. 2004; Gilbert et al. 2020). It is also known that changes in head size and eye position may be related to dietary conduct, especially in terms of size and diet type (Haas et al. 2010; Kentao and Jearranaiprepame 2021).
Differences in the patterns of variation found in D.anale and D.petenense once again demonstrated that the phenotypic expression of morphological characteristics is a product of the interaction among diverse biological, environmental, and historic processes (Michel et al. 2008; Franssen et al. 2013; Bracciali et al. 2016). Thus, there is a need to continue conducting studies of the fish communities within the Grijalva and Usumacinta basins. As previously demonstrated, northern Middle America is an extraordinary natural laboratory for evaluating and understanding past and recent processes that have promoted morphological diversification among diverse species of Neotropical fish.
Conclusion
There is little intraspecific variation in body shape in D.anale and D.petenense. However, we uncovered statistically significant differences between specimens of both species from the Grijalva basins and the upper Usumacinta. Morphological differentiation was based on body height, head length, pelvic fin position, and anal fin length. Nevertheless, even though variation was observed for the same attributes in both species, the direction and magnitude differed. Since the variation in these morphological attributes seems to be related to biological, environmental, and geographic factors, it could serve to define ecotypes. For both species, morphological differences among specimens from the Grijalva basin could be due to geographic isolation. Meanwhile, differentiation among D.petenense specimens from the upper Usumacinta appears to support the hypothesis regarding the existence of two lineages in the Usumacinta basin. Additionally, in D.petenense, differentiation was detected among specimens from the Miramar Lagoon. Notably, there is a need for further taxonomic and biogeographic studies of the ichthyofauna in northern Middle America to better comprehend their diversity and the processes related to their evolution—particularly in the Grijalva and Usumacinta basins, which possess some of the most interesting and complex fish communities of the Neotropics.
Acknowledgments
We thank the anonymous reviewers who provided helpful comments that improved the manuscript. AMC thanks to the Consejo Nacional de Ciencia y Tecnología (CONACYT) for the scholarship 297018 granted for the period of 2017–2020. Financial support for this study was received from the project “Conectividad y diversidad funcional de la Cuenca del río Usumacinta” (Fondo de investigación Científica y Desarrollo Tecnológico de El Colegio de la Frontera Sur, FID-784), coordinated by RRH. AFGA thanks to COFAA and EDI_IPN Programs, and SNI-CONACYT
ReferencesAguilar-ContrerasYGonzález-DíazAAMejíaORodiles-HernándezR (2021) Morphometric variation of Middle-American cichlids: Theraps–Paraneetroplus clade (Actinopterygii: Cichliformes: Cichlidae).51(4): 403–412. https://doi.org/10.3897/aiep.51.69363AguirreWEBellMA (2012) Twenty years of body shape evolution in a threespine stickleback population adapting to a lake environment.105(4): 817–831. https://doi.org/10.1111/j.1095-8312.2011.01825.xAlbertJSTagliacolloVADagostaF (2020) Diversification of Neotropical freshwater fishes.51(1): 27–53. https://doi.org/10.1146/annurev-ecolsys-011620-031032Álvarez-PliegoNSánchezAJFloridoRSalcedoMACruz-RamírezAKBarba-MaciasE (2021) Diversidad de peces en la reserva de la biosfera Pantanos de Centla. Ecosistemas Recursos Agropecuarios 8(I): e2713. https://doi.org/10.19136/era.a8nI.2713AraújoMSPerezSIMagazoniMJPetryAC (2014) Body size and allometric shape variation in the molly Poeciliavivipara along a gradient of salinity and predation.14(1): 251. https://doi.org/10.1186/s12862-014-0251-7Barrientos-VillalobosJSchmitter-SotoJJEspinosaAJ (2018) Several subspecies or phenotypic plasticity? A geometric morphometric and molecular analysis of variability of the Mayan cichlid Mayaherosurophthalmus in the Yucatan.106(2): 268–278. https://doi.org/10.1643/CI-17-657Beltrán-LópezRGGonzález-DíazAASoria-BarretoMGarduño-SánchezMAXochitla-CastrejónCRodiles-HernándezROrnelas-GarcíaCP (2021) Genetic diversity and structure of one of the most endangered freshwater fish species in Mexico: Tlalochildebrandi (Miller, 1950) and recognition of its evolutionarily significant units. PeerJ 9: e11952. https://doi.org/10.7717/peerj.11952BloomDDEganJP (2018) Systematics of Clupeiformes and testing for ecological limits on species richness in a trans-marine/freshwater clade. Neotropical Ichthyology 16(3): e180095. https://doi.org/10.1590/1982-0224-20180095BraccialiCGuzzoGGiacomaCDeamJMSaraG (2016) Fish functional traits are affected by hydrodynamics at small spatial scale.113: 116–123. https://doi.org/10.1016/j.marenvres.2015.12.002BraviRRuffiniMScaliciM (2013) Morphological variation in riverine cyprinids: A geometric morphometric contribution.80(4): 536–546. https://doi.org/10.1080/11250003.2013.829129Castro-AguirreJLEspinosaHSchmitter-SotoJJ (1999) Ictiofauna estuarina, lagunar y vicaria de México. Limusa, Mexico.DrydenLLMardiaKV (1998) Statistical shape analysis. John Wiley Sons, Chichester. England.ElíasDJMcMahanCDPillerKR (2021) Molecular data elucidate cryptic diversity within the widespread threadfin shad (Dorosomapetenense: Clupeidae) across the Nearctic and Northern Neotropics.849(1): 89–111. https://doi.org/10.1007/s10750-021-04713-8ElliottNGHaskardKKoslowJA (1995) Morphometric analysis of orange roughy (Hoplostethusatlanticus) of the continental slope of Southern Australia.46(2): 202–220. https://doi.org/10.1111/j.1095-8649.1995.tb05962.xFarréMTusetVMMaynouFRecasensLLombarteA (2016) Selection of landmarks and semilandmarks in fishes for geometric morphometric analyses: A comparative study based on analytical methods.80(2): 175–186. https://doi.org/10.3989/scimar.04280.15AFeilichKL (2016) Correlated evolution of body and fin morphology in the cichlid fishes.70(10): 2247–2267. https://doi.org/10.1111/evo.13021FranssenNR (2011) Anthropogenic habitat alteration induces rapid morphological divergence in a native stream fish.4(6): 791–804. https://doi.org/10.1111/j.1752-4571.2011.00200.xFranssenNRStewartLKSchaeferJF (2013) Morphological divergence and flow-induced phenotypic plasticity in a native fish from anthropogenically altered stream habitats.3(14): 4648–4657. https://doi.org/10.1002/ece3.842Garita-AlvaradoCAOrnelas-GarciaCP (2021) Parallel evolution of allometric trajectories of trophic morphology between sympatric morphs of Mesoamerican Astyanax (Characidae). Applied Sciences 11(17): e8020. https://doi.org/10.3390/app11178020Garita-AlvaradoCABarluengaMOrnelas-GarcíaCP (2018) Parallel evolution of morphs of Astyanax species (Teleostei: Characidae) in México and Central America.124(4): 706–717. https://doi.org/10.1093/biolinnean/bly082GilbertMCAkamaACoxFCCraigAR (2020) Rapid morphological change in multiple cichlid ecotypes following the damming of a major clearwater river in Brazil.13(10): 2754–2771. https://doi.org/10.1111/eva.13080Gómez-GonzálezAEVelázquez-VelázquezEAnzuetoCMMaza-CruzM (2015) Fishes of the Grijalva River basin of Mexico and Guatemala.11(5): 1726. https://doi.org/10.15560/11.1726Gómez-GonzálezAEÁlvarezFMatamorosWAVelázquez-VelázquezESchmitter-SotoJJGonzález-DíazAAMcMahanCD (2018) Redescription of Viejahartwegi (Taylor & Miller 1980) (Teleostei: Cichlidae) from the Grijalva River basin, Mexico and Guatemala, with description of a rheophilic morph.4375(3): 371–391. https://doi.org/10.11646/zootaxa.4375.3.5HaasTCBlumMJHeinsDC (2010) Morphological responses of a stream fish to water impoundment.6(6): 803–806. https://doi.org/10.1098/rsbl.2010.0401HammerØHarperDATRyanPD (2001) PAST: Paleontological statistics software package for education data analysis. Paleontologia Electronica 4(1): e4.HrbekTSeckingerJMeyerA (2007) A phylogenetic and biogeographic perspective on the evolution of poeciliid fishes.43(3): 986–998. https://doi.org/10.1016/j.ympev.2006.06.009HubbsCLLaglerKF (1947) Fishes of the Great Lakes region. Bulletin of the Cranbrook Institute of Science 26: 186 pp.Jiménez-PradoPAguirreWW (2021) Variación corporal paralela en peces de dos ríos costeros del Chocó ecuatoriano.69(1): 45–59. https://doi.org/10.15517/rbt.v69i1.41814JourdanJKrauseSTLazarVMZimmerCSommerCArias-RodriguezLKlausSRieschRPlathM (2016) Shared and unique patterns of phenotypic diversification along a stream gradient in two congeneric species. Scientific Reports 6(1): e38971. https://doi.org/10.1038/srep38971KassamDMizoiriSYamaokaK (2004) Interspecific variation of body shape and sexual dimorphism in three coexisting species of the genus Petrotilapia (Teleostei: Cichlidae) from Lake Malawi.51(3): 195–201. https://doi.org/10.1007/s10228-004-0215-9KentaoAJearranaiprepameP (2021) Ecomorphological diversification of some barbs and carps (Cyprininae, Cyprinidae) in the Lower Mekong Basin of Thailand. Zoology 143: 125830. https://doi.org/10.1016/j.zool.2020.125830KerschbaumerMSturmbauerC (2011) The utility of geometric morphometrics to elucidate pathways of cichlid fish evolution.2011: 1–8. https://doi.org/10.4061/2011/290245KlingenbergCP (2011) MorphoJ: An integrated software package for geometric morphometrics.11(2): 353–357. https://doi.org/10.1111/j.1755-0998.2010.02924.xLangerhansRB (2008) Predictability of phenotypic differentiation across flow regimes in fishes.48(6): 750–768. https://doi.org/10.1093/icb/icn092LangerhansRBLaymanCALangerhansAKDewittTJ (2003) Habitat-associated morphological divergence in two Neotropical fish species.80(4): 689–698. https://doi.org/10.1111/j.1095-8312.2003.00266.xLázaro-VázquezACastilloMMJarquín-SánchezACarrilloLCappsKA (2018) Temporal changes in the hydrology and nutrient concentrations of a large tropical river: Anthropogenic influence in the Lower Grijalva River, Mexico.34(7): 649–660. https://doi.org/10.1002/rra.3301LishchenkoFJonesJB (2021) Application of shape analyses to recording structures of marine organisms for stock discrimination and taxonomic purposes. Frontiers in Marine Science 8: e667183. https://doi.org/10.3389/fmars.2021.667183Macossay-CortezASánchezAJFloridoRHuidobroLMontalvo-UrgelH (2011) Historical and environmental distribution of ichthyofauna in the tropical wetland of Pantanos de Centla, southern Gulf of Mexico.41(3): 229–245. https://doi.org/10.3750/AIP2011.41.3.11MichelCHicksBJStöltingKNClarkeACStevensMITanaRMeyerAvan den HeuvelMR (2008) Distinct migratory and non–migratory ecotypes of an endemic New Zealand eleotrid Gobiomorphuscotidianus: Implications for the incipient speciation in island freshwater fish species.8(1): 49–63. https://doi.org/10.1186/1471-2148-8-49Ornelas-GarcíaCPCórdova-TapiaFZambranoLBermúdez-GonzálezMMercado-SilvaNMendoza-GarfiasBBautistaA (2018) Trophic specialization and morphological divergence between two sympatric species in Lake Catemaco, Mexico.8(10): 4867–4875. https://doi.org/10.1002/ece3.4042OufieroCEWhitlowKR (2016) The evolution of phenotypic plasticity in fish swimming.62(5): 475–488. https://doi.org/10.1093/cz/zow084PerazzoGXCorrêaFSalzburgerWGavaA (2019) Morphological differences between an artificial lentic and adjacent lotic environment in a characid specie.29(4): 935–949. https://doi.org/10.1007/s11160-019-09582-yPeris-TamayoAMDevineauOPræbelKKahilainenKKØstbyeK (2020) A brain and a head for a different habitat: Size variation in four morphs of Arctic charr (Salvelinusalpinus (L.)) in a deep oligotrophic lake.10(20): 11335–11351. https://doi.org/10.1002/ece3.6771Rodiles-HernándezRHendricksonDALundbergJGHumphriesJM (2005) Lacantuniaenigmatica (Teleostei: Siluriformes) a new and phylogenetically puzzling freshwater fish from Mesoamerica.1000(1): 1–24. https://doi.org/10.11646/zootaxa.1000.1.1RohlfFJ (2015) The tps series of software.26(1): 9–12. https://doi.org/10.4404/hystrix-26.1-11264RohlfFJSliceDE (1990) Extensions of the Procrustes method for the optimal superimposition of landmarks.39(1): 40–50. https://doi.org/10.2307/2992207RosenDE (1967) New poeciliid fishes from Guatemala, with comments on the origins of some South and Central American forms.2303: 1–15.RosenDE (1970) A new tetragonopterine characid fish from Guatemala.2435: 1–17.RosenDE (1979) Fishes from the uplands and intermontane basins of Guatemala: Revisionary studies and comparatсive geography.162(5): 267–376.SánchezAJÁlvarez-PliegoNEspinosa-PérezHFloridoRMacossay-CortezABarbaESalcedoMAGarrido-MoraA (2019) Species richness of urban and rural fish assemblages in the Grijalva Basin floodplain, southern Gulf of Mexico.43: 239–254. https://doi.org/10.26028/cybium/2019-433-005SantosABAraújoFG (2015) Evidence of morphological differences between Astyanaxbimaculatus (Actinopterygii: Characidae) from reaches above and below dams on a tropical river.98(1): 183–191. https://doi.org/10.1007/s10641-014-0248-5SheetsHD (2014) Integrated Morphometrics Package 8: 2014. [IMP]Soria-BarretoMGonzález-DíazACastillo-DomínguezAÁlvarez-PliegoNRodiles-HernándezR (2018) Diversidad íctica en la cuenca del Usumacinta, México. Revista Mexicana de Biodiversidad 89(Suppl. 2018): 100–117. https://doi.org/10.22201/ib.20078706e.2018.0.2462Terán-MartínezJRodiles-HernándezRGarduño-SánchezMAOrnelas-GarcíaCP (2021) Molecular characterization of the common snook, Centropomusundecimalis (Bloch, 1792) in the Usumacinta basin.13(8): 347. https://doi.org/10.3390/d13080347TripathySK (2020) Significance of traditional and advanced morphometry to fishery science.1(3): 153–166. https://doi.org/10.28991/HEF-2020-01-03-05Velázquez-VelázquezELópez-VilaJMGómez-GonzálezAERomero-BernyEILievano-TrujilloJLMatamorosWA (2016) Checklist of the continental fishes of the state of Chiapas, Mexico and their distribution.632: 99–120. https://doi.org/10.3897/zookeys.632.9747WhiteheadPJP (1985) Part 1 Chirocentridae, Clupeidae and Pristigasteridae. FAO Fisheries Synopsis No. 125, 303 pp.