Ulrike Krüger

Effective diagnosis of genetic disease by computational phenotype analysis of the disease-associated genome

Patients with genetic disease of unknown causes can be rapidly diagnosed by bioinformatic analysis of disease-associated DNA sequences and phenotype. Efficient Diagnosis of Genetic Disease We know which genes are mutated in almost 3000 inherited human diseases and have good descriptions of how these mutations affect the human phenotype. Now, Zemojtel et al. have coupled this knowledge with rapid sequencing of these genes in a group of 40 patients with undiagnosed genetic diseases. Bioinformatic matching of the patients’ clinical characteristics and their disease gene sequences to databases of current genetic and phenotype knowledge enabled the authors to successfully diagnose almost 30% of the patients. The process required only about 2 hours of a geneticists’ time. Zemojtel et al. have made their tools available to the community, enabling a fast straightforward process by which clinicians and patients can easily identify the genetic basis of inherited disease in certain people. Less than half of patients with suspected genetic disease receive a molecular diagnosis. We have therefore integrated next-generation sequencing (NGS), bioinformatics, and clinical data into an effective diagnostic workflow. We used variants in the 2741 established Mendelian disease genes [the disease-associated genome (DAG)] to develop a targeted enrichment DAG panel (7.1 Mb), which achieves a coverage of 20-fold or better for 98% of bases. Furthermore, we established a computational method [Phenotypic Interpretation of eXomes (PhenIX)] that evaluated and ranked variants based on pathogenicity and semantic similarity of patients’ phenotype described by Human Phenotype Ontology (HPO) terms to those of 3991 Mendelian diseases. In computer simulations, ranking genes based on the variant score put the true gene in first place less than 5% of the time; PhenIX placed the correct gene in first place more than 86% of the time. In a retrospective test of PhenIX on 52 patients with previously identified mutations and known diagnoses, the correct gene achieved a mean rank of 2.1. In a prospective study on 40 individuals without a diagnosis, PhenIX analysis enabled a diagnosis in 11 cases (28%, at a mean rank of 2.4). Thus, the NGS of the DAG followed by phenotype-driven bioinformatic analysis allows quick and effective differential diagnostics in medical genetics.

Mutations in PIGO, a member of the GPI-anchor-synthesis pathway, cause hyperphosphatasia with mental retardation

Hyperphosphatasia with mental retardation syndrome (HPMRS), an autosomal-recessive form of intellectual disability characterized by facial dysmorphism, seizures, brachytelephalangy, and persistent elevated serum alkaline phosphatase (hyperphosphatasia), was recently shown to be caused by mutations in PIGV, a member of the glycosylphosphatidylinositol (GPI)-anchor-synthesis pathway. However, not all individuals with HPMRS harbor mutations in this gene. By exome sequencing, we detected compound-heterozygous mutations in PIGO, a gene coding for a membrane protein of the same molecular pathway, in two siblings with HPMRS, and we then found by Sanger sequencing further mutations in another affected individual; these mutations cosegregated in the investigated families. The mutant transcripts are aberrantly spliced, decrease the membrane stability of the protein, or impair enzyme function such that GPI-anchor synthesis is affected and the level of GPI-anchored substrates localized at the cell surface is reduced. Our data identify PIGO as the second gene associated with HPMRS and suggest that a deficiency in GPI-anchor synthesis is the underlying molecular pathomechanism of HPMRS.

PGAP2 Mutations, Affecting the GPI-Anchor-Synthesis Pathway, Cause Hyperphosphatasia with Mental Retardation Syndrome

Recently, mutations in genes involved in the biosynthesis of the glycosylphosphatidylinositol (GPI) anchor have been identified in a new subclass of congenital disorders of glycosylation (CDGs) with a distinct spectrum of clinical features. To date, mutations have been identified in six genes (PIGA, PIGL, PIGM, PIGN, PIGO, and PIGV) encoding proteins in the GPI-anchor-synthesis pathway in individuals with severe neurological features, including seizures, muscular hypotonia, and intellectual disability. We developed a diagnostic gene panel for targeting all known genes encoding proteins in the GPI-anchor-synthesis pathway to screen individuals matching these features, and we detected three missense mutations in PGAP2, c.46C>T, c.380T>C, and c.479C>T, in two unrelated individuals with hyperphosphatasia with mental retardation syndrome (HPMRS). The mutations cosegregated in the investigated families. PGAP2 is involved in fatty-acid GPI-anchor remodeling, which occurs in the Golgi apparatus and is required for stable association between GPI-anchored proteins and the cell-surface membrane rafts. Transfection of the altered protein constructs, p.Arg16Trp (NP_001243169.1), p.Leu127Ser, and p.Thr160Ile, into PGAP2-null cells showed only partial restoration of GPI-anchored marker proteins, CD55 and CD59, on the cell surface. In this work, we show that an impairment of GPI-anchor remodeling also causes HPMRS and conclude that targeted sequencing of the genes encoding proteins in the GPI-anchor-synthesis pathway is an effective diagnostic approach for this subclass of CDGs.

Mutations in PGAP3 Impair GPI-Anchor Maturation, Causing a Subtype of Hyperphosphatasia with Mental Retardation

Glycosylphophatidylinositol (GPI)-anchored proteins play important roles in many biological processes, and mutations affecting proteins involved in the synthesis of the GPI anchor are reported to cause a wide spectrum of intellectual disabilities (IDs) with characteristic additional phenotypic features. Here, we describe a total of five individuals (from three unrelated families) in whom we identified mutations in PGAP3, encoding a protein that is involved in GPI-anchor maturation. Three siblings in a consanguineous Pakistani family presented with profound developmental delay, severe ID, no speech, psychomotor delay, and postnatal microcephaly. A combination of autozygosity mapping and exome sequencing identified a 13.8 Mb region harboring a homozygous c.275G>A (p.Gly92Asp) variant in PGAP3 region 17q11.2-q21.32. Subsequent testing showed elevated serum alkaline phosphatase (ALP), a GPI-anchored enzyme, in all three affected children. In two unrelated individuals in a cohort with developmental delay, ID, and elevated ALP, we identified compound-heterozygous variants c.439dupC (p.Leu147Profs(∗)16) and c.914A>G (p.Asp305Gly) and homozygous variant c.314C>G (p.Pro105Arg). The 1 bp duplication causes a frameshift and nonsense-mediated decay. Further evidence supporting pathogenicity of the missense mutations c.275G>A, c.314C>G, and c.914A>G was provided by the absence of the variants from ethnically matched controls, phylogenetic conservation, and functional studies on Chinese hamster ovary cell lines. Taken together with recent data on PGAP2, these results confirm the importance of the later GPI-anchor remodelling steps for normal neuronal development. Impairment of PGAP3 causes a subtype of hyperphosphatasia with ID, a congenital disorder of glycosylation that is also referred to as Mabry syndrome.

A case of paroxysmal nocturnal hemoglobinuria caused by a germline mutation and a somatic mutation in PIGT

To ascertain the genetic basis of a paroxysmal nocturnal hemoglobinuria (PNH) case without somatic mutations in PIGA, we performed deep next-generation sequencing on all exons of known genes of the glycosylphosphatidylinositol (GPI) anchor synthesis pathway. We identified a heterozygous germline splice site mutation in PIGT and a somatic 8-MB deletion in granulocytes affecting the other copy of PIGT. PIGA is essential for GPI anchor synthesis, whereas PIGT is essential for attachment of the preassembled GPI anchor to proteins. Although a single mutation event in the X-chromosomal gene PIGA is known to cause GPI-anchored protein deficiency, 2 such hits are required in the autosomal gene PIGT. Our data indicate that PNH can occur even in the presence of fully assembled GPI if its transfer to proteins is defective in hematopoietic stem cells.

Recurrent De Novo Mutations Affecting Residue Arg138 of Pyrroline-5-Carboxylate Synthase Cause a Progeroid Form of Autosomal-Dominant Cutis Laxa

Progeroid disorders overlapping with De Barsy syndrome (DBS) are collectively denoted as autosomal-recessive cutis laxa type 3 (ARCL3). They are caused by biallelic mutations in PYCR1 or ALDH18A1, encoding pyrroline-5-carboxylate reductase 1 and pyrroline-5-carboxylate synthase (P5CS), respectively, which both operate in the mitochondrial proline cycle. We report here on eight unrelated individuals born to non-consanguineous families clinically diagnosed with DBS or wrinkly skin syndrome. We found three heterozygous mutations in ALDH18A1 leading to amino acid substitutions of the same highly conserved residue, Arg138 in P5CS. A de novo origin was confirmed in all six probands for whom parental DNA was available. Using fibroblasts from affected individuals and heterologous overexpression, we found that the P5CS-p.Arg138Trp protein was stable and able to interact with wild-type P5CS but showed an altered sub-mitochondrial distribution. A reduced size upon native gel electrophoresis indicated an alteration of the structure or composition of P5CS mutant complex. Furthermore, we found that the mutant cells had a reduced P5CS enzymatic activity leading to a delayed proline accumulation. In summary, recurrent de novo mutations, affecting the highly conserved residue Arg138 of P5CS, cause an autosomal-dominant form of cutis laxa with progeroid features. Our data provide insights into the etiology of cutis laxa diseases and will have immediate impact on diagnostics and genetic counseling.

Delineation of PIGV mutation spectrum and associated phenotypes in hyperphosphatasia with mental retardation syndrome

Three different genes of the glycosylphosphatidylinositol anchor synthesis pathway, PIGV, PIGO, and PGAP2, have recently been implicated in hyperphosphatasia-mental retardation syndrome (HPMRS), also known as Mabry syndrome, a rare autosomal recessive form of intellectual disability. The aim of this study was to delineate the PIGV mutation spectrum as well as the associated phenotypic spectrum in a cohort of 16 individuals diagnosed with HPMRS on the basis of intellectual disability and elevated serum alkaline phosphate as minimal diagnostic criteria. All PIGV exons and intronic boundaries were sequenced in 16 individuals. Biallelic PIGV mutations were identified in 8 of 16 unrelated families with HPMRS. The most frequent mutation detected in about 80% of affected families including the cases reported here is the c.1022C>A PIGV mutation, which was found in both the homozygous as well as the heterozygous state. Four further mutations found in this study (c. 176T>G, c.53G>A, c.905T>C, and c.1405C>T) are novel. Our findings in the largest reported cohort to date significantly extend the range of reported manifestations associated with PIGV mutations and demonstrate that the severe end of the clinical spectrum presents as a multiple congenital malformation syndrome with a high frequency of Hirschsprung disease, vesicoureteral, and renal anomalies as well as anorectal malformations. PIGV mutations are the major cause of HPMRS, which displays a broad clinical variability regarding associated malformations and growth patterns. Severe developmental delays, particular facial anomalies, brachytelephalangy, and hyperphosphatasia are consistently found in PIGV-positive individuals.

The allele distribution in next-generation sequencing data sets is accurately described as the result of a stochastic branching process

With the availability of next-generation sequencing (NGS) technology, it is expected that sequence variants may be called on a genomic scale. Here, we demonstrate that a deeper understanding of the distribution of the variant call frequencies at heterozygous loci in NGS data sets is a prerequisite for sensitive variant detection. We model the crucial steps in an NGS protocol as a stochastic branching process and derive a mathematical framework for the expected distribution of alleles at heterozygous loci before measurement that is sequencing. We confirm our theoretical results by analyzing technical replicates of human exome data and demonstrate that the variance of allele frequencies at heterozygous loci is higher than expected by a simple binomial distribution. Due to this high variance, mutation callers relying on binomial distributed priors are less sensitive for heterozygous variants that deviate strongly from the expected mean frequency. Our results also indicate that error rates can be reduced to a greater degree by technical replicates than by increasing sequencing depth.

First description of a patient with Vici syndrome due to a mutation affecting the penultimate exon of EPG5 and review of the literature

Vici syndrome is a rare autosomal recessively inherited multisystem disorder characterized by agenesis of the corpus callosum, cataracts, cardiomyopathy, combined immunodeficiency, psychomotor delay, and hypopigmentation. Cullup et al. recently identified mutations in the gene EPG5 as the cause of Vici syndrome. EPG5 is involved in autophagy, an evolutionarily conserved lysosomal degradation process that is essential for cell homeostasis. Following the first description in 1988 by Vici et al., 24 other cases of Vici syndrome have been published with variable expression of the defining features. Here, we report on a further case of Vici syndrome with a homozygous truncating mutation of EPG5, identified by whole‐exome sequencing. The mutation in our patient is the first reported affecting the penultimate exon of EPG5 and presenting with typical clinical manifestations of Vici syndrome. Additionally, we present a detailed clinical analysis of Vici syndrome comprising all cases previously described in the literature.

Screening for single nucleotide variants, small indels and exon deletions with a next-generation sequencing based gene panel approach for Usher syndrome

Usher syndrome is an autosomal recessive disorder characterized both by deafness and blindness. For the three clinical subtypes of Usher syndrome causal mutations in altogether 12 genes and a modifier gene have been identified. Due to the genetic heterogeneity of Usher syndrome, the molecular analysis is predestined for a comprehensive and parallelized analysis of all known genes by next‐generation sequencing (NGS) approaches. We describe here the targeted enrichment and deep sequencing for exons of Usher genes and compare the costs and workload of this approach compared to Sanger sequencing. We also present a bioinformatics analysis pipeline that allows us to detect single‐nucleotide variants, short insertions and deletions, as well as copy number variations of one or more exons on the same sequence data. Additionally, we present a flexible in silico gene panel for the analysis of sequence variants, in which newly identified genes can easily be included. We applied this approach to a cohort of 44 Usher patients and detected biallelic pathogenic mutations in 35 individuals and monoallelic mutations in eight individuals of our cohort. Thirty‐nine of the sequence variants, including two heterozygous deletions comprising several exons of USH2A, have not been reported so far. Our NGS‐based approach allowed us to assess single‐nucleotide variants, small indels, and whole exon deletions in a single test. The described diagnostic approach is fast and cost‐effective with a high molecular diagnostic yield.