Single-cell analysis reveals transcriptional dynamics in healthy primary parathyroid tissue

Studies on human parathyroids are generally limited to hyperfunctioning glands owing to the difficulty in obtaining normal human tissue. We therefore obtained non-human primate (NHP) parathyroids to provide a suitable alternative for sequencing that would bear a close semblance to human organs. Single-cell RNA expression analysis of parathyroids from four healthy adult M. mulatta reveals a continuous trajectory of epithelial cell states. Pseudotime analysis based on transcriptomic signatures suggests a progression from GCM2 ^hi progenitors to mature parathyroid hormone (PTH)-expressing epithelial cells with increasing core mitochondrial transcript abundance along pseudotime. We sequenced, as a comparator, four histologically characterized hyperfunctioning human parathyroids with varying oxyphil and chief cell abundance and leveraged advanced computational techniques to highlight similarities and differences from non-human primate parathyroid expression dynamics. Predicted cell-cell communication analysis reveals abundant endothelial cell interactions in the parathyroid cell microenvironment in both human and NHP parathyroid glands. We show abundant RARRES2 transcripts in both human adenoma and normal primate parathyroid cells and use coimmunostaining to reveal high levels of RARRES2 protein (also known as chemerin) in PTH-expressing cells, which could indicate that RARRES2 plays an unrecognized role in parathyroid endocrine function. The data obtained are the first single-cell RNA transcriptome to characterize nondiseased parathyroid cell signatures and to show a transcriptomic progression of cell states within normal parathyroid glands, which can be used to better understand parathyroid cell biology.