ABSTRACT
BACKGROUND: Colorectal cancer (CRC) prognosis prediction is currently a major challenge. Epigenetic regulation has been widely reported for its role in cancer development. AIM: To construct a robust prognostic signature, we used developed and validated across datasets. METHODS: After constructing the signature, the prognostic value of the signature was evaluated in the TCGA cohort and six independent datasets (GSE17526, GSE17537, GSE33113, GSE37892, GSE39048 and GSE39582). The clinical, genomic and transcriptomic features related to the signature were identified. The correlations of the signature score with immune cell infiltration and cell-cell interactions were analyzed. The correlations between the signature score and the sensitivity to different drugs were also predicted. RESULTS: In the TCGA cohort, patients in the low-risk group according to the signature score had longer survival than those in the high-risk group, and this finding was validated in the validation datasets. The signature was a prognostic factor independent of age and sex and was correlated with stage and PD-1/PD-L1 expression. Area under the receiving operating characteristic curve was 0.72. Genomic association analyses revealed that samples from high-risk patients exhibited chromosomal instability. Transcriptomic analyses revealed that the signature score was significantly associated with multiple cellular pathways. Bulk RNA-seq and single-cell sequencing data revealed that the signature reflected differences in infiltrating immune cell-tumor cell interactions, especially for macrophages. The signature also predicted the putative drug sensitivity of CRC samples. CONCLUSION: The signature is a valuable biomarker for predicting CRC prognosis and reflects multiple features of CRC, especially macrophage infiltration in the microenvironment.
ABSTRACT
BACKGROUND: Early menarche is an established risk factor for breast cancer but its molecular contribution to tumor biology and prognosis remains unclear. METHODS: We profiled transcriptome-wide gene expression in breast tumors (N = 846) and tumor-adjacent normal tissues (N = 666) from women in the Nurses' Health Studies (NHS) to investigate whether early menarche (age < 12) is associated with tumor molecular and prognostic features in women with breast cancer. Multivariable linear regression and pathway analyses using competitive gene set enrichment analysis were conducted in both tumor and adjacent-normal tissue and externally validated in TCGA (N = 116). Subgroup analyses stratified on ER-status based on the tumor were also performed. PAM50 signatures were used for tumor molecular subtyping and to generate proliferation and risk of recurrence scores. We created a gene expression score using LASSO regression to capture early menarche based on 28 genes from FDR-significant pathways in breast tumor tissue in NHS and tested its association with 10-year disease-free survival in both NHS (N = 836) and METABRIC (N = 952). RESULTS: Early menarche was significantly associated with 369 individual genes in adjacent-normal tissues implicated in extracellular matrix, cell adhesion, and invasion (FDR ≤ 0.1). Early menarche was associated with upregulation of cancer hallmark pathways (18 significant pathways in tumor, 23 in tumor-adjacent normal, FDR ≤ 0.1) related to proliferation (e.g. Myc, PI3K/AKT/mTOR, cell cycle), oxidative stress (e.g. oxidative phosphorylation, unfolded protein response), and inflammation (e.g. pro-inflammatory cytokines IFN α and IFN γ ). Replication in TCGA confirmed these trends. Early menarche was associated with significantly higher PAM50 proliferation scores (ß = 0.082 [0.02-0.14]), odds of aggressive molecular tumor subtypes (basal-like, OR = 1.84 [1.18-2.85] and HER2-enriched, OR = 2.32 [1.46-3.69]), and PAM50 risk of recurrence score (ß = 4.81 [1.71-7.92]). Our NHS-derived early menarche gene expression signature was significantly associated with worse 10-year disease-free survival in METABRIC (N = 952, HR = 1.58 [1.10-2.25]). CONCLUSIONS: Early menarche is associated with more aggressive molecular tumor characteristics and its gene expression signature within tumors is associated with worse 10-year disease-free survival among women with breast cancer. As the age of onset of menarche continues to decline, understanding its relationship to breast tumor characteristics and prognosis may lead to novel secondary prevention strategies.
Subject(s)
Breast Neoplasms , Gene Expression Profiling , Menarche , Neoplasm Recurrence, Local , Transcriptome , Humans , Female , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Breast Neoplasms/mortality , Menarche/genetics , Neoplasm Recurrence, Local/genetics , Neoplasm Recurrence, Local/pathology , Middle Aged , Prognosis , Adult , Biomarkers, Tumor/genetics , Risk Factors , Gene Expression Regulation, Neoplastic , Age FactorsABSTRACT
Modern advanced radiotherapy techniques have improved the precision and accuracy of radiotherapy delivery, with resulting plans being highly personalised based on individual anatomy. Adaptation for individual tumour biology remains elusive. There is an unmet need for biomarkers of intrinsic radiosensitivity that can predict tumour response to radiation to facilitate individualised decision-making, dosing and treatment planning. Over the last few decades, the use of high throughput molecular biology technologies has led to an explosion of newly discovered cancer biomarkers. Gene expression signatures are now used routinely in clinic to aid decision-making regarding adjuvant systemic therapy. They have great potential as radiotherapy biomarkers. A previous systematic review published in 2015 reported only five studies of signatures evaluated for their ability to predict radiotherapy benefits in clinical cohorts. This updated systematic review encompasses the expanded number of studies reported in the last decade. An additional 27 studies were identified. In total, 22 distinct signatures were recognised (5 pre-2015, 17 post-2015). Seventeen signatures were 'radiosensitivity' signatures and five were breast cancer prognostic signatures aiming to identify patients at an increased risk of local recurrence and therefore were more likely to benefit from adjuvant radiation. Most signatures (15/22) had not progressed beyond the discovery phase of development, with no suitable validated clinical-grade assay for application. Very few signatures (4/17 'radiosensitivity' signatures) had undergone any laboratory-based biological validation of their ability to predict tumour radiosensitivity. No signatures have been assessed prospectively in a phase III biomarker-led trial to date and none are recommended for routine use in clinical guidelines. A phase III prospective evaluation is ongoing for two breast cancer prognostic signatures. The most promising radiosensitivity signature remains the radiosensitivity index (RSI), which is used to calculate a genomic adjusted radiation dose (GARD). There is an ongoing phase II prospective biomarker-led study of RSI/GARD in triple negative breast cancer. The results of these trials are eagerly anticipated over the coming years. Future work in this area should focus on (1) robust biological validation; (2) building biobanks alongside large radiotherapy randomised controlled trials with dose variance (to demonstrate an interaction between radiosensitivity signature and dose); (3) a validation of clinical-grade cost-effective assays that are deliverable within current healthcare infrastructure; and (4) an integration with biomarkers of other determinants of radiation response.
ABSTRACT
The prognostication of survival trajectories in multiple myeloma (MM) patients presents a substantial clinical challenge. Leveraging transcriptomic and clinical profiles from an expansive cohort of 2,088 MM patients, sourced from the Gene Expression Omnibus and The Cancer Genome Atlas repositories, we applied a sophisticated nested lasso regression technique to construct a prognostic model predicated on 28 gene pairings intrinsic to cell death pathways, thereby deriving a quantifiable risk stratification metric. Employing a threshold of 0.15, we dichotomized the MM samples into discrete high-risk and low-risk categories. Notably, the delineated high-risk cohort exhibited a statistically significant diminution in survival duration, a finding which consistently replicated across both training and external validation datasets. The prognostic acumen of our cell death signature was further corroborated by TIME ROC analyses, with the model demonstrating robust performance, evidenced by AUC metrics consistently surpassing the 0.6 benchmark across the evaluated arrays. Further analytical rigor was applied through multivariate COX regression analyses, which ratified the cell death risk model as an independent prognostic determinant. In an innovative stratagem, we amalgamated this risk stratification with the established International Staging System (ISS), culminating in the genesis of a novel, refined ISS categorization. This tripartite classification system was subjected to comparative analysis against extant prognostic models, whereupon it manifested superior predictive precision, as reflected by an elevated C-index. In summation, our endeavors have yielded a clinically viable gene pairing model predicated on cellular mortality, which, when synthesized with the ISS, engenders an augmented prognostic tool that exhibits pronounced predictive prowess in the context of multiple myeloma.
Subject(s)
Cell Death , Multiple Myeloma , Multiple Myeloma/pathology , Multiple Myeloma/genetics , Multiple Myeloma/mortality , Humans , Prognosis , Male , Female , Risk Assessment , Gene Expression Profiling , Middle Aged , Neoplasm Staging , Aged , Survival AnalysisABSTRACT
PURPOSE: Brain metastasis (BM) in colorectal cancer (CRC) is a rare event with poor prognosis. Apart from (K)RAS status and lung and bone metastasis no biomarkers exist to identify patients at risk. This study aimed to identify a gene expression signature associated with colorectal BM. METHODS: Three patient groups were formed: 1. CRC with brain metastasis (BRA), 2. exclusive liver metastasis (HEP) and, 3. non-metastatic disease (M0). RNA was extracted from primary tumors and mRNA expression was measured using a NanoString Panel (770 genes). Expression was confirmed by qPCR in a validation cohort. Statistical analyses including multivariate logistic regression followed by receiver operating characteristic (ROC) analysis were performed. RESULTS: EMILIN3, MTA1, SV2B, TMPRSS6, ACVR1C, NFAT5 and SMC3 were differentially expressed in BRA and HEP/M0 groups. In the validation cohort, differential NFAT5, ACVR1C and SMC3 expressions were confirmed. BRA patients showed highest NFAT5 levels compared to HEP/M0 groups (global p = 0.02). High ACVR1C expression was observed more frequently in the BRA group (42.9%) than in HEP (0%) and M0 (7.1%) groups (global p = 0.01). High SMC3 expressions were only detectable in the BRA group (global p = 0.003). Only patients with BM showed a combined high expression of NFAT5, ACVR1C or SMC3 as well as of all three genes. ROC analysis revealed a good prediction of brain metastasis by the three genes (area under the curve (AUC) = 0.78). CONCLUSIONS: The NFAT5, ACVR1C and SMC3 gene expression signature is associated with colorectal BM. Future studies should further investigate the importance of this biomarker signature.
Subject(s)
Biomarkers, Tumor , Brain Neoplasms , Colorectal Neoplasms , Humans , Brain Neoplasms/secondary , Brain Neoplasms/genetics , Colorectal Neoplasms/pathology , Colorectal Neoplasms/genetics , Male , Female , Middle Aged , Biomarkers, Tumor/genetics , Biomarkers, Tumor/metabolism , Aged , Transcriptome , Liver Neoplasms/secondary , Liver Neoplasms/genetics , Liver Neoplasms/metabolism , Prognosis , Gene Expression Profiling , ROC Curve , Adult , Gene Expression Regulation, NeoplasticABSTRACT
BACKGROUND/AIM: Patients with triple-negative breast cancer (TNBC) have a high rate of recurrence within 3 years of diagnosis and a high rate of death within 5 years compared to other subtypes. The number of clinical trials investigating various new agents and combination therapies has recently increased; however, current strategies benefit only a minority of patients. This study aimed to identify specific genes that predict patients at high risk of recurrence and the immune status of the tumor microenvironment at an early stage, thereby providing insight into potential therapeutic targets to improve clinical outcomes in TNBC patients. MATERIALS AND METHODS: We evaluated the prognostic significance of microarray mRNA expression of 20,603 genes in 233 TNBC patients from the METABRIC dataset and further validated the results using RNA-seq mRNA expression data in 143 TNBC patients from the GSE96058 dataset. RESULTS: Eighteen differentially expressed genes (AKNA, ARHGAP30, CA9, CD3D, CD3G, CD6, CXCR6, CYSLTR1, DOCK10, ENO1, FLT3LG, IFNG, IL2RB, LPXN, PRKCB, PVRIG, RASSF5, and STAT4) identified in both datasets were found to be reliable biomarkers for predicting TNBC recurrence and progression. Notably, the genes whose low expression was associated with increased risk of recurrence and death were immune-related genes, with significant differences in levels of immune cell infiltration in the tumor microenvironment between high- and low- expression groups. CONCLUSION: Genes reported herein may be effective biomarkers to identify TNBC patients who will and will not benefit from immunotherapy and may be particularly important genes for developing future treatment strategies, including immunotherapy.
Subject(s)
Neoplasm Recurrence, Local , Triple Negative Breast Neoplasms , Humans , Triple Negative Breast Neoplasms/genetics , Triple Negative Breast Neoplasms/immunology , Triple Negative Breast Neoplasms/pathology , Female , Neoplasm Recurrence, Local/genetics , Neoplasm Recurrence, Local/pathology , Neoplasm Recurrence, Local/immunology , Prognosis , Biomarkers, Tumor/genetics , Tumor Microenvironment/immunology , Tumor Microenvironment/genetics , Gene Expression Profiling , Gene Expression Regulation, Neoplastic , Transcriptome , Lymphocytes, Tumor-Infiltrating/immunology , Lymphocytes, Tumor-Infiltrating/metabolism , Middle AgedABSTRACT
OBJECTIVE: To propose a new gene expression signature that identifies endometrial disruptions independent of endometrial luteal phase timing and predicts if patients are at risk of endometrial failure. DESIGN: Multicentric, prospective study. SETTING: Reproductive medicine research department in a public hospital affiliated with private fertility clinics and a reproductive genetics laboratory. PATIENTS: Caucasian women (n = 281; 39.4 ± 4.8 years old with a body mass index of 22.9 ± 3.5 kg/m2) undergoing hormone replacement therapy between July 2018 and July 2021. Endometrial samples from 217 patients met RNA quality criteria for signature discovery and analysis. INTERVENTION(S): Endometrial biopsies collected in the mid-secretory phase. MAIN OUTCOME MEASURE(S): Endometrial luteal phase timing-corrected expression of 404 genes and reproductive outcomes of the first single embryo transfer (SET) after biopsy collection to identify prognostic biomarkers of endometrial failure. RESULTS: Removal of endometrial timing variation from gene expression data allowed patients to be stratified into poor (n = 137) or good (n = 49) endometrial prognosis groups on the basis of their clinical and transcriptomic profiles. Significant differences were found between endometrial prognosis groups in terms of reproductive rates: pregnancy (44.6% vs. 79.6%), live birth (25.6% vs. 77.6%), clinical miscarriage (22.2% vs. 2.6%), and biochemical miscarriage (20.4% vs. 0%). The relative risk of endometrial failure for patients predicted as a poor endometrial prognosis was 3.3 times higher than those with a good prognosis. The differences in gene expression between both profiles were proposed as a biomarker, coined the endometrial failure risk (EFR) signature. Poor prognosis profiles were characterized by 59 upregulated and 63 downregulated genes mainly involved in regulation (17.0%), metabolism (8.4%), immune response, and inflammation (7.8%). This EFR signature had a median accuracy of 0.92 (min = 0.88, max = 0.94), median sensitivity of 0.96 (min = 0.91, max = 0.98), and median specificity of 0.84 (min = 0.77, max = 0.88), positioning itself as a promising biomarker for endometrial evaluation. CONCLUSION(S): The EFR signature revealed a novel endometrial disruption, independent of endometrial luteal phase timing, present in 73.7% of patients. This EFR signature stratified patients into 2 significantly distinct and clinically relevant prognosis profiles providing opportunities for personalized therapy. Nevertheless, further validations are needed before implementing this gene signature as an artificial intelligence (AI)-based tool to reduce the risk of patients experiencing endometrial failure.
Subject(s)
Endometrium , Humans , Female , Endometrium/metabolism , Endometrium/drug effects , Endometrium/pathology , Adult , Prospective Studies , Pregnancy , Biomarkers/metabolism , Infertility, Female/genetics , Infertility, Female/therapy , Infertility, Female/diagnosis , Infertility, Female/metabolism , Pregnancy Rate , Hormone Replacement Therapy , Risk Assessment , Time Factors , Transcriptome , Embryo Transfer , Risk Factors , Luteal Phase/metabolism , Gene Expression Profiling , Predictive Value of TestsABSTRACT
PURPOSE: Medulloblastoma is the most common childhood malignant brain tumor and is a leading cause of cancer-related death in children. Recent transcriptional studies have shown that medulloblastomas comprise at least four molecular subgroups, each with distinct demographics, genetics, and clinical outcomes. Medulloblastoma subtyping has become critical for subgroup-specific therapies. The use of gene expression assays to determine the molecular subgroup of clinical specimens is a long-awaited application of molecular biology for this pediatric cancer. METHODS: In the current study, we established a medulloblastoma transcriptome database of 460 samples retrieved from three published datasets (GSE21140, GSE37382, and GSE37418). With this database, we identified a 23-gene signature that is significantly associated with the medulloblastoma subgroups and achieved a classification accuracy of 95.2%. RESULTS: The 23-gene signature was further validated in a long-term cohort of 142 Chinese medulloblastoma patients. The 23-gene signature classified 21 patients as WNT (15%), 41 as SHH (29%), 16 as Group 3 (11%), and 64 as Group 4 (45%). For patients of WNT, SHH, Group 3, and Group 4, 5-year overall-survival rate reached 80%, 62%, 27%, and 47%, respectively (p < 0.0001), meanwhile 5-year progression-free survival reached 80%, 52%, 27%, and 45%, respectively (p < 0.0001). Besides, SHH/TP53-mutant tumors were associated with worse prognosis compared with SHH/TP53 wild-type tumors and other subgroups. We demonstrated that subgroup assignments by the 23-gene signature and Northcott's NanoString assay were highly comparable with a concordance rate of 96.4%. CONCLUSIONS: In conclusion, we present a novel gene signature that is capable of accurately and reliably assigning FFPE medulloblastoma samples to their molecular subgroup, which may serve as an auxiliary tool for medulloblastoma subtyping in the clinic. Future incorporation of this gene signature into prospective clinical trials is warranted to further evaluate its clinical.
Subject(s)
Brain Neoplasms , Cerebellar Neoplasms , Medulloblastoma , Humans , Child , Medulloblastoma/diagnosis , Medulloblastoma/genetics , Transcriptome/genetics , Prospective Studies , Cerebellar Neoplasms/genetics , ChinaABSTRACT
Models based on risk stratification are increasingly reported for Diffuse large B cell lymphoma (DLBCL). Due to a rising interest in nomograms for cancer patients, we aimed to review and critically appraise prognostic models based on nomograms in DLBCL patients. A literature search in PubMed/Embase identified 59 articles that proposed prognostic models for DLBCL by combining parameters of interest (e.g., clinical, laboratory, immunohistochemical, and genetic) between January 2000 and 2024. Of them, 40 studies proposed different gene expression signatures and incorporated them into nomogram-based prognostic models. Although most studies assessed discrimination and calibration when developing the model, many lacked external validation. Current nomogram-based models for DLBCL are mainly developed from publicly available databases, lack external validation, and have no applicability in clinical practice. However, they may be helpful in individual patient counseling, although careful considerations should be made regarding model development due to possible limitations when choosing nomograms for prognostication.
Subject(s)
Lymphoma, Large B-Cell, Diffuse , Nomograms , Lymphoma, Large B-Cell, Diffuse/diagnosis , Lymphoma, Large B-Cell, Diffuse/genetics , Lymphoma, Large B-Cell, Diffuse/pathology , Humans , Prognosis , Risk Assessment/methodsABSTRACT
BACKGROUND: Although substantial efforts have been made to build molecular biomarkers to predict radiation sensitivity, the ability to accurately stratify the patients is still limited. In this study, we aim to leverage large-scale radiogenomics datasets to build genomic predictors of radiation response using the integral of the radiation dose-response curve. METHODS: Two radiogenomics datasets consisting of 511 and 60 cancer cell lines were utilized to develop genomic predictors of radiation sensitivity. The intrinsic radiation sensitivity, defined as the integral of the dose-response curve (AUC) was used as the radioresponse variable. The biological determinants driving AUC and SF2 were compared using pathway analysis. To build the predictive model, the largest and smallest datasets consisting of 511 and 60 cancer cell lines were used as the discovery and validation cohorts, respectively, with AUC as the response variable. RESULTS: Utilizing a compendium of three pathway databases, we illustrated that integral of the radiobiological model provides a more comprehensive characterization of molecular processes underpinning radioresponse compared to SF2. Furthermore, more pathways were found to be unique to AUC than SF2-30, 288 and 38 in KEGG, REACTOME and WIKIPATHWAYS, respectively. Also, the leading-edge genes driving the biological pathways using AUC were unique and different compared to SF2. With regards to radiation sensitivity gene signature, we obtained a concordance index of 0.65 and 0.61 on the discovery and validation cohorts, respectively. CONCLUSION: We developed an integrated framework that quantifies the impact of physical radiation dose and the biological effect of radiation therapy in interventional pre-clinical model systems. With the availability of more data in the future, the clinical potential of this signature can be assessed, which will eventually provide a framework to integrate genomics into biologically-driven precision radiation oncology.
Subject(s)
Neoplasms , Transcriptome , Humans , Radiation Tolerance/genetics , Neoplasms/genetics , Neoplasms/radiotherapy , Cell Line , BiomarkersABSTRACT
Isoniazid is a first-line drug in antitubercular therapy. Isoniazid is one of the most commonly used drugs that can cause liver injury or acute liver failure, leading to death or emergency liver transplantation. Therapeutic approaches for the prevention of isoniazid-induced liver injury are yet to be established. In this study, we identified the gene expression signature for isoniazid-induced liver injury using a public transcriptome dataset, focusing on the differences in susceptibility to isoniazid in various mouse strains. We predicted that lansoprazole is a potentially protective drug against isoniazid-induced liver injury using connectivity mapping and an adverse event reporting system. We confirmed the protective effects of lansoprazole against isoniazid-induced liver injury using zebrafish and patients' electronic health records. These results suggest that lansoprazole can ameliorate isoniazid-induced liver injury. The integrative approach used in this study may be applied to identify novel functions of clinical drugs, leading to drug repositioning.
ABSTRACT
The RATIONALE-307 study (ClinicalTrials.gov: NCT03594747) demonstrates prolonged progression-free survival (PFS) with first-line tislelizumab plus chemotherapy versus chemotherapy in advanced lung squamous cell carcinoma (LUSC; N = 360). Here we describe an immune-related gene expression signature (GES), composed of genes involved in both innate and adaptive immunity, that appears to differentiate tislelizumab plus chemotherapy PFS benefit versus chemotherapy. In contrast, a tislelizumab plus chemotherapy PFS benefit is observed regardless of programmed death ligand 1 (PD-L1) expression or tumor mutational burden (TMB). Genetic analysis reveals that NRF2 pathway activation is enriched in PD-L1positive and TMBhigh patients. NRF2 pathway activation is negatively associated with PFS, which affects efficacy outcomes associated with PD-L1 and TMB status, impairing their predictive potential. Mechanistic studies demonstrate that NRF2 directly mediates PD-L1 constitutive expression independent of adaptive PD-L1 regulation in LUSC. In summary, the GES is an immune signature that might identify LUSC patients likely to benefit from first-line tislelizumab plus chemotherapy.
Subject(s)
Carcinoma, Non-Small-Cell Lung , Carcinoma, Squamous Cell , Lung Neoplasms , Humans , B7-H1 Antigen/genetics , Biomarkers, Tumor/genetics , Carcinoma, Non-Small-Cell Lung/drug therapy , Carcinoma, Squamous Cell/drug therapy , Carcinoma, Squamous Cell/genetics , Carcinoma, Squamous Cell/pathology , Lung/pathology , Lung Neoplasms/drug therapy , Lung Neoplasms/genetics , Lung Neoplasms/pathology , Mutation , NF-E2-Related Factor 2/genetics , Programmed Cell Death 1 Receptor , Treatment Outcome , Tumor Microenvironment/geneticsABSTRACT
Epithelial-mesenchymal transition (EMT) is a cellular plasticity program critical for embryonic development and tissue regeneration, and aberrant EMT is associated with disease including cancer. The high degree of plasticity in the mammary epithelium is reflected in extensive heterogeneity among breast cancers. Here, we have analyzed RNA-sequencing data from three different mammary epithelial cell line-derived EMT models and identified a robust mammary EMT gene expression signature that separates breast cancers into distinct subgroups. Most strikingly, the basal-like breast cancers form two subgroups displaying partial-EMT and post-EMT gene expression patterns. We present evidence that key EMT-associated transcription factors play distinct roles at different stages of EMT in mammary epithelial cells.
ABSTRACT
Background: Diffuse large B-cell lymphoma (DLBCL) is a hematological malignancy representing one-third of non-Hodgkin's lymphoma cases. Notwithstanding immunotherapy in combination with chemotherapy (R-CHOP) is an effective therapeutic approach for DLBCL, a subset of patients encounters treatment resistance, leading to low survival rates. Thus, there is an urgent need to identify predictive biomarkers for DLBCL including the elderly population, which represents the fastest-growing segment of the population in Western countries. Methods: Gene expression profiles of n=414 DLBCL biopsies were retrieved from the public dataset GSE10846. Differentially expressed genes (DEGs) (fold change >1.4, p-value <0.05, n=387) have been clustered in responder and non-responder patient cohorts. An enrichment analysis has been performed on the top 30 up-regulated genes of responder and non-responder patients to identify the signatures involved in gene ontology (MSigDB). The more significantly up-regulated DEGs have been validated in our independent collection of formalin-fixed paraffin-embedded (FFPE) biopsy samples of elderly DLBCL patients, treated with R-CHOP as first-line therapy. Results: From the analysis of two independent cohorts of DLBCL patients emerged a gene signature able to predict the response to R-CHOP therapy. In detail, expression levels of EBF1, MYO6, CALR are associated with a significant worse overall survival. Conclusions: These results pave the way for a novel characterization of DLBCL biomarkers, aiding the stratification of responder versus non-responder patients.
Subject(s)
Lymphoma, Large B-Cell, Diffuse , Lymphoma, Non-Hodgkin , Humans , Aged , Antibodies, Monoclonal, Murine-Derived/therapeutic use , Lymphoma, Large B-Cell, Diffuse/drug therapy , Lymphoma, Large B-Cell, Diffuse/genetics , Lymphoma, Large B-Cell, Diffuse/pathology , Rituximab/therapeutic use , Lymphoma, Non-Hodgkin/drug therapy , Cyclophosphamide/therapeutic use , Vincristine/therapeutic use , Prednisone/therapeutic use , Doxorubicin/therapeutic use , Biomarkers , Trans-ActivatorsABSTRACT
Background: In recent years, with the development of transcriptome sequencing, the molecular characteristics of tumors are gradually revealed. Because of the complexity of tumor transcriptome, there is a need to look for the molecular signatures which can be used to evaluate the tissue origin and cell stemness of tumors in order to promote the diagnosis and treatment of tumors. Methods: Tumor tissue-specific gene sets (TTSGs) consisting of 200 genes were selected using RNA expression data of 9,875 patients from 33 tumor types. t-distributed Stochastic Neighbor Embedding (t-SNE) was used for dimensionality reduction and visualization of TTSGs in each tumor type. To evaluate oncogenic dedifferentiation and loss of cell stemness, Euclidean distance from each sample to a human embryo single-cell RNA-seq dataset (GSE36552) of TTSGs was calculated as TTSGs index indicating dissimilarity of tumors and embryo. TTSGs index was evaluated for prognosis in each tumor type. Two published signature indexes, the mRNA signature index (mRNAsi) and CIBERSORT, were compared to assess the correlation between the TTSGs index with cell stemness and immune microenvironment. Finally, the difference of prognosis, immune microenvironment and radiotherapy outcomes were compared between patients with high and low TTSGs index. Results: In this study, all 33 tumor types in The Cancer Genome Atlas (TCGA) were embedded into isolated clusters by t-SNE and confirmed by k-nearest neighbors (kNN) algorithm. Clusters of squamous-cell carcinoma were adjacent to each other revealing similar histologic origin. Basal-like breast cancer was separated from luminal and HER-2-amplified subtypes and closed to squamous-cell carcinoma. TTSGs index was related to overall survival outcomes in cancers derived from liver, thyroid, brain, cervical and kidney. There was a positive correlation between mRNAsi and TTSGs index in pan-kidney and pan-neuronal cancers. Furthermore, cell fractions of M2 macrophages and total leukocytes increased in the group with higher TTSGs index. Patients with higher TTSGs index had longer overall survival time and less radiation therapy resistance compared to patients with lower TTSGs index. Conclusions: The signature of TTSGs is related to tumor expression features that distinguish tumors of different histologic origin using t-SNE. The signature also relates to prognosis of certain kinds of tumors.
ABSTRACT
Background: Periodontits (PD) and Alzheimer's disease (AD) are both associated with ageing and clinical studies increasingly evidence their association. However, specific mechanisms underlying this association remain undeciphered, and immune-related processes are purported to play a signifcant role. The accrual of publicly available transcriptomic datasets permits secondary analysis and the application of data-mining and bioinformatic tools for biological discovery. Aim: The present study aimed to leverage publicly available transcriptomic datasets and databases, and apply a series of bioinformatic analysis to identify a robust signature of immune-related signature of PD and AD linkage. Methods: We downloaded gene-expresssion data pertaining PD and AD and identified crosstalk genes. We constructed a protein-protein network analysis, applied immune cell enrichment analysis, and predicted crosstalk immune-related genes and infiltrating immune cells. Next, we applied consisent cluster analysis and performed immune cell bias analysis, followed by LASSO regression to select biomarker immune-related genes. Results: The results showed a 3 gene set comprising of DUSP14, F13A1 and SELE as a robust immune-related signature. Macrophages M2 and NKT, B-cells, CD4+ memory T-cells and CD8+ naive T-cells emerged as key immune cells linking PD with AD. Conclusion: Candidate immune-related biomarker genes and immune cells central to the assocation of PD with AD were identified, and merit investigation in experimental and clinical research.
ABSTRACT
Background and Objectives: This study aimed to investigate the causes of continuous deep fluctuations in the absolute lymphocyte count (ALC) in an untreated patient with Chronic Lymphocytic Leukemia (CLL), who has had a favorable prognosis since the time of diagnosis. Up until now, the patient has voluntarily chosen to adopt a predominantly vegetarian and fruitarian diet, along with prolonged periods of total fasting (ranging from 4 to 39 days) each year. Materials and Methods: For this purpose, we decided to analyze the whole transcriptome profiling of peripheral blood (PB) CD19+ cells from the patient (#1) at different time-points vs. the same cells of five other untreated CLL patients who followed a varied diet. Consequently, the CLL patients were categorized as follows: the 1st group comprised patient #1 at 20 different time-points (16 time-points during nutrition and 4 time-points during fasting), whereas the 2nd group included only one time point for each of the patients (#2, #3, #4, #5, and #6) as they followed a varied diet. We performed microarray experiments using a powerful tool, the Affymetrix Human Clariom™ D Pico Assay, to generate high-fidelity biomarker signatures. Statistical analysis was employed to identify differentially expressed genes and to perform sample clustering. Results: The lymphocytosis trend in patient #1 showed recurring fluctuations since the time of diagnosis. Interestingly, we observed that approximately 4-6 weeks after the conclusion of fasting periods, the absolute lymphocyte count was reduced by about half. The gene expression profiling analysis revealed that nine genes were statistically differently expressed between the 1st group and the 2nd group. Specifically, IGLC3, RPS26, CHPT1, and PCDH9 were under expressed in the 1st group compared to the 2nd group of CLL patients. Conversely, IGHV3-43, IGKV3D-20, PLEKHA1, CYBB, and GABRB2 were over-expressed in the 1st group when compared to the 2nd group of CLL patients. Furthermore, clustering analysis validated that all the samples from patient #1 clustered together, showing clear separation from the samples of the other CLL patients. Conclusions: This study unveiled a small gene expression signature consisting of nine genes that distinguished an untreated CLL patient who followed prolonged periods of total fasting, maintaining a gradual growth trend of lymphocytosis, compared to five untreated CLL patients with a varied diet. Future investigations focusing on patient #1 could potentially shed light on the role of prolonged periodic fasting and the implication of this specific gene signature in sustaining the lymphocytosis trend and the favorable course of the disease.
Subject(s)
Fasting , Leukemia, Lymphocytic, Chronic, B-Cell , Transcriptome , Humans , Male , Female , Middle Aged , Aged , Aged, 80 and over , Cluster Analysis , Diet, Vegetarian , Leukemia, Lymphocytic, Chronic, B-Cell/genetics , LymphocytosisABSTRACT
Background: In the vaccine era, individuals receive multiple vaccines in their lifetime. Host gene expression in response to antigenic stimulation is usually virus-specific; however, identifying shared pathways of host response across a wide spectrum of vaccine pathogens can shed light on the molecular mechanisms/components which can be targeted for the development of broad/universal therapeutics and vaccines. Method: We isolated PBMCs, monocytes, B cells, and CD8+ T cells from the peripheral blood of healthy donors, who received both seasonal influenza vaccine (within <1 year) and smallpox vaccine (within 1 - 4 years). Each of the purified cell populations was stimulated with either influenza virus or vaccinia virus. Differentially expressed genes (DEGs) relative to unstimulated controls were identified for each in vitro viral infection, as well as for both viral infections (shared DEGs). Pathway enrichment analysis was performed to associate identified DEGs with KEGG/biological pathways. Results: We identified 2,906, 3,888, 681, and 446 DEGs in PBMCs, monocytes, B cells, and CD8+ T cells, respectively, in response to influenza stimulation. Meanwhile, 97, 120, 20, and 10 DEGs were identified as gene signatures in PBMCs, monocytes, B cells, and CD8+ T cells, respectively, upon vaccinia stimulation. The majority of DEGs identified in PBMCs were also found in monocytes after either viral stimulation. Of the virus-specific DEGs, 55, 63, and 9 DEGs occurred in common in PBMCs, monocytes, and B cells, respectively, while no DEGs were shared in infected CD8+ T cells after influenza and vaccinia. Gene set enrichment analysis demonstrated that these shared DEGs were over-represented in innate signaling pathways, including cytokine-cytokine receptor interaction, viral protein interaction with cytokine and cytokine receptor, Toll-like receptor signaling, RIG-I-like receptor signaling pathways, cytosolic DNA-sensing pathways, and natural killer cell mediated cytotoxicity. Conclusion: Our results provide insights into virus-host interactions in different immune cells, as well as host defense mechanisms against viral stimulation. Our data also highlights the role of monocytes as a major cell population driving gene expression in ex vivo PBMCs in response to viral stimulation. The immune response signaling pathways identified in this study may provide specific targets for the development of novel virus-specific therapeutics and improved vaccines for vaccinia and influenza. Although influenza and vaccinia viruses have been selected in this study as pathogen models, this approach could be applicable to other pathogens.
Subject(s)
Influenza Vaccines , Influenza, Human , Vaccinia , Humans , Vaccinia virus/genetics , Influenza, Human/genetics , CD8-Positive T-Lymphocytes , Transcriptome , VaccinationABSTRACT
Asthma is a common chronic disease that is characterized by respiratory symptoms including cough, wheeze, shortness of breath, and chest tightness. The underlying mechanisms of this disease are not fully elucidated, so more research is needed to identify better therapeutic compounds and biomarkers to improve disease outcomes. In this present study, we used bioinformatics to analyze the gene expression of adult asthma in publicly available microarray datasets to identify putative therapeutic molecules for this disease. We first compared gene expression in healthy volunteers and adult asthma patients to obtain differentially expressed genes (DEGs) for further analysis. A final gene expression signature of 49 genes, including 34 upregulated and 15 downregulated genes, was obtained. Protein-protein interaction and hub analyses showed that 10 genes, including POSTN, CPA3, CCL26, SERPINB2, CLCA1, TPSAB1, TPSB2, MUC5B, BPIFA1, and CST1, may be hub genes. Then, the L1000CDS2 search engine was used for drug repurposing studies. The top approved drug candidate predicted to reverse the asthma gene signature was lovastatin. Clustergram results showed that lovastatin may perturb MUC5B expression. Moreover, molecular docking, molecular dynamics simulation, and computational alanine scanning results supported the notion that lovastatin may interact with MUC5B via key residues such as Thr80, Thr91, Leu93, and Gln105. In summary, by analyzing gene expression signatures, hub genes, and therapeutic perturbation, we show that lovastatin is an approved drug candidate that may have potential for treating adult asthma.
Subject(s)
Asthma , Gene Expression Profiling , Humans , Adult , Gene Expression Profiling/methods , Molecular Docking Simulation , Asthma/drug therapy , Asthma/genetics , Genes, Regulator , Computational Biology/methods , Lovastatin , Gene Regulatory Networks , Protein Interaction Maps/geneticsABSTRACT
BACKGROUND: Type 1 diabetes is a complex heterogenous autoimmune disease without therapeutic interventions available to prevent or reverse the disease. This study aimed to identify transcriptional changes associated with the disease progression in patients with recent-onset type 1 diabetes. METHODS: Whole-blood samples were collected as part of the INNODIA study at baseline and 12 months after diagnosis of type 1 diabetes. We used linear mixed-effects modelling on RNA-seq data to identify genes associated with age, sex, or disease progression. Cell-type proportions were estimated from the RNA-seq data using computational deconvolution. Associations to clinical variables were estimated using Pearson's or point-biserial correlation for continuous and dichotomous variables, respectively, using only complete pairs of observations. FINDINGS: We found that genes and pathways related to innate immunity were downregulated during the first year after diagnosis. Significant associations of the gene expression changes were found with ZnT8A autoantibody positivity. Rate of change in the expression of 16 genes between baseline and 12 months was found to predict the decline in C-peptide at 24 months. Interestingly and consistent with earlier reports, increased B cell levels and decreased neutrophil levels were associated with the rapid progression. INTERPRETATION: There is considerable individual variation in the rate of progression from appearance of type 1 diabetes-specific autoantibodies to clinical disease. Patient stratification and prediction of disease progression can help in developing more personalised therapeutic strategies for different disease endotypes. FUNDING: A full list of funding bodies can be found under Acknowledgments.