Inter- and intra-host viral diversity in a large seasonal DENV2 outbreak

Whole Genome Sequencing of DENV-2 isolated from Aedes aegypti mosquitoes in Esmeraldas, Ecuador. Genomic epidemiology of genotype III Southern Asian-American in the country

Ecuador is a tropical country reporting Dengue virus (DENV) outbreaks with areas of hyperendemic viral transmission. Entomo-virological surveillance and monitoring effort conducted in the Northwestern border province of Esmeraldas in April 2022, five pools of female Aedes aegypti mosquitoes from a rural community tested positive for DENV serotype 2 by RT-qPCR. One pool was sequenced by Illumina MiSeq, and it corresponded to genotype III Southern Asian-American. Comparison with other genomes revealed genetic similarity to a human DENV genome sequenced in 2021, also from Esmeraldas. Potential introduction events to the country could have originated from Colombia, considering the vicinity of the collection sites to the neighboring country and high human movement. The inclusion of genomic information complements entomo-virological surveillance, providing valuable insights into genetic variants. This contribution enhances our understanding of Dengue virus (DENV) epidemiology in rural areas and guides evidence-based decisions for surveillance and interventions.

Detection of coinfection with Primate Erythroparvovirus 1 and arboviruses (DENV, CHIKV and ZIKV) in individuals with acute febrile illness in the state of Rio Grande do Norte in 2016

BACKGROUND

Arthropod-borne viruses, known as arboviruses, pose substantial risks to global public health. Dengue (DENV), Chikungunya (CHIKV) and Zika (ZIKV) viruses stand out as significant concerns in Brazil and worldwide. Their overlapping clinical manifestations make accurate diagnosis a challenge, underscoring the need for reliable laboratory support. This study employs a comprehensive molecular diagnostic approach to track viral infections in individuals with acute febrile illness, a period marked by widespread outbreaks of DENV, CHIKV and ZIKV.

METHODS

Between January and August 2016, we received a total of 713 serum samples obtained from individuals with acute febrile illness, previously tested for DENV, CHIKV or ZIKV, with initial negative results, from LACEN-NATAL. Of the total 713 samples, 667 were from females (354 of them pregnant) and 46 from males. Molecular diagnosis was conducted using the Multiplex RT-qPCR technique for simultaneous detection of DENV, CHIKV and ZIKV. Additionally, we performed differential diagnosis by RT-qPCR for other viruses of the Flavivirus, Alphavirus Enterovirus genera and qPCR for Primate Erythroparvovirus 1 (B19V) species, in accordance with Ministry of Health guidelines.

RESULTS

Among the 713 cases, 78.2% tested positive for viral infections, including 48% with CHIKV viremia, 0.6% with DENV and 0.1% with ZIKV. Arboviral coinfections totaled 2.4%, including DENV-CHIKV (1.7%) and CHIKV-ZIKV (0.7%). Moreover, 8% exhibited B19V viremia. Simultaneous infections were identified in 17.5%, encompassing B19V-CHIKV (17.1%), B19V-DENV (0.1%), and B19V-ZIKV (0.3%) Triple infections were observed in 1.3% of cases with B19V-DENV-CHIKV (1%) and B19V-CHIKV-ZIKV (0.3%).

CONCLUSION

Molecular testing demonstrated high efficacy in diagnosing prevalent arboviruses and detecting multiple coinfections. This approach helps to elucidate etiologies for symptomatic cases, especially during arbovirus outbreaks, and aids comprehensive surveillance. Our findings underscore the importance of monitoring co-circulating pathogens, such as B19V, with implications for clinical management, particularly in pregnant individuals. This study enhances our understanding of arbovirus epidemiology and reinforces the critical role of molecular diagnosis in disease surveillance and control.

Intracellular Diversity of WNV within Circulating Avian Peripheral Blood Mononuclear Cells Reveals Host-Dependent Patterns of Polyinfection

Arthropod-borne virus (arbovirus) populations exist as mutant swarms that are maintained between arthropods and vertebrates. West Nile virus (WNV) population dynamics are host-dependent. In American crows, purifying selection is weak and population diversity is high compared to American robins, which have 100- to 1000-fold lower viremia. WNV passed in robins leads to fitness gains, whereas that passed in crows does not. Therefore, we tested the hypothesis that high crow viremia allows for higher genetic diversity within individual avian peripheral blood mononuclear cells (PBMCs), reasoning that this could have produced the previously observed host-specific differences in genetic diversity and fitness. Specifically, we infected cells and birds with a molecularly barcoded WNV and sequenced viral RNA from single cells to quantify the number of WNV barcodes in each. Our results demonstrate that the richness of WNV populations within crows far exceeds that in robins. Similarly, rare WNV variants were maintained by crows more frequently than by robins. Our results suggest that increased viremia in crows relative to robins leads to the maintenance of defective genomes and less prevalent variants, presumably through complementation. Our findings further suggest that weaker purifying selection in highly susceptible crows is attributable to this higher viremia, polyinfections and complementation.

Transgene-induced cell death following dengue-2 virus infection in Aedes aegypti

Dengue viruses (DENVs) are mosquito-borne flaviviruses causing millions of human infections each year and pose a challenge for public health systems worldwide. Aedes aegypti is the principal vector species transmitting DENVs to humans. Controlling Ae. aegypti is difficult due to the abundance of breeding sites and increasing insecticide resistance in the vector populations. Developing new vector control strategies is critical for decreasing the disease burden. One potential approach is genetically replacing Ae. aegypti populations with vector populations highly resistant to DENV transmission. Here, we focus on an alternative strategy for generating dengue 2 virus (DENV-2) resistance in genetically-modified Ae. aegypti in which the mosquitoes express an inactive form of Michelob_x (Mx), an antagonist of the Inhibitor of Apoptosis (IAP), to induce apoptosis in those cells in which actively replicating DENV-2 is present. The inactive form of Mx was flanked by the RRRRSAG cleavage motif, which was recognized by the NS2B/NS3 protease of the infecting DENV-2 thereby releasing and activating Mx which then induced apoptosis. Our transgenic strain exhibited a significantly higher mortality rate than the non-transgenic control when infected with DENV-2. We also transfected a DNA construct containing inactive Mx fused to eGFP into C6/36 mosquito cells and indirectly observed Mx activation on days 3 and 6 post-DENV-2 infections. There were clear signs that the viral NS2B/NS3 protease cleaved the transgene, thereby releasing Mx protein into the cytoplasm, as was confirmed by the detection of eGFP expression in infected cells. The present study represents proof of the concept that virus infection can be used to induce apoptosis in infected mosquito cells.

Intracellular diversity of WNV within circulating avian peripheral blood mononuclear cells reveals host-dependent patterns of polyinfection

Error-prone replication of RNA viruses generates the genetic diversity required for adaptation within rapidly changing environments. Thus, arthropod-borne virus (arbovirus) populations exist in nature as mutant swarms that are maintained between arthropods and vertebrates. Previous studies have demonstrated that West Nile virus (WNV) population dynamics are host dependent: In American crows, which experience extremely high viremia, purifying selection is weak and population diversity is high compared to American robins, which have 100 to 1000-fold lower viremia. WNV passed in robins experiences fitness gains, whereas that passed in crows does not. Therefore, we tested the hypothesis that high crow viremia allows higher genetic diversity within individual avian peripheral-blood mononuclear cells (PBMCs), reasoning that this could have produced the previously observed host-specific differences in genetic diversity and fitness. Specifically, we infected cells and birds with a novel, barcoded version of WNV and sequenced viral RNA from single cells to quantify the number of WNV barcodes that each contained. Our results demonstrate that the richness of WNV populations within crows far exceeds that in robins. Similarly, rare WNV variants were maintained by crows more frequently than by robins. Our results suggest that increased viremia in crows relative to robins leads to maintenance of defective genomes and less prevalent variants, presumably through complementation. Our findings further suggest that weaker purifying selection in highly susceptible crows is attributable to this higher viremia, polyinfections and complementation. These studies further document the role of particular, ecologically relevant hosts in shaping virus population structure.

Author Summary

WNV mutational diversity in vertebrates is species-dependent. In crows, low frequency variants are common, and viral populations are more diverse. In robins, fewer mutations become permanent fixtures of the overall viral population. We infected crows, robins and a chicken cell line with a genetically marked (barcoded) WNV. Higher levels of virus led to multiple unique WNV genomes infecting individual cells, even when a genotype was present at low levels in the input viral stock. Our findings suggest that higher levels of circulating virus in natural hosts allow less fit viruses to survive in RNA virus populations through complementation by more fit viruses. This is significant as it allows less represented and less fit viruses to be maintained at low levels until they potentially emerge when virus environments change. Overall our data reveal new insights on the relationships between host susceptibility to high viremia and virus evolution.

Contribution of phylogenetics to understanding the evolution and epidemiology of dengue virus

Dengue virus (DENV) is one of the most important arboviral pathogens in the tropics and subtropics, and nearly one‐third of the world's population is at risk of infection. The transmission of DENV involves a sylvatic cycle between nonhuman primates (NHP) and Aedes genus mosquitoes, and an endemic cycle between human hosts and predominantly Aedes aegypti. DENV belongs to the genus Flavivirus of the family Flaviviridae and consists of four antigenically distinct serotypes (DENV‐1‐4). Phylogenetic analyses of DENV have revealed its origin, epidemiology, and the drivers that determine its molecular evolution in nature. This review discusses how phylogenetic research has improved our understanding of DENV evolution and how it affects viral ecology and improved our ability to analyze and predict future DENV emergence.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Application of Next-Generation Sequencing to Reveal How Evolutionary Dynamics of Viral Population Shape Dengue Epidemiology

Dengue viral (DENV) infection results in a wide spectrum of clinical manifestations from asymptomatic, mild fever to severe hemorrhage diseases upon infection. Severe dengue is the leading cause of pediatric deaths and/or hospitalizations, which are a major public health burden in dengue-endemic or hyperendemic countries. Like other RNA viruses, DENV continues to evolve. Adaptive mutations are obscured by the major consensus sequence (so-called wild-type sequences) and can only be identified once they become the dominant viruses in the virus population, a process that can take months or years. Traditional surveillance systems still rely on Sanger consensus sequencing. However, with the recent advancement of high-throughput next-generation sequencing (NGS) technologies, the genome-wide investigation of virus population within-host and between-hosts becomes achievable. Thus, viral population sequencing by NGS can increase our understanding of the changing epidemiology and evolution of viral genomics at the molecular level. This review focuses on the studies within the recent decade utilizing NGS in different experimental and epidemiological settings to understand how the adaptive evolution of dengue variants shapes the dengue epidemic and disease severity through its transmission. We propose three types of studies that can be pursued in the future to enhance our surveillance for epidemic prediction and better medical management.

An Outbreak of Human Parvovirus B19 Hidden by Dengue Fever

BACKGROUND

Seasonal outbreaks of dengue often result in hundreds of dengue-suspected cases where a clinical diagnosis cannot be confirmed. Usually, during large outbreaks of dengue and other pathogens that can cause acute febrile illnesses, the search for secondary pathogens with similar disease outcomes is rare.

METHODS

Using total RNA sequencing and targeted diagnostic assays, we discovered an outbreak of parvovirus B19 in dengue-suspected patients that occurred from November 2013 to February 2014.

RESULTS

Of the 182 cases investigated, 63% were viremic for the B19 virus. Moreover, we found that >43% of infected patients had no serological evidence of prior infection. Parvovirus B19 is a typical childhood infection, yet we observed that 82% of the infected patients were adults. Additionally, we perceived that infected adults had significantly higher presentations of myalgia than in children. We also obtained viral protein (VP) 1/VP2 gene nucleotide sequences from 43 patients.

CONCLUSIONS

Our results support the utility of next-generation sequencing for symptomatic patients with unknown etiologies during seasonal outbreaks of dengue and other arborviruses. Our findings could improve the vigilance of hospitals and laboratories by raising awareness of co-circulating pathogens such as parvovirus B19 that may be hidden in plain sight.

The first characterization of a cystatin and a cathepsin L-like peptidase from Aedes aegypti and their possible role in DENV infection by the modulation of apoptosis

Recently, a salivary gland transcriptome study demonstrated that the transcripts of a putative cystatin gene (SeqID AAEL013287; Aacystatins) from Aedes aegypti were increased in DENV2-infected mosquitoes and that silencing of the Aacystatin gene resulted in an increase in DENV titres. In this work, Aacystatin was biochemically characterized; the purified recombinant inhibitor was able to inhibit typical cysteine proteases with a Ki in the nM range. Pulldown assays using Aag2 cell extracts identified a cathepsin L-like peptidase (AaCatL) as a possible target of Aacystatin. Purified recombinant AaCatL had an optimal pH of 5.0 and displayed a preference for Leu, Val and Phe residues at P2, which is common for other cathepsin L-like peptidases. Transcription analysis of Aacystatin and AaCatL in the salivary glands and midgut of DENV2-infected mosquitoes revealed a negative correlation between DENV2 titres and levels of the inhibitor and peptidase, suggesting their involvement in DENV2-mosquito interactions. Considering that apoptosis may play an important role during viral infections, the possible involvement of Aacystatin in staurosporine-induced apoptosis in Aag2 cells was investigated; the results showed higher expression of the inhibitor in treated cells; moreover, pre incubation with rAacystatin was able to increase Aag2 cell viability.

Transgenic Aedes aegypti Mosquitoes Transfer Genes into a Natural Population

In an attempt to control the mosquito-borne diseases yellow fever, dengue, chikungunya, and Zika fevers, a strain of transgenically modified Aedes aegypti mosquitoes containing a dominant lethal gene has been developed by a commercial company, Oxitec Ltd. If lethality is complete, releasing this strain should only reduce population size and not affect the genetics of the target populations. Approximately 450 thousand males of this strain were released each week for 27 months in Jacobina, Bahia, Brazil. We genotyped the release strain and the target Jacobina population before releases began for >21,000 single nucleotide polymorphisms (SNPs). Genetic sampling from the target population six, 12, and 27-30 months after releases commenced provides clear evidence that portions of the transgenic strain genome have been incorporated into the target population. Evidently, rare viable hybrid offspring between the release strain and the Jacobina population are sufficiently robust to be able to reproduce in nature. The release strain was developed using a strain originally from Cuba, then outcrossed to a Mexican population. Thus, Jacobina Ae. aegypti are now a mix of three populations. It is unclear how this may affect disease transmission or affect other efforts to control these dangerous vectors. These results highlight the importance of having in place a genetic monitoring program during such releases to detect un-anticipated outcomes.

Endless Forms: Within-Host Variation in the Structure of the West Nile Virus RNA Genome during Serial Passage in Bird Hosts

The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

RNA viruses are infamous for their high rates of mutation, which produce swarms of genetic variants within individual hosts. To date, analyses of intrahost genetic diversity have focused on the primary genome sequence. However, virus phenotypes are shaped not only by primary sequence but also by the secondary structures into which this sequence folds. Such structures enable viral replication, translation, and binding of small RNAs, yet within-host variation at the structural level has not been adequately explored. We characterized the structural diversity of the 5′ untranslated region (UTR) of populations of West Nile virus (WNV) that had been subject to five serial passages in triplicate in each of three bird species. Viral genomes were sampled from host serum samples at each passage (n = 45 populations) and subjected to next-generation sequencing. For populations derived from passages 1, 3, and 5 (n = 9 populations), we predicted the impact of each mutation occurring at a frequency of ≥1% on the secondary structure of the 5′ UTR. As expected, mutations in double-stranded (DS) regions of the 5′ UTR stem structures caused structural changes of significantly greater magnitude than did mutations in single-stranded (SS) regions. Despite the greater impact of mutations in DS regions, mutations in DS and SS regions occurred at similar frequencies, with no evidence of enhanced selection against mutation in DS regions. In contrast, mutations in two regions that mediate genome cyclization and thereby regulate replication and translation, the 5′ cyclization sequence and the UAR flanking stem (UFS), were suppressed in all three hosts.

IMPORTANCE The enzymes that copy RNA genomes lack proofreading, and viruses that possess RNA genomes, such as West Nile virus, rapidly diversify into swarms of mutant lineages within a host. Intrahost variation of the primary genomic sequence of RNA viruses has been studied extensively because the extent of this variation shapes key virus phenotypes. However, RNA genomes also form complex secondary structures based on within-genome nucleotide complementarity, which are critical regulators of the cyclization of the virus genome that is necessary for efficient replication and translation. We sought to characterize variation in these secondary structures within populations of West Nile virus during serial passage in three bird species. Our study indicates that the intrahost population of West Nile virus is a diverse assortment of RNA secondary structures that should be considered in future analyses of intrahost viral diversity, but some regions that are critical for genome cyclization are conserved within hosts. Besides potential impacts on viral replication, structural diversity can influence the efficacy of small RNA antiviral therapies.

Dengue virus transmission from live donor liver graft

Arboviral transmission through transplanted organs is rare. We report a highly probable case of dengue viral transmission during live donor liver transplantation. Fever with severe thrombocytopenia was observed in the donor and recipient within 6 and 9 days after transplantation, respectively. Dengue diagnosis was confirmed by testing blood and explant tissue from the donor and recipient using dengue-specific NAT (nucleic acid testing) and serology. Serology indicated the donor had secondary dengue infection that ran a mild course. However, the dengue illness in the recipient was severe and deteriorated rapidly, eventually proving fatal. The recipient's explant liver tissue tested negative for viral RNA indicative of a pretransplant naïve status. The prM-Envelope gene sequence analysis of the donor and recipient viral RNA identified a similar serotype (DENV1) with almost 100% sequence identity in the envelope region. Molecular phylogenetic analysis of donor and recipient viral envelope sequences with regional and local dengue strains further confirmed their molecular similarity, suggesting a probable donor-to-recipient transmission via organ transplantation. Screening of living donors for dengue virus may be considered in endemic regions.

Tracking dengue virus type 1 genetic diversity during lineage replacement in an hyperendemic area in Colombia

Dengue virus (DENV) is a flavivirus responsible for the most common and burdensome arthropod-borne viral disease of humans[1]. DENV evolution has been extensively studied on broad geographic and time scales, using sequences from a single gene[2,3]. It is believed that DENV evolution in humans is dominated primarily by purifying selection due to the constraint of maintaining fitness in both humans and mosquitoes[4,5]. Few studies have explored DENV evolutionary dynamics using whole genome sequences, nor have they explored changes in viral diversity that occur during intra-epidemic periods. We used deep sequencing of the viral coding region to characterize DENV-1 evolution in a Colombian population sampled during two high-prevalence dengue seasons in which serotype dominance shifted. Our data demonstrate patterns of strain extinction and replacement within DENV-1 as its prevalence waned and DENV-3 became established. A comparison of whole-genome versus single-gene-based phylogenetic analyses highlights an important difference in evolutionary patterns. We report a trend of higher nonsynonymous to synonymous diversity ratios among non-structural (NS) genes, and statistically significantly higher values among these ratios in the NS1 gene after DENV-1 strain replacement. These results suggest that positive selection could be driving DENV evolution within individual communities. Signals of positive selection coming from distinct samples may be drowned out when combining multiple regions with differing patterns of endemic transmission as commonly done by large-scale geo-temporal assessments. Here, we frame our findings within a small, local transmission history which aids significance. Moreover, these data suggest that the NS1 gene, rather than the E gene, may be a target of positive selection, although not mutually exclusive, and potentially useful sentinel of adaptive changes at the population level.

Application of deep sequencing methods for inferring viral population diversity

The first deep sequencing method was announced in 2005. Due to an increasing number of sequencing data and a reduction in the costs of each sequencing dataset, this innovative technique was soon applied to genetic investigations of viral genome diversity in various viruses, particularly RNA viruses. These deep sequencing findings documented viral epidemiology and evolution and provided high-resolution data on the genetic changes in viral populations. Here, we review deep sequencing platforms that have been applied in viral quasispecies studies. Further, we discuss recent deep sequencing studies on viral inter- and intrahost evolution, drug resistance, and humoral immune selection, especially in emerging and re-emerging viruses. Deep sequencing methods are becoming the standard for providing comprehensive results of viral population diversity, and their applications are discussed.

Impact of Dengue Virus Serotype 2 Strain Diversity on Serological Immune Responses to Dengue

Dengue is a mosquito-borne disease caused by four dengue virus serotypes (DENV1-4) that are loosely categorized by sequence commonalities and antibody recognition profiles. The highly variable envelope protein (E) that is prominently displayed on the surface of DENV is an essential component of vaccines currently under development, yet the impact of using single strains to represent each serotype in tetravalent vaccines has not been adequately studied. We synthesized chimeric E by replacing highly variable residues from a dengue virus serotype 2 vaccine strain (PUO-218) with those from 16 DENV2 lineages spanning 60 years of antigen evolution. Examining sera from human and rhesus macaques challenged with single strains of DENV2, antibody-E interactions were markedly inhibited or enhanced by residues mainly focused within a 480 Ų footprint displayed on the E backbone. The striking impact of E diversity on polyclonal immune responses suggests that frequent antigen updates may be necessary for vaccines to counter shifts in circulating strains.

Inter- and intra-host sequence diversity reveal the emergence of viral variants during an overwintering epidemic caused by dengue virus serotype 2 in southern Taiwan

Purifying selection during dengue viral infection has been suggested as the driving force of viral evolution and the higher complexity of the intra-host quasi-species is thought to offer an adaptive advantage for arboviruses as they cycle between arthropod and vertebrate hosts. However, very few studies have been performed to investigate the viral genetic changes within (intra-host) and between (inter-host) humans in a spatio-temporal scale. Viruses of different serotypes from various countries imported to Taiwan cause annual outbreaks. During 2001-2003, two consecutive outbreaks were caused by dengue virus serotype 2 (DENV-2) and resulted in a larger-scale epidemic with more severe dengue cases in the following year. Phylogenetic analyses showed that the viruses from both events were similar and related to the 2001 DENV-2 isolate from the Philippines. We comprehensively analyzed viral sequences from representative dengue patients and identified three consensus genetic variants, group Ia, Ib and II, with different spatio-temporal population dynamics. The phylodynamic analysis suggested group Ib variants, characterized by lower genetic diversity, transmission rate, and intra-host variant numbers, might play the role of maintenance variants. The residential locations among the patients infected by group Ib variants were in the outer rim of case clusters throughout the 2001-2003 period whereas group Ia and II variants were located in the centers of case clusters, suggesting that group Ib viruses might serve as "sheltered overwintering" variants in an undefined ecological niche. Further deep sequencing of the viral envelope (E) gene directly from individual patient serum samples confirmed the emergence of variants belonging to three quasi-species (group Ia, Ib, and II) and the ancestral role of the viral variants in the latter phase of the 2001 outbreak contributed to the later, larger-scale epidemic beginning in 2002. These findings enhanced our understanding of increasing epidemic severity over time in the same epidemic area. It also highlights the importance of combining phylodynamic and deep sequencing analysis as surveillance tools for detecting dynamic changes in viral variants, particularly searching for and monitoring any specific viral subpopulation. Such subpopulations might have selection advantages in both fitness and transmissibility leading to increased epidemic severity.

Understanding dengue virus evolution to support epidemic surveillance and counter-measure development

Dengue virus (DENV) causes a profound burden of morbidity and mortality, and its global burden is rising due to the co-circulation of four divergent DENV serotypes in the ecological context of globalization, travel, climate change, urbanization, and expansion of the geographic range of the Ae.aegypti and Ae.albopictus vectors. Understanding DENV evolution offers valuable opportunities to enhance surveillance and response to DENV epidemics via advances in RNA virus sequencing, bioinformatics, phylogenetic and other computational biology methods. Here we provide a scoping overview of the evolution and molecular epidemiology of DENV and the range of ways that evolutionary analyses can be applied as a public health tool against this arboviral pathogen.

Serum from dengue virus-infected patients with and without plasma leakage differentially affects endothelial cells barrier function in vitro

Background

Although most of cases of dengue infections are asymptomatic or mild symptomatic some individuals present warning signs progressing to severe dengue in which plasma leakage is a hallmark.

Methodology/Principal findings

The present study used Electric Cell-substrate Impedance Sensing (ECIS^®) which allows for electrical monitoring of cellular barrier function measuring changes in Transendothelial Electric Resistance (TEER) to investigate the parameters associated with dengue induced leakage. Three groups of individuals were tested: dengue-positives with plasma leakage (leakage), dengue-positives without plasma leakage (no leakage), and dengue-negatives (control). Data show that TEER values of human umbilical vein endothelial cells (HUVECs) was significantly lower after incubation with serum from subjects of the leakage group in comparison to the no leakage or control groups. The serum levels of CXCL1, EGF, eotaxin, IFN-γ, sCD40L, and platelets were significantly decreased in the leakage group, while IL-10, IL-6, and IP-10 levels were significantly increased. We also found a strong correlation between TEER values and augmented levels of IP-10, GM-CSF, IL-1α, and IL-8, as well as decreased levels of CXCL1 and platelets.

Conclusions/Significance

The present work shows that the magnitude of the immune response contributes to the adverse plasma leakage outcomes in patients and that serum components are important mediators of changes in endothelial homeostasis during dengue infections. In particular, the increased levels of IP-10 and the decreased levels of CXCL1 and platelets seem to play a significant role in the disruption of vascular endothelium associated with leakage outcomes after DENV infection. These findings may have important implications for both diagnostic and therapeutic approaches to predict and mitigate vascular permeabilization in those experiencing the most severe clinical disease outcomes after dengue infection.

Targeted full-genome amplification and sequencing of dengue virus types 1-4 from South America

We report optimized workflows for full-genome sequencing of dengue viruses (DENVs) 1-4. Based on alignments of publicly available complete genomes we modified and expanded existing primers sets to amplify DENV genotypes that were previously difficult or impossible to sequence. We also report improvements to streamline laboratory handling, including a dual amplification strategy for easy and difficult to sequence "high-copy" and "low-copy" templates, respectively, and normalization of PCR cycling conditions across serotypes. High-copy templates can be sequenced following amplification of as few as 5 overlapping segments covering the complete viral genome, whereas low-copy templates can be sequenced following amplification of no more than 10 overlapping segments of smaller size. These changes have been validated using a balanced set of wild-type DENV genomes (11 of DENV1, 14 of DENV2, 13 of DENV3 and 7 of DENV4) derived from human serum samples collected throughout South America over the past 15 years. The changes described enable generation of complete DENV genomes from wild-type samples without the need for viral enrichment via passaging through laboratory cell lines. This should facilitate quick and cost-effective generation of DENV full-genome sequences of the type needed for accurate epidemiological surveillance and thorough evolutionary studies of wild-type DENVs.

Increasing Clinical Severity during a Dengue Virus Type 3 Cuban Epidemic: Deep Sequencing of Evolving Viral Populations

During the dengue virus type 3 (DENV-3) epidemic that occurred in Havana in 2001 to 2002, severe disease was associated with the infection sequence DENV-1 followed by DENV-3 (DENV-1/DENV-3), while the sequence DENV-2/DENV-3 was associated with mild/asymptomatic infections. To determine the role of the virus in the increasing severity demonstrated during the epidemic, serum samples collected at different time points were studied. A total of 22 full-length sequences were obtained using a deep-sequencing approach. Bayesian phylogenetic analysis of consensus sequences revealed that two DENV-3 lineages were circulating in Havana at that time, both grouped within genotype III. The predominant lineage is closely related to Peruvian and Ecuadorian strains, while the minor lineage is related to Venezuelan strains. According to consensus sequences, relatively few nonsynonymous mutations were observed; only one was fixed during the epidemic at position 4380 in the NS2B gene. Intrahost genetic analysis indicated that a significant minor population was selected and became predominant toward the end of the epidemic. In conclusion, greater variability was detected during the epidemic's progression in terms of significant minority variants, particularly in the nonstructural genes. An increasing trend of genetic diversity toward the end of the epidemic was observed only for synonymous variant allele rates, with higher variability in secondary cases. Remarkably, significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in the structural proteins premembrane (PrM) and envelope (E). Therefore, the dynamic of evolving viral populations in the context of heterotypic antibodies could be related to the increasing clinical severity observed during the epidemic.

IMPORTANCE

Based on the evidence that DENV fitness is context dependent, our research has focused on the study of viral factors associated with intraepidemic increasing severity in a unique epidemiological setting. Here, we investigated the intrahost genetic diversity in acute human samples collected at different time points during the DENV-3 epidemic that occurred in Cuba in 2001 to 2002 using a deep-sequencing approach. We concluded that greater variability in significant minor populations occurred as the epidemic progressed, particularly in the nonstructural genes, with higher variability observed in secondary infection cases. Remarkably, for the first time significant intrahost genetic variation was demonstrated within the same patient during the course of secondary infection with DENV-1/DENV-3, including changes in structural proteins. These findings indicate that high-resolution approaches are needed to unravel molecular mechanisms involved in dengue pathogenesis.

Intra-epidemic evolutionary dynamics of a Dengue virus type 1 population reveal mutant spectra that correlate with disease transmission

Dengue virus (DENV) is currently the most prevalent mosquito-borne viral pathogen. DENVs naturally exist as highly heterogeneous populations. Even though the descriptions on DENV diversity are plentiful, only a few studies have narrated the dynamics of intra-epidemic virus diversity at a fine scale. Such accounts are important to decipher the reciprocal relationship between viral evolutionary dynamics and disease transmission that shape dengue epidemiology. In the current study, we present a micro-scale genetic analysis of a monophyletic lineage of DENV-1 genotype III (epidemic lineage) detected from November 2012 to May 2014. The lineage was involved in an unprecedented dengue epidemic in Singapore during 2013–2014. Our findings showed that the epidemic lineage was an ensemble of mutants (variants) originated from an initial mixed viral population. The composition of mutant spectrum was dynamic and positively correlated with case load. The close interaction between viral evolution and transmission intensity indicated that tracking genetic diversity through time is potentially a useful tool to infer DENV transmission dynamics and thereby, to assess the epidemic risk in a disease control perspective. Moreover, such information is salient to understand the viral basis of clinical outcome and immune response variations that is imperative to effective vaccine design.

Genomic approaches for understanding dengue: insights from the virus, vector, and host

The incidence and geographic range of dengue have increased dramatically in recent decades. Climate change, rapid urbanization and increased global travel have facilitated the spread of both efficient mosquito vectors and the four dengue virus serotypes between population centers. At the same time, significant advances in genomics approaches have provided insights into host–pathogen interactions, immunogenetics, and viral evolution in both humans and mosquitoes. Here, we review these advances and the innovative treatment and control strategies that they are inspiring.

West Nile Virus Population Structure, Injury, and Interferon-Stimulated Gene Expression in the Brain From a Fatal Case of Encephalitis

Background. West Nile virus (WNV) infection in humans can result in severe, acute encephalitis typically involving subcortical gray matter brain regions. West Nile virus replication within specific human brain regions from a human case of acute encephalitis has not been studied.

Methods. We describe a fatal case of WNV encephalitis in which we obtained tissue from specific brain regions at autopsy to evaluate viral-host interactions using next-generation sequencing and immunohistochemistry analysis.

Results. We found that WNV populations in the injured subcortical brain regions exhibited increased amino acid variation and increased expression of specific interferon genes compared with cortical tissues despite similar viral burden.

Conclusions. These observational, patient-based data suggest that neuronal injury and the strength of viral selection pressure may be associated with the level of the innate immune responses. Further studies in human and animal models evaluating the role of innate immune responses on injury patterns and viral selection pressure are needed.

Nature and Extent of Genetic Diversity of Dengue Viruses Determined by 454 Pyrosequencing

Dengue virus (DENV) populations are characteristically highly diverse. Regular lineage extinction and replacement is an important dynamic DENV feature, and most DENV lineage turnover events are associated with increased incidence of disease. The role of genetic diversity in DENV lineage extinctions is not understood. We investigated the nature and extent of genetic diversity in the envelope (E) gene of DENV serotype 1 representing different lineages histories. A region of the DENV genome spanning the E gene was amplified and sequenced by Roche/454 pyrosequencing. The pyrosequencing results identified distinct sub-populations (haplotypes) for each DENV-1 E gene. A phylogenetic tree was constructed with the consensus DENV-1 E gene nucleotide sequences, and the sequences of each constructed haplotype showed that the haplotypes segregated with the Sanger consensus sequence of the population from which they were drawn. Haplotypes determined through pyrosequencing identified a recombinant DENV genome that could not be identified through Sanger sequencing. Nucleotide level sequence diversities of DENV-1 populations determined from SNP analysis were very low, estimated from 0.009–0.01. There were also no stop codon, frameshift or non-frameshift mutations observed in the E genes of any lineage. No significant correlations between the accumulation of deleterious mutations or increasing genetic diversity and lineage extinction were observed (p>0.5). Although our hypothesis that accumulation of deleterious mutations over time led to the extinction and replacement of DENV lineages was ultimately not supported by the data, our data does highlight the significant technical issues that must be resolved in the way in which population diversity is measured for DENV and other viruses. The results provide an insight into the within-population genetic structure and diversity of DENV-1 populations.

Antibodies are not required to a protective immune response against dengue virus elicited in a mouse encephalitis model

Generating neutralizing antibodies have been considered a prerequisite to control dengue virus (DENV) infection. However, T lymphocytes have also been shown to be important in a protective immune state. In order to investigate the contribution of both humoral and cellular immune responses in DENV immunity, we used an experimental model in which a non-lethal DENV2 strain (ACS46) is used to intracranially prime Balb/C mice which develop protective immunity against a lethal DENV2 strain (JHA1). Primed mice generated envelope-specific antibodies and CD8(+) T cell responses targeting mainly non-structural proteins. Immune sera from protected mice did not confer passive protection to naïve mice challenged with the JHA1 strain. In contrast, depletion of CD4(+) and CD8(+) T lymphocytes significantly reduced survival of ACS46-primed mice challenged with the JHA1 strain. Collectively, results presented in this study show that a cellular immune response targeting non-structural proteins are a promising way in vaccine development against dengue.

Complete Genome Sequence of an Atypical Dengue Virus Serotype 2 Lineage Isolated in Brazil

Here, we report the complete polyprotein sequence of a dengue virus 2 strain isolated in Brazil. This virus belongs to the American genotype and has the ability to cause neurovirulence in immunocompetent adult mice. The data presented here may help understand the genetic determinants responsible for neurovirulence.

Error correction and statistical analyses for intra-host comparisons of feline immunodeficiency virus diversity from high-throughput sequencing data

Background

Infection with feline immunodeficiency virus (FIV) causes an immunosuppressive disease whose consequences are less severe if cats are co-infected with an attenuated FIV strain (PLV). We use virus diversity measurements, which reflect replication ability and the virus response to various conditions, to test whether diversity of virulent FIV in lymphoid tissues is altered in the presence of PLV. Our data consisted of the 3′ half of the FIV genome from three tissues of animals infected with FIV alone, or with FIV and PLV, sequenced by 454 technology.

Results

Since rare variants dominate virus populations, we had to carefully distinguish sequence variation from errors due to experimental protocols and sequencing. We considered an exponential-normal convolution model used for background correction of microarray data, and modified it to formulate an error correction approach for minor allele frequencies derived from high-throughput sequencing. Similar to accounting for over-dispersion in counts, this accounts for error-inflated variability in frequencies – and quite effectively reproduces empirically observed distributions. After obtaining error-corrected minor allele frequencies, we applied ANalysis Of VAriance (ANOVA) based on a linear mixed model and found that conserved sites and transition frequencies in FIV genes differ among tissues of dual and single infected cats. Furthermore, analysis of minor allele frequencies at individual FIV genome sites revealed 242 sites significantly affected by infection status (dual vs. single) or infection status by tissue interaction. All together, our results demonstrated a decrease in FIV diversity in bone marrow in the presence of PLV. Importantly, these effects were weakened or undetectable when error correction was performed with other approaches (thresholding of minor allele frequencies; probabilistic clustering of reads). We also queried the data for cytidine deaminase activity on the viral genome, which causes an asymmetric increase in G to A substitutions, but found no evidence for this host defense strategy.

Conclusions

Our error correction approach for minor allele frequencies (more sensitive and computationally efficient than other algorithms) and our statistical treatment of variation (ANOVA) were critical for effective use of high-throughput sequencing data in understanding viral diversity. We found that co-infection with PLV shifts FIV diversity from bone marrow to lymph node and spleen.

Electronic supplementary material

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

Lessons from Reverse Vaccinology for viral vaccine design

Although almost 15 years have passed since the birthdate of Reverse Vaccinology (RV), there are very limited applications of this approach to viral vaccines discovery. Undeniably, RV presents a series of advantages as it can virtually identify all potential antigens coded by a genome, irrespective of their abundance, phase of expression and immunogenicity. Additionally, it can be applied to all pathogens, including those that cannot be grown in vitro. In this review we summarize the few examples of RV application to viruses, in particular the Herpesviridae, and report the advantage and limitations of this approach. Next we focus on the novel approaches and additional technologies to vaccine development including structure based approach (Structural Vaccinology [SV]), synthetic biology and some examples of their application in the development of viral vaccines.

High temporal and spatial diversity in marine RNA viruses implies that they have an important role in mortality and structuring plankton communities

Viruses in the order Picornavirales infect eukaryotes, and are widely distributed in coastal waters. Amplicon deep-sequencing of the RNA dependent RNA polymerase (RdRp) revealed diverse and highly uneven communities of picorna-like viruses in the coastal waters of British Columbia (BC), Canada. Almost 300 000 pyrosequence reads revealed 145 operational taxonomic units (OTUs) based on 95% sequence similarity at the amino-acid level. Each sample had between 24 and 71 OTUs and there was little overlap among samples. Phylogenetic analysis revealed that some clades of OTUs were only found at one site; whereas, other clades included OTUs from all sites. Since most of these OTUs are likely from viruses that infect eukaryotic phytoplankton, and viral isolates infecting phytoplankton are strain-specific; each OTU probably arose from the lysis of a specific phytoplankton taxon. Moreover, the patchiness in OTU distribution, and the high turnover of viruses in the mixed layer, implies continuous infection and lysis by RNA viruses of a diverse array of eukaryotic phytoplankton taxa. Hence, these viruses are likely important elements structuring the phytoplankton community, and play a significant role in nutrient cycling and energy transfer.

Restriction of Francisella novicida genetic diversity during infection of the vector midgut

The genetic diversity of pathogens, and interactions between genotypes, can strongly influence pathogen phenotypes such as transmissibility and virulence. For vector-borne pathogens, both mammalian hosts and arthropod vectors may limit pathogen genotypic diversity (number of unique genotypes circulating in an area) by preventing infection or transmission of particular genotypes. Mammalian hosts often act as “ecological filters” for pathogen diversity, where novel variants are frequently eliminated because of stochastic events or fitness costs. However, whether vectors can serve a similar role in limiting pathogen diversity is less clear. Here we show using Francisella novicida and a natural tick vector of Francisella spp. (Dermacentor andersoni), that the tick vector acted as a stronger ecological filter for pathogen diversity compared to the mammalian host. When both mice and ticks were exposed to mixtures of F. novicida genotypes, significantly fewer genotypes co-colonized ticks compared to mice. In both ticks and mice, increased genotypic diversity negatively affected the recovery of available genotypes. Competition among genotypes contributed to the reduction of diversity during infection of the tick midgut, as genotypes not recovered from tick midguts during mixed genotype infections were recovered from tick midguts during individual genotype infection. Mediated by stochastic and selective forces, pathogen genotype diversity was markedly reduced in the tick. We incorporated our experimental results into a model to demonstrate how vector population dynamics, especially vector-to-host ratio, strongly affected pathogen genotypic diversity in a population over time. Understanding pathogen genotypic population dynamics will aid in identification of the variables that most strongly affect pathogen transmission and disease ecology.

Co-infection, the presence of multiple genotypes of the same pathogen species within an infected individual, is common. Genotype diversity, defined as the number of unique genotypes, and the interaction between genotypes, can strongly influence virulence and pathogen transmission. Understanding how genotypic diversity affects transmission of pathogens that naturally cycle among disparate hosts, such as vector-borne pathogens, is especially important as the capacity of the host and vector to sustain genotypic diversity may differ. To address this, we exposed Dermacentor andersoni ticks, via infected mice, to variably diverse populations of Francisella novicida genotypes. Interestingly, we found that ticks served as greater ecological filters for genotypic diversity compared to mice. This loss in genotypic diversity was due to both stochastic and selective forces. Based on these data and a model, we determined that high numbers of ticks in an environment support high genotypic diversity, while genotypic diversity will be lost rapidly in environments with low tick numbers. Together, these results provide evidence that vector population dynamics, vector-to-host ratios, and competition among pathogen genotypes play critical roles in the maintenance of pathogen genotypic diversity.

Serum angiopoietin-2 and soluble VEGF receptor 2 are surrogate markers for plasma leakage in patients with acute dengue virus infection

BACKGROUND

Endothelial cell dysfunction is believed to play an important role in the pathogenesis of plasma leakage in patients with acute dengue virus (DENV) infection. Several factors, produced by activated endothelial cells, have been associated with plasma leakage or severe disease in patients with infectious diseases.

OBJECTIVES

The aim of this study was to investigate which of these markers could serve as a surrogate marker for the occurrence of plasma leakage in patients with acute DENV infection.

STUDY DESIGN

A case-control study was performed in patients with acute DENV infection in Santos, Brazil. Plasma leakage was detected with X-ray and/or ultrasound examination at admission. Serum levels of soluble endoglin, endothelin-1, angiopoietin-2, VEGF, soluble VEGFR-2, MMP-2, MMP-9, TIMP-1 and TIMP-2 were determined using commercially available ELISAs.

RESULTS

Increased levels of angiopoietin-2, endothelin-1 and MMP-2 and decreased levels of soluble VEGFR-2 were significantly associated with the occurrence of plasma leakage. An unsupervised cluster analysis confirmed that angiopoietin-2 and soluble VEGFR-2 were strongly associated with clinical apparent vascular leakage.

CONCLUSION

Angiopoietin-2 and soluble VEGFR-2 can serve as surrogate markers for the occurrence of plasma leakage in patients with acute DENV infection.

Origin and distribution of divergent dengue virus: novel database construction and phylogenetic analyses

Dengue virus (DENV), a mosquito-borne agent that exists as four serotypes (DENV-1–4), induces dengue illness. DENV has a positive-sense, ssRNA genome of approximately 11 kb that encodes a capsid protein, a premembrane protein and an envelope glycoprotein, in addition to seven nonstructural proteins. These individual genes show sequence variations that can be analyzed phylogenetically to yield several genotypes within each serotype. Here, the sequences of individual DENV genes were collected and used to construct a novel DENV database. This database was then used to characterize the evolution of individual genotypes in several countries. Interestingly, the database provided evidence for recombination between two or three different genotypes to yield new genotypes. This novel database will be available on the internet and is expected to be highly useful for dengue genetic studies, including phylogenetic analyses.