Genome-wide analysis of gene expression in adult Anopheles gambiae

A potential link between aromatics-induced oviposition repellency behaviors and specific odorant receptor of Aedes albopictus

BACKGROUND

The Asian tiger mosquito, Aedes albopictus, is a competent vector for the spread of several viral arboviruses including dengue, chikungunya, and Zika. Several vital mosquito behaviors linked to survival and reproduction are primarily dependent on a sophisticated olfactory system for semiochemical perception. However, a limited number of studies has hampered our understanding of the relationship between the A. albopictus acute olfactory system and the complex chemical world.

RESULTS

Here, we performed a qRT-PCR assay on antennae from A. albopictus of differing sex, age and physiological states, and found that AalbOr10 was enriched in blood-fed female mosquitoes. We then undertook single sensillum recording to de-orphan AalbOr10 using a panel of physiologically and behaviorally relevant odorants in a Drosophila 'empty neuron' system. The results indicated that AalbOr10 was activated by seven aromatic compounds, all of which hampered egg-laying in blood-fed mosquitoes. Furthermore, using a post-RNA interference oviposition assay, we found that reducing the transcript level of AalbOr10 affected repellent activity mediated by 2-ethylphenol at low concentrations (10-4 vol/vol). Computational modeling and molecular docking studies suggested that hydrogen bonds to Y68 and Y150 mediated the interaction of 2-ethylphenol with AalbOr10.

CONCLUSION

We reveal a potential link between aromatics-induced oviposition repellency behaviors and a specific odorant receptor in A. albopictus. Our findings provide a foundation for identifying active semiochemicals for the monitoring or controlling of mosquito populations. © 2024 Society of Chemical Industry.

Development of circulating isolates of Plasmodium falciparum is accelerated in Anopheles vectors with reduced reproductive output

Anopheles gambiae and its sibling species Anopheles coluzzii are the most efficient vectors of the malaria parasite Plasmodium falciparum. When females of these species feed on an infected human host, oogenesis and parasite development proceed concurrently, but interactions between these processes are not fully understood. Using multiple natural P. falciparum isolates from Burkina Faso, we show that in both vectors, impairing steroid hormone signaling to disrupt oogenesis leads to accelerated oocyst growth and in a manner that appears to depend on both parasite and mosquito genotype. Consistently, we find that egg numbers are negatively linked to oocyst size, a metric for the rate of oocyst development. Oocyst growth rates are also strongly accelerated in females that are in a pre-gravid state, i.e. that fail to develop eggs after an initial blood meal. Overall, these findings advance our understanding of mosquito-parasite interactions that influence P. falciparum development in malaria-endemic regions.

Age determination of Chrysomya megacephala (Diptera: Calliphoridae) using lifespan patterns, gene expression, and pteridine concentration under constant and variable temperatures

Chrysomya megacephala (Fabricius, 1794) (Diptera: Calliphoridae), is a blowfly species widely studied in medical, veterinary, and entomological research. Our study examined the impact of constant (15, 20, 25, 30, and 35 °C) and variable (ranging from 21.0 to 25.4 °C, with an average of 23.31 °C) temperatures on the development and larval body length of C. megacephala. Additionally, we analyzed the age of the adult C. megacephala through pteridine content and related metabolic genes analysis. Our findings revealed three distinct growth patterns: isomorphen diagram, isomegalen diagram, and thermal accumulated models. At constant temperatures of 15, 20, 25, 30, and 35 °C, egg-hatching times were 44.5 ± 8.9, 26.7 ± 4.6, 12.6 ± 1.1, 11.0 ± 1.0, and 9.9 ± 1.9 h, respectively, while it was 15.3 ± 5.9 h at variable temperatures. The total development times from oviposition to adult eclosion in C. megacephala required 858.1 ± 69.2, 362.3 ± 5.9, 289.6 ± 17.8, 207.3 ± 9.3, and 184.7 ± 12.1 h at constant temperatures of 15, 20, 25, 30, and 35 °C, respectively. This duration was extended to 282.0 ± 64.1 h under variable temperatures. However, no significant differences were found in hatching times and the total developmental durations between 25 °C and variable temperatures. A developmental threshold temperature (D0) of 9.90 ± 0.77 °C and a thermal summation constant (K) of 4244.0 ± 347.0° hours were ascertained. Pteridine content patterns varied significantly across constant temperatures, but not between 25 °C and variable temperatures. Sex and temperature were identified as the primary factors influencing pteridine levels in the head of C. megacephala. Gene expression associated with pteridine metabolism decreased following adult eclosion, matching with increased pteridine concentration. Further investigations are needed to explore the use of pteridine cofactors for age-grading adult necrophagous flies. These findings provide valuable insights into the lifespan of C. megacephala, thereby offering valuable groundwork for forthcoming investigations and PMImin determination.

GeCoNet-Tool: a software package for gene co-expression network construction and analysis

BACKGROUND

Network analysis is a powerful tool for studying gene regulation and identifying biological processes associated with gene function. However, constructing gene co-expression networks can be a challenging task, particularly when dealing with a large number of missing values.

RESULTS

We introduce GeCoNet-Tool, an integrated gene co-expression network construction and analysis tool. The tool comprises two main parts: network construction and network analysis. In the network construction part, GeCoNet-Tool offers users various options for processing gene co-expression data derived from diverse technologies. The output of the tool is an edge list with the option of weights associated with each link. In network analysis part, the user can produce a table that includes several network properties such as communities, cores, and centrality measures. With GeCoNet-Tool, users can explore and gain insights into the complex interactions between genes.

The structure of AgamOBP5 in complex with the natural insect repellents Carvacrol and Thymol: Crystallographic, fluorescence and thermodynamic binding studies

Among several proteins participating in the olfactory perception process of insects, Odorant Binding Proteins (OBPs) are today considered valid targets for the discovery of compounds that interfere with their host-detection behavior. The 3D structures of Anopheles gambiae mosquito AgamOBP1 in complex with the known synthetic repellents DEET and Icaridin have provided valuable information on the structural characteristics that govern their selective binding. However, no structure of a plant-derived repellent bound to an OBP has been available until now. Herein, we present the novel three-dimensional crystal structures of AgamOBP5 in complex with two natural phenolic monoterpenoid repellents, Carvacrol and Thymol, and the MPD molecule. Structural analysis revealed that both monoterpenoids occupy a binding site (Site-1) by adopting two alternative conformations. An additional Carvacrol was also bound to a secondary site (Site-2) near the central cavity entrance. A protein-ligand hydrogen-bond network supplemented by van der Waals interactions spans the entire binding cavity, bridging α4, α6, and α3 helices and stabilizing the overall structure. Fluorescence competition and Differential Scanning Calorimetry experiments verified the presence of two binding sites and the stabilization effect on AgamOBP5. While Carvacrol and Thymol bind to Site-1 with equal affinity in the submicromolar range, they exhibit a significantly lower and distinct binding capacity for Site-2 with Kd's of ~7 μΜ and ~18 μΜ, respectively. Finally, a comparison of AgamOBP5 complexes with the AgamOBP4-Indole structure revealed that variations of ligand-interacting aminoacids such as A109T, I72M, A112L, and A105T cause two structurally similar and homologous proteins to display different binding specificities.

A male steroid controls female sexual behaviour in the malaria mosquito

Insects, unlike vertebrates, are widely believed to lack male-biased sex steroid hormones¹. In the malaria mosquito Anopheles gambiae, the ecdysteroid 20-hydroxyecdysone (20E) appears to have evolved to both control egg development when synthesized by females² and to induce mating refractoriness when sexually transferred by males³. Because egg development and mating are essential reproductive traits, understanding how Anopheles females integrate these hormonal signals can spur the design of new malaria control programs. Here we reveal that these reproductive functions are regulated by distinct sex steroids through a sophisticated network of ecdysteroid-activating/inactivating enzymes. We identify a male-specific oxidized ecdysteroid, 3-dehydro-20E (3D20E), which safeguards paternity by turning off female sexual receptivity following its sexual transfer and activation by dephosphorylation. Notably, 3D20E transfer also induces expression of a reproductive gene that preserves egg development during Plasmodium infection, ensuring fitness of infected females. Female-derived 20E does not trigger sexual refractoriness but instead licenses oviposition in mated individuals once a 20E-inhibiting kinase is repressed. Identifying this male-specific insect steroid hormone and its roles in regulating female sexual receptivity, fertility and interactions with Plasmodium parasites suggests the possibility for reducing the reproductive success of malaria-transmitting mosquitoes.

A global [Formula: see text] gene co-expression network constructed from hundreds of experimental conditions with missing values

BACKGROUND

Gene co-expression networks (GCNs) can be used to determine gene regulation and attribute gene function to biological processes. Different high throughput technologies, including one and two-channel microarrays and RNA-sequencing, allow evaluating thousands of gene expression data simultaneously, but these methodologies provide results that cannot be directly compared. Thus, it is complex to analyze co-expression relations between genes, especially when there are missing values arising for experimental reasons. Networks are a helpful tool for studying gene co-expression, where nodes represent genes and edges represent co-expression of pairs of genes.

RESULTS

In this paper, we establish a method for constructing a gene co-expression network for the Anopheles gambiae transcriptome from 257 unique studies obtained with different methodologies and experimental designs. We introduce the sliding threshold approach to select node pairs with high Pearson correlation coefficients. The resulting network, which we name AgGCN1.0, is robust to random removal of conditions and has similar characteristics to small-world and scale-free networks. Analysis of network sub-graphs revealed that the core is largely comprised of genes that encode components of the mitochondrial respiratory chain and the ribosome, while different communities are enriched for genes involved in distinct biological processes.

CONCLUSION

Analysis of the network reveals that both the architecture of the core sub-network and the network communities are based on gene function, supporting the power of the proposed method for GCN construction. Application of network science methodology reveals that the overall network structure is driven to maximize the integration of essential cellular functions, possibly allowing the flexibility to add novel functions.

Beyond the eye: Kynurenine pathway impairment causes midgut homeostasis dysfunction and survival and reproductive costs in blood-feeding mosquitoes

Insect ommochrome biosynthesis pathways metabolize tryptophan to generate eye-color pigments and wild-type alleles of pathway genes are useful phenotypic markers in transgenesis studies. Pleiotropic effects of mutations in some genes exert a load on both survival and reproductive success in blood-feeding species. Here, we investigated the challenges imposed on mosquitoes by the increase of tryptophan metabolites resulting from blood meal digestion and the impact of disruptions of the ommochrome biosynthesis pathway. Female mosquitoes with spontaneous and induced mutations in the orthologs of the genes encoding kynurenine hydroxylase in Aedes aegypti, Anopheles stephensi and Culex quinquefasciatus exhibited impaired survival and reproductive phenotypes that varied in type and severity among the species. A compromised midgut permeability barrier function was also observed in An. stephensi. Surprisingly, mutant mosquitoes displayed an increase in microbiota compared to controls that was not accompanied by a general induction of immune genes. Antibiotic treatment rescued some deleterious traits implicating a role for the kynurenine pathway (KP) in midgut homeostasis. Supplemental xanthurenic acid, a KP end-product, rescued lethality and limited microbiota proliferation in Ae. aegypti. These data implicate the KP in the regulation of the host/microbiota interface. These pleiotropic effects on mosquito physiology are important in the development of genetic strategies targeting vector mosquitoes.

Functional constraints on insect immune system components govern their evolutionary trajectories

Roles of constraints in shaping evolutionary outcomes are often considered in the contexts of developmental biology and population genetics, in terms of capacities to generate new variants and how selection limits or promotes consequent phenotypic changes. Comparative genomics also recognizes the role of constraints, in terms of shaping evolution of gene and genome architectures, sequence evolutionary rates, and gene gains or losses, as well as on molecular phenotypes. Characterizing patterns of genomic change where putative functions and interactions of system components are relatively well described offers opportunities to explore whether genes with similar roles exhibit similar evolutionary trajectories. Using insect immunity as our test case system, we hypothesize that characterizing gene evolutionary histories can define distinct dynamics associated with different functional roles. We develop metrics that quantify gene evolutionary histories, employ these to characterize evolutionary features of immune gene repertoires, and explore relationships between gene family evolutionary profiles and their roles in immunity to understand how different constraints may relate to distinct dynamics. We identified three main axes of evolutionary trajectories characterized by gene duplication and synteny, maintenance/stability and sequence conservation, and loss and sequence divergence, highlighting similar and contrasting patterns across these axes amongst subsets of immune genes. Our results suggest that where and how genes participate in immune responses limit the range of possible evolutionary scenarios they exhibit. The test case study system of insect immunity highlights the potential of applying comparative genomics approaches to characterize how functional constraints on different components of biological systems govern their evolutionary trajectories.

Resistance to a CRISPR-based gene drive at an evolutionarily conserved site is revealed by mimicking genotype fixation

CRISPR-based homing gene drives can be designed to disrupt essential genes whilst biasing their own inheritance, leading to suppression of mosquito populations in the laboratory. This class of gene drives relies on CRISPR-Cas9 cleavage of a target sequence and copying ('homing') therein of the gene drive element from the homologous chromosome. However, target site mutations that are resistant to cleavage yet maintain the function of the essential gene are expected to be strongly selected for. Targeting functionally constrained regions where mutations are not easily tolerated should lower the probability of resistance. Evolutionary conservation at the sequence level is often a reliable indicator of functional constraint, though the actual level of underlying constraint between one conserved sequence and another can vary widely. Here we generated a novel adult lethal gene drive (ALGD) in the malaria vector Anopheles gambiae, targeting an ultra-conserved target site in a haplosufficient essential gene (AGAP029113) required during mosquito development, which fulfils many of the criteria for the target of a population suppression gene drive. We then designed a selection regime to experimentally assess the likelihood of generation and subsequent selection of gene drive resistant mutations at its target site. We simulated, in a caged population, a scenario where the gene drive was approaching fixation, where selection for resistance is expected to be strongest. Continuous sampling of the target locus revealed that a single, restorative, in-frame nucleotide substitution was selected. Our findings show that ultra-conservation alone need not be predictive of a site that is refractory to target site resistance. Our strategy to evaluate resistance in vivo could help to validate candidate gene drive targets for their resilience to resistance and help to improve predictions of the invasion dynamics of gene drives in field populations.

Use of transcriptional age grading technique to determine the chronological age of Sri Lankan Aedes aegypti and Aedes albopictus females

BACKGROUND

Aedes aegypti and Ae. albopictus are important vectors of human diseases such as dengue, chikungunya, and zika. In Sri Lanka, they have been responsible for transmitting dengue virus. One of the most important parameters influencing the likelihood of arbovirus transmission is the age structure of the mosquito population. However, mosquito age is difficult to measure with accuracy. This study aims to construct multivariate calibration models using the transcriptional abundance of three age-responsive genes: Ae15848 (calcium-binding protein), Ae8505 (structural component of cuticle), and Ae4274 (fizzy cell cycle/cell division cycle 20).

METHODS

The transcriptional age-grading technique was applied to determine the chronological age of Ae. aegypti and Ae. albopictus female mosquito populations from Sri Lanka using the age-responsive genes Ae15848, Ae8505, and Ae4274. Furthermore, Ae. aegypti samples obtained from colonies reared at two temperatures (23 and 27 °C) were used to investigate the influence of temperature on this age-grading technique. Expression levels of these three genes were quantified using reverse transcription qualitative PCR (qRT-PCR), and results were normalized against the housekeeping gene ribosomal gene S17 (RpS17).

RESULTS

The expression of Ae15848 and Ae8505 decreased with the age of mosquitoes and showed the most significant and consistent change while expression of Ae4274 increased with age. The multivariate calibration models showed > 80% correlation between expression of these age-responsive genes and the age of female mosquitoes at both temperatures. At 27 °C the accuracy of age predictions using the models was 2.19 (± 1.66) days and 2.58 (± 2.06) days for Ae. aegypti and Ae. albopictus females, respectively. The accuracy of the model for Ae. aegypti at 23 °C was 3.42 (± 2.74) days.

CONCLUSIONS

An adult rearing temperature difference of 4 °C (23-27 °C) did not significantly affect the age predictions. The calibration models created during this study could be successfully used to estimate the age of wild Ae. aegypti and Ae. albopictus mosquitoes from Sri Lanka.

Review and Meta-Analysis of the Evidence for Choosing between Specific Pyrethroids for Programmatic Purposes

Simple Summary

A group of insecticides, called pyrethroids, are the main strategy for controlling the mosquito vectors of malaria. Pyrethroids are used in all insecticide-treated bednets, and many indoor residual spray programmes (in which insecticides are sprayed on the interior walls of houses). There are different types of pyrethroids within the class (e.g., deltamethrin and permethrin). Across the world, mosquitoes are showing signs of resistance to the pyrethroids, such as reduced mortality following contact. However, it is unclear if this resistance is uniform across the pyrethroid class (i.e., if a mosquito is resistant to deltamethrin, whether it is resistant to permethrin at the same level). In addition, it is not known if switching between different pyrethroids can be used to effectively maintain mosquito control when resistance to a single pyrethroid has been detected. This review examined the evidence from molecular studies, resistance testing from laboratory and field data, and mosquito behavioural assays to answer these questions. The evidence suggested that in areas where pyrethroid resistance exists, different mortality seen between the pyrethroids is not necessarily indicative of an operationally relevant difference in control performance, and there is no reason to rotate between common pyrethroids (i.e., deltamethrin, permethrin, and alpha-cypermethrin) as an insecticide resistance management strategy.

Abstract

Pyrethroid resistance is widespread in malaria vectors. However, differential mortality in discriminating dose assays to different pyrethroids is often observed in wild populations. When this occurs, it is unclear if this differential mortality should be interpreted as an indication of differential levels of susceptibility within the pyrethroid class, and if so, if countries should consider selecting one specific pyrethroid for programmatic use over another. A review of evidence from molecular studies, resistance testing with laboratory colonies and wild populations, and mosquito behavioural assays were conducted to answer these questions. Evidence suggested that in areas where pyrethroid resistance exists, different results in insecticide susceptibility assays with specific pyrethroids currently in common use (deltamethrin, permethrin, α-cypermethrin, and λ-cyhalothrin) are not necessarily indicative of an operationally relevant difference in potential performance. Consequently, it is not advisable to use rotation between these pyrethroids as an insecticide-resistance management strategy. Less commonly used pyrethroids (bifenthrin and etofenprox) may have sufficiently different modes of action, though further work is needed to examine how this may apply to insecticide resistance management.

Unravelling population structure heterogeneity within the genome of the malaria vector Anopheles gambiae

Background

Whole genome re-sequencing provides powerful data for population genomic studies, allowing robust inferences of population structure, gene flow and evolutionary history. For the major malaria vector in Africa, Anopheles gambiae, other genetic aspects such as selection and adaptation are also important. In the present study, we explore population genetic variation from genome-wide sequencing of 765 An. gambiae and An. coluzzii specimens collected from across Africa. We used t-SNE, a recently popularized dimensionality reduction method, to create a 2D-map of An. gambiae and An. coluzzii genes that reflect their population structure similarities.

Results

The map allows intuitive navigation among genes distributed throughout the so-called “mainland” and numerous surrounding “island-like” gene clusters. These gene clusters of various sizes correspond predominantly to low recombination genomic regions such as inversions and centromeres, and also to recent selective sweeps. Because this mosquito species complex has been studied extensively, we were able to support our interpretations with previously published findings. Several novel observations and hypotheses are also made, including selective sweeps and a multi-locus selection event in Guinea-Bissau, a known intense hybridization zone between An. gambiae and An. coluzzii.

Conclusions

Our results present a rich dataset that could be utilized in functional investigations aiming to shed light onto An. gambiae s.l genome evolution and eventual speciation. In addition, the methodology presented here can be used to further characterize other species not so well studied as An. gambiae, shortening the time required to progress from field sampling to the identification of genes and genomic regions under unique evolutionary processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07722-y.

RNA-Seq analysis of blood meal induced gene-expression changes in Aedes aegypti ovaries

Background

Transmission of pathogens by vector mosquitoes is intrinsically linked with mosquito’s reproductive strategy because anautogenous mosquitoes require vertebrate blood to develop a batch of eggs. Each cycle of egg maturation is tightly linked with the intake of a fresh blood meal for most species. Mosquitoes that acquire pathogens during the first blood feeding can transmit the pathogens to susceptible hosts during subsequent blood feeding and also vertically to the next generation via infected eggs. Large-scale gene-expression changes occur following each blood meal in various tissues, including ovaries. Here we analyzed mosquito ovary transcriptome following a blood meal at three different time points to investigate blood-meal induced changes in gene expression in mosquito ovaries.

Results

We collected ovaries from Aedes aegypti that received a sugar meal or a blood meal on days 3, 10 and 20 post blood meal for transcriptome analysis. Over 4000 genes responded differentially following ingestion of a blood meal on day 3, and 660 and 780 genes on days 10 and 20, respectively. Proteins encoded by differentially expressed genes (DEGs) on day 3 include odorant binding proteins (OBPs), defense-specific proteins, and cytochrome P450 detoxification enzymes. In addition, we identified 580 long non-coding RNAs that are differentially expressed at three time points. Gene ontology analysis indicated that genes involved in peptidase activity, oxidoreductase activity, extracellular space, and hydrolase activity, among others were enriched on day 3. Although most of the DEGs returned to the nonsignificant level compared to the sugar-fed mosquito ovaries following oviposition on days 10 and 20, there remained differences in the gene expression pattern in sugar-fed and blood-fed mosquitoes.

Conclusions

Enrichment of OBPs following blood meal ingestion suggests that these genes may have other functions besides being part of the olfactory system. The enrichment of immune-specific genes and cytochrome P450 genes indicates that ovaries become well prepared to protect their germ line from any pathogens that may accompany the blood meal or from environmental contamination during oviposition, and to deal with the detrimental effects of toxic metabolites.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07551-z.

Transcriptomic analysis reveals pronounced changes in gene expression due to sub-lethal pyrethroid exposure and ageing in insecticide resistance Anopheles coluzzii

BACKGROUND

Malaria control is heavily reliant on the use of insecticides that target and kill the adult female Anopheline vector. The intensive use of insecticides of the pyrethroid class has led to widespread resistance in mosquito populations. The intensity of pyrethroid resistance in some settings in Africa means mosquitoes can contact bednets treated with this insecticide class multiple times with minimal mortality effects. Furthermore, both ageing and diel cycle have been shown to have large impacts on the resistance phenotype. Together, these traits may affect other aspects of vector biology controlling the vectorial capacity or fitness of the mosquito.

RESULTS

Here we show that sublethal exposure of a highly resistant Anopheles coluzzii population originally from Burkina Faso to the pyrethroid deltamethrin results in large and sustained changes to transcript expression. We identify five clear patterns in the data showing changes to transcripts relating to: DNA repair, respiration, translation, behaviour and oxioreductase processes. Further, we highlight differential regulation of transcripts from detoxification families previously linked with insecticide resistance, in addition to clear down-regulation of the oxidative phosphorylation pathway both indicative of changes in metabolism post-exposure. Finally, we show that both ageing and diel cycle have major effects on known insecticide resistance related transcripts.

CONCLUSION

Sub-lethal pyrethroid exposure, ageing and the diel cycle results in large-scale changes in the transcriptome of the major malaria vector Anopheles coluzzii. Our data strongly supports further phenotypic studies on how transcriptional changes such as reduced expression of the oxidative phosphorylation pathway or pyrethroid induced changes to redox state might impact key mosquito traits, such as vectorial capacity and life history traits.

Sex-Biased Gene Expression and Evolution in the Cerebrum and Syrinx of Chinese Hwamei (Garrulax canorus)

It is common that males and females display sexual dimorphisms, which usually result from sex-biased gene expression. Chinese hwamei (Garrulax canorus) is a good model for studying sex-biased gene expression because the song between the sexes is quite different. In this study, we analyze cerebrum and syrinx sex-biased gene expression and evolution using the de novo assembled Chinese hwamei transcriptome. In both the cerebrum and syrinx, our study revealed that most female-biased genes were actively expressed in females only, while most male-biased genes were actively expressed in both sexes. In addition, both male- and female-biased genes were enriched on the putative Z chromosome, suggesting the existence of sexually antagonistic genes and the insufficient dosage compensation of the Z-linked genes. We also identified a 9 Mb sex linkage region on the putative 4A chromosome which enriched more than 20% of female-biased genes. Resultantly, male-biased genes in both tissues had significantly higher Ka/Ks and effective number of codons (ENCs) than unbiased genes, and this suggested that male-biased genes which exhibit accelerated divergence may have resulted from positive selection. Taken together, our results initially revealed the reasons for the differences in singing behavior between males and females of Chinese hwamei.

Time Flies-Age Grading of Adult Flies for the Estimation of the Post-Mortem Interval

The estimation of the minimum time since death is one of the main applications of forensic entomology. This can be done by calculating the age of the immature stage of necrophagous flies developing on the corpse, which is confined to approximately 2–4 weeks, depending on temperature and species of the first colonizing wave of flies. Adding the age of the adult flies developed on the dead body could extend this time frame up to several weeks when the body is in a building or closed premise. However, the techniques for accurately estimating the age of adult flies are still in their beginning stages or not sufficiently validated. Here we review the current state of the art of analysing the aging of flies by evaluating the ovarian development, the amount of pteridine in the eyes, the degree of wing damage, the modification of their cuticular hydrocarbon patterns, and the increasing number of growth layers in the cuticula. New approaches, including the use of age specific molecular profiles based on the levels of gene and protein expression and the application of near infrared spectroscopy, are introduced, and the forensic relevance of these methods is discussed.

The Phlebotomus papatasi systemic transcriptional response to trypanosomatid-contaminated blood does not differ from the non-infected blood meal

Background

Leishmaniasis, caused by parasites of the genus Leishmania, is a disease that affects up to 8 million people worldwide. Parasites are transmitted to human and animal hosts through the bite of an infected sand fly. Novel strategies for disease control require a better understanding of the key step for transmission, namely the establishment of infection inside the fly.

Methods

The aim of this work was to identify sand fly systemic transcriptomic signatures associated with Leishmania infection. We used next generation sequencing to describe the transcriptome of whole Phlebotomus papatasi sand flies when fed with blood alone (control) or with blood containing one of three trypanosomatids: Leishmania major, L. donovani and Herpetomonas muscarum, the latter being a parasite not transmitted to humans.

Results

Of the trypanosomatids studied, only L. major was able to successfully establish an infection in the host P. papatasi. However, the transcriptional signatures observed after each parasite-contaminated blood meal were not specific to success or failure of a specific infection and they did not differ from each other. The transcriptional signatures were also indistinguishable after a non-contaminated blood meal.

Conclusions

The results imply that sand flies perceive Leishmania as just one feature of their microbiome landscape and that any strategy to tackle transmission should focus on the response towards the blood meal rather than parasite establishment. Alternatively, Leishmania could suppress host responses. These results will generate new thinking around the concept of stopping transmission by controlling the parasite inside the insect.

CRISPR/Cas9 knockout of female-biased genes AeAct-4 or myo-fem in Ae. aegypti results in a flightless phenotype in female, but not male mosquitoes

Aedes aegypti is a vector of dengue, chikungunya, and Zika viruses. Current vector control strategies such as community engagement, source reduction, and insecticides have not been sufficient to prevent viral outbreaks. Thus, interest in novel strategies involving genetic engineering is growing. Female mosquitoes rely on flight to mate with males and obtain a bloodmeal from a host. We hypothesized that knockout of genes specifically expressed in female mosquitoes associated with the indirect flight muscles would result in a flightless female mosquito. Using CRISPR-Cas9 we generated loss-of-function mutations in several genes hypothesized to control flight in mosquitoes, including actin (AeAct-4) and myosin (myo-fem) genes expressed specifically in the female flight muscle. Genetic knockout of these genes resulted in 100% flightless females, with homozygous males able to fly, mate, and produce offspring, albeit at a reduced rate when compared to wild type males. Interestingly, we found that while AeAct-4 was haplosufficient, with most heterozygous individuals capable of flight, this was not the case for myo-fem, where about half of individuals carrying only one intact copy could not fly. These findings lay the groundwork for developing novel mechanisms of controlling Ae. aegypti populations, and our results suggest that this mechanism could be applicable to other vector species of mosquito.

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.

Analysis of blood-induced Anopheles gambiae midgut proteins and sexual stage Plasmodium falciparum interaction reveals mosquito genes important for malaria transmission

Plasmodium invasion of mosquito midguts is a mandatory step for malaria transmission. The roles of mosquito midgut proteins and parasite interaction during malaria transmission are not clear. This study aims to identify mosquito midgut proteins that interact with and affect P. falciparum invasion. Based on gene expression profiles and protein sequences, 76 mosquito secretory proteins that are highly expressed in midguts and up-regulated by blood meals were chosen for analysis. About 61 candidate genes were successfully cloned from Anopheles gambiae and expressed in insect cells. ELISA analysis showed that 25 of the insect cell-expressed recombinant mosquito proteins interacted with the P. falciparum-infected cell lysates. Indirect immunofluorescence assays confirmed 17 of them interacted with sexual stage parasites significantly stronger than asexual stage parasites. Knockdown assays found that seven candidate genes significantly changed mosquitoes' susceptibility to P. falciparum. Four of them (AGAP006268, AGAP002848, AGAP006972, and AGAP002851) played a protective function against parasite invasion, and the other three (AGAP008138, FREP1, and HPX15) facilitated P. falciparum transmission to mosquitoes. Notably, AGAP008138 is a unique gene that only exists in Anopheline mosquitoes. These gene products are ideal targets to block malaria transmission.

Dynamic miRNA-mRNA interactions coordinate gene expression in adult Anopheles gambiae

microRNAs (miRNAs) are increasingly recognized as important regulators of many biological processes in mosquitoes, vectors of numerous devastating infectious diseases. Identification of bona fide targets remains the bottleneck for functional studies of miRNAs. In this study, we used CLEAR-CLIP assays to systematically analyze miRNA-mRNA interactions in adult female Anopheles gambiae mosquitoes. Thousands of miRNA-target pairs were captured after direct ligation of the miRNA and its cognate target transcript in endogenous Argonaute–miRNA–mRNA complexes. Using two interactions detected in this manner, miR-309-SIX4 and let-7-kr-h1, we demonstrated the reliability of this experimental approach in identifying in vivo gene regulation by miRNAs. The miRNA-mRNA interaction dataset provided an invaluable opportunity to decipher targeting rules of mosquito miRNAs. Enriched motifs in the diverse targets of each miRNA indicated that the majority of mosquito miRNAs rely on seed-based canonical target recognition, while noncanonical miRNA binding sites are widespread and often contain motifs complementary to the central or 3’ ends of miRNAs. The time-lapse study of miRNA-target interactomes in adult female mosquitoes revealed dynamic miRNA regulation of gene expression in response to varying nutritional sources and physiological demands. Interestingly, some miRNAs exhibited flexibility to use distinct sequences at different stages for target recognition. Furthermore, many miRNA-mRNA interactions displayed stage-specific patterns, especially for those genes involved in metabolism, suggesting that miRNAs play critical roles in precise control of gene expression to cope with enormous physiological demands associated with egg production. The global mapping of miRNA-target interactions contributes to our understanding of miRNA targeting specificity in non-model organisms. It also provides a roadmap for additional studies focused on regulatory functions of miRNAs in Anopheles gambiae.

Metazoan miRNAs typically bind to partially complementary sites in their target mRNAs. The interactions between miRNAs and target RNAs are generally stage-specific and context-dependent. Thus, identification of authentic miRNA targets remains a big challenge. Target identification is even more difficult in mosquitoes where miRNA-mRNA pairing rules are poorly characterized. Using an experimental approach, this study captures thousands of endogenous miRNA-target interactions in female mosquitoes at several critical stages during adult reproduction. Analyses of the target sequences reveal how individual miRNAs accomplish their target recognition in mosquitoes. Interestingly, many mosquito miRNAs exhibit flexibility to use distinct sequences at different stages to pair with their targets, greatly altering target selectivity and expanding target repertoire of miRNAs. Drastic changes in mRNA abundance have been previously reported when adult female mosquitoes attend to varying nutritional sources and physiological demands. The temporal patterns of miRNA-target interactions obtained in this study provide new insights into the roles of miRNAs in tightly controlled gene expression associated with blood-feeding and mosquito oogenesis.

Artificial Light at Night Influences Clock-Gene Expression, Activity, and Fecundity in the Mosquito Culex pipiens f. molestus

Light is an important environmental cue, and exposure to artificial light at night (ALAN) may disrupt organismal physiology and behavior. We investigated whether ALAN led to changes in clock-gene expression, diel activity patterns, and fecundity in laboratory populations of the mosquito Culex pipiens f. molestus (Diptera, Culicidae), a species that occurs in urban areas and is thus regularly exposed to ALAN. Populations were kept under 16hours (h):8h light:dark cycles or were subjected to an additional 3.5 h of light (100–300 lx) in the evenings. ALAN induced significant changes in expression in all genes studied, either alone (period) or as an interaction with time (timeless, cryptochrome2, Clock, cycle). Changes were sex-specific: period was down-regulated in both sexes, cycle was up-regulated in females, and Clock was down-regulated in males. ALAN-exposed mosquitoes were less active during the extra-light phase, but exposed females were more active later in the night. ALAN-exposed females also produced smaller and fewer eggs. Our findings indicate a sex-specific impact of ALAN on the physiology and behavior of Culex pipiens f. molestus and that changes in clock-gene expression, activity, and fecundity may be linked.

The contribution of gut bacteria to insecticide resistance and the life histories of the major malaria vector Anopheles arabiensis (Diptera: Culicidae)

The gut microbiota of mosquitoes is a crucial determinant of their fitness. As such, the biology of the gut microbiota of Anopheles arabiensis, a major malaria vector of Southern Africa, was investigated. Two laboratory strains of An. arabiensis were used; SENN, an insecticide susceptible strain, and SENN-DDT, a resistant strain. The strains were supplemented with either non-commensal bacteria or antibiotics via a sucrose source to sterilize the gut. The strains were fed the broad-spectrum bactericidal antibiotic gentamicin, or a preferentially gram-positive bactericidal (vancomycin), gram-negative bactericidal (streptomycin) or broad-spectrum bacteriostatic (erythromycin), either by sugar supplementation or by artificially-spiked blood-meal. The effects on adult mosquito longevity and insecticide resistance phenotype were assessed. Bacteria from the midgut of both strains were characterised by MALDI-TOF mass spectroscopy. Bactericidal antibiotics increased longevity in SENN-DDT. Bacterial supplementation increased insecticide tolerance. Antibiotic supplementation via sugar decreased tolerance to the insecticides deltamethrin and malathion. Blood-supplemented vancomycin decreased insecticide resistance, while gentamicin and streptomycin increased resistance. SENN showed a greater gut bacterial diversity than SENN-DDT, with both strains dominated by Gram-negative bacteria. This study suggests a crucial role for bacteria in An. arabiensis life history, and that gut microflora play variable roles in insecticide resistant and susceptible mosquitoes.

Small RNA-Seq Analysis Reveals miRNA Expression Dynamics Across Tissues in the Malaria Vector, Anopheles gambiae

Malaria continues to be a major global health problem, where disease transmission is deeply linked to the repeated blood feeding nature of the anautogenous mosquito. Given the tight link between blood feeding and disease transmission, understanding basic biology behind mosquito physiology is a requirement for developing effective vector-borne disease control strategies. In the mosquito, numerous loss of function studies with notable phenotypes demonstrate microRNAs (miRNAs) play significant roles in mosquito physiology. While the field appreciates the importance of a handful of miRNAs, we still need global mosquito tissue miRNA transcriptome studies. To address this need, our goal was to determine the miRNA transcriptome for multiple tissues of the pre-vitellogenic mosquito. To this end, by using small RNA-Seq analysis, we determined miRNA transcriptomes in tissues critical for mosquito reproduction and immunity including (i) fat body-abdominal wall enriched tissues, (ii) midguts, (iii) ovaries, and (iv) remaining tissues comprised of the head and thorax. We found numerous examples of miRNAs exhibiting pan-tissue high- or low- expression, tissue exclusion, and tissue enrichment. We also updated and consolidated the miRNA catalog and provided a detailed genome architecture map for the malaria vector, Anopheles gambiae. This study aims to build a foundation for future research on how miRNAs and potentially other small RNAs regulate mosquito physiology as it relates to vector-borne disease transmission.

Steroid Hormone Function Controls Non-competitive Plasmodium Development in Anopheles

Transmission of malaria parasites occurs when a female Anopheles mosquito feeds on an infected host to acquire nutrients for egg development. How parasites are affected by oogenetic processes, principally orchestrated by the steroid hormone 20-hydroxyecdysone (20E), remains largely unknown. Here we show that Plasmodium falciparum development is intimately but not competitively linked to processes shaping Anopheles gambiae reproduction. We unveil a 20E-mediated positive correlation between egg and oocyst numbers; impairing oogenesis by multiple 20E manipulations decreases parasite intensities. These manipulations, however, accelerate Plasmodium growth rates, allowing sporozoites to become infectious sooner. Parasites exploit mosquito lipids for faster growth, but they do so without further affecting egg development. These results suggest that P. falciparum has adopted a non-competitive evolutionary strategy of resource exploitation to optimize transmission while minimizing fitness costs to its mosquito vector. Our findings have profound implications for currently proposed control strategies aimed at suppressing mosquito populations.

Sexual Transcription Differences in Brachymeria lasus (Hymenoptera: Chalcididae), a Pupal Parasitoid Species of Lymantria dispar (Lepidoptera: Lymantriidae)

Sex differences in gene expression have been extensively documented, but little is known about these differences in parasitoid species that are widely applied to control pests. Brachymeria lasus is a solitary parasitoid species and has been evaluated as a potential candidate for release to control Lymantria dispar. In this study, gender differences in B. lasus were investigated using Illumina-based transcriptomic analysis. The resulting 37,453 unigene annotations provided a large amount of useful data for molecular studies of B. lasus. A total of 1416 differentially expressed genes were identified between females and males, and the majority of the sex-biased genes were female biased. Gene Ontology (GO) and Pathway enrichment analyses showed that (1) the functional categories DNA replication, fatty acid biosynthesis, and metabolism were enhanced in females and that (2) the only pathway enriched in males was phototransduction, while the GO subcategories enriched in males were those involved in membrane and ion transport. In addition, thirteen genes involving transient receptor potential (TRP) channels were annotated in B. lasus. We further explored and discussed the functions of TRPs in sensory signaling of light and temperature. In general, this study provides new molecular insights into the biological and sexually dimorphic traits of parasitoids, which may improve the application of these insects to the biological control of pests.

Identification, characterization and expression analysis of Anopheles stephensi double peroxidase

Peroxidases catalyze the reduction of peroxides and that, in turn, oxidize various substrates. They have been widely reported to play an important role in mosquito innate immunity against various pathogens. Here, we have characterized double heme peroxidase (AsDBLOX) gene from the Indian malaria vector Anopheles stephensi. It is a true ortholog of An. gambiae DBLOX. This 4209 bp AsDBLOX gene encodes for a protein of 1402 amino acids that has two duplicated peroxidase domains, domain I (from amino acid 61 to 527) and domain II (from amino acid 714 to 1252). The first domain has only substrate binding sites and lacks all other motifs of a functional heme peroxidase (e.g. heme binding site, calcium binding site and homodimer interface). Instead, it has two integrin binding motifs-LDV (Leu-Asp-Val) and RGD (Arg-Gly-Asp). The second peroxidase domain, however, has all the features of a complete heme peroxidase along with an integrin binding motif LDI (Leu-Asp-Ile). Thus, AsDBLOX gene is a unique type of peroxinectin as these groups of proteins are characterized by integrin binding motifs along with a heme peroxidase domain. We also observed that the AsDBLOX gene is expressed in all the life cycle stages of mosquito and is highly induced in the pupal stage of development which indicates its possible role in development.

Identification, Molecular Characterization, and In Silico Structural Analysis of Carboxypeptidase B2 of Anopheles stephensi

Malaria is a vector-borne infectious disease that is considered a priority of the World Health Organization due to its enormous impacts on global health. Plasmodium spp. (Haemosporida: Plasmodiidae), Anopheles spp. (Diptera: Culicidae), and a suitable host are the key elements for malaria transmission. To disrupt the parasitic life cycle of malaria or prevent its transmission, these three key elements should be targeted by effective control strategies. Development of vaccines that interrupt malaria transmission is one of the solutions that has been recommended to the countries that aim to eliminate malaria. With respect to the important role of Anopheles stephensi in malaria transmission and involvement of Anopheles carboxypeptidase B1 in sexual parasite development, we characterized the second member of cpb gene family (cpbAs2) of An. Stephensi to provide some basic information and evaluate significance of cpbAs2's role in complementing sexual plasmodium development role of cpbAs1. The cpbAs2 mRNA sequence was characterized by 3' and 5' RACE and the structural features of its coded protein were studied by in silico modeling. The coding sequence and gene structure of cpbAs2 were determined empirically and compared with the in silico predictions from the An. stephensi genome sequencing project. Furthermore, homology modeling revealed that its structure is very similar to the structurally important domains of procarboxypeptidase B2 in humans. This study provides basic molecular and structural information about another member of the cpb gene family of An. stephensi. The reported results are informative and necessary for evaluation of the role of this gene in sexual parasite development by future studies.

Functional Consequences of the Evolution of Matrimony, a Meiosis-Specific Inhibitor of Polo Kinase

Meiosis is a defining characteristic of eukaryotes, believed to have evolved only once, over one billion years ago. While the general progression of meiotic events is conserved across multiple diverse organisms, the specific pathways and proteins involved can be highly divergent, even within species from the same genus. Here we investigate the rapid evolution of Matrimony (Mtrm), a female meiosis-specific regulator of Polo kinase (Polo) in Drosophila. Mtrm physically interacts with Polo and is required to restrict the activity of Polo during meiosis. Despite Mtrm’s critical role in meiosis, sequence conservation within the genus Drosophila is poor. To explore the functional significance of this rapid divergence, we expressed Mtrm proteins from 12 different Drosophila species in the Drosophila melanogaster female germline. Distantly related Mtrm homologs are able to both physically interact with D. melanogaster Polo and rescue the meiotic defects seen in mtrm mutants. However, these distant homologs are not properly degraded after the completion of meiosis. Rather, they continue to inhibit Polo function in the early embryo, resulting in dominant maternal-effect lethality. We show that the ability of Mtrm to be properly degraded, and thus release Polo, is partially due to residues or motifs found within Mtrm’s least-conserved regions. We hypothesize that, while Mtrm regions critical for its meiotic function are under strong purifying selection, changes that occurred in its unconserved regions may have been advantageous, potentially by affecting the timing or duration of meiosis and/or the early embryonic divisions.

The Evolution and Metamorphosis of Arthropod Proteomics and Genomics

This article presents an overview of the development of techniques for analyzing cuticular proteins (CPs), their transcripts, and their genes over the past 50 years based primarily on experience in the laboratory of J.H. Willis. It emphasizes changes in the kind of data that can be gathered and how such data provided insights into the molecular underpinnings of insect metamorphosis and cuticle structure. It describes the techniques that allowed visualization of the location of CPs at both the anatomical and intracuticular levels and measurement of the appearance and deployment of transcripts from CP genes as well as what was learned from genomic and transcriptomic data. Most of the early work was done with the cecropia silkmoth, Hyalophora cecropia, and later work was with Anopheles gambiae.

How micronutrients influence the physiology of mosquitoes

Micronutrients or non-energetic nutrients (NEN) are needed in reduced amounts, but are essential for many mosquito physiological processes that influence biological traits from vector competence to reproductive capacity. The NEN include amino acids (AA), vitamins, salts, metals and sterols. Free AA plays critical roles controlling most physiological processes, from digestion to reproduction. Particularly proline connects metabolic pathways in energy production, flight physiology and ammonia detoxification. Metal, in particular iron and calcium, salts, sterol and vitamin homeostasis are critical for cell signaling, respiration, metabolism and reproduction. Micronutrient homeostasis influence the symbiotic relationships with microorganisms, having important implications in mosquitoes' nutrition, physiology and behavior, as well as in mosquito immunity and vector competence.

Association of microRNAs with Argonaute proteins in the malaria mosquito Anopheles gambiae after blood ingestion

Drastic changes in gene expression occur after adult female mosquitoes take a blood meal and use the nutrients for egg maturation. A growing body of evidence indicates that microRNAs (miRNAs) contribute to this tightly controlled tissue- and stage-specific gene expression. To investigate the role of miRNAs, we monitored miRNA expression in the mosquito Anopheles gambiae during the 72-h period immediately after blood feeding. We also measured the association of miRNAs with Argonaute 1 (Ago1) and Argonaute 2 (Ago2) to assess the functional status of individual miRNA species. Overall, 173 mature miRNAs were precipitated with Ago1 and Ago2, including 12 new miRNAs, the orthologs of which are found thus far only in other Anopheles species. Ago1 is the predominant carrier of miRNAs in Anopheles gambiae. The abundance and Ago loading of most of the mature miRNAs were relatively stable after blood ingestion. However, miRNAs of the miR-309/286/2944 cluster were considerably upregulated after blood feeding. Injection of the specific antagomir for miR-309 resulted in smaller developing oocytes and ultimately fewer eggs. In addition, the Ago association of some miRNAs was not proportional to their cellular abundance, suggesting that integration of miRNAs into the Ago complexes is regulated by additional mechanisms.

The impact of metagenomic interplay on the mosquito redox homeostasis

Mosquitoes are exposed to oxidative challenges throughout their life cycle. The primary challenge comes from a blood meal. The blood digestion turns the midgut into an oxidative environment, which imposes pressure not only on mosquito fecundity and other physiological traits but also on the microbiota in the midgut. During evolution, mosquitoes have developed numerous oxidative defense mechanisms to maintain redox homeostasis in the midgut. In addition to antioxidants, SOD, catalase, and glutathione system, sufficient supply of the reducing agent, NADPH, is vital for a successful defense against oxidative stress. Increasing evidence indicates that in response to oxidative stress, cells reconfigure metabolic pathways to increase the generation of NADPH through NADP-reducing networks including the pentose phosphate pathway and others. The microbial homeostasis is critical for the functional contributions to various host phenotypes. The symbiotic microbiota is regulated largely by the Duox-ROS pathway in Drosophila. In mosquitoes, Duox-ROS pathway, heme-mediated signaling, antimicrobial peptide production and C-type lectins work in concert to maintain the dynamic microbial community in the midgut. Microbial mechanisms against oxidative stress in this context are not well understood. Emerging evidence that microbial metabolites trigger host oxidative response warrants further study on the metagenomic interplay in an oxidative environment like mosquito gut ecosystem. Besides the classical Drosophila model, hematophagous insects like mosquitoes provide an alternative model system to study redox homeostasis in a symbiotic metagenomic context.

Evolutionary origin and function of NOX4-art, an arthropod specific NADPH oxidase

BACKGROUND

NADPH oxidases (NOX) are ROS producing enzymes that perform essential roles in cell physiology, including cell signaling and antimicrobial defense. This gene family is present in most eukaryotes, suggesting a common ancestor. To date, only a limited number of phylogenetic studies of metazoan NOXes have been performed, with few arthropod genes. In arthropods, only NOX5 and DUOX genes have been found and a gene called NOXm was found in mosquitoes but its origin and function has not been examined. In this study, we analyzed the evolution of this gene family in arthropods. A thorough search of genomes and transcriptomes was performed enabling us to browse most branches of arthropod phylogeny.

RESULTS

We have found that the subfamilies NOX5 and DUOX are present in all arthropod groups. We also show that a NOX gene, closely related to NOX4 and previously found only in mosquitoes (NOXm), can also be found in other taxonomic groups, leading us to rename it as NOX4-art. Although the accessory protein p22-phox, essential for NOX1-4 activation, was not found in any of the arthropods studied, NOX4-art of Aedes aegypti encodes an active protein that produces H₂O₂. Although NOX4-art has been lost in a number of arthropod lineages, it has all the domains and many signature residues and motifs necessary for ROS production and, when silenced, H₂O₂ production is considerably diminished in A. aegypti cells.

CONCLUSIONS

Combining all bioinformatic analyses and laboratory work we have reached interesting conclusions regarding arthropod NOX gene family evolution. NOX5 and DUOX are present in all arthropod lineages but it seems that a NOX2-like gene was lost in the ancestral lineage leading to Ecdysozoa. The NOX4-art gene originated from a NOX4-like ancestor and is functional. Although no p22-phox was observed in arthropods, there was no evidence of neo-functionalization and this gene probably produces H₂O₂ as in other metazoan NOX4 genes. Although functional and present in the genomes of many species, NOX4-art was lost in a number of arthropod lineages.

The Anopheles FBN9 immune factor mediates Plasmodium species-specific defense through transgenic fat body expression

Mosquitoes have a multifaceted innate immune system that is actively engaged in warding off various pathogens, including the protozoan malaria parasite Plasmodium. Various immune signaling pathways and effectors have been shown to mediate a certain degree of defense specificity against different Plasmodium species. A key pattern recognition receptor of the Anopheles gambiae immune system is the fibrinogen domain-containing immunolectin FBN9, which has been shown to be transcriptonally induced by Plasmodium infection, and to mediate defense against both rodent and human malaria parasites and bacteria. Here we have further studied the defense specificity of FBN9 using a transgenic approach, in which FBN9 is overexpressed in the fat body tissue after a blood meal through a vitellogenin promoter. Interestingly, the Vg-FBN9 transgenic mosquitoes showed increased resistance only to the rodent parasite P. berghei, and not to the human parasite P. falciparum, pointing to differences in the mosquito's defense mechanisms against the two parasite species. The Vg-FBN9 transgenic mosquitoes were also more resistant to infection with both Gram-positive and Gram-negative bacteria and showed increased longevity when infected with P. berghei. Our study points to the importance of both experimentally depleting and enriching candidate anti-Plasmodium effectors in functional studies in order to ascertain their suitability for the development of transgenic mosquito-based malaria control strategies.

Transgenic Expression of the Anti-parasitic Factor TEP1 in the Malaria Mosquito Anopheles gambiae

Mosquitoes genetically engineered to be resistant to Plasmodium parasites represent a promising novel approach in the fight against malaria. The insect immune system itself is a source of anti-parasitic genes potentially exploitable for transgenic designs. The Anopheles gambiae thioester containing protein 1 (TEP1) is a potent anti-parasitic protein. TEP1 is secreted and circulates in the mosquito hemolymph, where its activated cleaved form binds and eliminates malaria parasites. Here we investigated whether TEP1 can be used to create malaria resistant mosquitoes. Using a GFP reporter transgene, we determined that the fat body is the main site of TEP1 expression. We generated transgenic mosquitoes that express TEP1r, a potent refractory allele of TEP1, in the fat body and examined the activity of the transgenic protein in wild-type or TEP1 mutant genetic backgrounds. Transgenic TEP1r rescued loss-of-function mutations, but did not increase parasite resistance in the presence of a wild-type susceptible allele. Consistent with previous reports, TEP1 protein expressed from the transgene in the fat body was taken up by hemocytes upon a challenge with injected bacteria. Furthermore, although maturation of transgenic TEP1 into the cleaved form was impaired in one of the TEP1 mutant lines, it was still sufficient to reduce parasite numbers and induce parasite melanization. We also report here the first use of Transcription Activator Like Effectors (TALEs) in Anopheles gambiae to stimulate expression of endogenous TEP1. We found that artificial elevation of TEP1 expression remains moderate in vivo and that enhancement of endogenous TEP1 expression did not result in increased resistance to Plasmodium. Taken together, our results reveal the difficulty of artificially influencing TEP1-mediated Plasmodium resistance, and contribute to further our understanding of the molecular mechanisms underlying mosquito resistance to Plasmodium parasites.

RNA Interference for Mosquito and Mosquito-Borne Disease Control

RNA interference (RNAi) is a powerful tool to silence endogenous mosquito and mosquito-borne pathogen genes in vivo. As the number of studies utilizing RNAi in basic research grows, so too does the arsenal of physiological targets that can be developed into products that interrupt mosquito life cycles and behaviors and, thereby, relieve the burden of mosquitoes on human health and well-being. As this technology becomes more viable for use in beneficial and pest insect management in agricultural settings, it is exciting to consider its role in public health entomology. Existing and burgeoning strategies for insecticide delivery could be adapted to function as RNAi trigger delivery systems and thereby expedite transformation of RNAi from the lab to the field for mosquito control. Taken together, development of RNAi-based vector and pathogen management techniques & strategies are within reach. That said, tools for successful RNAi design, studies exploring RNAi in the context of vector control, and studies demonstrating field efficacy of RNAi trigger delivery have yet to be honed and/or developed for mosquito control.

Larval diet affects mosquito development and permissiveness to Plasmodium infection

The larval stages of malaria vector mosquitoes develop in water pools, feeding mostly on microorganisms and environmental detritus. Richness in the nutrient supply to larvae influences the development and metabolism of larvae and adults. Here, we investigated the effects of larval diet on the development, microbiota content and permissiveness to Plasmodium of Anopheles coluzzii. We tested three fish diets often used to rear mosquitoes in the laboratory, including two pelleted diets, Dr. Clarke’s Pool Pellets and Nishikoi Fish Pellets, and one flaked diet, Tetramin Fish-Flakes. Larvae grow and develop faster and produce bigger adults when feeding on both types of pellets compared with flakes. This correlates with a higher microbiota load in pellet-fed larvae, in agreement with the known positive effect of the microbiota on mosquito development. Larval diet also significantly influences the prevalence and intensity of Plasmodium berghei infection in adults, whereby Nishikoi Fish Pellets-fed larvae develop into adults that are highly permissive to parasites and survive longer after infection. This correlates with a lower amount of Enterobacteriaceae in the midgut microbiota. Together, our results shed light on the influence of larval feeding on mosquito development, microbiota and vector competence; they also provide useful data for mosquito rearing.

1.45 Å resolution structure of SRPN18 from the malaria vector Anopheles gambiae

Serine protease inhibitors (serpins) in insects function within development, wound healing and immunity. The genome of the African malaria vector, Anopheles gambiae, encodes 23 distinct serpin proteins, several of which are implicated in disease-relevant physiological responses. A. gambiae serpin 18 (SRPN18) was previously categorized as non-inhibitory based on the sequence of its reactive-center loop (RCL), a region responsible for targeting and initiating protease inhibition. The crystal structure of A. gambiae SRPN18 was determined to a resolution of 1.45 Å, including nearly the entire RCL in one of the two molecules in the asymmetric unit. The structure reveals that the SRPN18 RCL is extremely short and constricted, a feature associated with noncanonical inhibitors or non-inhibitory serpin superfamily members. Furthermore, the SRPN18 RCL does not contain a suitable protease target site and contains a large number of prolines. The SRPN18 structure therefore reveals a unique RCL architecture among the highly conserved serpin fold.

The Genetic Basis of Host Preference and Resting Behavior in the Major African Malaria Vector, Anopheles arabiensis

Malaria transmission is dependent on the propensity of Anopheles mosquitoes to bite humans (anthropophily) instead of other dead end hosts. Recent increases in the usage of Long Lasting Insecticide Treated Nets (LLINs) in Africa have been associated with reductions in highly anthropophilic and endophilic vectors such as Anopheles gambiae s.s., leaving species with a broader host range, such as Anopheles arabiensis, as the most prominent remaining source of transmission in many settings. An. arabiensis appears to be more of a generalist in terms of its host choice and resting behavior, which may be due to phenotypic plasticity and/or segregating allelic variation. To investigate the genetic basis of host choice and resting behavior in An. arabiensis we sequenced the genomes of 23 human-fed and 25 cattle-fed mosquitoes collected both in-doors and out-doors in the Kilombero Valley, Tanzania. We identified a total of 4,820,851 SNPs, which were used to conduct the first genome-wide estimates of "SNP heritability" for host choice and resting behavior in this species. A genetic component was detected for host choice (human vs cow fed; permuted P = 0.002), but there was no evidence of a genetic component for resting behavior (indoors versus outside; permuted P = 0.465). A principal component analysis (PCA) segregated individuals based on genomic variation into three groups which were characterized by differences at the 2Rb and/or 3Ra paracentromeric chromosome inversions. There was a non-random distribution of cattle-fed mosquitoes between the PCA clusters, suggesting that alleles linked to the 2Rb and/or 3Ra inversions may influence host choice. Using a novel inversion genotyping assay, we detected a significant enrichment of the standard arrangement (non-inverted) of 3Ra among cattle-fed mosquitoes (N = 129) versus all non-cattle-fed individuals (N = 234; χ2, p = 0.007). Thus, tracking the frequency of the 3Ra in An. arabiensis populations may be of use to infer selection on host choice behavior within these vector populations; possibly in response to vector control. Controlled host-choice assays are needed to discern whether the observed genetic component has a direct relationship with innate host preference. A better understanding of the genetic basis for host feeding behavior in An. arabiensis may also open avenues for novel vector control strategies based on driving genes for zoophily into wild mosquito populations.

Characterization of the target of ivermectin, the glutamate-gated chloride channel, from Anopheles gambiae

The use of insecticide-treated nets and indoor residual insecticides targeting adult mosquito vectors is a key element in malaria control programs. However, mosquito resistance to the insecticides used in these applications threatens malaria control efforts. Recently, the mass drug administration of ivermectin (IVM) has been shown to kill Anopheles gambiae mosquitoes and disrupt Plasmodium falciparum transmission in the field. We cloned the molecular target of IVM from A. gambiae, the glutamate-gated chloride channel (AgGluCl), and characterized its transcriptional patterns, protein expression and functional responses to glutamate and IVM. AgGluCl cloning revealed an unpredicted fourth splice isoform as well as a novel exon and splice site. The predicted gene products contained heterogeneity in the N-terminal extracellular domain and the intracellular loop region. Responses to glutamate and IVM were measured using two-electrode voltage clamp on Xenopus laevis oocytes expressing AgGluCl. IVM induced non-persistent currents in AgGluCl-a1 and did not potentiate glutamate responses. In contrast, AgGluCl-b was insensitive to IVM, suggesting that the AgGluCl gene could produce IVM-sensitive and -insensitive homomultimers from alternative splicing. AgGluCl isoform-specific transcripts were measured across tissues, ages, blood feeding status and sex, and were found to be differentially transcribed across these physiological variables. Lastly, we stained adult, female A. gambiae for GluCl expression. The channel was expressed in the antenna, Johnston's organ, supraesophageal ganglion and thoracic ganglia. In summary, we have characterized the first GluCl from a mosquito, A. gambiae, and described its unique activity and expression with respect to it as the target of the insecticide IVM.

CLIPB8 is part of the prophenoloxidase activation system in Anopheles gambiae mosquitoes

In insects and other arthropods the formation of eumelanin (melanization) is a broad spectrum and potent immune response that is used to encapsulate and kill invading pathogens. This immune response is regulated by the activation of prophenoxidase (proPO), which is controlled by proteinase cascades and its serpin inhibitors, together forming the proPO activation system. While the molecular composition of these protease cascades are well understood in insect model systems, major knowledge gaps remain in mosquitoes. Recently, a regulatory unit of melanization in Anopheles gambiae was documented, comprised of the inhibitory serpin-clip-serine proteinase, CLIPB9 and its inhibitor serpin-2 (SRPN2). Partial reversion of SRPN2 phenotypes in melanotic tumor formation and adult survival by SRPN2/CLIPB9 double knockdown suggested other target proteinases of SRPN2 in regulating melanization. Here we report that CLIPB8 supplements the SRPN2/CLIPB9 regulatory unit in controlling melanization in An. gambiae. As with CLIPB9, knockdown of CLIPB8 partially reversed the pleiotropic phenotype induced by SRPN2 silencing with regards to adult survival and melanotic tumor formation. Recombinant SRPN2 protein formed an SDS-stable protein complex with activated recombinant CLIPB8, however did not efficiently inhibit CLIPB8 activity in vitro. CLIPB8 did not directly activate proPO in vitro nor was it able to cleave and activate proCLIPB9. Nevertheless, epistasis analysis using RNAi placed CLIPB8 and CLIPB9 in the same pathway leading to melanization, suggesting that CLIPB8 either acts further upstream of CLIPB9 or is required for activation of a yet to be identified serine proteinase homolog. Taken together, this study identifies CLIPB8 as an additional player in proPO activation cascade and highlights the complexity of the proteinase network that regulates melanization in An. gambiae.

Plasticity in Insect Olfaction: To Smell or Not to Smell?

In insects, olfaction plays a crucial role in many behavioral contexts, such as locating food, sexual partners, and oviposition sites. To successfully perform such behaviors, insects must respond to chemical stimuli at the right moment. Insects modulate their olfactory system according to their physiological state upon interaction with their environment. Here, we review the plasticity of behavioral responses to different odor types according to age, feeding state, circadian rhythm, and mating status. We also summarize what is known about the underlying neural and endocrinological mechanisms, from peripheral detection to central nervous integration, and cover neuromodulation from the molecular to the behavioral level. We describe forms of olfactory plasticity that have contributed to the evolutionary success of insects and have provided them with remarkable tools to adapt to their ever-changing environment.

A comprehensive transcriptomic view of renal function in the malaria vector, Anopheles gambiae

Renal function is essential to maintain homeostasis. This is particularly significant for insects that undergo complete metamorphosis; larval mosquitoes must survive a freshwater habitat whereas adults are terrestrial, and mature females must maintain ion and fluid homeostasis after blood feeding. To investigate the physiological adaptations required for successful development to adulthood, we studied the Malpighian tubule transcriptome of Anopheles gambiae using Affymetrix arrays. We assessed transcription under several conditions; as third instar larvae, as adult males fed on sugar, as adult females fed on sugar, and adult females after a blood meal. In addition to providing the most detailed transcriptomic data to date on the Anopheles Malpighian tubules, the data provide unique information on the renal adaptations required for the switch from freshwater to terrestrial habitats, on gender differences, and on the contrast between nectar-feeding and haematophagy. We found clear differences associated with ontogenetic change in lifestyle, gender and diet, particularly in the neuropeptide receptors that control fluid secretion, and the water and ion transporters that impact volume and composition. These data were also combined with transcriptomics from the Drosophila melanogaster tubule, allowing meta-analysis of the genes which underpin tubule function across Diptera. To further investigate renal conservation across species we selected four D. melanogaster genes with orthologues highly enriched in the Anopheles tubules, and generated RNAi knockdown flies. Three of these genes proved essential, showing conservation of critical functions across 150 million years of phylogenetic separation. This extensive data-set is available as an online resource, MozTubules.org, and could potentially be mined for novel insecticide targets that can impact this critical organ in this pest species.

Increased production of piRNAs from euchromatic clusters and genes in Anopheles gambiae compared with Drosophila melanogaster

BACKGROUND

Specific genomic loci, termed Piwi-interacting RNA (piRNA) clusters, manufacture piRNAs that serve as guides for the inactivation of complementary transposable elements (TEs). The piRNA pathway has been accurately detailed in Drosophila melanogaster, while it remains poorly examined in other insects. This pathway is increasingly recognized as critical for germline development and reproduction. Understanding of the piRNA functions in mosquitoes could offer an opportunity for disease vector control by the reduction of their reproductive potential.

RESULTS

To analyze the similarities and differences in this pathway between Drosophila and mosquito, we performed an in-depth analysis of the genomic loci producing piRNAs and their targets in the African malaria vector Anopheles gambiae. We identified 187 piRNA clusters in the An. gambiae genome and 155 piRNA clusters in the D. melanogaster genome. We demonstrate that many more piRNA clusters in the mosquito compared with the fruit fly are uni-directionally transcribed and are located outside pericentromeric heterochromatin. About 11 % of the An. gambiae piRNA population map to gene transcripts. This is a noticeable increase compared with the ~6 % of the piRNA population mapped to genes in D. melanogaster. A subset of the piRNA-enriched genes in An. gambiae has functions related to reproduction and development. At least 24 and 65 % of the mapped piRNAs correspond to genomic TE sequences in An. gambiae and D. melanogaster, respectively. DNA transposons and non-LTR retrotransposons are more abundant in An. gambiae, while LTR retrotransposons are more abundant in D. melanogaster. Yet, piRNAs predominantly target LTR retrotransposons in both species, which may point to a distinct feature of these elements compared to the other classes of TEs concerning their silencing by the piRNA pathway.

CONCLUSIONS

Here, we demonstrate that piRNA-producing loci have more ubiquitous distribution in the An. gambiae genome than in the genome of D. melanogaster. Also, protein-coding genes have an increased role in production of piRNAs in the germline of this mosquito. Genes involved in germline and embryonic development of An. gambiae generate a substantial portion of piRNAs, suggesting a role of the piRNA pathway in the epigenetic regulation of the reproductive processes in the African malaria vector.

Age and prior blood feeding of Anopheles gambiae influences their susceptibility and gene expression patterns to ivermectin-containing blood meals

Background

Ivermectin has been proposed as a novel malaria transmission control tool based on its insecticidal properties and unique route of acquisition through human blood. To maximize ivermectin’s effect and identify potential resistance/tolerance mechanisms, it is important to understand its effect on mosquito physiology and potential to shift mosquito population age-structure. We therefore investigated ivermectin susceptibility and gene expression changes in several age groups of female Anopheles gambiae mosquitoes.

Methods

The effect of aging on ivermectin susceptibility was analyzed in three age groups (2, 6, and 14-days) of colonized female Anopheles gambiae mosquitoes using standard survivorship assays. Gene expression patterns were then analyzed by transcriptome sequencing on an Illumina HiSeq 2500 platform. RT-qPCR was used to validate transcriptional changes and also to examine expression in a different, colonized strain and in wild mosquitoes, both of which blood fed naturally on an ivermectin-treated person.

Results

Mosquitoes of different ages and blood meal history died at different frequencies after ingesting ivermectin. Mortality was lowest in 2-day old mosquitoes exposed on their first blood meal and highest in 6-day old mosquitoes exposed on their second blood meal. Twenty-four hours following ivermectin ingestion, 101 and 187 genes were differentially-expressed relative to control blood-fed, in 2 and 6-day groups, respectively. Transcription patterns of select genes were similar in membrane-fed, colonized, and naturally-fed wild vectors. Transcripts from several unexpected functional classes were highly up-regulated, including Niemann-Pick Type C (NPC) genes, peritrophic matrix-associated genes, and immune-response genes, and these exhibited different transcription patterns between age groups, which may explain the observed susceptibility differences. Niemann-Pick Type 2 genes were the most highly up-regulated transcripts after ivermectin ingestion (up to 160 fold) and comparing phylogeny to transcriptional patterns revealed that NPCs have rapidly evolved and separate members respond to either blood meals or to ivermectin.

Conclusion

We present evidence of increased ivermectin susceptibility in older An. gambiae mosquitoes that had previously bloodfed. Differential expression analysis suggests complex midgut interactions resulting from ivermectin ingestion that likely involve blood meal digestion physiological responses, midgut microflora, and innate immune responses. Thus, the transcription of certain gene families is consistently affected by ivermectin ingestion, and may provide important clues to ivermectin’s broad effects on malaria vectors. These findings contribute to the growing understanding of ivermectin’s potential as a transmission control tool.

Electronic supplementary material

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

Association mapping by pooled sequencing identifies TOLL 11 as a protective factor against Plasmodium falciparum in Anopheles gambiae

Background

The genome-wide association study (GWAS) techniques that have been used for genetic mapping in other organisms have not been successfully applied to mosquitoes, which have genetic characteristics of high nucleotide diversity, low linkage disequilibrium, and complex population stratification that render population-based GWAS essentially unfeasible at realistic sample size and marker density.

Methods

We designed a novel mapping strategy for the mosquito system that combines the power of linkage mapping with the resolution afforded by genetic association. We established founder colonies from West Africa, controlled for diversity, linkage disequilibrium and population stratification. Colonies were challenged by feeding on the infectious stage of the human malaria parasite, Plasmodium falciparum, mosquitoes were phenotyped for parasite load, and DNA pools for phenotypically similar mosquitoes were Illumina sequenced. Phenotype-genotype mapping was carried out in two stages, coarse and fine.

Results

In the first mapping stage, pooled sequences were analysed genome-wide for intervals displaying relativereduction in diversity between phenotype pools, and candidate genomic loci were identified for influence upon parasite infection levels. In the second mapping stage, focused genotyping of SNPs from the first mapping stage was carried out in unpooled individual mosquitoes and replicates. The second stage confirmed significant SNPs in a locus encoding two Toll-family proteins. RNAi-mediated gene silencing and infection challenge revealed that TOLL 11 protects mosquitoes against P. falciparum infection.

Conclusions

We present an efficient and cost-effective method for genetic mapping using natural variation segregating in defined recent Anopheles founder colonies, and demonstrate its applicability for mapping in a complex non-model genome. This approach is a practical and preferred alternative to population-based GWAS for first-pass mapping of phenotypes in Anopheles. This design should facilitate mapping of other traits involved in physiology, epidemiology, and behaviour.

Electronic supplementary material

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

Temporal Gene Expression Profiles of Pre Blood-Fed Adult Females Immediately Following Eclosion in the Southern House Mosquito Culex Quinquefasciatus

Prior to acquisition of the first host blood meal, the anautogenous mosquito Culex quinquefasciatus requires a period of time in order to prepare for the blood feeding and, later, vitellogenesis. In the current study, we conducted whole transcriptome analyses of adult female Culex mosquitoes to identify genes that may be necessary for both taking of the blood meal, and processing of the blood meal in adult female mosquitoes Cx. quinquefasciatus. We examined temporal expression of genes for the periods of post eclosion and prior to the female freely taking a blood meal. We further evaluated the temporal expression of certain genes for the periods after the taking of a blood meal to identify genes that may be necessary for both the taking of the blood meal, and the processing of the blood meal. We found that adult females required a minimum of 48 h post-eclosion before they freely took their first blood meal. We hypothesized that gene expression signatures were altered in the mosquitoes before blood feeding in preparation for the acquisition of the blood meal through changes in multiple gene expression. To identify the genes involved in the acquisition of blood feeding, we quantified the gene expression levels of adult female Cx. quinquefasciatus using RNA Seq throughout a pre-blooding period from 2 to 72 h post eclosion at 12 h intervals. A total of 325 genes were determined to be differentially-expressed throughout the pre-blooding period, with the majority of differentially-expressed genes occurring between the 2 h and 12 h post-eclosion time points. Among the up-regulated genes were salivary proteins, cytochrome P450s, odorant-binding proteins, and proteases, while the majority of the down-regulated genes were hypothetical or cuticular genes. In addition, Trypsin was found to be up-regulated immediately following blood feeding, while trypsin and chymotrypsin were up-regulated at 48h and 60h post blood-feeding, respectively, suggesting that these proteases are likely involved in the digestion of the blood meal. Overall, this study reviewed multiple genes that might be involved in the adult female competency for blood meal acquisition in mosquitoes.

Regulation of Gene Expression Patterns in Mosquito Reproduction

In multicellular organisms, development, growth and reproduction require coordinated expression of numerous functional and regulatory genes. Insects, in addition to being the most speciose animal group with enormous biological and economical significance, represent outstanding model organisms for studying regulation of synchronized gene expression due to their rapid development and reproduction. Disease-transmitting female mosquitoes have adapted uniquely for ingestion and utilization of the huge blood meal required for swift reproductive events to complete egg development within a 72-h period. We investigated the network of regulatory factors mediating sequential gene expression in the fat body, a multifunctional organ analogous to the vertebrate liver and adipose tissue, of the female Aedes aegypti mosquito. Transcriptomic and bioinformatics analyses revealed that ~7500 transcripts are differentially expressed in four sequential waves during the 72-h reproductive period. A combination of RNA-interference gene-silencing and in-vitro organ culture identified the major regulators for each of these waves. Amino acids (AAs) regulate the first wave of gene activation between 3 h and 12 h post-blood meal (PBM). During the second wave, between 12 h and 36 h, most genes are highly upregulated by a synergistic action of AAs, 20-hydroxyecdysone (20E) and the Ecdysone-Receptor (EcR). Between 36 h and 48 h, the third wave of gene activation—regulated mainly by HR3—occurs. Juvenile Hormone (JH) and its receptor Methoprene-Tolerant (Met) are major regulators for the final wave between 48 h and 72 h. Each of these key regulators also has repressive effects on one or more gene sets. Our study provides a better understanding of the complexity of the regulatory mechanisms related to temporal coordination of gene expression during reproduction. We have detected the novel function of 20E/EcR responsible for transcriptional repression. This study also reveals the previously unidentified large-scale effects of HR3 and JH/Met on transcriptional regulation during the termination of vitellogenesis and remodeling of the fat body.

In addition to being vectors of devastating human diseases, mosquitoes represent outstanding model organisms for studying regulatory mechanisms of differential gene expression due to their rapid reproductive cycles. About 7500 transcripts are differentially expressed in four sequential waves during the 72-h reproductive period in the fat body, a critical reproductive organ. The major regulators for these waves of gene expression are the two very important insect hormones, 20-hydroxyecdysone (20E) and Juvenile hormone (JH), their respective receptors Ecdysone Receptor (EcR) and Methoprene-Tolerant (Met), amino acids and the orphan nuclear receptor HR3. These key regulators are responsible for activation and repression of co-regulated gene sets, at different time points, within the 72-h reproductive period. Importantly, this study, apart from providing an insight into the regulatory complexity involved in the temporal coordination of gene expression, also reveals the previously unidentified roles of 20E/EcR, JH/Met and HR3 during the 72-h period post blood meal.