Loss-of-function in RBBP5 results in a syndromic neurodevelopmental disorder associated with microcephaly.

PURPOSE: Epigenetic dysregulation has been associated with many inherited disorders. RBBP5 (HGNC:9888) encodes a core member of the protein complex that methylates histone 3 lysine-4 (H3K4) and has not been implicated in human disease. METHODS: We identify five unrelated individuals with de novo heterozygous variants in RBBP5. Three nonsense/frameshift and two missense variants were identified in probands with neurodevelopmental symptoms including global developmental delay, intellectual disability, microcephaly, and short stature. Here, we investigate the pathogenicity of the variants through protein structural analysis and transgenic Drosophila models. RESULTS: Both missense p.(T232I) and p.(E296D) variants affect evolutionarily conserved amino acids located at the interface between RBBP5 and the nucleosome. In Drosophila, overexpression analysis identifies partial loss-of-function mechanisms when the variants are expressed using the fly Rbbp5 or human RBBP5 cDNA. Loss of Rbbp5 leads to a reduction in brain size. The human reference or variant transgenes fail to rescue this loss and expression of either missense variant in an Rbbp5 null background results in a less severe microcephaly phenotype than the human reference, indicating both missense variants are partial loss-of-function alleles. CONCLUSION: Haploinsufficiency of RBBP5 observed through de novo null and hypomorphic loss-of-function variants is associated with a syndromic neurodevelopmental disorder.

Reanalysis of RNA sequencing data ends diagnostic odyssey and expands the phenotypic spectrum of congenital titinopathy.

Although next-generation sequencing has enabled diagnoses for many patients with Mendelian disorders, the majority remain undiagnosed. Here, we present a sibling pair who were clinically diagnosed with Escobar syndrome, however targeted gene testing was negative. Exome sequencing (ES), and later genome sequencing (GS), revealed compound heterozygous TTN variants in both siblings, a maternally inherited frameshift variant [(NM_133378.4):c.36812del; p.(Asp12271Valfs*10)], and a paternally inherited missense variant [(NM_133378.4):c.12322G > A; p.(Asp4108Asn)]. This result was considered nondiagnostic due to poor clinical fit and limited pathogenicity evidence for the missense variant of uncertain significance (VUS). Following initial nondiagnostic RNA sequencing (RNAseq) on muscle and further pursuit of other variants detected on the ES/GS, a reanalysis of noncanonical splice sites in the muscle transcriptome identified an out-of-frame exon retraction in TTN, near the known VUS. Interim literature included reports of patients with similar TTN variants who had phenotypic concordance with the siblings, and a diagnosis of a congenital titinopathy was given 4 years after the TTN variants had been initially reported. This report highlights the value of reanalysis of RNAseq with a different approach, expands the phenotypic spectrum of congenital titinopathy and also illustrates how a perceived phenotypic mismatch, and failure to consider known variants, can result in a prolongation of the diagnostic journey.

Disseminated Coccidioidomycosis Treated with Interferon-γ and Dupilumab.

We describe a case of life-threatening disseminated coccidioidomycosis in a previously healthy child. Like most patients with disseminated coccidioidomycosis, this child had no genomic evidence of any known, rare immune disease. However, comprehensive immunologic testing showed exaggerated production of interleukin-4 and reduced production of interferon-γ. Supplementation of antifungal agents with interferon-γ treatment slowed disease progression, and the addition of interleukin-4 and interleukin-13 blockade with dupilumab resulted in rapid resolution of the patient's clinical symptoms. This report shows that blocking of type 2 immune responses can treat infection. This immunomodulatory approach could be used to enhance immune clearance of refractory fungal, mycobacterial, and viral infections. (Supported by the Jeffrey Modell Foundation and the National Institutes of Health.).

Mosaic de novo SNRPN gene variant associated with Prader-Willi syndrome.

BACKGROUND: Prader-Willi syndrome (PWS) is an imprinting disorder caused by the absence of paternal expressed genes in the Prader-Willi critical region (PWCR) on chromosome 15q11.2-q13. Three molecular mechanisms have been known to cause PWS, including a deletion in the PWCR, uniparental disomy 15 and imprinting defects. RESULTS: We report the first case of PWS associated with a single-nucleotide SNRPN variant in a 10-year-old girl presenting with clinical features consistent with PWS, including infantile hypotonia and feeding difficulty, developmental delay with cognitive impairment, excessive eating with central obesity, sleep disturbances, skin picking and related behaviour issues. Whole-exome sequencing revealed a de novo mosaic nonsense variant of the SNRPN gene (c.73C>T, p.R25X) in 10% of DNA isolated from buccal cells and 19% of DNA from patient-derived lymphoblast cells. DNA methylation study did not detect an abnormal methylation pattern in the SNRPN locus. Parental origin studies showed a paternal source of an intronic single-nucleotide polymorphism within the locus in proximity to the SNRPN variant. CONCLUSIONS: This is the first report that provides evidence of a de novo point mutation of paternal origin in SNRPN as a new disease-causing mechanism for PWS. This finding suggests that gene sequencing should be considered as part of the diagnostic workup in patients with clinical suspicion of PWS.

Genetic and functional analysis of a Pacific hagfish opioid system.

The actions of endogenous opioids and nociceptin/orphanin FQ are mediated by four homologous G protein-coupled receptors that constitute the opioid receptor family. However, little is known about opioid systems in cyclostomes (living jawless fish) and how opioid systems might have evolved from invertebrates. Here, we leveraged de novo transcriptome and low-coverage whole-genome assembly in the Pacific hagfish (Eptatretus stoutii) to identify and characterize the first full-length coding sequence for a functional opioid receptor in a cyclostome. Additionally, we define two novel endogenous opioid precursors in this species that predict several novel opioid peptides. Bioinformatic analysis shows no closely related opioid receptor genes in invertebrates with regard either to the genomic organization or to conserved opioid receptor-specific sequences that are common in all vertebrates. Furthermore, no proteins analogous to vertebrate opioid precursors could be identified by genomic searches despite previous claims of protein or RNA-derived sequences in several invertebrate species. The presence of an expressed orthologous receptor and opioid precursors in the Pacific hagfish confirms that a functional opioid system was likely present in the common ancestor of all extant vertebrates some 550 million years ago, earlier than all previous authenticated accounts. We discuss the premise that the cyclostome and vertebrate opioid systems evolved from invertebrate systems concerned with antimicrobial defense and speculate that the high concentrations of opioid precursors in tissues such as the testes, gut, and activated immune cells are key remnants of this evolutionary role.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases.

PURPOSE: Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful. METHODS: We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols. RESULTS: We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases. CONCLUSION: The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

GATAD2B-associated neurodevelopmental disorder (GAND): clinical and molecular insights into a NuRD-related disorder.

PURPOSE: Determination of genotypic/phenotypic features of GATAD2B-associated neurodevelopmental disorder (GAND). METHODS: Fifty GAND subjects were evaluated to determine consistent genotypic/phenotypic features. Immunoprecipitation assays utilizing in vitro transcription-translation products were used to evaluate GATAD2B missense variants' ability to interact with binding partners within the nucleosome remodeling and deacetylase (NuRD) complex. RESULTS: Subjects had clinical findings that included macrocephaly, hypotonia, intellectual disability, neonatal feeding issues, polyhydramnios, apraxia of speech, epilepsy, and bicuspid aortic valves. Forty-one novelGATAD2B variants were identified with multiple variant types (nonsense, truncating frameshift, splice-site variants, deletions, and missense). Seven subjects were identified with missense variants that localized within two conserved region domains (CR1 or CR2) of the GATAD2B protein. Immunoprecipitation assays revealed several of these missense variants disrupted GATAD2B interactions with its NuRD complex binding partners. CONCLUSIONS: A consistent GAND phenotype was caused by a range of genetic variants in GATAD2B that include loss-of-function and missense subtypes. Missense variants were present in conserved region domains that disrupted assembly of NuRD complex proteins. GAND's clinical phenotype had substantial clinical overlap with other disorders associated with the NuRD complex that involve CHD3 and CHD4, with clinical features of hypotonia, intellectual disability, cardiac defects, childhood apraxia of speech, and macrocephaly.

Diagnostic utility of transcriptome sequencing for rare Mendelian diseases.

PURPOSE: We investigated the value of transcriptome sequencing (RNAseq) in ascertaining the consequence of DNA variants on RNA transcripts to improve the diagnostic rate from exome or genome sequencing for undiagnosed Mendelian diseases spanning a wide spectrum of clinical indications. METHODS: From 234 subjects referred to the Undiagnosed Diseases Network, University of California-Los Angeles clinical site between July 2014 and August 2018, 113 were enrolled for high likelihood of having rare undiagnosed, suspected genetic conditions despite thorough prior clinical evaluation. Exome or genome sequencing and RNAseq were performed, and RNAseq data was integrated with genome sequencing data for DNA variant interpretation genome-wide. RESULTS: The molecular diagnostic rate by exome or genome sequencing was 31%. Integration of RNAseq with genome sequencing resulted in an additional seven cases with clear diagnosis of a known genetic disease. Thus, the overall molecular diagnostic rate was 38%, and 18% of all genetic diagnoses returned required RNAseq to determine variant causality. CONCLUSION: In this rare disease cohort with a wide spectrum of undiagnosed, suspected genetic conditions, RNAseq analysis increased the molecular diagnostic rate above that possible with genome sequencing analysis alone even without availability of the most appropriate tissue type to assess.

Correction: GATAD2B-associated neurodevelopmental disorder (GAND): clinical and molecular insights into a NuRD-related disorder.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genetic screening of a large series of North American sporadic and familial frontotemporal dementia cases.

INTRODUCTION: The Advancing Research and Treatment for Frontotemporal Lobar Degeneration (ARTFL) and Longitudinal Evaluation of Familial Frontotemporal Dementia Subjects (LEFFTDS) consortia are two closely connected studies, involving multiple North American centers that evaluate both sporadic and familial frontotemporal dementia (FTD) participants and study longitudinal changes. METHODS: We screened the major dementia-associated genes in 302 sporadic and 390 familial (symptomatic or at-risk) participants enrolled in these studies. RESULTS: Among the sporadic patients, 16 (5.3%) carried chromosome 9 open reading frame 72 (C9orf72), microtubule-associated protein tau (MAPT), and progranulin (GRN) pathogenic variants, whereas in the familial series we identified 207 carriers from 146 families. Of interest, one patient was found to carry a homozygous C9orf72 expansion, while another carried both a C9orf72 expansion and a GRN pathogenic variant. We also identified likely pathogenic variants in the TAR DNA binding protein (TARDBP), presenilin 1 (PSEN1), and valosin containing protein (VCP) genes, and a subset of variants of unknown significance in other rare FTD genes. DISCUSSION: Our study reports the genetic characterization of a large FTD series and supports an unbiased sequencing screen, irrespective of clinical presentation or family history.

Genome-wide survey of copy number variants finds MAPT duplications in progressive supranuclear palsy.

BACKGROUND: Progressive supranuclear palsy is a neurodegenerative tauopathy manifesting clinically as a progressive akinetic-rigid syndrome. In this study, we sought to identify genetic variants influencing PSP susceptibility through a genome-wide association analysis of a cohort of well-characterized patients who had participated in the Neuroprotection and Natural History in Parkinson Plus Syndromes and Blood Brain Barrier in Parkinson Plus Syndromes studies. METHODS: We genotyped single-nucleotide polymorphisms in 283 PSP cases from the United Kingdom, Germany, and France and compared these with genotypes from 4472 controls. Copy number variants were identified from genotyping data. RESULTS: We observed associations on chromosome 17 within or close to the MAPT gene and explored the genetic architecture at this locus. We confirmed the previously reported association of rs1768208 in the MOBP gene (P = 3.29 × 10-13 ) and rs1411478 in STX6 (P = 3.45 × 10-10 ). The population-attributable risk from the MAPT, MOBP, and STX6 single-nucleotide polymorphisms was found to be 0.37, 0.26, and 0.08, respectively. In addition, we found 2 instances of copy number variants spanning the MAPT gene in patients with PSP. These copy number variants include tau but few other genes within the chromosome 17 haplotype region, providing additional support for the direct pathogenicity of MAPT in PSP. CONCLUSIONS: Clinicians should also be aware of MAPT duplication as a possible genetic cause of PSP, especially in patients presenting with young age at onset. © 2019 International Parkinson and Movement Disorder Society.

Genetic screen in a large series of patients with primary progressive aphasia.

INTRODUCTION: Primary progressive aphasia (PPA) is a neurological syndrome, associated with both frontotemporal dementia and Alzheimer's disease, in which progressive language impairment emerges as the most salient clinical feature during the initial stages of disease. METHODS: We screened the main genes associated with Alzheimer's disease and frontotemporal dementia for pathogenic and risk variants in a cohort of 403 PPA cases. RESULTS: In this case series study, 14 (3.5%) cases carried (likely) pathogenic variants: four C9orf72 expansions, nine GRN, and one TARDBP mutation. Rare risk variants, TREM2 R47H and MAPT A152T, were associated with a three- to seven-fold increase in risk for PPA. DISCUSSION: Our results show that while pathogenic variants within the most common dementia genes were rarely associated with PPA, these were found almost exclusively in GRN and C9orf72, suggesting that PPA is more TDP43- than tau-related in our series. This is consistent with the finding that PPA frequency in dominantly inherited dementias is the highest in kindreds with GRN variants.

Evidence for a role of the rare p.A152T variant in MAPT in increasing the risk for FTD-spectrum and Alzheimer's diseases.

Rare mutations in the gene encoding for tau (MAPT, microtubule-associated protein tau) cause frontotemporal dementia-spectrum (FTD-s) disorders, including FTD, progressive supranuclear palsy (PSP) and corticobasal syndrome, and a common extended haplotype spanning across the MAPT locus is associated with increased risk of PSP and Parkinson's disease. We identified a rare tau variant (p.A152T) in a patient with a clinical diagnosis of PSP and assessed its frequency in multiple independent series of patients with neurodegenerative conditions and controls, in a total of 15 369 subjects. Tau p.A152T significantly increases the risk for both FTD-s (n = 2139, OR = 3.0, CI: 1.6-5.6, P = 0.0005) and Alzheimer's disease (AD) (n = 3345, OR = 2.3, CI: 1.3-4.2, P = 0.004) compared with 9047 controls. Functionally, p.A152T (i) decreases the binding of tau to microtubules and therefore promotes microtubule assembly less efficiently; and (ii) reduces the tendency to form abnormal fibers. However, there is a pronounced increase in the formation of tau oligomers. Importantly, these findings suggest that other regions of the tau protein may be crucial in regulating normal function, as the p.A152 residue is distal to the domains considered responsible for microtubule interactions or aggregation. These data provide both the first genetic evidence and functional studies supporting the role of MAPT p.A152T as a rare risk factor for both FTD-s and AD and the concept that rare variants can increase the risk for relatively common, complex neurodegenerative diseases, but since no clear significance threshold for rare genetic variation has been established, some caution is warranted until the findings are further replicated.

Mutations in SLC20A2 are a major cause of familial idiopathic basal ganglia calcification.

Familial idiopathic basal ganglia calcification (IBGC) or Fahr's disease is a rare neurodegenerative disorder characterized by calcium deposits in the basal ganglia and other brain regions, which is associated with neuropsychiatric and motor symptoms. Familial IBGC is genetically heterogeneous and typically transmitted in an autosomal dominant fashion. We performed a mutational analysis of SLC20A2, the first gene found to cause IBGC, to assess its genetic contribution to familial IBGC. We recruited 218 subjects from 29 IBGC-affected families of varied ancestry and collected medical history, neurological exam, and head CT scans to characterize each patient's disease status. We screened our patient cohort for mutations in SLC20A2. Twelve novel (nonsense, deletions, missense, and splice site) potentially pathogenic variants, one synonymous variant, and one previously reported mutation were identified in 13 families. Variants predicted to be deleterious cosegregated with disease in five families. Three families showed nonsegregation with clinical disease of such variants, but retrospective review of clinical and neuroimaging data strongly suggested previous misclassification. Overall, mutations in SLC20A2 account for as many as 41% of our familial IBGC cases. Our screen in a large series expands the catalog of SLC20A2 mutations identified to date and demonstrates that mutations in SLC20A2 are a major cause of familial IBGC. Non-perfect segregation patterns of predicted deleterious variants highlight the challenges of phenotypic assessment in this condition with highly variable clinical presentation.

Familial cortical myoclonus with a mutation in NOL3.

OBJECTIVE: Myoclonus is characterized by sudden, brief involuntary movements, and its presence is debilitating. We identified a family suffering from adult onset, cortical myoclonus without associated seizures. We performed clinical, electrophysiological, and genetic studies to define this phenotype. METHODS: A large, 4-generation family with a history of myoclonus underwent careful questioning, examination, and electrophysiological testing. Thirty-five family members donated blood samples for genetic analysis, which included single nucleotide polymorphism mapping, microsatellite linkage, targeted massively parallel sequencing, and Sanger sequencing. In silico and in vitro experiments were performed to investigate functional significance of the mutation. RESULTS: We identified 11 members of a Canadian Mennonite family suffering from adult onset, slowly progressive, disabling, multifocal myoclonus. Somatosensory evoked potentials indicated a cortical origin of the myoclonus. There were no associated seizures. Some severely affected individuals developed signs of progressive cerebellar ataxia of variable severity late in the course of their illness. The phenotype was inherited in an autosomal dominant fashion. We demonstrated linkage to chromosome 16q21-22.1. We then sequenced all coding sequence in the critical region, identifying only a single cosegregating, novel, nonsynonymous mutation, which resides in the gene NOL3. Furthermore, this mutation was found to alter post-translational modification of NOL3 protein in vitro. INTERPRETATION: We propose that familial cortical myoclonus is a novel movement disorder that may be caused by mutation in NOL3. Further investigation of the role of NOL3 in neuronal physiology may shed light on neuronal membrane hyperexcitability and pathophysiology of myoclonus and related disorders.

Whole-Exome Sequencing Identifies Homozygote Nonsense Variants in LMOD2 Gene Causing Infantile Dilated Cardiomyopathy.

As an essential component of the sarcomere, actin thin filament stems from the Z-disk extend toward the middle of the sarcomere and overlaps with myosin thick filaments. Elongation of the cardiac thin filament is essential for normal sarcomere maturation and heart function. This process is regulated by the actin-binding proteins Leiomodins (LMODs), among which LMOD2 has recently been identified as a key regulator of thin filament elongation to reach a mature length. Few reports have implicated homozygous loss of function variants of LMOD2 in neonatal dilated cardiomyopathy (DCM) associated with thin filament shortening. We present the fifth case of DCM due to biallelic variants in the LMOD2 gene and the second case with the c.1193G>A (p.W398*) nonsense variant identified by whole-exome sequencing. The proband is a 4-month male infant of Hispanic descent with advanced heart failure. Consistent with previous reports, a myocardial biopsy exhibited remarkably short thin filaments. However, compared to other cases of identical or similar biallelic variants, the patient presented here has an unusually late onset of cardiomyopathy during infancy. Herein, we present the phenotypic and histological features of this variant, confirm the pathogenic impact on protein expression and sarcomere structure, and discuss the current knowledge of LMOD2-related cardiomyopathy.

Full-length isoform sequencing for resolving the molecular basis of Charcot-Marie-Tooth 2A.

Objectives: Transcript sequencing of patient derived samples has been shown to improve the diagnostic yield for solving cases of likely Mendelian disorders, yet the added benefit of full-length long-read transcript sequencing is largely unexplored. Methods: We applied short-read and full-length isoform cDNA sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis. Results: We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts a branch point critical for intron 6 spicing. Full-length long-read isoform cDNA sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates five distinct altered splicing transcripts. All five altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 protein levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A). Discussion: This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Full-length Isoform Sequencing for Resolving the Molecular Basis of Charcot-Marie-Tooth 2A.

Objectives: Transcript sequencing of patient-derived samples has been shown to improve the diagnostic yield for solving cases of suspected Mendelian conditions, yet the added benefit of full-length long-read transcript sequencing is largely unexplored. Methods: We applied short-read and full-length transcript sequencing and mitochondrial functional studies to a patient-derived fibroblast cell line from an individual with neuropathy that previously lacked a molecular diagnosis. Results: We identified an intronic homozygous MFN2 c.600-31T>G variant that disrupts the branch point critical for intron 6 splicing. Full-length long-read isoform complementary DNA (cDNA) sequencing after treatment with a nonsense-mediated mRNA decay (NMD) inhibitor revealed that this variant creates 5 distinct altered splicing transcripts. All 5 altered splicing transcripts have disrupted open reading frames and are subject to NMD. Furthermore, a patient-derived fibroblast line demonstrated abnormal lipid droplet formation, consistent with MFN2 dysfunction. Although correctly spliced full-length MFN2 transcripts are still produced, this branch point variant results in deficient MFN2 levels and autosomal recessive Charcot-Marie-Tooth disease, axonal, type 2A (CMT2A). Discussion: This case highlights the utility of full-length isoform sequencing for characterizing the molecular mechanism of undiagnosed rare diseases and expands our understanding of the genetic basis for CMT2A.

Mutations in rare ataxia genes are uncommon causes of sporadic cerebellar ataxia.

BACKGROUND: Sporadic-onset ataxia is common in a tertiary care setting but a significant percentage remains unidentified despite extensive evaluation. Rare genetic ataxias, reported only in specific populations or families, may contribute to a percentage of sporadic ataxia. METHODS: Patients with adult-onset sporadic ataxia, who tested negative for common genetic ataxias (SCA1, SCA2, SCA3, SCA6, SCA7, and/or Friedreich ataxia), were evaluated using a stratified screening approach for variants in 7 rare ataxia genes. RESULTS: We screened patients for published mutations in SYNE1 (n = 80) and TGM6 (n = 118), copy number variations in LMNB1 (n = 40) and SETX (n = 11), sequence variants in SACS (n = 39) and PDYN (n = 119), and the pentanucleotide insertion of spinocerebellar ataxia type 31 (n = 101). Overall, we identified 1 patient with a LMNB1 duplication, 1 patient with a PDYN variant, and 1 compound SACS heterozygote, including a novel variant. CONCLUSIONS: The rare genetic ataxias examined here do not significantly contribute to sporadic cerebellar ataxia in our tertiary care population.

Synaptic processes and immune-related pathways implicated in Tourette syndrome.

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.