A developmental lineage-based gene co-expression network for mouse pancreatic β-cells reveals a role for Zfp800 in pancreas development

Deletion of Ascl1 in pancreatic β-cells improves insulin secretion, promotes parasympathetic innervation, and attenuates dedifferentiation during metabolic stress

OBJECTIVE

ASCL1, a pioneer transcription factor, is essential for neural cell differentiation and function. Previous studies have shown that Ascl1 expression is increased in pancreatic β-cells lacking functional KATP channels or after feeding of a high fat diet (HFD) suggesting that it may contribute to the metabolic stress response of β-cells.

METHODS

We generated β-cell-specific Ascl1 knockout mice (Ascl1βKO) and assessed their glucose homeostasis, islet morphology and gene expression after feeding either a normal diet or HFD for 12 weeks, or in combination with a genetic disruption of Abcc8, an essential KATP channel component.

RESULTS

Ascl1 expression is increased in response to both a HFD and membrane depolarization and requires CREB-dependent Ca2+ signaling. No differences in glucose homeostasis or islet morphology were observed in Ascl1βKO mice fed a normal diet or in the absence of KATP channels. However, male Ascl1βKO mice fed a HFD exhibited decreased blood glucose levels, improved glucose tolerance, and increased β-cell proliferation. Bulk RNA-seq analysis of islets from Ascl1βKO mice from three studied conditions showed alterations in genes associated with the secretory function. HFD-fed Ascl1βKO mice showed the most extensive changes with increased expression of genes necessary for glucose sensing, insulin secretion and β-cell proliferation, and a decrease in genes associated with β-cell dysfunction, inflammation and dedifferentiation. HFD-fed Ascl1βKO mice also displayed increased expression of parasympathetic neural markers and cholinergic receptors that was accompanied by increased insulin secretion in response to acetylcholine and an increase in islet innervation.

CONCLUSIONS

Ascl1 expression is induced by stimuli that cause Ca2+-signaling to the nucleus and contributes in a multifactorial manner to the loss of β-cell function by promoting the expression of genes associated with cellular dedifferentiation, attenuating β-cells proliferation, suppressing acetylcholine sensitivity, and repressing parasympathetic innervation of islets. Thus, the removal of Ascl1 from β-cells improves their function in response to metabolic stress.

Unbiased transcription factor CRISPR screen identifies ZNF800 as master repressor of enteroendocrine differentiation

Enteroendocrine cells (EECs) are hormone-producing cells residing in the epithelium of stomach, small intestine (SI), and colon. EECs regulate aspects of metabolic activity, including insulin levels, satiety, gastrointestinal secretion, and motility. The generation of different EEC lineages is not completely understood. In this work, we report a CRISPR knockout screen of the entire repertoire of transcription factors (TFs) in adult human SI organoids to identify dominant TFs controlling EEC differentiation. We discovered ZNF800 as a master repressor for endocrine lineage commitment, which particularly restricts enterochromaffin cell differentiation by directly controlling an endocrine TF network centered on PAX4. Thus, organoid models allow unbiased functional CRISPR screens for genes that program cell fate.

Caloric restriction promotes beta cell longevity and delays aging and senescence by enhancing cell identity and homeostasis mechanisms

Caloric restriction (CR) extends organismal lifespan and health span by improving glucose homeostasis mechanisms. How CR affects organellar structure and function of pancreatic beta cells over the lifetime of the animal remains unknown. Here, we used single nucleus transcriptomics to show that CR increases the expression of genes for beta cell identity, protein processing, and organelle homeostasis. Gene regulatory network analysis link this transcriptional phenotype to transcription factors involved in beta cell identity (Mafa) and homeostasis (Atf6). Imaging metabolomics further demonstrates that CR beta cells are more energetically competent. In fact, high-resolution light and electron microscopy indicates that CR reduces beta cell mitophagy and increases mitochondria mass, increasing mitochondrial ATP generation. Finally, we show that long-term CR delays the onset of beta cell aging and senescence to promote longevity by reducing beta cell turnover. Therefore, CR could be a feasible approach to preserve compromised beta cells during aging and diabetes.

Caloric restriction promotes beta cell longevity and delays aging and senescence by enhancing cell identity and homeostasis mechanisms

Caloric restriction (CR) extends organismal lifespan and health span by improving glucose homeostasis mechanisms. How CR affects organellar structure and function of pancreatic beta cells over the lifetime of the animal remains unknown. Here, we used single nucleus transcriptomics to show that CR increases the expression of genes for beta cell identity, protein processing, and organelle homeostasis. Gene regulatory network analysis link this transcriptional phenotype to transcription factors involved in beta cell identity (Mafa) and homeostasis (Atf6). Imaging metabolomics further demonstrates that CR beta cells are more energetically competent. In fact, high-resolution light and electron microscopy indicates that CR reduces beta cell mitophagy and increases mitochondria mass, increasing mitochondrial ATP generation. Finally, we show that long-term CR delays the onset of beta cell aging and senescence to promote longevity by reducing beta cell turnover. Therefore, CR could be a feasible approach to preserve compromised beta cells during aging and diabetes.

Single-Cell RNA Sequencing of Sox17-Expressing Lineages Reveals Distinct Gene Regulatory Networks and Dynamic Developmental Trajectories

During early embryogenesis, the transcription factor SOX17 contributes to hepato-pancreato-biliary system formation and vascular-hematopoietic emergence. To better understand Sox17 function in the developing endoderm and endothelium, we developed a dual-color temporal lineage-tracing strategy in mice combined with single-cell RNA sequencing to analyze 6934 cells from Sox17-expressing lineages at embryonic days 9.0-9.5. Our analyses showed 19 distinct cellular clusters combined from all 3 germ layers. Differential gene expression, trajectory and RNA-velocity analyses of endothelial cells revealed a heterogenous population of uncommitted and specialized endothelial subtypes, including 2 hemogenic populations that arise from different origins. Similarly, analyses of posterior foregut endoderm revealed subsets of hepatic, pancreatic, and biliary progenitors with overlapping developmental potency. Calculated gene-regulatory networks predict gene regulons that are dominated by cell type-specific transcription factors unique to each lineage. Vastly different Sox17 regulons found in endoderm versus endothelial cells support the differential interactions of SOX17 with other regulatory factors thereby enabling lineage-specific regulatory actions.

ZFP92, a KRAB domain zinc finger protein enriched in pancreatic islets, binds to B1/Alu SINE transposable elements and regulates retroelements and genes

Repressive KRAB domain-containing zinc-finger proteins (KRAB-ZFPs) are abundant in mammalian genomes and contribute both to the silencing of transposable elements (TEs) and to the regulation of developmental stage- and cell type-specific gene expression. Here we describe studies of zinc finger protein 92 (Zfp92), an X-linked KRAB-ZFP that is highly expressed in pancreatic islets of adult mice, by analyzing global Zfp92 knockout (KO) mice. Physiological, transcriptomic and genome-wide chromatin binding studies indicate that the principal function of ZFP92 in mice is to bind to and suppress the activity of B1/Alu type of SINE elements and modulate the activity of surrounding genomic entities. Deletion of Zfp92 leads to changes in expression of select LINE and LTR retroelements and genes located in the vicinity of ZFP92-bound chromatin. The absence of Zfp92 leads to altered expression of specific genes in islets, adipose and muscle that result in modest sex-specific alterations in blood glucose homeostasis, body mass and fat accumulation. In islets, Zfp92 influences blood glucose concentration in postnatal mice via transcriptional effects on Mafb, whereas in adipose and muscle, it regulates Acacb, a rate-limiting enzyme in fatty acid metabolism. In the absence of Zfp92, a novel TE-Capn11 fusion transcript is overexpressed in islets and several other tissues due to de-repression of an IAPez TE adjacent to ZFP92-bound SINE elements in intron 3 of the Capn11 gene. Together, these studies show that ZFP92 functions both to repress specific TEs and to regulate the transcription of specific genes in discrete tissues.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease

A major fraction of loci identified by genome-wide association studies (GWASs) lead to alterations in alternative splicing, but interpretation of how such alterations impact proteins is hindered by the technical limitations of short-read RNA-seq, which cannot directly link splicing events to full-length transcript or protein isoforms. Long-read RNA-seq represents a powerful tool to define and quantify transcript isoforms, and recently, infer protein isoform existence. Here we present a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H 4 PP ≥ 0.75). We generated deep coverage PacBio long-read RNA-seq data (N=∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were novel. By casting the colocalized sQTLs directly onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Using these data, we created one of the first proteome-scale resources defining full-length isoforms impacted by colocalized sQTLs. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense mediated decay (NMD) and 190 that potentially resulted in the expression of new protein isoforms. Finally, we identified colocalizing sQTLs in TPM2 for splice junctions between two mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization. We expect our approach to be widely generalizable across diverse clinical traits and accelerate system-scale analyses of protein isoform activities modulated by GWAS loci.

Genome-Wide Analysis on Transcriptome and Methylome in Prevention of Mammary Tumor Induced by Early Life Combined Botanicals

Breast cancer (BC) is the most common malignancy and the second leading cause of cancer death among women in the United States. The consumption of natural dietary components such as broccoli sprouts (BSp) and green tea polyphenols (GTPs) has demonstrated exciting potential in reducing the risk of BC through the regulation of epigenetic mechanisms. However, little is known about their impacts on reversing epigenomic aberrations that are centrally involved in the initiation and progression of BC. Previously, we have determined the efficacy of combined BSp and GTPs treatment on the inhibition of the growth of a mammary tumor in a transgenic Her2/neu mouse model. We sought to extend our previous study to identify universal biomarkers that represent common mechanistic changes among different mouse models in response to this dietary regime by including a new transgenic mouse model, C3(1)-SV40 TAg (SV40). As a result, we identified novel target genes that were differentially expressed and methylated in response to dietary botanicals when administered singly (BSp and GTPs) and in combination (BSp + GTPs) in both mouse models. We discovered more differentially expressed and methylated genes in the combination treatment group compared to the singly administered groups. Subsequently, several biological pathways related to epigenetic regulations were identified in response to the combination treatment. Furthermore, when compared to the BSp and GTPs treatment alone, the combinatorial treatment showed a more significant impact on the regulation of the epigenetic modifier activities involved in DNA methylation and histone modifications. Our study provides key insights about the impact of the combined administration of BSp and GTPs on BC using a multi-omics analysis, suggesting a combinatorial approach is more efficacious in preventing and inhibiting BC by impacting key tumor-related genes at transcriptomic and methylomic levels. Our findings could be further extrapolated as a comprehensive source for understanding the epigenetic modifications that are associated with the effects of these dietary botanicals on BC prevention.

Identification of Transcription Factor Networks during Mouse Hindlimb Development

Mammalian hindlimb development involves a variety of cells and the regulation of spatiotemporal molecular events, but regulatory networks of transcription factors contributing to hindlimb morphogenesis are not well understood. Here, we identified transcription factor networks during mouse hindlimb morphology establishment through transcriptome analysis. We used four stages of embryonic hindlimb transcription profiles acquired from the Gene Expression Omnibus database (GSE30138), including E10.5, E11.5, E12.5 and E13.5, to construct a gene network using Weighted Gene Co-expression Network Analysis (WGCNA), and defined seven stage-associated modules. After filtering 7625 hub genes, we further prioritized 555 transcription factors with AnimalTFDB3.0. Gene ontology enrichment showed that transcription factors of different modules were enriched in muscle tissue development, connective tissue development, embryonic organ development, skeletal system morphogenesis, pattern specification process and urogenital system development separately. Six regulatory networks were constructed with key transcription factors, which contribute to the development of different tissues. Knockdown of four transcription factors from regulatory networks, including Sox9, Twist1, Snai2 and Klf4, showed that the expression of limb-development-related genes was also inhibited, which indicated the crucial role of transcription factor networks in hindlimb development.

Deletion of Jazf1 gene causes early growth retardation and insulin resistance in mice

Single-nucleotide polymorphisms in the human juxtaposed with another zinc finger protein 1 (JAZF1) gene have repeatedly been associated with both type 2 diabetes (T2D) and height in multiple genome-wide association studies (GWAS); however, the mechanism by which JAZF1 causes these traits is not yet known. To investigate the possible functional role of JAZF1 in growth and glucose metabolism in vivo, we generated Jazf1 knockout (KO) mice and examined body composition and insulin sensitivity both in young and adult mice by using ¹H-nuclear magnetic resonance and hyperinsulinemic-euglycemic clamp techniques. Plasma concentrations of insulin-like growth factor 1 (IGF-1) were reduced in both young and adult Jazf1 KO mice, and young Jazf1 KO mice were shorter in stature than age-matched wild-type mice. Young Jazf1 KO mice manifested reduced fat mass, whereas adult Jazf1 KO mice manifested increased fat mass and reductions in lean body mass associated with increased plasma growth hormone (GH) concentrations. Adult Jazf1 KO manifested muscle insulin resistance that was further exacerbated by high-fat diet feeding. Gene set enrichment analysis in Jazf1 KO liver identified the hepatocyte hepatic nuclear factor 4 alpha (HNF4α), which was decreased in Jazf1 KO liver and in JAZF1 knockdown cells. Moreover, GH-induced IGF-1 expression was inhibited by JAZF1 knockdown in human hepatocytes. Taken together these results demonstrate that reduction of JAZF1 leads to early growth retardation and late onset insulin resistance in vivo which may be mediated through alterations in the GH-IGF-1 axis and HNF4α.

Mechanisms Underlying the Expansion and Functional Maturation of β-Cells in Newborns: Impact of the Nutritional Environment

The functional maturation of insulin-secreting β-cells is initiated before birth and is completed in early postnatal life. This process has a critical impact on the acquisition of an adequate functional β-cell mass and on the capacity to meet and adapt to insulin needs later in life. Many cellular pathways playing a role in postnatal β-cell development have already been identified. However, single-cell transcriptomic and proteomic analyses continue to reveal new players contributing to the acquisition of β-cell identity. In this review, we provide an updated picture of the mechanisms governing postnatal β-cell mass expansion and the transition of insulin-secreting cells from an immature to a mature state. We then highlight the contribution of the environment to β-cell maturation and discuss the adverse impact of an in utero and neonatal environment characterized by calorie and fat overload or by protein deficiency and undernutrition. Inappropriate nutrition early in life constitutes a risk factor for developing diabetes in adulthood and can affect the β-cells of the offspring over two generations. A better understanding of these events occurring in the neonatal period will help developing better strategies to produce functional β-cells and to design novel therapeutic approaches for the prevention and treatment of diabetes.

Insm1, Neurod1, and Pax6 promote murine pancreatic endocrine cell development through overlapping yet distinct RNA transcription and splicing programs

Insm1, Neurod1, and Pax6 are essential for the formation and function of pancreatic endocrine cells. Here, we report comparative immunohistochemical, transcriptomic, functional enrichment, and RNA splicing analyses of these genes using gene knock-out mice. Quantitative immunohistochemical analysis confirmed that elimination of each of these three factors variably impairs the proliferation, survival, and differentiation of endocrine cells. Transcriptomic analysis revealed that each factor contributes uniquely to the transcriptome although their effects were overlapping. Functional enrichment analysis revealed that genes downregulated by the elimination of Insm1, Neurod1, and Pax6 are commonly involved in mRNA metabolism, chromatin organization, secretion, and cell cycle regulation, and upregulated genes are associated with protein degradation, autophagy, and apoptotic process. Elimination of Insm1, Neurod1, and Pax6 impaired expression of many RNA-binding proteins thereby altering RNA splicing events, including for Syt14 and Snap25, two genes required for insulin secretion. All three factors are necessary for normal splicing of Syt14, and both Insm1 and Pax6 are necessary for the processing of Snap25. Collectively, these data provide new insights into how Insm1, Neurod1, and Pax6 contribute to the formation of functional pancreatic endocrine cells.