Development and verification of an immune-related gene pairs prognostic signature in ovarian cancer

MUL1 identified as mitochondria-linked biomarker promoting cisplatin resistance in OC cells

Ovarian cancer (OC) ranks among the prevalent tumors affecting the female reproductive system. The aim of this study was to evaluate mitochondria-associated platinum resistance genes using organoid models. Univariate Cox regression, LASSO and multivariate Cox regression analyses were performed on The Cancer Genome Atlas (TCGA) database to construct 2-gene prognostic signature (MUL1 and SSBP1), and GSE26712 dataset was used for external validation. In addition, the relationship between MUL1 and platinum resistance was examined by organoid culture, lentiviral transduction, CCK8 assay, and Western blot. The results showed that patients in the high-risk group exhibited significantly worse OS (P = 0.002, P = 0.017). Drug sensitivity analysis revealed that platinum resistance increased with the upregulation of MUL1 expression (Cor = 0.5154, P = 0.02). Our experimental findings demonstrated that knockout of the MUL1 gene significantly increased apoptosis and enhanced the sensitivity of the OC cell line A2780 to cisplatin. Through this study, we have provided strong evidence for further research on prognostic risk factors and individualized treatment in OC patients, and provided new insights into addressing platinum resistance in OC.

An immune-related gene pair signature predicts the prognosis and immunotherapeutic response in glioblastoma

Background

Glioblastoma (GBM) has the feature of aggressive growth and high rates of recurrence. Immunotherapy was not included in standard therapy for GBM due to lacking the predictive biomarkers. In the present study, we performed an immune-related gene pair (IRGP) signature to predict the prognosis and immunotherapy response of GBM.

Methods

A total of 160 GBM patients from TCGA were included. ssGSEA was conducted to evaluate the immune infiltration level. Univariate Cox, LASSO regression analysis, ROC analysis, and Kaplan-Meier survival analysis were applied to construct and evaluate the risk model. Moreover, the association between immune infiltration and the risk score was assessed. Finally, the expression of immune checkpoints between different risk groups was explored.

Results

According to the normal/tumor, high-/low-immunity group, we identified 125 differentially expressed immune-related genes. Subsequently, a prognostic model including 22 IRGPs was established. The area under the ROC curve to predict 1, 3, and 5-year was 0.811, 0.958, and 0.99 respectively. According to the optimal cut-off value of the 3-year ROC curve, patients were classified into high- and low-risk groups. The Kaplan-Meier analysis result indicated that patients in the low-risk group have longer survival time. The risk score was an independent prognostic predictor (P < 0.001). Moreover, PDCD1 was positively correlated with the risk score (P < 0.01). We also found that patients with high PDCD1 expression had worse survival.

Conclusions

The IRGP signature was built to predict the prognosis of GBM patients. This signature can serve as a tool to predict the response to immunotherapy in GBM.

Based on the prognosis model of immunogenes, the prognosis model was constructed to predict the invasion of immune genes and immune cells related to primary liver cancer and its experimental validation

Background

Primary liver cancer (PLC) is a prevalent malignancy of the digestive system characterized by insidious symptom onset and a generally poor prognosis. Recent studies have highlighted a significant correlation between the initiation and prognosis of liver cancer and the immune function of PLC patients.

Purpose

Revealing the expression of PLC-related immune genes and the characteristics of immune cell infiltration provides assistance for the analysis of clinical pathological parameters and prognosis of PLC patients.

Methods

PLC-related differentially expressed genes (DEGs) with a median absolute deviation (MAD > 0.5) were identified from TCGA and GEO databases. These DEGs were intersected with immune-related genes (IRGs) from the ImmPort database to obtain PLC-related IRGs. The method of constructing a prognostic model through immune-related gene pairs (IRGPs) is used to obtain IRGPs and conduct the selection of central immune genes. The central immune genes obtained from the selection of IRGPs are validated in PLC. Subsequently, the relative proportions of 22 types of immune cells in different immune risk groups are evaluated, and the differential characteristics of PLC-related immune cells are verified through animal experiments.

Results

Through database screening and the construction of an IRGP prognosis model, 84 pairs of IRGPs (P < 0.001) were ultimately obtained. Analysis of these 84 IRGPs revealed 11 central immune genes related to PLC, showing differential expression in liver cancer tissues compared to normal liver tissues. Results from the CiberSort platform indicate differential expression of immune cells such as naive B cells, macrophages, and neutrophils in different immune risk groups. Animal experiments demonstrated altered immune cell proportions in H22 tumor-bearing mice, validating findings from peripheral blood and spleen homogenate analyses.

Conclusion

Our study successfully predicted and validated PLC-related IRGs and immune cells, suggesting their potential as prognostic indicators and therapeutic targets for PLC.

Personalization of Therapy in High-Grade Serous Tubo-Ovarian Cancer-The Possibility or the Necessity?

High-grade serous tubo-ovarian cancer (HGSTOC) is the most lethal tumor of the female genital tract. The foregoing therapy consists of cytoreduction followed by standard platinum/taxane chemotherapy; alternatively, for primary unresectable tumors, neo-adjuvant platinum/taxane chemotherapy followed by delayed interval cytoreduction. In patients with suboptimal surgery or advanced disease, different forms of targeted therapy have been accepted or tested in clinical trials. Studies on HGSTOC discovered its genetic and proteomic heterogeneity, epigenetic regulation, and the role of the tumor microenvironment. These findings turned attention to the fact that there are several distinct primary tumor subtypes of HGSTOC and the unique biology of primary, metastatic, and recurrent tumors may result in a differential drug response. This results in both chemo-refractoriness of some primary tumors and, what is significantly more frequent and destructive, secondary chemo-resistance of metastatic and recurrent HGSTOC tumors. Treatment possibilities for platinum-resistant disease include several chemotherapeutics with moderate activity and different targeted drugs with difficult tolerable effects. Therefore, the question appears as to why different subtypes of ovarian cancer are predominantly treated based on the same therapeutic schemes and not in an individualized way, adjusted to the biology of a specific tumor subtype and temporal moment of the disease. The paper reviews the genomic, mutational, and epigenetic signatures of HGSTOC subtypes and the tumor microenvironment. The clinical trials on personalized therapy and the overall results of a new, comprehensive approach to personalized therapy for ovarian cancer have been presented and discussed.

Novel immune-related gene signature for risk stratification and prognosis prediction in ovarian cancer

BACKGROUND

The immune system played a multifaceted role in ovarian cancer (OC) and was a significant mediator of ovarian carcinogenesis. Various immune cells and immune gene products played an integrated role in ovarian cancer (OC) progression, proved the significance of the immune microenvironment in prognosis. Therefore, we aimed to establish and validate an immune gene prognostic signature for OC patients' prognosis prediction.

METHODS

Differently expressed Immune-related genes (DEIRGs) were identified in 428 OC and 77 normal ovary tissue specimens from 9 independent GEO datasets. The Cancer Genome Atlas (TCGA) cohort was used as a training cohort, Univariate Cox analysis was used to identify prognostic DEIRGs in TCGA cohort. Then, an immune gene-based risk model for prognosis prediction was constructed using the LASSO regression analysis, and validated the accuracy and stability of the model in 374 and 93 OC patients in TCGA training cohort and International Cancer Genome Consortium (ICGC) validation cohort respectively. Finally, the correlation among risk score model, clinicopathological parameters, and immune cell infiltration were analyzed.

RESULTS

Five DEIRGs were identified to establish the immune gene signature and divided OC patients into the low- and high-risk groups. In TCGA and ICGC datasets, patients in the low-risk group showed a substantially higher survival rate than high-risk group. Receiver operating characteristic (ROC) curves, t-distributed stochastic neighbor embedding (t-SNE) analysis and principal component analysis (PCA) showed the good performance of the risk model. Clinicopathological correlation analysis proved the risk score model could serve as an independent prognostic factor in 2 independent datasets.

CONCLUSIONS

The prognostic model based on immune-related genes can function as a superior prognostic indicator for OC patients, which could provide evidence for individualized treatment and clinical decision making.

A novel defined risk signature of endoplasmic reticulum stress-related genes for predicting the prognosis and immune infiltration status of ovarian cancer

Endoplasmic reticulum (ER) stress, as an emerging hallmark feature of cancer, has a considerable impact on cell proliferation, metastasis, invasion, and chemotherapy resistance. Ovarian cancer (OvCa) is one of the leading causes of cancer-related mortality across the world due to the late stage of disease at diagnosis. Studies have explored the influence of ER stress on OvCa in recent years, while the predictive role of ER stress-related genes in OvCa prognosis remains unexplored. Here, we enrolled 552 cases of ER stress-related genes involved in OvCa from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts for the screening of prognosis-related genes. The least absolute shrinkage and selection operator (LASSO) regression was applied to establish an ER stress-related risk signature based on the TCGA cohort. A seven-gene signature revealed a favorable predictive efficacy for the TCGA, International Cancer Genome Consortium (ICGC), and another GEO cohort (P<0.001, P<0.001, and P=0.04, respectively). Moreover, functional annotation indicated that this signature was enriched in cellular response and senescence, cytokines interaction, as well as multiple immune-associated terms. The immune infiltration profiles further delineated an immunologic unresponsive status in the high-risk group. In conclusion, ER stress-related genes are vital factors predicting the prognosis of OvCa, and possess great application potential in the clinic.

Overexpression of LILRA2 indicated poor prognosis of ovarian carcinoma: A new potential biomarker and therapeutic target

OBJECTIVE

This study aimed to assess the role of leukocyte immunoglobulin-like receptor A2 (LILRA2) in ovarian carcinoma (OC) oncogenesis and prognosis.

MATERIALS AND METHODS

Using the Cancer Genome Atlas, Genotype-Tissue Expression, and Gene Expression Omnibus databases, the association between clinicopathological profiles and LILRA2 expression was investigated using logistic regression analysis. Kaplan-Meier analysis, Cox regression analysis, and column plots predicted the clinical outcomes of patients with OC and determine the predictive value of LILRA2. The biological functions of LILRA2 were assessed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. We used single-sample Gene Set Enrichment Analysis to investigate the relationship between immune cell infiltration and LILRA2 expression.

RESULTS

LILRA2 expression in OC tumors was significantly higher than in normal tissue (P < 0.05). The high LILRA2 expression in OC was correlated with lymphatic invasion (P = 0.014). The results showed consistency indices of 0.611 [95% confidence interval (CI), 0.572-0.649] and 0.623 (95% CI, 0.584-0.663) for the overall and disease-specific survival nomograms, respectively. Cox regression analysis showed that LILRA2 was an independent risk factor for overall survival (hazard ratio [HR], 1.511; P = 0.002) and disease-specific survival (HR, 1.537; P = 0.003). Functional annotation revealed enrichment with immunoglobulin-corresponding pathways when LILRA2 expression was high.

CONCLUSION

By evaluating gene expression profiles, we demonstrated that LILRA2 has considerable potential to act as a therapeutic target and prognostic biomarker in OC.

Identification of an Individualized Prognostic Biomarker for Serous Ovarian Cancer: A Qualitative Model

Serous ovarian cancer is the most common type of ovarian epithelial cancer and usually has a poor prognosis. The objective of this study was to construct an individualized prognostic model for predicting overall survival in serous ovarian cancer. Based on the relative expression orderings (Ea > Eb/Ea ≤ Eb) of gene pairs closely associated with serous ovarian prognosis, we tried constructing a potential individualized qualitative biomarker by the greedy algorithm and evaluated the performance in independent validation datasets. We constructed a prognostic biomarker consisting of 20 gene pairs (SOV-P20). The overall survival between high- and low-risk groups stratified by SOV-P20 was statistically significantly different in the training and independent validation datasets from other platforms (p < 0.05, Wilcoxon test). The average area under the curve (AUC) values of the training and three validation datasets were 0.756, 0.590, 0.630, and 0.680, respectively. The distribution of most immune cells between high- and low-risk groups was quite different (p < 0.001, Wilcoxon test). The low-risk patients tended to show significantly better tumor response to chemotherapy than the high-risk patients (p < 0.05, Fisher’s exact test). SOV-P20 achieved the highest mean index of concordance (C-index) (0.624) compared with the other seven existing prognostic signatures (ranging from 0.511 to 0.619). SOV-P20 is a promising prognostic biomarker for serous ovarian cancer, which will be applicable for clinical predictive risk assessment.

Integration of Transcriptome and Epigenome to Identify and Develop Prognostic Markers for Ovarian Cancer

DNA methylation is a widely researched epigenetic modification. It is associated with the occurrence and development of cancer and has helped evaluate patients' prognoses. However, most existing DNA methylation prognosis models have not simultaneously considered the changes of the downstream transcriptome. Methods. The RNA-Sequencing data and DNA methylation omics data of ovarian cancer patients were downloaded from The Cancer Genome Atlas (TCGA) database. The Consensus Cluster Plus algorithm was used to construct the methylated molecular subtypes of the ovary. Lasso regression was employed to build a multi-gene signature. An independent data set was applied to verify the prognostic value of the signature. The Gene Set Variation Analysis (GSVA) was used to carry out the enrichment analysis of the pathways linked to the gene signature. The IMvigor 210 cohort was used to explore the predictive efficacy of the gene signature for immunotherapy response. Results. We distinguished ovarian cancer samples into two subtypes with different prognosis, based on the omics data of DNA methylation. Differentially expressed genes and enrichment analysis among subtypes indicated that DNA methylation was related to fatty acid metabolism and the extracellular matrix (ECM)-receptor. Furthermore, we constructed an 8-gene signature, which proved to be efficient and stable in predicting prognostics in ovarian cancer patients with different data sets and distinctive pathological characteristics. Finally, the 8-gene signature could predict patients' responses to immunotherapy. The polymerase chain reaction experiment was further used to verify the expression of 8 genes. Conclusion. We analyzed the prognostic value of the related genes of methylation in ovarian cancer. The 8-gene signature predicted the prognosis and immunotherapy response of ovarian cancer patients well and is expected to be valuable in clinical application.

A risk model of gene signatures for predicting platinum response and survival in ovarian cancer

Background

Ovarian cancer (OC) is the deadliest tumor in the female reproductive tract. And increased resistance to platinum-based chemotherapy represents the major obstacle in the treatment of OC currently. Robust and accurate gene expression models are crucial tools in distinguishing platinum therapy response and evaluating the prognosis of OC patients.

Methods

In this study, 230 samples from The Cancer Genome Atlas (TCGA) OV dataset were subjected to mRNA expression profiling, single nucleotide polymorphism (SNP), and copy number variation (CNV) analysis comprehensively to screen out the differentially expressed genes (DEGs). An SVM classifier and a prognostic model were constructed using the Random Forest algorithm and LASSO Cox regression model respectively via R. The Gene Expression Omnibus (GEO) database was applied as the validation set.

Results

Forty-eight differentially expressed genes (DEGs) were figured out through integrated analysis of gene expression, single nucleotide polymorphism (SNP), and copy number variation (CNV) data. A 10-gene classifier was constructed which could discriminate platinum-sensitive samples precisely with an AUC of 0.971 in the training set and of 0.926 in the GEO dataset (GSE638855). In addition, 8 optimal genes were further selected to construct the prognostic risk model whose predictions were consistent with the actual survival outcomes in the training cohort (p = 9.613e-05) and validated in GSE638855 (p = 0.04862). PNLDC1, SLC5A1, and SYNM were then identified as hub genes that were associated with both platinum response status and prognosis, which was further validated by the Fudan University Shanghai cancer center (FUSCC) cohort.

Conclusion

These findings reveal a specific risk model that could serve as effective biomarkers to identify patients’ platinum response status and predict survival outcomes for OC patients. PNLDC1, SLC5A1, and SYNM are the hub genes that may serve as potential biomarkers in OC treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-022-00969-3.

Prognostic immunologic signatures in epithelial ovarian cancer

Epithelial Ovarian Cancer (EOC) is a deadly gynecologic malignancy in which patients frequently develop recurrent disease following initial platinum-taxane chemotherapy. Analogous to many other cancer subtypes, EOC clinical trials have centered upon immunotherapeutic approaches, most notably programmed cell death 1 (PD-1) inhibitors. While response rates to these immunotherapies in EOC patients have been low, evidence suggests that ovarian tumors are immunogenic and that immune-related genomic profiles can serve as prognostic markers. This review will discuss recent advances in the development of immune-based prognostic signatures in EOC that predict patient clinical outcomes, as well as emphasize specific research areas that need to be addressed to drive this field forward.

Identification of an Innate Immune-Related Prognostic Signature in Early-Stage Lung Squamous Cell Carcinoma

Background

Early-stage lung squamous cell carcinoma (LUSC) progression is accompanied by changes in immune microenvironments and the expression of immune-related genes (IRGs). Identifying innate IRGs associated with prognosis may improve treatment and reveal new immunotherapeutic targets.

Methods

Gene expression profiles and clinical data of early-stage LUSC patients were obtained from the Gene Expression Omnibus and The Cancer Genome Atlas databases and IRGs from the InnateDB database. Univariate and multivariate Cox regression and LASSO regression analyses were performed to identify an innate IRG signature model prognostic in patients with early-stage LUSC. The predictive ability of this model was assessed by time-dependent receiver operator characteristic curve analysis, with the independence of the model-determined risk score assessed by univariate and multivariate Cox regression analyses. Overall survival (OS) in early-stage LUSC patients was assessed using a nomogram and decision curve analysis (DCA). Functional and biological pathways were determined by gene set enrichment analysis, and differences in biological functions and immune microenvironments between the high- and low-risk groups were assessed by ESTIMATE and the CIBERSORT algorithm.

Results

A signature involving six IRGs (SREBF2, GP2, BMX, NR1H4, DDX41, and GOPC) was prognostic of OS. Samples were divided into high- and low-risk groups based on median risk scores. OS was significantly shorter in the high-risk than in the low-risk group in the training (P < 0.001), GEO validation (P = 0.00021) and TCGA validation (P = 0.034) cohorts. Multivariate Cox regression analysis showed that risk score was an independent risk factor for OS, with the combination of risk score and T stage being optimally predictive of clinical benefit. GSEA, ESTIMATE, and the CIBERSORT algorithm showed that immune cell infiltration was higher and immune-related pathways were more strongly expressed in the low-risk group.

Conclusion

A signature that includes these six innate IRGs may predict prognosis in patients with early-stage LUSC.

OLFML2B Is a Robust Prognostic Biomarker in Bladder Cancer Through Genome-Wide Screening: A Study Based on Seven Cohorts

Background

Bladder cancer lacks useful and robust prognostic markers to stratify patients at risk. Our study is to identify a robust prognostic marker for bladder cancer.

Methods

The transcriptome and clinical data of bladder cancer were downloaded from multiple databases. We searched for genes with robust prognosis by Kaplan-Meier analysis of the whole genome. CIBERSORT and TIMER algorithm was used to calculate the degree of immune cell infiltration.

Results

We identified OLFML2B as a robust prognostic marker for bladder cancer in five cohorts. Kaplan-Meier analysis showed that patients with a high level of OLFML2B expression had a poor prognosis. The expression of OLFML2B increased with the increase of stage and grade. We found that patients with high expression of OLFML2B still had a poor prognosis in two small bladder cancer cohorts. OLFML2B also has the prognostic ability in ten other tumors, and the prognosis is poor in high expression. The correlation analysis between OLFML2B and immune cells showed that it was positively correlated with the degree of macrophage infiltration and highly co-expressed with tumor-associated macrophage markers. Finally, the Wound-healing assay and Colony formation assay results showed that the migration and proliferation ability of bladder cancer cell lines decreased after the knockdown of OLFML2B.

Conclusions

In summary, OLFML2B is a robust risk prognostic marker, and it can help patients with bladder cancer improve individualized treatment.

An Oxidative Stress-Related Gene Pair (CCNB1/PKD1), Competitive Endogenous RNAs, and Immune-Infiltration Patterns Potentially Regulate Intervertebral Disc Degeneration Development

Oxidative stress (OS) irreversibly affects the pathogenesis of intervertebral disc degeneration (IDD). Certain non-coding RNAs act as competitive endogenous RNAs (ceRNAs) that regulate IDD progression. Analyzing the signatures of oxidative stress-related gene (OSRG) pairs and regulatory ceRNA mechanisms and immune-infiltration patterns associated with IDD may enable researchers to distinguish IDD and reveal the underlying mechanisms. In this study, OSRGs were downloaded and identified using the Gene Expression Omnibus database. Functional-enrichment analysis revealed the involvement of oxidative stress-related pathways and processes, and a ceRNA network was generated. Differentially expressed oxidative stress-related genes (De-OSRGs) were used to construct De-OSRG pairs, which were screened, and candidate De-OSRG pairs were identified. Immune cell-related gene pairs were selected via immune-infiltration analysis. A potential long non-coding RNA-microRNA-mRNA axis was determined, and clinical values were assessed. Eighteen De-OSRGs were identified that were primarily related to intricate signal-transduction pathways, apoptosis-related biological processes, and multiple kinase-related molecular functions. A ceRNA network consisting of 653 long non-coding RNA-microRNA links and 42 mRNA-miRNA links was constructed. Three candidate De-OSRG pairs were screened out from 13 De-OSRG pairs. The abundances of resting memory CD4+ T cells, resting dendritic cells, and CD8+ T cells differed between the control and IDD groups. CD8+ T cell infiltration correlated negatively with cyclin B1 (CCNB1) expression and positively with protein kinase D1 (PKD1) expression. CCNB1-PKD1 was the only pair that was differentially expressed in IDD, was correlated with CD8+ T cells, and displayed better predictive accuracy compared to individual genes. The PKD1-miR-20b-5p-AP000797 and CCNB1-miR-212-3p-AC079834 axes may regulate IDD. Our findings indicate that the OSRG pair CCNB1-PKD1, which regulates oxidative stress during IDD development, is a robust signature for identifying IDD. This OSRG pair and increased infiltration of CD8+ T cells, which play important roles in IDD, were functionally associated. Thus, the OSRG pair CCNB1-PKD1 is promising target for treating IDD.

Identification of a glycolysis-related gene signature for survival prediction of ovarian cancer patients

Background

Ovarian cancer (OV) is deemed the most lethal gynecological cancer in women. The aim of this study was to construct an effective gene prognostic model for predicting overall survival (OS) in patients with OV.

Methods

The expression profiles of glycolysis‐related genes (GRGs) and clinical data of patients with OV were extracted from The Cancer Genome Atlas (TCGA) database. Univariate, multivariate, and least absolute shrinkage and selection operator Cox regression analyses were conducted, and a prognostic signature based on GRGs was constructed. The predictive ability of the signature was analyzed using training and test sets.

Results

A gene risk signature based on nine GRGs (ISG20, CITED2, PYGB, IRS2, ANGPTL4, TGFBI, LHX9, PC, and DDIT4) was identified to predict the survival outcome of patients with OV. The signature showed a good prognostic ability for OV, particularly high‐grade OV, in the TCGA dataset, with areas under the curve (AUC) of 0.709 and 0.762 for 3‐ and 5‐year survival, respectively. Similar results were found in the test sets, and the AUCs of 3‐, 5‐year OS were 0.714 and 0.772 in the combined test set. And our signature was an independent prognostic factor. Moreover, a nomogram combining the prediction model and clinical factors was developed.

Conclusion

Our study established a nine‐GRG risk model and nomogram to better predict OS in patients with OV. The risk model represents a promising and independent prognostic predictor for patients with OV. Moreover, our study on GRGs could offer guidance for the elucidation of underlying mechanisms in future studies.

Development and verification of an immune-related gene pairs prognostic signature in ovarian cancer

Ovarian cancer (OV) is the most common gynaecological cancer worldwide. Immunotherapy has recently been proven to be an effective treatment strategy. The work here attempts to produce a prognostic immune-related gene pair (IRGP) signature to estimate OV patient survival. The Gene Expression Omnibus (GEO) and Cancer Genome Atlas (TCGA) databases provided the genetic expression profiles and clinical data of OV patients. Based on the InnateDB database and the least absolute shrinkage and selection operator (LASSO) regression model, we first identified a 17-IRGP signature associated with survival. The average area under the curve (AUC) values of the training, validation, and all TCGA sets were 0.869, 0.712, and 0.778, respectively. The 17-IRGP signature noticeably split patients into high- and low-risk groups with different prognostic outcomes. As suggested by a functional study, some biological pathways, including the Toll-like receptor and chemokine signalling pathways, were significantly negatively correlated with risk scores; however, pathways such as the p53 and apoptosis signalling pathways had a positive correlation. Moreover, tumour stage III, IV, grade G1/G2, and G3/G4 samples had significant differences in risk scores. In conclusion, an effective 17-IRGP signature was produced to predict prognostic outcomes in OV, providing new insights into immunological biomarkers.