Unraveling the importance of functionally extreme tadpole types to functional diversity: a case study in temperate montane streams
Frontiers in Zoology volume 20, Article number: 7 (2023)
Functional diversity is important to maintain ecosystem functioning. Species with different ecomorphological traits may display distinct functional roles in ecosystems. Accordingly, functionally extreme species are more important as they can exhibit specific strategies. However, little is known about the distribution patterns of functionally extreme species at a local scale and whether the prior extinction of extreme species can cause significant effects on functional diversity. In addition, no empirical studies have been conducted on the microhabitat determinants of extreme species to maintain the functional diversity.
This study collected 1470 tadpoles belonging to 6 families and 20 anuran species. These species were subsequently divided into 65 functional entities based on their developmental stages to incorporate intraspecific traits variability. As a result, we detected seven extreme functional entities, accounting for 10.7% of the total number of entities. Moreover, the prior extinction of extreme entities can lead to a significant decrease in functional diversity compared with the random extinction of entities. Microhabitat variables such as conductivity, water depth, and current velocity determined the distribution of extreme entities.
Although the functionally extreme entities only represented a small proportion of the total number of tadpoles, they played irreplaceable roles in maintaining functional diversity. Their extinction may induce high functional vulnerability in tadpole communities. Therefore, anuran species with extreme tadpole traits need to be projected for amphibian conservation.
Biodiversity has a positive relationship with ecosystem functioning and has attracted the attention of ecologists for a long time [1, 2]. However, most previous studies only focused on the taxonomic facet of biodiversity (e.g., species richness), despite biodiversity contains multidimensional information on functional and phylogenetic facets [3, 4]. In recent decades, with the development of functional ecology, increasingly studies have started to quantify the functional dimensions of biodiversity (i.e., functional diversity) using the functional characteristics of species in a particular region [5, 6]. As a result, functional diversity is considered to be more sensitive and predictive in quantifying the relationship between biodiversity and ecosystem functioning . For instance, the increased fish functional diversity can promote material circulation and energy flow in a tropical stream in Água Azeda Creek, Bahia, Brazil . On the other hand, for birds in some northern hemisphere countries such as China, Russia, and Finland, the decrease in their functional diversity strongly threatens the ecosystem functioning of forests due to human activity and seed dispersal .
Among all the organisms, amphibians have rich diversity and contribute to many key ecosystem processes, such as promoting material exchange, regulating nutrient cycling and energy flow, and controlling pests . This is especially true for their larval tadpoles, which play critical functional roles in freshwater ecosystems, especially in montane streams . For instance, tadpole excretion and digestion behaviors can promote nutrient circulation and litter decomposition in streams . Therefore, their consumption activities and interactions can strongly affect stream food web structures . Previous tadpole functional studies usually selected some ecomorphological traits as their functional traits [14, 15]. Although these traits do not reflect the fundamental functions of tadpoles in ecosystems, they can also reflect critical ecological functions through food acquisition and locomotion to regulate tadpole contributions to ecosystem processes . Accordingly, studies have been conducted to investigate the responses of tadpole functional traits to perturbations  and to explore the spatial patterns of tadpole functional diversity [18,19,20]. However, little attention has been paid to the distinct ecological strategies of different tadpoles.
The distribution patterns of traits of organisms reflect their distinct ecological strategies . As the trait distribution is often unimodal in animals , species with the core distribution are generalists, while species distributed in tails can be considered as the specialists . Although only a few species had extreme traits in communities , the importance of functionally extreme species to functional diversity is still unclear. Moreover, it is widely recognized that species traits are strongly affected by environmental conditions (e.g., plants: , fish: , and amphibian tadpoles: ). However, empirical evidence related to the environmental determinants of the distribution of extreme species is still limited.
In the present study, we used anuran tadpoles in montane streams of a temperate forest as models to explore the importance of functionally extreme species to functional diversity. Specifically, we aim to (1) identify the extreme functional tadpole types in the studied montane streams, (2) investigate whether the prior extinction of extreme tadpole types can cause significant effects on functional diversity indices than the random extinction of tadpoles, and (3) quantify the microhabitat determinants of the distribution of extreme tadpole types (Fig. 1).
Materials and methods
This study was conducted in the montane streams of Mount Emei (29°16′–29°43′ N, 103°11′–103°37′ E) located on the southwestern edge of the Sichuan Basin, China. The highest altitude of this mountain is 3099 m, and the altitude span is about 2600 m . This area belongs to a typical subtropical climate. However, a large elevational gradient induced four distinct climate zones from low to high elevations. Specifically, it is the subtropical climate below 1500 m, which was covered by evergreen broad-leaved forest. The warm temperate climate is between 1500 and 2100 m, and the main vegetation cover is evergreen, deciduous broad-leaved mixed forest. Middle temperate climate can be detected between 2100 and 2500 m, with broad-leaved and coniferous forests dominating this area. In addition, it is the sub-cold zone climate with cold temperate coniferous forest above 2500 m . Abundant vegetation and freshwater resources provide various microhabitats for amphibians, and 35 amphibian species were recorded in this area based on literatures .
Field-work was conducted in 18 streams, covering an elevational gradient from 485 to 2865 m (Additional file 1: Table S1). These streams contained diverse microhabitats, providing suitable habitats for different tadpoles. Only 300 m of each stream was selected and sampled, as high waterfalls or cliffs can usually block the streams. Specifically, tadpoles in each stream were sampled between 20:30 h and 23:00 h at night using dip nets (opening diameter: 40 cm; depth: 35 cm; mesh size: 3 mm; RBF326, Renniaofei, China) from May to June in 2018 and 2019, respectively, with one stream being sampled per night. The developmental stages of sampled tadpoles were determined according to Gosner . For individuals that can be directly identified to species by morphological characteristics (e.g., Leptobrachium boringii and Quasipaa boulengeri), we collected at least 15 individuals of each species at the different developmental stages and brought them back to the laboratory for further analyses. The rest individuals were counted and released back to their original habitats. For individuals that cannot be identified to species directly by external morphology (e.g., Megophrys omeimontis and Megophrys shapingehsis), all the sampled individuals were brought back to the laboratory and identified to species using DNA barcoding analyses based on mitochondrial 16S rRNA gene fragments .
We also recorded 12 microhabitat variables based on previous studies showing that they can potentially determine tadpole distribution and diversity [5, 31, 32]. These variables included conductivity (Con), water depth (Wd), river width (Rw), dissolved oxygen (Do), water pH (pH), water temperature (Wt), total nitrogen (TN), total phosphorus (TP), current velocity (Cv), ammonia nitrogen (NH4-N), chlorophyll a (Chl. a), and substrate type (Sub). In addition, a previous study has provided details of the measurement approaches .
Functional traits and functional entities
Based on previous studies [15, 33], we selected 11 functional traits related to the two key ecological functions (i.e., food acquisition and locomotion) of tadpoles displayed in aquatic ecosystems. Specifically, there are four functional traits related to food acquisition (i.e., oral disk shape, oral disk position, eye position, and gut length), and six functional traits are related to locomotion (i.e., tail shape, tail position, tail throttling, trunk bending shape, spiracle position, and body section shape). In addition, mass is related to both food acquisition and locomotion. These traits were calculated as the ratios of 11 morphological traits (Additional file 1: Table S2). Details of the measurement of these traits can be found in a previous report .
In the past two decades, increasing studies have reported the importance of intraspecific traits variability in functional ecological studies [34,35,36], and they suggested that species should be divided into different functional entities according to their life stages (e.g., [6, 37]). Therefore, individuals within tadpole of each anuran species were divided into five functional entities based on their developmental stages following the description in a previous study . These stages included functional entity 1 (Gosner stages 1–25; embryonic stage), functional entity 2 (Gosner stages 26–30; developing period of limb bud), functional entity 3 (Gosner stages 31–39; developing period of toes), functional entity 4 (Gosner stages 40–41; preparation period of metamorphosis), and functional entity 5 (Gosner stages 42–46; metamorphosis period ).
First, we conducted a principal component analysis (PCA) based on scaled 11 functional traits of all tadpole functional entities. The first four principal components had their eigenvalues > 1, which were selected to construct the functional traits space. They occupied 76.74% of the initial inertia in trait values (PC1 = 27.25%, PC2 = 23.29%, PC3 = 15.74%, and PC4 = 10.46%, respectively). Following Su et al. , the 0.5% of entities with the lowest values and the 0.5% of entities with the highest values on each of the four PCA axes were defined as the extreme functional entities (EFE), and the remaining were defined as the general functional entities (GFE).
All the entities were gathered together as the community (i.e., the regional pool), as amphibians could move freely in the mountain. The functional diversity of the community was assessed based on three indices, including functional richness (FRic), functional originality (FOri), and functional specialization (FSpe). We selected these indices because their calculations involved in extreme functional entities . Specifically, functional richness indicated the volume of functional space occupied by all entities . Functional originality referred to a metric of original isolation when entities occupied functional space in a community, quantifying how altered functional redundancy between entities [39, 40]. Functional specialization reflected the specialization of the entities in the functional space, showing how the abundance of extreme entities increased [24, 39, 41, 42]. We generated null models to explore whether the prior extinction of EFE/GFE can cause significant effects on functional diversity indices than the random extinction of entities from the regional pool. Finally, we designed an appropriate randomization procedure to test the changes in functional diversity indices with the extinct process of EFE/GFE (see Additional file 1: Appendix Figure S1 for details). The null models randomized the tadpole assemblages after EFE/GFE had been extinct. The randomization process was carried out 999 times for each index, and the P values were calculated to indicate the statistical significance. In addition, standardized effect sizes (SES) were calculated as the differences between the observed functional diversity values and the random functional diversity values standardized by the standard deviation.
Finally, we used generalized linear models (GLMs) to explore the microhabitat determinants of the richness and the relative abundance of EFE, respectively. Specifically, we used a Poisson distribution with a logarithmic link function in the models, and different models were compared by deleting variables one by one from the global model. Because of the small sample size, the best-fitted models were selected based on the AICc values . Moreover, model averaging was performed if needed . Furthermore, the significant variables were identified when P < 0.05. After that, hierarchical partitioning analyses were conducted to quantify the relative contribution of each microhabitat variable in the best model to the richness and to the relative abundance of EFE, respectively.
All calculations and statistical analyses were performed in R 4.0.3 . PCA analyses were conducted using vegan package . GLMs were performed using the lme4 package . Hierarchical partitioning was performed using hier.part package . The null models were performed by picante package .
Functional entities with extreme traits
Within two years of sampling, we collected 1496 tadpole individuals, which belonged to 20 anuran species from 6 families (Table 1). These individuals were subsequently divided into 65 functional entities based on their developmental stages, with each functional entity being measured for 23.0 ± 47.2 (mean ± SD) individuals. Specifically, the common entities were tadpole of Leptobrachium boringii 1, Oreolalax omeimontis 1, and Megophrys omeimontis 3, accounting for 20.7%, 12.4%, and 9.5% of the total number of individuals, respectively. The rare entities were tadpole of Oreolalax major 1, Leptobrachella oshanensis 4, and Fejervarya multistriata 5, accounting for 0.06%, 0.12%, and 0.17% of the total number of individuals, respectively.
In the functional space, PC1 was mainly contributed by functional traits related to locomotion, with higher values indicating stronger mobility. PC2 was mainly contributed by functional traits related to food acquisition, and the higher values in PC2 demonstrated that tadpoles preferred to prey on the surface of the water. Moreover, these tadpoles also exhibited strong detection but weak digestion ability. In addition, PC3 and PC4 were mainly driven by spiracle position and mass, respectively. Based on the distribution in the four PC axes, seven functional entities were identified as EFE, accounting for 10.7% of the total number of functional entities. Specifically, these entities were tadpole of L. boringii 5, Amolops chunganensis 4, F. multistriata 5, Megophrys minor 5, O. omeimontis 5, L. oshanensis 2, and Rana omeimontis 1 (Fig. 2). Furthermore, they filled 9.75%, 3.29%, 14.25%, and 19.38% of the range of PC1, PC2, PC3, and PC4, respectively (Additional file 1: Figure S2).
Effects of tadpole extinction on functional diversity
Overall, the prior extinction of EFE had significant effects on functional diversity compared to those calculated from the random extinction of functional entities (Fig. 3). Specifically, functional richness could significantly decrease after prior extinction of one EFE (SES = − 2.087, P = 0.028). However, there was no significant effects on functional specialization (SES = − 1.624, P = 0.81) and functional originality (SES = − 1.497, P = 0.058). All the functional diversity indices for the prior extinction of two to seven EFE were significantly lower than those calculated based on random extinction (Additional file 1: Table S4; Fig. 3). In contrast, the prior extinction of one to seven GFE had no significant effects on all the functional diversity indices (Additional file 1: Table S4; Fig. 3).
Environmental determinants of the distribution of EFE
In the best model constructed by GLMs, six microhabitat factors (i.e., elevation, conductivity, water depth, current velocity, total phosphorus, and substrate type) significantly impacted the distribution of EFE richness (Table 2 and Additional file 1: Table S3). Specifically, conductivity and substrate type were significantly and positively correlated with EFE richness (P < 0.05), while water depth, elevation, current velocity, and total phosphorus had negative and significant correlations with EFE richness (P < 0.05). The results of hierarchical partitioning showed that total phosphorus (25.5%), current velocity (25.3%), and conductivity (24.4%) were the three most important factors contributing to the variance of EFE richness. Moreover, water depth and substrate type contributed 11.6% and 10.6% to the variance of EFE richness, respectively. In addition, the independent contribution of elevation was only 2.5% (Fig. 4A).
In terms of the distribution of the relative biomass of EFE, there were four significant microhabitat factors (i.e., conductivity, water depth, current velocity, and total phosphorus; Table 3 and Additional file 1: Table S3) in the best model constructed by GLMs. Specifically, conductivity had a significant and positive correlation with the relative biomass of EFE (P < 0.05), while water depth, total phosphorus, and current velocity were significantly and negatively correlated with it (P < 0.05). Hierarchical partitioning analyses indicated that the independent contribution of conductivity to the variance of EFE relative biomass was 44.0%, followed by water depth (20.7%), current velocity (19.7%), and total phosphorus (15.6%; Fig. 4B).
The present study revealed that seven of the 65 anuran tadpole functional entities we sampled were identified as EFE. Although these EFEs accounted for a small proportion of the total number of tadpole entities, they had an important contribution to each PC axis. Interestingly, these functional entities were either in the early developmental stage (e.g., R. omeimontis 1 and L. oshanensis 2) or in the metamorphosis period (e.g., L. boringii 5, F. multistriata 5, M. minor 5, and O. omeimontis 5), exhibiting distinct eco-morphology from other entities. Specifically, O. omeimontis 5, M. minor 5, and A. chunganensis 4 were distributed in the functional space with positive values for PC1, but negative values for PC2. Our results were compensated with previous findings showing that the locomotion of these tadpoles was high acceleration and/or maneuverability and low magnitude of vertebral curvature [41, 50]. This is probably because they were close to completing the metamorphosis and thus exhibited more vital locomotion ability. These tadpoles also had higher gut length, lower oral disk shape, and eye position values, indicating that the digestion ability was poor, and their prey detection was more focused on the water surface. This is because these tadpoles were in the metamorphosis period, and their functional traits exhibited significant changes to adapt living from water to land . Moreover, F. multistriata 5, R. omeimontis 1, and L. boringii 5 were distributed in the functional space with negative values for PC1 and PC2. Accordingly, the locomotion of these tadpoles was low acceleration and/or high magnitude of vertebral curvature [15, 50, 52]. This may be because these tadpoles either had poor locomotion ability (i.e., R. omeimontis 1) or relied more on the limbs to move (i.e., F. multistriata 5 and L. boringii 5). In addition, L. oshanensis 2 was distributed in the functional space with negative values for PC1, but positive values for PC2. Therefore, these tadpoles relied more on the vertebral column to swim, and they can detect prey easily on both sides of the body. Consequently, they can swim and prey for food agilely .
The results of the null models indicated that the prior extinction of EFE can cause a significant decrease in tadpole functional richness. These results showed that despite EFE was rare in a community, they occupied most of the functional space. Therefore, our results were consistent with previous studies showing that the functional richness of a community was mainly supported by EFE [39, 42]. Moreover, the extinction of EFE in the community can cause some vacant functional space, leading ecological resources not to be effectively utilized in the ecosystems . Indeed, montane streams were fragile ecosystems, and their energy input mainly relied on leaf litter decomposition . The vacancy of ecological space caused by the extinction of EFE would break the ecological balance and thus decrease the stability of ecosystems. More importantly, the prior extinction of tadpole EFE can also induce a significant decrease in functional specialization and functional originality, causing an increase in functional redundancy and functional vulnerability in the community . Consequently, these communities can become less resistant and resilient when facing human disturbances [39, 56]. The increase in functional redundancy is also related to the weak and asymptotic relationship between species richness and ecosystem functioning because of the increased functional homogenization [56, 57]. Furthermore, high homogenization can reduce spatial diversity and potentially cut down the ecosystem functioning, thus jeopardizing ecosystem services .
The results of GLMs showed that EFE tended to live in water conditions with high conductivity, shallow water, low current velocity, and low phosphorus. These environments are usually related to the high abundance of food resources [52, 59, 60], thus can provide suitable habitats for their growth (e.g., R. omeimontis 1 and L. oshanensis 2). Previous studies indicated that the adults of R. omeimontis and L. oshanensis preferred to spawn in still and clean water bodies . Therefore, the selection of adult breeding sites may also determine the distribution of tadpole EFE. Moreover, these water conditions can also benefit the growth of aquatic plants, associated with high habitat heterogeneity . Other EFE, such as F. multistriata 5, L. boringii 5, M. minor 5, and O. omeimontis 5 prefer to live in these places as they can find abundant shelters during metamorphosis .
Overall, this study identified the anuran species with extreme tadpole traits in streams of Mount Emei. These species exhibited specific ecomorphological traits and occupied most of the functional space. More importantly, although EFE accounted for only a small proportion (10.7%) of the total number of entities, they played irreplaceable roles in maintaining the functional diversity of the community. Accordingly, EFE should be better concerned and protected to preserve tadpole diversity. Furthermore, since EFE reflects the functional vulnerability of ecosystems, future research can explore how the extinction of EFE will affect ecosystem functioning. In addition, since the abundance of EFE was not considered in the present study, functional diversity indices incorporating entity traits and abundance could be applied to better understand the effects of EFE on biodiversity. In addition, the EFE distribution was strongly determined by microhabitat variables such as conductivity, water depth, and current velocity. Therefore, microhabitat protection should be the priority to maintain the richness and abundance of EFE.
Availability of data and materials
Data are available from the authors upon reasonable request.
Hector A, Joshi J, Lawler S, Spehn EM, Wilby A. Conservation implications of the link between biodiversity and ecosystem functioning. Oecologia. 2001;129:624–8.
Kardol P, Fanin N, Wardle DA. Long-term effects of species loss on community properties across contrasting ecosystems. Nature. 2018;557:710–13.
Devictor V, Mouillot D, Meynard C, Jiguet F, Thuiller W, Mouquet N. Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world. Ecol Lett. 2010;13:1030–40.
Lin L, Deng WD, Huang XX, Kang B. Fish taxonomic, functional, and phylogenetic diversity and their vulnerabilities in the largest river in southeastern China. Ecol Evol. 2021;11:11533–48.
Sun Z, Zhao C, Zhu W, Zhu W, Feng J, Su S, Zhao T. Microhabitat features determine the tadpole diversity in mountainous streams. Ecol Ind. 2021;126:107647.
Zhao T, Villeger S, Cucherousset J. Accounting for intraspecific diversity when examining relationships between non-native species and functional diversity. Oecologia. 2019;189:171–83.
Villeger S, Mason NWH, Mouillot D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology. 2008;89:2290–301.
Leduc A, da Silva EM, Rosenfeld JS. Effects of species vs. functional diversity: understanding the roles of complementarity and competition on ecosystem function in a tropical stream fish assemblage. Ecol Ind. 2015;48:627–35.
Matuoka MA, Benchimol M, JMd A-R, Morante-Filho JC. Effects of anthropogenic disturbances on bird functional diversity: a global meta-analysis.Ecological Indicators. 2020;116.
Whiles MR, Lips KR, Pringle CM, Kilham SS, Bixby RJ, Brenes R, Connelly S, Colon-Gaud JC, Hunte-Brown M, Huryn AD, et al. The effects of amphibian population declines on the structure and function of neotropical stream ecosystems. Front Ecol Environ. 2006;4:27–34.
Schmidt K, Pearson RG, Alford RA, Puschendorf R. Tadpole species have variable roles in litter breakdown, sediment removal, and nutrient cycling in a tropical stream. Freshw Sci. 2019;38:103–12.
Ramamonjisoa N, Natuhara Y. Contrasting effects of functionally distinct tadpole species on nutrient cycling and litter breakdown in a tropical rainforest stream. Freshw Biol. 2018;63:202–13.
Schmidt K, Blanchette ML, Pearson RG, Alford RA, Davis AM. Trophic roles of tadpoles in tropical australian streams. Freshw Biol. 2017;62:1929–41.
Jordani MX, Mouquet N, Casatti L, Menin M, Rossa-Feres DD, Albert CH. Intraspecific and interspecific trait variability in tadpole meta-communities from the brazilian Atlantic rainforest. Ecol Evol. 2019;9:4025–37.
Zhao T, Li C, Wang XY, Xie F, Jiang JP. Unraveling the relative contribution of inter- and intrapopulation functional variability in wild populations of a tadpole species. Ecol Evol. 2017;7:4726–34.
Iwai N, Pearson RG, Alford RA. Shredder-tadpole facilitation of leaf litter decomposition in a tropical stream. Freshw Biol. 2009;54:2573–80.
Zhao T, Wang XY, Wang XG, Wang SS, Chen YH, Jiang JP. Effects of urea on behavior and functional traits of Asiatic toad (Bufo gargarizans) tadpoles. Aquat Ecol. 2019;53:9–19.
Leao-Pires TA, Luiz AM, Sawaya RJ. The complex roles of space and environment in structuring functional, taxonomic and phylogenetic beta diversity of frogs in the Atlantic Forest. PLoS ONE. 2018;13:e0196066.
Strauss A, Reeve E, Randrianiaina R-D, Vences M, Glos J. The world’s richest tadpole communities show functional redundancy and low functional diversity: ecological data on Madagascar’s stream-dwelling amphibian larvae. BMC Ecol. 2010;10:12.
Sun Z, Zhu W, Zhu W, Zhao C, Liao C, Zou B, Xu D, Fan W, Su S, Jiang J, Zhao T. Spatiotemporal patterns of anuran functional diversity in temperate montane forests. Zoological Res. 2021;42:412–6.
Takatsu K, Kishida O. An offensive predator phenotype selects for an amplified defensive phenotype in its prey. Evol Ecol. 2012;27:1–11.
Jamil T, Ter Braak CJ. Generalized linear mixed models can detect unimodal species-environment relationships. PeerJ. 2013;1:95.
Verberk W, van der Velde G, Esselink H. Explaining abundance-occupancy relationships in specialists and generalists: a case study on aquatic macroinvertebrates in standing waters. J Anim Ecol. 2010;79:589–601.
Su G, Villéger S, Brosse S, Leprieur F. Morphological diversity of freshwater fishes differs between realms, but morphologically extreme species are widespread. Glob Ecol Biogeogr. 2019;28:211–21.
van der Merwe S, Greve M, Olivier B, le Roux PC. Testing the role of functional trait expression in plant-plant facilitation. Funct Ecol. 2021;35:255–65.
Zhao T, Grenouillet G, Pool T, Tudesque L, Cucherousset J. Environmental determinants of fish community structure in gravel pit lakes. Ecol Freshw Fish. 2016;25:412–21.
Bai H, Xu Q, Shi S, Zeng W, Tan L, Peng P. Radial growth of Abies fabri response to climate change along an altitudinal transect in the Mount Emei. J Cent South Univ Forestry Technol. 2020;40:128–38.
Tang CQ. Forest vegetation as related to climate and soil conditions at varying altitudes on a humid subtropical mountain, Mount Emei, Sichuan, China. Ecol Res. 2005;21:174–80.
Wang XY, Yang SN, Guo CP, Tang K, Jiang JP, Hu JH. Amphibian diversity and conservation along an elevational gradient on Mount Emei, southwestern China. Amphibian & Reptile Conservation. 2020;14:46–56.
Gosner KL. A simplified table for staging Anuran embryos and larvae with notes on ldentification. Herpetologists’ Leag. 1960;16:183–90.
Simpkins CA, Shuker JD, Lollback GW, Castley JG, Hero JM. Environmental variables associated with the distribution and occupancy of habitat specialist tadpoles in naturally acidic, oligotrophic waterbodies. Austral Ecol. 2014;39:95–105.
Thomas A, Das S, Manish K. Influence of stream habitat variables on distribution and abundance of tadpoles of the endangered purple frog, Nasikabatrachus sahyadrensis (Anura: Nasikabatrachidae). J Asia-Pacific Biodivers. 2019;12:144–51.
Xu D, Wu X, Zhao C, Zhu W, Zhao T, Su S. Functional approach reveals niche differentiation between two coexisting tadpoles: Leptobrachium boringii and Quasipaa boulengeri. Acta Ecol Sin. 2020;40:43–50.
Jung V, Violle C, Mondy C, Hoffmann L, Muller S. Intraspecific variability and trait-based community assembly. J Ecol. 2010;98:1134–40.
Manna LR, Villeger S, Rezende CF, Mazzoni R. High intraspecific variability in morphology and diet in tropical stream fish communities. Ecol Freshw Fish. 2019;28:41–52.
Zhao T, Villeger S, Lek S, Cucherousset J. High intraspecific variability in the functional niche of a predator is associated with ontogenetic shift and individual specialization. Ecol Evol. 2014;4:4649–57.
Rudolf VHW. A multivariate approach reveals diversity of ontogenetic niche shifts across taxonomic and functional groups. Freshw Biol. 2020;65:745–56.
Zhao E. An introduction of embryo and development stages of tadpoles. Bull Biol. 1990;1:13–5.
Mouillot D, Graham NAJ, Villeger S, Mason NWH, Bellwood DR. A functional approach reveals community responses to disturbances. Trends Ecol Evol. 2013;28:167–77.
Mouillot D, Culioli JM, Pelletier D, Tomasini JA. Do we protect biological originality in protected areas? A new index and an application to the Bonifacio Strait Natural Reserve. Biol Conserv. 2008;141:1569–80.
Eidietis L. The tactile-stimulated startle response of tadpoles: Acceleration performance and its relationship to the anatomy of wood frog (Rana sylvatica), bullfrog (Rana catesbeiana), and American toad (Bufo americanus) tadpoles. J Exp Zool A Comp Exp Biol. 2006;305(4):348–62.
Villeger S, Miranda JR, Hernandez DF, Mouillot D. Contrasting changes in taxonomic vs. functional diversity of tropical fish communities after habitat degradation. Ecol Appl. 2010;20:1512–22.
Anderson DR, Burnham KP. Avoiding pitfalls when using information-theoretic methods. J Wildl Manage. 2002;66:912–8.
Anderson DR. Model based inference in the life sciences: a primer on evidence. New York: Springer; 2008.
R Core Team. : R: A Language and Environment for Statistical Computing. 2022.
Oksanen J. Multivariate analysis of ecological communities in R: vegan tutorial. r package version 2010.
Bates D, Machler M, Bolker BM, Walker SC. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Walsh CJ, Papas PJ, Crowther D, Yoo J. Stormwater drainage pipes as a threat to a stream-dwelling amphipod of conservation significance, Austrogammarus australis, in southeastern Australia. Biodivers Conserv. 2004;13:781–93.
Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, Blomberg SP, Webb CO. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.
Azizi E, Landberg T, Wassersug RJ. Vertebral function during tadpole locomotion. Zoology. 2007;110:290–7.
Jiang Y, Wang J, Hu S, Jia Q. Breeding technology of tadpole stage of Rana spinosa David. Hebei Fishery. 2013;7:34–5.
Kopp K, Wachlevski M, Eterovick PC. Environmental complexity reduces tadpole predation by water bugs. Can J Zool. 2006;84:136–40.
Fei L, Ye C, Jiang J. Colored atlas of chinese amphibians and their distributions. ChengDu: Sichuan Science and Technology Press; 2012.
Mason NWH, Mouillot D, Lee WG, Wilson JB. Functional richness, functional evenness and functional divergence: the primary components of functional diversity. Oikos. 2005;111:112–8.
Stoler AB, Relyea RA. Reviewing the role of plant litter inputs to forested wetland ecosystems: leafing through the literature. Ecol Monogr. 2020;90:23.
Mouillot D, Villeger S, Parravicini V, Kulbicki M, Arias-Gonzalez JE, Bender M, Chabanet P, Floeter SR, Friedlander A, Vigliola L, Bellwood DR. Functional over-redundancy and high functional vulnerability in global fish faunas on tropical reefs. Proc Natl Acad Sci. 2014;111:13757–62.
Parravicini V, Villeger S, McClanahan TR, Arias-Gonzalez JE, Bellwood DR, Belmaker J, Chabanet P, Floeter SR, Friedlander AM, Guilhaumon F, et al. Global mismatch between species richness and vulnerability of reef fish assemblages. Ecol Lett. 2014;17:1101–10.
Clavel J, Julliard R, Devictor V. Worldwide decline of specialist species: toward a global functional homogenization? Front Ecol Environ. 2011;9:222–8.
de Melo LSO, Goncalves-Souza T, Garey MV, de Cerqueira D. Tadpole species richness within lentic and lotic microhabitats: an interactive influence of environmental and spatial factors. Herpetological J. 2017;27:339–45.
Peltzer PM, Lajmanovich RC. Anuran tadpole assemblages in riparian areas of the middle Parana River, Argentina. Biodivers Conserv. 2004;13:1833–42.
Bogota-Gregory JD, Lima FCT, Correa SB, Silva-Oliveira C, Jenkins DG, Ribeiro FR, Lovejoy NR, Reis RE, Crampton WGR. Biogeochemical water type influences community composition, species richness, and biomass in megadiverse amazonian fish assemblages. Sci Rep. 2020;10:15.
We are grateful to the Editor Sebastian Steinfartz and two anonymous Reviewers who provided insightful comments on an earlier version of the manuscript. We thank Janak Raj Khatiwada, Quangui Dan, Weiyi Song, Dan Xu, Da Kang, Jiwei Xia, Xiaoyi Wang, Shengnan Yang, Xuecheng Guo, Wenyan Zhang, Yinmeng Hou, Luyao Xiao, and Chunpeng Guo for their help in the field. We also thank Jin Zhou, Yinmeng Hou, and Shengchao Shi for laboratory assistance.
This study was supported by the Second Tibetan Plateau Scientific Expedition and Research Program (STEP), Grant No. 2019QZKK0501, the National Natural Science Foundation of China (31700353), and China Biodiversity Observation Networks (Sino BON). The Animal Care and Use Committee of Chengdu Institute of Biology, CAS provided full approval for this study (Number: CIB2011032201). Fieldwork was approved by the Management Office of the Emei Mountain Management Committee (project number: EMS-201304006).
Ethics approval and consent to participate
The animal study protocol was approved by the Institutional Ethics Committee of Animal Ethical and Welfare Committee of Chengdu Institute of Biology, Chinese Academy of Sciences (CIB20171105).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lan, J., Sun, Z., Feng, J. et al. Unraveling the importance of functionally extreme tadpole types to functional diversity: a case study in temperate montane streams. Front Zool 20, 7 (2023). https://doi.org/10.1186/s12983-023-00485-0
- Functional vulnerability
- Ecomorphological traits
- Functionally extreme species
- Microhabitat variables