Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease

BACKGROUND

X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique.

RESULTS

In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis.

CONCLUSIONS

This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Cre/lox regulated conditional rescue and inactivation with zebrafish UFlip alleles generated by CRISPR-Cas9 targeted integration

The ability to regulate gene activity spatially and temporally is essential to investigate cell type-specific gene function during development and in postembryonic processes and disease models. The Cre/lox system has been widely used for performing cell and tissue-specific conditional analysis of gene function in zebrafish. However, simple and efficient methods for isolation of stable, Cre/lox regulated zebrafish alleles are lacking. Here we applied our GeneWeld CRISPR-Cas9 targeted integration strategy to generate floxed alleles that provide robust conditional inactivation and rescue. A universal targeting vector, UFlip, with sites for cloning short homology arms flanking a floxed 2A-mRFP gene trap, was integrated into an intron in rbbp4 and rb1. rbbp4^off and rb1^off integration alleles resulted in strong mRFP expression, >99% reduction of endogenous gene expression, and recapitulated known indel loss of function phenotypes. Introduction of Cre led to stable inversion of the floxed cassette, loss of mRFP expression, and phenotypic rescue. rbbp4^on and rb1^on integration alleles did not cause phenotypes in combination with a loss of function mutation. Addition of Cre led to conditional inactivation by stable inversion of the cassette, gene trapping and mRFP expression, and the expected mutant phenotype. Neural progenitor Cre drivers were used for conditional inactivation and phenotypic rescue to showcase how this approach can be used in specific cell populations. Together these results validate a simplified approach for efficient isolation of Cre/lox responsive conditional alleles in zebrafish. Our strategy provides a new toolkit for generating genetic mosaics and represents a significant advance in zebrafish genetics.

Pathogenic SPTBN1 variants cause an autosomal dominant neurodevelopmental syndrome

SPTBN1 encodes βII-spectrin, the ubiquitously expressed β-spectrin that forms micrometer-scale networks associated with plasma membranes. Mice deficient in neuronal βII-spectrin have defects in cortical organization, developmental delay and behavioral deficiencies. These phenotypes, while less severe, are observed in haploinsufficient animals, suggesting that individuals carrying heterozygous SPTBN1 variants may also show measurable compromise of neural development and function. Here we identify heterozygous SPTBN1 variants in 29 individuals with developmental, language and motor delays; mild to severe intellectual disability; autistic features; seizures; behavioral and movement abnormalities; hypotonia; and variable dysmorphic facial features. We show that these SPTBN1 variants lead to effects that affect βII-spectrin stability, disrupt binding to key molecular partners, and disturb cytoskeleton organization and dynamics. Our studies define SPTBN1 variants as the genetic basis of a neurodevelopmental syndrome, expand the set of spectrinopathies affecting the brain and underscore the critical role of βII-spectrin in the central nervous system.

CSNK2B: A broad spectrum of neurodevelopmental disability and epilepsy severity

CSNK2B has recently been implicated as a disease gene for neurodevelopmental disability (NDD) and epilepsy. Information about developmental outcomes has been limited by the young age and short follow-up for many of the previously reported cases, and further delineation of the spectrum of associated phenotypes is needed. We present 25 new patients with variants in CSNK2B and refine the associated NDD and epilepsy phenotypes. CSNK2B variants were identified by research or clinical exome sequencing, and investigators from different centers were connected via GeneMatcher. Most individuals had developmental delay and generalized epilepsy with onset in the first 2 years. However, we found a broad spectrum of phenotypic severity, ranging from early normal development with pharmacoresponsive seizures to profound intellectual disability with intractable epilepsy and recurrent refractory status epilepticus. These findings suggest that CSNK2B should be considered in the diagnostic evaluation of patients with a broad range of NDD with treatable or intractable seizures.

SPEN haploinsufficiency causes a neurodevelopmental disorder overlapping proximal 1p36 deletion syndrome with an episignature of X chromosomes in females

Deletion 1p36 (del1p36) syndrome is the most common human disorder resulting from a terminal autosomal deletion. This condition is molecularly and clinically heterogeneous. Deletions involving two non-overlapping regions, known as the distal (telomeric) and proximal (centromeric) critical regions, are sufficient to cause the majority of the recurrent clinical features, although with different facial features and dysmorphisms. SPEN encodes a transcriptional repressor commonly deleted in proximal del1p36 syndrome and is located centromeric to the proximal 1p36 critical region. Here, we used clinical data from 34 individuals with truncating variants in SPEN to define a neurodevelopmental disorder presenting with features that overlap considerably with those of proximal del1p36 syndrome. The clinical profile of this disease includes developmental delay/intellectual disability, autism spectrum disorder, anxiety, aggressive behavior, attention deficit disorder, hypotonia, brain and spine anomalies, congenital heart defects, high/narrow palate, facial dysmorphisms, and obesity/increased BMI, especially in females. SPEN also emerges as a relevant gene for del1p36 syndrome by co-expression analyses. Finally, we show that haploinsufficiency of SPEN is associated with a distinctive DNA methylation episignature of the X chromosome in affected females, providing further evidence of a specific contribution of the protein to the epigenetic control of this chromosome, and a paradigm of an X chromosome-specific episignature that classifies syndromic traits. We conclude that SPEN is required for multiple developmental processes and SPEN haploinsufficiency is a major contributor to a disorder associated with deletions centromeric to the previously established 1p36 critical regions.

Impact of integrated translational research on clinical exome sequencing

PURPOSE

Exome sequencing often identifies pathogenic genetic variants in patients with undiagnosed diseases. Nevertheless, frequent findings of variants of uncertain significance necessitate additional efforts to establish causality before reaching a conclusive diagnosis. To provide comprehensive genomic testing to patients with undiagnosed disease, we established an Individualized Medicine Clinic, which offered clinical exome testing and included a Translational Omics Program (TOP) that provided variant curation, research activities, or research exome sequencing.

METHODS

From 2012 to 2018, 1101 unselected patients with undiagnosed diseases received exome testing. Outcomes were reviewed to assess impact of the TOP and patient characteristics on diagnostic rates through descriptive and multivariate analyses.

RESULTS

The overall diagnostic yield was 24.9% (274 of 1101 patients), with 174 (15.8% of 1101) diagnosed on the basis of clinical exome sequencing alone. Four hundred twenty-three patients with nondiagnostic or without access to clinical exome sequencing were evaluated by the TOP, with 100 (9% of 1101) patients receiving a diagnosis, accounting for 36.5% of the diagnostic yield. The identification of a genetic diagnosis was influenced by the age at time of testing and the disease phenotype of the patient.

CONCLUSION

Integration of translational research activities into clinical practice of a tertiary medical center can significantly increase the diagnostic yield of patients with undiagnosed disease.

SPECC1L regulates palate development downstream of IRF6

SPECC1L mutations have been identified in patients with rare atypical orofacial clefts and with syndromic cleft lip and/or palate (CL/P). These mutations cluster in the second coiled-coil and calponin homology domains of SPECC1L and severely affect the ability of SPECC1L to associate with microtubules. We previously showed that gene-trap knockout of Specc1l in mouse results in early embryonic lethality. We now present a truncation mutant mouse allele, Specc1lΔC510, that results in perinatal lethality. Specc1lΔC510/ΔC510 homozygotes showed abnormal palate rugae but did not show cleft palate. However, when crossed with a gene-trap allele, Specc1lcGT/ΔC510 compound heterozygotes showed a palate elevation delay with incompletely penetrant cleft palate. Specc1lcGT/ΔC510 embryos exhibit transient oral epithelial adhesions at E13.5, which may delay shelf elevation. Consistent with oral adhesions, we show periderm layer abnormalities, including ectopic apical expression of adherens junction markers, similar to Irf6 hypomorphic mutants and Arhgap29 heterozygotes. Indeed, SPECC1L expression is drastically reduced in Irf6 mutant palatal shelves. Finally, we wanted to determine if SPECC1L deficiency also contributed to non-syndromic (ns) CL/P. We sequenced 62 Caucasian, 89 Filipino, 90 Ethiopian, 90 Nigerian and 95 Japanese patients with nsCL/P and identified three rare coding variants (p.Ala86Thr, p.Met91Iso and p.Arg546Gln) in six individuals. These variants reside outside of SPECC1L coiled-coil domains and result in milder functional defects than variants associated with syndromic clefting. Together, our data indicate that palate elevation is sensitive to deficiency of SPECC1L dosage and function and that SPECC1L cytoskeletal protein functions downstream of IRF6 in palatogenesis.

Loss-of-function mutations in UDP-Glucose 6-Dehydrogenase cause recessive developmental epileptic encephalopathy

Developmental epileptic encephalopathies are devastating disorders characterized by intractable epileptic seizures and developmental delay. Here, we report an allelic series of germline recessive mutations in UGDH in 36 cases from 25 families presenting with epileptic encephalopathy with developmental delay and hypotonia. UGDH encodes an oxidoreductase that converts UDP-glucose to UDP-glucuronic acid, a key component of specific proteoglycans and glycolipids. Consistent with being loss-of-function alleles, we show using patients' primary fibroblasts and biochemical assays, that these mutations either impair UGDH stability, oligomerization, or enzymatic activity. In vitro, patient-derived cerebral organoids are smaller with a reduced number of proliferating neuronal progenitors while mutant ugdh zebrafish do not phenocopy the human disease. Our study defines UGDH as a key player for the production of extracellular matrix components that are essential for human brain development. Based on the incidence of variants observed, UGDH mutations are likely to be a frequent cause of recessive epileptic encephalopathy.

A tailored approach to fusion transcript identification increases diagnosis of rare inherited disease

BACKGROUND

RNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data.

METHODS

We describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization.

RESULTS

We demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance.

CONCLUSIONS

The approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.