PurposePrevious studies that identify putative genes associated with diabetic retinopathy are only focusing on specific clinical stages, thus resulting genes are not necessarily reflective of disease progression.This study identified genes associated with the severity level of diabetic retinopathy using the likelihood-ratio test (LRT) and ordinal logistic regression (OLR) model, as well as to profile immune and retinal cell landscape in progressive diabetic retinopathy using a machine learning deconvolution approach.MethodsThis study used a published transcriptomic dataset (GSE160306) from macular regions of donors with different degrees of diabetic retinopathy (10 healthy controls, 10 cases of diabetes, 9 cases of nonproliferative diabetic retinopathy, and 10 cases of proliferative diabetic retinopathy or combined with diabetic macular edema). LRT and OLR models were applied to identify severity-associated genes. In addition, CIBERSORTx was used to estimate proportional changes of immune and retinal cells in progressive diabetic retinopathy.ResultsBy controlling for gender and age using LRT and OLR, 50 genes were identified to be significantly increased in expression with the severity of diabetic retinopathy. Functional enrichment analyses suggested these severity-associated genes are related to inflammation and immune responses. CCND1 and FCGR2B are further identified as key regulators to interact with many other severity-associated genes and are crucial to inflammation. Deconvolution analyses demonstrated that the proportions of memory B cells, M2 macrophages, and Müller glia were significantly increased with the progression of diabetic retinopathy.ConclusionsThese findings demonstrate that deep analyses of transcriptomic data can advance our understanding of progressive ocular diseases, such as diabetic retinopathy, by applying LRT and OLR models as well as bulk gene expression deconvolution.