Genomic Analysis of Detoxification Supergene Families in the Mosquito Anopheles sinensis

The role of P450 enzymes in malaria and other vector-borne infectious diseases

Vector-borne infectious diseases are still an important global health problem. Malaria is the most important among them, mainly pediatric, life-threatening disease. Malaria and other vector-borne disorders caused by parasites, bacteria, and viruses have a strong impact on public health and significant economic costs. Most vector-borne diseases could be prevented by vector control, with attention to the ecological and biodiversity conservation aspects. Chemical control with pesticides and insecticides is widely used as a measure of prevention although increasing resistance to insecticides is a serious issue in vector control. Metabolic resistance is the most common mechanism and poses a big challenge. Insect enzyme systems, including monooxygenase CYP P450 enzymes, are employed by vectors mainly to metabolize insecticides thus causing resistance. The discovery and application of natural specific inhibitors/blockers of vector P450 enzymes as synergists for commonly used pesticides will contribute to the "greening" of insecticides. Besides vector CYPs, host CYP enzymes could also be exploited to fight against vector-borne diseases: using mostly their detoxifying properties and involvement in the immune response. Here, we review published research data on P450 enzymes from all players in vector-borne infections, that is, pathogens, vectors, and hosts, regarding the potential role of CYPs in disease. We discuss strategies on how to exploit cytochromes P450 in vector-borne disease control.

Unravelling the genome of the brackish water malaria vector Anopheles aquasalis

Malaria is a severe public health problem in several developing tropical and subtropical countries. Anopheles aquasalis is the primary coastal malaria vector in Central and South America and the Caribbean Islands, and it has the peculiar feature of living in water with large changes in salinity. Recent research has recognised An. aquasalis as an important model for studying the interactions of murine and human Plasmodium parasites. This study presents the complete genome of An. aquasalis and offers insights into its evolution and physiology. The genome is similar in size and gene content to other Neotropical anophelines, with 162 Mb and 12,446 protein-coding genes. There are 1387 single-copy orthologs at the Diptera level (eg. An. gambiae, An. darlingi and Drosophila melanogaster). An. aquasalis diverged from An. darlingi, the primary malaria vector in inland South America, nearly 20 million years ago. Proteins related to ion transport and metabolism belong to the most abundant gene families with 660 genes. We identified gene families relevant to osmosis control (e.g., aquaporins, vacuolar-ATPases, Na+/K+-ATPases, and carbonic anhydrases). Evolutionary analysis suggests that all osmotic regulation genes are under strong purifying selection. We also observed low copy number variation in insecticide resistance and immunity-related genes for all known classical pathways. The data provided by this study offers candidate genes for further studies of parasite-vector interactions and for studies on how anophelines of brackish water deal with the high fluctuation in water salinity. We also established data and insights supporting An. aquasalis as an emerging Neotropical malaria vector model for genetic and molecular studies.

Glutathione S-transferase (GST) genes and their function associated with pyrethroid resistance in the malaria vector Anopheles sinensis

BACKGROUND

Glutathione S-transferases (GSTs), a multifunctional protein family, are involved in insecticide resistance. However, a systematic analysis of GSTs in Anopheles sinensis, an important vector for malaria transmission, is lacking. In this study, we investigated the diversity and characteristics of GST genes, and analyzed their expression patterns and functions associated with insecticide resistance in this species.

RESULTS

We identified 32 putative cytosolic and three putative microsomal GST genes in the An. sinensis genome. Transcriptome analysis showed that GSTs were highly expressed in larvae, and mainly expressed in the antennae, midgut and Malpighian tubules of adults. In addition, we found that GSTd2 and GSTe2 were significantly upregulated in four An. sinensis pyrethroid-resistant field populations. Furthermore, silencing of GSTd2 and GSTe2 significantly increased the susceptibility of An. sinensis to deltamethrin, and recombinant GSTd2 and GSTe2 exhibited high enzymatic activity in the metabolism of 1-chloro-2,4-dinitrobenzene and dichlorodiphenyltrichloroethane (DDT).

CONCLUSION

These results showed that GSTs are involved in the development of insecticide resistance in An. sinensis through transcriptional overexpression and enzymatic metabolization, facilitating our understanding of insecticide resistance in insects. © 2022 Society of Chemical Industry.

Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236 M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolize pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets cotreated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid cotreated LLINs across African countries.

Host-pathogen interaction between Asian citrus psyllid and entomopathogenic fungus (Cordyceps fumosorosea) is regulated by modulations in gene expression, enzymatic activity and HLB-bacterial population of the host

The host-pathogen interaction has been explored by several investigations, but the impact of fungal pathogens against insect resistance is still ambiguous. Therefore, we assessed the enzymatic activity and defense-related gene expression of Asian citrus psyllid (ACP) nymphal and adult populations on Huanglongbing-diseased citrus plants under the attack of Cordyceps fumosorosea. Overall, five enzymes viz. superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), glutathione S-transferase (GST), carboxylesterase (CarE), and four genes, namely SOD, 16S, CYP4C68, CYP4BD1, were selected for respective observations from ACP populations. Enzymatic activity of four enzymes (SOD, POD, GST, CarE) was significantly decreased after 5-days post-treatment (dpt) and 3-dpt fungal exposure in fungal treated ACP adult and nymphal populations, respectively, whereas the activity of CAT was boosted substantially post-treatment time schedule. Besides, we recorded drastic fluctuations in the expression of CYP4 genes among fungal treated ACP populations. After 24 hours post-treatment (hpt), expression of both CYP4 genes was boosted in fungal treated populations than controlled populations (adult and nymph). After 3-dpt, however, the expression of CYP4 genes was declined in the given populations. Likewise, fungal attack deteriorated the resistance of adult and nymphal of ACP population, as SOD expression was down-regulated in fungal-treated adult and nymphs after 5-dpt and 3-dpt exposure, respectively. Moreover, bacterial expression via the 16S gene was significantly increased in fungal-treated adult and nymphal ACP populations with increasing post-treatment time. Overall, our data illustrate that the fungal application disrupted the insect defense system. The expression of these genes and enzymes suppress the immune function of adult and nymphal ACP populations. As it is reported first time that the applications of C. fumosorosea against ACP reduce insect resistance by interfering with the CYP4 and SOD system. Therefore, we propose new strategies to discover the role of certain toxic compounds from fungus, which can reduce insect resistance, focusing on resistance-related genes and enzymes.

Chromosome-level genome assembly of the mirid predator Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), an important natural enemy in the rice ecosystem

Though the genomes of many rice herbivorous pests have recently been well characterized, little is known about the genome of their natural enemies. Here, by using the Illumina and PacBio platforms, we sequenced and assembled the whole genome of the mirid species Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), which is an economically and ecologically important natural enemy in the rice ecosystem acting as a dominant predator for planthoppers and leafhoppers in the field. Through Hi-C scaffolding, 1615 scaffolds with a total size of 338.08 Mb were successfully anchored onto 13 chromosomes. The assembled genome size was 345.75 Mb with a final scaffold N50 of 27.58 Mb. Approximately 107.51 Mb of sequences accounting for 31.10% of the genome were identified as repeat elements, and 14,644 protein-coding genes were annotated. Phylogenetic analysis showed that C. lividipennis clustered with other Hemipteran species and diverged from Apolygus lucorum about 66.7 million years ago. Gene families related to detoxification, environmental adaptation and digestion were analysed comparatively with other Hemipteran species, but no significant expansion or contraction was found in C. lividipennis. We also observed male meiosis in C. lividipennis, which showed a typical post-reduction of sex chromosomes and a karyotype of 2n = 22 + XY. As the first natural-enemy genome in the rice ecosystem, the genomic resource of C. lividipennis not only expands our understanding of the multitrophic interactions (host plant-prey-predator), but also provides a genomic basis for better understanding this dominant predator and therefore promotes sustainable rice pest management and food grain production.

Anopheles metabolic proteins in malaria transmission, prevention and control: a review

The increasing resistance to currently available insecticides in the malaria vector, Anopheles mosquitoes, hampers their use as an effective vector control strategy for the prevention of malaria transmission. Therefore, there is need for new insecticides and/or alternative vector control strategies, the development of which relies on the identification of possible targets in Anopheles. Some known and promising targets for the prevention or control of malaria transmission exist among Anopheles metabolic proteins. This review aims to elucidate the current and potential contribution of Anopheles metabolic proteins to malaria transmission and control. Highlighted are the roles of metabolic proteins as insecticide targets, in blood digestion and immune response as well as their contribution to insecticide resistance and Plasmodium parasite development. Furthermore, strategies by which these metabolic proteins can be utilized for vector control are described. Inhibitors of Anopheles metabolic proteins that are designed based on target specificity can yield insecticides with no significant toxicity to non-target species. These metabolic modulators combined with each other or with synergists, sterilants, and transmission-blocking agents in a single product, can yield potent malaria intervention strategies. These combinations can provide multiple means of controlling the vector. Also, they can help to slow down the development of insecticide resistance. Moreover, some metabolic proteins can be modulated for mosquito population replacement or suppression strategies, which will significantly help to curb malaria transmission.

Genome of the webworm Hyphantria cunea unveils genetic adaptations supporting its rapid invasion and spread

BACKGROUND

The fall webworm Hyphantria cunea is an invasive and polyphagous defoliator pest that feeds on nearly any type of deciduous tree worldwide. The silk web of H. cunea aids its aggregating behavior, provides thermal regulation and is regarded as one of causes for its rapid spread. In addition, both chemosensory and detoxification genes are vital for host adaptation in insects.

RESULTS

Here, a high-quality genome of H. cunea was obtained. Silk-web-related genes were identified from the genome, and successful silencing of the silk protein gene HcunFib-H resulted in a significant decrease in silk web shelter production. The CAFE analysis showed that some chemosensory and detoxification gene families, such as CSPs, CCEs, GSTs and UGTs, were expanded. A transcriptome analysis using the newly sequenced H. cunea genome showed that most chemosensory genes were specifically expressed in the antennae, while most detoxification genes were highly expressed during the feeding peak. Moreover, we found that many nutrient-related genes and one detoxification gene, HcunP450 (CYP306A1), were under significant positive selection, suggesting a crucial role of these genes in host adaptation in H. cunea. At the metagenomic level, several microbial communities in H. cunea gut and their metabolic pathways might be beneficial to H. cunea for nutrient metabolism and detoxification, and might also contribute to its host adaptation.

CONCLUSIONS

These findings explain the host and environmental adaptations of H. cunea at the genetic level and provide partial evidence for the cause of its rapid invasion and potential gene targets for innovative pest management strategies.

Proteomic Investigation on Anopheles gambiae in Burkina Faso Related to Insecticide Pressures from Different Climatic Regions

In Sub-Saharan Africa, An. gambiae sensu lato (s.l.) Giles 190, largely contributes to malaria transmission. Therefore, the authors carry out a proteomic analysis to compare its metabolic state, depending on different pesticide pressures by selecting areas with/without cotton crops. The proteomes data are available via ProteomeXchange with identifier PXD016300. From a total of 1.182 identified proteins, 648 are retained for further statistical analysis and are attributed to biological functions, the most important of which being energy metabolism (120 proteins) followed by translation-biogenesis (74), cytoskeleton (71), stress response (62), biosynthetic process (60), signalling (44), cellular respiration (38), cell redox homeostasis (25), DNA processing (17), pheromone binding (10), protein folding (9), RNA processing (9), other proteins (26) and unknown functions (83). In the Sudano-Sahelian region, 421 (91.3%) proteins are found in samples from areas both with and without cotton crops. By contrast, in the Sahelian region, only 271 (55.0%) are common to both crop areas, and 233 proteins are up-regulated from the cotton area. The focus is placed on proteins with putative roles in insecticide resistance, according to literature. This study provides the first whole-body proteomic characterisation of An. gambiae s.l. in Burkina Faso, as a framework to strengthen vector control strategies.

Characterisation of GST genes from the Hyphantria cunea and their response to the oxidative stress caused by the infection of Hyphantria cunea nucleopolyhedrovirus (HcNPV)

The fall webworm, Hyphantria cunea (Drury) (Lepidoptera: Noctuidae), is a major pest found in forests. In this study, the effects of Hyphantria cunea nucleopolyhedrovirus (HcNPV) infection on the transcription levels and activities of glutathione S-transferases (GSTs) in H. cunea were determined. In the present study, 18 GST family genes were identified from the H. cunea transcriptome dataset by using bioinformatic analyses. These GST genes were classified into cytosolic (15 genes) and microsomal (three genes) classes. The 15 cytosolic GST genes belonged to four different subclasses (epsilon, sigma and delta). The all GST genes, especially GSTe4, showed high expression levels in egg and 1st~4th instar larval stage while their low expression levels in 5th~7th instar larvae using real-time quantitative PCR analysis. However, the expression levels of the 18 GST genes were varied after exposure to sublethal doses of HcNPV. The expression levels of most GSTs were downregulated and upregulated at low and high concentrations of HcNPV, respectively. The corresponding total GST activities also showed similar patterns. In H. cunea, changes in the expression levels and enzymatic activities of GSTs after exposure to HcNPV indicated that they may have important functions in the defense against HcNPV, and the stress, which may be reflected by the high GST enzymatic activities.

De novo transcriptome analysis and identification of genes associated with immunity, detoxification and energy metabolism from the fat body of the tephritid gall fly, Procecidochares utilis

The fat body, a multifunctional organ analogous to the liver and fat tissue of vertebrates, plays an important role in insect life cycles. The fat body is involved in protein storage, energy metabolism, elimination of xenobiotics, and production of immunity regulator-like proteins. However, the molecular mechanism of the fat body's physiological functions in the tephritid stem gall-forming fly, Procecidochares utilis, are still unknown. In this study, we performed transcriptome analysis of the fat body of P. utilis using Illumina sequencing technology. In total, 3.71 G of clean reads were obtained and assembled into 30,559 unigenes, with an average length of 539 bp. Among those unigenes, 21,439 (70.16%) were annotated based on sequence similarity to proteins in NCBI's non-redundant protein sequence database (Nr). Sequences were also compared to NCBI's non-redundant nucleotide sequence database (Nt), a manually curated and reviewed protein sequence database (SwissProt), and KEGG and gene ontology annotations were applied to better understand the functions of these unigenes. A comparative analysis was performed to identify unigenes related to detoxification, immunity and energy metabolism. Many unigenes involved in detoxification were identified, including 50 unigenes of putative cytochrome P450s (P450s), 18 of glutathione S-transferases (GSTs), 35 of carboxylesterases (CarEs) and 26 of ATP-binding cassette (ABC) transporters. Many unigenes related to immunity were identified, including 17 putative serpin genes, five peptidoglycan recognition proteins (PGRPs) and four lysozyme genes. In addition, unigenes potentially involved in energy metabolism, including 18 lipase genes, five fatty acid synthase (FAS) genes and six elongases of very long chain fatty acid (ELOVL) genes, were identified. This transcriptome improves our genetic understanding of P. utilis and the identification of a numerous transcripts in the fat body of P. utilis offer a series of valuable molecular resources for future studies on the functions of these genes.

UDP-glycosyltransferase genes and their association and mutations associated with pyrethroid resistance in Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

UDP-glycosyltransferase (UGT) is an important biotransformation superfamily of enzymes. They catalyze the transfer of glycosyl residues from activated nucleotide sugars to acceptor hydrophobic molecules, and function in several physiological processes, including detoxification, olfaction, cuticle formation, pigmentation. The diversity, classification, scaffold location, characteristics, phylogenetics, and evolution of the superfamily of genes at whole genome level, and their association and mutations associated with pyrethroid resistance are still little known.

METHODS

The present study identified UGT genes in Anopheles sinensis genome, classified UGT genes in An. sinensis, Anopheles gambiae, Aedes aegypti and Drosophila melanogaster genomes, and analysed the scaffold location, characteristics, phylogenetics, and evolution of An. sinensis UGT genes using bioinformatics methods. The present study also identified the UGTs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and the mutations associated with pyrethroid resistance with genome re-sequencing in An. sinensis.

RESULTS

There are 30 putative UGTs in An. sinensis genome, which are classified into 12 families (UGT301, UGT302, UGT306, UGT308, UGT309, UGT310, UGT313, UGT314, UGT315, UGT36, UGT49, UGT50) and further into 23 sub-families. The UGT308 is significantly expanded in gene number compared with other families. A total of 119 UGTs from An. sinensis, An. gambiae, Aedes aegypti and Drosophila melanogaster genomes are classified into 19 families, of which seven are specific for three mosquito species and seven are specific for Drosophila melanogaster. The UGT308 and UGT302 are proposed to main families involved in pyrethroid resistance. The AsUGT308D3 is proposed to be the essential UGT gene for the participation in biotransformation in pyrethroid detoxification process, which is possibly regulated by eight SNPs in its 3' flanking region. The UGT302A3 is also associated with pyrethroid resistance, and four amino acid mutations in its coding sequences might enhance its catalytic activity and further result in higher insecticide resistance.

CONCLUSIONS

This study provides the diversity, phylogenetics and evolution of UGT genes, and potential UGT members and mutations involved in pyrethroid resistance in An. sinensis, and lays an important basis for the better understanding and further research on UGT function in defense against insecticide stress.

Genome-wide identification, characterization and evolution of cuticular protein genes in the malaria vector Anopheles sinensis (Diptera: Culicidae)

Thirteen cuticular protein (CP) families have been recognized in arthropods. In this study, 250 Anopheles sinensis CP genes were identified and named based on genome and transcriptome sequences. They were classified into 10 families based on motifs and phylogenetic analyses. In 11 other insect species, nine had CP numbers > 150 while Apis mellifera and Tribolium castaneum had CP numbers less than 52. The CPs of eight species occupied > 1.4% of the total genomic gene number, whereas in three species the CPs occupied < 1%. The phylogenies for each CP family in An. sinensis were constructed and discussed. The 250 CPs each had 1-8 exons with 144 CPs (57.6%) having two exons. The intron length ranged from 66-3888 bp with 174 introns (54.0%) being 66-100 bp long. Except for two CPs on two contigs, 248 CPs were mapped onto 28 scaffolds with 136 genes (54.4%) restricted to five scaffolds. A total of 107 CPs were clustered and located at 27 loci. The CPR family had the conserved motif GSYSLVEPDGTVRTV. The RR-1 subfamily had an additional 21 amino acid (aa) motifs with the YVADENGF sequence that is common in insects. The RR-2 subfamily had an additional 50 aa motifs with two additional regions RDGDVVKG and G-x(3)-VV. A comparison with 115 orthologous counterparts of An. gambiae CPs suggested purifying selection for all of these genes. This study provides basic information useful for further studies on biological functions of An. sinensis CPs as well as for comparative genomics of insect CPs.

De novo transcriptomic analysis of the alimentary tract of the tephritid gall fly, Procecidochares utilis

The tephritid gall fly, Procecidochares utilis, is an important obligate parasitic insect of the malignant weed Eupatorium adenophorum which biosynthesizes toxic secondary metabolites. Insect alimentary tracts secrete several enzymes that are used for detoxification, including cytochrome P450s, glutathione S-transferases, and carboxylesterases. To explore the adaptation of P. utilis to its toxic host plant, E. adenophorum at molecular level, we sequenced the transcriptome of the alimentary tract of P. utilis using Illumina sequencing. Sequencing and de novo assembly yielded 62,443 high-quality contigs with an average length of 604 bp that were further assembled into 45,985 unigenes with an average length of 674 bp and an N50 of 983 bp. Among the unigenes, 30,430 (66.17%) were annotated by alignment against the NCBI non-redundant protein (Nr) database, while 16,700 (36.32%), 16,267 (35.37%), and 11,530 (25.07%) were assigned functions using the Clusters of Orthologous Groups (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) databases, respectively. Using the comprehensive transcriptome data set, we manually identified several important gene families likely to be involved in the detoxification of toxic compounds including 21 unigenes within the glutathione S-transferase (GST) family, 22 unigenes within the cytochrome P450 (P450) family, and 16 unigenes within the carboxylesterase (CarE) family. Quantitative PCR was used to verify eight, six, and two genes of GSTs, P450s, and CarEs, respectively, in different P. utilis tissues and at different developmental stages. The detoxification enzyme genes were mainly expressed in the foregut and midgut. Moreover, the unigenes were higher expressed in the larvae, pupae, and 3-day adults, while they were expressed at lower levels in eggs. These transcriptomic data provide a valuable molecular resource for better understanding the function of the P. utilis alimentary canal. These identified genes could be pinpoints to address the molecular mechanisms of P. utilis interacting with toxic plant host.

ASGDB: a specialised genomic resource for interpreting Anopheles sinensis insecticide resistance

BACKGROUND

Anopheles sinensis is an important malaria vector in Southeast Asia. The widespread emergence of insecticide resistance in this mosquito species poses a serious threat to the efficacy of malaria control measures, particularly in China. Recently, the whole-genome sequencing and de novo assembly of An. sinensis (China strain) has been finished. A series of insecticide-resistant studies in An. sinensis have also been reported. There is a growing need to integrate these valuable data to provide a comprehensive database for further studies on insecticide-resistant management of An. sinensis.

RESULTS

A bioinformatics database named An. sinensis genome database (ASGDB) was built. In addition to being a searchable database of published An. sinensis genome sequences and annotation, ASGDB provides in-depth analytical platforms for further understanding of the genomic and genetic data, including visualization of genomic data, orthologous relationship analysis, GO analysis, pathway analysis, expression analysis and resistance-related gene analysis. Moreover, ASGDB provides a panoramic view of insecticide resistance studies in An. sinensis in China. In total, 551 insecticide-resistant phenotypic and genotypic reports on An. sinensis distributed in Chinese malaria-endemic areas since the mid-1980s have been collected, manually edited in the same format and integrated into OpenLayers map-based interface, which allows the international community to assess and exploit the high volume of scattered data much easier. The database has been given the URL: http://www.asgdb.org /.

CONCLUSIONS

ASGDB was built to help users mine data from the genome sequence of An. sinensis easily and effectively, especially with its advantages in insecticide resistance surveillance and control.

Carbonic anhydrase II confers resistance to deltamethrin in Culex pipiens pallens

Carbonic anhydrases (CAs) are metabolic enzymes that regulate the physiological equilibrium in a variety of organisms. In this study, we investigated the function of CA II in Culex pipiens pallens using real-time quantitative polymerase chain reaction and double-stranded RNA injection. CA II transcripts were more abundant in the deltamethrin-susceptible strain than the deltamethrin-resistant strain. The activities of metabolic enzymes increased when the CA II expression was silenced. These findings suggest CA II regulates deltamethrin resistance by altering metabolic enzyme activity, and could serve as a potential genetic marker for monitoring deltamethrin resistance in mosquitoes.

Biology, Bionomics and Molecular Biology of Anopheles sinensis Wiedemann 1828 (Diptera: Culicidae), Main Malaria Vector in China

China has set a goal to eliminate all malaria in the country by 2020, but it is unclear if current understanding of malaria vectors and transmission is sufficient to achieve this objective. Anopheles sinensis is the most widespread malaria vector specie in China, which is also responsible for vivax malaria outbreak in central China. We reviewed literature from 1954 to 2016 on An. sinensis with emphasis on biology, bionomics, and molecular biology. A total of 538 references were relevant and included. An. sienesis occurs in 29 Chinese provinces. Temperature can affect most life-history parameters. Most An. sinensis are zoophilic, but sometimes they are facultatively anthropophilic. Sporozoite analysis demonstrated An. sinensis efficacy on Plasmodium vivax transmission. An. sinensis was not stringently refractory to P. falciparum under experimental conditions, however, sporozoite was not found in salivary glands of field collected An. sinensis. The literature on An. sienesis biology and bionomics was abundant, but molecular studies, such as gene functions and mechanisms, were limited. Only 12 molecules (genes, proteins or enzymes) have been studied. In addition, there were considerable untapped omics resources for potential vector control tools. Existing information on An. sienesis could serve as a baseline for advanced research on biology, bionomics and genetics relevant to vector control strategies.

piRNA-3878 targets P450 (CpCYP307B1) to regulate pyrethroid resistance in Culex pipiens pallens

Piwi-interacting RNAs (piRNAs) are a novel class of noncoding single-strand RNAs. They play an important role in the germ cell maintenance, brain development, epigenetic regulation of cancer, and antiviral function. However, little is known about the relationship between the piRNAs and insecticide resistance in mosquitoes. In this study, we reported that piRNA-3878 was related with pyrethroid resistance in Culex pipiens pallens. The expression level of piRNA-3878 was lower in both laboratory and field-collected deltamethrin-resistant (DR) strains. Overexpression of piRNA-3878 increased the susceptibility of the DR strain, while inhibiting the expression of piRNA-3878 in DS strain made the mosquitoes more resistant to deltamethrin. Furthermore, we identified that CpCYP307B1 was the target of piRNA-3878. The mosquito mortality rate was increased after downregulating the expression of CpCYP307B1. These findings revealed that piRNA-3878 could target CpCYP307B1 to regulate pyrethroid resistance.

Insecticide resistance in malaria vectors along the Thailand-Myanmar border

BACKGROUND

There is a paucity of data about the susceptibility status of malaria vectors to Public Health insecticides along the Thailand-Myanmar border. This lack of data is a limitation to guide malaria vector-control in this region. The aim of this study was to assess the susceptibility status of malaria vectors to deltamethrin, permethrin and DDT and to validate a simple molecular assay for the detection of knock-down resistance (kdr) mutations in the study area.

METHODS

Anopheles mosquitoes were collected in four sentinel villages during August and November 2014 and July 2015 using human landing catch and cow bait collection methods. WHO susceptibility tests were carried out to measure the mortality and knock-down rates of female mosquitoes to deltamethrin (0.05%), permethrin (0.75%) and DDT (4%). DNA sequencing of a fragment of the voltage-gated sodium channel gene was carried out to identify knock-down resistance (kdr) mutations at position 1014 in mosquitoes surviving exposure to insecticides.

RESULTS

A total of 6295 Anopheles belonging to ten different species were bioassayed. Resistance or suspected resistance to pyrethroids was detected in An. barbirostris (s.l.) (72 and 84% mortality to deltamethrin (n = 504) and permethrin (n = 493) respectively), An. hyrcanus (s.l.) (33 and 48% mortality to deltamethrin (n = 172) and permethrin (n = 154), respectively), An. jamesii (87% mortality to deltamethrin, n = 111), An. maculatus (s.l.) (85 and 97% mortality to deltamethrin (n = 280) and permethrin (n = 264), respectively), An. minimus (s.l.) (92% mortality, n = 370) and An. vagus (75 and 95% mortality to deltamethrin (n =148) and permethrin (n = 178), respectively). Resistance or suspected resistance to DDT was detected in An. barbirostris (s.l.) (74% mortality, n = 435), An. hyrcanus (s.l.) (57% mortality, n = 91) and An. vagus (97% mortality, n = 133). The L1014S kdr mutation at both heterozygous and homozygous state was detected only in An. peditaeniatus (Hyrcanus Group).

CONCLUSION

Resistance to pyrethroids is present along the Thailand-Myanmar border, and it represents a threat for malaria vector control. Further investigations are needed to better understand the molecular basis of insecticide resistance in malaria vectors in this area.

A neural network prediction of environmental determinants of Anopheles sinensis knockdown resistance mutation to pyrethroids in China

Selection pressure caused by long-term intensive use of insecticides is the key driving force in resistance development. Additional parameters such as environmental conditions may affect both the mosquito response to insecticides and the selection of resistance mechanisms. In this context, we analyzed the environmental determinants of kdr prevalence in Anopheles sinensis across China. We collected kdr frequency from 48 sites across central and southern China, together with key environmental factors including long-term climatic data, topographic features, main crops, and land cover types. Trend surface analysis found that the distribution of kdr frequency can be partitioned into three regions, namely central China (kdr frequency >80%), western China (kdr frequency varies from 0% to 60%), and southern China (kdr frequency <10%). Seven predictor variables were selected based on a radial basis function neural network model. A multilayer perceptron (MLP) network model revealed that the number of crops in a year was the most important predictor for the kdr mutation rate. Topography, long-term mean climate and land cover all contributed to the kdr mutation rate. The observed mean kdr frequency was 53.0% and the MLP network model-predicted mean was 52.6%, a 0.1% relative error. Predicted kdr frequencies closely matched the observed values. The model explained 92% of the total variance in kdr frequency. The results indicated that kdr was associated with the intensity of pesticide usage. Crop cultivation information, together with environmental factors, may well predict the spatial heterogeneity of kdr mutations in An. sinensis in China.

Identification of glutathione S-transferase genes in Leptinotarsa decemlineata and their expression patterns under stress of three insecticides

Glutathione S-transferases (GSTs) is a family of multifunctional enzymes that are involved in detoxification of poisonous compounds. In the present paper, the Leptinotarsa decemlineata genome and transcriptome dataset were mined and 30 GST genes were identified. These GSTs belonged to cytosolic (29 genes) and microsomal (1 gene) classes. Among them 3 GSTs (LdGSTe2, LdGSTs4, and LdGSTo3) possessed splice variants. Of the 29 cytosolic LdGSTs, 3, 10, 5, 4, 4, and 1 members were classified as delta, epsilon, omega, sigma, theta, and zeta subclasses respectively, along with 2 unclassified genes. Phylogenetic analysis suggest that epsilon, omega and sigma subclasses appear to undergo species-specific bloom. Moreover, most epsilon, omega and sigma GSTs are tandemly arranged in three chromosome scaffolds. To find GST candidates involving in insecticide detoxification, we tested the mRNA levels of 20 GST transcripts under stress of cyhalothrin, fipronil or endosulfan. Out of them, LdGSTe2a, LdGSTe2b, LdGSTo5 and LdGSTt1 were significantly overexpressed after exposure to each of the three insecticides. Two other genes were respectively upregulated after cyhalothrin (LdGSTe10 and LdGSTu2) or endosulfan (LdGSTd1 and LdGSTu2) treatment. The diversified expression responses to insecticide exposure suggest that the LdGSTs may depend on a functionally complex system to detoxify different classes of insecticides. In addition, our findings provide a base for a better understanding of the evolution of insecticide resistance, and functional research on specific GST genes.