Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2017.8367 mmr-17-03-3473 Articles Bioinformatic analysis of gene expression profiling of intracranial aneurysm BoLijuan 1 * WeiBo 2 * WangZhanfeng 2 LiChaohui 2 GaoZheng 3 MiaoZhuang 2 Department of Infection, China-Japan Union Hospital of Jilin University, Changchun, Jilin 130033, P.R. China Department of Neurosurgery, China-Japan Union Hospital of Jilin University, Changchun, Jilin 130033, P.R. China Department of Neurosurgery, Dandong First Hospital, Dandong, Liaoning 118000, P.R. China Correspondence to: Dr Zhuang Miao, Department of Neurosurgery, China-Japan Union Hospital of Jilin University, 126 Xiantai Street, Changchun, Jilin 130033, P.R. China, E-mail: miaozhzhzhuuuu@hotmail.com

Contributed equally

 032018 29122017 17 3 3473 3480 19052016 04042017 Copyright: © Bo 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.

Intracranial aneurysm (IA) is a severe clinical condition of primary concern and currently, there is no effective therapeutic reagent. The present study aimed to investigate the molecular mechanism of IA via bioinformatic analysis. Various gene expression profiles (GSE26969) were downloaded from the Gene Expression Omnibus database, including 3 IA and 3 normal superficial temporal artery samples. Firstly, the limma package in R language was used to identify differentially expressed genes (DEGs; P-value <0.01 and |log2 FC|≥1). Secondly, the database for annotation, visualization and integrated discovery software was utilized to perform pathway and functional enrichment analyses (false discovery rate ≤0.05). Finally, protein-protein interaction (PPI) network and sub-network clustering analyses were performed using the biomolecular interaction network database and ClusterONE software, respectively. Following this, a transcription factor regulatory network was identified from the PPI network. A total of 1,124 DEGs were identified, of which 989 were upregulated and 135 downregulated. Pathway and functional enrichment analyses revealed that the DEGs primarily participated in RNA splicing, functioning of the spliceosome, RNA processing and the mRNA metabolic process. Following PPI network analysis, 1 hepatocyte nuclear factor (HNF) 4A (transcription factor)-centered regulatory network and 5 DEG-centered sub-networks were identified. On analysis of the transcription factor regulatory network, 6 transcription factors (HNF6, HNF4A, E2F4, YY1, H4 and H31T) and a regulatory pathway (HNF6-HNF4-E2F4) were identified. The results of the present study suggest that activating transcription factor-5, Jun proto-oncogene, activator protein-1 transcription factor subunit, HNF6, HNF4 and E2F4 may participate in IA progression via vascular smooth muscle cell apoptosis, inflammation, vessel wall remodeling and damage and the tumor necrosis factor-β signaling pathway. However, further experimental studies are required to validate these predictions.

 intracranial aneurysm differentially expressed genes functional enrichment analysis protein-protein interaction regulatory network Introduction

Intracranial aneurysm (IA) is a cerebrovascular disorder characterized by a regional ballooning of intracranial arteries. IAs are classified into saccular aneurysms, fusiform aneurysms, and micro-aneurysms according to their size and shape. IA is a frequently occurring condition with a prevalence of 2–5% in the global population, and ~0.7–1.9% of cases transform into subarachnoid hemorrhage (SAH) (1). SAH, which constitutes 1–7% of all strokes, has high mortality and morbidity rates, even though diagnosis and treatment with modern medical technologies are available (2). Currently, imaging methods including digital subtraction angiography, computed tomography (CT) and CT angiography are the primary IA diagnostic methods, however there is no available preventive treatment of IA prior to rupture, except for surgical procedures and endovascular therapy. However, these treatment methods often result in further complications. The requirement to develop novel therapeutic strategies for IA is of primary concern, and understanding the underlying molecular mechanism of its initiation will aid this research.

Recently, numerous studies have been conducted to identify the pathogenesis of IA formation and rupture. Based on affymetrix microarray data from 3 unruptured IAs and a control superficial temporal artery, Li et al (3) demonstrated that upregulated genes are extracellular matrix (ECM)-associated genes and downregulated genes are enriched in the immune/inflammation response. Based on data obtained from 8 ruptured IAs, 6 unruptured IAs and 5 control intracranial arteries, Pera et al (4) demonstrated that upregulated genes are enriched in the immune/inflammatory response and downregulated genes are enriched in the muscle system and cell adhesion. Based on data from 11 ruptured IAs and 8 unruptured IAs, Kurki et al (5) demonstrated that upregulated genes are enriched in chemotaxis, leukocyte migration, oxidative stress, vascular remodeling, ECM degradation and the response to turbulent blood flow. In addition, Yagi et al (6) demonstrated that phosphodiesterase-4 is involved in inflammatory diseases and blocking it reduces macrophage migration and inhibits the progression of cerebral aneurysms (CA). Aoki et al (7) revealed that v-ets avian erythroblastosis virus E26 oncogene homolog EST-1 is over-expressed in vascular smooth muscle cells (VSMCs) in CA walls, and thus promotes CA progression.

The present study firstly downloaded gene expression profiling data. The raw data were then analyzed to identify differentially expressed genes (DEGs) between IA and normal superficial temporal artery tissue. Following this, pathway and functional enrichment analyses were conducted to investigate the functions of the identified DEGs. A protein-protein interaction (PPI) network was then constructed and visualized, from which a transcription factor regulatory network was identified.

 Materials and methods <sec> <title>Obtaining and preprocessing of microarray data

Gene expression profiles of IA and normal superficial temporal artery samples were obtained from the National Center of Biotechnology Information (NCBI) Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo/) database. The access number is GSE26969, and this dataset included a total of 6 samples; 3 IA and 3 normal superficial temporal artery samples. The platform used was the Affymetrix Human Genome U133 Plus 2.0 Array (GPL570; Affymetrix; Thermo Fisher Scientific, Inc., Waltham, MA, USA). Firstly, the original data at the probe symbol level were converted into expression values at the gene symbol level. For each sample, the expression values of all probes mapped to a given gene symbol were averaged. Then, missing data were inserted and median data normalization performed by using the robust multichip averaging (RMA) method (8).

 DEG screening and clustering analysis

The limma package (http://www.bioconductor.org/packages/2.9/bioc/html/limma.html) (9) in R language with multiple testing correction based on the Benjamini & Hochberg method (10) was employed to identify DEGs between IA and normal superficial temporal artery samples. P-value <0.01 and |log2 fold change (FC) |≥1 were used as the thresholds for this analysis. Then, the pheatmap package (https://cran.rproject.org/web/packages/pheatmap/index.html) (11) in R language was used to cluster the identified DEGs.

 Enrichment analysis of DEGs

To study DEGs at a functional level, GO (Gene Ontology) functional enrichment analysis (12) and KEGG (Kyoto Encyclopedia of Gene and Genomes) pathway enrichment (13) analysis were performed using the online biological tool, database for the annotation, visualization and integrated discovery (DAVID version 6.7; https://david-d.ncifcrf.gov/) (14). DAVID software has been extensively used to identify biological processes involving a given list of genes. In the present study, fold change discovery (FDR) ≤0.05 was set as the cut-off criterion for enrichment analysis.

 PPI network construction

PPIs are important for all biological processes. The present study constructed a PPI network among DEGs based on the biomolecular interaction network (BIND) database (http://binddb.org) (15). BIND is available for querying, viewing and submitting records regarding molecular interactions, complexes and pathway information. Additionally, ClusterONE software (http://www.paccanarolab.org/cluster-one/) (16) was used to perform sub-network clustering analysis and then GO functional enrichment analysis was then performed for DEGs in each sub-network.

 Results <sec> <title>Data preprocessing and DEG screening

Data prior to and following normalization are presented in Fig. 1. Following normalization, gene expression values of differing samples ranged almost uniformly, which indicated that all of the 6 chips were available for further analysis. Following DEG screening using the limma package, a total of 1124 DEGs (|log2 FC|≥1 and P-value ≤0.01) were identified, of which 989 DEGs were upregulated and 135 DEGs downregulated (Fig. 2A). In addition, the top 5 upregulated and top 5 downregulated DEGs are presented in Table I.

 Clustering analysis

A clustering analysis was performed using the pheatmap package in R language and 1124 DEGs were divided into two major clusters (Fig. 2B). The DEGs in the upper cluster were significantly downregulated, whereas the DEGs in the lower cluster were significantly downregulated. The separation of the clusters into two differing categories indicated that the identified DEGs were significantly differentially expressed between IA and normal superficial temporal artery samples.

 Enrichment analysis of DEGs

KEGG pathway and GO functional enrichment analyses demonstrated that 1124 DEGs were enriched in 1 KEGG pathway (hsa03040: spliceosome, FDR≤0.05) and 20 GO functions (FDR≤0.05; Fig. 3). These GO functions were primarily associated with RNA splicing and the spliceosome, RNA processing and the mRNA metabolic process.

 PPI network construction

A PPI network was then constructed using the BIND database (Fig. 4), in which a hepatocyte nuclear factor (HNF) 4A (transcription factor)-centered regulatory network was identified. In addition, a sub-network clustering analysis was performed using ClusterONE software and 5 DEG-centered sub-networks were identified (Fig. 5A-E) and significantly enriched in 3 GO functions (Table II). The central DEGs observed in the sub-networks were origin recognition complex subunit 2 (ORC2L), activating transcription factor 5 (ATF5), MAF BZIP transcription factor K (MAFK) and Jun proto-oncogene, activator protein-1 transcription factor subunit (JUN).

 Transcription factor regulatory network construction

The interactions between DEGs and transcription factors were identified from the PPI network, and presented in Fig. 6. A total of 6 transcription factors (HNF6, HNF4A, E2F4, YY1, H4 and H31T) and 24 DEGs (5 upregulated and 19 downregulated genes) were involved in the transcription factor regulatory network. It was observed that HNF4A exhibited the ability to regulate the DEGs, including prolylcarboxypeptidase (PRCP) and caspase-8 and Fas associated via death domain-like apoptosis regulator (CFLAR).

 Discussion

IA is a severe clinical condition and as of yet, there is no current effective therapy. In the present study, bioinformatic analyses were conducted to investigate the potential molecular mechanism underlying the occurrence of IA. Consequently, 1,124 DEGs were identified between IA and normal superficial temporal artery tissue. The DEGs primarily participated in RNA splicing and the spliceosome, RNA processing, and mRNA metabolic processing. Following this, a PPI network was constructed and 5 DEG-centered sub-networks were revealed, involving 4 DEGs (ORC2L, ATF5, MAFK and JUN). Furthermore, a transcription factor regulatory network was abstracted, and 6 transcription factors (HNF6, HNF4A, E2F4, YY1, H4 and H31T) and 24 DEGs (including PRCP and CFLAR) were involved.

From the results of the sub-pathway clustering analysis, the PPI network was divided into 5 DEG-centered sub-networks, involving 4 DEGs (ORC2L, ATF5, MAFK and JUN). Various DEGs identified included ATF5 (sub-network 2), which is a stress-response transcription factor and is highly expressed in various tumors (17). CCAAT/enhancer binding protein β (CEBPB) interacts with ATF5 and encodes a basic-region leucine zipper transcription factor (bZIP). CEBPB participates in the regulation of the immune and inflammatory response, and thus may participate in the pathogenesis of IA (18,19). Therefore, ATF5 may participate in IA via its interaction with CEBPB. Furthermore, JUN (sub-network 4 and 5) is a bZIP transcription factor involved in numerous cellular processes, such as cell growth and apoptosis (20), and tumor necrosis factor (TNF)-α has been demonstrated to stimulate prolonged activation of JUN and interleukin (IL)-1β gene expression (21). TNF-α is known to activate matrix metalloproteinase which is important in IA pathology (22,23). In addition, Moriwaki et al (24) demonstrated that IL-1β is important for the progression of IA via the induction of VSMC apoptosis. Therefore, JUN may participate in the progression of IA by responding to TNF-α and IL-1β.

The HNF4A-centered transcription factor regulatory network led to the identification of the regulatory pathway HNF6-HNF4A-E2F4. HNF6 encodes a transcription factor that has been demonstrated to influence a variety of cellular processes, including cell proliferation, cell-matrix adhesion and inflammation (25). The pathological pathway for IA involves endothelial dysfunction/injury, the mounting inflammatory response, VSMC phenotypic modulation, extra-cellular matrix remodeling, subsequent cell death and vessel wall degeneration (26). Therefore, HNF6 may participate in IA via induction of inflammation. HNF4A encodes a transcription factor that participates in lipid metabolic processes. Frösen et al (27) demonstrated that accumulation of lipids and oxidation in the saccular intracranial aneurysm (sIA) wall are associated with its degeneration, resulting in fatal SAH (27). Therefore, HNF4A may participate in the progression of IA via regulation of vessel wall remodeling and damage. Furthermore, E2F4 encodes E2F transcription factor 4 which participates in the cell cycle and tumor suppression. E2F4 participates in the transforming growth factor (TGF)-β signaling pathway which is involved in cell differentiation, cell growth and apoptosis. It has previously been reported that the TGF-β signaling pathway is important in IA (28). Therefore, E2F4 may be involved in IA progression via its participation in the TGF-β signaling pathway.

Of the DEGs regulated by HNF4A, PRCP encodes an enzyme which cleaves COOH-terminal amino acids linked to proline, in peptides including angiotension II (29). Angiotensin II, a part of the renin-angiotensin system (RAS), is a potent vasoconstrictor and pro-inflammatory stimulant. Shoja et al (30) demonstrated that RAS is associated with the pathogenesis of IA. Therefore, PRCP was predicted to participate in IA progression via the regulation of RAS. In addition, CFLAR encodes the caspase-8 and Fas associated via death domain-like apoptosis regulator (31), involved in the TNF-α/nuclear factor (NF)-κB signaling pathway. NF-κB recruits macrophages into the vessel walls involved in IA, initiating inflammation and participating in aneurysm formation (32). Therefore, CFLAR was hypothesized to participate in IA progression via the TNF-α/NF-κB signaling pathway.

In conclusion, the results of the present study revealed that inflammation may be an important factor in AI development. In order to elucidate the pathology of IA, gene expression profiling data of 6 samples were downloaded and analyzed. Following a sub-network clustering analysis, 5 sub-networks were identified, involving 4 DEGs (ORC2L, ATF5, MAFK and JUN). ATF5 and JUN served crucial roles in the immune and inflammatory response and VSMC apoptosis during AI development. In addition, the regulatory pathway HNF6-HNF4-E2F4 was identified as an important mechanism involved in the development of AI via the regulation of inflammation, the TGF-β signaling pathway, the renin-angiotensin system and the TNF-α/NF-κB signaling pathway. However, these results were derived from in silico analysis; thus, further experimental studies are required to verify the results of the present study.

 References AlgVSSofatRHouldenHWerringDJGenetic risk factors for intracranial aneurysms: A meta-analysis in more than 116,000 individualsNeurology8021542165201310.1212/WNL.0b013e318295d751237335523716358 BakkeSJLindegaardKFSubarachnoid haemorrhage-diagnosis and managementTidsskr Nor Laegeforen127107410782007(In Norwegian)17457396 LiLYangXJiangFDustingGJWuZTranscriptome-wide characterization of gene expression associated with unruptured intracranial aneurysmsEur Neurol62330337200910.1159/00023691119752560 PeraJKorostynskiMKrzyszkowskiTCzopekJSlowikADziedzicTPiechotaMStachuraKMoskalaMPrzewlockiRSzczudlikAGene expression profiles in human ruptured and unruptured intracranial aneurysms: What is the role of inflammation?Stroke41224231201010.1161/STROKEAHA.109.56200920044533 KurkiMIHakkinenSKFrosenJTulamoRvon und zu FraunbergMWongGTrompGNiemeläMHernesniemiJJääskeläinenJEYlä-HerttualaSUpregulated signaling pathways in ruptured human saccular intracranial aneurysm wall: An emerging regulative role of Toll-like receptor signaling and nuclear factor-κB, hypoxia-inducible factor-1A and ETS transcription factorsNeurosurgery6816661675201110.1227/NEU.0b013e318210f001 YagiKTadaYKitazatoKTTamuraTSatomiJNagahiroSIbudilast inhibits cerebral aneurysms by down-regulating inflammation-related molecules in the vascular wall of ratsNeurosurgery66551559201010.1227/01.NEU.0000365771.89576.7720124930 AokiTKataokaHNishimuraMIshibashiRMorishitaRMiyamotoSEts-1 promotes the progression of cerebral aneurysm by inducing the expression of MCP-1 in vascular smooth muscle cellsGene Ther1711171123201010.1038/gt.2010.6020428211 IrizarryRAHobbsBCollinFBeazer-BarclayYDAntonellisKJScherfUSpeedTPExploration, normalization and summaries of high density oligonucleotide array probe level dataBiostatistics4249264200310.1093/biostatistics/4.2.24912925520 GentlemanRCCareyVJBatesDMBolstadBDettlingMDudoitSEllisBGautierLGeYGentryJBioconductor: Open software development for computational biology and bioinformaticsGenome Biol5R80200410.1186/gb-2004-5-10-r8015461798545600 ReinerAYekutieliDBenjaminiYIdentifying differentially expressed genes using false discovery rate controlling proceduresBioinformatics19368375200310.1093/bioinformatics/btf87712584122 KoldeRPheatmap: Pretty heatmaps. R package version 0.612012 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 Huang daWShermanBTLempickiRASystematic and integrative analysis of large gene lists using DAVID bioinformatics resourcesNat Protoc44457200910.1038/nprot.2008.21119131956 BaderGDBetelDHogueCWBIND: The biomolecular interaction network databaseNucleic Acids Res31248250200310.1093/nar/gkg05612519993165503 NepuszTYuHPaccanaroADetecting overlapping protein complexes in protein-protein interaction networksNat Methods9471472201210.1038/nmeth.1938224264913543700 GreeneLALeeHYAngelastroJMThe transcription factor ATF5: Role in neurodevelopment and neural tumorsJ Neurochem1081122200910.1111/j.1471-4159.2008.05749.x19046351 SchoenenHHuberASondaNZimmermannSJantschJLepeniesBBronteVLangRDifferential control of mincle-dependent cord factor recognition and macrophage responses by the transcription factors C/EBPβ and HIF1αJ Immunol201410.4049/jimmunol.130159325156364 ChalouhiNAliMSJabbourPMTjoumakarisSIGonzalezLFRosenwasserRHKochWJDumontASBiology of intracranial aneurysms: Role of inflammationJ Cereb Blood Flow Metab3216591676201210.1038/jcbfm.2012.84227813303434628 WisdomRJohnsonRSMooreCc-Jun regulates cell cycle progression and apoptosis by distinct mechanismsEMBO J18188197199910.1093/emboj/18.1.18898780621171114 BrennerDAO'HaraMAngelPChojkierMKarinMProlonged activation of jun and collagenase genes by tumour necrosis factor-alphaNature337661663198910.1038/337661a02537468 PolavarapuRGongoraMCWinklesJAYepesMTumor necrosis factor-like weak inducer of apoptosis increases the permeability of the neurovascular unit through nuclear factor-kappa B pathway activationJ Neurosci251009410100200510.1523/JNEUROSCI.3382-05.200516267216 ThomsonEMWilliamsAYaukCLVincentROverexpression of tumor necrosis factor-α in the lungs alters immune response, matrix remodeling and repair and maintenance pathwaysAm J Pathol18014131430201210.1016/j.ajpath.2011.12.02022322299 MoriwakiTTakagiYSadamasaNAokiTNozakiKHashimotoNImpaired progression of cerebral aneurysms in interleukin-1beta-deficient miceStroke37900905200610.1161/01.STR.0000204028.39783.d916439700 WangKHoltermanAXPathophysiologic role of hepatocyte nuclear factor 6Cell Signal24916201210.1016/j.cellsig.2011.08.00921893194 JayaramanTPagetAShinYSLiXMayerJChaudhryHNiimiYSilaneMBerensteinATNF-alpha-mediated inflammation in cerebral aneurysms: A potential link to growth and ruptureVasc Health Risk Manag4805817200810.2147/VHRM.S2700190659972597764 FrösenJTulamoRHeikuraTSammalkorpiSNiemeläMHernesniemiJLevonenALHörkköSYlä-HerttualaSLipid accumulation, lipid oxidation and low plasma levels of acquired antibodies against oxidized lipids associate with degeneration and rupture of the intracranial aneurysm wallActa Neuropathol Commun171201310.1186/2051-5960-1-71242526583893371 RuigrokYMTanSMedicJRinkelGJWijmengaCGenes involved in the transforming growth factor beta signalling pathway and the risk of intracranial aneurysmsJ Neurol Neurosurg Psychiatry79722724200810.1136/jnnp.2007.12804118487557 TanFMorrisPWSkidgelRAErdosEGSequencing and cloning of human prolylcarboxypeptidase (angiotensinase C). Similarity to both serine carboxypeptidase and prolylendopeptidase familiesJ Biol Chem268166311663819938344943 ShojaMMAgutterPSTubbsRSPaynerTDGhabiliKCohen-GadolAAThe role of the renin-angiotensin system in the pathogenesis of intracranial aneurysmsJ Renin Angiotensin Aldosterone Syst12262273201110.1177/147032031038784521270100 DuttonAYoungLSMurrayPGThe role of cellular FLICE inhibitory protein (c-FLIP) in the pathogenesis and treatment of cancerExpert Opin Ther Targets102735200610.1517/14728222.10.1.2716441226 AokiTKataokaHMorimotoMNozakiKHashimotoNMacrophage-derived matrix metalloproteinase-2 and −9 promote the progression of cerebral aneurysms in ratsStroke38162169200710.1161/01.STR.0000252129.18605.c817122420

Boxplot of data prior to and following normalizationn.

Volcano plot and heatmap of differentially expressed genes. (A) Volcano plot and (B) heatmap of differentially expressed genes following clustering analysis.

Enrichment analysis of differentially expressed genes. FDR, false discovery rate; GO, gene ontology.

Protein-protein interaction network. A representative figure demonstrating the large scale and connectivity of the protein-protein interaction network. Red squares represent DEGs and blue circles represent genes that interact with the DEGs. Red arrows represent interactions between transcription factors and target genes and blue lines represent interactions between genes. DEGS, differentially expressed genes.

Sub-networks of differentially expressed genes extracted from the protein-protein interaction network. Sub-networks (A) 1, (B) 2, (C) 3, (D) 4 and (E) 5. Red squares represent DEGs, blue circles represent genes that interact with the DEGs and blue lines represent interactions between genes. DEGs, differentially expressed genes.

Transcription factor regulatory network. Red squares represent upregulated DEGs, green squares represent downregulated DEGs and blue circles represent genes that interact with the DEGs. Red arrows represent interactions between transcription factors and target genes and blue lines represent interactions between genes. DEGs, differentially expressed genes.

Top 5 upregulated and downregulated differentially expressed genes.

Gene Average expression value Log2 FC P-value Expression alteration
PLN     8.622363 −5.62222 3.47×10−5 Downregulated
ADH1C     8.253548 −5.56344 4.51×10−3 Downregulated
PLN     7.481133 −5.18628 1.00×10−3 Downregulated
MYL9     9.273575 −4.81707 3.61×10−3 Downregulated
SORBS1     8.440836 −4.69688 4.44×10−3 Downregulated
MMP16     7.516438 2.166793 1.17×10−3 Upregulated
SOX4     9.123515 2.227376 9.69×10−4 Upregulated
NUFIP2     8.381299 2.294707 9.57×10−3 Upregulated
TWIST1     6.786085 2.636303 1.01×10−5 Upregulated
COL5A2 10.29836     3.242704 2.12×10−4 Upregulated

FC, fold change.

Sub-networks from protein-protein interaction network.

Sub-network 1 Sub-network 2 Sub-network 3 Sub-network 4 Sub-network 5
Nodes 8 14 8 15 15
Density 0.679 0.615 0.607 0.6 0.6
Quality 0.792 0.560 0.708 0.525 0.496
P-value 4.552×10−4 5.069×10−4 5.367×10−4 1.00×10−3 3.00×10−3
Genes MCM10 CEBPE CDK6 MAFK MAFK
MCM2 DDIT3 MCM10 JUN JUN
CDC6 CEBPG MCM2 DDIT3 DDIT3
ORC2L CEBPB CDC6 CEBPG CEBPG
ORC1 DBP ORC2L MAFG MAFG
ORC5 CEBPD ORC1 CREB3 NFAC2
CDC5L ATF5 ORC5 FOSL1 FOSL1
ORC4 CEBPA ORC4 FOS FOS
CREB3 FOSL2 FOSL2
HLF ATF3 ATF3
ATF3 JUNB JUNB
BATF3 BATF3 BATF3
ATF7 ATF7 ATF7
BATF ATF2 ATF2
BACH2 BATF
Top significant GO:0006260~ GO:0046983~ GO:0006260~ GO:0043565~ GO:0043565~
GO term DNA replication protein dimerization DNA replication sequence-specific sequence-specific
activity DNA binding DNA binding

GO, gene ontology.