Viral Landscape Research Articles

On the importance of negative controls in viral landscape phylogeography.

Phylogeographic reconstructions are becoming an established procedure to evaluate the factors that could impact virus spread. While a discrete phylogeographic approach can be used to test predictors of transition rates among discrete locations, alternative continuous phylogeographic reconstructions can also be exploited to investigate the impact of underlying environmental layers on the dispersal velocity of a virus. The two approaches are complementary tools for studying pathogens' spread, but in both cases, care must be taken to avoid misinterpretations. Here, we analyse rice yellow mottle virus (RYMV) sequence data from West and East Africa to illustrate how both approaches can be used to study the impact of environmental factors on the virus’ dispersal frequency and velocity. While it was previously reported that host connectivity was a major determinant of RYMV spread, we show that this was a false positive result due to the lack of appropriate negative controls. We also discuss and compare the phylodynamic tools currently available for investigating the impact of environmental factors on virus spread.

Within-patient mutation frequencies reveal fitness costs of CpG dinucleotides and drastic amino acid changes in HIV.

HIV has a high mutation rate, which contributes to its ability to evolve quickly. However, we know little about the fitness costs of individual HIV mutations in vivo, their distribution and the different factors shaping the viral fitness landscape. We calculated the mean frequency of transition mutations at 870 sites of the pol gene in 160 patients, allowing us to determine the cost of these mutations. As expected, we found high costs for non-synonymous and nonsense mutations as compared to synonymous mutations. In addition, we found that non-synonymous mutations that lead to drastic amino acid changes are twice as costly as those that do not and mutations that create new CpG dinucleotides are also twice as costly as those that do not. We also found that G→A and C→T mutations are more costly than A→G mutations. We anticipate that our new in vivo frequency-based approach will provide insights into the fitness landscape and evolvability of not only HIV, but a variety of microbes.

Metagenomic analysis of Varroa-free Australian honey bees (Apis mellifera) shows a diverse Picornavirales virome.

The viral landscape of the honey bee (Apismellifera) has changed as a consequence of the global spread of the parasitic mite Varroa destructor and accompanying virulent strains of the iflavirus deformed wing virus (DWV), which the mite vectors. The presence of DWV in honey bee populations is known to influence the occurrence of other viruses, suggesting that the current known virome of A. mellifera may be undercharacterized. Here we tested this hypothesis by examining the honey bee virome in Australia, which is uniquely free of parasitic mites or DWV. Using a high-throughput sequencing (HTS) approach, we examined the RNA virome from nine pools of A. mellifera across Australia. In addition to previously reported honey bee viruses, several other insect viruses were detected, including strains related to aphid lethal paralysis virus (ALPV) and Rhopalosiphum padi virus (RhPV), which have recently been identified as infecting honey bees in the USA, as well as several other viruses recently found in Drosophila spp. A further 42 putative novel insect virus genomes spanning the order Picornavirales were assembled, which significantly increases the known viral diversity in A. mellifera. Among these novel genomes, we identified several that were similar (but different) to key A. mellifera viruses, such as DWV, that warrant further investigation. We propose that A. mellifera may be preferentially infected with viruses of the order Picornavirales and that a diverse population of these viruses may be representative of a Varroa-free landscape.

VIII. Chronic liver diseases after eradiation of viral hepatitis-Clinical landscape of NASH and alcoholic liver disease-

VIII. Chronic liver diseases after eradiation of viral hepatitis-Clinical landscape of NASH and alcoholic liver disease-

Evidence of Varroa-mediated deformed wing virus spillover in Hawaii.

Varroa destructor, a parasitic mite of honey bees, is also a vector for viral diseases. The mite displays high host specificity and requires access to colonies of Apis spp. to complete its lifecycle. In contrast, the Deformed Wing Virus (DWV), one of the many viruses transmitted by V. destructor, appears to have a much broader host range. Previous studies have detected DWV in a variety of insect groups that are not directly parasitized by the mite. In this study, we take advantage of the discrete distribution of the Varroa mite in the Hawaiian archipelago to compare DWV prevalence on non-Apis flower visitors, and test whether Varroa presence is linked to a "viral spillover". We selected two islands with different viral landscapes: Oahu, where V. destructor has been present since 2007, and Maui, where the mite is absent. We sampled individuals of Apis mellifera, Ceratina smaragdula, Polistes aurifer, and Polistes exclamens, to assess and compare the DWV prevalence in the Hymenoptera community of the two islands. The results indicated that, as expected, honey bee colonies on Oahu have much higher incidence of DWV compared to Maui. Correspondingly, DWV was detected on the Non-Apis Hymenoptera collected from Oahu, but was absent in the species examined on Maui. The study sites selected shared a similar geography, climate, and insect fauna, but differed in the presence of the Varroa mite, suggesting an indirect, but significant, increase on DWV prevalence in the Hymenoptera community on mite-infected islands.

Absence of deformed wing virus and Varroa destructor in Australia provides unique perspectives on honeybee viral landscapes and colony losses

Honeybee (Apis mellifera) health is threatened globally by the complex interaction of multiple stressors, including the parasitic mite Varroa destructor and a number of pathogenic viruses. Australia provides a unique opportunity to study this pathogenic viral landscape in the absence of V. destructor. We analysed 1,240A. mellifera colonies across Australia by reverse transcription-polymerase chain reaction (RT-PCR) and next-generation sequencing (NGS). Five viruses were prevalent: black queen cell virus (BQCV), sacbrood virus (SBV), Israeli acute paralysis virus (IAPV) and the Lake Sinai viruses (LSV1 and LSV2), of which the latter three were detected for the first time in Australia. We also showed several viruses were absent in our sampling, including deformed wing virus (DWV) and slow bee paralysis virus (SBPV). Our findings highlight that viruses can be highly prevalent in A. mellifera populations independently of V. destructor. Placing these results in an international context, our results support the hypothesis that the co-pathogenic interaction of V. destructor and DWV is a key driver of increased colony losses, but additional stressors such as pesticides, poor nutrition, etc. may enable more severe and frequent colony losses to occur.

A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations.

ABSTRACTUnderstanding the diversity and consequences of viruses present in honey bees is critical for maintaining pollinator health and managing the spread of disease. The viral landscape of honey bees (Apis mellifera) has changed dramatically since the emergence of the parasitic mite Varroa destructor, which increased the spread of virulent variants of viruses such as deformed wing virus. Previous genomic studies have focused on colonies suffering from infections by Varroa and virulent viruses, which could mask other viral species present in honey bees, resulting in a distorted view of viral diversity. To capture the viral diversity within colonies that are exposed to mites but do not suffer the ultimate consequences of the infestation, we examined populations of honey bees that have evolved naturally or have been selected for resistance to Varroa. This analysis revealed seven novel viruses isolated from honey bees sampled globally, including the first identification of negative-sense RNA viruses in honey bees. Notably, two rhabdoviruses were present in three geographically diverse locations and were also present in Varroa mites parasitizing the bees. To characterize the antiviral response, we performed deep sequencing of small RNA populations in honey bees and mites. This provided evidence of a Dicer-mediated immune response in honey bees, while the viral small RNA profile in Varroa mites was novel and distinct from the response observed in bees. Overall, we show that viral diversity in honey bee colonies is greater than previously thought, which encourages additional studies of the bee virome on a global scale and which may ultimately improve disease management.IMPORTANCE Honey bee populations have become increasingly susceptible to colony losses due to pathogenic viruses spread by parasitic Varroa mites. To date, 24 viruses have been described in honey bees, with most belonging to the order Picornavirales. Collapsing Varroa-infected colonies are often overwhelmed with high levels of picornaviruses. To examine the underlying viral diversity in honey bees, we employed viral metatranscriptomics analyses on three geographically diverse Varroa-resistant populations from Europe, Africa, and the Pacific. We describe seven novel viruses from a range of diverse viral families, including two viruses that are present in all three locations. In honey bees, small RNA sequences indicate that these viruses are processed by Dicer and the RNA interference pathway, whereas Varroa mites produce strikingly novel small RNA patterns. This work increases the number and diversity of known honey bee viruses and will ultimately contribute to improved disease management in our most important agricultural pollinator.

KSHV encoded ORF59 modulates histone arginine methylation of the viral genome to promote viral reactivation.

Kaposi’s sarcoma associated herpesvirus (KSHV) persists in a highly-ordered chromatin structure inside latently infected cells with the majority of the viral genome having repressive marks. However, upon reactivation the viral chromatin landscape changes into ‘open’ chromatin through the involvement of lysine demethylases and methyltransferases. Besides methylation of lysine residues of histone H3, arginine methylation of histone H4 plays an important role in controlling the compactness of the chromatin. Symmetric methylation of histone H4 at arginine 3 (H4R3me2s) negatively affects the methylation of histone H3 at lysine 4 (H3K4me3), an active epigenetic mark deposited on the viral chromatin during reactivation. We identified a novel binding partner to KSHV viral DNA processivity factor, ORF59-a protein arginine methyl transferase 5 (PRMT5). PRMT5 is an arginine methyltransferase that dimethylates arginine 3 (R3) of histone H4 in a symmetric manner, one hallmark of condensed chromatin. Our ChIP-seq data of symmetrically methylated H4 arginine 3 showed a significant decrease in H4R3me2s on the viral genome of reactivated cells as compared to the latent cells. Reduction in arginine methylation correlated with the binding of ORF59 on the viral chromatin and disruption of PRMT5 from its adapter protein, COPR5 (cooperator of PRMT5). Binding of PRMT5 through COPR5 is important for symmetric methylation of H4R3 and the expression of ORF59 competitively reduces the association of PRMT5 with COPR5, leading to a reduction in PRMT5 mediated arginine methylation. This ultimately resulted in a reduced level of symmetrically methylated H4R3 and increased levels of H3K4me3 marks, contributing to the formation of an open chromatin for transcription and DNA replication. Depletion of PRMT5 levels led to a decrease in symmetric methylation and increase in viral gene transcription confirming the role of PRMT5 in viral reactivation. In conclusion, ORF59 modulates histone-modifying enzymes to alter the chromatin structure during lytic reactivation.

RNA-binding protein CPEB1 remodels host and viral RNA landscapes.

Host and virus interactions at the post-transcriptional level are critical for infection but remain poorly understood. Human cytomegalovirus (HCMV) is a prevalent herpesvirus family member that causes severe complications in immunocompromised patients and newborns. Here, we perform comprehensive transcriptome-wide analyses revealing that HCMV infection results in widespread alternative splicing (AS), shorter 3′-untranslated regions (3′UTRs) and polyA tail lengthening in host genes. The host RNA binding protein cytoplasmic polyadenylation element binding protein 1 (CPEB1) is highly induced upon infection and ectopic expression of CPEB1 in non-infected cells recapitulates infection-related post-transcriptional changes. CPEB1 is also required for polyA-tail lengthening of viral RNAs important for productive infection. Strikingly, depletion of CPEB1 reverses infection-related cytopathology and post-transcriptional changes, and decreases productive HCMV titers. Host RNA processing is also altered in herpes simplex virus-2 (HSV-2) infected cells, indicating a common theme among herpesvirus infections. Our work is a starting point for therapeutic targeting of host RNA binding proteins in herpesvirus infections.

Hepatitis Delta Virus Detected in Salivary Glands of Sjögren's Syndrome Patients and Recapitulates a Sjögren's Syndrome-Like Phenotype in Vivo.

Background:Low-level, chronic viral infections have been suspect in the development of select autoimmune diseases, including primary Sjögren's syndrome (pSS). Multiple studies have shown stimulation of antiviral response pathways in pSS tissues suggestive of a viral infection. Yet, with this data in hand, a causal link between a viral infection and development of pSS had not been identified. Therefore, a study was designed to further define the viral landscape within pSS-affected salivary gland tissue to identify potential viral-mediated triggers in the pathogenesis of this autoimmune disease.Methods:A viral microarray was utilized to measure viral transcripts present in salivary gland tissue from patients diagnosed with pSS compared to healthy controls. Murine models of salivary gland localized HDV antigen expression were developed to evaluate the capacity of a chronic HDV signature to trigger the development of a pSS-like phenotype.Results:Through this analysis, two distinct viral profiles were identified, including the increased presence of hepatitis delta virus (HDV) in 50% of pSS patients evaluated. Presence of HDV antigen and sequence were confirmed in minor salivary gland tissue. Patients with elevated HDV levels in salivary gland tissue were negative for detectible hepatitis B virus (HBV) surface antigen and antibodies to HBV or HDV. Expression of HDV antigens in vivo resulted in reduced stimulated saliva flow, increase in focal lymphocytic infiltrates, and development of autoantibodies.Conclusion:Identification of HDV in pSS patients and induction of a complete pSS-like phenotype in vivo provides further support of a viral-mediated etiopathology in the development of pSS.

Distinct Viral and Mutational Spectrum of Endemic Burkitt Lymphoma

Endemic Burkitt lymphoma (eBL) is primarily found in children in equatorial regions and represents the first historical example of a virus-associated human malignancy. Although Epstein-Barr virus (EBV) infection and MYC translocations are hallmarks of the disease, it is unclear whether other factors may contribute to its development. We performed RNA-Seq on 20 eBL cases from Uganda and showed that the mutational and viral landscape of eBL is more complex than previously reported. First, we found the presence of other herpesviridae family members in 8 cases (40%), in particular human herpesvirus 5 and human herpesvirus 8 and confirmed their presence by immunohistochemistry in the adjacent non-neoplastic tissue. Second, we identified a distinct latency program in EBV involving lytic genes in association with TCF3 activity. Third, by comparing the eBL mutational landscape with published data on sporadic Burkitt lymphoma (sBL), we detected lower frequencies of mutations in MYC, ID3, TCF3 and TP53, and a higher frequency of mutation in ARID1A in eBL samples. Recurrent mutations in two genes not previously associated with eBL were identified in 20% of tumors: RHOA and cyclin F (CCNF). We also observed that polyviral samples showed lower numbers of somatic mutations in common altered genes in comparison to sBL specimens, suggesting dual mechanisms of transformation, mutation versus virus driven in sBL and eBL respectively.

Modularity enhances the rate of evolution in a rugged fitness landscape

Biological systems are modular, and this modularity affects the evolution of biological systems over time and in different environments. We here develop a theory for the dynamics of evolution in a rugged, modular fitness landscape. We show analytically how horizontal gene transfer couples to the modularity in the system and leads to more rapid rates of evolution at short times. The model, in general, analytically demonstrates a selective pressure for the prevalence of modularity in biology. We use this model to show how the evolution of the influenza virus is affected by the modularity of the proteins that are recognized by the human immune system. Approximately 25% of the observed rate of fitness increase of the virus could be ascribed to a modular viral landscape.

On the front line: quantitative virus dynamics in honeybee (Apis mellifera L.) colonies along a new expansion front of the parasite Varroa destructor.

Over the past fifty years, annual honeybee (Apis mellifera) colony losses have been steadily increasing worldwide. These losses have occurred in parallel with the global spread of the honeybee parasite Varroa destructor. Indeed, Varroa mite infestations are considered to be a key explanatory factor for the widespread increase in annual honeybee colony mortality. The host-parasite relationship between honeybees and Varroa is complicated by the mite's close association with a range of honeybee viral pathogens. The 10-year history of the expanding front of Varroa infestation in New Zealand offered a rare opportunity to assess the dynamic quantitative and qualitative changes in honeybee viral landscapes in response to the arrival, spread and level of Varroa infestation. We studied the impact of de novo infestation of bee colonies by Varroa on the prevalence and titres of seven well-characterised honeybee viruses in both bees and mites, using a large-scale molecular ecology approach. We also examined the effect of the number of years since Varroa arrival on honeybee and mite viral titres. The dynamic shifts in the viral titres of black queen cell virus and Kashmir bee virus mirrored the patterns of change in Varroa infestation rates along the Varroa expansion front. The deformed wing virus (DWV) titres in bees continued to increase with Varroa infestation history, despite dropping infestation rates, which could be linked to increasing DWV titres in the mites. This suggests that the DWV titres in mites, perhaps boosted by virus replication, may be a major factor in maintaining the DWV epidemic after initial establishment. Both positive and negative associations were identified for several pairs of viruses, in response to the arrival of Varroa. These findings provide important new insights into the role of the parasitic mite Varroa destructor in influencing the viral landscape that affects honeybee colonies.

Locking and blocking the viral landscape of an alphavirus with neutralizing antibodies.

Alphaviruses are serious, sometimes lethal human pathogens that belong to the family Togaviridae. The structures of human Venezuelan equine encephalitis virus (VEEV), an alphavirus, in complex with two strongly neutralizing antibody Fab fragments (F5 and 3B4C-4) have been determined using a combination of cryo-electron microscopy and homology modeling. We characterize these monoclonal antibody Fab fragments, which are known to abrogate VEEV infectivity by binding to the E2 (envelope) surface glycoprotein. Both of these antibody Fab fragments cross-link the surface E2 glycoproteins and therefore probably inhibit infectivity by blocking the conformational changes that are required for making the virus fusogenic. The F5 Fab fragment cross-links E2 proteins within one trimeric spike, whereas the 3B4C-4 Fab fragment cross-links E2 proteins from neighboring spikes. Furthermore, F5 probably blocks the receptor-binding site, whereas 3B4C-4 sterically hinders the exposure of the fusion loop at the end of the E2 B-domain. Alphaviral infections are transmitted mainly by mosquitoes. Venezuelan equine encephalitis virus (VEEV) is an alphavirus with a wide distribution across the globe. No effective vaccines exist for alphaviral infections. Therefore, a better understanding of VEEV and its associated neutralizing antibodies will help with the development of effective drugs and vaccines.

Discovery and evolution of bunyavirids in arctic phantom midges and ancient bunyavirid-like sequences in insect genomes.

Bunyaviridae is a large family of RNA viruses chiefly comprised of vertebrate and plant pathogens. We discovered novel bunyavirids that are approximately equally divergent from each of the five known genera. We characterized novel genome sequences for two bunyavirids, namely, Kigluaik phantom virus (KIGV), from tundra-native phantom midges (Chaoborus), and Nome phantom virus (NOMV), from tundra-invading phantom midges, and demonstrated that these bunyavirid-like sequences belong to an infectious virus by passaging KIGV in mosquito cell culture, although the infection does not seem to be well sustained beyond a few passages. Virus and host gene sequences from individuals collected on opposite ends of North America, a region spanning 4,000 km, support a long-term, vertically transmitted infection of KIGV in Chaoborus trivittatus. KIGV-like sequences ranging from single genes to full genomes are present in transcriptomes and genomes of insects belonging to six taxonomic orders, suggesting an ancient association of this clade with insect hosts. In Drosophila, endogenous virus genes have been coopted, forming an orthologous tandem gene family that has been maintained by selection during the radiation of the host genus. Our findings indicate that bunyavirid-host interactions in nonbloodsucking arthropods have been much more extensive than previously thought. Very little is known about the viral diversity in polar freshwater ponds, and perhaps less is known about the effects that climate-induced habitat changes in these regions will have on virus-host interactions in the coming years. Our results show that at the tundra-boreal boundary, a hidden viral landscape is being altered as infected boreal phantom midges colonize tundra ponds. Likewise, relatively little is known of the deeper evolutionary history of bunyavirids that has led to the stark lifestyle contrasts between some genera. The discovery of this novel bunyavirid group suggests that ancient and highly divergent bunyavirid lineages remain undetected in nature and may offer fresh insight into host reservoirs, potential sources of emerging disease, and major lifestyle shifts in the evolutionary history of viruses in the family Bunyaviridae.

Human immune responses to dengue virus infection: lessons learned from prospective cohort studies.

Dengue virus (DENV) continues to spread globally and is a major cause of morbidity and mortality. Currently, there is no antiviral treatment to diminish severe illness or a vaccine to induce protection from infection. An effective dengue vaccine that protects against all four DENV serotypes is a high priority for endemic countries and several candidates are in development by various United States Federal Agencies and private pharmaceutical companies. Challenges faced by dengue vaccine developers include creating tetravalent formulations that provide tetravalent protection, the lack of a correlate of protective immunity, a changing viral landscape as DENV evolves, and a complex viral-host pathogenesis that can result in a spectrum of illness from subclinical infection to severe hemorrhagic fever. There have been a number of long-term prospective studies on DENV transmission and dengue severity that have provided invaluable information on DENV epidemiology and pathogenesis of this disease. In this section, we will review the critical lessons learned from these studies and their application for dengue vaccine development.

A Viral Genome Landscape of RNA Polyadenylation from KSHV Latent to Lytic Infection

RNA polyadenylation (pA) is one of the major steps in regulation of gene expression at the posttranscriptional level. In this report, a genome landscape of pA sites of viral transcripts in B lymphocytes with Kaposi sarcoma-associated herpesvirus (KSHV) infection was constructed using a modified PA-seq strategy. We identified 67 unique pA sites, of which 55 could be assigned for expression of annotated ∼90 KSHV genes. Among the assigned pA sites, twenty are for expression of individual single genes and the rest for multiple genes (average 2.7 genes per pA site) in cluster-gene loci of the genome. A few novel viral pA sites that could not be assigned to any known KSHV genes are often positioned in the antisense strand to ORF8, ORF21, ORF34, K8 and ORF50, and their associated antisense mRNAs to ORF21, ORF34 and K8 could be verified by 3′RACE. The usage of each mapped pA site correlates to its peak size, the larger (broad and wide) peak size, the more usage and thus, the higher expression of the pA site-associated gene(s). Similar to mammalian transcripts, KSHV RNA polyadenylation employs two major poly(A) signals, AAUAAA and AUUAAA, and is regulated by conservation of cis-elements flanking the mapped pA sites. Moreover, we found two or more alternative pA sites downstream of ORF54, K2 (vIL6), K9 (vIRF1), K10.5 (vIRF3), K11 (vIRF2), K12 (Kaposin A), T1.5, and PAN genes and experimentally validated the alternative polyadenylation for the expression of KSHV ORF54, K11, and T1.5 transcripts. Together, our data provide not only a comprehensive pA site landscape for understanding KSHV genome structure and gene expression, but also the first evidence of alternative polyadenylation as another layer of posttranscriptional regulation in viral gene expression.

Heterosubtypic cross‐protection induced by whole inactivated influenza virus vaccine in mice: influence of the route of vaccine administration

Development of influenza vaccines capable of inducing broad protection against different virus subtypes is necessary given the ever-changing viral genetic landscape. Previously, we showed that vaccination with whole inactivated virus (WIV) induces heterosubtypic protection against lethal virus infection in mice. Whole inactivated virus-induced cross-protection was found to be mediated primarily by flu-specific CD8+ T cells. As it has been demonstrated that the route of vaccine administration strongly influences both the quantity and quality of vaccine-induced immunity, in this study, we determined which route of WIV administration induces optimal heterosubtypic cross-protection. We compared the magnitude of the immune response and heterosubtypic protection against lethal A/PR/8/34 (H1N1) infection after subcutaneous (SC), intramuscular (IM), and intranasal (IN) vaccination with A/NIBRG-14 (H5N1) WIV. Subcutaneous and IM administration was superior to IN administration of influenza WIV in terms of flu-specific CD8+ T-cell induction and protection of mice against lethal heterosubtypic challenge. Surprisingly, despite the very low flu-specific CD8+ T-cell responses detected in IN-vaccinated mice, these animals were partially protected, most likely due to cross-reactive IgA antibodies. The results of this study show that the magnitude of WIV-induced flu-specific CD8+ T-cell activity depends on the applied vaccination route. We conclude that parenteral administration of WIV vaccine, in particular IM injection, is superior to IN vaccine delivery for the induction of heterosubtypic cross-protection and generally appears to elicit stronger immune responses than mucosal vaccination with WIV.

The role of deformed wing virus in the initial collapse of varroa infested honey bee colonies in the UK

SummaryThe mite Varroa destructor has been associated with the collapse of millions of Apis mellifera honey bee colonies world-wide. During the past decade, a large body of research has revealed various interactions between varroa, the honey bee and various viral pathogens. One pathogen in particular, deformed wing virus (DWV), has emerged as the key pathogen involved in colony collapse. As varroa has permanently changed the viral landscape in which honey bees exist, we present a large body of data on the effects of DWV during the initial phase of varroa infestation in the UK during 1998. This provides baseline data for further comparative studies. We carried out DWV transmission studies, and observed the effects of DWV on bee longevity. As the ELISA technique used in these studies had a detection limit of ~107 viral particles per bee, only high viral (overt) titres were detected. During the initial phase of varroa establishment, DWV was detected in 0.6% of non-infested sealed brood, but in 89% of sealed brood invaded by varroa. Once DWV was introduced into the bee's haemolymph via mite feeding on either pupae or adults, an overt virus infection was rapidly produced in 3–4 days. In sealed brood the presence of varroa was fatal for 21% of the brood, caused wing deformity in some emerging adults and significantly reduced longevity as an adult. However, adult bees that became infected after they had emerged, did not develop wing deformity nor show any reduced longevity, but acted as reservoirs of DWV infection.

From HIV protein sequences to viral fitness landscapes: a new paradigm for in silico vaccine design

Background An inexpensive prophylactic vaccine offers the best hope to curb the HIV/AIDS epidemic gripping sub-Saharan Africa. Systematic means to guide the design of an effective immunogen for this, and other, infectious diseases are not available. What is required is a method to chart the peaks and valleys of viral fitness as a function of amino acid sequence. An efficacious vaccine would eject the virus from the high fitness peaks, and drive it into the valleys where its compromised fitness impairs its ability to replicate and inflict damage to the host. Methods Appealing to spin glass models in statistical physics, we present a novel approach to translate viral sequence databases into landscapes of viral fitness. These inferred models furnish a quantitative description of viral replicative capacity as a function of amino acid sequence. We illustrate this approach in the development of landscapes for the proteins of HIV-1 clade B Gag. Results In comparisons to experimental and clinical data, our inferred landscapes demonstrate excellent agreement with: 1) in vitro replicative fitness measurements, 2) clinically observed high-fitness circulating viral strains, 3) documented HLA associated CTL escape mutations, and 4) intrahost temporal adaptation pathways revealed by deep sequencing. These favorable comparisons support our landscapes as reflections of intrinsic viral fitness. We illustrate the value of such descriptions in the computational design of a CTL Gag immunogen. Conclusion We present a novel methodology to translate viral sequence data into quantitative landscapes of viral fitness. In an application to HIV-1 Gag, we illustrate excellent agreement of our model predictions with experimental and clinical data, and demonstrate a powerful new approach for HIV immunogen design. We anticipate that this approach may represent a heretofore unprecedented means to synthesize fitness landscapes for diverse pathogens, and may provide the basis for the design of improved prophylactic and therapeutic strategies.