A Web-based database of genetic association studies in cutaneous melanoma enhanced with network-driven data exploration tools

Abstract

Introduction

Although a small fraction of mainly familial patients with cutaneous melanoma (CM) carry highly penetrant gene mutations ( 1–5 ), i.e. mutations in CDKN2A and CDK4 , the more common sporadic form of CM is likely caused by the complex interplay of environmental and multiple genetic risk factors that exert moderate risk effects ( 1 ). Up to now, an increasing number of genetic association studies including candidate-gene and genome-wide association studies (GWAS) have been published that report on novel CM risk genes or attempt to validate previously reported genetic risk factors in CM ( 1 ). Owing to the enormous amount of partly contradictory information derived from those studies, the evaluation and interpretation of the genetic predisposition of CM is not a trivial task ( 6 ). To effectively collect and analyze all available information related to CM, we have implemented an online database named MelGene that provides a systematic and in-depth qualitative and quantitative catalog of genetic association studies in CM. Our database includes random-effects meta- analysis results of eligible polymorphisms that highlight the most compelling CM risk loci ( 7 ).

We conducted a systematic update in the MelGene database by including detailed summaries of all recently published association studies and by performing meta- analyses in all eligible polymorphisms that have been investigated in multiple studies to provide a summary effect for the association of each single-nucleotide polymorphism (SNP) to CM risk. The epidemiological validity of nominally significant meta-analysis results was assessed using the ‘Venice’ criteria suggested by the Human Genome Epidemiology Network ( 8 ). In this study, we present the new design of the front end and the back end of the database, providing the user with a more functional interface, easy-to-handle queries and embedded tools that facilitate the visualization and the exploration of the molecular relationship networks of putative genetic risk factors of CM. MelGene database provides a systematic and comprehensive overview of genetic association studies (both candidate-gene and GWAS) focusing exclusively on CM. Besides the embedded tools for data searching using keywords, MelGene database provides tools for automated meta- analysis of the collected data and allows for construction of networks using the available polymorphisms in the database. These features make MelGene database unique compared with databases such as GWAScentral ( http://www.gwascentral.org/ ) or ‘The catalogue of Published Genome-Wide Association studies’ from the National Human Genome Research Institute that act as a compilation of summary-level findings and allow for searches between loci derived from GWA studies for various outcomes ( http://www.genome.gov/gwastudies/ ).

Material and methods

Search strategy, data collection and meta-analysis

For the continuous curation of the MelGene database, we performed systematic literature searches for peer-reviewed genetic association studies on CM using PubMed ( http://www.ncbi.nlm.nih.gov/pubmed ), the Human Genome and Epidemiology Network Navigator ( http://hugenavigator.net ) and the Melanoma Molecular Maps Project ( http://www.mmmp.org/MMMP ). The last search was conducted on 31 August 2013. The search strategy has been described in detail elsewhere ( 7 ). The current version of MelGene includes 192 publications that fulfilled our inclusion criteria [outlined in ref. ( 7 )] and that report on 1114 polymorphisms across 280 genes. For each biallelic polymorphism included in the database with data of at least four independent case-control data sets (n = 79), a random- effects meta-analysis was calculated based on the DerSimonian and Laird model ( 7 ). In addition, between-study heterogeneity was quantified by the I ² metric. Forest plots of the respective meta-analysis results were created using the ‘ metafor ’ package ( 9 ) in the R programming language ( http://www.r-project.org/ ). Users are able to download those forest plots in high resolution ( Figure 1 B and C). The database curation is supported by an experienced team, which includes clinicians, biologists, bioinformaticians, biostatisticians and genetic epidemiologists.

Figure 1.

( A ) Updated MelGene database search engine. Users are able to retrieve information available on MelGene based on keywords such as the gene name, polymorphism name, chromosome, first author of a publication, year of publication, ethnicity and the country of origin of study populations. ( B ) Polymorphism overview page and meta-analysis of polymorphism rs1042602 as an example. All publications that were included in MelGene and assessed rs1042602 in their data sets are listed in a sortable interactive table, and cross-links to the corresponding publications indexed in PubMed are provided in the database. ( C ) Forest plot of rs1042602 displays study-specific results as well as the summary OR , 95% CI and heterogeneity estimate.

Open in new tabDownload slide

Database construction

The database scheme was created using MySQL (version 5.5.27, http://www.mysql.com/ ) and comprises the following fields: entry id, gene symbol, chromosome, location, study name consisting of first author name and year of publication, ethnicity, population and polymorphism name (where applicable, the official NCBI ’s rs identifiers have been used, http://www.ncbi.nlm.nih.gov/snp/ ). In addition, fields were also supplemented by the following information where available: number of melanoma cases, number of controls, significance assessment by the authors of each publication, the minor and the major allele name based on genetic data available in the respective publication, minor–minor, minor–major and major–major genotype counts per study population, the allele frequency in CM cases as well as in control subjects, the additive odds ratio ( OR ) and 95% confidence interval (CI) limits if provided in the respective publication.

Web application interface

The updated publicly available version of MelGene enables users to search the database based on a variety of parameters. More precisely, the database can be searched by gene name, polymorphism name, chromosome, name of the first author of a publication, the year of publication, the minimal number of cases per population, the geographical origin of study populations and by using a free text keyword search field. The updated Web-application search engine was implemented using html , PHP ( http://php.net/ ) and MySQL ( http://www.mysql.com/ ) queries ( Figure 1 A).

Moreover, for each polymorphism included in the MelGene , the updated database version provides links to several other genetic databases, and thus facilitates the retrieval of additional information on specific polymorphisms and genes of interest, i.e. NCBI ’s dbSNP ( http://www.ncbi.nlm.nih.gov/projects/SNP/ ), the International HapMap project’s database ( http://hapmap.ncbi.nlm.nih.gov/ ), the Ensembl browser ( www.ensembl.org ), the SNPedia Web site ( http://www.snpedia.com ) and GWAS Central ( https://www.gwascentral.org/ ).

Embedded tools to visualize and explore molecular relationship networks

Pairs of polymorphisms and/or their corresponding genes that have been included in MelGene can be assessed for 13 different types of relationship. Each type of relationship drives to a different network representation where the nodes are either polymorphisms or genes and the edges represent the existence and the strength of the relationship between two nodes. Network representations of the MelGene data were created and uploaded on the MelGene Web server. The software used for the construction of each network and the type of the networks are provided in Table 1 .

The list of all 280 genes included in MelGene has been used as input to Cytoscape ( http://www.cytoscape.org/ ) ( 10 ), an open-source software for integration, visualization and analysis of biological networks. Specifically, the GeneMANIA plug-in ( http://www.genemania.org/ ) ( 11 ) of Cytoscape has been used. This plug-in retrieves a list of genes that are related to the input genes based on a large set of functional data, including (i) co-expression, where two genes are linked if their expression levels are similar across conditions in a gene expression study, (ii) protein–protein interactions, (iii) genetic interaction, where two genes are functionally linked if perturbations of one gene has an impact to the second, (iv) shared protein domains, where two genes are linked if they have the same protein domain, (v) co-localization, where genes are linked if they are expressed in the same tissue and (vi) pathways, where two genes are linked if they are part of the same pathway [see ( 11 ) for more details]. With the use of GeneMANIA , new members of a pathway or a complex interaction can be highlighted.

Moreover, two additional relationship networks have been constructed. The first network was a protein–protein interaction network created by means of the STRING database ( http://string-db.org/ ) ( 12 ), whereas a second network was created based on the question of whether the input genes have been reported or predicted as being involved in CM pathophysiology in the literature, according to SABiosciences Gene Network Central ( 13 ).

In addition, all 1114 polymorphisms included in MelGene were submitted to SNAP ( http://www.broad institute.org/mpg/snap/ ) ( 14 ), yielding another six sets of networks. SNAP provides pair-wise SNP’s calculations of input against the ‘1000 Genome Pilot 1’ SNP data set based on phased genotype data from the International HapMap Project ( http://hapmap.ncbi.nlm.nih.gov/ ) ( 15 ). Using the SNAP tool with the 1114 SNPs, five additional molecular relationship networks were created using as network relationship measure (i) the recombination rate in centimorgans per million bases, (ii) the genetic map distance (that is, the distance from the query SNP to the proxy SNP in centimorgans), (iii) the genetic map position (that is, the position of the SNP on the genetic map for this chromosome in centimorgans), (iv) D’ [measure of linkage disequilibrium ( LD ) normalized to allele frequency] and (v) the r² correlation coefficient [see ( 14 ) for more details].

Both SNP and gene interaction network visualization tools were implemented in MelGene by means of the jQuery ( http://jquery.com/ ) and ‘ sigma.js ’ ( http://sigmajs.org/ ) JavaScript libraries. Figure 2 shows an example of an interaction network based on genes’ co-expression properties. Networks can be created on a random or circular layout. The node’s size is proportional to the number of first neighbors that are interconnected to the specific node. A larger size signifies a larger number of interconnections. Users of MelGene can locate and highlight a specific node interactively, either by using a drop-down box or by clicking on the node. The color and the weight of the edge between two interconnected nodes are an indicator of how closely these two nodes are related. White color corresponds to low, whereas cyan to high linkage weight. In addition, when a node is selected, node-specific network features [see ( 16 ) for detailed description] of their first neighbors are also provided, including the degree (i.e. number of adjacent edges), closeness (which measures how many steps are required to access every other node from a given node), betweenness (i.e. the shortest path from one node to another), and eigenvector centrality (a natural extension of degree centrality measuring the importance of a node if it is linked to by other important nodes). All the aforementioned features were apriori calculated by means of the ‘ igraph ’ ( http://cran.r-project.org/web/packages/igraph/ ) ( 16 ) R package tool ( http://www.r-project.org/ ) applied to each created relationship network. For each network, users are also able to apply network visualization filtering based on the mean value and the standard deviation calculated by the entire selected feature vector. This filtering can be performed either locally by selecting a node or globally to the entire network and can reveal ‘hub’ nodes in a more robust approach ( Figure 3 ).

Figure 2.

( A ) Dynamic gene interaction network based on protein–protein interactions using the gene TP53 as an example. First neighborhood interactions of TP53 are highlighted against the entire protein–protein interaction network. In the present example, edges (interconnections) with weight <20% of the maximum weight value were omitted, and node (gene) labels with >10 interactions are displayed. ( B ) List of the calculated node-specific network features corresponding to node TP53 is provided.

Open in new tabDownload slide

Figure 3.

Screenshot of the gene network created using GeneMANIA physical interactions, filtered with degree ≥ μ + σ, were μ and σ are the mean value and standard deviation of the entire degree feature vector, respectively.

Open in new tabDownload slide

Results and discussion

After a comprehensive data collection and systematic meta-analyses of CM association studies in MelGene , 20 genes showed genome-wide significant ( P  < 5 × 10 ⁻⁸ ) evidence for association with CM risk ( MelGene top genes , Table 2 ). In addition, summary OR s and 95% CIs, heterogeneity as well as summary estimates on exclusion of the first published study were calculated and the respective plots have been made publicly available in the database. Moreover, the new network visualization tools developed in the updated MelGene database provides additional information about interconnections among genes or SNPs that may be indirectly related to CM. In the following result sections, we depict some genes and polymorphisms highlighted from the network-driven analyses that were performed in the MelGene environment.

Table 2

Open in new tab

.Ranked list of the most significantly associated CM genes according to MelGene .

Network features were calculated on the relationship network produced by the STRING platform. In the present example, we have performed gene ranking for each network feature (of note, the pseudogene AFG3L1 was not found in the STRING platform and thus excluded from this list). The genes that rank among the top 50 of the respective network features are highlighted in yellow.

Table 2

Open in new tab

.Ranked list of the most significantly associated CM genes according to MelGene .

Network features were calculated on the relationship network produced by the STRING platform. In the present example, we have performed gene ranking for each network feature (of note, the pseudogene AFG3L1 was not found in the STRING platform and thus excluded from this list). The genes that rank among the top 50 of the respective network features are highlighted in yellow.

The gene co-expression network reconstructed by the GeneMANIA plug-in revealed that the majority of MelGene top genes (18 of 20) are interconnected with at least one edge, presenting a high level of co-expression between these genes. In Figure 4 , we see the corresponding co-expression matrix for these genes.

Figure 4.

The co-expression matrix for the MelGene top genes that showed genome-wide significant ( P  < 5x10 ⁻⁸ ) evidence for association with CM risk.

Open in new tabDownload slide

Furthermore, in the same gene co-expression network, we investigated which genes act as ‘hubs’ and assessed their relationship with the MelGene top genes . We calculated four network properties per node: degree, closeness, betweenness and eigenvector centrality. In the sequel, we calculated their average values (μ) and the corresponding standard deviations (σ). For each property, we considered as ‘hubs’ the genes that present a property value above μ + σ, and finally we merged the ‘hub genes’ from the four properties in one list. Afterward, we examined whether genes from Table 2 are ‘hub genes’ or are they connected with other ‘hub genes’. We found that the genes MITF , TYR and PARP1 of Table 2 are ‘hub genes’ and that the genes CDKN2A , OCA2 , MX2 , VDR , MITF and CASP8 are interconnected with >10 ‘hub genes’ as shown in Table 3 .

Table 3.

Open in new tab

Ranked list of the most significantly associated CM genes according to MelGene along with the number and the names of the interconnected ‘hub genes’ in the corresponding co-expression network

Yellow highlighted are those genes that act also as ‘hubs’. * AFG3L1 is reported as pseudogene.

Table 3.

Open in new tab

Ranked list of the most significantly associated CM genes according to MelGene along with the number and the names of the interconnected ‘hub genes’ in the corresponding co-expression network

Yellow highlighted are those genes that act also as ‘hubs’. * AFG3L1 is reported as pseudogene.

The current update of MelGene database includes 20 genes that showed genome-wide significant ( P  < 5 × 10 ⁻⁸ ) evidence for association with CM risk compared with 12 genes described in the previous field synopsis ( 7 ). Investigating the relationship between the previous and the current MelGene top gene list, we mapped the genes on the network constructed by GeneMANIA with all the available relationship types (co-expression, co-localization, genetic interactions, physical interactions, pathways and share protein domain). The fully connected (i.e. nodes without a connection are absent) subnetwork of these genes is shown in Figure 5 . The significant genes of the previous update are visualized with orange diamonds; the significant genes of the current field-synopsis are visualized in yellow whereas the genes required for a fully connected subnetwork are highlighted in blue ( TP53 , PMEL , ITGA7 and ITGA9 ; Figure 5 ).

Figure 5.

The fully connected subnetwork of the genes from two successive field synopses of MelGene that showed genome-wide significant evidence for association with CM risk. Orange, blue and yellow corresponding to the significant genes of the previous MelGene version, the significant genes of the current MelGene analysis and the genes required for a fully connected subnetwork, respectively.

Open in new tabDownload slide

We found that the genes TYRP1 , TYR , MITF and ERCC5 displayed in Figure 5 are ‘important connectors’ in the network as they interact with 3, 5, 6 and 3 ‘significant’ genes that were significant using random-effects meta- analysis. Investigating these genes further, we found that they have been implicated in CM in the literature and/or by several other databases. More specifically, p53 has been suggested to be a major player suppressing progression from nevi to melanoma ( 17 , 18 ). Aside from its role in pigmentation, PMEL (SILV) encodes antigenic epitopes that are recognized by multiple melanoma diagnostic antibodies including HMB-45, currently one of the most commonly used melanocytic markers for clinical melanoma diagnosis in humans ( 19 ). Furthermore, for PMEL , MalaCards, a database of human maladies and their annotations ( www.malacards.org ) ( 20 ), as well as the DISEASES database (Disease-gene associations mined from literature) developed by the University of Copenhagen ( http://diseases.jensenlab.org ) rank CM as the first disease related for PMEL . Finally, proteins encoded by ITGA7 and ITGA9 genes belong to the integrin alpha chain family. Changes in integrins expression have been reported during the malignant progression of many tumors, and much evidence exists implicating their involvement in CM metastasis ( 21 ). Malacards also implicates ITGA7 and ITGA9 in CM (eighth and third rank of gene-related diseases, respectively).

Another important gene in the co-expression network is MLANA . MLANA , a key network element, was found to be co-expressed with 8 ( TYRP1 , CDKN2A , OCA2 , SLC45A2 , PLA2G6 , MX2 , MITF and TYR ) from the 20 MelGene top genes of the current MelGene gene list. CM is ranking first among the diseases related to MLANA according to the MalaCards and DISEASES databases, but its role as a risk gene for CM in MelGene currently remains unclear due to lack of sufficient association data.

In addition, ITGB2 was found to act as a ‘hub gene’ by all centrality measures and is co-expressed with 5 ( CDKN2A , MX2 , VDR , MITF and CASP8 ) from the 20 MelGene top genes . Moreover, CM is ranking sixth among the diseases related to ITGB2 according to Malacards.

Finally, regarding HLA-A , CM is ranking third among the related diseases according to MalaCards and the DISEASES database. HLA-A is co-expressed with 2 ( MX2 and VDR ) of the 20 MelGene top genes and it acts as a ‘hub gene’ according to three important centrality metrics (Degree, Closeness and Eigenvector centrality).

In a similar way, we examined the networks of single-nucleotide polymorphisms (SNPs) based on their DPrime value with proxy SNPs. Specifically, we found the SNPs acting as ‘hubs’ in these networks and investigated their relationship with the top 20 SNPs that correspond to the 20 MelGene top genes . Again, we calculated the four network centrality measures. The SNPs that present at least one network property value above μ + σ are considered to act as ‘hub’ SNPs. We merged the ‘hub’ SNPs derived from the four different metrics resulting in 607 unique ‘hub’ SNPs, eight of which belong to the top 20 SNPs of MelGene (rs401681, rs1408799, rs2218220, rs6001027, rs11263498, rs1393350, rs1801516 and rs3219090— Table 4 , highlighted in yellow). The SNP4Disease database ( http://snp4disease.mpi-bn.mpg. de/ ), which was developed by the Max Planck Institute for Heart and Lung Research and provides information on diseases linked to SNPs by literature-mining techniques from various sources, was used to find which of the 607 ‘hub’ SNPs have been implicated in CM. In all, 72 ‘hub’ SNPs have been linked to CM, seven of them belong to the top 20 SNPs. Furthermore, we examined whether the top 20 SNPs are connected with ‘hub’ SNPs in the DPrime network (rs14961795, located in MITF , was not found in SNP DPrime network and thus excluded from this list). Five of the top 20 SNPs interact with at least one ‘hub’ SNP (rs408799, rs2218220, rs6001027, rs1126349 and rs1393350) as shown in Table 4 .

Table 4.

Open in new tab

Ranked list of the most significantly associated CM SNPs according to MelGene along with the number and the names of the interconnected ‘hub SNPs’ in the corresponding DPrime SNP network

Yellow highlighted are those SNPs that act also as ‘hubs’. * Indicates that the specific SNP does not exist on the DPrime SNP network.

Table 4.

Open in new tab

Ranked list of the most significantly associated CM SNPs according to MelGene along with the number and the names of the interconnected ‘hub SNPs’ in the corresponding DPrime SNP network

Yellow highlighted are those SNPs that act also as ‘hubs’. * Indicates that the specific SNP does not exist on the DPrime SNP network.

When exploring the 72 ‘hub’ SNPs that were implicated in CM by the SNP4Disease database in the polymorphism networks based on SNP DPrime values, we observed that four SNPs [rs2218220 (located in MTAP ), rs4636294 (located in MTAP ), rs854145 (located in GRM1 ) and rs935053 (located Near MTAP )] acted as ‘hubs’ based on three centrality measures (degree, closeness and betweenness centrality). We observed that one of these SNPs, i.e. rs854145 in GRM1 , is the most central polymorphism (see Figure 6 A). On the other hand, rs2218220, which is located in MTAP , belongs to the 20 most significant SNPs in the MelGene meta-analysis results. Finally, rs4636294 and rs935053 were also found as ‘hub’ SNPs from the three centrality metrics. As shown in Figure 6 B, they are also highly intercorrelated with rs2218220.

Figure 6.

Screenshots of SNP networks created based on the DPrime SNP network of ( A ) rs854145 with no visualization filtering, ( B ) rs935053 with Degree ≥ μ + 3σ.

Open in new tabDownload slide

To highlight the most important ‘hub’ SNPs, we applied the same procedure but we considered as ‘hubs’ the SNPs that present at least one centrality measure above μ + 3σ. We obtained 256 ‘hub’ SNPs for the 4 centrality measures of which 2 are from the 20 most significant SNPs in the MelGene meta-analysis results (rs11263498 [located in CCND1 ] and rs2218220 [located in MTAP ]). From the 256 ‘hub’ SNPs, 21 are reported in the ‘SNP4Disease’ database in relation to CM. Strikingly, the 21 SNPs derived with more rigorous centrality filtering (at least one centrality measure above μ + 3σ) included again the same central polymorphisms (rs854145, rs2218220, rs4636294 and rs935053) found previously in the 72 ‘hub’ SNPs that were implicated in CM by the SNP4Disease database with at least one centrality measure above μ + σ.

In summary, by using several pipelines network-driven data exploration tools implemented in MelGene , we have identified seven genes from the MelGene database (none of which had sufficient association data to perform a meta-analysis in MelGene ) and four polymorphisms (of which three have not been meta-analyzed in MelGene because of lack of association data) that may play a potential role in CM pathophysiology ( Table 5 ).

Conclusion and future work

We present here the new re-designed version of the MelGene database empowered with network analysis tools that has led to the identification of known and several new promising genes implicated in melanoma pathophysiology. MelGene serves as a comprehensive reference repository for genetic association data in CM and can provide further insights in the predisposition of CM by systems biology approaches.

In the future, it is planned to further expanding the MelGene database by curating our data set with new small- and large-scale genetic association studies in regular time intervals. In addition, our database can be expanded so as to include and integrate more data in the future including somatic mutations in melanoma from publicly available databases such as the COSMIC database ( 22 ), as well as more and of different types of pre-compiled relationship networks.

Funding

This research has been co-financed by the European Union (European Social Fund – ESF) and Greek national funds through the Operational Program ‘Education and Lifelong Learning’ of the National Strategic Reference Framework (NSRF) - Research Funding Program: Aristeia I - 1094. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflict of interest . None declared.

References

Author notes