This article has Open Peer Review reports available.
Combinations of newly confirmed Glioma-Associated loci link regions on chromosomes 1 and 9 to increased disease risk
© Yang et al; licensee BioMed Central Ltd. 2011
Received: 26 July 2010
Accepted: 9 August 2011
Published: 9 August 2011
Glioblastoma multiforme (GBM) tends to occur between the ages of 45 and 70. This relatively early onset and its poor prognosis make the impact of GBM on public health far greater than would be suggested by its relatively low frequency. Tissue and blood samples have now been collected for a number of populations, and predisposing alleles have been sought by several different genome-wide association (GWA) studies. The Cancer Genome Atlas (TCGA) at NIH has also collected a considerable amount of data. Because of the low concordance between the results obtained using different populations, only 14 predisposing single nucleotide polymorphism (SNP) candidates in five genomic regions have been replicated in two or more studies. The purpose of this paper is to present an improved approach to biomarker identification.
Association analysis was performed with control of population stratifications using the EIGENSTRAT package, under the null hypothesis of "no association between GBM and control SNP genotypes," based on an additive inheritance model. Genes that are strongly correlated with identified SNPs were determined by linkage disequilibrium (LD) or expression quantitative trait locus (eQTL) analysis. A new approach that combines meta-analysis and pathway enrichment analysis identified additional genes.
(i) A meta-analysis of SNP data from TCGA and the Adult Glioma Study identifies 12 predisposing SNP candidates, seven of which are reported for the first time. These SNPs fall in five genomic regions (5p15.33, 9p21.3, 1p21.2, 3q26.2 and 7p15.3), three of which have not been previously reported. (ii) 25 genes are strongly correlated with these 12 SNPs, eight of which are known to be cancer-associated. (iii) The relative risk for GBM is highest for risk allele combinations on chromosomes 1 and 9. (iv) A combined meta-analysis/pathway analysis identified an additional four genes. All of these have been identified as cancer-related, but have not been previously associated with glioma. (v) Some SNPs that do not occur reproducibly across populations are in reproducible (invariant) pathways, suggesting that they affect the same biological process, and that population discordance can be partially resolved by evaluating processes rather than genes.
We have uncovered 29 glioma-associated gene candidates; 12 of them known to be cancer related (p = 1. 4 × 10-6), providing additional statistical support for the relevance of the new candidates. This additional information on risk loci is potentially important for identifying Caucasian individuals at risk for glioma, and for assessing relative risk.
Determining the molecular changes that underlie phenotypic distinctions is a major thrust of cell biology. More specifically, identifying the precise DNA alterations in the genes and regulatory regions that underlie predisposition, initiation and progression of tumors is a central theme of biomedical research. Understanding the molecular changes associated with initiation and progression requires tissue samples from the tumor itself which are often difficult to obtain, as well as from a suitable control population. On the other hand, understanding molecular changes associated with predisposition requires only genomic DNA (e.g., from white blood cells) from the target and control populations, which can be obtained relatively readily. In this manuscript we focus on the latter, since that is where the preponderance of available information is. The methods can, however, be easily extended to the study of somatic genomic associations as control tissue samples from the brain become available.
Identifying predisposition involves (i) identification of the approximate genomic location of a change correlating with phenotypic distinction, which is usually done by finding correlative single nucleotide polymorphisms (SNPs), followed by (ii) the identification of genes or promoters in strong linkage disequilibrium with the SNPs, i.e. those that are coinherited. Following this is (iii) a search for mechanisms, such as point mutations, deletions, and translocations, which can be carried out by sequencing identified genomic regions in a sufficiently large number of samples from affected and control populations.
Here we focus on identifying regions and genes that predispose to glioblastoma multiforme (GBM) (i.e. (i) and (ii), above) by conducting a genome-wide association (GWA) study. A number of such studies have already been carried out for GBM and, as is typical for such studies, very few genes have been consistently identified across different populations [1, 2].
The Cancer Genome Atlas (TCGA) samples
Normal European American blood samples (n = 1366) were downloaded from the Illumina iControlDB (iControls) as the comparison group. After screening, 1306 control samples remained. Because the iControl has more than 3000 samples, we treat it as a background population for our TCGA analysis and use it as a control group (that is also used in AGS). Its treatment as a background population means that we would not expect the results to change with another sample from the same background population. As a partial check of this statement, we divided the glioma samples from TCGA into 2 independent samples, conducting a GWA analysis on each using the same iControl population. If the use of a single control created bias, we'd expect overlapping results. In fact the results have no overlap SNP at the 10-6 significance level (Additional file 1).
Association analysis was performed using the EIGENSTRAT package, under the null hypothesis of "no association between the GBM and control SNP genotypes" based on an additive inheritance model . The significance threshold p for association needed to be set stringently to allow for the large number of multiple tests in the GWA study. In particular, if we require that the probability of 1 or more false positives be less than 0.05, we must have 1-e -Np ≤ 0.05 or p ≤ 0.93 × 10-7 ≈ 10-7 with N, the total number of SNPs examined, taken as 550,000. At this level, few if any SNPs will be detected for typical glioma population sizes. The alternative is to accept a less stringent p-value, and to eliminate false discoveries by seeking confirmation in an independent study.
where we take sample-specific weights (W i ) proportional to the square root of the "total number of individuals" and Z i = F -1(1-p i ); F -1 (.) is an inverse standard normal CDF (cumulative distribution function). The false discovery rate is estimated as the fused probability multiplied by the total number of SNPs, which was 300,000 in the AGS study.
Concordant SNPs recovered from TCGA and AGS data, and associated genes
Genes correlated with SNP (R2)
We use the coefficient of determination, R 2, to identify the genes with strong eQTL correlation with the SNPs, identified by meta-analysis. R 2 was calculated based on the correlation between gene expression level and SNP genotype, using the "SCAN: SNP and CNV Annotation Database" resource (http://scan.bsd.uchicago.edu/newinterface/index.html). Genes/SNP combinations with R 2 greater than 0.8 were considered to be strongly correlated.
The odds ratio for glioma susceptibility is defined and calculated as follows. For any pair or triplet of the above-mentioned risk alleles, define:
n 11: number of individuals with the given risk allele combination in the glioma sample
n 12: number of individuals without the risk allele combination in the glioma sample
n 21: number of individuals with the risk allele combination in the control sample
n 22: number of individuals without the risk allele combination in the control sample
Pairwise and triplet SNP combinations with odds ratios greater than 3
Identification of Associated Pathways and Genes
Tumor initiation is associated with alterations in physiological processes that involve sets of genes and allelic variants in any of several such genes, not all revealed in a single population.
The standard method for identifying altered processes is a pathway enrichment analysis, which can be carried out using a single population . In this case, pathways would be identified by showing that the number of significant SNPs/genes occurring in a particular pathway is above chance expectation. The procedure that we describe here extends this methodology to multiple populations. The assignment of a SNP/gene to a particular pathway using a single population is done using a significance threshold which is loose enough to allow multiple assignments from that population, but not stringent enough for an acceptable FDR in the single population. The FDR is brought down to an acceptable level, as described below, when both populations assign the same gene(s) to the same pathway.
The procedure begins as follows: (1) Identify SNPs having a p-value < 10-3 in either of the populations. (2) Identify genes that include these SNPs. (3) Assign the genes thus obtained to KEGG  pathways.
In any given pathway, genes identified by the AGS SNPs are generally different from those identified by our TCGA SNPs. Because the p-value is not stringent (p < 0.001), there is a reasonable chance that a number of SNPs, and therefore pathways, are false positives. We reduced the likelihood of false positives by determining whether, for a given pathway, the number of genes that are common to the two datasets is greater than expected by chance. In a particular pathway, if n 1 genes are identified by TCGA data, and n 2 by AGS (all at nominal p < 0.001), we calculate the hypergeometric probability of finding at least n genes (the observed number) common to the two sets. If we set an FDR = 0.05, we eliminate all pathways for which the hypergeometric probability exceeds 0.05 divided by the number of pathways to which genes were assigned.
Significant SNP candidates
Wrensch et al.  screened the Adult Glioma Study (AGS) population using p = 10-6 and inferred 13 SNPs, 3 of which were confirmed in the Mayo Clinic population at a multiple hypothesis corrected p-value of 0.0038 (0.05/13). In a similar fashion, we validated 4 of these 13 SNPs in the TCGA dataset: rs2736100 at chromosome 5p15.33, rs1412829 at 9p21.3, and rs7530361 and rs501700 at 1p21.2 (Table 1, boldface). The last two SNPs are reported as validations for the first time. The above confirmation-based method of correction used for multiple hypotheses is an approximation which gives primacy to one of the populations. The difficulty can be illustrated by first screening on the TCGA data set rather than on the AGS dataset. Using p = 10-6, two SNPs are found for the TCGA screen: rs11840214 in gene EFNB2 at 13q33 (p = 1.08 × 10-7) and rs1909486 at 8q24 (p = 4.3 × 10-7), and neither of these can be confirmed on the AGS population.
Joint rather than sequential analysis of data from two or more populations can increase the power to detect genetic associations . In particular, using equation (1) in Methods, we identified 12 SNPs (Table 1), confirmed by AGS and TCGA at an FDR < 0.05, one of which was previously confirmed by Wrensch et al.  and Shete et al. . Of the 11 remaining SNPs, 4 have been reported by Shete et al., and the remaining 7 SNPs are reported for the first time here. The 12 SNPs are distributed over five genomic regions, including chromsomes 5q15.33, 9q21.3, 1p21.2, 3q26.2 and 7p15.3. Two of these, 5q15.33 and 9q21.3, have been reported in previous studies [1, 2]. The 12 SNPs are in strong linkage disequilibrium with 25 genes, 8 of which are known to be associated with cancer (indicated in boldface in Table 1). An additional SNP of interest is rs12341266 at 9q32, which has an FDR of 0.06 and is in the glioma associated gene, RGS3 (Additional file 2).
Odds ratio for combinatory SNPs
The analysis of TCGA and AGS identifies 12 SNPs, 7 of which have not been previously reported. One of the implications of identifying these additional SNPs is that the number of associated gene groups that can be used to estimate the odds ratio increases combinatorially. Consequently we can expect higher prognostic reliability for individuals possessing a combination of risk alleles, although at some loss of population coverage. We consider here all combinations of two and three SNPs, while constraining our choices to SNPs that are more than a mega-base (Mb) apart, in order to minimize redundant (disequilibrated) information. Specifically, the 12 SNPs are divided into 5 groups based on their chromosomal locations. Chromosome 1 has 5 SNPs clustered together within 1 Mb, and chromosome 9 has 4 SNPs within 1 Mb around genes CDKN2A, 2B. The remaining 3 SNPs are located on chromosomes 3, 5, and 7.
Genes Identified by Conserved Pathway Analysis
Pathways that contain significant SNPs (p < 10-3) inferred from both AGS and TCGA samples
Purine metabolism (p = 3.50E-04)**
Small cell lung cancer (p = 4.35E-04) **
Non-small cell lung cancer (p = 2.6E-04) **
Neuroactive ligand-receptor interaction (p = 8.00E-04) **
Vascular smooth muscle contraction (p = 3.48E-04) **
Gap junction (p = 1.30E-04) **
Long-term depression (p = 6.95E-04) **
Olfactory transduction (p = 3.47E-04) **
Axon guidance (p = 3.91E-04) **
Pathways in cancer (p = 2.13E-03)
Colorectal cancer (p = 8.69E-05) **
Focal adhesion (p = 1.95E-03)
ECM-receptor interaction (p = 8.69E-04) **
Cell adhesion molecules (CAMs) (p = 1.74E-04) **
Regulation of actin cytoskeleton (p = 1.56E-03)
Hypertrophic cardiomyopathy (HCM) (p = 1.22E-03)
Arrhythmogenic right ventricular cardiomyopathy (ARVC) (p = 9.12E-04) **
Dilated cardiomyopathy (p = 1.04E-03)
Genes strongly correlated with SNP candidates
Eight of the 25 genes we identified are known to be cancer-associated. These include TERT [1, 9, 10], SLC6A18 , CLPTM1L [9, 10], CDKN2A/B [1, 11, 12], SASS6 , ITGB8 , and MACC1  (Table 1). Five of the genes, TERT, SLC6A18, CLPTM1L, and CDKN2A/B, have previously been shown to be associated with glioma by other GWA studies.
Our combined GWA/pathway analysis predicts four additional genes that are identified in the literature as cancer-related. We therefore predict 29 glioma-associated genes, 12 of them known from previous studies to be cancer-related. It is useful to evaluate the probability that as many as 12 cancer-related genes in a set of 29 would be found by chance. If we use the fraction of OMIM genes that are cancer-related as an estimate of the background frequency of cancer genes in the disease gene population, the probability that 29 genes include 12 cancer-associated genes by chance is 1.4 × 10-6. The fraction of OMIM genes that are cancer-related is 0.1 (750 cancer-associated genes in 7381 OMIM genes).
Each of the 8 cancer-related genes listed above plays one or more key roles in processes known to be altered during tumor initiation and development . For example, MACC1 is a growth pathway regulator influencing angiogenesis and processes related to metastasis ; CDKN2A is a well-studied cell cycle regulator  and a known tumor suppressor whose loss results in a diminished ability to regulate growth and predisposition to cancer ; ITGB8 has been implicated in activities related to metastasis, including adhesion and migration ; and the telomerase enzyme (TERT) is linked to unlimited replication . It is worth noting that CDKN2A, 2B are in strong linkage disequilibrium with rs1412829 at 9p21.3, which has now been identified in 3 independent studies and should therefore be considered an unusually high confidence cancer gene marker.
As noted in the Background, the identification of disease-associated SNPs by GWA studies tends to have low concordance when different populations are compared. An example is the above-mentioned verification in the Mayo Clinic population of only 3 out of the 13 SNPs identified in the AGS population. Our own results identified 12 SNPs based on a meta-analysis of TCGA and AGS. These include 5 of 14 glioma-associated SNPs that had already been verified by at least 2 other studies. There are many reasons for this, including possible differences in ethnicity and other less obvious stratifiers, and differences in sample size and background populations. We believe, however, that an additional reason of an entirely different nature is present which, roughly stated, is that different genes contribute to the same phenotype.
Recent evidence based in clinical epidemiology, computational genomics and various model systems [17–21] suggests that disease phenotypes emerge from dysfunction of one or more components of a functionally coherent gene module. Importantly, the dysfunctions (such as mutations in a gene or its promoter, post-translational modifications, multiple copy number variations, translocations) need not be the same, nor need they be in the same gene, in different individuals with the same disease phenotype. Since alterations in different genes in the same functional module (e.g. a pathway) can lead to the same dysfunction, low reproducibility rates in GWA studies -- aimed at identifying DNA loci having inherited alterations that predispose to complex disorders -- are not surprising. We might, however, expect greater concordance between populations if aberrant functional modules of genes were compared.
In fact we do find that some SNPs not invariant across populations appear in pathways that are invariant, suggesting that the SNPs, however different, are linked to genes that contribute to the same biological process, and that population discordance of SNPs can be partially resolved by evaluating processes rather than genes. Among the 406 AGS SNPs having p-value < 0.001, only 3% (12 out of 406) are invariant; i.e. have significant fused p-values. This contrasts with 19% of the non-invariant SNPs (75 out of the remaining 394) being in one of the 13 invariant pathways. In addition, as noted above, five genes, each appearing in two or more of these pathways, are invariant across the two populations.
The crucial role of dysregulated signalling pathways in cancer development, and the frequent therapeutic targeting of kinases, makes association with the cyclic GMP-dependent protein kinase PRKG1 a plausible finding. This appears to be the first time this particular kinase has been associated with human cancer.
DCC, the deleted colorectal carcinoma gene, has (as its name implies) been well studied and has been implicated in a number of cancers in addition to colorectal cancer . DCC is an axon guidance receptor that responds to netrin-1 , and is a component of a pathway implicated in the regulation of angiogenesis, cell survival, apoptosis, and cell positioning and migration , all of which adds to the biological plausibility of this finding.
The GABA-A receptor gene family encodes the major inhibitory neurotransmitter receptors in the central nervous system; changes in GABA-A receptor function have been implicated in diseases as diverse as alcoholism, epilepsy, schizophrenia, autism and Alzheimer's disease. Over 19 different subunit genes code for the pentameric GABA-A receptors, pharmacologically distinct forms of inhibitory neurotransmitter receptors that change their expression levels either during development or in disease . Because of this complexity, elucidating molecular details of its role in many diseases is slow. The ligand GABA may be involved in metastatic prostate cancer , making association with aberrant forms of the receptor plausible, and the receptor itself has in fact been implicated in prostate cancer , and at least one of its subunits has been found to stimulate the development of pancreatic cancer . This brief context, and the fact that glioma is also a disease of the brain, suggest that the association we find is biologically plausible.
In fact the finding that all 4 genes are plausibly associated with cancer may itself be significant. If we use the conservative estimate that two of the four (FHIT and DCC) are taken as cancer-associated because of their inclusion in GAD, we obtain 0.0015 as the probability of chance occurrence. We therefore expect that although the procedure we outline is perhaps more heuristic than rigorous, it yields highly suggestive candidates, which should be of interest to glioma researchers.
It is evident that the greater the number of risk allele candidates, the greater the chance of identifying a potentially predisposed individual if all such genes are assayed. In addition, however, increasing knowledge of risk alleles affects not just the chance of identification, but the reliability, which will also increase to some extent. The reason for this is related to the chance of observing a combination of risk alleles in the same individual, and the fact that the number of combinations grows as a 2 r , where r is the number of known risk alleles. As an example we calculated the odds for all doublet and triplet combinations for the 12 SNPs obtained by a meta-analysis of the AGS and TCGA populations, and list those with odds ratio > 3 in Table S3. Perhaps the most interesting result is the extremely high odds ratio (OR = 5) associated with the two alleles rs1412829 and rs4977756, both of which are in a cyclin dependent kinase.
In this paper we demonstrate that a meta-analysis of SNP data from TCGA which have not been previously analyzed, together with data from the Adult Glioma Study, identify 12 glioma-associated SNPs which represent 5 genomic regions -- 5p15.33, 9p21.3, 1p21.2, 3q26.2 and 7p15.3 -- three of these not previously reported. Of the 12 SNPs identified in this study, 5 have been previously reported and verified by other studies [1, 2], while the remaining 7 are novel candidates. Eight genes known to be cancer-associated are included in, or are in strong linkage disequilibrium with, one or more of these SNPs. An additional 4 genes (PRKG1, FHIT, GABRG3 and DCC ) are identified by a combined pathway enrichment GWA analysis. In all we obtain 29 genes, of which 12 are known to be cancer related.
We have shown that the greatest relative risk occurs when risk alleles are present on chromosomes 1 and 9. Finally based on an analysis of the processes in which candidate genes are involved, we suggest that processes rather than genes are the most informative way to compare different populations, and have shown that such comparisons reduce discordance by approximately 19%.
This work was supported by TCGA grant CA135882-02S1.
- Wrensch M, Jenkins RB, Chang JS, Yeh RF, Xiao Y, Decker PA, Ballman KV, Berger M, Buckner JC, Chang S, et al: Variants in the CDKN2B and RTEL1 regions are associated with high-grade glioma susceptibility. Nature genetics. 2009, 41 (8): 905-908. 10.1038/ng.408.View ArticlePubMedPubMed CentralGoogle Scholar
- Shete S, Hosking FJ, Robertson LB, Dobbins SE, Sanson M, Malmer B, Simon M, Marie Y, Boisselier B, Delattre JY, et al: Genome-wide association study identifies five susceptibility loci for glioma. Nature genetics. 2009, 41 (8): 899-904. 10.1038/ng.407.View ArticlePubMedPubMed CentralGoogle Scholar
- Network CGAR: Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008, 455 (7216): 1061-1068. 10.1038/nature07385.View ArticleGoogle Scholar
- Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nature genetics. 2006, 38 (8): 904-909. 10.1038/ng1847.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, et al: American journal of human genetics. 2007, 81 (3): 559-575. 10.1086/519795.View ArticlePubMedPubMed CentralGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanehisa M: The KEGG database. Novartis Foundation symposium. 2002, 247: 91-101. discussion 101-103, 119-128, 244-152View ArticlePubMedGoogle Scholar
- Skol AD, Scott LJ, Abecasis GR, Boehnke M: Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nature genetics. 2006, 38 (2): 209-213. 10.1038/ng1706.View ArticlePubMedGoogle Scholar
- Kang JU, Koo SH, Kwon KC, Park JW, Kim JM: Gain at chromosomal region 5p15.33, containing TERT, is the most frequent genetic event in early stages of non-small cell lung cancer. Cancer genetics and cytogenetics. 2008, 182 (1): 1-11. 10.1016/j.cancergencyto.2007.12.004.View ArticlePubMedGoogle Scholar
- Stacey SN, Sulem P, Masson G, Gudjonsson SA, Thorleifsson G, Jakobsdottir M, Sigurdsson A, Gudbjartsson DF, Sigurgeirsson B, Benediktsdottir KR, et al: New common variants affecting susceptibility to basal cell carcinoma. Nature genetics. 2009, 41 (8): 909-914. 10.1038/ng.412.View ArticlePubMedPubMed CentralGoogle Scholar
- Bisio A, Nasti S, Jordan JJ, Gargiulo S, Pastorino L, Provenzani A, Quattrone A, Queirolo P, Bianchi-Scarra G, Ghiorzo P, et al: Functional analysis of CDKN2A/p16INK4a 5'-UTR variants predisposing to melanoma. Human molecular genetics. 2010, 19 (8): 1479-1491. 10.1093/hmg/ddq022.View ArticlePubMedGoogle Scholar
- Sherr CJ: Cancer cell cycles. Science (New York, NY. 1996, 274 (5293): 1672-1677. 10.1126/science.274.5293.1672.View ArticleGoogle Scholar
- Leidel S, Delattre M, Cerutti L, Baumer K, Gonczy P: SAS-6 defines a protein family required for centrosome duplication in C. elegans and in human cells. Nature cell biology. 2005, 7 (2): 115-125. 10.1038/ncb1220.View ArticlePubMedGoogle Scholar
- Culhane AC, Quackenbush J: Confounding effects in "A six-gene signature predicting breast cancer lung metastasis". Cancer research. 2009, 69 (18): 7480-7485. 10.1158/0008-5472.CAN-08-3350.View ArticlePubMedPubMed CentralGoogle Scholar
- Boardman LA: Overexpression of MACC1 leads to downstream activation of HGF/MET and potentiates metastasis and recurrence of colorectal cancer. Genome medicine. 2009, 1 (4): 36.View ArticlePubMedPubMed CentralGoogle Scholar
- Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2000, 100 (1): 57-70. 10.1016/S0092-8674(00)81683-9.View ArticlePubMedGoogle Scholar
- Rzhetsky A, Wajngurt D, Park N, Zheng T: Probing genetic overlap among complex human phenotypes. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (28): 11694-11699. 10.1073/pnas.0704820104.View ArticlePubMedPubMed CentralGoogle Scholar
- Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (21): 8685-8690. 10.1073/pnas.0701361104.View ArticlePubMedPubMed CentralGoogle Scholar
- Park J, Lee DS, Christakis NA, Barabasi AL: The impact of cellular networks on disease comorbidity. Molecular systems biology. 2009, 5: 262.View ArticlePubMedPubMed CentralGoogle Scholar
- Linghu B, Snitkin ES, Hu Z, Xia Y, Delisi C: Genome-wide prioritization of disease genes and identification of disease-disease associations from an integrated human functional linkage network. Genome biology. 2009, 10 (9): R91-10.1186/gb-2009-10-9-r91.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee DS, Park J, Kay KA, Christakis NA, Oltvai ZN, Barabasi AL: The implications of human metabolic network topology for disease comorbidity. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (29): 9880-9885. 10.1073/pnas.0802208105.View ArticlePubMedPubMed CentralGoogle Scholar
- McKusick VA: Mendelian Inheritance in Man. A Catalog of Human Genes and Genetic Disorders. 1998, Johns Hopkins University Press, 12Google Scholar
- Becker KG, Barnes KC, Bright TJ, Wang SA: The genetic association database. Nature genetics. 2004, 36 (5): 431-432. 10.1038/ng0504-431.View ArticlePubMedGoogle Scholar
- Huebner K, Garrison PN, Barnes LD, Croce CM: The role of the FHIT/FRA3B locus in cancer. Annual review of genetics. 1998, 32: 7-31. 10.1146/annurev.genet.32.1.7.View ArticlePubMedGoogle Scholar
- Frank S, Muller J, Plaschke J, Hahn M, Hampl J, Hampl M, Pistorius S, Schackert G, Schackert HK: The putative tumor suppressor gene FHIT at 3p14.2 is rarely affected by loss of heterozygosity in primary human brain tumors. Cancer research. 1997, 57 (13): 2638-2641.PubMedGoogle Scholar
- Reyes-Mugica M, Rieger-Christ K, Ohgaki H, Ekstrand BC, Helie M, Kleinman G, Yahanda A, Fearon ER, Kleihues P, Reale MA: Loss of DCC expression and glioma progression. Cancer research. 1997, 57 (3): 382-386.PubMedGoogle Scholar
- Bernet A, Mehlen P: Dependence receptors: when apoptosis controls tumor progression. Bulletin du cancer. 2007, 94 (4): E12-17.PubMedGoogle Scholar
- Hung JH, Whitfield TW, Yang TH, Hu Z, Weng Z, Delisi C: Identification of functional modules that correlate with phenotypic difference: the influence of network topology. Genome biology. 11 (2): R23.Google Scholar
- Whiting PJ, Bonnert TP, McKernan RM, Farrar S, Le Bourdelles B, Heavens RP, Smith DW, Hewson L, Rigby MR, Sirinathsinghji DJ, et al: Molecular and functional diversity of the expanding GABA-A receptor gene family. Annals of the New York Academy of Sciences. 1999, 868: 645-653. 10.1111/j.1749-6632.1999.tb11341.x.View ArticlePubMedGoogle Scholar
- Azuma H, Inamoto T, Sakamoto T, Kiyama S, Ubai T, Shinohara Y, Maemura K, Tsuji M, Segawa N, Masuda H, et al: Gamma-aminobutyric acid as a promoting factor of cancer metastasis; induction of matrix metalloproteinase production is potentially its underlying mechanism. Cancer research. 2003, 63 (23): 8090-8096.PubMedGoogle Scholar
- Abdul M, McCray SD, Hoosein NM: Expression of gamma-aminobutyric acid receptor (subtype A) in prostate cancer. Acta oncologica (Stockholm, Sweden). 2008, 47 (8): 1546-1550. 10.1080/02841860801961265.View ArticleGoogle Scholar
- Takehara A, Hosokawa M, Eguchi H, Ohigashi H, Ishikawa O, Nakamura Y, Nakagawa H: Gamma-aminobutyric acid (GABA) stimulates pancreatic cancer growth through overexpressing GABAA receptor pi subunit. Cancer research. 2007, 67 (20): 9704-9712. 10.1158/0008-5472.CAN-07-2099.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/4/63/prepub
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.