Compartment-specific total RNA profile of Hippocampal and Cortical cells from Mesial Temporal Lobe Epilepsy tissue

Human tissue

Human pharmacoresistant mTLE patient hippocampus (HC) and neocortical (Cx) tissue was obtained after surgery, as described previously (Kan et al., 2012; Vangoor et al., 2019) at the University Medical Center Utrecht. Informed consent was obtained from all patients for use of tissue and clinical information for research purposes and all procedures were approved by the Review board of University Medical Center Utrecht (Utrecht, The Netherlands). Postmortem human HC and Cx tissue was obtained from the Netherlands Brain Bank (Amsterdam, The Netherlands). A written informed consent was obtained from all donors for brain autopsy and for use of the material and clinical information for research purposes. Use of postmortem tissue and clinical information for research purposes was approved by the medical ethics board of the Amsterdam University Medical Center (Amsterdam, The Netherlands). For HC, a total of five HC mTLE (no hippocampal sclerosis (mTLE non-HS) (male and female)) and five control HC (male and female) samples were used. For Cx, six mTLE (non-HS (male and female) and six HS International League Against Epilepsy (ILAE) Type 1 (male and female) or mTLE+HS) and six control samples were used (Supplementary Table S1).

Nuclear-cytoplasmic fractionation of human hippocampus and cortex

Before fractionation, frozen brain samples were powdered in dry ice using a mortar pestle. Approximately 500 mg of powdered tissue was resuspended in Hypotonic Lysis Buffer (HLB) (10 mM Tris (pH 7.4), 3 mM CaCl₂, 2 mM MgCl₂, 1% Nonidet P-40 (Sigma), Protease inhibitor (Roche), 60 U Superase-In/mL (ThermoFisher Scientific)) for 10 minutes (min) on ice to lyse: homogenization and lysis were enhanced with the help of 30-80 strokes in a tightly fitting glass Douncer. The cytoplasmic fraction was collected following a gentle centrifugation (500 g) for 5 min at 4°C, and frozen in Qiazol (Qiagen) until RNA purification. The semi-pure nuclear fraction (pellet) was resuspended in 2 mL Sucrose Buffer I (0.32 M Sucrose, 5 mM CaCl₂, 3 mM Mg(CH₃COO)₂, 0.1 mM EDTA, 10 mM Tris (pH 8.0), 0.5% Nonidet P-40, 1 mM DTT, 25 U Superase-In/mL) and placed on top of two thicker sucrose cushions: Sucrose Buffer II (1.6 M Sucrose, 3 mM Mg(CH₃COO)₂, 0.1 mM EDTA, 10 mM Tris (pH 8.0), 1 mM DTT, Protease Inhibitor, 35 U Superase-In/mL), and a Sucrose Buffer III (2.0 M Sucrose, 3 mM Mg(CH₃COO)₂, 0.1 mM EDTA, 10 mM Tris (pH 8.0), 1 mM DTT, 35 U Superase-In/mL, Protease Inhibitor), respectively, in ultracentrifugation tubes (Beckman Coulter). Samples were centrifuged using a SW41.4 rotor at 30000 g at 4°C for 45 min. Pure nuclei were collected and stored in Qiazol (Qiagen) at −80°C until RNA extraction was performed.

RNA purification

Total RNA was extracted from fractionated human samples using miRNeasy Mini Kit (Qiagen) according to the manufacturer’s instructions and DNaseI treated using RNase free-DNase (Catalogue#79254, Qiagen) to eliminate genomic DNA contamination. When required, RNA was further purified using RNA Clean-Up and Concentration Kit (Catalogue#23600, Norgen Biotek Corporation). Finally, RNA integrity (RIN) value was estimated in an Agilent 2100 Bioanalyzer (Agilent Technologies) (Supplementary Table S1).

Quality control fractionation procedure

The efficiency of the fractionation protocol was estimated at the protein level by western blotting (WB) and at the RNA level by quantitative PCR (RT-qPCR). For WB analysis, a fraction of cytoplasmic and nuclear samples were collected in RIPA buffer (50 mM Tris, pH.7.5, 150 mM NaCl, 0.5% NP-40, 0.5% NaDoc, 1% Triton, Protease inhibitors in MilliQ). Nuclear sample was lysed by passing through 27 g needle (BD Plastipak) for five times and by sonication. 12 µl of the cytosol and 6 µl of nuclear protein samples were resuspended in 4x NuPAGE LDS Sample buffer (ThermoFisher Scientific) and heated at 95°C for 5 min. The samples were separated in SDS-PAGE gels (8%) for approximately 1 h at 160 V, until the dye ran out of the gel. Separated proteins were transferred onto nitrocellulose blotting membrane (GE Healthcare) for 1 h at 100 V under cold conditions. Then blots were blocked for 1 h at RT in 5% milk powder in 1xTBS-Tween. Blots were incubated overnight at 4°C with antibodies against the nuclear-enriched protein fibrillarin (1:1000, Catalogue# 2639S, Cell signalling) or the cytoplasmic-enriched protein β-tubulin (1:2000, Catalogue# T5168, Sigma) in blocking buffer. Blots were stained with peroxidase-conjugated secondary antibodies (1:30000 for HRP-anti-mouse; 1:50000 for HRP-anti-rabbit) for 1 h at RT in 1xTBS-Tween, and signal was detected by incubating blots with Pierce ECL substrate (Thermo Fisher Scientific). Images were acquired using a FluorChem M imaging system (Protein Simple).

RNA levels of nuclear-enriched markers (NEAT1, GOMAFU, pre-GAPDH) and cytoplasmic-enriched markers designed to amplify an intron spanning amplicon (GAPDH) were measured by RT-qPCR, both in cytoplasmic and nuclear preparations. Briefly, 100 ng of total RNA reaction was reverse transcribed using the SuperScript III First-Strand Synthesis System (ThermoFisher Scientific) with random hexamers, according to the manufacturer’s instructions. RT-qPCR quantification of mRNA was conducted using FastStart Universal SYBR Green Master (Rox) (Roche) for all samples, in triplicates, in a QuantStudio(tm) 6 Flex Real-Time PCR System (ThermoFisher Scientific). Real-time reactions were carried out as follows: pre-denaturation at 95°C for 10 min, followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. Melting curves were generated for the final PCR products by decreasing the temperature to 60°C for 1 min followed by an increase in temperature to 95°C. Primer sequences used for RT–qPCR reactions are listed in Supplementary Table S3.

Total RNA library preparation and sequencing

A total of 2 µg of RNA per sample in a 5 µl volume was rRNA-depleted using the Ribo-Zero Magnetic Kit (human/mouse/rat; Illumina) before library preparation. Stranded, paired-end sequencing libraries were prepared using ScriptSeq v2 (Epicentre), after which samples were subjected to quality control in the 2100 Bioanalyzer (Agilent). Sequencing was performed as paired end 100 bp reads on an Illumina HiSeq 4000 sequencer at the Beijing Genomics Institute (BGI), Hongkong, China.

Read mapping and differential expression analysis

Read mapping and differential expression analysis was performed as described previously (Vangoor et al., 2019). Briefly, following trimming of low-quality bases and adapter sequences with FASTQ-MCF (version 0.0.13), processed reads were mapped to the hg19/GCh37 reference human genome (iGenomes) with TopHat2 (version 2.0.13) (Kim et al., 2013). For the alignments of the total RNA-seq data, the ‘fr-secondstrand’ option was chosen. Mapped counts for each gene were summarised using the python script htseq-count (Anders et al., 2015). Alignment statistics were obtained from Bam files generated using Bowtie (Langmead et al., 2009). All raw and processed RNA-seq data are deposited at NCBI Gene Expression Omnibus (GEO) with reference number GSEXXXX.

Mapped counts were visualized using a principal component analysis (PCA; R version 4.0.3, package ggfortify). For differential expression analysis, genes and transcripts count data was analysed for differential expression in R using the Bioconductor package EdgeR version 3.12.1 (Robinson et al., 2010), with the trimmed mean of M-values (TMM) normalisation method (Robinson and Oshlack, 2010). Gene expression levels were corrected for batch and gender effects by including the series of sequencing rounds, sexes of controls and patients to a generalised linear model. Using the Benjamini-Hochberg false discovery rate (FDR) adjusted P-values for multiple testing were calculated and only genes with an FDR < 0.05 were considered significantly differentially expressed (DE). Visualisation of data was performed in R using the ggplot2 library (version 2.1.0). Using the Bioconductor heatmap package in R, gene expression heatmaps with hierarchical clustering of expression profiles were created. To estimate proportion of protein coding RNAs, total RNA-seq reads were pseudoaligned to the human protein coding transcriptome (Homo_sapiens.GRCh38.cds.all.fa from the Ensembl FTP site; (Yates et al., 2019)) using the Kallisto algorithm (Bray et al., 2016). To set the number of bootstraps to 100, the extra argument -b 100 was passed to the function.

Validating the expression of differentially expressed (DE) genes by RT-qPCR

RNA from an independent set of mTLE-HS and control samples were used to validate DE genes. Reverse transcription (RT) was performed with the High-Capacity RNA-to-cDNA Kit (Catalogue# 4387406, Thermo Fisher) according to manufacturer’s instructions. Briefly, 2x RT Buffer Mix and 20x RT Enzyme Mix were mixed with RNA (20 ng) and nuclease-free water in a final reaction volume of 20 μl. Reactions were then incubated in a thermal cycler with the following program: 60 min at 37°C and 5 min at 95°C, with the infinite hold set at 4°C. For RT-qPCR, reaction mixes were prepared according to manufacturer’s instructions. Briefly, 1x reaction was prepared by mixing 2x Universal Master Mix (Catalogue # 4440040, Thermo Fisher), 20x Taqman Gene Expression Assay, and nuclease-free water in a final volume of 9 μl. Next, 1 μl of cDNA was added to the reaction mix. Finally, 384-well RT-qPCR plates were spun down for 5 min in a centrifuge at 5000 g. The RT-qPCR reactions were performed in a QuantStudioTM 6 Flex Real-Time PCR System, using the following program: 10 min at 95°C and 40x (15 s at 95°C + 1 min at 60°C). Pre-designed gene expression assays were obtained from ThermoFisher Scientific: BTG2 (Hs00198887_m1), IL1ß (Hs01555410_m1), FAM126A (Hs00261521_m1), WISP1 (Hs05047584_s1) and HLA-C (Hs05049205_g1). GAPDH (Hs03929097_g1) was used as endogenous control. Every sample was run in triplicate and fold changes were estimated using the 2^-ΔCt method. Statistical significance between conditions was estimated by performing Mann-Whitney tests using Graphpad Prism software (Version 9.2.0).

Pathway analysis using Reactome

Pathway enrichment analysis was performed using the Reactome pathway analysis tool (https://reactome.org), by applying the hypergeometric distribution for P value comparison (Fabregat et al., 2017). Analysis was performed by using the DE gene list containing 174 genes for the cytoplasm and 337 genes for the nuclear compartment of the HC, 99 genes for the cytoplasm and 138 genes for the nuclear compartment of Cx mTLE non-HS and 223 genes for the cytoplasm and 409 genes for the nuclear compartment of Cx mTLE+HS. Enrichment of pathways was considered significant with a corrected P value of < 0.05 (Benjamin-Hochberg).

Total RNA profile of individual cell compartments

To understand the subcellular RNA profile of human mTLE patient hippocampal and neo-cortical tissue, total RNA-seq was performed on 1) five mTLE non-HS HC samples and compared to five control HC, and 2) six mTLE non-HS, six mTLE+HS Cx samples were compared separately to six control Cx samples (Figure 1A). Nuclear and cytoplasmic fractions of human brain tissue from patients and controls were successfully isolated by lysing tissue in hypotonic lysis buffer followed with ultracentrifugation. The purity of nucleus samples had evidently improved after ultracentrifugation (Figure 1B). Fractionation success was examined at both the protein and RNA level. The cytoplasmic protein β-Tubulin was only observed in the cytoplasmic samples, whereas the nucleus-enriched protein Fibrillarin was observed in nuclear samples (Figure 1C). At the RNA level, the cytoplasmic marker gene GAPDH was enriched in cytoplasmic samples, whereas in nuclear samples several nucleus-specific lncRNAs were found (GOMAFU, MALAT1, NEAT1) in addition to GAPDH intronic transcripts (Figure 1D).

• Download figure
• Open in new tab

Figure 1: Nuclear fractionation and RNA-seq of mTLE brain tissue.

(A) Summary legend of the cohorts (cortical and hippocampal tissue) used for the RNA-seq experiment. For hippocampus (HC), a total of five mTLE non-HS HC samples and five control HC samples were used. For Cx, six mTLE non-HS, six mTLE+HS and six control samples were used. Similar color codes indicate that both HC and Cx tissue were obtained from the same individual control or mTLE non-HS patient.

(B) Overview of the cell fractionation procedure for the isolation of RNA from nuclear and cytoplasmic fractions. DAPI (4',6-diamidino-2-phenylindole). Scale bars = 20 µm.

(C) Western blot analysis of the lysates of the nuclear and cytoplasmic fractions. Fibrillarin and β-Tubulin were used as nuclear and cytoplasmic markers, respectively. Consistent with its reported subcellular localisation, fibrillarin was specifically detected in the nuclear fraction, whereas anti-β-Tubulin antibody showed signal only in cytoplasmic fractions. Nuc: nucleus, Cyt: cytoplasm.

(D) Bar plot showing the enrichment of transcripts in the nuclear and cytoplasmic fractions. As expected, RT-qPCR detected nuclear enrichment for GOMAFU, MALAT1, NEAT1 and GAPDH int (an intron retained variant of GAPDH). Higher levels of the spliced form of GAPDH were reported in the cytoplasmic fraction. Enrichment is represented as a percentage of the distribution (% distribution) of a marker across the nuclear and cytoplasmic compartments.

Total RNA-seq was performed and an average of 61 million high-quality reads per cytoplasmic and 65 million reads per nuclear samples were obtained for Cx samples. For HC samples, 51 and 58 million reads in cytoplasm and nuclear samples, respectively, were obtained. Of these, 59% reads were uniquely mapped to the reference human genome in cytoplasmic Cx samples, and 62% in the nuclear Cx samples. For HC samples, 63% of reads mapped to the genome in cytoplasmic fractions and 66% in nuclear samples. We observed that approximately 36% of reads from cytoplasmic samples were aligned to exons, whereas in the nuclear samples approximately 15% of reads were aligned to exons. We hypothesize that the lower level of alignment in nucleus could be due to the high level of intronic or intergenic sequences present in this compartment (Supplementary Table S2). Due to the presence of nascent transcripts, an increased alignment of nuclear samples with intronic and intergenic regions has also been shown in recent studies performing total RNA-seq from human brain nuclear and cytosolic fractions (Zaghlool et al., 2021). Principal component analysis (PCA) showed a clear separation between cytoplasmic and nuclear samples, indicating differential RNA profiles between different cell compartments (Figure 2A, D, E). It is, however, important to note that the separation of control and mTLE samples was less clear (for both compartments (nuclear and cytoplasmic) and tissue types) (HC – Figure 2B, 2C, Cx – Figure 2D, 2E).

• Download figure
• Open in new tab

Figure 2: Principal component analysis (PCA) plots of mRNA raw count data.

PCA plots showing variance captured by the first two PCs along the respective axes.

(A) Cytoplasmic and nuclear samples are clustered separately for HC control versus mTLE non-HS group. However, individual comparison of HC control versus mTLE non-HS group for each compartment did not result in separate clusters either in Cyt (Cytoplasmic) (B) or Nuc (Nuclear) samples (C). For Cx, separate clustering of cytoplasmic and nuclear samples was observed for control versus mTLE non-HS (D) and control versus mTLE+HS groups (E). (F) Similarly, pre- and postmortem treatment groups were separated on basis of compartments.

Postmortem effect on gene expression

Postmortem delay is a common feature of control tissue samples obtained from brain banks (Ferrer et al., 2008). It is often not clear to which extent postmortem delay in control tissues will affect transcript levels and affect differential gene expression analysis when compared with samples collected without any postmortem delay. To estimate this potential confounding effect, resected Cx tissue obtained directly from the surgery room was subjected to postmortem delay treatment. Briefly, tissue was cut into two halves, where one half (pre-) was stored immediately on dry ice and then at −80°C until processing. The second half (post-) was incubated at 4°C in saline for 24 h to mimic postmortem conditions observed for tissue samples collected from brain banks. Cytoplasmic and nuclear fractions were collected from the pre- and postmortem samples as described above. RNA-seq of the pre- and postmortem samples did not show strong differences in the overall coverage and mapping of reads (Supplementary Table S2). PCA plots of the pre- and postmortem-treated samples showed an overall similar RNA profile in cytoplasmic and nuclear fractions (Figure 2F). Finally, no major changes in RNA expression were found in the same sample before and after postmortem delay, either in cytoplasmic or nuclear fractions (Supplemental Figure S1A). In comparison, significant differences in gene expression were observed in Cx samples between control and mTLE+HS or mTLE non-HS (Supplemental Figure S1B-S1C). Together, these data show that postmortem delay is unlikely to cause major differences in the transcriptomic profiles of control tissue.

• Download figure
• Open in new tab

Supplementary Figure S1:

(A) Volcano plots showing overall gene expression changes observed between Pre- versus Postmortem-treated nuclear and cytoplasmic samples. (B) Volcano plots for Cx Control versus mTLE+HS Cytoplasm and Nucleus. (C) Volcano plots for Cx Controls versus mTLE non-HS Cytoplasm and nucleus. Blue colored dots indicate significantly downregulated genes and the red colored dots indicate significantly upregulated genes.

Differential gene expression and pathway analysis

Several transcripts were found to be differentially expressed in cytoplasmic and nuclear compartments (Figure 3A-B). Several of the differentially expressed genes (DEGs) were unique to a given compartment, whereas other transcripts were altered in both the nuclear and cytoplasmic compartments, in mTLE non-HS HC nucleus and cytoplasm (Figure 3C), mTLE non-HS Cx nucleus and cytoplasm (Figure 3D) and mTLE+HS Cx nucleus and cytoplasm (Figure 3E). A few transcripts commonly upregulated in hippocampal cytoplasmic and nuclear samples (IL1β, FAM126A, BTG2, WISP1) were confirmed using RT-qPCR in an independent set of unfractionated patient and control samples. A significant upregulation of the selected transcripts in mTLE non-HS samples compared to control samples was found, confirming the RNA-seq analysis: IL1β (FC = 12.09, n = 7, p = 0.0011), FAM126A (FC = 2.00, n = 7, p = 0.0011), BTG2 (FC = 3.05, n = 7, p = 0.0005), WISP1 (FC = 6.08, n = 7, p = 0.004) (Figure 3F). In contrast, HLA-C was specifically downregulated in the nuclear compartment by RNA-seq analysis, but was unchanged using RT-qPCR analysis of unfractionated tissue; HLA-C (FC = 1.14, n = 7, p = 0.7103) (Figure 3F).

• Download figure
• Open in new tab

Figure 3: Differentially expressed genes and experimental validation.

Heatmaps of differentially expressed mRNAs across the 10 hippocampal samples (five controls and five mTLE non-HS patients) for cytoplasmic (A) and nuclear (B) fractions. Hierarchical clustering is shown at the top. HC: hippocampus, Controls, mTLE non-HS, Cyt: cytoplasm, Nuc: nucleus. Venn diagrams showing overlap of the number of differentially expressed genes (DEGs) between nuclear (Nuc) and cytoplasmic (Cyt) compartments in HC mTLE non-HS (C), Cx mTLE non-HS (D) and Cx mTLE+HS (E). Color code: red indicates upregulated genes and blue indicates downregulated genes. (F) RT-qPCR validation of a selection of DEGs in an independent set of HC samples, comparing expression levels between controls and mTLE non-HS. GAPDH was used for normalization. n = 7 mTLE non-HS and 7 control samples. Each colored dot is a sample. Error bars represent standard error of mean. Mann-Whitney tests. ***p < 0.001, **p < 0.01, n.s. non-significant.

In order to identify pathways that are enriched in the cytoplasm and nucleus based on DEGs identified by RNA-seq, reactome pathway enrichment analysis was performed. Pathways enriched in HC mTLE non-HS cytoplasmic samples were related to Neuronal system (Upregulated transcripts, entities = 17 out of 487, p = 5.22e-05, FDR = 0.021), Adrenoceptors (Upregulated transcripts, entities = 3 out of 11, p = 3.01e-04, FDR = 0.061), Transmission across Chemical Synapses (Upregulated transcripts, entities = 12 out of 341, p = 6.35e-04, FDR = 0.065), GABA receptor activation (Upregulated transcripts, entities = 5 out of 67, p = 0.001, FDR = 0.092), Potassium channels (Upregulated transcripts, entities = 6 out of 107, p = 0.002, FDR = 0.092), Activation of AMPA receptor (Upregulated transcripts, entities = 2 out of 7, p = 0.003, FDR = 0.127), Caspase mediated cleavage of cytoskeletal proteins (Downregulated transcripts, entities = 3 out of 12, p = 3.61e-06, FDR = 5.12e-04) and Apoptotic cleavage of cellular proteins (Downregulated transcripts, entities = 3 out of 38, p = 1.10e-04, FDR = 0.005). Pathways derived from nuclear samples included NGF-stimulated transcription (Upregulated transcripts, entities = 8 out of 56, p = 1.93e-05, FDR = 0.016), Signaling by Tyrosine kinase receptors (NTRK1) (Upregulated transcripts, entities = 9 out of 143, p = 0.002, FDR = 0.542), Endosomal/vacuolar pathway (Downregulated transcripts, entities = 14 out of 82, p = 4.11e-15, FDR = 1.27e-12), Antigen presentation (Downregulated transcripts, entities = 14 out of 102, p = 7.76e-14, FDR = 1.20e-11), and Interferon alpha/beta signaling (Downregulated transcripts, entities = 15 out of 186, p = 1.72e-11, FDR = 1.77e-09) (Supplementary file 2). No common pathways were found when comparing HC and Cx mTLE non-HS conditions. Importantly, pathways identified for Cx mTLE+HS in nuclear and cytoplasmic samples were commonly enriched for Interleukin-4 and Interleukin-13 signaling (Nuc: entities = 18 out of 211, p = 1.96e-04, FDR = 0.087. Cyt: entities = 16 out of 211, p = 1.88e-06, FDR = 0.001) hinting at activation of immune response genes which is a pathological signature of mTLE+HS tissue (Supplementary file 3).