Functional characterisation of the amyotrophic lateral sclerosis risk locus GPX3/TNIP1

Background: Amyotrophic lateral sclerosis (ALS) is a complex, late-onset, neurodegenerative disease with a genetic contribution to disease liability. Genome-wide association studies (GWAS) have identified ten risk loci to date, including the TNIP1/GPX3 locus on chromosome five. Given association analysis data alone cannot determine the most plausible risk gene for this locus, we undertook a comprehensive suite of in silico, in vivo and in vitro studies to address this. // Methods: The Functional Mapping and Annotation (FUMA) pipeline and five tools (conditional and joint analysis (GCTA-COJO), Stratified Linkage Disequilibrium Score Regression (S-LDSC), Polygenic Priority Scoring (PoPS), Summary-based Mendelian Randomisation (SMR-HEIDI) and transcriptome-wide association study (TWAS) analyses) were used to perform bioinformatic integration of GWAS data (Ncases = 20,806, Ncontrols = 59,804) with ‘omics reference datasets including the blood (eQTLgen consortium N = 31,684) and brain (N = 2581). This was followed up by specific expression studies in ALS case-control cohorts (microarray Ntotal = 942, protein Ntotal = 300) and gene knockdown (KD) studies of human neuronal iPSC cells and zebrafish-morpholinos (MO). // Results: SMR analyses implicated both TNIP1 and GPX3 (p < 1.15 × 10−6), but there was no simple SNP/expression relationship. Integrating multiple datasets using PoPS supported GPX3 but not TNIP1. In vivo expression analyses from blood in ALS cases identified that lower GPX3 expression correlated with a more progressed disease (ALS functional rating score, p = 5.5 × 10−3, adjusted R2 = 0.042, Beffect = 27.4 ± 13.3 ng/ml/ALSFRS unit) with microarray and protein data suggesting lower expression with risk allele (recessive model p = 0.06, p = 0.02 respectively). Validation in vivo indicated gpx3 KD caused significant motor deficits in zebrafish-MO (mean difference vs. control ± 95% CI, vs. control, swim distance = 112 ± 28 mm, time = 1.29 ± 0.59 s, speed = 32.0 ± 2.53 mm/s, respectively, p for all < 0.0001), which were rescued with gpx3 expression, with no phenotype identified with tnip1 KD or gpx3 overexpression. // Conclusions: These results support GPX3 as a lead ALS risk gene in this locus, with more data needed to confirm/reject a role for TNIP1. This has implications for understanding disease mechanisms (GPX3 acts in the same pathway as SOD1, a well-established ALS-associated gene) and identifying new therapeutic approaches. Few previous examples of in-depth investigations of risk loci in ALS exist and a similar approach could be applied to investigate future expected GWAS findings

Functional characterisation of the amyotrophic lateral sclerosis risk locus GPX3/TNIP1

Background Amyotrophic lateral sclerosis (ALS) is a complex, late-onset, neurodegenerative disease with a genetic contribution to disease liability. Genome-wide association studies (GWAS) have identified ten risk loci to date, including the TNIP1/GPX3 locus on chromosome five. Given association analysis data alone cannot determine the most plausible risk gene for this locus, we undertook a comprehensive suite of in silico, in vivo and in vitro studies to address this. Methods The Functional Mapping and Annotation (FUMA) pipeline and five tools (conditional and joint analysis (GCTA-COJO), Stratified Linkage Disequilibrium Score Regression (S-LDSC), Polygenic Priority Scoring (PoPS), Summary-based Mendelian Randomisation (SMR-HEIDI) and transcriptome-wide association study (TWAS) analyses) were used to perform bioinformatic integration of GWAS data (Ncases = 20,806, Ncontrols = 59,804) with ‘omics reference datasets including the blood (eQTLgen consortium N = 31,684) and brain (N = 2581). This was followed up by specific expression studies in ALS case-control cohorts (microarray Ntotal = 942, protein Ntotal = 300) and gene knockdown (KD) studies of human neuronal iPSC cells and zebrafish-morpholinos (MO). Results SMR analyses implicated both TNIP1 and GPX3 (p < 1.15 × 10−6), but there was no simple SNP/expression relationship. Integrating multiple datasets using PoPS supported GPX3 but not TNIP1. In vivo expression analyses from blood in ALS cases identified that lower GPX3 expression correlated with a more progressed disease (ALS functional rating score, p = 5.5 × 10−3, adjusted R2 = 0.042, Beffect = 27.4 ± 13.3 ng/ml/ALSFRS unit) with microarray and protein data suggesting lower expression with risk allele (recessive model p = 0.06, p = 0.02 respectively). Validation in vivo indicated gpx3 KD caused significant motor deficits in zebrafish-MO (mean difference vs. control ± 95% CI, vs. control, swim distance = 112 ± 28 mm, time = 1.29 ± 0.59 s, speed = 32.0 ± 2.53 mm/s, respectively, p for all < 0.0001), which were rescued with gpx3 expression, with no phenotype identified with tnip1 KD or gpx3 overexpression. Conclusions These results support GPX3 as a lead ALS risk gene in this locus, with more data needed to confirm/reject a role for TNIP1. This has implications for understanding disease mechanisms (GPX3 acts in the same pathway as SOD1, a well-established ALS-associated gene) and identifying new therapeutic approaches. Few previous examples of in-depth investigations of risk loci in ALS exist and a similar approach could be applied to investigate future expected GWAS findings

Using WormBase: A Genome Biology Resource for Caenorhabditis elegans and Related Nematodes

WormBase (www.wormbase.org) provides the nematode research community with a centralized database for information pertaining to nematode genes and genomes. As more nematode genome sequences are becoming available and as richer data sets are published, WormBase strives to maintain updated information, displays, and services to facilitate efficient access to and understanding of the knowledge generated by the published nematode genetics literature. This chapter aims to provide an explanation of how to use basic features of WormBase, new features, and some commonly used tools and data queries. Explanations of the curated data and step-by-step instructions of how to access the data via the WormBase website and available data mining tools are provided

Miniaturized Embryo Array for Automated Trapping, Immobilization and Microperfusion of Zebrafish Embryos

Zebrafish (Danio rerio) has recently emerged as a powerful experimental model in drug discovery and environmental toxicology. Drug discovery screens performed on zebrafish embryos mirror with a high level of accuracy the tests usually performed on mammalian animal models, and fish embryo toxicity assay (FET) is one of the most promising alternative approaches to acute ecotoxicity testing with adult fish. Notwithstanding this, automated in-situ analysis of zebrafish embryos is still deeply in its infancy. This is mostly due to the inherent limitations of conventional techniques and the fact that metazoan organisms are not easily susceptible to laboratory automation. In this work, we describe the development of an innovative miniaturized chip-based device for the in-situ analysis of zebrafish embryos. We present evidence that automatic, hydrodynamic positioning, trapping and long-term immobilization of single embryos inside the microfluidic chips can be combined with time-lapse imaging to provide real-time developmental analysis. Our platform, fabricated using biocompatible polymer molding technology, enables rapid trapping of embryos in low shear stress zones, uniform drug microperfusion and high-resolution imaging without the need of manual embryo handling at various developmental stages. The device provides a highly controllable fluidic microenvironment and post-analysis eleuthero-embryo stage recovery. Throughout the incubation, the position of individual embryos is registered. Importantly, we also for first time show that microfluidic embryo array technology can be effectively used for the analysis of anti-angiogenic compounds using transgenic zebrafish line (fli1a:EGFP). The work provides a new rationale for rapid and automated manipulation and analysis of developing zebrafish embryos at a large scale

Automated High-Content Live Animal Drug Screening Using C. elegans Expressing the Aggregation Prone Serpin α1-antitrypsin Z

The development of preclinical models amenable to live animal bioactive compound screening is an attractive approach to discovering effective pharmacological therapies for disorders caused by misfolded and aggregation-prone proteins. In general, however, live animal drug screening is labor and resource intensive, and has been hampered by the lack of robust assay designs and high throughput work-flows. Based on their small size, tissue transparency and ease of cultivation, the use of C. elegans should obviate many of the technical impediments associated with live animal drug screening. Moreover, their genetic tractability and accomplished record for providing insights into the molecular and cellular basis of human disease, should make C. elegans an ideal model system for in vivo drug discovery campaigns. The goal of this study was to determine whether C. elegans could be adapted to high-throughput and high-content drug screening strategies analogous to those developed for cell-based systems. Using transgenic animals expressing fluorescently-tagged proteins, we first developed a high-quality, high-throughput work-flow utilizing an automated fluorescence microscopy platform with integrated image acquisition and data analysis modules to qualitatively assess different biological processes including, growth, tissue development, cell viability and autophagy. We next adapted this technology to conduct a small molecule screen and identified compounds that altered the intracellular accumulation of the human aggregation prone mutant that causes liver disease in α1-antitrypsin deficiency. This study provides powerful validation for advancement in preclinical drug discovery campaigns by screening live C. elegans modeling α1-antitrypsin deficiency and other complex disease phenotypes on high-content imaging platforms

Development and Validation of an Automated High-Throughput System for Zebrafish In Vivo Screenings

The zebrafish is a vertebrate model compatible with the paradigms of drug discovery. The small size and transparency of zebrafish embryos make them amenable for the automation necessary in high-throughput screenings. We have developed an automated high-throughput platform for in vivo chemical screenings on zebrafish embryos that includes automated methods for embryo dispensation, compound delivery, incubation, imaging and analysis of the results. At present, two different assays to detect cardiotoxic compounds and angiogenesis inhibitors can be automatically run in the platform, showing the versatility of the system. A validation of these two assays with known positive and negative compounds, as well as a screening for the detection of unknown anti-angiogenic compounds, have been successfully carried out in the system developed. We present a totally automated platform that allows for high-throughput screenings in a vertebrate organism

Bioenergetic status modulates motor neuron vulnerability and pathogenesis in a zebrafish model of spinal muscular atrophy

Degeneration and loss of lower motor neurons is the major pathological hallmark of spinal muscular atrophy (SMA), resulting from low levels of ubiquitously-expressed survival motor neuron (SMN) protein. One remarkable, yet unresolved, feature of SMA is that not all motor neurons are equally affected, with some populations displaying a robust resistance to the disease. Here, we demonstrate that selective vulnerability of distinct motor neuron pools arises from fundamental modifications to their basal molecular profiles. Comparative gene expression profiling of motor neurons innervating the extensor digitorum longus (disease-resistant), gastrocnemius (intermediate vulnerability), and tibialis anterior (vulnerable) muscles in mice revealed that disease susceptibility correlates strongly with a modified bioenergetic profile. Targeting of identified bioenergetic pathways by enhancing mitochondrial biogenesis rescued motor axon defects in SMA zebrafish. Moreover, targeting of a single bioenergetic protein, phosphoglycerate kinase 1 (Pgk1), was found to modulate motor neuron vulnerability in vivo. Knockdown of pgk1 alone was sufficient to partially mimic the SMA phenotype in wild-type zebrafish. Conversely, Pgk1 overexpression, or treatment with terazosin (an FDA-approved small molecule that binds and activates Pgk1), rescued motor axon phenotypes in SMA zebrafish. We conclude that global bioenergetics pathways can be therapeutically manipulated to ameliorate SMA motor neuron phenotypes in vivo

Motor neuron-expressed microRNAs 218 and their enhancers are nested within introns of Slit2/3 genes

miR218-1 and miR218-2 are embedded in introns of SLIT2 and SLIT3, respectively, an arrangement conserved throughout vertebrate genomes. Both miR218 genes are predicted to be transcribed in the same orientation as their host genes and were assumed to be spliced from Slit2/3 primary transcripts. In zebrafish miR218 is active in cranial nerve motor nuclei and spinal cord motor neurons, while slit2 and slit3 are expressed predominantly in the midline. This differential expression pattern suggested independent regulation of miR218 genes by distinct enhancers. We tested conserved noncoding elements for regulatory activity by reporter gene transgenesis in zebrafish. Two human enhancers, 76 kb and 130 kb distant from miR218-2, were identified that drove GFP expression in zebrafish in an almost complete miR218 expression pattern. In the zebrafish slit3 locus, two enhancers with identical activity were discovered. In human SLIT2 one enhancer 52 kb upstream of miR218-1 drove an expression pattern very similar to the enhancers of miR218-2. This establishes that miR218-1/-2 regulatory units are nested within SLIT2/3 and that they are duplicates of an ancestral single locus. Due to the strong activity of the enhancers, unique transgenic lines were created that facilitate morphological and gene functional genetic experiments in motor neurons

Developmental suppression of schizophrenia-associated miR-137 alters sensorimotor function in zebrafish

The neurodevelopmentally regulated microRNA miR-137 was strongly implicated as risk locus for schizophrenia in the most recent genome wide association study coordinated by the Psychiatric Genome Consortium (PGC). This molecule is highly conserved in vertebrates enabling the investigation of its function in the developing zebrafish. We utilized this model system to achieve overexpression and suppression of miR-137, both transiently and stably through transgenesis. While miR-137 overexpression was not associated with an observable specific phenotype, downregulation by antisense morpholino and/or transgenic expression of miR-sponge RNA induced significant impairment of both embryonic and larval touch-sensitivity without compromising overall anatomical development. We observed miR-137 expression and activity in sensory neurons including Rohon–Beard neurons and dorsal root ganglia, two neuronal cell types that confer touch-sensitivity in normal zebrafish, suggesting a role of these cell types in the observed phenotype. The lack of obvious anatomical or histological pathology in these cells, however, suggested that subtle axonal network defects or a change in synaptic function and neural connectivity might be responsible for the behavioral phenotype rather than a change in the cellular morphology or neuroanatomy