In our study, a clinical factor model to classify RA-ILD stages was developed based on KL-6, and then achieved an AUC of 0.752 in the testing cohorts. This evidence concerns the gene MUC1 and rheumatoid arthritis.