Comprehensive Molecular Portraits of Invasive Lobular Breast Cancer.

Invasive lobular carcinoma (ILC) is the second most prevalent histologic subtype of invasive breast cancer. Here, we comprehensively profiled 817 breast tumors, including 127 ILC, 490 ductal (IDC), and 88 mixed IDC/ILC. Besides E-cadherin loss, the best known ILC genetic hallmark, we identified mutations targeting PTEN, TBX3, and FOXA1 as ILC enriched features. PTEN loss associated with increased AKT phosphorylation, which was highest in ILC among all breast cancer subtypes. Spatially clustered FOXA1 mutations correlated with increased FOXA1 expression and activity. Conversely, GATA3 mutations and high expression characterized luminal A IDC, suggesting differential modulation of ER activity in ILC and IDC. Proliferation and immune-related signatures determined three ILC transcriptional subtypes associated with survival differences. Mixed IDC/ILC cases were molecularly classified as ILC-like and IDC-like revealing no true hybrid features. This multidimensional molecular atlas sheds new light on the genetic bases of ILC and provides potential clinical options.

Radiogenomic Associations Clear Cell Renal Cell Carcinoma: An Exploratory Study.

INTRODUCTION: This study investigates how quantitative texture analysis can be used to non-invasively identify novel radiogenomic correlations with clear cell renal cell carcinoma (ccRCC) biomarkers. METHODS: The Cancer Genome Atlas-Kidney Renal Clear Cell Carcinoma open-source database was used to identify 190 sets of patient genomic data that had corresponding multiphase contrast-enhanced CT images in The Cancer Imaging Archive. 2,824 radiomic features spanning fifteen texture families were extracted from CT images using a custom-built MATLAB software package. Robust radiomic features with strong inter-scanner reproducibility were selected. Random forest, AdaBoost, and elastic net machine learning (ML) algorithms evaluated the ability of the selected radiomic features to predict the presence of 12 clinically relevant molecular biomarkers identified from the literature. ML analysis was repeated with cases stratified by stage (I/II vs. III/IV) and grade (1/2 vs. 3/4). 10-fold cross validation was used to evaluate model performance. RESULTS: Before stratification by tumor grade and stage, radiomics predicted the presence of several biomarkers with weak discrimination (AUC 0.60-0.68). Once stratified, radiomics predicted KDM5C, SETD2, PBRM1, and mTOR mutation status with acceptable to excellent predictive discrimination (AUC ranges from 0.70 to 0.86). CONCLUSIONS: Radiomic texture analysis can potentially identify a variety of clinically relevant biomarkers in patients with ccRCC and may have a prognostic implication.

Florida-California Cancer Research, Education and Engagement (CaRE2) Health Equity Center: Structure, Innovations, and Initial Outcomes.

INTRODUCTION: The Florida-California Cancer Research, Education, and Engagement (CaRE2) Health Equity Center is a triad partnership committed to increasing institutional capacity for cancer disparity research, the diversity of the cancer workforce, and community empowerment. This article provides an overview of the structure, process innovations, and initial outcomes from the first 4 years of the CaRE2 triad partnership. METHODS: CaRE2 serves diverse populations in Florida and California using a "molecule to the community and back" model. We prioritize research on the complex intersection of biological, environmental, and social determinants health, working together with scientific and health disparities communities, sharing expertise across institutions, bidirectional training, and community outreach. Partnership progress and outcomes were assessed using mixed methods and four Program Steering Committee meetings. RESULTS: Research capacity was increased through development of a Living Repository of 81 cancer model systems from minority patients for novel cancer drug development. CaRE2 funded 15 scientific projects resulting in 38 publications. Workforce diversity entailed supporting 94 cancer trainees (92 URM) and 34 ESIs (32 URM) who coauthored 313 CaRE2-related publications and received 48 grants. Community empowerment was promoted via outreaching to more than 3000 individuals, training 145 community cancer advocates (including 28 Community Scientist Advocates), and publishing 10 community reports. CaRE2 members and trainees together have published 639 articles, received 61 grants, and 57 awards. CONCLUSION: The CaRE2 partnership has achieved its initial aims. Infrastructure for translational cancer research was expanded at one partner institution, and cancer disparities research was expanded at the two cancer centers.

TENET 2.0: Identification of key transcriptional regulators and enhancers in lung adenocarcinoma.

Lung cancer is the leading cause of cancer-related death and lung adenocarcinoma is its most common subtype. Although genetic alterations have been identified as drivers in subsets of lung adenocarcinoma, they do not fully explain tumor development. Epigenetic alterations have been implicated in the pathogenesis of tumors. To identify epigenetic alterations driving lung adenocarcinoma, we used an improved version of the Tracing Enhancer Networks using Epigenetic Traits method (TENET 2.0) in primary normal lung and lung adenocarcinoma cells. We found over 32,000 enhancers that appear differentially activated between normal lung and lung adenocarcinoma. Among the identified transcriptional regulators inactivated in lung adenocarcinoma vs. normal lung, NKX2-1 was linked to a large number of silenced enhancers. Among the activated transcriptional regulators identified, CENPA, FOXM1, and MYBL2 were linked to numerous cancer-specific enhancers. High expression of CENPA, FOXM1, and MYBL2 is particularly observed in a subgroup of lung adenocarcinomas and is associated with poor patient survival. Notably, CENPA, FOXM1, and MYBL2 are also key regulators of cancer-specific enhancers in breast adenocarcinoma of the basal subtype, but they are associated with distinct sets of activated enhancers. We identified individual lung adenocarcinoma enhancers linked to CENPA, FOXM1, or MYBL2 that were associated with poor patient survival. Knockdown experiments of FOXM1 and MYBL2 suggest that these factors regulate genes involved in controlling cell cycle progression and cell division. For example, we found that expression of TK1, a potential target gene of a MYBL2-linked enhancer, is associated with poor patient survival. Identification and characterization of key transcriptional regulators and associated enhancers in lung adenocarcinoma provides important insights into the deregulation of lung adenocarcinoma epigenomes, highlighting novel potential targets for clinical intervention.

CT-based radiomics stratification of tumor grade and TNM stage of clear cell renal cell carcinoma.

OBJECTIVES: To evaluate the utility of CT-based radiomics signatures in discriminating low-grade (grades 1-2) clear cell renal cell carcinomas (ccRCC) from high-grade (grades 3-4) and low TNM stage (stages I-II) ccRCC from high TNM stage (stages III-IV). METHODS: A total of 587 subjects (mean age 60.2 years ± 12.2; range 22-88.7 years) with ccRCC were included. A total of 255 tumors were high grade and 153 were high stage. For each subject, one dominant tumor was delineated as the region of interest (ROI). Our institutional radiomics pipeline was then used to extract 2824 radiomics features across 12 texture families from the manually segmented volumes of interest. Separate iterations of the machine learning models using all extracted features (full model) as well as only a subset of previously identified robust metrics (robust model) were developed. Variable of importance (VOI) analysis was performed using the out-of-bag Gini index to identify the top 10 radiomics metrics driving each classifier. Model performance was reported using area under the receiver operating curve (AUC). RESULTS: The highest AUC to distinguish between low- and high-grade ccRCC was 0.70 (95% CI 0.62-0.78) and the highest AUC to distinguish between low- and high-stage ccRCC was 0.80 (95% CI 0.74-0.86). Comparable AUCs of 0.73 (95% CI 0.65-0.8) and 0.77 (95% CI 0.7-0.84) were reported using the robust model for grade and stage classification, respectively. VOI analysis revealed the importance of neighborhood operation-based methods, including GLCM, GLDM, and GLRLM, in driving the performance of the robust models for both grade and stage classification. CONCLUSION: Post-validation, CT-based radiomics signatures may prove to be useful tools to assess ccRCC grade and stage and could potentially add to current prognostic models. Multiphase CT-based radiomics signatures have potential to serve as a non-invasive stratification schema for distinguishing between low- and high-grade as well as low- and high-stage ccRCC. KEY POINTS: • Radiomics signatures derived from clinical multiphase CT images were able to stratify low- from high-grade ccRCC, with an AUC of 0.70 (95% CI 0.62-0.78). • Radiomics signatures derived from multiphase CT images yielded discriminative power to stratify low from high TNM stage in ccRCC, with an AUC of 0.80 (95% CI 0.74-0.86). • Models created using only robust radiomics features achieved comparable AUCs of 0.73 (95% CI 0.65-0.80) and 0.77 (95% CI 0.70-0.84) to the model with all radiomics features in classifying ccRCC grade and stage, respectively.

ZFX acts as a transcriptional activator in multiple types of human tumors by binding downstream of transcription start sites at the majority of CpG island promoters.

High expression of the transcription factor ZFX is correlated with proliferation, tumorigenesis, and patient survival in multiple types of human cancers. However, the mechanism by which ZFX influences transcriptional regulation has not been determined. We performed ChIP-seq in four cancer cell lines (representing kidney, colon, prostate, and breast cancers) to identify ZFX binding sites throughout the human genome. We identified ~9,000 ZFX binding sites and found that the majority of the sites are in CpG island promoters. Moreover, genes with promoters bound by ZFX are expressed at higher levels than genes with promoters not bound by ZFX. To determine if ZFX contributes to regulation of the promoters to which it is bound, we performed RNA-seq analysis after knockdown of ZFX by siRNA in prostate and breast cancer cells. Many genes with promoters bound by ZFX were downregulated upon ZFX knockdown, supporting the hypothesis that ZFX acts as a transcriptional activator. Surprisingly, ZFX binds at +240 bp downstream of the TSS of the responsive promoters. Using Nucleosome Occupancy and Methylome Sequencing (NOMe-seq), we show that ZFX binds between the open chromatin region at the TSS and the first downstream nucleosome, suggesting that ZFX may play a critical role in promoter architecture. We have also shown that a closely related zinc finger protein ZNF711 has a similar binding pattern at CpG island promoters, but ZNF711 may play a subordinate role to ZFX. This functional characterization of ZFX provides important new insights into transcription, chromatin structure, and the regulation of the cancer transcriptome.

Epigenetic plasticity potentiates a rapid cyclical shift to and from an aggressive cancer phenotype.

Highly tumorigenic, drug-resistant cancer stem-like cells drive cancer progression. These aggressive cells can arise repeatedly from bulk tumor cells independently of mutational events, suggesting an epigenetic mechanism. To test this possibility, we studied bladder cancer cells as they cyclically shifted to and from a cancer stem-like phenotype, and we discovered that these two states exhibit distinct DNA methylation and chromatin accessibility. Most differential chromatin accessibility was independent of methylation and affected the expression of driver genes such as E2F3, a cell cycle regulator associated with aggressive bladder cancer. Cancer stem-like cells exhibited increased E2F3 promoter accessibility and increased E2F3 expression that drove cell migration, invasiveness and drug resistance. Epigenetic interference using a DNA methylation inhibitor blocked the transition to a cancer stem-like state and reduced E2F3 expression. Our findings indicate that epigenetic plasticity plays a key role in the transition to and from an aggressive, drug-resistant phenotype.

Cell-type-specific enrichment of risk-associated regulatory elements at ovarian cancer susceptibility loci.

Understanding the regulatory landscape of the human genome is a central question in complex trait genetics. Most single-nucleotide polymorphisms (SNPs) associated with cancer risk lie in non-protein-coding regions, implicating regulatory DNA elements as functional targets of susceptibility variants. Here, we describe genome-wide annotation of regions of open chromatin and histone modification in fallopian tube and ovarian surface epithelial cells (FTSECs, OSECs), the debated cellular origins of high-grade serous ovarian cancers (HGSOCs) and in endometriosis epithelial cells (EECs), the likely precursor of clear cell ovarian carcinomas (CCOCs). The regulatory architecture of these cell types was compared with normal human mammary epithelial cells and LNCaP prostate cancer cells. We observed similar positional patterns of global enhancer signatures across the three different ovarian cancer precursor cell types, and evidence of tissue-specific regulatory signatures compared to non-gynecological cell types. We found significant enrichment for risk-associated SNPs intersecting regulatory biofeatures at 17 known HGSOC susceptibility loci in FTSECs (P = 3.8 × 10(-30)), OSECs (P = 2.4 × 10(-23)) and HMECs (P = 6.7 × 10(-15)) but not for EECs (P = 0.45) or LNCaP cells (P = 0.88). Hierarchical clustering of risk SNPs conditioned on the six different cell types indicates FTSECs and OSECs are highly related (96% of samples using multi-scale bootstrapping) suggesting both cell types may be precursors of HGSOC. These data represent the first description of regulatory catalogues of normal precursor cells for different ovarian cancer subtypes, and provide unique insights into the tissue specific regulatory variation with respect to the likely functional targets of germline genetic susceptibility variants for ovarian cancer.

Comprehensive functional annotation of 77 prostate cancer risk loci.

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations--we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at r(2) ≥ 0.88%. 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium (r(2) = 0.91) with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process.

A comprehensive examination of breast cancer risk loci in African American women.

Genome-wide association studies have identified 73 breast cancer risk variants mainly in European populations. Given considerable differences in linkage disequilibrium structure between populations of European and African ancestry, the known risk variants may not be informative for risk in African ancestry populations. In a previous fine-mapping investigation of 19 breast cancer loci, we were able to identify SNPs in four regions that better captured risk associations in African American women. In this study of breast cancer in African American women (3016 cases, 2745 controls), we tested an additional 54 novel breast cancer risk variants. Thirty-eight variants (70%) were found to have an association with breast cancer in the same direction as previously reported, with eight (15%) replicating at P < 0.05. Through fine-mapping, in three regions (1q32, 3p24, 10q25), we identified variants that better captured associations with overall breast cancer or estrogen receptor positive disease. We also observed suggestive associations with variants (at P < 5 × 10(-6)) in three separate regions (6q25, 14q13, 22q12) that may represent novel risk variants. Directional consistency of association observed for ∼65-70% of currently known genetic variants for breast cancer in women of African ancestry implies a shared functional common variant at most loci. To validate and enhance the spectrum of alleles that define associations at the known breast cancer risk loci, as well as genome-wide, will require even larger collaborative efforts in women of African ancestry.

Nucleosome positioning and histone modifications define relationships between regulatory elements and nearby gene expression in breast epithelial cells.

BACKGROUND: The precise nature of how cell type specific chromatin structures at enhancer sites affect gene expression is largely unknown. Here we identified cell type specific enhancers coupled with gene expression in two different types of breast epithelial cells, HMEC (normal breast epithelial cells) and MDAMB231 (triple negative breast cancer cell line). RESULTS: Enhancers were defined by modified neighboring histones [using chromatin immunoprecipitation followed by sequencing (ChIP-seq)] and nucleosome depletion [using formaldehyde-assisted isolation of regulatory elements followed by sequencing (FAIRE-seq)]. Histone modifications at enhancers were related to the expression levels of nearby genes up to 750 kb away. These expression levels were correlated with enhancer status (poised or active), defined by surrounding histone marks. Furthermore, about fifty percent of poised and active enhancers contained nucleosome-depleted regions. We also identified response element motifs enriched at these enhancer sites that revealed key transcription factors (e.g. TP63) likely involved in regulating breast epithelial enhancer-mediated gene expression. By utilizing expression data, potential target genes of more than 600 active enhancers were identified. These genes were involved in proteolysis, epidermis development, cell adhesion, mitosis, cell cycle, and DNA replication. CONCLUSIONS: These findings facilitate the understanding of epigenetic regulation specifically, such as the relationships between regulatory elements and gene expression and generally, how breast epithelial cellular phenotypes are determined by cell type specific enhancers.

A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11.

Genome-wide association studies (GWAS) of breast cancer defined by hormone receptor status have revealed loci contributing to susceptibility of estrogen receptor (ER)-negative subtypes. To identify additional genetic variants for ER-negative breast cancer, we conducted the largest meta-analysis of ER-negative disease to date, comprising 4754 ER-negative cases and 31 663 controls from three GWAS: NCI Breast and Prostate Cancer Cohort Consortium (BPC3) (2188 ER-negative cases; 25 519 controls of European ancestry), Triple Negative Breast Cancer Consortium (TNBCC) (1562 triple negative cases; 3399 controls of European ancestry) and African American Breast Cancer Consortium (AABC) (1004 ER-negative cases; 2745 controls). We performed in silico replication of 86 SNPs at P ≤ 1 × 10(-5) in an additional 11 209 breast cancer cases (946 with ER-negative disease) and 16 057 controls of Japanese, Latino and European ancestry. We identified two novel loci for breast cancer at 20q11 and 6q14. SNP rs2284378 at 20q11 was associated with ER-negative breast cancer (combined two-stage OR = 1.16; P = 1.1 × 10(-8)) but showed a weaker association with overall breast cancer (OR = 1.08, P = 1.3 × 10(-6)) based on 17 869 cases and 43 745 controls and no association with ER-positive disease (OR = 1.01, P = 0.67) based on 9965 cases and 22 902 controls. Similarly, rs17530068 at 6q14 was associated with breast cancer (OR = 1.12; P = 1.1 × 10(-9)), and with both ER-positive (OR = 1.09; P = 1.5 × 10(-5)) and ER-negative (OR = 1.16, P = 2.5 × 10(-7)) disease. We also confirmed three known loci associated with ER-negative (19p13) and both ER-negative and ER-positive breast cancer (6q25 and 12p11). Our results highlight the value of large-scale collaborative studies to identify novel breast cancer risk loci.

FunciSNP: an R/bioconductor tool integrating functional non-coding data sets with genetic association studies to identify candidate regulatory SNPs.

UNLABELLED: Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance. AVAILABILITY: FunciSNP is available from Bioconductor (bioconductor.org).

A genome-wide association study of breast cancer in women of African ancestry.

Genome-wide association studies (GWAS) in diverse populations are needed to reveal variants that are more common and/or limited to defined populations. We conducted a GWAS of breast cancer in women of African ancestry, with genotyping of >1,000,000 SNPs in 3,153 African American cases and 2,831 controls, and replication testing of the top 66 associations in an additional 3,607 breast cancer cases and 11,330 controls of African ancestry. Two of the 66 SNPs replicated (p < 0.05) in stage 2, which reached statistical significance levels of 10(-6) and 10(-5) in the stage 1 and 2 combined analysis (rs4322600 at chromosome 14q31: OR = 1.18, p = 4.3 × 10(-6); rs10510333 at chromosome 3p26: OR = 1.15, p = 1.5 × 10(-5)). These suggestive risk loci have not been identified in previous GWAS in other populations and will need to be examined in additional samples. Identification of novel risk variants for breast cancer in women of African ancestry will demand testing of a substantially larger set of markers from stage 1 in a larger replication sample.

Genomic, epigenomic, and transcriptomic signatures of prostate cancer between African American and European American patients.

Prostate cancer is the second most common cancer in men in the United States, and racial disparities are greatly observed in the disease. Specifically, African American (AA) patients have 60% higher incidence and mortality rates, in addition to higher grade and stage prostate tumors, than European American (EA) patients. In order to narrow the gap between clinical outcomes for these two populations, genetic and molecular signatures contributing to this disparity have been characterized. Over the past decade, profiles of prostate tumor samples from different ethnic groups have been developed using molecular and functional assays coupled with next generation sequencing or microarrays. Comparative genome-wide analyses of genomic, epigenomic, and transcriptomic profiles from prostate tumor samples have uncovered potential race-specific mutations, copy number alterations, DNA methylation, and gene expression patterns. In this study, we reviewed over 20 published studies that examined the aforementioned molecular contributions to racial disparities in AA and EA prostate cancer patients. The reviewed genomic studies revealed mutations, deletions, amplifications, duplications, or fusion genes differentially enriched in AA patients relative to EA patients. Commonly reported genomic alterations included mutations or copy number alterations of FOXA1, KMT2D, SPOP, MYC, PTEN, TP53, ZFHX3, and the TMPRSS2-ERG fusion. The reviewed epigenomic studies identified that CpG sites near the promoters of PMEPA1, RARB, SNRPN, and TIMP3 genes were differentially methylated between AA and EA patients. Lastly, the reviewed transcriptomic studies identified genes (e.g. CCL4, CHRM3, CRYBB2, CXCR4, GALR1, GSTM3, SPINK1) and signaling pathways dysregulated between AA and EA patients. The most frequently found dysregulated pathways were involved in immune and inflammatory responses and neuroactive ligand signaling. Overall, we observed that the genomic, epigenomic, and transcriptomic alterations evaluated between AA and EA prostate cancer patients varied between studies, highlighting the impact of using different methods and sample sizes. The reported genomic, epigenomic, and transcriptomic alterations do not only uncover molecular mechanisms of tumorigenesis but also provide researchers and clinicians valuable resources to identify novel biomarkers and treatment modalities to improve the disparity of clinical outcomes between AA and EA patients.

m6A epitranscriptome analysis reveals differentially methylated transcripts that drive early chemoresistance in bladder cancer.

N 6-Methyladenosine (m6A) RNA modifications dynamically regulate messenger RNA processing, differentiation and cell fate. Given these functions, we hypothesized that m6A modifications play a role in the transition to chemoresistance. To test this, we took an agnostic discovery approach anchored directly to chemoresistance rather than to any particular m6A effector protein. Specifically, we used methyl-RNA immunoprecipitation followed by sequencing (MeRIP-seq) in parallel with RNA sequencing to identify gene transcripts that were both differentially methylated and differentially expressed between cisplatin-sensitive and cisplatin-resistant bladder cancer (BC) cells. We filtered and prioritized these genes using clinical and functional database tools, and then validated several of the top candidates via targeted quantitative polymerase chain reaction (qPCR) and MeRIP-PCR. In cisplatin-resistant cells, SLC7A11 transcripts had decreased methylation associated with decreased m6A reader YTHDF3 binding, prolonged RNA stability, and increased RNA and protein levels, leading to reduced ferroptosis and increased survival. Consistent with this, cisplatin-sensitive BC cell lines and patient-derived organoids exposed to cisplatin for as little as 48 h exhibited similar mechanisms of SLC7A11 upregulation and chemoresistance, trends that were also reflected in public cancer survival databases. Collectively, these findings highlight epitranscriptomic plasticity as a mechanism of rapid chemoresistance and a potential therapeutic target.

Validation of Myc-Associated Protein X (MAX) regulation in growth hormone secreting and nonfunctional pituitary adenoma.

INTRODUCTION: Many patients with growth hormone-secreting pituitary adenoma (GHPA) fail to achieve biochemical remission, warranting investigation into epigenetic and molecular signatures associated with tumorigenesis and hormonal secretion. Prior work exploring the DNA methylome showed Myc-Associated Protein X (MAX), a transcription factor involved in cell cycle regulation, was differentially methylated between GHPA and nonfunctional pituitary adenoma (NFPA). We aimed to validate the differential DNA methylation and related MAX protein expression profiles between NFPA and GHPA. METHODS: DNA methylation levels were measured in 52 surgically resected tumors (37 NFPA, 15 GHPA) at ~100,000 known MAX binding sites derived using ChIP-seq analysis from ENCODE. Findings were correlated with MAX protein expression using a constructed tissue microarray (TMA). Gene ontology analysis was performed to explore downstream genetic and signaling pathways regulated by MAX. RESULTS: GHPA had more hypomethylation events across all known MAX binding sites. Of binding sites defined using ChIP-seq analysis, 1,551 sites had significantly different methylation patterns between the two cohorts; 432 occurred near promoter regions potentially regulated by MAX, including promoters of TNF and MMP9. Gene ontology analysis suggested enrichment in genes involved in oxygen response, immune system regulation, and cell proliferation. Thirteen MAX binding sites were within coding regions of genes. GHPA demonstrated significantly increased expression of MAX protein compared to NFPA. CONCLUSION: GHPA have significantly different DNA methylation and downstream protein expression levels of MAX compared to NFPA. These differences may influence mechanisms involved with cellular proliferation, tumor invasion and hormonal secretion.

Phosphorylation and stabilization of EZH2 by DCAF1/VprBP trigger aberrant gene silencing in colon cancer.

Our recent work has shown that DCAF1 (also known as VprBP) is overexpressed in colon cancer and phosphorylates histone H2AT120 to drive epigenetic gene inactivation and oncogenic transformation. We have extended these observations by investigating whether DCAF1 also phosphorylates non-histone proteins as an additional mechanism linking its kinase activity to colon cancer development. We now demonstrate that DCAF1 phosphorylates EZH2 at T367 to augment its nuclear stabilization and enzymatic activity in colon cancer cells. Consistent with this mechanistic role, DCAF1-mediated EZH2 phosphorylation leads to elevated levels of H3K27me3 and altered expression of growth regulatory genes in cancer cells. Furthermore, our preclinical studies using organoid and xenograft models revealed that EZH2 requires phosphorylation for its oncogenic function, which may have therapeutic implications for gene reactivation in colon cancer cells. Together, our data define a mechanism underlying DCAF1-driven colonic tumorigenesis by linking DCAF1-mediated EZH2 phosphorylation to EZH2 stability that is crucial for establishing H3K27me3 and gene silencing program.

Cell-free DNA Methylation as a Predictive Biomarker of Response to Neoadjuvant Chemotherapy for Patients with Muscle-invasive Bladder Cancer in SWOG S1314.

BACKGROUND: Neoadjuvant chemotherapy (NAC) is the standard of care in muscle-invasive bladder cancer (MIBC). However, treatment is intense, and the overall benefit is small, necessitating effective biomarkers to identify patients who will benefit most. OBJECTIVE: To characterize cell-free DNA (cfDNA) methylation in patients receiving NAC in SWOG S1314, a prospective cooperative group trial, and to correlate the methylation signatures with pathologic response at radical cystectomy. DESIGN, SETTING, AND PARTICIPANTS: SWOG S1314 is a prospective cooperative group trial for patients with MIBC (cT2-T4aN0M0, ≥5 mm of viable tumor), with a primary objective of evaluating the coexpression extrapolation (COXEN) gene expression signature as a predictor of NAC response, defined as achieving pT0N0 or ≤pT1N0 at radical cystectomy. For the current exploratory analysis, blood samples were collected prospectively from 72 patients in S1314 before and during NAC, and plasma cfDNA methylation was measured using the Infinium MethylationEPIC BeadChip array. INTERVENTION: No additional interventions besides plasma collection. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Differential methylation between pathologic responders (≤pT1N0) and nonresponders was analyzed, and a classifier predictive of treatment response was generated using the Random Forest machine learning algorithm. RESULTS AND LIMITATIONS: Using prechemotherapy plasma cfDNA, we developed a methylation-based response score (mR-score) predictive of pathologic response. Plasma samples collected after the first cycle of NAC yielded mR-scores with similar predictive ability. Furthermore, we used cfDNA methylation data to calculate the circulating bladder DNA fraction, which had a modest but independent predictive ability for treatment response. In a model combining mR-score and circulating bladder DNA fraction, we correctly predicted pathologic response in 79% of patients based on their plasma collected at baseline and after one cycle of chemotherapy. Limitations of this study included a limited sample size and relatively low circulating bladder DNA levels. CONCLUSIONS: Our study provides the proof of concept that cfDNA methylation can be used to generate classifiers of NAC response in bladder cancer patients. PATIENT SUMMARY: In this exploratory analysis of S1314, we demonstrated that cell-free DNA methylation can be profiled to generate biomarker signatures associated with neoadjuvant chemotherapy response. With validation in additional cohorts, this minimally invasive approach may be used to predict chemotherapy response in locally advanced bladder cancer and perhaps also in metastatic disease.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry.

BACKGROUND: Genetic factors play an important role in prostate cancer (PCa) susceptibility. OBJECTIVE: To discover common genetic variants contributing to the risk of PCa in men of African ancestry. DESIGN, SETTING, AND PARTICIPANTS: We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness. RESULTS AND LIMITATIONS: Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4). CONCLUSIONS: This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry. PATIENT SUMMARY: In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.