The size and shape of organs is species specific, and even in species in which organ size is strongly influenced by environmental cues, such as nutrition or temperature, it follows defined rules. Therefore, mechanisms must exist to ensure a tight control of organ size within a given species, while being flexible enough to allow for the evolution of different organ sizes in different species. We combined computational modeling and quantitative measurements to analyze growth control in the Drosophila eye disc. We find that the area growth rate declines inversely proportional to the increasing total eye disc area. We identify two growth laws that are consistent with the growth data and that would explain the extraordinary robustness and evolutionary plasticity of the growth process and thus of the final adult eye size. These two growth laws correspond to very different control mechanisms and we discuss how each of these laws constrains the set of candidate biological mechanisms for growth control in the Drosophila eye disc.

Patterning and growth are linked during early development and have to be tightly controlled to result in a functional tissue or organ. During the development of the Drosophila eye, this linkage is particularly clear: the growth of the eye primordium mainly results from proliferating cells ahead of the morphogenetic furrow (MF), a moving signaling wave that sweeps across the tissue from the posterior to the anterior side, that induces proliferating cells anterior to it to differentiate and become cell cycle quiescent in its wake. Therefore, final eye disc size depends on the proliferation rate of undifferentiated cells and on the speed with which the MF sweeps across the eye disc. We developed a spatio-temporal model of the growing eye disc based on the regulatory interactions controlled by the signals Decapentaplegic (Dpp), Hedgehog (Hh) and the transcription factor Homothorax (Hth) and explored how the signaling patterns affect the movement of the MF and impact on eye disc growth. We used published and new quantitative data to parameterize the model. In particular, two crucial parameter values, the degradation rate of Hth and the diffusion coefficient of Hh, were measured. The model is able to reproduce the linear movement of the MF and the termination of growth of the primordium. We further show that the model can explain several mutant phenotypes, but fails to reproduce the previously observed scaling of the Dpp gradient in the anterior compartment.

A fundamental question in developmental biology is how organ size is controlled. We have previously shown that the area growth rate in the Drosophila eye primordium declines inversely proportionally to the increase in its area. How the observed reduction in the growth rate is achieved is unknown. Here, we explore the dilution of the cytokine Unpaired (Upd) as a possible candidate mechanism. In the developing eye, upd expression is transient, ceasing at the time when the morphogenetic furrow first emerges. We confirm experimentally that the diffusion and stability of the JAK/STAT ligand Upd are sufficient to control eye disc growth via a dilution mechanism. We further show that sequestration of Upd by ectopic expression of an inactive form of the receptor Domeless (Dome) results in a substantially lower growth rate, but the area growth rate still declines inversely proportionally to the area increase. This growth rate-to-area relationship is no longer observed when Upd dilution is prevented by the continuous, ectopic expression of Upd. We conclude that a mechanism based on the dilution of the growth modulator Upd can explain how growth termination is controlled in the eye disc.

Gene regulatory networks (GRNs) represent the molecular interactions that govern the behavior of cells in tissues during development. The building and analysis of GRNs require quantitative information on gene expression from tissues. Laser Scanning Confocal Microscopy (LSCM) is commonly used to obtain such information, where immunofluorescence signal can be used as a correlate of gene expression or protein levels. However, a critical step for the extraction of this information is the segmentation of LSCM digital images. Popular segmentation algorithms are frequently based on watershed methods. Here we present an algorithm for the 3D segmentation of nuclei from LSCM (x,y,z) image stacks based on regional merging and graph contractions. This algorithm outperforms watershed methods, especially when the density of images along the z-axis is low and there is a high nuclear signal crowding. In addition, it reduces the parameterization since no filter is needed in order to prevent signal noise side effects (e.g. oversegmentation). Based on this algorithm, we developed an application (iFLIC, immunoFLuorescence Imaging Cytometry tool) for the Java Virtual Machine (JVM). The application supports basic operations for reading, writing and filtering 8-bit depth multicolor TIFF image formats, including indexed file directories (IFD), which are provided by the Java Advanced Imaging (JAI) library. It also provides with basic 3D rendering and ROI specification that make extensive use of the Java3D library. iFLIC is also a plugin based application powered by the Java Plugin Platform (JPF), so each specific operation is declared as a unique command associated to one plugin and linked to a common interface. Results from segmentation can be exported both as TIFF images and a descriptive file format (iFLIC format).

Dipterans show a striking range of eye sizes, shapes, and functional specializations. Their eye is of the compound type, the most frequent eye architecture in nature. The development of this compound eye has been most studied in Drosophila melanogaster. The early development of the Drosophila eye is under the control of a gene regulatory network of transcription factors and signaling molecules called the retinal determination gene network (RDGN). Nodes in this network have been found to be involved not only in the development of different eye types in invertebrates and vertebrates, but also of other organs. Here we have analyzed the network properties in detail. First, we have generated quantitative expression profiles for a number of the key RDGN transcription factors, at a single-cell resolution. With these profiles, and applying a correlation analysis, we revisited several of the links in the RDGN. Our study uncovers a new link, that we confirm experimentally, between the transcription factors Hth/Meis1 and Optix/Six3 and indicates that, at least during the period of eye differentiation, positive feedback regulation from Eya and Dac on the Pax6 gene Ey is not operating. From this revised RDGN we derive a simplified gene network that we model mathematically. This network integrates three basic motifs: a coherent feedforward loop, a toggle-switch and a positive autoregulation which, together with the input from the Dpp/BMP2 signaling molecule, recapitulate the gene expression profiles obtained experimentally, while ensuring a robust transition from progenitor cells into retinal precursors.