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). Here, TNF is linked to interstitial lung disease.