- Research article
- Open Access
- Open Peer Review
Similar gene expression profiles of sporadic, PGL2-, and SDHD-linked paragangliomas suggest a common pathway to tumorigenesis
BMC Medical Genomicsvolume 2, Article number: 25 (2009)
Paragangliomas of the head and neck are highly vascular and usually clinically benign tumors arising in the paraganglia of the autonomic nervous system. A significant number of cases (10–50%) are proven to be familial. Multiple genes encoding subunits of the mitochondrial succinate-dehydrogenase (SDH) complex are associated with hereditary paraganglioma: SDHB, SDHC and SDHD. Furthermore, a hereditary paraganglioma family has been identified with linkage to the PGL2 locus on 11q13. No SDH genes are known to be located in the 11q13 region, and the exact gene defect has not yet been identified in this family.
We have performed a RNA expression microarray study in sporadic, SDHD- and PGL2-linked head and neck paragangliomas in order to identify potential differences in gene expression leading to tumorigenesis in these genetically defined paraganglioma subgroups. We have focused our analysis on pathways and functional gene-groups that are known to be associated with SDH function and paraganglioma tumorigenesis, i.e. metabolism, hypoxia, and angiogenesis related pathways. We also evaluated gene clusters of interest on chromosome 11 (i.e. the PGL2 locus on 11q13 and the imprinted region 11p15).
We found remarkable similarity in overall gene expression profiles of SDHD -linked, PGL2-linked and sporadic paraganglioma. The supervised analysis on pathways implicated in PGL tumor formation also did not reveal significant differences in gene expression between these paraganglioma subgroups. Moreover, we were not able to detect differences in gene-expression of chromosome 11 regions of interest (i.e. 11q23, 11q13, 11p15).
The similarity in gene-expression profiles suggests that PGL2, like SDHD, is involved in the functionality of the SDH complex, and that tumor formation in these subgroups involves the same pathways as in SDH linked paragangliomas. We were not able to clarify the exact identity of PGL2 on 11q13. The lack of differential gene-expression of chromosome 11 genes might indicate that chromosome 11 loss, as demonstrated in SDHD-linked paragangliomas, is an important feature in the formation of paragangliomas regardless of their genetic background.
Paragangliomas are tumors originating in cells of neural crest origin in the extra-adrenal paraganglia associated with the autonomic nervous system. Most paragangliomas arise in the parasympathetic paraganglia of the head and neck region, but they can also arise in the parasympathetic paraganglia of the mediastinum or in the orthosympathetic para-aortic and retroperitoneal paraganglia. They are highly vascular and usually characterized by an indolent, non-invasive growth pattern. Most cases are sporadic, but a significant number (10–50%) have been shown to be familial. Mutations in 3 of the 4 genes encoding subunits of succinate dehydrogenase (SDH, complex II in the mitochondrial respiratory chain) have been implicated in the familial forms of the disease: SDHB, SDHC, and SDHD [1–3]. In our population, the majority of hereditary paraganglioma cases are associated with two founder mutations in the SDHD gene on 11q23. In addition to these SDH related cases, another hereditary paraganglioma family has been identified with linkage to a region on 11q13, the PGL2 locus. No mitochondrial complex II genes, including SDHA, are located in the 11q13 region, and the identity and function of the PGL2 gene are yet unknown. Mutations in SDHB, SDHC and SDHD are also implicated in the formation of phaeochromocytomas, tumors arising in cells derived from the neural crest in the adrenal medulla [6–8]. In PGL2 -linked cases no association with phaeochromocytoma formation has been found to date. A recent genome-wide expression study of phaeochromocytomas identified two distinct clusters: one containing SDH- and VHL-associated phaeochromocytomas and another containing MEN2- and NF1-associated phaeochromocytomas, while both clusters contained sporadic cases. The cluster containing VHL and SDH associated phaeochromocytomas was characterized by a transcription signature of reduced oxidoreductase activity and increased angiogenesis and hypoxia. In order to gain further insight into PGL2 function and identity, we have performed a gene expression study evaluating gene expression in head and neck paragangliomas of different genetic backgrounds: SDHD -linked, PGL2 -linked and sporadic cases without a mutation in the SDHB, SDHC or SDHD gene. In addition to a supervised gene-based analysis, a supervised pathway-based analysis was performed, evaluating differences in gene-expression for predefined pathways and functional gene groups. We evaluated in more detail gene groups that are known to be associated with SDH function and paraganglioma-or phaeochromocytoma formation, i.e. metabolism, cell cycle, hypoxia, and angiogenesis related pathways. In addition, we evaluated the gene sets that differentiate the SDH/VHL – from the NF1/MEN2-associated phaeochromocytoma cluster in the aforementioned phaeochromocytoma gene expression study, using our dataset. Finally, gene clusters located within or close to the PGL2 locus on 11q13, the SDHD locus on 11q23, and the imprinted 11p15 region were assessed. The latter region has previously been implicated in SDHD -linked paraganglioma formation. The results of both gene- and pathway- based analyses show remarkable similarity in the gene-expression profiles of SDHD -linked, PGL2 -linked and sporadic paragangliomas, suggesting that paraganglioma formation involves the same mechanisms and pathways in these paraganglioma subgroups.
Samples from head and neck paragangliomas were obtained from the tissue banks of the department of Pathology at the Leiden University Medical Center (LUMC) (all sporadic and SDHD related cases and one PGL2 -linked case) or the University Medical Center (UMC) St. Radboud (all but one PGL2 -linked cases). All specimens were handled according to the ethical guidelines, as described in the Code for Proper Secondary Use of Human Tissue in the Netherlands of the Dutch Federation of Medical Scientific Societies (FEDERA). Diagnosis of paraganglioma was confirmed by histology in all cases. All paragangliomas were carotid body tumors arising in the carotid bifurcation in the neck. No malignant paragangliomas were included in the study. Eighteen paraganglioma cases were selected: 7 cases with a known D92Y founder mutation in the SDHD gene, 6 cases from the family with significant linkage tot the PGL2 locus on 11q13, and 5 sporadic cases. The latter were defined as 'sporadic' because mutation scanning of SDHB, SDHC, and SDHD was negative, while the family histories of these cases were negative for HN-paraganglioma or any of the other clinical stigmata that would suggest the involvement of VHL, NF1 or the RET gene.
SDHB, SDHC, and SDHD genes were scanned for the presence of mutations at the laboratory for DNA diagnostics at the LUMC. All exonic regions of these genes were tested by direct sequencing using the Sanger method on an ABI 3177 Genetic Analyzer, starting with the exon containing the known Dutch founder mutations in SDHD followed by exons that had previously been found to contain pathogenic mutations in SDHD, SDHB, and SDHC (in that order) in the Dutch population[4, 11]. If that remained negative, scanning was completed by analyzing the remainder of exons of these genes. More recently, the sporadic, mutation-negative cases were also examined by MLPA for the presence of large deletions in SDHB, SDHC, and SDHD. MLPA was carried out with the P226 MLPA kit, containing probes for all exons and the promoter of each of these genes (27 different probes), according to the MRC Holland protocol.
RNA isolation and microarray hybridization
Tissue samples were snap frozen in liquid nitrogen and stored at -70°C. An experienced pathologist (PCWH) estimated the tumor percentage of the samples. Only samples with a tumor percentage of more than 70% were included in the study. Sample preparation was performed according to the Affymetrix protocol (Affymetrix, Inc., Santa Clara, CA). In brief, 30 5 μm sections were taken from each frozen tissue sample and total RNA was extracted using Trizol (Life Technologies, Inc., Rockville, MD), and purified using RNeasy columns according to the manufacturers protocol (Qiagen, Valencia, CA). A minimum of 10 μg of total RNA was used to synthesize cDNA with the Superscript Choice system (Life Technologies, Rockville, MD). First strand cDNA synthesis was performed with T7- (dT24) oligomer primer, followed by second strand synthesis using T4 DNA polymerase. The resultant was purified using Phase Lock Gel and precipitated in ethanol. Synthesis of biotine labeled cRNA was performed using the BioArray HighYield Transcript Labeling Kit (Enzo Diagnostics, Inc., Farmingdale, NY) according to the protocol of the manufacturer. In vitro transcription (IVT) reactions took place at 37°C for 4,5 hours. The labeled cRNA was purified using RNeasy columns (Qiagen, Valencia, CA) and fragmented in fragmentation buffer at 94°C for 35 minutes. Fragmented cRNA prepared from each individual sample was then transferred to a specialized Affymetrix hybridization centre (Leiden Genome Technology Centre, LGTC). Here the samples were hybridized according to the manufacturers protocol in a concentration of 0,5 μg/μl to a human GeneChip U95A-v2 (Affymetrix), containing approximately 8500 probe sets. The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus (GEO), and are accessible through GEO Series accession number GSE12921.
Sample size calculation
Sample size calculations were performed according to the method described by Pounds and Cheng.
Normalization and expression analysis
Acquisition and quantification of array images was performed using the MAS software package (Affymetrix). All arrays were normalized with gcrma normalization using the R statistical software package available from Bioconductor [17–19].
Unsupervised clustering analysis
The R package 'Linear Models for Microarray Data' (LIMMA) was used for the assessment of differential expression of individual genes between paraganglioma subgroups. Overall gene-expression differences between paraganglioma subgroups were evaluated with the 'global test' designed by J.J. Goeman using the R package 'global test' available on Bioconductor[18, 19, 21]. In order to evaluate subtle differences between paraganglioma subgroups, we analyzed all pathways in the Catalog of Human Gene Sets v2.0, containing 1687 gene sets, available from the Broad Institute as part of their publicly accessible Gene Set Enrichment Analysis (GSEA) software package [22, 23]. Instead of the statistical method used in the GSEA software, we used the global test developed by Goeman et al., because the latter tends to have more power to detect gene sets with small effect sizes [24–26]. Specific attention was paid to the gene sets that were significantly represented in SDH-linked phaeochromocytomas in a recent gene-expression study by Dahia et al.. Next, we applied the gene set that differentiated SDHD from MEN2-associated phaeochromocytomas in the aforementioned study to our data using the global test[9, 21]. Furthermore, we performed a pathway based analysis using the global test on manually curated gene sets, focusing specifically on pathways involved in processes or conditions that are known or assumed to play a role in paraganglioma formation, i.e. proliferation, survival, apoptosis, cell cycle regulation, metabolism and hypoxia, based on pathways described in literature and the publicly available pathway databases KEGG and Biocarta [27–29]. In addition to the evaluation of functionally related genes we also performed the global test on some topographically related gene groups on chromosome 11, i.e. the PGL2 minimal haplotype on 11q13, the SDHD region on 11q23, and 11p15, an imprinted region that has been implicated in SDHD -linked paraganglioma and phaeochromocytoma formation[10, 21]. In all, 264 manually curated pathways and functionally related gene sets were tested. All tests, both for genes and pathways, were corrected for multiple testing based on the false discovery rate (FDR) criterion, using the method of Benjamini and Hochberg.
Due to the rarity of PGL2-linked paragangliomas, sample sizes in this study are inevitably limited. In all, 21 samples were hybridized including 3 duplicates. Four samples (1 SDHD-linked sample, 2 PGL2-linked samples and 1 duplicate experiment) were excluded because of poor RNA or hybridization quality, leaving 15 different tumors in the analysis (5 sporadic, 6 SDHD-linked and 4 PGL2-linked samples) (Table 1).
Sample size calculation
Calculations showed that with this sample set and assuming that at least 30 to 35 genes are truly differentially expressed between subgroups with a fold change of 2.0 or more, at least 10 differentially expressed genes would be detected with a false discovery rate of 0.1.
Two-way hierarchical clustering of SDHD-linked, PGL2 -linked and sporadic paragangliomas revealed no clear clusters. No grouping according to genetic background was found (fig. 1). In fact, overall gene expression was very similar in all paraganglioma samples, with high correlation coefficients for overall gene-expression between all tumors irrespective of genetic background.
Using the LIMMA analysis, we did not find individual genes that are significantly differentially expressed between sporadic, SDHD and PGL2 -linked paragangliomas. The global test did not reveal significant differences in overall gene expression between paraganglioma subgroups. Using all 1687 functional gene sets from the Catalog of Human Gene Sets incorporated in the GSEA software, analysis with the global test revealed no significant differences in gene expression between SDHD and PGL2 -linked tumors, SDHD and sporadic tumors, or PGL2 and sporadic tumors for any gene set when corrected for multiple testing. In a recent phaeochromocytoma gene-expression study, several gene sets from the Catalog of Human Gene Sets were found to be significantly represented in SDH-associated phaeochromocytomas. These gene sets comprise microtubule activity, oxidoreductase activity, HIF1α, angiogenesis, proteasome degradation, electron transport chain, CCR3, collagen and glutathione metabolism. In our study, no significant differential expression between sporadic, SDHD and PGL2 -linked paragangliomas was found for these gene sets. Dahia et al. also identified a gene set differentiating SDH- from MEN2-associated phaeochromocytomas. This gene set contained 400 probes, encoding 288 different annotated genes. 212 of these 288 genes were also represented on the Affymetrix U95A chip used in this study. No significant differential expression between sporadic, SDHD and PGL2 -linked head and neck paragangliomas was observed for this gene set (data not shown). Next, we performed the global test on manually selected pathways assumed to play a role in paraganglioma formation, i.e. proliferation-, survival-, apoptosis, cell cycle regulation-, metabolism- and hypoxia related pathways. In all, 264 pathways and functional gene sets were tested. No significant differential expression was observed for any of these gene sets between the paraganglioma subgroups (data partially shown in fig. 2). Last, we performed a more detailed evaluation of genes located on chromosome 11 loci of interest (11q23, 11q13 and 11p15). This analysis also did not reveal significant differences between paraganglioma subgroups (data partially shown in fig. 3).
In our gene-expression analysis of sporadic, SDHD – and PGL2 -linked paragangliomas of the head and neck, no significant differences in gene-expression profile were observed between these genetically defined paraganglioma subgroups. Instead, we found considerable similarity between PGL2 -, SDHD – and sporadic tumor samples in both unsupervised and supervised analyses (fig. 1, 2 and 3). This correlates well with the observation that sporadic as well as SDHD -linked and PGL2 linked paragangliomas of the head and neck share important clinical characteristics like the age of onset of symptoms, the indolent growth pattern, and a usually benign behaviour of the tumor, although multiple paragangliomas are less often observed in sporadic cases [31–34]. Furthermore, all head and neck paraganglioma subtypes share the typical histological architecture of the 'zellballen', groups of neoplastic chief cells surrounded by sustentacular cells[35, 36].
In a recent gene-expression study by Dahia et al. of sporadic, SDHB-, SDHD -, VHL-, MEN2- and NF1-associated phaeochromocytomas, two phaeochromocytoma clusters were identified: a cluster containing VHL-, and SDH- linked tumors and another containing MEN2- and NF1-linked tumors. Gene set enrichment analysis showed that microtubule activity, oxidoreductase activity, HIF1α, angiogenesis, proteasome degeneration, electron transport chain, chemokine CCR3, collagen and glutathione metabolism gene sets were significantly represented in the gene expression signature of SDH-linked phaeochromocytomas. In our study, the GSEA pathway-based supervised analysis of sporadic, PGL2 – and SDHD -linked paragangliomas did not reveal significant differences between the subgroups for all GSEA gene sets, including the gene sets that characterized SDH tumors in the study by Dahia et al.. The authors also identified a gene set that differentiated SDH-linked tumors from MEN2-linked phaeochromocytomas. When applying this differentiating gene set to our dataset, significant differences in gene expression could not be found. These findings suggest that all paraganglioma subgroups in our study share the characteristics that defined the SDH-linked tumors in the study by Dahia et al., i.e. a signature of hypoxia, reduced oxidoreductase, and increased angiogenesis. Further characterization of the gene expression profiles of head and neck paragangliomas would require comparison with normal paraganglionic tissue. However, due to the microscopic size of normal paraganglia and their close anatomical relations with essential nerves and blood vessels it is not feasible to acquire this in sufficient quantity and quality to reliably perform RNA-based tests such gene-expression microarrays. In the present study, more detailed analysis of manually selected pathways and functional gene sets that are assumed to play a role in paraganglioma formation, i.e. processes of metabolism, angiogenesis and hypoxia as well as proliferation, survival, apoptosis and cell cycle related pathways also did not reveal significant differential expression between sporadic, SDHD -linked and PGL2 -linked paragangliomas. A striking finding is that there is no significant differential expression of SDH genes between paraganglioma subgroups. This is in agreement with prior observations of SDHB suppression and enhanced expression of SDHA in sporadic, SDHD – and PGL2 -associated tumors[9, 37]. Of further interest is the observed similar gene expression between all paraganglioma subgroups for HIF1α and HIF1α downstream target genes (fig. 2). HIF1α and HIF1α downstream target genes have been shown to be upregulated in SDH-linked tumors[9, 38–40]. The mechanism of HIF1α induction in tumors with SDH mutations has recently been shown to be succinate accumulation resulting from loss of SDH function, leading to inhibition of HIF-α-prolyl hydroxylases and thus to elevated HIF1α activity[39, 41]. The transcription factor HIF1α regulates a host of genes that are involved in proliferation and survival, angiogenesis and glucose metabolism, and the elevated HIF1α activity or pseudo-hypoxic drive is thought to be the basic mechanism of tumorigenesis in SDH-linked paragangliomas[39, 42, 43]. It has been demonstrated that in PGL2 – linked tumors SDH function is disrupted, as it is in SDHD -linked paragangliomas. PGL2 – and SDHD -linked tumors also appear to share the features of increased HIF1α activity and upregulation of HIF1α targets that results from SDH inactivity[9, 37, 41]. These findings may hold important clues for the function of the yet unidentified PGL2 gene on 11q13, as a defect in the yet unidentified PGL2 gene seems to have consequences similar to a mutation in the SDHD gene. No mitochondrial complex II genes are known to be located in the 11q13 region, but the PGL2 gene could affect SDH function by interfering with SDH assembly, transport or insertion into the mitochondrial membrane, or encode a cofactor that is essential for proper SDH function. Alternatively, PGL2 gene function could be more directly associated with HIF1α stability and thus constitute the pseudo-hypoxic drive that leads to paraganglioma formation. We did not find significant differences in expression between paraganglioma subgroups for the PGL2 minimal haplotype on 11q13, and further research to clarify the exact PGL2 identity is currently ongoing.
Another important clinical feature shared by both SDHD -and PGL2 -linked tumors is the remarkable parent-of origin dependent inheritance of disease. Inheritance of paraganglioma occurs in an autosomal dominant way only when paternally transmitted, while no phenotype develops after maternal transmission[44, 45]. Previously, we demonstrated that in SDHD -linked head and neck paragangliomas and phaeochromocytomas this exclusive paternal transmission of the disease is caused by consistent loss of the entire maternal chromosome 11. We hypothesized that selective loss of an as yet unidentified, imprinted gene on the 11p15 region drives this selective chromosome loss, and may also be important in the formation of non- SDHD linked paraganglioma. In line with this hypothesis, recently H19, a paternally imprinted gene on 11p15, has been put forward as the tumor suppressor gene responsible for the parent-of-origin dependent inheritance in SDHD -linked head and neck paragangliomas. In the present study, supervised analysis of all chromosome 11 probe sets on the array, as well as more detailed analysis of genes on chromosome 11p15, 11q23 (location of the SDHD gene) and 11q13 (location of the PGL2 locus), did not show significant expression differences between sporadic, PGL2 – and SDHD -linked tumors (fig. 3). It is possible that this result reflects the loss of chromosome 11 in all these paraganglioma subgroups. As the relation between chromosome loss and gene expression alterations is complex, we must interpret the observed lack of gene expression differences between these groups cautiously in this context. It has been shown previously that all SDHD-linked HN-paragangliomas show loss of the entire copy of the wildtype maternal chromosome 11, and the same applies to PGL2 -linked paragangliomas. Partial or entire chromosome 11 loss has also been observed in sporadic paragangliomas, although only in 2 out of 9 cases. Chromosome 11 loss could thus be an important step in paraganglioma formation irrespective of the genetic background.
In this study of sporadic, SDHD – and PGL2 -linked paragangliomas of the head and neck, we have found very similar gene-expression profiles for all three genetic subgroups. This correlates well with observations of comparable histopathology and clinical behaviour. More detailed analysis of gene sets that have previously been shown to characterize SDH linked tumors, as well as pathways known to be implicated in SDH linked paraganglioma formation, show no differential gene-expression for these paraganglioma subgroups. This suggests that a defect in the yet unidentified PGL2 gene, like a mutation in the SDHD gene, disrupts normal SDH function. Further gene-expression analysis of the PGL2 locus on 11q13 in this study did not reveal the PGL2 identity. The lack of differential gene-expression of chromosome 11 genes between the paraganglioma subgroups might further indicate that chromosome 11 loss, as demonstrated in SDHD-linked paragangliomas, is an important feature in the formation of a paraganglioma regardless of the genetic background.
Baysal BE, Ferrell RE, Willett-Brozick JE, Lawrence EC, Myssiorek D, Bosch A, et al: Mutations in SDHD, a mitochondrial complex II gene, in hereditary paraganglioma. Science. 2000, 287: 848-851. 10.1126/science.287.5454.848.
Astuti D, Latif F, Dallol A, Dahia PLM, Douglas F, George E, et al: Gene mutations in the succinate dehydrogenase subunit SDHB cause susceptibility to familial pheochromocytoma and to familial paraganglioma. Am J Hum Genet. 2001, 69: 49-54. 10.1086/321282.
Niemann S, Muller U: Mutations in SDHC cause autosomal dominant paraganglioma, type 3. Nat Genet. 2000, 26: 268-270. 10.1038/81551.
Taschner PE, Jansen JC, Baysal BE, Bosch A, Rosenberg EH, Brocker-Vriends AH, et al: Nearly all hereditary paragangliomas in the Netherlands are caused by two founder mutations in the SDHD gene. Genes Chromosomes Cancer. 2001, 31: 274-281. 10.1002/gcc.1144.
Mariman EC, Van Beersum SE, Cremers CW, Van Baars FM, Ropers HH: Analysis of a 2nd Family with Hereditary Nonchromaffin Paragangliomas Locates the Underlying Gene at the Proximal Region of Chromosome-11Q. Hum Genet. 1993, 91: 357-361. 10.1007/BF00217356.
Van Houtum WH, Corssmit EP, Douwes Dekker PB, Jansen JC, van der Mey AG, Bröcker-Vriends AH, et al: Increased prevalence of catecholamine excess and phaeochromocytomas in a well-defined Dutch population with SDHD-linked head and neck paragangliomas. Eur J Endocrinol. 2005, 52: 87-94. 10.1530/eje.1.01833.
Astuti D, Latif F, Dallol A, Dahia PLM, Douglas F, George E, et al: Gene Mutations in the Succinate Dehydrogenase Subunit SDHB Cause Susceptibility to Familial Pheochromocytoma and to Familial Paraganglioma. Am J Hum Genet. 2001, 69: 49-54. 10.1086/321282.
Astuti D, Hart-Holden N, Latif F, Lalloo F, Black GC, Lim C, et al: Genetic analysis of mitochondrial complex II subunits SDHD, SDHB and SDHC in paraganglioma and phaeochromocytoma susceptibility. Clin Endocrinol (Oxf). 2003, 59 (6): 728-33. 10.1046/j.1365-2265.2003.01914.x.
Dahia PL, Ross KN, Wright ME, Hayashida CY, Santagata S, Barontini M, et al: A HIF1 alpha regulatory loop links hypoxia and mitochondrial signals in pheochromocytomas. Plos Genet. 2005, 1: 72-80. 10.1371/journal.pgen.0010008.
Hensen EF, Jordanova ES, van Minderhout IJ, Hogendoorn PC, Taschner PE, Mey van der AG, et al: Somatic loss of maternal chromosome 11 causes parent-of-origin-dependent inheritance in SDHD-linked paraganglioma and phaeochromocytoma families. Oncogene. 2004, 23: 4076-4083. 10.1038/sj.onc.1207591.
Bayley JP, van Minderhout I, Weiss MM, Jansen JC, Oomen PHN, Menko FH, et al: Mutation analysis of SDHB and SDHC: novel germline mutations in sporadic head and neck paraganglioma and familial paraganglioma and/or pheochromocytoma. BMC Med Genet. 2006, 7: 1-10.1186/1471-2350-7-1.
Bayley JP, Grimbergen AE, Van Bunderen PA, Van der Wielen M, Kunst HP, Lenders JW, et al: The first Dutch SDHB founder deletion in paraganglioma – pheochromocytoma patients. BMC Med Genet. 2009, 10: 34-10.1186/1471-2350-10-34.
MRC Holland website. [http://www.mrc-holland.com]
Affymetrix website – manuals. [http://www.affymetrix.com/support/technical/manuals.affx]
NCBI Gene Expression Omnibus (GEO) website. [http://www.ncbi.nlm.nih.gov/geo]
Pounds S, Cheng C: Sample size determination for the false discovery rate. Bioinformatics. 2005, 21: 4263-4271. 10.1093/bioinformatics/bti699.
Wu ZJ, Irizarry RA, Gentleman R, Martinez-Murillo F, Spencer F: A model-based background adjustment for oligonucleotide expression arrays. J Am Stat Assoc. 2004, 99: 909-917. 10.1198/016214504000000683.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80. Epub 2004 Sep 15
Bioconductor website. [http://www.bioconductor.org]
Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3.
Goeman JJ, Geer van de SA, de Kort F, van Houwelingen HC: A global test for groups of genes: testing association with a clinical outcome. Bioinformatics. 2004, 20: 93-99. 10.1093/bioinformatics/btg382.
GSEA website of the Broad institute. [http://www.broad.mit.edu/gsea]
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al: PGC-1 alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34: 267-273. 10.1038/ng1180.
Goeman JJ, Buhlmann P: Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007, 23: 980-987. 10.1093/bioinformatics/btm051.
Manoli T, Gretz N, Grone HJ, Kenzelmann M, Eils R, Brors B: Group testing for pathway analysis improves comparability of different microarray datasets. Bioinformatics. 2006, 22: 2500-2506. 10.1093/bioinformatics/btl424.
Goeman JJ, Geer van de SA, van Houwelingen HC: Testing against a high dimensional alternative. J Royal Stat Soc B Stat Meth. 2006, 68: 477-493. 10.1111/j.1467-9868.2006.00551.x.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research. 1999, 27: 29-34. 10.1093/nar/27.1.29.
KEGG website: Kyoto Encyclopedia of Genes and Genomes. [http://www.genome.jp/kegg]
Biocarta website. [http://www.biocarta.com]
Benjamini Y, Hochberg Y: Controlling the False Discovery Rate – A Practical and Powerful Approach to Multiple Testing. J Royal Stat Soc B. 1995, 57: 289-300.
Vanbaars FM, Cremers CWRJ, Vandenbroek P, Veldman JE: Familiar Non-Chromaffinic Paragangliomas (Glomus Tumors) – Clinical and Genetic-Aspects (abridged). Acta Otolaryngol. 1981, 91: 589-593. 10.3109/00016488109138545.
Vandermey AGL, Frijns JHM, Cornelisse CJ, Brons EN, Vandulken H, Terpstra HL, et al: Does Intervention Improve the Natural Course of Glomus Tumors – A Series of 108 Patients Seen in A 32-Year Period. Ann Otol Rhinol Laryngol. 1992, 101: 635-642.
Benn DE, Gimenez-Roqueplo AP, Reilly JR, Bertherat J, Burgess J, Byth K, et al: Clinical presentation and penetrance of pheochromocytoma/paraganglioma syndromes. J Clin Endocrinol Metab. 2006, 91: 827-836. 10.1210/jc.2005-1862.
Neumann HP, Pawlu C, Peczkowska M, Bausch B, McWhinney SR, Muresan M, et al: Distinct clinical features of paraganglioma syndromes associated with SDHB and SDHD gene mutations. JAMA. 2004, 292: 943-951. 10.1001/jama.292.8.943.
Lack EE, Cubilla AL, Woodruff JM: Paragangliomas of the Head and Neck Region – Pathologic-Study of Tumors from 71 Patients. Hum Pathol. 1979, 10: 191-218. 10.1016/S0046-8177(79)80008-8.
Dekker PB, Corver WE, Hogendoom PC, Mey van der AG, Cornelisse CJ: Multiparameter DNA flow-sorting demonstrates diploidy and SDHD wild-type gene retention in the sustentacular cell compartment of head and neck paragangliomas: chief cells are the only neoplastic component. J Pathol. 2004, 202: 456-462. 10.1002/path.1535.
Dekker PB, Hogendoorn PC, Kuipers-Dijkshoorn N, Prins FA, van Duinen SG, Taschner PE, et al: SDHD mutations in head and neck paragangliomas result in destabilization of complex II in the mitochondrial respiratory chain with loss of enzymatic activity and abnormal mitochondrial morphology. J Pathol. 2003, 201: 480-486. 10.1002/path.1461.
Pollard PJ, El Bahrawy M, Poulsom R, Elia G, Killick P, Kelly G, et al: Expression of HIF-1 alpha, HIF-2 alpha (EPAS1), and their target genes in paraganglioma and pheochromocytoma with VHL and SDH mutations. J Clin Endocrinol Metab. 2006, 91: 4593-4598. 10.1210/jc.2006-0920.
Pollard PJ, Briere JJ, Alam NA, Barwell J, Barclay E, Wortham NC, et al: Accumulation of Krebs cycle intermediates and over-expression of HIF1 alpha in tumours which result from germline FH and SDH mutations. Hum Mol Genet. 2005, 14: 2231-2239. 10.1093/hmg/ddi227.
Gimenez-Roqueplo AP, Favier J, Rustin P, Mourad JJ, Plouin PF, Corvol P, et al: The R22X mutation of the SDHD gene in hereditary paraganglioma abolishes the enzymatic activity of complex II in the mitochondrial respiratory chain and activates the hypoxia pathway. Am J Hum Genet. 2001, 69: 1186-1197. 10.1086/324413.
Selak MA, Armour SM, MacKenzie ED, Boulahbel H, Watson DG, Mansfield KD, et al: Succinate links TCA cycle dysfunction to oncogenesis by inhibiting HIF-alpha prolyl hydroxylase. Cancer Cell. 2005, 7: 77-85. 10.1016/j.ccr.2004.11.022.
Hirota K, Semenza GL: Regulation of angiogenesis by hypoxia-inducible factor 1. Crit Rev Oncol Hematol. 2006, 59: 15-26. 10.1016/j.critrevonc.2005.12.003.
Choi KS, Bae MK, Jeong JW, Moon HE, Kim KW: Hypoxia-induced angiogenesis during carcinogenesis. J Biochem Mol Biol. 2003, 36: 120-127.
Struycken PM, Cremers CW, Mariman EC, Joosten FB, Bleker RJ: Glomus tumours and genomic imprinting: Influence of inheritance along the paternal or maternal line. Clin Otolaryngol. 1997, 22: 71-76. 10.1046/j.1365-2273.1997.00884.x.
van der mey AG, Maaswinkelmooy PD, Cornelisse CJ, Schmidt PH, van de kamp JJ: Genomic Imprinting in Hereditary Glomus Tumors – Evidence for New Genetic Theory. Lancet. 1989, 2: 1291-1294. 10.1016/S0140-6736(89)91908-9.
Pigny P, Vincent A, Cardot BC, Bertrand M, de Montpreville VT, Crepin M, et al: Paraganglioma after maternal transmission of a succinate dehydrogenase gene mutation. J Clin Endocrinol Metab. 2008, 93: 1609-1615. 10.1210/jc.2007-1989.
Dannenberg H, de Krijger RR, Zhao JM, Speel EJ, Saremaslani P, Dinjens WN , et al: Differential loss of chromosome 11q in familial and sporadic parasympathetic paragangliomas detected by comparative genomic hybridization. Am J Pathol. 2001, 158: 1937-1942.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/25/prepub
This research was financially supported by a grant from the 'De Drie Lichten' foundation, the Netherlands.
The authors declare that they have no competing interests.
EFH participated in the design of the study, performed the microarray experiments and drafted the manuscript. JJG and JO performed the statistical analysis and helped to draft the manuscript. CWRJC and PCWH helped with designing the study and drafting the manuscript. AGLM coordinated the study and helped to draft the manuscript. CJC conceived of the study, participated in its coordination and helped to draft the manuscript. All authors read and approved the final manuscript.