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.

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.

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.

Clinical Interpretation of Variants from Next-Generation Sequencing: The 2016 Scientific Meeting of the Human Genome Variation Society.

The 2016 scientific meeting of the Human Genome Variation Society (HGVS; http://www.hgvs.org) was held on the 20th of May in Barcelona, Spain, with the theme of "Clinical Interpretation of Variants from Next-Generation Sequencing."

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.

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.

The Matchmaker Exchange: a platform for rare disease gene discovery.

There are few better examples of the need for data sharing than in the rare disease community, where patients, physicians, and researchers must search for "the needle in a haystack" to uncover rare, novel causes of disease within the genome. Impeding the pace of discovery has been the existence of many small siloed datasets within individual research or clinical laboratory databases and/or disease-specific organizations, hoping for serendipitous occasions when two distant investigators happen to learn they have a rare phenotype in common and can "match" these cases to build evidence for causality. However, serendipity has never proven to be a reliable or scalable approach in science. As such, the Matchmaker Exchange (MME) was launched to provide a robust and systematic approach to rare disease gene discovery through the creation of a federated network connecting databases of genotypes and rare phenotypes using a common application programming interface (API). The core building blocks of the MME have been defined and assembled. Three MME services have now been connected through the API and are available for community use. Additional databases that support internal matching are anticipated to join the MME network as it continues to grow.

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.

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").

In search of triallelism in Bardet-Biedl syndrome.

Bardet-Biedl syndrome (BBS) is a model disease for ciliopathy in humans. The remarkable genetic heterogeneity that characterizes this disease is consistent with accumulating data on the interaction between the proteins encoded by the 14 BBS genes identified to date. Previous reports suggested that such interaction may also extend to instances of oligogenic inheritance in the form of triallelism which defies the long held view of BBS as an autosomal recessive disease. In order to investigate the magnitude of triallelism in BBS, we conducted a comprehensive analysis of all 14 BBS genes as well as the CCDC28B-modifier gene in a cohort of 29 BBS families, most of which are multiplex. Two in trans mutations in a BBS gene were identified in each of these families for a total of 20 mutations including 12 that are novel. In no instance did we observe two mutations in unaffected members of a given family, or observe the presence of a third allele that convincingly acted as a modifier of penetrance and supported the triallelic model of BBS. In addition to presenting a comprehensive genotype/phenotype overview of a large set of BBS mutations, including the occurrence of nonsyndromic retinitis pigmentosa in a family with a novel BBS9 mutation, our study argues in favor of straightforward autosomal recessive BBS in most cases.

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.

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.

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.

Sdhd and SDHD/H19 knockout mice do not develop paraganglioma or pheochromocytoma.

BACKGROUND: Mitochondrial succinate dehydrogenase (SDH) is a component of both the tricarboxylic acid cycle and the electron transport chain. Mutations of SDHD, the first protein of intermediary metabolism shown to be involved in tumorigenesis, lead to the human tumors paraganglioma (PGL) and pheochromocytoma (PC). SDHD is remarkable in showing an 'imprinted' tumor suppressor phenotype. Mutations of SDHD show a very high penetrance in man and we postulated that knockout of Sdhd would lead to the development of PGL/PC, probably in aged mice. METHODOLOGY/PRINCIPAL FINDINGS: We generated a conventional knockout of Sdhd in the mouse, removing the entire third exon. We also crossed this mouse with a knockout of H19, a postulated imprinted modifier gene of Sdhd tumorigenesis, to evaluate if loss of these genes together would lead to the initiation or enhancement of tumor development. Homozygous knockout of Sdhd results in embryonic lethality. No paraganglioma or other tumor development was seen in Sdhd KO mice followed for their entire lifespan, in sharp contrast to the highly penetrant phenotype in humans. Heterozygous Sdhd KO mice did not show hyperplasia of paraganglioma-related tissues such as the carotid body or of the adrenal medulla, or any genotype-related pathology, with similar body and organ weights to wildtype mice. A cohort of Sdhd/H19 KO mice developed several cases of profound cardiac hypertrophy, but showed no evidence of PGL/PC. CONCLUSIONS: Knockout of Sdhd in the mouse does not result in a disease phenotype. H19 may not be an initiator of PGL/PC tumorigenesis.

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.