Utilizing genomic signatures to gain insights into the dynamics of SARS-CoV-2 through Machine and Deep Learning techniques

被引:0
作者
Ahmed M. A. Elsherbini
Amr Hassan Elkholy
Youssef M. Fadel
Gleb Goussarov
Ahmed Mohamed Elshal
Mohamed El-Hadidi
Mohamed Mysara
机构
[1] Nile University,Bioinformatics Group, Center for Informatics Science, School of Information Technology and Computer Science
[2] Belgian Nuclear Research Centre (SCK•CEN),Microbiology Unit
来源
BMC Bioinformatics | / 25卷
关键词
SARS-CoV-2; Genomic signature; Di nucleotide frequency; Tri nucleotide frequency; GenoSig; Deep Learning; Machine Learning; Random Forest;
D O I
暂无
中图分类号
学科分类号
摘要
The global spread of the SARS-CoV-2 pandemic, originating in Wuhan, China, has had profound consequences on both health and the economy. Traditional alignment-based phylogenetic tree methods for tracking epidemic dynamics demand substantial computational power due to the growing number of sequenced strains. Consequently, there is a pressing need for an alignment-free approach to characterize these strains and monitor the dynamics of various variants. In this work, we introduce a swift and straightforward tool named GenoSig, implemented in C++. The tool exploits the Di and Tri nucleotide frequency signatures to delineate the taxonomic lineages of SARS-CoV-2 by employing diverse machine learning (ML) and deep learning (DL) models. Our approach achieved a tenfold cross-validation accuracy of 87.88% (± 0.013) for DL and 86.37% (± 0.0009) for Random Forest (RF) model, surpassing the performance of other ML models. Validation using an additional unexposed dataset yielded comparable results. Despite variations in architectures between DL and RF, it was observed that later clades, specifically GRA, GRY, and GK, exhibited superior performance compared to earlier clades G and GH. As for the continental origin of the virus, both DL and RF models exhibited lower performance than in predicting clades. However, both models demonstrated relatively higher accuracy for Europe, North America, and South America compared to other continents, with DL outperforming RF. Both models consistently demonstrated a preference for cytosine and guanine over adenine and thymine in both clade and continental analyses, in both Di and Tri nucleotide frequencies signatures. Our findings suggest that GenoSig provides a straightforward approach to address taxonomic, epidemiological, and biological inquiries, utilizing a reductive method applicable not only to SARS-CoV-2 but also to similar research questions in an alignment-free context.
引用
收藏
相关论文
共 183 条
[31]  
Giorgi FM(2008)Scikit-learn: machine learning in Python J Stat Softw 6 90-undefined
[32]  
Hamed SM(2014)Building predictive models in R using the caret package Bioinformatics 585 1474-undefined
[33]  
Elkhatib WF(2011)circlize implements and enhances circular visualization in R Biometrics 9 764-undefined
[34]  
Khairalla AS(2020)ggplot2: elegant graphics for data analysis by Wickham, H Nature 50 66-undefined
[35]  
Noreddin AM(2007)Array programming with NumPy Comput Sci Eng 27 66-undefined
[36]  
Zhou Y(2018)Matplotlib: a 2D graphics environment Nat Genet 6 93-undefined
[37]  
Zhi H(2011)Functional classification of long non-coding RNAs by k-mer content Bioinformatics 6 186-undefined
[38]  
Teng Y(2022)A fast, lock-free approach for efficient parallel counting of occurrences of k-mers bioRxiv 1 1494-undefined
[39]  
Randhawa GS(2018)MerCat2: a versatile k-mer counter and diversity estimator for database-independent property analysis obtained from omics data Bioinformatics 18 66-undefined
[40]  
Soltysiak MPM(2018)NextStrain: real-time tracking of pathogen evolution Annu Rev Biomed Data Sci 94 15-undefined