Doug E. Brackney

C6/36 Aedes albopictus Cells Have a Dysfunctional Antiviral RNA Interference Response

Mosquitoes rely on RNA interference (RNAi) as their primary defense against viral infections. To this end, the combination of RNAi and invertebrate cell culture systems has become an invaluable tool in studying virus-vector interactions. Nevertheless, a recent study failed to detect an active RNAi response to West Nile virus (WNV) infection in C6/36 (Aedes albopictus) cells, a mosquito cell line frequently used to study arthropod-borne viruses (arboviruses). Therefore, we sought to determine if WNV actively evades the hosts RNAi response or if C6/36 cells have a dysfunctional RNAi pathway. C6/36 and Drosophila melanogaster S2 cells were infected with WNV (Flaviviridae), Sindbis virus (SINV, Togaviridae) and La Crosse virus (LACV, Bunyaviridae) and total RNA recovered from cell lysates. Small RNA (sRNA) libraries were constructed and subjected to high-throughput sequencing. In S2 cells, virus-derived small interfering RNAs (viRNAs) from all three viruses were predominantly 21 nt in length, a hallmark of the RNAi pathway. However, in C6/36 cells, viRNAs were primarily 17 nt in length from WNV infected cells and 26–27 nt in length in SINV and LACV infected cells. Furthermore, the origin (positive or negative viral strand) and distribution (position along viral genome) of S2 cell generated viRNA populations was consistent with previously published studies, but the profile of sRNAs isolated from C6/36 cells was altered. In total, these results suggest that C6/36 cells lack a functional antiviral RNAi response. These findings are analogous to the type-I interferon deficiency described in Vero (African green monkey kidney) cells and suggest that C6/36 cells may fail to accurately model mosquito-arbovirus interactions at the molecular level.

RNAi Targeting of West Nile Virus in Mosquito Midguts Promotes Virus Diversification

West Nile virus (WNV) exists in nature as a genetically diverse population of competing genomes. This high genetic diversity and concomitant adaptive plasticity has facilitated the rapid adaptation of WNV to North American transmission cycles and contributed to its explosive spread throughout the New World. WNV is maintained in nature in a transmission cycle between mosquitoes and birds, with intrahost genetic diversity highest in mosquitoes. The mechanistic basis for this increase in genetic diversity in mosquitoes is poorly understood. To determine whether the high mutational diversity of WNV in mosquitoes is driven by RNA interference (RNAi), we characterized the RNAi response to WNV in the midguts of orally exposed Culex pipiens quinquefasciatus using high-throughput, massively parallel sequencing and estimated viral genetic diversity. Our data demonstrate that WNV infection in orally exposed vector mosquitoes induces the RNAi pathway and that regions of the WNV genome that are more intensely targeted by RNAi are more likely to contain point mutations compared to weakly targeted regions. These results suggest that, under natural conditions, positive selection of WNV within mosquitoes is stronger in regions highly targeted by the host RNAi response. Further, they provide a mechanistic basis for the relative importance of mosquitoes in driving WNV diversification.

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.

Highly Sensitive and Specific Detection of Rare Variants in Mixed Viral Populations from Massively Parallel Sequence Data

Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing methods and obtain the high sequence coverage needed to detect rare genetic variants (<1%) within an infected host, and to assay variants without prior knowledge. Critical to interpreting deep sequence data sets is the ability to distinguish biological variants from process errors with high sensitivity and specificity. To address this challenge, we describe V-Phaser, an algorithm able to recognize rare biological variants in mixed populations. V-Phaser uses covariation (i.e. phasing) between observed variants to increase sensitivity and an expectation maximization algorithm that iteratively recalibrates base quality scores to increase specificity. Overall, V-Phaser achieved >97% sensitivity and >97% specificity on control read sets. On data derived from a patient after four years of HIV-1 infection, V-Phaser detected 2,015 variants across the ∼10 kb genome, including 603 rare variants (<1% frequency) detected only using phase information. V-Phaser identified variants at frequencies down to 0.2%, comparable to the detection threshold of allele-specific PCR, a method that requires prior knowledge of the variants. The high sensitivity and specificity of V-Phaser enables identifying and tracking changes in low frequency variants in mixed populations such as RNA viruses.

Expression profiling and comparative analyses of seven midgut serine proteases from the yellow fever mosquito, Aedes aegypti.

Aedes aegypti utilizes blood for energy production, egg maturation and replenishment of maternal reserves. The principle midgut enzymes responsible for bloodmeal digestion are endoproteolytic serine-type proteases within the S1.A subfamily. While there are hundreds of serine protease-like genes in the A. aegypti genome, only five are known to be expressed in the midgut. We describe the cloning, sequencing and expression profiling of seven additional serine proteases and provide a genomic and phylogenetic assessment of these findings. Of the seven genes, four are constitutively expressed and three are transcriptionally induced upon blood feeding. The amount of transcriptional induction is strongly correlated among these genes. Alignments reveal that, in general, the conserved catalytic triad, active site and accessory catalytic residues are maintained in these genes and phylogenetic analysis shows that these genes fall within three distinct clades; trypsins, chymotrypsins and serine collagenases. Interestingly, a previously described trypsin consistently arose with other serine collagenases in phylogenetic analyses. These results suggest that multiple gene duplications have arisen within the S1.A subfamily of midgut serine proteases and/or that A. aegypti has evolved an array of proteases with a broad range of substrate specificities for rapid, efficient digestion of bloodmeals.

West Nile Virus Genetic Diversity is Maintained during Transmission by Culex pipiens quinquefasciatus Mosquitoes

Due to error-prone replication, RNA viruses exist within hosts as a heterogeneous population of non-identical, but related viral variants. These populations may undergo bottlenecks during transmission that stochastically reduce variability leading to fitness declines. Such bottlenecks have been documented for several single-host RNA viruses, but their role in the population biology of obligate two-host viruses such as arthropod-borne viruses (arboviruses) in vivo is unclear, but of central importance in understanding arbovirus persistence and emergence. Therefore, we tracked the composition of West Nile virus (WNV; Flaviviridae, Flavivirus) populations during infection of the vector mosquito, Culex pipiens quinquefasciatus to determine whether WNV populations undergo bottlenecks during transmission by this host. Quantitative, qualitative and phylogenetic analyses of WNV sequences in mosquito midguts, hemolymph and saliva failed to document reductions in genetic diversity during mosquito infection. Further, migration analysis of individual viral variants revealed that while there was some evidence of compartmentalization, anatomical barriers do not impose genetic bottlenecks on WNV populations. Together, these data suggest that the complexity of WNV populations are not significantly diminished during the extrinsic incubation period of mosquitoes.

Evaluation of ivermectin mass drug administration for malaria transmission control across different West African environments

BackgroundMass drug administration (MDA) of ivermectin to humans for control and elimination of filarial parasites can kill biting malaria vectors and lead to Plasmodium transmission reduction. This study examines the degree and duration of mosquitocidal effects resulting from single MDAs conducted in three different West African countries, and the subsequent reductions in parity and Plasmodium sporozoite rates.MethodsIndoor-resting, blood-fed and outdoor host-seeking Anopheles spp. were captured on days surrounding MDAs from 2008–2013 in Senegalese, Liberian and Burkinabé villages. Mortality was assessed on a portion of the indoor collection, and parity status was determined on host-seeking mosquitoes. The effect of MDA was then analysed against the time relative to the MDA, the distributed drugs and environmental variables.ResultsAnopheles gambiae survivorship was reduced by 33.9% for one week following MDA and parity rates were significantly reduced for more than two weeks after the MDAs. Sporozoite rates were significantly reduced by >77% for two weeks following the MDAs in treatment villages despite occurring in the middle of intense transmission seasons. These observed effects were consistent across three different West African transmission dynamics.ConclusionsThese data provide a comprehensive and crucial evidence base for the significant reduction in malaria transmission following single ivermectin MDAs across diverse field sites. Despite the limited duration of transmission reduction, these results support the hypothesis that repeated MDAs with optimal timing could help sustainably control malaria as well as filarial transmission.

Population variation of West Nile virus confers a host-specific fitness benefit in mosquitoes.

West Nile virus is similar to most other RNA viruses in that it exists in nature as a genetically diverse population. However, the role of this genetic diversity within natural transmission cycles and its importance to virus perpetuation remain poorly understood. Therefore, we determined whether highly genetically diverse populations are more fit compared to less genetically diverse WNV populations. Specifically, we generated three WNV populations that varied in their genetic diversity and evaluated their fitness relative to genetically marked control WNV in vivo in Culex quinquefasciatus mosquitoes and chickens. Our results demonstrate that high genetic diversity leads to fitness gains in vector mosquitoes, but not chickens.

West African Anopheles gambiae mosquitoes harbor a taxonomically diverse virome including new insect-specific flaviviruses, mononegaviruses, and totiviruses.

Anopheles gambiae are a major vector of malaria in sub-Saharan Africa. Viruses that naturally infect these mosquitoes may impact their physiology and ability to transmit pathogens. We therefore used metagenomics sequencing to search for viruses in adult Anopheles mosquitoes collected from Liberia, Senegal, and Burkina Faso. We identified a number of virus and virus-like sequences from mosquito midgut contents, including 14 coding-complete genome segments and 26 partial sequences. The coding-complete sequences define new viruses in the order Mononegavirales, and the families Flaviviridae, and Totiviridae. The identification of a flavivirus infecting Anopheles mosquitoes broadens our understanding of the evolution and host range of this virus family. This study increases our understanding of virus diversity in general, begins to define the virome of a medically important vector in its natural setting, and lays groundwork for future studies examining the potential impact of these viruses on anopheles biology and disease transmission.

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.