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

Sources of Human Immunodeficiency Virus Infections Among Men Who Have Sex With Men With a Migration Background: A Viral Phylogenetic Case Study in Amsterdam, The Netherlands

BACKGROUND

Men and women with a migration background comprise an increasing proportion of incident human immunodeficiency virus (HIV) cases across Western Europe.

METHODS

To characterize sources of transmission in local transmission chains, we used partial HIV consensus sequences with linked demographic and clinical data from the opt-out AIDS Therapy Evaluation in the Netherlands (ATHENA) cohort of people with HIV in the Netherlands and identified phylogenetically and epidemiologically possible HIV transmission pairs in Amsterdam. We interpreted these in the context of estimated infection dates, and quantified population-level sources of transmission to foreign-born and Dutch-born Amsterdam men who have sex with men (MSM) within Amsterdam transmission chains.

RESULTS

We estimate that Dutch-born MSM were the predominant sources of infections among all Amsterdam MSM who acquired their infection locally in 2010-2021, and among almost all foreign-born Amsterdam MSM subpopulations. Stratifying by 2-year intervals indicated time trends in transmission dynamics, with a majority of infections originating from foreign-born MSM since 2016, although uncertainty ranges remained wide.

CONCLUSIONS

Native-born MSM have predominantly driven HIV transmissions in Amsterdam in 2010-2021. However, in the context of rapidly declining incidence in Amsterdam, the contribution from foreign-born MSM living in Amsterdam is increasing, with some evidence that most local transmissions have been from foreign-born Amsterdam MSM since 2016.

Direct whole-genome sequencing of HIV-1 for clinical drug-resistance analysis and public health surveillance

BACKGROUND

Human Immunodeficiency virus type 1 (HIV-1) remains a significant global health threat partly due to its ability to develop resistance to anti-retroviral therapies. HIV-1 genotype and drug resistance analysis of the polymerase (pol) sequence is a mainstay of its clinical and public health management. However, as new treatments and resistances evolve, analysis methods must change accordingly. In this study, we outline the development and implementation of a direct whole-genome sequencing approach (dWGS) using probe-capture target-enrichment for HIV-1 genotype and drug resistance analysis.

METHODS

We implemented dWGS and performed parallel pol Sanger sequencing for clinical samples, followed by comparative genotype and drug-resistance analysis. These HIV-1 WGS sequences were also utilised for a novel partitioned phylogenetic analysis.

RESULTS

Optimised nucleic acid extraction and DNAse I treatment significantly increased HIV-1 whole-genome coverage and depth, and improved recovery of high-quality genomes from low viral load clinical samples, enabling routine sequencing of viral loads as low as 1000 copies/mL. Overall, dWGS was robust, accurate and more sensitive for detecting low-frequency variants at drug-resistance sites compared to Sanger sequencing. Analysis of multiple sequence regions improved phylogenetic reconstruction for recombinant HIV-1 sequences compared to analysis of pol sequence alone.

CONCLUSIONS

These findings demonstrate dWGS enhances HIV-1 drug-resistance analysis by quantitative variant detection and improves reconstruction of HIV-1 phylogenies compared to traditional pol sequencing. This work supports that HIV-1 dWGS is a viable option to replace Sanger sequencing for clinical and public health applications.

MARVEL-minimising the emergence and dissemination of HIV-1 drug resistance in Portuguese-speaking African Countries (PALOP): low-cost portable NGS platform for HIV-1 surveillance in Africa

BACKGROUND

HIV-1 infections remain a global public health concern. Scaled-up antiretroviral treatment (ART) is crucial for reducing morbidity and mortality related to HIV/AIDS. The emergence of drug-resistance mutations (DRMs) compromises viral suppression and contributes to the continued HIV-1 transmission. Several reports indicate a recent increase in acquired (ADR) and transmitted (TDR) drug resistance in Africa, probably linked to the lack of implementation of HIV drug resistance (HIVDR) testing and suboptimal treatment adherence. Herein, we will develop a low-cost protocol using third-generation sequencing (Oxford Nanopore Technology) for HIV-1 surveillance in Portuguese-speaking African Countries - PALOP [Angola (AO), Cape Verde (CV), Mozambique (MZ), and Sao Tome & Principe (STP)].

METHODS

This is a multicentric cross-sectional study that includes around 600 adult patients newly diagnosed with HIV-1 in the PALOP. An epidemiological questionnaire previously validated by our research team will be used to collect sociodemographic and clinical data. Also, whole blood samples will be collected and the plasma samples will be subjected to drug resistance testing using an in-house low-cost NGS protocol. Data analysis will involve bioinformatics, biostatistics and machine learning techniques to generate accurate and up-to-date information about HIV-1 genetic diversity, ADR and TDR.

DISCUSSION

The implementation of this low-cost NGS platform for HIV-1 surveillance in the PALOP will allow: (i) to increase DRM surveillance capacity in resource-limited settings; (ii) to understand the pattern and determinants of dissemination of resistant HIV-1 strains; and (iii) to promote the development of technical and scientific skills of African researchers for genomic surveillance of viral pathogens and bioinformatics analysis. These objectives will contribute to reinforcing the capacity to combat HIV infection in Africa by optimizing the selection of ART regimens, improving viral suppression, and reducing ADR or TDR prevalence in PALOPs, with relevant implications for public health.

HIV-1 diversity and pre-treatment drug resistance in the era of integrase inhibitor among newly diagnosed ART-naïve adult patients in Luanda, Angola

The surveillance of drug resistance in the HIV-1 naïve population remains critical to optimizing the effectiveness of antiretroviral therapy (ART), mainly in the era of integrase strand transfer inhibitor (INSTI) regimens. Currently, there is no data regarding resistance to INSTI in Angola since Dolutegravir-DTG was included in the first-line ART regimen. Herein, we investigated the HIV-1 genetic diversity and pretreatment drug resistance (PDR) profile against nucleoside/tide reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), protease inhibitors (PIs), and INSTIs, using a next-generation sequencing (NGS) approach with MinION, established to track and survey DRMs in Angola. This was a cross-sectional study comprising 48 newly HIV-diagnosed patients from Luanda, Angola, screened between March 2022 and May 2023. PR, RT, and IN fragments were sequenced for drug resistance and molecular transmission cluster analysis. A total of 45 out of the 48 plasma samples were successfully sequenced. Of these, 10/45 (22.2%) presented PDR to PIs/NRTIs/NNRTIs. Major mutations for NRTIs (2.2%), NNRTIs (20%), PIs (2.2%), and accessory mutations against INSTIs (13.3%) were detected. No major mutations against INSTIs were detected. M41L (2%) and I85V (2%) mutations were detected for NRTI and PI, respectively. K103N (7%), Y181C (7%), and K101E (7%) mutations were frequently observed in NNRTI. The L74M (9%) accessory mutation was frequently observed in the INSTI class. HIV-1 pure subtypes C (33%), F1 (17%), G (15%), A1 (10%), H (6%), and D (4%), CRF01_AG (4%) were observed, while about 10% were recombinant strains. About 31% of detected HIV-1C sequences were in clusters, suggesting small-scale local transmission chains. No major mutations against integrase inhibitors were detected, supporting the continued use of INSTI in the country. Further studies assessing the HIV-1 epidemiology in the era of INSTI-based ART regimens are needed in Angola.

UnCoVar: a reproducible and scalable workflow for transparent and robust virus variant calling and lineage assignment using SARS-CoV-2 as an example

BACKGROUND

At a global scale, the SARS-CoV-2 virus did not remain in its initial genotype for a long period of time, with the first global reports of variants of concern (VOCs) in late 2020. Subsequently, genome sequencing has become an indispensable tool for characterizing the ongoing pandemic, particularly for typing SARS-CoV-2 samples obtained from patients or environmental surveillance. For such SARS-CoV-2 typing, various in vitro and in silico workflows exist, yet to date, no systematic cross-platform validation has been reported.

RESULTS

In this work, we present the first comprehensive cross-platform evaluation and validation of in silico SARS-CoV-2 typing workflows. The evaluation relies on a dataset of 54 patient-derived samples sequenced with several different in vitro approaches on all relevant state-of-the-art sequencing platforms. Moreover, we present UnCoVar, a robust, production-grade reproducible SARS-CoV-2 typing workflow that outperforms all other tested approaches in terms of precision and recall.

CONCLUSIONS

In many ways, the SARS-CoV-2 pandemic has accelerated the development of techniques and analytical approaches. We believe that this can serve as a blueprint for dealing with future pandemics. Accordingly, UnCoVar is easily generalizable towards other viral pathogens and future pandemics. The fully automated workflow assembles virus genomes from patient samples, identifies existing lineages, and provides high-resolution insights into individual mutations. UnCoVar includes extensive quality control and automatically generates interactive visual reports. UnCoVar is implemented as a Snakemake workflow. The open-source code is available under a BSD 2-clause license at github.com/IKIM-Essen/uncovar.

Applying Next-Generation Sequencing to Track HIV-1 Drug Resistance Mutations Circulating in Portugal

BACKGROUND

The global scale-up of antiretroviral treatment (ART) offers significant health benefits by suppressing HIV-1 replication and increasing CD4 cell counts. However, incomplete viral suppression poses a potential threat for the emergence of drug resistance mutations (DRMs), limiting ART options, and increasing HIV transmission.

OBJECTIVE

We investigated the patterns of transmitted drug resistance (TDR) and acquired drug resistance (ADR) among HIV-1 patients in Portugal.

METHODS

Data were obtained from 1050 HIV-1 patient samples submitted for HIV drug resistance (HIVDR) testing from January 2022 to June 2023. Evaluation of DRM affecting viral susceptibility to nucleoside/tide reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), protease inhibitors (PIs), and integrase strand transfer inhibitors (INSTIs) was performed using an NGS technology, the Vela Diagnostics Sentosa SQ HIV-1 Genotyping Assay.

RESULTS

About 71% of patients were ART naïve and 29% were experienced. Overall, 20% presented with any DRM. The prevalence of TDR and ADR was 12.6% and 41.1%, respectively. M184V, T215S, and M41L mutations for NRTI, K103N for NNRTI, and M46I/L for PIs were frequent in naïve and treated patients. E138K and R263K mutations against INSTIs were more frequent in naïve than treated patients. TDR and ADR to INSTIs were 0.3% and 7%, respectively. Patients aged 50 or over (OR: 1.81, p = 0.015), originating from Portuguese-speaking African countries (PALOPs) (OR: 1.55, p = 0.050), HIV-1 subtype G (OR: 1.78, p = 0.010), and with CD4 < 200 cells/mm3 (OR: 1.70, p = 0.043) were more likely to present with DRMs, while the males (OR: 0.63, p = 0.003) with a viral load between 4.1 to 5.0 Log10 (OR: 0.55, p = 0.003) or greater than 5.0 Log10 (OR: 0.52, p < 0.001), had lower chances of presenting with DRMs.

CONCLUSIONS

We present the first evidence on TDR and ADR to INSTI regimens in followed up patients presenting for healthcare in Portugal. We observed low levels of TDR to INSTIs among ART-naïve and moderate levels in ART-exposed patients. Regimens containing PIs could be an alternative second line in patients with intermediate or high-level drug resistance, especially against second-generation INSTIs (dolutegravir, bictegravir, and cabotegravir).

CNCA aligns small annotated genomes

BACKGROUND

To explore the evolutionary history of sequences, a sequence alignment is a first and necessary step, and its quality is crucial. In the context of the study of the proximal origins of SARS-CoV-2 coronavirus, we wanted to construct an alignment of genomes closely related to SARS-CoV-2 using both coding and non-coding sequences. To our knowledge, there is no tool that can be used to construct this type of alignment, which motivated the creation of CNCA.

RESULTS

CNCA is a web tool that aligns annotated genomes from GenBank files. It generates a nucleotide alignment that is then updated based on the protein sequence alignment. The output final nucleotide alignment matches the protein alignment and guarantees no frameshift. CNCA was designed to align closely related small genome sequences up to 50 kb (typically viruses) for which the gene order is conserved.

CONCLUSIONS

CNCA constructs multiple alignments of small genomes by integrating both coding and non-coding sequences. This preserves regions traditionally ignored in conventional back-translation methods, such as non-coding regions.

The Role of Late Presenters in HIV-1 Transmission Clusters in Europe

BACKGROUND

Investigating the role of late presenters (LPs) in HIV-1 transmission is important, as they can contribute to the onward spread of HIV-1 virus before diagnosis, when they are not aware of their HIV status.

OBJECTIVE

To characterize individuals living with HIV-1 followed up in Europe infected with subtypes A, B, and G and to compare transmission clusters (TC) in LP vs. non-late presenter (NLP) populations.

METHODS

Information from a convenience sample of 2679 individuals living with HIV-1 was collected from the EuResist Integrated Database between 2008 and 2019. Maximum likelihood (ML) phylogenies were constructed using FastTree. Transmission clusters were identified using Cluster Picker. Statistical analyses were performed using R.

RESULTS

2437 (91.0%) sequences were from subtype B, 168 (6.3%) from subtype A, and 74 (2.8%) from subtype G. The median age was 39 y/o (IQR: 31.0-47.0) and 85.2% of individuals were males. The main transmission route was via homosexual (MSM) contact (60.1%) and 85.0% originated from Western Europe. In total, 54.7% of individuals were classified as LPs and 41.7% of individuals were inside TCs. In subtype A, individuals in TCs were more frequently males and natives with a recent infection. For subtype B, individuals in TCs were more frequently individuals with MSM transmission route and with a recent infection. For subtype G, individuals in TCs were those with a recent infection. When analyzing cluster size, we found that LPs more frequently belonged to small clusters (<8 individuals), particularly dual clusters (2 individuals).

CONCLUSION

LP individuals are more present either outside or in small clusters, indicating a limited role of late presentation to HIV-1 transmission.

Multiple sequence alignment based on deep reinforcement learning with self-attention and positional encoding

MOTIVATION

Multiple sequence alignment (MSA) is one of the hotspots of current research and is commonly used in sequence analysis scenarios. However, there is no lasting solution for MSA because it is a Nondeterministic Polynomially complete problem, and the existing methods still have room to improve the accuracy.

RESULTS

We propose Deep reinforcement learning with Positional encoding and self-Attention for MSA, based on deep reinforcement learning, to enhance the accuracy of the alignment Specifically, inspired by the translation technique in natural language processing, we introduce self-attention and positional encoding to improve accuracy and reliability. Firstly, positional encoding encodes the position of the sequence to prevent the loss of nucleotide position information. Secondly, the self-attention model is used to extract the key features of the sequence. Then input the features into a multi-layer perceptron, which can calculate the insertion position of the gap according to the features. In addition, a novel reinforcement learning environment is designed to convert the classic progressive alignment into progressive column alignment, gradually generating each column's sub-alignment. Finally, merge the sub-alignment into the complete alignment. Extensive experiments based on several datasets validate our method's effectiveness for MSA, outperforming some state-of-the-art methods in terms of the Sum-of-pairs and Column scores.

AVAILABILITY AND IMPLEMENTATION

The process is implemented in Python and available as open-source software from https://github.com/ZhangLab312/DPAMSA.

Investigation of the individual genetic evolution of SARS-CoV-2 in a small cluster during the rapid spread of the BF.5 lineage in Tokyo, Japan

There has been a decreasing trend in new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases and fatalities worldwide. The virus has been evolving, indicating the potential emergence of new variants and uncertainties. These challenges necessitate continued efforts in disease control and mitigation strategies. We investigated a small cluster of SARS-CoV-2 Omicron variant infections containing a common set of genomic mutations, which provided a valuable model for investigating the transmission mechanism of genetic alterations. We conducted a study at a medical center in Japan during the Omicron surge (sub-lineage BA.5), sequencing the entire SARS-CoV-2 genomes from infected individuals and evaluating the phylogenetic tree and haplotype network among the variants. We compared the mutations present in each strain within the BA.5 strain, TKYnat2317, which was first identified in Tokyo, Japan. From June 29th to July 4th 2022, nine healthcare workers (HCWs) tested positive for SARS-CoV-2 by real-time PCR. During the same period, five patients also tested positive by real-time PCR. Whole genome sequencing revealed that the infected patients belonged to either the isolated BA.2 or BA.5 sub-lineage, while the healthcare worker infections were classified as BF.5. The phylogenetic tree and haplotype network clearly showed the specificity and similarity of the HCW cluster. We identified 12 common mutations in the cluster, including I110V in nonstructural protein 4 (nsp4), A1020S in the Spike protein, and H47Y in ORF7a, compared to the BA.5 reference. Additionally, one case had the extra nucleotide-deletion mutation I27* in ORF10, and low frequencies of genetic alterations were also found in certain instances. The results of genome sequencing showed that the nine HCWs shared a set of genetic mutations, indicating transmission within the cluster. Minor mutations observed in five HCW individuals suggested the emergence of new virus variants. Five amino acid substitutions occurred in nsp3, which could potentially affect virus replication or immune escape. Intra-host evolution also generated additional mutations. The cluster exhibited a mild disease course, with individuals in this case, recovering without requiring any medical treatments. Further investigation is needed to understand the relationship between the genetic evolution of the virus and the symptoms.

Use of whole genome sequencing to identify low-frequency mutations in SARS-CoV-2 patients treated with remdesivir

BACKGROUND

Remdesivir (RDV) has been shown to reduce hospitalization and mortality in COVID-19 patients. Resistance mutations caused by RDV are rare and have been predominantly reported in patients who are on prolonged therapy and immunocompromised. We investigate the effects of RDV treatment on intra-host SARS-CoV-2 diversity and low-frequency mutations in moderately ill hospitalized COVID-19 patients and compare them to patients without RDV treatment.

METHODS

From March 2020 to April 2022, sequential collections of nasopharyngeal and mid-turbinate swabs were obtained from 14 patients with and 30 patients without RDV treatment. Demographic and clinical data on all patients were reviewed. A total of 109 samples were sequenced and mutation analyses were performed.

RESULTS

Previously reported drug resistant mutations in nsp12 were not identified during short courses of RDV therapy. In genes encoding and surrounding the replication complex (nsp6-nsp14), low-frequency minority variants were detected in 7/14 (50%) and 18/30 (60%) patients with and without RDV treatment, respectively. We did not detect significant differences in within-host diversity and positive selection between the RDV-treated and untreated groups.

CONCLUSIONS

Minimal intra-host variability and stochastic low-frequency variants detected in moderately ill patients suggests little selective pressure in patients receiving short courses of RDV. The barrier to RDV resistance is high in patients with moderate disease. Patients undergoing short regimens of RDV therapy should continue to be monitored.

Reconstruction of SARS-CoV-2 outbreaks in a primary school using epidemiological and genomic data

Mathematical modelling studies have shown that repetitive screening can be used to mitigate SARS-CoV-2 transmission in primary schools while keeping schools open. However, not much is known about how transmission progresses within schools and whether there is a risk of importation to households. During the academic year 2020-2021, a prospective surveillance study using repetitive screening was conducted in a primary school and associated households in Liège (Belgium). SARS-CoV-2 screening was performed via throat washing either once or twice a week. We used genomic and epidemiological data to reconstruct the observed school outbreaks using two different models. The outbreaker2 model combines information on the generation time and contact patterns with a model of sequence evolution. For comparison we also used SCOTTI, a phylogenetic model based on the structured coalescent. In addition, we performed a simulation study to investigate how the accuracy of estimated positivity rates in a school depends on the proportion of a school that is sampled in a repetitive screening strategy. We found no difference in SARS-CoV-2 positivity between children and adults and children were not more often asymptomatic compared to adults. Both models for outbreak reconstruction revealed that transmission occurred mainly within the school environment. Uncertainty in outbreak reconstruction was lowest when including genomic as well as epidemiological data. We found that observed weekly positivity rates are a good approximation to the true weekly positivity rate, especially in children, even when only 25% of the school population is sampled. These results indicate that, in addition to reducing infections as shown in modelling studies, repetitive screening in school settings can lead to a better understanding of the extent of transmission in schools during a pandemic and importation risk at the community level.

COVID Variants, Villain and Victory: A Bioinformatics Perspective

The SARS-CoV-2 virus, a novel member of the Coronaviridae family, is responsible for the viral infection known as Coronavirus Disease 2019 (COVID-19). In response to the urgent and critical need for rapid detection, diagnosis, analysis, interpretation, and treatment of COVID-19, a wide variety of bioinformatics tools have been developed. Given the virulence of SARS-CoV-2, it is crucial to explore the pathophysiology of the virus. We intend to examine how bioinformatics, in conjunction with next-generation sequencing techniques, can be leveraged to improve current diagnostic tools and streamline vaccine development for emerging SARS-CoV-2 variants. We also emphasize how bioinformatics, in general, can contribute to critical areas of biomedicine, including clinical diagnostics, SARS-CoV-2 genomic surveillance and its evolution, identification of potential drug targets, and development of therapeutic strategies. Currently, state-of-the-art bioinformatics tools have helped overcome technical obstacles with respect to genomic surveillance and have assisted in rapid detection, diagnosis, and delivering precise treatment to individuals on time.

Mapping the Dynamics of Contemporary PRRSV-2 Evolution and Its Emergence and Spreading Hotspots in the U.S. Using Phylogeography

The repeated emergence of new genetic variants of PRRSV-2, the virus that causes porcine reproductive and respiratory syndrome (PRRS), reflects its rapid evolution and the failure of previous control efforts. Understanding spatiotemporal heterogeneity in variant emergence and spread is critical for future outbreak prevention. Here, we investigate how the pace of evolution varies across time and space, identify the origins of sub-lineage emergence, and map the patterns of the inter-regional spread of PRRSV-2 Lineage 1 (L1)-the current dominant lineage in the U.S. We performed comparative phylogeographic analyses on subsets of 19,395 viral ORF5 sequences collected across the U.S. and Canada between 1991 and 2021. The discrete trait analysis of multiple spatiotemporally stratified sampled sets (n = 500 each) was used to infer the ancestral geographic region and dispersion of each sub-lineage. The robustness of the results was compared to that of other modeling methods and subsampling strategies. Generally, the spatial spread and population dynamics varied across sub-lineages, time, and space. The Upper Midwest was a main spreading hotspot for multiple sub-lineages, e.g., L1C and L1F, though one of the most recent emergence events (L1A(2)) spread outwards from the east. An understanding of historical patterns of emergence and spread can be used to strategize disease control and the containment of emerging variants.

Genetic Characterization of a Novel Equus caballus Papillomavirus Isolated from a Thoroughbred Mare

Papillomaviruses (PVs) are small, non-enveloped viruses, ubiquitous across the animal kingdom. PVs induce diverse forms of infection, such as cutaneous papillomas, genital papillomatosis, and carcinomas. During a survey on the fertility status of a mare, a novel Equus caballus PV (EcPV) has been identified using Next Generation Sequencing, and it was further confirmed with genome-walking PCR and Sanger sequencing. The complete circular genome 7607 bp long shares 67% average percentage of identity with EcPV9, EcPV2, EcPV1, and EcPV6, justifying a new classification as Equus caballus PV 10 (EcPV10). All EcPV genes are conserved in EcPV10, and phylogenetic analysis indicates that EcPV10 is closely related to EcPV9 and EcPV2, genus Dyoiota 1. A preliminary EcPV10 genoprevalence study, carried out on 216 horses using Real Time PCRs, suggested a low incidence of this isolate (3.7%) compared to EcPVs of the same genus such as EcPV2 and EcPV9 in the same horse population. We hypothesize a transmission mechanism different from the one observed in the closely related EcPV9 and EcPV2 that particularly infect Thoroughbreds. This horse breed is usually submitted to natural mating, thus indicating a possible sexual diffusion. No differences were detected for breeds in terms of susceptibility to EcPV10. Further studies are needed to investigate the molecular mechanisms behind the host and EcPV10 infection to explain the reduced viral spread.

VPrimer: A Method of Designing and Updating Primer and Probe With High Variant Coverage for RNA Virus Detection

Fatal infectious diseases caused by RNA viruses, such as COVID-19, have emerged around the world. RT-PCR is widely employed for virus detection, and its accuracy depends on the primers and probes since RT-PCR can detect a virus only when the primers and probes bind to the target gene of the virus. Most of primer design methods are for a single host and so require a great deal of effort to design for RNA virus detection, including homology tests among the host and all the viruses for the host using BLAST-like tools. Furthermore, they do not consider variant sequences, which are very common in viruses. In this study, we describe VPrimer, a method of designing high-quality primer-probe sets for RNA viruses. VPrimer can find primer-probe sets that cover more than 95% of the variants of a target virus but do not cover any sequences of other viruses or the host. With VPrimer, we found 381,698,582 primer-probe sets for 3,104 RNA viruses. Multiplex PCR assays using the top 2 primer-probe sets suggested by VPrimer usually cover 100% of variants. To address the rapid changes in viral genomes, VPrimer finds the best and up-to-date primer-probe sets incrementally against the most recently reported variants.

Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission

Wildlife reservoirs of broad-host-range viruses have the potential to enable evolution of viral variants that can emerge to infect humans. In North America, there is phylogenomic evidence of continual transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from humans to white-tailed deer (Odocoileus virginianus) through unknown means, but no evidence of transmission from deer to humans. We carried out an observational surveillance study in Ontario, Canada during November and December 2021 (n = 300 deer) and identified a highly divergent lineage of SARS-CoV-2 in white-tailed deer (B.1.641). This lineage is one of the most divergent SARS-CoV-2 lineages identified so far, with 76 mutations (including 37 previously associated with non-human mammalian hosts). From a set of five complete and two partial deer-derived viral genomes we applied phylogenomic, recombination, selection and mutation spectrum analyses, which provided evidence for evolution and transmission in deer and a shared ancestry with mink-derived virus. Our analysis also revealed an epidemiologically linked human infection. Taken together, our findings provide evidence for sustained evolution of SARS-CoV-2 in white-tailed deer and of deer-to-human transmission.

Estimating the potential to prevent locally acquired HIV infections in a UNAIDS Fast-Track City, Amsterdam

Background

More than 300 cities including the city of Amsterdam in the Netherlands have joined the UNAIDS Fast-Track Cities initiative, committing to accelerate their HIV response and end the AIDS epidemic in cities by 2030. To support this commitment, we aimed to estimate the number and proportion of Amsterdam HIV infections that originated within the city, from Amsterdam residents. We also aimed to estimate the proportion of recent HIV infections during the 5-year period 2014-2018 in Amsterdam that remained undiagnosed.

Methods

We located diagnosed HIV infections in Amsterdam using postcode data (PC4) at time of registration in the ATHENA observational HIV cohort, and used HIV sequence data to reconstruct phylogeographically distinct, partially observed Amsterdam transmission chains. Individual-level infection times were estimated from biomarker data, and used to date the phylogenetically observed transmission chains as well as to estimate undiagnosed proportions among recent infections. A Bayesian Negative Binomial branching process model was used to estimate the number, size, and growth of the unobserved Amsterdam transmission chains from the partially observed phylogenetic data.

Results

Between 1 January 2014 and 1 May 2019, there were 846 HIV diagnoses in Amsterdam residents, of whom 516 (61%) were estimated to have been infected in 2014-2018. The rate of new Amsterdam diagnoses since 2014 (104 per 100,000) remained higher than the national rates excluding Amsterdam (24 per 100,000), and in this sense Amsterdam remained a HIV hotspot in the Netherlands. An estimated 14% [12-16%] of infections in Amsterdan MSM in 2014-2018 remained undiagnosed by 1 May 2019, and 41% [35-48%] in Amsterdam heterosexuals, with variation by region of birth. An estimated 67% [60-74%] of Amsterdam MSM infections in 2014-2018 had an Amsterdam resident as source, and 56% [41-70%] in Amsterdam heterosexuals, with heterogeneity by region of birth. Of the locally acquired infections, an estimated 43% [37-49%] were in foreign-born MSM, 41% [35-47%] in Dutch-born MSM, 10% [6-18%] in foreign-born heterosexuals, and 5% [2-9%] in Dutch-born heterosexuals. We estimate the majority of Amsterdam MSM infections in 2014-2018 originated in transmission chains that pre-existed by 2014.

Conclusions

This combined phylogenetic, epidemiologic, and modelling analysis in the UNAIDS Fast-Track City Amsterdam indicates that there remains considerable potential to prevent HIV infections among Amsterdam residents through city-level interventions. The burden of locally acquired infection remains concentrated in MSM, and both Dutch-born and foreign-born MSM would likely benefit most from intensified city-level interventions.

Funding

This study received funding as part of the H-TEAM initiative from Aidsfonds (project number P29701). The H-TEAM initiative is being supported by Aidsfonds (grant number: 2013169, P29701, P60803), Stichting Amsterdam Dinner Foundation, Bristol-Myers Squibb International Corp. (study number: AI424-541), Gilead Sciences Europe Ltd (grant number: PA-HIV-PREP-16-0024), Gilead Sciences (protocol numbers: CO-NL-276-4222, CO-US-276-1712, CO-NL-985-6195), and M.A.C AIDS Fund.

Genomic surveillance of Rift Valley fever virus: from sequencing to lineage assignment

Genetic evolution of Rift Valley fever virus (RVFV) in Africa has been shaped mainly by environmental changes such as abnormal rainfall patterns and climate change that has occurred over the last few decades. These gradual environmental changes are believed to have effected gene migration from macro (geographical) to micro (reassortment) levels. Presently, 15 lineages of RVFV have been identified to be circulating within the Sub-Saharan Africa. International trade in livestock and movement of mosquitoes are thought to be responsible for the outbreaks occurring outside endemic or enzootic regions. Virus spillover events contribute to outbreaks as was demonstrated by the largest epidemic of 1977 in Egypt. Genomic surveillance of the virus evolution is crucial in developing intervention strategies. Therefore, we have developed a computational tool for rapidly classifying and assigning lineages of the RVFV isolates. The computational method is presented both as a command line tool and a web application hosted at https://www.genomedetective.com/app/typingtool/rvfv/. Validation of the tool has been performed on a large dataset using glycoprotein gene (Gn) and whole genome sequences of the Large (L), Medium (M) and Small (S) segments of the RVFV retrieved from the National Center for Biotechnology Information (NCBI) GenBank database. Using the Gn nucleotide sequences, the RVFV typing tool was able to correctly classify all 234 RVFV sequences at species level with 100% specificity, sensitivity and accuracy. All the sequences in lineages A (n = 10), B (n = 1), C (n = 88), D (n = 1), E (n = 3), F (n = 2), G (n = 2), H (n = 105), I (n = 2), J (n = 1), K (n = 4), L (n = 8), M (n = 1), N (n = 5) and O (n = 1) were also correctly classified at phylogenetic level. Lineage assignment using whole RVFV genome sequences (L, M and S-segments) did not achieve 100% specificity, sensitivity and accuracy for all the sequences analyzed. We further tested our tool using genomic data that we generated by sequencing 5 samples collected following a recent RVF outbreak in Kenya. All the 5 samples were assigned lineage C by both the partial (Gn) and whole genome sequence classifiers. The tool is useful in tracing the origin of outbreaks and supporting surveillance efforts.

Availability: https://github.com/ajodeh-juma/rvfvtyping

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

Identification of conserved regions from 230,163 SARS-CoV-2 genomes and their use in diagnostic PCR primer design

BACKGROUND

As the rapidly evolving characteristic of SARS-CoV-2 could result in false negative diagnosis, the use of as much sequence data as possible is key to the identification of conserved viral sequences. However, multiple alignment of massive genome sequences is computationally intensive.

OBJECTIVE

To extract conserved sequences from SARS-CoV-2 genomes for the design of diagnostic PCR primers using a bioinformatics approach that can handle massive genomic sequences efficiently.

METHODS

A total of 230,163 full-length viral genomes were retrieved from the NCBI SARS-CoV-2 Resources and GISAID EpiCoV database. This number was reduced to 14.11% following removal of 5'-/3'-untranslated regions and sequence dereplication. Fast, reference-based, multiple sequence alignments identified conserved sequences and specific primer sets were designed against these regions using a conventional tool. Primer sets chosen among the candidates were evaluated by in silico PCR and RT-qPCR.

RESULTS

Out of 17 conserved sequences (totaling 4.3 kb), two primer sets targeting the nsp2 and ORF3a genes were picked that exhibited > 99.9% in silico amplification coverage against the original dataset (230,163 genomes) when a 5% mismatch between the primers and target was allowed. In addition, the primer sets successfully detected nine SARS-CoV-2 variant RNA samples (Alpha, Beta, Gamma, Delta, Epsilon, Zeta, Eta, Iota, and Kappa) in experimental RT-qPCR validations.

CONCLUSION

In addition to the RdRp, E, N, and S genes that are targeted commonly, our approach can be used to identify novel primer targets in SARS-CoV-2 and should be a priority strategy in the event of novel SARS-CoV-2 variants or other pandemic outbreaks.

Differential patterns of postmigration HIV-1 infection acquisition among Portuguese immigrants of different geographical origins

OBJECTIVE

To investigate the dynamics of phylogenetic transmission clusters involving immigrants of Portuguese Speaking Countries living in Portugal.

DESIGN/METHODS

We included genomic sequences, sociodemographic and clinical data from 772 HIV migrants followed in Portugal between 2001 and 2017. To reconstruct HIV-1 transmission clusters, we applied phylogenetic inference from 16 454 patients: 772 migrants, 2973 Portuguese and 12 709 global controls linked to demographic and clinical data. Transmission clusters were defined using: clusters with SH greater than 90% (phylogenetic support), genetic distance less than 3.5% and clusters that included greater than 66% of patients from one specific geographic origin compared with the total of sequences within the cluster. Logistic regression was performed to assess factors associated with clustering.

RESULTS

Three hundred and six (39.6%) of migrants were included in transmission clusters. This proportion differed substantially by region of origin [Brazil 54% vs. Portuguese Speaking African Countries (PALOPs) 36%, P < 0.0001] and HIV-1 infecting subtype (B 52%, 43% subtype G and 32% CRF02_AG, P < 0.001). Belonging to a transmission cluster was independently associated with treatment-naive patients, CD4+ greater than 500, with recent calendar years of sampling, origin from PALOPs and with seroconversion. Among Brazilian migrants - mainly infected with subtype B - 40.6% were infected by Portuguese. Among migrants from PALOPs - mainly infected with subtypes G and CFR02_AG - the transmission occurred predominantly within the migrants' community (53 and 80%, respectively).

CONCLUSION

The acquisition of infection among immigrants living in Portugal differs according to the country of origin. These results can contribute to monitor the HIV epidemic and prevent new HIV infections among migrants.

HIV-1-Transmitted Drug Resistance and Transmission Clusters in Newly Diagnosed Patients in Portugal Between 2014 and 2019

Objective

To describe and analyze transmitted drug resistance (TDR) between 2014 and 2019 in newly infected patients with HIV-1 in Portugal and to characterize its transmission networks.

Methods

Clinical, socioepidemiological, and risk behavior data were collected from 820 newly diagnosed patients in Portugal between September 2014 and December 2019. The sequences obtained from drug resistance testing were used for subtyping, TDR determination, and transmission cluster (TC) analyses.

Results

In Portugal, the overall prevalence of TDR between 2014 and 2019 was 11.0%. TDR presented a decreasing trend from 16.7% in 2014 to 9.2% in 2016 (p _for-trend = 0.114). Multivariate analysis indicated that TDR was significantly associated with transmission route (MSM presented a lower probability of presenting TDR when compared to heterosexual contact) and with subtype (subtype C presented significantly more TDR when compared to subtype B). TC analysis corroborated that the heterosexual risk group presented a higher proportion of TDR in TCs when compared to MSMs. Among subtype A1, TDR reached 16.6% in heterosexuals, followed by 14.2% in patients infected with subtype B and 9.4% in patients infected with subtype G.

Conclusion

Our molecular epidemiology approach indicates that the HIV-1 epidemic in Portugal is changing among risk group populations, with heterosexuals showing increasing levels of HIV-1 transmission and TDR. Prevention measures for this subpopulation should be reinforced.

A survey on the algorithm and development of multiple sequence alignment

Multiple sequence alignment (MSA) is an essential cornerstone in bioinformatics, which can reveal the potential information in biological sequences, such as function, evolution and structure. MSA is widely used in many bioinformatics scenarios, such as phylogenetic analysis, protein analysis and genomic analysis. However, MSA faces new challenges with the gradual increase in sequence scale and the increasing demand for alignment accuracy. Therefore, developing an efficient and accurate strategy for MSA has become one of the research hotspots in bioinformatics. In this work, we mainly summarize the algorithms for MSA and its applications in bioinformatics. To provide a structured and clear perspective, we systematically introduce MSA's knowledge, including background, database, metric and benchmark. Besides, we list the most common applications of MSA in the field of bioinformatics, including database searching, phylogenetic analysis, genomic analysis, metagenomic analysis and protein analysis. Furthermore, we categorize and analyze classical and state-of-the-art algorithms, divided into progressive alignment, iterative algorithm, heuristics, machine learning and divide-and-conquer. Moreover, we also discuss the challenges and opportunities of MSA in bioinformatics. Our work provides a comprehensive survey of MSA applications and their relevant algorithms. It could bring valuable insights for researchers to contribute their knowledge to MSA and relevant studies.

Evolution of Viral Pathogens Follows a Linear Order

Although lessons have been learned from previous severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) outbreaks, the rapid evolution of the viruses means that future outbreaks of a much larger scale are possible, as shown by the current coronavirus disease 2019 (COVID-19) outbreak. Therefore, it is necessary to better understand the evolution of coronaviruses as well as viruses in general. This study reports a comparative analysis of the amino acid usage within several key viral families and genera that are prone to triggering outbreaks, including coronavirus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2], SARS-CoV, MERS-CoV, human coronavirus-HKU1 [HCoV-HKU1], HCoV-OC43, HCoV-NL63, and HCoV-229E), influenza A (H1N1 and H3N2), flavivirus (dengue virus serotypes 1 to 4 and Zika) and ebolavirus (Zaire, Sudan, and Bundibugyo ebolavirus). Our analysis reveals that the distribution of amino acid usage in the viral genome is constrained to follow a linear order, and the distribution remains closely related to the viral species within the family or genus. This constraint can be adapted to predict viral mutations and future variants of concern. By studying previous SARS and MERS outbreaks, we have adapted this naturally occurring pattern to determine that although pangolin plays a role in the outbreak of COVID-19, it may not be the sole agent as an intermediate animal. In addition to this study, our findings contribute to the understanding of viral mutations for subsequent development of vaccines and toward developing a model to determine the source of the outbreak. IMPORTANCE This study reports a comparative analysis of amino acid usage within several key viral genera that are prone to triggering outbreaks. Interestingly, there is evidence that the amino acid usage within the viral genomes is not random but in a linear order.

Many but small HIV-1 non-B transmission chains in the Netherlands

OBJECTIVE

The aim of this study was to investigate introductions and spread of different HIV-1 subtypes in the Netherlands.

DESIGN

We identified distinct HIV-1 transmission chains in the Netherlands within the global epidemic context through viral phylogenetic analysis of partial HIV-1 polymerase sequences from individuals enrolled in the ATHENA national HIV cohort of all persons in care since 1996, and publicly available international background sequences.

METHODS

Viral lineages circulating in the Netherlands were identified through maximum parsimony phylogeographic analysis. The proportion of HIV-1 infections acquired in-country among heterosexuals and MSM was estimated from phylogenetically observed, national transmission chains using a branching process model that accounts for incomplete sampling.

RESULTS

As of 1 January 2019, 2589 (24%) of 10 971 (41%) HIV-1 sequenced individuals in ATHENA had non-B subtypes (A1, C, D, F, G) or circulating recombinant forms (CRF01AE, CRF02AG, CRF06-cpx). The 1588 heterosexuals were in 1224, and 536 MSM in 270 phylogenetically observed transmission chains. After adjustments for incomplete sampling, most heterosexual (75%) and MSM (76%) transmission chains were estimated to include only the individual introducing the virus (size = 1). Onward transmission occurred mostly in chains size 2-5 amongst heterosexuals (62%) and in chains size at least 10 amongst MSM (64%). Considering some chains originated in-country from other risk-groups, 40% (95% confidence interval: 36-44) of non-B-infected heterosexuals and 62% (95% confidence interval: 49-73) of MSM-acquired infection in-country.

CONCLUSION

Although most HIV-1 non-B introductions showed no or very little onward transmission, a considerable proportion of non-B infections amongst both heterosexuals and MSM in the Netherlands have been acquired in-country.

OUP accepted manuscript

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.

Multiple Occurrences of a 168-Nucleotide Deletion in SARS-CoV-2 ORF8, Unnoticed by Standard Amplicon Sequencing and Variant Calling Pipelines

Genomic surveillance of the SARS-CoV-2 pandemic is crucial and mainly achieved by amplicon sequencing protocols. Overlapping tiled-amplicons are generated to establish contiguous SARS-CoV-2 genome sequences, which enable the precise resolution of infection chains and outbreaks. We investigated a SARS-CoV-2 outbreak in a local hospital and used nanopore sequencing with a modified ARTIC protocol employing 1200 bp long amplicons. We detected a long deletion of 168 nucleotides in the ORF8 gene in 76 samples from the hospital outbreak. This deletion is difficult to identify with the classical amplicon sequencing procedures since it removes two amplicon primer-binding sites. We analyzed public SARS-CoV-2 sequences and sequencing read data from ENA and identified the same deletion in over 100 genomes belonging to different lineages of SARS-CoV-2, pointing to a mutation hotspot or to positive selection. In almost all cases, the deletion was not represented in the virus genome sequence after consensus building. Additionally, further database searches point to other deletions in the ORF8 coding region that have never been reported by the standard data analysis pipelines. These findings and the fact that ORF8 is especially prone to deletions, make a clear case for the urgent necessity of public availability of the raw data for this and other large deletions that might change the physiology of the virus towards endemism.

Limited HIV-1 Subtype C nef 3'PPT Variation in Combination Antiretroviral Therapy Naïve and Experienced People Living with HIV in Botswana

Dolutegravir (DTG) is a potent anti-HIV drug that is used to treat HIV globally. There have been reports of mutations in the HIV-1 3'-polypurine tract (3'PPT) of the nef gene, contributing to DTG failure; however, there are limited 'real-world' data on this. In addition, there is a knowledge gap on the variability of 3'PPT residues in patients receiving combination antiretroviral therapy (cART) with and without viral load (VL) suppression. HIV-1 subtype C (HIV-1C) whole-genome sequences from cART naïve and experienced individuals were generated using next-generation sequencing. The nef gene sequences were trimmed from the generated whole-genome sequences using standard bioinformatics tools. In addition, we generated separate integrase and nef gene sequences by Sanger sequencing of plasma samples from individuals with virologic failure (VF) while on a DTG/raltegravir (RAL)-based cART. Analysis of 3'PPT residues was performed, and comparison of proportions computed using Pearson's chi-square test with p-values < 0.05 was considered statistically significant. A total of 6009 HIV-1C full genome sequences were generated and had a median log10 HIV-1 VL (Q1, Q3) copies/mL of 1.60 (1.60, 2.60). A total of 12 matching integrase and nef gene sequences from therapy-experienced participants failing DTG/ RAL-based cART were generated. HIV-1C 3'PPT nef gene sequences from therapy-experienced patients failing DTG cART (n = 12), cART naïve individuals (n = 1263), and individuals on cART with and without virological suppression (n = 4696) all had a highly conserved 3'PPT motif with no statistically significant differences identified. Our study confirms the high conservation of the HIV-1 nef gene 3'PPT motif in 'real-world' patients and showed no differences in the motif according to VL suppression or INSTI-based cART failure. Future studies should explore other HIV-1 regions outside of the pol gene for associations with DTG failure.

SARS-CoV-2 variant evolution in the United States: High accumulation of viral mutations over time likely through serial Founder Events and mutational bursts

Since the first case of COVID-19 in December 2019 in Wuhan, China, SARS-CoV-2 has spread worldwide and within a year and a half has caused 3.56 million deaths globally. With dramatically increasing infection numbers, and the arrival of new variants with increased infectivity, tracking the evolution of its genome is crucial for effectively controlling the pandemic and informing vaccine platform development. Our study explores evolution of SARS-CoV-2 in a representative cohort of sequences covering the entire genome in the United States, through all of 2020 and early 2021. Strikingly, we detected many accumulating Single Nucleotide Variations (SNVs) encoding amino acid changes in the SARS-CoV-2 genome, with a pattern indicative of RNA editing enzymes as major mutators of SARS-CoV-2 genomes. We report three major variants through October of 2020. These revealed 14 key mutations that were found in various combinations among 14 distinct predominant signatures. These signatures likely represent evolutionary lineages of SARS-CoV-2 in the U.S. and reveal clues to its evolution such as a mutational burst in the summer of 2020 likely leading to a homegrown new variant, and a trend towards higher mutational load among viral isolates, but with occasional mutation loss. The last quartile of 2020 revealed a concerning accumulation of mostly novel low frequency replacement mutations in the Spike protein, and a hypermutable glutamine residue near the putative furin cleavage site. Finally, end of the year data and 2021 revealed the gradual increase to prevalence of known variants of concern, particularly B.1.1.7, that have acquired additional Spike mutations. Overall, our results suggest that predominant viral genomes are dynamically evolving over time, with periods of mutational bursts and unabated mutation accumulation. This high level of existing variation, even at low frequencies and especially in the Spike-encoding region may become problematic when super-spreader events, akin to serial Founder Events in evolution, drive these rare mutations to prominence.

A Novel SARS-CoV-2 Viral Sequence Bioinformatic Pipeline Has Found Genetic Evidence That the Viral 3' Untranslated Region (UTR) Is Evolving and Generating Increased Viral Diversity

An unprecedented amount of SARS-CoV-2 sequencing has been performed, however, novel bioinformatic tools to cope with and process these large datasets is needed. Here, we have devised a bioinformatic pipeline that inputs SARS-CoV-2 genome sequencing in FASTA/FASTQ format and outputs a single Variant Calling Format file that can be processed to obtain variant annotations and perform downstream population genetic testing. As proof of concept, we have analyzed over 229,000 SARS-CoV-2 viral sequences up until November 30, 2020. We have identified over 39,000 variants worldwide with increased polymorphisms, spanning the ORF3a gene as well as the 3' untranslated (UTR) regions, specifically in the conserved stem loop region of SARS-CoV-2 which is accumulating greater observed viral diversity relative to chance variation. Our analysis pipeline has also discovered the existence of SARS-CoV-2 hypermutation with low frequency (less than in 2% of genomes) likely arising through host immune responses and not due to sequencing errors. Among annotated non-sense variants with a population frequency over 1%, recurrent inactivation of the ORF8 gene was found. This was found to be present in the newly identified B.1.1.7 SARS-CoV-2 lineage that originated in the United Kingdom. Almost all VOC-containing genomes possess one stop codon in ORF8 gene (Q27∗), however, 13% of these genomes also contains another stop codon (K68∗), suggesting that ORF8 loss does not interfere with SARS-CoV-2 spread and may play a role in its increased virulence. We have developed this computational pipeline to assist researchers in the rapid analysis and characterization of SARS-CoV-2 variation.

ViralMSA: massively scalable reference-guided multiple sequence alignment of viral genomes

MOTIVATION

In molecular epidemiology, the identification of clusters of transmissions typically requires the alignment of viral genomic sequence data. However, existing methods of multiple sequence alignment (MSA) scale poorly with respect to the number of sequences.

RESULTS

ViralMSA is a user-friendly reference-guided MSA tool that leverages the algorithmic techniques of read mappers to enable the MSA of ultra-large viral genome datasets. It scales linearly with the number of sequences, and it is able to align tens of thousands of full viral genomes in seconds. However, alignments produced by ViralMSA omit insertions with respect to the reference genome.

AVAILABILITY AND IMPLEMENTATION

ViralMSA is freely available at https://github.com/niemasd/ViralMSA as an open-source software project.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computational strategies to combat COVID-19: useful tools to accelerate SARS-CoV-2 and coronavirus research

SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) is a novel virus of the family Coronaviridae. The virus causes the infectious disease COVID-19. The biology of coronaviruses has been studied for many years. However, bioinformatics tools designed explicitly for SARS-CoV-2 have only recently been developed as a rapid reaction to the need for fast detection, understanding and treatment of COVID-19. To control the ongoing COVID-19 pandemic, it is of utmost importance to get insight into the evolution and pathogenesis of the virus. In this review, we cover bioinformatics workflows and tools for the routine detection of SARS-CoV-2 infection, the reliable analysis of sequencing data, the tracking of the COVID-19 pandemic and evaluation of containment measures, the study of coronavirus evolution, the discovery of potential drug targets and development of therapeutic strategies. For each tool, we briefly describe its use case and how it advances research specifically for SARS-CoV-2. All tools are free to use and available online, either through web applications or public code repositories. Contact:evbc@unj-jena.de.

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent "paralog" gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.IMPORTANCE Orf8 is one of the most puzzling genes in the SARS lineage of coronaviruses, including SARS-CoV-2. Using sophisticated sequence comparisons, we confirm its origins from Orf7a, another gene in the lineage that appears as more conserved, compared to Orf8. Orf7a is a potential immune antagonist of known structure, while a deletion of Orf8 was shown to decrease the severity of the infection in a cohort study. The subtle sequence similarities imply that Orf8 has the same immunoglobulin-like fold as Orf7a, confirmed by structure determination. We characterize the subgroups of this superfamily and demonstrate the highly idiosyncratic divergence patterns during the evolution of the virus.

Isolating SARS-CoV-2 Strains From Countries in the Same Meridian: Genome Evolutionary Analysis

BACKGROUND

COVID-19, caused by the novel SARS-CoV-2, is considered the most threatening respiratory infection in the world, with over 40 million people infected and over 0.934 million related deaths reported worldwide. It is speculated that epidemiological and clinical features of COVID-19 may differ across countries or continents. Genomic comparison of 48,635 SARS-CoV-2 genomes has shown that the average number of mutations per sample was 7.23, and most SARS-CoV-2 strains belong to one of 3 clades characterized by geographic and genomic specificity: Europe, Asia, and North America.

OBJECTIVE

The aim of this study was to compare the genomes of SARS-CoV-2 strains isolated from Italy, Sweden, and Congo, that is, 3 different countries in the same meridian (longitude) but with different climate conditions, and from Brazil (as an outgroup country), to analyze similarities or differences in patterns of possible evolutionary pressure signatures in their genomes.

METHODS

We obtained data from the Global Initiative on Sharing All Influenza Data repository by sampling all genomes available on that date. Using HyPhy, we achieved the recombination analysis by genetic algorithm recombination detection method, trimming, removal of the stop codons, and phylogenetic tree and mixed effects model of evolution analyses. We also performed secondary structure prediction analysis for both sequences (mutated and wild-type) and "disorder" and "transmembrane" analyses of the protein. We analyzed both protein structures with an ab initio approach to predict their ontologies and 3D structures.

RESULTS

Evolutionary analysis revealed that codon 9628 is under episodic selective pressure for all SARS-CoV-2 strains isolated from the 4 countries, suggesting it is a key site for virus evolution. Codon 9628 encodes the P0DTD3 (Y14_SARS2) uncharacterized protein 14. Further investigation showed that the codon mutation was responsible for helical modification in the secondary structure. The codon was positioned in the more ordered region of the gene (41-59) and near to the area acting as the transmembrane (54-67), suggesting its involvement in the attachment phase of the virus. The predicted protein structures of both wild-type and mutated P0DTD3 confirmed the importance of the codon to define the protein structure. Moreover, ontological analysis of the protein emphasized that the mutation enhances the binding probability.

CONCLUSIONS

Our results suggest that RNA secondary structure may be affected and, consequently, the protein product changes T (threonine) to G (glycine) in position 50 of the protein. This position is located close to the predicted transmembrane region. Mutation analysis revealed that the change from G (glycine) to D (aspartic acid) may confer a new function to the protein-binding activity, which in turn may be responsible for attaching the virus to human eukaryotic cells. These findings can help design in vitro experiments and possibly facilitate a vaccine design and successful antiviral strategies.

Field-Adapted Full Genome Sequencing of Peste-Des-Petits-Ruminants Virus Using Nanopore Sequencing

Peste-des-petits-ruminants virus (PPRV) is currently the focus of a control and eradication program. Full genome sequencing has the opportunity to become a powerful tool in the eradication program by improving molecular epidemiology and the study of viral evolution. PPRV is prevalent in many resource-constrained areas, with long distances to laboratory facilities, which can lack the correct equipment for high-throughput sequencing. Here we present a protocol for near full or full genome sequencing of PPRV. The use of a portable miniPCR and MinION brings the laboratory to the field and in addition makes the production of a full genome possible within 24 h of sampling. The protocol has been successfully used on virus isolates from cell cultures and field isolates from tissue samples of naturally infected goats.

RNA genome conservation and secondary structure in SARS-CoV-2 and SARS-related viruses: a first look

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Here we present a first look at RNA sequence conservation and structural propensities in the SARS-CoV-2 genome. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nt as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences subsequently reported from the COVID-19 outbreak, and we present a curated list of 30 "SARS-related-conserved" regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 "SARS-CoV-2-conserved-structured" regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the extended 5' UTR, frameshifting stimulation element, and 3' UTR. Lastly, we predict regions of the SARS-CoV-2 viral genome that have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 "SARS-CoV-2-conserved-unstructured" genomic regions may be most easily accessible by hybridization in primer-based diagnostic strategies.

Mat_peptide: comprehensive annotation of mature peptides from polyproteins in five virus families

MOTIVATION

Sequence repositories have few well-annotated virus mature peptide sequences. Therefore post-translational proteolytic processing of polyproteins into mature peptides (MPs) has been performed in silico, with a new computational method, for over 200 species in 5 pathogenic virus families (Caliciviridae, Coronaviridae, Flaviviridae, Picornaviridae and Togaviridae).

RESULTS

Using pairwise alignment with reference sequences, MPs have been annotated and their sequences made available for search, analysis and download. At publication the method had produced 156 216 sequences, a large portion of the protein sequences now available in https://www.viprbrc.org. It represents a new and comprehensive mature peptide collection.

AVAILABILITY AND IMPLEMENTATION

The data are available at the Virus Pathogen Resource https://www.viprbrc.org, and the software at https://github.com/VirusBRC/vipr_mat_peptide.

In Search of Covariates of HIV-1 Subtype B Spread in the United States-A Cautionary Tale of Large-Scale Bayesian Phylogeography

Infections with HIV-1 group M subtype B viruses account for the majority of the HIV epidemic in the Western world. Phylogeographic studies have placed the introduction of subtype B in the United States in New York around 1970, where it grew into a major source of spread. Currently, it is estimated that over one million people are living with HIV in the US and that most are infected with subtype B variants. Here, we aim to identify the drivers of HIV-1 subtype B dispersal in the United States by analyzing a collection of 23,588 pol sequences, collected for drug resistance testing from 45 states during 2004-2011. To this end, we introduce a workflow to reduce this large collection of data to more computationally-manageable sample sizes and apply the BEAST framework to test which covariates associate with the spread of HIV-1 across state borders. Our results show that we are able to consistently identify certain predictors of spread under reasonable run times across datasets of up to 10,000 sequences. However, the general lack of phylogenetic structure and the high uncertainty associated with HIV trees make it difficult to interpret the epidemiological relevance of the drivers of spread we are able to identify. While the workflow we present here could be applied to other virus datasets of a similar scale, the characteristic star-like shape of HIV-1 phylogenies poses a serious obstacle to reconstructing a detailed evolutionary and spatial history for HIV-1 subtype B in the US.

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.

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.

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.

Increasing importance of European lineages in seeding the hepatitis C virus subtype 1a epidemic in Spain

BackgroundReducing the burden of the hepatitis C virus (HCV) requires large-scale deployment of intervention programmes, which can be informed by the dynamic pattern of HCV spread. In Spain, ongoing transmission of HCV is mostly fuelled by people who inject drugs (PWID) infected with subtype 1a (HCV1a).AimOur aim was to map how infections spread within and between populations, which could help formulate more effective intervention programmes to halt the HCV1a epidemic in Spain.MethodsEpidemiological links between HCV1a viruses from a convenience sample of 283 patients in Spain, mostly PWID, collected between 2014 and 2016, and 1,317, 1,291 and 1,009 samples collected abroad between 1989 and 2016 were reconstructed using sequences covering the NS3, NS5A and NS5B genes. To efficiently do so, fast maximum likelihood-based tree estimation was coupled to a flexible Bayesian discrete phylogeographic inference method.ResultsThe transmission network structure of the Spanish HCV1a epidemic was shaped by continuous seeding of HCV1a into Spain, almost exclusively from North America and European countries. The latter became increasingly relevant and have dominated in recent times. Export from Spain to other countries in Europe was also strongly supported, although Spain was a net sink for European HCV1a lineages. Spatial reconstructions showed that the epidemic in Spain is diffuse, without large, dominant within-country networks.ConclusionTo boost the effectiveness of local intervention efforts, concerted supra-national strategies to control HCV1a transmission are needed, with a strong focus on the most important drivers of ongoing transmission, i.e. PWID and other high-risk populations.

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.