Abstract
Vertebrate sex determination and differentiation are coordinated by the activations and maintenance of reproductive transcriptional-regulatory networks (TRNs). There is considerable interest in studying the conserved design principles and functions of reproductive TRNs given that their intricate regulation is susceptible to disruption by gene mutations or exposures to exogenous endocrine disrupting chemicals (or EDCs). In this manuscript, the Boolean rules describing reproductive TRNs in humans, mice, and zebrafish, were represented as a pseudo-stoichiometric matrix model. This model mathematically described the interactions of 35 transcription factors with 21 sex determination and differentiation genes across the three species. The in silico approach of Extreme Pathway (ExPa) analysis was used to predict the extent of TRN gene activations subject to the species-specific transcriptomics data, from across various developmental life-stages. A goal of this work was to identify conserved and functional reproductive TRNs across the three species. ExPa analyses predicted the sex differentiation genes, DHH, DMRT1, and AR, to be highly active in male humans, mice, and zebrafish. Whereas FOXL2 was the most active gene in female humans and mice; and CYP19A1A in female zebrafish. These results agree with the expectation that regardless of a lack of sex determination genes in zebrafish, the TRNs responsible for canalizing male vs. female sexual differentiation are conserved with mammalian taxa. ExPa analysis therefore provides a framework with which to study the TRNs that influence the development of sexual phenotypes. And the in silico predicted conservation of sex differentiation TRNs between mammals and zebrafish identifies the piscine species as an effective in vivo model to study mammalian reproductive systems under normal or perturbed pathologies.