Transcriptome analysis in a humanized mouse model of familial dysautonomia reveals tissue-specific gene expression disruption in the peripheral nervous system.

Familial dysautonomia (FD) is a rare recessive neurodevelopmental disease caused by a splice mutation in the Elongator acetyltransferase complex subunit 1 (ELP1) gene. This mutation results in a tissue-specific reduction of ELP1 protein, with the lowest levels in the central and peripheral nervous systems (CNS and PNS, respectively). FD patients exhibit complex neurological phenotypes due to the loss of sensory and autonomic neurons. Disease symptoms include decreased pain and temperature perception, impaired or absent myotatic reflexes, proprioceptive ataxia, and progressive retinal degeneration. While the involvement of the PNS in FD pathogenesis has been clearly recognized, the underlying mechanisms responsible for the preferential neuronal loss remain unknown. In this study, we aimed to elucidate the molecular mechanisms underlying FD by conducting a comprehensive transcriptome analysis of neuronal tissues from the phenotypic mouse model TgFD9; Elp1Δ20/flox. This mouse recapitulates the same tissue-specific ELP1 mis-splicing observed in patients while modeling many of the disease manifestations. Comparison of FD and control transcriptomes from dorsal root ganglion (DRG), trigeminal ganglion (TG), medulla (MED), cortex, and spinal cord (SC) showed significantly more differentially expressed genes (DEGs) in the PNS than the CNS. We then identified genes that were tightly co-expressed and functionally dependent on the level of full-length ELP1 transcript. These genes, defined as ELP1 dose-responsive genes, were combined with the DEGs to generate tissue-specific dysregulated FD signature genes and networks. Within the PNS networks, we observed direct connections between Elp1 and genes involved in tRNA synthesis and genes related to amine metabolism and synaptic signaling. Importantly, transcriptomic dysregulation in PNS tissues exhibited enrichment for neuronal subtype markers associated with peptidergic nociceptors and myelinated sensory neurons, which are known to be affected in FD. In summary, this study has identified critical tissue-specific gene networks underlying the etiology of FD and provides new insights into the molecular basis of the disease.

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries.

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

Genetic risk variants in New Yorkers of Puerto Rican and Dominican Republic heritage with Parkinson's disease.

There is a paucity of genetic characterization in people with Parkinson's disease (PD) of Latino and Afro-Caribbean descent. Screening LRRK2 and GBA variants in 32 New Yorkers of Puerto Rican ethnicity with PD and in 119 non-Hispanic-non-Jewish European PD cases revealed that Puerto Rican participants were more likely to harbor the LRRK2-p.G2019S variant (15.6% vs. 4.2%, respectively). Additionally, whole exome sequencing of twelve Puerto Rican and Dominican PD participants was performed as an exploratory study.

Transcriptome analysis in a humanized mouse model of familial dysautonomia reveals tissue-specific gene expression disruption in the peripheral nervous system.

Familial dysautonomia (FD) is a rare recessive neurodevelopmental disease caused by a splice mutation in the Elongator acetyltransferase complex subunit 1 ( ELP1 ) gene. This mutation results in a tissue-specific reduction of ELP1 protein, with the lowest levels in the central and peripheral nervous systems (CNS and PNS, respectively). FD patients exhibit complex neurological phenotypes due to the loss of sensory and autonomic neurons. Disease symptoms include decreased pain and temperature perception, impaired or absent myotatic reflexes, proprioceptive ataxia, and progressive retinal degeneration. While the involvement of the PNS in FD pathogenesis has been clearly recognized, the underlying mechanisms responsible for the preferential neuronal loss remain unknown. In this study, we aimed to elucidate the molecular mechanisms underlying FD by conducting a comprehensive transcriptome analysis of neuronal tissues from the phenotypic mouse model TgFD9 ; Elp1 Δ 20/flox . This mouse recapitulates the same tissue-specific ELP1 mis-splicing observed in patients while modeling many of the disease manifestations. Comparison of FD and control transcriptomes from dorsal root ganglion (DRG), trigeminal ganglion (TG), medulla (MED), cortex, and spinal cord (SC) showed significantly more differentially expressed genes (DEGs) in the PNS than the CNS. We then identified genes that were tightly co-expressed and functionally dependent on the level of full-length ELP1 transcript. These genes, defined as ELP1 dose-responsive genes, were combined with the DEGs to generate tissue-specific dysregulated FD signature genes and networks. Within the PNS networks, we observed direct connections between Elp1 and genes involved in tRNA synthesis and genes related to amine metabolism and synaptic signaling. Importantly, transcriptomic dysregulation in PNS tissues exhibited enrichment for neuronal subtype markers associated with peptidergic nociceptors and myelinated sensory neurons, which are known to be affected in FD. In summary, this study has identified critical tissue-specific gene networks underlying the etiology of FD and provides new insights into the molecular basis of the disease.

A marker chromosome in psychosis identifies glycine decarboxylase (GLDC) as a novel regulator of neuronal and synaptic function in the hippocampus.

The biological significance of a small supernumerary marker chromosome that results in dosage alterations to chromosome 9p24.1, including triplication of the GLDC gene encoding glycine decarboxylase, in two patients with psychosis is unclear. In an allelic series of copy number variant mouse models, we identify that triplication of Gldc reduces extracellular glycine levels as determined by optical fluorescence resonance energy transfer (FRET) in dentate gyrus (DG) but not in CA1, suppresses long-term potentiation (LTP) in mPP-DG synapses but not in CA3-CA1 synapses, reduces the activity of biochemical pathways implicated in schizophrenia and mitochondrial bioenergetics, and displays deficits in prepulse inhibition, startle habituation, latent inhibition, working memory, sociability and social preference. Our results thus provide a link between a genomic copy number variation, biochemical, cellular and behavioral phenotypes, and further demonstrate that GLDC negatively regulates long-term synaptic plasticity at specific hippocampal synapses, possibly contributing to the development of neuropsychiatric disorders.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Statistical and functional convergence of common and rare genetic influences on autism at chromosome 16p.

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism's common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models.

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

Tissue-specific and repeat length-dependent somatic instability of the X-linked dystonia parkinsonism-associated CCCTCT repeat.

X-linked dystonia-parkinsonism (XDP) is a progressive adult-onset neurodegenerative disorder caused by insertion of a SINE-VNTR-Alu (SVA) retrotransposon in the TAF1 gene. The SVA retrotransposon contains a CCCTCT hexameric repeat tract of variable length, whose length is inversely correlated with age at onset. This places XDP in a broader class of repeat expansion diseases, characterized by the instability of their causative repeat mutations. Here, we observe similar inverse correlations between CCCTCT repeat length with age at onset and age at death and no obvious correlation with disease duration. To gain insight into repeat instability in XDP we performed comprehensive quantitative analyses of somatic instability of the XDP CCCTCT repeat in blood and in seventeen brain regions from affected males. Our findings reveal repeat length-dependent and expansion-based instability of the XDP CCCTCT repeat, with greater levels of expansion in brain than in blood. The brain exhibits regional-specific patterns of instability that are broadly similar across individuals, with cerebellum exhibiting low instability and cortical regions exhibiting relatively high instability. The spectrum of somatic instability in the brain includes a high proportion of moderate repeat length changes of up to 5 repeats, as well as expansions of ~ 20- > 100 repeats and contractions of ~ 20-40 repeats at lower frequencies. Comparison with HTT CAG repeat instability in postmortem Huntington's disease brains reveals similar brain region-specific profiles, indicating common trans-acting factors that contribute to the instability of both repeats. Analyses in XDP brains of expansion of a different SVA-associated CCCTCT located in the LIPG gene, and not known to be disease-associated, reveals repeat length-dependent expansion at overall lower levels relative to the XDP CCCTCT repeat, suggesting that expansion propensity may be modified by local chromatin structure. Together, the data support a role for repeat length-dependent somatic expansion in the process(es) driving the onset of XDP and prompt further investigation into repeat dynamics and the relationship to disease.

Can Machine Learning Models Predict Asparaginase-associated Pancreatitis in Childhood Acute Lymphoblastic Leukemia.

Asparaginase-associated pancreatitis (AAP) frequently affects children treated for acute lymphoblastic leukemia (ALL) causing severe acute and persisting complications. Known risk factors such as asparaginase dosing, older age and single nucleotide polymorphisms (SNPs) have insufficient odds ratios to allow personalized asparaginase therapy. In this study, we explored machine learning strategies for prediction of individual AAP risk. We integrated information on age, sex, and SNPs based on Illumina Omni2.5exome-8 arrays of patients with childhood ALL (N=1564, 244 with AAP 1.0 to 17.9 yo) from 10 international ALL consortia into machine learning models including regression, random forest, AdaBoost and artificial neural networks. A model with only age and sex had area under the receiver operating characteristic curve (ROC-AUC) of 0.62. Inclusion of 6 pancreatitis candidate gene SNPs or 4 validated pancreatitis SNPs boosted ROC-AUC somewhat (0.67) while 30 SNPs, identified through our AAP genome-wide association study cohort, boosted performance (0.80). Most predictive features included rs10273639 (PRSS1-PRSS2), rs10436957 (CTRC), rs13228878 (PRSS1/PRSS2), rs1505495 (GALNTL6), rs4655107 (EPHB2) and age (1 to 7 y). Second AAP following asparaginase re-exposure was predicted with ROC-AUC: 0.65. The machine learning models assist individual-level risk assessment of AAP for future prevention trials, and may legitimize asparaginase re-exposure when AAP risk is predicted to be low.

Dystonia-specific mutations in THAP1 alter transcription of genes associated with neurodevelopment and myelin.

Dystonia is a neurologic disorder associated with an increasingly large number of genetic variants in many genes, resulting in characteristic disturbances in volitional movement. Dissecting the relationships between these mutations and their functional outcomes is critical in understanding the pathways that drive dystonia pathogenesis. Here we established a pipeline for characterizing an allelic series of dystonia-specific mutations. We used this strategy to investigate the molecular consequences of genetic variation in THAP1, which encodes a transcription factor linked to neural differentiation. Multiple pathogenic mutations associated with dystonia cluster within distinct THAP1 functional domains and are predicted to alter DNA-binding properties and/or protein interactions differently, yet the relative impact of these varied changes on molecular signatures and neural deficits is unclear. To determine the effects of these mutations on THAP1 transcriptional activity, we engineered an allelic series of eight alterations in a common induced pluripotent stem cell background and differentiated these lines into a panel of near-isogenic neural stem cells (n = 94 lines). Transcriptome profiling followed by joint analysis of the most robust signatures across mutations identified a convergent pattern of dysregulated genes functionally related to neurodevelopment, lysosomal lipid metabolism, and myelin. On the basis of these observations, we examined mice bearing Thap1-disruptive alleles and detected significant changes in myelin gene expression and reduction of myelin structural integrity relative to control mice. These results suggest that deficits in neurodevelopment and myelination are common consequences of dystonia-associated THAP1 mutations and highlight the potential role of neuron-glial interactions in the pathogenesis of dystonia.

Physiological Characterization and Transcriptomic Properties of GnRH Neurons Derived From Human Stem Cells.

Gonadotropin-releasing hormone (GnRH) neurons in the hypothalamus play a key role in the regulation of reproductive function. In this study, we sought an efficient method for generating GnRH neurons from human embryonic and induced pluripotent stem cells (hESC and hiPSC, respectively). First, we found that exposure of primitive neuroepithelial cells, rather than neuroprogenitor cells, to fibroblast growth factor 8 (FGF8), was more effective in generating GnRH neurons. Second, addition of kisspeptin to FGF8 further increased the efficiency rates of GnRH neurogeneration. Third, we generated a fluorescent marker mCherry labeled human embryonic GnRH cell line (mCh-hESC) using a CRISPR-Cas9 targeting approach. Fourth, we examined physiological characteristics of GnRH (mCh-hESC) neurons: similar to GnRH neurons in vivo, they released the GnRH peptide in a pulsatile manner at ~60 min intervals; GnRH release increased in response to high potassium, kisspeptin, estradiol, and neurokinin B challenges; and injection of depolarizing current induced action potentials. Finally, we characterized developmental changes in transcriptomes of GnRH neurons using hESC, hiPSC, and mCh-hESC. The developmental pattern of transcriptomes was remarkably similar among the 3 cell lines. Collectively, human stem cell-derived GnRH neurons will be an important tool for establishing disease models to understand diseases, such as idiopathic hypothalamic hypogonadism, and testing contraceptive drugs.

Expanding the Genotypic Spectrum of Congenital Sensory and Autonomic Neuropathies Using Whole-Exome Sequencing.

OBJECTIVE: To test the hypothesis that many patients presenting with congenital insensitivity to pain have lesser known or unidentified mutations not captured by conventional genetic panels, we performed whole-exome sequencing in a cohort of well-characterized patients with a clinical diagnosis of congenital hereditary sensory and autonomic neuropathy with unrevealing conventional genetic testing. METHODS: We performed whole-exome sequencing (WES) in 13 patients with congenital impaired or absent sensation to pain and temperature with no identified molecular diagnosis from a conventional genetic panel. Patients underwent a comprehensive phenotypic assessment including autonomic function testing, and neurologic and ophthalmologic examinations. RESULTS: We identified known or likely pathogenic genetic causes of congenital insensitivity to pain in all 13 patients, spanning 9 genes, the vast majority of which were inherited in an autosomal recessive manner. These included known pathogenic variants (3 patients harboring mutations in TECPR2 and SCN11A), suspected pathogenic variants in genes described to cause congenital sensory and autonomic syndromes (7 patients harboring variants in NGF, LIFR, SCN9A, and PRDM12), and likely pathogenic variants in novel genes (4 patients harboring variants in SMPDL3A, PLEKHN1, and SCN10A). CONCLUSIONS: Our results expand the genetic landscape of congenital sensory and autonomic neuropathies. Further validation of some identified variants should confirm their pathogenicity. WES should be clinically considered to expedite diagnosis, reduce laboratory investigations, and guide enrollment in future gene therapy trials.

DYT-TUBB4A (DYT4 Dystonia): New Clinical and Genetic Observations.

OBJECTIVE: To report 4 novel TUBB4A mutations leading to laryngeal and cervical dystonia with frequent generalization. METHODS: We screened 4 families including a total of 11 definitely affected members with a clinical picture resembling the original description. RESULTS: Four novel variants in the TUBB4A gene have been identified: D295N, R46M, Q424H, and R121W. In silico modeling showed that all variants have characteristics similar to R2G. The variants segregate with the disease in 3 of the families with evidence of incomplete penetrance in 2 of them. All 4 variants would be classified as likely pathogenic. The clinical picture particularly included laryngeal dystonia (often the site of onset), associated with cervical and upper limb dystonia and frequent generalization. Laryngeal dystonia was extremely prevalent (>90%) both in the original cases and in this case series. The hobby horse gait was evident in only 1 patient in this case series. CONCLUSIONS: Our interpretation is that laryngeal involvement is a hallmark feature of DYT-TUBB4A. Nevertheless, TUBB4A mutations remain an exceedingly rare cause of laryngeal or other isolated dystonia.

Isolated dystonia: clinical and genetic updates.

Four genes associated with isolated dystonia are currently well replicated and validated. DYT-THAP1 manifests as young-onset generalized dystonia with predominant craniocervical symptoms; and is associated with mostly deleterious missense variation in the THAP1 gene. De novo and inherited missense and protein truncating variation in GNAL as well as primarily missense variation in ANO3 cause isolated focal and/or segmental dystonia with preference for the upper half of the body and older ages at onset. The GAG deletion in TOR1A is associated with generalized dystonia with onset in childhood in the lower limbs. Rare variation in these genes causes monogenic sporadic and inherited forms of isolated dystonia; common variation may confer risk and imply that dystonia is a polygenic trait in a subset of cases. Although candidate gene screens have been successful in the past in detecting gene-disease associations, recent application of whole-genome and whole-exome sequencing methods enable unbiased capture of all genetic variation that may explain the phenotype. However, careful variant-level evaluation is necessary in every case, even in genes that have previously been associated with disease. We review the genetic architecture and phenotype of DYT-THAP1, DYT-GNAL, DYT-ANO3, and DYT-TOR1A by collecting case reports from the literature and performing variant classification using pathogenicity criteria.

New gene discoveries highlight functional convergence in autism and related neurodevelopmental disorders.

Over the last two years, remarkable gene discovery efforts have implicated disruption of pathways involving gene regulatory functions and neuronal processes in autism spectrum disorder (ASD), and more broadly defined neurodevelopmental disorders (NDDs). Functional studies in the developing brain and across cell types demonstrate that the spatiotemporal expression patterns of many of these genes coalesce on subnetworks with distinct developmental trajectories. Here, we review the convergent biological processes derived from gene discovery and functional genomics in ASD and NDD from 2018-2020. We further probe the mechanistic insights that suggest these frequently perturbed pathways are interconnected and, ultimately, converge on specific functional deficits in human neurodevelopment.

Trypsin-encoding PRSS1-PRSS2 variations influence the risk of asparaginase-associated pancreatitis in children with acute lymphoblastic leukemia: a Ponte di Legno toxicity working group report.

Asparaginase-associated pancreatitis is a life-threatening toxicity to childhood acute lymphoblastic leukemia treatment. To elucidate genetic predisposition and asparaginase-associated pancreatitis pathogenesis, ten trial groups contributed remission samples from patients aged 1.0-17.9 years treated for acute lymphoblastic leukemia between 2000 and 2016. Cases (n=244) were defined by the presence of at least two of the following criteria: (i) abdominal pain; (ii) levels of pancreatic enzymes ≥3 × upper normal limit; and (iii) imaging compatible with pancreatitis. Controls (n=1320) completed intended asparaginase therapy, with 78% receiving ≥8 injections of pegylated-asparaginase, without developing asparaginase-associated pancreatitis. rs62228256 on 20q13.2 showed the strongest association with the development of asparaginase-associated pancreatitis (odds ratio=3.75; P=5.2×10-8). Moreover, rs13228878 (OR=0.61; P=7.1×10-6) and rs10273639 (OR=0.62; P=1.1×10-5) on 7q34 showed significant association with the risk of asparaginase-associated pancreatitis. A Dana Farber Cancer Institute ALL Consortium cohort consisting of patients treated on protocols between 1987 and 2004 (controls=285, cases=33), and the Children's Oncology Group AALL0232 cohort (controls=2653, cases=76) were available as replication cohorts for the 20q13.2 and 7q34 variants, respectively. While rs62228256 was not validated as a risk factor (P=0.77), both rs13228878 (P=0.03) and rs10273639 (P=0.04) were. rs13228878 and rs10273639 are in high linkage disequilibrium (r2=0.94) and associated with elevated expression of the PRSS1 gene, which encodes for trypsinogen, and are known risk variants for alcohol-associated and sporadic pancreatitis in adults. Intra-pancreatic trypsinogen cleavage to proteolytic trypsin induces autodigestion and pancreatitis. In conclusion, this study finds a shared genetic predisposition between asparaginase-associated pancreatitis and non-asparaginase-associated pancreatitis, and targeting the trypsinogen activation pathway may enable identification of effective interventions for asparaginase-associated pancreatitis.