The V-ATPase is the main regulator of intra-organellar acidification. Assembly of this complex has extensively been studied in yeast, while limited knowledge exists for man. We identified 11 male patients with hemizygous missense mutations in ATP6AP1, encoding accessory protein Ac45 of the V-ATPase. Homology detection at the level of sequence profiles indicated Ac45 as the long-sought human homologue of yeast V-ATPase assembly factor Voa1. Processed wild-type Ac45, but not its disease mutants, restored V-ATPase-dependent growth in Voa1 mutant yeast. Patients display an immunodeficiency phenotype associated with hypogammaglobulinemia, hepatopathy and a spectrum of neurocognitive abnormalities. Ac45 in human brain is present as the common, processed ∼40-kDa form, while liver shows a 62-kDa intact protein, and B-cells a 50-kDa isoform. Our work unmasks Ac45 as the functional ortholog of yeast V-ATPase assembly factor Voa1 and reveals a novel link of tissue-specific V-ATPase assembly with immunoglobulin production and cognitive function.

Retinitis pigmentosa (RP) is a group of inherited progressive retinal dystrophies (RD) characterized by photoreceptor degeneration. RP is highly heterogeneous both clinically and genetically, which complicates the identification of causative genes and mutations. Targeted next-generation sequencing (NGS) has been demonstrated to be an effective strategy for the detection of mutations in RP. In our study, an in-house gene panel comprising 75 known RP genes was used to analyze a cohort of 47 unrelated Spanish families pre-classified as autosomal recessive or isolated RP. Disease-causing mutations were found in 27 out of 47 cases achieving a mutation detection rate of 57.4%. In total, 33 pathogenic mutations were identified, 20 of which were novel mutations (60.6%). Furthermore, not only single nucleotide variations but also copy-number variations, including three large deletions in the USH2A and EYS genes, were identified. Finally seven out of 27 families, displaying mutations in the ABCA4, RP1, RP2 and USH2A genes, could be genetically or clinically reclassified. These results demonstrate the potential of our panel-based NGS strategy in RP diagnosis.

Mutations in more than a hundred genes have been reported to cause X-linked recessive intellectual disability (ID) mainly in males. In contrast, the number of identified X-linked genes in which de novo mutations specifically cause ID in females is limited. Here, we report 17 females with de novo loss-of-function mutations in USP9X, encoding a highly conserved deubiquitinating enzyme. The females in our study have a specific phenotype that includes ID/developmental delay (DD), characteristic facial features, short stature, and distinct congenital malformations comprising choanal atresia, anal abnormalities, post-axial polydactyly, heart defects, hypomastia, cleft palate/bifid uvula, progressive scoliosis, and structural brain abnormalities. Four females from our cohort were identified by targeted genetic testing because their phenotype was suggestive for USP9X mutations. In several females, pigment changes along Blaschko lines and body asymmetry were observed, which is probably related to differential (escape from) X-inactivation between tissues. Expression studies on both mRNA and protein level in affected-female-derived fibroblasts showed significant reduction of USP9X level, confirming the loss-of-function effect of the identified mutations. Given that some features of affected females are also reported in known ciliopathy syndromes, we examined the role of USP9X in the primary cilium and found that endogenous USP9X localizes along the length of the ciliary axoneme, indicating that its loss of function could indeed disrupt cilium-regulated processes. Absence of dysregulated ciliary parameters in affected female-derived fibroblasts, however, points toward spatiotemporal specificity of ciliary USP9X (dys-)function.

With the widespread adoption of next generation sequencing technologies by the genetics community and the rapid decrease in costs per base, exome sequencing has become a standard within the repertoire of genetic experiments for both research and diagnostics. Although bioinformatics now offers standard solutions for the analysis of exome sequencing data, many challenges still remain; especially the increasing scale at which exome data are now being generated has given rise to novel challenges in how to efficiently store, analyze and interpret exome data of this magnitude. In this review we discuss some of the recent developments in bioinformatics for exome sequencing and the directions that this is taking us to. With these developments, exome sequencing is paving the way for the next big challenge, the application of whole genome sequencing.

For next-generation sequencing technologies, sufficient base-pair coverage is the foremost requirement for the reliable detection of genomic variants. We investigated whether whole-genome sequencing (WGS) platforms offer improved coverage of coding regions compared with whole-exome sequencing (WES) platforms, and compared single-base coverage for a large set of exome and genome samples. We find that WES platforms have improved considerably in the last years, but at comparable sequencing depth, WGS outperforms WES in terms of covered coding regions. At higher sequencing depth (95x-160x), WES successfully captures 95% of the coding regions with a minimal coverage of 20x, compared with 98% for WGS at 87-fold coverage. Three different assessments of sequence coverage bias showed consistent biases for WES but not for WGS. We found no clear differences for the technologies concerning their ability to achieve complete coverage of 2,759 clinically relevant genes. We show that WES performs comparable to WGS in terms of covered bases if sequenced at two to three times higher coverage. This does, however, go at the cost of substantially more sequencing biases in WES approaches. Our findings will guide laboratories to make an informed decision on which sequencing platform and coverage to choose.

Over a relatively short period of time, the clinical geneticist's "toolbox" has been expanded by machine-learning algorithms for image analysis, which can be applied to the task of syndrome identification on the basis of facial photographs, but these technologies harbor potential beyond the recognition of established phenotypes. Here, we comprehensively characterized two individuals with a hitherto unknown genetic disorder caused by the same de novo mutation in LEMD2 (c.1436C>T;p.Ser479Phe), the gene which encodes the nuclear envelope protein LEM domain-containing protein 2 (LEMD2). Despite different ages and ethnic backgrounds, both individuals share a progeria-like facial phenotype and a distinct combination of physical and neurologic anomalies, such as growth retardation; hypoplastic jaws crowded with multiple supernumerary, yet unerupted, teeth; and cerebellar intention tremor. Immunofluorescence analyses of patient fibroblasts revealed mutation-induced disturbance of nuclear architecture, recapitulating previously published data in LEMD2-deficient cell lines, and additional experiments suggested mislocalization of mutant LEMD2 protein within the nuclear lamina. Computational analysis of facial features with two different deep neural networks showed phenotypic proximity to other nuclear envelopathies. One of the algorithms, when trained to recognize syndromic similarity (rather than specific syndromes) in an unsupervised approach, clustered both individuals closely together, providing hypothesis-free hints for a common genetic etiology. We show that a recurrent de novo mutation in LEMD2 causes a nuclear envelopathy whose prognosis in adolescence is relatively good in comparison to that of classical Hutchinson-Gilford progeria syndrome, and we suggest that the application of artificial intelligence to the analysis of patient images can facilitate the discovery of new genetic disorders.

We report on four families affected by a clinical presentation of complex hereditary spastic paraplegia (HSP) due to recessive mutations in DDHD2, encoding one of the three mammalian intracellular phospholipases A(1) (iPLA(1)). The core phenotype of this HSP syndrome consists of very early-onset (<2 years) spastic paraplegia, intellectual disability, and a specific pattern of brain abnormalities on cerebral imaging. An essential role for DDHD2 in the human CNS, and perhaps more specifically in synaptic functioning, is supported by a reduced number of active zones at synaptic terminals in Ddhd-knockdown Drosophila models. All identified mutations affect the protein's DDHD domain, which is vital for its phospholipase activity. In line with the function of DDHD2 in lipid metabolism and its role in the CNS, an abnormal lipid peak indicating accumulation of lipids was detected with cerebral magnetic resonance spectroscopy, which provides an applicable diagnostic biomarker that can distinguish the DDHD2 phenotype from other complex HSP phenotypes. We show that mutations in DDHD2 cause a specific complex HSP subtype (SPG54), thereby linking a member of the PLA(1) family to human neurologic disease.

We have characterized a novel autosomal recessive Crouzon-like craniosynostosis syndrome in a 12-affected member family from Antakya, Turkey, the presenting features of which include: multiple suture synostosis, midface hypoplasia, variable degree of exophthalmos, relative prognathism, a beaked nose, and conductive hearing loss. Homozygosity mapping followed by targeted next-generation sequencing identified a c.479+6T>G mutation in the interleukin 11 receptor alpha gene (IL11RA) on chromosome 9p21. This donor splice-site mutation leads to a high percentage of aberrant IL11RA mRNA transcripts in an affected individual and altered mRNA splicing determined by in vitro exon trapping. An extended IL11RA mutation screen was performed in a cohort of 79 patients with an initial clinical diagnosis of Crouzon syndrome, pansynostosis, or unclassified syndromic craniosynostosis. We identified mutations segregating with the disease in five families: a German patient of Turkish origin and a Turkish family with three affected sibs all of whom were homozygous for the previously identified IL11RA c.479+6T>G mutation; a family with pansynostosis with compound heterozygous missense mutations, p.Pro200Thr and p.Arg237Pro; and two further Turkish families with Crouzon-like syndrome carrying the homozygous nonsense mutations p.Tyr232* and p.Arg292*. Using transient coexpression in HEK293T and COS7 cells, we demonstrated dramatically reduced IL11-mediated STAT3 phosphorylation for all mutations. Immunofluorescence analysis of mouse Il11ra demonstrated specific protein expression in cranial mesenchyme which was localized around the coronal suture tips and in the lambdoidal suture. In situ hybridization analysis of adult zebrafish also detected zfil11ra expression in the coronal suture between the overlapping frontal and parietal plates. This study demonstrates that mutations in the IL11RA gene cause an autosomal recessive Crouzon-like craniosynostosis.

Whole exomes of patients with a genetic disorder are nowadays routinely sequenced but interpretation of the identified genetic variants remains a major challenge. The increased availability of population-based human genetic variation has given rise to measures of genetic tolerance that have been used, for example, to predict disease-causing genes in neurodevelopmental disorders. Here, we investigated whether combining variant information from homologous protein domains can improve variant interpretation. For this purpose, we developed a framework that maps population variation and known pathogenic mutations onto 2,750 "meta-domains." These meta-domains consist of 30,853 homologous Pfam protein domain instances that cover 36% of all human protein coding sequences. We find that genetic tolerance is consistent across protein domain homologues, and that patterns of genetic tolerance faithfully mimic patterns of evolutionary conservation. Furthermore, for a significant fraction (68%) of the meta-domains high-frequency population variation re-occurs at the same positions across domain homologues more often than expected. In addition, we observe that the presence of pathogenic missense variants at an aligned homologous domain position is often paired with the absence of population variation and vice versa. The use of these meta-domains can improve the interpretation of genetic variation.

Cantú syndrome is a rare disorder characterized by congenital hypertrichosis, neonatal macrosomia, a distinct osteochondrodysplasia, and cardiomegaly. Using an exome-sequencing approach applied to one proband-parent trio and three unrelated single cases, we identified heterozygous mutations in ABCC9 in all probands. With the inclusion of the remaining cohort of ten individuals with Cantú syndrome, a total of eleven mutations in ABCC9 were found. The de novo occurrence in all six simplex cases in our cohort substantiates the presence of a dominant disease mechanism. All mutations were missense, and several mutations affect Arg1154. This mutation hot spot lies within the second type 1 transmembrane region of this ATP-binding cassette transporter protein, which may suggest an activating mutation. ABCC9 encodes the sulfonylurea receptor (SUR) that forms ATP-sensitive potassium channels (K(ATP) channels) originally shown in cardiac, skeletal, and smooth muscle. Previously, loss-of-function mutations in this gene have been associated with idiopathic dilated cardiomyopathy type 10 (CMD10). These findings identify the genetic basis of Cantú syndrome and suggest that this is a new member of the potassium channelopathies.

Intellectual disability (ID) disorders are genetically and phenotypically highly heterogeneous and present a major challenge in clinical genetics and medicine. Although many genes involved in ID have been identified, the etiology is unknown in most affected individuals. Moreover, the function of most genes associated with ID remains poorly characterized. Evidence is accumulating that the control of gene transcription through epigenetic modification of chromatin structure in neurons has an important role in cognitive processes and in the etiology of ID. However, our understanding of the key molecular players and mechanisms in this process is highly fragmentary. Here, we identify a chromatin-modification module that underlies a recognizable form of ID, the Kleefstra syndrome phenotypic spectrum (KSS). In a cohort of KSS individuals without mutations in EHMT1 (the only gene known to be disrupted in KSS until now), we identified de novo mutations in four genes, MBD5, MLL3, SMARCB1, and NR1I3, all of which encode epigenetic regulators. Using Drosophila, we demonstrate that MBD5, MLL3, and NR1I3 cooperate with EHMT1, whereas SMARCB1 is known to directly interact with MLL3. We propose a highly conserved epigenetic network that underlies cognition in health and disease. This network should allow the design of strategies to treat the growing group of ID pathologies that are caused by epigenetic defects.

A significant component of the variation in cognitive disability that is observed in Duchenne muscular dystrophy (DMD) is known to be under genetic regulation. In this study we report correlations between standardised measures of intelligence and mutational class, mutation size, mutation location and the involvement of dystrophin isoforms.

The implementation of whole-exome sequencing in clinical diagnostics has generated a need for functional evaluation of genetic variants. In the field of inborn errors of metabolism (IEM), a diverse spectrum of targeted biochemical assays is employed to analyze a limited amount of metabolites. We now present a single-platform, high-resolution liquid chromatography quadrupole time of flight (LC-QTOF) method that can be applied for holistic metabolic profiling in plasma of individual IEM-suspected patients. This method, which we termed "next-generation metabolic screening" (NGMS), can detect >10,000 features in each sample. In the NGMS workflow, features identified in patient and control samples are aligned using the "various forms of chromatography mass spectrometry (XCMS)" software package. Subsequently, all features are annotated using the Human Metabolome Database, and statistical testing is performed to identify significantly perturbed metabolite concentrations in a patient sample compared with controls. We propose three main modalities to analyze complex, untargeted metabolomics data. First, a targeted evaluation can be done based on identified genetic variants of uncertain significance in metabolic pathways. Second, we developed a panel of IEM-related metabolites to filter untargeted metabolomics data. Based on this IEM-panel approach, we provided the correct diagnosis for 42 of 46 IEMs. As a last modality, metabolomics data can be analyzed in an untargeted setting, which we term "open the metabolome" analysis. This approach identifies potential novel biomarkers in known IEMs and leads to identification of biomarkers for as yet unknown IEMs. We are convinced that NGMS is the way forward in laboratory diagnostics of IEMs.

Intellectual disability (ID) affects approximately 1%-3% of humans with a gender bias toward males. Previous studies have identified mutations in more than 100 genes on the X chromosome in males with ID, but there is less evidence for de novo mutations on the X chromosome causing ID in females. In this study we present 35 unique deleterious de novo mutations in DDX3X identified by whole exome sequencing in 38 females with ID and various other features including hypotonia, movement disorders, behavior problems, corpus callosum hypoplasia, and epilepsy. Based on our findings, mutations in DDX3X are one of the more common causes of ID, accounting for 1%-3% of unexplained ID in females. Although no de novo DDX3X mutations were identified in males, we present three families with segregating missense mutations in DDX3X, suggestive of an X-linked recessive inheritance pattern. In these families, all males with the DDX3X variant had ID, whereas carrier females were unaffected. To explore the pathogenic mechanisms accounting for the differences in disease transmission and phenotype between affected females and affected males with DDX3X missense variants, we used canonical Wnt defects in zebrafish as a surrogate measure of DDX3X function in vivo. We demonstrate a consistent loss-of-function effect of all tested de novo mutations on the Wnt pathway, and we further show a differential effect by gender. The differential activity possibly reflects a dose-dependent effect of DDX3X expression in the context of functional mosaic females versus one-copy males, which reflects the complex biological nature of DDX3X mutations.

Noncanonical splice-site mutations are an important cause of inherited diseases. Based on in vitro and stem-cell-based studies, some splice-site variants show a stronger splice defect than expected based on their predicted effects, suggesting that other sequence motifs influence the outcome. We investigated whether splice defects due to human-inherited-disease-associated variants in noncanonical splice-site sequences in ABCA4, DMD, and TMC1 could be rescued by strengthening the splice site on the other side of the exon. Noncanonical 5'- and 3'-splice-site variants were selected. Rescue variants were introduced based on an increase in predicted splice-site strength, and the effects of these variants were analyzed using in vitro splice assays in HEK293T cells. Exon skipping due to five variants in noncanonical splice sites of exons in ABCA4, DMD, and TMC1 could be partially or completely rescued by increasing the predicted strengths of the other splice site of the same exon. We named this mechanism "splicing interdependency", and it is likely based on exon recognition by splicing machinery. Awareness of this interdependency is of importance in the classification of noncanonical splice-site variants associated with disease and may open new opportunities for treatments.

Mutations in retina-specific ATP-binding cassette transporter 4 (ABCA4) are responsible for over 95% of cases of Stargardt disease (STGD), as well as a minor proportion of retinitis pigmentosa (RP) and cone-rod dystrophy cases (CRD). Since the knowledge of the genetic causes of inherited retinal diseases (IRDs) in Poland is still scarce, the purpose of this study was to identify pathogenic ABCA4 variants in a subgroup of Polish IRD patients. We recruited 67 families with IRDs as a part of a larger study. The patients were screened with next generation sequencing using a molecular inversion probes (MIPs)-based technique targeting 108 genes involved in the pathogenesis of IRDs. All identified mutations were validated and their familial segregation was tested using Sanger sequencing. In the case of the most frequent complex allele, consisting of two variants in exon 12 and 21, familial segregation was tested using restriction fragment length polymorphism (RFLP). The most prevalent variant, a complex change c.[1622T>C;3113C>T], p.[Leu541Pro;Ala1038Val], was found in this cohort in 54% of all solved ABCA4-associated disorder cases, which is the highest frequency reported thus far. Additionally, we identified nine families displaying a pseudo-dominant mode of inheritance, indicating a high frequency of pathogenic variants within this population.

Leber congenital amaurosis (LCA) can be caused by mutations in more than 20 different genes. One of these, RPE65, encodes a protein essential for the visual cycle that is expressed in retinal pigment epithelium cells. In this work, we describe the generation and characterization of the human iPSC line SCTCi16-A. This hiPSC line was generated from peripheral blood mononuclear cells (PBMCs) from a patient affected with LCA caused by the homozygous c.11+5G>A variant in the RPE65 gene. Reprograming was conducted using episomal vectors containing OCT3/4, SOX2, KLF4, L-MYC, and LIN28.

Recurrent somatic variants in SPOP are cancer specific; endometrial and prostate cancers result from gain-of-function and dominant-negative effects toward BET proteins, respectively. By using clinical exome sequencing, we identified six de novo pathogenic missense variants in SPOP in seven individuals with developmental delay and/or intellectual disability, facial dysmorphisms, and congenital anomalies. Two individuals shared craniofacial dysmorphisms, including congenital microcephaly, that were strikingly different from those of the other five individuals, who had (relative) macrocephaly and hypertelorism. We measured the effect of SPOP variants on BET protein amounts in human Ishikawa endometrial cancer cells and patient-derived cell lines because we hypothesized that variants would lead to functional divergent effects on BET proteins. The de novo variants c.362G>A (p.Arg121Gln) and c. 430G>A (p.Asp144Asn), identified in the first two individuals, resulted in a gain of function, and conversely, the c.73A>G (p.Thr25Ala), c.248A>G (p.Tyr83Cys), c.395G>T (p.Gly132Val), and c.412C>T (p.Arg138Cys) variants resulted in a dominant-negative effect. Our findings suggest that these opposite functional effects caused by the variants in SPOP result in two distinct and clinically recognizable syndromic forms of intellectual disability with contrasting craniofacial dysmorphisms.

Joubert syndrome (JBTS) is a genetically heterogeneous group of neurodevelopmental syndromes caused by primary cilia dysfunction. Usually the neurological presentation starts with abnormal neonatal breathing followed by muscular hypotonia, psychomotor delay, and cerebellar ataxia. Cerebral MRI shows mid- and hindbrain anomalies including the molar tooth sign. We report a male patient with atypical presentation of Joubert syndrome type 23, thus expanding the phenotype.

Members of the chromodomain-helicase-DNA binding (CHD) protein family are chromatin remodelers implicated in human pathologies, with CHD6 being one of its least studied members. We discovered a de novo CHD6 missense mutation in a patient clinically presenting the rare Hallermann-Streiff syndrome (HSS). We used genome editing to generate isogenic iPSC lines and model HSS in relevant cell types. By combining genomics with functional in vivo and in vitro assays, we show that CHD6 binds a cohort of autophagy and stress response genes across cell types. The HSS mutation affects CHD6 protein folding and impairs its ability to recruit co-remodelers in response to DNA damage or autophagy stimulation. This leads to accumulation of DNA damage burden and senescence-like phenotypes. We therefore uncovered a molecular mechanism explaining HSS onset via chromatin control of autophagic flux and genotoxic stress surveillance.