Fallopian tube single cell analysis reveals myeloid cell alterations in high-grade serous ovarian cancer.

Most high-grade serous ovarian cancers (HGSCs) likely initiate from fallopian tube (FT) epithelia. While epithelial subtypes have been characterized using single-cell RNA-sequencing (scRNA-Seq), heterogeneity of other compartments and their involvement in tumor progression are poorly defined. Integrated analysis of human FT scRNA-Seq and HGSC-related tissues, including tumors, revealed greater immune and stromal transcriptional diversity than previously reported. We identified abundant monocytes in FTs across two independent cohorts. The ratio of macrophages to monocytes is similar between benign FTs, ovaries, and adjacent normal tissues but significantly greater in tumors. FT-defined monocyte and macrophage signatures, cell-cell communication, and gene set enrichment analyses identified monocyte- and macrophage-specific interactions and functional pathways in paired tumors and adjacent normal tissues. Further reanalysis of HGSC scRNA-Seq identified monocyte and macrophage subsets associated with neoadjuvant chemotherapy. Taken together, our work provides data that an altered FT myeloid cell composition could inform the discovery of early detection markers for HGSC.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

Rewiring of master transcription factor cistromes during high-grade serous ovarian cancer development.

The transcription factors MECOM, PAX8, SOX17 and WT1 are candidate master regulators of high-grade serous 'ovarian' cancer (HGSC), yet their cooperative role in the hypothesized tissue of origin, the fallopian tube secretory epithelium (FTSEC) is unknown. We generated 26 epigenome (CUT&TAG, CUT&RUN, ATAC-seq and HiC) data sets and 24 profiles of RNA-seq transcription factor knock-down followed by RNA sequencing in FTSEC and HGSC models to define binding sites and gene sets regulated by these factors in cis and trans. This revealed that MECOM, PAX8, SOX17 and WT1 are lineage-enriched, super-enhancer associated master regulators whose cooperative DNA-binding patterns and target genes are re-wired during tumor development. All four TFs were indispensable for HGSC clonogenicity and survival but only depletion of PAX8 and WT1 impaired FTSEC cell survival. These four TFs were pharmacologically inhibited by transcriptional inhibitors only in HGSCs but not in FTSECs. Collectively, our data highlights that tumor-specific epigenetic remodeling is tightly related to MECOM, PAX8, SOX17 and WT1 activity and these transcription factors are targetable in a tumor-specific manner through transcriptional inhibitors.

Single-cell transcriptomic analysis of endometriosis.

Endometriosis is a common condition in women that causes chronic pain and infertility and is associated with an elevated risk of ovarian cancer. We profiled transcriptomes of >370,000 individual cells from endometriomas (n = 8), endometriosis (n = 28), eutopic endometrium (n = 10), unaffected ovary (n = 4) and endometriosis-free peritoneum (n = 4), generating a cellular atlas of endometrial-type epithelial cells, stromal cells and microenvironmental cell populations across tissue sites. Cellular and molecular signatures of endometrial-type epithelium and stroma differed across tissue types, suggesting a role for cellular restructuring and transcriptional reprogramming in the disease. Epithelium, stroma and proximal mesothelial cells of endometriomas showed dysregulation of pro-inflammatory pathways and upregulation of complement proteins. Somatic ARID1A mutation in epithelial cells was associated with upregulation of pro-angiogenic and pro-lymphangiogenic factors and remodeling of the endothelial cell compartment, with enrichment of lymphatic endothelial cells. Finally, signatures of ciliated epithelial cells were enriched in ovarian cancers, reinforcing epidemiologic associations between these two diseases.

circCsnk1g3- and circAnkib1-regulated interferon responses in sarcoma promote tumorigenesis by shaping the immune microenvironment.

Exonic circular RNAs (circRNAs) produce predominantly non-coding RNA species that have been recently profiled in many tumors. However, their functional contribution to cancer progression is still poorly understood. Here, we identify the circRNAs expressed in soft tissue sarcoma cells and explore how the circRNAs regulate sarcoma growth in vivo. We show that circCsnk1g3 and circAnkib1 promote tumor growth by shaping a pro-tumorigenic microenvironment, possibly due to their capabilities to regulate tumor-promoting elements extrinsic to the tumor cells. Accordingly, circCsnk1g3 and circAnkib1 can control the expression of interferon-related genes and pro-inflammatory factors in the sarcoma cells, thus directing immune cell recruitment into the tumor mass, and hence their activation. Mechanistically, circRNAs may repress pro-inflammatory elements by buffering activation of the pathways mediated by RIG-I, the cytosolic viral RNA sensor. The current findings suggest that the targeting of specific circRNAs could augment the efficacy of tumor and immune response to mainstay therapies.

chromMAGMA: regulatory element-centric interrogation of risk variants.

Candidate causal risk variants from genome-wide association studies reside almost exclusively in noncoding regions of the genome and innovative approaches are necessary to understand their biological function. Multi-marker analysis of genomic annotation (MAGMA) is a widely used program that nominates candidate risk genes by mapping single-nucleotide polymorphism summary statistics from genome-wide association studies to gene bodies. We augmented MAGMA to create chromatin-MAGMA (chromMAGMA), a method to nominate candidate risk genes based on the presence of risk variants within noncoding regulatory elements (REs). We applied chromMAGMA to a genetic susceptibility dataset for epithelial ovarian cancer (EOC), a rare gynecologic malignancy characterized by high mortality. This identified 155 unique candidate EOC risk genes across five EOC histotypes; 83% (105/127) of high-grade serous ovarian cancer risk genes had not previously been implicated in this EOC histotype. Risk genes nominated by chromMAGMA converged on mRNA splicing and transcriptional dysregulation pathways. chromMAGMA is a pipeline that nominates candidate risk genes through a gene regulation-focused approach and helps interpret the biological mechanism of noncoding risk variants for complex diseases.

Predicting master transcription factors from pan-cancer expression data.

Critical developmental "master transcription factors" (MTFs) can be subverted during tumorigenesis to control oncogenic transcriptional programs. Current approaches to identifying MTFs rely on ChIP-seq data, which is unavailable for many cancers. We developed the CaCTS (Cancer Core Transcription factor Specificity) algorithm to prioritize candidate MTFs using pan-cancer RNA sequencing data. CaCTS identified candidate MTFs across 34 tumor types and 140 subtypes including predictions for cancer types/subtypes for which MTFs are unknown, including e.g. PAX8, SOX17, and MECOM as candidates in ovarian cancer (OvCa). In OvCa cells, consistent with known MTF properties, these factors are required for viability, lie proximal to superenhancers, co-occupy regulatory elements globally, co-bind loci encoding OvCa biomarkers, and are sensitive to pharmacologic inhibition of transcription. Our predictions of MTFs, especially for tumor types with limited understanding of transcriptional drivers, pave the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

Single-cell transcriptomics identifies gene expression networks driving differentiation and tumorigenesis in the human fallopian tube.

The human fallopian tube harbors the cell of origin for the majority of high-grade serous "ovarian" cancers (HGSCs), but its cellular composition, particularly the epithelial component, is poorly characterized. We perform single-cell transcriptomic profiling of around 53,000 individual cells from 12 primary fallopian specimens to map their major cell types. We identify 10 epithelial subpopulations with diverse transcriptional programs. Based on transcriptional signatures, we reconstruct a trajectory whereby secretory cells differentiate into ciliated cells via a RUNX3high intermediate. Computational deconvolution of advanced HGSCs identifies the "early secretory" population as a likely precursor state for the majority of HGSCs. Its signature comprises both epithelial and mesenchymal features and is enriched in mesenchymal-type HGSCs (p = 6.7 × 10-27), a group known to have particularly poor prognoses. This cellular and molecular compendium of the human fallopian tube in cancer-free women is expected to advance our understanding of the earliest stages of fallopian epithelial neoplasia.

Non-coding somatic mutations converge on the PAX8 pathway in ovarian cancer.

The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineage-specific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

A Study of High-Grade Serous Ovarian Cancer Origins Implicates the SOX18 Transcription Factor in Tumor Development.

Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (padj < 8 × 10-4). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (padj < 2 × 10-4), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10-16). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs.

Super-Enhancer-Associated LncRNA UCA1 Interacts Directly with AMOT to Activate YAP Target Genes in Epithelial Ovarian Cancer.

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of tumorigenesis, and yet their mechanistic roles remain challenging to characterize. Here, we integrate functional proteomics with lncRNA-interactome profiling to characterize Urothelial Cancer Associated 1 (UCA1), a candidate driver of ovarian cancer development. Reverse phase protein array (RPPA) analysis indicates that UCA1 activates transcription coactivator YAP and its target genes. In vivo RNA antisense purification (iRAP) of UCA1 interacting proteins identified angiomotin (AMOT), a known YAP regulator, as a direct binding partner. Loss-of-function experiments show that AMOT mediates YAP activation by UCA1, as UCA1 enhances the AMOT-YAP interaction to promote YAP dephosphorylation and nuclear translocation. Together, we characterize UCA1 as a lncRNA regulator of Hippo-YAP signaling and highlight the UCA1-AMOT-YAP signaling axis in ovarian cancer development.

A transcriptome-wide association study of high-grade serous epithelial ovarian cancer identifies new susceptibility genes and splice variants.

We sought to identify susceptibility genes for high-grade serous ovarian cancer (HGSOC) by performing a transcriptome-wide association study of gene expression and splice junction usage in HGSOC-relevant tissue types (N = 2,169) and the largest genome-wide association study available for HGSOC (N = 13,037 cases and 40,941 controls). We identified 25 transcriptome-wide association study significant genes, 7 at the junction level only, including LRRC46 at 19q21.32, (P = 1 × 10-9), CHMP4C at 8q21 (P = 2 × 10-11) and a PRC1 junction at 15q26 (P = 7 × 10-9). In vitro assays for CHMP4C showed that the associated variant induces allele-specific exon inclusion (P = 0.0024). Functional screens in HGSOC cell lines found evidence of essentiality for three of the new genes we identified: HAUS6, KANSL1 and PRC1, with the latter comparable to MYC. Our study implicates at least one target gene for 6 out of 13 distinct genome-wide association study regions, identifying 23 new candidate susceptibility genes for HGSOC.

RNA editing in nascent RNA affects pre-mRNA splicing.

In eukaryotes, nascent RNA transcripts undergo an intricate series of RNA processing steps to achieve mRNA maturation. RNA editing and alternative splicing are two major RNA processing steps that can introduce significant modifications to the final gene products. By tackling these processes in isolation, recent studies have enabled substantial progress in understanding their global RNA targets and regulatory pathways. However, the interplay between individual steps of RNA processing, an essential aspect of gene regulation, remains poorly understood. By sequencing the RNA of different subcellular fractions, we examined the timing of adenosine-to-inosine (A-to-I) RNA editing and its impact on alternative splicing. We observed that >95% A-to-I RNA editing events occurred in the chromatin-associated RNA prior to polyadenylation. We report about 500 editing sites in the 3' acceptor sequences that can alter splicing of the associated exons. These exons are highly conserved during evolution and reside in genes with important cellular function. Furthermore, we identified a second class of exons whose splicing is likely modulated by RNA secondary structures that are recognized by the RNA editing machinery. The genome-wide analyses, supported by experimental validations, revealed remarkable interplay between RNA editing and splicing and expanded the repertoire of functional RNA editing sites.

Genomic analysis of ADAR1 binding and its involvement in multiple RNA processing pathways.

Adenosine deaminases acting on RNA (ADARs) are the primary factors underlying adenosine to inosine (A-to-I) editing in metazoans. Here we report the first global study of ADAR1-RNA interaction in human cells using CLIP-seq. A large number of CLIP sites are observed in Alu repeats, consistent with ADAR1's function in RNA editing. Surprisingly, thousands of other CLIP sites are located in non-Alu regions, revealing functional and biophysical targets of ADAR1 in the regulation of alternative 3' UTR usage and miRNA biogenesis. We observe that binding of ADAR1 to 3' UTRs precludes binding by other factors, causing 3' UTR lengthening. Similarly, ADAR1 interacts with DROSHA and DGCR8 in the nucleus and possibly out-competes DGCR8 in primary miRNA binding, which enhances mature miRNA expression. These functions are dependent on ADAR1's editing activity, at least for a subset of targets. Our study unfolds a broad landscape of the functional roles of ADAR1.

Latency-associated nuclear antigen of Kaposi sarcoma-associated herpesvirus promotes angiogenesis through targeting notch signaling effector Hey1.

Notch signaling has been implicated in the pathogenesis of Kaposi sarcoma. Kaposi sarcoma is an angioproliferative neoplasm that originates from Kaposi sarcoma-associated herpesvirus (KSHV) infection. Previously, we showed that the KSHV LANA protein can stabilize intracellular Notch in KSHV-infected tumor cells and promote cell proliferation. However, whether Notch signaling functions in pathologic angiogenesis of Kaposi sarcoma remains largely unknown. Hey1, an essential downstream effector of the Notch signaling pathway, has been demonstrated to play a fundamental role in vascular development. In the present study, we performed whole transcriptome, paired-end sequencing on three patient-matched clinical Kaposi sarcoma specimens and their corresponding adjacent stroma samples, with an average depth of 42 million reads per sample. Dll4, Hey1, and HeyL displayed significant upregulation in Kaposi sarcoma. Further verification based on immunohistochemistry analysis demonstrated that Hey1 was indeed highly expressed in Kaposi sarcoma lesions. Using the Matrigel plug assay, we showed that downregulation of Hey1 and γ-secretase inhibitor treatment caused dramatic reduction in the formation of new blood vessels in mice. Interestingly, LANA was responsible for the elevated level of Hey1 through inhibition of its degradation. Importantly, Hey1 stabilized by LANA promoted the neoplastic vasculature. Taken together, our data suggest that hijacking of the proangiogenic property of Hey1 by LANA is an important strategy utilized by KSHV to achieve pathologic angiogenesis and that Hey1 is a potential therapeutic target in Kaposi sarcoma.

MicroRNAs and unusual small RNAs discovered in Kaposi's sarcoma-associated herpesvirus virions.

It is widely held that any given virus uses only one type of nucleic acid for genetic information storage. However, this consensus has been challenged slightly by several recent studies showing that many RNA species are present within a range of DNA viruses that include Kaposi's sarcoma-associated herpesvirus (KSHV). RNAs extracted from purified DNA virus particles exhibit great diversity in terms of length, abundance, temporal expression, cellular localization, and coding capacity during viral infection. In addition to known RNA species, the current study showed that small regulatory RNAs were present in KSHV virions. A large number of viral and cellular microRNAs (miRNAs), as well as unusual small RNAs (usRNAs), were detected in KSHV virions by using deep sequencing. Both viral and host miRNAs detected in small RNAs extracted from KSHV virions were further shown to colocalize with KSHV virions directly by in situ hybridization (ISH)-electron microscopy (EM) (ISH-EM). Some of these miRNAs were differentially present in the host cells and KSHV virions, suggesting that they are not randomly present in KSHV virions. The virional miRNAs could be transported into host cells, and they are biologically functional during de novo viral infection. Our study revealed miRNAs and usRNAs as a novel group of components in KSHV virions.

Kaposi's sarcoma-associated herpesvirus-encoded microRNA miR-K12-11 attenuates transforming growth factor beta signaling through suppression of SMAD5.

Kaposi's sarcoma-associated herpesvirus (KSHV) encodes 12 pre-microRNAs (pre-miRNAs). Current studies have shown that these miRNAs are involved in regulation of viral and host gene expression, implicating a role in the maintenance of viral latency and suppression of antiviral innate immunity. However, the functions of these miRNAs remain largely unknown. On the basis of the sequence homology between oncogenic miR-155 and KSHV-encoded miR-K12-11, we hypothesized that miR-K12-11 could attenuate transforming growth factor ß (TGF-ß) signaling, facilitating viral infection and tumorigenesis. In the present study, we demonstrated that ectopic expression of miR-K12-11 in Ramos, a TGF-ß-sensitive cell line, downregulated TGF-ß signaling and facilitated cell proliferation upon TGF-ß treatment by directly targeting SMAD5, an important mediator in TGF-ß signaling. In addition, the downregulation of SMAD5 by miR-K12-11 was further confirmed in a de novo KSHV infection system or latently infected KSHV-positive B-lymphoma cell lines. More importantly, repression of miR-K12-11 by a specific sponge inhibitor restored the expression of SMAD5 in both de novo-infected and latently infected cells. Finally, we found that restoration of SMAD5, in addition to the TGF-ß type II receptor, which was epigenetically silenced by the latent viral protein latency-associated nuclear antigen, sensitized BC3 cells to the cytostatic effect of TGF-ß signaling. Taken together, our findings highlight a novel mechanism in which miR-K12-11 downregulates TGF-ß signaling and suggest that viral miRNAs and proteins may exert a dichotomy regulation in virus-induced oncogenesis by targeting the same signaling pathway.

Kaposi's sarcoma-associated herpesvirus-encoded latency-associated nuclear antigen reduces interleukin-8 expression in endothelial cells and impairs neutrophil chemotaxis by degrading nuclear p65.

Latency-associated nuclear antigen 1 (LANA-1) of Kaposi's sarcoma-associated herpesvirus (KSHV) is the major viral latent protein and functions as a multifaceted protein. Here, we report that LANA-1 attenuates the endothelial response to tumor necrosis factor alpha (TNF-α) stimulation and inhibits consequent neutrophil chemotaxis. Reporter assays showed that LANA-1 constantly repressed nuclear factor (NF)-κB transactivity upon TNF-α stimulation. We also found that LANA-1 decreased nuclear p65 protein levels through enhancement of polyubiquitylation-mediated p65 degradation and that an elongin B/elongin C-cullin 5-LANA-1-p65 complex assembled by LANA-1 was responsible for this enhanced p65 degradation. In telomerase-immortalized human umbilical vein endothelial cells, LANA-1 was demonstrated to repress interleukin-8 expression, which was involved in neutrophil recruitment to the inflammatory site. Through an in vitro transmigration assay, we determined a suppressive effect of LANA-1 on neutrophil chemotaxis. Our work suggests that KSHV LANA-1 is a negative modulator of acute inflammation and sheds light on a new mechanism by which KSHV during the latent life cycle evades the host innate immune response.