The impact of inversions across 33,924 families with rare disease from a national genome sequencing project.

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

Recurrent de novo missense variants across multiple histone H4 genes underlie a neurodevelopmental syndrome.

Chromatin is essentially an array of nucleosomes, each of which consists of the DNA double-stranded fiber wrapped around a histone octamer. This organization supports cellular processes such as DNA replication, DNA transcription, and DNA repair in all eukaryotes. Human histone H4 is encoded by fourteen canonical histone H4 genes, all differing at the nucleotide level but encoding an invariant protein. Here, we present a cohort of 29 subjects with de novo missense variants in six H4 genes (H4C3, H4C4, H4C5, H4C6, H4C9, and H4C11) identified by whole-exome sequencing and matchmaking. All individuals present with neurodevelopmental features of intellectual disability and motor and/or gross developmental delay, while non-neurological features are more variable. Ten amino acids are affected, six recurrently, and are all located within the H4 core or C-terminal tail. These variants cluster to specific regions of the core H4 globular domain, where protein-protein interactions occur with either other histone subunits or histone chaperones. Functional consequences of the identified variants were evaluated in zebrafish embryos, which displayed abnormal general development, defective head organs, and reduced body axis length, providing compelling evidence for the causality of the reported disorder(s). While multiple developmental syndromes have been linked to chromatin-associated factors, missense-bearing histone variants (e.g., H3 oncohistones) are only recently emerging as a major cause of pathogenicity. Our findings establish a broader involvement of H4 variants in developmental syndromes.

Craniosynostosis, inner ear, and renal anomalies in a child with complete loss of SPRY1 (sprouty homolog 1) function.

INTRODUCTION: SPRY1 encodes protein sprouty homolog 1 (Spry-1), a negative regulator of receptor tyrosine kinase signalling. Null mutant mice display kidney/urinary tract abnormalities and altered size of the skull; complete loss-of-function of Spry-1 in humans has not been reported. METHODS: Analysis of whole-genome sequencing data from individuals with craniosynostosis enrolled in the 100,000 Genomes Project identified a likely pathogenic variant within SPRY1. Reverse-transcriptase PCR and western blot analysis were used to investigate the effect of the variant on SPRY1 mRNA and protein, in lymphoblastoid cell lines from the patient and both parents. RESULTS: A nonsense variant in SPRY1, encoding p.(Leu27*), was confirmed to be heterozygous in the unaffected parents and homozygous in the child. The child's phenotype, which included sagittal craniosynostosis, subcutaneous cystic lesions overlying the lambdoid sutures, hearing loss associated with bilateral cochlear and vestibular dysplasia and a unilateral renal cyst, overlapped the features reported in Spry1-/- null mice. Functional studies supported escape from nonsense-mediated decay, but western blot analysis demonstrated complete absence of full-length protein in the affected child and a marked reduction in both parents. CONCLUSION: This is the first report of complete loss of Spry-1 function in humans, associated with abnormalities of the cranial sutures, inner ear, and kidneys.

De Novo SOX6 Variants Cause a Neurodevelopmental Syndrome Associated with ADHD, Craniosynostosis, and Osteochondromas.

SOX6 belongs to a family of 20 SRY-related HMG-box-containing (SOX) genes that encode transcription factors controlling cell fate and differentiation in many developmental and adult processes. For SOX6, these processes include, but are not limited to, neurogenesis and skeletogenesis. Variants in half of the SOX genes have been shown to cause severe developmental and adult syndromes, referred to as SOXopathies. We here provide evidence that SOX6 variants also cause a SOXopathy. Using clinical and genetic data, we identify 19 individuals harboring various types of SOX6 alterations and exhibiting developmental delay and/or intellectual disability; the individuals are from 17 unrelated families. Additional, inconstant features include attention-deficit/hyperactivity disorder (ADHD), autism, mild facial dysmorphism, craniosynostosis, and multiple osteochondromas. All variants are heterozygous. Fourteen are de novo, one is inherited from a mosaic father, and four offspring from two families have a paternally inherited variant. Intragenic microdeletions, balanced structural rearrangements, frameshifts, and nonsense variants are predicted to inactivate the SOX6 variant allele. Four missense variants occur in residues and protein regions highly conserved evolutionarily. These variants are not detected in the gnomAD control cohort, and the amino acid substitutions are predicted to be damaging. Two of these variants are located in the HMG domain and abolish SOX6 transcriptional activity in vitro. No clear genotype-phenotype correlations are found. Taken together, these findings concur that SOX6 haploinsufficiency leads to a neurodevelopmental SOXopathy that often includes ADHD and abnormal skeletal and other features.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance.

PURPOSE: Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis. METHODS: Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins. RESULTS: Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis. CONCLUSION: This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

Unexpected role of SIX1 variants in craniosynostosis: expanding the phenotype of SIX1-related disorders.

BACKGROUND: Pathogenic heterozygous SIX1 variants (predominantly missense) occur in branchio-otic syndrome (BOS), but an association with craniosynostosis has not been reported. METHODS: We investigated probands with craniosynostosis of unknown cause using whole exome/genome (n=628) or RNA (n=386) sequencing, and performed targeted resequencing of SIX1 in 615 additional patients. Expression of SIX1 protein in embryonic cranial sutures was examined in the Six1nLacZ/+ reporter mouse. RESULTS: From 1629 unrelated cases with craniosynostosis we identified seven different SIX1 variants (three missense, including two de novo mutations, and four nonsense, one of which was also present in an affected twin). Compared with population data, enrichment of SIX1 loss-of-function variants was highly significant (p=0.00003). All individuals with craniosynostosis had sagittal suture fusion; additionally four had bilambdoid synostosis. Associated BOS features were often attenuated; some carrier relatives appeared non-penetrant. SIX1 is expressed in a layer basal to the calvaria, likely corresponding to the dura mater, and in the mid-sagittal mesenchyme. CONCLUSION: Craniosynostosis is associated with heterozygous SIX1 variants, with possible enrichment of loss-of-function variants compared with classical BOS. We recommend screening of SIX1 in craniosynostosis, particularly when sagittal±lambdoid synostosis and/or any BOS phenotypes are present. These findings highlight the role of SIX1 in cranial suture homeostasis.

De Novo Missense Substitutions in the Gene Encoding CDK8, a Regulator of the Mediator Complex, Cause a Syndromic Developmental Disorder.

The Mediator is an evolutionarily conserved, multi-subunit complex that regulates multiple steps of transcription. Mediator activity is regulated by the reversible association of a four-subunit module comprising CDK8 or CDK19 kinases, together with cyclin C, MED12 or MED12L, and MED13 or MED13L. Mutations in MED12, MED13, and MED13L were previously identified in syndromic developmental disorders with overlapping phenotypes. Here, we report CDK8 mutations (located at 13q12.13) that cause a phenotypically related disorder. Using whole-exome or whole-genome sequencing, and by international collaboration, we identified eight different heterozygous missense CDK8 substitutions, including 10 shown to have arisen de novo, in 12 unrelated subjects; a recurrent mutation, c.185C>T (p.Ser62Leu), was present in five individuals. All predicted substitutions localize to the ATP-binding pocket of the kinase domain. Affected individuals have overlapping phenotypes characterized by hypotonia, mild to moderate intellectual disability, behavioral disorders, and variable facial dysmorphism. Congenital heart disease occurred in six subjects; additional features present in multiple individuals included agenesis of the corpus callosum, ano-rectal malformations, seizures, and hearing or visual impairments. To evaluate the functional impact of the mutations, we measured phosphorylation at STAT1-Ser727, a known CDK8 substrate, in a CDK8 and CDK19 CRISPR double-knockout cell line transfected with wild-type (WT) or mutant CDK8 constructs. These experiments demonstrated a reduction in STAT1 phosphorylation by all mutants, in most cases to a similar extent as in a kinase-dead control. We conclude that missense mutations in CDK8 cause a developmental disorder that has phenotypic similarity to syndromes associated with mutations in other subunits of the Mediator kinase module, indicating probable overlap in pathogenic mechanisms.

The p190 RhoGAPs, ARHGAP35, and ARHGAP5 are implicated in GnRH neuronal development: Evidence from patients with idiopathic hypogonadotropic hypogonadism, zebrafish, and in vitro GAP activity assay.

PURPOSE: The study aimed to identify novel genes for idiopathic hypogonadotropic hypogonadism (IHH). METHODS: A cohort of 1387 probands with IHH underwent exome sequencing and de novo, familial, and cohort-wide investigations. Functional studies were performed on 2 p190 Rho GTPase-activating proteins (p190 RhoGAP), ARHGAP35 and ARHGAP5, which involved in vivo modeling in larval zebrafish and an in vitro p190A-GAP activity assay. RESULTS: Rare protein-truncating variants (PTVs; n = 5) and missense variants in the RhoGAP domain (n = 7) in ARHGAP35 were identified in IHH cases (rare variant enrichment: PTV [unadjusted P = 3.1E-06] and missense [adjusted P = 4.9E-03] vs controls). Zebrafish modeling using gnrh3:egfp phenotype assessment showed that mutant larvae with deficient arhgap35a, the predominant ARHGAP35 paralog in the zebrafish brain, display decreased GnRH3-GFP+ neuronal area, a readout for IHH. In vitro GAP activity studies showed that 1 rare missense variant [ARHGAP35 p.(Arg1284Trp)] had decreased GAP activity. Rare PTVs (n = 2) also were discovered in ARHGAP5, a paralog of ARHGAP35; however, arhgap5 zebrafish mutants did not display significant GnRH3-GFP+ abnormalities. CONCLUSION: This study identified ARHGAP35 as a new autosomal dominant genetic driver for IHH and ARHGAP5 as a candidate gene for IHH. These observations suggest a novel role for the p190 RhoGAP proteins in GnRH neuronal development and integrity.

Dissection of contiguous gene effects for deletions around ERF on chromosome 19.

Heterozygous intragenic loss-of-function mutations of ERF, encoding an ETS transcription factor, were previously reported to cause a novel craniosynostosis syndrome, suggesting that ERF is haploinsufficient. We describe six families harboring heterozygous deletions including, or near to, ERF, of which four were characterized by whole-genome sequencing and two by chromosomal microarray. Based on the severity of associated intellectual disability (ID), we identify three categories of ERF-associated deletions. The smallest (32 kb) and only inherited deletion included two additional centromeric genes and was not associated with ID. Three larger deletions (264-314 kb) that included at least five further centromeric genes were associated with moderate ID, suggesting that deletion of one or more of these five genes causes ID. The individual with the most severe ID had a more telomerically extending deletion, including CIC, a known ID gene. Children found to harbor ERF deletions should be referred for craniofacial assessment, to exclude occult raised intracranial pressure.

TAOK1 is associated with neurodevelopmental disorder and essential for neuronal maturation and cortical development.

Thousand and one amino-acid kinase 1 (TAOK1) is a MAP3K protein kinase, regulating different mitogen-activated protein kinase pathways, thereby modulating a multitude of processes in the cell. Given the recent finding of TAOK1 involvement in neurodevelopmental disorders (NDDs), we investigated the role of TAOK1 in neuronal function and collected a cohort of 23 individuals with mostly de novo variants in TAOK1 to further define the associated NDD. Here, we provide evidence for an important role for TAOK1 in neuronal function, showing that altered TAOK1 expression levels in the embryonic mouse brain affect neural migration in vivo, as well as neuronal maturation in vitro. The molecular spectrum of the identified TAOK1 variants comprises largely truncating and nonsense variants, but also missense variants, for which we provide evidence that they can have a loss of function or dominant-negative effect on TAOK1, expanding the potential underlying causative mechanisms resulting in NDD. Taken together, our data indicate that TAOK1 activity needs to be properly controlled for normal neuronal function and that TAOK1 dysregulation leads to a neurodevelopmental disorder mainly comprising similar facial features, developmental delay/intellectual disability and/or variable learning or behavioral problems, muscular hypotonia, infant feeding difficulties, and growth problems.

A de novo substitution in BCL11B leads to loss of interaction with transcriptional complexes and craniosynostosis.

Craniosynostosis, the premature ossification of cranial sutures, is a developmental disorder of the skull vault, occurring in approximately 1 in 2250 births. The causes are heterogeneous, with a monogenic basis identified in ~25% of patients. Using whole-genome sequencing, we identified a novel, de novo variant in BCL11B, c.7C>A, encoding an R3S substitution (p.R3S), in a male patient with coronal suture synostosis. BCL11B is a transcription factor that interacts directly with the nucleosome remodelling and deacetylation complex (NuRD) and polycomb-related complex 2 (PRC2) through the invariant proteins RBBP4 and RBBP7. The p.R3S substitution occurs within a conserved amino-terminal motif (RRKQxxP) of BCL11B and reduces interaction with both transcriptional complexes. Equilibrium binding studies and molecular dynamics simulations show that the p.R3S substitution disrupts ionic coordination between BCL11B and the RBBP4-MTA1 complex, a subassembly of the NuRD complex, and increases the conformational flexibility of Arg-4, Lys-5 and Gln-6 of BCL11B. These alterations collectively reduce the affinity of BCL11B p.R3S for the RBBP4-MTA1 complex by nearly an order of magnitude. We generated a mouse model of the BCL11B p.R3S substitution using a CRISPR-Cas9-based approach, and we report herein that these mice exhibit craniosynostosis of the coronal suture, as well as other cranial sutures. This finding provides strong evidence that the BCL11B p.R3S substitution is causally associated with craniosynostosis and confirms an important role for BCL11B in the maintenance of cranial suture patency.

Mutations in the BAF-Complex Subunit DPF2 Are Associated with Coffin-Siris Syndrome.

Variants affecting the function of different subunits of the BAF chromatin-remodelling complex lead to various neurodevelopmental syndromes, including Coffin-Siris syndrome. Furthermore, variants in proteins containing PHD fingers, motifs recognizing specific histone tail modifications, have been associated with several neurological and developmental-delay disorders. Here, we report eight heterozygous de novo variants (one frameshift, two splice site, and five missense) in the gene encoding the BAF complex subunit double plant homeodomain finger 2 (DPF2). Affected individuals share common clinical features described in individuals with Coffin-Siris syndrome, including coarse facial features, global developmental delay, intellectual disability, speech impairment, and hypoplasia of fingernails and toenails. All variants occur within the highly conserved PHD1 and PHD2 motifs. Moreover, missense variants are situated close to zinc binding sites and are predicted to disrupt these sites. Pull-down assays of recombinant proteins and histone peptides revealed that a subset of the identified missense variants abolish or impaire DPF2 binding to unmodified and modified H3 histone tails. These results suggest an impairment of PHD finger structural integrity and cohesion and most likely an aberrant recognition of histone modifications. Furthermore, the overexpression of these variants in HEK293 and COS7 cell lines was associated with the formation of nuclear aggregates and the recruitment of both wild-type DPF2 and BRG1 to these aggregates. Expression analysis of truncating variants found in the affected individuals indicated that the aberrant transcripts escape nonsense-mediated decay. Altogether, we provide compelling evidence that de novo variants in DPF2 cause Coffin-Siris syndrome and propose a dominant-negative mechanism of pathogenicity.

De Novo Variants in the F-Box Protein FBXO11 in 20 Individuals with a Variable Neurodevelopmental Disorder.

Next-generation sequencing combined with international data sharing has enormously facilitated identification of new disease-associated genes and mutations. This is particularly true for genetically extremely heterogeneous entities such as neurodevelopmental disorders (NDDs). Through exome sequencing and world-wide collaborations, we identified and assembled 20 individuals with de novo variants in FBXO11. They present with mild to severe developmental delay associated with a range of features including short (4/20) or tall (2/20) stature, obesity (5/20), microcephaly (4/19) or macrocephaly (2/19), behavioral problems (17/20), seizures (5/20), cleft lip or palate or bifid uvula (3/20), and minor skeletal anomalies. FBXO11 encodes a member of the F-Box protein family, constituting a subunit of an E3-ubiquitin ligase complex. This complex is involved in ubiquitination and proteasomal degradation and thus in controlling critical biological processes by regulating protein turnover. The identified de novo aberrations comprise two large deletions, ten likely gene disrupting variants, and eight missense variants distributed throughout FBXO11. Structural modeling for missense variants located in the CASH or the Zinc-finger UBR domains suggests destabilization of the protein. This, in combination with the observed spectrum and localization of identified variants and the lack of apparent genotype-phenotype correlations, is compatible with loss of function or haploinsufficiency as an underlying mechanism. We implicate de novo missense and likely gene disrupting variants in FBXO11 in a neurodevelopmental disorder with variable intellectual disability and various other features.

Evaluating the performance of a clinical genome sequencing program for diagnosis of rare genetic disease, seen through the lens of craniosynostosis.

PURPOSE: Genome sequencing (GS) for diagnosis of rare genetic disease is being introduced into the clinic, but the complexity of the data poses challenges for developing pipelines with high diagnostic sensitivity. We evaluated the performance of the Genomics England 100,000 Genomes Project (100kGP) panel-based pipelines, using craniosynostosis as a test disease. METHODS: GS data from 114 probands with craniosynostosis and their relatives (314 samples), negative on routine genetic testing, were scrutinized by a specialized research team, and diagnoses compared with those made by 100kGP. RESULTS: Sixteen likely pathogenic/pathogenic variants were identified by 100kGP. Eighteen additional likely pathogenic/pathogenic variants were identified by the research team, indicating that for craniosynostosis, 100kGP panels had a diagnostic sensitivity of only 47%. Measures that could have augmented diagnoses were improved calling of existing panel genes (+18% sensitivity), review of updated panels (+12%), comprehensive analysis of de novo small variants (+29%), and copy-number/structural variants (+9%). Recent NHS England recommendations that partially incorporate these measures should achieve 85% overall sensitivity (+38%). CONCLUSION: GS identified likely pathogenic/pathogenic variants in 29.8% of previously undiagnosed patients with craniosynostosis. This demonstrates the value of research analysis and the importance of continually improving algorithms to maximize the potential of clinical GS.

amplimap: a versatile tool to process and analyze targeted NGS data.

SUMMARY: amplimap is a command-line tool to automate the processing and analysis of data from targeted next-generation sequencing experiments with PCR-based amplicons or capture-based enrichment systems. From raw sequencing reads, amplimap generates output such as read alignments, annotated variant calls, target coverage statistics and variant allele counts and frequencies for each target base pair. In addition to its focus on user-friendliness and reproducibility, amplimap supports advanced features such as consensus base calling for read families based on unique molecular identifiers and filtering false positive variant calls caused by amplification of off-target loci. AVAILABILITY AND IMPLEMENTATION: amplimap is available as a free Python package under the open-source Apache 2.0 License. Documentation, source code and installation instructions are available at https://github.com/koelling/amplimap.

Correction: SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism nearBMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure.

Identification of mobile retrocopies during genetic testing: Consequences for routine diagnosis.

Human retrocopies, that is messenger RNA transcripts benefitting from the long interspersed element 1 machinery for retrotransposition, may have specific consequences for genomic testing. Next genetration sequencing (NGS) techniques allow the detection of such mobile elements but they may be misinterpreted as genomic duplications or be totally overlooked. We report eight observations of retrocopies detected during diagnostic NGS analyses of targeted gene panels, exome, or genome sequencing. For seven cases, while an exons-only copy number gain was called, read alignment inspection revealed a depth of coverage shift at every exon-intron junction where indels were also systematically called. Moreover, aberrant chimeric read pairs spanned entire introns or were paired with another locus for terminal exons. The 8th retrocopy was present in the reference genome and thus showed a normal NGS profile. We emphasize the existence of retrocopies and strategies to accurately detect them at a glance during genetic testing and discuss pitfalls for genetic testing.

Disruption of TWIST1 translation by 5' UTR variants in Saethre-Chotzen syndrome.

Saethre-Chotzen syndrome (SCS), one of the most common forms of syndromic craniosynostosis (premature fusion of the cranial sutures), results from haploinsufficiency of TWIST1, caused by deletions of the entire gene or loss-of-function variants within the coding region. To determine whether non-coding variants also contribute to SCS, we screened 14 genetically undiagnosed SCS patients using targeted capture sequencing, and identified novel single nucleotide variants (SNVs) in the 5' untranslated region (UTR) of TWIST1 in two unrelated SCS cases. We show experimentally that these variants, which create translation start sites in the TWIST1 leader sequence, reduce translation from the main open reading frame (mORF). This is the first demonstration that non-coding SNVs of TWIST1 can cause SCS, and highlights the importance of screening the 5' UTR in clinically diagnosed SCS patients without a coding mutation. Similar 5' UTR variants, particularly of haploinsufficient genes, may represent an under-ascertained cause of monogenic disease.

Junctophilin-1 is a modifier gene of GDAP1-related Charcot-Marie-Tooth disease.

Mutations in the GDAP1 gene cause different forms of Charcot-Marie-Tooth (CMT) disease, and the primary clinical expression of this disease is markedly variable in the dominant inheritance form (CMT type 2K; CMT2K), in which carriers of the GDAP1 p.R120W mutation can display a wide range of clinical severity. We investigated the JPH1 gene as a genetic modifier of clinical expression variability because junctophilin-1 (JPH1) is a good positional and functional candidate. We demonstrated that the JPH1-GDAP1 cluster forms a paralogon and is conserved in vertebrates. Moreover, both proteins play a role in Ca(2+) homeostasis, and we demonstrated that JPH1 is able to restore the store-operated Ca(2+) entry (SOCE) activity in GDAP1-silenced cells. After the mutational screening of JPH1 in a series of 24 CMT2K subjects who harbour the GDAP1 p.R120W mutation, we characterized the JPH1 p.R213P mutation in one patient with a more severe clinical picture. JPH1(p.R213P) cannot rescue the SOCE response in GDAP1-silenced cells. We observed that JPH1 colocalizes with STIM1, which is the activator of SOCE, in endoplasmic reticulum-plasma membrane puncta structures during Ca(2+) release in a GDAP1-dependent manner. However, when GDAP1(p.R120W) is expressed, JPH1 seems to be retained in mitochondria. We also established that the combination of GDAP1(p.R120W) and JPH1(p.R213P) dramatically reduces SOCE activity, mimicking the effect observed in GDAP1 knock-down cells. In summary, we conclude that JPH1 and GDAP1 share a common pathway and depend on each other; therefore, JPH1 can contribute to the phenotypical consequences of GDAP1 mutations.