Multi-omics analysis of early molecular mechanisms of type 1 diabetes

Type 1 diabetes (T1D) is a complicated autoimmune disease with largely unknown disease mechanisms. The diagnosis is preceded by a long asymptomatic period of autoimmune activity in the insulin-producing pancreatic islets. Currently the only clinical markers used for T1D prediction are islet autoantibodies, which are a sign of already-broken immune tolerance. The focus of this dissertation is on the early asymptomatic period preceding seroconversion to islet autoantibody positivity. The genetic risk of type 1 diabetes has been thoroughly mapped in genome-wide association studies, but environmental factors and molecular mechanisms that mediate the risk are less well understood. According to the hygiene hypothesis, the risk of immune-mediated disorders is increased by the lack of exposure to pathogens in modern environments. Within a study on the hygiene hypothesis, we compared umbilical cord blood gene expression patterns between children born in environments with contrasting standards of living and type 1 diabetes incidences (Finland, Russia, and Estonia). The differentially expressed genes were associated with innate immunity and immune maturation. Our results suggest that the environment influences the immune system development already in-utero. Furthermore, we analyzed genome-wide DNA methylation and gene expression profiles in samples collected prospectively from Finnish children and newborn infants at risk of type 1 diabetes. Bisulfite sequencing analysis did not show any association of neonatal DNA methylation with later progression to T1D. However, antiviral type I interferon response in early childhood was found to be a risk factor of T1D. This transcriptomic signature was detectable in the peripheral blood already before islet autoantibodies, and the main observations were confirmed in an independent German study. These results contributed to the hypothesis that virus infections might play a role in T1D. Additionally, this dissertation contributed to transcriptomic and epigenomic data analysis workflows. Simple probe-level analysis of exon array data was shown to improve the reproducibility, specificity, and sensitivity of detected differential exon inclusion events. Type 1 error rate was markedly reduced by permutation-based significance assessment of differential methylation in bisulfite sequencing studies.Tyypin 1 diabeteksen varhaisten molekulaaristen mekanismien multiomiikka-analyysi Tyypin 1 diabetes (T1D) on autoimmuunitauti, jonka taustalla olevista mekanismeista tiedetään vähän. Diagnoosia edeltää pitkä oireeton jakso, jonka aikana insuliinia tuottaviin beetasoluihin kohdistuva autoimmuunireaktio etenee haiman saarekkeissa. Tämä väitöskirjatutkimus keskittyy T1D:n varhaiseen oireettomaan ajanjaksoon, joka edeltää serokonversiota autovasta-ainepositiiviseksi. Tyypin 1 diabeteksen geneettiset riskitekijät on kartoitettu perusteellisesti genominlaajuisissa assosiaatiotutkimuksissa, mutta ympäristön riskitekijöistä ja riskiä välittävistä molekyylimekanismeista tiedetään vähemmän. Hygieniahypoteesin mukaan vähäinen altistuminen taudinaiheuttajille lisää immuunijärjestelmän häiriöiden riskiä. Hygieniahypoteesiin liittyvässä osatyössä vertasimme hygienian ja T1D:n ilmaantuvuuden suhteen erilaisissa ympäristöissä (Suomi, Venäjä ja Viro) syntyneiden lasten napaveren geeniekpressioprofiileja. Erilaisesti ekspressoituneet geenit liittyivät synnynnäiseen immuniteettiin ja immuunijärjestelmän maturaatioon. Näiden tulosten perusteella ympäristö saattaa vaikuttaa immuunijärjestelmän kehitykseen jo raskauden aikana. Genominlaajuista DNA-metylaatiota ja geeniekspressiota analysoitiin näytteistä, jotka oli kerätty laajassa suomalaisessa seurantatutkimuksessa T1D:n riskiryhmään kuuluvilta lapsilta ja vastasyntyneiltä. Bisulfiittisekvensointianalyysin perusteella vastasyntyneen DNA-metylaation ja lapsuuden aikana kehittyvän T1D:n välillä ei ollut yhteyttä. Sen sijaan RNA:n tasolla havaittava viruksiin kohdistuva tyypin 1 interferonivaste varhaislapsuudessa todettiin T1D:n riskitekijäksi. Tämä havainto tehtiin perifeerisestä verestä jo ennen saarekevasta-aineiden ilmaantumista, ja päähavainnot vahvistettiin saksalaisessa tutkimuksessa. Nämä tulokset vahvistivat hypoteesia, jonka mukaan virukset voivat vaikuttaa T1D:n puhkeamiseen. T1D-tutkimuksen ohella tämä väitöskirjatyö kehitti transkriptomiikkaan ja epigenomiikkaan sopivia analyysimenetelmiä. Eksonimikrosirujen koetintasoisen analyysin todettiin parantavan toistettavuutta, sensitiivisyyttä ja tarkkuutta vaihtoehtoisen silmukoinniin kartoittamisessa. Tilastollisen merkitsevyyden permutaatiopohjainen analyysi vähensi tyypin 1 virhettä bisulfiittisekvensointidatan analyysissa

The tumor and plasma cytokine profiles of renal cell carcinoma patients

Renal cell carcinoma (RCC) accounts for 90% of all renal cancers and is considered highly immunogenic. Although many studies have reported the circulating peripheral cytokine profiles, the signatures between the tumor tissue and matching healthy adjacent renal tissue counterparts have not been explored. We aimed to comprehensively investigate the cytokine landscape of RCC tumors and its correlation between the amount and phenotype of the tumor infiltrating lymphocytes (TILs). We analyzed the secretion of 42 cytokines from the tumor (n = 46), adjacent healthy kidney tissues (n = 23) and matching plasma samples (n = 33) with a Luminex-based assay. We further explored the differences between the tissue types, as well as correlated the findings with clinical data and detailed immunophenotyping of the TILs. Using an unsupervised clustering approach, we observed distinct differences in the cytokine profiles between the tumor and adjacent renal tissue samples. The tumor samples clustered into three distinct profiles based on the cytokine expressions: high (52.2% of the tumors), intermediate (26.1%), and low (21.7%). Most of the tumor cytokines positively correlated with each other, except for IL-8 that showed no correlation with any of the measured cytokine expressions. Furthermore, the quantity of lymphocytes in the tumor samples analyzed with flow cytometry positively correlated with the chemokine-family of cytokines, CXCL10 (IP-10) and CXCL9 (MIG). No significant correlations were found between the tumor and matching plasma cytokines, suggesting that circulating cytokines poorly mirror the tumor cytokine environment. Our study highlights distinct cytokine profiles in the RCC tumor microenvironment and provides insights to potential biomarkers for the treatment of RCC.Peer reviewe

The tumor and plasma cytokine profiles of renal cell carcinoma patients

Renal cell carcinoma (RCC) accounts for 90% of all renal cancers and is considered highly immunogenic. Although many studies have reported the circulating peripheral cytokine profiles, the signatures between the tumor tissue and matching healthy adjacent renal tissue counterparts have not been explored. We aimed to comprehensively investigate the cytokine landscape of RCC tumors and its correlation between the amount and phenotype of the tumor infiltrating lymphocytes (TILs). We analyzed the secretion of 42 cytokines from the tumor (n = 46), adjacent healthy kidney tissues (n = 23) and matching plasma samples (n = 33) with a Luminex-based assay. We further explored the differences between the tissue types, as well as correlated the findings with clinical data and detailed immunophenotyping of the TILs. Using an unsupervised clustering approach, we observed distinct differences in the cytokine profiles between the tumor and adjacent renal tissue samples. The tumor samples clustered into three distinct profiles based on the cytokine expressions: high (52.2% of the tumors), intermediate (26.1%), and low (21.7%). Most of the tumor cytokines positively correlated with each other, except for IL-8 that showed no correlation with any of the measured cytokine expressions. Furthermore, the quantity of lymphocytes in the tumor samples analyzed with flow cytometry positively correlated with the chemokine-family of cytokines, CXCL10 (IP-10) and CXCL9 (MIG). No significant correlations were found between the tumor and matching plasma cytokines, suggesting that circulating cytokines poorly mirror the tumor cytokine environment. Our study highlights distinct cytokine profiles in the RCC tumor microenvironment and provides insights to potential biomarkers for the treatment of RCC

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.Peer reviewe

Characterization and non-parametric modeling of the developing serum proteome during infancy and early childhood

Children develop rapidly during the first years of life, and understanding the sources and associated levels of variation in the serum proteome is important when using serum proteins as markers for childhood diseases. The aim of this study was to establish a reference model for the evolution of a healthy serum proteome during early childhood. Label-free quantitative proteomics analyses were performed for 103 longitudinal serum samples collected from 15 children at birth and between the ages of 3-36 months. A flexible Gaussian process-based probabilistic modelling framework was developed to evaluate the effects of different variables, including age, living environment and individual variation, on the longitudinal expression profiles of 266 reliably identified and quantified serum proteins. Age was the most dominant factor influencing approximately half of the studied proteins, and the most prominent age-associated changes were observed already during the first year of life. High inter-individual variability was also observed for multiple proteins. These data provide important details on the maturing serum proteome during early life, and evaluate how patterns detected in cord blood are conserved in the first years of life. Additionally, our novel modelling approach provides a statistical framework to detect associations between covariates and non-linear time series data

Characterization and non-parametric modeling of the developing serum proteome during infancy and early childhood

Children develop rapidly during the first years of life, and understanding the sources and associated levels of variation in the serum proteome is important when using serum proteins as markers for childhood diseases. The aim of this study was to establish a reference model for the evolution of a healthy serum proteome during early childhood. Label-free quantitative proteomics analyses were performed for 103 longitudinal serum samples collected from 15 children at birth and between the ages of 3–36 months. A flexible Gaussian process-based probabilistic modelling framework was developed to evaluate the effects of different variables, including age, living environment and individual variation, on the longitudinal expression profiles of 266 reliably identified and quantified serum proteins. Age was the most dominant factor influencing approximately half of the studied proteins, and the most prominent age-associated changes were observed already during the first year of life. High inter-individual variability was also observed for multiple proteins. These data provide important details on the maturing serum proteome during early life, and evaluate how patterns detected in cord blood are conserved in the first years of life. Additionally, our novel modelling approach provides a statistical framework to detect associations between covariates and non-linear time series data

Permutation-based significance analysis reduces the type 1 error rate in bisulfite 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.Peer reviewe

Genome-wide Analysis of STAT3-Mediated Transcription during Early Human Th17 Cell Differentiation

The development of therapeutic strategies to combat immune-associated diseases requires the molecular mechanisms of human Th17 cell differentiation to be fully identified and understood. To investigate transcriptional control of Th17 cell differentiation, we used primary human CD4+ T cells in small interfering RNA (siRNA)-mediated gene silencing and chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq) to identify both the early direct and indirect targets of STAT3. The integrated dataset presented in this study confirms that STAT3 is critical for transcriptional regulation of early human Th17 cell differentiation. Additionally, we found that a number of SNPs from loci associated with immune-mediated disorders were located at sites where STAT3 binds to induce Th17 cell specification. Importantly, introduction of such SNPs alters STAT3 binding in DNA affinity precipitation assays. Overall, our study provides important insights for modulating Th17-mediated pathogenic immune responses in humans.</p