Life cycle transcriptome of the malaria mosquito Anopheles gambiae and comparison with the fruitfly Drosophila melanogaster

Chemical Probes to Interrogate the Extreme Environment of Mosquito Larval Guts

Mosquito control methods are vital to curtail the spread of life-threatening illnesses, such as dengue fever, malaria, and yellow fever. Vector control technologies must be selective to minimize deleterious effects on our ecosystem. Successful methods that control mosquito larva populations utilize the uniquely high alkaline nature of the midgut. Here, we present novel protected triazabutadienes (pTBD) that are deprotected under basic conditions of the larval midgut, releasing an aryl diazonium ion (ADI) that results in protein modification. The probes contain a bioorthogonal terminal alkyne handle, enabling a selective Cu-click reaction with an azidofluorophore for quantification by SDS PAGE and visualization using fluorescence microscopy. A control TBD, unable to release an ADI, did not label the midgut. We envision our chemical probes will aid in the development of new selective mosquito control methods, thus preventing the spread of mosquito-borne illnesses with minimal impact on other organisms in the ecosystem.

Chemical probes to interrogate the extreme environment of mosquito larval guts

Mosquito control methods are vital for the spread of life-threatening illnesses such as dengue fever, malaria, and yellow fever. Vector control technologies must be selective to minimize deleterious effects to our ecosystem. Successful methods that control mosquito larva populations utilize the uniquely high alkaline nature of the midgut. Here, we present novel protected triazabutadienes (pTBD) which are deprotected under basic conditions of the larval midgut, releasing an aryl diazonium ion (ADI) that results in protein modification. The probes contain a bioorthogonal terminal alkyne handle, enabling a selective Cu-click reaction with an azido-fluorophore for quantification by SDS PAGE and visualization using fluorescence microscopy. A control TBD, unable to release an ADI, did not label the midgut. We envision our chemical probes will aid in the development of new selective mosquito control methods thus preventing the spread of mosquito-borne illnesses with minimal impact on other organisms in the ecosystem.

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.

Pipiserpin, a Culex factor Xa inhibitor, affects female reproductive capacity and serves as a potential target for mosquito control

BACKGROUND

Mosquito control is still the main prevention and control measure for numerous mosquito-borne diseases causing millions of deaths each year. New strategies for mosquito control are in demand. Proteases play an important role in mosquito physiology, therefore this study explored the inhibition of a serpin (serine protease inhibitor) in mosquitoes and its effect on reproductive capacity.

RESULTS

A factor Xa inhibitor homolog (named Pipiserpin) was amplified and identified in Culex pipiens pallens mosquitoes. We expressed a recombinant Pipiserpin protein in vitro against which a mouse antiserum was generated. We found that female mosquitoes expressed more Pipiserpin protein than male mosquitoes. After mating, female mosquitoes were fed with blood mixed with different amounts of antisera and results showed that consumption of Pipiserpin impeded ovary development and decreased eggs hatching rates compared to that of the pre-immune serum group.

CONCLUSION

We identified a Culex mosquito factor Xa inhibitor, Pipiserpin, which affects female reproductive potential. Our results suggest that Pipiserpin may be a novel target for mosquito population control. The conclusions from our study on Cx. pipiens pallens might serve as a reference for the development of control measures for other mosquitoes as well. © 2022 Society of Chemical Industry.

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.

Chromosome-level genome assembly, annotation, and phylogenomics of the gooseneck barnacle Pollicipes pollicipes

BACKGROUND

The barnacles are a group of >2,000 species that have fascinated biologists, including Darwin, for centuries. Their lifestyles are extremely diverse, from free-swimming larvae to sessile adults, and even root-like endoparasites. Barnacles also cause hundreds of millions of dollars of losses annually due to biofouling. However, genomic resources for crustaceans, and barnacles in particular, are lacking.

RESULTS

Using 62× Pacific Biosciences coverage, 189× Illumina whole-genome sequencing coverage, 203× HiC coverage, and 69× CHi-C coverage, we produced a chromosome-level genome assembly of the gooseneck barnacle Pollicipes pollicipes. The P. pollicipes genome is 770 Mb long and its assembly is one of the most contiguous and complete crustacean genomes available, with a scaffold N50 of 47 Mb and 90.5% of the BUSCO Arthropoda gene set. Using the genome annotation produced here along with transcriptomes of 13 other barnacle species, we completed phylogenomic analyses on a nearly 2 million amino acid alignment. Contrary to previous studies, our phylogenies suggest that the Pollicipedomorpha is monophyletic and sister to the Balanomorpha, which alters our understanding of barnacle larval evolution and suggests homoplasy in a number of naupliar characters. We also compared transcriptomes of P. pollicipes nauplius larvae and adults and found that nearly one-half of the genes in the genome are differentially expressed, highlighting the vastly different transcriptomes of larvae and adult gooseneck barnacles. Annotation of the genes with KEGG and GO terms reveals that these stages exhibit many differences including cuticle binding, chitin binding, microtubule motor activity, and membrane adhesion.

CONCLUSION

This study provides high-quality genomic resources for a key group of crustaceans. This is especially valuable given the roles P. pollicipes plays in European fisheries, as a sentinel species for coastal ecosystems, and as a model for studying barnacle adhesion as well as its key position in the barnacle tree of life. A combination of genomic, phylogenetic, and transcriptomic analyses here provides valuable insights into the evolution and development of barnacles.

Genome and transcriptome analysis of the beet armyworm Spodoptera exigua reveals targets for pest control

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies.

Spatial and morphological reorganization of endosymbiosis during metamorphosis accommodates adult metabolic requirements in a weevil

Bacterial intracellular symbiosis (endosymbiosis) is widespread in nature and impacts many biological processes. In holometabolous symbiotic insects, metamorphosis entails a complete and abrupt internal reorganization that creates a constraint for endosymbiont transmission from larvae to adults. To assess how endosymbiosis copes-and potentially evolves-throughout this major host-tissue reorganization, we used the association between the cereal weevil Sitophilus oryzae and the bacterium Sodalis pierantonius as a model system. S. pierantonius are contained inside specialized host cells, the bacteriocytes, that group into an organ, the bacteriome. Cereal weevils require metabolic inputs from their endosymbiont, particularly during adult cuticle synthesis, when endosymbiont load increases dramatically. By combining dual RNA-sequencing analyses and cell imaging, we show that the larval bacteriome dissociates at the onset of metamorphosis and releases bacteriocytes that undergo endosymbiosis-dependent transcriptomic changes affecting cell motility, cell adhesion, and cytoskeleton organization. Remarkably, bacteriocytes turn into spindle cells and migrate along the midgut epithelium, thereby conveying endosymbionts to midgut sites where future mesenteric caeca will develop. Concomitantly, endosymbiont genes encoding a type III secretion system and a flagellum apparatus are transiently up-regulated while endosymbionts infect putative stem cells and enter their nuclei. Infected cells then turn into new differentiated bacteriocytes and form multiple new bacteriomes in adults. These findings show that endosymbiosis reorganization in a holometabolous insect relies on a synchronized host-symbiont molecular and cellular "choreography" and illustrates an adaptive feature that promotes bacteriome multiplication to match increased metabolic requirements in emerging adults.

Genetic engineering and bacterial pathogenesis against the vectorial capacity of mosquitoes

Mosquitoes are the main vector of multiple diseases worldwide and transmit viral (malaria, chikungunya, encephalitis, yellow fever, as well as dengue fever), as well as bacterial diseases (tularemia). To manage the outbreak of mosquito populations, various management programs include the application of chemicals, followed by biological and genetic control. Here we aimed to focus on the role of bacterial pathogenesis and molecular tactics for the management of mosquitoes and their vectorial capacity. Bacterial pathogenesis and molecular manipulations have a substantial impact on the biology of mosquitoes, and both strategies change the gene expression and regulation of disease vectors. The strategy for genetic modification is also proved to be excellent for the management of mosquitoes, which halt the development of population via incompatibility of different sex. Therefore, the purpose of the present discussion is to illustrate the impact of both approaches against the vectorial capacity of mosquitoes. Moreover, it could be helpful to understand the relationship of insect-pathogen and to manage various insect vectors as well as diseases.

Development of the transcriptome for a sediment ecotoxicological model species, Chironomus dilutus

Chironomus dilutus is a prominent model species in conventional sediment toxicity testing and sediment contamination diagnosis. However, lack of genomic data significantly limited its application in identifying toxicological mode of action (MOA) and molecular biomarkers of toxicants. Here the transcriptome of C. dilutus in full life span and both sexes (1st, 2nd, 3rd and 4th instar larvae, pupae, and adults) were developed and temporal gene expression across adjacent life stages were investigated to understand the regulation of development. Furthermore, transcriptional response of Midges (the 4th instar larvae) exposed to chemicals of different MOAs (CdCl₂, nonylphenol and triclosan) were profiled based on the reference transcriptome. Consequently, a complete transcriptome of 31132 unigenes with N50 of 3117bp, covering 98.8% of the arthropod single-copy orthologs were assembled. While 364 genes were differentially expressed among adjacent larval stages, 7142 and 2127 of transcripts were significantly changed for the transition of larvae-pupae and pupae-adults, respectively. Finally, chemical-specific gene expression profile were identified in the midges, showed its potential in classifying distinct contaminants. Overall, the comprehensive transcriptome of C. dilutus developed here could not only facilitate the mechanistic understanding of environmental toxicants during critical life stage of aquatic insects, but also provide molecular diagnostic tools in sediment ecotoxicology.

The mosquito melanization response requires hierarchical activation of non-catalytic clip domain serine protease homologs

Serine protease cascades regulate important insect immune responses namely melanization and Toll pathway activation. An important component of these cascades are clip-domain serine protease homologs (cSPHs), which are non-catalytic, but essential for activating the enzyme prophenoloxidase (PPO) in the melanization response during septic infections. The activation of cSPHs requires their proteolytic cleavage, yet factors that control their activation and the complexity of their interactions within these cascades remain unclear. Here, we report the identification of CLIPA28 as a novel immune-related cSPH in the malaria vector Anopheles gambiae. Functional genetic analysis using RNA interference (RNAi) revealed that CLIPA28 is essential for the melanization of Plasmodium berghei parasites in refractory mosquitoes, and for mosquito resistance to fungal infections. We further show, using combined biochemical and genetic approaches, that CLIPA28 is member of a network of at least four cSPHs, whereby members are activated in a hierarchical manner following septic infections. Depletion of the complement-like protein TEP1 abolished the activation of this network after septic infections, whereas, depletion of the serine protease inhibitor 2 (SRPN2) triggered enhanced network activation, even in naïve mosquitoes, culminating in a dramatic reduction in cSPHs hemolymph levels, which paralleled that of PPO. Our data suggest that cSPHs are engaged in complex and multilayered interactions within serine protease cascades that regulate melanization, and identify TEP1 and SRPN2 as two master regulators of the cSPH network.

Melanization is a spectacular immune response of insects that culminates in the deposition of melanin on microbial surfaces leading to their death. Despite more than half a century of investigation of this immune reaction in several insect models, important knowledge gaps, specifically concerning the regulation of this response, remain, especially that its output is tightly regulated by complex genetic and biochemical interactions between members of structurally diverse immune gene families. Here, we identify a new hemolymph protein that is involved in the melanization of malaria parasites and show that it forms, together with other proteins of the same gene family, a network, whereby members are activated in an orderly manner during mosquito septic infections. Furthemore, we show that this network is controlled by two structurally distinct immune proteins which act as master regulators of its activation. Dissecting the complexity of the regulatory mechanisms mediating insect melanization responses, may facilitate fine manipulation of mosquito immunity to reduce vectorial capacity.

Genome-wide analysis of developmental stage-specific transcriptome in Bradysia odoriphaga

Bradysia odoriphaga is a serious pest of the Chinese chive; however, detailed information regarding the developmental stage-specific gene expression patterns of B. odoriphaga is not yet available. In this study, RNA sequencing (RNA-seq) was performed to determine the gene expression patterns of developmental stages including the eggs, second instar larvae, fourth instar larvae, pupae, and adults of B. odoriphaga. Analysis of 15 samples revealed an average of 89.56% of the clean reads could be mapped onto the assembled UniGene database. Cluster tree analysis showed that the expression patterns were stage-specific and that samples of the second and fourth instar larvae clustered in one group, while those of eggs, pupae, and adults clustered in another group. Differential expression unigenes (DEUs) for sequential developmental stages were between 3314 and 10,632. A total of 1910-7756 DEUs of sequential developmental stages were assigned into 45-56 gene ontology categories and 1165-3845 DEUs were mapped into Kyoto Encyclopedia of Genes and Genomes pathways. The expression of DEUs related to growth and development showed that hormone receptors highly expressed in the pupal stage, while chitinases were highly expressed in the larval stage. The results of quantitative real time polymerase chain reaction (qRT-PCR) and RNA-seq expression agreed well for 12 growth- and development-related unigenes. This study identified DEUs for sequential developmental stages of B. odoriphaga. Gene Ontology classifications and KEGG pathway identification of DEUs not only provide information useful for understanding insect growth and development but also for exploring novel approaches to control B. odoriphaga.

Growth-disrupting Murraya koenigii leaf extracts on Anopheles gambiae larvae and identification of associated candidate bioactive constituents

Plant-based constituents have been proposed as eco-friendly alternatives to synthetic insecticides for control of mosquito vectors of malaria. In this study, we first screened the effects of methanolic leaf extracts of curry tree (Murraya koenigii) growing in tropical (Mombasa, Malindi) and semi-arid (Kibwezi, and Makindu) ecological zones of Kenya on third instar An. gambiae s.s. larvae. Extracts of the plant from the semi-arid region, and particularly from Kibwezi, led to high mortality of the larvae. Bioassay-guided fractionation of the methanolic extract of the leaves of the plants from Kibwezi was then undertaken and the most active fraction (20 fold more potent than the crude extract) was then analyzed by Liquid chromatography quadruple time of flight coupled with mass spectrometry (LC-QtoF-MS) and a number of constituents were identified, including a major alkaloid constituent, Neplanocin A (5). Exposure of the third instar larvae to a sub-lethal dose (4.43 ppm) of this fraction over 7-day periods induced gross morphogenetic abnormalities in the larvae, with reduced locomotion, and delayed pupation. Moreover, the few adults that emerged from some pupae failed to fly from the water surface, unlike in the untreated control group. These results demonstrate subtle growth-disrupting effects of the phytochemical blend from M. koenigii leaves on aquatic stages An. gambiae mosquito. The study lays down some useful groundwork for the downstream development of phytochemical blends that can be evaluated for integration into eco-friendly control of An. gambiae vector population targeting the often overlooked but important immature stages of the malaria vector.

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.

Differential gene transcription across the life cycle in Daphnia magna using a new all genome custom-made microarray

Background

Unravelling the link between genes and environment across the life cycle is a challenging goal that requires model organisms with well-characterized life-cycles, ecological interactions in nature, tractability in the laboratory, and available genomic tools. Very few well-studied invertebrate model species meet these requirements, being the waterflea Daphnia magna one of them. Here we report a full genome transcription profiling of D. magna during its life-cycle. The study was performed using a new microarray platform designed from the complete set of gene models representing the whole transcribed genome of D. magna.

Results

Up to 93% of the existing 41,317 D. magna gene models showed differential transcription patterns across the developmental stages of D. magna, 59% of which were functionally annotated. Embryos showed the highest number of unique transcribed genes, mainly related to DNA, RNA, and ribosome biogenesis, likely related to cellular proliferation and morphogenesis of the several body organs. Adult females showed an enrichment of transcripts for genes involved in reproductive processes. These female-specific transcripts were essentially absent in males, whose transcriptome was enriched in specific genes of male sexual differentiation genes, like doublesex.

Conclusion

Our results define major characteristics of transcriptional programs involved in the life-cycle, differentiate males and females, and show that large scale gene-transcription data collected in whole animals can be used to identify genes involved in specific biological and biochemical processes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4725-7) contains supplementary material, which is available to authorized users.

Comparative transcriptomics reveal developmental turning points during embryogenesis of a hemimetabolous insect, the damselfly Ischnura elegans

Identifying transcriptional changes during embryogenesis is of crucial importance for unravelling evolutionary, molecular and cellular mechanisms that underpin patterning and morphogenesis. However, comparative studies focusing on early/embryonic stages during insect development are limited to a few taxa. Drosophila melanogaster is the paradigm for insect development, whereas comparative transcriptomic studies of embryonic stages of hemimetabolous insects are completely lacking. We reconstructed the first comparative transcriptome covering the daily embryonic developmental progression of the blue-tailed damselfly Ischnura elegans (Odonata), an ancient hemimetabolous representative. We identified a "core" set of 6,794 transcripts - shared by all embryonic stages - which are mainly involved in anatomical structure development and cellular nitrogen compound metabolic processes. We further used weighted gene co-expression network analysis to identify transcriptional changes during Odonata embryogenesis. Based on these analyses distinct clusters of transcriptional active sequences could be revealed, indicating that embryos at different development stages have their own transcriptomic profile according to the developmental events and leading to sequential reprogramming of metabolic and developmental genes. Interestingly, a major change in transcriptionally active sequences is correlated with katatrepsis (revolution) during mid-embryogenesis, a 180° rotation of the embryo within the egg and specific to hemimetabolous insects.

Anopheles gambiae larvae mount stronger immune responses against bacterial infection than adults: evidence of adaptive decoupling in mosquitoes

BACKGROUND

The immune system of adult mosquitoes has received significant attention because of the ability of females to vector disease-causing pathogens while ingesting blood meals. However, few studies have focused on the immune system of larvae, which, we hypothesize, is highly robust due to the high density and diversity of microorganisms that larvae encounter in their aquatic environments and the strong selection pressures at work in the larval stage to ensure survival to reproductive maturity. Here, we surveyed a broad range of cellular and humoral immune parameters in larvae of the malaria mosquito, Anopheles gambiae, and compared their potency to that of newly-emerged adults and older adults.

RESULTS

We found that larvae kill bacteria in their hemocoel with equal or greater efficiency compared to newly-emerged adults, and that antibacterial ability declines further with adult age, indicative of senescence. This phenotype correlates with more circulating hemocytes and a differing spatial arrangement of sessile hemocytes in larvae relative to adults, as well as with the individual hemocytes of adults carrying a greater phagocytic burden. The hemolymph of larvae also possesses markedly stronger antibacterial lytic and melanization activity than the hemolymph of adults. Finally, infection induces a stronger transcriptional upregulation of immunity genes in larvae than in adults, including differences in the immunity genes that are regulated.

CONCLUSIONS

These results demonstrate that immunity is strongest in larvae and declines after metamorphosis and with adult age, and suggest that adaptive decoupling, or the independent evolution of larval and adult traits made possible by metamorphosis, has occurred in the mosquito lineage.

The Anopheles gambiae transcriptome - a turning point for malaria control

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

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.

Transcriptomic differences between euryhaline and stenohaline malaria vector sibling species in response to salinity stress

Evolution of osmoregulatory systems is a key factor in the transition of species between fresh- and saltwater habitats. Anopheles coluzzii and Anopheles merus are stenohaline and euryhaline malaria vector mosquitoes belonging to a larger group of sibling species, the Anopheles gambiae complex, which radiated in Africa within the last 2 million years. Comparative ecological genomics of these vector species can provide insight into the mechanisms that permitted the rapid radiation of this species complex into habitats of contrasting salinity. Here, we use RNA-Seq to investigate gene expression differences between An. coluzzii and An. merus after briefly exposing both young and old larval instars of each species to either saltwater (SW) or freshwater (FW). Our study aims to identify candidate genes and pathways responsible for the greater SW tolerance of An. merus. Our results are congruent with the ability of gene induction to mediate salinity tolerance, with both species showing increasing amounts of differential gene expression between SW and FW as salt concentrations increase. Besides ion transporters such as AgAE2 that may serve as effectors for osmoregulation, we also find mitogen-activated protein kinases that may serve in a phosphorylation signalling pathway responding to salinity, and report potential cross-talk between the mosquito immune response and osmoregulation. This study provides a key step towards applying the growing molecular knowledge of these malaria vectors to improve understanding of their ecological tolerances and habitat occupancy.

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

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

Long non-coding RNA discovery across the genus anopheles reveals conserved secondary structures within and beyond the Gambiae complex

Background

Long non-coding RNAs (lncRNAs) have been defined as mRNA-like transcripts longer than 200 nucleotides that lack significant protein-coding potential, and many of them constitute scaffolds for ribonucleoprotein complexes with critical roles in epigenetic regulation. Various lncRNAs have been implicated in the modulation of chromatin structure, transcriptional and post-transcriptional gene regulation, and regulation of genomic stability in mammals, Caenorhabditis elegans, and Drosophila melanogaster. The purpose of this study is to identify the lncRNA landscape in the malaria vector An. gambiae and assess the evolutionary conservation of lncRNAs and their secondary structures across the Anopheles genus.

Results

Using deep RNA sequencing of multiple Anopheles gambiae life stages, we have identified 2,949 lncRNAs and more than 300 previously unannotated putative protein-coding genes. The lncRNAs exhibit differential expression profiles across life stages and adult genders. We find that across the genus Anopheles, lncRNAs display much lower sequence conservation than protein-coding genes. Additionally, we find that lncRNA secondary structure is highly conserved within the Gambiae complex, but diverges rapidly across the rest of the genus Anopheles.

Conclusions

This study offers one of the first lncRNA secondary structure analyses in vector insects. Our description of lncRNAs in An. gambiae offers the most comprehensive genome-wide insights to date into lncRNAs in this vector mosquito, and defines a set of potential targets for the development of vector-based interventions that may further curb the human malaria burden in disease-endemic countries.

Electronic supplementary material

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

Deciphering life history transcriptomes in different environments

We compared whole transcriptome variation in six pre-adult stages and seven adult female ages in two populations of cactophilic Drosophila mojavensis reared on two host plants to understand how differences in gene expression influence standing life history variation. We used singular value decomposition (SVD) to identify dominant trajectories of life cycle gene expression variation, performed pairwise comparisons of stage and age differences in gene expression across the life cycle, identified when genes exhibited maximum levels of life cycle gene expression, and assessed population and host cactus effects on gene expression. Life cycle SVD analysis returned four significant components of transcriptional variation, revealing functional enrichment of genes responsible for growth, metabolic function, sensory perception, neural function, translation and ageing. Host cactus effects on female gene expression revealed population- and stage-specific differences, including significant host plant effects on larval metabolism and development, as well as adult neurotransmitter binding and courtship behaviour gene expression levels. In 3- to 6-day-old virgin females, significant upregulation of genes associated with meiosis and oogenesis was accompanied by downregulation of genes associated with somatic maintenance, evidence for a life history trade-off. The transcriptome of D. mojavensis reared in natural environments throughout its life cycle revealed core developmental transitions and genome-wide influences on life history variation in natural populations.

Dissecting the organ specificity of insecticide resistance candidate genes in Anopheles gambiae: known and novel candidate genes

Background

The elevated expression of enzymes with insecticide metabolism activity can lead to high levels of insecticide resistance in the malaria vector, Anopheles gambiae. In this study, adult female mosquitoes from an insecticide susceptible and resistant strain were dissected into four different body parts. RNA from each of these samples was used in microarray analysis to determine the enrichment patterns of the key detoxification gene families within the mosquito and to identify additional candidate insecticide resistance genes that may have been overlooked in previous experiments on whole organisms.

Results

A general enrichment in the transcription of genes from the four major detoxification gene families (carboxylesterases, glutathione transferases, UDP glucornyltransferases and cytochrome P450s) was observed in the midgut and malpighian tubules. Yet the subset of P450 genes that have previously been implicated in insecticide resistance in An gambiae, show a surprisingly varied profile of tissue enrichment, confirmed by qPCR and, for three candidates, by immunostaining. A stringent selection process was used to define a list of 105 genes that are significantly (p ≤0.001) over expressed in body parts from the resistant versus susceptible strain. Over half of these, including all the cytochrome P450s on this list, were identified in previous whole organism comparisons between the strains, but several new candidates were detected, notably from comparisons of the transcriptomes from dissected abdomen integuments.

Conclusions

The use of RNA extracted from the whole organism to identify candidate insecticide resistance genes has a risk of missing candidates if key genes responsible for the phenotype have restricted expression within the body and/or are over expression only in certain tissues. However, as transcription of genes implicated in metabolic resistance to insecticides is not enriched in any one single organ, comparison of the transcriptome of individual dissected body parts cannot be recommended as a preferred means to identify new candidate insecticide resistant genes. Instead the rich data set on in vivo sites of transcription should be consulted when designing follow up qPCR validation steps, or for screening known candidates in field populations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1018) contains supplementary material, which is available to authorized users.

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.

Insights into the epigenomic landscape of the human malaria vector Anopheles gambiae

The epigenome of the human malaria vector Anopheles gambiae was characterized in midgut cells by mapping the distribution and levels of two post-translational histone modifications, H3K27ac and H3K27me3. These histone profiles were then correlated with levels of gene expression obtained by RNA-seq. Analysis of the transcriptome of A. gambiae midguts and salivary glands led to the discovery of 13,898 new transcripts not present in the most recent genome assembly. A subset of these transcripts is differentially expressed between midgut and salivary glands. The enrichment profiles of H3K27ac and H3K27me3 are mutually exclusive and associate with high and low levels of transcription, respectively. This distribution agrees with previous findings in Drosophila showing association of these two histone modifications with either active or inactive transcriptional states, including Polycomb-associated domains in silenced genes. This study provides a mosquito epigenomics platform for future comparative studies in other mosquito species, opening future investigations into the role of epigenetic processes in vector-borne systems of medical and economic importance.

Gene expression divergence between malaria vector sibling species Anopheles gambiae and An. coluzzii from rural and urban Yaoundé Cameroon

Divergent selection based on aquatic larval ecology is a likely factor in the recent isolation of two broadly sympatric and morphologically identical African mosquito species, the malaria vectors Anopheles gambiae and An. coluzzii. Population-based genome scans have revealed numerous candidate regions of recent positive selection, but have provided few clues as to the genetic mechanisms underlying behavioural and physiological divergence between the two species, phenotypes which themselves remain obscure. To uncover possible genetic mechanisms, we compared global transcriptional profiles of natural and experimental populations using gene-based microarrays. Larvae were sampled as second and fourth instars from natural populations in and around the city of Yaoundé, capital of Cameroon, where the two species segregate along a gradient of urbanization. Functional enrichment analysis of differentially expressed genes revealed that An. coluzzii--the species that breeds in more stable, biotically complex and potentially polluted urban water bodies--overexpresses genes implicated in detoxification and immunity relative to An. gambiae, which breeds in more ephemeral and relatively depauperate pools and puddles in suburbs and rural areas. Moreover, our data suggest that such overexpression by An. coluzzii is not a transient result of induction by xenobiotics in the larval habitat, but an inherent and presumably adaptive response to repeatedly encountered environmental stressors. Finally, we find no significant overlap between the differentially expressed loci and previously identified genomic regions of recent positive selection, suggesting that transcriptome divergence is regulated by trans-acting factors rather than cis-acting elements.

Transcriptome analysis of Anopheles stephensi embryo using expressed sequence tags

Germ band retraction (GBR) stage is one of the important stages during insect development. It is associated with an extensive epithelial morphogenesis and may also be pivotal in generation of morphological diversity in insects. Despite its importance, only a handful of studies report the transcriptome repertoire of this stage in insects. Here, we report generation, annotation and analysis of ESTs from the embryonic stage (16-22 h post fertilization) of laboratoryreared Anopheles stephensi mosquitoes. A total of 1002 contigs were obtained upon clustering of 1140 high-quality ESTs, which demonstrates an astonishingly low transcript redundancy (12.1 percent). Putative functions were assigned only to 213 contigs (21 percent), comprising mainly of transcripts encoding protein synthesis machinery. Approximately 78 percent of the transcripts remain uncharacterized, illustrating a lack of sequence information about the genes expressed in the embryonic stages of mosquitoes. This study highlights several novel transcripts, which apart from insect development, may significantly contribute to the essential biological complexity underlying insect viability in adverse environments. Nonetheless, the generated sequence information from this work provides a comprehensive resource for genome annotation, microarray development, phylogenetic analysis and other molecular biology applications in entomology.

The CLIP-domain serine protease homolog SPCLIP1 regulates complement recruitment to microbial surfaces in the malaria mosquito Anopheles gambiae

The complement C3-like protein TEP1 of the mosquito Anopheles gambiae is required for defense against malaria parasites and bacteria. Two forms of TEP1 are present in the mosquito hemolymph, the full-length TEP1-F and the proteolytically processed TEP1_cut that is part of a complex including the leucine-rich repeat proteins LRIM1 and APL1C. Here we show that the non-catalytic serine protease SPCLIP1 is a key regulator of the complement-like pathway. SPCLIP1 is required for accumulation of TEP1 on microbial surfaces, a reaction that leads to lysis of malaria parasites or triggers activation of a cascade culminating with melanization of malaria parasites and bacteria. We also demonstrate that the two forms of TEP1 have distinct roles in the complement-like pathway and provide the first evidence for a complement convertase-like cascade in insects analogous to that in vertebrates. Our findings establish that core principles of complement activation are conserved throughout the evolution of animals.

Mosquitoes are vectors of numerous human diseases including malaria. Disease transmission requires that microbes overcome the robust mosquito immune system. In the African malaria mosquito, the TEP1 protein that is homologous to mammalian complement factor C3 is shown to play a central role in mosquito immunity to malaria parasites and bacteria. In this study, we report that another mosquito protein belonging to a class of non-catalytic enzymes that are specific to arthropods is a core component of the mosquito complement-like immune pathway. We found that this new protein, named SPCLIP1, regulates the accumulation of TEP1 on malaria parasites and bacteria, and show that this can lead to distinct defense reactions including lysis and melanization of the pathogen. This work is valuable because it reveals novel insight into the regulation of mosquito complement on microbial surfaces such as those of the malaria parasites. Unraveling the molecular mechanisms regulating these defense responses may ultimately lead to the design of novel disease blocking strategies in the vector.

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.

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.

Biochemical properties, expression profiles, and tissue localization of orthologous acetylcholinesterase-2 in the mosquito, Anopheles gambiae

Acetylcholinesterases (AChEs) catalyze the hydrolysis of acetylcholine, a neurotransmitter for cholinergic neurotransmission in animals. Most insects studied so far possess two AChE genes: ace-1 paralogous and ace-2 orthologous to Drosophila melanogaster ace. We characterized the catalytic domain of Anopheles gambiae AChE1 in a previous study (Jiang et al., 2009) and report here biochemical properties of A. gambiae AChE2 expressed in Sf9 cells. An unknown protease in the expression system cleaved the recombinant AChE2 next to Arg(110), yielding two non-covalently associated polypeptides. A mixture of the intact and cleaved AChE2 had a specific activity of 72.3 U/mg, much lower than that of A. gambiae AChE1 (523 U/mg). The order of V(max)/K(M) values for the model substrates was acetylthiocholine > propionylthiocholine ≈ acetyl-(β-methyl)thiocholine > butyrylthiocholine. The IC(50)'s for eserine, carbaryl, BW284C51, paraoxon and malaoxon were 1.32, 13.6, 26.8, 192 and 294 nM, respectively. A. gambiae AChE2 bound eserine and carbaryl stronger than paraoxon and malaoxon, whereas eserine and malaoxon modified the active site Ser(232) faster than carbaryl or paraoxon did. Consequently, the k(i)'s were 1.173, 0.245, 0.029 and 0.018 μM(-1)min(-1) for eserine, carbaryl, paraoxon and malaoxon, respectively. Quantitative polymerase chain reactions showed a similar pattern of ace-1 and ace-2 expression. Their mRNAs were abundant in early embryos, greatly decreased in late embryos, larvae, pupae, and pharate adult, and became abundant again in adults. Both transcripts were higher in head and abdomen than thorax of adults and higher in male than female mosquitoes. Transcript levels of ace-1 were 1.9- to 361.8-fold higher than those of ace-2, depending on developmental stages and body parts. Cross-reacting polyclonal antibodies detected AChEs in adult brains, thoracic ganglia, and genital/rectal area. Activity assays, immunoblotting, and tandem mass spectrometric analysis indicated that A. gambiae AChE1 is responsible for most of acetylthiocholine hydrolysis in the head extracts. Taken together, these data indicate that A. gambiae AChE2 may play a less significant role than AChE1 does in the mosquito nervous system.

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.

Post-embryonic transcriptomes of the prawn Macrobrachium rosenbergii: multigenic succession through metamorphosis

Like many metazoans, the freshwater prawn Macrobrachium rosenbergii begins its post-embryonic life with a set of morphologically distinct planktonic larval stages, followed by a benthic post-larval stage during which the maturing organism differs from the larvae both ecologically and physiologically. Understanding of the molecular basis underlying morphogenesis in crustaceans is limited to the observation that methyl farnesoate, the non-epoxidated form of the insect juvenile hormone, acts as the active crustacean juvenoid. Molt steroids were also linked to morphogenesis and several other molecular pathways, such as Hedgehog and Wnt, are known to underlie morphogenesis in all metazoans examined and, as such, are thought to do the same in crustaceans. Using next generation sequencing, we deep-sequenced the transcriptomes of several larval and post-larval stages. De novo assembly, followed by bioinformatics analysis, revealed that many novel transcripts are over-expressed in either larvae- or post-larvae-stage prawn, shedding light on the molecular basis underlying M. rosenbergii metamorphosis. Fast larval molting rates and periodic morphological changes were reflected in over-expression of transcripts annotated to the cell cycle, DNA replication and morphogenic pathways (i.e., Hedgehog and Wnt). Further characterization of transcripts assigned to morphogenic pathways by real-time RT-PCR reconfirmed their over-expression in larvae, albeit with a more complex expression pattern when examined in the individual developmental stages. The expression level of an orthologue of cytochrome P450, 15A1, known to epoxidize methyl farnesoate in insects, was increased in the late larval and early post-larval stages, in accordance with the role of methyl farnesoate in crustacean metamorphosis. This study exemplifies the applicability of a high-throughput sequencing approach for studying complex traits, including metamorphosis, providing new insight into this unexplored area of crustacean research.

Nature, nurture and evolution of intra-species variation in mosquito arbovirus transmission competence

Mosquitoes vary in their competence or ability to transmit arthropod-borne viruses (arboviruses). Many arboviruses cause disease in humans and animals. Identifying the environmental and genetic causes of variation in mosquito competence for arboviruses is one of the great challenges in public health. Progress identifying genetic (nature) and environmental (nurture) factors influencing mosquito competence for arboviruses is reviewed. There is great complexity in the various traits that comprise mosquito competence. The complex interactions between environmental and genetic factors controlling these traits and the factors shaping variation in Nature are largely unknown. The norms of reaction of specific genes influencing competence, their distributions in natural populations and the effects of genetic polymorphism on phenotypic variation need to be determined. Mechanisms influencing competence are not likely due to natural selection because of the direct effects of the arbovirus on mosquito fitness. More likely the traits for mosquito competence for arboviruses are the effects of adaptations for other functions of these competence mechanisms. Determining these other functions is essential to understand the evolution and distributions of competence for arboviruses. This information is needed to assess risk from mosquito-borne disease, predict new mosquito-arbovirus systems, and provide novel strategies to mitigate mosquito-borne arbovirus transmission.

Transcriptome of the adult female malaria mosquito vector Anopheles albimanus

BACKGROUND

Human Malaria is transmitted by mosquitoes of the genus Anopheles. Transmission is a complex phenomenon involving biological and environmental factors of humans, parasites and mosquitoes. Among more than 500 anopheline species, only a few species from different branches of the mosquito evolutionary tree transmit malaria, suggesting that their vectorial capacity has evolved independently. Anopheles albimanus (subgenus Nyssorhynchus) is an important malaria vector in the Americas. The divergence time between Anopheles gambiae, the main malaria vector in Africa, and the Neotropical vectors has been estimated to be 100 My. To better understand the biological basis of malaria transmission and to develop novel and effective means of vector control, there is a need to explore the mosquito biology beyond the An. gambiae complex.

RESULTS

We sequenced the transcriptome of the An. albimanus adult female. By combining Sanger, 454 and Illumina sequences from cDNA libraries derived from the midgut, cuticular fat body, dorsal vessel, salivary gland and whole body, we generated a single, high-quality assembly containing 16,669 transcripts, 92% of which mapped to the An. darlingi genome and covered 90% of the core eukaryotic genome. Bidirectional comparisons between the An. gambiae, An. darlingi and An. albimanus predicted proteomes allowed the identification of 3,772 putative orthologs. More than half of the transcripts had a match to proteins in other insect vectors and had an InterPro annotation. We identified several protein families that may be relevant to the study of Plasmodium-mosquito interaction. An open source transcript annotation browser called GDAV (Genome-Delinked Annotation Viewer) was developed to facilitate public access to the data generated by this and future transcriptome projects.

CONCLUSIONS

We have explored the adult female transcriptome of one important New World malaria vector, An. albimanus. We identified protein-coding transcripts involved in biological processes that may be relevant to the Plasmodium lifecycle and can serve as the starting point for searching targets for novel control strategies. Our data increase the available genomic information regarding An. albimanus several hundred-fold, and will facilitate molecular research in medical entomology, evolutionary biology, genomics and proteomics of anopheline mosquito vectors. The data reported in this manuscript is accessible to the community via the VectorBase website (http://www.vectorbase.org/Other/AdditionalOrganisms/).

Identification of early zygotic genes in the yellow fever mosquito Aedes aegypti and discovery of a motif involved in early zygotic genome activation

During early embryogenesis the zygotic genome is transcriptionally silent and all mRNAs present are of maternal origin. The maternal-zygotic transition marks the time over which embryogenesis changes its dependence from maternal RNAs to zygotically transcribed RNAs. Here we present the first systematic investigation of early zygotic genes (EZGs) in a mosquito species and focus on genes involved in the onset of transcription during 2–4 hr. We used transcriptome sequencing to identify the “pure” (without maternal expression) EZGs by analyzing transcripts from four embryonic time ranges of 0–2, 2–4, 4–8, and 8–12 hr, which includes the time of cellular blastoderm formation and up to the start of gastrulation. Blast of 16,789 annotated transcripts vs. the transcriptome reads revealed evidence for 63 (P<0.001) and 143 (P<0.05) nonmaternally derived transcripts having a significant increase in expression at 2–4 hr. One third of the 63 EZG transcripts do not have predicted introns compared to 10% of all Ae. aegypti genes. We have confirmed by RT-PCR that zygotic transcription starts as early as 2–3 hours. A degenerate motif VBRGGTA was found to be overrepresented in the upstream sequences of the identified EZGs using a motif identification software called SCOPE. We find evidence for homology between this motif and the TAGteam motif found in Drosophila that has been implicated in EZG activation. A 38 bp sequence in the proximal upstream sequence of a kinesin light chain EZG (KLC2.1) contains two copies of the mosquito motif. This sequence was shown to support EZG transcription by luciferase reporter assays performed on injected early embryos, and confers early zygotic activity to a heterologous promoter from a divergent mosquito species. The results of these studies are consistent with the model of early zygotic genome activation via transcriptional activators, similar to what has been found recently in Drosophila.

Transcription profiling associated with life cycle of Anopheles gambiae

Complex biological events occur during the developmental process of the mosquito Anopheles gambiae (Giles). Using cDNA expression microarrays, the expression patterns of 13,440 clones representing 8,664 unique transcripts were revealed from six different developmental stages: early larvae (late third instar/early fourth instar), late larvae (late fourth instar), early pupae (< 30 min after pupation), late pupae (after tanning), and adult female and male mosquitoes (24 h postemergence). After microarray analysis, 560 unique transcripts were identified to show at least a fourfold up- or down-regulation in at least one developmental stage. Based on the expression patterns, these gene products were clustered into 13 groups. In total, eight genes were analyzed by quantitative real-time polymerase chain reaction to validate microarray results. Among 560 unique transcripts, 446 contigs were assigned to respective genes from the An. gambiae genome. The expression patterns and annotations of the genes in the 13 groups are discussed in the context of development including metabolism, transport, protein synthesis and degradation, cellular processes, cellular communication, intra- or extra-cellular architecture maintenance, response to stress or immune-related defense, and spermatogenesis.

An expression map for Anopheles gambiae

BACKGROUND

Quantitative transcriptome data for the malaria-transmitting mosquito Anopheles gambiae covers a broad range of biological and experimental conditions, including development, blood feeding and infection. Web-based summaries of differential expression for individual genes with respect to these conditions are a useful tool for the biologist, but they lack the context that a visualisation of all genes with respect to all conditions would give. For most organisms, including A. gambiae, such a systems-level view of gene expression is not yet available.

RESULTS

We have clustered microarray-based gene-averaged expression values, available from VectorBase, for 10194 genes over 93 experimental conditions using a self-organizing map. Map regions corresponding to known biological events, such as egg production, are revealed. Many individual gene clusters (nodes) on the map are highly enriched in biological and molecular functions, such as protein synthesis, protein degradation and DNA replication. Gene families, such as odorant binding proteins, can be classified into distinct functional groups based on their expression and evolutionary history. Immunity-related genes are non-randomly distributed in several distinct regions on the map, and are generally distant from genes with house-keeping roles. Each immunity-rich region appears to represent a distinct biological context for pathogen recognition and clearance (e.g. the humoral and gut epithelial responses). Several immunity gene families, such as peptidoglycan recognition proteins (PGRPs) and defensins, appear to be specialised for these distinct roles, while three genes with physically interacting protein products (LRIM1/APL1C/TEP1) are found in close proximity.

CONCLUSIONS

The map provides the first genome-scale, multi-experiment overview of gene expression in A. gambiae and should also be useful at the gene-level for investigating potential interactions. A web interface is available through the VectorBase website http://www.vectorbase.org/. It is regularly updated as new experimental data becomes available.

Comparative genomic analysis of Drosophila melanogaster and vector mosquito developmental genes

Genome sequencing projects have presented the opportunity for analysis of developmental genes in three vector mosquito species: Aedes aegypti, Culex quinquefasciatus, and Anopheles gambiae. A comparative genomic analysis of developmental genes in Drosophila melanogaster and these three important vectors of human disease was performed in this investigation. While the study was comprehensive, special emphasis centered on genes that 1) are components of developmental signaling pathways, 2) regulate fundamental developmental processes, 3) are critical for the development of tissues of vector importance, 4) function in developmental processes known to have diverged within insects, and 5) encode microRNAs (miRNAs) that regulate developmental transcripts in Drosophila. While most fruit fly developmental genes are conserved in the three vector mosquito species, several genes known to be critical for Drosophila development were not identified in one or more mosquito genomes. In other cases, mosquito lineage-specific gene gains with respect to D. melanogaster were noted. Sequence analyses also revealed that numerous repetitive sequences are a common structural feature of Drosophila and mosquito developmental genes. Finally, analysis of predicted miRNA binding sites in fruit fly and mosquito developmental genes suggests that the repertoire of developmental genes targeted by miRNAs is species-specific. The results of this study provide insight into the evolution of developmental genes and processes in dipterans and other arthropods, serve as a resource for those pursuing analysis of mosquito development, and will promote the design and refinement of functional analysis experiments.

A comprehensive gene expression atlas of sex- and tissue-specificity in the malaria vector, Anopheles gambiae

Background

The mosquito, Anopheles gambiae, is the primary vector of human malaria, a disease responsible for millions of deaths each year. To improve strategies for controlling transmission of the causative parasite, Plasmodium falciparum, we require a thorough understanding of the developmental mechanisms, physiological processes and evolutionary pressures affecting life-history traits in the mosquito. Identifying genes expressed in particular tissues or involved in specific biological processes is an essential part of this process.

Results

In this study, we present transcription profiles for ~82% of annotated Anopheles genes in dissected adult male and female tissues. The sensitivity afforded by examining dissected tissues found gene activity in an additional 20% of the genome that is undetected when using whole-animal samples. The somatic and reproductive tissues we examined each displayed patterns of sexually dimorphic and tissue-specific expression. By comparing expression profiles with Drosophila melanogaster we also assessed which genes are well conserved within the Diptera versus those that are more recently evolved.

Conclusions

Our expression atlas and associated publicly available database, the MozAtlas (http://www.tissue-atlas.org), provides information on the relative strength and specificity of gene expression in several somatic and reproductive tissues, isolated from a single strain grown under uniform conditions. The data will serve as a reference for other mosquito researchers by providing a simple method for identifying where genes are expressed in the adult, however, in addition our resource will also provide insights into the evolutionary diversity associated with gene expression levels among species.

The RNA-Seq approach to studying the expression of mosquito mitochondrial genes

In this study, we used extensive expressed sequence tag evidence obtained through 454 and Solexa next-generation sequencing to explore mtDNA transcription in male and female first instar larvae of Aedes aegypti and adults of Aedes aegypti, Anopheles gambiae, and Anopheles quadrimaculatus. Relative abundances of individual transcripts differed considerably within each sample, consistent with the differential stability of messenger RNA species. Large differences were also observed between species and between larval and adult stages; however, the male and female larval samples were remarkably similar. Quantitative PCR analysis of selected genes, cox1, l-rRNA and nd5, in larvae and adults of Ae. aegypti and in An. gambiae adults was consistent with the RNA-Seq-based quantification of expression. Finally, the absence of a conserved mtDNA region involved in transcriptional control in other dipterans suggests that mosquitoes have evolved a distinct mechanism of regulation of gene expression in the mitochondrion.

A phylogenetically based transcriptome age index mirrors ontogenetic divergence patterns

Parallels between phylogeny and ontogeny have been discussed for almost two centuries, and a number of theories have been proposed to explain such patterns. Especially elusive is the phylotypic stage, a phase during development where species within a phylum are particularly similar to each other. Although this has formerly been interpreted as a recapitulation of phylogeny, it is now thought to reflect an ontogenetic progression phase, where strong constraints on developmental regulation and gene interactions exist. Several studies have shown that genes expressed during this stage evolve at a slower rate, but it has so far not been possible to derive an unequivocal molecular signature associated with this stage. Here we use a combination of phylostratigraphy and stage-specific gene expression data to generate a cumulative index that reflects the evolutionary age of the transcriptome at given ontogenetic stages. Using zebrafish ontogeny and adult development as a model, we find that the phylotypic stage does indeed express the oldest transcriptome set and that younger sets are expressed during early and late development, thus faithfully mirroring the hourglass model of morphological divergence. Reproductively active animals show the youngest transcriptome, with major differences between males and females. Notably, ageing animals express increasingly older genes. Comparisons with similar data sets from flies and nematodes show that this pattern occurs across phyla. Our results indicate that an old transcriptome marks the phylotypic phase and that phylogenetic differences at other ontogenetic stages correlate with the expression of newly evolved genes.

Transcriptomic and functional analysis of the Anopheles gambiae salivary gland in relation to blood feeding

Background

The Anopheles gambiae salivary glands play a major role in malaria transmission and express a variety of bioactive components that facilitate blood-feeding by preventing platelet aggregation, blood clotting, vasodilatation, and inflammatory and other reactions at the probing site on the vertebrate host.

Results

We have performed a global transcriptome analysis of the A. gambiae salivary gland response to blood-feeding, to identify candidate genes that are involved in hematophagy. A total of 4,978 genes were found to be transcribed in this tissue. A comparison of salivary gland transcriptomes prior to and after blood-feeding identified 52 and 41 transcripts that were significantly up-regulated and down-regulated, respectively. Ten genes were further selected to assess their role in the blood-feeding process using RNAi-mediated gene silencing methodology. Depletion of the salivary gland genes encoding D7L2, anophelin, peroxidase, the SG2 precursor, and a 5'nucleotidase gene significantly increased probing time of A. gambiae mosquitoes and thereby their capacity to blood-feed.

Conclusions

The salivary gland transcriptome comprises approximately 38% of the total mosquito transcriptome and a small proportion of it is dynamically changing already at two hours in response to blood feeding. A better understanding of the salivary gland transcriptome and its function can contribute to the development of pathogen transmission control strategies and the identification of medically relevant bioactive compounds.

Evolutionary analysis of the kinesin light chain genes in the yellow fever mosquito Aedes aegypti: gene duplication as a source for novel early zygotic genes

Background

The maternal zygotic transition marks the time at which transcription from the zygotic genome is initiated and a subset of maternal RNAs are progressively degraded in the developing embryo. A number of early zygotic genes have been identified in Drosophila melanogaster and comparisons to sequenced mosquito genomes suggest that some of these early zygotic genes such as bottleneck are fast-evolving or subject to turnover in dipteran insects. One objective of this study is to identify early zygotic genes from the yellow fever mosquito Aedes aegypti to study their evolution. We are also interested in obtaining early zygotic promoters that will direct transgene expression in the early embryo as part of a Medea gene drive system.

Results

Two novel early zygotic kinesin light chain genes we call AaKLC2.1 and AaKLC2.2 were identified by transcriptome sequencing of Aedes aegypti embryos at various time points. These two genes have 98% nucleotide and amino acid identity in their coding regions and show transcription confined to the early zygotic stage according to gene-specific RT-PCR analysis. These AaKLC2 genes have a paralogous gene (AaKLC1) in Ae. aegypti. Phylogenetic inference shows that an ortholog to the AaKLC2 genes is only found in the sequenced genome of Culex quinquefasciatus. In contrast, AaKLC1 gene orthologs are found in all three sequenced mosquito species including Anopheles gambiae. There is only one KLC gene in D. melanogaster and other sequenced holometabolous insects that appears to be similar to AaKLC1. Unlike AaKLC2, AaKLC1 is expressed in all life stages and tissues tested, which is consistent with the expression pattern of the An. gambiae and D. melanogaster KLC genes. Phylogenetic inference also suggests that AaKLC2 genes and their likely C. quinquefasciatus ortholog are fast-evolving genes relative to the highly conserved AaKLC1-like paralogs. Embryonic injection of a luciferase reporter under the control of a 1 kb fragment upstream of the AaKLC2.1 start codon shows promoter activity at least as early as 3 hours in the developing Ae. aegypti embryo. The AaKLC2.1 promoter activity reached ~1600 fold over the negative control at 5 hr after egg deposition.

Conclusions

Transcriptome profiling by use of high throughput sequencing technologies has proven to be a valuable method for the identification and discovery of early and transient zygotic genes. The evolutionary investigation of the KLC gene family reveals that duplication is a source for the evolution of new genes that play a role in the dynamic process of early embryonic development. AaKLC2.1 may provide a promoter for early zygotic-specific transgene expression, which is a key component of the Medea gene drive system.

Developmental plasticity and the evolution of animal complex life cycles

Metazoan life cycles can be complex in different ways. A number of diverse phenotypes and reproductive events can sequentially occur along the cycle, and at certain stages a variety of developmental and reproductive options can be available to the animal, the choice among which depends on a combination of organismal and environmental conditions. We hypothesize that a diversity of phenotypes arranged in developmental sequence throughout an animal's life cycle may have evolved by genetic assimilation of alternative phenotypes originally triggered by environmental cues. This is supported by similarities between the developmental mechanisms mediating phenotype change and alternative phenotype determination during ontogeny and the common ecological condition that favour both forms of phenotypic variation. The comparison of transcription profiles from different developmental stages throughout a complex life cycle with those from alternative phenotypes in closely related polyphenic animals is expected to offer critical evidence upon which to evaluate our hypothesis.

Analysis of expression in the Anopheles gambiae developing testes reveals rapidly evolving lineage-specific genes in mosquitoes

Background

Male mosquitoes do not feed on blood and are not involved in delivery of pathogens to humans. Consequently, they are seldom the subjects of research, which results in a very poor understanding of their biology. To gain insights into male developmental processes we sought to identify genes transcribed exclusively in the reproductive tissues of male Anopheles gambiae pupae.

Results

Using a cDNA subtraction strategy, five male-specifically or highly male-biased expressed genes were isolated, four of which remain unannotated in the An. gambiae genome. Spatial and temporal expression patterns suggest that each of these genes is involved in the mid-late stages of spermatogenesis. Their sequences are rapidly evolving; however, two genes possess clear homologs in a wide range of taxa and one of these probably acts in a sperm motility control mechanism conserved in many organisms, including humans. The other three genes have no match to sequences from non-mosquito taxa, thus can be regarded as orphans. RNA in situ hybridization demonstrated that one of the orphans is transcribed in spermatids, which suggests its involvement in sperm maturation. Two other orphans have unknown functions. Expression analysis of orthologs of all five genes indicated that male-biased transcription was not conserved in the majority of cases in Aedes and Culex.

Conclusion

Discovery of testis-expressed orphan genes in mosquitoes opens new prospects for the development of innovative control methods. The orphan encoded proteins may represent unique targets of selective anti-mosquito sterilizing agents that will not affect non-target organisms.

The cellular geometry of growth drives the amino acid economy of Caenorhabditis elegans

The nematode Caenorhabditis elegans grows largely by increases in cell size. As a consequence of this, the surface: volume ratio of its cells must decline in the course of postembryonic growth. Here we use transcriptomic and metabolomic data to show that this change in geometry can explain a variety of phenomena during growth, including: (i) changes in the relative expression levels of cytoplasmic and membrane proteins; (ii) changes in the relative usage of the twenty amino acids in expressed proteins, as estimated by changes in the transcriptome; and (iii) changes in metabolite pools of free amino acids. We expect these relations to be universal in single cells and in whole multicellular organisms that grow largely by increases in cell size, but not those that grow by cell proliferation.