Identification of novel prognostic alternative splicing signature in papillary renal cell carcinoma

Papillary renal cell carcinoma (pRCC) is a heterogeneous disease containing multifocal or solitary tumors with an aggressive phenotype. Increasing evidence has indicated the involvement of aberrant splicing variants in renal cell cancer, while systematic profiling of aberrant alternative splicing (AS) in pRCC was lacking and largely unknown. In the current study, comprehensive profiling of AS events were performed based on the integration of pRCC cohort from the Cancer Genome Atlas database and SpliceSeq software. With rigorous screening and univariate Cox analysis, a total of 2077 prognoses AS events from 1642 parent genes were identified. Then, stepwise least absolute shrinkage and selection operator method-penalized Cox regression analyses with 10-fold cross-validation followed by multivariate Cox regression were used to construct the prognostic AS signatures within each AS type. And a final 21 AS event-based signature was proposed which showed potent prognostic capability in stratifying patients into low- and high-risk subgroups (P < .0001). Furthermore, time-dependent receiver operating characteristics curves confirmed that the final AS signature was effective and robust in predicting overall survival for pRCC patients with the area under the curve above 0.9 from 1 to 5 years. In addition, splicing correlation network was built to uncover the potential regulatory pattern among prognostic splicing factors and candidate AS events. Besides, gene set enrichment analysis revealed the involvement of these candidates AS events in tumor-related pathways including extracellular matrix organization, oxidative phosphorylation, and P53 signaling pathways. Taken together, our results could contribute to elucidating the underlying mechanism of AS in the oncogenesis process and broaden the novel field of prognostic and clinical application of molecule-targeted approaches in pRCC.     

A novel DNA damage and repair‐related gene signature to improve predictive capacity of overall survival for patients with gliomas

Gliomas, as the most lethal and malignant brain tumours in adults, remain a major challenge worldwide. DNA damage and repair‐related genes (DDRRGs) appear to play a significant role in gliomas, but the studies of DDRRGs are still insufficient. Herein, we systematically explored and analysed 1547 DDRRGs in 938 glioma samples from TCGA and CGGA datasets. Using least absolute shrinkage and selection operator (LASSO) Cox regression analysis, we identified a 16‐DDRRG signature, characterized by high‐risk and low‐risk patterns. This risk model harbours robust predictive capability for overall survival of glioma patients. We found the high‐risk score is strongly associated with well‐known malignant features of gliomas, such as the mesenchymal subtype, IDH‐wildtype, 1p/19q non‐codeletion and MGMT promoter unmethylated status. In addition, we found that the high‐risk score is also linked with multiple oncogenic pathways and therapeutic resistance. Significantly, we found the high‐risk group has higher enrichment of immunosuppressive cells (M2‐type macrophages, Tregs and MDSCs) and immune inhibition biomarkers (PD‐1, PD‐L1 and CTLA‐4). Lastly, we proved that SMC4, which has the highest positive regression coefficient in our risk model, is strongly linked with malignant progression and TMZ resistance of gliomas in a E2F1‐dependent manner.     

Genome‐wide methylomic analyses identify prognostic epigenetic signature in lower grade glioma

Glioma is the most malignant and aggressive type of brain tumour with high heterogeneity and mortality. Although some clinicopathological factors have been identified as prognostic biomarkers, the individual variants and risk stratification in patients with lower grade glioma (LGG) have not been fully elucidated. The primary aim of this study was to identify an efficient DNA methylation combination biomarker for risk stratification and prognosis in LGG. We conducted a retrospective cohort study by analysing whole genome DNA methylation data of 646 patients with LGG from the TCGA and GEO database. Cox proportional hazard analysis was carried out to screen and construct biomarker model that predicted overall survival (OS). The Kaplan‐Meier survival curves and time‐dependent ROC were constructed to prove the efficiency of the signature. Then, another independent cohort was used to further validate the finding. A two‐CpG site DNA methylation signature was identified by multivariate Cox proportional hazard analysis. Further analysis indicated that the signature was an independent survival predictor from other clinical factors and exhibited higher predictive accuracy compared with known biomarkers. This signature was significantly correlated with immune‐checkpoint blockade, immunotherapy‐related signatures and ferroptosis regulator genes. The expression pattern and functional analysis showed that these two genes corresponding with two methylation sites contained in the model were correlated with immune infiltration level, and involved in MAPK and Rap1 signalling pathway. The signature may contribute to improve the risk stratification of patients and provide a more accurate assessment for precision medicine in the clinic.     

Efficacy and safety of combination PD‐1/PD‐L1 checkpoint inhibitors for malignant solid tumours: A systematic review

Treatment of multiple malignant solid tumours with programmed death (PD)‐1/PD ligand (PD‐L) 1 inhibitors has been reported. However, the efficacy and immune adverse effects of combination therapies are controversial. This meta‐analysis was performed with PubMed, Web of Science, Medline, EMBASE and Cochrane Library from their inception until January 2020. Random‐effect model was adopted because of relatively high heterogeneity. We also calculated hazard ratio (HR) of progression‐free survival (PFS), overall survival (OS) and risk ratio (RR) of adverse events (AEs), the incidence of grade 3‐5 AEs by tumour subgroup, therapeutic schedules and therapy lines. Nineteen articles were selected using the search strategy for meta‐analysis. Combined PD‐1/PD‐L1 inhibitors prolonged OS and PFS (HR 0.72, P < 0.001) and (HR 0.66, P < 0.001). In addition, incidence of all‐grade and grade 3‐5 AEs was not significant in the two subgroup analyses (HR 1.01, P = 0.31) and (HR 1.10, P = 0.07), respectively. Our meta‐analysis indicated that combination therapy with PD‐1/PD‐L1 inhibitors had greater clinical benefits and adverse events were not increased significantly.     

The Proneural Molecular Signature Is Enriched in Oligodendrogliomas and Predicts Improved Survival among Diffuse Gliomas

The Cancer Genome Atlas Project (TCGA) has produced an extensive collection of ‘-omic’ data on glioblastoma (GBM), resulting in several key insights on expression signatures. Despite the richness of TCGA GBM data, the absence of lower grade gliomas in this data set prevents analysis genes related to progression and the uncovering of predictive signatures. A complementary dataset exists in the form of the NCI Repository for Molecular Brain Neoplasia Data (Rembrandt), which contains molecular and clinical data for diffuse gliomas across the full spectrum of histologic class and grade. Here we present an investigation of the significance of the TCGA consortium''s expression classification when applied to Rembrandt gliomas. We demonstrate that the proneural signature predicts improved clinical outcome among 176 Rembrandt gliomas that includes all histologies and grades, including GBMs (log rank test p = 1.16e-6), but also among 75 grade II and grade III samples (p = 2.65e-4). This gene expression signature was enriched in tumors with oligodendroglioma histology and also predicted improved survival in this tumor type (n = 43, p = 1.25e-4). Thus, expression signatures identified in the TCGA analysis of GBMs also have intrinsic prognostic value for lower grade oligodendrogliomas, and likely represent important differences in tumor biology with implications for treatment and therapy. Integrated DNA and RNA analysis of low-grade and high-grade proneural gliomas identified increased expression and gene amplification of several genes including GLIS3, TGFB2, TNC, AURKA, and VEGFA in proneural GBMs, with corresponding loss of DLL3 and HEY2. Pathway analysis highlights the importance of the Notch and Hedgehog pathways in the proneural subtype. This demonstrates that the expression signatures identified in the TCGA analysis of GBMs also have intrinsic prognostic value for low-grade oligodendrogliomas, and likely represent important differences in tumor biology with implications for treatment and therapy.     

Alternative splicing coupled mRNA decay shapes the temperature‐dependent transcriptome

Mammalian body temperature oscillates with the time of the day and is altered in diverse pathological conditions. We recently identified a body temperature‐sensitive thermometer‐like kinase, which alters SR protein phosphorylation and thereby globally controls alternative splicing (AS). AS can generate unproductive variants which are recognized and degraded by diverse mRNA decay pathways—including nonsense‐mediated decay (NMD). Here we show extensive coupling of body temperature‐controlled AS to mRNA decay, leading to global control of temperature‐dependent gene expression (GE). Temperature‐controlled, decay‐inducing splicing events are evolutionarily conserved and pervasively found within RNA‐binding proteins, including most SR proteins. AS‐coupled poison exon inclusion is essential for rhythmic GE of SR proteins and has a global role in establishing temperature‐dependent rhythmic GE profiles, both in mammals under circadian body temperature cycles and in plants in response to ambient temperature changes. Together, these data identify body temperature‐driven AS‐coupled mRNA decay as an evolutionary ancient, core clock‐independent mechanism to generate rhythmic GE.     

Prognostic significance of gamma‐glutamyl transpeptidase to albumin ratio in patients with intrahepatic cholangiocarcinoma after hepatectomy

Inflammation has been reported to play an important role in tumour progression and prognosis. In this study, we evaluated the prognostic significance of γ‐glutamyl transpeptidase (GGT) to albumin ratio (GAR) in patients with intrahepatic cholangiocarcinoma (ICC) after hepatectomy. We retrospectively analysed 650 ICC patients underwent hepatectomy at three Chinese medical centres between January 2009 and September 2017. Patients were classified into derivation cohort (n = 509) and validation cohort (n = 141). Receiver operating characteristic (ROC) curve was used to determine the optimal cut‐off value for GAR. Survival curve and cox regression analysis were applied to assess the prognostic power of GAR. The prognostic accuracy of GAR was compared with other variables by ROC curve. The optimal cut‐off value for GAR was 1.3655. Preoperative high GAR was closely related to tumour number, lymph node invasion and GGT. The survival curve of derivation and validation cohorts showed that patients in the high GAR group had significantly shorter overall survival (OS) and disease‐free survival (DFS) than patients in the low GAR group. Multivariate analysis in the derivation cohort confirmed that GAR was an independent prognostic factor for survival outcomes. Moreover, the ROC curve revealed that GAR had better predictive accuracy than other variables. High GAR predicted poor OS and DFS in ICC patients after hepatectomy. GAR may be a novel, simple and effective prognostic marker for ICC patients.     

Classifying behavior from short‐interval biologging data: An example with GPS tracking of birds

1. Recent advances in digital data collection have spurred accumulation of immense quantities of data that have potential to lead to remarkable ecological insight, but that also present analytic challenges. In the case of biologging data from birds, common analytical approaches to classifying movement behaviors are largely inappropriate for these massive data sets.
2. We apply a framework for using K‐means clustering to classify bird behavior using points from short time interval GPS tracks. K‐means clustering is a well‐known and computationally efficient statistical tool that has been used in animal movement studies primarily for clustering segments of consecutive points. To illustrate the utility of our approach, we apply K‐means clustering to six focal variables derived from GPS data collected at 1–11 s intervals from free‐flying bald eagles (Haliaeetus leucocephalus) throughout the state of Iowa, USA. We illustrate how these data can be used to identify behaviors and life‐stage‐ and age‐related variation in behavior.
3. After filtering for data quality, the K‐means algorithm identified four clusters in >2 million GPS telemetry data points. These four clusters corresponded to three movement states: ascending, flapping, and gliding flight; and one non‐moving state: perching. Mapping these states illustrated how they corresponded tightly to expectations derived from natural history observations; for example, long periods of ascending flight were often followed by long gliding descents, birds alternated between flapping and gliding flight.
4. The K‐means clustering approach we applied is both an efficient and effective mechanism to classify and interpret short‐interval biologging data to understand movement behaviors. Furthermore, because it can apply to an abundance of very short, irregular, and high‐dimensional movement data, it provides insight into small‐scale variation in behavior that would not be possible with many other analytical approaches.

     

Identification of survival-related alternative splicing signatures in acute myeloid leukemia

Aberrant RNA alternative splicing (AS) variants play critical roles in tumorigenesis and prognosis in human cancers. Here, we conducted a comprehensive profiling of aberrant AS events in acute myeloid leukemia (AML). RNA AS profile, including seven AS types, and the percent spliced in (PSI) value for each patient were generated by SpliceSeq using RNA-seq data from TCGA. Univariate followed by multivariate Cox regression analysis were used to identify survival-related AS events and develop the AS signatures. A nomogram was developed, and its predictive efficacy was assessed. About 27,892 AS events and 3,178 events were associated with overall survival (OS) after strict filtering. Parent genes of survival-associated AS events were mainly enriched in leukemia-associated processes including chromatin modification, autophagy, and T-cell receptor signaling pathway. The 10 AS signature based on seven types of AS events showed better efficacy in predicting OS of patients than those built on a single AS event type. The area under curve (AUC) value of the 10 AS signature for 3-year OS was 0.91. Gene set enrichment analysis (GSEA) confirmed that these survival-related AS events contribute to AML progression. Moreover, the nomogram showed good predictive performance for patient''s prognosis. Finally, the correlation network of AS variants with splicing factor genes found potential important regulatory genes in AML. The present study presented a systematic analysis of survival-related AS events and developed AS signatures for predicting the patient’s survival. Further studies are needed to validate the signatures in independent AML cohorts and might provide a promising perspective for developing therapeutic targets.     

Prognostic power of a lipid metabolism gene panel for diffuse gliomas

Lipid metabolism reprogramming plays important role in cell growth, proliferation, angiogenesis and invasion in cancers. However, the diverse lipid metabolism programmes and prognostic value during glioma progression remain unclear. Here, the lipid metabolism‐related genes were profiled using RNA sequencing data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) database. Gene ontology (GO) and gene set enrichment analysis (GSEA) found that glioblastoma (GBM) mainly exhibited enrichment of glycosphingolipid metabolic progress, whereas lower grade gliomas (LGGs) showed enrichment of phosphatidylinositol metabolic progress. According to the differential genes of lipid metabolism between LGG and GBM, we developed a nine‐gene set using Cox proportional hazards model with elastic net penalty, and the CGGA cohort was used for validation data set. Survival analysis revealed that the obtained gene set could differentiate the outcome of low‐ and high‐risk patients in both cohorts. Meanwhile, multivariate Cox regression analysis indicated that this signature was a significantly independent prognostic factor in diffuse gliomas. Gene ontology and GSEA showed that high‐risk cases were associated with phenotypes of cell division and immune response. Collectively, our findings provided a new sight on lipid metabolism in diffuse gliomas.     

Pyroptosis‐related molecular classification and immune microenvironment infiltration in breast cancer: A novel therapeutic target

The underlying role of pyroptosis in breast cancer (BC) remains unknown. Herein, we investigated the correlations of 33 pyroptosis‐related genes (PRGs) with immune checkpoints and immune cell infiltrations in BC patients based on The Cancer Genome Atlas cohort (n = 996) and Gene Expression Omnibus cohort (n = 3,262). Enrichment analysis revealed that these PRGs mainly functioned in pyroptosis, inflammasomes and regulation of autophagy pathway. Four prognostic independent PRGs (CASP9, TIRAP, GSDMC and IL18) were identified. Then, cluster 1/2 was recognized using consensus clustering for these four PRGs. Patients from cluster 1 had a favourable prognosis and diverse immune cell infiltrations. A nomogram was developed based on age, TNM stage, tumour subtype and pyroptosis score. Patients with the high‐risk group exhibited worse 5‐year OS, and the result was consistent in the external cohort. Additionally, high‐risk group patients were associated with downregulated immune checkpoint expression. Further analysis suggested that the high‐risk group patients were associated with a higher IC50 of paclitaxel, doxorubicin, cisplatin, methotrexate and vinorelbine. In summarizing, the pyroptosis score‐based nomogram might serve as an independent prognostic predictor and could guide medication for chemotherapy. Additionally, it may bring novel insight into the regulation of tumour immune microenvironment in BC and help to achieve precision immunotherapy.     

SQLE facilitates the pancreatic cancer progression via the lncRNA‐TTN‐AS1/miR‐133b/SQLE axis

Studies have shown that SQLE is highly expressed in a variety of tumours and promotes tumour progression. However, the role of SQLE in pancreatic cancer (PC) has not been reported. Here, we aim to study the role and molecular mechanism of SQLE in PC. Immunohistochemistry and functional experiments showed that SQLE was highly expressed in PC tissues and promoted the proliferation and invasion of PC cells. Terbinafine, an inhibitor of SQLE, inhibited this effect. In order to further study the upstream mechanism that regulates SQLE, we used bioinformatics technology to lock miR‐133b and lncRNA‐TTN‐AS. In situ hybridization was used to detect the expression of miR‐133b and lncRNA‐TTN‐AS1 in PC tissues. The luciferase reporter gene experiment was used to confirm the binding of miR‐133b and lncRNA‐TTN‐AS1. The results showed that miR‐133b was down‐regulated in PC tissues and negatively correlated with the expression of SQLE. LncRNA‐TTN‐AS1 was upregulated in pancreatic cancer tissues and positively correlated with the expression of SQLE. Luciferase gene reporter gene analysis confirmed lncRNA‐TTN‐AS1 directly binded to miR‐133b. Therefore, we propose that targeting the lncRNA‐TTN‐AS1/miR‐133b/SQLE axis is expected to provide new ideas for the clinical treatment of PC patients.     

Platelet biomarkers for a descending cognitive function: A proteomic approach

Memory loss is the most common clinical sign in Alzheimer''s disease (AD); thus, searching for peripheral biomarkers to predict cognitive decline is promising for early diagnosis of AD. As platelets share similarities to neuron biology, it may serve as a peripheral matrix for biomarkers of neurological disorders. Here, we conducted a comprehensive and in‐depth platelet proteomic analysis using TMT‐LC‐MS/MS in the populations with mild cognitive impairment (MCI, MMSE = 18–23), severe cognitive impairments (AD, MMSE = 2–17), and the age‐/sex‐matched normal cognition controls (MMSE = 29–30). A total of 360 differential proteins were detected in MCI and AD patients compared with the controls. These differential proteins were involved in multiple KEGG pathways, including AD, AMP‐activated protein kinase (AMPK) pathway, telomerase RNA localization, platelet activation, and complement activation. By correlation analysis with MMSE score, three positively correlated pathways and two negatively correlated pathways were identified to be closely related to cognitive decline in MCI and AD patients. Partial least squares discriminant analysis (PLS‐DA) showed that changes of nine proteins, including PHB, UQCRH, CD63, GP1BA, FINC, RAP1A, ITPR1/2, and ADAM10 could effectively distinguish the cognitively impaired patients from the controls. Further machine learning analysis revealed that a combination of four decreased platelet proteins, that is, PHB, UQCRH, GP1BA, and FINC, was most promising for predicting cognitive decline in MCI and AD patients. Taken together, our data provide a set of platelet biomarkers for predicting cognitive decline which may be applied for the early screening of AD.     

Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma

There is growing evidence that alternative splicing (AS) plays an important role in cancer development. However, a comprehensive analysis of AS signatures in kidney renal clear cell carcinoma (KIRC) is lacking and urgently needed. It remains unclear whether AS acts as diagnostic biomarkers in predicting the prognosis of KIRC patients. In the work, gene expression and clinical data of KIRC were obtained from The Cancer Genome Atlas (TCGA), and profiles of AS events were downloaded from the SpliceSeq database. The RNA sequence/AS data and clinical information were integrated, and we conducted the Cox regression analysis to screen survival-related AS events and messenger RNAs (mRNAs). Correlation between prognostic AS events and gene expression were analyzed using the Pearson correlation coefficient. Protein-protein interaction analysis was conducted for the prognostic AS-related genes, and a potential regulatory network was built using Cytoscape (version 3.6.1). Meanwhile, functional enrichment analysis was conducted. A prognostic risk score model is then established based on seven hub genes (KRT222, LENG8, APOB, SLC3A1, SCD5, AQP1, and ADRA1A) that have high performance in the risk classification of KIRC patients. A total 46,415 AS events including 10,601 genes in 537 patients with KIRC were identified. In univariate Cox regression analysis, 13,362 survival associated AS events and 8,694 survival-specific mRNAs were detected. Common 3,105 genes were screen by overlapping 13,362 survival associated AS events and 8,694 survival-specific mRNAs. The Pearson correlation analysis suggested that 13 genes were significantly correlated with AS events (Pearson correlation coefficient >0.8 or <−0.8). Then, We conducted multivariate Cox regression analyses to select the potential prognostic AS genes. Seven genes were identified to be significantly related to OS. A prognostic model based on seven genes was constructed. The area under the ROC curve was 0.767. In the current study, a robust prognostic prediction model was constructed for KIRC patients, and the findings revealed that the AS events could act as potential prognostic biomarkers for KIRC.