RESUMEN
Organogenesis is the phase of embryonic development leading to the formation of fully functional organs. In the case of the thyroid, organogenesis starts from the endoderm and generates a multitude of closely packed independent spherical follicular units surrounded by a dense network of capillaries. Follicular organisation is unique and essential for thyroid function, i.e. thyroid hormone production. Previous in vivo studies showed that, besides their nutritive function, endothelial cells play a central role during thyroid gland morphogenesis. However, the precise mechanisms and biological parameters controlling the transformation of the multi-layered thyroid epithelial primordium into a multitude of single-layered follicles are mostly unknown. Animal studies used to improve understanding of organogenesis are costly and time-consuming, with recognised limitations. Here, we developed and used a 2-D vertex model of thyroid growth, angiogenesis and folliculogenesis, within the open-source Chaste framework. Our in silico model, based on in vivo images, correctly simulates the differential growth and proliferation of central and peripheral epithelial cells, as well as the morphogen-driven migration of endothelial cells, consistently with our experimental data. Our simulations further showed that reduced epithelial cell adhesion was critical to allow endothelial invasion and fission of the multi-layered epithelial mass. Finally, our model also allowed epithelial cell polarisation and follicular lumen formation by endothelial cell abundance and proximity. Our study illustrates how constant discussion between theoretical and experimental approaches can help us to better understand the roles of cellular movement, adhesion and polarisation during thyroid embryonic development. We anticipate that the use of in silico models like the one we describe can push forward the fields of developmental biology and regenerative medicine.