Mutalyzer 2: next generation HGVS nomenclature checker.

MOTIVATION: Unambiguous variant descriptions are of utmost importance in clinical genetic diagnostics, scientific literature and genetic databases. The Human Genome Variation Society (HGVS) publishes a comprehensive set of guidelines on how variants should be correctly and unambiguously described. We present the implementation of the Mutalyzer 2 tool suite, designed to automatically apply the HGVS guidelines so users do not have to deal with the HGVS intricacies explicitly to check and correct their variant descriptions. RESULTS: Mutalyzer is profusely used by the community, having processed over 133 million descriptions since its launch. Over a five year period, Mutalyzer reported a correct input in ∼50% of cases. In 41% of the cases either a syntactic or semantic error was identified and for ∼7% of cases, Mutalyzer was able to automatically correct the description. AVAILABILITY AND IMPLEMENTATION: Mutalyzer is an Open Source project under the GNU Affero General Public License. The source code is available on GitHub (https://github.com/mutalyzer/mutalyzer) and a running instance is available at: https://mutalyzer.nl.

Synonymous Somatic Variants in Human Cancer Are Not Infamous: A Plea for Full Disclosure in Databases and Publications.

Single-nucleotide variants (SNVs) are the most frequent genetic changes found in human cancer. Most driver alterations are missense and nonsense variants localized in the coding region of cancer genes. Unbiased cancer genome sequencing shows that synonymous SNVs (sSNVs) can be found clustered in the coding regions of several cancer oncogenes or tumor suppressor genes suggesting purifying selection. sSNVs are currently underestimated, as they are usually discarded during analysis. Furthermore, several public databases do not display sSNVs, which can lead to analytical bias and the false assumption that this mutational event is uncommon. Recent progress in our understanding of the deleterious consequences of these sSNVs for RNA stability and protein translation shows that they can act as strong drivers of cancer, as demonstrated for several cancer genes such as TP53 or BCL2L12. It is therefore essential that sSNVs be properly reported and analyzed in order to provide an accurate picture of the genetic landscape of the cancer genome.

HGVS Recommendations for the Description of Sequence Variants: 2016 Update.

The consistent and unambiguous description of sequence variants is essential to report and exchange information on the analysis of a genome. In particular, DNA diagnostics critically depends on accurate and standardized description and sharing of the variants detected. The sequence variant nomenclature system proposed in 2000 by the Human Genome Variation Society has been widely adopted and has developed into an internationally accepted standard. The recommendations are currently commissioned through a Sequence Variant Description Working Group (SVD-WG) operating under the auspices of three international organizations: the Human Genome Variation Society (HGVS), the Human Variome Project (HVP), and the Human Genome Organization (HUGO). Requests for modifications and extensions go through the SVD-WG following a standard procedure including a community consultation step. Version numbers are assigned to the nomenclature system to allow users to specify the version used in their variant descriptions. Here, we present the current recommendations, HGVS version 15.11, and briefly summarize the changes that were made since the 2000 publication. Most focus has been on removing inconsistencies and tightening definitions allowing automatic data processing. An extensive version of the recommendations is available online, at http://www.HGVS.org/varnomen.

An efficient algorithm for the extraction of HGVS variant descriptions from sequences.

MOTIVATION: Unambiguous sequence variant descriptions are important in reporting the outcome of clinical diagnostic DNA tests. The standard nomenclature of the Human Genome Variation Society (HGVS) describes the observed variant sequence relative to a given reference sequence. We propose an efficient algorithm for the extraction of HGVS descriptions from two sequences with three main requirements in mind: minimizing the length of the resulting descriptions, minimizing the computation time and keeping the unambiguous descriptions biologically meaningful. RESULTS: Our algorithm is able to compute the HGVS descriptions of complete chromosomes or other large DNA strings in a reasonable amount of computation time and its resulting descriptions are relatively small. Additional applications include updating of gene variant database contents and reference sequence liftovers. AVAILABILITY: The algorithm is accessible as an experimental service in the Mutalyzer program suite (https://mutalyzer.nl). The C++ source code and Python interface are accessible at: https://github.com/mutalyzer/description-extractor. CONTACT: j.k.vis@lumc.nl.

A germline chromothripsis event stably segregating in 11 individuals through three generations.

PURPOSE: Parentally transmitted germ-line chromothripsis (G-CTH) has been identified in only a few cases. Most of these rearrangements were stably transmitted, in an unbalanced form, from a healthy mother to her child with congenital abnormalities probably caused by de novo copy-number changes of dosage sensitive genes. We describe a G-CTH transmitted through three generations in 11 healthy carriers. METHODS: Conventional cytogenetic analysis, mate-pair sequencing, and polymerase chain reaction (PCR) were used to identify the chromosome rearrangement and characterize the breakpoints in all three generations. RESULTS: We identified an apparently balanced translocation t(3;5), later shown to be a G-CTH, in all individuals of a three-generation family. The G-CTH stably segregated without occurrence of additional rearrangements; however, several spontaneous abortions were reported, possibly due to unbalanced transmission. Although seven protein-coding genes are interrupted, no clinical features can be definitively attributed to the affected genes. However, it can be speculated that truncation of one of these genes, encoding ataxia-telangiectasia and Rad3-related protein kinase (ATR), a key component of the DNA damage response, may be related to G-CTH formation. CONCLUSION: G-CTH rearrangements are not always associated with abnormal phenotypes and may be misinterpreted as balanced two-way translocations, suggesting that G-CTH is an underdiagnosed phenomenon.Genet Med 18 5, 494-500.

Recommendations for analyzing and reporting TP53 gene variants in the high-throughput sequencing era.

The architecture of TP53, the most frequently mutated gene in human cancer, is more complex than previously thought. Using TP53 variants as clinical biomarkers to predict response to treatment or patient outcome requires an unequivocal and standardized procedure toward a definitive strategy for the clinical evaluation of variants to provide maximum diagnostic sensitivity and specificity. An intronic promoter and two novel exons have been identified resulting in the expression of multiple transcripts and protein isoforms. These regions are additional targets for mutation events impairing the tumor suppressive activity of TP53. Reassessment of variants located in these regions is needed to refine their prognostic value in many malignancies. We recommend using the stable Locus Reference Genomic reference sequence for detailed and unequivocal reports and annotations of germ line and somatic alterations on all TP53 transcripts and protein isoforms according to the recommendations of the Human Genome Variation Society. This novel and comprehensive description framework will generate standardized data that are easy to understand, analyze, and exchange across various cancer variant databases. Based on the statistical analysis of more than 45,000 variants in the latest version of the UMD TP53 database, we also provide a classification of their functional effects ("pathogenicity").

Curating gene variant databases (LSDBs): toward a universal standard.

Gene variant databases or Locus-Specific DataBases (LSDBs) are used to collect and display information on sequence variants on a gene-by-gene basis. Their most frequent use is in relation to DNA-based diagnostics, giving clinicians and scientists easy access to an up-to-date overview of all gene variants identified worldwide and whether they influence the function of the gene ("pathogenic or not"). While literature on gene variant databases is extensive, little has been published on the process of database curation itself. Based on our extensive experience as LSDB curators and our contributions to database curation courses, we discuss the subject of database curation. We describe the tasks involved, the steps to take, and the issues that might occur. Our overview is a first step toward establishing overall guidelines for database curation and ultimately covers one aspect of establishing quality-assured gene variant databases.

A formalized description of the standard human variant nomenclature in Extended Backus-Naur Form.

BACKGROUND: The use of a standard human sequence variant nomenclature is advocated by the Human Genome Variation Society in order to unambiguously describe genetic variants in databases and literature. There is a clear need for tools that allow the mining of data about human sequence variants and their functional consequences from databases and literature. Existing text mining focuses on the recognition of protein variants and their effects. The recognition of variants at the DNA and RNA levels is essential for dissemination of variant data for diagnostic purposes. Development of new tools is hampered by the complexity of the current nomenclature, which requires processing at the character level to recognize the specific syntactic constructs used in variant descriptions. RESULTS: We approached the gene variant nomenclature as a scientific sublanguage and created two formal descriptions of the syntax in Extended Backus-Naur Form: one at the DNA-RNA level and one at the protein level. To ensure compatibility to older versions of the human sequence variant nomenclature, previously recommended variant description formats have been included. The first grammar versions were designed to help build variant description handling in the Alamut mutation interpretation software. The DNA and RNA level descriptions were then updated and used to construct the context-free parser of the Mutalyzer 2 sequence variant nomenclature checker, which has already been used to check more than one million variant descriptions. CONCLUSIONS: The Extended Backus-Naur Form provided an overview of the full complexity of the syntax of the sequence variant nomenclature, which remained hidden in the textual format and the division of the recommendations across the DNA, RNA and protein sections of the Human Genome Variation Society nomenclature website (http://www.hgvs.org/mutnomen/). This insight into the syntax of the nomenclature could be used to design detailed and clear rules for software development. The Mutalyzer 2 parser demonstrated that it facilitated decomposition of complex variant descriptions into their individual parts. The Extended Backus-Naur Form or parts of it can be used or modified by adding rules, allowing the development of specific sequence variant text mining tools and other programs, which can generate or handle sequence variant descriptions.

Describing structural changes by extending HGVS sequence variation nomenclature.

New technologies allow rapid discovery of novel sequence variants among which those involving complex structural rearrangements. The description of such complex variants challenges the existing standard sequence variation nomenclature of the Human Genome Variation Society (HGVS, http://www.hgvs.org/mutnomen), because this mainly focuses on simple variants. Here, we suggest several extensions of the HGVS nomenclature guidelines to facilitate unambiguous description of complex sequence variants at the DNA level. These include: (1) nesting to support description of changes within inversions and duplications, and (2) composite changes to support concatenation of inserted sequences. The advantage of these additions is that inversions and duplications with small differences and more complex variants can be described without reverting to the less informative indel description. In addition, they should provide sufficient flexibility and consistency, thereby limiting alternative interpretations and ambiguous descriptions. The specifications should allow easy implementation in sequence variant nomenclature checkers (e.g., Mutalyzer, http://www.mutalyzer.nl/). We are extending the functionality of Mutalyzer to incorporate the latest version of the HGVS sequence variation nomenclature guidelines.

LOVD v.2.0: the next generation in gene variant databases.

Locus-Specific DataBases (LSDBs) store information on gene sequence variation associated with human phenotypes and are frequently used as a reference by researchers and clinicians. We developed the Leiden Open-source Variation Database (LOVD) as a platform-independent Web-based LSDB-in-a-Box package. LOVD was designed to be easy to set up and maintain and follows the Human Genome Variation Society (HGVS) recommendations. Here we describe LOVD v.2.0, which adds enhanced flexibility and functionality and has the capacity to store sequence variants in multiple genes per patient. To reduce redundancy, patient and sequence variant data are stored in separate tables. Tables are linked to generate connections between sequence variant data for each gene and every patient. The dynamic structure allows database managers to add custom columns. The database structure supports fast queries and allows storage of sequence variants from high-throughput sequence analysis, as demonstrated by the X-chromosomal Mental Retardation LOVD installation. LOVD contains measures to ensure database security from unauthorized access. Currently, the LOVD Website (http://www.LOVD.nl/) lists 71 public LOVD installations hosting 3,294 gene variant databases with 199,000 variants in 84,000 patients. To promote LSDB standardization and thereby database interoperability, we offer free server space and help to establish an LSDB on our Leiden server.

Deep sequencing to reveal new variants in pooled DNA samples.

We evaluated massive parallel sequencing and long-range PCR (LRP) for rare variant detection and allele frequency estimation in pooled DNA samples. Exons 2 to 16 of the MUTYH gene were analyzed in breast cancer patients with Illumina's (Solexa) technology. From a pool of 287 genomic DNA samples we generated a single LRP product, while the same LRP was performed on 88 individual samples and the resulting products then pooled. Concentrations of constituent samples were measured with fluorimetry for genomic DNA and high-resolution melting curve analysis (HR-MCA) for LRP products. Illumina sequencing results were compared to Sanger sequencing data of individual samples. Correlation between allele frequencies detected by both methods was poor in the first pool, presumably because the genomic samples amplified unequally in the LRP, due to DNA quality variability. In contrast, allele frequencies correlated well in the second pool, in which all expected alleles at a frequency of 1% and higher were reliably detected, plus the majority of singletons (0.6% allele frequency). We describe custom bioinformatics and statistics to optimize detection of rare variants and to estimate required sequencing depth. Our results provide directions for designing high-throughput analyses of candidate genes.

Caenorhabditis elegans as a model for lysosomal storage disorders.

The nematode Caenorhabditis elegans is the simplest animal model available to study human disease. In this review, the worm homologues for the 58 human genes involved in lysosomal storage disorders and for 105 human genes associated with lysosomal function have been compiled. Most human genes had at least one worm homologue. In addition, the phenotypes of 147 mutants, in which these genes have been disrupted or knocked down, have been summarized and discussed. The phenotypic spectrum of worm models of lysosomal storage disorders varies from lethality to none obvious, with a large variety of intermediate phenotypes. The genetic power of C. elegans provides a means to identify genes involved in specific processes with relative ease. The overview of potential lysosomal phenotypes presented here might be used as a starting point for the phenotypic characterization of newly developed knock-out models or for the design of genetic screens selecting for loss or gain of suitable knock-out model phenotypes. Screens for genes involved in lysosomal biogenesis and function have been performed successfully resulting in the cup and glo mutants, but screens involving subtle phenotypes are likely to be difficult.

Improving sequence variant descriptions in mutation databases and literature using the Mutalyzer sequence variation nomenclature checker.

Unambiguous and correct sequence variant descriptions are of utmost importance, not in the least since mistakes and uncertainties may lead to undesired errors in clinical diagnosis. We developed the Mutation Analyzer (Mutalyzer) sequence variation nomenclature checker (www.lovd.nl/mutalyzer; last accessed 13 September 2007) for automated analysis and correction of sequence variant descriptions using reference sequences from any organism. Mutalyzer handles most variation types: substitution, deletion, duplication, insertion, indel, and splice-site changes following current recommendations of the Human Genome Variation Society (HGVS). Input is a GenBank accession number or an uploaded reference sequence file in GenBank format with user-modified annotation, an HGNC gene symbol, and the variant (single or in a batch file). Mutalyzer generates variant descriptions at DNA level, the level of all annotated transcripts and the deduced outcome at protein level. To validate Mutalyzer's performance and to investigate the sequence variant description quality in locus-specific mutation databases (LSDBs), more than 11,000 variants in the PAH, BIC BRCA2, and HbVar databases were analyzed, showing that 87%, 25%, and 38%, respectively, were error-free and following the recommendations. Low recognition rates in BIC and HbVar (38% and 51%, respectively) were due to lack of a well-annotated genomic reference sequence (HbVar) or noncompliance to the guidelines (BRCA2). Provided with well-annotated genomic reference sequences, Mutalyzer is very effective for the curation of newly discovered sequence variation descriptions and existing LSDB data. Mutalyzer will be linked to the Leiden Open source Variation Database (LOVD) (www.LOVD.nl; last accessed 13 September 2007) and is the first module of a sequence variant effect prediction package.

Characterizing pathogenic processes in Batten disease: use of small eukaryotic model systems.

The neuronal ceroid lipofuscinoses (NCLs) are neurodegenerative disorders. Nevertheless, small model organisms, including those lacking a nervous system, have proven invaluable in the study of mechanisms that underlie the disease and in studying the functions of the conserved proteins associated to each disease. From the single-celled yeast, Saccharomyces cerevisiae and Schizosaccharomyces pombe, to the worm, Caenorhabditis elegans and the fruitfly, Drosophila melanogaster, biochemical and, in particular, genetic studies on these organisms have provided insight into the NCLs.

Recommended Guidelines for Validation, Quality Control, and Reporting of TP53 Variants in Clinical Practice.

Accurate assessment of TP53 gene status in sporadic tumors and in the germline of individuals at high risk of cancer due to Li-Fraumeni Syndrome (LFS) has important clinical implications for diagnosis, surveillance, and therapy. Genomic data from more than 20,000 cancer genomes provide a wealth of information on cancer gene alterations and have confirmed TP53 as the most commonly mutated gene in human cancer. Analysis of a database of 70,000 TP53 variants reveals that the two newly discovered exons of the gene, exons 9ß and 9γ, generated by alternative splicing, are the targets of inactivating mutation events in breast, liver, and head and neck tumors. Furthermore, germline rearrange-ments in intron 1 of TP53 are associated with LFS and are frequently observed in sporadic osteosarcoma. In this context of constantly growing genomic data, we discuss how screening strategies must be improved when assessing TP53 status in clinical samples. Finally, we discuss how TP53 alterations should be described by using accurate nomenclature to avoid confusion in scientific and clinical reports. Cancer Res; 77(6); 1250-60. ©2017 AACR.

Mutation analysis of SDHB and SDHC: novel germline mutations in sporadic head and neck paraganglioma and familial paraganglioma and/or pheochromocytoma.

BACKGROUND: Germline mutations of the SDHD, SDHB and SDHC genes, encoding three of the four subunits of succinate dehydrogenase, are a major cause of hereditary paraganglioma and pheochromocytoma, and demonstrate that these genes are classic tumor suppressors. Succinate dehydrogenase is a heterotetrameric protein complex and a component of both the Krebs cycle and the mitochondrial respiratory chain (succinate:ubiquinone oxidoreductase or complex II). METHODS: Using conformation sensitive gel electrophoresis (CSGE) and direct DNA sequencing to analyse genomic DNA from peripheral blood lymphocytes, here we describe the mutation analysis of the SDHB and SDHC genes in 37 patients with sporadic (i.e. no known family history) head and neck paraganglioma and five pheochromocytoma and/or paraganglioma families. RESULTS: Two sporadic patients were found to have a SDHB splice site mutation in intron 4, c.423+1G>A, which produces a mis-spliced transcript with a 54 nucleotide deletion, resulting in an 18 amino acid in-frame deletion. A third patient was found to carry the c.214C>T (p.Arg72Cys) missense mutation in exon 4 of SDHC, which is situated in a highly conserved protein motif that constitutes the quinone-binding site of the succinate: ubiquinone oxidoreductase (SQR) complex in E. coli. Together with our previous results, we found 27 germline mutations of SDH genes in 95 cases (28%) of sporadic head and neck paraganglioma. In addition all index patients of five families showing hereditary pheochromocytoma-paraganglioma were found to carry germline mutations of SDHB: four of which were novel, c.343C>T (p.Arg115X), c.141G>A (p.Trp47X), c.281G>A (p.Arg94Lys), and c.653G>C (p.Trp218Ser), and one reported previously, c.136C>T, p.Arg46X. CONCLUSION: In conclusion, these data indicate that germline mutations of SDHB and SDHC play a minor role in sporadic head and neck paraganglioma and further underline the importance of germline SDHB mutations in cases of familial pheochromocytoma-paraganglioma.

Somatic loss of maternal chromosome 11 causes parent-of-origin-dependent inheritance in SDHD-linked paraganglioma and phaeochromocytoma families.

Germline mutations in succinate dehydrogenase subunits B, C and D (SDHB, SDHC and SDHD), genes encoding subunits of mitochondrial complex II, cause hereditary paragangliomas and phaeochromocytomas. In SDHB (1p36)- and SDHC (1q21)-linked families, disease inheritance is autosomal dominant. In SDHD (11q23)-linked families, the disease phenotype is expressed only upon paternal transmission of the mutation, consistent with maternal imprinting. However, SDHD shows biallelic expression in brain, kidney and lymphoid tissues (Baysal et al., 2000). Moreover, consistent loss of the wild-type (wt) maternal allele in SDHD-linked tumours suggests expression of the maternal SDHD allele in normal paraganglia. Here we demonstrate exclusive loss of the entire maternal chromosome 11 in SDHD-linked paragangliomas and phaeochromocytomas, suggesting that combined loss of the wt SDHD allele and maternal 11p region is essential for tumorigenesis. We hypothesize that this is driven by selective loss of one or more imprinted genes in the 11p15 region. In paternally, but not in maternally derived SDHD mutation carriers, this can be achieved by a single event, that is, non-disjunctional loss of the maternal chromosome 11. Thus, the exclusive paternal transmission of the disease can be explained by a somatic genetic mechanism targeting both the SDHD gene on 11q23 and a paternally imprinted gene on 11p15.5, rather than imprinting of SDHD.

LOVD: easy creation of a locus-specific sequence variation database using an "LSDB-in-a-box" approach.

The completion of the human genome project has initiated, as well as provided the basis for, the collection and study of all sequence variation between individuals. Direct access to up-to-date information on sequence variation is currently provided most efficiently through web-based, gene-centered, locus-specific databases (LSDBs). We have developed the Leiden Open (source) Variation Database (LOVD) software approaching the "LSDB-in-a-Box" idea for the easy creation and maintenance of a fully web-based gene sequence variation database. LOVD is platform-independent and uses PHP and MySQL open source software only. The basic gene-centered and modular design of the database follows the recommendations of the Human Genome Variation Society (HGVS) and focuses on the collection and display of DNA sequence variations. With minimal effort, the LOVD platform is extendable with clinical data. The open set-up should both facilitate and promote functional extension with scripts written by the community. The LOVD software is freely available from the Leiden Muscular Dystrophy pages (www.DMD.nl/LOVD/). To promote the use of LOVD, we currently offer curators the possibility to set up an LSDB on our Leiden server.

The SDH mutation database: an online resource for succinate dehydrogenase sequence variants involved in pheochromocytoma, paraganglioma and mitochondrial complex II deficiency.

BACKGROUND: The SDHA, SDHB, SDHC and SDHD genes encode the subunits of succinate dehydrogenase (succinate: ubiquinone oxidoreductase), a component of both the Krebs cycle and the mitochondrial respiratory chain. SDHA, a flavoprotein and SDHB, an iron-sulfur protein together constitute the catalytic domain, while SDHC and SDHD encode membrane anchors that allow the complex to participate in the respiratory chain as complex II. Germline mutations of SDHD and SDHB are a major cause of the hereditary forms of the tumors paraganglioma and pheochromocytoma. The largest subunit, SDHA, is mutated in patients with Leigh syndrome and late-onset optic atrophy, but has not as yet been identified as a factor in hereditary cancer. DESCRIPTION: The SDH mutation database is based on the recently described Leiden Open (source) Variation Database (LOVD) system. The variants currently described in the database were extracted from the published literature and in some cases annotated to conform to current mutation nomenclature. Researchers can also directly submit new sequence variants online. Since the identification of SDHD, SDHC, and SDHB as classic tumor suppressor genes in 2000 and 2001, studies from research groups around the world have identified a total of 120 variants. Here we introduce all reported paraganglioma and pheochromocytoma related sequence variations in these genes, in addition to all reported mutations of SDHA. The database is now accessible online. CONCLUSION: The SDH mutation database offers a valuable tool and resource for clinicians involved in the treatment of patients with paraganglioma-pheochromocytoma, clinical geneticists needing an overview of current knowledge, and geneticists and other researchers needing a solid foundation for further exploration of both these tumor syndromes and SDHA-related phenotypes.

Hereditary paraganglioma due to the SDHD M1I mutation in a second Chinese family: a founder effect?

OBJECTIVES/HYPOTHESIS: Hereditary paraganglioma is a rare condition that is inherited in an autosomal-dominant fashion. Four distinct loci have been associated with hereditary paraganglioma, including the SDHD, SDHC, and SDHB genes and a locus at 11q13. The SDHD, SDHC, and SDHB genes code for subunits of succinate dehydrogenase, which forms part of the mitochondrial respiratory chain. SDHD mutations are widely distributed along the gene with no apparent hot spots, although a founder effect has been described in the Dutch population. METHODS: Following a prior report of the SDHD M1I mutation in an Australian Chinese family, a second Chinese family with the same mutation is reported. The proband developed bilateral head and neck paragangliomas at age 34 years and a functioning adrenal pheochromocytoma and two extra-adrenal abdominal paragangliomas 7 years later. His brother had unilateral head and neck paraganglioma at age 39 years. Given the multicentricity of the proband's tumor and the familial clustering of paragangliomas, a clinical diagnosis of hereditary paraganglioma was made, and the proband was tested for a mutation in the SDHD gene. RESULTS: The proband was found to be heterozygous for the SDHD MII mutation that removes the start codon, and his brother subsequently tested positive for the same mutation. The family is not related to the Australian Chinese family. CONCLUSION: The finding suggests the possibility of a founder effect in the Chinese population and warrants further investigation.