Narrowing the diagnostic gap: Genomes, episignatures, long-read sequencing, and health economic analyses in an exome-negative intellectual disability cohort

PURPOSE

Genome sequencing (GS)-specific diagnostic rates in prospective tightly ascertained exome sequencing (ES)-negative intellectual disability (ID) cohorts have not been reported extensively.

METHODS

ES, GS, epigenetic signatures, and long-read sequencing diagnoses were assessed in 74 trios with at least moderate ID.

RESULTS

The ES diagnostic yield was 42 of 74 (57%). GS diagnoses were made in 9 of 32 (28%) ES-unresolved families. Repeated ES with a contemporary pipeline on the GS-diagnosed families identified 8 of 9 single-nucleotide variations/copy-number variations undetected in older ES, confirming a GS-unique diagnostic rate of 1 in 32 (3%). Episignatures contributed diagnostic information in 9% with GS corroboration in 1 of 32 (3%) and diagnostic clues in 2 of 32 (6%). A genetic etiology for ID was detected in 51 of 74 (69%) families. Twelve candidate disease genes were identified. Contemporary ES followed by GS cost US$4976 (95% CI: $3704; $6969) per diagnosis and first-line GS at a cost of $7062 (95% CI: $6210; $8475) per diagnosis.

CONCLUSION

Performing GS only in ID trios would be cost equivalent to ES if GS were available at $2435, about a 60% reduction from current prices. This study demonstrates that first-line GS achieves higher diagnostic rate than contemporary ES but at a higher cost.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

De novo variants in the non-coding spliceosomal snRNA gene RNU4-2 are a frequent cause of syndromic neurodevelopmental disorders

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes 1 . Increasingly, large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a novel syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome 2 . We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 119 individuals with NDD. The vast majority of individuals (77.3%) have the same highly recurrent single base-pair insertion (n.64_65insT). We estimate that variants in this region explain 0.41% of individuals with NDD. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to its contiguous counterpart RNU4-1 and other U4 homologs, supporting RNU4-2 's role as the primary U4 transcript in the brain. Overall, this work underscores the importance of non-coding genes in rare disorders. It will provide a diagnosis to thousands of individuals with NDD worldwide and pave the way for the development of effective treatments for these individuals.

Mono and biallelic variants in HCN2 cause severe neurodevelopmental disorders

Hyperpolarization activated Cyclic Nucleotide (HCN) gated channels are crucial for various neurophysiological functions, including learning and sensory functions, and their dysfunction are responsible for brain disorders, such as epilepsy. To date, HCN2 variants have only been associated with mild epilepsy and recently, one monoallelic missense variant has been linked to developmental and epileptic encephalopathy. Here, we expand the phenotypic spectrum of HCN2- related disorders by describing twenty-one additional individuals from fifteen unrelated families carrying HCN2 variants. Seventeen individuals had developmental delay/intellectual disability (DD/ID), two had borderline DD/ID, and one had borderline DD. Ten individuals had epilepsy with DD/ID, with median age of onset of 10 months, and one had epilepsy with normal development. Molecular diagnosis identified thirteen different pathogenic HCN2 variants, including eleven missense variants affecting highly conserved amino acids, one frameshift variant, and one in-frame deletion. Seven variants were monoallelic of which five occurred de novo, one was not maternally inherited, one was inherited from a father with mild learning disabilities, and one was of unknown inheritance. The remaining six variants were biallelic, with four homozygous and two compound heterozygous variants. Functional studies using two-electrode voltage-clamp recordings in Xenopus laevis oocytes were performed on three monoallelic variants, p.(Arg324His), p.(Ala363Val), and p.(Met374Leu), and three biallelic variants, p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp). The p.(Arg324His) variant induced a strong increase of HCN2 conductance, while p.(Ala363Val) and p.(Met374Leu) displayed dominant negative effects, leading to a partial loss of HCN2 channel function. By confocal imaging, we found that the p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp) pathogenic variants impaired membrane trafficking, resulting in a complete loss of HCN2 elicited currents in Xenopus oocytes. Structural 3D-analysis in depolarized and hyperpolarized states of HCN2 channels, revealed that the pathogenic variants p.(His205Gln), p.(Ser409Leu), p.(Arg324Cys), p.(Asn369Ser) and p.(Gly460Asp) modify molecular interactions altering HCN2 function. Taken together, our data broadens the clinical spectrum associated with HCN2 variants, and disclose that HCN2 is involved in developmental encephalopathy with or without epilepsy.

Biallellic variants in CACNA1S cause fetal akinesia sequence, progressive hydrops and stillbirth

Fetal arthrogryposis is a well-recognised ultrasonographic phenotype, caused by both genetic, maternal and extrinsic factors. When present with fetal growth restriction, pulmonary hypoplasia and multiple joint contractures, it is often referred to as fetal akinesia deformation sequence (FADS). Historically, elucidating genetic causes of arthryogryposis/FADS has been challenging; there are now more than 150 genes known to cause arthrogryposis through myopathic, neuromuscular and metabolic pathways affecting fetal movement. FADS is associated with over 400 medical conditions making prenatal diagnosis challenging. Here we present a case of FADS diagnosed at 19 weeks gestation with progression to severe fetal hydrops and stillbirth at 26-weeks gestation. Initial investigations including combined first trimester screening, TORCH (infection) screen and chromosomal microarray were normal. Trio whole exome sequencing (WES) detected compound heterozygous likely pathogenic CACNA1S gene variants associated with autosomal dominant (AD) and autosomal recessive (AR) congenital myopathy and FADS. To our knowledge, this is the first prenatal diagnosis of this condition.

Diagnostic Utility of Genome-wide DNA Methylation Analysis in Genetically Unsolved Developmental and Epileptic Encephalopathies and Refinement of a CHD2 Episignature

Sequence-based genetic testing currently identifies causative genetic variants in ∼50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. Rare epigenetic variations ("epivariants") can drive disease by modulating gene expression at single loci, whereas genome-wide DNA methylation changes can result in distinct "episignature" biomarkers for monogenic disorders in a growing number of rare diseases. Here, we interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 516 individuals with genetically unsolved DEEs who had previously undergone extensive genetic testing. We identified rare differentially methylated regions (DMRs) and explanatory episignatures to discover causative and candidate genetic etiologies in 10 individuals. We then used long-read sequencing to identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and two copy number variants. We also identify pathogenic sequence variants associated with episignatures; some had been missed by previous exome sequencing. Although most DEE genes lack known episignatures, the increase in diagnostic yield for DNA methylation analysis in DEEs is comparable to the added yield of genome sequencing. Finally, we refine an episignature for CHD2 using an 850K methylation array which was further refined at higher CpG resolution using bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate genetic causes as ∼2% (10/516) for unsolved DEE cases.

Integration of EpiSign, facial phenotyping, and likelihood ratio interpretation of clinical abnormalities in the re-classification of an ARID1B missense variant

Heterozygous ARID1B variants result in Coffin-Siris syndrome. Features may include hypoplastic nails, slow growth, characteristic facial features, hypotonia, hypertrichosis, and sparse scalp hair. Most reported cases are due to ARID1B loss of function variants. We report a boy with developmental delay, feeding difficulties, aspiration, recurrent respiratory infections, slow growth, and hypotonia without a clinical diagnosis, where a previously unreported ARID1B missense variant was classified as a variant of uncertain significance. The pathogenicity of this variant was refined through combined methodologies including genome-wide methylation signature analysis (EpiSign), Machine Learning (ML) facial phenotyping, and LIRICAL. Trio exome sequencing and EpiSign were performed. ML facial phenotyping compared facial images using FaceMatch and GestaltMatcher to syndrome-specific libraries to prioritize the trio exome bioinformatic pipeline gene list output. Phenotype-driven variant prioritization was performed with LIRICAL. A de novo heterozygous missense variant, ARID1B p.(Tyr1268His), was reported as a variant of uncertain significance. The ACMG classification was refined to likely pathogenic by a supportive methylation signature, ML facial phenotyping, and prioritization through LIRICAL. The ARID1B genotype-phenotype has been expanded through an extended analysis of missense variation through genome-wide methylation signatures, ML facial phenotyping, and likelihood-ratio gene prioritization.

Stretch-activated ion channel TMEM63B associates with developmental and epileptic encephalopathies and progressive neurodegeneration

By converting physical forces into electrical signals or triggering intracellular cascades, stretch-activated ion channels allow the cell to respond to osmotic and mechanical stress. Knowledge of the pathophysiological mechanisms underlying associations of stretch-activated ion channels with human disease is limited. Here, we describe 17 unrelated individuals with severe early-onset developmental and epileptic encephalopathy (DEE), intellectual disability, and severe motor and cortical visual impairment associated with progressive neurodegenerative brain changes carrying ten distinct heterozygous variants of TMEM63B, encoding for a highly conserved stretch-activated ion channel. The variants occurred de novo in 16/17 individuals for whom parental DNA was available and either missense, including the recurrent p.Val44Met in 7/17 individuals, or in-frame, all affecting conserved residues located in transmembrane regions of the protein. In 12 individuals, hematological abnormalities co-occurred, such as macrocytosis and hemolysis, requiring blood transfusions in some. We modeled six variants (p.Val44Met, p.Arg433His, p.Thr481Asn, p.Gly580Ser, p.Arg660Thr, and p.Phe697Leu), each affecting a distinct transmembrane domain of the channel, in transfected Neuro2a cells and demonstrated inward leak cation currents across the mutated channel even in isotonic conditions, while the response to hypo-osmotic challenge was impaired, as were the Ca2+ transients generated under hypo-osmotic stimulation. Ectopic expression of the p.Val44Met and p.Gly580Cys variants in Drosophila resulted in early death. TMEM63B-associated DEE represents a recognizable clinicopathological entity in which altered cation conductivity results in a severe neurological phenotype with progressive brain damage and early-onset epilepsy associated with hematological abnormalities in most individuals.

Identification of Regulatory Molecular "Hot Spots" for LH/PLOD Collagen Glycosyltransferase Activity

Hydroxylysine glycosylations are post-translational modifications (PTMs) essential for the maturation and homeostasis of fibrillar and non-fibrillar collagen molecules. The multifunctional collagen lysyl hydroxylase 3 (LH3/PLOD3) and the collagen galactosyltransferase GLT25D1 are the human enzymes that have been identified as being responsible for the glycosylation of collagen lysines, although a precise description of the contribution of each enzyme to these essential PTMs has not yet been provided in the literature. LH3/PLOD3 is thought to be capable of performing two chemically distinct collagen glycosyltransferase reactions using the same catalytic site: an inverting beta-1,O-galactosylation of hydroxylysines (Gal-T) and a retaining alpha-1,2-glucosylation of galactosyl hydroxylysines (Glc-T). In this work, we have combined indirect luminescence-based assays with direct mass spectrometry-based assays and molecular structure studies to demonstrate that LH3/PLOD3 only has Glc-T activity and that GLT25D1 only has Gal-T activity. Structure-guided mutagenesis confirmed that the Glc-T activity is defined by key residues in the first-shell environment of the glycosyltransferase catalytic site as well as by long-range contributions from residues within the same glycosyltransferase (GT) domain. By solving the molecular structures and characterizing the interactions and solving the molecular structures of human LH3/PLOD3 in complex with different UDP-sugar analogs, we show how these studies could provide insights for LH3/PLOD3 glycosyltransferase inhibitor development. Collectively, our data provide new tools for the direct investigation of collagen hydroxylysine PTMs and a comprehensive overview of the complex network of shapes, charges, and interactions that enable LH3/PLOD3 glycosyltransferase activities, expanding the molecular framework and facilitating an improved understanding and manipulation of glycosyltransferase functions in biomedical applications.

The clinical and molecular spectrum of the KDM6B-related neurodevelopmental disorder

De novo variants are a leading cause of neurodevelopmental disorders (NDDs), but because every monogenic NDD is different and usually extremely rare, it remains a major challenge to understand the complete phenotype and genotype spectrum of any morbid gene. According to OMIM, heterozygous variants in KDM6B cause "neurodevelopmental disorder with coarse facies and mild distal skeletal abnormalities." Here, by examining the molecular and clinical spectrum of 85 reported individuals with mostly de novo (likely) pathogenic KDM6B variants, we demonstrate that this description is inaccurate and potentially misleading. Cognitive deficits are seen consistently in all individuals, but the overall phenotype is highly variable. Notably, coarse facies and distal skeletal anomalies, as defined by OMIM, are rare in this expanded cohort while other features are unexpectedly common (e.g., hypotonia, psychosis, etc.). Using 3D protein structure analysis and an innovative dual Drosophila gain-of-function assay, we demonstrated a disruptive effect of 11 missense/in-frame indels located in or near the enzymatic JmJC or Zn-containing domain of KDM6B. Consistent with the role of KDM6B in human cognition, we demonstrated a role for the Drosophila KDM6B ortholog in memory and behavior. Taken together, we accurately define the broad clinical spectrum of the KDM6B-related NDD, introduce an innovative functional testing paradigm for the assessment of KDM6B variants, and demonstrate a conserved role for KDM6B in cognition and behavior. Our study demonstrates the critical importance of international collaboration, sharing of clinical data, and rigorous functional analysis of genetic variants to ensure correct disease diagnosis for rare disorders.

De novo MCM6 variants in neurodevelopmental disorders: a recognizable phenotype related to zinc binding residues

The minichromosome maintenance (MCM) complex acts as a DNA helicase during DNA replication, and thereby regulates cell cycle progression and proliferation. In addition, MCM-complex components localize to centrosomes and play an independent role in ciliogenesis. Pathogenic variants in genes coding for MCM components and other DNA replication factors have been linked to growth and developmental disorders as Meier-Gorlin syndrome and Seckel syndrome. Trio exome/genome sequencing identified the same de novo MCM6 missense variant p.(Cys158Tyr) in two unrelated individuals that presented with overlapping phenotypes consisting of intra-uterine growth retardation, short stature, congenital microcephaly, endocrine features, developmental delay and urogenital anomalies. The identified variant affects a zinc binding cysteine in the MCM6 zinc finger signature. This domain, and specifically cysteine residues, are essential for MCM-complex dimerization and the induction of helicase activity, suggesting a deleterious effect of this variant on DNA replication. Fibroblasts derived from the two affected individuals showed defects both in ciliogenesis and cell proliferation. We additionally traced three unrelated individuals with de novo MCM6 variants in the oligonucleotide binding (OB)-fold domain, presenting with variable (neuro)developmental features including autism spectrum disorder, developmental delay, and epilepsy. Taken together, our findings implicate de novo MCM6 variants in neurodevelopmental disorders. The clinical features and functional defects related to the zinc binding residue resemble those observed in syndromes related to other MCM components and DNA replication factors, while de novo OB-fold domain missense variants may be associated with more variable neurodevelopmental phenotypes. These data encourage consideration of MCM6 variants in the diagnostic arsenal of NDD.

A mutational hotspot in AMOTL1 defines a new syndrome of orofacial clefting, cardiac anomalies, and tall stature

AMOTL1 encodes angiomotin-like protein 1, an actin-binding protein that regulates cell polarity, adhesion, and migration. The role of AMOTL1 in human disease is equivocal. We report a large cohort of individuals harboring heterozygous AMOTL1 variants and define a core phenotype of orofacial clefting, congenital heart disease, tall stature, auricular anomalies, and gastrointestinal manifestations in individuals with variants in AMOTL1 affecting amino acids 157-161, a functionally undefined but highly conserved region. Three individuals with AMOTL1 variants outside this region are also described who had variable presentations with orofacial clefting and multi-organ disease. Our case cohort suggests that heterozygous missense variants in AMOTL1, most commonly affecting amino acid residues 157-161, define a new orofacial clefting syndrome, and indicates an important functional role for this undefined region.

Dominant-negative variants in CBX1 cause a neurodevelopmental disorder

PURPOSE

This study aimed to establish variants in CBX1, encoding heterochromatin protein 1β (HP1β), as a cause of a novel syndromic neurodevelopmental disorder.

METHODS

Patients with CBX1 variants were identified, and clinician researchers were connected using GeneMatcher and physician referrals. Clinical histories were collected from each patient. To investigate the pathogenicity of identified variants, we performed in vitro cellular assays, neurobehavioral and cytological analyses of neuronal cells obtained from newly generated Cbx1 mutant mouse lines.

RESULTS

In three unrelated individuals with developmental delay, hypotonia, and autistic features, we identified heterozygous de novo variants in CBX1. The identified variants were in the chromodomain, the functional domain of HP1β, which mediates interactions with chromatin. Cbx1 chromodomain mutant mice displayed increased latency-to-peak response, suggesting the possibility of synaptic delay or myelination deficits. Cytological and chromatin immunoprecipitation experiments confirmed the reduction of mutant HP1β binding to heterochromatin, while HP1β interactome analysis demonstrated that the majority of HP1β-interacting proteins remained unchanged between the wild-type and mutant HP1β.

CONCLUSION

These collective findings confirm the role of CBX1 in developmental disabilities through the disruption of HP1β chromatin binding during neurocognitive development. As HP1β forms homodimers and heterodimers, mutant HP1β likely sequesters wild-type HP1β and other HP1 proteins, exerting dominant-negative effects.

Biallelic variants in TUBGCP6 result in microcephaly and chorioretinopathy 1: Report of four cases and a literature review

Autosomal recessive microcephaly and chorioretinopathy-1 (MCCRP1) is a rare Mendelian disorder resulting from biallelic loss of function variants in Tubulin-Gamma Complex Associated Protein 6 (TUBGCP6, MIM#610053). Clinical features of this disorder include microcephaly, cognitive impairment, dysmorphic features, and variable ophthalmological anomalies including chorioretinopathy. Microcephaly can be recognized prenatally and visual impairment becomes evident during the first year of life. The clinical presentation resembles the findings in some acquired conditions such as congenital toxoplasmosis and cytomegalovirus infections; thus, it is important to recognize and diagnose this syndrome in view of its impact on patient health management and familial reproductive plans. To date, only seven molecularly confirmed patients from five unrelated families have been reported. We report an additional four unrelated patients with TUBGCP6 variants including one prenatal diagnosis and review the clinical phenotypes and genotypes of all the known cases. This report expands the molecular and phenotypic spectrum of TUBGCP6 and includes additional prenatal findings associated with MCCRP1.

Bi-allelic TTI1 variants cause an autosomal-recessive neurodevelopmental disorder with microcephaly

Telomere maintenance 2 (TELO2), Tel2 interacting protein 2 (TTI2), and Tel2 interacting protein 1 (TTI1) are the three components of the conserved Triple T (TTT) complex that modulates activity of phosphatidylinositol 3-kinase-related protein kinases (PIKKs), including mTOR, ATM, and ATR, by regulating the assembly of mTOR complex 1 (mTORC1). The TTT complex is essential for the expression, maturation, and stability of ATM and ATR in response to DNA damage. TELO2- and TTI2-related bi-allelic autosomal-recessive (AR) encephalopathies have been described in individuals with moderate to severe intellectual disability (ID), short stature, postnatal microcephaly, and a movement disorder (in the case of variants within TELO2). We present clinical, genomic, and functional data from 11 individuals in 9 unrelated families with bi-allelic variants in TTI1. All present with ID, and most with microcephaly, short stature, and a movement disorder. Functional studies performed in HEK293T cell lines and fibroblasts and lymphoblastoid cells derived from 4 unrelated individuals showed impairment of the TTT complex and of mTOR pathway activity which is improved by treatment with Rapamycin. Our data delineate a TTI1-related neurodevelopmental disorder and expand the group of disorders related to the TTT complex.

Australian Genomics: Outcomes of a 5-year national program to accelerate the integration of genomics in healthcare

Australian Genomics is a national collaborative partnership of more than 100 organizations piloting a whole-of-system approach to integrating genomics into healthcare, based on federation principles. In the first five years of operation, Australian Genomics has evaluated the outcomes of genomic testing in more than 5,200 individuals across 19 rare disease and cancer flagship studies. Comprehensive analyses of the health economic, policy, ethical, legal, implementation and workforce implications of incorporating genomics in the Australian context have informed evidence-based change in policy and practice, resulting in national government funding and equity of access for a range of genomic tests. Simultaneously, Australian Genomics has built national skills, infrastructure, policy, and data resources to enable effective data sharing to drive discovery research and support improvements in clinical genomic delivery.

De novo putative loss-of-function variants in TAF4 are associated with a neuro-developmental disorder

TATA-binding protein associated factor 4 (TAF4) is a subunit of the Transcription Factor IID (TFIID) complex, a central player in transcription initiation. Other members of this multimeric complex have been implicated previously as monogenic disease genes in human developmental disorders. TAF4 has not been described to date as a monogenic disease gene. We here present a cohort of eight individuals, each carrying de novo putative loss-of-function (pLoF) variants in TAF4 and expressing phenotypes consistent with a neuro-developmental disorder (NDD). Common features include intellectual disability, abnormal behavior, and facial dysmorphisms. We propose TAF4 as a novel dominant disease gene for NDD, and coin this novel disorder "TAF4-related NDD" (T4NDD). We place T4NDD in the context of other disorders related to TFIID subunits, revealing shared features of T4NDD with other TAF-opathies.

The role of exome sequencing in childhood interstitial or diffuse lung disease

BACKGROUND

Children's interstitial and diffuse lung disease (chILD) is a complex heterogeneous group of lung disorders. Gene panel approaches have a reported diagnostic yield of ~ 12%. No data currently exist using trio exome sequencing as the standard diagnostic modality. We assessed the diagnostic utility of using trio exome sequencing in chILD. We prospectively enrolled children meeting specified clinical criteria between 2016 and 2020 from 16 Australian hospitals. Exome sequencing was performed with analysis of an initial gene panel followed by trio exome analysis. A subset of critically ill infants underwent ultra-rapid trio exome sequencing as first-line test.

RESULTS

36 patients [median (range) age 0.34 years (0.02-11.46); 11F] were recruited from multiple States and Territories. Five patients had clinically significant likely pathogenic/pathogenic variants (RARB, RPL15, CTCF, RFXANK, TBX4) and one patient had a variant of uncertain significance (VIP) suspected to contribute to their clinical phenotype, with VIP being a novel gene candidate.

CONCLUSIONS

Trio exomes (6/36; 16.7%) had a better diagnostic rate than gene panel (1/36; 2.8%), due to the ability to consider a broader range of underlying conditions. However, the aetiology of chILD in most cases remained undetermined, likely reflecting the interplay between low penetrant genetic and environmental factors.

Heterozygous NOTCH1 variants cause CNS immune activation and microangiopathy

NOTCH1 belongs to the NOTCH family of proteins that regulate cell fate and inflammatory responses. Somatic and germline NOTCH1 variants have been implicated in cancer, Adams-Oliver syndrome and cardiovascular defects. We describe seven unrelated patients grouped by the presence of leukoencephalopathy with calcifications and heterozygous de novo gain-of-function variants in NOTCH1. Immunologic profiling showed upregulated CSF IP-10, a cytokine secreted downstream of NOTCH1 signaling. Autopsy revealed extensive leukoencephalopathy and microangiopathy with vascular calcifications. This evidence implicates that heterozygous gain-of-function variants in NOTCH1 lead to a chronic CNS inflammatory response resulting in a calcifying microangiopathy with leukoencephalopathy. This article is protected by copyright. All rights reserved.

Functional correlation of genome-wide DNA methylation profiles in genetic neurodevelopmental disorders

An expanding range of genetic syndromes are characterized by genome-wide disruptions in DNA methylation profiles referred to as episignatures. Episignatures are distinct, highly sensitive and specific biomarkers that have recently been applied in clinical diagnosis of genetic syndromes. Episignatures are contained within the broader disorder-specific genome-wide DNA methylation changes which can share significant overlap amongst different conditions. In this study we performed functional genomic assessment and comparison of disorder-specific and overlapping genome-wide DNA methylation changes related to 65 genetic syndromes with previously described episignatures. We demonstrate evidence of disorder-specific and recurring genome-wide differentially methylated probes (DMPs) and regions (DMRs). The overall distribution of DMPs and DMRs across the majority of the neurodevelopmental genetic syndromes analyzed showed substantial enrichment in gene promoters and CpG islands, and under-representation of the more variable intergenic regions. Analysis showed significant enrichment of the DMPs and DMRs in gene pathways and processes related to neurodevelopment, including neurogenesis, synaptic signaling and synaptic transmission. This study expands beyond the diagnostic utility of DNA methylation episignatures by demonstrating correlation between the function of the mutated genes and the consequent genomic DNA methylation profiles as a key functional element in the molecular etiology of genetic neurodevelopmental disorders. This article is protected by copyright. All rights reserved.

Prenatal phenotyping: A community effort to enhance the Human Phenotype Ontology

Technological advances in both genome sequencing and prenatal imaging are increasing our ability to accurately recognize and diagnose Mendelian conditions prenatally. Phenotype-driven early genetic diagnosis of fetal genetic disease can help to strategize treatment options and clinical preventive measures during the perinatal period, to plan in utero therapies, and to inform parental decision-making. Fetal phenotypes of genetic diseases are often unique and at present are not well understood; more comprehensive knowledge about prenatal phenotypes and computational resources have an enormous potential to improve diagnostics and translational research. The Human Phenotype Ontology (HPO) has been widely used to support diagnostics and translational research in human genetics. To better support prenatal usage, the HPO consortium conducted a series of workshops with a group of domain experts in a variety of medical specialties, diagnostic techniques, as well as diseases and phenotypes related to prenatal medicine, including perinatal pathology, musculoskeletal anomalies, neurology, medical genetics, hydrops fetalis, craniofacial malformations, cardiology, neonatal-perinatal medicine, fetal medicine, placental pathology, prenatal imaging, and bioinformatics. We expanded the representation of prenatal phenotypes in HPO by adding 95 new phenotype terms under the Abnormality of prenatal development or birth (HP:0001197) grouping term, and revised definitions, synonyms, and disease annotations for most of the 152 terms that existed before the beginning of this effort. The expansion of prenatal phenotypes in HPO will support phenotype-driven prenatal exome and genome sequencing for precision genetic diagnostics of rare diseases to support prenatal care.

Germline variants in tumor suppressor FBXW7 lead to impaired ubiquitination and a neurodevelopmental syndrome

Neurodevelopmental disorders are highly heterogenous conditions resulting from abnormalities of brain architecture and/or function. FBXW7 (F-box and WD-repeat-domain-containing 7), a recognized developmental regulator and tumor suppressor, has been shown to regulate cell-cycle progression and cell growth and survival by targeting substrates including CYCLIN E1/2 and NOTCH for degradation via the ubiquitin proteasome system. We used a genotype-first approach and global data-sharing platforms to identify 35 individuals harboring de novo and inherited FBXW7 germline monoallelic chromosomal deletions and nonsense, frameshift, splice-site, and missense variants associated with a neurodevelopmental syndrome. The FBXW7 neurodevelopmental syndrome is distinguished by global developmental delay, borderline to severe intellectual disability, hypotonia, and gastrointestinal issues. Brain imaging detailed variable underlying structural abnormalities affecting the cerebellum, corpus collosum, and white matter. A crystal-structure model of FBXW7 predicted that missense variants were clustered at the substrate-binding surface of the WD40 domain and that these might reduce FBXW7 substrate binding affinity. Expression of recombinant FBXW7 missense variants in cultured cells demonstrated impaired CYCLIN E1 and CYCLIN E2 turnover. Pan-neuronal knockdown of the Drosophila ortholog, archipelago, impaired learning and neuronal function. Collectively, the data presented herein provide compelling evidence of an F-Box protein-related, phenotypically variable neurodevelopmental disorder associated with monoallelic variants in FBXW7.

Neonatal-lethal dilated cardiomyopathy due to a homozygous LMOD2 donor splice-site variant

Dilated cardiomyopathy (DCM) is characterized by cardiac enlargement and impaired ventricular contractility leading to heart failure. A single report identified variants in leiomodin-2 (LMOD2) as a cause of neonatally-lethal DCM. Here, we describe two siblings with DCM who died shortly after birth due to heart failure. Exome sequencing identified a homozygous LMOD2 variant in both siblings, (GRCh38)chr7:g.123656237G > A; NM_207163.2:c.273 + 1G > A, ablating the donor 5' splice-site of intron-1. Pre-mRNA splicing studies and western blot analysis on cDNA derived from proband cardiac tissue, MyoD-transduced proband skin fibroblasts and HEK293 cells transfected with LMOD2 gene constructs established variant-associated absence of canonically spliced LMOD2 mRNA and full-length LMOD2 protein. Immunostaining of proband heart tissue unveiled abnormally short actin-thin filaments. Our data are consistent with LMOD2 c.273 + 1G > A abolishing/reducing LMOD2 transcript expression by: (1) variant-associated perturbation in initiation of transcription due to ablation of the intron-1 donor; and/or (2) degradation of aberrant LMOD2 transcripts (resulting from use of alternative transcription start-sites or cryptic splice-sites) by nonsense-mediated decay. LMOD2 expression is critical for life and the absence of LMOD2 is associated with thin filament shortening and severe cardiac contractile dysfunction. This study describes the first splice-site variant in LMOD2 and confirms the role of LMOD2 variants in DCM.

Novel diagnostic DNA methylation episignatures expand and refine the epigenetic landscapes of Mendelian disorders

Overlapping clinical phenotypes and an expanding breadth and complexity of genomic associations are a growing challenge in the diagnosis and clinical management of Mendelian disorders. The functional consequences and clinical impacts of genomic variation may involve unique, disorder-specific, genomic DNA methylation episignatures. In this study, we describe 19 novel episignature disorders and compare the findings alongside 38 previously established episignatures for a total of 57 episignatures associated with 65 genetic syndromes. We demonstrate increasing resolution and specificity ranging from protein complex, gene, sub-gene, protein domain, and even single nucleotide-level Mendelian episignatures. We show the power of multiclass modeling to develop highly accurate and disease-specific diagnostic classifiers. This study significantly expands the number and spectrum of disorders with detectable DNA methylation episignatures, improves the clinical diagnostic capabilities through the resolution of unsolved cases and the reclassification of variants of unknown clinical significance, and provides further insight into the molecular etiology of Mendelian conditions.

Standardized practices for RNA diagnostics using clinically accessible specimens reclassifies 75% of putative splicing variants

PURPOSE

Genetic variants causing aberrant premessenger RNA splicing are increasingly being recognized as causal variants in genetic disorders. In this study, we devise standardized practices for polymerase chain reaction (PCR)-based RNA diagnostics using clinically accessible specimens (blood, fibroblasts, urothelia, biopsy).

METHODS

A total of 74 families with diverse monogenic conditions (31% prenatal-congenital onset, 47% early childhood, and 22% teenage-adult onset) were triaged into PCR-based RNA testing, with comparative RNA sequencing for 19 cases.

RESULTS

Informative RNA assay data were obtained for 96% of cases, enabling variant reclassification for 75% variants that can be used for genetic counseling (71%), to inform clinical care (32%) and prenatal counseling (41%). Variant-associated mis-splicing was highly reproducible for 28 cases with samples from ≥2 affected individuals or heterozygotes and 10 cases with ≥2 biospecimens. PCR amplicons encompassing another segregated heterozygous variant was vital for clinical interpretation of 22 of 79 variants to phase RNA splicing events and discern complete from partial mis-splicing.

CONCLUSION

RNA diagnostics enabled provision of a genetic diagnosis for 64% of recruited cases. PCR-based RNA diagnostics has capacity to analyze 81.3% of clinically significant genes, with long amplicons providing an advantage over RNA sequencing to phase RNA splicing events. The Australasian Consortium for RNA Diagnostics (SpliceACORD) provide clinically-endorsed, standardized protocols and recommendations for interpreting RNA assay data.

Clinically Responsive Genomic Analysis Pipelines: Elements to Improve Detection Rate and Efficiency

Massively parallel sequencing has markedly improved mendelian diagnostic rates. This study assessed the effects of custom alterations to a diagnostic genomic bioinformatic pipeline in response to clinical need and derived practice recommendations relative to diagnostic rates and efficiency. The Genomic Annotation and Interpretation Application (GAIA) bioinformatics pipeline was designed to detect panel, exome, and genome sample integrity and prioritize gene variants in mendelian disorders. Reanalysis of selected negative cases was performed after improvements to the pipeline. GAIA improvements and their effect on sensitivity are described, including addition of a PubMed search for gene-disease associations not in the Online Mendelian Inheritance of Man database, inclusion of a process for calling low-quality variants (known as QPatch), and gene symbol nomenclature consistency checking. The new pipeline increased the diagnostic rate and reduced staff costs, resulting in a saving of US$844.34 per additional diagnosis. Recommendations for genomic analysis pipeline requirements are summarized. Clinically responsive bioinformatics pipeline improvements increase diagnostic sensitivity and increase cost-effectiveness.

Different types of disease-causing noncoding variants revealed by genomic and gene expression analyses in families with X-linked intellectual disability

The pioneering discovery research of X-linked intellectual disability (XLID) genes has benefitted thousands of individuals worldwide; however, approximately 30% of XLID families still remain unresolved. We postulated that noncoding variants that affect gene regulation or splicing may account for the lack of a genetic diagnosis in some cases. Detecting pathogenic, gene-regulatory variants with the same sensitivity and specificity as structural and coding variants is a major challenge for Mendelian disorders. Here, we describe three pedigrees with suggestive XLID where distinctive phenotypes associated with known genes guided the identification of three different noncoding variants. We used comprehensive structural, single-nucleotide, and repeat expansion analyses of genome sequencing. RNA-Seq from patient-derived cell lines, reverse-transcription polymerase chain reactions, Western blots, and reporter gene assays were used to confirm the functional effect of three fundamentally different classes of pathogenic noncoding variants: a retrotransposon insertion, a novel intronic splice donor, and a canonical splice variant of an untranslated exon. In one family, we excluded a rare coding variant in ARX, a known XLID gene, in favor of a regulatory noncoding variant in OFD1 that correlated with the clinical phenotype. Our results underscore the value of genomic research on unresolved XLID families to aid novel, pathogenic noncoding variant discovery.

Multisystem inflammation and susceptibility to viral infections in human ZNFX1 deficiency

BACKGROUND

Recognition of viral nucleic acids is one of the primary triggers for a type I interferon-mediated antiviral immune response. Inborn errors of type I interferon immunity can be associated with increased inflammation and/or increased susceptibility to viral infections as a result of dysbalanced interferon production. NFX1-type zinc finger-containing 1 (ZNFX1) is an interferon-stimulated double-stranded RNA sensor that restricts the replication of RNA viruses in mice. The role of ZNFX1 in the human immune response is not known.

OBJECTIVE

We studied 15 patients from 8 families with an autosomal recessive immunodeficiency characterized by severe infections by both RNA and DNA viruses and virally triggered inflammatory episodes with hemophagocytic lymphohistiocytosis-like disease, early-onset seizures, and renal and lung disease.

METHODS

Whole exome sequencing was performed on 13 patients from 8 families. We investigated the transcriptome, posttranscriptional regulation of interferon-stimulated genes (ISGs) and predisposition to viral infections in primary cells from patients and controls stimulated with synthetic double-stranded nucleic acids.

RESULTS

Deleterious homozygous and compound heterozygous ZNFX1 variants were identified in all 13 patients. Stimulation of patient-derived primary cells with synthetic double-stranded nucleic acids was associated with a deregulated pattern of expression of ISGs and alterations in the half-life of the mRNA of ISGs and also associated with poorer clearance of viral infections by monocytes.

CONCLUSION

ZNFX1 is an important regulator of the response to double-stranded nucleic acids stimuli following viral infections. ZNFX1 deficiency predisposes to severe viral infections and a multisystem inflammatory disease.

Pathogenic variants in CDH11 impair cell adhesion and cause Teebi hypertelorism syndrome

Teebi hypertelorism syndrome (THS; OMIM 145420) is a rare craniofacial disorder characterized by hypertelorism, prominent forehead, short nose with broad or depressed nasal root. Some cases of THS have been attributed to SPECC1L variants. Homozygous variants in CDH11 truncating the transmembrane and intracellular domains have been implicated in Elsahy-Waters syndrome (EWS; OMIM 211380) with hypertelorism. We report THS due to CDH11 heterozygous missense variants on 19 subjects from 9 families. All affected residues in the extracellular region of Cadherin-11 (CHD11) are highly conserved across vertebrate species and classical cadherins. Six of the variants that cluster around the EC2-EC3 and EC3-EC4 linker regions are predicted to affect Ca²⁺ binding that is required for cadherin stability. Two of the additional variants [c.164G > C, p.(Trp55Ser) and c.418G > A, p.(Glu140Lys)] are also notable as they are predicted to directly affect trans-homodimer formation. Immunohistochemical study demonstrates that CDH11 is strongly expressed in human facial mesenchyme. Using multiple functional assays, we show that five variants from the EC1, EC2-EC3 linker, and EC3 regions significantly reduced the cell-substrate trans adhesion activity and one variant from EC3-EC4 linker results in changes in cell morphology, focal adhesion, and migration, suggesting dominant negative effect. Characteristic features in this cohort included depressed nasal root, cardiac and umbilical defects. These features distinguished this phenotype from that seen in SPECC1L-related hypertelorism syndrome and CDH11-related EWS. Our results demonstrate heterozygous variants in CDH11, which decrease cell-cell adhesion and increase cell migratory behavior, cause a form of THS, as termed CDH11-related THS.

ClinSV: clinical grade structural and copy number variant detection from whole genome sequencing data

Whole genome sequencing (WGS) has the potential to outperform clinical microarrays for the detection of structural variants (SV) including copy number variants (CNVs), but has been challenged by high false positive rates. Here we present ClinSV, a WGS based SV integration, annotation, prioritization, and visualization framework, which identified 99.8% of simulated pathogenic ClinVar CNVs > 10 kb and 11/11 pathogenic variants from matched microarrays. The false positive rate was low (1.5-4.5%) and reproducibility high (95-99%). In clinical practice, ClinSV identified reportable variants in 22 of 485 patients (4.7%) of which 35-63% were not detectable by current clinical microarray designs. ClinSV is available at https://github.com/KCCG/ClinSV .

Diagnostic Yield of Whole Genome Sequencing After Nondiagnostic Exome Sequencing or Gene Panel in Developmental and Epileptic Encephalopathies

OBJECTIVE

To assess the benefits and limitations of whole genome sequencing (WGS) compared to exome sequencing (ES) or multigene panel (MGP) in the molecular diagnosis of developmental and epileptic encephalopathies (DEE).

METHODS

We performed WGS of 30 comprehensively phenotyped DEE patient trios that were undiagnosed after first-tier testing, including chromosomal microarray and either research ES (n = 15) or diagnostic MGP (n = 15).

RESULTS

Eight diagnoses were made in the 15 individuals who received prior ES (53%): 3 individuals had complex structural variants; 5 had ES-detectable variants, which now had additional evidence for pathogenicity. Eleven diagnoses were made in the 15 MGP-negative individuals (68%); the majority (n = 10) involved genes not included in the panel, particularly in individuals with postneonatal onset of seizures and those with more complex presentations including movement disorders, dysmorphic features, or multiorgan involvement. A total of 42% of diagnoses were autosomal recessive or X-chromosome linked.

CONCLUSION

WGS was able to improve diagnostic yield over ES primarily through the detection of complex structural variants (n = 3). The higher diagnostic yield was otherwise better attributed to the power of re-analysis rather than inherent advantages of the WGS platform. Additional research is required to assist in the assessment of pathogenicity of novel noncoding and complex structural variants and further improve diagnostic yield for patients with DEE and other neurogenetic disorders.

The expanding LARS2 phenotypic spectrum: HLASA, Perrault syndrome with leukodystrophy, and mitochondrial myopathy

LARS2 variants are associated with Perrault syndrome, characterized by premature ovarian failure and hearing loss, and with an infantile lethal multisystem disorder: Hydrops, lactic acidosis, sideroblastic anemia (HLASA) in one individual. Recently we reported LARS2 deafness with (ovario) leukodystrophy. Here we describe five patients with a range of phenotypes, in whom we identified biallelic LARS2 variants: three patients with a HLASA-like phenotype, an individual with Perrault syndrome whose affected siblings also had leukodystrophy, and an individual with a reversible mitochondrial myopathy, lactic acidosis, and developmental delay. Three HLASA cases from two unrelated families were identified. All were males with genital anomalies. Two survived multisystem disease in the neonatal period; both have developmental delay and hearing loss. A 55-year old male with deafness has not displayed neurological symptoms while his female siblings with Perrault syndrome developed leukodystrophy and died in their 30s. Analysis of muscle from a child with a reversible myopathy showed reduced LARS2 and mitochondrial complex I levels, and an unusual form of degeneration. Analysis of recombinant LARS2 variant proteins showed they had reduced aminoacylation efficiency, with HLASA-associated variants having the most severe effect. A broad phenotypic spectrum should be considered in association with LARS2 variants.

Expansion of the phenotypic spectrum of de novo missense variants in kinesin family member 1A (KIF1A)

Defects in the motor domain of kinesin family member 1A (KIF1A), a neuron-specific ATP-dependent anterograde axonal transporter of synaptic cargo, are well-recognized to cause a spectrum of neurological conditions, commonly known as KIF1A-associated neurological disorders (KAND). Here, we report one mutation-negative female with classic Rett syndrome (RTT) harboring a de novo heterozygous novel variant [NP_001230937.1:p.(Asp248Glu)] in the highly conserved motor domain of KIF1A. In addition, three individuals with severe neurodevelopmental disorder along with clinical features overlapping with KAND are also reported carrying de novo heterozygous novel [NP_001230937.1:p.(Cys92Arg) and p.(Pro305Leu)] or previously reported [NP_001230937.1:p.(Thr99Met)] variants in KIF1A. In silico tools predicted these variants to be likely pathogenic, and 3D molecular modeling predicted defective ATP hydrolysis and/or microtubule binding. Using the neurite tip accumulation assay, we demonstrated that all novel KIF1A variants significantly reduced the ability of the motor domain of KIF1A to accumulate along the neurite lengths of differentiated SH-SY5Y cells. In vitro microtubule gliding assays showed significantly reduced velocities for the variant p.(Asp248Glu) and reduced microtubule binding for the p.(Cys92Arg) and p.(Pro305Leu) variants, suggesting a decreased ability of KIF1A to move along microtubules. Thus, this study further expanded the phenotypic characteristics of KAND individuals with pathogenic variants in the KIF1A motor domain to include clinical features commonly seen in RTT individuals.

Gene selection for the Australian Reproductive Genetic Carrier Screening Project ("Mackenzie's Mission")

Reproductive genetic carrier screening aims to offer couples information about their chance of having children with certain autosomal recessive and X-linked genetic conditions. We developed a gene list for use in "Mackenzie's Mission", a research project in which 10,000 couples will undergo screening. Criteria for selecting genes were: the condition should be life-limiting or disabling, with childhood onset, such that couples would be likely to take steps to avoid having an affected child; and/or be one for which early diagnosis and intervention would substantially change outcome. Strong evidence for gene-phenotype relationship was required. Candidate genes were identified from OMIM and via review of 23 commercial and published gene lists. Genes were reviewed by 16 clinical geneticists using a standard operating procedure, in a process overseen by a multidisciplinary committee which included clinical geneticists, genetic counselors, an ethicist, a parent of a child with a genetic condition and scientists from diagnostic and research backgrounds. 1300 genes met criteria. Genes associated with non-syndromic deafness and non-syndromic differences of sex development were not included. Our experience has highlighted that gene selection for a carrier screening panel needs to be a dynamic process with ongoing review and refinement.

Feasibility of Ultra-Rapid Exome Sequencing in Critically Ill Infants and Children With Suspected Monogenic Conditions in the Australian Public Health Care System

IMPORTANCE

Widespread adoption of rapid genomic testing in pediatric critical care requires robust clinical and laboratory pathways that provide equitable and consistent service across health care systems.

OBJECTIVE

To prospectively evaluate the performance of a multicenter network for ultra-rapid genomic diagnosis in a public health care system.

DESIGN, SETTING, AND PARTICIPANTS

Descriptive feasibility study of critically ill pediatric patients with suspected monogenic conditions treated at 12 Australian hospitals between March 2018 and February 2019, with data collected to May 2019. A formal implementation strategy emphasizing communication and feedback, standardized processes, coordination, distributed leadership, and collective learning was used to facilitate adoption.

EXPOSURES

Ultra-rapid exome sequencing.

MAIN OUTCOMES AND MEASURES

The primary outcome was time from sample receipt to ultra-rapid exome sequencing report. The secondary outcomes were the molecular diagnostic yield, the change in clinical management after the ultra-rapid exome sequencing report, the time from hospital admission to the laboratory report, and the proportion of laboratory reports returned prior to death or hospital discharge.

RESULTS

The study population included 108 patients with a median age of 28 days (range, 0 days to 17 years); 34% were female; and 57% were from neonatal intensive care units, 33% were from pediatric intensive care units, and 9% were from other hospital wards. The mean time from sample receipt to ultra-rapid exome sequencing report was 3.3 days (95% CI, 3.2-3.5 days) and the median time was 3 days (range, 2-7 days). The mean time from hospital admission to ultra-rapid exome sequencing report was 17.5 days (95% CI, 14.6-21.1 days) and 93 reports (86%) were issued prior to death or hospital discharge. A molecular diagnosis was established in 55 patients (51%). Eleven diagnoses (20%) resulted from using the following approaches to augment standard exome sequencing analysis: mitochondrial genome sequencing analysis, exome sequencing-based copy number analysis, use of international databases to identify novel gene-disease associations, and additional phenotyping and RNA analysis. In 42 of 55 patients (76%) with a molecular diagnosis and 6 of 53 patients (11%) without a molecular diagnosis, the ultra-rapid exome sequencing result was considered as having influenced clinical management. Targeted treatments were initiated in 12 patients (11%), treatment was redirected toward palliative care in 14 patients (13%), and surveillance for specific complications was initiated in 19 patients (18%).

CONCLUSIONS AND RELEVANCE

This study suggests feasibility of ultra-rapid genomic testing in critically ill pediatric patients with suspected monogenic conditions in the Australian public health care system. However, further research is needed to understand the clinical value of such testing, and the generalizability of the findings to other health care settings.

Biallelic Mutations in MTPAP Associated with a Lethal Encephalopathy

BACKGROUND

A homozygous founder mutation in MTPAP/TENT6, encoding mitochondrial poly(A) polymerase (MTPAP), was first reported in six individuals of Old Order Amish descent demonstrating an early-onset, progressive spastic ataxia with optic atrophy and learning difficulties. MTPAP contributes to the regulation of mitochondrial gene expression through the polyadenylation of mitochondrially encoded mRNAs. Mitochondrial mRNAs with severely truncated poly(A) tails were observed in affected individuals, and mitochondrial protein expression was altered.

OBJECTIVE

To determine the genetic basis of a perinatal encephalopathy associated with stereotyped neuroimaging and infantile death in three patients from two unrelated families.

METHODS

Whole-exome sequencing was performed in two unrelated patients and the unaffected parents of one of these individuals. Variants and familial segregation were confirmed by Sanger sequencing. Polyadenylation of mitochondrial transcripts and de novo synthesis of mitochondrial proteins were assessed in patient's fibroblasts.

RESULTS

Compound heterozygous p.Ile428Thr and p.Arg523Trp substitutions in MTPAP were recorded in two affected siblings from one family, and a homozygous p.Ile385Phe missense variant identified in a further affected child from a second sibship. Mitochondrial poly(A) tail analysis demonstrated shorter posttranscriptional additions to the mitochondrial transcripts, as well as an altered expression of mitochondrial proteins in the fibroblasts of the two siblings compared with healthy controls.

CONCLUSION

Mutations in MTPAP likely cause an autosomal recessive perinatal encephalopathy with lethality in the first year of life.

CDH1-related blepharocheilodontic syndrome is associated with diffuse gastric cancer risk

We report the first case of diffuse gastric cancer in an individual with familial blepharocheilodontic syndrome (BCD) due to a germline CDH1 likely pathogenic variant. To date, other BCD affected relatives are nonpenetrant for diffuse gastric cancer posing challenges to counseling regarding gastric and breast cancer surveillance, and preventative total gastrectomy.

Impacts of genomics on the health and social costs of intellectual disability

BACKGROUND

This study provides an integrated assessment of the economic and social impacts of genomic sequencing for the detection of monogenic disorders resulting in intellectual disability (ID).

METHODS

Multiple knowledge bases were cross-referenced and analysed to compile a reference list of monogenic disorders associated with ID. Multiple literature searches were used to quantify the health and social costs for the care of people with ID. Health and social expenditures and the current cost of whole-exome sequencing and whole-genome sequencing were quantified in relation to the more common causes of ID and their impact on lifespan.

RESULTS

On average, individuals with ID incur annual costs in terms of health costs, disability support, lost income and other social costs of US$172 000, accumulating to many millions of dollars over a lifetime.

CONCLUSION

The diagnosis of monogenic disorders through genomic testing provides the opportunity to improve the diagnosis and management, and to reduce the costs of ID through informed reproductive decisions, reductions in unproductive diagnostic tests and increasingly targeted therapies.

Missense variants in TAF1 and developmental phenotypes: challenges of determining pathogenicity

We recently described a new neurodevelopmental syndrome (TAF1/MRXS33 intellectual disability syndrome) (MIM# 300966) caused by pathogenic variants involving the X-linked gene TAF1, which participates in RNA polymerase II transcription. The initial study reported eleven families, and the syndrome was defined as presenting early in life with hypotonia, facial dysmorphia, and developmental delay that evolved into intellectual disability (ID) and/or autism spectrum disorder (ASD). We have now identified an additional 27 families through a genotype-first approach. Familial segregation analysis, clinical phenotyping, and bioinformatics were capitalized on to assess potential variant pathogenicity, and molecular modelling was performed for those variants falling within structurally characterized domains of TAF1. A novel phenotypic clustering approach was also applied, in which the phenotypes of affected individuals were classified using 51 standardized Human Phenotype Ontology (HPO) terms. Phenotypes associated with TAF1 variants show considerable pleiotropy and clinical variability, but prominent among previously unreported effects were brain morphological abnormalities, seizures, hearing loss, and heart malformations. Our allelic series broadens the phenotypic spectrum of TAF1/MRXS33 intellectual disability syndrome and the range of TAF1 molecular defects in humans. It also illustrates the challenges for determining the pathogenicity of inherited missense variants, particularly for genes mapping to chromosome X. This article is protected by copyright. All rights reserved.

Seave: a comprehensive web platform for storing and interrogating human genomic variation

Motivation

Genome sequencing has had a remarkable impact on our ability to study the effects of human genetic variation, however, variant interpretation remains the major bottleneck. Understanding the potential impact of variants, including structural variants, requires extensive annotation from disparate sources of knowledge, and in silico prediction algorithms.

Results

We introduce Seave, an intuitive web platform that enables all types of variants to be securely stored, annotated and filtered. Variants are annotated with allele frequencies and pathogenicity assessments from many popular databases and in silico pathogenicity prediction scores. Seave enables filtering of variants with specific inheritance patterns, including somatic variants, by quality, allele frequencies and gene lists which can be curated and saved. Seave was made for whole genome data and is capable of storing and querying copy number and structural variants.

Availability and implementation

To demo Seave with public data, see https://www.seave.bio. Source code is available at http://code.seave.bio and extensive documentation is available at http://documentation.seave.bio. Seave can be locally installed on an Apache server with PHP and MySQL, or we provide an Amazon Machine Image for quick deployment. For commercial and clinical diagnostic licensing, contact the corresponding author.

Supplementary information

Supplementary data are available at Bioinformatics online.

Fetal diagnosis of Mowat-Wilson syndrome by whole exome sequencing

Mowat-Wilson syndrome (MWS) is a complex genetic disorder associated with heterozygous variation in ZEB2. It is mainly characterized by moderate-to-severe intellectual disability, facial dysmorphism, epilepsy, and various malformations including Hirschsprung disease, corpus callosum anomalies, and congenital heart defects. It is rarely diagnosed prenatally and there is limited information available on the prenatal phenotype associated with MWS. Here we report the detection of a heterozygous de novo nonsense variant in ZEB2 by whole exome sequencing in a fetus with microphthalmia in addition to cardiac defects and typical MWS facial dysmorphism. As the prenatal phenotypic spectrum of MWS expands, the routine addition of fetal genomic testing particularly in the presence of multiple malformations will increase both the sensitivity and specificity of prenatal diagnostics.

Mutations in GDF11 and the extracellular antagonist, Follistatin, as a likely cause of Mendelian forms of orofacial clefting in humans

Cleft lip with or without cleft palate (CL/P) is generally viewed as a complex trait with multiple genetic and environmental contributions. In 70% of cases, CL/P presents as an isolated feature and/or deemed nonsyndromic. In the remaining 30%, CL/P is associated with multisystem phenotypes or clinically recognizable syndromes, many with a monogenic basis. Here we report the identification, via exome sequencing, of likely pathogenic variants in two genes that encode interacting proteins previously only linked to orofacial clefting in mouse models. A variant in GDF11 (encoding growth differentiation factor 11), predicting a p.(Arg298Gln) substitution at the Furin protease cleavage site, was identified in one family that segregated with CL/P and both rib and vertebral hypersegmentation, mirroring that seen in Gdf11 knockout mice. In the second family in which CL/P was the only phenotype, a mutation in FST (encoding the GDF11 antagonist, Follistatin) was identified that is predicted to result in a p.(Cys56Tyr) substitution in the region that binds GDF11. Functional assays demonstrated a significant impact of the specific mutated amino acids on FST and GDF11 function and, together with embryonic expression data, provide strong evidence for the importance of GDF11 and Follistatin in the regulation of human orofacial development.

Glibenclamide treatment in a Cantú syndrome patient with a pathogenic ABCC9 gain-of-function variant: Initial experience

Cantú syndrome (CS), characterized by hypertrichosis, distinctive facial features, and complex cardiovascular abnormalities, is caused by pathogenic variants in ABCC9 and KCNJ8 genes. These genes encode gain‐of‐function mutations in the regulatory (SUR2) and pore‐forming (Kir6.1) subunits of K_ATP channels, respectively, suggesting that channel‐blocking sulfonylureas could be a viable therapy. Here we report a neonate with CS, carrying a heterozygous ABCC9 variant (c.3347G>A, p.Arg1116His), born prematurely at 32 weeks gestation. Initial echocardiogram revealed a large patent ductus arteriosus (PDA), and high pulmonary pressures with enlarged right ventricle. He initially received surfactant and continuous positive airway pressure ventilation and was invasively ventilated for 4 weeks, until PDA ligation. After surgery, he still had ongoing bilevel positive airway pressure (BiPAP) requirement, but was subsequently weaned to nocturnal BiPAP. He was treated for pulmonary hypertension with Sildenafil, but failed to make further clinical improvement. A therapeutic glibenclamide trial was commenced in week 11 (initial dose of 0.05 mg^–1 kg^–1 day^–1 in two divided doses). After 1 week of treatment, he began to tolerate time off BiPAP when awake, and edema improved. Glibenclamide was well tolerated, and the dose was slowly increased to 0.15 mg⁻¹ kg⁻¹day⁻¹ over the next 12 weeks. Mild transient hypoglycemia was observed, but there was no cardiovascular dysfunction. Confirmation of therapeutic benefit will require studies of more CS patients but, based on this limited experience, consideration should be given to glibenclamide as CS therapy, although problems associated with prematurity, and complications of hypoglycemia, might limit outcome in critically ill neonates with CS.

Pathogenic variants in PLOD3 result in a Stickler syndrome-like connective tissue disorder with vascular complications

Background

Pathogenic PLOD3 variants cause a connective tissue disorder (CTD) that has been described rarely. We further characterise this CTD and propose a clinical diagnostic label to improve recognition and diagnosis of PLOD3-related disease.

Methods

Reported PLOD3 phenotypes were compared with known CTDs utilising data from three further individuals from a consanguineous family with a homozygous PLOD3 c.809C>T; p.(Pro270Leu) variant. PLOD3 mRNA expression in the developing embryo was analysed for tissue-specific localisation. Mouse microarray expression data were assessed for phylogenetic gene expression similarities across CTDs with overlapping clinical features.

Results

Key clinical features included ocular abnormalities with risk for retinal detachment, sensorineural hearing loss, reduced palmar creases, finger contractures, prominent knees, scoliosis, low bone mineral density, recognisable craniofacial dysmorphisms, developmental delay and risk for vascular dissection. Collated clinical features showed most overlap with Stickler syndrome with variable features of Ehlers-Danlos syndrome (EDS) and epidermolysis bullosa (EB). Human lysyl hydroxylase 3/PLOD3 expression was localised to the developing cochlea, eyes, skin, forelimbs, heart and cartilage, mirroring the clinical phenotype of this disorder.

Conclusion

These data are consistent with pathogenic variants in PLOD3 resulting in a clinically distinct Stickler-like syndrome with vascular complications and variable features of EDS and EB. Early identification of PLOD3 variants would improve monitoring for comorbidities and may avoid serious adverse ocular and vascular outcomes.