Gene signature associated with neuro-endocrine activity predicting prognosis of pancreatic carcinoma

Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows

BACKGROUND

Heat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses.

RESULTS

A total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively.

CONCLUSIONS

Physiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response.

Identification of genomic biomarkers and their pathway crosstalks for deciphering mechanistic links in glioblastoma

Glioblastoma is a grade IV pernicious neoplasm occurring in the supratentorial region of brain. As its causes are largely unknown, it is essential to understand its dynamics at the molecular level. This necessitates the identification of better diagnostic and prognostic molecular candidates. Blood-based liquid biopsies are emerging as a novel tool for cancer biomarker discovery, guiding the treatment and improving its early detection based on their tumour origin. There exist previous studies focusing on the identification of tumour-based biomarkers for glioblastoma. However, these biomarkers inadequately represent the underlying pathological state and incompletely illustrate the tumour because of non-recursive nature of this approach to monitor the disease. Also, contrary to the tumour biopsies, liquid biopsies are non-invasive and can be performed at any interval during the disease span to surveil the disease. Therefore, in this study, a unique dataset of blood-based liquid biopsies obtained primarily from tumour-educated blood platelets (TEP) is utilised. This RNA-seq data from ArrayExpress is acquired comprising human cohort with 39 glioblastoma subjects and 43 healthy subjects. Canonical and machine learning approaches are applied for identification of the genomic biomarkers for glioblastoma and their crosstalks. In our study, 97 genes appeared enriched in 7 oncogenic pathways (RAF-MAPK, P53, PRC2-EZH2, YAP conserved, MEK-MAPK, ErbB2 and STK33 signalling pathways) using GSEA, out of which 17 have been identified participating actively in crosstalks. Using PCA, 42 genes are found enriched in 7 pathways (cytoplasmic ribosomal proteins, translation factors, electron transport chain, ribosome, Huntington's disease, primary immunodeficiency pathways, and interferon type I signalling pathway) harbouring tumour when altered, out of which 25 actively participate in crosstalks. All the 14 pathways foster well-known cancer hallmarks and the identified DEGs can serve as genomic biomarkers, not only for the diagnosis and prognosis of Glioblastoma but also in providing a molecular foothold for oncogenic decision making in order to fathom the disease dynamics. Moreover, SNP analysis for the identified DEGs is performed to investigate their roles in disease dynamics in an elaborated manner. These results suggest that TEPs are capable of providing disease insights just like tumour cells with an advantage of being extracted anytime during the course of disease in order to monitor it.

Proteomic Analysis of Prostate Cancer FFPE Samples Reveals Markers of Disease Progression and Aggressiveness

Simple Summary

Prostate cancer (PCa) is the second most frequently diagnosed type of cancer in men. The lack of tools for accurate risk assessment is causing over-treatment of men with indolent PCa but also delayed detection of metastatic disease and thus high mortality. The aim of our study was to identify proteins related to the progression and aggressiveness of PCa that could serve as potential biomarkers for better risk stratification. To this end, we performed proteomic analysis of Formalin Fixed Paraffin Embedded (FFPE) prostate tissue specimens (n = 86) and compared them based on grade groups and biochemical recurrence status. Based on the valuable data generated by these comparisons, we have selected seven proteins (NMP1, UQCRH, HSPA9, MRPL3, VCAN, SERBP1, HSPE1) as common denominators of PCa aggressiveness and persistence that could potentially be used for the development of risk assessment tools. Notably, our observations are largely validated by transcriptomics data and literature.

Abstract

Prostate cancer (PCa) is the second most common cancer in men. Diagnosis and risk assessment are widely based on serum Prostate Specific Antigen (PSA) and biopsy, which might not represent the exact degree of PCa risk. Towards the discovery of biomarkers for better patient stratification, we performed proteomic analysis of Formalin Fixed Paraffin Embedded (FFPE) prostate tissue specimens using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). Comparative analysis of 86 PCa samples among grade groups 1–5 identified 301 significantly altered proteins. Additional analysis based on biochemical recurrence (BCR; BCR+ n = 14, BCR- n = 51) revealed 197 significantly altered proteins that indicate disease persistence. Filtering the overlapping proteins of these analyses, seven proteins (NPM1, UQCRH, HSPA9, MRPL3, VCAN, SERBP1, HSPE1) had increased expression in advanced grades and in BCR+/BCR- and may play a critical role in PCa aggressiveness. Notably, all seven proteins were significantly associated with progression in Prostate Cancer Transcriptome Atles (PCTA) and NPM1NPM1, UQCRH, and VCAN were further validated in The Cancer Genome Atlas (TCGA), where they were upregulated in BCR+/BCR-. UQCRH levels were also associated with poorer 5-year survival. Our study provides valuable insights into the key regulators of PCa progression and aggressiveness. The seven selected proteins could be used for the development of risk assessment tools.

Identification of a Three-Glycolysis-Related lncRNA Signature Correlated With Prognosis and Metastasis in Clear Cell Renal Cell Carcinoma

The clear cell renal cell carcinoma (ccRCC) is not only a malignant disease but also an energy metabolic disease, we aimed to identify a novel prognostic model based on glycolysis-related long non-coding RNA (lncRNAs) and explore its mechanisms. With the use of Pearson correlation analysis between the glycolysis-related differentially expressed genes and lncRNAs from The Cancer Genome Atlas (TCGA) dataset, we identified three glycolysis-related lncRNAs and successfully constructed a prognostic model based on their expression. The diagnostic efficacy and the clinically predictive capacity of the signature were evaluated by univariate and multivariate Cox analyses, Kaplan–Meier survival analysis, and principal component analysis (PCA). The glycolysis-related lncRNA signature was constructed based on the expressions of AC009084.1, AC156455.1, and LINC00342. Patients were grouped into high- or low-risk groups according to risk score demonstrated significant differences in overall survival (OS) period, which were validated by patients with ccRCC from the International Cancer Genome Consortium (ICGC) database. Univariate Cox analyses, multivariate Cox analyses, and constructed nomogram-confirmed risk score based on our signature were independent prognosis predictors. The CIBERSORT algorithms demonstrated significant correlations between three-glycolysis-related lncRNAs and the tumor microenvironment (TME) components. Functional enrichment analysis demonstrated potential pathways and processes correlated with the risk model. Clinical samples validated expression levels of three-glycolysis-related lncRNAs, and LINC00342 demonstrated the most significant aberrant expression. in vitro, the general overexpression of LINC00342 was detected in ccRCC cells. After silencing LINC00342, the aberrant glycolytic levels and migration abilities in 786-O cells were decreased significantly, which might be explained by suppressed Wnt/β-catenin signaling pathway and reversed Epithelial mesenchymal transformation (EMT) process. Collectively, our research identified a novel three-glycolysis-related lncRNA signature as a promising model for generating accurate prognoses for patients with ccRCC, and silencing lncRNA LINC00342 from the signature could partly inhibit the glycolysis level and migration of ccRCC cells.

Probabilistic ratiocination of hepatocellular carcinoma after resection: evaluation of expected to be promising approaches

Background

Precise prediction of survival after treatment is of great importance for patients with diseases with high mortality. RNA sequencing data and deep learning (DL) methods are expected to become promising approaches in the development of prediction models in the future. We aimed to evaluate the optimal covariates and methodology for patients with hepatocellular carcinoma (HCC) undergoing surgical resection.

Methods

The Cox proportional hazards regression model and the DL approach were used to develop prediction models incorporating clinical, genetic, and combined clinical and genetic variables for survival prediction in patients with HCC after resection. A total of 1,114 patients and 184 patients were enrolled in the present study from 2,163 and 601 patients from Eastern Hepatobiliary Surgery Hospital and Renji Hospital, respectively. The models were internally validated through random sampling and externally validated in clinical cohorts. Between-model comparisons were carried out in terms of the integrated discrimination improvement and net reclassification index.

Results

The Cox and DL clinical models were developed by adopting 7 independent prognostic factors (total bilirubin, prothrombin time, tumor size, tumor number, lymph node metastasis, and vascular invasion) and 22 clinical factors, respectively. Both the Cox clinical model and the DL clinical model showed excellent performances in the derivation [area under the curve (AUC): 0.75 vs. 0.77] and validation (AUC: 0.83 vs. 0.80) sets. The derived Cox genetic model with 6 significant prognostic genes was not as effective as the DL approach involving 686 genes. A combined clinical and genetic approach modified the performances of both the Cox and DL models. The integrated discrimination improvement and net reclassification index of the DL clinical model were generally better than those of the Cox clinical model.

Conclusions

Our Cox clinical model sufficiently provided precise survival prediction in patients with HCC after resection. It may serve as an accurate and cost-effective tool for predicting survival in such patients.

Identification of a Novel Glycolysis-Related Gene Signature for Predicting Breast Cancer Survival

To identify a glycolysis-related gene signature for the evaluation of prognosis in patients with breast cancer, we analyzed the data of a training set from TCGA database and four validation cohorts from the GEO and ICGC databases which included 1,632 patients with breast cancer. We conducted GSEA, univariate Cox regression, LASSO, and multiple Cox regression analysis. Finally, an 11-gene signature related to glycolysis for predicting survival in patients with breast cancer was developed. And Kaplan–Meier analysis and ROC analyses suggested that the signature showed a good prognostic ability for BC in the TCGA, ICGC, and GEO datasets. The analyses of univariate Cox regression and multivariate Cox regression revealed that it’s an important prognostic factor independent of multiple clinical features. Moreover, a prognostic nomogram, combining the gene signature and clinical characteristics of patients, was constructed. These findings provide insights into the identification of breast cancer patients with a poor prognosis.

Abnormal Expression of Mitochondrial Ribosomal Proteins and Their Encoding Genes with Cell Apoptosis and Diseases

Mammalian mitochondrial ribosomes translate 13 proteins encoded by mitochondrial genes, all of which play roles in the mitochondrial respiratory chain. After a long period of reconstruction, mitochondrial ribosomes are the most protein-rich ribosomes. Mitochondrial ribosomal proteins (MRPs) are encoded by nuclear genes, synthesized in the cytoplasm and then, transported to the mitochondria to be assembled into mitochondrial ribosomes. MRPs not only play a role in mitochondrial oxidative phosphorylation (OXPHOS). Moreover, they participate in the regulation of cell state as apoptosis inducing factors. Abnormal expressions of MRPs will lead to mitochondrial metabolism disorder, cell dysfunction, etc. Many researches have demonstrated the abnormal expression of MRPs in various tumors. This paper reviews the basic structure of mitochondrial ribosome, focuses on the structure and function of MRPs, and their relationships with cell apoptosis and diseases. It provides a reference for the study of the function of MRPs and the disease diagnosis and treatment.

Using ESTIMATE algorithm to establish an 8-mRNA signature prognosis prediction system and identify immunocyte infiltration-related genes in Pancreatic adenocarcinoma

OBJECTIVE

The tumour microenvironment is one of the significant factors driving the carcinogenesis of Pancreatic adenocarcinoma (PAAD). However, the underlying mechanism of how the tumour microenvironment impacts the prognosis of PAAD is not completely clear.

RESULTS

The transcriptome and clinical data of 182 PAAD program cases were downloaded from the TCGA database. Three hundred thirty-three differentially expressed genes (DEGs) between high and low stromal groups and 314 DEGs between high and low immune score groups were identified using ESTIMATE score. Based on the 203 genes differentially expressed simultaneously in two score-related comparisons, we established an 8-mRNA signature to evaluate the prognosis of PAAD patients. Kaplan-Meier curves showed significantly worse survival for patients with high-risk scores in both the training and validation groups. The risk score was an independent prognostic factor and had a high predictive value for the prognosis of patients with PAAD. By searching the TCGA database, we showed that CA9, CXCL9, and GIMAP7 from the 8-mRNA signature were associated with the infiltration levels of immunocytes by regulating FOXO1 expression in PAAD.

CONCLUSIONS

Unlike traditional methods of screening for differential genes in cancer and healthy tissues, we constructed a novel 8-mRNA signature to predict the prognosis of PAAD patients by applying ESTIMATE scoring to RNA-seq-based transcriptome data. Most importantly, we identified CA9, CXCL9, and GIMAP7 from the above eight genes as regulators of immunocyte infiltration by adjusting the expression of FOXO1 in PAAD. Thus, CA9, CXCL9, and GIMAP7 might be the ideal targets of immune therapy of PAAD.

METHODS

ESTIMATE scoring was used to determine the stromal and immune scores of transcriptome datasets downloaded from the TCGA database. An mRNA-based prognostic signature was built for the training cohort via the LASSO Cox regression model. The signature was verified using a validation cohort. Kaplan-Meier curves and log-rank analysis were used to identify survival differences. Western blot analysis and RT-qPCR analysis were carried out to analyze the expression of specific proteins and mRNAs. IHC was performed to assess the protein levels of Forkhead box-O 1 (FOXO1), Carbonic anhydrase 9 (CA9), C-X-C motif chemokine ligand 9 (CXCL9), and GTPase, IMAP family member 7 (GIMAP7) in the tissue microarray of PAAD.

Glycolysis-Based Genes Associated with the Clinical Outcome of Pancreatic Ductal Adenocarcinoma Identified by The Cancer Genome Atlas Data Analysis

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadly tumors in digestive tract tumors. Although there has been advancement in PDAC treatment, its prognosis still remains unsatisfactory, mainly because of dismal diagnosis. This article aims to develop new prognostic factors related to energy metabolism in PDAC and to use these genes for novel risk stratification. Hundred fifty messenger RNA (mRNA) expression profiles and clinicopathological data of PDAC were downloaded from The Cancer Genome Atlas dataset. The glycolysis pathway was the significant pathway based on the gene set enrichment analysis. We chose the glycolysis pathway-related 176 genes for further analysis. Multivariate Cox regression analysis and forward stepwise Cox regression model established a novel three-gene glycolytic signature (including MET, B3GNT3, and SPAG4) for PDAC patients' prognosis prediction. All 150 patients were classified into two groups by the median risk score. High-risk group had a worse outcome compared to the low-risk group. The risk score was also significantly correlated with age and radiotherapy. A nomogram, including the glycolytic gene signature, has shown some clinical net benefit for overall survival prediction. We also validated the validity and reliability in the Puleo dataset. This novel gene expression signature may be involved in the pathophysiology and used for risk stratification and prognosis prediction in PDAC.

Gene signature associated with neuro-endocrine activity predicting prognosis of pancreatic carcinoma

Background

Genomic analysis is the promising tool to clear understanding of the tumorigenesis and guide molecular classification for pancreatic cancer. Our purpose was to develop a critical predictive model for prognosis in pancreatic carcinoma, based on the genomic data.

Methods

The online The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets were queried as training and validation cohorts for comprehensive bioinformatic analysis. We applied Lasso and multivariate Cox regression to shrink genes and construct predictive model.

Results

A four genes model (DNAH10: HR = 0.71, 95% CI = 0.57–0.88, HSBP1L1: HR = 1.51, 95% CI = 1.18–1.92, KIAA0513: HR = 0.69, 95% CI = 0.50–0.96, and MRPL3: HR = 3.73, 95% CI = 2.03–6.86), was proposed and validated. The C‐index was 0.73 (95% CI: 0.7–0.77). Patients in high‐risk and low‐risk group, stratified by model, suffered significantly different overall survival time (15.1 vs. 49.3 months, p < 0.0001 in TCGA; 423 vs. 618 days, p = 0.038 in ICGC). Taken clinical parameters into consideration, the risk‐score was independent marker in clinical subpopulation. To explore the molecular mechanisms, 579 differential expression genes (DEG) in two groups were identified by edgeR. Functional enrichment of DEG indicated neuro‐endocrine activity was the potential mechanism for the discrepant prognosis.

Conclusion

A specific four genes signature with the ability to predicted survival of pancreatic carcinoma was generated, which may indicate the connection between neuro‐endocrine activity and patients’ prognosis.