Gene identification for ocular congenital cranial motor neuron disorders using human sequencing, zebrafish screening, and protein binding microarrays.

Purpose: To functionally evaluate novel human sequence-derived candidate genes and variants for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods: Through exome and genome sequencing of a genetically unsolved human oCCDD cohort, we previously identified variants in 80 strong candidate genes. Here, we further prioritized a subset of these (43 human genes, 57 zebrafish genes) using a G0 CRISPR/Cas9-based knockout assay in zebrafish and generated F2 germline mutants for seventeen. We tested the functionality of variants of uncertain significance in known and novel candidate transcription factor-encoding genes through protein binding microarrays. Results: We first demonstrated the feasibility of the G0 screen by targeting known oCCDD genes phox2a and mafba . 70-90% of gene-targeted G0 zebrafish embryos recapitulated germline homozygous null-equivalent phenotypes. Using this approach, we then identified three novel candidate oCCDD genes ( SEMA3F , OLIG2, and FRMD4B ) with putative contributions to human and zebrafish cranial motor development. In addition, protein binding microarrays demonstrated reduced or abolished DNA binding of human variants of uncertain significance in known and novel sequence-derived transcription factors PHOX2A (p.(Trp137Cys)), MAFB (p.(Glu223Lys)), and OLIG2 (p.(Arg156Leu)). Conclusions: This study nominates three strong novel candidate oCCDD genes ( SEMA3F , OLIG2, and FRMD4B ) and supports the functionality and putative pathogenicity of transcription factor candidate variants PHOX2A p.(Trp137Cys), MAFB p.(Glu223Lys), and OLIG2 p.(Arg156Leu). Our findings support that G0 loss-of-function screening in zebrafish can be coupled with human sequence analysis and protein binding microarrays to aid in prioritizing oCCDD candidate genes/variants.

Protein-extending ACTN2 frameshift variants cause variable myopathy phenotypes by protein aggregation.

OBJECTIVE: The objective of the study is to characterize the pathomechanisms underlying actininopathies. Distal myopathies are a group of rare, inherited muscular disorders characterized by progressive loss of muscle fibers that begin in the distal parts of arms and legs. Recently, variants in a new disease gene, ACTN2, have been shown to cause distal myopathy. ACTN2, a gene previously only associated with cardiomyopathies, encodes alpha-actinin-2, a protein expressed in both cardiac and skeletal sarcomeres. The primary function of alpha-actinin-2 is to link actin and titin to the sarcomere Z-disk. New ACTN2 variants are continuously discovered; however, the clinical significance of many variants remains unknown. Thus, lack of clear genotype-phenotype correlations in ACTN2-related diseases, actininopathies, persists. METHODS: Functional characterization in C2C12 cell model of several ACTN2 variants is conducted, including frameshift and missense variants associated with dominant and recessive actininopathies. We assess the genotype-phenotype correlations of actininopathies using clinical data from several patients carrying these variants. RESULTS: The results show that the missense variants associated with a recessive form of actininopathy do not cause detectable alpha-actinin-2 aggregates in the cell model. Conversely, dominant frameshift variants causing a protein extension do form alpha-actinin-2 aggregates. INTERPRETATION: The results suggest that alpha-actinin-2 aggregation is the disease mechanism underlying some dominant actininopathies, and thus, we recommend that protein-extending frameshift variants in ACTN2 should be classified as pathogenic. However, this mechanism is likely elicited by only a limited number of variants. Alternative functional characterization methods should be explored to further investigate other molecular mechanisms underlying actininopathies.

Early B-cell transcription factor-2 defect as a novel cause of lipodystrophy: disruption of the adipose tissue character and integrity.

We report a novel cause of partial lipodystrophy associated with early B cell factor 2 (EBF2) nonsense variant (EBF2 8:26033143 C>A, c.493G>T, p.E165X) in a patient with an atypical form of partial lipodystrophy. The patient presented with progressive adipose tissue loss and metabolic deterioration at pre-pubertal age. In vitro and in vivo disease modeling demonstrates that the EBF2 variant impairs adipogenesis, causing excess accumulation of undifferentiated CD34+ cells, extracellular matrix proteins, and inflammatory myeloid cells in subcutaneous adipose tissues. Thus, this EBF2 p.E165X variant disrupts adipose tissue structure and function, leading to the development of partial lipodystrophy syndrome.

CIAO1 loss of function causes a neuromuscular disorder with compromise of nucleocytoplasmic Fe-S enzymes.

Cytoplasmic and nuclear iron-sulfur (Fe-S) enzymes that are essential for genome maintenance and replication depend on the cytoplasmic Fe-S assembly (CIA) machinery for cluster acquisition. The core of the CIA machinery consists of a complex of CIAO1, MMS19 and FAM96B. The physiological consequences of loss of function in the components of the CIA pathway have thus far remained uncharacterized. Our study revealed that patients with biallelic loss of function in CIAO1 developed proximal and axial muscle weakness, fluctuating creatine kinase elevation, and respiratory insufficiency. In addition, they presented with CNS symptoms including learning difficulties and neurobehavioral comorbidities, along with iron deposition in deep brain nuclei, mild normocytic to macrocytic anemia, and gastrointestinal symptoms. Mutational analysis revealed reduced stability of the variants compared with WT CIAO1. Functional assays demonstrated failure of the variants identified in patients to recruit Fe-S recipient proteins, resulting in compromised activities of DNA helicases, polymerases, and repair enzymes that rely on the CIA complex to acquire their Fe-S cofactors. Lentivirus-mediated restoration of CIAO1 expression reversed all patient-derived cellular abnormalities. Our study identifies CIAO1 as a human disease gene and provides insights into the broader implications of the cytosolic Fe-S assembly pathway in human health and disease.

Pathogenic variants in KMT2C result in a neurodevelopmental disorder distinct from Kleefstra and Kabuki syndromes.

Trithorax-related H3K4 methyltransferases, KMT2C and KMT2D, are critical epigenetic modifiers. Haploinsufficiency of KMT2C was only recently recognized as a cause of neurodevelopmental disorder (NDD), so the clinical and molecular spectrums of the KMT2C-related NDD (now designated as Kleefstra syndrome 2) are largely unknown. We ascertained 98 individuals with rare KMT2C variants, including 75 with protein-truncating variants (PTVs). Notably, ∼15% of KMT2C PTVs were inherited. Although the most highly expressed KMT2C transcript consists of only the last four exons, pathogenic PTVs were found in almost all the exons of this large gene. KMT2C variant interpretation can be challenging due to segmental duplications and clonal hematopoesis-induced artifacts. Using samples from 27 affected individuals, divided into discovery and validation cohorts, we generated a moderate strength disorder-specific KMT2C DNA methylation (DNAm) signature and demonstrate its utility in classifying non-truncating variants. Based on 81 individuals with pathogenic/likely pathogenic variants, we demonstrate that the KMT2C-related NDD is characterized by developmental delay, intellectual disability, behavioral and psychiatric problems, hypotonia, seizures, short stature, and other comorbidities. The facial module of PhenoScore, applied to photographs of 34 affected individuals, reveals that the KMT2C-related facial gestalt is significantly different from the general NDD population. Finally, using PhenoScore and DNAm signatures, we demonstrate that the KMT2C-related NDD is clinically and epigenetically distinct from Kleefstra and Kabuki syndromes. Overall, we define the clinical features, molecular spectrum, and DNAm signature of the KMT2C-related NDD and demonstrate they are distinct from Kleefstra and Kabuki syndromes highlighting the need to rename this condition.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

PURPOSE: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). METHODS: We coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. RESULTS: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. CONCLUSION: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

Dominant stop-loss HNRNPA1 variants in juvenile-onset myopathy.

INTRODUCTION/AIMS: Heterogeneous nuclear ribonucleoprotein A1 is involved in nucleic acid homeostatic functions. The encoding gene HNRNPA1 has been associated with several neuromuscular disorders including an amyotrophic lateral sclerosis-like phenotype, distal hereditary motor neuropathy, multisystem proteinopathy, and various myopathies. We report two unrelated individuals with monoallelic stop loss variants affecting the same codon of HNRNPA1. METHODS: Two individuals with unsolved juvenile-onset myopathy were enrolled under approved institutional protocols. Phenotype data were collected and genetic analyses were performed, including whole-exome sequencing (WES). RESULTS: The two probands (MNOT002-01 and K1440-01) showed a similar onset of slowly progressive extremity and facial weakness in early adolescence. K1440-01 presented with facial weakness, winged scapula, elevated serum creatine kinase (CK) levels, and mild neck weakness. MNOT002-01 also exhibited elevated CK levels along with facial weakness, cardiomyopathy, respiratory dysfunction, pectus excavatum, a mildly rigid spine, and loss of ambulation. On quadriceps muscle biopsy, K1440-01 displayed rounded myofibers, mild variation in fiber diameter, and type 2 fiber hypertrophy, while MNOT002-01 displayed rimmed vacuoles. Monoallelic stop-loss variants in HNRNPA1 were identified for both probands: c.1119A>C p.*373Tyrext*6 (K1440-01) and c.1118A>C p.*373Serext*6 (MNOT002-01) affect the same codon and are both predicted to lead to the addition of six amino acids before termination at an alternative stop codon. DISCUSSION: Both stop-loss variants in our probands are likely pathogenic. Our findings contribute to the disease characterization of pathogenic variants in HNRNPA1. This gene should be screened in clinical diagnostic testing of unsolved cases of sporadic or dominant juvenile-onset myopathy.

Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

De novo TLK1 and MDM1 mutations in a patient with a neurodevelopmental disorder and immunodeficiency.

The Tousled-like kinases 1 and 2 (TLK1/TLK2) regulate DNA replication, repair and chromatin maintenance. TLK2 variants underlie the neurodevelopmental disorder (NDD) 'Intellectual Disability, Autosomal Dominant 57' (MRD57), characterized by intellectual disability and microcephaly. Several TLK1 variants have been reported in NDDs but their functional significance is unknown. A male patient presenting with ID, seizures, global developmental delay, hypothyroidism, and primary immunodeficiency was determined to have a heterozygous TLK1 variant (c.1435C>G, p.Q479E), as well as a mutation in MDM1 (c.1197dupT, p.K400∗). Cells expressing TLK1 p.Q479E exhibited reduced cytokine responses and elevated DNA damage, but not increased radiation sensitivity or DNA repair defects. The TLK1 p.Q479E variant impaired kinase activity but not proximal protein interactions. Our study provides the first functional characterization of NDD-associated TLK1 variants and suggests that, such as TLK2, TLK1 variants may impact development in multiple tissues and should be considered in the diagnosis of rare NDDs.

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Mono and biallelic variants in HCN2 cause severe neurodevelopmental disorders.

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

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

Purpose: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods: We coupled phenotyping with exome or genome sequencing of 467 pedigrees with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. Results: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. Conclusion: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

A metabolic signature for NADSYN1-dependent congenital NAD deficiency disorder.

Nicotinamide adenine dinucleotide (NAD) is essential for embryonic development. To date, biallelic loss-of-function variants in 3 genes encoding nonredundant enzymes of the NAD de novo synthesis pathway - KYNU, HAAO, and NADSYN1 - have been identified in humans with congenital malformations defined as congenital NAD deficiency disorder (CNDD). Here, we identified 13 further individuals with biallelic NADSYN1 variants predicted to be damaging, and phenotypes ranging from multiple severe malformations to the complete absence of malformation. Enzymatic assessment of variant deleteriousness in vitro revealed protein domain-specific perturbation, complemented by protein structure modeling in silico. We reproduced NADSYN1-dependent CNDD in mice and assessed various maternal NAD precursor supplementation strategies to prevent adverse pregnancy outcomes. While for Nadsyn1+/- mothers, any B3 vitamer was suitable to raise NAD, preventing embryo loss and malformation, Nadsyn1-/- mothers required supplementation with amidated NAD precursors (nicotinamide or nicotinamide mononucleotide) bypassing their metabolic block. The circulatory NAD metabolome in mice and humans before and after NAD precursor supplementation revealed a consistent metabolic signature with utility for patient identification. Our data collectively improve clinical diagnostics of NADSYN1-dependent CNDD, provide guidance for the therapeutic prevention of CNDD, and suggest an ongoing need to maintain NAD levels via amidated NAD precursor supplementation after birth.

Diagnosing missed cases of spinal muscular atrophy in genome, exome, and panel sequencing datasets.

Spinal muscular atrophy (SMA) is a genetic disorder that causes progressive degeneration of lower motor neurons and the subsequent loss of muscle function throughout the body. It is the second most common recessive disorder in individuals of European descent and is present in all populations. Accurate tools exist for diagnosing SMA from genome sequencing data. However, there are no publicly available tools for GRCh38-aligned data from panel or exome sequencing assays which continue to be used as first line tests for neuromuscular disorders. This deficiency creates a critical gap in our ability to diagnose SMA in large existing rare disease cohorts, as well as newly sequenced exome and panel datasets. We therefore developed and extensively validated a new tool - SMA Finder - that can diagnose SMA not only in genome, but also exome and panel sequencing samples aligned to GRCh37, GRCh38, or T2T-CHM13. It works by evaluating aligned reads that overlap the c.840 position of SMN1 and SMN2 in order to detect the most common molecular causes of SMA. We applied SMA Finder to 16,626 exomes and 3,911 genomes from heterogeneous rare disease cohorts sequenced at the Broad Institute Center for Mendelian Genomics as well as 1,157 exomes and 8,762 panel sequencing samples from Tartu University Hospital. SMA Finder correctly identified all 16 known SMA cases and reported nine novel diagnoses which have since been confirmed by clinical testing, with another four novel diagnoses undergoing validation. Notably, out of the 29 total SMA positive cases, 23 had an initial clinical diagnosis of muscular dystrophy, congenital myasthenic syndrome, or myopathy. This underscored the frequency with which SMA can be misdiagnosed as other neuromuscular disorders and confirmed the utility of using SMA Finder to reanalyze phenotypically diverse neuromuscular disease cohorts. Finally, we evaluated SMA Finder on 198,868 individuals that had both exome and genome sequencing data within the UK Biobank (UKBB) and found that SMA Finder's overall false positive rate was less than 1 / 200,000 exome samples, and its positive predictive value (PPV) was 97%. We also observed 100% concordance between UKBB exome and genome calls. This analysis showed that, even though it is located within a segmental duplication, the most common causal variant for SMA can be detected with comparable accuracy to monogenic disease variants in non-repetitive regions. Additionally, the high PPV demonstrated by SMA Finder, the existence of treatment options for SMA in which early diagnosis is imperative for therapeutic benefit, as well as widespread availability of clinical confirmatory testing for SMA, warrants the addition of SMN1 to the ACMG list of genes with reportable secondary findings after genome and exome sequencing.

Rare ACTN2 Frameshift Variants Resulting in Protein Extension Cause Distal Myopathy and Hypertrophic Cardiomyopathy through Protein Aggregation.

Distal myopathies are a group of rare, inherited muscular disorders characterized by progressive loss of muscle fibers that begins in the distal parts of arms and legs. Recently, variants in a new disease gene, ACTN2 , have been shown to cause distal myopathy. ACTN2 , a gene previously only associated with cardiomyopathies, encodes alpha-actinin-2, a protein expressed in both cardiac and skeletal sarcomeres. The primary function of alpha-actinin-2 is to link actin and titin to the sarcomere Z-disk. New ACTN2 variants are continuously discovered, however, the clinical significance of many variants remains unknown. Thus, lack of clear genotype-phenotype correlations in ACTN2 -related diseases, actininopathies, persists. Objective: The objective of the study is to characterize the pathomechanisms underlying actininopathies. Methods: Functional characterization in C2C12 cell models of several ACTN2 variants is conducted, including frameshift and missense variants associated with dominant actininopathies. We assess the genotype-phenotype correlations of actininopathies using clinical data from several patients carrying these variants. Results: The results show that the missense variants associated with a recessive form of actininopathy do not cause detectable alpha-actinin-2 aggregates in the cell model. Conversely, dominant frameshift variants causing a protein extension do produce alpha-actinin-2 aggregates. Interpretation: The results suggest that alpha-actinin-2 aggregation is the disease mechanism underlying some dominant actininopathies, and thus we recommend that protein-extending frameshift variants in ACTN2 should be classified as pathogenic. However, this mechanism is likely elicited by only a limited number of variants. Alternative functional characterization methods should be explored to further investigate other molecular mechanisms underlying actininopathies.

Erythromelalgia caused by the missense mutation p.Arg220Pro in an alternatively spliced exon of SCN9A (NaV1.7).

Erythromelalgia (EM), is a familial pain syndrome characterized by episodic 'burning' pain, warmth, and erythema. EM is caused by monoallelic variants in SCN9A, which encodes the voltage-gated sodium channel (NaV) NaV1.7. Over 25 different SCN9A mutations attributed to EM have been described to date, all identified in the SCN9A transcript utilizing exon 6N. Here we report a novel SCN9A missense variant identified in seven related individuals with stereotypic episodes of bilateral lower limb pain presenting in childhood. The variant, XM_011511617.3:c.659G>C;p.(Arg220Pro), resides in the exon 6A of SCN9A, an exon previously shown to be selectively incorporated by developmentally regulated alternative splicing. The mutation is located in the voltage-sensing S4 segment of domain I, which is important for regulating channel activation. Functional analysis showed the p.Arg220Pro mutation altered voltage-dependent activation and delayed channel inactivation, consistent with a NaV1.7 gain-of-function molecular phenotype. These results demonstrate that alternatively spliced isoforms of SCN9A should be included in all genomic testing of EM.

Exome copy number variant detection, analysis and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and with new innovative methods can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the GREGoR consortium. Each family's CNV data was analyzed using the seqr platform and candidate CNVs classified using the 2020 ACMG/ClinGen CNV interpretation standards. We developed additional evidence criteria to address situations not covered by the current standards. The addition of CNV calling to exome analysis identified causal CNVs for 173 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb with estimates that 44% would not have been detected by standard chromosomal microarrays. The causal CNVs consisted of 141 deletions, 15 duplications, 4 suspected complex structural variants (SVs), 3 insertions and 10 complex SVs, the latter two groups being identified by orthogonal validation methods. We interpreted 153 CNVs as likely pathogenic/pathogenic and 20 CNVs as high interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Identification of a de novo mutation in TLK1 associated with a neurodevelopmental disorder and immunodeficiency.

Background: The Tousled-like kinases 1 and 2 (TLK1/TLK2) regulate DNA replication, repair and chromatin maintenance. TLK2 variants are associated with 'Intellectual Disability, Autosomal Dominant 57' (MRD57), a neurodevelopmental disorder (NDD) characterized by intellectual disability (ID), autism spectrum disorder (ASD) and microcephaly. Several TLK1 variants have been reported in NDDs but their functional significance is unknown. Methods: A male patient presenting with ID, seizures, global developmental delay, hypothyroidism, and primary immunodeficiency was determined to have a novel, heterozygous variant in TLK1 (c.1435C>G, p.Q479E) by genome sequencing (GS). Single cell gel electrophoresis, western blot, flow cytometry and RNA-seq were performed in patient-derived lymphoblast cell lines. In silico, biochemical and proteomic analysis were used to determine the functional impact of the p.Q479E variant and previously reported NDD-associated TLK1 variant, p.M566T. Results: Transcriptome sequencing in patient-derived cells confirmed expression of TLK1 transcripts carrying the p.Q479E variant and revealed alterations in genes involved in class switch recombination and cytokine signaling. Cells expressing the p.Q479E variant exhibited reduced cytokine responses and higher levels of spontaneous DNA damage but not increased sensitivity to radiation or DNA repair defects. The p.Q479E and p.M566T variants impaired kinase activity but did not strongly alter localization or proximal protein interactions. Conclusion: Our study provides the first functional characterization of TLK1 variants associated with NDDs and suggests potential involvement in central nervous system and immune system development. Our results indicate that, like TLK2 variants, TLK1 variants may impact development in multiple tissues and should be considered in the diagnosis of rare NDDs.