Eur J Endocrinol
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


DOI: 10.1530/eje.1.01864
European Journal of Endocrinology, Vol 152, Issue 3, 459-470
Copyright © 2005 by European Society of Endocrinology
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Forsberg, L.
Right arrow Articles by Larsson, C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Forsberg, L.
Right arrow Articles by Larsson, C.

EXPERIMENTAL STUDY

Distinction in gene expression profiles demonstrated in parathyroid adenomas by high-density oligoarray technology

Lars Forsberg1,2, Erik Björck1, Jamileh Hashemi1, Jan Zedenius3, Anders Höög4, Lars-Ove Farnebo2, Mark Reimers5 and Catharina Larsson1

1 Department of Molecular Medicine, Karolinska University Hospital CMM L8:01, SE-17176 Stockholm, Sweden, 2 Department of Surgical Sciences, Karolinska University Hospital P9:03, SE-17176 Stockholm, Sweden, 3 Department of Surgery, Center for Metabolism and Endocrinology, Karolinska University Hospital K54, SE-14186 Stockholm, Sweden, 4 Department of Oncology and Pathology, Karolinska University Hospital P1:02, SE-17176 Stockholm, Sweden and 5 Center for Genomics and Bioinformatics, Karolinska Institutet, SE-17177 Stockholm, Sweden

(Correspondence should be addressed to L Forsberg or C Larsson, Department of Molecular Medicine, Karolinska University Hospital CMM L8:01, SE-17176 Stockholm, Sweden; Email: Lars.Forsberg{at}cmm.ki.se or Catharina.Larsson{at}cmm.ki.se)


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Objectives: Somatic deletion of chromosome 11q13 is the most frequent genetic aberration in parathyroid adenoma. To gain further insight into the genetic etiology of parathyroid tumor development, we examined a comprehensive gene expression profile of parathyroid adenomas and normal parathyroid tissues. The results were then evaluated with respect to differences between adenomas and normal parathyroid tissue, and to the presence of loss of heterozygosity (LOH) in chromosomal region 11q13.

Design and methods: Sporadic parathyroid adenomas and normal parathyroids were hybridized against HG-U95Av2 oligonucleotide arrays (Affymetrix) containing a total of 12 625 probe sets. Quantitative real-time PCR (QRT-PCR) was performed in a larger series of parathyroid adenomas, in order to con-firm the microarray results.

Results: Cyclin D1 and c-Jun showed increased expression in adenomas vs normal parathyroids by microarray analysis and QRT-PCR, suggesting an oncogenic role of these genes in parathyroid tumor development. At unsupervised hierarchical clustering, the adenomas fell into two groups: Group I adenomas were characterized by 11q13 LOH, while Group II adenomas lacked this abnormality. In addition, a t-test analysis identified largely non-overlapping genes with differential expression in the tumors subgroups; e.g. in Group I tumors the putative oncogene ENC 1 was found highly over-expressed vs Group II adenomas.

Conclusions: The microarray analyses revealed partly distinctive and partly common expression profiles in parathyroid adenomas with and without 11q13 LOH. In addition, approximately half of the under-expressed genes were mapped to chromosome 11, in agreement with a dose effect following loss of this chromosome.


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Primary hyperparathyroidism (PHPT) is the second most common endocrinopathy following diabetes. In the typical form, a benign adenoma overproduces parathyroid hormone (PTH) leading to hypercalcemia and symptoms following therefrom. Parathyroid tumors are also recognized in a variety of less common conditions, in which one or more glands can be involved and the tumors may be benign or malignant. In addition to the more common sporadic cases, various familial forms have also been recognized, such as multiple endocrine neoplasia type 1 (MEN 1), the hyperparathyroidism-jaw tumor syndrome (HPT-JT) and familial isolated hyperparathyroidism (FIHP). To date, a number of parathyroid tumor-related genes have been identified. For the hereditary disorders, these include germline mutations of the MEN1 gene in MEN 1 and FIHP, of CASR in familial hypocalciuric hypercalcemia variants and of HRPT2 in HPT-JT and FIHP (15).

In sporadic forms of the disease, inactivation of HRPT2 is preferentially found in the malignant cases while MEN1 inactivation is mainly found in the benign form (3, 68). Furthermore, the Cyclin D1 gene as well as several genes related to calcium homeostasis have been suggested to be involved in the development of PHPT (911). Somatic deletion of chromosome 11 is the most frequent genetic aberration in parathyroid adenoma, and has been demonstrated by both comparative genomic hybridization (CGH) and loss of heterozygosity (LOH) techniques (1214). This alteration is seen in a third of the tumors, and generally represents a quite specific chromosomal alteration accompanied by inactivating mutation of the MEN1 gene in the other allele (6, 15). In addition, losses of chromosomes 1p, 6q, 9p, 11p, 13q and 15q as well as gains of chromosomes 7, 16p and 19p are also demonstrated in sporadic hyperparathyroidism (12, 13, 16).

In the present study, we have characterized the global gene expression profiles in a series of sporadic parathyroid adenomas, in an attempt to obtain an improved picture of the genetic etiology behind parathyroid tumor development.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Clinical cases

Seventeen parathyroid tumors (T1–T17) were collected from patients treated surgically for PHPT at the Karolinska University Hospital (Table 1Go). Informed consent was obtained from each patient and the local ethics committee approved the study of the clinical material. All patients had a sporadic form of the disease without a personal or family history of MEN 1, HPT-JT or other familial forms of PHPT. The tumors were classified histopathologically as uniglandular parathyroid adenomas according to WHO criteria (17). Representative sections were analysed to ensure that the samples used for DNA and RNA extractions contained a sufficient proportion of tumor cells for the analysis (i.e. > 70%). Tumors T1–T8 were analysed by microarray, LOH and CGH, while tumors T1–T17 were all analysed by real-time PCR. The two biopsies from normal parathyroid glands (N1 and N2) were obtained from patients treated for nodular thyroid lesions when parathyroid glands had to be excised and reimplanted into the sternomastoid muscle for anatomical reasons. Both these patients had normal levels of PTH, serum calcium and thyroid hormones preoperatively. They were both verified as histopathologically normal, and were used as reference samples in the microarray and real-time PCR analysis.


View this table:
[in this window]
[in a new window]
 
Table 1 Clinical data for the 17 cases with sporadic uniglandular parathyroid adenoma in this study.
 
Isolation of DNA and RNA

All tissue samples were collected immediately after surgery and stored at –70 °C until use. Total RNA was extracted from all 17 tumors (T1–T17) and the two normal parathyroid biopsies (N1 and N2) using the TRIzol reagent (Invitrogen, Carlsbad, CA, USA), followed by purification with an Rneasy mini kit (Qiagen, Hilden, Germany). The RNA purity and quality was confirmed by analysis in an Agilent 2100 Bioanalyser (Agilent Technnologies, Inc, Paolo Alto, CA, USA). For eight of the adenomas (T1–T8), DNA was isolated in parallel with the RNA during the TRIzol extraction. Peripheral blood leukocytes were used as a source of constitutional DNA in the microsatellite analysis (T1–T8).

LOH analyses

The eight matched blood and tumor DNA samples from T1–T8 were genotyped for three microsatellite markers located at the MEN1 locus in 11q13: cen-(D11S4946/MEN1)- D11S4938-D11S4937-tel. For T2 and T3, LOH analyses of markers from the MEN1 region have been previously published (15). The markers were analysed using fluorescent detection in an ABI 377 laser fluorescent sequencer and the GeneScan 3.1 software (Applied Biosystems, Perkin Elmer Corp., Foster City, CA, USA).

Microarray experiments and data analyses

For eight of the tumors (T1–T8) and the two normal biopsies (N1 and N2), cDNA was synthesized from total RNA using a T-7 linked oligodT primer, and then used for cRNA synthesis with biotinylated UTP and CTP. After fragmenting, the labeled RNA was hybridized to HG-U95Av2 oligonucleotide arrays (Affymetrix Inc., Santa Clara, CA, USA) according to the protocol recommended by the supplier. After hybridization, the arrays were washed and stained with streptavidin–phycoerythrin (Molecular Probes, Inc., Eugene, OR, USA) in an Affymetrix fluidics station. The arrays were then scanned in an Affymetrix scanner and the expression values for each probe set (represented by 16–20 selected oligomers) were estimated using the Affymetrix Microarray Suite Software (MAS v5.0). To minimize discrepancies resulting from individual variations regarding sample preparations, hybridization conditions, staining or different oligoarray lots, all probe sets were scaled from each array to a target intensity value of 100 (using Affymetrix MAS v5.0), thus allowing interarray comparisons. Scaling factors (range 2.28–5.26) varied within the acceptable range i.e. <threefold. 3'/5' ratios for GAPDH and ß-actin were confirmed to be within acceptable limits (< threefold), and BioB spike controls were found to be present on all chips. In addition, the internal controls BioC, BioD and CreX were present in increasing intensity.

By using the statistical detection algorithm in Affymetrix MAS v5.0, we were able to separate the transcripts of particular probe sets that are reliably detected (present call), from the transcripts below the threshold of detection and therefore considered as absent. From the 12 626 probe sets on the array, control probes and genes scored as marginal or absent (not detected) were excluded, leaving 4436 present probe sets for the subsequent analyses.

Comparison files were created using the batch analysis tool in the MAS v5.0. These included files for each of the parathyroid adenomas (T1–T8), which were then individually compared with each of the normal parathyroids (N1and N2), thus yielding 16 comparison files. To identify differentially expressed genes between parathyroid adenomas and normal parathyroids, two different algorithms from the MAS v5.0 tool were applied. The change algorithm reported the probe sets to be unchanged, increased, decreased or marginal. The signal log ratio algorithm instead calculated an estimated log value for each probe set in each sample. The algorithms are described at http://www.affymetrix.com/support/technical/technotes/statistical_reference_guide.pdf. For both algorithms the signal values were imported into the GeneSpring 5.0 software tool (Silicon Genetics, Redwood City, CA, USA) to sort out genes with increased or decreased expression. Unsupervised hierarchical clustering analysis of the 4436 present probe sets was performed using Gene-Spring 5.0 software with the minimum distance set to 0.001 and the separation ratio set to 0.95. The comparisons of mean expression levels were compared in the two groups of tumors identified by clustering, i.e. T1–T4 with 11q13 LOH (Group I) and T5–T8 without 11q13 LOH (Group II). The comparisons were carried out using the expression values estimated by dChip analysis (www.dchip.org) (18). This was done to obtain more sensitive estimates of the low-abundance genes than were available through MAS v5.0. To avoid effects related to the wide variation at the upper end of the distribution of values, the logarithmic values were used. t-statistics were computed for each gene between the two groups of tumors identified by clustering. The false positive rate was estimated at several different thresholds for the t-score, by permutating the sample labels and re-computing all gene-wise t-scores. The t-score threshold was set so that the estimated number of false positives was half the number of genes with t-scores exceeding that threshold (i.e. the ‘false discovery rate’ was one half).

CGH

CGH was carried out on tumors T1–T8 using methods previously described (13). For T2 and T3 the results have been previously published (13). A minimum of ten three-color digital images (DAPI, FITC and Spectrum Red fluorescence) were captured from each hybridization using an Axioplan 2 epifluorescence microscope (Carl Zeiss Jena GmbH, Jena, Germany) equipped with a Sensys charged-coupled-device camera (Photometrics, Tucson, AZ, USA) and analysed with the isis/CGH software (Metasystem, Altlussheim, Germany). Relative DNA sequence copy number changes were detected along the length of all chromosomes in each metaphase spread. Green to red ratios < 0.8 were scored as losses, ratios > 1.2 as gains and ratios > 2.0 as high-level amplifications. Heterochromatic regions in the centromeric and paracentromeric parts of some chromosomes, the Y chromosome, the short arm of the acrocentric chromosomes as well as the region next to the telomeres were not included in the evaluation. All tumor samples were hybridized and analysed on two different occasions, with concordant results.

Quantitative real-time PCR (QRT-PCR)

Total RNA was isolated from T1–T17 and N1 and N2 as described above. To avoid the presence of genomic DNA all samples were treated with DNAase after RNA extraction. Ten micrograms of RNA were reverse transcribed into cDNA using the TaqMan reverse transcription Kit (Applied Biosystems). Reverse transcription was performed in a thermal cycler at 25 °C for 10 min and 37 °C for 120 min. For five of the six genes analysed by TaqMan real-time PCR, Assays-on-Demand kits (Applied Biosystems) were available including Cyclin D1 (assay identification (ID) is Hs00233498_m1), c-Jun (Hs00277190_s1), ß-actin (Hs99999903_m1), TBP (Hs99999910_m1) and 18S rRNA (Hs99999901_s1). For each single sample, the master-mix included: 12.5 µl TaqMan 2 x Universal Master-Mix (Applied Biosystems), 8.25 µl RNAse-free water and 1.25 µl 20 x Assays-on-Demand gene expression assay mix. For the 36B4 gene (public ID M17885 [GenBank] ) no Assays-on-Demand are available, while the following previously published primers and probe were used (19): forward primer: 5'-AGATGCAGCAGATCCGCAT-3'; reverse primer: 5'-ATATGAGG AGCAGTTTCTCCAG-3'. The corresponding TaqMan probe 5'-AGGCTGTGGTGC-TGATGGGCAAGAAC-3' was purchased from MedProbe (catalogue number OL-0361-FAMDQ02). For 36B4, the TaqMan PCR Master-Mix was constituted of 12.5 µl TaqMan 2 x Universal Master-Mix, 7 µl RNAse free water, 1.5 µl primers (10 µM) and 1 µl 36B4 probemix (5 µM;) per sample to be analysed. For each sample and gene, 22 µl aliquots of the MasterMix and 3 µl (~15 ng) cDNA were transferred in duplicates to adjacent wells of a MicroAmp Optical 96-well reaction plate (Applied Biosystems) and covered with Optical Caps (Applied Biosystems). The plate was then placed in the ABI Prism 7700 Sequence Detector (Applied Biosystems) and the run commenced. The conditions for every TaqMan PCR reaction were: 2 min at 50 °C (stage 1), 10 min at 95 °C (stage 2) and 40 cycles of 15 s at 95 °C followed by 60 s at 60 °C (stage 3). When the run was completed, output data were analysed using the SDS software and exported into a Microsoft Excel file for analysis. The 36B4 housekeeping gene was used as an endogenous control to normalize the expression level of target genes. Standard curves for 36B4 and target genes were constructed using 1:10 serial dilutions of control cDNA prepared from the human neuroblastoma cell-line SK-NAS.

In the fold-change analysis, the RNA level was standardized by correcting for the 36B4 level in the same sample. For all experimental samples (T1–T17), the quantity was determined from the standard curve and divided by the quantity of the mean normal quantity (N1/N2). Thus, the normal quantity becomes the 1 x sample, and all other quantities are expressed as a relative n-fold difference.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Expression profiling was performed on eight parathyroid adenomas (T1–T8) and two normal parathyroid tissues (N1 and N2) by hybridization to high-density oligonucleotide arrays. Some of the array results were then verified by QRT-PCR in an extended series of 17 sporadic parathyroid adenomas (T1–T17). The clinical details of the 17 cases with PHPT are listed in Table 1Go. The findings and expression analyses are illustrated in Fig. 1Go and Tables 2Go–5Go.



View larger version (44K):
[in this window]
[in a new window]
 
Figure 1 (A) Unsupervised hierarchical clustering of the eight parathyroid adenomas (T1–T8) and two normal parathyroids (N1 and N2) analyzed by microarray. All probe sets that were present in at least eight of the ten arrays, i.e. 4436, were used to create the hierarchical tree. The dendogram was based on MAS v5.0 estimates. The hierarchical tree revealed two distinct cluster subgroups, which were found to consist of the four tumors with 11q13 LOH (Group I) and the four with no 11q13 LOH (Group II). All eight adenomas were distinctly separated from the two normal parathyroids (N1 and N2). The highlighted genes showed significantly altered expression in Group I vs Group II adenomas based on dChip estimates (Table 5Go). Red indicates relatively higher and green indicates relatively lower gene expression. (B) Chromosomal location of the 85 genes showing significantly (P < 0.0032) differential mean expression between adenoma Groups I and II. The t-tests were performed by comparing the mean expression levels between Group I and Group II using the dChip software. Genes detected with increased expression in Group I compared with Group II are shown by green circles to the right of the chromosomes, whereas those with decreased expression are shown by red circles to the left. (C and D) QRT-PCR analysis of the RNA expression levels of (C) Cyclin D1 and (D) C-Jun in normal parathyroid (Norm) and parathyroid adenomas (T1–T17). The results presented represent the mean of duplicate PCR assays and after normalization against the housekeeping gene 36B4. Norm shows the average expression level from the two normal parathyroids analysed (N1 and N2). The numbers above the bars are relative quantities.

 

View this table:
[in this window]
[in a new window]
 
Table 2 Genes with increased expression in parathyroid adenomas (T1–T8) vs normal parathyroids (N1 and N2).
 

View this table:
[in this window]
[in a new window]
 
Table 5 Genes showing significantly differential mean expression levels in parathyroid adenomas with 11q13 LOH (Group I) vs adenomas without 11q13 LOH (Group II).
 
Comparison of parathyroid adenomas and normal parathyroid

A total of 4436 of the 12 625 probe sets represented on the array gave a present call in at least eight of the ten arrays. This set of present genes was considered as informative and used in all comparisons in this study. However, for the analysis of the genes with decreased expression in the adenomas, the 4533 genes detected as present in both normal parathyroids were included in the analysis.

Differentially expressed genes between parathyroid adenomas (T1–T8) and normal parathyroids (N1 and N2) were identified by pair-wise comparison using change and signal log ratio algorithms in the Affymetrix MAS v5.0. The computing results from both analyses were then imported into the GeneSpring software tool, in order to identify probe sets with increased or decreased expression. A total of 224 probe sets showed increased expression in the adenomas vs normal parathyroids using at least one of the algorithms. The 35 probe sets (representing 31 genes) that were independently identified with increased expression in the adenomas using both algorithms are listed in Table 2Go. Several of the over-expressed genes are putative oncogenes involved in cell cycle control or transcriptional regulation, e.g. the Cyclin D1 gene and the c-Jun protooncogene. Similarly, 361 probe sets were detected with decreased expression in the adenomas by at least one algorithm. Eighty of these (representing 75 genes) were identified as decreased by both algorithms (Table 4Go). Noteworthy, CACNB2 and CXCL12 involved in calcium transport and homeostasis were present among these genes with decreased expression.


View this table:
[in this window]
[in a new window]
 
Table 4 Genes with decreased expression in parathyroid adenomas (T1–T8) vs normal parathyroid (N1 and N2).
 
Expression of parathyroid-related genes

In addition, we investigated the microarray expression of genes selected because of their demonstrated or implicated relation to parathyroid functions (Table 3Go). The calcium-sensing receptor gene (CASR) showed a distinct pattern of decreased expression in all eight adenomas. The mean expression level was 0.61 (range 0.43–0.80) as compared with the mean expression of the two normal parathyroids (1.00). Similarly, a reduced expression was noted for the vitamin D receptor gene (VDR), with a mean reduction from 1.0 in the normal parathyroids to 0.47 (range 0.20–0.64) in the eight adenomas.


View this table:
[in this window]
[in a new window]
 
Table 3 Relative expression levels of selected genes with parathyroid-related functions in adenomas as compared with normal parathyroid.
 
Different expression profiles in parathyroid adenomas with and without 11q13 LOH

Similarities in the global gene expression patterns of normal and pathological parathyroids were studied by creating a hierarchical tree. The unsupervised hierarchical clustering analysis of the present probe sets (4436) revealed that the eight tumor samples fell into two distinct cluster groups (Fig. 1AGo). Group I was found to consist of the four tumors with 11q13 LOH (T1–T4) and Group II of the four tumors without 11q13 LOH (T5–T8). The expression profiles between the tumors in the two groups were then compared by direct t-test comparisons of the means between the two groups. The t-tests revealed 85 genes that showed significantly (P < 0.0032) differential mean expression between the two groups (Table 5Go). The results obtained from the two independent approaches (MAS v5.0: Fig. 1Go and dChip: Table 5Go) support the distinction in expression between Group I and Group II adenomas.

As is evident from Table 5Go, approximately 50% of the genes with decreased expression in Group I vs Group II were located on chromosome 11 within and outside the q13 region (Fig. 1BGo). This indicated that the lower expression of some genes in Group I could be a result of low DNA dose, i.e. a reduction in DNA copies. In addition, no gene located on chromosome 11 was detected with increased expression when comparing the adenomas in Group I vs Group II. In order to detect the chromosome 11 status in the adenomas used for array analyses we complemented the LOH analysis with CGH. Of the eight adenomas hybridized onto microrarrays, tumors T1–T4 were by LOH performance identified with a deletion in 11q13 while T5–T8 had retention of the same chromosomal region. CGH results confirmed that the entire chromosome 11 was lost in three (T1–T3) of the four tumors in Group I. In T4, no loss of 11 was detected by CGH, suggesting that tumor T4 carries a balanced alteration or a small deletion of chromsome 11q13 beyond the resolution of CGH. The four adenomas in Group II (T5–T8) all showed retention of chromosome 11 by CGH. In addition, losses were detected in 1p32-ter (T4), 6cen-q22 (T4), 15q (T3, T8), 21q (T3, T8) and 22q (T8), while trisomies were observed for chromosomes 6, 8, 14 and 16 (T7).

QRT-PCR analyses of selected genes

Two genes (Cyclin D1 and c-Jun), detected with increased expression in tumors vs normal parathyroid by the array analysis, were selected for QRT-PCR on a larger series of sporadic parathyroid tumors, including 17 adenomas T1–T17 (Table 1Go) and two normal parathyroid biopsies (N1 and N2). In these analyses, the quantity of the tumor cDNA was determined from the standard curve and divided by the quantity of the normal parathyroid. The 36B4 housekeeping gene was used as an endogenous control to normalize the expression level of target genes. Cyclin D1 thus showed more than twofold increased expression in 4/17 (24%) of the adenomas compared with the average expression of the two normal parathyroid samples (Fig. 1CGo). Similarly, c-Jun showed more than a twofold increased expression in 6/17 (35%) of the adenomas (Fig. 1DGo).


    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
In approaching the molecular genetic pathways of parathyroid tumor development we have characterized global expression profiles in parathyroid adenomas and normal parathyroid tissues, using the oligonucleotide microarray technology. Although the number of samples studied was limited, several significant results were obtained. We could identify genes that were differentially expressed in tumors vs normal tissue. We also investigated the expression profiles of selected genes, previously implicated in parathyroid-related functions, and observed altered expression profiles for some of them. In addition, we detected different expression patterns in parathyroid adenoma subgroups.

Comparison of parathyroid adenomas and normal parathyroid

The result of over-expressed genes in adenomas vs normal parathyroid identifies several putative oncogenes implicated in cell growth and transcriptional regulation, such as the Cyclin D1 gene and the c-Jun protooncogene (Table 2Go). Cyclin D1 over-expression is in agreement with the protein over-expression previously detected by immunohistochemistry in parathyroid adenomas (20, 21). A pericentric inversion with breakpoints at 11q13 and 11p15 is suggested to bring the Cyclin D1 oncogene under the influence of the 5'-regulatory sequences of the PTH gene (9), which results in dramatic over-expression of Cyclin D1. While this rearrangement has been described in a few tumors only, Cyclin D1 protein over-expression has been reported in 18–40% of adenomas by immunohistochemistry studies. Cyclin D1 promotes the G1–S-phase transition, probably by activating a cyclin-dependent kinase, but may also interact with other proteins and influence other pathways related to cell proliferation (22). In our study, over-expression of Cyclin D1 and c-Jun in the adenomas vs normal parathyroid tissue was confirmed by real-time PCR, which represents an accurate method for mRNA expression quantification. In the extended series of 17 adenomas, Cyclin D1 was over-expressed (>twofold) in 24% of the cases, which is in good agreement with the reported immunohistochemistry results (20, 21). Notably, in the four adenomas with Cyclin D1 over-expression, the level was largely increased compared with the normal parathyroid (4- to 36-fold), suggesting its relevance for tumor development.

The c-Jun gene, located in chromosomal region 1p32, has previously been shown to induce oncogenic transformation, probably due to several mechanisms: c-Jun is upregulated in many carcinomas and one suggested function (e.g. in primary hepatocytes and in liver tumors) is to antagonize the proapoptotic activity of p53 (23). Furthermore, c-Jun is a member of the same family of Jun transcription factors as JunD (24). Menin, encoded by the MEN1 gene in 11q13, interacts with the AP-1 transcription factor leading to a repression of the JunD-mediated transcription (25). Interestingly, menin was recently shown to uncouple JunD and c-Jun phosporylation from MAP kinase activation (26). Given the role of MEN1 in parathyroid tumor development, it is tempting to speculate that c-Jun could also play a role in this tumor type. In our array analysis, over-expression of c-Jun was confirmed by two independent probe sets (Table 2Go). Furthermore, real-time PCR showed a >twofold increased expression in 6/17 of the adenomas (35%) as compared with normal parathyroid.

QRT-PCR analyses was performed on two of the most interesting over-expressed genes and correlated well with the array results, thus supporting the accuracy of the array methodology applied. In the real-time PCR, the normalization to a housekeeping gene is currently the most acceptable method to correct for minor variations due to differences in RNA amounts or in efficiencies of reverse transcription. An ideal housekeeping gene should be expressed at an approximately constant level among different tissues and between, for example, tumor and normal samples. We found that the 36B4 housekeeping gene (coding for the acidic ribosomal phosphoprotein) is suitable for studies of normal and adenomatous parathyroid tissues.

Expression of parathyroid-related genes

In the array analyses, Group I adenomas showed lower levels of PTH expression as compared with Group II adenomas. Since the PTH gene is located in 11p15, the lower expression might be explained by a reduced dose effect following loss of a chromosome 11 homologue in Group I adenomas. Two receptors of main interest for parathyroid tumorigenesis are CASR and VDR. Both these receptors showed a distinct pattern of under-expression in all adenomas as compared with the normal parathyroids (Table 3Go). Both CASR and VDR participate in the calcium homeostasis controlled by the parathyroid gland, and it is suggested that the reduced expression of CASR and VDR has a major impact on the altered calcium homeostasis in PHPT. Our array data regarding CASR expression (median 63%, range 43–80%) correlate well with the mRNA in situ data previously published which also showed reduced expression in adenomas vs normal parathyroid (reduced to median 64%, range 41–98%) (10). Furthermore, in our study the VDR expression was reduced in the adenomas (to median 47%, range 20–64%) compared with the expression in normal parathyroid tissue. A similar pattern of VDR under-expression in parathyroid adenomas was also reported by Carling et al. (27). The consistent reduced levels of both VDR and CASR expression are of interest given the recent demonstration that the CASR gene promoters have vitamin D response elements (VDREs) that confer transcriptional responsiveness to 1,25-dihydroxyvitamin D levels (28).

Expression profiles in parathyroid adenomas with and without 11q13 LOH

In order to display potential parathyroid adenoma subgroups, we performed an ‘unsupervised’ clustering analysis, in which no external input was used to guide the analysis process. We used the most common unsupervised analysis method, i.e. hierarchical cluster analysis (29). This method is used to search for patterns without any expectation concerning the number or type of groups that creates the dendogram, as illustrated in Fig. 1AGo. The resulting hierarchical tree demonstrated that the adenomas clustered together in two individual groups. Group I consisted of the four adenomas with 11q13 LOH and Group II of the four adenomas without 11q13 LOH (Fig. 1AGo). In addition, several genes were detected with a significantly (P < 0.0032) distinct expression profile in Group I and II (Table 5Go). Among the genes over-expressed in Group I can be noted ENC1 (5q12–q13.3), which showed 20-fold increase compared with Group II. This may reflect an important downstream effect. ENC1 is over-expressed in primary colon cancer, neuroblastoma and haircell leukemia (30). Functional in vitro assays suggest that the p53 inducible ENC1 is regulated by the ß-catenin/Tcf pathway, and that its over-expression stimulates proliferation and suppresses differentiation. The present series of tumors is too limited in number to draw conclusions about clinical phenotypes and expression profiles. However, it has been suggested that sporadic and MEN 1-associated parathyroid tumors with LOH of the MEN1 locus in 11q13 are larger than tumors without this abnormality (14).

Approximately 50% of the under-expressed genes in Group I are located on chromosome 11 suggesting an underlying gene dose effect (Table 5Go and Fig. 1BGo). This hypothesis was substantiated by the demonstration of loss of the entire chromosome 11 by CGH. The results suggest that adenomas in Group I have partly similar and partly distinguishing expression profiles as compared with the Group II adenomas. The distinct expression profiles would suggest that adenomas with 11q13 LOH partly affect different pathways as compared with the adenomas without LOH. Nevertheless, another set of genes showed similarly differential expression in the entire adenoma group as compared with normal parathyroid. This may result from late downstream effects in common for tumor development of most adenomas, e.g. over-expression of Cyclin D1 and c-Jun. The common alterations could also represent secondary effects of the tumor status or altered calcium homeostasis. The reduced expression observed for, for example CASR and VDR, would preferentially be placed in this last category. PHPT is clearly a heterogenous disease both regarding clinical development and molecular genetic alterations. Although the biological effect of the altered expression pattern is difficult to survey, it can be speculated that partly different molecular pathways could lead to benign tumor development in the parathyroid gland.

In summary, we have shown that genetic subsets in parathyroid adenomas distinguished by the 11q13 LOH status are reflected in the gene expression profiles, but that commonly altered gene expression patterns are also evident. In the future, detailed analysis at the protein level will help us to understand the impact on tumorigenesis of the genes identified in this study.


    Acknowledgements
 
The authors are grateful to Lisa Åhnfalk for excellent technical assistance in handling the tumor samples. The study was financially supported by grants from Torsten and Ragnar Söderberg Foundations, Swedish Cancer Foundation, Stockholm Cancer Society, Gustav V Jubilee Foundation, Vera and Emil Cornell Foundation, Milton Foundation, Stockholm County Council, Fredrik and Ingrid Thuring Foundataion and Wallenberg Consortium North.


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 

    1. Chandrasekharappa SC, Guru SC, Manickam P, Olufemi SE, Collins FS, Emmert-Buck MR, Debelenko LV, Zhuang Z, Lubensky IA, Liotta LA, Crabtree JS, Wang Y, Roe BA, Weisemann J, Boguski MS, Agarwal SK, Kester MB, Kim YS, Heppner C, Dong Q, Spiegel AM, Burns AL & Marx SJ. Positional cloning of the gene for multiple endocrine neoplasia-type 1. Science 1997 276 404–407.[Abstract/Free Full Text]

    2. Teh BT, Esapa CT, Houlston R, Grandell U, Farnebo F, Nordenskjöld M, Pearce CJ, Carmichael D, Larsson C & Harris PE. A family with isolated hyperparathyroidism segregating a missense MEN1 mutation and showing loss of the wild-type alleles in the parathyroid tumors. American Journal of Human Genetics 1998 63 1544–1549.[CrossRef][ISI][Medline]

    3. Carpten JD, Robbins CM, Villablanca A, Forsberg L, Presciuttini S, Bailey-Wilson J, Simonds WF, Gillanders EM, Kennedy AM, Chen JD, Agarwal SK, Sood R, Jones MP, Moses TY, Haven C, Petillo D, Leotlela PD, Harding B, Cameron D, Pannett AA, Höög A, Heath H, James-Newton LA, Robinson B, Zarbo RJ, Cavaco BM, Wassif W, Perrier ND, Rosen IB, Kristoffersson U, Turnpenny PD, Farnebo L-O, Besser GM, Jackson CE, Morreau H, Trent JM, Thakker RV, Marx SJ, Teh BT, Larsson C & Hobbs MR. HRPT2, encoding parafibromin, is mutated in hyperparathyroidism-jaw tumor syndrome. Nature Genetics 2002 32 676–680.[CrossRef][ISI][Medline]

    4. Carling T, Szabo E, Bai M, Ridefelt P, Westin G, Gustavsson P, Trivedi S, Hellman P, Brown EM, Dahl N & Rastad J. Familial hypercalcemia and hypercalciuria caused by a novel mutation in the cytoplasmic tail of the calcium receptor. Journal of Clinical Endocrinology and Metabolism 2000 85 2042–2047.[Abstract/Free Full Text]

    5. Villablanca A, Calender A, Forsberg L, Höög A, Cheng JD, Petillo D, Bauters C, Kahnoski K, Ebeling T, Salmela P, Richardson AL, Delbridge L, Meyrier A, Proye C, Carpten JD, Teh BT, Robinson BG & Larsson C. Germline and de novo mutations in the HRPT2 tumour suppressor gene in familial isolated hyperparathyroidism (FIHP). Journal of Medical Genetics 2004 41 e32.[Free Full Text]

    6. Heppner C, Kester MB, Agarwal SK, Debelenko LV, Emmert-Buck MR, Guru SC, Manickam P, Olufemi SE, Skarulis MC, Doppman JL, Alexander RH, Kim YS, Saggar SK, Lubensky IA, Zhuang Z, Liotta LA, Chandrasekharappa SC, Collins FS, Spiegel AM, Burns AL & Marx SJ. Somatic mutation of the MEN1 gene in parathyroid tumours. Nature Genetics 1997 16 375–378.[CrossRef][ISI][Medline]

    7. Howell VM, Haven CJ, Kahnoski K, Khoo SK, Petillo D, Chen J, Fleuren GJ, Robinson BG, Delbridge LW, Philips J, Nelson AE, Krause U, Hammje K, Dralle H, Hoang-Vu C, Gimm O, Marsh DJ, Morreau H & Teh BT. HRPT2 mutations are associated with malignancy in sporadic parathyroid tumours. Journal of Medical Genetics 2003 40 657–663.[Abstract/Free Full Text]

    8. Shattuck TM, Välimäki S, Obara T, Gaz RD, Clark OH, Shoback D, Wierman ME, Tojo K, Robbins CM, Carpten JD, Farnebo L-O, Larsson C & Arnold A. Somatic and germ-line mutations of the HRPT2 gene in sporadic parathyroid carcinoma. New England Journal of Medicine 2003 349 1722–1729.[Abstract/Free Full Text]

    9. Arnold A, Kim HG, Gaz RD, Eddy RL, Fukushima Y, Byers MG, Shows TB & Kronenberg HM. Molecular cloning and chromosomal mapping of DNA rearranged with the parathyroid hormone gene in a parathyroid adenoma. Journal of Clinical Investigation 1989 83 2034–2040.

    10. Farnebo F, Enberg U, Grimelius L, Bäckdahl M, Schalling M, Larsson C & Farnebo L-O. Tumor-specific decreased expression of calcium sensing receptor messenger ribonucleic acid in sporadic primary hyperparathyroidism. Journal of Clinical Endocrinology and Metabolism 1997 82 3481–3486.[Abstract/Free Full Text]

    11. Carling T, Kindmark A, Hellman P, Lundgren E, Ljunghall S, Rastad J, Åkerström G & Melhus H. Vitamin D receptor genotypes in primary hyperparathyroidism. Nature Medicine 1995 1 1309–1311.[CrossRef][ISI][Medline]

    12. Palanisamy N, Imanishi Y, Rao PH, Tahara H, Chaganti RS & Arnold A. Novel chromosomal abnormalities identified by comparative genomic hybridization in parathyroid adenomas. Journal of Clinical Endocrinology and Metabolism 1998 83 1766–1770.[Abstract/Free Full Text]

    13. Farnebo F, Kytölä S, Teh BT, Dwight T, Wong FK, Höög A, Elvius M, Wassif WS, Thompson NW, Farnebo L-O, Sandelin K & Larsson C. Alternative genetic pathways in parathyroid tumorigenesis. Journal of Clinical Endocrinology and Metabolism 1999 84 3775–3780.[Abstract/Free Full Text]

    14. Friedman E, De Marco L, Gejman PV, Norton JA, Bale AE, Aurbach GD, Spiegel AM & Marx SJ. Allelic loss from chromosome 11 in parathyroid tumors. Cancer Research 1992 52 6804–6809.[Abstract/Free Full Text]

    15. Farnebo F, Teh BT, Kytölä S, Svensson A, Phelan C, Sandelin K, Thompson NW, Höög A, Weber G, Farnebo L-O & Larsson C. Alterations of the MEN1 gene in sporadic parathyroid tumors. Journal of Clinical Endocrinology and Metabolism 1998 83 2627–2630.[Abstract/Free Full Text]

    16. Tahara H, Smith AP, Gas RD, Cryns VL & Arnold A. Genomic localization of novel candidate tumor suppressor gene loci in human parathyroid adenomas. Cancer Research 1996 56 599–605.[Abstract/Free Full Text]

    17. Solcia E, Klöppel G & Sobin LH. Histological classification of tumours of the parathyroid gland. In Histological Typing of Endocrine Tumors, edn 2, pp 48–55, Heidelberg: Springer-Verlag, 2000.

    18. Li C & Wong WH. Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. PNAS 2001 98 31–36.[Abstract/Free Full Text]

    19. Juneja HS, Lee S, Thomazy V, Shipley G & Davies PJ. Acute activation of gp130 gene expression in bone marrow stromal cells by contact with myeloma-derived lymphoblastic cell line ARH77 cell membranes. Journal of Interferon and Cytokine Research 2001 21 157–166.[ISI][Medline]

    20. Hemmer S, Wasenius VM, Haglund C, Zhu Y, Knuutila S, Franssila K & Joensuu H. Deletion of 11q23 and cyclin D1 over-expression are frequent aberrations in parathyroid adenomas. American Journal of Pathology 2001 158 1355–1362.[Abstract/Free Full Text]

    21. Hsi ED, Zukerberg LR, Yang WI & Arnold A. Cyclin D1/PRAD1 expression in parathyroid adenomas: an immunohistochemical study. Journal of Clinical Endocrinology and Metabolism 1996 81 1736–1739.[Abstract]

    22. Arnold A. Genetic basis of endocrine disease 5. Molecular genetics of parathyroid gland neoplasia. Journal of Clinical Endocrinology and Metabolism 1993 77 1108–1112.[CrossRef][ISI][Medline]

    23. Eferl R, Ricci R, Kenner L, Zenz R, David JP, Rath M & Wagner EF. Liver tumor development. c-Jun antagonizes the proapoptotic activity of p53. Cell 2003 112 181–192.[CrossRef][ISI][Medline]

    24. Mechta-Grigoriou F, Gerald D & Yaniv M. The mammalian Jun proteins: redundancy and specificity. Oncogene 2001 20 2378–2389.[CrossRef][ISI][Medline]

    25. Agarwal SK, Guru SC, Heppner C, Erdos MR, Collins RM, Park SY, Saggar S, Chandrasekharappa SC, Collins FS, Spiegel AM, Marx SJ & Burns AL. Menin interacts with the AP1 transcription factor JunD and represses JunD-activated transcription. Cell 1999 96 143–152.[CrossRef][ISI][Medline]

    26. Gallo A, Cuozzo C, Esposito I, Maggiolini M, Bonofiglio D, Vivacqua A, Garramone M, Weiss C, Bohmann D & Musti AM. Menin uncouples Elk-1, JunD and c-Jun phosphorylation from MAP kinase activation. Oncogene 2002 21 6434–6445.[CrossRef][ISI][Medline]

    27. Carling T, Rastad J, Szabo E, Westin G & Åkerström G. Reduced parathyroid vitamin D receptor messenger ribonucleic acid levels in primary and secondary hyperparathyroidism. Journal of Clinical Endocrinology and Metabolism 2000 85 2000–2003.[Abstract/Free Full Text]

    28. Canaff L & Hendy GN. Human calcium-sensing receptor gene. Vitamin D response elements in promoters P1 and P2 confer transcriptional responsiveness to 1,25-dihydroxyvitamin D. Journal of Biological Chemistry 2002; 277: 30337–30350. Epub 32002 May 30329.[Abstract/Free Full Text]

    29. Eisen MB, Spellman PT, Brown PO & Botstein D. Cluster analysis and display of genome-wide expression patterns. PNAS 1998 95 14863–14868.[Abstract/Free Full Text]

    30. Fujita M, Furukawa Y, Tsunoda T, Tanaka T, Ogawa M & Nakamura Y. Up-regulation of the ectodermal-neural cortex 1 (ENC1) gene, a downstream target of the beta-catenin/T-cell factor complex, in colorectal carcinomas. Cancer Research 2001 61 7722–7726.[Abstract/Free Full Text]


Received 25 August 2004
Accepted 14 December 2004




This article has been cited by other articles:


Home page
Am. J. Pathol.Home page
B. Setzer, D. Lebrecht, and U. A. Walker
Pyrimidine Nucleoside Depletion Sensitizes to the Mitochondrial Hepatotoxicity of the Reverse Transcriptase Inhibitor Stavudine
Am. J. Pathol., March 1, 2008; 172(3): 681 - 690.
[Abstract] [Full Text] [PDF]


Home page
Endocr Relat CancerHome page
C Juhlin, C Larsson, T Yakoleva, I Leibiger, B Leibiger, A Alimov, G Weber, A Hoog, and A Villablanca
Loss of parafibromin expression in a subset of parathyroid adenomas.
Endocr. Relat. Cancer, June 1, 2006; 13(2): 509 - 523.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Forsberg, L.
Right arrow Articles by Larsson, C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Forsberg, L.
Right arrow Articles by Larsson, C.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS