Treatment of triple-negative breast cancer (TNBC) is challenging because of its “COLD” tumor immunosuppressive microenvironment (TIME). Here, we present a hydrogel-mediated localized delivery of a combination of docetaxel (DTX) and carboplatin (CPT) (called DTX-CPT-Gel therapy) that ensured enhanced anticancer effect and tumor regression on multiple murine syngeneic and xenograft tumor models. DTX-CPT-Gel therapy modulated the TIME by an increase of antitumorigenic M1 macrophages, attenuation of myeloid-derived suppressor cells, and increase of granzyme B + CD8 + T cells. DTX-CPT-Gel therapy elevated ceramide levels in tumor tissues that activated the protein kinase R (PKR)–like endoplasmic reticulum kinase (PERK)–mediated unfolded protein response (UPR). This UPR-mediated activation of apoptotic cell death led to release of damage-associated molecular patterns, thereby activating the immunogenic cell death that could even clear the metastatic tumors. This study provides a promising hydrogel-mediated platform for DTX-CPT therapy that induces tumor regression and effective immune modulation and, therefore, can be explored further for treatment of TNBC.
Single nucleotide polymorphisms (SNPs) impacting the alternative splicing (AS) process (sQTLs) or isoform expression (iso-eQTL) are implicated as important cancer regulatory elements. To find the sQTL and iso-eQTL, we retrieved prostate cancer (PrCa) tissue RNA-seq and genotype data originating from 385 PrCa European patients from The Cancer Genome Atlas. We conducted RNA-seq analysis with isoform-based and splice event-based approaches. The MatrixEQTL was used to identify PrCa-associated sQTLs and iso-eQTLs. The overlap between sQTL and iso-eQTL with GWAS loci and those that are differentially expressed between cancer and normal tissue were identified. The cis-acting associations (FDR < 0.05) for PrCa-risk SNPs identified 42, 123, and 90 PrCa-associated cassette exons, intron retention, and mRNA isoforms belonging to 25, 95, and 83 genes, respectively; while assessment of trans-acting association (FDR < 0.05) yielded 59, 65, and 196 PrCa-associated cassette exons, intron retention and mRNA isoforms belonging to 35, 55, and 181 genes, respectively. The results suggest that functional PrCa-associated SNPs can play a role in PrCa genesis by making an important contribution to the dysregulation of AS and, consequently, impacting the expression of the mRNA isoforms.
We performed multi-omic profiling of epidermal keratinocytes, precancerous actinic keratoses, and squamous cell carcinomas to understand the molecular transitions during skin carcinogenesis. Single-cell mutational analyses of normal skin cells showed that most keratinocytes have remarkably low mutation burdens, despite decades of sun exposure, however keratinocytes with TP53 or NOTCH1 mutations had substantially higher mutation burdens. These observations suggest that wild-type keratinocytes (i.e. without pathogenic mutations) are able to withstand high dosages of cumulative UV radiation, but certain pathogenic mutations break these adaptive mechanisms, priming keratinocytes for transformation by increasing their mutation rate. Mutational profiling of squamous cell carcinomas adjacent to actinic keratoses revealed TERT promoter and CDKN2A mutations emerging in actinic keratoses, whereas additional mutations inactivating ARID2 and activating the MAPK-pathway delineated the transition to squamous cell carcinomas. Surprisingly, actinic keratoses were often not related to their neighboring squamous cell carcinoma, indicating that collisions of unrelated neoplasms are common in the skin. Spatial variation in gene expression patterns was common in both tumor and immune cells, with high expression of checkpoint molecules at the invasive front of tumors. In conclusion, this study catalogues the key events during the evolution of cutaneous squamous cell carcinoma.
Alternative splicing (AS) is a key post-transcriptional modification that helps in increasing protein diversity. Almost 90% of the protein-coding genes in humans are known to undergo AS and code for different transcripts. Some transcripts are associated with diseases such as breast cancer, lung cancer and glioblastoma. Hence, these transcripts can serve as novel therapeutic and prognostic targets for drug discovery. Herein, we have developed a pipeline, Finding Alternative Splicing Events (FASE), as the R package that includes modules to determine the structure and concentration of transcripts using differential AS. To predict the correct structure of expressed transcripts in given conditions, FASE combines the AS events with the information of exons, introns and junctions using graph theory. The estimated concentration of predicted transcripts is reported as the relative expression in terms of log2CPM. Using FASE, we were able to identify several unique transcripts of EMILIN1 and SLK genes in the TCGA-BRCA data, which were validated using RT-PCR. The experimental study demonstrated consistent results, which signify the high accuracy and precision of the developed methods. In conclusion, the developed pipeline, FASE, can efficiently predict novel transcripts that are missed in general transcript-level differential expression analysis. It can be applied selectively from a single gene to simple or complex genome even in multiple experimental conditions for the identification of differential AS-based biomarkers, prognostic targets and novel therapeutics.