REGγ is associated with multiple oncogenic pathways in human cancers

Background Recent studies suggest a role of the proteasome activator, REGγ, in cancer progression. Since there are limited numbers of known REGγ targets, it is not known which cancers and pathways are associated with REGγ. Methods REGγ protein expressions in four different cancers were investigated by immunohistochemistry (IHC) analysis. Following NCBI Gene Expression Omnibus (GEO) database search, microarray platform validation, differential expressions of REGγ in corresponding cancers were statistically analyzed. Genes highly correlated with REGγ were defined based on Pearson's correlation coefficient. Functional links were estimated by Ingenuity Core analysis. Finally, validation was performed by RT-PCR analysis in established cancer cell lines and IHC in human colon cancer tissues Results Here, we demonstrate overexpression of REGγ in four different cancer types by micro-tissue array analysis. Using meta-analysis of publicly available microarray databases and biological studies, we verified elevated REGγ gene expression in the four types of cancers and identified genes significantly correlated with REGγ expression, including genes in p53, Myc pathways, and multiple other cancer-related pathways. The predicted correlations were largely consistent with quantitative RT-PCR analysis. Conclusions This study provides us novel insights in REGγ gene expression profiles and its link to multiple cancer-related pathways in cancers. Our results indicate potentially important pathogenic roles of REGγ in multiple cancer types and implicate REGγ as a putative cancer marker.


Background
REGγ, also known as PA28gamma, 11Sgamma, or PSME3, was first identified as Ki antigen, a nuclear protein targeted by autoantibodies found in sera of patients with systemic lupus erythematosus [1]. It is a member of the 11S family of proteasomal activators that have the ability to stimulate the proteolytic activity of the 20S core proteasome independent of ubiquitination and ATP [1]. Accumulating evidence suggests REGγ is involved in cancer progression [2]. REGγ has been reported to be overexpressed in colorectal cancer [3] and thyroid cancer [2], and is involved in cancer development [2,[4][5][6]. It is unknown, however, whether REGγ is involved in additional cancers. REGγ is known to degrade both oncogenic and tumor suppressing proteins such as SRC-3, HCV core protein, PTTG1, p21, p16, p19, and p53. In this study we try to understand expression profiles of REGγ in multiple cancer types and correlations of REGγ with known cancer or cancer related pathways.
Microarray assays have been widely adopted in cancer marker exploration and expression profiling of tumor genes [3,4]. Microarray studies have contributed valuable information to our understanding of cancer by identifying biomarkers and enabling classification of tumor subtypes [5][6][7][8].
In this study, we focused on thyroid cancer, colon cancer, liver cancer and lung cancer since the first two cancers were reported with over-expression of REGγ [3,9] and the other two are among the list of the most malicious cancers. We analyzed REGγ expression in cancer tissue arrays by using publicly available microarray data from NCBI GEO database. We acquired datasets and integrated the analyzed results across different datasets and cancer types to characterize a general REGγ expression pattern in four different cancer types by comparing human cancer versus normal tissues. We set clear criteria along with quality controls for dataset screening and normalization, which allowed us to carry out comprehensive datasetbased meta-analysis across differing cancers. A set of genes highly correlated with REGγ expression were identified and validated by RT-PCR to identify putative functional interactions associated with REGγ.

Methods
Cell types and cell culture A549, HepG2, and HCT116 cells were purchased from ATCC and maintained at Cell Culture Core at the Department of Cell Biology, BCM. The human thyroid carcinoma cell line ARO was kindly provided by Dr. Adel El-Naggar at the University of Texas M.D. Anderson Cancer Center. The ARO cell line was authenticated at Genotyping Center of John's Hopkins University. The shN and shR stable cell lines were generated in ARO, A549, and HCT116 by introducing retroviral shRNA vectors specific for REGγ or a control vector from OriGene (Rockville, MD). ARO cells were cultured in 1640 supplemented with 10% fetal bovine growth serum (GIBCO). All other cells were cultured under standard conditions described by the ATCC.

Immunohistochemical assay
IHC analysis was performed to analyze REGγ expression of protein level in several human cancers including lung, colon, thyroid and liver cancer. Sections were deparaffinised and rehydrated. The slides were then heated in a 100°C water bath for 30 minutes in a 0.01 M citrate buffer solution at pH 6.0, and cooled to room temperature. After quenching the endogenous peroxidase activity with 0.3% H2O2 (in absolute methanol) for 10 minutes, the sections were treated for 10 minutes at room temperature with the serum albumin (HISTOSTAIN-PLUS DAB kit) to block non-specific staining. Duplicate sections were incubated overnight in 4°C with the primary specific antibodies. Slides were then incubated for 10 minutes with biotinylated anti-rab-bit IgG (DAB kit) for REGγ recognition. The sections were incubated with the HRP for 10 minutes. Finally, the sections were counterstained with Mayer's haematoxylin.

Preliminary datasets collection
Microarray expression profiles were obtained from Gene Expression Ominibus of National Center of Biotechnology Institute [10]. All datasets in this study were published within the past 5 years (2004-2009) and following the Minimum Information about a Microarray Experiment (MIAME) guidelines, including 49 datasets, with 3,832 samples containing 16, 15, 11, and 7 datasets from colon, liver, lung and thyroid cancer, respectively (Additional file 1: Table S2). The following preliminary datasets were retrieved: a) primary tumors, carcinoma and adenoma along with normal controls in each tissue: primary colon cancer samples including early onset colorectal carcinoma, colon tumor and adenoma; primary hepatocelluar carcinoma (HCC); lung cancer including non-small cell lung cancer, adenocarcinoma, and squamous cell carcinoma; thyroid cancer samples including papillary thyroid carcinomas (PTC), anaplastic thyroid carcinoma (ATC), follicular carcinomas (FC) and follicular adenomas. b) noncancer diseases originated from colon, liver, lung and thyroid tissues, including inflammatory bowel disease (IBD), Crohn's disease (CD), ulcerative colitis (UC), HCV cirrhosis, HCV-induced dysplasia, pneumonia, and follicular goiter. c) different stages of cancers with a stage 0 tissue or healthy or distant adjacent tissues as control.
The following datasets/samples were excluded: 1) datasets with no contorl tissue; 2) datasets without REGγ probe/probe set included in platform; 3) datasets without corresponding publication; 4) datasets with samples collected in time courses; 5) datasets without gene symbol annotation for the probes by the Human Gene Nomenclature (HGNC) guidlines; 6) datasets without REGγ data in the microarray platform. The preliminary sample screening yielded 23 datasets (n = 1,070) for differentially expressed gene analysis (Additional file 2: Table S3). Of these output, there were 15 cancer datasets (n = 614), 2 non-cancer diseases datasets of (n = 86), and 6 datasets (n = 399) containing both cancer (n = 260) and disease samples (n = 139).

Dataset based expression analysis
Microarray datasets (described above) were analyzed by GEOquery and Limma packages in R http://www.r-project.org/ as described previously [11][12][13]. First, all raw data were downloaded from GEO and mono-channel data were normalized using MAS5.0. Samples in each dataset were grouped into three classes, namely cancer, non-cancer disease and normal samples. The log2 ratio values of disease group versus normal group were calculated based on the normalized data. For all two-channel datasets, log2 transformed expression ratios were calculated and used in all subsequent analyses. Two-sample paired t-test [14] were carried out between cancer vs. non-cancer diseases and cancer vs. normal following statistic analysis as described [15,16]. Internal quality controls were set up for each dataset by validating the statistical significance of specific genes with what was reported in relevant publications. Two-sample comparisons were statistically analyzed for all 21 cancer datasets containing 874 cancer samples and 625 paired normal samples. An additional two-sample comparison was performed with the 8 non-cancer disease datasets, including 196 non-cancer disease samples and 174 normal control samples (Additional file 2: Table S3).

Correlation analysis
The Pearson's Correlation Coefficient (PCC) was used as a measure of correlation between REGγ and its potentially related genes based on 13 datasets, 4 from liver (n = 299) and 3 from each of lung (n = 164), colon (n = 77), and thyroid (n = 126) respectively. Pearson Correlation analysis was conducted using R [15,17] on datasets with significant overexpression of REGγ. PCC of REGγ with each gene in each dataset was calculated. Genes whose expression correlated with REGγ in each dataset were ranked based on their p-value. In order to produce at least 600 candidates in each datasets for subsequently selection, we used a cutoff of 0.001. The top 20%, 15%, 10%, and 5% genes were selected from thyroid, colon, liver, and, lung cancer datasets respectively. All subsequent selections and analyses were based on these genes referred to as REGγ correlated genes.
Genes were selected from all REGγ correlated genes based on cancer type except for the initial pilot testing (a PCC cutoff as ± 0.6 [18]). Our criteria were that each gene was present in at least 2 datasets, according to binomial distribution (p < 0.05), in one cancer type and the cutoff of PCC in one cancer type was set to ± 0.6. Genes that fulfill these criteria were considered as highly-correlated with REGγ and used for downstream pathway analysis.

Pathway and network analysis
Genes highly correlated with REGγ were analyzed with the IPA (Ingenuity Pathway Analysis) system: http:// www.ingenuity.com. With core analysis, all qRT-PCR validated REGγ-correlated genes were mapped and then analyzed using Ingenuity Knowledge Base (genes only) to yield bio-function pathway annotation and networks showing direct and indirect relationships between genes and molecules.
To calculate the composition of REGγ-correlated genes/pathways in cancers, results from Ingenuity pathway analysis were grouped into three clusters: cancer pathways, cancer related pathways, and other pathways. These pathway clusters were grouped based on the following characterization: 1) cancer pathways included bio-function of cancer, tumor or tumorigenesis, neoplasia, carcinoma or adenocarcinoma, lymphoma and sarcoma; 2) cancer related pathways included a) pathways related to cell cycle with following bio-function category: mitosis or mitotic, G2/M/S phase, cell division, check point, and arresting; b) cell growth pathways involved in: survival, growth and proliferation; c) cell death pathways with bio-function of apoptosis and death. 3) Other pathways: all the rest of the pathways not included in cancer or cancer-related pathway clusters.
Genes highly correlated with REGγ were also searched against the KEGG pathways database http://www.genome.jp/kegg/tool/colorpathway.html to highlight and augment the published graphical pathways analyzed by Ingenuity. Protein-protein interaction network analysis was performed by checking REGγ highly-correlated genes in the STRING database http://string-db.org) [19]. To make the network concise, genes with connections equal or greater than 3 were selected.

PCR validation
Confirmatory qRT-PCR was performed on randomly selected REGγ correlated genes. Fifteen genes were selected from the REGγ-correlated genes and an additional fifteen genes highly correlated with REGγ expression were selected for qRT-PCR. RT-PCR experiments were carried out in cells originated from colon, liver, lung and thyroid cancer.
RNA interference was performed in HepG2 cells to knock down REG. HepG2 was seeded in 6-well plate at 60% confluence overnight and was transfected with 10 nM siRNA along with lipofection 2000. Cells were harvested 72 hours later for RNA extraction and qRT-PCR analysis.

Statistic analysis
Weighted student t-test for two-sample with unequal variance was used to calculate statistics and p values for IHC tissue array [16]. Two-tailed student's t-test was used in microarray expression analysis and Fisher's Z transformation was used to adjust p value. A p-value of less than 0.05 was defined as "significant" for all statistic analysis involved in expression analysis. Datasets in which REGγ is highly significant (p < 0.001) [18] were selected for subsequent correlation analysis. In correlation analysis, Pearson's correlation coefficient was set with a cutoff PCC ± 0.6 [18] and binomial coefficient was used based on datasets number in each cancer type to selection REGγ highly related genes. Pathways with a p-value less than 0.01 were chosen to be studied in Ingenuity core analysis.

REGg protein is highly expressed in multiple cancers
To understand whether REGγ is a tumor-associated protein, we examined REGγ expression levels in multiple human carcinomas. IHC experiment was performed using tissue arrays containing 92 cases of primary lung cancer, 48 colon cancers, 49 thyroid cancers, and 206 liver cancer samples along with corresponding normal tissues, all arranged in duplicates ( Table 1, Additional file 3: Table S1). The expression of REGγ in cancer samples was scored double-blindly by comparing with normal tissues or adjacent non-cancer tissues which have no positive staining or low levels of REGγ staining ( Figure 1 and Table 1). The scored REGγ expression is consistent for most of the duplicate samples and a representative scored result was shown in Additional file 3: Table S1. The overall rate of REGγ overexpression (a combination of ++ and +++ staining) in different carcinoma is higher than 50%. We observed a statistically significant increase in the number of late-stage cancers with the highest REGγ expression (+++), such as in stage III of adenocarcinoma and squamous cell carcinoma (Additional file 3: Table  S1). Our results provide the first evidence for an association of REGγ with primary human lung carcinoma and liver cancer, substantiating previous observations that REGγ is increased in colon and thyroid cancers.

Integrated analysis of microarray datasets revealed overexpression of REGg in selective cancers
Overexpression of REGγ protein in four different human cancers prompted us to investigate whether elevation of REGγ is regulated at the mRNA level. We searched GEO database by keywords and identified 49 datasets (Additional file 1: Table S2), of which 23 were qualified for expression analysis in this study (Figure 2A). Significantly higher REGγ expression (p < 0.05) was observed in 67% of cancer datasets (14 of 21, n = 597) when compared with normal tissues ( Figure 2B). Consistently, our comparative analysis of control vs. non-cancer diseases [20][21][22], revealed that most of the non-cancer datasets had no significant differences in REGγ expression (66%, n = 117). On the contrary, only small percent of cancer datasets (7 of 21, n = 228) had no significant elevation in REGγ levels ( Figure 2B), indicating potential association of REGγ in the development of these cancers. Cancer type based analysis indicated an increase of REGγ in 60-83% of cancer datasets ( Figure 2C), concordant with our IHC studies (Figure 1).
A detailed analyses of pathologically classified, stagespecific cancers and non-cancer diseases were executed using dataset GSE6764 [23], GSE4183 [22], GSE6339 [21] and GSE7670 [20], which originate from liver, colon, thyroid, and lung respectively and disclosed detailed cancer stage information (Additional file 2: Table S3). The representative REGγ expression patterns in the four cancers, non-cancer diseases, normal controls, and some staged cancer samples are illustrated in Figure 3. Statistically significant elevation of REGγ gene expression in cancers ranged from 1.25 to 2.43 fold-change (overexpression cut-off: fold change > 1), consistent with our IHC result in corresponding cancer tissue arrays ( Figure 1). In liver cancer samples, REGγ appeared considerable up-regulation, consistent with the original publication where potential cancer biomarkers were linked with stepwise carcinogenic process [23]. The results of stage-associated increases of REGγ in advanced liver cancers in Figure 3A is consistent with our observation of higher REGγ staining in advanced cancer types. For most of the non-cancer datasets, the p-values of disease classes, such as IBD & CD, showed no significant changes in REGγ expression. In conclusion, expression of REGγ is significantly increased in multiple human cancer types and likely involved in late-stage cancers.   Correlation analysis links REGg to large numbers of cancer related genes/pathways To explore potential mechanisms of REGγ in cancer development, we further investigated genes whose expression is highly related to REGγ expression in the four cancer types profiled. A statistical meta-analysis based on Pearson correlation coefficient (PCC) was conducted on the defined (REGγ differential expressed) datasets (Additional file 4: Table S4). The correlation between REGγ and every other gene in these datasets were calculated and evaluated statistically. With the assumption that a high absolute PCC value would reflect a potentially close relation to REGγ functionally, only genes bearing high PCC scores were selected for subsequent studies.
To estimate that the approach we used in our analysis could indeed generate meaningful results, we first set up a PCC cutoff value +/-0.6 in at least one dataset for a pilot test. Since previous study [24,25] has demonstrated REGγ mediated regulation of p53, we examined if p53 targets can be identified among REGγ highly correlated genes. A total of 29 published genes in p53 signaling pathway were identified as significantly correlated with REGγ (Additional file 5: Table S5), indicating that our normalized datasets contain valuable information required for further analysis.
By more stringent PCC-value based criteria, we screened genes highly-correlated with REGγ and identified a total of 588 genes, with 521 positively correlated, 99 negatively correlated, and 31 being both negatively and positively correlated (Additional file 6: Table S6). Among these genes strongly correlated with REGγ, 467 were identified from colon cancer, while 75, 53, and 25 genes were from lung, thyroid and liver cancer respectively. Interestingly, multiple cancers shared significant amount of these genes. Based on all calculated results, there were 32% genes in two cancers, 43% genes in three cancers and 21% genes in four cancers simultaneously (Additional file 7: Figure  S1A). The PCC between REGγ and all other genes in each dataset are shown in Figure 4A. The range of the PCC plot reflects more positive-correlation points than negative ones, suggesting that more genes are positively correlated with REGγ expression. To understand functional diversities of the genes computationally correlated with REGγ, we performed Ingenuity pathway analysis of the 588 genes. Our analysis displayed that all mapped genes were functionally annotated into 500 pathways in which 207 were statistically significant (P <0.01). Among the 207 pathways analyzed, 20 cancer pathways (10%), 102 cancer related pathways (49%), and 85 other pathways (41%) were classified ( Figure 4B). The top 15 pathways based on statistic significance (p-values) are shown in Additional file 7: S1B. Cancer related pathways were composed of 86 pathways related to cell cycle regulation, 9 pathways in apoptosis, and 7 pathways in cell growth (Additional file 6: Table S6). Due to the presence of subset of genes in multiple pathways, such as Myc (in 37 function annotation pathways), HSP90AB1 (in 11 pathways), ILF2 and ODC1 (in 6 pathways respectively), the number of genes in cancer, cancer-related and other pathways were 163, 168, and 218 respectively with overlaps indicated in Figure 4C. Detailed information of pathway analysis is included in Additional file 6: Table  S6. Based on Ingenuity analysis of cancers with REGγ overexpression, our results indicate that over 50% of REGγ highly-correlated genes/pathways are cancer or cancer-related.
We also validated our pathway analysis of REGγ-correlated genes by applying all 588 REGγ highly correlated genes to KEGG pathway annotation. The results were consistent with Ingenuity analysis whereby cell cycle and cancer pathways were ranked among the top (Additional file 8: Table S7). Based on these annotation analyses, we discovered that REGγ is linked to large numbers of cancer To illustrate potential links between the gene products among the REGγ correlated genes, we conducted further analysis of protein-protein interaction (PPI) network using STRING [19], which is an online database of known and predicted protein interactions. This generated network (Additional file 9: Figure S2) integrated information from experimental repositories, computational prediction and published collections, and showed their interaction with default parameters. PPI network revealed potential interactions among five clusters of REGγ correlated gene products, including those in metabolic pathways, proteasome pathways, cell cycle related pathways, DNA repair pathways, and tRNA biosynthesis pathways. These results provide additional information for future study of cellular function of REGγ as well as its regulation.

Confirmatory analysis of REGg correlated genes from bioinformatic analysis
Our computational analysis indicated strong correlation of REGγ to genes regulated by p53 and in cancer related pathways. To validate our bioinformatics-based predictions, we selected 30 genes for expression analysis using specific cancer cell lines. In addition to genes associated with p53 pathways, we selected two representative genes from each of the top cancer/cancer related pathways, metabolic pathways as well as those appeared in KEGG and Ingenuity network analysis (Table 2). We made use of stable cell lines constitutively expressing a control shRNA (shN) or a REGγ specific shRNA (shR). Three pairs of shRNA expressing cell lines were originated from lung, colon, and thyroid (A549, HCT116, and ARO, Additional file 10: Figure S3). The REGγ knockdown in HepG2 liver cancer cell lines was generated by introducing synthetic siRNA against REGγ [26]. The significant differences of REGγ expression between control (shN) and REGγ knockdown in each pair of the cell lines allowed us to examine the relative levels of genes predicted to be highly correlated to REGγ (Additional file 10: Figure S3). Real time RT-PCR was performed with specific primers for each gene in multiple pairs of cell lines or in a specific pair of cells dependent upon cancer-type specific bioinformatics data (Additional file 11: Table S8). The relative expression ratios of each specific gene in shN and shR cancer cell lines were calculated from average of three independent experiments and the data with statistically significant changes (p <0.05) between shN and shR cells were plotted ( Figure 5A and 5B). The expression of most genes (66.7%, n = 20) in at least one cell line with or without REGγ knockdown was consistent with predicted correlation to REGγ levels ( Table 2).
All genes validated by RT-PCR were applied into Ingenuity system for core analysis. Network information showing the link among these REGγ correlated genes was displayed in Figure 6A. This analysis placed Myc as the hub of the interaction network and prompted us to perform further analysis on its biological significance. Interestingly, analysis of 11 colon cancer samples suggested significant positive correlation between cMyc and REGγ ( Figure 6B and 6C). The p53 target, PTEN, was also observed in this network analyses, reinforcing the close correlation among REGγ, p53 pathway, and other cancer related pathways.

Discussion
REGγ-proteasome system represents an emerging pathway recently recognized to be involved in cancer development. This study provides further links between REGγ and multiple cancer related pathways by a combination of bioinfomatic analysis and molecular biological approach.
To our knowledge, this is the first computational study so far in REGγ association with multiple cancers. We are also the first to demonstrate high expression of REGγ in lung and liver cancers despite that overexpression of REGγ in thyroid and colon cancer were reported [6,27]. Tissue array analyses of four different human cancers, including lung, colon, thyroid, and liver cancers, revealed significant increase of REGγ protein in over 50% of these cancer samples. Bioinformatic analysis of human microarray gene expression profiles indicates that REGγ gene expression is  Figure 5 Results of confirmatory qRT-PCR on selected genes. (A) Polar plot of qRT-PCR validation results of selected genes highly correlated with REGγ. Following qRT-PCR analysis of specific genes in REGγ knockdown or control cells, data averaged from three independent experiments were converted into relative fold changes. Fold change values greater than one was shown as positive correlation (brown) and values less than one represented negative correlation (blue). Genes were arranged in different theta (θ) and radius represent the qRT-PCR fold change values from 0 to 3. Only data consistent with prediction were shown (see Table 2). (B) Representative qPCR validation experiments. Quantitative RT-PCR was performed in paired cancer cell lines (shN and shR) with differential levels of REGγ expression. Subsets of qRT-PCR results showed the actual expression differences of a specific gene in these cell lines (results were average from three independent experiments). The relative PCR values ( Table 2) in shN (usually normalized as 1) were divided by the values in shR cells to yield fold changes as shown in Figure 5A. A value greater than 1 indicates a positive correlation with REGγ. also increased in most of these human cancers, providing new evidence that REGγ-proteasome pathway may be involved in the development of multiple cancers.
Computational analysis of datasets from thyroid cancer with thyroid non-cancer disease and liver cancer with HCC clinical stage information indicated a potential correlation of gradual increase of REGγ level with cancer stages ( Figure 3A). Although the sample size and numbers are relatively small, the results suggest a potential of REGγ as a prognostic cancer marker and hinted some molecular mechanisms linking REGγ to development of cancers toward later stages or malignancy.
Our meta-analysis disclosed significant correlation between REGγ and many genes in cancer and cancer related pathways from ingenuity analysis, including colorectal cancer, lung carcinoma, sarcoma, lymphoma, tumorigenesis, cell division and apoptosis related pathways etc. Importantly, genes downstream of the previously identified REGγ regulated proteins [24,25], p53, was found highly correlated with REGγ expression. Despite that p53 mutation in different cancer may complicate the correlation status of its downstream target genes with REGγ, the overall high correlation values strongly support the previous finding that REGγ-mediated regulation of p53 may play an important role in cancer development. Annotation analysis indicated significant correlation of REGγ with many different proteasome components, suggesting that REGγ may be elevated and function together with other proteasome complexes. Recent studies have indicated important roles of the ubiquitin-mediated protein degradation pathway in cancer and application of proteasome inhibitors as a promising anti-cancer therapeutic approach [28]. An additional interesting finding is that numerous REGγ-correlated genes are involved in metabolism, particularly in energy metabolism (Additional file 8: Table S7). The link between cell metabolism and cancers has been well documented [29]. The result concurs a recent notion that cancer cell metabolisms are controlled by oncogenes and tumor suppressor genes [30].
The mathematical approaches of bioinformatics used in this study are quite standard [18]. To minimize the chance of acquiring false-positive results and ensure that most of the strong candidate genes are selected, we set rigorous criteria for all studies performed. The significance of our computational analysis has been underlined by laboratory validation experiments in 4 pairs of cancer cell lines in which differential expression of REGγ was created in the same background to facilitate correlation studies. Results from quantitative analysis of selected genes were largely consistent with predicted correlations, suggesting powerful combination of bioinformatics and molecular biological studies in disclosing potentially novel functions of REGγ-proteasome in cancer progression. It is not unlikely that REGγ could serve as a cancer marker, particularly for cancers with aggressive behavior.
Given that REGγ mainly functions as a proteasome activator to induce protein degradation, the biological links between REGγ and its correlated genes may reflect a result of direct or indirect regulation on transcription. Ingenuity analysis of validated gene network led our attention to the correlation between REGγ and Myc gene. Previous research documented overexpression of both genes in colorectal cancers [31][32][33]. Coincident with the largest amount of datasets and highest REGγ differential expression from colon cancers, the positive correlation between REGγ and Myc was validated specifically in HTC116-shN and -shR cells. Since Myc functions as a transcription factor, we searched REGγ promoter and found numerous Myc binding sites within 1.5 KB upstream REGγ transcriptional initiation site (data not shown). Yet we could not exclude the possibility that REGγ may target a negative regulator of Myc for degradation. Further experiments will be carried out to understand the molecular detail of these hypotheses. It is likely that elevated expression of Myc in certain cancer cells is one of the potential mechanisms contributing to higher expression of REGγ.

Conclusions
This study provides REGγ expression profiles based on computational analysis of published microarray datasets and laboratory experiments on cancer samples. Data analysis links REGγ to multiple cancer-related pathways.
Our results indicate potentially important roles of REGγ in multiple cancer types and implicate REGγ as a putative cancer marker.

Funding
This work was supported by National Institutes of Health (1R01CA131914) and Norman Hackerman