Integrative multi-omic cancer profiling reveals DNA methylation patterns associated with therapeutic vulnerability and cell-of-origin.

DNA methylation plays a critical role in establishing and maintaining cellular identity. However, it is frequently dysregulated during tumor development and is closely intertwined with other genetic alterations. Here, we leveraged multi-omic profiling of 687 tumors and matched non-involved adjacent tissues from the kidney, brain, pancreas, lung, head and neck, and endometrium to identify aberrant methylation associated with RNA and protein abundance changes and build a Pan-Cancer catalog. We uncovered lineage-specific epigenetic drivers including hypomethylated FGFR2 in endometrial cancer. We showed that hypermethylated STAT5A is associated with pervasive regulon downregulation and immune cell depletion, suggesting that epigenetic regulation of STAT5A expression constitutes a molecular switch for immunosuppression in squamous tumors. We further demonstrated that methylation subtype-enrichment information can explain cell-of-origin, intra-tumor heterogeneity, and tumor phenotypes. Overall, we identified cis-acting DNA methylation events that drive transcriptional and translational changes, shedding light on the tumor's epigenetic landscape and the role of its cell-of-origin.

VISTA is an activating receptor in human monocytes.

As indicated by its name, V-domain Ig suppressor of T cell activation (VISTA) is thought to serve primarily as an inhibitory protein that limits immune responses. VISTA antibodies can dampen the effects of several concomitantly elicited activation signals, including TCR and TLR activation, but it is currently unclear if VISTA agonism could singly affect immune cell biology. In this study, we discovered two novel VISTA antibodies and characterized their effects on human peripheral blood mononuclear cells by scRNA/CITE-seq. Both antibodies appeared to agonize VISTA in an Fc-functional manner to elicit transcriptional and functional changes in monocytes consistent with activation. We also used pentameric VISTA to identify Syndecan-2 and several heparan sulfate proteoglycan synthesis genes as novel regulators of VISTA interactions with monocytic cells, adding further evidence of bidirectional signaling. Together, our study highlights several novel aspects of VISTA biology that have yet to be uncovered in myeloid cells and serves as a foundation for future research.

Co-evolution of tumor and immune cells during progression of multiple myeloma.

Multiple myeloma (MM) is characterized by the uncontrolled proliferation of plasma cells. Despite recent treatment advances, it is still incurable as disease progression is not fully understood. To investigate MM and its immune environment, we apply single cell RNA and linked-read whole genome sequencing to profile 29 longitudinal samples at different disease stages from 14 patients. Here, we collect 17,267 plasma cells and 57,719 immune cells, discovering patient-specific plasma cell profiles and immune cell expression changes. Patients with the same genetic alterations tend to have both plasma cells and immune cells clustered together. By integrating bulk genomics and single cell mapping, we track plasma cell subpopulations across disease stages and find three patterns: stability (from precancer to diagnosis), and gain or loss (from diagnosis to relapse). In multiple patients, we detect "B cell-featured" plasma cell subpopulations that cluster closely with B cells, implicating their cell of origin. We validate AP-1 complex differential expression (JUN and FOS) in plasma cell subpopulations using CyTOF-based protein assays, and integrated analysis of single-cell RNA and CyTOF data reveals AP-1 downstream targets (IL6 and IL1B) potentially leading to inflammation regulation. Our work represents a longitudinal investigation for tumor and microenvironment during MM progression and paves the way for expanding treatment options.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma.

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

Identification of murine gammaherpesvirus 68 miRNA-mRNA hybrids reveals miRNA target conservation among gammaherpesviruses including host translation and protein modification machinery.

Gammaherpesviruses, including the human pathogens Epstein-Barr virus (EBV) and Kaposi's sarcoma-associated herpesvirus (KSHV), establish lifelong latent infection in B cells and are associated with a variety of tumors. In addition to protein coding genes, these viruses encode numerous microRNAs (miRNAs) within their genomes. While putative host targets of EBV and KSHV miRNAs have been previously identified, the specific functions of these miRNAs during in vivo infection are largely unknown. Murine gammaherpesvirus 68 (MHV68) is a natural pathogen of rodents that is genetically related to both EBV and KSHV, and thus serves as an excellent model for the study of EBV and KSHV genetic elements such as miRNAs in the context of infection and disease. However, the specific targets of MHV68 miRNAs remain completely unknown. Using a technique known as qCLASH (quick crosslinking, ligation, and sequencing of hybrids), we have now identified thousands of Ago-associated, direct miRNA-mRNA interactions during lytic infection, latent infection and reactivation from latency. Validating this approach, detailed molecular analyses of specific interactions demonstrated repression of numerous host mRNA targets of MHV68 miRNAs, including Arid1a, Ctsl, Ifitm3 and Phc3. Notably, of the 1,505 MHV68 miRNA-host mRNA targets identified in B cells, 86% were shared with either EBV or KSHV, and 64% were shared among all three viruses, demonstrating significant conservation of gammaherpesvirus miRNA targeting. Pathway analysis of MHV68 miRNA targets further revealed enrichment of cellular pathways involved in protein synthesis and protein modification, including eIF2 Signaling, mTOR signaling and protein ubiquitination, pathways also enriched for targets of EBV and KSHV miRNAs. These findings provide substantial new information about specific targets of MHV68 miRNAs and shed important light on likely conserved functions of gammaherpesvirus miRNAs.

Computational analysis of ribonomics datasets identifies long non-coding RNA targets of γ-herpesviral miRNAs.

Ribonomics experiments involving crosslinking and immuno-precipitation (CLIP) of Ago proteins have expanded the understanding of the miRNA targetome of several organisms. These techniques, collectively referred to as CLIP-seq, have been applied to identifying the mRNA targets of miRNAs expressed by Kaposi's Sarcoma-associated herpes virus (KSHV) and Epstein-Barr virus (EBV). However, these studies focused on identifying only those RNA targets of KSHV and EBV miRNAs that are known to encode proteins. Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are also targeted by miRNAs. In this study, we performed a systematic re-analysis of published datasets from KSHV- and EBV-driven cancers. We used CLIP-seq data from lymphoma cells or EBV-transformed B cells, and a crosslinking, ligation and sequencing of hybrids dataset from KSHV-infected endothelial cells, to identify novel lncRNA targets of viral miRNAs. Here, we catalog the lncRNA targetome of KSHV and EBV miRNAs, and provide a detailed in silico analysis of lncRNA-miRNA binding interactions. Viral miRNAs target several hundred lncRNAs, including a subset previously shown to be aberrantly expressed in human malignancies. In addition, we identified thousands of lncRNAs to be putative targets of human miRNAs, suggesting that miRNA-lncRNA interactions broadly contribute to the regulation of gene expression.

Modified Cross-Linking, Ligation, and Sequencing of Hybrids (qCLASH) Identifies Kaposi's Sarcoma-Associated Herpesvirus MicroRNA Targets in Endothelial Cells.

Kaposi's sarcoma (KS) tumors are derived from endothelial cells and express Kaposi's sarcoma-associated herpesvirus (KSHV) microRNAs (miRNAs). Although miRNA targets have been identified in B cell lymphoma-derived cells and epithelial cells, little has been done to characterize the KSHV miRNA targetome in endothelial cells. A recent innovation in the identification of miRNA targetomes, cross-linking, ligation, and sequencing of hybrids (CLASH), unambiguously identifies miRNAs and their targets by ligating the two species while both species are still bound within the RNA-induced silencing complex (RISC). We developed a streamlined quick CLASH (qCLASH) protocol that requires a lower cell input than the original method and therefore has the potential to be used on patient biopsy samples. Additionally, we developed a fast-growing, KSHV-negative endothelial cell line derived from telomerase-immortalized vein endothelial long-term culture (TIVE-LTC) cells. qCLASH was performed on uninfected cells and cells infected with either wild-type KSHV or a mutant virus lacking miR-K12-11/11*. More than 1,400 cellular targets of KSHV miRNAs were identified. Many of the targets identified by qCLASH lacked a canonical seed sequence match. Additionally, most target regions in mRNAs originated from the coding DNA sequence (CDS) rather than the 3' untranslated region (UTR). This set of genes includes some that were previously identified in B cells and some new genes that warrant further study. Pathway analysis of endothelial cell targets showed enrichment in cell cycle control, apoptosis, and glycolysis pathways, among others. Characterization of these new targets and the functional consequences of their repression will be important in furthering our understanding of the role of KSHV miRNAs in oncogenesis.IMPORTANCE KS lesions consist of endothelial cells latently infected with KSHV. Cells that make up these lesions express KSHV miRNAs. Identification of the targets of KSHV miRNAs will help us understand their role in viral oncogenesis. The cross-linking and sequencing of hybrids (CLASH) protocol is a method for unambiguously identifying miRNA targetomes. We developed a streamlined version of CLASH, called quick CLASH (qCLASH). qCLASH requires a lower initial input of cells than for its parent protocol. Additionally, a new fast-growing KSHV-negative endothelial cell line, named TIVE-EX-LTC cells, was established. qCLASH was performed on TIVE-EX-LTC cells latently infected with wild-type (WT) KSHV or a mutant virus lacking miR-K12-11/11*. A number of novel targets of KSHV miRNAs were identified, including targets of miR-K12-11, the ortholog of the cellular oncogenic miRNA (oncomiR) miR-155. Many of the miRNA targets were involved in processes related to oncogenesis, such as glycolysis, apoptosis, and cell cycle control.

microRNA dependent and independent deregulation of long non-coding RNAs by an oncogenic herpesvirus.

Kaposi's sarcoma (KS) is a highly prevalent cancer in AIDS patients, especially in sub-Saharan Africa. Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiological agent of KS and other cancers like Primary Effusion Lymphoma (PEL). In KS and PEL, all tumors harbor latent KSHV episomes and express latency-associated viral proteins and microRNAs (miRNAs). The exact molecular mechanisms by which latent KSHV drives tumorigenesis are not completely understood. Recent developments have highlighted the importance of aberrant long non-coding RNA (lncRNA) expression in cancer. Deregulation of lncRNAs by miRNAs is a newly described phenomenon. We hypothesized that KSHV-encoded miRNAs deregulate human lncRNAs to drive tumorigenesis. We performed lncRNA expression profiling of endothelial cells infected with wt and miRNA-deleted KSHV and identified 126 lncRNAs as putative viral miRNA targets. Here we show that KSHV deregulates host lncRNAs in both a miRNA-dependent fashion by direct interaction and in a miRNA-independent fashion through latency-associated proteins. Several lncRNAs that were previously implicated in cancer, including MEG3, ANRIL and UCA1, are deregulated by KSHV. Our results also demonstrate that KSHV-mediated UCA1 deregulation contributes to increased proliferation and migration of endothelial cells.