Gene misexpression in a Smoc2+ve/ Sox2-low population in juvenile Prop1-mutant pituitary gland
Mutations in the pituitary specific transcription factor Prophet of Pit-1 (PROP1) are the most common genetic etiology of combined pituitary hormone deficiency (CPHD). CPHD is associated with short stature, attributable to growth hormone deficiency and/or thyroid stimulating hormone deficiency, as well as hypothyroidism and infertility. Pathogenic lesions impair pituitary development and differentiation of endocrine cells. We performed single-cell RNA sequencing of pituitary cells from a wild-type and a Prop1-mutant P4 female to elucidate population-specific differential gene expression. We observed a Smoc2+ve population that expressed low Sox2, which trajectory analyses suggest are a transitional cell state as stem cells differentiate into endocrine cells. We also detected ectopic expression of Sox21 in these cells in the Prop1df/df mutant. Prop1-mutant mice are known to overexpress Pou3f4, which we now show to be also enriched in this Smoc2+ve population. We sought to elucidate the role of Pou3f4 during pituitary development and to determine the contributions of Pou3f4 upregulation to pituitary disease by utilizing double-mutant mice lacking both Prop1 and Pou3f4. However, our data showed that Pou3f4 is not required for normal pituitary development and function. Double mutants further demonstrated that the upregulation of Pou3f4 was not causative for the overexpression of Sox21. These data indicate loss of Pou3f4 is not a potential cause of CPHD, and further studies may investigate the functional consequence of upregulation of Pou3f4 and Sox21, if any, in the novel Smoc2+ve cell population. Overall design: ScRNAseq RNA libraries were prepared from one P4 female Prop1+/+ pituitary gland and one P4 female Prop1df/df pituitary gland from the same litter. Single-cell RNA libraries were prepared by the University of Michigan's Advanced Genomics Core on the 10X Genomics Chromium 3' Gene Expression v3 platform following the manufacturer's instructions and sequenced on the NovaSeq6000 platform with an S4 chip. Samples were analyzed similar to our previous study 16. Samples were demultiplexed and fastq files aligned using Cellranger Single Cell Software Suite 3.1.0. Sequence reads were aligned to a modified mm10 containing extended annotation for the Prop1 3' untranslated region 17. Binary sequence alignment maps were converted into loom files using velocyto.py 18 and imported into R for primary analyses using Seurat 3.1.2 and 5.1.0 19,20. Core Seurat functions were used to import loom files, integrate samples, generate uniform manifold approximation and projection (UMAP), calculate clusters, and calculate differential gene expression. Expression plots used Nebulosa 1.12.1 21. Plots and images generated in R or Excel were compiled in Adobe Photoshop 2024.