Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian 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.

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.

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.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

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.

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.

The genome of a Late Pleistocene human from a Clovis burial site in western Montana.

Clovis, with its distinctive biface, blade and osseous technologies, is the oldest widespread archaeological complex defined in North America, dating from 11,100 to 10,700 (14)C years before present (bp) (13,000 to 12,600 calendar years bp). Nearly 50 years of archaeological research point to the Clovis complex as having developed south of the North American ice sheets from an ancestral technology. However, both the origins and the genetic legacy of the people who manufactured Clovis tools remain under debate. It is generally believed that these people ultimately derived from Asia and were directly related to contemporary Native Americans. An alternative, Solutrean, hypothesis posits that the Clovis predecessors emigrated from southwestern Europe during the Last Glacial Maximum. Here we report the genome sequence of a male infant (Anzick-1) recovered from the Anzick burial site in western Montana. The human bones date to 10,705 ± 35 (14)C years bp (approximately 12,707-12,556 calendar years bp) and were directly associated with Clovis tools. We sequenced the genome to an average depth of 14.4× and show that the gene flow from the Siberian Upper Palaeolithic Mal'ta population into Native American ancestors is also shared by the Anzick-1 individual and thus happened before 12,600 years bp. We also show that the Anzick-1 individual is more closely related to all indigenous American populations than to any other group. Our data are compatible with the hypothesis that Anzick-1 belonged to a population directly ancestral to many contemporary Native Americans. Finally, we find evidence of a deep divergence in Native American populations that predates the Anzick-1 individual.

Disease onset in X-linked dystonia-parkinsonism correlates with expansion of a hexameric repeat within an SVA retrotransposon in TAF1.

X-linked dystonia-parkinsonism (XDP) is a neurodegenerative disease associated with an antisense insertion of a SINE-VNTR-Alu (SVA)-type retrotransposon within an intron of TAF1 This unique insertion coincides with six additional noncoding sequence changes in TAF1, the gene that encodes TATA-binding protein-associated factor-1, which appear to be inherited together as an identical haplotype in all reported cases. Here we examined the sequence of this SVA in XDP patients (n = 140) and detected polymorphic variation in the length of a hexanucleotide repeat domain, (CCCTCT)n The number of repeats in these cases ranged from 35 to 52 and showed a highly significant inverse correlation with age at disease onset. Because other SVAs exhibit intrinsic promoter activity that depends in part on the hexameric domain, we assayed the transcriptional regulatory effects of varying hexameric lengths found in the unique XDP SVA retrotransposon using luciferase reporter constructs. When inserted sense or antisense to the luciferase reading frame, the XDP variants repressed or enhanced transcription, respectively, to an extent that appeared to vary with length of the hexamer. Further in silico analysis of this SVA sequence revealed multiple motifs predicted to form G-quadruplexes, with the greatest potential detected for the hexameric repeat domain. These data directly link sequence variation within the XDP-specific SVA sequence to phenotypic variability in clinical disease manifestation and provide insight into potential mechanisms by which this intronic retroelement may induce transcriptional interference in TAF1 expression.

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.

Global gene expression profiling of brown to white adipose tissue transformation in sheep reveals novel transcriptional components linked to adipose remodeling.

BACKGROUND: Large mammals are capable of thermoregulation shortly after birth due to the presence of brown adipose tissue (BAT). The majority of BAT disappears after birth and is replaced by white adipose tissue (WAT). RESULTS: We analyzed the postnatal transformation of adipose in sheep with a time course study of the perirenal adipose depot. We observed changes in tissue morphology, gene expression and metabolism within the first two weeks of postnatal life consistent with the expected transition from BAT to WAT. The transformation was characterized by massively decreased mitochondrial abundance and down-regulation of gene expression related to mitochondrial function and oxidative phosphorylation. Global gene expression profiling demonstrated that the time points grouped into three phases: a brown adipose phase, a transition phase and a white adipose phase. Between the brown adipose and the transition phase 170 genes were differentially expressed, and 717 genes were differentially expressed between the transition and the white adipose phase. Thirty-eight genes were shared among the two sets of differentially expressed genes. We identified a number of regulated transcription factors, including NR1H3, MYC, KLF4, ESR1, RELA and BCL6, which were linked to the overall changes in gene expression during the adipose tissue remodeling. Finally, the perirenal adipose tissue expressed both brown and brite/beige adipocyte marker genes at birth, the expression of which changed substantially over time. CONCLUSIONS: Using global gene expression profiling of the postnatal BAT to WAT transformation in sheep, we provide novel insight into adipose tissue plasticity in a large mammal, including identification of novel transcriptional components linked to adipose tissue remodeling. Moreover, our data set provides a useful resource for further studies in adipose tissue plasticity.

Transcriptome profiling of brown adipose tissue during cold exposure reveals extensive regulation of glucose metabolism.

We applied digital gene expression profiling to determine the transcriptome of brown and white adipose tissues (BAT and WAT, respectively) during cold exposure. Male C57BL/6J mice were exposed to cold for 2 or 4 days. A notable induction of genes related to glucose uptake, glycolysis, glycogen metabolism, and the pentose phosphate pathway was observed in BAT from cold-exposed animals. In addition, glycerol-3-phosphate dehydrogenase 1 expression was induced in BAT from cold-challenged mice, suggesting increased synthesis of glycerol from glucose. Similarly, expression of lactate dehydrogenases was induced by cold in BAT. Pyruvate dehydrogenase kinase 2 (Pdk2) and Pdk4 were expressed at significantly higher levels in BAT than in WAT, and Pdk2 was induced in BAT by cold. Of notice, only a subset of the changes detected in BAT was observed in WAT. Based on changes in gene expression during cold exposure, we propose a model for the intermediary glucose metabolism in activated BAT: 1) fluxes through glycolysis and the pentose phosphate pathway are induced, the latter providing reducing equivalents for de novo fatty acid synthesis; 2) glycerol synthesis from glucose is increased, facilitating triacylglycerol synthesis/fatty acid re-esterification; 3) glycogen turnover and lactate production are increased; and 4) entry of glucose carbon into the tricarboxylic acid cycle is restricted by PDK2 and PDK4. In summary, our results demonstrate extensive and diverse gene expression changes related to glucose handling in activated BAT.

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.

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.

TIMP-1 increases expression and phosphorylation of proteins associated with drug resistance in breast cancer cells.

Tissue inhibitor of metalloproteinase 1 (TIMP-1) is a protein with a potential biological role in drug resistance. To elucidate the unknown molecular mechanisms underlying the association between high TIMP-1 levels and increased chemotherapy resistance, we employed SILAC-based quantitative mass spectrometry to analyze global proteome and phosphoproteome differences of MCF-7 breast cancer cells expressing high or low levels of TIMP-1. In TIMP-1 high expressing cells, 312 proteins and 452 phosphorylation sites were up-regulated. Among these were the cancer drug targets topoisomerase 1, 2A, and 2B, which may explain the resistance phenotype to topoisomerase inhibitors that was observed in cells with high TIMP-1 levels. Pathway analysis showed an enrichment of proteins from functional categories such as apoptosis, cell cycle, DNA repair, transcription factors, drug targets and proteins associated with drug resistance or sensitivity, and drug transportation. The NetworKIN algorithm predicted the protein kinases CK2a, CDK1, PLK1, and ATM as likely candidates involved in the hyperphosphorylation of the topoisomerases. Up-regulation of protein and/or phosphorylation levels of topoisomerases in TIMP-1 high expressing cells may be part of the mechanisms by which TIMP-1 confers resistance to treatment with the widely used topoisomerase inhibitors in breast and colorectal cancer.

TIMP1 overexpression mediates resistance of MCF-7 human breast cancer cells to fulvestrant and down-regulates progesterone receptor expression.

High levels of Tissue Inhibitor of Metalloproteinases-1 (TIMP1) are associated with poor prognosis, reduced response to chemotherapy, and, potentially, also poor response to endocrine therapy in breast cancer patients. Our objective was to further investigate the hypothesis that TIMP1 is associated with endocrine sensitivity. We established a panel of 11 MCF-7 subclones with a wide range of TIMP1 mRNA and protein expression levels. Cells with high expression of TIMP1 versus low TIMP1 displayed significantly reduced sensitivity to the antiestrogen fulvestrant (ICI 182,780, Faslodex®), while TIMP1 levels did not influence the sensitivity to 4-hydroxytamoxifen. An inverse correlation between expression of the progesterone receptor and TIMP1 was found, but TIMP1 levels did not correlate with estrogen receptor levels or growth-promoting effects of estrogen (estradiol, E2). Additionally, the effects of fulvestrant, 4-hydroxytamoxifen, or estrogen on estrogen receptor expression were not associated with TIMP1 levels. Gene expression analyses revealed associations between expression of TIMP1 and genes involved in metabolic pathways, epidermal growth factor receptor 1/cancer signaling pathways, and cell cycle. Gene and protein expression analyses showed no general defects in estrogen receptor signaling except from lack of progesterone receptor expression and estrogen inducibility in clones with high TIMP1. The present study suggests a relation between high expression level of TIMP1 and loss of progesterone receptor expression combined with fulvestrant resistance. Our findings in vitro may have clinical implications as the data suggest that high tumor levels of TIMP1 may be a predictive biomarker for reduced response to fulvestrant.

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.

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.