Corey L. Campbell

Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics.

Closing the Vector Circle The genome sequence of Culex quinquefasciatus offers a representative of the third major genus of mosquito disease vectors for comparative analysis. In a major international effort, Arensburger et al. (p. 86) uncovered divergences in the C. quinquefasciatus genome compared with the representatives of the other two genera Aedes aegypti and Anopheles gambiae. The main difference noted is the expansion of numbers of genes, particularly for immunity, oxidoreductive functions, and digestive enzymes, which may reflect specific aspects of the Culex life cycle. Bartholomay et al. (p. 88) explored infection-response genes in Culex in more depth and uncovered 500 immune response-related genes, similar to the numbers seen in Aedes, but fewer than seen in Anopheles or the fruit fly Drosophila melanogaster. The higher numbers of genes were attributed partly to expansions in those encoding serpins, C-type lectins, and fibrinogen-related proteins, consistent with greater immune surveillance and associated signaling needed to monitor the dangers of breeding in polluted, urbanized environments. Transcriptome analysis confirmed that inoculation with unfamiliar bacteria prompted strong immune responses in Culex. The worm and virus pathogens that the mosquitoes transmit naturally provoked little immune activation, however, suggesting that tolerance has evolved to any damage caused by replication of the pathogens in the insects. The genome of a third mosquito species reveals distinctions related to vector capacities and habitat preferences. Culex quinquefasciatus (the southern house mosquito) is an important mosquito vector of viruses such as West Nile virus and St. Louis encephalitis virus, as well as of nematodes that cause lymphatic filariasis. C. quinquefasciatus is one species within the Culex pipiens species complex and can be found throughout tropical and temperate climates of the world. The ability of C. quinquefasciatus to take blood meals from birds, livestock, and humans contributes to its ability to vector pathogens between species. Here, we describe the genomic sequence of C. quinquefasciatus: Its repertoire of 18,883 protein-coding genes is 22% larger than that of Aedes aegypti and 52% larger than that of Anopheles gambiae with multiple gene-family expansions, including olfactory and gustatory receptors, salivary gland genes, and genes associated with xenobiotic detoxification.

Aedes aegypti uses RNA interference in defense against Sindbis virus infection

BackgroundRNA interference (RNAi) is an important anti-viral defense mechanism. The Aedes aegypti genome encodes RNAi component orthologs, however, most populations of this mosquito are readily infected by, and subsequently transmit flaviviruses and alphaviruses. The goal of this study was to use Ae. aegypti as a model system to determine how the mosquitos anti-viral RNAi pathway interacts with recombinant Sindbis virus (SINV; family Togaviridae, genus Alphavirus).ResultsSINV (TR339-eGFP) (+) strand RNA, infectious virus titers and infection rates transiently increased in mosquitoes following dsRNA injection to cognate Ago2, Dcr2, or TSN mRNAs. Detection of SINV RNA-derived small RNAs at 2 and 7 days post-infection in non-silenced mosquitoes provided important confirmation of RNAi pathway activity. Two different recombinant SINV viruses (MRE16-eGFP and TR339-eGFP) with significant differences in infection kinetics were used to delineate vector/virus interactions in the midgut. We show virus-dependent effects on RNAi component transcript and protein levels during infection. Monitoring midgut Ago2, Dcr2, and TSN transcript levels during infection revealed that only TSN transcripts were significantly increased in midguts over blood-fed controls. Ago2 protein levels were depleted immediately following a non-infectious bloodmeal and varied during SINV infection in a virus-dependent manner.ConclusionWe show that silencing RNAi components in Ae. aegypti results in transient increases in SINV replication. Furthermore, Ae. aegypti RNAi is active during SINV infection as indicated by production of virus-specific siRNAs. Lastly, the RNAi response varies in a virus-dependent manner. These data define important features of RNAi anti-viral defense in Ae. aegypti.

Comparison of dengue virus type 2-specific small RNAs from RNA interference-competent and -incompetent mosquito cells.

The exogenous RNA interference (RNAi) pathway is an important antiviral defense against arboviruses in mosquitoes, and virus-specific small interfering (si)RNAs are key components of this pathway. Understanding the biogenesis of siRNAs in mosquitoes could have important ramifications in using RNAi to control arbovirus transmission. Using deep sequencing technology, we characterized dengue virus type 2 (DENV2)-specific small RNAs produced during infection of Aedes aegypti mosquitoes and A. aegypti Aag2 cell cultures and compared them to those produced in the C6/36 Aedes albopictus cell line. We show that the size and mixed polarity of virus-specific small RNAs from DENV-infected A. aegypti cells indicate that they are products of Dicer-2 (Dcr2) cleavage of long dsRNA, whereas C6/36 cells generate DENV2-specific small RNAs that are longer and predominantly positive polarity, suggesting that they originate from a different small RNA pathway. Examination of virus-specific small RNAs after infection of the two mosquito cell lines with the insect-only flavivirus cell fusing agent virus (CFAV) corroborated these findings. An in vitro assay also showed that Aag2 A. aegypti cells are capable of siRNA production, while C6/36 A. albopictus cells exhibit inefficient Dcr2 cleavage of long dsRNA. Defective expression or function of Dcr2, the key initiator of the RNAi pathway, might explain the comparatively robust growth of arthropod-borne viruses in the C6/36 cell line, which has been used frequently as a surrogate for studying molecular interactions between arboviruses and cells of their mosquito hosts.

Comparative genomics of small RNA regulatory pathway components in vector mosquitoes

BackgroundSmall RNA regulatory pathways (SRRPs) control key aspects of development and anti-viral defense in metazoans. Members of the Argonaute family of catalytic enzymes degrade target RNAs in each of these pathways. SRRPs include the microRNA, small interfering RNA (siRNA) and PIWI-type gene silencing pathways. Mosquitoes generate viral siRNAs when infected with RNA arboviruses. However, in some mosquitoes, arboviruses survive antiviral RNA interference (RNAi) and are transmitted via mosquito bite to a subsequent host. Increased knowledge of these pathways and functional components should increase understanding of the limitations of anti-viral defense in vector mosquitoes. To do this, we compared the genomic structure of SRRP components across three mosquito species and three major small RNA pathways.ResultsThe Ae. aegypti, An. gambiae and Cx. pipiens genomes encode putative orthologs for all major components of the miRNA, siRNA, and piRNA pathways. Ae. aegypti and Cx. pipiens have undergone expansion of Argonaute and PIWI subfamily genes. Phylogenetic analyses were performed for these protein families. In addition, sequence pattern recognition algorithms MEME, MDScan and Weeder were used to identify upstream regulatory motifs for all SRRP components. Statistical analyses confirmed enrichment of species-specific and pathway-specific cis-elements over the rest of the genome.ConclusionAnalysis of Argonaute and PIWI subfamily genes suggests that the small regulatory RNA pathways of the major arbovirus vectors, Ae. aegypti and Cx. pipiens, are evolving faster than those of the malaria vector An. gambiae and D. melanogaster. Further, protein and genomic features suggest functional differences between subclasses of PIWI proteins and provide a basis for future analyses. Common UCR elements among SRRP components indicate that 1) key components from the miRNA, siRNA, and piRNA pathways contain NF-kappaB-related and Broad complex transcription factor binding sites, 2) purifying selection has occurred to maintain common pathway-specific elements across mosquito species and 3) species-specific differences in upstream elements suggest that there may be differences in regulatory control among mosquito species. Implications for arbovirus vector competence in mosquitoes are discussed.

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

BackgroundSmall 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.ResultsIn 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.ConclusionsThese 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.

Pathogenomics of Culex quinquefasciatus and Meta-Analysis of Infection Responses to Diverse Pathogens

Closing the Vector Circle The genome sequence of Culex quinquefasciatus offers a representative of the third major genus of mosquito disease vectors for comparative analysis. In a major international effort, Arensburger et al. (p. 86) uncovered divergences in the C. quinquefasciatus genome compared with the representatives of the other two genera Aedes aegypti and Anopheles gambiae. The main difference noted is the expansion of numbers of genes, particularly for immunity, oxidoreductive functions, and digestive enzymes, which may reflect specific aspects of the Culex life cycle. Bartholomay et al. (p. 88) explored infection-response genes in Culex in more depth and uncovered 500 immune response-related genes, similar to the numbers seen in Aedes, but fewer than seen in Anopheles or the fruit fly Drosophila melanogaster. The higher numbers of genes were attributed partly to expansions in those encoding serpins, C-type lectins, and fibrinogen-related proteins, consistent with greater immune surveillance and associated signaling needed to monitor the dangers of breeding in polluted, urbanized environments. Transcriptome analysis confirmed that inoculation with unfamiliar bacteria prompted strong immune responses in Culex. The worm and virus pathogens that the mosquitoes transmit naturally provoked little immune activation, however, suggesting that tolerance has evolved to any damage caused by replication of the pathogens in the insects. The genome of a third mosquito species reveals distinctions related to vector capacities and habitat preferences. The mosquito Culex quinquefasciatus poses a substantial threat to human and veterinary health as a primary vector of West Nile virus (WNV), the filarial worm Wuchereria bancrofti, and an avian malaria parasite. Comparative phylogenomics revealed an expanded canonical C. quinquefasciatus immune gene repertoire compared with those of Aedes aegypti and Anopheles gambiae. Transcriptomic analysis of C. quinquefasciatus genes responsive to WNV, W. bancrofti, and non-native bacteria facilitated an unprecedented meta-analysis of 25 vector-pathogen interactions involving arboviruses, filarial worms, bacteria, and malaria parasites, revealing common and distinct responses to these pathogen types in three mosquito genera. Our findings provide support for the hypothesis that mosquito-borne pathogens have evolved to evade innate immune responses in three vector mosquito species of major medical importance.

Complex Modulation of the Aedes aegypti Transcriptome in Response to Dengue Virus Infection

Dengue fever is the most important arboviral disease world-wide, with Aedes aegypti being the major vector. Interactions between the mosquito host and dengue viruses (DENV) are complex and vector competence varies among geographically-distinct Ae. aegypti populations. Additionally, dengue is caused by four antigenically-distinct viral serotypes (DENV1–4), each with multiple genotypes. Each virus genotype interacts differently with vertebrate and invertebrate hosts. Analyses of alterations in mosquito transcriptional profiles during DENV infection are expected to provide the basis for identifying networks of genes involved in responses to viruses and contribute to the molecular-genetic understanding of vector competence. In addition, this knowledge is anticipated to support the development of novel disease-control strategies. RNA-seq technology was used to assess genome-wide changes in transcript abundance at 1, 4 and 14 days following DENV2 infection in carcasses, midguts and salivary glands of the Ae. aegypti Chetumal strain. DENV2 affected the expression of 397 Ae. aegypti genes, most of which were down-regulated by viral infection. Differential accumulation of transcripts was mainly tissue- and time-specific. Comparisons of our data with other published reports reveal conservation of functional classes, but limited concordance of specific mosquito genes responsive to DENV2 infection. These results indicate the necessity of additional studies of mosquito-DENV interactions, specifically those focused on recently-derived mosquito strains with multiple dengue virus serotypes and genotypes.

MicroRNA levels are modulated in Aedes aegypti after exposure to Dengue-2

To define microRNA (miRNA) involvement during arbovirus infection of Aedes aegypti, we mined deep sequencing libraries of Dengue type 2 (DENV2)‐exposed mosquitoes. Three biological replicates for each timepoint [2, 4 and 9 days post‐exposure (dpe)] and treatment group allowed us to remove the outliers associated with sample‐to‐sample variability. Using edgeR (R Bioconductor), designed for use with replicate deep sequencing data, we determined the log fold‐change (logFC) of miRNA levels [18–23 nucleotides (nt)]. The number of significantly modulated miRNAs increased from ≤5 at 2 and 4 dpe to 23 unique miRNAs by 9 dpe. Putative miRNA targets were predicted by aligning miRNAs to the transcriptome, and the list was reduced to include the intersection of hits found using the Miranda, PITA, and TargetScan algorithms. To further reduce false‐positives, putative targets were validated by cross‐checking them with mRNAs reported in recent DENV2 host response transcriptome reports; 4076 targets were identified. Of these, 464 gene targets have predicted miRNA‐binding sites in 3′ untranslated regions. Context‐specific target functional groups include proteins involved in transport, transcriptional regulation, mitochondrial function, chromatin modification and signal transduction processes known to be required for viral replication and dissemination. The miRNA response is placed in context with other vector host response studies by comparing the predicted targets with those of transcriptome studies. Together, these data are consistent with the hypothesis that profound and persistent changes to gene expression occur in DENV2‐exposed mosquitoes.

Strain variation in the transcriptome of the dengue fever vector, aedes aegypti

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

Experimental Evolution of an RNA Virus in Wild Birds: Evidence for Host-Dependent Impacts on Population Structure and Competitive Fitness

Within hosts, RNA viruses form populations that are genetically and phenotypically complex. Heterogeneity in RNA virus genomes arises due to error-prone replication and is reduced by stochastic and selective mechanisms that are incompletely understood. Defining how natural selection shapes RNA virus populations is critical because it can inform treatment paradigms and enhance control efforts. We allowed West Nile virus (WNV) to replicate in wild-caught American crows, house sparrows and American robins to assess how natural selection shapes RNA virus populations in ecologically relevant hosts that differ in susceptibility to virus-induced mortality. After five sequential passages in each bird species, we examined the phenotype and population diversity of WNV through fitness competition assays and next generation sequencing. We demonstrate that fitness gains occur in a species-specific manner, with the greatest replicative fitness gains in robin-passaged WNV and the least in WNV passaged in crows. Sequencing data revealed that intrahost WNV populations were strongly influenced by purifying selection and the overall complexity of the viral populations was similar among passaged hosts. However, the selective pressures that control WNV populations seem to be bird species-dependent. Specifically, crow-passaged WNV populations contained the most unique mutations (~1.7× more than sparrows, ~3.4× more than robins) and defective genomes (~1.4× greater than sparrows, ~2.7× greater than robins), but the lowest average mutation frequency (about equal to sparrows, ~2.6× lower than robins). Therefore, our data suggest that WNV replication in the most disease-susceptible bird species is positively associated with virus mutational tolerance, likely via complementation, and negatively associated with the strength of selection. These differences in genetic composition most likely have distinct phenotypic consequences for the virus populations. Taken together, these results reveal important insights into how different hosts may contribute to the emergence of RNA viruses.