Glutathione S-transferase pi 1 (GSTP1) is lowly expressed in normal prostate luminal cells and becomes induced in most proliferative inflammatory atrophy (PIA) lesions. GSTP1 becomes silenced in prostatic intraepithelial neoplasia (PIN) and prostate adenocarcinoma (CaP) via cytosine-phospho-guanine (CpG) island promoter hypermethylation. However, GSTP1 methylation patterns in PIA and PIN, and their relationship to patterns in CaP are poorly understood. We used bisulfite genomic sequencing to examine patterns of GSTP1 promoter CpG island methylation in laser capture microdissected benign, PIA, PIN, and CaP regions from 32 subjects that underwent radical prostatectomy. We analyzed 908 sequence clones across 24 normal epithelium, 37 PIA, 18 PIN, and 23 CaP regions, allowing assessment of 34,863 CpG sites with allelic phasing. Normal and PIA lesions were mostly unmethylated with 0.52 and 1.3% of total CpG sites methylated, respectively. PIN and CaP lesions had greater methylation with 24% and 51% of total CpG sites methylated, respectively. The degree of GSTP1 methylation showed progression from PIA << PIN < CaP. PIN lesions showed more partial methylation compared with CaP lesions. Partially methylated lesions were enriched for methylation changes at AP1 and SP1 transcription factor binding sites. These results demonstrate that methylation density in the GSTP1 CpG island in PIN was intermediate relative to that in normal prostate epithelium/PIA and CaP lesions. These results are consistent with gradual spreading of DNA methylation centered at the SP1/AP1 transcription factor binding sites in precursor lesions, with subsequent spreading of methylation across the entire CpG island in transition to CaP.

Prevention Relevance:

DNA hypermethylation at the GSTP1 promoter progressively spreads from being unmethylated in normal prostate to intermediate levels in precursor lesions to extensive methylation in cancer. This molecular progression of GSTP1 promoter methylation patterns in early prostate carcinogenesis could be useful for identification and interception of prostate cancer precursors.

Initiation of cancer involves multiple genetic events characterized by chromosomal translocations, deletions, amplifications, and point mutations of critical genes (1). In addition to these genetic changes, abnormal epigenetic alterations, including changes in patterns of DNA methylation, have been observed in a wide spectrum of human cancers, including widespread genomic hypomethylation and simultaneous regional increases in DNA methylation (2, 3). This aberrant somatic DNA methylation of cytosine-phospho-guanine (CpG) dinucleotides within upstream regulatory regions of genes is generally associated with a compacted chromatin structure, and accompanying repression of transcription (4). A promising area of strength for DNA methylation in the clinic is in the area of molecular diagnostics and early detection (5). The hypermethylation of promoter CpG islands of tumor-suppressor genes is a common alteration in cancer cells and leads to the transcriptional inactivation and loss of normal cellular functions (6). In many cases hypermethylation of the CpG island in gene regulatory regions has been correlated with a loss of gene expression; thus, it is likely that DNA methylation provides an alternate pathway to gene deletion or mutation for the loss of tumor suppressor gene function.

A large body of evidence suggests that the major precursor lesion to invasive prostatic adenocarcinoma (CaP) is high-grade prostatic intraepithelial neoplasia (HGPIN; refs. 7–10). This hypothesis is supported by findings of proximity to CaP, morphological changes and molecular alterations shared between the two lesions, and two-fold increased prevalence and extent of HGPIN in prostates with carcinoma vs. those without. In addition, regions of atrophic epithelial cells associated with inflammatory infiltrates, known as proliferative inflammatory atrophy (PIA), frequently merge with HGPIN lesions and harbor some of the hallmark genome alterations in prostate cancer and prostatic intraepithelial neoplasia (PIN; refs. 11, 12). Thus, PIA lesions have been proposed as risk factor lesions and potential intermediates towards prostate cancer development (7, 11, 13). Given limited understanding of the specific patterns of methylation from benign to invasive prostate cancer and its intermediate steps, understanding the patterns of DNA methylation at CpG islands in both PIA and PIN lesions is of importance.

The CpG island promoter region spanning the glutathione S-transferase (GSTP1) gene is methylated in approximately 90% of prostatic adenocarcinomas. GSTP1 encodes an enzyme involved in detoxification and protection from oxidants and carcinogens by conjugation to glutathione (14). GSTP1 is expressed at high levels in normal Basal cells and is expressed at very low levels in normal luminal cells. However, it appears highly induced in PIA, in which many cells express high levels of GSTP1, apparently in response to increased chemical stress. In contrast, HGPIN and CaP are often devoid of GSTP1 protein, and often harbor methylation of the CpG island at the GSTP1 promoter (4, 10, 11, 15, 16).

As hypermethylation in the 5′ region and promoter of GSTP1 appears to be a frequent and early event in prostate cancer development, we wished to examine the specific changes in the pattern of methylation spanning the gene promoter in CaP and candidate precursor lesions, HGPIN and PIA. We used bisulfite genomic sequencing to determine the methylation status of 39 cytosine residues present spanning the GSTP1 gene promoter. Specifically, we wished to address whether the hypermethylation (1) proceeded via intermediate molecular steps (2), was enriched at specific regions in the promoter, such as transcription factor binding sites, and (3) presented with inter- or intra-tumor heterogeneity across lesions within patients.

Prostate tissue samples

Frozen or formalin-fixed, paraffin-embedded (FFPE) radical prostatectomy specimens were selected from 32 patients undergoing radical prostatectomy at The Johns Hopkins Hospital in which PIA and PIN lesions were readily discernible histologically using hematoxylin and eosin (H&E) staining and encompassed regions large enough to be amenable to laser capture microdissection (LCM). No subjects were excluded from our study, and male subjects were eligible for the study if they were undergoing radical prostatectomy for nonmetastatic prostate cancer. Also, to avoid the potential that PIN lesions represented intraductal/intra-acinar spread of carcinoma, most of them were isolated from tissue blocks either not containing carcinoma or in which the PIN lesions were “away” from carcinoma by at least 2 mm. Also, because HGPIN can be difficult to diagnose on frozen sections, most of the PIN lesions analyzed were from FFPE sections (12/18). The PIN lesions that we did use from frozen tissues were clearly diagnosable as HGPIN on those sections. The use of tissues was approved by the Johns Hopkins University School of Medicine Institutional Review Board under a waiver of consent and performed in accordance with the ethical guidelines set for by the U.S. Common Rule.

LCM and isolation of genomic DNA

Three serial 6-μm sections were stained with hematoxylin. Matched regions of normal epithelium (4 of the regions of normal were from transition zone in histological BPH nodules), PIA, HGPIN, and CaP were obtained by LCM using the Veritas system (Arcturus Engineering Inc., Mountain View, CA) and the CapSure HS LCM Caps (Arcturus Engineering Inc.). Approximately 1,000 to 3,000 cells were obtained in most cases by LCM. DNA was isolated with a standard phenol/chloroform protocol. DNA quantification was carried out by real-time PCR of the Beta-globin gene to ensure that DNA was of ample quality for reliable quantitative PCR amplification (11). Some of the sections were stained with antibodies against GSTP1 for IHC.

Bisulfite genomic sequencing

Two nanograms of genomic DNA was bisulfite converted using the EZ DNA methylation kit (Zymo Research, Orange, CA) and eluted in 11 μL of TE buffer, pH 7.4. Portions (1 μL) of the DNA were amplified by PCR using primers that would bind to a CpG-free sequence from the MYOD1 gene only after complete bisulfite conversion of all cytosines, thus allowing assessment of the efficacy of bisulfite conversion. MYOD1 primer set was 5′-CCAACTCCAAATCCCCTCTCTAT-3′ (forward), and 5′-TGATTAATTTAGATTGGGTTTAGAGAAGGA-3′ (reverse). 9 μL of the bisulfite converted DNA was subjected to nested PCR for GSTP1 alleles. The outer PCR was carried out in 40 μL reactions containing 10x Platinum Taq buffer (Invitrogen, Carlsbad, CA), 1.5 U Platinum Taq (Invitrogen), 250 μmol/L each dNTPs, 1.5 mmol/L MgCl2, 0.25 μg/μL BSA, 2 μL dimethyl sulfoxide, 400 nmol/L forward primer and 400 nmol/L reverse primer. Cycling conditions were 95°C for 3 minutes, 40 cycles of 95°C for 30 seconds, 50°C for 45 seconds and 72°C for 45 seconds, followed by a 7-minute extension step at 72°C. The outer PCR products amplified using primer #1 and #2 were purified and concentrated using the DNA Clean and Concentrator 5 kit (Zymo Research, Orange, CA) and eluted in 20 μL of TE buffer. Aliquots (4 μL) of the cleaned PCR product used in the inner PCR amplification using the #3 and #4 primers in a 40 μL reaction mixture for 20 cycles. Primers amplifying GSTP1 CpG islands without bias to methylation patterns were #1 forward primer, 5′-GTTGGTTTTATGTTGGGAGTTTTGAG-3′; #2 reverse primer, 5′-ATCCTCTTCCTACTATCTATTTACTCCCTAA-3′; #3 forward primer, 5′-GTTGGGAGTTTTGAGTTTTATTTT-3′; and #4 reverse primer, 5′-ACTATCTATTTACTCCCTAAACCCC-3′. Nested PCR products were gel purified and cloned into pCR2.1-TOPO vectors (Invitrogen, Carlsbad, CA). Plasmid DNA from 10 positive clones was sequenced using an ABI 3730 DNA Analyzer (Applied Biosystems, Foster City, CA). Bisulfite sequencing data were analyzed using a custom Java program called DNAMethylMap, which facilitates the analysis of Sanger bisulfite sequencing clones with virtually bisulfite converted reference amplicon sequences (17). Clones that showed conversion of >95% of cytosines outside of a CpG context were considered valid. DNA from the prostate cancer cell line LNCaP (RRID:CVCL_0395), which is known to be completely methylated at all CpGs in the region analyzed (18), was used as a positive control for methylation, and normal human WBC DNA (Novagen, Madison, WI) was used as a negative control.

Statistical assessment

The overall extent of GSTP1 gene promoter methylation across the four tissue types (normal, PIA, HGPIN, and CaP) were assessed with the Mann–Whitney U test. The extent of GSTP1 CpG dinucleotide methylation on an individual clone/allele across all 39 CpGs analyzed for that allele was classified as negative (<10%), mild (10%–25%), moderate (25%–50%) and high (>50%). Differences in frequency of GSTP1 hypermethylation between normal, atrophy, HGPIN, and CaP were assessed with the χ2 test, Fisher exact test, and Mann–Whitney U test. For further interactive analysis and visualization of the data, we used a custom Java program called DNAMethylMap v1.2 (17).

Data availability

The data generated in this study are available upon request from the corresponding authors.

LCM and bisulfite sequencing of the GSTP1 promoter CpG island

The GSTP1 gene is approximately 4 kb in length, comprises 7 exons and 6 introns (refs. 19, 20; Fig. 1A). Primers amplifying the GSTP1 CpG island without bias to methylation patterns allowed investigation into the methylation status of 39 CpG dinucleotides between −250 to +81 upstream of the GSTP1 promoter region [GenBank accession number AY324387, (GRCh38, Chr11:67,583,561–67,583,896); Fig. 1A].

To obtain genomic DNA from prostate tissues enriched for epithelial cells from normal prostate, PIA lesions, HGPIN, and prostate adenocarcinoma (CaP), LCM was undertaken using radical prostatectomy specimens from 102 lesions across 32 patients (Supplementary Table S1); representative images of PIA and PIN lesions before and after LCM are shown in Fig. 1B. Using nested amplification of bisulfite-treated DNA followed by cloning of individual PCR products, sequencing for 5-meC versus C was performed on individual clones/alleles from GSTP1 promoter amplicons. A total of 212 clones from normal epithelium (24 LCM regions), 327 clones from PIA (37 LCM regions), 167 clones from HGPIN (18 LCM regions), and 202 clones from CaP (23 LCM regions) were sequenced, giving a sum total of 34,863 CpG sites analyzed. The promoter region of GSTP1 in the LNCaP human CaP cell line was almost completely methylated (5,644/5,770 or 98%), while white blood cells showed virtually no CpG methylation (24/5,670 or 0.42%), as expected.

GSTP1 promoter CpG island methylation patterns in normal glands, PIA, PIN, and invasive cancer lesions

Overall, a gradual increase in the average methylation level was evident across normal, PIA, PIN, and cancer regions (Fig. 1C). The difference in average percent methylation of the GSTP1 promoter was significant between normal prostate and HGPIN (P < 0.0001; Mann–Whitney U test), and between PIA and HGPIN (P < 0.0001; Mann–Whitney U test; Fig. 1C).

We next examined the patterns of methylation at each individual allele in each tissue type (Fig. 2). Normal prostate epithelium was almost completely unmethylated (42/8,079 CpGs methylated or 0.52%; Fig. 2A), and all clones had less than 10% methylation (Supplementary Fig. S1). In PIA lesions, very sparse GSTP1 promoter methylation was evident (164/12,521 or 1.31%; Fig. 2B) with only 7 of 37 lesions, representing 18 of 327 clones, showing any evidence of methylation (Supplementary Tables S2 and S3). One of the PIA lesions showed one clone with moderate methylation. In addition, 6 of 7 PIA lesions with >10% methylated CpGs were either next to or merging with PIN or atypical glands suspicious for, but not diagnostic of, cancer (e.g., often referred to as “atyp” or atypical small acinar proliferation [ASAP]), compared with 10 of 30 PIA lesions with <10% methylated CpGs (P = 0.0118; χ2 test). The pattern of methylation of GSTP1 alleles in PIA lesions showed only sparse methylation, usually having less dense methylation changes compared with PIN and CaP lesions (Fig. 2).

HGPIN lesions showed an intermediate level of GSTP1 promoter methylation between normal/PIA and CaP. HGPIN had 1,534 of 6,461 (24%) CpGs methylated (Fig. 2C). 28.7% of PIN alleles had a mild or moderate level of methylation that was significantly different from the 5.5% in PIA (P = 6.9E-13; χ2 test) and 11.4% in cancer (P = 2.6E-05; χ2 test). Many individual alleles in HGPIN showed an intermediate density of methylation compared with normal/PIA and to CaP (Fig. 2C). CaP lesions showed an increase in CpG methylation level and density over that in HGPIN, with 4,002 of 7,802 (51.3%) CpG sites showing methylation and 124 of 202 individual GSTP1 alleles (61.4%) exhibiting moderate and high density of methylation. The increase in highly methylated GSTP1 promoter alleles in CaP versus HGPIN was significant (P = 3.5E-06; χ2 test).

GSTP1 promoter CpG sites were more commonly methylated near the transcription start site (OR = 3.94; P < 2.2E-16; Fisher test) and transcription factor binding sites for AP1 and SP1 (OR = 2.04; P = 7.8E-12; Fisher test) compared with methylation at CpG sites in the rest of the CpG island (Fig. 2; Supplementary Fig. S2). Six proximally located CpG sites (between −4 and +38 of the transcription start site) could be considered “hot spots” for somatic methylation changes, appearing methylated in GSTP1 with either sparse, intermediate, or dense methylation patterns from atrophy, HGPIN, and CaP (Fig. 2). There was also a smaller upstream region (−197 to −176) with recurrently methylated CpGs not correlating with any known transcription factor binding sites.

Allele level intra- and inter-patient heterogeneity in GSTP1 promoter CpG island methylation

To highlight the patterns of intra- and inter-individual heterogeneity, we carefully examined four cases in which genomic DNA was isolated from laser capture microdissected normal, PIA, PIN, and invasive cancer regions (Figs. 3,45). In all four of these patients, we observed very minimal methylation in all individual alleles in both the normal and PIA lesions. In comparing patients 21964 and 1618, both cases had highly dense methylation across all the GSTP1 alleles in CaP and showed intermediate GSTP1 methylation in the PIN lesions (Fig. 3). However, the PIN lesion in patient 21964 had alleles with heterogeneity in the GSTP1 promoter methylation patterns (Fig. 3A), but patient 1618 had nearly identical methylation across all 10 clones from the PIN lesion (Fig. 3B).

In examining the GSTP1 alleles for CaP lesions in Patient 28126, an averaged methylation view signifies an intermediate level of methylation (Fig. 4). However, there were clearly two distinct groups of alleles, 5 clones with nearly complete methylation and the other four without any methylation (Fig. 4). Unlike the previous two patients (Fig. 3), the PIN lesion from this subject had a similar methylation pattern to the CaP lesion (Fig. 4). This pattern of having two distinct fully methylated or unmethylated allele subsets was observed in a total of 3 of 23 CaP lesions.

In patient 26559, it is notable that three different cancer lesions harbored three distinct methylation patterns (Fig. 5). One cancer lesion had clones with full methylation across the GSTP1 promoter, the next one had clones with intermediate levels of methylation, and the last cancer lesion had no methylation across all of the clones. This highlights the degree of heterogeneity there can be with the GSTP1 CpG island promoter hypermethylation patterns within different tumor sites in a given patient.

The expression of GSTP1 protein as detected by IHC in patient 21964 (Fig. 6) demonstrates that lack of methylation in the GSTP1 promoter (Fig. 3A) corresponds with GSTP1 expression in normal and PIA lesions (Fig. 6AD). In addition, the presence of hypermethylation (Fig. 3A) corresponds with silencing of GSTP1 in luminal cells, with GSTP1 expression present only in stromal cells in the cancer lesion (Fig. 6G and H). The PIN allele-level distribution shows 3 alleles without and 6 alleles with methylation (Fig. 3A), likely corresponding with the distinct regions of GSTP1-positive and GSTP1-negative expression respectively (Fig. 6E and F). In a different patient, a similar pattern of GSTP1 expression heterogeneity corresponding with GSTP1 methylation heterogeneity is observed in an HGPIN lesion and in a PIA lesion, in which a few regions of low levels of GSTP1 expression are demonstrated (Supplementary Fig. S3).

Somatic CpG island hypermethylation, which results in gene silencing, accompanies human prostate cancer pathogenesis, and is often evident in early stages of disease development, including cancer precursor lesions (4). Intraepithelial neoplasia refers to epithelial lesions in many organ sites that exhibit cytological features of carcinoma but fail to invade through basement membrane structures and are thought to represent precursors to cancer, and have been proposed to represent conditions suitable for therapeutic intervention in pursuit of cancer prevention or interception (21). In the prostate, HGPIN is an intraepithelial neoplasia lesion that is traditionally thought to serve as a precursor to many CaP lesions. In this study, the levels and density of GSTP1 CpG island hypermethylation in HGPIN lesions were often intermediate between normal prostate epithelium and CaP. This observation provides molecular evidence that HGPIN lesions often represent an intermediate step in the pathogenesis of CaP. Interestingly, the levels of GSTP1 CpG island hypermethylation in PIA lesions showed a trend to be higher than normal, but significantly less dense than HGPIN.

Methylation-specific primers (MSP), and related techniques, are used to detect methylation of CpG dinucleotides at specific genome sites. MSP assays have reported a limit of detection below 0.1% methylated alleles (22). Using MSP, CpG methylation levels are directly assessed only at primer annealing sites, so using different MSP primer sets, the frequency of GSTP1 promoter methylation levels in different prostate lesions from different prostate cancer cases has appeared variable throughout different studies (11). The different primer sets are prone to false negative results depending on the distribution of methylated CpGs across the GSTP1 promoter (4, 11, 18). The bisulfite genomic sequencing assay strategy used in this study was capable of providing the pattern and phasing of 5-meCpG along individual alleles. In this strategy, bisulfite treatment of genomic DNA is combined with the sequencing of individual product clones generated by PCR amplification (23). The method also permits monitoring of conversion efficiency for unmethylated cytosine bases; PCR product clones from a bisulfite genomic sequencing assay can be discarded if inadequately converted. The major limitation of the bisulfite genomic sequencing approach appears to be DNA fragmentation during the bisulfite reaction, with as much as 90% ± 6% of input DNA subject to some level of degradation during bisulfite treatment (24). Hence, although it is possible to amplify sequences from as few as 10 copies of bisulfite-treated genomes in theory, the consistency and quality of the bisulfite genomic sequencing results tend to be markedly better if more DNA is used. For our study, we collected as much as 2 ng of genomic DNA and were successfully able to recover bisulfite converted GSTP1 promoter sequences as PCR clones. However, even though we sequenced 10 clones per sample for the study, the PCR products could have been derived from as few as 1 or 2 DNA molecules, resulting in stochastic clonal amplification (25). Thus, observing identical methylation patterns across clones (for example, the PIN lesion from patient 28126, Fig. 3B) may either indicate a true homogeneity in the methylation pattern in that sample, or could be due to clonal amplification of a single allele. Despite this limitation, the insights generated by examining the density of CpG methylation on individual GSTP1 alleles were very informative: among all GSTP1 alleles sequenced, the density of methylation in genomic DNA from HGPIN was clearly intermediate between normal prostate/PIA and CaP.

GSTP1 methylation serves a significant role as a candidate biomarker for prostate cancer. In particular, GSTP1 methylation is assessed with methylation-specific PCR (MSP), as described above, in the ConfirmMDX assay as part of a 3-gene panel in detection of occult prostate cancer on patient's previously-biopsied, negative tissue (26). On the basis of the MATLOC trial that demonstrated clinical utility of MSP-based methylation of the GSTP1 promoter to detect occult prostate cancer in histopathologically negative biopsies, the primer sets in these assays are likely targeting CpGs around 65–85 and 105–125 bases upstream of the transcription start site (27, 28).

Given the patterns of methylated CpGs across the GSTP1 promoter observed in this study in those regions, it is possible that certain CaP lesions are missed due to lack of dense methylation at those sites in an MSP-based assay. In addition, given the methylation distribution immediately upstream to the AP1 binding site seen in the PIN lesions in this study, an MSP-based assay may detect PIN in the absence of cancer being present and hence would result in false positives for certain patient samples. However, these PIN lesions may also be missed, fortuitously so, due to the limitations of the MSP-based assay, similar to the situation for cancer lesions described above.

Recent evidence suggests that at times, lesions morphologically identifiable as HGPIN may also arise as retrograde spread from invasive carcinomas that invade otherwise benign glands and mimic the morphology of HGPIN (29, 30). We refer to such lesions as post-invasive intraepithelial carcinomas (PIC). This is analogous to intraductal carcinoma of the prostate, which represents an aggressive form of prostate cancer that is presumed to arise through such retrograde invasion of invasive cancer into preexisting benign glandular structures. Therefore, it is important to develop methods to determine whether a given lesion with the appearance morphologically of HGPIN, represents de novo HGPIN or represents PIC mimicking PIN. We speculate that DNA methylation changes in cancer cells offer the potential to distinguish whether HGPIN lesions represent CaP precursors or CaP itself (e.g., PIC), because somatic CpG dinucleotide methylation patterns at specific genome sites can be sparse or dense, allowing the stepwise molecular progression from initially unmethylated (normal prostate) to ultimately extensively methylated (CaP) alleles, transiting through intermediate methylation in precursor lesions. In this study, while we did not provide any direct evidence of such cases, it is possible that because CaP cells can grow along established prostate glands and duct structures, some fraction of HGPIN lesions could in fact represent CaP mimicking intraepithelial neoplasia (29, 30). The identical GSTP1 promoter CpG island methylation patterns between CaP and HGPIN lesions in the same patient (Fig. 4) has the possibility to represent such cases of retrograde spread of invasive carcinomas into prostate glands, or it may represent a case in which the PIN had “progressed” in methylation and gave rise to this invasive carcinoma. Future studies from our group and others will work to address these unanswered questions by performing comprehensive analyses on such samples with whole-genome sequencing, transcriptome sequencing, and genome-wide methylation sequencing techniques. This will allow characterization of the clonal relationships of such lesions and assess whether PIC lesions more commonly harbor the type of dense methylation associated with invasive cancer compared with bona fide PIN precursor lesions.

Among GSTP1 alleles from prostate cancer precursor lesions subjected to bisulfite genomic sequencing, an apparent “hot spot” for CpG methylation surrounding the transcription start site and the AP1/SP1 binding sites was identified. In an analysis of de novo CpG methylation in transfected GSTP1 promoters in prostate cancer cells, Song and colleagues observed a spread of cytosine methylation with serial passage of the cells if (i) an initial 5-meCpG was present to “seed” the allele, and (ii) transcription from the allele was silenced (31). In their experiments, CpGs immediately downstream of the transcription start site (and near the “hot spot” seen in the current study) most commonly showed de novo methylation. A seed and spread model may be able to account for the apparent increase in the CpG methylation density in GSTP1 alleles during prostatic carcinogenesis. Rare GSTP1 CpG dinucleotides might become methylated spontaneously or in response to stress in normal prostate epithelial cells/atrophy cells, and then additional methylated CpG dinucleotides might accumulate in the setting of transcriptional silencing. Perhaps, CpG methylation at the “hot spot” itself contributes to inhibition of transcription, serving to drive both seeding and spreading of a wave of cytosine methylation. The intrapatient heterogeneity of GSTP1 promoter CpG island methylation of CaP lesions (Fig. 5), may represent three events in the spectrum of such a seed and spread model. The observed cancer lesion devoid of GSTP1 promoter CpG methylation could represent an early cancer event in which no seeding event had yet occurred, while the cancer with intermediate levels of methylation contains cells soon after initial seeding and spread of CpG methylation around the “hotspot” area (Fig. 5). The cancer lesion with complete GSTP1 promoter methylation may indicate complete progression once the CpG methylation has spread throughout the entire GSTP1 promoter CpG island. The findings from this study raise new questions as to whether other densely methylated genes in prostate cancer show similar “seeding and spreading” patterns in prostate cancer precursor lesions.

W.G. Nelson reports grants from NIH during the conduct of the study. A.M. De Marzo reports grants from Janssen; personal fees from Cepheid; and personal fees from Merck Inc. outside the submitted work. S. Yegnasubramanian reports grants from NIH, Department of Defense CDMRP PCRP, Irving A. Hansen Memorial Foundation, Commonwealth Foundation, Bristol-Myers Squibb; and grants from Prostate Cancer Foundation during the conduct of the study; grants and personal fees from Cepheid; other support from Digital Harmonic, Brahm Astra Therapeutics; and grants from Janssen outside the submitted work; in addition, S. Yegnasubramanian has a patent for Epigenetic Tests for Prostate Cancer pending, a patent for COMPARE-MS issued and licensed to ClonTech, a patent for Agents for Reversing Epigenetic Silencing of Genes issued, and a patent for DNA Methylation Markers for Metastatic Prostate Cancer issued. No disclosures were reported by the other authors.

H. Gupta: Conceptualization, data curation, formal analysis, investigation, visualization, methodology, writing–original draft, writing–review and editing. H. Inoue: Conceptualization, data curation, formal analysis, investigation, visualization, methodology, writing–review and editing. Y. Nakai: Investigation, writing–review and editing. M. Nakayama: Investigation, writing–review and editing. T. Jones: Investigation, writing–review and editing. J.L. Hicks: Investigation, writing–review and editing. B. Kumar: Software, visualization, writing–review and editing. M. Gurel: Investigation, writing–review and editing. W.G. Nelson: Conceptualization, resources, funding acquisition, writing–original draft, writing–review and editing. A.M. De Marzo: Conceptualization, resources, data curation, formal analysis, supervision, funding acquisition, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing. S. Yegnasubramanian: Conceptualization, resources, data curation, software, formal analysis, supervision, funding acquisition, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing.

This work was funded in part by NIH/NCI grants U54CA274370 (to A.M. De Marzo, S. Yegnasubramanian), U01CA196390 (to A.M. De Marzo, S. Yegnasubramanian), P50CA058236 (to W.G. Nelson, A.M. De Marzo, S. Yegnasubramanian), P30CA006973 (to W.G. Nelson, A.M. De Marzo, S. Yegnasubramanian), the DOD CDMRP PCRP grant W81XWH-21–1-0295 (to W.G. Nelson, A.M. De Marzo, S. Yegnasubramanian), the Irving A. Hansen Memorial Foundation (to S. Yegnasubramanian), the Commonwealth Foundation for Cancer Research Foundation (to W.G. Nelson, S. Yegnasubramanian), and the Prostate Cancer Foundation (to W.G. Nelson, A.M. De Marzo, S. Yegnasubramanian). We thank Christina DeStefano Shields and Brenna Hairston for critical reading of the manuscript.

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Cancer Prevention Research Online (http://cancerprevres.aacrjournals.org/).

1.
Knudson
AG
.
Genetics of human cancer
.
Annu Rev Genet
1986
;
20
:
231
51
.
2.
Jones
PA
,
Rideout
WM
3rd
,
Shen
JC
,
Spruck
CH
,
Tsai
YC
.
Methylation, mutation, and cancer
.
Bioessays
1992
;
14
:
33
6
.
3.
Schmutte
C
,
Yang
AS
,
Nguyen
TT
,
Beart
RW
,
Jones
PA
.
Mechanisms for the involvement of DNA methylation in colon carcinogenesis
.
Cancer Res
1996
;
56
:
2375
81
.
4.
Nelson
WG
,
De Marzo
AM
,
Yegnasubramanian
S
.
Minireview: epigenetic alterations in human prostate cancers
.
Endocrinology
2009
;
150
:
3991
4002
.
5.
Laird
PW
.
The power and the promise of DNA methylation markers
.
Nat Rev Cancer
2003
;
3
:
253
66
.
6.
Esteller
M
.
Cancer epigenomics: DNA methylomes and histone-modification maps
.
Nat Rev Genet
2007
;
8
:
286
98
.
7.
De Marzo
AM
,
Platz
EA
,
Sutcliffe
S
,
Xu
J
,
Grönberg
H
,
Drake
CG
, et al
.
Inflammation in prostate carcinogenesis
.
Nat Rev Cancer
2007
;
7
:
256
69
.
8.
Brawer
MK
.
Androgen deprivation therapy: a cornerstone in the treatment of advanced prostate cancer
.
Rev Urol
2004
;
6
:
S3
9
.
9.
Bostwick
DG
,
Cheng
L
.
Precursors of prostate cancer
.
Histopathology
2012
;
60
:
4
27
.
10.
Trabzonlu
L
,
Kulac
I
,
Zheng
Q
,
Hicks
JL
,
Haffner
MC
,
Nelson
WG
, et al
.
Molecular pathology of high-grade prostatic intraepithelial neoplasia: challenges and opportunities
.
Cold Spring Harb Perspect Med
2019
;
9
:
a030403
.
11.
Nakayama
M
,
Bennett
CJ
,
Hicks
JL
,
Epstein
JI
,
Platz
EA
,
Nelson
WG
, et al
.
Hypermethylation of the human glutathione S-transferase-pi gene (GSTP1) CpG island is present in a subset of proliferative inflammatory atrophy lesions but not in normal or hyperplastic epithelium of the prostate: a detailed study using laser-capture microdissection
.
Am J Pathol
2003
;
163
:
923
33
.
12.
Shrestha
E
,
Coulter
JB
,
Guzman
W
,
Ozbek
B
,
Hess
MM
,
Mummert
L
, et al
.
Oncogenic gene fusions in nonneoplastic precursors as evidence that bacterial infection can initiate prostate cancer
.
Proc Natl Acad Sci USA
2021
;
118
:
e2018976118
.
13.
Sfanos
KS
,
De Marzo
AM
.
Prostate cancer and inflammation: the evidence
.
Histopathology
2012
;
60
:
199
215
.
14.
Rushmore
TH
,
Pickett
CB
.
Glutathione S-transferases, structure, regulation, and therapeutic implications
.
J Biol Chem
1993
;
268
:
11475
8
.
15.
Brooks
JD
,
Weinstein
M
,
Lin
X
,
Sun
Y
,
Pin
SS
,
Bova
GS
, et al
.
CG island methylation changes near the GSTP1 gene in prostatic intraepithelial neoplasia
.
Cancer Epidemiol Biomarkers Prev
1998
;
7
:
531
6
.
16.
Lin
X
,
Asgari
K
,
Putzi
MJ
,
Gage
WR
,
Yu
X
,
Cornblatt
BS
, et al
.
Reversal of GSTP1 CpG island hypermethylation and reactivation of pi-class glutathione S-transferase (GSTP1) expression in human prostate cancer cells by treatment with procainamide
.
Cancer Res
2001
;
61
:
8611
6
.
17.
Esopi
D
,
Graham
MK
,
Brosnan-Cashman
JA
,
Meyers
J
,
Vaghasia
A
,
Gupta
A
, et al
.
Pervasive promoter hypermethylation of silenced TERT alleles in human cancers
.
Cell Oncol
2020
;
43
:
847
61
.
18.
Yegnasubramanian
S
,
Lin
X
,
Haffner
MC
,
DeMarzo
AM
,
Nelson
WG
.
Combination of methylated-DNA precipitation and methylation-sensitive restriction enzymes (COMPARE-MS) for the rapid, sensitive, and quantitative detection of DNA methylation
.
Nucleic Acids Res
2006
;
34
:
e19
.
19.
Cowell
IG
,
Dixon
KH
,
Pemble
SE
,
Ketterer
B
,
Taylor
JB
.
The structure of the human glutathione S-transferase π gene
.
Biochem J
1988
;
255
:
79
83
.
20.
Morrow
CS
,
Cowan
KH
,
Goldsmith
ME
.
Structure of the human genomic glutathione S-transferase-pi gene
.
Gene
1989
;
75
:
3
11
.
21.
O'Shaughnessy
JA
,
Kelloff
GJ
,
Gordon
GB
,
Dannenberg
AJ
,
Hong
WK
,
Fabian
CJ
, et al
.
Treatment and prevention of intraepithelial neoplasia: an important target for accelerated new agent development
.
Clin Cancer Res
2002
;
8
:
314
46
.
22.
Herman
JG
,
Graff
JR
,
Myöhänen
S
,
Nelkin
BD
,
Baylin
SB
.
Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands
.
Proc Natl Acad Sci USA
1996
;
93
:
9821
6
.
23.
Frommer
M
,
McDonald
LE
,
Millar
DS
,
Collis
CM
,
Watt
F
,
Grigg
GW
, et al
.
A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands
.
Proc Natl Acad Sci USA
1992
;
89
:
1827
31
.
24.
Grunau
C
,
Clark
SJ
,
Rosenthal
A
.
Bisulfite genomic sequencing: systematic investigation of critical experimental parameters
.
Nucleic Acids Res
2001
;
29
:
E65
5
.
25.
Warnecke
PM
,
Mann
JR
,
Frommer
M
,
Clark
SJ
.
Bisulfite sequencing in preimplantation embryos: DNA methylation profile of the upstream region of the mouse imprintedH19Gene
.
Genomics
1998
;
51
:
182
90
.
26.
Partin
AW
,
Van Neste
L
,
Klein
EA
,
Marks
LS
,
Gee
JR
,
Troyer
DA
, et al
.
Clinical validation of an epigenetic assay to predict negative histopathological results in repeat prostate biopsies
.
J Urol
2014
;
192
:
1081
7
.
27.
Stewart
GD
,
Van Neste
L
,
Delvenne
P
,
Delrée
P
,
Delga
A
,
McNeill
SA
, et al
.
Clinical utility of an epigenetic assay to detect occult prostate cancer in histopathologically negative biopsies: results of the MATLOC study
.
J Urol
2013
;
189
:
1110
6
.
28.
Van Neste
L
,
Bigley
J
,
Toll
A
,
Otto
G
,
Clark
J
,
Delrée
P
, et al
.
A tissue biopsy-based epigenetic multiplex PCR assay for prostate cancer detection
.
BMC Urol
2012
;
12
:
16
.
29.
De Marzo
AM
,
Haffner
MC
,
Lotan
TL
,
Yegnasubramanian
S
,
Nelson
WG
.
Premalignancy in prostate cancer: rethinking what we know
.
Cancer Prev Res
2016
;
9
:
648
56
.
30.
Haffner
MC
,
Weier
C
,
Xu
MM
,
Vaghasia
A
,
Gürel
B
,
Gümüşkaya
B
, et al
.
Molecular evidence that invasive adenocarcinoma can mimic prostatic intraepithelial neoplasia (PIN) and intraductal carcinoma through retrograde glandular colonization
.
J Pathol
2016
;
238
:
31
41
.
31.
Song
JZ
,
Stirzaker
C
,
Harrison
J
,
Melki
JR
,
Clark
SJ
.
Hypermethylation trigger of the glutathione-S-transferase gene (GSTP1) in prostate cancer cells
.
Oncogene
2002
;
21
:
1048
61
.