A Comparison of Structural Variant Calling from Short-Read and Nanopore-Based Whole-Genome Sequencing Using Optical Genome Mapping as a Benchmark.

The identification of structural variants (SVs) in genomic data represents an ongoing challenge because of difficulties in reliable SV calling leading to reduced sensitivity and specificity. We prepared high-quality DNA from 9 parent-child trios, who had previously undergone short-read whole-genome sequencing (Illumina platform) as part of the Genomics England 100,000 Genomes Project. We reanalysed the genomes using both Bionano optical genome mapping (OGM; 8 probands and one trio) and Nanopore long-read sequencing (Oxford Nanopore Technologies [ONT] platform; all samples). To establish a "truth" dataset, we asked whether rare proband SV calls (n = 234) made by the Bionano Access (version 1.6.1)/Solve software (version 3.6.1_11162020) could be verified by individual visualisation using the Integrative Genomics Viewer with either or both of the Illumina and ONT raw sequence. Of these, 222 calls were verified, indicating that Bionano OGM calls have high precision (positive predictive value 95%). We then asked what proportion of the 222 true Bionano SVs had been identified by SV callers in the other two datasets. In the Illumina dataset, sensitivity varied according to variant type, being high for deletions (115/134; 86%) but poor for insertions (13/58; 22%). In the ONT dataset, sensitivity was generally poor using the original Sniffles variant caller (48% overall) but improved substantially with use of Sniffles2 (36/40; 90% and 17/23; 74% for deletions and insertions, respectively). In summary, we show that the precision of OGM is very high. In addition, when applying the Sniffles2 caller, the sensitivity of SV calling using ONT long-read sequence data outperforms Illumina sequencing for most SV types.

Reassessing the association: Evaluation of a polyalanine deletion variant of RUNX2 in non-syndromic sagittal and metopic craniosynostosis.

The RUNT-related transcription factor RUNX2 plays a critical role in osteoblast differentiation, and alterations to gene dosage cause distinct craniofacial anomalies. Uniquely amongst the RUNT-related family, vertebrate RUNX2 encodes a polyglutamine/polyalanine repeat (Gln23-Glu-Ala17 in humans), with the length of the polyalanine component completely conserved in great apes. Surprisingly, a frequent 6-amino acid deletion polymorphism, p.(Ala84_Ala89)del, occurs in humans (termed 11A allele), and a previous association study (Cuellar et al. Bone 137:115395;2020) reported that the 11A variant was significantly more frequent in non-syndromic sagittal craniosynostosis (nsSag; allele frequency [AF] = 0.156; 95% confidence interval [CI] 0.126-0.189) compared to non-syndromic metopic craniosynostosis (nsMet; AF = 0.068; 95% CI 0.045-0.098). However, the gnomAD v.2.1.1 control population used by Cuellar et al. did not display Hardy-Weinberg equilibrium, hampering interpretation. To re-examine this association, we genotyped the RUNX2 11A polymorphism in 225 individuals with sporadic nsSag as parent-child trios and 164 singletons with sporadic nsMet, restricting our analysis to individuals of European ancestry. We compared observed allele frequencies to the non-transmitted alleles in the parent-child trios, and to the genome sequencing data from gnomAD v.4, which display Hardy-Weinberg equilibrium. Observed AFs (and 95% CI) were 0.076 (0.053-0.104) in nsSag and 0.082 (0.055-0.118) in nsMet, compared with 0.062 (0.042-0.089) in non-transmitted parental alleles and 0.065 (0.063-0.067) in gnomAD v.4.0.0 non-Finnish European control genomes. In summary, we observed a non-significant excess, compared to gnomAD data, of 11A alleles in both nsSag (relative risk 1.18, 95% CI 0.83-1.67) and nsMet (relative risk 1.29, 95% CI 0.87-1.92), but we did not replicate the much higher excess of RUNX2 11A alleles in nsSag previously reported (p = 0.0001).

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.

Review of Recurrently Mutated Genes in Craniosynostosis Supports Expansion of Diagnostic Gene Panels.

Craniosynostosis, the premature fusion of the cranial sutures, affects ~1 in 2000 children. Although many patients with a genetically determined cause harbor a variant in one of just seven genes or have a chromosomal abnormality, over 60 genes are known to be recurrently mutated, thus comprising a long tail of rarer diagnoses. Genome sequencing for the diagnosis of rare diseases is increasingly used in clinical settings, but analysis of the data is labor intensive and involves a trade-off between achieving high sensitivity or high precision. PanelApp, a crowd-sourced disease-focused set of gene panels, was designed to enable prioritization of variants in known disease genes for a given pathology, allowing enhanced identification of true-positives. For heterogeneous disorders like craniosynostosis, these panels must be regularly updated to ensure that diagnoses are not being missed. We provide a systematic review of genetic literature on craniosynostosis over the last 5 years, including additional results from resequencing a 42-gene panel in 617 affected individuals. We identify 16 genes (representing a 25% uplift) that should be added to the list of bona fide craniosynostosis disease genes and discuss the insights that these new genes provide into pathophysiological mechanisms of craniosynostosis.

Obesity-Associated GNAS Mutations and the Melanocortin Pathway.

BACKGROUND: GNAS encodes the Gαs (stimulatory G-protein alpha subunit) protein, which mediates G protein-coupled receptor (GPCR) signaling. GNAS mutations cause developmental delay, short stature, and skeletal abnormalities in a syndrome called Albright's hereditary osteodystrophy. Because of imprinting, mutations on the maternal allele also cause obesity and hormone resistance (pseudohypoparathyroidism). METHODS: We performed exome sequencing and targeted resequencing in 2548 children who presented with severe obesity, and we unexpectedly identified 22 GNAS mutation carriers. We investigated whether the effect of GNAS mutations on melanocortin 4 receptor (MC4R) signaling explains the obesity and whether the variable clinical spectrum in patients might be explained by the results of molecular assays. RESULTS: Almost all GNAS mutations impaired MC4R signaling. A total of 6 of 11 patients who were 12 to 18 years of age had reduced growth. In these patients, mutations disrupted growth hormone-releasing hormone receptor signaling, but growth was unaffected in carriers of mutations that did not affect this signaling pathway (mean standard-deviation score for height, -0.90 vs. 0.75, respectively; P = 0.02). Only 1 of 10 patients who reached final height before or during the study had short stature. GNAS mutations that impaired thyrotropin receptor signaling were associated with developmental delay and with higher thyrotropin levels (mean [±SD], 8.4±4.7 mIU per liter) than those in 340 severely obese children who did not have GNAS mutations (3.9±2.6 mIU per liter; P = 0.004). CONCLUSIONS: Because pathogenic mutations may manifest with obesity alone, screening of children with severe obesity for GNAS deficiency may allow early diagnosis, improving clinical outcomes, and melanocortin agonists may aid in weight loss. GNAS mutations that are identified by means of unbiased genetic testing differentially affect GPCR signaling pathways that contribute to clinical heterogeneity. Monogenic diseases are clinically more variable than their classic descriptions suggest. (Funded by Wellcome and others.).

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.

Biallelic P4HTM variants associated with HIDEA syndrome and mitochondrial respiratory chain complex I deficiency.

We report a patient with profound congenital hypotonia, central hypoventilation, poor visual behaviour with retinal hypopigmentation, and significantly decreased mitochondrial respiratory chain complex I activity in muscle, who died at 7 months of age having made minimal developmental progress. Biallelic predicted truncating P4HTM variants were identified following trio whole-genome sequencing, consistent with a diagnosis of hypotonia, hypoventilation, intellectual disability, dysautonomia, epilepsy and eye abnormalities (HIDEA) syndrome. Very few patients with HIDEA syndrome have been reported previously and mitochondrial abnormalities were observed in three of four previous cases who had a muscle biopsy, suggesting the possibility that HIDEA syndrome represents a primary mitochondrial disorder. P4HTM encodes a transmembrane prolyl 4-hydroxylase with putative targets including hypoxia inducible factors, RNA polymerase II and activating transcription factor 4, which has been implicated in the integrated stress response observed in cell and animal models of mitochondrial disease, and may explain the mitochondrial dysfunction observed in HIDEA syndrome.

Major skull manifestations of skeletal dysplasias - pictorial essay.

Skeletal dysplasias are a large group of rare conditions with widely heterogeneous manifestations and a reputation for being diagnostically difficult. Involvement of the brain and craniovertebral junction are features familiar to the paediatric neuroradiologist. Involvement of the skull itself represents an area of overlap between the domains of the neuroradiologist and the skeletal dysplasia radiologist. In this pictorial essay, we review the principal skull manifestations of skeletal dysplasias as they present to the neuroradiologist.

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.

A genome-wide association study implicates the BMP7 locus as a risk factor for nonsyndromic metopic craniosynostosis.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10-8): rs781716 (P = 4.71 × 10-9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10-8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10-9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10-8, OR = 0.45; P = 3.31 × 10-8, OR = 0.45; P = 1.09 × 10-8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10-8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10-6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821-55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.

Cenani-Lenz syndactyly in siblings with a novel homozygous LRP4 mutation and recurrent hypoglycaemia.

Cenani-Lenz syndactyly (CLS) is a rare autosomal recessive syndrome characterized by disorganized oligosyndactyly of upper and lower limbs as well as radioulnar synostosis. Structural renal abnormalities are also common. We report two affected brothers, born to orthodox Jewish parents, in whom we found a novel homozygous missense variant c.4910G>A; p.(Cys1637Tyr) in LRP4 situated in an EGF-like domain between the fourth beta-propeller and transmembrane domains. Both brothers have had recurrent ketotic hypoglycaemia which has not been associated previously. We present 3D computed tomographic imaging illustrating the limb abnormalities in detail.

Juvenile Paget's disease with compound heterozygous mutations in TNFRSF11B presenting with recurrent clavicular fractures and a mild skeletal phenotype.

Juvenile Paget's disease (JPD) is a rare recessively-inherited bone dysplasia. The great majority of cases described to date have had homozygous mutations in TNFRSF11B, the gene encoding osteoprotegerin. We describe a boy who presented with recurrent clavicular fractures following minor trauma (8 fractures from age 2 to 11). He was of normal height and despite mild lateral bowing of the thighs and anterior bowing of the shins he remained physically active. Abnormal modelling was noted in ribs and humeri on clavicular radiographs, and a skeletal survey at the age of 7 showed generalised diaphyseal expansion of the long bones with thickening of the periosteal and endosteal surfaces of the cortices. On biochemical evaluation, serum alkaline phosphatase was noted to be persistently elevated. The diagnosis of JPD was confirmed by the finding of compound heterozygous mutations in TNFRSF11B: a maternally-inherited A>G missense mutation at position 1 of the first amino acid codon (previously reported) and a paternally-inherited splice acceptor site mutation in intron 3 at a highly conserved position (not previously reported). Bioinformatics analysis suggested both mutations were disease-causing. Compound heterozygote mutations in TNFRSF11B causing JPD have been previously reported only once - in a boy who also had a relatively mild skeletal phenotype. The milder features may lead to delay in diagnosis and diagnostic confusion with other entities, but the extraskeletal features of JPD may nonetheless develop.

Biallelic human ITCH variants causing a multisystem disease with dysmorphic features: A second report.

We report a 23 year old female with biallelic truncating variants in the ITCH (Itchy E3 Ubiquitin protein ligase, mouse homolog of; OMIM60649) gene associated with marked short stature, severe early onset chronic lung disease resembling asthma, dysmorphic facial features, and symmetrical camptodactyly of the fingers but normal intellect. The condition has only been reported once previously (Lohr et al., American Journal of Human Genetics, 2010, 86, 447-453) in 10 children from an Old Order Amish family found to have a homozygous frameshift truncating variant in association with failure to thrive, chronic lung disease, motor and cognitive delay, and variable autoimmune diseases including autoimmune hepatitis, enteropathy, hypothyroidism, and diabetes. The condition is listed in OMIM as Autoimmune disease, Multisystem with Facial Dysmorphism (OMIM613385). The clinical course as well as the dysmorphic facial and limb features overlap closely with our patient. We believe the triad of marked syndromic short stature, chronic lung disease, and dysmorphism (with or without cognitive impairment and wider autoimmune involvement) is distinctive.

ERF-related craniosynostosis: The phenotypic and developmental profile of a new craniosynostosis syndrome.

Mutations in the ERF gene, coding for ETS2 repressor factor, a member of the ETS family of transcription factors cause a recently recognized syndromic form of craniosynostosis (CRS4) with facial dysmorphism, Chiari-1 malformation, speech and language delay, and learning difficulties and/or behavioral problems. The overall prevalence of ERF mutations in patients with syndromic craniosynostosis is around 2%, and 0.7% in clinically nonsyndromic craniosynostosis. Here, we present findings from 16 unrelated probands with ERF-related craniosynostosis, with additional data from 20 family members sharing the mutations. Most of the probands exhibited multisutural (including pan-) synostosis but a pattern involving the sagittal and lambdoid sutures (Mercedes-Benz pattern) predominated. Importantly the craniosynostosis was often postnatal in onset, insidious and progressive with subtle effects on head morphology resulting in a median age at presentation of 42 months among the probands and, in some instances, permanent visual impairment due to unsuspected raised intracranial pressure (ICP). Facial dysmorphism (exhibited by all of the probands and many of the affected relatives) took the form of orbital hypertelorism, mild exorbitism and malar hypoplasia resembling Crouzon syndrome but, importantly, a Class I occlusal relationship. Speech delay, poor gross and/or fine motor control, hyperactivity and poor concentration were common. Cranial vault surgery for raised ICP and/or Chiari-1 malformation was expected when multisutural synostosis was observed. Variable expressivity and nonpenetrance among genetically affected relatives was encountered. These observations form the most complete phenotypic and developmental profile of this recently identified craniosynostosis syndrome yet described and have important implications for surgical intervention and follow-up.

A New Multisystem Disorder Caused by the Gαs Mutation p.F376V.

CONTEXT: The α subunit of the stimulatory G protein (Gαs) links numerous receptors to adenylyl cyclase. Gαs, encoded by GNAS, is expressed predominantly from the maternal allele in certain tissues. Thus, maternal heterozygous loss-of-function mutations cause hormonal resistance, as in pseudohypoparathyroidism type Ia, whereas somatic gain-of-function mutations cause hormone-independent endocrine stimulation, as in McCune-Albright syndrome. OBJECTIVE: We report two unrelated boys presenting with a new combination of clinical findings that suggest both gain and loss of Gαs function. DESIGN AND SETTING: Clinical features were studied and sequencing of GNAS was performed. Signaling capacities of wild-type and mutant Gαs were determined in the presence of different G protein-coupled receptors (GPCRs) under basal and agonist-stimulated conditions. RESULTS: Both unrelated patients presented with unexplained hyponatremia in infancy, followed by severe early onset gonadotrophin-independent precocious puberty and skeletal abnormalities. An identical heterozygous de novo variant (c.1136T>G; p.F376V) was found on the maternal GNAS allele in both patients; this resulted in a clinical phenotype that differed from known Gαs-related diseases and suggested gain of function at the vasopressin 2 receptor (V2R) and lutropin/choriogonadotropin receptor (LHCGR), yet increased serum PTH concentrations indicative of impaired proximal tubular PTH1 receptor (PTH1R) function. In vitro studies demonstrated that Gαs-F376V enhanced ligand-independent signaling at the PTH1R, LHCGR, and V2R and, at the same time, blunted ligand-dependent responses. Structural homology modeling suggested mutation-induced modifications at the C-terminal α5 helix of Gαs that are relevant for interaction with GPCRs and signal transduction. CONCLUSIONS: The Gαs p.F376V mutation causes a previously unrecognized multisystem disorder.

Natural history and genotype-phenotype correlations in 72 individuals with SATB2-associated syndrome.

SATB2-associated syndrome (SAS) is an autosomal dominant disorder characterized by significant neurodevelopmental disabilities with limited to absent speech, behavioral issues, and craniofacial anomalies. Previous studies have largely been restricted to case reports and small series without in-depth phenotypic characterization or genotype-phenotype correlations. Seventy two study participants were identified as part of the SAS clinical registry. Individuals with a molecularly confirmed diagnosis of SAS were referred after clinical diagnostic testing. In this series we present the most comprehensive phenotypic and genotypic characterization of SAS to date, including prevalence of each clinical feature, neurodevelopmental milestones, and when available, patient management. We confirm that the most distinctive features are neurodevelopmental delay with invariably severely limited speech, abnormalities of the palate (cleft or high-arched), dental anomalies (crowding, macrodontia, abnormal shape), and behavioral issues with or without bone or brain anomalies. This comprehensive clinical characterization will help clinicians with the diagnosis, counseling and management of SAS and help provide families with anticipatory guidance.

Autosomal dominant frontometaphyseal dysplasia: Delineation of the clinical phenotype.

Frontometaphyseal dysplasia (FMD) is caused by gain-of-function mutations in the X-linked gene FLNA in approximately 50% of patients. Recently we characterized an autosomal dominant form of FMD (AD-FMD) caused by mutations in MAP3K7, which accounts for the condition in the majority of patients who lack a FLNA mutation. We previously also described a patient with a de novo variant in TAB2, which we hypothesized was causative of another form of AD-FMD. In this study, a cohort of 20 individuals with AD-FMD is clinically evaluated. This cohort consists of 15 individuals with the recently described, recurrent mutation (c.1454C>T) in MAP3K7, as well as three individuals with missense mutations that result in substitutions in the N-terminal kinase domain of TGFß-activated kinase 1 (TAK1), encoded by MAP3K7. Additionally, two individuals have missense variants in the gene TAB2, which encodes a protein with a close functional relationship to TAK1, TAK1-associated binding protein 2 (TAB2). Although the X-linked and autosomal dominant forms of FMD are very similar, there are distinctions to be made between the two conditions. Individuals with AD-FMD have characteristic facial features, and are more likely to be deaf, have scoliosis and cervical fusions, and have a cleft palate. Furthermore, there are features only found in AD-FMD in our review of the literature including valgus deformity of the feet and predisposition to keloid scarring. Finally, intellectual disability is present in a small number of subjects with AD-FMD but has not been described in association with X-linked FMD.

Diagnostic value of exome and whole genome sequencing in craniosynostosis.

BACKGROUND: Craniosynostosis, the premature fusion of one or more cranial sutures, occurs in ∼1 in 2250 births, either in isolation or as part of a syndrome. Mutations in at least 57 genes have been associated with craniosynostosis, but only a minority of these are included in routine laboratory genetic testing. METHODS: We used exome or whole genome sequencing to seek a genetic cause in a cohort of 40 subjects with craniosynostosis, selected by clinical or molecular geneticists as being high-priority cases, and in whom prior clinically driven genetic testing had been negative. RESULTS: We identified likely associated mutations in 15 patients (37.5%), involving 14 different genes. All genes were mutated in single families, except for IL11RA (two families). We classified the other positive diagnoses as follows: commonly mutated craniosynostosis genes with atypical presentation (EFNB1, TWIST1); other core craniosynostosis genes (CDC45, MSX2, ZIC1); genes for which mutations are only rarely associated with craniosynostosis (FBN1, HUWE1, KRAS, STAT3); and known disease genes for which a causal relationship with craniosynostosis is currently unknown (AHDC1, NTRK2). In two further families, likely novel disease genes are currently undergoing functional validation. In 5 of the 15 positive cases, the (previously unanticipated) molecular diagnosis had immediate, actionable consequences for either genetic or medical management (mutations in EFNB1, FBN1, KRAS, NTRK2, STAT3). CONCLUSIONS: This substantial genetic heterogeneity, and the multiple actionable mutations identified, emphasises the benefits of exome/whole genome sequencing to identify causal mutations in craniosynostosis cases for which routine clinical testing has yielded negative results.

Mutations in PROSC Disrupt Cellular Pyridoxal Phosphate Homeostasis and Cause Vitamin-B6-Dependent Epilepsy.

Pyridoxal 5'-phosphate (PLP), the active form of vitamin B6, functions as a cofactor in humans for more than 140 enzymes, many of which are involved in neurotransmitter synthesis and degradation. A deficiency of PLP can present, therefore, as seizures and other symptoms that are treatable with PLP and/or pyridoxine. Deficiency of PLP in the brain can be caused by inborn errors affecting B6 vitamer metabolism or by inactivation of PLP, which can occur when compounds accumulate as a result of inborn errors of other pathways or when small molecules are ingested. Whole-exome sequencing of two children from a consanguineous family with pyridoxine-dependent epilepsy revealed a homozygous nonsense mutation in proline synthetase co-transcribed homolog (bacterial), PROSC, which encodes a PLP-binding protein of hitherto unknown function. Subsequent sequencing of 29 unrelated indivduals with pyridoxine-responsive epilepsy identified four additional children with biallelic PROSC mutations. Pre-treatment cerebrospinal fluid samples showed low PLP concentrations and evidence of reduced activity of PLP-dependent enzymes. However, cultured fibroblasts showed excessive PLP accumulation. An E.coli mutant lacking the PROSC homolog (ΔYggS) is pyridoxine sensitive; complementation with human PROSC restored growth whereas hPROSC encoding p.Leu175Pro, p.Arg241Gln, and p.Ser78Ter did not. PLP, a highly reactive aldehyde, poses a problem for cells, which is how to supply enough PLP for apoenzymes while maintaining free PLP concentrations low enough to avoid unwanted reactions with other important cellular nucleophiles. Although the mechanism involved is not fully understood, our studies suggest that PROSC is involved in intracellular homeostatic regulation of PLP, supplying this cofactor to apoenzymes while minimizing any toxic side reactions.