A Strategy to Quantify Myofibroblast Activation on a Continuous Spectrum

Abstract

Myofibroblasts are a highly secretory and contractile phenotype most commonly identified by the de novo expression and assembly of alpha-smooth muscle actin stress fibers. Traditionally, this activation process has been thought of as a binary process, with cells being labeled as 'activated' or 'quiescent (non-activated)'. More recently, this view has been expanded to consider activation on a continuous spectrum. However, there is no established method to quantify the position of a cell on this spectrum, and as a result, the binary labeling system is still widely used. While transcriptomic analyses provide a continuous measure of myofibroblast markers, a faster and more facile screening method is needed. To this end, we utilized optical microscopy and machine learning methods to quantify myofibroblast activation on a spectrum. We first measured size and shape features of over 1,000 individual cardiac fibroblasts and found that these features provide enough information to predict activation state, on the binary scale, with 94% accuracy as compared to manual classification. We next performed dimensionality reduction techniques on these features to create a continuous scale of activation. Importantly, this new classification system captures a range of fibroblast activation states, but still possesses inherent bias due to choice of morphological features. Thus, we next used self-supervised machine learning to create a second continuous labeling system free from biases associated with the manually measured features. Lastly, we compared our findings for mechanically activated cardiac fibroblasts to a distribution of cell phenotypes generated from transcriptomic data using single-cell RNA sequencing. Altogether, these results demonstrate a continuous spectrum of activation from fibroblast to myofibroblast and provide a strategy to quantify the position of a cell on that spectrum.