Single-cell transcriptomics identifies divergent developmental lineage trajectories during human pituitary development

A complex multisystem disorder including hypopituitarism and hypoparathyroidism, associated with mutation in the gene encoding fatty acid synthase (FASN)

Whole exome sequencing performed on a male patient with a complex phenotype including short stature associated with hypopituitarism, sensorineural deafness, hypoparathyroidism, retinal dystrophy, and developmental delay revealed a novel de novo variant in FASN (p.Ala2132Val), encoding fatty acid synthase. The patient failed to respond to growth-promoting treatment, only reaching a height of 128.3 cm (-6.98 SDS) at 24.7 years of age, and was prepubertal with a delayed bone age (13.6 years). Subsequent metabolic investigations demonstrated high triglyceride concentrations throughout an 18 h fast with a failure to increase 3-hydroxybutyrate, suggesting a defect in fatty acid oxidation or ketone body synthesis. Human embryonic brain analysis revealed FASN expression in the diencephalon, hypothalamus and Rathke's pouch. Following the labelling of glucose with carbon-13 (C13) in cultured fibroblasts, mass spectrometry data revealed that more C13-glucose was incorporated into de novo synthesised palmitic acid in controls compared to patient cells, suggesting reduced fatty acid synthesis in the patient. Our data suggest that the FASN p.Ala2132Val variant is associated with a complex phenotype including hypothalamo-pituitary dysfunction, consistent with previous studies showing that rodent neural/progenitor brain stem cells are governed by Fasn-dependent de novo lipogenesis (fatty acid synthesis) for proliferation.

Gonadotrophs have a dual origin, with most derived from early postnatal pituitary stem cells

Gonadotrophs are the essential pituitary endocrine cells for reproduction. They produce both luteinizing (LH) and follicle-stimulating (FSH) hormones that act on the gonads to promote germ cell maturation and steroidogenesis. Their secretion is controlled by the hypothalamic gonadotrophin-releasing hormone (GnRH), and gonadal steroid feedback. Gonadotrophs first appear in the embryonic pituitary, along with other endocrine cell types, and all expand after birth. While gonadotrophs may display heterogeneity in their response to GnRH, they appear, at least transcriptionally, as a homogenous population. The pituitary also contains a population of stem cells (SCs), whose contribution to postnatal growth is unclear, in part because endocrine cells maintain the ability to proliferate. Here we show an unsuspected dual origin of the murine adult gonadotroph population, with most gonadotrophs originating from postnatal pituitary stem cells starting early postnatally and up to puberty, while embryonic gonadotrophs are maintained. We further demonstrate that postnatal gonadotroph differentiation happens independently of gonadal signals and is not affected by impairment of GnRH signalling. The division of gonadotrophs based on separate origins has implications for our understanding of the establishment and regulation of reproductive functions, both in health and in disease.

Single-cell technology for cell-based drug delivery and pharmaceutical research

Leveraging the capacity to precisely manipulate and analyze individual cells, single-cell technology has rapidly become an indispensable tool in the advancement of cell-based drug delivery systems and innovative cell therapies. This technology offers powerful means to address cellular heterogeneity and significantly enhance therapeutic efficacy. Recent breakthroughs in techniques such as single-cell electroporation, mechanical perforation, and encapsulation, particularly when integrated with microfluidics and bioelectronics, have led to remarkable improvements in drug delivery efficiency, reductions in cytotoxicity, and more precise targeting of therapeutic effects. Moreover, single-cell analyses, including advanced sequencing and high-resolution sensing, offer profound insights into complex disease mechanisms, the development of drug resistance, and the intricate processes of stem cell differentiation. This review summarizes the most significant applications of these single-cell technologies, highlighting their impact on the landscape of modern biomedicine. Furthermore, it provides a forward-looking perspective on future research directions aimed at further optimizing drug delivery strategies and enhancing therapeutic outcomes in the treatment of various diseases.

Transcriptomic sex differences in early human fetal brain development

The influence of sex chromosomes and sex hormones on early human brain development is poorly understood. We therefore undertook transcriptomic analysis of 46,XY and 46,XX human brain cortex samples (n = 64) at four different time points between 7.5 and 17 weeks post conception (wpc), in two independent studies. This developmental period encompasses the onset of testicular testosterone secretion in the 46,XY fetus (8wpc). We show differences in sex chromosome gene expression including X-inactivation genes (XIST, TSIX) in 46,XX samples; core Y chromosome genes (n = 18) in 46,XY samples; and two Y chromosome brain specific genes, PCDH11Y and RP11-424G14.1. PCDH11Y (protocadherin11 Y-linked) regulates excitatory neurons; this gene is unique to humans and is implicated in language development. RP11-424G14.1 is a long non-coding RNA. Fewer differences in sex hormone pathway-related genes are seen. The androgen receptor (AR, NR3C4) shows cortex expression in both sexes, which decreases with age. Global cortical sex hormone effects are not seen, but more localized AR mechanisms may be important with time (e.g., hypothalamus). Taken together, our data suggest that limited but potentially important sex differences occur during early human fetal brain development.

Old players and new insights: unraveling the role of RNA-binding proteins in brain tumors

The human genome harbors >1,600 evolutionarily conserved RNA-binding proteins (RBPs), with extensive multi-omics investigations documenting their pervasive dysregulation in malignancies ranging from glioblastoma to melanoma. These RBPs are integral to the complex regulatory networks governing hallmark cancer processes. Recent studies have investigated the multifaceted contributions of RBPs to tumorigenesis, tumor metabolism, the tumor-immune microenvironment, and resistance to therapy. This complexity is further compounded by the intricate regulation of RNA function at various levels by RBPs, as well as the post-translational modifications of RBPs, which improve their functional capacity. Moreover, numerous RBP-based therapeutics have emerged, each underpinned by distinct molecular mechanisms that extend from genomic analysis to the interference of RBPs' function. This review aims to provide a comprehensive overview of the recent progress in the meticulous roles of RBPs in brain tumors and to explore potential therapeutic interventions targeting these RBPs, complemented by a discussion of innovative techniques emerging in this research field. Advances in deciphering RNA-RBP interactomes and refining targeted therapeutic strategies are revealing the transformative potential of RBP-centric approaches in brain tumor treatment, establishing them as pivotal agents for overcoming current clinical challenges.

Human Pituitary Organoids: Transcriptional Landscape Deciphered by scRNA-Seq and Stereo-Seq, with Insights into SOX3's Role in Pituitary Development

The 3D human pituitary organoid represents a promising laboratory model for investigating human pituitary diseases. Nonetheless, this technology is still in its nascent stage, with uncertainties regarding the cellular composition, intercellular interactions, and spatial distribution of the human pituitary organoids. To address these gaps, the culture conditions are systematically adjusted and the efficiency of induced pluripotent stem cells' (iPSCs') differentiation into pituitary organoids is successfully improved, achieving results comparable to or exceeding those of previous studies. Additionally, single-cell RNA-sequencing (scRNA-seq) and stereomics sequencing (Stereo-seq) are performed on the pituitary organoids for the first time, and unveil the diverse cell clusters, intricate intercellular interactions, and spatial information within the organoids. Furthermore, the SOX3 gene interference impedes the iPSCs' differentiation into pituitary organoids, thereby highlighting the potential of pituitary organoids as an ideal experimental model. Altogether, the research provides an optimized protocol for the human pituitary organoid culture and a valuable transcriptomic dataset for future explorations, laying the foundation for subsequent research in the field of pituitary organoids or pituitary diseases.

Expanding the Spectrum of Endocrine Abnormalities Associated With SOX11-related Disorders

CONTEXT

SOX11 variants cause Coffin-Siris syndrome, characterized by developmental delay, hypogonadotropic hypogonadism, and skeletal and facial defects.

OBJECTIVE

To examine the contribution of SOX11 variants to the pathogenesis of idiopathic hypogonadotropic hypogonadism (IHH), a disorder caused by hypothalamic GnRH deficiency.

SETTING

The Reproductive Endocrine Unit and the Pediatric Endocrinology Division, Massachusetts General Hospital.

PATIENTS OR OTHER PARTICIPANTS

A cohort of 1810 unrelated IHH probands.

INTERVENTIONS

Exome sequencing data from the entire cohort were examined for SOX11 rare single nucleotide variants (SNVs) (minor allele frequency in the gnomAD database <0.1%). Rare SOX11 variant association testing was performed between the IHH and gnomAD population. Phenotyping of individuals harboring pathogenic/likely pathogenic SNVs (determined by the American College of Medical Genetics criteria) was performed.

MAIN OUTCOMES/RESULTS

Four pathogenic SOX11 SNVs were identified in 5 IHH probands. The IHH cohort was enriched for SOX11 protein truncating SNVs (frameshift/nonsense) across the entire protein (2 SNVs in 3 IHH cases [p.S303X (de novo); p.S345Afs*13]; P = .0004981) and for SOX11 missense SNVs within the SOX11 high-mobility group domain (2 SNVs in 2 IHH cases p.G84D [de novo]; p.P114S; P = .00313922). The phenotypic spectrum of SOX11 variant carriers revealed additional endocrine defects including anosmic and normosmic forms of IHH, GH deficiency, pituitary and hypothalamic structural defects, and hypothyroidism. A pathogenic SOX11 SNV was also identified in a patient with functional hypogonadotropic hypogonadism (p.R100Q). Coffin-Siris syndrome-associated features were present in 4/5 probands.

CONCLUSION

Deleterious SOX11 variants cause IHH and other pituitary hormone deficiencies, suggesting that the human SOX11-associated disorder may stem from both hypothalamic and pituitary level defects.

Transcriptome Complexity Disentangled: A Regulatory Molecules Approach

Transcription factors (TFs) and microRNAs (miRNAs) are fundamental regulators of gene expression, cell state, and biological processes. This study investigated whether a small subset of TFs and miRNAs could accurately predict genome-wide gene expression. We analyzed 8895 samples across 31 cancer types from The Cancer Genome Atlas and identified 28 miRNA and 28 TF clusters using unsupervised learning. Medoids of these clusters could differentiate tissues of origin with 92.8% accuracy, demonstrating their biological relevance. We developed Tissue-Agnostic and Tissue-Aware models to predict 20,000 gene expressions using the 56 selected medoid miRNAs and TFs. The Tissue-Aware model attained an R2 of 0.70 by incorporating tissue-specific information. Despite measuring only 1/400th of the transcriptome, the prediction accuracy was comparable to that achieved by the 1000 landmark genes. This suggests the transcriptome has an intrinsically low-dimensional structure that can be captured by a few regulatory molecules. Our approach could enable cheaper transcriptome assays and analysis of low-quality samples. It also provides insights into genes that are heavily regulated by miRNAs/TFs versus alternative mechanisms. However, model transportability was impacted by dataset discrepancies, especially in miRNA distribution. Overall, this study demonstrates the potential of a biology-guided approach for robust transcriptome representation.

Immune-developmental processes contribute to schizophrenia risk: insights from a genetic overlap study with height

BACKGROUND

Shorter stature has been phenotypically linked to increased prevalence of schizophrenia (SCZ), but the nature of this association is unknown.

METHODS

Using genome-wide genetic data, we studied the SCZ-height relationship on a genetic level. Applying novel genetic methods and tools, we analyzed gene-sets, tissue-types, cell-types, local genetic correlation, conditional genetic analyses, and fine-mapping of effector-genes to scrutinize the SCZ-height relationship.

RESULTS

We identified 142 genes statistically associated with both SCZ and height and found enrichment in 3 functional gene-sets. Genetic annotations implicated the pituitary and specifically mesenchymal stem cells for height and thyrotropic cells for SCZ. While the global SCZ-height genetic correlation was nonsignificant, 9 genomic regions showed robust local genetic correlations (7 negative, 6 in the MHC-region). The shared genetic signal for SCZ and height within the 6 MHC-regions was partially explained by mutual genetic overlap with white blood cell count, particularly lymphocytes. Fine-mapping prioritized 3 shared effector-genes (GIGYF2, HLA-C, and LIN28B) involved in immune response sensitivity and development of immune and pituitary cell-types.

CONCLUSIONS

Overall, our findings suggest an involvement during height-development of thyrotropic cells and immune response sensitivity contributing towards risk of SCZ.

CoupleVAE: coupled variational autoencoders for predicting perturbational single-cell RNA sequencing data

With the rapid advances in single-cell sequencing technology, it is now feasible to conduct in-depth genetic analysis in individual cells. Study on the dynamics of single cells in response to perturbations is of great significance for understanding the functions and behaviors of living organisms. However, the acquisition of post-perturbation cellular states via biological experiments is frequently cost-prohibitive. Predicting the single-cell perturbation responses poses a critical challenge in the field of computational biology. In this work, we propose a novel deep learning method called coupled variational autoencoders (CoupleVAE), devised to predict the postperturbation single-cell RNA-Seq data. CoupleVAE is composed of two coupled VAEs connected by a coupler, initially extracting latent features for controlled and perturbed cells via two encoders, subsequently engaging in mutual translation within the latent space through two nonlinear mappings via a coupler, and ultimately generating controlled and perturbed data by two separate decoders to process the encoded and translated features. CoupleVAE facilitates a more intricate state transformation of single cells within the latent space. Experiments in three real datasets on infection, stimulation and cross-species prediction show that CoupleVAE surpasses the existing comparative models in effectively predicting single-cell RNA-seq data for perturbed cells, achieving superior accuracy.

Single-cell transcriptome atlas of male mouse pituitary across postnatal life highlighting its stem cell landscape

The pituitary represents the master gland governing the endocrine system. We constructed a single-cell (sc) transcriptomic atlas of male mouse endocrine pituitary by incorporating existing and new data, spanning important postnatal ages in both healthy and injured condition. We demonstrate strong applicability of this new atlas to unravel pituitary (patho)biology by focusing on its stem cells and investigating their complex identity (unveiling stem cell markers) and niche (pinpointing regulatory factors). Importantly, we functionally validated transcriptomic findings using pituitary stem cell organoids, revealing roles for Krüppel-like transcription factor 5 (KLF5), activator protein-1 (AP-1) complex and epidermal growth factor (EGF) pathways in pituitary stem cell regulation. Our investigation substantiated changes in stem cell dynamics during aging, reinforcing the inflammatory/immune nature in elderly pituitary and stem cells. Finally, we show translatability of mouse atlas-based findings to humans, particularly regarding aging-associated profile. This pituitary sc map is a valuable tool to unravel pituitary (patho)biology.

Splicing diversity enhances the molecular classification of pituitary neuroendocrine tumors

Pituitary neuroendocrine tumors (PitNETs) are one of the most common intracranial tumors with diverse clinical manifestations. Current pathological classification systems rely primarily on histological hormone staining and transcription factors (TFs) expression. While effective in identifying three major lineages, molecular characteristics based on hormones and TFs lack sufficient resolution to fully capture the complex tumor heterogeneity. Transcriptional diversity by alternative splicing (AS) offered additional insight to address this challenge. Here, we perform bulk and full-length single-cell RNA sequencing to comprehensively investigate AS dysregulation across all PitNET lineages. We reveal pervasive splicing dysregulations that better depict tumor heterogeneity. Additionally, we delineate fundamental splicing heterogeneity at single-cell resolution, confirming bulk findings and refining splicing dysregulation varying among tumor cell types. Notably, we effectively distinguish the silent corticotroph subtype and define a distinct TPIT lineage subtype, which is associated with worse clinical outcomes and increased splicing abnormalities driven by altered ESRP1 expression. In conclusion, our results characterize the subtype specific AS landscape in PitNETs, enhancing the understanding of the PitNETs subtyping.

Update on regulation of GHRH and its actions on GH secretion in health and disease

This review focuses on our current understanding of how growth hormone releasing hormone (GHRH): 1) stimulates GH release and synthesis from pituitary growth hormone (GH)-producing cells (somatotropes), 2) drives somatotrope proliferation, 3) is negatively regulated by somatostatin (SST), GH and IGF1, 4) is altered throughout lifespan and in response to metabolic challenges, and 5) analogues can be used clinically to treat conditions of GH excess or deficiency. Although a large body of early work provides an underpinning for our current understanding of GHRH, this review specifically highlights more recent work that was made possible by state-of-the-art analytical tools, receptor-specific agonists and antagonists, high-resolution in vivo and ex vivo imaging and the development of tissue (cell) -specific ablation mouse models, to paint a more detailed picture of the regulation and actions of GHRH.

Co-profiling of single-cell gene expression and chromatin landscapes in stickleback pituitary

The pituitary gland is a key endocrine gland with various physiological functions including metabolism, growth, and reproduction. It comprises several distinct cell populations that release multiple polypeptide hormones. Although the major endocrine cell types are conserved across taxa, the regulatory mechanisms of gene expression and chromatin organization in specific cell types remain poorly understood. Here, we performed simultaneous profiling of the transcriptome and chromatin landscapes in the pituitary cells of the three-spined stickleback (Gasterosteus aculeatus), which represents a good model for investigating the genetic mechanisms underlying adaptive evolution. We obtained pairwise gene expression and chromatin profiles for 5184 cells under short- and long-day conditions. Using three independent clustering analyses, we identified 16 distinct cell clusters and validated their consistency. These results advance our understanding of the regulatory dynamics in the pituitary gland and provide a reference for future research on comparative physiology and evolutionary biology.

Comparative evaluation of ACetic - MEthanol high salt dissociation approach for single-cell transcriptomics of frozen human tissues

Current dissociation methods for solid tissues in scRNA-seq studies do not guarantee intact single-cell isolation, especially for sensitive and complex human endocrine tissues. Most studies rely on enzymatic dissociation of fresh samples or nuclei isolation from frozen samples. Dissociating whole intact cells from fresh-frozen samples, commonly collected by biobanks, remains a challenge. Here, we utilized the acetic-methanol dissociation approach (ACME) to capture transcriptional profiles of individual cells from fresh-frozen tissue samples. This method combines acetic acid-based dissociation and methanol-based fixation. In our study, we optimized this approach for human endocrine tissue samples for the first time. We incorporated a high-salt washing buffer instead of the standard PBS to stabilize RNA and prevent RNases reactivation during rehydration. We have designated this optimized protocol as ACME HS (ACetic acid-MEthanol High Salt). This technique aims to preserve cell morphology and RNA integrity, minimizing transcriptome changes and providing a more accurate representation of mature mRNA. We compared the ability of enzymatic, ACME HS, and nuclei isolation methods to preserve major cell types, gene expression, and standard quality parameters across 41 tissue samples. Our results demonstrated that ACME HS effectively dissociates and fixes cells, preserving cell morphology and high RNA integrity. This makes ACME HS a valuable alternative for scRNA-seq protocols involving challenging tissues where obtaining a live cell suspension is difficult or disruptive.

Modeling corticotroph deficiency with pituitary organoids supports the functional role of NFKB2 in human pituitary differentiation

Deficient Anterior pituitary with common Variable Immune Deficiency (DAVID) syndrome results from NFKB2 heterozygous mutations, causing adrenocorticotropic hormone deficiency (ACTHD) and primary hypogammaglobulinemia. While NFKB signaling plays a crucial role in the immune system, its connection to endocrine symptoms is unclear. We established a human disease model to investigate the role of NFKB2 in pituitary development by creating pituitary organoids from CRISPR/Cas9-edited human induced pluripotent stem cells (hiPSCs). Introducing homozygous TBX19K146R/K146R missense pathogenic variant in hiPSC, an allele found in congenital isolated ACTHD, led to a strong reduction of corticotrophs number in pituitary organoids. Then, we characterized the development of organoids harboring NFKB2D865G/D865G mutations found in DAVID patients. NFKB2D865G/D865G mutation acted at different levels of development with mutant organoids displaying changes in the expression of genes involved on pituitary progenitor generation (HESX1, PITX1, LHX3), hypothalamic secreted factors (BMP4, FGF8, FGF10), epithelial-to-mesenchymal transition, lineage precursors development (TBX19, POU1F1) and corticotrophs terminal differentiation (PCSK1, POMC), and showed drastic reduction in the number of corticotrophs. Our results provide strong evidence for the direct role of NFKB2 mutations in the endocrine phenotype observed in patients leading to a new classification of a NFKB2 variant of previously unknown clinical significance as pathogenic in pituitary development.

Recent Advances in the Understanding of Gonadotrope Lineage Differentiation in the Developing Pituitary

BACKGROUND

The pituitary gland is a vital endocrine organ regulating body homoeostasis through six hormone-secreting cell types. Among these, pituitary gonadotrope cells are essential for reproductive function. Throughout pituitary ontogenesis, gonadotrope cells differentiate in a stepwise process, involving both morphogenic cues and transcription factors, which drives specification of progenitor cells into specialised endocrine cells. It is crucial to understand the mechanisms underlying gonadotrope differentiation, as developmental defects and abnormalities in this process can lead to many reproductive pathologies.

SUMMARY

This review offers a detailed overview of the latest advances in gonadotrope cell differentiation. We addressed this question with a specific focus on three important aspects of gonadotrope differentiation: the identification of the progenitor population giving rise to gonadotrope cells, the early mechanisms that initiate Nr5a1 expression and thus gonadotrope fate commitment, and finally, the mechanisms driving the formation of physical and functional gonadotrope networks.

KEY MESSAGES

Overall, this review aimed to provide new insights into three aspects of the gonadotrope differentiation process by reconsidering pioneering studies in the light of data gained from latest technological developments. Firstly, we re-investigated the long debated developmental trajectory of pituitary gonadotrope cells. Secondly, we reported new regulatory mechanisms of Nr5a1 expression, focusing on the involvement of ERα. Finally, we highlighted the molecular and cellular mechanisms driving gonadotrope network formation during embryogenesis, a process that seems essential for regulation of gonadotrope activity.

BACKGROUND

The pituitary gland is a vital endocrine organ regulating body homoeostasis through six hormone-secreting cell types. Among these, pituitary gonadotrope cells are essential for reproductive function. Throughout pituitary ontogenesis, gonadotrope cells differentiate in a stepwise process, involving both morphogenic cues and transcription factors, which drives specification of progenitor cells into specialised endocrine cells. It is crucial to understand the mechanisms underlying gonadotrope differentiation, as developmental defects and abnormalities in this process can lead to many reproductive pathologies.

SUMMARY

This review offers a detailed overview of the latest advances in gonadotrope cell differentiation. We addressed this question with a specific focus on three important aspects of gonadotrope differentiation: the identification of the progenitor population giving rise to gonadotrope cells, the early mechanisms that initiate Nr5a1 expression and thus gonadotrope fate commitment, and finally, the mechanisms driving the formation of physical and functional gonadotrope networks.

KEY MESSAGES

Overall, this review aimed to provide new insights into three aspects of the gonadotrope differentiation process by reconsidering pioneering studies in the light of data gained from latest technological developments. Firstly, we re-investigated the long debated developmental trajectory of pituitary gonadotrope cells. Secondly, we reported new regulatory mechanisms of Nr5a1 expression, focusing on the involvement of ERα. Finally, we highlighted the molecular and cellular mechanisms driving gonadotrope network formation during embryogenesis, a process that seems essential for regulation of gonadotrope activity.

Fertility induction associated with improved peripheral reproductive parameters in male Prop1df/df mice subjected to GH and levothyroxine replacement

Objective

The aim of this study was to characterize the parameters of reproductive anatomy and pituitary hormone expression levels in ames dwarf mice (Prop1 df/df).

Materials and methods

Male Prop1df/df mice aged 30 days received daily intraperitoneal injections of recombinant human GH and levothyroxine three times weekly for 60 days. The sexual maturation of these animals was compared with that of their wild-type (Prop +/+) and untreated (Prop1df/df ) siblings.

Results

The Prop1df/df treated group developed sexual maturation 2 weeks later than the Prop +/+ group and presented an increase in testicular weight, complete spermatogenesis, and enhanced LH and FSH expression. The Prop1df/df untreated group had low FSH expression and no offspring; most animals in this group did not develop sexual maturation during the study period.

Conclusion

Replacement with GH and levothyroxine appeared to play a crucial role in restoring peripheral reproductive parameters and increasing pituitary hormone expression in Prop1df/df mice.

Joint representation and visualization of derailed cell states with Decipher

Biological insights often depend on comparing conditions such as disease and health, yet we lack effective computational tools for integrating single-cell genomics data across conditions or characterizing transitions from normal to deviant cell states. Here, we present Decipher, a deep generative model that characterizes derailed cell-state trajectories. Decipher jointly models and visualizes gene expression and cell state from normal and perturbed single-cell RNA-seq data, revealing shared and disrupted dynamics. We demonstrate its superior performance across diverse contexts, including in pancreatitis with oncogene mutation, acute myeloid leukemia, and gastric cancer.

Comprehensive mapping of somatotroph pituitary neuroendocrine tumour heterogeneity using spatial and single-cell transcriptomics

BACKGROUND

Pituitary neuroendocrine tumours (PitNETs) are common intracranial tumours that are highly heterogeneous with unpredictable growth patterns. The driver genes and mechanisms that are crucial for tumour progression in somatotroph PitNETs are poorly understood.

METHODS

In this study, we performed integrative spatial transcriptomics (ST) and single-cell RNA sequencing (scRNA-seq) analysis on somatotroph tumours and normal pituitary samples to comprehensively characterize the differences in cellular characteristics.

RESULTS

By analyzing combined copy number variations (CNVs), tumour tissues were divided into two regions, which included the CNVhigh and CNVlow areas. The protumour genes DLK1 and RCN1 were highly expressed in the CNVhigh area, which might be related to tumour progression and could be targeted for precision therapy. We also found that the transforming growth factor beta signalling pathway participated in tumour progression and identified heterogeneity in the expression profiles of key genes. We assessed the intertumoral and intratumoral heterogeneity in somatotroph PitNETs and emphasized the importance of individualized treatment.

CONCLUSION

In summary, we visualized the cellular distribution and transcriptional differences in normal pituitary and somatotroph PitNETs by ST and scRNA-seq for the first time. This study provides a strong theoretical foundation to comprehensively understand the crucial mechanisms involved in tumour progression and develop new strategies to treat somatotroph PitNETs.

KEY POINTS

The first-ever visualization of cellular distributions in normal and tumor pituitary tissues. The inter- and intra-tumoral transcriptomic heterogeneity of somatotroph PitNETs was comprehensively revealed. Identification of potential protumor factors and critical signaling pathways, opening new avenues for therapeutic intervention.

GnRH-Gonadotropes Interactions Revealed by Pituitary Single-cell Transcriptomics in Zebrafish

GnRH governs reproduction by regulating pituitary gonadotropins. Unlike most vertebrates, gnrh-/- zebrafish are fertile. To elucidate the role of the hypophysiotropic-Gnrh3 and other mechanisms regulating pituitary gonadotropes, we profiled the gene expression of all individual pituitary cells of wild-type and gnrh3-/- adult female zebrafish. The single-cell RNA sequencing showed that LH and FSH gonadotropes express the 2 gonadotropin beta subunits with a ratio of 140:1 (lhb:fshb) and 4:1 (fshb:lhb), respectively. Lh gonadotropes predominantly express genes encoding receptors for GnRH (gnrhr2), thyroid hormone, estrogen, and steroidogenic factor 1. No GnRH receptor transcript was enriched in FSH gonadotropes. Instead, cholecystokinin receptor-b and galanin receptor-1b transcripts were enriched in these cells. The loss of the Gnrh3 gene in gnrh3-/- zebrafish resulted in downregulation of fshb in LH gonadotropes and upregulation of pituitary hormones like TSH, GH, prolactin, and proopiomelanocortin-a. Likewise, targeted chemogenetic ablation of Gnrh3 neurons led to a decrease in the number of fshb+, lhb + and fshb+/lhb + cells. Our studies suggest that Gnrh3 directly acts on LH gonadotropes through Gnrhr2, but the outcome of this interaction is still unknown. Gnrh3 also regulates fshb expression in both gonadotropes, most likely via a non-GnRH receptor route. Altogether, while LH secretion and synthesis are likely regulated in a GnRH-independent manner, Gnrh3 seems to play a role in the cellular organization of the pituitary. Moreover, the coexpression of lhb and fshb in both gonadotropes provides a possible explanation as to why gnrh3-/- zebrafish are fertile.

Long Noncoding RNAs Expressed in Mouse Pituitary Development and Mature Hormone-Producing Cells

Mammalian genomes contain thousands of genes for long noncoding RNA (lncRNAs), some of which have been shown to affect protein coding gene expression through diverse mechanisms. The lncRNA transcripts are longer than 200 nucleotides and are often capped, spliced, and polyadenylated, but not translated into protein. Nuclear lncRNAs can modify chromatin structure and transcription in trans or cis by interacting with the DNA, forming R-loops, and recruiting regulatory proteins. Not much is known about the role of lncRNA in pituitary gland differentiation and function. We mined transcriptome data from mouse pituitary glands collected at embryonic days 12.5 and 14.5 and identified over 200 different lncRNA transcripts. To develop a research resource for the study of lncRNA, we used pituitary cre transgenes to tag pituitary cell types in adult mice with fluorescent markers, and enriched for thyrotropes, gonadotropes, and somatotropes using fluorescence-activated cell sorting. We determined the transcriptome of each cell population using RNA sequencing and mined the data for lncRNA. We detected hundreds of lncRNAs in adult pituitary cells; a few were located immediately nearby genes that encode pituitary hormones or lineage-specific transcription factors. The location of these lncRNAs suggests the possibility of a cis-acting regulatory role in pituitary development or function, and we observe coordinated expression of 2 of them with their putative target genes in transgenic mice. This research resource sets the foundation for examining the actions of lncRNAs on their putative target genes and determining whether they have roles during development and in response to physiological demand.

Morphological identification and distribution comparison of telocytes in pituitary gland between normal and cryptorchid yaks

BACKGROUND

Telocytes (TCs) is a novel type of interstitial cells in many mammals organs, which participate in the organizational metabolism, mechanical support, immunomodulation and other aspects. The aim of this study was to explore the organizational chemical characteristics of TCs in pituitary gland and their changes in cryptorchid yaks.

METHODS

Transmission electron microscopy (TEM), toluidine blue staining, immunofluorescence, qRT-PCR, and Western blotting may enable us to understand TCs distribution characteristics and biological functions.

RESULT

TEM confirmed the presence of TCs in the pituitary gland with small bodies and moniliform telopodes (Tps). The Tps extending out from the cell body to the peri-sinusoidal vessels spaces, the number of Tps is closely related to the morphology of the nucleus. The most obvious changes of TCs in the pituitary gland of cryptorchid yaks is the Tps are relatively shorter and decreased secretory vesicles. H.E. and toluidine blue staining revealed that TCs not only distributed between the sinusoidal blood vessels and the glandular cell clusters, but also present on the surface of vascular endothelial cells. The co-expression of TCs biomarkers, such as Vimentin/CD34, CD117/CD34 and α-SMA/CD34, were evaluated by immunofluorescence to further determine the phenotypic characteristics of TCs. Besides, we analyzed the mRNA and protein expression of these biomarkers to determine the characteristics of TCs changes and possible biological roles. Both the mRNA and protein expression of CD117 were significantly higher in the pituitary gland of cryptorchid yaks than in the normal (p < 0.01), the protein expression of CD34 in the cryptorchid yaks was significantly higher than the normal (p < 0.01). There were no significant difference in mRNA expression of Vimentin and α-SMA (p>0.05), while the protein expression were significantly increased in the normal yaks (p < 0.05).

CONCLUSIONS

In summary, this study reports for the first time that the biological characteristics of TCs in yak pituitary gland. Although there is no significant change in the distribution characteristics, the changes in biological features of TCs in cryptorchid yaks are clear, suggesting that TCs participated in alteration in the local microenvironment of the pituitary gland. Therefore, our study provides clues for further investigating the role of TCs in the pituitary gland during the occurrence of cryptorchidism in yaks.

The adipokine profile in the plasma and anterior pituitary of pigs during the estrous cycle

Adipokines play crucial roles in both reproductive and energy metabolic processes. This study aimed to compare the hormonal plasma profile of adiponectin, apelin, vaspin, chemerin, resistin, visfatin, and adipolin, and the expression of their receptors in the anterior pituitary (AP) between normal-weight Large White (LW) and fat Meishan (MS) pigs during different phases of the estrous cycle. We measured adipokine levels in the plasma and assessed their gene expression in the AP. We used Pearson's correlation analysis to examine potential links between adipokines levels, their receptors, and metabolic parameters (body weight; backfat thickness) and reproductive parameters (pituitary weight; age at puberty; levels of gonadotropins, steroid hormones; and gene expression of gonadotropin-releasing hormone receptor and gonadotropins in AP). The plasma levels of the evaluated adipokines fluctuated with phase and breed, except for visfatin and adipolin. Moreover, adipokine expression in AP varied significantly between breeds and estrous cycle phases, except for resistin receptor CAP1. Notably, we observed a positive correlation between plasma levels of adiponectin and its transcript in the AP only in MS pigs. Apelin gene expression correlated negatively with its receptor in MS, while we observed a breed-dependent correlation between chemerin gene expression and its receptor CMKLR1. We identified significant positive or negative correlations between adipokines or their receptor levels in plasma and AP as well as metabolic or reproductive parameters, depending on the breed. In conclusion, we have demonstrated breed-specific and estrous cycle-dependent regulation of adipokines in AP, underscoring their potential impact on metabolic and reproductive processes in swine.

Single nucleus/cell RNA-seq of the chicken hypothalamic-pituitary-ovarian axis offers new insights into the molecular regulatory mechanisms of ovarian development

The hypothalamic-pituitary-ovarian (HPO) axis represents a central neuroendocrine network essential for reproductive function. Despite its critical role, the intrinsic heterogeneity within the HPO axis across vertebrates and the complex intercellular interactions remain poorly defined. This study provides the first comprehensive, unbiased, cell type-specific molecular profiling of all three components of the HPO axis in adult Lohmann layers and Liangshan Yanying chickens. Within the hypothalamus, pituitary, and ovary, seven, 12, and 13 distinct cell types were identified, respectively. Results indicated that the pituitary adenylate cyclase activating polypeptide (PACAP), follicle-stimulating hormone (FSH), and prolactin (PRL) signaling pathways may modulate the synthesis and secretion of gonadotropin-releasing hormone (GnRH), FSH, and luteinizing hormone (LH) within the hypothalamus and pituitary. In the ovary, interactions between granulosa cells and oocytes involved the KIT, CD99, LIFR, FN1, and ANGPTL signaling pathways, which collectively regulate follicular maturation. The SEMA4 signaling pathway emerged as a critical mediator across all three tissues of the HPO axis. Additionally, gene expression analysis revealed that relaxin 3 (RLN3), gastrin-releasing peptide (GRP), and cocaine- and amphetamine regulated transcripts (CART, also known as CARTPT) may function as novel endocrine hormones, influencing the HPO axis through autocrine, paracrine, and endocrine pathways. Comparative analyses between Lohmann layers and Liangshan Yanying chickens demonstrated higher expression levels of GRP, RLN3, CARTPT, LHCGR, FSHR, and GRPR in the ovaries of Lohmann layers, potentially contributing to their superior reproductive performance. In conclusion, this study provides a detailed molecular characterization of the HPO axis, offering novel insights into the regulatory mechanisms underlying reproductive biology.

SMIM1 absence is associated with reduced energy expenditure and excess weight

BACKGROUND

Obesity rates have nearly tripled in the past 50 years, and by 2030 more than 1 billion individuals worldwide are projected to be obese. This creates a significant economic strain due to the associated non-communicable diseases. The root cause is an energy expenditure imbalance, owing to an interplay of lifestyle, environmental, and genetic factors. Obesity has a polygenic genetic architecture; however, single genetic variants with large effect size are etiological in a minority of cases. These variants allowed the discovery of novel genes and biology relevant to weight regulation and ultimately led to the development of novel specific treatments.

METHODS

We used a case-control approach to determine metabolic differences between individuals homozygous for a loss-of-function genetic variant in the small integral membrane protein 1 (SMIM1) and the general population, leveraging data from five cohorts. Metabolic characterization of SMIM1-/- individuals was performed using plasma biochemistry, calorimetric chamber, and DXA scan.

FINDINGS

We found that individuals homozygous for a loss-of-function genetic variant in SMIM1 gene, underlying the blood group Vel, display excess body weight, dyslipidemia, altered leptin to adiponectin ratio, increased liver enzymes, and lower thyroid hormone levels. This was accompanied by a reduction in resting energy expenditure.

CONCLUSION

This research identified a novel genetic predisposition to being overweight or obese. It highlights the need to investigate the genetic causes of obesity to select the most appropriate treatment given the large cost disparity between them.

FUNDING

This work was funded by the National Institute of Health Research, British Heart Foundation, and NHS Blood and Transplant.

Diatomite-incorporated hierarchical scaffolds for osteochondral regeneration

Osteochondral regeneration involves the highly challenging and complex reconstruction of cartilage and subchondral bone. Silicon (Si) ions play a crucial role in bone development. Current research on Si ions mainly focuses on bone repair, by using silicate bioceramics with complex ion compositions. However, it is unclear whether the Si ions have important effect on cartilage regeneration. Developing a scaffold that solely releases Si ions to simultaneously promote subchondral bone repair and stimulate cartilage regeneration is critically important. Diatomite (DE) is a natural diatomaceous sediment that can stably release Si ions, known for its abundant availability, low cost, and environmental friendliness. Herein, a hierarchical osteochondral repair scaffold is uniquely designed by incorporating gradient DE into GelMA hydrogel. The adding DE microparticles provides a specific Si source for controlled Si ions release, which not only promotes osteogenic differentiation of rBMSCs (rabbit bone marrow mesenchymal stem cells) but also enhances proliferation and maturation of chondrocytes. Moreover, DE-incorporated hierarchical scaffolds significantly promoted the regeneration of cartilage and subchondral bone. The study suggests the significant role of Si ions in promoting cartilage regeneration and solidifies their foundational role in enhancing bone repair. Furthermore, it offers an economic and eco-friendly strategy for developing high value-added osteochondral regenerative bioscaffolds from low-value ocean natural materials.

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.

scRNA sequencing technology for PitNET studies

Pituitary neuroendocrine tumors (PitNETs) are common, most likely benign tumors with complex clinical characteristics related to hormone hypersecretion and/or growing sellar tumor mass. PitNET types are classified according to their expression of specific transcriptional factors (TFs) and hormone secretion levels. Some types show aggressive, invasive, and reoccurrence behavior. Current research is being conducted to understand the molecular mechanisms regulating these high-heterogeneous neoplasms originating from adenohypophysis, and single-cell RNA sequencing (scRNA-seq) technology is now playing an essential role in these studies due to its remarkable resolution at the single-cell level. This review describes recent studies on human PitNETs performed with scRNA-seq technology, highlighting the potential of this approach in revealing these tumor pathologies, behavior, and regulatory mechanisms.

Obesity disrupts the pituitary-hepatic UPR communication leading to NAFLD progression

Obesity alters levels of pituitary hormones that govern hepatic immune-metabolic homeostasis, dysregulation of which leads to nonalcoholic fatty liver disease (NAFLD). However, the impact of obesity on intra-pituitary homeostasis is largely unknown. Here, we uncovered a blunted unfolded protein response (UPR) but elevated inflammatory signatures in pituitary glands of obese mice and humans. Furthermore, we found that obesity inflames the pituitary gland, leading to impaired pituitary inositol-requiring enzyme 1α (IRE1α)-X-box-binding protein 1 (XBP1) UPR branch, which is essential for protecting against pituitary endocrine defects and NAFLD progression. Intriguingly, pituitary IRE1-deletion resulted in hypothyroidism and suppressed the thyroid hormone receptor B (THRB)-mediated activation of Xbp1 in the liver. Conversely, activation of the hepatic THRB-XBP1 axis improved NAFLD in mice with pituitary UPR defect. Our study provides the first evidence and mechanism of obesity-induced intra-pituitary cellular defects and the pathophysiological role of pituitary-liver UPR communication in NAFLD progression.

Bone mineral density in adults growth hormone deficiency with different ages of onset: a real-world retrospective study

PURPOSE

Bone mineral density (BMD) impairment is one of the critical factors for long-term quality of life in adults growth hormone deficiency (AGHD). This study aims to investigate the annual changes in BMD in AGHD patients with different ages of onset and to identify predicting factors that influence BMD.

METHODS

AGHD patients (n = 160) with available data for 4 years follow-up from a major tertiary medical center in China were retrospectively included (110 [68.8%] childhood-onset, 119 [74.4%] male). BMD of the axial bone (including total hip, neck of femur, and L1-4) derived from dual X-ray absorptiometry and final height were investigated at the first visit, 12 months, 24 months, 36 months, and 48 months thereafter. Low BMD was defined as Z-score ≤ -2.

RESULTS

The prevalence of low BMD was 30.0% at baseline and 12.5% at 4 years of follow-up. The CO AGHD group presented a significantly lower BMD than the AO AGHD group at the baseline (P = 0.009). In contrast, the CO AGHD group had significantly greater median annual BMD change than the AO AGHD group (0.044 vs. -0.0003 g/cm2/year in L1-4, P < 0.001), indicating a significant difference in the overall BMD trend between CO and AO groups. Childhood-onset (odds ratio [OR] 0.326, P = 0.012), low serum testosterone (OR 0.847; P = 0.004) and FT4 (OR 0.595; P = 0.039) level were independent risk factors for BMD loss.

CONCLUSION

The annual changes of BMD show a different pattern in AGHD patients with varying ages of onset. Patients with CO AGHD have a lower bone mass, and in general, appropriate replacement therapy is necessary for long-term bone health in AGHD patients.

scCDC: a computational method for gene-specific contamination detection and correction in single-cell and single-nucleus RNA-seq data

In droplet-based single-cell and single-nucleus RNA-seq assays, systematic contamination of ambient RNA molecules biases the quantification of gene expression levels. Existing methods correct the contamination for all genes globally. However, there lacks specific evaluation of correction efficacy for varying contamination levels. Here, we show that DecontX and CellBender under-correct highly contaminating genes, while SoupX and scAR over-correct lowly/non-contaminating genes. Here, we develop scCDC as the first method to detect the contamination-causing genes and only correct expression levels of these genes, some of which are cell-type markers. Compared with existing decontamination methods, scCDC excels in decontaminating highly contaminating genes while avoiding over-correction of other genes.

Cardiopulmonary progenitors facilitate cardiac repair via exosomal transfer of miR-27b-3p targeting the SIK1-CREB1 axis

Ischemic heart disease, especially myocardial infarction (MI), is one of the leading causes of death worldwide, and desperately needs effective treatments, such as cell therapy. Cardiopulmonary progenitors (CPPs) are stem cells for both heart and lung, but their repairing role in damaged heart is still unknown. Here, we obtained CPPs from E9.5 mouse embryos, maintained their stemness while expanding, and identified their characteristics by scRNA-seq, flow cytometry, quantitative reverse transcription-polymerase chain reaction, and differentiation assays. Moreover, we employed mouse MI model to investigate whether CPPs could repair the injured heart. Our data identified that CPPs exhibit hybrid fibroblastic, endothelial, and mesenchymal state, and they could differentiate into cell lineages within the cardiopulmonary system. Moreover, intramyocardial injection of CPPs improves cardiac function through CPPs exosomes (CPPs-Exo) by promotion of cardiomyocytic proliferation and vascularization. To uncover the underlying mechanism, we used miRNA-seq, bulk RNA-seq, and bioinformatic approaches, and found the highly expressed miR-27b-3p in CPPs-Exo and its target gene Sik1, which can influence the transcriptional activity of CREB1. Therefore, we postulate that CPPs facilitate cardiac repair partially through the SIK1-CREB1 axis via exosomal miR-27b-3p. Our study offers a novel insight into the role of CPPs-Exo in heart repair and highlights the potential of CPPs-Exo as a promising therapeutic strategy for MI.

Single-Cell Transcriptional Profile Construction of Rat Pituitary Glands before and after Sexual Maturation and Identification of Novel Marker Spp1 in Gonadotropes

The mammalian pituitary gland drives highly conserved physiological processes such as somatic cell growth, pubertal transformation, fertility, and metabolism by secreting a variety of hormones. Recently, single-cell transcriptomics techniques have been used in pituitary gland research. However, more studies have focused on adult pituitary gland tissues from different species or different sexes, and no research has yet resolved cellular differences in pituitary gland tissue before and after sexual maturation. Here, we identified a total of 15 cell clusters and constructed single-cell transcriptional profiles of rats before and after sexual maturation. Furthermore, focusing on the gonadotrope cluster, 106 genes were found to be differentially expressed before and after sexual maturation. It was verified that Spp1, which is specifically expressed in gonadotrope cells, could serve as a novel marker for this cell cluster and has a promotional effect on the synthesis and secretion of follicle-stimulating hormone. The results provide a new resource for further resolving the regulatory mechanism of pituitary gland development and pituitary hormone synthesis and secretion.

Single-cell transcriptomics reveal distinct immune-infiltrating phenotypes and macrophage-tumor interaction axes among different lineages of pituitary neuroendocrine tumors

BACKGROUND

Pituitary neuroendocrine tumors (PitNETs) are common gland neoplasms demonstrating distinctive transcription factors. Although the role of immune cells in PitNETs has been widely recognized, the precise immunological environment and its control over tumor cells are poorly understood.

METHODS

The heterogeneity, spatial distribution, and clinical significance of macrophages in PitNETs were analyzed using single-cell RNA sequencing (scRNA-seq), bulk RNA-seq, spatial transcriptomics, immunohistochemistry, and multiplexed quantitative immunofluorescence (QIF). Cell viability, cell apoptosis assays, and in vivo subcutaneous xenograft experiments have confirmed that INHBA-ACVR1B influences the process of tumor cell apoptosis.

RESULTS

The present study evaluated scRNA-seq data from 23 PitNET samples categorized into 3 primary lineages. The objective was to explore the diversity of tumors and the composition of immune cells across these lineages. Analyzed data from scRNA-seq and 365 bulk RNA sequencing samples conducted in-house revealed the presence of three unique subtypes of tumor immune microenvironment (TIME) in PitNETs. These subtypes were characterized by varying levels of immune infiltration, ranging from low to intermediate to high. In addition, the NR5A1 lineage is primarily associated with the subtype characterized by limited infiltration of immune cells. Tumor-associated macrophages (TAMs) expressing CX3CR1+, C1Q+, and GPNMB+ showed enhanced contact with tumor cells expressing NR5A1 + , TBX19+, and POU1F1+, respectively. This emphasizes the distinct interaction axes between TAMs and tumor cells based on their lineage. Moreover, the connection between CX3CR1+ macrophages and tumor cells via INHBA-ACVR1B regulates tumor cell apoptosis.

CONCLUSIONS

In summary, the different subtypes of TIME and the interaction between TAM and tumor cells offer valuable insights into the control of TIME that affects the development of PitNET. These findings can be utilized as prospective targets for therapeutic interventions.

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

We present a gene-level regulatory model, single-cell ATAC + RNA linking (SCARlink), which predicts single-cell gene expression and links enhancers to target genes using multi-ome (scRNA-seq and scATAC-seq co-assay) sequencing data. The approach uses regularized Poisson regression on tile-level accessibility data to jointly model all regulatory effects at a gene locus, avoiding the limitations of pairwise gene-peak correlations and dependence on peak calling. SCARlink outperformed existing gene scoring methods for imputing gene expression from chromatin accessibility across high-coverage multi-ome datasets while giving comparable to improved performance on low-coverage datasets. Shapley value analysis on trained models identified cell-type-specific gene enhancers that are validated by promoter capture Hi-C and are 11× to 15× and 5× to 12× enriched in fine-mapped eQTLs and fine-mapped genome-wide association study (GWAS) variants, respectively. We further show that SCARlink-predicted and observed gene expression vectors provide a robust way to compute a chromatin potential vector field to enable developmental trajectory analysis.

Comprehensive transcriptomic analysis revealing the regulatory dynamics and networks of the pituitary-testis axis in sheep across developmental stages

Spermatogenesis is a complex process intricately regulated by the hypothalamic-pituitary-testis (HPT) axis. However, research on the regulatory factors governing the HPT axis remains limited. This study addresses this gap by conducting a comprehensive analysis of transcriptomes from the pituitary and testis tissues across various developmental stages, encompassing embryonic day (E120), neonatal period (P0), pre-puberty (P90), and post-puberty day (P270). Utilizing edgeR and WGCNA, we identified stage-specific genes in both the pituitary and testis throughout the four developmental stages. Notably, 380, 242, 34, and 479 stage-specific genes were identified in the pituitary, while 886, 297, 201, and 3,678 genes were identified in the testis. Subsequent analyses unveiled associations between these stage-specific genes and crucial pathways such as the cAMP signaling pathway, GnRH secretion, and male gamete generation. Furthermore, leveraging single-cell data from the pituitary and testis, we identified some signaling pathways involving BMP, HGF, IGF, and TGF-β, highlighting mutual regulation between the pituitary and testis at different developmental stages. This study sheds light on the pivotal role of the pituitary-testis axis in the reproductive process of sheep across four distinct developmental stages. Additionally, it delves into the intricate regulatory networks governing reproduction, offering novel insights into the dynamics of the pituitary-testis axis within the reproductive system.

Erythropoietin-producing hepatocellular receptor B6 is highly expressed in non-functioning pituitary neuroendocrine tumors and its expression correlates with tumor size

BACKGROUND

Erythropoietin-producing hepatocellular (EPH) receptors are the largest known family of receptor tyrosine kinases characterized in humans. These proteins are involved in tissue organization, synaptic plasticity, vascular development and the progression of various diseases including cancer. The Erythropoietin-producing hepatocellular receptor tyrosine kinase member EphB6 is a pseudokinase which has not attracted an equivalent amount of interest as its enzymatically-active counterparts. The aim of this study was to assess the expression of EphB6 in pituitary tumors.

METHODS AND RESULTS

Human normal pituitaries and pituitary tumors were examined for EphB6 mRNA expression using real-time PCR and for EphB6 protein by immunohistochemistry and Western blotting. EphB6 was highly expressed in non-functioning pituitary neuroendocrine tumors (NF-PitNETs) versus the normal pituitary and GH-secreting PitNETs. EphB6 mRNA expression was correlated with tumor size.

CONCLUSIONS

Our results suggest EphB6 aberrant expression in NF-PitNETs. Future studies are warranted to determine the role and significance of EphB6 in NF-PitNETs tumorigenesis.

Pituitary stem cells: past, present and future perspectives

Pituitary cells that express the transcription factor SOX2 are stem cells because they can self-renew and differentiate into multiple pituitary hormone-producing cell types as organoids. Wounding and physiological challenges can activate pituitary stem cells, but cell numbers are not fully restored, and the ability to mobilize stem cells decreases with increasing age. The basis of these limitations is still unknown. The regulation of stem cell quiescence and activation involves many different signalling pathways, including those mediated by WNT, Hippo and several cytokines; more research is needed to understand the interactions between these pathways. Pituitary organoids can be formed from human or mouse embryonic stem cells, or from human induced pluripotent stem cells. Human pituitary organoid transplantation is sufficient to induce corticosterone release in hypophysectomized mice, raising the possibility of therapeutic applications. Today, pituitary organoids have the potential to assess the role of individual genes and genetic variants on hormone production ex vivo, providing an important tool for the advancement of exciting frontiers in pituitary stem cell biology and pituitary organogenesis. In this article, we provide an overview of notable discoveries in pituitary stem cell function and highlight important areas for future research.

Pathogenesis, clinical features, and treatment of plurihormonal pituitary adenoma

Plurihormonal pituitary adenoma (PPA) is a type of pituitary tumor capable of producing two or more hormones and usually presents as an aggressive, large adenoma. As yet, its pathogenesis remains unclear. This is the first study to systematically summarize the underlying pathogenesis of PPA. The pathogenesis is related to plurihormonal primordial stem cells, co-transcription factors, hormone co-expression, differential gene expression, and cell transdifferentiation. We conducted a literature review of PPA and analyzed its clinical characteristics. We found that the average age of patients with PPA was approximately 40 years, and most showed only one clinical symptom. The most common manifestation was acromegaly. Currently, PPA is treated with surgical resection. However, recent studies suggest that immunotherapy may be a potentially effective treatment.

Single-cell transcriptomic analysis reveals tumor cell heterogeneity and immune microenvironment features of pituitary neuroendocrine tumors

BACKGROUND

Pituitary neuroendocrine tumors (PitNETs) are one of the most common types of intracranial tumors. Currently, the cellular characteristics of normal pituitary and various other types of PitNETs are still not completely understood.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) on 4 normal samples and 24 PitNET samples for comprehensive bioinformatics analysis. Findings regarding the function of PBK in the aggressive tumor cells were validated by siRNA knockdown, overexpression, and transwell experiments.

RESULTS

We first constructed a reference cell atlas of the human pituitary. Subsequent scRNA-seq analysis of PitNET samples, representing major tumor subtypes, shed light on the intrinsic cellular heterogeneities of the tumor cells and tumor microenvironment (TME). We found that the expression of hormone-encoding genes defined the major variations of the PIT1-lineage tumor cell transcriptomic heterogeneities. A sub-population of TPIT-lineage tumor cells highly expressing GZMK suggested a novel subtype of corticotroph tumors. In immune cells, we found two clusters of tumor-associated macrophages, which were both highly enriched in PitNETs but with distinct functional characteristics. In PitNETs, the stress response pathway was significantly activated in T cells. While a majority of these tumors are benign, our study unveils a common existence of aggressive tumor cells in the studied samples, which highly express a set of malignant signature genes. The following functional experiments confirmed the oncogenic role of selected up-regulated genes. The over-expression of PBK could promote both tumor cell proliferation and migration, and it was also significantly associated with poor prognosis in PitNET patients.

CONCLUSIONS

Our data and analysis manifested the basic cell types in the normal pituitary and inherent heterogeneity of PitNETs, identified several features of the tumor immune microenvironments, and found a novel epithelial cell sub-population with aggressive signatures across all the studied cases.

Stem cell-associated transcription factors in non-functioning pituitary neuroendocrine tumours

Background: Cells with stem cell features have been described in pituitary neuroendocrine tumours (PitNETs). Transcription factors SOX2 and SOX9 are stem cell-associated markers while the pituitary progenitor marker PROP1 is involved in anterior pituitary development. We characterised the presence of these markers known to be present in the human pituitary in non-functioning (NF) PitNETs. Methods: We investigated the pituitary transcription factors SOX2, SOX9 and PROP1 by immunohistochemistry (IHC) (N = 125) and RT-qPCR (N = 78) in a retrospective cohort of clinically NF-PitNETs. The markers were scored based on the percentage of immunolabeled cells. IHC staining scores were compared to reintervention rates for the whole cohort, and to expression of FSH, LH or ER in gonadotroph NF-PitNETs. Results: Most tumours showed no or few cells positive for SOX2, SOX9 and PROP1. More patients with SOX2-negative tumours went through reintervention (40 % vs 19 %, p = 0.03). SOX2, SOX9 and PROP1 staining correlated positively to each other (SOX2 and SOX9 rs = 0.666, SOX2 and PROP1 rs = 0.704, SOX9 and PROP1 rs = 0.570, and p < 0.001 for all). In gonadotroph NF-PitNETs, staining for SOX2 and PROP1 was positively associated to FSHβ staining (p < 0.001 for both). Staining for SOX2, SOX9 and PROP1 was positively associated with gene expression of Estrogen Receptor 1 (ESR1) (p < 0.001, p = 0.004 and p < 0.001) and IHC staining for ERα (p = 0.001, p = 0.03 and p = 0.05, respectively). Conclusion: SOX2, SOX9 and PROP1 were present at low levels in NF-PitNETs. Absence of SOX2 staining was associated with a higher reintervention rate. The stem cell markers correlated positively with markers of gonadotroph differentiation in gonadotroph NF-PitNETs. SOX2 and SOX9 were frequently coexpressed and showed positivity in intratumoural cells with epithelial features, however without coexpression of pituitary transcription factors.

Neural plate progenitors give rise to both anterior and posterior pituitary cells

The pituitary is the master neuroendocrine gland, which regulates body homeostasis. It consists of the anterior pituitary/adenohypophysis harboring hormones producing cells and the posterior pituitary/neurohypophysis, which relays the passage of hormones from the brain to the periphery. It is accepted that the adenohypophysis originates from the oral ectoderm (Rathke's pouch), whereas the neural ectoderm contributes to the neurohypophysis. Single-cell transcriptomics of the zebrafish pituitary showed that cyp26b1-positive astroglial pituicytes of the neurohypophysis and prop1-positive adenohypophyseal progenitors expressed common markers implying lineage relatedness. Genetic tracing identifies that, in contrast to the prevailing dogma, neural plate precursors of zebrafish (her4.3+) and mouse (Sox1+) contribute to both neurohypophyseal and a subset of adenohypophyseal cells. Pituicyte-derived retinoic-acid-degrading enzyme Cyp26b1 fine-tunes differentiation of prop1+ progenitors into hormone-producing cells. These results challenge the notion that adenohypophyseal cells are exclusively derived from non-neural ectoderm and demonstrate that crosstalk between neuro- and adeno-hypophyseal cells affects differentiation of pituitary cells.

Primary Pituitary Carcinoids Do Not Exist: A Reappraisal in the Era of Pituitary Neuroendocrine Tumours

The World Health Organization classification of pituitary tumours, published in 2022, supported a change in the terminology from "pituitary adenoma" to "pituitary neuroendocrine tumour" (PitNET). The neuroendocrine cells represent an integral part of the diffuse neuroendocrine system, including, among others, thyroid C cells, the parathyroid chief cells, and the anterior pituitary. Normal and neoplastic adenohypophyseal neuroendocrine cells have light microscopic, ultrastructural features and an immunoprofile compatible with the neuroendocrine cells and neuroendocrine tumours from other organs. Moreover, neuroendocrine cells of pituitary origin express transcription factors which indicate their cell-lineage origin. Thus, pituitary tumours are now considered as a continuum with other neuroendocrine tumours. PitNETs may occasionally be aggressive. In this context, the term "pituitary carcinoid" has no specific meaning: it either represents a PitNET, or a metastasis to the pituitary gland of a neuroendocrine tumour (NET). An accurate pathological evaluation, combined where necessary with functional radionuclide imaging, can define the origin of the tumour. We recommend that clinicians liaise with patient groups to understand the terminology to define primary tumours of adenohypophyseal cells. It is incumbent upon the responsible clinician to explain the use of the word "tumour" in a given clinical context.

Single-cell transcriptomic analysis reveals rich pituitary-Immune interactions under systemic inflammation

The pituitary represents an essential hub in the hypothalamus-pituitary-adrenal (HPA) axis. Pituitary hormone-producing cells (HPCs) release several hormones to regulate fundamental bodily functions under normal and stressful conditions. It is well established that the pituitary endocrine gland modulates the immune system by releasing adrenocorticotropic hormone (ACTH) in response to neuronal activation in the hypothalamus. However, it remains unclear how systemic inflammation regulates the transcriptomic profiles of pituitary HPCs. Here, we performed single-cell RNA-sequencing (scRNA-seq) of the mouse pituitary and revealed that upon inflammation, all major pituitary HPCs respond robustly in a cell type-specific manner, with corticotropes displaying the strongest reaction. Systemic inflammation also led to the production and release of noncanonical bioactive molecules, including Nptx2 by corticotropes, to modulate immune homeostasis. Meanwhile, HPCs up-regulated the gene expression of chemokines that facilitated the communication between the HPCs and immune cells. Together, our study reveals extensive interactions between the pituitary and immune system, suggesting multifaceted roles of the pituitary in mediating the effects of inflammation on many aspects of body physiology.

SOX9-positive pituitary stem cells differ according to their position in the gland and maintenance of their progeny depends on context

Stem cell (SC) differentiation and maintenance of resultant progeny underlie cell turnover in many organs, but it is difficult to pinpoint the contribution of either process. In the pituitary, a central regulator of endocrine axes, adult SCs undergo activation after target organ ablation, providing a well-characterized paradigm to study an adaptative response in a multi-organ system. Here, we used single-cell technologies to characterize SC heterogeneity and mobilization together with lineage tracing. We show that SC differentiation occurs more frequently than thought previously. In adaptative conditions, differentiation increases and is more diverse than demonstrated by the lineage tracing experiments. Detailed examination of SC progeny suggests that maintenance of selected nascent cells underlies SC output, highlighting a trophic role for the microenvironment. Analyses of cell trajectories further predict pathways and potential regulators. Our model provides a valuable system to study the influence of evolving states on the mechanisms of SC mobilization.

Characterization of Somatotrope Cell Expansion in Response to GHRH in the Neonatal Mouse Pituitary

In humans and mice, loss-of-function mutations in growth hormone-releasing hormone receptor (GHRHR) cause isolated GH deficiency. The mutant GHRHR mouse model, GhrhrLit/Lit (LIT), exhibits loss of serum GH, but also fewer somatotropes. However, how loss of GHRH signaling affects expansion of stem and progenitor cells giving rise to GH-producing cells is unknown. LIT mice and wild-type littermates were examined for differences in proliferation and gene expression of pituitary lineage markers by quantitative reverse transcription polymerase chain reaction and immunohistochemistry at postnatal day 5 (p5) and 5 weeks. At p5, the LIT mouse shows a global decrease in pituitary proliferation measured by proliferation marker Ki67 and phospho-histone H3. This proliferative defect is seen in a pituitary cell expressing POU1F1 with or without GH. SOX9-positive progenitors show no changes in proliferation in p5 LIT mice. Additionally, the other POU1F1 lineage cells are not decreased in number; rather, we observe an increase in lactotrope cell population as well as messenger RNA for Tshb and Prl. In the 5-week LIT pituitary, the proliferative deficit in POU1F1-expressing cells observed neonatally persists, while the number and proliferative proportion of SOX9 cells do not appear changed. Treatment of cultured pituitary explants with GHRH promotes proliferation of POU1F1-expressing cells, but not GH-positive cells, in a mitogen-activated protein kinase-dependent manner. These findings indicate that hypothalamic GHRH targets proliferation of a POU1F1-positive cell, targeted to the somatotrope lineage, to fine tune their numbers.

A Single-Cell Atlas of an Early Mongolian Sheep Embryo

Cell types have been established during organogenesis based on early mouse embryos. However, our understanding of cell types and molecular mechanisms in the early embryo development of Mongolian sheep has been hampered. This study presents the first comprehensive single-cell transcriptomic characterization at E16 in Ujumqin sheep and Hulunbuir short-tailed sheep. Thirteen major cell types were identified at E16 in Ujumqin sheep, and eight major cell types were identified at E16 in Hulunbuir short-tailed sheep. Function enrichment analysis showed that several pathways were significantly enriched in the TGF-beta signaling pathway, the Hippo signaling pathway, the platelet activation pathway, the riboflavin metabolism pathway, the Wnt signaling pathway, regulation of the actin cytoskeleton, and the insulin signaling pathway in the notochord cluster. Glutathione metabolism, glyoxylate, and dicarboxylate metabolism, the citrate cycle, thyroid hormone synthesis, pyruvate metabolism, cysteine and methionine metabolism, thermogenesis, and the VEGF signaling pathway were significantly enriched in the spinal cord cluster. Steroid biosynthesis, riboflavin metabolism, the cell cycle, the Hippo signaling pathway, the Hedgehog signaling pathway, the FoxO signaling pathway, the JAK-STAT signaling pathway, and the Wnt signaling pathway were significantly enriched in the paraxial mesoderm cluster. The notochord cluster, spinal cord cluster, and paraxial mesoderm cluster were found to be highly associated with tail development. Pseudo-time analysis demonstrated that the mesenchyme can translate to the notochord in Ujumqin sheep. Molecular assays revealed that the Hippo signaling pathway was enriched in Ujumqin sheep. This comprehensive single-cell map revealed previously unrecognized signaling pathways that will further our understanding of the mechanism of short-tailed sheep formation.

RNA-seq analysis identifies age-dependent changes in expression of mRNAs - encoding N-glycosylation pathway enzymes in mouse gonadotropes

Follicle-stimulating hormone (FSH) is a glycoprotein that is assembled as a heterodimer of α/β subunits in gonadotropes. Each subunit contains two N-glycan chains. Our previous in vivo genetic studies identified that at least one N-glycan chain must be present on the FSHβ subunit for efficient FSH dimer assembly and secretion. Moreover, macroheterogeneity observed uniquely on human FSHβ results in ratiometric changes in age-specific FSH glycoforms, particularly during menopausal transition. Despite the recognition of many prominent roles of sugars on FSH including dimer assembly and secretion, serum half-life, receptor binding and signal transduction, the N-glycosylation machinery in gonadotropes has never been defined. Here, we used a mouse model in which gonadotropes are GFP-labeled in vivo and achieved rapid purification of GFP+ gonadotropes from pituitaries of female mice at reproductively young, middle, and old ages. We identified by RNA-seq analysis 52 mRNAs encoding N-glycosylation pathway enzymes expressed in 3- and 8-10-month-old mouse gonadotropes. We hierarchically mapped and localized the enzymes to distinct subcellular organelles within the N-glycosylation biosynthetic pathway. Of the 52 mRNAs, we found 27 mRNAs are differentially expressed between the 3- and 8-10-month old mice. We subsequently selected 8 mRNAs which showed varying changes in expression for confirmation of abundance in vivo via qPCR analysis, using more expanded aging time points with distinct 8-month and 14-month age groups. Real time qPCR analysis indicated dynamic changes in expression of N-glycosylation pathway enzyme-encoding mRNAs across the life span. Notably, computational analysis predicted the promoters of genes encoding these 8 mRNAs contain multiple high probability binding sites for estrogen receptor-1 and progesterone receptor. Collectively, our studies define the N-glycome and identify age-specific dynamic changes in mRNAs encoding N-glycosylation pathway enzymes in mouse gonadotropes. Our studies suggest the age-related decline in ovarian steroids may regulate expression of N-glycosylation enzymes in mouse gonadotropes and explain the age-related N-glycosylation shift previously observed on human FSHβ subunit in pituitaries of women.