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.

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.

Proteome-centric studies, although have identified numerous lncRNA-encoded polypeptides, lack differential expression analysis of lncRNA-peptidome across primary tissues, cell lines and cancer states. We established a computational-proteogenomic workflow involving re-processing of publicly available LC-MS/MS data, which facilitated the identification of tissue-specific and universally expressed (UExp) lncRNA-polypeptides across 14 primary human tissues and 11 cell lines. The utility of lncRNA-peptidome as cancer-biomarkers was investigated by re-processing LC-MS/MS data from 92 colon-adenocarcinoma (COAD) and 30 normal colon-epithelium tissues. Intriguingly, a significant upregulation of five lncRNA UExp-polypeptides in COAD tissues was observed. Furthermore, clustering of the UExp-polypeptides led to the classification of COAD patients that coincided with the clinical stratification, underlining the prognostic potential of the UExp-polypeptides. Lastly, we identified differential abundance of the UExp-polypeptides in the plasma of prostate-cancer patients highlighting their potential as plasma-biomarker. The analysis of lncRNA-peptidome may pave the way to identify effective tissue/plasma biomarkers for different cancer types.

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.

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.

Transcription factors (TFs) are spatially and temporally regulated during gut organ specification. Although accumulating evidence shows aberrant reactivation of developmental programs in cancer, little is known about how TFs drive lineage specification in development and cancer. We first defined gastrointestinal tissue-specific chromatin accessibility and gene expression during development, identifying the dynamic epigenetic regulation of SOX family of TFs. We revealed that Sox2 is not only essential for gastric specification, by maintaining chromatin accessibility at forestomach lineage loci, but also sufficient to promote forestomach/esophageal transformation upon Cdx2 deletion. By comparing our gastrointestinal lineage-specific transcriptome to human gastrointestinal cancer data, we found that stomach and intestinal lineage-specific programs are reactivated in Sox2high /Sox9high and Cdx2high cancers, respectively. By analyzing mice deleted for both Sox2 and Sox9, we revealed their potentially redundant roles in both gastric development and cancer, highlighting the importance of developmental lineage programs reactivated by gastrointestinal TFs in cancer.

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.

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.

Neuroendocrine prostate cancer (NEPC), a lethal form of the disease, is characterized by loss of androgen receptor (AR) signaling during neuroendocrine transdifferentiation, which results in resistance to AR-targeted therapy. Clinically, genomically and epigenetically, NEPC resembles other types of poorly differentiated neuroendocrine tumors (NETs). Through pan-NET analyses, we identified ONECUT2 as a candidate master transcriptional regulator of poorly differentiated NETs. ONECUT2 ectopic expression in prostate adenocarcinoma synergizes with hypoxia to suppress androgen signaling and induce neuroendocrine plasticity. ONEUCT2 drives tumor aggressiveness in NEPC, partially through regulating hypoxia signaling and tumor hypoxia. Specifically, ONECUT2 activates SMAD3, which regulates hypoxia signaling through modulating HIF1α chromatin-binding, leading NEPC to exhibit higher degrees of hypoxia compared to prostate adenocarcinomas. Treatment with hypoxia-activated prodrug TH-302 potently reduces NEPC tumor growth. Collectively, these results highlight the synergy between ONECUT2 and hypoxia in driving NEPC, and emphasize the potential of hypoxia-directed therapy for NEPC patients.

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.

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.

Prostate cancer is the second most commonly diagnosed malignancy among men worldwide. Recurrently mutated in primary and metastatic prostate tumors, FOXA1 encodes a pioneer transcription factor involved in disease onset and progression through both androgen receptor-dependent and androgen receptor-independent mechanisms. Despite its oncogenic properties however, the regulation of FOXA1 expression remains unknown. Here, we identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic single-nucleotide variants in primary prostate tumors. We find that deletion and repression of these cis-regulatory elements significantly decreases FOXA1 expression and prostate cancer cell growth. Six of the ten single-nucleotide variants mapping to FOXA1 regulatory plexus significantly alter the transactivation potential of cis-regulatory elements by modulating the binding of transcription factors. Collectively, our results identify cis-regulatory elements within the FOXA1 plexus mutated in primary prostate tumors as potential targets for therapeutic intervention.

Thousands of noncoding somatic single-nucleotide variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1, and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators' tumor cistromes, exemplified at the 8q24 locus harboring both risk variants and somatic SNVs in cis-regulatory elements upregulating MYC expression. However, Massively Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual cis-regulatory elements. Instead, similar to inherited risk variants, SNVs accumulate in cistromes of master transcription regulators required for prostate cancer development.

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.

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.

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.