@article{WangorschButtMarketal.2011,
  author    = {Wangorsch, Gaby and Butt, Elke and Mark, Regina and Hubertus, Katharina and Geiger, J{\"o}rg and Dandekar, Thomas and Dittrich, Marcus},
  title     = {Time-resolved in silico modeling of fine-tuned cAMP signaling in platelets: feedback loops, titrated phosphorylations and pharmacological modulation},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-69145},
  year      = {2011},
  abstract  = {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.},
  subject      = {Vasodilatator-stimuliertes Phosphoprotein},
  language  = {en}
}
@article{UrbanRemmeleDittrichetal.2020,
  author    = {Urban, Lara and Remmele, Christian W. and Dittrich, Marcus and Schwarz, Roland F. and M{\"u}ller, Tobias},
  title     = {covRNA: discovering covariate associations in large-scale gene expression data},
  series = {BMC Reserach Notes},
  volume    = {13},
  journal   = {BMC Reserach Notes},
  doi       = {10.1186/s13104-020-04946-1},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-229258},
  year      = {2020},
  abstract  = {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.},
  language  = {en}
}
@article{StaigerCadotKooteretal.2012,
  author    = {Staiger, Christine and Cadot, Sidney and Kooter, Raul and Dittrich, Marcus and M{\"u}ller, Tobias and Klau, Gunnar W. and Wessels, Lodewyk F. A.},
  title     = {A Critical Evaluation of Network and Pathway-Based Classifiers for Outcome Prediction in Breast Cancer},
  series = {PLoS One},
  volume    = {7},
  journal   = {PLoS One},
  number    = {4},
  doi       = {10.1371/journal.pone.0034796},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-131323},
  pages     = {e34796},
  year      = {2012},
  abstract  = {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.},
  language  = {en}
}
@article{SchneiderDittrichBoecketal.2016,
  author    = {Schneider, Eberhard and Dittrich, Marcus and B{\"o}ck, Julia and Nanda, Indrajit and M{\"u}ller, Tobias and Seidmann, Larissa and Tralau, Tim and Galetzka, Danuta and El Hajj, Nady and Haaf, Thomas},
  title     = {CpG sites with continuously increasing or decreasing methylation from early to late human fetal brain development},
  series = {Gene},
  volume    = {592},
  journal   = {Gene},
  number    = {1},
  doi       = {10.1016/j.gene.2016.07.058},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-186936},
  pages     = {110-118},
  year      = {2016},
  abstract  = {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.},
  language  = {en}
}
@article{MortonFliesserDittrichetal.2014,
  author    = {Morton, Charles Oliver and Fliesser, Mirjam and Dittrich, Marcus and M{\"u}ller, Tobias and Bauer, Ruth and Kneitz, Susanne and Hope, William and Rogers, Thomas Richard and Einsele, Hermann and L{\"o}ffler, J{\"u}rgen},
  title     = {Gene Expression Profiles of Human Dendritic Cells Interacting with Aspergillus fumigatus in a Bilayer Model of the Alveolar Epithelium/Endothelium Interface},
  doi       = {10.1371/journal.pone.0098279},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-112893},
  year      = {2014},
  abstract  = {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.},
  language  = {en}
}
@article{MaierhoferFlunkertOshimaetal.2019,
  author    = {Maierhofer, Anna and Flunkert, Julia and Oshima, Junko and Martin, George M. and Poot, Martin and Nanda, Indrajit and Dittrich, Marcus and M{\"u}ller, Tobias and Haaf, Thomas},
  title     = {Epigenetic signatures of Werner syndrome occur early in life and are distinct from normal epigenetic aging processes},
  series = {Aging Cell},
  volume    = {18},
  journal   = {Aging Cell},
  doi       = {10.1111/acel.12995},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-202733},
  pages     = {e12995},
  year      = {2019},
  abstract  = {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.},
  language  = {en}
}
@article{MaierhoferFlunkertDittrichetal.2017,
  author    = {Maierhofer, Anna and Flunkert, Julia and Dittrich, Marcus and M{\"u}ller, Tobias and Schindler, Detlev and Nanda, Indrajit and Haaf, Thomas},
  title     = {Analysis of global DNA methylation changes in primary human fibroblasts in the early phase following X-ray irradiation},
  series = {PLoS ONE},
  volume    = {12},
  journal   = {PLoS ONE},
  number    = {5},
  doi       = {10.1371/journal.pone.0177442},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-170895},
  pages     = {e0177442},
  year      = {2017},
  abstract  = {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.},
  language  = {en}
}
@article{LutherBrandtVylkovaetal.2023,
  author    = {Luther, Christian H. and Brandt, Philipp and Vylkova, Slavena and Dandekar, Thomas and M{\"u}ller, Tobias and Dittrich, Marcus},
  title     = {Integrated analysis of SR-like protein kinases Sky1 and Sky2 links signaling networks with transcriptional regulation in Candida albicans},
  series = {Frontiers in Cellular and Infection Microbiology},
  volume    = {13},
  journal   = {Frontiers in Cellular and Infection Microbiology},
  issn      = {2235-2988},
  doi       = {10.3389/fcimb.2023.1108235},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-311771},
  year      = {2023},
  abstract  = {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.},
  language  = {en}
}
@article{KuhnGrippFliederetal.2015,
  author    = {Kuhn, Joachim and Gripp, Tatjana and Flieder, Tobias and Dittrich, Marcus and Hendig, Doris and Busse, Jessica and Knabbe, Cornelius and Birschmann, Ingvild},
  title     = {UPLC-MRM Mass Spectrometry Method for Measurement of the Coagulation Inhibitors Dabigatran and Rivaroxaban in Human Plasma and Its Comparison with Functional Assays},
  series = {PLOS ONE},
  volume    = {10},
  journal   = {PLOS ONE},
  number    = {12},
  doi       = {10.1371/journal.pone.0145478},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-136023},
  pages     = {e0145478},
  year      = {2015},
  abstract  = {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.},
  language  = {en}
}
@article{KlughammerDittrichBlometal.2017,
  author    = {Klughammer, Johanna and Dittrich, Marcus and Blom, Jochen and Mitesser, Vera and Vogel, Ulrich and Frosch, Matthias and Goesmann, Alexander and M{\"u}ller, Tobias and Schoen, Christoph},
  title     = {Comparative genome sequencing reveals within-host genetic changes in Neisseria meningitidis during invasive disease},
  series = {PLoS ONE},
  volume    = {12},
  journal   = {PLoS ONE},
  number    = {1},
  doi       = {10.1371/journal.pone.0169892},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-159547},
  pages     = {e0169892},
  year      = {2017},
  abstract  = {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.},
  language  = {en}
}