Microarray and RNAi analysis of P450s in Anopheles gambiae male and female steroidogenic tissues: CYP307A1 is required for ecdysteroid synthesis

Pyrethroid resistance and gene expression profile of a new resistant An. gambiae colony from Uganda reveals multiple resistance mechanisms and overexpression of Glutathione-S-Transferases linked to survival of PBO-pyrethroid combination

Background

The effectiveness of long-lasting insecticidal nets (LLINs) are being threatened by growing resistance to pyrethroids. To restore their efficacy, a synergist, piperonyl butoxide (PBO) which inhibits cytochrome P450s has been incorporated into pyrethroid treated nets. A trial of PBO-LLINs was conducted in Uganda from 2017 and we attempted to characterize mechanisms of resistance that could impact intervention efficacy.

Methods

We established an Anopheles gambiae s.s colony in 2018 using female mosquitoes collected from Busia district in eastern Uganda. We first assessed the phenotypic resistance profile of this colony using WHO tube and net assays using a deltamethrin dose-response approach. The Busia colony was screened for known resistance markers and RT-qPCR targeting 15 genes previously associated with insecticide resistance was performed.

Results

The Busia colony had very high resistance to deltamethrin, permethrin and DDT. In addition, the colony had moderate resistance to alpha-cypermethrin and lambda-cyhalothrin but were fully susceptible to bendiocarb and fenitrothion. Exposure to PBO in combination with permethrin and deltamethrin resulted in higher mortality rates in both net and tube assays, with a higher mortality observed in net assays than tube assays. The kdr marker, Vgsc-995S was at very high frequency (91.7-98.9%) whilst the metabolic markers Coeae1d and Cyp4j5-L43F were at very low (1.3% - 11.5%) and moderate (39.5% - 44.7%) frequencies respectively. Our analysis showed that gene expression pattern in mosquitoes exposed to deltamethrin, permethrin or DDT only were similar in comparison to the susceptible strain and there was significant overexpression of cytochrome P450s, glutathione-s-transferases (GSTs) and carboxyl esterases (COEs). However, mosquitoes exposed to both PBO and pyrethroid strikingly and significantly only overexpressed closely related GSTs compared to unexposed mosquitoes while major cytochrome P450s were underexpressed.

Conclusions

The high levels of pyrethroid resistance observed in Busia appears associated with a wide range of metabolic gene families.

Time-series analysis of transcriptomic changes due to permethrin exposure reveals that Aedes aegypti undergoes detoxification metabolism over 24 h

Insecticide resistance is a multifaceted response and an issue across taxa. Aedes aegypti, the mosquito that vectors Zika, dengue, chikungunya, and yellow fever, demonstrates high levels of pyrethroid resistance across the globe, presenting a challenge to public health officials. To examine the transcriptomic shifts across time after exposure to permethrin, a 3'Tag-Seq analysis was employed on samples 6, 10, and 24 h after exposure along with controls. Differential expression analysis revealed significant shifts in detoxifying enzymes and various energy-producing metabolic processes. These findings indicate significant alterations in gene expression associated with key energy mobilization pathways within the system. These changes encompass a coordinated response involving lipolysis, beta-oxidation, and the citric acid cycle, required for the production of energetic molecules such as ATP, NADH, NADPH, and FADH. These findings highlight a complex interplay of metabolic processes that may have broader implications for understanding insect physiology and response to environmental stimuli. Among the upregulated detoxifying enzymes are cytochrome P450s, glutathione s-transferases and peroxidases, and ATP-binding cassette transporters. Additionally, eight heat shock genes or genes with heat shock domains exhibit the highest fold change across time. Twenty-four hours after exposure, samples indicate a global downregulation of these processes, though principal component analysis suggests lasting signatures of the response. Understanding the recovery response to insecticide exposure provides information on possible new genetic and synergist targets to explore.

Functional Characterization of Drosophila melanogaster CYP6A8 Fatty Acid Hydroxylase

Genomic analysis indicated that the genome of Drosophila melanogaster contains more than 80 cytochrome P450 genes. To date, the enzymatic activity of these P450s has not been extensively studied. Here, the biochemical properties of CYP6A8 were characterized. CYP6A8 was cloned into the pCW vector, and its recombinant enzyme was expressed in Escherichia coli and purified using Ni²⁺-nitrilotriacetate affinity chromatography. Its expression level was approximately 130 nmol per liter of culture. Purified CYP6A8 exhibited a low-spin state in the absolute spectra of the ferric forms. Binding titration analysis indicated that lauric acid and capric acid produced type І spectral changes, with K_d values 28 ± 4 and 144 ± 20 μM, respectively. Ultra-performance liquid chromatography-mass spectrometry analysis showed that the oxidation reaction of lauric acid produced (ω-1)-hydroxylated lauric acid as a major product and ω-hydroxy-lauric acid as a minor product. Steady-state kinetic analysis of lauric acid hydroxylation yielded a k_cat value of 0.038 ± 0.002 min^-1 and a K_m value of 10 ± 2 μM. In addition, capric acid hydroxylation of CYP6A8 yielded kinetic parameters with a k_cat value of 0.135 ± 0.007 min^-1 and a K_m value of 21 ± 4 μM. Because of the importance of various lipids as carbon sources, the metabolic analysis of fatty acids using CYP6A8 in this study can provide an understanding of the biochemical roles of P450 enzymes in many insects, including Drosophila melanogaster.

Diallyl trisulfide reduced the reproductive capacity of male Sitotroga cerealella via the regulation of juvenile and ecdysone hormones

Environmental pollution and resistance in animals are major concerns for the application of synthetic pesticides. Diallyl trisulfide (DAT), an active compound in garlic essential oil, is a novel tool for active and safe control of agricultural insect pests. In this study, we analysed the effects of DAT (0.01 μL/L) on the protein content in male reproductive tissues (accessory glands, ejaculatory ducts, and testis), and juvenile hormone (JH) and ecdysone titres in a highly detrimental pest of stored products, Sitotroga cerealella. Evaluation of the expression profile of JH and ecdysone pathway-related genes in various tissues indicated that the accessory gland protein and ecdysone titres were markedly decreased after DAT fumigation, whereas the testis protein content and JH titre were increased. However, the protein content of the ejaculatory ducts remained unchanged between the treated and control groups. Further investigation revealed that DAT disrupted the mRNA expression of key enzymes involved in JH and ecdysone pathways. While increased mRNA levels of juvenile hormone acid O-methyltransferase (JHMAT) and Kruppel homologue 1 (Kr-h1) were observed after 4 and 7 h of DAT fumigation, the levels of juvenile hormone epoxide hydrolase (JHEH) were substantially reduced 3 h post-fumigation. mRNA levels of the ecdysone-responsive gene, FTZF1, and cytochrome P450 enzyme, CYP315A1, were notably decreased at 7 h and 4 h, respectively, post-fumigation, whereas CYP314A1 and CYP302A1 mRNA levels decreased after 3 h and 4 h, respectively. While DAT fumigation disrupted sperm number in the testis, ejaculatory ducts, and seminal vesicles, topical application of the 20-hydroxyecdysone (20E) analogue also lowered sperm number in the ejaculatory ducts. Topical application of methoprene, a JH analogue, increased the protein content in the testes, but not in the accessory glands or ejaculatory ducts. However, the survival rate was not affected by the topical application of methoprene or 20E. These data suggest that DAT regulates JH and ecdysone via its molecular pathway genes and modulates endocrine secretion during the male reproductive process.

Transcriptional Dynamics Induced by Diapause Hormone in the Silkworm, Bombyx mori

Diapause is a form of dormancy that organisms use to adapt to extreme environments by exhibiting developmental arrest. In the silkworm, Bombyx mori, diapause is thought to be elicited by diapause hormone (DH) signaling, which consists of interactions between DH and the DH receptor (DHR). However, the steps downstream of the DH signaling pathway are largely unknown. In the present study, we directly injected synthesized DH into the female pupae of a multivoltine, non-diapausing strain at 36 h after pupation. We found that the mRNA level of DHR declined at 4 h and recovered at 12 h after the injection of DH. Thus, we sequenced the transcriptome of the ovaries at 4 h and 12 h after the injection of DH. We identified 60 and 221 differentially expressed genes at 4 h and 12 h after the injection, respectively. All DEGs were identified, relating to 20E-related genes, JH-related genes, cellular detoxification, ribosomal proteins, lipid metabolism, and epigenetic modifications. Eleven genes were selected from the above categories to verify the transcriptome data. The qRT-PCR and RNA-Seq expression patterns of the genes were consistent, which indicated the authenticity and reliability of the transcriptome data. This study dramatically expands upon our knowledge of gene expression variation at the early phase of DH release.

Silencing of cytochrome P450 gene CYP321A1 effects tannin detoxification and metabolism in Spodoptera litura

Cytochrome P450 monooxygenase (P450 or CYP) plays an important role in the metabolism of insecticides and plant allelochemicals by insects. CYP321B1, a novel Spodoptera litura P450 gene, was identified and characterized. CYP321B1 contains a 1488 bp open reading frame (ORF) that encodes a 495 amino acid protein. In fourth instar larvae, the highest CYP321B1 expression levels were found in the midgut and fat body. In the tannin feeding test, tannin can significantly induce the expression of CYP321B1 in the midgut and fat body of 4th instar larvae. To verify the function of CYP321B1, RNA interference and metabolome analysis were performed. The results showed that silencing CYP321B1 significantly reduced the rate of weight gain under tannin induction. Metabolome analysis showed silencing affected 47 different metabolites, mainly involved in secondary metabolite biosynthesis and amino acid metabolism, including amino acids, lipid fatty acids, organic acids and their derivatives. Henoxyacetic acid and cysteamine are the most highly regulated metabolites, respectively. These findings demonstrate that CYP321B1 plays an important role in tannin detoxification and metabolism. Functional knowledge about metabolite detoxification genes in this major herbivorous insect pest can provide new insights into this biological process and provide new targets for agricultural pest control.

Genome-Wide Selective Signature Analysis Revealed Insecticide Resistance Mechanisms in Cydia pomonella

Simple Summary

The codling moth, Cydia pomonella, is a quarantine pest that causes extensive damage to many important pome fruits. To control this pest, insecticides are frequently used, leading to the development of resistance. In this study, we analyzed resequencing data of two resistant and one susceptible strains of codling moth, detecting the positively selected genes under the insecticide selective pressure. Coupled with transcriptome data, we discussed the potential role in insecticide resistance of these positively selected genes. Our results identified eight genes including CYP6b2, CYP307a1, 5-hydroxytryptaminereceptor, cuticle protein, and acetylcholinesterase, which are potentially involved in cross-resistance to azinphos-methyl and deltamethrin. Overall, our finding indicated that the insecticide resistance mechanism in C. pomonella is a complex physiological and biochemical process.

Abstract

The codling moth, Cydia pomonella L. (Lepidoptera, Tortricidae), is a serious invasive pest of pome fruits. Currently, C. pomonella management mainly relies on the application of insecticides, which have driven the development of resistance in the insect. Understanding the genetic mechanisms of insecticide resistance is of great significance for developing new pest resistance management techniques and formulating effective resistance management strategies. Using existing genome resequencing data, we performed selective sweep analysis by comparing two resistant strains and one susceptible strain of the insect pest and identified seven genes, among which, two (glycine receptor and glutamate receptor) were under strong insecticide selection, suggesting their functional importance in insecticide resistance. We also found that eight genes including CYP6B2, CYP307a1, 5-hydroxytryptamine receptor, cuticle protein, and acetylcholinesterase, are potentially involved in cross-resistance to azinphos-methyl and deltamethrin. Moreover, among several P450s identified as positively selected genes, CYP6B2, CYP4C1, and CYP4d2 showed the highest expression level in larva compared to other stages tested, and CYP6B2 also showed the highest expression level in midgut, supporting the roles they may play in insecticide metabolism. Our results provide several potential genes that can be studied further to advance understanding of complexity of insecticide resistance mechanisms in C. pomonella.

RNA interference shows that Spook, the precursor gene of 20-hydroxyecdysone (20E), regulates the molting of Macrobrachium nipponense

The aim of this study was to explore the function of the Mn-Spook gene, which was found in the ovary transcriptome of the Oriental river prawn (Macrobrachium nipponense). The Spook gene, which is the precursor gene of 20-hydroxyecdysone (20E), plays an important role in the process of molting in many arthropods, but its function in M. nipponense is unclear. We cloned the full-length Mn-Spook gene from the ovary of M. nipponense and found that it had the same conserved domains as the P450 gene of the Halloween family of genes. The Mn-Spook gene was highly expressed in ovary and gill tissue during the breeding period. During ovarian development, Mn-spook gene expression was highest at the nearly-ripe stage, and it also was highly expressed in the zoea developmental stage. Cellular localization analysis showed that Mn-Spook signals accumulated in the cytoplasmic membrane and nucleus of oocytes. Finally, we used RNA interference to evaluate the function of the Mn-Spook gene. Compared with the control group, in vivo injection of Mn-Spook dsRNA effectively downregulated the expression of Mn-Spook and the content of 20E. The molting frequency of M. nipponense in the experimental group also was significantly inhibited. These results demonstrated that the Mn-Spook gene played an important role in the molting process of M. nipponense.

Understanding molt control switches: Transcriptomic and expression analysis of the genes involved in ecdysteroidogenesis and cholesterol uptake pathways in the Y-organ of the blue crab, Callinectes sapidus

The crustacean molting process is regulated by an interplay of hormones produced by the eyestalk ganglia and Y-organs (YO). Molt-inhibiting hormone and crustacean hyperglycemic hormone released by the sinus gland of the eyestalk ganglia (EG) inhibit the synthesis and secretion of ecdysteroid by the YO, hence regulating hemolymph levels during the molt cycle. The purpose of this study is to investigate the ecdysteroidogenesis pathway, specifically genes linked to changes in ecdysteroid levels occurring at early premolt (ePM). To this end, a reference transcriptome based on YO, EG, and hepatopancreas was de novo assembled. Two genes (cholesterol 7-desaturase Neverland and cytochrome p450 307a1-like Spook) involved in ecdysteroidogenesis were identified from the YO transcriptome using sequence comparisons and transcript abundance. Two other candidates, Hormone receptor 4 and probable cytochrome p450 49a1 potentially involved in ecdysteroidogenesis were also identified. Since cholesterol is the ecdysteroid precursor, a putative cholesterol carrier (Apolipoprotein D-like) was also examined to understand if cholesterol uptake coincided with the increase in the ecdysteroid levels at the ePM stage. The expression level changes of the five candidate genes in the YO were compared between intermolt (IM) and induced ePM (iePM) stages using transcriptomic analysis. Expression analysis using qPCR were carried out at IM, iePM, and normal ePM. The increase in Spook and Neverland expression in the YO at the ePM was accompanied by a concomitant rise in ecdysteroid levels. The data obtained from iePM stage were congruent with those obtained from the normal ePM stage of intact control animals. The present findings support the role of Halloween genes in the ecdysteroidogenesis and molt cycle in the blue crab, Callinectes sapidus.

Comparative transcriptome analysis of different developmental stage of Bactrocera minax (Diptera: Tephritidae): Implication of the molecular basis of its obligatory diapause induction

The Chinese citrus fruit fly, Bactrocera minax (Enderlein), a devastating pest in citrus, has an obligatory diapause at the pupal stage. Although the physiology of obligatory diapause has been well described, the molecular mechanisms underlying the process remain unknow. Here we investigated the molecular mechanisms of obligatory diapause induction in B. minax using high-throughput RNA-Seq data from second-instar larva (2L), third-instar larva (3L) and pupa (P) stages. A total of 116,402 unigenes were obtained, of which 54,781 unigenes were successfully annotated in public databases, and the differentially expressed genes in the 3L vs 2L, P vs 2L, and P vs 3L comparisons were identified. The cluster co-expression patterns of the differentially expressed genes revealed that significantly differentially-expressed genes in the pupal stage were predicted to be related to diapause induction. All differentially expressed genes were investigated by GO functional and KEGG pathway analysis, and the results showed that genes involved in processes such as 20-hydroxyecdysone (20E) biosynthesis, cell cycle and metabolic pathways are likely related to obligatory diapause induction in B. minax. These results provide important information on the transcriptome of the Chinese citrus fruit fly that can be used for further functional studies as well as contributing to our understanding of the molecular basis of obligatory diapause induction and suggesting potential molecular targets for the control of this pest.

20-Hydroxyecdysone (20E) signaling as a promising target for the chemical control of malaria vectors

With the rapid development and spread of resistance to insecticides among anopheline malaria vectors, the efficacy of current World Health Organization (WHO)-approved insecticides targeting these vectors is under threat. This has led to the development of novel interventions, including improved and enhanced insecticide formulations with new targets or synergists or with added sterilants and/or antimalarials, among others. To date, several studies in mosquitoes have revealed that the 20-hydroxyecdysone (20E) signaling pathway regulates both vector abundance and competence, two parameters that influence malaria transmission. Therefore, insecticides which target 20E signaling (e.g. methoxyfenozide and halofenozide) may be an asset for malaria vector control. While such insecticides are already commercially available for lepidopteran and coleopteran pests, they still need to be approved by the WHO for malaria vector control programs. Until recently, chemicals targeting 20E signaling were considered to be insect growth regulators, and their effect was mostly studied against immature mosquito stages. However, in the last few years, promising results have been obtained by applying methoxyfenozide or halofenozide (two compounds that boost 20E signaling) to Anopheles populations at different phases of their life-cycle. In addition, preliminary studies suggest that methoxyfenozide resistance is unstable, causing the insects substantial fitness costs, thereby potentially circumventing one of the biggest challenges faced by current vector control efforts. In this review, we first describe the 20E signaling pathway in mosquitoes and then summarize the mechanisms whereby 20E signaling regulates the physiological processes associated with vector competence and vector abundance. Finally, we discuss the potential of using chemicals targeting 20E signaling to control malaria vectors.

Inhibition of Asaia in Adult Mosquitoes Causes Male-Specific Mortality and Diverse Transcriptome Changes

Mosquitoes can transmit many infectious diseases, such as malaria, dengue, Zika, yellow fever, and lymphatic filariasis. Current mosquito control strategies are failing to reduce the severity of outbreaks that still cause high human morbidity and mortality worldwide. Great expectations have been placed on genetic control methods. Among other methods, genetic modification of the bacteria colonizing different mosquito species and expressing anti-pathogen molecules may represent an innovative tool to combat mosquito-borne diseases. Nevertheless, this emerging approach, known as paratransgenesis, requires a detailed understanding of the mosquito microbiota and an accurate characterization of selected bacteria candidates. The acetic acid bacteria Asaia is a promising candidate for paratransgenic approaches. We have previously reported that Asaia symbionts play a beneficial role in the normal development of Anopheles mosquito larvae, but no study has yet investigated the role(s) of Asaia in adult mosquito biology. Here we report evidence on how treatment with a highly specific anti-Asaia monoclonal antibody impacts the survival and physiology of adult Anopheles stephensi mosquitoes. Our findings offer useful insight on the role of Asaia in several physiological systems of adult mosquitoes, where the influence differs between males and females.

Regulation of Hormone-Related Genes in Ericerus pela (Hemiptera: Coccidae) for Dimorphic Metamorphosis

Insect hormones regulate metamorphosis including that leading to sexual dimorphism. Using RNA-Seq, we discovered that the second-instar male larva (SM) of the white wax insect, Ericerus pela, have 5,968 and 8,620 differentially expressed transcripts compared with the second-instar female larva (SF) and the first-instar male larva (FM), respectively. The expression levels of genes involved in the apoptosis of old tissues and the reconstruction of new ones in the SM significantly enhanced, while the SF mainly has enhanced expression levels of anabolic genes such as chitin. We predicted that the second-instar larvae are the developmental origin of sexual dimorphic metamorphosis. Meanwhile, in the juvenile hormone (JH) metabolic pathway, CYP15A1 and JH esterase (JHE) are differentially expressed; and in the 20-hydroxyecdysone (20E) metabolic pathway, CYP307A1, CYP314A1, and CYP18A1 are differentially expressed. In the SM, the expression levels of CYP307A1 and CYP314A1 are significantly increased, whereas the expression level of CYP18A1 is significantly decreased; in the SF, the expression levels of the above genes are opposite to that of the SM. Expression trends of RNA-seq is consistent with the expression level of qRT-PCR, and seven of them are highly correlated (R ≥ 0.610) and four are moderately correlated (0.588 ≥ R ≥ 0.542).

Dual role of the Anopheles coluzzii Venus Kinase Receptor in both larval growth and immunity

Vector-borne diseases and especially malaria are responsible for more than half million deaths annually. The increase of insecticide resistance in wild populations of Anopheles malaria vectors emphasises the need for novel vector control strategies as well as for identifying novel vector targets. Venus kinase receptors (VKRs) constitute a Receptor Tyrosine Kinase (RTK) family only found in invertebrates. In this study we functionally characterized Anopheles VKR in the Gambiae complex member, Anophelescoluzzii. Results showed that Anopheles VKR can be activated by L-amino acids, with L-arginine as the most potent agonist. VKR was not required for the fecundity of A. coluzzii, in contrast to reports from other insects, but VKR function is required in both Anopheles males and females for development of larval progeny. Anopheles VKR function is also required for protection against infection by Plasmodium parasites, thus identifying a novel linkage between reproduction and immunity in Anopheles. The insect specificity of VKRs as well as the essential function for reproduction and immunity suggest that Anopheles VKR could be a potentially druggable target for novel vector control strategies.

Evolution of gene expression levels in the male reproductive organs of Anopheles mosquitoes

Modifications in gene expression determine many of the phenotypic differentiations between closely related species. This is particularly evident in reproductive tissues, where evolution of genes is more rapid, facilitating the appearance of distinct reproductive characteristics which may lead to species isolation and phenotypic variation. Large-scale, comparative analyses of transcript expression levels have been limited until recently by lack of inter-species data mining solutions. Here, by combining expression normalisation across lineages, multivariate statistical analysis, evolutionary rate, and protein-protein interaction analysis, we investigate ortholog transcripts in the male accessory glands and testes across five closely related species in the Anopheles gambiae complex. We first demonstrate that the differentiation by transcript expression is consistent with the known Anopheles phylogeny. Then, through clustering, we discover groups of transcripts with tissue-dependent expression patterns conserved across lineages, or lineage-dependent patterns conserved across tissues. The strongest associations with reproductive function, transcriptional regulatory networks, protein-protein subnetworks, and evolutionary rate are found for the groups of transcripts featuring large expression differences in lineage or tissue-conserved patterns.

Relevance of estrogen-related receptor gene and ecdysone receptor gene in adult testis of the cricket Teleogryllus emma (Orthoptera: Gryllidae)

Estrogen-related receptor gene (ERR) and ecdysone receptor gene (EcR) belong to the nuclear receptor gene superfamily, both of which are associated with the regulation of insect reproductive development. However, the relationship between ERR and EcR and whether ERR participates in the 20E signal pathway during male reproduction are unclear. In this paper, adult male crickets Teleogryllus emma Ohmschi & Matsumura were divided into the experimental group, negative group, and control group. Crickets of the experimental group were injected with TeERR or TeEcR-dsRNA, and those in the negative group received EGFP-dsRNA. The efficiency of TeERR and TeEcR-RNAi was detected in the experimental group. Furthermore, the transcription level, morphological characteristics as well as weight were analyzed in the TeERR or TeEcR knocked-down testis. Results showed that the expression level of TeERR or TeEcR was significantly down-regulated (P < 0.05) when treated with 2000 ng TeERR or TeEcR-dsRNA for 48 h. The expression level of TeERR could be down-regulated (P < 0.05) using TeEcR-RNAi and vice versa. TeERR and TeEcR-RNAi caused morphological changes in testes, but they had no obvious effect on weight (P > 0.05). These results indicate that TeERR and TeEcR are intimately related to each other. In addition, TeERR may be involved in the 20E signal pathway and maintain the function of adult cricket testis.

RNA interference in mosquito: understanding immune responses, double-stranded RNA delivery systems and potential applications in vector control

RNA interference (RNAi) refers to the process of post-transcriptional silencing of cellular mRNA by the application of double-stranded RNA (dsRNA). RNAi strategies have been widely employed to regulate gene expression in plants and animals including insects. With the availability of the full genome sequences of major vector mosquitoes, RNAi has been increasingly used to conduct genetic studies of human pathogens in mosquito vectors and to study the evolution of insecticide resistance in mosquitoes. This review summarizes the recent progress in our understanding of mosquito-pathogen interactions using RNAi and various methods of dsRNA delivery in mosquitoes at different stages. We also discuss potential applications of this technology to develop novel tools for vector control.

The Anopheles gambiae transcriptome - a turning point for malaria control

Mosquitoes are important vectors of several pathogens and thereby contribute to the spread of diseases, with social, economic and public health impacts. Amongst the approximately 450 species of Anopheles, about 60 are recognized as vectors of human malaria, the most important parasitic disease. In Africa, Anopheles gambiae is the main malaria vector mosquito. Current malaria control strategies are largely focused on drugs and vector control measures such as insecticides and bed-nets. Improvement of current, and the development of new, mosquito-targeted malaria control methods rely on a better understanding of mosquito vector biology. An organism's transcriptome is a reflection of its physiological state and transcriptomic analyses of different conditions that are relevant to mosquito vector competence can therefore yield important information. Transcriptomic analyses have contributed significant information on processes such as blood-feeding parasite-vector interaction, insecticide resistance, and tissue- and stage-specific gene regulation, thereby facilitating the path towards the development of new malaria control methods. Here, we discuss the main applications of transcriptomic analyses in An. gambiae that have led to a better understanding of mosquito vector competence.

Transcriptomic, cellular and life-history responses of Daphnia magna chronically exposed to benzotriazoles: Endocrine-disrupting potential and molting effects

Benzotriazoles (BZTs) are ubiquitous aquatic contaminants used in a wide range of industrial and domestic applications from aircraft deicers to dishwasher tablets. Acute toxicity has been reported in aquatic organisms for some of the BZTs but their mode of action remains unknown. The objectives of this study were to evaluate the transcriptomic response of D. magna exposed to sublethal doses of 1H-benzotriazole (BTR), 5-methyl-1H-benzotriazole (5MeBTR) and 5-chloro-1H-benzotriazole (5ClBTR) using RNA-sequencing and quantitative real-time PCR. Cellular and life-history endpoints (survival, number of neonates, growth) were also investigated. Significant effects on the molting frequency were observed after 21-d exposure to 5MeBTR and 5ClBTR. No effects on molting frequency were observed for BTR but RNA-seq results indicated that this BZT induced the up-regulation of genes coding for cuticular proteins, which could have compensated the molting disruption. Molting in cladocerans is actively controlled by ecdysteroid hormones. Complementary short-term temporal analysis (4- and 8-d exposure) of the transcription of genes related to molting and hormone-mediated processes indicated that the three compounds had specific modes of action. BTR induced the transcription of genes involved in 20-hydroxyecdysone synthesis, which suggests pro-ecdysteroid properties. 5ClBTR exposure induced protein activity and transcriptional levels of chitinase enzymes, associated with an impact on ecdysteroid signaling pathways, which could explain the decrease in molt frequency. Finally, 5MeBTR seemed to increase molt frequency through epigenetic processes. Overall, results suggested that molting effects observed at the physiological level could be linked to endocrine regulation impacts of BZTs at the molecular level.

The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate 'genomic island of speciation'

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a 'genomic island of divergence' spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X-island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X-island despite a lack of autosomal divergence; (ii) individuals with multiple-recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination-driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.

RNA interference-mediated knockdown of the Halloween gene Spookiest (CYP307B1) impedes adult eclosion in the western tarnished plant bug, Lygus hesperus

Ecdysteroids play a critical role in coordinating insect growth, development and reproduction. A suite of cytochrome P450 monooxygenases coded by what are collectively termed Halloween genes mediate ecdysteroid biosynthesis. In this study, we describe cloning and RNA interference (RNAi)-mediated knockdown of the CYP307B1 Halloween gene (Spookiest) in the western tarnished plant bug, Lygus hesperus. Transcripts for Ly. hesperus Spookiest (LhSpot) were amplified from all life stages and correlated well with timing of the pre-moult ecdysteroid pulse. In adults, LhSpot was amplified from heads of both genders as well as female reproductive tissues. Heterologous expression of a LhSpot fluorescent chimera in cultured insect cells co-localized with a fluorescent marker of the endoplasmic reticulum/secretory pathway. RNAi-mediated knockdown of LhSpot in fifth instars reduced expression of ecdysone-responsive genes E74 and E75, and prevented adult development. This developmental defect was rescued following application of exogenous 20-hydroxyecdysone but not exogenous 7-dehydrocholesterol. The unequivocal RNAi effects on Ly. hesperus development and the phenotypic rescue by 20-hydroxyecdysone are causal proof of the involvement of LhSpot in ecdysteroid biosynthesis and related developmental processes, and may provide an avenue for development of new control measures against Ly. hesperus.

Calcium influx enhances neuropeptide activation of ecdysteroid hormone production by mosquito ovaries

A critical step in mosquito reproduction is the ingestion of a blood meal from a vertebrate host. In mosquitoes like Aedes aegypti, blood feeding stimulates the release of ovary ecdysteroidogenic hormone (OEH) and insulin-like peptide 3 (ILP3). This induces the ovaries to produce ecdysteroid hormone (ECD), which then drives egg maturation. In many immature insects, prothoracicotropic hormone (PTTH) stimulates the prothoracic glands to produce ECD that directs molting and metamorphosis. The receptors for OEH, ILP3 and PTTH are different receptor tyrosine kinases with OEH and ILP3 signaling converging downstream in the insulin pathway and PTTH activating the mitogen-activated protein kinase pathway. Calcium (Ca(2+)) flux and cAMP have also been implicated in PTTH signaling, but the role of Ca(2+) in OEH, ILP3, and cAMP signaling in ovaries is unknown. Here, we assessed whether Ca(2+) flux affects OEH, ILP3, and cAMP activity in A. aegypti ovaries and also asked whether PTTH stimulated ovaries to produce ECD. Results indicated that Ca(2+) flux enhanced but was not essential for OEH or ILP3 activity, whereas cAMP signaling was dependent on Ca(2+) flux. Recombinant PTTH from Bombyx mori fully activated ECD production by B. mori PTGs, but exhibited no activity toward A. aegypti ovaries. Recombinant PTTH from A. aegypti also failed to stimulate either B. mori PTGs or A. aegypti ovaries to produce ECD. We discuss the implications of these results in the context of mosquito reproduction and ECD biosynthesis by insects generally.

Early transcriptional events linked to induction of diapause revealed by RNAseq in larvae of drosophilid fly, Chymomyza costata

Background

Diapause is a developmental alternative to direct ontogeny in many invertebrates. Its primary adaptive meaning is to secure survival over unfavourable seasons in a state of developmental arrest usually accompanied by metabolic suppression and enhanced tolerance to environmental stressors. During photoperiodically triggered diapause of insects, the ontogeny is centrally turned off under hormonal control, the molecular details of this transition being poorly understood. Using RNAseq technology, we characterized transcription profiles associated with photoperiodic diapause induction in the larvae of the drosophilid fly Chymomyza costata with the goal of identifying candidate genes and processes linked to upstream regulatory events that eventually lead to a complex phenotypic change.

Results

Short day photoperiod triggering diapause was associated to inhibition of 20-hydroxy ecdysone (20-HE) signalling during the photoperiod-sensitive stage of C. costata larval development. The mRNA levels of several key genes involved in 20-HE biosynthesis, perception, and signalling were significantly downregulated under short days. Hormonal change was translated into downregulation of a series of other transcripts with broad influence on gene expression, protein translation, alternative histone marking by methylation and alternative splicing. These changes probably resulted in blockade of direct development and deep restructuring of metabolic pathways indicated by differential expression of genes involved in cell cycle regulation, metabolism, detoxification, redox balance, protection against oxidative stress, cuticle formation and synthesis of larval storage proteins. This highly complex alteration of gene transcription was expressed already during first extended night, within the first four hours after the change of the photoperiodic signal from long days to short days. We validated our RNAseq differential gene expression results in an independent qRT-PCR experiment involving wild-type (photoperiodic) and NPD-mutant (non-photoperiodic) strains of C. costata.

Conclusions

Our study revealed several strong candidate genes for follow-up functional studies. Candidate genes code for upstream regulators of a complex change of gene expression, which leads to phenotypic switch from direct ontogeny to larval diapause.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1907-4) contains supplementary material, which is available to authorized users.

Redefining metamorphosis in spiny lobsters: molecular analysis of the phyllosoma to puerulus transition in Sagmariasus verreauxi

The molecular understanding of crustacean metamorphosis is hindered by small sized individuals and inability to accurately define molt stages. We used the spiny lobster Sagmariasus verreauxi where the large, transparent larvae enable accurate tracing of the transition from a leaf-shaped phyllosoma to an intermediate larval-juvenile phase (puerulus). Transcriptomic analysis of larvae at well-defined stages prior to, during, and following this transition show that the phyllosoma-puerulus metamorphic transition is accompanied by vast transcriptomic changes exceeding 25% of the transcriptome. Notably, genes previously identified as regulating metamorphosis in other crustaceans do not fluctuate during this transition but in the later, morphologically-subtle puerulus-juvenile transition, indicating that the dramatic phyllosoma-puerulus morphological shift relies on a different, yet to be identified metamorphic mechanism. We examined the change in expression of domains and gene families, with focus on several key genes. Our research implies that the separation in molecular triggering systems between the phyllosoma-puerulus and puerulus-juvenile transitions might have enabled the extension of the oceanic phase in spiny lobsters. Study of similar transitions, where metamorphosis is uncoupled from the transition into the benthic juvenile form, in other commercially important crustacean groups might show common features to point on the evolutionary advantage of this two staged regulation.

Rice stripe virus counters reduced fecundity in its insect vector by modifying insect physiology, primary endosymbionts and feeding behavior

Virus-vector relationships can be complex and diverse as a result of long-term coevolution. Understanding these interactions is crucial for disease and vector management. Rice stripe virus (RSV) is known to be transovarially transmitted within its vector, Laodelphax striatellus, and causes serious rice stripe disease. In RSV-infected L. striatellus, we found contrasting changes in vector fecundity, physiology, primary endosymbionts (i.e. yeast-like symbionts, YLS) and feeding behavior that can interact to affect the spread of RSV. RSV-infected L. striatellus exhibited a significant decrease in fecundity that could lead a reduction of viruliferous individuals in populations. As a potential response to this loss, RSV infection also significantly shortened nymphal stage duration, which can strengthen RSV vertical circulation in L. striatellus populations and promote RSV spreading by adult migration and dispersal. Down-regulated JHAMT and up-regulated CYP307A1 in the juvenile hormone and ecdysteroid pathways, respectively, were linked to accelerated development. RSV-infected adults were also found to have higher body weight in conjunction with increased YLS abundance. Furthermore, prolonged host plant phloem exposure to salivation by RSV-infected adults should further enhance RSV horizontal transmission. Our study highlights potential strategies of RSV in enhancing its transmission, and provides new insights into the complexity of virus-vector interactions.