Nagy Ádám, Ligeti Balázs, Szebeni János, Pongor Sándor, Gyrffy Balázs
Department of Bioinformatics, Semmelweis University, u 7-9, Tűzoltó, Budapest H-1094, Hungary.
TTK Momentum Cancer Biomarker Research Group, 2, Magyar tudósok körútja, Budapest H-1117, Hungary.
Database (Oxford). 2021 May 8;2021. doi: 10.1093/database/baab020.
Numerous studies demonstrate frequent mutations in the genome of SARS-CoV-2. Our goal was to statistically link mutations to severe disease outcome. We used an automated machine learning approach where 1594 viral genomes with available clinical follow-up data were used as the training set (797 'severe' and 797 'mild'). The best algorithm, based on random forest classification combined with the LASSO feature selection algorithm, was employed to the training set to link mutation signatures and outcome. The performance of the final model was estimated by repeated, stratified, 10-fold cross validation (CV) and then adjusted for multiple testing with Bootstrap Bias Corrected CV. We identified 26 protein and Untranslated Region (UTR) mutations significantly linked to severe outcome. The best classification algorithm uses a mutation signature of 22 mutations as well as the patient's age as the input and shows high classification efficiency with an area under the curve (AUC) of 0.94 [confidence interval (CI): [0.912, 0.962]] and a prediction accuracy of 87% (CI: [0.830, 0.903]). Finally, we established an online platform (https://covidoutcome.com/) that is capable to use a viral sequence and the patient's age as the input and provides a percentage estimation of disease severity. We demonstrate a statistical association between mutation signatures of SARS-CoV-2 and severe outcome of COVID-19. The established analysis platform enables a real-time analysis of new viral genomes.
大量研究表明,严重急性呼吸综合征冠状病毒2(SARS-CoV-2)的基因组中存在频繁突变。我们的目标是从统计学角度将突变与严重疾病结局联系起来。我们采用了一种自动化机器学习方法,将1594个具有可用临床随访数据的病毒基因组用作训练集(797个“重症”和797个“轻症”)。基于随机森林分类与套索特征选择算法相结合的最佳算法应用于训练集,以关联突变特征与结局。通过重复分层10折交叉验证(CV)评估最终模型的性能,然后使用Bootstrap偏差校正CV对多重检验进行调整。我们鉴定出26个与严重结局显著相关的蛋白质和非翻译区(UTR)突变。最佳分类算法使用22个突变的突变特征以及患者年龄作为输入,曲线下面积(AUC)为0.94[置信区间(CI):[0.912,