Distinct cellular immune responses in children en route to type 1 diabetes with different first-appearing autoantibodies.

Previous studies have revealed heterogeneity in the progression to clinical type 1 diabetes in children who develop islet-specific antibodies either to insulin (IAA) or glutamic acid decarboxylase (GADA) as the first autoantibodies. Here, we test the hypothesis that children who later develop clinical disease have different early immune responses, depending on the type of the first autoantibody to appear (GADA-first or IAA-first). We use mass cytometry for deep immune profiling of peripheral blood mononuclear cell samples longitudinally collected from children who later progressed to clinical disease (IAA-first, GADA-first, ≥2 autoantibodies first groups) and matched for age, sex, and HLA controls who did not, as part of the Type 1 Diabetes Prediction and Prevention study. We identify differences in immune cell composition of children who later develop disease depending on the type of autoantibodies that appear first. Notably, we observe an increase in CD161 expression in natural killer cells of children with ≥2 autoantibodies and validate this in an independent cohort. The results highlight the importance of endotype-specific analyses and are likely to contribute to our understanding of pathogenic mechanisms underlying type 1 diabetes development.

Long noncoding RNA LIRIL2R modulates FOXP3 levels and suppressive function of human CD4+ regulatory T cells by regulating IL2RA.

Regulatory T cells (Tregs) are central in controlling immune responses, and dysregulation of their function can lead to autoimmune disorders or cancer. Despite extensive studies on Tregs, the basis of epigenetic regulation of human Treg development and function is incompletely understood. Long intergenic noncoding RNAs (lincRNA)s are important for shaping and maintaining the epigenetic landscape in different cell types. In this study, we identified a gene on the chromosome 6p25.3 locus, encoding a lincRNA, that was up-regulated during early differentiation of human Tregs. The lincRNA regulated the expression of interleukin-2 receptor alpha (IL2RA), and we named it the lincRNA regulator of IL2RA (LIRIL2R). Through transcriptomics, epigenomics, and proteomics analysis of LIRIL2R-deficient Tregs, coupled with global profiling of LIRIL2R binding sites using chromatin isolation by RNA purification, followed by sequencing, we identified IL2RA as a target of LIRIL2R. This nuclear lincRNA binds upstream of the IL2RA locus and regulates its epigenetic landscape and transcription. CRISPR-mediated deletion of the LIRIL2R-bound region at the IL2RA locus resulted in reduced IL2RA expression. Notably, LIRIL2R deficiency led to reduced expression of Treg-signature genes (e.g., FOXP3, CTLA4, and PDCD1), upregulation of genes associated with effector T cells (e.g., SATB1 and GATA3), and loss of Treg-mediated suppression.

A proximal enhancer regulates RORA expression during early human Th17 cell differentiation.

Gene regulatory elements, such as enhancers, greatly influence cell identity by tuning the transcriptional activity of specific cell types. Dynamics of enhancer landscape during early human Th17 cell differentiation remains incompletely understood. Leveraging ATAC-seq-based profiling of chromatin accessibility and comprehensive analysis of key histone marks, we identified a repertoire of enhancers that potentially exert control over the fate specification of Th17 cells. We found 23 SNPs associated with autoimmune diseases within Th17-enhancers that precisely overlapped with the binding sites of transcription factors actively engaged in T-cell functions. Among the Th17-specific enhancers, we identified an enhancer in the intron of RORA and demonstrated that this enhancer positively regulates RORA transcription. Moreover, CRISPR-Cas9-mediated deletion of a transcription factor binding site-rich region within the identified RORA enhancer confirmed its role in regulating RORA transcription. These findings provide insights into the potential mechanism by which the RORA enhancer orchestrates Th17 differentiation.

PIM kinases regulate early human Th17 cell differentiation.

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

HIC1 interacts with FOXP3 multi protein complex: Novel pleiotropic mechanisms to regulate human regulatory T cell differentiation and function.

Transcriptional repressor, hypermethylated in cancer 1 (HIC1) participates in a range of important biological processes, such as tumor repression, immune suppression, embryonic development and epigenetic gene regulation. Further to these, we previously demonstrated that HIC1 provides a significant contribution to the function and development of regulatory T (Treg) cells. However, the mechanism by which it regulates these processes was not apparent. To address this question, we used affinity-purification mass spectrometry to characterize the HIC1 interactome in human Treg cells. Altogether 61 high-confidence interactors were identified, including IKZF3, which is a key transcription factor in the development of Treg cells. The biological processes associated with these interacting proteins include protein transport, mRNA processing, non-coding (ncRNA) transcription and RNA metabolism. The results revealed that HIC1 is part of a FOXP3-RUNX1-CBFB protein complex that regulates Treg signature genes thus improving our understanding of HIC1 function during early Treg cell differentiation.

Serum APOC1 levels are decreased in young autoantibody positive children who rapidly progress to type 1 diabetes.

Better understanding of the early events in the development of type 1 diabetes is needed to improve prediction and monitoring of the disease progression during the substantially heterogeneous presymptomatic period of the beta cell damaging process. To address this concern, we used mass spectrometry-based proteomics to analyse longitudinal pre-onset plasma sample series from children positive for multiple islet autoantibodies who had rapidly progressed to type 1 diabetes before 4 years of age (n = 10) and compared these with similar measurements from matched children who were either positive for a single autoantibody (n = 10) or autoantibody negative (n = 10). Following statistical analysis of the longitudinal data, targeted serum proteomics was used to verify 11 proteins putatively associated with the disease development in a similar yet independent and larger cohort of children who progressed to the disease within 5 years of age (n = 31) and matched autoantibody negative children (n = 31). These data reiterated extensive age-related trends for protein levels in young children. Further, these analyses demonstrated that the serum levels of two peptides unique for apolipoprotein C1 (APOC1) were decreased after the appearance of the first islet autoantibody and remained relatively less abundant in children who progressed to type 1 diabetes, in comparison to autoantibody negative children.

Circulating metabolic signatures of rapid and slow progression to type 1 diabetes in islet autoantibody-positive children.

Aims/hypothesis: Appearance of multiple islet cell autoantibodies in early life is indicative of future progression to overt type 1 diabetes, however, at varying rates. Here, we aimed to study whether distinct metabolic patterns could be identified in rapid progressors (RP, disease manifestation within 18 months after the initial seroconversion to autoantibody positivity) vs. slow progressors (SP, disease manifestation at 60 months or later from the appearance of the first autoantibody). Methods: Longitudinal samples were collected from RP (n=25) and SP (n=41) groups at the ages of 3, 6, 12, 18, 24, or ≥ 36 months. We performed a comprehensive metabolomics study, analyzing both polar metabolites and lipids. The sample series included a total of 239 samples for lipidomics and 213 for polar metabolites. Results: We observed that metabolites mediated by gut microbiome, such as those involved in tryptophan metabolism, were the main discriminators between RP and SP. The study identified specific circulating molecules and pathways, including amino acid (threonine), sugar derivatives (hexose), and quinic acid that may define rapid vs. slow progression to type 1 diabetes. However, the circulating lipidome did not appear to play a major role in differentiating between RP and SP. Conclusion/interpretation: Our study suggests that a distinct metabolic profile is linked with the type 1 diabetes progression. The identification of specific metabolites and pathways that differentiate RP from SP may have implications for early intervention strategies to delay the development of type 1 diabetes.

Targeted serum proteomics of longitudinal samples from newly diagnosed youth with type 1 diabetes distinguishes markers of disease and C-peptide trajectory.

AIMS/HYPOTHESIS: There is a growing need for markers that could help indicate the decline in beta cell function and recognise the need and efficacy of intervention in type 1 diabetes. Measurements of suitably selected serum markers could potentially provide a non-invasive and easily applicable solution to this challenge. Accordingly, we evaluated a broad panel of proteins previously associated with type 1 diabetes in serum from newly diagnosed individuals during the first year from diagnosis. To uncover associations with beta cell function, comparisons were made between these targeted proteomics measurements and changes in fasting C-peptide levels. To further distinguish proteins linked with the disease status, comparisons were made with measurements of the protein targets in age- and sex-matched autoantibody-negative unaffected family members (UFMs). METHODS: Selected reaction monitoring (SRM) mass spectrometry analyses of serum, targeting 85 type 1 diabetes-associated proteins, were made. Sera from individuals diagnosed under 18 years (n=86) were drawn within 6 weeks of diagnosis and at 3, 6 and 12 months afterwards (288 samples in total). The SRM data were compared with fasting C-peptide/glucose data, which was interpreted as a measure of beta cell function. The protein data were further compared with cross-sectional SRM measurements from UFMs (n=194). RESULTS: Eleven proteins had statistically significant associations with fasting C-peptide/glucose. Of these, apolipoprotein L1 and glutathione peroxidase 3 (GPX3) displayed the strongest positive and inverse associations, respectively. Changes in GPX3 levels during the first year after diagnosis indicated future fasting C-peptide/glucose levels. In addition, differences in the levels of 13 proteins were observed between the individuals with type 1 diabetes and the matched UFMs. These included GPX3, transthyretin, prothrombin, apolipoprotein C1 and members of the IGF family. CONCLUSIONS/INTERPRETATION: The association of several targeted proteins with fasting C-peptide/glucose levels in the first year after diagnosis suggests their connection with the underlying changes accompanying alterations in beta cell function in type 1 diabetes. Moreover, the direction of change in GPX3 during the first year was indicative of subsequent fasting C-peptide/glucose levels, and supports further investigation of this and other serum protein measurements in future studies of beta cell function in type 1 diabetes.

Exhausted T cells hijacking the cancer-immunity cycle: Assets and liabilities.

T cell exhaustion is an alternative differentiation path of T cells, sometimes described as a dysfunction. During the last decade, insights of T cell exhaustion acting as a bottle neck in the field of cancer immunotherapy have undoubtedly provoked attention. One of the main drivers of T cell exhaustion is prolonged antigen presentation, a prerequisite in the cancer-immunity cycle. The umbrella term "T cell exhaustion" comprises various stages of T cell functionalities, describing the dynamic, one-way exhaustion process. Together these qualities of T cells at the exhaustion continuum can enable tumor clearance, but if the exhaustion acquired timeframe is exceeded, tumor cells have increased possibilities of escaping immune system surveillance. This could be considered a tipping point where exhausted T cells switch from an asset to a liability. In this review, the contrary role of exhausted T cells is discussed.

Gene expression signature predicts rate of type 1 diabetes progression.

BACKGROUND: Type 1 diabetes is a complex heterogenous autoimmune disease without therapeutic interventions available to prevent or reverse the disease. This study aimed to identify transcriptional changes associated with the disease progression in patients with recent-onset type 1 diabetes. METHODS: Whole-blood samples were collected as part of the INNODIA study at baseline and 12 months after diagnosis of type 1 diabetes. We used linear mixed-effects modelling on RNA-seq data to identify genes associated with age, sex, or disease progression. Cell-type proportions were estimated from the RNA-seq data using computational deconvolution. Associations to clinical variables were estimated using Pearson's or point-biserial correlation for continuous and dichotomous variables, respectively, using only complete pairs of observations. FINDINGS: We found that genes and pathways related to innate immunity were downregulated during the first year after diagnosis. Significant associations of the gene expression changes were found with ZnT8A autoantibody positivity. Rate of change in the expression of 16 genes between baseline and 12 months was found to predict the decline in C-peptide at 24 months. Interestingly and consistent with earlier reports, increased B cell levels and decreased neutrophil levels were associated with the rapid progression. INTERPRETATION: There is considerable individual variation in the rate of progression from appearance of type 1 diabetes-specific autoantibodies to clinical disease. Patient stratification and prediction of disease progression can help in developing more personalised therapeutic strategies for different disease endotypes. FUNDING: A full list of funding bodies can be found under Acknowledgments.

Profiling of peripheral blood B-cell transcriptome in children who developed coeliac disease in a prospective study.

Background: In coeliac disease (CoD), the role of B-cells has mainly been considered to be production of antibodies. The functional role of B-cells has not been analysed extensively in CoD. Methods: We conducted a study to characterize gene expression in B-cells from children developing CoD early in life using samples collected before and at the diagnosis of the disease. Blood samples were collected from children at risk at 12, 18, 24 and 36 months of age. RNA from peripheral blood CD19+ cells was sequenced and differential gene expression was analysed using R package Limma. Findings: Overall, we found one gene, HNRNPL, modestly downregulated in all patients (logFC -0·7; q = 0·09), and several others downregulated in those diagnosed with CoD already by the age of 2 years. Interpretation: The data highlight the role of B-cells in CoD development. The role of HNRPL in suppressing enteroviral replication suggests that the predisposing factor for both CoD and enteroviral infections is the low level of HNRNPL expression. Funding: EU FP7 grant no. 202063, EU Regional Developmental Fund and research grant PRG712, The Academy of Finland Centre of Excellence in Molecular Systems Immunology and Physiology Research (SyMMyS) 2012-2017, grant no. 250114) and, AoF Personalized Medicine Program (grant no. 292482), AoF grants 292335, 294337, 319280, 31444, 319280, 329277, 331790) and grants from the Sigrid Jusélius Foundation (SJF).

Long Intergenic Noncoding RNA MIAT as a Regulator of Human Th17 Cell Differentiation.

T helper 17 (Th17) cells protect against fungal and bacterial infections and are implicated in autoimmunity. Several long intergenic noncoding RNAs (lincRNA) are induced during Th17 differentiation, however, their contribution to Th17 differentiation is poorly understood. We aimed to characterize the function of the lincRNA Myocardial Infarction Associated Transcript (MIAT) during early human Th17 cell differentiation. We found MIAT to be upregulated early after induction of human Th17 cell differentiation along with an increase in the chromatin accessibility at the gene locus. STAT3, a key regulator of Th17 differentiation, directly bound to the MIAT promoter and induced its expression during the early stages of Th17 cell differentiation. MIAT resides in the nucleus and regulates the expression of several key Th17 genes, including IL17A, IL17F, CCR6 and CXCL13, possibly by altering the chromatin accessibility of key loci, including IL17A locus. Further, MIAT regulates the expression of protein kinase C alpha (PKCα), an upstream regulator of IL17A. A reanalysis of published single-cell RNA-seq data showed that MIAT was expressed in T cells from the synovium of RA patients. Our results demonstrate that MIAT contributes to human Th17 differentiation by upregulating several genes implicated in Th17 differentiation. High MIAT expression in T cells of RA patient synovia suggests a possible role of MIAT in Th17 mediated autoimmune pathologies.

Biomarkers of nanomaterials hazard from multi-layer data.

There is an urgent need to apply effective, data-driven approaches to reliably predict engineered nanomaterial (ENM) toxicity. Here we introduce a predictive computational framework based on the molecular and phenotypic effects of a large panel of ENMs across multiple in vitro and in vivo models. Our methodology allows for the grouping of ENMs based on multi-omics approaches combined with robust toxicity tests. Importantly, we identify mRNA-based toxicity markers and extensively replicate them in multiple independent datasets. We find that models based on combinations of omics-derived features and material intrinsic properties display significantly improved predictive accuracy as compared to physicochemical properties alone.

Umbilical cord blood DNA methylation in children who later develop type 1 diabetes.

AIMS/HYPOTHESIS: Distinct DNA methylation patterns have recently been observed to precede type 1 diabetes in whole blood collected from young children. Our aim was to determine whether perinatal DNA methylation is associated with later progression to type 1 diabetes. METHODS: Reduced representation bisulphite sequencing (RRBS) analysis was performed on umbilical cord blood samples collected within the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) Study. Children later diagnosed with type 1 diabetes and/or who tested positive for multiple islet autoantibodies (n = 43) were compared with control individuals (n = 79) who remained autoantibody-negative throughout the DIPP follow-up until 15 years of age. Potential confounding factors related to the pregnancy and the mother were included in the analysis. RESULTS: No differences in the umbilical cord blood methylation patterns were observed between the cases and controls at a false discovery rate <0.05. CONCLUSIONS/INTERPRETATION: Based on our results, differences between children who progress to type 1 diabetes and those who remain healthy throughout childhood are not yet present in the perinatal DNA methylome. However, we cannot exclude the possibility that such differences would be found in a larger dataset.

A systematic comparison of FOSL1, FOSL2 and BATF-mediated transcriptional regulation during early human Th17 differentiation.

Th17 cells are essential for protection against extracellular pathogens, but their aberrant activity can cause autoimmunity. Molecular mechanisms that dictate Th17 cell-differentiation have been extensively studied using mouse models. However, species-specific differences underscore the need to validate these findings in human. Here, we characterized the human-specific roles of three AP-1 transcription factors, FOSL1, FOSL2 and BATF, during early stages of Th17 differentiation. Our results demonstrate that FOSL1 and FOSL2 co-repress Th17 fate-specification, whereas BATF promotes the Th17 lineage. Strikingly, FOSL1 was found to play different roles in human and mouse. Genome-wide binding analysis indicated that FOSL1, FOSL2 and BATF share occupancy over regulatory regions of genes involved in Th17 lineage commitment. These AP-1 factors also share their protein interacting partners, which suggests mechanisms for their functional interplay. Our study further reveals that the genomic binding sites of FOSL1, FOSL2 and BATF harbour hundreds of autoimmune disease-linked SNPs. We show that many of these SNPs alter the ability of these transcription factors to bind DNA. Our findings thus provide critical insights into AP-1-mediated regulation of human Th17-fate and associated pathologies.

Pancreas Whole Tissue Transcriptomics Highlights the Role of the Exocrine Pancreas in Patients With Recently Diagnosed Type 1 Diabetes.

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.

Permutation-based significance analysis reduces the type 1 error rate in bisulphite sequencing data analysis of human umbilical cord blood samples.

DNA methylation patterns are largely established in-utero and might mediate the impacts of in-utero conditions on later health outcomes. Associations between perinatal DNA methylation marks and pregnancy-related variables, such as maternal age and gestational weight gain, have been earlier studied with methylation microarrays, which typically cover less than 2% of human CpG sites. To detect such associations outside these regions, we chose the bisulphite sequencing approach. We collected and curated clinical data on 200 newborn infants; whose umbilical cord blood samples were analysed with the reduced representation bisulphite sequencing (RRBS) method. A generalized linear mixed-effects model was fit for each high coverage CpG site, followed by spatial and multiple testing adjustment of P values to identify differentially methylated cytosines (DMCs) and regions (DMRs) associated with clinical variables, such as maternal age, mode of delivery, and birth weight. Type 1 error rate was then evaluated with a permutation analysis. We discovered a strong inflation of spatially adjusted P values through the permutation analysis, which we then applied for empirical type 1 error control. The inflation of P values was caused by a common method for spatial adjustment and DMR detection, implemented in tools comb-p and RADMeth. Based on empirically estimated significance thresholds, very little differential methylation was associated with any of the studied clinical variables, other than sex. With this analysis workflow, the sex-associated differentially methylated regions were highly reproducible across studies, technologies, and statistical models.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age.

AIMS/HYPOTHESIS: Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies. METHODS: Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4+ T cell, CD8+ T cell and CD4-CD8- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing. RESULTS: We identified 79, 56 and 45 differentially methylated regions in CD4+ T cells, CD8+ T cells and CD4-CD8- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4+ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4+ T cells. CONCLUSIONS/INTERPRETATION: These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

Persistent coxsackievirus B1 infection triggers extensive changes in the transcriptome of human pancreatic ductal cells.

Enteroviruses, particularly the group B coxsackieviruses (CVBs), have been associated with the development of type 1 diabetes. Several CVB serotypes establish chronic infections in human cells in vivo and in vitro. However, the mechanisms leading to enterovirus persistency and, possibly, beta cell autoimmunity are not fully understood. We established a carrier-state-type persistent infection model in human pancreatic cell line PANC-1 using two distinct CVB1 strains and profiled the infection-induced changes in cellular transcriptome. In the current study, we observed clear changes in the gene expression of factors associated with the pancreatic microenvironment, the secretory pathway, and lysosomal biogenesis during persistent CVB1 infections. Moreover, we found that the antiviral response pathways were activated differently by the two CVB1 strains. Overall, our study reveals extensive transcriptional responses in persistently CVB1-infected pancreatic cells with strong opposite but also common changes between the two strains.