ABSTRACT
Long-read RNA sequencing (lrRNA-seq) produces detailed information about full-length transcripts, including novel and sample-specific isoforms. Furthermore, there is an opportunity to call variants directly from lrRNA-seq data. However, most state-of-the-art variant callers have been developed for genomic DNA. Here, there are two objectives: first, we perform a mini-benchmark on GATK, DeepVariant, Clair3, and NanoCaller primarily on PacBio Iso-Seq, data, but also on Nanopore and Illumina RNA-seq data; second, we propose a pipeline to process spliced-alignment files, making them suitable for variant calling with DNA-based callers. With such manipulations, high calling performance can be achieved using DeepVariant on Iso-seq data.
Subject(s)
High-Throughput Nucleotide Sequencing , RNA , Sequence Analysis, RNA , RNA-Seq , Exome SequencingABSTRACT
Biological psychiatry scholarship on posttraumatic stress disorder (PTSD) is making strides with new omics technologies. In this context, there is growing recognition that noncoding RNAs are vital for the regulation of gene and protein expression. Long noncoding RNAs (lncRNAs) can modulate splicing, influence RNA editing, messenger RNA (mRNA) stability, translation activation, and microRNA-mRNA interactions, are highly abundant in the brain, and have been implicated in neurodevelopmental disorders. The largest subclass of lncRNAs is long intergenic noncoding RNAs (lincRNAs). We report on lincRNAs and their predicted mRNA targets associated with fear extinction induced by co-administration of D-cycloserine and behavioral fear extinction in a PTSD animal model. Forty-three differentially expressed lincRNAs and 190 differentially expressed mRNAs were found to be associated with fear extinction. Eight lincRNAs were predicted to interact with and regulate 108 of these mRNAs, while seven lincRNAs were predicted to interact with 22 of their pre-mRNA transcripts. Based on the functions of their target mRNAs, we inferred that these lincRNAs bind to nucleotides, ribonucleotides, and proteins; subsequently influence nervous system development, morphology, and immune system functioning; and could be associated with nervous system and mental health disorders. We found the quantitative trait loci that overlapped with fear extinction-related lincRNAs included traits such as serum corticosterone level, neuroinflammation, anxiety, stress, and despair-related responses. To the best of our knowledge, this is the first study to identify lincRNAs and their RNA targets with a putative role in transcriptional regulation during fear extinction in the context of an animal model of PTSD.
Subject(s)
Emotions , Gene Expression Profiling , Mental Disorders/etiology , Mental Disorders/psychology , Transcriptome , Animals , Computational Biology/methods , Cycloserine/adverse effects , Disease Models, Animal , Disease Susceptibility , Extinction, Psychological , Fear , Gene Expression Regulation , Gene Ontology , Male , Quantitative Trait Loci , RNA, Long Noncoding/genetics , RNA, Messenger/genetics , Rats , Stress Disorders, Post-Traumatic/etiology , Stress Disorders, Post-Traumatic/psychologyABSTRACT
In response to cortical stroke and unilateral corticospinal tract degeneration, compensatory sprouting of spared corticospinal fibers is associated with recovery of skilled movement in rodents. To date, little is known about the molecular mechanisms orchestrating this spontaneous rewiring. In this study, we provide insights into the molecular changes in the spinal cord tissue after large ischemic cortical injury in adult female mice, with a focus on factors that might influence the reinnervation process by contralesional corticospinal neurons. We mapped the area of cervical gray matter reinnervation by sprouting contralesional corticospinal axons after unilateral photothrombotic stroke of the motor cortex in mice using anterograde tracing. The mRNA profile of this reinnervation area was analyzed using whole-genome sequencing to identify differentially expressed genes at selected time points during the recovery process. Bioinformatic analysis revealed two phases of processes: early after stroke (4-7 d post-injury), the spinal transcriptome is characterized by inflammatory processes, including phagocytic processes as well as complement cascade activation. Microglia are specifically activated in the denervated corticospinal projection fields in this early phase. In a later phase (28-42 d post-injury), biological processes include tissue repair pathways with upregulated genes related to neurite outgrowth. Thus, the stroke-denervated spinal gray matter, in particular its intermediate laminae, represents a growth-promoting environment for sprouting corticospinal fibers originating from the contralesional motor cortex. This dataset provides a solid starting point for future studies addressing key elements of the post-stroke recovery process, with the goal to improve neuroregenerative treatment options for stroke patients.SIGNIFICANCE STATEMENT We show that the molecular changes in the spinal cord target tissue of the stroke-affected corticospinal tract are mainly defined by two phases: an early inflammatory phase during which microglia are specifically activated in the target area of reinnervating corticospinal motor neurons; and a late phase during which growth-promoting factors are upregulated which can influence the sprouting response, arborization, and synapse formation. By defining for the first time the endogenous molecular machinery in the stroke-denervated cervical spinal gray matter with a focus on promotors of axon growth through the growth-inhibitory adult CNS, this study will serve as a basis to address novel neuroregenerative treatment options for chronic stroke patients.
