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.

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.

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.

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.

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.

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.

Polymorphism of μ-Opioid Receptor Gene (OPRM1:c.118A>G ) Does Not Protect Against Opioid-induced Respiratory Depression despite Reduced Analgesic Response

Background

The effect of a single nucleotide polymorphism of the mu-opioid receptor at nucleotide position 118 (OPRM1:c.118A>G) was investigated on morphine-6-glucuronide (M6G)-induced analgesia and respiratory depression in a group of healthy volunteers.

Methods

Sixteen subjects of either sex received 0.4 mg/kg (n = 8) or 0.6 mg/kg M6G (n = 8). At regular time intervals, the isocapnic acute hypoxic ventilatory response, pain tolerance (derived from a transcutaneous electrical acute pain model), and arterial blood samples were obtained. Data acquisition continued for 14 h after drug infusion. Population pharmacokinetic-pharmacodynamic sigmoid Emax models were applied to the respiratory and pain data. All collected data were analyzed using the statistical program NONMEM (San Francisco, CA).

Results

Four of the subjects were OPRM1:c.118GA heterozygotes, and the remainder of the subjects were OPRM1:c.118AA homozygotes. M6G analgesia: In contrast to analgesic responses in OPRM1:c.118AA homozygotes, responses were small and inconsistent in OPRM1:c.118GA heterozygotes and best described by the function Effect(t) = baseline (P < 0.01 vs. OPRM1:c.118AA homozygotes). Emax and C50 values in heterozygotes equaled 0.55 +/- 0.18 (or a 55% increase in current above baseline) and 161 +/- 42 ng/ml, respectively. M6G-induced respiratory depression: For the acute hypoxic response, neither Emax nor C50 (value = 282 +/- 72 ng/ml) differed between genotypes.

Conclusions

The data indicate that the OPRM1:c.118A>G polymorphism affects opioid analgesic and respiratory effects differentially. Despite reduced analgesic responses to M6G the OPRM1:c.118A>G single-nucleotide polymorphism does not protect against the toxic effects of the tested opioid. However, some caution in the interpretation of the data is needed because of the small sample size. Further studies are needed to explore the link between this polymorphism and respiratory/analgesic responses beyond the small human sample. In OPRM1:c.118AA homozygotes, the potency parameters differed by a factor of 2 for analgesic versus respiratory effect. In this respect, M6G differs favorably from morphine.

[From gene to disease; from CLN1, CLN2 and CLN3 to neuronal ceroid lipofuscinosis]

The neuronal ceroid lipofuscinoses (NCL) are worldwide the most common lysosomal storage disorders of childhood. Clinical features often include progressive visual impairment, seizures, psychomotor deterioration, dementia, and premature death. Most NCL cases are caused by mutations in the CLN1, CLN2 and CLN3 genes, which play an essential role in lysosomal protein degradation. Laboratory diagnostics for a patient suspected of NCL should start with enzyme analysis in the case of INCL and LINCL and investigation of lymphocyte vacuolisation for JNCL. Diagnosis at the protein level is not available for JNCL, but CLN3 mutation analysis is possible. The carrier status of healthy relatives in families with known mutations in either CLN1, CLN2, CLN3 or CLN6 can be determined with certainty by mutation analysis.

Deletion of the Caenorhabditis elegans homologues of the CLN3 gene, involved in human juvenile neuronal ceroid lipofuscinosis, causes a mild progeric phenotype

The CLN3 gene is involved in juvenile neuronal ceroid lipofuscinosis (JNCL), or Batten-Spielmeyer-Vogt disease, a severe hereditary neurodegenerative lysosomal storage disorder characterized by progressive disease pathology, with loss of vision as the first symptom. Another characteristic of JNCL is the lysosomal accumulation of autofluorescent lipopigments, forming fingerprint storage patterns visible by electron microscopy. The function of the CLN3 protein is still unknown, although the evolutionarily conserved CLN3 protein is being functionally analysed using different experimental models. We have explored the potential of the nematode Caenorhabditis elegans as a model for Batten disease in order to bridge the gap between the unicellular yeast and very complex mouse JNCL models. C. elegans has three genes homologous to CLN3, for each of which deletion mutants were isolated. Cln-3.1 deletion mutants have a decreased lifespan, and cln-3.2 deletion mutants a decreased brood size. However, the neuronal or movement defects and aberrant lipopigment distribution or accumulation observed in JNCL were not found in the worms. To detect possible redundancy, single deletion mutants were crossed to obtain double and triple mutants, which were viable but showed no JNCL-specific defects. The cln-3 triple mutants show a more prominent decrease in lifespan and brood size, the latter most conspicuously at the end of the egg-laying period, suggesting premature ageing. To focus our functional analysis we examined the C. elegans cln-3 expression patterns, using promoter-GFP (green fluorescent protein) gene fusions. Fluorescence patterns suggest cln-3.1 expression in the intestine, cln-3.2 expression in the hypoderm, and cln-3.3 expression in intestinal muscle, male-specific posterior muscle and hypoderm. Further life stage- and tissue-specific analysis of the processes causing the phenotype of the cln-3 triple mutants may provide more information about the function of the cln-3 protein and contribute to a better understanding of the basic processes affected in Batten disease patients.

Increased prevalence of catecholamine excess and phaeochromocytomas in a well-defined Dutch population with SDHD-linked head and neck paragangliomas

OBJECTIVE

The aim of this study was to identify the prevalence of catecholamine excess and phaeochromocytomas in a well-defined population of people with hereditary head and neck paragangliomas.

METHODS

We studied in a prospective follow-up protocol all consecutive patients referred to the Department of Endocrinology, Leiden University Medical Center, Leiden, The Netherlands with documented head and neck paragangliomas and either a positive family history for paragangliomas or a proven SDHD gene mutation. Initial analysis included medical history, physical examination and the measurement of excretion of catecholamines in two 24-h urine collections. In the case of documented catecholamine excess iodinated meta-iodobenzylguanidine (123I-MIBG) scintigraphy and magnetic resonance imaging were done.

RESULTS

Between 1988 and 2003, 40 consecutive patients (20 male and 20 female) with documented head and neck paragangliomas were screened. Biochemical screening revealed urinary catecholamine excess in 15 patients (37.5%). In nine of these 15 patients a lesion was found by 123I-MIBG scintigraphy. Exact localization by magnetic resonance imaging revealed phaeochromocytomas in seven of the 15 patients. One of the nine patients had an extra-adrenal paraganglioma. Histopathological examination in a subset of tumors displayed loss of heterozygosity of the wild-type SDHD allele in all cases.

CONCLUSIONS

The prevalence of catecholamine excess (37.5%) and phaeochromocytomas (20.0%) is high in patients with familial head and neck paragangliomas. Therefore, patients with hereditary head and neck paragangliomas require lifelong follow up by biochemical testing for catecholamine excess.

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.

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.

Pharmacokinetic-pharmacodynamic modeling of morphine-6-glucuronide-induced analgesia in healthy volunteers: absence of sex differences

BACKGROUND

Morphine-6-glucuronide (M6G) is a metabolite of morphine and a micro-opioid agonist. To quantify the potency and speed of onset-offset of M6G and explore putative sex dependency, the authors studied the pharmacokinetics and pharmacodynamics of M6G in volunteers using a placebo-controlled, randomized, double-blind study design.

METHODS

Ten men and 10 women received 0.3 mg/kg intravenous M6G and placebo (two thirds of the dose as bolus, one third as a continuous infusion over 1 h) on separate occasions. For 7 h, pain tolerance was measured using gradually increasing transcutaneous electrical stimulation, and blood samples were obtained. A population pharmacokinetic (inhibitory sigmoid Emax)-pharmacodynamic analysis was used to analyze M6G-induced changes in tolerated stimulus intensity. The improvement in model fits by inclusion of covariate sex was tested for significance. P values less than 0.01 were considered significant. Taking into account previous morphine data, a predictive pharmacokinetic-pharmacodynamic model was constructed to determine the contribution of M6G to morphine analgesia.

RESULTS

M6G concentrations did not differ between men and women. M6G caused analgesia significantly greater than that observed with placebo (P < 0.01). The M6G analgesia data were well described by the pharmacokinetic-pharmacodynamic model. The M6G effect site concentration causing a 25% increase in current (C25) was 275 +/- 135 nm (population estimate +/- SE), the blood effect site equilibration half-life was 6.2 +/- 3.3 h, and the steepness parameter was 0.71 +/- 0.18. Intersubject variability was 167% for C25 and 218% for the effect half-life. None of the model parameters showed sex dependency.

CONCLUSIONS

A cumulative dose of 0.3 mg/kg M6G, given over 1 h, produces long-term analgesia greater than that observed with placebo, with equal dynamics (potency and speed of onset-offset) in men and women. Possible causes for the great intersubject response variability, such as genetic polymorphism of the micro-opioid receptor and placebo-related phenomena, are discussed. The predictive pharmacokinetic-pharmacodynamic model was applied successfully and was used to estimate M6G analgesia after morphine in patients with normal and impaired renal function.

SDHD mutations in head and neck paragangliomas result in destabilization of complex II in the mitochondrial respiratory chain with loss of enzymatic activity and abnormal mitochondrial morphology

Hereditary head and neck paragangliomas are tumours associated with the autonomic nervous system. Recently, mutations in genes coding for subunits of mitochondrial complex II, succinate-ubiquinone-oxidoreductase (SDHB, SDHC, and SDHD), have been identified in the majority of hereditary tumours and a number of isolated cases. In addition, a fourth locus, PGL2, has been mapped to chromosome 11q13 in an isolated family. In order to characterize phenotypic effects of these mutations, the present study investigated the immunohistochemical expression of the catalytic subunits of complex II (flavoprotein and iron protein), SDH enzyme activity, and mitochondrial morphology in a series of 22 head and neck paragangliomas. These included 11 SDHD-, one SDHB-, two PGL2-linked tumours, and eight sporadic tumours. In the majority of the tumours (approximately 90%), the enzyme-histochemical SDH reaction was negative and immunohistochemistry of catalytic subunits of complex II showed reduced expression of iron protein and enhanced expression of flavoprotein. Ultrastructural examination revealed elevated numbers of tightly packed mitochondria with abnormal morphology in SDHD-linked and sporadic tumours. Immuno-electron microscopy showed localization of the flavoprotein on the remnants of the mitochondrial inner membranes, whereas virtually no signal for the iron protein was detected. These results indicate that the function of mitochondrial complex II is compromised in the majority of head and neck paragangliomas.

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.

[From gene to disease; from SDHD, a defect in the respiratory chain, to paragangliomas and pheochromocytomas]

Hereditary paragangliomas are rare benign tumours arising from neuroectodermal tissue in the head and neck region. In families with paraganglioma, occasionally adrenal and extra-adrenal pheochromocytomas are found. Paragangliomas, adrenal and extra-adrenal pheochromocytomas may be caused by mutations in the SDHB, SDHC and SDHD genes encoding different subunits of mitochondrial respiratory chain complex II. Most paraganglioma cases in the Netherlands are caused by SDHD mutations. Presymptomatic DNA diagnosis is available for families with paragangliomas caused by SDHD mutations.