Integrative analysis of COL6A3 in lupus nephritis: insights from single-cell transcriptomics and proteomics

Lupus nephritis (LN), a severe complication of systemic lupus erythematosus (SLE), presents significant challenges in patient management and treatment outcomes. The identification of novel LN-related biomarkers and therapeutic targets is critical to enhancing treatment outcomes and prognosis for patients. In this study, we analyzed single-cell expression data from LN (n=21) and healthy controls (n=3). A total of 143 differentially expressed genes were identified between the LN and control groups. Then, proteomics analysis of LN patients (n=9) and control (SLE patients without LN, n=11) revealed 55 differentially expressed genes among patients with LN and control group. We further utilizes protein-protein interaction network and functional enrichment analyses to elucidate the pivotal role of COL6A3 in key signaling pathways. Its diagnostic value is evaluate through its correlation with disease progression and renal function metrics, as well as Receiver Operating Characteristic Curve (ROC) analysis. Additionally, immunohistochemistry and qPCR experiments were performed to validate the expression of COL6A3 in LN. By comparison of single-cell and proteomics data, we discovered that COL6A3 is significantly upregulated, highlighting it as a critical biomarker of LN. Our findings emphasize the substantial involvement of COL6A3 in the pathogenesis of LN, particularly noting its expression in mesangial cells. Through comprehensive protein-protein interaction network and functional enrichment analyses, we uncovered the pivotal role of COL6A3 in key signaling pathways including integrin-mediated signaling pathways, collagen-activated signaling pathways, and ECM-receptor interaction, suggesting potential therapeutic targets. The diagnostic utility is confirmed by its correlation with disease progression and renal function metrics of the glomerular filtration rate. ROC analysis further validates the diagnostic value of COL6A3, with the area under the ROC values of 0.879 in the in-house cohort, and 0.802 and 0.915 in tubular and glomerular external cohort samples, respectively. Furthermore, immunohistochemistry and qPCR experiments were consistent with those obtained from the single-cell RNA sequencing and proteomics studies. These results proved that COL6A3 is a promising biomarker and therapeutic target, advancing personalized medicine strategies for LN.