help button home button Endocrine Society JCEM
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Journal of Clinical Endocrinology & Metabolism , doi:10.1210/jc.2004-1337
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Jacques, C.
Right arrow Articles by Reynier, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jacques, C.
Right arrow Articles by Reynier, P.
Related Collections
Right arrow Thyroid
Right arrow Endocrine Oncology
The Journal of Clinical Endocrinology & Metabolism Vol. 90, No. 4 2314-2320
Copyright © 2005 by The Endocrine Society

Two-Step Differential Expression Analysis Reveals a New Set of Genes Involved in Thyroid Oncocytic Tumors

Caroline Jacques, Olivier Baris, Delphine Prunier-Mirebeau, Frédérique Savagner, Patrice Rodien, Vincent Rohmer, Brigitte Franc, Serge Guyetant, Yves Malthiery and Pascal Reynier

Institut National de la Santé et de la Recherche Médicale, E0018, Laboratoire de Biochimie et Biologie Moléculaire, Centre Hospitalier Universitaire (C.J., O.B., D.P.-M., F.S., P.R., V.R., Y.M., P.R.), F-49033 Angers, France; Laboratoire d’Anatomie Pathologique, Hôpital Ambroise Paré (B.F.), F-92104 Boulogne, France; and Laboratoire d’Anatomie Pathologique, Centre Hospitalier Régional Universitaire (S.G.), F-37044 Tours, France

Address all correspondence and requests for reprints to: Dr. Caroline Jacques, Institut National de la Santé et de la Recherche Médicale, E0018, Laboratoire de Biochimie et Biologie Moléculaire, Centre Hospitalier Universitaire d’Angers, 4 rue Larrey, F-49033 Angers, France. E-mail: address: caroline.jacques{at}med.univ-angers.fr.


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Thyroid oncocytic adenomas are a class of tumors characterized by the presence of abundant mitochondria. We performed a differential display RT-PCR analysis on two oncocytic adenomas and their paired controls. We then carried out a microarray analysis using the 460 selected, differentially expressed clones on four other oncocytomas and their paired controls. Thirty genes, 12 encoded by mitochondrial DNA and 18 nuclear-encoded, were overexpressed by a factor of at least 2 in the tumors compared with the controls. Seven of the 18 nuclear-encoded genes are involved in protein metabolism: DKFZP434I116, B3GTL, SNX19, RP42, SENP1, UBE2D3, and the CTSB gene, which is known to be particularly deregulated in most thyroid tumors. Other genes are implicated in signal transduction (ITGAV) or tumorigenesis (AF1q). Immunohistochemistry allowed us to confirm overexpression of the ITGAV and CTSB genes at the protein level and showed a marked relocation of the CTSB protein. We confirmed the overexpression of the AF1q oncogene in 56% of 18 oncocytic tumors by quantitative RT-PCR analysis, which attested to the heterogeneity of these tumors. Our results show an increased expression of genes involved in protein metabolism in oncocytoma, the significance of which requires investigation.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
ONCOCYTIC TUMORS, WHICH may develop in various tissues of the human body, are frequently located in the thyroid. Thyroid oncocytomas, which occur in the sporadic or familial form (1), are characterized by a specific, granular cytoplasmic eosinophilia, reflecting the abundance of mitochondria in most of the cells. The proliferation of mitochondria is related to the increased activity of the respiratory chain complexes in the tumors compared with normal controls (2). Until now, in the normal thyroid, no specific behavior of the mitochondria has been described compared with other tissues (3). The causes of the greater frequency of oncocytomas in the thyroid than in other tissues and the link between mitochondrial proliferation in oncocytes and tumorigenesis remain unknown.

A feedback mechanism of mitochondrial proliferation in response to some mitochondrial defect, similar to that observed in various pathologies involving mitochondrial disorders (4), has been postulated for thyroid oncocytoma (5, 6). Recently, a compensatory mechanism was proposed to explain mitochondrial proliferation in renal oncocytoma (7). A defect in the enzymatic activity and protein content of complex I was specifically found in renal oncocytoma as well as in the vicinity of the tumor. In thyroid oncocytomas, no defects of oxidative phosphorylation or any significant mutations in mitochondrial DNA (mtDNA) have been discovered to date (5, 8). However, we reported a deficit in ATP production in thyroid oncocytic tumors as well as in a thyroid oncocytic cellular model (XTC.UC1) (5, 9). The slight uncoupling was associated with the overexpression of uncoupling protein 2 (UCP2). Transcriptional microarray profiling demonstrated the coordinated up-regulation of mitochondrial respiratory chain genes as well as the coordinated regulation of most of the metabolic pathways involved in oxidative metabolism, e.g. overexpression of the genes involved in cytoplasmic glycolysis, the tricarboxylic cycle, and the respiratory chain, and down-regulation of lactate dehydrogenase A. These results suggest that oncocytomas, in contrast with most solid tumors, preferentially use an oxidative metabolism to produce ATP (6). We also found that a ubiquitous pathway of factors involved in the coordination of mitochondrial biogenesis was strongly up-regulated in thyroid oncocytoma (10).

The differential display, an RT-PCR-derived technique developed by Liang and Pardee (11), offers the advantage of revealing the differential expression of identified and unidentified genes. This method, based on the use of short arbitrary primers, allows the selection of reproducible, differentially expressed fragments in two sets of samples, controls and tests. Overexpressed, underexpressed, and isoexpressed fragments can all be visualized in a single experiment. However, although the differential display is considered to be a powerful technique (12), it may produce up to 80% false-positive results (13). Thus, differential expression of the fragments needs to be confirmed using an adequate independent method, such as macroarray analysis.

In the case of thyroid oncocytoma, very few gene expression profiles are currently available (6, 14). To reveal new genes specifically deregulated in oncocytoma, we screened about 60% of the entire transcriptome belonging to two thyroid oncocytic adenomas using the differential display RT-PCR technique. We confirmed the systematic up-regulation of 30 genes in these two samples by the macroarray analysis of four other thyroid oncocytic adenomas. Finally, quantitative RT-PCR and immunohistochemistry were used to validate the overexpression of several genes and proteins.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Biological material

Two sporadic tumors, diagnosed, according to the World Health Organization classification (15), as a follicular adenoma with Hurtle cells and an oncocytic tumor of uncertain malignant potential, and two paired normal thyroid tissue samples were chosen to perform the differential display analysis. For the macroarray analysis, we used four other sporadic thyroid tumors, diagnosed as follicular oncocytic adenomas, and four paired normal thyroid tissues. For the quantitative RT-PCR analysis, 18 thyroid tumors, consisting of 13 oncocytic adenomas and five oncocytic carcinomas, were chosen. Normal thyroid samples were taken distant from the tumor and histologically tested as being appropriate for controls. All samples, which were rendered anonymous (i.e. with patient identifiers deleted before the study), were deep-frozen in liquid nitrogen immediately after surgery and conserved at –80 C.

For the immunohistochemistry study, we used 11 formalin-fixed, paraffin-embedded, thyroid oncocytic tumor samples and their normal counterparts: five oncocytic follicular carcinomas and six oncocytic follicular adenomas.

RNA preparation

RNA was isolated from frozen tissues using the guanidinium isothiocyanate procedure (TRIzol reagent, Invitrogen Life Technologies, Inc., Gaithersburg, MD). Quantification, degradation, and DNA contamination of RNA were assessed using an RNA 6000 Nano Assay (Agilent Technologies, Palo Alto, CA) following the manufacturer’s procedure.

Differential display analysis

The differential display analysis was performed using the RNAimage sets G501, G502, and G503 (RNAimage kit, GeneHunter Corp., Nashville, TN) (11) according to the manufacturer’s instructions. Briefly, these three sets allow 72 PCR amplifications per sample, using 24 different arbitrary primers. Statistically, this 72 PCR per sample analysis is estimated to cover between 56 and 71% of the cell transcriptome (see www.genhunter.com/support/index.html). The RTs and the PCR amplifications were performed on a PTC 200-MJ Research Thermocycler (MJ Research, Watertown, MA). The PCR amplifications of oncocytoma and paired normal thyroid samples obtained with the same primer pair were loaded in consecutive wells to visually select differentially expressed bands (Fig. 1Go). The bands were cut and eluted from the acrylamide gels, and reamplified using the original set of primers. The sizes of all the reamplified PCR products were checked on a 2.5% agarose gel with the DNA Molecular Weight Marker IX (Roche, Mannheim, Germany).



View larger version (120K):
[in this window]
[in a new window]
 
FIG. 1. Differential display gel. Examples of selected bands from tumors (T) and paired control thyroid tissues (C) are highlighted: three underexpressed bands (white ovals) and three overexpressed bands (white rectangles). Each four-lane group corresponds to one of the 72 primer pair amplifications used for the two thyroid oncocytic adenomas and control tissues analyzed.

 
Cloning and sequencing

After reamplification, the PCR products were cloned into the pCR 2.1-TOPO vector (TOPO TA Cloning kit, Invitrogen Life Technologies, Carlsbad, CA) following the manufacturer’s protocol. Plasmid DNA was extracted from three overnight subcultures of each cloned fragment, with the QIA Prep Spin Miniprep kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The sequencing reactions were performed on 50-fmol plasmid DNA using the two insert flanking primers, the M13 reverse primer and the M13 forward primer (TOPO TA Cloning kit), and the CEQ DTCS-Quick Start kit (Beckman Coulter, Fullerton, CA) following the manufacturer’s recommendations. The sequencing reactions were analyzed on the capillary automated sequencer CEQ 8000 Genetic Analysis System (Beckman Coulter).

Macroarrays

The plasmids containing the 460 inserts obtained by differential display and two control cDNA clones (UCP2 and ANT2) were PCR-amplified in 96-well plates using the two M13 primers mentioned above. The amplification conditions were 95 C for 10 min, followed by 35 cycles at 95 C for 30 sec, 48 C for 30 sec, and 72 C for 30 sec, with a final step at 72 C for 10 min. Several wells containing DNA-free water were used as PCR controls. The PCR amplification products were spotted on Hybond TM-N+ membranes optimized for nucleic acid transfer (Amersham Biosciences, Piscataway, NJ) using a Microgrid II spotter (Biorobotics, Cambridge, UK). After denaturation (0.5 M NaOH and 1.5 M NaCl) and DNA fixation [2 h at 80 C, and UV cross-linking (UV Stratalinker 1800; Stratagene, La Jolla, CA)], the membranes were kept at room temperature until hybridization. The first hybridization was performed with an oligonucleotide corresponding to a small portion of the cloning plasmid (5'-CCCTATAGTGAGTCGTATTA-3') to quantify the amount of DNA available for hybridization. One microgram of the oligonucleotide was labeled with [{gamma}-33P]ATP using the T4 polynucleotide kinase kit (Forward Labeling Reaction, Invitrogen Life Technologies, Inc., Groningen, The Netherlands) following the manufacturer’s protocol. The labeled oligonucleotide was then purified on MicroSpin G-25 columns (Amersham Biosciences).

The membranes were prehybridized for 4 h at 42 C in 10 ml 5x standard saline citrate (Eurobio, Les Ulis, France), 5x Denhardt’s, 0.5% sodium dodecyl sulfate, and 100 µg/ml sonicated salmon sperm DNA (Roche) buffer. After the addition of 400,000 cpm of the labeled plasmid oligonucleotide to the buffer, hybridization was performed for an additional 2 h at 42 C. After a washing step in a 2x standard saline citrate/0.1% sodium dodecyl sulfate buffer, the membranes were dried and exposed on a Super Resolution Cyclone Storage Phosphor Screen (Packard-PerkinElmer Life Sciences, Boston, MA; Fig. 2AGo). The signals were analyzed on a Cyclone Storage Phosphor System (Packard-Perkin-Elmer Life Sciences) and quantified with the Quantarray software (Packard-PerkinElmer Life Sciences).



View larger version (39K):
[in this window]
[in a new window]
 
FIG. 2. A macroarray set. A, Pattern of the spotted fragments after hybridization with the [{gamma}-33P]ATP-labeled plasmid oligonucleotide, used to correct the cDNA probe signal measurements. B, Pattern after hybridization with [{alpha}-33P]dCTP-labeled total cDNA of the control thyroid of one patient. C, Pattern after hybridization with [{alpha}-33P]dCTP-labeled total cDNA of the tumor of the same patient. Several overexpressed bands are observed in the tumor; three of these are indicated as examples (arrows).

 
After stripping the membranes, hybridizations were performed with the [{alpha}-33P]deoxy (d)-CTP-labeled total cDNA of the thyroid tumors and the normal thyroid samples. Each sample was hybridized on an individual array. Five micrograms of total RNA were used for the RT. After addition of 8 µg deoxythymidine25 oligonucleotide (Promega Corp., Madison, WI), the mix was incubated 8 min at 68 C, and for 30 min while decreasing the temperature from 68 to 42 C. After the addition of 40 U RNasin ribonuclease inhibitor (Promega Corp.), 5x buffer, 6.7 mM dithiothreitol (Superscript II kit, Invitrogen Life Technologies, Inc.), 0.4 mM dATP, dCTP, dGTP, 2.4 µM dCTP, 30 µCi dCTP ([{alpha}-33P]-labeled; PerkinElmer Life Sciences), and 200 U SuperScript II, the incubation was restarted for 1 h at 42 C. At this step, 200 U of the enzyme were added, and the incubation was prolonged for 1 h. The reaction was stopped by the addition of 0.33% sodium dodecyl sulfate, 13 mM EDTA, and 0.3 M NaOH and incubation for 30 min at 68 C, followed by 15 min at room temperature. Before purification on an AutoSeq G-50 column (Amersham Biosciences), 0.25 M Tris and 0.15 mM HCl were added. The purified labeled cDNA was then mixed with 2 µg dA80 oligonucleotide, denatured for 5 min at 100 C, and incubated for 2 h 30 min at 65 C in 1 ml hybridization buffer.

The stripped membranes were prehybridized in 10 ml prehybridization solution (see above) for 6 h at 68 C. The total cDNA-labeled probe was then added, and the incubation was performed for 48 h at 68 C. After the washing steps, the membranes were exposed, and the signals were detected as described above (Fig. 2Go, B and C).

The total cDNA probe signal measurements were first corrected by means of the amount of spotted DNA measured by hybridization of the plasmid oligonucleotide as previously described (16). For each spot, the final measurement corresponded to the ratio of the sample signal to the oligonucleotide signal. Spots exhibiting weak plasmid hybridization (small amounts of spotted DNA) were considered as missing values. Spots with expression similar to the background were withdrawn from analysis. Normalization between arrays was obtained by dividing each gene measurement by the median expression value of the array. For each gene, we then calculated the O/T ratio, defined as the ratio between gene expression in the tumor vs. that in normal thyroid tissue. Only the genes for which the mean cut-off differential expression value for the four oncocytomas vs. normal thyroid samples was confirmed (i.e. O/T ≥ 2) were retained for additional analysis.

In silico analysis

The FASTA format of the insert sequences was compared with the databanks GenBank, dbEST (www.ncbi.nlm.nih.gov/), and UCSC Genome Browser (www.genome.ucsc.edu/cgi-bin/, July 2003 version) for identification. The best matching sequences were considered to be identical (e-value: 0–7 E-36; mean e-value, 3.04 E-37; median, 1E-64). For each unidentified confirmed sequence, we checked a putative gene prediction with a BLAT alignment of the sequence with the UCSC Genome Browser databank sequences. For each of the predicted genes, the corresponding protein sequence was analyzed with the Mitoprot package (http://ihg.gsf.de/ihg/) and the TargetP package (www.cbs.dtu.dk/services/TargetP) to predict a putative N-terminal mitochondrial targeting signal (see Table 2Go).


View this table:
[in this window]
[in a new window]
 
TABLE 2. MitoProt and TargetP values of a putative mitochondrial N-terminal targeting signal of the proteins or potential proteins confirmed differentially expressed genes in four oncocytomas vs. control tissues

 
Quantitative PCR

The expression level for the AF1q gene was explored by quantitative real-time PCR analysis using the actin gene as a reference, as previously described (10). The AF1q standard was obtained using the forward primer 5'-CCATCTTTGGAACACGCCAG-3' and the reverse primer 5'-TTTCTCCTGGTCTGCTGCAG-3'. The actin standard was obtained using the forward primer 5'-CGACATGGAGAAAATCTGGC-3' and the reverse primer 5'-AGGTCCAGACGCAGGATGG-3'.

Immunohistochemistry

The 11 formalin-fixed, paraffin-embedded, thyroid tumor samples (five oncocytic follicular carcinomas and six oncocytic follicular adenomas) and their normal counterparts were used for tissue array construction, as described by Kononen et al. (17). A triplet of tumor tissue and a triplet of the normal counterpart were arrayed for each sample. Immunostaining was performed using the standard avidin-biotin peroxidase technique. The primary antibodies were anticathepsin B (human liver; dilution, 1:400; Calbiochem, San Diego, CA) and antihuman CD51/61 (integrin {alpha}Vß3-vitronectin receptor (ITGAV); clone 23C6; dilution, 1:10; Leinco Technologies, St. Louis, MO). Diaminobenzidine was used as the chromogen, and hematoxylin was used as the nuclear counterstain. For negative controls, the primary antibody was either omitted or replaced by a suitable concentration of normal IgG of the same species.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
We used a two-step differential expression analysis method to study six thyroid oncocytic adenomas and their paired normal tissues. After visualizing approximately 10,000 bands per sample, we selected and cloned 100 underexpressed and 360 overexpressed bands in the tumor (Fig. 1Go). Because the differential display technique is known to produce a high rate of false-positive reactions (up to 80%) (13), we designed a nylon membrane macroarray assay to confirm the differential expression of clones among the 460 bands identified by differential display analysis of the two tumors. We spotted the ANT2 and UCP2 genes to act as positive controls on the membranes. The microarray measurement of mRNA expression was made in four other pairs of thyroid oncocytic adenomas and their normal paired controls (Fig. 2Go). We selected the clones differentially expressed in all four thyroid oncocytic adenomas. This highly selective, two-step approach enabled us to confirm 12 underexpressed and 49 overexpressed clones in the tumor among the 460 differentially expressed clones identified. Among these differentially expressed clones, the O/T ratio ranged from 0.71–6.52. The two control genes spotted on the macroarrays, i.e. the ANT2 and UCP2 genes, were confirmed as being overexpressed, with overexpression factors of 1.98 and 1.77, respectively, according to our previous observations (5, 6). For additional analysis, we chose to focus on the clones that showed at least a 2-fold differential expression compared with controls, although some of the excluded genes could be of interest (for example, the pendrin gene, which was down-regulated by a factor of 0.79). The cut-off value of an O/T of 2 or greater highlighted 30 overexpressed genes.

The 30 clones selected had factors of overexpression in the tumor ranging from 2.06–6.52 compared with normal thyroid tissue. Table 1Go recapitulates the identification information for each of these sequences. Twenty-five of the 30 sequences (83%) corresponded to well-known genes, whereas five did not match any previously identified genes. Twelve of the 30 sequences (40%) were mtDNA-encoded genes (11 of 12) or a mtDNA sequence (one of 12). Ten of these sequences corresponded to genes coding for subunits of the four respiratory complexes. As shown in Table 1Go, many of these genes are among the most up-regulated genes in our study. However, the overexpression factor of the mtDNA-encoded genes in oncocytoma ranges from 2.15–6.52 (mean, 3.15). One of the remaining two clones corresponded to the 16S rRNA gene. The last sequence was identical to part of the RNA primer sequence, located in the mtDNA control region, which is necessary for initiation of mtDNA replication.


View this table:
[in this window]
[in a new window]
 
TABLE 1. Overexpressed genes, with O/T ratio of 2 or more, selected by differential display and confirmed by macroarrays in six thyroid oncocytic adenomas compared to the normal thyroid paired tissues

 
In addition to the 12 mtDNA genes, 18 genes or sequences (60%) were nuclear-encoded and overexpressed by factors varying from 2.06–3.62 (mean, 2.6). Two genes were found to encode transcription factors (ZMYND12 and TCF8). Seven of the 18 genes were involved at multiple levels of protein metabolism: mRNA transport (DKFZP434I116), posttranslational modifications (B3GTL), protein trafficking (SNX19), ubiquitination/proteasome (RP42 protein) (18), and proteolysis (cathepsin B, SENP1, and UBE2D3). One gene was involved in iron metabolism (FTH1). Other genes corresponded to proteins involved in cell adhesion/signal transduction and regulation (ITGAV and TCF8), and one gene was found to encode a cell membrane protein (EPB41L4A) of unknown function.

Figure 3Go shows the results of the immunostaining with the anticathepsin B and the anti-CD51/61 (ITGAV protein) antibodies of control samples and thyroid oncocytomas (six adenomas and five carcinomas). All tumors studied presented an overexpression of the two proteins compared with the control tissues. For the cathepsin B protein, cellular localization of the signal varied from one case to another. However, in the adenomas as well as the carcinomas, the signal was cytoplasmic and frequently concentrated at the apical border in the form of droplets of different sizes. For the integrin {alpha}Vß3 protein (CD51/61), a dense cytoplasmic signal was found in adenomas (Fig. 3Go) as well as in carcinomas. The presence of droplets and a lining at the apical border was rarely observed.



View larger version (120K):
[in this window]
[in a new window]
 
FIG. 3. Immunohistochemistry staining of normal thyroid tissues, oncocytic follicular adenomas, and oncocytic follicular carcinomas for anticathepsin B and antihuman CD51/61 (integrin {alpha}Vß3) antibodies (magnification, x400). Cathepsin B and CD51/61, overexpression in either oncocytic adenomas and carcinomas (brown cytoplasmic staining). For the cathepsin B protein, droplets are often located in the vicinity of the apical border of the cells (black arrows).

 
One sequence of the overexpressed clones corresponded to the AF1q gene (ALL1-fused gene from chromosome 1q) located in 1q21.3. We explored 18 oncocytic tumors (13 adenomas and five carcinomas) and their paired control tissues for the expression level of this gene by quantitative RT-PCR analysis. AF1q was overexpressed in 10 (56%) of the oncocytomas (seven adenomas and three carcinomas) vs. normal thyroid tissue, with a mean overexpression factor of 1.87 ± 0.54. Seven of the 13 adenomas (54%) showed a mean overexpression factor of 1.97 ± 0.6. Both of these results were significant using the Wilcoxon test (P < 0.05). The eight remaining oncocytomas presented isoexpression in three cases (17%) or underexpression in five cases (28%) compared with normal thyroid tissue.

Finally, in silico analysis allowed a gene prediction for three of the five isolated sequences that did not correspond to any known genes. One of these three sequences, 2q33.3, showed 100% identity (with respect to 59 bp) with the 3' extremity of the alternative variant gDec03 of the NDUFS1 gene. (The protein NDUFS1 is a nuclear-encoded component of the mitochondrial respiratory chain complex I.) The second sequence, 7p21.1, showed a 30% identity (with respect to 75 amino acids) with the elongation factor 1{alpha} of Macrosiphoniella ludovicianae involved in cytoplasmic translation (eIF1{alpha}; e-value, 1.6). The gene LOC115294 was found to be similar to the hypothetical protein FLJ0883 involved in posttranslational modification (protein L-isoaspartate O-methyltransferase). These three predicted proteins (2q33.3, 7p21.1, and LOC115294) displayed a high probability of an N-terminal mitochondrial targeting signal by analysis with the Mitoprot and TargetP targeting prediction packages (Table 2Go). We did not find any particular homology concerning the unknown sequences located in 2q12.1 and 14q13.2.


    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
The differential gene expression analysis coupled with the macroarray analysis of six thyroid oncocytic adenomas and corresponding controls was aimed at the identification of genes implicated in mitochondrial biogenesis and thyroid oncocytic tumorigenesis. We identified 30 genes overexpressed at least 2-fold in the thyroid oncocytic tumors. This result is in accordance with a previous microarray study that showed mostly up-regulated genes in oncocytomas (5).

The 30 overexpressed sequences included 10 of the 13 mtDNA-encoded genes corresponding to 13 subunits of four of the five respiratory chain complexes. The mtDNA-encoded genes are known to be up-regulated in the context of oncocytomas, which are characterized by increased mitochondrial biogenesis (5, 6, 19). The up-regulation of ND5, ND2, and UCP2 mRNA has been demonstrated with the RT-PCR method in 22 thyroid oncocytomas compared with the paired control tissues (5). Our results showed that nearly all mtRNA genes are up-regulated. The transcription of mtDNA is polycistronic; each strand of the molecule is entirely transcribed before cleavage of the different rRNA, transfer RNA, and mRNA sequences. The surprising variation in the level of expression of the 10 mRNAs may be explained by differences in posttranscriptional processing and in the half-life of mitochondrial mRNA.

Because the differential display technique uses polyT primer (HT11V) at the RT step, the up-regulation of the mitochondrial 16S rRNA, which is also mtDNA-encoded, was rather unexpected. However, it has been shown that 16S rRNA is polyadenylated by up to 10 residues (20), which is sufficient to allow RT.

We also isolated a fragment that matches the RNA replication primer sequence of the mtDNA control region. The transcription of this region produces an RNA sequence that, after cleavage, provides an RNA primer involved in the initiation of mtDNA replication (21). The overexpression of this sequence is of particular interest, because it is probably a key regulator of the balance between the two tightly coupled processes of replication and transcription of the mitochondrial genome. It has been shown in thyroid oncocytomas that this balance differs from that in normal thyroid tissues; the mRNA level is particularly high considering the mtDNA content (10). Although the RT of this nonpolyadenylated sequence is singular, its effective and reproducible overexpression may be involved in the uncoupling of replication and transcription.

In addition to the mtDNA-encoded genes, 13 known nuclear-encoded genes were also up-regulated in oncocytomas. Seven of these genes are associated with the cellular protein turnover at various levels of protein metabolism. In particular, the cathepsin B gene is the most overexpressed nuclear gene in our differential expression study, and its overexpression was confirmed at the protein level. The cathepsin B gene codes for a lysosomal cysteine endopeptidase located in the apical membrane of thyroid epithelial cells (22). This protein is involved in processing of the thyroglobulin hormone in the normal thyroid and is well known to be up-regulated at the transcriptional, translational, and posttranslational levels in human thyroid tumors (6, 23). It has therefore been suggested that this protein may be a reliable marker of the progression from the premalignant to the malignant state in the development of thyroid tumors (24, 25). The major dislocation of the protein at the apical border of the cells in all the tumors explored suggests an abnormal lysosome location or abnormal protein targeting or processing in oncocytoma.

Genes associated with other functions, such as transcription regulation, cell adhesion, and signal transduction, were also well represented. Among these genes, several have been shown to be involved in tumorigenesis or the cell cycle, e.g. TCF8 (transcription factor 8), ZMYND12 (zinc finger MYND domain containing 12 proteins), and ITGAV (integrin {alpha}V). ITGAV encodes the integrin {alpha}-chain V protein involved in cell adhesion (cell to cell junctions) and signal transduction. Integrins are known to be involved in various cancer processes, such as invasion, angiogenesis, and metastasis (26). The overexpression of genes involved in cell adhesion and proteolysis has been described in thyroid oncocytomas at a similar level of overexpression (6).

Of particular interest was the overexpression of the nuclear-encoded gene AF1q, which codes for a 9-kDa protein of unknown function with no similarity to any other protein. The AF1q gene was previously implicated in a translocation t(1;11)(q21;q23) responsible for cases of leukemia (27). The up-regulation of this gene has been shown to be greater in leukemic cell lines than in normal hemopoietic tissues (28). Moreover, the expression of this gene in pediatric acute myeloid leukemia may correlate with the clinical outcome of the patients (29). Our results, displaying no differences in the overexpression of the gene between adenoma and carcinoma oncocytic thyroid tumors, suggest that overexpression of the AF1q oncogene may be implicated in thyroid tumorigenesis. Finally, we revealed three putative genes encoding potential mitochondrial proteins, in which the functions have to be explored.

The set of 30 genes, overexpressed by a factor of at least 2, detected in thyroid oncocytic tumors reveals a great up-regulation of mtDNA genes and of the genes implicated in protein metabolism (seven known nuclear-encoded genes, two in silico predicted mitochondrial nuclear-encoded genes, and the mtDNA-encoded 16S RNA). These results show that mitochondrial proliferation is clearly associated with intense protein turnover and raise an important question. Is this increased protein turnover related to the oncocytoma tumorigenesis? We could postulate that increased protein turnover is necessary for the intense mitochondrial proliferation or that saturation of the cytoplasm by mitochondria could disturb cell trafficking and processing of the proteins, as suggested by relocation of the cathepsin B protein in tumors.


    Acknowledgments
 
We are grateful to Aurélie Paquereau, Dominique Couturier, Florence Fatih, and Laurence Preisser for their technical support and advice, and to Kanaya Malkani and Nicola Jordan for their critical reading of the manuscript.


    Footnotes
 
This work was supported by grants from the University Hospital of Angers (PHRC 01–10), Institut National de la Santé et de la Recherche Médicale, and the University of Angers.

First Published Online December 28, 2004

Abbreviations: mtDNA, Mutations in mitochondrial DNA; UCP2, uncoupling protein 2.

Received July 9, 2004.

Accepted December 20, 2004.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. Katoh R, Harach HR, Williams ED 1998 Solitary, multiple, and familial oxyphil tumours of the thyroid gland. J Pathol 186:292–299[CrossRef][Medline]
  2. Ebner D, Rodel G, Pavenstaedt I, Haferkamp O 1991 Functional and molecular analysis of mitochondria in thyroid oncocytoma. Virchows Arch B Cell Pathol 60:139–144[Medline]
  3. Lamy FM, Rodesch FR, Dumont JE 1967 Action of thyrotropin on thyroid energetic metabolism. VI. Regulation of mitochondrial respiration. Exp Cell Res 46:518–532[CrossRef][Medline]
  4. Joseph AM, Rungi AA, Robinson BH, Hood DA 2004 Compensatory responses of protein import and transcription factor expression in mitochondrial DNA defects. Am J Physiol 286:C867–C875
  5. Savagner F, Franc B, Guyetant S, Rodien P, Reynier P, Malthiery Y 2001 Defective mitochondrial ATP synthesis in oxyphilic thyroid tumors. J Clin Endocrinol Metab 86:4920–4925[Abstract/Free Full Text]
  6. Baris O, Savagner F, Nasser V, Loriod B, Granjeaud S, Guyetant S, Franc B, Rodien P, Rohmer V, Bertucci F, Birnbaum D, Malthiery Y, Reynier P, Houlgatte R 2004 Transcriptional profiling reveals coordinated up-regulation of oxidative metabolism genes in thyroid oncocytic tumors. J Clin Endocrinol Metab 89:994–1005[Abstract/Free Full Text]
  7. Simonnet H, Demont J, Pfeiffer K, Guenaneche L, Bouvier R, Brandt U, Schagger H, Godinot C 2003 Mitochondrial complex I is deficient in renal oncocytomas. Carcinogenesis 24:1461–1466[Abstract/Free Full Text]
  8. Yeh JJ, Lunetta KL, van Orsouw NJ, Moore Jr FD, Mutter GL, Vijg J, Dahia PL, Eng C 2000 Somatic mitochondrial DNA (mtDNA) mutations in papillary thyroid carcinomas and differential mtDNA sequence variants in cases with thyroid tumours. Oncogene 19:2060–2066[CrossRef][Medline]
  9. Zielke A, Tezelman S, Jossart GH, Wong M, Siperstein AE, Duh QY, Clark OH 1998 Establishment of a highly differentiated thyroid cancer cell line of Hurthle cell origin. Thyroid 8:475–483[Medline]
  10. Savagner F, Mirebeau D, Jacques C, Guyetant S, Morgan C, Franc B, Reynier P, Malthiery Y 2003 PGC-1-related coactivator and targets are upregulated in thyroid oncocytoma. Biochem Biophys Res Commun 310:779–784[CrossRef][Medline]
  11. Liang P, Pardee AB 1992 Differential display of eukaryotic messenger RNA by means of the polymerase chain reaction. Science 257:967–971[Abstract/Free Full Text]
  12. Broude NE 2002 Differential display in the time of microarrays. Expert Rev Mol Diagn 2:209–216[CrossRef][Medline]
  13. Janssen YM, Driscoll KE, Timblin CR, Hassenbein D, Mossman BT 1998 Modulation of mitochondrial gene expression in pulmonary epithelial cells exposed to oxidants. Environ Health Perspect 106(Suppl 5):1191–1195
  14. Mirebeau-Prunier D, Guyetant S, Rodien P, Franc B, Baris O, Rohmer V, Reynier P, Tourmen Y, Malthiery Y, Savagner F 2004 Decreased expression of thyrotropin receptor gene suggests a high-risk subgroup for oncocytic adenoma. Eur J Endocrinol 150:269–276[Abstract]
  15. DeLellis R, Lloyd R, Heitz P, Eng C 2004 World Health Organization classification of tumours: pathology and genetics of tumours of endocrine organs. Geneva: World Health Organization
  16. Bernard K, Auphan N, Granjeaud S, Victorero G, Schmitt-Verhulst AM, Jordan BR, Nguyen C 1996 Multiplex messenger assay: simultaneous, quantitative measurement of expression of many genes in the context of T cell activation. Nucleic Acids Res 24:1435–1442[Abstract/Free Full Text]
  17. Kononen J, Bubendorf L, Kallioniemi A, Barlund M, Schraml P, Leighton S, Torhorst J, Mihatsch MJ, Sauter G, Kallioniemi OP 1998 Tissue microarrays for high-throughput molecular profiling of tumor specimens. Nat Med 4:844–847[CrossRef][Medline]
  18. Mas C, Bourgeois F, Bulfone A, Levacher B, Mugnier C, Simonneau M 2000 Cloning and expression analysis of a novel gene, RP42, mapping to an autism susceptibility locus on 6q16. Genomics 65:70–74[CrossRef][Medline]
  19. Heddi A, Faure-Vigny H, Wallace DC, Stepien G 1996 Coordinate expression of nuclear and mitochondrial genes involved in energy production in carcinoma and oncocytoma. Biochim Biophys Acta 1316:203–209[Medline]
  20. Dubin DT, Montoya J, Timko KD, Attardi G 1982 Sequence analysis and precise mapping of the 3' ends of HeLa cell mitochondrial ribosomal RNAs. J Mol Biol 157:1–19[CrossRef][Medline]
  21. Lee DY, Clayton DA 1996 Properties of a primer RNA-DNA hybrid at the mouse mitochondrial DNA leading-strand origin of replication. J Biol Chem 271:24262–24269[Abstract/Free Full Text]
  22. Linke M, Herzog V, Brix K 2002 Trafficking of lysosomal cathepsin B-green fluorescent protein to the surface of thyroid epithelial cells involves the endosomal/lysosomal compartment. J Cell Sci 115:4877–4889
  23. Srisomsap C, Subhasitanont P, Otto A, Mueller EC, Punyarit P, Wittmann-Liebold B, Svasti J 2002 Detection of cathepsin B up-regulation in neoplastic thyroid tissues by proteomic analysis. Proteomics 2:706–712[CrossRef][Medline]
  24. Shuja S, Murnane MJ 1996 Marked increases in cathepsin B and L activities distinguish papillary carcinoma of the thyroid from normal thyroid or thyroid with non-neoplastic disease. Int J Cancer 66:420–426[CrossRef][Medline]
  25. Kusunoki T, Nishida S, Nakano T, Funasaka K, Kimoto S, Murata K, Tomura T 1995 Study on cathepsin B activity in human thyroid tumors. Auris Nasus Larynx 22:43–48[Medline]
  26. Schwartz MA, Ginsberg MH 2002 Networks and crosstalk: integrin signalling spreads. Nat Cell Biol 4:E65–E68
  27. Tse W, Zhu W, Chen HS, Cohen A 1995 A novel gene, AF1q, fused to MLL in t(1;11) (q21;q23), is specifically expressed in leukemic and immature hematopoietic cells. Blood 85:650–656[Abstract/Free Full Text]
  28. So CW, Ma SK, Wan TS, Chan GC, Ha SY, Chan LC 2000 Analysis of MLL-derived transcripts in infant acute monocytic leukemia with a complex translocation (1;11:4)(q21;q23;p16). Cancer Genet Cytogenet 117:24–27
  29. Tse W MS, Alonzo TA, Stirewalt DL, Gerbing RB, Woods WG, Appelbaum FR, Radich JP 2004 Elevated expression of the AF1q gene, an MLL fusion partner, is an independent adverse prognostic factor in pediatric acute myeloid leukemia. Blood 104:3058–3063[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Molecular Cancer TherapeuticsHome page
N. N. Co, W. P. Tsang, T. W.L. Wong, H. H. Cheung, T. Y. Tsang, S. K. Kong, and T. T. Kwok
Oncogene AF1q enhances doxorubicin-induced apoptosis through BAD-mediated mitochondrial apoptotic pathway
Mol. Cancer Ther., October 1, 2008; 7(10): 3160 - 3168.
[Abstract] [Full Text] [PDF]


Home page
J. Clin. Endocrinol. Metab.Home page
K. Krause, M. Eszlinger, O. Gimm, S. Karger, C. Engelhardt, H. Dralle, and D. Fuhrer
TFF3-Based Candidate Gene Discrimination of Benign and Malignant Thyroid Tumors in a Region with Borderline Iodine Deficiency
J. Clin. Endocrinol. Metab., April 1, 2008; 93(4): 1390 - 1393.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
T. Bawa-Khalfe, J. Cheng, Z. Wang, and E. T. H. Yeh
Induction of the SUMO-specific Protease 1 Transcription by the Androgen Receptor in Prostate Cancer Cells
J. Biol. Chem., December 28, 2007; 282(52): 37341 - 37349.
[Abstract] [Full Text] [PDF]


Home page
Endocr Relat CancerHome page
K. Fujarewicz, M. Jarzab, M. Eszlinger, K. Krohn, R. Paschke, M. Oczko-Wojciechowska, M. Wiench, A. Kukulska, B. Jarzab, and A. Swierniak
A multi-gene approach to differentiate papillary thyroid carcinoma from benign lesions: gene selection using support vector machines with bootstrapping
Endocr. Relat. Cancer, September 1, 2007; 14(3): 809 - 826.
[Abstract] [Full Text] [PDF]


Home page
Endocr. Rev.Home page
M. Eszlinger, K. Krohn, A. Kukulska, B. Jarzab, and R. Paschke
Perspectives and Limitations of Microarray-Based Gene Expression Profiling of Thyroid Tumors
Endocr. Rev., May 1, 2007; 28(3): 322 - 338.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
K. Kozma, J. J. Keusch, B. Hegemann, K. B. Luther, D. Klein, D. Hess, R. S. Haltiwanger, and J. Hofsteenge
Identification and Characterization of abeta1,3-Glucosyltransferase That Synthesizes the Glc-beta1,3-Fuc Disaccharide on Thrombospondin Type 1 Repeats
J. Biol. Chem., December 1, 2006; 281(48): 36742 - 36751.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Jacques, C.
Right arrow Articles by Reynier, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jacques, C.
Right arrow Articles by Reynier, P.
Related Collections
Right arrow Thyroid
Right arrow Endocrine Oncology


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals