Gene Expression Analysis of the Pre-Diabetic Pancreas to Identify Pathogenic Mechanisms and Biomarkers of Type 1 Diabetes

Predictive biomarkers of rapidly developing insulin deficiency in children with type 1 diabetes

INTRODUCTION

The rate of progression to complete insulin deficiency varies greatly in type 1 diabetes. This constitutes a challenge, especially when randomizing patients in intervention trials aiming to preserve beta cell function. This study aimed to identify biomarkers predictive of either a rapid or slow disease progression in children with new-onset type 1 diabetes.

RESEARCH DESIGN AND METHODS

A retrospective, longitudinal cohort study of children (<18 years) with type 1 diabetes (N=46) was included at diagnosis and followed until complete insulinopenia (C-peptide <0.03 nmol/L). Children were grouped into rapid progressors (n=20, loss within 30 months) and slow progressors (n=26). A sex-matched control group of healthy children (N=45) of similar age was included for comparison. Multiple biomarkers were assessed by proximity extension assay (PEA) at baseline and follow-up.

RESULTS

At baseline, rapid progressors had lower C-peptide and higher autoantibody levels than slow. Three biomarkers were higher in the rapid group: carbonic anhydrase 9, corticosteroid 11-beta-dehydrogenase isozyme 1, and tumor necrosis factor receptor superfamily member 21. In a linear mixed model, 25 proteins changed over time, irrespective of group. One protein, a coxsackievirus B-adenovirus receptor (CAR) increased over time in rapid progressors. Eighty-one proteins differed between type 1 diabetes and healthy controls. Principal component analysis could not distinguish between rapid, slow, and healthy controls.

CONCLUSIONS

Despite differences in individual proteins, the combination of multiple biomarkers analyzed by PEA could not distinguish the rate of progression in children with new-onset type 1 diabetes. Only one marker was altered significantly when considering both time and group effects, namely CAR, which increased significantly over time in the rapid group. Nevertheless, we did find some markers that may be useful in predicting the decline of the C-peptide. Moreover, these could potentially be important for understanding type 1 diabetes pathogenesis.

GWAS for systemic sclerosis identifies six novel susceptibility loci including one in the Fcγ receptor region

Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

Lactiplantibacillus plantarum and Saussurea costus as Therapeutic Agents against a Diabetic Rat Model-Approaches to Investigate Pharmacophore Modeling of Human IkB Kinase and Molecular Interaction with Dehydrocostus Lactone of Saussurea costus

Lactic acid bacteria is well-known as a vital strategy to alleviate or prevent diabetes. Similarly, the plant Saussurea costus (Falc) Lipsch is a preventive power against diabetes. Here, we aimed to determine whether lactic acid bacteria or Saussurea costus is more effective in treating a diabetic rat model in a comparative study manner. An in vivo experiment was conducted to test the therapeutic activity of Lactiplantibacillus plantarum (MW719476.1) and S. costus plants against an alloxan-induced diabetic rat model. Molecular, biochemical, and histological analyses were investigated to evaluate the therapeutic characteristics of different treatments. The high dose of S. costus revealed the best downregulated expression for the IKBKB, IKBKG, NfkB1, IL-17A, IL-6, IL-17F, IL-1β, TNF-α, TRAF6, and MAPK genes compared to Lactiplantibacillus plantarum and the control groups. The downregulation of IKBKB by S. costus could be attributed to dehydrocostus lactone as an active compound with proposed antidiabetic activity. So, we performed another pharmacophore modeling analysis to test the possible interaction between human IkB kinase beta protein and dehydrocostus lactone as an antidiabetic drug. Molecular docking and MD simulation data confirmed the interaction between human IkB kinase beta protein and dehydrocostus lactone as a possible drug. The target genes are important in regulating type 2 diabetes mellitus signaling, lipid and atherosclerosis signaling, NF-κB signaling, and IL-17 signaling pathways. In conclusion, the S. costus plant could be a promising source of novel therapeutic agents for treating diabetes and its complications. Dehydrocostus lactone caused the ameliorative effect of S. costus by its interaction with human IkB kinase beta protein. Further, future studies could be conducted to find the clinical efficacy of dehydrocostus lactone.

An Early Islet Transcriptional Signature Is Associated With Local Inflammation in Autoimmune Diabetes

Identifying the early islet cellular processes of autoimmune type 1 diabetes (T1D) in humans is challenging given the absence of symptoms during this period and the inaccessibility of the pancreas for sampling. In this article, we study temporal events in pancreatic islets in LEW.1WR1 rats, in which autoimmune diabetes can be induced with virus infection, by performing transcriptional analysis of islets harvested during the prediabetic period. Single-cell RNA-sequencing and differential expression analyses of islets from prediabetic rats reveal subsets of β- and α-cells under stress as evidenced by heightened expression, over time, of a transcriptional signature characterized by interferon-stimulated genes, chemokines including Cxcl10, major histocompatibility class I, and genes for the ubiquitin-proteasome system. Mononuclear phagocytes show increased expression of inflammatory markers. RNA-in situ hybridization of rat pancreatic tissue defines the spatial distribution of Cxcl10+ β- and α-cells and their association with CD8+ T cell infiltration, a hallmark of insulitis and islet destruction. Our studies define early islet transcriptional events during immune cell recruitment to islets and reveal spatial associations between stressed β- and α-cells and immune cells. Insights into such early processes can assist in the development of therapeutic and prevention strategies for T1D.

A discovery-based proteomics approach identifies protein disulphide isomerase (PDIA1) as a biomarker of β cell stress in type 1 diabetes

BACKGROUND

Stress responses within the β cell have been linked with both increased β cell death and accelerated immune activation in type 1 diabetes (T1D). At present, information on the timing and scope of these responses as well as disease-related changes in islet β cell protein expression during T1D development is lacking.

METHODS

Data independent acquisition-mass spectrometry was performed on islets collected longitudinally from NOD mice and NOD-SCID mice rendered diabetic through T cell adoptive transfer.

FINDINGS

In islets collected from female NOD mice at 10, 12, and 14 weeks of age, we found a time-restricted upregulation of proteins involved in stress mitigation and maintenance of β cell function, followed by loss of expression of protective proteins that heralded diabetes onset. EIF2 signalling and the unfolded protein response, mTOR signalling, mitochondrial function, and oxidative phosphorylation were commonly modulated pathways in both NOD mice and NOD-SCID mice rendered acutely diabetic by T cell adoptive transfer. Protein disulphide isomerase A1 (PDIA1) was upregulated in NOD islets and pancreatic sections from human organ donors with autoantibody positivity or T1D. Moreover, PDIA1 plasma levels were increased in pre-diabetic NOD mice and in the serum of children with recent-onset T1D compared to non-diabetic controls.

INTERPRETATION

We identified a core set of modulated pathways across distinct mouse models of T1D and identified PDIA1 as a potential human biomarker of β cell stress in T1D.

FUNDING

NIH (R01DK093954, DK127308, U01DK127786, UC4DK104166, R01DK060581, R01GM118470, and 5T32DK101001-09). VA Merit Award I01BX001733. JDRF (2-SRA-2019-834-S-B, 2-SRA-2018-493-A-B, 3-PDF-20016-199-A-N, 5-CDA-2022-1176-A-N, and 3-PDF-2017-385-A-N).

The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.

Effect of Inulin Intervention on Metabolic Control and Methylation of INS and IRS1 Genes in Patients with Type 2 Diabetes Mellitus

BACKGROUND AND AIMS

Currently, treatments are being sought to improve the control of type II diabetes mellitus (T2DM), and inulin has been shown to be effective in reducing glucose levels and other metabolic control parameters. These effects on metabolic control may be associated with changes in the epigenetic modulation of genes of the insulin pathway. Therefore, our objective is to determine the effect of agave inulin in metabolic control parameters and in INS and IRS1 genes' methylation in T2DM patients.

METHODS

This was a longitudinal experimental study with 67 Mexican participants who received an intervention of inulin agave (10 g daily) for 2 months. The methylation of the INS and IRS1 genes was determined by MSP.

RESULTS

For the INS gene, we found a significant decrease in the proportions of T2DM patients with methylated DNA after inulin intervention (p = 0.0001). In contrast, the difference in the proportions of the unmethylated IRS1 gene before and after the inulin intervention was not significant (p = 0.79). On the other hand, we observed changes in the number of T2DM patients' recommended categories for metabolic control depending on the methylation of INS and IRS1 genes before and after treatment with inulin.

CONCLUSION

For the first time, we report the modification in the methylation of two genes, INS and IRS1, of the insulin pathway and provide information on the possible relevant role of epigenetics as a key factor in positive changes in metabolic control parameters by inulin intake in T2DM patients.

Retinal Transcriptome and Cellular Landscape in Relation to the Progression of Diabetic Retinopathy

Purpose

Previous studies that identify putative genes associated with diabetic retinopathy are only focusing on specific clinical stages, thus resulting genes are not necessarily reflective of disease progression.

This study identified genes associated with the severity level of diabetic retinopathy using the likelihood-ratio test (LRT) and ordinal logistic regression (OLR) model, as well as to profile immune and retinal cell landscape in progressive diabetic retinopathy using a machine learning deconvolution approach.

Methods

This study used a published transcriptomic dataset (GSE160306) from macular regions of donors with different degrees of diabetic retinopathy (10 healthy controls, 10 cases of diabetes, 9 cases of nonproliferative diabetic retinopathy, and 10 cases of proliferative diabetic retinopathy or combined with diabetic macular edema). LRT and OLR models were applied to identify severity-associated genes. In addition, CIBERSORTx was used to estimate proportional changes of immune and retinal cells in progressive diabetic retinopathy.

Results

By controlling for gender and age using LRT and OLR, 50 genes were identified to be significantly increased in expression with the severity of diabetic retinopathy. Functional enrichment analyses suggested these severity-associated genes are related to inflammation and immune responses. CCND1 and FCGR2B are further identified as key regulators to interact with many other severity-associated genes and are crucial to inflammation. Deconvolution analyses demonstrated that the proportions of memory B cells, M2 macrophages, and Müller glia were significantly increased with the progression of diabetic retinopathy.

Conclusions

These findings demonstrate that deep analyses of transcriptomic data can advance our understanding of progressive ocular diseases, such as diabetic retinopathy, by applying LRT and OLR models as well as bulk gene expression deconvolution.

Highly perturbed genes and hub genes associated with type 2 diabetes in different tissues of adult humans: a bioinformatics analytic workflow

Type 2 diabetes (T2D) has a complex etiology which is not yet fully elucidated. The identification of gene perturbations and hub genes of T2D may deepen our understanding of its genetic basis. We aimed to identify highly perturbed genes and hub genes associated with T2D via an extensive bioinformatics analytic workflow consisting of five steps: systematic review of Gene Expression Omnibus and associated literature; identification and classification of differentially expressed genes (DEGs); identification of highly perturbed genes via meta-analysis; identification of hub genes via network analysis; and downstream analysis of highly perturbed genes and hub genes. Three meta-analytic strategies, random effects model, vote-counting approach, and p value combining approach, were applied. Hub genes were defined as those nodes having above-average betweenness, closeness, and degree in the network. Downstream analyses included gene ontologies, Kyoto Encyclopedia of Genes and Genomes pathways, metabolomics, COVID-19-related gene sets, and Genotype-Tissue Expression profiles. Analysis of 27 eligible microarrays identified 6284 DEGs (4592 downregulated and 1692 upregulated) in four tissue types. Tissue-specific gene expression was significantly greater than tissue non-specific (shared) gene expression. Analyses revealed 79 highly perturbed genes and 28 hub genes. Downstream analyses identified enrichments of shared genes with certain other diabetes phenotypes; insulin synthesis and action-related pathways and metabolomics; mechanistic associations with apoptosis and immunity-related pathways; COVID-19-related gene sets; and cell types demonstrating over- and under-expression of marker genes of T2D. Our approach provided valuable insights on T2D pathogenesis and pathophysiological manifestations. Broader utility of this pipeline beyond T2D is envisaged.

Identification of Novel Disease-Relevant Genes and Pathways in the Pathogenesis of Type 1 Diabetes: A Potential Defect in Pancreatic Iron Homeostasis

Multiple pathways contribute to the pathophysiological development of type 1 diabetes (T1D); however, the exact mechanisms involved are unclear. We performed differential gene expression analysis in pancreatic islets of NOD mice versus age-matched congenic NOD.B10 controls to identify genes that may contribute to disease pathogenesis. Novel genes related to extracellular matrix development and glucagon and insulin signaling/secretion were changed in NOD mice during early inflammation. During "respective" insulitis, the expression of genes encoding multiple chemosensory olfactory receptors were upregulated, and during "destructive" insulitis, the expression of genes involved in antimicrobial defense and iron homeostasis were downregulated. Islet inflammation reduced the expression of Hamp that encodes hepcidin. Hepcidin is expressed in β-cells and serves as the key regulator of iron homeostasis. We showed that Hamp and hepcidin levels were lower, while iron levels were higher in the pancreas of 12-week-old NOD versus NOD.B10 mice, suggesting that a loss of iron homeostasis may occur in the islets during the onset of "destructive" insulitis. Interestingly, we showed that the severity of NOD disease correlates with dietary iron intake. NOD mice maintained on low-iron diets had a lower incidence of hyperglycemia, while those maintained on high-iron diets had an earlier onset and higher incidence of disease, suggesting that high iron exposure combined with a loss of pancreatic iron homeostasis may exacerbate NOD disease. This mechanism may explain the link seen between high iron exposure and the increased risk for T1D in humans.

Honey polyphenolic fraction inhibits cyclooxygenase-2 expression via upregulation of microRNA-26a-5p expression in pancreatic islets

Objectives

Honey total polyphenolic fraction (HTPF) is reported to have anti-disease potential, however the role of HTPF in the regulation of microRNAs (miRNAs) has never been investigated. This study was undertaken to investigate the potential of HTPF against inflammation via regulation of miRNAs in pancreatic islets of Langerhans.

Methods

Pancreatic islets were isolated from C57BL/6 mice and HTPF was purified from honey. Bioinformatics algorithms were used to determine miRNA target genes. Expression of miRNA and mRNA was determined using their specific taqman assays. Pairing between miRNA and 3′ untranslated region (3′UTR) of mRNA was confirmed using luciferase reporter clone containing the 3′UTR of mRNA sequences and results were verified by transfection of mouse pancreatic β-cell line Min6 with miRNA inhibitors.

Results

The data showed that mmu-miR-26a-5p is a direct regulator of cyclooxygenase-2 (COX-2) expression and HTPF inhibits COX-2 expression or prostaglandin E2 (PGE2) production via up-regulating mmu-miR-26a-5p expression. Transfection of islets with anti-miR-26a-5p significantly enhanced COX-2 expression and PGE2 production ( p < .01), while HTPF treatment significantly inhibited anti-miR-26a-5p transfection-induced COX-2 expression or PGE2 production ( p < .05). These findings were further verified in pancreatic β-cells Min6. Moreover, the data also determined that HTPF also inhibits glucose-induced nuclear transcription factor (NF)-κB activity.

Conclusion

HTPF suppresses glucose-induced PGE2 production and activation of NF-κB via negative regulation of COX-2 and mmu-miR26a-5p. These novel pharmacological actions of HTPF on glucose-stimulated pancreatic islets provide new suggestions that HTPF or HTPF-derived compounds inhibit glucose induced inflammation in pancreas by up-regulating the expression of microRNAs.

Genome-Wide Identification of N6-Methyladenosine Associated SNPs as Potential Functional Variants for Type 1 Diabetes

OBJECTIVES

N6-methyladenosine (m6A) is essential in the regulation of the immune system, but the role that its single nucleotide polymorphisms (SNPs) play in the pathogenesis of type 1 diabetes (T1D) remains unknown. This study demonstrated the association between genetic variants in m6A regulators and T1D risk based on a case-control study in a Chinese population.

METHODS

The tagging SNPs in m6A regulators were genotyped in 1005 autoantibody-positive patients with T1D and 1257 controls using the Illumina Human OmniZhongHua-8 platform. Islet-specific autoantibodies were examined by radioimmunoprecipitation in all the patients. The mixed-meal glucose tolerance test was performed on 355 newly diagnosed patients to evaluate their residual islet function. The functional annotations for the identified SNPs were performed in silico. Using 102 samples from a whole-genome expression microarray, key signaling pathways associated with m6A regulators in T1D were comprehendingly evaluated.

RESULTS

Under the additive model, we observed three tag SNPs in the noncoding region of the PRRC2A (rs2260051, rs3130623) and YTHDC2 (rs1862315) gene are associated with T1D risk. Although no association was found between these SNPs and islet function, patients carrying risk variants had a higher positive rate for ZnT8A, GADA, and IA-2A. Further analyses showed that rs2260051[T] was associated with increased expression of PRRC2A mRNA (P = 7.0E-13), and PRRC2A mRNA was significantly higher in peripheral blood mononuclear cell samples from patients with T1D compared to normal samples (P = 0.022). Enrichment analyses indicated that increased PRRC2A expression engages in the most significant hallmarks of cytokine-cytokine receptor interaction, cell adhesion and chemotaxis, and neurotransmitter regulation pathways. The potential role of increased PRRC2A in disrupting immune homeostasis is through the PI3K/AKT pathway and neuro-immune interactions.

CONCLUSION

This study found intronic variants in PRRC2A and YTHDC2 associated with T1D risk in a Chinese Han population. PRRC2A rs2260051[T] may be implicated in unbalanced immune homeostasis by affecting the expression of PRRC2A mRNA. These findings enriched our understanding of m6A regulators and their intronic SNPs that underlie the pathogenesis of T1D.

Novel insights into immunohistochemical analysis for diagnosing serous neoplasm of the pancreas: aquaporin 1, stereocilin, and transmembrane protein 255B

AIMS

Serous (cystic) neoplasm (SCN) of the pancreas is generally benign, and surgical treatment is recommended in only a limited number of cases. To avoid unnecessary surgery, an accurate diagnosis of SCN is essential. In the present study, we aimed to identify new immunohistochemical markers with which to distinguish SCN from other tumours.

METHODS AND RESULTS

We compared the comprehensive gene expression profiles of SCN with those of normal pancreas and pancreatic ductal adenocarcinoma (PDAC). We selected the candidate molecules that were up-regulated in SCN, were minimally expressed or unexpressed in PDAC, and had specific and available antibodies suitable for immunohistochemistry, and then analysed their immunohistochemical expression in various tumours. We selected aquaporin 1 (AQP1), stereocilin (STRC), fibroblast growth factor receptor 3 (FGFR3), and transmembrane protein 255B (TMEM255B), which were diffusely expressed in SCN cells in 79%, 100%, 100% and 100% of SCN cases. AQP1 was not expressed in other tumours, except in 20% of mucinous cystic neoplasms (MCNs) and 19% of PDACs. STRC was rarely expressed in MCNs, neuroendocrine neoplasms (NENs), and PDACs. FGFR3 was expressed in 31% of intraductal papillary mucinous neoplasms (IPMNs), 50% of intraductal oncocytic papillary neoplasms, 40% of NENs, 30% of acinar cell carcinomas, 40% of solid pseudopapillary neoplasms, and 52% of PDACs. TMEM255B was not expressed in the other tumours, except in 50% of MCNs, 80% of gastric-subtype IPMNs, and 29% of PDACs. All antigens were usually expressed in a small proportion of cells when they were positive in tumours other than SCN.

CONCLUSIONS

These findings indicate that AQP1 and STRC, and potentially TMEM255B, may act as SCN markers.