VectorBase: a home for invertebrate vectors of human pathogens

Dishevelled Has Anti-Viral Activity in Rift Valley Fever Virus Infected Aedes aegypti

Mosquitoes in the genera Aedes and Culex are vectors of Rift Valley fever virus (RVFV), which emerges in periodic epidemics in Africa and Saudi Arabia. Factors that influence the transmission dynamics of RVFV are not well characterized. To address this, we interrogated mosquito host-signaling responses through analysis of differentially expressed genes (DEGs) in two mosquito species with marked differences in RVFV vector competence: Aedes aegypti (Aae, low competence) and Culex tarsalis (Cxt, high competence). Mosquito-host transcripts related to three different signaling pathways were investigated. Selected genes from the Wingless (Wg, WNT-beta-catenin) pathway, which is a conserved regulator of cell proliferation and differentiation, were assessed. One of these, dishevelled (DSH), differentially regulates progression/inhibition of the WNT and JNK (c-Jun N-terminal Kinase) pathways. A negative regulator of the JNK-signaling pathway, puckered, was also assessed. Lastly, Janus kinase/signal transducers and activators of transcription (JAK-STAT) are important for innate immunity; in this context, we tested domeless levels. Here, individual Aae and Cxt were exposed to RVFV MP-12 via oral bloodmeals and held for 14 days. Robust decreases in DEGs in both Aae and Cxt were observed. In particular, Aae DSH expression, but not Cxt DSH, was correlated to the presence/absence of viral RNA at 14 days post-challenge (dpc). Moreover, there was an inverse relationship between the viral copy number and aaeDSH expression. DSH silencing resulted in increased viral copy numbers compared to controls at 3 dpc, consistent with a role for aaeDSH in antiviral immunity. Analysis of cis-regulatory regions for the genes of interest revealed clues to upstream regulation of these pathways.

The Trypanosoma brucei MISP family of invariant proteins is co-expressed with BARP as triple helical bundle structures on the surface of salivary gland forms, but is dispensable for parasite development within the tsetse vector

Trypanosoma brucei spp. develop into mammalian-infectious metacyclic trypomastigotes inside tsetse salivary glands. Besides acquiring a variant surface glycoprotein (VSG) coat, little is known about the metacyclic expression of invariant surface antigens. Proteomic analyses of saliva from T. brucei-infected tsetse flies identified, in addition to VSG and Brucei Alanine-Rich Protein (BARP) peptides, a family of glycosylphosphatidylinositol (GPI)-anchored surface proteins herein named as Metacyclic Invariant Surface Proteins (MISP) because of its predominant expression on the surface of metacyclic trypomastigotes. The MISP family is encoded by five paralog genes with >80% protein identity, which are exclusively expressed by salivary gland stages of the parasite and peak in metacyclic stage, as shown by confocal microscopy and immuno-high resolution scanning electron microscopy. Crystallographic analysis of a MISP isoform (MISP360) and a high confidence model of BARP revealed a triple helical bundle architecture commonly found in other trypanosome surface proteins. Molecular modelling combined with live fluorescent microscopy suggests that MISP N-termini are potentially extended above the metacyclic VSG coat, and thus could be tested as a transmission-blocking vaccine target. However, vaccination with recombinant MISP360 isoform did not protect mice against a T. brucei infectious tsetse bite. Lastly, both CRISPR-Cas9-driven knock out and RNAi knock down of all MISP paralogues suggest they are not essential for parasite development in the tsetse vector. We suggest MISP may be relevant during trypanosome transmission or establishment in the vertebrate's skin.

Spatiotemporal Environmental Drivers of Eastern Equine Encephalitis Virus in Central Florida: Towards a Predictive Model for a Lethal Disease

Eastern equine encephalitis virus (EEEV) is a mosquito-borne virus that affects humans and horses, with a high case fatality rate in both species. The virus can be transmitted by several mosquito species and maintained in multiple reservoir hosts, making EEEV dynamics difficult to anticipate. In this study, we identified spatial and temporal factors that precede EEEV detection using sentinel chicken surveillance data from Orange County, Florida, from 2003 to 2017. We first examined the land cover and mosquito species composition associated with sentinel chicken sites. We then fit distributed lag nonlinear models of EEEV detection at the county scale, using monthly temperature, precipitation, and Southern Oscillation Index values, and at the sentinel flock-scale, using remotely sensed temperature and wetness indicators. We found positive associations between the percent wooded wetlands and the count of EEEV detections. We found Culiseta melanura (Diptera: Culicidae) were more abundant at positive sites in winter and summer, but Coquillettidia perturbans (Walker) were more abundant at positive sites in spring. In the county-wide model, precipitation, temperature, and Southern Oscillation Index values at lags of two, nine, and twelve months were significant, respectively, while temperature and wetness were significant at lags of eight and six months in the flock-specific models.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Aedes aegypti miRNA-33 modulates permethrin induced toxicity by regulating VGSC transcripts

Aedes aegypti is a major vector of Zika, dengue, and other arboviruses. Permethrin adulticidal spraying, which targets the voltage-gated sodium channel (VGSC), is commonly done to reduce local mosquito populations and protect humans from exposure to arbovirus pathogens transmitted by this dangerous pest. Permethrin resistance, however, is a growing problem and understanding its underlying molecular basis may identify avenues to combat it. We identified a single G:C polymorphism in pre-miR-33 that was genetically associated with permethrin resistance; resulting isoforms had structural differences that may affect DICER-1/pre-miRNA processing rates. We then assessed the effects of overexpression of pre-miR-33 isoforms on permethrin toxicological phenotypes, VGSC transcript abundance and protein levels for two genetically related mosquito strains. One strain had its naturally high permethrin resistance levels maintained by periodic treatment, and the other was released from selection. VGSC protein levels were lower in the permethrin resistant strain than in the related permethrin-susceptible strain. Overexpression of the G-pre-miR-33 isoform reduced VGSC expression levels in both strains. To further elucidate changes in gene expression associated with permethrin resistance, exome-capture gDNA deep sequencing, genetic association mapping and subsequent gene set enrichment analysis revealed that transport genes, in particular, were selected in resistant versus susceptible mosquitoes. Collectively, these data indicate that miR-33 regulates VGSC expression as part of a nuanced system of neuronal regulation that contributes to a network of heritable features determining permethrin resistance.

ASDB: A comprehensive omics database for Anopheles sinensis

Anopheles sinensis is a key disease vector for human malaria and parasitic diseases such as malayan filariasis, and it is considered to be one of the most important malaria vectors in China and Southeast Asia. As high-throughput sequencing and assembly technology are widely used in An. sinensis, a lot of omics data have been generated, and abundant genome, mRNA transcriptome, miRNA transcriptome and resequencing results have been accumulated. In addition, lots of valuable morphological images and publications have been produced with the in-depth studies on An. sinensis. However, the increased quantity, variety, and structure complexity of the omics data create inconveniences for researchers to use and manage this information. We have built an An. sinensis omics database (ASDB, http://asdb.jungleran.com/) - a comprehensive and integrated database to promote scientific research on An. sinensis. Docker was used to deploy a development environment and Drupal to build ASDB. ASDB provides a Blast tool to do sequence alignment of genome sequence, gene sequence and protein sequence of An. sinensis. It also offers JBrowse (a next-generation genome visualization and analysis web platform) to facilitate researchers visualize the gene structure, non-coding RNA (include miRNA, snRNA, tRNA and so on) structure and genomic variation sites as desired. ASDB has integrated various latest omics data of An. Sinensis, including de novo genome and its annotation data, genome variation data (such as SNP and InDel), transcriptome and its expression value, miRNA expression value and miRNA-mRNA interaction, metagenomes. The database has also included the morphological images of different developmental stages and tissues, and important literatures associated with An. sinensis. ASDB provides a user-friendly search and displays pages. The integration of these resources will contribute to the study of basic biology and functional genome of An. sinensis.

Gene co-expression network analysis of Trypanosoma brucei in tsetse fly vector

Background

Trypanosoma brucei species are motile protozoan parasites that are cyclically transmitted by tsetse fly (genus Glossina) causing human sleeping sickness and nagana in livestock in sub-Saharan Africa. African trypanosomes display digenetic life cycle stages in the tsetse fly vector and in their mammalian host. Experimental work on insect-stage trypanosomes is challenging because of the difficulty in setting up successful in vitro cultures. Therefore, there is limited knowledge on the trypanosome biology during its development in the tsetse fly. Consequently, this limits the development of new strategies for blocking parasite transmission in the tsetse fly.

Methods

In this study, RNA-Seq data of insect-stage trypanosomes were used to construct a T. brucei gene co-expression network using the weighted gene co-expression analysis (WGCNA) method. The study identified significant enriched modules for genes that play key roles during the parasite’s development in tsetse fly. Furthermore, potential 3′ untranslated region (UTR) regulatory elements for genes that clustered in the same module were identified using the Finding Informative Regulatory Elements (FIRE) tool.

Results

A fraction of gene modules (12 out of 27 modules) in the constructed network were found to be enriched in functional roles associated with the cell division, protein biosynthesis, mitochondrion, and cell surface. Additionally, 12 hub genes encoding proteins such as RNA-binding protein 6 (RBP6), arginine kinase 1 (AK1), brucei alanine-rich protein (BARP), among others, were identified for the 12 significantly enriched gene modules. In addition, the potential regulatory elements located in the 3′ untranslated regions of genes within the same module were predicted.

Conclusions

The constructed gene co-expression network provides a useful resource for network-based data mining to identify candidate genes for functional studies. This will enhance understanding of the molecular mechanisms that underlie important biological processes during parasite’s development in tsetse fly. Ultimately, these findings will be key in the identification of potential molecular targets for disease control.

Vgsc-interacting proteins are genetically associated with pyrethroid resistance in Aedes aegypti

Association mapping of factors that condition pyrethroid resistance in Aedes aegypti has consistently identified genes in multiple functional groups. Toward better understanding of the mechanisms involved, we examined high throughput sequencing data (HTS) from two Aedes aegypti aegypti collections from Merida, Yucatan, Mexico treated with either permethrin or deltamethrin. Exome capture enrichment for coding regions and the AaegL5 annotation were used to identify genes statistically associated with resistance. The frequencies of single nucleotide polymorphisms (SNPs) were compared between resistant and susceptible mosquito pools using a contingency χ² analysis. The -log₁₀(χ²p value) was calculated at each SNP site, with a weighted average determined from all sites in each gene. Genes with -log₁₀(χ²p value) ≥ 4.0 and present among all 3 treatment groups were subjected to gene set enrichment analysis (GSEA). We found that several functional groups were enriched compared to all coding genes. These categories were transport, signal transduction and metabolism, in order from highest to lowest statistical significance. Strikingly, 21 genes with demonstrated association to synaptic function were identified. In the high association group (n = 1,053 genes), several genes were identified that also genetically or physically interact with the voltage-gated sodium channel (VGSC). These genes were eg., CHARLATAN (CHL), a transcriptional regulator, several ankyrin-domain proteins, PUMILIO (PUM), a translational repressor, and NEDD4 (E3 ubiquitin-protein ligase). There were 13 genes that ranked among the top 10%: these included VGSC; CINGULIN, a predicted neuronal gap junction protein, and the aedine ortholog of NERVY (NVY), a transcriptional regulator. Silencing of CHL and NVY followed by standard permethrin bottle bioassays validated their association with permethrin resistance. Importantly, VGSC levels were also reduced about 50% in chl- or nvy-dsRNA treated mosquitoes. These results are consistent with the contribution of a variety of neuronal pathways to pyrethroid resistance in Ae. aegypti.

Symbiotic microbes affect the expression of male reproductive genes in Glossina m. morsitans

Background

Tsetse flies (Diptera, Glossinidae) display unique reproductive biology traits. Females reproduce through adenotrophic viviparity, nourishing the growing larva into their modified uterus until parturition. Males transfer their sperm and seminal fluid, produced by both testes and male accessory glands, in a spermatophore capsule transiently formed within the female reproductive tract upon mating. Both sexes are obligate blood feeders and have evolved tight relationships with endosymbionts, already shown to provide essential nutrients lacking in their diet. However, the partnership between tsetse and its symbionts has so far been investigated, at the molecular, genomic and metabolomics level, only in females, whereas the roles of microbiota in male reproduction are still unexplored.

Results

Here we begin unravelling the impact of microbiota on Glossina m. morsitans (G. morsitans) male reproductive biology by generating transcriptomes from the reproductive tissues of males deprived of their endosymbionts (aposymbiotic) via maternal antibiotic treatment and dietary supplementation. We then compared the transcriptional profiles of genes expressed in the male reproductive tract of normal and these aposymbiotic flies. We showed that microbiota removal impacts several male reproductive genes by depressing the activity of genes in the male accessory glands (MAGs), including sequences encoding seminal fluid proteins, and increasing expression of genes in the testes. In the MAGs, in particular, the expression of genes related to mating, immunity and seminal fluid components’ synthesis is reduced. In the testes, the absence of symbionts activates genes involved in the metabolic apparatus at the basis of male reproduction, including sperm production, motility and function.

Conclusions

Our findings mirrored the complementary roles male accessory glands and testes play in supporting male reproduction and open new avenues for disentangling the interplay between male insects and endosymbionts. From an applied perspective, unravelling the metabolic and functional relationships between tsetse symbionts and male reproductive physiology will provide fundamental information useful to understanding the biology underlying improved male reproductive success in tsetse. This information is of particular importance in the context of tsetse population control via Sterile Insect Technique (SIT) and its impact on trypanosomiasis transmission.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1289-2) contains supplementary material, which is available to authorized users.

Unbiased classification of mosquito blood cells by single-cell genomics and high-content imaging

Mosquito blood cells are immune cells that help control infection by vector-borne pathogens. Despite their importance, little is known about mosquito blood cell biology beyond morphological and functional criteria used for their classification. Here, we combined the power of single-cell RNA sequencing, high-content imaging flow cytometry, and single-molecule RNA hybridization to analyze a subset of blood cells of the malaria mosquito Anopheles gambiae By demonstrating that blood cells express nearly half of the mosquito transcriptome, our dataset represents an unprecedented view into their transcriptional program. Analyses of differentially expressed genes identified transcriptional signatures of two cell types and provide insights into the current classification of these cells. We further demonstrate the active transfer of a cellular marker between blood cells that may confound their identification. We propose that cell-to-cell exchange may contribute to cellular diversity and functional plasticity seen across biological systems.

Unravelling the relationship between the tsetse fly and its obligate symbiont Wigglesworthia: transcriptomic and metabolomic landscapes reveal highly integrated physiological networks

Insects with restricted diets rely on obligate microbes to fulfil nutritional requirements essential for biological function. Tsetse flies, vectors of African trypanosome parasites, feed exclusively on vertebrate blood and harbour the obligate endosymbiont Wigglesworthia glossinidia. Without Wigglesworthia, tsetse are unable to reproduce. These symbionts are sheltered within specialized cells (bacteriocytes) that form the midgut-associated bacteriome organ. To decipher the core functions of this symbiosis essential for tsetse's survival, we performed dual-RNA-seq analysis of the bacteriome, coupled with metabolomic analysis of bacteriome and haemolymph collected from normal and symbiont-cured (sterile) females. Bacteriocytes produce immune regulatory peptidoglycan recognition protein (pgrp-lb) that protects Wigglesworthia, and a multivitamin transporter (smvt) that can aid in nutrient dissemination. Wigglesworthia overexpress a molecular chaperone (GroEL) to augment their translational/transport machinery and biosynthesize an abundance of B vitamins (specifically B₁-, B₂-, B₃- and B₆-associated metabolites) to supplement the host's nutritionally deficient diet. The absence of Wigglesworthia's contributions disrupts multiple metabolic pathways impacting carbohydrate and amino acid metabolism. These disruptions affect the dependent downstream processes of nucleotide biosynthesis and metabolism and biosynthesis of S-adenosyl methionine (SAM), an essential cofactor. This holistic fundamental knowledge of the symbiotic dialogue highlights new biological targets for the development of innovative vector control methods.

Alternative patterns of sex chromosome differentiation in Aedes aegypti (L)

Background

Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex F_ST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings.

Results

Intersex standardized variance (F_ST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with F_ST values ≥ 0.100 and/or with increased expected heterozygosity (H_exp, one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex F_ST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex F_ST genes proximal to the M locus, whereas SenAae had inter-sex F_ST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased H_exp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation.

Conclusions

The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae.

Electronic supplementary material

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

Identification of cis-regulatory sequences reveals potential participation of lola and Deaf1 transcription factors in Anopheles gambiae innate immune response

The innate immune response of Anopheles gambiae involves the transcriptional upregulation of effector genes. Therefore, the cis-regulatory sequences and their cognate binding factors play essential roles in the mosquito’s immune response. However, the genetic control of the mosquito’s innate immune response is not yet fully understood. To gain further insight on the elements, the factors and the potential mechanisms involved, an open chromatin profiling was carried out on A. gambiae-derived immune-responsive cells. Here, we report the identification of cis-regulatory sites, immunity-related transcription factor binding sites, and cis-regulatory modules. A de novo motif discovery carried out on this set of cis-regulatory sequences identified immunity-related motifs and cis-regulatory modules. These modules contain motifs that are similar to binding sites for REL-, STAT-, lola- and Deaf1-type transcription factors. Sequence motifs similar to the binding sites for GAGA were found within a cis-regulatory module, together with immunity-related transcription factor binding sites. The presence of Deaf1- and lola-type binding sites, along with REL- and STAT-type binding sites, suggests that the immunity function of these two factors could have been conserved both in Drosophila and Anopheles gambiae.

Exon-Enriched Libraries Reveal Large Genic Differences Between Aedes aegypti from Senegal, West Africa, and Populations Outside Africa

Aedes aegypti is one of the most studied mosquito species, and the principal vector of several arboviruses pathogenic to humans. Recently failure to oviposit, low fecundity, and poor egg-to-adult survival were observed when Ae. aegypti from Senegal (SenAae) West Africa were crossed with Ae. aegypti (Aaa) from outside of Africa, and in SenAae intercrosses. Fluorescent in situ hybridization analyses indicated rearrangements on chromosome 1, and pericentric inversions on chromosomes 2 and 3. Herein, high throughput sequencing (HTS) of exon-enriched libraries was used to compare chromosome-wide genetic diversity among Aaa collections from rural Thailand and Mexico, a sylvatic collection from southeastern Senegal (PK10), and an urban collection from western Senegal (Kaolack). Sex-specific polymorphisms were analyzed in Thailand and PK10 to assess genetic differences between sexes. Expected heterozygosity was greatest in SenAae F_ST distributions of 15,735 genes among all six pairwise comparisons of the four collections indicated that Mexican and Thailand collections are genetically similar, while F_ST distributions between PK10 and Kaolack were distinct. All four comparisons of SenAae with Aaa indicated extreme differentiation. F_ST was uniform between sexes across all chromosomes in Thailand, but were different, especially on the sex autosome 1, in PK10. These patterns correlate with the reproductive isolation noted earlier. We hypothesize that cryptic Ae. aegypti taxa may exist in West Africa, and the large genic differences between Aaa and SenAae detected in the present study have accumulated over a long period following the evolution of chromosome rearrangements in allopatric populations that subsequently cause reproductive isolation when these populations became sympatric.

Protein Bioinformatics Databases and Resources

Many publicly available data repositories and resources have been developed to support protein-related information management, data-driven hypothesis generation, and biological knowledge discovery. To help researchers quickly find the appropriate protein-related informatics resources, we present a comprehensive review (with categorization and description) of major protein bioinformatics databases in this chapter. We also discuss the challenges and opportunities for developing next-generation protein bioinformatics databases and resources to support data integration and data analytics in the Big Data era.

A review on host-pathogen interactions: classification and prediction

The research on host-pathogen interactions is an ever-emerging and evolving field. Every other day a new pathogen gets discovered, along with comes the challenge of its prevention and cure. As the intelligent human always vies for prevention, which is better than cure, understanding the mechanisms of host-pathogen interactions gets prior importance. There are many mechanisms involved from the pathogen as well as the host sides while an interaction happens. It is a vis-a-vis fight of the counter genes and proteins from both sides. Who wins depends on whether a host gets an infection or not. Moreover, a higher level of complexity arises when the pathogens evolve and become resistant to a host's defense mechanisms. Such pathogens pose serious challenges for treatment. The entire human population is in danger of such long-lasting persistent infections. Some of these infections even increase the rate of mortality. Hence there is an immediate emergency to understand how the pathogens interact with their host for successful invasion. It may lead to discovery of appropriate preventive measures, and the development of rational therapeutic measures and medication against such infections and diseases. This review, a state-of-the-art updated scenario of host-pathogen interaction research, has been done by keeping in mind this urgency. It covers the biological and computational aspects of host-pathogen interactions, classification of the methods by which the pathogens interact with their hosts, different machine learning techniques for prediction of host-pathogen interactions, and future scopes of this research field.

A generic schema and data collection forms applicable to diverse entomological studies of mosquitoes

BACKGROUND

Standardized schemas, databases, and public data repositories are needed for the studies of malaria vectors that encompass a remarkably diverse array of designs and rapidly generate large data volumes, often in resource-limited tropical settings lacking specialized software or informatics support.

RESULTS

Data from the majority of mosquito studies conformed to a generic schema, with data collection forms recording the experimental design, sorting of collections, details of sample pooling or subdivision, and additional observations. Generically applicable forms with standardized attribute definitions enabled rigorous, consistent data and sample management with generic software and minimal expertise. Forms use now includes 20 experiments, 8 projects, and 15 users at 3 research and control institutes in 3 African countries, resulting in 11 peer-reviewed publications.

CONCLUSION

We have designed generic data schema that can be used to develop paper or electronic based data collection forms depending on the availability of resources. We have developed paper-based data collection forms that can be used to collect data from majority of entomological studies across multiple study areas using standardized data formats. Data recorded on these forms with standardized formats can be entered and linked with any relational database software. These informatics tools are recommended because they ensure that medical entomologists save time, improve data quality, and data collected and shared across multiple studies is in standardized formats hence increasing research outputs.

Silencing the odorant receptor co-receptor RproOrco affects the physiology and behavior of the Chagas disease vector Rhodnius prolixus

Olfaction is one of the main sensory modalities that allow insects to interpret their environment. Several proteins, including odorant-binding proteins (OBPs) and odorant receptors (ORs), are involved in this process. Odorant receptors are ion channels formed by a binding unit OR and an odorant receptor co-receptor (Orco). The main goal of this study was to characterize the Orco gene of Rhodnius prolixus (RproOrco) and to infer its biological functions using gene silencing. The full-length RproOrco gene sequence was downloaded from VectorBase. This gene has 7 introns and is located in the genome SuperContig GL563069: 1,017,713-1,023,165. RproOrco encodes a protein of 473 amino acids, with predicted 7 transmembrane domains, and is highly expressed in the antennae during all R. prolixus developmental stages. The RNAi technique effectively silenced RproOrco, reducing the gene's expression by approximately 73%. Interestingly, the effect of gene silencing persisted for more than 100 days, indicating a prolonged effect of dsRNA that was maintained even after molting. The phenotypic effects of silencing involved the following: (1) loss of the ability to find a vertebrate host in a timely manner, (2) decreased ingested blood volume, (3) delayed and decreased molt rate, (4) increased mortality rate, and (5) decreased egg laying. Our data strongly suggest that dsOrco disrupts R. prolixus host-finding behavior, which is further reflected in the blood ingestion, molting, mortality, and egg laying data. This study clearly demonstrates that Orco is an excellent target for controlling triatomine populations. Thus, the data presented here open new possibilities for the control of vector-borne diseases.

"Salivary gland cellular architecture in the Asian malaria vector mosquito Anopheles stephensi"

Background

Anopheles mosquitoes are vectors for malaria, a disease with continued grave outcomes for human health. Transmission of malaria from mosquitoes to humans occurs by parasite passage through the salivary glands (SGs). Previous studies of mosquito SG architecture have been limited in scope and detail.

Methods

We developed a simple, optimized protocol for fluorescence staining using dyes and/or antibodies to interrogate cellular architecture in Anopheles stephensi adult SGs. We used common biological dyes, antibodies to well-conserved structural and organellar markers, and antibodies against Anopheles salivary proteins to visualize many individual SGs at high resolution by confocal microscopy.

Results

These analyses confirmed morphological features previously described using electron microscopy and uncovered a high degree of individual variation in SG structure. Our studies provide evidence for two alternative models for the origin of the salivary duct, the structure facilitating parasite transport out of SGs. We compare SG cellular architecture in An. stephensi and Drosophila melanogaster, a fellow Dipteran whose adult SGs are nearly completely unstudied, and find many conserved features despite divergence in overall form and function. Anopheles salivary proteins previously observed at the basement membrane were localized either in SG cells, secretory cavities, or the SG lumen. Our studies also revealed a population of cells with characteristics consistent with regenerative cells, similar to muscle satellite cells or midgut regenerative cells.

Conclusions

This work serves as a foundation for linking Anopheles stephensi SG cellular architecture to function and as a basis for generating and evaluating tools aimed at preventing malaria transmission at the level of mosquito SGs.

Electronic supplementary material

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

Dengue virus models based on mice as experimental hosts

Dengue virus (DENV) causes dengue fever, a widely distributed endemic disease transmitted by mosquitoes. The complex interaction of DENV with the human immune system has complicated the development of an effective vaccine. This may be attributed, at least in part, to the lack of a suitable animal model capable to reproduce symptoms observed in humans. Mouse models are simple but usually rely on host-adapted virus strains or immunodeficient mouse lineages. Recent evidences indicated that some natural DENV strains are capable to infect immunocompetent mice. In addition, humanized mouse lineages can more faithfully reproduce some of the symptoms observed in humans. Such experimental models are valuable tools for the study of DENV biology.

Anopheles gambiae salivary protein expression modulated by wild Plasmodium falciparum infection: highlighting of new antigenic peptides as candidates of An. gambiae bites

Background

Malaria is the major parasitic disease worldwide caused by Plasmodium infection. The objective of integrated malaria control programs is to decrease malaria transmission, which needs specific tools to be accurately assessed. In areas where the transmission is low or has been substantially reduced, new complementary tools have to be developed to improve surveillance. A recent approach, based on the human antibody response to Anopheles salivary proteins, has been shown to be efficient in evaluating human exposure to Anopheles bites. The aim of the present study was to identify new An. gambiae salivary proteins as potential candidate biomarkers of human exposure to P. falciparum-infective bites.

Methods

Experimental infections of An. gambiae by wild P. falciparum were carried out in semi-field conditions. Then a proteomic approach, combining 2D-DIGE and mass spectrometry, was used to identify the overexpressed salivary proteins in infected salivary glands compared to uninfected An. gambiae controls. Subsequently, a peptide design of each potential candidate was performed in silico and their antigenicity was tested by an epitope-mapping technique using blood from individuals exposed to Anopheles bites.

Results

Five salivary proteins (gSG6, gSG1b, TRIO, SG5 and long form D7) were overexpressed in the infected salivary glands. Eighteen peptides were designed from these proteins and were found antigenic in children exposed to the Anopheles bites. Moreover, the results showed that the presence of wild P. falciparum in salivary glands modulates the expression of several salivary proteins and also appeared to induce post-translational modifications.

Conclusions

This study is, to our knowledge, the first that compares the sialome of An. gambiae both infected and not infected by wild P. falciparum, making it possible to mimic the natural conditions of infection. This is a first step toward a better understanding of the close interactions between the parasite and the salivary gland of mosquitoes. In addition, these results open the way to define biomarkers of infective bites of Anopheles, which could, in the future, improve the estimation of malaria transmission and the evaluation of malaria vector control tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-014-0599-y) contains supplementary material, which is available to authorized users.

iPathCons and iPathDB: an improved insect pathway construction tool and the database

Insects are one of the most successful animal groups on earth. Some insects, such as the silkworm and honeybee, are beneficial to humans, whereas others are notorious pests of crops. At present, the genomes of 38 insects have been sequenced and made publically available. In addition, the transcriptomes of dozens of insects have been sequenced. As gene data rapidly accumulate, constructing the pathway of molecular interactions becomes increasingly important for entomological research. Here, we developed an improved tool, iPathCons, for knowledge-based construction of pathways from the transcriptomes or the official gene sets of genomes. Considering the high evolution diversity in insects, iPathCons uses a voting system for Kyoto Encyclopedia of Genes and Genomes Orthology assignment. Both stand-alone software and a web server of iPathCons are provided. Using iPathCons, we constructed the pathways of molecular interactions of 52 insects, including 37 genome-sequenced and 15 transcriptome-sequenced ones. These pathways are available in the iPathDB, which provides searches, web server, data downloads, etc. This database will be highly useful for the insect research community.

Database URL:http://ento.njau.edu.cn/ipath/

Dual RNA-seq of parasite and host reveals gene expression dynamics during filarial worm-mosquito interactions

BACKGROUND

Parasite biology, by its very nature, cannot be understood without integrating it with that of the host, nor can the host response be adequately explained without considering the activity of the parasite. However, due to experimental limitations, molecular studies of parasite-host systems have been predominantly one-sided investigations focusing on either of the partners involved. Here, we conducted a dual RNA-seq time course analysis of filarial worm parasite and host mosquito to better understand the parasite processes underlying development in and interaction with the host tissue, from the establishment of infection to the development of infective-stage larva.

METHODOLOGY/PRINCIPAL FINDINGS

Using the Brugia malayi-Aedes aegypti system, we report parasite gene transcription dynamics, which exhibited a highly ordered developmental program consisting of a series of cyclical and state-transitioning temporal patterns. In addition, we contextualized these parasite data in relation to the concurrent dynamics of the host transcriptome. Comparative analyses using uninfected tissues and different host strains revealed the influence of parasite development on host gene transcription as well as the influence of the host environment on parasite gene transcription. We also critically evaluated the life-cycle transcriptome of B. malayi by comparing developmental stages in the mosquito relative to those in the mammalian host, providing insight into gene expression changes underpinning the mosquito-borne parasitic lifestyle of this heteroxenous parasite.

CONCLUSIONS/SIGNIFICANCE

The data presented herein provide the research community with information to design wet lab experiments and select candidates for future study to more fully dissect the whole set of molecular interactions of both organisms in this mosquito-filarial worm symbiotic relationship. Furthermore, characterization of the transcriptional program over the complete life cycle of the parasite, including stages within the mosquito, could help devise novel targets for control strategies.

A novel highly divergent protein family identified from a viviparous insect by RNA-seq analysis: a potential target for tsetse fly-specific abortifacients

In tsetse flies, nutrients for intrauterine larval development are synthesized by the modified accessory gland (milk gland) and provided in mother's milk during lactation. Interference with at least two milk proteins has been shown to extend larval development and reduce fecundity. The goal of this study was to perform a comprehensive characterization of tsetse milk proteins using lactation-specific transcriptome/milk proteome analyses and to define functional role(s) for the milk proteins during lactation. Differential analysis of RNA-seq data from lactating and dry (non-lactating) females revealed enrichment of transcripts coding for protein synthesis machinery, lipid metabolism and secretory proteins during lactation. Among the genes induced during lactation were those encoding the previously identified milk proteins (milk gland proteins 1-3, transferrin and acid sphingomyelinase 1) and seven new genes (mgp4-10). The genes encoding mgp2-10 are organized on a 40 kb syntenic block in the tsetse genome, have similar exon-intron arrangements, and share regions of amino acid sequence similarity. Expression of mgp2-10 is female-specific and high during milk secretion. While knockdown of a single mgp failed to reduce fecundity, simultaneous knockdown of multiple variants reduced milk protein levels and lowered fecundity. The genomic localization, gene structure similarities, and functional redundancy of MGP2-10 suggest that they constitute a novel highly divergent protein family. Our data indicates that MGP2-10 function both as the primary amino acid resource for the developing larva and in the maintenance of milk homeostasis, similar to the function of the mammalian casein family of milk proteins. This study underscores the dynamic nature of the lactation cycle and identifies a novel family of lactation-specific proteins, unique to Glossina sp., that are essential to larval development. The specificity of MGP2-10 to tsetse and their critical role during lactation suggests that these proteins may be an excellent target for tsetse-specific population control approaches.

Silencing of P-glycoprotein increases mortality in temephos-treated Aedes aegypti larvae

Re-emergence of vector-borne diseases such as dengue and yellow fever, which are both transmitted by the Aedes aegypti mosquito, has been correlated with insecticide resistance. P-glycoproteins (P-gps) are ATP-dependent efflux pumps that are involved in the transport of substrates across membranes. Some of these proteins have been implicated in multidrug resistance (MDR). In this study, we identified a putative P-glycoprotein in the Ae. aegypti database based on its significantly high identity with Anopheles gambiae, Culex quinquefasciatus, Drosophila melanogaster and human P-gps. The basal ATPase activity of ATP-binding cassette transporters in larvae was significantly increased in the presence of MDR modulators (verapamil and quinidine). An eightfold increase in Ae. aegypti P-gp (AaegP-gp) gene expression was detected in temephos-treated larvae as determined by quantitative PCR. To analyse the potential role of AaegP-gp in insecticide efflux, a temephos larvicide assay was performed in the presence of verapamil. The results showed an increase of 24% in temephos toxicity, which is in agreement with the efflux reversing effect. RNA interference (RNAi)-mediated silencing of the AaegP-gp gene caused a significant increase in temephos toxicity (57%). In conclusion, we have demonstrated for the first time in insects that insecticide-induced P-gp expression can be involved in the modulation of insecticide efflux.

Transgenic mosquitoes expressing a phospholipase A(2) gene have a fitness advantage when fed Plasmodium falciparum-infected blood

Background

Genetically modified mosquitoes have been proposed as an alternative strategy to reduce the heavy burden of malaria. In recent years, several proof-of-principle experiments have been performed that validate the idea that mosquitoes can be genetically modified to become refractory to malaria parasite development.

Results

We have created two transgenic lines of Anophelesstephensi, a natural vector of Plasmodium falciparum, which constitutively secrete a catalytically inactive phospholipase A₂ (mPLA₂) into the midgut lumen to interfere with Plasmodium ookinete invasion. Our experiments show that both transgenic lines expressing mPLA₂ significantly impair the development of rodent malaria parasites, but only one line impairs the development of human malaria parasites. In addition, when fed on malaria-infected blood, mosquitoes from both transgenic lines are more fecund than non-transgenic mosquitoes. Consistent with these observations, cage experiments with mixed populations of transgenic and non-transgenic mosquitoes show that the percentage of transgenic mosquitoes increases when maintained on Plasmodium-infected blood.

Conclusions

Our results suggest that the expression of an anti-Plasmodium effector gene gives transgenic mosquitoes a fitness advantage when fed malaria-infected blood. These findings have important implications for future applications of transgenic mosquito technology in malaria control.

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.

Identification of novel arthropod vector G protein-coupled receptors

BACKGROUND

The control of vector-borne diseases, such as malaria, dengue fever, and typhus fever is often achieved with the use of insecticides. Unfortunately, insecticide resistance is becoming common among different vector species. There are currently no chemical alternatives to these insecticides because new human-safe classes of molecules have yet to be brought to the vector-control market. The identification of novel targets offer opportunities for rational design of new chemistries to control vector populations. One target family, G protein-coupled receptors (GPCRs), has remained relatively under explored in terms of insecticide development.

METHODS

A novel classifier, Ensemble*, for vector GPCRs was developed. Ensemble* was validated and compared to existing classifiers using a set of all known GPCRs from Aedes aegypti, Anopheles gambiae, Apis Mellifera, Drosophila melanogaster, Homo sapiens, and Pediculus humanus. Predictions for unidentified sequences from Ae. aegypti, An. gambiae, and Pe. humanus were validated. Quantitative RT-PCR expression analysis was performed on previously-known and newly discovered Ae. aegypti GPCR genes.

RESULTS

We present a new analysis of GPCRs in the genomes of Ae, aegypti, a vector of dengue fever, An. gambiae, a primary vector of Plasmodium falciparum that causes malaria, and Pe. humanus, a vector of epidemic typhus fever, using a novel GPCR classifier, Ensemble*, designed for insect vector species. We identified 30 additional putative GPCRs, 19 of which we validated. Expression of the newly discovered Ae. aegypti GPCR genes was confirmed via quantitative RT-PCR.

CONCLUSION

A novel GPCR classifier for insect vectors, Ensemble*, was developed and GPCR predictions were validated. Ensemble* and the validation pipeline were applied to the genomes of three insect vectors (Ae. aegypti, An. gambiae, and Pe. humanus), resulting in the identification of 52 GPCRs not previously identified, of which 11 are predicted GPCRs, and 19 are predicted and confirmed GPCRs.

ASGARD: an open-access database of annotated transcriptomes for emerging model arthropod species

The increased throughput and decreased cost of next-generation sequencing (NGS) have shifted the bottleneck genomic research from sequencing to annotation, analysis and accessibility. This is particularly challenging for research communities working on organisms that lack the basic infrastructure of a sequenced genome, or an efficient way to utilize whatever sequence data may be available. Here we present a new database, the Assembled Searchable Giant Arthropod Read Database (ASGARD). This database is a repository and search engine for transcriptomic data from arthropods that are of high interest to multiple research communities but currently lack sequenced genomes. We demonstrate the functionality and utility of ASGARD using de novo assembled transcriptomes from the milkweed bug Oncopeltus fasciatus, the cricket Gryllus bimaculatus and the amphipod crustacean Parhyale hawaiensis. We have annotated these transcriptomes to assign putative orthology, coding region determination, protein domain identification and Gene Ontology (GO) term annotation to all possible assembly products. ASGARD allows users to search all assemblies by orthology annotation, GO term annotation or Basic Local Alignment Search Tool. User-friendly features of ASGARD include search term auto-completion suggestions based on database content, the ability to download assembly product sequences in FASTA format, direct links to NCBI data for predicted orthologs and graphical representation of the location of protein domains and matches to similar sequences from the NCBI non-redundant database. ASGARD will be a useful repository for transcriptome data from future NGS studies on these and other emerging model arthropods, regardless of sequencing platform, assembly or annotation status. This database thus provides easy, one-stop access to multi-species annotated transcriptome information. We anticipate that this database will be useful for members of multiple research communities, including developmental biology, physiology, evolutionary biology, ecology, comparative genomics and phylogenomics.

Database URL:asgard.rc.fas.harvard.edu

RNAi in Arthropods: Insight into the Machinery and Applications for Understanding the Pathogen-Vector Interface

The availability of genome sequencing data in combination with knowledge of expressed genes via transcriptome and proteome data has greatly advanced our understanding of arthropod vectors of disease. Not only have we gained insight into vector biology, but also into their respective vector-pathogen interactions. By combining the strengths of postgenomic databases and reverse genetic approaches such as RNAi, the numbers of available drug and vaccine targets, as well as number of transgenes for subsequent transgenic or paratransgenic approaches, have expanded. These are now paving the way for in-field control strategies of vectors and their pathogens. Basic scientific questions, such as understanding the basic components of the vector RNAi machinery, is vital, as this allows for the transfer of basic RNAi machinery components into RNAi-deficient vectors, thereby expanding the genetic toolbox of these RNAi-deficient vectors and pathogens. In this review, we focus on the current knowledge of arthropod vector RNAi machinery and the impact of RNAi on understanding vector biology and vector-pathogen interactions for which vector genomic data is available on VectorBase.

Gal4-based enhancer-trapping in the malaria mosquito Anopheles stephensi

Transposon-based forward and reverse genetic technologies will contribute greatly to ongoing efforts to study mosquito functional genomics. A piggyBac transposon-based enhancer-trap system was developed that functions efficiently in the human malaria vector, Anopheles stephensi. The system consists of six transgenic lines of Anopheles stephensi, each with a single piggyBac-Gal4 element in a unique genomic location; six lines with a single piggyBac-UAStdTomato element; and two lines, each with a single Minos element containing the piggyBac-transposase gene under the regulatory control of the hsp70 promoter from Drosophila melanogaster. Enhancer detection depended upon the efficient remobilization of piggyBac-Gal4 transposons, which contain the yeast transcription factor gene Gal4 under the regulatory control of a basal promoter. Gal4 expression was detected through the expression of the fluorescent protein gene tdTomato under the regulatory control of a promoter with Gal4-binding UAS elements. From five genetic screens for larval- and adult-specific enhancers, 314 progeny were recovered from 24,250 total progeny (1.3%) with unique patterns of tdTomato expression arising from the influence of an enhancer. The frequency of piggyBac remobilization and enhancer detection was 2.5- to 3-fold higher in female germ lines compared with male germ lines. A small collection of enhancer-trap lines are described in which Gal4 expression occurred in adult female salivary glands, midgut, and fat body, either singly or in combination. These three tissues play critical roles during the infection of Anopheles stephensi by malaria-causing Plasmodium parasites. This system and the lines generated using it will be valuable resources to ongoing mosquito functional genomics efforts.

Novel RNAi-mediated approach to G protein-coupled receptor deorphanization: proof of principle and characterization of a planarian 5-HT receptor

G protein-coupled receptors (GPCRs) represent the largest known superfamily of membrane proteins extending throughout the Metazoa. There exists ample motivation to elucidate the functional properties of GPCRs given their role in signal transduction and their prominence as drug targets. In many target organisms, these efforts are hampered by the unreliable nature of heterologous receptor expression platforms. We validate and describe an alternative loss-of-function approach for ascertaining the ligand and G protein coupling properties of GPCRs in their native cell membrane environment. Our efforts are focused on the phylum Platyhelminthes, given the heavy health burden exacted by pathogenic flatworms, as well as the role of free-living flatworms as model organisms for the study of developmental biology. RNA interference (RNAi) was used in conjunction with a biochemical endpoint assay to monitor cAMP modulation in response to the translational suppression of individual receptors. As proof of principle, this approach was used to confirm the neuropeptide GYIRFamide as the cognate ligand for the planarian neuropeptide receptor GtNPR-1, while revealing its endogenous coupling to Gα_i/o. The method was then extended to deorphanize a novel Gα_s-coupled planarian serotonin receptor, DtSER-1. A bioinformatics protocol guided the selection of receptor candidates mediating 5-HT-evoked responses. These results provide functional data on a neurotransmitter central to flatworm biology, while establishing the great potential of an RNAi-based deorphanization protocol. Future work can help optimize and adapt this protocol for higher-throughput platforms as well as other phyla.

Juvenile hormone (JH) esterase of the mosquito Culex quinquefasciatus is not a target of the JH analog insecticide methoprene

Juvenile hormones (JHs) are essential sesquiterpenes that control insect development and reproduction. JH analog (JHA) insecticides such as methoprene are compounds that mimic the structure and/or biological activity of JH. In this study we obtained a full-length cDNA, cqjhe, from the southern house mosquito Culex quinquefasciatus that encodes CqJHE, an esterase that selectively metabolizes JH. Unlike other recombinant esterases that have been identified from dipteran insects, CqJHE hydrolyzed JH with specificity constant (k_cat/K_M ratio) and V_max values that are common among JH esterases (JHEs). CqJHE showed picomolar sensitivity to OTFP, a JHE-selective inhibitor, but more than 1000-fold lower sensitivity to DFP, a general esterase inhibitor. To our surprise, CqJHE did not metabolize the isopropyl ester of methoprene even when 25 pmol of methoprene was incubated with an amount of CqJHE that was sufficient to hydrolyze 7,200 pmol of JH to JH acid under the same assay conditions. In competition assays in which both JH and methoprene were available to CqJHE, methoprene did not show any inhibitory effects on the JH hydrolysis rate even when methoprene was present in the assay at a 10-fold higher concentration relative to JH. Our findings indicated that JHE is not a molecular target of methoprene. Our findings also do not support the hypothesis that methoprene functions in part by inhibiting the action of JHE.

Mosquito genomics: progress and challenges

The whole-genome sequencing of mosquitoes has facilitated our understanding of fundamental biological processes at their basic molecular levels and holds potential for application to mosquito control and prevention of mosquito-borne disease transmission. Draft genome sequences are available for Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. Collectively, these represent the major vectors of African malaria, dengue fever and yellow fever viruses, and lymphatic filariasis, respectively. Rapid advances in genome technologies have revealed detailed information on genome architecture as well as phenotype-specific transcriptomics and proteomics. These resources allow for detailed comparative analyses within and across populations as well as species. Next-generation sequencing technologies will likely promote a proliferation of genome sequences for additional mosquito species as well as for individual insects. Here we review the current status of genome research in mosquitoes and identify potential areas for further investigations.

BioMart Central Portal: an open database network for the biological community

BioMart Central Portal is a first of its kind, community-driven effort to provide unified access to dozens of biological databases spanning genomics, proteomics, model organisms, cancer data, ontology information and more. Anybody can contribute an independently maintained resource to the Central Portal, allowing it to be exposed to and shared with the research community, and linking it with the other resources in the portal. Users can take advantage of the common interface to quickly utilize different sources without learning a new system for each. The system also simplifies cross-database searches that might otherwise require several complicated steps. Several integrated tools streamline common tasks, such as converting between ID formats and retrieving sequences. The combination of a wide variety of databases, an easy-to-use interface, robust programmatic access and the array of tools make Central Portal a one-stop shop for biological data querying. Here, we describe the structure of Central Portal and show example queries to demonstrate its capabilities.

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.

An automated homology-based approach for identifying transposable elements

Background

Transposable elements (TEs) are mobile sequences found in nearly all eukaryotic genomes. They have the ability to move and replicate within a genome, often influencing genome evolution and gene expression. The identification of TEs is an important part of every genome project. The number of sequenced genomes is rapidly rising, and the need to identify TEs within them is also growing. The ability to do this automatically and effectively in a manner similar to the methods used for genes is of increasing importance. There exist many difficulties in identifying TEs, including their tendency to degrade over time and that many do not adhere to a conserved structure. In this work, we describe a homology-based approach for the automatic identification of high-quality consensus TEs, aimed for use in the analysis of newly sequenced genomes.

Results

We describe a homology-based approach for the automatic identification of TEs in genomes. Our modular approach is dependent on a thorough and high-quality library of representative TEs. The implementation of the approach, named TESeeker, is BLAST-based, but also makes use of the CAP3 assembly program and the ClustalW2 multiple sequence alignment tool, as well as numerous BioPerl scripts. We apply our approach to newly sequenced genomes and successfully identify consensus TEs that are up to 99% identical to manually annotated TEs.

Conclusions

While TEs are known to be a major force in the evolution of genomes, the automatic identification of TEs in genomes is far from mature. In particular, there is a lack of automated homology-based approaches that produce high-quality TEs. Our approach is able to generate high-quality consensus TE sequences automatically, requiring the user to only provide a few basic parameters. This approach is intentionally modular, allowing researchers to use components separately or iteratively. Our approach is most effective for TEs with intact reading frames. The implementation, TESeeker, is available for download as a virtual appliance, while the library of representative TEs is available as a separate download.

Comparative genomics of the odorant-binding and chemosensory protein gene families across the Arthropoda: origin and evolutionary history of the chemosensory system

Chemoreception is a biological process essential for the survival of animals, as it allows the recognition of important volatile cues for the detection of food, egg-laying substrates, mates, or predators, among other purposes. Furthermore, its role in pheromone detection may contribute to evolutionary processes, such as reproductive isolation and speciation. This key role in several vital biological processes makes chemoreception a particularly interesting system for studying the role of natural selection in molecular adaptation. Two major gene families are involved in the perireceptor events of the chemosensory system: the odorant-binding protein (OBP) and chemosensory protein (CSP) families. Here, we have conducted an exhaustive comparative genomic analysis of these gene families in 20 Arthropoda species. We show that the evolution of the OBP and CSP gene families is highly dynamic, with a high number of gains and losses of genes, pseudogenes, and independent origins of subfamilies. Taken together, our data clearly support the birth-and-death model for the evolution of these gene families with an overall high gene turnover rate. Moreover, we show that the genome organization of the two families is significantly more clustered than expected by chance and, more important, that this pattern appears to be actively maintained across the Drosophila phylogeny. Finally, we suggest the homologous nature of the OBP and CSP gene families, dating back their most recent common ancestor after the terrestrialization of Arthropoda (380--450 Ma) and we propose a scenario for the origin and diversification of these families.

Small RNA profiling of Dengue virus-mosquito interactions implicates the PIWI RNA pathway in anti-viral defense

BACKGROUND

Small RNA (sRNA) regulatory pathways (SRRPs) are important to anti-viral defence in mosquitoes. To identify critical features of the virus infection process in Dengue serotype 2 (DENV2)-infected Ae. aegypti, we deep-sequenced small non-coding RNAs. Triplicate biological replicates were used so that rigorous statistical metrics could be applied.

RESULTS

In addition to virus-derived siRNAs (20-23 nts) previously reported for other arbovirus-infected mosquitoes, we show that PIWI pathway sRNAs (piRNAs) (24-30 nts) and unusually small RNAs (usRNAs) (13-19 nts) are produced in DENV-infected mosquitoes. We demonstrate that a major catalytic enzyme of the siRNA pathway, Argonaute 2 (Ago2), co-migrates with a ~1 megadalton complex in adults prior to bloodfeeding. sRNAs were cloned and sequenced from Ago2 immunoprecipitations. Viral sRNA patterns change over the course of infection. Host sRNAs were mapped to the published aedine transcriptome and subjected to analysis using edgeR (Bioconductor). We found that sRNA profiles are altered early in DENV2 infection, and mRNA targets from mitochondrial, transcription/translation, and transport functional categories are affected. Moreover, small non-coding RNAs (ncRNAs), such as tRNAs, spliceosomal U RNAs, and snoRNAs are highly enriched in DENV-infected samples at 2 and 4 dpi.

CONCLUSIONS

These data implicate the PIWI pathway in anti-viral defense. Changes to host sRNA profiles indicate that specific cellular processes are affected during DENV infection, such as mitochondrial function and ncRNA levels. Together, these data provide important progress in understanding the DENV2 infection process in Ae. aegypti.

Hyperexpression of the X chromosome in both sexes results in extensive female bias of X-linked genes in the flour beetle

A genome's ability to produce two separate sexually dimorphic phenotypes is an intriguing biological mystery. Microarray-based studies of a handful of model systems suggest that much of the mystery can be explained by sex-biased gene expression evolved in response to sexually antagonistic selection. We present the first whole-genome study of sex-biased expression in the red flour beetle, Tribolium castaneum. Tribolium is a model for the largest eukaryotic order, Coleoptera, and we show that in whole-body adults, ∼20% of the transcriptome is differentially regulated between the sexes. Among T. castaneum, Drosophila melanogaster, and Anopheles gambiae, we identify 416 1:1:1 orthologs with conserved sex-biased expression. Overrepresented functional categories among sex-biased genes are primarily those involved in gamete production and development. The genomic distribution of sex-biased genes in T. castaneum is distinctly nonrandom, with the strongest deficit of male-biased genes on the X chromosome (9 of 793) of any species studied to date. Tribolium also shows a significant enrichment of X-linked female-biased genes (408 of 793). Our analyses suggest that the extensive female bias of Tribolium X chromosome gene expression is due to hyperexpression of X-linked genes in both males and females. We propose that the overexpression of X chromosomes in females is an evolutionary side effect of the need to dosage compensate in males and that mechanisms to reduce female X chromosome gene expression to autosomal levels are sufficient but imperfect.

Mobility properties of the Hermes transposable element in transgenic lines of Aedes aegypti

The Hermes transposable element has been used to genetically transform a wide range of insect species, including the mosquito, Aedes aegypti, a vector of several important human pathogens. Hermes integrations into the mosquito germline are characterized by the non-canonical integration of the transposon and flanking plasmid and, once integrated, Hermes is stable in the presence of its transposase. In an effort to improve the post-integration mobility of Hermes in the germline of Ae. aegypti, a transgenic helper Mos1 construct expressing Hermes transposase under the control of a testis-specific promoter was crossed to a separate transgenic strain containing a target Hermes transposon. In less than 1% of the approximately 1,500 progeny from jumpstarter lines analyzed, evidence of putative Hermes germline remobilizations were detected. These recovered transposition events occur through an aberrant mechanism and provide insight into the non-canonical cut-and-paste transposition of Hermes in the germ line of Ae. aegypti.

Alternative exon usage in the single CPT1 gene of Drosophila generates functional diversity in the kinetic properties of the enzyme: differential expression of alternatively spliced variants in Drosophila tissues

The Drosophila melanogaster genome contains only one CPT1 gene (Jackson, V. N., Cameron, J. M., Zammit, V. A., and Price, N. T. (1999) Biochem. J. 341, 483-489). We have now extended our original observation to all insect genomes that have been sequenced, suggesting that a single CPT1 gene is a universal feature of insect genomes. We hypothesized that insects may be able to generate kinetically distinct variants by alternative splicing of their single CPT1 gene. Analysis of the insect genomes revealed that (a) the single CPT1 gene in each and every insect genome contains two alternative exons and (ii) in all cases, the putative alternative splicing site occurs within a small region corresponding to 21 amino acid residues that are known to be essential for the binding of substrates and of malonyl-CoA in mammalian CPT1A. We performed PCR analyses of mRNA from different Drosophila tissues; both of the anticipated splice variants of CPT1 mRNA were found to be expressed in all of the tissues tested (both in larvae and adults), with the expression level for one of the splice variants being significantly different between flight muscle and the fat body of adult Drosophila. Heterologous expression of the full-length cDNAs corresponding to the two putative variants of Drosophila CPT1 in the yeast Pichia pastoris revealed two important differences between the properties of the two variants: (i) their affinity (K(0.5)) for one of the substrates, palmitoyl-CoA, differed by 5-fold, and (ii) the sensitivity to inhibition by malonyl-CoA at fixed, higher palmitoyl-CoA concentrations was 2-fold different and associated with different kinetics of inhibition. These data indicate that alternative splicing that specifically affects a structurally crucial region of the protein is an important mechanism through which functional diversity of CPT1 kinetics is generated from the single gene that occurs in insects.

Transcriptome analysis of reproductive tissue and intrauterine developmental stages of the tsetse fly (Glossina morsitans morsitans)

Background

Tsetse flies, vectors of African trypanosomes, undergo viviparous reproduction (the deposition of live offspring). This reproductive strategy results in a large maternal investment and the deposition of a small number of progeny during a female's lifespan. The reproductive biology of tsetse has been studied on a physiological level; however the molecular analysis of tsetse reproduction requires deeper investigation. To build a foundation from which to base molecular studies of tsetse reproduction, a cDNA library was generated from female tsetse (Glossina morsitans morsitans) reproductive tissues and the intrauterine developmental stages. 3438 expressed sequence tags were sequenced and analyzed.

Results

Analysis of a nonredundant catalogue of 1391 contigs resulted in 520 predicted proteins. 475 of these proteins were full length. We predict that 412 of these represent cytoplasmic proteins while 57 are secreted. Comparison of these proteins with other tissue specific tsetse cDNA libraries (salivary gland, fat body/milk gland, and midgut) identified 51 that are unique to the reproductive/immature cDNA library. 11 unique proteins were homologus to uncharacterized putative proteins within the NR database suggesting the identification of novel genes associated with reproductive functions in other insects (hypothetical conserved). The analysis also yielded seven putative proteins without significant homology to sequences present in the public database (unknown genes). These proteins may represent unique functions associated with tsetse's viviparous reproductive cycle. RT-PCR analysis of hypothetical conserved and unknown contigs was performed to determine basic tissue and stage specificity of the expression of these genes.

Conclusion

This paper identifies 51 putative proteins specific to a tsetse reproductive/immature EST library. 11 of these proteins correspond to hypothetical conserved genes and 7 proteins are tsetse specific.

Design and utilization of epitope-based databases and predictive tools

In the last decade, significant progress has been made in expanding the scope and depth of publicly available immunological databases and online analysis resources, which have become an integral part of the repertoire of tools available to the scientific community for basic and applied research. Herein, we present a general overview of different resources and databases currently available. Because of our association with the Immune Epitope Database and Analysis Resource, this resource is reviewed in more detail. Our review includes aspects such as the development of formal ontologies and the type and breadth of analytical tools available to predict epitopes and analyze immune epitope data. A common feature of immunological databases is the requirement to host large amounts of data extracted from disparate sources. Accordingly, we discuss and review processes to curate the immunological literature, as well as examples of how the curated data can be used to generate a meta-analysis of the epitope knowledge currently available for diseases of worldwide concern, such as influenza and malaria. Finally, we review the impact of immunological databases, by analyzing their usage and citations, and by categorizing the type of citations. Taken together, the results highlight the growing impact and utility of immunological databases for the scientific community.

Two frequenins in Drosophila: unveiling the evolutionary history of an unusual neuronal calcium sensor (NCS) duplication

BACKGROUND

Drosophila frequenin (Frq), the homolog of the mammalian neuronal calcium sensor-1 (NCS-1), is a high affinity calcium-binding protein with ubiquitous expression in the nervous system. This protein has an important role in the regulation of neurotransmitter release per synapse, axonal growth and bouton formation. In D. melanogaster, frequenin is encoded by two genes (frq1 and frq2), a very unexpected feature in the Frq/NCS-1 subfamily. These genes are located in tandem in the same genomic region, and their products are 95% identical in their amino acid sequence, clearly indicating their recent origin by gene duplication. Here, we have investigated the factors involved in this unusual feature by examining the molecular evolution of the two frq genes in Drosophila and the evolutionary dynamics of NCS family in a large set of bilaterian species.

RESULTS

Surprisingly, we have found no amino acid replacements fixed across the twelve Drosophila species surveyed. In contrast, synonymous substitutions have been prevalent in the evolution of the coding region of frq1 and frq2, indicating the presence of strong functional constraints following gene duplication. Despite that, we have detected that significant evolutionary rate acceleration had occurred in Frq1 in early times from the duplication, in which positive selection (likely promoting functional diversification) had probably an important role. The analysis of sequence conservation and DNA topology at the non-coding regions of both genes has allowed the identification of DNA regions candidates to be cis-regulatory elements. The results reveal a possible mechanism of regulatory diversification between frq1 and frq2.

CONCLUSIONS

The presence of two frequenins in Drosophila and the rapid accumulation of amino acid substitutions after gene duplication are very unusual features in the evolution of the Frq/NCS-1 subfamily. Here we show that the action of positive selection in concordance with some extent of regulatory diversification might explain these findings. Selected amino acid substitutions in Frq1 likely contributed to the functional divergence between the two duplicates, which, in turn, should have diverged in their regulation by ecdysone-induced early genes.

Next generation transcriptomes for next generation genomes using est2assembly

Background

The decreasing costs of capillary-based Sanger sequencing and next generation technologies, such as 454 pyrosequencing, have prompted an explosion of transcriptome projects in non-model species, where even shallow sequencing of transcriptomes can now be used to examine a range of research questions. This rapid growth in data has outstripped the ability of researchers working on non-model species to analyze and mine transcriptome data efficiently.

Results

Here we present a semi-automated platform 'est2assembly' that processes raw sequence data from Sanger or 454 sequencing into a hybrid de-novo assembly, annotates it and produces GMOD compatible output, including a SeqFeature database suitable for GBrowse. Users are able to parameterize assembler variables, judge assembly quality and determine the optimal assembly for their specific needs. We used est2assembly to process Drosophila and Bicyclus public Sanger EST data and then compared them to published 454 data as well as eight new insect transcriptome collections.

Conclusions

Analysis of such a wide variety of data allows us to understand how these new technologies can assist EST project design. We determine that assembler parameterization is as essential as standardized methods to judge the output of ESTs projects. Further, even shallow sequencing using 454 produces sufficient data to be of wide use to the community. est2assembly is an important tool to assist manual curation for gene models, an important resource in their own right but especially for species which are due to acquire a genome project using Next Generation Sequencing.

Discovery of Plasmodium modulators by genome-wide analysis of circulating hemocytes in Anopheles gambiae

Insect hemocytes mediate important cellular immune responses including phagocytosis and encapsulation and also secrete immune factors such as opsonins, melanization factors, and antimicrobial peptides. However, the molecular composition of these important immune cells has not been elucidated in depth, because of their scarcity in the circulating hemolymph, their adhesion to multiple tissues and the lack of primary culture methods to produce sufficient material for a genome-wide analysis. In this study, we report a genome-wide molecular characterization of circulating hemocytes collected from the hemolymph of adult female Anopheles gambiae mosquitoes--the major mosquito vector of human malaria in subSaharan Africa. Their molecular profile identified 1,485 transcripts with enriched expression in these cells, and many of these genes belong to innate immune gene families. This hemocyte-specific transcriptome is compared to those of Drosophila melanogaster and two other mosquitoes, Aedes aegypti and Armigeres subalbatus. We report the identification of two genes as ubiquitous hemocyte markers and several others as hemocyte subpopulation markers. We assess, via an RNAi screen, the roles in development of Plasmodium berghei of 63 genes expressed in hemocytes and provide a molecular comparison of the transcriptome of these cells during malaria infection.