A Somatically Acquired Enhancer of the Androgen Receptor Is a Noncoding Driver in Advanced Prostate Cancer.

Increased androgen receptor (AR) activity drives therapeutic resistance in advanced prostate cancer. The most common resistance mechanism is amplification of this locus presumably targeting the AR gene. Here, we identify and characterize a somatically acquired AR enhancer located 650 kb centromeric to the AR. Systematic perturbation of this enhancer using genome editing decreased proliferation by suppressing AR levels. Insertion of an additional copy of this region sufficed to increase proliferation under low androgen conditions and to decrease sensitivity to enzalutamide. Epigenetic data generated in localized prostate tumors and benign specimens support the notion that this region is a developmental enhancer. Collectively, these observations underscore the importance of epigenomic profiling in primary specimens and the value of deploying genome editing to functionally characterize noncoding elements. More broadly, this work identifies a therapeutic vulnerability for targeting the AR and emphasizes the importance of regulatory elements as highly recurrent oncogenic drivers.

Targetable leukaemia dependency on noncanonical PI3Kγ signalling.

Phosphoinositide-3-kinase-γ (PI3Kγ) is implicated as a target to repolarize tumour-associated macrophages and promote antitumour immune responses in solid cancers1-4. However, cancer cell-intrinsic roles of PI3Kγ are unclear. Here, by integrating unbiased genome-wide CRISPR interference screening with functional analyses across acute leukaemias, we define a selective dependency on the PI3Kγ complex in a high-risk subset that includes myeloid, lymphoid and dendritic lineages. This dependency is characterized by innate inflammatory signalling and activation of phosphoinositide 3-kinase regulatory subunit 5 (PIK3R5), which encodes a regulatory subunit of PI3Kγ5 and stabilizes the active enzymatic complex. We identify p21 (RAC1)-activated kinase 1 (PAK1) as a noncanonical substrate of PI3Kγ that mediates this cell-intrinsic dependency and find that dephosphorylation of PAK1 by PI3Kγ inhibition impairs mitochondrial oxidative phosphorylation. Treatment with the selective PI3Kγ inhibitor eganelisib is effective in leukaemias with activated PIK3R5. In addition, the combination of eganelisib and cytarabine prolongs survival over either agent alone, even in patient-derived leukaemia xenografts with low baseline PIK3R5 expression, as residual leukaemia cells after cytarabine treatment have elevated G protein-coupled purinergic receptor activity and PAK1 phosphorylation. Together, our study reveals a targetable dependency on PI3Kγ-PAK1 signalling that is amenable to near-term evaluation in patients with acute leukaemia.

HIF1A employs CDK8-mediator to stimulate RNAPII elongation in response to hypoxia.

The transcription factor HIF1A is a key mediator of the cellular response to hypoxia. Despite the importance of HIF1A in homeostasis and various pathologies, little is known about how it regulates RNA polymerase II (RNAPII). We report here that HIF1A employs a specific variant of the Mediator complex to stimulate RNAPII elongation. The Mediator-associated kinase CDK8, but not the paralog CDK19, is required for induction of many HIF1A target genes. HIF1A induces binding of CDK8-Mediator and the super elongation complex (SEC), containing AFF4 and CDK9, to alleviate RNAPII pausing. CDK8 is dispensable for HIF1A chromatin binding and histone acetylation, but it is essential for binding of SEC and RNAPII elongation. Global analysis of active RNAPII reveals that hypoxia-inducible genes are paused and active prior to their induction. Our results provide a mechanistic link between HIF1A and CDK8, two potent oncogenes, in the cellular response to hypoxia.

Cancer-stromal cell interactions in breast cancer brain metastases induce glycocalyx-mediated resistance to HER2-targeting therapies.

Brain metastatic breast cancer is particularly lethal largely due to therapeutic resistance. Almost half of the patients with metastatic HER2-positive breast cancer develop brain metastases, representing a major clinical challenge. We previously described that cancer-associated fibroblasts are an important source of resistance in primary tumors. Here, we report that breast cancer brain metastasis stromal cell interactions in 3D cocultures induce therapeutic resistance to HER2-targeting agents, particularly to the small molecule inhibitor of HER2/EGFR neratinib. We investigated the underlying mechanisms using a synthetic Notch reporter system enabling the sorting of cancer cells that directly interact with stromal cells. We identified mucins and bulky glycoprotein synthesis as top-up-regulated genes and pathways by comparing the gene expression and chromatin profiles of stroma-contact and no-contact cancer cells before and after neratinib treatment. Glycoprotein gene signatures were also enriched in human brain metastases compared to primary tumors. We confirmed increased glycocalyx surrounding cocultures by immunofluorescence and showed that mucinase treatment increased sensitivity to neratinib by enabling a more efficient inhibition of EGFR/HER2 signaling in cancer cells. Overexpression of truncated MUC1 lacking the intracellular domain as a model of increased glycocalyx-induced resistance to neratinib both in cell culture and in experimental brain metastases in immunodeficient mice. Our results highlight the importance of glycoproteins as a resistance mechanism to HER2-targeting therapies in breast cancer brain metastases.

MIRA: joint regulatory modeling of multimodal expression and chromatin accessibility in single cells.

Rigorously comparing gene expression and chromatin accessibility in the same single cells could illuminate the logic of how coupling or decoupling of these mechanisms regulates fate commitment. Here we present MIRA, probabilistic multimodal models for integrated regulatory analysis, a comprehensive methodology that systematically contrasts transcription and accessibility to infer the regulatory circuitry driving cells along cell state trajectories. MIRA leverages topic modeling of cell states and regulatory potential modeling of individual gene loci. MIRA thereby represents cell states in an efficient and interpretable latent space, infers high-fidelity cell state trees, determines key regulators of fate decisions at branch points and exposes the variable influence of local accessibility on transcription at distinct loci. Applied to epidermal differentiation and embryonic brain development from two different multimodal platforms, MIRA revealed that early developmental genes were tightly regulated by local chromatin landscape whereas terminal fate genes were titrated without requiring extensive chromatin remodeling.

An Enhancer-Driven Stem Cell-Like Program Mediated by SOX9 Blocks Intestinal Differentiation in Colorectal Cancer.

BACKGROUND AND AIMS: Genomic alterations that encourage stem cell activity and hinder proper maturation are central to the development of colorectal cancer (CRC). Key molecular mediators that promote these malignant properties require further elucidation to galvanize translational advances. We therefore aimed to characterize a key factor that blocks intestinal differentiation, define its transcriptional and epigenetic program, and provide preclinical evidence for therapeutic targeting in CRC. METHODS: Intestinal tissue from transgenic mice and patients were analyzed by means of histopathology and immunostaining. Human CRC cells and neoplastic murine organoids were genetically manipulated for functional studies. Gene expression profiling was obtained through RNA sequencing. Histone modifications and transcription factor binding were determined with the use of chromatin immunoprecipitation sequencing. RESULTS: We demonstrate that SRY-box transcription factor 9 (SOX9) promotes CRC by activating a stem cell-like program that hinders intestinal differentiation. Intestinal adenomas and colorectal adenocarcinomas from mouse models and patients demonstrate ectopic and elevated expression of SOX9. Functional experiments indicate a requirement for SOX9 in human CRC cell lines and engineered neoplastic organoids. Disrupting SOX9 activity impairs primary CRC tumor growth by inducing intestinal differentiation. By binding to genome wide enhancers, SOX9 directly activates genes associated with Paneth and stem cell activity, including prominin 1 (PROM1). SOX9 up-regulates PROM1 via a Wnt-responsive intronic enhancer. A pentaspan transmembrane protein, PROM1 uses its first intracellular domain to support stem cell signaling, at least in part through SOX9, reinforcing a PROM1-SOX9 positive feedback loop. CONCLUSIONS: These studies establish SOX9 as a central regulator of an enhancer-driven stem cell-like program and carry important implications for developing therapeutics directed at overcoming differentiation defects in CRC.

Subependymal giant cell astrocytomas are characterized by mTORC1 hyperactivation, a very low somatic mutation rate, and a unique gene expression profile.

Subependymal giant-cell astrocytomas (SEGAs) are slow-growing brain tumors that are a hallmark feature seen in 5-10% of patients with Tuberous Sclerosis Complex (TSC). Though histologically benign, they can cause serious neurologic symptoms, leading to death if untreated. SEGAs consistently show biallelic loss of TSC1 or TSC2. Herein, we aimed to define other somatic events beyond TSC1/TSC2 loss and identify potential transcriptional drivers that contribute to SEGA formation. Paired tumor-normal whole-exome sequencing was performed on 21 resected SEGAs from 20 TSC patients. Pathogenic variants in TSC1/TSC2 were identified in 19/21 (90%) SEGAs. Copy neutral loss of heterozygosity (size range: 2.2-46 Mb) was seen in 76% (16/21) of SEGAs (44% chr9q and 56% chr16p). An average of 1.4 other somatic variants (range 0-7) per tumor were identified, unlikely of pathogenic significance. Whole transcriptome RNA-sequencing analyses revealed 190 common differentially expressed genes in SEGA (n = 16, 13 from a prior study) in pairwise comparison to each of: low grade diffuse gliomas (n = 530) and glioblastoma (n = 171) from The Cancer Genome Atlas (TCGA) consortium, ganglioglioma (n = 10), TSC cortical tubers (n = 15), and multiple normal tissues. Among these, homeobox transcription factors (TFs) HMX3, HMX2, VAX1, SIX3; and TFs IRF6 and EOMES were all expressed >12-fold higher in SEGAs (FDR/q-value < 0.05). Immunohistochemistry supported the specificity of IRF6, VAX1, SIX3 for SEGAs in comparison to other tumor entities and normal brain. We conclude that SEGAs have an extremely low somatic mutation rate, suggesting that TSC1/TSC2 loss is sufficient to drive tumor growth. The unique and highly expressed SEGA-specific TFs likely reflect the neuroepithelial cell of origin, and may also contribute to the transcriptional and epigenetic state that enables SEGA growth following two-hit loss of TSC1 or TSC2 and mTORC1 activation.

Embryonic transcription factor SOX9 drives breast cancer endocrine resistance.

The estrogen receptor (ER) drives the growth of most luminal breast cancers and is the primary target of endocrine therapy. Although ER blockade with drugs such as tamoxifen is very effective, a major clinical limitation is the development of endocrine resistance especially in the setting of metastatic disease. Preclinical and clinical observations suggest that even following the development of endocrine resistance, ER signaling continues to exert a pivotal role in tumor progression in the majority of cases. Through the analysis of the ER cistrome in tamoxifen-resistant breast cancer cells, we have uncovered a role for an RUNX2-ER complex that stimulates the transcription of a set of genes, including most notably the stem cell factor SOX9, that promote proliferation and a metastatic phenotype. We show that up-regulation of SOX9 is sufficient to cause relative endocrine resistance. The gain of SOX9 as an ER-regulated gene associated with tamoxifen resistance was validated in a unique set of clinical samples supporting the need for the development of improved ER antagonists.

VIPER: Visualization Pipeline for RNA-seq, a Snakemake workflow for efficient and complete RNA-seq analysis.

BACKGROUND: RNA sequencing has become a ubiquitous technology used throughout life sciences as an effective method of measuring RNA abundance quantitatively in tissues and cells. The increase in use of RNA-seq technology has led to the continuous development of new tools for every step of analysis from alignment to downstream pathway analysis. However, effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts. RESULTS: Using the workflow management system Snakemake we have developed a user friendly, fast, efficient, and comprehensive pipeline for RNA-seq analysis. VIPER (Visualization Pipeline for RNA-seq analysis) is an analysis workflow that combines some of the most popular tools to take RNA-seq analysis from raw sequencing data, through alignment and quality control, into downstream differential expression and pathway analysis. VIPER has been created in a modular fashion to allow for the rapid incorporation of new tools to expand the capabilities. This capacity has already been exploited to include very recently developed tools that explore immune infiltrate and T-cell CDR (Complementarity-Determining Regions) reconstruction abilities. The pipeline has been conveniently packaged such that minimal computational skills are required to download and install the dozens of software packages that VIPER uses. CONCLUSIONS: VIPER is a comprehensive solution that performs most standard RNA-seq analyses quickly and effectively with a built-in capacity for customization and expansion.

ChiLin: a comprehensive ChIP-seq and DNase-seq quality control and analysis pipeline.

BACKGROUND: Transcription factor binding, histone modification, and chromatin accessibility studies are important approaches to understanding the biology of gene regulation. ChIP-seq and DNase-seq have become the standard techniques for studying protein-DNA interactions and chromatin accessibility respectively, and comprehensive quality control (QC) and analysis tools are critical to extracting the most value from these assay types. Although many analysis and QC tools have been reported, few combine ChIP-seq and DNase-seq data analysis and quality control in a unified framework with a comprehensive and unbiased reference of data quality metrics. RESULTS: ChiLin is a computational pipeline that automates the quality control and data analyses of ChIP-seq and DNase-seq data. It is developed using a flexible and modular software framework that can be easily extended and modified. ChiLin is ideal for batch processing of many datasets and is well suited for large collaborative projects involving ChIP-seq and DNase-seq from different designs. ChiLin generates comprehensive quality control reports that include comparisons with historical data derived from over 23,677 public ChIP-seq and DNase-seq samples (11,265 datasets) from eight literature-based classified categories. To the best of our knowledge, this atlas represents the most comprehensive ChIP-seq and DNase-seq related quality metric resource currently available. These historical metrics provide useful heuristic quality references for experiment across all commonly used assay types. Using representative datasets, we demonstrate the versatility of the pipeline by applying it to different assay types of ChIP-seq data. The pipeline software is available open source at https://github.com/cfce/chilin . CONCLUSION: ChiLin is a scalable and powerful tool to process large batches of ChIP-seq and DNase-seq datasets. The analysis output and quality metrics have been structured into user-friendly directories and reports. We have successfully compiled 23,677 profiles into a comprehensive quality atlas with fine classification for users.

CistromeFinder for ChIP-seq and DNase-seq data reuse.

SUMMARY: Chromatin immunoprecipitation and DNase I hypersensitivity assays with high-throughput sequencing have greatly accelerated the understanding of transcriptional and epigenetic regulation, although data reuse for the community of experimental biologists has been challenging. We created a data portal CistromeFinder that can help query, evaluate and visualize publicly available Chromatin immunoprecipitation and DNase I hypersensitivity assays with high-throughput sequencing data in human and mouse. The database currently contains 6378 samples over 4391 datasets, 313 factors and 102 cell lines or cell populations. Each dataset has gone through a consistent analysis and quality control pipeline; therefore, users could evaluate the overall quality of each dataset before examining binding sites near their genes of interest. CistromeFinder is integrated with UCSC genome browser for visualization, Primer3Plus for ChIP-qPCR primer design and CistromeMap for submitting newly available datasets. It also allows users to leave comments to facilitate data evaluation and update. AVAILABILITY: http://cistrome.org/finder. CONTACT: xsliu@jimmy.harvard.edu or henry_long@dfci.harvard.edu.

scATAnno: Automated Cell Type Annotation for single-cell ATAC Sequencing Data.

Recent advances in single-cell epigenomic techniques have created a growing demand for scATAC-seq analysis. One key analysis task is to determine cell type identity based on the epigenetic data. We introduce scATAnno, a python package designed to automatically annotate scATAC-seq data using large-scale scATAC-seq reference atlases. This workflow generates the reference atlases from publicly available datasets enabling accurate cell type annotation by integrating query data with reference atlases, without the use of scRNA-seq data. To enhance annotation accuracy, we have incorporated KNN-based and weighted distance-based uncertainty scores to effectively detect cell populations within the query data that are distinct from all cell types in the reference data. We compare and benchmark scATAnno against 7 other published approaches for cell annotation and show superior performance in multiple data sets and metrics. We showcase the utility of scATAnno across multiple datasets, including peripheral blood mononuclear cell (PBMC), Triple Negative Breast Cancer (TNBC), and basal cell carcinoma (BCC), and demonstrate that scATAnno accurately annotates cell types across conditions. Overall, scATAnno is a useful tool for scATAC-seq reference building and cell type annotation in scATAC-seq data and can aid in the interpretation of new scATAC-seq datasets in complex biological systems.

Loss of Kmt2c or Kmt2d drives brain metastasis via KDM6A-dependent upregulation of MMP3.

KMT2C and KMT2D, encoding histone H3 lysine 4 methyltransferases, are among the most commonly mutated genes in triple-negative breast cancer (TNBC). However, how these mutations may shape epigenomic and transcriptomic landscapes to promote tumorigenesis is largely unknown. Here we describe that deletion of Kmt2c or Kmt2d in non-metastatic murine models of TNBC drives metastasis, especially to the brain. Global chromatin profiling and chromatin immunoprecipitation followed by sequencing revealed altered H3K4me1, H3K27ac and H3K27me3 chromatin marks in knockout cells and demonstrated enhanced binding of the H3K27me3 lysine demethylase KDM6A, which significantly correlated with gene expression. We identified Mmp3 as being commonly upregulated via epigenetic mechanisms in both knockout models. Consistent with these findings, samples from patients with KMT2C-mutant TNBC have higher MMP3 levels. Downregulation or pharmacological inhibition of KDM6A diminished Mmp3 upregulation induced by the loss of histone-lysine N-methyltransferase 2 (KMT2) and prevented brain metastasis similar to direct downregulation of Mmp3. Taken together, we identified the KDM6A-matrix metalloproteinase 3 axis as a key mediator of KMT2C/D loss-driven metastasis in TNBC.

Transcription factor dynamics, oscillation, and functions in human enteroendocrine cell differentiation.

Enteroendocrine cells (EECs) secrete serotonin (enterochromaffin [EC] cells) or specific peptide hormones (non-EC cells) that serve vital metabolic functions. The basis for terminal EEC diversity remains obscure. By forcing activity of the transcription factor (TF) NEUROG3 in 2D cultures of human intestinal stem cells, we replicated physiologic EEC differentiation and examined transcriptional and cis-regulatory dynamics that culminate in discrete cell types. Abundant EEC precursors expressed stage-specific genes and TFs. Before expressing pre-terminal NEUROD1, post-mitotic precursors oscillated between transcriptionally distinct ASCL1+ and HES6hi cell states. Loss of either factor accelerated EEC differentiation substantially and disrupted EEC individuality; ASCL1 or NEUROD1 deficiency had opposing consequences on EC and non-EC cell features. These TFs mainly bind cis-elements that are accessible in undifferentiated stem cells, and they tailor subsequent expression of TF combinations that underlie discrete EEC identities. Thus, early TF oscillations retard EEC maturation to enable accurate diversity within a medically important cell lineage.

Transcription factor dynamics, oscillation, and functions in human enteroendocrine cell differentiation.

Enteroendocrine cells (EECs), which secrete serotonin (enterochromaffin cells, EC) or a dominant peptide hormone, serve vital physiologic functions. As with any adult human lineage, the basis for terminal cell diversity remains obscure. We replicated human EEC differentiation in vitro , mapped transcriptional and chromatin dynamics that culminate in discrete cell types, and studied abundant EEC precursors expressing selected transcription factors (TFs) and gene programs. Before expressing the pre-terminal factor NEUROD1, non-replicating precursors oscillated between epigenetically similar but transcriptionally distinct ASCL1 + and HES6 hi cell states. Loss of either factor substantially accelerated EEC differentiation and disrupted EEC individuality; ASCL1 or NEUROD1 deficiency had opposing consequences on EC and hormone-producing cell features. Expressed late in EEC differentiation, the latter TFs mainly bind cis -elements that are accessible in undifferentiated stem cells and tailor the subsequent expression of TF combinations that specify EEC types. Thus, TF oscillations retard EEC maturation to enable accurate EEC diversification.

Increased chromatin accessibility drives transition to androgen receptor splice variant dependence in castration-resistant prostate cancer.

Androgen receptor (AR) splice variants, of which ARv7 is the most common, are increased in prostate cancer (PC) that develops resistance to androgen signaling inhibitor drugs, but the extent to which these variants drive AR activity, and whether they have novel functions or dependencies, remain to be determined. We generated a subline of VCaP PC cells (VCaP16) that is resistant to the AR inhibitor enzalutamide (ENZ) and found that AR activity was independent of the full-length AR (ARfl), despite its continued high-level expression, and was instead driven by ARv7. The ARv7 cistrome and transcriptome in VCaP16 cells mirrored that of the ARfl in VCaP cells, although ARv7 chromatin binding was weaker, and strong ARv7 binding sites correlated with higher affinity ARfl binding sites across multiple models and clinical samples. Notably, although ARv7 expression in VCaP cells increased rapidly in response to ENZ, there was a long lag before it gained chromatin binding and transcriptional activity. This lag was associated with an increase in chromatin accessibility, with the AR and nuclear factor I (NFI) motifs being most enriched at these more accessible sites. Moreover, the transcriptional effects of combined NFIB and NFIX knockdown versus ARv7 knockdown were highly correlated. These findings indicate that ARv7 can drive the AR program, but that its activity is dependent on adaptations that increase chromatin accessibility to enhance its intrinsically weak chromatin binding.

Mammalian SWI/SNF complex activity regulates POU2F3 and constitutes a targetable dependency in small cell lung cancer.

Small cell lung cancers (SCLC) are comprised of heterogeneous subtypes marked by lineage-specific transcription factors, including ASCL1, NEUROD1, and POU2F3. POU2F3-positive SCLC, ∼12% of all cases, are uniquely dependent on POU2F3 itself; as such, approaches to attenuate POU2F3 expression may represent new therapeutic opportunities. Here using genome-scale screens for regulators of POU2F3 expression and SCLC proliferation, we define mSWI/SNF complexes, including non-canonical BAF (ncBAF) complexes, as top dependencies specific to POU2F3-positive SCLC. Notably, clinical-grade pharmacologic mSWI/SNF inhibition attenuates proliferation of all POU2F3-positive SCLCs, while disruption of ncBAF via BRD9 degradation is uniquely effective in pure non-neuroendocrine POU2F3-SCLCs. mSWI/SNF maintains accessibility over gene loci central to POU2F3-mediated gene regulatory networks. Finally, chemical targeting of SMARCA4/2 mSWI/SNF ATPases and BRD9 decrease POU2F3-SCLC tumor growth and increase survival in vivo . Taken together, these results characterize mSWI/SNF-mediated global governance of the POU2F3 oncogenic program and suggest mSWI/SNF inhibition as a therapeutic strategy for SCLC.

Thio-2 Inhibits Key Signaling Pathways Required for the Development and Progression of Castration-resistant Prostate Cancer.

Therapies that abrogate persistent androgen receptor (AR) signaling in castration-resistant prostate cancer (CRPC) remain an unmet clinical need. The N-terminal domain of the AR that drives transcriptional activity in CRPC remains a challenging therapeutic target. Herein we demonstrate that BCL-2-associated athanogene-1 (BAG-1) mRNA is highly expressed and associates with signaling pathways, including AR signaling, that are implicated in the development and progression of CRPC. In addition, interrogation of geometric and physiochemical properties of the BAG domain of BAG-1 isoforms identifies it to be a tractable but challenging drug target. Furthermore, through BAG-1 isoform mouse knockout studies, we confirm that BAG-1 isoforms regulate hormone physiology and that therapies targeting the BAG domain will be associated with limited "on-target" toxicity. Importantly, the postulated inhibitor of BAG-1 isoforms, Thio-2, suppressed AR signaling and other important pathways implicated in the development and progression of CRPC to reduce the growth of treatment-resistant prostate cancer cell lines and patient-derived models. However, the mechanism by which Thio-2 elicits the observed phenotype needs further elucidation as the genomic abrogation of BAG-1 isoforms was unable to recapitulate the Thio-2-mediated phenotype. Overall, these data support the interrogation of related compounds with improved drug-like properties as a novel therapeutic approach in CRPC, and further highlight the clinical potential of treatments that block persistent AR signaling which are currently undergoing clinical evaluation in CRPC.

Universal NicE-Seq: A Simple and Quick Method for Accessible Chromatin Detection in Fixed Cells.

Genome-wide accessible chromatin sequencing and identification has enabled deciphering the epigenetic information encoded in chromatin, revealing accessible promoters, enhancers, nucleosome positioning, transcription factor occupancy, and other chromosomal protein binding. The starting biological materials are often fixed using formaldehyde crosslinking. Here, we describe accessible chromatin library preparation from low numbers of formaldehyde-crosslinked cells using a modified nick translation method, where a nicking enzyme nicks one strand of DNA and DNA polymerase incorporates biotin-conjugated dATP, dCTP, and methyl-dCTP. Once the DNA is labeled, it can be isolated for NGS library preparation. We termed this method as universal NicE-seq (nicking enzyme-assisted sequencing). We also demonstrate a single tube method that enables direct NGS library preparation from low cell numbers without DNA purification. Furthermore, we demonstrated universal NicE-seq on FFPE tissue section sample.

Tutorial: integrative computational analysis of bulk RNA-sequencing data to characterize tumor immunity using RIMA.

RNA-sequencing (RNA-seq) has become an increasingly cost-effective technique for molecular profiling and immune characterization of tumors. In the past decade, many computational tools have been developed to characterize tumor immunity from gene expression data. However, the analysis of large-scale RNA-seq data requires bioinformatics proficiency, large computational resources and cancer genomics and immunology knowledge. In this tutorial, we provide an overview of computational analysis of bulk RNA-seq data for immune characterization of tumors and introduce commonly used computational tools with relevance to cancer immunology and immunotherapy. These tools have diverse functions such as evaluation of expression signatures, estimation of immune infiltration, inference of the immune repertoire, prediction of immunotherapy response, neoantigen detection and microbiome quantification. We describe the RNA-seq IMmune Analysis (RIMA) pipeline integrating many of these tools to streamline RNA-seq analysis. We also developed a comprehensive and user-friendly guide in the form of a GitBook with text and video demos to assist users in analyzing bulk RNA-seq data for immune characterization at both individual sample and cohort levels by using RIMA.