Gene expression profiling of puberty-associated genes reveals abundant tissue and sex-specific changes across postnatal development.

The timing of human puberty is highly variable, sexually dimorphic, and associated with adverse health outcomes. Over 20 genes carrying rare mutations have been identified in known pubertal disorders, many of which encode critical components of the hypothalamic-pituitary-gonadal (HPG) axis. Recent genome-wide association studies (GWAS) have identified more than 100 candidate genes at loci associated with age at menarche or voice breaking in males. We know little about the spatial, temporal or postnatal expression patterns of the majority of these puberty-associated genes. Using a high-throughput and sensitive microfluidic quantitative PCR strategy, we profiled the gene expression patterns of the mouse orthologs of 178 puberty-associated genes in male and female mouse HPG axis tissues, the pineal gland, and the liver at five postnatal ages spanning the pubertal transition. The most dynamic gene expression changes were observed prior to puberty in all tissues. We detected known and novel tissue-enhanced gene expression patterns, with the hypothalamus expressing the largest number of the puberty-associated genes. Notably, over 40 puberty-associated genes in the pituitary gland showed sex-biased gene expression, most of which occurred peri-puberty. These sex-biased genes included the orthologs of candidate genes at GWAS loci that show sex-discordant effects on pubertal timing. Our findings provide new insight into the expression of puberty-associated genes and support the possibility that the pituitary plays a role in determining sex differences in the timing of puberty.