Next Article in Journal
Effects of Elevated CO2 Concentration on Host Adaptability and Chlorantraniliprole Susceptibility in Spodoptera frugiperda
Previous Article in Journal
Transcriptomic Identification and Expression Profile Analysis of Odorant-Degrading Enzymes from the Asian Corn Borer Moth, Ostrinia furnacalis
Previous Article in Special Issue
Temporal Evaluation of Insecticide Resistance in Populations of the Major Arboviral Vector Aedes Aegypti from Northern Nigeria
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Vivo Glutathione S-Transferases Superfamily Proteome Analysis: An Insight into Aedes albopictus Mosquitoes upon Acute Xenobiotic Challenges

by
Siti Nasuha Hamzah
1,*,
Silas Wintuma Avicor
2,
Zazali Alias
3,
Sarah Abdul Razak
3,
Siti Khadijah Mohd Bakhori
4,
Ting Chuan Hsieh
1,
Nurin Nazifa Syanizam
1 and
Salinah Abdul Farouk
1
1
School of Biological Sciences, Universiti Sains Malaysia, George Town 11800, Penang, Malaysia
2
Entomology Division, Cocoa Research Institute of Ghana, New Tafo-Akim P.O. Box 8, Ghana
3
Institute of Biological Sciences, Faculty of Science, Universiti Malaya, Kuala Lumpur 50603, Malaysia
4
School of Physics, Universiti Sains Malaysia, George Town 11800, Penang, Malaysia
*
Author to whom correspondence should be addressed.
Insects 2022, 13(11), 1028; https://doi.org/10.3390/insects13111028
Submission received: 20 August 2022 / Revised: 29 October 2022 / Accepted: 1 November 2022 / Published: 7 November 2022
(This article belongs to the Collection Pesticide Chemistry and Toxicology)

Abstract

:

Simple Summary

Glutathione S-transferases (GSTs) are a superfamily of enzymes that play a crucial part in phase II detoxification reactions. In many cases, GSTs have been implicated in toxicological challenges by association only (i.e., detected using model substrates or due to an increase in GST activity) and not through the identification of the individual peptides involved. In this study, the aim was to identify the GSTs induced by acute insecticide challenges up to the isoform or peptide level and the detoxification mechanisms involved in Aedes albopictus larvae. Upon permethrin and malathion inductions, the GST-based peptides were identified, and the functional characteristics as well as protein-protein interactions were analyzed along with the quantification of the abundance. Twenty-four GST peptide isoforms were identified under seven different classes with a highly significant protein-protein interaction between GST groups and within other detoxification metabolism-related peptides (p-value < 1.0 × 10−16). Overall, our findings indicate that out of seven GST classes, only Delta and Epsilon GSTs were mainly involved in the detoxification mechanism via the direct glutathione metabolism or the sequestration of the insecticides, as well as providing the protection of sulfhydryl groups in enzymes against oxidative stress induced by insecticide exposure.

Abstract

In this study, the induction of glutathione S-transferase (GST) enzymatic activities in Aedes albopictus under 24 h of xenobiotic challenges was investigated. From LCMS analysis, 23 GST isoforms were identified under Delta, Epsilon, Sigma, Zeta, Omega, and Iota classes, together with one GSTX1-1 isoform, in both treated and untreated samples. Using STRING 11.5, the functional enrichment network of Gene Ontology (GO) analysis, the identified peptides were found to be involved in the glutathione metabolic biological process (GO:0006749, p-value: 1.93 × 10−29), and the molecular functions involved are due to glutathione transferase (GO:0016848, p-value: 2.92 × 10−8) aside from carbon-halide lyase activity (GO:004364, p-value: 1.21 × 10−31). The Protein-Protein Interaction (PPI) network (STRING 11.5) showed significant interactions within the GST superfamily and some of the GST classes interacted with other proteins among the input domain of the identified peptides (p-value < 1.0 × 10−16). In TMT labeling for the quantification of peptide abundance, isoforms from Delta (GSTD1-2, GSTD1-3, GSTD1-4) and Epsilon (GSTE3-1, GSTE4-2) were found to be overexpressed (between 1.5-fold and 2-fold changes). In the PPI analysis, 12 common enriched pathways of Kyoto Encyclopedia of Genes and Genomes (KEGG) were found to be intercorrelated with the identified GSTs at PPI enrichment p-value < 1.0 × 10−16. Overall, this study indicates that distinct GST enzymes, which were identified up to their specific protein isoforms, are involved in the metabolic mechanisms underlying xenobiotic stress.

1. Introduction

Glutathione S-transferases (GSTs: E.C.2.5.1.18) represent a major group of detoxification enzymes, comprised of numerous hydrophobic toxic compounds, and are multifunctional dimeric proteins crucial in drug metabolism, xenobiotic and endobiotic biotransformation, intracellular transport, the biosynthesis of hormones, and the metabolism of xenobiotics [1,2,3,4,5]. Among these, GSTs play a predominant role in providing protection against electrophiles and oxidative stress products [5].
The mammalian GSTs have been categorized into eight classes which are Alpha, Mu, Pi, Theta, Sigma, Zeta, Kappa, and the microsomal class [5,6]. Specific interest has been shown towards GSTs due to their role in metabolic enzyme-based insecticide resistance as various insects have demonstrated elevated GST enzymatic activity [6,7]. Studies have emphasized the action of insect GSTs in insecticide resistance, particularly dichlorodiphenyltrichloroethane (DDT) and pyrethroid resistance [2,4,8,9,10]. The role of GSTs in the dehydrochlorination of DDT is the primary route of its detoxification and serves as the most common DDT resistance mechanism in mosquitoes [11,12].
A total of six cytosolic GST classes have been identified in insects including mosquitoes which are Delta, Sigma, Epsilon, Omega, Theta, and Zeta [13,14,15,16]. Among these classes, the Delta and Epsilon classes are apparently the primary classes involved in the metabolic detoxification of xenobiotics and have been identified to be responsible for the propagation of insecticide resistance [6,10,13,15]. Penilla et al. [17] reported that there is an increase in GST activities under continuous pyrethroid selection in field strains of Anopheles albimanus. The metabolic resistance to DDT in Anopheles spp. larvae and adults is primarily conferred by an increase in GST activity [15,18].
The species of interest in this study is Aedes albopictus, which acts as the secondary dengue vector in Malaysia as well as being the pathogen carrier of various other viruses [19]. Dengue poses a serious health threat in Malaysia, and the study of the species is, therefore, pertinent to reducing the burden of the disease. Aedes mosquito management is the sole strategy that can be used since a stable dengue vaccine is yet to be established [20,21]. Utilizing pesticides, eliminating larval breeding grounds, and using spatial repellents are the primary methods for controlling Ae. albopictus [22,23,24]. Unfortunately, mosquitoes metabolize chemicals using their defense mechanisms as a result of the continuous, intensive, and broad application of insecticides, which eventually results in resistance to the insecticides, a rapidly emerging phenomenon [22,25]. The selection pressure of insecticides and the inheritability of resistance among generations of vector populations are some other aspects that undermine the efforts to control vectors [10,26].
Insects evolve xenobiotic resistance primarily through improving the metabolic efficiency of their detoxification systems and/or decreasing the sensitivity of the xenobiotic target site [26]. Metabolic-based resistance occurs when the amount of insecticide that reaches the target site is reduced due to an enhanced rate of xenobiotic metabolism caused by the actions of enzymes [27]. Resistance involving the enzymatic activities of metabolic enzymes is of concern due to the fact that the enzymatic activities are not constant and differ according to different strains from different locations which are caused by various factors. For integrated vector control programs to be successful yet sustainable, it is essential to understand the underlying mechanisms of pesticide resistance [24]. A biochemical assay provides information regarding the total amount of protein and the catalytic behaviors of the target enzyme. Proteomic analysis on the other hand enables the identification and quantification of the different expression levels of the enzymes. Considering insecticides from the pyrethroid and organophosphate groups are frequently employed in programs to control mosquitoes, the intended outcomes of this study were to investigate the effect of 24 h of acute insecticide treatments on Ae. albopictus in terms of GST enzymatic activities in order to identify the toxicological challenge mechanisms in Ae. albopictus larvae. Subsequently, the different classes and isoforms of GST enzyme proteomes, as well as their expression levels based on their peptide abundance upon the acute 24 h treatment with malathion and permethrin, were observed.

2. Materials and Methods

2.1. Mosquito Strains

An established 284th generation of a susceptible strain of Ae. albopictus was obtained from the Vector Control Research Unit (VCRU), Universiti Sains Malaysia. The 4th instar larvae of the susceptible Ae. albopictus were used as the acute insecticide-treated strains and the untreated control strains.

2.2. Determination of Toxicity Parameters and Insecticide Treatment

The susceptible VCRU strain was subjected to the standard World Health Organization (WHO) larval bioassay procedure [28,29] to determine the sub-lethal concentrations (LC50) of malathion (96.0% purity) and permethrin (96.9% purity). Batches of 25 early 4th instar larvae of the laboratory susceptible strain were transferred to a 500 mL plastic cup containing a mixture of seasoned water and insecticide solution with a final volume of 250 mL. Ethanol was used as the solvent to dilute the insecticides. Seasoned water containing 1% ethanol and larvae summing up to 250 mL was used as control. A range of concentrations in four replicates against the yielded mortality ranging from 0% to 100% within 24 h was plotted using Probit analysis in SPSS software version 24 [30] to obtain the LC50 value [29]. For the acute treatment, the susceptible larvae were exposed to the determined malathion and permethrin LC50 for a period of 24 h. Every treated larva that survived was collected and used in further analysis.

2.3. Protein Determination and Enzyme Assays

The early 4th instar larvae of the susceptible strains and larvae which survived acute treatment with the insecticides for 24 h were individually homogenized in 200 μL of seasoned water on ice and centrifuged at 14,000 rpm, 4 °C for 30 s, and the supernatant was used as an enzyme source in the assays. At least 100 individual replicates were used in the assays. The assays were performed in a 96-well microplate on ice and the absorbance (optical density [OD] values) was measured on the microtiter plate reader (Thermo Fischer Scientific) with Magellan data analysis software. The total protein content (mg) was determined using the Bradford protein assay with bovine serum albumin (BSA) as a standard to normalize the activities for protein concentration [31]. The GST activities toward the model substrates were conducted according to Hemingway [32] and Hamzah and Alias [33]. The GST conjugation activities with 1-Chloro-2,4-dinitrobenzene were determined at 340 nm. The specific activity was calculated in μmol/min/mg protein by following Beer’s Law, A= εcl where A is the absorbance; ε is the extinction coefficient; c = concentration; and l is the path length. Then the total enzyme activity was expressed as μmol/min at 25 °C [34].

2.4. Trypsin Digestion and Tandem Mass Tag (TMT) Labelling

A tandem mass tag (TMT) was used as a chemical label for mass spectrometry (MS)-based identification and to quantify the identified peptides. Briefly, the supernatant of homogenized susceptible, acute permethrin-treated and malathion-treated larvae were lyophilized and re-suspended in a Radioimmunoprecipitation assay (RIPA) buffer. The samples were concentrated by precipitation using the Total Protein Precipitation Kit (ITSIBIO, USA). Each sample was assayed for protein content using the total protein assay kit (ITSI Biosciences). The dried protein pellet was re-suspended in a lysis buffer. An equal amount of protein (20–100 µg) from each sample was reduced with tris(2-carboxyethyl) phosphine (TCEP) and alkylated with iodoacetamide. The precipitated samples were then re-suspended in 100 mM triethylammonium bicarbonate (TEAB), and trypsin digestion was performed overnight. The digested samples were individually labeled with 3 different TMT reagents (6 plex amine reactive tags kit) according to the manufacturer’s protocol (Thermo Scientific). The fractions were desalted by using C18 zip tips and dried down.

2.5. Liquid Chromatography Mass Spectrometry

Liquid Chromatography-Mass Spectrometry (LCMS) was performed according to the manufacturer’s manual (Thermo Fisher) by using an Orbitrap Fusion mass spectrometer (Thermo Fisher). It was equipped with a 100 µ × 20 mm C18 100 Å 5U reverse-phase trap for online desalting and a 75 µ × 150 mm C18 200 Å 3U reverse-phase column for peptide separation. Subsequently, 1 μL of the peptide mix was loaded onto the reverse-phase column at an isocratic flow rate of 300 mL/min and a gradient of 0.1% formic acid (A) and 100% acetonitrile (B). The eluents gradient used was 5% to 35% B for 60 min. The high-mass accuracy MS full spectra data were obtained in the data-dependent mode with the Orbitrap MS (OTMS). The resolution was set to 120,000 at MS1 with a scan range of 310–1800 m/z, an automatic gain control (AGC) target of 400,000, and a maximum injection time of 50 milliseconds. The method consisted of a 3 s top-speed mode where precursors were selected for a maximum 3 s cycle. Only precursors with an assigned monoisotopic m/z and a charge state of 2–7 were further analyzed for MS2. All precursors were filtered using a 20 s dynamic exclusion window and an intensity threshold of 5000. The MS2 spectra were analyzed with an ion trap MS (ITMS) set with a rapid scan rate at a resolving power of 60,000, an AGC target of 100, a 1.6 m/z isolation window, and a maximum injection time of 250 milliseconds. The precursors were fragmented by collision-induced dissociation (CID) and high-energy collision dissociation (HCD) at a normalized collision energy of 30% and 28%, respectively. The top 10 fragmented ions from MS2 were fragmented again to produce the MS3 spectrum to increase the TMT reporter ion population. Every MS3 precursor population was fragmented by HCD at 65% collision energy.

2.6. Analytical Statistics

The data obtained in WHO bioassays were statistically analyzed by using the Probit analysis computer program SPSS version 24 [30] to determine LC50 values [29]. Adjustments were made using Abbott’s formula [34] if the control mortality ranged from 5% to 20%. If the pupation exceeded 10%, the test was rejected [29]. The LC50s of permethrin and malathion obtained for the control strain were used further as the dose for the acute-insecticide-treatment effect study. The mortality percentages and enzymatic activities of GST were tested for normality as well as variance homogeneity by using Komolgorov–Smirnov and Levene’s tests, respectively. The GST enzymatic activity data were presented as mean ± standard deviation (S.D.). One sample t-test (T-test) was applied to compare the total protein content and particular enzyme expression levels between each acute insecticide-treated strain with the susceptible strain only.

2.7. Querying Domain and Interaction Information

For proteomic analysis, the GST raw data files were searched against the most recent database for Aedes downloaded from Uniprot by using the MudPIT option in Proteome Discoverer 2.1 (Thermo Scientific) and the Sequest HT search algorithm. Only peptides identified with high confidence were considered and used for protein identification. Trypsin was the selected enzyme allowing for up to two missed cleavages per peptide. Carbamidomethylation of cysteine, N-terminal TMT3-plex, and Lysine TMT3-plex were used as a static modification whereas the oxidation of Methionine and the TMT3-plex of Threonine were used as variable modifications. The group peptide abundances were generated using Proteome Discoverer 2.1 software coupled with the Sequest HT search algorithm. After the Peptide and Protein Quantifier node classified the peptide groups, it calculated protein abundances for the samples as the simple summation of its associated peptide group abundances. Precursor ion area detection presented the average peak areas of the top N unique peptides for each protein, usually set to 3 by default. The node creates the Abundance Counts column on the Proteins page, which shows the number of peptide group abundances used for calculating abundances. Finally, the node calculated the ratios (fold-change) by dividing the abundance values of the associated sample groups (susceptible, malathion-treated, and permethrin-treated samples). Thus, a TMT ratio ≥ 1.2 was classified as upregulated, and ≤0.80 was classified as downregulated. A fold-change from 0.81 to 1.19 was considered a moderate-to-no change. Ratios above 100 exceeded the maximum allowed threshold that was reasonably expected.

2.8. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and Protein-Protein Interaction (PPI) Analysis

Uniprot (https://www.uniprot.org/ accessed on 24 May 2022) was used to obtain information on the identified protein domains. Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Protein-Protein Interaction (PPI) network were integrated using STRING analyses. The interaction between the identified proteins was determined by using STRING 11.5 at an interaction coefficient of 0.4 as the minimum standard. Protein interactions with confidence coefficients larger than 0.4 (medium confidence) were regarded as relevant because it was the most often employed analytic criterion in STRING. Twenty-four identified protein accessions generated from the MudPIT option in Proteome Discoverer 2.1 (Section 2.7) were used in the STRING analysis. The accessions are J9HXZ8, J9HHL7, Q17MA9, A0A1S4EXN6, Q6PTY1, A0A0N8ES64, Q17MB8, Q16SH6, A0A023ENG1, Q16SH7, Q170C6, Q5PY78, Q170C9, A0A0P6IV26, Q170C7, Q16P53, J9E9C0, A0A1S4G560, A0A0P6J0T5, Q1HQK1, Q16P79, Q16P80, Q0C791, and Q16NL9. The biological information contained in the genome, such as the biological processes, molecular data, and cellular components, was interpreted using Gene Ontology (GO). By connecting the function of gene products with metabolic pathways, the Kyoto Encyclopedia of Genes and Genomes (KEGG) enables the understanding of biological systems at the molecular level.

3. Results

3.1. Identification of GST Isoforms

This study aimed to identify and analyze the differential expression of each isoform or class of GST enzymes. The detection of all targeted enzymes was conducted by blasting the mass spectrometry data obtained with the Aedes genus database. The data obtained were further analyzed using UniProtKB up to the GST enzyme classes and each isoform was designated as its proposed identification label for easy analysis and reference. The identification of GSTs was carried out according to the International Union of Biochemistry and Molecular Biology (IUBMB). The details of the isoforms identified are summarized and presented in Table 1. In total, 24 GST isoforms were successfully identified in up to 6 different classes, which were Delta, Epsilon, Sigma, Zeta, Omega, and Iota. However, from the analysis, one unidentified gene/isoform which did not belong to any of the current GST classes was detected and designated as GSTX1-1 (24.817 kDa). For the Delta class, four isoforms of the Delta1 class (GSTD1-1 to GSTD1-4) and three isoforms of the Delta4 class (GSTD4-1 to GSTD4-3) with a molecular weight ranging from 23.146 kDa to 24.855 kDa were identified. Two isoforms of the Delta6 class (GSTD6-1 and GSTD6-2) and one isoform of the Delta11 class (GSTD11-1) were also discovered with a higher molecular weight of 28.215 kDa, 28.148 kDa, and 25.879 kDa, respectively.
For the classification of Epsilon, one of Epsilon3 (GSTE3-1), two isoforms of Epsilon4 (GSTE4-1 and GSTE4-2), one isoform of Epsilon5 (GSTE5-1), and one isoform of Epsilon6 (GSTE6-1) were discovered with a molecular weight range from 24.635 kDa to 27.264 kDa. Only one isoform of the Iota class was identified with a molecular weight of 26.177 kDa and was designated as GSTI-1. For the Omega class, only the Omega1 class was detected with four isoforms designated as GSTO1-1 to GSTO1-4 with a molecular weight ranging from 28.551 kDa to 29.618 kDa. A similar pattern was observed for the Sigma class whereby only the Sigma1 class was found with two isoforms (GSTS1-1 and GSTS1-2) and exhibited a molecular weight of 23.249 kDa and 23.216 kDa. One isoform of the Zeta1 class was identified and designated as GSTZ1-1 (26.39 kDa). Overall, a close range of 23.146 kDa to 29.618 kDa in molecular weight as well as a pI value of 5.17 to 8.51 was observed among all the classes of GST detected in the in vivo homogenized samples.

3.2. Protein-Protein Interaction (PPI) Network Analysis of GST Families

To predict the interacting proteins and the outcomes of the in vivo biological analysis, Ae. albopictus GST protein in the listed classes (Table 1) was used as the model and applied to the STRING 11.5 tool. Based on the results from 13 GSTs outputs, 10 functional partners were identified in the network analyses for GSTs, listed in Table 2 with their scores and functions. The computerized output of the 13 GSTs was reflected in the in vivo analysis earlier.
Furthermore, the enriched pathways of the biological process, the molecular function, and the KEGG obtained by PPI analysis for the identified GSTs are listed in Table 3 and Table 4 with their p-values. The K-means algorithm was used for protein clustering in 13 different colors based on the different GST classes that contributed to the analysis (Figure 1). The result of the biological processes and molecular functioning in GO showed that the protein peptides are involved in several processes. The results obtained a total of 47 terms, containing 27 biological processes and 20 molecular functions (Table 3). The glutathione metabolic process and the peptide metabolic process had the higher enrichment level among the biological processes. In terms of molecular functions, the peptides were mainly associated with Glutathione transferase activity. Catalytic activity, Carbon-halide lyase activity, Glutathione hydrolase activity, and Peptidyltransferase activity also showed a correlation between the molecular functions and the identified isoforms corresponding to the significant p-value (the smaller the p, the higher the enrichment level).
The results of the KEGG analysis of GST proteins identified 12 pathways involved in all 13 GST classes (Table 4). The enrichment results showed that the glutathione metabolism and metabolic pathways were the highest in the KEGG analysis. However, the GSTs were also associated with the drug metabolism of cytochrome P450 and other enzymes as well as the arachidonic acid metabolism. On the basis of the GO and KEGG results, our initial hypothesis was that the GSTs play an important role in xenobiotic metabolism and interacted with other enzymes in the molecule to maximize the metabolism process in the organisms.

3.3. Effect of Acute Insecticide Treatment on GSTs Activities

The present study used permethrin (pyrethroid) and malathion (organophosphate), which are from the common insecticide classes used in mosquito control strategies. Table 5 shows the calculated sub-lethal concentration (LC50) values obtained, which were used in acute insecticide treatments. It indicated that permethrin was more effective against Ae. albopictus larvae followed by malathion insecticide with LC50 values of 0.023 mg/L and 0.099 mg/L, respectively. The determined LC50s were used to investigate the effect of acute insecticide treatment on GST activities and moreover on the change in the abundance of each GST peptide. Table 6 shows the result of a 1-sample t-test of the total protein content (µg) of the larvae upon 24 h of acute treatment with malathion and permethrin. The total protein content of the larvae exposed to 24 h of acute permethrin and malathion treatments increased significantly (p < 0.05) after 24 h, at 5.3 × 10−3 ± 0.006 mg (1.36-fold) and 4.3 × 10−3 ± 0.002 mg (1.10-fold), respectively, compared to the value of the susceptible strain. The total GST activity in the larvae after 24 h of permethrin and malathion treatment increased up to 13.97 × 10−3 ± 0.01 μmol/min (1.84-fold) and 10.75 × 10−3 ± 0.02 μmol/min (1.41-fold), respectively, and were found to be significant (p < 0.05) upon comparison to the control strain. Similarly, the specific activity of GSTs was 2.635 ± 0.08 µmol/min/mg (1.35-fold) and 2.502 ± 0.050 µmol/min/mg (1.29-fold) for larvae under 24 h of permethrin and malathion treatments, respectively, compared to the control strain (p < 0.05). Conclusively, a varying increment in total GST enzyme activities and specific GST activities was detected after a 24 h period of acute malathion and permethrin treatments.

3.4. Differential Analysis of the Peptide Abundance under Xenobiotic Challenge

The fold-change in the peptide abundance of various isoforms of GSTs is summarized in Table 7. The elevated peptide abundance fold-change was detected from isoforms GSTD1-2, GSTD1-3, and GSTD1-4, with a fold range of 1.92-fold to 2.26-fold, followed by GSTE3-1 and GSTE4-2 under permethrin and malathion treatments. Meanwhile, the fold-change in the peptide abundance of GSTE5-1 and GSTE6-1 after 24 h of acute permethrin treatment decreased by 0.66-fold each. Only GSTS1-2 was found to decrease the peptide abundance fold-change (0.75-fold) under malathion treatment.

4. Discussion

4.1. Effect of the GST Enzymatic Activities of Aedes albopictus Larvae under 24 h of Acute Insecticide Challenges

The acute treatment of early 4th instar larvae treated for 24 h with sub-lethal concentrations of permethrin and malathion was conducted to observe the induction effect of the enzyme’s activities in the presence of toxicological challenges. Sub-lethal doses were selected for this study and defined as acute treatment because this dose has been tested to yield an optimum level of 50% of surviving larvae to be used for further analysis. A pilot study with lesser doses yielded less satisfactory enzymatic activities while higher doses resulted in a very limited number of surviving larvae and enzymatic activities which were lower than the control strain.
In insects, the metabolism of xenobiotics is one of the predominant insecticide detoxification mechanisms [35]. The biosynthesis of enzymes is suggested to occur because of a direct reaction upon exposure to xenobiotics by enzymes that hydrolyze, oxidize, and conjugate xenobiotics into less toxic, more water-soluble products thereby causing the products to be easily excreted [36,37]. The acute insecticide treatment of mosquitoes induces different levels of metabolic enzyme production which cumulatively increase the overall total protein content in the mosquito. The overall total protein fold-change of larvae post-acute malathion and permethrin treatments showed a significant increment (p < 0.05) upon comparison with the control strain, respectively. It has been suggested that the increase in protein content may have an impact on the increased activity of enzymes under toxicant stress, which is connected to an oxidative stress state in induced larvae. According to Devonshire et al. [38], metabolic defense mechanisms can be accelerated by the overproduction of GST content followed by an increase in its activity. Moreover, in Yan et al. [39], GSTs have been reported to be present in high concentrations in Apis cerana after being challenged by an extreme oxidative stress agent. The introduction of oxidative stress agents, namely insecticides such as permethrin and malathion, may lead to antioxidant defense by the initiation of the repairing activity of the damaged secondary product generated by reactive oxygen species (ROS) [33,40,41]. The determination of the total and specific enzyme activities is crucial to observe the effect of acute insecticide exposure on the rate of the metabolic enzyme activities of the test mosquitoes. In samples treated with acute malathion for 24 h, the difference in the fold between the total and specific enzyme activities of GSTs was found to be significantly elevated (p < 0.05) compared to the control strain. The elevation of the fold-change in the enzymatic activity of specific GSTs upon acute malathion treatment for 24 h is justifiable because it has a wide range of substrate specificity to different insecticide classes including organophosphates [42]. Upon comparison to the control strain, the elevation of the total and specific enzyme activity of GSTs (p < 0.05) was detected in the susceptible larvae after a 24 h acute permethrin treatment. The same pattern has been recorded in previous studies whereby there was a considerable increment in the enzymatic activities of GSTs in the field strain as well as the laboratory-selected Ae. aegypti mosquitoes, which implied its possible involvement in pyrethroid metabolism and resistance [43,44,45]. Apart from that, the results obtained are in line with several earlier studies on the involvement of GST enzyme activity in insecticide resistance including pyrethroid resistance [2,4,9,10]. Exposure to insecticides activates the enzymatic activities of metabolic enzymes. Upon prolonged exposure to insecticides, the enzyme production increases for the mosquito to detoxify more insecticides to survive despite the exposure to the toxicity of the insecticides. These populations of surviving mosquitoes will eventually propagate and become resistant to insecticides and hamper vector control programs.

4.2. Proteomic Study of Aedes albopictus Larvae upon 24 h of Acute Insecticide Treatments

Proteomics is regarded as the detailed study of proteomes which are the protein complements of the genome [46]. This approach allows for a more comprehensive study of the structure of proteins, while elucidating the diverse functions of, as well as the interactions between proteins [47]. This branch of study has also been used to analyze mosquitoes in an effort to identify and characterize the expression of the proteins of various species, tissues, and physiological stages, as well as the response to parasites, environmental conditions, and xenobiotic challenges [47,48].
Even though elevations of induced GST enzymatic activities were detected in the larvae upon 24 h of acute permethrin and malathion treatments, it was imperative to identify the presence and peptide abundance of the specific GST isoforms to determine their involvement in the detoxification of xenobiotics. Hence, a proteomic study was conducted using the crude homogenized Ae. albopictus larvae which had been subjected to 24 h of acute insecticide treatments. As the identification and the peptide abundance detection were completed up to the peptide and amino acid level, this proteomic data verified the GST enzymatic activity data from the previous sections.

4.3. Identification of GST Enzyme Classes up to the Isoform and Their Response towards Acute Insecticide Treatments

A broad range of classes and isoforms of GST enzymes were obtained after further analysis using UniProtKB. The target detoxification enzyme, GST, is the enzyme substantially responsible for the metabolism of all four major classes of insecticides [49,50,51]. The capability of GST enzymes to metabolize insecticides is of huge concern due to the fact that the elevated rate of their enzymatic activities has a tendency to result in the emergence of tolerance and eventually resistance to insecticides which threatens vector control. From the analysis, different GST isoforms were identified in this study from the control and the acute permethrin and malathion treatments, which validated their presence in the larvae of Ae. albopictus.
The classification of GSTs of the same classes was made if there was more than 40% identity of amino acid sequence and other properties, namely immunological properties, phylogenetic criteria, tertiary structure, chromosomal locations, and the capacity to form heterodimers [1,3,5,15]. To the best of our knowledge, only studies by Hamzah and Alias [16,33] and Hamzah et al. [52] characterized, identified, and elucidated the effects of chemical challenges on GST enzymes in Ae. albopictus. In the study of Hamzah et al. [52], the protein expression was determined and quantified by analyzing the peak of the gel-expressed for each visible GST spot. The TMT labeling MS-based approach was used in the present work as an essential step to increase the number of the identified proteins. As a result, out of the six classes of insect GSTs, five were successfully detected in this study and they were identified up to their specific isoforms, and one isoform from the class Iota as well as one GSTX1 were also identified.
In this study, via protein-protein interaction analysis using STRING 11.5, all 7 GST classes were further annotated into 13 domains as shown in Figure 2. The network shows more significant intra-network interactions than expected. This means that the proteins have more interactions among themselves than what would be expected for a random set of proteins of the same size and degree distribution drawn from the genome. Such an enrichment indicates that the proteins are at least partially biologically connected as a group. In the functional enrichment network, from the Gene Ontology analysis, the 13 domains were identified under the glutathione metabolic biological process (GO:0006749, p-value: 1.93 × 10−29), and the molecular functions involved are due to glutathione transferase (GO: 0016848, p-value: 2.92 × 10−8) and carbon-halide lyase activity (GO:004364, p-value: 1.21 × 10−31). The domains were located in the subcellular cytoplasm (GO:0005737, p-value: 0.00015) and highly significantly involved in the glutathione metabolism (p-value: 5.39 × 10−25), metabolic pathways (p-value: 1.29 × 10−14), drug metabolism corresponding with other enzymes (p-value: 2.08 × 10−19), drug metabolism corresponding with cytochrome P450 (p-value: 8.55 × 10−24), and the metabolism of xenobiotics by cytochrome P450 (p-value: 8.55 × 10−24) with the KEGG Pathway IDs aag00480, aag01100, aag00983, aag00980, and aag00980, respectively.
The analysis further showed that the xenobiotic metabolism involved interactions not only among the domains but the interaction enrichment involved a total of 33 proteins around the input domains (Figure 1). In the analysis, STRING 11.5 then further added (by default) other occurring proteins around the network that showed expected connectivity enrichment with an average local clustering coefficient of 0.779 and a PPI enrichment p-value < 1.0 × 10−16. Generally, proteins do not work independently; this discovery lends credence to the basic proposition that proteins mediate all biological processes through interactions with other proteins, the formation of protein complexes, and the development of metabolic and signaling pathways [53,54]. Networks in biology can be established using PPI data. The special and significant nodes in the biological network may be located and the relationships between these nodes can be further understood by a thorough investigation of the network [55,56].
Furthermore, of all the isoforms identified in the samples that underwent 24 h of acute permethrin treatment, the peptide abundance of isoforms from Delta class GSTs (GSTD1-2, GSTD1-3, and GSTD1-4) were overexpressed (more than 2-fold), followed by peptide abundance of isoforms from Epsilon class GSTs (GSTE4-2). A similar but lower fold-change in the peptide abundance of isoforms was also detected in the samples that underwent 24 h of acute malathion treatment, with the highest overexpression in the peptide abundance in isoforms from the GST Delta class (GSTD1-2, GSTD1-3, and GSTD1-4 with a GO term of GO:0006749) followed by isoforms from Epsilon class which are GSTE3-1 (GO:0006749) and GSTE4-2 (GO:0016740; GO:0006749). These results showed a major involvement of Delta and Epsilon GSTs in the metabolism of insecticides. The information regarding the contribution of insect-specific Epsilon and Delta GSTs, which are the largest classes of cytosolic GSTs to insecticide detoxifications as well as resistance, has been well documented [7,13,57,58]. From their GO terms, these five isoforms have a biological and molecular function that involves the glutathione metabolic process and transferase activity. As in the chemical structure, permethrin is a synthetic insecticide that contains a cyclopropane ring which is substituted with a 2,2-dichlorovinyl and gem-dimethyl groups while malathion contains phosphate groups in the chemical structure (PubChem, NCBI). GSTE4-2 was found to be involved in the molecular function of transferase activity where the mode of action to counter the toxicological effect is by the catalytic activity of the methyl group in permethrin and the phosphate group in malathion from the donor to the acceptor to produce water-soluble metabolites that are more readily excreted (GO:0016740). Moreover, the main mode of action is by performing the glutathione metabolic biological process (GO:0006749) which acts as an antioxidant in the protection of sulfhydryl groups in enzymes and other important functional proteins, hence it could also be protective against substances such as permethrin and malathion which are known as oxidative stress inducers [52].
Moreover, overexpression in the peptide abundance of different GST isoforms suggests the involvement and induction of these isoforms upon a 24 h acute permethrin and malathion treatment and validates the outcomes attained from the biochemical assay, whereby the fold-change of the total as well as specific enzymatic activity of GSTs were also elevated in larvae that underwent 24 h of acute permethrin and malathion treatment. The GST, as the major phase II detoxification enzyme, plays a main role as a defense enzyme in the mosquito by catalyzing the conjugation process of GSH to the exposed toxic insecticides [13]. Apart from that, there have been reports about the direct involvement of GSTs in oxidative stress and in one of the mechanisms in the metabolism of insecticides [39,57,59,60]. Hence, GST overexpression would definitely increase tolerance to cytotoxicity [1].
Upon 24 h of acute permethrin treatment, the Zeta class of GSTs (GSTZ1-1) exhibited an overexpression in the fold-change of peptide abundances which implies its involvement in detoxification activities. Blackburn et al. [61] reported that a deficiency of glutathione transferase Zeta causes oxidative stress and the activation of antioxidant response pathways. There was an overall reduction in the fold-change of the peptide abundance from Sigma class (GSTS1-1 and GSTS1-2) for both the larvae that underwent 24 h of acute permethrin treatment and those that underwent 24 h of acute malathion treatment. Omega class GSTs (isoforms: GSTO1-1 to GSTO1-4) from both the larvae that underwent 24 h of acute permethrin treatment and those that underwent 24 h of acute malathion treatment showed a fold change anywhere from negligible to slightly higher compared to the control. The pattern of fold-changes suggests that these classes of GST enzymes (Sigma and Omega) did not contribute to the detoxification activities. In Figure 1, the GSTS1, GSTO1, and GSTZ1 nodes interacted significantly with proteins outside the protein group (the non-GST-related proteins). It should be noted that the involvement of the Omega and Sigma classes of GSTs in housekeeping has also been reported. A few studies indicate that apart from the Omega and Sigma classes of GST, the Zeta class also has a relatively broad taxonomic distribution and serves an essential role in housekeeping, and this information, therefore, dismisses its involvement in the metabolism of xenobiotics [62,63]. Collectively, from this study, the involvement of different classes and isoforms of GST enzymes was detected and identified in Ae. albopictus larvae upon 24 h of acute permethrin and malathion treatment. This study enables the pinpointing of GST enzymes that potentially underpin the development of resistance in response to pesticide challenges. Invariably, this will help in the early detection and prediction of resistance development in field strains of the vector.

5. Conclusions

The present study provides some biochemical and proteomic information on the targeted detoxification enzymes (GSTs) in susceptible insecticide-treated strains. Furthermore, the total protein content, total enzymatic activity, and specific enzymatic activity of GSTs were found to be significantly elevated after 24 h of acute permethrin and malathion treatment on larvae. In the proteomic study, the occurrence of multiple isoforms of GSTs in samples that underwent 24 h of acute insecticide treatment with varying fold-change in their peptide abundance was identified. Various isoforms from five classes (Delta, Epsilon, Omega, Sigma, Iota, and Zeta) as well as one unidentified GST, summing up to 24 GST isoforms, were identified. Cumulatively, the results of the peptide abundance quantitative study reflected the results of the enzymatic activities from the biochemical analysis. The in vivo peptide abundance correlated with its functional analysis (p-value < 1.0 × 10−16), suggesting that the mechanism of specific insecticide (permethrin and malathion) detoxification involved specific GST isoenzymes (Delta and Epsilon GSTs) via specific metabolic strategies (direct glutathione metabolism, the sequestration of the insecticides, or providing protection for sulfhydryl groups). The results obtained from this research provide more knowledge about the current resistance status and enzymatic activity of GST enzymes in Ae. albopictus mosquitoes upon acute and xenobiotic challenges. Knowledge of the resistance status and GST enzymatic activity in the vector is key to understanding the patterns of GST-mediated resistance for effective resistance management. The differentially expressed proteins could be used for monitoring and surveillance purposes to detect incipient pesticide resistance for early management in Ae. albopictus.

Author Contributions

Conceptualization: S.N.H., S.W.A. and Z.A.; Data curation, S.N.H., S.K.M.B. and S.A.F.; Formal analysis, S.N.H., S.K.M.B. and S.A.R.; Funding acquisition: S.N.H.; Investigation, T.C.H., N.N.S. and S.A.F.; Methodology, S.N.H., S.W.A., S.A.R., S.K.M.B. and S.A.F.; Project administration, S.N.H.; Resources, Z.A.; Supervision, S.N.H., S.W.A. and Z.A.; Validation, S.A.R., S.K.M.B., T.C.H. and N.N.S., Visualization, S.A.R.; Writing—original draft, S.N.H. and S.A.F.; Writing—review and editing, S.W.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundamental Research Grant Scheme (FRGS), Ministry of Higher Education Malaysia. Grant number FRGS/1/2018/STG04/USM/02/7 and account code 203.PBIOLOGI.6711639.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Data is contained within the article except spectra data which is available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ketterman, A.J.; Saiwasang, C.; Wongsantichon, J. Insect glutathione S-transferases. Drug Metab. Rev. 2011, 43, 253–265. [Google Scholar] [CrossRef] [PubMed]
  2. Enayati, A.A.; Ranson, H.; Hemingway, J. Insect glutathione s-transferases and insecticide resistance. Insect Mol. Biol. 2005, 14, 3–8. [Google Scholar] [CrossRef] [Green Version]
  3. Ranson, H.; Hemingway, J. Mosquito glutathione s-transferases. Methods Enzymol. 2005, 401, 226–241. [Google Scholar]
  4. Vontas, J.G.; Small, G.J.; Hemingway, J. Glutathione S-transferases as antioxidant defence agents confer pyrethroid resistance in Nilaparvata Lugens. Biochem. J. 2001, 357, 65–72. [Google Scholar] [CrossRef] [PubMed]
  5. Hayes, J.D.; Pulford, D.J. The glutathione s-transferase supergene family: Regulation of GST and the contribution of the isoenzymes to cancer chemoprotection and drug resistance part I. Crit. Rev. Biochem. Mol. Biol. 1995, 30, 445–520. [Google Scholar] [CrossRef]
  6. Ranson, H.; Collins, F.; Hemingway, J. The role of alternative mRNA splicing in generating heterogeneity within the Anopheles gambiae class I glutathione s-transferase family. Proc. Natl. Acad. Sci. USA 1998, 95, 14284–14289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Che-Mendoza, A.; Penilla, R.P.; Rodriguez, D.A. Insecticide resistance and glutathione S-transferases in mosquitoes: A review. Afr. J. Biotechnol. 2009, 8, 1386–1397. [Google Scholar]
  8. Ranson, H.; Prapanthadara, L.A.; Hemingway, J. Cloning and characterization of two glutathione S-transferases from a DDT-resistant strain of Anopheles gambiae. Biochem. J. 1997, 324, 97–102. [Google Scholar] [CrossRef] [Green Version]
  9. Dou, W.; Wu, S.; Hassan, M.W.; Wang, J.J. Purification and biochemical characterization of glutathione S-transferases from three strains of Liposcelis bostrychophila Badonnel (Psocoptera: Liposcelididae): Implication of insecticide resistance. Pestic. Biochem. Physiol. 2009, 94, 10–14. [Google Scholar] [CrossRef]
  10. Hemingway, J.; Ranson, H. Insecticide resistance in insect vectors of human disease. Annu. Rev. Entomol. 2000, 45, 371–391. [Google Scholar] [CrossRef]
  11. Hemingway, J. The molecular basis of two contrasting metabolic mechanisms of insecticide resistance. Insect Biochem. Mol. Biol. 2000, 30, 1009–1015. [Google Scholar] [CrossRef]
  12. Lewis, A.D.; Hayes, J.D.; Wolf, C.R. Glutathione and glutathione-dependent enzymes in ovarian adenocarcinoma cell lines derived from a patient before and after the onset of drug resistance: Intrinsic differences and cell cycle effects. Carcinogenesis 1988, 9, 1283–1287. [Google Scholar] [CrossRef] [PubMed]
  13. Ranson, H.; Rossiter, L.; Ortelli, F.; Jensen, B.; Wang, X.; Roth, C.W.; Collins, F.H.; Hemingway, J. Identification of a novel class of insect glutathione S-transferases involved in resistance to DDT in the malaria vector Anopheles Gambiae. Biochem. J. 2001, 359, 295–304. [Google Scholar] [CrossRef] [PubMed]
  14. Ranson, H.; Claudianos, C.; Ortelli, F.; Abgrall, C.; Hemingway, J.; Sharakhova, M.V.; Unger, M.F.; Collins, F.H.; Feyereisen, R. Evolution of supergene families associated with insecticide resistance. Science 2002, 298, 179–181. [Google Scholar] [CrossRef] [PubMed]
  15. Ding, Y.; Hawkes, N.; Meredith, J.; Eggleton, P.; Hemingway, J.; Ranson, H. Characterization of the promoters of Epsilon glutathione S-transferases in the mosquito Anopheles gambiae and their response to oxidative stress. Biochem. J. 2005, 387, 879–888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hamzah, S.N.; Alias, Z. Purification, expression and partial characterisation of glutathione S-transferases (GSTs) from three different strains of Aedes albopictus (Diptera: Culicidae). Trop. Biomed. 2016, 33, 335–347. [Google Scholar]
  17. Penilla, R.P.; Rodriguez, A.D.; Hemingway, J.; Torres, J.L.; Solis, F.; Rodriguez, M.H. Changes in glutathione s-transferase activity in DDT resistant natural Mexican populations of Anopheles albimanus under different insecticide resistance management strategies. Pestic. Biochem. Physiol. 2006, 86, 63–71. [Google Scholar] [CrossRef]
  18. Hemingway, J.; Hawkes, N.; Prapanthadara, L.A.; Jayawardenal, K.I.; Ranson, H. The role of gene splicing, gene amplification and regulation in mosquito insecticide resistance. Philos. Trans. R. Soc. B Biol. Sci. 1998, 353, 1695–1699. [Google Scholar] [CrossRef]
  19. Caminade, C.; Medlock, J.M.; Ducheyne, E.; McIntyre, K.M.; Leach, S.; Baylis, M.; Morse, A.P. Suitability of European climate for the Asian tiger mosquito Aedes albopictus: Recent trends and future scenarios. J. R. Soc. Interface 2012, 9, 2708–2717. [Google Scholar] [CrossRef] [Green Version]
  20. Vanlerberghe, V.; Toledo, M.E.; Rodriguez, M.; Gomez, D.; Baly, A.; Benitez, J.R.; Van Der Stuyft, P. Community involvement in dengue vector control: Cluster randomised trial. Br. Med. J. 2009, 338, b1959. [Google Scholar] [CrossRef] [Green Version]
  21. Ab Hamid, N.; Noor, S.N.M.; Saadatian-Elahi, M.; Isa, N.R.; Rodzay, R.M.; Ruslan, B.M.; Omar, T.; Norsham, M.I.M.; Amanzuri, N.H.; Abd Khalil, N.; et al. Residual spray for the control of Aedes vectors in dengue outbreak residential areas. Adv. Entomol. 2019, 7, 105. [Google Scholar] [CrossRef] [Green Version]
  22. Wan-Norafikah, O.; Nazni, W.A.; Lee, H.L.; Zainol-Ariffin, P.; Sofian-Azirun, M. Susceptibility of Aedes albopictus Skuse (Diptera: Culicidae) to permethrin in Kuala Lumpur, Malaysia. Asian Biomed. 2013, 7, 51–62. [Google Scholar]
  23. Koou, S.Y.; Chong, C.S.; Vythilingam, I.; Ng, L.C.; Lee, C.Y. Pyrethroid resistance in Aedes aegypti larvae (Diptera: Culicidae) from Singapore. J. Med. Entomol. 2014, 51, 1–12. [Google Scholar] [CrossRef] [Green Version]
  24. Grigoraki, L.; Lagnel, J.; Kioulos, I.; Kampouraki, A.; Morou, E.; Labbé, P.; Weill, M.; Vontas, J. Transcriptome profiling and genetic study reveal amplified carboxylesterase genes implicated in temephos resistance, in the Asian tiger mosquito Aedes albopictus. PLoS Negl. Trop. Dis. 2015, 9, e0003771. [Google Scholar] [CrossRef] [PubMed]
  25. Sokhna, C.; Ndiath, M.O.; Rogier, C. The changes in mosquito vector behaviour and the emerging resistance to insecticides will challenge the decline of malaria. Clin. Microbiol. Infect. 2013, 19, 902–907. [Google Scholar] [CrossRef] [Green Version]
  26. Li, X.; Schuler, M.A.; Berenbaum, M.R. Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu. Rev. Entomol. 2007, 52, 231–253. [Google Scholar] [CrossRef]
  27. Ranson, H.; N’Guessan, R.; Lines, J.; Moiroux, N.; Nkuni, Z.; Corbel, V. Pyrethroid resistance in African anopheline mosquitoes: What are the implications for malaria control? Trends Parasitol. 2011, 27, 91–98. [Google Scholar] [CrossRef]
  28. World Health Organisation Instructions for Determining the Susceptibility or Resistance of Mosquito Larvae to Insecticides. WHO/VBC/81.807. 1981. Available online: https://apps.who.int/iris/bitstream/handle/10665/69615/WHO_VBC_81.807_eng.pdf (accessed on 8 March 2018).
  29. World Health Organization Guidelines for Laboratory and Field Testing of Mosquito Larvicides. WHO/CDS/WHOPES/GCDPP/2005.13. 2005. Available online: https://apps.who.int/iris/bitstream/handle/10665/69101/WHO_CDS_WHOPES_GCDPP_2005.13.pdf (accessed on 1 August 2018).
  30. IBM SPSS Statistics for Windows, Version 24.0; IBM Corp: Armonk, NY, USA, 2016.
  31. Bradford, M.M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar] [CrossRef]
  32. World Health Organisation Manual for Techniques to Detect Insecticide Resistance Mechanism (Field and Laboratory Manual). WHO/CDC/CPC/MAL/98.6. World Health Organization: Geneva, Switzerland, 1998. Available online: https://apps.who.int/iris/handle/10665/83780 (accessed on 3 May 2018).
  33. Hamzah, S.N.; Alias, Z. Developmental expression and oxidative stress induction of proteome of glutathione S-transferases in Aedes albopictus (Diptera: Culicidae). J. Asia-Pac. Entomol. 2016, 19, 869–875. [Google Scholar] [CrossRef]
  34. Abbott, W.S. A method of computing the effectiveness of an insecticide. J. Econ. Entomol. 1925, 18, 265–267. [Google Scholar] [CrossRef]
  35. Marriel, N.B.; Tomé, H.V.V.; Guedes, R.C.N.; Martins, G.F. Deltamethrin-mediated survival, behaviour, and ooenocyte morphology of insecticide-susceptible and resistant yellow fever mosquitoes (Aedes aegypti). Acta Trop. 2016, 158, 88–96. [Google Scholar] [CrossRef] [PubMed]
  36. Brogdon, W.G.; Janet, C.; McAllister, J.C.; Vulule, J. Heme peroxidase activity measured in single mosquitoes identifies individuals expressing an elevated oxidase for insecticide resistance. J. Am. Mosq. Control Assoc. 1997, 13, 233–237. [Google Scholar] [PubMed]
  37. Poupardin, R.; Reynaud, S.; Strode, C.; Ranson, H.; Vontas, J.; David, J.P. Cross-induction of detoxification genes by environmental xenobiotics and insecticides in the mosquito Aedes aegypti: Impact on larval tolerance to chemical insecticides. Insect Biochem. Mol. Biol. 2008, 38, 540–551. [Google Scholar] [CrossRef] [PubMed]
  38. Devonshire, A.L.; Field, L.M.; Williams, M.S. Molecular Biology of Insecticide Resistance; Academic Press: London, UK, 1992. [Google Scholar]
  39. Yan, H.; Jia, H.; Gao, X.; Xu, B. Identification, genomic organization, and oxidative stress response of a sigma class glutathione s-transferase gene (AccGSTS1) in the honey bee, Apis Cerana Cerana. Cell Stress Chaperones 2013, 18, 415–462. [Google Scholar] [CrossRef]
  40. Wang, X.; Martínez, M.A.; Dai, M.; Chen, D.; Ares, I.; Romero, A.; Castellano, V.; Martínez, M.; Rodríguez, J.L.; Martínez-Larrañaga, M.R.; et al. Permethrin-induced oxidative stress and toxicity and metabolism. A review. Environ. Res. 2016, 149, 86–104. [Google Scholar] [CrossRef]
  41. Ullah, S.; Li, Z.; Hasan, Z.; Khan, S.U.; Fahad, S. Malathion induced oxidative stress leads to histopathological and biochemical toxicity in the liver of rohu (Labeo rohita, Hamilton) at acute concentration. Ecotoxicol. Environ. Saf. 2018, 161, 270–280. [Google Scholar] [CrossRef]
  42. Panini, M.; Manicardi, G.C.; Moores, G.D.; Mazzoni, E. An overview of the main pathways of metabolic resistance in insects. Invertebr. Surviv. J. 2016, 13, 326–335. [Google Scholar]
  43. Montella, I.R.; Martins, A.J.; Viana-Medeiros, P.F.; Lima, J.B.P.; Braga, I.A.; Valle, D. Insecticide resistance mechanisms of Brazilian Aedes aegypti populations from 2001 to 2004. Am. J. Trop. Med. Hyg. 2007, 77, 467–477. [Google Scholar] [CrossRef] [Green Version]
  44. Jagadeshwaran, U.; Vijayan, V.A. Biochemical characterization of deltamethrin resistance in a laboratory-selected strain of Aedes aegypti. Parasitol. Res. 2009, 104, 1431–1438. [Google Scholar] [CrossRef]
  45. Amelia-Yap, Z.H.; Chen, C.D.; Sofian-Azirun, M.; Low, V.L. Pyrethroid resistance in the dengue vector Aedes aegypti in the Southeast Asia: Present situation in the prospects for management. Parasites Vectors 2018, 11, 332. [Google Scholar] [CrossRef] [Green Version]
  46. Jeffrey, S.S.; Lønning, P.E.; Hillner, B.E. Genomics-based prognosis and therapeutic prediction in breast cancer. J. Natl. Compr. Cancer Netw. 2005, 3, 291–300. [Google Scholar] [CrossRef] [PubMed]
  47. Sehrawat, N.; Gakhar, S.K. Mosquito proteomics: Present and future prospective. Res. Biotechnol. 2014, 5, 25–33. [Google Scholar]
  48. Rosilawati, R.; Nabila, R.; Siti Futri Farahininajua, F.; Nazni, W.A.; Lee, H.L. A preliminary proteomic study of permethrin resistant and susceptible Aedes aegypti (L.). Trop. Biomed. 2019, 36, 855–865. [Google Scholar] [PubMed]
  49. Müller, P.; Warr, E.; Stevenson, B.J.; Pignatelli, P.M.; Morgan, J.C.; Steven, A.; Yawson, A.E.; Mitchell, S.N.; Ranson, H.; Hemingway, J.; et al. Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genet. 2008, 4, e1000286. [Google Scholar] [CrossRef] [Green Version]
  50. Ngoagouni, C.; Kamgang, B.; Brengues, C.; Yahouedo, G.; Paupy, C.; Nakouné, E.; Kazanji, M.; Chandre, F. Susceptibility profile and metabolic mechanisms involved in Aedes aegypti and Aedes albopictus resistant to DDT and deltamethrin in the Central African Republic. Parasites Vectors 2016, 9, 599. [Google Scholar] [CrossRef] [Green Version]
  51. Paiva, M.H.; Lovin, D.D.; Mori, A.; Melo-Santos, M.A.; Severson, D.W.; Ayres, C.F. Identification of a major quantitative trait locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes Aegypti. Genomics 2016, 107, 40–48. [Google Scholar] [CrossRef]
  52. Hamzah, S.N.; Farouk, S.; Alias, Z. Isoenzymes of Aedes albopictus (Diptera: Culicidae) Glutathione S-transferases: Isolation and expression after acute insecticide treatment. Pestic. Biochem. Physiol. 2019, 153, 116–121. [Google Scholar] [CrossRef]
  53. Ito, T.; Tashiro, K.; Muta, S.; Ozawa, R.; Chiba, T.; Nishizawa, M.; Sakaki, Y. Toward a protein-protein interaction map of the budding yeast: A comprehensive system to examine two-hybrid interactions in all possible combinations between the yeast proteins. Proc. Natl. Acad. Sci. USA 2000, 97, 1143–1147. [Google Scholar] [CrossRef] [Green Version]
  54. Liu, Q.; Zheng, J.; Sun, W.; Huo, Y.; Zhang, L.; Hao, P.; Zhuang, M. A proximity-tagging system to identify membrane protein-protein interactions. Nat. Methods 2018, 15, 715–722. [Google Scholar] [CrossRef]
  55. Xin, S.; Zhang, W. Construction and analysis of the protein–protein interaction network for the olfactory system of the silkworm Bombyx mori. Arch. Insect Biochem. Physiol. 2020, 105, e21737. [Google Scholar] [CrossRef]
  56. Zhang, G.; Zhang, W. Protein-protein interaction network analysis of insecticide resistance molecular mechanism in Drosophila melanogaster. Arch. Insect Biochem. Physiol. 2019, 100, e21523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Lumjuan, N.; Stevenson, B.J.; Prapanthadara, L.A.; Somboon, P.; Brophy, P.M.; Loftus, B.J.; Severson, D.W.; Ranson, H. The Aedes aegypti glutathione s-transferase family. Insect Biochem. Mol. Biol. 2007, 37, 1026–1035. [Google Scholar] [CrossRef] [PubMed]
  58. Wongsantichon, J.; Robinson, R.C.; Ketterman, A.J. Epsilon glutathione transferases possess a unique class-conserved subunit interface motif that directly interacts with glutathione in the active site. Biosci. Rep. 2015, 35, e00272. [Google Scholar] [CrossRef] [PubMed]
  59. Alias, Z.; Clark, A. Studies on the glutathione s-transferase proteome of adult Drosophila melanogaster: Responsiveness to chemical challenge. Proteomics 2007, 7, 2618–3628. [Google Scholar] [CrossRef] [PubMed]
  60. Nardini, L.; Christian, R.N.; Coetzer, N.; Koekemoer, L.L. DDT and pyrethroid resistance in Anopheles arabiensis from South Africa. Parasites Vectors 2013, 6, 229–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Blackburn, A.C.; Matthaei, K.I.; Lim, C.; Taylor, M.C.; Cappello, J.Y.; Hayes, J.D.; Anders, M.W.; Board, P.G. Deficiency of glutathione transferase zeta causes oxidative stress and activation of antioxidant response pathways. Mol. Pharmacol. 2005, 69, 650–657. [Google Scholar] [CrossRef]
  62. Board, P.G.; Coggan, M.; Chelvanayagam, G.; Easteal, S.; Jermin, L.S.; Schulte, G.K.; Danley, D.E.; Hoth, L.R.; Griffor, M.C.; Kamath, A.V.; et al. Identification, characterization, and crystal structure of the Omega class glutathione s-transferases. J. Biol. Chem. 2000, 275, 24798–24806. [Google Scholar] [CrossRef] [Green Version]
  63. Wildenburg, G.; Liebau, E.; Henkle-Dührsen, K. Onchocerca volvulus: Ultrastructural localization of two glutathione s-transferases. Exp. Parasitol. 1998, 88, 34–42. [Google Scholar] [CrossRef]
Figure 1. The interactive network view of Aedes GSTs with predicted proteins of protein-protein interactions analysis using STRING 11.5 tool at PPI enrichment p-value < 1.0 × 10−16. Network nodes are proteins, and the edges represent the predicted functional associations. The K-means algorithm was used to cluster the proteins in 13 different GST classes. Every colored node corresponds to a cluster. Inter-cluster edges are represented by dashed lines.
Figure 1. The interactive network view of Aedes GSTs with predicted proteins of protein-protein interactions analysis using STRING 11.5 tool at PPI enrichment p-value < 1.0 × 10−16. Network nodes are proteins, and the edges represent the predicted functional associations. The K-means algorithm was used to cluster the proteins in 13 different GST classes. Every colored node corresponds to a cluster. Inter-cluster edges are represented by dashed lines.
Insects 13 01028 g001
Figure 2. The interactive network view among Aedes GSTs of predicted protein-protein interaction analysis using STRING 11.5 tool. No. of nodes: 13; no of edge: 12; average node degree: 11.7; average local clustering coefficient: 0.977; PPI enrichment p-value < 1.0 × 10−16.
Figure 2. The interactive network view among Aedes GSTs of predicted protein-protein interaction analysis using STRING 11.5 tool. No. of nodes: 13; no of edge: 12; average node degree: 11.7; average local clustering coefficient: 0.977; PPI enrichment p-value < 1.0 × 10−16.
Insects 13 01028 g002
Table 1. In vivo glutathione S-transferase isoforms identified in susceptible, permethrin-treated, and malathion-treated Aedes albopictus.
Table 1. In vivo glutathione S-transferase isoforms identified in susceptible, permethrin-treated, and malathion-treated Aedes albopictus.
AccessionIdentified ProteinClassMW (kDa)pI ValuePeptide
No.
AA CoverageProposed Isoform Identification
J9HXZ8AAEL001061-PC GN = GSTD1Delta123.8455.968211GSTD1-1
J9HHL7AAEL001061-PBDelta123.7786.077209GSTD1-2
Q17MA9AAEL001061-PA GN = GSTD1Delta124.7555.442219GSTD1-3
A0A1S4EXN6Glutathione S-transferase GN = 5568355Delta124.5385.92216GSTD1-4
Q6PTY1Glutathione S-transferase OX = 7159Delta423.1467.942205GSTD4-1
A0A0N8ES64Glutathione S-transferase GSTD4 OX = 7159Delta424.8556.141218GSTD4-2
Q17MB8AAEL001054-PA GN = GSTD4Delta424.0856.141211GSTD4-3
Q16SH6AAEL010591-PA GN = GSTD6Delta628.2155.741249GSTD6-1
A0A023ENG1Glutathione S-transferase e2 OX = 7160Delta628.1485.541249GSTD6-2
Q16SH7AAEL010582-PA GN = GSTD11Delta1125.8795.953222GSTD11-1
Q170C6AAEL007947-PA GN = GSTE3Epsilon324.8246.541222GSTE3-1
Q5PY78AAEL007962-PA GN = GSTe4Epsilon425.037.113224GSTE4-1
A0A0P6IV26Putative glutathione S-transferase e4Epsilon427.2646.812244GSTE4-2
Q170C9AAEL007964-PA GN = GSTE5Epsilon524.6355.21221GSTE5-1
Q170C7AAEL007946-PA GN = GSTE6Epsilon624.7385.972220GSTE6-1
Q16P53AAEL011752-PA GN = GSTI1Iota126.1776.151231GSTI1-1
J9E9C0AAEL017085-PA GN = GSTO1Omega128.5897.036248GSTO1-1
A0A1S4G560Glutathione S-transferase GN = 23687505Omega129.6188.516257GSTO1-2
A0A0P6J0T5Glutathione S-transferase OX = 7159Omega129.5728.516257GSTO1-3
Q1HQK1Glutathione S-transferaseOX = 7159Omega128.6087.426248GSTO1-4
Q16P79AAEL011741-PA GN = GSTS1Sigma123.2495.29203GSTS1-1
Q16P80AAEL011741-PB GN = GSTS1Sigma123.2165.172203GSTS1-2
Q0C791AAEL000092-PA GN = GSTX1GSTX124.8175.784218GSTX1-1
Q16NL9AAEL011934-PA GN = GSTZ1Zeta126.396.921233GSTZ1-1
MW: molecular weight; kDa: kilodalton, pI: isoelectric point value; AA: amino acid. All identified peptides were validated based on q-value (Percolator® algorithm) with less than 5% false discovery rate.
Table 2. Characteristics of input proteins and functional partners in GST-class protein-protein interactions from STRING 11.5 analysis.
Table 2. Characteristics of input proteins and functional partners in GST-class protein-protein interactions from STRING 11.5 analysis.
Input ProteinInformation
GSTD1Glutathione s-transferase 1 isoform ×1; Glutathione S-transferase (GSTD1); Glutathione transferase (219 aa)
GSTE4Glutathione S-transferase e4; Glutathione transferase (224 aa)
GSTI1Glutathione S-transferase 1; Belongs to the GST superfamily (231 aa)
GSTD4Glutathione S-transferase 1; AAEL001054-PA; Glutathione transferase (211 aa)
GSTD6Glutathione S-transferase 1; Glutathione transferase; Belongs to the GST superfamily (249 aa)
GSTE3Glutathione S-transferase 1; Belongs to the GST superfamily (222 aa)
GSTX1Glutathione S-transferase d4; Glutathione transferase (218 aa)
GSTS1Prostaglandin-H2 D-isomerase/glutathione transferase; Glutathione transferase (203 aa)
GSTE5Glutathione S-transferase 1; Belongs to the GST superfamily (221 aa)
GSTZ1Probable maleylacetoacetate isomerase 2 isoform x1; Belongs to the GST superfamily (233 aa)
GSTD11Glutathione S-transferase 1-1; Belongs to the GST superfamily (222 aa)
GSTE6Glutathione S-transferase 1; Glutathione transferase (220 aa)
GSTO1Pyrimidodiazepine synthase; Belongs to the GST superfamily (248 aa)
Predicted Functional Partner
AccessionInformationScore
AAEL013637Homogentisate 1,2-dioxygenase0.999
AAEL011973Fumarylacetoacetate hydrolase0.998
AAEL009154Glutathione synthetase isoform x1; Glutathione synthetase0.968
AAEL004088Aldose reductase isoform x1; Aldo-keto reductase0.958
AAEL004102Aldehyde reductase; Aldo-keto reductase0.958
AAEL013637Homogentisate 1,2-dioxygenase0.999
AAEL011973Fumarylacetoacetate hydrolase0.998
AAEL009154Glutathione synthetase isoform x1; Glutathione synthetase0.968
AAEL004088Aldose reductase isoform x1; Aldo-keto reductase0.958
AAEL004102Aldehyde reductase; Aldo-keto reductase0.958
Table 3. Characteristics of biological process and molecular function of protein-protein interaction of GSTs with metabolic detoxification related protein STRING 11.5.
Table 3. Characteristics of biological process and molecular function of protein-protein interaction of GSTs with metabolic detoxification related protein STRING 11.5.
Biological Process (GO)
GO-TermPathway DescriptionCount in Networkp-Values
GO:0000413Protein peptidyl-prolyl isomerization3 of 290.0134
GO:0006457Protein folding8 of 1342.39 × 10−6
GO:0006518Peptide metabolic process21 of 4711.25 × 10−17
GO:0006559L-phenylalanine catabolic process3 of 50.00028
GO:0006570Tyrosine metabolic process4 of 134.58 × 10−5
GO:0006572Tyrosine catabolic process3 of 30.00014
GO:0006749Glutathione metabolic process20 of 431.91 × 10−34
GO:0006750Glutathione biosynthetic process3 of 40.00020
GO:0006751Glutathione catabolic process4 of 123.76 × 10−5
GO:0006807Nitrogen compound metabolic process29 of 43032.34 × 10−5
GO:0008152Metabolic process34 of 58331.22 × 10−5
GO:0009074Aromatic amino acid family catabolic process5 of 125.50 × 10−7
GO:0009987Cellular process39 of 87010.00012
GO:0034605Cellular response to heat4 of 390.00016
GO:0034641Cellular nitrogen compound metabolic process23 of 20980.0011
GO:0034605Cellular response to heat4 of 394.25 × 10−7
GO:0043043Peptide biosynthetic process7 of 3660.0093
GO:0043171Peptide catabolic process5 of 540.00017
GO:0044237Cellular metabolic process31 of 44552.26 × 10−6
GO:0044248Cellular catabolic process10 of 8850.0164
GO:0050821Protein stabilization4 of 470.0020
GO:0071704Organic substance metabolic process29 of 51310.00063
GO:0071722Detoxification of arsenic-containing substance2 of 20.0072
GO:0098754Detoxification4 of 710.0080
GO:1901564Organonitrogen compound metabolic process28 of 31832.23 × 10−7
GO:1901565Organonitrogen compound catabolic process10 of 5970.00081
GO:1901606Alpha-amino acid catabolic process5 of 600.00024
Molecular Function (GO)
GO:0000048Peptidyltransferase activity4 of 101.57 × 10−5
GO:0003755Peptidyl-prolyl cis-trans isomerase activity3 of 300.0089
GO:0003824Catalytic activity33 of 47212.22 × 10−7
GO:00038684-hydroxyphenylpyruvate dioxygenase activity2 of 20.0049
GO:0004032alditol:NADP+ 1-oxidoreductase activity2 of 110.0373
GO:0004364Glutathione transferase activity13 of 304.08 × 10−21
GO:0008144Drug binding3 of 160.0027
GO:0016018Cyclosporin a binding2 of 90.0283
GO:0016702Oxidoreductase activity, acting on single donors with incorporation of molecular oxygen, incorporation of two atoms of oxygen3 of 160.0027
GO:0016740Transferase activity17 of 13851.57 × 10−5
GO:0016823Hydrolase activity, acting on acid carbon-carbon bonds, in ketonic substances2 of 20.0049
GO:0016848Carbon-halide lyase activity 4.96 × 10−6
GO:0016853Isomerase activity6 of 1080.00012
GO:0016859Cis-trans isomerase activity4 of 320.00044
GO:0031072Heat shock protein binding3 of 420.0199
GO:0033218Amide binding5 of 1380.0049
GO:0036374Glutathione hydrolase activity4 of 101.57 × 10−5
GO:0042277Peptide binding4 of 1050.0180
GO:0050220Prostaglandin-E synthase activity2 of 30.0069
GO:0051879Hsp90 protein binding2 of 100.0327
Count In Network: The first number indicates how many proteins in the network are annotated with a particular term. The second number indicates how many proteins in total (in the network and in the background) have this term assigned; p-value: the false discovery rate which corrected for multiple testing within each category using the Benjamini–Hochberg procedure to measure how significant the enrichment is.
Table 4. Characteristics of KEGG pathways of protein-protein interaction of GSTs with STRING 11.5.
Table 4. Characteristics of KEGG pathways of protein-protein interaction of GSTs with STRING 11.5.
KEGG Pathway
Pathway IDPathway DescriptionCount in Networkp-Values
aag00051Fructose and mannose metabolism2 of 230.0260
aag00052Galactose metabolism2 of 280.0343
aag00270Cysteine and methionine metabolism3 of 270.0011
aag00350Tyrosine metabolism5 of 273.81 × 10−7
aag00430Taurine and hypotaurine metabolism4 of 73.02 × 10−7
aag00480Glutathione metabolism21 of 742.66 × 10−34
aag00590Arachidonic acid metabolism8 of 174.27 × 10−14
aag00790Folate biosynthesis3 of 330.0017
aag00980Metabolism of xenobiotics by cytochrome P45011 of 533.76 × 10−16
aag00982Drug metabolism—cytochrome P45011 of 523.76 × 10−16
aag00983Drug metabolism—other enzymes10 of 807.93 × 10−13
aag01100Metabolic pathways30 of 10395.51 × 10−24
Count In Network: The first number indicates how many proteins in the network are annotated with a particular term. The second number indicates how many proteins in total (in the network and in the background) have this term assigned; p-value: the false discovery rate which corrected for multiple testing within each category using the Benjamini–Hochberg procedure to measure how significant the enrichment is.
Table 5. Aedes albopictus mortality under malathion and permethrin treatments along with the calculated sub-lethal dosage (LC50).
Table 5. Aedes albopictus mortality under malathion and permethrin treatments along with the calculated sub-lethal dosage (LC50).
InsecticideConcentration
(mg/L)
Mortality
(%)
LC50
(mg/L)
95% Confidence Limit (mg/L)Slope ± SD
Lower BoundUpper Bound
Malathion0.0240
0.05315
0.07230
0.115600.0990.0900.1083.692 ± 0.228
0.2594
0.40100
Permethrin0.0020
0.00610
0.018400.0230.0190.0273.482 + 0.223
0.03465
0.12095
0.203100
The LC50 values were calculated statistically from the data collected from the WHO bioassays using the Probit analysis software SPSS version 24.
Table 6. Mean of total protein contents (mean ± SD), total GST enzyme activities (mean ± SD), and specific GST enzyme activities (mean ± SD) in Aedes albopictus larvae upon 24 h of acute treatment by permethrin and malathion sub-lethal dosages (LC50).
Table 6. Mean of total protein contents (mean ± SD), total GST enzyme activities (mean ± SD), and specific GST enzyme activities (mean ± SD) in Aedes albopictus larvae upon 24 h of acute treatment by permethrin and malathion sub-lethal dosages (LC50).
TreatmentTreatment DurationTotal Protein Contents (mg)PFTotal Activity
(μmol/min)
TAFSpecific Activity
(μmol/min/mg)
SAF
PermethrinControl
24 h
3.9 × 10−3 ± 0.001
5.3 × 10−3 ± 0.006 *
1.00
1.36
7.59 × 10−3 ± 0.01
13.97 × 10−3 ± 0.01 *
1.00
1.84
1.947 ± 0.08
2.635 ± 0.08 *
1.00
1.35
MalathionControl
24 h
3.9 × 10−3 ± 0.001
4.3 × 10−3 ± 0.002 *
1.00
1.10
7.59 × 10−3 ± 0.04
10.75 × 10−3 ± 0.02*
1.00
1.41
1.947 ± 0.08
2.502 ± 0.05 *
1.00
1.29
Mean of total protein content, total activity, and specific activity with an asterisk (*) shows a significant difference at p = 0.05. The t-test was performed to compare the total protein contents as well as the enzyme activities between each acute insecticide-treated strain with susceptible strain only. SD: standard deviation, calculations were made per individual. PF: total protein content fold-change. TAF: total GST enzymatic activity fold-change. SAF: specific GST enzymatic activity fold-change. Fold-change calculations were calculated by dividing the values of each strain by its control counterpart. n = 100.
Table 7. Peptide abundance response of glutathione s-transferase isoforms in Aedes albopictus larvae upon 24 h of acute permethrin and malathion treatments.
Table 7. Peptide abundance response of glutathione s-transferase isoforms in Aedes albopictus larvae upon 24 h of acute permethrin and malathion treatments.
Peptide AbundanceFold-Change
NoProposed IdentificationControlPermethrinMalathionPermethrinMalathion
1GSTD1-1155.8180.2155.31.161.00
2GSTD1-248.5107.3102.42.212.11
3GSTD1-327.662.553.12.261.92
4GSTD1-427.662.553.12.261.92
5GSTD4-1162.8176.5188.81.081.16
6GSTD4-278.595.890.11.221.15
7GSTD4-378.595.890.11.221.15
8GSTD6-1173.3174252.61.001.46
9GSTD6-2173.3174252.61.001.46
10GSTD11-1147.3126131.40.860.89
11GSTE3-1148192227.21.301.54
12GSTE4-1198.3217.7183.91.100.93
13GSTE4-2143.1230.1226.71.611.58
14GSTE5-1175.6116.1191.80.661.09
15GSTE6-1175.6116.1191.80.661.09
16GSTI1-1198.9224.4176.71.130.89
17GSTO1-147.64753.81.001.13
18GSTO1-247.64753.81.001.13
19GSTO1-347.64753.81.001.13
20GSTO1-447.64753.81.001.13
21GSTS1-1165.5149.41390.900.84
22GSTS1-2166.4136.9124.60.820.75
23GSTX1-187.8109.4110.91.251.26
24GSTZ1-116.420.119.41.231.18
GST: Glutathione S-transferase, D: Delta, E: Epsilon, I: Iota, O: Omega, S: Sigma, Z: Zeta. Fold-change of peptide abundance in acute permethrin and malathion treatment was obtained by dividing its peptide abundance value by the peptide abundance value of the control (susceptible).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hamzah, S.N.; Avicor, S.W.; Alias, Z.; Razak, S.A.; Bakhori, S.K.M.; Hsieh, T.C.; Syanizam, N.N.; Farouk, S.A. In Vivo Glutathione S-Transferases Superfamily Proteome Analysis: An Insight into Aedes albopictus Mosquitoes upon Acute Xenobiotic Challenges. Insects 2022, 13, 1028. https://doi.org/10.3390/insects13111028

AMA Style

Hamzah SN, Avicor SW, Alias Z, Razak SA, Bakhori SKM, Hsieh TC, Syanizam NN, Farouk SA. In Vivo Glutathione S-Transferases Superfamily Proteome Analysis: An Insight into Aedes albopictus Mosquitoes upon Acute Xenobiotic Challenges. Insects. 2022; 13(11):1028. https://doi.org/10.3390/insects13111028

Chicago/Turabian Style

Hamzah, Siti Nasuha, Silas Wintuma Avicor, Zazali Alias, Sarah Abdul Razak, Siti Khadijah Mohd Bakhori, Ting Chuan Hsieh, Nurin Nazifa Syanizam, and Salinah Abdul Farouk. 2022. "In Vivo Glutathione S-Transferases Superfamily Proteome Analysis: An Insight into Aedes albopictus Mosquitoes upon Acute Xenobiotic Challenges" Insects 13, no. 11: 1028. https://doi.org/10.3390/insects13111028

APA Style

Hamzah, S. N., Avicor, S. W., Alias, Z., Razak, S. A., Bakhori, S. K. M., Hsieh, T. C., Syanizam, N. N., & Farouk, S. A. (2022). In Vivo Glutathione S-Transferases Superfamily Proteome Analysis: An Insight into Aedes albopictus Mosquitoes upon Acute Xenobiotic Challenges. Insects, 13(11), 1028. https://doi.org/10.3390/insects13111028

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop