ORIGINAL RESEARCH article

Front. Mol. Neurosci., 28 April 2023

Sec. Brain Disease Mechanisms

Volume 16 - 2023 | https://doi.org/10.3389/fnmol.2023.1166875

mRNA and lncRNA co-expression network in mice of acute intracerebral hemorrhage

  • 1. Department of Integrated Traditional Chinese and Western Medicine, Institute of Integrative Medicine, Xiangya Hospital, Central South University, Changsha, Hunan, China

  • 2. National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, Hunan, China

Article metrics

View details

13

Citations

4,2k

Views

1,3k

Downloads

Abstract

Background:

Intracerebral hemorrhage (ICH) is a severe subtype of stroke lacking effective pharmacological targets. Long noncoding RNA (lncRNA) has been confirmed to participate in the pathophysiological progress of various neurological disorders. However, how lncRNA affects ICH outcomes in the acute phase is not completely clear. In this study, we aimed to reveal the relationship of lncRNA-miRNA-mRNA following ICH.

Method:

We conducted the autologous blood injection ICH model and extracted total RNAs on day 7. Microarray scanning was used to obtain mRNA and lncRNA profiles, which were validated by RT-qPCR. GO/KEGG analysis of differentially expressed mRNAs was performed using the Metascape platform. We calculated the Pearson correlation coefficients (PCCs) of lncRNA-mRNA for co-expression network construction. A competitive endogenous (Ce-RNA) network was established based on DIANALncBase and miRDB database. Finally, the Ce-RNA network was visualized and analyzed by Cytoscape.

Results:

In total, 570 differentially expressed mRNAs and 313 differentially expressed lncRNAs were identified (FC ≥ 2 and value of p <0.05). The function of differentially expressed mRNAs was mainly enriched in immune response, inflammation, apoptosis, ferroptosis, and other typical pathways. The lncRNA-mRNA co-expression network contained 57 nodes (21 lncRNAs and 36 mRNAs) and 38 lncRNA-mRNA pairs. The ce-RNA network was generated with 303 nodes (29 lncRNAs, 163 mRNAs, and 111 miRNAs) and 906 edges. Three hub clusters were selected to indicate the most significant lncRNA-miRNA-mRNA interactions.

Conclusion:

Our study suggests that the top differentially expressed RNA molecules may be the biomarker of acute ICH. Furthermore, the hub lncRNA-mRNA pairs and lncRNA-miRNA-mRNA correlations may provide new clues for ICH treatment.

1. Introduction

Stroke remains one of the leading causes of death worldwide and is associated with an extremely high rate of disability (Organization, G.W.H., 2022). Intracerebral hemorrhage (ICH) accounts for approximately 27.9% of all strokes, with a higher mortality than that of ischemic stroke (GBD 2019 Stroke Collaborators, 2021). Especially the mortality is as high as 30–40% in acute and subacute phases (Greenberg et al., 2022). Although surgical therapies for evacuating hematoma and alleviating the compressive effect reduce the mortality in ICH acute phase, there is no significant benefit for the long-term prognosis (Zhao et al., 2020). Apart from the drugs targeting the risk factors, there is still a lack of effective drugs for ICH patients. Therefore, the novel therapeutic targets for acute ICH need further exploration.

Long noncoding RNA (lncRNA) is a type of non-protein-coding RNA with transcripts of more than 200 nucleotides (Mattick et al., 2023). According to the lncRNA location with protein-coding genes (PCGs), lncRNA is classified into intergenic lncRNA, antisense lncRNA, and intronic lncRNA (Jarroux et al., 2017). LncRNA evolves rapidly and rarely has an ortholog in another species (Mattick et al., 2023). Although their poor conservation has confounded efforts to predict lncRNA functions across species, it has also highlighted potential functional areas (Ransohoff et al., 2018). Over the last decade, increasing attention has been paid to its functions. LncRNA contributes to translational regulation, epigenetic modifications, post-translational modifications, and protein/mRNA stability by interacting with DNA, RNA, and protein (Bridges et al., 2021). Recent studies have even found that some lncRNA encodes small peptides less than 100 aa (Choi et al., 2019). The lncRNA functions are associated with their subcellular localization and transcriptional loci relative to the modulated gene (Kopp and Mendell, 2018). Based on the biological functions previously described, lncRNA participates in multiple neurological disorders development. Accumulating evidence has implicated lncRNA dysregulation in neurodegenerative disorders, including Alzheimer’s disease (Canseco-Rodriguez et al., 2022), Parkinson’s disease (Na et al., 2022), and Huntington’s disease (Tan et al., 2021).

In ischemic stroke, some lncRNAs play the role of anti-inflammation and/or anti-apoptosis, such as MALAT1 (Zhang et al., 2017) and Nespas (Deng et al., 2019). Whereas others produce an effect of pro-inflammation and pro-inflammatory factors releasing, such as NEAT1 (Ni et al., 2020) and Gm4419 (Wen et al., 2017). In addition, lncRNA H19 (Wang et al., 2017), MEG3 (Li et al., 2020), and RMST (Sun et al., 2019) have been proven to regulate M1 microglia polarization. Although there are currently some studies on the mechanism of lncRNA in stroke, most of them merely focus on ischemic stroke. So, the mechanisms of biological functions induced by lncRNA in ICH remain scarce.

lncRNA regulation of gene transcription is also influenced by other non-coding RNA, such as microRNA (miRNA). LncRNA may competitively bind to miRNA to conduct its functions, termed competitive endogenous RNA (ce-RNA). As the sponges of miRNA, lncRNA negatively regulates the target transcripts with similar miRNA response elements (MREs). In ischemic stroke, SNHG12/LNC04080 promotes angiogenesis and blocks inflammation by interacting with miR-199a (Long et al., 2018). MALAT1 targets miR-375 and miR-145 to aggravate neurological injury through increasing reactive oxygen species (ROS) and inhibiting Aquaporin 4 (AQP4), respectively (Zhang et al., 2017; Wang et al., 2020). Currently, there is little focus on the ce-RNA in acute ICH. It is suggested that lncRNA Mtss1 promotes secondary inflammation along with miR-709 (Chen et al., 2020). To dig into the underlying relationship among mRNA, lncRNA, and miRNA post-ICH, a comprehensive ce-RNA network should be revealed.

In this study, we examined the differentially expressed lncRNA and mRNA profiles in the acute phase of autologous blood-induced ICH animals. After validation of microarray results, we detected the functions of the differentially expressed mRNAs through GO/KEGG enrichment analysis. Furthermore, the lncRNA-mRNA co-expression network was established according to Pearson correlation coefficients. Moreover, to investigate the complicated correlation and interaction between lncRNA and mRNA, the ce-RNA network was carried out. The flow chart of this research is shown in Figure 1.

Figure 1

Figure 1

The flow chart of this study.

2. Materials and methods

2.1. Animals and ICH model

Studies were carried out on 8-week-old male C57BL/6 J mice obtained from the Laboratory Animal Centre of Central South University (CSU) (No. 2020sydw0929). All animals were housed in five per cage under stable environmental conditions (room temperature at approximately 25°C, relative humidity of the room at approximately 50%, 12 h light–dark cycle) with access to full water and food. All experimental procedures were approved by the Institutional Animal Care and Use Committee of CSU.

The ICH model was induced using autologous blood injection according to the published protocol (Rynkowski et al., 2008). Mice were randomly divided into sham and ICH groups. They were anesthetized by 0.3% pentobarbital sodium (22 mL/kg, i.p.). Autologous blood was obtained from the distal end of the mouse tail using a microinjector. For ICH induction, 15 μl fresh autologous blood was injected into the right globus pallidus at a rate of 1.5 μL min − 1. After injecting, the syringe was left in place for an additional 20 min to prevent backflow of autologous blood. The injection coordinates were 2 mm to the right, 0.5 mm posterior to Bregma, and 4 mm below the cortical surface. The sham group was injected with the same volume of 0.9% saline instead of autologous blood in the same site.

2.2. Behavioral tests

Behavioral tests were conducted on days 1, 3, and 7 after ICH induction. The modified neurological severity score (mNSS) and foot-fault test were used to evaluate the neurological injury (Li et al., 2019). In brief, mNSS assesses motor, sensory, and reflex abilities, and the foot-fault test evaluates the forelimbs dysfunction of mice. A higher score of mNSS indicates a more severe neurological function deficit. Meanwhile, data from the foot-fault test were presented as a percentage of foot faults (number of foot faults for left forelimb) per total number of steps.

2.3. Histological staining and quantify

The 3 to 4-μm coronal paraffin brain sections were stained following the instructions of the H&E (Servicebio, G1001&G1004) and Nissl’s staining kits (Servicebio, G1032), respectively. The cross-sectional hematoma areas were measured by ImageJ software in reference to a published study (Fei et al., 2022). To calculate the survival neurons, four peri-hematoma fields were photographed from each section at 200 magnifications. Within each filed, a fixed-size hematoma-free area was selected to count survival neurons. The number of survival neurons in the fixed-size area was calculated for statistical analysis.

2.4. mRNA and lncRNA microarray

Ten globus pallidus tissues (five from the sham group and five from the ICH group) were used to extract total RNA to detect microarray. The quantity and quality of total RNA were validated by NanoDrop ND-1000. Before microarray scanning, total RNA was purified by removing rRNAs (mRNA-ONLY™ Eukaryotic mRNA Isolation Kit, Epicenter) to obtain mRNAs. Next, the samples were amplified and transcribed into fluorescent cRNA (Arraystar Flash RNA Labeling Kit, Arraystar) and purified by RNeasy Mini Kit (Qiagen) subsequently. Arraystar Mouse lncRNA Microarray (V4.0, including about 22,692 mRNAs and 37,949 lncRNAs) was designed to profile mRNA along with lncRNA systematically. Agilent DNA Microarray Scanner (part number G2505C) was used for hybridized array scanning. The acquired array images were analyzed by Agilent Feature Extraction software (version 11.0.1.1). GeneSpring GX v12.1 software package (Agilent Technologies) performed quantile normalization and subsequent data processing. Differential expression of mRNAs and lncRNAs between the two groups were identified by applying a filter of fold change (FC) ≥ 2 and value of p <0.05.

2.5. LncRNA homology query

The sequence of lncRNA was obtained by referring to the NCBI and UCSC databases. Then, the sequence was aligned to the genome of human by the BLAT section of the UCSC database. The nucleotide blast tool of NCBI was used to validate the homology of lncRNA between mouse and human.

2.6. Real-time quantitative PCR

According to the manufacturer’s instructions, total RNA was extracted from globus pallidus tissue using a Trizol reagent (Takara, 9,108). Reverse transcription of mRNA/lncRNA used HiScript II Q RT SuperMix for qPCR (Vazyme, R223) as instructed by the manufacturer. Differential expression of mRNAs/lncRNAs was determined by RT-qPCR with ChamQTM Universal SYBR® qPCR Master Mix (Vazyme, Q711) following the manufacturer’s instructions, on a CFX Connect™ Real-time system (BIO-RAD, United States). Briefly, the cycling settings were as follows: incubation at 95°C for 5 min, followed by 40 cycles of 95°C for 10 s, 60°C for 30 s, and the instrument default melting curve acquisition procedures. GAPDH was used as the standardized control. Specific primers are listed in Table 1. The results of RT-qPCR were analyzed by the 2−ΔΔCt method.

Table 1

Name Primers (forward and primer)
GAPDH F: 5’ GGCAAATTCAACGGCACAGTCAAG 3′
R: 5’ TCGCTCCTGGAAGATGGTGATGG 3′
Saa3 F: 5’ AGCTGACCAGTTTGCCAATGAGTG 3′
R: 5’ ACCCAGTAGTTGCTCCTCTTCTCG 3′
Col10a1 F: 5’ ATGCCGCTTGTCAGTGCTAACC 3′
R: 5’ GGGTCGTAATGCTGCTGCCTATTG 3′
Mmp13 F: 5’ ACAGTTGACAGGCTCCGAGAAATG 3’
R: 5’ CCACATCAGGCACTCCACATCTTG 3’
Oxt F: 5’ ACCATCACCTACAGCGGATCTCAG 3’
F: 5’ GTCAGAGCCAGTAAGCCAAGCAG 3’
Pchhb6 F: 5’ GCAGCACGACTCTTCACAGGATG 3’
R: 5’ CGAACAGCAGCACAGACAGGAG 3’
Resp18 F: 5’ GTGAGGAGCAAGCAAGAGGAGAAAC 3’
R: 5’ GTGGTAAAGCATCGGTCCCTGTTC 3’
ENSMUST00000161503 F: 5’ TCTCAACTGCCCTTCACTCTCTTAG 3’
R: 5’ CATCCTTACCTCTCCTCCACTTCC 3’
AK041738 F: 5’ TCTCTTATCTGCCCCTGGACTTGTG 3’
R: 5’ CAAACCACTCACTCCCACTTCCTG 3’
ENSMUST00000210375 F: 5’ GCCACAACCAGCGCAGAC 3’
R: 5’ TCCTAGTCGTGCTTGAGAGTGAG 3’
ENSMUST00000126225 F: 5’ TGACCTCCCGGATGTTCTCT 3’
R:5’ TTTGTTTTGCAGGTCACGGC 3’
ENSMUST00000152254 F: 5’ GGTGGAAACATCAGGGCAGAAG 3’
R: 5’ GTTGGGAAGAAGGTCACTAGGAATG 3’
AF177020 F: 5’ AGTAGCCTAGTCTAGCACCCAACC 3’
R: 5’ CTGAGGGCGTGTTCTGACTTATAGC 3’

RT-qPCR primers of candidate mRNAs and lncRNAs.

2.7. Go and KEGG analysis

The roles differentially expressed mRNAs played in biological pathways were applied by GO and KEGG analysis. Metascape database was used to perform GO (biological process, cellular component, molecular function) and KEGG-enriched pathways (Zhou et al., 2019).

2.8. Establishment of lncRNA-mRNA co-expression network

Co-expression analysis was performed using the Sangerbox tools,1 which is a free online platform for data analysis. Pearson correlation coefficients (PCCs) were calculated to present the correlation between differentially expressed mRNAs and lncRNAs. LncRNA-mRNA pairs of co-expression were defined by |PCCs| > 0.99 and value of p <0.05. Cytoscape (version v3.9.1) was used for visualization of the co-expression network.

Table 2

Transcription ID Gene symbol Fold change (FC) value of p Regulation
NM_011315 Saa3 168.1430479 0.000232152 up
NM_009925 Col10a1 86.4408001 0.000152025 up
NM_009690 Cd5l 51.5293354 7.95593E-05 up
NM_008607 Mmp13 43.8708593 0.000143571 up
NM_008599 Cxcl9 31.6133097 0.000279797 up
NM_001123362 Prdm12 5.1310422 0.03269006 down
NM_011025 Oxt 2.6697481 0.042354851 down
NM_013833 Rax 2.5607242 0.008107894 down
NM_009049 Resp18 2.5185397 0.023800052 down
NM_153079 Nmur2 2.4707521 0.027802601 down

Top 5 up-regulated and down-regulated mRNAs post-ICH.

2.9. Construction of ce-RNA network

We selected DIANALncBaseV.32 to establish the lncRNA-miRNA interaction linkage. mRNA-miRNA interaction was predicted through miRDB-MicroRNA Target Prediction Database3 by setting the screening target score over 80. Then, the data of lncRNA-miRNA and mRNA-miRNA were integrated. Cytoscape (version v3.9.1) was used for visualization of the ce-RNA network. MCODE plugin of Cytoscape screened hub clusters according to filter criteria (degree cutoff: 0.2, node score cutoff: 0.2, k-core: 2, and max. Depth: 100).

2.10. Statistical analysis

GraphPad Prism 9 (GraphPad Software, Inc., La Jolla, CA, United States) was used for statistical analysis. Statistical differences were analyzed by t-test. Data from behavioral tests were analyzed by a two-way ANOVA test followed by Bonferroni’s multiple comparisons. The results were expressed as the mean ± SEM. A value of p of <0.05 was defined as statistically significant.

3. Results

3.1. Neurological deficits and histological alteration post-ICH

To evaluate the alteration of neurological function of mice after ICH, mNSS and foot-fault rates were recorded on days 1, 3, and 7 after the operation (n = 7–10) (Figure 2A). On the first day after ICH, the mNSS scores and foot-fault rates significantly increased compared with the sham group (Figures 2B,C). These neurological deficits persisted until day 7 when the tissues were harvested. Histological staining was conducted on brain tissues obtained on day 7 post-ICH. H&E staining revealed the hematoma area after stereotaxic injected autologous blood at the globus pallidus on day 7 (n = 5) (Figures 2D,E). In Nissl’s staining, injured neurons show up as dark neurons characterized by massive shrinkage and abnormal basophilia (Ooigawa et al., 2006). Therefore, we used the number of survival neurons to reflect the severity of neurological damage. The number of survival neurons decreased notably in contrast to dark neurons after ICH (n = 4) (Figures 2F,G).

Figure 2

Figure 2

ICH-induced neurological deficits and histology alterations. (A) Schedule of behavioral tests and tissue harvest. The mNSS (B) and foot-fault test (C) were measured for the two groups at three time points (day 1, day 3, day 7) post-ICH. (n = 7–10 mice in each group) Representative H&E staining (D) and Nissl staining (F) images on day 7 post-ICH. (E) Compared with the sham group, the ICH group exists significant hematoma. (n = 5 mice for each group) (G) The number of survival neurons in the ICH group was significantly decreased compared with that in the sham group. (n = 4 mice for each group) Data are presented as mean ± SEM, *p<0.05, **p<0.01, ****p<0.001.

3.2. Differential expression of mRNAs and lncRNAs in ICH

In order to investigate the alteration of mRNAs and lncRNAs in ICH, we used Arraystar Mouse lncRNA Microarray, including 22,692 mRNAs and 37,949 lncRNAs. Before microarray scanning, the quantity and quality of RNAs and labeling efficiency were validated by NanoDrop ND-1000 and specific activity (Supplementary Tables S1, S2). With an FC ≥ 2 and value of p <0.05 filtering, heat maps and volcano maps for differentially expressed mRNAs and lncRNAs were achieved (Figure 3). Five hundred and seventy mRNAs (including 543 up-regulated transcripts and 27 down-regulated transcripts) and 313 lncRNAs (including 273 up-regulated transcripts and 40 down-regulated transcripts) were identified between sham and ICH groups (n = 5). The heatmaps were also performed to visualize the cluster analysis of differentially expressed mRNAs and lncRNAs (Figures 3A,B). The top 5 up/down-regulated mRNAs and lncRNAs were listed in Tables 2, 3.

Figure 3

Figure 3

Differential expression of mRNAs and lncRNAs in sham and ICH group. Heat maps of differentially expressed mRNAs (A) and lncRNAs (B). Cluster analysis was carried out with rows. Volcano maps of differentially expressed mRNAs (C) and lncRNAs (D). Differentially expressed mRNAs and lncRNAs were according to the cutoff of log |FC| ≥ 2 and value of p <0.05. The up-regulated transcripts were marked in red, and the down-regulated transcripts in blue.

Table 3

GeneID Gene symbol Fold change (FC) p-value Regulation
ENSMUSG00000085786 Gm15987 13.6022041 0.000351216 up
AK143011 AK143011 12.7551382 6.7097E-06 up
ENSMUSG00000075010 AW112010 11.819704 0.000498888 up
ENSMUSG00000047415 Gpr68 10.1714873 3.46477E-05 up
ENSMUSG00000023349 Clec4n 10.0060345 0.000126743 up
AK030988 AK030988 6.0191623 0.002926642 down
ENSMUSG00000097620 4921514A10Rik 2.8131091 0.026730509 down
ENSMUSG00000087694 A530058N18Rik 2.5930208 0.044473003 down
AF177020 AF177020 2.4809815 0.008580349 down
ENSMUSG00000099696 2900052N01Rik 2.4651633 0.006742548 down

Top 5 up-regulated and down-regulated lncRNAs post-ICH.

Table 4

GO Category Description Count Log10(P)
GO:0008009 MF Chemokine activity 14 −12.03
GO:0038187 MF Pattern recognition receptor activity 11 −9.96
GO:0042287 MF MHC protein binding 14 −9.43
GO:0004857 MF Enzyme inhibitor activity 35 −9.12
GO:0140375 MF Immune receptor activity 19 −8.63
GO:0032103 BP Positive regulation of response to external stimulus 75 −37.57
GO:0006954 BP Inflammatory response 74 −33.6
GO:0002683 BP Negative regulation of immune system process 62 −26.76
GO:0050778 BP Positive regulation of immune response 75 −26.49
GO:0001775 BP Cell activation 69 −21.4
GO:0009897 CC External side of plasma membrane 69 −22.64
GO:0000323 CC Lytic vacuole 51 −14.43
GO:0045121 CC Membrane raft 34 −9.05
GO:0098797 CC Plasma membrane protein complex 41 −7.88
GO:0005581 CC Collagen trimer 12 −5.95

Top 5 GO terms enriched by differentially expressed mRNA post-ICH.

To verify the accuracy of microarray, 6 mRNAs (Saa3, Col10a1, Mmp13, Oxt, Pcdhb6 and Resp18) and 6 lncRNAs (ENSMUST00000161503, AK041738, ENSMUST00000210375, ENSMUST00000126225, ENSMUST00000152254, and AF177020) were selected for RT-qPCR. The 6 validated lncRNAs were confirmed a curtained conservation with human (Supplementary Table S3). Consistent with the results of array scanning, RT-qPCR confirmed that the expressions of Saa3, Col10a1, Mmp13, ENSMUST00000161503, AK041738, ENSMUST00000210375, and ENSMUST00000126225 were up-regulated post-ICH (n = 4–5) (Figures 4A–C,G). Meanwhile, the expressions of Oxt, Pcdhb6, Resp18, AK030988, ENSMUST00000152254, and AF177020 were lower in the ICH group than that of the sham group (n = 4–5) (Figures 4D–F,K). The RT-qPCR results supported the reliability of the microarray.

Figure 4

Figure 4

RT-qPCR confirmation of differentially expressed mRNAs and lncRNAs. The relative expression of Saa3 (A), Col10a1 (B), Mmp13 (C), ENSMUST00000161503 (G), AK041738 (H), ENSMUST00000210375 (I), and ENSMUST00000126225 (J) was upregulated in ICH group. The relative expression of Oxt (D), Pcdhb6 (E), Resp18 (F), ENSMUST00000152254 (K), and AF177020 (L) was downregulated in the ICH group. (n = 4–5 mice in each group) Data are presented as mean ± SEM, *p<0.05, **p<0.01.

3.3. Functional enrichment of differentially expressed mRNAs

Metascape database was used to analyze the functions of 570 differentially expressed mRNAs. GO analysis was performed in three dimensions involving biological process (BP), cellular compound (CC), and molecular function (MF). In total, 1,046 GO items were enriched, including 856 BP terms, 76 CC terms, and 114 MF terms. A subset of enriched terms was rendered as a network plot and colored by cluster-ID. Simultaneously, the plots (with a similarity >0.3) were connected by edges. Finally, we selected 20 clusters of best value of ps to be visualized by the network (Figure 5A). The bar graph showed the top 5 terms in BP, CC, and MF, respectively (Figure 5B) and listed in Table 4. BP was mainly involved in chemokine activity, pattern recognition, MHC protein binding, enzyme inhibitor activity, and immune receptor. CC was mainly focused on the positive regulation of response to external stimulus, inflammatory response, negative regulation of immune system process, positive regulation of immune response, and cell activation. Meanwhile, MF was enriched in the external side of plasma membrane, lytic vacuole, membrane raft, plasma membrane protein complex, and collagen trimer.

Figure 5

Figure 5

GO enrichment analysis of differentially expressed mRNA post-ICH. (A) Network of enriched GO terms. The nodes were colored by cluster-ID. If similarity >0.3, the nodes were connected by edges. (B) The bar graph of the top 5 GO terms in BP, CC, and MF, respectively.

KEGG pathway analysis was also conducted based on the Metascape platform. There were 78 pathways enriched with differentially expressed mRNAs. Ranked by lg value of p from the smallest to the largest, the top 20 pathways were shown in a bubble graph (Figure 6) and listed in Table 5, mainly containing cytokine-cytokine receptor interaction, hematopoietic cell lineage, TNF signaling pathway, NF-κB signaling pathway, apoptosis, natural killer cell mediated cytotoxicity, NOD-like receptor signaling pathway, lysosome, p53 signaling pathway, ferroptosis, and platelet activation.

Figure 6

Figure 6

Bubble graph representing KEGG pathways from the enrichment analysis of differentially expressed mRNA post-ICH. The top 20 enriched terms were shown ranked by value of p, and the size of bubbles represents the gene counts corresponding to the pathway.

Table 5

GO Description Count Log10(p) Log10(q)
mmu05332 Graft-versus-host disease 18 −13.75 −11.21
mmu04060 Cytokine-cytokine receptor interaction 34 −12.62 −10.78
mmu05133 Pertussis 18 −12.09 −10.41
mmu05323 Rheumatoid arthritis 18 −11.13 −9.62
mmu04640 Hematopoietic cell lineage 17 −9.57 −8.31
mmu05417 Lipid and atherosclerosis 24 −8.66 −7.43
mmu04668 TNF signaling pathway 17 −8.29 −7.12
mmu04380 Osteoclast differentiation 18 −8.26 −7.12
mmu04064 NF-kappa B signaling pathway 16 −7.92 −6.85
mmu05202 Transcriptional misregulation in cancer 23 −7.67 −6.63
mmu04210 Apoptosis 17 −7.06 −6.03
mmu04650 Natural killer cell mediated cytotoxicity 16 −7.03 −6.01
mmu05164 Influenza A 19 −6.9 −5.9
mmu05144 Malaria 11 −6.71 −5.73
mmu04621 NOD-like receptor signaling pathway 20 −6.12 −5.2
mmu04142 Lysosome 15 −5.64 −4.75
mmu04110 Cell cycle 14 −5.35 −4.48
mmu04115 p53 signaling pathway 10 −4.8 −3.94
mmu04216 Ferroptosis 7 −4.18 −3.34
mmu04611 Platelet activation 12 −4.04 −3.23

Top 20 KEGG terms enriched by differentially expressed mRNA post-ICH.

3.4. LncRNA-mRNA co-expression network

To investigate the relationship between differentially expressed lncRNAs and mRNAs, the Pearson correlation coefficients were calculated. If the absolute value of PCCs >0.99 and the value of p <0.05, we defined the lncRNA-mRNA pairs as co-expression. The co-expression network contained 57 nodes (21 lncRNAs and 36 mRNAs) and 38 connections (Figure 7). The top 10 co-expression lncRNA-mRNA pairs were ENSMUSG00000015451-Tnfaip8l2, ENSMUSG00000099930-Gm10134, ENSMUSG00000026656-Ctss, ENSMUSG00000020914-Top2a, KY468055-Parp9, ENSMUSG00000023992-Tmem273, ENSMUSG00000090272-Gmip, AK143011-Ms4a14, AK041738-Cybrd1 and AK041738-Csf2rb2 (Table 6).

Figure 7

Figure 7

The lncRNA-mRNA co-expression network post-ICH. The red circles represent mRNAs, and the green triangles represent lncRNAs.

Table 6

lncRNA mRNA p-value PCCs
ENSMUSG00000015451 Tnfaip8l2 1.92E-10 0.997425788
ENSMUSG00000099930 Gm10134 2.23E-10 0.997327088
ENSMUSG00000026656 Ctss 8.71E-10 0.996239743
ENSMUSG00000020914 Top2a 2.23E-09 0.99523903
KY468055 Parp9 4.14E-09 0.994444226
ENSMUSG00000023992 Tmem273 4.54E-09 0.994315768
ENSMUSG00000090272 Gmip 4.86E-09 0.994217239
AK143011 Ms4a14 5.96E-09 0.993914269
AK041738 Cybrd1 7.00E-09 0.993664268
AK041738 Csf2rb2 7.04E-09 0.993654819
ENSMUSG00000020914 Vim 9.16E-09 0.993221134
ENSMUSG00000002111 Slamf7 1.01E-08 0.99305264
ENSMUSG00000022901 Cd86 1.47E-08 0.992374707
ENSMUSG00000020914 Vim 1.47E-08 0.992362848
ENSMUSG00000015451 Ptpn18 1.62E-08 0.992180803
ENSMUSG00000054717 Cdkn2b 1.68E-08 0.992106452
ENSMUSG00000015451 Mafb 1.72E-08 0.992057591
ENSMUSG00000028581 Sat1 1.93E-08 0.991830154
ENSMUSG00000015451 Fcgr2b 1.94E-08 0.991814958
ENSMUSG00000047415 Slamf6 2.04E-08 0.99171702
ENSMUSG00000026656 C1qc 2.27E-08 0.991492486
ENSMUSG00000050708 Tmc4 2.29E-08 0.991472954
ENSMUSG00000026304 Gfap 2.36E-08 0.991409206
ENSMUSG00000020914 Abcg3 2.52E-08 0.991267119
ENSMUSG00000022901 Fcgr4 2.54E-08 0.991245292
ENSMUSG00000048852 Gm12185 2.68E-08 0.991130117
ENSMUSG00000020914 C4b 2.90E-08 0.990951017
AK033727 Tlr2 2.94E-08 0.990920626
AK079793 Myo1g 2.94E-08 0.990920012
ENSMUSG00000015451 Lyl1 3.01E-08 0.9908701
ENSMUSG00000109251 Siglec1 3.35E-08 0.990620722
ENSMUSG00000086150 Tnfsf11 3.50E-08 0.990514182
AK041738 Tnfsf11 3.57E-08 0.990467798
AK041738 Clec4n 3.84E-08 0.990295043
ENSMUSG00000020914 Ftl1 3.93E-08 0.990235332
ENSMUSG00000015451 Plcg2 4.11E-08 0.990127434
AK143011 Tnfaip2 4.20E-08 0.990072078
ENSMUSG00000020914 Arhgdib 4.25E-08 0.990042744

LncRNA-mRNA co-expression pairs. (|PCCs| > 0.99).

3.5. LncRNA-miRNA-mRNA (ce-RNA) network

For the purpose of exploring further functions based on the mechanism of competitive endogenous RNA (ce-RNA), we predicted the miRNAs which were correlated with differentially expressed lncRNAs and mRNAs. We defined lncRNA-miRNA pairs and mRNA-miRNA pairs with target scores >80 as closely correlated. A total of 570 differentially expressed mRNAs were correlated with 1,514 miRNAs. Additionally, 242 miRNAs were predicted by 313 differentially expressed lncRNAs based on the DIANALncBaseV.3 database. Totally, we defined 490 mRNA-miRNA pairs and 416 lncRNA-miRNA pairs. We combined with miRNAs predicted by mRNAs and lncRNAs to establish the ce-RNA network (Figure 8A). The ce-RNA network was including 303 nodes (29 lncRNAs, 163 mRNAs, and 111 miRNAs) and 906 edges. Then, we used MCODE, a plugin module of Cytoscape, to analyze the ce-RNA network. Three clusters were screened based on the filter criteria (degree cutoff: 0.2, node score cutoff: 0.2, k-core: 2, and max. Depth: 100) (Figures 8B–D). Cluster1 was composed with 2 mRNAs (Col1a2 and Tgfbr1) and 5 mircoRNAs (mmu-let-7f-5p, mmu-let-7c-5p, mmu-let-7i-5p, mmu-let-7 g-5p and mmu-miR-98-5p) (Figure 8B). In Cluster2 (Has2, Bcl2l11, ENSMUSG00000099757, ENSMUSG00000075010, ENSMUSG00000109644, mmu-miR-106b-5p, mmu-miR-30a-5p, mmu-let-7a-5p, mmu-let-7e-5p) (Figure 8C), there were 7 mRNA-miRNA-lncRNA pairs be screened comparing with 4 pairs in Cluster3 (Kif23, ENSMUSG00000093594, ENSMUSG00000091575, mmu-miR-106a-5p, mmu-miR-93-5p, mmu-miR-20b-5p, mmu-miR-15b-5p, mmu-miR-18a-5p) (Figure 8D).

Figure 8

Figure 8

The lncRNA-miRNA-mRNA interactions post-ICH. (A) The ce-RNA network post-ICH. The hub clusters of the ce-RNA network: are Cluster1 (B), Cluster2 (C), and Cluster3 (D). The yellow circles represent mRNAs, green rhombuses represent miRNAs, and red hexagons represent lncRNAs.

4. Discussion

In our previous study, we investigated the relationship between mRNA and lncRNA in the ICH model 21 days (chronic phase) post-ICH (Hanjin et al., 2018). In this study, to understand the pathophysiological changes during acute and subacute phases of intracerebral hemorrhage, we used microarray to obtain the differential expression of mRNAs and lncRNAs between sham and ICH groups. Afterward, we validated the accuracy of results from the microarray relying on RT-qPCR. The relative expression of 3 mRNAs and 4 lncRNAs (Saa3, Col10a1, Mmp13, ENSMUST00000161503, AK041738, ENSMUST00000210375, and ENSMUST00000126225) were up-regulated in contrast to the other 3 mRNAs and 2 lncRNAs (Oxt, Pcdhb6, Resp18, ENSMUST00000152254 and AF177020). These were consistent with the results of microarray, demonstrating that the results of microarray were reliable. To investigate the hiding function of differentially expressed mRNAs and lncRNAs, GO/KEGG was performed. The pathways and functions enriched by GO/KEGG mainly focused on the pathophysiological processes related to immune response, inflammatory factors, inflammatory cell activation, and apoptosis. In addition, the results of GO/KEGG were also involved in ferroptosis, platelet activation, and biological processes of membrane. The regulations of mRNA transcription and translation by lncRNA have attracted increasing attention. Besides, we are also concerned about the ternary regulatory relationship of mRNA-miRNA-lncRNA. Therefore, to gain further insight into the association between differentially expressed mRNAs and lncRNAs, we constructed the lncRNA-mRNA co-expression network and the ce-RNA network. Thirty-eight lncRNA-mRNA pairs were focused on in the co-expression network. Three hub clusters were selected from the complicated ce-RNA network.

In our study, the top five up-regulated mRNAs were Saa3, Col10a1, Cd5l, Mmp13, and Cxcl19. To date, none of these five mRNAs have been investigated in ICH. Saa3 is the target gene of C/EBPβ, which is up-regulated in cerebral ischemia (Hanjin et al., 2018). However, the exact function of Saa3 in stroke is still unclear. There is a study speculates that Saa3 plays a role in the recruitment and adhesion of monocytes (Han et al., 2007). Under the hypoxic condition, Col10a1 was up-regulated in certain types of cells, including human retinal microvascular endothelial cells, umbilical vein endothelial cells, RPE19 cells, and RF/6A cells (Lv et al., 2022). It indicates that Col10a1 may be related to microvascular formation and angiogenesis. Danger-associated molecular pattern (DAMP) is a new phase secondary to the primary injury, such as cell death. The sterile inflammation occurs following DAMP. It was found that CD5 antigen-like (Cd5l) administration attenuates sterile inflammation and that Cd5l deficiency courses severe neurological damage with higher mortality (Maehara et al., 2021). The up-regulation of Cd5l and Cxcl19 in the study may be a compensatory elevation after ICH. As a family of diverse zinc endopeptidases, matrix metalloproteinases (Mmps) play a major role in the physiology and pathophysiology of the central nervous system (CNS) (Rosell and Lo, 2008). The baseline of plasma Mmp13 has been proven to be strongly associated with the severity of the ischemic stroke and indicates poorer outcomes (Ma et al., 2016). Similar to the top 5 up-regulated mRNAs, most of the top 5 down-regulated mRNAs are still uninvestigated in ICH. Prdm12 was confirmed to contribute to nervous system development (Kinameri et al., 2008). Oxytocin (OXT) is the ligand of the Oxytocin receptor (OXTR), which is located in microvascular endothelial cells. Exogenous OXT administration benefits cerebral ischemia by reducing ROS and neuron apoptosis (Stary et al., 2019). Rax also plays a significant role in forebrain development and neuronal regeneration (Muranishi et al., 2012; Stary et al., 2019). Moreover, Resp18 knockdown mice have poorer locomotor activity in Parkinson’s disease (Su et al., 2018). It may prompt Resp18 to be involved in CNS function. Nmur2 (neuromedin U receptor 2), the receptor of neuromedin U (NMU), is identified predominantly centrally. The role of receptor-ligand binding mainly involves central energy homeostasis, thermoregulation, circadian rhythms, and stress response (Malendowicz and Rucinski, 2021). In brief, the top five up-regulated mRNAs are mainly associated with inflammation and angiogenesis. While most of the top 5 down-regulated mRNAs are related to CNS development, neuroprotection, and neurodegeneration. The specific roles of these genes in the pathophysiological process of ICH remain to be further investigated.

The essential biological function and pathways enriched by GO/KEGG are immune responses and inflammation. The CNS injury after ICH includes primary injury and secondary injury. The extravasated blood from ruptured vessels causes hematoma formation within a few hours (Tschoe et al., 2020). The resident glial cells are suggested to be the earliest inflammatory cells in response to the primary injury. Additionally, neuroinflammation is a significant progress of secondary injury. Neuroinflammation is triggered by thrombin, Hemoglobin (Hb), Hemin, iron, and other stimulator released by blood and dead cells (Zhou et al., 2014). Besides, brain edema aggravates damage to neurons and glial cells, causing complex inflammatory cascades. Therefore, inflammation and immune response extend almost through the entire process of ICH. Notably, biological processes of GO analysis enriched for two opposing signaling pathways simultaneously, negative regulation of immune system progress and positive regulation of immune response. Activation of immune response is a double-edged sword. Taking microglia activation as an example, it has both positive and negative effects on ICH outcomes. In the hyperacute and acute phase of ICH, microglia activation produces a mass of pro-inflammatory factors leading to an inflammatory waterfall reaction, eventually promotes blood–brain barrier (BBB) injury, edema, and neurologic deficit. However, M2, a phenotype of microglia, improves CNS recovery by phagocytosis and secreting anti-inflammatory factors like IL-10 and CD36 (Liu et al., 2021). To sum up, the organism exerts pro-inflammatory effects in parallel with anti-inflammation in an attempt to achieve some balance after ICH.

GO/KEGG analysis indicates the functions of differentially expressed mRNAs associated with apoptosis and ferroptosis. It has been described previously that cell death improves DAMP and drives inflammatory cascades. Various types of cell death exert during ICH, including apoptosis, autophagocytosis, necrosis, necroptosis, and pyroptosis. Apoptosis is the programmed death of cells. Complex factors post-stroke lead to apoptosis, which is concerned with ATP depletion, misfolded proteins, calcium influx, free radicals, inflammatory response, α-Syn aggregation, and so on (Sekerdag et al., 2018). Red blood cells released from the ruptured vessels contain an abundance of Hb. Microglia cells phagocytize Hb, which is then metabolized into ferrous iron by ROS and lipid peroxidation. The pathological products trigger the onset of ferroptosis jointly (Wan et al., 2019).

Some typical pathways of stroke were also enriched by KEGG analysis. TNF signaling pathway has a compact relationship with microglia activation and polarization (Zhou et al., 2014). The NF-κB signaling pathway is activated within a few minutes after ICH and persists for more than a week. Accumulating evidence demonstrates that NF-κB activation leads to a large amount of pro-inflammatory cytokines release, especially TNF-α and IL-1β (Sivandzade et al., 2019). NOD-like thermal receptor protein domain associated protein 3 (NLRP3), belonging to the NOD-like receptor family, is confirmed to be a vital contributor to ICH. NLRP3 increases in microglia cells remarkably in ICH and down-regulation of NLRP3 attenuates inflammation after ICH (Yang et al., 2015). The signaling pathway of p53 contributes to apoptotic programs closely after stroke. After the stroke onset, p53 is activated immediately. Then, the activated p53 translocates to the outer membrane of mitochondria and into nuclei. On the outer membrane of mitochondria, p53 binding with B-cell lymphoma-2 (BCL-2) leads to cytosol release for the apoptotic program. In nuclei, p53 interacts with its responsive elements to initiate the transcription of apoptosis (Almeida et al., 2021). The molecules upstream or downstream of these signaling pathways may be the potential targets for the treatment of ICH.

A growing number of researchers show solicitude for the effects of non-coding RNA in pro-ICH and post-ICH. It is reported that lncRNA MEG3 participates in reducing oxidative stress and inflammation by targeting miR-181b in ICH (Xie B. et al., 2021). Apart from elevating the inflammatory level, lncRNA Blnc1 is able to mediate the permeability through the PPAR-γ/SIRT6/FoxO3 pathway in ICH (Xie L. et al., 2021). In our ce-RNA network, three significant clusters were identified. Some interactions of them have been proven to be involved in the pathological processes of a multiplicity of diseases. For instance, let-7f-5p regulated osteogenesis and bone formation by targeting TGFBR1 to ameliorate osteoporosis (Shen et al., 2019). Moreover, TGFBR1 also contributed to pulmonary fibrosis and non-small cell lung cancer (NSLC) by interacting with let-7i-5p and miR-98-5p, respectively (Jiang et al., 2019; Xu C. et al., 2022). The miR-98-5p-TGFBR1 interaction could be regulated by lncRNA Linc00511 which participated in the recurrence of NSLC (Li et al., 2022). Zhong and colleagues confirmed that miR-98-5p promoted COL1A2 inducing extracellular matrix deposition and contributing to asthma pathogenesis in human (Zhong et al., 2022). Additionally, miR-30a-5p was detected to down-regulate BCL2L11 depending on P65 for aggravating multiple myeloma (Xie et al., 2022). In Cluster 2, miR-106b-5p is a hub node. Chen and colleagues found that lncRNA H19 inhibits ICH injuries in miR-106b-5p dependent manner (Chen et al., 2021). However, the interactions among miR-106b-5p, Has2, Bcl2l11, ENSMUSG00000099757 and ENSMUSG00000075010 need further detected. While H19 is regarded as an accomplice of brain damage in another study (Mao et al., 2022). Similarly, MALAT1 inhibits inflammation and oxidative stress in ICH, as in ischemic stroke (Xu Z. et al., 2022). Through multiple post-transcriptional mechanisms, lncRNA can act as miRNA regulators or effectors, ultimately regulating the expression of mRNA in ICH. Nonetheless, the journey of finding out how lncRNA-miRNA-mRNA is involved in the physiological and pathological mechanisms of ICH is still far from over.

Compared to other studies focusing on the function of lncRNA in ICH (Cao et al., 2020; Chen et al., 2020; Xie B. et al., 2021), we employed the autologous blood injection model rather than the collagenase-induced model. The autologous blood injection model and collagenase-induced model are the most frequent and stable ICH models (Ma et al., 2011). Although neither of the two ICH models can perfectly mimic the pathophysiological situation of the disease, they have their unique advantages in simulating some processes of ICH. In the autologous blood injection model, blood enters the brain parenchyma from a single needle in a short time, which mimics the initial hemorrhage occurrence. Furthermore, the volume of the hematoma can be controlled (relating to the injected blood volume) (Kirkman et al., 2011). To investigate edema and inflammation, which is the most significant process in the acute and subacute phases, we chose the autologous blood injection model to conduct this research.

The study has several limitations that need to be addressed. Differentially expressed mRNAs and lncRNAs did not undergo functional experiments to determine their exact roles in ICH. Furthermore, the expression of differentially expressed RNAs at other phases was not investigated.

In summary, our study identified hundreds of differentially expressed mRNAs and lncRNAs in the acute ICH mouse model. Top RNA molecules of them may be the biomarker of ICH. Moreover, we analyzed the underlying relationships of lncRNA-mRNA and lncRNA-miRNA-mRNA, which may offer a potential target for ICH treatment.

Funding

This study was funded by grants from the National Natural Science Foundation of China (no’s. 81874425, 81973665, and 82174259); the Hunan Provincial Key Research and Development Program (2022SK2015); the Excellent Clinical Talent Training Project of SATCM (2022-543313); the Leadership Training Project of HNATCM (2022–24); the Shennong Youth Scholar Project of Hunan (2022). The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author. The microarray data presented in the study are publicly available in the NCBI GEO database. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi accession number GSE227683.

Ethics statement

The animal study was reviewed and approved by Institutional Animal Care and Use Committee of CSU.

Author contributions

ZY: investigation, validation, visualization, and writing—original draft. EH: investigation, validation, writing—review and editing. YC: investigation and validation. WZ: investigation. QC: visualization and software. TL: investigation. ZL: investigation. YW: supervision and writing—review and editing. TT: resources, data curation, funding acquisition, supervision, and writing—review and editing. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2023.1166875/full#supplementary-material

References

  • 1

    Almeida A. Sánchez-Morán I. Rodríguez C. (2021). Mitochondrial-nuclear p53 trafficking controls neuronal susceptibility in stroke. IUBMB Life73, 582–591. doi: 10.1002/iub.2453

  • 2

    Bridges M. C. Daulagala A. C. Kourtidis A. (2021). LNCcation: lncRNA localization and function. J. Cell Biol.220:e202009045. doi: 10.1083/jcb.202009045

  • 3

    Canseco-Rodriguez A. Masola V. Aliperti V. Meseguer-Beltran M. Donizetti A. Sanchez-Perez A. M. (2022). Long non-coding RNAs, extracellular vesicles and inflammation in Alzheimer's disease. Int. J. Mol. Sci.23:13171. doi: 10.3390/ijms232113171

  • 4

    Cao F. Guo Y. Zhang Q. Fan Y. Liu Q. Song J. et al . (2020). Integration of transcriptome resequencing and quantitative proteomics analyses of collagenase VII-induced intracerebral hemorrhage in mice. Front. Genet.11:551065. doi: 10.3389/fgene.2020.551065

  • 5

    Chen B. Wang H. Lv C. Mao C. Cui Y. (2021). Long non-coding RNA H19 protects against intracerebral hemorrhage injuries via regulating microRNA-106b-5p/acyl-CoA synthetase long chain family member 4 axis. Bioengineered12, 4004–4015. doi: 10.1080/21655979.2021.1951070

  • 6

    Chen J. X. Wang Y. P. Zhang X. Li G. X. Zheng K. Duan C. Z. (2020). lncRNA Mtss1 promotes inflammatory responses and secondary brain injury after intracerebral hemorrhage by targeting miR-709 in mice. Brain Res. Bull.162, 20–29. doi: 10.1016/j.brainresbull.2020.04.017

  • 7

    Choi S. W. Kim H. W. Nam J. W. (2019). The small peptide world in long noncoding RNAs. Brief. Bioinform.20, 1853–1864. doi: 10.1093/bib/bby055

  • 8

    Deng Y. Chen D. Wang L. Gao F. Jin B. Lv H. et al . (2019). Silencing of Long noncoding RNA Nespas aggravates microglial cell death and Neuroinflammation in ischemic stroke. Stroke50, 1850–1858. doi: 10.1161/strokeaha.118.023376

  • 9

    Fei X. Dou Y. N. Wang L. Wu X. Huan Y. Wu S. et al . (2022). Homer1 promotes the conversion of A1 astrocytes to A2 astrocytes and improves the recovery of transgenic mice after intracerebral hemorrhage. J. Neuroinflammation19:67. doi: 10.1186/s12974-022-02428-8

  • 10

    GBD 2019 Stroke Collaborators 2021 Stroke Collaborators (2021). Global, regional, and national burden of stroke and its risk factors 1990-2019: a systematic analysis for the global burden of disease study 2019. Lancet Neurol.20, 795–820. doi: 10.1016/s1474-4422(21)00252-0

  • 11

    Greenberg S. M. Ziai W. C. Cordonnier C. Dowlatshahi D. Francis B. Goldstein J. N. et al . (2022). 2022 guideline for the Management of Patients with Spontaneous Intracerebral Hemorrhage: a guideline from the American Heart Association/American Stroke Association. Stroke53, e282–e361. doi: 10.1161/str.0000000000000407

  • 12

    Han C. Y. Subramanian S. Chan C. K. Omer M. Chiba T. Wight T. N. et al . (2007). Adipocyte-derived serum amyloid A3 and hyaluronan play a role in monocyte recruitment and adhesion. Diabetes56, 2260–2273. doi: 10.2337/db07-0218

  • 13

    Hanjin C. Tao L. Pengfei L. Ali Y. Huajun Z. Jiekun L. et al . (2018). Altered Long noncoding RNA and messenger RNA expression in experimental intracerebral hemorrhage -a preliminary study. Cell. Physiol. Biochem.45, 1284–1301. doi: 10.1159/000487464

  • 14

    Jarroux J. Morillon A. Pinskaya M. (2017). History, discovery, and classification of lncRNAs. Adv. Exp. Med. Biol.1008, 1–46. doi: 10.1007/978-981-10-5203-3_1

  • 15

    Jiang F. Yu Q. Chu Y. Zhu X. Lu W. Liu Q. et al . (2019). MicroRNA-98-5p inhibits proliferation and metastasis in non-small cell lung cancer by targeting TGFBR1. Int. J. Oncol.54, 128–138. doi: 10.3892/ijo.2018.4610

  • 16

    Kinameri E. Inoue T. Aruga J. Imayoshi I. Kageyama R. Shimogori T. et al . (2008). Prdm proto-oncogene transcription factor family expression and interaction with the notch-Hes pathway in mouse neurogenesis. PLoS One3:e3859. doi: 10.1371/journal.pone.0003859

  • 17

    Kirkman M. A. Allan S. M. Parry-Jones A. R. (2011). Experimental intracerebral hemorrhage: avoiding pitfalls in translational research. J. Cereb. Blood Flow Metab.31, 2135–2151. doi: 10.1038/jcbfm.2011.124

  • 18

    Kopp F. Mendell J. T. (2018). Functional classification and experimental dissection of Long noncoding RNAs. Cells172, 393–407. doi: 10.1016/j.cell.2018.01.011

  • 19

    Li C. Li Z. Yi H. Liu Z. (2022). IncRNA Linc00511 upregulation elevates TGFBR1 and participates in the postoperative distant recurrence of non-small-cell lung cancer by targeting miR-98-5p. Crit. Rev. Eukaryot. Gene Expr.32, 1–10. doi: 10.1615/CritRevEukaryotGeneExpr.2022039498

  • 20

    Li T. Luo Y. Zhang P. Guo S. Sun H. Yan D. et al . (2020). LncRNA MEG3 regulates microglial polarization through KLF4 to affect cerebral ischemia-reperfusion injury. J. Appl. Physiol. (1985)129, 1460–1467. doi: 10.1152/japplphysiol.00433.2020

  • 21

    Li P. Tang T. Liu T. Zhou J. Cui H. He Z. et al . (2019). Systematic analysis of tRNA-derived small RNAs reveals novel potential therapeutic targets of traditional Chinese medicine (Buyang-Huanwu-decoction) on intracerebral hemorrhage. Int. J. Biol. Sci.15, 895–908. doi: 10.7150/ijbs.29744

  • 22

    Liu J. Liu L. Wang X. Jiang R. Bai Q. Wang G. (2021). Microglia: a double-edged sword in intracerebral hemorrhage from basic mechanisms to clinical research. Front. Immunol.12:675660. doi: 10.3389/fimmu.2021.675660

  • 23

    Long F. Q. Su Q. J. Zhou J. X. Wang D. S. Li P. X. Zeng C. S. et al . (2018). LncRNA SNHG12 ameliorates brain microvascular endothelial cell injury by targeting miR-199a. Neural Regen. Res.13, 1919–1926. doi: 10.4103/1673-5374.238717

  • 24

    Lv D. Chen D. Wang Z. Cui Z. Ma J. H. Ji S. et al . (2022). COL10A1 is a novel factor in the development of choroidal neovascularization. Microvasc. Res.139:104239. doi: 10.1016/j.mvr.2021.104239

  • 25

    Ma Q. Khatibi N. H. Chen H. Tang J. Zhang J. H. (2011). History of preclinical models of intracerebral hemorrhage. Acta Neurochir. Suppl.111, 3–8. doi: 10.1007/978-3-7091-0693-8_1

  • 26

    Ma F. Rodriguez S. Buxo X. Morancho A. Riba-Llena I. Carrera A. et al . (2016). Plasma matrix metalloproteinases in patients with stroke during intensive rehabilitation therapy. Arch. Phys. Med. Rehabil.97, 1832–1840. doi: 10.1016/j.apmr.2016.06.007

  • 27

    Maehara N. Taniguchi K. Okuno A. Ando H. Hirota A. Li Z. et al . (2021). AIM/CD5L attenuates DAMPs in the injured brain and thereby ameliorates ischemic stroke. Cell Rep.36:109693. doi: 10.1016/j.celrep.2021.109693

  • 28

    Malendowicz L. K. Rucinski M. (2021). Neuromedins NMU and NMS: an updated overview of their functions. Front. Endocrinol. (Lausanne)12:713961. doi: 10.3389/fendo.2021.713961

  • 29

    Mao S. Huang H. Chen X. (2022). lncRNA H19 aggravates brain injury in rats following experimental intracerebral hemorrhage via NF-κB pathway. Comput. Math. Methods Med.2022, 3017312–3017317. doi: 10.1155/2022/3017312

  • 30

    Mattick J. S. Amaral P. P. Carninci P. Carpenter S. Chang H. Y. Chen L. L. et al . (2023). Long non-coding RNAs: definitions, functions, challenges and recommendations. Nat. Rev. Mol. Cell Biol. doi: 10.1038/s41580-022-00566-8

  • 31

    Muranishi Y. Terada K. Furukawa T. (2012). An essential role for Rax in retina and neuroendocrine system development. Develop. Growth Differ.54, 341–348. doi: 10.1111/j.1440-169X.2012.01337.x

  • 32

    Na C. Wen-Wen C. Li W. Ao-Jia Z. Ting W. (2022). Significant role of Long non-coding RNAs in Parkinson's disease. Curr. Pharm. Des.28, 3085–3094. doi: 10.2174/1381612828666220922110551

  • 33

    Ni X. Su Q. Xia W. Zhang Y. Jia K. Su Z. et al . (2020). Knockdown lncRNA NEAT1 regulates the activation of microglia and reduces AKT signaling and neuronal apoptosis after cerebral ischemic reperfusion. Sci. Rep.10:19658. doi: 10.1038/s41598-020-71411-1

  • 34

    Ooigawa H. Nawashiro H. Fukui S. Otani N. Osumi A. Toyooka T. et al . (2006). The fate of Nissl-stained dark neurons following traumatic brain injury in rats: difference between neocortex and hippocampus regarding survival rate. Acta Neuropathol.112, 471–481. doi: 10.1007/s00401-006-0108-2

  • 35

    Organization, G. W. H . (2022). World health statistics 2022: monitoring health for the SDGs, sustainable development goals [online] Available at: https://www.who.int/news/item/20-05-2022-world-health-statistics-2022 (Accessed 20 May 2022).

  • 36

    Ransohoff J. D. Wei Y. Khavari P. A. (2018). The functions and unique features of long intergenic non-coding RNA. Nat. Rev. Mol. Cell Biol.19, 143–157. doi: 10.1038/nrm.2017.104

  • 37

    Rosell A. Lo E. H. (2008). Multiphasic roles for matrix metalloproteinases after stroke. Curr. Opin. Pharmacol.8, 82–89. doi: 10.1016/j.coph.2007.12.001

  • 38

    Rynkowski M. A. Kim G. H. Komotar R. J. Otten M. L. Ducruet A. F. Zacharia B. E. et al . (2008). A mouse model of intracerebral hemorrhage using autologous blood infusion. Nat. Protoc.3, 122–128. doi: 10.1038/nprot.2007.513

  • 39

    Sekerdag E. Solaroglu I. Gursoy-Ozdemir Y. (2018). Cell death mechanisms in stroke and novel molecular and cellular treatment options. Curr. Neuropharmacol.16, 1396–1415. doi: 10.2174/1570159x16666180302115544

  • 40

    Shen G. Y. Ren H. Shang Q. Zhao W. H. Zhang Z. D. Yu X. et al . (2019). Let-7f-5p regulates TGFBR1 in glucocorticoid-inhibited osteoblast differentiation and ameliorates glucocorticoid-induced bone loss. Int. J. Biol. Sci.15, 2182–2197. doi: 10.7150/ijbs.33490

  • 41

    Sivandzade F. Prasad S. Bhalerao A. Cucullo L. (2019). NRF2 and NF-Ò›B interplay in cerebrovascular and neurodegenerative disorders: molecular mechanisms and possible therapeutic approaches. Redox Biol.21:101059. doi: 10.1016/j.redox.2018.11.017

  • 42

    Stary C. M. Xu L. Voloboueva L. A. Alcántara-Hernández M. Arvola O. J. Idoyaga J. et al . (2019). Nursing markedly protects postpartum mice from stroke: associated central and peripheral Neuroimmune changes and a role for oxytocin. Front. Neurosci.13:609. doi: 10.3389/fnins.2019.00609

  • 43

    Su J. Wang H. Yang Y. Wang J. Li H. Huang D. et al . (2018). RESP18 deficiency has protective effects in dopaminergic neurons in an MPTP mouse model of Parkinson's disease. Neurochem. Int.118, 195–204. doi: 10.1016/j.neuint.2018.06.010

  • 44

    Sun X. L. Wang Z. L. Wu Q. Jin S. Q. Yao J. Cheng H. (2019). LncRNA RMST activates TAK1-mediated NF-κB signaling and promotes activation of microglial cells via competitively binding with hnRNPK. IUBMB Life71, 1785–1793. doi: 10.1002/iub.2125

  • 45

    Tan X. Liu Y. Liu Y. Zhang T. Cong S. (2021). Dysregulation of long non-coding RNAs and their mechanisms in Huntington's disease. J. Neurosci. Res.99, 2074–2090. doi: 10.1002/jnr.24825

  • 46

    Tschoe C. Bushnell C. D. Duncan P. W. Alexander-Miller M. A. Wolfe S. Q. (2020). Neuroinflammation after intracerebral hemorrhage and potential therapeutic targets. J Stroke22, 29–46. doi: 10.5853/jos.2019.02236

  • 47

    Wan J. Ren H. Wang J. (2019). Iron toxicity, lipid peroxidation and ferroptosis after intracerebral haemorrhage. Stroke Vasc. Neurol.4, 93–95. doi: 10.1136/svn-2018-000205

  • 48

    Wang J. Zhao H. Fan Z. Li G. Ma Q. Tao Z. et al . (2017). Long noncoding RNA H19 promotes Neuroinflammation in ischemic stroke by driving histone deacetylase 1-dependent M1 microglial polarization. Stroke48, 2211–2221. doi: 10.1161/strokeaha.117.017387

  • 49

    Wang H. Zheng X. Jin J. Zheng L. Guan T. Huo Y. et al . (2020). LncRNA MALAT1 silencing protects against cerebral ischemia-reperfusion injury through miR-145 to regulate AQP4. J. Biomed. Sci.27:40. doi: 10.1186/s12929-020-00635-0

  • 50

    Wen Y. Yu Y. Fu X. (2017). LncRNA Gm4419 contributes to OGD/R injury of cerebral microglial cells via IκB phosphorylation and NF-κB activation. Biochem. Biophys. Res. Commun.487, 923–929. doi: 10.1016/j.bbrc.2017.05.005

  • 51

    Xie B. Qiao M. Xuan J. (2021). lncRNA MEG3 downregulation relieves intracerebral hemorrhage by inhibiting oxidative stress and inflammation in an miR-181b-dependent manner. Med. Sci. Monit.27:e929435. doi: 10.12659/msm.929435

  • 52

    Xie L. Wang Y. Chen Z. (2021). LncRNA Blnc1 mediates the permeability and inflammatory response of cerebral hemorrhage by regulating the PPAR-γ/SIRT6/FoxO3 pathway. Life Sci.267:118942. doi: 10.1016/j.lfs.2020.118942

  • 53

    Xie L. Wei J. Gao Z. Huang H. Ju S. Wang X. (2022). Significance of a tumor microenvironment-mediated P65-miR-30a-5p-BCL2L11 amplification loop in multiple myeloma. Exp. Cell Res.415:113113. doi: 10.1016/j.yexcr.2022.113113

  • 54

    Xu C. Hou L. Zhao J. Wang Y. Jiang F. Jiang Q. et al . (2022). Exosomal let-7i-5p from three-dimensional cultured human umbilical cord mesenchymal stem cells inhibits fibroblast activation in silicosis through targeting TGFBR1. Ecotoxicol. Environ. Saf.233:113302. doi: 10.1016/j.ecoenv.2022.113302

  • 55

    Xu Z. Zhao B. Mao J. Sun Z. (2022). Knockdown of long noncoding RNA metastasis-associated lung adenocarcinoma transcript 1 protects against intracerebral hemorrhage through microRNA-146a-mediated inhibition of inflammation and oxidative stress. Bioengineered13, 3969–3980. doi: 10.1080/21655979.2022.2031401

  • 56

    Yang Z. Zhong L. Xian R. Yuan B. (2015). MicroRNA-223 regulates inflammation and brain injury via feedback to NLRP3 inflammasome after intracerebral hemorrhage. Mol. Immunol.65, 267–276. doi: 10.1016/j.molimm.2014.12.018

  • 57

    Zhang X. Tang X. Liu K. Hamblin M. H. Yin K. J. (2017). Long noncoding RNA Malat1 regulates cerebrovascular pathologies in ischemic stroke. J. Neurosci.37, 1797–1806. doi: 10.1523/jneurosci.3389-16.2017

  • 58

    Zhao W. Wu C. Stone C. Ding Y. Ji X. (2020). Treatment of intracerebral hemorrhage: current approaches and future directions. J. Neurol. Sci.416:117020. doi: 10.1016/j.jns.2020.117020

  • 59

    Zhong Z. G. Dong C. P. Guo X. H. Chen J. Zhu L. P. Zhang M. (2022). Long noncoding RNA ANRIL up-regulates CCND1 via sponging miR-98-5p to promote TGF-β1-induced human airway smooth muscle cell proliferation, migration, and extracellular matrix deposition. Kaohsiung J. Med. Sci.38, 633–642. doi: 10.1002/kjm2.12538

  • 60

    Zhou Y. Wang Y. Wang J. Anne Stetler R. Yang Q. W. (2014). Inflammation in intracerebral hemorrhage: from mechanisms to clinical translation. Prog. Neurobiol.115, 25–44. doi: 10.1016/j.pneurobio.2013.11.003

  • 61

    Zhou Y. Zhou B. Pache L. Chang M. Khodabakhshi A. H. Tanaseichuk O. et al . (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun.10:1523. doi: 10.1038/s41467-019-09234-6

Summary

Keywords

mRNA, lncRNA, co-expression network, competitive endogenous RNA, intracerebral hemorrhage

Citation

Yu Z, Hu E, Cai Y, Zhu W, Chen Q, Li T, Li Z, Wang Y and Tang T (2023) mRNA and lncRNA co-expression network in mice of acute intracerebral hemorrhage. Front. Mol. Neurosci. 16:1166875. doi: 10.3389/fnmol.2023.1166875

Received

23 February 2023

Accepted

05 April 2023

Published

28 April 2023

Volume

16 - 2023

Edited by

Xiao Xiao, Fudan University, China

Reviewed by

Zhi Zheng, Memorial Sloan Kettering Cancer Center, United States; Lingjie Sang, University of Texas Southwestern Medical Center, United States; Hao Wu, Exelixis, United States

Updates

Copyright

*Correspondence: Tao Tang,

†These authors have contributed equally to this work and share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics