Skip to main content

3DVizSNP: a tool for rapidly visualizing missense mutations identified in high throughput experiments in iCn3D

Abstract

Background

High throughput experiments in cancer and other areas of genomic research identify large numbers of sequence variants that need to be evaluated for phenotypic impact. While many tools exist to score the likely impact of single nucleotide polymorphisms (SNPs) based on sequence alone, the three-dimensional structural environment is essential for understanding the biological impact of a nonsynonymous mutation.

Results

We present a program, 3DVizSNP, that enables the rapid visualization of nonsynonymous missense mutations extracted from a variant caller format file using the web-based iCn3D visualization platform. The program, written in Python, leverages REST APIs and can be run locally without installing any other software or databases, or from a webserver hosted by the National Cancer Institute. It automatically selects the appropriate experimental structure from the Protein Data Bank, if available, or the predicted structure from the AlphaFold database, enabling users to rapidly screen SNPs based on their local structural environment. 3DVizSNP leverages iCn3D annotations and its structural analysis functions to assess changes in structural contacts associated with mutations.

Conclusions

This tool enables researchers to efficiently make use of 3D structural information to prioritize mutations for further computational and experimental impact assessment. The program is available as a webserver at https://analysistools.cancer.gov/3dvizsnp or as a standalone python program at https://github.com/CBIIT-CGBB/3DVizSNP.

Peer Review reports

Background

High-throughput sequencing experiments generate large numbers of variant calls. Though there are many types of variants, here we focus on nonsynonymous missense mutations affecting protein coding regions. There are databases with annotations of known SNPs, such as dbSNP [1], ClinVar [2], and COSMIC [3], but the volume of data generated by ongoing experiments means that researchers need tools to evaluate the effects of SNPs that, for most, are not present in these databases. The computational challenge is in predicting which variants may impact biological function [4,5,6,7,8,9]. There are numerous tools available that predict the possible deleterious impact of a mutation, most of which use one-dimensional sequence only [8, 10, 11]. The three-dimensional structural context of variants provides critical information for assessing the potential impact of an amino acid mutation. While the number of experimentally solved 3D structures is still relatively small compared to the number of sequences, the success of deep-learning based structure predictors such as AlphaFold [12] and others [13, 14] have greatly increased the number of proteins with high-quality structural information. The challenge is to make use of all this information in efficient ways. It is typically cumbersome and time consuming for non-structural biologists to identify the protein where a given mutation occurs, determine if the mutation is in a known PDB structure or not, load the PDB structure (if available) or a predicted model, bring up the mutation in a 3D viewer where one can visualize the local structural context for the mutation, and analyze differences due to a mutation in structural terms. Some tools make use of 3D information to predict the effect of a mutation, such as Missense3D [15], HOPE [16], or StructMAn [17], but these tools are generally not set up to operate efficiently on a large number of mutations. They also may require more advanced structural biology and biophysical knowledge to operate and/or interpret the results, while non-structural biologists need to quickly and efficiently make use of the available 3D structural information to help assess the likely impact of a nonsynonymous mutation. Below we describe a tool called 3DVizSNP designed to fill this gap (Fig. 1).

Fig. 1
figure 1

Illustration depicting the role 3DVizSNP fills in the process of evaluating missense mutations. Many programs provide impact predictions for large numbers of SNPs based on sequence (purple box), and there are algorithms that assess the impact of individual mutations using 3D structure (blue box), but 3DVizSNP allows the user to easily visualize a moderate number of mutations in 3D to facilitate prioritization of mutations for further study

Implementation

3DVizSNP is a Python program that, combined with the iCn3D structure viewing and analysis platform [18, 19], allows the user to quickly process a VCF file and produce a table with Ensembl Gene IDs, Gene symbols, UniProt IDs, PDB IDs (if applicable), SIFT [20] and PolyPhen [21] predictions and scores, and a link that can open iCn3D highlighting the mutations mapped in 3D on either the relevant PDB structure or the AlphaFold predicted structure. The user can quickly go through tens to hundreds of nonsynonymous mutations identified in high-throughput experiments and visualize the structure and structural contacts to inform and prioritize further studies, without requiring extensive structural biology expertise or heavy computational resources, and without installing a standalone 3D structure viewer or modeling software. It leverages REST APIs including that of iCn3D, which enables the rapid retrieval of experimental structures from the PDB or predictions from the AlphaFold database, and display of contacts in 1D/2D/3D formats, all through a saveable and shareable URL [18].

3DVizSNP workflow

3DVizSNP reads in a bgzipped VCF file with an associated tabix index and extracts the variants, with the option to only select variants from an input list of HGNC gene symbols. It then submits the variants to the VEP server [11] via a REST API and gets back the Ensembl Gene IDs, Gene symbols, SwissProt IDs, amino acid mutation, and SIFT and PolyPhen predictions and scores. It uses the PDB REST API to identify the PDB ID of the X-ray crystal or electron micrograph with the highest resolution, or the NMR structure in the absence of X-ray or EM structures, for each amino acid mutation (correcting for mismatches between UniProt and PDB numbering). It also avoids structures with engineered mutations at the mutation site. If an experimental structure is not available, it loads the AlphaFold prediction (for sequences under 2700 amino acids in length) color-coded by prediction confidence. The output is an HTML table and a comma separated value (csv) file, which includes iCn3D links for each mutation. (Fig. 2) The links can be clicked to open the PDB structure or AlphaFold predicted structure in iCn3D, with the mutations added as separate tracks in the 1D viewer for SIFT and PolyPhen predictions. The iCn3D ‘scap interactions’ command is used to highlight the wild-type and mutant amino acid side chains (which can be toggled) in iCn3D. iCn3D has the capability to add multiple annotations to the sequence viewer, enabling the user to quickly determine if the variant is in a known functional domain, matches known mutations in dbSNP or ClinVar, etc. The csv file can be loaded into a spreadsheet program, enabling the user to record notes about variants, sort them, store the results, etc.

Fig. 2
figure 2

The 3DVizSNP Workflow

In addition to the python script, we have provided a website where the user can upload a VCF file and view the sortable, filterable output table along with an embedded version of iCn3D (Fig. 3). Clicking on a row brings up the mutation in iCn3D above the table, while clicking on the iCn3D button will open the mutation in a separate window with the full version of iCn3D.

Fig. 3
figure 3

The web interface to 3DVizSNP. The current row is highlighted in green and shown in the embedded iCn3D interface. Clicking on a different row loads that mutation in the iCn3D viewer. The table is sortable and filterable, and clicking on the ‘View in iCn3D’ button opens the mutation in a new full page iCn3D

Results

Known KRAS mutations

To illustrate the ability to find and visualize known mutations in the KRAS oncogene, we ran 3DVizSNP against a colorectal adenocarcinoma (COAD) vcf file from a single patient from the TCGA database with the flag ‘-g KRAS’ to select only variants in the KRAS gene. This produced an output table with two mutations, the well-known G12V mutation and a 2nd mutation, D154N. Figure 4 depicts these mutations in iCn3D. The G12V mutation is predicted to be ‘deleterious_low_confidence’ by SIFT (with a score of 0) and is not called in this instance by PolyPhen. (Note that the PDB ID 7rov has the lowest resolution structure, so it is used for the D154N mutation, but it is an engineered G12D mutant, so 3DVizSNP discards it in favor of 7vvb for the G12V mutation.) The impact of the G12V mutant is not obvious upon initial inspection, but it is adjacent to the GTP binding site and involves replacing a glycine in a turn with a valine, which is likely to cause steric hindrances. (Fig. 4A) The D154N mutation is also predicted to be ‘deleterious_low_confidence’ by SIFT with a score of 0 and is not called by PolyPhen. The mutation is subtle, only appearing to disrupt a single hydrogen bond to a water molecule (Fig. 4B). However, KRAS activity depends upon dimerization and D154 helps stabilize the alpha-alpha form of the KRAS dimer [22]. It has been reported that a D154Q mutation affects dimerization and has a growth inhibitory effect on oncogenic versions of KRAS, [23] though the impact of mutating this residue on dimerization is disputed [24].

Fig. 4
figure 4

Mutations in the KRAS protein. A The G12V mutation (WT, left; G12V, right). The SIFT score is zero and the prediction is ‘deleterious_low_confidence’. The active site magnesium ion (green sphere) and nonhydrolyzable ATP analog GNP are seen near the mutated residue. B The D154N mutation (WT, left; D154N, right). The SIFT score is zero and the prediction is ‘deleterious_low_confidence’. The red arrows highlight the hydrogen bond to a water molecule that is removed by the mutation

TCGA VCF example

To illustrate the utility of 3DVizSNP in screening mutations of unknown effect, we took a randomly selected VCF file from a Glioblastoma patient in the TCGA database, which had 27,526 variants in it. We submitted it to the OpenCRAVAT server [10] to filter the list down to a more manageable size, by using the filters coding = yes, polyphen-2 HumVar rank score 0.8, and SIFT prediction ‘damaging’. This produced a list of 901 variants, which we ran through 3DVizSNP. We looked at each variant in iCn3D and investigated variants that met certain criteria for further inspection. These criteria included being in a PDB structure or in a Confident or Very High region of an AlphaFold prediction, not being on the surface of the structure, and having significant changes in the structure, such as changes in size/polarity/charge, changes in contacts, etc. Low confidence areas of AlphaFold predictions are typically not well-ordered and little inference can be drawn regarding the effects of mutations in those regions. We avoided surface mutations because it is harder to predict what effects they might have on protein structure and function. Obviously, they can affect binding to other proteins but in the absence of knowledge about these interfaces one cannot make quick determinations of the impact of mutations on binding. We then looked at the UniProt page for information on the protein function, known mutations, etc. We also submitted the mutation to the Missense3D server [15] and the SAAMBE server [25] to obtain structure-based predictions of mutational impact. The 44th mutation in the output list was a L194R mutation in UniProt ID Q5TCH4, a Cytochrome P450 4A protein. (Fig. 5) The Missense3D server predicts it causes structural damage, as it introduces a buried hydrophilic and charged residue for a buried hydrophobic residue. The SAAMBE server predicts the mutation is disruptive with a ΔΔG of 0.31 kcal/mol (that is, the free energy of folding of the mutant protein is predicted to be 0.31 kcal/mol higher than that of the wild-type). The UniProt entry has a link to the DisGenNET database (entry 284541) which links mutations to diseases. DisGenNET lists liver carcinoma as a disease linked to CYP4A22, and links to PubMed ID 30069903 [26], which indicates that higher CYP4A22 expression in hepatocellular carcinoma correlates with better disease prognosis. Obviously hepatocellular carcinoma differs from glioblastoma; the point is not to make a strong claim that this particular mutation contributes to glioblastoma progression, but rather to point out that there is evidence linking mutations that damage CYP4A22 function, as the L194R mutation is predicted to, to cancer progression, indicating further study is warranted. This process illustrates that one can use tools such as OpenCRAVAT to filter a larger set of variants, then use 3DVizSNP to further sort and prioritize variants based on the available 3D structural information. One can then use other tools and databases such as Missense3D, SAAMBE, and UniProt to derive further support for the potential functional relevance of a given mutation. Continuing this process could winnow down the list of variants to around 50 for further in-depth computational assessment and experimental characterization.

Fig. 5
figure 5

A closeup view of the wildtype leucine 194 (gray) and the mutation to arginine (red) in the AlphaFold prediction for UniProt ID Q5TCH4, a Cytochrome P450 4A protein. The mutation was identified in a vcf file from TCGA (see text for details). The blue color of the other residues indicates that this is a very high confidence region for the AlphaFold structure prediction

Discussion

There are other tools that allow one to visualize missense mutations in a 3D context. (Additional file 1) PhyreRisk from the Sternberg lab does similar things to 3DVizSNP [27]. It accepts more formats for submitting variants than 3DVizSNP, also uses the VEP server, displays the list of variants in a table, and provides an interactive 3D interface for visualizing the mutation. The main advantages of 3DVizSNP are the ability to toggle between the mutant and wild-type residues, the depiction of intramolecular contacts for the wild-type and mutant, the use of Alphafold models instead of Phyre models, and the fact that the structure viewer is embedded in the webpage with the output table. Phyre has been used for many years to successfully model protein structures, but it does not have the same accuracy or breadth of coverage that AlphaFold2 does [28]. MuPit from the Karchin lab also displays mutations on 3D structures along with annotations [29, 30]. Like PhyreRisk, it does not allow toggling between the mutant and wild-type residue and does not show intramolecular contacts, and it is limited to mapping mutations onto PDB structures only. VIVID is a recently developed tool designed to accomplish many of the same tasks as 3DVizSNP [31]. It requires the user to input a gene sequence, a VCF file with mutations, a GFF file, and a PDB file (which can be obtained from an Alphafold prediction). It allows visualization of nonsynonymous mutations on a 3D structure, along with a variety of metrics including both population level information and mutational impact calculations. It appears to be limited to evaluating one protein structure at a time, and does not allow toggling between the wild-type and mutant sidechains.

3DVizSNP provides complementary features to these existing powerful tools. Generally speaking, these other tools are designed from the perspective of viewing multiple mutations on a particular protein structure, whereas 3DVizSNP is designed to from the perspective of viewing mutations from a particular VCF file. The existing tools provide less ability to screen large numbers of mutations found in a VCF file than 3DVizSNP while providing more in-depth information about each mutation or protein.

While 3DVizSNP provides significant novel capabilities, it does have limitations. Due to querying REST APIs it is not particularly fast. It takes 3 min, 21 s to produce output for 1000 variants on a MacBook Pro laptop with 2.3 GHz 8-Core Intel Core i9 processor and can take closer to 30 min for longer input VCF files. Nonetheless the running time is small compared to the time required to visually assess tens to hundreds of mutations, and is comparable to the running times of the servers discussed above. 3DVizSNP requires filtering of VCF files (which may have millions of variants) in advance using tools like vcftools [32], bcftools [33], or OpenCRAVAT [10]), as it is not practical to use with more than approximately 1000 variants in the output.

Conclusions

Here we present a tool, 3DVizSNP, that facilitates the rapid assessment of phenotypic impact of missense mutations extracted from VCF files by enabling the user to visualize multiple mutations quickly in the iCn3D web-based structure and sequence analysis platform. This tool provides structural context not present in sequence-based impact prediction programs but enables the viewing of multiple mutations more rapidly and easily than other 3D impact assessment programs. This tool will make it easier for researchers to prioritize mutations for further study, a critical bottleneck in modern high-throughput experiments.

Availability and requirements

Project name: 3DVizSNP. Project home page: https://analysistools.cancer.gov/3dvizsnp. (Source code: https://github.com/CBIIT-CGBB/3DVizSNP). Operating system(s): Platform independent. Programming language: Python. Other requirements: Pysam v. 0.19.1, Pandas library. License: MIT License.

Availability of data and materials

The python code presented here is publicly available at https://github.com/CBIIT-CGBB/3DVizSNP. The webserver is publicly available at https://analysistools.cancer.gov/3dvizsnp. Example VCF files used are publicly available at the Github, at TCGA, or from the authors upon request.

Abbreviations

SNP:

Single nucleotide polymorphism

PDB:

Protein data bank

VCF:

Variant caller format

REST:

Representational state transfer

API:

Application programming interface

References

  1. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Landrum MJ, Lee JM, Benson M, Brown G, Chao C, Chitipiralla S, Gu B, Hart J, Hoffman D, Hoover J, et al. ClinVar: public archive of interpretations of clinically relevant variants. Nucleic Acids Res. 2016;44(D1):D862-868.

    Article  CAS  PubMed  Google Scholar 

  3. Bamford S, Dawson E, Forbes S, Clements J, Pettett R, Dogan A, Flanagan A, Teague J, Futreal PA, Stratton MR, et al. The COSMIC (catalogue of somatic mutations in cancer) database and website. Br J Cancer. 2004;91(2):355–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. George Priya Doss C, Sudandiradoss C, Rajasekaran R, Choudhury P, Sinha P, Hota P, Batra UP, Rao S. Applications of computational algorithm tools to identify functional SNPs. Funct Integr Genom. 2008;8(4):309–16.

    Article  CAS  Google Scholar 

  5. Espinosa O, Mitsopoulos K, Hakas J, Pearl F, Zvelebil M. Deriving a mutation index of carcinogenicity using protein structure and protein interfaces. PLoS ONE. 2014;9(1):e84598.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Gerasimavicius L, Livesey BJ, Marsh JA. Loss-of-function, gain-of-function and dominant-negative mutations have profoundly different effects on protein structure. Nat Commun. 2022;13(1):3895.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Gong S, Worth CL, Cheng TM, Blundell TL. Meet me halfway: when genomics meets structural bioinformatics. J Cardiovasc Transl Res. 2011;4(3):281–303.

    Article  PubMed  Google Scholar 

  8. Hu Z, Yu C, Furutsuki M, Andreoletti G, Ly M, Hoskins R, Adhikari AN, Brenner SE. VIPdb, a genetic variant impact predictor database. Hum Mutat. 2019;40(9):1202–14.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Glusman G, Rose PW, Prlic A, Dougherty J, Duarte JM, Hoffman AS, Barton GJ, Bendixen E, Bergquist T, Bock C, et al. Mapping genetic variations to three-dimensional protein structures to enhance variant interpretation: a proposed framework. Genome Med. 2017;9(1):113.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Pagel KA, Kim R, Moad K, Busby B, Zheng L, Tokheim C, Ryan M, Karchin R. Integrated informatics analysis of cancer-related variants. JCO Clin Cancer Inform. 2020;4:310–7.

    Article  PubMed  Google Scholar 

  11. McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, Flicek P, Cunningham F. The Ensembl variant effect predictor. Genome Biol. 2016;17(1):122.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Varadi M, Anyango S, Deshpande M, Nair S, Natassia C, Yordanova G, Yuan D, Stroe O, Wood G, Laydon A, et al. AlphaFold Protein structure database: massively expanding the structural coverage of protein-sequence space with high-accuracy models. Nucleic Acids Res. 2022;50(D1):D439–44.

    Article  CAS  PubMed  Google Scholar 

  13. Chowdhury R, Bouatta N, Biswas S, Floristean C, Kharkar A, Roy K, Rochereau C, Ahdritz G, Zhang J, Church GM, et al. Single-sequence protein structure prediction using a language model and deep learning. Nat Biotechnol. 2022;40(11):1617–23.

    Article  CAS  PubMed  Google Scholar 

  14. Lin Z, Akin H, Rao R, Hie B, Zhu Z, Lu W, Smetanin N, Verkuil R, Kabeli O, Shmueli Y et al. Evolutionary-scale prediction of atomic level protein structure with a language model. Science. 2023;379(6637):1123–30.

    Article  Google Scholar 

  15. Ittisoponpisan S, Islam SA, Khanna T, Alhuzimi E, David A, Sternberg MJE. Can predicted protein 3D structures provide reliable insights into whether missense variants are disease associated? J Mol Biol. 2019;431(11):2197–212.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Venselaar H, Te Beek TA, Kuipers RK, Hekkelman ML, Vriend G. Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces. BMC Bioinf. 2010;11:548.

    Article  Google Scholar 

  17. Gress A, Ramensky V, Buch J, Keller A, Kalinina OV. StructMAn: annotation of single-nucleotide polymorphisms in the structural context. Nucleic Acids Res. 2016;44(W1):W463-468.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wang J, Youkharibache P, Marchler-Bauer A, Lanczycki C, Zhang D, Lu S, Madej T, Marchler GH, Cheng T, Chong LC, et al. iCn3D: from web-based 3D viewer to structural analysis tool in batch mode. Front Mol Biosci. 2022;9:831740.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wang J, Youkharibache P, Zhang D, Lanczycki CJ, Geer RC, Madej T, Phan L, Ward M, Lu S, Marchler GH, et al. iCn3D, a web-based 3D viewer for sharing 1D/2D/3D representations of biomolecular structures. Bioinformatics. 2020;36(1):131–5.

    Article  CAS  PubMed  Google Scholar 

  20. Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31(13):3812–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;7:20.

    Google Scholar 

  22. Lee KY, Enomoto M, Gebregiworgis T, Gasmi-Seabrook GMC, Ikura M, Marshall CB. Oncogenic KRAS G12D mutation promotes dimerization through a second, phosphatidylserine-dependent interface: a model for KRAS oligomerization. Chem Sci. 2021;12(38):12827–37.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ambrogio C, Kohler J, Zhou ZW, Wang H, Paranal R, Li J, Capelletti M, Caffarra C, Li S, Lv Q, et al. KRAS dimerization impacts MEK inhibitor sensitivity and oncogenic activity of mutant KRAS. Cell. 2018;172(4):857–68.

    Article  CAS  PubMed  Google Scholar 

  24. Grozavu I, Stuart S, Lyakisheva A, Yao Z, Pathmanathan S, Ohh M, Stagljar I. D154Q mutation does not alter KRAS dimerization. J Mol Biol. 2022;434(2):167392.

    Article  CAS  PubMed  Google Scholar 

  25. Pahari S, Li G, Murthy AK, Liang S, Fragoza R, Yu H, Alexov E. SAAMBE-3D: predicting effect of mutations on protein-protein interactions. Int J Mol Sci. 2020;21:7.

    Article  Google Scholar 

  26. Eun HS, Cho SY, Lee BS, Kim S, Song IS, Chun K, Oh CH, Yeo MK, Kim SH, Kim KH. Cytochrome P450 4A11 expression in tumor cells: a favorable prognostic factor for hepatocellular carcinoma patients. J Gastroenterol Hepatol. 2019;34(1):224–33.

    Article  CAS  PubMed  Google Scholar 

  27. Ofoegbu TC, David A, Kelley LA, Mezulis S, Islam SA, Mersmann SF, Stromich L, Vakser IA, Houlston RS, Sternberg MJE. PhyreRisk: a dynamic web application to bridge genomics, proteomics and 3D structural data to guide interpretation of human genetic variants. J Mol Biol. 2019;431(13):2460–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kryshtafovych A, Schwede T, Topf M, Fidelis K, Moult J. Critical assessment of methods of protein structure prediction (CASP)-Round XIV. Proteins. 2021;89(12):1607–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Niknafs N, Kim D, Kim R, Diekhans M, Ryan M, Stenson PD, Cooper DN, Karchin R. MuPIT interactive: webserver for mapping variant positions to annotated, interactive 3D structures. Hum Genet. 2013;132(11):1235–43.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Tokheim C, Bhattacharya R, Niknafs N, Gygax DM, Kim R, Ryan M, Masica DL, Karchin R. Exome-scale discovery of hotspot mutation regions in human cancer using 3D protein structure. Cancer Res. 2016;76(13):3719–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Tichkule S, Myung Y, Naung MT, Ansell BRE, Guy AJ, Srivastava N, Mehra S, Caccio SM, Mueller I, Barry AE, et al. VIVID: a web application for variant interpretation and visualization in multi-dimensional analyses. Mol Biol Evol. 2022;39:9.

    Article  Google Scholar 

  32. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whitwham A, Keane T, McCarthy SA, Davies RM, et al. Twelve years of SAMtools and BCFtools. Gigascience. 2021;10(2):giab008.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This project was initiated at a Hackathon organized by two of the authors (JW and PY) held at ISMB2022. MS was a team leader for the project during the Hackathon and MW was a team member and set up the initial GitHub repository for the project. The authors would like to thank the other team members Bonface Onyango, and Pranavathiyani G for their help in getting the project started. The authors also would like to thank the other members of the CGBB for testing the script and the website, and for providing valuable feedback.

Funding

Open Access funding provided by the National Institutes of Health (NIH).

Author information

Authors and Affiliations

Authors

Contributions

MS, DM, JW, and PY conceived of the project. MW, SR, and MS wrote the python code. BC and BP built the webserver. MS wrote the manuscript. MW, MS, and PY edited the manuscript. All authors read and approved of the final manuscript.

Corresponding author

Correspondence to Michael Sierk.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

A table summarizing the features of 3DVizSNP, PhyreRisk, MuPit, and VIVID.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sierk, M., Ratnayake, S., Wagle, M.M. et al. 3DVizSNP: a tool for rapidly visualizing missense mutations identified in high throughput experiments in iCn3D. BMC Bioinformatics 24, 244 (2023). https://doi.org/10.1186/s12859-023-05370-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12859-023-05370-5

Keywords