Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2018.9000 mmr-18-01-0333 Articles Bioinformatics analysis of microarray data to reveal the pathogenesis of brain ischemia HeJiaxuan 1 GaoYa 2 WuGang 1 LeiXiaoming 1 ZhangYong 1 PanWeikang 2 YuHui 2 Department of Anesthesia, Second Affiliated Hospital of Xi'an Jiaotong University, Xi'an, Shaanxi 710004, P.R. China Department of Pediatric Surgery, Second Affiliated Hospital of Xi'an Jiaotong University, Xi'an, Shaanxi 710004, P.R. China Correspondence to: Dr Ya Gao, Department of Pediatric Surgery, Second Affiliated Hospital of Xi'an Jiaotong University, 157 XiWu Road, Xi'an, Shaanxi 710004, P.R. China, E-mail: yatan27weihe@163.com 072018 09052018 18 1 333 341 01112017 17042018 Copyright: © He et al. 2018 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Brain ischemia leads to energy depletion, mitochondrial dysfunction and neuronal cell death. The present study was designed to identify key genes and pathways associated with brain ischemia. The gene expression profile GSE52001, including 3 normal brain samples and 3 cerebral ischemia samples, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified using the limma package. Then functional and pathway enrichment analyses were performed by the MATHT tool. Protein-protein interaction (PPI) network, module selection and microRNA (miRNA)-target gene network were constructed utilizing Cytoscape software. A total of 488 DEGs were identified (including 281 upregulated and 207 downregulated genes). In the PPI network, Rac family small GTPase 2 (RAC2) had higher degrees. RAC2 was significantly enriched in the FcγR-mediated phagocytosis pathway. miR-29A/B/C had a higher degree in the miRNA-target gene network. Insulin like growth factor 1 (Igf1) was identified as the target gene for miR-29A/B/C. RAC2 may function in brain ischemia through mediating the FcγR-mediated phagocytosis pathway. Meanwhile, miR-29A/B/C and their targets gene Igf1 may serve important roles in the development and progression of brain ischemia.

 brain ischemia differentially-expressed genes protein-protein interaction network microRNA-target gene network Introduction

Ischemic stroke, the third leading cause of death, leads to neuronal cell death by necrosis or apoptosis, mitochondrial dysfunction, energy depletion, and its complications such as coma, and hemiplegia (1–3). The study shows that the incidence of ischemic stroke decrease over time among men, but it is stable among women (4). It is report that the incidence of ischemic stroke is 170/100 thousand in adult women (5,6), and it is 212/100 thousand in men. Notably, the incidence of ischemic stroke is 91.3/100 thousand-263.1/100 thousand in China (7). Usually, the middle cerebral artery is related to ischemic stroke (8). Despite some clot lysing drugs have applied to ameliorate cerebral ischemic according to clinical experience, the treatment efficacy is limited by the narrow therapeutic safety and time window (1,9). In addition, reperfusion also aggravates brain injury, such as neuronal apoptosis, reactive oxygen species overproduction, and neuro-inflammation (10,11). Thus, it is critical to explore the novel therapeutic agents and targets of ischemic stroke.

Currently, numerous studies involved the pathophysiological of ischemic stroke are performed. For example, the Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/STAT3) signaling pathway is suggested to play a vital role in central nervous system (12). In this pathway, hypoxia preconditioning (13), rhEPO (14), IL-6 (15) can activate JAK2-STAT3 pathway and promote neurological recovery (13,16). In addition, the suppressor of cytokine signaling (SOCS) family of proteins, including SOCS1 and SOCS3, can suppress cytokine activity by interacting with JAK (17), indicating JAK2/STAT3 pathway is associated with cerebral ischemia reperfusion injury. Moreover, studies show that Rac family small GTPase 2 (Rac2), a well-studied small GTPase, has the effect on hematopoietic and endothelial cell integrin and immunoreceptor signaling (18,19). Joshi et al demonstrate that Rac2 is related to macrophage autonomous process, which can control tumor growth (20). However, the relationship between Rac2 and cerebral ischemia need to be further investigated.

MicroRNA (miRNA), a small non-coding RNA molecule, has important functions in the RNA silencing and post-transcriptional regulation of gene expression (21). The study shows that miRNAs are involved in neuro protection, ischemia, and injury (21). Previous study indicated that miR-29a had the protective effect on reperfusion injury by targeting a pro-apoptotic family member (22). However, potential gene markers related to brain ischemia based on gene or miRNA expression remains unclear.

The GSE52001 is obtained on Agilent Array platform and firstly analyzed by Lai et al (23). However, based on the huge information of gene expression profile, the data about the role of potential gene markers in cerebral ischemia are limited. In the present study, a bioinformatics study was performed based on the microarray data deposited by Lai et al (23). On the basis of differentially expressed gene (DEGs) between sham brain samples (sham group) and cerebral ischemia brain samples (ischemia brain group), the function and pathways analyses were investigated. Protein-protein interaction (PPI) network analysis was also conducted. Then the analysis for potential miRNA-target regulation in the process of glioma was performed. We expected to explore a detailed mechanism of transcriptional regulation in the cerebral ischemia, and provide a novel strategy for cerebral ischemia therapy.

 Materials and methods <sec> <title>Microarray data

The gene expression profiling GSE52001 was downloaded from the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) database (24), which was based on the platform of GPL14746 Agilent-028282 Whole Rat Genome Microarray 4×44 K V3.0 expression beadchip. The organism of this dataset was Rattus norvegicus, including 3 normal brain samples (Sham brain, SB, SB1, SB2, SB3) and 3 brain samples of cerebral ischemia (ischemia brain, IB, IB1, IB2, IB3) (23). The IB samples were collected as follows: Male Sprague Dawley rats (220±20 g; 7–8 weeks old) were anesthetized with 10% chloral hydrate (3 ml/kg). Then a silicone-coated nylon monofilament was inserted from the left common carotid artery to the origin of the middle cerebral artery. After 2 h of occlusion, reperfusion was induced by withdrawing the filament. Sham animals were operated on in the same manner except that the middle cerebral artery was not occluded. The animal experiments were conducted in accordance with the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals. The study of Lai et al was approved by the Animal Care and Use Committee of Fujian University of Traditional Chinese Medicine (23).

 Pretreatment and differential analysis

The robust multi-array average (RMA) method in limma package (http://www.bioconductor.org/packages/2.9/bioc/html/limma.html) (25) was applied to preprocess the raw CEL data by performing background correction, data normalization, conversion of original data, and quartile data normalization. Then the DEGs were identified by the non-paired t-test in limma package (25). Here, the adjusted P-value <0.05 and |log fold change (FC) |≥1 were set as the threshold value. Finally, the heat map for DEGs was generated via the Pheatmap package (https://CRAN.R-project.org/package=pheatmap) (26) in R (version 3.3.2).

 Functional and pathway enrichment analyses

Gene Ontology (GO) (http://www.geneontology.org) analysis (27) is used for analyzing the functions of a large number of genes. The Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.ad.jp/kegg/) (28,29) pathway is the major recognized pathway-related database which contains varieties of biochemical pathways (29). Multifaceted Analysis Tool for Human Transcriptome (MATHT) (www.biocloudservice.com) was used to perform GO term and KEGG pathway enrichment analyses for the DEGs. The setting of cut-off value was P-value <0.05.

 PPI network and module analyses

The Search Tool for the Retrieval of Interacting Genes (STRING) (version 10.0) (30) (http://www.string-db.org/) is an online database providing experimental and predicted PPI information. Here, the STRING database (30) was applied to analyze the PPIs among the proteins encoded by the DEGs. The parameter was set as medium confidence score >0.4. Then the PPI networks for the upregulated genes and the downregulated genes were separately visualized by Cytoscape software (version 3.2.0) (http://www.cytoscape.org/) (31), and node degrees were determined. In addition, the significant modules were obtained using the MCODE plug-in (http://apps.cytoscape.org/apps/mcode) (32) in Cytoscape software. Additionally, the KEGG pathway enrichment analysis for nodes in the significant modules was performed using MATHT tool.

 miRNA-target gene regulatory network analyses

With the discovery of RNA interference (RNAi), the function of noncoding RNAs in gene expression and regulation was widely focused (33). MiRNAs regulate the expression of genes by interacting with their target genes at the post transcription stage (33). In the present study, the miRNAs associated with DEGs were searched utilizing WebGestal (http://www.webgestalt.org/option.php) (34,35) online tool, and miRNA-DEG regulatory network was visualized by Cytoscape software (31).

 Results <sec> <title>DEGs and clusters

In total, 488 DEGs were identified, including 281 upregulated and 207 downregulated DEGs. Thereafter, the 488 DEGs and 6 samples were clustered, and DEGs could well differentiate the IB samples from the SB controls (Fig. 1).

 Functional and pathway enrichment analyses

The enriched GO terms and KEGG pathways for DEGs were identified. According to the P-values (ascending sort), the top 5 enriched terms were exhibited in Fig. 2. The upregulated genes were significantly enriched in the functions of aging (BP, P=1.69×10−7), extracellular exosome (CC, P=3.58×10−11), extracellular space (CC, P=5.70×10−14), protein homodimerization activity (MF, P=6.23×10−4), and identical protein binding (MF, P=3.80×10−3) (Fig. 2A). While the downregulated genes were dramatically enriched in the functions of response to drug (BP, P=3.20×10−3), extracellular space (CC, P=1.25×10−2), transcriptional activator activity (MF, P=1.03×10−4), and RNA polymerase II core promoter proximal region sequence-specific binding (MF, P=1.03×10−4) (Fig. 2B). Besides, the significantly enriched KEGG pathways were presented in Table I. For the upregulated genes, the significantly enriched KEGG pathways mainly include staphylococcus aureus infection (pathway, P=8.15×10−10), pertussis (pathway, P=2.64×10−7), and complement and coagulation cascades (pathway, P=2.36×10−6). There were only 2 significantly enriched KEGG pathways for the downregulated genes, which include African trypanosomiasis and Aldosterone synthesis (pathway, P=4.70×10−3), and secretion (pathway, P=3.96×10−2) (Table I).

 PPI network and module analyses

The PPI network with 295 nodes and 827 edges was constructed (Fig. 3). Upregulated gene with higher node degree were Rac2, angiotensinogen (Agt), integrin β2 (Itgb2), protein tyrosine phosphatase, receptor type, C (Ptprc), protein tyrosine phosphatase, non-receptor type 6 (Ptpn6), and hematopoietic cell kinase (Hck). Downregulated genes with higher degrees were Fos, Hras, and Junb. The degree of top 20 DEGs was listed in the Table II. In this study, 3 significant modules were obtained by MCODE plug-in, which included module A (10 nodes and 45 edges), module B (12 nodes and 33 edges), and module C (5 nodes and 10 edges) (Fig. 4).

In addition, 2 KEGG pathways were significantly enriched in module A, including neuroactive ligand-receptor interaction (P=9.62×10−4), and inflammatory mediator regulation of TRP channels (P=3.03×10−3). Meanwhile, 8 KEGG pathways were dramatically enriched in module B, such as FcγR-mediated phagocytosis (P=1.60×10−6), B cell receptor signaling pathway (P=5.40×10−5), and natural killer cell mediated cytotoxicity (P=1.54×10−4) (Table III). However, no pathways were enriched for the nodes in module C.

 miRNA-target regulatory network analysis

A total of 58 DEGs, 13 miRNAs, and 128 edges were contained in the miRNA-DEG regulatory network (Fig. 5). The nodes with top 10 degrees were listed in Table IV. Among them, the degrees of miR-29A, miR-29B and miR-29C were higher than other miRNAs in the miRNA-target regulatory network. Besides, target genes with higher degrees such as Mycn, Plcb1, Igf1 were shown in Fig. 5.

 Discussion

In the present study, a total of 488 DEGs were identified, including 281 upregulated and 207 downregulated DEGs. Rac2, with higher degree in the PPI network, was associated with FcγR-mediated phagocytosis pathway. In the miRNA-target gene network, the degrees of miR-29A, miR-29B and miR-29C were higher than other miRNAs. Notably, the target gene Igf1 was regulated by miR-29A, miR-29B and miR-29C.

Rac2, a member of Rac sub-class 3 proteins (including Rac1, Rac2 and Rac3), is a well-studied small GTPase (36,37). Among sub-class proteins, there are 92% sequence identity between Rac1 and Rac2, 83% identity between Rac2 and Rac3, and 77% identity between Rac1 and Rac3 (18). Rac2 is only expressed in hematopoietic and endothelial cells, while Rac1 and Rac3 are comprehensively expressed in mammalian systems (36–38). Rac2, the hematopoietic specific GTPase, plays an obligate role in endothelial integrin signaling and the postnatal neovascularization response in vivo (19). Additionally, some studies demonstrated that Rac2 regulates FcγR-mediated phagocytosis (39–41). Consistent with Yang et al (42), this study also shows that RAC2 is related to FcγR-mediated phagocytosis pathway in brain ischemia by bioinformatics methods. Yang et al (42) point that pathological nerve pain may be related to immune dysfunctions. Here, our results showed that RAC2 gene and FcγR-mediated phagocytosis pathway may have vital effect on the progress of pathological nerve pain in nervous system. Therefore, RAC2 may play an important role in the development of brain ischemia by mediating FcγR-mediated phagocytosis pathway.

In recent years, Kriegel et al reveal that the miR-29 family, including miR-29A, miR-29B-1, miR-29B-2 and miR-29C (43), is found to be enriched in astrocytes (44). Previous study also shows that miR-29 family is downregulated in cortex (45), but is upregulated in hippocampus after focal ischemia (46). miR-29A/B-1 is reported in Alzheimer's disease (47). Ouyang et al uncover that miR-29A is significantly upregulated in astrocytes, and regulates ischemic injury by BH3-only protein PUMA (22). The study shows that miR-29B loss at the infarct site is an important contributor to stroke lesion by 12-lipoxygenase pathway (48). Downregulated miR-29C promotes ischemic brain damage by its target gene DNMT3a. REST, an upstream transcriptional controller of miR-29C, can impede miR-29C downregulation and ischemic neuronal death by reducing REST induction (49). In this study, the degrees of miR-29A, miR-29B and miR-29C were higher than other miRNAs in the miRNA-target regulatory network. These results suggest that miR-29A/B/C may be novel biomarkers for the protection of brain ischemia injury. Intriguingly, these results were also in accordance with previous studies (22,48,49). In addition, our results indicated that miR-29A/B/C regulated the upregulated target gene Igf1 in brain ischemia. Nicholas et al find that Igf1 and Il1RAP are direct targets of miR-29 in biliary atresia (50). Therefore, we speculated that Igf1 targeted by miR-29A/B/C may have significant effect in brain ischemia. However, the detailed regulatory relationship between Igf1 and miR-29A/B/C in brain ischemia is not validated.

In conclusion, this study indicated that RAC2 may function in brain ischemia through the FcγR-mediated phagocytosis pathway. Meanwhile, miR-29A/B/C and their target gene Igf1 may have critical roles in brain ischemia. This study provides new insights into the molecular mechanisms for the progression of brain ischemia and suggests directions for future study. However, it is essential for verifying these results by the experiment in the future.

 Acknowledgements

Not applicable.

 Funding

This study was supported by the National Natural Science Foundation of China (grant no. 81270435).

 Availability of data and materials

The datasets used and/or analyzed during the present study are available from the corresponding author on reasonable request.

 Authors' contributions

JH and YG conceived and designed the research. YZ acquired the data. WP, GW and XL analyzed and interpreted the data. JH and HY performed the statistical analysis. YG obtained the funding. JH drafted the manuscript. YG, GW and XL revised the manuscript. All authors have read and approved the final manuscript.

 Ethics approval and consent to participate

The animal experiments were conducted in accordance with the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals. The study of Lai et al was approved by the Animal Care and Use Committee of Fujian University of Traditional Chinese Medicine.

 Consent for publication

Not applicable.

 Competing interests

The authors declare that they have no competing interest.

 References MuresanuDFBuzoianuAFlorianSIvon WildTTowards a roadmap in brain protection and recoveryJ Cell Mol Med1628612871201210.1111/j.1582-4934.2012.01605.x228635214393716 HankeyGJAnticoagulant therapy for patients with ischaemic strokeNat Rev Neurol8319328201210.1038/nrneurol.2012.7722565207 KlimkiewiczPKubsikAWoldańska-OkońskaMNDT-Bobath method used in the rehabilitation of patients with a history of ischemic strokeWiad Lek651021072012(In Polish)23289255 MalkiNKoupilIElorantaSWeibullCETiikkajaSIngelssonESparénPTemporal trends in incidence of myocardial infarction and ischemic stroke by socioeconomic position in Sweden 1987–2010PLoS One9e105279201410.1371/journal.pone.0105279251709194149372 SherzaiAZMaHHornrossPCancholaAJVoutsinasJWilleyJZGuYScarmeasNSherzaiDBernsteinLAbstract MP85: Mediterranean diet and incidence of stroke in the California teachers studyCirculation131AMP85AMP852015 KhouryJCKisselaBSucharewHAlwellKMoomawCWooDFlahertyMAdeoyeOKhatriPFerioliSAbstract 145: Is the rate of ischemic stroke incidence changing differentially over time for women and men? In: American-Heart-Association/American-Stroke-Association InternationalStroke44Suppl 1A1452013 TuXSEpidemiological studies of acute ischemic strokeChin J Clin Neurosci245945992016 LiLLiHLiMCurcumin protects against cerebral ischemia-reperfusion injury by activating JAK2/STAT3 signaling pathway in ratsInt J Clin Exp Med814985149912015266289814658870 YepesMRousselBDAliCVivienDTissue-type plasminogen activator in the ischemic brain: More than a thrombolyticTrends Neurosci324855200910.1016/j.tins.2008.09.00618963068 YangYJiangSDongYFanCZhaoLYangXLiJDiSYueLLiangGMelatonin prevents cell death and mitochondrial dysfunction via a SIRT1-dependent mechanism during ischemic-stroke in miceJ Pineal Res586170201510.1111/jpi.1219325401748 PalenciaGMedranoJÁOrtiz-PlataAFarfánDJSoteloJSánchezATrejo-SolísCAnti-apoptotic, anti-oxidant, and anti-inflammatory effects of thalidomide on cerebral ischemia/reperfusion injury in ratsJ Neurol Sci3517887201510.1016/j.jns.2015.02.04325818676 NicolasCSAmiciMBortolottoZADohertyACsabaZFafouriADournaudPGressensPCollingridgeGLPeineauSThe role of JAK-STAT signaling within the CNSJAKSTAT2e229252013240587893670265 WangGZhouDWangCGaoYZhouQQianGDecosterMAHypoxic preconditioning suppresses group III secreted phospholipase A2-induced apoptosis via JAK2-STAT3 activation in cortical neuronsJ Neurochem11410391048201020492356 ZhouTFYuJGRecombinant human erythropoietin attenuates neuronal apoptosis and cognitive defects via JAK2/STAT3 signaling in experimental endotoxemiaJ Surg Res183304312201310.1016/j.jss.2012.11.03523236988 FengQWangYIYangYNeuroprotective effect of interleukin-6 in a rat model of cerebral ischemiaExp Ther Med916951701201510.3892/etm.2015.2363261368794471691 ZhuHZouLTianJDuGGaoYSMND-309, a novel derivative of salvianolic acid B, protects rat brains ischemia and reperfusion injury by targeting the JAK2/STAT3 pathwayEur J Pharmacol7142331201310.1016/j.ejphar.2013.05.04323764464 YasukawaHMisawaHSakamotoHMasuharaMSasakiAWakiokaTOhtsukaSImaizumiTMatsudaTIhleJNThe JAK-binding protein JAB inhibits Janus tyrosine kinase activity through binding in the activation loopEMBO J1813091320199910.1093/emboj/18.5.1309100645971171221 PradipDPengXDurdenDLRac2 specificity in macrophage integrin signaling: Potential role for Syk kinaseJ Biol Chem2784166114669200310.1074/jbc.M30649120012917394 DePPengQTraktuevDOLiWYoderMCMarchKLDurdenDLExpression of RAC2 in endothelial cells is required for the postnatal neovascular responseExp Cell Res315248263200910.1016/j.yexcr.2009.04.012191232682767303 JoshiSSinghARZulcicMBaoLMesserKIdekerTDutkowskiJDurdenDLRac2 controls tumor growth, metastasis and M1-M2 macrophage differentiation in vivoPLoS One9e95893201410.1371/journal.pone.0095893247703464000195 SaugstadJAMicroRNAs as effectors of brain function with roles in ischemia and injury, neuroprotection, and neurodegenerationJ Cereb Blood Flow Metab3015641576201010.1038/jcbfm.2010.101206066862932764 OuyangYBXuLLuYSunXYueSXiongXXGiffardRGAstrocyte-enriched miR-29a targets PUMA and reduces neuronal vulnerability to forebrain ischemiaGlia6117841794201310.1002/glia.2255624038396 LaiWZhengZZhangXWeiYChuKBrownJHongGChenLSalidroside-mediated neuroprotection is associated with induction of early growth response genes (Egrs) across a wide therapeutic windowNeurotox Res28108121201510.1007/s12640-015-9529-925911293 BarrettTSuzekTOTroupDBWilhiteSENgauWCLedouxPRudnevDLashAEFujibuchiWEdgarRNCBI GEO: Mining millions of expression profiles-database and toolsNucleic Acids Res33D562D566200510.1093/nar/gki02215608262 SmythGKLimma: Linear Models for Microarray DataBioinformatics and Computational Biology Solutions Using R and BioconductorGentlemanRCareyVJHuberWIrizarryRADudoitSSpringerNew York, NY397420200510.1007/0-387-29362-0_23 KoldeRPheatmap: Pretty HeatmapsR PackageVersion 0.7. 7CRAN Repository2012 AshburnerMBallCABlakeJABotsteinDButlerHCherryJMDavisAPDolinskiKDwightSSEppigJTGene ontology: Tool for the unification of biology. The Gene Ontology ConsortiumNat Genet252529200010.1038/75556108026513037419 KanehisaMGotoSKEGG: Kyoto Encyclopedia of Genes and GenomesNucleic Acids Res282730200010.1093/nar/28.1.2710592173102409 DuJYuanZMaZSongJXieXChenYKEGG-PATH: Kyoto encyclopedia of genes and genomes-based pathway analysis using a path analysis modelMol Biosyst1024412447201410.1039/C4MB00287C24994036 SzklarczykDFranceschiniAWyderSForslundKHellerDHuerta-CepasJSimonovicMRothASantosATsafouKPSTRING v10: PRotein-protein interaction networks, integrated over the tree of lifeNucleic Acids Res43(Database Issue)D447D452201510.1093/nar/gku100325352553 ShannonPMarkielAOzierOBaligaNSWangJTRamageDAminNSchwikowskiBIdekerTCytoscape: A software environment for integrated models of biomolecular interaction networksGenome Res1324982504200310.1101/gr.123930314597658403769 BandettiniWPKellmanPManciniCBookerOJVasuSLeungSWWilsonJRShanbhagSMChenMYAraiAEMultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: A clinical validation studyJ Cardiovasc Magn Reson1483201210.1186/1532-429X-14-83231993623552709 Slezak-ProchazkaIDurmusSKroesenBJvan den BergAMicroRNAs, macrocontrol: Regulation of miRNA processingRNA1610871095201010.1261/rna.1804410204239802874160 WangJDuncanDShiZZhangBWEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013Nucleic Acids Res41W77W83201310.1093/nar/gkt439237032153692109 ZhangBKirovSSnoddyJWebGestalt: An integrated system for exploring gene sets in various biological contextsNucleic Acids Res33W741W748200510.1093/nar/gki475159805751160236 DidsburyJWeberRFBokochGMEvansTSnydermanRrac, a novel ras-related family of proteins that are botulinum toxin substratesJ Biol Chem264163781638219892674130 HaatajaLGroffenJHeisterkampNCharacterization of RAC3, a novel member of the Rho familyJ Biol Chem2722038420388199710.1074/jbc.272.33.203849252344 DiekmannDNobesCDBurbeloPDAboAHallARac GTPase interacts with GAPs and target proteins through multiple effector sitesEMBO J145297530519957489719394639 UeyamaTEtoMKamiKTatsunoTKobayashiTShiraiYLennartzMRTakeyaRSumimotoHSaitoNIsoform-specific membrane targeting mechanism of Rac during Fc gamma R-mediated phagocytosis: Positive charge-dependent and independent targeting mechanism of Rac to the phagosomeJ Immunol17523812390200510.4049/jimmunol.175.4.238116081809 YamauchiALiSJMarchalCCDinauerMCRac2-null macrophages have defects in superoxide production and Fc gamma-R-mediated phagocytosisBlood10017722002 PatelJCHallACaronEVav regulates activation of Rac but not Cdc42 during FcgammaR-mediated phagocytosisMol Biol Cell1312151226200210.1091/mbc.02-01-000211950933102263 YangYKLuXBWangYHYangMMJiangDMIdentification crucial genes in peripheral neuropathic pain induced by spared nerve injuryEur Rev Med Pharmacol Sci1821522159201425070821 KriegelAJLiuYFangYDingXLiangMThe miR-29 family: Genomics, cell biology, and relevance to renal and cardiovascular injuryPhysiol Genomics44237244201210.1152/physiolgenomics.00141.2011222146003289120 SmirnovaLGräfeASeilerASchumacherSNitschRWulczynFGRegulation of miRNA expression during neural cell specificationEur J Neurosci2114691477200510.1111/j.1460-9568.2005.03978.x15845075 DharapABowenKPlaceRLiLCVemugantiRTransient focal ischemia induces extensive temporal changes in rat cerebral microRNAomeJ Cereb Blood Flow Metab29675687200910.1038/jcbfm.2008.157191421922743462 YuanYWangJYXuLYCaiRChenZLuoBYMicroRNA expression changes in the hippocampi of rats subjected to global ischemiaJ Clin Neurosci17774778201010.1016/j.jocn.2009.10.00920080409 ShioyaMObayashiSTabunokiHArimaKSaitoYIshidaTSatohJAberrant microRNA expression in the brains of neurodegenerative diseases: miR-29a decreased in Alzheimer disease brains targets neurone navigator 3Neuropathol Appl Neurobiol36320330201010.1111/j.1365-2990.2010.01076.x20202123 KhannaSRinkCGhoorkhanianRGnyawaliSHeigelMWijesingheDSChalfantCEChanYCBanerjeeJHuangYLoss of miR-29b following acute ischemic stroke contributes to neural cell death and infarct sizeJ Cereb Blood Flow Metab3311971206201310.1038/jcbfm.2013.68236329683734770 PandiGNakkaVPDharapARoopraAVemugantiRMicroRNA miR-29c downregulation leading to de-repression of its target DNA methyltransferase 3a promotes ischemic brain damagePLoS One8e58039201310.1371/journal.pone.0058039235164283596316 HandNJHornerAMMasterZRBoatengLALeGuenCUvaydovaMFriedmanJRMicroRNA profiling identifies miR-29 as a regulator of disease-associated pathways in experimental biliary atresiaJ Pediatr Gastroenterol Nutr54186192201210.1097/MPG.0b013e318244148b221670213264748

Cluster analysis of DEGs. DEGs, differentially expressed genes; IB, ischemia brain samples; SB, sham brain samples. Cluster analysis was performed at both gene level (vertical) and sample level (horizontal).

Top 5 GO terms enriched separately. GO terms separated according to (A) upregulated and (B) downregulated genes. MF, Molecular function; BP, Biological processes; CC, Cellular components; GO, gene ontology; The horizontal axis represents the enriched GO terms The vertical axis represents the count of enriched DEGs.

Protein-protein interaction network constructed for the DEGs. The pink circle and the blue square represent upregulated genes and downregulated genes, respectively. DEGs, differentially-expressed genes.

Module analyses of protein-protein interaction network. (A) The pink circles represent upregulated DEGs, and the blue squares represent downregulated DEGs. (B) The pink circles represent upregulated DEGs, and the blue square represents downregulated DEGs. (C) The blue squares represent downregulated DEGs. DEGs, differentially-expressed genes.

miRNA-target regulatory network. The pink circle and the blue square represent upregulated genes and downregulated genes, respectively. The green triangle represents miRNA. Arrow represents the miRNA-target relationship. DEGs, differentially-expressed genes; miRNAs, microRNAs.

KEGG pathways significantly enriched by DEGs.

A, Upregulated genes
Pathway ID Pathway name Count P-value Genes
rno05150 Staphylococcus aureus infection 12 8.15×10−10 C1QA, C1QB, C5AR1, FCGR2B, C3, LOC498276, C1R, ITGB2, C2, C1S, FCGR3A, C1QC
rno05133 Pertussis 11 2.64×10−7 C1QA, C1QB, C3, PYCARD, SERPING1, C1R, ITGB2, C2, C1S, C1QC, CD14
rno04610 Complement and coagulation cascades 10 2.36×10−6 C1QA, C1QB, A2M, C5AR1, C3, SERPING1, C1R, C2, C1S, C1QC
rno04650 Natural killer cell mediated cytotoxicity 10 3.05×10−5 CD48, PTPN6, RAC2, FCER1G, ITGB2, VAV2, FCGR3A, IFNGR1, HCST, TYROBP
rno05140 Leishmaniasis 7 1.13×10−3 PTPN6, CYBA, C3, LOC498276, ITGB2, FCGR3A, IFNGR1
rno04145 Phagosome 11 1.41×10−3 RT1-A2, CYBA, RT1-A1, FCGR2B, C3, LOC498276, C1R, ITGB2, CTSS, FCGR3A, CD14
rno05322 Systemic lupus erythematosus 9 1.44×10−3 C1QA, C1QB, C3, LOC498276, C1R, C2, C1S, FCGR3A, C1QC
rno05152 Tuberculosis 10 3.01×10−3 LSP1, FCGR2B, C3, LOC498276, FCER1G, ITGB2, CTSS, FCGR3A, IFNGR1, CD14
rno04666 FcγR-mediated phagocytosis 6 1.28×10−2 PTPRC, RAC2, FCGR2B, HCK, LOC498276, VAV2
rno04142 Lysosome 7 1.81×10−2 CTSZ, GUSB, LGMN, CTSE, CTSC, CTSS, CD63
rno04380 Osteoclast differentiation 7 1.94×10−2 CYBA, FCGR2B, LOC498276, TREM2, FCGR3A, IFNGR1, TYROBP
rno04611 Platelet activation 7 2.37×10−2 P2RY12, ORAI1, TBXAS1, FERMT3, LOC498276, COL3A1, FCER1G
rno00860 Porphyrin and chlorophyll metabolism 4 2.94×10−2 GUSB, HMOX1, HEPH, CP
rno05146 Amoebiasis 6 3.70×10−2 ARG1, COL3A1, SERPINB1A, ITGB2, SERPINB1B, CD14
rno04670 Leukocyte transendothelial migration 6 4.75×10−2 CYBA, RAC2, CLDN1, ITGB2, VAV2, MMP2
B, Downregulated genes
PathwayID Pathway name Count P-value Genes
rno05143 African trypanosomiasis 4 4.70×10−3 LOC689064, HBB-B1, PLCB1, HBB
rno04925 Aldosterone synthesis and secretion 4 3.96×10−2 CYP11B1, NR4A1, PLCB1, CACNA1S

KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes.

Degree of top 20 differentially-expressed genes in the protein-protein interaction network.

A, Upregulated
Gene Degree
Rac2 35
Agt 34
Tyrobp 30
Itgb2 30
Ptprc 28
Timp1 27
Ptpn6 24
Lgals3 23
Hck 21
Igf1 20
Jak3 20
Anxa1 20
Cd53 20
Dcn 18
C1qb 16
B, Downregulated
Gene Degree
Fos 34
Hras 32
Junb 20
Sst 19
Egr1 17

Enriched pathways for the nodes in module A and B.

Pathway ID Pathway name Count P-value Genes
MEA
  rno04080 Neuroactive ligand-receptor interaction 4 9.62×10−4 P2RY6, P2RY2, LPAR4, HTR2A
  rno04750 Inflammatory mediator regulation of TRP channels 3 3.13×10−3 P2RY2, PLCB1, HTR2A
MEB
  rno04666 FcγR-mediated phagocytosis 5 1.60×10−6 PTPRC, RAC2, FCGR2B, HCK, VAV2
  rno04662 B cell receptor signaling pathway 4 5.40×10−5 PTPN6, RAC2, FCGR2B, VAV2
  rno04650 Natural killer cell mediated cytotoxicity 4 1.54×10−4 PTPN6, RAC2, ITGB2, VAV2
  rno05150 Staphylococcus aureus infection 3 1.65×10−3 C5AR1, FCGR2B, ITGB2
  rno04660 T cell receptor signaling pathway 3 6.45×10−3 PTPN6, PTPRC, VAV2
  rno04670 Leukocyte transendothelial migration 3 7.91×10−3 RAC2, ITGB2, VAV2
  rno04062 Chemokine signaling pathway 3 1.67×10−2 RAC2, HCK, VAV2
  rno04810 Regulation of actin cytoskeleton 3 2.47×10−2 RAC2, ITGB2, VAV2

Degree of top 10 miRNAs in the miRNAs-target regulatory network.

miRNA Degree
miR-29A 15
miR-29B 15
miR-29C 15
miR-27A 12
miR-27B 12
miR-26A   9
miR-26B   9
miR-377   8
miR-202   8
miR-520G   8

miRNA, microRNA.