Analyses of Early ZIKV Genomes Are Consistent with Viral Spread from Northeast Brazil to the Americas.

The Americas, particularly Brazil, were greatly impacted by the widespread Zika virus (ZIKV) outbreak in 2015 and 2016. Efforts were made to implement genomic surveillance of ZIKV as part of the public health responses. The accuracy of spatiotemporal reconstructions of the epidemic spread relies on the unbiased sampling of the transmission process. In the early stages of the outbreak, we recruited patients exhibiting clinical symptoms of arbovirus-like infection from Salvador and Campo Formoso, Bahia, in Northeast Brazil. Between May 2015 and June 2016, we identified 21 cases of acute ZIKV infection and subsequently recovered 14 near full-length sequences using the amplicon tiling multiplex approach with nanopore sequencing. We performed a time-calibrated discrete phylogeographic analysis to trace the spread and migration history of the ZIKV. Our phylogenetic analysis supports a consistent relationship between ZIKV migration from Northeast to Southeast Brazil and its subsequent dissemination beyond Brazil. Additionally, our analysis provides insights into the migration of ZIKV from Brazil to Haiti and the role Brazil played in the spread of ZIKV to other countries, such as Singapore, the USA, and the Dominican Republic. The data generated by this study enhances our understanding of ZIKV dynamics and supports the existing knowledge, which can aid in future surveillance efforts against the virus.

Genome-wide diversity of Zika virus: Exploring spatio-temporal dynamics to guide a new nomenclature proposal.

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized 'African' and 'Asian' genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

Increasing Prevalence of HIV-1 Transmitted Drug Resistance in Portugal: Implications for First Line Treatment Recommendations.

INTRODUCTION: Treatment for All recommendations have allowed access to antiretroviral (ARV) treatment for an increasing number of patients. This minimizes the transmission of infection but can potentiate the risk of transmitted (TDR) and acquired drug resistance (ADR). OBJECTIVE: To study the trends of TDR and ADR in patients followed up in Portuguese hospitals between 2001 and 2017. METHODS: In total, 11,911 patients of the Portuguese REGA database were included. TDR was defined as the presence of one or more surveillance drug resistance mutation according to the WHO surveillance list. Genotypic resistance to ARV was evaluated with Stanford HIVdb v7.0. Patterns of TDR, ADR and the prevalence of mutations over time were analyzed using logistic regression. RESULTS AND DISCUSSION: The prevalence of TDR increased from 7.9% in 2003 to 13.1% in 2017 (p < 0.001). This was due to a significant increase in both resistance to nucleotide reverse transcriptase inhibitors (NRTIs) and non-nucleotide reverse transcriptase inhibitors (NNRTIs), from 5.6% to 6.7% (p = 0.002) and 2.9% to 8.9% (p < 0.001), respectively. TDR was associated with infection with subtype B, and with lower viral load levels (p < 0.05). The prevalence of ADR declined from 86.6% in 2001 to 51.0% in 2017 (p < 0.001), caused by decreasing drug resistance to all antiretroviral (ARV) classes (p < 0.001). CONCLUSIONS: While ADR has been decreasing since 2001, TDR has been increasing, reaching a value of 13.1% by the end of 2017. It is urgently necessary to develop public health programs to monitor the levels and patterns of TDR in newly diagnosed patients.

Drivers of HIV-1 transmission: The Portuguese case.

BACKGROUND: Portugal has one of the most severe HIV-1 epidemics in Western Europe. Two subtypes circulate in parallel since the beginning of the epidemic. Comparing their transmission patterns and its association with transmitted drug resistance (TDR) is important to pinpoint transmission hotspots and to develop evidence-based treatment guidelines. METHODS: Demographic, clinical and genomic data were collected from 3599 HIV-1 naive patients between 2001 and 2014. Sequences obtained from drug resistance testing were used for subtyping, TDR determination and transmission clusters (TC) analyses. RESULTS: In Portugal, transmission of subtype B was significantly associated with young males, while transmission of subtype G was associated with older heterosexuals. In Portuguese originated people, there was a decreasing trend both for prevalence of subtype G and for number of TCs in this subtype. The active TCs that were identified (i.e. clusters originated after 2008) were associated with subtype B-infected males residing in Lisbon. TDR was significantly different when comparing subtypes B (10.8% [9.5-12.2]) and G (7.6% [6.4-9.0]) (p = 0.001). DISCUSSION: TC analyses shows that, in Portugal, the subtype B epidemic is active and fueled by young male patients residing in Lisbon, while transmission of subtype G is decreasing. Despite similar treatment rates for both subtypes in Portugal, TDR is significantly different between subtypes.

Advances in Visualization Tools for Phylogenomic and Phylodynamic Studies of Viral Diseases.

Genomic and epidemiological monitoring have become an integral part of our response to emerging and ongoing epidemics of viral infectious diseases. Advances in high-throughput sequencing, including portable genomic sequencing at reduced costs and turnaround time, are paralleled by continuing developments in methodology to infer evolutionary histories (dynamics/patterns) and to identify factors driving viral spread in space and time. The traditionally static nature of visualizing phylogenetic trees that represent these evolutionary relationships/processes has also evolved, albeit perhaps at a slower rate. Advanced visualization tools with increased resolution assist in drawing conclusions from phylogenetic estimates and may even have potential to better inform public health and treatment decisions, but the design (and choice of what analyses are shown) is hindered by the complexity of information embedded within current phylogenetic models and the integration of available meta-data. In this review, we discuss visualization challenges for the interpretation and exploration of reconstructed histories of viral epidemics that arose from increasing volumes of sequence data and the wealth of additional data layers that can be integrated. We focus on solutions that address joint temporal and spatial visualization but also consider what the future may bring in terms of visualization and how this may become of value for the coming era of real-time digital pathogen surveillance, where actionable results and adequate intervention strategies need to be obtained within days.

A computational method for the identification of Dengue, Zika and Chikungunya virus species and genotypes.

In recent years, an increasing number of outbreaks of Dengue, Chikungunya and Zika viruses have been reported in Asia and the Americas. Monitoring virus genotype diversity is crucial to understand the emergence and spread of outbreaks, both aspects that are vital to develop effective prevention and treatment strategies. Hence, we developed an efficient method to classify virus sequences with respect to their species and sub-species (i.e. serotype and/or genotype). This tool provides an easy-to-use software implementation of this new method and was validated on a large dataset assessing the classification performance with respect to whole-genome sequences and partial-genome sequences. Available online: http://krisp.org.za/tools.php.

An Evolutionary Model-Based Approach To Quantify the Genetic Barrier to Drug Resistance in Fast-Evolving Viruses and Its Application to HIV-1 Subtypes and Integrase Inhibitors.

Viral pathogens causing global disease burdens are often characterized by high rates of evolutionary changes. The extensive viral diversity at baseline can shorten the time to escape from therapeutic or immune selective pressure and alter mutational pathways. The impact of genotypic background on the barrier to resistance can be difficult to capture, particularly for agents in experimental stages or that are recently approved or expanded into new patient populations. We developed an evolutionary model-based counting method to quickly quantify the population genetic potential to resistance and assess population differences. We demonstrate its applicability to HIV-1 integrase inhibitors, as their increasing use globally contrasts with limited availability of non-B subtype resistant sequence data and corresponding knowledge gap. A large sequence data set encompassing most prevailing HIV-1 subtypes and resistance-associated mutations of currently approved integrase inhibitors was investigated. A complex interplay between codon predominance, polymorphisms, and associated evolutionary costs resulted in a subtype-dependent varied genetic potential for 15 resistance mutations against integrase inhibitors. While we confirm the lower genetic barrier of subtype B for G140S, we convincingly discard a similar effect previously suggested for G140C. A supplementary analysis for HIV-1 reverse transcriptase inhibitors identified a lower genetic barrier for K65R in subtype C through differential codon usage not reported before. To aid evolutionary interpretations of genomic differences for antiviral strategies, we advanced existing counting methods with increased sensitivity to identify subtype dependencies of resistance emergence. Future applications include novel HIV-1 drug classes or vaccines, as well as other viral pathogens.

SANTA-SIM: simulating viral sequence evolution dynamics under selection and recombination.

Simulations are widely used to provide expectations and predictive distributions under known conditions against which to compare empirical data. Such simulations are also invaluable for testing and comparing the behaviour and power of inference methods. We describe SANTA-SIM, a software package to simulate the evolution of a population of gene sequences forwards through time. It models the underlying biological processes as discrete components: replication, recombination, point mutations, insertion-deletions, and selection under various fitness models and population size dynamics. The software is designed to be intuitive to work with for a wide range of users and executable in a cross-platform manner.

VIRULIGN: fast codon-correct alignment and annotation of viral genomes.

SUMMARY: Virus sequence data are an essential resource for reconstructing spatiotemporal dynamics of viral spread as well as to inform treatment and prevention strategies. However, the potential benefit of these applications critically depends on accurate and correctly annotated alignments of genetically heterogeneous data. VIRULIGN was built for fast codon-correct alignments of large datasets, with standardized and formalized genome annotation and various alignment export formats. AVAILABILITY AND IMPLEMENTATION: VIRULIGN is freely available at https://github.com/rega-cev/virulign as an open source software project. SUPPLEMENTARY INFORMATION: Supplementary data is available at Bioinformatics online.

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

[This corrects the article DOI: 10.1371/journal.pgen.1007420.].

Time to Harmonize Dengue Nomenclature and Classification.

Dengue virus (DENV) is estimated to cause 390 million infections per year worldwide. A quarter of these infections manifest clinically and are associated with a morbidity and mortality that put a significant burden on the affected regions. Reports of increased frequency, intensity, and extended geographical range of outbreaks highlight the virus's ongoing global spread. Persistent transmission in endemic areas and the emergence in territories formerly devoid of transmission have shaped DENV's current genetic diversity and divergence. This genetic layout is hierarchically organized in serotypes, genotypes, and sub-genotypic clades. While serotypes are well defined, the genotype nomenclature and classification system lack consistency, which complicates a broader analysis of their clinical and epidemiological characteristics. We identify five key challenges: (1) Currently, there is no formal definition of a DENV genotype; (2) Two different nomenclature systems are used in parallel, which causes significant confusion; (3) A standardized classification procedure is lacking so far; (4) No formal definition of sub-genotypic clades is in place; (5) There is no consensus on how to report antigenic diversity. Therefore, we believe that the time is right to re-evaluate DENV genetic diversity in an essential effort to provide harmonization across DENV studies.

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.

A genetic IFN/STAT1/FAS axis determines CD4 T stem cell memory levels and apoptosis in healthy controls and Adult T-cell Leukemia patients.

Adult T-cell leukemia (ATL) is an aggressive, chemotherapy-resistant CD4+CD25+ leukemia caused by HTLV-1 infection, which usually develops in a minority of patients several decades after infection. IFN + AZT combination therapy has shown clinical benefit in ATL, although its mechanism of action remains unclear. We have previously shown that an IFN-responsive FAS promoter polymorphism in a STAT1 binding site (rs1800682) is associated to ATL susceptibility and survival. Recently, CD4 T stem cell memory (TSCM) Fashi cells have been identified as the hierarchical cellular apex of ATL, but a possible link between FAS, apoptosis, proliferation and IFN response in ATL has not been studied. In this study, we found significant ex vivo antiproliferative, antiviral and immunomodulatory effects of IFN-α treatment in short-term culture of primary mononuclear cells from ATL patients (n = 25). Bayesian Network analysis allowed us to integrate ex vivo IFN-α response with clinical, genetic and immunological data from ATL patients, thereby revealing a central role for FAS -670 polymorphism and apoptosis in the coordinated mechanism of action of IFN-α. FAS genotype-dependence of IFN-induced apoptosis was experimentally validated in an independent cohort of healthy controls (n = 20). The same FAS -670 polymorphism also determined CD4 TSCM levels in a genome-wide twin study (p = 7 × 10-11, n = 460), confirming a genetic link between apoptosis and TSCM levels. Transcriptomic analysis and cell type deconvolution confirmed the FAS genotype/TSCM link and IFN-α-induced downregulation of CD4 TSCM-specific genes in ATL patient cells. In conclusion, ex vivo IFN-α treatment exerts a pleiotropic effect on primary ATL cells, with a genetic IFN/STAT1/Fas axis determining apoptosis vs. proliferation and underscoring the CD4 TSCM model of ATL leukemogenesis.

HIV-1 Infection in Cyprus, the Eastern Mediterranean European Frontier: A Densely Sampled Transmission Dynamics Analysis from 1986 to 2012.

Since HIV-1 treatment is increasingly considered an effective preventionstrategy, it is important to study local HIV-1 epidemics to formulate tailored preventionpolicies. The prevalence of HIV-1 in Cyprus was historically low until 2005. To investigatethe shift in epidemiological trends, we studied the transmission dynamics of HIV-1 in Cyprususing a densely sampled Cypriot HIV-1 transmission cohort that included 85 percent ofHIV-1-infected individuals linked to clinical care between 1986 and 2012 based on detailedclinical, epidemiological, behavioral and HIV-1 genetic information. Subtyping andtransmission cluster reconstruction were performed using maximum likelihood and Bayesianmethods, and the transmission chain network was linked to the clinical, epidemiological andbehavioral data. The results reveal that for the main HIV-1 subtype A1 and B sub-epidemics,young and drug-naïve HIV-1-infected individuals in Cyprus are driving the dynamics of thelocal HIV-1 epidemic. The results of this study provide a better understanding of thedynamics of the HIV-1 infection in Cyprus, which may impact the development of preventionstrategies. Furthermore, this methodology for analyzing densely sampled transmissiondynamics is applicable to other geographic regions to implement effective HIV-1 preventionstrategies in local settings.

The impact of HIV-1 within-host evolution on transmission dynamics.

The adaptive potential of HIV-1 is a vital mechanism to evade host immune responses and antiviral treatment. However, high evolutionary rates during persistent infection can impair transmission efficiency and alter disease progression in the new host, resulting in a delicate trade-off between within-host virulence and between-host infectiousness. This trade-off is visible in the disparity in evolutionary rates at within-host and between-host levels, and preferential transmission of ancestral donor viruses. Understanding the impact of within-host evolution for epidemiological studies is essential for the design of preventive and therapeutic measures. Herein, we review recent theoretical and experimental work that generated new insights into the complex link between within-host evolution and between-host fitness, revealing temporal and selective processes underlying the structure and dynamics of HIV-1 transmission.

Characterization of Nucleoside Reverse Transcriptase Inhibitor-Associated Mutations in the RNase H Region of HIV-1 Subtype C Infected Individuals.

The South African national treatment programme includes nucleoside reverse transcriptase inhibitors (NRTIs) in both first and second line highly active antiretroviral therapy regimens. Mutations in the RNase H domain have been associated with resistance to NRTIs but primarily in HIV-1 subtype B studies. Here, we investigated the prevalence and association of RNase H mutations with NRTI resistance in sequences from HIV-1 subtype C infected individuals. RNase H sequences from 112 NRTI treated but virologically failing individuals and 28 antiretroviral therapy (ART)-naive individuals were generated and analysed. In addition, sequences from 359 subtype C ART-naive sequences were downloaded from Los Alamos database to give a total of 387 sequences from ART-naive individuals for the analysis. Fisher's exact test was used to identify mutations and Bayesian network learning was applied to identify novel NRTI resistance mutation pathways in RNase H domain. The mutations A435L, S468A, T470S, L484I, A508S, Q509L, L517I, Q524E and E529D were more prevalent in sequences from treatment-experienced compared to antiretroviral treatment naive individuals, however, only the E529D mutation remained significant after correction for multiple comparison. Our findings suggest a potential interaction between E529D and NRTI-treatment; however, site-directed mutagenesis is needed to understand the impact of this RNase H mutation.

Epidemic dispersion of HIV and HCV in a population of co-infected Romanian injecting drug users.

Co-infections with HIV and HCV are very frequent among people who inject drugs (PWID). However, very few studies comparatively reconstructed the transmission patterns of both viruses in the same population. We have recruited 117 co-infected PWID during a recent HIV outbreak in Romania. Phylogenetic analyses were performed on HIV and HCV sequences in order to characterize and compare transmission dynamics of the two viruses. Three large HIV clusters (2 subtype F1 and one CRF14_BG) and thirteen smaller HCV transmission networks (genotypes 1a, 1b, 3a, 4a and 4d) were identified. Eighty (65%) patients were both in HIV and HCV transmission chains and 70 of those shared the same HIV and HCV cluster with at least one other patient. Molecular clock analysis indicated that all identified HIV clusters originated around 2006, while the origin of the different HCV clusters ranged between 1980 (genotype 1b) and 2011 (genotypes 3a and 4d). HCV infection preceded HIV infection in 80.3% of cases. Coincidental transmission of HIV and HCV was estimated to be rather low (19.65%) and associated with an outbreak among PWID during detention in the same penitentiary. This study has reconstructed and compared the dispersion of these two viruses in a PWID population.

PhyloGeoTool: interactively exploring large phylogenies in an epidemiological context.

MOTIVATION: Clinicians, health officials and researchers are interested in the epidemic spread of pathogens in both space and time to support the optimization of intervention measures and public health policies. Large sequence databases of virus sequences provide an interesting opportunity to study this spread through phylogenetic analysis. To infer knowledge from large phylogenetic trees, potentially encompassing tens of thousands of virus strains, an efficient method for data exploration is required. The clades that are visited during this exploration should be annotated with strain characteristics (e.g. transmission risk group, tropism, drug resistance profile) and their geographic context. RESULTS: PhyloGeoTool implements a visual method to explore large phylogenetic trees and to depict characteristics of strains and clades, including their geographic context, in an interactive way. PhyloGeoTool also provides the possibility to position new virus strains relative to the existing phylogenetic tree, allowing users to gain insight in the placement of such new strains without the need to perform a de novo reconstruction of the phylogeny. AVAILABILITY AND IMPLEMENTATION: https://github.com/rega-cev/phylogeotool (Freely available: open source software project). CONTACT: phylogeotool@kuleuven.be. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exploring resistance pathways for first-generation NS3/4A protease inhibitors boceprevir and telaprevir using Bayesian network learning.

Resistance-associated variants (RAVs) have been shown to influence treatment response to direct-acting antivirals (DAAs) and first generation NS3/4A protease inhibitors (PIs) in particular. Interpretation of hepatitis C virus (HCV) genotypic drug resistance remains a challenge, especially in patients who previously failed DAA therapy and need to be retreated with a second DAA based regimen. Bayesian network (BN) learning on HCV sequence data from patients treated with DAAs could provide insight in resistance pathways against PIs for HCV subtypes 1a and 1b, in a similar way as applied before for HIV. The publicly available 'Rega-BN' tool chain was developed to study associative analyses for various pathogens. Our first analysis, comparing sequences from PI-naïve and PI-experienced patients, determined that NS3 substitutions R155K and V36M arise with PI-exposure in HCV1a infected patients, and were defined as major and minor resistance-associated variants respectively. NS3 variant 174H was newly identified as potentially related to PI resistance. In a second analysis, NS3 sequences from PI-naïve patients who cleared the virus during PI therapy and from PI-naïve patients who failed PI therapy were compared, showing that NS3 baseline variant 67S predisposes to treatment-failure and variant 72I to treatment success. This approach has the potential to better characterize the role of more RAVs, if sufficient therapy annotated sequence data becomes available in curated public databases. In addition, polymorphisms present in baseline sequences that predispose patients to therapy failure can be identified using this approach.