LβT2 and αT3-1 are important, widely studied cell line models for the pituitary gonadotropes that were generated by targeted tumorigenesis in transgenic mice. LβT2 cells are more mature gonadotrope precursors than αT3-1 cells. Microsatellite authentication patterns, chromosomal characteristics, and their intercellular variation have not been reported. We performed microsatellite and cytogenetic analysis of both cell types at early passage numbers. Short tandem repeat (STR) profiling was consistent with a mixed C57BL/6J × BALB/cJ genetic background, with distinct patterns for each cell type. Spectral karyotyping in αT3-1 cells revealed cell-to-cell variation in chromosome composition and pseudodiploidy. In LβT2 cells, chromosome counting and karyotyping demonstrated pseudotriploidy and high chromosomal variation among cells. Chromosome copy number variation was confirmed by single-cell DNA sequencing. Chromosomal compositions were consistent with a male sex for αT3-1 and a female sex for LβT2 cells. Among LβT2 stocks used in multiple laboratories, we detected two genetically similar but distinguishable lines via STR authentication, LβT2a and LβT2b. The two lines differed in morphological appearance, with LβT2a having significantly smaller cell and nucleus areas. Analysis of immediate early gene and gonadotropin subunit gene expression revealed variations in basal expression and responses to continuous and pulsatile GnRH stimulation. LβT2a showed higher basal levels of Egr1, Fos, and Lhb but lower Fos induction. Fshb induction reached significance only in LβT2b cells. Our study highlights the heterogeneity in gonadotrope cell line genomes and provides reference STR authentication patterns that can be monitored to improve experimental reproducibility and facilitate comparisons of results within and across laboratories.

Despite their importance in tissue homeostasis and renewal, human pituitary stem cells (PSCs) are incompletely characterized. We describe a human single nucleus RNA-seq and ATAC-seq resource from pediatric, adult, and aged postmortem pituitaries (snpituitaryatlas.princeton.edu) and characterize cell-type-specific gene expression and chromatin accessibility programs for all major pituitary cell lineages. We identify uncommitted PSCs, committing progenitor cells, and sex differences. Pseudotime trajectory analysis indicates that early-life PSCs are distinct from the other age groups. Linear modeling of same-cell multiome data identifies regulatory domain accessibility sites and transcription factors that are significantly associated with gene expression in PSCs compared with other cell types and within PSCs. We identify distinct deterministic mechanisms that contribute to heterogeneous marker expression within PSCs. These findings characterize human stem cell lineages and reveal diverse mechanisms regulating key PSC genes and cell type identity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

To provide a multi-omics resource and investigate transcriptional regulatory mechanisms, we profile the transcriptome, chromatin accessibility, and methylation status of over 70,000 single nuclei (sn) from adult mouse pituitaries. Paired snRNAseq and snATACseq datasets from individual animals highlight a continuum between developmental epigenetically-encoded cell types and transcriptionally-determined transient cell states. Co-accessibility analysis-based identification of a putative Fshb cis-regulatory domain that overlaps the fertility-linked rs11031006 human polymorphism, followed by experimental validation illustrate the use of this resource for hypothesis generation. We also identify transcriptional and chromatin accessibility programs distinguishing each major cell type. Regulons, which are co-regulated gene sets sharing binding sites for a common transcription factor driver, recapitulate cell type clustering. We identify both cell type-specific and sex-specific regulons that are highly correlated with promoter accessibility, but not with methylation state, supporting the centrality of chromatin accessibility in shaping cell-defining transcriptional programs. The sn multi-omics atlas is accessible at snpituitaryatlas.princeton.edu.