Preserving biological heterogeneity with a permuted surrogate variable analysis for genomics batch correction

Neurobiology of attention-deficit hyperactivity disorder: historical challenges and emerging frontiers

Extensive investigations spanning multiple levels of inquiry, from genetic to behavioural studies, have sought to unravel the mechanistic foundations of attention-deficit hyperactivity disorder (ADHD), with the aspiration of developing efficacious treatments for this condition. Despite these efforts, the pathogenesis of ADHD remains elusive. In this Review, we reflect on what has been learned about ADHD while also providing a framework that may serve as a roadmap for future investigations. We emphasize that ADHD is a highly heterogeneous disorder with multiple aetiologies that necessitates a multifactorial dimensional phenotype, rather than a fixed dichotomous conceptualization. We highlight new findings that suggest a more brain-wide, 'global' view of the disorder, rather than the traditional localizationist framework, which asserts that a limited set of brain regions or networks underlie ADHD. Last, we underscore how underpowered studies that have aimed to associate neurobiology with ADHD phenotypes have long precluded the field from making progress. However, a new age of ADHD research with refined phenotypes, advanced methods, creative study designs and adequately powered investigations is beginning to put the field on a good footing. Indeed, the field is at a promising juncture to advance the neurobiological understanding of ADHD and fulfil the promise of clinical utility.

Assessing and mitigating batch effects in large-scale omics studies

Batch effects in omics data are notoriously common technical variations unrelated to study objectives, and may result in misleading outcomes if uncorrected, or hinder biomedical discovery if over-corrected. Assessing and mitigating batch effects is crucial for ensuring the reliability and reproducibility of omics data and minimizing the impact of technical variations on biological interpretation. In this review, we highlight the profound negative impact of batch effects and the urgent need to address this challenging problem in large-scale omics studies. We summarize potential sources of batch effects, current progress in evaluating and correcting them, and consortium efforts aiming to tackle them.

Thinking points for effective batch correction on biomedical data

Batch effects introduce significant variability into high-dimensional data, complicating accurate analysis and leading to potentially misleading conclusions if not adequately addressed. Despite technological and algorithmic advancements in biomedical research, effectively managing batch effects remains a complex challenge requiring comprehensive considerations. This paper underscores the necessity of a flexible and holistic approach for selecting batch effect correction algorithms (BECAs), advocating for proper BECA evaluations and consideration of artificial intelligence-based strategies. We also discuss key challenges in batch effect correction, including the importance of uncovering hidden batch factors and understanding the impact of design imbalance, missing values, and aggressive correction. Our aim is to provide researchers with a robust framework for effective batch effects management and enhancing the reliability of high-dimensional data analyses.

Exhaustive Exploitation of Nature-Inspired Computation for Cancer Screening in an Ensemble Manner

Accurate screening of cancer types is crucial for effective cancer detection and precise treatment selection. However, the association between gene expression profiles and tumors is often limited to a small number of biomarker genes. While computational methods using nature-inspired algorithms have shown promise in selecting predictive genes, existing techniques are limited by inefficient search and poor generalization across diverse datasets. This study presents a framework termed Evolutionary Optimized Diverse Ensemble Learning (EODE) to improve ensemble learning for cancer classification from gene expression data. The EODE methodology combines an intelligent grey wolf optimization algorithm for selective feature space reduction, guided random injection modeling for ensemble diversity enhancement, and subset model optimization for synergistic classifier combinations. Extensive experiments were conducted across 35 gene expression benchmark datasets encompassing varied cancer types. Results demonstrated that EODE obtained significantly improved screening accuracy over individual and conventionally aggregated models. The integrated optimization of advanced feature selection, directed specialized modeling, and cooperative classifier ensembles helps address key challenges in current nature-inspired approaches. This provides an effective framework for robust and generalized ensemble learning with gene expression biomarkers.

Characterizing mitochondrial features in osteoarthritis through integrative multi-omics and machine learning analysis

Purpose

Osteoarthritis (OA) stands as the most prevalent joint disorder. Mitochondrial dysfunction has been linked to the pathogenesis of OA. The main goal of this study is to uncover the pivotal role of mitochondria in the mechanisms driving OA development.

Materials and methods

We acquired seven bulk RNA-seq datasets from the Gene Expression Omnibus (GEO) database and examined the expression levels of differentially expressed genes related to mitochondria in OA. We utilized single-sample gene set enrichment analysis (ssGSEA), gene set enrichment analysis (GSEA), and weighted gene co-expression network analysis (WGCNA) analyses to explore the functional mechanisms associated with these genes. Seven machine learning algorithms were utilized to identify hub mitochondria-related genes and develop a predictive model. Further analyses included pathway enrichment, immune infiltration, gene-disease relationships, and mRNA-miRNA network construction based on these hub mitochondria-related genes. genome-wide association studies (GWAS) analysis was performed using the Gene Atlas database. GSEA, gene set variation analysis (GSVA), protein pathway analysis, and WGCNA were employed to investigate relevant pathways in subtypes. The Harmonizome database was employed to analyze the expression of hub mitochondria-related genes across various human tissues. Single-cell data analysis was conducted to examine patterns of gene expression distribution and pseudo-temporal changes. Additionally, The real-time polymerase chain reaction (RT-PCR) was used to validate the expression of these hub mitochondria-related genes.

Results

In OA, the mitochondria-related pathway was significantly activated. Nine hub mitochondria-related genes (SIRT4, DNAJC15, NFS1, FKBP8, SLC25A37, CARS2, MTHFD2, ETFDH, and PDK4) were identified. They constructed predictive models with good ability to predict OA. These genes are primarily associated with macrophages. Unsupervised consensus clustering identified two mitochondria-associated isoforms that are primarily associated with metabolism. Single-cell analysis showed that they were all expressed in single cells and varied with cell differentiation. RT-PCR showed that they were all significantly expressed in OA.

Conclusion

SIRT4, DNAJC15, NFS1, FKBP8, SLC25A37, CARS2, MTHFD2, ETFDH, and PDK4 are potential mitochondrial target genes for studying OA. The classification of mitochondria-associated isoforms could help to personalize treatment for OA patients.

Machine learning-derived immunosenescence index for predicting outcome and drug sensitivity in patients with skin cutaneous melanoma

The functions of immunosenescence are closely related to skin cutaneous melanoma (SKCM). The aim of this study is to uncover the characteristics of immunosenescence index (ISI) to identify novel biomarkers and potential targets for treatment. Firstly, integrated bioinformatics analysis was carried out to identify risk prognostic genes, and their expression and prognostic value were evaluated. Then, we used the computational algorithm to estimate ISI. Finally, the distribution characteristics and clinical significance of ISI in SKCM by using multi-omics analysis. Patients with a lower ISI had a favorable survival rate, lower chromosomal instability, lower somatic copy-number alterations, lower somatic mutations, higher immune infiltration, and sensitive to immunotherapy. The ISI exhibited robust, which was validated in multiple datasets. Besides, the ISI is more effective than other published signatures in predicting survival outcomes for patients with SKCM. Single-cell analysis revealed higher ISI was specifically expressed in monocytes, and correlates with the differentiation fate of monocytes in SKCM. Besides, individuals exhibiting elevated ISI levels could potentially receive advantages from chemotherapy, and promising compounds with the potential to target high ISI were recognized. The ISI model is a valuable tool in categorizing SKCM patients based on their prognosis, gene mutation signatures, and response to immunotherapy.

Research on the impact of pilot free trade zones on urban green development: A case study based on the Yangtze River Economic Belt in China

Green development is an important component of China's new development concept. Pilot Free Trade Zones (PFTZs), as "experimental fields" for promoting reform, deepening opening-up, and raising the level of an open economy, are important open areas for China to promote green development. However, existing related research is not extensive. This article takes PFTZs as quasi-natural experiments, with the Yangtze River Economic Belt (YREB) as the research area. Based on urban panel data from 2006 to 2020, using multi-period differences-in-differences and spatial differences-in-differences models, it explores the impact effects of PFTZs on urban green development and their potential mechanisms. The research findings indicate: (1) Overall, PFTZs significantly promote urban green development, with variations in impact effects due to different batches and locations of establishment. (2) Mechanism tests show that PFTZs mainly promote urban green development by stimulating technological innovation, industrial upgrading, and reducing government intervention. (3) From the perspective of spatial spillover effects, the establishment of PFTZs not only promotes the green development process in the host cities but also has a promoting effect on the green development of surrounding cities.

Identification of HDAC9 and ARRDC4 as potential biomarkers and targets for treatment of type 2 diabetes

We aimed to identify the key potential insulin resistance (IR)-related genes and investigate their correlation with immune cell infiltration in type 2 diabetes (T2D). The GSE78721 dataset (68 diabetic patients and 62 controls) was downloaded from the Gene Expression Omnibus database and utilized for single-sample gene set enrichment analysis. IR-related genes were obtained from the Comparative Toxicology Genetics Database, and the final IR-differentially expressed genes (DEGs) were screened by intersecting with the DEGs obtained from the GSE78721 datasets. Functional enrichment analysis was performed, and the networks of the target gene with microRNA, transcription factor, and drug were constructed. Hub genes were identified based on a protein-protein interaction network. Least absolute shrinkage and selection operator regression and Random Forest and Boruta analysis were combined to screen diagnostic biomarkers in T2D, which were validated using the GSE76894 (19 diabetic patients and 84 controls) and GSE9006 (12 diabetic patients and 24 controls) datasets. Quantitative real-time polymerase chain reaction was performed to validate the biomarker expression in IR mice and control mice. In addition, infiltration of immune cells in T2D and their correlation with the identified markers were computed using CIBERSORT. We identified differential immune gene set regulatory T-cells in the GSE78721 dataset, and T2D samples were assigned into three clusters based on immune infiltration. A total of 2094 IR-DEGs were primarily enriched in response to endoplasmic reticulum stress. Importantly, HDAC9 and ARRDC4 were identified as markers of T2D and associated with different levels of immune cell infiltration. HDAC9 mRNA level were higher in the IR mice than in control mice, while ARRDC4 showed the opposite trend. In summary, we discovered potential vital biomarkers that contribute to immune cell infiltration associated with IR, which offers a new sight of immunotherapy for T2D.

Integration of bulk RNA sequencing to reveal protein arginine methylation regulators have a good prognostic value in immunotherapy to treat lung adenocarcinoma

Background

Given the differential expression and biological functions of protein arginine methylation (PAM) regulators in lung adenocarcinoma (LUAD), it may be of great value in the diagnosis, prognosis, and treatment of LUAD. However, the expression and function of PAM regulators in LUAD and its relationship with prognosis are unclear.

Methods

8 datasets including 1798 LUAD patients were selected. During the bioinformatic study in LUAD, we performed (i) consensus clustering to identify clusters based on 9 PAM regulators related expression profile data, (ii) to identify hub genes between the 2 clusters, (iii) principal component analysis to construct a PAM.score based on above genes, and (iv) evaluation of the effect of PAM.score on the deconstruction of tumor microenvironment and guidance of immunotherapy.

Results

We identified two different clusters and a robust and clinically practicable prognostic scoring system. Meanwhile, a higher PAM.score subgroup showed poorer prognosis, and was validated by multiple cohorts. Its prognostic effect was validated by ROC (Receiver operating characteristic curve) curve and found to have a relatively good prediction efficacy. High PAM.score group exhibited lower immune score, which associated with an immunosuppressive microenvironment in LUAD. Finally, patients exhibiting a lower PAM.score presented noteworthy therapeutic benefits and clinical advantages.

Conclusion

Our PAM.score model can help clinicians to select personalized therapy for LUAD patients, and PAM.score may act a part in the development of LUAD.

Identification and analysis of diverse cell death patterns in diabetic kidney disease using microarray-based transcriptome profiling and single-nucleus RNA sequencing

BACKGROUND

Diabetic kidney disease (DKD) is the most lethal complication of diabetes. Diverse programmed cell death (PCD) has emerged as a crucial disease phenotype that has the potential to serve as an indicator of renal function decline and can be used as a target for researching drugs for DKD.

METHODS

Microarray-based transcriptome profiling and single-nucleus transcriptome sequencing (snRNA-seq) related to DKD were retrieved from the Gene Expression Omnibus (GEO) database. 13 PCD-related genes (including alkaliptosis, apoptosis, autophagy-dependent cell death, cuproptosis, disulfidptosis, entotic cell death, ferroptosis, lysosome-dependent cell death, necroptosis, netotic cell death, oxeiptosis, parthanatos, and pyroptosis) were obtained from various public databases and reviews. The gene set variation analysis (GSVA) analysis was used to explore the pathway activity of these 13 PCDs in DKD, and the pathway activity of these PCDs in different renal cells was studied based on DKD-related snRNA-seq data. To identify the core PCDs that play a significant role in DKD, we analyzed the relationships between different types of PCD and immune infiltration, fibrosis-related gene expression levels, glomerular filtration rate (GFR), and diagnostic efficiency in DKD. Using the Weighted Gene Co-expression Network Analysis (WGCNA) algorithm, we screened for core death genes among the core PCDs and constructed a cell death-related signature (CDS) risk score based on the Least Absolute Shrinkage and Selection Operator (LASSO). Finally, we validated the predictive performance of the CDS risk score in an independent validation set.

RESULTS

We identified 4 core PCD pathways, namely entotic cell death, apoptosis, necroptosis, and pyroptosis in DKD, and further applied the WGCNA algorithm to screen 4 core death genes (CASP1, CYBB, PLA2G4A, and CTSS) and constructed a CDS risk score based on these genes. The CDS risk score demonstrated high diagnostic efficiency for DKD patients, and those with higher scores had higher levels of immune cell infiltration and poorer GFR.

CONCLUSION

Our study sheds light on the fact that multiple PCDs contribute to the progression of DKD, highlighting potential therapeutic targets for treating this disease.

A disulfidptosis-related lncRNAs signature in hepatocellular carcinoma: prognostic prediction, tumor immune microenvironment and drug susceptibility

Disulfidptosis, a novel type of programmed cell death, has attracted researchers' attention worldwide. However, the role of disulfidptosis-related lncRNAs (DRLs) in liver hepatocellular carcinoma (LIHC) not yet been studied. We aimed to establish and validate a prognostic signature of DRLs and analyze tumor microenvironment (TME) and drug susceptibility in LIHC patients. RNA sequencing data, mutation data, and clinical data were obtained from the Cancer Genome Atlas Database (TCGA). Lasso algorithm and cox regression analysis were performed to identify a prognostic DRLs signature. Kaplan-Meier curves, principal component analysis (PCA), nomogram and calibration curve, function enrichment, TME, immune dysfunction and exclusion (TIDE), tumor mutation burden (TMB), and drug sensitivity analyses were analyzed. External datasets were used to validate the predictive value of DRLs. qRT-PCR was also used to validate the differential expression of the target lncRNAs in tissue samples and cell lines. We established a prognostic signature for the DRLs (MKLN1-AS and TMCC1-AS1) in LIHC. The signature could divide the LIHC patients into low- and high-risk groups, with the high-risk subgroup associated with a worse prognosis. We observed discrepancies in tumor-infiltrating immune cells, immune function, function enrichment, and TIDE between two risk groups. LIHC patients in the high-risk group were more sensitive to several chemotherapeutic drugs. External datasets, clinical tissue, and cell lines confirmed the expression of MKLN1-AS and TMCC1-AS1 were upregulated in LIHC and associated with a worse prognosis. The novel signature based on the two DRLs provide new insight into LIHC prognostic prediction, TME, and potential therapeutic strategies.

Identification of a 5-Gene Cuproptosis Signature Predicting the Prognosis for Colon Adenocarcinoma Based on WGCNA

Background: Colorectal cancer is a highly aggressive malignant tumor that primarily affects the digestive system. It is frequently diagnosed at an advanced stage. Cuproptosis is a copper-dependent form cell death mechanism, distinct from all other known pathways underlying cell death, tumor progression, prognosis, and immune response. Although the role of cuproptosis in colorectal cancer has been investigated over time, there is still an urgent need to explore new methods and insights to understand its potential function. Methods: The Gene Expression Omnibus and The Cancer Genome Atlas gene expression data were systematically explored to investigate the role of cuproptosis in colon adenocarcinoma. The weighted gene coexpression network analysis was used to construct a gene coexpression network and identify the critical module and cuproptosis-related genes correlated with colon adenocarcinoma prognosis. A cuproptosis-related genes prognostic signature for colon adenocarcinoma was identified and validated. To validate the identified gene signature, quantitative reverse transcription-polymerase chain reaction was performed. Cell proliferation assays were analyzed by CCK8 and cell cycle detection. In addition, reactive oxygen species assay was also analyzed. Results: Five hub cuproptosis-related genes (Dihydrolipoamide S-acetyltransferase, Cyclin-dependent kinase inhibitor 2A, ATOX1, VEGFA, and ULK1) were screened and a prognostic risk model for predicting overall survival was established based on these genes. The model was successfully tested in the validation cohort and the GEPIA database. Colon adenocarcinoma patients were categorized into high-risk and low-risk groups based on risk scores. The study revealed that patients with higher risk scores were more likely to have a poor prognosis. Moreover, Dihydrolipoamide S-acetyltransferase was a tumor suppressor gene that can induce cell death and affected the redox reactions in the colon cancer cell line. Conclusions: These findings suggest that the newly identified 5-gene signature may serve as a more reliable prognostic factor than clinical factors such as age and stage of disease. These findings offer a theoretical foundation for further investigation into potential cuproptosis-related biomarkers for predicting colon adenocarcinoma prognosis in the future.

Comprehensive analysis of the relationship between xanthine oxidoreductase activity and chronic kidney disease

Chronic kidney disease (CKD) is a common disease that seriously endangers human health. However, the potential relationship between xanthine oxidoreductase (XOR) activity and CKD remains unclear. In this study, we used clinical data, CKD datasets from the Gene Expression Omnibus database, and untargeted metabolomics to explain the relationship between XOR activity and CKD. First, XOR activity showed high correlation with the biomarkers of CKD, such as serum creatinine, blood urea nitrogen, uric acid, and estimated glomerular filtration rate. Then, we used least absolute shrinkage and selection operator logical regression algorithm and random forest algorithm to screen CKD molecular markers from differentially expressed genes, and the results of qRT-PCR of XDH, KOX-1, and ROMO1 were in accordance with the results of bioinformatics analyses. In addition, untargeted metabolomics analysis revealed that the purine metabolism pathway was significantly enriched in CKD patients in the simulated models of kidney fibrosis.

GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership

Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.

Machine learning-derived identification of tumor-infiltrating immune cell-related signature for improving prognosis and immunotherapy responses in patients with skin cutaneous melanoma

BACKGROUND

Immunoblockade therapy based on the PD-1 checkpoint has greatly improved the survival rate of patients with skin cutaneous melanoma (SKCM). However, existing anti-PD-1 therapeutic efficacy prediction markers often exhibit a poor situation of poor reliability in identifying potential beneficiary patients in clinical applications, and an ideal biomarker for precision medicine is urgently needed.

METHODS

10 multicenter cohorts including 4 SKCM cohorts and 6 immunotherapy cohorts were selected. Through the analysis of WGCNA, survival analysis, consensus clustering, we screened 36 prognostic genes. Then, ten machine learning algorithms were used to construct a machine learning-derived immune signature (MLDIS). Finally, the independent data sets (GSE22153, GSE54467, GSE59455, and in-house cohort) were used as the verification set, and the ROC index standard was used to evaluate the model.

RESULTS

Based on computing framework, we found that patients with high MLDIS had poor overall survival and has good prediction performance in all cohorts and in-house cohort. It is worth noting that MLDIS performs better in each data set than almost all models which from 51 prognostic signatures for SKCM. Meanwhile, high MLDIS have a positive prognostic impact on patients treated with anti-PD-1 immunotherapy by driving changes in the level of infiltration of immune cells in the tumor microenvironment. Additionally, patients suffering from SKCM with high MLDIS were more sensitive to immunotherapy.

CONCLUSIONS

Our study identified that MLDIS could provide new insights into the prognosis of SKCM and predict the immunotherapy response in patients with SKCM.

Predicting key gene related to immune infiltration and myofibroblast-like valve interstitial cells in patients with calcified aortic valve disease based on bioinformatics analysis

BACKGROUND

Calcified aortic valve disease (CAVD) is the most prevalent valvular disease that can be treated only through valve replacement. We aimed to explore potential biomarkers and the role of immune cell infiltration in CAVD progression through bioinformatics analysis.

METHODS

Differentially ex-pressed genes (DEGs) were screened out based on three microarray datasets: GSE12644, GSE51472 and GSE83453. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to evaluate gene expression differences. Machine learning algorithms and DEGs were used to screen key gene. We used CIBERSORT to evaluate the immune cell infiltration of CAVD and evaluated the correlation between the biomarkers and infiltrating immune cells. We also compared bioinformatics analysis results with the valve interstitial cells (VICs) gene expression in single-cell RNA sequencing.

RESULTS

Collagen triple helix repeat containing 1 (CTHRC1) was identified as the key gene of CAVD. We identified a cell subtype valve interstitial cells-fibroblast, which was closely associated with fibro-calcific progress of aortic valve. CTHRC1 highly expressed in the VIC subpopulation. Immune infiltration analysis demonstrated that mast cells, B cells, dendritic cells and eosinophils were involved in pathogenesis of CAVD. Correlation analysis demonstrated that CTHRC1 was correlated with mast cells mostly.

CONCLUSIONS

In summary, the study suggested that CTHRC1 was a key gene of CAVD and CTHRC1 might participate in the potential molecular pathways involved in the connection between infiltrating immune cells and myofibroblast phenotype VICs.

Comprehensive analysis of the progression mechanisms of CRPC and its inhibitor discovery based on machine learning algorithms

Background: Almost all patients treated with androgen deprivation therapy (ADT) eventually develop castration-resistant prostate cancer (CRPC). Our research aims to elucidate the potential biomarkers and molecular mechanisms that underlie the transformation of primary prostate cancer into CRPC. Methods: We collected three microarray datasets (GSE32269, GSE74367, and GSE66187) from the Gene Expression Omnibus (GEO) database for CRPC. Differentially expressed genes (DEGs) in CRPC were identified for further analyses, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). Weighted gene coexpression network analysis (WGCNA) and two machine learning algorithms were employed to identify potential biomarkers for CRPC. The diagnostic efficiency of the selected biomarkers was evaluated based on gene expression level and receiver operating characteristic (ROC) curve analyses. We conducted virtual screening of drugs using AutoDock Vina. In vitro experiments were performed using the Cell Counting Kit-8 (CCK-8) assay to evaluate the inhibitory effects of the drugs on CRPC cell viability. Scratch and transwell invasion assays were employed to assess the effects of the drugs on the migration and invasion abilities of prostate cancer cells. Results: Overall, a total of 719 DEGs, consisting of 513 upregulated and 206 downregulated genes, were identified. The biological functional enrichment analysis indicated that DEGs were mainly enriched in pathways related to the cell cycle and metabolism. CCNA2 and CKS2 were identified as promising biomarkers using a combination of WGCNA, LASSO logistic regression, SVM-RFE, and Venn diagram analyses. These potential biomarkers were further validated and exhibited a strong predictive ability. The results of the virtual screening revealed Aprepitant and Dolutegravir as the optimal targeted drugs for CCNA2 and CKS2, respectively. In vitro experiments demonstrated that both Aprepitant and Dolutegravir exerted significant inhibitory effects on CRPC cells (p < 0.05), with Aprepitant displaying a superior inhibitory effect compared to Dolutegravir. Discussion: The expression of CCNA2 and CKS2 increases with the progression of prostate cancer, which may be one of the driving factors for the progression of prostate cancer and can serve as diagnostic biomarkers and therapeutic targets for CRPC. Additionally, Aprepitant and Dolutegravir show potential as anti-tumor drugs for CRPC.

Identification and characterization of a novel molecular classification based on disulfidptosis-related genes to predict prognosis and immunotherapy efficacy in hepatocellular carcinoma

BACKGROUND

Disulfidptosis has been discovered as a mechanism of cell death mediating by SLC7A11. Nonetheless, little is known about the relationship between disulfidptosis-related genes (DRG) and hepatocellular carcinoma (HCC).

METHODS

7 datasets including 1,302 HCC patients and 62,530 cells were downloaded. We adopted consensus clustering algorithm to construct the consensus matrix and cluster the samples' DRG related expression profile data. Then, weighted gene co-expression network analysis (WGCNA) was conducted to identify hub gene modules associated with the identified clusters and determine the correlation between modules. A DRG.score was constructed based on genes through differential analysis and WGCNA of the 2 clusters.

RESULTS

Univariate and multivariate Cox regression analysis show that SLC7A11 and LRPPRC can be used as an independent factor in HCC. Then, two molecular subgroups with significantly different survival were identified based on 10 DRG. The cluster.A showed a worse prognosis, higher immune infiltration, and higher immune checkpoint expression. Then, by differential analysis and WGCNA of the 2 clusters, we identified 5 hub genes, and constructed a DRG.score. Univariate and multivariate Cox regression analysis show that DRG.score can be used as an independent factor to predict the prognosis in HCC. Furthermore, high DRG.score group had a worse prognosis, and was validated in TCGA-LIHC, LIRI-JP, GSE14520, GSE36376, and GSE76427. Preclinically, patients with higher DRG.score demonstrated significant immunotherapy therapeutic advantages and transcatheter arterial chemoembolization clinical benefits.

CONCLUSIONS

SLC7A11 and LRPPRC play an essential role in HCC prognosis prediction. The DRG.score might become useful biomarkers for novel therapeutic targets.

Genetic characterization of primary and metastatic high-grade serous ovarian cancer tumors reveals distinct features associated with survival

High-grade serous ovarian cancer (HGSC) is the most lethal histotype of ovarian cancer and the majority of cases present with metastasis and late-stage disease. Over the last few decades, the overall survival for patients has not significantly improved, and there are limited targeted treatment options. We aimed to better characterize the distinctions between primary and metastatic tumors based on short- or long-term survival. We characterized 39 matched primary and metastatic tumors by whole exome and RNA sequencing. Of these, 23 were short-term (ST) survivors (overall survival (OS) < 3.5 years) and 16 were long-term (LT) survivors (OS > 5 years). We compared somatic mutations, copy number alterations, mutational burden, differential gene expression, immune cell infiltration, and gene fusion predictions between the primary and metastatic tumors and between ST and LT survivor cohorts. There were few differences in RNA expression between paired primary and metastatic tumors, but significant differences between the transcriptomes of LT and ST survivors in both their primary and metastatic tumors. These findings will improve the understanding of the genetic variation in HGSC that exist between patients with different prognoses and better inform treatments by identifying new targets for drug development.

Characterizing the Immune Microenvironment and Neoantigen Landscape of Hürthle Cell Carcinoma to Identify Potential Immunologic Vulnerabilities

Hürthle cell carcinoma (HCC) is a rare type of thyroid cancer with high rates of distant metastasis and recurrence. Along with the scarcity of effective systemic therapies for HCC, these factors contribute to poor clinical outcomes. The immunologic features of HCC are poorly defined and response rates with immune checkpoint blockade have not been reported. A more comprehensive understanding of the immune landscape and factors that predict response to checkpoint inhibitors is needed. We performed RNA sequencing on 40 tumors to characterize the neoantigen landscape and immune microenvironment of HCC. We analyzed transcriptomic profiles, tumor-infiltrating immune cell populations, and measures of T-cell activation/dysfunction and correlated these to genetic features such as tumor mutation burden, neoantigen burden, mitochondrial mutations, and LOH from chromosomal uniparental disomy. Finally, immune profiles of patients with recurrence were compared with those of patients without recurrence. HCC tumors exhibited low levels of immune infiltration, with the more aggressive widely invasive phenotype associated with more immune depletion. There was a negative correlation between tumor mutation burden, neoantigen burden, programmed cell death ligand 1 (PD-L1) expression, and the immune infiltration score. HCC tumors that exhibited a global LOH from chromosomal uniparental disomy or haploidization had the lowest level of immune infiltration. HCC tumors that recurred displayed an immune-depleted microenvironment associated with global LOH and aerobic glycolysis. These findings offer new insights into the functional immune landscapes and immune microenvironment of HCC. Our data identify potential immunologic vulnerabilities for these understudied and often fatal cancers.

Significance

The immune landscape of HCC is poorly defined and response rates to immunotherapy have not been reported. The authors found the immune microenvironment in HCC to be depleted. This immunosuppression is associated with a global LOH from haploidization and uniparental disomy, resulting in whole chromosome losses across the genome.

Identification and analysis of cellular senescence-associated signatures in diabetic kidney disease by integrated bioinformatics analysis and machine learning

Background

Diabetic kidney disease (DKD) is a common complication of diabetes that is clinically characterized by progressive albuminuria due to glomerular destruction. The etiology of DKD is multifactorial, and numerous studies have demonstrated that cellular senescence plays a significant role in its pathogenesis, but the specific mechanism requires further investigation.

Methods

This study utilized 5 datasets comprising 144 renal samples from the Gene Expression Omnibus (GEO) database. We obtained cellular senescence-related pathways from the Molecular Signatures Database and evaluated the activity of senescence pathways in DKD patients using the Gene Set Enrichment Analysis (GSEA) algorithm. Furthermore, we identified module genes related to cellular senescence pathways through Weighted Gene Co-Expression Network Analysis (WGCNA) algorithm and used machine learning algorithms to screen for hub genes related to senescence. Subsequently, we constructed a cellular senescence-related signature (SRS) risk score based on hub genes using the Least Absolute Shrinkage and Selection Operator (LASSO), and verified mRNA levels of hub genes by RT-PCR in vivo. Finally, we validated the relationship between the SRS risk score and kidney function, as well as their association with mitochondrial function and immune infiltration.

Results

The activity of cellular senescence-related pathways was found to be elevated among DKD patients. Based on 5 hub genes (LIMA1, ZFP36, FOS, IGFBP6, CKB), a cellular senescence-related signature (SRS) was constructed and validated as a risk factor for renal function decline in DKD patients. Notably, patients with high SRS risk scores exhibited extensive inhibition of mitochondrial pathways and upregulation of immune cell infiltration.

Conclusion

Collectively, our findings demonstrated that cellular senescence is involved in the process of DKD, providing a novel strategy for treating DKD.

Transcriptomics for Clinical and Experimental Biology Research: Hang on a Seq

Sequencing the human genome empowers translational medicine, facilitating transcriptome-wide molecular diagnosis, pathway biology, and drug repositioning. Initially, microarrays are used to study the bulk transcriptome; but now short-read RNA sequencing (RNA-seq) predominates. Positioned as a superior technology, that makes the discovery of novel transcripts routine, most RNA-seq analyses are in fact modeled on the known transcriptome. Limitations of the RNA-seq methodology have emerged, while the design of, and the analysis strategies applied to, arrays have matured. An equitable comparison between these technologies is provided, highlighting advantages that modern arrays hold over RNA-seq. Array protocols more accurately quantify constitutively expressed protein coding genes across tissue replicates, and are more reliable for studying lower expressed genes. Arrays reveal long noncoding RNAs (lncRNA) are neither sparsely nor lower expressed than protein coding genes. Heterogeneous coverage of constitutively expressed genes observed with RNA-seq, undermines the validity and reproducibility of pathway analyses. The factors driving these observations, many of which are relevant to long-read or single-cell sequencing are discussed. As proposed herein, a reappreciation of bulk transcriptomic methods is required, including wider use of the modern high-density array data-to urgently revise existing anatomical RNA reference atlases and assist with more accurate study of lncRNAs.

Comprehensive analysis identifies ARHGEF6 as a potential prognostic and immunological biomarker in lung adenocarcinoma

Lung adenocarcinoma (LUAD), the most common histological type in lung cancer, is one of leading cancers with considerable morbidity/mortality worldwide. Treating LUAD remains an outstanding challenge due to the lack of early diagnosis and the poor therapeutic effects. Rac/Cdc42 guanine nucleotide exchange factor 6 (ARHGEF6), one of cytoskeletal regulators, exerts crucial biological functions in T cell migration. The potential biological role of ARHGEF6 in LUAD has yet to be established. Using multiple bioinformatics tools and statistical methods, we discovered that the mRNA and protein expression level of ARHGEF6 was significantly downregulated in tumor tissues comparing to normal controls. Moreover, ARHGEF6 presented high diagnostic value in LUAD patients (AUC = 0.949), and the patients with low ARHGEF6 expression had more somatic mutations and poor T stage, N stage, clinical prognosis. Experimental validation indicated that ARHGEF6 was low expressed in A549 and PC-9 cells comparing to the normal lung epithelial cells. The overexpression of ARHGEF6 remarkably attenuated the abilities of cell proliferation and colony formation. Furthermore, the immune landscape analysis in TME revealed that ARHGEF6 expression was positively associated with immune cell infiltration and immune checkpoints. Single-cell transcriptome analysis indicated that ARHGEF6 expression was also distributed in immune cell types in TME based on TISCH database. Additionally, differentially expressed genes (DEGs) and functional enrichment analyses uncovered that ARHGEF6 was involved in T cell activation. Finally, LUAD samples were classified two clusters based on DEGs for subgroups analysis. In summary, this study comprehensively uncovered that ARHGEF6 could be identified as a potential prognostic and immunological biomarker in lung adenocarcinoma.

A novel risk model based on the correlation between the expression of basement membrane genes and immune infiltration to predict the invasiveness of pituitary adenomas

Objective

Invasive pituitary adenomas (IPAs) are common tumors of the nervous system tumors for which invasive growth can lead to difficult total resection and a high recurrence rate. The basement membrane (BM) is a special type of extracellular matrix and plays an important role in the invasion of pituitary adenomas (PAs). The aim of this study was to develop a risk model for predicting the invasiveness of PAs by analyzing the correlation between the expression of BM genes and immune infiltration.

Methods

Four datasets, featuring samples IPAs and non-invasive pituitary adenomas (NIPAs), were obtained from the Gene Expression Omnibus database (GEO). R software was then used to identify differentially expressed genes (DEGs) and analyze their functional enrichment. Protein-protein interaction (PPI) network was used to screen BM genes, which were analyzed for immune infiltration; this led to the generation of a risk model based on the correlation between the expression of BM genes and immunity. A calibration curve and receiver operating characteristic (ROC) curve were used to evaluate and validate the model. Subsequently, the differential expression levels of BM genes between IPA and NIPA samples collected in surgery were verified by Quantitative Polymerase Chain Reaction (qPCR) and the prediction model was further evaluated. Finally, based on our analysis, we recommend potential drug targets for the treatment of IPAs.

Results

The merged dataset identified 248 DEGs that were mainly enriching in signal transduction, the extracellular matrix and channel activity. The PPI network identified 11 BM genes from the DEGs: SPARCL1, GPC3, LAMA1, SDC4, GPC4, ADAMTS8, LAMA2, LAMC3, SMOC1, LUM and THBS2. Based on the complex correlation between these 11 genes and immune infiltration, a risk model was established to predict PAs invasiveness. Calibration curve and ROC curve analysis (area under the curve [AUC]: 0.7886194) confirmed the good predictive ability of the model. The consistency between the qPCR results and the bioinformatics results confirmed the reliability of data mining.

Conclusion

Using a variety of bioinformatics methods, we developed a novel risk model to predict the probability of PAs invasion based on the correlation between 11 BM genes and immune infiltration. These findings may facilitate closer surveillance and early diagnosis to prevent or treat IPAs in patients and improve the clinical awareness of patients at high risk of IPAs.

Normalization of Large-Scale Transcriptome Data Using Heuristic Methods

In this study, we introduce an artificial intelligent method for addressing the batch effect of a transcriptome data. The method has several clear advantages in comparison with the alternative methods presently in use. Batch effect refers to the discrepancy in gene expression data series, measured under different conditions. While the data from the same batch (measurements performed under the same conditions) are compatible, combining various batches into 1 data set is problematic because of incompatible measurements. Therefore, it is necessary to perform correction of the combined data (normalization), before performing biological analysis. There are numerous methods attempting to correct data set for batch effect. These methods rely on various assumptions regarding the distribution of the measurements. Forcing the data elements into pre-supposed distribution can severely distort biological signals, thus leading to incorrect results and conclusions. As the discrepancy between the assumptions regarding the data distribution and the actual distribution is wider, the biases introduced by such “correction methods” are greater. We introduce a heuristic method to reduce batch effect. The method does not rely on any assumptions regarding the distribution and the behavior of data elements. Hence, it does not introduce any new biases in the process of correcting the batch effect. It strictly maintains the integrity of measurements within the original batches.

Transcriptome Analysis Reveals Modulation of Human Stem Cells from the Apical Papilla by Species Associated with Dental Root Canal Infection

Interaction of oral bacteria with stem cells from the apical papilla (SCAP) can negatively affect the success of regenerative endodontic treatment (RET). Through RNA-seq transcriptomic analysis, we studied the effect of the oral bacteria Fusobacterium nucleatum and Enterococcus faecalis, as well as their supernatants enriched by bacterial metabolites, on the osteo- and dentinogenic potential of SCAPs in vitro. We performed bulk RNA-seq, on the basis of which differential expression analysis (DEG) and gene ontology enrichment analysis (GO) were performed. DEG analysis showed that E. faecalis supernatant had the greatest effect on SCAPs, whereas F. nucleatum supernatant had the least effect (Tanimoto coefficient = 0.05). GO term enrichment analysis indicated that F. nucleatum upregulates the immune and inflammatory response of SCAPs, and E. faecalis suppresses cell proliferation and cell division processes. SCAP transcriptome profiles showed that under the influence of E. faecalis the upregulation of VEGFA, Runx2, and TBX3 genes occurred, which may negatively affect the SCAP's osteo- and odontogenic differentiation. F. nucleatum downregulates the expression of WDR5 and TBX2 and upregulates the expression of TBX3 and NFIL3 in SCAPs, the upregulation of which may be detrimental for SCAPs' differentiation potential. In conclusion, the present study shows that in vitro, F. nucleatum, E. faecalis, and their metabolites are capable of up- or downregulating the expression of genes that are necessary for dentinogenic and osteogenic processes to varying degrees, which eventually may result in unsuccessful RET outcomes. Transposition to the clinical context merits some reservations, which should be approached with caution.

Integrated bioinformatics-based identification of diagnostic markers in Alzheimer disease

Alzheimer disease (AD) is a progressive neurodegenerative disease resulting from the accumulation of extracellular amyloid beta (Aβ) and intracellular neurofibrillary tangles. There are currently no objective diagnostic measures for AD. The aim of this study was to identify potential diagnostic markers for AD and evaluate the role of immune cell infiltration in disease pathogenesis. AD expression profiling data for human hippocampus tissue (GSE48350 and GSE5281) were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified using R software and the Human Protein Atlas database was used to screen AD-related DEGs. We performed functional enrichment analysis and established a protein-protein interaction (PPI) network to identify disease-related hub DEGs. The fraction of infiltrating immune cells in samples was determined with the Microenvironment Cell Populations-counter method. The random forest algorithm was used to develop a prediction model and receiver operating characteristic (ROC) curve analysis was performed to validate the diagnostic utility of the candidate AD markers. The correlation between expression of the diagnostic markers and immune cell infiltration was also analyzed. A total of 107 AD-related DEGs were screened in this study, including 28 that were upregulated and 79 that were downregulated. The DEGs were enriched in the Gene Ontology terms GABAergic synapse, Morphine addiction, Nicotine addiction, Phagosome, and Synaptic vesicle cycle. We identified 10 disease-related hub genes and 20 candidate diagnostic genes. Synaptophysin (SYP) and regulator of G protein signaling 4 (RGS4) (area under the ROC curve = 0.909) were verified as potential diagnostic markers for AD in the GSE28146 validation dataset. Natural killer cells, B lineage cells, monocytic lineage cells, endothelial cells, and fibroblasts were found to be involved in AD; additionally, the expression levels of both SYP and RGS4 were negatively correlated with the infiltration of these immune cell types. These results suggest that SYP and RGS4 are potential diagnostic markers for AD and that immune cell infiltration plays an important role in AD development and progression.

Identification of key genes and immune cell infiltration in recurrent implantation failure: A study based on integrated analysis of multiple microarray studies

PROBLEM

Recurrent implantation failure (RIF) refers to a challenging topic in assisted reproductive technology (ART), the etiology of which may be attributed to impaired endometrial receptivity; however, the precise pathogenesis of RIF has not been thoroughly elucidated.

METHOD OF STUDY

Four RIF microarray datasets were obtained from the Gene Expression Omnibus database and integrated by the "sva" R package. The differentially expressed genes (DEGs) were analyzed using the "limma" package and then GO, KEGG, GSEA, and GSVA were applied to perform functional and pathway enrichment analysis. The immune cell infiltration in the RIF process was evaluated by the CIBERSORT algorithm. Finally, the hub genes were identified through the CytoHubba and subsequently verified using two items of external endometrial data.

RESULTS

236 genes were differentially expressed in the endometrium of the RIF group. Functional enrichment analysis demonstrated that the biological functions of DEGs were mainly correlated to the immune-related pathways, including immune response, TNF signaling pathway, complement and coagulation cascades. Among the immune cells, γδ T cells decreased significantly in the endometrium of RIF patients. In addition, the key DEGs such as PTGS2, FGB, MUC1, SST, VCAM1, MMP7, ERBB4, FOLR1, and C3 were screened and identified as the hub genes involved in the pathogenesis of RIF.

CONCLUSIONS

Abnormal immune response regulation of endometrium contributes to the occurrence of RIF, and γδ T cells may be the pivotal immune cells causing RIF. At the same time, the novel hub genes identified will provide effective targets for the prediction and therapy of RIF.

Identification and validation of P4HB as a novel autophagy-related biomarker in diabetic nephropathy

Diabetic nephropathy (DN), a frequent microvascular complication of diabetes, has been recognized as a primary cause of chronic kidney disease (CKD) and end-stage renal disease (ESRD). Previous studies found that autophagy of renal tubular epithelial cells plays an important role in DN pathogenesis. Our research aimed to investigate the differentially expressed autophagy-related genes (DEARGs) between DN and healthy renal tubule samples and identify a novel autophagy-related biomarker associated with tubulointerstitial injury in DN. In this study, gene expression profiles of renal tubules from 10 DN patients and 24 healthy controls in the GSE30122 dataset were analyzed, and 43 DEARGs were identified by bioinformatics analysis. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and correlation analysis were performed on DEARGs, and the hub gene prolyl 4-hydroxylase subunit beta (P4HB) was screened by protein–protein interaction and verified by utilizing other datasets and stimulating HK-2 cells under high glucose concentration. We found that the expression of P4HB in renal tubules was correlated with renal function. In summary, our research provided novel insights for comprehension of DN molecular mechanisms and identified P4HB as a novel autophagy-related biomarker of DN.

ATF3 as a potential diagnostic marker of early-stage osteoarthritis and its correlation with immune infiltration through bioinformatics analysis

AIMS

This study aimed, through bioinformatics analysis, to identify the potential diagnostic markers of osteoarthritis, and analyze the role of immune infiltration in synovial tissue.

METHODS

The gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified by R software. Functional enrichment analyses were performed and protein-protein interaction networks (PPI) were constructed. Then the hub genes were screened. Biomarkers with high value for the diagnosis of early osteoarthritis (OA) were validated by GEO datasets. Finally, the CIBERSORT algorithm was used to evaluate the immune infiltration between early-stage OA and end-stage OA, and the correlation between the diagnostic marker and infiltrating immune cells was analyzed.

RESULTS

A total of 88 DEGs were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that DEGs were significantly enriched in leucocyte migration and interleukin (IL)-17 signalling pathways. Disease ontology (DO) indicated that DEGs were mostly enriched in rheumatoid arthritis. Six hub genes including FosB proto-oncogene, AP-1 transcription factor subunit (FOSB); C-X-C motif chemokine ligand 2 (CXCL2); CXCL8; IL-6; Jun proto-oncogene, AP-1 transcription factor subunit (JUN); and Activating transcription factor 3 (ATF3) were identified and verified by GEO datasets. ATF3 (area under the curve = 0.975) turned out to be a potential biomarker for the diagnosis of early OA. Several infiltrating immune cells varied significantly between early-stage OA and end-stage OA, such as resting NK cells (p = 0.016), resting dendritic cells (p = 0.043), and plasma cells (p = 0.043). Additionally, ATF3 was significantly correlated with resting NK cells (p = 0.034), resting dendritic cells (p = 0.026), and regulatory T cells (Tregs, p = 0.018).

CONCLUSION

ATF3 may be a potential diagnostic marker for early diagnosis and treatment of OA, and immune cell infiltration provides new perspectives for understanding the mechanism during OA progression.Cite this article: Bone Joint Res 2022;11(9):679-689.

Identification of potential biomarkers of inflammation-related genes for ischemic cardiomyopathy

Objective

Inflammation plays an important role in the pathophysiology of ischemic cardiomyopathy (ICM). We aimed to identify potential biomarkers of inflammation-related genes for ICM and build a model based on the potential biomarkers for the diagnosis of ICM.

Materials and methods

The microarray datasets and RNA-Sequencing datasets of human ICM were downloaded from the Gene Expression Omnibus database. We integrated 8 microarray datasets via the SVA package to screen the differentially expressed genes (DEGs) between ICM and non-failing control samples, then the differentially expressed inflammation-related genes (DEIRGs) were identified. The least absolute shrinkage and selection operator, support vector machine recursive feature elimination, and random forest were utilized to screen the potential diagnostic biomarkers from the DEIRGs. The potential biomarkers were validated in the RNA-Sequencing datasets and the functional experiment of the ICM rat, respectively. A nomogram was established based on the potential biomarkers and evaluated via the area under the receiver operating characteristic curve (AUC), calibration curve, decision curve analysis (DCA), and Clinical impact curve (CIC).

Results

64 DEGs and 19 DEIRGs were identified, respectively. 5 potential biomarkers (SERPINA3, FCN3, PTN, CD163, and SCUBE2) were ultimately selected. The validation results showed that each of these five potential biomarkers showed good discriminant power for ICM, and their expression trends were consistent with the bioinformatics results. The results of AUC, calibration curve, DCA, and CIC showed that the nomogram demonstrated good performance, calibration, and clinical utility.

Conclusion

SERPINA3, FCN3, PTN, CD163, and SCUBE2 were identified as potential biomarkers associated with the inflammatory response to ICM. The proposed nomogram could potentially provide clinicians with a helpful tool to the diagnosis and treatment of ICM from an inflammatory perspective.

Inhibition of IL1R1 or CASP4 attenuates spinal cord injury through ameliorating NLRP3 inflammasome-induced pyroptosis

Spinal cord injury (SCI) is a devastating trauma characterized by serious neuroinflammation and permanent neurological dysfunction. However, the molecular mechanism of SCI remains unclear, and few effective medical therapies are available at present. In this study, multiple bioinformatics methods were used to screen out novel targets for SCI, and the mechanism of these candidates during the progression of neuroinflammation as well as the therapeutic effects were both verified in a rat model of traumatic SCI. As a result, CASP4, IGSF6 and IL1R1 were identified as the potential diagnostic and therapeutic targets for SCI by computational analysis, which were enriched in NF-κB and IL6-JAK-STATA3 signaling pathways. In the injured spinal cord, these three signatures were up-regulated and closely correlated with NLRP3 inflammasome formation and gasdermin D (GSDMD) -induced pyroptosis. Intrathecal injection of inhibitors of IL1R1 or CASP4 improved the functional recovery of SCI rats and decreased the expression of these targets and inflammasome component proteins, such as NLRP3 and GSDMD. This treatment also inhibited the pp65 activation into the nucleus and apoptosis progression. In conclusion, our findings of the three targets shed new light on the pathogenesis of SCI, and the use of immunosuppressive agents targeting these proteins exerted anti-inflammatory effects against spinal cord inflammation by inhibiting NF-kB and NLRP3 inflammasome activation, thus blocking GSDMD -induced pyroptosis and immune activation.

Systematic identification of key extracellular proteins as the potential biomarkers in lupus nephritis

Background

Lupus nephritis (LN) is the most common and severe clinical manifestation of systemic lupus erythematosus (SLE) with considerable morbidity/mortality and limited treatment options. Since kidney biopsy is a relative hysteretic indicator, it is indispensable to investigate potential biomarkers for early diagnosis and predicting clinical outcomes of LN patients. Extracellular proteins may become the promising biomarkers by the secretion into body fluid. Our study linked extracellular proteins with lupus nephritis to identify the emerging biomarkers.

Methods

The expression profiling data were acquired from the Gene Expression Omnibus (GEO) database. Meanwhile, the two gene lists encoding extracellular proteins were collected from the Human Protein Atlas (HPA) and UniProt database. Subsequently, the extracellular protein-differentially expressed genes (EP-DEGs) were screened out, and the key EP-DEGs were determined by MCODE, MCC, and Degree methods via the protein–protein interaction (PPI) network. The expression level, immune characteristics, and diagnostic value of these candidate biomarkers were investigated. Finally, the Nephroseq V5 tool was applied to evaluate the clinical significance of the key EP-DEGs.

Results

A total of 164 DEGs were acquired by comparing LN samples with healthy controls based on GSE32591 datasets. Then, 38 EP-DEGs were screened out through the intersection between DEGs and extracellular protein gene lists. Function enrichment analysis indicated that these EP-DEGs might participate in immune response and constitute the extracellular matrix. Four key EP-DEGs (LUM, TGFBI, COL1A2, and POSTN) were eventually identified as candidate biomarkers, and they were all overexpressed in LN samples. Except that LUM expression was negatively correlated with most of the immune regulatory genes, there was a positive correlation between the remaining three biomarkers and the immune regulatory genes. In addition, these biomarkers had high diagnostic value, especially the AUC value of the LUM–TGFBI combination which reached almost 1 (AUC = 0.973), demonstrating high accuracy in distinguishing LN from controls. Finally, we found a meaningful correlation of these biomarkers with sex, WHO class, and renal function such as glomerular filtration rate (GFR), serum creatinine level, and proteinuria.

Conclusion

In summary, our study comprehensively identified four key EP-DEGs exerting a vital role in LN diagnosis and pathogenesis and serving as promising therapeutic targets.

Necroptosis-Related Genes Signatures Identified Molecular Subtypes and Underlying Mechanisms in Hepatocellular Carcinoma

Background

Although emerging evidence supports the relationship between necroptosis (NEC) related genes and hepatocellular carcinoma (HCC), the contribution of these necroptosis-related genes to the development, prognosis, and immunotherapy of HCC is unclear.

Methods

The expression of genes and relevant clinical information were downloaded from TCGA-LIHC, LIRI-JP, GSE14520/NCI, GSE36376, GSE76427, GSE20140, GSE27150, and IMvigor210 datasets. Next, we used an unsupervised clustering method to assign the samples into phenotype clusters base on 15 necroptosis-related genes. Subsequently, we constructed a NEC score based on NEC phenotype-related prognostic genes to quantify the necroptosis related subtypes of individual patients.

Results

We divided the samples into the high and low NEC score groups, and the high NEC score showed a poor prognosis. Simultaneously, NEC score is an effective and stable model and had a good performance in predicting the prognosis of HCC patients. A high NEC score was characterized by activation of the stroma and increased levels of immune infiltration. A high NEC score was also related to low expression of immune checkpoint molecules (PD-1/PD-L1). Importantly, the established NEC score would contribute to predicting the response to anti-PD-1/L1 immunotherapy.

Conclusions

Our study provide a comprehensive analysis of necroptosis-related genes in HCC. Stratification based on the NEC score may enable HCC patients to benefit more from immunotherapy and help identify new cancer treatment strategies.

Integration of Genetic and Immune Infiltration Insights into Data Mining of Multiple Sclerosis Pathogenesis

Background

Multiple sclerosis (MS) is an immune-mediated demyelinating disease of the central nervous system. MS pathogenesis is closely related to the environment, genetic, and immune system, but the underlying interactions have not been clearly elucidated. This study aims to unveil the genetic basis and immune landscape of MS pathogenesis with bioinformatics.

Methods

Gene matrix was retrieved from the gene expression database NCBI-GEO. Then, bioinformatics was used to standardize the samples and obtain differentially expressed genes (DEGs). The protein-protein interaction network was constructed with DEGs on the STRING website. Cytohubba plug-in and MCODE plug-in were used to mine hub genes. Meanwhile, the CIBERSORTX algorithm was used to explore the characteristics of immune cell infiltration in MS brain tissues. Spearman correlation analysis was performed between genes and immune cells, and the correlation between genes and different types of brain tissues was also analyzed using the WGCNA method.

Results

A total of 90 samples from 2 datasets were included, and 882 DEGs and 10 hub genes closely related to MS were extracted. Functional enrichment analysis suggested the role of immune response in MS. Besides, CIBERSORTX algorithm results showed that MS brain tissues contained a variety of infiltrating immune cells. Correlation analysis suggested that the hub genes were highly relevant to chronic active white matter lesions. Certain hub genes played a role in the activation of immune cells such as macrophages and natural killer cells.

Conclusions

Our study shall provide guidance for the further study of the genetic basis and immune infiltration mechanism of MS.

Transcriptomic analysis reveals pathophysiological relationship between chronic obstructive pulmonary disease (COPD) and periodontitis

Background

The aim of this study was to detect potential crosstalk genes, pathways and immune cells between periodontitis and chronic obstructive pulmonary disease (COPD).

Methods

Chronic periodontitis (CP, GSE156993) and COPD (GSE42057, GSE94916) datasets were downloaded. Differential expressed genes (DEGs; p < 0.05) were assessed and screened for overlapping results, following functional pathway enrichment analyses (p < 0.05). The xCell method was used to assess immune cell infiltration relationship between CP and COPD. Features of the detected cross-talk genes were revealed using conventional Recursive Feature Elimination (RFE) algorithm in R project. Receiver-operating characteristic curves were applied to evaluate the predictive value of the genes. Furthermore, Pearson correlation analysis was performed on crosstalk markers and infiltrating immune cells in CP and COPD, respectively.

Results

A total of 904 DEGs of COPD and 763 DEGs of CP were acquired, showing 22 overlapping DEGs between the two diseases. Thereby 825 nodes and 923 edges were found in the related protein–protein-interaction network. Eight immune cell pairs were found to be highly correlated to both CP and COPD (|correlation coefficients |> 0.5 and p-value < 0.05). Most immune cells were differently expressed between COPD and CP. RFE identified three crosstalk genes, i.e. EPB41L4A-AS1, INSR and R3HDM1. In correlation analysis, INSR was positively correlated with Hepatocytes in CP (r = 0.6714, p = 0.01679) and COPD (r = 0.5209, p < 0.001). R3HDM was positively correlated with Th1 cells in CP (r = 0.6783, p = 0.0153) and COPD (r = 0.4120, p < 0.01).

Conclusion

EPB41L4A-AS1, INSR and R3HDM1 are potential crosstalk genes between COPD and periodontitis. R3HDM was positively correlated with Th1 cells in both diseases, while INSR was positively correlated with Hepatocytes in periodontitis and COPD, supporting a potential pathophysiological relationship between periodontitis and COPD.

Characterization of RNA modifications in gastric cancer to identify prognosis-relevant gene signatures

BACKGROUND

Most human genes have diverse transcript isoforms, which mainly arise from alternative cleavage and polyadenylation (APA) at 3' ends. N7-methylguanosine (m⁷ G) is also an essential epigenetic modification at the 5' end. However, the contribution of these two RNA modifications to the development, prognosis, regulation mechanisms, and drug sensitivity of gastric cancer (GC) is unclear.

METHODS

The expression data of 2412 patients were extracted from 12 cohorts and the RNA modification patterns of 20 marker genes were systematically identified into phenotypic clusters using the unsupervised clustering approach. Following that, we developed an RNA modification model (RMscore) to quantify each GC patient's RNA modification index. Finally, we examined the correlation between RMscore and clinical features such as survival outcomes, molecular subtypes identified by the Asian Cancer Research Group (ACRG), posttranscriptional regulation, and chemotherapeutic sensitivity in GC.

RESULTS

The samples were categorized into two groups on the basis of their RMscore: high and low. The group with a low RMscore had a bad prognosis. Moreover, the low RMscore was associated with KRAS, Hedgehog, EMT, and TGF-β signaling, whereas a high RMscore was related to abnormal cell cycle signaling pathway activation. The findings also revealed that the RMscore contributes to the regulation of the miRNA-mRNA network. Drug sensitivity analysis revealed that RMscore is associated with the response to some anticancer drugs.

CONCLUSIONS

The RMscore model has the potential to be a useful tool for prognosis prediction in patients with GC. A comprehensive investigation of APA-RNA and m⁷ G-RNA modifications may reveal novel insights into the epigenetics of GC and aid in the development of more effective treatment strategies.

Development and Validation of a Novel Gene Signature for Predicting the Prognosis of Idiopathic Pulmonary Fibrosis Based on Three Epithelial-Mesenchymal Transition and Immune-Related Genes

Background: Increasing evidence has revealed that epithelial–mesenchymal transition (EMT) and immunity play key roles in idiopathic pulmonary fibrosis (IPF). However, correlation between EMT and immune response and the prognostic significance of EMT in IPF remains unclear.

Methods: Two microarray expression profiling datasets (GSE70866 and GSE28221) were downloaded from the Gene Expression Omnibus (GEO) database. EMT- and immune-related genes were identified by gene set variation analysis (GSVA) and the Estimation of STromal and Immune cells in MAlignant Tumors using Expression data (ESTIMATE) algorithm. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to investigate the functions of these EMT- and immune-related genes. Cox and least absolute shrinkage and selection operator (LASSO) regression analyses were used to screen prognostic genes and establish a gene signature. Gene Set Enrichment Analysis (GSEA) and Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) were used to investigate the function of the EMT- and immune-related signatures and correlation between the EMT- and immune-related signatures and immune cell infiltration. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed to investigate the mRNA expression of genes in the EMT- and immune-related signatures.

Results: Functional enrichment analysis suggested that these genes were mainly involved in immune response. Moreover, the EMT- and immune-related signatures were constructed based on three EMT- and immune-related genes (IL1R2, S100A12, and CCL8), and the K–M and ROC curves presented that the signature could affect the prognosis of IPF patients and could predict the 1-, 2-, and 3-year survival well. Furthermore, a nomogram was developed based on the expression of IL1R2, S100A12, and CCL8, and the calibration curve showed that the nomogram could visually and accurately predict the 1-, 2-, 3-year survival of IPF patients. Finally, we further found that immune-related pathways were activated in the high-risk group of patients, and the EMT- and immune-related signatures were associated with NK cells activated, macrophages M0, dendritic cells resting, mast cells resting, and mast cells activated. qRT-PCR suggested that the mRNA expression of IL1R2, S100A12, and CCL8 was upregulated in whole blood of IPF patients compared with normal samples.

Conclusion: IL1R2, S100A12, and CCL8 might play key roles in IPF by regulating immune response and could be used as prognostic biomarkers of IPF.

Identification and Verification of Feature Biomarkers Associated With Immune Cells in Dilated Cardiomyopathy by Bioinformatics Analysis

Objective: To explore immune-related feature genes in patients with dilated cardiomyopathy (DCM).

Methods: Expression profiles from three datasets (GSE1145, GSE21610 and GSE21819) of human cardiac tissues of DCM and healthy controls were downloaded from the GEO database. After data preprocessing, differentially expressed genes (DEGs) were identified by the ‘limma’ package in R software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were then performed to identify biological functions of the DEGs. The compositional patterns of stromal and immune cells were estimated using xCell. Hub genes and functional modules were identified based on protein-protein interaction (PPI) network analysis by STRING webtool and Cytoscape application. Correlation analysis was performed between immune cell subtypes and hub genes. Hub genes with |correlation coefficient| > 0.5 and p value <0.05 were selected as feature biomarkers. A logistic regression model was constructed based on the selected biomarkers and validated in datasets GSE5406 and GSE57338.

Results: A total of 1,005 DEGs were identified. Functional enrichment analyses indicated that extracellular matrix remodeling and immune and inflammation disorder played important roles in the pathogenesis of DCM. Immune cells, including CD8⁺ T-cells, macrophages M1 and Th1 cells, were proved to be significantly changed in DCM patients by immune cell infiltration analysis. In the PPI network analysis, STAT3, IL6, CCL2, PIK3R1, ESR1, CCL5, IL17A, TLR2, BUB1B and MYC were identified as hub genes, among which CCL2, CCL5 and TLR2 were further screened as feature biomarkers by using hub genes and immune cells correlation analysis. A diagnosis model was successfully constructed by using the three biomarkers with area under the curve (AUC) scores 0.981, 0.867 and 0.946 in merged dataset, GSE5406 and GSE57338, respectively.

Conclusion: The present study identified three immune-related genes as diagnostic biomarkers for DCM, providing a novel perspective of immune and inflammatory response for the exploration of DCM molecular mechanisms.

Identification of Immune-Related Genes Concurrently Involved in Critical Illnesses Across Different Etiologies: A Data-Driven Analysis

Severe trauma and sepsis can lead to multiple organ dysfunction syndrome, which is a leading cause of death in intensive care units with mortality rates in excess of 50%. In addition to infection, the degree of immuno-inflammatory response also influences the outcome. The genomic changes observed after a variety of pathophysiological insults, such as trauma, sepsis, burns are similar, and consist of innate immune activation and adaptive immunity suppression. However, the characteristics of the shared mechanisms of aforementioned critical illnesses and the clinical relevance remain less explored. In the present study, we performed a data analysis to identify functional genes concurrently involved in critical illnesses across differing etiologies (trauma and sepsis derived from community-acquired pneumonia/abdominal source) and explored the shared signaling pathways these common genes involved in to gain insight into the underlying molecular mechanisms. A number of immune-related biological functions were found to be dysregulated in both trauma and sepsis in the present study, so we continued to identify immune-related common genes, profiled the immune cell proportion, and explored the relationships between them. The diagnostic and prognostic value of the immune-related common genes was also evaluated to address their potential clinical utilization as novel biomarkers. Notably, we identified a list of 14 immune-related genes concurrently dysregulated in trauma and sepsis showing favorable diagnostic value, among which S100P can predict prognosis of sepsis patients. Moreover, a spectrum of immune cell subsets including naïve B cells, CD8+ T cells, CD4+ memory resting T cells, activated NK cells, resting dendritic cells, plasma cells, Tregs, macrophages M0 and macrophages M1 was found to be concurrently dysregulated in both trauma and sepsis, and a close relation between above identified immune-related genes and immune cell subsets was observed. Our data-driven findings lay a foundation for future research to elucidate the pathophysiology regarding the aspect of inflammatory and immune response in critical illnesses, and suggest future studies focus on interpreting the function roles of the identified immune-related genes, as well as the reactive immune cell subsets.

Expression of CCL2, FOS, and JUN May Help to Distinguish Patients With IgA Nephropathy From Healthy Controls

Background

IgA nephropathy (IgAN), the most common type of glomerulonephritis worldwide, can only be diagnosed mainly by renal biopsy owing to lack of effective biomarkers. It is urgent to explore and identify the potential diagnostic biomarkers through assessing the gene expression profiles of patients with IgAN.

Methods

Two datasets were obtained from the Gene Expression Omnibus (GEO) database, including GSE115857 (55 IgAN, 7 living healthy donors) and GSE35487 (25 IgAN, 6 living healthy donors), then underwent differentially expressed genes (DEGs) and function enrichment analyses utilizing R packages. The common gene list was screened out between DEGs and immune-associated genes by Venn diagram, then performed gene-gene interaction, protein-protein interaction (PPI) and function enrichment analyses. Top three immune-associated hub genes were selected by Maximal Clique Centrality (MCC) method, then the expression and diagnostic value of these hub genes were determined. Consensus clustering algorithm was applied to conduct the unsupervised cluster analysis of the immune-associated hub gene list in IgAN. Finally, the Nephroseq V5 tool was applied to identify the expression level of CCL2, FOS, JUN in kidney diseases, as well as the correlation between CCL2, FOS, JUN expression and renal function in the patients with IgAN.

Results

A total of 129 DEGs were obtained through comparing IgAN with healthy controls via the GSE115857 and GSE35487 datasets. Then, we screened out 24 immune-associated IgAN DEGs. CCL2, JUN, and FOS were identified as the top three hub genes, and they were all remarkably downregulated in IgAN. More importantly, CCL2, JUN, and FOS had a high accuracy [area under the curve (AUC) reached almost 1] in predicting IgAN, which could easily distinguish between IgAN patients and healthy individuals. Three distinct subgroups of IgAN were determined based on 24 immune-associated DEGs, with significant differences in the expression of CCL2, JUN, and FOS genes. Finally, CCL2, FOS, JUN were manifested a meaningful association with proteinuria, glomerular filtration rate (GFR), and serum creatinine level.

Conclusion

In summary, our study comprehensively uncovers that CCL2, JUN, and FOS may function as promising biomarkers for diagnosis of IgAN.

Integrative Analyses of Genes Associated With Otologic Disorders in Turner Syndrome

Background: Loss or partial loss of one X chromosome induces Turner syndrome (TS) in females, causing major medical concerns, including otologic disorders. However, the underlying genetic pathophysiology of otologic disorders in TS is mostly unclear.

Methods: Ear-related genes of TS (TSEs) were identified by analyzing differentially expressed genes (DEGs) in two Gene Expression Omnibus (GEO)-derived expression profiles and ear-genes in the Comparative Toxicogenomic Database (CTD). Subsequently, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Disease Ontology (DO) analyses; Gene Set Enrichment Analysis (GSEA); and Gene Set Variation Analysis (GSVA) were adopted to study biological functions. Moreover, hub genes within the TSEs were identified by assessing protein-protein interaction (PPI), gene-microRNA, and gene-transcription factor (TF) networks. Drug-Gene Interaction Database (DGIdb) analysis was performed to predict molecular drugs for TS. Furthermore, three machine-learning analysis outcomes were comprehensively compared to explore optimal biomarkers of otologic disorders in TS. Finally, immune cell infiltration was analyzed.

Results: The TSEs included 30 significantly upregulated genes and 14 significantly downregulated genes. Enrichment analyses suggested that TSEs play crucial roles in inflammatory responses, phospholipid and glycerolipid metabolism, transcriptional processes, and epigenetic processes, such as histone acetylation, and their importance for inner ear development. Subsequently, we described three hub genes in the PPI network and confirmed their involvement in Wnt/β-catenin signaling pathway and immune cell regulation and roles in maintaining normal auditory function. We also constructed gene-microRNA and gene-TF networks. A novel biomarker (SLC25A6) of the pathogenesis of otologic disorders in TS was identified by comprehensive comparisons of three machine-learning analyses with the best predictive performance. Potential therapeutic agents in TS were predicted using the DGIdb. Immune cell infiltration analysis showed that TSEs are related to immune-infiltrating cells.

Conclusion: Overall, our findings have deepened the understanding of the pathophysiology of otologic disorders in TS and made contributions to present a promising biomarker and treatment targets for in-depth research.

Expression Patterns of Glycosylation Regulators Define Tumor Microenvironment and Immunotherapy in Gastric Cancer

Glycosylation (Glyc) is prevalently related to gastric cancer (GC) pathophysiology. However, studies on the relationship between glycosylation regulators and tumor microenvironment (TME) and immunotherapy of GC remain scarce. We extracted expression data of 1,956 patients with GC from eight cohorts and systematically characterized the glycosylation patterns of six marker genes into phenotype clusters using the unsupervised clustering method. Next, we constructed a Glyc. score to quantify the glycosylation index of each patient with GC. Finally, we analyzed the relationship between Glyc. score and clinical traits including molecular subtype, TME, and immunotherapy of GC. On the basis of prognostic glycosylation-related differentially expressed genes, we constructed the Glyc. score and divided the samples into the high– and low–Glyc. score groups. The high–Glyc. score group showed a poor prognosis and was validated in multiple cohorts. Functional enrichment analysis revealed that the high–Glyc. score group was enriched in metabolism-related pathways. Furthermore, the high–Glyc. score group was associated with the infiltration of immune cells. Importantly, the established Glyc. score would contribute to predicting the response to anti–PD-1/L1 immunotherapy. In conclusion, the Glyc. score is a potentially useful tool to predict the prognosis of GC. Comprehensive analysis of glycosylation may provide novel insights into the epigenetics of GC and improve treatment strategies.

Machine learning approaches for prediction of bipolar disorder based on biological, clinical and neuropsychological markers: a systematic review and meta-analysis

Applying machine learning (ML) to objective markers may overcome prognosis uncertainty due to the subjective nature of the diagnosis of bipolar disorder (BD). This PRISMA-compliant meta-analysis provides new systematic evidence of the BD classification accuracy reached by different markers and ML algorithms. We focused on neuroimaging, electrophysiological techniques, peripheral biomarkers, genetic data, neuropsychological or clinical measures, and multimodal approaches. PubMed, Embase and Scopus were searched through 3^rd December 2020. Meta-analyses were performed using random-effect models. Overall, 81 studies were included in this systematic review and 65 in the meta-analysis (11,336 participants, 3,903 BD). The overall pooled classification accuracy was 0.77 (95%CI[0.75;0.80]). Despite subgroup analyses for diagnostic comparison group, psychiatric disorders, marker, ML algorithm, and validation procedure were not significant, linear discriminant analysis significantly outperformed support vector machine for peripheral biomarkers (p=0.03). Sample size was inversely related to accuracy. Evidence of publication bias was detected. Ultimately, although ML reached a high accuracy in differentiating BD from other psychiatric disorders, best practices in methodology are needed for the advancement of future studies.

Identification and Verification of Diagnostic Biomarkers for Glomerular Injury in Diabetic Nephropathy Based on Machine Learning Algorithms

Diabetic nephropathy (DN) is regarded as the leading cause of end-stage renal disease worldwide and lacks novel therapeutic targets. To screen and verify special biomarkers for glomerular injury in patients with DN, fifteen datasets were retrieved from the Gene Expression Omnibus (GEO) database, correspondingly divided into training and testing cohorts and then merged. Using the limma package, 140 differentially expressed genes (DEGs) were screened out between 81 glomerular DN samples and 41 normal ones from the training cohort. With the help of the ConsensusClusterPlus and WGCNA packages, the 81 glomerular DN samples were distinctly divided into two subclusters, and two highly associated modules were identified. By using machine learning algorithms (LASSO, RF, and SVM-RFE) and the Venn diagram, two overlapping genes (PRKAR2B and TGFBI) were finally determined as potential biomarkers, which were further validated in external testing datasets and the HFD/STZ-induced mouse models. Based on the biomarkers, the diagnostic model was developed with reliable predictive ability for diabetic glomerular injury. Enrichment analyses indicated the apparent abnormal immune status in patients with DN, and the two biomarkers played an important role in the immune microenvironment. The identified biomarkers demonstrated a meaningful correlation between the immune cells' infiltration and renal function. In conclusion, two robust genes were identified as diagnostic biomarkers and may serve as potential targets for therapeutics of DN, which were closely associated with multiple immune cells.

An Robust Rank Aggregation and Least Absolute Shrinkage and Selection Operator Analysis of Novel Gene Signatures in Dilated Cardiomyopathy

Objective: Dilated cardiomyopathy (DCM) is a heart disease with high mortality characterized by progressive cardiac dilation and myocardial contractility reduction. The molecular signature of dilated cardiomyopathy remains to be defined. Hence, seeking potential biomarkers and therapeutic of DCM is urgent and necessary.

Methods: In this study, we utilized the Robust Rank Aggregation (RRA) method to integrate four eligible DCM microarray datasets from the GEO and identified a set of significant differentially expressed genes (DEGs) between dilated cardiomyopathy and non-heart failure. Moreover, LASSO analysis was carried out to clarify the diagnostic and DCM clinical features of these genes and identify dilated cardiomyopathy derived diagnostic signatures (DCMDDS).

Results: A total of 117 DEGs were identified across the four microarrays. Furthermore, GO analysis demonstrated that these DEGs were mainly enriched in the regulation of inflammatory response, the humoral immune response, the regulation of blood pressure and collagen–containing extracellular matrix. In addition, KEGG analysis revealed that DEGs were mainly enriched in diverse infected signaling pathways. Moreover, Gene set enrichment analysis revealed that immune and inflammatory biological processes such as adaptive immune response, cellular response to interferon and cardiac muscle contraction, dilated cardiomyopathy are significantly enriched in DCM. Moreover, Least absolute shrinkage and selection operator (LASSO) analyses of the 18 DCM-related genes developed a 7-gene signature predictive of DCM. This signature included ANKRD1, COL1A1, MYH6, PERELP, PRKACA, CDKN1A, and OMD. Interestingly, five of these seven genes have a correlation with left ventricular ejection fraction (LVEF) in DCM patients.

Conclusion: Our present study demonstrated that the signatures could be robust tools for predicting DCM in clinical practice. And may also be potential treatment targets for clinical implication in the future.

Bioinformatics analysis identifies immune-related gene signatures and subtypes in diabetic nephropathy

BACKGROUND

Systemic inflammation and immune response are involved in the pathogenesis of diabetic nephropathy (DN). However, the specific immune-associated signature during DN development is unclear. Our study aimed to reveal the roles of immune-related genes during DN progression.

METHODS

The GSE30529 and GSE30528 datasets were acquired from the Gene Expression Omnibus (GEO) database. Then, the intersection between differentially expressed genes (DEGs) and immune score-related genes (ISRGs) was screened. Subsequently, functional enrichment analyses were performed. The different immune phenotype-related subgroups were finally divided using unsupervised clustering. The core genes were identified by WGCNA and the protein-protein interaction (PPI) network. xCell algorithm was applied to assess the proportion of immune cell infiltration.

RESULTS

92 immune score-related DEGs (ISRDEGs) were identified, and these genes were enriched in inflammation- and immune-associated pathways. Furthermore, two distinct immune-associated subgroups (C1 and C2) were identified, and the C1 subgroup exhibited activated immune pathways and a higher percentage of immune cells compared to the C2 subgroup. Two core genes (LCK and HCK) were identified and all up-regulated in DN, and the expressions were verified using GSE30122, GSE142025, and GSE104954 datasets. GSEA indicated the core genes were mainly enriched in immune-related pathways. Correlation analysis indicated LCK and HCK expressions were positively correlated with aDC, CD4+ Tem, CD8+T cells, CD8+ Tem, and mast cells.

CONCLUSIONS

We identified two immune-related genes and two immune-associated subgroups, which might help to design more precise tailored immunotherapy for DN patients.

Comprehensive analysis identifies IFI16 as a novel signature associated with overall survival and immune infiltration of skin cutaneous melanoma

Background

Skin cutaneous melanoma (SKCM) is the most common skin tumor with high mortality. The unfavorable outcome of SKCM urges the discovery of prognostic biomarkers for accurate therapy. The present study aimed to explore novel prognosis-related signatures of SKCM and determine the significance of immune cell infiltration in this pathology.

Methods

Four gene expression profiles (GSE130244, GSE3189, GSE7553 and GSE46517) of SKCM and normal skin samples were retrieved from the GEO database. Differentially expressed genes (DEGs) were then screened, and the feature genes were identified by the LASSO regression and Boruta algorithm. Survival analysis was performed to filter the potential prognostic signature, and GEPIA was used for preliminary validation. The area under the receiver operating characteristic curve (AUC) was obtained to evaluate discriminatory ability. The Gene Set Variation Analysis (GSVA) was performed, and the composition of the immune cell infiltration in SKCM was estimated using CIBERSORT. At last, paraffin-embedded specimens of primary SKCM and normal skin tissues were collected, and the signature was validated by fluorescence in situ hybridization (FISH) and immunohistochemistry (IHC).

Results

Totally 823 DEGs and 16 feature genes were screened. IFI16 was identified as the signature associated with overall survival of SKCM with a great discriminatory ability (AUC > 0.9 for all datasets). GSVA noticed that IFI16 might be involved in apoptosis and ultraviolet response in SKCM, and immune cell infiltration of IFI16 was evaluated. At last, FISH and IHC both validated the differential expression of IFI16 in SKCM.

Conclusions

In conclusion, our comprehensive analysis identified IFI16 as a signature associated with overall survival and immune infiltration of SKCM, which may play a critical role in the occurrence and development of SKCM.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02409-6.

Integrated weighted gene co-expression network analysis uncovers STAT1(signal transducer and activator of transcription 1) and IFI44L (interferon-induced protein 44-like) as key genes in pulmonary arterial hypertension

Despite the multiple diagnostic and therapeutic strategies implemented in clinical practice, the mortality rate of patients with pulmonary arterial hypertension (PAH) remains high. Understanding the mechanisms and key genes involved could provide insight into the drivers of the pathogenesis of PAH. In this research, we aimed to examine the mechanisms underlying PAH and identify key genes with potential usefulness as clinical biomarkers of PAH and thereby establish therapeutic targets for PAH. The datasets GSE117261, GSE113439, and GSE53408 were downloaded from the Gene Expression Omnibus (GEOs) database. We used weighted gene coexpression network analysis (WGCNA) to identify networks and the most relevant modules in PAH. Functional enrichment analysis was performed for the selected clinically relevant modules. The least absolute shrinkage and selection operator (LASSO) was applied to identify key genes in lung samples from patients with PAH. The genes were validated in a monocrotaline-induced PAH rat model. Three clinically relevant modules were identified through average linkage hierarchical clustering. The genes in the clinically relevant modules were related to endothelial cell differentiation, inflammation, and autoimmunity. Seven genes were screened as key genes significantly associated with PAH. Interferon-induced protein 44-like (IFI44L) and signal transducer and activator of transcription 1 (STAT1) were expressed at higher levels in the lung tissues of the PAH rat model than in those of the controls. Our findings reveal the novel pathological mechanisms underlying PAH and indicate that STAT1 and IFI44L may represent potential therapeutic targets in PAH.