Search

Abstract It has been believed that immunosenescence plays a crucial role in tumorigenesis and cancer therapy. Nevertheless, there is still a lack of understanding regarding its role in determining clinical outcomes and therapy selection for gastric cancer patients, due to the lack of a feasible immunosenescence signature. Therefore, this research aims to develop a gene signature based on immunosenescence, which is used for stratification of gastric cancer. By integrative analysis of bulk transcriptome and single-cell data, we uncovered immunosenescence features in gastric cancer. Random forest algorithm was used to select hub genes and multivariate Cox algorithm was applied to construct a scoring system to evaluate the prognosis and the response to immunotherapy and chemotherapy. The Cancer Genome Atlas of Stomach Adenocarcinoma (TCGA-STAD) cohort was implemented as the training cohort and two independent cohorts from the Gene Expression Omnibus (GEO) database were used for validation. The model was further tested by our Fudan cohort. In this study, immunosenescence was identified as a hallmark of gastric cancer that is linked with transcriptomic features, genomic variations, and distinctive tumor microenvironment (TME). Four immunosenescence genes, including APOD, ADIPOR2, BRAF, and C3, were screened out to construct a gene signature for risk stratification. Higher risk scores indicated strong predictive power for poorer overall survival. Notably, the risk score signature could reliably predict response to chemotherapy and immunotherapy, with patients with high scores benefiting from immunotherapy and patients with low scores responding to chemotherapy. We report immunosenescence as a hitherto unheralded hallmark of gastric cancer that affects prognosis and treatment efficiency.

Abstract Background Colorectal cancer (CRC) is one of the most threatening tumors in the world, and chemotherapy remains dominant in the treatment of metastatic CRC (mCRC) patients. The purpose of this study was to develop a biomarker panel to predict the response of the first line chemotherapy in mCRC patients. Methods Totally 190 mCRC patients treated with FOLFOX or XEOLX chemotherapy in 3 different institutions were included. We extracted the plasma extracellular vesicle (EV) RNA, performed RNA sequencing, constructed a model and generated a signature through shrinking the number of variables by the random forest algorithm and the least absolute shrinkage and selection operator (LASSO) algorithm in the training cohort (n = 80). We validated it in an internal validation cohort (n = 62) and a prospective external validation cohort (n = 48). Results We established a signature consisted of 22 EV RNAs which could identify responders, and the area under the receiver operating characteristic curve (AUC) values was 0.986, 0.821, and 0.816 in the training, internal validation, and external validation cohort respectively. The signature could also identify the progression-free survival (PFS) and overall survival (OS). Besides, we constructed a 7-gene signature which could predict tumor response to first-line oxaliplatin-containing chemotherapy and simultaneously resistance to second-line irinotecan-containing chemotherapy. Conclusions The study was first to develop a signature of EV-derived RNAs to predict the response of the first line chemotherapy in mCRC with high accuracy using a non-invasive approach, indicating that the signature could help to select the optimal regimen for mCRC patients.

Abstract Background Esophageal squamous cell carcinoma (ESCC) is a highly heterogeneous cancer that lacks comprehensive understanding and effective treatment. Although multi-omics study has revealed features and underlying drivers of advanced ESCC, research on molecular characteristics of the early stage ESCC is quite limited. Materials and methods We presented characteristics of genomics and transcriptomics in 10 matched pairs of tumor and normal tissues of early ESCC patients in the China region. Results We identified the specific patterns of cancer gene mutations and copy number variations. We also found a dramatic change in the transcriptome, with more than 4,000 genes upregulated in cancer. Among them, more than one-third of HOX family genes were specifically and highly expressed in early ESCC samples of China and validated by RT-qPCR. Gene regulation network analysis indicated that alteration of Hox family genes promoted the proliferation and metabolism remodeling of early ESCC. Conclusions We characterized the genomic and transcriptomic landscape of 10 paired normal adjacent and early ESCC tissues in the China region, and provided a new perspective to understand the development of ESCC and insight into potential prevention and diagnostic targets for the management of early ESCC in China.

Abstract Chronic liver diseases usually developed through stepwise pathological transitions under the persistent risk factors. The molecular changes during liver transitions are pivotal to improve liver diagnostics and therapeutics yet still remain elusive. Cumulative large-scale liver transcriptomic studies have been revealing molecular landscape of various liver conditions at bulk and single-cell resolution, however, neither single experiment nor databases enabled thorough investigations of transcriptomic dynamics along the progression of liver diseases. Here we establish GepLiver, a longitudinal and multidimensional liver expression atlas integrating expression profiles of 2469 human bulk tissues, 492 mouse samples, 409,775 single cells from 347 human samples and 27 liver cell lines spanning 16 liver phenotypes with uniformed processing and annotating methods. Using GepLiver, we have demonstrated dynamic changes of gene expression, cell abundance and crosstalk harboring meaningful biological associations. GepLiver can be applied to explore the evolving expression patterns and transcriptomic features for genes and cell types respectively among liver phenotypes, assisting the investigation of liver transcriptomic dynamics and informing biomarkers and targets for liver diseases.

Background The heterogeneity of colorectal cancer (CRC) is the main cause of the disparity of drug sensitivity and the variability of prognosis. Pyroptosis is closely associated with the development and prognosis of various tumors, including CRC. Dividing CRC into distinct subgroups based on pyroptosis is a worthwhile topic for improving the precision treatment and prognosis prediction of CRC. Methods We classified patients into two clusters using the consensus clustering based on the pyroptosis-related genes (PRGs). Next, the prognostic signature was developed with LASSO regression analysis using the screened genes from differentially expressed genes (DEGs) by univariate and multivariate Cox analyses. According to the pyroptosis-related score (PR score) calculated with the signature, patients belonged to two groups with distinct prognosis. Moreover, we assessed the immune profile to explore the relationship between the signature and immunological characteristics. Two single cell sequencing databases were adopted for further exploration of tumor immune microenvironment (TME). In addition, we applied our own cohort and Drugbank to explore the correlation of the signature and clinical therapies. We also studied the expression of key genes by immunohistochemistry. Results The signature performed well in predicting the prognosis of CRC as the high area under curve (AUC) value demonstrated. Patients with a higher PR score had poorer prognosis and higher expression of immune checkpoints but more abundant infiltration of immune cells. Combining with the indicator of therapeutic analysis, they might benefit more from immune checkpoint blockade (ICB) and neo-adjuvant chemoradiotherapy (nCRT). Conclusion In conclusion, our study is based on genomics and transcriptomics to investigate the role of PRGs in CRC. We have established a prognostic signature and integrated single-cell data to study the relationship between the signature with the TME in CRC. Its clinical application in reliable prediction of prognosis and personalized treatment was validated by public and own sequencing cohort. It provided a new insight for the personalized treatment of CRC.

Abstract One major characteristic of tumor cells is the aberrant activation of epigenetic regulatory elements, which remodel the tumor transcriptome and ultimately promote cancer progression and drug resistance. However, the oncogenic functions and mechanisms of ovarian cancer (OC) remain elusive. Here, super‐enhancer (SE) regulatory elements that are aberrantly activated in OC are identified and it is found that SEs drive the relative specific expression of the transcription factor KLF5 in OC patients and poly(ADP‐ribose) polymerase inhibitor (PARPi)‐resistant patients. KLF5 expression is associated with poor outcomes in OC patients and can drive tumor progression in vitro and in vivo. Mechanistically, KLF5 forms a transcriptional complex with EHF and ELF3 and binds to the promoter region of RAD51 to enhance its transcription, strengthening the homologous recombination repair (HRR) pathway. Notably, the combination of suberoylanilide hydroxamic acid (SAHA) and olaparib significantly inhibits tumor growth and metastasis of PARPi‐resistant OC cells with high KLF5. In conclusion, it is discovered that SEs‐driven KLF5 is a key regulatory factor in OC progression and PARPi resistance; and potential therapeutic strategies for OC patients with PARPi resistance and high KLF5 are identified.

Abstract Aldolase A (ALDOA), a crucial glycolytic enzyme, is often aberrantly expressed in various types of cancer. Although ALDOA has been reported to play additional roles beyond its conventional enzymatic role, its nonmetabolic function and underlying mechanism in cancer progression remain elusive. Here, it is shown that ALDOA promotes liver cancer growth and metastasis by accelerating mRNA translation independent of its catalytic activity. Mechanistically, ALDOA interacted with insulin‐ like growth factor 2 mRNA‐binding protein 1 (IGF2BP1) to facilitate its binding to m6A‐modified eIF4G mRNA, thereby increasing eIF4G protein levels and subsequently enhancing overall protein biosynthesis in cells. Importantly, administration of GalNAc‐conjugated siRNA targeting ALDOA effectively slows the tumor growth of orthotopic xenografts. Collectively, these findings uncover a previously unappreciated nonmetabolic function of ALDOA in modulating mRNA translation and highlight the potential of specifically targeting ALDOA as a prospective therapeutic strategy in liver cancer.

Abstract Background Long noncoding RNAs (lncRNAs) have driven research focused on their effects as oncogenes or tumor suppressors involved in carcinogenesis. However, the functions and mechanisms of most lncRNAs in colorectal cancer (CRC) remain unclear. Methods The expression of DLGAP1-AS2 was assessed by quantitative RT-PCR in multiple CRC cohorts. The impacts of DLGAP1-AS2 on CRC growth and metastasis were evaluated by a series of in vitro and in vivo assays. Furthermore, the underlying mechanism of DLGAP1-AS2 in CRC was revealed by RNA pull down, RNA immunoprecipitation, RNA sequencing, luciferase assays, chromatin immunoprecipitation, and rescue experiments. Results We discovered that DLGAP1-AS2 promoted CRC tumorigenesis and metastasis by physically interacting with Elongin A (ELOA) and inhibiting its protein stability by promoting tripartite motif containing 21 (Trim21)-mediated ubiquitination modification and degradation of ELOA. In particular, we revealed that DLGAP1-AS2 decreases phospholysine phosphohistidine inorganic pyrophosphate phosphatase (LHPP) expression by inhibiting ELOA-mediated transcriptional activating of LHPP and thus blocking LHPP-dependent suppression of the AKT signaling pathway. In addition, we also demonstrated that DLGAP1-AS2 was bound and stabilized by cleavage and polyadenylation specificity factor (CPSF2) and cleavage stimulation factor (CSTF3). Conclusions The discovery of DLGAP1-AS2, a promising prognostic biomarker, reveals a new dimension into the molecular pathogenesis of CRC and provides a prospective treatment target for this disease.

Abstract Medullary thyroid carcinoma (MTC) is a rare neuroendocrine malignancy derived from parafollicular cells (C cells) of the thyroid. Here we presented a comprehensive multi-omics landscape of 102 MTCs through whole-exome sequencing, RNA sequencing, DNA methylation array, proteomic and phosphoproteomic profiling. Integrated analyses identified BRAF and NF1 as novel driver genes in addition to the well-characterized RET and RAS proto-oncogenes. Proteome-based stratification of MTCs revealed three molecularly heterogeneous subtypes named as: (1) Metabolic, (2) Basal and (3) Mesenchymal, which are distinct in genetic drivers, epigenetic modification profiles, clinicopathologic factors and clinical outcomes. Furthermore, we explored putative therapeutic targets of each proteomic subtype, and found that two tenascin family members TNC/TNXB might serve as potential prognostic biomarkers for MTC. Collectively, our study expands the knowledge of MTC biology and therapeutic vulnerabilities, which may serve as an important resource for future investigation on this malignancy.

Abstract Gastric cancer (GC) is the second cause of cancer-related death and metastasis is an important cause of death. Considering difficulties in searching for metastatic driver mutations, we tried a novel strategy here. We conducted an integrative genomic analysis on GC and identified early drivers lead to metastasis. Whole-exome sequencing (WES), transcriptomes sequencing and targeted-exome sequencing (TES) were performed on tumors and matched normal tissues from 432 Chinese GC patients, especially the comparative analysis between higher metastatic-potential (HMP) group with T1 stage and lymph-node metastasis, and lower metastatic-potential (LMP) group without lymph-nodes or distant metastasis. HMP group presented higher mutation load and heterogeneity, enrichment in immunosuppressive signaling, more immune cell infiltration than LMP group. An integrated mRNA-lncRNA signature based on differentially expressed genes was constructed and its prognostic value was better than traditional TNM stage. We identified 176 candidate prometastatic mutations by WES and selected 8 genes for following TES. Mutated TP53 and MADCAM1 were significantly associated with poor metastasis-free survival. We further demonstrated that mutated MADCAM1 could not only directly promote cancer cells migration, but also could trigger tumor metastasis by establishing immunosuppressive microenvironment, including promoting PD-L1-mediated immune escape and reprogramming tumor-associated macrophages by regulating CCL2 through Akt/mTOR axis. In conclusion, GCs with different metastatic-potential are distinguishable at the genetic level and we revealed a number of potential metastatic driver mutations. Driver mutations in early-onset metastatic GC could promote metastasis by establishing an immunosuppressive microenvironment. This study provided possibility for future target therapy of GC.

Identification of clinically applicable molecular subtypes of pancreatic ductal adenocarcinoma (PDAC) is crucial to improving patient outcomes. However, the traditional tissue-dependent transcriptional subtyping strategies are invasive and not amenable to routine clinical evaluation. In this study, we developed a circulating extracellular vesicle (cEV) long RNA (exLR)-based PDAC subtyping method and provided exLR-derived signatures for predicting immunogenic features and clinical outcomes in PDAC. We enrolled 426 individuals, among which 227 PDACs served as an internal cohort, 118 PDACs from two other medical centers served as an independent validation cohort, and 81 healthy individuals served as the control. ExLR sequencing was performed on all plasma samples. We found that PDAC could be categorized into three subtypes based on plasma exLR profiles. Each subpopulation showed its own molecular features and was associated with patient clinical prognosis. The immunocyte-derived cEV fractions were altered among PDAC subtypes and interconnected with tumor-infiltrating lymphocytes in cancerous tissue. Additionally, we found a significant concordance of immunoregulators between tissue and blood EVs, and we harvested potential PDAC therapeutic targets. Most importantly, we constructed a nine exLR-derived, tissue-applicable signature for prognostic assessment of PDAC. The circulating exLR-based features may offer an attractive platform for personalized treatment and predicting patient outcomes in multiple types of cancer.

Abstract A large number RNAs are enriched and stable in extracellular vesicles (EVs), and they can reflect their tissue origins and are suitable as liquid biopsy markers for cancer diagnosis and treatment efficacy prediction. In this study, we used extracellular vesicle long RNA (exLR) sequencing to characterize the plasma-derived exLRs from 112 breast cancer patients, 19 benign patients and 41 healthy participants. The different exLRs profiling was found between the breast cancer and non-cancer groups. Thus, we constructed a breast cancer diagnostic signature which showed high accuracy with an area under the curve (AUC) of 0.960 in the training cohort and 0.900 in the validation cohort. The signature was able to identify early stage BC (I/II) with an AUC of 0.940. Integrating the signature with breast imaging could increase the diagnosis accuracy for breast cancer patients. Moreover, we enrolled 58 patients who received neoadjuvant treatment and identified an exLR (exMSMO1), which could distinguish pathological complete response (pCR) patients from non-pCR with an AUC of 0.790. Silencing MSMO1 could significantly enhance the sensitivity of MDA-MB-231 cells to paclitaxel and doxorubicin through modulating mTORC1 signaling pathway. This study demonstrated the value of exLR profiling to provide potential biomarkers for early detection and treatment efficacy prediction of breast cancer.

Abstract PTEN is a crucial tumor suppressor and loss of PTEN protein is involved in various cancers. However, the detailed molecular mechanisms of PTEN loss in cancers remain elusive, especially the involvement of lncRNAs. Here, lncRNA RP11-295G20.2 is found to be significantly upregulated in hepatocellular carcinoma (HCC) and promotes the growth of liver cancer cells both in vitro and in vivo. Furthermore, RP11-295G20.2 inhibits autophagy in liver cancer cells. Interestingly, RP11-295G20.2 directly binds to the PTEN protein and leads to its degradation. RP11-295G20.2 expression is inversely correlated with PTEN protein expression in 82 TCGA/TCPA-LIHC samples. Surprisingly, RP11-295G20.2-induced PTEN degradation occurs through the lysosomal pathway instead of the proteasome pathway. RP11-295G20.2 binds to the N terminus of PTEN and facilitates the interaction of p62 with PTEN. Thus, PTEN is translocated into lysosomes and degraded. RP11-295G20.2 also influences AKT phosphorylation and forkhead box O 3a (FOXO3a) translocation into the nucleus, in turn regulating the transcription of autophagy-related genes. Collectively, RP11-295G20.2 directly binds to PTEN and enables its lysosomal degradation. This newly identified RP11-295G20.2/PTEN axis reveals an unexplored molecular mechanism regarding PTEN loss in liver cancer and might provide new therapeutic benefits for liver cancer patients.

Abstract Background Viral infections are prevalent in human cancers and they have great diagnostic and theranostic values in clinical practice. Recently, their potential of shaping the tumor immune microenvironment (TIME) has been related to the immunotherapy of human cancers. However, the landscape of viral expressions and immune status in human cancers remains incompletely understood. Methods We developed a next-generation sequencing (NGS)-based pipeline to detect viral sequences from the whole transcriptome and used machine learning algorithms to classify different TIME subtypes. Results We revealed a pan-cancer landscape of viral expressions in human cancers where 9 types of viruses were detected in 744 tumors of 25 cancer types. Viral infections showed different tissue tendencies and expression levels. Multi-omics analyses further revealed their distinct impacts on genomic, transcriptomic and immune responses. Epstein-Barr virus (EBV)-infected stomach adenocarcinoma (STAD) and Human Papillomavirus (HPV)-infected head and neck squamous cell carcinoma (HNSC) showed decreased genomic variations, significantly altered gene expressions, and effectively triggered anti-viral immune responses. We identified three TIME subtypes, in which the “Immune-Stimulation” subtype might be the promising candidate for immunotherapy. EBV-infected STAD and HPV-infected HNSC showed a higher frequency of the “Immune-Stimulation” subtype. Finally, we constructed the eVIIS pipeline to simultaneously evaluate viral infection and immune status in external datasets. Conclusions Viral infections are prevalent in human cancers and have distinct influences on hosts. EBV and HPV infections combined with the TIME subtype could be promising biomarkers of immunotherapy in STAD and HNSC, respectively. The eVIIS pipeline could be a practical tool to facilitate clinical practice and relevant studies.

Abstract Background Small nucleolar RNA host gene (SNHG) long noncoding RNAs (lncRNAs) are frequently dysregulated in human cancers and involved in tumorigenesis and progression. SNHG17 has been reported as a candidate oncogene in several cancer types, however, its regulatory role in colorectal cancer (CRC) is unclear. Methods SNHG17 expression in multiple CRC cohorts was assessed by RT-qPCR or bioinformatic analyses. Cell viability was evaluated using Cell Counting Kit-8 (CCK-8) and colony formation assays. Cell mobility and invasiveness were assessed by Transwell assays. Tumor xenograft and metastasis models were applied to confirm the effects of SNHG17 on CRC tumorigenesis and metastasis in vivo. Immunohistochemistry staining was used to measure protein expression in cancer tissues. RNA pull-down, RNA immunoprecipitation, chromatin immunoprecipitation, and dual luciferase assays were used to investigate the molecular mechanism of SNHG17 in CRC. Results Using multiple cohorts, we confirmed that SNHG17 is aberrantly upregulated in CRC and correlated with poor survival. In vitro and in vivo functional assays indicated that SNHG17 facilitates CRC proliferation and metastasis. SNHG17 impedes PES1 degradation by inhibiting Trim23-mediated ubiquitination of PES1. SNHG17 upregulates FOSL2 by sponging miR-339-5p, and FOSL2 transcription activates SNHG17 expression, uncovering a SNHG17-miR-339-5p-FOSL2-SNHG17 positive feedback loop. Conclusions We identified SNHG17 as an oncogenic lncRNA in CRC and identified abnormal upregulation of SNHG17 as a prognostic risk factor for CRC. Our mechanistic investigations demonstrated, for the first time, that SNHG17 promotes tumor growth and metastasis through two different regulatory mechanisms, SNHG17-Trim23-PES1 axis and SNHG17-miR-339-5p-FOSL2-SNHG17 positive feedback loop, which may be exploited for CRC therapy.

Background Despite the success of HER2-targeted therapy in achieving prolonged survival in approximately 50% of treated individuals, treatment resistance is still an important challenge for HER2+ breast cancer (BC) patients. The influence of both adaptive and innate immune responses on the therapeutic outcomes of HER2+BC patients has been extensively demonstrated.Methods Long non-coding RNAs expressed in non-pathological complete response (pCR) HER2 positive BC were screened and validated by RNA-seq. Survival analysis were made by Kaplan-Meier method. Cell death assay and proliferation assay were performed to confirm the phenotype of LINC00624. RT-qPCR and western blot were used to assay the IFN response. Xenograft mouse model were used for in vivo confirmation of anti-neu treatment resistance. RNA pull-down and immunoblot were used to confirm the interaction of ADAR1 and LINC00624. ADAR1 recombinant protein were purified from baculovirus expression system. B16-OVA cells were used to study antigen presentation both in vitro and in vivo. Flow cytometry was used to determine the tumor infiltrated immune cells of xenograft model. Antisense oligonucleotides (ASOs) were used for in vivo treatment.Results In this study, we found that LINC00624 blocked the antitumor effect of HER2- targeted therapy both in vitro and in vivo by inhibiting type I interferon (IFN) pathway activation. The double-stranded RNA-like structure of LINC00624 can bind and be edited by the adenosine (A) to inosine (I) RNA-editing enzyme adenosine deaminase RNA specific 1 (ADAR1), and this editing has been shown to release the growth inhibition and attenuate the innate immune response caused by the IFN response. Notably, LINC00624 promoted the stabilization of ADAR1 by inhibiting its ubiquitination-induced degradation triggered by β-TrCP. In contrast, LINC00624 inhibited major histocompatibility complex (MHC) class I antigen presentation and limited CD8+T cell infiltration in the cancer microenvironment, resulting in immune checkpoint blockade inhibition and anti-HER2 treatment resistance mediated through ADAR1.Conclusions In summary, these results suggest that LINC00624 is a cancer immunosuppressive lncRNA and targeting LINC00624 through ASOs in tumors expressing high levels of LINC00624 has great therapeutic potential in future clinical applications.

Background Prostate cancer (PCa) is a major type of cancer in man worldwide. Androgen deprivation therapy (ADT) and the next‐generation androgen receptor (AR) pathway inhibitors have acquired great success in treating PCa. However, patients treated with ADT or AR targeted therapy are inevitably developing into castration‐resistant prostate cancer (CRPC) or becoming drug resistance. The role of mRNA 5‐methylcytosine (m5C) modification in cancers is largely unknown. This study aimed to explore the role of the m5C methyltransferase NSUN2 in Prostate cancer (PCa). Methods The expression of NSUN2 and its clinicopathological impact were evaluated in PCa cohorts. The effect of NSUN2 on the biological characteristics of PCa cells was investigated on the basis of gain‐offunction and loss‐of‐function analyses. Subcutaneous models further uncovered the role of NSUN2 in tumor growth. Epi‐transcriptome assays with RNA bisulfite sequencing (RNA‐BisSeq) analysis and in vitro enzyme reaction assays were performed to validate the targeted effect of NSUN2 on AR. AR‐binding sites in the NSUN2 promoter were investigated by ChIP and luciferase assays to uncover the interplay between NSUN2 and AR signaling. RIP‐qPCR and EMSA methods were performed to confirm that YBX1 binds to AR m5C sites. Results NSUN2 is highly expressed in PCa and predicts poor outcome. NSUN2 plays roles as a PCa oncogene both in vitro and in vivo. Depletion of NSUN2 results in decreased expression and activities of AR, including AR‐V7. Mechanistically, NSUN2 posttranscriptionally stabilized AR by cluster m5C modification in a m5CYBX1‐dependent manner. Strikingly, treatment with enzalutamide, an effective AR inhibitor, reduces NSUN2 expression and decreases the m5C modification level in prostate cancer cells. Finally, we found that AR transcriptionally regulates NSUN2. Conclusion NSUN2 stabilizes AR mRNA through cluster 5‐methylcytosine modification and activates a positive feedback loop to promote prostate cancer.

Peritoneal implantation and lymph node metastasis have different driving mechanisms in ovarian cancer. Elucidating the underlying mechanism of lymph node metastasis is important for treatment outcomes. A new cell line, FDOVL, was established from a metastatic lymph node of a patient with primary platinum-resistant ovarian cancer and was then characterized. The effect of NOTCH1-p.C702fs mutation and NOTCH1 inhibitor on migration was evaluated in vitro and in vivo. Ten paired primary sites and metastatic lymph nodes were analyzed by RNA sequencing. The FDOVL cell line with serious karyotype abnormalities could be stably passaged and could be used to generated xenografts. NOTCH1-p.C702fs mutation was found exclusively in the FDOVL cell line and the metastatic lymph node. The mutation promoted migration and invasion in cell and animal models, and these effects were markedly repressed by the NOTCH inhibitor LY3039478. RNA sequencing confirmed CSF3 as the downstream effector of NOTCH1 mutation. Furthermore, the mutation was significantly more common in metastatic lymph nodes than in other peritoneal metastases in 10 paired samples (60% vs. 20%). The study revealed that NOTCH1 mutation is probably a driver of lymph node metastasis in ovarian cancer, which offers new ideas for the treatment of ovarian cancer lymph node metastasis with NOTCH inhibitors.

Abstract Tumor‐associated macrophages (TAMs) are one of the most abundant cell types in colorectal cancer (CRC) tumor microenvironment (TME). Recent studies observed complicated “cross‐talks” between cancer cells and macrophages in TME. However, the underlying mechanisms are still poorly elucidated. Here, PD‐L1 levels are very low in CRC cells but highly abundant in TAMs, and a specific PD‐L1+CD206+ macrophage subpopulation are identified, which is induced by tumor cells and associated with a poor prognosis. Mechanistic investigations reveal that CRC cells can secrete small extracellular vesicles (sEVs) taken up by macrophages that induce M2 like polarization and PD‐L1 expression, resulting in increased PD‐L1+CD206+ macrophage abundance and decreased T cell activity in CRC TME. sEV‐derived miR‐21‐5p and miR‐200a are identified as key signaling molecules mediating the regulatory effects of CRC on macrophages. Further studies reveal that CRC‐derived miR‐21‐5p and miR‐200a synergistically induces macrophage M2 like polarization and PD‐L1 expression by regulating the PTEN/AKT and SCOS1/STAT1 pathways, resulting in decreased CD8+ T cell activity and increased tumor growth. This study suggests that inhibiting the secretion of specific sEV‐miRNAs from CRC and targeting PD‐L1 in TAMs may serve as novel methods for CRC treatment as well as a sensitization method for anti‐PD‐L1 therapy in CRC.

Extracellular vesicles (EVs) are complex ecosystems that can be derived from all body cells and circulated in the body fluids. Characterizing the tissue-cellular source contributing to circulating EVs provides biological information about the cell or tissue of origin and their functional states. However, the relative proportion of tissue-cellular origin of circulating EVs in body fluid has not been thoroughly characterized. Here, we developed an approach for digital EVs quantification, called EV-origin, that enables enumerating of EVs tissue-cellular source contribution from plasma extracellular vesicles long RNA sequencing profiles. EV-origin was constructed by the input matrix of gene expression signatures and robust deconvolution algorithm, collectively used to separate the relative proportions of each tissue or cell type of interest. EV-origin respectively predicted the relative enrichment of seven types of hemopoietic cells and sixteen solid tissue subsets from exLR-seq profile. Using the EV-origin approach, we depicted an integrated landscape of the traceability system of plasma EVs for healthy individuals. We also compared the heterogenous tissue-cellular source components from plasma EVs samples with diverse disease status. Notably, the aberrant liver fraction could reflect the development and progression of hepatic disease. The liver fraction could also serve as a diagnostic indicator and effectively separate HCC patients from normal individuals. The EV-origin provides an approach to decipher the complex heterogeneity of tissue-cellular origin in circulating EVs. Our approach could inform the development of exLR-based applications for liquid biopsy.

Abstract Background Hypoxic tumors are refractory to DNA damage drugs. However, the underlying mechanism has yet to be elucidated. We aimed to identify lncRNAs that upregulated under hypoxia and their effects on colorectal cancer (CRC). Methods CRC cells were treated with 1% O2 to identify lncRNAs that upregulated under hypoxia. We integrated these lncRNAs with RNA-seq of 4 paired CRC tissues and TCGA data to get candidate lncRNAs. Multiple in vitro and in vivo assays were used to explore the role of LUCAT1 in CRC. Results We identified a hypoxia-induced lncRNA LUCAT1 that facilitated the growth of CRC cells and contributed to drug resistance of CRC cells both in vitro and in vivo. Mechanically, LUCAT1 interacts with polypyrimidine tract binding protein 1 (PTBP1) in CRC cells, facilitates the association of a set of DNA damage related genes with PTBP1, thus resulting in altered alternative splicing of these genes. Moreover, ectopic expression of PTBP1 in CRC cells with knockdown of LUCAT1 abrogated the effects induced by LUCAT1 knockdown. Chemotherapeutics drug combined with LUCAT1 knockdown via antisense oligonucleotides (ASO) would get a better outcome in vivo, compared with group treated with chemotherapeutic drug only. Notably, LUCAT1 is upregulated in CRC tissues, compared to adjacent normal tissues; and CRC patients with higher LUCAT1 have a worse prognosis and poorly responded to chemotherapy in the clinic. Conclusions Our data suggested CRC cells utilizes LUCAT1 to develop resistance to DNA damage drugs, and disrupting the LUCAT1/PTBP1 axis might be a promising therapeutic strategy for refractory hypoxic tumors.

BackgroundGastric cancer (GC) is a highly molecular heterogeneous tumor with poor prognosis. Epithelial-mesenchymal transition (EMT) process and cancer stem cells (CSCs) are reported to share common signaling pathways and cause poor prognosis in GC. Considering about the close relationship between these two processes, we aimed to establish a gene signature based on both processes to achieve better prognostic prediction in GC.MethodsThe gene signature was constructed by univariate Cox and the least absolute shrinkage and selection operator (LASSO) Cox regression analyses by using The Cancer Genome Atlas (TCGA) GC cohort. We performed enrichment analyses to explore the potential mechanisms of the gene signature. Kaplan-Meier analysis and time-dependent receiver operating characteristic (ROC) curves were implemented to assess its prognostic value in TCGA cohort. The prognostic value of gene signature on overall survival (OS), disease-free survival (DFS), and drug sensitivity was validated in different cohorts. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) validation of the prognostic value of gene signature for OS and DFS prediction was performed in the Fudan cohort.ResultsA prognostic signature including SERPINE1, EDIL3, RGS4, and MATN3 (SERM signature) was constructed to predict OS, DFS, and drug sensitivity in GC. Enrichment analyses illustrated that the gene signature has tight connection with the CSC and EMT processes in GC. Patients were divided into two groups based on the risk score obtained from the formula. The Kaplan-Meier analyses indicated high-risk group yielded significantly poor prognosis compared with low-risk group. Pearson’s correlation analysis indicated that the risk score was positively correlated with carboplatin and 5-fluorouracil IC50 of GC cell lines. Multivariate Cox regression analyses showed that the gene signature was an independent prognostic factor for predicting GC patients’ OS, DFS, and susceptibility to adjuvant chemotherapy.ConclusionsOur SERM prognostic signature is of great value for OS, DFS, and drug sensitivity prediction in GC, which may give guidance to the development of targeted therapy for CSC- and EMT-related gene in the future.

Abstract Background The majority of breast cancer patients die of metastasis rather than primary tumors, whereas the molecular mechanisms orchestrating cancer metastasis remains poorly understood. Long noncoding RNAs (lncRNA) have been shown to regulate cancer occurrence and progression. However, the lncRNAs that drive metastasis in cancer patients and their underlying mechanisms are still largely unknown. Methods lncRNAs highly expressed in metastatic lymph nodes were identified by microarray. Survival analysis were made by Kaplan-Meier method. Cell proliferation, migration, and invasion assay was performed to confirm the phenotype of LINC02273. Tail vein model and mammary fat pad model were used for in vivo study. RNA pull-down and RIP assay were used to confirm the interaction of hnRNPL and LINC02273. Chromatin isolation by RNA purification followed by sequencing (ChIRP-seq), RNA-seq, ChIP-seq, and luciferase reporter assay reveal hnRNPL-LINC02273 regulates AGR2. Antisense oligonucleotides were used for in vivo treatment. Results We identified a novel long noncoding RNA LINC02273, whose expression was significantly elevated in metastatic lesions compared to the primary tumors, by genetic screen of matched tumor samples. Increased LINC02273 promoted breast cancer metastasis in vitro and in vivo. We further showed that LINC02273 was stabilized by hnRNPL, a protein increased in metastatic lesions, in breast cancer cells. Mechanistically, hnRNPL-LINC02273 formed a complex which activated AGR2 transcription and promoted cancer metastasis. The recruitment of hnRNPL-LINC02273 complex to AGR2 promoter region epigenetically upregulated AGR2 by augmenting local H3K4me3 and H3K27ac levels. Combination of AGR2 and LINC02273 was an independent prognostic factor for predicting breast cancer patient survival. Moreover, our data revealed that LINC02273-targeting antisense oligonucleotides (ASO) substantially inhibited breast cancer metastasis in vivo. Conclusions Our findings uncover a key role of LINC02273-hnRNPL-AGR2 axis in breast cancer metastasis and provide potential novel therapeutic targets for metastatic breast cancer intervention.

Abstract Circular RNAs (circRNAs) are an intriguing class of widely prevalent endogenous RNAs, the vast majority of which have not been characterized functionally. Here, we identified a novel oncogenic circRNA originating from the back‐splicing of Exon2 and Exon3 of a tumor suppressor gene, ARHGAP35 (also known as P190‐A), termed as circARHGAP35. have observe that circARHGAP35 and linear ARHGAP35 have antithetical expression and functions. Interestingly, circARHGAP35 contains a 3867 nt long ORF with an m6A‐modified start codon and encodes a truncated protein comprising four FF domains and lacking the Rho GAP domain. Mechanistically, circARHGAP35 protein promotes cancer cell progression by interacting with TFII‐I protein in the nucleus. The RNA binding protein, HNRNPL, facilitates the formation of circARHGAP35. Clinically, circARHGAP35 is associated with poor survival in cancer patients. Our findings characterize an oncogenic circRNA and demonstrate a novel mechanism of oncogene activation in cancer by circRNA through the production of a truncated protein.

RNA splicing may generate different kinds of splice junctions, such as linear, back-splice and fusion junctions. Only a limited number of programs are available for detection and quantification of splice junctions. Here, we present Assembling Splice Junctions Analysis (ASJA), a software package that identifies and characterizes all splice junctions from high-throughput RNA sequencing (RNA-seq) data. ASJA processes assembled transcripts and chimeric alignments from the STAR aligner and StringTie assembler. ASJA provides the unique position and normalized expression level of each junction. Annotations and integrative analysis of the junctions enable additional filtering. It is also appropriate for the identification of novel junctions. ASJA is available at https://github.com/HuangLab-Fudan/ASJA. Keywords: RNA splicing, Splice junctions, RNA-seq, Circular RNA

Background: The prevalence of HER2 alterations in pan-cancer indicates a broader range of application of HER2-targeted therapies; however, biomarkers for such therapies are still insufficient and limited to breast cancer and gastric cancer. Methods: Using multi-omics data from The Cancer Genome Atlas (TCGA), the landscape of HER2 alterations was exhibited across 33 tumor types. A HER2 index was constructed using one-class logistic regression (OCLR). With the predictive value validated in GEO cohorts and pan-cancer cell lines, the index was then applied to evaluate the HER2-enriched expression pattern across TCGA pan-cancer types. Findings: Increased HER2 somatic copy number alterations (SCNAs) could be divided into two patterns, focal- or arm-level. The expression-based HER2 index successfully distinguished the HER2-enriched subtype from the others and provided a stable and superior performance in predicting the response to HER2-targeted therapies both in breast tumor tissue and pan-cancer cell lines. With frequencies varying from 12.0% to 0.9%, tumors including head and neck squamous tumors, gastrointestinal tumors, bladder cancer, lung cancer and uterine tumors exhibited high HER2 indices together with HER2 amplification or overexpression, which may be more suitable for HER2-targeted therapies. The BLCA.3 and HNSC.Basal were the most distinguishable subtypes within bladder cancer and head and neck cancer respectively by HER2 index, implying their potential benefits from HER2-targeted therapies. Interpretation: As a pan-cancer predictive biomarker of HER2-targeted therapies, the HER2 index could help identify potential candidates for such treatment in multiple tumor types by combining with HER2 multi-omics features. The discoveries of our study highlight the importance of incorporating transcriptional pattern into the assessment of HER2 status for better patient selection. Funding: The National Key Research and Development Program of China; Clinical Research and Cultivation Project of Shanghai ShenKang Hospital Development Center.

Summary: Hepatocellular carcinoma (HCC) initiation is characterized by stepwise accumulation of molecular alterations, during which the early events are largely unknown. Here, we presented a comprehensive genomic and transcriptomic landscape at stages of hepatitis, cirrhosis, and HCC by using a diethylnitrosamine-induced rat HCC model. We observed the early occurrence of gene instability and aberrant cancer associated signaling pathways in liver hepatitis. We further characterized the progressive molecular changes during hepatocarcinogenesis, wherein the intense rivalry between tumor-suppressive and oncogenic strengths occurred in cirrhosis stage. Despite the significant pathological difference, mutation signatures and expression landscape are highly similar between hepatitis and cirrhosis stages. Furthermore, we identified PI3K-Akt signaling pathway as a key pathway in the process of hepatocarcinogenesis through integrative analysis, and PIK3CD is a potential biomarker indicating HCC recurrence. The dynamic immune response during hepatocarcinogenesis, such as continuous decline of monocytes, suggests an immunological intervention strategy beyond chemoprevention for liver cancer.

African swine fever virus (ASFV) is the causative agent of the epidemic of African swine fever (ASF), with virulent strains having a mortality rate of up to 100% and presenting devastating impacts on animal farming. Since ASF was first reported in China in 2018, ASFV still exists and poses a potential threat to the current pig industry. Low-virulence and genotype I strains of ASFV have been reported in China, and the prevention and control of ASF is more complicated. Insufficient understanding of the interaction of ASFV with the host immune system hinders vaccine development. Physical barriers, nonspecific immune response and acquired immunity are the three barriers of the host against infection. To escape the innate immune response, ASFV invades monocytes/macrophages and dendritic cells, thereby inhibiting IFN expression, regulating cytokine expression and the body’s inflammatory response process. Meanwhile, in order to evade the adaptive immune response, ASFV inhibits antigen presentation, induces the production of non-neutralizing antibodies, and inhibits apoptosis. Recently, significant advances have been achieved in vaccine development around the world. Live attenuated vaccines (LAVs) based on artificially deleting specific virulence genes can achieve 100% homologous protection and partial heterologous protection. The key of subunit vaccines is identifying the combination of antigens that can effectively provide protection and selecting carriers that can effectively deliver the antigens. In this review, we introduce the epidemic trend of ASF and the impact on the pig industry, analyze the interaction mechanism between ASFV and the body’s immune system, and compare the current status of potential vaccines in order to provide a reference for the development of effective ASF vaccines.

Long intergenic non-coding RNAs have been linked to cancer development. Here the authors, using RNA-seq and genomic amplification data, identify lincRNAs deregulated in hepatocellular carcinoma and propose that Linc01138 is stabilized by IGF2BP1/3 in the cytoplasm, and binds and stabilizes the methyltransferase PRMT5 by preventing the association of PRMT5 to the E3 ubiquitin ligase CHIP.

Summary: The diversity and complexity of the cancer transcriptome may contain transcripts unique to the tumor environment. Here, we report a LIN28B variant, LIN28B-TST, which is specifically expressed in hepatocellular carcinoma (HCC) and many other cancer types. Expression of LIN28B-TST is associated with significantly poor prognosis in HCC patients. LIN28B-TST initiates from a de novo alternative transcription initiation site that harbors a strong promoter regulated by NFYA but not c-Myc. Demethylation of the LIN28B-TST promoter might be a prerequisite for its transcription and transcriptional regulation. LIN28B-TST encodes a protein isoform with additional N-terminal amino acids and is critical for cancer cell proliferation and tumorigenesis. Our findings reveal a mechanism of LIN28B activation in cancer and the potential utility of LIN28B-TST for clinical purposes.

Circular RNAs are formed from exon back-splicing, the significance of these endogenous RNAs is beginning to be unraveled. Here, the authors identify thousands of circular RNAs differentially expressed between normal and cancer tissues and show that an abundant circular RNA generated from HIPK3regulates cell growth.

Tumor cells increase their glucose consumption through aerobic glycolysis to manufacture the necessary biomass required for proliferation, commonly known as the Warburg effect. Accumulating evidences suggest that microRNAs (miRNAs) interact with their target genes and contribute to metabolic reprogramming in cancer cells. By integrating high-throughput screening data and the existing miRNA expression datasets, we explored the roles of candidate glycometabolism-regulating miRNAs in gastric cancer (GC). Subsequent investigation of the characterized miRNAs indicated that miR-129-5p inhibits glucose metabolism in GC cells. miRNA-129-5p directly targets the 3′-UTR of SLC2A3, thereby suppressing glucose consumption, lactate production, cellular ATP levels, and glucose uptake of GC cells. In addition, the PI3K-Akt and MAPK signaling pathways are involved in the effects of the miR-129-5p/SLC2A3 axis, regulating GC glucose metabolism and growth. These results reveal a novel role of the miR-129-5p/SLC2A3 axis in reprogramming the glycometabolism process in GC cells and indicate a potential therapeutic target for the treatment of this disease.

Although studies provide insights into the neurobiology of stress and depression, the exact molecular mechanisms underlying their pathologies remain largely unknown. Long non-coding RNA (lncRNA) has been implicated in brain functions and behavior. A potential link between lncRNA and psychiatric disorders has been proposed. However, it remains undetermined whether IncRNA regulation, in the brain, contributes to stress or depression pathologies. In this study, we used a valid animal model of depression-like symptoms; namely learned helplessness, RNA-seq, Gene Ontology and co-expression network analyses to profile the expression pattern of lncRNA and mRNA in the hippocampus of mice. We identified 6346 differentially expressed transcripts. Among them, 340 lncRNAs and 3559 protein coding mRNAs were differentially expressed in helpless mice in comparison with control and/or non-helpless mice (inescapable stress resilient mice). Gene Ontology and pathway enrichment analyses indicated that induction of helplessness altered expression of mRNAs enriched in fundamental biological functions implicated in stress/depression neurobiology such as synaptic, metabolic, cell survival and proliferation, developmental and chromatin modification functions. To explore the possible regulatory roles of the altered lncRNAs, we constructed co-expression networks composed of the lncRNAs and mRNAs. Among our differentially expressed lncRNAs, 17% showed significant correlation with genes. Functional co-expression analysis linked the identified lncRNAs to several cellular mechanisms implicated in stress/depression neurobiology. Importantly, 57% of the identified regulatory lncRNAs significantly correlated with 18 different synapse-related functions. Thus, the current study identifies for the first time distinct groups of lncRNAs regulated by induction of learned helplessness in the mouse brain. Our results suggest that lncRNA-directed regulatory mechanisms might contribute to stress-induced pathologies; in particular, to inescapable stress-induced synaptic modifications.

Protocol Summary Here, we describe targeted gene deletion and replacement in human cells via the CRISPR/Cas9 system using two guide RNAs. The system effectively generated targeted deletions of varied length, regardless of the transcriptional status of the target gene. It is notable that targeted gene deletions generated via CRISPR/Cas9 and two guide RNAs resulted in the formation of correct junctions at high efficiency. Moreover, in the presence of a homology repair donor, the CRISPR/Cas9 system could guide precise gene replacement. Our results illustrate that the CRISPR/Cas9 system can be used to precisely and effectively generate targeted deletions or gene replacement in human cells, which will facilitate characterization of functional domains in protein-coding genes as well as noncoding regulatory sequences in animal genomes.

The prokaryotic type II CRISPR/Cas9 system has been adapted to perform targeted genome editing in cells and model organisms. Here, we describe targeted gene deletion and replacement in human cells via the CRISPR/Cas9 system using two guide RNAs. The system effectively generated targeted deletions of varied length, regardless of the transcriptional status of the target gene. It is notable that targeted gene deletions generated via CRISPR/Cas9 and two guide RNAs resulted in the formation of correct junctions at high efficiency. Moreover, in the presence of a homology repair donor, the CRISPR/Cas9 system could guide precise gene replacement. Our results illustrate that the CRISPR/Cas9 system can be used to precisely and effectively generate targeted deletions or gene replacement in human cells, which will facilitate characterization of functional domains in protein-coding genes as well as noncoding regulatory sequences in animal genomes.

MicroRNAs (miRNAs) are small, single-stranded, non-coding RNAs that play pivotal roles in human cancer development and progression, such as tumor metastasis. Here, we identified the miRNAs that regulate hepatocellular carcinoma (HCC) cell migration by a high-throughput screening method using the classical wound-healing assay with time-lapse video microscopy and validation with a transwell migration assay. Eleven miRNAs (miR-134, -146b-3p, -188-3p, -525-3p, -661, -767-5p, -891a, -891b, -1244, -1247 and miR-1471) were found to promote or inhibit HCC cell migration. Further investigation revealed that miR-134 suppressed the invasion and metastasis of HCC cells in vitro and in vivo, and integrin beta 1 (ITGB1) was a direct and functional target gene of miR-134. Moreover, miR-134 inhibited the phosphorylation of focal adhesion kinase (FAK) and the activation of RhoA downstream of the ITGB1 pathway, thereby decreasing stress fiber formation and cell adhesion in HCC cells. In conclusion, we demonstrated that miR-134 is a novel metastasis suppressor in HCC and could be a potential therapeutic target for the treatment of HCC.

MicroRNAs (miRNAs) are a class of small, non-coding RNA molecules that are often found at chromosomal breakpoints and play a vital role in human cancer. Our previous study found that miR-550a, a frequently amplified miRNA on 7p14.3, was upregulated in hepatocellular carcinoma (HCC). However, the possible functions and molecular mechanisms of miR-550a in HCC remain unknown. In this study, gain-of-function and loss-of-function assays revealed that miR-550a markedly promoted HCC cell migration and invasion. In addition, we discovered that cytoplasmic polyadenylation element binding protein 4 (CPEB4) was a potential target of miR-550a in HCC. Further analyses showed that knockdown of CPEB4 expression significantly facilitated HCC cell migration and invasion, which phenocopied the effects of miR-550a on HCC cells. Moreover, a decrease in CPEB4 expression mediated miR-550a-induced liver cancer cell migration and invasion. Interestingly, CPEB4 is frequently downregulated in HCC, and its expression levels correlate with the overall survival of HCC patients. Together, these results suggested that this newly identified miR-550a-CPEB4 axis may be involved in HCC cell metastasis. Moreover, the expression levels of CPEB4 could be used to predict outcomes in HCC patients. Our findings provide novel potential targets for HCC therapy and prognosis.