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

Integrated analysis of ATAC-seq and RNA-seq reveals the transcriptional regulation network in SLE

BACKGROUND

CD4⁺ T cells have a vital role in the pathogenesis of systemic lupus erythematosus (SLE), abnormal gene expression in CD4⁺ T cells partly accounting for dysfunctional CD4⁺T cells. However, the underying regulatory mechanisms of abnormal gene expression in CD4⁺ T cells derived from SLE patients are not fully understood.

METHODS

The peripheral blood CD4⁺ T cells were acquired from 4 SLE patients and 4 matched healthy controls. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) was conducted to screen differentially accessible chromatin regions between SLE and normals, and motif prediction was used to identify potentially key transcription factors (TFs) involved in CD4⁺T dysfunction. RNA sequencing (RNA-seq) was performed to screen differentially expressed genes in SLE CD4⁺T cells. ATAC-seq and RNA-seq were integrated to further analyze the relationship between chromatin accessibility and gene expression. KEGG pathway enrichment analysis was to determine enriched pathways of interactions between all predicted TFs and differentially expressed genes (DEGs). Meanwhile, the expression changes of target genes followed by siRNA knockdown of the predicted TF were experimentally verified by qPCR. Finally, the H3K27ac modification levels of immune-related genes with open chromatin and up-regulated expression in SLE CD4⁺T cells was detected by ChIP-qPCR.

RESULTS

We identified 3067 differentially accessible regions (DARs) and 1292 DEGs. TF prediction and functional enrichment analyses showed the TF-gene interaction networks were enriched predominantly in T helper 17 (Th17) cell differentiation, the cell cycle and some signaling pathways. Top 5 TFs were predicted based on overlapping genes between the DAR-related genes and the DEGs: ZNF770, THAP11, ZBTB14, ETV1, POU3F1. Validation experiments indicated that the expression of TRIM25, CD163, BST2, IFIT5, IFITM3, OASL, TBX21, IL15RA and IL12RB2 was significantly downregulated in CD4⁺Tcells with ZNF770 knockdown. H3K27ac showed significantly higher levels in the promoter regions of KLF4 and MX2 in SLE CD4⁺ T cells.

CONCLUSION

These DARs associated with this disease may become targets for future treatment of SLE.

GC-Globulin/Vitamin D-Binding Protein Is Required for Pancreatic α-Cell Adaptation to Metabolic Stress

GC-globulin (GC), or vitamin D-binding protein, is a multifunctional protein involved in the transport of circulating vitamin 25(OH)D and fatty acids, as well as actin scavenging. In the pancreatic islets, the gene encoding GC, GC/Gc, is highly localized to glucagon-secreting α-cells. Despite this, the role of GC in α-cell function is poorly understood. We previously showed that GC is essential for α-cell morphology, electrical activity, and glucagon secretion. We now show that loss of GC exacerbates α-cell failure during metabolic stress. High-fat diet-fed GC-/- mice have basal hyperglucagonemia, which is associated with decreased α-cell size, impaired glucagon secretion and Ca2+ fluxes, and changes in glucose-dependent F-actin remodelling. Impairments in glucagon secretion can be rescued using exogenous GC to replenish α-cell GC levels, increase glucagon granule area, and restore the F-actin cytoskeleton. Lastly, GC levels decrease in α-cells of donors with type 2 diabetes, which is associated with changes in α-cell mass, morphology, and glucagon expression. Together, these data demonstrate an important role for GC in α-cell adaptation to metabolic stress.

Microvascular reconstructions in oral and maxillofacial surgery - Results of a survey among oral and maxillofacial surgeons in Germany, Austria, and Switzerland

This study aimed to evaluate the use of microvascular free flaps (MFF) in oral and maxillofacial surgery (OMFS) in Germany, Austria, and Switzerland. A dynamic online questionnaire, using 42-46 questions, was sent to OMF surgeons based in hospitals in Germany, Austria, and Switzerland. The questionnaire was evaluated internally and externally. Aside from general information, data were collected on organizational aspects, approaches, MFF types and frequency, presurgical planning, intraoperative procedures, perioperative medications, flap monitoring, and patient management. Participants mostly performed 30-40 MFF each year (11/53). Most stated that the COVID-19 pandemic did influence MFF frequency (25/53) to varying extents. Radial forearm flap was most frequently used (37/53), followed by ALT (5/53), and fibula flap (5/53). Primary reconstruction was performed by most participants (35/48). Irradiated bony transplants were mostly used for implant placement after 12 months (23/48). Most participants (38/48) used reconstruction plates, followed by miniplates (36/48), PSI reconstruction (31/48), and PSI miniplates (10/48). Regarding the postoperative use of anticoagulants, low-molecular-weight (37/48) and unfractioned heparins (15/48) were widely used, most often for 3-7 days (26/48). Clinical evaluation was mostly preferred for flap monitoring (47/48), usually every 2 h (34/48), for at least 48 h (19/48). Strong heterogeneity in MFF reconstructions in OMFS was found, especially regarding the timepoints of reconstruction, types of osteosynthesis, and postoperative MFF management. These findings provide the chance to further compare the different treatment algorithms regarding relevant MFF aspects, such as postoperative management. This could create evidence-based treatment algorithms that will further improve the clinical outcomes in MFF reconstructions.

Direct Reprogramming of Somatic Cells into Induced β-Cells: An Overview

The persistent shortage of insulin-producing islet mass or β-cells for transplantation in the ever-growing diabetic population worldwide is a matter of concern. To date, permanent cure to this medical complication is not available and soon after the establishment of lineage-specific reprogramming, direct β-cell reprogramming became a viable alternative for β-cell regeneration. Direct reprogramming is a straightforward and powerful technique that can provide an unlimited supply of cells by transdifferentiating terminally differentiated cells toward the desired cell type. This approach has been extensively used by multiple groups to reprogram non-β-cells toward insulin-producing β-cells. The β-cell identity has been achieved by various studies via ectopic expression of one or more pancreatic-specific transcription factors in somatic cells, bypassing the pluripotent state. This work highlights the importance of the direct reprogramming approaches (both integrative and non-integrative) in generating autologous β-cells for various applications. An in-depth understanding of the strategies and cell sources could prove beneficial for the efficient generation of integration-free functional insulin-producing β-cells for diabetic patients lacking endogenous β-cells.

RPS6 transcriptional modulation in neural tissues of Nauphoeta cinerea during streptozotocin-associated sugar metabolism impairment

The use of insects to model molecular events that characterize degenerative conditions was originally met with scepticism. However, the discovery of insect insulin-like peptides in the 1970's and the demonstration of evolutionary conservation of insulin-related signalling from insects to mammals have highlighted the importance and reduced cost of insect models in biomedical research. Here, we expand on our earlier described modelling of streptozotocin-induced brain glucose metabolic disruption in Nauphoeta cinerea, using RNA-sequencing analysis to study the transcriptional and genetic signatures of degeneration and stress signalling when glucose levels are elevated in the brain of the lobster cockroach. Nymphs were randomly divided into three groups: Control (0.8% NaCl), and two single streptozotocin injection doses (74 nmol and 740 nmol). The transcriptional analyses featured a dysregulation of 226 genes at high dose STZ treatment and 278 genes at the low dose. Our mRNA-sequencing data showed that ribosomal protein genes were the most upregulated genes at both 74 and 740 nmol STZ treatment. We therefore used RT-qPCR and relative transcriptional methods to validate our proposed mechanism of brain glucose toxicity-induced degeneration in Nauphoeta cinerea, which involved the upregulation of ribosomal proteins and rpS6 regulators (mTORC1, protein kinases, casein kinase 1 and Death-associated protein kinase), the upregulation of MAPK cascades (RAS, ERK, P38 and JNK), alongside the downregulation of the PI3K/AKT cascade. Taken together, this study highlights the remarkable opportunity for Nauphoeta cinerea use as an experimental organism in hyperglycaemia, degeneration, and stress signalling.

ATAC-Seq Analysis of Accessible Chromatin: From Experimental Steps to Data Analysis

Accessible chromatin often represents gene regulatory elements, including promoters and enhancers, essential for gene expression. Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is one of the most popular techniques to investigate chromatin accessibility across the genome. Here we describe, step by step, a series of optimized experimental methods and bioinformatics pipelines for ATAC-seq analysis. As an example, we present an analysis of murine spermatogenic cells: a method to isolate germ cells, a reaction step using Tn5 transposase to insert sequencing adapters into accessible DNA, a library preparation method for high-throughput sequencing, and bioinformatics analysis of sequencing data. Overall, we introduce a framework of ATAC-seq analysis that can be applied to any cell population to identify cell-type-specific gene regulatory elements and their cis-regulatory networks.

Understanding islet dysfunction in type 2 diabetes through multidimensional pancreatic phenotyping: The Human Pancreas Analysis Program

In this perspective, we provide an overview of a recently established National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) initiative, the Human Pancreas Analysis Program for Type 2 Diabetes (HPAP-T2D). This program is designed to define the molecular pathogenesis of islet dysfunction by studying human pancreatic tissue samples from organ donors with T2D. HPAP-T2D generates detailed datasets of physiological, histological, transcriptomic, epigenomic, and genomic information. Importantly, all data collected, generated, and analyzed by HPAP-T2D are made immediately and freely available through a centralized database, PANC-DB, thus providing a comprehensive data resource for the diabetes research community.

T3 and glucose increase expression of phosphoenolpyruvate carboxykinase (PCK1) leading to increased β-cell proliferation

OBJECTIVES

Thyroid hormone (T3) and high glucose concentrations are critical components of β-cell maturation and function. In the present study, we asked whether T3 and glucose signaling pathways coordinately regulate transcription of genes important for β-cell function and proliferation.

METHODS

RNA-seq analysis was performed on cadaveric human islets from five different donors in response to low and high glucose concentrations and in the presence or absence of T3. Gene expression was also studies in sorted human β-cells, mouse islets and Ins-1 cells by RT-qPCR. Silencing of the thyroid hormone receptors (THR) was conducted using lentiviruses. Proliferation was assessed by ki67 immunostaining in primary human/mouse islets. Chromatin immunoprecipitation and proximity ligation assay were preformed to validate interactions of ChREBP and THR.

RESULTS

We found glucose-mediated expression of carbohydrate response element binding protein alpha and beta (ChREBPα and ChREBPβ) mRNAs and their target genes are highly dependent on T3 concentrations in rodent and human β-cells. In β-cells, T3 and glucose coordinately regulate the expression of ChREBPβ and PCK1 (phosphoenolpyruvate carboxykinase-1) among other important genes for β-cell maturation. Additionally, we show the thyroid hormone receptor (THR) and ChREBP interact, and their relative response elements are located near to each other on mutually responsive genes. In FACS-sorted adult human β-cells, we found that high concentrations of glucose and T3 induced the expression of PCK1. Next, we show that overexpression of Pck1 together with dimethyl malate (DMM), a substrate precursor, significantly increased β-cell proliferation in human islets. Finally, using a Cre-Lox approach, we demonstrated that ChREBPβ contributes to Pck1-dependent β-cell proliferation in mouse β-cells.

CONCLUSIONS

We conclude that T3 and glucose act together to regulate ChREBPβ, leading to increased expression and activity of Pck1, and ultimately increased β-cell proliferation.

The early function of cortisol in liver during Aeromonas hydrophila infection: Dynamics of the transcriptome and accessible chromatin landscapes

In China, channel catfish (Ictalurus punctatus) is an important aquaculture species; however, haemorrhagic disease (Aeromonas hydrophila induced disease) in these fish has caused tremendous economic loss due to high morbidity and mass mortality in the breeding industry. The role of cortisol in bacterial diseases, particularly in the acute phase, remains unclear. In this study, liver transcriptome (RNA-seq) and chromatin accessibility (ATAC-seq) analyses were employed to investigate the early functional role of cortisol in Aeromonas hydrophila-stimulated responses. Our experiments confirmed that A. hydrophila infection can initially significantly increase serum cortisol levels at 1 h after infection. At this time point, the increased serum cortisol levels can significantly regulate A. hydrophila-regulated genes by affecting both transcriptome and chromatin accessibility. Cross-analysis of RNA-seq and ATAC-seq revealed that a certain gene group (92 target_DEGs) was regulated at an early time point by cortisol. KEGG enrichment analysis revealed that the top three pathways according to target_DEGs were cancer, glutathione metabolism, and the Notch signalling pathway. The protein-protein interaction analysis of target_DEGs revealed that they may be primarily involved in cell proliferation, CD8⁺ T cell function, glutathione synthesis, and activation of the NF-κB signalling pathway. This suggests that after the emergence of immune stress, the early regulation of cortisol is positive against the immune response. It is possible that in this situation, the animal is attempting to avoid dangerous situations and risks and then cope with the imbalance produced by the stressor to ultimately restore homeostasis. Our results will contribute to future research on fish and provide valuable insight regarding the mechanism of immune regulation by cortisol and the study of bacterial haemorrhagic disease in channel catfish.

Altered microvasculature in pancreatic islets from subjects with type 1 diabetes

AIMS

The transcriptome of different dissociated pancreatic islet cells has been described in enzymatically isolated islets in both health and disease. However, the isolation, culturing, and dissociation procedures likely affect the transcriptome profiles, distorting the biological conclusions. The aim of the current study was to characterize the cells of the islets of Langerhans from subjects with and without type 1 diabetes in a way that reflects the in vivo situation to the highest possible extent.

METHODS

Islets were excised using laser capture microdissection directly from frozen pancreatic tissue sections obtained from organ donors with (n = 7) and without (n = 8) type 1 diabetes. Transcriptome analysis of excised samples was performed using AmpliSeq. Consecutive pancreatic sections were used to estimate the proportion of beta-, alpha-, and delta cells using immunofluorescence and to examine the presence of CD31 positive endothelial regions using immunohistochemistry.

RESULTS

The proportion of beta cells in islets from subjects with type 1 diabetes was reduced to 0% according to both the histological and transcriptome data, and several alterations in the transcriptome were derived from the loss of beta cells. In total, 473 differentially expressed genes were found in the islets from subjects with type 1 diabetes. Functional enrichment analysis showed that several of the most upregulated gene sets were related to vasculature and angiogenesis, and histologically, vascular density was increased in subjects with type 1 diabetes. Downregulated in type 1 diabetes islets was the gene set epithelial mesenchymal transition.

CONCLUSION

A number of transcriptional alterations are present in islets from subjects with type 1 diabetes. In particular, several gene sets related to vasculature and angiogenesis are upregulated and there is an increased vascular density, suggesting an altered microvasculature in islets from subjects with type 1 diabetes. By studying pancreatic islets extracted directly from snap-frozen pancreatic tissue, this study reflects the in vivo situation to a high degree and gives important insights into islet pathophysiology in type 1 diabetes.

Circular RNAs in diabetes and its complications: Current knowledge and future prospects

A novel class of non-coding RNA transcripts called circular RNAs (circRNAs) have been the subject of significant recent studies. Accumulating evidence points that circRNAs play an important role in the cellular processes, inflammatory expression, and immune responses through sponging miRNA, binding, or translating in proteins. Studies have found that circRNAs are involved in the physiologic and pathologic processes of diabetes. There has been an increased focus on the relevance of between abnormal circRNA expression and the development and progression of various types of diabetes and diabetes-related diseases. These circRNAs not only serve as promising diagnostic and prognostic molecular biomarkers, but also have important biological roles in islet cells, diabetes, and its complications. In addition, many circRNA signaling pathways have been found to regulate the occurrence and development of diabetes. Here we comprehensively review and discuss recent advances in our understanding of the physiologic function and regulatory mechanisms of circRNAs on pancreatic islet cells, different subtypes in diabetes, and diabetic complications.

Sex Differences in the Molecular Programs of Pancreatic Cells Contribute to the Differential Risks of Type 2 Diabetes

Epidemiology studies demonstrate that women are at a significantly lower risk of developing type 2 diabetes (T2D) compared to men. However, the molecular basis of this risk difference is not well understood. In this study, we examined the sex differences in the genetic programs of pancreatic endocrine cells. We combined pancreas perifusion data and single-cell genomic data from our laboratory and from publicly available data sets to investigate multiple axes of the sex differences in the human pancreas at the single-cell type and single-cell level. We systematically compared female and male islet secretion function, gene expression program, and regulatory principles of pancreatic endocrine cells. The perifusion data indicate that female endocrine cells have a higher secretion capacity than male endocrine cells. Single-cell RNA-sequencing analysis suggests that endocrine cells in male controls have molecular signatures that resemble T2D. In addition, we identified genomic elements associated with genome-wide association study T2D loci to have differential accessibility between female and male delta cells. These genomic elements may play a sex-specific causal role in the pathogenesis of T2D. We provide molecular mechanisms that explain the differential risk of T2D between women and men. Knowledge gained from our study will accelerate the development of diagnostics and therapeutics in sex-aware precision medicine for diabetes.

Dynamics of transcriptome and chromatin accessibility revealed sequential regulation of potential transcription factors during the brown adipose tissue whitening in rabbits

Brown adipose tissue (BAT) represents a valuable target for treating obesity in humans. BAT losses of thermogenic capacity and gains a “white adipose tissue-like (WAT-like)” phenotype (BAT whitening) under thermoneutral environments, which could lead to potential low therapy responsiveness in BAT-based obesity treatments. However, the epigenetic mechanisms of BAT whitening remain largely unknown. In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). RNA-sequencing (RNA-seq) and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were performed to investigate transcriptome and chromatin accessibility of BATs at the four whitening stages, respectively. Our data showed that many genes and chromatin accessible regions (refer to as “peaks”) were identified as significantly changed during BAT whitening in rabbits. The BAT-selective genes downregulated while WAT-selective genes upregulated from D0 to Y2, and the de novo lipogenesis-related genes reached the highest expression levels at D85. Both the highly expressed genes and accessible regions in Y2 were significantly enriched in immune response-related signal pathways. Analysis of different relationships between peaks and their nearby genes found an increased proportion of the synchronous changes between chromatin accessibility and gene expression during BAT whitening. The synergistic changes between the chromatin accessibility of promoter and the gene expression were found in the key adipose genes. The upregulated genes which contained increased peaks were significantly enriched in the PI3K-Akt signaling pathway, steroid biosynthesis, TGF-beta signaling pathway, osteoclast differentiation, and dilated cardiomyopathy. Moreover, the footprinting analysis suggested that sequential regulation of potential transcription factors (TFs) mediated the loss of thermogenic phenotype and the gain of a WAT-like phenotype of BAT. In conclusion, our study provided the transcriptional and epigenetic frameworks for understanding BAT whitening in rabbits for the first time and might facilitate potential insights into BAT-based obesity treatments.

3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.

Decoding mechanism of action and sensitivity to drug candidates from integrated transcriptome and chromatin state

Omics-based technologies are driving major advances in precision medicine, but efforts are still required to consolidate their use in drug discovery. In this work, we exemplify the use of multi-omics to support the development of 3-chloropiperidines, a new class of candidate anticancer agents. Combined analyses of transcriptome and chromatin accessibility elucidated the mechanisms underlying sensitivity to test agents. Furthermore, we implemented a new versatile strategy for the integration of RNA- and ATAC-seq (Assay for Transposase-Accessible Chromatin) data, able to accelerate and extend the standalone analyses of distinct omic layers. This platform guided the construction of a perturbation-informed basal signature predicting cancer cell lines' sensitivity and to further direct compound development against specific tumor types. Overall, this approach offers a scalable pipeline to support the early phases of drug discovery, understanding of mechanisms, and potentially inform the positioning of therapeutics in the clinic.

Integrated Analysis of Transcriptome, microRNAs, and Chromatin Accessibility Revealed Potential Early B-Cell Factor1-Regulated Transcriptional Networks during the Early Development of Fetal Brown Adipose Tissues in Rabbits

In domestic mammals, cold stress decreases the survival rate of newborns and increases the cost of management. Brown adipose tissue (BAT) is the main thermogenic organ for cubs, and well-developed fetal BAT (FBAT) is beneficial for newborns to maintain core temperatures during the first several days of life. However, our knowledge of the epigenetic mechanisms during the early development of FBAT remains largely unknown. Rabbits (Oryctolagus cuniculus) are economically important domestic animals. In this study, a histological analysis shows that the tissue content, thermogenic capacity, and lipid content of FBAT dramatically increases from gestational day 21 (G21) to gestational day 24 (G24) in rabbits. RNA-seq, microRNA-seq (miRNA-seq), and the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) show that many genes, miRNAs, and chromatin-accessible regions (referred to as peaks) were identified as significantly changed from G21 to G24, respectively. The upregulated genes from G21 to G24 were significantly enriched in the mitochondrial metabolism and thermogenesis-related signal pathways. The integrated analysis of transcriptome and chromatin accessibility reveals that the peaks in the promoters have a more regulatory effect than peaks in other genomic elements on the expression of their nearby genes in FBATs. The upregulated genes that are associated with increased chromatin accessibility in the promoter regions are involved in the energy metabolism-related signaling pathways. The genes that have a greater tendency to be regulated by miRNAs than the chromatin accessibility in gene promoters are involved in the apelin, insulin, and endocytosis signaling pathways. Furthermore, genome-wide transcription factor (TF) footprinting analysis identifies early B-cell factor1 (EBF1) as playing a key role during early FBAT development. The carbon metabolism, citrate cycle, and PPAR signaling pathways are significantly enriched by the predicted EBF1-regulated cascade TF-network. In conclusion, our work provides a framework for understanding epigenetics regulatory mechanisms underlying the early development of FBAT and identifies potential TF involved in the early development of FBAT in rabbits.

Variant-to-gene-mapping analyses reveal a role for pancreatic islet cells in conferring genetic susceptibility to sleep-related traits

We investigated the potential role of sleep-trait associated genetic loci in conferring a degree of their effect via pancreatic α- and β-cells, given that both sleep disturbances and metabolic disorders, including type 2 diabetes and obesity, involve polygenic contributions and complex interactions. We determined genetic commonalities between sleep and metabolic disorders, conducting linkage disequilibrium genetic correlation analyses with publicly available GWAS summary statistics. Then we investigated possible enrichment of sleep-trait associated SNPs in promoter-interacting open chromatin regions within α- and β-cells, intersecting public GWAS reports with our own ATAC-seq and high-resolution promoter-focused Capture C data generated from both sorted human α-cells and an established human beta-cell line (EndoC-βH1). Finally, we identified putative effector genes physically interacting with sleep-trait associated variants in α- and EndoC-βH1cells running variant-to-gene mapping and establish pathways in which these genes are significantly involved. We observed that insomnia, short and long sleep-but not morningness-were significantly correlated with type 2 diabetes, obesity and other metabolic traits. Both the EndoC-βH1 and α-cells were enriched for insomnia loci (p = .01; p = .0076), short sleep loci (p = .017; p = .022) and morningness loci (p = 2.2 × 10-7; p = .0016), while the α-cells were also enriched for long sleep loci (p = .034). Utilizing our promoter contact data, we identified 63 putative effector genes in EndoC-βH1 and 76 putative effector genes in α-cells, with these genes showing significant enrichment for organonitrogen and organophosphate biosynthesis, phosphatidylinositol and phosphorylation, intracellular transport and signaling, stress responses and cell differentiation. Our data suggest that a subset of sleep-related loci confer their effects via cells in pancreatic islets.

Recent Progress in the Diagnosis and Management of Type 2 Diabetes Mellitus in the Era of COVID-19 and Single Cell Multi-Omics Technologies

Type 2 diabetes mellitus (T2DM) is one of the world's leading causes of death and life-threatening conditions. Therefore, we review the complex vicious circle of causes responsible for T2DM and risk factors such as the western diet, obesity, genetic predisposition, environmental factors, and SARS-CoV-2 infection. The prevalence and economic burden of T2DM on societal and healthcare systems are dissected. Recent progress on the diagnosis and clinical management of T2DM, including both non-pharmacological and latest pharmacological treatment regimens, are summarized. The treatment of T2DM is becoming more complex as new medications are approved. This review is focused on the non-insulin treatments of T2DM to reach optimal therapy beyond glycemic management. We review experimental and clinical findings of SARS-CoV-2 risks that are attributable to T2DM patients. Finally, we shed light on the recent single-cell-based technologies and multi-omics approaches that have reached breakthroughs in the understanding of the pathomechanism of T2DM.

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

Genome-Wide Analysis of Yeast Metabolic Cycle through Metabolic Network Models Reveals Superiority of Integrated ATAC-seq Data over RNA-seq Data

Saccharomyces cerevisiae undergoes robust oscillations to regulate its physiology for adaptation and survival under nutrient-limited conditions. Environmental cues can induce rhythmic metabolic alterations in order to facilitate the coordination of dynamic metabolic behaviors. Of such metabolic processes, the yeast metabolic cycle enables adaptation of the cells to varying nutritional status through oscillations in gene expression and metabolite production levels. In this process, yeast metabolism is altered between diverse cellular states based on changing oxygen consumption levels: quiescent (reductive charging [RC]), growth (oxidative [OX]), and proliferation (reductive building [RB]) phases. We characterized metabolic alterations during the yeast metabolic cycle using a variety of approaches. Gene expression levels are widely used for condition-specific metabolic simulations, whereas the use of epigenetic information in metabolic modeling is still limited despite the clear relationship between epigenetics and metabolism. This prompted us to investigate the contribution of epigenomic information to metabolic predictions for progression of the yeast metabolic cycle. In this regard, we determined altered pathways through the prediction of regulated reactions and corresponding model genes relying on differential chromatin accessibility levels. The predicted metabolic alterations were confirmed via data analysis and literature. We subsequently utilized RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data sets in the contextualization of the yeast model. The use of ATAC-seq data considerably enhanced the predictive capability of the model. To the best of our knowledge, this is the first attempt to use genome-wide chromatin accessibility data in metabolic modeling. The preliminary results showed that epigenomic data sets can pave the way for more accurate metabolic simulations.

IMPORTANCE Dynamic chromatin organization mediates the emergence of condition-specific phenotypes in eukaryotic organisms. Saccharomyces cerevisiae can alter its metabolic profile via regulation of genome accessibility and robust transcriptional oscillations under nutrient-limited conditions. Thus, both epigenetic information and transcriptomic information are crucial in the understanding of condition-specific metabolic behavior in this organism. Based on genome-wide alterations in chromatin accessibility and transcription, we investigated the yeast metabolic cycle, which is a remarkable example of coordinated and dynamic yeast behavior. In this regard, we assessed the use of ATAC-seq and RNA-seq data sets in condition-specific metabolic modeling. To our knowledge, this is the first attempt to use chromatin accessibility data in the reconstruction of context-specific metabolic models, despite the extensive use of transcriptomic data. As a result of comparative analyses, we propose that the incorporation of epigenetic information is a promising approach in the accurate prediction of metabolic dynamics.

The Nuclear DNA Sensor IFI16 Indiscriminately Binds to and Diminishes Accessibility of the HSV-1 Genome to Suppress Infection

Human cells identify invading pathogens and activate immune signaling pathways through a wide array of pattern recognition receptors, including DNA sensors. The interferon-inducible protein 16 (IFI16) is a nuclear DNA sensor that recognizes double-stranded DNA from a number of viral sources, including genomes of nuclear-replicating viruses. Among these is the prevalent human pathogen herpes simplex virus 1 (HSV-1). Upon binding to the HSV-1 DNA genome, IFI16 both induces antiviral cytokine expression and suppresses virus gene expression. Here, we used a multiomics approach of DNA sequencing techniques paired with targeted mass spectrometry to obtain an extensive view of the interaction between IFI16 and the HSV-1 genome and how this binding affects the viral DNA structure and protein expression. Through chromatin immunoaffinity purification coupled with next-generation DNA sequencing (ChIP-seq), we found that IFI16 binds to the HSV-1 genome in a sequence-independent manner while simultaneously exhibiting broad enrichment at two loci: UL30, the viral DNA polymerase gene, and US1 to US7. The assay for transposase-accessible chromatin with sequencing (ATAC-seq) revealed that these two regions are among the most accessible stretches of DNA on the genome, thereby facilitating IFI16 binding. Accessibility of the entire HSV-1 genome is elevated upon IFI16 knockout, indicating that expression of IFI16 globally induces chromatinization of viral DNA. Deletion of IFI16 also results in a global increase in the expression of HSV-1 proteins, as measured by parallel reaction monitoring-mass spectrometry of viral proteins representing 80% of the HSV-1 genome. Altogether, we demonstrate that IFI16 interacts with the HSV-1 genome in a sequence-independent manner, coordinating epigenetic silencing of the viral genome and decreasing protein expression and virus replication. IMPORTANCE Mammalian host defense against viral infection includes broad-acting cellular restriction factors, as well as effectors of intrinsic and innate immunity. IFI16 is a critical nuclear host defense factor and intrinsic immune protein involved in binding viral DNA genomes, thereby repressing the replication of nucleus-replicating viruses, including the human herpes simplex virus 1. What has remained unclear is where on the viral genome IFI16 binds and how binding affects both viral DNA structural accessibility and viral protein expression. Our study provides a global view of where and how a nuclear restriction factor of DNA viruses associates with viral genomes to exert antiviral functions during early stages of an acute virus infection. Our study can additionally serve as a systems-level model to evaluate nuclear DNA sensor interactions with viral genomes, as well as the antiviral outcomes of transcriptionally silencing pathogen-derived DNA.

The influence of phytochemicals on cell heterogeneity in chronic inflammation-associated diseases: the prospects of single cell sequencing

Chronic inflammation-associated diseases include, but is not limited to cardiovascular disease, cancer, obesity, diabetes, etc. Cell heterogeneity is a prerequisite for understanding the physiological and pathological development of cell metabolism, and its response to external stimuli. Recently, dietary habits based on phytochemicals became increasingly recognized to play a pivotal role in chronic inflammation. Phytochemicals can relieve chronic inflammation by regulating inflammatory cell differentiation and immune cell response, but the influence of phytochemicals on cell heterogeneity from in vitro and ex vivo studies cannot simulate the complexity of cell differentiation in vivo due to the differences in cell lines and extracellular environment. Therefore, there is no consensus on the regulation mechanism of phytochemicals on chronic diseases based on cell heterogeneity. The purpose of this review is to summarize cell heterogeneity in common chronic inflammation-associated diseases and trace the effects of phytochemicals on cell differentiation in chronic diseases development. More importantly, by discussing the problems and challenges which hinder the study of cell heterogeneity in recent nutritional assessment experiments, we propose new prospects based on the drawbacks of existing research to optimize the research on the regulation mechanism of phytochemicals on chronic diseases. The need to explore precise measurements of cell heterogeneity is a key pillar in understanding the influence of phytochemicals on certain diseases. In the future, deeper understanding of cell-to-cell variation and the impact of food components and their metabolites on cell function by single-cell genomics and epigenomics with the focus on individual differences will open new avenues for the next generation of health care.

Overcoming the Limitations of Stem Cell-Derived Beta Cells

Great advances in type 1 diabetes (T1D) and type 2 diabetes (T2D) treatment have been made to this day. However, modern diabetes therapy based on insulin injections and cadaveric islets transplantation has many disadvantages. That is why researchers are developing new methods to regenerate the pancreatic hormone-producing cells in vitro. The most promising approach is the generation of stem cell-derived beta cells that could provide an unlimited source of insulin-secreting cells. Recent studies provide methods to produce beta-like cell clusters that display glucose-stimulated insulin secretion—one of the key characteristics of the beta cell. However, in comparison with native beta cells, stem cell-derived beta cells do not undergo full functional maturation. In this paper we review the development and current state of various protocols, consider advantages, and propose ways to improve them. We examine molecular pathways, epigenetic modifications, intracellular components, and the microenvironment as a possible leverage to promote beta cell functional maturation. A possibility to create islet organoids from stem cell-derived components, as well as their encapsulation and further transplantation, is also examined. We try to combine modern research on beta cells and their crosstalk to create a holistic overview of developing insulin-secreting systems.

Harnessing changes in open chromatin determined by ATAC-seq to generate insulin-responsive reporter constructs

BACKGROUND

Gene regulation is critical for proper cellular function. Next-generation sequencing technology has revealed the presence of regulatory networks that regulate gene expression and essential cellular functions. Studies investigating the epigenome have begun to uncover the complex mechanisms regulating transcription. Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is quickly becoming the assay of choice for many epigenomic investigations. However, whether intervention-mediated changes in accessible chromatin determined by ATAC-seq can be harnessed to generate intervention-inducible reporter constructs has not been systematically assayed.

RESULTS

We used the insulin signaling pathway as a model to investigate chromatin regions and gene expression changes using ATAC- and RNA-seq in insulin-treated Drosophila S2 cells. We found correlations between ATAC- and RNA-seq data, especially when stratifying differentially-accessible chromatin regions by annotated feature type. In particular, our data demonstrated a weak but significant correlation between chromatin regions annotated to enhancers (1-2 kb from the transcription start site) and downstream gene expression. We cloned candidate enhancer regions upstream of luciferase and demonstrate insulin-inducibility of several of these reporters.

CONCLUSIONS

Insulin-induced chromatin accessibility determined by ATAC-seq reveals enhancer regions that drive insulin-inducible reporter gene expression.

Chromatin accessibility and transcriptome integrative analysis revealed AP-1-mediated genes potentially modulate histopathology features in psoriasis

Background

Psoriasis is a chronic and hyperproliferative skin disease featured by hyperkeratosis with parakeratosis, Munro micro-abscess, elongation of rete pegs, granulosa thinning, and lymphocyte infiltration. We previously profiled gene expression and chromatin accessibility of psoriatic skins by transcriptome sequencing and ATAC-seq. However, integrating both of these datasets to unravel gene expression regulation is lacking. Here, we integrated transcriptome and ATAC-seq of the same psoriatic and normal skin tissues, trying to leverage the potential role of chromatin accessibility and their function in histopathology features.

Results

By inducing binding and expression target analysis (BETA) algorithms, we explored the target prediction of transcription factors binding in 15 psoriatic and 19 control skins. BETA identified 408 upregulated genes (rank product < 0.01) and 133 downregulated genes linked with chromatin accessibility. We noticed that cumulative fraction of genes in upregulation group was statistically higher than background, while that of genes in downregulation group was not significant. KEGG pathway analysis showed that the upregulated 408 genes were enriched in TNF, NOD, and IL-17 signaling pathways. In addition, the motif module in BETA suggested the 57 upregulated genes are targeted by transcription factor AP-1, indicating that increased chromatin accessibility facilitated the binding of AP-1 to the target regions and further induced expression of relevant genes. Among these genes, SQLE, STRN, EIF4, and MYO1B expression was increased in patients with hyperkeratosis, parakeratosis, and acanthosis thickening.

Conclusions

In summary, with the advantage of BETA, we identified a series of genes that contribute to the disease pathogenesis, especially in modulating histopathology features, providing us with new clues in treating psoriasis.

Analysis of chromatin accessibility in p53 deficient spermatogonial stem cells for high frequency transformation into pluripotent state

OBJECTIVES

Spermatogonial stem cells (SSCs), the germline stem cells (GSCs) committed to spermatogenesis in niche, can transform into pluripotent state in long-term culture without introduction of exogenous factors, typically in p53 deficiency condition. As the guardian for genomic stability, p53 is associated with epigenetic alterations during SSCs transformation. However, the mechanism is still unknown, since complicated roles of p53 baffle our understanding of the regulating process.

MATERIALS AND METHODS

The chromatin accessibility and differentially expressed genes (DEGs) were analysed in p53+/+ and p53-/- SSCs using the Assay for Transposase-Accessible Chromatin with high-throughput Sequencing (ATAC-seq) and RNA-sequencing (RNA-seq), to explore the connection of p53 and cell fate at chromosomal level.

RESULTS

Several transcription factors (TFs), such as CTCF, SMAD3 and SOX2, were predicted as important factors mediating the transformation. Molecular evidence suggested that SMAD3 efficiently promoted pluripotency-associated gene expression both in fresh and long-term cultured SSCs. However, p53 knockout (KO) is insufficient to induce SMAD3 expression in SSCs.

CONCLUSIONS

These observations indicate that SMAD3 is a key factor for SSCs transformation, and an unknown event is required to activate SMAD3 as the prerequisite for SSCs reprogramming, which may occur in the long-term culture of SSCs. This study demonstrates the connection of p53 and pluripotency-associated factors, providing new insight for understanding the mechanisms of SSCs reprogramming and germline tumorigenesis.

Alterations in chromatin accessibility during osteoblast and adipocyte differentiation in human mesenchymal stem cells

Although differential expression of genes is apparent during the adipogenic/osteogenic differentiation of marrow mesenchymal stem cells (MSCs), it is not known whether this is associated with changes in chromosomal structure. In this study, we used ATAC-sequencing technology to observe variations in chromatin assembly during the early stages of MSC differentiation. This showed significant changes in the number and distribution of chromosome accessibility at different time points of adipogenic/osteogenic differentiation. Sequencing of differential peaks indicated alterations in transcription factor motifs involved in MSC differentiation. Gene Ontology (GO) and pathway analysis indicated that changes in biological function resulted from the alterations in chromatin accessibility. We then integrated ATAC-seq and RNA-seq and found that only a small proportion of the overlapping genes were screened out from ATAC-seq and RNA-seq overlapping. Through GO and pathway analysis of these overlapped genes, we not only observed some known biological functions related to adipogenic/osteogenic differentiation but also noticed some unusual biological clustering during MSC differentiation. In summary, our work not only presents the landscape of chromatin accessibility of MSC during differentiation but also helps to further our understanding of the underlying mechanisms of gene expression in these processes.

Mitogen Synergy: An Emerging Route to Boosting Human Beta Cell Proliferation

Decreased number and function of beta cells are a key aspect of diabetes mellitus (diabetes), a disease that remains an onerous global health problem. Means of restoring beta cell mass are urgently being sought as a potential cure for diabetes. Several strategies, such as de novo beta cell derivation via pluripotent stem cell differentiation or mature somatic cell transdifferentiation, have yielded promising results. Beta cell expansion is another promising strategy, rendered challenging by the very low proliferative capacity of beta cells. Many effective mitogens have been identified in rodents, but the vast majority do not have similar mitogenic effects in human beta cells. Extensive research has led to the identification of several human beta cell mitogens, but their efficacy and specificity remain insufficient. An approach based on the simultaneous application of several mitogens has recently emerged and can yield human beta cell proliferation rates of up to 8%. Here, we discuss recent advances in restoration of the beta cell population, focusing on mitogen synergy, and the contribution of RNA-sequencing (RNA-seq) to accelerating the elucidation of signaling pathways in proliferating beta cells and the discovery of novel mitogens. Together, these approaches have taken beta cell research up a level, bringing us closer to a cure for diabetes.

Transcriptional profile and chromatin accessibility in zebrafish melanocytes and melanoma tumors

Transcriptional and epigenetic characterization of melanocytes and melanoma cells isolated from their in vivo context promises to unveil key differences between these developmentally related normal and cancer cell populations. We therefore engineered an enhanced Danio rerio (zebrafish) melanoma model with fluorescently labeled melanocytes to allow for isolation of normal (wild type) and premalignant (BRAFV600E-mutant) populations for comparison to fully transformed BRAFV600E-mutant, p53 loss-of-function melanoma cells. Using fluorescence-activated cell sorting to isolate these populations, we performed high-quality RNA- and ATAC-seq on sorted zebrafish melanocytes vs. melanoma cells, which we provide as a resource here. Melanocytes had consistent transcriptional and accessibility profiles, as did melanoma cells. Comparing melanocytes and melanoma, we note 4128 differentially expressed genes and 56,936 differentially accessible regions with overall gene expression profiles analogous to human melanocytes and the pigmentation melanoma subtype. Combining the RNA- and ATAC-seq data surprisingly revealed that increased chromatin accessibility did not always correspond with increased gene expression, suggesting that though there is widespread dysregulation in chromatin accessibility in melanoma, there is a potentially more refined gene expression program driving cancerous melanoma. These data serve as a resource to identify candidate regulators of the normal vs. diseased states in a genetically controlled in vivo context.

Stem Cell-Derived β Cells: A Versatile Research Platform to Interrogate the Genetic Basis of β Cell Dysfunction

Pancreatic β cell dysfunction is a central component of diabetes progression. During the last decades, the genetic basis of several monogenic forms of diabetes has been recognized. Genome-wide association studies (GWAS) have also facilitated the identification of common genetic variants associated with an increased risk of diabetes. These studies highlight the importance of impaired β cell function in all forms of diabetes. However, how most of these risk variants confer disease risk, remains unanswered. Understanding the specific contribution of genetic variants and the precise role of their molecular effectors is the next step toward developing treatments that target β cell dysfunction in the era of personalized medicine. Protocols that allow derivation of β cells from pluripotent stem cells, represent a powerful research tool that allows modeling of human development and versatile experimental designs that can be used to shed some light on diabetes pathophysiology. This article reviews different models to study the genetic basis of β cell dysfunction, focusing on the recent advances made possible by stem cell applications in the field of diabetes research.

Identification of direct transcriptional targets of NFATC2 that promote β cell proliferation

The transcription factor NFATC2 induces β cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced β cell proliferation, suggesting the FOXP family works to regulate β cell proliferation in concert with NFATC2. NFATC2 induced β cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce β cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate β cell proliferation.

SWI/SNF (BAF) complexes: From framework to a functional role in endothelial mechanotransduction

Endothelial cells (ECs) are constantly subjected to an array of mechanical cues, especially shear stress, due to their luminal placement in the blood vessels. Blood flow can regulate various aspects of endothelial biology and pathophysiology by regulating the endothelial processes at the transcriptomic, proteomic, miRNomic, metabolomics, and epigenomic levels. ECs sense, respond, and adapt to altered blood flow patterns and shear profiles by specialized mechanisms of mechanosensing and mechanotransduction, resulting in qualitative and quantitative differences in their gene expression. Chromatin-regulatory proteins can regulate transcriptional activation by modifying the organization of nucleosomes at promoters, enhancers, silencers, insulators, and locus control regions. Recent research efforts have illustrated that SWI/SNF (SWItch/Sucrose Non-Fermentable) or BRG1/BRM-associated factor (BAF) complex regulates DNA accessibility and chromatin structure. Since the discovery, the gene-regulatory mechanisms of the BAF complex associated with chromatin remodeling have been intensively studied to investigate its role in diverse disease phenotypes. Thus far, it is evident that (1) the SWI/SNF complex broadly regulates the activity of transcriptional enhancers to control lineage-specific differentiation and (2) mutations in the BAF complex proteins lead to developmental disorders and cancers. It is unclear if blood flow can modulate the activity of SWI/SNF complex to regulate EC differentiation and reprogramming. This review emphasizes the integrative role of SWI/SNF complex from a structural and functional standpoint with a special reference to cardiovascular diseases (CVDs). The review also highlights how regulation of this complex by blood flow can lead to the discovery of new therapeutic interventions for the treatment of endothelial dysfunction in vascular diseases.

The Human Islet: Mini-Organ With Mega-Impact

This review focuses on the human pancreatic islet-including its structure, cell composition, development, function, and dysfunction. After providing a historical timeline of key discoveries about human islets over the past century, we describe new research approaches and technologies that are being used to study human islets and how these are providing insight into human islet physiology and pathophysiology. We also describe changes or adaptations in human islets in response to physiologic challenges such as pregnancy, aging, and insulin resistance and discuss islet changes in human diabetes of many forms. We outline current and future interventions being developed to protect, restore, or replace human islets. The review also highlights unresolved questions about human islets and proposes areas where additional research on human islets is needed.

A Functional Polymorphism in Accessible Chromatin Region Confers Risk of Non-Small Cell Lung Cancer in Chinese Population

Background

The disease-associated non-coding variants identified by genome-wide association studies (GWASs) were enriched in open chromatin regions (OCRs) and implicated in gene regulation. Genetic variants in OCRs thus may exert regulatory functions and contribute to non-small cell lung cancer (NSCLC) susceptibility.

Objective

To fine map potential functional variants in GWAS loci that contribute to NSCLC predisposition using chromatin accessibility and histone modification data and explore their functions by population study and biochemical experimental analyses.

Methods

We mapped the chromatin accessible regions of lung tissues using data of assay for transposase-accessible chromatin using sequencing (ATAC-seq) in The Cancer Genome Atlas (TCGA) and prioritized potential regulatory variants within lung cancer GWAS loci by aligning with histone signatures using data of chromatin immunoprecipitation assays followed by sequencing (ChIP-seq) in the Encyclopedia of DNA Elements (ENCODE). A two-stage case–control study with 1,830 cases and 2,001 controls was conducted to explore the associations between candidate variants and NSCLC risk in Chinese population. Bioinformatic annotations and biochemical experiments were performed to further reveal the potential functions of significant variants.

Results

Sixteen potential functional single-nucleotide polymorphisms (SNPs) were selected as candidates from bioinformatics analyses. Three variants out of the 16 candidate SNPs survived after genotyping in stage 1 case–control study, and only the results of SNP rs13064999 were successfully validated in the analyses of stage 2 case–control study. In combined analyses, rs13064999 was significantly associated with NSCLC risk [additive model; odds ratio (OR) = 1.17; 95%CI, 1.07–1.29; p = 0.001]. Functional annotations indicated its potential enhancer bioactivity, and dual-luciferase reporter assays revealed a significant increase in luciferase activity for the reconstructed plasmid with rs13064999 A allele, when compared to the one with wild-type G allele (p_A549 < 0.001, p_SK-MES-1 = 0.004). Further electrophoretic mobility shift assays (EMSA) and super-shift assays confirmed a stronger affinity of HP1γ for the binding motif containing SNP rs13064999 A allele.

Conclusion

These findings suggested that the functional variant rs13064999, identified by the integration of ATAC-seq and ChIP-seq data, contributes to the susceptibility of NSCLC by affecting HP1γ binding, while the exact biological mechanism awaits further exploration.

Profiling of chromatin accessibility identifies transcription factor binding sites across the genome of Aspergillus species

Background

The identification of open chromatin regions and transcription factor binding sites (TFBs) is an important step in understanding the regulation of gene expression in diverse species. ATAC-seq is a technique used for such purpose by providing high-resolution measurements of chromatin accessibility revealed through integration of Tn5 transposase. However, the existence of cell walls in filamentous fungi and associated difficulty in purifying nuclei have precluded the routine application of this technique, leading to a lack of experimentally determined and computationally inferred data on the identity of genome-wide cis-regulatory elements (CREs) and TFBs. In this study, we constructed an ATAC-seq platform suitable for filamentous fungi and generated ATAC-seq libraries of Aspergillus niger and Aspergillus oryzae grown under a variety of conditions.

Results

We applied the ATAC-seq assay for filamentous fungi to delineate the syntenic orthologue and differentially changed chromatin accessibility regions among different Aspergillus species, during different culture conditions, and among specific TF-deleted strains. The syntenic orthologues of accessible regions were responsible for the conservative functions across Aspergillus species, while regions differentially changed between culture conditions and TFs mutants drove differential gene expression programs. Importantly, we suggest criteria to determine TFBs through the analysis of unbalanced cleavage of distinct TF-bound DNA strands by Tn5 transposase. Based on this criterion, we constructed data libraries of the in vivo genomic footprint of A. niger under distinct conditions, and generated a database of novel transcription factor binding motifs through comparison of footprints in TF-deleted strains. Furthermore, we validated the novel TFBs in vivo through an artificial synthetic minimal promoter system.

Conclusions

We characterized the chromatin accessibility regions of filamentous fungi species, and identified a complete TFBs map by ATAC-seq, which provides valuable data for future analyses of transcriptional regulation in filamentous fungi.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01114-0.

Functional characterization of T2D-associated SNP effects on baseline and ER stress-responsive β cell transcriptional activation

Genome-wide association studies (GWAS) have linked single nucleotide polymorphisms (SNPs) at >250 loci in the human genome to type 2 diabetes (T2D) risk. For each locus, identifying the functional variant(s) among multiple SNPs in high linkage disequilibrium is critical to understand molecular mechanisms underlying T2D genetic risk. Using massively parallel reporter assays (MPRA), we test the cis-regulatory effects of SNPs associated with T2D and altered in vivo islet chromatin accessibility in MIN6 β cells under steady state and pathophysiologic endoplasmic reticulum (ER) stress conditions. We identify 1,982/6,621 (29.9%) SNP-containing elements that activate transcription in MIN6 and 879 SNP alleles that modulate MPRA activity. Multiple T2D-associated SNPs alter the activity of short interspersed nuclear element (SINE)-containing elements that are strongly induced by ER stress. We identify 220 functional variants at 104 T2D association signals, narrowing 54 signals to a single candidate SNP. Together, this study identifies elements driving β cell steady state and ER stress-responsive transcriptional activation, nominates causal T2D SNPs, and uncovers potential roles for repetitive elements in β cell transcriptional stress response and T2D genetics.

The Anna Karenina Model of β-Cell Maturation in Development and Their Dedifferentiation in Type 1 and Type 2 Diabetes

Loss of mature β-cell function and identity, or β-cell dedifferentiation, is seen in both type 1 and type 2 diabetes. Two competing models explain β-cell dedifferentiation in diabetes. In the first model, β-cells dedifferentiate in the reverse order of their developmental ontogeny. This model predicts that dedifferentiated β-cells resemble β-cell progenitors. In the second model, β-cell dedifferentiation depends on the type of diabetogenic stress. This model, which we call the "Anna Karenina" model, predicts that in each type of diabetes, β-cells dedifferentiate in their own way, depending on how their mature identity is disrupted by any particular diabetogenic stress. We directly tested the two models using a β-cell-specific lineage-tracing system coupled with RNA sequencing in mice. We constructed a multidimensional map of β-cell transcriptional trajectories during the normal course of β-cell postnatal development and during their dedifferentiation in models of both type 1 diabetes (NOD) and type 2 diabetes (BTBR-Lep^ob/ob ). Using this unbiased approach, we show here that despite some similarities between immature and dedifferentiated β-cells, β-cell dedifferentiation in the two mouse models is not a reversal of developmental ontogeny and is different between different types of diabetes.

Single cell transcriptomics reveal trans-differentiation of pancreatic beta cells following inactivation of the TFIID subunit Taf4

Regulation of gene expression involves a complex and dynamic dialogue between transcription factors, chromatin remodelling and modification complexes and the basal transcription machinery. To address the function of the Taf4 subunit of general transcription factor TFIID in the regulation of insulin signalling, it was inactivated in adult murine pancreatic beta cells. Taf4 inactivation impacted the expression of critical genes involved in beta-cell function leading to increased glycaemia, lowered plasma insulin levels and defective glucose-stimulated insulin secretion. One week after Taf4-loss, single-cell RNA-seq revealed cells with mixed beta cell, alpha and/or delta cell identities as well as a beta cell population trans-differentiating into alpha-like cells. Computational analysis of single-cell RNA-seq defines how known critical beta cell and alpha cell determinants may act in combination with additional transcription factors and the NuRF chromatin remodelling complex to promote beta cell trans-differentiation.

Using high-throughput multi-omics data to investigate structural balance in elementary gene regulatory network motifs

MOTIVATION

The simultaneous availability of ATAC-seq and RNA-seq experiments allows to obtain a more in-depth knowledge on the regulatory mechanisms occurring in gene regulatory networks. In this article, we highlight and analyze two novel aspects that leverage on the possibility of pairing RNA-seq and ATAC-seq data. Namely we investigate the causality of the relationships between transcription factors, chromatin and target genes and the internal consistency between the two omics, here measured in terms of structural balance in the sample correlations along elementary length-3 cycles.

RESULTS

We propose a framework that uses the a priori knowledge on the data to infer elementary causal regulatory motifs (namely chains and forks) in the network. It is based on the notions of conditional independence and partial correlation, and can be applied to both longitudinal and non-longitudinal data. Our analysis highlights a strong connection between the causal regulatory motifs that are selected by the data and the structural balance of the underlying sample correlation graphs: strikingly, >97% of the selected regulatory motifs belong to a balanced subgraph. This result shows that internal consistency, as measured by structural balance, is close to a necessary condition for 3-node regulatory motifs to satisfy causality rules.

AVAILABILITY AND IMPLEMENTATION

The analysis was carried out in MATLAB and the code can be found at https://github.com/albertozenere/Multi-omics-elementary-regulatory-motifs.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

AIAP: A Quality Control and Integrative Analysis Package to Improve ATAC-seq Data Analysis

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a technique widely used to investigate genome-wide chromatin accessibility. The recently published Omni-ATAC-seq protocol substantially improves the signal/noise ratio and reduces the input cell number. High-quality data are critical to ensure accurate analysis. Several tools have been developed for assessing sequencing quality and insertion size distribution for ATAC-seq data; however, key quality control (QC) metrics have not yet been established to accurately determine the quality of ATAC-seq data. Here, we optimized the analysis strategy for ATAC-seq and defined a series of QC metrics for ATAC-seq data, including reads under peak ratio (RUPr), background (BG), promoter enrichment (ProEn), subsampling enrichment (SubEn), and other measurements. We incorporated these QC tests into our recently developed ATAC-seq Integrative Analysis Package (AIAP) to provide a complete ATAC-seq analysis system, including quality assurance, improved peak calling, and downstream differential analysis. We demonstrated a significant improvement of sensitivity (20%–60%) in both peak calling and differential analysis by processing paired-end ATAC-seq datasets using AIAP. AIAP is compiled into Docker/Singularity, and it can be executed by one command line to generate a comprehensive QC report. We used ENCODE ATAC-seq data to benchmark and generate QC recommendations, and developed qATACViewer for the user-friendly interaction with the QC report. The software, source code, and documentation of AIAP are freely available at https://github.com/Zhang-lab/ATAC-seq_QC_analysis.

Successful ATAC-Seq From Snap-Frozen Equine Tissues

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

Runx1 shapes the chromatin landscape via a cascade of direct and indirect targets

Runt-related transcription factor 1 (Runx1) can act as both an activator and a repressor. Here we show that CRISPR-mediated deletion of Runx1 in mouse metanephric mesenchyme-derived mK4 cells results in large-scale genome-wide changes to chromatin accessibility and gene expression. Open chromatin regions near down-regulated loci enriched for Runx sites in mK4 cells lose chromatin accessibility in Runx1 knockout cells, despite remaining Runx2-bound. Unexpectedly, regions near upregulated genes are depleted of Runx sites and are instead enriched for Zeb transcription factor binding sites. Re-expressing Zeb2 in Runx1 knockout cells restores suppression, and CRISPR mediated deletion of Zeb1 and Zeb2 phenocopies the gained expression and chromatin accessibility changes seen in Runx1KO due in part to subsequent activation of factors like Grhl2. These data confirm that Runx1 activity is uniquely needed to maintain open chromatin at many loci, and demonstrate that Zeb proteins are required and sufficient to maintain Runx1-dependent genome-scale repression.

Runt-related transcription factor (Runx) 1 & 2 impact development and disease by activating or repressing transcription. In this manuscript we used genome editing tools to remove Runx1, and as expected, observed widespread changes in chromatin accessibility. Newly closed areas contained Runx1 binding sites and were enriched near genes whose expression depended on Runx1. Interestingly, this occurred despite continued binding of Runx2 to the same regions of DNA, which suggests that Runx2 is insufficient to maintain open chromatin and expression of Runx1 target genes in this cellular context. By contrast, newly opened chromatin regions, many near genes that were upregulated in Runx1 knockout cells, did not enrich for Runx1 binding sites. Instead, these regions were enriched for sites for the repressor Zeb proteins. We found that the loss of Zeb 1 & 2 expression, direct transcriptional targets of Runx1, resulted in the opening of chromatin and upregulation of genes residing near the newly open sites in Runx1 knockout cells. The same sites were also open and nearby genes expressed in edited Zeb1 and Zeb2 knockout cells. Among them were transcription factors, such as the Grhl2 gene, which in turn bind to and upregulate their target genes. Thus, the loss of a single transcription factor initiates a cascade of direct and indirect ramifications with likely negative effects on development and health.

Age-dependent transition from islet insulin hypersecretion to hyposecretion in mice with the long QT-syndrome loss-of-function mutation Kcnq1-A340V

Loss-of-function (LoF) mutations in KCNQ1, encoding the voltage-gated K⁺ channel K_v7.1, lead to long QT syndrome 1 (LQT1). LQT1 patients also present with post-prandial hyperinsulinemia and hypoglycaemia. In contrast, KCNQ1 polymorphisms are associated with diabetes, and LQTS patients have a higher prevalence of diabetes. We developed a mouse model with a LoF Kcnq1 mutation using CRISPR-Cas9 and hypothesized that this mouse model would display QT prolongation, increased glucose-stimulated insulin secretion and allow for interrogation of K_v7.1 function in islets. Mice were characterized by electrocardiography and oral glucose tolerance tests. Ex vivo, islet glucose-induced insulin release was measured, and beta-cell area quantified by immunohistochemistry. Homozygous mice had QT prolongation. Ex vivo, glucose-stimulated insulin release was increased in islets from homozygous mice at 12–14 weeks, while beta-cell area was reduced. Non-fasting blood glucose levels were decreased at this age. In follow-up studies 8–10 weeks later, beta-cell area was similar in all groups, while glucose-stimulated insulin secretion was now reduced in islets from hetero- and homozygous mice. Non-fasting blood glucose levels had normalized. These data suggest that K_v7.1 dysfunction is involved in a transition from hyper- to hyposecretion of insulin, potentially explaining the association with both hypoglycemia and hyperglycemia in LQT1 patients.

Pancreatic cell fate specification: insights into developmental mechanisms and their application for lineage reprogramming

Diabetes is a group of metabolic disorders, which results from insufficient functional pancreatic β-cell mass either due to the autoimmune destruction of insulin producing β-cells, or their death or de-differentiation as compensation for insulin resistance. The ability to reprogram cell types within close developmental proximity to β-cells offers a strategy to replenish β-cell mass and a future possible treatment of diabetes. Here, we review recent advances in the fields of pancreas development and lineage reprogramming. We also probe the possibility of using reprogrammed cells as an approach by which to further understand developmental mechanisms, in particular roadblocks to changing cell identity. Finally, we highlight fundamental challenges that need to be overcome to advance lineage reprogramming for generating pancreatic cells.

The landscape of chromatin accessibility in skeletal muscle during embryonic development in pigs

BACKGROUND

The development of skeletal muscle in pigs during the embryonic stage is precisely regulated by transcriptional mechanisms, which depend on chromatin accessibility. However, how chromatin accessibility plays a regulatory role during embryonic skeletal muscle development in pigs has not been reported. To gain insight into the landscape of chromatin accessibility and the associated genome-wide transcriptome during embryonic muscle development, we performed ATAC-seq and RNA-seq analyses of skeletal muscle from pig embryos at 45, 70 and 100 days post coitus (dpc).

RESULTS

In total, 21,638, 35,447 and 60,181 unique regions (or peaks) were found across the embryos at 45 dpc (LW45), 70 dpc (LW70) and 100 dpc (LW100), respectively. More than 91% of the peaks were annotated within - 1 kb to 100 bp of transcription start sites (TSSs). First, widespread increases in specific accessible chromatin regions (ACRs) from embryos at 45 to 100 dpc suggested that the regulatory mechanisms became increasingly complicated during embryonic development. Second, the findings from integrated ATAC-seq and RNA-seq analyses showed that not only the numbers but also the intensities of ACRs could control the expression of associated genes. Moreover, the motif screening of stage-specific ACRs revealed some transcription factors that regulate muscle development-related genes, such as MyoG, Mef2c, and Mef2d. Several potential transcriptional repressors, including E2F6, OTX2 and CTCF, were identified among the genes that exhibited different regulation trends between the ATAC-seq and RNA-seq data.

CONCLUSIONS

This work indicates that chromatin accessibility plays an important regulatory role in the embryonic muscle development of pigs and regulates the temporal and spatial expression patterns of key genes in muscle development by influencing the binding of transcription factors. Our results contribute to a better understanding of the regulatory dynamics of genes involved in pig embryonic skeletal muscle development.

PRDM1 decreases sensitivity of human NK cells to IL2-induced cell expansion by directly repressing CD25 (IL2RA)

IL2 receptor signaling is crucial for human NK cell activation and gain of effector functions. The molecular mechanisms involved in termination of IL2 activation are largely unknown in human NK cells. PR/SET domain 1 was previously reported to decrease cell growth and increase apoptosis in an IL2-dependent manner in malignant NK cell lines, suggesting the possibility of down-regulation of IL2 signaling pathway gene(s) through direct transcriptional repression. Using ChIP-Seq, we identified a PRDM1 binding site on the first intron of CD25 (IL2RA), which codes for the IL2 receptor subunit regulating sensitivity to IL2 signaling, in primary NK cells activated with engineered K562 cells or IL2. Ectopic expression of PRDM1 down-regulated CD25 expression at transcript and protein levels in two PRDM1 nonexpressing NK cell lines. shRNA-mediated knockdown of CD25 in two malignant NK cell lines led to progressive depletion of NK cells in low IL2 concentrations. By contrast, ectopic CD25 expression in primary human NK cells led to progressive increase in cell number in CD25-transduced cells in low IL2 concentrations. Altogether these results reveal a pivotal role of PRDM1 in inhibition of IL2-induced NK cell expansion through direct repression of CD25 in activated human NK cells. These observations provide additional support for the role of PRDM1 in attenuation of NK cell activation and growth, with implications on neoplastic transformation or NK cell function when it is deregulated.

Scrt1, a transcriptional regulator of β-cell proliferation identified by differential chromatin accessibility during islet maturation

Glucose-induced insulin secretion, a hallmark of mature β-cells, is achieved after birth and is preceded by a phase of intense proliferation. These events occurring in the neonatal period are decisive for establishing an appropriate functional β-cell mass that provides the required insulin throughout life. However, key regulators of gene expression involved in functional maturation of β-cells remain to be elucidated. Here, we addressed this issue by mapping open chromatin regions in newborn versus adult rat islets using the ATAC-seq assay. We obtained a genome-wide picture of chromatin accessible sites (~ 100,000) among which 20% were differentially accessible during maturation. An enrichment analysis of transcription factor binding sites identified a group of transcription factors that could explain these changes. Among them, Scrt1 was found to act as a transcriptional repressor and to control β-cell proliferation. Interestingly, Scrt1 expression was controlled by the transcriptional repressor RE-1 silencing transcription factor (REST) and was increased in an in vitro reprogramming system of pancreatic exocrine cells to β-like cells. Overall, this study led to the identification of several known and unforeseen key transcriptional events occurring during β-cell maturation. These findings will help defining new strategies to induce the functional maturation of surrogate insulin-producing cells.