| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology and Genetics |
1 Department of Genetics, Institute for Cancer Research and Departments of 2 Pathology and 3 Medical Oncology and Radiotherapy, The Norwegian Radium Hospital and University of Oslo; 4 Department of Chemical Toxicology, Division of Environmental Medicine, Norwegian Institute of Public Health, Oslo, Norway; 5 Biomedicum Biochip Center, Biomedicum, University of Helsinki, Helsinki, Finland; 6 Medical Biotechnology Group, VTT Technical Research Centre of Finland and University of Turku, Turku, Finland; and 7 Department of Biomedical Science, University of Sheffield, Sheffield, United Kingdom
Requests for reprints: Ragnhild A. Lothe, Department of Genetics, Institute for Cancer Research, The Norwegian Radium Hospital, N-0310 Oslo, Norway. Phone: 47-22934415; Fax: 47-22934440; E-mail: rlothe{at}radium.uio.no.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
|
To explore the transcriptional programs during germ cell tumorigenesis, we studied samples of each histologic subtype of TGCT, testicular parenchyma from healthy individuals, and samples taken from testis with IGCN. We also tested the RA-induced in vitro differentiation of the two embryonal carcinoma cell lines NTERA2 and 2102Ep. Particular attention was paid on dysregulation of gene categories known to be related to the early embryogenesis, such as the homeobox gene family, genes specifically expressed in embryonic stem cells, genes involved in the WNT, transforming growth factor-ß (TGF-ß), and Notch signaling pathways, and genes related to DNA methylation and imprinting.
| Materials and Methods |
|---|
|
|
|---|
|
Gene expression microarrays. We used Agilent Human 1A oligomicroarrays (GEO accession no. GPL885, Agilent Technologies, Palo Alto, CA) containing 60-mer DNA probes synthesized in situ in a 22k format. Of 19,061 spots, 18,086 are noncontrols and there are 17,086 unique transcript sequences, matching to 15,989 unique human genes.8
Probe preparations, hybridizations, image generation, and image analyses were done according to the manufacturer's protocol. Briefly, we obtained the total RNA fraction from ground tissue samples (in liquid N2) and cell culture pellets by using the Trizol reagent (Life Technologies, Rockville, MD). The RNA quality was evaluated by use of the Agilent 2100 Bioanalyzer (Agilent Technologies). Labeled cDNA was synthesized from 20 µg RNA (Fluorescent Direct Label kit, Agilent Technologies) in the presence of cyanine 3-dCTP for the test sample and cyanine 5-dCTP (Perkin-Elmer Life Sciences, Boston, MA) for the common reference, consisting of a pool of 10 human cell lines, including one from embryonal carcinoma (Universal Human Reference RNA, Stratagene, La Jolla, CA). Differentially labeled test and reference samples were mixed with Agilent control targets before hybridization onto the oligomicroarrays for 17 hours at 60°C in a rotating oven.
Data processing and statistics. The fluorescence intensities at the targets were detected by a laser confocal scanner (Agilent Technologies), and resulting images were processed using the Feature Extraction software, version 6.1.1.1 (Agilent Technologies). This included defining the spots, measuring intensities, flagging spots with inadequate measurements, subtracting local background, and LOWESS dye normalization. For spots that were not flagged as having inadequate measurements, ratios (sample over reference) of the processed intensities were used further. The measured ratios of all genes and samples were divided by the median of the ratios of the three normal samples before log2 transformation. This was done to facilitate the interpretation of the expression values, because now the sign of the value indicates relative up-regulation or down-regulation in the sample compared with the expression level in normal testis.
For the 60-mer DNA sequences that are represented more than once on the array (100 sequences are present in 10 times abundance), the medians were used for further analyses. The total gene expression matrix contained after this adjustment results for 17,086 unique sequences, of which 259 had technically inadequate measurements in >5 of the 29 samples and were excluded from the analyses. Thus, we ended up with 16,817 unique and valid transcript sequences represented by 15,756 individual gene symbols. For further analyses, missing values were estimated by the k-nearest neighbor imputation (k = 10).
Principal component analysis was done by input of the 1,575 genes that were altered in at least three samples to at least a 4-fold deviation from the median across the sample set (J-Express Pro version 2.0, MolMine, Bergen, Norway). Hierarchical clustering for improved visualization in Figs. 2 and 3 was done by average linkage of Euclidean distance measures (J-Express Pro).
|
|
In the first step of finding discriminator genes for the undifferentiated and pluripotent phenotype, a SAM analysis of embryonal carcinomas versus differentiated tissues was done (1.84% FDR; d = 0.965), resulting in 217 significant genes. These were again filtered, leaving only genes that were significant in the same direction from a second SAM analysis (10% FDR) comparing undifferentiated cell lines (the three samples from 2102Ep and the untreated sample from NTERA2) in one group and cell lines where differentiation were induced (the two treated NTERA2 samples) in the other.
Functional classes of gene products were assigned according to the Gene Ontology (GO) Consortium database.9 The total list of gene symbols with valid DNA microarray measurements and our lists of significant genes were along with the June 2004 version of the GO database uploaded into the GoMiner software (9). Here, the null hypothesis is that genes are evenly distributed among the GO categories represented by the genes on the microarray, and Fisher's exact tests are done to test the significance level for enrichment of genes within each GO category. Also for the pairwise comparisons of the self-curated gene lists (Table 2), Fisher's exact tests were used to estimate the statistical significance levels.
|
For GAL, cases with granular staining in the cytoplasm and/or the extracellular spaces were scored as positive. For POU5F1, tissue cores with staining in >5% of the relevant nuclei were scored as positive.
| Results |
|---|
|
|
|---|
To illustrate the relationship between expression profiles of the individual samples, principal component analysis was carried out (Fig. 1C). Here, samples from the same histologic subtypes appear in the vicinity of each other, demonstrating, as expected, that part of the variability in gene expression profiles of TGCTs is dependent on histologic type. A measure of the actual differences between the individual histologic subtypes is shown in the matrix of numbers of genes with expression medians deviating >3-fold between each pair of all combinations of the histologic subtypes (Fig. 1D). For example, there are 292 genes with a >3-fold higher expression value in seminomas than in embryonal carcinomas and vice versa and 182 genes that were more highly expressed by the embryonal carcinomas compared with seminomas. Interestingly, the pluripotent, undifferentiated cell lines have a similar number of genes differentially expressed from seminomas as they have from the embryonal carcinomas, the subtype that they are generally believed to model.
Histologic subtype-specific gene expression. We identified sets of differentially expressed genes within each of the histologic subgroups compared with the rest of the samples (individual SAM analyses; ref. 8). The top 10 up-regulated genes within each subgroup are shown in Fig. 2. From each of the seven histologic subtypes, 80 genes with the most significant up-regulation and the 20 most down-regulated genes were listed as interhistologic subtype-regulated genes (Supplementary Table S1). As 13 genes appeared for more than one histologic subtype, the complete set contained altogether 687 genes and not the expected 700. The GOs of the differentially expressed genes were explored, and several GO categories were significantly enriched within the interhistologic subtype-regulated gene list (Table 3; Supplementary Table S2). Significance levels of enrichment of genes belonging to embryogenesis-related gene categories were also calculated (Table 2). The 687 interhistologic subtype-regulated genes were particularly enriched among embryonic stem cellspecific genes (12 of 83, P = 2 x 104) identified by a DNA microarray study (12).
|
Undifferentiated and pluripotency-specific gene expression. To identify genes that are robustly differentially expressed between undifferentiated and pluripotent embryonal carcinoma and differentiated tissues or cell lines, a stepwise gene selection was implemented. First, SAM analysis was done, comparing the two extreme in vivo groups in terms of differentiation [i.e., the highly undifferentiated embryonal carcinomas (n = 5) versus the differentiated tissue types (n = 12; normal testis, yolk sac tumor, choriocarcinoma, and teratoma)]. This resulted in 176 genes that were most highly expressed in the embryonal carcinomas and 41 genes that were most highly expressed in the differentiated tissues (Table 4; Supplementary Table S5). The same analysis was also done, including only malignant tissue samples, and again, the same five genes, GAL, POU5F1, NANOG, DPPA4, and MT1H, appeared as the most significantly overexpressed in embryonal carcinoma (data not shown). The 217 (176 + 41) genes that were overall differentially expressed between the tissue samples were filtered further, leaving only those genes showing up-regulation or down-regulation in the in vitro models of TGCT differentiation (altered in the same direction after RA treatment for the pluripotent NTERA2 but remained unchanged in the relatively nullipotent 2102Ep). This left us with 68 genes, 58 that are characteristic of the undifferentiated embryonic stem celllike state and 10 that are characteristic of the differentiated derivatives (Fig. 3). Significantly enriched GO terms assigned to these 68 genes are listed in Table 5. The involvement of these genes in various gene function categories related to developmental biology are found from Table 2. Gene categories, such as DNA methylation genes, homeobox genes, and genes highly expressed in embryonic stem cells, were particularly abundant among the 68 genes specifically expressed in undifferentiated embryonal carcinomas.
|
|
Protein expression of GAL and POU5F1. The in situ protein expression of GAL and POU5F1, two genes identified as significantly overexpressed in embryonal carcinoma, were analyzed on tissue microarrays containing 510 tissue samples from normal testis, IGCN, and the various histologic subgroups of TGCT. Both GAL and POU5F1 were negative for all normal testis samples but were positive for 21% and 54% of the overall tumor samples, respectively, and both with significant differences between various histologic subtypes (Fig. 4). Interestingly, the frequency of POU5F1-positive seminomas (63%) was similar to that of embryonal carcinomas (65%), whereas for GAL this was different (8% positive seminomas and 34% positive embryonal carcinomas; P = 5 x 106).
|
| Discussion |
|---|
|
|
|---|
Histologic subtype-specific gene expression and parallels to early embryogenesis. The transcriptional profiling of individual histologic subtypes using SAM analyses revealed mostly novel genes with respect to testicular tumorigenesis. However, some previously known biological and clinical markers also appeared and add support to the importance of the resulting genes, reliability of the technology, and strength of the statistics [e.g.,
-fetoprotein (AFP) was highly expressed in yolk sac tumors, and subunits of the chorionic gonadotropin (CGB5 and CGA) were highly expressed in choriocarcinomas (Fig. 2)]. In addition, when comparing normal testis and IGCN (Supplementary Fig. S2), several genes that are expressed in the IGCN samples, such as KIT and the pluripotency markers NANOG and POU5F1, were also identified in a recent DNA microarray study focusing on IGCN (23). Among the novel genes to this study, DPPA4 was specifically expressed in IGCN compared with normal testis, and interestingly, the expression of these pluripotency-specific genes was also seen in seminomas (Fig. 3).
When we searched for gene expression patterns characterizing the pluripotent and undifferentiated phenotype of embryonal carcinomas, we did therefore not include IGCN and seminomas in any of the groups compared to avoid exclusion of the pluripotency-specific markers. The 68 resulting genes with specific expression pattern in the undifferentiated and pluripotent phenotype of embryonal carcinomas included 7 genes that were reported previously as enriched in embryonic stem cells (P = 1 x 107; Table 2; ref. 12), and 2 of these, NANOG and POU5F1, were also abundantly present in IGCN and seminomas. Both these genes are known as key regulators of pluripotency (2426), and IGCN and seminoma cells might therefore have similar pluripotent capabilities as embryonal carcinoma cells. Also on the protein level, as seen from immunohistochemistry on the tissue microarrays, POU5F1 was present in most IGCN and seminomas. This was not the case for GAL, which was highly associated to embryonal carcinomas compared with seminomas. However, a few seminomas were also positive for GAL, which lends support to the existence of a distinct subgroup of seminomas that resembles embryonal carcinomas on the molecular level.11
For the nonseminomas, morphologic parallels to early embryogenesis are described. In the same way that embryonal carcinomas share many characteristics with the inner cell mass of blastocysts (27), the choriocarcinomas parallel the trophoblastic and syncytiotrophoblastic cells of the placenta, yolk sac tumors parallel the endodermal differentiation giving rise to the yolk sac, and teratomas parallel the differentiation into somatic tissues of all three primary germ layers (2). Hence, transcriptional programs of the early embryogenesis are likely to be altered in testicular tumorigenesis. The enrichment of embryonic stem cellspecific genes in testicular tumorigenesis validates this morphologic link to early embryogenesis also on the gene transcription level. Hence, by presenting gene expression patterns specific to the histologic subtypes of human TGCT (Fig. 2), we also provide novel transcriptional information relevant to the parallel differentiation stages of human development.
Pluripotent embryonal carcinomas are considered to be the malignant counterparts of embryonic stem cells (28). Embryonic stem cells are again potential mediators to regenerative medicine, but as transplants may be rejected from the host's immune system, it would be advantageous if personalized embryonic stem cells could be made from human adult stem cells by dedifferentiation back to the embryonic stage without the need for cloning (29). The identified embryonal carcinomaspecific transcripts (as in Fig. 3) can potentially work as reprogramming factors that lead to dedifferentiation of adult cells into personalized embryonic stem cells. Human embryonic stem cells have thus far been considered pluripotent, rather than totipotent, because of their limited differentiation potency into extraembryonic membranes and tissue types (30). Recently, however, they were shown to have the capability to differentiate into trophoblast cells (31), an observation that may broaden the definition of their developmental potential. However, as TGCTs often contain tissue types with extraembryonic morphology (yolk sac tumor and choriocarcinoma), the current histologic subtype-specific gene expression profiles (Fig. 2) may indeed contain information that can be useful in the search for ways to make the embryonic stem cells differentiate into these lineages; e.g., by overexpression of genes specifically expressed in the TGCT subtype paralleling the desired lineage differentiation.
Gene categories involved in early embryogenesis. Knowledge from the field of embryogenesis may also be used to understand germ cell tumorigenesis, and we have identified subsets of genes that are also known to play a role in early embryonic development (Table 2). For example, all the four most significantly overexpressed genes in the undifferentiated embryonal carcinomas, POU5F1, GAL, NANOG, and DPPA4 (Table 4), are known as pluripotency-associated genes (12). For WNT signaling, we found seven components of this pathway within the interhistologic subtype-regulated genes list (Table 2). Except FRAT2, which was highest expressed in seminomas and embryonal carcinomas, these WNT signaling components (AES, DKK1, FZD4, TCF7L1, and TCF7L2) were highest expressed in yolk sac tumors and teratomas. Further, the three significantly altered target genes of this signaling pathway, FZD7, ID2, and, SLC2A1, were also up-regulated in the yolk sac tumors and teratomas. Thus, activation of the WNT signaling pathway may be a key step in the transition of embryonal carcinomas into differentiated nonseminomas.
High expression of the de novo DNA methyltransferase DNMT3B and its homologue involved in establishment of imprinting, DNMT3L, was significantly associated with the embryonal carcinoma subtype (Fig. 3; Tables 2 and 4). The third homologue, DNMT3A, had a similar expression profile not reaching statistical significance. This discovery of an up-regulated DNA methylation machinery in embryonal carcinomas provide a malignant parallel to the complete epigenetic reprogramming events taking place in the inner cell mass of the blastocysts (32, 33).
In vivo versus in vitro clues to potential drug targets. The RA-induced differentiation of embryonal carcinoma cell lines in this study is used to model the in vivo differentiation of embryonal carcinomas. The shared transcriptional changes between the in vivo and the in vitro systems (Fig. 3) are of special interest as they may be specific and/or necessary to make embryonal carcinoma differentiate and play biologically as well as clinically important roles. In fact, one of the genes sorted out by this approach, NALP7, was recently confirmed to play a functional oncogenic role in TGCT (34). In the treatment of TGCT, it could therefore be advantageous to repress the genes specifically expressed in the undifferentiated cells (e.g., GAL, POU5F1, NANOG, and DPPA4), as growth and malignant potential could be repressed by forcing the cells into terminal differentiation. Further intriguing, cancer drugs targeting genes that in the normal setting are specifically expressed during the blastocyst stage should in the postembryonal body be specific to tumor cells.
The current data could also facilitate cancer drug discovery by use of the genes always expressed in the differentiated tissues but not in the undifferentiated ones (e.g., MGC5576, TTC3, or C9orf154; Table 4), genes that may be used as general markers of differentiation. Their gene products can be used in drug screening programs where the agents causing expression of these selective markers in embryonal carcinoma cells are considered candidate cancer drugs.
In conclusion, we have presented the first comprehensive gene expression profiling of all histologic subtypes of TGCT as well as of precursor lesions and normal testis. This pinpointed specifically expressed genes within all analyzed subgroups. We also uncovered novel transcriptional information about tumorigenesis and its functional relation to normal development. The key genes described, as GAL and POU5F1, will be useful in molecularly assisted diagnosis as well as serve as potential drug targets.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Christine Pigott for the cell culturing work and Tuula Airaksinen for the demonstration of DNA microarray protocols.
| Footnotes |
|---|
8 TIGR Resourcerer 10.0 July 2004 Release (http://www.tigr.org/tigr-scripts/magic/r1.pl). ![]()
9 Gene Ontology Consortium downloadable database as of June 2004 (http://www.godatabase.org/dev/database/). ![]()
10 http://www.ncbi.nlm.nih.gov/geo/ ![]()
11 Hofer et al., unpublished data. ![]()
Received 1/17/05. Revised 3/15/05. Accepted 4/13/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Godmann, I. Gashaw, K. Eildermann, S. Schweyer, M. Bergmann, R.I. Skotheim, and R. Behr The pluripotency transcription factor Kruppel-like factor 4 is strongly expressed in intratubular germ cell neoplasia unclassified and seminoma Mol. Hum. Reprod., August 1, 2009; 15(8): 479 - 488. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. B. Sonne, K. Almstrup, M. Dalgaard, A. S. Juncker, D. Edsgard, L. Ruban, N. J. Harrison, C. Schwager, A. Abdollahi, P. E. Huber, et al. Analysis of Gene Expression Profiles of Microdissected Cell Populations Indicates that Testicular Carcinoma In situ Is an Arrested Gonocyte Cancer Res., June 15, 2009; 69(12): 5241 - 5250. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lalancette, C. Thibault, I. Bachand, N. Caron, and N. Bissonnette Transcriptome Analysis of Bull Semen with Extreme Nonreturn Rate: Use of Suppression-Subtractive Hybridization to Identify Functional Markers for Fertility Biol Reprod, April 1, 2008; 78(4): 618 - 635. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Y. Kim, M. K. Kee, S. A. Chong, and M. J. Nam Galanin Is Up-Regulated in Colon Adenocarcinoma Cancer Epidemiol. Biomarkers Prev., November 1, 2007; 16(11): 2373 - 2378. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Rajpert-De Meyts Developmental model for the pathogenesis of testicular carcinoma in situ: genetic and environmental aspects Hum. Reprod. Update, May 1, 2006; 12(3): 303 - 323. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. Korkola, J. Houldsworth, R. S.V. Chadalavada, A. B. Olshen, D. Dobrzynski, V. E. Reuter, G. J. Bosl, and R.S.K. Chaganti Down-Regulation of Stem Cell Genes, Including Those in a 200-kb Gene Cluster at 12p13.31, Is Associated with In vivo Differentiation of Human Male Germ Cell Tumors Cancer Res., January 15, 2006; 66(2): 820 - 827. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |