# Combining two-dimensional gel electrophoresis and metabolomic data in support of dry-season survival in the two main species of the malarial mosquito Anopheles gambiae

Kévin Hidalgo, Karine Mouline, W. Mamai, Natacha Foucreau, Rock K.

Dabiré, Alain Bouchereau, Frédéric Simard, D Renault

► To cite this version:

Kévin Hidalgo, Karine Mouline, W. Mamai, Natacha Foucreau, Rock K. Dabiré, et al. Combining two-dimensional gel electrophoresis and metabolomic data in support of dry-season survival in the two main species of the malarial mosquito Anopheles gambiae. Data in Brief, 2015, 5, pp.255-268. 10.1016/j.dib.2015.08.031. hal-01234296

HAL Id: hal-01234296

https://univ-rennes.hal.science/hal-01234296

Submitted on 29 Apr 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Data Article

# Combining two-dimensional gel electrophoresis and metabolomic data in support of dry-season survival in the two main species of the malarial mosquito *Anopheles gambiae* ☆

<figure><img src="image_2.png" /></figure>

CrossMark

Hidalgo K.<sup>a,b,*</sup>, Mouline K.<sup>b,c</sup>, Mamai W.<sup>b,c</sup>, Foucreau N.<sup>d</sup>, Dabiré K.R.<sup>c</sup>, Bouchereau A.<sup>e</sup>, Simard F.<sup>b</sup>, Renault D.<sup>a</sup>

ARTICLE INFO

Article history:

Received 24 June 2015

Received in revised form

25 August 2015

Accepted 25 August 2015

Available online 7 September 2015

Keywords:

2D electrophoresis

Amino acid

Metabolomic

Mosquito

Proteomic

Sugar

ABSTRACT

In dry savannahs of West-Africa, the malarial mosquitoes of the *Anopheles gambiae sensu stricto* complex annually survive the harsh desiccating conditions of the dry season. However, the physiological and biochemical mechanisms underlying how these mosquitoes survive such desiccating conditions are still undefined, and controversial. In this context, we provide the first work examining both proteomic and metabolomic changes in the two molecular forms of *A. gambiae* s.s (M and S forms) experimentally exposed to the rainy and dry season conditions as they experience in the field. Protein abundances of the mosquitoes were measured using a two-dimensional fluorescence difference gel electrophoresis (2D DIGE) coupled with a matrix-assisted laser desorption/ionisation-time of flight (MALDI-TOF) and tandem mass spectrometry (MS) for protein identification. These assays were

<footnote><sup>a</sup> Université de Rennes 1, UMR CNRS 6553 Ecobio, Campus de Beaulieu, 263 Avenue du Général Leclerc, CS 74205 35042 Rennes Cedex, France</footnote>

<footnote><sup>b</sup> Institut de Recherche pour le Développement (IRD), UMR IRD 224-CNRS 5290-Université de Montpellier MIVEGEC, 911 Avenue Agropolis, BP 64501, 34394 Montpellier cedex 5, France</footnote>

<footnote><sup>c</sup> Institut de Recherche en Sciences de la Santé (IRSS), Direction Régionale de l'Ouest (DRO), 399 Avenue de la Liberté 01, BP 545 Bobo-Dioulasso, Burkina Faso</footnote>

<footnote><sup>d</sup> Université Claude Bernard Lyon 1, UMR CNRS 5023 LEHNA, 43 Bd du 11 Novembre 1918, 69622 Villeurbanne Cedex, France</footnote>

<footnote><sup>e</sup> Université de Rennes 1, UMR INRA IGEPP, Campus de Beaulieu, 263 Avenue du Général Leclerc, CS 74205, 35042 Rennes Cedex, France</footnote>

<footnote><sup>∗</sup>For more details, please refer to: K. Hidalgo, K. Mouline, W. Mamai, N. Foucreau, K.R. Dabiré, A. Bouchereau, F. Simard and D. Renault, Novel insights into the metabolic and biochemical underpinnings assisting dry-season survival in female malaria mosquitoes of the *Anopheles gambiae* complex, Journal of Insect Physiology (2014). doi:10.1016/j.jinsphys.2014.07.003.</footnote>

<footnote>* Corresponding author.</footnote>

<footnote>E-mail address: kevin.hidalgo@gmx.com (K. Hidalgo).</footnote>
conducted by Applied Biomics (http://www.appliedbiomics.com, Applied Biomics, Inc. Hayward, CA, USA), and the mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifier PXD000294. The metabolomic analysis was conducted using both Acquity UPLC® system (for amino acid identification), and a gas-chromatography-mass spectrometry platform (for sugars identification). Metabolomic fingerprints were assessed in the University of Rennes 1, UMR CNRS 6553 EcoBio (France). A detailed interpretation of the obtained data can be found in Hidalgo et al. (2014) [1] (Journal of Insect Physiology (2014)).

© 2015 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license

(http://creativecommons.org/licenses/by/4.0/).

<table><caption>Specifications Table [please fill in right-hand column of the table below]</caption><tr><td>Subject area</td><td>Biology</td></tr><tr><td>More specific subject area</td><td>Ecophysiology; Proteomics; Medical Entomology</td></tr><tr><td>Type of data</td><td>Tables &amp; barplots</td></tr><tr><td>How data was acquired</td><td>Proteomic data were assessed using a two-dimensional fluorescence difference gel electrophoresis (2D DIGE) coupled with a matrix-assisted laser desorption/ionisation-time of flight (MALDI-TOF) and tandem mass spectrometry (MS) Metabolomic data were assessed using a gas-chromatography-mass spectrometry (GC-MS) and an Acquity UPLC® system (Waters Corporation, Milford, MA, USA) platforms</td></tr><tr><td>Data format</td><td>Raw, analyzed</td></tr><tr><td>Experimental factors</td><td>We compared the protein and metabolites (amino acids, polyols, sugars) in two malarial mosquito species form the Anopheles gambiae species complex (formerly known as the M and S molecular forms). The two species were reared under climatic conditions mimicking the rainy and dry season conditions as they experience naturally.</td></tr><tr><td>Experimental features</td><td>Two-dimensional fluorescence difference gel electrophoresis (2D DIGE); matrix-assisted laser desorption/ionisation-time of flight (MALDI-TOF) and tandem mass spectrometry (MS); gas-chromatography-mass spectrometry (GC-MS) and Acquity UPLC® platforms</td></tr><tr><td>Data source location</td><td>Both M and S mosquitoes were collected in Burkina-Faso (West Africa) in human dwellings of Bama (11°23'N, 04°24'W) and Soumoussou (11°01'N, 04°02'W), respectively. The metabolomic data were collected at the University of Rennes 1, UMR CNRS 6553 EcoBio (France). The proteomic data were obtained from Applied Biomics (USA).</td></tr><tr><td>Data accessibility</td><td>The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifier PXD000294.</td></tr></table>

Value of the data

* The use of 2D gel electrophoresis comparisons identified interesting changes in several protein abundances from one season to the next in mosquitoes;

* The comparisons of mosquito's metabolic fingerprints highlighted specific changes in their metabolite contents (amino acids and sugars) from one season to the next;
* By combining both 2D gel electrophoresis and metabolic fingerprint comparisons our approach improves the identification of physiological and biochemical changes in mosquitoes.

# Data, experimental design, materials and methods

## 1. Data

Mosquitoes from the *Anopheles gambiae* s.l. complex (both M and S molecular forms of *A. gambiae*) are widespread in West Africa where they seasonally experienced the harsh desiccating conditions of the dry season. How these mosquitoes survive the dry season and are engaged in exponential population growth as soon as the rainy season starts is still undefined and very controversial [2,3]. This ability to survive the dry season suggests high physiological plasticity in mosquitoes. In order to highlight the general physiological changes involved in mosquitoes at the onset of the dry season, we compared protein abundance and metabolic fingerprints in the two molecular forms of *A. gambiae* exposed to the climatic conditions of the rainy and the early dry season.

The M and S molecular forms of *A. gambiae* were recently described as two different species, named *Anopheles coluzzii* and *A. gambiae*, respectively [4,5]. Here, we are still referring to the M and S molecular forms, as in our published paper (Hidalgo et al. [1]).

## 2. Experimental design

Mosquitoes were reared from eggs to adults into four programmable climatic chambers (Sanyo MLR 315 H, Sanyo Electric Co., Osaka, Japan). Two climatic chambers were used to reproduce the climatic conditions of the rainy season, and two others to reproduce those of the dry season. Climatic conditions were programmed using the temperature and humidity cycles recorded in south-western Burkina-Faso (11°23′N, 04°24′W) with a Vantage Pro2 weather monitoring station (Weatherlink; Davis Instruments, Hayward, CA, U.S.A.). Climatic conditions were recorded in August 2010 (rainy season) and December 2010 (onset of the dry season), then hourly averaged to design 12 step cycles required to programme as closely as possible the natural daily climatic variations inside the climatic chambers [6]. Conditions inside the climatic chambers were tightly monitored throughout experiments using MicroLog Pro monitors (EC750, Davis Instruments, Hayward, CA, U.S.A.).

## 3. Proteomic assays

For each mosquito species (M and S molecular forms) and rearing conditions (rainy and dry season), proteomic samples consisted of a pool of 30 teneral (1 h-old adult) female mosquitoes. Samples were sent to the proteomics department of Applied Biomics (http://www.appliedbiomics.com, Applied Biomics, Inc. Hayward, CA, USA) for proteomic assays as described below.

### 3.1. Protein extraction and adjustment of the concentration level

A 2-D cell-lysis buffer (30 mM Tris-HCl, pH 8.8, containing 7 M urea, 2 M thiourea, and 4% 3-[(3-Cholamidopropyl) dimethylammonio]-1 propanesulfonate [CHAPS]) was used to extract the proteins in each sample of 30 pooled females. Then, the protein concentration of each sample was measured using a Bio-Rad protein assay method, and adjusted to 5 mg.ml⁻¹ using the 2-D cell-lysis buffer.

### 3.2. Samples labelling

For both species, we used the CyDye DIGE Fluor Cy3 (red) and Cy5 (green) saturation dyes (4 mL, 2 mM) to label samples from the rainy and the dry season rearing conditions, respectively. A Cy2 dye,
corresponding to an internal standard for the normalisation of protein abundances, was also used to label a 50:50 mixture of the two rearing conditions samples. The three CyDye DIGE Fluors were incubated with corresponding samples for 30 min on ice, under dark conditions, before adding a 1 µL volume of 10 mM lysine. A 2 × 2-D sample buffer (8 M urea, 4% CHAPS, 20 mg/mL DTT, 2% phar- malytes and a trace amount of bromophenol blue) was used to terminate the reaction, and the samples were stored at –80 °C until further use for 2D DIGE assays.

### 3.3. Two-dimensional fluorescence difference gel electrophoresis (2D DIGE)

First, 30 µg of both Cy3- and Cy5-labelled samples were mixed with a DeStreak™ rehydration solution (0.5% pharmalyte pH 3–10 in a 250 µL final volume; GE Healthcare, Waukesha, WI, USA). Then, samples were loaded onto an isoelectric focusing (IEF) strip (pH 4–10 linear range; GE Healthcare) with an Ettan IPGPhore II at standard conditions. After IEF, the strips were incubated in an equilibration buffer (50 mM Tris–HCl, pH 8.8, containing 6 M urea, 30% glycerol, 2% SDS, a trace amount of bromophenol blue and 10 mg.mL⁻¹ DTT) for 15 min, and rinsed in an equilibration buffer 2 (50 mM Tris–HCl, pH 8.8, containing 6 M urea, 30% glycerol, 2% SDS, a trace amount of bromophenol blue and 45 mg.mL⁻¹ iodoacetamide) for 10 min.

The electrophoresis was performed at 15 °C on a 12% SDS-PAGE gel. Note that each gel contained the three Cy3, Cy5 and Cy2 labelled samples. For each mosquito species, we performed three-replicate gels to compare the effects of rearing conditions on their protein expression (N=6 gels in total). A seventh gel containing the two species exposed to dry conditions and their 50:50 mixture was run. In this last gel, the M molecular form was labelled with Cy5 dye, whereas the S form was labelled with Cy3 dye.

### 3.4. Image analysis

Effect of rearing conditions on protein expression in both species and differences between species were assessed by overlaying images of the seven different gels. Note that the abundance of each protein spot was normalised using the seventh gel. The resulting 2D gels were scanned using a Typhoon Trio scanner (Amersham BioSciences) with excitation and emission wavelengths for Cy2-labelled (488/520 nm), Cy3-labelled (548/560 nm) and Cy5-labelled (641/660 nm) proteins using settings that match in similar relative fluorescence intensities for the Cy3- or Cy5-labelled samples (Fig. 1). Image analysis for intensity measurements of the different protein spots was performed using ImageQuantTL and DeCyder softwares (GE Healthcare). Then, images were subjected to in-gel analysis and cross-gel analysis using DeCyder v6.5 (GE-Healthcare) with a *P*-value < 0.05 according to Student's *t*-tests. Based on these tests, we retained a total of 109 spots, corresponding to proteins with molecular masses ranging from 14 to 150 kDa, and isoelectric points between 4 and 9. These spots showed significant variation in at least one species and/or one rearing condition.

### 3.5. Protein identification using matrix-assisted laser desorption/ionisation-time of flight (MALDI-TOF) and tandem mass spectrometry (MS)

Identification of the 109 protein spots of interest was performed by Applied Biomics, Inc. Spots were picked up from the gel using an Ettan spot picker (GE Healthcare), and washed twice with 25 mM ammonium bicarbonate and 50% acetonitrile to remove the remaining dyes. Then, proteins were digested in-gel at 37 °C using a modified porcine trypsin protease (Trypsin Gold, Promega) according the method described in Rosenfeld et al. [7]. The digested tryptic peptides were desalted using Zip-tip C18 (Millipore, Billerica, MA), mixed with 0.5 µL of α-cyano-4-hydroxycinnamic acid (CHCA) matrix and spotted into wells of a MALDI plate. Mass spectra of the peptides in each digested spot were obtained using a MALDI-TOF (MS) and a TOF/TOF (tandem MS/MS) equipment (AB Sciex). The MALDI-TOF mass spectra data were acquired in reflectron positive ion mode, averaging 2000 laser shots per spectrum. The TOF/TOF tandem MS fragmentation spectra were acquired for each sample, averaging 2000 laser shots per fragmentation spectrum on each of the 5-10 most abundant ions present in each sample (excluding trypsin autolytic peptides and other known background ions). Both
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Representative image of the two-dimensional fluorescence difference gel electrophoresis (2D DIGE); *A. gambiae* M (A), *A. gambiae* S (B). Females reared under the dry conditions were labelled with Cy3 (green), whereas females reared under the rainy conditions were labelled with Cy5 (red). From Hidalgo et al. [1].</figcaption></figure>
the peptide mass and the associated fragmentation spectra were submitted for database search using GPS Explorer software version 3.5 (Applied Biosystems) equipped with the MASCOT search engine ([http://www.matrixscience.com](http://www.matrixscience.com), Matrix science). This submission allows the identification of proteins from the National Center for Biotechnology Information non-redundant protein sequence database (NCBI nr), restricted to the *A. gambiae* species complex. The searches were performed without con- straining the protein molecular weight or isoelectric point, with carbamidomethylation and oxidation as variable modifications, and with one missed cleavage allowed in the search parameters. The highest protein scoring hit with a protein score confidence interval over 95% from the database search for each 2D gel spot was accepted as positive identification.

Based on an average difference of at least 1.3-fold (absolute value), only 39 spots were selected, which corresponded to 29 distinct proteins, for further identification. The ratio of the 39 spot expressions between the species and/or rearing conditions was obtained from in-gel analysis using DeCyder. Table 1 shows the identification of the 39 spots retained in *A. gambiae* based on this method.

### 3.6. Data analysis

Average-fold differences of the 39 spot abundances were first normalised (Logarithmic abundance of the protein contents from each 2D DIGE gel), then calculated for the four following comparisons: (i) rainy vs. dry condition in *A. gambiae* M; (ii) rainy vs. dry condition in *A. gambiae* S; (iii) *A. gambiae* M vs. *A. gambiae* S when exposed to rainy conditions; and (iv) *A. gambiae* M vs. *A. gambiae* S when exposed to dry conditions. For each ratio, differences in spot abundance were assessed using Student's *t*-tests on individual spots (*P*-value < 0.05).

Accordingly, we observed that 15 proteins exhibited significant variations between the two rearing conditions in 1-h old females *A. gambiae* M, among which 10 were more abundant when females grew under the dry conditions (Fig. 2A). In contrast, 14 proteins showed significant variations between the two rearing conditions in 1-h old females *A. gambiae* S, among which 8 were more abundant when females grew under dry conditions (Fig. 2B). Overall, nine proteins exhibited the same pattern of variations in both species: six proteins increased under dry conditions, whereas the three others decreased under these conditions (Fig. 2A and B).

Finally, 24 proteins showed significant variation between species whatever the rearing conditions. Six of the 24 proteins were always more abundant in females *A. gambiae* S (Fig. 3). For details about data interpretation and conclusion please report to Hidalgo et al. [1].

## 4. Metabolomic assays: amino acid and sugar concentrations

Metabolomic fingerprinting was conducted on both 1 h-old and 24 h-old adult females. For each mosquito species (M and S molecular forms) and rearing conditions (rainy and dry season), meta- bolomic samples consisted of a pool of five-six female specimens, so that the dry mass of each sample reached at least 1 mg. Five to nine replicates were used for the different experimental conditions.

### 4.1. Extraction of the metabolites

Samples were first freeze-dried (Lyovac™ GT3) for 72 h, before being weighed using a micro- balance (Mettler Toledo GmbH®, Greinfense, Switzerland, *d*=1 µg). Two 3 mm tungsten beads and a 1000 µL volume of methanol-chloroform (2:1) solution were added to each sample. Then, each sample was shaken at 30 Hz for 1.5 min (Retsch™ MM301, bead-beating, Retsch GbM, Haan, Germany). We added a 500 µL volume of ultrapure water. Samples were vortexed for homogenisation and centrifuged at 4000g for 10 min at 4 °C. Finally, we collected and transferred a 600 µL aliquot of the upper aqueous phase (containing amino acids, polyols and sugars) into clean microtubes. The 600 µL aliquot was vacuum-dried (Speed Vac Concentrator, Genevac Ltd., Ipswitch, England), and a volume of 600 µL of ultrapure water was added to the residual.
<caption>**Table 1**</caption>

List of the 39 protein spots displaying differential abundance (average fold difference > 1.3, absolute value; *P* < 0.05) between the two rearing conditions in 1-h old female *A. gambiae* M and S, or between the two species whatever the rearing conditions. Adapted from Hidalgo et al. [1].

<table><thead><tr><th>Spot no.<sup>a</sup></th><th>Protein name</th><th>Protein scores<sup>b</sup></th><th>No. Of matched peptides</th><th>Source species</th><th>Accession no.<sup>c</sup></th><th>Functions<sup>d</sup></th></tr></thead><tbody><tr><td>01</td><td>tropomyosin invertebrate</td><td>318</td><td>12</td><td>Culex quinquefasciatus</td><td>gil170056897</td><td>Muscle contraction</td></tr><tr><td>02</td><td>hexamerin 2 beta</td><td>41</td><td>2</td><td>Aedes aegypti</td><td>gil157119837</td><td>Storage of amino acid</td></tr><tr><td>03</td><td>arylphorin subunit alpha</td><td>134</td><td>4</td><td>Culex quinquefasciatus</td><td>gil170043201</td><td>Storage of amino acid / Constituent of Cuticle sclerotizing system</td></tr><tr><td>07</td><td>hexamerin A</td><td>284</td><td>15</td><td>Anopheles melas</td><td>gil3420171</td><td>Storage of amino acid</td></tr><tr><td>09</td><td>hypothetical protein AND_22551</td><td>153</td><td>15</td><td>Anopheles darlingi</td><td>gil312371166</td><td>Storage of amino acid</td></tr><tr><td>12</td><td>hexamerin 2 beta</td><td>70</td><td>4</td><td>Aedes aegypti</td><td>gil157110143</td><td>Storage of amino acid</td></tr><tr><td>13</td><td>myosin-1d</td><td>172</td><td>5</td><td>Culex quinquefasciatus</td><td>gil170029188</td><td>Muscle contraction</td></tr><tr><td>19</td><td>catalase</td><td>215</td><td>16</td><td>Anopheles gambiae</td><td>gil118638436</td><td>Response to oxidative stresses</td></tr><tr><td>21</td><td>hypothetical protein AND_22551</td><td>321</td><td>12</td><td>Anopheles darlingi</td><td>gil312371166</td><td>Storage of amino acid</td></tr><tr><td>24</td><td>AGAP005558-PA</td><td>1000</td><td>19</td><td>Anopheles gambiae</td><td>gil31213235</td><td>Proteolysis</td></tr><tr><td>25</td><td>AGAP005558-PA</td><td>570</td><td>16</td><td>Anopheles gambiae</td><td>gil31213235</td><td>Proteolysis</td></tr><tr><td>27</td><td>aspartate ammonia lyase</td><td>301</td><td>15</td><td>Aedes aegypti</td><td>gil157118058</td><td>Tricarboxylic acid cycle</td></tr><tr><td>37</td><td>myosin regulatory light chain 2 (mlc-2)</td><td>336</td><td>5</td><td>Aedes aegypti</td><td>gil157167683</td><td>Muscle contraction</td></tr><tr><td>38</td><td>myosin regulatory light chain 2 (mlc-2)</td><td>299</td><td>5</td><td>Aedes aegypti</td><td>gil157167683</td><td>Muscle contraction</td></tr><tr><td>39</td><td>myosin regulatory light chain 2 (mlc-2)</td><td>330</td><td>5</td><td>Aedes aegypti</td><td>gil157167683</td><td>Muscle contraction</td></tr><tr><td>45</td><td>AGAP007963-PA</td><td>601</td><td>16</td><td>Anopheles gambiae</td><td>gil118789564</td><td>Muscle contraction</td></tr><tr><td>52</td><td>ATP synthase beta subunit</td><td>250</td><td>12</td><td>Culex quinquefasciatus</td><td>gil170040305</td><td>ATP synthesis</td></tr><tr><td>55</td><td>AGAP011284-PA</td><td>304</td><td>11</td><td>Anopheles gambiae</td><td>gil118779554</td><td>Positive regulation of translational elongation and termination</td></tr><tr><td>56</td><td>NADH dehydrogenase</td><td>268</td><td>6</td><td>Culex quinquefasciatus</td><td>gil170037145</td><td>Respiratory chain process</td></tr><tr><td>57</td><td>myosin light chain 1, putative</td><td>202</td><td>4</td><td>Aedes aegypti</td><td>gil157167807</td><td>Muscle contraction</td></tr><tr><td>58</td><td>myosin light chain 1, putative</td><td>117</td><td>4</td><td>Aedes aegypti</td><td>gil157167807</td><td>Muscle contraction</td></tr><tr><td>59</td><td>myosin light chain 1, putative</td><td>156</td><td>5</td><td>Aedes aegypti</td><td>gil157167807</td><td>Muscle contraction</td></tr><tr><td>60</td><td>myosin light chain 1, putative</td><td>181</td><td>5</td><td>Aedes aegypti</td><td>gil157167807</td><td>Muscle contraction</td></tr><tr><td>64</td><td>AGAP010657-PA</td><td>165</td><td>8</td><td>Anopheles gambiae str. PEST</td><td>gil158289706</td><td>Oxygen transporter activity</td></tr><tr><td>65</td><td>*cuticular protein 117, RR-2 family (AGAP003379-PA)</td><td>709</td><td>9</td><td>Anopheles gambiae</td><td>gil158290052</td><td>Structural constituent of the rigid cuticle</td></tr><tr><td>67</td><td>*cuticular protein 70, RR-2 family (AGAP006283-PB)</td><td>608</td><td>9</td><td>Anopheles gambiae</td><td>gil158295676</td><td>Structural constituent of the rigid cuticle</td></tr><tr><td>68</td><td>hexamerin A</td><td>395</td><td>26</td><td>Anopheles gambiae</td><td>gil3420159</td><td>Storage of amino acid</td></tr><tr><td>33</td><td>myosin heavy chain, nonmuscle or smooth muscle</td><td>124</td><td>17</td><td>Aedes aegypti</td><td>gil157111095</td><td>Muscle contraction</td></tr></tbody></table>
<table><caption>**Table 1 (continued)**</caption><thead><tr><th>Spot no.<sup>a</sup></th><th>Protein name</th><th>Protein scores<sup>b</sup></th><th>No. Of matched peptides</th><th>Source species</th><th>Accession no.<sup>c</sup></th><th>Functions<sup>d</sup></th></tr></thead><tbody><tr><td>73</td><td>myosin heavy chain, nonmuscle or smooth muscle</td><td>498</td><td>36</td><td>Aedes aegypti</td><td>gil157110721</td><td>Muscle contraction</td></tr><tr><td>75</td><td>glycogen phosphorylase</td><td>635</td><td>21</td><td>Aedes aegypti</td><td>gil157108521</td><td>Glycogenolysis process</td></tr><tr><td>76</td><td>prophenoloxidase</td><td>955</td><td>38</td><td>Anopheles gambiae</td><td>gil3892092</td><td>Oxidation-reduction process</td></tr><tr><td>80</td><td>AGAP007393-PB</td><td>570</td><td>19</td><td>Anopheles gambiae str. PEST</td><td>gil118778070</td><td>Cell redox homoeostasis / Glycerol ether metabolic process</td></tr><tr><td>84</td><td>citrate synthase</td><td>58</td><td>12</td><td>Aedes aegypti</td><td>gil157133341</td><td>Tricarboxylic acid cycle</td></tr><tr><td>85</td><td>isocitrate dehydrogenase cytoplasmic</td><td>283</td><td>14</td><td>Culex quinquefasciatus</td><td>gil170028051</td><td>Tricarboxylic acid cycle</td></tr><tr><td>92</td><td>AGAP009218-PA</td><td>141</td><td>4</td><td>Anopheles gambiae str. PEST</td><td>gil118791868</td><td>Proteolysis</td></tr><tr><td>93</td><td>AGAP005423-PA</td><td>588</td><td>11</td><td>Anopheles gambiae str. PEST</td><td>gil118786445</td><td>Ubiquitin-dependent protein catabolic process</td></tr><tr><td>96</td><td>arginine kinase</td><td>503</td><td>14</td><td>Mimomyia luzonensis</td><td>gil284159531</td><td>Transferase / Kinase</td></tr><tr><td>99</td><td>glutathione transferase epsilon3</td><td>358</td><td>4</td><td>Anopheles dirus</td><td>gil74275399</td><td>Transferase</td></tr><tr><td>101</td><td>cuticular protein 123, RR-2 family (AGAP003385-PA)</td><td>814</td><td>9</td><td>Anopheles gambiae str. PEST</td><td>gil158290062</td><td>Structural constituent of the rigid cuticle</td></tr></tbody></table>

<footnote><sup>a</sup> Protein scores derived from Mascot algorithm, indicating identity or extensive homology (P < 0.05).</footnote>

<footnote><sup>b</sup> Protein accession numbers from the National Center for Biotechnology Information non-redundant (NCBI) database.</footnote>

<footnote><sup>c</sup> Positive values mean more expressed in "rs" (rainy season) conditions whereas negative values mean over expressed in "ods" (onset of the dry season) conditions.</footnote>

<footnote><sup>d</sup> Protein functions are checked using http://www.uniprot.org.</footnote>
A

An. gambiae M

<figure><img src="image_5.png" /></figure>

B

An. gambiae S

<figure><img src="image_6.png" /></figure>

**Fig. 2.** Differential fold-change of the 21 selected spot proteins between rainy and dry conditions in 1-h old females of *Anopheles gambiae* M (A) and S (B). Blue bars are positive values that represent spots that were more abundant in rainy conditions, whereas red bars are negative values that represent spots that were more abundant in dry conditions. Spots with "*" represent the nine protein spots exhibiting the same pattern of variations in both species. From Hidalgo et al. [1].

## 4.2. Identification and quantification of amino acids

Twenty µL of each sample previously extracted were diluted with 60 µL of ultrapure water. Then, a 5 µL volume of each diluted sample was transferred into a new microtube for derivatization. Derivatization was conducted according to the AccQTag ultra derivatization kit (Waters Corporation, Milford, MA, USA) procedures. After reconstituting the AccQTag Ultra Reagent Powder, we prepared a
<figure><img src="image_7.png" /><figcaption>**Fig. 3.** Differential fold-change of the 24 selected spot proteins between M and S forms reared under rainy (blue bars) and dry (red bars) conditions. Positive values represent spots that were more abundant in *Anopheles gambiae* M, whereas negative values represent spots that were more abundant in *Anopheles gambiae* S. Six proteins were more abundant in females *A. gambiae* S, whatever the experimental conditions, among which prophenoloxidase (spot #76), citrate synthase (spot #84), isocitrate dehydrogenase (spot #85), and glutathione transferase (spot #99) can be mentioned. The RR-2 cuticular protein 123 (spot #101) was also significantly more abundant under the dry conditions (Student *t*-test, *P* < 0.001), although its fold change was only 1.27. In addition, a phosphorylase (spot #75) was more abundant in *A. gambiae* S compared to *A. gambiae* M when the females were exposed to the dry conditions, but no differential expression was found under rainy conditions. Arginine kinase (spot #96) was always more abundant in *A. gambiae* M, whatever the experimental conditions. From Hidalgo et al. [1].</figcaption></figure>

100 pmol/µl calibration standard solution (Acide Hydrolase Standard). Then, a 10 µl volume of the calibration standard solution and a 20 µl volume of AccQTag Ultra Reagent were added to each sample. Samples were vortexed and heated for 10 min at 55 °C.

After derivatization, 1 µL of each sample was transferred in an Acquity UPLC® system (Waters Corporation, Milford, MA, USA) equipped with an Acquity UPLC® BEH C18 1.7 µm 2.1 × 100 mm² column heated at 55 °C, as described by Renaultet al. [8]. The derivatized amino acids were detected at 260 nm using a photo diode-array detector. Peaks were identified according to their retention time compared with a list of 22 amino acid commercial standards (Table 2). Concentration was calculated by comparison of each amino acid peak area with those of the individual external standards.

### 4.3. Identification and quantification of sugars

A 50 µL volume of each sample previously extracted was transferred into a clean glass vial and vacuum-dried. Then, residuals were resuspended in 50 µL of freshly prepared methoxyamine hydrochloride (Sigma-Aldrich, St. Louis, MO, USA) in pyridine (20 mg/mL). Samples were incubated under orbital shaking at 30 °C for 90 min. Following incubation, 50 µl volume of N-methyl-N-(tri-methylsilyl) trifluoroacetamide (MSTFA) were added for derivatization. Derivatization was conducted overnight at 37 °C.

To identify and quantify sugars, we used gas chromatography coupled with mass spectrometry (GC-MS). The GC-MS system consisted of a TriPlus autosampler, a Trace GC Ultra chromatograph and a Trace DSQII quadrupole mass spectrometer (Thermo Fischer Scientific Inc, Waltham, MA, USA). A 1 µL volume of the derivatized sample was injected into the GC-MS system using the split mode (25:1). The injector temperature was held at 260 °C. We programmed the GC-MS in order that the oven temperature remained 4 min at 70 °C after the injection, and then increased at a rate of 5 °C/min until it reached 300 °C, where it remained for 10 min. For elution, we used a 30 m fused silica column (TR5 MS, I.D. 0.25 mm, 95% dimethyl siloxane, 5% phenyl polysilphenylene-siloxane) and helium as the carrier gas at a rate of 1 mL/min.
<table><caption>**Table 2**
List of the 22 standard amino acids used for metabolomic assays.</caption><thead><tr><td>Name</td><td>abbreviation</td></tr></thead><tbody><tr><td>Alanine</td><td>Ala</td></tr><tr><td>Arginine</td><td>Arg</td></tr><tr><td>Asparagine</td><td>Asn</td></tr><tr><td>Aspartate</td><td>Asp</td></tr><tr><td>Cystein</td><td>Cys</td></tr><tr><td>GABA</td><td>GABA</td></tr><tr><td>Glutamate</td><td>Glu</td></tr><tr><td>Glutamine</td><td>Gln</td></tr><tr><td>Glycine</td><td>Gly</td></tr><tr><td>Histidine</td><td>His</td></tr><tr><td>Isoleucine</td><td>Ile</td></tr><tr><td>Leucine</td><td>Leu</td></tr><tr><td>Lysine</td><td>Lys</td></tr><tr><td>Methionine</td><td>Met</td></tr><tr><td>Ornithine</td><td>Orn</td></tr><tr><td>Phenylalanine</td><td>Phe</td></tr><tr><td>Proline</td><td>Pro</td></tr><tr><td>Serine</td><td>Ser</td></tr><tr><td>Threonine</td><td>Thr</td></tr><tr><td>Tryptophan</td><td>Trp</td></tr><tr><td>Tyrosine</td><td>Tyr</td></tr><tr><td>Valine</td><td>Val</td></tr></tbody></table>

Detection of each sugar was conducted by mass spectrometry using the electronic impact (EI) method. Briefly, the temperature of the ion source and mass spectrometry transfer line was set at 260 °C. All of the samples were run under the SIM mode (electron energy: –70 eV). Each sugar peak was annotated using both mass spectra (two specific ions), and the retention index specific to each compound. Randomised sample and standard sequences were established for the sample injection. Chromatograms were deconvoluted using the XCalibur v2.0.7 software (Thermo Fisher Scientific, San Jose, CA, USA). Standard samples consisted of the pure reference compounds (arabinose, fructose, galactose, glucose, ribose and trehalose) at 100, 200, 300, 500 and 1000 µM concentrations. Sugars concentrations were thus quantified using the quadratic calibration curves for each pure reference compound.

## 4.4. Data analysis

### 4.4.1. Amino acids

Among the 22 amino acids we listed in our external standards, 21 were correctly detected and quantified (i.e. cysteine did not provide reliable measurements, as this amino acid is unstable and is quickly transformed into cystine during derivatization). Significance of amino acid fingerprintings among species and rearing conditions was assessed using multivariate discriminant analyses and MANOVAs. Because false discovery rate can be induced by multiple comparisons, we used a Benjamini–Hochberg procedure to adjust all *P*-values.

Comparisons of the amino acid fingerprints of each experimental modality showed significant differences in the concentration of at least one amino acid between the two species, their ageing and their rearing conditions. Of particular interest, both species and especially M form, accumulated phenylalanine, tryptophan, tyrosine, and valine 1 h after emergence (Fig. 4A and B). S form specimens did not show fingerprint variation according to the rearing conditions 1 h after emergence, but significant changes were observed 24 h after emergence, with proline being accumulated under dry conditions (Fig. 4B). For details about data interpretation and conclusion, please report to Hidalgo et al. [1].
A

An. gambiae M

<figure><img src="image_8.png" /></figure>

B

An. gambiae S

<figure><img src="image_9.png" /><figcaption>**Fig. 4.** Left: Sample projections onto the first LDA discriminant plane for female *Anopheles gambiae* M (28 samples, A) and S (23 samples, B) species reared under dry (red squares) and rainy (blue dots) season conditions. Singular values refer to the ratio of the between-class and within-class inertia. Right: correlation circles depict the normalised relation (from -1 to 1) between each amino acid and LDA axes. From Hidalgo et al. [1].</figcaption></figure>

### 4.4.2. Sugars

Among the six sugars listed in our reference standard, only galactose and glucose were reliably quantified in our mosquito samples. Indeed, arabinose, fructose, ribose and trehalose were below the quantification limit of the equipment used. Significance of the galactose and glucose concentrations variation among species and rearing conditions was assessed using Kruskal–Wallis tests, followed by a Bonferroni procedure resulting in a decreased threshold of significance from P-value < 0.05 to P-value < 0.01.
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 5.** Mean (± s.e.) galactose (A and B) and glucose (C and D) levels (in n moles mg⁻¹ dry mass) in 1-h and 24-h old females of *Anopheles gambiae* M (A and C) and S (B and D). Blue bars correspond to the females reared under rainy season conditions and light red bars correspond to females reared under dry season conditions. Letters above the bars report significant differences after Bonferroni correction to account for multiple tests (*P* < 0.01). From Hidalgo et al. [1].</figcaption></figure>

Accordingly, both galactose and glucose concentrations were significantly higher 24 hours after emergence in the two species when mosquitoes were reared under the dry conditions (Fig. 5A–D). Moreover, glucose concentration was lower in M form specimens reared under the dry conditions (Fig. 5B and C). For details about data interpretation and conclusion please report to Hidalgo et al. [1].

**Acknowledgements**

These assays were supported by the Agence Nationale de la Recherche through grant ANR-08-MIEN-006 to FS. The authors would like to thank Moundai Tchonfienet, Sougrinoma Zoungrana and Boubakar Nikiema (IRSS, Bobo-Dioulasso) for their help in mosquito rearing and sample preparation. Authors would like to thank the PRIDE team for their support.
# Appendix A. Supplementary material

Supplementary data associated with this article can be found in the online version at http://dx.doi.org/10.1016/j.dib.2015.08.031.

## References

[1] K. Hidalgo, K. Mouline, W. Mamai, N. Foucreau, K.R. Dabiré, A. Bouchereau, et al., Novel insights into the metabolic and biochemical underpinnings assisting dry-season survival in female malaria mosquitoes of the *Anopheles gambiae* complex, *J. Insect Physiol.* (2014), http://dx.doi.org/10.1016/j.jinsphys.2014.07.003.

[2] T. Lehmann, A. Dao, A.S. Yaro, M. Diallo, S. Timbiné, D.L. Huestis, et al., Seasonal Variation in Spatial Distributions of *Anopheles gambiae* in a Sahelian Village: Evidence for Aestivation, *J. Med. Entomol* 51 (2014) 27–38. http://dx.doi.org/10.1603/ME13094.

[3] D.L. Huestis, T. Lehmann, Ecophysiology of *Anopheles gambiae* s.l.: Persistence in the Sahel, *Infect. Genet. Evol.* J. Mol. Epidemiol. Evol. Genet. Infect. Dis (2014), http://dx.doi.org/10.1016/j.meegid.2014.05.027.

[4] M. Coetzee, R.H. Hunt, R. Wilkerson, A.D. Torre, M.B. Coulibaly, N.J. Besansky, *Anopheles coluzzii* and *Anopheles amharicus*, new members of the *Anopheles gambiae* complex, *Zootaxa* (2013), http://dx.doi.org/10.11646/zootaxa.3619.3.2.

[5] Y. Lee, C.D. Marsden, L.C. Norris, T.C. Collier, B.J. Main, A. Fofana, et al., Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, *Anopheles gambiae*, *Proc. Natl. Acad. Sci* 110 (2013) 19854–19859.

[6] K. Mouline, W. Mamai, P. Agnew, M. Tchonfienet, C. Brengues, R. Dabire, et al., Physiology and development of the M and S molecular forms of *Anopheles gambiae* in Burkina Faso (West Africa), *Med. Vet. Entomol.* 26 (2012) 447–454. http://dx.doi.org/10.1111/j.1365-2915.2012.01018.x.

[7] J. Rosenfeld, J. Capdevielle, J.C. Guillemot, P. Ferrara, In-gel digestion of proteins for internal sequence analysis after one- or two-dimensional gel electrophoresis, *Anal. Biochem.* 203 (1992) 173–179.

[8] H. Renault, V. Roussel, A.E. Amrani, M. Arzel, D. Renault, A. Bouchereau, et al., The Arabidopsis pop2-1 mutant reveals the involvement of GABA transaminase in salt stress tolerance, *BMC Plant Biol.* 10 (20) (2010), http://dx.doi.org/10.1186/1471-2229-10-20.