About

Stephen Aaron Liebhaber, M.D., is an Emeritus Professor of Genetics at the University of Pennsylvania's Perelman School of Medicine. His research focuses on the roles of chromatin structure and epigenetic controls in eukaryotic gene activation, as well as the mechanisms of mRNA-protein interactions in the regulation of mRNA stability and expression. His work investigates how epigenetic modifications in histones and DNA influence gene expression profiles, particularly through long-range control of gene promoters, using the human Growth Hormone gene cluster as a model. Liebhaber's studies employ transgenic mouse models, in vitro chromatin analyses, and cell culture models to understand developmentally controlled gene regulation. He also explores post-transcriptional gene regulation, including the stabilization of mRNAs by RNA-binding proteins such as αCPs, and the pathways involved in mRNA decay and surveillance. His research extends to activation of placental hormone gene expression, especially within the human Growth Hormone locus, utilizing transgenic mice and stem cell models to identify chromatin-based determinants of tissue-specific and developmental gene expression. Liebhaber has contributed extensively to the understanding of gene regulation mechanisms in development and differentiation, and his laboratory maintains a focus on chromatin modifications, epigenetic controls, and mRNA stability, with ongoing projects involving transgenic models, stem cell cultures, and biochemical approaches.

Research topics

• Endocrinology
• Genetics
• Internal medicine
• Biology
• Cell biology

Selected publications

• Multi-omic profiling of pituitary thyrotropic cells and progenitors

  BMC Biology · 2021-04-15 · 8 citations

  articleOpen access

  BACKGROUND: The pituitary gland is a neuroendocrine organ containing diverse cell types specialized in secreting hormones that regulate physiology. Pituitary thyrotropes produce thyroid-stimulating hormone (TSH), a critical factor for growth and maintenance of metabolism. The transcription factors POU1F1 and GATA2 have been implicated in thyrotrope fate, but the transcriptomic and epigenomic landscapes of these neuroendocrine cells have not been characterized. The goal of this work was to discover transcriptional regulatory elements that drive thyrotrope fate. RESULTS: We identified the transcription factors and epigenomic changes in chromatin that are associated with differentiation of POU1F1-expressing progenitors into thyrotropes using cell lines that represent an undifferentiated Pou1f1 lineage progenitor (GHF-T1) and a committed thyrotrope line that produces TSH (TαT1). We compared RNA-seq, ATAC-seq, histone modification (H3K27Ac, H3K4Me1, and H3K27Me3), and POU1F1 binding in these cell lines. POU1F1 binding sites are commonly associated with bZIP transcription factor consensus binding sites in GHF-T1 cells and Helix-Turn-Helix (HTH) or basic Helix-Loop-Helix (bHLH) factors in TαT1 cells, suggesting that these classes of transcription factors may recruit or cooperate with POU1F1 binding at unique sites. We validated enhancer function of novel elements we mapped near Cga, Pitx1, Gata2, and Tshb by transfection in TαT1 cells. Finally, we confirmed that an enhancer element near Tshb can drive expression in thyrotropes of transgenic mice, and we demonstrate that GATA2 enhances Tshb expression through this element. CONCLUSION: These results extend the ENCODE multi-omic profiling approach to the pituitary gland, which should be valuable for understanding pituitary development and disease pathogenesis.

  PublisherOA PDFDOI

• RNA-Binding Proteins PCBP1 and PCBP2 Are Critical Determinants of Murine Erythropoiesis

  Molecular and Cellular Biology · 2021-06-28 · 18 citations

  articleOpen accessSenior author

  analyses of primary erythroid progenitors selectively depleted of these two RNA-binding proteins revealed that they mediate a combination of overlapping and isoform-specific impacts on hematopoietic lineage transcriptome, impacting both mRNA representation and exon splicing. These data lead us to conclude that PCBP1 and PCBP2 mediate functions critical to differentiation of the erythroid lineage.

  PublisherOA PDFDOI

• Dynamic changes in RNA–protein interactions and RNA secondary structure in mammalian erythropoiesis

  Life Science Alliance · 2021-07-27 · 7 citations

  articleOpen accessCorresponding

  Two features of eukaryotic RNA molecules that regulate their post-transcriptional fates are RNA secondary structure and RNA-binding protein (RBP) interaction sites. However, a comprehensive global overview of the dynamic nature of these sequence features during erythropoiesis has never been obtained. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approach to reveal the global landscape of RNA secondary structure and RBP-RNA interaction sites and the dynamics of these features during this important developmental process. We identify dynamic patterns of RNA secondary structure and RBP binding throughout the process and determine a set of corresponding protein-bound sequence motifs along with their dynamic structural and RBP-binding contexts. Finally, using these dynamically bound sequences, we identify a number of RBPs that have known and putative key functions in post-transcriptional regulation during mammalian erythropoiesis. In total, this global analysis reveals new post-transcriptional regulators of mammalian blood cell development.

  PublisherOA PDFDOI

• Single-cell transcriptomic analysis of adult mouse pituitary reveals sexual dimorphism and physiologic demand-induced cellular plasticity

  Protein & Cell · 2020 · 89 citations

  Senior authorCorresponding
  • Biology
  • Cell biology
  • Endocrinology

  The anterior pituitary gland drives highly conserved physiologic processes in mammalian species. These hormonally controlled processes are central to somatic growth, pubertal transformation, fertility, lactation, and metabolism. Current cellular models of mammalian anteiror pituitary, largely built on candidate gene based immuno-histochemical and mRNA analyses, suggest that each of the seven hormones synthesized by the pituitary is produced by a specific and exclusive cell lineage. However, emerging evidence suggests more complex relationship between hormone specificity and cell plasticity. Here we have applied massively parallel single-cell RNA sequencing (scRNA-seq), in conjunction with complementary imaging-based single-cell analyses of mRNAs and proteins, to systematically map both cell-type diversity and functional state heterogeneity in adult male and female mouse pituitaries at single-cell resolution and in the context of major physiologic demands. These quantitative single-cell analyses reveal sex-specific cell-type composition under normal pituitary homeostasis, identify an array of cells associated with complex complements of hormone-enrichment, and undercover non-hormone producing interstitial and supporting cell-types. Interestingly, we also identified a Pou1f1-expressing cell population that is characterized by a unique multi-hormone gene expression profile. In response to two well-defined physiologic stresses, dynamic shifts in cellular diversity and transcriptome profiles were observed for major hormone producing and the putative multi-hormone cells. These studies reveal unanticipated cellular complexity and plasticity in adult pituitary, and provide a rich resource for further validating and expanding our molecular understanding of pituitary gene expression programs and hormone production.

  PublisherOA PDFDOI

• Identification of pituitary thyrotrope signature genes and regulatory elements

  bioRxiv (Cold Spring Harbor Laboratory) · 2020-08-06

  preprintOpen access

  Abstract Pituitary thyrotropes are specialized cells that produce thyroid stimulating hormone (TSH), a critical factor for growth and maintenance of metabolism. The transcription factors POU1F1 and GATA2 have been implicated in thyrotrope fate and transcriptional regulation of the beta subunit of TSH, Tshb , but no transcriptomic or epigenomic analyses of these cells has been undertaken. The goal of this work was to discover key transcriptional regulatory elements that drive thyrotrope fate. We identified the transcription factors and epigenomic changes in chromatin that are associated with differentiation of POU1F1-expressing progenitors into thyrotropes, a process modeled by two cell lines: one that represents an early, undifferentiated Pou1f1 lineage progenitor (GHF-T1) and one that is a committed thyrotrope that produces TSH (TαT1). We generated and compared RNA-seq, ATAC-seq, histone modification (including H3K27Ac, H3K4Me1, and H3K27Me3), and transcription factor (POU1F1) binding in these two cell lines to identify regulatory elements and candidate transcriptional regulators. We identified POU1F1 binding sites that were unique to each cell line. POU1F1 binding sites are commonly associated with bZIP transcription factor consensus binding sites in GHF-T1 cells and Helix-Turn-Helix (HTH) or basic Helix-Loop-Helix (bHLH) factors in TαT1 cells, suggesting that these classes of transcription factors may recruit or cooperate with POU1F1 binding to unique sites. We validated enhancer function of novel elements we mapped near Cga, Pitx1, Gata2, and Tshb by transfection in TαT1 cells. Finally, we confirmed that an enhancer element near Tshb can drive expression in thyrotropes of transgenic mice, and we demonstrate that GATA2 enhances Tshb expression through this element. These results extend the ENCODE multi-omic profiling approach to an organ that is critical for growth and metabolism, which should be valuable for understanding pituitary development and disease pathogenesis.

  PublisherOA PDFDOI

• The Transcription Factor NR4A2 Plays an Essential Role in Driving Prolactin Expression in Female Pituitary Lactotropes

  Endocrinology · 2020-03-19 · 8 citations

  articleOpen accessSenior author

  Differentiation of the hormone-producing cells of the pituitary represents an informative model of cell fate determination. The generation and maintenance of 2 pituitary lineages, the growth hormone (GH)- producing somatotropes and the prolactin (PRL)- producing lactotropes, are dependent on the pituitary-specific transcription factor, POU1F1. While POU1F1 is expressed in both cell types, and plays a role in activation of both the Gh and Prl genes, expression of Gh and Prl is restricted to somatotropes and lactotropes, respectively. These observations imply the existence of additional factors that contribute to the somatotrope and lactotrope identities and their hormone expressions. Prior transcriptome analysis of primary somatotropes and lactotropes isolated from the mouse pituitary identified enrichment of a transcription factor, Nr4a2, in the lactotropes. Nr4a2 was shown in a cell culture model to bind the Prl promoter at a position adjacent to Pou1f1 and to synergize with Pou1f1 in driving Prl transcription. Here we demonstrate in vivo the role of Nr4a2 as an enhancer of Prl expression by conditional gene inactivation of the Nr4a2 gene in mouse lactotropes. We demonstrate that nuclear orphan receptor transcription factor (NR4A2) binding at the Prl promoter is dependent on actions of POU1F1; while POU1F1 is essential to loading polymerase (Pol) II on the Prl promoter, Nr4a2 plays a role in enhancing Pol II release into the Prl gene body. These studies establish an in vivo role of Nr4a2 in enhancing Prl expression in mouse lactotropes, explore its mechanism of action, and establish a system for further study of the lactotrope lineage in the pituitary.

  PublisherOA PDFDOI

• A cytosine-rich splice regulatory determinant enforces functional processing of the human α-globin gene transcript

  Blood · 2019-03-05 · 4 citations

  articleOpen accessSenior author

  Abstract The establishment of efficient and stable splicing patterns in terminally differentiated cells is critical to maintenance of specific functions throughout the lifespan of an organism. The human α-globin (hα-globin) gene contains 3 exons separated by 2 short introns. Naturally occurring α-thalassemia mutations that trigger aberrant splicing have revealed the presence of cryptic splice sites within the hα-globin gene transcript. How cognate (functional) splice sites are selectively used in lieu of these cryptic sites has remained unexplored. Here we demonstrate that the preferential selection of a cognate splice donor essential to functional splicing of the hα-globin transcript is dependent on the actions of an intronic cytosine (C)-rich splice regulatory determinant and its interacting polyC-binding proteins. Inactivation of this determinant by mutation of the C-rich element or by depletion of polyC-binding proteins triggers a dramatic shift in splice donor activity to an upstream, out-of-frame, cryptic donor. The essential role of the C-rich element in hα-globin gene expression is supported by its coevolution with the cryptic donor site in primate species. These data lead us to conclude that an intronic C-rich determinant enforces functional splicing of the hα-globin transcript, thus acting as an obligate determinant of hα-globin gene expression.

  PublisherOA PDFDOI

• Poly(C)-Binding Protein Pcbp2 Enables Differentiation of Definitive Erythropoiesis by Directing Functional Splicing of the Runx1 Transcript

  Molecular and Cellular Biology · 2018-06-05 · 18 citations

  articleSenior authorCorresponding

  Formation of the mammalian hematopoietic system is under a complex set of developmental controls. Here, we report that mouse embryos lacking the KH domain poly(C) binding protein, Pcbp2, are selectively deficient in the definitive erythroid lineage. Compared to wild-type controls, transcript splicing analysis of the Pcbp2−/− embryonic liver reveals accentuated exclusion of an exon (exon 6) that encodes a highly conserved transcriptional control segment of the hematopoietic master regulator, Runx1. Embryos rendered homozygous for a Runx1 locus lacking this cassette exon (Runx1ΔE6) effectively phenocopy the loss of the definitive erythroid lineage in Pcbp2−/− embryos. These data support a model in which enhancement of Runx1 cassette exon 6 inclusion by Pcbp2 serves a critical role in development of hematopoietic progenitors and constitutes a critical step in the developmental pathway of the definitive erythropoietic lineage.

  PublisherDOI

• Domain-focused CRISPR screen identifies HRI as a fetal hemoglobin regulator in human erythroid cells

  Science · 2018-07-19 · 153 citations

  article

  Increasing fetal hemoglobin (HbF) levels in adult red blood cells provides clinical benefit to patients with sickle cell disease and some forms of β-thalassemia. To identify potentially druggable HbF regulators in adult human erythroid cells, we employed a protein kinase domain-focused CRISPR-Cas9-based genetic screen with a newly optimized single-guide RNA scaffold. The screen uncovered the heme-regulated inhibitor HRI (also known as EIF2AK1), an erythroid-specific kinase that controls protein translation, as an HbF repressor. HRI depletion markedly increased HbF production in a specific manner and reduced sickling in cultured erythroid cells. Diminished expression of the HbF repressor BCL11A accounted in large part for the effects of HRI depletion. Taken together, these results suggest HRI as a potential therapeutic target for hemoglobinopathies.

  PublisherDOI

• Single cell transcriptomic analysis of the adult mouse pituitary reveals a novel multi-hormone cell cluster and physiologic demand-induced lineage plasticity

  bioRxiv (Cold Spring Harbor Laboratory) · 2018-11-22 · 8 citations

  preprintOpen accessSenior author

  Abstract The anterior pituitary gland drives a set of highly conserved physiologic processes in mammalian species. These hormonally-controlled processes are central to somatic growth, pubertal transformation, fertility, lactation, and metabolism. Current models, largely built upon candidate gene based immuno-histochemical and mRNA analyses, suggest that each of the seven hormones synthesized by the pituitary is produced by a specific and exclusive cell lineage. However, emerging evidence suggests more complex models of hormone specificity and cell plasticity. Here we have applied massively parallel single-cell RNA sequencing (scRNA-seq), in conjunction with a set of orthogonal mRNA and protein imaging studies, to systematically map the cellular composition of adult male and female mouse pituitaries at single-cell resolution and in the setting of major physiologic demands. These analyses reveal sex-specific cellular diversity associated with normal pituitary homeostasis, and identify an array of cells with complex complements of hormone-enrichment as well as a series of non-hormone producing interstitial and supporting cell lineages. These scRNA-seq studies identify a major cell population that is characterized by a unique multi-hormone gene expression profile. The detection of dynamic shifts in cellular representations and transcriptome profiles in response to two well-defined physiologic stresses suggests corresponding roles of a number of these clusters in cellular plasticity within the adult pituitary. These studies point to an unanticipated complexity and plasticity in pituitary cellular composition that expands upon current models and concepts of pituitary gene expression and hormone production.

  PublisherOA PDFDOI

Recent grants

• NIH Grant P60HL038632

  NIH · $32.1M · 2003

• Determinants of Human Growth Hormone Expression and Pituitary Cell Differentiation

  NIH · $1.6M · 2016–2020

• University of Pennsylvania Diabetes Research Center

  NIH · $17.0M · 1997–2027

• NIH Grant R01HD046737

  NIH · $4.1M · 2015

• NIH Grant R01AM033975

  NIH · $108k · 1987

Frequent coauthors

• Nancy E. Cooke

  155 shared

• Xinjun Ji

  University of Pennsylvania

  44 shared

• Jharna Ray

  University of Calcutta

  34 shared

• Janice Russell

  University of Pennsylvania

  30 shared

• James N. MacLeod

  University of Kentucky

  29 shared

• Susan P. Perrine

  Boston University

  25 shared

• Millicent Sutton

  25 shared

• Yelena Galperin

  Icahn School of Medicine at Mount Sinai

  25 shared