Jing Wang

Background/Objectives: Cancer cell migration is a hallmark of cancer and is associated with metastasis. While large-scale functional screens have identified regulators of migration, less is known about how intrinsic transcriptional heterogeneity drives highly migratory phenotypes within individual cancer models or whether these transcriptional changes are conserved across different models of varying tissues of origin. This study aims to define shared and cell line-specific transcriptional programs associated with cancer cell migration and analyze their relevance to patient datasets. Methods: Five cancer cell lines across three cancer types (breast, colorectal, and melanoma) were subjected to transwell-based migratory sorting to isolate highly and weakly migratory subpopulations. Bulk RNA sequencing, differential gene expression analysis, Gene Ontology (GO) enrichment, and upstream regulator prediction were performed. Public tumor datasets were analyzed to evaluate gene expression and its association with patient survival. Results: EVA1A was consistently upregulated in all highly migratory (HM) subpopulations. Multiple GO terms were enriched across all cell lines, often driven by distinct gene signatures, indicating convergence at the level of biological processes despite transcriptional divergence. TEAD4 was predicted as an upstream regulator, and increased TEAD4 nuclear localization was observed in four of the five HM subpopulations. EVA1A and TEAD4 expression were elevated in tumors relative to normal tissues, with cancer type-dependent survival outcomes. Conclusions: Migratory selection was accompanied by extensive transcriptional change within each model, yet across five models spanning three tissue types these changes converged on shared biological processes rather than shared genes. Migration-associated phenotypes may therefore be better defined by pathway-level than single-gene analyses, and the clinical relevance of regulatory nodes such as TEAD4 appears conditional on cancer type.

Background/Objectives: SCCOHT is an aggressive and often fatal cancer that belongs to the ~20% of cancers defined by mutations to subunits of the SWI/SNF chromatin remodeling complex. In SCCOHT, mutations to the SMARCA4 gene, which encodes the SWI/SNF ATPase BRG1, are sufficient to impair SWI/SNF function. This single genetic lesion leads to a cascade of events that promote tumorigenesis, some of which may involve the intersection of SWI/SNF with oncogenic pathways such as those regulated by the MYC oncogene. In SCCOHT tumors and other cancers marked by SWI/SNF subunit mutation, MYC target genes are recurrently activated, pointing to a relationship between SWI/SNF and MYC that has yet to be fully explored. Methods: In this study, we investigate the contribution of MYC to SCCOHT biology by performing a combination of chromatin binding and transcriptome assays in genetically engineered SCCOHT cell lines, with subsequent validation using patient tumor expression data. Results: We find that MYC binds to thousands of active promoters in the BIN-67 SCCOHT cell line and that the depletion of MYC results in a broad range of gene expression changes with a notable effect on the expression of genes related to DNA repair. We uncover an MYC-regulated DNA repair gene expression program in BIN-67 cells that is antagonized by BRG1 reintroduction. Finally, we identify a DNA repair gene signature that is upregulated in SCCOHT tumors and in tumors defined by loss of the SWI/SNF subunit SNF5. Conclusions: Collectively, these data implicate MYC as a robust regulator of DNA repair gene expression in SCCOHT and lay a foundation for future studies focused on interrogating the relationship between BRG1 and MYC.

Nascent RNA sequencing tracks primary transcriptional events, making it crucial for studying the immediate regulatory changes of genes and enhancers in response to both endogenous and exogenous stimuli. NRSA is a widely used tool for analyzing nascent transcriptomic data, enabling quantification of transcriptional changes at proximal promoters and gene bodies, estimation of pausing indices, identifying active enhancers, and establishing enhancer–target gene relationships. To improve its functionality and broaden its applicability to diverse organisms and complex study designs, we have developed an enhanced version, eNRSA. Key advancements include adaptive selection of major transcripts, support for any organism with known gene structures, compatibility with complex study designs, and identification of alternative transcription start and termination sites, as well as transcription readthrough events. Additionally, eNRSA achieves a ∼20-fold increase in analysis speed while significantly reducing memory usage. These enhancements make eNRSA a faster, more versatile, and more powerful tool for nascent transcriptome analysis. eNRSA is freely available at https://bioinfo.vanderbilt.edu/eNRSA/.

Mutations in the SWI/SNF chromatin remodeling complex occur in ~20% of cancers. In rhabdoid tumors defined by loss of the SWI/SNF subunit SMARCB1, dysregulation of enhancer-mediated gene expression is pivotal in driving oncogenesis. Enhancer dysregulation in this setting is tied to retention of the SWI/SNF ATPase BRG1—which becomes essential in the absence of SMARCB1—but precisely how BRG1 contributes to this process remains unknown. To characterize how BRG1 participates in chromatin remodeling and gene expression in SMARCB1-deficient cells, we performed a genome-wide characterization of the impact of BRG1 depletion in multiple rhabdoid tumor cell lines. We find that although BRG1-regulated open chromatin sites are distinct at the locus level, the biological characteristics of the loci are very similar, converging on a set of thematically related genes and pointing to the involvement of the AP-1 transcription factor. The open chromatin sites regulated by BRG1 colocalize with histone-marked enhancers and intriguingly include almost all super-enhancers, revealing that BRG1 plays a critical role in maintaining super-enhancer function in this setting. These studies can explain the essentiality of BRG1 to rhabdoid tumor cell identity and survival and implicate the involvement of AP-1 as a critical downstream effector of rhabdoid tumor cell transcriptional programs.

Motivation: Extracellular vesicles (EVs) are produced and released by most cells and are now recognized to play a role in intercellular communication through the delivery of molecular cargo, including proteins, lipids, and RNA. Small RNA sequencing (small RNA-seq) has been widely used to characterize the small RNA content in EVs. However, there is a lack of a systematic assessment of the quality, technical biases, RNA composition, and RNA biotypes enrichment for small RNA profiling of EVs across cell types, biofluids, and conditions. Methods: We collected and reanalyzed small RNA-seq datasets for 2756 samples from 83 studies involving 55 with EVs only and 28 with both EVs and matched donor cells. We assessed their quality by the total number of reads after adapter trimming, the overall alignment rate to the host and non-host genomes, and the proportional abundance of total small RNA and specific biotypes, such as miRNA, tRNA, rRNA, and Y RNA. Results: We found that EV extraction methods varied in their reproducibility in isolating small RNAs, with effects on small RNA composition. Comparing proportional abundances of RNA biotypes between EVs and matched donor cells, we discovered that rRNA and tRNA fragments were relatively enriched, but miRNAs and snoRNA were depleted in EVs. Except for the export of eight miRNAs being context-independent, the selective release of most miRNAs into EVs was study-specific. Conclusion: This work guides quality control and the selection of EV isolation methods and enhances the interpretation of small RNA contents and preferential loading in EVs.

Alternative splicing (AS) of pre-mRNAs is an essential process that regulates gene expression and functional diversity, yet the underlying regulatory mechanisms that control this process during cancer metastasis are not clear. Here, we uncovered that a long non-coding RNA (lncRNA), MIR99AHG, mediates regulation of AS to alter chromatin remodeler function and promote invadopodia formation in colorectal cancer (CRC). We found that MIR99AHG is highly expressed in metastatic CRC cells and patient tissue samples. In vitro and in vivo assays showed that MIR99AHG potently drove CRC cell motility, invasion, and metastasis. MIR99AHG bound to and stabilized the RNA splicing factor PTBP1, which cooperatively increased cassette exon inclusion of the chromatin remodeling gene SMARCA1. Mechanistically, MIR99AHG acted as an address label for PTBP1 to fine-tune its binding pattern on SMARCA1 pre-mRNA, thereby controlling the skipping-to-inclusion splicing switch of SMARCA1 and favouring the production of a long isoform (SMARCA1-L). Canonical SMARCA1 was found to suppress cell migration and invasiveness by inhibiting invadopodia formation, but SMARCA1-L was functionally inert. Clinicaldata revealed that MIR99AHG was positively correlated with PTBP1 and SMARCA1-L in human CRC specimens and predicted patient outcome. Furthermore, we showed that the crosstalk with cancer-associated fibroblasts triggers TGF-β/SAMD signalling in CRC cells and activates MIR99AHG expression. Our findings establish a novel mechanism for a lncRNA that interacts with PTBP1 to regulate AS process, providing a potential therapeutic target and predictive biomarker of CRC.

Slowly cycling/infrequently proliferating tumor cells present a clinical challenge due to their ability to evade treatment. Previous studies established that high levels of SOX2 in both fetal and tumor cells restrict cell proliferation and induce a slowly cycling state. However, the mechanisms through which elevated SOX2 levels inhibit tumor cell proliferation have not been identified. To identify common mechanisms through which SOX2 elevation restricts tumor cell proliferation, we initially performed RNA-seq using two diverse tumor cell types. SOX2 elevation in both cell types downregulated MYC target genes. Consistent with these findings, elevating SOX2 in five cell lines representing three different human cancer types decreased MYC expression. Importantly, the expression of a dominant-negative MYC variant, omomyc, recapitulated many of the effects of SOX2 on proliferation, cell cycle, gene expression, and biosynthetic activity. We also demonstrated that rescuing MYC activity in the context of elevated SOX2 induces cell death, indicating that the downregulation of MYC is a critical mechanistic step necessary to maintain survival in the slowly cycling state induced by elevated SOX2. Altogether, our findings uncover a novel SOX2:MYC signaling axis and provide important insights into the molecular mechanisms through which SOX2 elevation induces a slowly cycling proliferative state.

Colorectal cancer (CRC) is one of the most common malignancies worldwide with poor prognosis. Studies have showed that abnormal microRNA (miRNA) expression can affect CRC pathogenesis and development through targeting critical genes in cellular system. However, it is unclear about which miRNAs play central roles in CRC’s pathogenesis and how they interact with transcription factors (TFs) to regulate the cancer-related genes. To address this issue, we systematically explored the major regulation motifs, namely feed-forward loops (FFLs), that consist of miRNAs, TFs and CRC-related genes through the construction of a miRNA-TF regulatory network in CRC. First, we compiled CRC-related miRNAs, CRC-related genes, and human TFs from multiple data sources. Second, we identified 13,123 3-node FFLs including 25 miRNA-FFLs, 13,005 TF-FFLs and 93 composite-FFLs, and merged the 3-node FFLs to construct a CRC-related regulatory network. The network consists of three types of regulatory subnetworks (SNWs): miRNA-SNW, TF-SNW, and composite-SNW. To enhance the accuracy of the network, the results were filtered by using The Cancer Genome Atlas (TCGA) expression data in CRC, whereby we generated a core regulatory network consisting of 58 significant FFLs. We then applied a hub identification strategy to the significant FFLs and found 5 significant components, including two miRNAs (hsa-miR-25 and hsa-miR-31), two genes (ADAMTSL3 and AXIN1) and one TF (BRCA1). The follow up prognosis analysis indicated all of the 5 significant components having good prediction of overall survival of CRC patients. In summary, we generated a CRC-specific miRNA-TF regulatory network, which is helpful to understand the complex CRC regulatory mechanisms and guide clinical treatment. The discovered 5 regulators might have critical roles in CRC pathogenesis and warrant future investigation.

Transcription Factors (TFs), essential for many cellular processes, generally work coordinately to induce transcriptional change in response to internal and external signals. Disrupted cooperation between TFs, leading to dysregulation of target genes, contributes to the pathogenesis of many diseases, including cancer. Although the aberrant activation of individual TFs and the functional effects have been widely studied, the perturbation of TF cooperativity in cancer has rarely been explored. We used TF co-expression as proxy as cooperativity and performed a large-scale study on disrupted TF cooperation across seven cancer types. While the connectivity of downstream effectors, like metabolic genes and TF targets, were more or similarly disrupted than/with non-TFs, the cooperativity of TFs (upstream regulators) were consistently less disturbed in all studied cancer types. Highly coordinated TFs in normal, however, generally lost that cooperation in cancer. Although different types of cancer shared very few TF pairs with highly disrupted cooperation, the cooperativity of interferon regulatory factors (IRF) was highly disrupted in six cancer types. Specifically, the cooperativity of IRF8 was highly perturbed in lung cancer, which was further validated by two independent lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD) datasets. More interestingly, the cooperativity of IRF8 was markedly associated with tumor progression and even contributed to the patient survival independent of tumor stage. Our findings underscore the far more important role of TF cooperativity in tumorigenesis than previously appreciated. Disrupted cooperation of TFs provides potential clinical utility as prognostic markers for predicting the patient survival.

It is well known that carcinogenesis is in part dictated by dysregulated transcription events and signal pathways. Large-scale transcriptional profiling studies in each cancer type have reported aberrant gene expression associated with cancer development. However, common and specific patterns altered across cancer types, especially the contribution of transcriptional and post-transcriptional regulators, are rarely explored. Using transcriptional profiles from matched tumor and normal samples in the Cancer Genome Atlas pan-cancer dataset, we performed a comprehensive analysis on the altered expression across 9 cancer types, focusing on transcriptional and post-transcriptional regulators and cancer-related genes. As we expected, the transcription of cancer-related genes was significantly deregulated in tumor vs. normal across all cancer types. Surprisingly, the expression of RNA-binding proteins (RBPs), master regulators of post-transcriptional gene expression, was also significantly changed across most studied cancer types. Although the expression of RBPs was not as strongly deregulated as cancer-related genes, their direct interaction partners are enriched by cancer-related genes, suggesting the cascade regulation effect of RBPs. Integrating genetic and epigenetic profiles found that deregulated RBPs were frequently caused by genetic rather than epigenetic alterations. Furthermore, tissue-specific genes were under-expressed in tumor vs. normal across all cancer types except prostate cancer. Dysregulated transcription across cancer types reveals the importance of RBPs in carcinogenesis. The aberrant expression of RBPs is caused by genetic alterations and spreads their effect to cancer-related genes. In addition, disruption of tissue-specific genes contributes to the corresponding cancer pathology.