ABSTRACT
Most algorithms that are used to predict the effects of variants rely on evolutionary conservation. However, a majority of such techniques compute evolutionary conservation by solely using the alignment of multiple sequences while overlooking the evolutionary context of substitution events. We had introduced PHACT, a scoring-based pathogenicity predictor for missense mutations that can leverage phylogenetic trees, in our previous study. By building on this foundation, we now propose PHACTboost, a gradient boosting tree-based classifier that combines PHACT scores with information from multiple sequence alignments, phylogenetic trees, and ancestral reconstruction. By learning from data, PHACTboost outperforms PHACT. Furthermore, the results of comprehensive experiments on carefully constructed sets of variants demonstrated that PHACTboost can outperform 40 prevalent pathogenicity predictors reported in the dbNSFP, including conventional tools, metapredictors, and deep learning-based approaches as well as more recent tools such as AlphaMissense, EVE, and CPT-1. The superiority of PHACTboost over these methods was particularly evident in case of hard variants for which different pathogenicity predictors offered conflicting results. We provide predictions of 215 million amino acid alterations over 20,191 proteins. PHACTboost is available at https://github.com/CompGenomeLab/PHACTboost. PHACTboost can improve our understanding of genetic diseases and facilitate more accurate diagnoses.
Subject(s)
Mutation, Missense , Phylogeny , Humans , Software , Computational Biology/methods , Algorithms , Sequence AlignmentABSTRACT
This paper presents an evaluation of predictions submitted for the "HMBS" challenge, a component of the sixth round of the Critical Assessment of Genome Interpretation held in 2021. The challenge required participants to predict the effects of missense variants of the human HMBS gene on yeast growth. The HMBS enzyme, critical for the biosynthesis of heme in eukaryotic cells, is highly conserved among eukaryotes. Despite the application of a variety of algorithms and methods, the performance of predictors was relatively similar, with Kendall's tau correlation coefficients between predictions and experimental scores around 0.3 for a majority of submissions. Notably, the median correlation (≥ 0.34) observed among these predictors, especially the top predictions from different groups, was greater than the correlation observed between their predictions and the actual experimental results. Most predictors were moderately successful in distinguishing between deleterious and benign variants, as evidenced by an area under the receiver operating characteristic (ROC) curve (AUC) of approximately 0.7 respectively. Compared with the recent two rounds of CAGI competitions, we noticed more predictors outperformed the baseline predictor, which is solely based on the amino acid frequencies. Nevertheless, the overall accuracy of predictions is still far short of positive control, which is derived from experimental scores, indicating the necessity for considerable improvements in the field. The most inaccurately predicted variants in this round were associated with the insertion loop, which is absent in many orthologs, suggesting the predictors still heavily rely on the information from multiple sequence alignment.
ABSTRACT
Ultraviolet (UV) component of solar radiation impairs genome stability by inducing the formation of pyrimidine-pyrimidone (6-4) photoproducts [(6-4)PPs] in plant genomes. (6-4)PPs disrupt growth and development by interfering with transcription and DNA replication. To resist UV stress, plants employ both photoreactivation and nucleotide excision repair that excises oligonucleotide containing (6-4)PPs through two subpathways: global and transcription-coupled excision repair (TCR). Here, we analyzed the genome-wide excision repair-mediated repair of (6-4)PPs in Arabidopsis thaliana and found that (6-4)PPs can be repaired by TCR; however, the main subpathway to remove (6-4)PPs from the genome is global repair. Our analysis showed that open chromatin genome regions are more rapidly repaired than heterochromatin regions, and the repair level peaks at the promoter, transcription start site and transcription end site of genes. Our study revealed that the repair of (6-4)PP in plants showed a distinct genome-wide repair profile compared to the repair of other major UV-induced DNA lesion called cyclobutane pyrimidine dimers (CPDs).