Background: The starlet sea anemone Nematostella vectensis is a diploblastic cnidarian that expresses a set of conserved genes for gut formation during its early development. During the last decade, the spatial distribution of many of these genes has been visualized with RNA hybridization or protein immunolocalization techniques. However, due to N. vectensis' curved and changing morphology, quantification of these spatial data is problematic. A method is developed for two-dimensional gene expression quantification, which enables a numerical analysis and dynamic modeling of these spatial patterns.
Methods/result: In this work, first standardized gene expression profiles are generated from publicly available N. vectensis embryo images that display mRNA and/or protein distributions. Then, genes expressed during gut formation are clustered based on their expression profiles, and further grouped based on temporal appearance of their gene products in embryonic development. Representative expression profiles are manually selected from these clusters, and used as input for a simulation-based optimization scheme. This scheme iteratively fits simulated profiles to the selected profiles, leading to an optimized estimation of the model parameters. Finally, a preliminary gene regulatory network is derived from the optimized model parameters.
Outlook: While the focus of this study is N. vectensis, the approach outlined here is suitable for inferring gene regulatory networks in the embryonic development of any animal, thus allowing to comparatively study gene regulation of gut formation in silico across various species.