NODAL variants are associated with a continuum of laterality defects from simple D-transposition of the great arteries to heterotaxy

BACKGROUND

NODAL signaling plays a critical role in embryonic patterning and heart development in vertebrates. Genetic variants resulting in perturbations of the TGF-β/NODAL signaling pathway have reproducibly been shown to cause laterality defects in humans. To further explore this association and improve genetic diagnosis, the study aims to identify and characterize a broader range of NODAL variants in a large number of individuals with laterality defects.

METHODS

We re-analyzed a cohort of 321 proband-only exomes of individuals with clinically diagnosed laterality congenital heart disease (CHD) using family-based, rare variant genomic analyses. To this cohort we added 12 affected subjects with known NODAL variants and CHD from institutional research and clinical cohorts to investigate an allelic series. For those with candidate contributory variants, variant allele confirmation and segregation analysis were studied by Sanger sequencing in available family members. Array comparative genomic hybridization and droplet digital PCR were utilized for copy number variants (CNV) validation and characterization. We performed Human Phenotype Ontology (HPO)-based quantitative phenotypic analyses to dissect allele-specific phenotypic differences.

RESULTS

Missense, nonsense, splice site, indels, and/or structural variants of NODAL were identified as potential causes of heterotaxy and other laterality defects in 33 CHD cases. We describe a recurrent complex indel variant for which the nucleic acid secondary structure predictions implicate secondary structure mutagenesis as a possible mechanism for formation. We identified two CNV deletion alleles spanning NODAL in two unrelated CHD cases. Furthermore, 17 CHD individuals were found (16/17 with known Hispanic ancestry) to have the c.778G > A:p.G260R NODAL missense variant which we propose reclassification from variant of uncertain significance (VUS) to likely pathogenic. Quantitative HPO-based analyses of the observed clinical phenotype for all cases with p.G260R variation, including heterozygous, homozygous, and compound heterozygous cases, reveal clustering of individuals with biallelic variation. This finding provides evidence for a genotypic-phenotypic correlation and an allele-specific gene dosage model.

CONCLUSION

Our data further support a role for rare deleterious variants in NODAL as a cause for sporadic human laterality defects, expand the repertoire of observed anatomical complexity of potential cardiovascular anomalies, and implicate an allele specific gene dosage model.

Genomic Balancing Act: Deciphering DNA rearrangements in the Complex Chromosomal Aberration involving 5p15.2, 2q31.1 and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long read genome sequencing unveiled the presence of four distinct junctions in the father, who is carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndrome and to better understand genetic disease mechanisms.

Intrafamilial phenotypic heterogeneity in siblings with pseudohypoparathyroidism 1B due to maternal STX16 deletion

OBJECTIVES

Pseudohypoparathyroidism (PHP1B) is most commonly caused by epigenetic defects resulting in loss of methylation at the GNAS locus, although deletions of STX16 leading to GNAS methylation abnormalities have been previously reported. The phenotype of this disorder is variable and can include hormonal resistances and severe infantile obesity with hyperphagia. A possible time relationship between the onset of obesity and endocrinopathies has been previously reported but remains unclear. Understanding of the condition's natural history is limited, partly due to a scarcity of literature, especially in children.

CASE PRESENTATION

We report three siblings with autosomal dominant PHP1B caused by a deletion in STX16 who presented with early childhood onset PTH-resistance with normocalcemia with a progressive nature, accompanied by TSH-resistance and severe infantile obesity with hyperphagia in some, not all of the affected individuals.

CONCLUSIONS

PHP1B from a STX16 deletion displays intrafamilial phenotypic variation. It is a novel cause of severe infantile obesity, which is not typically included in commercially available gene panels but must be considered in the genetic work-up. Finally, it does not seem to have a clear time relationship between the onset of obesity and hormonal resistance.

RAB1A haploinsufficiency phenocopies the 2p14-p15 microdeletion and is associated with impaired neuronal differentiation

Hereditary spastic parapareses (HSPs) are clinically heterogeneous motor neuron diseases with variable age of onset and severity. Although variants in dozens of genes are implicated in HSPs, much of the genetic basis for pediatric-onset HSP remains unexplained. Here, we re-analyzed clinical exome-sequencing data from siblings with HSP of unknown genetic etiology and identified an inherited nonsense mutation (c.523C>T [p.Arg175Ter]) in the highly conserved RAB1A. The mutation is predicted to produce a truncated protein with an intact RAB GTPase domain but without two C-terminal cysteine residues required for proper subcellular protein localization. Additional RAB1A mutations, including two frameshift mutations and a mosaic missense mutation (c.83T>C [p.Leu28Pro]), were identified in three individuals with similar neurodevelopmental presentations. In rescue experiments, production of the full-length, but not the truncated, RAB1a rescued Golgi structure and cell proliferation in Rab1-depleted cells. In contrast, the missense-variant RAB1a disrupted Golgi structure despite intact Rab1 expression, suggesting a dominant-negative function of the mosaic missense mutation. Knock-down of RAB1A in cultured human embryonic stem cell-derived neurons resulted in impaired neuronal arborization. Finally, RAB1A is located within the 2p14-p15 microdeletion syndrome locus. The similar clinical presentations of individuals with RAB1A loss-of-function mutations and the 2p14-p15 microdeletion syndrome implicate loss of RAB1A in the pathogenesis of neurodevelopmental manifestations of this microdeletion syndrome. Our study identifies a RAB1A-related neurocognitive disorder with speech and motor delay, demonstrates an essential role for RAB1a in neuronal differentiation, and implicates RAB1A in the etiology of the neurodevelopmental sequelae associated with the 2p14-p15 microdeletion syndrome.

Chromosomal microarray analysis supplements exome sequencing to diagnose children with suspected inborn errors of immunity

Purpose

Though copy number variants (CNVs) have been suggested to play a significant role in inborn errors of immunity (IEI), the precise nature of this role remains largely unexplored. We sought to determine the diagnostic contribution of CNVs using genome-wide chromosomal microarray analysis (CMA) in children with IEI.

Methods

We performed exome sequencing (ES) and CMA for 332 unrelated pediatric probands referred for evaluation of IEI. The analysis included primary, secondary, and incidental findings.

Results

Of the 332 probands, 134 (40.4%) received molecular diagnoses. Of these, 116/134 (86.6%) were diagnosed by ES alone. An additional 15/134 (11.2%) were diagnosed by CMA alone, including two likely de novo changes. Three (2.2%) participants had diagnostic molecular findings from both ES and CMA, including two compound heterozygotes and one participant with two distinct diagnoses. Half of the participants with CMA contribution to diagnosis had CNVs in at least one non-immune gene, highlighting the clinical complexity of these cases. Overall, CMA contributed to 18/134 diagnoses (13.4%), increasing the overall diagnostic yield by 15.5% beyond ES alone.

Conclusion

Pairing ES and CMA can provide a comprehensive evaluation to clarify the complex factors that contribute to both immune and non-immune phenotypes. Such a combined approach to genetic testing helps untangle complex phenotypes, not only by clarifying the differential diagnosis, but in some cases by identifying multiple diagnoses contributing to the overall clinical presentation.

Expansion of the clinical and molecular spectrum of WWOX-related epileptic encephalopathy

WWOX biallelic loss-of-function pathogenic single nucleotide variants (SNVs) and copy number variants (CNVs) including exonic deletions and duplications cause WWOX-related epileptic encephalopathy (WOREE) syndrome. This disorder is characterized by refractory epilepsy, axial hypotonia, peripheral hypertonia, progressive microcephaly, and premature death. Here we report five patients with WWOX biallelic predicted null variants identified by exome sequencing (ES), genome sequencing (GS), and/or chromosomal microarray analysis (CMA). SNVs and intragenic deletions of one or more exons were commonly reported in WOREE syndrome patients which made the genetic diagnosis challenging and required a combination of different diagnostic technologies. These patients presented with severe, developmental and epileptic encephalopathy (DEE), and other cardinal features consistent with WOREE syndrome. This report expands the clinical phenotype associated with this condition, including failure to thrive in most patients and epilepsy that responded to a ketogenic diet in three patients. Dysmorphic features and abnormal prenatal findings were not commonly observed. Additionally, recurrent pancreatitis and sensorineural hearing loss each were observed in single patients. In summary, these phenotypic features broaden the clinical spectrum of WOREE syndrome.

Compound heterozygosity of a de novo submicroscopic deletion and an inherited frameshift pathogenic variant in the PKHD1 gene in a fetus with bilaterally enlarged and echogenic kidneys, enlarged abdomen and oligohydramnios

We present a fetus with bilaterally enlarged and echogenic kidneys. Prenatal testing detected compound heterozygosity for a 0.676 Mb de novo deletion and an inherited pathogenic variant in PKHD1. This is the first case of autosomal recessive polycystic kidney disease (ARPKD) with a prenatally detected disease-causing PKHD1 deletion.

Alternative polyadenylation alters protein dosage by switching between intronic and 3'UTR sites

Alternative polyadenylation (APA) creates distinct transcripts from the same gene by cleaving the pre-mRNA at poly(A) sites that can lie within the 3' untranslated region (3'UTR), introns, or exons. Most studies focus on APA within the 3'UTR; however, here, we show that CPSF6 insufficiency alters protein levels and causes a developmental syndrome by deregulating APA throughout the transcript. In neonatal humans and zebrafish larvae, CPSF6 insufficiency shifts poly(A) site usage between the 3'UTR and internal sites in a pathway-specific manner. Genes associated with neuronal function undergo mostly intronic APA, reducing their expression, while genes associated with heart and skeletal function mostly undergo 3'UTR APA and are up-regulated. This suggests that, under healthy conditions, cells toggle between internal and 3'UTR APA to modulate protein expression.

TCEAL1 loss-of-function results in an X-linked dominant neurodevelopmental syndrome and drives the neurological disease trait in Xq22.2 deletions

An Xq22.2 region upstream of PLP1 has been proposed to underly a neurological disease trait when deleted in 46,XX females. Deletion mapping revealed that heterozygous deletions encompassing the smallest region of overlap (SRO) spanning six Xq22.2 genes (BEX3, RAB40A, TCEAL4, TCEAL3, TCEAL1, and MORF4L2) associate with an early-onset neurological disease trait (EONDT) consisting of hypotonia, intellectual disability, neurobehavioral abnormalities, and dysmorphic facial features. None of the genes within the SRO have been associated with monogenic disease in OMIM. Through local and international collaborations facilitated by GeneMatcher and Matchmaker Exchange, we have identified and herein report seven de novo variants involving TCEAL1 in seven unrelated families: three hemizygous truncating alleles; one hemizygous missense allele; one heterozygous TCEAL1 full gene deletion; one heterozygous contiguous deletion of TCEAL1, TCEAL3, and TCEAL4; and one heterozygous frameshift variant allele. Variants were identified through exome or genome sequencing with trio analysis or through chromosomal microarray. Comparison with previously reported Xq22 deletions encompassing TCEAL1 identified a more-defined syndrome consisting of hypotonia, abnormal gait, developmental delay/intellectual disability especially affecting expressive language, autistic-like behavior, and mildly dysmorphic facial features. Additional features include strabismus, refractive errors, variable nystagmus, gastroesophageal reflux, constipation, dysmotility, recurrent infections, seizures, and structural brain anomalies. An additional maternally inherited hemizygous missense allele of uncertain significance was identified in a male with hypertonia and spasticity without syndromic features. These data provide evidence that TCEAL1 loss of function causes a neurological rare disease trait involving significant neurological impairment with features overlapping the EONDT phenotype in females with the Xq22 deletion.

Exome sequencing efficacy and phenotypic expansions involving esophageal atresia/tracheoesophageal fistula plus

Esophageal atresia/tracheoesophageal fistula (EA/TEF) is a life-threatening birth defect that often occurs with other major birth defects (EA/TEF+). Despite advances in genetic testing, a molecular diagnosis can only be made in a minority of EA/TEF+ cases. Here, we analyzed clinical exome sequencing data and data from the DECIPHER database to determine the efficacy of exome sequencing in cases of EA/TEF+ and to identify phenotypic expansions involving EA/TEF. Among 67 individuals with EA/TEF+ referred for clinical exome sequencing, a definitive or probable diagnosis was made in 11 cases for an efficacy rate of 16% (11/67). This efficacy rate is significantly lower than that reported for other major birth defects, suggesting that polygenic, multifactorial, epigenetic, and/or environmental factors may play a particularly important role in EA/TEF pathogenesis. Our cohort included individuals with pathogenic or likely pathogenic variants that affect TCF4 and its downstream target NRXN1, and FANCA, FANCB, and FANCC, which are associated with Fanconi anemia. These cases, previously published case reports, and comparisons to other EA/TEF genes made using a machine learning algorithm, provide evidence in support of a potential pathogenic role for these genes in the development of EA/TEF.

Clinical exome sequencing uncovers a high frequency of Mendelian disorders in infants with stroke: A retrospective analysis

Stroke causes significant disability and is a common cause of death worldwide. Previous studies have estimated that 1%-5% of stroke is attributable to monogenic etiologies. We set out to assess the utility of clinical exome sequencing (ES) in the evaluation of stroke. We retrospectively analyzed 124 individuals who received ES at the Baylor Genetics reference lab between 2012 and 2021 who had stroke as a major part of their reported phenotype. Ages ranged from 10 days to 69 years. 8.9% of the cohort received a diagnosis, including 25% of infants less than 1 year old; an additional 10.5% of the cohort received a probable diagnosis. We identified several syndromes that predispose to stroke such as COL4A1-related brain small vessel disease, homocystinuria caused by CBS mutation, POLG-related disorders, TTC19-linked mitochondrial disease, and RNASEH2A associated Aicardi-Goutieres syndrome. We also observed pathogenic variants in NSD1, PKHD1, HRAS, and ATP13A2, which are genes rarely associated with stroke. Although stroke is a complex phenotype with varying pathologies and risk factors, these results show that use of exome sequencing can be highly relevant in stroke, especially for those presenting <1 year of age.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation

BACKGROUND

The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases.

METHODS

Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype.

RESULTS

In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations.

CONCLUSIONS

Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Sequencing individual genomes with recurrent genomic disorder deletions: an approach to characterize genes for autosomal recessive rare disease traits

BACKGROUND

In medical genetics, discovery and characterization of disease trait contributory genes and alleles depends on genetic reasoning, study design, and patient ascertainment; we suggest a segmental haploid genetics approach to enhance gene discovery and molecular diagnostics.

METHODS

We constructed a genome-wide map for nonallelic homologous recombination (NAHR)-mediated recurrent genomic deletions and used this map to estimate population frequencies of NAHR deletions based on large-scale population cohorts and region-specific studies. We calculated recessive disease carrier burden using high-quality pathogenic or likely pathogenic variants from ClinVar and gnomAD. We developed a NIRD (NAHR deletion Impact to Recessive Disease) score for recessive disorders by quantifying the contribution of NAHR deletion to the overall allele load that enumerated all pairwise combinations of disease-causing alleles; we used a Punnett square approach based on an assumption of random mating. Literature mining was conducted to identify all reported patients with defects in a gene with a high NIRD score; meta-analysis was performed on these patients to estimate the representation of NAHR deletions in recessive traits from contemporary human genomics studies. Retrospective analyses of extant clinical exome sequencing (cES) were performed for novel rare recessive disease trait gene and allele discovery from individuals with NAHR deletions.

RESULTS

We present novel genomic insights regarding the genome-wide impact of NAHR recurrent segmental variants on recessive disease burden; we demonstrate the utility of NAHR recurrent deletions to enhance discovery in the challenging context of autosomal recessive (AR) traits and biallelic variation. Computational results demonstrate new mutations mediated by NAHR, involving recurrent deletions at 30 genomic regions, likely drive recessive disease burden for over 74% of loci within these segmental deletions or at least 2% of loci genome-wide. Meta-analyses on 170 literature-reported patients implicate that NAHR deletions are depleted from the ascertained pool of AR trait alleles. Exome reanalysis of personal genomes from subjects harboring recurrent deletions uncovered new disease-contributing variants in genes including COX10, ERCC6, PRRT2, and OTUD7A.

CONCLUSIONS

Our results demonstrate that genomic sequencing of personal genomes with NAHR deletions could dramatically improve allele and gene discovery and enhance clinical molecular diagnosis. Moreover, results suggest NAHR events could potentially enable human haploid genetic screens as an approach to experimental inquiry into disease biology.

Recurring germline mosaicism in a family due to reversion of an inherited derivative chromosome 8 from an 8;21 translocation with interstitial telomeric sequences

BACKGROUND

Mosaicism for chromosomal structural abnormalities, other than marker or ring chromosomes, is rarely inherited.

METHODS

We performed cytogenetics studies and breakpoint analyses on a family with transmission of mosaicism for a derivative chromosome 8 (der(8)), resulting from an unbalanced translocation between the long arms of chromosomes 8 and 21 over three generations.

RESULTS

The proband and his maternal half-sister had mosaicism for a der(8) cell line leading to trisomy of the distal 21q, and both had Down syndrome phenotypic features. Mosaicism for a cell line with the der(8) and a normal cell line was also detected in a maternal half-cousin. The der(8) was inherited from the maternal grandmother who had four abnormal cell lines containing the der(8), in addition to a normal cell line. One maternal half-aunt had the der(8) and an isodicentric chromosome 21 (idic(21)). Sequencing studies revealed microhomologies at the junctures of the der(8) and idic(21) in the half-aunt, suggesting a replicative mechanism in the rearrangement formation. Furthermore, interstitial telomeric sequences (ITS) were identified in the juncture between chromosomes 8 and 21 in the der(8).

CONCLUSION

Mosaicism in the proband, his half-sister and half-cousin resulting from loss of chromosome 21 material from the der(8) appears to be a postzygotic event due to the genomic instability of ITS and associated with selective growth advantage of normal cells. The reversion of the inherited der(8) to a normal chromosome 8 in this family resembles revertant mosaicism of point mutations. We propose that ITS could mediate recurring revertant mosaicism for some constitutional chromosomal structural abnormalities.

Clinical Exome Sequencing of 1000 Families with Complex Immune Phenotypes: Towards comprehensive genomic evaluations

BACKGROUND

Prospective genetic evaluation of patients at our referral research hospital presents clinical research challenges.

OBJECTIVE

This study sought not only a single-gene explanation for participants' immune-related presentations, but viewed each participant holistically, with the potential to have multiple genetic contributions to their immune-phenotype and other heritable comorbidities relevant to their presentation and health.

METHODS

We developed a program integrating exome sequencing, chromosomal microarray, phenotyping, results return with genetic counseling, and reanalysis in 1505 individuals from 1000 families with suspected or known inborn errors of immunity.

RESULTS

Probands were 50.8% female, 71.5% >18 years, and had diverse immune presentations. Overall, 327/1000 probands (32.7%) received 361 molecular diagnoses. These included 17 probands with diagnostic copy number variants, 32 probands with secondary findings, and 31 probands with multiple molecular diagnoses. Reanalysis added 22 molecular diagnoses, predominantly due to new disease-gene associations (9/22, 40.9%). One-quarter of the molecular diagnoses (92/361) did not involve immune-associated genes. Molecular diagnosis was correlated with younger age, male sex, and a higher number of organ systems involved. This program also facilitated the discovery of new gene-disease associations such as SASH3-related immunodeficiency. A review of treatment options and ClinGen actionability curations suggest that at least 251/361 (69.5%) of these molecular diagnoses could translate into >1 management option.

CONCLUSION

This program contributes to our understanding of the diagnostic and clinical utility whole exome analysis on a large scale.

CLINICAL

Implication: Comprehensive analysis of exome data has diagnostic and clinical utility for patients with suspected inborn errors of immunity.

LMOD2-related dilated cardiomyopathy presenting in late infancy

Leiomodin-2 (LMOD2) is an important regulator of the thin filament length, known to promote elongation of actin through polymerization at pointed ends. Mice with Lmod2 deficiency die around 3 weeks of age due to severe dilated cardiomyopathy (DCM), resulting from decreased heart contractility due to shorter thin filaments. To date, there have been three infants from two families reported with biallelic variants in LMOD2, presenting with perinatal onset DCM. Here, we describe a third family with a child harboring a previously described homozygous frameshift variant, c.1243_1244delCT (p.L415Vfs*108) with DCM, presenting later in infancy at 9 months of age. Family history was relevant for a sibling who died suddenly at 1 year of age after being diagnosed with cardiomegaly. LMOD2-related cardiomyopathy is a rare form of inherited cardiomyopathy resulting from thin filament length dysregulation and should be considered in genetic evaluation of newborns and infants with suspected autosomal recessive inheritance or sporadic early onset cardiomyopathy.

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.

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.

A rare description of pure partial trisomy of 16q12.2q24.3 and review of the literature

Trisomy 16 is the most common autosomal trisomy in humans, which is almost uniformly embryonic lethal. Partial trisomy 16 including a segment of the long arm of chromosome 16 is occasionally compatible with life and has been associated with severe congenital defects, growth retardation, and early lethality. Segmental trisomy of 16q is usually described concomitantly with partial monosomy of another chromosome, often resulting from a parental balanced translocation. Pure partial chromosome 16q trisomy is exceedingly rare. About nine children with 16q12→qter and 16q13→qter duplication have been reported in the literature, almost all described with monosomy of a second chromosome, and highlighting very few long-term survivors. A single individual with pure partial distal 16q12.1q23.3 duplication has been reported in an infant, underscoring complexities of genetic counseling and management, especially in view of life-limiting congenital anomalies in rare survivors. Here, we present a 12-month-old child with pure 16q12.2q24.3 trisomy, having continued morbidity related to pulmonary hypertension and chronic lung disease. The features of intrauterine growth retardation, facial dysmorphism, hypotonia, congenital heart defect, distal contractures, urogenital abnormalities, and hearing loss support the association with 16q partial trisomy, as in previous studies. This report expands our current understanding related to the survival of infants with large segmental aneusomy of the long arm of chromosome 16.

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.

Quantitative Assessment of Parental Somatic Mosaicism for Copy-Number Variant (CNV) Deletions

As genome sequencing methodologies have become more sensitive in detecting low-frequency rare-variant events, the link between post-zygotic mutagenesis and somatic mosaicism in the etiology of several human genetic conditions other than cancers has become more clear. Given that current clinical-genomics diagnostic methods have limited detection sensitivity for mosaic events, a copy-number variant (CNV) deletion inherited from a parent with low-level (<10%) mosaicism can be erroneously interpreted in the proband to represent a de novo germline event. Here, we describe three sensitive, precise, and cost-efficient methods that can quantitatively assess the potential degree of parental somatic mosaicism levels for CNV deletions: droplet digital PCR (ddPCR), PCR amplicon-based next-generation sequencing (NGS), and quantitative PCR. ddPCR using the EvaGreen fluorescent dye protocol can specifically quantify the deleted or non-deleted alleles by analyzing the number of droplets positive for a fluorescent signal for each event. PCR amplicon-based NGS assesses the allele frequencies of a heterozygous single-nucleotide polymorphism within a deletion region. The difference in number of reads between the two genotypes indicates the level of somatic mosaicism for the CNV deletion. Quantitative PCR can be applied where the relative quantity of the deletion junction-specific product represents the level of mosaicism. Clinical implementation of these quantitative variant-detection methods enables potentially more accurate assessment of disease recurrence risk in family-based genetic counseling, allowing couples to engage in more informed family planning. © 2020 by John Wiley & Sons, Inc. Basic Protocol: Droplet digital PCR (ddPCR) Alternate Protocol 1: PCR amplicon-based next-generation sequencing Alternate Protocol 2: Quantitative real-time PCR (qPCR).

Clinical characterization of individuals with the distal 1q21.1 microdeletion

Distal 1q21.1 microdeletions have shown highly variable clinical expressivity and incomplete penetrance, with affected individuals manifesting a broad spectrum of nonspecific features. The goals of this study were to better describe the phenotypic spectrum of patients with distal 1q21.1 microdeletions and to compare the clinical features among affected individuals. We performed a retrospective chart review of 47 individuals with distal 1q21.1 microdeletions tested at a large clinical genetic testing laboratory, with most patients being clinically evaluated in the same children's hospital. Health information such as growth charts, results of imaging studies, developmental history, and progress notes were collected. Statistical analysis was performed using Fisher's exact test to compare clinical features among study subjects. Common features in our cohort include microcephaly (51.2%), seizures (29.8%), developmental delay (74.5%), failure to thrive (FTT) (68.1%), dysmorphic features (63.8%), and a variety of congenital anomalies such as cardiac abnormalities (23.4%) and genitourinary abnormalities (19.1%). Compared to prior literature, we found that seizures, brain anomalies, and FTT were more prevalent among our study cohort. Females were more likely than males to have microcephaly (p = 0.0199) and cardiac abnormalities (p = 0.0018). Based on existing genome-wide clinical testing results, at least a quarter of the cohort had additional genetic findings that may impact the phenotype of the individual. Our study represents the largest cohort of distal 1q21.1 microdeletion carriers available in the literature thus far, and it further illustrates the wide spectrum of clinical manifestations among symptomatic individuals. These results may allow for improved genetic counseling and management of affected individuals. Future studies may help to elucidate the underlying molecular mechanisms impacting the phenotypic variability observed with this microdeletion.

Integrated sequencing and array comparative genomic hybridization in familial Parkinson disease

Objective

To determine how single nucleotide variants (SNVs) and copy number variants (CNVs) contribute to molecular diagnosis in familial Parkinson disease (PD), we integrated exome sequencing (ES) and genome-wide array-based comparative genomic hybridization (aCGH) and further probed CNV structure to reveal mutational mechanisms.

Methods

We performed ES on 110 subjects with PD and a positive family history; 99 subjects were also evaluated using genome-wide aCGH. We interrogated ES and aCGH data for pathogenic SNVs and CNVs at Mendelian PD gene loci. We confirmed SNVs via Sanger sequencing and further characterized CNVs with custom-designed high-density aCGH, droplet digital PCR, and breakpoint sequencing.

Results

Using ES, we discovered individuals with known pathogenic SNVs in GBA (p.Glu365Lys, p.Thr408Met, p.Asn409Ser, and p.Leu483Pro) and LRRK2 (p.Arg1441Gly and p.Gly2019Ser). Two subjects were each double heterozygotes for variants in GBA and LRRK2. Based on aCGH, we additionally discovered cases with an SNCA duplication and heterozygous intragenic GBA deletion. Five additional subjects harbored both SNVs (p.Asn52Metfs*29, p.Thr240Met, p.Pro437Leu, and p.Trp453*) and likely disrupting CNVs at the PRKN locus, consistent with compound heterozygosity. In nearly all cases, breakpoint sequencing revealed microhomology, a mutational signature consistent with CNV formation due to DNA replication errors.

Conclusions

Integrated ES and aCGH yielded a genetic diagnosis in 19.3% of our familial PD cohort. Our analyses highlight potential mechanisms for SNCA and PRKN CNV formation, uncover multilocus pathogenic variation, and identify novel SNVs and CNVs for further investigation as potential PD risk alleles.

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.

Parental somatic mosaicism for CNV deletions - A need for more sensitive and precise detection methods in clinical diagnostics settings

To further assess the scale and level of parental somatic mosaicism, we queried the CMA database at Baylor Genetics. We selected 50 unrelated families where clinically relevant apparent de novo CNV-deletions were found in the affected probands. Parental blood samples screening using deletion junction-specific PCR revealed four parents with somatic mosaicism. Droplet digital PCR (ddPCR), qPCR, and amplicon-based next-generation sequencing (NGS) were applied to validate these findings. Using ddPCR levels of mosaicism ranged from undetectable to 18.5%. Amplicon-based NGS and qPCR for the father with undetectable mosaicism was able to detect mosaicism at 0.39%. In one mother, ddPCR analysis revealed 15.6%, 10.6%, 8.2%, and undetectable levels of mosaicism in her blood, buccal cells, saliva, and urine samples, respectively. Our data suggest that more sensitive and precise methods, e.g. CNV junction-specific LR-PCR, ddPCR, or qPCR may allow for a more refined assessment of the potential disease recurrence risk for an identified variant.

GARS-related disease in infantile spinal muscular atrophy: Implications for diagnosis and treatment

The majority of patients with spinal muscular atrophy (SMA) identified to date harbor a biallelic exonic deletion of SMN1. However, there have been reports of SMA-like disorders that are independent of SMN1, including those due to pathogenic variants in the glycyl-tRNA synthetase gene (GARS1). We report three unrelated patients with de novo variants in GARS1 that are associated with infantile-onset SMA (iSMA). Patients were ascertained during inpatient hospital evaluations for complications of neuropathy. Evaluations were completed as indicated for clinical care and management and informed consent for publication was obtained. One newly identified, disease-associated GARS1 variant, identified in two out of three patients, was analyzed by functional studies in yeast complementation assays. Genomic analyses by exome and/or gene panel and SMN1 copy number analysis of three patients identified two previously undescribed de novo missense variants in GARS1 and excluded SMN1 as the causative gene. Functional studies in yeast revealed that one of the de novo GARS1 variants results in a loss-of-function effect, consistent with other pathogenic GARS1 alleles. In sum, the patients' clinical presentation, assessments of previously identified GARS1 variants and functional assays in yeast suggest that the GARS1 variants described here cause iSMA. GARS1 variants have been previously associated with Charcot-Marie-Tooth disease (CMT2D) and distal SMA type V (dSMAV). Our findings expand the allelic heterogeneity of GARS-associated disease and support that severe early-onset SMA can be caused by variants in this gene. Distinguishing the SMA phenotype caused by SMN1 variants from that due to pathogenic variants in other genes such as GARS1 significantly alters approaches to treatment.

Pathogenic variants in TNRC6B cause a genetic disorder characterised by developmental delay/intellectual disability and a spectrum of neurobehavioural phenotypes including autism and ADHD

BACKGROUND

Rare variants in hundreds of genes have been implicated in developmental delay (DD), intellectual disability (ID) and neurobehavioural phenotypes. TNRC6B encodes a protein important for RNA silencing. Heterozygous truncating variants have been reported in three patients from large cohorts with autism, but no full phenotypic characterisation was described.

METHODS

Clinical and molecular characterisation was performed on 17 patients with TNRC6B variants. Clinical data were obtained by retrospective chart review, parent interviews, direct patient interaction with providers and formal neuropsychological evaluation.

RESULTS

Clinical findings included DD/ID (17/17) (speech delay in 94% (16/17), fine motor delay in 82% (14/17) and gross motor delay in 71% (12/17) of subjects), autism or autistic traits (13/17), attention deficit and hyperactivity disorder (ADHD) (11/17), other behavioural problems (7/17) and musculoskeletal findings (12/17). Other congenital malformations or clinical findings were occasionally documented. The majority of patients exhibited some dysmorphic features but no recognisable gestalt was identified. 17 heterozygous TNRC6B variants were identified in 12 male and five female unrelated subjects by exome sequencing (14), a targeted panel (2) and a chromosomal microarray (1). The variants were nonsense (7), frameshift (5), splice site (2), intragenic deletions (2) and missense (1).

CONCLUSIONS

Variants in TNRC6B cause a novel genetic disorder characterised by recurrent neurocognitive and behavioural phenotypes featuring DD/ID, autism, ADHD and other behavioural abnormalities. Our data highly suggest that haploinsufficiency is the most likely pathogenic mechanism. TNRC6B should be added to the growing list of genes of the RNA-induced silencing complex associated with ID/DD, autism and ADHD.

BAZ2B haploinsufficiency as a cause of developmental delay, intellectual disability, and autism spectrum disorder

The bromodomain adjacent to zinc finger 2B gene (BAZ2B) encodes a protein involved in chromatin remodeling. Loss of BAZ2B function has been postulated to cause neurodevelopmental disorders. To determine whether BAZ2B deficiency is likely to contribute to the pathogenesis of these disorders, we performed bioinformatics analyses that demonstrated a high level of functional convergence during fetal cortical development between BAZ2B and genes known to cause autism spectrum disorder (ASD) and neurodevelopmental disorder. We also found an excess of de novo BAZ2B loss-of-function variants in exome sequencing data from previously published cohorts of individuals with neurodevelopmental disorders. We subsequently identified seven additional individuals with heterozygous deletions, stop-gain, or de novo missense variants affecting BAZ2B. All of these individuals have developmental delay (DD), intellectual disability (ID), and/or ASD. Taken together, our findings suggest that haploinsufficiency of BAZ2B causes a neurodevelopmental disorder, whose cardinal features include DD, ID, and ASD.

Loss of Oxidation Resistance 1, OXR1, Is Associated with an Autosomal-Recessive Neurological Disease with Cerebellar Atrophy and Lysosomal Dysfunction

We report an early-onset autosomal-recessive neurological disease with cerebellar atrophy and lysosomal dysfunction. We identified bi-allelic loss-of-function (LoF) variants in Oxidative Resistance 1 (OXR1) in five individuals from three families; these individuals presented with a history of severe global developmental delay, current intellectual disability, language delay, cerebellar atrophy, and seizures. While OXR1 is known to play a role in oxidative stress resistance, its molecular functions are not well established. OXR1 contains three conserved domains: LysM, GRAM, and TLDc. The gene encodes at least six transcripts, including some that only consist of the C-terminal TLDc domain. We utilized Drosophila to assess the phenotypes associated with loss of mustard (mtd), the fly homolog of OXR1. Strong LoF mutants exhibit late pupal lethality or pupal eclosion defects. Interestingly, although mtd encodes 26 transcripts, severe LoF and null mutations can be rescued by a single short human OXR1 cDNA that only contains the TLDc domain. Similar rescue is observed with the TLDc domain of NCOA7, another human homolog of mtd. Loss of mtd in neurons leads to massive cell loss, early death, and an accumulation of aberrant lysosomal structures, similar to what we observe in fibroblasts of affected individuals. Our data indicate that mtd and OXR1 are required for proper lysosomal function; this is consistent with observations that NCOA7 is required for lysosomal acidification.

Validation Studies for Single Circulating Trophoblast Genetic Testing as a Form of Noninvasive Prenatal Diagnosis

It has long been appreciated that genetic analysis of fetal or trophoblast cells in maternal blood could revolutionize prenatal diagnosis. We implemented a protocol for single circulating trophoblast (SCT) testing using positive selection by magnetic-activated cell sorting and single-cell low-coverage whole-genome sequencing to detect fetal aneuploidies and copy-number variants (CNVs) at ∼1 Mb resolution. In 95 validation cases, we identified on average 0.20 putative trophoblasts/mL, of which 55% were of high quality and scorable for both aneuploidy and CNVs. We emphasize the importance of analyzing individual cells because some cells are apoptotic, in S-phase, or otherwise of poor quality. When two or more high-quality trophoblast cells were available for singleton pregnancies, there was complete concordance between all trophoblasts unless there was evidence of confined placental mosaicism. SCT results were highly concordant with available clinical data from chorionic villus sampling (CVS) or amniocentesis procedures. Although determining the exact sensitivity and specificity will require more data, this study further supports the potential for SCT testing to become a diagnostic prenatal test.

Biallelic mutation of FBXL7 suggests a novel form of Hennekam syndrome

Hennekam lymphangiectasia-lymphedema syndrome is an autosomal recessive disorder characterized by congenital lymphedema, intestinal lymphangiectasia, facial dysmorphism, and variable intellectual disability. Known disease genes include CCBE1, FAT4, and ADAMTS3. In a patient with clinically diagnosed Hennekam syndrome but without mutations or copy-number changes in the three known disease genes, we identified a homozygous single-exon deletion affecting FBXL7. Specifically, exon 3, which encodes the F-box domain and several leucine-rich repeats of FBXL7, is eliminated. Our analyses of databases representing >100,000 control individuals failed to identify biallelic loss-of-function variants in FBXL7. Published studies in Drosophila indicate Fbxl7 interacts with Fat, of which human FAT4 is an ortholog, and mutation of either gene yields similar morphological consequences. These data suggest that FBXL7 may be the fourth gene for Hennekam syndrome, acting via a shared pathway with FAT4.

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.

Homozygous loss-of-function variants of TASP1, a gene encoding an activator of the histone methyltransferases KMT2A and KMT2D, cause a syndrome of developmental delay, happy demeanor, distinctive facial features, and congenital anomalies

We report four unrelated children with homozygous loss-of-function variants in TASP1 and an overlapping phenotype comprising developmental delay with hypotonia and microcephaly, feeding difficulties with failure-to-thrive, recurrent respiratory infections, cardiovascular malformations, cryptorchidism, happy demeanor, and distinctive facial features. Two children had a homozygous founder deletion encompassing exons 5-11 of TASP1, the third had a homozygous missense variant, c.701 C>T (p.Thr234Met), affecting the active site of the encoded enzyme, and the fourth had a homozygous nonsense variant, c.199 C>T (p.Arg67*). TASP1 encodes taspase 1 (TASP1), which is responsible for cleaving, thus activating, the lysine methyltransferases KMT2A and KMT2D, which are essential for histone methylation and transcription regulation. The consistency of the phenotype, the critical biological function of TASP1, the deleterious nature of the TASP1 variants, and the overlapping features with Wiedemann-Steiner and Kabuki syndromes respectively caused by pathogenic variants in KMT2A and KMT2D all support that TASP1 is a disease-related gene.

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.

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.

Clinical exome sequencing reveals locus heterogeneity and phenotypic variability of cohesinopathies

PURPOSE

Defects in the cohesin pathway are associated with cohesinopathies, notably Cornelia de Lange syndrome (CdLS). We aimed to delineate pathogenic variants in known and candidate cohesinopathy genes from a clinical exome perspective.

METHODS

We retrospectively studied patients referred for clinical exome sequencing (CES, N = 10,698). Patients with causative variants in novel or recently described cohesinopathy genes were enrolled for phenotypic characterization.

RESULTS

Pathogenic or likely pathogenic single-nucleotide and insertion/deletion variants (SNVs/indels) were identified in established disease genes including NIPBL (N = 5), SMC1A (N = 14), SMC3 (N = 4), RAD21 (N = 2), and HDAC8 (N = 8). The phenotypes in this genetically defined cohort skew towards the mild end of CdLS spectrum as compared with phenotype-driven cohorts. Candidate or recently reported cohesinopathy genes were supported by de novo SNVs/indels in STAG1 (N = 3), STAG2 (N = 5), PDS5A (N = 1), and WAPL (N = 1), and one inherited SNV in PDS5A. We also identified copy-number deletions affecting STAG1 (two de novo, one of unknown inheritance) and STAG2 (one of unknown inheritance). Patients with STAG1 and STAG2 variants presented with overlapping features yet without characteristic facial features of CdLS.

CONCLUSION

CES effectively identified disease-causing alleles at the mild end of the cohensinopathy spectrum and enabled characterization of candidate disease genes.

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.

[The value of lean nystagmus and sitting to supine positioning nystagmus in the diagnosis of horizontal semicircular canal benign paroxysmal positional vertigo]

Objective: To explore the value of 1ean nystagmus and sitting to supine positioning nystagmus in the diagnosis of horizontal semicircular canal benign paroxysmal positional vertigo. Method: One hundred cases of patients with definitive diagnosis of horizontal semicircular canal benign paroxysmal positional vertigo were tested by 1ean nystagmus and sitting to supine positioning nystagmus and supine roll test. The affected side was recorded according to the nystagmus direction. After diagnosis, they were treated with canalith repositioning procedure. The canalith repositioning procedure was made according to the supine roll test result when 1ean nystagmus and sitting to supine positioning nystagmuscannot be induced. Furthermore，the canalith repositioning procedure was made according to the1ean nystagmus and sitting to supine positioning nystagmuswhenthe supine roll test cannot diagnose. The patients with canalolithiasis were randomly divided into two groups when both 1ean nystagmus and sitting to supine positioning nystagmusand the supine roll test can be induced. One group was treated with canalith repositioning procedure based on results of supine roll test and the other group was treated based on the results of1ean nystagmus and sitting to supine positioning nystagmus. The detection rate and diagnostic coincidence rate of 1ean nystagmus and sitting to supine positioning nystagmus were calculated and the shortterm outcome were evaluated one day after treatment. Result: The detection rate of 1ean nystagmus and sitting to supine positioning nystagmus was 83% and the coincidence rate with the roll test was 90.1%，respectively. There was no significant statistical difference between the treatment effect according to lean nystagmus and sitting to supine positioning nystagmus and supine roll test. Conclusion: Lean nystagmus and sitting to supine positioning nystagmus cannot be used alone in the diagnosis ofhorizontal semicircular canal benign paroxysmal positional vertigo. They need to be combined with supine roll test. However, when the supine roll test is difficult to demonstrate the affected side, it can be used as an auxiliary diagnostic method. When the results of 1ean nystagmus and sitting to supine positioning nystagmusand the supine roll test are inconsistent, repeat the test and conduct a comprehensive assessment to avoid missed diagnosis and misdiagnosis..

Pathogenic variants in USP7 cause a neurodevelopmental disorder with speech delays, altered behavior, and neurologic anomalies

PURPOSE

Haploinsufficiency of USP7, located at chromosome 16p13.2, has recently been reported in seven individuals with neurodevelopmental phenotypes, including developmental delay/intellectual disability (DD/ID), autism spectrum disorder (ASD), seizures, and hypogonadism. Further, USP7 was identified to critically incorporate into the MAGEL2-USP7-TRIM27 (MUST), such that pathogenic variants in USP7 lead to altered endosomal F-actin polymerization and dysregulated protein recycling.

METHODS

We report 16 newly identified individuals with heterozygous USP7 variants, identified by genome or exome sequencing or by chromosome microarray analysis. Clinical features were evaluated by review of medical records. Additional clinical information was obtained on the seven previously reported individuals to fully elucidate the phenotypic expression associated with USP7 haploinsufficiency.

RESULTS

The clinical manifestations of these 23 individuals suggest a syndrome characterized by DD/ID, hypotonia, eye anomalies,feeding difficulties, GERD, behavioral anomalies, and ASD, and more specific phenotypes of speech delays including a nonverbal phenotype and abnormal brain magnetic resonance image findings including white matter changes based on neuroradiologic examination.

CONCLUSION

The consistency of clinical features among all individuals presented regardless of de novo USP7 variant type supports haploinsufficiency as a mechanism for pathogenesis and refines the clinical impact faced by affected individuals and caregivers.

PSTPIP1-associated myeloid-related proteinemia inflammatory syndrome: A rare cause of childhood neutropenia associated with systemic inflammation and hyperzincemia

Neutropenia in pediatric patients can be due to a variety of disorders. We describe two patients who underwent extensive evaluation over many years for arthralgias and moderate neutropenia of unclear etiology. Genetic testing identified a pathogenic variant in PSTPIP1 (proline-serine-threonine phosphatase-interacting protein 1) in both patients. Markedly elevated inflammatory markers and zinc levels confirmed the rare diagnosis of PSTPIP1-associated myeloid-related proteinemia inflammatory (PAMI) syndrome, tailoring treatment. Neutropenia is common in patients with PAMI syndrome. Unique mutations seen in PAMI syndrome may account for the specific phenotypic features of this disorder.

KCTD7 deficiency defines a distinct neurodegenerative disorder with a conserved autophagy-lysosome defect

OBJECTIVE

Several small case series identified KCTD7 mutations in patients with a rare autosomal recessive disorder designated progressive myoclonic epilepsy (EPM3) and neuronal ceroid lipofuscinosis (CLN14). Despite the name KCTD (potassium channel tetramerization domain), KCTD protein family members lack predicted channel domains. We sought to translate insight gained from yeast studies to uncover disease mechanisms associated with deficiencies in KCTD7 of unknown function.

METHODS

Novel KCTD7 variants in new and published patients were assessed for disease causality using genetic analyses, cell-based functional assays of patient fibroblasts and knockout yeast, and electron microscopy of patient samples.

RESULTS

Patients with KCTD7 mutations can exhibit movement disorders or developmental regression before seizure onset, and are distinguished from similar disorders by an earlier age of onset. Although most published KCTD7 patient variants were excluded from a genome sequence database of normal human variations, most newly identified patient variants are present in this database, potentially challenging disease causality. However, genetic analysis and impaired biochemical interactions with cullin 3 support a causal role for patient KCTD7 variants, suggesting deleterious alleles of KCTD7 and other rare disease variants may be underestimated. Both patient-derived fibroblasts and yeast lacking Whi2 with sequence similarity to KCTD7 have impaired autophagy consistent with brain pathology.

INTERPRETATION

Biallelic KCTD7 mutations define a neurodegenerative disorder with lipofuscin and lipid droplet accumulation but without defining features of neuronal ceroid lipofuscinosis or lysosomal storage disorders. KCTD7 deficiency appears to cause an underlying autophagy-lysosome defect conserved in yeast, thereby assigning a biological role for KCTD7. Ann Neurol 2018;84:774-788.

Amish nemaline myopathy and dilated cardiomyopathy caused by a homozygous contiguous gene deletion of TNNT1 and TNNI3 in a Mennonite child

Amish nemaline myopathy (ANM) is a severe congenital form of NM, known to be fatal in early childhood due to pulmonary insufficiency. Homozygous mutation in TNNT1 was originally ascertained in an Older Amish community in 2000. To date, only five reports with six pathogenic variants in TNNT1 have been described in both Amish and non-Amish families. Here, we describe a 16-month old female from a small Mennonite community from Mexico, presenting with congenital hypotonia and dilated cardiomyopathy, with a novel homozygous deletion of 19q13.42 of about 11 kb in size, encompassing TNNT1 and TNNI3. Cardiomyopathy has not been observed in association with ANM in previous reports. Conversely, homozygous mutation in TNNI3 have been described with dilated cardiomyopathy. Our report underscores the consideration of contiguous gene deletion in children with ANM who present with congenital hypotonia and cardiomyopathy. The report also expands the known spectrum of non-Amish related ANM mutations to include homozygous multi-exonic TNNT1 deletion.