|Year : 2019 | Volume
| Issue : 3 | Page : 218-222
A bioinformatics analysis of exosomal microRNAs released following mycobacterial infection
Shamila D Alipoor1, Ian M Adcock2, Gert Folkerts3, Johan Garssen4, Esmaeil Mortaz5
1 Department of Molecular Medicine, Institute of Medical Biotechnology, National Institute of Genetic Engineering and Biotechnology, Tehran, Iran
2 Priority Research Centre for Asthma and Respiratory Disease, Hunter Medical Research Institute, University of Newcastle, Newcastle, NSW, Australia; Airways Section, National Heart and Lung Institute, Imperial College London, London, United Kingdom
3 Division of Pharmacology, Faculty of Science, Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Utrecht, Netherlands
4 Division of Pharmacology, Faculty of Science, Utrecht Institute for Pharmaceutical Sciences, Utrecht University; Immunology Section, Nutricia Research, Utrecht, Netherlands
5 Division of Pharmacology, Faculty of Science, Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Utrecht, Netherlands; Clinical Tuberculosis and Epidemiology Research Center, National Research Institute of Tuberculosis and Lung Diseases, Shahid Beheshti University of Medical Sciences, Tehran, Iran
|Date of Web Publication||12-Sep-2019|
Dr Esmaeil Mortaz
Clinical Tuberculosis and Epidemiology Research Center, National Research Institute of Tuberculosis and Lung Diseases, Shahid Beheshti University of Medical Sciences, Tehran; Division of Pharmacology, Faculty of Science, Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Utrecht
Source of Support: None, Conflict of Interest: None
Background: Tuberculosis (TB) still remains a major health threat worldwide. The current TB diagnostics are suboptimal, and there is a high clinical need for identifying novel biomarkers of disease prevalence. Circulating exosomes have been currently attractive as novel biomarkers in a wide range of pathological conditions. Methods: In this study, we performed bioinformatics analysis on the downstream targets of a dysregulated microRNA (miRNA) cluster induced by Bacillus Calmette–Guerin infection of human macrophages to provide greater understanding of their potential roles in disease pathogenesis. Results: Our analysis demonstrated that these dysregulated miRNAs have central roles in the host metabolic and energy pathways. Conclusion: This suggests that the host miRNA network is perturbed by Mycobacterium to re-patterning host metabolism machinery to favor its intracellular survival. The dysregulated miRNAs can be delivered to local and distal cells by exosomes and thereby modulate their function.
Keywords: Clusters, exosomal, microRNAs, tuberculosis
|How to cite this article:|
Alipoor SD, Adcock IM, Folkerts G, Garssen J, Mortaz E. A bioinformatics analysis of exosomal microRNAs released following mycobacterial infection. Int J Mycobacteriol 2019;8:218-22
|How to cite this URL:|
Alipoor SD, Adcock IM, Folkerts G, Garssen J, Mortaz E. A bioinformatics analysis of exosomal microRNAs released following mycobacterial infection. Int J Mycobacteriol [serial online] 2019 [cited 2019 Nov 19];8:218-22. Available from: http://www.ijmyco.org/text.asp?2019/8/3/218/266496
| Introduction|| |
Tuberculosis (TB) is a potentially serious infection that still remains the most common cause of death from infection, despite the global advances in health and disease., According to a recent World Health Organization report, 9.6 million new TB cases and 1.5 million TB deaths occur each year., In addition, about one-third of the world's population has latent TB showing no overt symptoms of disease but acting as a source for spreading the infection. Since the early detection of TB is important in controlling and preventing infections from spreading, and because of the challenges with current TB diagnostics, the development of novel biomarkers is essential., Exosomes have been proposed as novel diagnostic biomarkers in a wide range of pathological conditions, such as cancers and infectious diseases.,,,
Exosomes are 30–100 nm vesicles containing functional molecules, such as nucleic acids (mRNA and non-coding RNAs, DNA), proteins, metabolites, and lipid mediators., Exosomes can be derived from nearly all cell types and are isolated from almost all human biofluids., Exosomes play a key role in cell–cell communication and transfer biological information by shuttling their cargo to either local or more distant cells and thereby modulate the function of the recipient cells.,, Exosomal contents have been used as disease signatures in various cancers ,, as well as in Alzheimer's disease.
MicroRNAs (miRNAs) are small 18–22 ntRNAs that have an important role in modulating gene expression and translation. miRNAs can affect most biological functions, and their dysregulation is associated with several pathologies. Functional miRNAs may be encapsulated inside exosomes, delivered to target cells causing effects on the recipient cell function by affecting their transcriptome and/or proteome.
miRNAs are involved in the regulation of inflammatory processes during Mycobacterium tuberculosis (Mtb) infection., Mtb infection triggers a range of physiological responses in infected cells, leading to host immune dysregulation and metabolic re-patterning., These modulations in host cell functions enable bacteria to sequester vital host factors to supply their food and energy requirements to enable intracellular survival. These processes may be controlled by subversion of host miRNA networks that are involved in the regulation of carbon, nitrogen, and lipid metabolism in the infected host cells.,
We have previously observed that infection of human monocyte-derived macrophages (MDMs) with Mycobacterium bovis Bacillus Calmette–Guerin (BCG) induced the secretion of a specific set of exosomal miRNAs, including mir-1224, mir-484, mir-1293, mir-423, and miR-96, which are involved in modulating key host metabolic and energy production pathways as well as regulating immunological and cell signaling events. In this study, we hypothesized that further bioinformatic analysis of the downstream targets of dysregulated miRNA cluster induced by BCG infection may provide further insight into their potential roles.
| Methods|| |
MicroRNAs network and enrichment analysis
The isolation and detection of exosomal miRNAs from BCG-exposed human MDMs have been described previously.
The target genes of the interested miRNAs were determined using miRTarBase  and microT-CDS algorithms. Protein–protein interaction (PPI) networks were constructed and analyzed by Cytoscape v. 3.4.0. Briefly, miRNA target genes were converted to seed proteins and PPI networks were constructed using GeneMANIA plugin and visualized and analyzed in Cytoscape.
GeneMANIA provides information on both the experimental and the predicted interactions based on their gene neighborhood, gene co-occurrence, gene fusions, gene co-expression, experiments, and literature mining.
The highly interconnected regions (clusters) were extracted using the Molecular COmplex DEtection (MCODE) plugin within Cytoscape. This clustering method isolates the dense proteins based on their local neighborhood density and distance from a local seed protein. Enrichment analysis was performed using EnrichR  and KEGG pathway  database for nodes within each cluster.
| Results|| |
Protein–protein interaction network of the candidate microRNA target genes and enrichment analysis
Differentially regulated miRNAs were called according to a log2-fold difference of greater than 2 and/or an adjusted P < 0.05. This gave a total of five miRNAs – miR-1224, miR-484, miR-1293, miR-423, and miR-96. To survey the interactions between the genes targeted by the identified dysregulated miRNAs, a PPI network was constructed. We first constructed an extended network of the seed proteins, their direct PPI neighbors, and their interactions. The extended network consisted of 96 nodes connected via 1022 edges. Extracting the highly connected regions from this network using MCODE identified five clusters [Figure 1]. The number of nodes in clusters 1–5 was 16, 12, 11, 10, and 15 with 90, 45, 27, 20, and 28 edges, respectively. The score of these subnetworks was 1 = 7.86; 2 = 5.81; 3 = 4.4; 4 = 4.22; and 5 = 3.57 [Table 1]. To understand the biological meaning behind the highly connected regions, enrichment analysis was performed on the nodes within each cluster. Using the EnrichR tool, the most significant gene ontology biological process terms identified within the clusters are provided in [Table 2]. The results showed that cluster 1 with the highest score of 7.8 contains proteins that participate in B-cell activation involved in immune response; cluster 2 was most associated with positive regulation of phosphatidylinositol 3-kinase signaling; cluster 3 with positive regulation of lipid metabolic process; cluster 4 with inositol phosphate biosynthetic process; and cluster 5 with the Fc receptor signaling pathway.
|Figure 1: Protein–protein interaction network of the dysregulated microRNA gene targets. Protein–protein interaction networks were constructed and the clusters contained within the highly connected regions in the network identified by Molecular COmplex DEtection plugin. (a) Cluster 1, (b) cluster 2, (c) cluster 3, (d) cluster 4, and (e) cluster 5. Nodes represent proteins and edges represent interactions between two proteins. Node size indicates the generality of each gene ontology term (larger = more specific). Node color indicates the P values of each gene ontology term in each intersection (darker = more significant)|
Click here to view
|Table 1: Protein–protein interaction data of intersection network identified by Molecular COmplex DEtection plugin|
Click here to view
|Table 2: Gene ontology biological process identified in the Molecular COmplex DEtection clusters using the EnrichR tool|
Click here to view
Pathway analysis utilizing KEGG, Reactome, and WikiPathways within EnrichR revealed that the target genes within each cluster were mostly involved in signaling pathways, metabolism, and immunological pathways [Table 3]. In particular, cluster 1 was most associated with endocytosis, bacterial invasion of epithelial cells, and EGF/EGFR signaling; cluster 2 with adherents junction, focal adhesion, and regulation of lipid metabolism; cluster 3 with N-glycan biosynthesis, proteasome degradation, and metabolic pathways; cluster 4 with the adaptive immune system, antigen processing, and ubiquitination and proteasome degradation; and cluster 5 with calcium signaling, focal adhesion, and B-cell receptor and interleukin-3 signaling.
|Table 3: Molecular COmplex DEtection clusters and their related pathways|
Click here to view
| Discussion|| |
In the current study, we assessed the potential regulatory effect of a cluster of dysregulated exosomal miRNAs released by human MDMs upon infection with Mycobacterium. Bioinformatic analysis demonstrated that these miRNAs regulate key metabolic and energy pathways including central carbon metabolism, fatty acids and sugar metabolism, bacterial invasion, and immunological-related pathways as well as cell signaling pathways.
This finding suggests that Mycobacterium may recruit host miRNAs to modulate host metabolic reprograming or re-patterning to favor its intracellular survival. These dysregulated miRNAs can be enclosed inside exosomes and released into the biological fluids to modulate the function of local and distal cells. Exosomes have been reported as novel biomarkers in several diseases, including acute myeloid leukemia, asthma, and sarcoidosis., In TB, the exosomal miRNA profile has recently been proposed as a diagnostic biomarker.
miRNAs generally exert their effects in the context of complex regulatory networks. Understanding the role of miRNAs in the modulation of pathological conditions depends on our knowledge of their functions in the context of these networks. Network analysis of miRNA incorporates known and predicted protein–protein interaction data, based on their gene neighborhood, gene fusions, gene co-occurrence, gene co-expression interactions determined from both laboratory and text mining analysis., In a given network, nodes are genes and the interactions between the nodes are defined as edges. Node degree shows the number of edges connected to a node, and nodes with high degree may indicate the key proteins with important biological functions.
Our miRNA network gave a total of 1022 interaction pairs over 96 nodes. Cluster 1 was the most connected with the highest score and as with two or the five other clusters contained nodes enriched for metabolism and infection-related pathways. This suggests that these miRNAs and their gene targets have extensive interactions that control the host response to Mycobacterium infection. Specific proteins targeted by this cluster include mitochondrial fission protein 1 [Figure 1] that is regulated by miR-484 and a number of lipid metabolizing proteins targeted by miR-1224. Furthermore, miR-425 has also been implicated in the control of various metabolic pathways and is associated with several metabolic disorders.
In addition, both miR-425 and miR-96 mediate insulin resistance with miR-96 regulating the expression of multiple genes that fine-tune insulin release., For example, miR-96 suppresses insulin secretion by increasing the level of granuphilin, an inhibitor of insulin exocytosis. miR-96 is strongly induced by fatty acids and downregulates the expression of the insulin receptor (INSR) and the INSR substrate 1. On the other hand, Mtb-infected macrophages show a disturbed sugar flux and a significant change in the intracellular levels of glucose, glycogen, NAD/NADP, and lactate., Mtb infection also causes increased aerobic glycolysis as a result of the pentose-phosphate shunt and elevated glucose uptake.,
In this small pilot study, we have obtained a comprehensive network analysis of the pathways affected by the dysregulated miRNAs found in exosomes from Mtb-infected human macrophages. We demonstrated that this cluster of miRNAs modulate not only key host immune functions but also re-pattern critical metabolic and energy production pathways. Since the nutrients required for intracellular survival of Mtb are restricted inside the host cell, it is likely that the bacteria sequester the host's cellular metabolism to optimize the chances of survival and growth. Our data suggest that this is accomplished, at least in part, through dysregulation of host miRNA networks that control the host's carbon and nitrogen metabolism. These Mtb-induced dysregulated miRNAs also impact on amino acid, pyrimidine, and purine nucleotide biosynthesis which are increased in Mtb-infected lung tissue.
Finally, the target genes of the altered miRNAs also were involved in cell membrane and communication pathways, such as cell junction, glycosaminoglycan biosynthesis, and heparan sulfate metabolism. This further emphasizes the ability of Mtb to alter the host cell membrane structure to support its intracellular growth and proliferation.
| Conclusion|| |
These findings emphasize the multifactorial effects of miRNAs on cellular metabolism and shed light on the host–pathogen interaction and modulation in cellular function following Mtb infection. The data highlight key PPI networks that have a potential pivotal role in host metabolic reprogramming to favor Mtb survival. Functional analysis of these effects is required to confirm their roles. It also may be possible to use exosomal dysregulated miRNAs as potential biomarkers for TB therapy monitoring and diagnosis of active TB.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Jia T, Fu H, Sun J, Zhang Y, Yang W, Li Y. Foxp3 expression in A549 cells is regulated by toll-like receptor 4 through nuclear factor-κB. Mol Med Rep 2012;6:167-72.
Rajaei E, Hadadi M, Madadi M, Aghajani J, Ahmad MM, Farnia P, et al
. Outdoor air pollution affects tuberculosis development based on geographical information system modeling. Biomed Biotechnol Res J 2018;2:39. [Full text]
Chockalingam AK, Hamed S, Goodwin DG, Rosenzweig BA, Pang E, Boyne MT 2nd
, et al.
The effect of oseltamivir on the disease progression of lethal influenza A virus infection: Plasma cytokine and miRNA responses in a mouse model. Dis Markers 2016;2016:9296457.
Kruh-Garcia NA, Wolfe LM, Dobos KM. Deciphering the role of exosomes in tuberculosis. Tuberculosis (Edinb) 2015;95:26-30.
Kiani A, Razavi F, Mortaz E, Emami H, Ghazali S, Anbardan AD, et al
. The prevalence of latent tuberculosis infection among Iranian sarcoidosis patients. Biomed Biotechnol Res J 2018;2:247. [Full text]
Velayati AA, Abeel T, Shea T, Konstantinovich Zhavnerko G, Birren B, Cassell GH, et al.
Populations of latent Mycobacterium tuberculosis
lack a cell wall: Isolation, visualization, and whole-genome characterization. Int J Mycobacteriol 2016;5:66-73. [Full text]
Velayati AA, Farnia P, Mozafari M, Mirsaeidi M. Nontuberculous mycobacteria isolation from clinical and environmental samples in Iran: Twenty years of surveillance. Biomed Res Int 2015;2015:254285.
Tsara V, Serasli E, Christaki P. Problems in diagnosis and treatment of tuberculosis infection. Hippokratia 2009;13:20-2.
Schorey JS, Harding CV. Extracellular vesicles and infectious diseases: New complexity to an old story. J Clin Invest 2016;126:1181-9.
Schorey JS, Dobos KM. Exosomes and Diagnostic Biomarkers. Google Patents; 2016.
Raposo G, Stoorvogel W. Extracellular vesicles: Exosomes, microvesicles, and friends. J Cell Biol 2013;200:373-83.
Frydrychowicz M, Kolecka-Bednarczyk A, Madejczyk M, Yasar S, Dworacki G. Exosomes – Structure, biogenesis and biological role in non-small-cell lung cancer. Scand J Immunol 2015;81:2-10.
Lin J, Li J, Huang B, Liu J, Chen X, Chen XM, et al.
Exosomes: Novel biomarkers for clinical diagnosis. ScientificWorldJournal 2015;2015:657086.
Kourembanas S. Exosomes: Vehicles of intercellular signaling, biomarkers, and vectors of cell therapy. Annu Rev Physiol 2015;77:13-27.
Zöller M. Exosomes in cancer disease. Methods Mol Biol 2016;1381:111-49.
Taylor DD, Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol 2008;110:13-21.
Nilsson J, Skog J, Nordstrand A, Baranov V, Mincheva-Nilsson L, Breakefield XO, et al.
Prostate cancer-derived urine exosomes: A novel approach to biomarkers for prostate cancer. Br J Cancer 2009;100:1603-7.
Yang TT, Liu CG, Gao SC, Zhang Y, Wang PC. The serum exosome derived microRNA-135a, -193b, and -384 were potential Alzheimer's disease biomarkers. Biomed Environ Sci 2018;31:87-96.
Alipoor SD, Adcock IM, Garssen J, Mortaz E, Varahram M, Mirsaeidi M, et al.
The roles of miRNAs as potential biomarkers in lung diseases. Eur J Pharmacol 2016;791:395-404.
Alipoor SD, Mortaz E, Garssen J, Movassaghi M, Mirsaeidi M, Adcock IM. Exosomes and exosomal miRNA in respiratory diseases. Mediators Inflamm 2016;2016:5628404.
Alipoor SD, Mortaz E, Tabarsi P, Farnia P, Mirsaeidi M, Garssen J, et al.
Calmette-Guerin (BCG) infection induces exosomal miRNA release by human macrophages. J Transl Med 2017;15:105.
Furci L, Schena E, Miotto P, Cirillo DM. Alteration of human macrophages microRNA expression profile upon infection with Mycobacterium tuberculosis
. Int J Mycobacteriol 2013;2:128-34. [Full text]
Moschos SA, Williams AE, Perry MM, Birrell MA, Belvisi MG, Lindsay MA. Expression profiling in vivo
demonstrates rapid changes in lung microRNA levels following lipopolysaccharide-induced inflammation but not in the anti-inflammatory action of glucocorticoids. BMC Genomics 2007;8:240.
Mehrotra P, Jamwal SV, Saquib N, Sinha N, Siddiqui Z, Manivel V, et al.
Pathogenicity of Mycobacterium tuberculosis
is expressed by regulating metabolic thresholds of the host macrophage. PLoS Pathog 2014;10:e1004265.
Eisenreich W, Heesemann J, Rudel T, Goebel W. Metabolic host responses to infection by intracellular bacterial pathogens. Front Cell Infect Microbiol 2013;3:24.
Smith I. Mycobacterium tuberculosis
pathogenesis and molecular determinants of virulence. Clin Microbiol Rev 2003;16:463-96.
Chou CH, Chang NW, Shrestha S, Hsu SD, Lin YL, Lee WH, et al.
MiRTarBase 2016: Updates to the experimentally validated miRNA-target interactions database. Nucleic Acids Res 2016;44:D239-47.
Paraskevopoulou MD, Georgakilas G, Kostoulas N, Vlachos IS, Vergoulis T, Reczko M, et al.
DIANA-microT web server v5.0: Service integration into miRNA functional analysis workflows. Nucleic Acids Res 2013;41:W169-73.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al.
Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res 2003;13:2498-504.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al.
Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res 2016;44:W90-7.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000;28:27-30.
Gurtan AM, Sharp PA. The role of miRNAs in regulating gene expression networks. J Mol Biol 2013;425:3582-600.
Hu JX, Thomas CE, Brunak S. Network biology concepts in complex disease comorbidities. Nat Rev Genet 2016;17:615-29.
Lovis P, Gattesco S, Regazzi R. Regulation of the expression of components of the exocytotic machinery of insulin-secreting cells by microRNAs. Biol Chem 2008;389:305-12.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al.
STRING v10: Protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res 2015;43:D447-52.
Wang K, Long B, Jiao JQ, Wang JX, Liu JP, Li Q, et al.
MiR-484 regulates mitochondrial network through targeting fis1. Nat Commun 2012;3:781.
Vickers KC, Sethupathy P, Baran-Gale J, Remaley AT. Complexity of microRNA function and the role of isomiRs in lipid homeostasis. J Lipid Res 2013;54:1182-91.
Barwari T, Skroblin P, Mayr M. When sweet turns salty: Glucose-induced suppression of atrial natriuretic peptide by microRNA-425. J Am Coll Cardiol 2016;67:813-6.
Chakraborty C, George Priya Doss C, Bandyopadhyay S. MiRNAs in insulin resistance and diabetes-associated pancreatic cancer: The 'minute and miracle' molecule moving as a monitor in the 'genomic galaxy'. Curr Drug Targets 2013;14:1110-7.
Yang WM, Min KH, Lee W. Induction of miR-96 by dietary saturated fatty acids exacerbates hepatic insulin resistance through the suppression of INSR and IRS-1. PLoS One 2016;11:e0169039.
Amila A, Acosta A, Sarmiento ME, Suraiya S, Zafarina Z, Panneerchelvam S, et al.
Sequence comparison of six human microRNAs genes between tuberculosis patients and healthy individuals. Int J Mycobacteriol 2015;4:341-6. [Full text]
Ndzi EN, Indu Viswanath AN, Adzemye NG, Tamgue O, Nsongka MV, Nair AS, et al.
Upregulated bovine tuberculosis microRNAs trigger oncogenic pathways: An In silico
perception. Int J Mycobacteriol 2019;8:70-4.
] [Full text]
Schuhmacher Neto R, Giacomelli IL, Schuller Nin C, da Silva Moreira J, Comaru Pasqualotto A, Marchiori E, et al.
High-resolution CT findings of pulmonary tuberculosis in liver transplant patients. Clin Radiol 2017;72:899.e9-8.99E+16.
Shin JH, Yang JY, Jeon BY, Yoon YJ, Cho SN, Kang YH, et al.
(1) H NMR-based metabolomic profiling in mice infected with Mycobacterium tuberculosis
. J Proteome Res 2011;10:2238-47.
[Table 1], [Table 2], [Table 3]