Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2017.7058 mmr-16-04-3777 Articles Integrated analysis of genome-wide gene expression and DNA methylation microarray of diffuse large B-cell lymphoma with TET mutations LiuPengfei 1 * JiangWenhua 2 * ZhaoJinkun 3 ZhangHuilai 1 Department of Lymphoma, Sino-US Center of Lymphoma and Leukemia, Tianjin Medical University Cancer Institute and Hospital, National Clinical Research Center for Cancer, Key Laboratory of Cancer Prevention and Therapy of Tianjin, Tianjin Clinical Research Center for Cancer, Tianjin 300162, P.R. China Department of Radiotherapy, The Second Hospital of Tianjin Medical University, Tianjin 300162, P.R. China Department of Radiology, Tianjin Medical University Cancer Institute and Hospital, National Clinical Research Center for Cancer, Key Laboratory of Cancer Prevention and Therapy, Tianjin 300162, P.R. China Correspondence to: Dr Huilai Zhang, Department of Lymphoma, Sino-US Center of Lymphoma and Leukemia, Tianjin Medical University Cancer Institute and Hospital, National Clinical Research Center for Cancer, Key Laboratory of Cancer Prevention and Therapy of Tianjin, Tianjin Clinical Research Center for Cancer, Huan-Hu-Xi Road, Ti-Yuan-Bei, He Xi, Tianjin 300162, P.R. China, E-mail: snowlpf@163.com

Contributed equally

 102017 21072017 16 4 3777 3782 24062016 20042017 Copyright: © Liu et al. 2017 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.

Diffuse large B-cell lymphoma (DLBCL), the most frequently occurring type of lymphoid malignancy, has been demonstrated to be associated with mutations of Ten-Eleven Translocation (TET). In order to explore the association between DLBCL and TET mutations, the present study analyzed the gene expression and methylation profiles in human DLBCL biopsy tissues with wildtype and mutated TET2. The microarray dataset GSE37365, containing two subseries: the genome-wide gene expression dataset GSE37362 and the DNA methylation microarray dataset GSE37363, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified using the limma package of R. Furthermore, differentially methylated sites and differentially methylated regions were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed via GO stats and GSEABase packages respectively. Finally, the Pathview package was used to construct the network of enriched pathways. A total of 198 DEGs (106 up- and 92 downregulated) were obtained. A total of 602 shared differentially methylated genes (DMGs) were identified according to differentially methylated levels. A total of 12 overlapping genes were identified in DEGs and DMGs. It was observed that 9 of the 12 overlapped genes were downregulated and hypermethylated, with 24 GO terms and one KEGG pathway significantly enriched. The results of the present study demonstrated that the genes cryptochrome circadian clock 1, zinc finger protein (ZNF) interacting with K protein 1, ZNF134, ZNF256 and ZNF615, which were hypermethylated and downregulated in DLBCL patients with TET2 mutations, were the key genes in the association between DLBCL and TET mutations. These genes may act as potential biomarkers for the diagnosis of DLBCL in the future.

 differentially expressed genes diffuse large B-cell lymphoma methylation Ten-Eleven Translocation Introduction

Diffuse large B-cell lymphoma (DLBCL) is the most frequently occurring type of lymphoid malignancy in adults worldwide. It arises from normal antigen-exposed B cells which have migrated to germinal centers in secondary lymphoid organs (1). DLBCL is a clinically heterogeneous disease and ~40% of patients may be cured by anthracycline-based chemotherapy (2). Therefore, the early diagnosis and treatment of DLBCL is critical.

A previous study indicated that specific genetic and epigenetic cancer characteristics were detectable in the plasma, serum and tissue of patients, and these were able to be used for the early detection and diagnosis of cancer (3). Ten-Eleven Translocation (TET) hydroxylases are able to convert 5-methylcytosine to 5-hydroxymethylcytosine and induce DNA demethylation (4). The use of a mutation screening assay and the genome-wide profiling has verified that TET2 mutations may result in the aberrant methylation of various genes (5). Microarray technology has contributed to the development of life sciences, making it possible to obtain large quantities of genes, microRNAs, or DNA methylated sites for further analysis of cancer at the molecular level (6,7). In addition, combining microarray data and statistical methods may be used to detect influential genes in cancers (8). Although TET mutations were considered to be associated with DLBCL, the mechanism of DNA methylation in DLBCL remains to be elucidated. Novel diagnostic targets and treatment strategies for this disease are urgently required and of primary concern.

The present study compared the data from the DNA methylation and gene expression profiles in TET2 mutated samples, with TET2 wild-type samples, and differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were identified. Following this, functional and pathway analyses were conducted and various key genes and pathways associated with DLBCL and TET2 mutations were identified.

 Materials and methods <sec> <title>Microarray data

The gene expression profile of GSE37365 was downloaded from the Gene Expression Omnibus (GEO) database (www.ncbi.nlm.nih.gov/geo/) and deposited by Asmar et al (5). It contained two subseries: the genome-wide gene expression dataset GSE37362 and the DNA methylation microarray dataset GSE37363. In GSE37362, there were 12 human DLBCL biopsy tissues with TET2 mutations and 19 TET2 wild-type DLBCL biopsy tissues, and the dataset was performed on the Illumina HumanMethylation450 BeadChip (GPL13534) platform (Illumina, San Diego, CA, USA). The gene expression dataset (GSE37363) contained 4 human DLBCL biopsy tissues with TET2 mutations and 5 TET2 wild-type DLBCL biopsy tissues and was detected with Affymetrix Human Gene 1.0 ST Array (GPL6244) platform (Affymetrix, Inc., Santa Clara, CA, USA).

 Data preprocessing

The original CEL gene microarray data were imported into R and the Affy package (http://bioconductor.org/packages/release/bioc/html/affy.html) was used for the background correction and robust multi array normalization (9). The probe-level symbols were converted into gene-level symbols. If multiple probes corresponded to one gene symbol, the average expression value of all probes was defined as the gene expression value. For the DNA methylation dataset, the IMA package (https://www.rforge.net/IMA) of R was used to normalize the methylation level in each methylation site. The methylation sites located in the X and Y chromosomes or SNPs were removed. P-value <0.05 was set as the cut-off criterion for the DNA-methylation matrix.

 Differential expression analysis

DEGs were identified by the limma package (http://bioconductor.org/packages/release/bioc/html/limma.html) (10) of R in TET2 mutated DLBCL samples compared with wild-type samples, according to the thresholds of Benjamini-Hochberg (BH) adjusted P-value <0.05 and |log2 (fold change)| >1. A total of 3 differentially methylated levels were identified. The first one was obtained by the P-value of a Student's t-test (P<0.05) in R (www.r-project.org) and the β value of the DNA methylation site in TET2 mutated samples compared with wild-type samples. The second one was determined by the gene regions, including the first exon area (EXON1), 1–200 bp upstream of the transcription start site (TSS200), 1–1500 bp upstream of the transcription start site (TSS1500), 3′-untranslated region (UTR) and 5′UTR. The third one was calculated by the locational data involving the CpG island, 1–200 bp upstream of the CpG island (N-CPGshore), 200–400 bp upstream of the CpG island (N-CPGshelf), 1–200 bp downstream of the CpG island (S-CPGshore) and 200–400 bp downstream of the CpG island (S-CPGshelf). Finally, with |∆β| >0.1 and P-value <0.05, the differentially methylated sites (DMSs) and DMRs were selected out in TET2 mutated samples and compared with wild-type samples.

The DMSs and DMRs were transformed to target genes termed differentially methylated genes (DMGs) using the annotation information, and the overlapping genes between the identified DMGs and DEGs were obtained.

 Functional enrichment analysis

The biological process (BP) of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for the aforementioned screened genes were identified using GOstats (http://www.bioconductor.org/packages/release/bioc/html/GOstats.html) (11) and GSEABase packages (http://www.bioconductor.org/packages/release/bioc/html/GSEABase.html) (12), respectively. The adjusted P-value <0.01 was selected as the cut-off criterion. The network of enriched pathways was constructed using the Pathview package (http://www.bioconductor.org/packages/release/bioc/html/pathview.html) (13).

 Results <sec> <title>Screening of DEGs and DMGs

The raw data were normalized in all samples and are presented in Fig. 1. A total of 198 DEGs (106 up- and 92 downregulated) were identified in TET2 mutated DLBCL samples compared with wild-type samples and the top 30 DEGs are listed in Table I. A total of 10493 DMSs, distributed over 3768 genes with |∆β| >0.1 and P-value <0.05 (group 1) were screened. A total of 1411 non-redundant genes (group 2) were selected from the DMRs (Fig. 2). As presented, the majority of DMRs were located in the gene body accounting for 26.5%, followed by TSS1500 with 20.2%. A further 1670 non-redundant genes (group 3) were identified by their location, of which 51.9% genes were located in CpG islands, 15.5% in N-shore, 12.1% in S-shore, 1.63% in N-shelf and 1.49% in S-shelf (Fig. 3). The results revealed that 602 DMGs were shared among the three groups. In addition, 12 overlapping genes, cryptochrome circadian clock (CRY)1, cytochrome B5 reductase (CYB5R) 2, doublecortin like kinase (DCLK) 2, fibronectin 1, glutathione peroxidase 3, short chain dehydrogenase/reductase family 42E, member 1 (SDR42E1), secreted frizzled related protein 2, spi-B transcription factor (SPIB), zinc finger protein (ZNF) interacting with K protein (ZIK) 1, ZNF134, ZNF256 and ZNF615 were identified in the 198 DEGs and 602 DMGs.

 Enriched GO terms and KEGG pathways of key target genes

It was demonstrated that 9 of the 12 overlapped genes were downregulated and hypermethylated, including CRY1, CYB5R2, DCLK2, SDR42E1, SPIB, ZIK1, ZNF134, ZNF256 and ZNF615. The heatmap of these 9 genes is presented in Fig. 4. Furthermore, using the cut-off criterion of P-value <0.01, 24 BPs of GO terms were obtained, including organic cyclic compound biosynthetic process (P=0.000131518), blue light signaling pathway (P=0.001212284) and cellular response to blue light (P=0.001212284; Table II). A significant KEGG pathway associated with circadian rhythm was obtained (Fig. 5).

 Discussion

TET2, one of the three proteins of the TET family, catalyzes the conversion of 5-methyl-cytosine to 5-hydroxymethyl-cytosine and promotes DNA demethylation (14). The present study, by using the data reported by Asmar et al (5), analyzed the DNA methylation and gene expression profiles in TET2 mutated samples compared with wild-type samples.

DNA methylation is one of the most common epigenetic modifications which may regulate gene expression in mammalian cells (15). Abnormal DNA methylation may lead to tumor formation and development (16,17). DLBCL represents the most frequently occurring B-cell malignancy with a high mortality rate, with a rapid progression and fatal outcome (18). It has previously been demonstrated that abnormal DNA methylation is associated with the pathogenesis and development of DLBCL. Shaknovich et al (19) reported that methylation of Src homology region 2 domain-containing phosphatase (SHP) 1 is frequently detected in plasma and tissues of DLBCL patients and therefore methylation of SHP1 may serve as a candidate biomarker for DLBCL. Additionally, TET mutations are typically observed in early hematopoietic progenitors (20) and it was hypothesized that they resulted in gene aberrant methylation in hematopoietic development (5). In the present study, 198 DEGs and 602 DMGs were identified and 9 genes CRY1, CYB5R2, DCLK2, SDR42E1, SPIB, ZIK1, ZNF134, ZNF256 and ZNF615 were revealed to be downregulated and hypermethylated in the TET2 mutated DLBCL samples compared with wild-type samples. The results demonstrated that CRY1, a gene associated with the circadian system (21), was significantly enriched in the circadian rhythm pathway. CRY1 and CRY2, interact with the period circadian clock (PER) proteins, and have a role in the negative limbs of the clock feedback loop (22). Circadian rhythms control a variety of metabolic and physiological functions, and abnormal circadian rhythms may be directly associated with the development of cancer (23). A previous study demonstrated that the downregulated clock gene PER2 in DLBCL patients was targeted by the transcription factor (CCAAT/enhancer-binding protein α and PER2 may act as a tumor suppressor gene (24). In addition, the epigenetic inactivation of aryl hydrocarbon receptor nuclear translocator-like protein (BMAL) 1, another circadian clock gene, has been revealed to contribute to the development of hematologic malignancies (25). Furthermore, CRY1 may modulate the ATR-mediated DNA damage checkpoint response by interacting with the timeless gene (26). The role of CRY1 as a key gene in the association between DLBCL and TET mutations has previously been investigated.

Furthermore, it was noted that ~50% of the aforementioned 9 genes encoded zinc finger proteins, including ZIK1, ZNF134, ZNF256 and ZNF615. Consistent with Asmar et al (5), ZIK1 was hypermethylated and downregulated in the results of the present study. Numerous studies have reported an association between zinc finger proteins and lymphoma or lymphatic disease. Bond et al (27) suggests that early hematopoietic zinc finger proteins are important in the regulation of human hematopoiesis and may be associated with the development of hematopoietic malignancies. Möröy et al (28) suggests that zinc finger protein growth factor independent 1 transcriptional repressor may maintain development and progression of lymphoid leukemia by suppressing the activation of the tumor suppressor p53. This suggests that different zinc finger proteins serve different roles in various types of cancer. Therefore, it was hypothesized that ZIK1, ZNF134, ZNF256 and ZNF615 may be involved in the development of DLBCL with TET2 mutations. Further validation of these genes is required.

In conclusion, the results of the present study suggested that the genes which were hypermethylated and downregulated in DLBCL biopsy tissues with TET2 mutations, exhibited important roles in the development of DLBCL. These genes, including CRY1, ZIK1, ZNF134, ZNF256 and ZNF615, may be affected by TET mutations and may act as potential biomarkers for the early diagnosis, therapy and prognosis of DLBCL. However, further experiments are required in order to validate the specific function of these genes in DLBCL.

 Acknowledgements

The present study was supported by the Municipal Science and Technology Commission of Tianjin (grant nos. 15ZLZLZF00440 and 16ZLZXZF00120), the Health Bureau Science and Technology Foundation of Tianjin (grant no. 2014KZ102) and the National Natural Science Foundation of China (grant nos. 81302082, 81272685, 31301151 and 81172355).

 References MontiSSavageKJKutokJLFeuerhakeFKurtinPMihmMWuBPasqualucciLNeubergDAguiarRCMolecular profiling of diffuse large B-cell lymphoma identifies robust subtypes including one characterized by host inflammatory responseBlood10518511861200510.1182/blood-2004-07-294715550490 CoiffierBDiffuse large cell lymphomaCurr Opin Oncol13325334200110.1097/00001622-200109000-0000311555708 TsangJCLoYMCirculating nucleic acids in plasma/serumPathology39197207200710.1080/0031302070123083117454749 TahilianiMKohKPShenYPastorWABandukwalaHBrudnoYAgarwalSIyerLMLiuDRAravindLRaoAConversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1Science324930935200910.1126/science.1170116193723912715015 AsmarFPunjVChristensenJPedersenMTPedersenANielsenABHotherCRalfkiaerUBrownPRalfkiaerEGenome-wide profiling identifies a DNA methylation signature that associates with TET2 mutations in diffuse large B-cell lymphomaHaematologica9819121920201310.3324/haematol.2013.088740238319203856967 SchumacherAKapranovPKaminskyZFlanaganJAssadzadehAYauPVirtanenCWinegardenNChengJGingerasTPetronisAMicroarray-based DNA methylation profiling: Technology and applicationsNucleic Acids Res34528542200610.1093/nar/gkj461164282481345696 SatoFTsuchiyaSTerasawaKTsujimotoGIntra-platform repeatability and inter-platform comparability of MicroRNA microarray technologyPLoS One4e5540200910.1371/journal.pone.0005540194367442677665 KhoshhaliMMahjubHSaidijamMPoorolajalJSoltanianARPredicting the survival time for diffuse large B-cell lymphoma using microarray dataJ Mol Genet Med6287292201210.4172/1747-0862.1000051231730133410377 IrizarryRAHobbsBCollinFBeazer-BarclayYDAntonellisKJScherfUSpeedTPExploration, normalization, and summaries of high density oligonucleotide array probe level dataBiostatistics4249264200310.1093/biostatistics/4.2.24912925520 da HuangWShermanBTLempickiRABioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene listsNucleic Acids Res37113200910.1093/nar/gkn92319033363 FalconSGentlemanRUsing GOstats to test gene lists for GO term associationBioinformatics23257258200710.1093/bioinformatics/btl56717098774 MericoDIsserlinRStuekerOEmiliABaderGDEnrichment map: A network-based method for gene-set enrichment visualization and interpretationPLoS One5e13984201010.1371/journal.pone.0013984210855932981572 LuoWBrouwerCPathview: An R/Bioconductor package for pathway-based data integration and visualizationBioinformatics2918301831201310.1093/bioinformatics/btt285237407503702256 SolaryEBernardOATefferiAFuksFVainchenkerWThe ten-eleven translocation-2 (TET2) gene in hematopoiesis and hematopoietic diseasesLeukemia28485496201410.1038/leu.2013.33724220273 OttoSPWalbotVDNA methylation in eukaryotes: Kinetics of demethylation and de novo methylation during the life cycleGenetics124429437199023073641203935 RobertsonKDDNA methylation, methyltransferases, and cancerOncogene2031393155200110.1038/sj.onc.120434111420731 JonesPABaylinSBThe fundamental role of epigenetic events in cancerNat Rev Genet3415428200212042769 PasqualucciLTrifonovVFabbriGMaJRossiDChiarenzaAWellsVAGrunnAMessinaMElliotOAnalysis of the coding genome of diffuse large B-cell lymphomaNat Genet43830837201110.1038/ng.892218045503297422 ShaknovichRGengHJohnsonNATsikitasLCerchiettiLGreallyJMGascoyneRDElementoOMelnickADNA methylation signatures define molecular subtypes of diffuse large B-cell lymphomaBlood116e81e89201010.1182/blood-2010-05-285320206108142993635 QuivoronCCouronnéLDella ValleVLopezCKPloIWagner-BallonODo CruzeiroMDelhommeauFArnulfBSternMHTET2 inactivation results in pleiotropic hematopoietic abnormalities in mouse and is a recurrent event during human lymphomagenesisCancer Cell202538201110.1016/j.ccr.2011.08.00621723201 KoCHTakahashiJSMolecular components of the mammalian circadian clockHum Mol Genet 15 Spec No2R271R277200610.1093/hmg/ddl207 KumeKZylkaMJSriramSShearmanLPWeaverDRJinXMaywoodESHastingsMHReppertSMmCRY1 and mCRY2 are essential components of the negative limb of the circadian clock feedback loopCell98193205199910.1016/S0092-8674(00)81014-410428031 SaharSSassone-CorsiPMetabolism and cancer: The circadian clock connectionNat Rev Cancer9886896200910.1038/nrc274719935677 ThoennissenNHThoennissenGBAbbassiSNabavi-NouisSSauerTDoanNBGerySMüller-TidowCSaidJWKoefflerHPTranscription factor CCAAT/enhancer-binding protein alpha and critical circadian clock downstream target gene PER2 are highly deregulated in diffuse large B-cell lymphomaLeuk Lymphoma5315771585201210.3109/10428194.2012.65879222260161 TaniguchiHFernándezAFSetiénFRoperoSBallestarEVillanuevaAYamamotoHImaiKShinomuraYEstellerMEpigenetic inactivation of the circadian clock gene BMAL1 in hematologic malignanciesCancer Res6984478454200910.1158/0008-5472.CAN-09-055119861541 KangTHLeemSHModulation of ATR-mediated DNA damage checkpoint response by cryptochrome 1Nucleic Acids Res4244274434201410.1093/nar/gku094244891203985666 BondHMMesuracaMCarboneEBonelliPAgostiVAmodioNDe RosaGDi NicolaMGianniAMMooreMAEarly hematopoietic zinc finger protein (EHZF), the human homolog to mouse Evi3, is highly expressed in primitive human hematopoietic cellsBlood10320622070200410.1182/blood-2003-07-238814630787 MöröyTThe zinc finger protein Gfi1 maintains development and progression of lymphoid leukemia by blocking the activation of the tumor suppressor p53Experimental Hematol42SupplS7201410.1016/j.exphem.2014.07.024

Boxplots for microarrays (A) prior to and (B) following application of normalization. Red boxes represent TET mutated diffuse large B-cell lymphoma samples. Blue boxes represent TET2 wild-type samples. TET, Ten-Eleven Translocation hydroxylases.

The distributions of differentially methylated genes based on the genetic regions. UTR, untranslated region; TSS, transcription start site.

The distributions of differentially methylated genes based on the locations of CpG islands.

Hierarchical clustering heatmap of CRY1, CYB5R2, DCLK2, SDR42E1, SPIB, ZIK1, ZNF134, ZNF256 and ZNF615. Each column corresponds to a single microarray and each row represents expression profile of a single gene. Red and blue represent the low and high values of the gene expression, respectively. Case-1, 2, 3, 4 represent human DLBCL biopsy tissue samples with TET mutations and control-1, 2, 3, 4, 5, represent TET2 wild-type DLBCL biopsy tissue samples. DLBCL, diffuse large B-cell lymphoma; TET, Ten-Eleven Translocation hydroxylases; CRY1 cryptochrome circadian clock 1; CYB5R2, cytochrome B5 reductase 2; DCLK2, doublecortin like kinase; SDR42E1, short chain dehydrogenase/reductase family 42E, member 1; SPIB, spi-B transcription factor; ZNF, zinc finger protein; ZIK1, ZNF interacting with K protein.

Pathway of circadian rhythm. The red gene is the significant gene involved in the pathway. CRY1 cryptochrome circadian clock 1; KEGG, kyoto encyclopedia of genes and genomes; CREB1, CAMP responsive element binding protein 1; MAPK, mitogen activated protein kinase; Fbx 13, F-box 13; Ror, RAR related orphan receptor; PER, period circadian clock; Cklε, casein kinase 1 ε; Cklδ, casein kinase 1δ; SCF, Skp1-Cul1-F-box protein; β-TrCP, β-transducin repeat containing E3 ubiquitin protein ligase.

The top 30 DEGs in TET2 mutated DLBCL samples compared with wild-type samples.

Gene LogFC P-value
ZIK1 −1.00204 5.43×10−5
SLC9A7 −2.12435 8.30×10−5
FAM129C −1.78254 0.000152
VCAN 1.888421 0.000272
ZNF256 −2.00847 0.00029
ZNF542P −1.48907 0.000365
RNU2-6P −2.03556 0.00037
KIAA0226L −3.08868 0.000435
VDR 1.206615 0.000439
IFT57 −1.36341 0.00056
APOBEC3B 2.477052 0.000594
ZNF667-AS1 −2.00769 0.000808
SERPINB9 1.405557 0.001202
SDR42E1 −1.28485 0.001228
BCL2 −1.23809 0.001266
KCNN3 1.915499 0.001271
ZNF350 −1.70697 0.001372
PRRG4 1.080011 0.001556
SLC7A11 2.276444 0.001604
CEP68 −1.00705 0.001643
SPP1 3.186619 0.001661
ZNF382 −1.46886 0.00197
EMP1 1.098811 0.001994
PLA2G2A 1.148008 0.002013
GPR183 −2.03511 0.002126
PDGFRA 2.247303 0.002385
CHPF 1.078177 0.002536
CRY1 −1.08059 0.00269
LZTFL1 −1.77782 0.002797
AKT3 −1.13853 0.002883

DEGs, differentially expressed genes; TET, Ten-Eleven Translocation; DLBCL, diffuse large B-cell lymphoma; LogFC, log (fold change).

Significant 24 GO terms for screened genes.

GO term Description P-value Genes
GO:1901362 Organic cyclic compound biosynthetic process 0.000131518 CRY1, SPIB, ZNF134, ZNF256, CYB5R2, SDR42E1, ZIK1, ZNF615
GO:0009785 Blue light signaling pathway 0.001212284 CRY1
GO:0071483 Cellular response to blue light 0.001212284 CRY1
GO:1901576 Organic substance biosynthetic process 0.001486849 CRY1, SPIB, ZNF134, ZNF256, CYB5R2, SDR42E1, ZIK1, ZNF615
GO:0009058 Biosynthetic process 0.00165506 CRY1, SPIB, ZNF134, ZNF256, CYB5R2, SDR42E1, ZIK1, ZNF615
GO:0009637 Response to blue light 0.001817937 CRY1
GO:1901360 Organic cyclic compound metabolic process 0.00276577 CRY1, SPIB, ZNF134, ZNF256, CYB5R2, SDR42E1, ZIK1, ZNF615
GO:0006355 Regulation of transcription, DNA-dependent 0.002884379 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:2001141 Regulation of RNA biosynthetic process 0.002991969 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0006694 Steroid biosynthetic process 0.003173823 CYB5R2, SDR42E1
GO:0051252 Regulation of RNA metabolic process 0.003375671 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:2000112 Regulation of cellular macromolecule biosynthetic process 0.004585358 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0006351 Transcription, DNA-dependent 0.004971651 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0010556 Regulation of macromolecule biosynthetic process 0.005245834 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0032774 RNA biosynthetic process 0.005738022 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0031326 Rregulation of cellular biosynthetic process 0.006573337 CRY1, SPIB, ZNF134, ZNF256, ZIK1,
GO:0009889 Regulation of biosynthetic process 0.006882303 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0010468 Regulation of gene expression 0.006955993 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0019219 Regulation of nucleobase-containing compound metabolic process 0.007355736 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0018298 Protein-chromophore linkage 0.007856527 CRY1
GO:0051171 Regulation of nitrogen compound metabolic process 0.008353762 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0034654 Nucleobase-containing compound biosynthetic process 0.008928906 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0018130 Heterocycle biosynthetic process 0.009820777 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615
GO:0019438 Aromatic compound biosynthetic process 0.009834588 CRY1, SPIB, ZNF134, ZNF256, ZIK1, ZNF615

GO, Gene Ontology.