Next Article in Journal
The Optimized γ-Globin Lentiviral Vector GGHI-mB-3D Leads to Nearly Therapeutic HbF Levels In Vitro in CD34+ Cells from Sickle Cell Disease Patients
Next Article in Special Issue
Correction: Martínez-Murcia et al. SARS-CoV-2 Variants Identification; A Fast and Affordable Strategy Based on Partial S-Gene Targeted PCR Sequencing. Viruses 2022, 14, 2588
Previous Article in Journal
High-Risk Human Papillomavirus Detection via Cobas® 4800 and REBA HPV-ID® Assays
Previous Article in Special Issue
Structural Insight into the Resistance of the SARS-CoV-2 Omicron BA.4 and BA.5 Variants to Cilgavimab
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes

1
Department of Microbiology-Immunology, Feinberg School of Medicine, Northwestern University, Chicago, IL 60611, USA
2
Department of Medicine, Division of Infectious Diseases, Feinberg School of Medicine, Northwestern University, Chicago, IL 60611, USA
3
Center for Pathogen Genomics and Microbial Evolution, Havey Institute for Global Health, Feinberg School of Medicine, Northwestern University, Chicago, IL 60611, USA
4
Biosciences Division, School of Medicine, University of California Riverside, Riverside, CA 92507, USA
5
Computational Biology Division, Department of Integrative Biomedical Sciences, University of Cape Town, Cape Town 7925, South Africa
*
Authors to whom correspondence should be addressed.
Viruses 2022, 14(12), 2714; https://doi.org/10.3390/v14122714
Submission received: 15 October 2022 / Revised: 14 November 2022 / Accepted: 21 November 2022 / Published: 5 December 2022
(This article belongs to the Special Issue SARS-CoV-2 Genomics)

Abstract

:
The unprecedented growth of publicly available SARS-CoV-2 genome sequence data has increased the demand for effective and accessible SARS-CoV-2 data analysis and visualization tools. The majority of the currently available tools either require computational expertise to deploy them or limit user input to preselected subsets of SARS-CoV-2 genomes. To address these limitations, we developed ViralVar, a publicly available, point-and-click webtool that gives users the freedom to investigate and visualize user-selected subsets of SARS-CoV-2 genomes obtained from the GISAID public database. ViralVar has two primary features that enable: (1) the visualization of the spatiotemporal dynamics of SARS-CoV-2 lineages and (2) a structural/functional analysis of genomic mutations. As proof-of-principle, ViralVar was used to explore the evolution of the SARS-CoV-2 pandemic in the USA in pediatric, adult, and elderly populations (n > 1.7 million genomes). Whereas the spatiotemporal dynamics of the variants did not differ between these age groups, several USA-specific sublineages arose relative to the rest of the world. Our development and utilization of ViralVar to provide insights on the evolution of SARS-CoV-2 in the USA demonstrates the importance of developing accessible tools to facilitate and accelerate the large-scale surveillance of circulating pathogens.

1. Introduction

Since the onset of the coronavirus disease 2019 (COVID-19) pandemic, the continued mutation and diversification of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in the repeated emergence of new “variants of concern” (VOCs) with increased infectivity, transmissibility, and/or immune evasion properties [1,2,3,4,5]. Each VOC has been defined by a distinct set of protein mutations (missense or nonsynonymous substitutions, in-frame insertions, and deletions) that confer unique functional properties [1,6,7,8,9,10,11,12]. For example, the Alpha (B.1.1.7*) VOC was defined using a set of nine Spike mutations (N501Y, A570D, D614G, P681H, T716I, S982A, D1118H, 69-70Δ, 144Δ) that increased infectivity [13], transmissibility [14], and resistance to monoclonal antibody therapeutics [15]. Especially within the Spike open reading frame, a greater proportion of missense compared to synonymous mutations is indicative of strong positive selection for Spike proteins with altered structure and function [16,17]. Continued SARS-CoV-2 genomic surveillance is essential to identify new emergent variants with novel phenotypic properties that may alter best practices in public health and clinical care.
The remarkable global scientific response to the COVID-19 pandemic has led to the generation of vast amounts of publicly available SARS-CoV-2 whole-genome sequence data. Worldwide, most genome sequences are deposited in the GISAID public database (gisaid.org) [18], and more than 13 million viral sequences from around the world have been deposited as of 12 September 2022. This massive and ongoing SARS-CoV-2 sequencing effort has provided a unique opportunity to study the virus’s evolution in exquisite detail. However, at the same time, the volume and diversity of available sequences exacerbates the complexity of the data analysis and calls for effective tools to allow researchers with little or no computational expertise to perform detailed analyses of relevant genomic data.
In part to address this problem, several web-based tools have been developed to facilitate the study of SARS-CoV-2 spatiotemporal dynamics, mutational frequency, and/or three-dimensional (3D) protein structures [19,20,21,22]. Though useful for gaining broad insights, these applications are often limited to the analysis of predetermined datasets with minimal user control, such as COVIDCG [23], outbreak info [24], covariants [25], 2019nCoVR [26], CoV-GLUE [27], and COG-UK [28]. However, even tools that allow the processing of user-defined data often accept a limited number of sequencing data such as covdb (limit = 100) [29], coronApp [30] (limit ~100 MB or ~3500), and VirusViz (limit = 50) [31]. In addition to the lack of options for large-scale data analysis, these tools have limited analytical features for the multilevel analysis and visualization of SARS-CoV-2 lineages and their mutations (e.g., spatiotemporal visualization of lineages, linear or 3D visualizations of mutations in the context of proteins and genomes).
Other tools and databases have been developed to study SARS-CoV-2 protein structures. One of these applications is SARS-CoV-2 3D [32] which provides tools for 3D structure predictions and energy calculations to evaluate targets and design new potential therapeutics. CoV3D [33] is a repository for 3D protein structures of SARS-CoV-2 and host antibodies. Neither tool provides information on mutational changes in the context of the 3D structures. Other webservers such as the GISAID [18], covariants [25], and COG-UK [28] provide limited 3D structural visualizations for only fixed sets of mutations (mostly clade-defining) and only for the Spike protein. To the best of our knowledge, there are currently two webservers that enable the visualization of mutations in the context of 3D protein structures for all SARS-CoV-2 proteins: Coronavirus3D [19] and COVID-3D [21]. However, both servers have a fixed list of SARS-CoV-2 genomes/mutations in their databases and lack the option to visualize mutations based on user-provided genomes. Although COVID-3D [21] provides an option for the user to input variants, this is limited to only one protein at a time and thus requires multiple file uploads to visualize genomic mutations in the context of more than one protein. Moreover, both servers lack the capacity to examine mutational patterns at selected time points in the pandemic, in selected geographical regions, and/or among specific lineages.
Here, we present a new, web-based software application “ViralVar” that incorporates user-selected genome data to visualize and study lineages over time by depicting the distribution of mutations at both the nucleotide and protein levels as well as providing the context of variants in the 3D structure of SARS-CoV-2 proteins. Protein visualizations provide detailed information on the functional protein domains and predicted B-cell epitopes. Additionally, ViralVar provides a currently unique feature among similar applications that allows for the binomial testing of protein mutations to identify potential over- and under-mutated proteins, k-means clustering of genomes based on protein mutations to expedite large-scale surveillance of new mutations, and the visualization of changes in the mutational patterns of the virus over selected date ranges, within defined geographical regions, and/or within or among lineages. A practical demonstration of the application of ViralVar is given here by examining the relative dynamics of the evolution of SARS-CoV-2 in the USA, using a total of 1,739,797 sequenced genomes collected in the USA between January 2020 and May 2022. The ViralVar webserver is freely available at http://viralvar.org/.

2. Methods

2.1. General Software Workflow

The ViralVar webtool is implemented in the R programming language using Shiny, an open-source R package for developing interactive web applications. Shiny implements layout features available in Bootstrap, an HTML 4.01/shiny-css 1.7.1/shiny-javascript 1.7.1 framework. To add more advanced content to ViralVar, the user interface was customized with HTML and Shiny’s HTML tag attributes, as well as custom cascading style sheets (CSS) and other R packages listed in the context of the relevant sections below. Briefly, SARS-CoV-2 genomic data retrieved from the GISAID [18] were used as input for ViralVar. The webapp is divided into two modules (Figure 1). In the first module, “Lineage Dynamics”, data are processed to depict the spatiotemporal dynamics of SARS-CoV-2 lineages and clades in the form of stacked bars, area plots, and pie charts. The second module “Mutational Analysis” visualizes mutation distributions along the SARS-CoV-2 genome and proteins (linear and 3D) and generates statistical analyses to identify over- and under-mutated proteins. Users can interact with the server to explore and compare the temporal dynamics of the lineages and mutations between the different sets of genomes and/or VOCs. Each module provides various control options, allowing users to customize analyses and view and export figures according to their requirements. The output files of ViralVar are either high-resolution figures (PDF, PNG) or tables (tsv).

2.2. Data Input

The GISAID is one of the largest global resources for sharing SARS-CoV-2 genome sequences and associated clinical and demographic information [18]. GISAID data are accessible to users through free registration via the GISAID website (https://gisaid.org/) accessed on 31 May 2022. The database provides genome consensus sequences, reference-based multiple sequence genome alignments, and lists of mutations for each genome with the associated lineage or clade designations in tabular format. Data to be downloaded from the database can be readily filtered to focus on dates of collection, specific geographical regions, or selected lineages or clades. ViralVar accepts input data from the GISAID in a tabular format that includes three sets of information for each genome: (1) PANGO lineage (users can opt to manually add Nextclade designations), (2) sample collection date, and (3) a list of protein mutations (denoted as “AA Substitutions” in the GISAID data; required for “Mutational Analysis”). Detailed guidance on retrieving GISAID data in the proper format for input into ViralVar is provided on the “Home” tab of ViralVar. A limited set of 3892 SARS-CoV-2 sequences collected through the Northwestern Medicine Healthcare (NMH) system in Chicago, IL, between February 2020 and May 2022, are included in the ViralVar webtool for example purposes and can be viewed by checking the “Visualize Example Data” checkbox in each module. The GISAID IDs are provided in Table S1.

2.3. Lineage Dynamics

The “Lineage Dynamics” module of ViralVar serves to provide tools for visualizing changing trends in SARS-CoV-2 lineages/clades over time using temporal abundances and geographical distributions. ViralVar uses the R package ggplot2 [34] to generate visualizations reflecting the trend of changes in the absolute and relative abundances of SARS-CoV-2 lineages over time in the input data set. After data input, the data are displayed in a tabular format in the “Data Overview” tab. Note that for this module, only collection date and PANGO lineage information are required. The “Area Charts” and “Bar Charts” tabs illustrate the dynamics of lineage distributions over user-specified date ranges. The “Geographical map” tab shows lineage distributions overlaid as pie charts on user-selected geographical maps for the world, the USA, or individual USA states and territories again over a user-specified date range. Geographical maps are drawn using the R package maps and scatterpie. The phylogenetic nomenclature option allows users to customize output data to use PANGO lineage, Nextclade clade, or World Health Organization-defined VOC nomenclature. Tables and customizable figures are downloadable in a portable document format (PDF).

2.4. Mutational Analysis

The “Mutational Analysis” module of ViralVar provides users with a suite of tools to visualize the genomic and structural context of SARS-CoV-2 mutations. The R package ggplot2 [34] is used to generate and annotate density plots. After data input, the data are displayed in tabular format in the “Data Overview” tab. Note that for this module, collection date, PANGO lineage, and amino acid (AA) substitution information are required. The “Genome Distribution” tab depicts mutation density among uploaded sequences across the SARS-CoV-2 genome. Briefly, the number of distinct mutation events at each genomic position or protein residue is determined relative to a reference sequence (NCBI: NC_045512.2) [35] and reported over a sliding 100 nucleotide window. Position counts are calculated separately for insertions, deletions, and substitutions. This method does not consider virus counts in its calculation (i.e., the number of uploaded genomes with a particular mutation) such that each mutation event is counted only once. This avoids potential biases in reporting mutational frequency due to unequal amplification or sequencing across the genome as well as bias sampling [16]. In the “Protein Distribution” tab, the frequencies of genomes (virus counts) with mutations at specific protein residues are visualized using the R package ggplot2 [34] and plotly R package (interactive visualization). Separate plots can be generated for all SARS-CoV-2 proteins, both structural and nonstructural. Protein domain boundaries are indicated as described in the literature [16,17]. The IEDB server (Bepipred Linear Epitope Prediction 2.0 at http://www.iedb.org/) [36] (accessed on 31 October 2021) was used to predict B-cell epitopes, which are indicated above the protein schematic. In the “3D Protein Structure” tab, the R library package r3dmol is used to visualize mutations in the context of 3D protein structures. The 3D coordinates were obtained from the Protein Data Bank (PDB) with PDB accession numbers provided for each structure [37]. For proteins with no available 3D structure, models, as predicted by AlphaFold, were used when available [38]. Alternatively, the positions of transmembrane helices for proteins with no available 3D structures were identified using the TMHMM 2.0 algorithm [39]. Lists of the top mutations along with their frequencies for each protein can be downloaded in the form of tab delimited tables. The 3D protein illustrations can be downloaded as portable network graphics (PNG) files. Each of the above tabs includes a date slider to allow users to restrict data to a specific date range and a “Select VOC/VOI” option to limit output to a specified VOC or VOI.
The above mutational analysis tabs are further complemented by two tabs for statistical analysis and k-means clustering. In the “Statistical Analysis” tab, ViralVar utilizes the binomial test to identify individual proteins within the uploaded dataset that have significantly different mutation frequencies. The method has been previously applied to identify significantly under- and over-mutated SARS-CoV-2 proteins [16,17]. Briefly, the arguments for the binomial test are the observed number of distinct protein mutations in a certain protein (the “number of successes”), the total number of distinct protein mutations in all SARS-CoV-2 proteins (the “number of trials”), and the length of a given protein divided by the length of all SARS-CoV-2 proteins (the “expected probability of success”). An example of binomial calculations is provided below. For more details, please refer to [16].
To simplify the calculations in this method, we hypothesize that each protein mutation is an independent event and that all SARS-CoV-2 proteins and all residues have the same probability of being mutated. Therefore, this method applies the binomial test to assess the null hypothesis: that protein mutations are distributed randomly across all SARS-CoV-2 proteins.
P ( M P ,   M T ) = ( M T M P ) P ( p ) M P ( 1 ( P ( p ) ) M T M P
MT = the total number of protein mutations observed for all proteins (for example, 325 mutations in user input data)
MP = the number of protein mutations in the target protein (for example, 66 mutations in Spike in user input data)
P(p) = length protein/length proteome (e.g., length Spike/total length = 1273/9930 = ~0.13)
P ( M P ,   M T ) = ( 325 66 ) 0.13 66 ( 1 ( 0.13 ) 325 66 = 0.00046
Based on the null hypothesis, we expect only 42 mutations in Spike, given that 66 out of the total 325 mutations identified in SARS-CoV-2 proteins are located in Spike, the length of the Spike protein is of 1273 amino acids, and the entire SARS-CoV-2 proteome is 9930 long. However, the binomial test p-value (0.00046) suggests rejection of the null hypothesis and indicates a significantly higher number of mutations in the Spike proteins compared to the background (entire proteome). ViralVar conducts the above calculation for user input data; therefore, MP, MT, and P(p) will be different for each input dataset. An option to exclude clade signature mutations is provided to avoid bias in the binomial test across highly divergent clades. ViralVar also provides control options to customize binomial test parameters, including the option to adjust the p-values for multiple comparisons. As above, the tab includes a date slider to allow users to restrict data to a specific date range and a “Select VOC/VOI” option to limit output to a specified VOC or VOI. A results table of the analysis can be downloaded as a tsv file.
In the “Genome Clustering” tab, ViralVar employs k-means clustering to facilitate rapid investigation of emerging clusters of genomes with specific protein mutation. As the selection of mutations in SARS-CoV-2 evolution has been shown to be largely impacted by positive selection, driven by changes in SARS-CoV-2 protein structures and functions [16,17], targeting protein mutations could cluster genomes relative to the phenotype. For instance, a common feature of SARS-CoV-2 genomes with the N501Y spike mutation (e.g., Alpha, Beta and Gamma strains) was enhanced infectivity and transmissibility over the previous variants [14].
The clustering of genomes based on pairwise distance-based methods is computationally intensive and might take days to run depending on the computational resources. The runtime for the first step of these approaches (the calculation of distance matrices for all pairs of genomes) increases exponentially with the increase in the number of genomes (Figure S3). In contrast, k-means clustering of SARS-CoV-2 genomes has been proposed in the recent literature as a rapid method to investigate emerging variants and tackle the computational challenges in large-data analysis [40,41]. Due to its simplicity and being computationally inexpensive, the k-means clustering of genomes, based on mutations in specific proteins, can be quickly and repeatedly run on large-scale genomic datasets (such as ~11.1 M SARS-CoV-2 genomes).
ViralVar uses k-means to group genomes-based on protein mutations. To avoid the effects of spurious mutations (e.g., due to sequencing or assembly errors), the clustering of the genomes is calculated only from protein mutations with a default minimum mutation frequency (MMF) of >0.005, although this cutoff is user-adjustable. To determine the optimal number of clusters, ViralVar repeats k-means clustering for numbers of clusters (determined based on the number of variables in the input file) and calculates the average silhouette width (ASW) index using the R package NbClust [42]. In the calculation of the ASW, ViralVar uses unique genomes (duplicated genomes with identical mutational patterns are removed) to make calculations less computationally expensive. However, the final clustering is applied to all of the genomes in the input data to produce counts of the genomes in each cluster. As with the previously described functions, VOC/VOI and date range are selectable. The protein selection option allows for targeting mutations along a protein of interest. Tables and customizable figures in a PDF format are downloadable.

2.5. Applying ViralVar to Assess Dynamics of SARS-CoV-2 Evolution

The following was retrieved from the GISAID for the collection dates between 1 January 2020 and 15 May 2022 (downloaded 31 May 2022): a total of 1,739,797 SARS-CoV-2 high-quality complete genome assemblies (the GISAID criteria, including N content <5%) from the USA for which patient age, collection date (day/month/year), and geographic location available. To study the dynamics of the SARS-CoV-2 evolution using ViralVar, sequenced samples were classified into three populations by age: children (0–18 years), adults (18–64 years), and elderly (65+ years) (Table 1). The list of the GISAID identifiers that compose each group is provided in Table S1. Sequence data for each age group were uploaded separately and analyzed using ViralVar. Mutation distributions were also compared between the SARS-CoV-2 genomes collected and sequenced for different age groups in the USA. Additionally, to show the application of ViralVar to smaller sample sets (i.e., from regions with lower rates of genome sequencing relative to the USA), we targeted 248 SARS-CoV-2 genome sequences collected in Nigeria between 15 December 2020 and 15 January 2021 and 90 sequences collected from Cape Town, South Africa between 15 December 2020 and 15 January 2021. The list of the GISAID identifiers is provided in Table S1. Each of the Cape Town and Nigeria samples were processed using ViralVar.

3. Results and Discussion

3.1. Spatiotemporal Dynamics of SARS-CoV-2 VOCs in the USA

The United States has experienced one of the world’s highest COVID-19 burdens during the pandemic, with a total of 86.4 M confirmed cases and 1.01 M deaths as of 31 May 2022. Whereas some reports detailing the evolution of the COVID-19 pandemic in select cities and states are available [43,44], there are few comprehensive reports at a national level. To demonstrate the capabilities of ViralVar, we downloaded all high-quality whole genome sequence data available in the GISAID on specimens collected in the USA between 1 January 2020 and 15 May 2022 (n = 1,739,797 SARS-CoV-2 sequences total). These data were sorted by age (children, adults, elderly) and uploaded into the ViralVar webtool for analysis.
The temporal dynamics of the VOCs across age groups in the USA were visualized using the “Area Chart” tab in the “Lineage Dynamics” module of ViralVar (Figure 2A). The results indicate that the dynamics of the VOCs were relatively similar for all age groups. The SARS-CoV-2 lineage B.1.1.7, designated by the WHO as “Alpha”, was the first named VOC and likely emerged in the United Kingdom (UK) in September 2020. Alpha rapidly displaced other circulating lineages in the USA and became one of the top circulating VOCs in the world in early 2021 [45]. The emergence of Alpha in the USA can be tracked back to November 2020 (Figure 2B), coincident with a spike of new cases and deaths between November 2020 and March 2021 (Figure S1). Using the date range feature to focus on dynamics during these months, Alpha emerges as the dominant variant at the tail-end of the surge in cases, suggesting that it was not responsible for the rise in cases but rather took over after the contraction in cases of the previous variant (Figure 2B).
Utilizing the “Geographical Map” feature in the “Lineage Dynamics” module, the distribution of VOCs collected between January 2020 and May 2022 was visualized for each age group by state (Figure 3A). The lineage distributions were similar across states between all age groups, with the Omicron and Delta VOCs making up the majority of cases, followed by Alpha (Figure 3A and Figure S2). In narrower timeframes, however, distinct spatiotemporal trends become more obvious. Using the date control feature, we adjusted this analysis to examine cases between November 2020 and March 2021 (Figure 3B). Whereas Alpha lineages made up the majority of cases in most states over this time period, region-specific trends emerged. For example, the Epsilon VOI was responsible for a substantial number of cases in southwestern states, whereas the Iota VOI was more prevalent in the northeast. Illinois specifically reported a substantial number of cases of the Gamma VOC that were not reflected in the neighboring states. These region-specific trends were consistent across age groups (Figure 3B).

3.2. Mutational Analysis of Alpha Variant Sublineages in the USA

We subset the USA data explained earlier (n = 1,739,797 SARS-CoV-2 sequences) to only include genomes assigned to the Alpha lineage (n = 140,100). Additionally, genomes assigned to the Alpha lineage collected from the entire world (n = 906,114 excluding the USA) were retrieved from the GISAID as of 31 May 2022. Using the ViralVar “Mutational Analysis” module, the mutation profile for the Alpha VOC in the USA was compared to specimens from other countries. All ages were grouped together for this analysis due to the relatively small sample size of the under 18 and over 65 populations compared to adults in this dataset. Using the “Protein Distribution” tab, we visualized the mutational frequency in the Alpha VOC sequences at sites across Spike and NSP12 in both the USA and in the rest of the world (Figure 4). Whereas the defining mutations of the Alpha VOC were universally present, a distinct subset of mutations were more prevalent in the USA, specifically the Spike mutation K1191N and the NSP12 mutation P227L. To further investigate these mutations, all genomes containing Spike K1191N (n = 51,713) and NSP12 P227L (n = 190,869) mutations were retrieved from the GISAID and uploaded into ViralVar. A majority of genomes with the Spike K1191N mutation were in the Alpha variant genomes (80.4%, 41,558 of 51,713 genomes), of which the vast majority came from the USA (93.5%, 38,837 of 41,558 genomes) (Figure 5A, top). Similarly, 169,314 genomes with the NSP12 P227L mutation were classified as Alpha variants (88.7%), of which 104,435 genomes were collected in the USA (61.6%) (Figure 5A, bottom).
The Spike K1191N and NSP12 P227L appear to be recurrent mutations that have emerged in several other VOCs (i.e., Delta, Omicron, and Gamma); however, there is a lack of evidence regarding their role in virus infectivity, transmissibility, and/or clinical outcomes. To gain insight into their possible functional roles, we examined the protein context of each mutation using the “3D Protein Structure” feature in ViralVar (Figure 5B). The NSP12 P227L mutation is located in the Nidovirus RdRp-associated nucleotidyl transferase (NiRAN) domain. Although it is surface-exposed, it is far from the RNA binding or enzymatic active site. That being said, a nearby mutation in the NiRAN domain, N198S, has been recently reported as a potential antiviral resistance mutation to the NSP12-targeting drug, remdesivir [46]. Given the high level of conservation among coronavirus RNA-dependent RNA polymerases (RdRps) [47] and the recurring, but infrequent, prevalence of this mutation, it may also be that the P227L mutation confers some selective benefit but at a fitness cost to the virus [46]. The Spike mutation K1191N is located in the S2 subunit in the heptad repeat 2 (HR2) subdomain of the Spike protein, which is involved in host cell membrane fusion and viral entry (Figure 5B) [48]. Other Spike protein mutations in the HR2 subdomain, such as V1176F, have been shown to augment the stability of Spike and have been associated with increased disease severity and mortality [49,50,51]. More studies are required to determine the functional consequences of both Spike K1191N and NSP12 P227L.

3.3. ViralVar K-Means Clustering Feature Identifies Subclusters of the Alpha Variant in the USA

To better understand the genomic context of these mutations, we used the “Clustering Analysis” feature in ViralVar to identify co-occurring groups of mutations. K-means clustering based on Euclidean distance was applied to all of the Alpha VOC sequences collected in the USA, using a minimum mutation frequency cutoff of 0.005 and with a focus on the Spike and NSP12 proteins. The clustering of the Alpha genomes based on Spike mutations resulted in three distinct clusters (Figure 5C), two of which were defined by the presence (cluster 3) or absence (cluster 1) of K1191N. A third cluster showed a minor presence of K1191N but concurrently lacked the S982A and/or T716I mutations (cluster 2). The clustering of the Alpha genomes based on NSP12 mutations identified two distinct clusters distinguished solely by the P227L mutation (Figure 5C). To determine if these clusters can also be identified using phylogenetic analysis, we examined these mutations using the Nextstrain webserver (Figure 5D). The time-resolved phylogenetic trees from Nextstrain suggest that the Spike K1191N mutation is monophyletic, whereas the P227L mutation arose in at least two distinct branches (Figure 5D). The k-means clustering is largely in accordance with the phylogenetic analysis but suggests that additional mutational information, including synonymous mutations and those that occur outside of the open reading frame of interest, capture additional information not accounted for in this approach.
One of the limitations of the phylogenetic tree-based analysis, clustering, and visualization of SARS-CoV-2 genomes and investigating protein mutations is the computational cost that multiplies with the number of available genomes. The majority of studies using phylogenetic trees to study SARS-CoV-2 variants of concern (VOCs), therefore, must rely on subsampling approaches [52,53]. The k-means-based clustering of SARS-CoV-2 genomes based on Euclidean distance is one way to overcome this challenge as the method calculates the distance of each datapoint to the centroid using pairwise distances instead, decreasing the computational cost of analyzing the additional sequences (Figure S3). Furthermore, the k-means clustering of genomes based on protein mutations can be leveraged to the group genomes in a way directly related to the phenotype [54,55]. The congruence between the approach taken by ViralVar (Figure 5C) and the phylogenetic analysis results (Figure 5D) support the potential use of k-means clustering for the rapid analysis of large genomic datasets to facilitate tracking emerging protein mutations using a generic clustering method. This method could also be readily adapted and applied to other viruses. That being said, this approach is not suitable for making specific evolutionary inferences and so can be considered complementary to traditional phylogenetic-tree-based methods and useful for initial analyses and hypothesis generation.

3.4. Significant Nonrandom Distribution of Mutations in SARS-CoV-2 Proteins

To explore the different mutational profiles in genomes collected from different age groups in the USA, we used the “Genome Distribution” feature of ViralVar to visualize the mutations in all of the collected specimens from the USA split by age group (Figure 6). Overall, the analyses of the mutation profiles of the SARS-CoV-2 genomes were relatively similar for the three age groups in the USA samples (Figure 6). Compared to structural and accessory proteins, nonstructural proteins seem to undergo a higher mutational constraint (Figure 6 and Table S2), consistent with the previous reports [16]. The slight variability in the mutational patterns between different data subsets could be partly attributed to the differences in the population size and sampling dates between regions and age groups.
One of the most noteworthy differences when comparing the results from the first year of the pandemic [16] and results obtained in this study is the increased frequency of protein indel events, especially the accumulation of insertions in the Spike NTD. This trend was consistent for samples collected across all age groups, though distinct deletion events appeared more prevalent in elderly populations (for example, in the NSP15 open reading frame, Figure 6). The increased frequency of recurrent indels and their nonrandom distribution is believed to be an adaptive response mechanism to elevated global herd immunity, resulting from vaccination, infection, or both [17,56,57]. Spike NTD indels could alter neutralizing epitopes in the region and are thought to result in reduced antibody protection against VOCs that harbor these indels [56].
Using the “Statistical Analysis” feature of ViralVar, we further identified significant accumulations of mutations in mostly the structural proteins of SARS-CoV-2 with two exceptions for the nonstructural proteins (NSP1 and NSP2). Of note, a higher concentration of mutations was observed in NSP1 (average odds ratio = 1.46, q-value = 0 across all age groups), NSP2 (average odds ratio = 1.3, q-value = 0 across all age groups), N (average odds ratio = 1.6, q-value = 0 across all age groups), NS6 (average odds ratio = 1.6, q-value = 0 across all age groups), NS7a (average odds ratio = 3.1, q-value = 0 across all age groups), NS7b (average odds ratio = 1.8, q-value = 0 across all age groups), NS8 (average odds ratio = 3.1, q-value = 0 across all age groups), and Spike (average odds ratio = 1.4, q-value = 0 across all age groups) (Table S2). All these proteins are involved in interactions with the host immune system [58,59,60]. Recurrent NSP1 substitutions and indels have been found to accumulate on the protein surface and near epitope regions [17] and are thought to adversely affect the host’s immune response and vaccine efficiency [61,62]. For instance, NSP1 Δ79-89 induces a lower IFN-I response in the infected Calu-3 cells [62], highlighting the biological importance of mutations in NSP1 and other nonstructural proteins. The significantly higher concentration of mutations in the specific proteins involved in host immune interactions, the emergence of new types of protein mutations (in-frame indels), and the expansion of mutations to new proteins or protein regions suggest the virus is evolving to combat the host immune system. Taken together, nonrandom distribution of the mutations in different SARS-CoV-2 proteins suggests proteins undergo different evolutionary pressures driven partly by the host immune system.

3.5. ViralVar Potential in Identifying Novel Variants in Small and Local Cohorts

Using ViralVar, we explored the evolution of 90 SARS-CoV-2 genome sequences collected in Cape Town, South Africa between 1 October 2020 and 30 November 2020 and 248 sequences collected in Nigeria between 15 December 2020 and 15 January 2021. Our analysis suggests the presence of two clusters based on the Spike protein mutations in the Cape Town samples (Figure S4A). Cluster 1 samples were all assigned to the Beta VOC (B.1.351) (Figure S4B) which was first described in South Africa [1]. We identified three distinct clusters of genomes based on Spike mutations in the Nigeria sample cohort (Figure S4C). The second cluster in this analysis corresponds to the Eta variant (B.1.525) (Figure S4D) which was identified as a variant local to West Africa [63]. These analyses show the potential of ViralVar for the analysis and tracking of mutations in small and regionally collected datasets.

4. Conclusions

The emergence of new variants of SARS-CoV-2 with higher transmissibility and enhanced immune evasion highlights the need for ongoing SARS-CoV-2 genomic surveillance. This work has been greatly facilitated by public sequence repositories such as the GISAID which contained available data for more than 11.1 M genome sequences as of 31 May 2022. At the same time, this vast amount of genomic data has increased the demand for more flexible and multilevel analysis platforms to help study the virus evolution. To complement and expand upon previously developed analysis tools, we created ViralVar, a webtool for visualizing and researching SARS-CoV-2 lineages and mutational patterns over time. We have shown that ViralVar can be deployed as a point-and-click tool to rapidly investigate the spatiotemporal evolution of large numbers of SARS-CoV-2 genomes. Overall, our findings utilizing ViralVar offer important insights into pathogen evolution dynamics and spread in the USA. This study demonstrates that ViralVar can be successfully used to study the evolution of SARS-CoV-2 and help in improving global COVID-19 mitigation plans as the pandemic continues to evolve.
As part of a larger project for facilitating the study of virus evolution and mutational patterns, the development of ViralVar will continue for the study of other viruses. Additional future work includes the addition of multiple data input options (i.e., consensus sequences or multiple sequence alignments) to facilitate users in analyzing their own data. The continued enrichment of the list of the structural and functional properties of SARS-CoV-2 and other viral proteins in ViralVar will also take place on a regular basis. The ViralVar databases will be updated at regular intervals based upon information provided for other viruses and updates in public databases for protein structural and functional properties. ViralVar complements current tools for studying the massive number of SARS-CoV-2 genomes and can provide a user-friendly platform for the multilevel study of SARS-CoV-2 evolution.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14122714/s1, Figure S1: Frequency of new cases and deaths in the USA; Figure S2: Area plots reflecting the absolute abundance of variants of concern and variants of interest collected in the USA over time; Figure S3: Runtime performance for K-means clustering. Figure S4: ViralVar Potential in Identifying Novel Variants in Small and Local Cohorts; Table S1: List of GISAID SARS-CoV-2 genome sequence IDs used in this paper; Table S2. Significant over and under-mutated SARS-CoV-2 proteins.

Author Contributions

Conceptualization, A.A.; methodology, A.A., L.J., A.I., R.L.-R., E.A.O. and A.G.; software, A.A.; data curation, A.A., L.M.S. and T.J.D.; writing—original draft preparation, A.A. and E.A.O.; writing—review and editing, A.A., E.A.O., L.M.S., R.L.-R. and J.F.H.; visualization, A.A., L.J. and A.I.; supervision, E.A.O. and J.F.H.; project administration, E.A.O. and J.F.H.; funding acquisition, E.A.O. and J.F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by grants (R21 AI163912 to Hultquist, U19 AI171110 to Hultquist, U19 AI135964 to Ozer) from the National Institutes of Health; a grant (to Lorenzo-Redondo) from the Northwestern University Havey Institute for Global Health; and a grant (to Ozer and Hultquist) from the Walder Foundation’s Chicago Coronavirus Assessment Network. This research was supported, in part, through the computational resources and staff contributions provided by the Quest high performance computing facility at Northwestern University, which is jointly supported by the Office of the Provost, the Office for Research, and Northwestern University Information Technology. The funding sources had no role in the study design, data collection, analysis, interpretation, or writing of the report.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The ViralVar webserver is freely accessible through http://viralvar.org/.

Acknowledgments

We gratefully acknowledge the authors from the originating laboratories and the submitting laboratories, who generated and shared via the GISAID genetic sequence data on which this research is based, as well as the structural biology groups for contributing their structures to the PDB. This research is supported by grants (R21 AI163912 to Hultquist, U19 AI171110 to Hultquist, U19 AI135964 to Ozer) from the National Institutes of Health; a grant (to Lorenzo-Redondo) from the Northwestern University Havey Institute for Global Health; and a grant (to Ozer and Hultquist) from the Walder Foundation’s Chicago Coronavirus Assessment Network. This research was supported, in part, through the computational resources and staff contributions provided by the Quest high performance computing facility at Northwestern University, which is jointly supported by the Office of the Provost, the Office for Research, and Northwestern University Information Technology. The funding sources had no role in the study design, data collection, analysis, interpretation, or writing of the report.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tegally, H.; Wilkinson, E.; Giovanetti, M.; Iranzadeh, A.; Fonseca, V.; Giandhari, J.; Doolabh, D.; Pillay, S.; San, E.J.; Msomi, N.; et al. Detection of a SARS-CoV-2 variant of concern in South Africa. Nature 2021, 592, 438–443. [Google Scholar] [CrossRef] [PubMed]
  2. Karim, S.S.A.; Karim, Q.A. Omicron SARS-CoV-2 variant: A new chapter in the COVID-19 pandemic. Lancet 2021, 398, 2126–2128. [Google Scholar] [CrossRef] [PubMed]
  3. Viana, R.; Moyo, S.; Amoako, D.G.; Tegally, H.; Scheepers, C.; Lessells, R.J.; Giandhari, J.; Wolter, N.; Everatt, J.; Rambaut, A. Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa. Nature 2022, 603, 679–686. [Google Scholar] [CrossRef] [PubMed]
  4. Madhi, S.A.; Baillie, V.; Cutland, C.L.; Voysey, M.; Koen, A.L.; Fairlie, L.; Padayachee, S.D.; Dheda, K.; Barnabas, S.L.; Bhorat, Q.E.; et al. Efficacy of the ChAdOx1 nCoV-19 CoV-19 Vaccine against the B.1.351 Variant. N. Engl. J. Med. 2021, 384, 1885–1898. [Google Scholar] [CrossRef] [PubMed]
  5. Jewell, B.L. Monitoring differences between the SARS-CoV-2 B.1.1.7 variant and other lineages. Lancet Public Health 2021, 6, e267–e268. [Google Scholar] [CrossRef]
  6. Jassat, W.; Mudara, C.; Ozougwu, L.; Tempia, S.; Blumberg, L.; Davies, M.-A.; Pillay, Y.; Carter, T.; Morewane, R.; Wolmarans, M. Difference in mortality among individuals admitted to hospital with COVID-19 during the first and second waves in South Africa: A cohort study. Lancet Glob. Health 2021, 9, e1216–e1225. [Google Scholar] [CrossRef]
  7. Edward, P.R.; Lorenzo-Redondo, R.; Reyna, M.E.; Simons, L.M.; Hultquist, J.F.; Patel, A.B.; Ozer, E.A.; Muller, W.J.; Heald-Sargent, T.; McHugh, M.; et al. Severity of illness caused by severe acute respiratory syndrome coronavirus 2 variants of concern in children: A single-center retrospective cohort study. J. Pediatr. Infect. Dis. 2022, 11, 440–447. [Google Scholar] [CrossRef]
  8. Duong, D. Alpha, Beta, Delta, Gamma: What’s important to know about SARS-CoV-2 variants of concern? Can. Med. Assoc. 2021, 193, E1059–E1060. [Google Scholar] [CrossRef]
  9. Khan, A.; Khan, T.; Ali, S.; Aftab, S.; Wang, Y.; Qiankun, W.; Khan, M.; Suleman, M.; Ali, S.; Heng, W. SARS-CoV-2 new variants: Characteristic features and impact on the efficacy of different vaccines. Biomed. Pharmacother. 2021, 143, 112176. [Google Scholar] [CrossRef]
  10. Wang, Y.; Chen, R.; Hu, F.; Lan, Y.; Yang, Z.; Zhan, C.; Shi, J.; Deng, X.; Jiang, M.; Zhong, S. Transmission, viral kinetics and clinical characteristics of the emergent SARS-CoV-2 Delta VOC in Guangzhou, China. EClinicalMedicine 2021, 40, 101129. [Google Scholar] [CrossRef]
  11. Tian, D.; Sun, Y.; Xu, H.; Ye, Q. The emergence and epidemic characteristics of the highly mutated SARS-CoV-2 Omicron variant. J. Med. Virol. 2022, 94, 2376–2383. [Google Scholar] [CrossRef] [PubMed]
  12. Post, L.A.; Lorenzo-Redondo, R. Omicron: Fewer adverse outcomes come with new dangers. Lancet 2022, 399, 1280–1281. [Google Scholar] [CrossRef] [PubMed]
  13. Frampton, D.; Rampling, T.; Cross, A.; Bailey, H.; Heaney, J.; Byott, M.; Scott, R.; Sconza, R.; Price, J.; Margaritis, M. Genomic characteristics and clinical effect of the emergent SARS-CoV-2 B. 1.1.7 lineage in London, UK: A whole-genome sequencing and hospital-based cohort study. Lancet Infect. Dis. 2021, 21, 1246–1256. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, Y.; Liu, J.; Plante, K.S.; Plante, J.A.; Xie, X.; Zhang, X.; Ku, Z.; An, Z.; Scharton, D.; Schindewolf, C. The N501Y spike substitution enhances SARS-CoV-2 infection and transmission. Nature 2022, 602, 294–299. [Google Scholar] [CrossRef]
  15. Lusvarghi, S.; Wang, W.; Herrup, R.; Neerukonda, S.N.; Vassell, R.; Bentley, L.; Eakin, A.E.; Erlandson, K.J.; Weiss, C.D. Key substitutions in the spike protein of SARS-CoV-2 variants can predict resistance to monoclonal antibodies, but other substitutions can modify the effects. J. Virol. 2022, 96, e01110–e01121. [Google Scholar] [CrossRef]
  16. Jaroszewski, L.; Iyer, M.; Alisoltani, A.; Sedova, M.; Godzik, A. The interplay of SARS-CoV-2 evolution and constraints imposed by the structure and functionality of its proteins. PLoS Comput. Biol. 2021, 17, e1009147. [Google Scholar] [CrossRef]
  17. Alisoltani, A.; Jaroszewski, L.; Iyer, M.; Iranzadeh, A.; Godzik, A. Increased frequency of indels in hypervariable regions of SARS-CoV-2 proteins–a possible signature of adaptive selection. Front. Genet. 2022, 13, 1019. [Google Scholar] [CrossRef]
  18. Shu, Y.; McCauley, J. GISAID: Global initiative on sharing all influenza data—From vision to reality. Eurosurveillance 2017, 22, 30494. [Google Scholar] [CrossRef] [Green Version]
  19. Sedova, M.; Jaroszewski, L.; Alisoltani, A.; Godzik, A. Coronavirus3D: 3D structural visualization of COVID-19 genomic divergence. Bioinformatics 2020, 36, 4360–4362. [Google Scholar] [CrossRef]
  20. Mercatelli, D.; Holding, A.N.; Giorgi, F.M. Web tools to fight pandemics: The COVID-19 experience. Brief. Bioinform. 2021, 22, 690–700. [Google Scholar] [CrossRef]
  21. Portelli, S.; Olshansky, M.; Rodrigues, C.H.; D’Souza, E.N.; Myung, Y.; Silk, M.; Alavi, A.; Pires, D.E.; Ascher, D.B. Exploring the structural distribution of genetic variation in SARS-CoV-2 with the COVID-3D online resource. Nat. Genet. 2020, 52, 999–1001. [Google Scholar] [CrossRef] [PubMed]
  22. Mei, L.-C.; Jin, Y.; Wang, Z.; Hao, G.-F.; Yang, G.-F. Web resources facilitate drug discovery in treatment of COVID-19. Drug Discov. Today 2021, 26, 2358–2366. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, A.T.; Altschuler, K.; Zhan, S.H.; Chan, Y.A.; Deverman, B.E. COVID-19 CG enables SARS-CoV-2 mutation and lineage tracking by locations and dates of interest. Elife 2021, 10, e63409. [Google Scholar] [CrossRef] [PubMed]
  24. Gangavarapu, K.; Latif, A.A.; Mullen, J.L.; Alkuzweny, M.; Hufbauer, E.; Tsueng, G.; Haag, E.; Zeller, M.; Aceves, C.M.; Zaiets, K. Outbreak. info genomic reports: Scalable and dynamic surveillance of SARS-CoV-2 variants and mutations. Res. Sq. 2022, rs. 3, rs-1723829. [Google Scholar]
  25. Hodcroft, E. CoVariants: SARS-CoV-2 Mutations and Variants of Interest. 2021. Available online: https://covariants.org (accessed on 4 October 2022).
  26. Lu, G.; Moriyama, N.E. 2019nCoVR—A comprehensive genomic resource for SARS-CoV-2 variant surveillance. Innovation 2021, 2, 100150. [Google Scholar] [CrossRef]
  27. Singer, J.; Gifford, R.; Cotten, M.; Robertson, D. CoV-GLUE: A Web Application for Tracking SARS-CoV-2 Genomic Variation. 2020. Available online: http://cov-glue.cvr.gla.ac.uk (accessed on 4 October 2022).
  28. Wright, D.W.; Harvey, W.T.; Hughes, J.; Cox, M.; Peacock, T.P.; Colquhoun, R.; Jackson, B.; Orton, R.; Nielsen, M.; Hsu, N.S. Tracking SARS-CoV-2 mutations and variants through the COG-UK-Mutation Explorer. Virus Evol. 2022, 8, veac023. [Google Scholar] [CrossRef]
  29. Tzou, P.L.; Tao, K.; Sahoo, M.K.; Pond, S.L.K.; Pinsky, B.A.; Shafer, R.W. Sierra SARS-CoV-2 sequence and antiviral resistance analysis program. 2022. J. Clin. Virol. 2022, 157, 105323. [Google Scholar] [CrossRef]
  30. Mercatelli, D.; Triboli, L.; Fornasari, E.; Ray, F.; Giorgi, F.M. Coronapp: A web application to annotate and monitor SARS-CoV-2 mutations. J. Med. Virol. 2021, 93, 3238–3245. [Google Scholar] [CrossRef]
  31. Bernasconi, A.; Gulino, A.; Alfonsi, T.; Canakoglu, A.; Pinoli, P.; Sandionigi, A.; Ceri, S. VirusViz: Comparative analysis and effective visualization of viral nucleotide and amino acid variants. Nucleic Acids Res. 2021, 49, e90. [Google Scholar] [CrossRef]
  32. Alsulami, A.F.; Thomas, S.E.; Jamasb, A.R.; Beaudoin, C.A.; Moghul, I.; Bannerman, B.; Copoiu, L.; Vedithi, S.C.; Torres, P.; Blundell, T.L. SARS-CoV-2 3D database: Understanding the coronavirus proteome and evaluating possible drug targets. Brief. Bioinform. 2021, 22, 769–780. [Google Scholar] [CrossRef]
  33. Gowthaman, R.; Guest, J.D.; Yin, R.; Adolf-Bryfogle, J.; Schief, W.R.; Pierce, B.G. CoV3D: A database of high resolution coronavirus protein structures. Nucleic Acids Res. 2021, 49, D282–D287. [Google Scholar] [CrossRef] [PubMed]
  34. Wickham, H. Data Analysis. In Ggplot2; Springer: Cham, Germany, 2016. [Google Scholar] [CrossRef]
  35. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Jespersen, M.C.; Peters, B.; Nielsen, M.; Marcatili, P. BepiPred-2.0: Improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res. 2017, 45, W24–W29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. DeepMind. Computational Predictions of Protein Structures Associated with COVID-19. Available online: https://www.deepmind.com/open-source/computational-predictions-of-protein-structures-associated-with-COVID-19 (accessed on 31 October 2021).
  39. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [Green Version]
  40. Hozumi, Y.; Wang, R.; Yin, C.; Wei, G.W. UMAP-assisted K-means clustering of large-scale SARS-CoV-2 mutation datasets. Comput. Biol. Med. 2021, 131, 104264. [Google Scholar] [CrossRef]
  41. Mullick, B.; Magar, R.; Jhunjhunwala, A.; Farimani, A.B. Understanding mutation hotspots for the SARS-CoV-2 spike protein using Shannon Entropy and K-means clustering. Comput. Biol. Med. 2021, 138, 104915. [Google Scholar] [CrossRef]
  42. Charrad, M.; Ghazzali, N.; Boiteau, V.; Niknafs, A. NbClust: An R package for determining the relevant number of clusters in a data set. J. Stat. Softw. 2014, 61, 1–36. [Google Scholar] [CrossRef] [Green Version]
  43. Deng, X.; Gu, W.; Federman, S.; Du Plessis, L.; Pybus, O.G.; Faria, N.R.; Wang, C.; Yu, G.; Bushnell, B.; Pan, C.-Y. Genomic surveillance reveals multiple introductions of SARS-CoV-2 into Northern California. Science 2020, 369, 582–587. [Google Scholar] [CrossRef]
  44. Russell, A.; O’Connor, C.; Lasek-Nesselquist, E.; Plitnick, J.; Kelly, J.P.; Lamson, D.M.; George, K.S. Spatiotemporal Analyses of 2 Co-Circulating SARS-CoV-2 Variants, New York State, USA. Emerg. Infect. Dis. 2022, 28, 650. [Google Scholar] [CrossRef]
  45. Alpert, T.; Brito, A.F.; Lasek-Nesselquist, E.; Rothman, J.; Valesano, A.L.; MacKay, M.J.; Petrone, M.E.; Breban, M.I.; Watkins, A.E.; Vogels, C.B. Early introductions and transmission of SARS-CoV-2 variant B. 1.1. 7 in the United States. Cell 2021, 184, 2595–2604. [Google Scholar] [CrossRef] [PubMed]
  46. Stevens, L.J.; Pruijssers, A.J.; Lee, H.W.; Gordon, C.J.; Tchesnokov, E.P.; Gribble, J.; George, A.S.; Hughes, T.M.; Lu, X.; Li, J. Mutations in the SARS-CoV-2 RNA dependent RNA polymerase confer resistance to remdesivir by distinct mechanisms. Sci. Transl. Med. 2022, 14, eabo0718. [Google Scholar] [CrossRef] [PubMed]
  47. Posthuma, C.C.; Te Velthuis, A.J.; Snijder, E.J. Nidovirus RNA polymerases: Complex enzymes handling exceptional RNA genomes. Virus Res. 2017, 234, 58–73. [Google Scholar] [CrossRef] [PubMed]
  48. Huang, Y.; Yang, C.; Xu, X.-F.; Xu, W.; Liu, S.-W. Structural and functional properties of SARS-CoV-2 spike protein: Potential antivirus drug development for COVID-19. Acta Pharmacol. Sin. 2020, 41, 1141–1149. [Google Scholar] [CrossRef] [PubMed]
  49. Nagy, Á.; Pongor, S.; Győrffy, B. Different mutations in SARS-CoV-2 associate with severe and mild outcome. Int. J. Antimicrob. Agents 2021, 57, 106272. [Google Scholar] [CrossRef]
  50. Farkas, C.; Mella, A.; Turgeon, M.; Haigh, J.J. 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. Front. Microbiol. 2021, 12, 665041. [Google Scholar] [CrossRef]
  51. Yang, K.; Wang, C.; White, K.I.; Pfuetzner, R.A.; Esquivies, L.; Brunger, A.T. Structural conservation among variants of the SARS-CoV-2 spike postfusion bundle. Proc. Natl. Acad. Sci. USA 2022, 119, e2119467119. [Google Scholar] [CrossRef]
  52. Tao, K.; Tzou, P.L.; Nouhin, J.; Gupta, R.K.; de Oliveira, T.; Kosakovsky Pond, S.L.; Fera, D.; Shafer, R.W. The biological and clinical significance of emerging SARS-CoV-2 variants. Nat. Rev. Genet. 2021, 22, 757–773. [Google Scholar] [CrossRef]
  53. Hadfield, J.; Megill, C.; Bell, S.M.; Huddleston, J.; Potter, B.; Callender, C.; Sagulenko, P.; Bedford, T.; Neher, R.A. Nextstrain: Real-time tracking of pathogen evolution. Bioinformatics 2018, 34, 4121–4123. [Google Scholar] [CrossRef] [Green Version]
  54. Du, P.; Ding, N.; Li, J.; Zhang, F.; Wang, Q.; Chen, Z.; Song, C.; Han, K.; Xie, W.; Liu, J.; et al. Genomic surveillance of COVID-19 cases in Beijing. Nat. Commun. 2020, 11, 5503. [Google Scholar] [CrossRef]
  55. Wang, R.; Chen, J.; Gao, K.; Hozumi, Y.; Yin, C.; Wei, G.W. Author Correction: Analysis of SARS-CoV-2 mutations in the United States suggests presence of four substrains and novel variants. Commun. Biol. 2021, 4, 300. [Google Scholar] [CrossRef] [PubMed]
  56. McCarthy, K.R.; Rennick, L.J.; Nambulli, S.; Robinson-McCarthy, L.R.; Bain, W.G.; Haidar, G.; Duprex, W.P. Recurrent deletions in the SARS-CoV-2 spike glycoprotein drive antibody escape. Science 2021, 371, 1139–1142. [Google Scholar] [CrossRef] [PubMed]
  57. Martin, D.P.; Weaver, S.; Tegally, H.; San, J.E.; Shank, S.D.; Wilkinson, E.; Lucaci, A.G.; Giandhari, J.; Naidoo, S.; Pillay, Y. The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages. Cell 2021, 184, 5189–5200.e7. [Google Scholar] [CrossRef] [PubMed]
  58. Lei, X.; Dong, X.; Ma, R.; Wang, W.; Xiao, X.; Tian, Z.; Wang, C.; Wang, Y.; Li, L.; Ren, L.; et al. Activation and evasion of type I interferon responses by SARS-CoV-2. Nat. Commun. 2020, 11, 3810. [Google Scholar] [CrossRef] [PubMed]
  59. Liang, T.; Cheng, M.; Teng, F.; Wang, H.; Deng, Y.; Zhang, J.; Qin, C.; Guo, S.; Zhao, H.; Yu, X. Proteome-wide epitope mapping identifies a resource of antibodies for SARS-CoV-2 detection and neutralization. Signal Transduct. Target. Ther. 2021, 6, 166. [Google Scholar] [CrossRef]
  60. Smith, C.C.; Olsen, K.S.; Gentry, K.M.; Sambade, M.; Beck, W.; Garness, J.; Entwistle, S.; Willis, C.; Vensko, S.; Woods, A. Landscape and selection of vaccine epitopes in SARS-CoV-2. Genome Med. 2021, 13, 101. [Google Scholar] [CrossRef]
  61. Mou, K.; Mukhtar, F.; Khan, M.T.; Darwish, D.B.; Peng, S.; Muhammad, S.; Al-Sehemi, A.G.; Wei, D.-Q. Emerging Mutations in Nsp1 of SARS-CoV-2 and Their Effect on the Structural Stability. Pathogens 2021, 10, 1285. [Google Scholar] [CrossRef]
  62. Lin, J.W.; Tang, C.; Wei, H.C.; Du, B.; Chen, C.; Wang, M.; Zhou, Y.; Yu, M.X.; Cheng, L.; Kuivanen, S.; et al. Genomic monitoring of SARS-CoV-2 uncovers an Nsp1 deletion variant that modulates type I interferon response. Cell Host Microbe 2021, 29, 489–502.e8. [Google Scholar] [CrossRef]
  63. Ozer, E.A.; Simons, L.M.; Adewumi, O.M.; Fowotade, A.A.; Omoruyi, E.C.; Adeniji, J.A.; Olayinka, O.A.; Dean, T.J.; Zayas, J.; Bhimalli, P.P. Multiple expansions of globally uncommon SARS-CoV-2 lineages in Nigeria. Nat. Commun. 2022, 13, 688. [Google Scholar] [CrossRef]
Figure 1. General workflow of ViralVar and its two main modules. Input data reflecting SARS-CoV-2 sequences of interest can be downloaded directly from the GISAID public repository. In the “Lineage Dynamics” module, the spatiotemporal dynamics of SARS-CoV-2 lineages and clades are represented in the form of stacked bars, area plots, and pie charts over user-selected timeframes and geographical areas. In the “Mutational Analysis” module, mutations are depicted in the context of the SARS-CoV-2 genome and relevant proteins (both primary sequence and 3D structural representations). This module also provides options to perform statistical analyses to identify over- and under mutated proteins over user-selected time periods and perform genome clustering within user- selected subsets. More details are available in the ViralVar User Manual.
Figure 1. General workflow of ViralVar and its two main modules. Input data reflecting SARS-CoV-2 sequences of interest can be downloaded directly from the GISAID public repository. In the “Lineage Dynamics” module, the spatiotemporal dynamics of SARS-CoV-2 lineages and clades are represented in the form of stacked bars, area plots, and pie charts over user-selected timeframes and geographical areas. In the “Mutational Analysis” module, mutations are depicted in the context of the SARS-CoV-2 genome and relevant proteins (both primary sequence and 3D structural representations). This module also provides options to perform statistical analyses to identify over- and under mutated proteins over user-selected time periods and perform genome clustering within user- selected subsets. More details are available in the ViralVar User Manual.
Viruses 14 02714 g001
Figure 2. Area plots reflecting (different colors represent variants of concern see legend for details) the relative abundance of variants of concern and variants of interest collected in the USA over time. (A) Frequency of indicated VOCs and VOIs over time in specimens collected between January 2020 and May 2022 in the USA (n = 1,739,797 sequences from GISAID as of 31 May 2022). (B) Frequency of indicated VOCs and VOIs over time in specimens collected between November 2020 and March 2021. Specimens were divided into three age groups: children (up to 18 years), adults (18–64 years), and the elderly (65 years or more). The number of sequences per age group is indicated above each plot. Each subset of genomes was processed separately using the ViralVar “Lineage Dynamics” module.
Figure 2. Area plots reflecting (different colors represent variants of concern see legend for details) the relative abundance of variants of concern and variants of interest collected in the USA over time. (A) Frequency of indicated VOCs and VOIs over time in specimens collected between January 2020 and May 2022 in the USA (n = 1,739,797 sequences from GISAID as of 31 May 2022). (B) Frequency of indicated VOCs and VOIs over time in specimens collected between November 2020 and March 2021. Specimens were divided into three age groups: children (up to 18 years), adults (18–64 years), and the elderly (65 years or more). The number of sequences per age group is indicated above each plot. Each subset of genomes was processed separately using the ViralVar “Lineage Dynamics” module.
Viruses 14 02714 g002
Figure 3. Distribution of SARS-CoV-2 variants of concern and interest by US state. Pie charts represent the proportion of SARS-CoV-2 VOIs and VOCs in each US state as reported to GISAID (as of 31 May 2022) (A) between January 2020 and May 2022 and (B) between November 2020 and March 2021. Specimens were divided into three age groups: children (up to 18 years), adults (18–64 years), and the elderly (65 years or more). The number of sequences per age group is indicated above each plot. The size of pie charts represents the relative frequency of sequenced data in each state. Each subset of genomes was visualized separately using the ViralVar “Geographical Map” feature.
Figure 3. Distribution of SARS-CoV-2 variants of concern and interest by US state. Pie charts represent the proportion of SARS-CoV-2 VOIs and VOCs in each US state as reported to GISAID (as of 31 May 2022) (A) between January 2020 and May 2022 and (B) between November 2020 and March 2021. Specimens were divided into three age groups: children (up to 18 years), adults (18–64 years), and the elderly (65 years or more). The number of sequences per age group is indicated above each plot. The size of pie charts represents the relative frequency of sequenced data in each state. Each subset of genomes was visualized separately using the ViralVar “Geographical Map” feature.
Viruses 14 02714 g003
Figure 4. Absolute frequency of mutations in SARS-CoV-2 Spike and NSP12 among Alpha VOCs. SARS-CoV-2 genome data for all sequences assigned to an Alpha variant lineage (B.1.1.7 and Q.*) from the USA (n = 140,100) and rest of the world (n = 906,114, the USA cases excluded) were retrieved from GISAID as of 31 May 2022. Plots represent the absolute frequency of mutations at each amino acid position across Spike (top) and NSP12 (bottom) in sequences from the USA (left) and rest of the world (right). Deletions (red), insertions (blue), and substitutions (green) are plotted in different colors at each position. Boundaries for protein domains of Spike and NSP12 proteins were obtained from [16,17]. Predicted B-cell epitopes are highlighted above in teal, as predicted by [36]. Each subset of genomes was visualized separately using the ViralVar “Protein Distribution” feature.
Figure 4. Absolute frequency of mutations in SARS-CoV-2 Spike and NSP12 among Alpha VOCs. SARS-CoV-2 genome data for all sequences assigned to an Alpha variant lineage (B.1.1.7 and Q.*) from the USA (n = 140,100) and rest of the world (n = 906,114, the USA cases excluded) were retrieved from GISAID as of 31 May 2022. Plots represent the absolute frequency of mutations at each amino acid position across Spike (top) and NSP12 (bottom) in sequences from the USA (left) and rest of the world (right). Deletions (red), insertions (blue), and substitutions (green) are plotted in different colors at each position. Boundaries for protein domains of Spike and NSP12 proteins were obtained from [16,17]. Predicted B-cell epitopes are highlighted above in teal, as predicted by [36]. Each subset of genomes was visualized separately using the ViralVar “Protein Distribution” feature.
Viruses 14 02714 g004
Figure 5. Analysis of Alpha VOC mutations predominantly found in USA specimens. (A) Relative frequency of Alpha VOC SARS-CoV-2 genomes harboring Spike K1191N (top) or NSP12 P227L (bottom) mutations. Calculations are based on GISAID data as of 31 May 2022. (B) Spike K1191N (top) or NSP12 P227L (bottom) mutations highlighted on available protein structures using the ViralVar “3D Protein Structure” feature. The Spike receptor-binding domain (RBD) and N-terminal domain (NTD) are colored in light blue and pink, respectively (D-I-TASSER model). NSP12 is colored in light blue with the NiRAN domain highlighted in pink (PDB: 6XEZ). (C) Euclidean distance-based k-means clustering of Alpha VOC SARS-CoV-2 genomes based on Spike and NSP12 mutations was performed using the “Genome Clustering” feature. Heatmaps represent the percent of genomes with a specific mutation within each cluster. Only protein mutations present in more than two thirds (70%) of genomes are shown here. (D) Time-resolved phylogenetic tree built by Nextstrain (https://nextstrain.org/ncov/gisaid/north-america/) accessed 14 October 2021 using a North America-focused subsampling between December 2020 and August 2021 (n = 399 sequences) visualized using R package ggtree. Yellow and hot pink branches and tips highlight genomes containing the Spike K1191N and NSP12 P227L mutations, respectively.
Figure 5. Analysis of Alpha VOC mutations predominantly found in USA specimens. (A) Relative frequency of Alpha VOC SARS-CoV-2 genomes harboring Spike K1191N (top) or NSP12 P227L (bottom) mutations. Calculations are based on GISAID data as of 31 May 2022. (B) Spike K1191N (top) or NSP12 P227L (bottom) mutations highlighted on available protein structures using the ViralVar “3D Protein Structure” feature. The Spike receptor-binding domain (RBD) and N-terminal domain (NTD) are colored in light blue and pink, respectively (D-I-TASSER model). NSP12 is colored in light blue with the NiRAN domain highlighted in pink (PDB: 6XEZ). (C) Euclidean distance-based k-means clustering of Alpha VOC SARS-CoV-2 genomes based on Spike and NSP12 mutations was performed using the “Genome Clustering” feature. Heatmaps represent the percent of genomes with a specific mutation within each cluster. Only protein mutations present in more than two thirds (70%) of genomes are shown here. (D) Time-resolved phylogenetic tree built by Nextstrain (https://nextstrain.org/ncov/gisaid/north-america/) accessed 14 October 2021 using a North America-focused subsampling between December 2020 and August 2021 (n = 399 sequences) visualized using R package ggtree. Yellow and hot pink branches and tips highlight genomes containing the Spike K1191N and NSP12 P227L mutations, respectively.
Viruses 14 02714 g005
Figure 6. Genomic distribution of SARS-CoV-2 mutations for three age groups. Each plot depicts the number of distinct protein mutations in a 100 nucleotide sliding window across the SARS-CoV-2 genome in specimens collected between January 2020 and May 2022 in the USA (n = 1,739,797 sequences from GISAID as of 31 May 2022). Sequences were divided into six groups based on the age of patients (children (up to 18 years), adults (18–64 years), and elderly (65 years or more)). The total number of sequences used per age group is indicated. Each subset of genomes was processed separately using the ViralVar “Lineage Dynamics” module.
Figure 6. Genomic distribution of SARS-CoV-2 mutations for three age groups. Each plot depicts the number of distinct protein mutations in a 100 nucleotide sliding window across the SARS-CoV-2 genome in specimens collected between January 2020 and May 2022 in the USA (n = 1,739,797 sequences from GISAID as of 31 May 2022). Sequences were divided into six groups based on the age of patients (children (up to 18 years), adults (18–64 years), and elderly (65 years or more)). The total number of sequences used per age group is indicated. Each subset of genomes was processed separately using the ViralVar “Lineage Dynamics” module.
Viruses 14 02714 g006
Table 1. Details of SARS-CoV-2 data used in this study. Data retrieved from GISAID and each of the three data subsets were separately analyzed using ViralVar.
Table 1. Details of SARS-CoV-2 data used in this study. Data retrieved from GISAID and each of the three data subsets were separately analyzed using ViralVar.
SequencesMean AgeMedian Age
Children (<18)282,10610.2210.5
Adults (18–65)1,287,05838.9237.5
Elderly (>65)170,63374.4272.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alisoltani, A.; Jaroszewski, L.; Godzik, A.; Iranzadeh, A.; Simons, L.M.; Dean, T.J.; Lorenzo-Redondo, R.; Hultquist, J.F.; Ozer, E.A. ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes. Viruses 2022, 14, 2714. https://doi.org/10.3390/v14122714

AMA Style

Alisoltani A, Jaroszewski L, Godzik A, Iranzadeh A, Simons LM, Dean TJ, Lorenzo-Redondo R, Hultquist JF, Ozer EA. ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes. Viruses. 2022; 14(12):2714. https://doi.org/10.3390/v14122714

Chicago/Turabian Style

Alisoltani, Arghavan, Lukasz Jaroszewski, Adam Godzik, Arash Iranzadeh, Lacy M. Simons, Taylor J. Dean, Ramon Lorenzo-Redondo, Judd F. Hultquist, and Egon A. Ozer. 2022. "ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes" Viruses 14, no. 12: 2714. https://doi.org/10.3390/v14122714

APA Style

Alisoltani, A., Jaroszewski, L., Godzik, A., Iranzadeh, A., Simons, L. M., Dean, T. J., Lorenzo-Redondo, R., Hultquist, J. F., & Ozer, E. A. (2022). ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes. Viruses, 14(12), 2714. https://doi.org/10.3390/v14122714

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop