MIXTURE: an improved algorithm for immune tumor microenvironment estimation based on gene expression data

RNA sequencing has proved to be an efficient high-throughput technique to robustly characterize the presence and quantity of RNA in tumor biopsies at a given time. Importantly, it can be used to computationally estimate the composition of the tumor immune infiltrate and to infer the immunological phenotypes of those cells. Given the significant impact of anti-cancer immunotherapies and the role of the associated immune tumor microenvironment (ITME) on its prognosis and therapy response, the estimation of the immune cell-type content in the tumor is crucial for designing effective strategies to understand and treat cancer. Current digital estimation of the ITME cell mixture content can be performed using different analytical tools. However, current methods tend to over-estimate the number of cell-types present in the sample, thus under-estimating true proportions, biasing the results. We developed MIXTURE, a noise-constrained recursive feature selection for support vector regression that overcomes such limitations. MIXTURE deconvolutes cell-type proportions of bulk tumor samples for both RNA microarray or RNA-Seq platforms from a leukocyte validated gene signature. We evaluated MIXTURE over simulated and benchmark data sets. It overcomes competitive methods in terms of accuracy on the true number of present cell-types and proportions estimates with increased robustness to estimation bias. It also shows superior robustness to collinearity problems. Finally, we investigated the human immune microenvironment of breast cancer, head and neck squamous cell carcinoma, and melanoma biopsies before and after anti-PD-1 immunotherapy treatment revealing associations to response to therapy which have not seen by previous methods.