Transcriptomic analysis of stage 1 versus advanced adult granulosa cell tumors

Comparative Tumor Microenvironment Analysis of Primary and Recurrent Ovarian Granulosa Cell Tumors

Adult-type granulosa cell tumors (aGCT) are rare ovarian sex cord tumors with few effective treatments for recurrent disease. The objective of this study was to characterize the tumor microenvironment (TME) of primary and recurrent aGCTs and to identify correlates of disease recurrence. Total RNA sequencing (RNA-seq) was performed on 24 pathologically confirmed, cryopreserved aGCT samples, including 8 primary and 16 recurrent tumors. After read alignment and quality-control filtering, DESeq2 was used to identify differentially expressed genes (DEG) between primary and recurrent tumors. Functional enrichment pathway analysis and gene set enrichment analysis was performed using "clusterProfiler" and "GSVA" R packages. TME composition was investigated through the analysis and integration of multiple published RNA-seq deconvolution algorithms. TME analysis results were externally validated using data from independent previously published RNA-seq datasets. A total of 31 DEGs were identified between primary and recurrent aGCTs. These included genes with known function in hormone signaling such as LHCGR and INSL3 (more abundant in primary tumors) and CYP19A1 (more abundant in recurrent tumors). Gene set enrichment analysis revealed that primarily immune-related and hormone-regulated gene sets expression was increased in recurrent tumors. Integrative TME analysis demonstrated statistically significant depletion of cancer-associated fibroblasts in recurrent tumors. This finding was confirmed in multiple independent datasets.

IMPLICATIONS

Recurrent aGCTs exhibit alterations in hormone pathway gene expression as well as decreased infiltration of cancer-associated fibroblasts, suggesting dual roles for hormonal signaling and TME remodeling underpinning disease relapse.

Transcriptomic Profiling of Gene Expression Associated with Granulosa Cell Tumor Development in a Mouse Model

Ovarian granulosa cell tumors (GCTs) are rare sex cord-stromal tumors, accounting for ~5% ovarian tumors. The etiology of GCTs remains poorly defined. Genetically engineered mouse models are potentially valuable for understanding the pathogenesis of GCTs. Mice harboring constitutively active TGFβ signaling (TGFBR1-CA) develop ovarian GCTs that phenocopy several hormonal and molecular characteristics of human GCTs. To determine molecular alterations in the ovary upon TGFβ signaling activation, we performed transcriptomic profiling of gene expression associated with GCT development using ovaries from 1-month-old TGFBR1-CA mice and age-matched controls. RNA-sequencing and bioinformatics analysis coupled with the validation of select target genes revealed dysregulations of multiple cellular events and signaling molecules/pathways. The differentially expressed genes are enriched not only for known GCT-related pathways and tumorigenic events but also for signaling events potentially mediated by neuroactive ligand-receptor interaction, relaxin signaling, insulin signaling, and complements in TGFBR1-CA ovaries. Additionally, a comparative analysis of our data in mice with genes dysregulated in human GCTs or granulosa cells overexpressing a mutant FOXL2, the genetic hallmark of adult GCTs, identified some common genes altered in both conditions. In summary, this study has revealed the molecular signature of ovarian GCTs in a mouse model that harbors the constitutive activation of TGFBR1. The findings may be further exploited to understand the pathogenesis of a class of poorly defined ovarian tumors.

The Role of Omics Approaches to Characterize Molecular Mechanisms of Rare Ovarian Cancers: Recent Advances and Future Perspectives

Rare ovarian cancers are ovarian cancers with an annual incidence of less than 6 cases per 100,000 women. They generally have a poor prognosis due to being delayed diagnosis and treatment. Exploration of molecular mechanisms in these cancers has been challenging due to their rarity and research efforts being fragmented across the world. Omics approaches can provide detailed molecular snapshots of the underlying mechanisms of these cancers. Omics approaches, including genomics, transcriptomics, proteomics, and metabolomics, can identify potential candidate biomarkers for diagnosis, prognosis, and screening of rare gynecological cancers and can aid in identifying therapeutic targets. The integration of multiple omics techniques using approaches such as proteogenomics can provide a detailed understanding of the molecular mechanisms of carcinogenesis and cancer progression. Further, omics approaches can provide clues towards developing immunotherapies, cancer recurrence, and drug resistance in tumors; and form a platform for personalized medicine. The current review focuses on the application of omics approaches and integrative biology to gain a better understanding of rare ovarian cancers.

FOXO1 mitigates the SMAD3/FOXL2C134W transcriptomic effect in a model of human adult granulosa cell tumor

BACKGROUND

Adult granulosa cell tumor (aGCT) is a rare type of stromal cell malignant cancer of the ovary characterized by elevated estrogen levels. aGCTs ubiquitously harbor a somatic mutation in FOXL2 gene, Cys134Trp (c.402C < G); however, the general molecular effect of this mutation and its putative pathogenic role in aGCT tumorigenesis is not completely understood. We previously studied the role of FOXL2C134W, its partner SMAD3 and its antagonist FOXO1 in cellular models of aGCT.

METHODS

In this work, seeking more comprehensive profiling of FOXL2C134W transcriptomic effects, we performed an RNA-seq analysis comparing the effect of FOXL2WT/SMAD3 and FOXL2C134W/SMAD3 overexpression in an established human GC line (HGrC1), which is not luteinized, and bears normal alleles of FOXL2.

RESULTS

Our data shows that FOXL2C134W/SMAD3 overexpression alters the expression of 717 genes. These genes include known and novel FOXL2 targets (TGFB2, SMARCA4, HSPG2, MKI67, NFKBIA) and are enriched for neoplastic pathways (Proteoglycans in Cancer, Chromatin remodeling, Apoptosis, Tissue Morphogenesis, Tyrosine Kinase Receptors). We additionally expressed the FOXL2 antagonistic Forkhead protein, FOXO1. Surprisingly, overexpression of FOXO1 mitigated 40% of the altered genome-wide effects specifically related to FOXL2C134W, suggesting it can be a new target for aGCT treatment.

CONCLUSIONS

Our transcriptomic data provide novel insights into potential genes (FOXO1 regulated) that could be used as biomarkers of efficacy in aGCT patients.

Ovarian sex cord-stromal tumors: an update on clinical features, molecular changes, and management

Sex cord stromal-tumors are rare tumors of the ovary that include numerous tumor subtypes of variable histological features and biological behavior. Surgery is the main therapeutic modality for the management of these tumors, while chemotherapy and hormonal therapy may be used in some patients with progressive and recurrent tumors. Several studies investigated molecular changes in the different tumor types. Understanding molecular changes underlying the development and progression of sex cord-stromal tumors provides valuable information for diagnostic and prognostic biomarkers and potential therapeutic targets for these tumors. In this review, we provide an update on the clinical presentation, molecular changes, and management of sex cord-stromal tumors.

Functional Profiling of FSH and Estradiol in Ovarian Granulosa Cell Tumors

Adult-type granulosa cell tumors (AGCTs) are sex-cord derived neoplasms with a propensity for late relapse. Hormonal modulators have been used empirically in the treatment of recurrent AGCT, albeit with limited success. To provide a more rigorous foundation for hormonal therapy in AGCT, we used a multimodal approach to characterize the expressions of key hormone biomarkers in 175 tumor specimens and 51 serum samples using RNA sequencing, immunohistochemistry, RNA in situ hybridization, quantitative PCR, and circulating biomarker analysis, and correlated these results with clinical data. We show that FSH receptor and estrogen receptor beta (ERβ) are highly expressed in the majority of AGCTs, whereas the expressions of estrogen receptor alpha (ERα) and G-protein coupled estrogen receptor 1 are less prominent. ERβ protein expression is further increased in recurrent tumors. Aromatase expression levels show high variability between tumors. None of the markers examined served as prognostic biomarkers for progression-free or overall survival. In functional experiments, we assessed the effects of FSH, estradiol (E2), and the aromatase inhibitor letrozole on AGCT cell viability using 2 in vitro models: KGN cells and primary cultures of AGCT cells. FSH increased cell viability in a subset of primary AGCT cells, whereas E2 had no effect on cell viability at physiological concentrations. Letrozole suppressed E2 production in AGCTs; however, it did not impact cell viability. We did not find preclinical evidence to support the clinical use of aromatase inhibitors in AGCT treatment, and thus randomized, prospective clinical studies are needed to clarify the role of hormonal treatments in AGCTs.

A modified Fox pentagon technique performed using a polytetrafluoroethylene sling in frontalis suspension to treat blepharophimosis syndrome

The purpose of this study was to evaluate the functional and cosmetic outcomes of a new modified Fox pentagon technique performed using polytetrafluoroethylene in frontalis suspension surgery for blepharophimosis syndrome. This retrospective study enrolled 40 patients diagnosed with blepharophimosis syndrome from March 2016 to October 2018. All patients underwent frontalis suspension using a new modified Fox pentagon technique. The functional and cosmetic outcomes were evaluated. After the operation, the mean palpebral fissure height increased from 2.68 to 6.93 in right eyes and from 2.73 to 6.98 in left eyes. The mean MRD1 increased from 0.53 to 3.76 in right eyes and from 0.50 to 3.78 in left eyes. While preoperative to postoperative differences were statistically significant (p < 0.01), there were no significant differences between right and left eyes either before or after the surgery (p > 0.01). All patients achieved good cosmetic results with an average score of 0.6. We have experimentally created a modified Fox pentagon technique performed using a polytetrafluoroethylene sling in a frontalis suspension to treat BPES; this approach yielded favorable cosmetic and functional outcomes.

Comparative transcriptome analysis of matched primary and distant metastatic ovarian carcinoma

BACKGROUND

High grade serous ovarian carcinoma (HGSOC) is the most common subtype of epithelial ovarian cancers (EOC) with poor prognosis. In most cases EOC is widely disseminated at the time of diagnosis. Despite the optimal cytoreductive surgery and chemotherapy most patients develop chemoresistance, and the 5-year overall survival being only 25-35%.

METHODS

Here we analyzed the gene expression profiles of 10 primary HGSOC tumors and 10 related omental metastases using RNA sequencing and identified 100 differentially expressed genes.

RESULTS

The differentially expressed genes were associated with decreased embryogenesis and vasculogenesis and increased cellular proliferation and organismal death. Top upstream regulators responsible for this gene signature were NR5A1, GATA4, FOXL2, TP53 and BMP7. A subset of these genes were highly expressed in the ovarian cancer among the cancer transcriptomes of The Cancer Genome Atlas. Importantly, the metastatic gene signature was suggestive of poor survival in TCGA data based on gene enrichment analysis.

CONCLUSION

By comparing the gene expression profiles of primary HGSOC tumors and their matched metastasis, we provide evidence that a signature of 100 genes is able to separate these two sample types and potentially predict patient survival. Our study identifies functional categories of genes and transcription factors that could play important roles in promoting metastases and serve as markers for cancer prognosis.

The transcriptional regulator CBX2 and ovarian function: A whole genome and whole transcriptome approach

The chromobox homolog 2 (CBX2) was found to be important for human testis development, but its role in the human ovary remains elusive. We conducted a genome-wide analysis based on DNA adenine methyltransferase identification (DamID) and RNA sequencing strategies to investigate CBX2 in the human granulosa cells. Functional analysis revealed that CBX2 was upstream of genes contributing to ovarian function like folliculogenesis and steroidogenesis (i.e. ESR1, NRG1, AKR1C1, PTGER2, BMP15, BMP2, FSHR and NTRK1/2). We identified CBX2 regulated genes associated with polycystic ovary syndrome (PCOS) such as TGFβ, MAP3K15 and DKK1, as well as genes implicated in premature ovarian failure (POF) (i.e. POF1B, BMP15 and HOXA13) and the pituitary deficiency (i.e. LHX4 and KISS1). Our study provided an excellent opportunity to identify genes surrounding CBX2 in the ovary and might contribute to the understanding of ovarian physiopathology causing infertility in women.

Combined PPARγ Activation and XIAP Inhibition as a Potential Therapeutic Strategy for Ovarian Granulosa Cell Tumors

Ovarian granulosa cell tumors (GCT) are characterized by indolent growth and late relapse. No therapeutic modalities aside from surgery have proven effective. We previously reported overexpression of the nuclear receptor, peroxisome proliferator-activated receptor-gamma (PPARγ), and constitutive activity of the NFκB and AP1 signaling pathways in GCT. PPARγ presents as a potential therapeutic target as it impedes proliferation and promotes terminal differentiation of granulosa cells. However, resistance to the actions of PPARγ is caused by NFκB transrepression in GCT-derived cell lines, KGN and COV434. We showed that abrogation of NFκB signaling in GCT cells enables PPARγ agonists to initiate apoptosis. In addition, we observed overexpression of an NFκB-induced gene, X-linked inhibitor of apoptosis protein (XIAP), in GCT and GCT-derived cells. XIAP is an attractive therapeutic target due to its role in inhibiting the apoptotic pathway. We investigated the antitumor effects of combined XIAP inhibition using Smac-mimetics and PPARγ activation using thiazolidinediones (TZD) in the GCT-derived cells. Transactivation assays revealed that NFκB transrepression of PPARγ can be relieved by NFκB or XIAP inhibition. Combined Smac-mimetic and TZD significantly induced apoptosis, reduced cell viability and proliferation in KGN cells in monolayer and 3D spheroid culture, and in GCT explant models. The Smac-mimetic and TZD cotreatment also delayed cell invasion, upregulated proapoptotic genes, and compromised cell metabolism in KGN cells. This study provides evidence that PPARγ and XIAP cotreatment has antineoplastic effects in GCT. As therapeutics that target these proteins are already in clinical or preclinical use, expedient translation to the clinic is possible.

FOXL2C134W-Induced CYP19 Expression via Cooperation With SMAD3 in HGrC1 Cells

Germline knockout studies in female mice demonstrated an essential role for forkhead box L2 (FOXL2) in early follicle development, whereas an inducible granulosa cell (GC)-specific deletion of Foxl2 in adults has shown ovary-to-testis somatic sex reprogramming. In women, over 120 different germline mutations in the FOXL2 gene have been shown to cause blepharophimosis/ptosis/epicantus inversus syndrome associated with or without primary ovarian insufficiency. By contrast, a single somatic mutation (FOXL2C134W) accounts for almost all adult-type GC tumors (aGCTs). To test the hypothesis that FOXL2C134W differentially regulates the expression of aGCT markers, we investigated the effect of FOXL2C134W on inhibin B and P450 aromatase expression using a recently established human GC line (HGrC1), which we now show to bear two normal alleles of FOXL2. Neither FOXL2wt nor FOXL2C134W regulate INHBB messenger RNA (mRNA) expression. However, FOXL2C134W selectively displays a 50-fold induction of CYP19 mRNA expression dependent upon activin A. Mechanistically, the CYP19 promoter is activated in a similar way by FOXL2C134W interaction with SMAD3, but not by FOXL2wt. SMAD2 had no effect. Moreover, FOXL2C134W interactions with SMAD3 and with the FOX binding element located at -199 bp upstream of the ATG initiation codon of CYP19 are more sustainable than FOXL2wt. Thus, FOXL2C134W potentiates CYP19 expression in HGrC1 cells via enhanced recruitment of SMAD3 to a proximal FOX binding element. These findings may explain the pathophysiology of estrogen excess in patients with aGCT.

Functional Analysis of a Novel FOXL2 Indel Mutation in Chinese Families with Blepharophimosis-Ptosis-Epicanthus Inversus Syndrome Type I

Background: Blepharophimosis-ptosis-epicanthus inversus syndrome (BPES) is an autosomal dominant disease with a low incidence rate. Indel mutations in the forkhead box L2 (FOXL2) gene cause two types of BPES that are distinguished by the presence (type I) or absence (type II) of premature ovarian failure (POF). The purpose of this study was to identify a possible deletion in FOXL2 in Chinese families with BPES and to clarify its relationship with POF.

Methods: An autosomal dominant Chinese BPES family with four generations was enrolled in this study. Peripheral venous blood was collected from all affected patients, and genomic DNA was extracted from leukocytes. The whole coding sequence and nearby 5' untranslated region (UTR) and 3'UTR of the FOXL2 gene were amplified using polymerase chain reaction (PCR) with three sets of overlapping primers, followed by sequencing analyses. The sequencing results were analysed using SeqMan software. Based on the patients' clinical manifestations and analysis of the identified indel mutation, we found that the mutation disturbed interactions between FOXL2 and the StAR gene. Furthermore, through subcellular localisation and functional studies, we observed significant mislocalisation of the mutant protein; the mutant protein was found in the cytoplasm, while the wild-type protein was found in the nucleus. Loss of function was confirmed by transcriptional activity assays, quantitative real-time PCR, and electrophoretic mobility shift assays.

Results: All affected patients presented with clinical features of BPES type I, including small palpebral fissures, ptosis, telecanthus, and epicanthus inversus with POF. A novel FOXL2 heterozygous indel mutation, c.19_95del, a 77-bp deletion that disrupts FOXL2 protein structure, was identified in all affected members of the family. In addition, this indel mutation significantly increased StAR mRNA expression by disrupting the ability of the FOXL2 protein to bind to the StAR promoter and act as a repressor of this gene.

Conclusions: A novel FOXL2 indel mutation was identified in Chinese families with BPES. Our results expand the spectrum of known FOXL2 mutations and provide additional insight into the structure-function relationships of the FOXL2 protein. Furthermore, this novel mutation resulted in the dysfunction of FOXL2 as a transcription factor, blocking its ability to bind to the promoter region of the StAR gene, resulting in POF in the affected patient.