Comparative genomics of the anopheline glutathione S-transferase epsilon cluster

Recombinant expression and characterization of GSTd3 from a resistant population of Anopheles arabiensis and comparison of DDTase activity with GSTe2

The development of insecticide resistance in malaria vectors is a challenge for the global effort to control and eradicate malaria. Glutathione S-transferases (GSTs) are multifunctional enzymes involved in the detoxification of many classes of insecticides. For mosquitoes, it is known that overexpression of an epsilon GST, GSTe2, confers resistance towards DDT and pyrethroids. In addition to GSTe2, consistent overexpression of a delta class GST, GSTd3, has been observed in insecticide resistant populations of different malaria vector species. However, the functional role of GSTd3 towards DDT resistance has not yet been investigated. Here, we recombinantly expressed both GSTe2 and GSTd3 from Anopheles arabiensis and compared their metabolic activities against DDT. Both AaGSTd3 and AaGSTe2 exhibited CDNB-conjugating and glutathione peroxidase activity and DDT metabolism was observed for both GSTs. However, the DDT dehydrochlorinase activity exhibited by AaGSTe2 was much higher than for AaGSTd3, and AaGSTe2 was also able to eliminate DDE although the metabolite could not be identified. Molecular modeling revealed subtle differences in the binding pocket of both enzymes and a better fit of DDT within the H-site of AaGSTe2. The overexpression but much lower DDT metabolic activity of AaGSTd3, might suggest that AaGSTd3 sequesters DDT. These findings highlight the complexity of insecticide resistance in the major malaria vectors and the difficulties associated with control of the vectors using DDT, which is still used for indoor residual spraying.

Tandem duplication of a genomic region encoding glutathione S-transferase epsilon-2 and -4 genes in DDT-resistant Anopheles stephensi strain from India

The glutathione S-transferases (GST) genes are a multigene family of enzymes involved in the metabolism of endogenous and xenobiotic compounds by catalysing the conjugation of the reduced form of glutathione to the substrate. The epsilon class of GST (GSTe), unique to arthropods, is known to be involved in the detoxification process of several classes of insecticides, and GSTe2 in particular is known to have DDT dehydrochlorinase activity. This communication reports a tandem duplication of a genomic region encoding GSTe2 and GSTe4 genes in a laboratory-colonized DDT-resistant Anopheles stephensi. We identified duplication breakpoints and the organization of gene duplication through Sanger sequencing performed on long-PCR products. Manual annotation of sequences revealed a tandemly-arrayed duplication of a 3.62 kb segment of GST epsilon gene clusters comprised of five genes: a partial GSTe1, GSTe2, GSTe2-pseudogene, GSTe4 and partial GSTe5, interconnected by a conserved 2.42 kb DNA insert segment major part of which is homologous to a genomic region located on a different chromosome. The tandemly duplicated array contained a total of two GSTe2 and three GSTe4 functional paralog genes. Read-depth coverage and split-read analysis of Illumina-based whole-genome sequence reads confirmed the presence of duplication in the corresponding region of the genome. The increased gene dose in mosquitoes as a result of the GSTe gene-duplication may be an adaptive process to increase levels of detoxifying enzymes to counter insecticide pressure.

Updated occurrence and bionomics of potential malaria vectors in Europe: a systematic review (2000-2021)

Despite the eradication of malaria across most European countries in the 1960s and 1970s, the anopheline vectors are still present. Most of the malaria cases that have been reported in Europe up to the present time have been infections acquired in endemic areas by travelers. However, the possibility of acquiring malaria by locally infected mosquitoes has been poorly investigated in Europe, despite autochthonous malaria cases having been occasionally reported in several European countries. Here we present an update on the occurrence of potential malaria vector species in Europe. Adopting a systematic review approach, we selected 288 papers published between 2000 and 2021 for inclusion in the review based on retrieval of accurate information on the following Anopheles species: An. atroparvus, An. hyrcanus sensu lato (s.l.), An. labranchiae, An. maculipennis sensu stricto (s.s.), An. messeae/daciae, An. sacharovi, An. superpictus and An. plumbeus. The distribution of these potential vector species across Europe is critically reviewed in relation to areas of major presence and principal bionomic features, including vector competence to Plasmodium. Additional information, such as geographical details, sampling approaches and species identification methods, are also reported. We compare the information on each species extracted from the most recent studies to comparable information reported from studies published in the early 2000s, with particular reference to the role of each species in malaria transmission before eradication. The picture that emerges from this review is that potential vector species are still widespread in Europe, with the largest diversity in the Mediterranean area, Italy in particular. Despite information on their vectorial capacity being fragmentary, the information retrieved suggests a re-definition of the relative importance of potential vector species, indicating An. hyrcanus s.l., An. labranchiae, An. plumbeus and An. sacharovi as potential vectors of higher importance, while An. messeae/daciae and An. maculipennis s.s. can be considered to be moderately important species. In contrast, An. atroparvus and An. superpictus should be considered as vectors of lower importance, particularly in relation to their low anthropophily. The presence of gaps in current knowledge of vectorial systems in Europe becomes evident in this review, not only in terms of vector competence but also in the definition of sampling approaches, highlighting the need for further research to adopt the appropriate surveillance system for each species.

Effects of elevated CO2 treatment of Populus davidiana × P. bolleana on growth and detoxifying enzymes in gypsy moth, Lymantria dispar

To date, elevated CO₂ concentrations in the environment caused by various human activities influence diverse areas of life, including the interactions between insects and plants. The Lymantria dispar is one of the most severely destructive pests, which further could inflict ecological and economical damage. In this experiment, one-year-old Populus davidiana × P. bolleana plants were grown in CO₂-enhanced environments for one month at three different CO₂ concentrations: 397 ppm (atmospheric CO₂ concentration), 550 ppm and 750 ppm (two predicted elevated CO₂ concentrations). The 3rd instar L. dispar larvae then fed on the treated poplar seedlings covered in a nylon bag. The L. dispar larvae fed on poplar seedling treated for 96 h showed the highest growth rate at all CO₂ concentrations. Enzymatic activity of treated larvae showed the highest GST and P450 activity at 750 ppm CO₂. The relative expressions of seven CYP and ten GST genes in L. dispar larvae were analyzed quantitatively using real-time RT-PCR, which the results were expressed variably. Compared to 397 ppm CO₂, the expression of CYP4L23 was down-regulated, while the expressions of other CYP genes were up-regulated. Meanwhile, only GSTo1 gene showed down-regulated at 48 h and 96 h in 750 ppm CO₂ treatment, while GST expression level for the other nine GST genes showed up-regulated at 48 h and 72 h. These results offer the insight into plant-insect interactions under global climate change and furthermore will provide essential information for strategic pest control based on biochemical and molecular levels changes in gypsy moths.

Insecticide resistance profiles in malaria vector populations from Sud-Kivu in the Democratic Republic of the Congo

BACKGROUND

Insecticide resistance has become a widespread problem causing a decline in the effectiveness of vector control tools in sub-Saharan Africa. In this situation, ongoing monitoring of vector susceptibility to insecticides is encouraged by the WHO to guide national malaria control programmes. Our study was conducted from April to November 2018 in Tchonka (Sud-Kivu, Democratic Republic of the Congo) and reported primary data on the resistance status of Anopheles funestus and Anopheles gambiae.

METHODS

Insecticide susceptibility bioassays were performed on wild populations of A. funestus and A. gambiae using WHO insecticide-impregnated papers at discriminating concentration. In addition, PCR was performed to identify mosquito species and to detect kdr and ace-1R mutations involved in insecticide resistance.

RESULTS

Bioassay results show resistance to all tested insecticides except pirimiphos-methyl, propoxur, fenitrothion and malathion with a mortality rate ranging from 95.48 to 99.86%. The addition of piperonyl butoxide (PBO) increased the susceptibility of vectors to deltamethrin and alpha-cypermethrin by exhibiting a mortality ranging from 91.50 to 95.86%. The kdr mutation was detected at high frequencies (approximately 0.98) within A. gambiae while ace-1R was not detected.

CONCLUSIONS

This study provides useful data on the insecticide resistance profiles of malaria vector populations to better manage vector control. Our results highlight that, despite the high level of resistance, organophosphorus compounds and pyrethroids + PBO remain effective against the vectors.

Polymorphisms in GSTE2 is associated with temephos resistance in Aedes aegypti

The glutathione S-transferases (GSTs) are enzymes involved in several distinct biological processes. In insects, the GSTs, especially delta and epsilon classes, play a key role in the metabolism of xenobiotics used to control insect populations. Here, we investigated its potential role in temephos resistance, examining the GSTE2 gene from susceptible (RecL) and resistant (RecR) strains of the mosquito Aedes aegypti, vector for several pathogenic arboviruses. Total GST enzymatic activity and the GSTE2 gene expression profile were evaluated, with the GSTE2 cDNA and genomic loci sequenced from both strains. Recombinant GSTE2 and mutants were produced in a heterologous expression system and assayed for enzyme kinetic parameters. These proteins also had their 3D structure predicted through molecular modeling. Our results showed that RecR has a profile of total GST enzymatic activity higher than RecL, with the expression of the GSTE2 gene in resistant larvae increasing six folds. Four exclusive RecR mutations were observed (L111S, I150V, E178A and A198E), which were absent in the laboratory susceptible strains. The enzymatic activity of the recombinant GSTE2 showed different kinetic parameters, with the GSTE2 RecR showing an enhanced ability to metabolize its substrate. The I150V mutation was shown to induce significant changes in catalytic parameters and a 3D modeling of GSTE2 mapped two of the RecR changes (L111S and I150V) near the enzyme's catalytic pocket, also implying an impact on its catalytic activity. Our results reinforce a potential role for GSTE2 in the metabolic resistance phenotype while contributing to the understanding of the molecular basis for the resistance mechanism.

An integrated approach to unravel a crucial structural property required for the function of the insect steroidogenic Halloween protein Noppera-bo

Ecdysteroids are the principal steroid hormones essential for insect development and physiology. In the last 18 years, several enzymes responsible for ecdysteroid biosynthesis encoded by Halloween genes were identified and genetically and biochemically characterized. However, the tertiary structures of these proteins have not yet been characterized. Here, we report the results of an integrated series of in silico, in vitro, and in vivo analyses of the Halloween GST protein Noppera-bo (Nobo). We determined crystal structures of Drosophila melanogaster Nobo (DmNobo) complexed with GSH and 17β-estradiol, a DmNobo inhibitor. 17β-Estradiol almost fully occupied the putative ligand-binding pocket and a prominent hydrogen bond formed between 17β-estradiol and Asp-113 of DmNobo. We found that Asp-113 is essential for 17β-estradiol–mediated inhibition of DmNobo enzymatic activity, as 17β-estradiol did not inhibit and physically interacted less with the D113A DmNobo variant. Asp-113 is highly conserved among Nobo proteins, but not among other GSTs, implying that this residue is important for endogenous Nobo function. Indeed, a homozygous nobo allele with the D113A substitution exhibited embryonic lethality and an undifferentiated cuticle structure, a phenocopy of complete loss-of-function nobo homozygotes. These results suggest that the nobo family of GST proteins has acquired a unique amino acid residue that appears to be essential for binding an endogenous sterol substrate to regulate ecdysteroid biosynthesis. To the best of our knowledge, ours is the first study describing the structural characteristics of insect steroidogenic Halloween proteins. Our findings provide insights relevant for applied entomology to develop insecticides that specifically inhibit ecdysteroid biosynthesis.

Identification of novel Glutathione S-Transferases epsilon 2 mutation in Anopheles maculipennis s.s. (Diptera: Culicidae)

Anopheles maculipennis complex comprises some important malaria vectors in Iran, Middle East, and Europ. The principal way to control of malaria remains on the use of chemical insecticides against its vectors because there is no vaccine for malaria prevention. Extensive use of organophosphate compounds has caused to emergence and distribution of insecticide resistance in Anopheles species in Asia. The current study aimed to the detection of three well-known amino acid substitutions (I114T, L119F, and F120L) in the Glutathione S-Transferases epsilon 2 (GSTe2) gene are associated with DDT and organophosphate insecticides resistance in an Anopheles maculipennis population collected from Iran. Adult samples of An. maculipennis were collected by hand and Total catch in Animal and Human Shelters from Azerbaijan-Gharbi and Zanjan provinces. Following morphological identification, DNA was extracted by YTA Genomic DNA Extraction Mini Kit for amplification of rDNA-ITS2 and GSTe2 fragments. ∼500 bp fragment was amplified using F rDNA-ITS2 and GSTe2 primers. rDNA-ITS2 sequence analysis showed 100% similarity with An. maculipennis. GSTe2 nucleotide sequence similarity within species was 99-100%, while, it was 95-96 % when compared with Anopheles sacharovi GSTe2 sequences available in GenBank. Amino acid sequence comparisons showed a novel amino acid substitution in N148D position with 15.79% frequency. The current study reports new GSTe2 amino acid substitution in An. maculipennis s.s., for the first time. The function of the mutation N148D and its association with resistance phenotype need to validate. However, the integration of these data into the malaria control program still remains a challenge.

Detection of alleles associated with resistance to chemical insecticide in the malaria vector Anopheles arabiensis in Santiago, Cabo Verde

BACKGROUND

Mosquitoes of the Anopheles gambiae complex are the main malaria vectors worldwide. Due to the lack of a vaccine to prevent malaria, the principal way to reduce the impact of this disease relies on the use of chemical insecticides to control its vectors. However, the intensive use of such compounds has led to the emergence of insecticide resistance in several Anopheles populations in Africa. This study aimed to investigate the presence of resistance alleles in an Anopheles arabiensis population from the City of Praia, capital of the Archipelago Cabo Verde, one of the countries on the World Health Organization list of countries that are on a path to eliminate local transmission of malaria.

METHODS

Larvae from the Anopheles genus were collected using a one-pint dipper in three areas of City of Praia. Larvae were fed and maintained until the emergence of adult mosquitoes, and these were morphologically identified. In addition, molecular identification was performed using IGS markers and all An. arabiensis samples were subjected to PCR to screen for mutations associated to resistance in the Ace-1, Na_v and GSTE2 genes.

RESULTS

From a total of 440 mosquitoes collected, 52.3% were morphologically identified as An. gambiae sensu lato (s.l.) and 46.7% as Anopheles pretoriensis. The molecular identification showed that 100% of the An. gambiae s.l. were An. arabiensis. The mutations G119S in the Ace-1 gene and L119F in the GSTE2 gene were screened but not found in any sample. However, sequencing analysis for GSTE2 revealed the presence of 37 haplotypes, 16 polymorphic sites and a high genetic diversity (π = 2.67). The L1014S mutation in the Na_v (voltage-gated sodium channel gene) was detected at a frequency of 7.3%.

CONCLUSION

This is the first study to investigate the circulation of insecticide resistance alleles in An. arabiensis from Cabo Verde. The circulation of the L1014S allele in the population of An. arabiensis in the city of Praia suggests that pyrethroid resistance may arise, be quickly selected, and may affect the process of malaria elimination in Cabo Verde. Molecular monitoring of resistance should continue in order to guide the development of strategies to be used in vector control in the study region.

Identification of a major Quantitative Trait Locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes aegypti

Organophosphate insecticides (OP) have extensively been used to control mosquitoes, such as the vector Aedes aegypti. Unfortunately, OP resistance has hampered control programs worldwide. We used Quantitative Trait Locus (QTL) mapping to evaluate temephos resistance in two F1 intercross populations derived from crosses between a resistant Ae. aegypti strain (RecR) and two susceptible strains (MoyoD and Red). A single major effect QTL was identified on chromosome 2 of both segregating populations, named rtt1 (resistance to temephos 1). Bioinformatics analyses identified a cluster of carboxylesterase genes (CCE) within the rtt1 interval. qRT-PCR demonstrated that different CCEs were up-regulated in F2 resistant individuals from both crosses. However, none exceeded the 2-fold expression. Primary mechanisms for temephos resistance may vary between Ae. aegypti populations, yet also appear to support previous findings suggesting that multiple linked esterase genes may contribute to temephos resistance in the RecR strain as well as other populations.

The silkworm glutathione S-transferase gene noppera-bo is required for ecdysteroid biosynthesis and larval development

Insect molting and metamorphosis are tightly controlled by ecdysteroids, which are important steroid hormones that are synthesized from dietary sterols in the prothoracic gland. One of the ecdysteroidogenic genes in the fruit fly Drosophila melanogaster is noppera-bo (nobo), also known as GSTe14, which encodes a member of the epsilon class of glutathione S-transferases. In D. melanogaster, nobo plays a crucial role in utilizing cholesterol via regulating its transport and/or metabolism in the prothoracic gland. However, it is still not known whether the orthologs of nobo from other insects are also involved in ecdysteroid biosynthesis via cholesterol transport and/or metabolism in the prothoracic gland. Here we report genetic evidence showing that the silkworm Bombyx mori ortholog of nobo (nobo-Bm; GSTe7) is essential for silkworm development. nobo-Bm is predominantly expressed in the prothoracic gland. To assess the functional importance of nobo-Bm, we generated a B. mori genetic mutant of nobo-Bm using TALEN-mediated genome editing. We show that loss of nobo-Bm function causes larval arrest and a glossy cuticle phenotype, which are rescued by the application of 20-hydroxyecdysone. Moreover, the prothoracic gland cells isolated from the nobo-Bm mutant exhibit an abnormal accumulation of 7-dehydrocholesterol, a cholesterol metabolite. These results suggest that the nobo family of glutathione S-transferases is essential for development and for the regulation of sterol utilization in the prothoracic gland in not only the Diptera but also the Lepidoptera. On the other hand, loss of nobo function mutants of D. melanogaster and B. mori abnormally accumulates different sterols, implying that the sterol utilization in the PG is somewhat different between these two insect species.

Glutathione S-Transferase (GST) Gene Diversity in the Crustacean Calanus finmarchicus--Contributors to Cellular Detoxification

Detoxification is a fundamental cellular stress defense mechanism, which allows an organism to survive or even thrive in the presence of environmental toxins and/or pollutants. The glutathione S-transferase (GST) superfamily is a set of enzymes involved in the detoxification process. This highly diverse protein superfamily is characterized by multiple gene duplications, with over 40 GST genes reported in some insects. However, less is known about the GST superfamily in marine organisms, including crustaceans. The availability of two de novo transcriptomes for the copepod, Calanus finmarchicus, provided an opportunity for an in depth study of the GST superfamily in a marine crustacean. The transcriptomes were searched for putative GST-encoding transcripts using known GST proteins from three arthropods as queries. The identified transcripts were then translated into proteins, analyzed for structural domains, and annotated using reciprocal BLAST analysis. Mining the two transcriptomes yielded a total of 41 predicted GST proteins belonging to the cytosolic, mitochondrial or microsomal classes. Phylogenetic analysis of the cytosolic GSTs validated their annotation into six different subclasses. The predicted proteins are likely to represent the products of distinct genes, suggesting that the diversity of GSTs in C. finmarchicus exceeds or rivals that described for insects. Analysis of relative gene expression in different developmental stages indicated low levels of GST expression in embryos, and relatively high expression in late copepodites and adult females for several cytosolic GSTs. A diverse diet and complex life history are factors that might be driving the multiplicity of GSTs in C. finmarchicus, as this copepod is commonly exposed to a variety of natural toxins. Hence, diversity in detoxification pathway proteins may well be key to their survival.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

A Halloween gene noppera-bo encodes a glutathione S-transferase essential for ecdysteroid biosynthesis via regulating the behaviour of cholesterol in Drosophila

In insects, the precise timing of moulting and metamorphosis is strictly guided by ecdysteroids that are synthesised from dietary cholesterol in the prothoracic gland (PG). In the past decade, several ecdysteroidogenic enzymes, some of which are encoded by the Halloween genes, have been identified and characterised. Here, we report a novel Halloween gene, noppera-bo (nobo), that encodes a member of the glutathione S-transferase family. nobo was identified as a gene that is predominantly expressed in the PG of the fruit fly Drosophila melanogaster. We generated a nobo knock-out mutant, which displayed embryonic lethality and a naked cuticle structure. These phenotypes are typical for Halloween mutants showing embryonic ecdysteroid deficiency. In addition, the PG-specific nobo knock-down larvae displayed an arrested phenotype and reduced 20-hydroxyecdysone (20E) titres. Importantly, both embryonic and larval phenotypes were rescued by the administration of 20E or cholesterol. We also confirm that PG cells in nobo loss-of-function larvae abnormally accumulate cholesterol. Considering that cholesterol is the most upstream material for ecdysteroid biosynthesis in the PG, our results raise the possibility that nobo plays a crucial role in regulating the behaviour of cholesterol in steroid biosynthesis in insects.

Parallel evolution or purifying selection, not introgression, explains similarity in the pyrethroid detoxification linked GSTE4 of Anopheles gambiae and An. arabiensis

Insecticide resistance is a major impediment to the control of vectors and pests of public health importance and is a strongly selected trait capable of rapid spread, sometimes even between closely related species. Elucidating the mechanisms generating insecticide resistance in mosquito vectors of disease, and understanding the spread of resistance within and between populations and species are vital for the development of robust resistance management strategies. Here, we studied the mechanisms of resistance in two sympatric members of the Anopheles gambiae species complex-the major vector of malaria in sub-Saharan Africa-to understand how resistance has developed and spread in eastern Uganda, a region with some of the highest levels of malaria. In eastern Uganda, where the mosquitoes Anopheles arabiensis and An. gambiae can be found sympatrically, low levels of hybrids (0.4 %) occur, offering a route for introgression of adaptively important variants between species. In independent microarray studies of insecticide resistance, Gste4, an insect-specific glutathione S-transferase, was among the most significantly up-regulated genes in both species. To test the hypothesis of interspecific introgression, we sequenced 2.3 kbp encompassing Gste4. Whilst this detailed sequencing ruled out introgression, we detected strong positive selection acting on Gste4. However, these sequences, followed by haplotype-specific qPCR, showed that the apparent up-regulation in An. arabiensis is a result of allelic variation across the microarray probe binding sites which artefactually elevates the gene expression signal. Thus, face-value acceptance of microarray data can be misleading and it is advisable to conduct a more detailed investigation of the causes and nature of such signal. The identification of positive selection acting on this locus led us to functionally express and characterise allelic variants of GSTE4. Although the in vitro data do not support a direct role for GSTE4 in metabolism, they do support a role for this enzyme in insecticide sequestration. Thus, the demonstration of a role for an up-regulated gene in metabolic resistance to insecticides should not be limited to simply whether it can metabolise insecticide; such a strict criterion would argue against the involvement of GSTE4 despite the weight of evidence to the contrary.