Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms.

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

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).

De Novo Variants in CNOT1, a Central Component of the CCR4-NOT Complex Involved in Gene Expression and RNA and Protein Stability, Cause Neurodevelopmental Delay.

CNOT1 is a member of the CCR4-NOT complex, which is a master regulator, orchestrating gene expression, RNA deadenylation, and protein ubiquitination. We report on 39 individuals with heterozygous de novo CNOT1 variants, including missense, splice site, and nonsense variants, who present with a clinical spectrum of intellectual disability, motor delay, speech delay, seizures, hypotonia, and behavioral problems. To link CNOT1 dysfunction to the neurodevelopmental phenotype observed, we generated variant-specific Drosophila models, which showed learning and memory defects upon CNOT1 knockdown. Introduction of human wild-type CNOT1 was able to rescue this phenotype, whereas mutants could not or only partially, supporting our hypothesis that CNOT1 impairment results in neurodevelopmental delay. Furthermore, the genetic interaction with autism-spectrum genes, such as ASH1L, DYRK1A, MED13, and SHANK3, was impaired in our Drosophila models. Molecular characterization of CNOT1 variants revealed normal CNOT1 expression levels, with both mutant and wild-type alleles expressed at similar levels. Analysis of protein-protein interactions with other members indicated that the CCR4-NOT complex remained intact. An integrated omics approach of patient-derived genomics and transcriptomics data suggested only minimal effects on endonucleolytic nonsense-mediated mRNA decay components, suggesting that de novo CNOT1 variants are likely haploinsufficient hypomorph or neomorph, rather than dominant negative. In summary, we provide strong evidence that de novo CNOT1 variants cause neurodevelopmental delay with a wide range of additional co-morbidities. Whereas the underlying pathophysiological mechanism warrants further analysis, our data demonstrate an essential and central role of the CCR4-NOT complex in human brain development.

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.

Missense Mutations in NKAP Cause a Disorder of Transcriptional Regulation Characterized by Marfanoid Habitus and Cognitive Impairment.

NKAP is a ubiquitously expressed nucleoplasmic protein that is currently known as a transcriptional regulatory molecule via its interaction with HDAC3 and spliceosomal proteins. Here, we report a disorder of transcriptional regulation due to missense mutations in the X chromosome gene, NKAP. These mutations are clustered in the C-terminal region of NKAP where NKAP interacts with HDAC3 and post-catalytic spliceosomal complex proteins. Consistent with a role for the C-terminal region of NKAP in embryogenesis, nkap mutant zebrafish with a C-terminally truncated NKAP demonstrate severe developmental defects. The clinical features of affected individuals are highly conserved and include developmental delay, hypotonia, joint contractures, behavioral abnormalities, Marfanoid habitus, and scoliosis. In affected cases, transcriptome analysis revealed the presence of a unique transcriptome signature, which is characterized by the downregulation of long genes with higher exon numbers. These observations indicate the critical role of NKAP in transcriptional regulation and demonstrate that perturbations of the C-terminal region lead to developmental defects in both humans and zebrafish.

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.

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.

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.

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.

GATAD2B-associated neurodevelopmental disorder (GAND): clinical and molecular insights into a NuRD-related disorder.

PURPOSE: Determination of genotypic/phenotypic features of GATAD2B-associated neurodevelopmental disorder (GAND). METHODS: Fifty GAND subjects were evaluated to determine consistent genotypic/phenotypic features. Immunoprecipitation assays utilizing in vitro transcription-translation products were used to evaluate GATAD2B missense variants' ability to interact with binding partners within the nucleosome remodeling and deacetylase (NuRD) complex. RESULTS: Subjects had clinical findings that included macrocephaly, hypotonia, intellectual disability, neonatal feeding issues, polyhydramnios, apraxia of speech, epilepsy, and bicuspid aortic valves. Forty-one novelGATAD2B variants were identified with multiple variant types (nonsense, truncating frameshift, splice-site variants, deletions, and missense). Seven subjects were identified with missense variants that localized within two conserved region domains (CR1 or CR2) of the GATAD2B protein. Immunoprecipitation assays revealed several of these missense variants disrupted GATAD2B interactions with its NuRD complex binding partners. CONCLUSIONS: A consistent GAND phenotype was caused by a range of genetic variants in GATAD2B that include loss-of-function and missense subtypes. Missense variants were present in conserved region domains that disrupted assembly of NuRD complex proteins. GAND's clinical phenotype had substantial clinical overlap with other disorders associated with the NuRD complex that involve CHD3 and CHD4, with clinical features of hypotonia, intellectual disability, cardiac defects, childhood apraxia of speech, and macrocephaly.

Correction: GATAD2B-associated neurodevelopmental disorder (GAND): clinical and molecular insights into a NuRD-related disorder.

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

BCL11A Haploinsufficiency Causes an Intellectual Disability Syndrome and Dysregulates Transcription.

Intellectual disability (ID) is a common condition with considerable genetic heterogeneity. Next-generation sequencing of large cohorts has identified an increasing number of genes implicated in ID, but their roles in neurodevelopment remain largely unexplored. Here we report an ID syndrome caused by de novo heterozygous missense, nonsense, and frameshift mutations in BCL11A, encoding a transcription factor that is a putative member of the BAF swi/snf chromatin-remodeling complex. Using a comprehensive integrated approach to ID disease modeling, involving human cellular analyses coupled to mouse behavioral, neuroanatomical, and molecular phenotyping, we provide multiple lines of functional evidence for phenotypic effects. The etiological missense variants cluster in the amino-terminal region of human BCL11A, and we demonstrate that they all disrupt its localization, dimerization, and transcriptional regulatory activity, consistent with a loss of function. We show that Bcl11a haploinsufficiency in mice causes impaired cognition, abnormal social behavior, and microcephaly in accordance with the human phenotype. Furthermore, we identify shared aberrant transcriptional profiles in the cortex and hippocampus of these mouse models. Thus, our work implicates BCL11A haploinsufficiency in neurodevelopmental disorders and defines additional targets regulated by this gene, with broad relevance for our understanding of ID and related syndromes.

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 cohort of 17 patients with kyphoscoliotic Ehlers-Danlos syndrome caused by biallelic mutations in FKBP14: expansion of the clinical and mutational spectrum and description of the natural history.

PurposeIn 2012 we reported in six individuals a clinical condition almost indistinguishable from PLOD1-kyphoscoliotic Ehlers-Danlos syndrome (PLOD1-kEDS), caused by biallelic mutations in FKBP14, and characterized by progressive kyphoscoliosis, myopathy, and hearing loss in addition to connective tissue abnormalities such as joint hypermobility and hyperelastic skin. FKBP14 is an ER-resident protein belonging to the family of FK506-binding peptidyl-prolyl cis-trans isomerases (PPIases); it catalyzes the folding of type III collagen and interacts with type III, type VI, and type X collagens. Only nine affected individuals have been reported to date.MethodsWe report on a cohort of 17 individuals with FKBP14-kEDS and the follow-up of three previously reported patients, and provide an extensive overview of the disorder and its natural history based on clinical, biochemical, and molecular genetics data.ResultsBased on the frequency of the clinical features of 23 patients from the present and previous cohorts, we define major and minor features of FKBP14-kEDS. We show that myopathy is confirmed by histology and muscle imaging only in some patients, and that hearing impairment is predominantly sensorineural and may not be present in all individuals.ConclusionOur data further support the extensive clinical overlap with PLOD1-kEDS and show that vascular complications are rare manifestations of FKBP14-kEDS.

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.

Mandibulofacial Dysostosis with Microcephaly: Mutation and Database Update.

Mandibulofacial dysostosis with microcephaly (MFDM) is a multiple malformation syndrome comprising microcephaly, craniofacial anomalies, hearing loss, dysmorphic features, and, in some cases, esophageal atresia. Haploinsufficiency of a spliceosomal GTPase, U5-116 kDa/EFTUD2, is responsible. Here, we review the molecular basis of MFDM in the 69 individuals described to date, and report mutations in 38 new individuals, bringing the total number of reported individuals to 107 individuals from 94 kindreds. Pathogenic EFTUD2 variants comprise 76 distinct mutations and seven microdeletions. Among point mutations, missense substitutions are infrequent (14 out of 76; 18%) relative to stop-gain (29 out of 76; 38%), and splicing (33 out of 76; 43%) mutations. Where known, mutation origin was de novo in 48 out of 64 individuals (75%), dominantly inherited in 12 out of 64 (19%), and due to proven germline mosaicism in four out of 64 (6%). Highly penetrant clinical features include, microcephaly, first and second arch craniofacial malformations, and hearing loss; esophageal atresia is present in an estimated ∼27%. Microcephaly is virtually universal in childhood, with some adults exhibiting late "catch-up" growth and normocephaly at maturity. Occasionally reported anomalies, include vestibular and ossicular malformations, reduced mouth opening, atrophy of cerebral white matter, structural brain malformations, and epibulbar dermoid. All reported EFTUD2 mutations can be found in the EFTUD2 mutation database (http://databases.lovd.nl/shared/genes/EFTUD2).

Biallelic mutations in CYP26B1: A differential diagnosis for Pfeiffer and Antley-Bixler syndromes.

Recently, a newly identified autosomal recessive skeletal dysplasia was described characterized by calvarial abnormalities (including cranium bifidum, coronal, and lambdoid synostosis), oligodactyly, femoral bowing, narrow thorax, small pelvic bones, and radiohumeral synostosis. In the two families described, a more severe phenotype led to in utero lethality in three siblings while in a single patient in a second family the phenotype was sufficiently mild to allow survival to 5 months of age. The disorder is caused by biallelic missense mutations in CYP26B1, which encodes for a cytochrome P450 enzyme responsible for the catabolism of retinoic acid in a temporally and spatially restricted fashion during embryonic development. Here, we provide the third family affected by the disorder and the first affected individual to survive beyond infancy. This woman homozygous for c.1303G>A; p.(Gly435Ser) in CYP26B1, which was associated with multisutural synostosis, radiohumeral synostosis, normal bone mineral density, and apparent intellectual disability, a phenotype with significant similarities to Antley-Bixler and Pfeiffer syndromes. © 2016 Wiley Periodicals, Inc.

PLA2G6, encoding a phospholipase A2, is mutated in neurodegenerative disorders with high brain iron.

Neurodegenerative disorders with high brain iron include Parkinson disease, Alzheimer disease and several childhood genetic disorders categorized as neuroaxonal dystrophies. We mapped a locus for infantile neuroaxonal dystrophy (INAD) and neurodegeneration with brain iron accumulation (NBIA) to chromosome 22q12-q13 and identified mutations in PLA2G6, encoding a calcium-independent group VI phospholipase A2, in NBIA, INAD and the related Karak syndrome. This discovery implicates phospholipases in the pathogenesis of neurodegenerative disorders with iron dyshomeostasis.

Personalized recurrence risk assessment following the birth of a child with a pathogenic de novo mutation.

Following the diagnosis of a paediatric disorder caused by an apparently de novo mutation, a recurrence risk of 1-2% is frequently quoted due to the possibility of parental germline mosaicism; but for any specific couple, this figure is usually incorrect. We present a systematic approach to providing individualized recurrence risk. By combining locus-specific sequencing of multiple tissues to detect occult mosaicism with long-read sequencing to determine the parent-of-origin of the mutation, we show that we can stratify the majority of couples into one of seven discrete categories associated with substantially different risks to future offspring. Among 58 families with a single affected offspring (representing 59 de novo mutations in 49 genes), the recurrence risk for 35 (59%) was decreased below 0.1%, but increased owing to parental mixed mosaicism for 5 (9%)-that could be quantified in semen for paternal cases (recurrence risks of 5.6-12.1%). Implementation of this strategy offers the prospect of driving a major transformation in the practice of genetic counselling.