Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity

Single-Cell RNA Sequencing and Its Applications in Pituitary Research

BACKGROUND

Mounting evidence underscores the significance of cellular diversity within the endocrine system and the intricate interplay between different cell types and tissues, essential for preserving physiological balance and influencing disease trajectories. The pituitary gland, a central player in the endocrine orchestra, exemplifies this complexity with its assortment of hormone-secreting and nonsecreting cells.

SUMMARY

The pituitary gland houses several types of cells responsible for hormone production, alongside nonsecretory cells like fibroblasts and endothelial cells, each playing a crucial role in the gland's function and regulatory mechanisms. Despite the acknowledged importance of these cellular interactions, the detailed mechanisms by which they contribute to pituitary gland physiology and pathology remain largely uncharted. The last decade has seen the emergence of groundbreaking technologies such as single-cell RNA sequencing, offering unprecedented insights into cellular heterogeneity and interactions. However, the application of this advanced tool in exploring the pituitary gland's complexities has been scant. This review provides an overview of this methodology, highlighting its strengths and limitations, and discusses future possibilities for employing it to deepen our understanding of the pituitary gland and its dysfunction in disease states.

KEY MESSAGE

Single-cell RNA sequencing technology offers an unprecedented means to study the heterogeneity and interactions of pituitary cells, though its application has been limited thus far. Further utilization of this tool will help uncover the complex physiological and pathological mechanisms of the pituitary, advancing research and treatment of pituitary diseases.

Transcriptomic Profiling of Lactotroph Pituitary Neuroendocrine Tumors via RNA Sequencing and Ingenuity Pathway Analysis

INTRODUCTION

Lactotroph pituitary neuroendocrine tumors (PitNETs) are common pituitary tumors, but their underlying molecular mechanisms remain unclear. This study aimed to investigate the transcriptomic landscape of lactotroph PitNETs and identify potential molecular mechanisms and therapeutic targets through RNA sequencing and ingenuity pathway analysis (IPA).

METHODS

Lactotroph PitNET tissues from five surgical cases without dopamine agonist treatment underwent RNA sequencing. Normal pituitary tissues from 3 patients served as controls. Differentially expressed genes (DEGs) were identified, and the functional pathways and gene networks were explored by IPA.

RESULTS

Transcriptome analysis revealed that lactotroph PitNETs had gene expression patterns that were distinct from normal pituitary tissues. We identified 1,172 upregulated DEGs, including nine long intergenic noncoding RNAs (lincRNAs) belonging to the top 30 DEGs. IPA of the upregulated DEGs showed that the estrogen receptor signaling, oxidative phosphorylation signaling, and EIF signaling were activated. In gene network analysis, key upstream regulators, such as EGR1, PRKACA, PITX2, CREB1, and JUND, may play critical roles in lactotroph PitNETs.

CONCLUSION

This study provides a comprehensive transcriptomic profile of lactotroph PitNETs and highlights the potential involvement of lincRNAs and specific signaling pathways in tumor pathogenesis. The identified upstream regulators may be potential therapeutic targets for future investigations.

Sex-Linked Growth Disorder and Aberrant Pituitary Gene Expression in Nestin-Cre-Mediated Egr1 Conditional Knockout Mice

Genes that regulate hormone release are essential for maintaining metabolism and energy balance. Egr1 encodes a transcription factor that regulates hormone production and release, and a decreased in growth hormones has been reported in Egr1 knockout mice. A reduction in growth hormones has also been observed in Nestin-Cre mice, a model frequently used to study the nervous system. Currently, it is unknown how Egr1 loss or the Nestin-Cre driver disrupt pituitary gene expression. Here, we compared the growth curves and pituitary gene expression profiles of Nestin-Cre-mediated Egr1 conditional knockout (Egr1cKO) mice with those of their controls. Reduced body weight was observed in both the Nestin-Cre and Egr1cKO mice, and the loss of Egr1 had a slightly more severe impact on female mice than on male mice. RNA-seq data analyses revealed that the sex-related differences were amplified in the Nestin-Cre-mediated Egr1 conditional knockout mice. Additionally, in the male mice, the influence of Egr1cKO on pituitary gene expression may be overridden by the Nestin-Cre driver. Differentially expressed genes associated with the Nestin-Cre driver were significantly enriched for genes related to growth factor activity and binding. Altogether, our results demonstrate that Nestin-Cre and the loss of Egr1 in the neuronal cell lineage have distinct impacts on pituitary gene expression in a sex-specific manner.

Practical Considerations for Complex Tissue Dissociation for Single-Cell Transcriptomics

Single-cell and single-nucleus RNA sequencing have revolutionized biomedical research, allowing analysis of complex tissues, identification of novel cell types, and mapping of development as well as disease states. Successful application of this technology critically relies on the dissociation of solid organs and tissues into high-quality single-cell (or nuclei) suspensions.In this chapter, we examine several key aspects of the tissue handling workflow that need to be considered when establishing an efficient tissue processing protocol for single-cell RNA sequencing (scRNA-seq). These include tissue collection, transport, and storage, as well as the choice of the dissociation conditions. We emphasize the importance of the tissue quality check and discuss the advantages (and potential limitations) of tissue cryopreservation. We provide practical tips and considerations on each of the steps of the processing workflow, and comment on how to maximize cell viability and integrity, which are critical for obtaining high-quality single-cell transcriptomic data.

A Simple Evolutionary Model of Genetic Robustness After Gene Duplication

When a dispensable gene is duplicated (referred to the ancestral dispensability denoted by O⁺), genetic buffering and duplicate compensation together maintain the duplicate redundancy, whereas duplicate compensation is the only mechanism when an essential gene is duplicated (referred to the ancestral essentiality denoted by O^-). To investigate these evolutionary scenarios of genetic robustness, I formulated a simple mixture model for analyzing duplicate pairs with one of the following states: double dispensable (DD), semi-dispensable (one dispensable one essential, DE), or double essential (EE). This model was applied to the yeast duplicate pairs from a whole-genome duplication (WGD) occurred about 100 million years ago (mya), and the mouse duplicate pairs from a WGD occurred about more than 500 mya. Both case studies revealed that the proportion of essentiality for those duplicates with ancestral essentiality [P_E(O^-)] was much higher than that for those with ancestral dispensability [P_E(O⁺)]. While it was negligible in the yeast duplicate pairs, P_E(O⁺) (about 20%) was shown statistically significant in the mouse duplicate pairs. These findings, together, support the hypothesis that both sub-functionalization and neo-functionalization may play some roles after gene duplication, though the former may be much faster than the later.

Single nucleus multi-omics regulatory landscape of the murine pituitary

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.

Cryopreservation of microglia enables single-cell RNA sequencing with minimal effects on disease-related gene expression patterns

Microglia play a key role in brain development, normal homeostasis, and neurodegenerative disorders. Single-cell technologies have led to important findings about microglia, with many animal model studies using single-cell RNA sequencing (scRNA-seq), whereas most human specimen studies using archived frozen brains for single-nucleus RNA sequencing (snRNA-seq). However, microglia compose a small proportion of the total brain tissue; snRNAseq depletes expression of microglia activation genes that characterize many diseases. Here we examine the use of purified, cryopreserved microglia for scRNA-seq. Comparison of scRNA-seq on paired fresh and cryopreserved microglia from rhesus monkeys revealed a high level of correlation of gene expression between the two conditions. Disease-related genes were relatively unaffected, but an increase in immediate-early gene expression was present in cryopreserved cells. Regardless, changes in immediate-early gene expression are still detectable. Cryopreservation of microglia is a suitable procedure for prospectively archiving samples.

Versatile workflow for cell type-resolved transcriptional and epigenetic profiles from cryopreserved human lung

Complexity of lung microenvironment and changes in cellular composition during disease make it exceptionally hard to understand molecular mechanisms driving development of chronic lung diseases. Although recent advances in cell type-resolved approaches hold great promise for studying complex diseases, their implementation relies on local access to fresh tissue, as traditional tissue storage methods do not allow viable cell isolation. To overcome these hurdles, we developed a versatile workflow that allows storage of lung tissue with high viability, permits thorough sample quality check before cell isolation, and befits sequencing-based profiling. We demonstrate that cryopreservation enables isolation of multiple cell types from both healthy and diseased lungs. Basal cells from cryopreserved airways retain their differentiation ability, indicating that cellular identity is not altered by cryopreservation. Importantly, using RNA sequencing and EPIC Array, we show that gene expression and DNA methylation signatures are preserved upon cryopreservation, emphasizing the suitability of our workflow for omics profiling of lung cells. Moreover, we obtained high-quality single-cell RNA-sequencing data of cells from cryopreserved human lungs, demonstrating that cryopreservation empowers single-cell approaches. Overall, thanks to its simplicity, our workflow is well suited for prospective tissue collection by academic collaborators and biobanks, opening worldwide access to viable human tissue.

Cell Atlas technologies and insights into tissue architecture

Since Robert Hooke first described the existence of 'cells' in 1665, scientists have sought to identify and further characterise these fundamental units of life. While our understanding of cell location, morphology and function has expanded greatly; our understanding of cell types and states at the molecular level, and how these function within tissue architecture, is still limited. A greater understanding of our cells could revolutionise basic biology and medicine. Atlasing initiatives like the Human Cell Atlas aim to identify all cell types at the molecular level, including their physical locations, and to make this reference data openly available to the scientific community. This is made possible by a recent technology revolution: both in single-cell molecular profiling, particularly single-cell RNA sequencing, and in spatially resolved methods for assessing gene and protein expression. Here, we review available and upcoming atlasing technologies, the biological insights gained to date and the promise of this field for the future.

Imaging of Endoplasmic Reticulum Ca2+ in the Intact Pituitary Gland of Transgenic Mice Expressing a Low Affinity Ca2+ Indicator

The adenohypophysis contains five secretory cell types (somatotrophs, lactotrophs, thyrotrophs, corticotrophs, and gonadotrophs), each secreting a different hormone, and controlled by different hypothalamic releasing hormones (HRHs). Exocytic secretion is regulated by cytosolic Ca²⁺ signals ([Ca²⁺]_C), which can be generated either by Ca²⁺ entry through the plasma membrane and/or by Ca²⁺ release from the endoplasmic reticulum (ER). In addition, Ca²⁺ entry signals can eventually be amplified by ER release via calcium-induced calcium release (CICR). We have investigated the contribution of ER Ca²⁺ release to the action of physiological agonists in pituitary gland. Changes of [Ca²⁺] in the ER ([Ca²⁺]_ER) were measured with the genetically encoded low-affinity Ca²⁺ sensor GAP3 targeted to the ER. We used a transgenic mouse strain that expressed erGAP3 driven by a ubiquitous promoter. Virtually all the pituitary cells were positive for the sensor. In order to mimick the physiological environment, intact pituitary glands or acute slices from the transgenic mouse were used to image [Ca²⁺]_ER. [Ca²⁺]_C was measured simultaneously with Rhod-2. Luteinizing hormone-releasing hormone (LHRH) or thyrotropin releasing hormone (TRH), two agonists known to elicit intracellular Ca²⁺ mobilization, provoked robust decreases of [Ca²⁺]_ER and concomitant rises of [Ca²⁺]_C. A smaller fraction of cells responded to thyrotropin releasing hormone (TRH). By contrast, depolarization with high K⁺ triggered a rise of [Ca²⁺]_C without a decrease of [Ca²⁺]_ER, indicating that the calcium-induced calcium-release (CICR) via ryanodine receptor amplification mechanism is not present in these cells. Our results show the potential of transgenic ER Ca²⁺ indicators as novel tools to explore intraorganellar Ca²⁺ dynamics in pituitary gland in situ.

scRNA-seq assessment of the human lung, spleen, and esophagus tissue stability after cold preservation

BACKGROUND

The Human Cell Atlas is a large international collaborative effort to map all cell types of the human body. Single-cell RNA sequencing can generate high-quality data for the delivery of such an atlas. However, delays between fresh sample collection and processing may lead to poor data and difficulties in experimental design.

RESULTS

This study assesses the effect of cold storage on fresh healthy spleen, esophagus, and lung from ≥ 5 donors over 72 h. We collect 240,000 high-quality single-cell transcriptomes with detailed cell type annotations and whole genome sequences of donors, enabling future eQTL studies. Our data provide a valuable resource for the study of these 3 organs and will allow cross-organ comparison of cell types. We see little effect of cold ischemic time on cell yield, total number of reads per cell, and other quality control metrics in any of the tissues within the first 24 h. However, we observe a decrease in the proportions of lung T cells at 72 h, higher percentage of mitochondrial reads, and increased contamination by background ambient RNA reads in the 72-h samples in the spleen, which is cell type specific.

CONCLUSIONS

In conclusion, we present robust protocols for tissue preservation for up to 24 h prior to scRNA-seq analysis. This greatly facilitates the logistics of sample collection for Human Cell Atlas or clinical studies since it increases the time frames for sample processing.

Using single-cell transcriptomics to understand functional states and interactions in microbial eukaryotes

Understanding the diversity and evolution of eukaryotic microorganisms remains one of the major challenges of modern biology. In recent years, we have advanced in the discovery and phylogenetic placement of new eukaryotic species and lineages, which in turn completely transformed our view on the eukaryotic tree of life. But we remain ignorant of the life cycles, physiology and cellular states of most of these microbial eukaryotes, as well as of their interactions with other organisms. Here, we discuss how high-throughput genome-wide gene expression analysis of eukaryotic single cells can shed light on protist biology. First, we review different single-cell transcriptomics methodologies with particular focus on microbial eukaryote applications. Then, we discuss single-cell gene expression analysis of protists in culture and what can be learnt from these approaches. Finally, we envision the application of single-cell transcriptomics to protist communities to interrogate not only community components, but also the gene expression signatures of distinct cellular and physiological states, as well as the transcriptional dynamics of interspecific interactions. Overall, we argue that single-cell transcriptomics can significantly contribute to our understanding of the biology of microbial eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.

Cell Type- and Sex-Dependent Transcriptome Profiles of Rat Anterior Pituitary Cells

Understanding the physiology and pathology of an organ composed of a variety of cell populations depends critically on genome-wide information on each cell type. Here, we report single-cell transcriptome profiling of over 6,800 freshly dispersed anterior pituitary cells from postpubertal male and female rats. Six pituitary-specific cell types were identified based on known marker genes and characterized: folliculostellate cells and hormone-producing corticotrophs, gonadotrophs, thyrotrophs, somatotrophs, and lactotrophs. Also identified were endothelial and blood cells from the pituitary capillary network. The expression of numerous developmental and neuroendocrine marker genes in both folliculostellate and hormone-producing cells supports that they have a common origin. For several genes, the validity of transcriptome analysis was confirmed by qRT-PCR and single cell immunocytochemistry. Folliculostellate cells exhibit impressive transcriptome diversity, indicating their major roles in production of endogenous ligands and detoxification enzymes, and organization of extracellular matrix. Transcriptome profiles of hormone-producing cells also indicate contributions toward those functions, while also clearly demonstrating their endocrine function. This survey highlights many novel genetic markers contributing to pituitary cell type identity, sexual dimorphism, and function, and points to relationships between hormone-producing and folliculostellate cells.

Follicle-Stimulating Hormone Glycobiology

FSH glycosylation varies in two functionally important aspects: microheterogeneity, resulting from oligosaccharide structure variation, and macroheterogeneity, arising from partial FSHβ subunit glycosylation. Although advances in mass spectrometry permit extensive characterization of FSH glycan populations, microheterogeneity remains difficult to illustrate, and comparisons between different studies are challenging because no standard format exists for rendering oligosaccharide structures. FSH microheterogeneity is illustrated using a consistent glycan diagram format to illustrate the large array of structures associated with one hormone. This is extended to commercially available recombinant FSH preparations, which exhibit greatly reduced microheterogeneity at three of four glycosylation sites. Macroheterogeneity is demonstrated by electrophoretic mobility shifts due to the absence of FSHβ glycans that can be assessed by Western blotting of immunopurified FSH. Initially, macroheterogeneity was hoped to matter more than microheterogeneity. However, it now appears that both forms of carbohydrate heterogeneity have to be taken into consideration. FSH glycosylation can reduce its apparent affinity for its cognate receptor by delaying initial interaction with the receptor and limiting access to all of the available binding sites. This is followed by impaired cellular signaling responses that may be related to reduced receptor occupancy or biased signaling. To resolve these alternatives, well-characterized FSH glycoform preparations are necessary.

Cytogenetic, Genomic, and Functional Characterization of Pituitary Gonadotrope Cell Lines

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.