Integration of ATAC-seq and RNA-seq identifies human alpha cell and beta cell signature genes

Mulberry twig (Sangzhi) alkaloids induce pancreatic α-to-β-cell phenotypic conversion in type 2 diabetic mice

BACKGROUND

Pancreatic β-cell regeneration represents a promising therapeutic strategy for diabetes, yet safe and effective treatments remain elusive.

PURPOSE

We aimed to investigate whether and how mulberry twig (Sangzhi) alkaloids (SZ-A), a newly approved anti-diabetic Chinese medicine, promoted β-cell regeneration.

METHODS

Diabetic db/db mice and pancreatic α-cell lineage-tracing mice were treated with SZ-A, acarbose, or vehicle daily via intragastric gavage. Blood glucose and plasma insulin levels were measured. The areas of islets, α-cells and β-cells were quantified. Cell transdifferentiation was assessed by double-immunostaining of glucagon or α-cell lineage-tracing marker with β-cell-specific markers. Mouse α-cells were incubated with SZ-A or its three main components, and the mRNA levels of cell transdifferentiation-related genes were detected. RNA-sequencing was performed to screen potential targets. The activities of five mitochondrial complexes were detected following treatment, and specific inhibitor was utilized to validate the involvement.

RESULTS

Both SZ-A and acarbose improved glycemic control, but only SZ-A enlarged islet and β-cell areas in the diabetic mice. SZ-A induced α-to-β-cell conversion, as indicated by glucagon and insulin double-immunostaining and confirmed by α-cell lineage-tracing. In cultured α-cells, SZ-A and its main component 1-deoxynojirimycin (DNJ) downregulated the expressions of α-cell-specific markers, while upregulated the expressions of β-cell-specific markers. DNJ-induced differentially expressed genes were enriched in the mitochondrial protein complex term. DNJ inhibited mitochondrial complex I activity, and the complex inhibitor induced α-to-β-cell conversion.

CONCLUSION

SZ-A, especially its main component DNJ, induces α-to-β-cell transdifferentiation via inhibiting mitochondrial complex I. Our finding provides a potential strategy for β-cell regeneration and diabetes treatment.

Exploring the functional, protective, and transcriptomic effects of GIP on cytokine-exposed human pancreatic islets and EndoC-βH5 cells

Immune-mediated beta-cell destruction and lack of alpha-cell responsiveness to hypoglycaemia are hallmarks of type 1 diabetes pathology. The incretin hormone glucose-dependent insulinotropic polypeptide (GIP) may hold therapeutic potential for type 1 diabetes due to its insulinotropic and glucagonotropic effects, as well as its cytoprotective effects shown in rodent beta cells. To further increase our understanding of GIP's effects on human beta cells, we here examined the functional, protective, and transcriptomic effects of GIP in human EndoC-βH5 beta cells and isolated human islets in the presence or absence of proinflammatory cytokines (interferon (IFN)-γ ± interleukin (IL)-1β) as a mimic of type 1 diabetes. GIP dose-dependently augmented glucose-stimulated insulin secretion from EndoC-βH5 cells and increased insulin and glucagon secretion from human islets at high and low glucose concentrations, respectively. The insulinotropic effect of GIP in EndoC-βH5 cells was abrogated by KN-93, an inhibitor of calcium/calmodulin-dependent protein kinase 2 (CaMK2). GIP did not prevent cytokine-induced apoptosis in EndoC-βH5 cells or human islets, and GIP did not protect against cytokine-induced functional impairment in EndoC-βH5 cells. GIP treatment of human islets for 24 h had no effects on the transcriptome and did not modulate cytokine-induced transcriptional changes. However, GIP augmented IL-1β + IFNγ-induced secretion of interleukin (IL)-10 and c-c motif chemokine ligand (CCL)-2 from human islets while decreasing the secretion of c-x-c motif chemokine ligand (CXCL)-8. In EndoC-βH5 cells, GIP reduced IFN-γ-induced secretion of tumor necrosis factor (TNF)-α, IL-2, IL-6, and IL-10 but increased the secretion of CXCL8, CCL2, CCL4, and CCL11. In conclusion, our results suggest that the insulinotropic effect of GIP is CaMK2-dependent. Furthermore, our findings indicate that GIP neither exerts cytoprotective effects against cytokines nor modulate the transcriptome of human islets. GIP may, however, exert selective modulatory effects on secreted inflammatory factors from cytokine-exposed beta cells and islets.

Genome-wide chromatin accessibility and selective signals of meat rabbits reveal key Cis-regulatory elements and variants during postnatal development of skeletal muscles in rabbits

BACKGROUND

The development of skeletal muscles is intricately modulated by multiple genetic factors and significantly impacts the economic value of meat rabbits. However, our knowledge of epigenetics in rabbit skeletal muscles remains largely unknown.

RESULTS

In this study, we collected leg skeletal muscles of rabbits and performed assays for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) to detect open chromatin across three developmental stages: birth (D1), weaning (D35), and adulthood (D75). A total of 126,959 accessible chromatin regions (ACRs) were identified across samples, and a broad increase and decrease in chromatin accessibility were found from D1 to D35 and D35 to D75, respectively. Integrative analysis of chromatin accessibility and transcriptome data revealed ACRs that were nearly closed at D1 but highly accessible at D35 and D75 were significantly enriched in skeletal muscle development. Cis-regulation analysis further revealed that genes dominated by enhancers mainly play roles in the neuron development of rabbit skeletal muscles. Moreover, the detection of selection signals of meat rabbits and the footprinting analysis of transcription factor at open chromatin revealed that both base transversion (Chr13:12144967 A-> G) and the dynamics of chromatin accessibility at the PRDM1 binding site might regulate ZSWIM5 during the development of skeletal muscles in rabbits.

CONCLUSIONS

Our study provided a category of potential cis-regulatory elements for understanding the development of skeletal muscles at the tissue level and might facilitate potential insights into growth regulation in rabbits.

Enhancer transcription profiling reveals an enhancer RNA-driven ferroptosis and new therapeutic opportunities in prostate cancer

Enhancer RNAs (eRNAs), a subclass of non-coding RNAs transcribed from enhancer regions, have emerged as critical regulators of gene expression; however, their functional roles in prostate cancer remain largely unexplored. In this study, we performed integrated chromatin accessibility and transcriptomic analyses using ATAC-seq and RNA-seq on twenty pairs of prostate cancer and matched benign tissues. By incorporating chromatin immunoprecipitation sequencing data, we identified a subset of differentially expressed eRNAs significantly associated with genes involved in prostate development and oncogenic signaling pathways. Among these, lactotransferrin-eRNA (LTFe) was markedly downregulated in prostate cancer tissues, with functional analyses revealing its tumor-suppressive role. Mechanistically, LTFe promotes the transcription of its target gene, lactotransferrin (LTF), by interacting with heterogeneous nuclear ribonucleoprotein F (HNRNPF) and facilitating enhancer-promoter chromatin interactions. Furthermore, we demonstrate that the LTFe-LTF axis facilitates ferroptosis by modulating iron transport. Notably, androgen receptor (AR) signaling disrupts LTFe-associated chromatin looping, leading to ferroptosis resistance. Therapeutically, co- administration of the AR inhibitor enzalutamide and the ferroptosis inducer RSL3 significantly suppressed tumor growth, offering a promising strategy for castration-resistant prostate cancer. Collectively, this study provides novel insights into the mechanistic role of eRNAs in prostate cancer, highlighting the LTFe-LTF axis as a critical epigenetic regulator and potential therapeutic target for improved treatment outcomes.

Unveiling Tissue-Specific RNA Landscapes in Mouse Organs During Fasting and Feeding Using Nanopore Direct RNA Sequencing

Understanding tissue-specific RNA landscapes is essential for uncovering the functional mechanisms of key organs in mammals. However, current knowledge remains limited, as short-read RNA sequencing-the predominant method for assessing gene expression-depends on incomplete gene annotations and struggles to resolve the diverse transcripts produced by genes. To address these limitations, an integrative approach combining nanopore direct RNA sequencing (DRS), ATAC-Seq, and short-read RNA-seq is used. This method enabled the analysis of RNA landscapes across major mouse organs under fasting and fed conditions, representing two extremes of the caloric cycle. This study uncovered tens of thousands of novel transcripts and identified hundreds of genes with tissue-specific expression, revealing additional layers of regulated pathways within each organ that conventional short-read RNA-seq cannot resolve. By profiling transcript expression across multiple organs under identical conditions, it is conducted comparative analyses exposing significant differences in transcript isoforms and regulations. Moreover, nanopore DRS revealed dynamic changes in poly(A) tail length and m6A modifications of transcripts, many regulated in a tissue-specific manner. These changes likely contribute to functional differentiation and metabolic specialization of various organs. Collectively, this findings reveal previously unrecognized layers of gene regulation, offering new insights into the metabolic basis of organ function.

Integrated analysis of ATAC-seq and RNA-seq reveals ADSCP2 regulates oxidative phosphorylation pathway in hypertrophic scar fibroblasts

The primary effector cells involved in the formation of hypertrophic scars are fibroblasts. A potential peptide, ADSCP2 (adipose-derived stem cell peptide 2, the peptide fragment of ALCAM protein), derived from adipose-derived stem cell-conditioned medium, has been identified as having the potential to mitigate hypertrophic scar formation by targeting pyruvate carboxylase. However, the underlying mechanisms remain incompletely understood. Whether ADSCP2 attenuates hypertrophic scar fibrosis at the transcription level remains unclear. Consequently, this study sought to elucidate the potential mechanism associated with ADSCP2 by examining genome-wide transcriptional alterations and changes in chromatin accessibility in fibroblasts. This was achieved through the integrated analysis of assay for transposase accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq). In the ADSCP2 treatment group, ATAC-seq identified a total of 7,805 differential peaks associated with 3,176 genes. RNA-seq analysis revealed 345 upregulated and 399 downregulated transcripts in the same group. A combined Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of both downregulated genes and close-ACRs (accessible chromatin regions) genes within the ADSCP2 treatment group indicated regulation of the oxidative phosphorylation pathway (OXPHOS) by ADSCP2. The amalgamation of ATAC-seq and RNA-seq data elucidates that two OXPHOS associated genes, namely COX6B1 (cytochrome c oxidase subunit 6B1) and NDUFA1 (NADH dehydrogenase (ubiquinone) alpha subcomplex-1), demonstrate significant downregulation in the presence of ADSCP2. Further analysis using the integrative genomics viewer indicates that the promoter regions of both COX6B1 and NDUFA1 exhibit a higher degree of closure in the ADSCP2 treatment group. Quantitative PCR analysis demonstrated that ADSCP2 treatment resulted in a reduction of COX6B1 and NDUFA1 mRNA expression levels. Furthermore, cellular ATP and lactic acid concentrations were diminished in the ADSCP2-treated group. Collectively, these findings suggest potential avenues for future research into the therapeutic application of the peptide ADSCP2 in the treatment of hypertrophic scars.

Mechanistic insights and approaches for beta cell regeneration

Diabetes is characterized by variable loss of insulin-producing beta cells, and new regenerative approaches to increasing the functional beta cell mass of patients hold promise for reversing disease progression. In this Review, we summarize recent chemical biology breakthroughs advancing our knowledge of beta cell regeneration. We present current chemical-based tools, sensors and mechanistic insights into pathways that can be targeted to enhance beta cell regeneration in model organisms. We group the pathways according to the cellular processes they affect, that is, proliferation, conversion of other mature cell types to beta cells and beta cell differentiation from progenitor-like populations. We also suggest assays for assessing the functionality of the regenerated beta cells. Although regeneration processes differ between animal models, such as zebrafish, mice and pigs, regenerative mechanisms identified in any one animal model may be translatable to humans. Overall, chemical biology-based approaches in beta cell regeneration give hope that specific molecular pathways can be targeted to enhance beta cell regeneration.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci, we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18, and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

ATAC-seq for Characterizing Host and Pathogen Genome Accessibility During Virus Infection

Chromatin regulation provides a mechanism through which cells dynamically and rapidly regulate their gene expression profiles, playing a pivotal role in diverse biological processes and disease states. The Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a method that enables genome-wide detection of accessible chromatin regions, providing information on nucleosome positioning and the epigenetic regulation of the chromatin structure. ATAC-seq has been used in various biological contexts, and several reports have demonstrated its application to studying infections with viral or bacterial pathogens. The ability to characterize changes in viral or bacterial genome accessibility during infections provides insights into both pathogen replication and host defense mechanisms. Viral genomes undergo dynamic changes in their structural landscape to facilitate replication and evade host immune responses. Additionally, host cells encode DNA sensors, which are specialized proteins that bind to viral genomes to initiate innate immune responses and sometimes, to suppress viral gene expression. ATAC-seq enables the systematic detection of key structural changes on the viral genome mediated by either viral or host proteins, offering mechanistic insights into virus-host interactions. Here, we describe an ATAC-seq method optimized for studying changes in chromatin accessibility in both host and viral genomes. We have previously applied this method to demonstrate a systematic decrease in the genome accessibility of herpes simplex virus type I (HSV-1) enabled by a host antiviral factor, the interferon-gamma inducible protein 16 (IFI16) during infection of human fibroblasts. This protocol can be adapted to various biological contexts involving the introduction of foreign DNA, making it a valuable tool for a broad range of research endeavors.

G6PC2 controls glucagon secretion by defining the set point for glucose in pancreatic α cells

Elevated glucagon concentrations have been reported in patients with type 2 diabetes (T2D). A critical role for α cell-intrinsic mechanisms in regulating glucagon secretion was previously established through genetic manipulation of the glycolytic enzyme glucokinase (GCK) in mice. Genetic variation at the glucose-6-phosphatase catalytic subunit 2 (G6PC2) locus, encoding an enzyme that opposes GCK, has been reproducibly associated with fasting blood glucose and hemoglobin A1c. Here, we found that trait-associated variants in the G6PC2 promoter are located in open chromatin not just in β but also in α cells and documented allele-specific G6PC2 expression of linked variants in human α cells. Using α cell-specific gene ablation of G6pc2 in mice, we showed that this gene plays a critical role in controlling glucose suppression of amino acid-stimulated glucagon secretion independent of alterations in insulin output, islet hormone content, or islet morphology, findings that we confirmed in primary human α cells. Collectively, our data demonstrate that G6PC2 affects glycemic control via its action in α cells and possibly suggest that G6PC2 inhibitors might help control blood glucose through a bihormonal mechanism.

Dynamic Chromatin Accessibility and Gene Expression Regulation During Maize Leaf Development

BACKGROUND/OBJECTIVES

Chromatin accessibility is closely associated with transcriptional regulation during maize (Zea mays) leaf development. However, its precise role in controlling gene expression at different developmental stages remains poorly understood. This study aimed to investigate the dynamics of chromatin accessibility and its influence on genome-wide gene expression during the BBCH_11, BBCH_13, and BBCH_17 stages of maize leaf development.

METHODS

Maize leaves were collected at the BBCH_11, BBCH_13, and BBCH_17 developmental stages, and chromatin accessibility was assessed using ATAC-seq. RNA-seq was performed to profile gene expression. Integrated analysis of ATAC-seq and RNA-seq data was conducted to elucidate the relationship between chromatin accessibility and transcriptional regulation.

RESULTS

A total of 46,808, 38,242, and 41,084 accessible chromatin regions (ACRs) were identified at the BBCH_11, BBCH_13, and BBCH_17 stages, respectively, with 23.4%, 12.2%, and 21.9% of these regions located near transcription start sites (TSSs). Integrated analyses revealed that both the number and intensity of ACRs significantly influence gene expression levels. Motif analysis identified key transcription factors associated with leaf development and potential transcriptional repressors among genes, showing divergent regulation patterns in ATAC-seq and RNA-seq datasets.

CONCLUSIONS

These findings demonstrate that chromatin accessibility plays a crucial role in regulating the spatial and temporal expression of key genes during maize leaf development by modulating transcription factor binding. This study provides novel insights into the regulatory mechanisms underlying maize leaf development, contributing to a deeper understanding of chromatin-mediated gene expression.

An Unbiased Approach to Identifying Cellular Reprogramming-Inducible Enhancers

Cellular reprogramming of somatic cells towards induced pluripotency is a multistep stochastic process mediated by the transcription factors Oct4, Sox2, Klf4 and c-Myc (OSKM), which orchestrate global epigenetic and transcriptional changes. We performed a large-scale analysis of integrated ChIP-seq, ATAC-seq and RNA-seq data and revealed the spatiotemporal highly dynamic pattern of OSKM DNA binding during reprogramming. We found that OSKM show distinct temporal patterns of binding to different classes of pluripotency-related enhancers. Genes involved in reprogramming are regulated by the coordinated activity of multiple enhancers, which are sequentially bound by OSKM for strict transcriptional control. Based on these findings, we developed an unbiased approach to identify Reprogramming-Inducible Enhancers (RIEs), constructed enhancer-traps and isolated cells undergoing reprogramming in real time. We used a representative RIE taken from the Upp1 gene fused to Gfp and isolated cells at different time-points during reprogramming and found that they have unique developmental capacities as they are reprogrammed with high efficiency due to their distinct molecular signatures. In conclusion, our experiments have led to the development of an unbiased method to identify and isolate reprogrammable cells in real time by exploiting the functional dynamics of OSKM, which can be used as efficient reprogramming biomarkers.

Role of the Egr2 Promoter Antisense RNA in Modulating the Schwann Cell Chromatin Landscape

Background: Schwann cells (SCs) and their plasticity contribute to the peripheral nervous system's capacity for nerve regeneration after injury. The Egr2/Krox20 promoter antisense RNA (Egr2-AS) recruits chromatin remodeling complexes to inhibit Egr2 transcription following peripheral nerve injury. Methods: RNA-seq and ATAC-seq were performed on control cells, Lenti-GFP-transduced cells, and cells overexpressing Egr2-AS (Lenti-AS). Egr2 AS-RNA was cloned into the pLVX-DsRed-Express2-N1 lentiviral expression vector (Clontech, Mountain View, CA, USA), and the levels of AS-RNA expression were determined. Ezh2 and Wdr5 were immunoprecipitated from rat SCs and RT-qPCR was performed against AS-Egr2 RNA. ChIP followed by DNA purification columns was used to perform qPCR for relevant promoters. Hi-C, HiC-DC+, R, Bioconductor, and TOBIAS were used for significant and differential loop analysis, identifications of COREs and CORE-promotor loops, comparisons of TF activity at promoter sites, and identification of site-specific TF footprints. OnTAD was used to detect TADs, and Juicer was used to identify A/B compartments. Results: Here we show that a Neuregulin-ErbB2/3 signaling axis mediates binding of the Egr2-AS to YY1Ser184 and regulates its expression. Egr2-AS modulates the chromatin accessibility of Schwann cells and interacts with two distinct histone modification complexes. It binds to EZH2 and WDR5 and enables targeting of H3K27me3 and H3K4me3 to promoters of Egr2 and C-JUN, respectively. Expression of the Egr2-AS results in reorganization of the global chromatin landscape and quantitative changes in the loop formation and contact frequency at domain boundaries exhibiting enrichment for AP-1 genes. In addition, the Egr2-AS induces changes in the hierarchical TADs and increases transcription factor binding scores on an inter-TAD loop between a super-enhancer regulatory hub and the promoter of mTOR. Conclusions: Our results show that Neuregulin-ErbB2/3-YY1 regulates the expression of Egr2-AS, which mediates remodeling of the chromatin landscape in Schwann cells.

LncRNA Snhg3 aggravates hepatic steatosis via PPARγ signaling

LncRNAs are involved in modulating the individual risk and the severity of progression in metabolic dysfunction-associated fatty liver disease (MASLD), but their precise roles remain largely unknown. This study aimed to investigate the role of lncRNA Snhg3 in the development and progression of MASLD, along with the underlying mechanisms. The result showed that Snhg3 was significantly downregulated in the liver of high-fat diet-induced obesity (DIO) mice. Notably, palmitic acid promoted the expression of Snhg3 and overexpression of Snhg3 increased lipid accumulation in primary hepatocytes. Furthermore, hepatocyte-specific Snhg3 deficiency decreased body and liver weight, alleviated hepatic steatosis and promoted hepatic fatty acid metabolism in DIO mice, whereas overexpression induced the opposite effect. Mechanistically, Snhg3 promoted the expression, stability and nuclear localization of SND1 protein via interacting with SND1, thereby inducing K63-linked ubiquitination modification of SND1. Moreover, Snhg3 decreased the H3K27me3 level and induced SND1-mediated chromatin loose remodeling, thus reducing H3K27me3 enrichment at the Pparg promoter and enhancing PPARγ expression. The administration of PPARγ antagonist T0070907 improved Snhg3-aggravated hepatic steatosis. Our study revealed a new signaling pathway, Snhg3/SND1/H3K27me3/PPARγ, responsible for mice MASLD and indicates that lncRNA-mediated epigenetic modification has a crucial role in the pathology of MASLD.

RUNX1 interacts with lncRNA SMANTIS to regulate monocytic cell functions

Monocytes, the circulating macrophage precursors, contribute to diseases like atherosclerosis and asthma. Long non-coding RNAs (lncRNAs) have been shown to modulate the phenotype and inflammatory capacity of monocytes. We previously discovered the lncRNA SMANTIS, which contributes to cellular phenotype expression by controlling BRG1 in mesenchymal cells. Here, we report that SMANTIS is particularly highly expressed in monocytes and lost during differentiation into macrophages. Moreover, different types of myeloid leukemia presented specific SMANTIS expression patterns. Interaction studies revealed that SMANTIS binds RUNX1, a transcription factor frequently mutated in AML, primarily through its Alu-element on the RUNT domain. RNA-seq after CRISPR/Cas9-mediated deletion of SMANTIS or RUNX1 revealed an association with cell adhesion and both limited the monocyte adhesion to endothelial cells. Mechanistically, SMANTIS KO reduced RUNX1 genomic binding and altered the interaction of RUNX1 with EP300 and CBFB. Collectively, SMANTIS interacts with RUNX1 and attenuates monocyte adhesion, which might limit monocyte vascular egress.

Implications of noncoding regulatory functions in the development of insulinomas

Insulinomas are rare neuroendocrine tumors arising from pancreatic β cells, characterized by aberrant proliferation and altered insulin secretion, leading to glucose homeostasis failure. With the aim of uncovering the role of noncoding regulatory regions and their aberrations in the development of these tumors, we coupled epigenetic and transcriptome profiling with whole-genome sequencing. As a result, we unraveled somatic mutations associated with changes in regulatory functions. Critically, these regions impact insulin secretion, tumor development, and epigenetic modifying genes, including polycomb complex components. Chromatin remodeling is apparent in insulinoma-selective domains shared across patients, containing a specific set of regulatory sequences dominated by the SOX17 binding motif. Moreover, many of these regions are H3K27me3 repressed in β cells, suggesting that tumoral transition involves derepression of polycomb-targeted domains. Our work provides a compendium of aberrant cis-regulatory elements affecting the function and fate of β cells in their progression to insulinomas and a framework to identify coding and noncoding driver mutations.

Short-term cold exposure induces persistent epigenomic memory in brown fat

Deficiency of the epigenome modulator histone deacetylase 3 (HDAC3) in brown adipose tissue (BAT) impairs the ability of mice to survive in near-freezing temperatures. Here, we report that short-term exposure to mild cold temperature (STEMCT: 15°C for 24 h) averted lethal hypothermia of mice lacking HDAC3 in BAT (HDAC3 BAT KO) exposed to 4°C. STEMCT restored the induction of the thermogenic coactivator PGC-1α along with UCP1 at 22°C, which is greatly impaired in HDAC3-deficient BAT, and deletion of either UCP1 or PGC-1α prevented the protective effect of STEMCT. Remarkably, this protection lasted for up to 7 days. Transcriptional activator C/EBPβ was induced by short-term cold exposure in mouse and human BAT and, uniquely, remained high for 7 days following STEMCT. Adeno-associated virus-mediated knockdown of BAT C/EBPβ in HDAC3 BAT KO mice erased the persistent memory of STEMCT, revealing the existence of a C/EBPβ-dependent and HDAC3-independent cold-adaptive epigenomic memory.

The chromatin accessibility and transcriptomic landscape of the aging mice cochlea and the identification of potential functional super-enhancers in age-related hearing loss

BACKGROUND

Presbycusis, also referred to as age-related hearing loss (ARHL), is a condition that results from the cumulative effects of aging on an individual's auditory capabilities. Given the limited understanding of epigenetic mechanisms in ARHL, our research focuses on alterations in chromatin-accessible regions.

METHODS

We employed assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) in conjunction with unique identifier (UID) mRNA-seq between young and aging cochleae, and conducted integrated analysis as well as motif/TF-gene prediction. Additionally, the essential role of super-enhancers (SEs) in the development of ARHL was identified by comparative analysis to previous research. Meanwhile, an ARHL mouse model and an aging mimic hair cell (HC) model were established with a comprehensive identification of senescence phenotypes to access the role of SEs in ARHL progression.

RESULTS

The control cochlear tissue exhibited greater chromatin accessibility than cochlear tissue affected by ARHL. Furthermore, the levels of histone 3 lysine 27 acetylation were significantly depressed in both aging cochlea and aging mimic HEI-OC1 cells, highlighting the essential role of SEs in the development of ARHL. The potential senescence-associated super-enhancers (SASEs) of ARHL were identified, most of which exhibited decreased chromatin accessibility. The majority of genes related to the SASEs showed obvious decreases in mRNA expression level in aging HCs and was noticeably altered following treatment with JQ1 (a commonly used SE inhibitor).

CONCLUSION

The chromatin accessibility in control cochlear tissue was higher than that in cochlear tissue affected by ARHL. Potential SEs involved in ARHL were identified, which might provide a basis for future therapeutics targeting SASEs related to ARHL.

Ca2+ signaling and metabolic stress-induced pancreatic β-cell failure

Early in the development of Type 2 diabetes (T2D), metabolic stress brought on by insulin resistance and nutrient overload causes β-cell hyperstimulation. Herein we summarize recent studies that have explored the premise that an increase in the intracellular Ca2+ concentration ([Ca2+]i), brought on by persistent metabolic stimulation of β-cells, causes β-cell dysfunction and failure by adversely affecting β-cell function, structure, and identity. This mini-review builds on several recent reviews that also describe how excess [Ca2+]i impairs β-cell function.

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

Update on the transdifferentiation of pancreatic cells into functional beta cells for treating diabetes

The increasing global prevalence and associated comorbidities need innovative approaches for type 2 diabetes mellitus (T2DM) prevention and treatment. Genetics contributes significantly to T2DM susceptibility, and genetic counseling is significant in detecting and informing people about the diabetic risk. T2DM is also intricately linked to overnutrition and obesity, and nutritional advising is beneficial to mitigate diabetic evolution. However, manipulating pancreatic cell plasticity and transdifferentiation could help beta cell regeneration and glucose homeostasis, effectively contributing to the antidiabetic fight. Targeted modulation of transcription factors is highlighted for their roles in various aspects of pancreatic cell differentiation and function, inducing non-beta cells' conversion into functional beta cells (responsive to glucose). In addition, pharmacological interventions targeting specific receptors and pathways might facilitate cell transdifferentiation aiming to maintain or increase beta cell mass and function. However, the mechanisms underlying cellular reprogramming are not yet well understood. The present review highlights the primary transcriptional factors in the endocrine pancreas, focusing on transdifferentiation as a primary mechanism. Therefore, islet cell reprogramming, converting one cell type to another and transforming non-beta cells into insulin-producing cells, depends, among others, on transcription factors. It is a promising fact that new transcription factors are discovered every day, and their actions on pancreatic islet cells are revealed. Exploring these pathways associated with pancreatic development and islet endocrine cell differentiation could unravel the molecular intricacies underlying transdifferentiation processes, exploring novel therapeutic strategies to treat diabetes. The medical use of this biotechnology is expected to be achievable within a short time.

Contribution of CENP-F to FOXM1-Mediated Discordant Centromere and Kinetochore Transcriptional Regulation

Proper chromosome segregation is required to ensure chromosomal stability. The centromere (CEN) is a unique chromatin domain defined by CENP-A and is responsible for recruiting the kinetochore (KT) during mitosis, ultimately regulating microtubule spindle attachment and mitotic checkpoint function. Upregulation of many CEN/KT genes is commonly observed in cancer. Here, we show that although FOXM1 occupies promoters of many CEN/KT genes with MYBL2, FOXM1 overexpression alone is insufficient to drive the FOXM1-correlated transcriptional program. CENP-F is canonically an outer kinetochore component; however, it functions with FOXM1 to coregulate G2/M transcription and proper chromosome segregation. Loss of CENP-F results in altered chromatin accessibility at G2/M genes and reduced FOXM1-MBB complex formation. We show that coordinated CENP-FFOXM1 transcriptional regulation is a cancer-specific function. We observe a small subset of CEN/KT genes including CENP-C, that are not regulated by FOXM1. Upregulation of CENP-C in the context of CENP-A overexpression leads to increased chromosome missegregation and cell death suggesting that escape of CENP-C from FOXM1 regulation is a cancer survival mechanism. Together, we show that FOXM1 and CENP-F coordinately regulate G2/M genes, and this coordination is specific to a subset of genes to allow for maintenance of chromosome instability levels and subsequent cell survival.

Genome-wide profiling of angiogenic cis-regulatory elements unravels cis-regulatory SNPs for vascular abnormality

Angiogenesis is extensively involved in embryonic development and requires complex regulation networks, whose defects can cause a variety of vascular abnormalities. Cis-regulatory elements control gene expression at all developmental stages, but they have not been studied or profiled in angiogenesis yet. In this study, we exploited public DNase-seq and RNA-seq datasets from a VEGFA-stimulated in vitro angiogenic model, and carried out an integrated analysis of the transcriptome and chromatin accessibility across the entire process. Totally, we generated a bank of 47,125 angiogenic cis-regulatory elements with promoter (marker by H3K4me3) and/or enhancer (marker by H3K27ac) activities. Motif enrichment analysis revealed that these angiogenic cis-regulatory elements interacted preferentially with ETS family TFs. With this tool, we performed an association study using our WES data of TAPVC and identified rs199530718 as a cis-regulatory SNP associated with disease risk. Altogether, this study generated a genome-wide bank of angiogenic cis-regulatory elements and illustrated its utility in identifying novel cis-regulatory SNPs for TAPVC, expanding new horizons of angiogenesis as well as vascular abnormality genetics.

Chromatin accessibility profiling reveals that human fibroblasts respond to mechanical stimulation in a cell-specific manner

Fibroblasts in the skin are highly heterogeneous, both in vivo and in vitro. One difference between follicular (dermal papilla fibroblasts [DP]) and interfollicular fibroblasts (papillary fibroblasts [PFi]) in vitro is their ability to differentiate in response to osteogenic media (OM), or mechanical stimulation. Here, we asked whether differences in the ability of DP and PFi to respond to differentiation stimuli are due to differences in chromatin accessibility. We performed chromatin accessibility and transcriptional profiling of DP and PFi in human skin, which arise from a common progenitor during development, yet display distinct characteristics in adult tissue and in vitro. We found that cells cultured in growth media had unique chromatin accessibility profiles; however, these profiles control similar functional networks. Upon introduction of a chemical perturbation (OM) to promote differentiation, we observed a divergence not only in the accessible chromatin signatures but also in the functional networks controlled by these signatures. The biggest divergence between DP and PFi was observed when we applied 2 perturbations to cells: growth in OM and mechanical stimulation (a shock wave [OMSW]). DP readily differentiate into bone in OMSW conditions, while PFi lack differentiation capability in vitro. In the DP we found a number of uniquely accessible promoters that controlled osteogenic interaction networks associated with bone and differentiation functions. Using ATAC-seq and RNA-seq we found that the combination of 2 stimuli (OMSW) could result in significant changes in chromatin accessibility associated with osteogenic differentiation, but only within the DP (capable of osteogenic differentiation). De novo motif analysis identified enrichment of motifs bound by the TEA domain (TEAD) family of transcription factors, and inter-cell comparisons (UpSet analysis) displayed large groups of genes to be unique to single cell types and conditions. Our results suggest that these 2 stimuli (OMSW) elicit cell-specific responses by modifying chromatin accessibility of osteogenic-related gene promoters.

Nuclear receptor corepressors non-canonically drive glucocorticoid receptor-dependent activation of hepatic gluconeogenesis

Nuclear receptor corepressors (NCoRs) function in multiprotein complexes containing histone deacetylase 3 (HDAC3) to alter transcriptional output primarily through repressive chromatin remodelling at target loci1-5. In the liver, loss of HDAC3 causes a marked hepatosteatosis largely because of de-repression of genes involved in lipid metabolism6,7; however, the individual roles and contribution of other complex members to hepatic and systemic metabolic regulation are unclear. Here we show that adult loss of both NCoR1 and NCoR2 (double knockout (KO)) in hepatocytes phenocopied the hepatomegalic fatty liver phenotype of HDAC3 KO. In addition, double KO livers exhibited a dramatic reduction in glycogen storage and gluconeogenic gene expression that was not observed with hepatic KO of individual NCoRs or HDAC3, resulting in profound fasting hypoglycaemia. This surprising HDAC3-independent activation function of NCoR1 and NCoR2 is due to an unexpected loss of chromatin accessibility on deletion of NCoRs that prevented glucocorticoid receptor binding and stimulatory effect on gluconeogenic genes. These studies reveal an unanticipated, non-canonical activation function of NCoRs that is required for metabolic health.

Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

Iterative assay for transposase-accessible chromatin by sequencing to isolate functionally relevant neuronal subtypes

The Drosophila brain contains tens of thousands of distinct cell types. Thousands of different transgenic lines reproducibly target specific neuron subsets, yet most still express in several cell types. Furthermore, most lines were developed without a priori knowledge of where the transgenes would be expressed. To aid in the development of cell type-specific tools for neuronal identification and manipulation, we developed an iterative assay for transposase-accessible chromatin (ATAC) approach. Open chromatin regions (OCRs) enriched in neurons, compared to whole bodies, drove transgene expression preferentially in subsets of neurons. A second round of ATAC-seq from these specific neuron subsets revealed additional enriched OCR2s that further restricted transgene expression within the chosen neuron subset. This approach allows for continued refinement of transgene expression, and we used it to identify neurons relevant for sleep behavior. Furthermore, this approach is widely applicable to other cell types and to other organisms.

Glucagon-producing α-cell transcriptional identity and reprogramming towards insulin production

β-Cell replacement by in situ reprogramming of non-β-cells is a promising diabetes therapy. Following the observation that near-total β-cell ablation in adult mice triggers the reprogramming of pancreatic α-, δ-, and γ-cells into insulin (INS)-producing cells, recent studies are delving deep into the mechanisms controlling adult α-cell identity. Systematic analyses of the α-cell transcriptome and epigenome have started to pinpoint features that could be crucial for maintaining α-cell identity. Using different transgenic and chemical approaches, significant advances have been made in reprogramming α-cells in vivo into INS-secreting cells in mice. The recent reprogramming of human α-cells in vitro is an important step forward that must now be complemented with a comprehensive molecular dissection of the mechanisms controlling α-cell identity.

Eomes-dependent mitochondrial regulation promotes survival of pathogenic CD4+ T cells during inflammation

The mechanisms whereby Eomes controls tissue accumulation of T cells and strengthens inflammation remain ill-defined. Here, we show that Eomes deletion in antigen-specific CD4+ T cells is sufficient to protect against central nervous system (CNS) inflammation. While Eomes is dispensable for the initial priming of CD4+ T cells, it is required for long-term maintenance of CNS-infiltrating CD4+ T cells. We reveal that the impact of Eomes on effector CD4+ T cell longevity is associated with sustained expression of multiple genes involved in mitochondrial organization and functions. Accordingly, epigenetic studies demonstrate that Eomes supports mitochondrial function by direct binding to either metabolism-associated genes or mitochondrial transcriptional modulators. Besides, the significance of these findings was confirmed in CD4+ T cells from healthy donors and multiple sclerosis patients. Together, our data reveal a new mechanism by which Eomes promotes severity and chronicity of inflammation via the enhancement of CD4+ T cell mitochondrial functions and resistance to stress-induced cell death.

TDP-43 chronic deficiency leads to dysregulation of transposable elements and gene expression by affecting R-loop and 5hmC crosstalk

TDP-43 is an RNA/DNA-binding protein that forms aggregates in various brain disorders. TDP-43 engages in many aspects of RNA metabolism, but its molecular roles in regulating genes and transposable elements (TEs) have not been extensively explored. Chronic TDP-43 knockdown impairs cell proliferation and cellular responses to DNA damage. At the molecular level, TDP-43 chronic deficiency affects gene expression either locally or distally by concomitantly altering the crosstalk between R-loops and 5-hydroxymethylcytosine (5hmC) in gene bodies and long-range enhancer/promoter interactions. Furthermore, TDP-43 knockdown induces substantial disease-relevant TE activation by influencing their R-loop and 5hmC homeostasis in a locus-specific manner. Together, our findings highlight the genomic roles of TDP-43 in modulating R-loop-5hmC coordination in coding genes, distal regulatory elements, and TEs, presenting a general and broad molecular mechanism underlying the contributions of proteinopathies to the etiology of neurodegenerative disorders.

Loss of the accessory chromosome converts a pathogenic tree-root fungus into a mutualistic endophyte

Some fungal accessory chromosomes (ACs) may contribute to virulence in plants. However, the mechanisms by which ACs determine specific traits associated with lifestyle transitions along a symbiotic continuum are not clear. Here we delineated the genetic divergence in two sympatric but considerably variable isolates (16B and 16W) of the poplar-associated fungus Stagonosporopsis rhizophilae. We identified a ∼0.6-Mb horizontally acquired AC in 16W that resulted in a mildly parasitic lifestyle in plants. Complete deletion of the AC (Δ16W) significantly altered the fungal phenotype. Specifically, Δ16W was morphologically more similar to 16B, showed enhanced melanization, and established beneficial interactions with poplar plants, thereby acting as a dark septate endophyte. RNA sequencing (RNA-seq) analysis showed that AC loss induced the upregulation of genes related to root colonization and biosynthesis of indole acetic acid and melanin. We observed that the AC maintained a more open status of chromatin across the genome, indicating an impressive remodeling of cis-regulatory elements upon AC loss, which potentially enhanced symbiotic effectiveness. We demonstrated that the symbiotic capacities were non-host-specific through comparable experiments on Triticum- and Arabidopsis-fungus associations. Furthermore, the three isolates generated symbiotic interactions with a nonvascular liverwort. In summary, our study suggests that the AC is a suppressor of symbiosis and provides insights into the underlying mechanisms of mutualism with vascular plants in the absence of traits encoded by the AC. We speculate that AC-situated effectors and other potential secreted molecules may have evolved to specifically target vascular plants and promote mild virulence.

Disrupted RNA editing in beta cells mimics early-stage type 1 diabetes

A major hypothesis for the etiology of type 1 diabetes (T1D) postulates initiation by viral infection, leading to double-stranded RNA (dsRNA)-mediated interferon response and inflammation; however, a causal virus has not been identified. Here, we use a mouse model, corroborated with human islet data, to demonstrate that endogenous dsRNA in beta cells can lead to a diabetogenic immune response, thus identifying a virus-independent mechanism for T1D initiation. We found that disruption of the RNA editing enzyme adenosine deaminases acting on RNA (ADAR) in beta cells triggers a massive interferon response, islet inflammation, and beta cell failure and destruction, with features bearing striking similarity to early-stage human T1D. Glycolysis via calcium enhances the interferon response, suggesting an actionable vicious cycle of inflammation and increased beta cell workload.

Chronic consumption of imbalance diets high in sucrose or fat induces abdominal obesity with different pattern of metabolic disturbances and lost in Langerhans cells population

AIM

Obesity is a worldwide health issue, associated with development of type 2 Diabetes Mellitus. The aim of this study is to analyze the effect of consumption of two hypercaloric diets on metabolic disturbance and beta cells damage.

MAIN METHODS

Male Wistar rats were subjected to twelve months consumption of three diets: a Control balanced diet (CTD, carbohydrates 58 %, proteins 29 %, lipids 13 %) and two hypercaloric diets, high in sucrose (HSD, carbohydrates 68 %, proteins 22 %, lipids 10 %) or high in fat (HFD, carbohydrates 31 %, proteins 14 %, lipids 55 %). Serum levels of glucose, triglycerides and free fatty acids were measured after zoometric parameters determination. Antioxidant enzymes activity and oxidative stress-marker were measured in pancreas tissue among histological analysis of Langerhans islets.

KEY FINDINGS

Although diets were hypercaloric, the amount of food consumed by rats decreased, resulting in an equal caloric consumption. The HSD induced hypertriglyceridemia and hyperglycemia with higher levels in free fatty acids (FFA, lipotoxicity); whereas HFD did not increased neither the triglycerides nor FFA, nevertheless the loss of islets' cell was larger. Both diets induced obesity with hyperglycemia and significant reduction in Langerhans islets size.

SIGNIFICANCE

Our results demonstrate that consumption of HSD induces more significant metabolic disturbances that HFD, although both generated pancreas damage; as well hypercaloric diet consumption is not indispensable to becoming obese; the chronic consumption of unbalanced diets (rich in carbohydrates or lipids) may lead to abdominal obesity with metabolic and functional disturbances, although the total amount of calories are similar.

An Insight into Vital Genes Responsible for β-cell Formation

The regulation of glucose homeostasis and insulin secretion by pancreatic β-cells, when disturbed, will result in diabetes mellitus. Replacement of dysfunctional or lost β-cells with fully functional ones can tackle the problem of β-cell generation in diabetes mellitus. Various pancreatic-specific genes are expressed during different stages of development, which have essential roles in pancreatogenesis and β-cell formation. These factors play a critical role in cellular-based studies like transdifferentiation or de-differentiation of somatic cells to multipotent or pluripotent stem cells and their differentiation into functional β-cells. This work gives an overview of crucial transcription factors expressed during various stages of pancreas development and their role in β-cell specification. In addition, it also provides a perspective on the underlying molecular mechanisms.

Bioactive Vitamins and Epigenetic Modifications in Diabetes: A Perspective

Diabetes is a complex metabolic disease that has been associated with epigenetic changes. External factors such as dietary patterns can induce an imbalance in the pools of micronutrients and macronutrients in the body. Consequently, bioactive vitamins may influence epigenetic mechanisms via several pathways: involvement in the control of gene expression, and in protein synthesis, by acting as coenzymes and co-factors in the metabolism of methyl groups or methylation of DNA and histones. Herein, we present a perspective on the relevance of bioactive vitamins in the epigenetic modifications that occur in diabetes.

Contribution of CENP-F to FOXM1-mediated discordant centromere and kinetochore transcriptional regulation

Proper chromosome segregation is required to ensure genomic and chromosomal stability. The centromere is a unique chromatin domain present throughout the cell cycle on each chromosome defined by the CENP-A nucleosome. Centromeres (CEN) are responsible for recruiting the kinetochore (KT) during mitosis, ultimately regulating spindle attachment and mitotic checkpoint function. Upregulation of many genes that encode the CEN/KT proteins is commonly observed in cancer. Here, we show although that FOXM1 occupies the promoters of many CEN/KT genes with MYBL2, occupancy is insufficient alone to drive the FOXM1 correlated transcriptional program. We show that CENP-F, a component of the outer kinetochore, functions with FOXM1 to coregulate G2/M transcription and proper chromosome segregation. Loss of CENP-F results in alteration of chromatin accessibility at G2/M genes, including CENP-A, and leads to reduced FOXM1-MBB complex formation. The FOXM1-CENP-F transcriptional coordination is a cancer-specific function. We observed that a few CEN/KT genes escape FOXM1 regulation such as CENP-C which when upregulated with CENP-A, leads to increased chromosome misegregation and cell death. Together, we show that the FOXM1 and CENP-F coordinately regulate G2/M gene expression, and this coordination is specific to a subset of genes to allow for proliferation and maintenance of chromosome stability for cancer cell survival.

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.

25-Hydroxyvitamin D Is Associated with Islet Homeostasis in Type-2 Diabetic Patients with Abdominal Obesity

OBJECTIVE

Islet α cells input is essential for insulin secretion from β cells. The present study aims to investigate the association between 25-hydroxyvitamin D [25(OH)D] and islet function homeostasis in type-2 diabetes (T2D) patients.

METHODS

A total of 4670 T2D patients from seven communities in Shanghai, China were enrolled. The anthropometric indices, biochemical parameters, serum 25(OH)D, and islet function [including C-peptide (C-p) and glucagon] were measured.

RESULTS

The fasting plasma glucose (FPG), glycated hemoglobin (HbA1c), glucagon, and C-p levels exhibited a significantly decreasing trend in T2D patients as the 25(OH)D levels increased. Next, the population was divided into two groups: abdominal obesity and non-abdominal obesity groups. After adjustment, the 25(OH)D level was found to be associated with HbA1c, glucagon, and homeostasis model assessment of β (HOMA-β) in the non-abdominal obesity group. There was a significant relationship between 25(OH)D and HbA1c, glucagon, HOMA-IR, baseline insulin or C-p in the abdominal obesity group. In the abdominal obesity group, the ordinary least squares (OLS) regression and quantile regression revealed that 25(OH) D was obviously associated with glucagon and fasting C-p levels. In the abdominal obesity group, the moderate analysis revealed a significant interaction effect of 25(OH)D and glucagon on C-p (P=0.0124). Furthermore, the conditional indirect effect of 25(OH)D on the glucagon/C-p ratio was significantly lower at 1 standard deviation (SD) below the mean (P=0.0002), and lower at the mean of the course of diabetes (P=0.0007).

CONCLUSION

25(OH)D was found to be negatively correlated to glucagon and C-p in T2D patients with abdominal obesity. The 25(OH)D influenced C-p in part by influencing glucagon. The effect of 25(OH)D on the glucagon/C-p ratio in T2D patients with abdominal obesity, in terms of islet homeostasis, is influenced by the course of diabetes.

Metabolic regulation of glucagon secretion

Since the discovery of glucagon 100 years ago, the hormone and the pancreatic islet alpha cells that produce it have remained enigmatic relative to insulin-producing beta cells. Canonically, alpha cells have been described in the context of glucagon's role in glucose metabolism in liver, with glucose as the primary nutrient signal regulating alpha cell function. However, current data reveal a more holistic model of metabolic signalling, involving glucagon-regulated metabolism of multiple nutrients by the liver and other tissues, including amino acids and lipids, providing reciprocal feedback to regulate glucagon secretion and even alpha cell mass. Here we describe how various nutrients are sensed, transported and metabolised in alpha cells, providing an integrative model for the metabolic regulation of glucagon secretion and action. Importantly, we discuss where these nutrient-sensing pathways intersect to regulate alpha cell function and highlight key areas for future research.

Not the second fiddle: α cell development, identity, and function in health and diabetes

Historic and emerging studies provide evidence for the deterioration of pancreatic α cell function and identity in diabetes mellitus. Increased access to human tissue and the availability of more sophisticated molecular technologies have identified key insights into how α cell function and identity are preserved in healthy conditions and how they become dysfunctional in response to stress. These studies have revealed evidence of impaired glucagon secretion, shifts in α cell electrophysiology, changes in α cell mass, dysregulation of α cell transcription, and α-to-β cell conversion prior to and during diabetes. In this review, we outline the current state of research on α cell identity in health and disease. Evidence in model organisms and humans suggests that in addition to β cell dysfunction, diabetes is associated with a fundamental dysregulation of α cell identity. Importantly, epigenetic studies have revealed that α cells retain more poised and open chromatin at key cell-specific and diabetes-dysregulated genes, supporting the model that the inherent epigenetic plasticity of α cells makes them susceptible to the transcriptional changes that potentiate the loss of identity and function seen in diabetes. Thus, additional research into the maintenance of α cell identity and function is critical to fully understanding diabetes. Furthermore, these studies suggest α cells could represent an alternative source of new β cells for diabetes treatment.

A genetically supported drug repurposing pipeline for diabetes treatment using electronic health records

BACKGROUND

The identification of new uses for existing drug therapies has the potential to identify treatments for comorbid conditions that have the added benefit of glycemic control while also providing a rapid, low-cost approach to drug (re)discovery.

METHODS

We developed and tested a genetically-informed drug-repurposing pipeline for diabetes management. This approach mapped genetically-predicted gene expression signals from the largest genome-wide association study for type 2 diabetes mellitus to drug targets using publicly available databases to identify drug-gene pairs. These drug-gene pairs were then validated using a two-step approach: 1) a self-controlled case-series (SCCS) using electronic health records from a discovery and replication population, and 2) Mendelian randomization (MR).

FINDINGS

After filtering on sample size, 20 candidate drug-gene pairs were validated and various medications demonstrated evidence of glycemic regulation including two anti-hypertensive classes: angiotensin-converting enzyme inhibitors as well as calcium channel blockers (CCBs). The CCBs demonstrated the strongest evidence of glycemic reduction in both validation approaches (SCCS HbA1c and glucose reduction: -0.11%, p = 0.01 and -0.85 mg/dL, p = 0.02, respectively; MR: OR = 0.84, 95% CI = 0.81, 0.87, p = 5.0 x 10-25).

INTERPRETATION

Our results support CCBs as a strong candidate medication for blood glucose reduction in addition to cardiovascular disease reduction. Further, these results support the adaptation of this approach for use in future drug-repurposing efforts for other conditions.

FUNDING

National Institutes of Health, Medical Research Council Integrative Epidemiology Unit at the University of Bristol, UK Medical Research Council, American Heart Association, and Department of Veterans Affairs (VA) Informatics and Computing Infrastructure and VA Cooperative Studies Program.

C/EBPα Confers Dependence to Fatty Acid Anabolic Pathways and Vulnerability to Lipid Oxidative Stress-Induced Ferroptosis in FLT3-Mutant Leukemia

Although transcription factor CCAAT-enhancer binding protein α (C/EBPα) is critical for normal and leukemic differentiation, its role in cell and metabolic homeostasis is largely unknown in cancer. Here, multiomics analyses uncovered a coordinated activation of C/EBPα and Fms-like tyrosine kinase 3 (FLT3) that increased lipid anabolism in vivo and in patients with FLT3-mutant acute myeloid leukemia (AML). Mechanistically, C/EBPα regulated the fatty acid synthase (FASN)-stearoyl-CoA desaturase (SCD) axis to promote fatty acid (FA) biosynthesis and desaturation. We further demonstrated that FLT3 or C/EBPα inactivation decreased monounsaturated FA incorporation to membrane phospholipids through SCD downregulation. Consequently, SCD inhibition enhanced susceptibility to lipid redox stress that was exploited by combining FLT3 and glutathione peroxidase 4 inhibition to trigger lipid oxidative stress, enhancing ferroptotic death of FLT3-mutant AML cells. Altogether, our study reveals a C/EBPα function in lipid homeostasis and adaptation to redox stress, and a previously unreported vulnerability of FLT3-mutant AML to ferroptosis with promising therapeutic application.

SIGNIFICANCE

FLT3 mutations are found in 30% of AML cases and are actionable by tyrosine kinase inhibitors. Here, we discovered that C/EBPα regulates FA biosynthesis and protection from lipid redox stress downstream mutant-FLT3 signaling, which confers a vulnerability to ferroptosis upon FLT3 inhibition with therapeutic potential in AML. This article is highlighted in the In This Issue feature, p. 1501.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Single-nucleus multi-omics of human stem cell-derived islets identifies deficiencies in lineage specification

Insulin-producing β cells created from human pluripotent stem cells have potential as a therapy for insulin-dependent diabetes, but human pluripotent stem cell-derived islets (SC-islets) still differ from their in vivo counterparts. To better understand the state of cell types within SC-islets and identify lineage specification deficiencies, we used single-nucleus multi-omic sequencing to analyse chromatin accessibility and transcriptional profiles of SC-islets and primary human islets. Here we provide an analysis that enabled the derivation of gene lists and activity for identifying each SC-islet cell type compared with primary islets. Within SC-islets, we found that the difference between β cells and awry enterochromaffin-like cells is a gradient of cell states rather than a stark difference in identity. Furthermore, transplantation of SC-islets in vivo improved cellular identities overtime, while long-term in vitro culture did not. Collectively, our results highlight the importance of chromatin and transcriptional landscapes during islet cell specification and maturation.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Chromatin accessibility differences between alpha, beta, and delta cells identifies common and cell type-specific enhancers

BACKGROUND

High throughput sequencing has enabled the interrogation of the transcriptomic landscape of glucagon-secreting alpha cells, insulin-secreting beta cells, and somatostatin-secreting delta cells. These approaches have furthered our understanding of expression patterns that define healthy or diseased islet cell types and helped explicate some of the intricacies between major islet cell crosstalk and glucose regulation. All three endocrine cell types derive from a common pancreatic progenitor, yet alpha and beta cells have partially opposing functions, and delta cells modulate and control insulin and glucagon release. While gene expression signatures that define and maintain cellular identity have been widely explored, the underlying epigenetic components are incompletely characterized and understood. However, chromatin accessibility and remodeling is a dynamic attribute that plays a critical role to determine and maintain cellular identity.

RESULTS

Here, we compare and contrast the chromatin landscape between mouse alpha, beta, and delta cells using ATAC-Seq to evaluate the significant differences in chromatin accessibility. The similarities and differences in chromatin accessibility between these related islet endocrine cells help define their fate in support of their distinct functional roles. We identify patterns that suggest that both alpha and delta cells are poised, but repressed, from becoming beta-like. We also identify patterns in differentially enriched chromatin that have transcription factor motifs preferentially associated with different regions of the genome. Finally, we not only confirm and visualize previously discovered common endocrine- and cell specific- enhancer regions across differentially enriched chromatin, but identify novel regions as well. We compiled our chromatin accessibility data in a freely accessible database of common endocrine- and cell specific-enhancer regions that can be navigated with minimal bioinformatics expertise.

CONCLUSIONS

Both alpha and delta cells appear poised, but repressed, from becoming beta cells in murine pancreatic islets. These data broadly support earlier findings on the plasticity in identity of non-beta cells under certain circumstances. Furthermore, differential chromatin accessibility shows preferentially enriched distal-intergenic regions in beta cells, when compared to either alpha or delta cells.

β Cell and Autophagy: What Do We Know?

Pancreatic β cells are central to glycemic regulation through insulin production. Studies show autophagy as an essential process in β cell function and fate. Autophagy is a catabolic cellular process that regulates cell homeostasis by recycling surplus or damaged cell components. Impaired autophagy results in β cell loss of function and apoptosis and, as a result, diabetes initiation and progress. It has been shown that in response to endoplasmic reticulum stress, inflammation, and high metabolic demands, autophagy affects β cell function, insulin synthesis, and secretion. This review highlights recent evidence regarding how autophagy can affect β cells' fate in the pathogenesis of diabetes. Furthermore, we discuss the role of important intrinsic and extrinsic autophagy modulators, which can lead to β cell failure.

Shedding new light on the role of ERAP1 in Type 1 diabetes: A perspective on disease management

Type 1 diabetes mellitus (T1D) is a multifactorial organ specific autoimmune disease which originates from the destruction of insulin-producing beta cells within the pancreatic islets by autoreactive CD8+ T lymphocytes. The autoimmune responses are raised against autoantigenic peptides presented in the context of the Major Histocompatibility Complex (MHC) class I molecules. Peptides are generated in the cytoplasm of the beta cell by degradation through the proteasome activity and other proteases. Proteolytic intermediate protein fragments are then vehicled into the endoplasmic reticulum (ER) by transporters associated with antigen processing TAP1 and TAP2. In the ER, Endoplasmic Reticulum Aminopeptidase 1 (ERAP1) and 2 (ERAP2) shape the intermediate proteins to produce the optimal peptide size for loading into the MHC class I molecules. Subsequently complexes are shuttled to the cell surface for antigen presentation. Genome Wide Association Studies (GWAS) have identified different SNPs of ERAP1 associated to several autoimmune diseases and in particular the T1D-related ERAP1 SNP rs30187 encoding for K528R ERAP1. An association between the ER stress and the increased exposure of beta cells to the immune system has been hypothesized to further contribute to the etiopathogenesis. In particular in a recent study by Thomaidou et al. 2020 (doi: https://doi.org/10.2337/db19-0984) the posttranscriptional regulation of ERAP1 is shown to shaping the recognition of the preproinsulin (PPI) signal peptide by cytotoxic T lymphocytes. In the light of foregoing ERAP1 inhibitors could potentially prevent the activation of epitope-specific autoimmune-promoting T cells and their cytokine production; further regulating ERAP1 expression at posttranscriptional level under stress conditions of the beta cells could help to reverse autoimmune process through limiting epitope-presentation to autoreactive T cells. In this article we provide a perspective on the role of ERAP1 as implicated in the pathogenesis of insulin-dependent diabetes mellitus by reviewing studies reported in literature and discussing our own experimental evidence.

The Landscape of Accessible Chromatin and Developmental Transcriptome Maps Reveal a Genetic Mechanism of Skeletal Muscle Development in Pigs

The epigenetic regulation mechanism of porcine skeletal muscle development relies on the openness of chromatin and is also precisely regulated by transcriptional machinery. However, fewer studies have exploited the temporal changes in gene expression and the landscape of accessible chromatin to reveal the underlying molecular mechanisms controlling muscle development. To address this, skeletal muscle biopsy samples were taken from Landrace pigs at days 0 (D0), 60 (D60), 120 (D120), and 180 (D180) after birth and were then analyzed using RNA-seq and ATAC-seq. The RNA-seq analysis identified 8554 effective differential genes, among which ACBD7, TMEM220, and ATP1A2 were identified as key genes related to the development of porcine skeletal muscle. Some potential cis-regulatory elements identified by ATAC-seq analysis contain binding sites for many transcription factors, including SP1 and EGR1, which are also the predicted transcription factors regulating the expression of ACBD7 genes. Moreover, the omics analyses revealed regulatory regions that become ectopically active after birth during porcine skeletal muscle development after birth and identified 151,245, 53,435, 30,494, and 40,911 peaks. The enriched functional elements are related to the cell cycle, muscle development, and lipid metabolism. In summary, comprehensive high-resolution gene expression maps were developed for the transcriptome and accessible chromatin during postnatal skeletal muscle development in pigs.

ISL1 controls pancreatic alpha cell fate and beta cell maturation

BACKGROUND

Glucose homeostasis is dependent on functional pancreatic α and ß cells. The mechanisms underlying the generation and maturation of these endocrine cells remain unclear.

RESULTS

We unravel the molecular mode of action of ISL1 in controlling α cell fate and the formation of functional ß cells in the pancreas. By combining transgenic mouse models, transcriptomic and epigenomic profiling, we uncover that elimination of Isl1 results in a diabetic phenotype with a complete loss of α cells, disrupted pancreatic islet architecture, downregulation of key ß-cell regulators and maturation markers of ß cells, and an enrichment in an intermediate endocrine progenitor transcriptomic profile.

CONCLUSIONS

Mechanistically, apart from the altered transcriptome of pancreatic endocrine cells, Isl1 elimination results in altered silencing H3K27me3 histone modifications in the promoter regions of genes that are essential for endocrine cell differentiation. Our results thus show that ISL1 transcriptionally and epigenetically controls α cell fate competence, and ß cell maturation, suggesting that ISL1 is a critical component for generating functional α and ß cells.