Abstract
CircularRNAs (circRNAs) are non-coding RNAs which compete for microRNA (miRNA) binding, influencing the abundance and stability of other RNA species. Herein we have investigated the effect of circRNAs on the mir200-ZEB1 feedback loop in relationship with the aggressiveness of human melanoma cells. We first compared the level of expression of key factors in the mir200-ZEB1 feedback loop in primary human melanoma cells compared with their matching metastatic one and found a correlation between the aggressiveness of the cells and the level of expression of ZEB1 and SNAI1. We also analyzed factors in the mir200-ZEB1 feedback loop, including circZEB1, during the phenotypic switching of human melanoma cells. Our results showed a correlation between the level of ZEB1 and SNAI1 and the fraction of cancer stem cells in the population. The level of circZEB1 was, however, consistently high during the entire phenotypic transformation. To understand this result we propose a mathematical model of the regulatory circuit. According to the model, the experimental observations can be explained by the presence of a back-splicing factor limiting circRNA production.