Satellite repeat RNA expression in epithelial ovarian cancer associates with a tumor-immunosuppressive phenotype.

Aberrant expression of viral-like repeat elements is a common feature of epithelial cancers, and the substantial diversity of repeat species provides a distinct view of the cancer transcriptome. Repeatome profiling across ovarian, pancreatic, and colorectal cell lines identifies distinct clustering independent of tissue origin that is seen with coding gene analysis. Deeper analysis of ovarian cancer cell lines demonstrated that human satellite II (HSATII) satellite repeat expression was highly associated with epithelial-mesenchymal transition (EMT) and anticorrelated with IFN-response genes indicative of a more aggressive phenotype. SATII expression - and its correlation with EMT and anticorrelation with IFN-response genes - was also found in ovarian cancer RNA-Seq data and was associated with significantly shorter survival in a second independent cohort of patients with ovarian cancer. Repeat RNAs were enriched in tumor-derived extracellular vesicles capable of stimulating monocyte-derived macrophages, demonstrating a mechanism that alters the tumor microenvironment with these viral-like sequences. Targeting of HSATII with antisense locked nucleic acids stimulated IFN response and induced MHC I expression in ovarian cancer cell lines, highlighting a potential strategy of modulating the repeatome to reestablish antitumor cell immune surveillance.

Reverse Transcriptase Inhibition Disrupts Repeat Element Life Cycle in Colorectal Cancer.

Altered RNA expression of repetitive sequences and retrotransposition are frequently seen in colorectal cancer, implicating a functional importance of repeat activity in cancer progression. We show the nucleoside reverse transcriptase inhibitor 3TC targets activities of these repeat elements in colorectal cancer preclinical models with a preferential effect in p53-mutant cell lines linked with direct binding of p53 to repeat elements. We translate these findings to a human phase II trial of single-agent 3TC treatment in metastatic colorectal cancer with demonstration of clinical benefit in 9 of 32 patients. Analysis of 3TC effects on colorectal cancer tumorspheres demonstrates accumulation of immunogenic RNA:DNA hybrids linked with induction of interferon response genes and DNA damage response. Epigenetic and DNA-damaging agents induce repeat RNAs and have enhanced cytotoxicity with 3TC. These findings identify a vulnerability in colorectal cancer by targeting the viral mimicry of repeat elements. SIGNIFICANCE: Colorectal cancers express abundant repeat elements that have a viral-like life cycle that can be therapeutically targeted with nucleoside reverse transcriptase inhibitors (NRTI) commonly used for viral diseases. NRTIs induce DNA damage and interferon response that provide a new anticancer therapeutic strategy. This article is highlighted in the In This Issue feature, p. 1397.

Temporal and spatial heterogeneity of host response to SARS-CoV-2 pulmonary infection.

The relationship of SARS-CoV-2 pulmonary infection and severity of disease is not fully understood. Here we show analysis of autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter-patient and intra-patient heterogeneity of pulmonary virus infection. There is a spectrum of high and low virus cases associated with duration of disease. High viral cases have high activation of interferon pathway genes and a predominant M1-like macrophage infiltrate. Low viral cases are more heterogeneous likely reflecting inherent patient differences in the evolution of host response, but there is consistent indication of pulmonary epithelial cell recovery based on napsin A immunohistochemistry and RNA expression of surfactant and mucin genes. Using a digital spatial profiling platform, we find the virus corresponds to distinct spatial expression of interferon response genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Temporal and Spatial Heterogeneity of Host Response to SARS-CoV-2 Pulmonary Infection.

The relationship of SARS-CoV-2 lung infection and severity of pulmonary disease is not fully understood. We analyzed autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter- and intra- patient heterogeneity of pulmonary virus infection. There was a spectrum of high and low virus cases that was associated with duration of disease and activation of interferon pathway genes. Using a digital spatial profiling platform, the virus corresponded to distinct spatial expression of interferon response genes and immune checkpoint genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Epithelial to mesenchymal plasticity and differential response to therapies in pancreatic ductal adenocarcinoma.

Transcriptional profiling has defined pancreatic ductal adenocarcinoma (PDAC) into distinct subtypes with the majority being classical epithelial (E) or quasi-mesenchymal (QM). Despite clear differences in clinical behavior, growing evidence indicates these subtypes exist on a continuum with features of both subtypes present and suggestive of interconverting cell states. Here, we investigated the impact of different therapies being evaluated in PDAC on the phenotypic spectrum of the E/QM state. We demonstrate using RNA-sequencing and RNA-in situ hybridization (RNA-ISH) that FOLFIRINOX combination chemotherapy induces a common shift of both E and QM PDAC toward a more QM state in cell lines and patient tumors. In contrast, Vitamin D, another drug under clinical investigation in PDAC, induces distinct transcriptional responses in each PDAC subtype, with augmentation of the baseline E and QM state. Importantly, this translates to functional changes that increase metastatic propensity in QM PDAC, but decrease dissemination in E PDAC in vivo models. These data exemplify the importance of both the initial E/QM subtype and the plasticity of E/QM states in PDAC in influencing response to therapy, which highlights their relevance in guiding clinical trials.

Stromal Microenvironment Shapes the Intratumoral Architecture of Pancreatic Cancer.

Single-cell technologies have described heterogeneity across tissues, but the spatial distribution and forces that drive single-cell phenotypes have not been well defined. Combining single-cell RNA and protein analytics in studying the role of stromal cancer-associated fibroblasts (CAFs) in modulating heterogeneity in pancreatic cancer (pancreatic ductal adenocarcinoma [PDAC]) model systems, we have identified significant single-cell population shifts toward invasive epithelial-to-mesenchymal transition (EMT) and proliferative (PRO) phenotypes linked with mitogen-activated protein kinase (MAPK) and signal transducer and activator of transcription 3 (STAT3) signaling. Using high-content digital imaging of RNA in situ hybridization in 195 PDAC tumors, we quantified these EMT and PRO subpopulations in 319,626 individual cancer cells that can be classified within the context of distinct tumor gland "units." Tumor gland typing provided an additional layer of intratumoral heterogeneity that was associated with differences in stromal abundance and clinical outcomes. This demonstrates the impact of the stroma in shaping tumor architecture by altering inherent patterns of tumor glands in human PDAC.

Engineered nanointerfaces for microfluidic isolation and molecular profiling of tumor-specific extracellular vesicles.

Extracellular vesicles (EVs) carry RNA, DNA, proteins, and lipids. Specifically, tumor-derived EVs have the potential to be utilized as disease-specific biomarkers. However, a lack of methods to isolate tumor-specific EVs has limited their use in clinical settings. Here we report a sensitive analytical microfluidic platform (EVHB-Chip) that enables tumor-specific EV-RNA isolation within 3 h. Using the EVHB-Chip, we achieve 94% tumor-EV specificity, a limit of detection of 100 EVs per µL, and a 10-fold increase in tumor RNA enrichment in comparison to other methods. Our approach allows for the subsequent release of captured tumor EVs, enabling downstream characterization and functional studies. Processing serum and plasma samples from glioblastoma multiforme (GBM) patients, we can detect the mutant EGFRvIII mRNA. Moreover, using next-generation RNA sequencing, we identify genes specific to GBM as well as transcripts that are hallmarks for the four genetic subtypes of the disease.

Whole blood stabilization for the microfluidic isolation and molecular characterization of circulating tumor cells.

Precise rare-cell technologies require the blood to be processed immediately or be stabilized with fixatives. Such restrictions limit the translation of circulating tumor cell (CTC)-based liquid biopsy assays that provide accurate molecular data in guiding clinical decisions. Here we describe a method to preserve whole blood in its minimally altered state by combining hypothermic preservation with targeted strategies that counter cooling-induced platelet activation. Using this method, whole blood preserved for up to 72 h can be readily processed for microfluidic sorting without compromising CTC yield and viability. The tumor cells retain high-quality intact RNA suitable for single-cell RT-qPCR as well as RNA-Seq, enabling the reliable detection of cancer-specific transcripts including the androgen-receptor splice variant 7 in a cohort of prostate cancer patients with an overall concordance of 92% between fresh and preserved blood. This work will serve as a springboard for the dissemination of diverse blood-based diagnostics.

Prediction of potent shRNAs with a sequential classification algorithm.

We present SplashRNA, a sequential classifier to predict potent microRNA-based short hairpin RNAs (shRNAs). Trained on published and novel data sets, SplashRNA outperforms previous algorithms and reliably predicts the most efficient shRNAs for a given gene. Combined with an optimized miR-E backbone, >90% of high-scoring SplashRNA predictions trigger >85% protein knockdown when expressed from a single genomic integration. SplashRNA can significantly improve the accuracy of loss-of-function genetics studies and facilitates the generation of compact shRNA libraries.

Enhanced Isolation and Release of Circulating Tumor Cells Using Nanoparticle Binding and Ligand Exchange in a Microfluidic Chip.

The detection of rare circulating tumor cells (CTCs) in the blood of cancer patients has the potential to be a powerful and noninvasive method for examining metastasis, evaluating prognosis, assessing tumor sensitivity to drugs, and monitoring therapeutic outcomes. In this study, we have developed an efficient strategy to isolate CTCs from the blood of breast cancer patients using a microfluidic immune-affinity approach. Additionally, to gain further access to these rare cells for downstream characterization, our strategy allows for easy detachment of the captured CTCs from the substrate without compromising cell viability or the ability to employ next generation RNA sequencing for the identification of specific breast cancer genes. To achieve this, a chemical ligand-exchange reaction was engineered to release cells attached to a gold nanoparticle coating bound to the surface of a herringbone microfluidic chip (NP-HBCTC-Chip). Compared to the use of the unmodified HBCTC-Chip, our approach provides several advantages, including enhanced capture efficiency and recovery of isolated CTCs.

The Functional Human C-Terminome.

All translated proteins end with a carboxylic acid commonly called the C-terminus. Many short functional sequences (minimotifs) are located on or immediately proximal to the C-terminus. However, information about the function of protein C-termini has not been consolidated into a single source. Here, we built a new "C-terminome" database and web system focused on human proteins. Approximately 3,600 C-termini in the human proteome have a minimotif with an established molecular function. To help evaluate the function of the remaining C-termini in the human proteome, we inferred minimotifs identified by experimentation in rodent cells, predicted minimotifs based upon consensus sequence matches, and predicted novel highly repetitive sequences in C-termini. Predictions can be ranked by enrichment scores or Gene Evolutionary Rate Profiling (GERP) scores, a measurement of evolutionary constraint. By searching for new anchored sequences on the last 10 amino acids of proteins in the human proteome with lengths between 3-10 residues and up to 5 degenerate positions in the consensus sequences, we have identified new consensus sequences that predict instances in the majority of human genes. All of this information is consolidated into a database that can be accessed through a C-terminome web system with search and browse functions for minimotifs and human proteins. A known consensus sequence-based predicted function is assigned to nearly half the proteins in the human proteome. Weblink: http://cterminome.bio-toolkit.com.

BRD4 Connects Enhancer Remodeling to Senescence Immune Surveillance.

UNLABELLED: Oncogene-induced senescence is a potent barrier to tumorigenesis that limits cellular expansion following certain oncogenic events. Senescent cells display a repressive chromatin configuration thought to stably silence proliferation-promoting genes while simultaneously activating an unusual form of immune surveillance involving a secretory program referred to as the senescence-associated secretory phenotype (SASP). Here, we demonstrate that senescence also involves a global remodeling of the enhancer landscape with recruitment of the chromatin reader BRD4 to newly activated super-enhancers adjacent to key SASP genes. Transcriptional profiling and functional studies indicate that BRD4 is required for the SASP and downstream paracrine signaling. Consequently, BRD4 inhibition disrupts immune cell-mediated targeting and elimination of premalignant senescent cells in vitro and in vivo Our results identify a critical role for BRD4-bound super-enhancers in senescence immune surveillance and in the proper execution of a tumor-suppressive program. SIGNIFICANCE: This study reveals how cells undergoing oncogene-induced senescence acquire a distinctive enhancer landscape that includes formation of super-enhancers adjacent to immune-modulatory genes required for paracrine immune activation. This process links BRD4 and super-enhancers to a tumor-suppressive immune surveillance program that can be disrupted by small molecule inhibitors of the bromo and extra terminal domain family of proteins. Cancer Discov; 6(6); 612-29. ©2016 AACR.See related commentary by Vizioli and Adams, p. 576This article is highlighted in the In This Issue feature, p. 561.

A tunable delivery platform to provide local chemotherapy for pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most devastating and painful cancers. It is often highly resistant to therapy owing to inherent chemoresistance and the desmoplastic response that creates a barrier of fibrous tissue preventing transport of chemotherapeutics into the tumor. The growth of the tumor in pancreatic cancer often leads to invasion of other organs and partial or complete biliary obstruction, inducing intense pain for patients and necessitating tumor resection or repeated stenting. Here, we have developed a delivery device to provide enhanced palliative therapy for pancreatic cancer patients by providing high concentrations of chemotherapeutic compounds locally at the tumor site. This treatment could reduce the need for repeated procedures in advanced PDAC patients to debulk the tumor mass or stent the obstructed bile duct. To facilitate clinical translation, we created the device out of currently approved materials and drugs. We engineered an implantable poly(lactic-co-glycolic)-based biodegradable device that is able to linearly release high doses of chemotherapeutic drugs for up to 60 days. We created five patient-derived PDAC cell lines and tested their sensitivity to approved chemotherapeutic compounds. These in vitro experiments showed that paclitaxel was the most effective single agent across all cell lines. We compared the efficacy of systemic and local paclitaxel therapy on the patient-derived cell lines in an orthotopic xenograft model in mice (PDX). In this model, we found up to a 12-fold increase in suppression of tumor growth by local therapy in comparison to systemic administration and reduce retention into off-target organs. Herein, we highlight the efficacy of a local therapeutic approach to overcome PDAC chemoresistance and reduce the need for repeated interventions and biliary obstruction by preventing local tumor growth. Our results underscore the urgent need for an implantable drug-eluting platform to deliver cytotoxic agents directly within the tumor mass as a novel therapeutic strategy for patients with pancreatic cancer.

GUIDE-seq enables genome-wide profiling of off-target cleavage by CRISPR-Cas nucleases.

CRISPR RNA-guided nucleases (RGNs) are widely used genome-editing reagents, but methods to delineate their genome-wide, off-target cleavage activities have been lacking. Here we describe an approach for global detection of DNA double-stranded breaks (DSBs) introduced by RGNs and potentially other nucleases. This method, called genome-wide, unbiased identification of DSBs enabled by sequencing (GUIDE-seq), relies on capture of double-stranded oligodeoxynucleotides into DSBs. Application of GUIDE-seq to 13 RGNs in two human cell lines revealed wide variability in RGN off-target activities and unappreciated characteristics of off-target sequences. The majority of identified sites were not detected by existing computational methods or chromatin immunoprecipitation sequencing (ChIP-seq). GUIDE-seq also identified RGN-independent genomic breakpoint 'hotspots'. Finally, GUIDE-seq revealed that truncated guide RNAs exhibit substantially reduced RGN-induced, off-target DSBs. Our experiments define the most rigorous framework for genome-wide identification of RGN off-target effects to date and provide a method for evaluating the safety of these nucleases before clinical use.

EWS-FLI1 utilizes divergent chromatin remodeling mechanisms to directly activate or repress enhancer elements in Ewing sarcoma.

The aberrant transcription factor EWS-FLI1 drives Ewing sarcoma, but its molecular function is not completely understood. We find that EWS-FLI1 reprograms gene regulatory circuits in Ewing sarcoma by directly inducing or repressing enhancers. At GGAA repeat elements, which lack evolutionary conservation and regulatory potential in other cell types, EWS-FLI1 multimers induce chromatin opening and create de novo enhancers that physically interact with target promoters. Conversely, EWS-FLI1 inactivates conserved enhancers containing canonical ETS motifs by displacing wild-type ETS transcription factors. These divergent chromatin-remodeling patterns repress tumor suppressors and mesenchymal lineage regulators while activating oncogenes and potential therapeutic targets, such as the kinase VRK1. Our findings demonstrate how EWS-FLI1 establishes an oncogenic regulatory program governing both tumor survival and differentiation.

Development of siRNA payloads to target KRAS-mutant cancer.

UNLABELLED: RNAi is a powerful tool for target identification and can lead to novel therapies for pharmacologically intractable targets such as KRAS. RNAi therapy must combine potent siRNA payloads with reliable in vivo delivery for efficient target inhibition. We used a functional "Sensor" assay to establish a library of potent siRNAs against RAS pathway genes and to show that they efficiently suppress their targets at low dose. This reduces off-target effects and enables combination gene knockdown. We administered Sensor siRNAs in vitro and in vivo and validated the delivery of KRAS siRNA alone and siRNA targeting the complete RAF effector node (A/B/CRAF) as promising strategies to treat KRAS-mutant colorectal cancer. We further demonstrate that improved therapeutic efficacy is achieved by formulating siRNA payloads that combine both single-gene siRNA and node-targeted siRNAs (KRAS + PIK3CA/B). The customizable nature of Sensor siRNA payloads offers a universal platform for the combination target identification and development of RNAi therapeutics. SIGNIFICANCE: To advance RNAi therapy for KRAS-mutant cancer, we developed a validated siRNA library against RAS pathway genes that enables combination gene silencing. Using an in vivo model for real-time siRNA delivery tracking, we show that siRNA-mediated inhibition of KRAS as well as RAF or PI3K combinations can impair KRAS-mutant colorectal cancer in xenograft models.

Dimeric CRISPR RNA-guided FokI nucleases for highly specific genome editing.

Monomeric CRISPR-Cas9 nucleases are widely used for targeted genome editing but can induce unwanted off-target mutations with high frequencies. Here we describe dimeric RNA-guided FokI nucleases (RFNs) that can recognize extended sequences and edit endogenous genes with high efficiencies in human cells. RFN cleavage activity depends strictly on the binding of two guide RNAs (gRNAs) to DNA with a defined spacing and orientation substantially reducing the likelihood that a suitable target site will occur more than once in the genome and therefore improving specificities relative to wild-type Cas9 monomers. RFNs guided by a single gRNA generally induce lower levels of unwanted mutations than matched monomeric Cas9 nickases. In addition, we describe a simple method for expressing multiple gRNAs bearing any 5' end nucleotide, which gives dimeric RFNs a broad targeting range. RFNs combine the ease of RNA-based targeting with the specificity enhancement inherent to dimerization and are likely to be useful in applications that require highly precise genome editing.

Non-cell-autonomous tumor suppression by p53.

The p53 tumor suppressor can restrict malignant transformation by triggering cell-autonomous programs of cell-cycle arrest or apoptosis. p53 also promotes cellular senescence, a tumor-suppressive program that involves stable cell-cycle arrest and secretion of factors that modify the tissue microenvironment. In the presence of chronic liver damage, we show that ablation of a p53-dependent senescence program in hepatic stellate cells increases liver fibrosis and cirrhosis associated with reduced survival and enhances the transformation of adjacent epithelial cells into hepatocellular carcinoma. p53-expressing senescent stellate cells release factors that skew macrophage polarization toward a tumor-inhibiting M1-state capable of attacking senescent cells in culture, whereas proliferating p53-deficient stellate cells secrete factors that stimulate polarization of macrophages into a tumor-promoting M2-state and enhance the proliferation of premalignant cells. Hence, p53 can act non-cell autonomously to suppress tumorigenesis by promoting an antitumor microenvironment, in part, through secreted factors that modulate macrophage function.

Secondary structure, a missing component of sequence-based minimotif definitions.

Minimotifs are short contiguous segments of proteins that have a known biological function. The hundreds of thousands of minimotifs discovered thus far are an important part of the theoretical understanding of the specificity of protein-protein interactions, posttranslational modifications, and signal transduction that occur in cells. However, a longstanding problem is that the different abstractions of the sequence definitions do not accurately capture the specificity, despite decades of effort by many labs. We present evidence that structure is an essential component of minimotif specificity, yet is not used in minimotif definitions. Our analysis of several known minimotifs as case studies, analysis of occurrences of minimotifs in structured and disordered regions of proteins, and review of the literature support a new model for minimotif definitions that includes sequence, structure, and function.

Efficient parallel and out of core algorithms for constructing large bi-directed de Bruijn graphs.

BACKGROUND: Assembling genomic sequences from a set of overlapping reads is one of the most fundamental problems in computational biology. Algorithms addressing the assembly problem fall into two broad categories - based on the data structures which they employ. The first class uses an overlap/string graph and the second type uses a de Bruijn graph. However with the recent advances in short read sequencing technology, de Bruijn graph based algorithms seem to play a vital role in practice. Efficient algorithms for building these massive de Bruijn graphs are very essential in large sequencing projects based on short reads. In an earlier work, an O(n/p) time parallel algorithm has been given for this problem. Here n is the size of the input and p is the number of processors. This algorithm enumerates all possible bi-directed edges which can overlap with a node and ends up generating Θ(nΣ) messages (Σ being the size of the alphabet). RESULTS: In this paper we present a Θ(n/p) time parallel algorithm with a communication complexity that is equal to that of parallel sorting and is not sensitive to Σ. The generality of our algorithm makes it very easy to extend it even to the out-of-core model and in this case it has an optimal I/O complexity of Θ(nlog(n/B)Blog(M/B)) (M being the main memory size and B being the size of the disk block). We demonstrate the scalability of our parallel algorithm on a SGI/Altix computer. A comparison of our algorithm with the previous approaches reveals that our algorithm is faster--both asymptotically and practically. We demonstrate the scalability of our sequential out-of-core algorithm by comparing it with the algorithm used by VELVET to build the bi-directed de Bruijn graph. Our experiments reveal that our algorithm can build the graph with a constant amount of memory, which clearly outperforms VELVET. We also provide efficient algorithms for the bi-directed chain compaction problem. CONCLUSIONS: The bi-directed de Bruijn graph is a fundamental data structure for any sequence assembly program based on Eulerian approach. Our algorithms for constructing Bi-directed de Bruijn graphs are efficient in parallel and out of core settings. These algorithms can be used in building large scale bi-directed de Bruijn graphs. Furthermore, our algorithms do not employ any all-to-all communications in a parallel setting and perform better than the prior algorithms. Finally our out-of-core algorithm is extremely memory efficient and can replace the existing graph construction algorithm in VELVET.