Next Article in Journal
Multi-Omics Reveals the Effect of Population Density on the Phenotype, Transcriptome and Metabolome of Mythimna separata
Previous Article in Journal
Comparison of Methyl Bromide and Ethyl Formate for Fumigation of Snail and Fly Pests of Imported Orchids
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 14
Issue 1
10.3390/insects14010067
insects-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

Environmental Factors Drive Chalcid Body Size Increases with Altitudinal Gradients for Two Hyper-Diverse Taxa

by
Ning Kang
1,2,
Hongying Hu
1,2,*,
Zengqian Huang
1,2,
Shungang Luo
1,2 and
Shuhan Guo
1,2
1
College of Life Science and Technology, Xinjiang University, Urumqi 830049, China
2
Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, Urumqi 830046, China
*
Author to whom correspondence should be addressed.
Insects 2023, 14(1), 67; https://doi.org/10.3390/insects14010067
Submission received: 1 November 2022 / Revised: 23 December 2022 / Accepted: 8 January 2023 / Published: 10 January 2023
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

The size of an organism is closely correlated with its physiological and ecological characteristics and strongly influences its fitness. Bergmann’s rule, originally widely applied in homeotherms, states that individuals living in a colder environment are larger than those living in a warmer region. Similar geographical patterns are also found in various groups of poikilotherms, but are still controversial, especially in Hymenoptera. In this investigation, we found a significant upward trend in body size with increasing elevation in two tiny groups of Chalcids (Pteromalidae and Eulophidae). The temperature and precipitation play a crucial role in their size variation. This result casts light on the environmental adaptation of parasitoids.

Abstract

Body size is the most essential feature that significantly correlates with insects’ longevity, fecundity, metabolic rate, and sex ratio. Numerous biogeographical rules have been proposed to illustrate the correlation between the body sizes of different taxa and corresponding geographical or environmental factors. Whether the minute and multifarious chalcids exhibit a similar geographical pattern is still little known. In this research, we analyzed morphological data from 2953 specimens worldwide, including the two most abundant and diverse taxa (Pteromalidae and Eulophidae), which are both composed of field-collected and BOLD system specimens. We examined forewing length as a surrogate of body size and analyzed the average size separately for males and females using two methods (species and assemblage-based method). To verify Bergmann’s rule, we included temperature, precipitation, wind speed and solar radiation as explanatory variables in a generalized linear model to analyze the causes of the size variation. We found that there was an increasing trend in the body size of Pteromalidae and Eulophidae with altitude. The optimal Akaike information criterion (AIC) models showed that larger sizes are significantly negatively correlated with temperature and positively correlated with precipitation, and the possible reasons for this variation are discussed and analyzed.

1. Introduction

Body size is one of the most remarkable characteristics of all organisms, since it is significantly correlated with the biological, physiological, and ecological traits involving longevity, fecundity, metabolic and reproduction rate, and population size [1,2,3]. Bergmann’s rule [4], a widely recognized ecogeographic rule, states that animals living in colder climates are larger than those living in a warmer regions. The geographical pattern has been broadly verified in mammals [5] and birds [6], which have broad ecological niches ranging from tropical to arctic, and their body size is positively correlated with latitude.
Bergmann’s rule was originally widely applied in homeotherms, then also verified in various poikilotherms, especially variations in arthropod sizes at both inter- and intraspecific levels along altitude or latitude gradients, but general patterns in the geographical trend of body size are still contentious. Indeed, several insect clades follow the converse Bergmann’s rule [7], such as Onthophagus spp.(Scarabaeidae) [8], Dichroplus Pratensis (Acrididae) [9], and Paropsis atomaria (Chrysomelidae) [10]. Conversely, the body size of Zygogramma bicolorata (Chrysomelidae) [11], weevils [12], moths [13], butterflies [14], stoneflies [15], and psyllids [16] generally increases with elevation. However, the body size of various insects lacks a clear relationship with altitude or latitude [17,18], and these contradictions are particularly apparent in Hymenoptera.
By summarizing the relevant literature (Table 1), we found that most of the investigations into Hymenoptera were focused on ants, bees, and other larger wasps. The geographical trends vary among species and different morphological parts, and some ant species even showed a nonlinear trend with elevation [19]. Relevant research on miniature taxa only involved Encyrtidae [20] and Mymaridae [21], and only conducted a simple comparison of body-size differences, so the knowledge of small-sized taxa is still limited. The forewing length and hind tibia length are two generally used surrogates for body length in laboratory experiments, but this method lacks validation in multiple wild species. Meanwhile, they are closely correlated with longevity, fecundity, reproductive duration, and oviposition index [22,23], which provide insights into the relationship between size and environmental factors.
The various geographical trends are generally caused by the combined effect of multiple environmental factors [24]. Based on the results of several models’ analyses, the explanatory mechanisms of these trends are usually correlated with temperature and precipitation in Diptera [25,26], Lepidoptera [13], Odonata [27,28], and Orthoptera [29]. Numerous hypotheses related to body-size control in parasitoids have been widely verified in laboratory experiments [30,31,32]. In consequence, it is critical to study these relationships under realistic natural conditions. As we all know, temperature significantly decreases with increasing elevation and is linearly correlated [33]; the most frequently verified hypotheses is the temperature-size rule [34,35], which states that a higher temperature may accelerate the development rate and result in adults with a smaller body size than those living in lower temperatures; therefore, insects living in lower temperatures tend to be relatively larger. Resources availability is also crucial to insect survival and development; insects tend to reduce their reproductive investment and increase energy storage to maintain survival under the stress of starvation [36]; thus, larger-sized insects have the potential to store more resources to resist that stress [37] and can also remain warm in cooler environments due to their lower surface-area-to-volume ratio. The wing-size polymorphism and variations are also closely related to altitude among flying insects. Wing loading and the wing aspect ratio are two common parameters for flight capability. The lower flight activity and larger wing size in a highland environment compensate for the thinner air encountered during flight [24]. However, it is still unknown whether these rules also apply to minute parasitic insects.
Table
Table 1. A review of the application of Bergmann’s rule in Hymenoptera insects based on the significant research. CB: Follow converse Bergmann’s rule, B: Follow Bergmann’s rule.
Table 1. A review of the application of Bergmann’s rule in Hymenoptera insects based on the significant research. CB: Follow converse Bergmann’s rule, B: Follow Bergmann’s rule.
TaxaGenus or SpeciesBergmann’s RuleData SourcesFactorsMeasured TreatReference
ApidaeBombusCBSpecimens from museumLatitudeThorax width[38]
ApidaeB. lucorum, B. magnus and B. cryptarumBSpecimens collected from Great BritainLatitudeThorax width[39]
ApidaeB. vancouverensisBSpecimens collected from the United StatesLatitudeInter-tegular span (ITS), Forewing area
Community-weighted mean of ITS, relative forewing length, relative hind tibia length, hair length at mesonotum, relative hair length at mesonotum		[40]
B. vosnesenskiiCBAltitude
Apoidea, Megachilidae, Andrenidae, Halictidae, ColletidaeMacrotera sp.,
Lasioglossum (Dialictus) sp., Lasioglossum (Lasioglossum) sp., Halictus sp.		B (only for Relative hair length)
CB		Specimens collected from north MexicoAltitude[41]
ApoideaAndrena, Dasypoda, Halictus, PanurginusBData from databaseLatitudeITS[42]
Bombus, Colletes, MelittaCB
ApoideaAt family levelB, species larger than 27.81 mg (dry weight)Specimens collected from Mediterranean; data from articles and museum collectionsTemperatureITS[43]
CB, species less than 27.81 mg
BraconidaeAt family levelCBLiterature recordedLatitudeBody length[44]
Ichneumonidae
BethylidaePristocera rufaCBLiterature recordedLatitudeBody length[20]
BraconinaeAt family levelCB
EncyrtidaeComperiella bifasciataCB
IchneumonidaeAt family levelB
Braconidae, IchneumonidaeAt family levelCB Specimens collected from Canada and data from BOLD systemLatitudeHind tibia lengths [45]
FormicidaeAnt species from 427 generaCBSpecimens collected from the United StatesAltitude/latitudeWeber’s length[46]
FormicidaeLeptothorax acervorumBSpecimens collected from AlpsAltitudeMaximum cephalic width[47]
FormicidaeAt family levelCBSpecimens collected from Australia, South AfricaAltitudeBody mass, critical thermal maximum (CTmax)[48]
FormicidaeLeptothorax acervorumBSpecimens collected throughout EuropeLatitudesThorax length[49]
MymaridaeAnaphes cultripennisBSpecimens collected from NorwayAltitudeBody length[21]
VespidaePolistes bahamensis, P. bellicosus, P. comanchus, P. dorsalis, P. exclamansBSpecimens from museumLatitudeBody length[50]
P. annularis, P. apachus, P. Carolina, P. flavus, P. fuscatus, P. metricus, P. dominula, P. aurifer, P. bahamensisCB
VespidaeAgelaia pallipesCBSpecimens collected from Santuario de Iguaque, ColombiaAltitudeHind femur length, mesosoma height, head width/length[51]
BForewing width/length
To address these gaps, we primarily focused on the geographical patterns of the forewing length of Pteromalidae and Eulophidae, based on a large-scale dataset in this analysis, which consists of field-collected samples from Northwestern China and data from the BOLD system. The main objectives of this research were to analyze the relationship between morphological characteristics and elevation; verify whether the body size increases with increasing elevation, as Bergmann’s rule states; investigate the dominant environmental factors contributing to this trend, such as temperature and precipitation; and offer alternative explanations for the relationship, including the temperature-size hypothesis and the resource’s availability hypothesis. We hope to provide insights that could improve the understanding of phenotypic plasticity in the future.

2. Materials and Methods

2.1. Study Area and Specimen Measurement

This study was based on a compiled dataset of 2953 specimens, which included 1406 Pteromalids (52 genera and 81 species) and 1547 Eulophids (54 genera and 92 species), including endemically common and cosmopolitan species. More detailed information and the composition of the dataset is shown in Table S1.
The field-collected specimens were primarily from Northwestern China (mainly covering the endemic habitats of Xinjiang, Gansu, Ningxia, and western Inner Mongolia). The area provides full-scale altitudinal gradient transect ranges from 220 m in the basin to 4336.9 m in the Altun mountain national nature reserve (Table S2). The specimens were collected using net-sweeping and yellow pan traps (malaise traps were used at two alpine stations) during July and August from 2017 to 2020, then labeled and identified definitely to the level of genus or species following a combination of classical taxonomic guides [52,53,54,55,56], and preserved in the Insect Collection of the College of Life Science and Technology, Xinjiang University, Urumqi, Xinjiang, China (ICXU). To increase the sample quantity and geographic coverage, we selected specimens from the BOLD system (http://www.boldsystems.org/. accessed on 1 February 2021), with a complete morphological graph (with scale bar), taxonomy, and collection site, to add to the dataset. The final dataset covers 217 collection sites worldwide (sites less than 10 km apart were merged into one site). We also extracted the GPS coordinates of each point for subsequent analysis. The measured specimens finally used for analysis were distributed across all the continents except the poles (Figure 1a).
We photographed each collected specimen using a Nikon SMZ25 stereomicroscope (Nikon, Japan) in conjunction with an NIS-Elements D 4.30.00 image-collecting system using 10× magnification to capture the entire body in a single image. To avoid deviation, for accurate measurements, and to increase the reliability of the data, any severely shriveled or damaged specimens were ignored and the average value of multiple measurements for each characteristic was calculated. We measured the body length (BL) from the top of the head to the end of the abdomen, excluding the ovipositor protrusion, the mesosoma length (MEL) from the anterior edge of the pronotum to the posterior margin of the propodeum, the forewing length (FL) from the tegula proximally to the outermost point of its distal end, and the hind tibia length (HTL) (Figure 1b). The forewing was spread on the microscope slides for accurate measurement. The wing area (WA) was measured by drawing the contour of the forewing, starting and ending at the tegula (Figure 1b). After summarizing the relevant literature, we found that the wing loading (WL) and wing aspect ratio (WAR) also have a certain relationship with altitude[13,57]. Since it is difficult to obtain accurate body mass data on parasitoids, the WL was replaced by MEL3/WA2 [58], and the WAR was evaluated by FL2/WA [59]. All the measurements were conducted using ImageJ/Fiji 1.46.

2.2. Environmental Data

Based on the previous predictions and hypotheses of the size-related Bergmann’s rule, wind speed has a certain effect on parasitoids’ flight ability [60], and the lower water-loss hypothesis for larger insects in drought climates [61,62], we only considered annual mean precipitation, temperature, wind speed and solar radiation to analyze their effects on the morphological traits. The environmental variables used in our study were obtained from WorldClim-Global Climate and Weather data (version 2.1 available at https://www.worldclim.org/data/worldclim21.html. accessed on 10 April 2022.), developed by Fick and Hijmans [63]. The four environmental factors for each site were extracted using the ‘surface tool’ in ArcGIS Spatial Analyst [64], based on the geographical coordinates of the specimens used in our study; the latter two variables were obtained by calculating the average value of 12 months of data; all the digital information was acquired at the spatial level resolution of 30 arcseconds.

2.3. Statistical Analysis

We identified the best morphological characteristics that can be used as a proxy for body length by using an OLS regression of seven morphological parameters from males and females of Pteromalidae and Eulophidae. After establishing that the FL is more closely related to the BL than other characteristics (Pteromalidae: Cor = 0.931; Eulophidae: Cor = 0.953), we analyzed the geographical patterns of FL with altitude for both sexes, using two methods at the family level: (1). In the species-based method, the arithmetic means of FL were calculated separately for each species per location, and species with fewer than three records were excluded; a total of 2134 (1233 Eulophidae and 901 Pteromalidae) specimens were used in this analysis. (2). For the assemblage-based method, the arithmetic means of FL were calculated for all specimens from the same location, ignoring the species; some sites with fewer than three records were excluded; 2706 (1420 Eulophidae and 1286 Pteromalidae) specimens were analyzed by this method. This method aims to reflect the complete situation of populations at a specific location and is more ecologically meaningful.
The generalized linear model (GLM) is more inclusive, allowing for dependent variables to be non-normally distributed. We used the FL of the different sexes as dependent variables in the GLM for the species-based and assemblage-based methods, and the corresponding environmental factors of each site were used as the explanatory variables. A gamma distribution of errors and a log link function were used in the GLM, while the stepAIC () function was applied to avoid collinearity and simplify the models according to the Akaike information criterion (AIC) values.
We analyzed the relationship between the FL and altitude for six dominant subfamilies with the largest sample sizes. In the genus-level analysis, we accessed the relationship between the altitude (dependent variable) and the FL of 16 widely spread genera (sample size: at least 30 specimens) from Pteromalidae and Eulophidae (independent variables), taking into account the interaction between sex and altitude. We indicated the correlation structure using the corAR () function to account for the spatial autocorrelation in the data and also used the AIC value to determine the best-fit model.
The models and statistical analyzes were performed in R (version 4.0.5, R Core Team, 2021) and Excel; the data visualization was conducted by using ggplot2, ggpubr, and ggpmisc; and the model results were checked by using the ANOVA function in the “car” package [65]. The R model code refers to the code in the related articles [25,26].

3. Results

The regression results of the FL and HTL on the BL for females and males of 51 Pteromalidae and 48 Eulophidae species confirmed the validity of the FL as the best surrogate for body size compared to the HTL (Figure 2, Figure S1). Meanwhile, the WL was considerably negatively correlated with the FL (Pteromalidae: cor = −0.639; Eulophidae: cor = −0.604), while there was no significant correlation between the WAR and the FL (Pteromalidae: cor = 0.024; Eulophidae: cor = −0.059). The FL increased significantly with altitude and decreased with temperature in both taxa for both sexes, and in both species and assemblage-based methods (Figure 3, Table 1).
The results of the species-based method revealed a more pronounced difference between sexes for both taxa; the mean FL of females was obviously larger than that of males, indicating that larger females dominated at all altitudinal ranges. Pteromalids are generally larger than Eulophids (Figure S2, all sexes combined; mean FL of Pteromalids = 1.51 mm, mean FL of Eulophids = 1.05 mm, t. test: p < 0.001); the FL size frequency of both taxa demonstrated significant right skewness (Figure S3), and the frequency peaks in both sexes transferred from left to right, indicating that the proportion of individuals with larger FL tends to be greater at a higher altitude, and this was remarkably evident in the males of both taxa. The slope values showed that the trend of FL variation with altitude was more pronounced in the assemblage-based method (Figure 3).
A stepwise best subset regression, including four environmental factors and sex, was conducted. The results of the GLM model on the FL from the two methods showed that it was negatively related to sex in both taxa (males have smaller FL, p < 0.01), and positively correlated with the annual mean precipitation (p < 0.01) and annual mean solar radiation, while negatively related to the annual mean temperature (p < 0.001) (Table 2). The annual mean wind speed was excluded and not included in any of the models with the lowest AIC. The temperature consistently provided the greatest independent explanation for the FL variation in both methods.
Five common subfamilies (Pteromalidae: Asaphinae, Miscogasterinae, Pteromalinae; Eulophidae: Entedoninae, and Tetrastichinae) showed an obvious trend (p < 0.05) of the FL increasing with altitude, except for Eulophinae (p > 0.05) (Figure S4), mainly due to the relatively large-sized Euplectrus, Cirrospilus, and Pnigalio generally being distributed at low altitudes (Figure 4). However, after removing the data of these genera, the FL of Eulophinae also demonstrated an increasing trend with altitude (p < 0.05). The altitudinal trend of Miscogasterinae and Entedoninae is particularly significant, and the variation in Halticoptera (p < 0.001) and Entedon (p < 0.001) made the greatest respective contributions (Figures S5 and S6). Based on the analysis of the genera with a sample size greater than 20 records, we found that most of the genera had a wide altitudinal range, between 200 and 4000 m; some alpine-dwelling genera, including Selderma, Callicarolynia, Stenomalina, Thinodytes, and Entedon, have a larger population at higher altitudes (≥3000 m), while some genera were only distributed at lower altitudes (≤1500 m), such as Norbanus, Philotrypesis, Dzhanokmenia, and Omphale.
We used the GLS model to validate Bergmann’s rule for 16 genera with widespread distribution and large populations (distributed in the altitudinal range of 200–4000 m and ≥ 30 records). The results indicated that the FL of 11 genera significantly increased with altitude among all 16 genera, while some genera demonstrated contrasting patterns (Table 3, Figures S5 and S6). The trend of Pteromalus was not significant (female: p = 0.244; male: p = 0.728), especially for its males, mostly due to their comparatively larger size in Pteromalidae, and the fact that there were also several larger species found at low altitudes. There was no significant correlation between the FL of Homoporus (female: p = 0.662; male: p = 0.707), Neochrysocharis (female: p = 0.411; male: p = 0.005), and Pediobius (female: p = 0.422; male: p = 0.281) and altitude in general. The FL of Chrysocharis (female: slope = −1.71; male: slope = −1.54) showed a reverse trend with increasing elevation. All the interactions between sex and altitude were not significant.

4. Discussion

Much of the previous research has documented the applicability of Bergmann’s rule in multiple insect groups at the species or family level, while no general pattern has been detected in parasitoids, even for the same group, and there has been little in-depth related research on minute parasitic wasps. In this study, we examined the geographical patterns in two hyper-diverse groups of Chalcidoidea and analyzed the dominant environmental factors responsible for this variation. We found that the FL of Pteromalidae and Eulophidae increased with elevation in both sexes at the family level (a pattern that follows Bergmann’s rule), with temperature and precipitation as the main environmental drivers.
Our results, based on the two different methods, might be explained by some of the hypotheses mentioned above. The average size of chalcids from these two taxa is larger at higher altitudes than at low altitudes. Temperature was obviously the most statistically explanatory variable in both taxa, explaining most of the variance in body size among other environmental factors in our study, and may play a more direct role. This phenomenon is also in line with numerous laboratory experiments, which show that Pnigalio soemius (Eulophidae: Eulophinae) raised at lower temperatures have larger FL [31], and are similar to Eretmocerus warrae (Hymenoptera: Aphelinidae). Larger insects usually have a smaller surface-area-to-volume ratio, which may help them maintain their thermal energy and reduce the water-loss rate [61], and also facilitate the storage of the energy needed to maintain normal activities [66]. Larger individuals are generally more fertile, which facilitates the reproduction of more offspring to maintain the population in the shorter reproductive periods in alpine regions. As stated in the starvation-resistance hypothesis, larger parasitoids are able to store more energy, such as lipids and glycogen, which are involved in egg production and other physiological mechanisms. The species’ diversity and coverage of flowering plants in the alpine area were significantly lower than those at lower altitudes, meaning that parasitoids need to restore more energy to cope with the extreme environment; thus, larger species are more suitable for survival in the alpine region. Precipitation also showed a high correlation with the FL variation; it is generally inversely proportional to elevation and is supposed to be negatively correlated with FL, as with temperature, while our results are contrary to the above speculation. Numerous investigations’ results about the effect of precipitation on body size are also controversial. Precipitation not only indirectly affects the growth of the host plants of parasitoids, and thus their population density, but also directly affects the water-loss rate of insects [67]. Therefore, precipitation may play an indirect role and interact with other environmental factors in complex ways. The geographical trend in size might be the result of interactions between various environmental factors. In addition to the factors considered above, the oxygen level (PO2) also declines significantly with altitude, and also has effects on the insect’s body size, combined with temperature [68]. The host identity and host plants also influence the parasitoid’s body size [69]. However, the various parasitoids analyzed in this study do not have specific host information, and the data on diet and DNA barcode also remain incomplete. This will also be a future research direction.
During the extensive field collection, we found brachypterous Pteromalidae and apterous Encyrtidae species in the alpine area (≥3500 m), while no similar species were found at low altitudes. We also found that the wing loading of Pteromalidae and Eulophidae was inversely proportional to the FL (Pteromalidae: cor = −0.639; Eulophidae: cor = −0.604) and decreased significantly (p < 0.001) with increasing altitude. A similar trend was also found in Drosophila [70], Bumbus [40], and moth [13], mainly because a lower wing loading is generally associated with a superior flight maneuverability, reduced energy being required for flight, and enhanced flight duration. Therefore, having this feature will help them live normally in the alpine environment. However, wing aspect ratio, also an important indicator of flight performance, did not significantly correlate with either body length or altitude in our research.
We found that the populations of Pahcyneuron, Halticoptera, and Selderma from the Pteromalidae, and Entedon, Neochrysocharis, and Diglyphus from the Eulophidae were highly dominant among other groups in the alpine region. Eleven out of sixteen genera followed Bergmann’s rule, and only the genera Pteromalus, Homoporus, Neochrysocharis, and Pediobius showed no significant trend, while Charysocharis showed a transverse trend. Pteromalus is one of the large-sized genera of Pteromalidae; the average FL of the Pteromalus specimens in our data is 1.782 mm(female), 1.35 mm(male), and the average FL of all the Pteromalidae is 1.652 mm(female), 1.422 mm(male). The larger individuals in this genus are also abundant at a lower altitude, resulting in an insignificant FL trend with elevation. Two female specimens of Chrysocharis pubicornis from Norway, found at elevations of 65 and 63 m, are larger than the average FL of the genus (female: 1.777 > 1.03 mm, male: 1.128 > 0.997 mm); the FL of the genus showed no significant trend with altitude when we removed these pieces of data. Such contradictory trends may also be caused by the interaction between various factors or there being less diversity in alpine species, and need further information on these genera to resolve this problem. Through the statistical analysis of numerous species, we also discovered some species with a wide distribution across altitudes, including Asaphes vulgaris, Pachyneuron korlense, Pachyneuron aphidis, Diglyphus isaea, Diaulinopsis arenaria, and Neochrysocharis formosus, etc., which will provide objects for further study of the adaptation mechanism.
In our research, the geographical trends in the size of two abundant chalcids were analyzed along a wide range of altitudinal gradients. The most seminal result is that the alpine region tends to harbor larger chalcids (following Bergmann’s rule) than the lower altitudes. Temperature and precipitation are the dominant environmental drivers of this variation, and the results indicate that temperature may be more directly related to body-size variation with altitudinal gradients than precipitation. Other hypotheses, such as resource availability, can also be used to explain the trend. The combination and interaction of these factors play an important role in the formation of body-size patterns. Additional information, such as genetic data (a test of whether body-size response to altitude shows a phylogenetic signal or there are positive selected genes) and more data on different species is needed to help further understand the mechanisms that cause body-size variations. Our results may provide a cornerstone for relevant future work.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects14010067/s1, Table S1: Overview of the numbers of analyzed specimens in Pteromalidae and Eulophidae. Figure S1: The correlation between seven morphological characteristics of both sexes in Pteromalidae(a) and Eulophidae(b); Figure S2: Frequency distribution histogram of the FL of Pteromalidae(a) and Eulophidae(b) from all elevation gradients; Figure S3: Histogram of the FL of Pteromalidae(a,b) and Eulophidae(c,d) from four elevation ranges; Figure S4: The relationship between forewing length and elevation in six subfamilies from Pteromalidae (a–c) and Eulophidae (d–f); Figure S5: The relationship between forewing length and elevation in eight genera(n ≥ 30 and with wide altitudinal range) of Pteromalidae; Figure S6: The relationship between forewing length and elevation in eight genera (n ≥ 30 and with wide altitudinal range) of Eulophidae. File S1: Forewing length data of Pteromalidae and Eulophidae which used for GLM and GLS model in R analysis.

Author Contributions

Conceptualization, N.K. and H.H.; methodology, N.K.; software, N.K.; validation, N.K.; formal analysis, N.K.; investigation, N.K., Z.H., S.L., and S.G.; resources, N.K., Z.H., S.L., and S.G.; data curation, N.K., Z.H., S.L., and S.G.; writing—original draft preparation, N.K.; writing—review and editing, N.K. and H.H.; visualization, N.K.; supervision, H.H.; project administration, H.H. and N.K.; funding acquisition, H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 31860612; Autonomous Region graduate scientific research innovation project, grant number XJ2019G025; and National Natural Science Foundation of China, grant number, 31672338.

Data Availability Statement

The data presented in this study are available on request from the corresponding author, and also available in supplementary material here.

Acknowledgments

We thank Bazhou Altun Mountain National Nature Reserve Administration for admission for scientific investigation. We sincerely thank the drivers Li Youyi and Bao gang for their favorable guidance, help, and concern during the scientific investigation in the Altun Mountain National Nature Reserve.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ollason, J.G.R.H. Peters 1986. The Ecological Implications of Body Size. Cambridge University Press, Cambridge. 329 Pages. ISBN 0-521-2886-x. J. Trop. Ecol. 1987, 3, 286–287. [Google Scholar] [CrossRef]
  2. Brown, J.H.; Gillooly, J.F.; Allen, A.P.; Savage, V.M.; West, G.B. Toward A Metabolic Theory of Ecology. Ecology 2004, 85, 1771–1789. [Google Scholar] [CrossRef]
  3. Chown, S.L.; Gaston, K.J. Body Size Variation in Insects: A Macroecological Perspective. Biol. Rev. 2010, 85, 139–169. [Google Scholar] [CrossRef] [PubMed]
  4. Bergmann, C. Über die Verhältnisse der Wärmeökonomie der Thiere zu Ihrer Grösse; Vandenhoeck und Ruprecht: Göttingen, Germany, 1847; Volume 1. [Google Scholar]
  5. Ashton, K.G.; Tracy, M.C.; de Queiroz, A. Is Bergmann’s Rule Valid for Mammals? Am. Nat. 2000, 156, 390–415. [Google Scholar] [CrossRef]
  6. Nwaogu, C.J.; Tieleman, B.I.; Bitrus, K.; Cresswell, W. Temperature and Aridity Determine Body Size Conformity to Bergmann’s Rule Independent of Latitudinal Differences in a Tropical Environment. J. Ornithol. 2018, 159, 1053–1062. [Google Scholar] [CrossRef] [Green Version]
  7. Blanckenhorn, W.U. Bergmann and Converse Bergmann Latitudinal Clines in Arthropods: Two Ends of a Continuum? Integr. Comp. Biol. 2004, 44, 413–424. [Google Scholar] [CrossRef] [Green Version]
  8. Salomão, R.P.; Arriaga-Jiménez, A.; Kohlmann, B. The Relationship between Altitudinal Gradients, Diversity, and Body Size in a Dung Beetle (Coleoptera: Scarabaeinae: Onthophagus) Model System. Can. J. Zool. 2021, 99, 33–43. [Google Scholar] [CrossRef]
  9. Bidau, C.J.; Martí, D.A. Clinal Variation of Body Size in Dichroplus Pratensis (Orthoptera: Acrididae): Inversion of Bergmann’s and Rensch’s Rules. Ann. Entomol. Soc. Am. 2007, 100, 850–860. [Google Scholar] [CrossRef]
  10. Schutze, M.K.; Clarke, A.R. Converse Bergmann Cline in a Eucalyptus Herbivore, Paropsis Atomaria Olivier (Coleoptera: Chrysomelidae): Phenotypic Plasticity or Local Adaptation? Glob. Ecol. Biogeogr. 2008, 17, 424–431. [Google Scholar] [CrossRef] [Green Version]
  11. Bhusal, D.R.; Patel, P.; Ghimire, K.C.; Bista, M.; Kumar, B. Size-Based Intraspecific Variations along Altitudinal Gradient and Follow of Bergmann’s Rule in Parthenium Beetle, Zygogramma Bicolorata Pallister. J. Asia-Pac. Entomol. 2019, 22, 1173–1179. [Google Scholar] [CrossRef]
  12. Chown, S.L.; Klok, C.J. Altitudinal Body Size Clines: Latitudinal Effects Associated with Changing Seasonality. Ecography 2003, 26, 445–455. [Google Scholar] [CrossRef]
  13. Brehm, G.; Zeuss, D.; Colwell, R.K. Moth Body Size Increases with Elevation along a Complete Tropical Elevational Gradient for Two Hyperdiverse Clades. Ecography 2019, 42, 632–642. [Google Scholar] [CrossRef] [Green Version]
  14. Xing, S.; Cheng, W.; Nakamura, A.; Tang, C.C.; Huang, S.; Odell, E.; Goodale, E.; Goodale, U.M.; Bonebrake, T.C. Elevational Clines in Morphological Traits of Subtropical and Tropical Butterfly Assemblages. Biol. J. Linn. Soc. 2018, 123, 506–517. [Google Scholar] [CrossRef] [Green Version]
  15. McCulloch, G.A.; Waters, J.M. Does Wing Reduction Influence the Relationship between Altitude and Insect Body Size? A Case Study Using New Zealand’s Diverse Stonefly Fauna. Ecol. Evol. 2018, 8, 953–960. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Pérez-Valencia, L.I.; Moya-Raygoza, G. Body Size Variation of Diaphorina citri (Hemiptera: Psyllidae) through an Elevation Gradient. Ann. Entomol. Soc. Am. 2015, 108, 800–806. [Google Scholar] [CrossRef]
  17. Beerli, N.; Bärtschi, F.; Ballesteros-Mejia, L.; Kitching, I.J.; Beck, J. How Has the Environment Shaped Geographical Patterns of Insect Body Sizes? A Test of Hypotheses Using Sphingid Moths. J. Biogeogr. 2019, 46, 1687–1698. [Google Scholar] [CrossRef]
  18. Horne, C.R.; Hirst, A.G.; Atkinson, D. Insect Temperature–Body Size Trends Common to Laboratory, Latitudinal and Seasonal Gradients Are Not Found across Altitudes. Funct. Ecol. 2018, 32, 948–957. [Google Scholar] [CrossRef] [Green Version]
  19. Sandoval-Arango, S.; Cárdenas Henao, H.; Montoya-Lerma, J. Divergence in Bergmann’s Clines: Elevational Variation and Heritability of Body Size in a Leaf-Cutting Ant. Insectes Soc. 2020, 67, 355–366. [Google Scholar] [CrossRef]
  20. Traynor, R.E.; Mayhew, P.J. A Comparative Study of Body Size and Clutch Size across the Parasitoid Hymenoptera. Oikos 2005, 109, 305–316. [Google Scholar] [CrossRef]
  21. Sveum, P.; Solem, J.O. Anaphes cultripennis Debauché, 1948 (Hymenoptera, Mymaridae) in Norway. Fauna Nor. 1980, 27, 17–18. [Google Scholar]
  22. Jervis, M.A.; Ferns, P.N.; Heimpel, G.E. Body Size and the Timing of Egg Production in Parasitoid Wasps: A Comparative Analysis. Funct. Ecol. 2003, 17, 375–383. [Google Scholar] [CrossRef]
  23. Sagarra, L.A.; Vincent, C.; Stewart, R.K. Body Size as an Indicator of Parasitoid Quality in Male and Female Anagyrus kamali (Hymenoptera: Encyrtidae). Bull. Entomol. Res. 2001, 91, 363–367. [Google Scholar] [CrossRef] [PubMed]
  24. Hodkinson, I.D. Terrestrial Insects along Elevation Gradients: Species and Community Responses to Altitude. Biol. Rev. 2005, 80, 489. [Google Scholar] [CrossRef] [Green Version]
  25. García-Navas, V.; Noguerales, V.; Cordero, P.J.; Ortego, J. Ecological Drivers of Body Size Evolution and Sexual Size Dimorphism in Short-Horned Grasshoppers (Orthoptera: Acrididae). J. Evol. Biol. 2017, 30, 1592–1608. [Google Scholar] [CrossRef] [Green Version]
  26. Baranov, V.; Jourdan, J.; Hunter-Moffatt, B.; Noori, S.; Schölderle, S.; Haug, J.T. Global Size Pattern in a Group of Important Ecological Indicators (Diptera, Chironomidae) Is Driven by Latitudinal Temperature Gradients. Insects 2021, 13, 34. [Google Scholar] [CrossRef]
  27. Wonglersak, R.; Fenberg, P.B.; Langdon, P.G.; Brooks, S.J.; Price, B.W. Temperature-body Size Responses in Insects: A Case Study of British Odonata. Ecol. Entomol. 2020, 45, 795–805. [Google Scholar] [CrossRef]
  28. Acquah-Lamptey, D.; Brändle, M.; Brandl, R.; Pinkert, S. Temperature-driven Color Lightness and Body Size Variation Scale to Local Assemblages of European Odonata but Are Modified by Propensity for Dispersal. Ecol. Evol. 2020, 10, 8936–8948. [Google Scholar] [CrossRef]
  29. Baroni, D.; Masoero, G. Complex Influence of Climate on the Distribution and Body Size of an Alpine Species. Insect Conserv. Divers. 2018, 11, 435–448. [Google Scholar] [CrossRef]
  30. Colinet, H.; Boivin, G.; Hance, T. Manipulation of Parasitoid Size Using the Temperature-Size Rule: Fitness Consequences. Oecologia 2007, 152, 425–433. [Google Scholar] [CrossRef]
  31. Bernardo, U.; Pedata, P.A.; Viggiani, G. Phenotypic Plasticity of Pigmentation and Morphometric Traits in Pnigalio soemius (Hymenoptera: Eulophidae). Bull. Entomol. Res. 2007, 97, 101–109. [Google Scholar] [CrossRef]
  32. Wang, T.; Keller, M.A. Larger Is Better in the Parasitoid Eretmocerus Warrae (Hymenoptera: Aphelinidae). Insects 2020, 11, 39. [Google Scholar] [CrossRef]
  33. Montgomery, K. Variation in Temperature With Altitude and Latitude. J. Geogr. 2006, 105, 133–135. [Google Scholar] [CrossRef]
  34. Angilletta, M.J., Jr.; Dunham, A.E. The Temperature-Size Rule in Ectotherms: Simple Evolutionary Explanations May Not Be General. Am. Nat. 2003, 162, 332–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Atkinson, D. Temperature and Organism Size—A Biological Law for Ectotherms? In Advances in Ecological Research; Elsevier: San Leandro, CA, USA, 1994; Volume 25, pp. 1–58. ISBN 978-0-12-013925-5. [Google Scholar]
  36. Zhang, D.-W.; Xiao, Z.-J.; Zeng, B.-P.; Li, K.; Tang, Y.-L. Insect Behavior and Physiological Adaptation Mechanisms Under Starvation Stress. Front. Physiol. 2019, 10, 163. [Google Scholar] [CrossRef]
  37. Gergs, A.; Jager, T. Body Size-Mediated Starvation Resistance in an Insect Predator. J. Anim. Ecol. 2014, 83, 758–768. [Google Scholar] [CrossRef] [PubMed]
  38. Ramírez-Delgado, V.H.; Sanabria-Urbán, S.; Serrano-Meneses, M.A.; Cueva del Castillo, R. The Converse to Bergmann’s Rule in Bumblebees, a Phylogenetic Approach. Ecol. Evol. 2016, 6, 6160–6169. [Google Scholar] [CrossRef] [PubMed]
  39. Scriven, J.J.; Whitehorn, P.R.; Goulson, D.; Tinsley, M.C. Bergmann’s Body Size Rule Operates in Facultatively Endothermic Insects: Evidence from a Complex of Cryptic Bumblebee Species. PLoS ONE 2016, 11, e0163307. [Google Scholar] [CrossRef]
  40. Lozier, J.D.; Parsons, Z.M.; Rachoki, L.; Jackson, J.M.; Pimsler, M.L.; Oyen, K.J.; Strange, J.; Dillon, M.E. Divergence in Body Mass, Wing Loading, and Population Structure Reveals Species-Specific and Potentially Adaptive Trait Variation Across Elevations in Montane Bumble Bees. Insect Syst. Divers. 2021, 5, 3. [Google Scholar] [CrossRef]
  41. Osorio-Canadas, S.; Flores-Hernández, N.; Sánchez-Ortiz, T.; Valiente-Banuet, A. Changes in Bee Functional Traits at Community and Intraspecific Levels along an Elevational Gradient in a Mexical-Type Scrubland. Oecologia 2022, 200, 145–158. [Google Scholar] [CrossRef]
  42. Gérard, M.; Vanderplanck, M.; Franzen, M.; Kuhlmann, M.; Potts, S.G.; Rasmont, P.; Schweiger, O.; Michez, D. Patterns of Size Variation in Bees at a Continental Scale: Does Bergmann’s Rule Apply? Oikos 2018, 127, 1095–1103. [Google Scholar] [CrossRef] [Green Version]
  43. Osorio-Canadas, S.; Arnan, X.; Rodrigo, A.; Torné-Noguera, A.; Molowny, R.; Bosch, J. Body Size Phenology in a Regional Bee Fauna: A Temporal Extension of Bergmann’s Rule. Ecol. Lett. 2016, 19, 1395–1402. [Google Scholar] [CrossRef]
  44. Santos, A.M.C.; Quicke, D.L.J. Large-Scale Diversity Patterns of Parasitoid Insects: Parasitoid Diversity Patterns. Entomol. Sci. 2011, 14, 371–382. [Google Scholar] [CrossRef]
  45. Eagalle, T.I. Latitudinal Effects on Diversity and Body Size: A Case Study with Parasitoid and Parasitic Wasps. Master’s Thesis, The University of Guelph, Guelph, ON, Canada, 2014. [Google Scholar]
  46. Geraghty, M.J.; Dunn, R.R.; Sanders, N.J. Body Size, Colony Size, and Range Size in Ants (Hymenoptera: Formicidae): Are Patterns along Elevational and Latitudinal Gradients Consistent with Bergmann’s Rule? Myrmecol. News 2007, 10, 51–58. [Google Scholar]
  47. Bernadou, A.; Römermann, C.; Gratiashvili, N.; Heinze, J. Body Size but Not Colony Size Increases with Altitude in the Holarctic Ant, Leptothorax acervorum: Bergmann’s Rule in a Holarctic Ant. Ecol. Entomol. 2016, 41, 733–736. [Google Scholar] [CrossRef]
  48. Nowrouzi, S.; Andersen, A.N.; Bishop, T.R.; Robson, S.K.A. Is Thermal Limitation the Primary Driver of Elevational Distributions? Not for Montane Rainforest Ants in the Australian Wet Tropics. Oecologia 2018, 188, 333–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Heinze, J.; Foitzik, S.; Fischer, B.; Wanke, T.; Kipyatkov, V.E. The Significance of Latitudinal Variation in Body Size in a Holarctic Ant, Leptothorax acervorum. Ecography 2003, 26, 349–355. [Google Scholar] [CrossRef] [Green Version]
  50. Miller, S.E.; Sheehan, M.J. Ecogeographical Patterns of Body Size Differ among North American Paper Wasp Species. Insectes Soc. 2021, 68, 109–122. [Google Scholar] [CrossRef]
  51. Rodríguez-Jimenez, A.; Sarmiento, C.E. Altitudinal Distribution and Body Resource Allocation in a High Mountain Social Wasp (Hymenoptera: Vespidae). Neotrop. Entomol. 2008, 37, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Gibson, G. Annotated Keys to the Genera of Nearctic Chalcidoidea (Hymenoptera); NRC Research Press: Ottawa, ON, Canada, 1997; ISBN 0-660-16669-0. [Google Scholar]
  53. Bouček, Z. Australasian Chalcidoidea (Hymenoptera). A Biosystematic Revision of Genera of Fourteen Families, with a Reclassification of Species; The Cambrian News Ltd.: Aberystwyth, Wales; CAB International: Wallingford, UK, 1988. [Google Scholar]
  54. Huang, D.-W.; Xiao, H. Fauna Sinica, Insecta, Volume 42, Hymenoptera: Pteromalidae; Science Press: Beijing, China, 2005. [Google Scholar]
  55. Bouček, Z.; Rasplus, J.Y. Ilustrated Key to West-Palearctic Genera of Pteromalidae; Institut National de la Recherche Agronomique: Paris, France, 1991; ISBN 2-7380-0343-5. [Google Scholar]
  56. De Vere Graham, M.W.R. The Pteromalidae of North-Western Europe (Hymenoptera: Chalcidoidea). Bull. Br. Mus. (Nat. Hist.) Entomol. Suppl. 1969, 16, 1–908. [Google Scholar] [CrossRef]
  57. Montejo-Kovacevich, G.; Salazar, P.A.; Smith, S.H.; Gavilanes, K.; Bacquet, C.N.; Chan, Y.F.; Jiggins, C.D.; Meier, J.I.; Nadeau, N.J. Genomics of Altitude-associated Wing Shape in Two Tropical Butterflies. Mol. Ecol. 2021, 30, 6387–6402. [Google Scholar] [CrossRef]
  58. Starmer, W.T.; Wolf, L.L. Causes of Variation in Wing Loading among Drosophila Species. Biol. J. Linn. Soc. 1989, 37, 247–261. [Google Scholar] [CrossRef]
  59. Bhat, S.S.; Sheridan, J.; Hourigan, K.; Thompson, M.C. Aspect Ratio Studies on Insect Wings. Phys. Fluids 2019, 31, 121301. [Google Scholar] [CrossRef]
  60. Sarig, A.; Ribak, G. To What Extent Can the Tiny Parasitoid Wasps, Eretmocerus Mundus, Fly Upwind? J. Appl. Entomol. 2021, 145, 660–674. [Google Scholar] [CrossRef]
  61. Kühsel, S.; Brückner, A.; Schmelzle, S.; Heethoff, M.; Blüthgen, N. Surface Area-Volume Ratios in Insects: Surface Area Allometry and Water Loss. Insect Sci. 2017, 24, 829–841. [Google Scholar] [CrossRef] [PubMed]
  62. Remmert, H. Body Size of Terrestrial Arthropods and Biomass of Their Populations in Relation to the Abiotic Parameters of Their Milieu. Oecologia 1981, 50, 12–13. [Google Scholar] [CrossRef]
  63. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km Spatial Resolution Climate Surfaces for Global Land Areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  64. ESRI. ArcGIS Desktop: Release 10.2.2; Environmental Systems Research Institute: Redlands, CA, USA, 2014. [Google Scholar]
  65. Fox, J.; Weisberg, S. An R Companion to Applied Regression, 3rd ed.; Sage: Thousand Oaks, CA, USA, 2018. [Google Scholar]
  66. Rivero, A.; West, S.A. The Physiological Costs of Being Small in a Parasitic Wasp. Evol. Ecol. Res. 2002, 4, 407–420. [Google Scholar]
  67. Addo-Bediako, A.; Chown, S.L.; Gaston, K.J. Revisiting Water Loss in Insects: A Large Scale View. J. Insect Physiol. 2001, 47, 1377–1388. [Google Scholar] [CrossRef]
  68. Dillon, M.E. Into Thin Air: Physiology and Evolution of Alpine Insects. Integr. Comp. Biol. 2006, 46, 49–61. [Google Scholar] [CrossRef] [Green Version]
  69. Symonds, M.R.E.; Elgar, M.A. The Evolution of Body Size, Antennal Size and Host Use in Parasitoid Wasps (Hymenoptera: Chalcidoidea): A Phylogenetic Comparative Analysis. PLoS ONE 2013, 8, e78297. [Google Scholar] [CrossRef]
  70. Gilchrist, G.W. Plastic and Genetic Variation in Wing Loading as a Function of Temperature Within and Among Parallel Clines in Drosophila Subobscura. Integr. Comp. Biol. 2004, 44, 461–470. [Google Scholar] [CrossRef] [PubMed]
Insects 14 00067 g001 550
Figure 1. (a) Geographic distribution of Pteromalidae and Eulophidae specimens used in this analysis. (b) Morphometric measurement illustrations. Upper: microslide specimen of male forewing of Sphegigaster stepicola and the scheme of wing length and wing area measurement. Lower: body of female Pachyneuron korlense and the scheme of body length, mesothorax length, tibia length, wing length, and wing area measurement.
Figure 1. (a) Geographic distribution of Pteromalidae and Eulophidae specimens used in this analysis. (b) Morphometric measurement illustrations. Upper: microslide specimen of male forewing of Sphegigaster stepicola and the scheme of wing length and wing area measurement. Lower: body of female Pachyneuron korlense and the scheme of body length, mesothorax length, tibia length, wing length, and wing area measurement.
Insects 14 00067 g001
Insects 14 00067 g002 550
Figure 2. Forewing length (FL) and hind tibia length (HTL) of Pteromalidae (n = 51 males and females of the same species) and Eulophidae (n = 48 males and females of the same species) in relation to body length (BL). (a) Correlation between FL and BL of Pteromalidae. (b) Correlation between FL and BL of Eulophidae. (c) Correlation between HTL and BL of Pteromalidae. (d) Correlation between HTL and BL of Eulophidae. All p < 0.001. Green circles: male; red circles: female. For detailed comparison of morphological characteristics, see Supplementary Material, Figure S1.
Figure 2. Forewing length (FL) and hind tibia length (HTL) of Pteromalidae (n = 51 males and females of the same species) and Eulophidae (n = 48 males and females of the same species) in relation to body length (BL). (a) Correlation between FL and BL of Pteromalidae. (b) Correlation between FL and BL of Eulophidae. (c) Correlation between HTL and BL of Pteromalidae. (d) Correlation between HTL and BL of Eulophidae. All p < 0.001. Green circles: male; red circles: female. For detailed comparison of morphological characteristics, see Supplementary Material, Figure S1.
Insects 14 00067 g002
Insects 14 00067 g003 550
Figure 3. Forewing length of Pteromalidae and Eulophidae in relation to altitude based on two methods: (a) species-based method, (b) assemblage-based method. Green circles: male; red circles: female. See Supplementary Material, Figures S2 and S3 for detailed frequency histograms of four different elevation ranges.
Figure 3. Forewing length of Pteromalidae and Eulophidae in relation to altitude based on two methods: (a) species-based method, (b) assemblage-based method. Green circles: male; red circles: female. See Supplementary Material, Figures S2 and S3 for detailed frequency histograms of four different elevation ranges.
Insects 14 00067 g003
Insects 14 00067 g004 550
Figure 4. Altitudinal distribution of Pteromalidae (a) and Eulophidae (b), and frequency of occurrence records by altitude. Box and whisker plots show the median (cross symbol), the 25th and 75th percentile (bottom and top of the box), the minimum and maximum values (lower and upper whiskers), and outliers (circles). The red line shows the frequency distribution of altitudinal values in the assemblage-based data.
Figure 4. Altitudinal distribution of Pteromalidae (a) and Eulophidae (b), and frequency of occurrence records by altitude. Box and whisker plots show the median (cross symbol), the 25th and 75th percentile (bottom and top of the box), the minimum and maximum values (lower and upper whiskers), and outliers (circles). The red line shows the frequency distribution of altitudinal values in the assemblage-based data.
Insects 14 00067 g004
Table
Table 2. The generalized linear model (GLM) regression of relationships between FL of Pteromalidae and Eulophidae and environmental variables. AMTEM: annual mean temperature (°C), AMPRE: annual mean precipitation (mm), AMSR: annual mean solar radiation (kJ/m2), AMWS: annual mean wind speed (m/s). (** </ =0.01; *** </ =0.001).
Table 2. The generalized linear model (GLM) regression of relationships between FL of Pteromalidae and Eulophidae and environmental variables. AMTEM: annual mean temperature (°C), AMPRE: annual mean precipitation (mm), AMSR: annual mean solar radiation (kJ/m2), AMWS: annual mean wind speed (m/s). (** </ =0.01; *** </ =0.001).
Species-Based MethodRegression CoefficientsSE95% Confidence Intervalp-Value
Pteromalidae
SexM−0.14310.0334−0.208, 0.07782.83 × 10−5 ***
AMTEM−0.02053.417 × 10−3−0.0291, −0.01583.59 × 10−8 ***
APRE1.449 × 10−46.262 × 10−51.253 × 10−4, 3.675 × 10−40.000116 ***
AMSR3.051 × 10−58.025 × 10−6−1.046 × 10−6, 3.003 × 10−50.0722
Eulophidae
SexM−0.1180.0321−0.018, −0.050.000262 **
AMTEM−0.017562.977 × 10−30.023, −0.01787.3 × 10−9 ***
APRE1.305 × 10−43.121 × 10−56.894 × 10−5, 1.921 × 10−43.5 × 10−5 ***
AMSR1.138 × 10−56.059 × 10−6−4.41 × 10−7, 2.32 × 10−50.061
Assemblage-based method
Pteromalidae
SexM−0.18860.0273−0.242, −0.1352.84 × 10−11 ***
AMTEM−0.02052.636 × 10−3−0.0256, −0.01541.33 × 10−13 ***
APRE1.499 × 10−42.773 × 10−58.945 × 10−5, 0.00023.27 × 10−7 ***
AMSR3.051 × 10−56.503 × 10−61.807 × 10−5, 4.298 × 10−54.13 × 10−6 ***
Eulophidae
SexM−0.15260.02768−0.2068, −0.09837.47 × 10−8 ***
AMTEM−0.01492.595 × 10−3−0.0200, −9.792 × 10−37.36 × 10−6 ***
APRE1.267 × 10−42.777 × 10−57.079 × 10−5, 1.8309 × 10−42.38 × 10−8 ***
AMSR1.136 × 10−56.031 × 10−6−4.274 × 10−7, 2.317 × 10−50.0607
Table
Table 3. Coefficient and ANOVA of the best fit GLS models on FL of 16 genera from Pteromalidae and Eulophidae (CE = coefficients, n = number of specimens). * denotes level of significance (* </ =0.05; ** </ =0.01; *** </ =0.001).
Table 3. Coefficient and ANOVA of the best fit GLS models on FL of 16 genera from Pteromalidae and Eulophidae (CE = coefficients, n = number of specimens). * denotes level of significance (* </ =0.05; ** </ =0.01; *** </ =0.001).
CEχ2p-Value CEχ2p-Value
Pteromalidae Eulophidae
Halticoptera (n = 221) Diglyphus (n = 98)
Altitude2.41 × 10−417.9192<0.0001 ***Altitude8.5 × 10−520.6508<0.0001 ***
SexM0.011.59680.2077SexM−0.0261.11480.2937
Altitude × Sex−1.5 × 10−62.37740.1245Altitude × Sex−1.27 × 10−50.10440.7474
Pachyneuron (n = 162) Diaulinopsis (n = 80)
Altitude4.28 × 10−519.1234<0.0001 ***Altitude6.85 × 10−553.949<0.0001 ***
SexM−3.9616.30630.0001 ***SexM−0.0130.326<0.0001 ***
Altitude × Sex1.08 × 10−32.06790.1524Altitude × Sex−4 × 10−510.1340.0021 **
Sphegigaster (n = 71) Hyssopus (n = 30)
Altitude1.81 × 10−516.20640.0001 ***Altitude8.16 × 10−539.0975<0.0001 ***
SexM−0.5816.17330.0001 ***SexM0.2271.08530.3071
Altitude × Sex1.45 × 10−43.71990.058Altitude × Sex5.43 × 10−51.53370.2266
Asaphes (n = 64) Chrysocharis
(n = 107)
Altitude2.13 × 10−570.0104 *Altitude−1.7 × 10−61.3540.2473
SexM−0.35827.5398<0.0001 ***SexM−0.0612.4330.1219
Altitude × Sex8.48 × 10−58.80150.0043 **Altitude × Sex1.6 × 10−60.0020.9641
Pteromalus (n = 145) Entedon (n = 150)
Altitude4.84 × 10−50.78180.3777Altitude2.02 × 10−4108.014<0.0001 ***
SexM−3.88536.9243<0.0001 ***SexM−0.0152.6820.1037
Altitude × Sex−4.9 × 10−51.41120.2364Altitude × Sex−1.39 × 10−50.0890.7656
Homoporus (n = 30) Neochrysocharis (n = 115)
Altitude7.69 × 10−53.6230.0681Altitude2.7 × 10−61.07250.3026
SexM−0.657873.692<0.0001 ***SexM−1.627.74980.0063 **
Altitude × Sex−1.24 × 10−41.8780.1823Altitude × Sex3.62 × 10−53.40630.0676
Selderma (n = 33) Pediobius (n = 30)
Altitude3.5 × 10−431.4179<0.0001 ***Altitude1.18 × 10−40.27450.6047
SexM0.11877.27940.0115 *SexM−0.3977.17930.0126 *
Altitude × Sex−1.15 × 10−41.4970.231Altitude × Sex8.37 × 10−50.12100.7308
Sphaeripalpus (n = 34) Aprostocetus (n = 356)
Altitude3.5 × 10−412.56130.0012 **Altitude6.23 × 10−520.992<0.0001 ***
SexM0.00481.82550.1859SexM−0.11719.163<0.0001 ***
Altitude × Sex−2.82 × 10−50.0670.7973Altitude × Sex−2.13 × 10−50.7390.3904
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

Kang, N.; Hu, H.; Huang, Z.; Luo, S.; Guo, S. Environmental Factors Drive Chalcid Body Size Increases with Altitudinal Gradients for Two Hyper-Diverse Taxa. Insects 2023, 14, 67. https://doi.org/10.3390/insects14010067

AMA Style

Kang N, Hu H, Huang Z, Luo S, Guo S. Environmental Factors Drive Chalcid Body Size Increases with Altitudinal Gradients for Two Hyper-Diverse Taxa. Insects. 2023; 14(1):67. https://doi.org/10.3390/insects14010067

Chicago/Turabian Style

Kang, Ning, Hongying Hu, Zengqian Huang, Shungang Luo, and Shuhan Guo. 2023. "Environmental Factors Drive Chalcid Body Size Increases with Altitudinal Gradients for Two Hyper-Diverse Taxa" Insects 14, no. 1: 67. https://doi.org/10.3390/insects14010067

APA Style

Kang, N., Hu, H., Huang, Z., Luo, S., & Guo, S. (2023). Environmental Factors Drive Chalcid Body Size Increases with Altitudinal Gradients for Two Hyper-Diverse Taxa. Insects, 14(1), 67. https://doi.org/10.3390/insects14010067

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