Next Article in Journal
Pomalidomide Ameliorates H2O2-Induced Oxidative Stress Injury and Cell Death in Rat Primary Cortical Neuronal Cultures by Inducing Anti-Oxidative and Anti-Apoptosis Effects
Next Article in Special Issue
Roles of Tristetraprolin in Tumorigenesis
Previous Article in Journal
Translational Research Opportunities Regarding Homologous Recombination in Ovarian Cancer
Previous Article in Special Issue
A G-Quadruplex Structure in the Promoter Region of CLIC4 Functions as a Regulatory Element for Gene Expression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

PR/SET Domain Family and Cancer: Novel Insights from The Cancer Genome Atlas

by
Anna Sorrentino
1,2,†,
Antonio Federico
2,3,†,
Monica Rienzo
4,
Patrizia Gazzerro
5,
Maurizio Bifulco
6,
Alfredo Ciccodicola
2,3,
Amelia Casamassimi
1,* and
Ciro Abbondanza
1,*
1
Department of Precision Medicine, University of Campania “Luigi Vanvitelli”, Via L. De Crecchio, 80138 Naples, Italy
2
Department of Science and Technology, University of Naples “Parthenope”, 80143 Naples, Italy
3
Institute of Genetics and Biophysics “Adriano Buzzati Traverso”, CNR, 80131 Naples, Italy
4
Department of Environmental, Biological, and Pharmaceutical Sciences and Technologies, University of Campania “Luigi Vanvitelli”, 81100 Caserta, Italy
5
Department of Pharmacy, University of Salerno, 84084 Salerno, Italy
6
Department of Molecular Medicine and Medical Biotechnologies, University of Naples “Federico II”, 80131 Naples, Italy
*
Authors to whom correspondence should be addressed.
These two authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(10), 3250; https://doi.org/10.3390/ijms19103250
Submission received: 30 September 2018 / Revised: 12 October 2018 / Accepted: 17 October 2018 / Published: 19 October 2018

Abstract

:
The PR/SET domain gene family (PRDM) encodes 19 different transcription factors that share a subtype of the SET domain [Su(var)3-9, enhancer-of-zeste and trithorax] known as the PRDF1-RIZ (PR) homology domain. This domain, with its potential methyltransferase activity, is followed by a variable number of zinc-finger motifs, which likely mediate protein–protein, protein–RNA, or protein–DNA interactions. Intriguingly, almost all PRDM family members express different isoforms, which likely play opposite roles in oncogenesis. Remarkably, several studies have described alterations in most of the family members in malignancies. Here, to obtain a pan-cancer overview of the genomic and transcriptomic alterations of PRDM genes, we reanalyzed the Exome- and RNA-Seq public datasets available at The Cancer Genome Atlas portal. Overall, PRDM2, PRDM3/MECOM, PRDM9, PRDM16 and ZFPM2/FOG2 were the most mutated genes with pan-cancer frequencies of protein-affecting mutations higher than 1%. Moreover, we observed heterogeneity in the mutation frequencies of these genes across tumors, with cancer types also reaching a value of about 20% of mutated samples for a specific PRDM gene. Of note, ZFPM1/FOG1 mutations occurred in 50% of adrenocortical carcinoma patients and were localized in a hotspot region. These findings, together with OncodriveCLUST results, suggest it could be putatively considered a cancer driver gene in this malignancy. Finally, transcriptome analysis from RNA-Seq data of paired samples revealed that transcription of PRDMs was significantly altered in several tumors. Specifically, PRDM12 and PRDM13 were largely overexpressed in many cancers whereas PRDM16 and ZFPM2/FOG2 were often downregulated. Some of these findings were also confirmed by real-time-PCR on primary tumors.

Graphical Abstract

1. Introduction

The positive regulatory domain (PRDM) gene family, a subfamily of Kruppel-like zinc finger gene products, currently includes 19 members in humans [1,2,3,4]. The protein products of this family share a conserved N-terminal PR (PRDI-BF1-RIZ1 homologous) domain, which is structurally and functionally similar to the catalytic SET domain that defines a large group of histone methyltransferases (HMTs) [5,6,7,8]. So far, enzymatic activity has been experimentally demonstrated for only a few family members; otherwise, PRDM proteins (PRDMs) lacking intrinsic enzymatic activity are able to recruit histone-modifying enzymes to mediate their function.
The PR domain is generally positioned at the protein N-terminal region and, with the exception of PRDM11, it is followed by repeated zinc fingers toward the C-terminus, potentially mediating sequence-specific DNA or RNA binding and protein–protein interactions [5,6,7,8,9] (Figure S1). Importantly, PRDMs have been established to tether transcription factors to target gene promoters by recognition of a specific DNA consensus sequence [10] or acting as non-DNA binding cofactors [11,12]. These features give PRDMs the ability to drive cell differentiation and to specify cell fate choice and, thus, contribute to many developmental processes [5,7,12,13,14].
A common characteristic of most PRDM genes is to express two main molecular variants, one lacking the PR domain (PR-minus isoform) but otherwise identical to the other PR-containing product (PR-plus isoform). These two isoforms, generated either by alternative splicing or alternative use of different promoters [7,8,15], play opposite roles, particularly in cancer. The full-length product PR-plus usually acts as a tumor suppressor, whereas the short isoform functions as an oncogene. This bivalent behavior has been tagged as ‘yin-yang’. The imbalance in favor of the PR-minus is found in many human malignancies and it can be due to inactivating mutations or silencing of the complete form and/or to increased expression of the PR-minus form [8].
PRDM1 and PRDM2 use alternative promoters to generate short isoforms lacking the PR domain, which show oncogenic properties. Increased levels of the short isoforms were reported in various cancer cell lines. A similar PR-less product was also described for PRDM3/MECOM (MDS1 and EVI1 complex locus), PRDM16, and PRDM6, thus suggesting that this ‘yin-yang’ expression pattern and its functional implications could be a hallmark of most, if not all, PRDMs [8].
Furthermore, several studies have described alterations (both mutations and/or gene expression changes) of most PRDMs in solid tumors and/or hematological malignancies [8]. For instance, frameshift mutations of microsatellite repeats within the PRDM2 coding region are frequent events in various cancers. A recent study has described a frameshift mutation in the C-terminal region of PRDM2, affecting the (A)9 repeat within exon 8, as a microsatellite indel driver hotspot and as a driver mutation in microsatellite instability (MSI) colorectal cancer [15,16]. Notably, a similar frameshift mutation was found to occur in a mononucleotide repeat (A7) of PRDM3/MECOM gene in this cancer type [17]. Intriguingly, recent findings also indicate that PRDM2 methyltransferase is required for BRCA1-mediated genome maintenance [15,18]. Moreover, a significant reduction of PRDM2 expression was observed in high-grade gliomas [19], and forced expression of PRDM2 in glioma cell lines inhibits cell proliferation and increases apoptosis. This evidence strongly suggests a possible tumor suppressive role for PRDM2 [19]. Interestingly, PRDM9 HMTase activity is essential for meiotic DNA double-strand break formation at its binding sites [20,21]. Moreover, both PRDM1 and PRDM5 negatively modulate WNT/β-catenin signaling, a pathway involved in the occurrence of several cancers, including glioma and colorectal cancer [22,23].
This evidence suggests that PRDMs are involved in human cancer through modulation of several processes, such as epigenetic modifications, genetic reprogramming, inflammation, and metabolic homeostasis.
To date, both mutations and altered expression have been reported for some PRDMs in specific cancer entities. However, our understanding of the role played by different PRDM family members in cancer is still limited by the lack of a systematic and comprehensive approach in deciphering the mutational status and the complete transcriptional profile of all the PRDMs across a large number of different cancer types.
Here, The Cancer Genome Atlas (TCGA) deposited exome and RNA-Seq data [24] were used to obtain a complete pan-cancer overview of the genomic and transcriptomic alterations for all PRDM genes across 31 distinct human cancer types.

2. Results

2.1. Mutational Profiling of PRDM Genes Across Human Cancers

To systematically identify somatic mutations within genes encoding PRDMs, we started with a mutational profiling of these genes across human cancers. We downloaded Exome-sequencing datasets from the TCGA web portal for 31 cancer types and about 11,000 patients. The number of samples for each cancer type is illustrated in Table S1 [25].
Overall, we identified 3995 point mutations, 180 deletions (39 in-frame and 141 frameshift), and 22 insertions (16 in-frame and 6 frameshift) affecting PRDM genes. Silent or synonymous mutations were 1531 (26.7% of total mutations) and ranged between 11% (PRDM6) and 41% (PRDM8) of the total mutations for each gene (Figure 1).
According to our reanalysis, the most mutated genes were PRDM2 (507 mutations; with 24% of silent mutations), PRDM3/MECOM (547 mutations; 22% silent), PRDM9 (899 mutations; 27% silent), PRDM16 (514 mutations; 31% silent), and ZFPM2/FOG2 (700 mutations; 24% silent). Non-sense mutations were more recurrent in PRDM5 (23), PRDM9 (60), and ZFPM2/FOG2 (34), whereas splice sites disrupting mutations were more frequently detected in PRDM3/MECOM (13), PRDM9 (19), PRDM10 (14) and PRDM12 (13) (Figure 1).
To measure the frequencies of somatic mutations for each PRDM gene across all tumor types, only non-synonymous mutations were considered. We observed heterogeneity in the mutation frequencies of these genes in the different tumor types. A global low mutation rate (from 0 to 8.2%) was found, except for PRDM3/MECOM, PRDM8, PRDM9, PRDM15, ZFPM1/FOG1, and ZFPM2/FOG2 (Table 1). In detail, PRDM8 and PRDM15 were mutated at low rates in most of the analyzed cancer types except PAAD where they were both frequently mutated (16.0% and 11.2%, respectively). PRDM3/MECOM was recurrently mutated in various cancer types, also reaching a value of 20.1% of mutated samples in SKCM. Similarly, PRDM9 was mutated with a high mutation rate in many cancer types, achieving values of 10.0% in UCEC, 14.2% in LUAD, and 15.4% in SKCM. Otherwise, ZFPM1/FOG1 was mutated at a low rate in a few cancer types, except in UCS (5.2%), COAD (6.6%), READ (9.4%), and ACC (50.5%). Finally, ZFPM2/FOG2 was frequently mutated at a high rate in various cancer types, reaching a value of 11.1% in LUAD and 16.5% in SKCM.
We visualized the mutation data in each tumor type by Oncostrip function (Supplementary file 1). Through this approach, we evaluated the percentage of samples with at least one mutated PRDM gene in each tumor type ranging from 1.02% (2/196) in LAML samples to 55.43% (51/92) in ACC samples. Furthermore, this function allowed us to visualize the mutation type affecting PRDM genes in each sample. Interestingly, ZFPM1/FOG1 revealed a high number of samples, especially in ACC, with “multi_hit” mutations (more than one mutation affecting the same gene in the same cancer sample). Specifically, we found 11/18 (61%) multi-hit mutations in COAD, 10/11 (90%) in READ, and 23/47 (48%) in ACC (see Supplementary file S1).
To distinguish between damaging and tolerated missense mutations, we carried out a variant effect predictor (VEP) analysis (Table S2). Missense mutations with a SIFT score ranging in the interval 0.0–0.05 and/or with a PolyPhen score in the interval 0.5–1 were considered as deleterious or probably damaging, respectively. As shown in Table 2, adding all the other deleterious somatic mutations (frameshift, in-frame deletions, stop gained and start lost mutations, splice site, UTR, and intron variants) to the deleterious missense mutations classified with the VEP analysis, we obtained the total number of deleterious mutations affecting each PRDM gene. Thus, we obtained the percentage of deleterious somatic mutations across the tumor samples. This number was ≥50% for MECOM/PRDM3 (52.7%), PRDM4 (55%), PRDM5 (54.8%), PRDM6 (58.5%), PRDM10 (55.2%), PRDM11 (54%), PRDM13 (51.7%), and PRDM16 (50%).
Additionally, to predict the potential functional effect of the identified PRDM somatic mutations on the affected proteins and to detect a possible mutation enrichment in some domains, we localized the deleterious missense mutations on the canonical protein isoform of each PRDM (Figure 2). Interestingly, a random sampling weighted on the size of the annotated protein domains demonstrated that somatic deleterious mutations were significantly enriched in the PR domain of PRDM1, PRDM5, PRDM6, PRDM8, PRDM9, PRDM12 and PRDM13 (p < 0.005).
Another important aspect of cancer genetic studies is the presence of possible recurrent and hotspot mutations. Figure 3 illustrates mutations in PRDM genes recurring in more than three tumor types. Interestingly, the frameshift mutation T/-→K678X, despite affecting PRDM3/MECOM in a region not containing known domains, was recurrent in different tumor types; similarly, also the missense mutation G/A→S237L occurred in a region without known domains but in many tumors. Otherwise, missense mutations affecting a Zn-finger domain and occurring in different tumors were observed for PRDM9, PRDM14, and PRDM16. Likewise, PRDM12 was frequently mutated in a splice donor site in a region coding for the PR domain whereas in different tumor types, ZFPM2/FOG2 was affected by the missense mutation C/T→R734C in a region without known domains. In addition, PRDM2 and PRDM15 revealed an in-frame deletion in various cancers and PRDM11 a frameshift mutation. Finally, ZFPM1/FOG1 showed several recurrent mutations; they all (frameshift mutations and in-frame deletions) hit a region without known domains (Figure 3).
Interestingly, all these mutations were particularly recurrent in ACC patients. In this cohort, ZFPM1/FOG1 also displayed five hotspot mutations, all localized in the same region outside the known domains (Figure 4a). To establish whether these hotspot mutations could have an impact on the ZFPM1/FOG1 structure, we utilized the I-TASSER web-tool to predict the tertiary structure of the annotated ZFPM1/FOG1 protein (Figure 4b) and proteins carrying the missense mutations and the in-frame deletions (Figure 4c–e). As illustrated, these mutations completely altered the structure of the canonical protein. Otherwise, the frameshift mutations E444X and P445X led to premature stop codons at the residues 669 and 796, respectively; both of the mutated proteins shared only the first 443 residues with the canonical protein whereas they changed in the 444–669 and 444–796 regions and missed respectively 337 and 210 residues at the C-terminal, which contains the last five zinc fingers of ZFPM1/FOG1.
Finally, to assess whether members of the PRDM family may be driver genes in a given cancer type, we used the OncodriveCLUST tool, which aims to identify genes whose mutations are biased towards a large spatial clustering. This method is based on the feature that cancer gene mutations frequently cluster in particular positions of the protein. Thus, mutations with a frequency higher than the background rate that tend to cluster in specific regions of protein-coding genes are likely to be driver genes. Based on the scores of this analysis, ZFPM1/FOG1 can be considered as a cancer driver for ACC (Figure 4f) and PRDM8 for PAAD (Figure S2).

2.2. Differentially Expressed PRDM Genes across Human Cancers

To evaluate whether the expression of PRDM genes is affected in human cancers, we took advantage of RNA-Seq datasets from paired samples (cancer vs. benign counterpart) available at the TCGA web portal. Globally, 585 patients across 21 cancer types were analyzed (Table S1). The gene expression profiles differed considerably between normal and tumor specimens, depending on the cancer type, as shown by the principal component analysis [26]. The results of gene expression profiling are summarized in Table S3 and Figure 5.
Data indicate that a large subset of PRDM genes is consistently deregulated across several cancer types. Particularly, a significant overexpression was measured for PRDM12 and PRDM13. On the other side, ZFPM2/FOG2, PRDM8, and PRDM16 were more often downregulated across tumors (Figure 5). Strong upregulation of almost all PRDM genes was measured in CHOL where 13/19 PRDM genes were overexpressed in tumor versus healthy counterparts. Among them, the most upregulated were MECOM/PRDM3 (FC = 12.64), PRDM5 (FC = 10.7) and PRDM16 (FC = 13.74). Conversely, the cancer types with the smallest number of deregulated PRDM genes were SARC with no PRDM genes deregulated, PAAD with 1 gene upregulated (MECOM/PRDM3, FC = 2.71), THYM with two genes strongly downregulated (PRDM1, FC = 0.13 and ZFPM2/FOG2, FC = 0.16), CESC with two genes strongly downregulated (PRDM8, FC = 0.07 and ZFPM2/FOG2, FC = 0.04) and one gene strongly upregulated (PRDM13 FC = 29.86), and, finally, ESCA with two genes strongly downregulated (PRDM16, FC = 0.37 and ZFPM2/FOG2, FC = 0.44) and one gene upregulated (PRDM1, FC = 2.28).

2.3. PRDM Expression in Human Primary Tumors

In an attempt to validate the findings obtained on RNA-Seq datasets, we assayed TissueScan cDNA panel arrays containing eight different tumors (breast, colon, kidney, liver, lung, ovary, prostate, and thyroid) [25].
Specifically, we analyzed the differential expression of PRDM3/MECOM, PRDM10, PRDM12, PRDM16, and ZFPM2/FOG2 genes (Figure S3). As illustrated in Figure 6, a general differential expression, even though not significant, was observed for both ZFPM2/FOG2 and PRDM16 in several tumor tissues (Figure 6a,b). PRDM3/MECOM was found to be significantly overexpressed in breast (p < 0.001), ovary (p < 0.001), and colon (p < 0.05) cancer specimens (Figure 6c). Similarly, PRDM10 was overexpressed at significant levels in breast (p < 0.001) and colon (p < 0.05) cancer specimens (Figure 6d). No statistical differences in the expression of all these genes were measured between tumor and healthy samples of the other cancer entities (Figure 6a–d). Otherwise, the PRDM12 gene was difficult to analyze. This gene was confirmed as being highly overexpressed in all the analyzed cancer tissues, except in ovarian cancer, as indicated by reanalysis of the TCGA RNA-Seq dataset. However, this gene was not expressed or was expressed at very low levels in all the healthy tissues used as controls; it was expressed in tumor specimens. For this reason, its relative gene expression could not be calculated using the 2−ΔΔCt method for all the analyzed tissues. Particularly, relative expression was measured only for thyroid, ovary, prostate, and liver (Figure 6e) whereas in breast, colon, kidney, and lung normal tissues, the amplification products were undetectable when observed by agarose gel electrophoresis analysis (Figure 6f). Of note, these cancer specimens showed measurable levels of PRDM12 (Figure 6e–f).

3. Discussion

In this study, we provide for the first time a systematic and comprehensive overview of both the mutational status and the expression profile of all the PRDM genes across a large number of different cancer types.
Recently, the availability of multi-omics datasets (such as TCGA) from human cancers, together with the development of advanced bioinformatics tools, represent a unique source to study human malignancies [24].
Our reanalysis of the TCGA Exome-sequencing datasets revealed PRDM2, PRDM3/MECOM, PRDM9, PRDM16 and ZFPM2/FOG2 as the most mutated genes. Heterogeneity in the mutation frequencies was observed in the different tumor types with higher mutation rates found for PRDM3/MECOM, PRDM8, PRDM9, PRDM15, ZFPM1/FOG1 and ZFPM2/FOG2 in specific cancers. Remarkably, VEP analysis indicated that the percentage of total deleterious mutations across the tumor samples was high for most genes. More interestingly, a random sampling weighted on the size of the annotated protein domains demonstrated that deleterious mutations were significantly enriched in the PR domain of PRDM1, PRDM5, PRDM6, PRDM8, PRDM9, PRDM12 and PRDM13. Frequent mutations disrupting the PR domain in tumor samples would be a mechanism for removing the tumor suppressor function of the PR-plus isoform in favor of the oncogenic PR-minus form.
A big challenge in cancer biology studies is distinguishing between mutations conferring a selective growth advantage to cancer cells (drivers) and those randomly accumulating and without significant effects on the oncogenic process (passengers) [27,28]. Many algorithms employing different approaches are now utilized to recognize driver genes although, when compared for their performance, all display both strengths and weaknesses [25,29,30,31]. Moreover, recent studies have highlighted the existence of “mini-driver” genes with weaker tumor-promoting effects, thus expanding the previous driver–passenger dichotomy to a continuous model [25,30,32]. Besides, a sub-classification has been proposed to differentiate “mut-driver genes”, usually altered by somatic gene mutations, from “epi-driver genes”, which are deregulated through epigenetic modifications but are not frequently mutated [25,27].
In this study, OncodriveCLUST analysis revealed two putative cancer mut-driver genes: PRDM8 in PAAD and ZFPM1/FOG1 in ACC. In the latter case, we found that the involved gene was mutated in a very high percentage (about 55%) of tumor samples. Additionally, a mutational hotspot region localized between the amino acid positions 444 and 447, outside the known domains, was recognizable. All these findings agree with the key parameters commonly used to discern drivers from passengers [29]. Notably, these hotspot mutations recurred also in other malignancies, such as COAD, KIRP, READ, STAD, and UCS, supporting an important role of this gene in carcinogenesis. Moreover, the finding of “multi_hit” mutations in ACC, as well as in COAD and READ, advises that this gene could function as a tumor suppressor gene. This is conceivable with the current knowledge about the role of ZFPM1/FOG1 in differentiation. Indeed, ZFPM1 is also known as a friend of GATA1 (FOG1) since it interacts with GATA1 and it is an essential cofactor for the transcription factor GATA-1 in erythroid and megakaryocytic differentiation. Reduced expression of ZFPM1/FOG1 was found in preleukemic progenitors of a mouse model of leukemia [33,34]. Besides, the high recurrence of these mutations together with results from 3D-modeling of the canonical and mutant ZFPM1/FOG1 proteins suggest that, although without known domains, this is a critical region for the protein. Interestingly, this region (particularly K431 residue) is evolutionary highly conserved (Figure S4).
It is conceivable that other PRDM genes may play a key role in the initiation and progression of specific or multiple tumor types, as also reported by previous literature [8,15,16,17]. Indeed, it is accepted that cancer driver genes are mainly involved in three core cellular processes: cell fate, cell survival, and genome maintenance. For instance, PRDM2 has a relevant role in all of them; specifically, it also participates to the formation of protein complexes involved in the DNA damage response and in genome maintenance [15]. Additionally, a recent study identified a driver PRDM2 mutation in MSI colorectal cancer [16]. However, our analysis has not considered this gene as a driver. We cannot exclude, among the explanations, the limitation of the utilized bioinformatics tool [30].
Our pan-cancer study has also identified PRDM12 as a possible epi-driver gene in multiple cancers. Of note, TCGA gene expression profiling of PRDM genes revealed significant overexpression of PRDM12 and PRDM13 in many tumor types whereas ZFPM2/FOG2, PRDM8, and PRDM16 were more often downregulated in tumor tissues. Our qRT-PCR analysis was not able to confirm all the results obtained through TCGA analysis. The main reason could be the utilization of unpaired samples in the validation through TissueScan cDNA panel arrays; otherwise, our analysis on TCGA was carried out on paired samples. In addition, we have analyzed a small number of samples by qRT-PCR compared to the huge number of cases from TCGA. Noteworthy is that when we measured the differential expression of PRDM12 in cDNA panel arrays, we found the expression of this gene only in cancer specimens but not in healthy samples of several tissues, suggesting that it could be putatively utilized as a biomarker in those malignancies. Our study represents the first analysis of all PRDM genes in pan-cancer; further studies using large cohorts are necessary to validate the most promising results, particularly for PRDM12. In addition, given the lack of literature data, we are aware that functional studies investigating the effect of altered expression both in vitro and in vivo are required to establish the possible impact of PRDM12 in cancer.
Altogether, our results can be useful for identifying a subset of relevant PRDMs that are frequently mutated and/or transcriptionally deregulated in certain tumor types. Functional studies on specific PRDM gene mutations should be accomplished to definitely prove their potential oncogenic role. Moreover, it would be interesting to investigate whether these mutations contribute to cancer progression and metastasis, as well as whether they correlate with prognosis and/or with drug response and resistance. The epigenetic changes underlying the altered gene expression observed in tumor samples should also be explored. In this context, the availability of novel multi-omics data integration tools and methods also offer the opportunity to further integrate our analysis of PRDM gene expression by a systematic pan-cancer study of the epigenetic marks in these genes [25,35,36].

4. Materials and Methods

4.1. TCGA Data Source Selection and Processing for Mutation Analysis

In this manuscript, we analyzed both whole Exome- and RNA-Seq data retrieved from publicly accessible repositories. Specifically, we retrieved the whole exome sequencing data from the GitHub R data package for pre-compiled somatic mutations from TCGA cohorts “TCGA mutations” and analyzed it using the Bioconductor package “maftools” [37].
The selection and nomenclature of PRDM genes were based on the HUGO Gene Nomenclature Committee [38]. To estimate the mutation enrichment within the PR domain of each of the PRDM proteins, we performed a random sampling iterated 1000 times weighted on the size of the annotated domains.
To assess whether one or more PRDM proteins could be considered as cancer driver genes based on the positional clustering of the variants in the selected human cancers, we used a re-implementation of the software OncodriveCLUST within the maftools package [39].
Three-dimensional (3D) modeling of the human canonical and mutant ZFPM1/FOG1 proteins was carried out using I-TASSER [26,40,41].

4.2. TCGA Data Source Selection and Processing for Expression Analysis

The RNA-Seq gene expression data were downloaded from TCGA [42]. The analysis of gene expression and the identification of differentially expressed genes were performed comparing the expression profiles of cancer vs. matched normal samples in a paired analysis. Therefore, expression data taken from human primary cancers for which healthy samples were not available were discarded. According to this criterion, 22 tumor entities were analyzed. To have a more robust differential expression analysis in paired samples, we applied generalized linear models implemented in the EdgeR Bioconductor package version 3.17.10. p-values adjustment was performed through the application of the false discovery rate (FDR) method. We considered differentially expressed genes with a logFC ≤−1 and logFC ≥1, and an FDR ≤ 0.01.

4.3. Real-Time RT-PCR Analysis

Quantitative real-time PCR (qRT-PCR) experiments were carried out on TissueScan Cancer Survey Panels, which contained cDNA samples from various normal and cancer tissues covering eight different tumors (breast, colon, kidney, liver, lung, ovary, prostate and thyroid) from independent patients diagnosed at various clinical disease stages and selected from mixed ages and genders. Tissue cDNAs of each array were synthesized from high-quality total RNAs of pathologist-verified tissues, normalized and validated with β-actin in two sequential qPCR analyses, and provided with clinical information and QC data [25].
To quantitatively determine the relative amount of PRDM3/MECOM, PRDM8, PRDM10, PRDM12, PRDM16, and ZFPM2/FOG2 RNAs, qRT-PCR was performed [25]. Primers were designed using Primer3Plus [43] and specificity was verified with the BLAST program and through in-silico PCR analysis by the UCSC Genome Browser [44].
The selected sequences of oligonucleotides forward (F) and reverse (R) were: PRDM3/MECOM-F 5′-AGTGGCAGTGACCTGGAAAC-3′; PRDM3/MECOM-R 5′-ACCGCAGTCTGCTCCTCTAA-3′; PRDM10-F 5′-CAGCACATTCGAAAGAAGCA-3′; PRDM10-R 5′-GCGTTCGGTAGTCTGTCGTT-3′; PRDM12-F 5′-GGGAGTCCTTACGCAACCTT-3′; PRDM12-R 5′-TTCCATTGTGCCTCCACTCT-3′; PRDM16-F 5′-ATGATGGACAAGGCAAAACC-3′; PRDM16-R 5′-GATGTGGGAGGTAGCAGAGG-3′; ZFPM2/FOG2-F 5′-GACAGTGCCCATCAGATTTC-3′; ZFPM2/FOG2-R 5′-GGGCAGGAATTCTTC CATTTT-3′.
The amplification products were also analyzed by agarose gel electrophoresis [45]. Data were normalized with β-actin gene provided with arrays. The relative gene expression was calculated using the 2−ΔΔCt method [25].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/10/3250/s1.

Author Contributions

Conceived and designed the analysis: C.A., M.R., A.C. (Alfredo Ciccodicola), and A.C. (Amelia Casamassimi) Analyzed the data and produced the report: A.F. and A.S. Contributed to the data analysis: A.F., A.S., M.R., and A.C. (Amelia Casamassimi) Interpreted and validated the results: C.A., M.R., and A.C. (Amelia Casamassimi) Wrote the paper: M.R., P.G., M.B., C.A., A.C. (Alfredo Ciccodicola), and A.C. (Amelia Casamassimi) All authors read and approved the manuscript.

Funding

This work was supported in part by ordinary funds from the University of Campania “Luigi Vanvitelli”.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

FCFold change
FDRFalse discovery rate
TCGAThe Cancer Genome Atlas
ACCAdrenocortical carcinoma
BLCABladder cancer
BRCABreast cancer
CESCCervical squamous cell carcinoma and endocervical adenocarcinoma
CHOLCholangiocarcinoma
COADColon adenocarcinoma
DLBCLymphoid neoplasm diffuse large B-cell lymphoma
ESCAEsophageal carcinoma
GBMGlioblastoma
HNSCHead and neck squamous cell carcinoma
KICHKidney chromophobe carcinoma
KIRCKidney renal clear cell carcinoma
KIRPKidney renal papillary cell carcinoma
LAMLAcute myeloid leukemia
LIHCLiver hepatocarcinoma
LUADLung adenocarcinoma
LUSCLung squamous cell carcinoma
OVOvarian cancer
PAADPancreas adenocarcinoma
PCPGPheochromocytoma and paraganglioma
PRADProstate adenocarcinoma
READRectum adenocarcinoma
SARCSarcoma
SKCMSkin cutaneous melanoma
STADStomach adenocarcinoma
TGCTTesticular germ cell tumors
THCAThyroid cancer
THYMThymoma
UCECUterine corpus endometrial carcinoma
UCSUterine carcinosarcoma
UVMUveal melanoma

References

  1. Fumasoni, I.; Meani, N.; Rambaldi, D.; Scafetta, G.; Alcalay, M.; Ciccarelli, F.D. Family expansion and gene rearrangements contributed to the functional specialization of PRDM genes in vertebrates. BMC. Evol. Biol. 2007, 7, 187. [Google Scholar] [CrossRef] [PubMed]
  2. Sun, X.J.; Xu, P.F.; Zhou, T.; Hu, M.; Fu, C.T.; Zhang, Y.; Jin, Y.; Chen, Y.; Chen, S.J.; Huang, Q.H.; et al. Genome-wide survey and developmental expression mapping of zebrafish SET domain-containing genes. PLoS ONE 2008, 3, e1499. [Google Scholar] [CrossRef] [PubMed]
  3. Clifton, M.K.; Westman, B.J.; Thong, S.Y.; O’Connell, M.R.; Webster, M.W.; Shepherd, N.E.; Quinlan, K.G.; Crossley, M.; Blobel, G.A.; Mackay, J.P. The identification and structure of an N-terminal PR domain show that FOG1 is a member of the PRDM family of proteins. PLoS ONE 2014, 9, e106011. [Google Scholar] [CrossRef] [PubMed]
  4. Vervoort, M.; Meulemeester, D.; Béhague, J.; Kerner, P. Evolution of Prdm Genes in Animals: Insights from Comparative Genomics. Mol. Biol. Evol. 2016, 33, 679–696. [Google Scholar] [CrossRef] [PubMed]
  5. Hohenauer, T.; Moore, A.W. The Prdm family: Expanding roles in stem cells and development. Development 2012, 139, 2267–2282. [Google Scholar] [CrossRef] [PubMed]
  6. Fog, C.K.; Galli, G.G.; Lund, A.H. PRDM proteins: Important players in differentiation and disease. BioEssays 2012, 34, 50–60. [Google Scholar] [CrossRef] [PubMed]
  7. Di Zazzo, E.; De Rosa, C.; Abbondanza, C.; Moncharmont, B. PRDM Proteins: Molecular Mechanisms in Signal Transduction and Transcriptional Regulation. Biology 2013, 2, 107–141. [Google Scholar] [CrossRef] [PubMed]
  8. Mzoughi, S.; Tan, Y.X.; Low, D.; Guccione, E. The role of PRDMs in cancer: One family, two sides. Curr. Opin. Genet. Dev. 2016, 36, 83–91. [Google Scholar] [CrossRef] [PubMed]
  9. Huang, S.; Shao, G.; Liu, L. The PR domain of the Rb-binding zinc finger protein RIZ1 is a protein binding interface and is related to the SET domain functioning in chromatin-mediated gene expression. J. Biol. Chem. 1998, 273, 15933–15939. [Google Scholar] [CrossRef] [PubMed]
  10. Ren, B.; Chee, K.J.; Kim, T.H.; Maniatis, T. PRDI-BF1/Blimp-1 repression is mediated by corepressors of the Groucho family of proteins. Genes Dev. 1999, 13, 125–137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Kajimura, S.; Seale, P.; Kubota, K.; Lunsford, E.; Frangioni, J.V.; Gygi, S.P.; Spiegelman, B.M. Initiation of myoblast to brown fat switch by a PRDM16-C/EBP-beta transcriptional complex. Nature 2009, 460, 1154–1158. [Google Scholar] [CrossRef] [PubMed]
  12. Seale, P.; Bjork, B.; Yang, W.; Kajimura, S.; Chin, S.; Kuang, S.; Scimè, A.; Devarakonda, S.; Conroe, H.M.; Erdjument-Bromage, H.; et al. PRDM16 controls a brown fat/skeletal muscle switch. Nature 2008, 454, 961–967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Okashita, N.; Suwa, Y.; Nishimura, O.; Sakashita, N.; Kadota, M.; Nagamatsu, G.; Kawaguchi, M.; Kashida, H.; Nakajima, A.; Tachibana, M.; et al. PRDM14 Drives OCT3/4 Recruitment via Active Demethylation in the Transition from Primed to Naive Pluripotency. Stem Cell Rep. 2016, 7, 1072–1086. [Google Scholar] [CrossRef] [PubMed]
  14. Chi, J.; Cohen, P. The Multifaceted Roles of PRDM16: Adipose Biology and Beyond. Trends Endocrinol. Metab. 2016, 27, 11–23. [Google Scholar] [CrossRef] [PubMed]
  15. Sorrentino, A.; Rienzo, M.; Ciccodicola, A.; Casamassimi, A.; Abbondanza, C. Human PRDM2: Structure, function and pathophysiology. Biochim. Biophys. Acta 2018, 1861, 657–671. [Google Scholar] [CrossRef] [PubMed]
  16. Maruvka, Y.E.; Mouw, K.W.; Karlic, R.; Parasuraman, P.; Kamburov, A.; Polak, P.; Haradhvala, N.J.; Hess, J.M.; Rheinbay, E.; Brody, Y.; et al. Analysis of somatic microsatellite indels identifies driver events in human tumors. Nat. Biotechnol. 2017, 35, 951–959. [Google Scholar] [CrossRef] [PubMed]
  17. Choi, E.J.; Kim, M.S.; Song, S.Y.; Yoo, N.J.; Lee, S.H. Intratumoral Heterogeneity of Frameshift Mutations in MECOM Gene is Frequent in Colorectal Cancers with High Microsatellite Instability. Pathol. Oncol. Res. 2017, 23, 145–149. [Google Scholar] [CrossRef] [PubMed]
  18. Khurana, S.; Kruhlak, M.J.; Kim, J.; Tran, A.D.; Liu, J.; Nyswaner, K.; Shi, L.; Jailwala, P.; Sung, M.H.; Hakim, O.; et al. A macrohistone variant links dynamic chromatin compaction to BRCA1-dependent genome maintenance. Cell Rep. 2014, 8, 1049–1062. [Google Scholar] [CrossRef] [PubMed]
  19. Zhang, C.; Zhu, Q.; He, H.; Jiang, L.; Qiang, Q.; Hu, L.; Hu, G.; Jiang, Y.; Ding, X.; Lu, Y. RIZ1: A potential tumor suppressor in glioma. BMC Cancer 2015, 15, 990. [Google Scholar] [CrossRef] [PubMed]
  20. Kang, R.; Zelazowski, M.J.; Cole, F. Missing the Mark: PRDM9-Dependent Methylation Is Required for Meiotic DSB Targeting. Mol. Cell 2018, 69, 725–727. [Google Scholar] [CrossRef]
  21. Diagouraga, B.; Clément, J.A.J.; Duret, L.; Kadlec, J.; de Massy, B.; Baudat, F. PRDM9 Methyltransferase Activity Is Essential for Meiotic DNA Double-Strand Break Formation at Its Binding Sites. Mol. Cell 2018, 69, 853–865. [Google Scholar] [CrossRef] [PubMed]
  22. Shu, X.; Geng, H.; Li, L.; Ying, J.; Ma, C.; Wang, Y.; Poon, F.F.; Wang, X.; Ying, Y.; Yeo, W.; et al. The Epigenetic Modifier PRDM5 Functions as a Tumor Suppressor through Modulating WNT/β-Catenin Signaling and Is Frequently Silenced in Multiple Tumors. PLoS ONE 2011, 6, e27346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Meani, N.; Pezzimenti, F.; Deflorian, G.; Mione, M.; Alcalay, M. The tumor suppressor PRDM5 regulates Wnt signaling at early stages of zebrafish development. PLoS ONE 2009, 4, e4273. [Google Scholar] [CrossRef] [PubMed]
  24. Weinstein, J.N.; Collisson, E.A.; Mills, G.B.; Shaw, K.R.; Ozenberger, B.A.; Ellrott, K.; Shmulevich, I.; Sander, C.; Stuart, J.M. The Cancer Genome Atlas Pan-Cancer analysis project. Nat. Genet. 2013, 45, 1113–1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Federico, A.; Rienzo, M.; Abbondanza, C.; Costa, V.; Ciccodicola, A.; Casamassimi, A. Pan-Cancer Mutational and Transcriptional Analysis of the Integrator Complex. Int. J. Mol. Sci. 2017, 18, 936. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9, 40. [Google Scholar] [CrossRef] [PubMed]
  27. Vogelstein, B.; Papadopoulos, N.; Velculescu, V.E.; Zhou, S.; Diaz, L.A., Jr.; Kinzler, K.W. Cancer genome landscapes. Science 2013, 339, 1546–1558. [Google Scholar] [CrossRef] [PubMed]
  28. Garraway, L.A.; Lander, E.S. Lessons from the cancer genome. Cell 2013, 153, 17–37. [Google Scholar] [CrossRef] [PubMed]
  29. Marx, V. Cancer genomes: Discerning drivers from passengers. Nat. Methods 2014, 11, 375–379. [Google Scholar] [CrossRef] [PubMed]
  30. Porta-Pardo, E.; Kamburov, A.; Tamborero, D.; Pons, T.; Grases, D.; Valencia, A.; Lopez-Bigas, N.; Getz, G.; Godzik, A. Comparison of algorithms for the detection of cancer drivers at subgene resolution. Nat. Methods 2017, 14, 782–788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Tokheim, C.J.; Papadopoulos, N.; Kinzler, K.W.; Vogelstein, B.; Karchin, R. Evaluating the evaluation of cancer driver genes. Proc. Natl. Acad. Sci. USA 2016, 113, 14330–14335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Castro-Giner, F.; Ratcliffe, P.; Tomlinson, I. The mini-driver model of polygenic cancer evolution. Nat. Rev. Cancer 2015, 15, 680–685. [Google Scholar] [CrossRef] [PubMed]
  33. Tsang, A.P.; Visvader, J.E.; Turner, C.A.; Fujiwara, Y.; Yu, C.; Weiss, M.J.; Crossley, M.; Orkin, S.H. FOG, a multitype zinc finger protein, acts as a cofactor for transcription factor GATA-1 in erythroid and megakaryocytic differentiation. Cell 1997, 90, 109–119. [Google Scholar] [CrossRef]
  34. Cai, Q.; Jeannet, R.; Hua, W.K.; Cook, G.J.; Zhang, B.; Qi, J.; Liu, H.; Li, L.; Chen, C.C.; Marcucci, G.; et al. CBFβ-SMMHC creates aberrant megakaryocyte-erythroid progenitors prone to leukemia initiation in mice. Blood 2016, 128, 1503–1515. [Google Scholar] [CrossRef] [PubMed]
  35. Ruffalo, M.; Koyutürk, M.; Sharan, R. Network-Based Integration of Disparate Omic Data to Identify “Silent Players” in Cancer. PLoS Comput. Biol. 2015, 11, e1004595. [Google Scholar] [CrossRef] [PubMed]
  36. Huang, S.; Chaudhary, K.; Garmire, L.X. More Is Better: Recent Progress in Multi-Omics Data Integration Methods. Front. Genet. 2017, 8, 84. [Google Scholar] [CrossRef] [PubMed]
  37. Mayakonda, A.; Koeffler, H.P. Maftools: Efficient analysis, visualization and summarization of MAF files from large-scale cohort based cancer studies. bioRxiv 2016. [Google Scholar] [CrossRef]
  38. HGNC. Available online: http://www.genenames.org (accessed on 19 October 2018).
  39. Tamborero, D.; Gonzalez-Perez, A.; Lopez-Bigas, N. OncodriveCLUST: Exploiting the positional clustering of somatic mutations to identify cancer genes. Bioinformatics 2013, 29, 2238–2244. [Google Scholar] [CrossRef] [PubMed]
  40. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef] [PubMed]
  41. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef] [PubMed]
  42. TCGA. Available online: https://tcga-data.nci.nih.gov/tcga/ (accessed on 25 September 2017).
  43. Primer3Plus. Available online: http://primer3plus.com/cgi-bin/dev/primer3plus.cgi (accessed on 12 October 2017).
  44. UCSC-Genome Browser. Available online: https://genome.ucsc.edu (accessed on 8 January 2018).
  45. Rienzo, M.; Schiano, C.; Casamassimi, A.; Grimaldi, V.; Infante, T.; Napoli, C. Identification of valid reference housekeeping genes for gene expression analysis in tumor neovascularization studies. Clin. Transl. Oncol. 2013, 15, 211–218. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Stacked histograms showing the number of different classes of somatic mutations affecting PRDM genes as reported in the Mutation Annotation Files across all analyzed cancer entities.
Figure 1. Stacked histograms showing the number of different classes of somatic mutations affecting PRDM genes as reported in the Mutation Annotation Files across all analyzed cancer entities.
Ijms 19 03250 g001
Figure 2. The stacked histogram represents the amount of deleterious somatic mutations affecting the different known domains of PRDM proteins. In detail, the missense mutations affecting the PR domain are reported in blue, the missense mutations affecting the Zn fingers are in red, and the missense mutations affecting other known domains (where present) are illustrated in green. The deleterious missense mutations not affecting known domains are shown in violet and the other classes of deleterious mutations (frameshift, in-frame deletions, stop gained and start lost mutations, splice site, UTR, and intron variants) are in orange.
Figure 2. The stacked histogram represents the amount of deleterious somatic mutations affecting the different known domains of PRDM proteins. In detail, the missense mutations affecting the PR domain are reported in blue, the missense mutations affecting the Zn fingers are in red, and the missense mutations affecting other known domains (where present) are illustrated in green. The deleterious missense mutations not affecting known domains are shown in violet and the other classes of deleterious mutations (frameshift, in-frame deletions, stop gained and start lost mutations, splice site, UTR, and intron variants) are in orange.
Ijms 19 03250 g002
Figure 3. Mutations in PRDM genes recurrent in more than three tumor types.
Figure 3. Mutations in PRDM genes recurrent in more than three tumor types.
Ijms 19 03250 g003
Figure 4. (a) ZFPM1/FOG1 hotspot mutations in ACC obtained by Lollipop plot visualization function. These hotspot mutations are all localized outside the PR domain (green cylinder) and Zinc fingers (blue dots) (be) I-TASSER predicted tertiary structures of the annotated ZFPM1/FOG1 canonical protein (b) and ZFPM1/FOG1 proteins carrying the mutations ELPA444-447A (c), LA446-447P (d), and K431N (e); the arrows show the mutated protein regions. (f) The scatter plot shows the results of the OncodriveCLUST algorithm analysis for ACC. The dimension of the dots is proportional to the number of clusters found in a certain gene, also indicated in the squared bracket. Specifically, in the ZFPM1/FOG1 locus, two mutation clusters were found (fdr < 2.87 × 10−6).
Figure 4. (a) ZFPM1/FOG1 hotspot mutations in ACC obtained by Lollipop plot visualization function. These hotspot mutations are all localized outside the PR domain (green cylinder) and Zinc fingers (blue dots) (be) I-TASSER predicted tertiary structures of the annotated ZFPM1/FOG1 canonical protein (b) and ZFPM1/FOG1 proteins carrying the mutations ELPA444-447A (c), LA446-447P (d), and K431N (e); the arrows show the mutated protein regions. (f) The scatter plot shows the results of the OncodriveCLUST algorithm analysis for ACC. The dimension of the dots is proportional to the number of clusters found in a certain gene, also indicated in the squared bracket. Specifically, in the ZFPM1/FOG1 locus, two mutation clusters were found (fdr < 2.87 × 10−6).
Ijms 19 03250 g004
Figure 5. The heatmap shows the expression profiles of PRDMs across the analyzed cancer types.
Figure 5. The heatmap shows the expression profiles of PRDMs across the analyzed cancer types.
Ijms 19 03250 g005
Figure 6. Relative expressions obtained by real-time PCR in different types of cancer tissues versus corresponding normal tissues (with arbitrary expression value equal to 1). The comparative threshold cycle (Ct) method was used with β-actin as the internal control. The results are expressed as the mean ± ES. The statistical significance of differences between experimental groups was calculated using the unpaired two-tailed Student’s t-test. (*) Results with a p-value < 0.05 were considered significant. (a) ZFPM2/FOG2; (b) PRDM16; (c) PRDM3/MECOM; (d) PRDM10; (e) PRDM12. (f) Representative samples analyzed for PRDM12 by agarose gel electrophoresis.
Figure 6. Relative expressions obtained by real-time PCR in different types of cancer tissues versus corresponding normal tissues (with arbitrary expression value equal to 1). The comparative threshold cycle (Ct) method was used with β-actin as the internal control. The results are expressed as the mean ± ES. The statistical significance of differences between experimental groups was calculated using the unpaired two-tailed Student’s t-test. (*) Results with a p-value < 0.05 were considered significant. (a) ZFPM2/FOG2; (b) PRDM16; (c) PRDM3/MECOM; (d) PRDM10; (e) PRDM12. (f) Representative samples analyzed for PRDM12 by agarose gel electrophoresis.
Ijms 19 03250 g006
Table 1. Frequency of patients carrying mutations in the PRDMs across the 31 analyzed tumors.
Table 1. Frequency of patients carrying mutations in the PRDMs across the 31 analyzed tumors.
CancersGenes
PRDM1PRDM2MECOMPRDM4PRDM5PRDM6PRDM7PRDM8PRDM9PRDM10PRDM11PRDM12PRDM13PRDM14PRDM15PRDM16ZNF408ZFPM1ZFPM2
ACC01.101.11.12.2003.21.1001.1001.1050.54.3
BLCA1.74.15.33.43.401.913.92.21.701.41.72.74.10.713.6
BRCA0.80.910.50.30.20.20.50.70.80.20.20.30.50.60.50.10.21.3
CESC0.5421.520.510.54.511.5000.53310.51
CHOL2.7005.402.7005.400000002.700
COAD2.26.65.51.52.62.90.71.15.54.41.10.41.51.82.24.41.56.62.6
DLBC8.222020028.20000026.1004.1
ESCA1.11.63.22.11.62.71.11.674.81.602.12.74.82.71.60.58.6
GBM0.60.62.80.31.100.303.31.90.60.80.30.31.90.60.600.3
HNSC1.32.12.50.81.300.60.67.21.11.10.80.81.90.92.30.60.90.8
KICH01.51.500001.56001.51.5001.5000
KIRC0.20.81.21.50.30.20.80110.20.30.70.51.20.50.50.21.2
KIRP0.32.400.30.30.70.30.30.710.70.31.411.71.41.40.30.3
LAML000000000.50000000.5000
LIHC1.63.22.91.81.10.80.50.82.41.61.800.81.32.42.10.803.9
LUAD1.93.72.11.61.901.91.214.231.61.10.94.42.34.20.90.211.1
LUSC2.83.4502.800.607.33.42.21.13.402.23.41.705
OV0.20.400000.401.300.2000.20.20001.3
PAAD0.54.32.71.11.602.1161.63.22.72.71.60.511.21.61.102.7
PCPG00.50000001.100000000.500.5
PRAD00.81.400.600.60.60.80.40.20.40.20.40.40.80.60.61.4
READ02.61.72.63.40.9002.61.70.90.9001.70.909.45.1
SARC1.11.91.91.10.82.70.80.431.51.10.40.81.53.430.80.82.3
SKCM5.34.220.11.33.602.31.515.43.631.32.53.63.67.81.90.816.5
STAD2.37.83.813012.35.84.31.50.32.82.534.81.815.8
TGCT02.50.6001.30.61.92.52.50.601.3001.301.30
THCA0.40.600.40.200.80.41.20.80.20.20.200.40.80.40.40.4
THYM2.43.23.22.402.43.20.83.21.60.80.80.82.42.43.20.80.81.6
UCEC47.25.62.44.420.81.6105.22.41.22.84.83.25.640.44.8
UCS5.21.7003.401.71.71.71.71.701.703.401.75.25.2
UVM01.201.21.2001.24.901.20001.20000
Table 2. Percentage of deleterious and tolerated mutations in the PRDMs across the analyzed tumor samples.
Table 2. Percentage of deleterious and tolerated mutations in the PRDMs across the analyzed tumor samples.
GenesDeleterious MutationsTotal Mutations% Deleterious Mutations
PRDM112027244.1
PRDM222350744.0
MECOM/PRDM328854752.7
PRDM48816055.0
PRDM513624854.8
PRDM6315358.5
PRDM73111427.2
PRDM85415435.1
PRDM940389944.8
PRDM1018533555.2
PRDM1110118754.0
PRDM12307241.7
PRDM138917251.7
PRDM1412526347.5
PRDM1516234147.5
PRDM1625751450.0
ZNF408/PRDM175515535.5
ZFPM1/FOG14310242.2
ZFPM2/FOG223470033.4

Share and Cite

MDPI and ACS Style

Sorrentino, A.; Federico, A.; Rienzo, M.; Gazzerro, P.; Bifulco, M.; Ciccodicola, A.; Casamassimi, A.; Abbondanza, C. PR/SET Domain Family and Cancer: Novel Insights from The Cancer Genome Atlas. Int. J. Mol. Sci. 2018, 19, 3250. https://doi.org/10.3390/ijms19103250

AMA Style

Sorrentino A, Federico A, Rienzo M, Gazzerro P, Bifulco M, Ciccodicola A, Casamassimi A, Abbondanza C. PR/SET Domain Family and Cancer: Novel Insights from The Cancer Genome Atlas. International Journal of Molecular Sciences. 2018; 19(10):3250. https://doi.org/10.3390/ijms19103250

Chicago/Turabian Style

Sorrentino, Anna, Antonio Federico, Monica Rienzo, Patrizia Gazzerro, Maurizio Bifulco, Alfredo Ciccodicola, Amelia Casamassimi, and Ciro Abbondanza. 2018. "PR/SET Domain Family and Cancer: Novel Insights from The Cancer Genome Atlas" International Journal of Molecular Sciences 19, no. 10: 3250. https://doi.org/10.3390/ijms19103250

APA Style

Sorrentino, A., Federico, A., Rienzo, M., Gazzerro, P., Bifulco, M., Ciccodicola, A., Casamassimi, A., & Abbondanza, C. (2018). PR/SET Domain Family and Cancer: Novel Insights from The Cancer Genome Atlas. International Journal of Molecular Sciences, 19(10), 3250. https://doi.org/10.3390/ijms19103250

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop