Abstract
RIZ1 displays strong tumor-suppressive activities, which has a potential histone methyltransferase activity. The objective of the study was to evaluate the level and the methylation status of RIZ1 and analyze its association with clinicopathological features and the histone in the pituitary adenomas. We found that RIZ1-positive cases were 11/50 and H-Scores 22.75 ± 11.83 in invasive pituitary adenomas and 26/53 and 66.3 ± 21.7 in non-invasive pituitary adenomas (χ2 = 8.182, p = 0.004). RIZ1 and C-myc showed the opposite trend in these cases. The methylation levels of RIZ1 were more than 50% in 30.4% (7/23) CpG sites through MALDI-TOF Mass array. There was significant difference (p < 0.01) in 4 CpG sites between invasive pituitary adenoma group and non-invasive pituitary adenoma group; furthermore, the relieved methylation levels of H3K4/H3K9 and enhanced methylation levels of H3K27 in the patients’ serum were found. Furthermore, there was statistic difference of H3K4 and H3K27 methylation between invasive pituitary adenoma and non-invasive pituitary adenoma group (p < 0.01). The average progression-free survival in high RIZ1 group was 52.63 ± 7.62 months and 26.06 ± 4.23 months in low RIZ1 group (p < 0.05). Promoter region methylation of RIZ1 may play an important role in the epigenetic silencing of RIZ1 expression in pituitary adenomas, which may translate into important diagnostic and therapeutic applications.
Introduction
Pituitary adenomas (PAs) account for 10%–15% of all intracranial neoplasms, and these tumors are invariably benign and exhibit features of differentiated pituitary cell function as well as premature proliferative arrest. 1 The mechanisms underlying pituitary tumorigenesis are not very clear. Altered expression of pituitary cell-cycle genes, activation of pituitary selective oncoproteins, or loss of pituitary suppressor factors may be associated with aberrant growth factor signaling in PAs. 2 Pituitary cells are among the few epithelial cell types that rarely undergo malignant transformation although several studies have been performed to identify the genetic cause of pituitary tumourigenesis. 3
In recent years, emerging evidence indicates that epigenetic modifications including deoxyribonucleic acid (DNA) methylation, histone modification, nucleosome remodeling, and ribonucleic acid (RNA)-mediated targeting, are an alternative force altering the expression level of genes involved in neoplastic development, including pituitary tumorigenesis.4,5 DNA methylation is a stable modification that leads to chromatin remodeling, resulting in transcriptional silencing without gene mutation. 6 Aberrant methylation patterns in DNA damage repair genes may serve as predictive, diagnostic, prognostic, and chemosensitive markers of human cancer. 7 Increasing evidence suggests, however, that most methylation changes are generated in a programmed manner and occur in a subpopulation of tissue cells during normal aging, probably predisposing them for tumorigenesis. 8 DNA methylation–mediated gene silencing is closely linked to the deacetylation and methylation of histones at key lysine residues. 9
Members of the PRDM family are characterized by an N terminal PR domain that is related to the SET methyltransferase domain and multiple zinc fingers that mediate sequence-specific DNA binding and protein–protein interactions. 10 RIZ1, also named PRDM2, is originally isolated as a molecule that is associated with the retinoblastoma protein Rb. 11 RIZ1 displays methyltransferase activity toward H3K9. 12 PRDM factors either act as direct histone methyltransferases or recruit a suite of histone-modifying enzymes to target promoters. 13
The comprehensive analysis of many histone modifications in the human genome demonstrates that the combination of several histone modified forms a modular pattern and influences transcriptional activations in a cooperative manner.14,15 The signals of methylations of H3K9 and H4K20 are high in silent genes and repressed chromatin domains. 16 H3K9 methylation can be mediated by several histone methyltransferases (HMTases). 12 Methylation of H3K27 is likely to spread over larger regions around transcription start sites of silent genes. 17
In this study, we carried out series of immunohistochemistry (IHC) to assess a scatter of data on RIZ1 expression in a large cohort of PAs and analyze the promoter methylation of RIZ1 and histone methylation in PAs.
Material and methods
Patients and tissue specimens
The files of the Department of Neurosurgery, the First Affiliated Hospital of Zhengzhou University, from January 2010 to July 2014 were searched for 103 PAs. Histologically, the tumors were classified into 23 prolactinomas (PRL), 23 GHomas, and 57 non-functional pituitary adenoma (NFPA) according to the 2007 World Health Organization histological classification of tumors of the PAs. 18 The following IPA diagnostic criteria were adopted: (1) Knosp classification grade IV tumors and Hardy classification invasive adenomas, (2) tumor cells confirmed via pathology as invading sellar bone or adjacent dura mater, and (3) tumor cells invading the sphenoid sinus cavity or peripheral vascular and nerve. Three normal pituitary specimens were from donors. The study was conducted according to an Institutional Review Board (IRB)–approved protocol (IRB of The First Affiliated Hospital of Zhengzhou University), and written informed consent was obtained from all patients for surgery and for research purposes (Table 1).
The clinicopathological data in 103 cases.
PRL: prolactinomas; NFPA: non-functional pituitary adenoma; SD: standard deviation; FSH: follicle-stimulating hormone; LH: luteinizing hormone; TSH: thyroid-stimulating hormone.
IHC
All slides were evaluated in advance using an H&E stain to assess tumor content and quality. The slides were placed in the Leica BOND-III purchased from Leica Biosystems, which is a fully automated, random, and continuous-access slide-staining system that processes IHC tests simultaneously. IHC protocol F was selected in the machine, and 3 min with ER1 (epitope retrieval) was set for heat-induced epitope retrieval (HIER) parameter. Bond™ Polymer Refine Detection (DS9800, Leica Biosystems) was used for detection of primary antibodies. The slides were scanned into digital pictures and expression was examined using Aperio AT2 (Leica Biosystems). Primary antibodies, anti-RIZ1 (1:1000, ab3791; Abcam) and C-myc (1:500, ab17767; Abcam) was used. The optimal titer of the primary Abs for the remainder was determined based on pre-experiment results. The results were calculated using Aperio AT2 (Leica Biosystems) with digital slide-viewing software.
The staining intensity was stratified on a scale of 0–3 (0 for no staining, 1 for weak, 2 for moderate, and 3 for strong). The H-score is obtained multiplying by the intensity of its stain and a constant adjusting the mean for the strongest staining (score = 1.0× (% weak) + 2.0× (% medium) + 3.0× (%strong)).
DNA extraction, purification, and bisulfite modification
Genomic DNA (gDNA) from frozen tissues was extracted by using a DNeasy kit (QIAGEN). All extracted gDNA was treated with sodium bisulfite (Sigma). Briefly, 2 µg gDNA was denatured by 5.5 µL of 3 mol fresh NaOH (final concentration 0.3 mol/L) for 10 min at 37°C; 30 µL of 10 mmol/L hydroquinone (Sigma) and 520 µL of 3 mol/L sodium bisulfite (pH 5.0) were added kept in the dark. The mixture was inverted, added to 200 µL liquid paraffin to prevent water evaporation and reagent oxidation, and then incubated at 50°C for 16 h. The modified DNA was purified using the Wizard DNA Clean-Up System (Promega). The purified DNA was treated again with NaHSO3 and precipitated. DNA was re-dissolved in 20 µL of TAE, 2 µL of which were subjected to polymerase chain reaction (PCR) amplification.
Methylation-specific PCR and sequencing
Three CpG islands were found using the EMBOSS Cpgplot. Methylation-specific primers were designed to cover 30 CpG dinucleotides numbered 3445 upstream of the transcription start site G (forward) and 4024 upstream of the transcription start site G (reverse) according equenom® EpiDesigner (#24). Primers specific for methylated DNA (forward 5′-ATTTAGGTTGGAGAGTAGTGGGGTA′; reverse 5′-TAATCCCAACACTTTAAAAAACCAA-3′) were added to the reaction and expected to generate 580 bp products. PCR conditions were 45 cycles of denaturation at 94°C for 20 s, annealing at 60°C for 30 s and 72°C for 1 min. Then, sequence PCR products with Sequenom® MassARRAY and specimens’ methylation were generated by the EpiTYPER software v1.0 (Sequenom).
Histone extraction and methyltransferase assays
Specimens were cut into small pieces and added TEB buffer (PBS containing 0.5% Triton X-100, 2 mM phenylmethane sulfonyl fluoride (PMSF), and 0.02% NaN3). The homogenized mixture was centrifuged at 3000 r/min for 5 min at 4°C. The pellet was resuspended with extraction buffer (0.5 N HCl, 10% glycerol), incubated on ice for 30 min and then centrifuged at 12,000 r/min for 5 min at 4°C. The supernatant was collected into a new vial, 8 volume of acetone was added, and incubated at −20°C overnight. Again centrifugation was carried out at 12,000 r/min for 5 min and the pellet was air-dried; the pellet was dissolved in distilled water and the protein concentration was quantified.
Histone methyltransferase activity was measured according to the Epigentek’s guide in https://www.epigentek.com/catalog/ (P-3024, P-3003, and P-3044).
Statistics
χ2, log rank analysis, and Fisher’s exact tests were used to determine the significance of categorical variables. One-way analysis of variance (ANOVA) test was applied to the examination of differential expression of methylation of H3K4, H3K9, and H3K27. All p values are two sided and 0.05 was applied as the significance level.
Results
Sample characteristics of the cohort
We selected 103 consecutive patients diagnosed as PAs who had primary tumor specimens. For 103 samples, sufficient archival material was present for standard H&E staining and IHC staining using anti-RIZ1 and C-myc antibody, including 54 men and 49 women with a median age of 45.1 years (range = 13–72 years). The tumor size ranged from 1.3 to 6.1 cm in greatest dimension (median size: 4.3 cm). Additional sample characteristics are given in Table 1. According to the Knosp classification and the intraoperation, all cases were divided into invasive pituitary adenoma (IPA) group (50 cases) and non-invasive pituitary adenoma (nIPA) group (53 cases).
The headache (45.6%, 47/103) is the most initial clinic symptom, and visual disturbance (38.9%, 40/103) is a common clinic symptom. The headache of IPA was 32/56 and 15/47 in nIPA (χ2 = 6.556, p = 0.01). The recurrence of IPA was 12/56 and 5/47 in nIPA (χ2 = 2.80, p = 0.094).
The RIZ1 and C-myc expression in PAs
IHC was conducted to investigate the expression pattern of RIZ1 and C-myc in IPA and nIPA. Immunoreactivity was observed primarily in cytoplasm of tumor cells, and IHC staining for RIZ1 and C-myc in representative samples of PAs tissues was shown in Figure 1. In this study, the positive case was defined as more than 10% of positive cells.

The expression of RIZ1 and c-myc in 103 cases: (a) Results of IHC (bar = 60 µm) and (b) results of western blot.
RIZ1-positive cases were 11/50 in IPA group and 26/53 in nIPA group (Table 2; χ2 = 8.182, p = 0.004). H-Scores were 22.75 ± 11.83 in IPA group and 66.3 ± 21.7 in nIPA group in Table 2 (p < 0.01). RIZ1 overexpression in malignant meningioma cells was associated with the downregulation of C-myc expression which played a role in cell-cycle progression, apoptosis, and cellular transformation. 19 In this study, C-myc-positive cases were 41/50 and 42/53, respectively. H-Scores were 265 ± 27.41 and 170 ± 41.72 as shown in Table 2 (t = 2.934 p < 0.01).
The RIZ1 and C-myc levels in 103 cases.
IPA: invasive pituitary adenoma; nIPA: non-invasive pituitary adenoma.
The RIZ1 methylation level in PAs
Epigenetic activation or inactivation of genes plays a critical role in many important human diseases, especially in cancer. A major mechanism of epigenetic inactivation of genes is methylation of CpG islands in genome DNA caused by DNA methyltransferases. We investigated the RIZ1 promoter methylation status in 12 PA samples and 3 normal pituitary samples with Sequenom’s EpiTYPER assay, and 23 CpG sites per sample were analyzed (total 245 CpG sites in 15 samples, Figure 2). About 30.4% (7/23) CpG sites of RIZ1 in PAs samples were defined as having a high-degree methylation compared to normal pituitary, with mean methylation levels above 50%. Furthermore, significant difference between methylation pattern of IPA group and nIPA group was detected in 4 CpG sites (8, 11, 12, and 27; p < 0.01).

The methylation CpG sites of RIZ1 promoter. Yellow clusters indicate 0% methylated, blue clusters indicate 100% methylated, color gradient between yellow and blue indicates methylation ranging from 0 to 100, red clusters indicate 0% methylated, and white clusters indicate CpG sites not analyzed. Highlight shows the hypermethylation CpG site.
The Methylation level of H3K4, H3K9 and H3K27 in PAs
Histone methyltransferases (HMTs) control or regulate DNA methylation through chromatin-dependent transcription repression or activation. Lysine methylation is associated to transcription activation or repression depending on the residue and degree of methylation. 20
To further understand the molecular function of RIZ1, we sought to detect the methylation of H3K4, H3K9, and H3K27 in 42 PAs samples and 3 normal pituitary samples. H3K4 mono-methylation was associated with silenced gene repression. Results of enzyme-linked immunosorbent assay (ELISA) showed the relieved mono-methylation H3K4 expression in IPA group (p < 0.01) but not in nIPA group (Figure 3). H3K9 methylation mediates heterochromatin formation by forming a binding site for HP1 and also participates in silencing gene expression at euchromatic sites. The methyltransferase activity of H3K9 in PAs was relieved compared to the normal pituitary (p < 0.01), and there was no statistic difference between IPA group and nIPA group. Increased global H3K27 methylation was also found be involved in some pathological processes such as cancer progress. There was enhanced global methylation H3K27 in IPA group (p < 0.01) compared to normal pituitary and nIPA group.

The methylation level of H3K4, H3K9, and H3K27 in normal pituitary and PA samples.
RIZ1 level and PAs’ recurrence and prognostic
Follow-up periods ranged from 2 to 8 years (mean = 4.8 years). The median expression level was used as the cutoff. There was no significant correlation between RIZ1 levels and age, gender, tumor size, or hormone serum levels. However, low RIZ1 levels were more frequently observed in IPA (χ2 = 9.356, p = 0.002) or recurrent tumors (χ2 = 5.500, p = 0.037) in Table 3. Log rank analysis showed that the average progression-free survival (PFS) in high RIZ1 group (4/51, 52.63 months) was longer than that in low RIZ1 group (13/52, 26.06 months; p < 0.05; Figure 4).
Relationship between RIZ1 H-Scores in 103 patients and various clinical parameters.
The median expression level was used as the cutoff: low RIZ1 H-Scores were defined as values below the 50th percentile of the 103 patients; values at or above 50th percentile were classified as high levels.

The average progression-free survival in 103 PA cases.
Discussion
Since the identification of RIZ1 as an RB-binding protein, several studies reported that the genomic region containing the RIZ1 gene (chromosome 1p36) is frequently deleted or rearranged in human tumors.11,21,22 In this study, we examined the association of RIZ1 expression with the clinical feature in 103 PA cases through tissue microarray (TMA) and found RIZ1 was associated with an increased risk of recurrence for adenomas.
The RIZ1 gene is a member of a superfamily of histone/protein methyltransferases. Buyse et al. first found RIZ1 by retinoblastoma (Rb) probes in 1995, while performing a functional screen for Rb, and observed that it was able to interact with Rb. 11 The gene is commonly inactivated in human cancers. Gene knock-out study has established RIZ1 as a tumor susceptibility gene in mice. The gene has potent tumor-suppressive activities in causing apoptosis, G2/M arrest, or both. 23 RIZ1 had a direct repressor function on C-myc gene expression. RIZ1 overexpression in malignant meningioma cells was associated with the downregulation of C-myc expression. 19 Moreover, C-myc, putative downstream targets of H3K9 methylation, may be involved in regulating RIZ1 tumor-suppressive effects. The reciprocal relationship between RIZ1 and C-myc was then validated in primary meningioma cells and human tumor samples. 24 In this study, the RIZ1-positive rate in nIPA is 22% and 49.1% in IPA group. In addition, the RIZ1 level was approximately two-fold lower in IPA group than in nIPA group through H-Scores (p < 0.01). Binary multivariate regression analysis revealed that low RIZ1 expression level was independently associated with PAs recurrence. The average PFS in high RIZ1 group was longer than that in low RIZ1 group (p < 0.05).
According to a genetic perspective, RIZ1 may be downregulated by chromosomal instability and microsatellite instability, as well as frame-shift mutations, point mutations, and heterozygote deficiency.25 –27 From an epigenetic perspective, the deactivation of RIZ1 may occur due to promoter methylation and histone acetylation. The RIZ1 promoter has been demonstrated to have the characteristics of a CpG island, which suggests that RIZ1 is a target of inactivation by epigenetic mechanisms. 28 In all, 42.6% prostate cancer cases and 69% gastric adenocarcinomas cases were reported to have RIZ1 methylation, and RIZ1 methylation was more frequent in patients with a high-grade malignancy.29,30 In this study, we explored the methylation status of RIZ1 among normal pituitary, nIPA, and IPA specimens using MALDI-TOF MS. After methylation quantification of the RIZ1, our results showed hypomethylation patterns for PAs compared with normal pituitary in 7 CpG sites; furthermore, there was statistically significant difference in 4 CpG sites between IPA group and nIPA group.
Dynamic changes in histone post-translational modifications can alter DNA-transcription factor interactions during many nuclear processes and may either accompany or precede transcriptional activation or repression. Transcription activation correlates with methylation of H3K4, together with histone acetylation (H3K9Ac). 31 Clinical data indicate that methylation of H3K4 may be elevated in breast, kidney, and colon cancers and correlates with a poor clinical outcome.32,33 By contrast, transcription repression often involves tri-methylation of lysine 27 of H3 (H3K27me3) and di- or tri-methylation of lysine 9 of H3 (H3K9me2/3), through the recruitment of repressive protein complexes. RIZ1 associates with >4400 promoters in G0 myoblasts, 55% of which are also marked with H3K9me2 and enriched for myogenic, cell cycle, and developmental regulators. Knockdown of RIZ1 alters histone methylation at key promoters such as Myogenin and CyclinA2 (CCNA2). 34 Epigenetic regulation by RIZ1 preserves key functions of the quiescent state, with implications for stem-cell self-renewal.
In this study, relieved methylation of H3K4 and enriched methylation of H3K27 were observed in IPA group and relieved methylation of H3K9 was observed between IPA group and nIPA group. Changes in expression of key histone modifications are found in PAs cases, which would be expected because tumor cells require instant changes in gene expression and chromatin structure in order to promote cell proliferation and tumor cell invasion. The results of the survival analyses of the individual markers reflect our expected results based on the cellular functions of RIZ1 and the respective histone modifications. We also found that nuclear expression of methylation on H3K4 and H4K27 has prognostic value in IPA cases. Methylation of H3K4 and H3K27 involves modifications found on gene promoter regions and is associated with activation of gene transcription, and lower expression of H3K4 and higher expression of H3K27 in tumors could lead to aberrant gene transcription, including genes required for cell survival, proliferation, and migration. For methylation of H3K9, literature shows conflicting results with respect to patient survival and prognosis, depending on the type of cancer. 35 Histone modification of H3K9 is associated with silencing of gene transcription. 36 In this study, there is no statistic difference between IPA group and nIPA group although the lower methylation of H3K9 in PAs.
Somatic gene mutations are uncommon in sporadic PAs, and recent studies lend support to epigenetic modification as a potential cause of tumorigenesis and tumor progression. We found that RIZ1 was silenced through DNA methylation and showed particular susceptibility to epigenetic modification, with abnormal DNA methylation in >50% of PA samples. And H3K4 and H3K27 methylation may play an important role regulating the tumor invasion.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by First Affiliated Hospital of Zhengzhou University Research and Development Funds.
