The clinical utility and diagnostic implementation of human subject cell transdifferentiation followed by RNA sequencing

RNA sequencing (RNA-seq) has recently been used in translational research settings to facilitate diagnoses of Mendelian disorders. A significant obstacle for clinical laboratories in adopting RNA-seq is the low or absent expression of a significant number of disease-associated genes/transcripts in clinically accessible samples. As this is especially problematic in neurological diseases, we developed a clinical diagnostic approach that enhanced the detection and evaluation of tissue-specific genes/transcripts through fibroblast-to-neuron cell transdifferentiation. The approach is designed specifically to suit clinical implementation, emphasizing simplicity, cost effectiveness, turnaround time, and reproducibility. For clinical validation, we generated induced neurons (iNeurons) from 71 individuals with primary neurological phenotypes recruited to the Undiagnosed Diseases Network. The overall diagnostic yield was 25.4%. Over a quarter of the diagnostic findings benefited from transdifferentiation and could not be achieved by fibroblast RNA-seq alone. This iNeuron transcriptomic approach can be effectively integrated into diagnostic whole-transcriptome evaluation of individuals with genetic disorders.

SPTSSA variants alter sphingolipid synthesis and cause a complex hereditary spastic paraplegia

Sphingolipids are a diverse family of lipids with critical structural and signalling functions in the mammalian nervous system, where they are abundant in myelin membranes. Serine palmitoyltransferase, the enzyme that catalyses the rate-limiting reaction of sphingolipid synthesis, is composed of multiple subunits including an activating subunit, SPTSSA. Sphingolipids are both essential and cytotoxic and their synthesis must therefore be tightly regulated. Key to the homeostatic regulation are the ORMDL proteins that are bound to serine palmitoyltransferase and mediate feedback inhibition of enzymatic activity when sphingolipid levels become excessive. Exome sequencing identified potential disease-causing variants in SPTSSA in three children presenting with a complex form of hereditary spastic paraplegia. The effect of these variants on the catalytic activity and homeostatic regulation of serine palmitoyltransferase was investigated in human embryonic kidney cells, patient fibroblasts and Drosophila. Our results showed that two different pathogenic variants in SPTSSA caused a hereditary spastic paraplegia resulting in progressive motor disturbance with variable sensorineural hearing loss and language/cognitive dysfunction in three individuals. The variants in SPTSSA impaired the negative regulation of serine palmitoyltransferase by ORMDLs leading to excessive sphingolipid synthesis based on biochemical studies and in vivo studies in Drosophila. These findings support the pathogenicity of the SPTSSA variants and point to excessive sphingolipid synthesis due to impaired homeostatic regulation of serine palmitoyltransferase as responsible for defects in early brain development and function.

Bi-allelic ATG4D variants are associated with a neurodevelopmental disorder characterized by speech and motor impairment

Autophagy regulates the degradation of damaged organelles and protein aggregates, and is critical for neuronal development, homeostasis, and maintenance, yet few neurodevelopmental disorders have been associated with pathogenic variants in genes encoding autophagy-related proteins. We report three individuals from two unrelated families with a neurodevelopmental disorder characterized by speech and motor impairment, and similar facial characteristics. Rare, conserved, bi-allelic variants were identified in ATG4D, encoding one of four ATG4 cysteine proteases important for autophagosome biogenesis, a hallmark of autophagy. Autophagosome biogenesis and induction of autophagy were intact in cells from affected individuals. However, studies evaluating the predominant substrate of ATG4D, GABARAPL1, demonstrated that three of the four ATG4D patient variants functionally impair ATG4D activity. GABARAPL1 is cleaved or "primed" by ATG4D and an in vitro GABARAPL1 priming assay revealed decreased priming activity for three of the four ATG4D variants. Furthermore, a rescue experiment performed in an ATG4 tetra knockout cell line, in which all four ATG4 isoforms were knocked out by gene editing, showed decreased GABARAPL1 priming activity for the two ATG4D missense variants located in the cysteine protease domain required for priming, suggesting that these variants impair the function of ATG4D. The clinical, bioinformatic, and functional data suggest that bi-allelic loss-of-function variants in ATG4D contribute to the pathogenesis of this syndromic neurodevelopmental disorder.

Continuing a search for a diagnosis: the impact of adolescence and family dynamics

The "diagnostic odyssey" describes the process those with undiagnosed conditions undergo to identify a diagnosis. Throughout this process, families of children with undiagnosed conditions have multiple opportunities to decide whether to continue or stop their search for a diagnosis and accept the lack of a diagnostic label. Previous studies identified factors motivating a family to begin searching, but there is limited information about the decision-making process in a prolonged search and how the affected child impacts a family's decision. This study aimed to understand how families of children with undiagnosed diseases decide whether to continue to pursue a diagnosis after standard clinical testing has failed. Parents who applied to the Undiagnosed Disease Network (UDN) at the National Institutes of Health (NIH) were recruited to participate in semi-structured interviews. The 2015 Supportive Care Needs model by Pelenstov, which defines critical needs in families with rare/undiagnosed diseases, provided a framework for interview guide development and transcript analysis (Pelentsov et al in Disabil Health J 8(4):475-491, 2015. https://doi.org/10.1016/J.DHJO.2015.03.009 ). A deductive, iterative coding approach was used to identify common unifying themes. Fourteen parents from 13 families were interviewed. The average child's age was 11 years (range 3-18) and an average 63% of their life had been spent searching for a diagnosis. Our analysis found that alignment or misalignment of parent and child needs impact the trajectory of the diagnostic search. When needs and desires align, reevaluation of a decision to pursue a diagnosis is limited. However, when there is conflict between parent and child desires, there is reevaluation, and often a pause, in the search. This tension is exacerbated when children are adolescents and attempting to balance their dependence on parents for medical care with a natural desire for independence. Our results provide novel insights into the roles of adolescents in the diagnostic odyssey. The tension between desired and realistic developmental outcomes for parents and adolescents impacts if, and how, the search for a diagnosis progresses.

Retrospective analysis of a clinical exome sequencing cohort reveals the mutational spectrum and identifies candidate disease-associated loci for BAFopathies

PURPOSE

BRG1/BRM-associated factor (BAF) complex is a chromatin remodeling complex that plays a critical role in gene regulation. Defects in the genes encoding BAF subunits lead to BAFopathies, a group of neurodevelopmental disorders with extensive locus and phenotypic heterogeneity.

METHODS

We retrospectively analyzed data from 16,243 patients referred for clinical exome sequencing (ES) with a focus on the BAF complex. We applied a genotype-first approach, combining predicted genic constraints to propose candidate BAFopathy genes.

RESULTS

We identified 127 patients carrying pathogenic variants, likely pathogenic variants, or de novo variants of unknown clinical significance in 11 known BAFopathy genes. Those include 34 patients molecularly diagnosed using ES reanalysis with new gene-disease evidence (n = 21) or variant reclassifications in known BAFopathy genes (n = 13). We also identified de novo or predicted loss-of-function variants in 4 candidate BAFopathy genes, including ACTL6A, BICRA (implicated in Coffin-Siris syndrome during this study), PBRM1, and SMARCC1.

CONCLUSION

We report the mutational spectrum of BAFopathies in an ES cohort. A genotype-driven and pathway-based reanalysis of ES data identified new evidence for candidate genes involved in BAFopathies. Further mechanistic and phenotypic characterization of additional patients are warranted to confirm their roles in human disease and to delineate their associated phenotypic spectrums.

De novo variants in H3-3A and H3-3B are associated with neurodevelopmental delay, dysmorphic features, and structural brain abnormalities

The histone H3 variant H3.3, encoded by two genes H3-3A and H3-3B, can replace canonical isoforms H3.1 and H3.2. H3.3 is important in chromatin compaction, early embryonic development, and lineage commitment. The role of H3.3 in somatic cancers has been studied extensively, but its association with a congenital disorder has emerged just recently. Here we report eleven de novo missense variants and one de novo stop-loss variant in H3-3A (n = 6) and H3-3B (n = 6) from Baylor Genetics exome cohort (n = 11) and Matchmaker Exchange (n = 1), of which detailed phenotyping was conducted for 10 individuals (H3-3A = 4 and H3-3B = 6) that showed major phenotypes including global developmental delay, short stature, failure to thrive, dysmorphic facial features, structural brain abnormalities, hypotonia, and visual impairment. Three variant constructs (p.R129H, p.M121I, and p.I52N) showed significant decrease in protein expression, while one variant (p.R41C) accumulated at greater levels than wild-type control. One H3.3 variant construct (p.R129H) was found to have stronger interaction with the chaperone death domain-associated protein 6.

Contribution of uniparental disomy in a clinical trio exome cohort of 2675 patients

Background

Uniparental disomy (UPD) is the inheritance of two homologous chromosomes from the same parent. UPD may result in clinical phenotypes when occurring on chromosomes with specific imprinting pattern, when leading to homozygosity of a deleterious recessive allele inherited from one carrier parent, or when associated with a mosaic aneuploidy. Due to the importance of UPD in genetic disease etiology, UPD analysis has started to be implemented in the context of exome sequencing (ES) or genome sequencing.

Methods

We developed an in‐house algorithm TRIPS (Trio Parentage/UPD Studies) to identify UPD events in trio ES cases. This method identifies regions with uniparental inheritance by utilizing the trio genotyping data obtained from the concurrent SNP array to delineate the parental origin of the SNPs in the proband.

Results

We identified 16 UPD events from 2675 ES trios. Among those, four events led to imprinting disorders, seven unmasked a pathogenic/likely pathogenic variant in a recessive disease gene, and two were consistent with a mosaic genome wide paternal UPD pattern. Twelve of these UPD events directly contributed to the molecular diagnosis of the patients.

Conclusion

Our study demonstrated the contribution of UPD to the molecular diagnosis in one clinical ES cohort, thus UPD analysis should be incorporated into routine clinical ES interpretation.

PPP3CA truncating variants clustered in the regulatory domain cause early-onset refractory epilepsy

PPP3CA encodes the catalytic subunit of calcineurin, a calcium-calmodulin-regulated serine-threonine phosphatase. Loss-of-function (LoF) variants in the catalytic domain have been associated with epilepsy, while gain-of-function (GoF) variants in the auto-inhibitory domain cause multiple congenital abnormalities. We herein report five new patients with de novo PPP3CA variants. Interestingly, the two frameshift variants in this study and the six truncating variants reported previously are all located within a 26-amino acid region in the regulatory domain (RD). Patients with a truncating variant had more severe earlier onset seizures compared to patients with a LoF missense variant, while autism spectrum disorder was a more frequent feature in the latter. Expression studies of a truncating variant showed apparent RNA expression from the mutant allele, but no detectable mutant protein. Our data suggest that PPP3CA truncating variants clustered in the RD, causing more severe early-onset refractory epilepsy and representing a type of variants distinct from LoF or GoF missense variants.

Heterozygous variants in SPTBN1 cause intellectual disability and autism

Spectrins are common components of cytoskeletons, binding to cytoskeletal elements and the plasma membrane, allowing proper localization of essential membrane proteins, signal transduction, and cellular scaffolding. Spectrins are assembled from α and β subunits, encoded by SPTA1 and SPTAN1 (α) and SPTB, SPTBN1, SPTBN2, SPTBN4, and SPTBN5 (β). Pathogenic variants in various spectrin genes are associated with erythroid cell disorders (SPTA1, SPTB) and neurologic disorders (SPTAN1, SPTBN2, and SPTBN4), but no phenotypes have been definitively associated with variants in SPTBN1 or SPTBN5. Through exome sequencing and case matching, we identified seven unrelated individuals with heterozygous SPTBN1 variants: two with de novo missense variants and five with predicted loss-of-function variants (found to be de novo in two, while one was inherited from a mother with a history of learning disabilities). Common features include global developmental delays, intellectual disability, and behavioral disturbances. Autistic features (4/6) and epilepsy (2/7) or abnormal electroencephalogram without overt seizures (1/7) were present in a subset. Identification of loss-of-function variants suggests a haploinsufficiency mechanism, but additional functional studies are required to fully elucidate disease pathogenesis. Our findings support the essential roles of SPTBN1 in human neurodevelopment and expand the knowledge of human spectrinopathy disorders.

Compound heterozygous KCTD7 variants in progressive myoclonus epilepsy

KCTD7 is a member of the potassium channel tetramerization domain-containing protein family and has been associated with progressive myoclonic epilepsy (PME), characterized by myoclonus, epilepsy, and neurological deterioration. Here we report four affected individuals from two unrelated families in which we identified KCTD7 compound heterozygous single nucleotide variants through exome sequencing. RNAseq was used to detect a non-annotated splicing junction created by a synonymous variant in the second family. Whole-cell patch-clamp analysis of neuroblastoma cells overexpressing the patients' variant alleles demonstrated aberrant potassium regulation. While all four patients experienced many of the common clinical features of PME, they also showed variable phenotypes not previously reported, including dysautonomia, brain pathology findings including a significantly reduced thalamus, and the lack of myoclonic seizures. To gain further insight into the pathogenesis of the disorder, zinc finger nucleases were used to generate kctd7 knockout zebrafish. Kctd7 homozygous mutants showed global dysregulation of gene expression and increased transcription of c-fos, which has previously been correlated with seizure activity in animal models. Together these findings expand the known phenotypic spectrum of KCTD7-associated PME, report a new animal model for future studies, and contribute valuable insights into the disease.

Low-level parental somatic mosaic SNVs in exomes from a large cohort of trios with diverse suspected Mendelian conditions

Purpose:

The goal of this study was to assess the scale of low-level parental mosaicism in exome sequencing (ES) databases.

Methods:

We analyzed approximately 2000 family trio ES datasets from the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) and Baylor Genetics (BG). Among apparent de novo single nucleotide variants (SNVs) identified in the affected probands, we selected rare unique variants with variant allele fraction (VAF) between 30-70% in the probands and lower than 10% in one of the parents.

Results:

Out of 102 candidate mosaic variants validated using amplicon-based NGS, droplet digital PCR, or blocker displacement amplification, 27 (26.4%) were confirmed to be low- (VAF between 1-10%) or very low- (VAF <1%) level mosaic. Detection precision in parental samples with two or more alternate reads was 63.6% (BHCMG) and 43.6% (BG). In nine investigated individuals, we observed variability of mosaic ratios among blood, saliva, fibroblast, buccal, hair, and urine samples.

Conclusion:

Our computational pipeline enables robust discrimination between true and false positive candidate mosaic variants and efficient detection of low-level mosaicism in ES samples. We confirm that the presence of two or more alternate reads in the parental sample is a reliable predictor of low-level parental somatic mosaicism.

CNVs cause autosomal recessive genetic diseases with or without involvement of SNV/indels

PURPOSE

Improved resolution of molecular diagnostic technologies enabled detection of smaller sized exonic level copy-number variants (CNVs). The contribution of CNVs to autosomal recessive (AR) conditions may be better recognized using a large clinical cohort.

METHODS

We retrospectively investigated the CNVs' contribution to AR conditions in cases subjected to chromosomal microarray analysis (CMA, N = ~70,000) and/or clinical exome sequencing (ES, N = ~12,000) at Baylor Genetics; most had pediatric onset neurodevelopmental disorders.

RESULTS

CNVs contributed to biallelic variations in 87 cases, including 81 singletons and three affected sibling pairs. Seventy cases had CNVs affecting both alleles, and 17 had a CNV and a single-nucleotide variant (SNV)/indel in trans. In total, 94.3% of AR-CNVs affected one gene; among these 41.4% were single-exon and 35.0% were multiexon partial-gene events. Sixty-nine percent of homozygous AR-CNVs were embedded in homozygous genomic intervals. Five cases had large deletions unmasking an SNV/indel on the intact allele for a recessive condition, resulting in multiple molecular diagnoses.

CONCLUSIONS

AR-CNVs are often smaller in size, transmitted through generations, and underrecognized due to limitations in clinical CNV detection methods. Our findings from a large clinical cohort emphasized integrated CNV and SNV/indel analyses for precise clinical and molecular diagnosis especially in the context of genomic disorders.

De novo EIF2AK1 and EIF2AK2 Variants Are Associated with Developmental Delay, Leukoencephalopathy, and Neurologic Decompensation

EIF2AK1 and EIF2AK2 encode members of the eukaryotic translation initiation factor 2 alpha kinase (EIF2AK) family that inhibits protein synthesis in response to physiologic stress conditions. EIF2AK2 is also involved in innate immune response and the regulation of signal transduction, apoptosis, cell proliferation, and differentiation. Despite these findings, human disorders associated with deleterious variants in EIF2AK1 and EIF2AK2 have not been reported. Here, we describe the identification of nine unrelated individuals with heterozygous de novo missense variants in EIF2AK1 (1/9) or EIF2AK2 (8/9). Features seen in these nine individuals include white matter alterations (9/9), developmental delay (9/9), impaired language (9/9), cognitive impairment (8/9), ataxia (6/9), dysarthria in probands with verbal ability (6/9), hypotonia (7/9), hypertonia (6/9), and involuntary movements (3/9). Individuals with EIF2AK2 variants also exhibit neurological regression in the setting of febrile illness or infection. We use mammalian cell lines and proband-derived fibroblasts to further confirm the pathogenicity of variants in these genes and found reduced kinase activity. EIF2AKs phosphorylate eukaryotic translation initiation factor 2 subunit 1 (EIF2S1, also known as EIF2α), which then inhibits EIF2B activity. Deleterious variants in genes encoding EIF2B proteins cause childhood ataxia with central nervous system hypomyelination/vanishing white matter (CACH/VWM), a leukodystrophy characterized by neurologic regression in the setting of febrile illness and other stressors. Our findings indicate that EIF2AK2 missense variants cause a neurodevelopmental syndrome that may share phenotypic and pathogenic mechanisms with CACH/VWM.

A Genocentric Approach to Discovery of Mendelian Disorders

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

Harmonizing Clinical Sequencing and Interpretation for the eMERGE III Network

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

A clinical survey of mosaic single nucleotide variants in disease-causing genes detected by exome sequencing

Background

Although mosaic variation has been known to cause disease for decades, high-throughput sequencing technologies with the analytical sensitivity to consistently detect variants at reduced allelic fractions have only recently emerged as routine clinical diagnostic tests. To date, few systematic analyses of mosaic variants detected by diagnostic exome sequencing for diverse clinical indications have been performed.

Methods

To investigate the frequency, type, allelic fraction, and phenotypic consequences of clinically relevant somatic mosaic single nucleotide variants (SNVs) and characteristics of the corresponding genes, we retrospectively queried reported mosaic variants from a cohort of ~ 12,000 samples submitted for clinical exome sequencing (ES) at Baylor Genetics.

Results

We found 120 mosaic variants involving 107 genes, including 80 mosaic SNVs in proband samples and 40 in parental/grandparental samples. Average mosaic alternate allele fraction (AAF) detected in autosomes and in X-linked disease genes in females was 18.2% compared with 34.8% in X-linked disease genes in males. Of these mosaic variants, 74 variants (61.7%) were classified as pathogenic or likely pathogenic and 46 (38.3%) as variants of uncertain significance. Mosaic variants occurred in disease genes associated with autosomal dominant (AD) or AD/autosomal recessive (AR) (67/120, 55.8%), X-linked (33/120, 27.5%), AD/somatic (10/120, 8.3%), and AR (8/120, 6.7%) inheritance. Of note, 1.7% (2/120) of variants were found in genes in which only somatic events have been described. Nine genes had recurrent mosaic events in unrelated individuals which accounted for 18.3% (22/120) of all detected mosaic variants in this study. The proband group was enriched for mosaicism affecting Ras signaling pathway genes.

Conclusions

In sum, an estimated 1.5% of all molecular diagnoses made in this cohort could be attributed to a mosaic variant detected in the proband, while parental mosaicism was identified in 0.3% of families analyzed. As ES design favors breadth over depth of coverage, this estimate of the prevalence of mosaic variants likely represents an underestimate of the total number of clinically relevant mosaic variants in our cohort.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0658-2) contains supplementary material, which is available to authorized users.

De novo missense variant in the GTPase effector domain (GED) of DNM1L leads to static encephalopathy and seizures

DNM1L encodes a GTPase of the dynamin superfamily, which plays a crucial role in mitochondrial and peroxisomal fission. Pathogenic variants affecting the middle domain and the GTPase domain of DNM1L have been implicated in encephalopathy because of defective mitochondrial and peroxisomal fission 1 (EMPF1, MIM #614388). Patients show variable phenotypes ranging from severe hypotonia leading to death in the neonatal period to developmental delay/regression, with or without seizures. Familial pathogenic variants in the GTPase domain have also been associated with isolated optic atrophy. We present a 27-yr-old woman with static encephalopathy, a history of seizures, and nystagmus, in whom a novel de novo heterozygous variant was detected in the GTPase effector domain (GED) of DNM1L (c.2072A>G, p.Tyr691Cys). Functional studies in Drosophila demonstrate large, abnormally distributed peroxisomes and mitochondria, an effect very similar to that of middle domain missense alleles observed in pediatric subjects with EMPF1. To our knowledge, not only is this the first report of a disease-causing variant in the GED domain in humans, but this is also the oldest living individual reported with EMPF1. Longitudinal data of this kind helps to expand our knowledge of the natural history of a growing list of DNM1L-related disorders.

Copy number variant and runs of homozygosity detection by microarrays enabled more precise molecular diagnoses in 11,020 clinical exome cases

Background

Exome sequencing (ES) has been successfully applied in clinical detection of single nucleotide variants (SNVs) and small indels. However, identification of copy number variants (CNVs) using ES data remains challenging. The purpose of this study is to understand the contribution of CNVs and copy neutral runs of homozygosity (ROH) in molecular diagnosis of patients referred for ES.

Methods

In a cohort of 11,020 consecutive ES patients, an Illumina SNP array analysis interrogating mostly coding SNPs was performed as a quality control (QC) measurement and for CNV/ROH detection. Among these patients, clinical chromosomal microarray analysis (CMA) was performed at Baylor Genetics (BG) on 3229 patients, either before, concurrently, or after ES. We retrospectively analyzed the findings from CMA and the QC array.

Results

The QC array can detect ~ 70% of pathogenic/likely pathogenic CNVs (PCNVs) detectable by CMA. Out of the 11,020 ES cases, the QC array identified PCNVs in 327 patients and uniparental disomy (UPD) disorder-related ROH in 10 patients. The overall PCNV/UPD detection rate was 5.9% in the 3229 ES patients who also had CMA at BG; PCNV/UPD detection rate was higher in concurrent ES and CMA than in ES with prior CMA (7.2% vs 4.6%). The PCNVs/UPD contributed to the molecular diagnoses in 17.4% (189/1089) of molecularly diagnosed ES cases with CMA and were estimated to contribute in 10.6% of all molecularly diagnosed ES cases. Dual diagnoses with both PCNVs and SNVs were detected in 38 patients. PCNVs affecting single recessive disorder genes in a compound heterozygous state with SNVs were detected in 4 patients, and homozygous deletions (mostly exonic deletions) were detected in 17 patients. A higher PCNV detection rate was observed for patients with syndromic phenotypes and/or cardiovascular abnormalities.

Conclusions

Our clinical genomics study demonstrates that detection of PCNV/UPD through the QC array or CMA increases ES diagnostic rate, provides more precise molecular diagnosis for dominant as well as recessive traits, and enables more complete genetic diagnoses in patients with dual or multiple molecular diagnoses. Concurrent ES and CMA using an array with exonic coverage for disease genes enables most effective detection of both CNVs and SNVs and therefore is recommended especially in time-sensitive clinical situations.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0639-5) contains supplementary material, which is available to authorized users.

A toolkit for genetics providers in follow-up of patients with non-diagnostic exome sequencing

There are approximately 7,000 rare diseases affecting 25-30 million Americans, with 80% estimated to have a genetic basis. This presents a challenge for genetics practitioners to determine appropriate testing, make accurate diagnoses, and conduct up-to-date patient management. Exome sequencing (ES) is a comprehensive diagnostic approach, but only 25%-41% of the patients receive a molecular diagnosis. The remaining three-fifths to three-quarters of patients undergoing ES remain undiagnosed. The Stanford Center for Undiagnosed Diseases (CUD), a clinical site of the Undiagnosed Diseases Network, evaluates patients with undiagnosed and rare diseases using a combination of methods including ES. Frequently these patients have non-diagnostic ES results, but strategic follow-up techniques identify diagnoses in a subset. We present techniques used at the CUD that can be adopted by genetics providers in clinical follow-up of cases where ES is non-diagnostic. Solved case examples illustrate different types of non-diagnostic results and the additional techniques that led to a diagnosis. Frequent approaches include segregation analysis, data reanalysis, genome sequencing, additional variant identification, careful phenotype-disease correlation, confirmatory testing, and case matching. We also discuss prioritization of cases for additional analyses.

Extracutaneous manifestations in phacomatosis cesioflammea and cesiomarmorata: Case series and literature review

Phacomatosis pigmentovascularis (PPV) comprises a family of rare conditions that feature vascular abnormalities and melanocytic lesions that can be solely cutaneous or multisystem in nature. Recently published work has demonstrated that both vascular and melanocytic abnormalities in PPV of the cesioflammea and cesiomarmorata subtypes can result from identical somatic mosaic activating mutations in the genes GNAQ and GNA11. Here, we present three new cases of PPV with features of the cesioflammea and/or cesiomarmorata subtypes and mosaic mutations in GNAQ or GNA11. To better understand the risk of potentially occult complications faced by such patients we additionally reviewed 176 cases published in the literature. We report the frequency of clinical findings, their patterns of co-occurrence as well as published recommendations for surveillance after diagnosis. Features assessed include: capillary malformation; dermal and ocular melanocytosis; glaucoma; limb asymmetry; venous malformations; and central nervous system (CNS) anomalies, such as ventriculomegaly and calcifications. We found that ocular findings are common in patients with phacomatosis cesioflammea and cesiomarmorata. Facial vascular involvement correlates with a higher risk of seizures (p = .0066). Our genetic results confirm the role of mosaic somatic mutations in GNAQ and GNA11 in phacomatosis cesioflammea and cesiomarmorata. Their clinical and molecular findings place these conditions on a clinical spectrum encompassing other GNAQ and GNA11 related disorders and inform recommendations for their management.

Atlas-CNV: a validated approach to call single-exon CNVs in the eMERGESeq gene panel

PURPOSE

To provide a validated method to confidently identify exon-containing copy-number variants (CNVs), with a low false discovery rate (FDR), in targeted sequencing data from a clinical laboratory with particular focus on single-exon CNVs.

METHODS

DNA sequence coverage data are normalized within each sample and subsequently exonic CNVs are identified in a batch of samples, when the target log₂ ratio of the sample to the batch median exceeds defined thresholds. The quality of exonic CNV calls is assessed by C-scores (Z-like scores) using thresholds derived from gold standard samples and simulation studies. We integrate an ExonQC threshold to lower FDR and compare performance with alternate software (VisCap).

RESULTS

Thirteen CNVs were used as a truth set to validate Atlas-CNV and compared with VisCap. We demonstrated FDR reduction in validation, simulation, and 10,926 eMERGESeq samples without sensitivity loss. Sixty-four multiexon and 29 single-exon CNVs with high C-scores were assessed by Multiplex Ligation-dependent Probe Amplification (MLPA).

CONCLUSION

Atlas-CNV is validated as a method to identify exonic CNVs in targeted sequencing data generated in the clinical laboratory. The ExonQC and C-score assignment can reduce FDR (identification of targets with high variance) and improve calling accuracy of single-exon CNVs respectively. We propose guidelines and criteria to identify high confidence single-exon CNVs.

Bi-allelic Variants in TONSL Cause SPONASTRIME Dysplasia and a Spectrum of Skeletal Dysplasia Phenotypes

SPONASTRIME dysplasia is an autosomal-recessive spondyloepimetaphyseal dysplasia characterized by spine (spondylar) abnormalities, midface hypoplasia with a depressed nasal bridge, metaphyseal striations, and disproportionate short stature. Scoliosis, coxa vara, childhood cataracts, short dental roots, and hypogammaglobulinemia have also been reported in this disorder. Although an autosomal-recessive inheritance pattern has been hypothesized, pathogenic variants in a specific gene have not been discovered in individuals with SPONASTRIME dysplasia. Here, we identified bi-allelic variants in TONSL, which encodes the Tonsoku-like DNA repair protein, in nine subjects (from eight families) with SPONASTRIME dysplasia, and four subjects (from three families) with short stature of varied severity and spondylometaphyseal dysplasia with or without immunologic and hematologic abnormalities, but no definitive metaphyseal striations at diagnosis. The finding of early embryonic lethality in a Tonsl-/- murine model and the discovery of reduced length, spinal abnormalities, reduced numbers of neutrophils, and early lethality in a tonsl-/- zebrafish model both support the hypomorphic nature of the identified TONSL variants. Moreover, functional studies revealed increased amounts of spontaneous replication fork stalling and chromosomal aberrations, as well as fewer camptothecin (CPT)-induced RAD51 foci in subject-derived cell lines. Importantly, these cellular defects were rescued upon re-expression of wild-type (WT) TONSL; this rescue is consistent with the hypothesis that hypomorphic TONSL variants are pathogenic. Overall, our studies in humans, mice, zebrafish, and subject-derived cell lines confirm that pathogenic variants in TONSL impair DNA replication and homologous recombination-dependent repair processes, and they lead to a spectrum of skeletal dysplasia phenotypes with numerous extra-skeletal manifestations.

De novo and inherited TCF20 pathogenic variants are associated with intellectual disability, dysmorphic features, hypotonia, and neurological impairments with similarities to Smith-Magenis syndrome

BACKGROUND

Neurodevelopmental disorders are genetically and phenotypically heterogeneous encompassing developmental delay (DD), intellectual disability (ID), autism spectrum disorders (ASDs), structural brain abnormalities, and neurological manifestations with variants in a large number of genes (hundreds) associated. To date, a few de novo mutations potentially disrupting TCF20 function in patients with ID, ASD, and hypotonia have been reported. TCF20 encodes a transcriptional co-regulator structurally related to RAI1, the dosage-sensitive gene responsible for Smith-Magenis syndrome (deletion/haploinsufficiency) and Potocki-Lupski syndrome (duplication/triplosensitivity).

METHODS

Genome-wide analyses by exome sequencing (ES) and chromosomal microarray analysis (CMA) identified individuals with heterozygous, likely damaging, loss-of-function alleles in TCF20. We implemented further molecular and clinical analyses to determine the inheritance of the pathogenic variant alleles and studied the spectrum of phenotypes.

RESULTS

We report 25 unique inactivating single nucleotide variants/indels (1 missense, 1 canonical splice-site variant, 18 frameshift, and 5 nonsense) and 4 deletions of TCF20. The pathogenic variants were detected in 32 patients and 4 affected parents from 31 unrelated families. Among cases with available parental samples, the variants were de novo in 20 instances and inherited from 4 symptomatic parents in 5, including in one set of monozygotic twins. Two pathogenic loss-of-function variants were recurrent in unrelated families. Patients presented with a phenotype characterized by developmental delay, intellectual disability, hypotonia, variable dysmorphic features, movement disorders, and sleep disturbances.

CONCLUSIONS

TCF20 pathogenic variants are associated with a novel syndrome manifesting clinical characteristics similar to those observed in Smith-Magenis syndrome. Together with previously described cases, the clinical entity of TCF20-associated neurodevelopmental disorders (TAND) emerges from a genotype-driven perspective.

Non-invasive prenatal sequencing for multiple Mendelian monogenic disorders using circulating cell-free fetal DNA

Current non-invasive prenatal screening is targeted toward the detection of chromosomal abnormalities in the fetus^1,2. However, screening for many dominant monogenic disorders associated with de novo mutations is not available, despite their relatively high incidence³. Here we report on the development and validation of, and early clinical experience with, a new approach for non-invasive prenatal sequencing for a panel of causative genes for frequent dominant monogenic diseases. Cell-free DNA (cfDNA) extracted from maternal plasma was barcoded, enriched, and then analyzed by next-generation sequencing (NGS) for targeted regions. Low-level fetal variants were identified by a statistical analysis adjusted for NGS read count and fetal fraction. Pathogenic or likely pathogenic variants were confirmed by a secondary amplicon-based test on cfDNA. Clinical tests were performed on 422 pregnancies with or without abnormal ultrasound findings or family history. Follow-up studies on cases with available outcome results confirmed 20 true-positive, 127 true-negative, zero false-positive, and zero-false negative results. The initial clinical study demonstrated that this non-invasive test can provide valuable molecular information for the detection of a wide spectrum of dominant monogenic diseases, complementing current screening for aneuploidies or carrier screening for recessive disorders.

Expanding the Spectrum of BAF-Related Disorders: De Novo Variants in SMARCC2 Cause a Syndrome with Intellectual Disability and Developmental Delay

SMARCC2 (BAF170) is one of the invariable core subunits of the ATP-dependent chromatin remodeling BAF (BRG1-associated factor) complex and plays a crucial role in embryogenesis and corticogenesis. Pathogenic variants in genes encoding other components of the BAF complex have been associated with intellectual disability syndromes. Despite its significant biological role, variants in SMARCC2 have not been directly associated with human disease previously. Using whole-exome sequencing and a web-based gene-matching program, we identified 15 individuals with variable degrees of neurodevelopmental delay and growth retardation harboring one of 13 heterozygous variants in SMARCC2, most of them novel and proven de novo. The clinical presentation overlaps with intellectual disability syndromes associated with other BAF subunits, such as Coffin-Siris and Nicolaides-Baraitser syndromes and includes prominent speech impairment, hypotonia, feeding difficulties, behavioral abnormalities, and dysmorphic features such as hypertrichosis, thick eyebrows, thin upper lip vermilion, and upturned nose. Nine out of the fifteen individuals harbor variants in the highly conserved SMARCC2 DNA-interacting domains (SANT and SWIRM) and present with a more severe phenotype. Two of these individuals present cardiac abnormalities. Transcriptomic analysis of fibroblasts from affected individuals highlights a group of differentially expressed genes with possible roles in regulation of neuronal development and function, namely H19, SCRG1, RELN, and CACNB4. Our findings suggest a novel SMARCC2-related syndrome that overlaps with neurodevelopmental disorders associated with variants in BAF-complex subunits.

Biallelic loss of function variants in PPP1R21 cause a neurodevelopmental syndrome with impaired endocytic function

Next-generation sequencing (NGS) has been instrumental in solving the genetic basis of rare inherited diseases, especially neurodevelopmental syndromes. However, functional workup is essential for precise phenotype definition and to understand the underlying disease mechanisms. Using whole exome (WES) and whole genome sequencing (WGS) in four independent families with hypotonia, neurodevelopmental delay, facial dysmorphism, loss of white matter, and thinning of the corpus callosum, we identified four previously unreported homozygous truncating PPP1R21 alleles: c.347delT p.(Ile116Lysfs*25), c.2170_2171insGGTA p.(Ile724Argfs*8), c.1607dupT p.(Leu536Phefs*7), c.2063delA p.(Lys688Serfs*26) and found that PPP1R21 was absent in fibroblasts of an affected individual, supporting the allele's loss of function effect. PPP1R21 function had not been studied except that a large scale affinity proteomics approach suggested an interaction with PIBF1 defective in Joubert syndrome. Our co-immunoprecipitation studies did not confirm this but in contrast defined the localization of PPP1R21 to the early endosome. Consistent with the subcellular expression pattern and the clinical phenotype exhibiting features of storage diseases, we found patient fibroblasts exhibited a delay in clearance of transferrin-488 while uptake was normal. In summary, we delineate a novel neurodevelopmental syndrome caused by biallelic PPP1R21 loss of function variants, and suggest a role of PPP1R21 within the endosomal sorting process or endosome maturation pathway.

Effect of Genetic Diagnosis on Patients with Previously Undiagnosed Disease

BACKGROUND

Many patients remain without a diagnosis despite extensive medical evaluation. The Undiagnosed Diseases Network (UDN) was established to apply a multidisciplinary model in the evaluation of the most challenging cases and to identify the biologic characteristics of newly discovered diseases. The UDN, which is funded by the National Institutes of Health, was formed in 2014 as a network of seven clinical sites, two sequencing cores, and a coordinating center. Later, a central biorepository, a metabolomics core, and a model organisms screening center were added.

METHODS

We evaluated patients who were referred to the UDN over a period of 20 months. The patients were required to have an undiagnosed condition despite thorough evaluation by a health care provider. We determined the rate of diagnosis among patients who subsequently had a complete evaluation, and we observed the effect of diagnosis on medical care.

RESULTS

A total of 1519 patients (53% female) were referred to the UDN, of whom 601 (40%) were accepted for evaluation. Of the accepted patients, 192 (32%) had previously undergone exome sequencing. Symptoms were neurologic in 40% of the applicants, musculoskeletal in 10%, immunologic in 7%, gastrointestinal in 7%, and rheumatologic in 6%. Of the 382 patients who had a complete evaluation, 132 received a diagnosis, yielding a rate of diagnosis of 35%. A total of 15 diagnoses (11%) were made by clinical review alone, and 98 (74%) were made by exome or genome sequencing. Of the diagnoses, 21% led to recommendations regarding changes in therapy, 37% led to changes in diagnostic testing, and 36% led to variant-specific genetic counseling. We defined 31 new syndromes.

CONCLUSIONS

The UDN established a diagnosis in 132 of the 382 patients who had a complete evaluation, yielding a rate of diagnosis of 35%. (Funded by the National Institutes of Health Common Fund.).

Clinical exome sequencing for fetuses with ultrasound abnormalities and a suspected Mendelian disorder

Background

Exome sequencing is now being incorporated into clinical care for pediatric and adult populations, but its integration into prenatal diagnosis has been more limited. One reason for this is the paucity of information about the clinical utility of exome sequencing in the prenatal setting.

Methods

We retrospectively reviewed indications, results, time to results (turnaround time, TAT), and impact of exome results for 146 consecutive “fetal exomes” performed in a clinical diagnostic laboratory between March 2012 and November 2017. We define a fetal exome as one performed on a sample obtained from a fetus or a product of conception with at least one structural anomaly detected by prenatal imaging or autopsy. Statistical comparisons were performed using Fisher’s exact test.

Results

Prenatal exome yielded an overall molecular diagnostic rate of 32% (n = 46/146). Of the 46 molecular diagnoses, 50% were autosomal dominant disorders (n = 23/46), 41% were autosomal recessive disorders (n = 19/46), and 9% were X-linked disorders (n = 4/46). The molecular diagnostic rate was highest for fetuses with anomalies affecting multiple organ systems and for fetuses with craniofacial anomalies. Out of 146 cases, a prenatal trio exome option designed for ongoing pregnancies was performed on 62 fetal specimens, resulting in a diagnostic yield of 35% with an average TAT of 14 days for initial reporting (excluding tissue culture time). The molecular diagnoses led to refined recurrence risk estimates, altered medical management, and informed reproductive planning for families.

Conclusion

Exome sequencing is a useful diagnostic tool when fetal structural anomalies suggest a genetic etiology, but other standard prenatal genetic tests did not provide a diagnosis.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0582-x) contains supplementary material, which is available to authorized users.

Phenotypic expansion in DDX3X - a common cause of intellectual disability in females

De novo variants in DDX3X account for 1–3% of unexplained intellectual disability (ID) cases and are amongst the most common causes of ID especially in females. Forty‐seven patients (44 females, 3 males) have been described. We identified 31 additional individuals carrying 29 unique DDX3X variants, including 30 postnatal individuals with complex clinical presentations of developmental delay or ID, and one fetus with abnormal ultrasound findings. Rare or novel phenotypes observed include respiratory problems, congenital heart disease, skeletal muscle mitochondrial DNA depletion, and late‐onset neurologic decline. Our findings expand the spectrum of DNA variants and phenotypes associated with DDX3X disorders.

De novo apparent loss-of-function mutations in PRR12 in three patients with intellectual disability and iris abnormalities

PRR12 encodes a proline-rich protein nuclear factor suspected to be involved in neural development. Its nuclear expression in fetal brains and in the vision system supports its role in brain and eye development more specifically. However, its function and potential role in human disease has not been determined. Recently, a de novo t(10;19) (q22.3;q13.33) translocation disrupting the PRR12 gene was detected in a girl with intellectual disability and neuropsychiatric alterations. Here we report on three unrelated patients with heterozygous de novo apparent loss-of-function mutations in PRR12 detected by clinical whole exome sequencing: c.1918G>T (p.Glu640*), c.4502_4505delTGCC (p.Leu1501Argfs*146) and c.903_909dup (p.Pro304Thrfs*46). All three patients had global developmental delay, intellectual disability, eye and vision abnormalities, dysmorphic features, and neuropsychiatric problems. Eye abnormalities were consistent among the three patients and consisted of stellate iris pattern and iris coloboma. Additional variable clinical features included hypotonia, skeletal abnormalities, sleeping problems, and behavioral issues such as autism and anxiety. In summary, we propose that haploinsufficiency of PRR12 is associated with this novel multisystem neurodevelopmental disorder.

Two de novo novel mutations in one SHANK3 allele in a patient with autism and moderate intellectual disability

SHANK3 encodes for a scaffolding protein that links neurotransmitter receptors to the cytoskeleton and is enriched in postsynaptic densities of excitatory synapses. Deletions or mutations in one copy of the SHANK3 gene cause Phelan-McDermid syndrome, also called 22q13.3 deletion syndrome, a neurodevelopmental disorder with common features including global developmental delay, absent to severely impaired language, autistic behavior, and minor dysmorphic features. By whole exome sequencing, we identified two de novo novel variants including one frameshift pathogenic variant and one missense variant of unknown significance in a 14-year-old boy with delayed motor milestones, delayed language acquisition, autism, intellectual disability, ataxia, progressively worsening spasticity of the lower extremities, dysmorphic features, short stature, microcephaly, failure to thrive, chronic constipation, intrauterine growth restriction, and bilateral inguinal hernias. Both changes are within the CpG island in exon 21, separated by a 375 bp sequence. Next generation sequencing of PCR products revealed that the two variants are most frequently associated with each other. Sanger sequencing of the cloned PCR products further confirmed that both changes were on a single allele. The clinical presentation in this individual is consistent with other patients with a truncating mutation in exon 21, suggesting that the missense change contributes none or minimally to the phenotypes. This is the first report of two de novo mutations in one SHANK3 allele.

Clinical and molecular characterization of de novo loss of function variants in HNRNPU

DNA alterations in the 1q43-q44 region are associated with syndromic neurodevelopmental disorders characterized by global developmental delay, intellectual disability, dysmorphic features, microcephaly, seizures, and agenesis of the corpus callosum. HNRNPU is located within the 1q43-q44 region and mutations in the gene have been reported in patients with early infantile epileptic encephalopathy. Here, we report on the clinical presentation of four patients with de novo heterozygous HNRNPU loss-of-function mutations detected by clinical whole exome sequencing: c.651_660del (p.Gly218Alafs*118), c.1089G>A (p.Trp363*), c.1714C>T (p.Arg572*), and c.2270_2271del (p.Pro757Argfs*7). All patients shared similar clinical features as previously reported including seizures, global developmental delay, intellectual disability, variable neurologic regression, behavior issues, and dysmorphic facial features. Features including heart defects and kidney abnormalities were not reported in our patients. These findings expands the clinical spectrum of HNRNPU-related disorder and shows that HNRNPU contributes to a subset of the clinical phenotypes associated with the contiguous 1q43-q44 deletion syndrome.

Renal cell carcinoma harboring somatic TSC2 mutations in a child with methylmalonic acidemia

Pediatric renal cell carcinoma (RCC) is a rare cancer that can be associated with inherited diseases including tuberous sclerosis complex (TSC) caused by germline mutations in TSC1 or TSC2. Somatic mutations in TSC1 and TSC2 have also been reported in adult RCC, which predict response to mTOR inhibitors. Here, we present the first case of RCC in a child with methylmalonic acidemia (MMA). Clinical whole exome sequencing of blood and tumor samples confirmed the diagnosis of MMA and revealed two somatic inactivating mutations in TSC2, suggesting the potential consideration of an mTOR inhibitor in the event of tumor recurrence.

Biallelic Variants in OTUD6B Cause an Intellectual Disability Syndrome Associated with Seizures and Dysmorphic Features

Ubiquitination is a posttranslational modification that regulates many cellular processes including protein degradation, intracellular trafficking, cell signaling, and protein-protein interactions. Deubiquitinating enzymes (DUBs), which reverse the process of ubiquitination, are important regulators of the ubiquitin system. OTUD6B encodes a member of the ovarian tumor domain (OTU)-containing subfamily of deubiquitinating enzymes. Herein, we report biallelic pathogenic variants in OTUD6B in 12 individuals from 6 independent families with an intellectual disability syndrome associated with seizures and dysmorphic features. In subjects with predicted loss-of-function alleles, additional features include global developmental delay, microcephaly, absent speech, hypotonia, growth retardation with prenatal onset, feeding difficulties, structural brain abnormalities, congenital malformations including congenital heart disease, and musculoskeletal features. Homozygous Otud6b knockout mice were subviable, smaller in size, and had congenital heart defects, consistent with the severity of loss-of-function variants in humans. Analysis of peripheral blood mononuclear cells from an affected subject showed reduced incorporation of 19S subunits into 26S proteasomes, decreased chymotrypsin-like activity, and accumulation of ubiquitin-protein conjugates. Our findings suggest a role for OTUD6B in proteasome function, establish that defective OTUD6B function underlies a multisystemic human disorder, and provide additional evidence for the emerging relationship between the ubiquitin system and human disease.

Lessons learned from additional research analyses of unsolved clinical exome cases

BACKGROUND

Given the rarity of most single-gene Mendelian disorders, concerted efforts of data exchange between clinical and scientific communities are critical to optimize molecular diagnosis and novel disease gene discovery.

METHODS

We designed and implemented protocols for the study of cases for which a plausible molecular diagnosis was not achieved in a clinical genomics diagnostic laboratory (i.e. unsolved clinical exomes). Such cases were recruited to a research laboratory for further analyses, in order to potentially: (1) accelerate novel disease gene discovery; (2) increase the molecular diagnostic yield of whole exome sequencing (WES); and (3) gain insight into the genetic mechanisms of disease. Pilot project data included 74 families, consisting mostly of parent-offspring trios. Analyses performed on a research basis employed both WES from additional family members and complementary bioinformatics approaches and protocols.

RESULTS

Analysis of all possible modes of Mendelian inheritance, focusing on both single nucleotide variants (SNV) and copy number variant (CNV) alleles, yielded a likely contributory variant in 36% (27/74) of cases. If one includes candidate genes with variants identified within a single family, a potential contributory variant was identified in a total of ~51% (38/74) of cases enrolled in this pilot study. The molecular diagnosis was achieved in 30/63 trios (47.6%). Besides this, the analysis workflow yielded evidence for pathogenic variants in disease-associated genes in 4/6 singleton cases (66.6%), 1/1 multiplex family involving three affected siblings, and 3/4 (75%) quartet families. Both the analytical pipeline and the collaborative efforts between the diagnostic and research laboratories provided insights that allowed recent disease gene discoveries (PURA, TANGO2, EMC1, GNB5, ATAD3A, and MIPEP) and increased the number of novel genes, defined in this study as genes identified in more than one family (DHX30 and EBF3).

CONCLUSION

An efficient genomics pipeline in which clinical sequencing in a diagnostic laboratory is followed by the detailed reanalysis of unsolved cases in a research environment, supplemented with WES data from additional family members, and subject to adjuvant bioinformatics analyses including relaxed variant filtering parameters in informatics pipelines, can enhance the molecular diagnostic yield and provide mechanistic insights into Mendelian disorders. Implementing these approaches requires collaborative clinical molecular diagnostic and research efforts.

Haploinsufficiency of the E3 ubiquitin-protein ligase gene TRIP12 causes intellectual disability with or without autism spectrum disorders, speech delay, and dysmorphic features

Impairment of ubiquitin-proteasome system activity involving ubiquitin ligase genes UBE3A, UBE3B, and HUWE1 and deubiquitinating enzyme genes USP7 and USP9X has been reported in patients with neurodevelopmental delays. To date, only a handful of single-nucleotide variants (SNVs) and copy-number variants (CNVs) involving TRIP12, encoding a member of the HECT domain E3 ubiquitin ligases family on chromosome 2q36.3 have been reported. Using chromosomal microarray analysis and whole-exome sequencing (WES), we have identified, respectively, five deletion CNVs and four inactivating SNVs (two frameshifts, one missense, and one splicing) in TRIP12. Seven of these variants were found to be de novo; parental studies could not be completed in two families. Quantitative PCR analyses of the splicing mutation showed a dramatically decreased level of TRIP12 mRNA in the proband compared to the family controls, indicating a loss-of-function mechanism. The shared clinical features include intellectual disability with or without autistic spectrum disorders, speech delay, and facial dysmorphism. Our findings demonstrate that E3 ubiquitin ligase TRIP12 plays an important role in nervous system development and function. The nine presented pathogenic variants further document that TRIP12 haploinsufficiency causes a childhood-onset neurodevelopmental disorder. Finally, our data enable expansion of the phenotypic spectrum of ubiquitin-proteasome dependent disorders.

A de novo mutation in the X-linked PAK3 gene is the underlying cause of intellectual disability and macrocephaly in monozygotic twins

Pathogenic variants in theP21 protein (Cdc42/Rac)-activated kinase 3gene (PAK3) lead to a rare non syndromic X-linked intellectual disability. The protein encoded by this gene forms an activated complex with GTP-bound RAS-like (P21), CDC2 and RAC1 proteins which then mediates a variety of cellular processes. So far, mutations in PAK3 gene have been reported in few families affected with intellectual disability associated with neurological manifestations such as speech defect, behavioral problem, brain structural abnormalities, microcephaly and cerebral palsy. In this study whole exome sequencing revealed a de novo likely pathogenic variant in PAK3 gene in monozygotic twins presented with intellectual disability, speech delay, behavioral problems and macrocephaly. Macrocephaly was noticed in our patients from birth at 35 weeks of gestation. This aspect of the phenotype has not been previously reported in other documented cases with pathogenic mutations in PAK3 gene. Our findings extend the phenotype of this disorder to include macrocephaly and offers further clues to the importance of the serine/threonine-protein kinase 3 (PAK3) protein in brain development and function.

Resolution of Disease Phenotypes Resulting from Multilocus Genomic Variation

BACKGROUND

Whole-exome sequencing can provide insight into the relationship between observed clinical phenotypes and underlying genotypes.

METHODS

We conducted a retrospective analysis of data from a series of 7374 consecutive unrelated patients who had been referred to a clinical diagnostic laboratory for whole-exome sequencing; our goal was to determine the frequency and clinical characteristics of patients for whom more than one molecular diagnosis was reported. The phenotypic similarity between molecularly diagnosed pairs of diseases was calculated with the use of terms from the Human Phenotype Ontology.

RESULTS

A molecular diagnosis was rendered for 2076 of 7374 patients (28.2%); among these patients, 101 (4.9%) had diagnoses that involved two or more disease loci. We also analyzed parental samples, when available, and found that de novo variants accounted for 67.8% (61 of 90) of pathogenic variants in autosomal dominant disease genes and 51.7% (15 of 29) of pathogenic variants in X-linked disease genes; both variants were de novo in 44.7% (17 of 38) of patients with two monoallelic variants. Causal copy-number variants were found in 12 patients (11.9%) with multiple diagnoses. Phenotypic similarity scores were significantly lower among patients in whom the phenotype resulted from two distinct mendelian disorders that affected different organ systems (50 patients) than among patients with disorders that had overlapping phenotypic features (30 patients) (median score, 0.21 vs. 0.36; P=1.77×10^-7).

CONCLUSIONS

In our study, we found multiple molecular diagnoses in 4.9% of cases in which whole-exome sequencing was informative. Our results show that structured clinical ontologies can be used to determine the degree of overlap between two mendelian diseases in the same patient; the diseases can be distinct or overlapping. Distinct disease phenotypes affect different organ systems, whereas overlapping disease phenotypes are more likely to be caused by two genes encoding proteins that interact within the same pathway. (Funded by the National Institutes of Health and the Ting Tsung and Wei Fong Chao Foundation.).

A Syndromic Neurodevelopmental Disorder Caused by De Novo Variants in EBF3

Early B cell factor 3 (EBF3) is a member of the highly evolutionarily conserved Collier/Olf/EBF (COE) family of transcription factors. Prior studies on invertebrate and vertebrate animals have shown that EBF3 homologs are essential for survival and that loss-of-function mutations are associated with a range of nervous system developmental defects, including perturbation of neuronal development and migration. Interestingly, aristaless-related homeobox (ARX), a homeobox-containing transcription factor critical for the regulation of nervous system development, transcriptionally represses EBF3 expression. However, human neurodevelopmental disorders related to EBF3 have not been reported. Here, we describe three individuals who are affected by global developmental delay, intellectual disability, and expressive speech disorder and carry de novo variants in EBF3. Associated features seen in these individuals include congenital hypotonia, structural CNS malformations, ataxia, and genitourinary abnormalities. The de novo variants affect a single conserved residue in a zinc finger motif crucial for DNA binding and are deleterious in a fly model. Our findings indicate that mutations in EBF3 cause a genetic neurodevelopmental syndrome and suggest that loss of EBF3 function might mediate a subset of neurologic phenotypes shared by ARX-related disorders, including intellectual disability, abnormal genitalia, and structural CNS malformations.

Primary immunodeficiency diseases: Genomic approaches delineate heterogeneous Mendelian disorders

BACKGROUND

Primary immunodeficiency diseases (PIDDs) are clinically and genetically heterogeneous disorders thus far associated with mutations in more than 300 genes. The clinical phenotypes derived from distinct genotypes can overlap. Genetic etiology can be a prognostic indicator of disease severity and can influence treatment decisions.

OBJECTIVE

We sought to investigate the ability of whole-exome screening methods to detect disease-causing variants in patients with PIDDs.

METHODS

Patients with PIDDs from 278 families from 22 countries were investigated by using whole-exome sequencing. Computational copy number variant (CNV) prediction pipelines and an exome-tiling chromosomal microarray were also applied to identify intragenic CNVs. Analytic approaches initially focused on 475 known or candidate PIDD genes but were nonexclusive and further tailored based on clinical data, family history, and immunophenotyping.

RESULTS

A likely molecular diagnosis was achieved in 110 (40%) unrelated probands. Clinical diagnosis was revised in about half (60/110) and management was directly altered in nearly a quarter (26/110) of families based on molecular findings. Twelve PIDD-causing CNVs were detected, including 7 smaller than 30 Kb that would not have been detected with conventional diagnostic CNV arrays.

CONCLUSION

This high-throughput genomic approach enabled detection of disease-related variants in unexpected genes; permitted detection of low-grade constitutional, somatic, and revertant mosaicism; and provided evidence of a mutational burden in mixed PIDD immunophenotypes.

Exome sequencing identifies de novo pathogenic variants in FBN1 and TRPS1 in a patient with a complex connective tissue phenotype

Here we describe a patient who presented with a history of congenital diaphragmatic hernia, inguinal hernia, and recurrent umbilical hernia. She also has joint laxity, hypotonia, and dysmorphic features. A unifying diagnosis was not identified based on her clinical phenotype. As part of her evaluation through the Undiagnosed Diseases Network, trio whole-exome sequencing was performed. Pathogenic variants in FBN1 and TRPS1 were identified as causing two distinct autosomal dominant conditions, each with de novo inheritance. Fibrillin 1 (FBN1) mutations are associated with Marfan syndrome and a spectrum of similar phenotypes. TRPS1 mutations are associated with trichorhinophalangeal syndrome types I and III. Features of both conditions are evident in the patient reported here. Discrepant features of the conditions (e.g., stature) and the young age of the patient may have made a clinical diagnosis more difficult in the absence of exome-wide genetic testing.

MIPEP recessive variants cause a syndrome of left ventricular non-compaction, hypotonia, and infantile death

BACKGROUND

Mitochondrial presequence proteases perform fundamental functions as they process about 70 % of all mitochondrial preproteins that are encoded in the nucleus and imported posttranslationally. The mitochondrial intermediate presequence protease MIP/Oct1, which carries out precursor processing, has not yet been established to have a role in human disease.

METHODS

Whole exome sequencing was performed on four unrelated probands with left ventricular non-compaction (LVNC), developmental delay (DD), seizures, and severe hypotonia. Proposed pathogenic variants were confirmed by Sanger sequencing or array comparative genomic hybridization. Functional analysis of the identified MIP variants was performed using the model organism Saccharomyces cerevisiae as the protein and its functions are highly conserved from yeast to human.

RESULTS

Biallelic single nucleotide variants (SNVs) or copy number variants (CNVs) in MIPEP, which encodes MIP, were present in all four probands, three of whom had infantile/childhood death. Two patients had compound heterozygous SNVs (p.L582R/p.L71Q and p.E602*/p.L306F) and one patient from a consanguineous family had a homozygous SNV (p.K343E). The fourth patient, identified through the GeneMatcher tool, a part of the Matchmaker Exchange Project, was found to have inherited a paternal SNV (p.H512D) and a maternal CNV (1.4-Mb deletion of 13q12.12) that includes MIPEP. All amino acids affected in the patients' missense variants are highly conserved from yeast to human and therefore S. cerevisiae was employed for functional analysis (for p.L71Q, p.L306F, and p.K343E). The mutations p.L339F (human p.L306F) and p.K376E (human p.K343E) resulted in a severe decrease of Oct1 protease activity and accumulation of non-processed Oct1 substrates and consequently impaired viability under respiratory growth conditions. The p.L83Q (human p.L71Q) failed to localize to the mitochondria.

CONCLUSIONS

Our findings reveal for the first time the role of the mitochondrial intermediate peptidase in human disease. Loss of MIP function results in a syndrome which consists of LVNC, DD, seizures, hypotonia, and cataracts. Our approach highlights the power of data exchange and the importance of an interrelationship between clinical and research efforts for disease gene discovery.