aeGEPUCI: a database of gene expression in the dengue vector mosquito, Aedes aegypti

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post-blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codon decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we discovered that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

Unified Aedes aegypti Protein Resource Database (UAAPRD): An Integrated High-Throughput In Silico Platform for Comprehensive Protein Structure Modeling and Functional Target Analysis to Enhance Vector Control Strategies

A comprehensive examination of Aedes aegypti's proteome to detect key proteins that can be targeted with small molecules can disrupt blood feeding and disease transmission. However, research currently only focuses on finding repellent-like compounds, limiting studies on identifying unexplored proteins in its proteome. High-throughput analysis generates vast amounts of data, raising concerns about accessibility and usability. Establishing a dedicated database is a solution, centralizing information on identified proteins, functions, and modeled structures for easy access and research. This study focuses on scrutinizing key proteins in A. aegypti, modeling their structures using RaptorX standalone tool, identification of druggable binding sites using BiteNet, validating the models via Ramachandran plot studies and refining them via 50-ns molecular dynamic simulations using Schrodinger Maestro. By analyzing ~ 18 k proteins in the proteome of A. aegypti in our previous studies, all proteins involved in the light and dark circadian rhythm of the mosquito, inclusive of proteins in blood feeding, metabolism, etc. were chosen for the current study. The outcome is UAAPRD, a unique repository housing information on hundreds of previously unmodeled and un-simulated mosquito proteins. This robust MYSQL database ( https://uaaprd.onrender.com/user ) houses data on 309 modeled & simulated proteins of A. aegypti. It allows users to obtain protein data, view evolutionary analysis data of the protein categories, visualize proteins of interest, and send request to screen against the pharmacophore models present in UAAPRD against ligand of interest. This study offers crucial insights for developing targeted studies, which will ultimately contribute to more effective vector control strategies.

Comprehensive Proteomics Analysis of the Hemolymph Composition of Sugar-Fed Aedes aegypti Female and Male Mosquitoes

In arthropods, hemolymph carries immune cells and solubilizes and transports nutrients, hormones, and other molecules that are involved in diverse physiological processes including immunity, metabolism, and reproduction. However, despite such physiological importance, little is known about its composition. We applied mass spectrometry-based label-free quantification approaches to study the proteome of hemolymph perfused from sugar-fed female and male Aedes aegypti mosquitoes. A total of 1403 proteins were identified, out of which 447 of them were predicted to be extracellular. In both sexes, almost half of these extracellular proteins were predicted to be involved in defense/immune response, and their relative abundances (based on their intensity-based absolute quantification, iBAQ) were 37.9 and 33.2%, respectively. Interestingly, among them, 102 serine proteases/serine protease-homologues were identified, with almost half of them containing CLIP regulatory domains. Moreover, proteins belonging to families classically described as chemoreceptors, such as odorant-binding proteins (OBPs) and chemosensory proteins (CSPs), were also highly abundant in the hemolymph of both sexes. Our data provide a comprehensive catalogue of A. aegypti hemolymph basal protein content, revealing numerous unexplored targets for future research on mosquito physiology and disease transmission. It also provides a reference for future studies on the effect of blood meal and infection on hemolymph composition.

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA (mRNA) codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codons' decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we identified that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

Ovariole-specific Yellow-g and Yellow-g2 proteins are required for fecundity and egg chorion rigidity in the red flour beetle, Tribolium castaneum

Most insects reproduce by laying eggs that have an eggshell/chorion secreted by follicle cells, which serves as a protective barrier for developing embryos. Thus, eggshell formation is vital for reproduction. Insect yellow family genes encode for secreted extracellular proteins that perform different, context-dependent functions in different tissues at various stages of development involving, for example, cuticle/eggshell coloration and morphology, molting, courtship behavior and embryo hatching. In this study we investigated the function of two of this family's genes, yellow-g (TcY-g) and yellow-g2 (TcY-g2), on the formation and morphology of the eggshell of the red flour beetle, Tribolium castaneum. Real-time PCR analysis revealed that both TcY-g and TcY-g2 were specifically expressed in the ovarioles of adult females. Loss of function produced by injection of double-stranded RNA (dsRNA) for either TcY-g or TcY-g2 gene resulted in failure of oviposition. There was no effect on maternal survival. Ovaries dissected from those dsRNA-treated females exhibited ovarioles containing not only developing oocytes but also mature eggs in their egg chambers. However, the ovulated eggs were collapsed and ruptured, resulting in swollen lateral oviducts and calyxes. TEM analysis showed that lateral oviducts were filled with electron-dense material, presumably from some cellular content leakage out of the collapsed eggs. In addition, morphological abnormalities in lateral oviduct epithelial cells and the tubular muscle sheath were evident. These results support the hypothesis that both TcY-g and TcY-g2 proteins are required for maintaining the rigidity and integrity of the chorion, which is critical for resistance to mechanical stress and/or rehydration during ovulation and egg activation in the oviducts of T. castaneum. Because Yellow-g and Yellow-g2 are highly conserved among insect species, both genes are potential targets for development of gene-based insect pest population control methods.

Genome-wide Association Study Reveals New Loci Associated With Pyrethroid Resistance in Aedes aegypti

Genome-wide association studies (GWAS) use genetic polymorphism across the genomes of individuals with distinct characteristics to identify genotype-phenotype associations. In mosquitoes, complex traits such as vector competence and insecticide resistance could benefit from GWAS. We used the Aedes aegypti 50k SNP chip to genotype populations with different levels of pyrethroid resistance from Northern Brazil. Pyrethroids are widely used worldwide to control mosquitoes and agricultural pests, and their intensive use led to the selection of resistance phenotypes in many insects including mosquitoes. For Ae. aegypti, resistance phenotypes are mainly associated with several mutations in the voltage-gated sodium channel, known as knockdown resistance (kdr). We phenotyped those populations with the WHO insecticide bioassay using deltamethrin impregnated papers, genotyped the kdr alleles using qPCR, and determined allele frequencies across the genome using the SNP chip. We identified single-nucleotide polymorphisms (SNPs) directly associated with resistance and one epistatic SNP pair. We also observed that the novel SNPs correlated with the known kdr genotypes, although on different chromosomes or not in close physical proximity to the voltage gated sodium channel gene. In addition, pairwise comparison of resistance and susceptible mosquitoes from each population revealed differentiated genomic regions not associated with pyrethroid resistance. These new bi-allelic markers can be used to genotype other populations along with kdr alleles to understand their worldwide distribution. The functional roles of the genes near the newly discovered SNPs require new studies to determine if they act synergistically with kdr alleles or reduce the fitness cost of maintaining resistant alleles.

Yellow-y Functions in Egg Melanization and Chorion Morphology of the Asian Tiger Mosquito, Aedes albopictus

The Asian tiger mosquito, Aedes albopictus, is one of the most serious public health pests, which can transmit various vector-borne diseases. Eggs from this mosquito species become dark black shortly after oviposition and exhibit high desiccation resistance. Some of the Yellow proteins that act as dopachrome conversion enzymes (DCEs) are involved in the tyrosine-mediated tanning (pigmentation and sclerotization) metabolic pathway that significantly accelerates melanization reactions in insects. In this research, we analyzed the function of one of the yellow genes, yellow-y (AalY-y), in eggshell/chorion melanization of Ae. albopictus eggs. Developmental and tissue-specific expression measured by real-time PCR showed that AalY-y transcripts were detected at all stages of development analyzed, with significantly higher levels in the ovaries from blood-fed adult females. Injection of double-stranded RNA for AalY-y (dsAalY-y) had no significant effect on fecundity. However, unlike dsEGFP-treated control eggs that become black by 2–3 h after oviposition (HAO), dsAalY-y eggs were yellow-brown at 2 HAO, and reddish-brown even at 48 HAO. dsEGFP eggs exhibited resistance to desiccation at 48 HAO, whereas approximately 50% of the dsAalY-y eggs collapsed when they were moved to a low humidity condition. In addition, TEM analysis revealed an abnormal morphology and ultrastructure of the outer-endochorion in the dsAalY-y eggs. These results support the hypothesis that AalY-y is involved in the tyrosine-induced melanin biosynthetic pathway, plays an important role in black melanization of the chorion and functions in conferring proper morphology of the outer-endochorion, a structure that is presumably required for egg desiccation resistance in Ae. albopictus.

Regulation of the antennal transcriptome of the dengue vector, Aedes aegypti, during the first gonotrophic cycle

Background

In the light of dengue being the fastest growing transmissible disease, there is a dire need to identify the mechanisms regulating the behaviour of the main vector Aedes aegypti. Disease transmission requires the female mosquito to acquire the pathogen from a blood meal during one gonotrophic cycle, and to pass it on in the next, and the capacity of the vector to maintain the disease relies on a sustained mosquito population.

Results

Using a comprehensive transcriptomic approach, we provide insight into the regulation of the odour-mediated host- and oviposition-seeking behaviours throughout the first gonotrophic cycle. We provide clear evidence that the age and state of the female affects antennal transcription differentially. Notably, the temporal- and state-dependent patterns of differential transcript abundance of chemosensory and neuromodulatory genes extends across families, and appears to be linked to concerted differential modulation by subsets of transcription factors.

Conclusions

By identifying these regulatory pathways, we provide a substrate for future studies targeting subsets of genes across disparate families involved in generating key vector behaviours, with the goal to develop novel vector control tools.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07336-w.

The SUMOylation pathway suppresses arbovirus replication in Aedes aegypti cells

Mosquitoes are responsible for the transmission of many clinically important arboviruses that cause significant levels of annual mortality and socioeconomic health burden worldwide. Deciphering the mechanisms by which mosquitoes modulate arbovirus infection is crucial to understand how viral-host interactions promote vector transmission and human disease. SUMOylation is a post-translational modification that leads to the covalent attachment of the Small Ubiquitin-like MOdifier (SUMO) protein to host factors, which in turn can modulate their stability, interaction networks, sub-cellular localisation, and biochemical function. While the SUMOylation pathway is known to play a key role in the regulation of host immune defences to virus infection in humans, the importance of this pathway during arbovirus infection in mosquito vectors, such as Aedes aegypti (Ae. aegypti), remains unknown. Here we characterise the sequence, structure, biochemical properties, and tissue-specific expression profiles of component proteins of the Ae. aegypti SUMOylation pathway. We demonstrate significant biochemical differences between Ae. aegypti and Homo sapiens SUMOylation pathways and identify cell-type specific patterns of SUMO expression in Ae. aegypti tissues known to support arbovirus replication. Importantly, depletion of core SUMOylation effector proteins (SUMO, Ubc9 and PIAS) in Ae. aegypti cells led to enhanced levels of arbovirus replication from three different families; Zika (Flaviviridae), Semliki Forest (Togaviridae), and Bunyamwera (Bunyaviridae) viruses. Our findings identify an important role for mosquito SUMOylation in the cellular restriction of arboviruses that may directly influence vector competence and transmission of clinically important arboviruses.

Half the world’s population is at risk of infection from arboviruses transmitted by mosquitoes. Deciphering the viral-host interactions that influence the outcome of arbovirus infection in mosquitoes is beneficial to the development of future vector control strategies to limit arbovirus transmission and viral emergence within the human population. Similar to humans, mosquitoes possess different immune pathways to limit the replication of arboviruses. While the Small Ubiquitin-like MOdifier (SUMO) pathway is known to play an important role in the regulation of immune defences to viral infection in humans, the influence of this pathway during arbovirus infection in mosquito cells is currently unknown. Here we define the conservation, biochemical activity, and tissue distribution of the core effector proteins of the Aedes aegypti SUMOylation pathway. We show that the mosquito SUMOylation pathway plays a broadly antiviral role against a wide range of clinically important arboviruses, including Zika, Semliki Forest, and Bunyamwera viruses. Our findings identify SUMOylation as an important component of the antiviral response to arbovirus infection in mosquito cells.

Limited influence of the microbiome on the transcriptional profile of female Aedes aegypti mosquitoes

The microbiome is an assemblage of microorganisms living in association with a multicellular host. Numerous studies have identified a role for the microbiome in host physiology, development, immunity, and behaviour. The generation of axenic (germ-free) and gnotobiotic model systems has been vital to dissecting the role of the microbiome in host biology. We have previously reported the generation of axenic Aedes aegypti mosquitoes, the primary vector of several human pathogenic viruses, including dengue virus and Zika virus. In order to better understand the influence of the microbiome on mosquitoes, we examined the transcriptomes of axenic and conventionally reared Ae. aegypti before and after a blood meal. Our results suggest that the microbiome has a much lower effect on the mosquito's gene expression than previously thought with only 170 genes influenced by the axenic state, while in contrast, blood meal status influenced 809 genes. The pattern of expression influenced by the microbiome is consistent with transient changes similar to infection rather than sweeping physiological changes. While the microbiome does seem to affect some pathways such as immune function and metabolism, our data suggest the microbiome is primarily serving a nutritional role in development with only minor effects in the adult.

Yellow-g and Yellow-g2 proteins are required for egg desiccation resistance and temporal pigmentation in the Asian tiger mosquito, Aedes albopictus

Eggs from Aedes mosquitoes exhibit desiccation resistance that helps them to survive and spread as human disease vectors throughout the world. Previous studies have suggested that eggshell/chorion melanization and/or serosal cuticle formation are important for desiccation resistance. In this study, using dsRNAs for target genes, we analyzed the functional importance of two ovary-specific yellow genes, AalY-g and AalY-g2, in the resistance to egg desiccation of the Asian tiger mosquito, Aedes albopictus, a species in which neither the timing of the melanization nor temporal development of the serosal cuticle is correlated with desiccation resistance. Injections of dsAalY-g, dsAalY-g2 or dsAalY-g/g2 (co-injection) into adult females have no effect on their fecundity. However, initial melanization is delayed by 1-2 h with the eggshells eventually becoming black similar to that observed in eggs from dsEGFP-injected control females. In addition, the shape of the eggs from dsAalY-g, -g2 and -g/g2-treated females is abnormally crescent-shaped and the outermost exochorion is more fragile and partially peeled off. dsEGFP control eggs, like those from the wild-type strain, acquire resistance to desiccation between 18 and 24 h after oviposition (HAO). In contrast, ~80% of the 24 HAO dsAalY-g and dsAalY-g2 eggs collapse when they are transferred to a low humidity environment. In addition, there is no electron-dense outer endochorion evident in either dsAalY-g or dsAalY-g2 eggs. These results support the hypothesis that AalY-g and AalY-g2 regulate the timing of eggshell darkening and are required for integrity of the exochorion as well as for rigidity, normal morphology and formation of the outer endochorion, a structure that apparently is critical for desiccation resistance of the Ae. albopictus egg.

Aedes aegypti Immune Response and Its Potential Impact on Dengue Virus Transmission

Dengue virus (DENV) transmission to human populations requires infection of vector mosquitoes as an essential component of the transmission process. DENV transmission leads to infections that range from asymptomatic to life-threatening pathologies, such as dengue hemorrhagic fever and dengue shock syndrome. Aedes aegypti is the principal vector of DENV, and its vector competence consists of the intrinsic factors, genes, molecules, and pathways that allow infection, replication, and dissemination of this virus throughout the cells of mosquito tissues. In the search for mosquito molecular targets to block DENV transmission, the effect of DENV infection on mosquitoes has been an important focus of research. In this study, we review the findings of research on the effect of DENV infection on mosquito tissue cells and the immunity pathways and molecules that are involved in this infection. We emphasize the relevance of recent findings to understand the relationship between Ae. aegypti immune response, vector competence, and DENV transmission to human hosts.

Thermoprotective adaptations are critical for arthropods feeding on warm-blooded hosts

Blood feeding in arthropods has evolved in multiple lineages. This feeding preference provides a source of ample proteins and lipids for egg production and survival, but ingestion of a large warm blood-meal can boost the arthropod's body temperature 15°-20°C within seconds to minutes. This represents one of, if not the most, rapid thermal change documented under a natural setting. Here, we describe mechanisms of thermoregulation and thermotolerance in arthropods during blood feeding. The ability to prevent blood-induced thermal damage is a fundamental physiological adaptation linked to the use of warm-blooded vertebrates as food sources. Specific functional and comparative studies have identified unique and divergent mechanisms that suppress or repair thermal stress during blood feeding. These mechanisms include countercurrent heat exchange, evaporative cooling, and upregulation of stress associated proteins.

In vitro and in silico Models to Study Mosquito-Borne Flavivirus Neuropathogenesis, Prevention, and Treatment

Mosquito-borne flaviviruses can cause disease in the nervous system, resulting in a significant burden of morbidity and mortality. Disease models are necessary to understand neuropathogenesis and identify potential therapeutics and vaccines. Non-human primates have been used extensively but present major challenges. Advances have also been made toward the development of humanized mouse models, but these models still do not fully represent human pathophysiology. Recent developments in stem cell technology and cell culture techniques have allowed the development of more physiologically relevant human cell-based models. In silico modeling has also allowed researchers to identify and predict transmission patterns and discover potential vaccine and therapeutic candidates. This review summarizes the research on in vitro and in silico models used to study three mosquito-borne flaviviruses that cause neurological disease in humans: West Nile, Dengue, and Zika. We also propose a roadmap for 21st century research on mosquito-borne flavivirus neuropathogenesis, prevention, and treatment.

Using genetic variation in Aedes aegypti to identify candidate anti-dengue virus genes

BACKGROUND

Transcriptomic profiling has generated extensive lists of genes that respond to viral infection in mosquitoes. These gene lists contain two types of genes; (1) those that are responsible for the insect's natural antiviral defense mechanisms, including some known innate immunity genes, and (2) genes whose change in expression may occur simply as a result of infection. As genetic modification tools for mosquitoes continue to improve, the opportunities to make refractory insects via allelic replacement or delivery of small RNAs that alter gene expression are expanding. Therefore, the ability to identify which genes in transcriptional profiles may have immune function has increasing value. Arboviruses encounter a range of mosquito tissues and physiologies as they traverse from the midgut to the salivary glands. While the midgut is well-studied as the primary tissue barrier, antiviral genes expressed in the subsequent tissues of the carcass offer additional candidates for second stage intervention in the mosquito body.

METHODS

Mosquito lines collected recently from field populations exhibit natural genetic variation for dengue virus susceptibility. We sought to use a modified full-sib breeding design to identify mosquito families that varied in their dengue viral load in their bodies post infection.

RESULTS

By delivering virus intrathoracically, we bypassed the midgut and focused on whole body responses in order to evaluate carcass-associated refractoriness. We tested 25 candidate genes selected for their appearance in multiple published transcriptional profiles and were able to identify 12 whose expression varied with susceptibility in the genetic families.

CONCLUSIONS

This method, using natural genetic variation, offers a simple means to screen and reduce candidate gene lists prior to carrying out more labor-intensive functional studies. The extracted RNA from the females across the families represents a storable resource that can be used to screen subsequent candidate genes in the future. The aspect of vector competence being assessed could be varied by focusing on different tissues or time points post infection.

Homologs of Human Dengue-Resistance Genes, FKBP1B and ATCAY, Confer Antiviral Resistance in Aedes aegypti Mosquitoes

Dengue virus (DENV) is transmitted by mosquitoes and is a major public health concern. The study of innate mosquito defense mechanisms against DENV have revealed crucial roles for the Toll, Imd, JAK-STAT, and RNAi pathways in mediating DENV in the mosquito. Often overlooked in such studies is the role of intrinsic cellular defense mechanisms that we hypothesize to work in concert with the classical immune pathways to affect organismal defense. Our understanding of the molecular interaction of DENV with mosquito host cells is limited, and we propose to expand upon the recent results from a genome-scale, small interfering RNA (siRNA)-based study that identified mammalian host proteins associated with resistance to dengue/West Nile virus (DENV/WNV) infection. The study identified 22 human DENV/WNV resistance genes (DVR), and we hypothesized that a subset would be functionally conserved in Aedes aegypti mosquitoes, imparting cellular defense against flaviviruses in this species. We identified 12 homologs of 22 human DVR genes in the Ae. aegypti genome. To evaluate their possible role in cellular resistance/antiviral defense against DENV, we used siRNA silencing targeted against each of the 12 homologs in an Ae. aegypti cell line (Aag2) infected with DENV2 and identified that silencing of the two candidates, AeFKBP1 and AeATCAY, homologs of human FKBP1B and ATCAY, were associated with a viral increase. We then used dsRNA to silence each of the two genes in adult mosquitoes to validate the observed antiviral functions in vivo. Depletion of AeFKBP1 or AeATCAY increased viral dissemination through the mosquito at 14 days post-infection. Our results demonstrated that AeFKBP1 and AeATCAY mediate resistance to DENV akin to what has been described for their homologs in humans. AeFKBP1 and AeATCAY provide a rare opportunity to elucidate a DENV-resistance mechanism that may be evolutionarily conserved between humans and Ae. aegypti.

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.

Modulation of mitochondrial metabolism as a biochemical trait in blood feeding organisms: the redox vampire hypothesis redux

Hematophagous organisms undergo remarkable metabolic changes during the blood digestion process, increasing fermentative glucose metabolism, and reducing respiratory rates, both consequence of functional mitochondrial remodeling. Here, we review the pathways involved in energy metabolism and mitochondrial functionality in a comparative framework across different hematophagous species, and consider how these processes regulate redox homeostasis during blood digestion. The trend across distinct species indicate that a switch in energy metabolism might represent an important defensive mechanism to avoid the potential harmful interaction of oxidants generated from aerobic energy metabolism with products derived from blood digestion. Indeed, in insect vectors, blood feeding transiently reduces respiratory rates and oxidant production, irrespective of tissue and insect model. On the other hand, a different scenario is observed in several unrelated parasite species when exposed to blood digestion products, as respiratory rates reduce and mitochondrial oxidant production increase. The emerging picture indicates that re-wiring of energy metabolism, through reduced mitochondrial function, culminates in improved tolerance to redox insults and seems to represent a key step for hematophagous organisms to cope with the overwhelming and potentially toxic blood meal.

Extensive Genetic Differentiation between Homomorphic Sex Chromosomes in the Mosquito Vector, Aedes aegypti

Mechanisms and evolutionary dynamics of sex-determination systems are of particular interest in insect vectors of human pathogens like mosquitoes because novel control strategies aim to convert pathogen-transmitting females into nonbiting males, or rely on accurate sexing for the release of sterile males. In Aedes aegypti, the main vector of dengue and Zika viruses, sex determination is governed by a dominant male-determining locus, previously thought to reside within a small, nonrecombining, sex-determining region (SDR) of an otherwise homomorphic sex chromosome. Here, we provide evidence that sex chromosomes in Ae. aegypti are genetically differentiated between males and females over a region much larger than the SDR. Our linkage mapping intercrosses failed to detect recombination between X and Y chromosomes over a 123-Mbp region (40% of their physical length) containing the SDR. This region of reduced male recombination overlapped with a smaller 63-Mbp region (20% of the physical length of the sex chromosomes) displaying high male–female genetic differentiation in unrelated wild populations from Brazil and Australia and in a reference laboratory strain originating from Africa. In addition, the sex-differentiated genomic region was associated with a significant excess of male-to-female heterozygosity and contained a small cluster of loci consistent with Y-specific null alleles. We demonstrate that genetic differentiation between sex chromosomes is sufficient to assign individuals to their correct sex with high accuracy. We also show how data on allele frequency differences between sexes can be used to estimate linkage disequilibrium between loci and the sex-determining locus. Our discovery of large-scale genetic differentiation between sex chromosomes in Ae. aegypti lays a new foundation for mapping and population genomic studies, as well as for mosquito control strategies targeting the sex-determination pathway.

Meta-Analysis of Aedes aegypti Expression Datasets: Comparing Virus Infection and Blood-Fed Transcriptomes to Identify Markers of Virus Presence

The mosquito Aedes aegypti (L.) is vector of several arboviruses including dengue, yellow fever, chikungunya, and more recently zika. Previous transcriptomic studies have been performed to elucidate altered pathways in response to viral infection. However, the intrinsic coupling between alimentation and infection were unappreciated in these studies. Feeding is required for the initial mosquito contact with the virus and these events are highly dependent. Addressing this relationship, we reinterrogated datasets of virus-infected mosquitoes with two different diet schemes (fed and unfed mosquitoes), evaluating the metabolic cross-talk during both processes. We constructed coexpression networks with the differentially expressed genes of these comparison: virus-infected versus blood-fed mosquitoes and virus-infected versus unfed mosquitoes. Our analysis identified one module with 110 genes that correlated with infection status (representing ~0.7% of the A. aegypti genome). Furthermore, we performed a machine-learning approach and summarized the infection status using only four genes (AAEL012128, AAEL014210, AAEL002477, and AAEL005350). While three of the four genes were annotated as hypothetical proteins, AAEL012128 gene is a membrane amino acid transporter correlated with viral envelope binding. This gene alone is able to discriminate all infected samples and thus should have a key role to discriminate viral infection in the A. aegypti mosquito. Moreover, validation using external datasets found this gene as differentially expressed in four transcriptomic experiments. Therefore, these genes may serve as a proxy of viral infection in the mosquito and the others 106 identified genes provides a framework to future studies.

The midgut transcriptome of Aedes aegypti fed with saline or protein meals containing chikungunya virus reveals genes potentially involved in viral midgut escape

BACKGROUND

The mosquito Aedes aegypti is the primary vector for medically important arthropod-borne viruses, including chikungunya virus (CHIKV). Following oral acquisition, an arbovirus has to persistently infect several organs in the mosquito before becoming transmissible to another vertebrate host. A major obstacle an arbovirus has to overcome during its infection cycle inside the mosquito is the midgut escape barrier, representing the exit mechanism arboviruses utilize when disseminating from the midgut. To understand the transcriptomic basis of midgut escape and to reveal genes involved in the process, we conducted a comparative transcriptomic analysis of midgut samples from mosquitoes which had received a saline meal (SM) or a protein meal (PM) (not) containing CHIKV.

RESULTS

CHIKV which was orally acquired by a mosquito along with a SM or PM productively infected the midgut epithelium and disseminated to secondary tissues. A total of 27 RNA-Seq libraries from midguts of mosquitoes that had received PM or SM (not) containing CHIKV at 1 and 2 days post-feeding were generated and sequenced. Fewer than 80 genes responded differentially to the presence of CHIKV in midguts of mosquitoes that had acquired the virus along with SM or PM. SM feeding induced differential expression (DE) of 479 genes at day 1 and 314 genes at day 2 when compared to midguts of sugarfed mosquitoes. By comparison, PM feeding induced 6029 DE genes at day 1 and 7368 genes at day 2. Twenty-three DE genes encoding trypsins, metalloproteinases, and serine-type endopeptidases were significantly upregulated in midguts of mosquitoes at day 1 following SM or PM ingestion. Two of these genes were Ae. aegypti late trypsin (AeLT) and serine collagenase 1 precursor (AeSP1). In vitro, recombinant AeLT showed strong matrix metalloproteinase activity whereas recombinant AeSP1 did not.

CONCLUSIONS

By substituting a bloodmeal for SM, we identified midgut-expressed genes not involved in blood or protein digestion. These included genes coding for trypsins, metalloproteinases, and serine-type endopeptidases, which could be involved in facilitating midgut escape for arboviruses in Ae. aegypti. The presence of CHIKV in any of the ingested meals had relatively minor effects on the overall gene expression profiles in midguts.

RNA-Seq Comparison of Larval and Adult Malpighian Tubules of the Yellow Fever Mosquito Aedes aegypti Reveals Life Stage-Specific Changes in Renal Function

Introduction: The life history of Aedes aegypti presents diverse challenges to its diuretic system. During the larval and pupal life stages mosquitoes are aquatic. With the emergence of the adult they become terrestrial. This shifts the organism within minutes from an aquatic environment to a terrestrial environment where dehydration has to be avoided. In addition, female mosquitoes take large blood meals, which present an entirely new set of challenges to salt and water homeostasis.

Methods: To determine differences in gene expression associated with these different life stages, we performed an RNA-seq analysis of the main diuretic tissue in A. aegypti, the Malpighian tubules. We compared transcript abundance in 4th instar larvae to that of adult females and analyzed the data with a focus on transcripts that encode proteins potentially involved in diuresis, like water and solute channels as well as ion transporters. We compared our results against the model of potassium- and sodium chloride excretion in the Malpighian tubules proposed by Hine et al. (2014), which involves at least eight ion transporters and a proton-pump.

Results: We found 3,421 of a total number of 17,478 (19.6%) unique transcripts with a P < 0.05 and at least a 2.5 fold change in expression levels between the two groups. We identified two novel transporter genes that are highly expressed in the adult Malpighian tubules, which have not previously been part of the transport model in this species and may play important roles in diuresis. We also identified candidates for hypothesized sodium and chloride channels. Detoxification genes were generally higher expressed in larvae.

Significance: This study represents the first comparison of Malpighian tubule transcriptomes between larval and adult A. aegypti mosquitoes, highlighting key differences in their renal systems that arise as they transform from an aquatic filter-feeding larval stage to a terrestrial, blood-feeding adult stage.

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.

Retention of duplicated long-wavelength opsins in mosquito lineages by positive selection and differential expression

Background

Opsins are light sensitive receptors associated with visual processes. Insects typically possess opsins that are stimulated by ultraviolet, short and long wavelength (LW) radiation. Six putative LW-sensitive opsins predicted in the yellow fever mosquito, Aedes aegypti and malaria mosquito, Anopheles gambiae, and eight in the southern house mosquito, Culex quinquefasciatus, suggest gene expansion in the Family Culicidae (mosquitoes) relative to other insects. Here we report the first detailed molecular and evolutionary analyses of LW opsins in three mosquito vectors, with a goal to understanding the molecular basis of opsin-mediated visual processes that could be exploited for mosquito control.

Results

Time of divergence estimates suggest that the mosquito LW opsins originated from 18 or 19 duplication events between 166.9/197.5 to 1.07/0.94 million years ago (MY) and that these likely occurred following the predicted divergence of the lineages Anophelinae and Culicinae 145–226 MY. Fitmodel analyses identified nine amino acid residues in the LW opsins that may be under positive selection. Of these, eight amino acids occur in the N and C termini and are shared among all three species, and one residue in TMIII was unique to culicine species. Alignment of 5′ non-coding regions revealed potential Conserved Non-coding Sequences (CNS) and transcription factor binding sites (TFBS) in seven pairs of LW opsin paralogs.

Conclusions

Our analyses suggest opsin gene duplication and residues possibly associated with spectral tuning of LW-sensitive photoreceptors. We explore two mechanisms - positive selection and differential expression mediated by regulatory units in CNS – that may have contributed to the retention of LW opsin genes in Culicinae and Anophelinae. We discuss the evolution of mosquito LW opsins in the context of major Earth events and possible adaptation of mosquitoes to LW-dominated photo environments, and implications for mosquito control strategies based on disrupting vision-mediated behaviors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0910-6) contains supplementary material, which is available to authorized users.

DsRNA-mediated targeting of ribosomal transcripts RPS6 and RPL26 induces long-lasting and significant reductions in fecundity of the vector Aedes aegypti

Ribosomal transcripts produce critical proteins that are involved in most cellular production processes. Targeting ribosomal transcripts has produced mortality in mites and ticks but the effect of ribosomal transcript knockdown has not been thoroughly examined in mosquitoes. We examine the effects of triggers targeting four ribosomal proteins (RP) transcripts. Although no significant mortality was observed after dsRNA microinjection and subsequent blood feeding, significant contrasts were observed on fecundity. Triggers targeting RPS6 and RPL26 effectively reduced gene expression but more importantly, reduced reproductive output by more than 96% and 91% at the first oviposition while triggers targeting RPL1 and RPS2 did not cause a reduction although gene expression was reduced. Significantly reduced fecundity continued through a second oviposition cycle in dsRPS6 and dsRPL26 cohorts, although the effect was not as strong. Relative gene expression levels confirmed specific transcript knockdown up to 20days post-injection in mosquitoes that did not oviposit or produced reduced clutch sizes. Dissections at 36h post-blood meal indicated defects in oocyte provisioning. The strong phenotype produced by dsRPS6 allowed us to examine the effects in various tissues as well as the dose response, trigger format, delivery method and trigger specificity in Aedes aegypti. Strong knockdown was observed in the abdomen and the ovaries. Greater than 50ng of dsRPS6 significantly reduced fecundity but not when delivered in a sugar meal or as an siRNA. Similar bioassays with mutated dsRPS6 triggers indicates that up to three mismatches per possible siRNA are still effective in reducing fecundity. These studies indicate that while active and effective triggers can be developed for vector species, the lack of an efficient delivery method is the biggest barrier to use as a potential control method.

Human C5a Protein Participates in the Mosquito Immune Response Against Dengue Virus

Dengue virus (DENV) is transmitted by Aedes spp mosquitoes during a bloodmeal uptake. The bloodmeal consists of host cells, immune factors, and possibly blood-borne pathogens, such as arboviruses. Human cells and immune-related factors, like the complement system, can remain active in the bloodmeal and may be able to interact with pathogens in the mosquito. Previous studies have shown that active complement proteins impact Plasmodium parasite viability in the Anopheles midgut. Thus, we investigated the effects of the human complement on DENV infection in the midgut of Aedes aegypti. Our findings indicate that mosquitoes receiving DENV mixed with normal non-inactivated human serum showed significantly lower viremia than those fed with heat-inactivated serum. This implies that human complement may act to limit DENV infection in the mosquito midgut. In addition, we found that human complement C5a protein was able to directly communicate with mosquito cells, affecting the cell antiviral response against DENV. Our results also show that human C5a protein is able to interact with several membrane-bound mosquito proteins. Together these results suggest an important role of human complement protein in DENV transmission.

Expression Profiles and RNAi Silencing of Inhibitor of Apoptosis Transcripts in Aedes, Anopheles, and Culex Mosquitoes (Diptera: Culicidae)

Effective mosquito control is vital to curtail the devastating health effects of many vectored diseases. RNA interference (RNAi)-mediated control of mosquitoes is an attractive alternative to conventional chemical pesticides. Previous studies have suggested that transcripts for inhibitors of apoptosis (IAPs) may be good RNAi targets. To revisit and extend previous reports, we examined the expression of Aedes aegypti (L.) IAPs (AaeIAPs) 1, 2, 5, 6, 9, and a viral IAP-associated factor (vIAF) as well as Anopheles quadrimaculatus Say and Culex quinquefasciatus Say IAP1 homologs (AquIAP1 and CquIAP1) in adult females. Expression profiles of IAPs suggested that some older female mosquitoes had significantly higher IAP mRNA levels when compared to the youngest ones. Minor differences in expression of AaeIAPs were observed in mosquitoes that imbibed a bloodmeal, but the majority of the time points (up to 48 h) were not significantly different. Although in vitro experiments with the Ae. aegypti Aag-2 cell line demonstrated that the various AaeIAPs could be effectively knocked down within one day after dsRNA treatment, only Aag-2 cells treated with dsIAP1 displayed apoptotic morphology. Gene silencing and mortality were also evaluated after topical application and microinjection of the same dsRNAs into female Ae. aegypti. In contrast to previous reports, topical administration of dsRNA against AaeIAP1 did not yield a significant reduction in gene expression or increased mortality. Knockdown of IAP1 and other IAPs by microinjection did not result in significant mortality. In toto, our findings suggest that IAPs may not be suitable RNAi targets for controlling adult mosquito populations.

Arbovirus-mosquito interactions: RNAi pathway

Arthropod-borne (arbo) viruses infect hematophagous arthropods (vectors) to maintain virus transmission between vertebrate hosts. The mosquito vector actively controls arbovirus infection to minimize its fitness costs. The RNA interference (RNAi) pathway is the major antiviral response vectors use to restrict arbovirus infections. We know this because depleting RNAi gene products profoundly impacts arbovirus replication, the antiviral RNAi pathway genes undergo positive, diversifying selection and arboviruses have evolved strategies to evade the vector's RNAi responses. The vector's RNAi defense and arbovirus countermeasures lead to an arms race that prevents potential virus-induced fitness costs yet maintains arbovirus infections needed for transmission. This review will discuss the latest findings in RNAi-arbovirus interactions in the model insect (Drosophila melanogaster) and in specific mosquito vectors.

Contrasting genetic structure between mitochondrial and nuclear markers in the dengue fever mosquito from Rio de Janeiro: implications for vector control

Dengue is the most prevalent global arboviral disease that affects over 300 million people every year. Brazil has the highest number of dengue cases in the world, with the most severe epidemics in the city of Rio de Janeiro (Rio). The effective control of dengue is critically dependent on the knowledge of population genetic structuring in the primary dengue vector, the mosquito Aedes aegypti. We analyzed mitochondrial and nuclear genomewide single nucleotide polymorphism markers generated via Restriction-site Associated DNA sequencing, as well as traditional microsatellite markers in Ae. aegypti from Rio. We found four divergent mitochondrial lineages and a strong spatial structuring of mitochondrial variation, in contrast to the overall nuclear homogeneity across Rio. Despite a low overall differentiation in the nuclear genome, we detected strong spatial structure for variation in over 20 genes that have a significantly altered expression in response to insecticides, xenobiotics, and pathogens, including the novel biocontrol agent Wolbachia. Our results indicate that high genetic diversity, spatially unconstrained admixing likely mediated by male dispersal, along with locally heterogeneous genetic variation that could affect insecticide resistance and mosquito vectorial capacity, set limits to the effectiveness of measures to control dengue fever in Rio.

Mitochondrial physiology in the major arbovirus vector Aedes aegypti: substrate preferences and sexual differences define respiratory capacity and superoxide production

Adult females of Aedes aegypti are facultative blood sucking insects and vectors of Dengue and yellow fever viruses. Insect dispersal plays a central role in disease transmission and the extremely high energy demand posed by flight is accomplished by a very efficient oxidative phosphorylation process, which take place within flight muscle mitochondria. These organelles play a central role in energy metabolism, interconnecting nutrient oxidation to ATP synthesis, but also represent an important site of cellular superoxide production. Given the importance of mitochondria to cell physiology, and the potential contributions of this organelle for A. aegypti biology and vectorial capacity, here, we conducted a systematic assessment of mitochondrial physiology in flight muscle of young adult A. aegypti fed exclusively with sugar. This was carried out by determining the activities of mitochondrial enzymes, the substrate preferences to sustain respiration, the mitochondrial bioenergetic efficiency and capacity, in both mitochondria-enriched preparations and mechanically permeabilized flight muscle in both sexes. We also determined the substrates preferences to promote mitochondrial superoxide generation and the main sites where it is produced within this organelle. We observed that respiration in A. aegypti mitochondria was essentially driven by complex I and glycerol 3 phosphate dehydrogenase substrates, which promoted distinct mitochondrial bioenergetic capacities, but with preserved efficiencies. Respiration mediated by proline oxidation in female mitochondria was strikingly higher than in males. Mitochondrial superoxide production was essentially mediated through proline and glycerol 3 phosphate oxidation, which took place at sites other than complex I. Finally, differences in mitochondrial superoxide production among sexes were only observed in male oxidizing glycerol 3 phosphate, exhibiting higher rates than in female. Together, these data represent a significant step towards the understanding of fundamental mitochondrial processes in A. aegypti, with potential implications for its physiology and vectorial capacity.

Exome and transcriptome sequencing of Aedes aegypti identifies a locus that confers resistance to Brugia malayi and alters the immune response

Many mosquito species are naturally polymorphic for their abilities to transmit parasites, a feature which is of great interest for controlling vector-borne disease. Aedes aegypti, the primary vector of dengue and yellow fever and a laboratory model for studying lymphatic filariasis, is genetically variable for its capacity to harbor the filarial nematode Brugia malayi. The genome of Ae. aegypti is large and repetitive, making genome resequencing difficult and expensive. We designed exome captures to target protein-coding regions of the genome, and used association mapping in a wild Kenyan population to identify a single, dominant, sex-linked locus underlying resistance. This falls in a region of the genome where a resistance locus was previously mapped in a line established in 1936, suggesting that this polymorphism has been maintained in the wild for the at least 80 years. We then crossed resistant and susceptible mosquitoes to place both alleles of the gene into a common genetic background, and used RNA-seq to measure the effect of this locus on gene expression. We found evidence for Toll, IMD, and JAK-STAT pathway activity in response to early stages of B. malayi infection when the parasites are beginning to die in the resistant genotype. We also found that resistant mosquitoes express anti-microbial peptides at the time of parasite-killing, and that this expression is suppressed in susceptible mosquitoes. Together, we have found that a single resistance locus leads to a higher immune response in resistant mosquitoes, and we identify genes in this region that may be responsible for this trait.

Two cathepsins B are responsible for the yolk protein hydrolysis in Culex quinquefasciatus

Despite the established role of Culex quinquefasciatus as a vector of various neurotropic viruses, such as the Rift Valley and West Nile viruses, as well as lymphatic filariasis, little is known regarding the organism’s reproductive physiology. As in other oviparous animals, vitellogenin, the most important source of nutrients for the embryo development, is digested by intracellular proteases. Using mass spectrometry, we have identified two cathepsin B homologues partially purified by self-proteolysis of Cx. quinquefasciatus total egg extract. The transcriptional profile of these two cathepsin B homologues was determined by quantitative RT-PCR, and the enzymatic activity associated with the peptidase was determined in ovaries after female engorgement. According to the VectorBase (vectorbase.org) annotation, both cathepsin B homologues shared approximately 66% identity in their amino acid sequences. The two cathepsin B genes are expressed simultaneously in the fat body of the vitellogenic females, and enzymatic activity was detected within the ovaries, suggesting an extra-ovarian origin. Similar to the transcriptional profile of vitellogenin, cathepsin B transcripts were shown to accumulate post-blood meal and reached their highest expression at 36 h PBM. However, while vitellogenin expression decreased drastically at 48 h PBM, the expression of the cathepsins increased until 84 h PBM, at which time the females of our colony were ready for oviposition. The similarity between their transcriptional profiles strongly suggests a role for the cathepsin B homologues in vitellin degradation.

Anopheles gambiae blood feeding initiates an anticipatory defense response to Plasmodium berghei

Mosquitoes have potent innate defense mechanisms that protect them from infection by diverse pathogens. Much remains unknown about how different pathogens are sensed and specific responses triggered. Leucine-Rich repeat IMmune proteins (LRIMs) are a mosquito-specific family of putative innate receptors. Although some LRIMs have been implicated in mosquito immune responses, the function of most family members is largely unknown. We screened Anopheles gambiae LRIMs by RNAi for effects on mosquito infection by rodent malaria and found that LRIM9 is a Plasmodium berghei antagonist with phenotypes distinct from family members LRIM1 and APL1C, which are key components of the mosquito complement-like pathway. LRIM9 transcript and protein levels are significantly increased after blood feeding but are unaffected by Plasmodium or midgut microbiota. Interestingly, LRIM9 in the hemolymph is strongly upregulated by direct injection of the ecdysteroid, 20-hydroxyecdysone. Our data suggest that LRIM9 may define a novel anti-Plasmodium immune defense mechanism triggered by blood feeding and that hormonal changes may alert the mosquito to bolster its defenses in anticipation of exposure to blood-borne pathogens.

Transcriptome sequencing and developmental regulation of gene expression in Anopheles aquasalis

BACKGROUND

Anopheles aquasalis is a major malaria vector in coastal areas of South and Central America where it breeds preferentially in brackish water. This species is very susceptible to Plasmodium vivax and it has been already incriminated as responsible vector in malaria outbreaks. There has been no high-throughput investigation into the sequencing of An. aquasalis genes, transcripts and proteins despite its epidemiological relevance. Here we describe the sequencing, assembly and annotation of the An. aquasalis transcriptome.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 419 thousand cDNA sequence reads, encompassing 164 million nucleotides, were assembled in 7544 contigs of ≥ 2 sequences, and 1999 singletons. The majority of the An. aquasalis transcripts encode proteins with their closest counterparts in another neotropical malaria vector, An. darlingi. Several analyses in different protein databases were used to annotate and predict the putative functions of the deduced An. aquasalis proteins. Larval and adult-specific transcripts were represented by 121 and 424 contig sequences, respectively. Fifty-one transcripts were only detected in blood-fed females. The data also reveal a list of transcripts up- or down-regulated in adult females after a blood meal. Transcripts associated with immunity, signaling networks and blood feeding and digestion are discussed.

CONCLUSIONS/SIGNIFICANCE

This study represents the first large-scale effort to sequence the transcriptome of An. aquasalis. It provides valuable information that will facilitate studies on the biology of this species and may lead to novel strategies to reduce malaria transmission on the South American continent. The An. aquasalis transcriptome is accessible at http://exon.niaid.nih.gov/transcriptome/An_aquasalis/Anaquexcel.xlsx.

Identification of life-stage and tissue-specific splice variants of an inward rectifying potassium (Kir) channel in the yellow fever mosquito Aedes aegypti

Inward-rectifier potassium (Kir) channels play key roles in nerve, muscle, and epithelial cells in mammals, but their physiological roles in insects remain to be determined. The yellow fever mosquito (Aedes aegypti) possesses five different genes encoding Kir channel subunits: Kir1, Kir2A, Kir2B, Kir2B', and Kir3. We have recently cloned and characterized the Kir1, Kir2B, and Kir3 cDNAs in the renal (Malpighian) tubules of adult female Ae. aegypti. Here we characterize the expression of the Kir2A gene in Ae. aegypti, which was not abundantly expressed in Malpighian tubules. We find that the 1) Kir2A gene is expressed primarily in the midgut and hindgut of adult female mosquitoes, and 2) Kir2A mRNAs are alternatively spliced into three distinct variants (Kir2A-a, -b, and -c). The deduced Kir2A proteins from these splice forms share a completely conserved transmembrane domain (a pore-forming domain flanked by two transmembrane-spanning segments), but possess novel NH2-terminal and/or COOH-terminal domains. Semi-quantitative RT-PCR analyses indicate that the splice variants exhibit both developmental- and tissue-specific expression. Lastly, we provide insights into the conservation of alternative splicing among the Kir2A genes of dipterans, which may add molecular diversity that compensates for the relatively limited number of Kir channel genes in insects compared to mammals.

Integrated proteomic and transcriptomic analysis of the Aedes aegypti eggshell

Background

Mosquito eggshells show remarkable diversity in physical properties and structure consistent with adaptations to the wide variety of environments exploited by these insects. We applied proteomic, transcriptomic, and hybridization in situ techniques to identify gene products and pathways that participate in the assembly of the Aedes aegypti eggshell. Aedes aegypti population density is low during cold and dry seasons and increases immediately after rainfall. The survival of embryos through unfavorable periods is a key factor in the persistence of their populations. The work described here supports integrated vector control approaches that target eggshell formation and result in Ae. aegypti drought-intolerant phenotypes for public health initiatives directed to reduce mosquito-borne diseases.

Results

A total of 130 proteins were identified from the combined mass spectrometric analyses of eggshell preparations.

Conclusions

Classification of proteins according to their known and putative functions revealed the complexity of the eggshell structure. Three novel Ae. aegypti vitelline membrane proteins were discovered. Odorant-binding and cysteine-rich proteins that may be structural components of the eggshell were identified. Enzymes with peroxidase, laccase and phenoloxidase activities also were identified, and their likely involvements in cross-linking reactions that stabilize the eggshell structure are discussed.

Targeted genome editing in Aedes aegypti using TALENs

The Culicine mosquito, Aedes aegypti, is both a major vector of arthropod-borne viruses (arboviruses) and a genetic model organism for arbovirus transmission. TALE nucleases (TALENs), a group of artificial enzymes capable of generating site-specific DNA lesions, consist of a non-specific FokI endonuclease cleavage domain fused to an engineered DNA binding domain specific to a target site. While TALENs have become an important tool for targeted gene disruption in a variety of organisms, application to the mosquito genome is a new approach. We recently described the use of TALENs to perform heritable genetic disruptions in A. aegypti. Here, we provide detailed methods that will allow other research laboratories to capitalize on the potential of this technology for understanding mosquito gene function. We describe target site selection, transient embryo-based assays to rapidly assess TALEN activity, embryonic microinjection and downstream screening steps to identify target site mutations.

Changes in transcript abundance for cuticular proteins and other genes three hours after a blood meal in Anopheles gambiae

Numerous studies have examined changes in transcript levels after Anopheles gambiae takes a blood meal. Marinotti et al. (2006) used microarrays and reported massive changes in transcript levels 3 h after feeding (BF3h) compared to non-blood fed (NBF). We were intrigued by the number of transcripts for structural cuticular proteins (CPs) that showed such major differences in levels and employed paired-end (50 bp) RNA-seq technology to compare whole body transcriptomes from 5-day-old females NBF and BF3h. We detected transcripts for the majority of CPs (164/243) but levels of only 12 were significantly altered by the blood meal. While relative transcript levels of NBF females were somewhat similar to the microarray data, there were major differences in BF3h animals, resulting in levels of many transcripts, both for CPs and other genes changing in the opposite direction. We compared our data also to other studies done with both microarrays and RNA-seq. Findings were consistent that a small number of CP genes have transcripts that persist even in 5-day-old adults. Some of these transcripts showed diurnal rhythms (Rund et al., 2013; Rinker et al., 2013). In situ hybridization revealed that transcripts for several of these CP genes were found exclusively or predominantly in the eye. Transcripts other than for CPs that changed in response to blood-feeding were predominantly expressed in midgut and Malpighian tubules. Even in these tissues, genes responsible for proteins with similar functions, such as immunity or digestion, responded differently, with transcript levels for some rising and others falling. These data demonstrate that genes coding for some CPs are dynamic in expression even in adults and that the response to a blood meal is rapid and precisely orchestrated.

Expression and accumulation of the two-domain odorant-binding protein AaegOBP45 in the ovaries of blood-fed Aedes aegypti

Background

Aedes aegypti mosquitoes are the main vectors of dengue viruses. Despite global efforts to reduce the prevalence of dengue using integrated vector management strategies, innovative alternatives are necessary to help prevent virus transmission. Detailed characterizations of Ae. aegypti genes and their products provide information about the biology of mosquitoes and may serve as foundations for the design of new vector control methods.

Findings

We studied the Ae. aegypti gene, AAEL010714, that encodes a two-domain odorant-binding protein, AaegOBP45. The predicted gene structure and sequence were validated, although single nucleotide polymorphisms were observed. Transcriptional and translational products accumulate in the ovaries of blood fed females and are not detected or are at low abundance in other tissues.

Conclusions

We validated the Ae. aegypti AAEL010714 gene sequence and characterized the expression profile of a two-domain OBP expressed in ovaries. We propose that AaegOBP45 function as a component of the mosquito eggshell.

The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector

Mosquitoes are vectors of a number of important human and animal diseases. The development of novel vector control strategies requires a thorough understanding of mosquito biology. To facilitate this, we used RNA-seq to identify novel genes and provide the first high-resolution view of the transcriptome throughout development and in response to blood feeding in a mosquito vector of human disease, Aedes aegypti, the primary vector for Dengue and yellow fever. We characterized mRNA expression at 34 distinct time points throughout Aedes development, including adult somatic and germline tissues, by using polyA+ RNA-seq. We identify a total of 14,238 novel new transcribed regions corresponding to 12,597 new loci, as well as many novel transcript isoforms of previously annotated genes. Altogether these results increase the annotated fraction of the transcribed genome into long polyA+ RNAs by more than twofold. We also identified a number of patterns of shared gene expression, as well as genes and/or exons expressed sex-specifically or sex-differentially. Expression profiles of small RNAs in ovaries, early embryos, testes, and adult male and female somatic tissues also were determined, resulting in the identification of 38 new Aedes-specific miRNAs, and ~291,000 small RNA new transcribed regions, many of which are likely to be endogenous small-interfering RNAs and Piwi-interacting RNAs. Genes of potential interest for transgene-based vector control strategies also are highlighted. Our data have been incorporated into a user-friendly genome browser located at www.Aedes.caltech.edu, with relevant links to Vectorbase (www.vectorbase.org)

Stage-specific transcription during development of Aedes aegypti

BACKGROUND

Aedes aegypti is the most important global vector of dengue virus infection in humans. Availability of the draft genome sequence of this mosquito provides unique opportunities to study different aspects of its biology, including identification of genes and pathways relevant to the developmental processes associated with transition across individual life stages. However, detailed knowledge of gene expression patterns pertaining to developmental stages of A. aegypti is largely lacking.

RESULTS

We performed custom cDNA microarray analyses to examine the expression patterns among six developmental stages: early larvae, late larvae, early pupae, late pupae, and adult male and female mosquitoes. Results revealed 1,551 differentially expressed transcripts (DETs) showing significant differences in levels of expression between these life stages. The data suggests that most of the differential expression occurs in a stage specific manner in A. aegypti. Based on hierarchical clustering of expression levels, correlated expression patterns of DETs were also observed among developmental stages. Weighted gene correlation network analysis revealed modular patterns of expression among the DETs. We observed that hydrolase activity, membrane, integral to membrane, DNA binding, translation, ribosome, nucleoside-triphosphatase activity, structural constituent of ribosome, ribonucleoprotein complex and receptor activity were among the top ten ranked GO (Gene Ontology) terms associated with DETs. Significant associations of DETs were also observed with specific KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway modules. Finally, comparisons with the previously reported developmental transcriptome of the malaria vector, Anopheles gambiae, indicated that gene expression patterns during developmental processes reflect both species-specific as well as common components of the two mosquito species.

CONCLUSIONS

Our study shows that genes involved in the developmental life cycle of A. aegypti are expressed in a highly stage-specific manner. This suggests that transcriptional events associated with transition through larval, pupal and adult stages are largely discrete.

miRNA genes of an invasive vector mosquito, Aedes albopictus

Aedes albopictus, a vector of Dengue and Chikungunya viruses, is a robust invasive species in both tropical and temperate environments. MicroRNAs (miRNAs) regulate gene expression and biological processes including embryonic development, innate immunity and infection. While a number of miRNAs have been discovered in some mosquitoes, no comprehensive effort has been made to characterize them from different developmental stages from a single species. Systematic analysis of miRNAs in Ae. albopictus will improve our understanding of its basic biology and inform novel strategies to prevent virus transmission. Between 10–14 million Illumina sequencing reads per sample were obtained from embryos, larvae, pupae, adult males, sugar-fed and blood-fed adult females. A total of 119 miRNA genes represented by 215 miRNA or miRNA star (miRNA*) sequences were identified, 15 of which are novel. Eleven, two, and two of the newly-discovered miRNA genes appear specific to Aedes, Culicinae, and Culicidae, respectively. A number of miRNAs accumulate predominantly in one or two developmental stages and the large number that showed differences in abundance following a blood meal likely are important in blood-induced mosquito biology. Gene Ontology (GO) analysis of the targets of all Ae. albopictus miRNAs provides a useful starting point for the study of their functions in mosquitoes. This study is the first systematic analysis of miRNAs based on deep-sequencing of small RNA samples of all developmental stages of a mosquito species. A number of miRNAs are related to specific physiological states, most notably, pre- and post-blood feeding. The distribution of lineage-specific miRNAs is consistent with mosquito phylogeny and the presence of a number of Aedes-specific miRNAs likely reflects the divergence between the Aedes and Culex genera.

Overlapping genes of Aedes aegypti: evolutionary implications from comparison with orthologs of Anopheles gambiae and other insects

Background

Although gene overlapping is a common feature of prokaryote and mitochondria genomes, such genes have also been identified in many eukaryotes. The overlapping genes in eukaryotes are extensively rearranged even between closely related species. In this study, we investigated retention and rearrangement of positionally overlapping genes between the mosquitoes Aedes aegypti (dengue virus vector) and Anopheles gambiae (malaria vector). The overlapping gene pairs of A. aegypti were further compared with orthologs of other selected insects to conduct several hypothesis driven investigations relating to the evolution and rearrangement of overlapping genes.

Results

The results show that as much as ~10% of the predicted genes of A. aegypti and A. gambiae are localized in positional overlapping manner. Furthermore, the study shows that differential abundance of introns and simple sequence repeats have significant association with positional rearrangement of overlapping genes between the two species. Gene expression analysis further suggests that antisense transcripts generated from the oppositely oriented overlapping genes are differentially regulated and may have important regulatory functions in these mosquitoes. Our data further shows that synonymous and non-synonymous mutations have differential but non-significant effect on overlapping localization of orthologous genes in other insect genomes.

Conclusion

Gene overlapping in insects may be a species-specific evolutionary process as evident from non-dependency of gene overlapping with species phylogeny. Based on the results, our study suggests that overlapping genes may have played an important role in genome evolution of insects.

Comparative genomics of odorant binding proteins in Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus

About 1 million people in the world die each year from diseases spread by mosquitoes, and understanding the mechanism of host identification by the mosquitoes through olfaction is at stake. The role of odorant binding proteins (OBPs) in the primary molecular events of olfaction in mosquitoes is becoming an important focus of biological research in this area. Here, we present a comprehensive comparative genomics study of OBPs in the three disease-transmitting mosquito species Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus starting with the identification of 110 new OBPs in these three genomes. We have characterized their genomic distribution and orthologous and phylogenetic relationships. The diversity and expansion observed with respect to the Aedes and Culex genomes suggests that the OBP gene family acquired functional diversity concurrently with functional constraints posed on these two species. Sequences with unique features have been characterized such as the "two-domain OBPs" (previously known as Atypical OBPs) and "MinusC OBPs" in mosquito genomes. The extensive comparative genomics featured in this work hence provides useful primary insights into the role of OBPs in the molecular adaptations of mosquito olfactory system and could provide more clues for the identification of potential targets for insect repellants and attractants.

The orthologue of the fruitfly sex behaviour gene fruitless in the mosquito Aedes aegypti: evolution of genomic organisation and alternative splicing

In Drosophila melanogaster the doublesex (dsx) and fruitless (fru) regulatory genes act at the bottom of the somatic sex determination pathway. Both are regulated via alternative splicing by an upstream female-specific TRA/TRA-2 complex, recognizing a common cis element. dsx controls somatic sexual differentiation of non-neural as well as of neural tissues. fru, on the other hand, expresses male-specific functions only in neural system where it is required to built the neural circuits underlying proper courtship behaviour. In the mosquito Aedes aegypti sex determination is different from Drosophila. The key male determiner M, which is located on one of a pair of homomorphic sex chromosomes, controls sex-specific splicing of the mosquito dsx orthologue. In this study we report the genomic organization and expression of the fru homologue in Ae. aegypti (Aeafru). We found that it is sex-specifically spliced suggesting that it is also under the control of the sex determination pathway. Comparative analyses between the Aeafru and Anopheles gambiae fru (Angfru) genomic loci revealed partial conservation of exon organization and extensive divergence of intron lengths. We find that Aeadsx and Aeafru share novel cis splicing regulatory elements conserved in the alternatively spliced regions. We propose that in Aedes aegypti sex-specific splicing of dsx and fru is most likely under the control of splicing regulatory factors which are different from TRA and TRA-2 found in other dipteran insects and discuss the potential use of fru and dsx for developing new genetic strategies in vector control.

The roles of serpins in mosquito immunology and physiology

In vector-borne diseases, the complex interplay between pathogen and its vector's immune system determines the outcome of infection and therefore disease transmission. Serpins have been shown in many animals to be key regulators of innate immune reactions. Their control over regulatory proteolytic cascades ultimately decides whether the recognition of a pathogen will lead to an appropriate immune response. In mosquitoes, serpins (SRPNs) regulate the activation of prophenoloxidase and thus melanization, contribute to malaria parasite lysis, and likely Toll pathway activation. Additionally, in culicine mosquitoes, SRPNs are able to regulate hemostasis in the vertebrate host, suggesting a crucial role during bloodfeeding. This review summarizes the annotation, transcriptional regulation, and current knowledge of SRPN function in the three mosquito species for which the complete genome sequence is available. Additionally, we give a brief overview of how SRPNs may be used to prevent transmission of vector-borne diseases.

SLC7 amino acid transporters of the yellow fever mosquito Aedes aegypti and their role in fat body TOR signaling and reproduction

BACKGROUND

An important function of the fat body in adult female mosquitoes is the conversion of blood meal derived amino acids (AA) into massive amounts of yolk protein precursors. A highly efficient transport mechanism for AAs across the plasma membrane of the fat body trophocytes is essential in order to deliver building blocks for the rapid synthesis of large amounts of these proteins. This mechanism consists in part of AA transporter proteins from the solute carrier family. These transporters have dual function; they function as transporters and participate in the nutrient signal transduction pathway that is activated in the fat body after a blood meal. In this study we focused on the solute carrier 7 family (SLC7), a family of AA transporters present in all metazoans that includes members with strong substrate specificity for cationic AAs.

METHODOLOGY/PRINCIPAL FINDINGS

We identified 11 putative SLC7 transporters in the genome sequence of Aedes aegypti. Phylogenetic analysis puts five of these in the cationic AA transporter subfamily (CAT) and six in the heterodimeric AA transporter (HAT) subfamily. All 11 A. aegypti SLC7 genes are expressed in adult females. Expression profiles are dynamic after a blood meal. We knocked down six fat body-expressed SLC7 transporters using RNAi and found that these 'knockdowns' reduced AA-induced TOR signaling. We also determined the effect these knockdowns had on the number of eggs deposited following a blood meal.

CONCLUSIONS/SIGNIFICANCE

Our analysis stresses the importance of SLC7 transporters in TOR signaling pathway and mosquito reproduction.

Strain Variation in the Transcriptome of the Dengue Fever Vector, Aedes aegypti

Studies of transcriptome dynamics provide a basis for understanding functional elements of the genome and the complexity of gene regulation. The dengue vector mosquito, Aedes aegypti, exhibits great adaptability to diverse ecological conditions, is phenotypically polymorphic, and shows variation in vectorial capacity to arboviruses. Previous genome sequencing showed richness in repetitive DNA and transposable elements that can contribute to genome plasticity. Population genetic studies revealed a varying degree of worldwide genetic polymorphism. However, the extent of functional genetic polymorphism across strains is unknown. The transcriptomes of three Ae. aegypti strains, Chetumal (CTM), Rexville D-Puerto Rico (Rex-D) and Liverpool (LVP), were compared. CTM is more susceptible than Rex- D to infection by dengue virus serotype 2. A total of 4188 transcripts exhibit either no or small variation (<2-fold) among sugar-fed samples of the three strains and between sugar- and blood-fed samples within each strain, corresponding most likely to genes encoding products necessary for vital functions. Transcripts enriched in blood-fed mosquitoes encode proteins associated with catalytic activities, molecular transport, metabolism of lipids, carbohydrates and amino acids, and functions related to blood digestion and the progression of the gonotropic cycle. Significant qualitative and quantitative differences were found in individual transcripts among strains including differential representation of paralogous gene products. The majority of immunity-associated transcripts decreased in accumulation after a bloodmeal and the results are discussed in relation to the different susceptibility of CTM and Rex-D mosquitoes to DENV2 infection.

The fat body transcriptomes of the yellow fever mosquito Aedes aegypti, pre- and post- blood meal

BACKGROUND

The fat body is the main organ of intermediary metabolism in insects and the principal source of hemolymph proteins. As part of our ongoing efforts to understand mosquito fat body physiology and to identify novel targets for insect control, we have conducted a transcriptome analysis of the fat body of Aedes aegypti before and in response to blood feeding.

RESULTS

We created two fat body non-normalized EST libraries, one from mosquito fat bodies non-blood fed (NBF) and another from mosquitoes 24 hrs post-blood meal (PBM). 454 pyrosequencing of the non-normalized libraries resulted in 204,578 useable reads from the NBF sample and 323,474 useable reads from the PBM sample. Alignment of reads to the existing reference Ae. aegypti transcript libraries for analysis of differential expression between NBF and PBM samples revealed 116,912 and 115,051 matches, respectively. De novo assembly of the reads from the NBF sample resulted in 15,456 contigs, and assembly of the reads from the PBM sample resulted in 15,010 contigs. Collectively, 123 novel transcripts were identified within these contigs. Prominently expressed transcripts in the NBF fat body library were represented by transcripts encoding ribosomal proteins. Thirty-five point four percent of all reads in the PBM library were represented by transcripts that encode yolk proteins. The most highly expressed were transcripts encoding members of the cathepsin b, vitellogenin, vitellogenic carboxypeptidase, and vitelline membrane protein families.

CONCLUSION

The two fat body transcriptomes were considerably different from each other in terms of transcript expression in terms of abundances of transcripts and genes expressed. They reflect the physiological shift of the pre-feeding fat body from a resting state to vitellogenic gene expression after feeding.