Widespread and Functional RNA Circularization in Localized Prostate Cancer.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Risk SNP-Mediated Promoter-Enhancer Switching Drives Prostate Cancer through lncRNA PCAT19.

The prostate cancer (PCa) risk-associated SNP rs11672691 is positively associated with aggressive disease at diagnosis. We showed that rs11672691 maps to the promoter of a short isoform of long noncoding RNA PCAT19 (PCAT19-short), which is in the third intron of the long isoform (PCAT19-long). The risk variant is associated with decreased and increased levels of PCAT19-short and PCAT19-long, respectively. Mechanistically, the risk SNP region is bifunctional with both promoter and enhancer activity. The risk variants of rs11672691 and its LD SNP rs887391 decrease binding of transcription factors NKX3.1 and YY1 to the promoter of PCAT19-short, resulting in weaker promoter but stronger enhancer activity that subsequently activates PCAT19-long. PCAT19-long interacts with HNRNPAB to activate a subset of cell-cycle genes associated with PCa progression, thereby promoting PCa tumor growth and metastasis. Taken together, these findings reveal a risk SNP-mediated promoter-enhancer switching mechanism underlying both initiation and progression of aggressive PCa.

CRISPR screen identifies genes that sensitize AML cells to double-negative T-cell therapy.

Acute myeloid leukemia (AML) remains a devastating disease in need of new therapies to improve patient survival. Targeted adoptive T-cell therapies have achieved impressive clinical outcomes in some B-cell leukemias and lymphomas but not in AML. Double-negative T cells (DNTs) effectively kill blast cells from the majority of AML patients and are now being tested in clinical trials. However, AML blasts obtained from ∼30% of patients show resistance to DNT-mediated cytotoxicity; the markers or mechanisms underlying this resistance have not been elucidated. Here, we used a targeted clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) screen to identify genes that cause susceptibility of AML cells to DNT therapy. Inactivation of the Spt-Ada-Gcn5-acetyltransferase (SAGA) deubiquitinating complex components sensitized AML cells to DNT-mediated cytotoxicity. In contrast, CD64 inactivation resulted in resistance to DNT-mediated cytotoxicity. Importantly, the level of CD64 expression correlated strongly with the sensitivity of AML cells to DNT treatment. Furthermore, the ectopic expression of CD64 overcame AML resistance to DNTs in vitro and in vivo. Altogether, our data demonstrate the utility of CRISPR/Cas9 screens to uncover mechanisms underlying the sensitivity to DNT therapy and suggest CD64 as a predictive marker for response in AML patients.

Genomic hallmarks of localized, non-indolent prostate cancer.

Prostate tumours are highly variable in their response to therapies, but clinically available prognostic factors can explain only a fraction of this heterogeneity. Here we analysed 200 whole-genome sequences and 277 additional whole-exome sequences from localized, non-indolent prostate tumours with similar clinical risk profiles, and carried out RNA and methylation analyses in a subset. These tumours had a paucity of clinically actionable single nucleotide variants, unlike those seen in metastatic disease. Rather, a significant proportion of tumours harboured recurrent non-coding aberrations, large-scale genomic rearrangements, and alterations in which an inversion repressed transcription within its boundaries. Local hypermutation events were frequent, and correlated with specific genomic profiles. Numerous molecular aberrations were prognostic for disease recurrence, including several DNA methylation events, and a signature comprised of these aberrations outperformed well-described prognostic biomarkers. We suggest that intensified treatment of genomically aggressive localized prostate cancer may improve cure rates.

CRISPR activation screens identify the SWI/SNF ATPases as suppressors of ferroptosis.

Ferroptosis is an iron-dependent cell death mechanism characterized by the accumulation of toxic lipid peroxides and cell membrane rupture. GPX4 (glutathione peroxidase 4) prevents ferroptosis by reducing these lipid peroxides into lipid alcohols. Ferroptosis induction by GPX4 inhibition has emerged as a vulnerability of cancer cells, highlighting the need to identify ferroptosis regulators that may be exploited therapeutically. Through genome-wide CRISPR activation screens, we identify the SWI/SNF (switch/sucrose non-fermentable) ATPases BRM (SMARCA2) and BRG1 (SMARCA4) as ferroptosis suppressors. Mechanistically, they bind to and increase chromatin accessibility at NRF2 target loci, thus boosting NRF2 transcriptional output to counter lipid peroxidation and confer resistance to GPX4 inhibition. We further demonstrate that the BRM/BRG1 ferroptosis connection can be leveraged to enhance the paralog dependency of BRG1 mutant cancer cells on BRM. Our data reveal ferroptosis induction as a potential avenue for broadening the efficacy of BRM degraders/inhibitors and define a specific genetic context for exploiting GPX4 dependency.

The N6-methyladenosine Epitranscriptomic Landscape of Lung Adenocarcinoma.

Comprehensive m6A epitranscriptome profiling of primary tumors remains largely uncharted. Here, we profiled the m6A epitranscriptome of 10 non-neoplastic lung (NL) tissues and 51 lung adenocarcinoma (LUAD) tumors, integrating the corresponding transcriptome, proteome and extensive clinical annotations. We identified distinct clusters and genes that were exclusively linked to disease progression through m6A modifications. In comparison with NL tissues, we identified 430 transcripts to be hypo-methylated and 222 to be hyper-methylated in tumors. Among these genes, EML4 emerged as a novel metastatic driver, displaying significant hyper-methylation in tumors. m6A modification promoted the translation of EML4, leading to its widespread overexpression in primary tumors. Functionally, EML4 modulated cytoskeleton dynamics through interacting with ARPC1A, enhancing lamellipodia formation, cellular motility, local invasion, and metastasis. Clinically, high EML4 protein abundance correlated with features of metastasis. METTL3 small molecule inhibitor markedly diminished both EML4 m6A and protein abundance, and efficiently suppressed lung metastases in vivo.

DNA methylation modulated genetic variant effect on gene transcriptional regulation.

BACKGROUND: Expression quantitative trait locus (eQTL) analysis has emerged as an important tool in elucidating the link between genetic variants and gene expression, thereby bridging the gap between risk SNPs and associated diseases. We recently identified and validated a specific case where the methylation of a CpG site influences the relationship between the genetic variant and gene expression. RESULTS: Here, to systematically evaluate this regulatory mechanism, we develop an extended eQTL mapping method, termed DNA methylation modulated eQTL (memo-eQTL). Applying this memo-eQTL mapping method to 128 normal prostate samples enables identification of 1063 memo-eQTLs, the majority of which are not recognized as conventional eQTLs in the same cohort. We observe that the methylation of the memo-eQTL CpG sites can either enhance or insulate the interaction between SNP and gene expression by altering CTCF-based chromatin 3D structure. CONCLUSIONS: This study demonstrates the prevalence of memo-eQTLs paving the way to identify novel causal genes for traits or diseases associated with genetic variations.

MYC reshapes CTCF-mediated chromatin architecture in prostate cancer.

MYC is a well characterized oncogenic transcription factor in prostate cancer, and CTCF is the main architectural protein of three-dimensional genome organization. However, the functional link between the two master regulators has not been reported. In this study, we find that MYC rewires prostate cancer chromatin architecture by interacting with CTCF protein. Through combining the H3K27ac, AR and CTCF HiChIP profiles with CRISPR deletion of a CTCF site upstream of MYC gene, we show that MYC activation leads to profound changes of CTCF-mediated chromatin looping. Mechanistically, MYC colocalizes with CTCF at a subset of genomic sites, and enhances CTCF occupancy at these loci. Consequently, the CTCF-mediated chromatin looping is potentiated by MYC activation, resulting in the disruption of enhancer-promoter looping at neuroendocrine lineage plasticity genes. Collectively, our findings define the function of MYC as a CTCF co-factor in three-dimensional genome organization.

An allosteric pan-TEAD inhibitor blocks oncogenic YAP/TAZ signaling and overcomes KRAS G12C inhibitor resistance.

The Hippo pathway is a key growth control pathway that is conserved across species. The downstream effectors of the Hippo pathway, YAP (Yes-associated protein) and TAZ (transcriptional coactivator with PDZ-binding motif), are frequently activated in cancers to drive proliferation and survival. Based on the premise that sustained interactions between YAP/TAZ and TEADs (transcriptional enhanced associate domain) are central to their transcriptional activities, we discovered a potent small-molecule inhibitor (SMI), GNE-7883, that allosterically blocks the interactions between YAP/TAZ and all human TEAD paralogs through binding to the TEAD lipid pocket. GNE-7883 effectively reduces chromatin accessibility specifically at TEAD motifs, suppresses cell proliferation in a variety of cell line models and achieves strong antitumor efficacy in vivo. Furthermore, we uncovered that GNE-7883 effectively overcomes both intrinsic and acquired resistance to KRAS (Kirsten rat sarcoma viral oncogene homolog) G12C inhibitors in diverse preclinical models through the inhibition of YAP/TAZ activation. Taken together, this work demonstrates the activities of TEAD SMIs in YAP/TAZ-dependent cancers and highlights their potential broad applications in precision oncology and therapy resistance.

Transposable elements are a significant contributor to tandem repeats in the human genome.

Sequence repeats are an important phenomenon in the human genome, playing important roles in genomic alteration often with phenotypic consequences. The two major types of repeat elements in the human genome are tandem repeats (TRs) including microsatellites, minisatellites, and satellites and transposable elements (TEs). So far, very little has been known about the relationship between these two types of repeats. In this study, we identified TRs that are derived from TEs either based on sequence similarity or overlapping genomic positions. We then analyzed the distribution of these TRs among TE families/subfamilies. Our study shows that at least 7,276 TRs or 23% of all minisatellites/satellites is derived from TEs, contributing ∼0.32% of the human genome. TRs seem to be generated more likely from younger/more active TEs, and once initiated they are expanded with time via local duplication of the repeat units. The currently postulated mechanisms for origin of TRs can explain only 6% of all TE-derived TRs, indicating the presence of one or more yet to be identified mechanisms for the initiation of such repeats. Our result suggests that TEs are contributing to genome expansion and alteration not only by transposition but also by generating tandem repeats.

Neonatal LTßR signaling is required for the accumulation of eosinophils in the inflamed adult mesenteric lymph node.

Although eosinophils are important contributors to mucosal immune responses, mechanisms that regulate their accumulation in mucosal-associated lymphoid tissues remain ill-defined. Combining bone marrow chimeras and pharmacological inhibition approaches, here we find that lymphotoxin-beta receptor (LTßR) signaling during the neonatal period is required for the accumulation of eosinophils in the mesenteric lymph nodes (MLN) during an enteric viral infection in adult male and female mice. We demonstrate that MLN stromal cells express genes that are important for eosinophil migration and survival, such as Ccl-11 (eotaxin-1), Ccl7, Ccl9, and Cxcl2, and that expression of most of these genes is downregulated as a consequence of neonatal LTßR blockade. We also find that neonatal LTßR signaling is required for the generation of a rotavirus-specific IgA antibody response in the adult MLN, but eosinophils are dispensable for this response. Collectively, our studies reveal a role for neonatal LTßR signaling in regulating eosinophil numbers in the adult MLN.

A noncoding single-nucleotide polymorphism at 8q24 drives IDH1-mutant glioma formation.

Establishing causal links between inherited polymorphisms and cancer risk is challenging. Here, we focus on the single-nucleotide polymorphism rs55705857, which confers a sixfold greater risk of isocitrate dehydrogenase (IDH)-mutant low-grade glioma (LGG). We reveal that rs55705857 itself is the causal variant and is associated with molecular pathways that drive LGG. Mechanistically, we show that rs55705857 resides within a brain-specific enhancer, where the risk allele disrupts OCT2/4 binding, allowing increased interaction with the Myc promoter and increased Myc expression. Mutating the orthologous mouse rs55705857 locus accelerated tumor development in an Idh1R132H-driven LGG mouse model from 472 to 172 days and increased penetrance from 30% to 75%. Our work reveals mechanisms of the heritable predisposition to lethal glioma in ~40% of LGG patients.

CRISPR screens identify cholesterol biosynthesis as a therapeutic target on stemness and drug resistance of colon cancer.

Cancer stem cells (CSCs) are responsible for tumor progression, recurrence, and drug resistance. To identify genetic vulnerabilities of colon cancer, we performed targeted CRISPR dropout screens comprising 657 Drugbank targets and 317 epigenetic regulators on two patient-derived colon CSC-enriched spheroids. Next-generation sequencing of pooled genomic DNAs isolated from surviving cells yielded therapeutic candidates. We unraveled 44 essential genes for colon CSC-enriched spheroids propagation, including key cholesterol biosynthetic genes (HMGCR, FDPS, and GGPS1). Cholesterol biosynthesis was induced in colon cancer tissues, especially CSC-enriched spheroids. The genetic and pharmacological inhibition of HMGCR/FDPS impaired self-renewal capacity and tumorigenic potential of the spheroid models in vitro and in vivo. Mechanistically, HMGCR or FDPS depletion impaired cancer stemness characteristics by activating TGF-ß signaling, which in turn downregulated expression of inhibitors of differentiation (ID) proteins, key regulators of cancer stemness. Cholesterol and geranylgeranyl diphosphate (GGPP) rescued the growth inhibitory and signaling effect of HMGCR/FDPS blockade, implying a direct role of these metabolites in modulating stemness. Finally, cholesterol biosynthesis inhibitors and 5-FU demonstrated antitumor synergy in colon CSC-enriched spheroids, tumor organoids, and xenografts. Taken together, our study unravels novel genetic vulnerabilities of colon CSC-enriched spheroids and suggests cholesterol biosynthesis as a potential target in conjunction with traditional chemotherapy for colon cancer treatment.

HNRNPM controls circRNA biogenesis and splicing fidelity to sustain cancer cell fitness.

High spliceosome activity is a dependency for cancer cells, making them more vulnerable to perturbation of the splicing machinery compared to normal cells. To identify splicing factors important for prostate cancer (PCa) fitness, we performed pooled shRNA screens in vitro and in vivo. Our screens identified heterogeneous nuclear ribonucleoprotein M (HNRNPM) as a regulator of PCa cell growth. RNA- and eCLIP-sequencing identified HNRNPM binding to transcripts of key homeostatic genes. HNRNPM binding to its targets prevents aberrant exon inclusion and backsplicing events. In both linear and circular mis-spliced transcripts, HNRNPM preferentially binds to GU-rich elements in long flanking proximal introns. Mimicry of HNRNPM-dependent linear-splicing events using splice-switching-antisense-oligonucleotides was sufficient to inhibit PCa cell growth. This suggests that PCa dependence on HNRNPM is likely a result of mis-splicing of key homeostatic coding and non-coding genes. Our results have further been confirmed in other solid tumors. Taken together, our data reveal a role for HNRNPM in supporting cancer cell fitness. Inhibition of HNRNPM activity is therefore a potential therapeutic strategy in suppressing growth of PCa and other solid tumors.

An androgen receptor switch underlies lineage infidelity in treatment-resistant prostate cancer.

Cancers adapt to increasingly potent targeted therapies by reprogramming their phenotype. Here we investigated such a phenomenon in prostate cancer, in which tumours can escape epithelial lineage confinement and transition to a high-plasticity state as an adaptive response to potent androgen receptor (AR) antagonism. We found that AR activity can be maintained as tumours adopt alternative lineage identities, with changes in chromatin architecture guiding AR transcriptional rerouting. The epigenetic regulator enhancer of zeste homologue 2 (EZH2) co-occupies the reprogrammed AR cistrome to transcriptionally modulate stem cell and neuronal gene networks-granting privileges associated with both fates. This function of EZH2 was associated with T350 phosphorylation and establishment of a non-canonical polycomb subcomplex. Our study provides mechanistic insights into the plasticity of the lineage-infidelity state governed by AR reprogramming that enabled us to redirect cell fate by modulating EZH2 and AR, highlighting the clinical potential of reversing resistance phenotypes.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer.

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

A computational approach for identification of epitopes in dengue virus envelope protein: a step towards designing a universal dengue vaccine targeting endemic regions.

A major problem in designing vaccine for the dengue virus has been the high antigenic variability in the envelope protein of different virus strains. In this study, a computational approach was adopted to identify a multi-epitope vaccine candidate against dengue virus that may be suitable for large populations in the dengue-endemic regions. Different bioinformatics tools were exploited that helped the identification of a conserved immunological hot-spot in the dengue envelope protein. The tools also rendered the prediction of immunogenicity and population coverage to the proposed 'in silico' vaccine candidate against dengue. A peptide region, spanning 19 amino acids, was identified in the envelope protein which found to be conserved in all four types of dengue viruses. Ten proteasomal cleavage sites were identified within the 19-mer conserved peptide sequence and a total of 8 overlapping putative cytotoxic T cell (CTL) epitopes were identified. The immunogenicity of these epitopes was evaluated in terms of their binding affinities to and dissociation half-time from respective human leukocyte antigen (HLA) molecules. The HLA allele frequencies were studied among populations in the dengue endemic regions and compared with respect to HLA restriction patterns of the overlapping epitopes. The cumulative population coverage for these epitopes as vaccine candidates was high ranging from approximately 80% to 92%. Structural analysis suggested that a 9-mer epitope fitted well into the peptide-binding groove of HLA-A*0201. In conclusion, the 19-mer epitope cluster was shown to have the potential for use as a vaccine candidate against dengue.

Comparative network clustering of direct repeats (DRs) and cas genes confirms the possibility of the horizontal transfer of CRISPR locus among bacteria.

CRISPRs are a diverse family of DNA repeat sequences that are widely distributed among archaea and bacteria. The CRISPR locus is usually composed of three key elements; direct repeats (DRs), spacer sequences and the cas genes. Although recent studies have suggested that spacers may be of extrachromosomal origin, the evolutionary origin of the other two elements of the CRISPR locus has remained unresolved. With the aim to elucidate the evolutionary origin and association of DRs and cas genes of the CRISPR locus, a comparative analysis of the evolutionary network clusters of DRs, cas1 and 16S rRNA genes sequences from 100 different bacteria was conducted. Significant matches of DR and cas1 gene clades imply that these CRISPR components are evolutionary closely linked and potentially evolving simultaneously as a whole locus. On the contrary, the prominent discordance between the CASS (DR and cas1) clades and the 16S rRNA clusters indicates that CRISPR locus has been transferred horizontally as a complete package. Sequence analysis also revealed that DR and cas1 genes are coevolving under analogous evolutionary pressure. This atypical evolutionary pattern also signifies the possibility of horizontal transfer event of CRISPR locus.

Integration of Online Omics-Data Resources for Cancer Research.

The manifestations of cancerous phenotypes necessitate alterations at different levels of information-flow from genome to proteome. The molecular alterations at different information processing levels serve as the basis for the cancer phenotype to emerge. To understand the underlying mechanisms that drive the acquisition of cancer hallmarks it is required to interrogate cancer cells using multiple levels of information flow represented by different omics - such as genomics, epigenomics, transcriptomics, and proteomics. The advantage of multi-omics data integration comes with a trade-off in the form of an added layer of complexity originating from inherently diverse types of omics-datasets that may pose a challenge to integrate the omics-data in a biologically meaningful manner. The plethora of cancer-specific online omics-data resources, if able to be integrated efficiently and systematically, may facilitate the generation of new biological insights for cancer research. In this review, we provide a comprehensive overview of the online single- and multi-omics resources that are dedicated to cancer. We catalog various online omics-data resources such as The Cancer Genome Atlas (TCGA) along with various TCGA-associated data portals and tools for multi-omics analysis and visualization, the International Cancer Genome Consortium (ICGC), Catalogue of Somatic Mutations in Cancer (COSMIC), The Pathology Atlas, Gene Expression Omnibus (GEO), and PRoteomics IDEntifications (PRIDE). By comparing the strengths and limitations of the respective online resources, we aim to highlight the current biological and technological challenges and possible strategies to overcome these challenges. We outline the available schemes for the integration of the multi-omics dimensions for stratifying cancer patients and biomarker prediction based on the integrated molecular-signatures of cancer. Finally, we propose the multi-omics driven systems-biology approaches to realize the potential of precision onco-medicine as the future of cancer research. We believe this systematic review will encourage scientists and clinicians worldwide to utilize the online resources to explore and integrate the available omics datasets that may provide a window of opportunity to generate new biological insights and contribute to the advancement of the field of cancer research.

Chromatin binding of FOXA1 is promoted by LSD1-mediated demethylation in prostate cancer.

FOXA1 functions as a pioneer transcription factor by facilitating the access to chromatin for steroid hormone receptors, such as androgen receptor and estrogen receptor1-4, but mechanisms regulating its binding to chromatin remain elusive. LSD1 (KDM1A) acts as a transcriptional repressor by demethylating mono/dimethylated histone H3 lysine 4 (H3K4me1/2)5,6, but also acts as a steroid hormone receptor coactivator through mechanisms that are unclear. Here we show, in prostate cancer cells, that LSD1 associates with FOXA1 and active enhancer markers, and that LSD1 inhibition globally disrupts FOXA1 chromatin binding. Mechanistically, we demonstrate that LSD1 positively regulates FOXA1 binding by demethylating lysine 270, adjacent to the wing2 region of the FOXA1 DNA-binding domain. Acting through FOXA1, LSD1 inhibition broadly disrupted androgen-receptor binding and its transcriptional output, and dramatically decreased prostate cancer growth alone and in synergy with androgen-receptor antagonist treatment in vivo. These mechanistic insights suggest new therapeutic strategies in steroid-driven cancers.