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.

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.

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).

Computational solutions for spatial transcriptomics

Transcriptome level expression data connected to the spatial organization of the cells and molecules would allow a comprehensive understanding of how gene expression is connected to the structure and function in the biological systems. The spatial transcriptomics platforms may soon provide such information. However, the current platforms still lack spatial resolution, capture only a fraction of the transcriptome heterogeneity, or lack the throughput for large scale studies. The strengths and weaknesses in current ST platforms and computational solutions need to be taken into account when planning spatial transcriptomics studies. The basis of the computational ST analysis is the solutions developed for single-cell RNA-sequencing data, with advancements taking into account the spatial connectedness of the transcriptomes. The scRNA-seq tools are modified for spatial transcriptomics or new solutions like deep learning-based joint analysis of expression, spatial, and image data are developed to extract biological information in the spatially resolved transcriptomes. The computational ST analysis can reveal remarkable biological insights into spatial patterns of gene expression, cell signaling, and cell type variations in connection with cell type-specific signaling and organization in complex tissues. This review covers the topics that help choosing the platform and computational solutions for spatial transcriptomics research. We focus on the currently available ST methods and platforms and their strengths and limitations. Of the computational solutions, we provide an overview of the analysis steps and tools used in the ST data analysis. The compatibility with the data types and the tools provided by the current ST analysis frameworks are summarized.

COVID-19-specific transcriptomic signature detectable in blood across multiple cohorts

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading across the world despite vast global vaccination efforts. Consequently, many studies have looked for potential human host factors and immune mechanisms associated with the disease. However, most studies have focused on comparing COVID-19 patients to healthy controls, while fewer have elucidated the specific host factors distinguishing COVID-19 from other infections. To discover genes specifically related to COVID-19, we reanalyzed transcriptome data from nine independent cohort studies, covering multiple infections, including COVID-19, influenza, seasonal coronaviruses, and bacterial pneumonia. The identified COVID-19-specific signature consisted of 149 genes, involving many signals previously associated with the disease, such as induction of a strong immunoglobulin response and hemostasis, as well as dysregulation of cell cycle-related processes. Additionally, potential new gene candidates related to COVID-19 were discovered. To facilitate exploration of the signature with respect to disease severity, disease progression, and different cell types, we also offer an online tool for easy visualization of the selected genes across multiple datasets at both bulk and single-cell levels.

Introducing untargeted data-independent acquisition for metaproteomics of complex microbial samples

Mass spectrometry-based metaproteomics is a relatively new field of research that enables the characterization of the functionality of microbiota. Recently, we demonstrated the applicability of data-independent acquisition (DIA) mass spectrometry to the analysis of complex metaproteomic samples. This allowed us to circumvent many of the drawbacks of the previously used data-dependent acquisition (DDA) mass spectrometry, mainly the limited reproducibility when analyzing samples with complex microbial composition. However, the DDA-assisted DIA approach still required additional DDA data on the samples to assist the analysis. Here, we introduce, for the first time, an untargeted DIA metaproteomics tool that does not require any DDA data, but instead generates a pseudospectral library directly from the DIA data. This reduces the amount of required mass spectrometry data to a single DIA run per sample. The new DIA-only metaproteomics approach is implemented as a new open-source software package named glaDIAtor, including a modern web-based graphical user interface to facilitate wide use of the tool by the community.

PhosPiR: an automated phosphoproteomic pipeline in R

Large-scale phosphoproteome profiling using mass spectrometry (MS) provides functional insight that is crucial for disease biology and drug discovery. However, extracting biological understanding from these data is an arduous task requiring multiple analysis platforms that are not adapted for automated high-dimensional data analysis. Here, we introduce an integrated pipeline that combines several R packages to extract high-level biological understanding from large-scale phosphoproteomic data by seamless integration with existing databases and knowledge resources. In a single run, PhosPiR provides data clean-up, fast data overview, multiple statistical testing, differential expression analysis, phosphosite annotation and translation across species, multilevel enrichment analyses, proteome-wide kinase activity and substrate mapping and network hub analysis. Data output includes graphical formats such as heatmap, box-, volcano- and circos-plots. This resource is designed to assist proteome-wide data mining of pathophysiological mechanism without a need for programming knowledge.

Improved risk prediction of chemotherapy-induced neutropenia-model development and validation with real-world data

BACKGROUND

The existing risk prediction models for chemotherapy-induced febrile neutropenia (FN) do not necessarily apply to real-life patients in different healthcare systems and the external validation of these models are often lacking. Our study evaluates whether a machine learning-based risk prediction model could outperform the previously introduced models, especially when validated against real-world patient data from another institution not used for model training.

METHODS

Using Turku University Hospital electronic medical records, we identified all patients who received chemotherapy for non-hematological cancer between the years 2010 and 2017 (N = 5879). An experimental surrogate endpoint was first-cycle neutropenic infection (NI), defined as grade IV neutropenia with serum C-reactive protein >10 mg/l. For predicting the risk of NI, a penalized regression model (Lasso) was developed. The model was externally validated in an independent dataset (N = 4594) from Tampere University Hospital.

RESULTS

Lasso model accurately predicted NI risk with good accuracy (AUROC 0.84). In the validation cohort, the Lasso model outperformed two previously introduced, widely approved models, with AUROC 0.75. The variables selected by Lasso included granulocyte colony-stimulating factor (G-CSF) use, cancer type, pre-treatment neutrophil and thrombocyte count, intravenous treatment regimen, and the planned dose intensity. The same model predicted also FN, with AUROC 0.77, supporting the validity of NI as an endpoint.

CONCLUSIONS

Our study demonstrates that real-world NI risk prediction can be improved with machine learning and that every difference in patient or treatment characteristics can have a significant impact on model performance. Here we outline a novel, externally validated approach which may hold potential to facilitate more targeted use of G-CSFs in the future.

Compressive stress-mediated p38 activation required for ERα + phenotype in breast cancer

Breast cancer is now globally the most frequent cancer and leading cause of women's death. Two thirds of breast cancers express the luminal estrogen receptor-positive (ERα + ) phenotype that is initially responsive to antihormonal therapies, but drug resistance emerges. A major barrier to the understanding of the ERα-pathway biology and therapeutic discoveries is the restricted repertoire of luminal ERα + breast cancer models. The ERα + phenotype is not stable in cultured cells for reasons not fully understood. We examine 400 patient-derived breast epithelial and breast cancer explant cultures (PDECs) grown in various three-dimensional matrix scaffolds, finding that ERα is primarily regulated by the matrix stiffness. Matrix stiffness upregulates the ERα signaling via stress-mediated p38 activation and H3K27me3-mediated epigenetic regulation. The finding that the matrix stiffness is a central cue to the ERα phenotype reveals a mechanobiological component in breast tissue hormonal signaling and enables the development of novel therapeutic interventions. Subject terms: ER-positive (ER + ), breast cancer, ex vivo model, preclinical model, PDEC, stiffness, p38 SAPK.

SimPhospho: a software tool enabling confident phosphosite assignment

Motivation

Mass spectrometry combined with enrichment strategies for phosphorylated peptides has been successfully employed for two decades to identify sites of phosphorylation. However, unambiguous phosphosite assignment is considered challenging. Given that site-specific phosphorylation events function as different molecular switches, validation of phosphorylation sites is of utmost importance. In our earlier study we developed a method based on simulated phosphopeptide spectral libraries, which enables highly sensitive and accurate phosphosite assignments. To promote more widespread use of this method, we here introduce a software implementation with improved usability and performance.

Results

We present SimPhospho, a fast and user-friendly tool for accurate simulation of phosphopeptide tandem mass spectra. Simulated phosphopeptide spectral libraries are used to validate and supplement database search results, with a goal to improve reliable phosphoproteome identification and reporting. The presented program can be easily used together with the Trans-Proteomic Pipeline and integrated in a phosphoproteomics data analysis workflow.

Availability and implementation

SimPhospho is open source and it is available for Windows, Linux and Mac operating systems. The software and its user’s manual with detailed description of data analysis as well as test data can be found at https://sourceforge.net/projects/simphospho/.

Supplementary information

Supplementary data are available at Bioinformatics online.

A comprehensive evaluation of popular proteomics software workflows for label-free proteome quantification and imputation

Label-free mass spectrometry (MS) has developed into an important tool applied in various fields of biological and life sciences. Several software exist to process the raw MS data into quantified protein abundances, including open source and commercial solutions. Each software includes a set of unique algorithms for different tasks of the MS data processing workflow. While many of these algorithms have been compared separately, a thorough and systematic evaluation of their overall performance is missing. Moreover, systematic information is lacking about the amount of missing values produced by the different proteomics software and the capabilities of different data imputation methods to account for them.In this study, we evaluated the performance of five popular quantitative label-free proteomics software workflows using four different spike-in data sets. Our extensive testing included the number of proteins quantified and the number of missing values produced by each workflow, the accuracy of detecting differential expression and logarithmic fold change and the effect of different imputation and filtering methods on the differential expression results. We found that the Progenesis software performed consistently well in the differential expression analysis and produced few missing values. The missing values produced by the other software decreased their performance, but this difference could be mitigated using proper data filtering or imputation methods. Among the imputation methods, we found that the local least squares (lls) regression imputation consistently increased the performance of the software in the differential expression analysis, and a combination of both data filtering and local least squares imputation increased performance the most in the tested data sets.

A Data Analysis Protocol for Quantitative Data-Independent Acquisition Proteomics

Data-independent acquisition (DIA) mode of mass spectrometry, such as the SWATH-MS technology, enables accurate and consistent measurement of proteins, which is crucial for comparative proteomics studies. However, there is lack of free and easy to implement data analysis protocols that can handle the different data processing steps from raw spectrum files to peptide intensity matrix and its downstream analysis. Here, we provide a data analysis protocol, named diatools, covering all these steps from spectral library building to differential expression analysis of DIA proteomics data. The data analysis tools used in this protocol are open source and the protocol is distributed at Docker Hub as a complete software environment that supports Linux, Windows, and macOS operating systems.

A systematic evaluation of normalization methods in quantitative label-free proteomics

To date, mass spectrometry (MS) data remain inherently biased as a result of reasons ranging from sample handling to differences caused by the instrumentation. Normalization is the process that aims to account for the bias and make samples more comparable. The selection of a proper normalization method is a pivotal task for the reliability of the downstream analysis and results. Many normalization methods commonly used in proteomics have been adapted from the DNA microarray techniques. Previous studies comparing normalization methods in proteomics have focused mainly on intragroup variation. In this study, several popular and widely used normalization methods representing different strategies in normalization are evaluated using three spike-in and one experimental mouse label-free proteomic data sets. The normalization methods are evaluated in terms of their ability to reduce variation between technical replicates, their effect on differential expression analysis and their effect on the estimation of logarithmic fold changes. Additionally, we examined whether normalizing the whole data globally or in segments for the differential expression analysis has an effect on the performance of the normalization methods. We found that variance stabilization normalization (Vsn) reduced variation the most between technical replicates in all examined data sets. Vsn also performed consistently well in the differential expression analysis. Linear regression normalization and local regression normalization performed also systematically well. Finally, we discuss the choice of a normalization method and some qualities of a suitable normalization method in the light of the results of our evaluation.

Phosphonormalizer: an R package for normalization of MS-based label-free phosphoproteomics

Motivation

Global centering-based normalization is a commonly used normalization approach in mass spectrometry-based label-free proteomics. It scales the peptide abundances to have the same median intensities, based on an assumption that the majority of abundances remain the same across the samples. However, especially in phosphoproteomics, this assumption can introduce bias, as the samples are enriched during sample preparation which can mask the underlying biological changes. To address this possible bias, phosphopeptides quantified in both enriched and non-enriched samples can be used to calculate factors that mitigate the bias.

Results

We present an R package phosphonormalizer for normalizing enriched samples in label-free mass spectrometry-based phosphoproteomics.

Availability and implementation

The phosphonormalizer package is freely available under GPL ( > =2) license from Bioconductor (https://bioconductor.org/packages/phosphonormalizer).

Contact

sohrab.saraei@utu.fi or laura.elo@utu.fi.

Supplementary information

Supplementary data are available at Bioinformatics online.

HIF prolyl hydroxylase PHD3 regulates translational machinery and glucose metabolism in clear cell renal cell carcinoma

Background

A key feature of clear cell renal cell carcinoma (ccRCC) is the inactivation of the von Hippel-Lindau tumour suppressor protein (pVHL) that leads to the activation of hypoxia-inducible factor (HIF) pathway also in well-oxygenated conditions. Important regulator of HIF-α, prolyl hydroxylase PHD3, is expressed in high amounts in ccRCC. Although several functions and downstream targets for PHD3 in cancer have been suggested, the role of elevated PHD3 expression in ccRCC is not clear.

Methods

To gain insight into the functions of high PHD3 expression in ccRCC, we used PHD3 knockdown by siRNA in 786-O cells under normoxic and hypoxic conditions and performed discovery mass spectrometry (LC-MS/MS) of the purified peptide samples. The LC-MS/MS results were analysed by label-free quantification of proteome data using a peptide-level expression-change averaging procedure and subsequent gene ontology enrichment analysis.

Results

Our data reveals an intriguingly widespread effect of PHD3 knockdown with 91 significantly regulated proteins. Under hypoxia, the response to PHD3 silencing was wider than under normoxia illustrated by both the number of regulated proteins and by the range of protein expression levels. The main cellular functions regulated by PHD3 expression were glucose metabolism, protein translation and messenger RNA (mRNA) processing. PHD3 silencing led to downregulation of most glycolytic enzymes from glucose transport to lactate production supported by the reduction in extracellular acidification and lactate production and increase in cellular oxygen consumption rate. Moreover, upregulation of mRNA processing-related proteins and alteration in a number of ribosomal proteins was seen as a response to PHD3 silencing. Further studies on upstream effectors of the translational machinery revealed a possible role for PHD3 in regulation of mTOR pathway signalling.

Conclusions

Our findings suggest crucial involvement of PHD3 in the maintenance of key cellular functions including glycolysis and protein synthesis in ccRCC.

Electronic supplementary material

The online version of this article (doi:10.1186/s40170-017-0167-y) contains supplementary material, which is available to authorized users.

ROTS: An R package for reproducibility-optimized statistical testing

Differential expression analysis is one of the most common types of analyses performed on various biological data (e.g. RNA-seq or mass spectrometry proteomics). It is the process that detects features, such as genes or proteins, showing statistically significant differences between the sample groups under comparison. A major challenge in the analysis is the choice of an appropriate test statistic, as different statistics have been shown to perform well in different datasets. To this end, the reproducibility-optimized test statistic (ROTS) adjusts a modified t-statistic according to the inherent properties of the data and provides a ranking of the features based on their statistical evidence for differential expression between two groups. ROTS has already been successfully applied in a range of different studies from transcriptomics to proteomics, showing competitive performance against other state-of-the-art methods. To promote its widespread use, we introduce here a Bioconductor R package for performing ROTS analysis conveniently on different types of omics data. To illustrate the benefits of ROTS in various applications, we present three case studies, involving proteomics and RNA-seq data from public repositories, including both bulk and single cell data. The package is freely available from Bioconductor (https://www.bioconductor.org/packages/ROTS).

Brief Isoflurane Anesthesia Produces Prominent Phosphoproteomic Changes in the Adult Mouse Hippocampus

Anesthetics are widely used in medical practice and experimental research, yet the neurobiological basis governing their effects remains obscure. We have here used quantitative phosphoproteomics to investigate the protein phosphorylation changes produced by a 30 min isoflurane anesthesia in the adult mouse hippocampus. Altogether 318 phosphorylation alterations in total of 237 proteins between sham and isoflurane anesthesia were identified. Many of the hit proteins represent primary pharmacological targets of anesthetics. However, findings also enlighten the role of several other proteins-implicated in various biological processes including neuronal excitability, brain energy homeostasis, synaptic plasticity and transmission, and microtubule function-as putative (secondary) targets of anesthetics. In particular, isoflurane increases glycogen synthase kinase-3β (GSK3β) phosphorylation at the inhibitory Ser(9) residue and regulates the phosphorylation of multiple proteins downstream and upstream of this promiscuous kinase that regulate diverse biological functions. Along with confirmatory Western blot data for GSK3β and p44/42-MAPK (mitogen-activated protein kinase; reduced phosphorylation of the activation loop), we observed increased phosphorylation of microtubule-associated protein 2 (MAP2) on residues (Thr(1620,1623)) that have been shown to render its dissociation from microtubules and alterations in microtubule stability. We further demonstrate that diverse anesthetics (sevoflurane, urethane, ketamine) produce essentially similar phosphorylation changes on GSK3β, p44/p42-MAPK, and MAP2 as observed with isoflurane. Altogether our study demonstrates the potential of quantitative phosphoproteomics to study the mechanisms of anesthetics (and other drugs) in the mammalian brain and reveals how already a relatively brief anesthesia produces pronounced phosphorylation changes in multiple proteins in the central nervous system.

Using Peptide-Level Proteomics Data for Detecting Differentially Expressed Proteins

The expression of proteins can be quantified in high-throughput means using different types of mass spectrometers. In recent years, there have emerged label-free methods for determining protein abundance. Although the expression is initially measured at the peptide level, a common approach is to combine the peptide-level measurements into protein-level values before differential expression analysis. However, this simple combination is prone to inconsistencies between peptides and may lose valuable information. To this end, we introduce here a method for detecting differentially expressed proteins by combining peptide-level expression-change statistics. Using controlled spike-in experiments, we show that the approach of averaging peptide-level expression changes yields more accurate lists of differentially expressed proteins than does the conventional protein-level approach. This is particularly true when there are only few replicate samples or the differences between the sample groups are small. The proposed technique is implemented in the Bioconductor package PECA, and it can be downloaded from http://www.bioconductor.org.

Hot spots, indicator taxa, complementarity and optimal networks of taiga

If hot spots for different taxa coincide, priority-setting surveys in a region could be carried out more cheaply by focusing on indicator taxa. Several previous studies show that hot spots of different taxa rarely coincide. However, in tropical areas indicator taxa may be used in selecting complementary networks to represent biodiversity as a whole. We studied beetles (Coleoptera), Heteroptera, polypores or bracket fungi (Polyporaceae) and vascular plants of old growth boreal taiga forests. Optimal networks for Heteroptera maximized the high overall species richness of beetles and vascular plants, but these networks were least favourable options for polypores. Polypores are an important group indicating the conservation value of old growth taiga forests. Random selection provided a better option. Thus, certain groups may function as good indicators for maximizing the overall species richness of some taxonomic groups, but all taxa should be examined separately.