Skip to main content

Expression analysis of genes and pathways associated with liver metastases of the uveal melanoma



Uveal melanoma is an aggressive cancer which has a high percentage metastasizing to the liver, with a worse prognosis. Identification of patients at high risk of metastases may provide information for early detection of metastases and treatment.


Expression profiling of ocular tumor tissues from 46 liver metastatic uveal melanoma samples and 45 non-metastatic uveal melanoma samples were got from GEO database. Bioinformatic analyses such as the Gene Oncology and Kyoto Encyclopedia of Genes and Genomes were used to identify genes and pathways specifically associated with liver metastases of the uveal melanoma.


A total of 1138 probes were differentially expressed in two group samples. All differential gene interactions in the Signal-Net were analyzed. Of them, 768 probes were up-regulated and 370 down-regulated. They mainly participated in 125 GO terms and 16 pathways. Of the genes differentially expressed between two group cancers, HTR2B, CHL1, the ZNF family, YWHAZ and FYN were the most significantly altered.


Bioinformatics may help excavate and analyze large amounts of data in microarrays by means of rigorous experimental planning, scientific statistical analysis and collection of complete data about liver metastases of uveal melanoma patients. In the present study, a novel differential gene expression pattern was constructed and advanced study will provide new targets for diagnosis and mechanism of uveal melanoma liver metastases.

Peer Review reports


Uveal melanoma is the most common primary intracellular tumor in adults with an estimated 5-year survival rate of 50%-70% [1]. About 50% of patients develop metastases within a median of 36 months, mostly to the liver, with a median survival of 6 months after metastases [2]. It seems important to identify high-risk patients at the time of the initial diagnosis for early detection and treatment of metastatic disease or for the administration of adjuvant therapy. Several clinical and histopathological features have been correlated with survival, including patient age (>60), anterior location of the tumor, tumor cell histology, largest diameter of the tumor, mitotic activity, and chromosome 3 monosomy. The most frequent chromosomal imbalances in uveal melanoma are loss of chromosome 3 and gains of 8q and 6p [3].

Despite the improvements in diagnosis and the development of more effective local therapies for primary tumors, the rate of metastatic death remains unchanged. Unfortunately, once uveal melanoma has spread to distant organs, the disease is largely resistant to currently available therapies [4]. Nowadays, many new prognostic factors such as Cytological Features, Standard Karyotyping, Fluorescence in situ Hybridization, Centromeric Probes, Single Nucleotide Polymorphism and Gene Expression Profiling were investigated. Following the technique advance and lower expense of gene expression microarray, it has become a useful tool for studying the development and progression of tumors owing to its high throughout, new genes which may affect metastases of uveal melanoma could be found.

Several gene expression profiling studies have identified two molecular classes strongly associated with metastatic risk [57]. Hepatocyte growth factor/scatter factor (HGF) [810], Insulin-like growth factor [11] and Stem cell factor [1214] receptors have been involved in metastatic progression of uveal melanoma. In addition, the chemokine receptor CXCR4 was recently related to liver homing of human uveal melanoma metastatic cells [15, 16]. However, little is known about the core genes and their potential mechanisms in liver metastases. The present study made use of bioinformatics method to analyze the data obtained from two public available datasets in combination with the clinical data about metastases of uveal melanoma patients in attempt to investigate the liver metastases-related genes. The different gene ontology and pathways would indicate the most important mechanisms and candidate genes in the process of liver metastases, and helpful in working out more specific and individualized target treatment regimens according to genetic characteristics of individual patients.


Tissue samples and clinical data

Tissue samples were obtained from two datasets in total 92 uveal melanomas after enucleation surgery upon approval of the institutional bioethics board. The patients in the two datasets were composed of 63 and 29 samples, respectively [17, 18]. All the samples were performed on Affymetrix Human Genome U133 Plus 2.0 Array. Both of two data series were accessible at NCBI GEO database, accession GSE22138 and GSE27831. Clinical, pathological, and molecular features of the tumors are presented in Table 1.

Table 1 Clinicopathologic characteristics of uveal melanoma patients

Significant differential gene analysis

All 92 tumors in two datasets were pooled and reanalyzed on Affymetrix Expression Console Software (Version 1.1). MAS5 was used to normalize the original data. One sample was removed after the normalized data filtered with Pearson’s Correlation (Additional file 1: Figure S1), and the remaining 91 samples were renormalized. Genes were standardized and interpreted functionally before comparison. Using t-test, SAM or RVM mode and the tumors with no metastases as the control group, the P value and the fold change were calculated for each differentially expressed gene. With a threshold of P value < 0.05 or FDR value < 0.05 and fold change > =1.5, related genes were picked out and the venn diagram was showed in Additional file 2: Figure S2. Based on our sample size and previous study [19], t-test result was chose for further analyze. Unsupervised hierarchical clustering was performed with Cluster using Pearson’s correlation distance metric and average linkage followed by visualization in Treeview [20].

Gene ontology (GO) analysis

Based on Gene Ontology Database (, the significant level of GOs of the liver metastases-related differentially expressed genes was analyzed by two-side Fisher’s exact test and χ 2 test using DAVID ( analysis [21]. The differential expression genes were analyzed independently according to up- and down-regulation of these genes. We computed P-values for all the differential expression genes in all GO categories, and the threshold of significance was defined as P-value < 0.05.

Pathway analysis

Based on KEGG ( database, the significant level of pathways of the liver metastases-related differentially expressed genes was analyzed by Pathway-Express [22, 23]. Significant differences from the expected were calculated with a two-sided binomial distribution. The numbers of genes corresponding to each pathway category among the differentially expressed genes was tallied and compared with the number of genes expected for each pathway category based on their representation on the Affymetrix Human Genome U133 plus 2.0 array. All signaling pathways were analyzed for the significance level, using gamma P < 0.05 as the threshold.

Signal-net analysis

Using java that allows users to build and analyze molecular networks, network maps were constructed. For instance, if there is confirmative evidence that two genes interact with each other, an interaction edge is assigned between the two genes. The considered evidence is the source of the interaction database from KEGG. Networks are stored and presented as graphs, where nodes are mainly genes (protein, compound, etc.) and edges represent relation types between the nodes, e.g. activation or phosphorylation. The graph nature of Networks raised our interest to investigate them with powerful tools implemented in R.

Data analysis

Numerical data were presented as means and standard deviation (± SD). Differences between means were analyzed using Student’s t test. All statistical analyses were performed using SPSS11.0 software (Chicago, IL).


Clinical characteristics of the two group samples

Comparing with the non-metastases uveal melanoma samples, the liver metastases group has no differences on patients’ age, gender or tumor diameter. The tumor thickness value was much larger in the liver metastases samples than in the non-metastases group. Though the tumor location, extrascleral extension, tumor cell type and chromosome 3 data were incomplete, the first two appeared no differences while the mixed subtype tumor and chromosome 3 monosomy showed positive correlation with liver metastases of uveal melanoma (Table 1).

Liver metastases-related differential expression genes

Using liver metastases as the demarcation, genes of non-metastases and liver metastases groups were compared and 1138 statistically significant differential expression probes were obtained. Of them, 768 probes were up-regulated in the non-metastases samples and 370 were down-regulated. Hierarchical clustering showed systematic variations in the expression of genes between the two groups (Figure 1). The results demonstrated these differential probes could clearly separate the two groups from the whole samples and have good consistency in the group. Different genes with the most obvious p value and fold change between the liver metastases and non-metastases group were listed (Tables 2 and 3).

Figure 1
figure 1

Unsupervised classification of uveal melanoma samples based on gene expression profiling. Classification of 91 uveal melanoma samples using the 1138-probe sets identified as differentially expressed between the 46 liver metastatic samples and the 45 non-metastatic samples. Expression data are depicted as a data matrix where each row represents a gene and each column represents a sample. Expression levels are depicted according to the color scale shown at the top. Red and green indicate expression levels, respectively, above and below the median. The magnitude of deviation from the median is represented by the color saturation.

Table 2 Most obviously dysregulated genes sorted by P value in liver metastatic uveal melanoma compare to non-metastatic tumors
Table 3 Most obviously dysregulated genes with fold change absolute value > 3 in liver Metastatic Uveal Melanoma compare to non-metastatic tumors

Significant GOs

GOs of the differential expression genes were statistically analyzed. It was found that the differential expression genes obtained from the Microarray mainly participated in 125 significant GOs. According to the enrichment list, the up-regulated differential expression genes mainly participated in 89 GOs including germ cell programmed cell death, developmental programmed cell death, germ cell migration, melanocyte differentiation (Table 4), and the down-regulated differential expression genes mainly participated in 36 GOs including embryonic skeletal system development, regulation of Rho protein signal transduction, regulation of cell morphogenesis (Table 5).

Table 4 Different gene significant upregulated GO
Table 5 Different gene significant downregulated GO

Significant Pathways

The pathways of liver metastatic uveal melanoma samples were analyzed according to the functions and interactions of the differential genes. By using Pathway-Express which contains both the up- and down- regulated differential genes in its analysis and the threshold of significance defined on the basis of gamma P-value < 0.05, 16 significant pathways were found. Many signaling pathways had been verified to relate to cell migration and invasion, including Phosphatidylinositol signaling system, Gap junction, Adherens junction (Figure 2).

Figure 2
figure 2

Histogram of signal pathways those were significantly different in liver metastatic and non-metastatic uveal melanoma showing. X axis, negative logarithm of the P value (-LgP); Y axis, the name of the pathway. The larger the -LgP, the smaller the P value.


According to the literature and experimental records in the databases, a diagram of the gene interaction network was drawn up based on the genes differentially expressed between liver metastases and non-metastases uveal melanoma (Figure 3). The total number of genes in the network was 297, and the particular relationships between them were listed (Additional file 3). In the network, cycle nodes represent genes, and edges between two nodes represent interactions between genes, which were quantified by degree. Degrees within the network which describe the number of single gene that regulates other genes represent the size of the cycle node. The higher the degree, the more central the gene occurs within the network. There were many more genes upregulated involved in the Signal-Net analysis, of which YWHAZ, PRKDC, and ESR1 were the three main overexpressed genes, while the three main underexpressed genes were FYN, TIAM1 and GNAI1.

Figure 3
figure 3

Signal transduction networks of liver metastatic related genes. Circle represents genes, red circle represents the upregulated gene, and blue circle represents the downregulated gene. Arrow represents the activation of (a); straight line represents combine; dotted line represents indirect effects; a represents activation; ex represents gene expression; b represents binding; ind represents indirect effects; inh represents inhibition; u represents ubiquination. For the most interesting nodes, the fold change and p value were added to the YWHAC and FYN.


The present study followed up microarray-based 91 uveal melanoma patients, of whom patients were not significantly different in age, gender and tumor diameter. More than 1100 probes were differentially dysregulated in liver metastases uveal melanoma in this study. Based on sorting the different genes by p value and fold change, numbers of genes with the largest difference individually were carried out. HTR2B is a type of serotonin receptor with both the extremely significant p value and the largest upregulated fold change in the two groups. Soll et al. had found serotonin can promote tumor growth in hepatocellular cancer [24] and this maybe an effective hint to us. CHL1 is a cell adhesion molecule with homology to L1CAM involved in the regulation of cell adhesion and migration. Current study had supported a significant role for CHL1 gene in the growth, migration and invasion of human cervical cancer cells [25]. MEGF10 and OVOS2 have the largest downregulated fold change and a significant p value. It is still lack of information about their roles both in tumor and ophthalmology. Zinc finger proteins are among the most abundant proteins which were involved in transcriptional activation, apoptosis regulation and protein folding. Some zinc finger protein had been detected related to melanoma [26]. However, there have been no reports of a differential expression of ZNF family members in liver metastases uveal melanoma. Comparing with the original papers of the expression microarray, we found PTP4A3/PRL3 still has significant difference while SDCBP dropped out with too little fold change. In our opinion, this is caused by the enlargement of sample number which may impact larger in the previously smaller sample group for GSE27831 has only 29 samples. For the sample number we integrate here is not so small, though these genes were still lack of information about their impact on uveal melanoma liver metastases, we have the reason to focus on them in further investigation.

The GO is widely recognized as the premier tool for the organization and functional annotation of molecular aspect [27]. GO-analysis was used to interpret each GO of differential expressed gene and analyzed it statistically. By using the criteria of P < 0.05, we obtained the significant GOs and genes involved in them. Guo et al. used GO-analysis to analyze miRNA microarray and found that miR-15b and miR-16 may be indispensable for apoptosis by targeting Bcl-2 [28]. GO terms about programmed cell death plays an important role in liver metastases, many studies had reported proteins such as TRAIL and TGF-β1 involved in cell apoptosis also affected cancer cell metastases [29, 30]. Ion transport related genes is little investigated in cancer cell metastases, but recently Lee et al. found monoamine carboxylate transporters could impact colon fibrosarcoma cell migration through regulating extracellular pH value [31]. This may be a new way investigating uveal melanoma liver metastases. Also the structural changes of intracellular organelles and cytoskeletal would of course act on cell migration and invasion. According to our results, programmed cell death, ion transport and cytoskeletal would affect cell environment to play roles in liver metastases. Furthermore, other biological process may also have their effects in uveal melanoma distant metastases.

Pathways analysis can show the distinct biological process and find significant pathways that differential expression genes participate in, based on which can we have a comprehensive understanding about interactions of genes, functions that they participate in and relations between up- and down-stream, and obtain genes involved in these significant pathways. Appearance of pathways about cell junction, melanogenesis and calcium signaling pathway confirm their concordance with GO terms and their critical role in liver metastases. Numerous studies had proved PI3K/Akt signaling pathway which belongs to phosphatidylinositol signaling pathway participated in different cancers’ liver metastases [3234]. Holder et al. had reported gap junction may affect cancer metastases since connections are made between the primary tumor cells and foreign host cells at the secondary metastatic site [35]. Found of epimorphin activating focal adhesion kinase/extracellular signal-regulated kinase/matrix metalloproteinase-9 axis to promote hepatocellular carcinoma invasion and metastases verified focal adhesion and ECM-receptor interaction’s role in liver metastases. For this, we have reasons to believe the other seemingly irrelevant pathways would have their functions in uveal melanoma metastases.

Investigating genes involved in significant GOs and pathways, 297 genes in common were found that may affect the liver metastases of uveal melanoma patients. YWHAZ and FYN were identified as key genes in uveal melanoma liver metastases and play crucial roles in cell proliferation, which may be related to the higher incidence of metastases in uveal melanoma. YWHAZ previously was thought as a reference gene in many cell lines, recently it was found to play a major role in YWHAZ/beta-Catenin Axis and promote Epithelial-Mesenchymal Transition and Lung Cancer Metastases [36, 37]. FYN was a member of the Src family of kinases and though lack of data about its role in uveal melanoma, it participated in many cancers metastases through different pathways [3840]. Based on these data, further studies of these genes’ expression and the protein functions of HTR2B, CHL1, the ZNF family, YWHAZ and FYN need to be performed in more samples using reverse transcriptase-polymerase chain reaction and western blotting; moreover, for the GO and pathway findings are based on the currently known database, the regulation of identified genes and protein functions may have some other ways.

The above results all suggest that differences in gene expression exist between liver metastases and non liver metastases uveal melanoma. These genes encode proteins involved in different GOs and signal pathways, the disruption of which can promote cancer metastases. Several genes, such as HTR2B, CHL1, the ZNF family, YWHAZ and FYN provide potential candidates for distinguishing between uveal melanoma whether contain liver metastases in the future. This distinction will aid in the diagnosis and prevention of uveal melanoma liver metastases, based on their different features. Therefore, our results may provide important referential merit for clinical investigation. Nevertheless, the genes and the related GOs and pathways identified here are required to be further dissected and confirmed in more patient samples by other clinic-related studies.


Bioinformatics may help excavate and analyze large amounts of data in microarrays by means of rigorous experimental planning, scientific statistical analysis and collection of complete data about liver metastases of uveal melanoma patients. In the present study, a novel differential gene expression pattern was constructed and advanced study will provide new targets for diagnosis and mechanism of uveal melanoma liver metastases.


  1. Shields JA, Shields CL, De Potter P, Singh AD: Diagnosis and treatment of uveal melanoma. Semin Oncol. 1996, 23 (6): 763-767.

    CAS  PubMed  Google Scholar 

  2. Sibbritt T, Patel HR, Preiss T: Mapping and significance of the mRNA methylome. Wiley Interdiscip Rev RNA. 2013, 4 (4): 397-422. 10.1002/wrna.1166.

    Article  CAS  PubMed  Google Scholar 

  3. Trolet J, Hupe P, Huon I, Lebigot I, Decraene C, Delattre O, Sastre-Garau X, Saule S, Thiery JP, Plancher C, Asselain B, Desjardins L, Mariani P, Piperno-Neumann S, Barillot E, Couturier J: Genomic profiling and identification of high-risk uveal melanoma by array CGH analysis of primary tumors and liver metastases. Invest Ophthalmol Vis Sci. 2009, 50 (6): 2572-2580. 10.1167/iovs.08-2296.

    Article  PubMed  Google Scholar 

  4. Augsburger JJ, Correa ZM, Shaikh AH: Effectiveness of treatments for metastatic uveal melanoma. Am J Ophthalmol. 2009, 148 (1): 119-127. 10.1016/j.ajo.2009.01.023.

    Article  PubMed  Google Scholar 

  5. Onken MD, Worley LA, Ehlers JP, Harbour JW: Gene expression profiling in uveal melanoma reveals two molecular classes and predicts metastatic death. Cancer Res. 2004, 64 (20): 7205-7209. 10.1158/0008-5472.CAN-04-1750.

    Article  CAS  PubMed  Google Scholar 

  6. Onken MD, Ehlers JP, Worley LA, Makita J, Yokota Y, Harbour JW: Functional gene expression analysis uncovers phenotypic switch in aggressive uveal melanomas. Cancer Res. 2006, 66 (9): 4602-4609. 10.1158/0008-5472.CAN-05-4196.

    Article  CAS  PubMed  Google Scholar 

  7. Petrausch U, Martus P, Tonnies H, Bechrakis NE, Lenze D, Wansel S, Hummel M, Bornfeld N, Thiel E, Foerster MH, Keilholz U: Significance of gene expression analysis in uveal melanoma in comparison to standard risk factors for risk assessment of subsequent metastases. Eye (Lond). 2008, 22 (8): 997-1007. 10.1038/sj.eye.6702779.

    Article  CAS  Google Scholar 

  8. Mallikarjuna K, Pushparaj V, Biswas J, Krishnakumar S: Expression of epidermal growth factor receptor, ezrin, hepatocyte growth factor, and c-Met in uveal melanoma: an immunohistochemical study. Curr Eye Res. 2007, 32 (3): 281-290. 10.1080/02713680601161220.

    Article  CAS  PubMed  Google Scholar 

  9. Di Cesare S, Marshall JC, Logan P, Antecka E, Faingold D, Maloney SC, Burnier MN: Expression and migratory analysis of 5 human uveal melanoma cell lines for CXCL12, CXCL8, CXCL1, and HGF. J Carcinog. 2007, 6: 2-

    PubMed  PubMed Central  Google Scholar 

  10. Peruzzi B, Bottaro DP: Targeting the c-met signaling pathway in cancer. Clin Cancer Res. 2006, 12 (12): 3657-3660. 10.1158/1078-0432.CCR-06-0818.

    Article  CAS  PubMed  Google Scholar 

  11. All-Ericsson C, Girnita L, Seregard S, Bartolazzi A, Jager MJ, Larsson O: Insulin-like growth factor-1 receptor in uveal melanoma: a predictor for metastatic disease and a potential therapeutic target. Invest Ophthalmol Vis Sci. 2002, 43 (1): 1-8.

    PubMed  Google Scholar 

  12. Mouriaux F, Kherrouche Z, Maurage CA, Demailly FX, Labalette P, Saule S: Expression of the c-kit receptor in choroidal melanomas. Melanoma Res. 2003, 13 (2): 161-166. 10.1097/00008390-200304000-00008.

    Article  CAS  PubMed  Google Scholar 

  13. Mouriaux F, Chahud F, Maurage CA, Malecaze F, Labalette P: Implication of stem cell factor in the proliferation of choroidal melanocytes. Exp Eye Res. 2001, 73 (2): 151-157. 10.1006/exer.2001.1017.

    Article  CAS  PubMed  Google Scholar 

  14. Lefevre G, Glotin AL, Calipel A, Mouriaux F, Tran T, Kherrouche Z, Maurage CA, Auclair C, Mascarelli F: Roles of stem cell factor/c-Kit and effects of Glivec/STI571 in human uveal melanoma cell tumorigenesis. J Biol Chem. 2004, 279 (30): 31769-31779. 10.1074/jbc.M403907200.

    Article  CAS  PubMed  Google Scholar 

  15. Scala S, Ierano C, Ottaiano A, Franco R, La Mura A, Liguori G, Mascolo M, Staibano S, Ascierto PA, Botti G, De Rosa G, Castello G: CXC chemokine receptor 4 is expressed in uveal malignant melanoma and correlates with the epithelioid-mixed cell type. Cancer Immunol Immunother. 2007, 56 (10): 1589-1595. 10.1007/s00262-007-0303-6.

    Article  CAS  PubMed  Google Scholar 

  16. Franco R, Botti G, Mascolo M, Loquercio G, Liguori G, Ilardi G, Losito S, La Mura A, Calemma R, Ierano C, Bryce J, D'Alterio C, Scala S: CXCR4-CXCL12 and VEGF correlate to uveal melanoma progression. Front Biosci (Elite Ed). 2010, 2: 13-21.

    Article  Google Scholar 

  17. Gangemi R, Mirisola V, Barisione G, Fabbi M, Brizzolara A, Lanza F, Mosci C, Salvi S, Gualco M, Truini M, Angelini G, Boccardo S, Cilli M, Airoldi I, Queirolo P, Jager MJ, Daga A, Pfeffer U, Ferrini S: Mda-9/syntenin is expressed in uveal melanoma and correlates with metastatic progression. Plos One. 2012, 7 (1): e29989-10.1371/journal.pone.0029989.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Laurent C, Valet F, Planque N, Silveri L, Maacha S, Anezo O, Hupe P, Plancher C, Reyes C, Albaud B, Rapinat A, Gentien D, Couturier J, Sastre-Garau X, Desjardins L, Thiery JP, Roman-Roman S, Asselain B, Barillot E, Piperno-Neumann S, Saule S: High PTP4A3 phosphatase expression correlates with metastatic risk in uveal melanoma patients. Cancer Res. 2011, 71 (3): 666-674. 10.1158/0008-5472.CAN-10-0605.

    Article  CAS  PubMed  Google Scholar 

  19. Jeanmougin M, de Reynies A, Marisa L, Paccard C, Nuel G, Guedj M: Should we abandon the t-test in the analysis of gene expression microarray data: a comparison of variance modeling strategies. Plos One. 2010, 5 (9): e12336-10.1371/journal.pone.0012336.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4 (5): 3-10.1186/gb-2003-4-5-p3.

    Article  Google Scholar 

  22. Khatri P, Draghici S, Ostermeier GC, Krawetz SA: Profiling gene expression using onto-express. Genomics. 2002, 79 (2): 266-270. 10.1006/geno.2002.6698.

    Article  CAS  PubMed  Google Scholar 

  23. Draghici S, Khatri P, Martins RP, Ostermeier GC, Krawetz SA: Global functional profiling of gene expression. Genomics. 2003, 81 (2): 98-104. 10.1016/S0888-7543(02)00021-6.

    Article  CAS  PubMed  Google Scholar 

  24. Soll C, Jang JH, Riener MO, Moritz W, Wild PJ, Graf R, Clavien PA: Serotonin promotes tumor growth in human hepatocellular cancer. Hepatology. 2010, 51 (4): 1244-1254. 10.1002/hep.23441.

    Article  CAS  PubMed  Google Scholar 

  25. Long MJ, Wu FX, Li P, Liu M, Li X, Tang H: MicroRNA-10a targets CHL1 and promotes cell growth, migration and invasion in human cervical cancer cells. Cancer Lett. 2012, 324 (2): 186-196. 10.1016/j.canlet.2012.05.022.

    Article  CAS  PubMed  Google Scholar 

  26. Yajima I, Kumasaka M, Thang ND, Yanagishita T, Ohgami N, Kallenberg D, Naito Y, Yoshikawa T, Sakashita N, Kato M: Zinc finger protein 28 as a novel melanoma-related molecule. J Dermatol Sci. 2009, 55 (1): 68-70. 10.1016/j.jdermsci.2009.02.010.

    Article  CAS  PubMed  Google Scholar 

  27. Lovering RC, Camon EB, Blake JA, Diehl AD: Access to immunology through the gene ontology. Immunology. 2008, 125 (2): 154-160. 10.1111/j.1365-2567.2008.02940.x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Guo CJ, Pan Q, Li DG, Sun H, Liu BW: miR-15b and miR-16 are implicated in activation of the rat hepatic stellate cell: an essential role for apoptosis. J Hepatol. 2009, 50 (4): 766-778. 10.1016/j.jhep.2008.11.025.

    Article  CAS  PubMed  Google Scholar 

  29. Yu SC, Leung TW, Lee KT, Wong LK: Angioplasty and stenting of intracranial atherosclerosis with the wingspan system: 1-year clinical and radiological outcome in a single Asian center. J Neurointerv Surg. 2014, 6: 96-102. 10.1136/neurintsurg-2012-010608.

    Article  PubMed  Google Scholar 

  30. Armstrong SR, Campbell CB, Richardson CL, Vickery RG, Tsuruda PR, Long DD, Hegde SS, Beattie DT: The in vivo pharmacodynamics of the novel opioid receptor antagonist, TD-1211, in models of opioid-induced gastrointestinal and CNS activity. Naunyn Schmiedebergs Arch Pharmacol. 2013, 386 (6): 471-478. 10.1007/s00210-013-0844-5.

    Article  CAS  PubMed  Google Scholar 

  31. Festa F, Steel J, Bian X, Labaer J: High-throughput cloning and expression library creation for functional proteomics. Proteomics. 2013, 13 (9): 1381-1399. 10.1002/pmic.201200456.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Lublin FD: MS as a gateway disease. J Neurol Sci. 2013, 333 (1-2): 73-75. 10.1016/j.jns.2013.02.014.

    Article  PubMed  Google Scholar 

  33. Grana RA: Electronic cigarettes: a new nicotine gateway?. J Adolesc Health. 2013, 52 (2): 135-136. 10.1016/j.jadohealth.2012.11.007.

    Article  PubMed  Google Scholar 

  34. Shukla S, Chandran S, Gadagkar R: Ovarian developmental variation in the primitively eusocial wasp ropalidia marginata suggests a gateway to worker ontogeny and the evolution of sociality. J Exp Biol. 2013, 216 (Pt 2): 181-187.

    Article  PubMed  Google Scholar 

  35. Mueller MB, Blunk T, Appel B, Maschke A, Goepferich A, Zellner J, Englert C, Prantl L, Kujat R, Nerlich M, Angele P: Insulin is essential for in vitro chondrogenesis of mesenchymal progenitor cells and influences chondrogenesis in a dose-dependent manner. Int Orthop. 2013, 37 (1): 153-158. 10.1007/s00264-012-1726-z.

    Article  PubMed  Google Scholar 

  36. Chen CH, Chuang SM, Yang MF, Liao JW, Yu SL, Chen JJ: A novel function of YWHAZ/beta-catenin axis in promoting epithelial-mesenchymal transition and lung cancer metastasis. Mol Cancer Res. 2012, 10 (10): 1319-1331. 10.1158/1541-7786.MCR-12-0189.

    Article  CAS  PubMed  Google Scholar 

  37. Chua SL, See Too WC, Khoo BY, Few LL: UBC and YWHAZ as suitable reference genes for accurate normalisation of gene expression using MCF7, HCT116 and HepG2 cell lines. Cytotechnology. 2011, 63 (6): 645-654. 10.1007/s10616-011-9383-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Wang Q, Qian J, Wang F, Ma Z: Cellular prion protein accelerates colorectal cancer metastasis via the Fyn-SP1-SATB1 axis. Oncol Rep. 2012, 28 (6): 2029-2034.

    CAS  PubMed  Google Scholar 

  39. Yadav V, Denning MF: Fyn is induced by Ras/PI3K/Akt signaling and is required for enhanced invasion/migration. Mol Carcinog. 2011, 50 (5): 346-352. 10.1002/mc.20716.

    Article  CAS  PubMed  Google Scholar 

  40. Kim AN, Jeon WK, Lim KH, Lee HY, Kim WJ, Kim BC: Fyn mediates transforming growth factor-beta1-induced down-regulation of E-cadherin in human A549 lung cancer cells. Biochem Biophys Res Commun. 2011, 407 (1): 181-184. 10.1016/j.bbrc.2011.02.134.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references


This study was supported by the 11th five-year plan of science and technology research projects of China (no. 2007BAI18B07).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Lei Chen or Jianhong Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YZ carried out the molecular genetic studies, participated in the gene alignment and drafted the manuscript. YY participated in the design of the study and performed the statistical analysis. LC participated in the gene alignment and its design and coordination. JZ conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Figure S1: Normalized microarray data by Pearson’s Correlation. Uveal melanoma samples microarray data were normalized by Pearson’s Correlation and one sample GSM685470 was removed. (TIFF 4 MB)


Additional file 2: Figure S2: Venn diagram for different methods of gene expression analysis. T-test, SAM and RVM mode were used to screen the differential expressed genes. (TIFF 354 KB)


Additional file 3: Gene feature and relationship in the liver metastatic related signal transduction network.(XLS 150 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Reprints and Permissions

About this article

Cite this article

Zhang, Y., Yang, Y., Chen, L. et al. Expression analysis of genes and pathways associated with liver metastases of the uveal melanoma. BMC Med Genet 15, 29 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: