Publications

Ischemic stroke triggers rapid transcriptional changes in the brain, including the induction of noncoding RNAs, which are well-established regulators of post-stroke pathophysiology. Among the numerous classes of noncoding RNAs, short interspersed nuclear element RNAs (SINE-RNAs) are transcribed by Pol III and reported to be upregulated in various paradigms of cellular stress. In the ischemic brain, Pol III-driven gene expression is not well-studied and the expression of SINE-RNAs is virtually unmapped. In the current study, we used a mouse model of transient focal ischemia to evaluate for the first time post-stroke SINE-RNA expression in the cerebral cortex on a genome-wide scale. We observed SINE-RNA induction as early as 0 to 3 h of reperfusion and peak expression at 6 h of reperfusion, with 335 SINE-RNAs induced at this timepoint as compared to sham controls. Many of these transcripts remained induced through later timepoints during the acute phase of reperfusion (24 h). Fluorescence in situ hybridization against the SINE-RNAs, combined with immunohistochemistry for cell-type markers, revealed that these RNAs are localized to the nuclei of post-ischemic neurons and microglia in the ipsilateral cortex and hippocampus in both males and females. Further, we found that SINE-RNA expression was recapitulated in vitro following oxygen-glucose deprivation in HT22 hippocampal neurons, showing that they are reproducibly expressed in neurons in both in vivo and in vitro models of ischemia. Together, this is the first study to map genome-wide SINE-RNA expression in the post-ischemic brain and reveals a new layer of the noncoding transcriptome that may play a role in the post-stroke pathophysiology.

BACKGROUND: Recent advances in high-throughput transcriptomics have revealed extensive gene expression changes during cerebral ischemia. However, the changes in 3-dimensional chromatin architecture and long-range chromatin looping between genes and distal regulatory elements that drive these gene expression changes remain virtually unexplored. In this study, we map the landscape of the dynamic alterations in topologically associating domains and chromatin loops in the ischemic cortex and evaluate their contributions to poststroke transcriptional changes.

METHODS: We used high-throughput chromatin conformation capture (Hi-C) to profile genome-wide changes of the 3-dimensional chromatin architecture in the cerebral cortex of adult mice following a 1-hour middle cerebral artery occlusion and 6 hours of reperfusion or sham treatment. We conducted RNA sequencing to identify the differentially expressed genes that are concomitantly altered in association with the stroke-induced topologically associating domains and loops.

RESULTS: We identified 293 altered topologically associating domains following stroke, harboring a total of 60 upregulated differentially expressed genes and 106 downregulated differentially expressed genes. Chromatin looping analysis identified 4323 differentially altered loops in response to stroke, of which 1583 had enriched interactions and 2741 had diminished interactions. These loci included previously unknown enhancer and silencer elements, which were linked to a total of 88 upregulated and 96 downregulated genes. Upregulated genes associated with altered topologically associating domains and loops were linked to endothelial and immune cell functions, suggesting a role for dynamic chromatin remodeling in the poststroke cerebrovascular and neuroimmune response. Conversely, downregulated genes were associated with neuronal and synaptic functions, suggesting a role in poststroke neuronal fate.

CONCLUSIONS: This study is the first to map changes in the 3-dimensional chromatin of the adult cerebral cortex following ischemic stroke. Our findings reveal novel regulatory elements directly associated with stroke-altered genes that may be involved in the regulation of these genes via dynamic changes in chromatin architecture and looping characteristics to fine-tune their expression.

Noncoding DNA undergoes widespread context-dependent transcription to produce noncoding RNAs. In recent decades, tremendous advances in genomics and transcriptomics have revealed important regulatory roles for noncoding DNA elements and the RNAs that they produce. Enhancers are one such element that are well-established drivers of gene expression changes in response to a variety of factors such as external stimuli, cellular responses, developmental cues, and disease states. They are known to act at long distances, interact with multiple target gene loci simultaneously, synergize with other enhancers, and associate with dynamic chromatin architectures to form a complex regulatory network. Recent advances in enhancer biology have revealed that upon activation, enhancers transcribe long noncoding RNAs, known as enhancer RNAs (eRNAs), that have been shown to play important roles in enhancer-mediated gene regulation and chromatin-modifying activities. In the brain, enhancer dysregulation and eRNA transcription has been reported in numerous disorders from acute injuries to chronic neurodegeneration. Because this is an emerging area, a comprehensive understanding of eRNA function has not yet been achieved in brain disorders; however, the findings to date have illuminated a role for eRNAs in activity-driven gene expression and phenotypic outcomes. In this review, we highlight the breadth of the current literature on eRNA biology in brain health and disease and discuss the challenges as well as focus areas and strategies for future in-depth research on eRNAs in brain health and disease.

Stressors such as hypoxia, hypothermia, and acute toxicity often result in widespread cell death. This study investigated the outcomes of Neuro-2a (N2a; mouse neuroblastoma) cells following a cryogenic storage failure that exposed them to a combination of these stressors over a period of approximately 24-30 hours. Remarkably, a small fraction of the cells survived the event, underwent a period of dormancy, and eventually recovered to a healthy state. To understand the underlying resilience mechanisms, we created a model to replicate the dewar failure event and examined changes in phenotype, transcriptomics, proteomics, and mitochondrial activity of the surviving cells during recovery. We found that the surviving cells initially displayed a stressed morphology with irregular membranes and a clustered apperance. They showed an increased expression of proteins related to DNA repair and chromatin modification pathways as well as heightened mitochondrial function shortly after the stress event. As recovery progressed, the stress-responsive pathways, mitochondrial activity, and growth rates normalized toward that of healthy controls, indicating a return to a stable baseline state. These findings suggest that an initial robust energetic state supports key stress-responsive and repair pathways at the early stages of recovery, facilitating cell survival and resiliency after extreme stress. This work provides valuable insights into cellular resilience mechanisms with potential implications for improving cell preservation and recovery in biomedical applications and developing therapeutic strategies for conditions involving cell damage and stress.

Enhancer-derived RNAs (eRNAs) are a new class of long noncoding RNA that have roles in modulating enhancer-mediated gene transcription, which ultimately influences phenotypic outcomes. We recently published the first study mapping genome-wide eRNA expression in the male mouse cortex during ischemic stroke and identified 77 eRNAs that were significantly altered following a 1 h middle cerebral artery occlusion (MCAO) and 6 h of reperfusion, as compared to sham controls. Knockdown of one such stroke-induced eRNA - eRNA_06347 - resulted in significantly larger infarcts, demonstrating a role for eRNA_06347 in modulating the post-stroke pathophysiology in males. In the current study, we applied quantitative real-time PCR to evaluate whether the 77 eRNAs identified in the male cortex also show altered expression in the post-stroke female cortex. Using age-matched and time-matched female mice, we found that only a subset of the 77 eRNAs were detected in the post-stroke female cortex. Of these, only a small fraction showed similar temporal expression characteristics as males, including eRNA_06347 which was highly induced in both sexes. Knockdown of eRNA_06347 in the female cortex resulted in significantly increased infarct volumes that were closely matched to those in males, indicating that eRNA_06347 modulates the post-stroke pathophysiology similarly in males and females. This suggests a common underlying role for eRNA_06347 in the two sexes. Overall, this is the first study to evaluate eRNA expression and perturbation in the female cortex during stroke, and present a comparative analysis between males and females. Our findings show that eRNAs have sex-dependent and sex-independent expression patterns that may be of significance to the pathophysiological responses to stroke in the two sexes.

Recent studies have reported widespread stimulus-dependent transcription of mammalian enhancers into noncoding enhancer RNAs (eRNAs), some of which have central roles in the enhancer-mediated induction of target genes and modulation of phenotypic outcomes during development and disease. In cerebral ischemia, the expression and functions of eRNAs are virtually unknown. Here, we applied genome-wide H3K27ac ChIP-seq and genome-wide RNA-seq to identify enhancer elements and stroke-induced eRNAs, respectively, in the mouse cerebral cortex during transient focal ischemia. Following a 1-h middle cerebral artery occlusion (MCAO) and 6 h of reperfusion, we identified 77 eRNAs that were significantly upregulated in stroke as compared to sham, of which 55 were exclusively expressed in stroke. The knockdown of two stroke-induced eRNAs in the mouse brain resulted in significantly larger infarct volumes as compared to controls, suggesting that these eRNAs are involved in the post-stroke neuroprotective response. A preliminary comparison of eRNA expression in the male versus female cortices revealed sex-dependent patterns that may underlie the physiological differences in response to stroke between the two sexes. Together, this study is the first to illuminate the eRNA landscape in the post-stroke cortex and demonstrate the significance of an eRNA in modulating post-stroke cortical brain damage.

Gene expression in cerebral ischemia has been a subject of intense investigations for several years. Studies utilizing probe-based high-throughput methodologies such as microarrays have contributed significantly to our existing knowledge but lacked the capacity to dissect the transcriptome in detail. Genome-wide RNA-sequencing (RNA-seq) enables comprehensive examinations of transcriptomes for attributes such as strandedness, alternative splicing, alternative transcription start/stop sites, and sequence composition, thus providing a very detailed account of gene expression. Leveraging this capability, we conducted an in-depth, genome-wide evaluation of the protein-coding transcriptome of the adult mouse cortex after transient focal ischemia at 6, 12, or 24 h of reperfusion using RNA-seq. We identified a total of 1007 transcripts at 6 h, 1878 transcripts at 12 h, and 1618 transcripts at 24 h of reperfusion that were significantly altered as compared to sham controls. With isoform-level resolution, we identified 23 splice variants arising from 23 genes that were novel mRNA isoforms. For a subset of genes, we detected reperfusion time-point-dependent splice isoform switching, indicating an expression and/or functional switch for these genes. Finally, for 286 genes across all three reperfusion time-points, we discovered multiple, distinct, simultaneously expressed and differentially altered isoforms per gene that were generated via alternative transcription start/stop sites. Of these, 165 isoforms derived from 109 genes were novel mRNAs. Together, our data unravel the protein-coding transcriptome of the cerebral cortex at an unprecedented depth to provide several new insights into the flexibility and complexity of stroke-related gene transcription and transcript organization.

Ischemic stroke is an acute brain injury with high mortality and disability rates worldwide. The pathophysiological effects of ischemic stroke are driven by a multitude of complex molecular and cellular interactions that ultimately result in brain damage and neurological dysfunction. The Human Genome Project revealed that the vast majority of the human genome (and mammalian genome in general) is transcribed into noncoding RNAs. These RNAs have several important roles in the molecular biology of the cell. Of these, the long noncoding RNAs are gaining particular importance in stroke biology. High-throughput analysis of gene expression using methodologies such as RNA-seq and microarrays have identified a number of aberrantly expressed lncRNAs in the post-stroke brain and blood in experimental models as well as in clinical samples. These expression changes exhibited distinct temporal and cell-type-dependent patterns. Many of these lncRNAs were shown to modulate molecular pathways that resulted in deleterious as well as neuroprotective outcomes in the post-stroke brain. In this review, we consolidate the latest data from the literature that elucidate the roles and functions of lncRNAs in ischemic stroke. We also summarize clinical studies identifying differential lncRNA expression changes between stroke patients and healthy individuals, and genetic variations in lncRNA loci that are correlated with an increased risk of stroke development.

Long noncoding RNAs (lncRNAs) play major roles in regulating gene expression in mammals, but are poorly understood in ischemic stroke. Using a mouse model of transient focal ischemia, we applied RNA-seq to evaluate for the first time the unbiased, genome-wide expression of lncRNAs as a function of reperfusion time in the cerebral cortex. Focal ischemia was induced in adult male C57BL/6 mice followed by reperfusion for 6, 12 or 24h. Total RNA from ipsilateral cortices was used for Illumina sequencing and reads were mapped to the mouse reference genome (GRCm38). Annotated and novel transcript isoforms were identified and differential expression between the groups was estimated. We observed that the baseline expression of lncRNAs in the healthy cortex was low, but many were highly altered after stroke. Very few of these altered lncRNAs were previously annotated. A total of 259 lncRNA isoforms at 6h, 378 isoforms at 12h, and 217 isoforms at 24h of reperfusion were differentially expressed versus sham controls. Of these, 213, 322 and 171 isoforms at 6, 12 and 24h of reperfusion, respectively, were novel lncRNAs. Reperfusion time-point-specific analyses revealed that the lncRNAs reached peak expression levels at 6h of reperfusion. Positional analysis of ischemia-responsive lncRNAs with respect to ischemia-responsive protein-coding genes identified potential gene-regulatory relationships. Overall, this work shows that transient focal ischemia induces widespread changes in the expression of lncRNAs in the mouse cortex with distinct reperfusion time-point-dependent expression characteristics that may underlie progression of the ischemic pathophysiology. The detection of hundreds of novel ischemia-responsive lncRNAs marks the discovery of new disease-related genomic regions in the adult cortex and may help identify novel opportunities for therapeutic targeting.

MicroRNAs (miRNAs) are small non-coding RNAs that are known to control mRNA translation. Most miRNAs are transcribed from specific genes with well-defined promoters located throughout the genome. The mechanisms that control miRNA expression under normal and pathological conditions are not yet understood clearly. Peroxisome proliferator-activated receptor (PPAR) γ is a ligand-activated transcription factor that is extensively distributed in the CNS. PPARγ activation induces neuroprotection by modulating genes that contain peroxisome proliferator response elements (PPREs) in their promoters. We presently evaluated if PPARγ modulates miRNA expression. When adult rats were treated with PPARγ agonist rosiglitazone, expression of 28 miRNAs altered significantly (12 up- and 16 down-regulated; 3-119 fold) in the cerebral cortex compared to vehicle-treated controls. In silico analysis showed 1-5 PPREs in the putative promoter regions (within 1 Kb upstream of the transcription start site) of these miRNA genes. Cotransfection with a PPARγ constitutively expressing vector significantly induced the miR-145 and miR-329 promoter vectors (each have four PPREs), which was curtailed by point mutations of PPREs in their promoters. Interestingly, the PPARγ promoter has binding sites for both these miRNAs and transfection with miR-329 mimic and miR-145 mimic induced the PPARγ expression. Thus, these studies show a cyclical induction of miRNAs and PPARγ, indicating that the pleiotropic beneficial effects of PPARγ agonists might be modulated in part by miRNAs and their down-stream mRNAs. We proposed that promoters of many microRNAs contain the binding sites for the transcription factor PPARγ. Activation of PPARγ modulates the expression of these microRNAs. Two such PPARγ-responsive microRNAs (miR-145 and miR-329) bind to PPARγ promoter to induce its expression. This indicates the presence of a feedback loop by which transcription factors and microRNAs can modulate each other.