Immune cells are highly heterogeneous and show diverse phenotypes, but the underlying mechanism remains to be elucidated. In this study, we proposed a theoretical framework for immune cell phenotypic classification based on gene plasticity, which herein refers to expressional change or variability in response to conditions. The system contains two core points. One is that the functional subsets of immune cells can be further divided into subdivisions based on their highly plastic genes, and the other is that loss of phenotype accompanies gain of phenotype during phenotypic conversion. The first point suggests phenotypic stratification or layerability according to gene plasticity, while the second point reveals expressional compatibility and mutual exclusion during the change in gene plasticity states. Abundant transcriptome data analysis in this study from both microarray and RNA sequencing in human CD4 and CD8 single-positive T cells, B cells, natural killer cells and monocytes supports the logical rationality and generality, as well as expansibility, across immune cells. A collection of thousands of known immunophenotypes reported in the literature further supports that highly plastic genes play an important role in maintaining immune cell phenotypes and reveals that the current classification model is compatible with the traditionally defined functional subsets. The system provides a new perspective to understand the characteristics of dynamic, diversified immune cell phenotypes and intrinsic regulation in the immune system. Moreover, the current substantial results based on plasticitomics analysis of bulk and single-cell sequencing data provide a useful resource for big-data–driven experimental studies and knowledge discoveries.

Background: Calcific aortic valve disease (CAVD) is a progressive fibrocalcific disease that can be treated only through valve replacement. This study aimed to determine the role of hub genes and immune cell infiltration in CAVD progression.

Methods: In this study, bioinformatics analysis was used to identify hub genes involved in CAVD. The datasets were downloaded from the Gene Expression Omnibus (GEO) database. Gene expression differences were evaluated via pathway and Gene Ontology analyses. Weighted gene co-expression network analysis (WGCNA) and differentially expressed genes were used to screen hub genes. The CIBERSORT algorithm was used to compare immune infiltration into the calcified aortic valve based on the hub genes between high- and low-expression groups. We also performed single-cell RNA sequencing based on six different human aortic valve leaflets. The expression of hub genes was identified in human and mouse samples through quantitative real-time polymerase chain reaction (qPCR), immunohistochemistry, immunofluorescence, and ELISA, and clinical features of the patients were investigated.

Results: In total, 454 differentially expressed genes were obtained from the GEO database. WGCNA was used to find 12 co-expression modules in the Array Express database, of which one hub module (brown module) was most correlated with CAVD. Two hub genes were identified after combining the differentially expressed genes S100A8 and S100A9. Regarding these genes, the immune infiltration profiles varied between high- and low-expression groups. Compared with that in the low hub gene expression group, the high hub gene expression group had a higher proportion of activated NK cells (p < 0.01) and M1 macrophages (p < 0.05). The expression of S100A8 and S100A9 was consistent with single-gene RNA sequencing results, confirming that the expression levels of these two hub genes are significantly upregulated in patients with CAVD (p < 0.01). Furthermore, these results were verified using mouse and human samples by performing immunofluorescence, immunohistochemistry, qPCR, and ELISA analyses. Finally, the localization of S100A8 and S100A9 in monocytes and macrophages was confirmed via immunofluorescence using human aortic valves.

Conclusion: These results demonstrate that S100A8 and S100A9 are two hub genes involved in CAVD, which might play an important role in its development through immune-related signaling pathways.