|Year : 2019 | Volume
| Issue : 1 | Page : 70-74
Upregulated bovine tuberculosis microRNAs Trigger oncogenic pathways: An In silico perception
Elvis Ndukong Ndzi1, Ambily Nath Indu Viswanath2, Nsadzetsen Gilbert Adzemye1, Ousman Tamgue3, Munji Victorine Nsongka1, Achuthsankar S Nair2, Céline Nguefeu Nkenfou4
1 Department of Animal Production and Fisheries, Laboratory of Animal Physiology and Health, Institute of Agricultural Research for Development, Bambui, Cameroon
2 Department of Computational Biology and Bioinformatics, University of Kerala, Thiruvananthapuram, Kerala, India
3 Division of Immunology and South Africa Medical Research Council, Immunology of Infectious Disease, Faculty of Health Sciences, University of Cape Town, Institute of Infectious Diseases and Molecular Medicine, South Africa; Department of Biochemistry, University of Douala, Cameroon
4 System Biology Laboratory, Chantal Biya International Reference Center, Yaounde; Department of Biology, Higher Teachers' Training College, University of Yaounde I, Cameroon
|Date of Web Publication||12-Mar-2019|
Elvis Ndukong Ndzi
Laboratory of Animal Physiology and Health, Institute of Agricultural Research for Development, P. O. Box 51, Irad Bambui
Source of Support: None, Conflict of Interest: None
Background/Objective: Although microRNA (miRNA)-directed regulation of bovine tuberculosis (bTB) has already been reported, very little is known about the incited pathways and genes. We profiled bTB-upregulated miRNAs through an in silico methodology. Methods: The data of upregulated miRNAs in bTB versus healthy controls were collected and clustered into three groups by their tissue specificity as follows: G1 (mammary gland-specific): bta-miR-146a; G2 (peripheral blood mononuclear cell-specific): bta-miR-155; and G3 (alveolar macrophage-specific): bta-miR-146a, bta-miR-155, bta-miR-142-5p, bta-miR-423-3p, bta-miR-21-5p, bta-miR-27a-3p, bta-miR-99b, bta-miR-147, bta-miR-223, and bta-let-7i. The miRNA–mRNA interaction network was defined by TargetScan. The gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways of these transcripts were examined. Results: The results illustrate the induction of pathways in cancer, highly enriched, and unanimous to all three gene sets (G1, G2, and G3). Mitogen-activated protein kinase and PI3K-Akt signaling were specific to G2 and G3 with fibroblast growth factors formed the key factors. Conclusion: The inferred cancer cascades denote a probable modulation of innate immune response in an infectious state. These baseline pictures could lay the ground for further substantive studies.
Keywords: Bovine tuberculosis, computational, microRNA, pathway prediction, target gene
|How to cite this article:|
Ndzi EN, Indu Viswanath AN, Adzemye NG, Tamgue O, Nsongka MV, Nair AS, Nkenfou CN. Upregulated bovine tuberculosis microRNAs Trigger oncogenic pathways: An In silico perception. Int J Mycobacteriol 2019;8:70-4
|How to cite this URL:|
Ndzi EN, Indu Viswanath AN, Adzemye NG, Tamgue O, Nsongka MV, Nair AS, Nkenfou CN. Upregulated bovine tuberculosis microRNAs Trigger oncogenic pathways: An In silico perception. Int J Mycobacteriol [serial online] 2019 [cited 2019 Mar 25];8:70-4. Available from: http://www.ijmyco.org/text.asp?2019/8/1/70/253967
| Introduction|| |
Human (Homo sapiens) and cattle (Bos taurus) microRNA (miRNAs) have been shown to be associated with tuberculosis (TB).,,,,,,,,,, However, little is known about the role played or pathways influenced by these miRNAs in the pathogenicity of TB.
The present study performed the functional profiling of miRNAs relatively upregulated in bovine TB compared to healthy controls using a computational approach. We compared the effector genes and functional pathways of miRNAs specific to the mammary gland, peripheral blood mononuclear cells (PBMCs), and alveolar macrophages induced at higher levels in the infectious state.
| Methods|| |
A systematic biomedical literature hunt in PubMed was made with the keywords “miRNA and bovine TB (bTB).” A parallel search was also conducted in Google with the same queries to calibrate the data mining. In total, 15 articles were tabulated; the retrieved substances were collated and sorted to retain miRNAs overexpressed in bTB as compared to non-TB infected (healthy controls) bovine. The miRNA sequences were downloaded from miRBase release 22 (www.mirbase.org), and alterations in the nomenclature were corrected by miRBase tracker (www.mirbasetracker.org).
Target microRNA prediction
The miRNA-targeting mRNAs were predicted by TargetScanHuman v7.2 (http://www.targetscan.org/vert_72/). This algorithm predicts biological targets of miRNAs by searching for the presence of conserved 8 mer, 7 mer, and 6 mer sites that match the seed region of each miRNA. It also identifies mismatches in the seed region that is compensated by conserved 3' pairing  and centered sites. The species was selected as “cow,” and the miRNA full name as found in miRBase was given as the input. Predictions were ranked based on the predicted efficacy of targeting as calculated using cumulative weighted context++ scores of the sites.
Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses
To unify the molecular, biological, cellular, and pathway level functionalities of tissue-specific bTB miRNAs, the predicted genes were submitted to Database for Annotation, Visualization, and Integrated Discovery Bioinformatics Resources 6.8 (http://david.abcc.ncifcrf.gov/home.jsp), a freely available integrated biological knowledge base and analytical tool. The gene list with the identifier “OFFICIAL_GENE_SYMBOL” was uploaded by selecting B. taurus as the background. The three gene ontology (GO) classifiers, namely, biological process (BP), molecular function (MF), and cellular components (CC) were applied. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the genes were also carried out. The statistical significance of each term/pathway was evaluated by setting the EASE score (modified Fisher's exact P value) to ≤ 0.05 and the respective functional class of the three most enriched pathways was recorded.
| Results|| |
Clustering of tissue-specific microRNAs
An extensive literature surfing was conducted for the 2010–2018 period to yield ten miRNAs reportedly upregulated in bTB with respect to healthy controls [Table 1]. These miRNAs were clustered by their tissue-specific induction into three groups – G1 (mammary gland-specific, one miRNA): bta-miR-146a; G2 (PBMC-specific, one miRNA): bta-miR-155; and G3 (alveolar macrophage-specific, ten miRNAs): bta-miR-146a, bta-miR-155, bta-miR-142-5p, bta-miR-423-3p, bta-miR-21-5p, bta-miR-27a-3p, bta-miR-99b, bta-miR-147, bta-miR-223, and bta-let-7i. All ten miRNAs were validated to be upregulated in the bovine alveolar macrophage (G3).
|Table 1: Retrieved microRNAs upregulated in bovine tuberculosis as compared to healthy controls|
Click here to view
MicroRNA-predicted target genes
The predicted gene counts for G1, G2, and G3 were 246, 498, and 3379, respectively. G3 had the highest number of gene records due to the participation of relatively higher number of miRNAs in the group.
Kyoto Encyclopedia of Genes and Genomes pathway enrichment
For G1, G2, and G3 genes, 57, 232, and 692 chart records (pathways) were characterized. The topmost enriched three channels by G1 were pathways in cancer, human T-cell leukemia virus type 1 infection, and T-cell receptor signaling [Table 2] while the G2-G3 gene sets augmented pathways in cancer, MAPK, and PI3K-Akt signaling [Table 3] and [Table 4]. Pathways in cancer were common to all three groups (G1, G2, and G3), while MAPK and PI3K-Akt signaling pathways were pertinent to G2 and G3 members. The fibroblast growth factors (FGFs) were identified as the core genes nourishing the prime signaling cascades in G2 and G3, while for G1, the key players were clusters of differentiation (CD) molecules.
|Table 2: Top 3 Kyoto Encyclopedia of Genes and Genomes functional pathways enriched by target genes of G1 microRNAs|
Click here to view
|Table 3: Top 3 Kyoto Encyclopedia of Genes and Genomes functional pathways enriched by target genes of G2 microRNAs|
Click here to view
|Table 4: Top 3 Kyoto Encyclopedia of Genes and Genomes functional pathways enriched by target genes of G3 microRNAs|
Click here to view
Gene ontology annotation
The GO terms for the gene groups were detailed. The BP of G1, G2, and G3 is marked as the regulation of transcription from RNA polymerase II. They bind to metal – zinc ions, adenosine triphosphate (ATP), and DNA (MF). They mostly found in the nucleus and cytoplasm (CC) [Table 5], [Table 6], [Table 7] for respective GO terms of G1, G2, and G3].
|Table 5: Top 3 gene ontology biological process, molecular function, and cellular components terms enriched by target genes of G1 microRNAs|
Click here to view
|Table 6: Top 3 gene ontology biological process, molecular function, and cellular components terms enriched by target genes of G2 microRNAs|
Click here to view
|Table 7: Top 3 gene ontology biological process, molecular function, and cellular components terms enriched by target genes of G3 microRNAs|
Click here to view
| Discussion|| |
This work is designed to disclose the pathogenic cascades and regulatory units (genes) enriched by ten upregulated bTB miRNAs. Our roadmap included grouping of miRNAs in a tissue-specific manner ensured by modeling the interaction network between referred miRNAs and their target genes and elucidation of signal transduction pathways.
Selection of a reliable target gene prediction tool is the first step in the workflow. Although a number of gene prediction algorithms (TargetScan, Diana Tools, TarBase, mirPath, miExTra, Miranda, PicTar, PITA, rna22, miRBridge, and miRSystem) are available, only TargetScan has the provision to predict miRNA target genes in B. taurus (bovine or cow). Most prediction algorithms follow the basic attributes of complementarity between seed region (miRNA) and target sites (mRNA), seed conservation, and miRNA/mRNA duplex energy.
The KEGG annotation accentuated that pathways in cancer were well-furnished by bTB miRNA target genes of all three tissue locations. These results were similar to a previous study in humans where tissue-specific miRNAs overexpressed in TB were involved in the regulation of multiple cancer pathways. Barh et al. recently disclosed the coexistence of TB and lung cancer through noncoding RNAs.
We identified MAPK and PI3K-Akt signaling as importantly enriched pathways common to bTB miRNAs of PBMC and alveoli macrophages with FGFs being the major enrichment factor. MAPK and PI3K-Akt pathways are usually triggered by the activation of Ras, a small GTPase, that plays a key role in fine-tuning the innate immune response. The MAPK pathway is involved in boosting the proliferation and survival of immune cells while PI3K works toward the activation and functioning of immune cells. Accidental activation of these pathways in a bacterial infection may, therefore, influence the state of innate immune system. Vegh et al. already corroborated the involvement of bovine alveolar macrophage upregulated miRNAs in innate immunity. All these expositions indicate that bTB contagion is established through a complex network of bTB miRNA-driven immune system alteration through different genes.
The GO structuring of the transcriptomic data revealed that upregulated bTB miRNAs mostly target transcription regulators from RNA polymerase II promoter with their prime MF being binding to metal ions, zinc ion, ATP, and DNA. They inferred to occupy the nucleus and cytoplasm.
MiRNAs regulate major signaling pathways and hence crucial to maintain the homeostasis. Aberrant expression of these molecules is indicative of pathological states of many infectious diseases. The potential application of miRNAs as clinical, diagnostic, prognostic, and therapeutic follow-up biomarkers for TB has extensively been explored in humans.,,,, However, very little is known about the role of miRNAs in bovine TB pathophysiology. Application of in vitro-in vivo methods to discern this enigma is cumbersome where in silico (computational) approaches form a good alternative stance. In silico platform generates large data sets which constitute an effective baseline for planned experiments. The present work intends to give a comprehensive understanding on the role of tissue-specific bTB transcripts in modulating the pathogenicity by inciting multiple signaling pathways. Our preliminary data need to be ascertained through realistic biological models.
The major limitation to present study is the use of a single miRNA target gene prediction algorithm (TargetScan). The practice of gene finding by two or more sophisticated tools could increase the accuracy and specificity of the results. Since TargetScan is the only tool with a provision for direct miRNA–mRNA interaction prediction in the species B. taurus, we pursued with this one.
| Conclusion|| |
This study demonstrated a direct link between bTB miRNAs and oncogenic cascades. Pathways in cancer were common to G1, G2, and G3 classes, whereas MAPK and PI3K-Akt signaling were specific to PBMC (G2) and alveolar macrophage (G3) miRNAs. FGFs were the key determinants in G2–G3 while CD molecules regulate the G1 functions. These genes could be novel therapeutic targets for anti-TB treatments. The GO classifiers disclosed that bTB miRNAs modulate the transcription from RNA polymerase II promoter (BP), target genes functionally involved in the binding interactions with ions and organic compounds (MF), and are confined to nucleus and cytoplasm (CC). The upregulated bTB transcripts presumed to complete the pathogenic circuit by governing specific molecular pathways. The oncogenic MAPK and PI3K-Akt signaling also connected to innate immune response. An addendum to this work with consolidated experimental proof could provide vital clues on the mechanistic and epigenetic roles of bTB miRNAs.
The authors would like to thank the Institute of Agricultural Research for Development (IRAD), Bambui, Cameroon, for providing the platform for the realization of this work.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Iannaccone M, Cosenza G, Pauciullo A, Garofalo F, Proroga YT, Capuano F, et al.
Milk microRNA-146a as a potential biomarker in bovine tuberculosis. J Dairy Res 2018;85:178-80.
Vegh P, Magee DA, Nalpas NC, Bryan K, McCabe MS, Browne JA, et al.
MicroRNA profiling of the bovine alveolar macrophage response to Mycobacterium bovis
infection suggests pathogen survival is enhanced by microRNA regulation of endocytosis and lysosome trafficking. Tuberculosis (Edinb) 2015;95:60-7.
Golby P, Villarreal-Ramos B, Dean G, Jones GJ, Vordermeier M. MicroRNA expression profiling of PPD-B stimulated PBMC from M. Bovis-challenged unvaccinated and BCG vaccinated cattle. Vaccine 2014;32:5839-44.
Vegh P, Foroushani AB, Magee DA, McCabe MS, Browne JA, Nalpas NC, et al.
Profiling microRNA expression in bovine alveolar macrophages using RNA-seq. Vet Immunol Immunopathol 2013;155:238-44.
Qi Y, Cui L, Ge Y, Shi Z, Zhao K, Guo X, et al.
Altered serum microRNAs as biomarkers for the early diagnosis of pulmonary tuberculosis infection. BMC Infect Dis 2012;12:384.
Yi Z, Fu Y, Ji R, Li R, Guan Z. Altered microRNA signatures in sputum of patients with active pulmonary tuberculosis. PLoS One 2012;7:e43184.
Fu Y, Yi Z, Wu X, Li J, Xu F. Circulating microRNAs in patients with active pulmonary tuberculosis. J Clin Microbiol 2011;49:4246-51.
Liu Y, Wang X, Jiang J, Cao Z, Yang B, Cheng X, et al.
Modulation of T cell cytokine production by miR-144* with elevated expression in patients with pulmonary tuberculosis. Mol Immunol 2011;48:1084-90.
Ndzi EN, Nkenfou CN, Mekue LM, Zentilin L, Tamgue O, Yone Pefura EW, et al
. MicroRNA hsa-miR-29a-3p is a plasma biomarker for the differential diagnosis and monitoring of tuberculosis. Tuberculosis 2018;114:69-76.
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]
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]
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. MiRBase: MicroRNA sequences, targets and gene nomenclature. Nucleic Acids Res 2006;34:D140-4.
Van Peer G, Lefever S, Anckaert J, Beckers A, Rihani A, Van Goethem A, et al.
MiRBase tracker: Keeping track of microRNA annotation changes. Database (Oxford) 2014;2014. pii: bau080.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. eLife 2015;4:e05005. eLife Lens view.
Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005;120:15-20.
Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 2009;19:92-105.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009;4:44-57.
Ndzi EN, Ambily Nath IV, Nair AS, Chellappan Biju V, Shidhi PR, Tamgue O, et al
. Bioinformatics probing of tissue-specifc miRNAs overexpressed in active tuberculosis reveals induction of multiple cancer pathways. J Appl Bioinforma Comput Biol 2018;7:2.
Barh D, Tiwari S, Kumavath RN, Azevedo V. Potential link between tuberculosis and lung cancer through non-coding RNAs. BioRxiv preprint 2017. [doi: 10.1101/188375].
Johnson DS, Chen YH. Ras family of small GTPases in immunity and inflammation. Curr Opin Pharmacol 2012;12:458-63.
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6], [Table 7]