Cuticular differences associated with aridity acclimation in African malaria vectors carrying alternative arrangements of inversion 2La
Parasites & Vectors volume 7, Article number: 176 (2014)
Principal malaria vectors in Africa, An. gambiae and An. coluzzii, share an inversion polymorphism on the left arm of chromosome 2 (2La/2L+a) that is distributed non-randomly in the environment. Genomic sequencing studies support the role of strong natural selection in maintaining steep clines in 2La inversion frequency along environmental gradients of aridity, and physiological studies have directly implicated 2La in heat and desiccation tolerance, but the precise genetic basis and the underlying behavioral and physiological mechanisms remain unknown. As the insect cuticle is the primary barrier to water loss, differences in cuticle thickness and/or epicuticular waterproofing associated with alternative 2La arrangements might help explain differences in desiccation resistance.
To test that hypothesis, two subcolonies of both An. gambiae and An. coluzzii were established that were fixed for alternative 2La arrangements (2La or 2L+a) on an otherwise homosequential and shared genetic background. Adult mosquitoes reared under controlled environmental conditions (benign or arid) for eight days post-eclosion were collected and analyzed. Measurements of cuticle thickness were made based on scanning electron microscopy, and cuticular hydrocarbon (CHC) composition was evaluated by gas chromatography–mass spectrometry.
After removing the allometric effects of body weight, differences in mean cuticle thickness were found between alternative 2La karyotypes, but not between alternative environments. Moreover, the thicker cuticle of the An. coluzzii 2La karyotype was contrary to the known higher rate of water loss of this karyotype relative to 2L+a. On the other hand, quantitative differences in individual CHCs and overall CHC profiles between alternative karyotypes and environmental conditions were consistent with expectation based on previous physiological studies.
Our results suggest that alternative arrangements of the 2La inversion are associated with differences in cuticle thickness and CHC composition, but that only CHC composition appears to be relevant for desiccation resistance. Differences in the CHC composition were consistent with previous findings of a lower rate of water loss for the 2L+a karyotype at eight days post-eclosion, suggesting that CHC composition is an important strategy for maintaining water balance in this genetic background, but not for 2La. Despite a higher rate of water loss at eight days, higher body water content of the 2La karyotype confers a level of desiccation resistance equivalent to that of the 2L+a karyotype.
A chromosomal inversion occurs when the chromosome breaks in two places and the intervening segment rotates 180° before rejoining. This results in a reversal of gene order without any change in gene content, aside from genes that may have been lesioned by the breakpoints themselves. Nevertheless, chromosomal inversions may have important phenotypic consequences, not least due to their effect on recombination. When evidence for chromosomal inversions was first documented in 1921, the potential for reduced recombination in individuals heterozygous for alternative gene arrangements was immediately hypothesized. It is now a tenet of classical genetics that a single crossover between inverted and uninverted arrangements within the inversion loop that forms during meiotic chromosomal pairing results in genetically unbalanced and inviable recombinant gametes, which are thus not recovered. It is this characteristic of effective suppression of recombination by inversions that has particularly interested evolutionary biologists, due to the potential for a chromosomal inversion event to capture between its breakpoints locally adapted alleles at multiple loci, and maintain their association against genetic exchange with deleterious immigrant alleles carried on an uninverted chromosome[2, 3]. In this way, inversions may facilitate adaptive evolution in heterogeneous environments[2, 4].
The notion that inversions are not neutral population genetic markers, but may be maintained by strong natural selection, arose from cytogenetic studies of Drosophila based on the banding pattern of giant larval salivary polytene chromosomes[5–7]. Species such as Drosophila pseudoobscura, D. subobscura, D. robusta and D. melanogaster are polymorphic for chromosomal inversions whose frequencies vary clinally along environmental gradients, cycle regularly with the seasons in concert with climatic variables, and shift in response to global warming[9, 10], suggesting that the adaptive benefit of alternative arrangements is dependent on local environmental conditions. One of the most compelling examples is inversion In(3R)P in D. melanogaster, which shows parallel clinal patterns along the eastern coasts of two different continents, from temperate regions where the inversion is rare, to tropical regions where its frequency is very high[11, 12]. Traits associated with inversion polymorphism in Drosophila include development time, body size, wing shape and thermal resistance, among others. Inversion polymorphism in a variety of other insects has also been associated with many of the same traits. Studies of inversion-trait associations in Anopheles mosquitoes have documented that wing shape in Anopheles funestus and resistance to heat, desiccation, and insecticides in An. gambiae and An. coluzzii[14–17] are linked to inversion polymorphisms.
Anopheles gambiae and An. coluzzii (formerly An. gambiae S and M molecular forms;) are principal vectors of human malaria in Africa, where 90% of the >600,000 deaths from malaria occur annually. Abundant inversion polymorphism in these species likely contributes to their successful exploitation of vast areas of sub-Saharan Africa, spanning diverse ecoclimatic zones from mesic rainforest to xeric savanna, but always in close association with humans[20, 21]. One of the oldest, most geographically widespread, and best-studied polymorphic inversion systems in these mosquitoes, 2La and its alternative 2L+a, is shared by common ancestry between An. gambiae and An. coluzzii. Frequencies of 2La are correlated with degree of aridity. This has been documented in West and Central Africa along environmental gradients of aridity from the mesic rainforest belt flanking the Gulf of Guinea where 2La is rare or absent, to xeric sahel savanna regions hundreds of kilometers to the north where 2La reaches fixation[23, 24]. Longitudinal studies in localities where 2La is polymorphic have demonstrated that its frequency fluctuates in a predictable pattern tied to annual rainfall variation between rainy and dry seasons[25, 26]. Even at the scale of a village, 2La is non-randomly distributed among mosquitoes resting indoors versus outdoors at night, a pattern attributed to the higher saturation deficit found inside dwellings where 2La is overrepresented. DNA sequencing studies of natural An. gambiae populations sampled along a cline of aridity in Cameroon support the notion that spatially varying selection is maintaining the 2La inversion polymorphism[23, 28]. Among the candidate genic targets of selection in 2La are cuticle protein (CPR) genes, many of which occur in a large cluster near the distal breakpoint[23, 28].
The strong correlation between 2La frequency and degree of aridity suggests that alleles within the 2La inversion may confer resistance to thermal and/or desiccation stress, both of which are environmental stressors routinely encountered by insects. Recent physiological studies have confirmed the association of 2La with thermal tolerance of An. coluzzii larvae and desiccation resistance of both An. coluzzii and An. gambiae adults[15, 17]. In particular, adult females homozygous for the 2La inversion were significantly more resistant to desiccation (measured as time to death in dry air), due to lower rates of water loss in newly emerged adults, and higher total body water content at four days post-emergence. Interestingly, by eight days post-emergence, there was no significant difference in survival between the karyotypes under desiccation stress. At this time point, the higher initial water content of 2La karyotypes was balanced by their higher rate of water loss relative to 2L+a—a pattern confirmed in eight-day old females acclimated to a lower humidity regimen. Prior acclimation to non-lethal desiccation stress improved desiccation resistance for females of both karyotypes, and although subsequent resistance to lethal desiccation stress was comparable between alternative karyotypes at eight days post-emergence, the orientation of the 2La rearrangement was associated with opposing effects on the storage of energy reserves in acclimated females. While 2La inverted homokaryotypes boosted their mass-specific glycogen content and reduced lipid stores, the converse was true for 2L+a homokaryotypes. It is noteworthy that increased storage of glycogen, utilized by Drosophila under desiccation stress, has been associated with higher dehydration tolerance in flies[29–31], consistent with its proposed role as an osmolyte for water retention. Taken together, these data suggest that genes inside the 2La inversion influence water balance and desiccation resistance, but also that resistance to desiccation is a complex trait, whose precise physiological basis and causal genes remain to be elucidated.
One of the main barriers to water loss in insects is the cuticle[32–35], which may be modified to enhance desiccation resistance. Increases in the amount of surface lipids (mainly cuticular hydrocarbons, CHCs) or changes in their chemical composition resulting in increased chain length, linearity, and saturation are the main means of minimizing cuticular transpiration in insects[32, 36]. Another non-mutually exclusive mechanism that may further reduce cuticular permeability is thickening of the cuticle. Increased cuticle thickness is associated with pyrethroid resistance in An. funestus, and structural cuticular proteins implicated in pyrethroid resistance in two other anophelines, An. gambiae and An. stephensi, are found in the endocuticle where they may contribute to cuticle thickening. In the Colorado potato beetle, transcripts encoding cuticular proteins were highly induced not only by insecticide treatment of resistant beetles but also by arid conditions, suggesting that increased cuticle thickness may be a common adaptive response to both environmental stresses.
A possible association of the 2La inversion system with cuticular characteristics of An. gambiae and An. coluzzii has not been investigated before now. The objective of this study was to assess karyotype-dependent effects on cuticle thickness as well as amount and composition of surface lipids in adults maintained under controlled environmental conditions.
Origin, husbandry and experimental treatments of mosquito colonies
Sub-lines of An. gambiae and An. coluzzii homokaryotypic for alternative arrangements of 2La (i.e., 2La/a and 2L+a/+a) were derived from parental colonies that were polymorphic for 2La. Apart from the alternative arrangements of 2La, the two sub-lines of each species that were employed for this study share the same genetic background. The An. coluzzii sub-lines SUCAM-2La and SUCAM-2L+a were derived in 2008 from a parental colony (SUCAM) established in 2005, as described previously. Using a similar approach, the homokaryotypic An. gambiae sub-lines NDKO-2La;2R+b and NDKO-2L+a;2R+b were derived from a parental colony (NDKO) established in 2008 from mosquitoes collected in Ndokayo, Cameroon. From the parental NDKO colony carrying 2La/2L+a and 2Rb/2R+b inversion polymorphisms, four homokaryotypic sub-lines were derived by molecular karyotyping of virgin males and females, as follows. On the morning of adult emergence, virgin females and males were separated into individual vials. DNA extracted from one rear leg of each mosquito was used for molecular karyotyping of 2La and 2Rb, performed in separate PCR reactions following White et al. and Lobo et al.. The second day, successfully karyotyped males and females were sorted into four population cages according to karyotype and allowed to mate: (1) 2La/a;2Rb/b, (2) 2La/a;2R+b/+b, (3) 2 L+a/+a;2Rb/b, and (4) 2L+a/+a;2R+b/+b. Females were offered blood meals to support egg production, and the progeny of each population cage were reared to adults. A random subset (N = 20) of the adult female progeny from each of the four populations was karyotyped by polytene chromosome analysis for each of three successive generations, to ensure the absence of inversion loops (indicative of heterokaryotypes). Experiments were started at least three generations after derivation of the NDKO-2La/2R+b and NDKO-2L+a/2R+b homokaryotypic sub-lines used in this study.
For routine colony husbandry, and for controlled experimental conditions referred to as “benign”, adult mosquitoes were held in insectaries maintained at 27°C and 85% relative humidity with a 12 h:12 h light:dark cycle and 1 h crepuscular (dawn and dusk) transitions. Non-lethal desiccation stress, referred to as “arid” experimental conditions, was imposed on adult females held in an insectary chamber maintained at 30°C and 60% RH. Immature stages of all sub-lines and both treatment groups were reared under common benign insectary conditions. Eggs were placed in plastic trays (27 × 16 × 6.5 cm) containing 2 L of water purified by reverse osmosis. Larvae were reared at equal density (200 per pan) and fed 180 mg of 2:1 finely ground tropical fish pellet:bakers yeast, every other day until day five, and daily thereafter until pupation. Pupae were transferred to 0.2 m3 emergence cages. Within 24 h of emergence, mosquitoes were sexed and virgin females were transferred to cages where they were held for eight days, either under “benign” or “arid” conditions. Adult females had free access to a 10% corn syrup solution supplied through saturated cotton balls provided fresh daily.
Cuticle thickness and total body dry mass measurements
At eight days post-emergence, at least 14 adult females from eight treatment groups (2 species x 2 karyotypes x 2 environmental conditions) were cold euthanized. Both mesothoracic legs of freshly killed mosquitoes were removed and preserved in 80% ethanol until they could be processed for scanning electron microscopy (SEM). The remaining carcass was desiccated overnight in an oven set to 70°C, and subsequent dry weight was determined using a microbalance (±0.2 μg, Mettler-Toledo, Columbus, OH, USA). Cuticle thickness was estimated by SEM from a cut through the left mesothoracic leg, following Wood et al.. First, the left leg was brushed gently with a small paintbrush to remove the scales. The first tarsal segment was cut at the midpoint with a new platinum coated razor blade. Both resulting leg fragments were mounted cut side up on carbon tape attached to an aluminum stub mount, and coated with 4 μm of iridium using a Cressington 208 HR sputter coater (Cressington Scientific Instruments, Watford, UK) in conjunction with the Cressington MTM 20 thickness monitor. Images of the cut of both leg fragments were taken with a FEI-Magellan 400 FESEM (FEI, Hillsboro, OR, USA), and the highest quality image chosen for measurement. If a high quality image of the left leg could not be obtained, the right leg was processed in the same manner.
Digital micrographs were analyzed using Zeiss AxioVision software. First, the inner and outer cuticle boundaries were delineated. The first two cuticle thickness measurements consisted of the length between the inner and outer boundary near the tendon and its mirror on the opposite side. Two additional thickness measurements were taken at points along the circumference midway between the first two measurements, and this process was repeated three more times to yield a maximum of 32 measurements (Figure 1). Sections containing cuticle imperfections that would impair measurements were bypassed, and images yielding less than 16 reliable measurements were discarded. Cuticle thickness for each female mosquito was reported as the mean value of 16–32 thickness measurements around the cuticle circumference.
Cuticular hydrocarbon (CHC) characteristics
Samples from a total of four treatments (2 karyotypes of An. coluzzii SUCAM x 2 environmental conditions) with five replicates each were analyzed on a coupled gas chromatography–mass spectrometer (GC-MS). Each replicate consisted of 10 virgin females. Total CHCs were extracted from a pool of 10 virgin females for each replicate by immersing the mosquitoes for 10 min in 200 μl of glass-distilled hexane (Fischer, ≥98.5% purity) in a 3 ml Wheaton® high recovery NextGen™ V-Vial attached with PTFE lined septa caps (Wheaton, Millville, NJ). Prior to their removal, mosquitoes were rinsed with an additional 150 μl of hexane. Hexane extracts were then evaporated under a gentle stream of charcoal filtered air using a Syntech Stimulus Controller CS-55 (Syntech, Kirchzarten, Germany). CHC extracts were reconstituted with 20 μl of hexane containing 10 ng/μl of octadecane as an internal standard (IS). Octadecane has widely been used as an IS in similar CHC analyses, and it was also chosen because its retention time did not overlap with any of the hydrocarbons in our biological samples.
Chemical analyses were performed on a 7890A GC system (Agilent Technologies, Santa Clara, CA) coupled with a 5975C Agilent Technologies mass spectrophotometer (inert XL MSD with a triple-Axis Detector). One microliter of an extract was injected in a split-splitless injector operating under splitless mode at 250°C in an Agilent HP-5 capillary column (30 m, 0.32 mm ID, 0.25 μm phase thickness). Helium (Ultra High Purity 5.0 Grade; Airgas, USA) was used as the carrier gas at a constant flow rate of 1 ml/min. The column was held isothermally at 70°C for 1 minute, then programmed to increase at a rate of 10°C per minute until 310°C, with a final hold of 10 minutes. MS was operated at 70 eV. Data recording and quantification was performed using Agilent MSD ChemStation software (E.02.02.1431). Initial chemical identity was determined using NIST 2011 MS library.
Synthetic hydrocarbons of known carbon chain lengths from C10 to C34 were run to determine Kovat’s Retention Index (RI). There were sixteen constituent peaks that consistently appeared in the chromatograms of each pool, and these were selected for further analyses since their amounts were well above the detection limit of the MS (as inferred from the IS amounts). Identification of these CHC was made in two ways: 1) comparison of the RI calculated for each of the 16 compounds relative to known RIs from Pherobase (http://www.pherobase.com); and 2) confirmation of retention times and their mass spectra compared with synthetic standards.
For each replicate, the area under each CHC peak was quantified and expressed as a proportional value relative to the sum of all peaks. Relative abundance values were Gaussian-transformed to normalize the data prior to analysis.
Statistical analyses were conducted using R v2.13.1 and SAS 9.2 (SAS Institute, Cary, NC).
Initial exploration of the relationship between total body dry weight and cuticle thickness was based on the Pearson’s product moment correlation by species and environmental condition (four groups). An allometric relationship between body weight and cuticle thickness meant that variation in body weight among individuals in a group would confound inferences about the direct contribution of karyotype to cuticle thickness. Accordingly, the normalization approach described by Lleonart et al. to remove allometric effects of body size was adapted for the general linear modeling framework applied here. Specifically, the initial general linear models, by species and condition, contained log-transformed cuticle thickness as the response variable, and log-transformed total dry body weight, karyotype, and their interaction as the independent variables. As the interaction terms from all four initial models were not statistically significant (all P-values from the t-tests of interaction were > 0.20), they were removed from the final models. From each of the four models, the least-squares means and 95% confidence intervals (CI) were obtained for log-transformed cuticle thickness by karyotype, which were then back-transformed to the linear scale (estimates known as geometric means) for reporting purposes. Statistical significance of the difference in standardized cuticle thickness between karyotypes for each species and condition (four comparisons) was determined from P-values based on t-statistics.
Statistical significance of differences in the amounts of individual compounds was assessed using ANOVA models that included karyotype, environmental condition, and their interaction as main (fixed) factors. Quantitative differences in overall CHC profiles were analyzed by principal component analysis (PCA). The first four PCs with eigenvalues greater than 1 explained 88% of the variance. Multivariate analysis of variance (MANOVA) with these first four PCs as response variables, and karyotype, environmental condition, and their interaction as factors was conducted to test for statistical significance of differences due to karyotype and environment.
Using the same SUCAM colony of An. coluzzii, Gray et al. showed that at eight days post-emergence, females homokaryotypic for the 2La inversion lost water at a significantly higher rate (standardized for body mass) than their 2L+a counterparts, despite an otherwise common genetic background. Although acclimation to arid environmental conditions reduced the rate of water loss relative to non-acclimated females regardless of karyotype, the difference between karyotypes—whereby female carriers of 2La lost water at a significantly higher rate—was consistent with or without acclimation. Accordingly, to assess whether changes in cuticle characteristics could be implicated in karyotype-based differences in standardized rates of water loss, we employed the same sex, age and experimental conditions as in Gray et al..
Positive allometry of cuticle thickness and total body dry weight
Adult female mosquitoes maintained for eight days from the time of eclosion under benign or arid insectary conditions were evaluated in two species, An. coluzzii (SUCAM) and An. gambiae (NDKO), to assess whether the orientation of the 2La chromosomal rearrangement impacted cuticle thickness. Total body dry weight was recorded together with cuticle thickness for the 14–18 individual females included in each of eight treatment groups (2 species x 2 karyotypes x 2 environmental conditions), in anticipation that body size would co-vary with cuticle thickness. Mean measurements of cuticle thickness and total body dry weight for each of the eight groups are provided in Table 1. Assessed on its own, total body dry weight was significantly higher for females acclimated to arid conditions relative to those maintained under benign conditions, and this trend applied for both species regardless of karyotype (P < 0.0001 for each karyotype/species; Figure 2). Moreover, initial tests of the expected positive relationship between dry weight and cuticle thickness for adult females maintained under the same environmental condition revealed a significant correlation in both species (Table 2). This relationship is illustrated in Figure 3 by logarithmic scatter plots of cuticle thickness by total body dry mass for each of the eight treatment groups. The t-tests of interaction between log(dry weight) and karyotype for the four species and condition combinations yielded non-significant P-values in all cases (Figure 3A-D), indicating that for a given species and environmental condition, the slopes of the regression lines for alternative karyotypes were indistinguishable statistically.
Inversion-specific effects on cuticle thickness
After removing the allometric effect of body size (see Methods), the resulting least-squares geometric means of cuticle thickness were compared between karyotypes, for each species and environmental condition (Figure 4). For An. gambiae NDKO, mean cuticle thickness between karyotypes was not statistically significant owing to strongly overlapping 95% confidence intervals. In contrast, for An. coluzzii SUCAM there was a statistically significant difference in cuticle thickness between 2L+a and 2La under benign conditions, and a marginally non-significant difference (P = 0.08) between the karyotypes under arid conditions (Figure 4). In this case, mean cuticle thickness was greater for the 2La karyotype under both benign and arid conditions.
Inversion-dependent effects on cuticular hydrocarbon composition
Epicuticular hydrocarbon composition was examined in An. coluzzii (SUCAM), using eight day-old virgin females of both (2La or 2L+a) karyotypes, maintained under benign or arid conditions. Regardless of karyotype or environmental condition, the same 16 peaks were consistently observed in all chromatograms, suggesting that these two factors had no qualitative effect on CHC composition (Figure 5). Hydrocarbon identification, based on mass spectra and retention times relative to known standards, revealed several alkanes, alkenes, and methyl-branched alkanes with base chain lengths ranging from 23 to 33 carbons; one compound (peak 14; KI 3158) could not be identified (Table 3). Peaks 4, 7, and 13 (heptacosane, nonacosane, and 13-methylhentriacontane, respectively) were the most abundant, accounting for about 61 to 68% of the total CHC profile (Figures 5 and6).
Assessment of quantitative differences in the 16 CHCs between treatment groups was conducted two ways, by compound and across the CHC profile. Compound-by-compound ANOVA tests revealed six CHCs whose relative abundance was significantly different owing either to karyotype or environmental condition, and an additional four CHCs were responsive to both factors (Figure 6), but no significant interaction was detected between karyotype and environment. Of the seven compounds responsive to karyotype (peaks 1, 5, 6, 7, 9, 10, 11), all save one (7, nonacosane) were more abundant in the 2La background. Interestingly, the relative proportion of nonacosane was significantly lower, while the relative proportions of two alkenes (peak 6, 1-nonacosene; and peak 11, 1-hentriacontene) and one methyl-branched alkane (peak 10, methyltriacontane) were significantly higher in the 2La background. Lower amounts of an unbranched alkane in favor of more alkenes and a branched alkane should disrupt lipid packing and lower the melting point relative to 2L+a. Seven compounds were responsive to environmental condition. Of those, five were unbranched alkanes that increased in proportion in arid-acclimated females (peaks 5, 7, 9, 12, 15). Two compounds decreased in response to aridity, an alkene (peak 6, 1-nonacosene) and the unknown (peak 14). Given that double bonds and other chain modifications disrupt lipid packing, and that the unknown compound decreased under arid conditions, it is tempting to speculate that this compound may be an alkene or a branched alkane.
Principal component analysis (PCA) was used to compare overall CHC profiles. The first three PCs accounted for 77.5% of the total variance (Figure 7, Table 4). The factor loadings (the correlation between individual compounds and the PC) above 0.25 were used to identify the compounds most influential for each of the first three PCs (Table 4), to guide interpretation of the biological significance of each vector. PC1 was weighted negatively by peaks 2, 4 and 7 (pentacosane, heptacosane, nonacosane), and positively by peaks 10, 11, 13, and 16 (three methyl-branched alkanes and an alkene). Thus, PC1 may describe trade-offs between shorter chain alkanes that were relatively more abundant in the 2L+a background—especially under arid conditions—and longer CHCs with chain modifications that were more abundant in the 2La background, especially under benign conditions (Figure 6). That PC1 describes a karyotype effect, particularly for 2L+a under arid conditions and 2La under benign conditions, is reflected in a plot of the mean PC1 scores by karyotype and environmental condition (Figure 8A). PC2 was negatively weighted by peaks 1, 3, 5, 6, and 8–10, including an alkene and two methyl-branched alkanes; all of the corresponding compounds were more abundant in 2La karyotypes relative to their 2L+a counterparts (Figure 6). Accordingly, PC2 seems to be a vector that describes variation in CHC profiles in response to karyotype differences, as suggested by Figure 8B. PC3 was negatively weighted by peaks 7, 9, 12, and 15, corresponding to the longer unbranched alkanes that were proportionally more abundant in arid-acclimated females. PC3 also was negatively weighted by peak 14, the unknown compound that decreased in abundance in response to aridity (Figure 6). Therefore, PC3 is a vector that explains variation in CHC profiles due to environmental conditions (Figure 8C). Multivariate analysis using the PCs with eigenvalues greater than 1 as response variables, and karyotype, environment and their interaction as fixed factors, indicated that karyotype and environment, but not their interaction, were significant influences (Table 5).
Climate is a major factor limiting the geographic and seasonal ranges of insects. Climatic variables, particularly those related to water availability, are effective predictors of the distribution of African malaria vectors in the Anopheles gambiae complex and the chromosomal inversion polymorphisms that segregate within them[49–52]. Physiological studies of An. gambiae and An. coluzzii, based on field isolates as well as laboratory colonies, suggest that the superior desiccation resistance afforded by alleles linked to inversion 2La and other inversions on chromosome 2 are major determinants responsible for their relative abundance under arid conditions[15, 17, 53]. Both species can mount a plastic response to desiccation stress, as acclimation to arid conditions increases their resistance. However, the differential desiccation resistance of alternative 2La karyotypes in an otherwise common genetic background, under controlled and identical environmental conditions, indicates that physiological differences between karyotypes have a genetic basis. Genomic scans of elevated sequence divergence between karyotypes suggested a possible role for cuticle protein genes[23, 28]. This study aimed to test whether alternative 2La karyotypes were associated with differences in cuticle characteristics that could contribute to differential desiccation resistance.
Mean cuticle thickness, when standardized for body weight, differed between alternative karyotypes, although the difference was statistically significant only in An. coluzzii maintained under benign conditions. Moreover, in both species, environmental condition—whether benign or arid—did not alter the relationship between cuticle thickness and karyotype. Recently, Vannini et al. studied the location of the structural cuticular protein CPF3, whose transcripts are significantly more abundant in virgin females of An. coluzzii relative to An. gambiae. The location of this protein in appendage exocuticle was deemed inconsistent with a role for CPF3 in pheromone display for mate recognition, as proposed previously[54, 55]. However, Vannini et al. suggested that its exocuticular location and higher transcript abundance in An. coluzzii might be related to the greater desiccation resistance of this species. If so, CPF3 may contribute either to a more densely packed or a thicker epicuticular layer in An. coluzzii. The CPF3 gene is not located inside the 2La inversion, and there is no a priori expectation that CPF3 should be differentially expressed between 2La karyotypes of An. coluzzii, but comparisons between species (by karyotype and condition) suggest that if anything, An. coluzzii has the thinner cuticle of the two species.
Within An. coluzzii, mean cuticle thickness was greater for 2La than 2L+a females at eight days post-eclosion. If a thicker cuticle were a better barrier to water loss, the expectation would be that females with thicker cuticles should lose water at a lower rate. However, physiological studies showed the contrary; based on the same An. coluzzii SUCAM colonies, it was found that eight day old An. coluzzii 2La females actually lost water at a significantly higher rate than their 2L+a counterparts with thinner cuticles. Taken together, these data suggest that for these anophelines, increased cuticle thickness per se does not appear to reduce water loss and is unlikely to be an adaptation to desiccation stress, although this trait has been implicated in resistance of Anopheles funestus to insecticides. A caveat to this conclusion is that cuticle thickness was only measured from an appendage; we cannot rule out the possibility that cuticle thickness elsewhere in the body increased in response to aridity.
Cuticular hydrocarbons are the dominant constituents of the waxy epicuticular layer that coats the insect cuticle, thought to be the main barrier to water loss. We were able to detect 16 CHCs in An. coluzzii females in common with previous studies of An. gambiae and An. coluzzii field populations[56, 57], but were unable to detect the dimethyl-branched CHCs (chain length C33-C45) observed in those studies despite similar methods. Based on the 16 CHCs detected, we found no qualitative differences in CHC composition between karyotypes or environmental conditions. However, we did find quantitative differences that were consistent with expectations. Exposure to arid conditions caused a decrease in the relative proportions of unsaturated CHCs and those with methyl-branched chains in favor of unbranched, saturated CHCs in both karyotypes, an effect expected to enable denser packing of the lipid layer, hence a higher melting point and superior waterproofing. In terms of karyotype-associated CHC profiles, the 2L+a karyotype was characterized by a greater proportion of saturated and unbranched hydrocarbon chains relative to the 2La karyotype, apparently lending 2L+a females better waterproofing under desiccation stress. Although this result seems counter-intuitive given the association of 2La (and not 2L+a) with aridity in natural populations, it is entirely consistent with previous physiological studies demonstrating that An. coluzzii 2La females lost water at a higher rate than their 2L+a counterparts, at eight days post-eclosion. Those same physiological studies found that at the vulnerable teneral stage when the cuticle is still hardening, 2La females lost water at a significantly lower rate than 2L+a females, potentially explaining the advantage of the 2La karyotype (more precisely, the advantage of causal alleles linked to 2La) in arid environments. In fact, the water loss rate for 2La females was essentially constant from adult eclosion to eight days post-eclosion; the difference in water loss rates between karyotypes was a function of the reduction of the rate of water loss in 2L+a females during that time-frame.
The challenge for insects in desiccating environments is to maintain their water balance. One important component of the overall strategy of water regulation, considered in this study, is to conserve water by minimizing its loss via the route responsible for the greatest body water loss, the cuticle. However, there are other components of insect water economy that were not considered here, including physiological adjustments following a blood meal, that operate in concert to determine the ultimate desiccation resistance phenotype. During the off-host stage (sensu), mosquitoes may increase their water intake by imbibing more nectar or other plant juices. Additionally, they may limit water loss by resting in humid microenvironments, and avoiding activity during the relatively hot, dry daytime hours. Water loss may also be limited by increasing body size, thereby reducing the surface-to-volume ratio—a strategy that both An. gambiae and An. coluzzii seem to have adopted in our experiments based on their increased dry mass under arid conditions (Figure 2). Increasing tolerance to desiccation is another non-exclusive strategy, a process that may be aided by heat shock proteins such as Hsp70 and Hsp90[59, 60]. In this regard, it is intriguing that a tandem trio of Hsp90 genes are found inside the 2La inversion, in a region of exceptionally high sequence divergence from 2L+a that has been implicated in the selective maintenance of the inversion. Not least, desiccation resistance may be improved by boosting glycogen content, as this molecule may be used to store water. Indeed, this may be one of the primary ways in which alleles within the 2La inversion contribute to desiccation resistance, as Gray et al. found that mass-specific body water content was significantly higher at four and eight days post-eclosion in 2La versus 2L+a females, a circumstance that presumably compensated for the significantly higher rate of water loss observed in eight day-old 2La females.
There remains much to be learned of the underlying molecular mechanisms governing the complex and sophisticated phenotype that is desiccation resistance. Yet ability to withstand desiccation is a primary determinant of the geographic, microgeographic, and seasonal distribution of malaria vectors. The difficulty of detecting even the most desiccation resistant An. coluzzii mosquitoes in the West African Sahel strongly testifies to the limits that aridity poses to these vectors[61, 62]. The differential effects of alternative inversion karyotypes on mosquito behaviors such as endophily, which can be interpreted in terms of optimal habitat choice with respect to degree of aridity, also have important impacts on the likelihood of vector-human contact and vector contact with indoor-based vector control. Thus, they affect vectorial capacity and the epidemiology of malaria. Importantly, although the 2La inversion is associated with superior desiccation resistance and its frequencies peak in An. gambiae and An. coluzzii populations living in arid climates, the same inversion is rare or absent in populations from mesic environments, implying fitness tradeoffs. Deeper understanding of water regulation and its limits should uncover the hidden costs as well as the benefits of water balance strategies, information that could be exploited for vector control.
Sturtevant AH: A Case of Rearrangement of Genes in Drosophila. Proc Natl Acad Sci U S A. 1921, 7: 235-237. 10.1073/pnas.7.8.235.
Kirkpatrick M, Barton N: Chromosome inversions, local adaptation and speciation. Genetics. 2006, 173: 419-434. 10.1534/genetics.105.047985.
Hoffmann AA, Rieseberg LH: Revisiting the impact of inversions in evolution: from population genetic markers to drivers of adaptive shifts and speciation?. Ann Rev Ecol Evol Syst. 2008, 39: 21-42. 10.1146/annurev.ecolsys.39.110707.173532.
Schaeffer SW: Selection in heterogeneous environments maintains the gene arrangement polymorphism of Drosophila pseudoobscura. Evolution. 2008, 62: 3082-3099. 10.1111/j.1558-5646.2008.00504.x.
Dobzhansky T: Genetics of the Evolutionary Process. 1970, New York: Columbia University Press
Powell JR: Progress and Prospects in Evolutionary Biology: The Drosophila Model. 1997, Oxford: Oxford University Press
Krimbas CB, Powell JR: Drosophila inversion polymorphism. 1992, London: CRC Press
Rodriguez-Trelles F: Seasonal cycles of allozyme-by-chromosomal-inversion gametic disequilibrium in Drosophila subobscura. Evolution. 2003, 57: 839-848.
Etges WJ, Arbuckle KL, Levitan M: Long-term frequency shifts in the chromosomal polymorphisms of Drosophila robusta in the Great Smoky Mountains. Biol J Linn Soc. 2006, 88: 131-141. 10.1111/j.1095-8312.2006.00617.x.
Balanya J, Oller JM, Huey RB, Gilchrist GW, Serra L: Global genetic change tracks global climate warming in Drosophila subobscura. Science. 2006, 313: 1773-1775. 10.1126/science.1131002.
Singh RS, Hickey DA, David J: Genetic differentiation between geographically distant populations of Drosophila melanogaster. Genetics. 1982, 101: 235-256.
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. 10.1007/s10709-006-9010-z.
Ayala D, Caro-Riano H, Dujardin JP, Rahola N, Simard F, Fontenille D: Chromosomal and environmental determinants of morphometric variation in natural populations of the malaria vector Anopheles funestus in Cameroon. Infect Genet Evol. 2011, 11: 940-947. 10.1016/j.meegid.2011.03.003.
Rocca KA, Gray EM, Costantini C, Besansky NJ: 2La chromosomal inversion enhances thermal tolerance of Anopheles gambiae larvae. Malar J. 2009, 8: 147-10.1186/1475-2875-8-147.
Fouet C, Gray E, Besansky NJ, Costantini C: Adaptation to aridity in the malaria mosquito Anopheles gambiae: chromosomal inversion polymorphism and body size influence resistance to desiccation. PLoS ONE. 2012, 7: e34841-10.1371/journal.pone.0034841.
Brooke BD, Hunt RH, Chandre F, Carnevale P, Coetzee M: Stable chromosomal inversion polymorphisms and insecticide resistance in the malaria vector mosquito Anopheles gambiae (Diptera: Culicidae). J Med Entomol. 2002, 39: 568-573. 10.1603/0022-2585-39.4.568.
Gray EM, Rocca KA, Costantini C, Besansky NJ: Inversion 2La is associated with enhanced desiccation resistance in Anopheles gambiae. Malar J. 2009, 8: 215-10.1186/1475-2875-8-215.
Coetzee M, Hunt RH, Wilkerson R, Della Torre A, Coulibaly MB, Besansky NJ: Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013, 3619: 246-274.
World Health Organisation: World Malaria Report: 2013.http://www.who.int/malaria/publications/world_malaria_report_2013/en/,
Della Torre A, Tu Z, Petrarca V: On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms. Insect Biochem Mol Biol. 2005, 35: 755-769. 10.1016/j.ibmb.2005.02.006.
Pombi M, Caputo B, Simard F, Di Deco MA, Coluzzi M, Della Torre A, Costantini C, Besansky NJ, Petrarca V: Chromosomal plasticity and evolutionary potential in the malaria vector Anopheles gambiae sensu stricto: insights from three decades of rare paracentric inversions. BMC Evol Biol. 2008, 8: 309-10.1186/1471-2148-8-309.
White BJ, Cheng C, Sangare D, Lobo NF, Collins FH, Besansky NJ: The population genomics of trans-specific inversion polymorphisms in Anopheles gambiae. Genetics. 2009, 183: 275-288. 10.1534/genetics.109.105817.
Cheng C, White BJ, Kamdem C, Mockaitis K, Costantini C, Hahn MW, Besansky NJ: Ecological genomics of Anopheles gambiae along a latitudinal cline: a population-resequencing approach. Genetics. 2012, 190: 1417-1432. 10.1534/genetics.111.137794.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 73: 483-497. 10.1016/0035-9203(79)90036-1.
Petrarca V, Sabatinelli G, Di Deco MA, Papakay M: The Anopheles gambiae complex in the Federal Islamic Republic of Comoros (Indian Ocean): some cytogenetic and biometric data. Parassitologia. 1990, 32: 371-380.
Rishikesh N, Di Deco MA, Petrarca V, Coluzzi M: Seasonal variations in indoor resting Anopheles gambiae and Anopheles arabiensis in Kaduna, Nigeria. Acta Trop. 1985, 42: 165-170.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Behavioural divergences between mosquitoes with different inversion karyotypes in polymorphic populations of the Anopheles gambiae complex. Nature. 1977, 266: 832-833. 10.1038/266832a0.
White BJ, Hahn MW, Pombi M, Cassone BJ, Lobo NF, Simard F, Besansky NJ: Localization of candidate regions maintaining a common polymorphic inversion (2La) in Anopheles gambiae. PLoS Genet. 2007, 3: e217-10.1371/journal.pgen.0030217.
Marron MT, Markow TA, Kain KJ, Gibbs AG: Effects of starvation and desiccation on energy metabolism in desert and mesic Drosophila. J Insect Physiol. 2003, 49: 261-270. 10.1016/S0022-1910(02)00287-1.
Gibbs AG, Chippindale AK, Rose MR: Physiological mechanisms of evolved desiccation resistance in Drosophila melanogaster. J Exp Biol. 1997, 200: 1821-1832.
Parkash R, Ranga P: Sex-specific divergence for adaptations to dehydration stress in Drosophila kikkawai. J Exp Biol. 2013, 216: 3301-3313. 10.1242/jeb.087650.
Benoit JB: Water management by dormant insects: comparisons between dehydration resistance during summer aestivation and winter diapause. Prog Mol Subcell Biol. 2010, 49: 209-229. 10.1007/978-3-642-02421-4_10.
Chown SL, Sorensen JG, Terblanche JS: Water loss in insects: an environmental change perspective. J Insect Physiol. 2011, 57: 1070-1084. 10.1016/j.jinsphys.2011.05.004.
Hadley NF: Water relations of terrestrial arthropods. 1994, San Diego: Academic Press
Benoit JB, Denlinger DL: Meeting the challenges of on-host and off-host water balance in blood-feeding arthropods. J Insect Physiol. 2010, 56: 1366-1376. 10.1016/j.jinsphys.2010.02.014.
Gibbs AG, Rajpurohit S: Cuticular lipids and water balance. Insect Hydrocarbons: Biology, Biochemistry, and Chemical Ecology. Edited by: Blomquist GJ, Bagnères A-G. 2010, Cambridge: Cambridge University Press, 100-120.
Wood O, Hanrahan S, Coetzee M, Koekemoer L, Brooke B: Cuticle thickening associated with pyrethroid resistance in the major malaria vector Anopheles funestus. Parasit Vectors. 2010, 3: 67-10.1186/1756-3305-3-67.
Awolola TS, Oduola OA, Strode C, Koekemoer LL, Brooke B, Ranson H: Evidence of multiple pyrethroid resistance mechanisms in the malaria vector Anopheles gambiae sensu stricto from Nigeria. Trans R Soc Trop Med Hyg. 2009, 103: 1139-1145. 10.1016/j.trstmh.2008.08.021.
Vontas J, David JP, Nikou D, Hemingway J, Christophides GK, Louis C, Ranson H: Transcriptional analysis of insecticide resistance in Anopheles stephensi using cross-species microarray hybridization. Insect Mol Biol. 2007, 16: 315-324. 10.1111/j.1365-2583.2007.00728.x.
Vannini L, Reed TW, Willis JH: Temporal and spatial expression of cuticular proteins of Anopheles gambiae implicated in insecticide resistance or differentiation of M/S incipient species. Parasit Vectors. 2014, 7: 24-10.1186/1756-3305-7-24.
Zhang J, Goyer C, Pelletier Y: Environmental stresses induce the expression of putative glycine-rich insect cuticular protein genes in adult Leptinotarsa decemlineata (Say). Insect Mol Biol. 2008, 17: 209-216. 10.1111/j.1365-2583.2008.00796.x.
Lee CE, Frost BW: Morphological stasis in the Eurytemora affinis species complex (Copepoda: Temoridae). Hydrobiologia. 2002, 480: 111-128. 10.1023/A:1021293203512.
White BJ, Santolamazza F, Kamau L, Pombi M, Grushko O, Mouline K, Brengues C, Guelbeogo W, Coulibaly M, Kayondo JK, Sharakhov I, Simard F, Petrarca V, Della Torre A, Besansky NJ: Molecular karyotyping of the 2La inversion in Anopheles gambiae. Am J Trop Med Hyg. 2007, 76: 334-339.
Lobo NF, Sangare DM, Regier AA, Reidenbach KR, Bretz DA, Sharakhova MV, Emrich SJ, Traore SF, Costantini C, Besansky NJ, Collins FH: Breakpoint structure of the Anopheles gambiae 2Rb chromosomal inversion. Malar J. 2010, 9: 293-10.1186/1475-2875-9-293.
Della Torre A: Polytene chromosome preparation from anopheline mosquitoes. Molecular Biology of Disease Vectors: A Methods Manual. Edited by: Crampton & Hall JM, Beard CB, Louis C. 1997, London: Chapman, 329-336.
Lleonart J, Salat J, Torres GJ: Removing allometric effects of body size in morphological analysis. J Theor Biol. 2000, 205: 85-93. 10.1006/jtbi.2000.2043.
Tabachnick B, Fidell L: Using multivariate statistics. 1989, New York: Harper Collins
Sutherst R, Maywald G, Skarratt D: Predicting insect distributions in a changed climate. Insects in a changing environment. Edited by: Harrington R, Stork N. 1995, London: Academic Press, 59-91.
Bayoh MN, Thomas CJ, Lindsay SW: Mapping distributions of chromosomal forms of Anopheles gambiae in West Africa using climate data. Med Vet Entomol. 2001, 15: 267-274. 10.1046/j.0269-283x.2001.00298.x.
Lindsay SW, Parson L, Thomas CJ: Mapping the ranges and relative abundance of the two principal African malaria vectors, Anopheles gambiae sensu stricto and An. arabiensis, using climate data. Proc R Soc Lond B Biol Sci. 1998, 265: 847-854. 10.1098/rspb.1998.0369.
Costantini C, Ayala D, Guelbeogo WM, Pombi M, Some CY, Bassole IHN, Ose K, Fotsing J-M, Sagnon NF, Fontenille D, Besansky NJ, Simard F: Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae. BMC Ecol. 2009, 9: 16-10.1186/1472-6785-9-16.
Simard F, Ayala D, Kamdem GC, Etouna J, Ose K, Fotsing J-M, Fontenille D, Besansky NJ, Costantini C: Ecological niche partitioning between the M and S molecular forms of Anopheles gambiae in Cameroon: the ecological side of speciation. BMC Ecol. 2009, 9: 17-10.1186/1472-6785-9-17.
Lee Y, Meneses CR, Fofana A, Lanzaro GC: Desiccation resistance among subpopulations of Anopheles gambiae s.s. from Selinkenyi, Mali. J Med Entomol. 2009, 46: 316-320. 10.1603/033.046.0216.
Cassone BJ, Mouline K, Hahn MW, White BJ, Pombi M, Simard F, Costantini C, Besansky NJ: Differential gene expression in incipient species of Anopheles gambiae. Mol Ecol. 2008, 17: 2491-2504. 10.1111/j.1365-294X.2008.03774.x.
Papandreou NC, Iconomidou VA, Willis JH, Hamodrakas SJ: A possible structural model of members of the CPF family of cuticular proteins implicating binding to components other than chitin. J Insect Physiol. 2010, 56: 1420-1426. 10.1016/j.jinsphys.2010.04.002.
Caputo B, Dani FR, Horne GL, N'Fale S, Diabate A, Turillazzi S, Coluzzi M, Costantini C, Priestman AA, Petrarca V, della Torre A:Comparative analysis of epicuticular lipid profiles of sympatric and allopatric field populations ofAnopheles gambiae s.s. molecular forms andAn. arabiensisfrom Burkina Faso (West Africa). Insect Biochem Mol Biol. 2007, 37: 389-398. 10.1016/j.ibmb.2007.01.002.
Caputo B, Dani FR, Horne GL, Petrarca V, Turillazzi S, Coluzzi M, Priestman AA, Della Torre A: Identification and composition of cuticular hydrocarbons of the major Afrotropical malaria vector Anopheles gambiae s.s. (Diptera: Culicidae): analysis of sexual dimorphism and age-related changes. J Mass Spectrom. 2005, 40: 1595-1604. 10.1002/jms.961.
Benoit JB, Denlinger DL: Suppression of water loss during adult diapause in the northern house mosquito, Culex pipiens. J Exp Biol. 2007, 210: 217-226. 10.1242/jeb.02630.
Benoit JB, Lopez-Martinez G, Phillips ZP, Patrick KR, Denlinger DL: Heat shock proteins contribute to mosquito dehydration tolerance. J Insect Physiol. 2010, 56: 151-156. 10.1016/j.jinsphys.2009.09.012.
Benoit JB, Lopez-Martinez G, Teets NM, Phillips SA, Denlinger DL: Responses of the bed bug, Cimex lectularius, to temperature extremes and dehydration: levels of tolerance, rapid cold hardening and expression of heat shock proteins. Med Vet Entomol. 2009, 23: 418-425. 10.1111/j.1365-2915.2009.00832.x.
Lehmann T, Dao A, Yaro AS, Adamou A, Kassogue Y, Diallo M, Sekou T, Coscaron-Arias C: Aestivation of the African malaria mosquito, Anopheles gambiae in the Sahel. Amer J Trop Med Hyg. 2010, 83: 601-606. 10.4269/ajtmh.2010.09-0779.
Adamou A, Dao A, Timbine S, Kassogue Y, Yaro AS, Diallo M, Traore SF, Huestis DL, Lehmann T: The contribution of aestivating mosquitoes to the persistence of Anopheles gambiae in the Sahel. Malar J. 2011, 10: 151-10.1186/1475-2875-10-151.
Powell JR, Petrarca V, Della Torre A, Caccone A, Coluzzi M: Population structure, speciation, and introgression in the Anopheles gambiae complex. Parassitologia. 1999, 41: 101-113.
We thank O. Wood for advice, and C. Costantini and J. H. Willis for helpful discussion and constructive comments on an earlier version of this manuscript. This work was supported by National Institutes of Health grant R01AI076584 (to NJB) and grants from Eck Institute of Global Health, University of Notre Dame (to ZS).
The authors declare no competing interests.
KR, NB and ZS conceived the study. CC selected the NDKO subcolonies of known karyotype. KR performed mosquito husbandry and treatments, and subsequent measurements of total dry body mass, cuticle thickness by SEM, and CHC composition by gas chromatography. KR, FL, CL and ZS analyzed results. KR, NB and ZS wrote the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Reidenbach, K.R., Cheng, C., Liu, F. et al. Cuticular differences associated with aridity acclimation in African malaria vectors carrying alternative arrangements of inversion 2La. Parasites Vectors 7, 176 (2014). https://doi.org/10.1186/1756-3305-7-176
- An. gambiae
- An. coluzzii
- Cuticular hydrocarbons
- Chromosomal inversion
- Desiccation resistance
- M and S molecular forms