1
|
Ascensión AM, Arrospide-Elgarresta M, Izeta A, Araúzo-Bravo MJ. NaviSE: superenhancer navigator integrating epigenomics signal algebra. BMC Bioinformatics 2017; 18:296. [PMID: 28587674 PMCID: PMC5461685 DOI: 10.1186/s12859-017-1698-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2016] [Accepted: 05/18/2017] [Indexed: 01/23/2023] Open
Abstract
Background Superenhancers are crucial structural genomic elements determining cell fate, and they are also involved in the determination of several diseases, such as cancer or neurodegeneration. Although there are pipelines which use independent pieces of software to predict the presence of superenhancers from genome-wide chromatin marks or DNA-interaction protein binding sites, there is not yet an integrated software tool that processes automatically algebra combinations of raw data sequencing into a comprehensive final annotated report of predicted superenhancers. Results We have developed NaviSE, a user-friendly streamlined tool which performs a fully-automated parallel processing of genome-wide epigenomics data from sequencing files into a final report, built with a comprehensive set of annotated files that are navigated through a graphic user interface dynamically generated by NaviSE. NaviSE also implements an ‘epigenomics signal algebra’ that allows the combination of multiple activation and repression epigenomics signals. NaviSE provides an interactive chromosomal landscaping of the locations of superenhancers, which can be navigated to obtain annotated information about superenhancer signal profile, associated genes, gene ontology enrichment analysis, motifs of transcription factor binding sites enriched in superenhancers, graphs of the metrics evaluating the superenhancers quality, protein-protein interaction networks and enriched metabolic pathways among other features. We have parallelised the most time-consuming tasks achieving a reduction up to 30% for a 15 CPUs machine. We have optimized the default parameters of NaviSE to facilitate its use. NaviSE allows different entry levels of data processing, from sra-fastq files to bed files; and unifies the processing of multiple replicates. NaviSE outperforms the more time-consuming processes required in a non-integrated pipeline. Alongside its high performance, NaviSE is able to provide biological insights, predicting cell type specific markers, such as SOX2 and ZIC3 in embryonic stem cells, CDK5R1 and REST in neurons and CD86 and TLR2 in monocytes. Conclusions NaviSE is a user-friendly streamlined solution for superenhancer analysis, annotation and navigation, requiring only basic computer and next generation sequencing knowledge. NaviSE binaries and documentation are available at: https://sourceforge.net/projects/navise-superenhancer/. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1698-5) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Alex M Ascensión
- Computational Biology and Systems Biomedicine, Biodonostia Health Research Institute, San Sebastián, 20014, Spain.,Tissue Engineering Laboratory, Bioengineering Area, Biodonostia Health Research Institute, San Sebastián, 20014, Spain.,Department of Biochemistry and Molecular Biology, University of the Basque Country, Leioa, 48940, Spain
| | - Mikel Arrospide-Elgarresta
- Computational Biology and Systems Biomedicine, Biodonostia Health Research Institute, San Sebastián, 20014, Spain
| | - Ander Izeta
- Tissue Engineering Laboratory, Bioengineering Area, Biodonostia Health Research Institute, San Sebastián, 20014, Spain.
| | - Marcos J Araúzo-Bravo
- Computational Biology and Systems Biomedicine, Biodonostia Health Research Institute, San Sebastián, 20014, Spain. .,IKERBASQUE, Basque Foundation for Science, Bilbao, 48013, Spain.
| |
Collapse
|
2
|
Wynants M, Vengethasamy L, Ronisz A, Meyns B, Delcroix M, Quarck R. NF-κB pathway is involved in CRP-induced effects on pulmonary arterial endothelial cells in chronic thromboembolic pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol 2013; 305:L934-42. [DOI: 10.1152/ajplung.00034.2013] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Chronic thromboembolic pulmonary hypertension (CTEPH) is characterized by thrombofibrotic obstruction of proximal pulmonary arteries. The cellular and molecular mechanisms underlying the pathogenesis remain incompletely understood, although we recently evidenced the potential involvement of the inflammatory marker C-reactive protein (CRP). We aimed to investigate the intracellular mechanisms induced by CRP in proximal pulmonary arterial endothelial cells (PAEC). PAEC were isolated from vascular material obtained during pulmonary endarterectomy. RNA was extracted from CRP-stimulated PAEC, and first-stand cDNA was generated. A RT2 profiler PCR Array was used to evaluate the expression of 84 key genes related to NF-κB-mediated signal transduction. CRP-induced NF-κB activation was studied. The effects of pyrrolidine-dithio-carbamate ammonium (PDTC), an inhibitor of the NF-κB pathway, were investigated on CRP-induced adhesion of monocytes to PAEC, adhesion molecule expression, endothelin-1 (ET-1), interleukin-6 (IL-6), and von Willebrand factor (vWF) secretion. Compared with nonstimulated PAEC, serotonin receptor 2B was downregulated by 25%, inhibitor of NF-κB kinase subunit epsilon (IKBKE) by 30%, and toll-like receptor-4 and -6 by 18 and 39%, respectively, in CRP-stimulated PAEC. The transcription factor FOS was threefold upregulated. CRP induced RelA/NF-κBp65 phosphorylation. PDTC dose dependently inhibited the adhesion of monocytes to CRP-stimulated PAEC. PDTC also inhibited the CRP-induced expression of ICAM-1 at the surface of PAEC. PDTC impaired the secretion of ET-1 by 18% and tended to inhibit the secretion of IL-6 by CRP-stimulated PAEC by 46%. PDTC did not inhibit the CRP-induced secretion of vWF. These results suggest an involvement of the NF-κB pathway in mediating different effects of CRP on proximal CTEPH-PAEC.
Collapse
Affiliation(s)
- Marijke Wynants
- Respiratory Division, Department of Clinical and Experimental Medicine, KU Leuven, Leuven, Belgium
| | - Leanda Vengethasamy
- Respiratory Division, Department of Clinical and Experimental Medicine, KU Leuven, Leuven, Belgium
| | - Alicja Ronisz
- Respiratory Division, Department of Clinical and Experimental Medicine, KU Leuven, Leuven, Belgium
| | - Bart Meyns
- Cardiac Surgery Department, University Hospitals of Leuven, Leuven, Belgium; and
| | - Marion Delcroix
- Respiratory Division, University Hospitals and Department of Clinical and Experimental Medicine, University of Leuven, Leuven, Belgium
| | - Rozenn Quarck
- Respiratory Division, Department of Clinical and Experimental Medicine, KU Leuven, Leuven, Belgium
| |
Collapse
|
3
|
Maiwald S, Zwetsloot PP, Sivapalaratnam S, Dallinga-Thie GM. Monocyte gene expression and coronary artery disease. Curr Opin Clin Nutr Metab Care 2013; 16:411-7. [PMID: 23739627 DOI: 10.1097/mco.0b013e32836236f9] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
PURPOSE OF REVIEW Despite current therapy, coronary artery disease (CAD) remains the major cause of morbidity and mortality worldwide. CAD is the consequence of a complex array of deranged metabolic processes including the immune system. In this context, monocytes and macrophages are indisputable players. Thus, monocyte gene expression analysis could be a powerful tool to provide new insights in the pathophysiology of CAD and improve identification of individuals at risk. We discuss current literature assessing monocyte gene expression and its association with CAD. RECENT FINDINGS Monocyte surface markers CD14 ⁺⁺and CD16⁺ have been established as biomarkers for increased cardiovascular disease risk in a large number of studies. More in-depth gene expression analysis identified several interesting genes, such as ABCA1, CD36 and MSR1 with an increased expression in circulating monocytes from patients with CAD. The results for CD36 were replicated in one other study. For ABCA1 and MSR1 conflicting data are published. SUMMARY Recent findings indicate that genetic differences exist in circulating monocytes of patients suffering from CAD, giving us more insights into the underlying mechanisms. However, larger studies are required to prove that monocytes' expression signature could serve as a marker for diagnostic purposes in the future.
Collapse
Affiliation(s)
- Stephanie Maiwald
- Department of Vascular Medicine, Academic Medical Center, Amsterdam, The Netherlands
| | | | | | | |
Collapse
|
4
|
Masud R, Shameer K, Dhar A, Ding K, Kullo IJ. Gene expression profiling of peripheral blood mononuclear cells in the setting of peripheral arterial disease. J Clin Bioinforma 2012; 2:6. [PMID: 22409835 PMCID: PMC3381689 DOI: 10.1186/2043-9113-2-6] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2011] [Accepted: 03/12/2012] [Indexed: 12/26/2022] Open
Abstract
BACKGROUND Peripheral arterial disease (PAD) is a relatively common manifestation of systemic atherosclerosis that leads to progressive narrowing of the lumen of leg arteries. Circulating monocytes are in contact with the arterial wall and can serve as reporters of vascular pathology in the setting of PAD. We performed gene expression analysis of peripheral blood mononuclear cells (PBMC) in patients with PAD and controls without PAD to identify differentially regulated genes. METHODS PAD was defined as an ankle brachial index (ABI) ≤0.9 (n = 19) while age and gender matched controls had an ABI > 1.0 (n = 18). Microarray analysis was performed using Affymetrix HG-U133 plus 2.0 gene chips and analyzed using GeneSpring GX 11.0. Gene expression data was normalized using Robust Multichip Analysis (RMA) normalization method, differential expression was defined as a fold change ≥1.5, followed by unpaired Mann-Whitney test (P < 0.05) and correction for multiple testing by Benjamini and Hochberg False Discovery Rate. Meta-analysis of differentially expressed genes was performed using an integrated bioinformatics pipeline with tools for enrichment analysis using Gene Ontology (GO) terms, pathway analysis using Kyoto Encyclopedia of Genes and Genomes (KEGG), molecular event enrichment using Reactome annotations and network analysis using Ingenuity Pathway Analysis suite. Extensive biocuration was also performed to understand the functional context of genes. RESULTS We identified 87 genes differentially expressed in the setting of PAD; 40 genes were upregulated and 47 genes were downregulated. We employed an integrated bioinformatics pipeline coupled with literature curation to characterize the functional coherence of differentially regulated genes. CONCLUSION Notably, upregulated genes mediate immune response, inflammation, apoptosis, stress response, phosphorylation, hemostasis, platelet activation and platelet aggregation. Downregulated genes included several genes from the zinc finger family that are involved in transcriptional regulation. These results provide insights into molecular mechanisms relevant to the pathophysiology of PAD.
Collapse
Affiliation(s)
- Rizwan Masud
- Division of Cardiovascular Diseases, Mayo Clinic, Rochester MN 55905, USA
| | - Khader Shameer
- Division of Cardiovascular Diseases, Mayo Clinic, Rochester MN 55905, USA
| | - Aparna Dhar
- Division of Cardiovascular Diseases, Mayo Clinic, Rochester MN 55905, USA
| | - Keyue Ding
- Division of Cardiovascular Diseases, Mayo Clinic, Rochester MN 55905, USA
| | - Iftikhar J Kullo
- Division of Cardiovascular Diseases, Mayo Clinic, Rochester MN 55905, USA
| |
Collapse
|
5
|
Differentially expressed genes in human peripheral blood as potential markers for statin response. J Mol Med (Berl) 2011; 90:201-11. [PMID: 21947165 DOI: 10.1007/s00109-011-0818-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2011] [Revised: 08/25/2011] [Accepted: 09/12/2011] [Indexed: 10/17/2022]
Abstract
There is a considerable inter-individual variation in response to statin therapy and one third of patients do not meet their treatment goals. We aimed to identify differentially expressed genes that might be involved in the effects of statin treatment and to suggest potential markers to guide statin therapy. Forty-six healthy Korean subjects received atorvastatin; their whole-genome expression profiles in peripheral blood were analyzed before and after atorvastatin administration in relation with changes in lipid profiles. The expression patterns of the differentially expressed genes were also compared with the data of familial hypercholesterolemia (FH) patients and controls. Pairwise comparison analyses revealed differentially expressed genes involved in diverse biological processes and molecular functions related with immune responses. Atorvastain mainly affected antigen binding, immune or inflammatory response including interleukin pathways. Similar expression patterns of the genes were observed in patients with FH and controls. The Charcol-Leyden crystal (CLC), CCR2, CX3CR1, LRRN3, FOS, LDLR, HLA-DRB1, ERMN, and TCN1 genes were significantly associated with cholesterol levels or statin response. Interestingly, the CLC gene, which was significantly altered by atorvastatin administration and differentially expressed between FH patients and controls, showed much bigger change in high-responsive group than in low-responsive group. We identified differentially expressed genes that might be involved in mechanisms underlying the known pleiotropic effects of atorvastatin, baseline cholesterol levels, and drug response. Our findings suggest CLC as a new candidate marker for statin response, and further validation is needed.
Collapse
|
6
|
Saliques S, Teyssier JR, Vergely C, Lorgis L, Lorin J, Donzel A, Sicard P, Berchoud J, Ragot S, Touzery C, Cottin Y, Rochette L, Zeller M. Smoking and FOS expression from blood leukocyte transcripts in patients with coronary artery disease. Atherosclerosis 2011; 219:931-6. [PMID: 21978922 DOI: 10.1016/j.atherosclerosis.2011.09.026] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/15/2011] [Revised: 08/12/2011] [Accepted: 09/11/2011] [Indexed: 10/17/2022]
Abstract
OBJECTIVE Analysis of the leukocyte transriptome, in particular the Finkel-Biskis-Jinkins Osteosarcoma (c-Fos) gene, which has a prominent role in inflammation, provides new insights into atherosclerosis mechanisms. Although smoking is a major risk factor, the links between smoking status and coronary artery disease (CAD) remains unclear. We aimed to analyze the relationship between smoking status and c-Fos expression in circulating leukocytes of patients with CAD. METHODS c-Fos expression was measured by RT-Q-PCR, from blood leukocytes of 239 consecutive patients after acute myocardial infarction (MI). The patients were asked about their smoking status and stratified into 3 groups: current smokers (CS) (N = 85), past smokers (PS) (N = 78) and never smokers (NS) (N = 76). RESULTS NS had a higher risk profile including hypertension, and CS were younger than PS and NS (-13 years and 17 years respectively). There was only a trend towards lower CRP levels in NS and PS than in CS. The mean c-Fos transcript level was slightly higher in CS than in PS and NS (0.924 vs. 0.908 and 0.861 AU, respectively; p = 0.005). By univariate analysis, neither age, nor sex, nor CRP nor white blood cell count was associated with c-Fos transcript levels. By multivariate analysis, CS (vs. PS + NS) was the strongest predictor of the c-Fos transcript level, (B = 0.042 ± 0.014, p = 0.003), even after adjustment for confounding factors (i.e. hypertension, chronic medication, family history of CAD, and prior MI). CONCLUSION Our work suggests that c-Fos transcript level in blood leukocyte could be considered a cumulative biomarker of smoking. As the c-Fos gene has been put forward as a new factor in the progression and severity of atherosclerosis, it could be considered a novel potential pathway of tobacco toxicity in coronary artery disease.
Collapse
Affiliation(s)
- Sebastien Saliques
- IFR 100 Santé-STIC, Laboratoire de Physiopathologie et Pharmacologie Cardiovasculaires Expérimentales, UFR Médecine, Université de Bourgogne EA2979, Dijon, France
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|