Search

Glycine is an obligatory co-agonist at excitatory NMDA receptors in the brain, especially in the dentate gyrus, which has been postulated to be crucial for the development of psychotic associations and memories with psychotic content. Drugs modulating glycine levels are in clinical development for improving cognition in schizophrenia. However, the functional relevance of the regulation of glycine metabolism by endogenous enzymes is unclear. Using a chromosome-engineered allelic series in mice, we report that a triplication of the gene encoding the glycine-catabolizing enzyme glycine decarboxylase (GLDC) - as found on a small supernumerary marker chromosome in patients with psychosis - reduces extracellular glycine levels as determined by optical fluorescence resonance energy transfer (FRET) in dentate gyrus (DG) and 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 schizophrenia-like behaviors which are in part known to be dependent on the activity of the dentate gyrus, e.g., prepulse inhibition, startle habituation, latent inhibition, working memory, sociability and social preference. Our results demonstrate that Gldc negatively regulates long-term synaptic plasticity in the dentate gyrus in mice, suggesting that an increase in GLDC copy number possibly contributes to the development of psychosis in humans., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2024. The Author(s).)

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., (© 2023. The Author(s).)

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., Competing Interests: Declaration of interests M.E.T. receives research funding and/or reagents from Levo Therapeutics, Microsoft Inc, and Illumina Inc. K.M. initiated employment with Tornado Bio during the writing of this manuscript and is at present an employee of the company. J.F.G. is a founder and member of the scientific advisory board of Triplet Therapeutics, Inc., and has been a paid consultant to Biogen, Inc., Pfizer, Inc., and Wave Biosciences, Inc., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

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., (© 2022. The Author(s).)

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., Competing Interests: Declaration of interests J.F.G. is a founder and member of the scientific advisory board of Triplet Therapeutics and has been a paid consultant to Biogen, Pfizer, and Wave Biosciences. M.E.T. receives research funding and/or reagents from Illumina, Levo Therapeutics, and Microsoft. E.M. receives research funding and reagents from PTC Therapeutics. The following co-authors are currently employed by for-profit companies or non-profit organizations: P.R. is founder and CEO of Kumuda; K.M. is currently employed by Tornado Bio; T.A. is a founder and CEO of Independent Data Lab and OmicsChart; and A.H. is currently employed by Bionano Genomics., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

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., (© 2022. The Author(s).)

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., Competing Interests: The authors declare no conflict of interest., (Copyright © 2021 The Author(s). Published by Wolters Kluwer Health, Inc.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Endocrine Society.)

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., (© 2020 American Academy of Neurology.)

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.

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., (Copyright © 2021 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.)

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., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

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 (r 2 =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., (Copyright© 2019 Ferrata Storti Foundation.)

Genomic association studies of common or rare protein-coding variation have established robust statistical approaches to account for multiple testing. Here we present a comparable framework to evaluate rare and de novo noncoding single-nucleotide variants, insertion/deletions, and all classes of structural variation from whole-genome sequencing (WGS). Integrating genomic annotations at the level of nucleotides, genes, and regulatory regions, we define 51,801 annotation categories. Analyses of 519 autism spectrum disorder families did not identify association with any categories after correction for 4,123 effective tests. Without appropriate correction, biologically plausible associations are observed in both cases and controls. Despite excluding previously identified gene-disrupting mutations, coding regions still exhibited the strongest associations. Thus, in autism, the contribution of de novo noncoding variation is probably modest in comparison to that of de novo coding variants. Robust results from future WGS studies will require large cohorts and comprehensive analytical strategies that consider the substantial multiple-testing burden.

Root-knot nematodes (RKNs, Meloidogyne incognita) are economically important endoparasites with a wide host range. We used a comprehensive transcriptomic approach to investigate the expression of both tomato and RKN genes in tomato roots at five infection time intervals from susceptible plants and two infection time intervals from resistant plants, grown under soil conditions. Differentially expressed genes during susceptible (1827, tomato; 462, RKN) and resistance (25, tomato; 160, RKN) interactions were identified. In susceptible responses, tomato genes involved in cell wall structure, development, primary and secondary metabolite, and defence signalling pathways, together with RKN genes involved in host parasitism, development and defence, are discussed. In resistance responses, tomato genes involved in secondary metabolite and hormone-mediated defence responses, together with RKN genes involved in starvation stress-induced apoptosis, are discussed. In addition, 40 novel differentially expressed RKN genes encoding secretory proteins were identified. Our findings provide novel insights into the temporal regulation of genes involved in various biological processes from tomato and RKN simultaneously during susceptible and resistance responses, and reveal the involvement of a complex network of biosynthetic pathways during disease development., (© 2017 BSPP AND JOHN WILEY & SONS LTD.)

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., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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., Competing Interests: The authors declare no conflict of interest., (Copyright © 2017 the Author(s). Published by PNAS.)

Gene disruption by clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) is highly efficient and relies on the error-prone non-homologous end-joining pathway. Conversely, precise gene editing requires homology-directed repair (HDR), which occurs at a lower frequency than non-homologous end-joining in mammalian cells. Here, by testing whether manipulation of DNA repair factors improves HDR efficacy, we show that transient ectopic co-expression of RAD52 and a dominant-negative form of tumour protein p53-binding protein 1 (dn53BP1) synergize to enable efficient HDR using a single-stranded oligonucleotide DNA donor template at multiple loci in human cells, including patient-derived induced pluripotent stem cells. Co-expression of RAD52 and dn53BP1 improves multiplexed HDR-mediated editing, whereas expression of RAD52 alone enhances HDR with Cas9 nickase. Our data show that the frequency of non-homologous end-joining-mediated double-strand break repair in the presence of these two factors is not suppressed and suggest that dn53BP1 competitively antagonizes 53BP1 to augment HDR in combination with RAD52. Importantly, co-expression of RAD52 and dn53BP1 does not alter Cas9 off-target activity. These findings support the use of RAD52 and dn53BP1 co-expression to overcome bottlenecks that limit HDR in precision genome editing.

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.

The present study aimed to identify genes exhibiting concomitant obesity-dependent changes in DNA methylation and gene expression in adipose tissues in the mouse using diet-induced obese (DIO) C57BL/6J and genetically obese ob/ob mice as models. Mature adipocytes were isolated from epididymal and inguinal adipose tissues of ob/ob and DIO C57BL/6J mice. DNA methylation was analyzed by MeDIP-sequencing and gene expression by microarray analysis. The majority of differentially methylated regions (DMRs) were hypomethylated in obese mice. Global methylation of long interspersed elements indicated that hypomethylation did not reflect methyl donor deficiency. In both DIO and ob/ob mice, we observed more obesity-associated methylation changes in epididymal than in inguinal adipocytes. Assignment of DMRs to promoter, exon, intron and intergenic regions demonstrated that DIO-induced changes in DNA methylation in C57BL/6J mice occurred primarily in exons, whereas inguinal adipocytes of ob/ob mice exhibited a higher enrichment of DMRs in promoter regions than in other regions of the genome, suggesting an influence of leptin on DNA methylation in inguinal adipocytes. We observed altered methylation and expression of 9 genes in epididymal adipocytes, including the known obesity-associated genes, Ehd2 and Kctd15, and a novel candidate gene, Irf8, possibly involved in immune type 1/type2 balance. The use of 2 obesity models enabled us to dissociate changes associated with high fat feeding from those associated with obesity per se. This information will be of value in future studies on the mechanisms governing the development of obesity and changes in adipocyte function associated with obesity.

There is increasing evidence from genome-wide association studies for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children, yet the effects of protein-coding variants on ALL risk have not been systematically evaluated. Here we show a missense variant in CDKN2A associated with the development of ALL at genome-wide significance (rs3731249, P=9.4 × 10(-23), odds ratio=2.23). Functional studies indicate that this hypomorphic variant results in reduced tumour suppressor function of p16(INK4A), increases the susceptibility to leukaemic transformation of haematopoietic progenitor cells, and is preferentially retained in ALL tumour cells. Resequencing the CDKN2A-CDKN2B locus in 2,407 childhood ALL cases reveals 19 additional putative functional germline variants. These results provide direct functional evidence for the influence of inherited genetic variation on ALL risk, highlighting the important and complex roles of CDKN2A-CDKN2B tumour suppressors in leukaemogenesis.

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.

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., (Copyright © 2015 the American Physiological Society.)

Building a population-specific catalogue of single nucleotide variants (SNVs), indels and structural variants (SVs) with frequencies, termed a national pan-genome, is critical for further advancing clinical and public health genetics in large cohorts. Here we report a Danish pan-genome obtained from sequencing 10 trios to high depth (50 × ). We report 536k novel SNVs and 283k novel short indels from mapping approaches and develop a population-wide de novo assembly approach to identify 132k novel indels larger than 10 nucleotides with low false discovery rates. We identify a higher proportion of indels and SVs than previous efforts showing the merits of high coverage and de novo assembly approaches. In addition, we use trio information to identify de novo mutations and use a probabilistic method to provide direct estimates of 1.27e-8 and 1.5e-9 per nucleotide per generation for SNVs and indels, respectively.

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.

Asthma exacerbations are among the most frequent causes of hospitalization during childhood, but the underlying mechanisms are poorly understood. We performed a genome-wide association study of a specific asthma phenotype characterized by recurrent, severe exacerbations occurring between 2 and 6 years of age in a total of 1,173 cases and 2,522 controls. Cases were identified from national health registries of hospitalization, and DNA was obtained from the Danish Neonatal Screening Biobank. We identified five loci with genome-wide significant association. Four of these, GSDMB, IL33, RAD50 and IL1RL1, were previously reported as asthma susceptibility loci, but the effect sizes for these loci in our cohort were considerably larger than in the previous genome-wide association studies of asthma. We also obtained strong evidence for a new susceptibility gene, CDHR3 (encoding cadherin-related family member 3), which is highly expressed in airway epithelium. These results demonstrate the strength of applying specific phenotyping in the search for asthma susceptibility genes.

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.

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.