Analysis on gene modular network reveals morphogen-directed development robustness in Drosophila

Genetic robustness is an important characteristic to tolerate genetic or nongenetic perturbations and ensure phenotypic stability. Morphogens, a type of evolutionarily conserved diffusible molecules, govern tissue patterns in a direction-dependent or concentration-dependent manner by differentially regulating downstream gene expression. However, whether the morphogen-directed gene regulatory network possesses genetic robustness remains elusive. In the present study, we collected 4217 morphogen-responsive genes along A-P axis of Drosophila wing discs from the RNA-seq data, and clustered them into 12 modules. By applying mathematical model to the measured data, we constructed a gene modular network (GMN) to decipher the module regulatory interactions and robustness in morphogen-directed development. The computational analyses on asymptotical dynamics of this GMN demonstrated that this morphogen-directed GMN is robust to tolerate a majority of genetic perturbations, which has been further validated by biological experiments. Furthermore, besides the genetic alterations, we further demonstrated that this morphogen-directed GMN can well tolerate nongenetic perturbations (Hh production changes) via computational analyses and experimental validation. Therefore, these findings clearly indicate that the morphogen-directed GMN is robust in response to perturbations and is important for Drosophila to ensure the proper tissue patterning in wing disc.