Explore the vast range of titles available.

Insertions and deletions (indels) in human genomes are associated with a wide range of phenotypes, including various clinical disorders. High-throughput, next generation sequencing (NGS) technologies enable the detection of short genetic variants, such as single nucleotide variants (SNVs) and indels. However, the variant calling accuracy for indels remains considerably lower than for SNVs. Here we present a comparative study of the performance of variant calling tools for indel calling, evaluated with a wide repertoire of NGS datasets. While there is no single optimal tool to suit all circumstances, our results demonstrate that the choice of variant calling tool greatly impacts the precision and recall of indel calling. Furthermore, to reliably detect indels, it is essential to choose NGS technologies that offer a long read length and high coverage coupled with specific variant calling tools.

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

Background The use of glucocorticoids has given contradictory results for treating acute respiratory distress syndrome (ARDS). The use of intravenous Interferon beta (IFN β) for the treatment of ARDS was recently tested in a phase III ARDS trial (INTEREST), in which more than half of the patients simultaneously received glucocorticoids. Trial results showed deleterious effects of glucocorticoids when administered together with IFN β, and therefore, we aimed at finding the reason behind this. Methods We first sequenced the genes encoding the IFN α/β receptor of the patients, who participated in the INTEREST study (ClinicalTrials.gov Identifier:  NCT02622724 , November 24, 2015) in which the patients were randomized to receive an intravenous injection of IFN β-1a (144 patients) or placebo (152 patients). Genetic background was analyzed against clinical outcome, concomitant medication, and pro-inflammatory cytokine levels. Thereafter, we tested the influence of the genetic background on IFN α/β receptor expression in lung organ cultures and whether, it has any effect on transcription factors STAT1 and STAT2 involved in IFN signaling. Results We found a novel disease association of a SNP rs9984273, which is situated in the interferon α/β receptor subunit 2 (IFNAR2) gene in an area corresponding to a binding motif of the glucocorticoid receptor (GR). The minor allele of SNP rs9984273 associates with higher IFNAR expression, more rapid decrease of IFN γ and interleukin-6 (IL-6) levels and better outcome in IFN β treated patients with ARDS, while the major allele associates with a poor outcome especially under concomitant IFN β and glucocorticoid treatment. Moreover, the minor allele of rs9984273 associates with a less severe form of coronavirus diseases (COVID-19) according to the COVID-19 Host Genetics Initiative database. Conclusions The distribution of this SNP within clinical study arms may explain the contradictory results of multiple ARDS studies and outcomes in COVID-19 concerning type I IFN signaling and glucocorticoids. Graphical abstract

We describe a new reproducibility-optimization method ROPECA for statistical analysis of proteomics data with a specific focus on the emerging data-independent acquisition (DIA) mass spectrometry technology. ROPECA optimizes the reproducibility of statistical testing on peptide-level and aggregates the peptide-level changes to determine differential protein-level expression. Using a ‘gold standard’ spike-in data and a hybrid proteome benchmark data we show the competitive performance of ROPECA over conventional protein-based analysis as well as state-of-the-art peptide-based tools especially in DIA data with consistent peptide measurements. Furthermore, we also demonstrate the improved accuracy of our method in clinical studies using proteomics data from a longitudinal human twin study.

The sequence contexts of genomic variants play important roles in understanding biological significances of variants and potential sequencing related variant calling issues. However, methods for assessing the diverse sequence contexts of genomic variants such as tandem repeats and unambiguous annotations have been limited. Herein, we describe the Variant Sequence Context Annotation Tool (VarSCAT) for annotating the sequence contexts of genomic variants, including breakpoint ambiguities, flanking bases of variants, wildtype/mutated DNA sequences, variant nomenclatures, distances between adjacent variants, tandem repeat regions, and custom annotation with user customizable options. Our analyses demonstrate that VarSCAT is more versatile and customizable than the currently available methods or strategies for annotating variants in short tandem repeat (STR) regions or insertions and deletions (indels) with breakpoint ambiguity. Variant sequence context annotations of high-confidence human variant sets with VarSCAT revealed that more than 75% of all human individual germline and clinically relevant indels have breakpoint ambiguities. Moreover, we illustrate that more than 80% of human individual germline small variants in STR regions are indels and that the sizes of these indels correlated with STR motif sizes. VarSCAT is available from https://github.com/elolab/VarSCAT .

Pathway analysis is a frequent step in studies involving gene or protein expression data, but most of the available pathway methods are designed for simple case versus control studies of two sample groups without further complexity. The few available methods allowing the pathway analysis of more complex study designs cannot use pathway structures or handle the situation where the variable of interest is not defined for all samples. Such scenarios are common in longitudinal studies with so long follow up time that healthy controls are required to identify the effect of normal aging apart from the effect of disease development, which is not defined for controls. To address the need, we introduce a new method for Pathway Analysis of Longitudinal data (PAL), which is suitable for complex study designs, such as longitudinal data. The main advantages of PAL are the use of pathway structures and the suitability of the approach for study settings beyond currently available tools. We demonstrate the performance of PAL with simulated data and three longitudinal datasets related to the early development of type 1 diabetes, which involve different study designs and only subtle biological signals, and include both transcriptomic and proteomic data. An R package implementing PAL is publicly available at https://github.com/elolab/PAL .

Frequent laboratory monitoring is recommended for early identification of toxicity when initiating conventional synthetic disease-modifying antirheumatic drugs (csDMARDs). We aimed at developing a risk prediction model to individualize laboratory testing at csDMARD initiation. We identified inflammatory joint disease patients (N = 1196) initiating a csDMARD in Turku University Hospital 2013–2019. Baseline and follow-up safety monitoring results were drawn from electronic health records. For rheumatoid arthritis patients, diagnoses and csDMARD initiation/cessation dates were manually confirmed. Primary endpoint was alanine transaminase (ALT) elevation of more than twice the upper limit of normal (ULN) within 6 months after treatment initiation. Computational models for predicting incident ALT elevations were developed using Lasso Cox proportional hazards regression with stable iterative variable selection (SIVS) and were internally validated against a randomly selected test cohort (1/3 of the data) that was not used for training the models. Primary endpoint was reached in 82 patients (6.9%). Among baseline variables, Lasso model with SIVS predicted subsequent ALT elevations of > 2 × ULN using higher ALT, csDMARD other than methotrexate or sulfasalazine and psoriatic arthritis diagnosis as important predictors, with a concordance index of 0.71 in the test cohort. Respectively, at first follow-up, in addition to baseline ALT and psoriatic arthritis diagnosis, also ALT change from baseline was identified as an important predictor resulting in a test concordance index of 0.72. Our computational model predicts ALT elevations after the first follow-up test with good accuracy and can help in optimizing individual testing frequency.

Background Detection of copy number variations (CNVs) from high-throughput next-generation whole-genome sequencing (WGS) data has become a widely used research method during the recent years. However, only a little is known about the applicability of the developed algorithms to ultra-low-coverage (0.0005–0.8×) data that is used in various research and clinical applications, such as digital karyotyping and single-cell CNV detection. Result Here, the performance of six popular read-depth based CNV detection algorithms (BIC-seq2, Canvas, CNVnator, FREEC, HMMcopy, and QDNAseq) was studied using ultra-low-coverage WGS data. Real-world array- and karyotyping kit-based validation were used as a benchmark in the evaluation. Additionally, ultra-low-coverage WGS data was simulated to investigate the ability of the algorithms to identify CNVs in the sex chromosomes and the theoretical minimum coverage at which these tools can accurately function. Our results suggest that while all the methods were able to detect large CNVs, many methods were susceptible to producing false positives when smaller CNVs (< 2 Mbp) were detected. There was also significant variability in their ability to identify CNVs in the sex chromosomes. Overall, BIC-seq2 was found to be the best method in terms of statistical performance. However, its significant drawback was by far the slowest runtime among the methods (> 3 h) compared with FREEC (~ 3 min), which we considered the second-best method. Conclusions Our comparative analysis demonstrates that CNV detection from ultra-low-coverage WGS data can be a highly accurate method for the detection of large copy number variations when their length is in millions of base pairs. These findings facilitate applications that utilize ultra-low-coverage CNV detection.

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. Previous studies have reported heterogeneity in the progression to clinical type 1 diabetes in children who develop either insulin- or glutamic acid decarboxylase-specific antibodies as their first autoantibodies. Here, the authors show that children who later develop disease have distinct characteristics in early immune responses, which are dependent on the type of autoantibodies that appear first.

This study aims to develop and validate a modeling framework to predict long-term weight change on the basis of self-reported weight data. The aim is to enable focusing resources of health systems on individuals that are at risk of not achieving their goals in weight loss interventions, which would help both health professionals and the individuals in weight loss management. The weight loss prediction models were built on 327 participants, aged 21–78, from a Finnish weight coaching cohort, with at least 9 months of self-reported follow-up weight data during weight loss intervention. With these data, we used six machine learning methods to predict weight loss after 9 months and selected the best performing models for implementation as modeling framework. We trained the models to predict either three classes of weight change (weight loss, insufficient weight loss, weight gain) or five classes (high/moderate/insufficient weight loss, high/low weight gain). Finally, the prediction accuracy was validated with an independent cohort of overweight UK adults (n = 184). Of the six tested modeling approaches, logistic regression performed the best. Most three-class prediction models achieved prediction accuracy of > 50% already with half a month of data and up to 97% with 8 months. The five-class prediction models achieved accuracies from 39% (0.5 months) to 89% (8 months). Our approach provides an accurate prediction method for long-term weight loss, with potential for easier and more efficient management of weight loss interventions in the future. A web application is available: https://elolab.shinyapps.io/WeightChangePredictor/ . The trial is registered at clinicaltrials.gov/ct2/show/NCT04019249 (Clinical Trials Identifier NCT04019249), first posted on 15/07/2019.