Next Article in Journal
Phylogenomic Insights into the Phylogeography of Halophila baillonii Asch.
Next Article in Special Issue
Are Intermittent Rivers in the Karst Mediterranean Region of the Balkans Suitable as Mayfly Habitats?
Previous Article in Journal
Topographical Based Significance of Sap-Sucking Heteropteran in European Wheat Cultivations: A Systematic Review
Previous Article in Special Issue
Trophic Structure of Macrozoobenthos in Permanent Streams in the Eastern Balkans
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 15
Issue 1
10.3390/d15010110
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity and Composition of Caddisflies (Insecta: Trichoptera) along an Elevation Gradient in Southeastern Mexico

by
María Razo-González
1,
Rodolfo Novelo-Gutiérrez
2,
Gabriela Castaño-Meneses
1,* and
Juan Márquez
3
1
Unidad Multidisciplinaria de Docencia e Investigación, Facultad de Ciencias, Campus Juriquilla, Universidad Nacional Autónoma de México, Juriquilla 76230, Mexico
2
Red de Biodiversidad y Sistemática, Instituto de Ecología A. C., Xalapa 91073, Mexico
3
Centro de Investigaciones Biológicas, Universidad Autónoma del Estado de Hidalgo, Hidalgo 42184, Mexico
*
Author to whom correspondence should be addressed.
Diversity 2023, 15(1), 110; https://doi.org/10.3390/d15010110
Submission received: 7 December 2022 / Revised: 8 January 2023 / Accepted: 10 January 2023 / Published: 13 January 2023
(This article belongs to the Special Issue Ecology, Diversity and Evolution of Aquatic Macroinvertebrates)
Browse Figures
Review Reports Versions Notes

Abstract

:
We analyzed the alpha and beta components of male Trichoptera diversity along an altitudinal gradient ranging from 136 to 2432 m a.s.l. by studying five assemblages of adult male caddisflies in five streams of the Sierra de Juárez, Oaxaca, Mexico. Samples were collected with light traps in the “nortes” (February) and rainy (October) seasons of 2016, recording air temperature and relative humidity. The assemblages showed an abundance pattern clearly associated with the rainy season, and a high dissimilarity in their species composition. Species richness and abundance were significantly higher at lower altitudinal levels. The highest alpha diversity was recorded in the Flor de Azahar stream; beta diversity increased with altitude and made an important contribution to gamma diversity. These results suggest that temporal and altitudinal variations affect environmental variables, which influence the structure of Trichoptera assemblages. Temperature was the primary factor determining spatial and temporal distribution patterns of caddisfly assemblages in the study area. However, in the rainy season, this factor explained a smaller proportion of the observed changes. Therefore, it is necessary to evaluate other environmental variables and explore other ecosystems to provide additional evidence on the factors affecting the diversity patterns of the Trichoptera in the Sierra de Juárez.

1. Introduction

In aquatic macroinvertebrate communities, caddisflies are particularly interesting as they form the most diverse order among exclusively aquatic insects. Their larvae display a wide range of feeding habits and are food for different organisms, thus contributing to trophic dynamics and energy flow between aquatic and the adjacent terrestrial systems [1,2,3]. In addition, they have been widely used in water quality biological assessment and monitoring [2,4].
Ecological diversity results from the interactions among groups of species, the relationships between species and the environment, and the evolutionary and biogeographical history of the areas and the biota inhabiting them [5]. To infer the contribution of all these elements, it is necessary to study alpha, beta, and gamma diversity [5]. Alpha diversity is defined as the variety and abundance of species at a specific study point [6] as a result of the ecological relationships in a particular area [5]. Beta diversity represents the replacement of species between communities located in different localities [7] due to the interaction between environmental, spatial, and historical factors [8]. Total beta diversity (βcc) can be divided into two elements: replacement (β-3), which involves the substitution of some species by others between assemblages, and nesting (βrich), which reflects the loss of species from the richest to the poorest assemblage [9]. Gamma diversity is defined by the species richness of a set of localities within a landscape [5].
Several studies conducted in Europe and North America have analyzed altitudinal distribution patterns of caddisflies [10,11,12,13]. Studies of this sort have recently been conducted in Central and South American countries [14,15,16,17,18], but none had been performed in Mexico to date.
In the mid-20th century, Janzen [19] noted that species richness is higher in tropical areas than in temperate areas, particularly in mountainous regions, because seasonal changes in temperature act as a selection factor for tighter thermal tolerances for organisms located along a tropical elevation gradient. Thermal gradation prevents dispersal and results in higher reproductive isolation and speciation rates [19]. This proposal prompted several studies on the subject [20,21], which generally confirmed Janzen’s prediction.
According to Jacobsen [22], the study of macroinvertebrate richness patterns along altitudinal gradients has shown variable results, and no general pattern has been detected. If such a pattern exists, it may be masked by methodological complexities, the intrinsic characteristics of the system studied, or anthropogenic effects [22]. However, temperature is one of the key factors in structuring river fauna [23]. In terms of diversity, it has been observed that alpha diversity can be positively [15] or negatively correlated [14,18] with altitude, and, in some cases, low diversity values are observed at the ends of the gradient, being highest at intermediate sites [10,12]. Beta diversity increases with altitude [18] and may have a greater contribution to gamma diversity than alpha diversity [24].
The Sierra de Juárez is a mountain range located in the north-central part of the state of Oaxaca, Mexico. Due to its geographic location and geological history, this area shows a high environmental heterogeneity, with the predominance of temperate forests, evergreen tropical forests, and mountain cloud forests. The presence of the tree Oreomunnea mexicana (Standl.) J. F. Leroy (Juglandaceae), a relict species that emerged in the late Miocene about 23 mya [25], has been recorded in the area, which is evidence of its age. This region is home to numerous species of plants and animals and several endemic species [26]. The highly uneven relief is associated with a complex hydrological network comprising a wide range of aquatic environments with the potential to host a large variety of caddisflies, among other taxa.
Although caddisflies in Mexico have been continuously studied over the past 45 years, information is still scarce due to the lack of specialists and because most studies have focused on the systematics and taxonomy of the order [27], while little work has been conducted on the ecology of caddisfly communities. The objective of the present study was to examine the diversity of caddisflies in the Sierra de Juárez eastern slope by analyzing five assemblages located at different altitudes, in relation to temporal variations and some environmental factors. According to the current knowledge on the Trichoptera, we expect that: (i) species richness will be higher at study sites located at low altitudes and lower at higher sites; (ii) alpha diversity will decrease as altitude increases, considering that the study sites are located in a tropical region; (iii) beta diversity will increase along the altitudinal gradient, and (iv) temperature will play a key role on the modeling of the assemblages studied.

2. Materials and Methods

2.1. Study Area

Sierra de Juárez is located in the north of the state of Oaxaca, in the district of Ixtlán de Juárez, delimited by the geographic coordinates 17°46′ N, 96°30′ W and 17°02′ N, 96°16′ W (Figure 1). Temperate climate predominates in the region, although there is a vast mosaic of climates on a local scale, ranging from semi-warm humid and sub-humid to semi-cold and temperate humid in the highest parts [28]. These are associated with a gradient of temperature and particular vegetation types, where pine forests predominate, followed by high evergreen tropical forests and mountain cloud forests, among others [29]. The mean annual temperature ranges between 12 °C and 18 °C, but it may drop to less than 12 °C at higher altitudes, and frost events may occur. Precipitation varies from 800 mm to 1200 mm per year, with a summer rainfall pattern. The five collection sites correspond to different streams distributed along an altitudinal gradient [22,23] ranging from 136 m a.s.l. (Soyolapam stream) to 2432 m a.s.l. (Cascadas de Comaltepec) (Table 1). These streams are located within the Papaloapan River hydrological basin (RH28) and in the subbasins of the Valle Nacional River [Soyolapam stream (SS), Bobo stream (BS), Cascadas de Comaltepec (CC)] and the Playa Vicente River [Yegu-Iag-Lagu stream (YS), Flor de Azahar stream (FA)] [30].

2.2. Field Data Collection and Insect Identification

To record the diversity of caddisflies in Sierra de Juárez, two sampling events were carried out during 2016; one during the “nortes” season (February) and the other at the end of the rainy season (October). Sampling was conducted by three persons over two consecutive days in each stream. Since many (but not all) adult caddisflies are of nocturnal or crepuscular habits and are also attracted by light, specimens were collected using two black-light traps. Each trap consisted of a 2.2 m × 1.5 m white sheet secured with ropes, illuminated with a 15-watt blacklight lamp and an energy source [32]. Light traps were placed on the banks of each of the five streams, with a minimum separation of 25 m between them, and operated for four hours from twilight. Caddisflies caught in the blacklight trap at each site were collected at 30-min intervals (considered sample units) between 18:30 and 22:30 during the winter and from 18:00 to 22:00 in the autumn. An amount of 32 samples were collected per site and 160 samples for the entire study area per season. In addition, some caddisflies were sacrificed in lethal chambers with ethyl acetate [33] to build a reference collection of pin-mounted specimens [32]. The remaining caddisflies collected were preserved in 96% alcohol bottles for subsequent identification. Similarly, air temperature and relative humidity were recorded at the same time intervals (every 30 min) using a HOBO Pro v2 thermo-hygrometer. Systematic sampling was carried out with four replicates obtained by placing two light traps on two consecutive days at each collection site. Additionally, an entomological net was used for one hour during daylight hours to capture diurnal caddisflies, such as those in the family Xiphocentronidae [34]. These were collected by sweeping the riverbank vegetation between 12:00 and 13:00 h.
For species identification, the genital structures of adult males were processed according to the technique described by Blahnik and Holzenthal [32]. The specimens collected were properly curated and deposited in the Colección Nacional de Insectos (National Insect Collection) of the Instituto de Biología (Institute of Biology) at Universidad Nacional Autónoma de México (CNIN-IBUNAM).

2.3. Data Analysis

The effects of the site and the collection season on air temperature and relative humidity were determined through a two-way ANOVA and a multiple correlation analysis to determine the direction of the relationship using the STATISTICA Ver. 6 program [35]. The overall sex ratio and the one for the most abundant species were calculated, and Chi2 values were obtained to demonstrate the significance of the differences observed. Species presence-absence matrices were constructed for each site, each season, and the entire region; these matrices were analyzed with Clench’s parametric model and the Chao2 non-parametric model to evaluate sampling efficiency and estimate the expected richness. These analyses were carried out with the programs EstimateS ver 9.1.0 [36] and STATISTICA ver. 6 [35]. Model performance was evaluated according to Jiménez-Valverde and Hortal [37] in the case of the Clench model and to González-Oreja et al. [38] for the Chao2 model. Sampling units were defined as 30-min sampling intervals. The similarities in species composition between the five assemblages and the two seasons were assessed by constructing a presence-absence matrix and a dendrogram using the unweighted pair group method with arithmetic mean (UPGMA) and performing a similarity analysis (ANOSIM) with the Jaccard index. The similarity between groups was statistically tested (p < 0.05) using 1000 permutations [39]. These analyses were carried out with the BioDiversity Pro ver. 2 program [40]. The different diversity components were analyzed by recording the number of species (richness) and estimating the “true diversity” using the exponential of the Shannon index [41] corresponding to order 1 (alpha). The unit of measurement is the effective number of species representing the diversity of a hypothetical community in which all species have exactly the same abundance [41], allowing direct comparison of the diversity of different assemblages. Only male abundances were considered for calculating diversity indices, since there is insufficient taxonomic information to classify females and assign them to the different species, because the species level taxonomy of the order is mainly based on structures of the adult male genitalia [2]. Beta diversity between assemblages and their components was calculated using the betapart package in R [42], and beta diversity between seasons was calculated using the SPADE program [43]. True gamma diversity uses multiplicative diversity partitioning proposed by Whittaker [7] and was calculated as the product of mean alpha diversity by total beta diversity [44]. A two-way ANOVA was performed to determine the effects of the site and the collection season on the total caddisfly richness and abundance, and a multiple correlation analysis to determine the direction of the relationship using the STATISTICA ver. 6 program [35]. To understand the relationship between temperature, relative humidity, and caddisfly assemblages, discriminant analyses were performed by season and taxonomic level (family, genus, and species) using the STATISTICA ver. 6 program [35]. For this and the previous analysis, the data were normalized by a √x + 0.5 correction for temperature and abundance data and an arcsine√(x/100) correction for relative humidity [45].

3. Results

3.1. Environmental Factors

The minimum and maximum air temperatures were recorded in CC (2432 m a.s.l.) during the “nortes” season (1.4 °C) and in SS (136 m a.s.l.) during the rainy season (29.2 °C), respectively. The lowest relative humidity was recorded in SS during the rainy season (68.7%) and the highest in CC during the rainy season (98%) (Table 2).
According to the two-way ANOVA, the site, the month of collection, and the interaction between both factors exerted a significant effect on temperature (site: F = 940.104, p < 0.005; month of collection: F = 2280.001, p < 0.005; interaction: F = 157.401, p < 0.005) and relative humidity (site: F = 25.204, p < 0.005; month of collection: F = 232.801, p < 0.005; interaction: F = 56.301, p < 0.005, respectively). Altitude was negatively correlated with temperature (r2 = −0.43, p < 0.001) and positively correlated with relative humidity (r2 = 0.22, p < 0.001).

3.2. Richness and Composition

From the analysis of 320 samples from the five streams, a total of 19,479 adult caddisflies were collected, belonging to 13 families, 37 genera, and 121 species (Table 3). Eighteen new species for the science were recorded: three in the genus Lepidostoma Rambur, two in Hydropsyche Pictet, two in Ochrotrichia Mosely, two in Metrichia Ross, one each in the genera Helicopsyche Siebold, Ithytrichia Eaton, Rhyacopsyche Müller, Zumatrichia Mosely, Oecetis McLachlan, Triaenodes McLachlan, and Cnodocentron Schmid, and the species Atopsyche bifurcata Razo-González and Novelo-Gutiérrez and A. jujmi Razo-González and Novelo-Gutiérrez. SS presented 71 species, while BS, YS, and FA showed lower richness, with 44, 40, and 28 species, respectively. CC showed only eight species. There was a negative and significant correlation between altitude and the number of species (r = −0.94, p < 0.005), while the correlation between altitude and abundance was not significant (r = −0.76, p = 0.14).
From the active sampling with the entomological net, we also collected the species Ithytrichia sp. in FA, Leucotrichia melleopicta Mosely in SS, and Oxyethira desadorna Moulton and Harris in YS and FA.
The geographic distribution of the species widely varied across sampling sites. At the landscape scale, 65.25% of the species were unique to a particular site (Table 4); most were recorded in SS, and ten or fewer species were unique to any other streams. No species were found in all streams, but five were found in four streams, and thirty-four (28.81%) occurred in two or three streams.
The most abundant species were Smicridea radula and Nectopsyche gemmoides. The families with the largest number of species unique to any sampling site were Hydroptilidae (20) and Hydropsychidae (10), in accordance with their high abundance. Seven species showed a wide distribution in Sierra de Juárez: Atopsyche dampfi, Hydropsyche toschiae, Hydropscyhe sp. 1, Smicridea dithyra, Oecetis disjuncta, O. metlacensis, and Wormaldia barbai. In contrast, we collected only a single specimen of the following species: Protoptila bicornuta, P. resolda, Helicopsyche dampfi, Atopsyche huenga, Leptonema championi, Smicridea dampfi, Metrichia circulatrix, Metrichia sp. 1, Metrichia. sp. 2, Neotrichia exicoma, Ochrotrichia sp. 1, Lepidostoma zapoteca, Nectosyche pavida, Marilia flexuosa, Cernotina calcea, C. chiapaneca, Polycentropus casicus, Polyplectropus carolae, Cnodocentron sp., and Xiphocentron aureum. Thirty-eight additional species were rare; i.e., between two and five specimens were collected over the entire sampling.
The overall sex ratio was significantly skewed toward females (60.24%, Chi2 = 788.36, p < 0.05). However, it was skewed toward males for some of the most abundant species, i.e., Smicridea radula (55.58%, Chi2 = 43.22, p < 0.05), S. lacanha (73.56%; Chi2 = 287.89, p < 0.05), and Oecetis disjuncta (74.24%, Chi2 = 175.96, p < 0.05). On the other hand, it was biased toward females for other abundant species, such as Nectopsyche gemmoides (85.10%, Chi2 = 1681.06, p < 0.05), S. dithyra (61.81%, Chi2 = 110.34, p < 0.05), Ochrotrichia pacifica (63.83%, Chi2 = 118.11, p < 0.05), Anchitrichia spangleri (58.81%, Chi2 = 15.86, p < 0.05), and Ochrotrichia sp. 2 (84.38%, Chi2 = 343.84, p < 0.05).
Regarding seasonality, fifty-seven species were recorded in the two seasons, only eight in the “nortes” season, and only fifty-eight in the rainy season. During the “nortes” season, 1697 caddisflies were collected, classified into 10 families, 28 genera, and 65 species; in the rainy season, 17,782 adults were collected, belonging to 13 families, 37 genera, and 115 species. At SS, 13 families, 34 genera, and 71 species were recorded; the dominant species were Nectopsyche gemmoides, Smicridea radula, and S. dithyra, representing 51.18% of the total abundance at the site. At BS, we recorded 44 species belonging to 11 families and 26 genera; Smicridea radula and Leptonema pinotepa were the dominant species, accounting for 56.58% of the abundance. At YS, forty species belonging to nine families and twenty-three genera were recorded; the dominant species were Chimarra mexicana and Leucotrichia sp., representing 38.23% of the abundance at the site. The 28 species recorded at FA belonged to 10 families and 18 genera; the dominant species were Ochrotrichia sp. 2 and Hydropsyche toschiae, accounting for 60.41% of the total. At CC, we collected fifty-four individuals classified into five families, six genera, and eight species; the dominant species were Lepidostoma sp. 3 and Wormaldia noveloi, which account for 64.81% of the total abundance at the site.
The dendrogram constructed with the Jaccard index showed that caddisfly assemblages were particular for each altitude (Figure 2). Cascadas de Comaltepec (CCll, CCn)—the site located at the highest altitude—was separated from the rest of the streams, which were clustered in two groups. On the one hand, samples of FA and YS were separated; on the other, samples of BS and YSr were separated, and, finally, the samples of SS, located at the lower part of the gradient. These results reveal that the composition of caddisfly assemblages is unique to each stream and season of the year, although the similarity analysis (ANOSIM) showed no significant differences between groups (R = 0.90, p > 0.05).
The richness estimated with the Clench and Chao2 models was very similar (133 and 134 species, respectively), so the 118 collected species represent 88.72% and 88.06% of the theoretical richness, respectively. The data showed a better fit to the Clench model (R2 ≥ 0.99, Figure 3), indicating that the inventory was still incomplete (slope greater than 0.01) as the proportion of diversity sampled at each site varied between 44.50% (FA) and 71.98% (SS) in the “nortes” season, and between 67.23% (YS) and 90.89% (SS) in the rainy season (Table 5).
In contrast to the number of species, true alpha diversity per site was 2.57-fold higher at FA (12.82 effective species) than at CC, where the lowest diversity was recorded (4.99 effective species, Figure 4).
True alpha diversity was generally lower in the “nortes” season (15.26 effective species) than in the rainy season (17.05 effective species). At each stream, the effective number of species was higher in the “nortes” than in the rainy season, except for SS and CC, which showed the opposite trend (Table 4, Figure 5). The highest true alpha diversity was recorded at FA during the “nortes” season (13.65 effective species) and the lowest at CC during the rainy season (5.05 effective species), disregarding that only a single specimen was collected at CC during the “nortes” season. True alpha diversity was similar to the expected value in most cases, except for FA in the “nortes” season (Figure 5).
Beta diversity significantly increased with altitude (r2 = 0.99, p < 0.005). The mean βcc values varied between 83% (BS and YS) and 94% (FA). The inverse relationship between taxonomic richness and altitude is characterized by the replacement of species (β-3). Dissimilarity percentages between assemblages ranged from 65% to 100%, while nesting explained between 3% (SS) and 13% (YS) of the total beta diversity. The most dissimilar assemblages were SS and BS relative to CC (Table 6).
Altitudinally closer assemblages were more similar; however, a significant replacement of species was recorded mainly between the CC and the rest of the assemblages in both seasons (Table 7).
The mean alpha diversity was 8.57 effective species, and beta diversity was 2.01 effective communities; thus, the gamma diversity of caddisflies in the Sierra de Juárez was 17.22 effective species.
The two-way ANOVA showed a significant effect of the site and the collection season on total caddisfly richness and abundance, and on the abundance of dominant species (Table 8). The season showed a positive correlation with species richness, total abundance, and abundances of the dominant species, while the site showed a negative correlation with all of the above variables (Table 9).
The discriminant analyses performed on families, genera, and species in relation to environmental variables and seasonality showed a clear differentiation between sites mainly determined by temperature (Figure 6). In the “nortes” season, the assemblages of families (Figure 6a), genera (Figure 6c), and species (Figure 6e) in CC are separated from the groups YS + FA and SS + BS. In the rainy season, assemblages are divided into two groups, CC + FA (sites at the highest altitudes) and SS + BS + YS (Figure 6b,d,f). For assemblages of families and genera in the “nortes” season, temperature was more strongly correlated with axis 1 (90.24% of the variance, eigenvalue = 29.90, and 86.26% of the variance, eigenvalue = 32.39, respectively), while relative humidity showed a stronger correlation with axis 2 (5.16% of the variance, eigenvalue = 1.71 and 6.47% of the variance, eigenvalue = 2.43, respectively). In contrast, for the species assemblage in the “nortes” season and for all assemblages in the rainy season, both temperature and relative humidity were more strongly correlated with axis 1, explaining between 74.28% (genera in the rainy season) and 88.64% (species in the “nortes” season) of the variance.

4. Discussion

This study evaluated the alpha, beta, and gamma components of male caddisfly diversity in relation to air temperature and relative humidity along an altitudinal gradient in Sierra de Juárez, Oaxaca. As far as we know, this is the first study of its kind conducted in Mexico and provides evidence to elucidate the ecological factors that shape caddisfly assemblages in mountain streams.
As expected, air temperature and relative humidity varied with altitude. Our results showed that temperature decreases and relative humidity increases with altitude [13,15].
These results suggest that caddisfly richness and diversity in Sierra de Juárez are high, as the 121 species reported here represent 20.86% of caddisfly diversity in Mexico and 58.17% in Oaxaca [46,47,48,49]. It should be noted that we collected 51 species (33.11%) that had not previously been reported for the state. Regarding estimates of richness by site, the results suggest that between 9.11% and 55.5% of the Trichoptera remains to be discovered, and a more extensive sampling in such a heterogeneous region as the Sierra de Juárez may significantly increase the number of species recorded. Other taxa in the area for which high richness has been reported are papilionid lepidopterans [50], some angiosperm groups [51], Crassulaceae of the genus Sedum [52], and mosses [53]. The above demonstrates the high species richness of caddisflies in Mexico and the desirability for further studies focused on various aspects of this group.
Hydropsychidae and Leptoceridae were dominant taxa, confirming that these are the most abundant families in tropical regions [15,54]. The observed changes in sex ratio may reflect the sex ratio of larvae in streams or may be due to the sex selectivity of the blacklight trap [55]. This study did not record data on larvae, so we cannot determine the cause.
The results showed a clear variation of caddisfly assemblages related to the altitudinal gradient, thus supporting our first hypothesis that Trichoptera richness would decrease with increasing altitude [13,18,22]; our results are in contrast with the findings of Malicky and Chantaramongkol [10] and Kučinić et al. [12], who recorded a higher richness concentrated in intermediate sites along the altitudinal range. Similarly, the dendrogram showed that the assemblages studied, although segregated according to altitude, were not structured along the spatial gradient, which may be due to the fact that other factors, such as the dispersal capacity, might have a greater influence on the structuring of assemblages [56]. In addition, the similarity analysis by Wigger et al. [13] showed the segregation of sites with significant differences, which was not observed here.
In contrast to the second hypothesis, alpha diversity showed an oscillatory trend, with the highest diversity in FA (Figure 3), one of the highest-altitude sites (1967 m a.s.l.). With regards to temporal variation, it was observed that alpha diversity increased with altitude in the “nortes” season, except for CC, where a single specimen was collected. In contrast, alpha diversity was higher at intermediate altitudes in the rainy season, as reported for a Thai stream running across an altitudinal gradient similar to the one studied here [10]. Species-rich assemblages such as SS, BS, and YS (J’ = 0.18, J’ = 0.16, and J’ = 0.28, respectively) had lower diversity values—mainly driven by the presence of highly dominant species—than FA, which showed a lower richness and a more even distribution of species abundances (J’ = 0.44). The best example of the first case was SS, where 38.64% of the individuals collected corresponded to the species Nectopsyche gemmoides and Smicridea radula.
The assemblages observed at the lowest elevations were shaped similarly to the SS assemblage. The CC assemblage showed a low richness, but it is highly important for the regional fauna because most of its species are not shared with the other assemblages studied and are most likely endemic to Mexico. This locality is situated in the highest part of the subbasin, forming part of a stream head. Despite its low alpha diversity, CC is critical to regional diversity because of its unique spatial characteristics [57]. In several works carried out on stream headwaters [58,59,60,61], up to 50% of all the taxa are unique to a particular site, as also observed in this study.
As for beta diversity, the results are also consistent with those reported by Jacobsen [62] and Castro et al. [18], where beta diversity increases with altitude in terms of greater species turnover between the lower and the higher sites of the gradient. The greatest contribution of beta diversity to the regional diversity of the Sierra de Juárez represents an addition of species to the community as a result of the environmental variation between sites, a pattern recorded in other tropical mountains [16,18,62] and vascular plant assemblages [63]. Differences in vegetation along the altitudinal gradient may translate into habitat complexity since the structure of streams is more complex at lower altitudes, allowing for the presence of more diverse assemblages [64]. Similarly, the higher diversity observed during the rainy season can be explained by the greater availability of resources at that time.
Beta diversity significantly contributed to gamma diversity, maybe due to the limited dispersal capacity of caddisflies [65,66], so the distance between study sites may have caused the variation in assemblage composition and a high level of endemism [19]. The overall beta diversity values (βcc) significantly contributed to the replacement of species (β-3), which, at the same time, is an expression of the number of unique species of a site. In contrast, the dissimilarity due to the species richness in each assemblage had a smaller contribution relative to species replacement. The pattern mentioned above is not new, as Hernández-Rodríguez and Aguirre-Hidalgo [53], when studying bryophyte assemblages in Santiago Comaltepec, reported similar contributions for the β-3 and βrich elements of beta diversity. Other studies have theorized that the high beta diversity at high-altitude sites may be a consequence of the lower dispersal capacity of insects caused by ecological and historical factors, including cold weather, wind, or impassable topographic barriers [57,62,65]; also, beta diversity increases as we approach the tropics [9].
Our results indicate that temperature is a driver of the spatial distribution of caddisfly assemblages, since the discriminant analysis showed that this parameter has a strong influence on them, regardless of the taxonomic level; furthermore, it explains a large percentage of the observed variance (more than 75.25% in all cases). In this sense, it has been argued that high temperatures favor high mutation rates and shorter generations, thus fostering evolution and the speciation process [23,67]. This pattern is dissimilar between seasons, probably because the factors associated with temperature that have the greatest influence on assemblages differ in each season; however, the segregation of sites according to temperature and altitude is evident. These results are consistent with several studies showing the preponderance of temperature over other parameters, such as the composition of the vegetation cover, the percentage of shade, and organic matter content in deep sediments, among others, in structuring caddisflies assemblages [13,68,69].
Additional important aspects shaping assemblages are historical factors. Of these, biogeographical history apparently has great influence because the Sierra de Juárez is located in the Zona de Transición Mexicana (Mexican Transition Zone), which in turn is located in the transition between the biogeographical provinces of Sierra Madre del Sur, Veracruzana, and Cuenca del Balsas [70]. A reflection of its biogeographical complexity are the eight generalized tracks crossing the region [71] and integrating species from different biotic components. The location of Sierra de Juárez enhances the expansion of the distribution ranges of species and promotes speciation events [72], which explains to some extent the high beta diversity recorded.
To improve the current knowledge of Trichoptera diversity in the Sierra de Juárez and the factors that shape its assemblages, it will be necessary to explore other habitats within the region and incorporate additional environmental factors into the analysis, such as water conditions and physicochemical parameters in relation to immature stage assemblages, which may be directly or indirectly defining diversity patterns at different scales, as suggested by other studies [8,18]. The diversity of adult caddisfly assemblages should also be analyzed in the dry season, as well as the diversity of juvenile assemblages. Nevertheless, the results reported here provide a solid argument for promoting protection and management actions that support the conservation of the streams in the Sierra de Juárez and their enormous biodiversity.
Trichoptera is one of the most diverse orders of aquatic insects and is widely distributed in both tropical and temperate environments. However, its diversity patterns have not been studied in depth in Mexico, and little is known about its geographic distribution patterns and ecology. The results reported here highlight the importance of these studies, not only for contributing to the species inventory of the Sierra de Juárez, but also for the information they provide about the importance of conserving the aquatic ecosystems in the region. Furthermore, in order to understand the actual diversity of caddisflies at the species level, considerable taxonomic research is needed to associate and describe the females and larvae of species currently only identified as males.

Author Contributions

Conceptualization and methodology, R.N.-G., G.C.-M., J.M. and M.R.-G.; data curation, software and formal analysis, G.C.-M. and M.R.-G.; sample collection, R.N.-G. and M.R.-G.; writing—original draft preparation, M.R.-G.; review and editing, R.N.-G., G.C.-M. and J.M.; project administration, R.N.-G.; funding acquisition, R.N.-G. and M.R.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Consejo Nacional de Ciencia y Tecnología (Mexico’s National Council of Science and Technology; CONACyT) through the Posgrado en Ciencias Biológicas at Universidad Nacional Autónoma de México. MRG was supported by grant number 326085.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

MRG. thanks José Antonio Gómez-Anaya (Instituto de Ecología, A. C.) for his support during fieldwork; José Arturo Casasola-González (Instituto de Estudios Ambientales, Universidad de la Sierra Juárez); José G. Palacios-Vargas (Facultad de Ciencias, Universidad Nacional Autónoma de México), Alicia Callejas-Chavero (Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional) and Alicia Rojas-Ascencio (Facultad de Ciencias, Universidad Nacional Autónoma de México) for facilitating laboratory work; Joaquín Bueno-Soria and Rafael Barba-Álvarez (Instituto de Biología, Universidad Nacional Autónoma de México) for their assistance in the identification of Trichoptera. Thanks also to the authorities of Santiago Comaltepec, San Martin Soyolapam, San Juan Tepanzacoalco, La Esperanza, Santiago Cuasimulco, and Santiago Laxopa for providing access to their land. To Posgrado en Ciencias Biológicas of Universidad Nacional Autónoma de México for the financial support during the Ph. D. studies of the first author and to CONACyT for the doctoral fellowship 326085. To INECOL A. C. for logistical support. María Elena Sánchez-Salazar translated the manuscript into English. Abel Ibáñez Huerta made the map show collecting sites. Suggestions of two anonymous reviewers improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wiggins, G.B. Larvae of the North America Caddisfly Genera (Trichoptera), 2nd ed.; University of Toronto Press: Toronto, ON, Canada, 1996; pp. 3–25. [Google Scholar]
  2. Holzenthal, R.W.; Blahnik, R.J.; Prather, A.L.; Kjer, K.M. Order Trichoptera Kirby, 1813 (Insecta), Caddisflies. Zootaxa 2007, 1668, 639–698. [Google Scholar] [CrossRef]
  3. Morse, J.C.; Frandsen, P.B.; Graf, W.; Thomas, J.A. Diversity and ecosystem services of Trichoptera. Insects 2019, 10, 125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ward, J.V. Aquatic Insect Ecology. 1. Biology and Habitat, 1st ed.; John Wiley and Sons, Inc.: New York, NY, USA, 1992; pp. 25–28. [Google Scholar]
  5. Halffter, G.; Moreno, C.E. Significado biológico de las diversidad alfa, beta y gamma. In Sobre Diversidad Biológica: El Significado de las Diversidades Alfa, Beta y Gamma; Halffter, G., Soberón, J., Koleff, P., Melic, A., Eds.; Zaragoza: m3m-Monografías Tercer Milenio; Sociedad Entomológica Aragonesa, Comisión Nacional para el uso y Conocimiento de la Biodiversidad, Grupo Diversitas y Consejo Nacional de Ciencia y Tecnología: Zaragoza, Spain, 2005; Volume 4, pp. 5–18. [Google Scholar]
  6. Maurer, B.A.; McGill, B.J. Measurement of species diversity. In Biological Diversity Frontiers in Measurement and Assessment, 1st ed.; Magurran, A.E., McGill, B.J., Eds.; Oxford University Press Inc.: New York, NY, USA, 2011; pp. 55–65. [Google Scholar]
  7. Whittaker, R.H. Vegetation of the Siskiyou Mountains, Oregon and California. Ecol. Monogr. 1960, 30, 279–338. [Google Scholar] [CrossRef]
  8. Breda, M.; Restello, R.M.; Giovenardi, R.; Vizzotto, A.P.; Soares, B.; Hepp, L.U. Alpha and beta diversities of Trichoptera (Insecta) assemblages in natural and rural subtropical streams. Acta Limnol. Bras. 2020, 32, e14. [Google Scholar] [CrossRef]
  9. Baselga, A.; Gómez-Rodríguez, C. Diversidad alfa, beta y gamma: ¿cómo medimos las diferencias entre comunidades biológicas? NACC 2019, 26, 39–45. [Google Scholar]
  10. Malicky, H.; Chantaramongkol, P. The altitudinal distribution of Trichoptera species in Mae Klang catchment on Doi Inthanon, northern Thailand: Stream zonation and cool- and warm-adapted groups. Rev. Hydrobiol. Trop. 1993, 26, 279–291. [Google Scholar]
  11. Ruiz-García, A.; Herrera-Grao, A.F.; Ferreras-Romero, M. Distribution of Trichoptera communities in the Hozgarganta catchment (Los Alcornocales Natural Park, SW Spain). Int. Rev. Hydrobiol. 2006, 91, 71–85. [Google Scholar] [CrossRef]
  12. Kučinić, M.; Vučković, I.; Kutnjak, H.; Jelaska, L.S.; Marguš, D. Diversity, distribution, ecology and biogeography of caddisflies (Insecta: Trichoptera) in the Krka river (National Park “Krka”, Croatia). Zoosymposia 2011, 5, 255–268. [Google Scholar] [CrossRef] [Green Version]
  13. Wigger, F.W.; Schmidlin, L.; Nagel, P.; von Fumetti, S. Macroinvertebrate assemblages of natural spring along an altitudinal gradient in the Bernese Alps, Switzerland. Ann. Limnol. 2015, 51, 237–247. [Google Scholar] [CrossRef] [Green Version]
  14. Miserendino, M.L. Distribución altitudinal de especies de Trichoptera en un sistema fluvial en Patagonia. Ecol. Austral 1999, 9, 28–34. [Google Scholar]
  15. Henriques-Oliveira, A.L.; Nessimian, J.L. Aquatic macroinvertebrate diversity and composition in streams along an altitudinal gradient in Southeastern Brazil. Biota Neotrop. 2010, 10, 115–128. [Google Scholar] [CrossRef] [Green Version]
  16. Scheibler, E.E.; Claps, M.C.; Roig-Juñent, S.A. Temporal and altitudinal variations in benthic macroinvertebrates assemblages in an Andean river basin of Argentina. J. Limnol. 2014, 73, 76–92. [Google Scholar] [CrossRef] [Green Version]
  17. Henriques-Oliveira, A.L.; Nessimian, J.L.; Baptista, D. Diversity and composition of Trichoptera (Insecta) larvae assemblages in streams with different environmental conditions at Serra da Bocaina, Southeastern Brazil. Acta Limnol. Bras. 2015, 27, 394–410. [Google Scholar] [CrossRef]
  18. Castro, D.M.; Callisto, M.; Solar, R.R.C.; Macedo, D.R.; Fernades, G.W. Beta diversity of aquatic invertebrates increases along an altitudinal gradient in a Neotropical mountain. Biotropica 2019, 51, 399–411. [Google Scholar] [CrossRef]
  19. Janzen, D.H. Why mountain passes are higher in the tropics. Am. Nat. 1967, 101, 233–249. [Google Scholar] [CrossRef]
  20. Graham, C.H.; Carnaval, A.C.; Cadena, C.D.; Zamudio, K.R.; Roberts, T.E.; Parra, J.L.; McCain, C.M.; Bowie, R.C.K.; Moritz, C.; Baines, S.B.; et al. The origin and maintenance of montane diversity: Integrating evolutionary and ecological processes. Ecography 2014, 37, 711–719. [Google Scholar] [CrossRef]
  21. Polato, N.R.; Gill, B.A.; Shah, A.A.; Gray, M.M.; Casner, K.L.; Barthelet, A.; Messer, P.W.; Simmons, M.P.; Guayasamin, J.M.; Encalada, A.C.; et al. Narrow thermal tolerance and low dispersal drive higher speciation in tropical mountains. Proc. Natl. Acad. Sci. USA 2018, 115, 12471–12476. [Google Scholar] [CrossRef] [Green Version]
  22. Jacobsen, D. Contrasting patters in local and zonal family richness of stream invertebrates along an Andean altitudinal gradient. Freshw. Biol. 2004, 49, 1293–1305. [Google Scholar] [CrossRef]
  23. Jacobsen, D.; Schultz, R.; Encalada, A. Structure and diversity of stream invertebrate assemblages: The influence of temperature with altitude and latitude. Freshw. Biol. 1997, 38, 247–261. [Google Scholar] [CrossRef] [Green Version]
  24. Magliozzi, C.; Meyer, A.; Usseglio, P.; Robertson, A.; Grabowski, R.C. Investigating invertebrate biodiversity around large Wood: Taxonomic vs functional metrics. Aquat. Sci. 2020, 82, 69. [Google Scholar] [CrossRef]
  25. Rzedowski, J.; Palacios, R. El bosque de Engelhardtia (Oreomunnea) mexicana en la región de la Chinantla, Oaxaca, México, una reliquia del Cenozoico. Bol. Soc. Bot. México 1977, 36, 93–123. [Google Scholar] [CrossRef] [Green Version]
  26. Morrone, J.J. Biogeographic regionalization of the Sierra Madre del Sur province, Mexico. Rev. Mex. Biodivers. 2017, 88, 710–714. [Google Scholar] [CrossRef]
  27. Bueno-Soria, J.; Barba-Alvarez, R. New species of Plectropsyche Ross 1947 (Trichoptera: Hydropsychidae: Hydropsychinae). Zootaxa 2015, 4040, 421–432. [Google Scholar] [CrossRef]
  28. Trejo, I. Clima. In Biodiversidad de Oaxaca; García-Mendoza, A.J., Ordoñez, M.J., Briones-Salas, M., Eds.; Instituto de Biología, UNAM-Fondo Oaxaqueño para la Conservación de la Naturaleza-World Wildlife Fund: México City, México, 2004; Volume 1, pp. 67–85. [Google Scholar]
  29. Arriaga, L.; Espinoza-Rodríguez, J.M.; Aguilar-Zúñiga, C.; Martínez-Romero, E.; Gómez-Mendoza, L.; Loa, E. (Eds.) Regiones Terrestres Prioritarias de México; Comisión Nacional para el Conocimiento y Uso de la Biodiversidad: México City, México, 2000. [Google Scholar]
  30. Instituto Nacional de Estadística, Geografía e Informática (INEGI). Hidrología. 2020. Available online: https://www.inegi.org.mx/temas/hidrologia/ (accessed on 10 October 2021).
  31. García, E. Modificaciones al Sistema Climático de Köppen, 5th ed.; Instituto de Geografía-Universidad Nacional Autónoma de México: Ciudad de México, Mexico, 2004; pp. 19–74. [Google Scholar]
  32. Blahnik, R.J.; Holzenthal, R.W. Collection and curation of Trichoptera, with emphasis on pinned material. Nectopsyche Neotrop. Trichoptera Newsl. 2004, 1, 8–22. [Google Scholar]
  33. Triplehorn, C.A.; Johnson, N.F. Borror and DeLong’s Introduction to the Study of Insects, 7th ed.; Thomson Brooks/Cole: Belmont, MA, USA, 2005; pp. 745–775. [Google Scholar]
  34. Schmid, F. La famille des Xiphocentronides (Trichoptera: Annulipalpia). Mem. Entomol. Soc. Can. 1982, 121, 3–127. [Google Scholar] [CrossRef]
  35. StatSoft, Inc. Statistical User Guide; Complete Statistical System Statsoft: Oklahoma, OK, USA, 1995. [Google Scholar]
  36. Colwell, R.K. Estimates: Statistical Estimation of Species Richness and Shared Species from Samples (Version 9.1.0) User’s Guide and Application. 2013. Available online: http://viceroy.eeb.uconn.edu/estimates/EstimateSPages/EstSUsersGuide/EstimateSUsersGuide.htm (accessed on 1 September 2020).
  37. Jiménez-Valverde, A.; Hortal, J. Las curvas de acumulación de especies y la necesidad de evaluarla la calidad de los inventarios biológicos. Rev. Iber. Aracnol. 2003, 8, 151–161. [Google Scholar]
  38. González-Oreja, J.A.; De la Fuente-Díaz, A.A.; Hernández-Santín, L.D.; Buzo-Franco, D.; Bonache-Regidor, C. Evaluación de los estimadores no paramétricos de la riqueza de especies. Un ejemplo con las aves en áreas verdes de la ciudad de Puebla, México. Anim. Biodivers. Conserv. 2010, 33, 31–45. [Google Scholar] [CrossRef]
  39. Clarke, K.R.; Green, R.H. Statistical design and analysis for a ‘biological effects’ study. Mar. Ecol. Prog. Ser. 1988, 4, 213–226. [Google Scholar] [CrossRef]
  40. McAleece, N.; Gage, J.D.G.; Lambshead, P.J.D.; Paterson, G.L.J. BioDiversity Professional Statistics Analysis Software. Jointly Developed by the Scottish Association for Marine Science and the Natural History Museum London. 1997. Available online: https://www.sams.ac.uk/science/outputs/ (accessed on 26 November 2020).
  41. Jost, L. Entropy and Diversity. Oikos 2006, 113, 363–375. [Google Scholar] [CrossRef]
  42. Baselga, A.; Orme, C.D.L.; Villéger, S.; De Bortoli, J.; Leprieur, F. Betapart: Partitioning Beta Diversity into Turnover and Nestedness Components. R Package Version 1.4. 2017. Available online: http://CRAN.R-project.org/package=betapart (accessed on 26 August 2021).
  43. Chao, A.; Shen, T.J. Program SPADE (Species Prediction and Diversity Estimation). Program and User Guide. 2010. Available online: http://chao.stat.nthu.edu.tw/wordpress/software-download/ (accessed on 26 August 2021).
  44. Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 2010, 19, 134–143. [Google Scholar] [CrossRef]
  45. Zar, J.H. Biostatystical Analysis, 2nd ed.; Prentice-Hall: Upper Saddle River, NJ, USA, 1984; pp. 286–296. [Google Scholar]
  46. Bueno-Soria, J. Guía de Identificación Ilustrada de los Géneros de Larvas de Insectos del Orden Trichoptera de México, 1st ed.; Instituto de Biología-Universidad Nacional Autónoma de México: México City, Mexico, 2010; pp. 49–234. [Google Scholar]
  47. Razo-González, M. Caddisflies (Insecta: Trichoptera) from Santa Catarina Lachatao, Oaxaca, México: New species, new geographical records, and checklist, México. Zootaxa 2018, 438, 22–40. [Google Scholar] [CrossRef] [PubMed]
  48. Razo-González, M.; Castaño-Meneses, G.; Novelo-Gutiérrez, R.; Márquez, J. Preliminary evaluation of the nocturnal flight of caddisflies (Insecta: Trichoptera) in a temperate forest in Oaxaca, Mexico. Aquat. Insects 2020, 41, 339–355. [Google Scholar] [CrossRef]
  49. Razo-González, G.; Novelo-Gutiérrez, R.; Castaño-Meneses, G.; Márquez, J. Atopsyche Banks (Trichoptera: Hydrobiosidae) from Mexico: New species, redescription, and identification key. Stud. Neotrop. Fauna Environ. 2021. [Google Scholar] [CrossRef]
  50. Luis-Martínez, A.; Vargas-Fernández, I.; Llorente-Bousquets, J. Lepidopterofauna de Oaxaca I: Distribución y fenología de los Papilionoidea de la Sierra de Juárez; Publicaciones Especiales del Museo de Zoología: México City, Mexico, 1991; Volume 3, pp. 1–119. [Google Scholar]
  51. Rodríguez, A.; Castro-Castro, A.; Vargas-Amado, G.; Vargas-Ponce, O.; Zamora-Tavares, P.; González-Gallegos, J.; Carrillo-Reyes, P.; Anguiano-Constante, M.; Carrasco-Ortiz, M.; García-Martínez, M.; et al. Richness, geographic distribution patterns, and areas of endemism of selected angiosperm groups in Mexico. J. Syst. Evol. 2018, 56, 537–549. [Google Scholar] [CrossRef] [Green Version]
  52. Aragón-Parada, J.; Carrillo-Reyes, P.; Rodríguez, A.; Munguía-Lino, G. Diversidad y distribución del género Sedum (Crassulaceae) en la Sierra Madre del Sur, México. Rev. Mex. Biodivers. 2019, 90, 1–17. [Google Scholar] [CrossRef]
  53. Hernández-Rodríguez, E.; Aguirre-Hidalgo, V. Diversidad de musgos del bosque nublado de la Sierra de Juárez, Oaxaca, México. Acta Bot. Mex. 2020, 127, 1–14. [Google Scholar] [CrossRef] [Green Version]
  54. Flint, O.S.; Holzenthal, R.W.; Harris, S.C. Catalog of the Neotropical Caddisflies (Insecta: Trichoptera); Ohio Biological Survey: Columbus, OH, USA, 1999; pp. 12–14. [Google Scholar]
  55. Kimura, G.; Inoue, E.; Hirabayashi, K. Seasonal abundance of adult caddisfly (Trichoptera) in the middle reaches of the Shinano River in Central Japan. In Proceedings of the 6th International Conference on Urban Pests, Budapest, Hungary, 13–16 July 2008; National Education Centre-Press Ltd.: Veszprem, Hungary, 2008; pp. 259–266. [Google Scholar]
  56. Thompson, R.; Townsend, C. A truce with neutral theory: Local deterministic factors, species traits and dispersal limitation together determine patterns of diversity in stream invertebrates. J. Anim. Ecol. 2006, 75, 476–484. [Google Scholar] [CrossRef]
  57. Clarke, A.; Mac Nally, R.; Bond, N.; Lake, P.S. Macroinvertebrate diversity in headwater streams: A review. Freshw. Biol. 2008, 53, 1707–1721. [Google Scholar] [CrossRef]
  58. Monaghan, M.T.; Robinson, C.T.; Spaak, P.; Ward, J.V. Macroinvertebrate diversity in fragmented Alpine streams: Implications for freshwater conservation. Aquat. Sci. 2005, 67, 454–464. [Google Scholar] [CrossRef]
  59. Danehy, R.J.; Chan, S.S.; Lester, G.T.; Langshaw, R.B.; Turner, T.R. Periphyton and macroinvertebrate assemblage structure in headwaters bordered by mature, thinned, and clearcut douglas-fir stands. For. Sci. 2007, 53, 294–307. [Google Scholar] [CrossRef]
  60. Frady, C.; Johnson, S.; Li, J. Stream macroinvertebrate community responses as legacies of forest harvest at the H.J. Andrews experimental forest, Oregon. For. Sci. 2007, 53, 281–293. [Google Scholar] [CrossRef]
  61. Moldenke, A.R.; Linden, C.V. Effects of clearcutting and riparian buffers on the yield of adult aquatic macroinvertebrates from headwater streams. For. Sci. 2007, 53, 308–319. [Google Scholar] [CrossRef]
  62. Jacobsen, D. Altitudinal changes in diversity of macroinvertebrates from small streams in the Ecuadorian Andes. Arch. Hydrobiol. 2003, 158, 145–167. [Google Scholar] [CrossRef]
  63. Zacarías-Eslava, Y.; del Castillo, R.F. Comunidades vegetales templadas de la Sierra Juárez Oaxaca: Pisos altitudinales y sus posibles implicaciones ante el cambio climático. Bol. Soc. Bot. México 2010, 87, 13–28. [Google Scholar] [CrossRef] [Green Version]
  64. Cava, M.B.; Corronca, J.A.; Echeverría, A.J. Diversidad alfa y beta de los artrópodos en diferentes ambientes del Parque Nacional Los Cardones, Salta (Argentina). Rev. Biol. Trop. 2013, 61, 1785–1798. [Google Scholar] [CrossRef]
  65. Hepp, L.U.; Melo, A.S. Dissimilarity of stream insect assemblages: Effects of multiple scales and spatial distances. Hydrobiologia 2013, 703, 239–246. [Google Scholar] [CrossRef]
  66. Curry, C.J.; Baird, D.J. Habitat type and dispersal ability influence spatial structuring of larval Odonata and Trichoptera assemblages. Freshw. Biol. 2015, 60, 2142–2155. [Google Scholar] [CrossRef]
  67. Rohde, K. Latitudinal gradients in species diversity: The search for the primary cause. Oikos 1992, 65, 514–527. [Google Scholar] [CrossRef] [Green Version]
  68. Mora, A.; Szivák, I.; Deák, C.; Boda, R.; Csabai, Z.; Sály, P.; Takács, P.; Erős, T.; Bíró, P. Environmental factors influencing distribution of EPT assemblages in streams of Lake Balaton’s catchment area, Hungary. Zoosymposia 2011, 5, 360–371. [Google Scholar] [CrossRef]
  69. De Mendoza, G.; Ventura, M.; Catalan, J. Environmental factors prevail over dispersal constraints in determining the distribution and assembly of Trichoptera species in mountain lakes. Ecol. Evol. 2015, 5, 2518–2532. [Google Scholar] [CrossRef] [Green Version]
  70. Morrone, J.J. Regionalización biogeográfica y evolución biótica de México: Encrucijada de la biodiversidad del Nuevo Mundo. Rev. Mex. Biodiver. 2019, 9, 1–68. [Google Scholar] [CrossRef] [Green Version]
  71. Razo-González, M.; Márquez, J.; Castaño-Meneses, G.; Novelo-Gutiérrez, R. La complejidad biogeográfica de la Sierra de Juárez, Oaxaca, México, revelada a través del análisis de parsimonia de endemismos de especies de tricópteros (Insecta: Trichoptera). Rev. Mex. Biodivers. 2021, 92, e923808. [Google Scholar] [CrossRef]
  72. Halffter, G. La zona de transición mexicana y la megadiversidad de México: Del marco histórico a la riqueza actual. Dugesiana 2017, 24, 77–89. [Google Scholar] [CrossRef]
Diversity 15 00110 g001 550
Figure 1. Geographic location of the Sierra de Juárez, Oaxaca, Mexico, and five study sites. SS = Soyolapam stream; BS = Bobo stream; CC = Cascadas de Comaltepec; YS = Yegu-Iag-Lagu stream; FA = Flor de Azahar stream. Modified from Google Earth.
Figure 1. Geographic location of the Sierra de Juárez, Oaxaca, Mexico, and five study sites. SS = Soyolapam stream; BS = Bobo stream; CC = Cascadas de Comaltepec; YS = Yegu-Iag-Lagu stream; FA = Flor de Azahar stream. Modified from Google Earth.
Diversity 15 00110 g001
Diversity 15 00110 g002 550
Figure 2. Cluster analysis based on the Jaccard index and UPGMA representing the faunistic relationships between the rivers sampled from the Sierra de Juárez, Oaxaca, Mexico.
Figure 2. Cluster analysis based on the Jaccard index and UPGMA representing the faunistic relationships between the rivers sampled from the Sierra de Juárez, Oaxaca, Mexico.
Diversity 15 00110 g002
Diversity 15 00110 g003 550
Figure 3. Species accumulation curve of Trichoptera from the Sierra de Juárez, Oaxaca, Mexico.
Figure 3. Species accumulation curve of Trichoptera from the Sierra de Juárez, Oaxaca, Mexico.
Diversity 15 00110 g003
Diversity 15 00110 g004 550
Figure 4. True alpha diversity of Trichoptera in the Sierra de Juárez, Oaxaca, Mexico.
Figure 4. True alpha diversity of Trichoptera in the Sierra de Juárez, Oaxaca, Mexico.
Diversity 15 00110 g004
Diversity 15 00110 g005 550
Figure 5. Temporal variation in observed and expected true diversity of Trichoptera in the Sierra de Juárez, Oaxaca, Mexico.
Figure 5. Temporal variation in observed and expected true diversity of Trichoptera in the Sierra de Juárez, Oaxaca, Mexico.
Diversity 15 00110 g005
Diversity 15 00110 g006 550
Figure 6. Discriminant analysis based on the abundance of different taxa and different seasons in relation to the temperature and relative humidity in the Sierra de Juárez, Oaxaca, Mexico. Axes 1 and 2. (a) families “nortes”, (b) families rains, (c) genera “nortes”, (d) genera rains, (e) species “nortes”, (f) species rains.
Figure 6. Discriminant analysis based on the abundance of different taxa and different seasons in relation to the temperature and relative humidity in the Sierra de Juárez, Oaxaca, Mexico. Axes 1 and 2. (a) families “nortes”, (b) families rains, (c) genera “nortes”, (d) genera rains, (e) species “nortes”, (f) species rains.
Diversity 15 00110 g006
Table
Table 1. Study sites in the Sierra de Juárez, Oaxaca, Mexico. MAP = mean annual precipitation; MAT = mean annual temperature. Climate according to García notation [31].
Table 1. Study sites in the Sierra de Juárez, Oaxaca, Mexico. MAP = mean annual precipitation; MAT = mean annual temperature. Climate according to García notation [31].
StreamGeographic LocationAltitude (msnm)MAP (mm)MAT (°C)Vegetation TypeClimate
Soyolapam (SS)17°41′47″ N
96°16′54″ W		1363000–350024.9Evergreen rainforestAw
Bobo (BS)17°39′23″ N
96°23′43″ W		4602500–300030Evergreen rainforestAf
Yegu-Iag-Lagu (YS)17°24′18″ N
96°23′11″ W		13201200–230021Cloud forest(A)c(m)
Flor de Azahar (FA)17°13′21″ N
96°19′39″ W		1967800–120022.5Oak-pine forestCw
Cascadas Comaltepec (CC)17°35′15″ N
96°29′35″ W		24322000–25009Pine forestCw
Table
Table 2. Mean, standard error (±SE), minimum and maximum values of air temperature (T °C) and relative humidity (RH) of the five study sites in the Sierra de Juárez, Oaxaca, Mexico. n = “nortes”, r = rainy. N = 50.
Table 2. Mean, standard error (±SE), minimum and maximum values of air temperature (T °C) and relative humidity (RH) of the five study sites in the Sierra de Juárez, Oaxaca, Mexico. n = “nortes”, r = rainy. N = 50.
Stream SSBSYSFACC
Seasonnrnrnrnrnr
T °C
(±SE) 		16.74 (0.13)21.90 (0.24)15.73 (0.04)19.71 (0.14)10.31 (0.07)19.42 (0.33)9.43 (0.26)14.89 (0.37)2.26 (0.13)11.77 (0.11)
Minimum15.6120.2215.4118.729.4617.086.0812.941.4010.88
Maximum18.1829.2416.3923.6211.6624.7714.1920.893.4914.65
RH%
(±SE) 		90.87 (0.15)88.60 (0.85)91.62 (0.10)91.49 (0.49)91.35 (0.09)92.13 (0.26)80.71 (0.95)93.19 (0.36)84.84 (0.87)94.45 (0.53)
Minimum88.3368.7489.8780.8789.6587.2169.2581.8972.3287.33
Maximum92.5392.7492.7293.7592.2694.3091.4796.7390.3998.00
Table
Table 3. Richness and composition of adult Trichoptera assemblages in the Sierra de Juárez, Oaxaca, Mexico. Taxa are alphabetically arranged. x = “nortes”, * = rainy, + = new species for the science; in bold the new record for state of Oaxaca; Freq = capture frequency (percentage of catch of each species in the ten samplings).
Table 3. Richness and composition of adult Trichoptera assemblages in the Sierra de Juárez, Oaxaca, Mexico. Taxa are alphabetically arranged. x = “nortes”, * = rainy, + = new species for the science; in bold the new record for state of Oaxaca; Freq = capture frequency (percentage of catch of each species in the ten samplings).
Species/AssemblageSSBSYSFACCFreq
Calamoceratidae
  Phylloicus aeneus (Hagen, 1861)** * 30
Glossosomatidae
  Culoptila jamapa Bueno-Soria and Santiago-Fragoso, 1996 x *x * 40
  Culoptila pararusia Blahnik and Holzenthal, 2006x * 20
  Mortoniella falcicula Blahnik and Holzenthal, 2008 * 10
  Mortoniella meralda (Mosely, 1954)x * 20
   Mortoniella mexicana Blahnik and Holzenthal, 2008*x 20
  Protoptila bicornuta Flint, 1963* 10
  Protoptila ixtala Mosely, 1937x * 20
   Protoptila liqua Mosely, 1954x * 20
  Protoptila resolda Mosely, 1937* 10
  Protoptila rota Mosely, 1937* 10
   Protoptila spangleri Flint, 1967x * 20
Helicopsychidae
  Cochliopsyche vazquezae Flint, 1986* 10
  Helicopsyche borealis (Hagen, 1861)x *** 40
   Helicopsyche dampfi Ross, 1956 * 10
  Helicopsyche planata Ross, 1956xx* 30
  Helicopsyche tuxtlensis Bueno-Soria, 1983x** 30
   Helicopsyche sp.+ * 10
Hydrobiosidae
   Atopsyche bifurcata Razo-González and Novelo-Gutiérrez, 2021 *10
  Atopsyche calopta Ross and King, 1952 x x 20
  Atopsyche dampfi Ross and King, 1952xxx *x * 60
  Atopsyche erigia Ross, 1947* 10
   Atopsyche huenga Flint, 1974 * 10
   Atopsyche japoda Ross and King, 1952 xx 20
   Atopsyche jujmi Razo-González and Novelo-Gutiérrez, 2021 *10
   Atopsyche majada Ross, 1947x ** 30
  Atopsyche pilcomayo Schmid, 1989 *10
Hydropsychidae
  Cheumatopsyche gelita Denning, 1952x 10
  Diplectrona chiapensis Flint, 1967 xxx **50
  Hydropsyche toschiae (Denning, 1965)xx *x *x * 70
   Hydropsyche sp. 1+ x *x **50
   Hydropsyche sp. 2+x ** 30
  Leptonema albovirens (Walker, 1852)x * 20
  Leptonema championi Mosely, 1933 x 10
  Leptonema crassum Ulmer, 1905* 10
  Leptonema pinotepa Bueno-Soria et al., 2001 x * 20
  Leptonema simulans mayanum Flint et al., 1987 x * 20
  Plectropsyche hoogstraali Ross, 1947x *x * 40
   Smicridea bulara Flint and Denning, 1989* 10
   Smicridea dampfi Flint, 1974 * 10
   Smicridea dithyra Flint, 1974x *x ** 50
   Smicridea lachanha Bueno-Soria and Hamilton, 1986x * 20
  Smicridea radula Flint, 1974x *x ** 50
  Smicridea signata (Banks, 1903)* 10
Hydroptilidae
   Anchitrichia spangleri Flint, 1970x * 20
   Costatrichia lodora Mosely, 1937x * 20
  Hydroptila denza Ross, 1948* 10
  Hydroptila mexicana Mosely, 1937x * 20
  Hydroptila misolha Bueno-Soria, 1984* 10
  Ithytrichia mexicana Harris and Contreras-Ramos, 1989 ** 20
  Leucotrichia imitator Flint, 1970 x 10
   Mayatrichia rualda Mosely, 1937 x ** 30
   Metrichia cf. circulatrix Bueno-Soria 2002 * 10
   Metrichia crenula Bueno-Soria 2002 x 10
  Metrichia yavesia Bueno-Soria, 2002 x * 20
   Metrichia sp. 1+ * 10
   Metrichia sp. 2+* 10
   Neotrichia exicoma (Mosely, 1937)* 10
   Neotrichia xicana (Mosely, 1937)*x 20
  Neotrichia cf. maria Bueno and Hamilton, 1986 x 10
   Ochrotrichia pacifica Flint, 1972x * 20
  Ochrotrichia yavesia Bueno and Holzenthal, 2004 *x 20
  Ochrotrichia yetla Bueno-Soria, 2009*x * 30
   Ochrotrichia sp. 1+ * 10
   Ochrotrichia sp. 2+ ** 20
  Oxyethira azteca (Mosely, 1937)xx 20
   Oxyethira tica Holzenthal and Harris, 1992* 10
   Rhyacopsyche mexicana (Flint, 1967)* 10
   Rhyacopsyche sp.+ * 10
   Zumatrichia filosa Mosely, 1937* 10
   Zumatrichia multisetosa Flint, 1970* 10
   Zumatrichia sp.+*x 20
Lepidostomatidae
  Lepidostoma bakeri Flint, 1965 x * 20
  Lepidostoma frontale (Banks, 1901) xx * 30
  Lepidostoma zapoteca Razo-González, 2018 x 10
   Lepidostoma sp. 1+ x *x 30
   Lepidostoma sp. 2+ x *20
   Lepidostoma sp. 3+ ** 20
Leptoceridae
   Nectopsyche dorsalis (Banks, 1901)x * 20
  Nectopsyche gemmoides Flint, 1981x *x ** 50
   Nectopsyche pavida (Hagen, 1861)x 10
   Oecetis disjuncta (Banks, 1920)x *x * 40
   Oecetis metlacensis Bueno-Soria, 1981x *x *** 60
   Oecetis silviae Bueno-Soria, 1981 ** 20
  Oecetis sordida Blahnik and Holzenthal, 2014 x * 20
   Oecetis sp.+* 10
   Triaenodes sp.+ *10
Odontoceridae
  Marilia flexuosa Ulmer, 1905* 10
Philopotamidae
  Chimarra acuta Ross, 1959* 10
  Chimarra barrettae (Banks, 1900)x ** 30
  Chimarra bicolor (Banks, 1901)*** 30
  Chimarra dentosa Ross, 1948* 10
  Chimarra laguna Ross, 1951* 10
  Chimarra mexicana (Banks, 1900) xx ** 40
  Chimarra oaxaca Blahnik, 1998 * 10
  Chimarra ortiziana Flint, 1967* 10
  Chimarra pelaezi Bueno-Soria, 1985*x ** 40
  Chimarra ridleyi (Denning, 1941)* 10
  Chimarra setosa Ross, 1959* 10
  Wormaldia barbai Muñoz-Quesada and Holzenthal, 2015*x *** 50
  Wormaldia noveloi Razo-González, 2018 * *20
  Wormaldia planae Ross and king, 1956x 10
  Wormaldia tarasca Bueno-Soria and Holzenthal, 1986 * 10
Polycentropodidae
   Cernotina calcea Ross, 1938* 10
   Cernotina chiapaneca Bueno-Soria, 2010* 10
  Polycentropus aliciae Barba-Álvarez and Bueno-Soria, 2005 * 10
  Polycentropus aztecus Flint, 1967 x * 20
  Polycentropus casicus Denning and Sykora, 1966 * 10
   Polycentropus mayanus Flint, 1981x ** 30
  Polycentropus mexicanus (Banks, 1901) x * 20
  Polycentropus veracruzensis Flint, 1981 ** 20
   Polyplectropus carolae Bueno-Soria, 1990 * 10
  Polyplectropus charlesi (Ross, 1941)* 10
Psychomyiidae
  Tinodes sp.x * 20
Xiphocentronidae
   Cnodocentron sp.+* 10
   Xiphocentron aureum Flint, 1967 * 10
Table
Table 4. Observed, estimated, and exclusive richness, and observed and estimated true diversity of Trichoptera in an altitudinal gradient from the Sierra de Juárez, Oaxaca, Mexico. In parentheses the percentage of richness and diversity obtained in the field with respect to those estimated.
Table 4. Observed, estimated, and exclusive richness, and observed and estimated true diversity of Trichoptera in an altitudinal gradient from the Sierra de Juárez, Oaxaca, Mexico. In parentheses the percentage of richness and diversity obtained in the field with respect to those estimated.
RichnessTrue Diversity
Site♂ AbundanceObs.ClenchChao2Excl.Obs.Est.
SSn5043447 (72.34)39 (87.18)13(11.02)6.606.90 (95.65)
BSn1912740 (67.50)69 (39.13)4(3.39)10.7612.08 (89.07)
YSn251224 (50.00)13 (92.31)1(0.85)10.8014.57 (74.12)
Fan291635 (45.71)30 (53.33)6(5.08)13.6523.13 (59.01)
CCn11--1(0.85)1-
SSr50386167 (91.04)66 (92.42)41(34.75)11.6811.77 (99.24)
BSr11383238 (94.12)62 (51.61)6(5.08)4.864.99 (97.39)
YSr5273755 (67.27)55 (67.27)8(6.78)10.4110.94 (95.16)
Far2832530 (83.33)29 (86.21)10(8.47)10.8511.43 (94.93)
CCr44811 (72.73)8 (100)5(4.24)5.055.50 (91.82)
SJ7402118133 (88.72)134 (88.06)77(65.25)17.2219.45 (88.53)
Table
Table 5. Recorded species richness (S), parameters a (rate of increase of new species), b (slope of the curve), and predictions of Clench’s species accumulation model by season in Trichoptera assemblages from the Sierra de Juárez, Oaxaca, Mexico.
Table 5. Recorded species richness (S), parameters a (rate of increase of new species), b (slope of the curve), and predictions of Clench’s species accumulation model by season in Trichoptera assemblages from the Sierra de Juárez, Oaxaca, Mexico.
Clench
StreamSnaba/bR2Slope%
SSn34323.760.0847.240.990.3071.98
BSn27322.420.0640.050.990.2867.42
YSn12320.7600.0324.250.990.1949.49
FAn16320.890.0235.950.990.2844.50
CCn116------
SSr613216.540.2567.120.990.2190.89
BSr32324.790.1238.930.990.2082.21
YSr37323.270.0655.030.990.3967.23
FAr25323.720.1230.860.990.1681.01
CCr8320.980.0911.190.990.0771.50
SJtotal1183202.440.02133.410.990.9787.70
Table
Table 6. Global beta diversity values (βcc), species turnover (β-3), and differences in richness (βrich) among the Trichoptera assemblages from the Sierra de Juárez, Oaxaca, Mexico.
Table 6. Global beta diversity values (βcc), species turnover (β-3), and differences in richness (βrich) among the Trichoptera assemblages from the Sierra de Juárez, Oaxaca, Mexico.
BSYSFACCMean
SSβcc0.720.850.921.000.87
β-30.590.780.861.000.81
βrich0.130.070.070.000.07
BSβcc 0.650.870.980.83
β-3 0.630.830.930.80
βrich 0.010.040.050.03
YSβcc 0.720.930.83
β-3 0.630.770.70
βrich 0.090.170.13
FAβcc 0.940.94
β-3 0.860.86
βrich 0.080.08
Table
Table 7. Values of total beta diversity between Trichoptera assemblages and between “nortes” and rainy seasons in the Sierra de Juárez, Oaxaca, Mexico.
Table 7. Values of total beta diversity between Trichoptera assemblages and between “nortes” and rainy seasons in the Sierra de Juárez, Oaxaca, Mexico.
SSnBSnYSnFAnCCnSSrBSrYSrFArCCr
SSn0.000.930.990.991.000.720.930.980.991.00
BSn 0.000.930.891.000.910.930.780.870.99
YSn 0.000.321.001.000.990.670.740.99
FAn 0.001.001.000.980.610.520.95
CCn 0.001.001.001.001.001.00
SSr 0.000.410.880.991.00
BSr 0.000.900.971.00
YSr 0.000.220.99
FAr 0.000.99
CCr 0.00
Table
Table 8. Values of F for two-way ANOVA test to evaluate season, site, and interaction (season X site) effects on total richness, total abundance of Trichoptera, and on the two dominant species, in the Sierra de Juárez, Oaxaca, Mexico. * p < 0.001; gl season = 1, gl site = 4, gl interaction = 4.
Table 8. Values of F for two-way ANOVA test to evaluate season, site, and interaction (season X site) effects on total richness, total abundance of Trichoptera, and on the two dominant species, in the Sierra de Juárez, Oaxaca, Mexico. * p < 0.001; gl season = 1, gl site = 4, gl interaction = 4.
Total RichnessTotal AbundanceNectopsyche gemmoidesSmicridea radula
Season256.40 *147.65 *37.39 *70.74 *
Site127.55 *74.55 *132.08 *38.84 *
Interaction16.79 *24.11 *33.26 *27.97 *
Table
Table 9. Multiple correlation coefficients between season, site, total richness, total abundance of Trichoptera, and the two dominant species, in the Sierra de Juárez, Oaxaca, Mexico. gl = degrees free, * significative differences.
Table 9. Multiple correlation coefficients between season, site, total richness, total abundance of Trichoptera, and the two dominant species, in the Sierra de Juárez, Oaxaca, Mexico. gl = degrees free, * significative differences.
Total RichnessTotal AbundanceNectopsyche gemmoidesSmicridea radula
Seasonr2 = 0.47 *
gl = 2.32
p = 0.00		r2 = 0.42 *
gl = 2.32
p = 0.00		r2 = 0.19 *
gl = 2.32
p = 0.00		r2 = 0.33 *
gl = 2.32
p = 0.00
Siter2 = −0.43 *
gl = 2.32
p = 0.00		r2 = −0.43 *
gl = 2.32
p = 0.00		r2 = −0.52 *
gl = 2.32
p = 0.00		r2 = −0.44 *
gl = 2.32
p = 0.00
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Razo-González, M.; Novelo-Gutiérrez, R.; Castaño-Meneses, G.; Márquez, J. Diversity and Composition of Caddisflies (Insecta: Trichoptera) along an Elevation Gradient in Southeastern Mexico. Diversity 2023, 15, 110. https://doi.org/10.3390/d15010110

AMA Style

Razo-González M, Novelo-Gutiérrez R, Castaño-Meneses G, Márquez J. Diversity and Composition of Caddisflies (Insecta: Trichoptera) along an Elevation Gradient in Southeastern Mexico. Diversity. 2023; 15(1):110. https://doi.org/10.3390/d15010110

Chicago/Turabian Style

Razo-González, María, Rodolfo Novelo-Gutiérrez, Gabriela Castaño-Meneses, and Juan Márquez. 2023. "Diversity and Composition of Caddisflies (Insecta: Trichoptera) along an Elevation Gradient in Southeastern Mexico" Diversity 15, no. 1: 110. https://doi.org/10.3390/d15010110

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop