A Boolean network model of human gonadal sex determination

Background: Gonadal sex determination (GSD) in humans is a complex biological process that takes place in early stages of embryonic development when the bipotential gonadal primordium (BGP) differentiates towards testes or ovaries. This decision is directed by one of two distinct pathways embedded in a GSD network activated in a population of coelomic epithelial cells, the Sertoli progenitor cells (SPC) and the granulosa progenitor cells (GPC). In males, the pathway is activated when the Sex-Determining Region Y (SRY) gene starts to be expressed, whereas in females the WNT4/beta-catenin pathway promotes the differentiation of the GPCs towards ovaries. The interactions and dynamics of the elements that constitute the GSD network are poorly understood, thus our group is interested in inferring the general architecture of this network as well as modeling the dynamic behavior of a set of genes associated to this process under wild-type and mutant conditions. Methods: We reconstructed the regulatory network of GSD with a set of genes directly associated with the process of differentiation from SPC and GPC towards Sertoli and granulosa cells, respectively. These genes are experimentally well-characterized and the effects of their deficiency have been clinically reported. We modeled this GSD network as a synchronous Boolean network model (BNM) and characterized its attractors under wild-type and mutant conditions. Results: Three attractors with a clear biological meaning were found; one of them corresponding to the currently known gene expression pattern of Sertoli cells, the second correlating to the granulosa cells and, the third resembling a disgenetic gonad. Conclusions: The BNM of GSD that we present summarizes the experimental data on the pathways for Sertoli and granulosa establishment and sheds light on the overall behavior of a population of cells that differentiate within the developing gonad. With this model we propose a set of regulatory interactions needed to activate either the SRY or the WNT4/beta-catenin pathway as well as their downstream targets, which are critical for further sex differentiation. In addition, we observed a pattern of altered regulatory interactions and their dynamics that lead to some disorders of sex development (DSD).