Integrative analysis of ultra-deep RNA-seq reveals alternative promoter usage as a mechanism of activating oncogenic programmes during prostate cancer progression.

Transcription factor (TF) proteins regulate gene activity by binding to regulatory regions, most importantly at gene promoters. Many genes have alternative promoters (APs) bound by distinct TFs. The role of differential TF activity at APs during tumour development is poorly understood. Here we show, using deep RNA sequencing in 274 biopsies of benign prostate tissue, localized prostate tumours and metastatic castration-resistant prostate cancer, that AP usage increases as tumours progress and APs are responsible for a disproportionate amount of tumour transcriptional activity. Expression of the androgen receptor (AR), the key driver of prostate tumour activity, is correlated with elevated AP usage. We identified AR, FOXA1 and MYC as potential drivers of AP activation. DNA methylation is a likely mechanism for AP activation during tumour progression and lineage plasticity. Our data suggest that prostate tumours activate APs to magnify the transcriptional impact of tumour drivers, including AR and MYC.

Single cell-transcriptomic analysis informs the lncRNA landscape in metastatic castration resistant prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal form of prostate cancer. Although long-noncoding RNAs (lncRNAs) have been implicated in mCRPC, past studies have relied on bulk sequencing methods with low depth and lack of single-cell resolution. Hence, we performed a lncRNA-focused analysis of single-cell RNA-sequencing data (n = 14) from mCRPC biopsies followed by integration with bulk multi-omic datasets. This yielded 389 cell-enriched lncRNAs in prostate cancer cells and the tumor microenvironment (TME). These lncRNAs demonstrated enrichment with regulatory elements and exhibited alterations during prostate cancer progression. Prostate-lncRNAs were correlated with AR mutational status and response to treatment with enzalutamide, while TME-lncRNAs were associated with RB1 deletions and poor prognosis. Finally, lncRNAs identified between prostate adenocarcinomas and neuroendocrine tumors exhibited distinct expression and methylation profiles. Our findings demonstrate the ability of single-cell analysis to refine our understanding of lncRNAs in mCRPC and serve as a resource for future mechanistic studies.

Genomic and epigenomic analysis of plasma cell-free DNA identifies stemness features associated with worse survival in AR -altered lethal prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) resistant to androgen receptor (AR)-targeted agents is often lethal. Unfortunately, biomarkers for this deadly disease remain under investigation, and underpinning mechanisms are ill-understood. Here, we applied deep sequencing to ∼100 mCRPC patients prior to the initiation of first-line AR-targeted therapy, which detected AR /enhancer alterations in over a third of patients, which correlated with lethality. To delve into the mechanism underlying why these patients with cell-free AR /enhancer alterations developed more lethal prostate cancer, we next performed genome-wide cell-free DNA epigenomics. Strikingly, we found that binding sites for transcription factors associated with developmental stemness were nucleosomally more accessible. These results were corroborated using cell-free DNA methylation data, as well as tumor RNA sequencing from a held-out cohort of mCRPC patients. Thus, we validated the importance of AR /enhancer alterations as a prognostic biomarker in lethal mCRPC, and showed that the underlying mechanism for lethality involves reprogramming developmental states toward increased stemness.

Single-cell transcriptomics reveals long noncoding RNAs associated with tumor biology and the microenvironment in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is highly heterogeneous and lethal. Long noncoding RNAs (lncRNAs) are an important class of genes regulating tumorigenesis and progression. Prior bulk transcriptomic studies in PDAC have revealed the dysregulation of lncRNAs but lack single-cell resolution to distinguish lncRNAs in tumor-intrinsic biology and the tumor microenvironment (TME). We analyzed single-cell transcriptome data from 73 multiregion samples in 21 PDAC patients to evaluate lncRNAs associated with intratumoral heterogeneity and the TME in PDAC. We found 111 cell-specific lncRNAs that reflected tumor, immune and stromal cell contributions, associated with outcomes, and validated across orthogonal datasets. Single-cell analysis of tumor cells revealed lncRNAs associated with TP53 mutations and FOLFIRINOX treatment that were obscured in bulk tumor analysis. Lastly, tumor subcluster analysis revealed widespread intratumor heterogeneity and intratumoral lncRNAs associated with cancer hallmarks and tumor processes such as angiogenesis, epithelial-mesenchymal transition, metabolism and immune signaling. Intratumoral subclusters and lncRNAs were validated across six datasets and showed clinically relevant associations with patient outcomes. Our study provides the first comprehensive assessment of the lncRNA landscape in PDAC using single-cell transcriptomic data and can serve as a resource, PDACLncDB (accessible at https://www.maherlab.com/pdaclncdb-overview), to guide future functional studies.

PACT: a pipeline for analysis of circulating tumor DNA.

MOTIVATION: Detection of genomic alterations in circulating tumor DNA (ctDNA) is currently used for active clinical monitoring of cancer progression and treatment response. While methods for analysis of small mutations are more developed, strategies for detecting structural variants (SVs) in ctDNA are limited. Additionally, reproducibly calling small-scale mutations, copy number alterations, and SVs in ctDNA is challenging due to the lack to unified tools for these different classes of variants. RESULTS: We developed a unified pipeline for the analysis of ctDNA [Pipeline for the Analysis of ctDNA (PACT)] that accurately detects SVs and consistently outperformed similar tools when applied to simulated, cell line, and clinical data. We provide PACT in the form of a Common Workflow Language pipeline which can be run by popular workflow management systems in high-performance computing environments. AVAILABILITY AND IMPLEMENTATION: PACT is freely available at https://github.com/ChrisMaherLab/PACT.

Chemerin triggers migration of a CD8 T cell subset with natural killer cell functions.

The recruitment of cells with effector functions into the tumor microenvironment holds potential for delaying cancer progression. We show that subsets of human CD28-effector CD8 T cells, CCR7- CD45RO+ effector memory, and CCR7- CD45RO- effector memory RA phenotypes, express the chemerin receptor CMKLR1 and bind chemerin via the receptor. CMKLR1-expressing human CD8 effector memory T cells present gene, protein, and cytotoxic features of NK cells. Active chemerin promotes chemotaxis of CMKLR1-expressing CD8 effector memory cells and triggers activation of the α4ß1 integrin. In an experimental prostate tumor mouse model, chemerin expression is downregulated in the tumor microenvironment, which is associated with few tumor-infiltrating CD8+ T cells, while forced overexpression of chemerin by mouse prostate cancer cells leads to an accumulation of intra-tumor CD8+ T cells. Furthermore, α4 integrin blockade abrogated the chemerin-dependent recruitment of CD8+ T effector memory cells into implanted prostate tumors in vivo. The results identify a role for chemerin:CMKLR1 in defining a specialized NK-like CD8 T cell, and suggest the use of chemerin-dependent modalities to target effector CMKLR1-expressing T cells to the tumor microenvironment for immunotherapeutic purposes.

Characterization of cell-type specific circular RNAs associated with colorectal cancer metastasis.

Colorectal cancer (CRC) is the most common gastrointestinal malignancy and a leading cause of cancer deaths in the United States. More than half of CRC patients develop metastatic disease (mCRC) with an average 5-year survival rate of 13%. Circular RNAs (circRNAs) have recently emerged as important tumorigenesis regulators; however, their role in mCRC progression remains poorly characterized. Further, little is known about their cell-type specificity to elucidate their functions in the tumor microenvironment (TME). To address this, we performed total RNA sequencing (RNA-seq) on 30 matched normal, primary and metastatic samples from 14 mCRC patients. Additionally, five CRC cell lines were sequenced to construct a circRNA catalog in CRC. We detected 47 869 circRNAs, with 51% previously unannotated in CRC and 14% novel candidates when compared to existing circRNA databases. We identified 362 circRNAs differentially expressed in primary and/or metastatic tissues, termed circular RNAs associated with metastasis (CRAMS). We performed cell-type deconvolution using published single-cell RNA-seq datasets and applied a non-negative least squares statistical model to estimate cell-type specific circRNA expression. This predicted 667 circRNAs as exclusively expressed in a single cell type. Collectively, this serves as a valuable resource, TMECircDB (accessible at https://www.maherlab.com/tmecircdb-overview), for functional characterization of circRNAs in mCRC, specifically in the TME.

LINC00355 regulates p27KIP expression by binding to MENIN to induce proliferation in late-stage relapse breast cancer.

Late-stage relapse (LSR) in patients with breast cancer (BC) occurs more than five years and up to 10 years after initial treatment and has less than 30% 5-year relative survival rate. Long non-coding RNAs (lncRNAs) play important roles in BC yet have not been studied in LSR BC. Here, we identify 1127 lncRNAs differentially expressed in LSR BC via transcriptome sequencing and analysis of 72 early-stage and 24 LSR BC patient tumors. Decreasing expression of the most up-regulated lncRNA, LINC00355, in BC and MCF7 long-term estrogen deprived cell lines decreases cellular invasion and proliferation. Subsequent mechanistic studies show that LINC00355 binds to MENIN and changes occupancy at the CDKN1B promoter to decrease p27Kip. In summary, this is a key study discovering lncRNAs in LSR BC and LINC00355 association with epigenetic regulation and proliferation in BC.

Loss of Long Noncoding RNA NXTAR in Prostate Cancer Augments Androgen Receptor Expression and Enzalutamide Resistance.

Androgen receptor (AR) signaling continues to play a dominant role in all stages of prostate cancer, including castration-resistant prostate cancers (CRPC) that have developed resistance to second generation AR antagonists such as enzalutamide. In this study, we identified a long noncoding RNA (lncRNA), NXTAR (LOC105373241) that is located convergent with the AR gene and is repressed in human prostate tumors and cell lines. NXTAR bound upstream of the AR promoter and promoted EZH2 recruitment, causing significant loss of AR (and AR-V7) expression. Paradoxically, AR bound the NXTAR promoter, and inhibition of AR by the ACK1/TNK2 small molecule inhibitor (R)-9b excluded AR from the NXTAR promoter. The histone acetyltransferase GCN5 bound and deposited H3K14 acetylation marks, enhancing NXTAR expression. Application of an oligonucleotide derived from NXTAR exon 5 (NXTAR-N5) suppressed AR/AR-V7 expression and prostate cancer cell proliferation, indicating the translational relevance of the negative regulation of AR. In addition, pharmacologic restoration of NXTAR using (R)-9b abrogated enzalutamide-resistant prostate xenograft tumor growth. Overall, this study uncovers a positive feedback loop, wherein NXTAR acts as a novel prostate tumor-suppressing lncRNA by inhibiting AR/AR-V7 expression, which in turn upregulates NXTAR levels, compromising enzalutamide-resistant prostate cancer. The restoration of NXTAR could serve as a new therapeutic modality for patients who have acquired resistance to second generation AR antagonists. SIGNIFICANCE: This study identifies NXTAR as a tumor suppressive lncRNA that can epigenetically downregulate AR/AR-V7 expression and provides a therapeutic strategy to reinstate NXTAR expression for treating recurrent CRPC.

Long noncoding RNAs in cancer metastasis.

Metastasis is a major contributor to cancer-associated deaths. It is characterized by a multistep process that occurs through the acquisition of molecular and phenotypic changes enabling cancer cells from a primary tumour to disseminate and colonize at distant organ sites. Over the past decade, the discovery and characterization of long noncoding RNAs (lncRNAs) have revealed the diversity of their regulatory roles, including key contributions throughout the metastatic cascade. Here, we review how lncRNAs promote metastasis by functioning in discrete pro-metastatic steps including the epithelial-mesenchymal transition, invasion and migration and organotrophic colonization, and by influencing the metastatic tumour microenvironment, often by interacting within ribonucleoprotein complexes or directly with other nucleic acid entities. We discuss well-characterized lncRNAs with in vivo phenotypes and highlight mechanistic commonalities such as convergence with the TGFß-ZEB1/ZEB2 axis or the nuclear factor-κB pathway, in addition to lncRNAs with controversial mechanisms and the influence of methodologies on mechanistic interpretation. Furthermore, some lncRNAs can help identify tumours with increased metastatic risk and spur novel therapeutic strategies, with several lncRNAs having shown potential as novel targets for antisense oligonucleotide therapy in animal models. In addition to well-characterized examples of lncRNAs functioning in metastasis, we discuss controversies and ongoing challenges in lncRNA biology. Finally, we present areas for future study for this rapidly evolving field.

Targeted Therapy to ß3 Integrin Reduces Chemoresistance in Breast Cancer Bone Metastases.

Breast cancer bone metastases are common and incurable. Tumoral integrin ß3 (ß3) expression is induced through interaction with the bone microenvironment. Although ß3 is known to promote bone colonization, its functional role during therapy of established bone metastases is not known. We found increased numbers of ß3+ tumor cells in murine bone metastases after docetaxel chemotherapy. ß3+ tumor cells were present in 97% of post-neoadjuvant chemotherapy triple-negative breast cancer patient samples (n = 38). High tumoral ß3 expression was associated with worse outcomes in both pre- and postchemotherapy triple-negative breast cancer groups. Genetic deletion of tumoral ß3 had minimal effect in vitro, but significantly enhanced in vivo docetaxel activity, particularly in the bone. Rescue experiments confirmed that this effect required intact ß3 signaling. Ultrastructural, transcriptomic, and functional analyses revealed an alternative metabolic response to chemotherapy in ß3-expressing cells characterized by enhanced oxygen consumption, reactive oxygen species generation, and protein production. We identified mTORC1 as a candidate for therapeutic targeting of this ß3-mediated, chemotherapy-induced metabolic response. mTORC1 inhibition in combination with docetaxel synergistically attenuated murine bone metastases. Furthermore, micelle nanoparticle delivery of mTORC1 inhibitor to cells expressing activated αvß3 integrins enhanced docetaxel efficacy in bone metastases. Taken together, we show that ß3 integrin induction by the bone microenvironment promotes resistance to chemotherapy through an altered metabolic response that can be defused by combination with αvß3-targeted mTORC1 inhibitor nanotherapy. Our work demonstrates the importance of the metastatic microenvironment when designing treatments and presents new, bone-specific strategies for enhancing chemotherapeutic efficacy.

SV-HotSpot: detection and visualization of hotspots targeted by structural variants associated with gene expression.

Whole genome sequencing (WGS) has enabled the discovery of genomic structural variants (SVs), including those targeting intergenic and intronic non-coding regions that eluded previous exome focused strategies. However, the field currently lacks an automated tool that analyzes SV candidates to identify recurrent SVs and their targeted sites (hotspot regions), visualizes these genomic events within the context of various functional elements, and evaluates their potential effect on gene expression. To address this, we developed SV-HotSpot, an automated tool that integrates SV candidates, copy number alterations, gene expression, and genome annotations (e.g. gene and regulatory elements) to discover, annotate, and visualize recurrent SVs and their targeted hotspot regions that may affect gene expression. We applied SV-HotSpot to WGS and matched transcriptome data from metastatic castration resistant prostate cancer patients and rediscovered recurrent SVs targeting coding and non-coding functional elements known to promote prostate cancer progression and metastasis. SV-HotSpot provides a valuable resource to integrate SVs, gene expression, and genome annotations for discovering biologically relevant SVs altering coding and non-coding genome. SV-HotSpot is available at https://github.com/ChrisMaherLab/SV-HotSpot .

Cell-free DNA alterations in the AR enhancer and locus predict resistance to AR-directed therapy in patients with metastatic prostate cancer.

PURPOSE: Cell-free DNA (cfDNA) and circulating tumor cell (CTC) based liquid biopsies have emerged as potential tools to predict responses to androgen receptor (AR)-directed therapy in metastatic prostate cancer. However, due to complex mechanisms and incomplete understanding of genomic events involved in metastatic prostate cancer resistance, current assays (e.g. CTC AR-V7) demonstrate low sensitivity and remain underutilized. The recent discovery of AR enhancer amplification in >80% of metastatic patients and its association with disease resistance presents an opportunity to improve upon current assays. We hypothesized that tracking AR/enhancer genomic alterations in plasma cfDNA would detect resistance with high sensitivity and specificity. METHODS: We developed a targeted sequencing and analysis method as part of a new assay called Enhancer and neighboring loci of Androgen Receptor Sequencing (EnhanceAR-Seq). We applied EnhanceAR-Seq to plasma collected from 40 patients with metastatic prostate cancer treated with AR-directed therapy to monitor AR/enhancer genomic alterations and correlate these events with therapy resistance, progression-free survival (PFS) and overall survival (OS). RESULTS: EnhanceAR-Seq identified genomic alterations in the AR/enhancer locus in 45% of cases, including a 40% rate of AR enhancer amplification. Patients with AR/enhancer alterations had significantly worse PFS and OS than those without (6-month PFS: 30% vs. 71%, P=0.0002; 6-month OS: 59% vs. 100%, P=0.0015). AR/enhancer alterations in plasma cfDNA detected 18 of 23 resistant cases (78%) and outperformed the CTC AR-V7 assay which was also run on a subset of patients. CONCLUSION: cfDNA-based AR locus alterations, including of the enhancer, are strongly associated with resistance to AR-directed therapy and significantly worse survival. cfDNA analysis using EnhanceAR-Seq may enable more precise risk stratification and personalized therapeutic approaches for metastatic prostate cancer.

Pan-cancer proteogenomic analysis reveals long and circular noncoding RNAs encoding peptides.

Recent studies show that annotated long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) encode for stable, functional peptides that contribute to human development and disease. To systematically discover lncRNAs and circRNAs encoding peptides, we performed a comprehensive integrative analysis of mass spectrometry-based proteomic and transcriptomic sequencing data from >900 patients across nine cancer types. This enabled us to identify 19,871 novel peptides derived from 8,903 lncRNAs. Further, we exploited open reading frames overlapping the backspliced region of circRNAs to identify 3,238 peptides that are uniquely derived from 2,834 circRNAs and not their corresponding linear RNAs. Collectively, our pan-cancer proteogenomic analysis will serve as a resource for evaluating the coding potential of lncRNAs and circRNAs that could aid future mechanistic studies exploring their function in cancer.

The DNA methylation landscape of advanced prostate cancer.

Although DNA methylation is a key regulator of gene expression, the comprehensive methylation landscape of metastatic cancer has never been defined. Through whole-genome bisulfite sequencing paired with deep whole-genome and transcriptome sequencing of 100 castration-resistant prostate metastases, we discovered alterations affecting driver genes that were detectable only with integrated whole-genome approaches. Notably, we observed that 22% of tumors exhibited a novel epigenomic subtype associated with hypermethylation and somatic mutations in TET2, DNMT3B, IDH1 and BRAF. We also identified intergenic regions where methylation is associated with RNA expression of the oncogenic driver genes AR, MYC and ERG. Finally, we showed that differential methylation during progression preferentially occurs at somatic mutational hotspots and putative regulatory regions. This study is a large integrated study of whole-genome, whole-methylome and whole-transcriptome sequencing in metastatic cancer that provides a comprehensive overview of the important regulatory role of methylation in metastatic castration-resistant prostate cancer.

The clonal evolution of metastatic colorectal cancer.

Tumor heterogeneity and evolution drive treatment resistance in metastatic colorectal cancer (mCRC). Patient-derived xenografts (PDXs) can model mCRC biology; however, their ability to accurately mimic human tumor heterogeneity is unclear. Current genomic studies in mCRC have limited scope and lack matched PDXs. Therefore, the landscape of tumor heterogeneity and its impact on the evolution of metastasis and PDXs remain undefined. We performed whole-genome, deep exome, and targeted validation sequencing of multiple primary regions, matched distant metastases, and PDXs from 11 patients with mCRC. We observed intricate clonal heterogeneity and evolution affecting metastasis dissemination and PDX clonal selection. Metastasis formation followed both monoclonal and polyclonal seeding models. In four cases, metastasis-seeding clones were not identified in any primary region, consistent with a metastasis-seeding-metastasis model. PDXs underrepresented the subclonal heterogeneity of parental tumors. These suggest that single sample tumor sequencing and current PDX models may be insufficient to guide precision medicine.

Long non-coding RNA RAMS11 promotes metastatic colorectal cancer progression.

Colorectal cancer (CRC) is the most common gastrointestinal malignancy in the U.S.A. and approximately 50% of patients develop metastatic disease (mCRC). Despite our understanding of long non-coding RNAs (lncRNAs) in primary colon cancer, their role in mCRC and treatment resistance remains poorly characterized. Therefore, through transcriptome sequencing of normal, primary, and distant mCRC tissues we find 148 differentially expressed RNAs Associated with Metastasis (RAMS). We prioritize RAMS11 due to its association with poor disease-free survival and promotion of aggressive phenotypes in vitro and in vivo. A FDA-approved drug high-throughput viability assay shows that elevated RAMS11 expression increases resistance to topoisomerase inhibitors. Subsequent experiments demonstrate RAMS11-dependent recruitment of Chromobox protein 4 (CBX4) transcriptionally activates Topoisomerase II alpha (TOP2α). Overall, recent clinical trials using topoisomerase inhibitors coupled with our findings of RAMS11-dependent regulation of TOP2α supports the potential use of RAMS11 as a biomarker and therapeutic target for mCRC.

Long non-coding RNA LCAL62 / LINC00261 is associated with lung adenocarcinoma prognosis.

BACKGROUND: More than half of non-small cell lung cancer (NSCLC) patients present with metastatic disease at initial diagnosis with an estimated five-year survival rate of ~5%. Despite advances in understanding primary lung cancer oncogenesis metastatic disease remains poorly characterized. Recent studies demonstrate important roles of long non-coding RNAs (lncRNAs) in tumor physiology and as prognostic markers. Therefore, we present the first transcriptome analysis to identify lncRNAs altered in metastatic lung adenocarcinoma leading to the discovery and characterization of the lncRNA LCAL62 as a prognostic biomarker. PATIENTS AND METHODS: RNA-Seq, microarray, nanoString expression, and clinical data from 1,116 LUAD patients across six independent cohorts and 83 LUAD cell lines were used to discover and evaluate the survival association of metastasis associated lncRNAs. Coexpression and gene set enrichment analyses were used to establish gene regulatory networks and implicate metastasis associated lncRNAs in specific biological processes. RESULTS: Our integrative analysis discovered LCAL62 as the most down-regulated lncRNA in metastasis. Further low LCAL62 expression promoted aggressive phenotypes and regulated genes associated with metastasis (such as metastasis repressor FOXA2). Low LCAL62 expression corresponded to poor overall patient survival across five independent lung adenocarcinoma cohorts (n = 881) including our own nanoString validation cohort. CONCLUSION: We discovered that LCAL62 was down-regulated in lung cancer progression to promote invasive phenotypes, and lower expression was significantly associated with poor patient outcome and aggressive lung adenocarcinoma.

Heparin-based hydrogel scaffolding alters the transcriptomic profile and increases the chemoresistance of MDA-MB-231 triple-negative breast cancer cells.

The tumor microenvironment plays a critical role in the proliferation and chemoresistance of cancer cells. Growth factors (GFs) are known to interact with the extracellular matrix (ECM) via heparin binding sites, and these associations influence cell behavior. In the present study, we demonstrate the ability to define signals presented by the scaffold by pre-mixing growth factors, such as epidermal growth factor, into the heparin-based (HP-B) hydrogel prior to gelation. In the 3D biomimetic microenvironment, breast cancer cells formed spheroids within 24 hours of initial seeding. Despite higher number of proliferating cells in 2D cultures, 3D spheroids exhibited a higher degree of chemoresistance after 72 hours. Further, our RNA sequencing results highlighted the phenotypic changes influenced by solid-phase GF presentation. Wnt/ß-catenin and TGF-ß signaling were upregulated in the cells grown in the hydrogel, while apoptosis, IL2-STAT5 and PI3K-AKT-mTOR signaling were downregulated. With emerging technologies for precision medicine in cancer, this nature of fine-tuning the microenvironment is paramount for cultivation and downstream characterization of primary cancer cells and rare circulating tumor cells (CTCs), and effective screening of chemotherapeutic agents.

Diverse AR Gene Rearrangements Mediate Resistance to Androgen Receptor Inhibitors in Metastatic Prostate Cancer.

PURPOSE: Prostate cancer is the second leading cause of male cancer deaths. Castration-resistant prostate cancer (CRPC) is a lethal stage of the disease that emerges when endocrine therapies are no longer effective at suppressing activity of the androgen receptor (AR) transcription factor. The purpose of this study was to identify genomic mechanisms that contribute to the development and progression of CRPC. EXPERIMENTAL DESIGN: We used whole-genome and targeted DNA-sequencing approaches to identify mechanisms underlying CRPC in an aggregate cohort of 272 prostate cancer patients. We analyzed structural rearrangements at the genome-wide level and carried out a detailed structural rearrangement analysis of the AR locus. We used genome engineering to perform experimental modeling of AR gene rearrangements and long-read RNA sequencing to analyze effects on expression of AR and truncated AR variants (AR-V). RESULTS: AR was among the most frequently rearranged genes in CRPC tumors. AR gene rearrangements promoted expression of diverse AR-V species. AR gene rearrangements occurring in the context of AR amplification correlated with AR overexpression. Cell lines with experimentally derived AR gene rearrangements displayed high expression of tumor-specific AR-Vs and were resistant to endocrine therapies, including the AR antagonist enzalutamide. CONCLUSIONS: AR gene rearrangements are an important mechanism of resistance to endocrine therapies in CRPC.