Chen, Zhihua

Abstract Introduction: 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 to identify lncRNA expression patterns, provide insight into their cellular contexts and functions, and assess their relevance to HGSOC pathogenesis and genetic susceptibility. Methods: We performed RNA sequencing on 220 primary HGSOCs and 116 benign epithelia (micro-dissected fallopian tube, ovarian surface, and inclusion cyst epithelia), and combined samples with 428 HGSOCs from TCGA, 60 HGSOCs from a prior study, and 180 bulk ovary tissues from GTEx. Reads were processed with a uniform bioinformatic and quality control pipeline, combined and batch corrected. Normalized TPM expression values were compared using the Wilcoxon rank sum test with >2-fold change and FDR <1% considered statistically significant. We computed weighted gene co-expression networks for each tissue type and tested for enrichment of differentially expressed genes (DEGs) and gene ontologies within co-expression modules. Results: Transcriptomes included 27,700 expressed genes (8,202 lncRNAs) across 706 HGSOCs, 180 bulk ovary, and 88 ovarian epithelia tissues. Most genes were expressed in all tissues, though 5% each showed HGSOC- or normal tissue-specific expression which comprised ≥50% lncRNA. Comparing HGSOCs to normal tissues 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.4 × 10−39), SNHG1 (P=3.0 × 10−27), and HAGLR (P=2.0 × 10−24) at the HOXD risk locus. Lower proportions of lncRNA were seen in the HGSOC co-expression modules compared to those in normal tissues and less (3% vs. 15%) HGSOC modules were considered strongly lncRNA regulated (>50% lncRNA in module). Tissue-specific DElncRNA tended toward upregulation in HGSOC compared to ovarian epithelia with enriched modules associated with cell cycle regulation (P=1.5 × 10−38); and toward downregulation in HGSOC compared to bulk ovary with enriched modules associated with chemokine signaling/response (P=1.7 × 10−113). Conclusions: The lower proportion of lncRNA involved in HGSOC co-expression networks indicates a marked loss of regulation by lncRNA. Preliminary interrogation of 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 support for previously nominated risk genes. eQTL analysis has recently been completed and integration with GWAS are underway to elucidate lncRNA and lncRNA-regulated susceptibility genes and will be presented at the conference. Citation Format: Brett M. Reid, Ann Chen, Zhihua Chen, Florian A. Karreth, Peter Kanetsky, Jennifer B. Permuth, Ozlen Saglam, Jamie K. Teer, Xiaoqing Yu, Simon Gayther, Ellen Goode, Paul Pharoah, Thomas A. Sellers, Kate Lawrenson, Jonathan Tyrer. Patterns of lncRNA-regulated gene expression in high-grade serous ovarian carcinomas [abstract]. In: Proceedings of the AACR Special Conference on Ovarian Cancer; 2023 Oct 5-7; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_2):Abstract nr A004.

AbstractBackground:Germline DNA copy number variation (CNV) is a ubiquitous source of genetic variation and remains largely unexplored in association with epithelial ovarian cancer (EOC) risk.Methods:CNV was quantified in the DNA of approximately 3,500 cases and controls genotyped with the Illumina 610k and HumanOmni2.5M arrays. We performed a genome-wide association study of common (>1%) CNV regions (CNVRs) with EOC and high-grade serous (HGSOC) risk and, using The Cancer Genome Atlas (TCGA), performed in silico analyses of tumor-gene expression.Results:Three CNVRs were associated (P < 0.01) with EOC risk: two large (∼100 kb) regions within the 610k set and one small (<5 kb) region with the higher resolution 2.5M data. Large CNVRs included a duplication at LILRA6 (OR = 2.57; P = 0.001) and a deletion at CYP2A7 (OR = 1.90; P = 0.007) that were strongly associated with HGSOC risk (OR = 3.02; P = 8.98 × 10−5). Somatic CYP2A7 alterations correlated with EGLN2 expression in tumors (P = 2.94 × 10−47). An intronic ERBB4/HER4 deletion was associated with reduced EOC risk (OR = 0.33; P = 9.5 × 10−2), and somatic deletions correlated with ERBB4 downregulation (P = 7.05 × 10−5). Five CNVRs were associated with HGSOC, including two reduced-risk deletions: one at 1p36.33 (OR = 0.28; P = 0.001) that correlated with lower CDKIIA expression in TCGA tumors (P = 2.7 × 10−7), and another at 8p21.2 (OR = 0.52; P = 0.002) that was present somatically where it correlated with lower GNRH1 expression (P = 5.9 × 10−5).Conclusions:Though CNV appears to not contribute largely to EOC susceptibility, a number of low-to-common frequency variants may influence the risk of EOC and tumor-gene expression.Impact:Further research on CNV and EOC susceptibility is warranted, particularly with CNVs estimated from high-density arrays.

Abstract IL-17 is known to play important roles in immune and inflammatory disease, such as in asthma, but its functions in allergic airway inflammation are still controversial, and the molecular mechanisms mediating these functions remain unclear. Increased production of eosinophils in bone marrow and their emergence in the airway have been linked to the onset and progression of allergic asthma. In this study, we investigated the effects of exogenous IL-17 on allergic airway inflammation and explored the underlying molecular mechanisms through eosinophil generation. Exogenous IL-17 significantly attenuated the features of allergic inflammation induced by ovalbumin in mice. It inhibited eosinophil differentiation both in vivo and in vitro, accompanied by down-regulated expression of CC chemokine receptor 3, GATA binding protein 1 (GATA-1), and GATA binding protein 2 (GATA-2), as well as reduced formation of common myeloid progenitors and eosinophil progenitors, but without influencing eosinophil apoptosis. IL-17 also significantly decreased the number of eosinophils in IL-5–transgenic mice, although it notably increased the levels of IL-3, IL-5, and granulocyte/macrophage colony–stimulating factor. In addition, IL-17 had little effect on secretion of the inflammatory cytokines by eosinophils. Neutralization of endogenous IL-17 significantly augmented eosinophil recruitment in the airways. Together, these findings suggest that exogenous IL-17 protects against allergic airway inflammation, most likely through inhibition of the eosinophil differentiation in bone marrow.

Abstract Background: Circadian rhythms of biological processes are regulated by endogenous clock genes and clock-controlled genes. Aberrant expression of circadian clock genes may have important consequences on the transactivation of downstream targets that control the cell cycle and cellular proliferation potentially promoting carcinogenesis. Animal models indicate that several circadian rhythm genes are expressed in the ovaries, where they influence and are influenced by estrous cycles. The goal of the current study was to examine the association of circadian gene variants and epithelial ovarian cancer (EOC) risk. Methods: Thirty-one SNPs from five circadian genes (i.e., ARNTL, CRY2, KLF10, NPAS2, PER3, TIMELESS) were genotyped in 14,736 cases and 23,448 control women of European ancestry from 43 studies in the Ovarian Cancer Association Consortium (OCAC), on a custom Illumina iSelect designed for the Collaborative Oncological Gene-Environment Study (COGS). Both invasive cancers combined and the four main histological subtypes (serous [n=8,372], endometroid [n=2,068], clear cell [n=1,025] and mucinous [n=943] were analyzed, SNP analyses were conducted using unconditional logistic regression under a log-additive model. All analyses were adjusted for study site and population substructure. Results: Eleven SNPs were found to be associated with ovarian cancer risk. The SNPs most associated with serous cancer risk were KLF10 rs2513928 (OR=.95 p=6.1x10−3), rs2513927 (OR=1.05 p=6.1x10−3), rs2511703 (OR=1.05 p=9.8x10−3), rs3191333 (OR=1.05, p=1.5x10−2) and NPAS2 rs13012930 (OR=.95, p=3.5x10−2). Endometroid cancer risk was significantly associated with ARNTL SNPs rs10732458 (OR=1.3, p=1.0x10−2) and rs7117836 (OR=1.2, p=4.6x10−2), while clear cell cancer risk was associated with ARNTL rs1562438 (OR=.88, p=1.5x10−2), rs1026071 (OR=.89, p=1.8x10−2), and rs3816360 (OR=.91, p=4.8x10−2) and KLF10 rs2388232 (OR=1.1, p=3.3x10−2). No variants were significantly associated with mucinous cancer risk. Four SNPs in KLF10 were associated with cancer invasiveness; they were rs2513928 (OR=.95, p=1.8x10−3), rs3191333 (OR=1.04, p=1.4x10−2), rs2513927 (OR=1.04, p=1.9x10−2), and rs2511703 (OR=1.03, p=2.8x10−2). Conclusions: Data from the current study suggest that polymorphisms in circadian genes ARNTL, KLF10, and NPAS2 are significantly associated with ovarian cancer histopathologic subtypes and invasiveness. These findings merit further investigation and replication. Funding: R01 CA149429 Citation Format: Heather Jim, Jonathan Tyrer, Hui-Yi Lin, Gang Han, Xiaotao Qu, Ellen L. Goode, Zhihua Chen, Ya-Yu Tsai, Julie M. Cunningham, Edward Iversen, Susan Ramus, Andrew Berchuck, Joellen Schildkraut, Alvaro Monteiro, Simon Gayther, Steven A. Narod, Thomas A. Sellers, Paul Pharoah, Catherine M. Phelan. Variation in circadian rhythm genes influence epithelial ovarian cancer risk and invasiveness. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 4850. doi:10.1158/1538-7445.AM2013-4850 Note: This abstract was not presented at the AACR Annual Meeting 2013 because the presenter was unable to attend.

Abstract Recent GWAS have identified and confirmed associations between several SNPs in the gene desert located on chromosome 8q24 with risk of cancers of the breast, prostate, glioma, bladder and colorectum. Evidence supporting an association with ovarian cancer is limited and inconsistent. Here we evaluate preliminary results in the 8q24 region from stage I of an ongoing invasive epithelial ovarian cancer GWAS among 1864 cases (1096 serous cases) and 1912 controls. All subjects were self-reported non-Hispanic non-Jewish Caucasian, and were genotyped with the Illumina 610quad Chip. Samples and SNPs with call rates less than 95% were excluded (9.6% and 9.9%, respectively) as well as 128 subjects with ambiguous gender, unresolved identical genotypes and less than 80% European ancestry were also excluded. Our most significant hit in the region was rs6651252 located at 8q24.21 with p = 0.000125 for all ovarian cancer, and 1.02*10−5 for serous cases only. This region is located approximately 849 kb downstream of the previous 8q24.21 findings on other cancers, and was not in linkage disequilibrium with any of the haplotype blocks reported previously. We further examined the block structure of this region with eight genotyped SNPs located within 129.61 to 129.64 MB (rs10088218, rs1516973, rs1516975, rs16903097, rs10098821, rs6651252, rs1561928, and rs1561927) and performed haplotype analysis to evaluate the association with ovarian cancer risk. One haplotype (AAGCAGGG, estimated frequency=0.11) was significantly associated with decreased ovarian cancer risk with an odds ratio of 0.67 (all cases, 95% CI: 0.78-0.91, p=0.0026, empirical p= 0.0020). This region contains no known genes, but it is located upstream and overlaps with a known copy number variation (CNV) reported by Redon et al. (2006). Previous published literature on array CGH and CNV results from TCGA (The Cancer Genome Atlas) also indicate that the 8q24.21 region is associated with frequent gains/amplification in ovarian tumors. Combined analysis of SNPs and CNV identified in this region is underway to elucidate the joint effect of these two types of genetic variation on ovarian cancer risk. Our result indicates a novel locus on 8q24.21 is associated with epithelial ovarian cancer susceptibility. Haplotype analysis reveals one specific haplotype that is associated with decreased ovarian cancer risk compared to the most common haplotype. Given previous 8q24 findings from other cancers and evidence from publicly available resources, fine-mapping and resequencing of this region is warranted. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 4724.