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.

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.

A form of muscular dystrophy associated with pathogenic variants in JAG2.

JAG2 encodes the Notch ligand Jagged2. The conserved Notch signaling pathway contributes to the development and homeostasis of multiple tissues, including skeletal muscle. We studied an international cohort of 23 individuals with genetically unsolved muscular dystrophy from 13 unrelated families. Whole-exome sequencing identified rare homozygous or compound heterozygous JAG2 variants in all 13 families. The identified bi-allelic variants include 10 missense variants that disrupt highly conserved amino acids, a nonsense variant, two frameshift variants, an in-frame deletion, and a microdeletion encompassing JAG2. Onset of muscle weakness occurred from infancy to young adulthood. Serum creatine kinase (CK) levels were normal or mildly elevated. Muscle histology was primarily dystrophic. MRI of the lower extremities revealed a distinct, slightly asymmetric pattern of muscle involvement with cores of preserved and affected muscles in quadriceps and tibialis anterior, in some cases resembling patterns seen in POGLUT1-associated muscular dystrophy. Transcriptome analysis of muscle tissue from two participants suggested misregulation of genes involved in myogenesis, including PAX7. In complementary studies, Jag2 downregulation in murine myoblasts led to downregulation of multiple components of the Notch pathway, including Megf10. Investigations in Drosophila suggested an interaction between Serrate and Drpr, the fly orthologs of JAG1/JAG2 and MEGF10, respectively. In silico analysis predicted that many Jagged2 missense variants are associated with structural changes and protein misfolding. In summary, we describe a muscular dystrophy associated with pathogenic variants in JAG2 and evidence suggests a disease mechanism related to Notch pathway dysfunction.

Variants in DTNA cause a mild, dominantly inherited muscular dystrophy.

DTNA encodes α-dystrobrevin, a component of the macromolecular dystrophin-glycoprotein complex (DGC) that binds to dystrophin/utrophin and α-syntrophin. Mice lacking α-dystrobrevin have a muscular dystrophy phenotype, but variants in DTNA have not previously been associated with human skeletal muscle disease. We present 12 individuals from four unrelated families with two different monoallelic DTNA variants affecting the coiled-coil domain of α-dystrobrevin. The five affected individuals from family A harbor a c.1585G > A; p.Glu529Lys variant, while the recurrent c.1567_1587del; p.Gln523_Glu529del DTNA variant was identified in the other three families (family B: four affected individuals, family C: one affected individual, and family D: two affected individuals). Myalgia and exercise intolerance, with variable ages of onset, were reported in 10 of 12 affected individuals. Proximal lower limb weakness with onset in the first decade of life was noted in three individuals. Persistent elevations of serum creatine kinase (CK) levels were detected in 11 of 12 affected individuals, 1 of whom had an episode of rhabdomyolysis at 20 years of age. Autism spectrum disorder or learning disabilities were reported in four individuals with the c.1567_1587 deletion. Muscle biopsies in eight affected individuals showed mixed myopathic and dystrophic findings, characterized by fiber size variability, internalized nuclei, and slightly increased extracellular connective tissue and inflammation. Immunofluorescence analysis of biopsies from five affected individuals showed reduced α-dystrobrevin immunoreactivity and variably reduced immunoreactivity of other DGC proteins: dystrophin, α, ß, δ and γ-sarcoglycans, and α and ß-dystroglycans. The DTNA deletion disrupted an interaction between α-dystrobrevin and syntrophin. Specific variants in the coiled-coil domain of DTNA cause skeletal muscle disease with variable penetrance. Affected individuals show a spectrum of clinical manifestations, with severity ranging from hyperCKemia, myalgias, and exercise intolerance to childhood-onset proximal muscle weakness. Our findings expand the molecular etiologies of both muscular dystrophy and paucisymptomatic hyperCKemia, to now include monoallelic DTNA variants as a novel cause of skeletal muscle disease in humans.

Fatal Case of Chronic Jamestown Canyon Virus Encephalitis Diagnosed by Metagenomic Sequencing in Patient Receiving Rituximab.

A 56-year-old man receiving rituximab who had months of neurologic symptoms was found to have Jamestown Canyon virus in cerebrospinal fluid by clinical metagenomic sequencing. The patient died, and postmortem examination revealed extensive neuropathologic abnormalities. Deep sequencing enabled detailed characterization of viral genomes from the cerebrospinal fluid, cerebellum, and cerebral cortex.

Rapid Detection of Powassan Virus in a Patient With Encephalitis by Metagenomic Sequencing.

We describe a patient with severe and progressive encephalitis of unknown etiology. We performed rapid metagenomic sequencing from cerebrospinal fluid and identified Powassan virus, an emerging tick-borne flavivirus that has been increasingly detected in the United States.

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.

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.

Expert Panel Curation of 31 Genes in Relation to Limb Girdle Muscular Dystrophy.

Introduction: Limb girdle muscular dystrophies (LGMDs) are a group of genetically heterogeneous autosomal conditions with some degree of phenotypic homogeneity. LGMD is defined as having onset >2 years of age with progressive proximal weakness, elevated serum creatine kinase levels and dystrophic features on muscle biopsy. Advances in massively parallel sequencing have led to a surge in genes linked to LGMD. Methods: The ClinGen Muscular Dystrophies and Myopathies gene curation expert panel (MDM GCEP, formerly Limb Girdle Muscular Dystrophy GCEP) convened to evaluate the strength of evidence supporting gene-disease relationships (GDR) using the ClinGen gene-disease clinical validity framework to evaluate 31 genes implicated in LGMD. Results: The GDR was exclusively LGMD for 17 genes, whereas an additional 14 genes were related to a broader phenotype encompassing congenital weakness. Four genes (CAPN3, COL6A1, COL6A2, COL6A3) were split into two separate disease entities, based on each displaying both dominant and recessive inheritance patterns, resulting in curation of 35 GDRs. Of these, 30 (86%) were classified as Definitive, 4 (11%) as Moderate and 1 (3%) as Limited. Two genes, POMGNT1 and DAG1, though definitively related to myopathy, currently have insufficient evidence to support a relationship specifically with LGMD. Conclusions: The expert-reviewed assertions on the clinical validity of genes implicated in LGMDs form an invaluable resource for clinicians and molecular geneticists. We encourage the global neuromuscular community to publish case-level data that help clarify disputed or novel LGMD associations.

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.

Novel syndromic neurodevelopmental disorder caused by de novo deletion of CHASERR, a long noncoding RNA.

Genes encoding long non-coding RNAs (lncRNAs) comprise a large fraction of the human genome, yet haploinsufficiency of a lncRNA has not been shown to cause a Mendelian disease. CHASERR is a highly conserved human lncRNA adjacent to CHD2-a coding gene in which de novo loss-of-function variants cause developmental and epileptic encephalopathy. Here we report three unrelated individuals each harboring an ultra-rare heterozygous de novo deletion in the CHASERR locus. We report similarities in severe developmental delay, facial dysmorphisms, and cerebral dysmyelination in these individuals, distinguishing them from the phenotypic spectrum of CHD2 haploinsufficiency. We demonstrate reduced CHASERR mRNA expression and corresponding increased CHD2 mRNA and protein in whole blood and patient-derived cell lines-specifically increased expression of the CHD2 allele in cis with the CHASERR deletion, as predicted from a prior mouse model of Chaserr haploinsufficiency. We show for the first time that de novo structural variants facilitated by Alu-mediated non-allelic homologous recombination led to deletion of a non-coding element (the lncRNA CHASERR) to cause a rare syndromic neurodevelopmental disorder. We also demonstrate that CHD2 has bidirectional dosage sensitivity in human disease. This work highlights the need to carefully evaluate other lncRNAs, particularly those upstream of genes associated with Mendelian disorders.

Genome and RNA sequencing boost neuromuscular diagnoses to 62% from 34% with exome sequencing alone.

OBJECTIVE: Most families with heritable neuromuscular disorders do not receive a molecular diagnosis. Here we evaluate diagnostic utility of exome, genome, RNA sequencing, and protein studies and provide evidence-based recommendations for their integration into practice. METHODS: In total, 247 families with suspected monogenic neuromuscular disorders who remained without a genetic diagnosis after standard diagnostic investigations underwent research-led massively parallel sequencing: neuromuscular disorder gene panel, exome, genome, and/or RNA sequencing to identify causal variants. Protein and RNA studies were also deployed when required. RESULTS: Integration of exome sequencing and auxiliary genome, RNA and/or protein studies identified causal or likely causal variants in 62% (152 out of 247) of families. Exome sequencing alone informed 55% (83 out of 152) of diagnoses, with remaining diagnoses (45%; 69 out of 152) requiring genome sequencing, RNA and/or protein studies to identify variants and/or support pathogenicity. Arrestingly, novel disease genes accounted for <4% (6 out of 152) of diagnoses while 36.2% of solved families (55 out of 152) harbored at least one splice-altering or structural variant in a known neuromuscular disorder gene. We posit that contemporary neuromuscular disorder gene-panel sequencing could likely provide 66% (100 out of 152) of our diagnoses today. INTERPRETATION: Our results emphasize thorough clinical phenotyping to enable deep scrutiny of all rare genetic variation in phenotypically consistent genes. Post-exome auxiliary investigations extended our diagnostic yield by 81% overall (34-62%). We present a diagnostic algorithm that details deployment of genomic and auxiliary investigations to obtain these diagnoses today most effectively. We hope this provides a practical guide for clinicians as they gain greater access to clinical genome and transcriptome sequencing.

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

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

A homozygous mutation in the tight-junction protein JAM3 causes hemorrhagic destruction of the brain, subependymal calcification, and congenital cataracts.

The tight junction, or zonula occludens, is a specialized cell-cell junction that regulates epithelial and endothelial permeability, and it is an essential component of the blood-brain barrier in the cerebrovascular endothelium. In addition to functioning as a diffusion barrier, tight junctions are also involved in signal transduction. In this study, we identified a homozygous mutation in the tight-junction protein gene JAM3 in a large consanguineous family from the United Arab Emirates. Some members of this family had a rare autosomal-recessive syndrome characterized by severe hemorrhagic destruction of the brain, subependymal calcification, and congenital cataracts. Their clinical presentation overlaps with some reported cases of pseudo-TORCH syndrome as well as with cases involving mutations in occludin, another component of the tight-junction complex. However, massive intracranial hemorrhage distinguishes these patients from others. Homozygosity mapping identified the disease locus in this family on chromosome 11q25 with a maximum multipoint LOD score of 6.15. Sequence analysis of genes in the candidate interval uncovered a mutation in the canonical splice-donor site of intron 5 of JAM3. RT-PCR analysis of a patient lymphoblast cell line confirmed abnormal splicing, leading to a frameshift mutation with early termination. JAM3 is known to be present in vascular endothelium, although its roles in cerebral vasculature have not been implicated. Our results suggest that JAM3 is essential for maintaining the integrity of the cerebrovascular endothelium as well as for normal lens development in humans.

Mutations in ARFGEF2 implicate vesicle trafficking in neural progenitor proliferation and migration in the human cerebral cortex.

Disruption of human neural precursor proliferation can give rise to a small brain (microcephaly), and failure of neurons to migrate properly can lead to an abnormal arrest of cerebral cortical neurons in proliferative zones near the lateral ventricles (periventricular heterotopia). Here we show that an autosomal recessive condition characterized by microcephaly and periventricular heterotopia maps to chromosome 20 and is caused by mutations in the gene ADP-ribosylation factor guanine nucleotide-exchange factor-2 (ARFGEF2). By northern-blot analysis, we found that mouse Arfgef2 mRNA levels are highest during embryonic periods of ongoing neuronal proliferation and migration, and by in situ hybridization, we found that the mRNA is widely distributed throughout the embryonic central nervous system (CNS). ARFGEF2 encodes the large (>200 kDa) brefeldin A (BFA)-inhibited GEF2 protein (BIG2), which is required for vesicle and membrane trafficking from the trans-Golgi network (TGN). Inhibition of BIG2 by BFA, or by a dominant negative ARFGEF2 cDNA, decreases cell proliferation in vitro, suggesting a cell-autonomous regulation of neural expansion. Inhibition of BIG2 also disturbed the intracellular localization of such molecules as E-cadherin and beta-catenin by preventing their transport from the Golgi apparatus to the cell surface. Our findings show that vesicle trafficking is an important regulator of proliferation and migration during human cerebral cortical development.

Transcriptome and Genome Analysis Uncovers a DMD Structural Variant: A Case Report.

Objective: Duchenne muscular dystrophy (DMD) is caused by pathogenic variants in the dystrophin gene (DMD). Hypermethylated CGG expansions within DIP2B 5' UTR are associated with an intellectual development disorder. Here, we demonstrate the diagnostic utility of genomic short-read sequencing (SRS) and transcriptome sequencing to identify a novel DMD structural variant (SV) and a DIP2B CGG expansion in a patient with DMD for whom conventional diagnostic testing failed to yield a genetic diagnosis. Methods: We performed genomic SRS, skeletal muscle transcriptome sequencing, and targeted programmable long-read sequencing (LRS). Results: The proband had a typical DMD clinical presentation, autism spectrum disorder (ASD), and dystrophinopathy on muscle biopsy. Transcriptome analysis identified 6 aberrantly expressed genes; DMD and DIP2B were the strongest underexpression and overexpression outliers, respectively. Genomic SRS identified a 216 kb paracentric inversion (NC_000023.11: g.33162217-33378800) overlapping 2 DMD promoters. ExpansionHunter indicated an expansion of 109 CGG repeats within the 5' UTR of DIP2B. Targeted genomic LRS confirmed the SV and genotyped the DIP2B repeat expansion as 270 CGG repeats. Discussion: Here, transcriptome data heavily guided genomic analysis to resolve a complex DMD inversion and a DIP2B repeat expansion. Longitudinal follow-up will be important for clarifying the clinical significance of the DIP2B genotype.

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.

Exome Sequencing and the Identification of New Genes and Shared Mechanisms in Polymicrogyria.

Importance: Polymicrogyria is the most commonly diagnosed cortical malformation and is associated with neurodevelopmental sequelae including epilepsy, motor abnormalities, and cognitive deficits. Polymicrogyria frequently co-occurs with other brain malformations or as part of syndromic diseases. Past studies of polymicrogyria have defined heterogeneous genetic and nongenetic causes but have explained only a small fraction of cases. Objective: To survey germline genetic causes of polymicrogyria in a large cohort and to consider novel polymicrogyria gene associations. Design, Setting, and Participants: This genetic association study analyzed panel sequencing and exome sequencing of accrued DNA samples from a retrospective cohort of families with members with polymicrogyria. Samples were accrued over more than 20 years (1994 to 2020), and sequencing occurred in 2 stages: panel sequencing (June 2015 to January 2016) and whole-exome sequencing (September 2019 to March 2020). Individuals seen at multiple clinical sites for neurological complaints found to have polymicrogyria on neuroimaging, then referred to the research team by evaluating clinicians, were included in the study. Targeted next-generation sequencing and/or exome sequencing were performed on probands (and available parents and siblings) from 284 families with individuals who had isolated polymicrogyria or polymicrogyria as part of a clinical syndrome and no genetic diagnosis at time of referral from clinic, with sequencing from 275 families passing quality control. Main Outcomes and Measures: The number of families in whom genetic sequencing yielded a molecular diagnosis that explained the polymicrogyria in the family. Secondarily, the relative frequency of different genetic causes of polymicrogyria and whether specific genetic causes were associated with co-occurring head size changes were also analyzed. Results: In 32.7% (90 of 275) of polymicrogyria-affected families, genetic variants were identified that provided satisfactory molecular explanations. Known genes most frequently implicated by polymicrogyria-associated variants in this cohort were PIK3R2, TUBB2B, COL4A1, and SCN3A. Six candidate novel polymicrogyria genes were identified or confirmed: de novo missense variants in PANX1, QRICH1, and SCN2A and compound heterozygous variants in TMEM161B, KIF26A, and MAN2C1, each with consistent genotype-phenotype relationships in multiple families. Conclusions and Relevance: This study's findings reveal a higher than previously recognized rate of identifiable genetic causes, specifically of channelopathies, in individuals with polymicrogyria and support the utility of exome sequencing for families affected with polymicrogyria.