Abstract
Although type 1 diabetes (T1D) is primarily a disease of the pancreatic beta-cells, understanding of the disease-associated alterations in the whole pancreas could be important for the improved treatment or the prevention of the disease. We have characterized the whole-pancreas gene expression of patients with recently diagnosed T1D from the Diabetes Virus Detection (DiViD) study and non-diabetic controls. Furthermore, another parallel dataset of the whole pancreas and an additional dataset from the laser-captured pancreatic islets of the DiViD patients and non-diabetic organ donors were analyzed together with the original dataset to confirm the results and to get further insights into the potential disease-associated differences between the exocrine and the endocrine pancreas. First, higher expression of the core acinar cell genes, encoding for digestive enzymes, was detected in the whole pancreas of the DiViD patients when compared to non-diabetic controls. Second, In the pancreatic islets, upregulation of immune and inflammation related genes was observed in the DiViD patients when compared to non-diabetic controls, in line with earlier publications, while an opposite trend was observed for several immune and inflammation related genes at the whole pancreas tissue level. Third, strong downregulation of the regenerating gene family (REG) genes, linked to pancreatic islet growth and regeneration, was observed in the exocrine acinar cell dominated whole-pancreas data of the DiViD patients when compared with the non-diabetic controls. Fourth, analysis of unique features in the transcriptomes of each DiViD patient compared with the other DiViD patients, revealed elevated expression of central antiviral immune response genes in the whole-pancreas samples, but not in the pancreatic islets, of one DiViD patient. This difference in the extent of antiviral gene expression suggests different statuses of infection in the pancreas at the time of sampling between the DiViD patients, who were all enterovirus VP1+ in the islets by immunohistochemistry based on earlier studies. The observed features, indicating differences in the function, status and interplay between the exocrine and the endocrine pancreas of recent onset T1D patients, highlight the importance of studying both compartments for better understanding of the molecular mechanisms of T1D.
Introduction
Type 1 diabetes (T1D) is an autoimmune disease characterized by the progressing destruction of pancreatic β-cells and the related loss of insulin production, blood glucose control and finally extreme hyperglycemia and ketoacidosis in untreated patients (–). The exact cause of T1D remains unknown with both genetic (, ) and environmental factors suspected to play a role in the disease etiology. The disease progression is assumed to be initiated by inflammation of the pancreas and the islets of Langerhans (pancreatic islets) containing β-cells.
In the Diabetes Virus Detection (DiViD) study, samples from the pancreatic tissue of six living adults with recently diagnosed T1D were collected using pancreatic tail resection by laparoscopy (). The insulin secretion capability of the isolated pancreatic islets was shown to be recoverable when removed from the diabetic environment, and all transcripts in the human insulin pathway were present, but less expressed in the DiViD patients when compared to the nondiabetic controls (). To investigate the potential link between virus infections and T1D, indicators of virus infection in the samples have been explored with multiple approaches, and the presence of enteroviruses has been shown in the pancreatic islets of all six DiViD patients ().
T1D is a disease of the whole pancreas, not just the pancreatic islets, although most of the studies into the etiology of T1D have focused on the endocrine tissue of the pancreas. The pancreas of individuals with T1D have been observed to be smaller in size compared to the pancreases of patients with type 2 diabetes (T2D) or control subjects (, ). In addition to decreased volume, other anomalies, including various histological changes in the exocrine pancreas of diabetic patients, have been discovered (). Also pancreatic exocrine dysfunction or pancreatic insufficiency are frequent for T1D patients (, ). Furthermore, preserved beta cell insulin secretion has been connected with preserved acinar cell function and volume (, ). Finally, infiltration of the exocrine pancreas by immune cells in connection with T1D has been noticed in several studies (, ).
There has been an increasing body of evidence linking viral infections to the occurrence or facilitation of T1D (). While multiple types of different viruses have been associated with T1D (, ), the common enteroviruses, especially the group B Coxsackieviruses, have been proposed as strong candidates for promoting T1D (, ). Enteroviruses possess a strong affinity for the pancreatic islets (), and the presence of low grade enterovirus infection based on the detection of enterovirus capsid protein VP1 in the pancreatic islets has been demonstrated in several studies, including the DiViD study (–). Compared to an acute infection, in persistent infection the virus may be present at extremely low numbers and the replication of the virus might be decreased, rendering the virus difficult to detect (). Recently, studies utilizing highly sensitive fluorescent in situ hybridization method to detect enteroviral RNA in pancreas tissue samples, have shown the presence of the virus both in the islets and in the exocrine pancreas of T1D patients and at-risk individuals positive for T1D autoantibodies (, ). The progression of enterovirus infection and the possible facilitation of T1D by the infection are likely influenced by the virus itself, as well as host responses and host factors (e.g. gender, age, dietary factors, infection timing, genetics) (). Generally, viral agents have been observed to be able to regulate host responses, and to downregulate the inflammatory responses in human cells (). There are also clear differences in the immunogenicity of different viruses, including different CVB1 strains (–). In addition, several cell surface receptors for enteroviruses are recognized, and the receptor usage seems to differ to some extent between distinct members of the virus family (). This will likely result in differences in the cell populations that are most susceptible for infection as well as the outcomes of infection.
To investigate the role of both the exocrine and the endocrine pancreas in early T1D, we characterized total gene expression in the whole pancreas tissue of the DiViD study patients with recently diagnosed T1D and compared it with the pancreas tissue transcriptomes of non-diabetic organ donors. To further characterize the contribution of different pancreas compartments on the observed changes, gene expression in the whole pancreas was compared to the gene expression in the pancreatic islets of the DiViD patients. Additionally, a gene expression dataset produced by an independent laboratory from pancreas whole-tissue samples of the same DiViD patients was used to verify part of the results. The combined analysis of these three gene expression datasets allowed for a comprehensive exploration of gene expression and systemic changes in the pancreas related to recently diagnosed T1D. Finally, in addition to investigating gene expression in the exocrine and endocrine pancreas, the generated datasets also enabled the exploration of potential enterovirus related signals in the whole pancreas and in the pancreatic islets of the DiViD patients.
Materials and Methods
Pancreas tissue samples from recent onset T1D patients participating in the Diabetes Virus Detection (DiViD) study () were analyzed. The samples were collected by pancreatic tail resections performed 3–9 weeks after the diagnosis of T1D. The DiViD study was approved by the Norwegian Government´s Regional Ethics Committee, and written informed consent was obtained from all patients.
Whole Pancreas Transcriptome 1
RNA was isolated from pancreas whole tissue samples preserved in RNAlater, using RNeasy Plus Mini kit (Qiagen). Samples from the DiViD cases 2-6, (two females and three males, aged 24-35 years) () were prepared and analyzed. Based on the quality control results, RNA from the pancreas tissue sample of DiViD case 2 was degraded (Supplementary Figure 1), and therefore an adjusted protocol was used for RNA-sequencing analysis of this sample. The adjusted protocol for DiViD case 2 was provided by Illumina TruSeq Stranded Total RNA Reference Guide (Document # 1000000040499 v00), where during library preparation of the sample, the RNA fragmentation time was reduced from 8 min to 6 min. Three commercial human pancreas total RNA preparates were used as controls (#AM7954, Ambion; #540023, Agilent Technologies; #636577, Clontech) (Supplementary Table 1). RNA-sequencing was performed using Illumina TruSeq® Stranded total RNA protocol and Illumina HiSeq 3000 instrument with TruSeq v4 sequencing. Paired-end sequencing with 2 x 150 bp read length with 6 bp index run was used.
Whole Pancreas Transcriptome 2
Independent total RNA-sequencing was performed for separate pancreas whole tissue pieces from the DiViD patients (cases 1-6) in the Oslo University Hospital (). The sequencing was done using Illumina TruSeq® Stranded total RNA protocol and Illumina HiSeq 2000 instrument with paired-end sequencing and 2 x 100 bp read length.
Pancreatic Islet Transcriptome
The laser-capture microarray data was generated for the DiViD cases 2-6 and for 18 Network for Pancreatic Organ donors with Diabetes (nPOD) non-diabetic controls and normalized as described before ().pt?>
Data Analyses
Cell type proportions in the whole pancreas transcriptome 1 dataset were estimated using the online deconvolution tool CIBERSORT () with a signature matrix constructed from pancreatic single-cell data () under the accession ID E-MTAB-5061 in ArrayExpress. The utilized signature matrix is available as Supplementary Data 1. For exploring differential expression between the DiViD cases and the controls, the reproducibility optimized test statistic (ROTS) (), shown to perform well with RNA-Seq data (), was used. ROTS false discovery rate (FDR) 0.05 and a fold change (FC) threshold of two were used to define the differentially expressed (DE) genes. To explore the enrichment of gene ontology (GO) biological processes (BP) among the DE genes, the Database for Annotation, Visualization and Integrated Discovery (DAVID) () version 6.8 was used. The GO FAT terms were considered, filtering out the broadest GO term annotation categories. Only GO BP terms with FDR ≤ 0.05 were considered as enriched. Protein-protein interactions among the gene products of interest were queried from the STRING functional protein association networks database () using both the predicted and known interactions. Only high confidence interactions (confidence score ≥ 0.7) were considered. To compare gene expression between the DiViD cases, gene expression in each dataset was z-score transformed using the case samples in each dataset only.
To ensure comparability, all datasets were filtered and preprocessed as similarly as possible. For details related to preprocessing and data analysis of the datasets, see Supplementary Methods.
Results
Recent Onset T1D Associated Gene Expression Profile Can be Identified in the Whole Pancreas Tissue
We performed RNA-sequencing for pancreas tissue samples collected from five recent onset T1D patients, as well as for three commercial human pancreas tissue RNA preparates of non-diabetic adults. To get a rough estimation on the contribution of different pancreatic cell types on total tissue transcriptomes of each sample, cell type deconvolution was performed using the CIBERSORT tool (). Acinar cells are clearly the most abundant cell type in the pancreas, which was also reflected in the cell type deconvolution results of the current RNA-sequencing data. While the contribution of acinar cells on the RNA-sequencing data was estimated to be over 90%, the estimated contribution of endocrine cells was only few percent for all the samples with no difference in proportions between the cases and the controls. (Figure 1A)
Figure 1
The pancreas tissue transcriptomes of individuals with recent onset of T1D were clearly separated from the control samples on the first two principal components of the Principal Component Analysis (PCA) plot (Figure 1B). Furthermore, the differential expression analysis of the transcriptomics data revealed T1D-associated differences in the expression of 365 genes (Figure 1C and Supplementary Table 2).
T1D Influences the Expression of Genes in the Exocrine and the Endocrine Pancreas
Study on the cell type enriched genes in the pancreas showed concurrent T1D associated differences in the beta-cell and acinar-cell gene expression (Figure 2A). Also for markers of the pancreatic ductal cells, KRT19 and SPP1, a trend towards higher expression in the control samples was observed, but this difference was not statistically significant. No T1D associated differences were observed in the expression of the alpha-cell and delta-cell markers GCG and SST, respectively.
Figure 2
The loss of insulin-producing beta-cells in the pancreatic islets is the core pathogenic feature in T1D. Also in the current study, INS mRNA levels were lower in every tissue sample from individuals with recent onset T1D when compared with the control samples, although the statistical significance remained slightly above the used FDR threshold of 0.05 (Figure 2A). Furthermore, five other genes shown to be enriched in the pancreatic beta cells by single cell RNA-sequencing (ADCYAP1, G6PC2, NPTX2, PFKFB2, SYT13) (, ) were significantly downregulated in the tissue samples of individuals with recent onset T1D. In addition to impaired insulin production, these results show a broader reduction in pancreatic beta-cell functions shortly after the diagnosis of T1D.
The core function of the acinar cells of the exocrine pancreas is the production of digestive enzymes. Transcripts encoding for several digestive enzymes including AMY2A, AMY2B, CELA3A and CELA2B were significantly upregulated in the pancreas tissue samples of T1D patients when compared with the control samples (Figure 2A). A similar trend, although not statistically significant, was also observed for PNLIP and PRSS1 often used as markers of acinar cells. Segerstolpe et al. have utilized single cell RNA-sequencing data to further define two distinct acinar cell subpopulations: the inflammatory cluster characterized by the increased expression of inflammatory genes and the enzymatic cluster characterized by the elevated levels of key acinar genes encoding for digestive enzymes (). Interestingly, many of the differentially expressed genes between the T1D and control pancreatic samples in the current study were also distinguishing the suggested acinar cell subpopulations from each other. Among the overlapping genes, almost all enzymatic cluster associated genes were upregulated in the pancreas tissue samples of recent onset T1D patients. On the other hand, the overlapping genes associated with the inflammatory acinar cell cluster were downregulated in the pancreas tissue samples of recent onset T1D patients, potentially reflecting a T1D associated change in balance between the different acinar cell subpopulations (Figure 2B). Finally, REG3A has been shown to mark a subpopulation of acinar cells located close to the islets of Langerhans (), and it has been suggested to support pancreatic islet function (). Interestingly, REG3A was the most strongly downregulated gene in the pancreas tissue samples of individuals with recent onset T1D (Figure 2C and Supplementary Table 2). Additionally REG1B and REG3G, also members of the Reg-family genes, were clearly downregulated in the pancreas tissue samples of individuals with recent onset T1D (Figure 2C).
Inflammation and Immune Response Associated Genes Are Downregulated in the Pancreas Tissue of Recent Onset T1D Patients
Based on the functional gene ontology (GO) enrichment analysis of the differentially expressed genes, significant T1D associated differences were observed in the expression of genes associated with extracellular matrix organization, cell migration and cell adhesion (Supplementary Figure 2), potentially reflecting changes in cell-cell contacts and communication. These changes might also influence the ability of enteroviruses and other viruses to enter the cells.
In addition, several genes associated with inflammation and immune responses were differentially expressed in the pancreas tissue samples and pancreatic islets between the DiViD patients and non-diabetic controls (Figure 3A and Supplementary Figure 2, Supplementary Table 3). Interestingly, these differentially expressed genes were almost exclusive for one tissue type. In line with earlier studies (), the upregulation of HLA class I genes and several interferon response genes was observed in the pancreatic islets of the DiViD patients when compared with non-diabetic controls. On the other hand, in the pancreas whole-tissue samples, several genes involved in the early steps of complement activation (e.g. C1qB, C1qC, C1R, C1S, C3), Fc gamma receptors (FCGR2A, FCGR3A), and Tumor necrosis factor receptors (TNFRSF1B, TNFRSF10B, TNFRSF12a) were significantly downregulated in the DiViD patients when compared with non-diabetic controls (Figures 3B, C). Some of these downregulated immune and inflammatory genes, including LCN2, CLDN2, MHC class II genes (HLA-DQB1 and HLA-DRB1) and some complement genes, were also among the genes distinguishing the inflammatory acinar cell cluster from the enzymatic acinar cell cluster in (), supporting the potential contribution of the exocrine pancreas on the changes observed in the pancreas total tissue samples. Also, three transcription factors with central roles in the activation of immune and inflammatory responses (IRF1, IRF8 and STAT3) were significantly downregulated in the whole pancreas of T1D patients. Intriguingly, the most strongly upregulated gene in the pancreas tissue samples of recent onset T1D patients was SIGLEC11 (Figure 3C), which is an immune inhibitory SIGLEC receptor.
Figure 3
Antiviral Immune Response Signature in the Exocrine Pancreas of One T1D Patient
In addition to the T1D associated differences in gene expression in the pancreas, our aim was to understand the differences between the individual recent onset T1D patients in the current study. For that, we studied the top 100 most highly expressed transcripts in each DiViD case relative to any other DiViD case, and named these gene sets as individual signatures. As further validation, similar analysis was performed on an independent RNA-sequencing dataset from distinct pancreas tissue samples collected at the same time from the same individuals. Genes that were included in the individual signatures in both datasets were considered as reliable observations. An individual signature of 45 genes was discovered in both RNA-sequencing datasets for DiViD case 6 (Figure 4A). Only in this signature of DiViD case 6, strong known interactions between the signature transcripts were detected (Figure 4B), indicating a presence of a group of tightly related features. This signature included several interconnected transcripts with central roles in antiviral interferon responses, including EIF2AK2, MX1, STAT1 and ISG15 (Figure 4B), suggesting the presence of viral infection in the pancreas of DiViD case 6. Interestingly, when these results were compared with microarray transcriptomics data collected from the laser capture microdissected islets of the same individuals, no such clear trend was observed (Figure 4C). Furthermore, the expression of the antiviral signature genes in DiViD case 6 was clearly higher than in the control pancreas tissue samples (Supplementary Figure 3). In conclusion, our data indicates strong contribution from the exocrine pancreas, rather than the pancreatic islets, on the antiviral signature observed in the current study, which is specific for DiViD case 6.
Figure 4
Discussion
Although T1D is primarily a disease of the pancreatic beta-cells, understanding of the disease-associated alterations in the whole pancreas could be important for the improved treatment or prevention of the disease. In this study, our aim was to characterize changes in pancreas gene expression that are associated with the early stages of human T1D soon after the diagnosis of the disease. Our pancreas tissue RNA-sequencing results showed concurrent T1D associated differences in the expression of both beta-cell and acinar cell linked transcripts.
T1D has been associated with reduced pancreas size (, , ) and exocrine dysfunction [summarized in ()]. The reduction in pancreas weight and volume in T1D could occur even before beta-cell loss and disease diagnosis (, ), and might be one of the early T1D predisposing changes in the pancreas. Curiously, recent studies have shown preserved acinar cell size () and density () in the pancreases of donors with T1D, implying that the reduced pancreas volume in T1D could be due to a reduced number of acinar cells rather than acinar atrophy. Interestingly, in the current study, higher expression of core acinar cell genes encoding for digestive enzymes was observed in the pancreas tissue samples from the recent onset T1D patients than from the controls. Increased expression of genes encoding for digestive enzymes might actually reflect an attempt of the remaining acinar cells to compensate for the reduced functional capacity associated with a reduced pancreas volume, and could even be linked to a shift in the balance between different acinar cell subpopulations, as suggested by Segerstolpe et al. (). Although our cell type deconvolution analysis did not reveal any major differences in the proportion of acinar cells between the DiViD patient samples and the controls, further studies at single cell level would be required to confirm this as the accuracy of computational deconvolution is lower than that of the corresponding experimental approaches ().
In addition to the subpopulation of acinar cells expressing high levels of genes encoding for secretory digestive enzymes, Segerstolpe et al. suggested the presence of another acinar cell subpopulation characterized by the elevated expression levels of inflammatory and immune associated genes (). Interestingly, a clear observation in the current pancreas tissue RNA-sequencing data was the downregulation of several of these genes, including LCN2, CLDN2 and MHC class II genes, in the pancreas of the recent onset T1D patients. The potentially reduced numbers of acinar cells expressing inflammation and immune regulatory genes could early on contribute to making the exocrine pancreas less capable of responding to immune challenges, such as virus infections. Still, the regulation of immune and inflammatory processes occurs at multiple levels, and other factors are also likely to influence this. For example, the inhibitory SIGLEC receptors dampen the inflammatory and immune responses, but they are also thought to have a role in the maintenance of self-tolerance (). The inhibitory SIGLEC receptor, SIGLEC11, known to be expressed in the pancreas (), was the most strongly downregulated gene in the pancreas of the recent onset T1D patients, and in earlier studies it has been associated with neuroinflammation (, ). Interestingly, while downregulated in the whole pancreas, SIGLEC11 was not detected as significantly regulated in the pancreatic islets of the DiViD patients when compared to non-diabetic controls. Moreover, the downregulation of immune and inflammatory genes in the whole pancreas differs from earlier reports from the pancreatic islets and beta-cells, showing T1D associated upregulation of such genes, including for example MHC class I and II genes, and interferon response related STAT1 and MX1 in (, ). This suggests largely differing expression patterns for immune and inflammatory genes between the endocrine and the exocrine pancreas. However, although the global expression of immune and inflammatory genes was downregulated in the pancreas tissue samples of the individuals with recent onset T1D in the current data, local infiltration of immune cells and activation of immune responses might still be present also in these samples.
The interplay between the exocrine and the endocrine pancreas is thought to have an important role in maintaining pancreatic homeostasis (). Reg-family genes are thought to promote pancreatic islet growth and regeneration (), and REG3A was discovered as a marker of acinar cells close to the pancreatic islets (). Strong downregulation of REG3A, and also REG3G and REG1B in the pancreas of T1D patients in the current study could therefore reflect impaired interplay between the exocrine and the endocrine pancreas. Furthermore, recent reports on intermediate cell types containing simultaneously both exocrine and endocrine features () or mixed characteristics of different endocrine cell types (), further complicate the classification and functional understanding of distinct cell types. Therefore, broader studies at the level of single cells are required to understand the presence and relationships between the different cell subtypes in the pancreas in T1D. However, in addition to being related to pancreatic islet growth, regeneration and the acinar cell type, members of the third subclass of the Reg family, have been observed to be associated with host defenses and hyperglycemia and it has been suggested that bacteria, injury, interleukins () or hyperglycemic conditions () can regulate or inhibit the expression of REG3A, further complicating the picture for the possible underlying cause for the observed gene expression.
Enteroviruses are thought to be one of the potential triggers for the development of T1D. Several studies have shown the presence of enteroviruses in the pancreatic beta-cells of T1D patients more often than in control populations (, , ). In addition, recent studies have shown the detection of enterovirus RNA also in the exocrine pancreas of T1D and non-diabetic individuals (, ). Krogvold et al. explored the presence of enteroviruses in the pancreas of the DiViD patients using immunohistochemistry (IHC), reverse transcription polymerase chain reaction (RT-PCR) and RNA-sequencing (). While no viral genome was detected in the pancreas tissue samples using RNA-sequencing, sensitive RT-PCR detected enterovirus RNA in the pancreas tissue in one of them, the DiViD case 6. On the other hand, pancreatic islets were positive for the enterovirus VP1 protein in IHC in all the six DiViD cases and culture supernatants of isolated islets were positive for enterovirus RNA by RT-PCR in the DiViD cases 2, 4, 5 and 6. Finally, an earlier study showed overexpression of several interferon stimulated genes in the insulitic islets of the DiViD patients, when compared to non-insulitic islets from organ donors without diabetes (). Our results, dominated by the exocrine pancreas, showed clearly higher expression of central antiviral immune response genes in the pancreas tissue sample of DiViD case 6 compared to all the other DiViD cases, and are thereby well in line with the earlier data on the presence of enteroviruses (). Many of the antiviral immune response genes upregulated in the pancreas tissue sample of DiViD case 6 were also included in an enterovirus induced blood gene expression signature defined in our earlier study (), supporting the link to enterovirus infection. On the other hand, studying the same set of central antiviral immune response genes in microarray data collected from the pancreatic islets of the same individuals showed no clear differences between the DiViD cases, indicating similar infection status in their pancreatic islets. These results highlight the importance of studying both the exocrine and the endocrine pancreas at an individual level when aiming to understand the role of enteroviruses in the development of T1D.
In our whole pancreas transcriptome 1, clear downregulation of genes and processes related to extracellular matrix organization, cell migration, cell adhesion, cell motility and cell surface receptors was observed (Supplementary Figure 2). Similar downregulation of genes and processes associated with cell-cell contacts, extracellular matrix and cell surface receptors was observed in a recent study by us and others as a result of a carrier-state-type persistent infection by two Group B coxsackievirus (CVB) strains in a pancreatic cell line (). Moreover, many cell surface receptors are known for enteroviruses including members of the integrin and immunoglobulin superfamilies typically involved in cell adhesion (). In earlier studies over two decades ago, it was observed that the human enterovirus coxsackievirus A21 binds with high affinity to the intercellular adhesion molecule 1 (ICAM1) (, ) for cell entry, a phenomenon already known earlier in relation to several rhinoviruses. Persistent infection by coxsackie B viruses of human microvascular endothelial cells have been observed to result in increased expression of adhesion molecules, including ICAM1, especially shortly after the infection (<30-40 days after infection) (), fluctuating to lower levels afterwards. More generally, the interplay between viral proteins and host proteins, namely cellular receptors and proteins, largely determines the attachment, entry and internalization of the virus into the host cells (). Interestingly, among other cell adhesion molecules, ICAM1 was also observed to be downregulated in the recently diagnosed DiViD patients when compared to the non-diabetic controls in our whole pancreas data (Supplementary Table 2). Such downregulation of cell adhesion and cell surface receptors observed in this study and as a response to persistent carrier state CVB infection by (), can reflect the hosts attempt in controlling the infection after a prolonged infectious state. However, even though there is a large amount of evidence linking enteroviruses and viruses in general to cell surface molecules related to cell adhesion, the exact mechanisms of virus-host interplay regarding these molecules are unknown and require further research.
Due to the invasiveness of pancreas tissue sampling, there are clear limitations in accessing appropriate control samples for the unique DiViD study of living young adults with recently diagnosed T1D. In the current study, three commercial pancreas total tissue controls from different sources and with differing cadaveric donor characteristics were used with an aim to minimize systematic bias due to the limitations with the control samples. Overall, our data suggests that the development of T1D is associated with co-occurring changes in the exocrine and the endocrine pancreas, and might involve an imbalance and impaired communication between different cell populations in the pancreas. In the future, broader comparisons across different levels of data from deep bulk omics analyses to single cell profiling and imaging analyses with spatial information could provide valuable knowledge for deeper understanding on the pathogenesis of T1D. Such cross-platform studies could also better take into account the heterogeneity of pancreas tissue and its effects on the unique readouts of each platform.
Funding
The DiViD study was funded by South-Eastern Norway Regional Health Authority (Grant to KD-J), The Novo Nordisk Foundation (Grant to KD-J) and through PEVNET. The present work was financially supported by the European Commission (Persistent Virus Infection in Diabetes Network [PEVNET] Frame Programme 7, contract number 261441). RL and LE groups are also supported by InFLAMES Flagship Programme of the Academy of Finland (decision number: 337530). NL was supported by the Academy of Finland decision no. 287423. RL received funding from the Academy of Finland (grants 292335, 294337, 319280, 31444, 319280, 329277, 331790), Business Finland and by grants from the JDRF, the Sigrid Jusélius Foundation (SJF), Jane and Aatos Erkko Foundation, the Novo Nordisk Foundation, Finnish Diabetes Foundation and the Finnish Cancer Foundation. LE reports grants from the European Research Council ERC (677943), European Union’s Horizon 2020 research and innovation programme (955321), Academy of Finland (296801, 310561, 314443, 329278, 335434 and 335611), and Sigrid Juselius Foundation, during the conduct of the study. CM and IG were supported by National Institutes of Health, UC4 DK104155 and the Juvenile Diabetes Research Foundation, JDRF 47-2013-520. Our research was also supported by the University of Turku Graduate School (UTUGS).
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this article are not readily available because of the sensitive nature of the data and possible high risks associated with patient confidentiality. Requests to access the datasets should be directed to knut.dahl-jorgensen@medisin.uio.no.
Ethics statement
DiViD study was approved by the Norwegian Government´s Regional Ethics Committee, and written informed consent was obtained from all patients. The patients/participants provided their written informed consent to participate in this study.
Author contributions
TV participated in designing the study, drafted the manuscript with NL and MJ, processed all the datasets included in the study, performed the data analysis and interpreted the results. NL participated in designing the study, drafted the manuscript with TV and MJ, analyzed the data and interpreted the results. MJ participated in drafting the manuscript with TV and NL and performed the cell type deconvolution analysis. LK was responsible for clinical coordination, the recruitment of patients, data collection, the whole transcriptome sequencing of the whole tissue transcriptome 2 together with ME and KD-J and participated in interpretation of the results. ME was responsible for the whole transcriptome sequencing of the whole tissue transcriptome 2 together with LK and KD-J. HK and ST participated in designing the study, performed the whole transcriptome sequencing analysis of the whole tissue transcriptome 1 and participated in interpreting the results. CM was responsible for the analysis of the islet transcriptome data together with IG. IG was responsible for generating the islet transcriptome data together with CM and participated in interpreting the results. SO was responsible for generating the pancreatic tissue samples for the whole tissue transcriptome 1 together with HH. HH was responsible for generating the pancreatic tissue samples for the whole tissue transcriptome 1 together with SO and participated in interpreting the results. KD-J had the initial idea and design of the DiViD study, participated in study design, funding, regulatory issues, international collaboration, data collection and interpreting the results. LE was responsible for the study design, supervised the performed data analysis, and participated in interpreting the results and drafting the manuscript. RL was responsible for the study design, supervised the study, and participated in interpreting the results and drafting the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors thank specialist nurse Trine Roald (Pediatric Department, Oslo University Hospital) for invaluable efforts in coordination of the study, the nurses and physicians at the local hospitals for providing contact with the patients, and the patients who participated in this study. A special thanks to Professor Bjørn Edwin and Professor Trond Buanes at the Intervention Centre and Department of Gastrosurgery, Oslo University Hospital for performing the minimal pancreatic tail resection providing the tissues of the DiViD cases. We would like to thank Sarita Heinonen and Marjo Hakkarainen for their excellent technical assistance. We acknowledge the Turku Bioscience Centre’s core facility, the Finnish Functional Genomics Centre (FFGC) supported by Biocenter Finland, for their assistance. We acknowledge the Finnish Centre for Scientific Computing (CSC) and Elixir Finland for computational resources.
Conflict of interest
Author HH is employed by Fimlab Laboratories Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.861985/full#supplementary-material
References
1
EisenbarthGS. Type I Diabetes Mellitus. A Chronic Autoimmune Dis N Engl J Med (1986) 314(21):1360–8. doi: 10.1056/NEJM198605223142106
2
WarshauerJTBluestoneJAAndersonMS. New Frontiers in the Treatment of Type 1 Diabetes. Cell Metab (2020) 31(1):46–61. doi: 10.1016/j.cmet.2019.11.017
3
InselRADunneJLAtkinsonMAChiangJLDabeleaDGottliebPAet al. Staging Presymptomatic Type 1 Diabetes: A Scientific Statement of JDRF, the Endocrine Society, and the American Diabetes Association. Diabetes Care (2015) 38(10):1964–74. doi: 10.2337/dc15-1419
4
RedondoMJJeffreyJFainPREisenbarthGSOrbanT. Concordance for Islet Autoimmunity Among Monozygotic Twins. N Engl J Med (2008) 359(26):2849–50. doi: 10.1056/NEJMc0805398
5
BonifacioEZieglerAG. Advances in the Prediction and Natural History of Type 1 Diabetes. Endocrinol Metab Clin North Am (2010) 39(3):513–25. doi: 10.1016/j.ecl.2010.05.007
6
KrogvoldLEdwinBBuanesTLudvigssonJKorsgrenOHyötyHet al. Pancreatic Biopsy by Minimal Tail Resection in Live Adult Patients at the Onset of Type 1 Diabetes: Experiences From the DiViD Study. Diabetologia (2014) 57(4):841–3. doi: 10.1007/s00125-013-3155-y
7
KrogvoldLSkogOSundstromGEdwinBBuanesTHanssenKFet al. Function of Isolated Pancreatic Islets From Patients at Onset of Type 1 Diabetes: Insulin Secretion Can Be Restored After Some Days in a Nondiabetogenic Environment In Vitro: Results From the DiViD Study. Diabetes (2015) 64(7):2506–12. doi: 10.2337/db14-1911
8
KrogvoldLEdwinBBuanesTFriskGSkogOAnagandulaMet al. Detection of a Low-Grade Enteroviral Infection in the Islets of Langerhans of Living Patients Newly Diagnosed With Type 1 Diabetes. Diabetes (2015) 64(5):1682–7. doi: 10.2337/db14-1370
9
Alexandre-HeymannLMalloneRBoitardCScharfmannRLargerE. Structure and Function of the Exocrine Pancreas in Patients With Type 1 Diabetes. Rev Endocr Metab Disord (2019) 20(2):129–49. doi: 10.1007/s11154-019-09501-3
10
Rodriguez-CalvoTEkwallOAmirianNZapardiel-GonzaloJvon HerrathMG. Increased Immune Cell Infiltration of the Exocrine Pancreas: A Possible Contribution to the Pathogenesis of Type 1 Diabetes. Diabetes (2014) 63(11):3880–90. doi: 10.2337/db14-0549
11
KondrashovaANurminenNLehtonenJHyötyMToppariJllonenJet al. Exocrine Pancreas Function Decreases During the Progression of the Beta-Cell Damaging Process in Young Prediabetic Children. Pediatr Diabetes (2018) 19(3):398–402. doi: 10.1111/pedi.12592
12
FoulisAKStewartJA. The Pancreas in Recent-Onset Type 1 (Insulin-Dependent) Diabetes Mellitus: Insulin Content of Islets, Insulitis and Associated Changes in the Exocrine Acinar Tissue. Diabetologia (1984) 26(6):456–61. doi: 10.1007/bf00262221
13
ValleAGiamporcaroGMScaviniMStabiliniAGroganPEleonoraBet al. Reduction of Circulating Neutrophils Precedes and Accompanies Type 1 Diabetes. Diabetes (2013) 62(6):2072–7. doi: 10.2337/db12-1345
14
Rodriguez-CalvoT. Enterovirus Infection and Type 1 Diabetes: Unraveling the Crime Scene. Clin Exp Immunol (2019) 195(1):15–24. doi: 10.1111/cei.13223
15
FilippiCMvon HerrathMG. Viral Trigger for Type 1 Diabetes: Pros and Cons. Diabetes (2008) 57(11):2863–71. doi: 10.2337/db07-1023
16
JunH-SYoonJ-W. A New Look at Viruses in Type 1 Diabetes. Diabetes Metab Res Rev (2003) 19(1):8–31. doi: 10.1002/dmrr.337
17
LaitinenOHHonkanenHPakkanenOOikarinenSHankaniemiMMHutalaHet al. Coxsackievirus B1 Is Associated With Induction of β-Cell Autoimmunity That Portends Type 1 Diabetes. Diabetes (2014) 63(2):446–55. doi: 10.2337/db13-0619
18
OikarinenSTauriainenSHoberDLucasBVazeouASioofy-KhojineAet al. Virus Antibody Survey in Different European Populations Indicates Risk Association Between Coxsackievirus B1 and Type 1 Diabetes. Diabetes (2014) 63(2):655–62. doi: 10.2337/db13-0620
19
YlipaastoPKlingelKLindbergAMOtonkoskiTKandolfRHoviTet al. Enterovirus Infection in Human Pancreatic Islet Cells, Islet Tropism In Vivo and Receptor Involvement in Cultured Islet Beta Cells. Diabetologia (2004) 47(2):225–39. doi: 10.1007/s00125-003-1297-z
20
TauriainenSOikarinenSOikarinenMHyötyH. Enteroviruses in the Pathogenesis of Type 1 Diabetes. Semin Immunopathol (2011) 33(1):45–55. doi: 10.1007/s00281-010-0207-y
21
BusseNParoniFRichardsonSJLaihoJEOikarinenMFriskGet al. Detection and Localization of Viral Infection in the Pancreas of Patients With Type 1 Diabetes Using Short Fluorescently-Labelled Oligonucleotide Probes. Oncotarget (2017) 8(8):12620–36. doi: 10.18632/oncotarget.14896
22
GeravandiSRichardsonSPuglieseAMaedlerK. Localization of Enteroviral RNA Within the Pancreas in Donors With T1D and T1D-Associated Autoantibodies. Cell Rep Med (2021) 2(8):100371. doi: 10.1016/j.xcrm.2021.100371
23
LiQCohenJI. Epstein-Barr Virus and the Human Leukocyte Antigen Complex. Curr Clin Microbiol Rep (2019) 6(3):175–81. doi: 10.1007/s40588-019-00120-9
24
LietzénNHirvonenKHonkimaaABuchacherTLaihoJEOikarinenSet al. Coxsackievirus B Persistence Modifies the Proteome and the Secretome of Pancreatic Ductal Cells. iScience (2019) 19:340–57. doi: 10.1016/j.isci.2019.07.040
25
HämäläinenSNurminenNAhlforsHOikarinenSSioofy-KhojineABFriskGet al. Coxsackievirus B1 Reveals Strain Specific Differences in Plasmacytoid Dendritic Cell Mediated Immunogenicity. J Med Virol (2014) 86(8):1412–20. doi: 10.1002/jmv.23903
26
BuchacherTHonkimaaAVälikangasTLietzénNHirvonenMKLaihoJEet al. Persistent Coxsackievirus B1 Infection Triggers Extensive Changes in the Transcriptome of Human Pancreatic Ductal Cells. iScience (2022) 25(1):103653. doi: 10.1016/j.isci.2021.103653
27
RichardsonSJRodriguez-CalvoTGerlingICMathewsCEKaddisJSRussellMAet al. Islet Cell Hyperexpression of HLA Class I Antigens: A Defining Feature in Type 1 Diabetes. Diabetologia (2016) 59(11):2448–58. doi: 10.1007/s00125-016-4067-4
28
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust Enumeration of Cell Subsets From Tissue Expression Profiles. Nat Methods (2015) 12(5):453–7. doi: 10.1038/nmeth.3337
29
SegerstolpeÅPalasantzaAEliassonPAnderssonEMAndréassonACSunXet al. Single-Cell Transcriptome Profiling of Human Pancreatic Islets in Health and Type 2 Diabetes. Cell Metab (2016) 24(4):593–607. doi: 10.1016/j.cmet.2016.08.020
30
EloLFilénSLahesmaaRAittokallioT. Reproducibility-Optimized Test Statistic for Ranking Genes in Microarray Studies. IEEE/ACM Trans Comput Biol Bioinform (2008) 5(3):423–31. doi: 10.1109/tcbb.2007.1078
31
SeyednasrollahFRantanenKJaakkolaPEloLL. ROTS: Reproducible RNA-Seq Biomarker Detector—Prognostic Markers for Clear Cell Renal Cell Cancer. Nucleic Acids Res (2015) 44(1):e1–1. doi: 10.1093/nar/gkv806
32
HuangDWShermanBTLempickiRA. Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat Protoc (2009) 4(1):44–57. doi: 10.1038/nprot.2008.211
33
SzklarczykDGableALLyonDJungeAWyderSHuerta-CepasJet al. STRING V11: Protein-Protein Association Networks With Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res (2019) 47(D1):D607–13. doi: 10.1093/nar/gky1131
34
MuraroMJDharmadhikariGGrünDGroenNDielenTJansenEet al. A Single-Cell Transcriptome Atlas of the Human Pancreas. Cell Syst (2016) 3(4):385–94.e3. doi: 10.1016/j.cels.2016.09.002
35
RickelsMRNorrisAWHullRL. A Tale of Two Pancreases: Exocrine Pathology and Endocrine Dysfunction. Diabetologia (2020) 63(10):2030–9. doi: 10.1007/s00125-020-05210-8
36
KramerAGreenJPollardJJTugendreichS. Causal Analysis Approaches in Ingenuity Pathway Analysis. Bioinformatics (2014) 30(4):523–30. doi: 10.1093/bioinformatics/btt703
37
ShannonPMarkielAOzierOBaligaNSWangJTRamageDet al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res (2003) 13(11):2498–504. doi: 10.1101/gr.1239303
38
LundbergMKrogvoldLKuricEDahl-JorgensenKSkogO. Expression of Interferon-Stimulated Genes in Insulitic Pancreatic Islets of Patients Recently Diagnosed With Type 1 Diabetes. Diabetes (2016) 65(10):3104–10. doi: 10.2337/db16-0616
39
DonchevaNTMorrisJHGorodkinJJensenLJ. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J Proteome Res (2019) 18(2):623–32. doi: 10.1021/acs.jproteome.8b00702
40
Campbell-ThompsonMWasserfallCMontgomeryELAtkinsonMAKaddisJS. Pancreas Organ Weight in Individuals With Disease-Associated Autoantibodies at Risk for Type 1 Diabetes. JAMA (2012) 308(22):2337–9. doi: 10.1001/jama.2012.15008
41
Campbell-ThompsonMLFilippSLGrajoJRNambanBBeeglesRMiddlebrooksEHet al. Relative Pancreas Volume Is Reduced in First-Degree Relatives of Patients With Type 1 Diabetes. Diabetes Care (2019) 42(2):281–7. doi: 10.2337/dc18-1512
42
WrightJJSaundersDCDaiCPoffenbergerGCairnsBSerrezeDVet al. Decreased Pancreatic Acinar Cell Number in Type 1 Diabetes. Diabetologia (2020) 63(7):1418–23. doi: 10.1007/s00125-020-05155-y
43
GranlundLHedinAWahlhütterMSeironPKorsgrenOSkogOet al. Histological and Transcriptional Characterization of the Pancreatic Acinar Tissue in Type 1 Diabetes. BMJ Open Diabetes Res Care (2021) 9(1). doi: 10.1136/bmjdrc-2020-002076
44
Avila CobosFVandesompeleJMestdaghPDe PreterK. Computational Deconvolution of Transcriptomics Data From Mixed Cell Populations. Bioinformatics (2018) 34(11):1969–79. doi: 10.1093/bioinformatics/bty019
45
MacauleyMSCrockerPRPaulsonJC. Siglec-Mediated Regulation of Immune Cell Function in Disease. Nat Rev Immunol (2014) 14(10):653–66. doi: 10.1038/nri3737
46
AngataTKerrSCGreavesDRVarkiNMCrockerPRVarkiA. Cloning and Characterization of Human Siglec-11. A Recently Evolved Signaling Molecule That can Interact With SHP-1 and SHP-2 and is Expressed by Tissue Macrophages, Including Brain Microglia. J Biol Chem (2002) 277(27):24466–74. doi: 10.1074/jbc.M202833200
47
RussellMARedickSDBlodgettDMRichardsonSJLeetePKrogvoldLet al. HLA Class II Antigen Processing and Presentation Pathway Components Demonstrated by Transcriptome and Protein Analyses of Islet β-Cells From Donors With Type 1 Diabetes. Diabetes (2019) 68(5):988–1001. doi: 10.2337/db18-0686
48
OkamotoH. The Reg Gene Family and Reg Proteins: With Special Attention to the Regeneration of Pancreatic Beta-Cells. J Hepatobil Pancreat Surg (1999) 6(3):254–62. doi: 10.1007/s005340050115
49
de BoerPPirozziNMWoltersAHGKuipersJKusmartsevaIAtkinsonMAet al. Large-Scale Electron Microscopy Database for Human Type 1 Diabetes. Nat Commun (2020) 11(1):2475. doi: 10.1038/s41467-020-16287-5
50
SeironPWibergAKuricEKrogvoldLJahnsenFLDahl-JørgensenKet al. Characterisation of the Endocrine Pancreas in Type 1 Diabetes: Islet Size is Maintained But Islet Number is Markedly Reduced. J Pathol Clin Res (2019) 5(4):248–55. doi: 10.1002/cjp2.140
51
WangLQuanYZhuYXieXWangZWangLet al. The Regenerating Protein 3A: A Crucial Molecular With Dual Roles in Cancer. Mol Biol Rep (2022) 49(2):1491–500. doi: 10.1007/s11033-021-06904-x
52
WuYQuanYLiuYLiuKLiHJiangZet al. Hyperglycaemia Inhibits REG3A Expression to Exacerbate TLR3-Mediated Skin Inflammation in Diabetes. Nat Commun (2016) 7(1):13393. doi: 10.1038/ncomms13393
53
RichardsonSJWillcoxABoneAJFoulisAKMorganNG. The Prevalence of Enteroviral Capsid Protein Vp1 Immunostaining in Pancreatic Islets In Human Type 1 Diabetes. Diabetologia (2009) 52(6):1143–51. doi: 10.1007/s00125-009-1276-0
54
RichardsonSJLeetePBoneAJFoulisAKMorganNG. Expression of the Enteroviral Capsid Protein VP1 in the Islet Cells of Patients With Type 1 Diabetes is Associated With Induction of Protein Kinase R and Downregulation of Mcl-1. Diabetologia (2013) 56(1):185–93. doi: 10.1007/s00125-012-2745-4
55
LietzenNAnLTTJaakkolaMKKallionpääHOikarinenSMykkänenJet al. Enterovirus-Associated Changes in Blood Transcriptomic Profiles of Children With Genetic Susceptibility to Type 1 Diabetes. Diabetologia (2018) 61(2):381–8. doi: 10.1007/s00125-017-4460-7
56
ShafrenDRDorahyDJInghamRABurnsGFBarryRD. Coxsackievirus A21 Binds to Decay-Accelerating Factor But Requires Intercellular Adhesion Molecule 1 for Cell Entry. J Virol (1997) 71(6):4736–43. doi: 10.1128/JVI.71.6.4736-4743.1997
57
ShafrenDRDorahyDJGreiveSJBurnsGFBarryRD. Mouse Cells Expressing Human Intercellular Adhesion Molecule-1 Are Susceptible to Infection by Coxsackievirus A21. J Virol (1997) 71(1):785–9. doi: 10.1128/JVI.71.1.785-789.1997
58
ZanoneMMFavaroEConaldiPGGreeningJBotteliAPerinPCet al. Persistent Infection of Human Microvascular Endothelial Cells by Coxsackie B Viruses Induces Increased Expression of Adhesion Molecules. J Immunol (2003) 171(1):438–46. doi: 10.4049/jimmunol.171.1.438
59
OwinoCOChuJJH. Recent Advances on the Role of Host Factors During non-Poliovirus Enteroviral Infections. J BioMed Sci (2019) 26(1):47. doi: 10.1186/s12929-019-0540-y
Summary
Keywords
DiViD, exocrine pancreas, gene expression, pancreatic islet, pancreas, transcriptomics, type 1 diabetes
Citation
Välikangas T, Lietzén N, Jaakkola MK, Krogvold L, Eike MC, Kallionpää H, Tuomela S, Mathews C, Gerling IC, Oikarinen S, Hyöty H, Dahl-Jorgensen K, Elo LL and Lahesmaa R (2022) Pancreas Whole Tissue Transcriptomics Highlights the Role of the Exocrine Pancreas in Patients With Recently Diagnosed Type 1 Diabetes. Front. Endocrinol. 13:861985. doi: 10.3389/fendo.2022.861985
Received
25 January 2022
Accepted
09 March 2022
Published
13 April 2022
Volume
13 - 2022
Edited by
Sangeeta Dhawan, City of Hope National Medical Center, United States
Reviewed by
Joachim Størling, Steno Diabetes Center Copenhagen (SDCC), Denmark; Jon D. Piganelli, University of Pittsburgh, United States
Updates
Copyright
© 2022 Välikangas, Lietzén, Jaakkola, Krogvold, Eike, Kallionpää, Tuomela, Mathews, Gerling, Oikarinen, Hyöty, Dahl-Jorgensen, Elo and Lahesmaa.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Riitta Lahesmaa, rilahes@utu.fi; Laura L. Elo, laura.elo@utu.fi
†These authors have contributed equally to this work
This article was submitted to Diabetes: Molecular Mechanisms, a section of the journal Frontiers in Endocrinology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.