EpiGEN: an epistasis simulation pipeline

Abstract Summary Simulated data are crucial for evaluating epistasis detection tools in genome-wide association studies. Existing simulators are limited, as they do not account for linkage disequilibrium (LD), support limited interaction models of single nucleotide polymorphisms (SNPs) and only dichotomous phenotypes or depend on proprietary software. In contrast, EpiGEN supports SNP interactions of arbitrary order, produces realistic LD patterns and generates both categorical and quantitative phenotypes. Availability and implementation EpiGEN is implemented in Python 3 and is freely available at https://github.com/baumbachlab/epigen. Supplementary information Supplementary data are available at Bioinformatics online

Alternative splicing impacts microRNA regulation within coding regions

MicroRNAs (miRNAs) are small non-coding RNA molecules that regulate post-transcriptional gene expression by binding to specific target sites. Approximately 95% of human multi-exon genes can be spliced alternatively, which enables the production of functionally diverse transcripts and proteins from a single gene. In complex diseases, such as cancer, gene but also miRNA dysregulation plays a significant role. According to most studies miRNAs preferably bind to 3’-untranslated regions of mRNA. However, through alternative splicing, transcripts might lose exons harboring miRNA target sites and, hence, become unresponsive to miRNA regulation. To check this hypothesis, we studied the role of miRNA target sites in both coding and noncoding regions using six cancer data sets from The Cancer Genome Atlas (TCGA). First, we predicted miRNA target sites on mRNAs from their sequence using TarPmiR. For our analysis, we focused on miRNAs whose expression was negatively correlated with gene expression (as evidence for active regulation) as well as genes that were at least moderately expressed and showed evidence of alternative splicing. We chose different subsets of transcripts to differentiate the effects of target sites in different gene regions. To check whether alternative splicing interferes with miRNA regulation, we trained linear regression models to predict miRNA expression from transcript expression. Using nested models, we compared the predictive power of transcripts with miRNA target sites to that of transcripts without target sites in the investigated gene region. For all six cancer data sets and all subsets, models containing transcripts with target sites predicted miRNA abundance significantly better. We conclude that alternative splicing does interfere with miRNA regulation by skipping exons with miRNA target sites within the coding region.Book of abstract: 4th Belgrade Bioinformatics Conference, June 19-23, 202

From protein-protein to isoform-isoform interactions: the toolkit to map alternative splicing to interactome

Alternative splicing (AS) can impact protein structure and lead to protein-protein interaction (PPI) rewiring. Available PPI networks neglect alternative splicing isoforms: as interactions might happen only between a subset of isoforms, the PPI network contains both false-positive and false-negative interactions. Since it is not feasible to validate all isoform-isoform interactions experimentally, we present a set of tools to investigate AS on a network level: DIGGER to map splicing to the PPI network, as well as NEASE and Spycone to evaluate the functional consequences of network rewiring. DIGGER (https://exbio.wzw.tum.de/digger) integrates PPIs, domain-domain, and residuelevel interactions - the structures that might be spliced in or out and result in interaction gain or loss. Users can explore possible rewiring for an isoform or exon of interest and extract relevant subnetworks. NEASE (https://github.com/louadi/NEASE) identifies pathways that are significantly affected by network rewiring. NEASE extends classic gene set enrichment analysis by considering isoform-specific interactions affecting pathways. Spycone (https://github.com/yollct/spycone) addresses the time-course changes in AS. It searches for isoforms that demonstrate similar temporal splicing patterns and reflect the splicing co-regulation. Spycone further integrates gene set, network, and splicing-aware NEASE enrichment. Overall, we offer a splicing-focused network analysis toolkit that allows for studying the mechanistic consequences of AS.Book of abstract: 4th Belgrade Bioinformatics Conference, June 19-23, 202

New Enhancing MRI Lesions Associate with IL-17, Neutrophil Degranulation and Integrin Microparticles: Multi-Omics Combined with Frequent MRI in Multiple Sclerosis

Background: Blood–barrier (BBB) breakdown and active inflammation are hallmarks of relapsing multiple sclerosis (RMS), but the molecular events contributing to the development of new lesions are not well explored. Leaky endothelial junctions are associated with increased production of endothelial-derived extracellular microvesicles (EVs) and result in the entry of circulating immune cells into the brain. MRI with intravenous gadolinium (Gd) can visualize acute blood–barrier disruption as the initial event of the evolution of new lesions. Methods: Here, weekly MRI with Gd was combined with proteomics, multiplex immunoassay, and endothelial stress-optimized EV array to identify early markers related to BBB disruption. Five patients with RMS with no disease-modifying treatment were monitored weekly using high-resolution 3T MRI scanning with intravenous gadolinium (Gd) for 8 weeks. Patients were then divided into three groups (low, medium, or high MRI activity) defined by the number of new, total, and maximally enhancing Gd-enhancing lesions and the number of new FLAIR lesions. Plasma samples taken at each MRI were analyzed for protein biomarkers of inflammation by quantitative proteomics, and cytokines using multiplex immunoassays. EVs were characterized with an optimized endothelial stress EV array based on exosome surface protein markers for the detection of soluble secreted EVs. Results: Proteomics analysis of plasma yielded quantitative information on 208 proteins at each patient time point (n = 40). We observed the highest number of unique dysregulated proteins (DEPs) and the highest functional enrichment in the low vs. high MRI activity comparison. Complement activation and complement/coagulation cascade were also strongly overrepresented in the low vs. high MRI activity comparison. Activation of the alternative complement pathway, pathways of blood coagulation, extracellular matrix organization, and the regulation of TLR and IGF transport were unique for the low vs. high MRI activity comparison as well, with these pathways being overrepresented in the patient with high MRI activity. Principal component analysis indicated the individuality of plasma profiles in patients. IL-17 was upregulated at all time points during 8 weeks in patients with high vs. low MRI activity. Hierarchical clustering of soluble markers in the plasma indicated that all four MRI outcomes clustered together with IL-17, IL-12p70, and IL-1β. MRI outcomes also showed clustering with EV markers CD62E/P, MIC A/B, ICAM-1, and CD42A. The combined cluster of these cytokines, EV markers, and MRI outcomes clustered also with IL-12p40 and IL-7. All four MRI outcomes correlated positively with levels of IL-17 (p < 0.001, respectively), and EV-ICAM-1 (p < 0.0003, respectively). IL-1β levels positively correlated with the number of new Gd-enhancing lesions (p < 0.01), new FLAIR lesions (p < 0.001), and total number of Gd-enhancing lesions (p < 0.05). IL-6 levels positively correlated with the number of new FLAIR lesions (p < 0.05). Random Forests and linear mixed models identified IL-17, CCL17/TARC, CCL3/MIP-1α, and TNF-α as composite biomarkers predicting new lesion evolution. Conclusions: Combination of serial frequent MRI with proteome, neuroinflammation markers, and protein array data of EVs enabled assessment of temporal changes in inflammation and endothelial dysfunction in RMS related to the evolution of new and enhancing lesions. Particularly, the Th17 pathway and IL-1β clustered and correlated with new lesions and Gd enhancement, indicating their importance in BBB disruption and initiating acute brain inflammation in MS. In addition to the Th17 pathway, abundant protein changes between MRI activity groups suggested the role of EVs and the coagulation system along with innate immune responses including acute phase proteins, complement components, and neutrophil degranulation.Background: Blood-barrier (BBB) breakdown and active inflammation are hallmarks of relapsing multiple sclerosis (RMS), but the molecular events contributing to the development of new lesions are not well explored. Leaky endothelial junctions are associated with increased production of endothelial-derived extracellular microvesicles (EVs) and result in the entry of circulating immune cells into the brain. MRI with intravenous gadolinium (Gd) can visualize acute blood-barrier disruption as the initial event of the evolution of new lesions.Methods: Here, weekly MRI with Gd was combined with proteomics, multiplex immunoassay, and endothelial stress-optimized EV array to identify early markers related to BBB disruption. Five patients with RMS with no disease-modifying treatment were monitored weekly using high-resolution 3T MRI scanning with intravenous gadolinium (Gd) for 8 weeks. Patients were then divided into three groups (low, medium, or high MRI activity) defined by the number of new, total, and maximally enhancing Gd-enhancing lesions and the number of new FLAIR lesions. Plasma samples taken at each MRI were analyzed for protein biomarkers of inflammation by quantitative proteomics, and cytokines using multiplex immunoassays. EVs were characterized with an optimized endothelial stress EV array based on exosome surface protein markers for the detection of soluble secreted EVs.Results: Proteomics analysis of plasma yielded quantitative information on 208 proteins at each patient time point (n = 40). We observed the highest number of unique dysregulated proteins (DEPs) and the highest functional enrichment in the low vs. high MRI activity comparison. Complement activation and complement/coagulation cascade were also strongly overrepresented in the low vs. high MRI activity comparison. Activation of the alternative complement pathway, pathways of blood coagulation, extracellular matrix organization, and the regulation of TLR and IGF transport were unique for the low vs. high MRI activity comparison as well, with these pathways being overrepresented in the patient with high MRI activity. Principal component analysis indicated the individuality of plasma profiles in patients. IL-17 was upregulated at all time points during 8 weeks in patients with high vs. low MRI activity. Hierarchical clustering of soluble markers in the plasma indicated that all four MRI outcomes clustered together with IL-17, IL-12p70, and IL-1β. MRI outcomes also showed clustering with EV markers CD62E/P, MIC A/B, ICAM-1, and CD42A. The combined cluster of these cytokines, EV markers, and MRI outcomes clustered also with IL-12p40 and IL-7. All four MRI outcomes correlated positively with levels of IL-17 (p < 0.001, respectively), and EV-ICAM-1 (p < 0.0003, respectively). IL-1β levels positively correlated with the number of new Gd-enhancing lesions (p < 0.01), new FLAIR lesions (p < 0.001), and total number of Gd-enhancing lesions (p < 0.05). IL-6 levels positively correlated with the number of new FLAIR lesions (p < 0.05). Random Forests and linear mixed models identified IL-17, CCL17/TARC, CCL3/MIP-1α, and TNF-α as composite biomarkers predicting new lesion evolution.Conclusions: Combination of serial frequent MRI with proteome, neuroinflammation markers, and protein array data of EVs enabled assessment of temporal changes in inflammation and endothelial dysfunction in RMS related to the evolution of new and enhancing lesions. Particularly, the Th17 pathway and IL-1β clustered and correlated with new lesions and Gd enhancement, indicating their importance in BBB disruption and initiating acute brain inflammation in MS. In addition to the Th17 pathway, abundant protein changes between MRI activity groups suggested the role of EVs and the coagulation system along with innate immune responses including acute phase proteins, complement components, and neutrophil degranulation

sPLINK : a hybrid federated tool as a robust alternative to meta-analysis in genome-wide association studies

Meta-analysis has been established as an effective approach to combining summary statistics of several genome-wide association studies (GWAS). However, the accuracy of meta-analysis can be attenuated in the presence of cross-study heterogeneity. We present sPLINK, a hybrid federated and user-friendly tool, which performs privacy-aware GWAS on distributed datasets while preserving the accuracy of the results. sPLINK is robust against heterogeneous distributions of data across cohorts while meta-analysis considerably loses accuracy in such scenarios. sPLINK achieves practical runtime and acceptable network usage for chi-square and linear/logistic regression tests.Peer reviewe

Helicobacter pylori infection associates with fecal microbiota composition and diversity

Helicobacter (H.) pylori is the most important cause for peptic ulcer disease and a risk factor for gastric carcinoma. How colonization with H. pylori affects the intestinal microbiota composition in humans is unknown. We investigated the association of H. pylori infection with intestinal microbiota composition in the population-based cohort Study-of-Health-in-Pomerania (SHIP)-TREND. Anti-H. pylori serology and H. pylori stool antigen tests were used to determine the H. pylori infection status. The fecal microbiota composition of 212 H. pylori positive subjects and 212 matched negative control individuals was assessed using 16S rRNA gene sequencing. H. pylori infection was found to be significantly associated with fecal microbiota alterations and a general increase in fecal microbial diversity. In infected individuals, the H. pylori stool antigen load determined a larger portion of the microbial variation than age or sex. The highest H. pylori stool antigen loads were associated with a putatively harmful microbiota composition. This study demonstrates profound alterations in human fecal microbiota of H. pylori infected individuals. While the increased microbiota diversity associated with H. pylori infection as well as changes in abundance of specific genera could be considered to be beneficial, others may be associated with adverse health effects, reflecting the complex relationship between H. pylori and its human host

Plasma proteome and metabolome characterization of an experimental human thyrotoxicosis model.

BACKGROUND: Determinations of thyrotropin (TSH) and free thyroxine (FT4) represent the gold standard in evaluation of thyroid function. To screen for novel peripheral biomarkers of thyroid function and to characterize FT4-associated physiological signatures in human plasma we used an untargeted OMICS approach in a thyrotoxicosis model. METHODS: A sample of 16 healthy young men were treated with levothyroxine for 8 weeks and plasma was sampled before the intake was started as well as at two points during treatment and after its completion, respectively. Mass spectrometry-derived metabolite and protein levels were related to FT4 serum concentrations using mixed-effect linear regression models in a robust setting. To compile a molecular signature discriminating between thyrotoxicosis and euthyroidism, a random forest was trained and validated in a two-stage cross-validation procedure. RESULTS: Despite the absence of obvious clinical symptoms, mass spectrometry analyses detected 65 metabolites and 63 proteins exhibiting significant associations with serum FT4. A subset of 15 molecules allowed a robust and good prediction of thyroid hormone function (AUC = 0.86) without prior information on TSH or FT4. Main FT4-associated signatures indicated increased resting energy expenditure, augmented defense against systemic oxidative stress, decreased lipoprotein particle levels, and increased levels of complement system proteins and coagulation factors. Further association findings question the reliability of kidney function assessment under hyperthyroid conditions and suggest a link between hyperthyroidism and cardiovascular diseases via increased dimethylarginine levels. CONCLUSION: Our results emphasize the power of untargeted OMICs approaches to detect novel pathways of thyroid hormone action. Furthermore, beyond TSH and FT4, we demonstrated the potential of such analyses to identify new molecular signatures for diagnosis and treatment of thyroid disorders. This study was registered at the German Clinical Trials Register (DRKS) [DRKS00011275] on the 16th of November 2016

Multi-omic signature of body weight change: results from a population-based cohort study

BACKGROUND: Excess body weight is a major risk factor for cardiometabolic diseases. The complex molecular mechanisms of body weight change-induced metabolic perturbations are not fully understood. Specifically, in-depth molecular characterization of long-term body weight change in the general population is lacking. Here, we pursued a multi-omic approach to comprehensively study metabolic consequences of body weight change during a seven-year follow-up in a large prospective study. METHODS: We used data from the population-based Cooperative Health Research in the Region of Augsburg (KORA) S4/F4 cohort. At follow-up (F4), two-platform serum metabolomics and whole blood gene expression measurements were obtained for 1,631 and 689 participants, respectively. Using weighted correlation network analysis, omics data were clustered into modules of closely connected molecules, followed by the formation of a partial correlation network from the modules. Association of the omics modules with previous annual percentage weight change was then determined using linear models. In addition, we performed pathway enrichment analyses, stability analyses, and assessed the relation of the omics modules with clinical traits. RESULTS: Four metabolite and two gene expression modules were significantly and stably associated with body weight change (P-values ranging from 1.9 × 10−4 to 1.2 × 10−24). The four metabolite modules covered major branches of metabolism, with VLDL, LDL and large HDL subclasses, triglycerides, branched-chain amino acids and markers of energy metabolism among the main representative molecules. One gene expression module suggests a role of weight change in red blood cell development. The other gene expression module largely overlaps with the lipid-leukocyte (LL) module previously reported to interact with serum metabolites, for which we identify additional co-expressed genes. The omics modules were interrelated and showed cross-sectional associations with clinical traits. Moreover, weight gain and weight loss showed largely opposing associations with the omics modules. CONCLUSIONS: Long-term weight change in the general population globally associates with serum metabolite concentrations. An integrated metabolomics and transcriptomics approach improved the understanding of molecular mechanisms underlying the association of weight gain with changes in lipid and amino acid metabolism, insulin sensitivity, mitochondrial function as well as blood cell development and function

Meta-analysis of human genome-microbiome association studies: the MiBioGen consortium initiative

Abstract Background In recent years, human microbiota, especially gut microbiota, have emerged as an important yet complex trait influencing human metabolism, immunology, and diseases. Many studies are investigating the forces underlying the observed variation, including the human genetic variants that shape human microbiota. Several preliminary genome-wide association studies (GWAS) have been completed, but more are necessary to achieve a fuller picture. Results Here, we announce the MiBioGen consortium initiative, which has assembled 18 population-level cohorts and some 19,000 participants. Its aim is to generate new knowledge for the rapidly developing field of microbiota research. Each cohort has surveyed the gut microbiome via 16S rRNA sequencing and genotyped their participants with full-genome SNP arrays. We have standardized the analytical pipelines for both the microbiota phenotypes and genotypes, and all the data have been processed using identical approaches. Our analysis of microbiome composition shows that we can reduce the potential artifacts introduced by technical differences in generating microbiota data. We are now in the process of benchmarking the association tests and performing meta-analyses of genome-wide associations. All pipeline and summary statistics results will be shared using public data repositories. Conclusion We present the largest consortium to date devoted to microbiota-GWAS. We have adapted our analytical pipelines to suit multi-cohort analyses and expect to gain insight into host-microbiota cross-talk at the genome-wide level. And, as an open consortium, we invite more cohorts to join us (by contacting one of the corresponding authors) and to follow the analytical pipeline we have developed

The AIMe registry for artificial intelligence in biomedical research

peer reviewe