Transcriptome of dorsal root ganglia caudal to a spinal cord injury with modulated behavioral activity

Spinal cord injury (SCI) is a devastating clinical condition resulting in significant disabilities. Apart from local injury within the spinal cord, SCI patients develop a myriad of complications including multi-organ dysfunction. Some of the dysfunctions may be directly or indirectly related to the sensory neurons of the dorsal root ganglia (DRG), which signal to both the spinal cord and the peripheral organs. After SCI, some classes of DRG neurons exhibit sensitization and undergo axonal sprouting both peripherally and centrally. Such physiological and anatomical re-organization after SCI contributes to both adaptive and maladaptive plasticity processes, which may be modulated by activity and exercise. In this study, we collected comprehensive gene expression data in whole DRG below the levels of the injury to compare the effects of SCI with and without two different forms of exercise in rats.


Background & Summary
Injury to the spinal cord leads to many complications, ranging from the more observable such as loss of motor control, to the less observable, such as multi-organ dysfunction, chronic pain, and autonomic dysregulation. With reduced or absent descending control (pathways originating in the brain that project onto the spinal cord 1 , leading to what is known as supraspinal input) following SCI, spinal circuitry becomes highly reliant on information coming from the periphery and carried by the sensory neurons 2 . This makes understanding the response of sensory neurons to SCI and treatments of paramount importance. A growing body of work indicates that although the vast majority of spinal sensory neurons housed in DRG do not appear to be directly injured by SCI, alterations in sensory neuron anatomy and physiology not only contribute to life-threatening conditions that develop after SCI but also contribute to some forms of recovered functions. Below the level of injury, some sensory neurons undergo intra-spinal axonal sprouting which provides new input to spinal circuitry 3,4 . In one example, sprouting of nociceptive axons with concomitant SCI-induced loss of descending controls allows sensory neurons unchecked access to sensory, autonomic, and motor circuitry, which contributes to chronic pain, autonomic dysreflexia, and impaired locomotion [5][6][7] . However, such plasticity can also be beneficial in enhancing outcomes of locomotor training, where inputs from skin and muscle are major drivers [8][9][10] . In addition to inducing intra-spinal axonal sprouting, SCI also leads to persistent hyperexcitability and spontaneous activity in some neuronal cell bodies and peripheral branches of DRG neurons 11 , which could be a major contributor to developing chronic pain and autonomic dysregulation.
Previous studies have shown spontaneous improvement in locomotor function following sufficient "in-cage activity" after SCI [12][13][14][15] . Therefore, as a control, this study analyzed DRG from rats housed in locomotor-restricting (tiny) cages where in-cage activity was reduced by 75-80% for both SCI and naïve groups (DSKM, unpublished observations). This type of housing mimics the clinical situation in which patients are bed-and wheelchairbound for weeks or months following injury. Spinal cord injured rats, all of which were housed in tiny cages, were divided into two injury groups: contusion SCI (CONT SCI), which resulted in some sparing of supraspinal input, and complete spinal cord transection (TX SCI), which resulted in the total loss of descending, supraspinal input. In addition, subsets of CONT SCI animals were randomly assigned to two exercise training paradigms: swimming (SWIM) and shallow water walking (SWW). These models of enhanced activity/exercise training have beneficial effects in terms of recovery of function and complement our previous study of transcriptional responses in liver [14][15][16][17][18] . Further, several clinical and pre-clinical studies suggest that physical activity and exercise can be effective in enhancing recovery of function and attenuating SCI-induced neuropathic pain and autonomic dysregulation in both humans and rodents 7,[19][20][21][22] . Exercise training ameliorates chronic pain in rodents and is associated with reducing SCI-induced intra-spinal axonal sprouting of subsets of afferents 7,20 . Exercise can also influence sensory neuron function and transcription in intact, nerve-injured, and SCI conditions 7,23,24 . Intrinsic changes in gene expression could underlie SCI-induced neuronal plasticity (sprouting and sensitization) in DRG neurons and the prevention/reversal with exercise.
These data were collected to evaluate the transcriptomic changes in DRG below the level of SCI and to better understand how this profile is altered with exercise training. Because the study was designed to identify broad effects at the transcriptional level, as opposed to organ/tissue/segment-specific effects, bilateral DRG from many segments below the injury were pooled within animals. Our technical validation indicates that the data are of high quality with gene expression occurring in expected gene regions. UCSC Genome Browser expression tracks are available with the raw data to facilitate exploration of the samples. The data described in this paper are offered to the scientific community as a valuable resource for future investigations into the molecular changes in DRG occurring in SCI-induced pathologies. . We chose to work exclusively with female adult rats to control for transcriptomic responses related to age and gender. Also, in our experience, female rats recover more quickly from surgery and are more motivated to exercise, providing better outcome measures for our experimental purposes. All rats were housed initially in standard cages and maintained in a continuous 12 h-light/dark cycle. Tap water and a standard rodent diet were available to all rats ad libitum. experimental design and SCI. The experimental workflow is displayed in Fig. 1 with the experimental timeline for SCI with and without exercise, given in greater detail in Fig. 2 www.nature.com/scientificdata www.nature.com/scientificdata/ were randomly assigned to five groups: naïve (No SCI, 4 replicates), T2 transection injury (TX SCI, 4 replicates), T2 contusion injury (CONT SCI, 6 replicates), T2 contusion injury followed by swimming as exercise (CONT SCI + SWIM, 5 replicates), T2 contusion injury followed by shallow water walking as exercise (CONT SCI + SWW, 5 replicates). Throughout the study, rats were doubly-housed with animals in the same experimental group.
All rats were initially gentled and acclimatized to testing and exercise facilities for two weeks. After this period, animals were anaesthetized with ketamine (50 mg/kg)/xylazine (0.024 mg/kg)/acepromazine (0.005 mg/kg) cocktail (IP) and glycopyrolate (0.08 mg/kg, IM) was given prior to SCI surgeries. For all injury groups (both TX SCI and CONT SCI), a dorsal midline incision was made in the superficial muscle overlying the T1-T3 vertebrae. A single level laminectomy was performed on the T2 vertebra; dura mater was kept intact. For animals in the TX SCI group, complete transection of the spinal cord at T2 was performed with a scalpel. Animals in the CONT SCI group received a moderately-severe contusion injury (25 g-cm) at the T2 spinal cord level using the NYU Impactor 25,26 . Injured animals were monitored on heating pads until they recovered from the anaesthesia. Rats were then double housed in standard cages with ALPHA-dri ® bedding (Shepherd's TM Specialty Paper, Milford, NJ). Post operative care consisted of prophylactic antibiotics (gentamicin sulphate, 20 mg/kg, SC), given for 7 consecutive days. For pain management injured animals received buprenorphine (0.03 mg/kg, SC) twice daily for three days and as needed thereafter. Animals were also given 5 ml boluses of lactated Ringer's solution for three days (and as needed for hydration thereafter). Manual bladder expression was performed thrice daily for each animal until reflexive voiding was re-established.
Three days after injury, all animals (naïve included) were doubly-housed with animals in the same experimental group in tiny cages (7.5′′ × 8.5′′ × 8′′) rather than standard cages (22′′ × 12.5′′ × 8′′) to restrict their overground locomotion for the duration of the study. If mortality resulted in an odd number of animals, then an injured animal not in the study was used as a cage mate. For animals in the exercise training groups, swimming or shallow water walking began 2.5 weeks post-injury and was conducted 5 consecutive days a week for 10 weeks. Animals exercised for a total of 30 minutes each day with 15 minutes of exercise in the morning and 15 minutes in the afternoon. The morning and afternoon sessions were separated by a minimum of one hour. Each 15 minute session consisted of three 5 minute periods of exercise with breaks between the periods lasting approximately 20-25 minutes. tissue collection and RNA extraction. Animals were sacrificed with a ketamine overdose at 8.5, 11.5 or 13.5 weeks post-surgery depending on experimental condition (see Fig. 2). Hearts were arrested in diastole with an injection of 3 M KCl. Animals were perfused with PBS supplemented with 20% RNAlater (Thermo Fisher Scientific, Waltham, MA). For each animal, T11-L4 DRG were retrieved from both sides of the spinal column and collected together in QIAzol (Qiagen, Germantown, MD). Total RNA from each sample was isolated with RNeasy Lipid Tissue Mini Kit (Qiagen).
Library preparation and sequencing. 1 ug of total RNA samples were used for poly A enrichment. First and second strands were synthesised followed by 3′ end adenylation. Samples were barcoded with Illumina TrueSeq adapters. 1.8 pM of barcoded library was denatured, and sequencing was performed on the University of Louisville Center for Genetics and Molecular Medicine (CGeMM) Illumina NextSeq 500 using the NextSeq 500/550 1X75 cycle high output kit.   Fig. 2 Experimental Paradigm. Rats were gentled for two weeks. After two weeks, a moderate-severe T2 contusion (Groups 3, 4 and 5) or complete T2 transection (Group 2) was performed. After 3~4 days, rats in all groups were transferred to tiny cages. Two weeks post-transfer, a subset of rats with contusion SCI were given swimming or shallow water walking exercise training for 10 weeks (5 days per week). All groups were euthanized after designated time points for terminal assessments.
www.nature.com/scientificdata www.nature.com/scientificdata/ RNA-seq data analysis. Sequencing produced nearly a billion single end reads across the 24 samples.
The vast majority of read lengths fell between 74 and 76 bases across all samples. The quality of the reads was assessed using FastQC v.0.10.1 27 , which indicated no sequence trimming was necessary. The sequences were directly aligned to the Rattus norvegicus reference genome assembly (Rn6) using the STAR aligner v.2.6 28 . Read counts were extracted with HTSeq version 0.10.0 29 , using Rn6 Ensembl annotations (build 93) to locate gene regions 30 . Prior to extracting read counts, mitochondrial genes were removed from the annotation file to reduce non-relevant variation in future analyses. This produced an annotation file testing 24,613 regions. Prior to examining expression, raw read counts were normalized using DESeq2's default method, relative log expression (RLE) 31,32 . UCSC Genome Browser tracks were created to easily explore expression across the tested sites 33 . The tracks were created using guidelines and utilities available on the UCSC Genome Browser website. Much of this involved converting sequencing alignment files in BAM format to BigWig format, required for use on the browser.
A principal component analysis (PCA) was performed using the plotPCA function in DESeq2. Prior to the analysis, a variance stabilizing transformation was performed on the raw read count matrix to achieve similar variance at all levels of mean expression across the samples. This minimizes the effect of highly varying genes on the spread of points in the PCA plot. The PCA plot was generated using the R programming language package ggplot2 34 .

Data Records
The data were submitted to NCBI Gene Expression Omnibus (GSE125630) 35 . This GEO project includes raw data in Fastq format, raw HTSeq 29 counts, and UCSC Genome Browser 33 tracks in bigwig format for all samples. This dataset is part of a larger study measuring the systemic transcriptional response to spinal cord injury, including liver 16 and soleus muscle, all of which are included as part of a GEO superseries (GSE129704) 36 .

technical Validation
RNA metrics: Sequencing generated 22.8 to 54.7 million reads/sample with a mean of 40.7 million and standard deviation of 8 million. Table 1 displays the number of raw reads successfully aligned for each of the samples. The alignment rate for uniquely mapped and multi-mapped reads combined ranged from 98 to 98.6 percent with a mean of 98.4 across the 24 samples. Quality assessment. Raw sequencing data was input to FastQC for quality assessment. All samples were deemed of high quality. In Fig. 3a, the Phred quality score per base is displayed for a representative sample from each experimental group. With the exception of the last base in some samples, the 25 th percentile of quality scores is at or above a Phred score of 30, reflecting 99.9 percent accuracy in base calling. When scores at the 25 th percentile drop below this point on the last base, they remain in the reasonable quality category. A gradual drop in quality toward the end of the sequenced read is inherent to Illumina's approach to sequencing by synthesis 37 . www.nature.com/scientificdata www.nature.com/scientificdata/ Gene expression variation of biological replicates. PCA was performed to assess the within-and between-group variation of the samples. The samples were separated into two distinct comparison sets, one appropriate for an examination of the effect of exercise on SCI and one appropriate for an examination of injury severity. Figure 3b displays the PCA plot for samples comparing CONT SCI alone with CONT SCI + SWIM/ SWW. Samples within each experimental group cluster together, and all experimental groups lie in different regions of the graph. Sixty-six percent of the variance is explained by the separation of CONT SCI + SWIM/ SWW samples and No SCI/CONT SCI samples. Figure 3c displays the PCA plot for samples with different types of injury. Once again, samples within each experimental group cluster together, and all experimental groups are located in different regions of the graph. Sixty-four percent of the variance is explained by the distinction between TX SCI samples, the more severe SCI, and No SCI/CONT SCI samples. www.nature.com/scientificdata www.nature.com/scientificdata/ Gene expression was confirmed in SCI samples for genes relevant to sensory functions of particular interest. Genes associated with ion channels, solute carriers, and axonogenesis were selected from Ensembl (version 93) and Rat Genome Database gene function descriptions for rat 30 . Genes identified as pain markers were selected from the Pain Research Forum 38 . In Fig. 4, mean expression for the CONT SCI samples is displayed for the highest expressed genes in each functional category. Mean expression for No SCI is included as a comparison. UCSC Genome Browser tracks for the CONT SCI samples positioned at Calca (CGRPα) are shown in Fig. 5. Gene expression appears to be consistent for both the transcript variants at all the locations across the six CONT SCI samples, suggesting strong intra-group consistency. potential batch effects. The period between injury and tissue collection varied from 8.5 to 13.5 weeks, depending on experimental condition. This was intentional and resulted from an attempt to balance the requirements of our experimental design and the well-being of the animals. The animals undergoing a complete spinal cord transection required extensive care to ensure their well-being. Given that physiological measures are known to stabilize by four weeks after complete transection 39 , it was deemed acceptable to collect tissue at 8.5 weeks, over 4 weeks beyond this "stabilized" time point. The exercised animals required an initial period of time for introduction to the exercise facility followed by 10 weeks of exercise training, resulting in the 13.5 weeks time point