TY - JOUR A1 - Wangorsch, Gaby A1 - Butt, Elke A1 - Mark, Regina A1 - Hubertus, Katharina A1 - Geiger, Jörg A1 - Dandekar, Thomas A1 - Dittrich, Marcus T1 - Time-resolved in silico modeling of fine-tuned cAMP signaling in platelets: feedback loops, titrated phosphorylations and pharmacological modulation N2 - Background: Hemostasis is a critical and active function of the blood mediated by platelets. Therefore, the prevention of pathological platelet aggregation is of great importance as well as of pharmaceutical and medical interest. Endogenous platelet inhibition is predominantly based on cyclic nucleotides (cAMP, cGMP) elevation and subsequent cyclic nucleotide-dependent protein kinase (PKA, PKG) activation. In turn, platelet phosphodiesterases (PDEs) and protein phosphatases counterbalance their activity. This main inhibitory pathway in human platelets is crucial for countervailing unwanted platelet activation. Consequently, the regulators of cyclic nucleotide signaling are of particular interest to pharmacology and therapeutics of atherothrombosis. Modeling of pharmacodynamics allows understanding this intricate signaling and supports the precise description of these pivotal targets for pharmacological modulation. Results: We modeled dynamically concentration-dependent responses of pathway effectors (inhibitors, activators, drug combinations) to cyclic nucleotide signaling as well as to downstream signaling events and verified resulting model predictions by experimental data. Experiments with various cAMP affecting compounds including antiplatelet drugs and their combinations revealed a high fidelity, fine-tuned cAMP signaling in platelets without crosstalk to the cGMP pathway. The model and the data provide evidence for two independent feedback loops: PKA, which is activated by elevated cAMP levels in the platelet, subsequently inhibits adenylyl cyclase (AC) but as well activates PDE3. By multi-experiment fitting, we established a comprehensive dynamic model with one predictive, optimized and validated set of parameters. Different pharmacological conditions (inhibition, activation, drug combinations, permanent and transient perturbations) are successfully tested and simulated, including statistical validation and sensitivity analysis. Downstream cyclic nucleotide signaling events target different phosphorylation sites for cAMP- and cGMP-dependent protein kinases (PKA, PKG) in the vasodilator-stimulated phosphoprotein (VASP). VASP phosphorylation as well as cAMP levels resulting from different drug strengths and combined stimulants were quantitatively modeled. These predictions were again experimentally validated. High sensitivity of the signaling pathway at low concentrations is involved in a fine-tuned balance as well as stable activation of this inhibitory cyclic nucleotide pathway. Conclusions: On the basis of experimental data, literature mining and database screening we established a dynamic in silico model of cyclic nucleotide signaling and probed its signaling sensitivity. Thoroughly validated, it successfully predicts drug combination effects on platelet function, including synergism, antagonism and regulatory loops. KW - Vasodilatator-stimuliertes Phosphoprotein KW - VASP KW - cyclic nucleotide signaling KW - silico model Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-69145 ER - TY - JOUR A1 - Urban, Lara A1 - Remmele, Christian W. A1 - Dittrich, Marcus A1 - Schwarz, Roland F. A1 - Müller, Tobias T1 - covRNA: discovering covariate associations in large-scale gene expression data JF - BMC Reserach Notes N2 - Objective The biological interpretation of gene expression measurements is a challenging task. While ordination methods are routinely used to identify clusters of samples or co-expressed genes, these methods do not take sample or gene annotations into account. We aim to provide a tool that allows users of all backgrounds to assess and visualize the intrinsic correlation structure of complex annotated gene expression data and discover the covariates that jointly affect expression patterns. Results The Bioconductor package covRNA provides a convenient and fast interface for testing and visualizing complex relationships between sample and gene covariates mediated by gene expression data in an entirely unsupervised setting. The relationships between sample and gene covariates are tested by statistical permutation tests and visualized by ordination. The methods are inspired by the fourthcorner and RLQ analyses used in ecological research for the analysis of species abundance data, that we modified to make them suitable for the distributional characteristics of both, RNA-Seq read counts and microarray intensities, and to provide a high-performance parallelized implementation for the analysis of large-scale gene expression data on multi-core computational systems. CovRNA provides additional modules for unsupervised gene filtering and plotting functions to ensure a smooth and coherent analysis workflow. KW - Multivariate analysis KW - Fourthcorner analysis KW - RLQ analysis KW - Transcriptomics KW - High-throughput data KW - Visualization KW - Ordination methods KW - RNA-Seq analysis KW - Microarray analysis Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-229258 VL - 13 ER - TY - JOUR A1 - Staiger, Christine A1 - Cadot, Sidney A1 - Kooter, Raul A1 - Dittrich, Marcus A1 - Müller, Tobias A1 - Klau, Gunnar W. A1 - Wessels, Lodewyk F. A. T1 - A Critical Evaluation of Network and Pathway-Based Classifiers for Outcome Prediction in Breast Cancer JF - PLoS One N2 - Recently, several classifiers that combine primary tumor data, like gene expression data, and secondary data sources, such as protein-protein interaction networks, have been proposed for predicting outcome in breast cancer. In these approaches, new composite features are typically constructed by aggregating the expression levels of several genes. The secondary data sources are employed to guide this aggregation. Although many studies claim that these approaches improve classification performance over single genes classifiers, the gain in performance is difficult to assess. This stems mainly from the fact that different breast cancer data sets and validation procedures are employed to assess the performance. Here we address these issues by employing a large cohort of six breast cancer data sets as benchmark set and by performing an unbiased evaluation of the classification accuracies of the different approaches. Contrary to previous claims, we find that composite feature classifiers do not outperform simple single genes classifiers. We investigate the effect of (1) the number of selected features; (2) the specific gene set from which features are selected; (3) the size of the training set and (4) the heterogeneity of the data set on the performance of composite feature and single genes classifiers. Strikingly, we find that randomization of secondary data sources, which destroys all biological information in these sources, does not result in a deterioration in performance of composite feature classifiers. Finally, we show that when a proper correction for gene set size is performed, the stability of single genes sets is similar to the stability of composite feature sets. Based on these results there is currently no reason to prefer prognostic classifiers based on composite features over single genes classifiers for predicting outcome in breast cancer. KW - modules KW - protein-interaction networks KW - expression signature KW - classification KW - set KW - metastasis KW - stability KW - survival KW - database KW - markers Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-131323 VL - 7 IS - 4 ER - TY - JOUR A1 - Schneider, Eberhard A1 - Dittrich, Marcus A1 - Böck, Julia A1 - Nanda, Indrajit A1 - Müller, Tobias A1 - Seidmann, Larissa A1 - Tralau, Tim A1 - Galetzka, Danuta A1 - El Hajj, Nady A1 - Haaf, Thomas T1 - CpG sites with continuously increasing or decreasing methylation from early to late human fetal brain development JF - Gene N2 - Normal human brain development is dependent on highly dynamic epigenetic processes for spatial and temporal gene regulation. Recent work identified wide-spread changes in DNA methylation during fetal brain development. We profiled CpG methylation in frontal cortex of 27 fetuses from gestational weeks 12-42, using Illumina 450K methylation arrays. Sites showing genome-wide significant correlation with gestational age were compared to a publicly available data set from gestational weeks 3-26. Altogether, we identified 2016 matching developmentally regulated differentially methylated positions (m-dDMPs): 1767 m-dDMPs were hypermethylated and 1149 hypomethylated during fetal development. M-dDMPs are underrepresented in CpG islands and gene promoters, and enriched in gene bodies. They appear to cluster in certain chromosome regions. M-dDMPs are significantly enriched in autism-associated genes and CpGs. Our results promote the idea that reduced methylation dynamics during fetal brain development may predispose to autism. In addition, m-dDMPs are enriched in genes with human-specific brain expression patterns and/or histone modifications. Collectively, we defined a subset of dDMPs exhibiting constant methylation changes from early to late pregnancy. The same epigenetic mechanisms involving methylation changes in cis-regulatory regions may have been adopted for human brain evolution and ontogeny. KW - Autism spectrum disorders KW - DNA methylation KW - Genome KW - Autism KW - Frontal cortex KW - Human prefrontal cortex KW - Gene-expression KW - Schizophrenia KW - Patterns KW - Transcription KW - Epigenetics KW - Environment KW - Fetal brain development KW - DNA methylation dynamics KW - Methylome Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-186936 VL - 592 IS - 1 ER - TY - JOUR A1 - Morton, Charles Oliver A1 - Fliesser, Mirjam A1 - Dittrich, Marcus A1 - Müller, Tobias A1 - Bauer, Ruth A1 - Kneitz, Susanne A1 - Hope, William A1 - Rogers, Thomas Richard A1 - Einsele, Hermann A1 - Löffler, Jürgen T1 - Gene Expression Profiles of Human Dendritic Cells Interacting with Aspergillus fumigatus in a Bilayer Model of the Alveolar Epithelium/Endothelium Interface N2 - The initial stages of the interaction between the host and Aspergillus fumigatus at the alveolar surface of the human lung are critical in the establishment of aspergillosis. Using an in vitro bilayer model of the alveolus, including both the epithelium (human lung adenocarcinoma epithelial cell line, A549) and endothelium (human pulmonary artery epithelial cells, HPAEC) on transwell membranes, it was possible to closely replicate the in vivo conditions. Two distinct sub-groups of dendritic cells (DC), monocyte-derived DC (moDC) and myeloid DC (mDC), were included in the model to examine immune responses to fungal infection at the alveolar surface. RNA in high quantity and quality was extracted from the cell layers on the transwell membrane to allow gene expression analysis using tailored custom-made microarrays, containing probes for 117 immune-relevant genes. This microarray data indicated minimal induction of immune gene expression in A549 alveolar epithelial cells in response to germ tubes of A. fumigatus. In contrast, the addition of DC to the system greatly increased the number of differentially expressed immune genes. moDC exhibited increased expression of genes including CLEC7A, CD209 and CCL18 in the absence of A. fumigatus compared to mDC. In the presence of A. fumigatus, both DC subgroups exhibited up-regulation of genes identified in previous studies as being associated with the exposure of DC to A. fumigatus and exhibiting chemotactic properties for neutrophils, including CXCL2, CXCL5, CCL20, and IL1B. This model closely approximated the human alveolus allowing for an analysis of the host pathogen interface that complements existing animal models of IA. KW - aspergillus fumigatus KW - gene expression KW - immune receptors KW - immune response KW - denritic cells KW - B cell receptors KW - gene regulation KW - RNA extraction Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-112893 ER - TY - JOUR A1 - Maierhofer, Anna A1 - Flunkert, Julia A1 - Oshima, Junko A1 - Martin, George M. A1 - Poot, Martin A1 - Nanda, Indrajit A1 - Dittrich, Marcus A1 - Müller, Tobias A1 - Haaf, Thomas T1 - Epigenetic signatures of Werner syndrome occur early in life and are distinct from normal epigenetic aging processes JF - Aging Cell N2 - Werner Syndrome (WS) is an adult‐onset segmental progeroid syndrome. Bisulfite pyrosequencing of repetitive DNA families revealed comparable blood DNA methylation levels between classical (18 WRN‐mutant) or atypical WS (3 LMNA‐mutant and 3 POLD1‐mutant) patients and age‐ and sex‐matched controls. WS was not associated with either age‐related accelerated global losses of ALU, LINE1, and α‐satellite DNA methylations or gains of rDNA methylation. Single CpG methylation was analyzed with Infinium MethylationEPIC arrays. In a correspondence analysis, atypical WS samples clustered together with the controls and were clearly separated from classical WS, consistent with distinct epigenetic pathologies. In classical WS, we identified 659 differentially methylated regions (DMRs) comprising 3,656 CpG sites and 613 RefSeq genes. The top DMR was located in the HOXA4 promoter. Additional DMR genes included LMNA, POLD1, and 132 genes which have been reported to be differentially expressed in WRN‐mutant/depleted cells. DMRs were enriched in genes with molecular functions linked to transcription factor activity and sequence‐specific DNA binding to promoters transcribed by RNA polymerase II. We propose that transcriptional misregulation of downstream genes by the absence of WRN protein contributes to the variable premature aging phenotypes of WS. There were no CpG sites showing significant differences in DNA methylation changes with age between WS patients and controls. Genes with both WS‐ and age‐related methylation changes exhibited a constant offset of methylation between WRN‐mutant patients and controls across the entire analyzed age range. WS‐specific epigenetic signatures occur early in life and do not simply reflect an acceleration of normal epigenetic aging processes. KW - (classical and atypical) Werner syndrome KW - bisulfite pyrosequencing KW - methylation array KW - premature aging KW - segmental progeria KW - transcription deficiency Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-202733 VL - 18 ER - TY - JOUR A1 - Maierhofer, Anna A1 - Flunkert, Julia A1 - Dittrich, Marcus A1 - Müller, Tobias A1 - Schindler, Detlev A1 - Nanda, Indrajit A1 - Haaf, Thomas T1 - Analysis of global DNA methylation changes in primary human fibroblasts in the early phase following X-ray irradiation JF - PLoS ONE N2 - Epigenetic alterations may contribute to the generation of cancer cells in a multi-step process of tumorigenesis following irradiation of normal body cells. Primary human fibroblasts with intact cell cycle checkpoints were used as a model to test whether X-ray irradiation with 2 and 4 Gray induces direct epigenetic effects (within the first cell cycle) in the exposed cells. ELISA-based fluorometric assays were consistent with slightly reduced global DNA methylation and hydroxymethylation, however the observed between-group differences were usually not significant. Similarly, bisulfite pyrosequencing of interspersed LINE-1 repeats and centromeric α-satellite DNA did not detect significant methylation differences between irradiated and non-irradiated cultures. Methylation of interspersed ALU repeats appeared to be slightly increased (one percentage point; p = 0.01) at 6 h after irradiation with 4 Gy. Single-cell analysis showed comparable variations in repeat methylation among individual cells in both irradiated and control cultures. Radiation-induced changes in global repeat methylation, if any, were much smaller than methylation variation between different fibroblast strains. Interestingly, α-satellite DNA methylation positively correlated with gestational age. Finally, 450K methylation arrays mainly targeting genes and CpG islands were used for global DNA methylation analysis. There were no detectable methylation differences in genic (promoter, 5' UTR, first exon, gene body, 3' UTR) and intergenic regions between irradiated and control fibroblast cultures. Although we cannot exclude minor effects, i.e. on individual CpG sites, collectively our data suggest that global DNA methylation remains rather stable in irradiated normal body cells in the early phase of DNA damage response. KW - DNA methylation KW - fibroblasts KW - methylation KW - alu elements KW - DNA damage KW - epigenetics KW - cancer treatment Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-170895 VL - 12 IS - 5 ER - TY - JOUR A1 - Luther, Christian H. A1 - Brandt, Philipp A1 - Vylkova, Slavena A1 - Dandekar, Thomas A1 - Müller, Tobias A1 - Dittrich, Marcus T1 - Integrated analysis of SR-like protein kinases Sky1 and Sky2 links signaling networks with transcriptional regulation in Candida albicans JF - Frontiers in Cellular and Infection Microbiology N2 - Fungal infections are a major global health burden where Candida albicans is among the most common fungal pathogen in humans and is a common cause of invasive candidiasis. Fungal phenotypes, such as those related to morphology, proliferation and virulence are mainly driven by gene expression, which is primarily regulated by kinase signaling cascades. Serine-arginine (SR) protein kinases are highly conserved among eukaryotes and are involved in major transcriptional processes in human and S. cerevisiae. Candida albicans harbors two SR protein kinases, while Sky2 is important for metabolic adaptation, Sky1 has similar functions as in S. cerevisiae. To investigate the role of these SR kinases for the regulation of transcriptional responses in C. albicans, we performed RNA sequencing of sky1Δ and sky2Δ and integrated a comprehensive phosphoproteome dataset of these mutants. Using a Systems Biology approach, we study transcriptional regulation in the context of kinase signaling networks. Transcriptomic enrichment analysis indicates that pathways involved in the regulation of gene expression are downregulated and mitochondrial processes are upregulated in sky1Δ. In sky2Δ, primarily metabolic processes are affected, especially for arginine, and we observed that arginine-induced hyphae formation is impaired in sky2Δ. In addition, our analysis identifies several transcription factors as potential drivers of the transcriptional response. Among these, a core set is shared between both kinase knockouts, but it appears to regulate different subsets of target genes. To elucidate these diverse regulatory patterns, we created network modules by integrating the data of site-specific protein phosphorylation and gene expression with kinase-substrate predictions and protein-protein interactions. These integrated signaling modules reveal shared parts but also highlight specific patterns characteristic for each kinase. Interestingly, the modules contain many proteins involved in fungal morphogenesis and stress response. Accordingly, experimental phenotyping shows a higher resistance to Hygromycin B for sky1Δ. Thus, our study demonstrates that a combination of computational approaches with integration of experimental data can offer a new systems biological perspective on the complex network of signaling and transcription. With that, the investigation of the interface between signaling and transcriptional regulation in C. albicans provides a deeper insight into how cellular mechanisms can shape the phenotype. KW - sky kinases KW - kinase signaling KW - network analysis KW - transcriptome KW - transcriptional regulation KW - phosphoproteome KW - Candida albicans Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-311771 SN - 2235-2988 VL - 13 ER - TY - JOUR A1 - Kuhn, Joachim A1 - Gripp, Tatjana A1 - Flieder, Tobias A1 - Dittrich, Marcus A1 - Hendig, Doris A1 - Busse, Jessica A1 - Knabbe, Cornelius A1 - Birschmann, Ingvild T1 - UPLC-MRM Mass Spectrometry Method for Measurement of the Coagulation Inhibitors Dabigatran and Rivaroxaban in Human Plasma and Its Comparison with Functional Assays JF - PLOS ONE N2 - Introduction The fast, precise, and accurate measurement of the new generation of oral anticoagulants such as dabigatran and rivaroxaban in patients' plasma my provide important information in different clinical circumstances such as in the case of suspicion of overdose, when patients switch from existing oral anticoagulant, in patients with hepatic or renal impairment, by concomitant use of interaction drugs, or to assess anticoagulant concentration in patients' blood before major surgery. Methods Here, we describe a quick and precise method to measure the coagulation inhibitors dabigatran and rivaroxaban using ultra-performance liquid chromatography electrospray ionization-tandem mass spectrometry in multiple reactions monitoring (MRM) mode (UPLC-MRM MS). Internal standards (ISs) were added to the sample and after protein precipitation; the sample was separated on a reverse phase column. After ionization of the analytes the ions were detected using electrospray ionization-tandem mass spectrometry. Run time was 2.5 minutes per injection. Ion suppression was characterized by means of post-column infusion. Results The calibration curves of dabigatran and rivaroxaban were linear over the working range between 0.8 and 800 mu g/L (r > 0.99). Limits of detection (LOD) in the plasma matrix were 0.21 mu g/L for dabigatran and 0.34 mu g/L for rivaroxaban, and lower limits of quantification (LLOQ) in the plasma matrix were 0.46 mu g/L for dabigatran and 0.54 mu g/L for rivaroxaban. The intraassay coefficients of variation (CVs) for dabigatran and rivaroxaban were < 4% and 6%; respectively, the interassay CVs were < 6% for dabigatran and < 9% for rivaroxaban. Inaccuracy was < 5% for both substances. The mean recovery was 104.5% (range 83.8-113.0%) for dabigatran and 87.0%(range 73.6-105.4%) for rivaroxaban. No significant ion suppressions were detected at the elution times of dabigatran or rivaroxaban. Both coagulation inhibitors were stable in citrate plasma at -20 degrees C, 4 degrees C and even at RT for at least one week. A method comparison between our UPLC-MRM MS method, the commercially available automated Direct Thrombin Inhibitor assay (DTI assay) for dabigatran measurement from CoaChrom Diagnostica, as well as the automated anti-Xa assay for rivaroxaban measurement from Chromogenix both performed by ACL-TOP showed a high degree of correlation. However, UPLC-MRM MS measurement of dabigatran and rivaroxaban has a much better selectivity than classical functional assays measuring activities of various coagulation factors which are susceptible to interference by other coagulant drugs. Conclusions Overall, we developed and validated a sensitive and specific UPLC-MRM MS assay for the quick and specific measurement of dabigatran and rivaroxaban in human plasma. KW - LC-MS/MS KW - validation KW - serum KW - quantification KW - apixaban KW - diagnostic accuracy KW - performance liquid chromatography KW - factor XA inhibitor KW - direct oral anticoagulants KW - direct thrombin inhibitor Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-136023 VL - 10 IS - 12 ER - TY - JOUR A1 - Klughammer, Johanna A1 - Dittrich, Marcus A1 - Blom, Jochen A1 - Mitesser, Vera A1 - Vogel, Ulrich A1 - Frosch, Matthias A1 - Goesmann, Alexander A1 - Müller, Tobias A1 - Schoen, Christoph T1 - Comparative genome sequencing reveals within-host genetic changes in Neisseria meningitidis during invasive disease JF - PLoS ONE N2 - Some members of the physiological human microbiome occasionally cause life-threatening disease even in immunocompetent individuals. A prime example of such a commensal pathogen is Neisseria meningitidis, which normally resides in the human nasopharynx but is also a leading cause of sepsis and epidemic meningitis. Using N. meningitidis as model organism, we tested the hypothesis that virulence of commensal pathogens is a consequence of within host evolution and selection of invasive variants due to mutations at contingency genes, a mechanism called phase variation. In line with the hypothesis that phase variation evolved as an adaptation to colonize diverse hosts, computational comparisons of all 27 to date completely sequenced and annotated meningococcal genomes retrieved from public databases showed that contingency genes are indeed enriched for genes involved in host interactions. To assess within-host genetic changes in meningococci, we further used ultra-deep whole-genome sequencing of throat-blood strain pairs isolated from four patients suffering from invasive meningococcal disease. We detected up to three mutations per strain pair, affecting predominantly contingency genes involved in type IV pilus biogenesis. However, there was not a single (set) of mutation(s) that could invariably be found in all four pairs of strains. Phenotypic assays further showed that these genetic changes were generally not associated with increased serum resistance, higher fitness in human blood ex vivo or differences in the interaction with human epithelial and endothelial cells in vitro. In conclusion, we hypothesize that virulence of meningococci results from accidental emergence of invasive variants during carriage and without within host evolution of invasive phenotypes during disease progression in vivo. KW - blood KW - comparative genomics KW - throat KW - genetic loci KW - Neisseria meningitidis KW - genomic libraries KW - genome sequencing KW - sequence assembly tools Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-159547 VL - 12 IS - 1 ER -