Differential transcript usage in the Parkinson's disease brain

Nigral transcriptomic profiles in Engrailed-1 hemizygous mouse models of Parkinson's disease reveal upregulation of oxidative phosphorylation-related genes associated with delayed dopaminergic neurodegeneration

Introduction

Parkinson's disease (PD) is the second most common neurodegenerative disorder, increasing both in terms of prevalence and incidence. To date, only symptomatic treatment is available, highlighting the need to increase knowledge on disease etiology in order to develop new therapeutic strategies. Hemizygosity for the gene Engrailed-1 (En1), encoding a conserved transcription factor essential for the programming, survival, and maintenance of midbrain dopaminergic neurons, leads to progressive nigrostriatal degeneration, motor impairment and depressive-like behavior in SwissOF1 (OF1-En1+/-). The neurodegenerative phenotype is, however, absent in C57Bl/6j (C57-En1+/-) mice. En1+/- mice are thus highly relevant tools to identify genetic factors underlying PD susceptibility.

Methods

Transcriptome profiles were defined by RNAseq in microdissected substantia nigra from 1-week old OF1, OF1- En1+/-, C57 and C57- En1+/- male mice. Differentially expressed genes (DEGs) were analyzed for functional enrichment. Neurodegeneration was assessed in 4- and 16-week old mice by histology.

Results

Nigrostriatal neurodegeneration was manifested in OF1- En1+/- mice by increased dopaminergic striatal axonal swellings from 4 to 16 weeks and decreased number of dopaminergic neurons in the SNpc at 16 weeks compared to OF1. In contrast, C57- En1+/- mice had no significant increase in axonal swellings or cell loss in SNpc at 16 weeks. Transcriptomic analyses identified 198 DEGs between OF1- En1+/- and OF1 mice but only 52 DEGs between C57- En1+/- and C57 mice. Enrichment analysis of DEGs revealed that the neuroprotective phenotype of C57- En1+/- mice was associated with a higher expression of oxidative phosphorylation-related genes compared to both C57 and OF1- En1+/- mice.

Discussion

Our results suggest that increased expression of genes encoding mitochondrial proteins before the onset of neurodegeneration is associated with increased resistance to PD-like nigrostriatal neurodegeneration. This highlights the importance of genetic background in PD models, how different strains can be used to model clinical and sub-clinical pathologies and provides insights to gene expression mechanisms associated with PD susceptibility and progression.

Comprehensive benchmark of differential transcript usage analysis for static and dynamic conditions

RNA sequencing offers unique insights into transcriptome diversity, and a plethora of tools have been developed to analyze alternative splicing. One important task is to detect changes in the relative transcript abundance in differential transcript usage (DTU) analysis. The choice of the right analysis tool is non-trivial and depends on experimental factors such as the availability of single- or paired-end and bulk or single-cell data. To help users select the most promising tool for their task, we performed a comprehensive benchmark of DTU detection tools. We cover a wide array of experimental settings, using simulated bulk and single-cell RNA-seq data as well as real transcriptomics datasets, including time-series data. Our results suggest that DEXSeq, edgeR, and LimmaDS are better choices for paired-end data, while DSGseq and DEXSeq can be used for single-end data. In single-cell simulation settings, we showed that satuRn performs better than DTUrtle. In addition, we showed that Spycone is optimal for time series DTU/IS analysis based on the evidence provided using GO terms enrichment analysis.

Mutations in DARS2 result in global dysregulation of mRNA metabolism and splicing

Leukoencephalopathy with brainstem and spinal cord involvement and lactate elevation (LBSL) is a rare neurological disorder caused by the mutations in the DARS2 gene, which encodes the mitochondrial aspartyl-tRNA synthetase. The objective of this study was to understand the impact of DARS2 mutations on cell processes through evaluation of LBSL patient stem cell derived cerebral organoids and neurons. We generated human cerebral organoids (hCOs) from induced pluripotent stem cells (iPSCs) of seven LBSL patients and three healthy controls using an unguided protocol. Single cells from 70-day-old hCOs were subjected to SMART-seq2 sequencing and bioinformatic analysis to acquire high-resolution gene and transcript expression datasets. Global gene expression analysis demonstrated dysregulation of a number of genes involved in mRNA metabolism and splicing processes within LBSL hCOs. Importantly, there were distinct and divergent gene expression profiles based on the nature of the DARS2 mutation. At the transcript level, pervasive differential transcript usage and differential spliced exon events that are involved in protein translation and metabolism were identified in LBSL hCOs. Single-cell analysis of DARS2 (exon 3) showed that some LBSL cells exclusively express transcripts lacking exon 3, indicating that not all LBSL cells can benefit from the "leaky" nature common to splice site mutations. At the gene- and transcript-level, we uncovered that dysregulated RNA splicing, protein translation and metabolism may underlie at least some of the pathophysiological mechanisms in LBSL. To confirm hCO findings, iPSC-derived neurons (iNs) were generated by overexpressing Neurogenin 2 using lentiviral vector to study neuronal growth, splicing of DARS2 exon 3 and DARS2 protein expression. Live cell imaging revealed neuronal growth defects of LBSL iNs, which was consistent with the finding of downregulated expression of genes related to neuronal differentiation in LBSL hCOs. DARS2 protein was downregulated in iNs compared to iPSCs, caused by increased exclusion of exon 3. The scope and complexity of our data imply that DARS2 is potentially involved in transcription regulation beyond its canonical role of aminoacylation. Nevertheless, our work highlights transcript-level dysregulation as a critical, and relatively unexplored, mechanism linking genetic data with neurodegenerative disorders.

Mutations in DARS2 result in global dysregulation of mRNA metabolism and splicing

Leukoencephalopathy with brain stem and spinal cord involvement and lactate elevation (LBSL) is a rare neurological disorder caused by the mutations in the DARS2 gene, which encodes the mitochondrial aspartyl-tRNA synthetase. The objective of this study was to understand the impact of DARS2 mutations on cell processes through evaluation of LBSL patient stem cell derived cerebral organoids and neurons. We generated human cerebral organoids (hCOs) from induced pluripotent stem cells (iPSCs) of seven LBSL patients and three healthy controls using an unguided protocol. Single cells from 70-day-old hCOs underwent SMART-seq2 sequencing and multiple bioinformatic analysis tools were applied to high-resolution gene and transcript expression analyses. To confirm hCO findings, iPSC-derived neurons (iNs) were generated by overexpressing Neurogenin 2 using lentiviral vector to study neuronal growth, splicing of DARS2 exon 3 and DARS2 protein expression. Global gene expression analysis demonstrated dysregulation of a number of genes involved in mRNA metabolism and splicing processes within LBSL hCOs. Importantly, there were distinct and divergent gene expression profiles based on the nature of the DARS2 mutation. At the transcript level, pervasive differential transcript usage and differential spliced exon events that are involved in protein translation and metabolism were identified in LBSL hCOs. Single-cell analysis of DARS2 (exon 3) showed that some LBSL cells exclusively express transcripts lacking exon 3, indicating that not all LBSL cells can benefit from the "leaky" nature common to splice site mutations. Live cell imaging revealed neuronal growth defects of LBSL iNs, which was consistent with the finding of downregulated expression of genes related to neuronal differentiation in LBSL hCOs. DARS2 protein was downregulated in iNs compared to iPSCs, caused by increased exclusion of exon 3. At the gene- and transcript-level, we uncovered that dysregulated RNA splicing, protein translation and metabolism may underlie at least some of the pathophysiological mechanisms in LBSL. The scope and complexity of our data imply that DARS2 is potentially involved in transcription regulation beyond its canonical role of aminoacylation. Nevertheless, our work highlights transcript-level dysregulation as a critical, and relatively unexplored, mechanism linking genetic data with neurodegenerative disorders.

Transcriptome analyses in infertile men reveal germ cell-specific expression and splicing patterns

The process of spermatogenesis-when germ cells differentiate into sperm-is tightly regulated, and misregulation in gene expression is likely to be involved in the physiopathology of male infertility. The testis is one of the most transcriptionally rich tissues; nevertheless, the specific gene expression changes occurring during spermatogenesis are not fully understood. To better understand gene expression during spermatogenesis, we generated germ cell-specific whole transcriptome profiles by systematically comparing testicular transcriptomes from tissues in which spermatogenesis is arrested at successive steps of germ cell differentiation. In these comparisons, we found thousands of differentially expressed genes between successive germ cell types of infertility patients. We demonstrate our analyses' potential to identify novel highly germ cell-specific markers (TSPY4 and LUZP4 for spermatogonia; HMGB4 for round spermatids) and identified putatively misregulated genes in male infertility (RWDD2A, CCDC183, CNNM1, SERF1B). Apart from these, we found thousands of genes showing germ cell-specific isoforms (including SOX15, SPATA4, SYCP3, MKI67). Our approach and dataset can help elucidate genetic and transcriptional causes for male infertility.

The landscape of isoform switches in sepsis: a multicenter cohort study

Sepsis is caused by an uncontrolled inflammatory response, whose underlying mechanisms are not fully understood. It is well known that the majority of human genes can be expressed as alternative isoforms. While isoform switching is implicated in many diseases and is particularly prominent in cancer, it has never been reported in the context of sepsis. Patients presented to the emergency department of three tertiary care hospitals from January 2020 to December 2020 were enrolled. Clinical variables and genome-wide transcriptome of peripheral blood mononuclear cells (PBMC) were obtained. Isoform switching analysis were performed to identify significant isoform switches and relevant biological consequences. A total of 48 subjects with sepsis, involving 42 survivors and 6 non-survivors, admitted to the emergency department of three tertiary care hospitals were enrolled in this study. PBMCs were extracted for RNA sequencing (RNA-seq). Patients (n = 4) with mild stroke or acute coronary syndrome without infection were enrolled in this study as controls. The most frequent functional changes resulting from isoform switching were changes affecting the open reading frame, protein domains and intron retention. Many genes without differences in gene expression showed significant isoform switching. Many genes with significant isoform switches ([Formula: see text]> 0.1) were associated with higher mortality risk, including PIGS, CASP3, LITAF, HBB and RUVBL2. The study for the first time described the landscape of isoform switching in sepsis, including differentially expressed isoform fractions between patients with and without sepsis and survivors and nonsurvivors. The biological consequences of isoform switching, including protein domain loss, signal peptide gain, and intron retention, were identified.

Long-Read RNA Sequencing Identifies Polyadenylation Elongation and Differential Transcript Usage of Host Transcripts During SARS-CoV-2 In Vitro Infection

Better methods to interrogate host-pathogen interactions during Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections are imperative to help understand and prevent this disease. Here we implemented RNA-sequencing (RNA-seq) using Oxford Nanopore Technologies (ONT) long-reads to measure differential host gene expression, transcript polyadenylation and isoform usage within various epithelial cell lines permissive and non-permissive for SARS-CoV-2 infection. SARS-CoV-2-infected and mock-infected Vero (African green monkey kidney epithelial cells), Calu-3 (human lung adenocarcinoma epithelial cells), Caco-2 (human colorectal adenocarcinoma epithelial cells) and A549 (human lung carcinoma epithelial cells) were analyzed over time (0, 2, 24, 48 hours). Differential polyadenylation was found to occur in both infected Calu-3 and Vero cells during a late time point (48 hpi), with Gene Ontology (GO) terms such as viral transcription and translation shown to be significantly enriched in Calu-3 data. Poly(A) tails showed increased lengths in the majority of the differentially polyadenylated transcripts in Calu-3 and Vero cell lines (up to ~101 nt in mean poly(A) length, padj = 0.029). Of these genes, ribosomal protein genes such as RPS4X and RPS6 also showed downregulation in expression levels, suggesting the importance of ribosomal protein genes during infection. Furthermore, differential transcript usage was identified in Caco-2, Calu-3 and Vero cells, including transcripts of genes such as GSDMB and KPNA2, which have previously been implicated in SARS-CoV-2 infections. Overall, these results highlight the potential role of differential polyadenylation and transcript usage in host immune response or viral manipulation of host mechanisms during infection, and therefore, showcase the value of long-read sequencing in identifying less-explored host responses to disease.

Systematic analysis of inflammation and pain pathways in a mouse model of gout

Gout is a prevalent and painful inflammatory arthritis, and its global burden continues to rise. Intense pain induced by gout attacks is a major complication of gout. However, systematic studies of gout inflammation and pain are lacking. Using a monosodium urate (MSU) crystal-induced gout model, we performed genome-wide transcriptome analysis of the inflamed ankle joint, dorsal root ganglion (DRG), and spinal cord of gouty mice. Our results revealed important transcriptional changes, including highly elevated inflammation and broad activation of immune pathways in both the joint and the nervous system, in gouty mice. Integrated analysis showed that there was a remarkable overlap between our RNAseq and human genome-wide association study (GWAS) of gout; for example, the risk gene, stanniocalcin-1 (STC1) showed significant upregulation in all three tissues. Interestingly, when compared to the transcriptomes of human osteoarthritis (OA) and rheumatoid arthritis (RA) joint tissues, we identified significant upregulation of cAMP/cyclic nucleotide-mediated signaling shared between gouty mice and human OA with high knee pain, which may provide excellent drug targets to relieve gout pain. Furthermore, we investigated the common and distinct transcriptomic features of gouty, inflammatory pain, and neuropathic pain mouse models in their DRG and spinal cord tissues. Moreover, we discovered distinct sets of genes with significant differential alternative splicing or differential transcript usage in each tissue, which were largely not detected by conventional differential gene expression analysis approaches. Based on these results, our study provided a more accurate and comprehensive depiction of transcriptomic alterations related to gout inflammation and pain.

Characterisation of the Function of a SINE-VNTR-Alu Retrotransposon to Modulate Isoform Expression at the MAPT Locus

SINE-VNTR-Alu retrotransposons represent one class of transposable elements which contribute to the regulation and evolution of the primate genome and have the potential to be involved in genetic instability and disease progression. However, these polymorphic elements have not been extensively analysed when addressing the missing heritability of neurodegenerative diseases, including Parkinson’s disease (PD) and amyotrophic lateral sclerosis (ALS). SVA_67, a retrotransposon insertion polymorphism, is located in a 1.8 Mb region of high linkage disequilibrium, called the MAPT locus, which is known to contribute to increased risk of developing PD, frontotemporal dementia and other tauopathies. To investigate the role of SVA_67 in directing differential gene expression at this locus, we characterised the impact of SVA_67 allele dosage on isoform expression of several genes in the MAPT locus using the datasets from both the Parkinson’s Progression Markers Initiative and New York Genome Center Consortium Target ALS cohort. The Parkinson’s data was from gene expression in the blood and the ALS data from a variety of CNS regions and allowed us to demonstrate that SVA_67 presence or absence correlated with both isoform- and tissue-specific expression of multiple genes at this locus. This study highlights the importance of addressing SVA polymorphism in disease genetics to gain insight into a better understanding of the role of these regulatory domains to a variety of neurodegenerative diseases.

Transcriptomic analysis of frontotemporal lobar degeneration with TDP-43 pathology reveals cellular alterations across multiple brain regions

Frontotemporal lobar degeneration (FTLD) is a group of heterogeneous neurodegenerative disorders affecting the frontal and temporal lobes of the brain. Nuclear loss and cytoplasmic aggregation of the RNA-binding protein TDP-43 represents the major FTLD pathology, known as FTLD-TDP. To date, there is no effective treatment for FTLD-TDP due to an incomplete understanding of the molecular mechanisms underlying disease development. Here we compared postmortem tissue RNA-seq transcriptomes from the frontal cortex, temporal cortex, and cerebellum between 28 controls and 30 FTLD-TDP patients to profile changes in cell-type composition, gene expression and transcript usage. We observed downregulation of neuronal markers in all three regions of the brain, accompanied by upregulation of microglia, astrocytes, and oligodendrocytes, as well as endothelial cells and pericytes, suggesting shifts in both immune activation and within the vasculature. We validate our estimates of neuronal loss using neuropathological atrophy scores and show that neuronal loss in the cortex can be mainly attributed to excitatory neurons, and that increases in microglial and endothelial cell expression are highly correlated with neuronal loss. All our analyses identified a strong involvement of the cerebellum in the neurodegenerative process of FTLD-TDP. Altogether, our data provides a detailed landscape of gene expression alterations to help unravel relevant disease mechanisms in FTLD.

Transcriptomic analysis of frontotemporal lobar degeneration with TDP-43 pathology reveals cellular alterations across multiple brain regions

Frontotemporal lobar degeneration (FTLD) is a group of heterogeneous neurodegenerative disorders affecting the frontal and temporal lobes of the brain. Nuclear loss and cytoplasmic aggregation of the RNA-binding protein TDP-43 represents the major FTLD pathology, known as FTLD-TDP. To date, there is no effective treatment for FTLD-TDP due to an incomplete understanding of the molecular mechanisms underlying disease development. Here we compared postmortem tissue RNA-seq transcriptomes from the frontal cortex, temporal cortex, and cerebellum between 28 controls and 30 FTLD-TDP patients to profile changes in cell-type composition, gene expression and transcript usage. We observed downregulation of neuronal markers in all three regions of the brain, accompanied by upregulation of microglia, astrocytes, and oligodendrocytes, as well as endothelial cells and pericytes, suggesting shifts in both immune activation and within the vasculature. We validate our estimates of neuronal loss using neuropathological atrophy scores and show that neuronal loss in the cortex can be mainly attributed to excitatory neurons, and that increases in microglial and endothelial cell expression are highly correlated with neuronal loss. All our analyses identified a strong involvement of the cerebellum in the neurodegenerative process of FTLD-TDP. Altogether, our data provides a detailed landscape of gene expression alterations to help unravel relevant disease mechanisms in FTLD.

Machine Learning Identifies Six Genetic Variants and Alterations in the Heart Atrial Appendage as Key Contributors to PD Risk Predictivity

Parkinson's disease (PD) is a complex neurodegenerative disease with a range of causes and clinical presentations. Over 76 genetic loci (comprising 90 SNPs) have been associated with PD by the most recent GWAS meta-analysis. Most of these PD-associated variants are located in non-coding regions of the genome and it is difficult to understand what they are doing and how they contribute to the aetiology of PD. We hypothesised that PD-associated genetic variants modulate disease risk through tissue-specific expression quantitative trait loci (eQTL) effects. We developed and validated a machine learning approach that integrated tissue-specific eQTL data on known PD-associated genetic variants with PD case and control genotypes from the Wellcome Trust Case Control Consortium. In so doing, our analysis ranked the tissue-specific transcription effects for PD-associated genetic variants and estimated their relative contributions to PD risk. We identified roles for SNPs that are connected with INPP5P, CNTN1, GBA and SNCA in PD. Ranking the variants and tissue-specific eQTL effects contributing most to the machine learning model suggested a key role in the risk of developing PD for two variants (rs7617877 and rs6808178) and eQTL associated transcriptional changes of EAF1-AS1 within the heart atrial appendage. Similarly, effects associated with eQTLs located within the Brain Cerebellum were also recognized to confer major PD risk. These findings were replicated in two additional, independent cohorts (the UK Biobank, and NeuroX) and thus warrant further mechanistic investigations to determine if these transcriptional changes could act as early contributors to PD risk and disease development.

RNA Modifications and RNA Metabolism in Neurological Disease Pathogenesis

The intrinsic cellular heterogeneity and molecular complexity of the mammalian nervous system relies substantially on the dynamic nature and spatiotemporal patterning of gene expression. These features of gene expression are achieved in part through mechanisms involving various epigenetic processes such as DNA methylation, post-translational histone modifications, and non-coding RNA activity, amongst others. In concert, another regulatory layer by which RNA bases and sugar residues are chemically modified enhances neuronal transcriptome complexity. Similar RNA modifications in other systems collectively constitute the cellular epitranscriptome that integrates and impacts various physiological processes. The epitranscriptome is dynamic and is reshaped constantly to regulate vital processes such as development, differentiation and stress responses. Perturbations of the epitranscriptome can lead to various pathogenic conditions, including cancer, cardiovascular abnormalities and neurological diseases. Recent advances in next-generation sequencing technologies have enabled us to identify and locate modified bases/sugars on different RNA species. These RNA modifications modulate the stability, transport and, most importantly, translation of RNA. In this review, we discuss the formation and functions of some frequently observed RNA modifications—including methylations of adenine and cytosine bases, and isomerization of uridine to pseudouridine—at various layers of RNA metabolism, together with their contributions to abnormal physiological conditions that can lead to various neurodevelopmental and neurological disorders.