Skip to main content
eLife logoLink to eLife
. 2022 Dec 6;11:e80859. doi: 10.7554/eLife.80859

Desiccation resistance differences in Drosophila species can be largely explained by variations in cuticular hydrocarbons

Zinan Wang 1,2,, Joseph P Receveur 1,2,3, Jian Pu 1,4, Haosu Cong 1, Cole Richards 1, Muxuan Liang 5, Henry Chung 1,2,
Editors: Lynn M Riddiford6, K VijayRaghavan7
PMCID: PMC9757832  PMID: 36473178

Abstract

Maintaining water balance is a universal challenge for organisms living in terrestrial environments, especially for insects, which have essential roles in our ecosystem. Although the high surface area to volume ratio in insects makes them vulnerable to water loss, insects have evolved different levels of desiccation resistance to adapt to diverse environments. To withstand desiccation, insects use a lipid layer called cuticular hydrocarbons (CHCs) to reduce water evaporation from the body surface. It has long been hypothesized that the water-proofing capability of this CHC layer, which can confer different levels of desiccation resistance, depends on its chemical composition. However, it is unknown which CHC components are important contributors to desiccation resistance and how these components can determine differences in desiccation resistance. In this study, we used machine-learning algorithms, correlation analyses, and synthetic CHCs to investigate how different CHC components affect desiccation resistance in 50 Drosophila and related species. We showed that desiccation resistance differences across these species can be largely explained by variation in CHC composition. In particular, length variation in a subset of CHCs, the methyl-branched CHCs (mbCHCs), is a key determinant of desiccation resistance. There is also a significant correlation between the evolution of longer mbCHCs and higher desiccation resistance in these species. Given that CHCs are almost ubiquitous in insects, we suggest that evolutionary changes in insect CHC components can be a general mechanism for the evolution of desiccation resistance and adaptation to diverse and changing environments.

Research organism: D. melanogaster, None

Introduction

Maintaining water balance is a universal challenge for organisms living in terrestrial environments where water constantly evaporates from the body surface (Hadley, 1994). Organisms such as insects have a high surface area to volume ratio due to their small sizes, rendering them vulnerable to water loss and desiccation (Gibbs, 2002b; Gibbs and Rajpurohit, 2010; Kühsel et al., 2017). Studies using closely related insect species collected worldwide showed that phylogenetically related insect species can have very different levels of desiccation resistance and occupy very distinct habitats, while insects that are phylogenetically distant but dwelling in similar habitats could have similar levels of desiccation resistance (Kellermann et al., 2012; Li et al., 2022; Menzel et al., 2017; Rane et al., 2019). This suggests that extant species have evolved different levels of desiccation resistance to survive in their different habitats. However, as climate change accelerates the expansion of dryland (Huang et al., 2016) and changes aridity in many areas across the globe (Sherwood and Fu, 2014; Shi et al., 2021), it is less understood how insect species, which are integral parts of our ecosystems, can evolve higher levels of desiccation resistance to adapt to the more arid environments.

To determine how insects can adapt to desiccation stress, understanding how insects conserve water and prevent desiccation is crucial. In insects, cuticular water loss is the leading cause of desiccation (Chown et al., 2011; Gibbs and Matzkin, 2001; Wang et al., 2021). In an extreme example, cuticular water transpiration has been found to account for 97% of increased water loss in queens of the harvester ant Pogonomyrmex barbatus (Johnson and Gibbs, 2004). To conserve water and prevent desiccation, a general mechanism in insects is the use of a lipid layer on the epicuticle, named cuticular hydrocarbons (CHCs) (Gibbs et al., 2003; Gibbs and Matzkin, 2001). This CHC layer, which can contain more than 100 different compounds on the same individual, provides a hydrophobic barrier against evaporative water loss through the cuticle (Blomquist and Ginzel, 2021). This was first demonstrated when the physical or chemical removal of this layer using abrasive dust and various solvents resulted in increased water loss in various insect species in experiments published almost eight decades ago (Wigglesworth, 1945). In recent years, after the identification of a cytochrome P450 decarbonylase responsible for the synthesis of insect CHCs (Qiu et al., 2012), genetic manipulations of this gene in the fruitfly Drosophila melanogaster as well as in aphids and cockroaches, led to almost complete loss of CHCs and significant decrease in desiccation resistance (Chen et al., 2016; Chen et al., 2019; Qiu et al., 2012).

Variations in the composition of this CHC layer have been suggested to contribute to intraspecific variation in desiccation resistance in D. melanogaster (Rouault et al., 2004), the Mediterranean dung beetle Onthophagus taurus (Leeson et al., 2020), and the Argentine ant Linepithema humile (Buellesbach et al., 2018). This is supported by studies manipulating CHC composition either chemically or genetically, in different Drosophila species, producing different levels of desiccation resistance (Chiang et al., 2016; Ferveur et al., 2018; Chung et al., 2014; Koto et al., 2019; Savage et al., 2021). Together, these studies suggest that different components in the CHC layer can influence desiccation resistance.

The varying ability of these CHC components to prevent desiccation depends on its chemical structure, which in turn determines its melting temperatures (Gibbs and Pomonis, 1995; Gibbs, 1998; Gibbs, 2002a). The melting temperature of the hydrocarbon is positively correlated with its water-proofing properties and contribution to desiccation resistance (Blomquist and Ginzel, 2021; Gibbs and Rajpurohit, 2010). In Drosophila and most other insects, CHCs range in lengths between approximately 21 and 50 carbons, and consist of linear alkanes (n-alkanes), alkenes (monoenes and dienes), and methyl-branched alkanes (mbCHCs) (Blomquist and Ginzel, 2021). Among them, n-alkanes have the highest melting temperature, followed by mbCHCs , monoenes, and dienes. Increase in CHC chain length can also increase melting temperature and potentially leads to higher desiccation resistance (Gibbs and Pomonis, 1995). This is consistent with laboratory selection experiments in D. melanogaster selecting for increased desiccation resistance resulted in longer carbon-chain CHCs in the desiccation selected flies than the control flies (Gibbs et al., 1997; Kwan and Rundle, 2010). Longer CHC chain length is also correlated with climatic factors such as higher temperature and lower precipitation (Jezovit et al., 2017; Rouault et al., 2004). These climatic factors are associated with desiccation resistance in Drosophila species (Hoffmann, 2010; Hoffmann and Weeks, 2007; Kellermann et al., 2018).

While the studies mentioned above showed strong evidence that CHCs and variations in CHC compositions are associated with desiccation resistance in insects, no study has investigated the extend that CHC variation can determine desiccation resistance and whether we can identify the important CHC components that underlie variation in desiccation resistance across different species. Several large-scale studies produced large datasets in determining the evolution of desiccation resistance across closely related Drosophila species (Kellermann et al., 2018; Kellermann et al., 2012; Matzkin et al., 2009) and ant species (Bujan et al., 2016; Hood and Tschinkel, 1990). Other studies produced large datasets in CHC variation across species (Khallaf et al., 2021; Menzel et al., 2017; Nunes et al., 2017; Van Oystaeyen et al., 2014) and focused on the communication aspects of CHCs. However, variations in CHCs and differences in desiccation resistance have not been experimentally connected in a phylogenetic framework. Analyzing CHC profiles and desiccation resistance across closely related species measured under the similar conditions can determine important CHC components that could predict various levels of desiccation resistance. Understanding the mechanisms underlying desiccation resistance and examining the evolution of these traits is key to predicting species’ responses in facing drier environments that could result from future climate change (Chown et al., 2010).

In this study, we used a cohort-based experimental design of 50 Drosophila and related species, and experimentally determined their CHC compositions and desiccation resistance under similar conditions. Using a random forest machine-learning algorithm, we built decision trees connecting CHC components and desiccation resistance in the 50 species and tested their correlation. First, we determined that CHC variations can largely explain differences in desiccation resistance across these species. Second, we identified a subset of CHCs, the methyl-branched CHCs (mbCHCs), as the most important CHC components that can determine desiccation resistance across these species. Third, we determined that the evolution of longer mbCHCs is also significantly correlated with evolution of higher desiccation resistance in Drosophila species. Taken together, our study showed that mbCHCs are a key predictive component of desiccation resistance in Drosophila species and suggests that evolution in CHC components can be a general mechanism for higher desiccation resistance and adaptation to more arid environments.

Results

A cohort-based experimental design for determining correlations between CHCs and desiccation resistance

To investigate how CHC variation across an evolutionary gradient affects desiccation resistance, for our experiments, we selected 46 Drosophila species representing both the Sophophora subgenus and the Drosophila subgenus, as well as three Scaptodrosophila species and one Chymomyza species (Figure 1—figure supplement 1). As the assays for CHCs and desiccation resistance cannot be performed on the same individuals, we used a cohort-based experimental design (Figure 1) where subsets of each cohort (5–6 per species) were used for gas chromatography–mass spectrometry (GC–MS) determination of CHC profiles, desiccation resistance assays, and measuring body weight of the F1 progeny for each sex. Therefore, all measurements were performed at the cohort level.

Figure 1. A cohort-based experimental design.

A cohort-based experimental design can facilitate determining the correlation between desiccation resistance, cuticular hydrocarbons (CHCs), and body weight. Five to six cohorts of each species were established.

Figure 1.

Figure 1—figure supplement 1. The phylogeny of 50 Drosophila and related species used in this study.

Figure 1—figure supplement 1.

Fifty species representing 46 Drosophila species (18 from Drosophila subgenus, 28 from Sophophora subgenus) and 4 species (3 Scaptodrosophila species and 1 Chymomyza species) as outgroup were used in our study. The 46 Drosophila species represent some of the main groups in the Drosophila genus: (1) melanogaster group, (2) obscura group, (3) willistoni group, (4) nasuta group, (5) virilis group, and (6) repleta group.
Figure 1—figure supplement 2. Chromatogram of the authentic standard n-alkane mixture which contains 34 n-alkanes from C7 to C40 at the same concentration.

Figure 1—figure supplement 2.

The authentic standard n-alkane mixture (C7–C40; Supelco 49,452U) was obtained from Sigma-Aldrich (St. Louis, MO). C7–C14 were excluded in the chromatogram due to the solvent delay in the gas chromatography–mass spectrometry (GC–MS) program. The integration of peak areas on longer cuticular hydrocarbons (CHCs) is biased to have lower values. To address this issue, we corrected the peak areas of CHCs based on the carbon-chain lengths using the integration from this standard n-alkane mixture.

Desiccation resistance varied across these 50 species with desert-dwelling species from the repleta group showing the highest desiccation resistance (e.g., D. mojavensis males: 58.2 ± 4.7 hr) and species from melanogaster group showing the lowest desiccation resistance (e.g., D. mauritiana males: 2.7 ± 0.2 hr) (Figure 2, Figure 2—figure supplement 1), consistent with the findings of previous research (Kellermann et al., 2012; Matzkin et al., 2009). GC–MS analyses of the CHCs in these 50 species detected five types of CHCs with different carbon-chain lengths and quantities, including linear alkanes (n-alkanes), methyl-branched alkanes (mbCHCs; the methyl branch is on the second carbon), monoenes (with one double bond), dienes (with two double bonds), and trienes (with three double bonds) (Figure 2—figure supplement 2). n-Alkanes, which have the highest melting temperatures, were only present in species from the melanogaster group, some of which have the lowest desiccation resistance among the species tested in this study. This suggests that n-alkanes may not have a general contribution to desiccation resistance across Drosophila species. Trienes, which have the lowest melting temperatures, were only present in eight species in low-to-moderate quantities. The other three types of CHCs, the mbCHCs, monoenes, and dienes, were observed in most species tested in our study.

Figure 2. Desiccation resistance and cuticular hydrocarbon (CHC) composition in 46 Drosophila species and 4 outgroup species.

Desiccation resistance and CHC composition of each species were plotted together. Males and females are shown separately. The shading intensity of each species represents the number of hours of desiccation resistance, while the size and color of each circle represent the type of CHC and its quantity. n-Alkanes are only present in some species from the melanogaster group.

Figure 2.

Figure 2—figure supplement 1. Desiccation resistance and body weight in 46 Drosophila species and 4 outgroup species.

Figure 2—figure supplement 1.

Desiccation resistance and body weight are shown in separate graphs.
Figure 2—figure supplement 2. The melting temperature of cuticular hydrocarbons (CHCs) is determined by methyl group, double bonds, and carbon-chain length.

Figure 2—figure supplement 2.

Five types of CHCs are detected in this study. Due to their chemical structures, n-alkanes have the highest melting temperatures followed by mbCHCs, monoenes, dienes, and trienes. Increasing length of the carbon-chain can also increase melting temperature. The positions of double bonds in this figure are for illustrative purposes and were not determined in this study.

Higher relative quantities of CHCs are not primarily responsible for higher desiccation resistance

We first sought to determine if higher desiccation resistance in flies could be due to having higher amounts of CHCs. One caveat in the analysis between absolute CHC quantity and desiccation resistance is a possible correlation between CHC quantities and body size: species with larger body size may possess higher quantities of CHCs. This can lead to potential biases when directly determining correlations between CHC quantities and desiccation resistance. Species with higher body weight can also have higher absolute water content and a lower surface area to volume ratio due to their larger size, which may also lead to higher desiccation resistance (Wang et al., 2021).

Our initial analyses showed a significant positive correlation (Females: r = 0.4, p < 0.001; Males: r = 0.5, p < 0.001) between body weight and desiccation (Figure 3A). We also found a significant positive correlation between body weight and total amount of CHCs (Female: r = 0.7, p < 0.001; Male: r = 0.7, p < 0.001) (Figure 3B), suggesting that variation in the body weight (or size) of different species may be a confound in determining the relationship between CHC quantity and desiccation resistance. To take this confound into consideration, we normalized the total CHC quantity by dividing by the body weight for each species. This shows that CHCs account for 0.02–0.5% of the total body weight of different species (Figure 3—figure supplement 1). Analyses between the normalized CHC quantity and desiccation resistance showed no correlation between these two variables in males (Male: p = 0.1) and positive correlation in females (Female: r = 0.1, p = 0.03) (Figure 3C). The low value of r (0.1) indicates a weak correlation between CHC quantities and desiccation resistance in females. This suggests that having higher amounts of CHCs only have a limited contribution to higher desiccation resistance.

Figure 3. Higher amounts of cuticular hydrocarbons (CHCs) do not contribute to higher desiccation resistance.

(A) Body weight is positively correlated with desiccation resistance (Females: r = 0.4, p < 0.001, Males: r = 0.5, p < 0.001). (B) Total amount of CHCs is correlated with having higher body weight (Females: r = 0.7, p < 0.001, Males: r = 0.7, p < 0.001). (C) A weak positive correlation between desiccation resistance and CHCs as a percentage of body weight in females, while no correlation in males (Females: r = 0.1, p = 0.03, Males: p = 0.1). All correlation analyses were conducted using Pearson’s method.

Figure 3.

Figure 3—figure supplement 1. Cuticular hydrocarbon (CHC) quantity as a percentage of body weight.

Figure 3—figure supplement 1.

As Drosophila species largely vary in size, we normalized the total CHCs as a percentage of body weight. CHCs account from 0.02% to 0.5% of the total body weight of different species.

mbCHCs are important determinants of desiccation resistance

As total CHC quantity is not a main contributor to desiccation resistance in our study, we hypothesized that the composition of CHC profiles may be important for desiccation resistance. The CHC layer of each Drosophila species is composed of different CHC components, each with different levels of abundance, so we used ‘beta diversity’ to represent this variation in the CHC composition of our sampled species, similar to what is used in microbiome analysis. This allows us to test the correlation between the compositional variation and desiccation resistance for the sampled Drosophila species. We then applied a permutational multivariate analysis of variance (PERMANOVA) to test whether the composition of CHCs differed across desiccation resistance (Anderson et al., 2006). Our results showed that the beta diversity of CHC compositions differs significantly across the increasing desiccation resistance (r2 = 0.1, p < 0.001; Figure 4—figure supplement 1), suggesting that CHC composition is important for desiccation resistance. We further sought to identify individual components of CHC profiles that are important determinants of desiccation resistance. To test if variation in CHC components can be used to determine desiccation resistance, we applied a random forest regression method that uses decision trees to connect the variables of interest (Liaw and Wiener, 2002; Svetnik et al., 2003), for example, CHC composition and desiccation resistance, and identified the importance of individual CHC components to predict desiccation resistance. In this analysis, we treated the CHC profiles of each species and sex as an individual dataset, giving us 100 datasets (50 species × 2 sexes) with five to six individual CHC profiles each. We then correlated the decision trees that generated from the CHC composition with desiccation resistance. The prediction process can identify key CHC components that are important contributors to desiccation resistance.

Random forest modeling of CHC composition was able to explain 85.5% of the variation in desiccation resistance (out of bag estimate: root mean square error ‘RMSE’ = 4.5) (Figure 4A). Models built with a 70:30 training/test split (n = 382, n = 164) performed similar to models built with out of bag estimate (RMSE = 5.4, Figure 4—figure supplement 2). Four mbCHCs were identified as having the highest contribution to predicting desiccation resistance in the regression model (listed in decreasing importance: 2MeC30, 2MeC28, 2MeC32, and 2MeC26), while many other CHCs did not substantially contribute to predicting desiccation resistance (Figure 4B). These results suggest that CHCs, and mbCHCs in particular, are important predictors of desiccation resistance in Drosophila species.

Figure 4. Cuticular hydrocarbon (CHC) composition can be used to predict desiccation resistance.

(A) Random Forest regression modeling of CHC abundance was able to explain 85.5% of the variation in time to desiccation with a root mean square error (RMSE) of 4.5. (B) The abundance of four mbCHCs, 2MeC30, 2MeC28, 2MeC32, and 2MeC26, has the highest importance to the desiccation resistance in the random forest regression model, while most of CHCs have less contribution to the accuracy of the model for desiccation resistance.

Figure 4.

Figure 4—figure supplement 1. The composition of cuticular hydrocarbons (CHCs) significantly differed across the increasing desiccation resistance.

Figure 4—figure supplement 1.

The non-metric multidimensional scaling (NMDS) plot showing the beta diversity of CHCs differs significantly across the increasing desiccation resistance (permutational multivariate analysis of variance [PERMANOVA]: r2 = 0.1, p < 0.001), suggesting CHC composition could affect desiccation resistance in the 50 Drosophila and selected species.
Figure 4—figure supplement 2. Cross-validation of the random forest regression model has a similar performance to the regression model using the full dataset.

Figure 4—figure supplement 2.

A 70:30 training/test split in the dataset (n = 382, n = 164) was used to evaluate the performance of the random forest regression model. The out of bag estimate of the root mean square error (RMSE) in the split dataset is 5.4, which is similar to the RMSE using the full dataset (RNMSE = 4.5).

Longer mbCHCs are correlated with higher desiccation resistance

We further sought to determine how mbCHCs contribute to desiccation resistance. As each species produces different combinations of these mbCHCs that are correlated with each other (Figure 5—figure supplement 1), we did not perform correlation analyses on these mbCHCs with a single regression model due to multicollinearity issues (Belsley et al., 2005). Instead, we performed correlation analyses between each of the most important CHCs identified from the random forest modeling, 2MeC26, 2MeC28, 2MeC30, and 2MeC32, and desiccation resistance. We found that the quantity of 2MeC26 negatively correlates with desiccation resistance in males (r = −0.4, p < 0.001) but no significant correlation was identified in females (p = 0.1) (Figure 5A). Significant positive correlations were observed between quantities of the longer mbCHCs and desiccation resistance for females (2MeC28: r = 0.4, p < 0.001; 2MeC30: r = 0.2, p = 0.009; 2MeC32: r = 0.3, p < 0.001) and males (2MeC28: r = 0.4, p < 0.001; 2MeC32: r = 0.4, p < 0.001) (Figure 5B–D). However, no significant correlation was observed between quantity of 2MeC30 and desiccation resistance in males (Figure 5C). These results suggested the production of longer mbCHCs could play a role in higher desiccation resistance.

Figure 5. Longer mbCHCs are associated with higher desiccation resistance.

(A) Quantity of 2MeC26 is negatively correlated with desiccation resistance in males, but no correlation in females (Females: p = 0.1, Males: r = −0.4, p < 0.001). (B) Quantity of 2MeC28 is positively correlated with desiccation resistance (Females: r = 0.4, p < 0.001, Males: r = 0.4, p < 0.001). (C) Quantity of 2MeC30 is positively correlated with desiccation resistance in females, but no correlation in males (Females: r = 0.2, p = 0.009, Males: p = 0.2). (D) Quantity of 2MeC32 is positively correlated with desiccation resistance (Females: r = 0.3, p < 0.001, Males: r = 0.4, p < 0.001). Correlations between the quantities of each mbCHC and desiccation resistance were determined using Pearson’s method. (E) D. melanogaster coated with longer mbCHCs have higher desiccation resistance. Coating of synthetic mbCHCs on D. melanogaster showed that 2MeC26 does not influence desiccation resistance. 2MeC30 increases desiccation resistance in female D. melanogaster while both 2MeC28 and 2MeC30 in increases desiccation resistance in male D. melanogaster. The bold horizontal line within each box plot corresponds to the median value, the box length to the interquartile range, and the lines emanating from the box (whiskers) extend to the smallest and largest observations. One-way analysis of variance (ANOVA) showed significant differences between D. melanogaster flies coated with different mbCHCs (Female: F(3,189) = 73.1, p < 0.001, Male: F(3,191) = 23.7, p < 0.001). Post hoc comparison was conducted using Tukey’s method. ***p < 0.001; n.s. refers to not significant.

Figure 5—source data 1. Desiccation resistance of Drosophila melanogaster attP40 flies when coated with individual mbCHC of 2MeC26, 2MeC28, and 2MeC30.
Figure 5—source data 2. Desiccation resistance of Drosophila melanogaster attP40 flies when coated with individual n-alkane of C23, C25, C27, C29, and C31.
Figure 5—source data 3. Desiccation resistance of Drosophila melanogaster CHC- flies when coated with individual CHC of 2MeC26, 2MeC28, 2MeC30, C23, C25, C27, C29, and C31.

Figure 5.

Figure 5—figure supplement 1. Correlation between different mbCHCs in 50 Drosophila and related species.

Figure 5—figure supplement 1.

Correlation between different mbCHCs was determined using Pearson’s method. Significant correlations between each pairs of mbCHCs were labeled with stars following with the correlation coefficients. The lines on the diagonals represent the distribution of the data for each mbCHC. **p < 0.01; ***p < 0.001.
Figure 5—figure supplement 2. The length of coated mbCHCs on D. melanogaster attP40 flies is positively correlated with desiccation resistance in both females and males.

Figure 5—figure supplement 2.

Pearson’s correlation: Female, r = 0.7, p < 0.001; Male, r = 0.5, p < 0.001.
Figure 5—figure supplement 3. Desiccation resistance of D. melanogaster attP40 flies coated with individual n-alkanes (C23, C25, C27, C29, and C31).

Figure 5—figure supplement 3.

One-way analysis of variance (ANOVA) was applied to test the differences in desiccation resistance between the coating of different n-alkanes (Female: F(5, 277) = 17.1, p < 0.001, Male: F(5, 287) = 8.0, p < 0.001), following with post hoc comparison using Tukey’s method at alpha = 0.05. *p < 0.05; **p < 0.01; ***p < 0.001; n.s. refers to not significant. The bold horizontal line within each box plot corresponds to the median value, the box length to the interquartile range, and the lines emanating from the box (whiskers) extend to the smallest and largest observations. Numbers on top of the boxplots are the mean in each treatment.
Figure 5—figure supplement 4. The length of coated n-alkanes on D. melanogaster attP40 flies is positively correlated with desiccation resistance in both females and males.

Figure 5—figure supplement 4.

Pearson’s correlation: Female, r = 0.4, p < 0.001; Male, r = 0.1, p = 0.02.
Figure 5—figure supplement 5. Desiccation resistance of CHC- D. melanogaster flies coated with individual mbCHCs (2MeC26, 2MeC28, 2MeC30) or n-alkanes (C23, C25, C27, C29, and C31).

Figure 5—figure supplement 5.

One-way analysis of variance (ANOVA) was applied to test the difference of desiccation resistance between the coating of different mbCHCs and n-alkanes (Female: F(8, 417) = 11.2, p < 0.001, Male: F(8, 408) = 7.1, p < 0.001), following with post hoc comparison using Tukey’s method at alpha = 0.05. *p < 0.05; **p < 0.01; ***p < 0.001. The bold horizontal line within each box plot corresponds to the median value, the box length to the interquartile range, and the lines emanating from the box (whiskers) extend to the smallest and largest observations. Numbers on top of the boxplots are the mean in each treatment.
Figure 5—figure supplement 6. The length of coated mbCHCs on CHC- flies has a weak positive correlation with desiccation resistance in females but did not correlate with the levels of desiccation resistance in males.

Figure 5—figure supplement 6.

Pearson’s correlation: Female, r = 0.2, p = 0.004; Male, p = 0.9.
Figure 5—figure supplement 7. The length of coated n-alkanes on CHC- flies did not correlate with desiccation resistance.

Figure 5—figure supplement 7.

Pearson’s correlation: Female, p = 0.4; Male, p = 1.0.
Figure 5—figure supplement 8. Pathway for the synthesis of branched cuticular hydrocarbons (mbCHCs) and linear CHCs (n-alkanes, monoenes, and dienes).

Figure 5—figure supplement 8.

All CHCs are synthesized via the fatty acyl-CoA synthesis pathway from acetyl-CoA. The pathway splits early due to the action of two fatty acid synthetases (FAS) into the branched CHC pathway or the linear CHC pathway. The precursors of each pathway are modified by synthesis enzymes such as elongases (Elo), reductases (FAR), and the terminal P450 decarbonylase (CYP4G) into CHCs. For linear CHCs, additional enzymes such as desaturases (Desat) incorporate double bonds during synthesis. The diagram was modified from Chung and Carroll, 2015.
Figure 5—figure supplement 9. Gas chromatography–mass spectrometry (GC–MS) chromatograms of 5′mFAS-GAL4, UAS-Cyp4g1-RNAi, and the F1 offspring from 5′mFAS-GAL4 × UAS-Cyp4g1-RNAi.

Figure 5—figure supplement 9.

CHC- D. melanogaster flies are generated by crossing 5′mFAS-GAL4 and UAS-Cyp4g1-RNAi, resulting in adult oenocyte specific knockdown of Cyp4g1 and loss of cuticular hydrocarbon (CHC) production. IS = internal standard spiked in for quantification (C26). cVA = cis-vaccenyl acetate, a lipid produced in the male ejaculatory bulb not affected by the RNAi knockdown of Cyp4g1 in the oenocytes.

Longer saturated CHCs can confer higher desiccation resistance in wild-type D. melanogaster

To investigate the above results suggesting that longer mbCHCs could underlie higher desiccation resistance in Drosophila species, we performed desiccation assays on D. melanogaster (attP40 strain) coated with synthetic mbCHCs of different lengths (2MeC26, 2MeC28, and 2MeC30). The D. melanogaster attP40 strain is a common lab strain used for transgenesis with a wild-type CHC profile. While previous studies showed that a mixture of 2MeC26, 2MeC28, and 2MeC30 can slightly increase desiccation resistance in D. melanogaster flies without CHCs, experiments were not performed on individual mbCHCs for desiccation resistance (Krupp et al., 2020). Our desiccation assays showed that the coating of 2MeC26 did not increase desiccation resistance in both sexes (post hoc comparison with Tukey’s method followed by one-way analysis of variance (ANOVA). Female: p = 0.07; Male: p = 1.0), while 2MeC30 significantly increased desiccation resistance in female D. melanogaster (t = 11.2, p < 0.001), and both 2MeC28 (t = 4.6, p < 0.001) and 2MeC30 (t = 6.9, p < 0.001) significantly increased desiccation resistance in male D. melanogaster (Figure 5E). There is a positive correlation with longer mbCHC length coated and higher desiccation resistance (Figure 5—figure supplement 2), suggesting that having higher quantities of longer mbCHCs can lead to higher desiccation resistance. As a control, we repeated the same experiment with synthetic n-alkanes (C23, C25, C27, C29, and C31). Coating with n-alkanes on the attP40 strain led to increases in desiccation resistance, which is positively correlated with carbon-chain length, similar to the mbCHCs experiments, suggesting that longer chain CHCs, with higher melting temperatures, can confer higher desiccation resistance when coated on the attP40 strain (Figure 5—figure supplement 3, Figure 5—figure supplement 4).

The coating of synthetic mbCHCs and n-alkanes on the attP40 strain showed that adding more CHCs to a strain with a wild-type CHC profile can increase desiccation resistance but did not show if these individual CHCs are sufficient to confer higher desiccation resistance. When these coating experiments (three mbCHCs, five n-alkanes) were individually repeated on D. melanogaster flies without CHCs (CHC-), we observed that while all mbCHCs and n-alkanes are able to confer higher desiccation resistance when coated on female CHC- flies, only a subset of these CHCs are able to confer resistance on male CHC- flies (Figure 5—figure supplement 5). In addition, only a positive correlation between longer mbCHCs and higher desiccation was observed in female CHC- flies. There is no correlation between mbCHC length and higher desiccation in male CHC- flies and no correlation between n-alkane length and higher desiccation resistance in female and male CHC- flies (Figure 5—figure supplement 6, Figure 5—figure supplement 7). We noted that coating mbCHCs and n-alkanes only produced a small increase in desiccation resistance in CHC- flies, compared to coating these CHCs on the attP40 strain. This is similar to a previous study by Krupp et al., 2020 where the authors showed similar small increases in desiccation resistance when only n-alkanes or mbCHCs are coated on CHC- flies, compared to the larger increase in desiccation resistance when wild-type extracts of CHCs were coated. Taken together, these results suggest that CHCs with higher melting temperatures such as longer mbCHCs or n-alkanes contribute to higher desiccation resistance when in a blend of CHCs but not sufficient to confer high desiccation resistance by themselves. This is consistent with the hypotheses from multiple papers that the cuticular lipid layer works as a blend or a mixture rather than as individual components (Gibbs, 1998; Menzel et al., 2019; Wigglesworth, 1945). Although we showed that both longer mbCHCs and n-alkanes can confer higher desiccation resistance when coated on D. melanogaster with a wild-type profile, as n-alkanes are mostly absent in the Drosophila genus, our experiments suggest a causal relationship between longer mbCHCs and higher desiccation resistance in Drosophila.

Evolution of mbCHCs underlies variation in desiccation resistance in Drosophila species

We have shown that mbCHCs are important determinants of desiccation resistance in Drosophila species and experimentally coating insects with longer mbCHCs confers higher desiccation resistance. To further determine whether the evolution of mbCHCs underlies the variation in desiccation resistance in Drosophila species, we first examined the evolutionary trajectory of mbCHCs and then tested the correlation between mbCHCs and desiccation resistance with phylogenetic correction. We mapped the evolution of mbCHC composition in the Drosophila genus using ancestral trait reconstruction and tested the phylogenetic signal using Pagel’s λ (Freckleton et al., 2002; Goolsby et al., 2017). Ranging from 0 to 1, Pagel’s λ measures the extent to which the variance of the trait can be explained by the phylogeny and therefore determines the degree of association between the trait evolution and the phylogeny (Lynch, 1991). If λ is closer to 1, the evolution of the trait has a higher association with the phylogeny. Phylogenetic signals in mbCHCs were detected for both females and males with λ = 0.75 and 0.83, respectively (Supplementary file 1a). This suggests a moderate to strong association between mbCHC evolution and phylogeny. The derived ancestral state for mbCHCs in the Drosophila genus suggest that 2MeC28 and 2MeC30 are the major mbCHCs (Figure 6—figure supplement 1). The evolutionary trajectory of mbCHC composition showed repeated evolution in the lengths of mbCHCs, as well as their quantities. Independent evolution of longer mbCHCs was observed in four clades of species, including the willistoni, nasuta, and repleta groups, and the ananassae subgroup (including D. ananassae, D. pseudoananassae, and D. bipectinata), while shorter mbCHCs were mainly observed in several species in the melanogaster subgroup (e.g., D. melanogaster, D. simulans, and D. teissieri) (Figure 2, Figure 6—figure supplement 1).

Desiccation resistance also has a strong association with the phylogeny of Drosophila species (Supplementary file 1b; Kellermann et al., 2018; Kellermann et al., 2012). We sought to determine if the variation of desiccation resistance could be explained by the evolution of mbCHCs in these species. We used a Phylogenetic Generalized Linear Square (PGLS) model to test for the correlation between mbCHCs and desiccation resistance when controlling the effects from the phylogenetic relationship of the species (Grafen, 1989; Mundry, 2014). Since mbCHCs are produced in a linear pathway with each species producing several different mbCHCs, we selected the mbCHC with the longest carbon-chain length in the species as the proxy to represent the lengths of mbCHCs in each species. We incorporated the quantity, length, and their interaction in the PGLS model for the longest mbCHC. The interaction term between the length and quantity of the longest mbCHC in the model can determine how the two variables combinatorically affect desiccation resistance. PGLS modeling showed, after correcting for phylogenetic effects, the higher quantity and longer length of the longest mbCHC both affect desiccation resistance (interaction term, Female: t = 3.5, p < 0.001; Male: t = 2.2, p = 0.03) (Supplementary file 1c; Figure 6). This suggests that the synthesis of mbCHCs with longer carbon-chain lengths could be a common mechanism underlying the evolution of higher desiccation resistance.

Figure 6. Evolution of longer mbCHCs is significantly correlated with evolution of higher desiccation resistance.

Patterns in the normalized quantities of the longest mbCHCs and desiccation resistance for females (top) and males (bottom) were listed across the phylogeny of the 50 Drosophila and related species. Phylogenetic Generalized Linear Square (PGLS) analysis between the longest mbCHCs and desiccation resistance showed both the higher quantity and longer length of the longest mbCHC affect desiccation resistance (interaction term, Female: t = 3.5, p < 0.001; Male: t = 2.2, p = 0.03).

Figure 6—source data 1. Normalized quantities of mbCHCs, the length and normalized quantities of the longest mbCHC in each species.

Figure 6.

Figure 6—figure supplement 1. Ancestral state reconstruction of mbCHC for the Drosophila genus.

Figure 6—figure supplement 1.

Bar plots on the right were the derived ancestral state of female and male mbCHC profiles for the common ancestor of the Drosophila genus (A), the Sophophora subgenus (B), the Drosophila subgenus (C), the melanogaster group (D), the willistoni group (E), the nasuta group (F), the virilis group (G), the repleta group (H), the clade of D. melanogaster, D. simulans, and D. mauritiana (I), and the clade of D. mojavensis and D. arizonae (J). 2MeC28 and 2MeC30 were the major cuticular hydrocarbons (CHCs) in the derived ancestral state (A).

Discussion

Reducing evaporative water loss through the cuticle using a layer of CHCs is one of the most important evolutionary innovations in insects that allows many species to survive and thrive in diverse and arid habitats. We showed that CHC composition can account for 85.5% of the variation in desiccation resistance in the 50 Drosophila and related species in this study. This suggests that CHC composition may be highly predictive of desiccation resistance. Algorithmic ranking in importance of CHC components using a random forest machine-learning model showed that mbCHCs have the highest contribution to determining desiccation resistance. Importantly, higher amounts of longer mbCHCs are important in desiccation resistance. This is consistent with previous studies showing that mbCHCs have higher melting temperatures than the other commonly present types of CHCs (monoenes and dienes) and longer CHC leads to higher melting temperature (Gibbs and Pomonis, 1995). Our overall analyses support the conclusion that the evolution of longer mbCHCs is correlated with higher desiccation resistance in our dataset, but we observed minor differences between males and females (Figures 3 and 5). In particular, while higher quantities of 2MeC28 and 2MeC32 are both positively correlated with higher desiccation resistance in both males and females (Figure 5B, D), there are differences between the sexes regarding 2MeC26 and 2MeC30 (Figure 5A, C). There could be several underlying reasons for this. Firstly, across our dataset, CHCs are sexually dimorphic in many species, leading to overall differences in the correlation analyses between the sexes. Secondly, other physiological differences, such as size differences, may contribute to the differences between males and females. In many species such as D. melanogaster, females are larger than males and this may contribute to the higher basal desiccation resistance (when CHCs are taken out of the equation) of female CHC- flies (3.6 ± 0.1 hr) compared to that of male CHC- flies (2.1 ± 0.1 hr) (Figure 5E). While these are minor differences in the overall correlation analyses between males and females, the synthetic CHC coating experiments showed that coating of longer mbCHCs confers higher desiccation resistance in both males and females.

mbCHCs and desiccation resistance

Previous studies showed support for mbCHCs in desiccation resistance in Drosophila. RNAi knockdown of a methyl-branched specific fatty acid synthase (mFAS) in D. serrata eliminates almost all mbCHCs and leads to significant decrease in desiccation resistance which could be partially rescued with synthetic mbCHCs (Chung et al., 2014) . Its sibling species, the rainforest desiccation sensitive D. birchii, have only trace amounts of mbCHCs (Howard et al., 2003) and is unable to evolve desiccation resistance despite strong laboratory selection over many generations (Hoffmann et al., 2003). The presence of mbCHCs in most Drosophila and related species leads to another question: Is this a general rule for insects using mbCHCs to develop desiccation resistance? A high proportion of very-long-chain mbCHCs with a single methyl branch have been reported in many insect species from different taxa that dwell in desert environments, such as the desert tenebrionid beetle Eleodes armata (Hadley, 1977), the desert ants Cataglyphis niger and P. barbatus (Johnson and Gibbs, 2004; Soroker and Hefetz, 2000), and the desert locust Schistocerca gregaria (Heifetz et al., 1998). Although a small proportion of n-alkanes was also reported in some of these species, this suggests the use of mbCHCs to minimize water evaporation may be a common mechanism for insects in extremely dry environments.

Ancestral trait reconstruction showed that the derived ancestral status for mbCHCs in the Drosophila genus is 2MeC28 and 2MeC30 as the major mbCHCs (Figure 6—figure supplement 1). As there are many extant species that do not have 2MeC30, this suggests that the ability to synthesize longer mbCHCs may be lost during evolution, especially for species that do not inhabit arid environments and have low levels of desiccation resistance such as D. melanogaster and D. biarmipes. On the contrary, for species that dwell in extremely arid environments such as D. mojavensis and D. arizonae, higher quantities of 2MeC30 and an even longer mbCHC, 2MeC32, are observed. The observation is consistent with the results from our PGLS modeling, showing significant correlation between having higher quantity of longer mbCHCs and higher desiccation resistance in the Drosophila genus.

Constraints in evolution of CHCs and desiccation resistance

Why are n-alkanes, which have the highest melting temperature and potentially better water-proofing properties, not prevalent in species in our studies? A hypothesis that could explain this is the competition on the common precursors between the synthesis of linear unsaturated CHCs (monoenes and dienes) and n-alkanes. In insects, the synthesis of all CHCs start from the fatty acyl CoA synthesis pathway. The pathway is split where a cytosolic fatty acid synthase (cFAS) will synthesize all the linear CHCs (n-alkanes, monoenes, dienes), while mFAS will synthesize mbCHCs (Figure 5—figure supplement 8; Chung and Carroll, 2015; Holze et al., 2021). The pathway suggests that the synthesis of n-alkanes, monoenes, and dienes are in competition. As many monoenes and dienes function in Drosophila as contact pheromones in regulating different behaviors such as mating and aggression (Blomquist and Ginzel, 2021; Chertemps et al., 2007; Chertemps et al., 2006; Chung and Carroll, 2015; Krupp et al., 2008; Wang et al., 2011), the synthesis for these CHCs may compete with the synthesis of n-alkanes. In contrast, the synthesis of mbCHCs is in another part of the pathway. This suggests that there are potential constraints in the roles of n-alkanes in desiccation resistance in Drosophila due to competition with CHCs which function as signaling molecules, while the synthesis of mbCHCs is relatively unaffected by other products of the pathway (Figure 5—figure supplement 8; Chung and Carroll, 2015). We speculate that in Drosophila and related species, the use of mbCHCs for modulating desiccation resistance may avoid the conflict between surviving desiccation and chemical signaling.

If there are no biosynthetic constraints to the synthesis of mbCHCs, would the evolution of longer mbCHCs be a general mechanism in the evolution of higher desiccation for species adapting to more arid climates? We suggest that this may be the case for many species, but for some species that use the shorter mbCHCs such as 2MeC26 and 2MeC28 as signaling molecules, such as D. serrata (Chenoweth and Blows, 2005; Chung et al., 2014) and the longhorned beetle, Mallodon dasystomus (Spikes et al., 2010), natural selection and sexual selection may have opposing effects on the lengths of these mbCHCs. In D. serrata, having higher amounts of 2MeC26 leads to higher mating success in male D. serrata (Chung et al., 2014). Northern populations of this species produce lower amounts of 2MeC26 but higher amounts of 2MeC30, due to a polymorphism in a fatty acyl-CoA reductase gene that shifts to an overall production of longer CHCs. As CHCs are sexually monomorphic in D. serrata, this polymorphism led to lower mating success in male northern flies and higher desiccation resistance in female northern flies compared to their common counterparts (Rusuwa et al., 2022). The sexually antagonistic effect due to the functions of the different mbCHCs may affect their evolutionary response to desiccation in many species.

CHC variation and desiccation resistance across insects

How much of what we learnt in this study is applicable broadly to natural populations of Drosophila and other insect species? The use of fly strains from the stock center reared on a common cornmeal artificial medium and the same abiotic factors rather than wild-caught flies allowed us to standardize this large-scale experiment. Previous studies have shown that diet (Fedina et al., 2012), temperature (Rajpurohit et al., 2021), and mating status (Everaerts et al., 2010) can affect CHC profiles in D. melanogaster, therefore the CHC profiles that we obtained in our studies may not be representative of the CHC variations present in some of these species in their natural habitats. However, these differences in CHC profiles are mostly quantitative, that is, different quantities of the same CHCs, rather than different CHCs. Similarly, CHC variations between different natural populations of the same Drosophila species are also largely quantitative than qualitative (Dembeck et al., 2015; Rusuwa et al., 2022). Therefore, stock center strains would be representative of the qualitative CHC variations between species and our conclusions we reach in this study may be extended to wild populations. The use of standardized conditions also precludes the selection of species that require specialized media and conditions (e.g., species from the Hawaiian radiation and mycophagus species). While these are some of the limitations of our study, the standardized conditions allow us to reduce the number of confounds in our analyses for the correlation between CHC variation and desiccation resistance. As our study focused on a broad sampling of Drosophila and related species from different habitats spanning millions of years of evolution, the overall interpretation of the data and the broad conclusion that we reached regarding the correlation between the evolution of mbCHCs and desiccation resistance is likely to extend to natural populations of the Drosophila genus.

Constraints in the evolution of CHCs due to their roles in many insect species as contact pheromones are likely to affect their evolution in response to desiccation resistance. In this study, we showed that the evolution of mbCHCs are major determinants of desiccation resistance across Drosophila species, as n-alkanes shared the same branch in the fatty acyl-CoA pathway as monoenes and dienes (Figure 5—figure supplement 8), many of which function as contact pheromones in Drosophila (Ferveur, 2005). In other insect orders such as Lepidoptera, where volatile non-CHC pheromones are predominantly used for chemical communication (Ando et al., 2004), constraints on each class of CHCs may be different from Drosophila. Likewise, across eusocial hymenopterans (Van Oystaeyen et al., 2014), saturated CHCs such as n-alkanes and mbCHCs are used as queen pheromones and are likely to be under constraints in the evolution of desiccation resistance in these species. Finally, in some rare examples such as dragonflies where very-long-chain methyl ketones and aldehydes are the major cuticular wax components (Futahashi et al., 2019), the evolution of CHCs may not be a major factor in the evolution of desiccation resistance. While we limit our conclusions that the synthesis of mbCHCs with longer carbon-chain lengths is a common mechanism underlying the evolution of higher desiccation resistance to Drosophila, the evolution of cuticular components that modulate the melting temperature of the cuticular lipid layer can reduce water loss and increase desiccation resistance. We suggest that the evolution of CHC components in insects can be largely predictive of desiccation resistance between closely related species.

Materials and methods

Drosophila species and strains

In this study, 46 Drosophila species, as well as three Scaptodrosophila species and one Chymomyza species were either obtained from the National Drosophila Species Stock Center (NDSSC) or were gifts from various colleagues. Details are listed in Supplementary file 1d. The D. melanogaster y w; attP40 strain and UAS-Cyp4g1 RNAi strain were obtained from the Perrimon lab and the Bloomington Drosophila Stock Center, respectively. The 5′mFAS-GAL4 strain was constructed by cloning the oenocyte enhancer fragment of the mFAS/FASN2 (CG3524) gene into the GAL4 vector pBPGUw (Addgene #17575) using primers 5′CG3524-TOPO-F (5′-CACCCCGCGGCGTGTTATTGAACC-3′) and 5′CG3524-TOPO-R (5′-CTTGTTGCGCAGACAGACTG-3′), before injecting into the D. melanogaster strain 51C and integrated into the genome using the PhiC31 integrase system. All species were reared on standard cornmeal medium (Flystuff 66-121 Nutri-Fly Bloomington Formulation). The phylogenetic relationship of all 50 species in this study was adapted from Finet et al., 2021. CHC- D. melanogaster were generated by crossing the 5′mFAS-GAL4 driver line (which expresses GAL4 in adult oenocytes) with the UAS-Cyp4g1 RNAi strain (Figure 5—figure supplement 9).

Experimental design

To investigate the contribution of CHCs to desiccation resistance, a cohort-based design was used. For each species, five to six cohorts were established and three measurements were conducted for each sex of the F1 progeny, including desiccation resistance, CHCs, and body weight (Figure 1). Each cohort was treated as a biological replicate. Each cohort in each species was established by pooling five females and five males on a standard cornmeal medium in the environmental chambers set at 25°C and a 12L:12D photoperiod. To maximize the food and spatial availability and minimize competition between the F1 progenies (Mueller, 1988), the parent flies in each cohort were transferred to fresh food after 5 days. F1 progeny flies were collected daily following emergence, separated by sex, and maintained on fresh cornmeal medium. All flies used for measurements were 4- to 5-day old.

Desiccation resistance assays

Desiccation resistance assays were performed in a randomized and blinded manner with the setup consistent with a previously published protocol (Figure 1; Chung et al., 2014). Briefly, in each cohort, 10 adults of the same sex were subjected to the assay setup containing 10 g of silica gel (S7500-1KG, Sigma-Aldrich, St. Louis, MO). After the assays were assembled, all the setups were randomly arranged with a number assigned. Mortality of the flies was recorded hourly after 2 hr. For each cohort, one vial was scored and the average time in hours until all flies died was recorded as desiccation resistance. Desiccation resistance assays for the 50 Drosophila and related species were conducted at 25°C and with a 12L:12D photoperiod, while the desiccation assays for the CHC coating experiments were conducted at 27°C (GAL4/UAS) and a 12L:12D photoperiod.

CHC analyses

CHCs were extracted and analyzed using GC–MS following previously published protocols (Lamb et al., 2020; Savage et al., 2021). Five flies of the same sex (4- to 5-day old) from each cohort were soaked for 10 min in 200 µl hexane containing hexacosane (C26; 25 ng/µl) as an internal standard. Extracts were directly analyzed by GC–MS (7890A, Agilent Technologies Inc, Santa Clara, CA) using a DB-17ht column (Agilent Technologies Inc, Santa Clara, CA). To identify the CHC composition, we first compared retention times and mass spectra to an authentic standard mixture (C7–C40) (Supelco 49,452U, Sigma-Aldrich, St. Louis, MO) with CHC samples. The types of CHCs, including methyl-branched alkanes, monoenes, dienes, and trienes, were then identified by a combination of their specific fragment ions on the side of functional groups (methyl branch or double bonds), retention times relative to linear hydrocarbon standards, and the m/z value of the molecular ion. The position of methyl branch in mbCHCs was determined using the protocol described in Carlson et al., 1998, while the position of double bonds was not determined in this study. Each CHC peak was quantified using its comparison with the peak area of the internal standard (C26) and represented as nanogram per fly (ng/fly). Because we observed a biased integration of peak areas on longer CHCs in running the standard mixture (Figure 1—figure supplement 2), we corrected the peak areas of CHCs based on the carbon-chain lengths using the integration from the standard mixture.

Body weight measurement

Body weight was incorporated in the data collection and analysis. The body weight was determined as the difference between the Eppendorf tube containing 5–10 flies of the same sex and the same empty Eppendorf tube.

Coating of synthetic compounds

To coat each of synthetic mbCHCs and n-alkanes on D. melanogaster, including 2MeC26, 2MeC28, 2MeC30, C23, C25, C27, C29, and C31, we first added 200 µl hexane containing 300 ng/µl of each CHC in a 2-ml glass vial and then used a nitrogen evaporator (BT1603 G-Biosciences) to evaporate the hexane with only mbCHCs precipitate at the bottom of the vial. For the control group, 200 µl hexane was added. After the hexane was evaporated, 10 flies of the same sex were transferred into each vial, following shaking on a Vortex for 20 s on, 20 s off, and 20 s on, following previous protocols (Billeter et al., 2009; Chung et al., 2014). The flies were then directly subjected to desiccation assays. Synthetic mbCHCs were kindly provided by Dr. Jocelyn Millar (University of California, Riverside) and synthetic n-alkanes were ordered from Sigma-Aldrich (C23 #91447-1G, C25 #286931-1G, C27 #51559-250MG, C29 #74156-250MG, and C31 #51529-250MG).

Statistics

All analyses were conducted in R (Version 4.1). Correlation analyses were conducted with Pearson’s method using ‘cor.test’ function. The dependent variables were log-transformed to better conform with assumptions of normality. Variance in CHC beta diversity across desiccation resistance was determined using PERMANOVA using the ‘vegan’ package (Anderson et al., 2006). The random forest regression analysis was performed using the ‘ranger’ and ‘randomForest’ packages (Liaw and Wiener, 2002; Wright and Ziegler, 2015). The random forest regression models were built using both Out Of Bag estimate and test/training sets (70:30 split). To determine how useful each CHC variable is in the prediction of desiccation resistance in the random forest regression analysis, the importance of top predictor CHCs was quantified using permutation importance (Altmann et al., 2010). Desiccation resistance in flies coated with different mbCHCs was determined using one-way ANOVA at alpha = 0.05. Post hoc comparison was further conducted using Tukey’s method. The ancestral trait reconstruction and estimation of phylogenetic signal, Pagel’s λ, for mbCHC composition and desiccation resistance were determined using ‘Phytools’, ‘Picante’, and ‘Rphylopars’ packages (Goolsby et al., 2017; Kembel et al., 2010; Revell, 2012). The PGLS analysis was conducted using generalized least squares fit by maximum likelihood and the covariance structure between species was used under a Brownian motion process of evolution. The PGLS analyses were conducted using ‘GLS’ function in ‘ape’ package (Paradis and Schliep, 2019).

Acknowledgements

We thank Dr. Jocelyn Millar for guidance and advice, as well as the synthetic compounds used in this study. We also thank Ye Ma, Elaina Giannetti, Taylor Hori, and Zhuo Chen for technical assistance, Yuzhang Shan for assistance with figure visualization, and the National Drosophila Species Stock Center for fly stocks of these different Drosophila species. This work is supported by a National Science Foundation grant (2054773) to H Chung.

Funding Statement

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Contributor Information

Zinan Wang, Email: wangzina@msu.edu.

Henry Chung, Email: hwchung@msu.edu.

Lynn M Riddiford, University of Washington, United States.

K VijayRaghavan, National Centre for Biological Sciences, Tata Institute of Fundamental Research, India.

Funding Information

This paper was supported by the following grant:

  • National Science Foundation 2054773 to Henry Chung.

Additional information

Competing interests

No competing interests declared.

Author contributions

Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing.

Software, Formal analysis, Validation, Methodology, Writing - review and editing.

Data curation, Validation, Investigation.

Data curation, Investigation.

Data curation, Investigation.

Formal analysis, Methodology, Writing - review and editing.

Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing.

Additional files

Supplementary file 1. Statistical output from phylogenetic comparative analyses and the list of species used in this study.

(a) Phylogenetic signals for mbCHCs in Drosophila species. (b) Phylogenetic signals for desiccation resistance in Drosophila species. (c) Summary of the Phylogenetic Generalized Linear Square (PGLS) models between the longest mbCHCs and desiccation resistance for females and males in 50 Drosophila and related species. (d) List of species used in this study.

elife-80859-supp1.docx (71.5KB, docx)
MDAR checklist
Source data 1. Data of normalized cuticular hydrocarbon (CHC) quantities, body weight, averaged desiccation resistance, total CHC quantity, and percentage of CHC quantities per body weight.
elife-80859-data1.xlsx (168.7KB, xlsx)
Source code 1. Code for Pearson’s correlation analysis and analysis of variance (ANOVA) tests in Figure 3 and Figure 5.
elife-80859-code1.zip (23.4KB, zip)
Source code 2. Code for NMDS plot and random forest analysis for Figure 4.
elife-80859-code2.zip (8.4KB, zip)
Source code 3. Code for Phylogenetic Generalized Linear Square (PGLS) analyses in Figure 6.
elife-80859-code3.zip (5.4KB, zip)

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting source data file. Code used is uploaded as source codes 1-3.

References

  1. Altmann A, Toloşi L, Sander O, Lengauer T. Permutation importance: a corrected feature importance measure. Bioinformatics. 2010;26:1340–1347. doi: 10.1093/bioinformatics/btq134. [DOI] [PubMed] [Google Scholar]
  2. Anderson MJ, Ellingsen KE, McArdle BH. Multivariate dispersion as a measure of beta diversity. Ecology Letters. 2006;9:683–693. doi: 10.1111/j.1461-0248.2006.00926.x. [DOI] [PubMed] [Google Scholar]
  3. Ando T, Inomata S, Yamamoto M. In: The Chemistry of Pheromones and Other Semiochemicals I. Schulz S, editor. Springer; 2004. Lepidopteran sex pheromones; pp. 51–96. [DOI] [PubMed] [Google Scholar]
  4. Belsley DA, Kuh E, Welsch RE. Regression diagnostics: Identifying influential data and sources of collinearity. John Wiley & Sons; 2005. [DOI] [Google Scholar]
  5. Billeter JC, Atallah J, Krupp JJ, Millar JG, Levine JD. Specialized cells tag sexual and species identity in Drosophila melanogaster. Nature. 2009;461:987–991. doi: 10.1038/nature08495. [DOI] [PubMed] [Google Scholar]
  6. Blomquist GJ, Ginzel MD. Chemical ecology, biochemistry, and molecular biology of insect hydrocarbons. Annual Review of Entomology. 2021;66:45–60. doi: 10.1146/annurev-ento-031620-071754. [DOI] [PubMed] [Google Scholar]
  7. Buellesbach J, Whyte BA, Cash E, Gibson JD, Scheckel KJ, Sandidge R, Tsutsui ND. Desiccation resistance and micro-climate adaptation: cuticular hydrocarbon signatures of different Argentine ant supercolonies across California. Journal of Chemical Ecology. 2018;44:1101–1114. doi: 10.1007/s10886-018-1029-y. [DOI] [PubMed] [Google Scholar]
  8. Bujan J, Yanoviak SP, Kaspari M. Desiccation resistance in tropical insects: causes and mechanisms underlying variability in a Panama ant community. Ecology and Evolution. 2016;6:6282–6291. doi: 10.1002/ece3.2355. [DOI] [PMC free article] [PubMed] [Google Scholar]
  9. Carlson DA, Bernier UR, Sutton BD. Elution patterns from capillary GC for methyl-branched alkanes. Journal of Chemical Ecology. 1998;24:1845–1865. doi: 10.1023/A:1022311701355. [DOI] [Google Scholar]
  10. Chen N, Fan Y-L, Bai Y, Li X-D, Zhang Z-F, Liu T-X. Cytochrome P450 gene, CYP4G51, modulates hydrocarbon production in the pea aphid, Acyrthosiphon pisum. Insect Biochemistry and Molecular Biology. 2016;76:84–94. doi: 10.1016/j.ibmb.2016.07.006. [DOI] [PubMed] [Google Scholar]
  11. Chen N, Pei XJ, Li S, Fan YL, Liu TX. Involvement of integument‐rich CYP4G19 in hydrocarbon biosynthesis and cuticular penetration resistance in Blattella germanica (L.) Pest Management Science. 2019;76:215–226. doi: 10.1002/ps.5499. [DOI] [PubMed] [Google Scholar]
  12. Chenoweth SF, Blows MW. Contrasting mutual sexual selection on homologous signal traits in Drosophila serrata. The American Naturalist. 2005;165:281–289. doi: 10.1086/427271. [DOI] [PubMed] [Google Scholar]
  13. Chertemps T, Duportets L, Labeur C, Ueyama M, Wicker-Thomas C. A female-specific desaturase gene responsible for diene hydrocarbon biosynthesis and courtship behaviour in Drosophila melanogaster. Insect Molecular Biology. 2006;15:465–473. doi: 10.1111/j.1365-2583.2006.00658.x. [DOI] [PubMed] [Google Scholar]
  14. Chertemps T, Duportets L, Labeur C, Ueda R, Takahashi K, Saigo K, Wicker-Thomas C. A female-biased expressed elongase involved in long-chain hydrocarbon biosynthesis and courtship behavior in Drosophila melanogaster. PNAS. 2007;104:4273–4278. doi: 10.1073/pnas.0608142104. [DOI] [PMC free article] [PubMed] [Google Scholar]
  15. Chiang YN, Tan KJ, Chung H, Lavrynenko O, Shevchenko A, Yew JY. Steroid hormone signaling is essential for pheromone production and oenocyte survival. PLOS Genetics. 2016;12:e1006126. doi: 10.1371/journal.pgen.1006126. [DOI] [PMC free article] [PubMed] [Google Scholar]
  16. Chown SL, Hoffmann AA, Kristensen TN, Angilletta MJ, Stenseth NC, Pertoldi C. Adapting to climate change: a perspective from evolutionary physiology. Climate Research. 2010;43:3–15. doi: 10.3354/cr00879. [DOI] [Google Scholar]
  17. Chown SL, Sørensen JG, Terblanche JS. Water loss in insects: an environmental change perspective. Journal of Insect Physiology. 2011;57:1070–1084. doi: 10.1016/j.jinsphys.2011.05.004. [DOI] [PubMed] [Google Scholar]
  18. Chung H, Loehlin DW, Dufour HD, Vaccarro K, Millar JG, Carroll SB. A single gene affects both ecological divergence and mate choice in Drosophila. Science. 2014;343:1148–1151. doi: 10.1126/science.1249998. [DOI] [PubMed] [Google Scholar]
  19. Chung H, Carroll SB. Wax, sex and the origin of species: dual roles of insect cuticular hydrocarbons in adaptation and mating. BioEssays. 2015;37:822–830. doi: 10.1002/bies.201500014. [DOI] [PMC free article] [PubMed] [Google Scholar]
  20. Dembeck LM, Böröczky K, Huang W, Schal C, Anholt RRH, Mackay TFC. Genetic architecture of natural variation in cuticular hydrocarbon composition in Drosophila melanogaster. eLife. 2015;4:e09861. doi: 10.7554/eLife.09861. [DOI] [PMC free article] [PubMed] [Google Scholar]
  21. Everaerts C, Farine JP, Cobb M, Ferveur JF. Drosophila cuticular hydrocarbons revisited: mating status alters cuticular profiles. PLOS ONE. 2010;5:e9607. doi: 10.1371/journal.pone.0009607. [DOI] [PMC free article] [PubMed] [Google Scholar]
  22. Fedina TY, Kuo TH, Dreisewerd K, Dierick HA, Yew JY, Pletcher SD. Dietary effects on cuticular hydrocarbons and sexual attractiveness in Drosophila. PLOS ONE. 2012;7:e49799. doi: 10.1371/journal.pone.0049799. [DOI] [PMC free article] [PubMed] [Google Scholar]
  23. Ferveur JF. Cuticular hydrocarbons: their evolution and roles in Drosophila pheromonal communication. Behavior Genetics. 2005;35:279–295. doi: 10.1007/s10519-005-3220-5. [DOI] [PubMed] [Google Scholar]
  24. Ferveur JF, Cortot J, Rihani K, Cobb M, Everaerts C. Desiccation resistance: effect of cuticular hydrocarbons and water content in Drosophila melanogaster adults. PeerJ. 2018;6:e4318. doi: 10.7717/peerj.4318. [DOI] [PMC free article] [PubMed] [Google Scholar]
  25. Finet C, Kassner VA, Carvalho AB, Chung H, Day JP, Day S, Delaney EK, De Ré FC, Dufour HD, Dupim E, Izumitani HF, Gautério TB, Justen J, Katoh T, Kopp A, Koshikawa S, Longdon B, Loreto EL, Nunes MDS, Raja KKB, Rebeiz M, Ritchie MG, Saakyan G, Sneddon T, Teramoto M, Tyukmaeva V, Vanderlinde T, Wey EE, Werner T, Williams TM, Robe LJ, Toda MJ, Marlétaz F. DrosoPhyla: resources for drosophilid phylogeny and systematics. Genome Biology and Evolution. 2021;13:evab179. doi: 10.1093/gbe/evab179. [DOI] [PMC free article] [PubMed] [Google Scholar]
  26. Freckleton RP, Harvey PH, Pagel M. Phylogenetic analysis and comparative data: a test and review of evidence. The American Naturalist. 2002;160:712–726. doi: 10.1086/343873. [DOI] [PubMed] [Google Scholar]
  27. Futahashi R, Yamahama Y, Kawaguchi M, Mori N, Ishii D, Okude G, Hirai Y, Kawahara-Miki R, Yoshitake K, Yajima S, Hariyama T, Fukatsu T. Molecular basis of wax-based color change and UV reflection in dragonflies. eLife. 2019;8:e43045. doi: 10.7554/eLife.43045. [DOI] [PMC free article] [PubMed] [Google Scholar]
  28. Gibbs AG, Pomonis JG. Physical properties of insect cuticular hydrocarbons: the effects of chain length, methyl-branching and unsaturation. Comparative Biochemistry and Physiology Part B. 1995;112:243–249. doi: 10.1016/0305-0491(95)00081-X. [DOI] [Google Scholar]
  29. Gibbs AG, Chippindale AK, Rose MR. Physiological mechanisms of evolved desiccation resistance in Drosophila melanogaster. The Journal of Experimental Biology. 1997;200:1821–1832. doi: 10.1242/jeb.200.12.1821. [DOI] [PubMed] [Google Scholar]
  30. Gibbs AG. Water-proofing properties of cuticular lipids. American Zoologist. 1998;38:471–482. doi: 10.1093/icb/38.3.471. [DOI] [Google Scholar]
  31. Gibbs AG, Matzkin LM. Evolution of water balance in the genus Drosophila. The Journal of Experimental Biology. 2001;204:2331–2338. doi: 10.1242/jeb.204.13.2331. [DOI] [PubMed] [Google Scholar]
  32. Gibbs AG. Lipid melting and cuticular permeability: new insights into an old problem. Journal of Insect Physiology. 2002a;48:391–400. doi: 10.1016/s0022-1910(02)00059-8. [DOI] [PubMed] [Google Scholar]
  33. Gibbs AG. Water balance in desert Drosophila: lessons from non-charismatic microfauna. Comparative Biochemistry and Physiology. Part A, Molecular & Integrative Physiology. 2002b;133:781–789. doi: 10.1016/s1095-6433(02)00208-8. [DOI] [PubMed] [Google Scholar]
  34. Gibbs AG, Fukuzato F, Matzkin LM. Evolution of water conservation mechanisms in Drosophila. The Journal of Experimental Biology. 2003;206:1183–1192. doi: 10.1242/jeb.00233. [DOI] [PubMed] [Google Scholar]
  35. Gibbs AG, Rajpurohit S. In: Insect Hydrocarbons: Biology, Biochemistry, and Chemical Ecology. Gibbs AG, editor. Cambridge University Press; 2010. Cuticular lipids and water balance; pp. 100–120. [DOI] [Google Scholar]
  36. Goolsby EW, Bruggeman J, Ané C, Fitzjohn R. Rphylopars: fast multivariate phylogenetic comparative methods for missing data and within‐species variation. Methods in Ecology and Evolution. 2017;8:22–27. doi: 10.1111/2041-210X.12612. [DOI] [Google Scholar]
  37. Grafen A. The phylogenetic regression. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences. 1989;326:119–157. doi: 10.1098/rstb.1989.0106. [DOI] [PubMed] [Google Scholar]
  38. Hadley NF. Epicuticular lipids of the desert tenebrionid beetle, Eleodes armata: seasonal and acclimatory effects on composition. Insect Biochemistry. 1977;7:277–283. doi: 10.1016/0020-1790(77)90025-7. [DOI] [Google Scholar]
  39. Hadley NF. Water Relations of Terrestrial Arthropods. Academic Press; 1994. [Google Scholar]
  40. Heifetz Y, Miloslavski I, Aizenshtat Z, Applebaum SW. Cuticular surface hydrocarbons of desert locust nymphs, Schistocerca gregaria, and their effect on phase behavior. Journal of Chemical Ecology. 1998;24:1033–1047. doi: 10.1023/A:1022302519373. [DOI] [Google Scholar]
  41. Hoffmann AA, Hallas RJ, Dean JA, Schiffer M. Low potential for climatic stress adaptation in a rainforest Drosophila species. Science. 2003;301:100–102. doi: 10.1126/science.1084296. [DOI] [PubMed] [Google Scholar]
  42. Hoffmann AA, Weeks AR. Climatic selection on genes and traits after a 100 year-old invasion: a critical look at the temperate-tropical clines in Drosophila melanogaster from eastern Australia. Genetica. 2007;129:133–147. doi: 10.1007/s10709-006-9010-z. [DOI] [PubMed] [Google Scholar]
  43. Hoffmann AA. Physiological climatic limits in Drosophila: patterns and implications. The Journal of Experimental Biology. 2010;213:870–880. doi: 10.1242/jeb.037630. [DOI] [PubMed] [Google Scholar]
  44. Holze H, Schrader L, Buellesbach J. Advances in deciphering the genetic basis of insect cuticular hydrocarbon biosynthesis and variation. Heredity. 2021;126:219–234. doi: 10.1038/s41437-020-00380-y. [DOI] [PMC free article] [PubMed] [Google Scholar]
  45. Hood WG, Tschinkel WR. Desiccation resistance in arboreal and terrestrial ants. Physiological Entomology. 1990;15:23–35. doi: 10.1111/j.1365-3032.1990.tb00489.x. [DOI] [Google Scholar]
  46. Howard RW, Jackson LL, Banse H, Blows MW. Cuticular hydrocarbons of Drosophila birchii and D. serrata: identification and role in mate choice in D. serrata. Journal of Chemical Ecology. 2003;29:961–976. doi: 10.1023/a:1022992002239. [DOI] [PubMed] [Google Scholar]
  47. Huang J, Yu H, Guan X, Wang G, Guo R. Accelerated dryland expansion under climate change. Nature Climate Change. 2016;6:166–171. doi: 10.1038/nclimate2837. [DOI] [Google Scholar]
  48. Jezovit JA, Levine JD, Schneider J. Phylogeny, environment and sexual communication across the Drosophila genus. The Journal of Experimental Biology. 2017;220:42–52. doi: 10.1242/jeb.143008. [DOI] [PubMed] [Google Scholar]
  49. Johnson RA, Gibbs AG. Effect of mating stage on water balance, cuticular hydrocarbons and metabolism in the desert harvester ant, Pogonomyrmex barbatus. Journal of Insect Physiology. 2004;50:943–953. doi: 10.1016/j.jinsphys.2004.07.006. [DOI] [PubMed] [Google Scholar]
  50. Kellermann V, Loeschcke V, Hoffmann AA, Kristensen TN, Fløjgaard C, David JR, Svenning JC, Overgaard J. Phylogenetic constraints in key functional traits behind species’ climate niches: patterns of desiccation and cold resistance across 95 Drosophila species. Evolution; International Journal of Organic Evolution. 2012;66:3377–3389. doi: 10.1111/j.1558-5646.2012.01685.x. [DOI] [PubMed] [Google Scholar]
  51. Kellermann V, Hoffmann AA, Overgaard J, Loeschcke V, Sgrò CM. Plasticity for desiccation tolerance across Drosophila species is affected by phylogeny and climate in complex ways. Proceedings. Biological Sciences. 2018;285:20180048. doi: 10.1098/rspb.2018.0048. [DOI] [PMC free article] [PubMed] [Google Scholar]
  52. 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–1464. doi: 10.1093/bioinformatics/btq166. [DOI] [PubMed] [Google Scholar]
  53. Khallaf MA, Cui R, Weißflog J, Erdogmus M, Svatoš A, Dweck HKM, Valenzano DR, Hansson BS, Knaden M. Large-scale characterization of sex pheromone communication systems in Drosophila. Nature Communications. 2021;12:4165. doi: 10.1038/s41467-021-24395-z. [DOI] [PMC free article] [PubMed] [Google Scholar]
  54. Koto A, Motoyama N, Tahara H, McGregor S, Moriyama M, Okabe T, Miura M, Keller L. Oxytocin/vasopressin-like peptide inotocin regulates cuticular hydrocarbon synthesis and water balancing in ants. PNAS. 2019;116:5597–5606. doi: 10.1073/pnas.1817788116. [DOI] [PMC free article] [PubMed] [Google Scholar]
  55. Krupp JJ, Kent C, Billeter JC, Azanchi R, So AKC, Schonfeld JA, Smith BP, Lucas C, Levine JD. Social experience modifies pheromone expression and mating behavior in male Drosophila melanogaster. Current Biology. 2008;18:1373–1383. doi: 10.1016/j.cub.2008.07.089. [DOI] [PubMed] [Google Scholar]
  56. Krupp JJ, Nayal K, Wong A, Millar JG, Levine JD. Desiccation resistance is an adaptive life-history trait dependent upon cuticular hydrocarbons, and influenced by mating status and temperature in D. melanogaster. Journal of Insect Physiology. 2020;121:103990. doi: 10.1016/j.jinsphys.2019.103990. [DOI] [PubMed] [Google Scholar]
  57. Kühsel S, Brückner A, Schmelzle S, Heethoff M, Blüthgen N. Surface area-volume ratios in insects. Insect Science. 2017;24:829–841. doi: 10.1111/1744-7917.12362. [DOI] [PubMed] [Google Scholar]
  58. Kwan L, Rundle HD. Adaptation to desiccation fails to generate pre- and postmating isolation in replicate Drosophila melanogaster laboratory populations. Evolution; International Journal of Organic Evolution. 2010;64:710–723. doi: 10.1111/j.1558-5646.2009.00864.x. [DOI] [PubMed] [Google Scholar]
  59. Lamb AM, Wang Z, Simmer P, Chung H, Wittkopp PJ. Ebony affects pigmentation divergence and cuticular hydrocarbons in Drosophila americana and D. novamexicana. Frontiers in Ecology and Evolution. 2020;8:184. doi: 10.3389/fevo.2020.00184. [DOI] [PMC free article] [PubMed] [Google Scholar]
  60. Leeson SA, Kennington WJ, Evans TA, Simmons LW. Phenotypic plasticity but no adaptive divergence in cuticular hydrocarbons and desiccation resistance among translocated populations of dung beetles. Evolutionary Ecology. 2020;34:929–944. doi: 10.1007/s10682-020-10074-w. [DOI] [Google Scholar]
  61. Li F, Rane RV, Luria V, Xiong Z, Chen J, Li Z, Catullo RA, Griffin PC, Schiffer M, Pearce S, Lee SF, McElroy K, Stocker A, Shirriffs J, Cockerell F, Coppin C, Sgrò CM, Karger A, Cain JW, Weber JA, Santpere G, Kirschner MW, Hoffmann AA, Oakeshott JG, Zhang G. Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation. Molecular Ecology Resources. 2022;22:1559–1581. doi: 10.1111/1755-0998.13561. [DOI] [PMC free article] [PubMed] [Google Scholar]
  62. Liaw A, Wiener M. Classification and regression by randomforest. R News. 2002;2:18–22. [Google Scholar]
  63. Lynch M. Methods for the analysis of comparative data in evolutionary biology. Evolution; International Journal of Organic Evolution. 1991;45:1065–1080. doi: 10.1111/j.1558-5646.1991.tb04375.x. [DOI] [PubMed] [Google Scholar]
  64. Matzkin LM, Watts TD, Markow TA. Evolution of stress resistance in Drosophila: interspecific variation in tolerance to desiccation and starvation. Functional Ecology. 2009;23:521–527. doi: 10.1111/j.1365-2435.2008.01533.x. [DOI] [Google Scholar]
  65. Menzel F, Blaimer BB, Schmitt T. How do cuticular hydrocarbons evolve? physiological constraints and climatic and biotic selection pressures act on a complex functional trait. Proceedings. Biological Sciences. 2017;284:20161727. doi: 10.1098/rspb.2016.1727. [DOI] [PMC free article] [PubMed] [Google Scholar]
  66. Menzel F, Morsbach S, Martens JH, Räder P, Hadjaje S, Poizat M, Abou B. Communication versus waterproofing: the physics of insect cuticular hydrocarbons. The Journal of Experimental Biology. 2019;222:jeb210807. doi: 10.1242/jeb.210807. [DOI] [PubMed] [Google Scholar]
  67. Mueller LD. Density-dependent population growth and natural selection in food-limited environments: the Drosophila model. The American Naturalist. 1988;132:786–809. doi: 10.1086/284890. [DOI] [Google Scholar]
  68. Mundry R. In: Modern Phylogenetic Comparative Methods and Their Application in Evolutionary Biology. Mundry R, editor. Springer; 2014. Statistical issues and assumptions of phylogenetic generalized least squares; pp. 131–153. [DOI] [Google Scholar]
  69. Nunes TM, Oldroyd BP, Elias LG, Mateus S, Turatti IC, Lopes NP. Evolution of queen cuticular hydrocarbons and worker reproduction in stingless bees. Nature Ecology & Evolution. 2017;1:1–3. doi: 10.1038/s41559-017-0185. [DOI] [Google Scholar]
  70. Paradis E, Schliep K. Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2019;35:526–528. doi: 10.1093/bioinformatics/bty633. [DOI] [PubMed] [Google Scholar]
  71. Qiu Y, Tittiger C, Wicker-Thomas C, Le Goff G, Young S, Wajnberg E, Fricaux T, Taquet N, Blomquist GJ, Feyereisen R. An insect-specific P450 oxidative decarbonylase for cuticular hydrocarbon biosynthesis. PNAS. 2012;109:14858–14863. doi: 10.1073/pnas.1208650109. [DOI] [PMC free article] [PubMed] [Google Scholar]
  72. Rajpurohit S, Vrkoslav V, Hanus R, Gibbs AG, Cvačka J, Schmidt PS. Post-eclosion temperature effects on insect cuticular hydrocarbon profiles. Ecology and Evolution. 2021;11:352–364. doi: 10.1002/ece3.7050. [DOI] [PMC free article] [PubMed] [Google Scholar]
  73. Rane RV, Pearce SL, Li F, Coppin C, Schiffer M, Shirriffs J, Sgrò CM, Griffin PC, Zhang G, Lee SF, Hoffmann AA, Oakeshott JG. Genomic changes associated with adaptation to arid environments in cactophilic Drosophila species. BMC Genomics. 2019;20:52. doi: 10.1186/s12864-018-5413-3. [DOI] [PMC free article] [PubMed] [Google Scholar]
  74. Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things) Methods in Ecology and Evolution. 2012;3:217–223. doi: 10.1111/j.2041-210X.2011.00169.x. [DOI] [Google Scholar]
  75. Rouault JD, Marican C, Wicker-Thomas C, Jallon JM. In: Drosophila Melanogaster, Drosophila Simulans: So Similar, so Different. Capy P, editor. Springer; 2004. Relations between cuticular hydrocarbon (HC) polymorphism, resistance against desiccation and breeding temperature; a model for HC evolution in D; pp. 195–212. [DOI] [PubMed] [Google Scholar]
  76. Rusuwa BB, Chung H, Allen SL, Frentiu FD, Chenoweth SF. Natural variation at a single gene generates sexual antagonism across fitness components in Drosophila. Current Biology. 2022;32:3161–3169. doi: 10.1016/j.cub.2022.05.038. [DOI] [PubMed] [Google Scholar]
  77. Savage BA, Wang Z, Chung H, Masten S, Grieshop M. An ozonolysis based method and applications for the non-lethal modification of insect cuticular hydrocarbons. Journal of Chemical Ecology. 2021;47:628–641. doi: 10.1007/s10886-021-01285-8. [DOI] [PubMed] [Google Scholar]
  78. Sherwood S, Fu Q. A drier future? Science. 2014;343:737–739. doi: 10.1126/science.1247620. [DOI] [PubMed] [Google Scholar]
  79. Shi H, Tian H, Lange S, Yang J, Pan S, Fu B, Reyer CPO. Terrestrial biodiversity threatened by increasing global aridity velocity under high-level warming. PNAS. 2021;118:e2015552118. doi: 10.1073/pnas.2015552118. [DOI] [PMC free article] [PubMed] [Google Scholar]
  80. Soroker V, Hefetz A. Hydrocarbon site of synthesis and circulation in the desert ant Cataglyphis niger. Journal of Insect Physiology. 2000;46:1097–1102. doi: 10.1016/s0022-1910(99)00219-x. [DOI] [PubMed] [Google Scholar]
  81. Spikes AE, Paschen MA, Millar JG, Moreira JA, Hamel PB, Schiff NM, Ginzel MD. First contact pheromone identified for a longhorned beetle (Coleoptera: Cerambycidae) in the subfamily Prioninae. Journal of Chemical Ecology. 2010;36:943–954. doi: 10.1007/s10886-010-9837-8. [DOI] [PubMed] [Google Scholar]
  82. Svetnik V, Liaw A, Tong C, Culberson JC, Sheridan RP, Feuston BP. Random forest: a classification and regression tool for compound classification and QSAR modeling. Journal of Chemical Information and Computer Sciences. 2003;43:1947–1958. doi: 10.1021/ci034160g. [DOI] [PubMed] [Google Scholar]
  83. Van Oystaeyen A, Oliveira RC, Holman L, van Zweden JS, Romero C, Oi CA, d’Ettorre P, Khalesi M, Billen J, Wäckers F, Millar JG, Wenseleers T. Conserved class of Queen pheromones stops social insect workers from reproducing. Science. 2014;343:287–290. doi: 10.1126/science.1244899. [DOI] [PubMed] [Google Scholar]
  84. Wang L, Han X, Mehren J, Hiroi M, Billeter JC, Miyamoto T, Amrein H, Levine JD, Anderson DJ. Hierarchical chemosensory regulation of male-male social interactions in Drosophila. Nature Neuroscience. 2011;14:757–762. doi: 10.1038/nn.2800. [DOI] [PMC free article] [PubMed] [Google Scholar]
  85. Wang Y, Ferveur JF, Moussian B. Eco-genetics of desiccation resistance in Drosophila. Biological Reviews of the Cambridge Philosophical Society. 2021;96:1421–1440. doi: 10.1111/brv.12709. [DOI] [PubMed] [Google Scholar]
  86. Wigglesworth VB. Transpiration through the cuticle of insects. Journal of Experimental Biology. 1945;21:97–114. doi: 10.1242/jeb.21.3-4.97. [DOI] [Google Scholar]
  87. Wright MN, Ziegler A. Ranger: A Fast Implementation of Random Forests for High Dimensional Data in C++ and R. arXiv. 2015 doi: 10.1093/bioinformatics/btq257. https://arxiv.org/abs/1508.04409 [DOI] [PMC free article] [PubMed]

Editor's evaluation

Lynn M Riddiford 1

These studies have presented convincing evidence that dessication resistance in Drosophila species is conferred primarily by long, methyl-branched cuticular hydrocarbons. These fundamental findings add significantly to our understanding of how Drosophila species have evolved to adapt and survive in different environments. This study will be of interest to insect physiologists and ecologists as well as evolutionary biologists.

Decision letter

Editor: Lynn M Riddiford1

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Desiccation resistance differences in Drosophila species can be largely explained by variations in cuticular hydrocarbons" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and K VijayRaghavan as the Senior Editor. The reviewers have opted to remain anonymous.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) Perform control coating experiments and "replacement" experiments as outlined by Reviewer 1 in order to validate your hypothesis that the longer methyl-branched CHCs are critical to desiccation resistance.

2) Discuss the differences seen between males and females and the role that the methyl-branched cuticular pheromones may be playing.

3) Discuss the implications of using lab-reared flies for this analysis rather than wild-caught flies for the overall interpretation of the data.

4) Add a short discussion of how this study on Drosophila species may relate to desiccation resistance in other insects – mention similarities and differences.

Reviewer #1 (Recommendations for the authors):

A central issue in this work is to lift the correlative nature of the impressive data to the level of causativity. The coating experiment performed is, to me, insufficient to do so. I propose "relacement" and control experiments to support this point.

Despite this, there is another crucial point that needs to be addressed more convincingly: there are quite a few pieces of evidence published that mbCHC may be a major type of CHCs that is involved in desiccation resistance. So, the findings here are not really new. Admittedly, the amount of data/species is impressive; nevertheless, the statstical data (relatively low r-values) are not doubtlessly clarifying this matter. Indeed, this is a complex problem (trade-off of the CHC composition), and there won't be any simple solution to it. I will happy to discuss this with the authors.

Reviewer #2 (Recommendations for the authors):

This was a well-organized manuscript, and the authors do not go beyond their data in reaching their conclusions.

Reviewer #3 (Recommendations for the authors):

I have several suggestions for the authors as indicated below about the general scope of the study and interpretation of how this study could relate to other insects. The way this paper is written, it is most relevant to studies on Drosophila and may not pertain to other insects.

An aspect that is not discussed is the fact that many insects have internally branched methyl groups. It is surprising to me that these Drosophila species only have 2 methyl branched hydrocarbons. Although 2meCHC does have the highest melting temperature. Internally branched methyl groups have lower melting temperatures and would thus not contribute as much to desiccation. The discussion could include the contribution of other lipids to desiccation rates in other insects.

Line 66 – The authors describe the general mechanism in conserving water and preventing desiccation as a lipid layer found on the epicuticle and they named it cuticular hydrocarbons. This is not accurate because there are other types of lipids that could be found on the surface in addition to hydrocarbons. In some insects the majority of the lipids are hydrocarbons but other insects have few hydrocarbons but other types of lipids. The term cuticular lipids refers to a broader category of lipids and not just to hydrocarbons.

Line 121 – In addition to other types of lipids that could be found on the cuticular surface not all insects would have methyl branched hydrocarbons. Therefore this might not be a general mechanism to evolve higher levels of desiccation resistance across Insecta.

Line 258 – Compare D. immigrans vs. D. mojavensis. Both apparently have similar quantities of 2meC32 yet D. immigrans is not very desiccation resistant. Indicating that other factors are involved and is not due to just 2meC32.

Line 375 – Were body weight measurements take on dead or alive flies?

Materials and methods – Statistics

I cannot comment on the statistics as I am not familiar with most of the statistical packages used in the study.

eLife. 2022 Dec 6;11:e80859. doi: 10.7554/eLife.80859.sa2

Author response


Essential revisions:

1) Perform control coating experiments and "replacement" experiments as outlined by Reviewer 1 in order to validate your hypothesis that the longer methyl-branched CHCs are critical to desiccation resistance.

Reviewer 1 has suggested two different sets of experiments to be performed on flies without CHCs (CHC- flies). First, “control” experiments which coats “CHCs with the same chain lengths but without branched methyl groups” i.e., n-alkanes of different lengths on CHC- flies. Second, “replacement” experiments which coats mbCHCs of different lengths on CHC- flies. To address the reviewer’s comment, we performed the “control” experiments on both WT and CHC- flies and the “replacement” experiments on the CHC- flies. In summary, we performed the following experiments and these new data are added in a new Figure 5—figure supplement 2:

Initial submission

[WT +mbCHCs]

Resubmission

[WT +mbCHCs from initial submission], [WT +n-alkanes], [CHC- +mbCHCs], [CHC- +n-alkanes]

a) In our original submission, we showed that the coating of longer mbCHCs, 2MeC28 and 2MeC30, on wildtype (WT) D. melanogaster flies can increase desiccation resistance, while the shorter 2MeC26 does not have a significant effect (Figure 5E). We now provide additional analyses testing the correlation between mbCHC length and desiccation resistance and show that the significant increase in desiccation resistance is positively correlated with longer mbCHCs (Figure 5—figure supplement 2), suggesting that longer mbCHCs are able to confer higher desiccation resistance in WT flies. We obtained similar results with the “control” experiment, showing that the coating WT D. melanogaster flies with alkanes (C23, C25, C27, C29, C31) is able to increase desiccation resistance which is positively correlated with carbon chain length (Figure 5—figure supplement 2, Figure 5—figure supplement 3). This is not surprising as alkanes are saturated CHCs with high melting temperatures.

b) In the new “replacement” experiments, we showed that the coating of each of the mbCHCs, 2MeC26, 2MeC28, and 2MeC30, is able to significantly increase desiccation resistance in CHC- flies of both sexes (Figure 5—figure supplement 5). We also showed that the coating of all of the alkanes (C23, C25, C27, C29, C31) is able to significantly increase desiccation resistance in female CHC- flies, but only C31 is able to significantly increase desiccation resistance in male CHC- flies (Figure 5—figure supplement 5). However, desiccation resistance of these CHC- flies are very low and the coating of these synthetic CHCs did not increase desiccation resistance by much, albeit statistically significantly. Only a weak positive correlation between longer mbCHCs and higher desiccation was observed in female CHC- flies. There is no correlation between mbCHC length and higher desiccation in male CHC- flies and no correlation between alkane length and higher desiccation in female and male CHC- flies (Figure 5—figure supplement 6). This is similar with Krupp et al. 2020 where the authors showed coating CHC- with either only alkanes or mbCHCs only produce small increases in desiccation resistance, while coating CHC- flies with full WT extracts is able to result in a much larger increase in desiccation resistance. Both Krupp et al. and our observations are consistent with previous hypotheses that the CHC layer work as a blend rather than as individual components (Gibbs, 1998; Menzel et al., 2019; Wigglesworth, 1945).

Although we showed longer chain n-alkanes are also positively correlated with higher desiccation resistance when perfumed on WT flies, these CHCs are not widely found in Drosophila beyond the melanogaster group, thus were not ranked highly in our random forest analysis. In other words, we again suggest that longer chain n-alkanes are able to confer higher desiccation resistance, but do not contribute to the evolution of higher desiccation resistance in Drosophila. Therefore, we reiterate the conclusion that longer mbCHCs are critical for desiccation resistance across Drosophila but will also add the caveat that mbCHCs alone (in the absence of all the other CHCs in the profile of a Drosophila fly) do not confer high desiccation resistance.

We have modified added new results manuscript (Lines 222 to 250). In addition, we have modified the Materials and methods to include the new experiments (Lines 410 to 419) that uses a slightly different but highly effective method in generating CHC- flies.

References

Gibbs, A.G. (1998). Water-proofing properties of cuticular lipids. American Zoologist 38, 471-482.

Krupp, J.J., Nayal, K., Wong, A., Millar, J.G., and Levine, J.D. (2020). Desiccation resistance is an adaptive life-history trait dependent upon cuticular hydrocarbons, and influenced by mating status and temperature in D. melanogaster. Journal of Insect Physiology 121, 103990.

Menzel, F., Morsbach, S., Martens, J.H., Räder, P., Hadjaje, S., Poizat, M., and Abou, B. (2019). Communication versus waterproofing: the physics of insect cuticular hydrocarbons. Journal of Experimental Biology 222.

Wigglesworth, V.B. (1945). Transpiration through the cuticle of insects. Journal of Experimental Biology 21, 97-114.

2) Discuss the differences seen between males and females and the role that the methyl-branched cuticular pheromones may be playing.

To address this comment, we have split the comment into two parts: (a) discuss the differences seen between males and females and (b) discuss the role that the methyl-branched cuticular pheromones may be playing.

(a) We have added the following section to the Discussion section (Lines 297 to 310) of our manuscript to discuss the differences seen between males and females:

“Our overall analyses support the conclusion that the evolution of longer mbCHCs is correlated with higher desiccation resistance in our dataset, but we observed minor differences between males and females (Figures 3 and 5). In particular, while higher quantities of 2MeC28 and 2MeC32 are both positively correlated with higher desiccation resistance in both males and females (Figure 5B, D), there are differences between the sexes regarding 2MeC26 and 2MeC30 (Figure 5A, C). There could be several underlying reasons for this. Firstly, across our dataset, CHCs are sexually dimorphic in many species, leading to overall differences in the correlation analyses between the sexes. Secondly, other physiological differences, such as size differences, may contribute to the differences between males and females. In many species such as D. melanogaster, females are larger than males and this may contribute to the higher basal desiccation resistance (when CHCs are taken out of the equation) of female CHC- flies (3.6 ± 0.1 h) compared to that of male CHC- flies (2.1 ± 0.1 h) (Figure 5E). While these are minor differences in the overall correlation analyses between males and females, the synthetic CHC coating experiments showed that coating of longer mbCHCs confers higher desiccation resistance in both males and females.”

(b) We have expanded the discussion regarding sexual antagonism in D. serrata to discuss the roles that methyl-branched cuticular pheromones could be playing in affecting the differences in the evolution of mbCHCs and desiccation resistance (Lines 359 to 367):

“In D. serrata, having higher amounts of 2MeC26 leads to higher mating success in male D. serrata (Chung et al., 2014). Northern populations of this species produce lower amounts of 2MeC26 but higher amounts 2MeC30, due to a polymorphism in a fatty acyl-CoA reductase gene that shifts to an overall production of longer CHCs. As CHCs are sexually monomorphic in D. serrata, this polymorphism led to lower mating success in male northern flies and higher desiccation resistance in female northern flies compared to their common counterparts (Rusuwa et al., 2022). The sexually antagonistic effect due to the functions of the different mbCHCs may affect their evolutionary response to desiccation in many species.”

References

Chung, H., Loehlin, D.W., Dufour, H.D., Vaccarro, K., Millar, J.G., and Carroll, S.B. (2014). A single gene affects both ecological divergence and mate choice in Drosophila. Science 343, 1148-1151.

Rusuwa, B.B., Chung, H., Allen, S.L., Frentiu, F.D., and Chenoweth, S.F. (2022). Natural variation at a single gene generates sexual antagonism across fitness components in Drosophila. Current Biology 32, 3161-3169.

(3) Discuss the implications of using lab-reared flies for this analysis rather than wild-caught flies for the overall interpretation of the data.

We agree with the reviewer’s comment that the CHC composition found in lab-reared flies may not be reflective of CHC variation of individual species in the field. In published studies, almost all the CHC variation between different populations of a single species are quantitative differences, rather than qualitative differences in CHCs. In addition, as our study focused on a broad sampling of Drosophila and related species from different habitats spanning millions of years of evolution, rather than individual species, the overall interpretation of the data and the broad conclusion that we reached regarding the correlation between the evolution of mbCHCs and desiccation resistance is likely to extend to natural populations or wild-caught flies of the Drosophila genus. We have acknowledged the reviewer’s comment in the paper and added discussion points regarding this in the text (Lines 369 to 387):

“How much of what we learnt in this study is applicable broadly to natural populations of Drosophila and other insect species? The use of fly strains from the stock center reared on a common cornmeal artificial medium and the same abiotic factors rather than wild-caught flies allowed us to standardize this large-scale experiment. Previous studies have shown that diet (Fedina et al., 2012), temperature (Rajpurohit et al., 2021), and mating status (Everaerts et al., 2010) can affect CHC profiles in D. melanogaster, therefore the CHC profiles that we obtained in our studies may not be representative of the CHC variations present in some of these species in their natural habitats. However, these differences in CHC profiles are mostly quantitative, i.e., different quantities of the same CHCs, rather than different CHCs. Similarly, CHC variations between different natural populations of the same Drosophila species are also largely quantitative than qualitative (Dembeck et al., 2015; Rusuwa et al., 2022). Therefore, stock center strains would be representative of the qualitative CHC variations between species and our conclusions we reach in this study may be extended to wild populations. The use of standardized conditions also preclude the selection of species that require specialized media and conditions (e.g., species from the Hawaiian radiation and mycophagus species). While these are some of the limitations of our study, the standardized conditions allow us to reduce the number of confounds in our analyses for the correlation between CHC variation and desiccation resistance. As our study focused on a broad sampling of Drosophila and related species from different habitats spanning millions of years of evolution, the overall interpretation of the data and the broad conclusion that we reached regarding the correlation between the evolution of mbCHCs and desiccation resistance is likely to extend to natural populations of the Drosophila genus.”

(4) Add a short discussion of how this study on Drosophila species may relate to desiccation resistance in other insects – mention similarities and differences.

We have added a new section (Line 388 to 404) at the end of the discussion to discuss how our study may relate to desiccation resistance in other insects including some differences and similarities:

“Constraints in the evolution of CHCs due to their roles in many insect species as contact pheromones are likely to affect their evolution in response to desiccation resistance. In this study, we showed that the evolution of mbCHCs are major determinants of desiccation resistance across Drosophila species, as n-alkanes shared the same branch in the fatty acyl-CoA pathway as monoenes and dienes (Figure S8), many of which function as contact pheromones in Drosophila (Ferveur, 2005). In other insect orders such as Lepidoptera, where volatile non-CHC pheromones are predominantly used for chemical communication (Ando et al., 2004), constraints on each class of CHCs may be different from Drosophila. Likewise, across eusocial hymenopterans (Van Oystaeyen et al., 2014), saturated CHCs such as n-alkanes and mbCHCs are used as queen pheromones and are likely to be under constraints in the evolution of desiccation resistance in these species. Finally, in some rare examples such as dragonflies where very long-chain methyl ketones and aldehydes are the major cuticular wax components (Futahashi et al., 2019), the evolution of CHCs may not be a major factor in the evolution of desiccation resistance. While we limit our conclusions that the synthesis of mbCHCs with longer carbon-chain lengths is a common mechanism underlying the evolution of higher desiccation resistance to Drosophila, the evolution of cuticular components that modulate the melting temperature of the cuticular lipid layer can reduce water loss and increase desiccation resistance. We suggest that the evolution of CHC components in insects can be largely predictive of desiccation resistance between closely related species.”

Reviewer #1 (Recommendations for the authors):

A central issue in this work is to lift the correlative nature of the impressive data to the level of causativity. The coating experiment performed is, to me, insufficient to do so. I propose "relacement" and control experiments to support this point.

Despite this, there is another crucial point that needs to be addressed more convincingly: there are quite a few pieces of evidence published that mbCHC may be a major type of CHCs that is involved in desiccation resistance. So, the findings here are not really new. Admittedly, the amount of data/species is impressive; nevertheless, the statstical data (relatively low r-values) are not doubtlessly clarifying this matter. Indeed, this is a complex problem (trade-off of the CHC composition), and there won't be any simple solution to it. I will happy to discuss this with the authors.

For many of the data in the Figures 3 and 5, we use the Pearson correlation test to identify correlations between our variables. The validity of the Pearson correlation test does not rely on the linear assumptions. When the relationship is close to be linear, the Pearson correlation test will have typically higher power than other nonparametric tests, i.e., it is capable of detecting possible monotone relationship close to a linear relationship with fewer samples. However, we agree with the reviewer that the slopes might look misleading. Therefore, we have removed the lines/slopes from the figures.

Reviewer #3 (Recommendations for the authors):

I have several suggestions for the authors as indicated below about the general scope of the study and interpretation of how this study could relate to other insects. The way this paper is written, it is most relevant to studies on Drosophila and may not pertain to other insects.

An aspect that is not discussed is the fact that many insects have internally branched methyl groups. It is surprising to me that these Drosophila species only have 2 methyl branched hydrocarbons. Although 2meCHC does have the highest melting temperature. Internally branched methyl groups have lower melting temperatures and would thus not contribute as much to desiccation. The discussion could include the contribution of other lipids to desiccation rates in other insects.

Line 66 – The authors describe the general mechanism in conserving water and preventing desiccation as a lipid layer found on the epicuticle and they named it cuticular hydrocarbons. This is not accurate because there are other types of lipids that could be found on the surface in addition to hydrocarbons. In some insects the majority of the lipids are hydrocarbons but other insects have few hydrocarbons but other types of lipids. The term cuticular lipids refers to a broader category of lipids and not just to hydrocarbons.

Line 121 – In addition to other types of lipids that could be found on the cuticular surface not all insects would have methyl branched hydrocarbons. Therefore this might not be a general mechanism to evolve higher levels of desiccation resistance across Insecta.

Line 258 – Compare D. immigrans vs. D. mojavensis. Both apparently have similar quantities of 2meC32 yet D. immigrans is not very desiccation resistant. Indicating that other factors are involved and is not due to just 2meC32.

The reviewer is correct. While our model showed that CHC composition was able to explain 85.5% of the variation in desiccation resistance, it does not attribute all the variation in desiccation resistance to CHC composition. We agree suggest there are other physiological factors that also can contribute to desiccation resistance across the Drosophila genus.

Line 375 – Were body weight measurements take on dead or alive flies?

The body weight measurements were taken on alive flies.

Associated Data

    This section collects any data citations, data availability statements, or supplementary materials included in this article.

    Supplementary Materials

    Figure 5—source data 1. Desiccation resistance of Drosophila melanogaster attP40 flies when coated with individual mbCHC of 2MeC26, 2MeC28, and 2MeC30.
    Figure 5—source data 2. Desiccation resistance of Drosophila melanogaster attP40 flies when coated with individual n-alkane of C23, C25, C27, C29, and C31.
    Figure 5—source data 3. Desiccation resistance of Drosophila melanogaster CHC- flies when coated with individual CHC of 2MeC26, 2MeC28, 2MeC30, C23, C25, C27, C29, and C31.
    Figure 6—source data 1. Normalized quantities of mbCHCs, the length and normalized quantities of the longest mbCHC in each species.
    Supplementary file 1. Statistical output from phylogenetic comparative analyses and the list of species used in this study.

    (a) Phylogenetic signals for mbCHCs in Drosophila species. (b) Phylogenetic signals for desiccation resistance in Drosophila species. (c) Summary of the Phylogenetic Generalized Linear Square (PGLS) models between the longest mbCHCs and desiccation resistance for females and males in 50 Drosophila and related species. (d) List of species used in this study.

    elife-80859-supp1.docx (71.5KB, docx)
    MDAR checklist
    Source data 1. Data of normalized cuticular hydrocarbon (CHC) quantities, body weight, averaged desiccation resistance, total CHC quantity, and percentage of CHC quantities per body weight.
    elife-80859-data1.xlsx (168.7KB, xlsx)
    Source code 1. Code for Pearson’s correlation analysis and analysis of variance (ANOVA) tests in Figure 3 and Figure 5.
    elife-80859-code1.zip (23.4KB, zip)
    Source code 2. Code for NMDS plot and random forest analysis for Figure 4.
    elife-80859-code2.zip (8.4KB, zip)
    Source code 3. Code for Phylogenetic Generalized Linear Square (PGLS) analyses in Figure 6.
    elife-80859-code3.zip (5.4KB, zip)

    Data Availability Statement

    All data generated or analyzed during this study are included in the manuscript and supporting source data file. Code used is uploaded as source codes 1-3.


    Articles from eLife are provided here courtesy of eLife Sciences Publications, Ltd

    RESOURCES