Batch Correction Research Articles

Development of a Spectral Flow Cytometry Analysis Pipeline for High-dimensional Immune Cell Characterization.

Flow cytometry is used for immune cell analysis for cell composition and function. Spectral flow cytometry allows for high-dimensional analysis of immune cells, overcoming limitations of conventional flow cytometry. However, analyzing data from large Ab panels is challenging using traditional biaxial gating strategies. We present, to our knowledge, a novel analysis pipeline to improve analysis of spectral flow cytometry. We employ this method to identify rare T cell populations in aging. We isolated splenocytes from young (2-3 mo old) and aged (18-19 mo old) female C57BL/6N mice and then stained these with a panel of 20 fluorescently labeled Abs. We performed spectral flow cytometry and then data processing and analysis using Python within a Jupyter Notebook environment to perform dimensionality reduction, batch correction, unsupervised clustering, and differential expression analysis. Our analysis of 3,776,804 T cells from 11 spleens revealed 35 distinct T cell clusters identified by surface marker expression. We observed significant differences between young and aged mice, with clusters enriched in one age group over the other. Naive, effector memory, and central memory CD8+ and CD4+ T cell subsets exhibited age-associated changes in abundance and marker expression. We also demonstrate the utility of our pipeline in a human PBMC dataset that used a 50-fluorescent color panel. By leveraging high-dimensional analysis methods, we provide insights into the immune aging process. This approach offers a robust and easily implemented analysis pipeline for spectral flow cytometry data that may facilitate the discovery of novel therapeutic targets for age-related immune dysfunction.

ScatterX： A software for fast processing of high-throughput small-angle scattering data

Abstract Scattering experiments become increasingly popular in modern scientific research, including the areas of materials, biology, chemistry, physics, etc. Besides, various types of scattering facilities have been developed recently, such as lab-based X-ray scattering equipment, national synchrotron facilities and large neutron facilities. These above-mentioned trends bring up fast-increasing data amounts of scattering data, as well as different scattering types (X-ray, neutron, laser and even microwaves). This struggles with data analysis software of the last generation, where data are usually analyzed assuming certain first-principle models and given data formats. To help researchers process and analyze scattering data more efficiently, we developed a general and model-free scattering data analysis software based on matrix operation, which has the unique advantage of high throughput scattering data processing, analysis and visualization. To maximize generality and efficiency, data processing is performed based on a three-dimensional matrix, where scattering curves are saved as matrices or vectors, rather than the traditional definition of paired values. It can not only realize image batch processing, background subtraction and correction, but also analyze data according to scattering theory and model, such as radius of gyration, fractal dimension and other physical quantities. In the aspect of visualization, the software allows the modify the color maps of two-dimensional scattering images and the gradual color variation of one-dimensional curves to suit efficient data communications. In all, this new software can work as a standalone platform for researchers to process, analyze and visualize scattering data from different research facilities without considering different file types or formats. All codes in this manuscript are open-sourced and can be easily implemented in matrix-based software, such as MATLAB, Python and Igor.

Batch effects correction in scRNA-seq based on biological-noise decoupling autoencoder and central-cross loss

Technical or biologically irrelevant differences caused by different experiments, times, or sequencing platforms can generate batch effects that mask the true biological information. Therefore, batch effects are typically removed when analyzing single-cell RNA sequencing (scRNA-seq) datasets for downstream tasks. Existing batch correction methods usually mitigate batch effects by reducing the data from different batches to a lower dimensional space before clustering, potentially leading to the loss of rare cell types. To address this problem, we introduce a novel single-cell data batch effect correction model using Biological-noise Decoupling Autoencoder (BDA) and Central-cross Loss termed BDACL. The model initially reconstructs raw data using an auto-encoder and conducts preliminary clustering. We then construct a similarity matrix and a hierarchical clustering tree to delineate relationships within and between different batches. Finally, we introduce a Central-cross Loss (CL). This loss leverages cross-entropy loss to prompt the model to better distinguish between different cluster labels. Additionally, it employs the Central Loss to encourage samples to form more compact clusters in the embedding space, thereby enhancing the consistency and interpretability of clustering results to mitigate differences between different batches. The primary innovation of this model lies in reconstructing data with an auto-encoder and gradually merging smaller clusters into larger ones using a hierarchical clustering tree. By using reallocated cluster labels as training labels and employing the Central-cross Loss, the model effectively eliminates batch effects in an unsupervised manner. Compared to current methods, BDACL can mitigate batch effects without losing rare cell types.

DGSIST: Clustering spatial transcriptome data based on deep graph structure Infomax

Although spatial transcriptomics data provide valuable insights into gene expression profiles and the spatial structure of tissues, most studies rely solely on gene expression information, underutilizing the spatial data. To fully leverage the potential of spatial transcriptomics and graph neural networks, the DGSI (Deep Graph Structure Infomax) model is proposed. This innovative graph data processing model uses graph convolutional neural networks and employs an unsupervised learning approach. It maximizes the mutual information between graph-level and node-level representations, emphasizing flexible sampling and aggregation of nodes and their neighbors. This effectively captures and incorporates local information from nodes into the overall graph structure. Additionally, this paper developed the DGSIST framework, an unsupervised cell clustering method that integrates the DGSI model, SVD dimensionality reduction algorithm, and k-means++ clustering algorithm. This aims to identify cell types accurately. DGSIST fully uses spatial transcriptomics data and outperforms existing methods in accuracy. Demonstrations of DGSIST’s capability across various tissue types and technological platforms have shown its effectiveness in accurately identifying spatial domains in multiple tissue sections. Compared to other spatial clustering methods, DGSIST excels in cell clustering and effectively eliminates batch effects without needing batch correction. DGSIST excels in spatial clustering analysis, spatial variation identification, and differential gene expression detection and directly applies to graph analysis tasks, such as node classification, link prediction, or graph clustering. Anticipation lies in the contribution of the DGSIST framework to a deeper understanding of the spatial organizational structures of diseases such as cancer.

7864 Accurate Machine Learning-Based Diagnosis of At-Risk MASH and MASLD Subtypes Using Categorical Gradient Boosting and Select Clinical, Biochemical, and Metabolomic Measurements: Building Highly Robust Models With Few Variables Through a Multinational, Multi-Center, Biopsy-Proven Cohort

Abstract Disclosure: K. Stefanakis: None. E. Axarloglou: None. C.S. Mantzoros: None. Metabolic dysfunction-associated steatotic liver disease (MASLD) is characterized by intrahepatic lipid accumulation, may precipitate steatohepatitis (MASH) and liver fibrosis, and is a significant cardiometabolic risk factor. Liver biopsy, though invasive and costly, is the diagnostic gold standard, thus creating an unmet clinical need for accurate and reliable non-invasive diagnostic methods to diagnose MASLD, particularly the combination of MASH and significant fibrosis, classified as “at-risk MASH” and known to markedly increase morbidity and mortality. We aimed to leverage novel machine-learning (ML) algorithms to create lightweight but highly precise models that may accurately diagnose at-risk MASH and related histological outcomes through a top-down approach. This is a study of 443 participants from two Gastroenterology-Hepatology Clinics (Greece, Australia) and one Metabolic Surgery Clinic (Italy), including healthy controls and MASLD patients across the entire spectrum of the disease. All participants underwent liver biopsies, anthropometric, standard biochemical and additional hormonal measurements, and global serum metabolomics using ultra-high-performance liquid chromatography-tandem mass spectrometry, followed by robust batch correction and quality control, eventually identifying 687 known and 152 unknown metabolites. The dataset was randomly divided into a discovery (n=353) and a validation cohort (n=90) ensuring a balanced distribution of at-risk MASH and clinic. Our main ML pipeline implemented Categorical Gradient Boosting Machines using the CatBoost package for R. We pre-specified domain knowledge variables and implemented recursive feature elimination, among additional techniques, to build models through an iterative, algorithm-based process. The models were tested in the validation and entire cohort using k-fold and leave-one-out cross-validation with random resampling, and further assessed on relevant subgroups (MASLD status, obesity, diabetes, clinic, country). The main ML model consists of three routine clinical and biochemical variables, supplemented by two aminoacid metabolites (full model to be shown), and can detect the presence of at-risk MASH with a mean AUROC of 93% and mean accuracy of 92% in the validation cohort, and AUCS>93% with very high sensitivities and specificities per J-statistic in all subgroups. We likewise present ML models for the detection of MASLD, MASH, and fibrosis grades, with AUCs of over 90% in the validation cohort and all subgroups. Our models significantly outperform all applicable biomarker-based indices per DeLong’s tests. Further enhancement with lipidomic and proteomic variables will fully delineate the circulating biochemical snapshot of MASLD and likely lead to even more accurate, similarly structured models. Presentation: 6/1/2024

BatchFLEX: feature-level equalization of X-batch.

Integrative analysis of heterogeneous expression data remains challenging due to variations in platform, RNA quality, sample processing, and other unknown technical effects. Selecting the approach for removing unwanted batch effects can be a time-consuming and tedious process, especially for more biologically focused investigators. Here, we present BatchFLEX, a Shiny app that can facilitate visualization and correction of batch effects using several established methods. BatchFLEX can visualize the variance contribution of a factor before and after correction. As an example, we have analyzed ImmGen microarray data and enhanced its expression signals that distinguishes each immune cell type. Moreover, our analysis revealed the impact of the batch correction in altering the gene expression rank and single-sample GSEA pathway scores in immune cell types, highlighting the importance of real-time assessment of the batch correction for optimal downstream analysis. Our tool is available through Github https://github.com/shawlab-moffitt/BATCH-FLEX-ShinyApp with an online example on Shiny.io https://shawlab-moffitt.shinyapps.io/batch_flex/.

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.

A Framework for Quality Control in Quantitative Proteomics.

A thorough evaluation of the quality, reproducibility, and variability of bottom-up proteomics data is necessary at every stage of a workflow, from planning to analysis. We share vignettes applying adaptable quality control (QC) measures to assess sample preparation, system function, and quantitative analysis. System suitability samples are repeatedly measured longitudinally with targeted methods, and we share examples where they are used on three instrument platforms to identify severe system failures and track function over months to years. Internal QCs incorporated at the protein and peptide levels allow our team to assess sample preparation issues and to differentiate system failures from sample-specific issues. External QC samples prepared alongside our experimental samples are used to verify the consistency and quantitative potential of our results during batch correction and normalization before assessing biological phenotypes. We combine these controls with rapid analysis (Skyline), longitudinal QC metrics (AutoQC), and server-based data deposition (PanoramaWeb). We propose that this integrated approach to QC is a useful starting point for groups to facilitate rapid quality control assessment to ensure that valuable instrument time is used to collect the best quality data possible. Data are available on Panorama Public and ProteomeXchange under the identifier PXD051318.

ImmCellTyper facilitates systematic mass cytometry data analysis for deep immune profiling.

Mass cytometry is a cutting-edge high-dimensional technology for profiling marker expression at the single-cell level, advancing clinical research in immune monitoring. Nevertheless, the vast data generated by cytometry by time-of-flight (CyTOF) poses a significant analytical challenge. To address this, we describe ImmCellTyper (https://github.com/JingAnyaSun/ImmCellTyper), a novel toolkit for CyTOF data analysis. This framework incorporates BinaryClust, an in-house developed semi-supervised clustering tool that automatically identifies main cell types. BinaryClust outperforms existing clustering tools in accuracy and speed, as shown in benchmarks with two datasets of approximately 4 million cells, matching the precision of manual gating by human experts. Furthermore, ImmCellTyper offers various visualisation and analytical tools, spanning from quality control to differential analysis, tailored to users' specific needs for a comprehensive CyTOF data analysis solution. The workflow includes five key steps: (1) batch effect evaluation and correction, (2) data quality control and pre-processing, (3) main cell lineage characterisation and quantification, (4) in-depth investigation of specific cell types; and (5) differential analysis of cell abundance and functional marker expression across study groups. Overall, ImmCellTyper combines expert biological knowledge in a semi-supervised approach to accurately deconvolute well-defined main cell lineages, while maintaining the potential of unsupervised methods to discover novel cell subsets, thus facilitating high-dimensional immune profiling.

Higher-order correction of persistent batch effects in correlation networks.

Systems biology analyses often use correlations in gene expression profiles to infer co-expression networks that are then used as input for gene regulatory network inference or to identify functional modules of co-expressed or putatively co-regulated genes. While systematic biases, including batch effects, are known to induce spurious associations and confound differential gene expression analyses (DE), the impact of batch effects on gene co-expression has not been fully explored. Methods have been developed to adjust expression values, ensuring conditional independence of mean and variance from batch or other covariates for each gene, resulting in improved fidelity of DE analysis. However, such adjustments do not address the potential for spurious differential co-expression (DC) between groups. Consequently, uncorrected, artifactual DC can skew the correlation structure, leading to the identification of false, non-biological associations, even when the input data are corrected using standard batch correction. In this work, we demonstrate the persistence of confounders in covariance after standard batch correction using synthetic and real-world gene expression data examples. We then introduce Co-expression Batch Reduction Adjustment (COBRA), a method for computing a batch-corrected gene co-expression matrix based on estimating a conditional covariance matrix. COBRA estimates a reduced set of parameters expressing the co-expression matrix as a function of the sample covariates, allowing control for continuous and categorical covariates. COBRA is computationally efficient, leveraging the inherently modular structure of genomic data to estimate accurate gene regulatory associations and facilitate functional analysis for high-dimensional genomic data. COBRA is available under the GLP3 open source license in R and Python in netZoo (https://netzoo.github.io).

Identification of a Novel Activated NK-Associated Gene Score Associated with Diagnosis and Biological Therapy Response in Ulcerative Colitis

Introduction: Natural killer (NK) cells are associated with the pathogenesis of ulcerative colitis (UC); however, their precise contributions remain unclear. The present study aimed to investigate the diagnostic value of the activated NK-associated gene (ANAG) score in UC and evaluate its predictive value in response to biological therapy. Methods: Bulk RNA-seq and scRNA-seq datasets were obtained from the Gene Expression Omnibus (GEO) and Single Cell Portal (SCP) databases. In the bulk RNA-seq, differentially expressed genes (DEGs) were screened by the “Batch correction” and “Robust rank aggregation” (RRA) methods. The immune infiltration landscape was estimated using single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT. DEGs that correlated with activated NK cells were identified as activated NK-associated genes (ANAGs). Protein-protein interaction (PPI) analysis and least absolute shrinkage and selection operator (LASSO) regression were used to screen key ANAGs and establish an ANAG score. The expression levels of the four key ANAGs were validated in human samples by real-time quantitative polymerase chain reaction (RT-qPCR) and immunofluorescence. The potential therapeutic drugs for UC were identified using the DSigDB database. Through scRNA-seq data analysis, the cell scores based on the ANAGs were calculated by “AddModuleScore” and “AUCell.” Results: Immune infiltration analysis revealed a higher abundance of activated NK cells in noninflamed UC tissues (ssGSEA, p < 0.001; CIBERSORT, p < 0.01). Fifty-four DEGs correlated with activated NK cells were identified as ANAGs. The ANAG score was established using four key ANAGs (SELP, TIMP1, MMP7, and ABCG2). The ANAG scores were significantly higher in inflamed tissues (p < 0.001) and in biological therapy nonresponders (NR) tissues before treatment (golimumab, p < 0.05; ustekinumab, p < 0.001). The ANAG score demonstrated an excellent diagnostic value (AUC = 0.979). Patients with higher ANAG scores before treatment were more likely to experience a lack of response to golimumab or ustekinumab (golimumab, p < 0.05; ustekinumab, p < 0.001). Conclusion: This study established a novel ANAG score with the ability to precisely diagnose UC and distinguish the efficacy of biological treatment.

A framework for quality control in quantitative proteomics.

A thorough evaluation of the quality, reproducibility, and variability of bottom-up proteomics data is necessary at every stage of a workflow from planning to analysis. We share vignettes applying adaptable quality control (QC) measures to assess sample preparation, system function, and quantitative analysis. System suitability samples are repeatedly measured longitudinally with targeted methods, and we share examples where they are used on three instrument platforms to identify severe system failures and track function over months to years. Internal QCs incorporated at protein and peptide-level allow our team to assess sample preparation issues and to differentiate system failures from sample-specific issues. External QC samples prepared alongside our experimental samples are used to verify the consistency and quantitative potential of our results during batch correction and normalization before assessing biological phenotypes. We combine these controls with rapid analysis (Skyline), longitudinal QC metrics (AutoQC), and server-based data deposition (PanoramaWeb). We propose that this integrated approach to QC is a useful starting point for groups to facilitate rapid quality control assessment to ensure that valuable instrument time is used to collect the best quality data possible. Data are available on Panorama Public and on ProteomeXchange under the identifier PXD051318.

Evaluating batch correction methods for image-based cell profiling

High-throughput image-based profiling platforms are powerful technologies capable of collecting data from billions of cells exposed to thousands of perturbations in a time- and cost-effective manner. Therefore, image-based profiling data has been increasingly used for diverse biological applications, such as predicting drug mechanism of action or gene function. However, batch effects severely limit community-wide efforts to integrate and interpret image-based profiling data collected across different laboratories and equipment. To address this problem, we benchmark ten high-performing single-cell RNA sequencing (scRNA-seq) batch correction techniques, representing diverse approaches, using a newly released Cell Painting dataset, JUMP. We focus on five scenarios with varying complexity, ranging from batches prepared in a single lab over time to batches imaged using different microscopes in multiple labs. We find that Harmony and Seurat RPCA are noteworthy, consistently ranking among the top three methods for all tested scenarios while maintaining computational efficiency. Our proposed framework, benchmark, and metrics can be used to assess new batch correction methods in the future. This work paves the way for improvements that enable the community to make the best use of public Cell Painting data for scientific discovery.

ITGB2 related to immune cell infiltration as a potential therapeutic target of inflammatory bowel disease using bioinformatics and functional research.

Inflammatory bowel disease (IBD) is a chronic systemic inflammatory condition regarded as a major risk factor for colitis-associated cancer. However, the underlying mechanisms of IBD remain unclear. First, five GSE data sets available in GEO were used to perform 'batch correction' and Robust Rank Aggregation (RRA) to identify differentially expressed genes (DEGs). Candidate molecules were identified using CytoHubba, and their diagnostic effectiveness was predicted. The CIBERSORT algorithm evaluated the immune cell infiltration in the intestinal epithelial tissues of patients with IBD and controls. Immune cell infiltration in the IBD and control groups was determined using the least absolute shrinkage selection operator algorithm and Cox regression analysis. Finally, a total of 51 DEGs were screened, and nine hub genes were identified using CytoHubba and Cytoscape. GSE87466 and GSE193677 were used as extra data set to validate the expression of the nine hub genes. CD4-naïve T cells, gamma-delta T cells, M1 macrophages and resting dendritic cells (DCs) are the main immune cell infiltrates in patients with IBD. Signal transducer and activator of transcription 1, CCR5 and integrin subunit beta 2 (ITGB2) were significantly upregulated in the IBD mouse model, and suppression of ITGB2 expression alleviated IBD inflammation in mice. Additionally, the expression of ITGB2 was upregulated in IBD-associated colorectal cancer (CRC). The silence of ITGB2 suppressed cell proliferation and tumour growth invitro and invivo. ITGB2 resting DCs may provide a therapeutic strategy for IBD, and ITGB2 may be a potential diagnostic marker for IBD-associated CRC.

Identification of diagnostic biomarkers of rheumatoid arthritis based on machine learning-assisted comprehensive bioinformatics and its correlation with immune cells

BackgroundRheumatoid arthritis (RA) is a chronic systemic autoimmune disease characterized by inflammatory cell infiltration, which can lead to chronic disability, joint destruction and loss of function. At present, the pathogenesis of RA is still unclear. The purpose of this study is to explore the potential biomarkers and immune molecular mechanisms of rheumatoid arthritis through machine learning-assisted bioinformatics analysis, in order to provide reference for the early diagnosis and treatment of RA disease. MethodsRA gene chips were screened from the public gene GEO database, and batch correction of different groups of RA gene chips was performed using Strawberry Perl. DEGs were obtained using the limma package of R software, and functional enrichment analysis such as gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), disease ontology (DO), and gene set (GSEA) were performed. Three machine learning methods, least absolute shrinkage and selection operator regression (LASSO), support vector machine recursive feature elimination (SVM-RFE) and random forest tree (Random Forest), were used to identify potential biomarkers of RA. The validation group data set was used to verify and further confirm its expression and diagnostic value. In addition, CIBERSORT algorithm was used to evaluate the infiltration of immune cells in RA and control samples, and the correlation between confirmed RA diagnostic biomarkers and immune cells was analyzed. ResultsThrough feature screening, 79 key DEGs were obtained, mainly involving virus response, Parkinson's pathway, dermatitis and cell junction components. A total of 29 hub genes were screened by LASSO regression, 34 hub genes were screened by SVM-RFE, and 39 hub genes were screened by Random Forest. Combined with the three algorithms, a total of 12 hub genes were obtained. Through the expression and diagnostic value verification in the validation group data set, 7 genes that can be used as diagnostic biomarkers for RA were preliminarily confirmed. At the same time, the correlation analysis of immune cells found that γδT cells, CD4+ memory activated T cells, activated dendritic cells and other immune cells were positively correlated with multiple RA diagnostic biomarkers, CD4+ naive T cells, regulatory T cells and other immune cells were negatively correlated with multiple RA diagnostic biomarkers. ConclusionsThe results of novel characteristic gene analysis of RA showed that KYNU, EVI2A, CD52, C1QB, BATF, AIM2 and NDC80 had good diagnostic and clinical value for the diagnosis of RA, and were closely related to immune cells. Therefore, these seven DEGs may become new diagnostic markers and immunotherapy markers for RA.

Comprehensive assessment of machine learning methods for diagnosing gastrointestinal diseases through whole metagenome sequencing data

ABSTRACT The gut microbiome, linked significantly to host diseases, offers potential for disease diagnosis through machine learning (ML) pipelines. These pipelines, crucial in modeling diseases using high-dimensional microbiome data, involve selecting profile modalities, data preprocessing techniques, and classification algorithms, each impacting the model accuracy and generalizability. Despite whole metagenome shotgun sequencing (WMS) gaining popularity for human gut microbiome profiling, a consensus on the optimal methods for ML pipelines in disease diagnosis using WMS data remains elusive. Addressing this gap, we comprehensively evaluated ML methods for diagnosing Crohn’s disease and colorectal cancer, using 2,553 fecal WMS samples from 21 case-control studies. Our study uncovered crucial insights: gut-specific, species-level taxonomic features proved to be the most effective for profiling; batch correction was not consistently beneficial for model performance; compositional data transformations markedly improved the models; and while nonlinear ensemble classification algorithms typically offered superior performance, linear models with proper regularization were found to be more effective for diseases that are linearly separable based on microbiome data. An optimal ML pipeline, integrating the most effective methods, was validated for generalizability using holdout data. This research offers practical guidelines for constructing reliable disease diagnostic ML models with fecal WMS data.

Exploration and verification a 13-gene diagnostic framework for ulcerative colitis across multiple platforms via machine learning algorithms

Ulcerative colitis (UC) is a chronic inflammatory bowel disease with intricate pathogenesis and varied presentation. Accurate diagnostic tools are imperative to detect and manage UC. This study sought to construct a robust diagnostic model using gene expression profiles and to identify key genes that differentiate UC patients from healthy controls. Gene expression profiles from eight cohorts, encompassing a total of 335 UC patients and 129 healthy controls, were analyzed. A total of 7530 gene sets were computed using the GSEA method. Subsequent batch correction, PCA plots, and intersection analysis identified crucial pathways and genes. Machine learning, incorporating 101 algorithm combinations, was employed to develop diagnostic models. Verification was done using four external cohorts, adding depth to the sample repertoire. Evaluation of immune cell infiltration was undertaken through single-sample GSEA. All statistical analyses were conducted using R (Version: 4.2.2), with significance set at a P value below 0.05. Employing the GSEA method, 7530 gene sets were computed. From this, 19 intersecting pathways were discerned to be consistently upregulated across all cohorts, which pertained to cell adhesion, development, metabolism, immune response, and protein regulation. This corresponded to 83 unique genes. Machine learning insights culminated in the LASSO regression model, which outperformed others with an average AUC of 0.942. This model's efficacy was further ratified across four external cohorts, with AUC values ranging from 0.694 to 0.873 and significant Kappa statistics indicating its predictive accuracy. The LASSO logistic regression model highlighted 13 genes, with LCN2, ASS1, and IRAK3 emerging as pivotal. Notably, LCN2 showcased significantly heightened expression in active UC patients compared to both non-active patients and healthy controls (P < 0.05). Investigations into the correlation between these genes and immune cell infiltration in UC highlighted activated dendritic cells, with statistically significant positive correlations noted for LCN2 and IRAK3 across multiple datasets. Through comprehensive gene expression analysis and machine learning, a potent LASSO-based diagnostic model for UC was developed. Genes such as LCN2, ASS1, and IRAK3 hold potential as both diagnostic markers and therapeutic targets, offering a promising direction for future UC research and clinical application.

Macrophages as determinants and regulators of systemic sclerosis-related interstitial lung disease

BackgroundInterstitial lung disease (ILD) is the primary cause of mortality in systemic sclerosis (SSc), an autoimmune disease characterized by tissue fibrosis. SSc-related ILD (SSc-ILD) occurs more frequently in females aged 30–55 years, whereas idiopathic pulmonary fibrosis (IPF) is more prevalent in males aged 60–75 years. SSc-ILD occurs earlier than IPF and progresses rapidly. FCN1, FABP4, and SPP1 macrophages are involved in the pathogenesis of lung fibrosis; SPP1 macrophages demonstrate upregulated expression in both SSc-ILD and IPF. To identify the differences between SSc-ILD and IPF using single-cell analysis, clarify their distinct pathogeneses, and propose directions for prevention and treatment.MethodsWe performed single-cell RNA sequencing on NCBI Gene Expression Omnibus (GEO) databases GSE159354 and GSE212109, and analyzed lung tissue samples across healthy controls, IPF, and SSc-ILD. The primary measures were the filtered genes integrated with batch correction and annotated cell types for distinguishing patients with SSc-ILD from healthy controls. We proposed an SSc-ILD pathogenesis using cell–cell interaction inferences, and predicted transcription factors regulating target genes using SCENIC. Drug target prediction of the TF gene was performed using Drug Bank Online.ResultsA subset of macrophages activates the MAPK signaling pathway under oxidative stress. Owing to the lack of inhibitory feedback from ANNEXIN and the autoimmune characteristics, this leads to an earlier onset of lung fibrosis compared to IPF. During initial lung injury, fibroblasts begin to activate the IL6 pathway under the influence of SPP1 alveolar macrophages, but IL6 appears unrelated to other inflammatory and immune cells. This may explain why tocilizumab (an anti-IL6-receptor antibody) only preserves lung function in patients with early SSc-ILD. Finally, we identified BCLAF1 and NFE2L2 as influencers of MAPK activation in macrophages. Metformin downregulates NFE2L2 and could serve as a repurposed drug candidate.ConclusionsSPP1 alveolar macrophages play a role in the profibrotic activity of IPF and SSc-ILD. However, SSc-ILD is influenced by autoimmunity and oxidative stress, leading to the continuous activation of MAPK in macrophages. This may result in an earlier onset of lung fibrosis than in IPF. Such differences could serve as potential research directions for early prevention and treatment.

Development of a Spectral Flow Cytometry Analysis Pipeline for High-Dimensional Immune Cell Characterization.

Flow cytometry is a widely used technique for immune cell analysis, offering insights into cell composition and function. Spectral flow cytometry allows for high-dimensional analysis of immune cells, overcoming limitations of conventional flow cytometry. However, analyzing data from large antibody panels can be challenging using traditional bi-axial gating strategies. Here, we present a novel analysis pipeline designed to improve analysis of spectral flow cytometry. We employ this method to identify rare T cell populations in aging. We isolated splenocytes from young (2-3 months) and aged (18-19 months) female mice then stained these with a panel of 20 fluorescently labeled antibodies. Spectral flow cytometry was performed, followed by data processing and analysis using Python within a Jupyter Notebook environment to perform batch correction, unsupervised clustering, dimensionality reduction, and differential expression analysis. Our analysis of 3,776,804 T cells from 11 spleens revealed 34 distinct T cell clusters identified by surface marker expression. We observed significant differences between young and aged mice, with certain clusters enriched in one age group over the other. Naïve, effector memory, and central memory CD8+ and CD4+ T cell subsets exhibited age-associated changes in abundance and marker expression. Additionally, γδ T cell clusters showed differential abundance between age groups. By leveraging high-dimensional analysis methods borrowed from single-cell RNA sequencing analysis, we identified age-related differences in T cell subsets, providing insights into the immune aging process. This approach offers a robust, free, and easily implemented analysis pipeline for spectral flow cytometry data that may facilitate the discovery of novel therapeutic targets for age-related immune dysfunction.

Identification of potential biomarkers and immunoregulatory mechanisms of rheumatoid arthritis based on multichip co-analysis of GEO database

To identify the biomarkers for early rheumatoid arthritis (RA) diagnosis and explore the possible immune regulatory mechanisms. The differentially expressed genesin RA were screened and functionally annotated using the limma, RRA, batch correction, and clusterProfiler. The protein-protein interaction network was retrieved from the STRING database, and Cytoscape 3.8.0 and GeneMANIA were used to select the key genes and predicting their interaction mechanisms. ROC curves was used to validate the accuracy of diagnostic models based on the key genes. The disease-specific immune cells were selected via machine learning, and their correlation with the key genes were analyzed using Corrplot package. Biological functions of the key genes were explored using GSEA method. The expression of STAT1 was investigated in the synovial tissue of rats with collagen-induced arthritis (CIA). We identified 9 core key genes in RA (CD3G, CD8A, SYK, LCK, IL2RG, STAT1, CCR5, ITGB2, and ITGAL), which regulate synovial inflammation primarily through cytokines-related pathways. ROC curve analysis showed a high predictive accuracy of the 9 core genes, among which STAT1 had the highest AUC (0.909). Correlation analysis revealed strong correlations of CD3G, ITGAL, LCK, CD8A, and STAT1 with disease-specific immune cells, and STAT1 showed the strongest correlation with M1-type macrophages (R=0.68, P=2.9e-08). The synovial tissues of the ankle joints of CIA rats showed high expressions of STAT1 and p-STAT1 with significant differential expression of STAT1 between the nucleus and the cytoplasm of the synovial fibroblasts. The protein expressions of p-STAT1 and STAT1 in the cell nuclei were significantly reduced after treatment. CD3G, CD8A, SYK, LCK, IL2RG, STAT1, CCR5, ITGB2, and ITGAL may serve as biomarkers for early diagnosis of RA. Gene-immune cell pathways such as CD3G/CD8A/LCK-γδ T cells, ITGAL-Tfh cells, and STAT1-M1-type macrophages may be closely related with the development of RA.