The divergence of distinct cell populations from multipotent progenitors is poorly understood, particularly in vivo. The gonad is an ideal place to study this process, because it originates as a bipotential primordium where multiple distinct lineages acquire sex-specific fates as the organ differentiates as a testis or an ovary. To gain a more detailed understanding of the process of gonadal differentiation at the level of the individual cell populations, we conducted microarrays on sorted cells from XX and XY mouse gonads at three time points spanning the period when the gonadal cells transition from sexually undifferentiated progenitors to their respective sex-specific fates. We analyzed supporting cells, interstitial/stromal cells, germ cells, and endothelial cells. This work identified genes specifically depleted and enriched in each lineage as it underwent sex-specific differentiation. We determined that the sexually undifferentiated germ cell and supporting cell progenitors showed lineage priming. We found that germ cell progenitors were primed with a bias toward the male fate. In contrast, supporting cells were primed with a female bias, indicative of the robust repression program involved in the commitment to XY supporting cell fate. This study provides a molecular explanation reconciling the female default and balanced models of sex determination and represents a rich resource for the field. More importantly, it yields new insights into the mechanisms by which different cell types in a single organ adopt their respective fates.

Complex patterns of cell-type-specific gene expression are thought to be achieved by combinatorial binding of transcription factors (TFs) to sequence elements in regulatory regions. Predicting cell-type-specific expression in mammals has been hindered by the oftentimes unknown location of distal regulatory regions. To alleviate this bottleneck, we used DNase-seq data from 19 diverse human cell types to identify proximal and distal regulatory elements at genome-wide scale. Matched expression data allowed us to separate genes into classes of cell-type-specific up-regulated, down-regulated, and constitutively expressed genes. CG dinucleotide content and DNA accessibility in the promoters of these three classes of genes displayed substantial differences, highlighting the importance of including these aspects in modeling gene expression. We associated DNase I hypersensitive sites (DHSs) with genes, and trained classifiers for different expression patterns. TF sequence motif matches in DHSs provided a strong performance improvement in predicting gene expression over the typical baseline approach of using proximal promoter sequences. In particular, we achieved competitive performance when discriminating up-regulated genes from different cell types or genes up-and down-regulated under the same conditions. We identified previously known and new candidate cell-type-specific regulators. The models generated testable predictions of activating or repressive functions of regulators. DNase I footprints for these regulators were indicative of their direct binding to DNA. In summary, we successfully used information of open chromatin obtained by a single assay, DNase-seq, to address the problem of predicting cell-type-specific gene expression in mammalian organisms directly from regulatory sequence.

Temperature dependent sex determination (TSD) is the process by which the environmental temperature experienced during embryogenesis influences the sex of an organism, as in the red-eared slider turtle Trachemys scripta elegans. In accord with current paradigms of vertebrate sex determination, temperature is believed to exert its effects on sexual development in T. scripta entirely within the middle third of development, when the gonad is forming. However, whether temperature regulates the transcriptome in T. scripta early embryos in a manner that could influence secondary sex characteristics or establish a pro-male or pro-female environment has not been investigated. In addition, apart from a handful of candidate genes, very little is known about potential similarities between the expression cascade during TSD and the genetic cascade that drives mammalian sex determination. Here, we conducted an unbiased transcriptome-wide analysis of the effects of male- and female-promoting temperatures on the turtle embryo prior to gonad formation, and on the gonad during the temperature sensitive period. We found sexually dimorphic expression reflecting differences in steroidogenic enzymes and brain development prior to gonad formation. Within the gonad, we mapped a cascade of differential expression similar to the genetic cascade established in mammals. Using a Hidden Markov Model based clustering approach, we identified groups of genes that show heterochronic shifts between M. musculus and T. scripta. We propose a model in which multiple factors influenced by temperature accumulate during early gonadogenesis, and converge on the antagonistic regulation of aromatase to canalize sex determination near the end of the temperature sensitive window of development.