Species richness and trait composition of butterfly assemblages change along an altitudinal gradient

Community metabarcoding reveals the relative role of environmental filtering and spatial processes in metacommunity dynamics of soil microarthropods across a mosaic of montane forests

Disentangling the relative role of environmental filtering and spatial processes in driving metacommunity structure across mountainous regions remains challenging, as the way we quantify spatial connectivity in topographically and environmentally heterogeneous landscapes can influence our perception of which process predominates. More empirical data sets are required to account for taxon- and context-dependency, but relevant research in understudied areas is often compromised by the taxonomic impediment. Here we used haplotype-level community DNA metabarcoding, enabled by stringent filtering of amplicon sequence variants (ASVs), to characterize metacommunity structure of soil microarthropod assemblages across a mosaic of five forest habitats on the Troodos mountain range in Cyprus. We found similar β diversity patterns at ASV and species (OTU, operational taxonomic unit) levels, which pointed to a primary role of habitat filtering resulting in the existence of largely distinct metacommunities linked to different forest types. Within-habitat turnover was correlated to topoclimatic heterogeneity, again emphasizing the role of environmental filtering. However, when integrating landscape matrix information for the highly fragmented Quercus alnifolia habitat, we also detected a major role of spatial isolation determined by patch connectivity, indicating that stochastic and niche-based processes synergistically govern community assembly. Alpha diversity patterns varied between ASV and OTU levels, with OTU richness decreasing with elevation and ASV richness following a longitudinal gradient, potentially reflecting a decline of genetic diversity eastwards due to historical pressures. Our study demonstrates the utility of haplotype-level community metabarcoding for characterizing metacommunity structure of complex assemblages and improving our understanding of biodiversity dynamics across mountainous landscapes worldwide.

Alpine butterflies want to fly high: Species and communities shift upwards faster than their host plants

Despite sometimes strong codependencies of insect herbivores and plants, the responses of individual taxa to accelerating climate change are typically studied in isolation. For this reason, biotic interactions that potentially limit species in tracking their preferred climatic niches are ignored. Here, we chose butterflies as a prominent representative of herbivorous insects to investigate the impacts of temperature changes and their larval host plant distributions along a 1.4-km elevational gradient in the German Alps. Following a sampling protocol of 2009, we revisited 33 grassland plots in 2019 over an entire growing season. We quantified changes in butterfly abundance and richness by repeated transect walks on each plot and disentangled the direct and indirect effects of locally assessed temperature, site management, and larval and adult food resource availability on these patterns. Additionally, we determined elevational range shifts of butterflies and host plants at both the community and species level. Comparing the two sampled years (2009 and 2019), we found a severe decline in butterfly abundance and a clear upward shift of butterflies along the elevational gradient. We detected shifts in the peak of species richness, community composition, and at the species level, whereby mountainous species shifted particularly strongly. In contrast, host plants showed barely any change, neither in connection with species richness nor individual species shifts. Further, temperature and host plant richness were the main drivers of butterfly richness, with change in temperature best explaining the change in richness over time. We concluded that host plants were not yet hindering butterfly species and communities from shifting upwards. However, the mismatch between butterfly and host plant shifts might become a problem for this very close plant-herbivore relationship, especially toward higher elevations, if butterflies fail to adapt to new host plants. Further, our results support the value of conserving traditional extensive pasture use as a promoter of host plant and, hence, butterfly richness.

Diversity of Lacewings (Neuroptera) in an Altitudinal Gradient of the Tacaná Volcano, Southern Mexico

Neuroptera is an order of insects with a moderate diversity of species numbers yet a high between-family morphological diversity, which has a significant ecological role as a predator. However, there are few studies focused on describing changes in species diversity along environmental gradients. We evaluated changes in the alpha and beta diversity of species and the higher taxa in Neuroptera communities in the Tacaná Volcano in southern Mexico. Five sites each at different altitudes were studied through systematic annual sampling. The taxonomic and phylogenetic alpha diversity were analyzed, as well as the beta diversity and its components, species turnover and nestedness. The alpha diversity had two trends: (1) decreased standardized richness and taxonomic distinctness with increasing altitude, and (2) increased estimated richness and species diversity at intermediate altitudes. The highest turnover values for species, as well as for supra-specific taxa, were recorded at sites with lower altitudes. The highest total beta diversity value was recorded at elevations above 3000 m, whereas the highest number of species and supra-specific taxa were observed at sites between 600 and 2000 m, with an evident decrease above 3000 m. The type of vegetation and environmental conditions may be influencing the decrease in diversity toward higher elevations, which could explain the niche specialization of Neuroptera species to particular sites within the gradient. These results highlight the need to study the environmental factors and their effects on species composition along an elevation gradient.

Species richness, abundance, distributional pattern and trait composition of butterfly assemblage change along an altitudinal gradient in the Gulmarg region of Jammu & Kashmir, India

Despite enormous diversity, abundance, and role in ecosystem processes, little is known about how butterflies differ across altitudinal gradients. For this, butterfly communities were investigated along an altitudinal gradient of 2700–3200 m a.s.l, along the Gulmarg region of Jammu & Kashmir, India. We aimed to determine how the altitudinal gradient and environmental factors affect the butterfly diversity and abundance. Our findings indicate that species richness and diversity are mainly affected by the synergism between climate and vegetation. Alpha diversity indices showed that butterfly communities were more diverse at lower elevations and declined significantly with increase in elevation. Overall, butterfly abundance and diversity is stronger at lower elevations and gradually keep dropping towards higher elevations because floristic diversity decreased on which butterflies rely for survival and propagation. A total of 2023 individuals of butterflies were recorded belonging to 40 species, represented by 27 genera and 05 families. Six survey sites (S I- S VI) were assessed for butterfly diversity from 2018 to 2020 in the Gulmarg region of Jammu & Kashmir. Across the survey, Nymphalidae was the most dominant family represented by 16 genera and 23 species, while Papilionidae and Hesperiidae were least dominant represented by 01 genera and 01 species each. Among the six collection sites selected, Site I was most dominant, represented by 16 genera and 21 species, while Site VI was least dominant, represented by 04 genera and 04 species.

Plant traits mediate the effects of climate on phytophagous beetle diversity on Mt. Kilimanjaro

Patterns of insect diversity along elevational gradients are well described in ecology. However, it remains little tested how variation in the quantity, quality, and diversity of food resources influence these patterns. Here we analyzed the direct and indirect effects of climate, food quantity (estimated by net primary productivity), quality (variation in the specific leaf area index, leaf nitrogen to phosphorus and leaf carbon to nitrogen ratio), and food diversity (diversity of leaf traits) on the species richness of phytophagous beetles along the broad elevation and land use gradients of Mt. Kilimanjaro, Tanzania. We sampled beetles at 65 study sites located in both natural and anthropogenic habitats, ranging from 866 to 4,550 m asl. We used path analysis to unravel the direct and indirect effects of predictor variables on species richness. In total, 3,154 phytophagous beetles representing 19 families and 304 morphospecies were collected. We found that the species richness of phytophagous beetles was bimodally distributed along the elevation gradient with peaks at the lowest (˜866 m asl) and upper mid-elevations (˜3,200 m asl) and sharply declined at higher elevations. Path analysis revealed temperature- and climate-driven changes in primary productivity and leaf trait diversity to be the best predictors of changes in the species richness of phytophagous beetles. Species richness increased with increases in mean annual temperature, primary productivity, and with increases in the diversity of leaf traits of local ecosystems. Our study demonstrates that, apart from temperature, the quantity and diversity of food resources play a major role in shaping diversity gradients of phytophagous insects. Drivers of global change, leading to a change of leaf traits and causing reductions in plant diversity and productivity, may consequently reduce the diversity of herbivore assemblages.

Species richness and beta diversity patterns of multiple taxa along an elevational gradient in pastured grasslands in the European Alps

To understand how diversity is distributed in space is a fundamental aim for optimizing future species and community conservation. We examined in parallel species richness and beta diversity components of nine taxonomic groups along a finite space, represented by pastured grasslands along an elevational gradient. Beta diversity, which is assumed to bridge local alpha diversity to regional gamma diversity was partitioned into the two components turnover and nestedness and analyzed at two levels: from the lowest elevation to all other elevations, and between neighboring elevations. Species richness of vascular plants, butterflies, beetles, spiders and earthworms showed a hump-shaped relationship with increasing elevation, while it decreased linearly for grasshoppers and ants, but increased for lichens and bryophytes. For most of the groups, turnover increased with increasing elevational distance along the gradient while nestedness decreased. With regard to step-wise beta diversity, rates of turnover or nestedness did not change notably between neighboring steps for the majority of groups. Our results support the assumption that species communities occupying the same habitat significantly change along elevation, however transition seems to happen continuously and is not detectable between neighboring steps. Our findings, rather than delineating levels of major diversity losses, indicate that conservation actions targeting at a preventive protection for species and their environment in mountainous regions require the consideration of entire spatial settings.

Diatom Species Richness in Swiss Springs Increases with Habitat Complexity and Elevation

Understanding the drivers of species richness gradients is a central challenge of ecological and biodiversity research in freshwater science. Species richness along elevational gradients reveals a great variety of patterns. Here, we investigate elevational changes in species richness and turnover between microhabitats in near-natural spring habitats across Switzerland. Species richness was determined for 175 subsamples from 71 near-natural springs, and Poisson regression was applied between species richness and environmental predictors. Compositional turnover was calculated between the different microhabitats within single springs using the Jaccard index based on observed species and the Chao index based on estimated species numbers. In total, 539 diatom species were identified. Species richness increased monotonically with elevation. Habitat diversity and elevation explaining some of the species richness per site. The Jaccard index for the measured compositional turnover showed a mean similarity of 70% between microhabitats within springs, whereas the Chao index which accounts for sampling artefacts estimated a turnover of only 37%. Thus, the commonly applied method of counting 500 valves led to an undersampling of the rare species and might need to be reconsidered when assessing diatom biodiversity.

Butterfly phenology in Mediterranean mountains using space-for-time substitution

Inferring species' responses to climate change in the absence of long-term time series data is a challenge, but can be achieved by substituting space for time. For example, thermal elevational gradients represent suitable proxies to study phenological responses to warming. We used butterfly data from two Mediterranean mountain areas to test whether mean dates of appearance of communities and individual species show a delay with increasing altitude, and an accompanying shortening in the duration of flight periods. We found a 14-day delay in the mean date of appearance per kilometer increase in altitude for butterfly communities overall, and an average 23-day shift for 26 selected species, alongside average summer temperature lapse rates of 3°C per km. At higher elevations, there was a shortening of the flight period for the community of 3 days/km, with an 8.8-day average decline per km for individual species. Rates of phenological delay differed significantly between the two mountain ranges, although this did not seem to result from the respective temperature lapse rates. These results suggest that climate warming could lead to advanced and lengthened flight periods for Mediterranean mountain butterfly communities. However, although multivoltine species showed the expected response of delayed and shortened flight periods at higher elevations, univoltine species showed more pronounced delays in terms of species appearance. Hence, while projections of overall community responses to climate change may benefit from space-for-time substitutions, understanding species-specific responses to local features of habitat and climate may be needed to accurately predict the effects of climate change on phenology.

Leaf traits mediate changes in invertebrate herbivory along broad environmental gradients on Mt. Kilimanjaro, Tanzania

Temperature, primary productivity, plant functional traits, and herbivore abundances are considered key predictors of leaf herbivory but their direct and indirect contributions to community-level herbivory are not well understood along broad climatic gradients. Here, we determined elevational herbivory patterns and used a path analytical approach to disentangle the direct and indirect effects of climate, land use, net primary productivity (NPP), herbivore abundance, and plant functional traits on community-level invertebrate herbivory along the extensive elevational and land use gradients at Mt. Kilimanjaro, Tanzania. We recorded standing leaf herbivory caused by leaf chewers, leaf miners and leaf gallers on 55 study sites distributed in natural and anthropogenic habitats along a 3,060 m elevation gradient. We related the total community-level herbivory to climate (temperature and precipitation), NPP, plant functional traits (specific leaf area, leaf carbon-to-nitrogen [CN] ratio and leaf nitrogen-to-phosphorus [NP] ratio) and herbivore abundances. Leaf herbivory ranged from 5% to 11% along the elevation gradient. Total leaf herbivory showed unimodal pattern in natural habitats but a strongly contrasting bimodal pattern in anthropogenic habitats. We also detected some variation in the patterns of leaf herbivory along environmental gradients across feeding guilds with leaf chewers being responsible for a disproportionally large part of herbivory. Path analyses indicated that the variation in leaf herbivory was mainly driven by changes in leaf CN and NP ratios which were closely linked to changes in NPP in natural habitats. Similarly, patterns of leaf herbivory in anthropogenic habitats were best explained by variation in leaf CN ratios and a negative effect of land use. Our study elucidates the strong role of leaf nutrient stoichiometry and its linkages to climate and NPP for explaining the variation in leaf herbivory along broad climatic gradients. Furthermore, the study suggests that climatic changes and nutrient inputs in the course of land use change may alter leaf herbivory and consequently energy and nutrient fluxes in terrestrial habitats.

Resource availability drives trait composition of butterfly assemblages

How species respond to environmental change is a fundamental question in ecology and species traits can help to tackle this question. In this study, we analyze how the functional structure of species assemblages changes with selected environmental variables along an elevational gradient. In particular, we used species traits of local butterfly communities (body size, voltinism, overwintering stages, and host specificity) in a national nature reserve in China to assess the impacts of temperature, net primary productivity, and land use. Our results show that productivity, measured as NDVI, had a stronger influence on the functional community structure of butterflies than temperature. Within the butterfly assemblages, net primary productivity mainly affected body size and supported few but large species. Length of vegetation period demonstrated dominating effects on the functional structure of local butterfly assemblages. However, an observed increase in dietary generalists with longer vegetation periods contradicted expectations based on niche breadth hypothesis, that more stable conditions should favor specialists. Furthermore, the general positive impact of vegetation period on species abundances differed considerably among functional groups. Only the group containing species hibernating as egg decreased with the length of vegetation period. Our results suggest that trait associations are instructive to explain environment-herbivore relationships, that resource availability can predominantly influence the functional composition of herbivore assemblages, and that conservation priority should be given to specialist butterfly species overwintering as egg, especially in the face of global warming.

Trait-based ecology of terrestrial arthropods

In focusing on how organisms' generalizable functional properties (traits) interact mechanistically with environments across spatial scales and levels of biological organization, trait-based approaches provide a powerful framework for attaining synthesis, generality and prediction. Trait-based research has considerably improved understanding of the assembly, structure and functioning of plant communities. Further advances in ecology may be achieved by exploring the trait-environment relationships of non-sessile, heterotrophic organisms such as terrestrial arthropods, which are geographically ubiquitous, ecologically diverse, and often important functional components of ecosystems. Trait-based studies and trait databases have recently been compiled for groups such as ants, bees, beetles, butterflies, spiders and many others; however, the explicit justification, conceptual framework, and primary-evidence base for the burgeoning field of 'terrestrial arthropod trait-based ecology' have not been well established. Consequently, there is some confusion over the scope and relevance of this field, as well as a tendency for studies to overlook important assumptions of the trait-based approach. Here we aim to provide a broad and accessible overview of the trait-based ecology of terrestrial arthropods. We first define and illustrate foundational concepts in trait-based ecology with respect to terrestrial arthropods, and justify the application of trait-based approaches to the study of their ecology. Next, we review studies in community ecology where trait-based approaches have been used to elucidate how assembly processes for terrestrial arthropod communities are influenced by niche filtering along environmental gradients (e.g. climatic, structural, and land-use gradients) and by abiotic and biotic disturbances (e.g. fire, floods, and biological invasions). We also review studies in ecosystem ecology where trait-based approaches have been used to investigate biodiversity-ecosystem function relationships: how the functional diversity of arthropod communities relates to a host of ecosystem functions and services that they mediate, such as decomposition, pollination and predation. We then suggest how future work can address fundamental assumptions and limitations by investigating trait functionality and the effects of intraspecific variation, assessing the potential for sampling methods to bias the traits and trait values observed, and enhancing the quality and consolidation of trait information in databases. A roadmap to guide observational trait-based studies is also presented. Lastly, we highlight new areas where trait-based studies on terrestrial arthropods are well positioned to advance ecological understanding and application. These include examining the roles of competitive, non-competitive and (multi-)trophic interactions in shaping coexistence, and macro-scaling trait-environment relationships to explain and predict patterns in biodiversity and ecosystem functions across space and time. We hope this review will spur and guide future applications of the trait-based framework to advance ecological insights from the most diverse eukaryotic organisms on Earth.

Elevational gradient and human effects on butterfly species richness in the French Alps

We examined how butterfly species richness is affected by human impact and elevation, and how species ranges are distributed along the elevational gradient (200–2700 m) in the Isère Department (French Alps). A total of 35,724 butterfly observations gathered in summer (May–September) between 1995 and 2015 were analyzed. The number of estimated species per 100‐m elevational band was fitted to the elevational gradient using a generalized additive model. Estimations were also performed on a 500 m × 500 m grid at low altitude (200–500 m) to test for the human impact on species richness using generalized least squares regression models. Each species elevational range was plotted against the elevational gradient. Butterfly richness along the elevational gradient first increased (200–500 m) to reach a maximum of 150 species at 700 m and then remained nearly constant till a sharp decrease after 1900 m, suggesting that after some temperature threshold, only few specialized species can survive. At low elevation, urbanization and arable lands had a strongly negative impact on butterfly diversity, which was buffered by a positive effect of permanent crops. Butterfly diversity is exceptionally high (185 species) in this alpine department that represents less than 5% of the French territory and yet holds more than 70% of all the Rhopalocera species recorded in France. Both climate and habitat shape the distribution of species, with a negative effect of anthropization at low altitude and strong climatic constraints at high altitude.

Change in Anopheles richness and composition in response to artificial flooding during the creation of the Jirau hydroelectric dam in Porto Velho, Brazil

BACKGROUND

Anopheles mosquitoes are the only vectors of human malaria. Anopheles species use standing water as breeding sites. Human activities, like the creation of an artificial lake during the implementation of hydroelectric power plants, lead to changes in environmental characteristics and, therefore, may changes the species richness and composition of Anopheles mosquitoes. The aim of the present study was to verify whether or not there is an association between the artificial flooding resulting from the construction of the Jirau hydroelectric power plant, and the richness and composition of anophelines.

METHODS

Mosquitoes samples were obtained monthly from the Jirau hydroelectric power plant area located at Porto Velho, Rondônia State, using Human Landing Catch (06:00-10:00 PM). Mosquitoes collected were transported to Laboratório de Entomologia Médica FIOCRUZ-RO where they were identified until species using dichotomous key.

RESULTS

A total of 6347 anophelines belonging to eight different species were collected. The anophelines species richness was significantly lower during the first flooding stage. Differences in anophelines species composition were found when comparing the first flooding stage with the other stages. Furthermore, the mean number of Anopheles darlingi, the main vector of malaria in the region, increases during the first and the third flooding stages.

CONCLUSIONS

The continual monitoring of these vectors during the late operational phase may be useful in order to understand how anophelines will behave in this area.

Life History Traits Reflect Changes in Mediterranean Butterfly Communities Due to Forest Encroachment

The biodiversity of the Southern Balkans, part of the Mediterranean global biodiversity hot-spot, is threatened by land use intensification and abandonment, the latter causing forest encroachment of formerly open habitats. We investigated the impact of forest encroachment on butterfly species richness, community species composition and the representation of life history traits by repeated seasonal visits of 150 one-hectare sites in five separate regions in three countries—Greece, Bulgaria, and the Republic of Macedonia (FYROM—the Former Yugoslav Republic of Macedonia)— 10 replicates for each habitat type of grasslands, open formations and scrub forest within each region. Grasslands and open formations sites hosted in average more species and more red-listed species than scrub forest, while no pattern was found for numbers of Mediterranean species. As shown by ordination analyses, each of the three habitat types hosted distinct butterfly communities, with Mediterranean species inclining either towards grasslands or open formations. Analysing the representation of life history traits revealed that successional development from grasslands and open formations towards scrub forest shifts the community composition towards species overwintering in earlier stages, having fewer generations per year, and inhabiting large European or Eurosiberian (e.g. northern) ranges; it decreases the representation of Mediterranean endemics. The loss of grasslands and semi-open formations due to forest encroachment thus threatens exactly the species that should be the focus of conservation attention in the Mediterranean region, and innovative conservation actions to prevent ongoing forest encroachment are badly needed.