SCODiA Report by MLBI Lab

Single-Cell Transcriptomic Profiling of Peripheral Blood in Idiopathic Pulmonary Fibrosis: Insights into Cellular Dynamics and Disease Progression

This report details the single-cell RNA sequencing analysis of peripheral blood from Idiopathic Pulmonary Fibrosis (IPF) patients and controls, revealing significant immune cell dysregulation. Key findings include a marked increase in monocyte populations in progressive IPF, accompanied by widespread pro-inflammatory and metabolic reprogramming. T cell subsets exhibit altered activation states, trafficking patterns, and signs of exhaustion, while distinct cell-cell interaction networks and surfaceome marker expressions further differentiate stable and progressive disease, highlighting systemic immune activation and potential novel therapeutic targets.

Contents

  1. Dataset overview
  2. UMAP Visualization of Single-Cell Transcriptomic Data from IPF Patients and Controls
  3. Major Cell Type Score Distribution and UMAP Visualization
  4. Marker Gene Expression Profiles Across Peripheral Blood Cell Subsets
  5. Peripheral Blood Minor Cell Type Population Analysis in Idiopathic Pulmonary Fibrosis
  6. Peripheral Blood T Cell Subset Composition in Idiopathic Pulmonary Fibrosis
  7. Peripheral Blood T cell and ILC Subsets Show Differential Proportions in Idiopathic Pulmonary Fibrosis
  8. Population Bar Plot of Monocyte Subsets
  9. Monocyte Population Dynamics in Idiopathic Pulmonary Fibrosis Progression
  10. Cell-Cell Interaction Patterns in Idiopathic Pulmonary Fibrosis Progression
  11. Condition-Specific Cell-Cell Interaction Patterns in Idiopathic Pulmonary Fibrosis
  12. T cell CD8+ Surfaceome Marker Expression in Idiopathic Pulmonary Fibrosis
  13. T cell CD4+ Surfaceome Marker Expression in Idiopathic Pulmonary Fibrosis
  14. Gene Ontology (GSA) Analysis of Monocytes in Idiopathic Pulmonary Fibrosis
  15. Gene Set Enrichment Analysis Reveals Systemic Immune Activation and Metabolic Reprogramming in Idiopathic Pulmonary Fibrosis
  16. Discussion
  17. Query List

0. Dataset overview

Dataset Summary

Precomputed Results

1. UMAP Visualization of Single-Cell Transcriptomic Data from IPF Patients and Controls

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Uniform Manifold Approximation and Projection (UMAP) plots, which are commonly used dimensionality reduction techniques to visualize single-cell RNA sequencing data. These UMAPs display the relationships between cells based on their gene expression profiles, colored by different metadata categories: condition (control, progressive_ipf, stable_ipf), sample, celltype_major, celltype_minor, and celltype_subset. The purpose of these visualizations is to assess the overall structure of the dataset, evaluate the quality of cell type annotations, identify potential batch effects, and observe condition-specific cellular distributions.

Visual Summary

Biological Interpretation

The UMAP visualizations provide a comprehensive overview of the cellular landscape in human blood samples from individuals with Idiopathic Pulmonary Fibrosis (IPF) and healthy controls.

  1. Disease-Associated Cellular Alterations: The condition UMAP suggests that while there is some overlap, cells from 'control' individuals tend to occupy distinct regions of the UMAP compared to 'progressive_ipf' and 'stable_ipf' samples. This implies that IPF, regardless of its progression state, induces systemic changes in blood cell gene expression profiles, leading to a modified cellular landscape. The intermixing of 'progressive_ipf' and 'stable_ipf' in several regions, yet some distinct enrichment for 'progressive_ipf' in others, might indicate shared underlying disease mechanisms but also unique features associated with disease progression.
  2. Robust Cell Type Annotation: The clear delineation of celltype_major, celltype_minor, and celltype_subset populations across the UMAPs demonstrates the success of the cell type annotation pipeline. Major immune cell lineages (T cells, B cells, Myeloid cells) form distinct clusters, which further resolve into functionally specialized minor cell types and highly granular subsets. This robust annotation is critical for downstream differential expression and pathway analysis, ensuring that comparisons are made between bona fide cell populations. The relatively low proportion and scattered distribution of 'unassigned' cells suggest that most cells have been confidently identified.
  3. Absence of Strong Batch Effects: The intermixing of sample origins across the UMAP indicates that the observed cellular structure is primarily driven by biological variation (cell types, conditions) rather than technical artifacts from individual sample processing. This is a crucial finding for the reliability of subsequent comparative analyses.
  4. Potential for Cell-Type Specific Disease Signatures: Given the distinct separation of cell types and the condition-specific clustering, these UMAPs provide a foundation for investigating how specific immune cell populations or their subsets are altered in IPF. For example, one could explore if certain T cell subsets (e.g., T cell (CD8 Exhausted) or T cell (Th17)) are disproportionately enriched or transcriptionally modified in progressive IPF compared to stable IPF or controls, which could shed light on disease pathogenesis.

Clinical or Translational Implications

These UMAP plots serve as an essential quality control and exploratory step in single-cell analysis.

2. Major Cell Type Score Distribution and UMAP Visualization

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the distribution of major cell types across a Uniform Manifold Approximation and Projection (UMAP) embedding of single-cell RNA sequencing data from human blood. Each plot displays the HiCAT_major_score for a specific major cell type, indicating the enrichment of marker genes associated with that cell identity. The final plot shows the assigned celltype_major annotations, serving as a reference to assess the concordance between computationally derived scores and established annotations. The primary goal is to evaluate the quality and coherence of cell type annotations and the underlying clustering structure.

Visual Summary

The UMAP embedding reveals a well-structured organization of cells, with distinct clusters reflecting different major cell lineages.

This strong agreement indicates that the major cell type annotations are robust and well-supported by their characteristic gene expression signatures.

Biological Interpretation

The UMAP visualization, complemented by major cell type scores, provides a clear representation of the cellular landscape in human peripheral blood. The distinct clustering and high concordance between HiCAT_major_score and celltype_major annotations for key immune populations (T cells, B cells, Plasma cells, Myeloid cells) underscore the robust resolution of different immune lineages based on their gene expression profiles. This clear separation reflects fundamental differences in cellular function and origin within the hematopoietic system.

The identification of MDSC scores within the broader Myeloid compartment suggests the presence of these immunosuppressive cells, which are known to be elevated in various disease states, including cancer and chronic inflammatory conditions like IPF. The low abundance of other cell types such as Mast cells, Erythroid cells, Endothelial cells, and Stromal cells aligns with their typical distribution; these cells are not expected to be major constituents of healthy peripheral blood. Their minor presence might indicate rare circulating populations, specific disease-related phenomena (e.g., circulating endothelial cells in vascular injury), or minimal contamination during sample processing.

Annotation Notes

3. Marker Gene Expression Profiles Across Peripheral Blood Cell Subsets

Report figure

[Analysis Visualization Results]...

Analysis Overview

This dot plot visualizes the expression of selected marker genes across various peripheral blood cell subsets (celltype_subset) identified in the single-cell RNA-seq data. The purpose of this analysis is to validate the assigned cell identities by examining whether the expression patterns of known canonical markers align with the cell type annotations. Each dot represents a specific gene in a given cell subset: its size indicates the fraction of cells within that group expressing the gene (non-zero expression percentage), and its color intensity (from light to dark red) reflects the mean expression level of the gene in that group. The bar plot on the right displays the total number of cells assigned to each subset.

Visual Summary

The dot plot clearly segregates distinct clusters of marker gene expression corresponding to different cell lineages and their respective subsets.

Biological Interpretation

The marker gene expression patterns generally provide robust support for the celltype_subset annotations.

Myeloid Cells:

T Cell Subsets:

Annotation Notes

The comprehensive and largely specific marker gene expression patterns observed across the various celltype_subset populations strongly support the quality and accuracy of the cell type annotations. The clear delineation of major lineages and their finer subsets based on canonical markers indicates a well-annotated dataset, which is crucial for subsequent downstream analyses related to Idiopathic Pulmonary Fibrosis. The identification of common markers and shared functional genes across related cell types also aligns with known biological relationships. This robust annotation serves as a reliable foundation for investigating disease-associated changes in IPF.

4. Peripheral Blood Minor Cell Type Population Analysis in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot of minor cell types across 'control', 'progressive_ipf', and 'stable_ipf' conditions, derived from single-cell RNA sequencing data of human blood. Each bar represents an individual sample, with stacked segments indicating the proportional representation of different minor cell types within that sample. This visualization aims to identify shifts in immune cell composition in the peripheral blood associated with Idiopathic Pulmonary Fibrosis (IPF) and its progression states.

Visual Summary

The stacked bar plots reveal distinct patterns in the relative abundance of peripheral blood minor cell types when comparing control individuals to those with stable and progressive IPF.

Biological Interpretation

The observed shifts in peripheral blood cell populations provide insights into the systemic immune alterations in IPF.

Clinical or Translational Implications

These findings suggest that peripheral blood cell composition can serve as a systemic indicator of disease activity and immune dysregulation in IPF.

5. Peripheral Blood T Cell Subset Composition in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents a population bar plot illustrating the relative proportions of T cell subsets, Natural Killer (NK) cells, and Innate Lymphoid Cells (ILCs) within the peripheral blood mononuclear cells (PBMCs) of individuals across three conditions: control, progressive idiopathic pulmonary fibrosis (IPF), and stable IPF. Each bar represents a single patient sample, showing the stacked proportions of different immune cell subsets. The goal is to identify potential shifts in T cell-related immune cell populations associated with IPF and its progression states.

Visual Summary

The visualization displays the proportional composition of T cell subsets and related immune cells (NK cells, ILCs) for each sample, grouped by condition: control, progressive_ipf, and stable_ipf.

Biological Interpretation

The observed shifts in peripheral blood immune cell populations provide insights into systemic immune responses in IPF.

Clinical or Translational Implications

These findings suggest that specific alterations in the peripheral immune cell landscape are associated with Idiopathic Pulmonary Fibrosis.

6. Peripheral Blood T cell and ILC Subsets Show Differential Proportions in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportions of various T cell and innate lymphoid cell (ILC) subsets in the peripheral blood of individuals with Idiopathic Pulmonary Fibrosis (IPF), categorized into stable IPF, progressive IPF, and control groups. Box plots illustrate the distribution of cell type proportions for each condition, and statistical comparisons (p-values) are presented to highlight significant differences between groups (p < 0.1). This helps to identify immune cell shifts in the circulation that may be associated with different stages of IPF.

Visual Summary

The box plots reveal several statistically significant differences in T cell and ILC subset proportions across the studied conditions:

Biological Interpretation

These findings highlight distinct shifts in peripheral immune cell composition associated with IPF and its progression, suggesting potential roles for these cell types in disease pathogenesis.

Clinical or Translational Implications

These differential immune cell proportions in the peripheral blood could serve as potential biomarkers for IPF diagnosis, monitoring disease progression, or predicting response to therapy. For example, the elevated Treg and Th17 cell proportions in progressive IPF suggest a dysregulated immune environment that could be targeted therapeutically. Furthermore, understanding the systemic reduction of ILC2s might provide insights into their trafficking dynamics or functional alterations in IPF, which could inform novel immunomodulatory strategies. However, further studies are required to establish the precise functional implications of these changes and their utility as clinical biomarkers.

References

  1. Th17 Cells in Fibrosis: PubMed search for "Th17 cells fibrosis" https://pubmed.ncbi.nlm.nih.gov/?term=Th17+cells+fibrosis
  2. B cells in IPF: PubMed search for "IPF B cell pathogenesis" https://pubmed.ncbi.nlm.nih.gov/?term=IPF+B+cell+pathogenesis
  3. ILC2 in Fibrosis: PubMed search for "ILC2 fibrosis" https://pubmed.ncbi.nlm.nih.gov/?term=ILC2+fibrosis

7. Population Bar Plot of Monocyte Subsets

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis visualizes the population distribution of Monocyte subsets across individual samples from different conditions: control, progressive idiopathic pulmonary fibrosis (IPF), and stable IPF. The goal was to examine whether there are shifts in the proportions of different Monocyte subsets within each sample and across disease states. The plot shows the relative proportions (percentage) of cell types categorized at the celltype_subset level, specifically focusing on cells initially identified as 'Monocyte' from the celltype_minor annotation.

Visual Summary

The bar plot displays the percentage of identified Monocyte subsets for each sample within the three conditions: control, progressive_ipf, and stable_ipf. Across all samples and all three conditions, a single bar labeled "Monocyte" is present for each sample, and it consistently reaches 100% of the population. There are no other distinct subsets visible for Monocytes in any of the samples or conditions.

Biological Interpretation

The consistent 100% representation of 'Monocyte' as the sole subset for all samples across control, progressive_ipf, and stable_ipf conditions indicates that, within this dataset's annotation scheme, Monocytes (as defined by celltype_minor) are not further resolved into finer, distinct subsets at the celltype_subset level.

While Monocytes are known to be a heterogeneous population comprising classical, intermediate, and non-classical monocytes, as well as being precursors to macrophages and dendritic cells, the current celltype_subset annotation does not provide this granularity for the cells classified as celltype_minor: Monocyte. Therefore, based on the provided annotations, there are no observable differences in the proportions of Monocyte *subsets* that can be attributed to IPF progression or stability.

Annotation Notes

This visualization primarily serves as a check on the granularity of the cell type annotations for Monocytes within this dataset. The result suggests that:

  1. The celltype_minor category 'Monocyte' is directly mapped to the celltype_subset category 'Monocyte'.
  2. If specific functional or phenotypical Monocyte subsets (e.g., CD14+ classical monocytes, CD16+ non-classical monocytes) are of interest for IPF research, deeper sub-clustering or a different annotation strategy would be required to identify and quantify these populations. The current annotation structure does not allow for differential analysis of Monocyte subsets.

8. Monocyte Population Dynamics in Idiopathic Pulmonary Fibrosis Progression

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the proportion of Monocyte cells in peripheral blood across different conditions: control, stable Idiopathic Pulmonary Fibrosis (IPF), and progressive IPF. The box plot visualizes the distribution of Monocyte cell proportions for each condition, along with statistical comparisons to identify significant differences. The goal is to understand how the relative abundance of monocytes changes with IPF status and progression.

Visual Summary

The box plot illustrates the Celltype proportion of Monocytes (labeled "Mono") across the three conditions: 'control', 'stable_ipf', and 'progressive_ipf'.

Statistical Significance:

Biological Interpretation

The observed increase in peripheral blood Monocyte proportion in IPF patients, particularly in the progressive stage, suggests a heightened systemic inflammatory and immune response. Monocytes are critical innate immune cells that differentiate into macrophages and dendritic cells, playing a central role in inflammation, tissue repair, and fibrosis.

Clinical or Translational Implications

The findings suggest that the proportion of circulating monocytes could serve as a valuable biomarker for monitoring IPF progression.

9. Cell-Cell Interaction Patterns in Idiopathic Pulmonary Fibrosis Progression

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns within peripheral blood mononuclear cells across different clinical conditions: control, progressive_ipf, and stable_ipf. The provided dot plots visualize the top 80 ligand-receptor interactions, highlighting their significance (dot size, -log10(p-value)) and interaction strength (dot color, log2(mean expression)) for various cell-cell pairs, using celltype_minor annotations. The goal is to identify distinct communicative landscapes associated with Idiopathic Pulmonary Fibrosis (IPF) and its progression.

Visual Summary

The three dot plots display the most significant cell-cell interactions for control, progressive_ipf, and stable_ipf conditions.

IPF-Associated Changes (Progressive & Stable vs. Control):

Differences between Progressive and Stable IPF:

Biological Interpretation

The observed cell-cell interaction patterns in peripheral blood reflect systemic immune dysregulation associated with Idiopathic Pulmonary Fibrosis.

Heightened Inflammatory and Chemotactic Signals in IPF:

Clinical or Translational Implications

The observed cell-cell interaction patterns have significant clinical and translational implications for IPF.

Therapeutic Target Prioritization:

10. Condition-Specific Cell-Cell Interaction Patterns in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates cell-cell interaction (CCI) patterns involving major blood cell types (T cell, Myeloid cell, Plasma cell, Megakaryocytic cell, B cell) across different Idiopathic Pulmonary Fibrosis (IPF) conditions: control, progressive_ipf, and stable_ipf. The dot plot visualizes the strength and statistical significance of specific ligand-receptor interactions between defined cell pairs, allowing for a comparison of communication networks in the peripheral blood during different disease stages of IPF.

Visual Summary

The dot plot displays cell-cell interactions, with conditions arranged as major columns (control, progressive_ipf, stable_ipf) and individual samples as rows within each condition. Each dot represents a unique ligand-receptor pair interacting between two specified cell types (e.g., ProstaglandinE2_byPTGES3_PTGER4 -- Monocyte | T cell CD8+).

Key Observations from the Plot:

  1. Condition-Specific Patterns: There are clear differences in CCI patterns across the three conditions. While some interactions are present in all conditions, their strength and significance vary considerably.
  2. Control Condition (C27-C40): This group exhibits a broad range of strong and significant interactions, particularly involving ProstaglandinE2-receptor pairs, PLAUR-integrin, and SIRPG-CD47 interactions, predominantly between Monocytes and T cells or NK cells. This suggests a robust basal immune communication network.
  3. Progressive IPF Condition (P01-P12): This condition shows a mixed pattern. Many interactions observed in the control group persist, often with comparable strength, but some appear dampened. Notably, there's a strong presence of co-stimulatory interactions like CD86-CD28 and CD40LG-CD40 involving T cells and B cells, as well as FASLG-FAS involving T cells and NK cells. PLAUR-integrin interactions also remain prominent.
  4. Stable IPF Condition (S14-S26): This group displays a distinct shift. While some general interactions are reduced, a unique cluster of highly significant Monocyte-Monocyte and Monocyte-T cell interactions emerges towards the right side of the plot. These include specific ProstaglandinE2-receptor pairs (e.g., PTGER2-Monocyte|Monocyte, PTGER4-Monocyte|Monocyte), CD47-SIRPB1 complex, NAMPT-NOX2 complex, HLA-C_KIR2DL1, and LGALS9-P4HB. This suggests a different inflammatory or regulatory immune environment in stable disease.

Biological Interpretation

The observed condition-specific CCI patterns provide insight into the systemic immune landscape in different stages of IPF.

Prostaglandin E2 (PGE2) Signaling:

Cell Adhesion and Migration:

T and B Cell Co-stimulation in Progressive IPF:

Immune Regulation and Monocyte Interactions in Stable IPF:

Mechanism:

Clinical or Translational Implications

The identification of condition-specific CCI patterns in peripheral blood offers potential avenues for understanding IPF pathogenesis and developing targeted therapies.

Therapeutic Targets:

11. T cell CD8+ Surfaceome Marker Expression in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the surfaceome marker expression profiles of CD8+ T cells in human blood samples from individuals with progressive Idiopathic Pulmonary Fibrosis (IPF), stable IPF, and healthy controls. The objective was to identify surface markers that are differentially expressed in CD8+ T cells across these conditions, comparing each IPF condition against the control group. The plot_markers_and_expression_dot tool was used to visualize these markers, focusing on surfaceome-specific genes with significant fold changes and expression scores. Up to 50 markers were considered per condition, enabling a detailed look into potential disease-associated phenotypes of CD8+ T cells.

Visual Summary

The dot plot displays the expression of selected surfaceome markers in CD8+ T cells across individual samples, grouped by condition (control, progressive_ipf, stable_ipf). The size of each dot represents the fraction of cells within that sample group expressing the marker, while the color intensity indicates the mean expression level. Samples are on the y-axis, and markers are on the x-axis, organized by the condition for which they were identified as differentially expressed against the control group. The bar plot on the right shows the total number of CD8+ T cells analyzed for each sample.

Key visual observations include:

Biological Interpretation

The differential expression of surfaceome markers on CD8+ T cells in the blood of IPF patients suggests significant alterations in their immune functions, migratory capacities, and activation states.

  1. SELL (CD62L) Upregulation:
  1. S1PR1 Upregulation:
  1. HLA-DRB5 Expression on CD8+ T cells:
  1. CD320 (Transcobalamin Receptor) Upregulation:
  1. PTGER2 (EP2 Receptor) Expression:

These observations indicate a systemic alteration of CD8+ T cell biology in IPF patients, characterized by changes in their trafficking molecules (SELL, S1PR1), activation state (HLA-DRB5), and metabolic profiles (CD320).

Clinical or Translational Implications

The identified differentially expressed surfaceome markers in CD8+ T cells offer potential insights into IPF pathogenesis and highlight avenues for diagnostic and therapeutic development.

  1. Biomarker Potential: The distinct expression patterns of SELL, S1PR1, HLA-DRB5, and CD320 could serve as circulating biomarkers for IPF diagnosis, monitoring disease progression, or distinguishing between stable and progressive forms. Flow cytometry-based assays could easily detect these surface markers on peripheral blood CD8+ T cells.
  2. Therapeutic Targets:

These findings necessitate further functional studies to ascertain the precise roles of these altered CD8+ T cell populations in IPF pathogenesis and to validate their utility as diagnostic tools or therapeutic targets.

12. T cell CD4+ Surfaceome Marker Expression in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis investigates the expression of surfaceome markers on CD4+ T cells from human blood, comparing patients with progressive Idiopathic Pulmonary Fibrosis (IPF) and stable IPF to a control group. The dot plot visualizes the mean expression level (color intensity) and the fraction of cells expressing a given marker (dot size) for each condition. Markers were selected based on differential expression criteria comparing each IPF condition against the reference control condition. The goal is to identify unique or enriched surfaceome markers that could serve as biological indicators or potential therapeutic targets for IPF.

Visual Summary

The dot plot effectively highlights condition-specific surfaceome markers for CD4+ T cells. Key observations include:

Biological Interpretation

The identified surfaceome markers on CD4+ T cells provide insights into altered immune cell function, trafficking, and signaling pathways in IPF, particularly in its progressive form.

Immune Cell Adhesion and Trafficking:

T Cell Activation and Co-signaling:

Cytokine and Prostaglandin Signaling:

Clinical or Translational Implications

The identified surfaceome markers on CD4+ T cells offer potential avenues for diagnostic and therapeutic interventions in IPF.

Therapeutic Targets for Immune Modulation:

These surfaceome markers provide a rich set of candidates for further functional validation and exploration as diagnostic tools or therapeutic targets to address the immune dysregulation and fibrosis characteristic of Idiopathic Pulmonary Fibrosis.

13. Gene Ontology (GSA) Analysis of Monocytes in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Analysis (GSA) results for Monocytes, leveraging Gene Ontology (GO) terms. The goal is to identify biological pathways and processes that are significantly enriched based on up-regulated differentially expressed genes (DEGs) in Monocytes across different Idiopathic Pulmonary Fibrosis (IPF) disease conditions, compared to control or other disease states. The analysis specifically focuses on three comparisons:

  1. control_vs_others: Monocytes from control individuals compared to both progressive_ipf and stable_ipf groups. This highlights pathways up-regulated in healthy control monocytes relative to disease states.
  2. progressive_ipf_vs_control: Monocytes from individuals with progressive IPF compared to control individuals. This identifies pathways up-regulated in monocytes during progressive disease.
  3. stable_ipf_vs_control: Monocytes from individuals with stable IPF compared to control individuals. This identifies pathways up-regulated in monocytes during stable disease.

The results are visualized as bar plots, showing the top enriched GO terms ranked by their statistical significance (-log(p-val) and -log(q-val)).

Visual Summary

GSA for Monocyte: control_vs_others

The first bar plot illustrates GO terms enriched in control Monocytes when compared to both progressive and stable IPF Monocytes. The most significantly enriched terms are primarily associated with neurodegenerative diseases (e.g., Amyotrophic lateral sclerosis, Parkinson disease, Huntington disease, Alzheimer disease, Pathways of neurodegeneration), and fundamental cellular processes such as Spliceosome, Proteasome, Protein processing in endoplasmic reticulum, Oxidative phosphorylation, Citrate cycle (TCA cycle), and RNA degradation. Some infection-related pathways like Epstein-Barr virus infection and Salmonella infection also appear. These terms, especially the core cellular machinery and metabolic pathways, suggest a baseline level of activity or homeostasis in healthy monocytes that is relatively higher or distinctly regulated compared to IPF disease states.

GSA for Monocyte: progressive_ipf_vs_control

The second bar plot highlights GO terms significantly up-regulated in Monocytes from progressive IPF patients compared to control. A striking feature is the strong enrichment of numerous infection-related pathways, with Coronavirus disease being the top hit, followed by Yersinia infection, Salmonella infection, Shigellosis, Human cytomegalovirus infection, and Leishmaniasis. Other highly enriched terms include basic cellular machinery and metabolic pathways like Ribosome, Endocytosis, Thermogenesis, Mitochondrion, and Lysosome. Key immune signaling pathways such as NOD-like receptor signaling pathway and Chemokine signaling pathway are also prominent. Osteoclast differentiation and processes related to cellular stress/dysfunction (e.g., Parkinson disease, Pathways of neurodegeneration, Oxidative phosphorylation) are also enriched.

GSA for Monocyte: stable_ipf_vs_control

The third bar plot shows GO terms up-regulated in Monocytes from stable IPF patients compared to control. Similar to progressive IPF, there is a substantial enrichment of infection-related pathways, including Salmonella infection, Yersinia infection, Human cytomegalovirus infection, Kaposi sarcoma-associated herpesvirus infection, Human immunodeficiency virus 1 infection, Leishmaniasis, Epstein-Barr virus infection, and Toxoplasmosis. Core innate immune signaling pathways are highly prominent: NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, NF-kappa B signaling pathway, Phagosome, and Chemokine signaling pathway. Notably, Lipid and atherosclerosis and PD-L1 expression and PD-1 checkpoint pathway in cancer are also highly enriched, alongside Autophagy and various neurotrophin/neurodegenerative pathways.

Biological Interpretation

The GSA results provide compelling insights into the functional state of circulating Monocytes in Idiopathic Pulmonary Fibrosis.

  1. Monocyte Activation and Innate Immunity in IPF: Both progressive and stable IPF monocytes exhibit a highly activated immune phenotype, characterized by the consistent and strong upregulation of numerous infection-related pathways and core innate immune signaling cascades (e.g., NOD-like receptor, Toll-like receptor, NF-kappa B, Chemokine signaling, Phagosome, Endocytosis). This indicates that monocytes in IPF patients are in a state of heightened alert, actively sensing and responding to perceived threats, which is a hallmark of chronic inflammatory diseases. This activation likely contributes to the ongoing inflammation and fibrotic processes observed in IPF.
  1. Metabolic Reprogramming in IPF Monocytes: Upregulation of pathways related to fundamental cellular machinery and metabolism (e.g., Ribosome, Mitochondrion, Lysosome, Oxidative phosphorylation, Thermogenesis, Lipid and atherosclerosis, Autophagy) in IPF monocytes suggests significant metabolic reprogramming. Activated immune cells, including monocytes, often undergo metabolic shifts to fuel their energetic demands for proliferation, cytokine production, and effector functions. For example, increased ribosomal activity indicates enhanced protein synthesis, while mitochondrial activity points to altered energy metabolism. This reprogramming can drive their pro-inflammatory and pro-fibrotic functions.
  2. Distinct Features in Progressive vs. Stable IPF:
  1. Baseline Monocyte Homeostasis in Control: In contrast, control monocytes display a relative enrichment of basic cellular maintenance and metabolic pathways, alongside terms related to neurodegenerative diseases. While the latter might reflect shared protein quality control or cellular stress response mechanisms, the overall profile for control suggests a focus on fundamental cellular housekeeping and metabolic efficiency compared to the actively inflammatory and reprogrammed state of IPF monocytes. The "infection" terms in control monocytes might represent a baseline surveillance state or responses to common environmental exposures.

Clinical or Translational Implications

14. Gene Set Enrichment Analysis Reveals Systemic Immune Activation and Metabolic Reprogramming in Idiopathic Pulmonary Fibrosis

Report figure

[Analysis Visualization Results]...

Analysis Overview

This analysis presents Gene Set Enrichment Analysis (GSEA) results for various cell types from peripheral blood, comparing Idiopathic Pulmonary Fibrosis (IPF) conditions (progressive_ipf, stable_ipf) against control samples. GSEA identifies biological pathways or gene sets that are statistically enriched or depleted in a given condition based on gene rankings. In the provided dot plot, each dot represents a specific pathway in a particular cell type and condition comparison. The color of the dot indicates the Normalized Enrichment Score (NES), where red signifies positive enrichment (pathway genes generally upregulated) and blue signifies negative enrichment (pathway genes generally downregulated). The size of the dot corresponds to the statistical significance (-log(p-value)), with larger dots indicating greater statistical confidence. This comprehensive visualization allows for an in-depth understanding of the systemic molecular alterations occurring in different immune cell populations during IPF, with a focus on comparisons against the healthy control state.

Visual Summary

The dot plot highlights widespread and distinct pathway enrichment patterns across various circulating immune cell types in IPF patients compared to controls.

Biological Interpretation

The GSEA results from peripheral blood provide critical insights into the systemic immunological and metabolic dysregulation associated with Idiopathic Pulmonary Fibrosis, offering a distinction between stable and progressive disease states.

Monocyte Dysregulation and Metabolic Reprogramming

Monocytes, as crucial innate immune cells, exhibit profound changes in IPF. The strong enrichment of inflammatory pathways such as "TNF signaling pathway", "Toll-like receptor signaling pathway", "JAK-STAT signaling pathway", "IL-17 signaling pathway", and "Chemokine signaling pathway" in progressive IPF monocytes signifies an escalated pro-inflammatory state. These pathways are pivotal for immune cell activation, cytokine production, and recruitment of other immune cells, collectively contributing to the chronic inflammation and tissue damage characteristic of IPF.

The contrasting enrichment of "Oxidative phosphorylation" in monocytes is particularly striking. Its downregulation in stable IPF monocytes and sharp upregulation in progressive IPF monocytes suggests a profound metabolic switch. This could reflect a transition from a less metabolically active or potentially immunosuppressive state in stable disease towards a highly energetic, activated pro-inflammatory state in progressive disease. Such metabolic reprogramming is known to dictate immune cell function and could play a role in monocyte-mediated fibrosis in the lung [[PubMed search: monocytes IPF metabolic reprogramming]](PubMed Search). The enrichment of "Necroptosis" pathway in progressive IPF monocytes further suggests that programmed necrotic cell death could contribute to sustained inflammation and tissue damage by releasing DAMPs.

T-cell Activation and Skewed Differentiation

Both CD4+ and CD8+ T cells demonstrate significant activation of critical immune signaling pathways. The enrichment of "Th1 and Th2 cell differentiation" and "Th17 cell differentiation" pathways in T cells from IPF patients, particularly those with progressive disease, indicates a skewed adaptive immune response. Th1 and Th17 cells are known pro-inflammatory T-helper subsets, while Th2 cells contribute to fibrosis through the production of profibrotic cytokines (e.g., IL-4, IL-13) [[PubMed search: Th1 Th2 Th17 cells IPF]](PubMed Search). The presence of these activated T-cell phenotypes in the circulation points to a systemic adaptive immune response that likely perpetuates lung inflammation and fibrosis.

NK Cell Responses

Natural Killer (NK) cells, essential for innate immunity and viral defense, show enrichment in "Natural killer cell mediated cytotoxicity", "TNF signaling pathway", and "Chemokine signaling pathway" in progressive IPF. This suggests an active cytotoxic and pro-inflammatory role for circulating NK cells during disease progression, potentially contributing to epithelial cell injury in the lung or engaging in complex interactions with other immune cells.

General Immune Activation Across Cell Types

The consistent enrichment of pathways such as "Cytokine-cytokine receptor interaction" across monocytes, T cells, and NK cells underscores a broad activation of the cytokine network, leading to sustained inflammation. Pathways like "Phagosome" in monocytes further highlight their active role in antigen presentation and clearance, which could be dysregulated in IPF.

Clinical or Translational Implications

The systemic immune and metabolic alterations identified through GSEA hold significant clinical and translational potential for Idiopathic Pulmonary Fibrosis.

15. Discussion

The comprehensive single-cell RNA sequencing analysis of peripheral blood provides critical insights into the systemic immune landscape in Idiopathic Pulmonary Fibrosis (IPF). UMAP visualizations revealed distinct cellular profiles distinguishing IPF patients from controls, suggesting systemic alterations in gene expression and cellular composition regardless of disease progression state. Notably, monocytes emerged as a central player, showing a significant increase in abundance in both stable and progressive IPF, with a particularly pronounced expansion in progressive disease. This monocyte expansion is accompanied by a profound metabolic reprogramming, evidenced by GSEA results showing downregulation of oxidative phosphorylation in stable IPF monocytes but a sharp upregulation in progressive IPF monocytes, indicative of a shift towards a highly energetic, pro-inflammatory state. Furthermore, GSA results highlighted the strong enrichment of numerous infection-related pathways and core innate immune signaling cascades (e.g., NOD-like receptor, Toll-like receptor, NF-κB, Chemokine signaling) in IPF monocytes, pointing to a state of heightened immune alert.

T cell populations also exhibited significant dysregulation. While overall CD4+ and CD8+ T cell proportions were reduced in IPF blood, within the T cell compartment, NK cells and CD8+ exhausted T cells showed relative increases. Box plot analyses further revealed an elevation of Tregs and Th17 cells in progressive IPF, alongside reductions in CD4+ Memory T cells, T follicular helper (Tfh) cells, Th2 cells, and ILC2s. Surfaceome marker analysis provided a deeper resolution, showing upregulation of trafficking molecules like SELL (CD62L) and S1PR1 on CD8+ T cells, suggesting altered homing. The presence of HLA-DRB5 on CD8+ T cells indicates chronic activation or an unusual T cell differentiation, while increased CD320 suggests metabolic adaptations. Similarly, CD4+ T cells in progressive IPF exhibited high expression of SELPLG (CD162), S1PR4, TRAT1, and PTGER2, collectively pointing to enhanced homing capacity, altered TCR signaling, and modified responsiveness to prostaglandin E2.

Cell-cell interaction (CCI) analyses underscored the importance of monocyte-centric communication networks. Key interactions like Zarachidonoylglycerol_by_DAGLA-FPR1 (involving monocytes, NK, and T cells) and CD40LG-CD40 (T cell-B cell/monocyte) were notably upregulated in IPF. Critically, the TGFB1-TGFBR1_complex interaction, particularly between monocytes and CD4+ T cells, was more pronounced in progressive IPF, reinforcing the central role of TGF-β in fibrosis progression. Condition-specific CCI patterns revealed unique signatures: progressive IPF showed heightened adaptive immune activity with interactions like CD86-CD28, while stable IPF featured distinct monocyte-monocyte interactions involving CD47-SIRPB1_complex, NAMPT-NOX2_complex, and the PD-L1 expression and PD-1 checkpoint pathway. This latter finding in stable IPF monocytes is particularly notable, suggesting a potential role for immune checkpoints in modulating disease stability or, conversely, contributing to immune evasion and inadequate resolution of inflammation. These systemic alterations collectively portray IPF as a disease driven by complex immune dysregulation, with circulating immune cells actively contributing to the inflammatory and pro-fibrotic milieu.

Hypotheses:

  1. Peripheral blood monocytes in progressive IPF undergo a metabolic switch, characterized by heightened oxidative phosphorylation, which fuels their increased pro-inflammatory and pro-fibrotic functions, thereby driving disease acceleration.
  2. The altered expression of homing receptors (SELL, SELPLG, S1PR1, S1PR4) and activation markers (HLA-DRB5, TRAT1) on circulating T cells in IPF reflects a dysregulated trafficking and activation state that contributes to chronic inflammation and impaired resolution of lung injury.
  3. The upregulation of the PD-L1 expression and PD-1 checkpoint pathway in monocytes of stable IPF patients signifies an active immune regulatory mechanism that helps to dampen inflammation and maintain disease stability, rather than promoting immune evasion.
  4. Persistent or recurrent viral/bacterial infections contribute to the systemic activation of innate immune pathways in IPF monocytes, potentially acting as triggers or exacerbating factors for disease progression.
  5. The dynamic and cell-type-specific alterations in prostaglandin E2 receptor expression (PTGER2) on T cells and its involvement in monocyte-monocyte interactions play a critical, context-dependent role in modulating immune responses and influencing IPF pathogenesis.
  6. Increased proportions of Tregs and Th17 cells in progressive IPF contribute to a dysregulated immune environment, where Tregs may be dysfunctional in suppressing fibrosis and Th17 cells actively promote inflammation and fibrogenesis.

Potential therapeutic targets:

  1. TGF-β signaling (TGFB1-TGFBR1 complex): The TGFB1-TGFBR1_complex interaction, particularly between monocytes and CD4+ T cells, is notably pronounced in progressive IPF. TGF-β is a master regulator of fibrosis and immunosuppression, making its pathway a critical driver of IPF pathogenesis. Evidence: Section 9 and 10 CCI analysis shows heightened TGFB1-TGFBR1_complex interactions in progressive IPF compared to controls and stable IPF. Validation: Evaluate small molecule inhibitors of TGF-β receptor (e.g., ALK5 inhibitors) in *in vitro* co-culture systems of IPF patient-derived monocytes and T cells to assess reduction in pro-fibrotic cytokine production and T cell activation. Test efficacy in relevant *in vivo* models of pulmonary fibrosis.
  2. FPR1 signaling (Zarachidonoylglycerol_by_DAGLA-FPR1): Upregulation of Zarachidonoylglycerol_by_DAGLA-FPR1 interactions in both stable and progressive IPF suggests an altered inflammatory microenvironment and enhanced immune cell activation and chemotaxis. FPR1 on phagocytes is crucial for immune cell activation and recruitment. Evidence: Section 9 CCI analysis demonstrates stronger Zarachidonoylglycerol_by_DAGLA-FPR1 interactions, particularly involving monocytes with NK and T cells, in IPF conditions. Validation: Test FPR1 antagonists in *in vitro* assays to block monocyte, NK cell, and T cell chemotaxis and activation induced by 2-arachidonoylglycerol (2-AG). Further, assess *in vivo* in fibrosis models for effects on immune cell infiltration and fibrotic lesion development.
  3. S1P receptors (S1PR1/S1PR4) on T cells: S1PR1 is upregulated on CD8+ T cells, and S1PR4 on CD4+ T cells, especially in progressive IPF. These receptors are crucial for lymphocyte egress from lymphoid organs and their trafficking, suggesting their modulation could impact T cell distribution and inflammatory processes. Evidence: Section 11 reports S1PR1 upregulation on CD8+ T cells. Section 12 shows S1PR4 upregulation on CD4+ T cells in IPF, particularly progressive IPF. Validation: Evaluate existing S1PR modulators (e.g., S1PR1 agonists like fingolimod, or specific S1PR4 modulators) in *in vitro* transwell migration assays using IPF patient-derived T cells. Assess their impact on T cell distribution, activation, and fibrosis progression in *in vivo* lung fibrosis models.
  4. PD-1/PD-L1 axis (PD-L1 on monocytes): The PD-L1 expression and PD-1 checkpoint pathway is uniquely and highly enriched in stable IPF monocytes. This suggests immune checkpoint mechanisms might be involved in maintaining disease stability or, conversely, contributing to immune dysfunction in IPF. Evidence: Section 13 GSA for monocytes comparing stable IPF vs. control explicitly highlights enrichment of PD-L1 expression and PD-1 checkpoint pathway in cancer. Validation: Investigate the effects of PD-1 or PD-L1 blocking antibodies (clinically available) on *ex vivo* IPF monocyte and T cell interactions, cytokine production, and T cell effector functions. Evaluate in relevant *in vivo* models for impact on immune activation and fibrotic responses, carefully considering the potential for unintended immune activation or exacerbation of inflammation.

Follow-up validation ideas:

  1. Validate altered cell population proportions (e.g., monocytes, Tregs, Th17, NK cells) and surface marker expression (e.g., SELL, S1PR1, HLA-DRB5, SELPLG, S1PR4) via flow cytometry or imaging mass cytometry in an independent cohort of IPF patients and controls.
  2. Conduct *in vitro* functional assays (e.g., chemotaxis, cytokine production, differentiation assays, cytotoxicity assays) using sorted monocytes and T cell subsets from IPF patients to confirm the functional consequences of altered gene expression and surface marker profiles.
  3. Perform spatial transcriptomics or multiplexed immunohistochemistry on lung tissue biopsies from IPF patients to correlate peripheral blood findings with localized immune cell infiltration, activation states, and cell-cell interaction patterns within the fibrotic lung microenvironment.
  4. Utilize targeted qPCR or Western blot to orthogonally validate the expression levels of key genes (e.g., SPP1, FPR1, CD40L, TGFB1, PD-L1, S1PR1, SELPLG) in specific sorted cell populations or lung tissue.
  5. Implement *in vitro* perturbation assays using small molecule inhibitors or blocking antibodies targeting identified pathways (e.g., S1PR modulators, TGF-β inhibitors, PD-1/PD-L1 blockers) to assess their impact on immune cell function, proliferation, differentiation, and interaction in IPF models.
  6. Conduct longitudinal cohort studies to track changes in peripheral blood cell populations, marker expression, and pathway activation over time in IPF patients, correlating these changes with disease progression, exacerbations, and response to existing anti-fibrotic therapies.

Limitations:

This study's reliance on peripheral blood single-cell RNA sequencing offers a systemic perspective but may not fully capture the intricate cellular and molecular events occurring within the primary affected lung tissue. The findings are predominantly correlative, necessitating further functional and mechanistic studies to establish causality. Cell type annotations, while robust, are based on predefined marker sets, and more granular sub-clustering might reveal additional subtle changes. Although batch effects appear minimal, inherent inter-individual variability within human cohorts warrants larger, independent validation studies. Finally, the proposed therapeutic targets are based on observational data and require rigorous preclinical and clinical validation before consideration for treatment.

16. Query List

  1. Show and save UMAPs for condition, sample, major cell type, minor cell type, and celltype_subset in two columns.
  2. Show major cell type scores on UMAP and save the result.
  3. Show and save a marker expression dot plot for celltype_subset. SET target_cell = None and var_group_rotation = 45. Leave all other arguments at their default values.
  4. Show and save a population bar plot of minor cell types.
  5. Show and save a population bar plot of T cell subsets.
  6. Show and save box plots of T cell subset populations with statistically significant differences between conditions, if any. Choose ncols based on the total number of panels.
  7. Show and save a population bar plot of Monocyte subsets.
  8. Show and save box plots of Monocyte subset populations with statistically significant differences between conditions, if any. Choose ncols based on the total number of panels.
  9. Show and save cell-cell interaction patterns. Select at most 80 cell-cell interactions per group.
  10. Find cell-cell interactions involving major blood cell types that differ significantly between conditions, show them as a dot plot, and save the result. Set max_n_items_per_group to 60. For each interaction identified under each condition, using the provided dataset context, describe the associated signaling pathways and mechanisms and, where relevant, its potential as a therapeutic target and related therapeutic agents.
  11. Show marker expression for T cell CD8+ cells in a dot plot, and save the result. Use only surfaceome markers, with at most 50 markers per condition. For each marker identified under each condition, using the provided dataset context, describe the associated signaling pathways and biological mechanisms within the blood and, where relevant, its potential as a therapeutic target for Idiopathic Pulmonary Fibrosis and related therapeutic agents.
  12. Show marker expressions for T cell CD4+ in a dot plot, and save the result. Use only surfaceome markers, with at most 40 markers per condition. For each marker identified under each condition, using the provided dataset context, describe the associated signaling pathways and biological mechanisms within the blood and, where relevant, its potential as a therapeutic target for Idiopathic Pulmonary Fibrosis and related therapeutic agents.
  13. Show Gene Ontology (GSA) analysis results for Monocytes as a bar plot, and save the result.
  14. Show Gene Set Enrichment Analysis results as a dot plot, and save it. Set the color map to RdBu_r, n_pws_to_show to 120, and target_cells to None.
↑ Top