Single-cell analyses reveal transient retinal progenitor cells in the ciliary margin of developing human retina

The emergence of retinal progenitor cells and differentiation to various retinal cell types represent fundamental processes during retinal development. Herein, we provide a comprehensive single cell characterisation of transcriptional and chromatin accessibility changes that underline retinal progenitor cell specification and differentiation over the course of human retinal development up to midgestation. Our lineage trajectory data demonstrate the presence of early retinal progenitors, which transit to late, and further to transient neurogenic progenitors, that give rise to all the retinal neurons. Combining single cell RNA-Seq with spatial transcriptomics of early eye samples, we demonstrate the transient presence of early retinal progenitors in the ciliary margin zone with decreasing occurrence from 8 post-conception week of human development. In retinal progenitor cells, we identified a significant enrichment for transcriptional enhanced associate domain transcription factor binding motifs, which when inhibited led to loss of cycling progenitors and retinal identity in pluripotent stem cell derived organoids.

Heterozygous loss-of-function SMC3 variants are associated with variable growth and developmental features

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia

BACKGROUND

Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus.

METHODS

Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region.

RESULTS

Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements.

CONCLUSION

Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

A human embryonic limb cell atlas resolved in space and time

Human limbs emerge during the fourth post-conception week as mesenchymal buds, which develop into fully formed limbs over the subsequent months1. This process is orchestrated by numerous temporally and spatially restricted gene expression programmes, making congenital alterations in phenotype common2. Decades of work with model organisms have defined the fundamental mechanisms underlying vertebrate limb development, but an in-depth characterization of this process in humans has yet to be performed. Here we detail human embryonic limb development across space and time using single-cell and spatial transcriptomics. We demonstrate extensive diversification of cells from a few multipotent progenitors to myriad differentiated cell states, including several novel cell populations. We uncover two waves of human muscle development, each characterized by different cell states regulated by separate gene expression programmes, and identify musculin (MSC) as a key transcriptional repressor maintaining muscle stem cell identity. Through assembly of multiple anatomically continuous spatial transcriptomic samples using VisiumStitcher, we map cells across a sagittal section of a whole fetal hindlimb. We reveal a clear anatomical segregation between genes linked to brachydactyly and polysyndactyly, and uncover transcriptionally and spatially distinct populations of the mesenchyme in the autopod. Finally, we perform single-cell RNA sequencing on mouse embryonic limbs to facilitate cross-species developmental comparison, finding substantial homology between the two species.

Heterozygous loss-of-function SMC3 variants are associated with variable and incompletely penetrant growth and developmental features

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 13 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated a milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, some instead having intriguing symptomatologies with rational biological links to SMC3 including bone marrow failure, acute myeloid leukemia, and Coats retinal vasculopathy. Analyses of transcriptomic and epigenetic data suggest that SMC3 pLoF variants reduce SMC3 expression but do not result in a blood DNA methylation signature clustering with that of CdLS, and that the global transcriptional signature of SMC3 loss is model-dependent. Our finding of substantial population-scale LoF intolerance in concert with variable penetrance in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multi-layered genomic data paired with careful phenotyping.

Characterization of an eye field-like state during optic vesicle organoid development

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Genomic Diagnosis of Rare Pediatric Disease in the United Kingdom and Ireland

BACKGROUND

Pediatric disorders include a range of highly penetrant, genetically heterogeneous conditions amenable to genomewide diagnostic approaches. Finding a molecular diagnosis is challenging but can have profound lifelong benefits.

METHODS

We conducted a large-scale sequencing study involving more than 13,500 families with probands with severe, probably monogenic, difficult-to-diagnose developmental disorders from 24 regional genetics services in the United Kingdom and Ireland. Standardized phenotypic data were collected, and exome sequencing and microarray analyses were performed to investigate novel genetic causes. We developed an iterative variant analysis pipeline and reported candidate variants to clinical teams for validation and diagnostic interpretation to inform communication with families. Multiple regression analyses were performed to evaluate factors affecting the probability of diagnosis.

RESULTS

A total of 13,449 probands were included in the analyses. On average, we reported 1.0 candidate variant per parent-offspring trio and 2.5 variants per singleton proband. Using clinical and computational approaches to variant classification, we made a diagnosis in approximately 41% of probands (5502 of 13,449). Of 3599 probands in trios who received a diagnosis by clinical assertion, approximately 76% had a pathogenic de novo variant. Another 22% of probands (2997 of 13,449) had variants of uncertain significance in genes that were strongly linked to monogenic developmental disorders. Recruitment in a parent-offspring trio had the largest effect on the probability of diagnosis (odds ratio, 4.70; 95% confidence interval [CI], 4.16 to 5.31). Probands were less likely to receive a diagnosis if they were born extremely prematurely (i.e., 22 to 27 weeks' gestation; odds ratio, 0.39; 95% CI, 0.22 to 0.68), had in utero exposure to antiepileptic medications (odds ratio, 0.44; 95% CI, 0.29 to 0.67), had mothers with diabetes (odds ratio, 0.52; 95% CI, 0.41 to 0.67), or were of African ancestry (odds ratio, 0.51; 95% CI, 0.31 to 0.78).

CONCLUSIONS

Among probands with severe, probably monogenic, difficult-to-diagnose developmental disorders, multimodal analysis of genomewide data had good diagnostic power, even after previous attempts at diagnosis. (Funded by the Health Innovation Challenge Fund and Wellcome Sanger Institute.).

IMPROVE-DD: Integrating multiple phenotype resources optimizes variant evaluation in genetically determined developmental disorders

Diagnosing rare developmental disorders using genome-wide sequencing data commonly necessitates review of multiple plausible candidate variants, often using ontologies of categorical clinical terms. We show that Integrating Multiple Phenotype Resources Optimizes Variant Evaluation in Developmental Disorders (IMPROVE-DD) by incorporating additional classes of data commonly available to clinicians and recorded in health records. In doing so, we quantify the distinct contributions of sex, growth, and development in addition to Human Phenotype Ontology (HPO) terms and demonstrate added value from these readily available information sources. We use likelihood ratios for nominal and quantitative data and propose a classifier for HPO terms in this framework. This Bayesian framework results in more robust diagnoses. Using data systematically collected in the Deciphering Developmental Disorders study, we considered 77 genes with pathogenic/likely pathogenic variants in ≥10 individuals. All genes showed at least a satisfactory prediction by receiver operating characteristic when testing on training data (AUC ≥ 0.6), and HPO terms were the best predictor for the majority of genes, though a minority (13/77) of genes were better predicted by other phenotypic data types. Overall, classifiers based upon multiple integrated phenotypic data sources performed better than those based upon any individual source, and importantly, integrated models produced notably fewer false positives. Finally, we show that IMPROVE-DD models with good predictive performance on cross-validation can be constructed from relatively few individuals. This suggests new strategies for candidate gene prioritization and highlights the value of systematic clinical data collection to support diagnostic programs.

Monoallelic variants resulting in substitutions of MAB21L1 Arg51 Cause Aniridia and microphthalmia

Classical aniridia is a congenital and progressive panocular disorder almost exclusively caused by heterozygous loss-of-function variants at the PAX6 locus. We report nine individuals from five families with severe aniridia and/or microphthalmia (with no detectable PAX6 mutation) with ultrarare monoallelic missense variants altering the Arg51 codon of MAB21L1. These mutations occurred de novo in 3/5 families, with the remaining families being compatible with autosomal dominant inheritance. Mice engineered to carry the p.Arg51Leu change showed a highly-penetrant optic disc anomaly in heterozygous animals with severe microphthalmia in homozygotes. Substitutions of the same codon (Arg51) in MAB21L2, a close homolog of MAB21L1, cause severe ocular and skeletal malformations in humans and mice. The predicted nucleotidyltransferase function of MAB21L1 could not be demonstrated using purified protein with a variety of nucleotide substrates and oligonucleotide activators. Induced expression of GFP-tagged wildtype and mutant MAB21L1 in human cells caused only modest transcriptional changes. Mass spectrometry of immunoprecipitated protein revealed that both mutant and wildtype MAB21L1 associate with transcription factors that are known regulators of PAX6 (MEIS1, MEIS2 and PBX1) and with poly(A) RNA binding proteins. Arg51 substitutions reduce the association of wild-type MAB21L1 with TBL1XR1, a component of the NCoR complex. We found limited evidence for mutation-specific interactions with MSI2/Musashi-2, an RNA-binding proteins with effects on many different developmental pathways. Given that biallelic loss-of-function variants in MAB21L1 result in a milder eye phenotype we suggest that Arg51-altering monoallelic variants most plausibly perturb eye development via a gain-of-function mechanism.

Robust Genetic Analysis of the X-Linked Anophthalmic (Ie) Mouse

Anophthalmia (missing eye) describes a failure of early embryonic ocular development. Mutations in a relatively small set of genes account for 75% of bilateral anophthalmia cases, yet 25% of families currently are left without a molecular diagnosis. Here, we report our experimental work that aimed to uncover the developmental and genetic basis of the anophthalmia characterising the X-linked Ie (eye-ear reduction) X-ray-induced allele in mouse that was first identified in 1947. Histological analysis of the embryonic phenotype showed failure of normal eye development after the optic vesicle stage with particularly severe malformation of the ventral retina. Linkage analysis mapped this mutation to a ~6 Mb region on the X chromosome. Short- and long-read whole-genome sequencing (WGS) of affected and unaffected male littermates confirmed the Ie linkage but identified no plausible causative variants or structural rearrangements. These analyses did reduce the critical candidate interval and revealed evidence of multiple variants within the ancestral DNA, although none were found that altered coding sequences or that were unique to Ie. To investigate early embryonic events at a genetic level, we then generated mouse ES cells derived from male Ie embryos and wild type littermates. RNA-seq and accessible chromatin sequencing (ATAC-seq) data generated from cultured optic vesicle organoids did not reveal any large differences in gene expression or accessibility of putative cis-regulatory elements between Ie and wild type. However, an unbiased TF-footprinting analysis of accessible chromatin regions did provide evidence of a genome-wide reduction in binding of transcription factors associated with ventral eye development in Ie, and evidence of an increase in binding of the Zic-family of transcription factors, including Zic3, which is located within the Ie-refined critical interval. We conclude that the refined Ie critical region at chrX: 56,145,000-58,385,000 contains multiple genetic variants that may be linked to altered cis regulation but does not contain a convincing causative mutation. Changes in the binding of key transcription factors to chromatin causing altered gene expression during development, possibly through a subtle mis-regulation of Zic3, presents a plausible cause for the anophthalmia phenotype observed in Ie, but further work is required to determine the precise causative allele and its genetic mechanism.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome

BACKGROUND

The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts.

METHODS

We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups.

RESULTS

We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants.

CONCLUSIONS

These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

Creation and evaluation of full-text literature-derived, feature-weighted disease models of genetically determined developmental disorders

There are >2500 different genetically determined developmental disorders (DD), which, as a group, show very high levels of both locus and allelic heterogeneity. This has led to the wide-spread use of evidence-based filtering of genome-wide sequence data as a diagnostic tool in DD. Determining whether the association of a filtered variant at a specific locus is a plausible explanation of the phenotype in the proband is crucial and commonly requires extensive manual literature review by both clinical scientists and clinicians. Access to a database of weighted clinical features extracted from rigorously curated literature would increase the efficiency of this process and facilitate the development of robust phenotypic similarity metrics. However, given the large and rapidly increasing volume of published information, conventional biocuration approaches are becoming impractical. Here, we present a scalable, automated method for the extraction of categorical phenotypic descriptors from the full-text literature. Papers identified through literature review were downloaded and parsed using the Cadmus custom retrieval package. Human Phenotype Ontology terms were extracted using MetaMap, with 76-84% precision and 65-73% recall. Mean terms per paper increased from 9 in title + abstract, to 68 using full text. We demonstrate that these literature-derived disease models plausibly reflect true disease expressivity more accurately than widely used manually curated models, through comparison with prospectively gathered data from the Deciphering Developmental Disorders study. The area under the curve for receiver operating characteristic (ROC) curves increased by 5-10% through the use of literature-derived models. This work shows that scalable automated literature curation increases performance and adds weight to the need for this strategy to be integrated into informatic variant analysis pipelines. Database URL: https://doi.org/10.1093/database/baac038.

Detecting cryptic clinically relevant structural variation in exome-sequencing data increases diagnostic yield for developmental disorders

Structural variation (SV) describes a broad class of genetic variation greater than 50 bp in size. SVs can cause a wide range of genetic diseases and are prevalent in rare developmental disorders (DDs). Individuals presenting with DDs are often referred for diagnostic testing with chromosomal microarrays (CMAs) to identify large copy-number variants (CNVs) and/or with single-gene, gene-panel, or exome sequencing (ES) to identify single-nucleotide variants, small insertions/deletions, and CNVs. However, individuals with pathogenic SVs undetectable by conventional analysis often remain undiagnosed. Consequently, we have developed the tool InDelible, which interrogates short-read sequencing data for split-read clusters characteristic of SV breakpoints. We applied InDelible to 13,438 probands with severe DDs recruited as part of the Deciphering Developmental Disorders (DDD) study and discovered 63 rare, damaging variants in genes previously associated with DDs missed by standard SNV, indel, or CNV discovery approaches. Clinical review of these 63 variants determined that about half (30/63) were plausibly pathogenic. InDelible was particularly effective at ascertaining variants between 21 and 500 bp in size and increased the total number of potentially pathogenic variants identified by DDD in this size range by 42.9%. Of particular interest were seven confirmed de novo variants in MECP2, which represent 35.0% of all de novo protein-truncating variants in MECP2 among DDD study participants. InDelible provides a framework for the discovery of pathogenic SVs that are most likely missed by standard analytical workflows and has the potential to improve the diagnostic yield of ES across a broad range of genetic diseases.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Identification and functional modelling of plausibly causative cis-regulatory variants in a highly-selected cohort with X-linked intellectual disability

Identifying causative variants in cis-regulatory elements (CRE) in neurodevelopmental disorders has proven challenging. We have used in vivo functional analyses to categorize rigorously filtered CRE variants in a clinical cohort that is plausibly enriched for causative CRE mutations: 48 unrelated males with a family history consistent with X-linked intellectual disability (XLID) in whom no detectable cause could be identified in the coding regions of the X chromosome (chrX). Targeted sequencing of all chrX CRE identified six rare variants in five affected individuals that altered conserved bases in CRE targeting known XLID genes and segregated appropriately in families. Two of these variants, FMR1CRE and TENM1CRE, showed consistent site- and stage-specific differences of enhancer function in the developing zebrafish brain using dual-color fluorescent reporter assay. Mouse models were created for both variants. In male mice Fmr1CRE induced alterations in neurodevelopmental Fmr1 expression, olfactory behavior and neurophysiological indicators of FMRP function. The absence of another likely causative variant on whole genome sequencing further supported FMR1CRE as the likely basis of the XLID in this family. Tenm1CRE mice showed no phenotypic anomalies. Following the release of gnomAD 2.1, reanalysis showed that TENM1CRE exceeded the maximum plausible population frequency of a XLID causative allele. Assigning causative status to any ultra-rare CRE variant remains problematic and requires disease-relevant in vivo functional data from multiple sources. The sequential and bespoke nature of such analyses renders them time-consuming and challenging to scale for routine clinical use.

Evaluating variants classified as pathogenic in ClinVar in the DDD Study

PURPOSE

Automated variant filtering is an essential part of diagnostic genome-wide sequencing but may generate false negative results. We sought to investigate whether some previously identified pathogenic variants may be being routinely excluded by standard variant filtering pipelines.

METHODS

We evaluated variants that were previously classified as pathogenic or likely pathogenic in ClinVar in known developmental disorder genes using exome sequence data from the Deciphering Developmental Disorders (DDD) study.

RESULTS

Of these ClinVar pathogenic variants, 3.6% were identified among 13,462 DDD probands, and 1134/1352 (83.9%) had already been independently communicated to clinicians using DDD variant filtering pipelines as plausibly pathogenic. The remaining 218 variants failed consequence, inheritance, or other automated variant filters. Following clinical review of these additional variants, we were able to identify 112 variants in 107 (0.8%) DDD probands as potential diagnoses.

CONCLUSION

Lower minor allele frequency (<0.0005%) and higher gold star review status in ClinVar (>1 star) are good predictors of a previously identified variant being plausibly diagnostic for developmental disorders. However, around half of previously identified pathogenic variants excluded by automated variant filtering did not appear to be disease-causing, underlining the continued need for clinical evaluation of candidate variants as part of the diagnostic process.

The contribution of X-linked coding variation to severe developmental disorders

Over 130 X-linked genes have been robustly associated with developmental disorders, and X-linked causes have been hypothesised to underlie the higher developmental disorder rates in males. Here, we evaluate the burden of X-linked coding variation in 11,044 developmental disorder patients, and find a similar rate of X-linked causes in males and females (6.0% and 6.9%, respectively), indicating that such variants do not account for the 1.4-fold male bias. We develop an improved strategy to detect X-linked developmental disorders and identify 23 significant genes, all of which were previously known, consistent with our inference that the vast majority of the X-linked burden is in known developmental disorder-associated genes. Importantly, we estimate that, in male probands, only 13% of inherited rare missense variants in known developmental disorder-associated genes are likely to be pathogenic. Our results demonstrate that statistical analysis of large datasets can refine our understanding of modes of inheritance for individual X-linked disorders.

Evidence for 28 genetic disorders discovered by combining healthcare and research data

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

Genomically Aided Diagnosis of Severe Developmental Disorders

Our ability to make accurate and specific genetic diagnoses in individuals with severe developmental disorders has been transformed by data derived from genomic sequencing technologies. These data reveal both the patterns and rates of different mutational mechanisms and identify regions of the human genome with fewer mutations than would be expected. In outbred populations, the most common identifiable cause of severe developmental disorders is de novo mutation affecting the coding region in one of approximately 500 different genes, almost universally showing constraint. Simply combining the location of a de novo genomic event with its predicted consequence on the gene product gives significant diagnostic power. Our knowledge of the diversity of phenotypic consequences associated with comparable diagnostic genotypes at each locus is improving. Computationally useful phenotype data will improve diagnostic interpretation of ultrarare genetic variants and, in the long run, indicate which specific embryonic processes have been perturbed.

Recurrent heterozygous PAX6 missense variants cause severe bilateral microphthalmia via predictable effects on DNA-protein interaction

PURPOSE

Most classical aniridia is caused by PAX6 haploinsufficiency. PAX6 missense variants can be hypomorphic or mimic haploinsufficiency. We hypothesized that missense variants also cause previously undescribed disease by altering the affinity and/or specificity of PAX6 genomic interactions.

METHODS

We screened PAX6 in 372 individuals with bilateral microphthalmia, anophthalmia, or coloboma (MAC) from the Medical Research Council Human Genetics Unit eye malformation cohort (HGUeye) and reviewed data from the Deciphering Developmental Disorders study. We performed cluster analysis on PAX6-associated ocular phenotypes by variant type and molecular modeling of the structural impact of 86 different PAX6 causative missense variants.

RESULTS

Eight different PAX6 missense variants were identified in 17 individuals (15 families) with MAC, accounting for 4% (15/372) of our cohort. Seven altered the paired domain (p.[Arg26Gln]x1, p.[Gly36Val]x1, p.[Arg38Trp]x2, p.[Arg38Gln]x1, p.[Gly51Arg]x2, p.[Ser54Arg]x2, p.[Asn124Lys]x5) and one the homeodomain (p.[Asn260Tyr]x1). p.Ser54Arg and p.Asn124Lys were exclusively associated with severe bilateral microphthalmia. MAC-associated variants were predicted to alter but not ablate DNA interaction, consistent with the electrophoretic mobility shifts observed using mutant paired domains with well-characterized PAX6-binding sites. We found no strong evidence for novel PAX6-associated extraocular disease.

CONCLUSION

Altering the affinity and specificity of PAX6-binding genome-wide provides a plausible mechanism for the worse-than-null effects of MAC-associated missense variants.

Microphthalmia, Linear Skin Defects, Callosal Agenesis, and Cleft Palate in a Patient with Deletion at Xp22.3p22.2

The authors describe the clinical findings observed in a Brazilian girl that are suggestive of microphthalmia and linear skin defects (MLS) also known as MIDAS syndrome (OMIM #309801). She also presented with short stature, agenesis of corpus callosum, cleft palate, enamel defects, and genitourinary anomalies, which are rarely reported within the clinical spectrum of MLS. The 11,5 Mb deletion in Xp22.3p22.2 observed in the patient includes the entire HCCS gene (responsible for the MLS phenotype) and also encompasses several other genes involved with behavioral phenotypes, craniofacial and central nervous system development such as MID1, NLGN4X, AMELX , ARHGAP6, and TBL1X. The whole clinical features of our proband possibly represents an unusual MLS syndromic phenotype caused by an Xp22.3p22.2 continuous gene deletion.

Contribution of retrotransposition to developmental disorders

Mobile genetic Elements (MEs) are segments of DNA which can copy themselves and other transcribed sequences through the process of retrotransposition (RT). In humans several disorders have been attributed to RT, but the role of RT in severe developmental disorders (DD) has not yet been explored. Here we identify RT-derived events in 9738 exome sequenced trios with DD-affected probands. We ascertain 9 de novo MEs, 4 of which are likely causative of the patient's symptoms (0.04%), as well as 2 de novo gene retroduplications. Beyond identifying likely diagnostic RT events, we estimate genome-wide germline ME mutation rate and selective constraint and demonstrate that coding RT events have signatures of purifying selection equivalent to those of truncating mutations. Overall, our analysis represents a comprehensive interrogation of the impact of retrotransposition on protein coding genes and a framework for future evolutionary and disease studies.

Enabling Global Clinical Collaborations on Identifiable Patient Data: The Minerva Initiative

The clinical utility of computational phenotyping for both genetic and rare diseases is increasingly appreciated; however, its true potential is yet to be fully realized. Alongside the growing clinical and research availability of sequencing technologies, precise deep and scalable phenotyping is required to serve unmet need in genetic and rare diseases. To improve the lives of individuals affected with rare diseases through deep phenotyping, global big data interrogation is necessary to aid our understanding of disease biology, assist diagnosis, and develop targeted treatment strategies. This includes the application of cutting-edge machine learning methods to image data. As with most digital tools employed in health care, there are ethical and data governance challenges associated with using identifiable personal image data. There are also risks with failing to deliver on the patient benefits of these new technologies, the biggest of which is posed by data siloing. The Minerva Initiative has been designed to enable the public good of deep phenotyping while mitigating these ethical risks. Its open structure, enabling collaboration and data sharing between individuals, clinicians, researchers and private enterprise, is key for delivering precision public health.

Cardiac Genetic Predisposition in Sudden Infant Death Syndrome

BACKGROUND

Sudden infant death syndrome (SIDS) is a leading cause of postneonatal mortality. Genetic heart diseases (GHDs) underlie some cases of SIDS.

OBJECTIVES

This study aimed to determine the spectrum and prevalence of GHD-associated mutations as a potential monogenic basis for SIDS.

METHODS

A cohort of 419 unrelated SIDS cases (257 male; average age 2.7 ± 1.9 months) underwent whole exome sequencing and a targeted analysis of 90 GHD-susceptibility genes. The yield of "potentially informative," ultra-rare variants (minor allele frequency <0.00005) in GHD-associated genes was assessed.

RESULTS

Overall, 53 of 419 (12.6%) SIDS cases had ≥1 "potentially informative," GHD-associated variant. The yield was 14.9% (21 of 141) for mixed-European ancestry cases and 11.5% (32 of 278) for European ancestry SIDS cases. Infants older than 4 months were more likely to host a "potentially informative" GHD-associated variant. There was significant overrepresentation of ultra-rare nonsynonymous variants in European SIDS cases (18 of 278 [6.5%]) versus European control subjects (30 of 973 [3.1%]; p = 0.013) when combining all 4 major cardiac channelopathy genes (KCNQ1, KCNH2, SCN5A, and RYR2). According to the American College of Medical Genetics guidelines, only 18 of 419 (4.3%) SIDS cases hosted a "pathogenic" or "likely pathogenic" variant.

CONCLUSIONS

Less than 15% of more than 400 SIDS cases had a "potentially informative" variant in a GHD-susceptibility gene, predominantly in the 4- to 12-month age group. Only 4.3% of cases possessed immediately clinically actionable variants. Consistent with previous studies, ultra-rare, nonsynonymous variants within the major cardiac channelopathy-associated genes were overrepresented in SIDS cases in infants of European ethnicity. These findings have major implications for the investigation of SIDS cases and families.

Flexible and scalable diagnostic filtering of genomic variants using G2P with Ensembl VEP

We aimed to develop an efficient, flexible and scalable approach to diagnostic genome-wide sequence analysis of genetically heterogeneous clinical presentations. Here we present G2P ( www.ebi.ac.uk/gene2phenotype ) as an online system to establish, curate and distribute datasets for diagnostic variant filtering via association of allelic requirement and mutational consequence at a defined locus with phenotypic terms, confidence level and evidence links. An extension to Ensembl Variant Effect Predictor (VEP), VEP-G2P was used to filter both disease-associated and control whole exome sequence (WES) with Developmental Disorders G2P (G2PDD; 2044 entries). VEP-G2PDD shows a sensitivity/precision of 97.3%/33% for de novo and 81.6%/22.7% for inherited pathogenic genotypes respectively. Many of the missing genotypes are likely false-positive pathogenic assignments. The expected number and discriminative features of background genotypes are defined using control WES. Using only human genetic data VEP-G2P performs well compared to other freely-available diagnostic systems and future phenotypic matching capabilities should further enhance performance.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Ca_v2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca²⁺ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Ca_v2.2 in normal human neurodevelopment.

NAA10 polyadenylation signal variants cause syndromic microphthalmia

BACKGROUND

A single variant in NAA10 (c.471+2T>A), the gene encoding N-acetyltransferase 10, has been associated with Lenz microphthalmia syndrome. In this study, we aimed to identify causative variants in families with syndromic X-linked microphthalmia.

METHODS

Three families, including 15 affected individuals with syndromic X-linked microphthalmia, underwent analyses including linkage analysis, exome sequencing and targeted gene sequencing. The consequences of two identified variants in NAA10 were evaluated using quantitative PCR and RNAseq.

RESULTS

Genetic linkage analysis in family 1 supported a candidate region on Xq27-q28, which included NAA10. Exome sequencing identified a hemizygous NAA10 polyadenylation signal (PAS) variant, chrX:153,195,397T>C, c.*43A>G, which segregated with the disease. Targeted sequencing of affected males from families 2 and 3 identified distinct NAA10 PAS variants, chrX:g.153,195,401T>C, c.*39A>G and chrX:g.153,195,400T>C, c.*40A>G. All three variants were absent from gnomAD. Quantitative PCR and RNAseq showed reduced NAA10 mRNA levels and abnormal 3' UTRs in affected individuals. Targeted sequencing of NAA10 in 376 additional affected individuals failed to identify variants in the PAS.

CONCLUSION

These data show that PAS variants are the most common variant type in NAA10-associated syndromic microphthalmia, suggesting reduced RNA is the molecular mechanism by which these alterations cause microphthalmia/anophthalmia. We reviewed recognised variants in PAS associated with Mendelian disorders and identified only 23 others, indicating that NAA10 harbours more than 10% of all known PAS variants. We hypothesise that PAS in other genes harbour unrecognised pathogenic variants associated with Mendelian disorders. The systematic interrogation of PAS could improve genetic testing yields.

ITPase deficiency causes a Martsolf-like syndrome with a lethal infantile dilated cardiomyopathy

Typical Martsolf syndrome is characterized by congenital cataracts, postnatal microcephaly, developmental delay, hypotonia, short stature and biallelic hypomorphic mutations in either RAB3GAP1 or RAB3GAP2. Genetic analysis of 85 unrelated "mutation negative" probands with Martsolf or Martsolf-like syndromes identified two individuals with different homozygous null mutations in ITPA, the gene encoding inosine triphosphate pyrophosphatase (ITPase). Both probands were from multiplex families with a consistent, lethal and highly distinctive disorder; a Martsolf-like syndrome with infantile-onset dilated cardiomyopathy. Severe ITPase-deficiency has been previously reported with infantile epileptic encephalopathy (MIM 616647). ITPase acts to prevent incorporation of inosine bases (rI/dI) into RNA and DNA. In Itpa-null cells dI was undetectable in genomic DNA. dI could be identified at a low level in mtDNA without detectable mitochondrial genome instability, mtDNA depletion or biochemical dysfunction of the mitochondria. rI accumulation was detectable in proband-derived lymphoblastoid RNA. In Itpa-null mouse embryos rI was detectable in the brain and kidney with the highest level seen in the embryonic heart (rI at 1 in 385 bases). Transcriptome and proteome analysis in mutant cells revealed no major differences with controls. The rate of transcription and the total amount of cellular RNA also appeared normal. rI accumulation in RNA-and by implication rI production-correlates with the severity of organ dysfunction in ITPase deficiency but the basis of the cellulopathy remains cryptic. While we cannot exclude cumulative minor effects, there are no major anomalies in the production, processing, stability and/or translation of mRNA.

Pathogenicity and selective constraint on variation near splice sites

Mutations that perturb normal pre-mRNA splicing are significant contributors to human disease. We used exome sequencing data from 7833 probands with developmental disorders (DDs) and their unaffected parents, as well as more than 60,000 aggregated exomes from the Exome Aggregation Consortium, to investigate selection around the splice sites and quantify the contribution of splicing mutations to DDs. Patterns of purifying selection, a deficit of variants in highly constrained genes in healthy subjects, and excess de novo mutations in patients highlighted particular positions within and around the consensus splice site of greater functional relevance. By using mutational burden analyses in this large cohort of proband–parent trios, we could estimate in an unbiased manner the relative contributions of mutations at canonical dinucleotides (73%) and flanking noncanonical positions (27%), and calculate the positive predictive value of pathogenicity for different classes of mutations. We identified 18 patients with likely diagnostic de novo mutations in dominant DD-associated genes at noncanonical positions in splice sites. We estimate 35%–40% of pathogenic variants in noncanonical splice site positions are missing from public databases.

Quantifying the contribution of recessive coding variation to developmental disorders

We estimated the genome-wide contribution of recessive coding variation in 6040 families from the Deciphering Developmental Disorders study. The proportion of cases attributable to recessive coding variants was 3.6% in patients of European ancestry, compared with 50% explained by de novo coding mutations. It was higher (31%) in patients with Pakistani ancestry, owing to elevated autozygosity. Half of this recessive burden is attributable to known genes. We identified two genes not previously associated with recessive developmental disorders, KDM5B and EIF3F, and functionally validated them with mouse and cellular models. Our results suggest that recessive coding variants account for a small fraction of currently undiagnosed nonconsanguineous individuals, and that the role of noncoding variants, incomplete penetrance, and polygenic mechanisms need further exploration.

Dysfunction of NaV1.4, a skeletal muscle voltage-gated sodium channel, in sudden infant death syndrome: a case-control study

BACKGROUND

Sudden infant death syndrome (SIDS) is the leading cause of post-neonatal infant death in high-income countries. Central respiratory system dysfunction seems to contribute to these deaths. Excitation that drives contraction of skeletal respiratory muscles is controlled by the sodium channel NaV1.4, which is encoded by the gene SCN4A. Variants in NaV1.4 that directly alter skeletal muscle excitability can cause myotonia, periodic paralysis, congenital myopathy, and myasthenic syndrome. SCN4A variants have also been found in infants with life-threatening apnoea and laryngospasm. We therefore hypothesised that rare, functionally disruptive SCN4A variants might be over-represented in infants who died from SIDS.

METHODS

We did a case-control study, including two consecutive cohorts that included 278 SIDS cases of European ancestry and 729 ethnically matched controls without a history of cardiovascular, respiratory, or neurological disease. We compared the frequency of rare variants in SCN4A between groups (minor allele frequency <0·00005 in the Exome Aggregation Consortium). We assessed biophysical characterisation of the variant channels using a heterologous expression system.

FINDINGS

Four (1·4%) of the 278 infants in the SIDS cohort had a rare functionally disruptive SCN4A variant compared with none (0%) of 729 ethnically matched controls (p=0·0057).

INTERPRETATION

Rare SCN4A variants that directly alter NaV1.4 function occur in infants who had died from SIDS. These variants are predicted to significantly alter muscle membrane excitability and compromise respiratory and laryngeal function. These findings indicate that dysfunction of muscle sodium channels is a potentially modifiable risk factor in a subset of infant sudden deaths.

FUNDING

UK Medical Research Council, the Wellcome Trust, National Institute for Health Research, the British Heart Foundation, Biotronik, Cardiac Risk in the Young, Higher Education Funding Council for England, Dravet Syndrome UK, the Epilepsy Society, the Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National Institutes of Health, and the Mayo Clinic Windland Smith Rice Comprehensive Sudden Cardiac Death Program.

NALCN Dysfunction as a Cause of Disordered Respiratory Rhythm With Central Apnea

The sodium leak channel nonselective protein (NALCN) is a regulator of the pacemaker neurons that are responsible for rhythmic behavior (including respiration), maintaining the resting membrane potential, and are required for action potential production. NALCN-null mice show early death associated with disrupted respiratory rhythms, characterized by frequent and profound apneas. We report 3 children (2 siblings) with compound heterozygous mutations in NALCN associated with developmental impairment, hypotonia, and central sleep-disordered breathing causing apneas. Supplemental oxygen normalized the respiratory rhythm. NALCN mutations have been previously reported to cause severe hypotonia, speech impairment, and cognitive delay as well as infantile neuroaxonal dystrophy and facial dysmorphism. Nonsynonymous changes in the 2 affected extracellular loops may be responsible for the deleterious effect on the stability of the respiratory rhythm. Although oxygen is known to be a stabilizer of respiratory rhythm in central apnea in children, its role in NALCN dysfunction requires further investigation.

De novo mutations in regulatory elements in neurodevelopmental disorders

We previously estimated that 42% of patients with severe developmental disorders carry pathogenic de novo mutations in coding sequences. The role of de novo mutations in regulatory elements affecting genes associated with developmental disorders, or other genes, has been essentially unexplored. We identified de novo mutations in three classes of putative regulatory elements in almost 8,000 patients with developmental disorders. Here we show that de novo mutations in highly evolutionarily conserved fetal brain-active elements are significantly and specifically enriched in neurodevelopmental disorders. We identified a significant twofold enrichment of recurrently mutated elements. We estimate that, genome-wide, 1-3% of patients without a diagnostic coding variant carry pathogenic de novo mutations in fetal brain-active regulatory elements and that only 0.15% of all possible mutations within highly conserved fetal brain-active elements cause neurodevelopmental disorders with a dominant mechanism. Our findings represent a robust estimate of the contribution of de novo mutations in regulatory elements to this genetically heterogeneous set of disorders, and emphasize the importance of combining functional and evolutionary evidence to identify regulatory causes of genetic disorders.

Publisher Correction: BRD4 interacts with NIPBL and BRD4 is mutated in a Cornelia de Lange-like syndrome

In the version of this article initially published, Wendy Bickmore and Madapura Pradeepa were incorrectly not indicated as corresponding authors. The error has been corrected in the HTML and PDF versions of the paper.

High Rate of Recurrent De Novo Mutations in Developmental and Epileptic Encephalopathies

Developmental and epileptic encephalopathy (DEE) is a group of conditions characterized by the co-occurrence of epilepsy and intellectual disability (ID), typically with developmental plateauing or regression associated with frequent epileptiform activity. The cause of DEE remains unknown in the majority of cases. We performed whole-genome sequencing (WGS) in 197 individuals with unexplained DEE and pharmaco-resistant seizures and in their unaffected parents. We focused our attention on de novo mutations (DNMs) and identified candidate genes containing such variants. We sought to identify additional subjects with DNMs in these genes by performing targeted sequencing in another series of individuals with DEE and by mining various sequencing datasets. We also performed meta-analyses to document enrichment of DNMs in candidate genes by leveraging our WGS dataset with those of several DEE and ID series. By combining these strategies, we were able to provide a causal link between DEE and the following genes: NTRK2, GABRB2, CLTC, DHDDS, NUS1, RAB11A, GABBR2, and SNAP25. Overall, we established a molecular diagnosis in 63/197 (32%) individuals in our WGS series. The main cause of DEE in these individuals was de novo point mutations (53/63 solved cases), followed by inherited mutations (6/63 solved cases) and de novo CNVs (4/63 solved cases). De novo missense variants explained a larger proportion of individuals in our series than in other series that were primarily ascertained because of ID. Moreover, these DNMs were more frequently recurrent than those identified in ID series. These observations indicate that the genetic landscape of DEE might be different from that of ID without epilepsy.

Genetic Interactions in Nonsyndromic Orofacial Clefts in Europe—EUROCRAN Study

Background

Nonsyndromic cleft lip with or without cleft palate (nsCL±P) and nonsyndromic cleft palate (nsCP) are caused by a combination of genetic and environmental risk factors. We investigated gene-environment and gene-gene joint effects in a large multicenter study of case-parent triads.

Methods

The nsCL±P or nsCP triads were recruited in 11 European countries between 2001 and 2005. We collected DNA samples from infants and from their mothers and fathers, and mothers completed a questionnaire on exposures, including smoking and folic acid supplement use during pregnancy. We used log-linear regression to estimate relative risks (RRs) and 95% confidence intervals (CIs) for associations between nsCL±P or nsCP and variants in MTHFR, MTHFD1, TGFA, SATB2, and MSX1, stratifying by environmental or genetic factors.

Results

We obtained genotype and exposure data for 728 nsCL±P triads and 292 nsCP triads. In male infants, there was no association between the mother's homozygous MSX1 p(CA) ∗4/∗4 genotype and nsCL±P (RR, 0.98; 95% CI, 0.63–1.54), but this maternal genotype resulted in a doubling of risk for female infants (RR, 2.21; 95% CI, 1.13–4.34). There was evidence suggestive of gene-gene joint-effects between MTHFR-TGFA for nsCP but not for nsCL±P.

Conclusion

Although we chose the genes and their variants and putative joint effects based on associations previously reported in the literature, we replicated few associations. These results do not provide evidence supporting associations between these genes and oral clefts in European populations, although gene-environment and gene-gene interactions could play a role in oral cleft etiology.

A recurrent de novo mutation in ACTG1 causes isolated ocular coloboma

Ocular coloboma (OC) is a defect in optic fissure closure and is a common cause of severe congenital visual impairment. Bilateral OC is primarily genetically determined and shows marked locus heterogeneity. Whole-exome sequencing (WES) was used to analyze 12 trios (child affected with OC and both unaffected parents). This identified de novo mutations in 10 different genes in eight probands. Three of these genes encoded proteins associated with actin cytoskeleton dynamics: ACTG1, TWF1, and LCP1. Proband-only WES identified a second unrelated individual with isolated OC carrying the same ACTG1 allele, encoding p.(Pro70Leu). Both individuals have normal neurodevelopment with no extra-ocular signs of Baraitser-Winter syndrome. We found this mutant protein to be incapable of incorporation into F-actin. The LCP1 and TWF1 variants each resulted in only minor disturbance of actin interactions, and no further plausibly causative variants were identified in these genes on resequencing 380 unrelated individuals with OC.

Novel PEX11B Mutations Extend the Peroxisome Biogenesis Disorder 14B Phenotypic Spectrum and Underscore Congenital Cataract as an Early Feature

Purpose

Peroxisomes perform complex metabolic and catabolic functions essential for normal growth and development. Mutations in 14 genes cause a spectrum of peroxisomal disease in humans. Most recently, PEX11B was associated with an atypical peroxisome biogenesis disorder (PBD) in a single individual. In this study, we identify further PEX11B cases and delineate associated phenotypes.

Methods

Probands from three families underwent next generation sequencing (NGS) for diagnosis of a multisystem developmental disorder. Autozygosity mapping was conducted in one affected sibling pair. ExomeDepth was used to identify copy number variants from NGS data and confirmed by dosage analysis. Biochemical profiling was used to investigate the metabolic signature of the condition.

Results

All patients presented with bilateral cataract at birth but the systemic phenotype was variable, including short stature, skeletal abnormalities, and dysmorphism-features not described in the original case. Next generation sequencing identified biallelic loss-of-function mutations in PEX11B as the underlying cause of disease in each case (PEX11B c.235C>T p.(Arg79Ter) homozygous; PEX11B c.136C>T p.(Arg46Ter) homozygous; PEX11B c.595C>T p.(Arg199Ter) heterozygous, PEX11B ex1-3 del heterozygous). Biochemical studies identified very low plasmalogens in one patient, whilst a mildly deranged very long chain fatty acid profile was found in another.

Conclusions

Our findings expand the phenotypic spectrum of the condition and underscore congenital cataract as the consistent primary presenting feature. We also find that biochemical measurements of peroxisome function may be disturbed in some cases. Furthermore, diagnosis by NGS is proficient and may circumvent the requirement for an invasive skin biopsy for disease identification from fibroblast cells.

Genotype-phenotype correlations in Cornelia de Lange syndrome: Behavioral characteristics and changes with age

Cornelia de Lange syndrome (CdLS) is a multisystem genetic disorder associated with unusual facial features, limb abnormalities, a wide range of health conditions, and intellectual disability. Mutations in five genes that encode (SMC1A, SMC3, RAD21) or regulate (NIPBL, HDAC8) the cohesin complex have been identified in up to 70% of individuals. Genetic cause remains unknown for a proportion of individuals. There is substantial heterogeneity in all aspects of CdLS but very little is known about what predicts phenotypic heterogeneity. In this study, we evaluated genotype-phenotype associations in 34 individuals with CdLS. Participants with NIPBL mutations had significantly lower self help skills and were less likely to have verbal skills relative to those who were negative for the NIPBL mutation. No significant differences were identified between the groups in relation to repetitive behavior, mood, interest and pleasure, challenging behavior, activity, impulsivity, and characteristics of autism spectrum disorder whilst controlling differences in self help skills. Significant correlations indicating lower mood, interest and pleasure, and increased insistence on sameness with older age were identified for those who were NIPBL mutation positive. The findings suggest similarities in the behavioral phenotype between those with and without the NIPBL mutation once differences in self help skills are controlled for. However, there may be subtle differences in the developmental trajectory of these behaviors according to genetic mutation status in CdLS.

Returning genome sequences to research participants: Policy and practice

Despite advances in genomic science stimulating an explosion of literature around returning health-related findings, the possibility of returning entire genome sequences to individual research participants has not been widely considered. Through direct involvement in large-scale translational genomics studies, we have identified a number of logistical challenges that would need to be overcome prior to returning individual genome sequence data, including verifying that the data belong to the requestor and providing appropriate informatics support. In addition, we identify a number of ethico-legal issues that require careful consideration, including returning data to family members, mitigating against unintended consequences, and ensuring appropriate governance. Finally, recognising that there is an opportunity cost to addressing these issues, we make some specific pragmatic suggestions for studies that are considering whether to share individual genomic datasets with individual study participants. If data are shared, research should be undertaken into the personal, familial and societal impact of receiving individual genome sequence data.

Clinical and molecular consequences of disease-associated de novo mutations in SATB2

PURPOSE

To characterize features associated with de novo mutations affecting SATB2 function in individuals ascertained on the basis of intellectual disability.

METHODS

Twenty previously unreported individuals with 19 different SATB2 mutations (11 loss-of-function and 8 missense variants) were studied. Fibroblasts were used to measure mutant protein production. Subcellular localization and mobility of wild-type and mutant SATB2 were assessed using fluorescently tagged protein.

RESULTS

Recurrent clinical features included neurodevelopmental impairment (19/19), absent/near absent speech (16/19), normal somatic growth (17/19), cleft palate (9/19), drooling (12/19), and dental anomalies (8/19). Six of eight missense variants clustered in the first CUT domain. Sibling recurrence due to gonadal mosaicism was seen in one family. A nonsense mutation in the last exon resulted in production of a truncated protein retaining all three DNA-binding domains. SATB2 nuclear mobility was mutation-dependent; p.Arg389Cys in CUT1 increased mobility and both p.Gly515Ser in CUT2 and p.Gln566Lys between CUT2 and HOX reduced mobility. The clinical features in individuals with missense variants were indistinguishable from those with loss of function.

CONCLUSION

SATB2 haploinsufficiency is a common cause of syndromic intellectual disability. When mutant SATB2 protein is produced, the protein appears functionally inactive with a disrupted pattern of chromatin or matrix association.Genet Med advance online publication 02 February 2017.

Clinical features associated with CTNNB1 de novo loss of function mutations in ten individuals

Loss of function mutations in CTNNB1 have been reported in individuals with intellectual disability [MIM #615075] associated with peripheral spasticity, microcephaly and central hypotonia, suggesting a recognisable phenotype associated with haploinsufficiency for this gene. Trio based whole exome sequencing via the Deciphering Developmental Disorders (DDD) study has identified eleven further individuals with de novo loss of function mutations in CTNNB1. Here we report detailed phenotypic information on ten of these. We confirm the features that have been previously described and further delineate the skin and hair findings, including fair skin and fair and sparse hair with unusual patterning.

SMCHD1 mutations associated with a rare muscular dystrophy can also cause isolated arhinia and Bosma arhinia microphthalmia syndrome

Arhinia, or absence of the nose, is a rare malformation of unknown etiology that is often accompanied by ocular and reproductive defects. Sequencing of 40 people with arhinia revealed that 84% of probands harbor a missense mutation localized to a constrained region of SMCHD1 encompassing the ATPase domain. SMCHD1 mutations cause facioscapulohumeral muscular dystrophy type 2 (FSHD2) via a trans-acting loss-of-function epigenetic mechanism. We discovered shared mutations and comparable DNA hypomethylation patterning between these distinct disorders. CRISPR/Cas9-mediated alteration of smchd1 in zebrafish yielded arhinia-relevant phenotypes. Transcriptome and protein analyses in arhinia probands and controls showed no differences in SMCHD1 mRNA or protein abundance but revealed regulatory changes in genes and pathways associated with craniofacial patterning. Mutations in SMCHD1 thus contribute to distinct phenotypic spectra, from craniofacial malformation and reproductive disorders to muscular dystrophy, which we speculate to be consistent with oligogenic mechanisms resulting in pleiotropic outcomes.

The RNA-binding landscape of RBM10 and its role in alternative splicing regulation in models of mouse early development

Mutations in the RNA-binding protein, RBM10, result in a human syndromic form of cleft palate, termed TARP syndrome. A role for RBM10 in alternative splicing regulation has been previously demonstrated in human cell lines. To uncover the cellular functions of RBM10 in a cell line that is relevant to the phenotype observed in TARP syndrome, we used iCLIP to identify its endogenous RNA targets in a mouse embryonic mandibular cell line. We observed that RBM10 binds to pre-mRNAs with significant enrichment in intronic regions, in agreement with a role for this protein in pre-mRNA splicing. In addition to protein-coding transcripts, RBM10 also binds to a variety of cellular RNAs, including non-coding RNAs, such as spliceosomal small nuclear RNAs, U2 and U12. RNA-seq was used to investigate changes in gene expression and alternative splicing in RBM10 KO mouse mandibular cells and also in mouse ES cells. We uncovered a role for RBM10 in the regulation of alternative splicing of common transcripts in both cell lines but also identified cell-type specific events. Importantly, those pre-mRNAs that display changes in alternative splicing also contain RBM10 iCLIP tags, suggesting a direct role of RBM10 in these events. Finally, we show that depletion of RBM10 in mouse ES cells leads to proliferation defects and to gross alterations in their differentiation potential. These results demonstrate a role for RBM10 in the regulation of alternative splicing in two cell models of mouse early development and suggests that mutations in RBM10 could lead to splicing changes that affect normal palate development and cause human disease.

Distinct genetic architectures for syndromic and nonsyndromic congenital heart defects identified by exome sequencing

Congenital heart defects (CHDs) have a neonatal incidence of 0.8-1% (refs. 1,2). Despite abundant examples of monogenic CHD in humans and mice, CHD has a low absolute sibling recurrence risk (∼2.7%), suggesting a considerable role for de novo mutations (DNMs) and/or incomplete penetrance. De novo protein-truncating variants (PTVs) have been shown to be enriched among the 10% of 'syndromic' patients with extra-cardiac manifestations. We exome sequenced 1,891 probands, including both syndromic CHD (S-CHD, n = 610) and nonsyndromic CHD (NS-CHD, n = 1,281). In S-CHD, we confirmed a significant enrichment of de novo PTVs but not inherited PTVs in known CHD-associated genes, consistent with recent findings. Conversely, in NS-CHD we observed significant enrichment of PTVs inherited from unaffected parents in CHD-associated genes. We identified three genome-wide significant S-CHD disorders caused by DNMs in CHD4, CDK13 and PRKD1. Our study finds evidence for distinct genetic architectures underlying the low sibling recurrence risk in S-CHD and NS-CHD.

A Restricted Repertoire of De Novo Mutations in ITPR1 Cause Gillespie Syndrome with Evidence for Dominant-Negative Effect

Gillespie syndrome (GS) is characterized by bilateral iris hypoplasia, congenital hypotonia, non-progressive ataxia, and progressive cerebellar atrophy. Trio-based exome sequencing identified de novo mutations in ITPR1 in three unrelated individuals with GS recruited to the Deciphering Developmental Disorders study. Whole-exome or targeted sequence analysis identified plausible disease-causing ITPR1 mutations in 10/10 additional GS-affected individuals. These ultra-rare protein-altering variants affected only three residues in ITPR1: Glu2094 missense (one de novo, one co-segregating), Gly2539 missense (five de novo, one inheritance uncertain), and Lys2596 in-frame deletion (four de novo). No clinical or radiological differences were evident between individuals with different mutations. ITPR1 encodes an inositol 1,4,5-triphosphate-responsive calcium channel. The homo-tetrameric structure has been solved by cryoelectron microscopy. Using estimations of the degree of structural change induced by known recessive- and dominant-negative mutations in other disease-associated multimeric channels, we developed a generalizable computational approach to indicate the likely mutational mechanism. This analysis supports a dominant-negative mechanism for GS variants in ITPR1. In GS-derived lymphoblastoid cell lines (LCLs), the proportion of ITPR1-positive cells using immunofluorescence was significantly higher in mutant than control LCLs, consistent with an abnormality of nuclear calcium signaling feedback control. Super-resolution imaging supports the existence of an ITPR1-lined nucleoplasmic reticulum. Mice with Itpr1 heterozygous null mutations showed no major iris defects. Purkinje cells of the cerebellum appear to be the most sensitive to impaired ITPR1 function in humans. Iris hypoplasia is likely to result from either complete loss of ITPR1 activity or structure-specific disruption of multimeric interactions.