Next Article in Journal
Identification and Candidate Gene Evaluation of a Large Fast Neutron-Induced Deletion Associated with a High-Oil Phenotype in Soybean Seeds
Next Article in Special Issue
Using Genetics to Investigate Relationships between Phenotypes: Application to Endometrial Cancer
Previous Article in Journal
Screening of CIMMYT and South Asian Bread Wheat Germplasm Reveals Marker–Trait Associations for Seedling Resistance to Septoria Nodorum Blotch
Previous Article in Special Issue
A Polygenic Risk Analysis for Identifying Ulcerative Colitis Patients with European Ancestry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Virus Variant Detection Based on the Optimal Natural Metric

1
Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China
2
Beijing Institute of Mathematical Sciences and Applications (Bimsa), Beijing 101408, China
*
Author to whom correspondence should be addressed.
Genes 2024, 15(7), 891; https://doi.org/10.3390/genes15070891
Submission received: 5 June 2024 / Revised: 3 July 2024 / Accepted: 5 July 2024 / Published: 7 July 2024
(This article belongs to the Special Issue Statistical Methods for Genetic Epidemiology)

Highlights

What are the main findings?
  • We introduce a new algorithm designed for the automatic detection of emerging virus variants.
  • The algorithm was tested on real datasets including SARS-CoV-2 and HIV-1, demonstrating nearly 100% precision in identification.
What is the implication of the main finding?
  • Our approach enables the efficient identification of new virus variants based solely on sequence data, eliminating the need for biologists to pinpoint key viral regions.
  • Our method pushes the boundaries of alignment-free techniques, expanding their application from classifying within known categories to recognizing new categories.

Abstract

:
The highly variable SARS-CoV-2 virus responsible for the COVID-19 pandemic frequently undergoes mutations, leading to the emergence of new variants that present novel threats to public health. The determination of these variants often relies on manual definition based on local sequence characteristics, resulting in delays in their detection relative to their actual emergence. In this study, we propose an algorithm for the automatic identification of novel variants. By leveraging the optimal natural metric for viruses based on an alignment-free perspective to measure distances between sequences, we devise a hypothesis testing framework to determine whether a given viral sequence belongs to a novel variant. Our method demonstrates high accuracy, achieving nearly 100% precision in identifying new variants of SARS-CoV-2 and HIV-1 as well as in detecting novel genera in Orthocoronavirinae. This approach holds promise for timely surveillance and management of emerging viral threats in the field of public health.

1. Introduction

Viruses are microorganisms that depend on host cells for replication, and are characterized by their high variability [1]. This trait is especially evident in RNA viruses such as Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) and Human Immunodeficiency Virus-1 (HIV-1) [2,3]. The rapid mutation rates of these viruses frequently lead to the emergence of new variants with distinct pathogenicity or transmissibility, posing significant challenges to public health efforts, as evidenced by the COVID-19 pandemic [4].
Therefore, the timely identification of new variants has become a crucial issue in biology. However, current processes for discovering and defining new variants often require manual intervention. After a new sequence is obtained, it is aligned with known sequences at biologically significant regions (e.g., segments corresponding to the spike protein) in order to identify mutation sites [5]. Subsequently, official organizations such as the World Health Organization (WHO) define variants based on key mutations [6,7]. For example, the Omicron variant of SARS-CoV-2 is characterized by a series of key mutations such as K417N, T478K, and N501Y [8]. This system proved vital during the COVID-19 pandemic; however, developing a more automated variant identification system remains a valuable goal.
Our theoretical framework differs from the aforementioned approaches. It is rooted in alignment-free principles, meaning that manual identification of important segments is not required and that there is no need to spend time on sequence alignment in order to pinpoint mutations. Instead, we assess whether a new sequence differs from known variants by analyzing its statistical characteristics. As sequencing technologies advance, the pool of known sequences is expanding, amplifying the importance of more efficient alignment-free methods [9,10,11].
Among the current alignment-free methods, the natural vector approach stands out. This approach utilizes statistical moments to transform sequences into feature vectors based on different k-mers and moment orders, and has garnered widespread attention [12,13,14]. It has demonstrated significant efficacy across such diverse fields as sequence classification, phylogenetic analysis, and chromosome fusion detection [15,16,17]. The optimal natural metric, an advancement of the natural vector method, fine tunes the integration of statistical information from various k-mers and moment orders [18]. It identifies the optimal combination of features through training on real data, thereby naturally defining dissimilarities between sequences. However, no previous applications of these metrics have addressed the discovery of new categories that do not belong to the existing classification system.
In this paper, we employ the optimal natural metric to measure the distances between sequences, then extend the metric to effectively assess the distances between individual sequences and known variants. We observe distinct distributional characteristics in the distances from sequences to their non-affiliated variants, allowing us to determine whether a sequence belongs to a given variant and if it originates from an unknown new variant through hypothesis testing. We validate this method on two significant viruses, SARS-CoV-2 and HIV-1, and further analyze different genera within the Orthocoronavirinae family at a higher classification level. These numerical experiments show very promising results, with minimal Type I and Type II errors, confirming the effectiveness of our method. This approach offers significant potential for the timely detection and management of emerging viral threats in public health.

2. Materials and Methods

2.1. Materials

We utilized three types of datasets in this study. For all datasets, sequences containing ambiguous letters and sequences with unidentified categories were removed. All sequences along with their information are available in the Github repository https://github.com/BobYHY/NewVariant (accessed on 29 June 2024).
SARS-CoV-2 Dataset: The first dataset was from GISAID (https://gisaid.org accessed on 30 June 2022), encompassing SARS-CoV-2 sequences up to 30 June 2022. We included sequences categorized as variants of concern (VOC) and variants of interest (VOI). Due to the imbalance in the number of sequences for different variants (e.g., 99,955 sequences for Omicron and only 3 for Theta), we excluded variants with fewer than ten sequences and randomly sampled 500 sequences for variants with more than 500 sequences. This resulted in a final dataset of 4917 sequences from twelve variants.
HIV-1 Dataset: The second dataset was sourced from the HIV Database (https://www.hiv.lanl.gov accessed on 8 April 2022), containing HIV-1 sequences up to 8 April 2022. Initially, we considered all non-recombinant subtype cDNA sequences available in the database and excluded subtypes with fewer than ten sequences, resulting in a dataset of 5549 sequences from six subtypes. Additionally, we included 117 recombinant viruses belonging to combinations of A+B, A+C, and B+C to further explore our method.
Orthocoronavirinae Dataset: The third dataset was sourced from NCBI (https://www.ncbi.nlm.nih.gov/ accessed on 17 November 2022), comprising Orthocoronavirinae sequences available up to 17 November 2022. Except for SARS-CoV-2, where we randomly selected 500 sequences due to the high number of available sequences, all Orthocoronavirinae sequences with lengths between 25,000 and 35,000 bases were included. This process resulted in a final dataset of 3169 sequences from four genera.

2.2. The Optimal Natural Metric

The optimal natural metric is used to measure the dissimilarity between sequences based on natural vectors. It leverages statistical moments to extract features from sequences, followed by training of the best weighted scheme for these features based on actual datasets [18].
The natural vector method is an alignment-free approach that transforms DNA sequences into vectors of moments [12]. Consider the sequence S = s 1 s 2 s n , where s i , α { A , T , C , G } ; we can define
w α ( s i ) = 1 , s i = α 0 , o t h e r w i s e .
Then, the j-th ordered element D j α of the natural vector can be defined as
D 0 α = n α = i = 1 n w α ( s i ) D 1 α = μ α = i = 1 n i n α w α ( s i ) D j α = i = 1 n ( i μ α ) j n α j 1 n j 1 w α ( s i ) ( j = 2 , 3 , 4 , ) ,
where n = n A + n T + n C + n G . The distance between j-th ordered elements of two natural vectors refers to the distance between their corresponding ( D j A , D j T , D j C , D j G ) .
The k-mer natural vector method extends the natural vector method [13]. A k-mer is a string composed of k nucleotides, resulting in 4 k possible k-mers, denoted as l 1 , , l 4 k . For a sequence S = s 1 s 2 s n , we can view it as a sequence consisting of n k + 1 k-mers, i.e.,  ( s 1 s k ) ( s n k + 1 s n ) . We can define the j-th ordered moments of k-mers similar to those of nucleotides, and thereby define the k-mer natural vector (if D 0 l i = 0 , then we set D 1 l i = D 2 l i = = 0 ).
Given two DNA sequences a and b, the Euclidean distance between the j-th ordered elements of their k-mer natural vectors is defined as d i s k j ( a , b ) . With an arbitrarily chosen weight matrix w k j , we can define a weighting metric as
D i s w ( a , b ) = k = 1 K j = 0 J w k j d i s k j ( a , b ) .
In this metric, different d i s k j are used to capture differences in various types of information within sequences. By integrating them, a more comprehensive measurement between sequences can be obtained. We aim to achieve the best classification performance of the metric under a set of weights, where the performance is measured by the 1NN accuracy [19]. Specifically, we seek to ensure that sequences and their nearest neighbors belong to the same class under the chosen metric. We can smooth the 1NN accuracy and use a gradient-based algorithm to calculate the optimal weights. In applications, we set K = 9 and J = 2 , and train these weights based on all complete virus reference sequences [18]. The specific optimal weight matrix are shown in the Table 1. In the subsequent text, we refer to the metric corresponding to these weights as Dis.

2.3. New Virus Detection Method

In the preceding section, we have defined the metric Dis between individual sequences. In the subsequent discussion, we extend this definition to encompass the distance from a single sequence to a set of sequences.
Considering that the known classification data of sequences may occasionally contain errors, we have devised a robust distance measure D i s R by computing the distance from a sequence to its R-th nearest neighbor in the set, which serves as the distance from the sequence to the set of sequences.
Assuming that all sequences can be divided into M categories, denoted as A 1 ,…, A M , we can define the following sets of numbers:
I o u t ( A i , A j ) : = { D i s R ( a i , A j ) | a i A i } ( i j ) I i n ( A i ) : = { D i s R ( a i , A i { a i } ) | a i A i } I o u t : = i j I o u t ( A i , A j ) I i n : = i = 1 M I i n ( A i )
where I o u t ( A i , A j ) reflects the distribution of D i s R from sequences in A i to the set A j and its union I o u t reflects the distance characteristics when sequences are not within a certain category.
Moreover, I i n ( A i ) reflects the distribution of D i s R from sequences in A i to the set composed of other sequences within the same category, while its union I i n reflects the distance characteristics when sequences are within a certain category.
In practical applications, significant differences are observed between the two sets of numbers I i n and I o u t . Specifically, they can be separated by a certain threshold B. Suppose that a new sequence a appears; if  D i s R ( a , A i ) > B , then a can be considered to not belonging to A i . If  min i = 1 , , M D i s R ( a , A i ) > B , then a does not belong to any known category, indicating that it belongs to a new category.
A straightforward choice for B is B = max I i n ; however, we aim for greater robustness of the algorithm to the data. Therefore, we consider the 99 % quantile of I i n , i.e., we select B such that B is just greater than 99 % of values in I i n .
The requirement for robustness is positively correlated with the amount of errors in the dataset. Therefore, for the selection of R, we associate it with the number of sequences in the largest category. We define R : = 1 + max i c a r d ( A i ) 1000 , where c a r d ( A i ) denotes the number of sequences in class A i and x means the nearest smaller integer for x.
We summarize our method in Algorithm 1. The first four steps can be performed in advance before the arrival of a new sequence, and these calculations need to be performed only once with multiple query sequences, ensuring high efficiency. The code can be found in both the Supplementary Materials and the Github repository https://github.com/BobYHY/NewVariant (accessed on 29 June 2024). In Figure 1, we present a flowchart to assist in understanding how we obtained our results from the original data.
We evaluate the performance of the algorithm from two perspectives. In statistical terms, hypothesis testing errors are divided into Type I errors and Type II errors [20]. Here, a Type I error refers to failing to identify sequences from a new category, while a Type II error refers to misidentifying sequences from an existing category as sequences from a new category. For a given dataset, we follow the leave-one-out approach to test the algorithm’s error rate. When assessing the Type I error rate, we examine each sequence v i in succession as a new sequence. Here, all sequences except those in the same category as v i are regarded as known sequences, and the probability of failing to detect the new virus is calculated. In evaluating the Type II error rate, we again consider each sequence in sequence as a new sequence. This time, all other known sequences, including those in the same category as v i , are considered known, and the probability of incorrectly categorizing the sequence as belonging to a new category is computed.
Algorithm 1 New variant detection method
Require: 
Sequences from M known categories A 1 , , A M . New sequence v.
1:
Compute pairwise distances for known sequences using the optimal natural metric.
2:
Determine R : = 1 + max i c a r d ( A i ) 1000 .
3:
Compute D i s R for all sequences with respect to their respective categories to obtain I i n .
4:
Determine B based on the 99 % quantile of I i n .
5:
Given a new sequence v, calculate d = min i = 1 , , M D i s R ( v , A i ) .
6:
if  d > B   then
7:
   Claim that v belongs to a new category.
8:
else
9:
   Assert that v does not belong to a new category and its classification result can be obtained from its nearest neighbors.
10:
end if

3. Results

3.1. SARS-CoV-2

The most significant variants of SARS-CoV-2 belong to the category of variants of concern (VOC), which includes Alpha, Beta, Gamma, Delta, and Omicron. Additionally, there are variants of interest (VOI) that are also noteworthy. After removing variants with insufficient sequence counts, seven VOIs remain: Epsilon, Zeta, Eta, Iota, Kappa, Lambda, and Mu. These twelve variants comprise the 4917 sequences that we analyze.
To begin with, we use the optimal natural metric to classify sequences by the 1NN algorithm with the leave-one-out strategy, achieving an accuracy of 99.9%. This high accuracy demonstrates that the optimal natural metric effectively reflects the relationships between SARS-CoV-2 RNA sequences.
Next, we illustrate that the optimal natural metric is a viable approach for detecting new variants. In Figure 2, we use kernel density estimation to show the distribution differences between I i n ( A l p h a ) and I o u t ( · , A l p h a ) [21]. It is evident that the distribution of D i s R ( · , A l p h a ) for sequences within variant Alpha falls within the threshold, forming a unimodal distribution. In contrast, D i s R ( · , A l p h a ) for sequences from other variants are outside this threshold. Figure 3 further illustrates the distributional differences between I i n and I o u t . The distinct separation between these two distributions indicates that our threshold can effectively distinguish whether a sequence belongs to a specific variant class.
We test Algorithm 1 for various SARS-CoV-2 variants using the previously described strategies for evaluating Type I and Type II errors. In the Type I error test, we assume that sequences from a given variant have never been discovered before and assess whether a newly appearing sequence from that variant can be correctly identified as new. In the Type II error test, we evaluate whether a sequence from known variants can be mistakenly identified as a new variant. As shown in Table 2, both types of errors are kept below 1%, which indicates a low rate of false positives and false negatives, yielding highly satisfactory results.

3.2. HIV-1

HIV-1 is categorized into a major group (M) and three minor groups (O, N, P). The major group is further subdivided into subtypes: A, B, C, D, E, F, G, H, I, J, K, L [22,23,24,25,26]. After excluding subtypes with insufficient sequence counts, six subtypes remain: A, B, C, D, F, G. These six variants comprise the 5549 sequences that we analyze.
We repeat the aforementioned analysis for the HIV-1 virus. Initially, we determine that the 1NN classification accuracy using the optimal natural metric is 99.8%, demonstrating a remarkably high level of accuracy. In Figure A1, significant disparities can be observed between I i n and I o u t , consistent with our previous findings. Subsequently, we evaluate Algorithm 1, recording Type I and Type II errors at 0.94% and 0.87%, respectively. These results indicate promising performance of the algorithm in classifying HIV-1 variants.
One notable characteristic of HIV-1 is the significant attention given to recombinant viruses, which exhibit properties of multiple subtypes without any single subtype being dominant (if one subtype is dominant, then it is classified as pure) [27]. We aim to investigate how these recombinant viruses perform with our algorithm. Specifically, we consider three combinations: A+B, A+C, and B+C, comprising 10, 38, and 69 sequences, respectively.
Generally speaking, for recombinant viruses of the X+Y type, there are four potential identification outcomes. The first outcome is when the virus is identified as either X or Y, indicating a closer similarity to one of the recombinant parents. This suggests that the virus should be assigned to the closer category. The second outcome occurs when the virus is recognized as within both X and Y, with distances to both variants falling within the boundary. This suggests similarly short distances to both parents. The third outcome is when the virus is identified as a new variant, with distances to all known variants falling outside the boundary. This indicates significant differences between the virus and both of its parents, suggesting that it should be regarded as a new variant instead of solely a recombinant. The fourth outcome occurs when none of the above three scenarios apply, such as being classified as Z type or X+Z type. Only this fourth scenario implies an error, as it contradicts the actual situation.
Out of the 117 recombinant sequences we analyzed, only three fell into the fourth category, representing misclassifications, with the remaining sequences providing valuable biological classification suggestions. This underscores the overall effectiveness of our method. In the following parts, we illustrate the assistance that our method can provide using specific examples.
For the A+B recombinant sequences, those identified as new sequences all belong to subtype A1B, while another subtype, 03_A6B, is mostly identified as belonging to both A and B. In Figure 4, we further illustrate this result using the points ( D i s R ( · , A ) , D i s R ( · , B ) ) for different subtypes. This suggests that A1B is more distantly related compared to our given dataset. For the A+C recombinant sequences, A2C is consistently recognized as a new variant, whereas A1C shows inconsistent results, indicating that A2C is a more internally homogeneous subtype and is more distantly related to its parent variants. For the B+C recombinant sequences, six out of eight subtypes are consistently recognized as C, implying that these subtypes are more similar to variant C. These results can provide valuable insights into different subtypes within a recombinant group.

3.3. Orthocoronavirinae

Coronaviridae is a family of positive-sense single-stranded RNA viruses that cause diseases in mammals and birds [28]. They are named for their characteristic crown-like appearance, attributed to the presence of spike proteins on their surface [29]. Within the Coronaviridae family, the Orthocoronavirinae subfamily represents the vast majority, comprising four genera: Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus [30].
In our previous analyses, we have concentrated on detecting new variants of specific viruses. Now, we expand our focus to include the identification of new genera. Initially, we calculate the 1NN classification accuracy at the genus level using the optimal natural metric, achieving a perfect score of 99.9%. In Figure A2, it can be noted that I i n and I o u t display even greater disparities compared to the previous analysis at the variant level. This is easily comprehensible, as the discrepancies between genera are more pronounced than those between variants.
Subsequently, we assess Algorithm 1, recording Type I and Type II errors at 0.03% and 0.98%, respectively. This outcome slightly surpasses the previous viral analyses, signifying the efficacy of our method at higher classification levels.

3.4. Time Complexity Analysis

Our approach incorporates an alignment-free perspective, significantly reducing its time complexity. We assume a total of N sequences with known labels, each with a length of O ( L ) , along with M known variants, where each variant comprises O ( N / M ) sequences alongside V new sequences with unknown labels. Given that the computational complexity of a single k-mer natural vector is O ( L ) , the overall time complexity of the algorithm is O ( N L + N 2 / M + V ( L + N ) ) , with three terms corresponding to calculating the natural vectors, calculating the boundary, and computation upon the arrival of new sequences, respectively.
When using the alignment method, although search algorithms such as BLAST can quickly identify candidate similar sequences, calculating the actual distance necessitates algorithms with a time complexity of O ( L 2 ) [31,32]. Moreover, the distance must be computed between every pair of sequences, while the alignment-free approach requires only a single embedding per sequence. Therefore, if the distance calculation component of our method were replaced with the alignment method, the overall algorithmic complexity would increase to O ( N 2 L 2 / M + V N L 2 ) , with two terms representing the boundary calculation and computation upon the arrival of new sequences, respectively. This is is much slower than alignment-free approaches.
To conduct a further comparison using a small sample, we randomly selected thirty sequences from each of the five variants within the VOC of SARS-CoV-2, compiling a dataset of 150 sequences. We then applied both the original method and the alignment-based approach to perform the same leave-one-out test as mentioned earlier. The time expenditure of the original method on a standard PC with CPU 2.3 GHz was 2 min and 15 s, and the Type I error and Type II error of the test were 2.00% and 1.33%, respectively. However, using the alignment-based method (employing the default parameters for global alignment in the Align.PairwiseAligner function from Biopython (version 1.7.9), with the additive inverse of the match score as the similarity metric) takes 12 h and 37 min, and the Type I error reaches 97%, indicating that the method is completely ineffective. It is worth noting that alignment methods offer various parameter options, and there may be certain configurations that can potentially make the method work. In terms of time efficiency alone, however, it is nearly impossible to process large datasets using this approach.

4. Discussion

In this paper, we have employed the optimal natural metric, which demonstrated exceptional performance across three distinct datasets, to measure the similarity between sequences. Using this metric, we have developed a hypothesis testing-based algorithm capable of quickly determining whether a sequence belongs to an unknown category. Our approach achieves nearly 100% accuracy in detecting new variants of SARS-CoV-2, new variants of HIV-1, and new genera within Orthocoronavirinae. Additionally, for HIV-1 our approach provides deeper insights into recombinant viruses by evaluating their relationships with their parental strains.
The significance of our method lies in two main aspects. First, in an era where pandemics are of increasing concern, the ability to rapidly detect new variants is crucial. Current methods for defining new variants rely heavily on manual selection, requiring the identification of key regions for alignment and specific mutations. Our method offers an automated approach to identifying new variants, reducing dependence on labor-intensive processes. In the process of infectious disease control and management, even after defining the known variants, there is always concern about the potential emergence of new threatening variants. Our method allows for the efficient detection of new variants based purely on sequence data, without the need to rely on biologists to identify critical regions of the virus. This capability is of significant importance.
Second, our method advances alignment-free techniques. Alignment-free methods which do not require time-consuming alignment processes and are not dependent on conserved regions are becoming increasingly important in the era of big data. Previous research has proposed various methods for extracting and comparing sequence features, with the optimal natural metric being one of the most effective. However, these methods typically only classify sequences when all categories are known, failing to address the presence of new categories. Our approach fills this gap by enabling the identification of sequences from previously unknown categories.
Certainly, our method has limitations that require further investigation in future studies. A primary issue is that our approach handles relatively complete sequences, and cannot determine whether a shorter unknown segment of a virus belongs to a new category. Therefore, while our method is very useful for disease control and management where complete sequences can be obtained, it is currently less effective in areas where complete sequences are not available. For example, the study of meta-viromes must often deal with a large number of mixed short sequences from different pathogens. If future research could extend our algorithm to make it effective on fragment sequences as well, this would be highly beneficial.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15070891/s1, Data S1: SARSCoV2.csv; Data S2: SARSCoV2.fasta; Data S3: HIV1.csv; Data S4: HIV1.fasta; Data S5: Orthocoronavirinae.csv; Data S6: Orthocoronavirinae.fasta; Code S1: metric.py; Code S2: algorithm1.py.

Author Contributions

Conceptualization, H.Y. and S.S.-T.Y.; methodology, H.Y. and S.S.-T.Y.; software, H.Y.; validation, H.Y.; formal analysis, H.Y.; investigation, H.Y.; resources, S.S.-T.Y.; data curation, H.Y.; writing—original draft preparation, H.Y.; writing—review and editing, S.S.-T.Y.; visualization, H.Y.; supervision, S.S.-T.Y.; project administration, S.S.-T.Y.; funding acquisition, S.S.-T.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by National Natural Science Foundation of China (NSFC) grant (12171275) and the Tsinghua University Education Foundation Fund (042202008).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The code and data are available in the Github repository https://github.com/BobYHY/NewVariant (accessed on 29 June 2024).

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study, in the collection, analysis, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
COVID-19Coronavirus Disease 2019
SARS-CoV-2Severe Acute Respiratory Syndrome Coronavirus-2
HIV-1Human Immunodeficiency Virus-1
WHOWorld Health Organization
VOCVariants Of Concern
VOIVariants Of Interest
NCBINational Center for Biotechnology Information
NNNearest Neighbor

Appendix A

Figure A1. Distribution of I i n and I o u t (HIV-1).
Figure A1. Distribution of I i n and I o u t (HIV-1).
Genes 15 00891 g0a1
Figure A2. Distribution of I i n and I o u t (Orthocoronavirinae).
Figure A2. Distribution of I i n and I o u t (Orthocoronavirinae).
Genes 15 00891 g0a2

References

  1. LaTourrette, K.; Garcia-Ruiz, H. Determinants of Virus Variation, Evolution, and Host Adaptation. Pathogens 2022, 11, 1039. [Google Scholar] [CrossRef]
  2. Uddin, M.; Mustafa, F.; Rizvi, T.A.; Loney, T.; Al Suwaidi, H.; Al-Marzouqi, A.H.H.; Kamal Eldin, A.; Alsabeeha, N.; Adrian, T.E.; Stefanini, C.; et al. SARS-CoV-2/COVID-19: Viral Genomics, Epidemiology, Vaccines, and Therapeutic Interventions. Viruses 2020, 12, 526. [Google Scholar] [CrossRef]
  3. Maartens, G.; Celum, C.; Lewin, S.R. HIV infection: Epidemiology, pathogenesis, treatment, and prevention. Lancet 2014, 384, 258–271. [Google Scholar] [CrossRef] [PubMed]
  4. Vulturar, D.-M.; Moacă, L.-Ș.; Neag, M.A.; Mitre, A.-O.; Alexescu, T.-G.; Gherman, D.; Făgărășan, I.; Chețan, I.M.; Gherman, C.D.; Melinte, O.-E.; et al. Delta Variant in the COVID-19 Pandemic: A Comparative Study on Clinical Outcomes Based on Vaccination Status. J. Pers. Med. 2024, 14, 358. [Google Scholar] [CrossRef]
  5. 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]
  6. Li, M.; Lou, F.; Fan, H. SARS-CoV-2 Variants of Concern Delta: A great challenge to prevention and control of COVID-19. Sig. Transduct. Target Ther. 2021, 6, 349. [Google Scholar] [CrossRef] [PubMed]
  7. Enhancing Response to Omicron SARS-CoV-2 Variant. Available online: https://www.who.int/publications/m/item/enhancing-readiness-for-omicron-(b.1.1.529)-technical-brief-and-priority-actions-for-member-states (accessed on 31 May 2024).
  8. Karim, S.S.A.; Karim, Q.A. Omicron SARS-CoV-2 variant: A new chapter in the COVID-19 pandemic. Lancet 2022, 399, 2126–2128. [Google Scholar]
  9. Zielezinski, A.; Vinga, S.; Almeida, J.S.; Karłowski, W.M. Alignment-free sequence comparison: Benefits, applications, and tools. Genome Biol. 2017, 18, 186. [Google Scholar] [CrossRef] [PubMed]
  10. Bonham-Carter, O.; Steele, J.; Bastola, D.R. Alignment-free genetic sequence comparisons: A review of recent approaches by word analysis. Brief. Bioinform. 2014, 15, 890–905. [Google Scholar] [CrossRef]
  11. Lu, Y.Y.; Tang, K.; Ren, J.; Fuhrman, J.A.; Waterman, M.S.; Sun, F. CAFE: ACcelerated Alignment-FrEe sequence analysis. Nucl. Acids Res. 2017, 45, W554–W559. [Google Scholar] [CrossRef] [PubMed]
  12. Deng, M.; Yu, C.; Liang, Q.; He, R.L.; Yau, S.S.T. A Novel Method of Characterizing Genetic Sequences: Genome Space with Biological Distance and Applications. PLoS ONE 2011, 6, e17293. [Google Scholar] [CrossRef]
  13. Wen, J.; Chan, R.H.; Yau, S.C.; He, R.L.; Yau, S.S.T. K-mer natural vector and its application to the phylogenetic analysis of genetic sequences. Gene 2014, 546, 25–34. [Google Scholar] [CrossRef] [PubMed]
  14. Yau, S.S.T.; Zhao, X.; Tian, K.; Yu, H. Mathematical Principles in Bioinformatics; Springer: Cham, Switzerland, 2023; pp. 91–144. [Google Scholar]
  15. Sun, N.; Pei, S.; He, L.; Yin, C.; He, R.L.; Yau, S.S.T. Geometric construction of viral genome space and its applications. Comput. Struct. Biotechnol. J. 2021, 19, 4226–4234. [Google Scholar] [CrossRef] [PubMed]
  16. Dong, R.; Pei, S.; Guan, M.; Yau, S.C.; Yin, C.; He, R.L.; Yau, S.S.T. Full Chromosomal Relationships Between Populations and the Origin of Humans. Front. Genet. 2022, 12, 828805. [Google Scholar] [CrossRef]
  17. Yu, H.; Yau, S.S.T. Automated recognition of chromosome fusion using an alignment-free natural vector method. Front. Genet. 2024, 15, 1364951. [Google Scholar] [CrossRef] [PubMed]
  18. Yu, H.; Yau, S.S.T. The optimal metric for viral genome space. Comput. Struct. Biotechnol. J. 2024, 23, 2083–2096. [Google Scholar] [CrossRef] [PubMed]
  19. Cover, T.M.; Hart, P.E. Nearest neighbor pattern classification. IEEE Trans. Inf. Theory 1967, 13, 21–27. [Google Scholar] [CrossRef]
  20. Dekking, F.M.; Kraaikamp, C.; Lopuhaä, H.P.; Meester, L.E. A Modern Introduction to Probability and Statistics: Understanding Why and How; Springer: London, UK, 2005; pp. 377–379. [Google Scholar]
  21. Weglarczyk, S. Kernel density estimation and its application. ITM Web Conf. 2018, 23, 00037. [Google Scholar] [CrossRef]
  22. Taylor, B.S.; Sobieszczyk, M.E.; McCutchan, F.E.; Hammer, S.M. The challenge of HIV-1 subtype diversity. N. Engl. J. Med. 2008, 358, 1590–1602. [Google Scholar] [CrossRef]
  23. D’arc, M.; Ayouba, A.; Esteban, A.; Learn, G.H.; Boué, V.; Liegeois, F.; Etienne, L.; Tagg, N.; Leendertz, F.H.; Boesch, C.; et al. Origin of the HIV-1 group O epidemic in western lowland gorillas. Proc. Natl. Acad. Sci. USA 2015, 112, E1343–E1352. [Google Scholar] [CrossRef]
  24. Mourez, T.; Simon, F.; Plantier, J.C. Non-m variants of human immunodeficiency virus type 1. Clin. Microbiol. Rev. 2013, 26, 448–461. [Google Scholar] [CrossRef] [PubMed]
  25. Plantier, J.C.; Leoz, M.; Dickerson, J.E.; De Oliveira, F.; Cordonnier, F.; Lemée, V.; Damond, F.; Robertson, D.L.; Simon, F. A new human immunodeficiency virus derived from gorillas. Nat. Med. 2009, 15, 871–872. [Google Scholar] [CrossRef] [PubMed]
  26. Hemelaar, J.; Gouws, E.; Ghys, P.D.; Osmanov, S. Global and regional distribution of HIV-1 genetic subtypes and recombinants in 2004. AIDS 2006, 20, W13–W23. [Google Scholar] [CrossRef] [PubMed]
  27. Smith, D.M.; Richman, D.D.; Little, S.J. HIV Superinfection. J. Infect. Dis. 2005, 192, 438–444. [Google Scholar] [CrossRef]
  28. Louten, J. Essential Human Virology, 2nd ed.; Academic Press: Boston, MA, USA, 2023; pp. 277–306. [Google Scholar]
  29. McBride, R.; van Zyl, M.; Fielding, B.C. The coronavirus nucleocapsid is a multifunctional protein. Viruses 2014, 6, 2991–3018. [Google Scholar] [CrossRef] [PubMed]
  30. Carstens, E.B. Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses. Arch. Virol. 2010, 155, 133–146. [Google Scholar] [CrossRef]
  31. Altschul, S.; Gish, W.; Miller, W.; Myers, E.; Lipman, D. Basic Local Aligment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  32. Needleman, S.B.; Wunsch, C.D. A general method applicable to the search for similarities in the amino acid sequence of two proteins. J. Mol. Biol. 1970, 48, 443–453. [Google Scholar] [CrossRef]
Figure 1. Flowchart of Algorithm 1.
Figure 1. Flowchart of Algorithm 1.
Genes 15 00891 g001
Figure 2. Distribution of I i n ( A l p h a ) and I o u t ( · , A l p h a ) (SARS-CoV-2).
Figure 2. Distribution of I i n ( A l p h a ) and I o u t ( · , A l p h a ) (SARS-CoV-2).
Genes 15 00891 g002
Figure 3. Distribution of I i n and I o u t (SARS-CoV-2).
Figure 3. Distribution of I i n and I o u t (SARS-CoV-2).
Genes 15 00891 g003
Figure 4. Scatter plot of D i s R distribution for four types of viruses: A, B, A1B, 03_A6B.
Figure 4. Scatter plot of D i s R distribution for four types of viruses: A, B, A1B, 03_A6B.
Genes 15 00891 g004
Table 1. The optimal weight matrix for viruses.
Table 1. The optimal weight matrix for viruses.
k-mer/Order012
1 3.5 × 10 2 3.9 × 10 2 3.0 × 10 1
2 3.1 × 10 2 1.4 × 10 2 8.5 × 10 2
3 2.1 × 10 1 9.7 × 10 3 3.7 × 10 2
4 4.9 × 10 1 1.8 × 10 3 2.6 × 10 3
5 7.6 × 10 1 1.3 × 10 3 5.4 × 10 3
6 1.0 × 10 0 1.8 × 10 3 1.6 × 10 3
7 8.4 × 10 1 1.1 × 10 3 1.1 × 10 4
8 7.1 × 10 1 1.0 × 10 3 9.3 × 10 4
9 1.8 × 10 1 8.1 × 10 4 5.1 × 10 6
Table 2. Two error rates for Algorithm 1 for three datasets.
Table 2. Two error rates for Algorithm 1 for three datasets.
DatasetType I ErrorType II Error
SARS-CoV-20.94%0.96%
HIV-10.94%0.87%
Orthocoronavirinae0.03%0.98%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, H.; Yau, S.S.-T. New Virus Variant Detection Based on the Optimal Natural Metric. Genes 2024, 15, 891. https://doi.org/10.3390/genes15070891

AMA Style

Yu H, Yau SS-T. New Virus Variant Detection Based on the Optimal Natural Metric. Genes. 2024; 15(7):891. https://doi.org/10.3390/genes15070891

Chicago/Turabian Style

Yu, Hongyu, and Stephen S.-T. Yau. 2024. "New Virus Variant Detection Based on the Optimal Natural Metric" Genes 15, no. 7: 891. https://doi.org/10.3390/genes15070891

APA Style

Yu, H., & Yau, S. S. -T. (2024). New Virus Variant Detection Based on the Optimal Natural Metric. Genes, 15(7), 891. https://doi.org/10.3390/genes15070891

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