Familial severe skeletal Class II malocclusion with gingival hyperplasia caused by a complex structural rearrangement at the KCNJ2-KCNJ16 locus.

The aim of this work was to identify the underlying genetic cause in a four-generation family segregating an unusual phenotype comprising a severe form of skeletal Class II malocclusion with gingival hyperplasia. SNP array identified a copy number gain on chromosome 1 (chr1); however, this chromosomal region did not segregate correctly in the extended family. Exome sequencing also failed to identify a candidate causative variant but highlighted co-segregating genetic markers on chr17 and chr19. Short- and long-read genome sequencing allowed us to pinpoint and characterize at nucleotide-level resolution a chromothripsis-like complex rearrangement (CR) inserted into the chr17 co-segregating region at the KCNJ2-SOX9 locus. The CR involved the gain of five different regions from chr1 that are shuffled, chained, and inserted as a single block (∼828 kb) at chr17q24.3. The inserted sequences contain craniofacial enhancers that are predicted to interact with KCNJ2/KCNJ16 through neo-topologically associating domain (TAD) formation to induce ectopic activation. Our findings suggest that the CR inserted at chr17q24.3 is the cause of the severe skeletal Class II malocclusion with gingival hyperplasia in this family and expands the panoply of phenotypes linked to variation at the KCNJ2-SOX9 locus. In addition, we highlight a previously overlooked potential role for misregulation of the KCNJ2/KCNJ16 genes in the pathomechanism of gingival hyperplasia associated with deletions and other rearrangements of the 17q24.2-q24.3 region (MIM 135400).

BTB domain mutations perturbing KCTD15 oligomerisation cause a distinctive frontonasal dysplasia syndrome.

INTRODUCTION: KCTD15 encodes an oligomeric BTB domain protein reported to inhibit neural crest formation through repression of Wnt/beta-catenin signalling, as well as transactivation by TFAP2. Heterozygous missense variants in the closely related paralogue KCTD1 cause scalp-ear-nipple syndrome. METHODS: Exome sequencing was performed on a two-generation family affected by a distinctive phenotype comprising a lipomatous frontonasal malformation, anosmia, cutis aplasia of the scalp and/or sparse hair, and congenital heart disease. Identification of a de novo missense substitution within KCTD15 led to targeted sequencing of DNA from a similarly affected sporadic patient, revealing a different missense mutation. Structural and biophysical analyses were performed to assess the effects of both amino acid substitutions on the KCTD15 protein. RESULTS: A heterozygous c.310G>C variant encoding p.(Asp104His) within the BTB domain of KCTD15 was identified in an affected father and daughter and segregated with the phenotype. In the sporadically affected patient, a de novo heterozygous c.263G>A variant encoding p.(Gly88Asp) was present in KCTD15. Both substitutions were found to perturb the pentameric assembly of the BTB domain. A crystal structure of the BTB domain variant p.(Gly88Asp) revealed a closed hexameric assembly, whereas biophysical analyses showed that the p.(Asp104His) substitution resulted in a monomeric BTB domain likely to be partially unfolded at physiological temperatures. CONCLUSION: BTB domain substitutions in KCTD1 and KCTD15 cause clinically overlapping phenotypes involving craniofacial abnormalities and cutis aplasia. The structural analyses demonstrate that missense substitutions act through a dominant negative mechanism by disrupting the higher order structure of the KCTD15 protein complex.

Germinal centers output clonally diverse plasma cell populations expressing high- and low-affinity antibodies.

Germinal centers (GCs) form in lymph nodes after immunization or infection to facilitate antibody affinity maturation and memory and plasma cell (PC) development. PC differentiation is thought to involve stringent selection for GC B cells expressing the highest-affinity antigen receptors, but how this plays out during complex polyclonal responses is unclear. We combine temporal lineage tracing with antibody characterization to gain a snapshot of PCs developing during influenza infection. GCs co-mature B cell clones with antibody affinities spanning multiple orders of magnitude; however, each generates PCs with similar efficiencies, including weak binders. Within lineages, PC selection is not restricted to variants with the highest-affinity antibodies. Differentiation is commonly associated with proliferative expansion to produce "nodes" of identical PCs. Immunization-induced GCs generate fewer PCs but still of low- and high-antibody affinities. We propose that generating low-affinity antibody PCs reflects an evolutionary compromise to facilitate diverse serum antibody responses.

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.

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.

CSynth: an interactive modelling and visualization tool for 3D chromatin structure.

MOTIVATION: The 3D structure of chromatin in the nucleus is important for gene expression and regulation. Chromosome conformation capture techniques, such as Hi-C, generate large amounts of data showing interaction points on the genome but these are hard to interpret using standard tools. RESULTS: We have developed CSynth, an interactive 3D genome browser and real-time chromatin restraint-based modeller to visualize models of any chromosome conformation capture (3C) data. Unlike other modelling systems, CSynth allows dynamic interaction with the modelling parameters to allow experimentation and effects on the model. It also allows comparison of models generated from data in different tissues/cell states and the results of third-party 3D modelling outputs. In addition, we include an option to view and manipulate these complicated structures using Virtual Reality (VR) so scientists can immerse themselves in the models for further understanding. This VR component has also proven to be a valuable teaching and a public engagement tool. AVAILABILITYAND IMPLEMENTATION: CSynth is web based and available to use at csynth.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Germinal Center B Cells Replace Their Antigen Receptors in Dark Zones and Fail Light Zone Entry when Immunoglobulin Gene Mutations are Damaging.

Adaptive immunity involves the development of bespoke antibodies in germinal centers (GCs) through immunoglobulin somatic hypermutation (SHM) in GC dark zones (DZs) and clonal selection in light zones (LZs). Accurate selection requires that cells fully replace surface B cell receptors (BCRs) following SHM, but whether this happens before LZ entry is not clear. We found that most GC B cells degrade pre-SHM receptors before leaving the DZ, and that B cells acquiring crippling mutations during SHM rarely reached the LZ. Instead, apoptosis was triggered preferentially in late G1, a stage wherein cells with functional BCRs re-entered cell cycle or reduced surface expression of the chemokine receptor CXCR4 to enable LZ migration. Ectopic expression of the anti-apoptotic gene Bcl2 was not sufficient for cells with damaging mutations to reach the LZ, suggesting that BCR-dependent cues may actively facilitate the transition. Thus, BCR replacement and pre-screening in DZs prevents the accumulation of clones with non-functional receptors and facilitates selection in the LZ.

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.

Sequencing of human genomes extracted from single cancer cells isolated in a valveless microfluidic device.

Sequencing the genomes of individual cells enables the direct determination of genetic heterogeneity amongst cells within a population. We have developed an injection-moulded valveless microfluidic device in which single cells from colorectal cancer derived cell lines (LS174T, LS180 and RKO) and fresh colorectal tumors have been individually trapped, their genomes extracted and prepared for sequencing using multiple displacement amplification (MDA). Ninety nine percent of the DNA sequences obtained mapped to a reference human genome, indicating that there was effectively no contamination of these samples from non-human sources. In addition, most of the reads are correctly paired, with a low percentage of singletons (0.17 ± 0.06%) and we obtain genome coverages approaching 90%. To achieve this high quality, our device design and process shows that amplification can be conducted in microliter volumes as long as the lysis is in sub-nanoliter volumes. Our data thus demonstrates that high quality whole genome sequencing of single cells can be achieved using a relatively simple, inexpensive and scalable device. Detection of genetic heterogeneity at the single cell level, as we have demonstrated for freshly obtained single cancer cells, could soon become available as a clinical tool to precisely match treatment with the properties of a patient's own tumor.

Sasquatch: predicting the impact of regulatory SNPs on transcription factor binding from cell- and tissue-specific DNase footprints.

In the era of genome-wide association studies (GWAS) and personalized medicine, predicting the impact of single nucleotide polymorphisms (SNPs) in regulatory elements is an important goal. Current approaches to determine the potential of regulatory SNPs depend on inadequate knowledge of cell-specific DNA binding motifs. Here, we present Sasquatch, a new computational approach that uses DNase footprint data to estimate and visualize the effects of noncoding variants on transcription factor binding. Sasquatch performs a comprehensive k-mer-based analysis of DNase footprints to determine any k-mer's potential for protein binding in a specific cell type and how this may be changed by sequence variants. Therefore, Sasquatch uses an unbiased approach, independent of known transcription factor binding sites and motifs. Sasquatch only requires a single DNase-seq data set per cell type, from any genotype, and produces consistent predictions from data generated by different experimental procedures and at different sequence depths. Here we demonstrate the effectiveness of Sasquatch using previously validated functional SNPs and benchmark its performance against existing approaches. Sasquatch is available as a versatile webtool incorporating publicly available data, including the human ENCODE collection. Thus, Sasquatch provides a powerful tool and repository for prioritizing likely regulatory SNPs in the noncoding genome.

A biallelic mutation in IL6ST encoding the GP130 co-receptor causes immunodeficiency and craniosynostosis.

Multiple cytokines, including interleukin 6 (IL-6), IL-11, IL-27, oncostatin M (OSM), and leukemia inhibitory factor (LIF), signal via the common GP130 cytokine receptor subunit. In this study, we describe a patient with a homozygous mutation of IL6ST (encoding GP130 p.N404Y) who presented with recurrent infections, eczema, bronchiectasis, high IgE, eosinophilia, defective B cell memory, and an impaired acute-phase response, as well as skeletal abnormalities including craniosynostosis. The p.N404Y missense substitution is associated with loss of IL-6, IL-11, IL-27, and OSM signaling but a largely intact LIF response. This study identifies a novel immunodeficiency with phenotypic similarities to STAT3 hyper-IgE syndrome caused by loss of function of GP130.

Localized TWIST1 and TWIST2 basic domain substitutions cause four distinct human diseases that can be modeled in Caenorhabditis elegans.

Twist transcription factors, members of the basic helix-loop-helix family, play crucial roles in mesoderm development in all animals. Humans have two paralogous genes, TWIST1 and TWIST2, and mutations in each gene have been identified in specific craniofacial disorders. Here, we describe a new clinical entity, Sweeney-Cox syndrome, associated with distinct de novo amino acid substitutions (p.Glu117Val and p.Glu117Gly) at a highly conserved glutamic acid residue located in the basic DNA binding domain of TWIST1, in two subjects with frontonasal dysplasia and additional malformations. Although about one hundred different TWIST1 mutations have been reported in patients with the dominant haploinsufficiency Saethre-Chotzen syndrome (typically associated with craniosynostosis), substitutions uniquely affecting the Glu117 codon were not observed previously. Recently, subjects with Barber-Say and Ablepharon-Macrostomia syndromes were found to harbor heterozygous missense substitutions in the paralogous glutamic acid residue in TWIST2 (p.Glu75Ala, p.Glu75Gln and p.Glu75Lys). To study systematically the effects of these substitutions in individual cells of the developing mesoderm, we engineered all five disease-associated alleles into the equivalent Glu29 residue encoded by hlh-8, the single Twist homolog present in Caenorhabditis elegans. This allelic series revealed that different substitutions exhibit graded severity, in terms of both gene expression and cellular phenotype, which we incorporate into a model explaining the various human disease phenotypes. The genetic analysis favors a predominantly dominant-negative mechanism for the action of amino acid substitutions at this highly conserved glutamic acid residue and illustrates the value of systematic mutagenesis of C. elegans for focused investigation of human disease processes.

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 CDC45, Encoding an Essential Component of the Pre-initiation Complex, Cause Meier-Gorlin Syndrome and Craniosynostosis.

DNA replication precisely duplicates the genome to ensure stable inheritance of genetic information. Impaired licensing of origins of replication during the G1 phase of the cell cycle has been implicated in Meier-Gorlin syndrome (MGS), a disorder defined by the triad of short stature, microtia, and a/hypoplastic patellae. Biallelic partial loss-of-function mutations in multiple components of the pre-replication complex (preRC; ORC1, ORC4, ORC6, CDT1, or CDC6) as well as de novo stabilizing mutations in the licensing inhibitor, GMNN, cause MGS. Here we report the identification of mutations in CDC45 in 15 affected individuals from 12 families with MGS and/or craniosynostosis. CDC45 encodes a component of both the pre-initiation (preIC) and CMG helicase complexes, required for initiation of DNA replication origin firing and ongoing DNA synthesis during S-phase itself, respectively, and hence is functionally distinct from previously identified MGS-associated genes. The phenotypes of affected individuals range from syndromic coronal craniosynostosis to severe growth restriction, fulfilling diagnostic criteria for Meier-Gorlin syndrome. All mutations identified were biallelic and included synonymous mutations altering splicing of physiological CDC45 transcripts, as well as amino acid substitutions expected to result in partial loss of function. Functionally, mutations reduce levels of full-length transcripts and protein in subject cells, consistent with partial loss of CDC45 function and a predicted limited rate of DNA replication and cell proliferation. Our findings therefore implicate the preIC as an additional protein complex involved in the etiology of MGS and connect the core cellular machinery of genome replication with growth, chondrogenesis, and cranial suture homeostasis.

Ubiquitin-mediated fluctuations in MHC class II facilitate efficient germinal center B cell responses.

Antibody affinity maturation occurs in germinal centers (GCs) through iterative rounds of somatic hypermutation and selection. Selection involves B cells competing for T cell help based on the amount of antigen they capture and present on their MHC class II (MHCII) proteins. How GC B cells are able to rapidly and repeatedly transition between mutating their B cell receptor genes and then being selected shortly after is not known. We report that MHCII surface levels and degradation are dynamically regulated in GC B cells. Through ectopic expression of a photoconvertible MHCII-mKikGR chimeric gene, we found that individual GC B cells differed in the rates of MHCII protein turnover. Fluctuations in surface MHCII levels were dependent on ubiquitination and the E3 ligase March1. Increases in March1 expression in centroblasts correlated with decreases in surface MHCII levels, whereas CD83 expression in centrocytes helped to stabilize MHCII at that stage. Defects in MHCII ubiquitination caused GC B cells to accumulate greater amounts of a specific peptide-MHCII (pMHCII), suggesting that MHCII turnover facilitates the replacement of old complexes. We propose that pMHCII complexes are periodically targeted for degradation in centroblasts to favor the presentation of recently acquired antigens, thereby promoting the fidelity and efficiency of selection.

Identification of Intragenic Exon Deletions and Duplication of TCF12 by Whole Genome or Targeted Sequencing as a Cause of TCF12-Related Craniosynostosis.

TCF12-related craniosynostosis can be caused by small heterozygous loss-of-function mutations in TCF12. Large intragenic rearrangements, however, have not been described yet. Here, we present the identification of four large rearrangements in TCF12 causing TCF12-related craniosynostosis. Whole-genome sequencing was applied on the DNA of 18 index cases with coronal synostosis and their family members (43 samples in total). The data were analyzed using an autosomal-dominant disease model. Structural variant analysis reported intragenic exon deletions (of sizes 84.9, 8.6, and 5.4 kb) in TCF12 in three different families. The results were confirmed by deletion-specific PCR and dideoxy-sequence analysis. Separately, targeted sequencing of the TCF12 genomic region in a patient with coronal synostosis identified a tandem duplication of 11.3 kb. The pathogenic effect of this duplication was confirmed by cDNA analysis. These findings indicate the importance of screening for larger rearrangements in patients suspected to have TCF12-related craniosynostosis.

A Recurrent Mosaic Mutation in SMO, Encoding the Hedgehog Signal Transducer Smoothened, Is the Major Cause of Curry-Jones Syndrome.

Curry-Jones syndrome (CJS) is a multisystem disorder characterized by patchy skin lesions, polysyndactyly, diverse cerebral malformations, unicoronal craniosynostosis, iris colobomas, microphthalmia, and intestinal malrotation with myofibromas or hamartomas. Cerebellar medulloblastoma has been described in a single affected individual; in another, biopsy of skin lesions showed features of trichoblastoma. The combination of asymmetric clinical features, patchy skin manifestations, and neoplastic association previously led to the suggestion that this could be a mosaic condition, possibly involving hedgehog (Hh) signaling. Here, we show that CJS is caused by recurrent somatic mosaicism for a nonsynonymous variant in SMO (c.1234C>T [p.Leu412Phe]), encoding smoothened (SMO), a G-protein-coupled receptor that transduces Hh signaling. We identified eight mutation-positive individuals (two of whom had not been reported previously) with highly similar phenotypes and demonstrated varying amounts of the mutant allele in different tissues. We present detailed findings from brain MRI in three mutation-positive individuals. Somatic SMO mutations that result in constitutive activation have been described in several tumors, including medulloblastoma, ameloblastoma, and basal cell carcinoma. Strikingly, the most common of these mutations is the identical nonsynonymous variant encoding p.Leu412Phe. Furthermore, this substitution has been shown to activate SMO in the absence of Hh signaling, providing an explanation for tumor development in CJS. This raises therapeutic possibilities for using recently generated Hh-pathway inhibitors. In summary, our work uncovers the major genetic cause of CJS and illustrates strategies for gene discovery in the context of low-level tissue-specific somatic mosaicism.