Chen, Zhihua

Abstract Purpose: Identifying dysregulated gene expression in ovarian cancers has been limited by a deficit of available normal tissues. Here, we generated the largest set of high-grade serous ovarian cancer (HGSOC) tumors with normal precursor tissues for transcriptome analyses. Methods: We performed RNA sequencing on 220 primary HGSOCs and 116 benign epithelia (micro-dissected fallopian tube, ovarian surface, and inclusion cysts), and combined samples with 428 HGSOCs from TCGA, 60 HGSOCs from a prior study, and 180 bulk ovary tissues from GTEx. Raw reads were processed with a uniform bioinformatic and quality control pipeline; combined data were batch corrected. We tested for differences in median normalized CPM expression values using the Wilcoxon rank sum test with >2-fold change and a false discovery rate <1% considered statistically significant. We also conducted weighted gene co-expression network analysis in each tissue type. The hypergeometric test was used for enrichment of differentially expressed genes (DEGs) and gene ontologies within modules. Results: Transcriptomes comprised 27,700 expressed genes (8,202 lncRNAs) across 706 HGSOCs, 180 bulk ovary, and 88 ovarian epithelia tissues. Most (~90%) genes were expressed in all tissues; 5% each showed HGSOC- or normal tissue-specific expression and ≥50% were lncRNA. Comparing HGSOCs to ovarian epithelia and to bulk ovary identified 11,804 DEGs with 4,522 lncRNAs (DElncRNA) of which ~50% were tissue-specific. DEGs included MUC16 and multiple GWAS/TWAS implicated susceptibility genes including RAD51, BRIP1, BNC2, TIPARP-AS1, PRC1, KANSL1, ANKLE1, CHMP4C, ESRP2, and CCNE1. The most highly expressed DElncRNA in HGSOC were upregulated RMRP (P=1.4x10-39), SNHG1 (P=3.0x10-27), and HAGLR (P=2.0x10-24) at the HOXD risk locus. The highest expressed DElncRNA in precursor tissues was the HGSOC-downregulated MEG3 (P=1.7x10-67). DEGs were enriched in HGSOC co-expression modules associated with immune response, cell motility/localization, cell cycle regulation, angiogenesis, reproductive development, transcriptional regulation, and metabolic processes. Tissue-specific DElncRNA tended toward upregulation compared to ovarian epithelia with enriched modules associated with cell cycle regulation (hub=BUB1B; top DElncRNA TRPM2-AS, P=1.5x10-38); and toward downregulation compared to bulk ovary with enriched modules associated with chemokine signaling/response (hub= GADD45B, top DElncRNA RP11-87P3.1, P=1.7x10-113). Conclusion: HGSOC-dysregulated lncRNA expression revealed tissue-specific differences that highlight unique biological pathways in precursor epithelia and the ovarian microenvironment that contribute to HGSOC pathogenesis. Our results provide additional evidence to support previously nominated risk genes. Integration with eQTL and GWAS are underway to further elucidate novel HGSOC susceptibility genes. Citation Format: Brett M. Reid, Ann Chen, Zhihua Chen, Florian A. Karreth, Peter Kanetsky, Jennifer B. Permuth, Ozlen Saglam, Jamie Teer, Xiaoqing Yu, Simon Gayther, Ellen Goode, Paul Pharoah, Thomas A. Sellers, Kate Lawrenson. Patterns of dysregulated coding and noncoding gene expression in high-grade serous ovarian carcinomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6065.

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni < 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P < 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

Abstract Background: Genome-wide association studies (GWAS) have identified multiple loci associated with epithelial ovarian cancer (EOC) susceptibility, but further progress requires integration of epidemiology and biology to illuminate true risk loci below genome-wide significance levels (P < 5 × 10−8). Most risk SNPs lie within non–protein-encoding regions, and we hypothesize that long noncoding RNA (lncRNA) genes are enriched at EOC risk regions and represent biologically relevant functional targets. Methods: Using imputed GWAS data from about 18,000 invasive EOC cases and 34,000 controls of European ancestry, the GENCODE (v19) lncRNA database was used to annotate SNPs from 13,442 lncRNAs for permutation-based enrichment analysis. Tumor expression quantitative trait locus (eQTL) analysis was performed for sub-genome-wide regions (1 × 10−5 > P > 5 × 10−8) overlapping lncRNAs. Results: Of 5,294 EOC-associated SNPs (P < 1.0 × 10−5), 1,464 (28%) mapped within 53 unique lncRNAs and an additional 3,484 (66%) SNPs were correlated (r2 > 0.2) with SNPs within 115 lncRNAs. EOC-associated SNPs comprised 130 independent regions, of which 72 (55%) overlapped with lncRNAs, representing a significant enrichment (P = 5.0 × 10−4) that was more pronounced among a subset of 5,401 lncRNAs with active epigenetic regulation in normal ovarian tissue. EOC-associated lncRNAs and their putative promoters and transcription factors were enriched for biologically relevant pathways and eQTL analysis identified five novel putative risk regions with allele-specific effects on lncRNA gene expression. Conclusions: lncRNAs are significantly enriched at EOC risk regions, suggesting a mechanistic role for lncRNAs in driving predisposition to EOC. Impact: lncRNAs represent key candidates for integrative epidemiologic and functional studies. Further research on their biologic role in ovarian cancer is indicated. Cancer Epidemiol Biomarkers Prev; 26(1); 116–25. ©2016 AACR.

Background The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study. Methods We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant. Results For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10−5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10−8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants. Conclusions This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

ABSTRACTEpithelial‐mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single‐nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome‐wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT‐related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive‐cancer patients and 23,447 controls. A P‐value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07–1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69–0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.

Abstract Background: Recent bioinformatic analyses and high-throughput sequencing efforts have revealed that RNA editing is a widespread post-transcriptional modification that affects non-coding repetitive elements throughout the genome. This modification, which involves the conversion of adenosine into inosine, is mediated by the action of the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on reports of dysfunctional ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in five ADAR genes (ADAD1, ADAR1, ADAR2, ADAR3, SND1) may contribute to EOC susceptibility, potentially by altering the expression of corresponding genes in ovarian tissues. Methods: Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 healthy controls of European ancestry represented in the international Collaborative Oncological Gene-environment Study, we evaluated the frequency of 7,133 SNPs in 5 RNA editing genes. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI) between genotypes and case status under a log-additive genetic model, with adjustment for the first five principal components representing European ancestry. Subgroup analysis was conducted for the four main histologic sub-types of EOC. We also aggregated SNP-level data and conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for the combined effect of common and rare variants (SKAT-CR). Results: Association analysis revealed a top risk-associated SNP (OR (95% CI) = 1.39 (1.17-1.64)) in ADAR3 with P = 1.3 x10−4, 17 additional SNPs in ADAR1, ADAR3, or SND1 with P = 10−4, and 51 SNPs with P = 10−3; this observed number of statistically-significant SNPs was greater than expected. When restricting to the 6,500 invasive serous adenocarcinomas, the magnitude of association strengthened for numerous SNPs. Exploratory analysis for the less common histologic subtypes (endometrioid (n = 1,439), mucinous (n = 696), and clear cell (n = 660)) revealed SNP-level associations (P&lt;10−3) unique to each sub-type. Gene-level analyses based on AML and SKAT-CR revealed that ADAR1 was globally associated (P&lt;0.10) with susceptibility to invasive EOC, with P = 0.024 and P = 0.023, respectively. Expression quantitative trait locus analysis is underway to evaluate associations between genotype and gene expression in EOC and precursor tissues using data from The Cancer Genome Atlas Project and other sources. Conclusions: Preliminary data from this large-scale collaboration suggest that germline variation involving RNA editing genes may influence susceptibility to EOC. This novel finding warrants further investigation, as inherited and acquired alterations affecting RNA editing may serve as important biological mechanisms that promote the development of EOC and possibly other malignancies. Citation Format: Jennifer Permuth-Wey, Brett M. Reid, Y. Ann Chen, Hui-Yi Lin, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Simon Gayther, Ellen Goode, Edwin Iversen, Leigh Pearce, Paul D. P. Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Thomas Sellers, on behalf of the Ovarian Cancer Association Consortium. Inherited variants affecting RNA editing may contribute to ovarian cancer susceptibility. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4635. doi:10.1158/1538-7445.AM2015-4635

Abstract Background: Cell-type-specific super enhancers, or large clusters of mediator-rich enhancers, are highly transcribed regions of the genome with key functions in the maintenance of cell identity and control. These regulatory elements are easily perturbed and enriched for disease-associated sequence variation. Identification of disease-associated super-enhancers and investigation of their sequence variants may help identify important, biologically relevant associations previously overlooked by genome-wide studies. Thus, we sought to identify variants within ovarian tissue specific super-enhancers and investigate their association with epithelial ovarian cancer (EOC) susceptibility. Methods: We utilized genotype data from a European population of ∼11,000 invasive EOC cases and 22,000 healthy controls from the Ovarian Cancer Association Consortium. A public catalog of 478 super-enhancers, identified by peak histone H3K27ac ChIP-seq signal in ovarian tissue, was used to identify 72,116 variants located within 344 super-enhancer regions. Unconditional logistic regression (log-additive model) was used to assess individual SNP-EOC risk associations. Additionally, we reviewed associations within 500 kb of super-enhancer regions to identify proximal associated variants. When LD structure indicated, we employed conditional analysis to determine independence of the super-enhancer signal. Gene-EOC risk associations were assessed with the Admixture Maximum Likelihood (AML) test and Sequence-Kernel Association Test (SKAT) for the combined effect of common and rare (MAF&lt;0.01) variants. Models were adjusted for population stratification using the first five principal components. To adjust for multiple comparisons, a false discovery rate of 10% was used. Results: The most significant SNP observed was rs147487657 (p = 7.7×10−6) on locus 2q37.3 in the HDAC4 gene and there were an additional 12 SNPs with p&lt;1.0×10−4, although not significant at 10% FDR. Gene-level analyses identified two super-enhancers significantly associated with EOC risk. An enhancer at the previously reported 2q31 locus, containing HOXD8 and HOXD-AS2, associated (SKAT p = 4.3×10−5; AML p&lt;0.001) with the top SNP located within the 5′UTR of HOXD8 transcription factor (rs847163 (T&gt;C); p&lt;2.08×10−5). Conditional analysis did not identify an independent association within the super-enhancer region. The second EOC-associated super-enhancer identified was at locus 1q32 (SKAT p = 0.00052, AML p = 0.001) containing the CSRP1 gene. There were 6 SNPs with p&lt;0.001 and 16 SNPs with p&lt;0.01, all were inversely associated with risk and consistent with the top hit at rs927901 T&gt;C (OR = 0.94; 95% CI: 0.90-0.97; p = 2.1×10−4). Associations varied by histologic subtype, possibly due to variance in disease relevant super-enhancer cell-type.Conclusion: These findings suggest germline variation within super-enhancer regions may contribute to EOC susceptibility Citation Format: Brett M. Reid, Jennifer Permuth-Wey, Y. Ann Chen, Hui-Yi Lin, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Simon Gayther, Ellen L. Goode, Edwin Iversen, Leigh Pearce, Paul Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Thomas Sellers, on behalf of the Ovarian Cancer Association Consortium. Variants within super-enhancer regulatory elements associate with epithelial ovarian cancer risk. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4634. doi:10.1158/1538-7445.AM2015-4634

Abstract Background: Genome wide association studies (GWAS) have identified 20 loci associated with epithelial ovarian cancer (EOC) risk. Additional EOC risk loci await identification, and biological approaches may represent a useful strategy for overcoming sample size limitations. The ENCODE project recently concluded that only ∼1.2% of the genome encodes proteins, but at least 20% exhibits biological function and over 80% exhibits biochemical indices of function. Post-GWAS follow-up studies have firmly established that some associations are due to inter-individual differences in non-coding RNAs (ncRNAs). To estimate the magnitude of impact of variants in long non-coding RNAs (lncRNAs), we assessed enrichment of EOC-associated SNPs in lncRNA regions by comparing the density of EOC-associated loci between lncRNA regions, non-lncRNA regions, and the whole genome. Methods: The study population included ∼18,000 invasive EOC cases and 34,000 healthy controls of European ancestry from the Ovarian Cancer Association Consortium with GWAS data imputed to 1000 Genomes Project (1KGP) density. Density was compared with two measures: average chromosome length required for one EOC-associated SNP (kb/locus) and average number of tested SNPs containing one EOC-associated SNP (SNPs/locus). Coordinates and annotation for 13,870 lncRNA genes were downloaded from the publicly available GENCODE (v19) database and were used to annotate SNPs from the 1KGP imputed GWAS data, yielding 13,192 lncRNA genes with biallelic variants imputed at rsquare ≥0.25. A SNP was considered EOC-associated if it reached a significance level of P&lt;10−5. Results: Genome-wide, 11.6% of the 15,123,646 tested SNPs fell within lncRNA regions. Of 5,294 EOC-associated SNPs that were identified, 27.7% were located within lncRNAs, a higher frequency than expected (&lt;v:shape style = “WIDTH: 16.5pt; HEIGHT: 18pt” id = _x0000_i1026 equationxml = ‘14pχ2′ type = “#_x0000_t75”&gt; = 9.35×10−22). Additionally, EOC-associated SNPs were located on average every 237.6 kb (or 1204 tested SNPs) in lncRNA regions versus every 657.7 kb (or 3489 SNPS) in non-lncRNA regions and every 541.5 kb (or 2857 SNPs) in the whole genome, clearly demonstrating an enrichment of EOC-associated loci in lncRNA regions. A total of 1,464 lncRNA SNPs from 53 unique lncRNAs associated with EOC risk (P&lt;10−5); 70 of the identified lncRNA SNPs had P-values between 10−20 and 10−32, and 221 had P-values between 10−10 and 10−15, representing 21 unique regions &gt; 500kb apart. Conclusions: These results justify efforts to integrate functional and population level genetic data on lncRNA regions to explore a new paradigm that addresses the problem of the missing heritability for cancer where restrictions in sample size and statistical power limit the ability to discover additional risk loci. Citation Format: Thomas A. Sellers, Brett M. Reid, Y. Ann Chen, Hui-Yi Lin, Edward Richards, Jamie Teer, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Ellen Goode, Edwin Iverson, Leigh Pearce, Paul Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Jin Cheng, Simon Gayther, Jennifer Permuth-Wey, on behalf of the Ovarian Cancer Association Consortium. Evidence that long non-coding RNA variants associate with epithelial ovarian cancer risk. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4633. doi:10.1158/1538-7445.AM2015-4633

Abstract Background: Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations. Methods: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls). Results: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P &lt; 0.05 and FDR &lt; 0.05). These results were replicated (P &lt; 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network. Conclusion: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development. Impact: Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization. Cancer Epidemiol Biomarkers Prev; 24(10); 1574–84. ©2015 AACR.