Research Article |
Corresponding author: César Camilo Julián-Caballero ( cjulianc1500@alumno.ipn.mx ) Academic editor: Monika Lipińska
© 2024 César Camilo Julián-Caballero, Emilio Martínez-Ramírez, Rosa María Gómez-Ugalde, Eufemia Cruz-Arenas.
This 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.
Citation:
Julián-Caballero CC, Martínez-Ramírez E, Gómez-Ugalde RM, Cruz-Arenas E (2024) Trophic ecology and niche overlap of two sympatric species of Rhamdia (Siluriformes, Heptapteridae) from northeast Oaxaca, Mexico. In: Lipińska M, Lopez-Selva MM, Sierra JM (Eds) Biodiversity research in Central America. Neotropical Biology and Conservation 19(2): 117-136. https://doi.org/10.3897/neotropical.19.e119908
|
The trophic ecology of fishes in the genus Rhamdia remains poorly reported. Here, we aim to describe the diet composition and infer the trophic levels, niche breadth, and niche overlap between Rhamdia guatemalensis and Rhamdia laticauda in northeast Oaxaca, Mexico. Fishes were sampled using an electrofishing device between 2016 and 2017. We calculated the Quotient index and the diet composition was analyzed using the percentage of the Index of Relative Importance (%IRI) to analyze possible ontogenetic and sexual differences on feeding ecology. We then calculated the trophic niche breadth using the Levins´ standardized index (BA). Trophic levels (TLk) of each species were calculated using the trophic index and trophic niche overlap using Pianka´s index (Ojk). The results indicate that both species are carnivorous and tend to prey on aquatic insects, with R. guatemalensis feeding on nine different orders of animals and R. laticauda displaying a broader diet spectrum, encompassing 19 orders. Our data showed ontogeny and sex differences in the trophic ecology of R. guatemalensis. They were secondary consumers (TLk > 3) and showed a specialized diet (BA < 0.4). We found evidence of trophic niche partitioning between adults (Ojk = 0.40) and females (Ojk = 0.42). Our results highlight the fundamental role of this species in the area as a carnivorous/insectivorous, predatory fishes, secondary consumers, and biological control for Diptera, Ephemeroptera, and Megaloptera. Both species minimize the trophic overlap through food resource partitioning mediated by fish body size. Additionally, adults and females consume different prey, which allows the coexistence of closely related species. Our findings regarding trophic ecology could be crucial for devising conservation and management plans for these particular Rhamdia species.
Coexistence, diet, feeding ecology, freshwater fishes, natural history, niche breadth
Ecological differences are essential to coexistence, and closely related, sympatric (co-occurring) species need to differ in at least one of the three dimensions of the ecological niche (i.e., temporal, spatial, and trophic) to avoid competition (
The genus Rhamdia Bleeker, 1858 (Siluriformes, Heptapteridae), is found in almost every freshwater environment in the Neotropical region from southern Mexico to central Argentina (
So far, seven species of Rhamdia are known from Mexico (
Previous research in Southern Mexico has reported that R. guatemalensis exhibits a diverse and opportunistic diet that encompasses aquatic invertebrates, detritus, fishes, terrestrial arthropods, and terrestrial plant remains (
To fill these knowledge gaps, we analyzed quantitatively the overall diet composition of R. guatemalensis and R. laticauda, and compared the dietary composition between size classes and sexes. We determine the niche breadth and trophic levels to quantify the degree of dietary specialization and to determine their positions within the food web, respectively. Finally, we assess the overlap in food items between the species to examine possible feeding mechanisms that allow them to coexist. We propose that differences in the consumption of prey items could be attributed to food resource partitioning between these coexisting species.
We sampled Rhamdia individuals between October-November 2015 and February-March 2016, using an electrofishing device in 21 sites in two dryland streams, the Río Salado and Río Grande, both tributaries of the Papaloapan River within the southern TCBR, northeast of Oaxaca State, Mexico (Fig.
In the laboratory, each individual was identified using identification keys (
Prey items were sorted using a stereomicroscope (Carl Zeiss Stemi 2000-C), identified to the lowest taxonomic level possible (Suppl. materials
To characterize the diet of both species, the Intestinal Quotient (IQ) was calculated following
We analyzed the diet composition of R. guatemalensis and R. laticauda (Suppl. material
Trophic niche breadth for each species, gender, and size-class (using %N) was calculated based on Levin´s niche breadth index (B) using the RInSp package (
Trophic levels (TLk) for each species were calculated using the trophic index (
We compared the diet between the species, genders, and size classes using Pianka´s Index (
Where Ojk y Okj = the niche overlap index between the species j (R. guatemalensis) an k (R. laticauda); pij = Proportion of the resource type i relative to the diet of species j; pik = Proportion of resource type i relative to the diet of species k; and n = Total number of resource categories. This index ranges between 0 (no food items in common) and 1 (complete overlap) (
We tested for differences in SL, W, and IL between the species using a Mann-Whitney U test. All identifiable prey items were later categorized into five functional feeding groups: collectors, predators, filterers, shredders, and scrapers (
We performed null model analyses to test whether the observed niche overlap differs from what would be expected by chance (algorithm = RA3, permutation = 1000) using the EcoSimR package (
We found significant difference in SL (U = 15748; P < 0.05), W (U = 15614; P < 0.05), and IL (U = 5047; P < 0.05) between the species. The mean IQ obtained (using all specimens) for R. guatemalensis was 0.69 and 0.59 for R. laticauda, with IQ < 1 for both sexual and ontogeny groups (Table
Species-prey network consumed by two Rhamdia species. The vertices sizes are proportional to the number of interactions according to %IRI data. P1 = Diptera, P2 = Coleoptera, P3 = Ephemeroptera, P4 = Trichoptera, P5= Hymenoptera, P6 = Megaloptera, P10 = Lepidoptera, P12 = Bassommatophora, and P18 = Cyprinodontiformes. For the rest of the prey orders see Table
The diet composition of both species varies as they grow (Table
Dietary difference in %IRI between sexes were marked in R. guatemalensis, but not in R. laticauda. Female specimens of R. guatemalensis consumed a greater proportion of Megaloptera and Cyprinodontiformes compared to males (Table
We present data on morphometric and trophic niche separation between R. guatemalensis and R. laticauda, which co-occur in the TCBR, supported by evidence of body measurements (SL, W, and IL) and analysis of stomach contents data. Rhamdia guatemalensis and R. laticauda differ significantly in standard length and body mass (Table
The mean IQ values (IQ < 1) recorded for both species in all samples (Table
Sample size and range of body measurements for morphometric comparison and trophic analysis in two Rhamdia species.
Species | Samples (n) | Range of standard length (mm) | Range of weight (g) |
---|---|---|---|
Rhamdia guatemalensis | Total (n = 124) | 40.56–213.05 | 0.87–105.01 |
Males (n = 32) | 59.08–183.27 | 2.55–61.34 | |
Females (n = 28) | 49.46–213.05 | 1.62–105.01 | |
Undifferentiated (n = 13) | 40.56–70.73 | 0.87–3.38 | |
Rhamdia laticauda | Total (n = 210) | 37.80–133.40 | 0.67–33.12 |
Males (n = 57) | 57.39–126.71 | 2.59–24.02 | |
Females (n = 52) | 56.46–133.40 | 2.08–33.12 | |
Undifferentiated (n = 4) | 44.00–50.53 | 1.00–4.10 |
Composition of the diet of Rhamdia guatemalensis and Rhamdia laticauda expressed as the percentage of the Index of Relative Importance (%IIR). T = Total diet, J = Juvenile, A = Adult, F = Female, and M = Male. No = total number of specimens with food in the stomach and intestine, IQ = Intestinal Quotient, and BA = Levins´ standardized index.
Prey source | T | J | A | M | F | T | J | A | M | F |
---|---|---|---|---|---|---|---|---|---|---|
Bassommatophora | 0.86 | 0.04 | 2.06 | 0.10 | 1.90 | 0.01 | – | 0.03 | 0.05 | – |
Veneroida | – | – | – | – | – | 0.01 | – | 0.01 | 0.01 | – |
Araneae | – | – | – | – | – | 0.01 | 0.08 | – | – | 0.05 |
Trombidiformes | – | – | – | – | – | 0.02 | 0.01 | 0.03 | 0.01 | 0.04 |
Podocopida | – | – | – | – | – | 0.30 | 0.02 | 0.89 | 1.15 | 0.04 |
Decapoda | – | – | – | – | – | 0.10 | – | 0.23 | 0.21 | 0.14 |
Coleoptera | 0.70 | 0.13 | 2.63 | 1.60 | 0.69 | 1.20 | 0.93 | 1.61 | 1.23 | 1.52 |
Diptera | 42.20 | 41.36 | 32.82 | 48.90 | 28.36 | 37.7 | 37.26 | 36.92 | 35.85 | 39.47 |
Ephemeroptera | 22.30 | 46.95 | 2.77 | 20.00 | 7.26 | 50.1 | 53.70 | 47.35 | 47.63 | 50.14 |
Hemiptera | 0.03 | 0.02 | 0.05 | – | 0.22 | 0.07 | 0.19 | 0.03 | 0.02 | 0.19 |
Hymenoptera | 0.07 | – | 0.40 | 0.30 | 0.04 | 0.040 | 0.04 | 0.03 | 0.02 | 0.08 |
Lepidoptera | – | – | – | – | – | 0.017 | 0.05 | 0.01 | 0.01 | 0.04 |
Megaloptera | 29.80 | 9.39 | 49.44 | 28.60 | 48.52 | 2.40 | 2.23 | 2.76 | 3.79 | 1.76 |
Odonata | – | – | – | – | – | 1.00 | 0.18 | 2.17 | 1.92 | 0.56 |
Plecoptera | – | – | – | – | – | 0.02 | 0.06 | 0.01 | 0.01 | 0.04 |
Trichoptera | 0.47 | 2.07 | 0.16 | 0.60 | 0.71 | 6.50 | 5.25 | 7.66 | 7.68 | 5.92 |
Anura | – | – | – | – | – | 0.006 | – | 0.02 | 0.07 | – |
Cyprinodontiformes | 3.40 | 0.05 | 9.66 | 0.10 | 12.30 | 0.014 | – | 0.05 | 0.07 | – |
Characiformes | – | – | – | – | – | 0.062 | – | 0.20 | 0.31 | – |
No | 72 | 38 | 34 | 32 | 28 | 112 | 51 | 61 | 56 | 52 |
IQ | 0.69 | 0.55 | 0.59 | 0.64 | 0.79 | 0.59 | 0.41 | 0.45 | 0.58 | 0.62 |
B A | 0.22 | 0.18 | 0.26 | 0.18 | 0.39 | 0.12 | 0.15 | 0.15 | 0.15 | 0.14 |
The diet composition of R. guatemalensis and R. laticauda in this area can be interpreted as aquatic invertivores, according to
Larvae of lower Diptera vary in length from a few millimeters to many centimeters, depending on the species (
The obtained trophic level values (TLk > 3) for R. guatemalensis and R. laticauda suggest that they are secondary consumers. These results are close to the values reported for the order Siluriformes (TLk = 3.42) (
Our analysis indicates that the diet composition of both species varies depending on size, with an increase in the proportions of fishes for R. guatemalensis and fishes and decapods for R. laticauda as they grow (Table
We found a dietary difference between sexes in R. guatemalensis, but not in R. laticauda. In the case of R. guatemalensis, this difference could be related to a broader range of body measurements (Table
Values of Levins´ standardized index (BA < 0.4) suggest that both species showed a tendency towards dietary specialization (
The Pianka index indicated a low to high dietary overlap, primarily driven by the consumption of Diptera and Ephemeroptera (Fig.
These two species show overlapping distributions in the Neotropical region (
Rhamdia guatemalensis and R. laticauda differ in standard length and weight, which could avoid direct overlap in food usage. Our result highlights the fundamental role of these species in the area as carnivores and biological controls for Diptera, Ephemeroptera, and Megaloptera. Our quantitative analysis of stomach contents allows us to classify the two species as aquatic invertivores and potential predators. The diet of both species varies in composition according to size and shows sexual differentiation for R. guatemalensis. Both species are secondary consumers and tend to specialize in consuming insects. We also provide evidence that adults and females consume different prey orders, enabling closely related species to coexist through trophic partitioning. This information is essential for understanding the tropic ecology of two Rhamdia species that coexist in the TCBR and could be crucial for developing conservation and management plans for these nocturnal predator fishes in a protected natural area of Mexico.
The authors are grateful to the authorities of the Tehuacán-Cuicatlán Biosphere Reserve for providing a research permit to carry out this investigation. The Comisión Nacional para el Conocimiento y Uso de la Biodiversidad and the Instituto Politécnico Nacional (grant SIP-2015-RE/054) supported this study throughout the project LI007. In addition, the Instituto Politécnico Nacional (México) provided funding for carry out this work (grant SIP 20164780). The first author is very grateful to the Master Scholarship granted by CONAHCYT (CVU 717180). We thank Miguel Calixto-Rojas, Jose Barrios-Gutiérrez, Rubí Marcos-Gómez, Carlos Pinacho-Pinacho, Gabriel Cruz-Ruiz, David Cruz, Lourdes Hernández, and Angel García for helping in field work. We are grateful to the anonymous reviewers for their valuable comments, which have greatly improved this manuscript.
The authors have declared that no competing interests exist.
No ethical statement was reported.
The Comisión Nacional para el Conocimiento y Uso de la Biodiversidad and the Instituto Politécnico Nacional (grant SIP-2015-RE/054) supported this study throughout the project LI007. In addition, the Instituto Politécnico Nacional (México) provided funding for carry out this work (grant SIP 20164780).
Conceptualization: CCJC, ECA, RMGU, EMR. Data curation: CCJC, EMR. Formal analysis: CCJC. Funding acquisition: EMR. Investigation: CCJC, ECA. Methodology: CCJC, EMR. Project administration: EMR. Resources: EMR. Software: CCJC. Supervision: ECA, EMR, RMGU. Validation: EMR, CCJC, RMGU, ECA. Visualization: RMGU, CCJC, EMR. Writing - original draft: CCJC. Writing - review and editing: RMGU, ECA, EMR, CCJC.
César Camilo Julián-Caballero https://orcid.org/0000-0002-1995-6649
Emilio Martínez-Ramírez https://orcid.org/0000-0002-0570-5013
Rosa María Gómez-Ugalde https://orcid.org/0000-0002-1710-6763
Eufemia Cruz-Arenas https://orcid.org/0000-0003-0646-4204
All of the data that support the findings of this study are available in the main text or Supplementary Information.
List of specimens with voucher codes, locality information, and body measurements
Data type: xlsx
Explanation note: List of specimens with voucher codes (CIDOAX), locality information (MEX), and body measurements: Total length, Standard length, and Body weight using in this manuscript taken from Collection of Continental Fish at the CIIDIR Unidad Oaxaca IPN, Mexico (OAX.PEC.122.0302).
Foods identified in the digestive tracts of Rhamdia guatemalensis
Data type: xlsx
Explanation note: Foods identified in the digestive tracts of Rhamdia guatemalensis (n = 72). AP = autochthonous prey, ALP = allochthonous prey, B = benthic prey; C = prey preferring the water column. Stage L = Larva, P = Pupa, N = Nymph, AD = Adult. J = juveniles, and Ind = Indeterminate prey.
Foods identified in the digestive tracts of Rhamdia laticauda
Data type: xlsx
Explanation note: Foods identified in the digestive tracts of Rhamdia laticauda (n = 112). AP = autochthonous prey, ALP = allochthonous prey, B = benthic prey; C = prey preferring the water column. Stage L = Larva, P = Pupae, N = Nymph, AD = Adult. J = Vertebrate juveniles, and Ind = Indeterminate prey.
Percentage of ocurrence, percentage number, and percent weight of stomach content data
Data type: xlsx
Explanation note: Percentage of ocurrence, percentage number, and percent weight of stomach content data of two Rhamdia species for calculated the Index of Relative Importance.