MUC1 and interstitial lung disease: Complementary least absolute shrinkage and selection operator (LASSO) modeling showed that age, KL-6, and TNF-α effectively differentiated pSS–ILD (ILA 2 + ILA3) from pSS without ILD (ILA 0 + ILA 1), with an area under the curve (AUC) of 0.883 (p value < 0.0001).