ABSTRACT:To understand genes, pathways, and networks related to puberty, we characterized the transcriptome of two tissues: the pituitary gland and ovaries. Samples were harvested from pre-and postpubertal Brahman heifers (same age group). Brahman heifers (Bos indicus) are older at puberty compared with Bos taurus, a productivity issue. With RNA sequencing, we identified differentially expressed (DEx) genes and important transcription factors (TF) and predicted coexpression networks. The number of DEx genes detected in the pituitary gland was 284 (P < 0.05), and VWC2L was the most DEx gene (fold change = 4.12, P = 0.01). The gene VWC2L promotes bone mineralization through transforming growth factor-β (TGFβ) signaling. Further studies of the link between bone mineralization and puberty could target VWC2L. In ovaries, 3,871 genes were DEx (P < 0.05). Four highly DEx genes were noteworthy for their function: SLC6A13 (a γ-aminobutyric acid [GABA] transporter), OXT (oxytocin), and NPY (neuropeptide Y) and its receptor NPY2R. These genes had higher ovarian expression in postpubertal heifers. The GABA and its receptors and transporters were expressed in the ovaries of many mammals, suggesting a role for this pathway beyond the brain. The OXT pathway has been known to influence the timing of puberty in rats, via modulation of GnRH. The effects of NPY at the hypothalamus, pituitary gland, and ovaries have been documented. Neuropeptide Y and its receptors are known factors in the release of GnRH, similar to OXT and GABA, although their roles in ovarian tissue are less clear. Pathways previously related to puberty such as TGFβ signaling (P = 6.71 × 10 −5 ), Wnt signaling (P = 4.1 × 10 −2 ), and peroxisome proliferator-activated receptor (PPAR) signaling (P = 4.84 × 10 −2 ) were enriched in our data set. Seven genes were identified as key TF in both tissues: HIC2, ZIC4, ZNF219, ZSCAN26, LHX1, OLIG1, and a novel gene. An ovarian subnetwork created with TF and significant ovarian DEx genes revealed five zinc fingers as regulators: ZNF507, ZNF12, ZNF512, ZNF184, and ZNF432. Recent work of hypothalamic gene expression also pointed to zinc fingers as TF for bovine puberty. Although some zinc fingers may be ubiquitously expressed, the identification of DEx genes in common across tissues points to key regulators of puberty. The hypothalamus and pituitary gland had eight DEx genes in common. The hypothalamus and ovaries had 89 DEx genes in common. The pituitary gland and ovaries had 48 DEx genes in common. Our study confirmed the complexity of puberty and suggested further investigation on genes that code zinc fingers.