Email updates

Keep up to date with the latest news and content from BMC Microbiology and BioMed Central.

Open Access Research article

Identification of latent tuberculosis infection-related microRNAs in human U937 macrophages expressing Mycobacterium tuberculosis Hsp16.3

Qing-Lin Meng13, Fei Liu2, Xing-Yuan Yang4, Xiao-Mei Liu13, Xia Zhang2, Chun Zhang1* and Zong-De Zhang2*

Author Affiliations

1 Suzhou Municipal Key Laboratory of Molecular Diagnostics and Therapeutics, Suzhou Institute of Biomedical Engineering and Technology, Chinese Academy of Sciences, Suzhou 215163, Jiangsu, China

2 Beijing Key Laboratory of Drug Resistance Tuberculosis, Beijing Tuberculosis and Thoracic Tumor Research Institute, Beijing Chest Hospital, Capital Medical University, Beijing 101149, China

3 Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun 130033, Jilin, China

4 College of Life Science and Technology, China Pharmaceutical University, Nanjing 210009, Jiangsu, China

For all author emails, please log on.

BMC Microbiology 2014, 14:37  doi:10.1186/1471-2180-14-37


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2180/14/37


Received:19 September 2013
Accepted:10 February 2014
Published:12 February 2014

© 2014 Meng et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Abstract

Background

Latent tuberculosis infection (LTBI) relies on a homeostasis of macrophages and Mycobacterium tuberculosis (Mtb). The small heat shock protein, Mtb Hsp16.3 (also known as latency-associated antigen), plays an important role in Mtb persistence within macrophages. However, the mechanism of LTBI remains elusive. The aim of this study was to delineate LTBI-related miRNA expression in U937 macrophages expressing Mtb Hsp16.3 protein. U937 macrophages were infected with an integrase-deficient Lentivirus vector to transiently express Mtb Hsp16.3, and green fluorescent protein (GFP) as a control. We used a microRNA (miRNA) microarray chip containing more than 1000 probes to identify the significant differentially expressed miRNAs in the infected U937 cells, and employed real-time quantitative polymerase chain reaction (qRT-PCR) for validation. Furthermore, we confirmed these candidate LTBI-related miRNAs in peripheral blood mononuclear cells from subjects with LTBI and in healthy control individuals. Functional annotation prediction of miRNA target genes and pathway enrichment analyses were used to explore the putative links between these miRNAs and LTBI.

Results

Analysis of the miRNA expression profile identified 149 miRNAs that were differentially expressed in U937 macrophages expressing Mtb Hsp16.3 compared with the control expressing GFP. The expression level of seven miRNAs (miR-424-5p, miR-493-5p, miR-296-5p, miR-27b-3p, miR-377-5p, miR-3680-5p, miR-191-5p) were validated by qRT-PCR. The expression level of four miRNAs (miR-424-5p, miR-27b-3p, miR-377-5p, miR-3680-5p) in the peripheral blood mononuclear cells samples from LTBI and healthy participants reflected the altered patterns observed in the microarray profile. The bioinformatic analyses suggest that the miRNAs may regulate Mtb latent infection by affecting the development of macrophage cells.

Conclusions

The results suggest that miRNA expression may play a considerable role in the pathogenesis of LTBI, and this would increase our understanding of the molecular basis of Hsp16.3-facilitated Mtb survival in macrophages.

Keywords:
microRNAs; Macrophages; Mycobacterium; Tuberculosis; Small heat shock protein; Latent tuberculosis infection

Background

In the latent tuberculosis infection (LTBI) state, the patient harbors the Mycobacterium tuberculosis (Mtb) bacilli in the body and is asymptomatic; no radiographic or bacteriological evidence of active tuberculosis is observed; however, the patients reveal immunological sensitization to Mtb-derived antigen proteins (e.g., ESAT6, CFP10, and Hsp16.3) [1]. The granuloma is thought to play a major role in maintaining latency and avoiding reactivation of Mtb, representing the intersection of innate and adaptive immunity. The hypoxic core of the granuloma is thought to induce a dormant state of Mtb. In this regard, in vitro studies have confirmed that Mtb dramatically upregulated dormancy survival regulon (DosR)-related genes, which are characteristic of nonreplicating persistence [2]. One of the most prominent of these is Rv2031c, which encodes the small heat shock protein Hsp16.3 (also known as α-crystalline related protein 1, or the 16 kDa antigen). Hsp16.3 constitutes one of the prominent antigens in the exponential phase after infection. It contains both T- and B-cell epitopes that contributes to enhance the cellular and humoral immune responses [3]. Interestingly, Hsp16.3 is maximally expressed during latency, play a role in facilitating the persistence of Mtb within macrophages [4]. Indeed, the functional versatility of macrophages is evident from their role in diverse biological processes, such as phagocytosis, inflammation, immunoregulation, differentiation, and metabolism [5]. Recent studies of these cells using system biology and a variety of -omics technologies in several disease models (e.g., atherosclerosis and metabolic disorders) suggest that they orchestrate crucial functions during homeostasis or pathogenesis in health/disease [6].

MicroRNAs (miRNAs) are endogenous, 22–25 nucleotide RNAs that play major regulatory roles in higher eukaryotes by targeting mRNAs for cleavage or translational repression. MiRNAs modulate the innate and adaptive immune responses to pathogens by affecting host immune cell differentiation and progression of diseases [7]. The clinical application of miRNAs as diagnostic or prognostic biomarkers has already been demonstrated in various types of cancers [8]. However, compared with their well-known role in cancer, the biological and diagnostic role of miRNAs in LTBI is still poorly understood. In the present study, we used U937 cell line as in vitro macrophage model, focused on the interaction between U937 macrophages and Mtb Hsp16.3, aiming to identify differentially expressed miRNAs in U937 macrophages. Our study intends to explore the potential function of miRNAs in the interaction of macrophages with Mtb Hsp16.3 and provide insights for investigating the role of macrophage homeostasis in LTBI.

Methods

Ethics statement and participants

The local ethics committee of the Beijing Tuberculosis and Thoracic Tumor Research Institute reviewed and approved the study. Written informed consent was obtained from participants before their enrollment in the study. Twenty clinical health care workers of Beijing Chest Hospital were recruited and all have history of close contact with active tuberculosis patient for more than two years. The four healthy controls were students of Suzhou Institute of Biomedical Engineering and Technology and had no history of contact with TB. Potential study participants were excluded if they had another infectious disease. The interferon gamma release assay (IGRA) (T-SPOT.TB, Oxford Immunuotec, Oxfordshire, UK) was used to distinguish the LTBI group from healthy control. Fourteen clinical health care worker participants were IGRA-positive and included as LTBI group while the four healthy control subjects were IGRA-negative.

PBMC samples preparation

Peripheral venous blood (10 ml) was drawn from each subject and PBMC samples were isolated by density gradient separation using Lympholyte-H, immediately mixed with TRIzol (1 ml) and frozen at -80°C until RNA were extracted.

Preparation of the IDLV and Infection

To obtain the Mtb Hsp16.3 expression vector pLVHsp-IRES-GFP, the encoding gene Rv2031c was amplified and cloned into the pLVX-IRES-GFP plasmid, and confirmed by sequencing. The Lenti-X HTX Packaging System (Integrase Deficient) (Clontech, Mountain View, CA, USA) was used to prepare the viral vector. The U937 cells were cultured in RPMI1640 medium (Gibco, Grand Island, NY, USA) containing 10% fetal bovine serum under 5% CO2 at 37°C, infected with viral IDLVs stock at 5:1 multiplicity of infection (MOI), refreshed with medium 6 h later and incubated for 64 h.

Western blot analysis

Briefly, U937 cells were infected with IDLVs (Hsp/GFP), and control IDLVs (GFP), respectively. After 64 h, the cells were collected and then heated for 5 min at 95°C in 1 × protein loading buffer containing β-mercaptoethanol, and cell extracts were separated on 12% SDS-PAGE gel and transferred to PVDF membranes. The membranes were blocked with 5% skimmed milk-TBST, incubated with polyclonal rabbit anti-Mtb Hsp16.3 antibodies (prepared in our laboratory) and β-actin antibodies (Solarbio, Beijing, China), and finally visualized via reaction with a chemiluminescent HRP substrate (Millipore, Billerica, MA, USA).

Microarray analyses of infected macrophages

KangCheng Biosciences (Shanghai, China) performed the miRNA profiling analysis. To determine the miRNA profiles for the two groups, total RNAs were purified using TRIzol (Invitrogen, Grand Island, NY, USA) and a miRNeasy mini kit (Qiagen, Shenzhen, China), labeled using the miRCURY™ Hy3™/Hy5™ Power labeling kit (Exiqon, Vedbaek, Denmark) and hybridized on the specific miRCURY™ LNA Array (v.18.0, Exiqon, Denmark) platform. The Exiqon miRCURY™ LNA Array (v.18.0) contains 2043 capture probes covering all human miRNAs, and could quantify genome-wide miRNA expression in the two groups. Images on the chip were scanned using an Axon GenePix 4000B microarray scanner (Axon Instruments, Foster City, CA, USA) and imported into GenePix Pro 6.0 software (Axon) for grid alignment and data extraction. MiRNAs with intensities >50 were used to calculate the normalization factor. Expression data were normalized using the median normalization. After normalization, average values of replicate spots of each miRNA were used for statistical analysis; differentially expressed miRNAs were identified through fold change filtering. Data are presented as means ± standard deviations. Analysis of variance tests or unpaired two-tailed Student t tests were used for statistical analysis. The data were regarded as significantly different at P < 0.05.

Reverse transcription and quantitative real time-polymerase chain reaction (qRT-PCR) validation

The total RNAs were extracted from each two groups of infected U937 macrophages and PBMC samples using a mirVana™ miRNA Isolation Kit (Ambion, Austin, TX, USA). cDNA was reverse transcribed from total RNAs using the miRcute miRNA cDNA first-strand synthesis kit (Tiangen, Beijing), according to the manufacturer’s instructions. Using U6/5S RNA as the endogenous reference for normalization, qRT-PCR assays were performed on an ABI 7500 Real-Time PCR System (Applied Biosystems, Foster, CA, USA) using the miRcute miRNA qPCR Detection kit (SYBR Green) (Tiangen, Beijing, China). The experiments were conducted in triplicate.

Pathway enrichment analyses

The predicted targets of the miRNAs were obtained from the TargetScan database [9], and the PITA database [10]. The intersections of the results obtained from these different software programs were regarded as the reliable target genes. The predicted miRNA target genes were analyzed for enriched KEGG pathways using the NCBI DAVID server ( http://david.abcc.ncifcrf.gov webcite) with default settings [11].

Results

U937 Macrophages expressed Mtb Hsp16.3 and GFP, respectively

To reduce the risk of insertional mutagenesis in U937 cells, the IDLV system was used to produce non integrative lentiviral vectors , which delivered the transgene into U937 macrophages for instantaneous expression. After 64 h, the fluorescence microscopy was used to detect the GFP expression, and GFP-positive cells was quantified by the flow cytometry. It was shown that the transformation efficiency of the test group of U937 cells (expressing Mtb Hsp16.3) was 73% (Figure  1A), and that of the control group (expressing GFP) was 82% (Figure  1B). Furthermore, to validate the expression of Mtb Hsp16.3 protein in the cells, western blot analysis was performed using anti-Mtb Hsp16.3 and the results demonstrated that Mtb Hsp16.3 was strongly expressed in the test group of U937 cells (Figure  1C).

thumbnailFigure 1. The integrase-deficient lentivirus vector (IDLV) transfected U937 cells with high efficiency and the cells expressed Mtb Hsp16.3. An IDLV delivered the transgene into U937 macrophages for instantaneous expression. The fluorescence microscopy and flow cytometry were used at 64 h after infection to detect GFP and analyse the transduction efficiency. A, the transduction efficiency of the test group of U937 cells (expressing Mtb Hsp16.3 and GFP) was 73%. B, the transduction efficiency of the control group (expressing GFP only) was 82%. C, western blot analysis with antibodies against Mtb Hsp16.3; β-actin was used as a loading control.

Expression profiles of miRNAs in U937 cells from the test group and the control group

To determine the miRNA profiles for the two groups, the Exiqon miRCURY™ LNA Array was employed to perform the 2043 miRNAs assay (1898 human and 145 human viral miRNAs represented in the Sanger miRBase v18.0). After normalization and unsupervised filtering (see Methods), the obtained average values for each miRNA spot were used for statistical analysis. Comparing the data from the two groups (test/control) and using fold change filtering (upregulated more than 2-fold and downregulated less than 0.5-fold ), total of 149 differentially expressed miRNAs was identified, of which 60 were upregulated (Table  1) and 89 were downregulated (Table  2). The P values for these 149 miRNAs were less than 0.05 in the test groups compared to results for the control groups.

Table 1. Summary of upregulated miRNAs

Table 2. Summary of downregulated miRNAs

qRT-PCR validation of candidate miRNA expression level

To validate the microarray findings, seven miRNAs were selected for qRT-PCR analysis. The miR-424-5p (previous ID: miR-424), miR-493-5p (previous ID: miR-493*), and miR-296-5p were reported as potential to discriminate between latent TB and healthy by the previous study [12], the other four miRNAs (miR-27b-3p, miR-377-5p, miR-3680-5p, miR-191-5p) were randomly selected. As shown in Figure  2A, the respective level of downregulated miR-27a-3p, miR-424-5p, and miR-493-5p in qRT-PCR results largely reflected the altered patterns of these selected miRNAs observed in the microarray profiles. In parallel, the levels of upregulated miR-296-5p, miR-377-5p, miR-3680-5p, and unchanged miR-191-5p were similar to the chip results as well (Figure  2B). Furthermore, to evaluated the relative expression level of the six differentially expressed miRNAs in LTBI group and healthy control, 14 LTBI subjects and four healthy control individuals were recruited for the qRT-PCR assay (Additional file 1: Table S1). As shown in Figure  3, the results of four miRNAs (miR-424-5p, miR-27a-3p, miR-377-5p, miR-3680-5p) recapitulated the microarray data, and the other two miRNAs (miR-493-5p and miR-296-5p) were not significant differentially expressed.

thumbnailFigure 2. Confirmation of miRNA expression profiles of the microarray by qPCR. After normalization to 1 in the control group (U937/GFP), the relative expressions of selected downregulated miRNAs (miR-27a-3p, miR-424-5p, and miR-496-5p) in the test group are shown in A; the relative expressions of upregulated miRNAs (miR-296-5p, miR-377-5p, and miR-3680-5p), and unchanged miR-191-5p in the test group are shown in B.

Additional file 1: Table S1. Characteristics of latent TB infection and healthy control participants used for qRT-PCR.

Format: DOCX Size: 18KB Download fileOpen Data

thumbnailFigure 3. qPCR validation of miRNA expression levels in samples from the latent tuberculosis infection (LTBI) group versus the healthy control group. Relative expressions of miR-424-5p, miR-496-5p, miR-27a-3p, miR-377-5p, and miR-3680-5p in LTBI and healthy samples. Statistical analysis was performed using the unpaired t-test. **P < 0.01, *P < 0.05, NS: not significant.

Bioinformatic exploratory analyses

Because only very small number of miRNA targets has been experimentally validated so far, we predicted the targets of these differentially expressed miRNAs using TargetScan and PITA software (Additional file 2: Table S2). The identified miRNAs were predicted to modulate 7044 target genes. We then used the NCBI DAVID server to identify the significantly enriched pathways involving the predicted target genes. As shown in Table  3, apart from cancer-associated pathways, the MAPK signaling, endocytosis, Wnt signaling, focal adhesion, axon guidance, and TGF-beta signaling pathways, which are related to differentiation, polarization, and versatility of macrophages, were significantly enriched. The results suggest that the miRNAs may regulate Mtb infection by affecting the development of immune cells.

Additional file 2: Table S2. Target genes of differently expressed miRNAs.

Format: XLSX Size: 2.7MB Download fileOpen Data

Table 3. Enriched pathways involving the predicted target genes

Discussion

The macrophage is the main replication niche of Mtb, despite the bactericidal characteristics and functions that this cell type normally has. The Mtb has evolved several strategies to reside and even replicate within the otherwise hostile environment of the macrophage, including the prevention of phagosome-lysosome fusion, inhibition of phagosomal maturation, and detoxification of the host’s stresses. Accordingly, the localization of Mtb inside the macrophage has been a matter of debate in recent years [13]. For a long time, an impermeable phagosome in the macrophage was thought to contain Mtb. However, recent evidence indicates that Mtb, as well as M. leprae, can escape its vacuole and reside in the host cell cytosol [14]. It is becoming clear that LTBI is not a static state with a homogenous population of non-replicating bacilli, but a constant endogenous Mtb reinfection process [15]. It is argued that both phagosomal maturation inhibition and escape from the phagosome are part of the survival strategies of Mtb. Owing to the lack of widely accepted animal or in vitro models to reflect the complex biology of LTBI, population-based studies are considered the best methods to investigate the intricate nature of LTBI. In some previous reports, human cell line U937 was used as in vitro model to investigate the molecular mechanism of Mtb during infection or persistence and its effect on the cell [16,17]. In this study, U937 cells expressing Hsp16.3 in the cytosol could partialy reflect the dynamic interplay of macrophages with dormant Mtb, which is necessary to prevent reactivation of the bacilli and development of active TB.

Indeed, some miRNAs that have been previously linked to carcinogenesis of different organs and tissues, such as miR-424-5p (previous ID: miR-424), miR-221-5p (previous ID: miR-221*), miR-675, miR-647, miR-125a-5p, miR-214-3p (previous ID: miR-214), miR-130b-3p (previous ID: miR-130b), miR-522-3p (previous ID: miR-522), and miR-16-5p (previous ID: miR-16) [18-21] were found to be up- or downregulated in our analysis. Forrest and colleagues [22] showed that induction of miR-424 (miR-424-5p) and miR-222 (miR-222-3p) promotes monocytic differentiation via combined regulation; both of these miRNAs were significantly downregulated in this analysis. Interestingly, miR-150-5p (previous ID: miR-150) has been shown to regulate the immune response and monocyte differentiation [23]; miR-150-5p was upregulated in our analysis. Conversely, miR-181a (miR-181a-5p) and miR-146a (miR-146a-5p), which have been proven to participate in the regulation of the adaptive immune responses, were 7- and 10-fold downregulated in our profiling data [24,25]. Furthermore, current research has demonstrated that miR-181a regulates inflammation responses in macrophages, and increased expression of miR-181a is strongly correlated with the expression of interleukin (IL)-1β, IL-6, and tumor necrosis factor alpha (TNFα) [26]. These results suggest that Hsp16.3 protein might be involved in blocking immunity against Mtb via miR-181a and miR-146a deregulation. In addition, Fu et al. demonstrated that miR-93*(miR-93-3p) was the most upregulated in active TB serum [27]; however, our analysis indicated that miR-93-3p was downregulated, making it a potential diagnostic marker to distinguish latent TB from active TB.

Although many target genes have been predicted by bioinformatic methods, the functions of most differentially expressed miRNAs remain unknown, and very few predicted target genes have been validated. More than half of the differentially expressed miRNAs did not find a target mRNA in either database; most of them were recently identified miRNAs. Bioinformatic exploratory provides a rapid analytic approach categorizing large amounts of genes into functionally related groups to thereby facilitate the uncovering of the biological content captured by transcriptomic profiling. KEGG pathway enrichment analyses further interpret the biological functions of these genes. The overrepresented pathways associated with glioma and basal cell carcinoma were enriched, which somewhat surprised. The formation of extra-pulmonary tuberculosis is caused by Mtb dissemination to other organs or tissue, such as the central nervous system (CNS) [28]. It is speculated that the occurrence of glioma might be related to tuberculoma in the CNS, and a tuberculoma-like granuloma is often misdiagnosed as a tumor [29]. This indicated that Mtb Hsp16.3 might be involved in carcinogenesis, which warrants further investigation.

Earlier studies, which used peripheral blood mononuclear cells (PBMCs) or whole blood cells to perform whole genome transcriptional profiling and miRNA profiling [27,30], described a number of candidate biomarkers that might function in active TB. Wang and collegues identified miRNAs that were differently expressed in latent TB versus healthy from the clinical PBMC samples [12], In present study, the microarray data and independent qRT-PCR results indicated that our in vitro model by used of U937 cells expressing Mtb Hsp16.3 protein has good repeatability. However, the weakness of the model is also obvious, it does not represent the real interaction of pathogen and host macrophage in vivo, it provided only mechanistic insights on the interaction between Mtb antigen and human cell line. Although the expressions of miR-424-5p (previous ID: miR-424), miR-27a-3p, miR-377-5p and miR-3680-5p were consistent in clinical PBMC samples, the small size of healthy controls weakened the statistical power. Our understanding the biology of latent tuberculosis as part of a broad range of responses that occur following infection with Mtb remains incomplete. Multiple factors are involved in this complex process. Herein, compared to previously studies, our experiments got more differentially expressed miRNAs since we focused on just whether the Mtb Hsp16.3 had great effects on the U937 macrophage cell. Furthermore, this model could also be used in the follow-up investigation of the miRNA candidates regulating the macrophage in chronic inflammatory response or other process correlated with LTBI.

Conclusions

Using miRNA expression profiling, we identified 149 differentially expressed miRNAs and validated that the transcription patterns of some miRNAs were consistent with previous reports. Our data provide evidences for the underlying biological processes involved in LTBI via the interaction between U937 macrophages and the Mtb Hsp16.3 protein. These findings provide an improved understanding of the link between miRNA homeostasis and LTBI. Further characterization of the pathogenetic roles of specific miRNAs and deciphering of the miRNA-controlled signaling regulatory network may help to enhance diagnosis and prevention of LTBI.

Supporting information

Microarray data submission for human arrays is MIAME-compliant. The chip data from this study have been deposited at NCBI Gene Expression Omnibus (GEO) database, and its accession number is GSE54630.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CZ and ZDZ conceived the study. QLM, FL, XYY, XML, and XZ carried out the experiments. QLM wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This work was supported by the National Major Project (grant 2013ZX10003003) and the Suzhou Science and Technology Project (grant ZXJ2012005).

References

  1. Tufariello JM, Chan J, Flynn JL: Latent tuberculosis: mechanisms of host and bacillus that contribute to persistent infection.

    Lancet Infect Dis 2003, 3:578-590. PubMed Abstract | Publisher Full Text OpenURL

  2. Yuan Y, Crane DD, Simpson RM, Zhu YQ, Hickey MJ, Sherman DR, Barry CE 3rd: The 16-kDa alpha-crystallin (Acr) protein of Mycobacterium tuberculosis is required for growth in macrophages.

    Proc Natl Acad Sci U S A 1998, 95:9578-9583. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Leyten EM, Lin MY, Franken KL, Friggen AH, Prins C, van Meijgaarden KE, Voskuil MI, Weldingh K, Andersen P, Schoolnik GK, et al.: Human T-cell responses to 25 novel antigens encoded by genes of the dormancy regulon of Mycobacterium tuberculosis.

    Microbes Infect 2006, 8:2052-2060. PubMed Abstract | Publisher Full Text OpenURL

  4. Yuan Y, Crane DD, Barry CE 3rd: Stationary phase-associated protein expression in Mycobacterium tuberculosis: function of the mycobacterial alpha-crystallin homolog.

    J Bacteriol 1996, 178:4484-4492. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Mueller P, Pieters J: Modulation of macrophage antimicrobial mechanisms by pathogenic mycobacteria.

    Immunobiology 2006, 211:549-556. PubMed Abstract | Publisher Full Text OpenURL

  6. Biswas SK, Chittezhath M, Shalova IN, Lim JY: Macrophage polarization and plasticity in health and disease.

    Immunol Res 2012, 53:11-24. PubMed Abstract | Publisher Full Text OpenURL

  7. Lodish HF, Zhou B, Liu G, Chen CZ: Micromanagement of the immune system by microRNAs.

    Nat Rev Immunol 2008, 8:120-130. PubMed Abstract | Publisher Full Text OpenURL

  8. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al.: MicroRNA expression profiles classify human cancers.

    Nature 2005, 435:834-838. PubMed Abstract | Publisher Full Text OpenURL

  9. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets.

    Cell 2003, 115:787-798. PubMed Abstract | Publisher Full Text OpenURL

  10. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E: The role of site accessibility in microRNA target recognition.

    Nat Genet 2007, 39:1278-1284. PubMed Abstract | Publisher Full Text OpenURL

  11. da Huang W, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists.

    Nucleic Acids Res 2007, 35:W169-W175. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Wang C, Yang S, Sun G, Tang X, Lu S, Neyrolles O, Gao Q: Comparative miRNA expression profiles in individuals with latent and active tuberculosis.

    PLoS One 2011, 6:e25832. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Welin A, Lerm M: Inside or outside the phagosome? The controversy of the intracellular localization of Mycobacterium tuberculosis.

    Tuberculosis (Edinb) 2012, 92:113-120. PubMed Abstract | Publisher Full Text OpenURL

  14. van der Wel N, Hava D, Houben D, Fluitsma D, van Zon M, Pierson J, Brenner M, Peters PJ: M. tuberculosis and M. leprae translocate from the phagolysosome to the cytosol in myeloid cells.

    Cell 2007, 129:1287-1298. PubMed Abstract | Publisher Full Text OpenURL

  15. Cardona PJ: A dynamic reinfection hypothesis of latent tuberculosis infection.

    Infection 2009, 37:80-86. PubMed Abstract | Publisher Full Text OpenURL

  16. McGarvey JA, Wagner D, Bermudez LE: Differential gene expression in mononuclear phagocytes infected with pathogenic and non-pathogenic mycobacteria.

    Clin Exp Immunol 2004, 136:490-500. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Samuel LP, Song CH, Wei J, Roberts EA, Dahl JL, Barry CE 3rd, Jo EK, Friedman RL: Expression, production and release of the Eis protein by Mycobacterium tuberculosis during infection of macrophages and its effect on cytokine secretion.

    Microbiology 2007, 153:529-540. PubMed Abstract | Publisher Full Text OpenURL

  18. Lui WO, Pourmand N, Patterson BK, Fire A: Patterns of known and novel small RNAs in human cervical cancer.

    Cancer Res 2007, 67:6031-6043. PubMed Abstract | Publisher Full Text OpenURL

  19. Garofalo M, Quintavalle C, Romano G, Croce CM, Condorelli G: miR221/222 in cancer: their role in tumor progression and response to therapy.

    Curr Mol Med 2012, 12:27-33. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Jiang L, Huang Q, Zhang S, Zhang Q, Chang J, Qiu X, Wang E: Hsa-miR-125a-3p and hsa-miR-125a-5p are downregulated in non-small cell lung cancer and have inverse effects on invasion and migration of lung cancer cells.

    BMC Cancer 2010, 10:318. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Finnerty JR, Wang WX, Hebert SS, Wilfred BR, Mao G, Nelson PT: The miR-15/107 group of microRNA genes: evolutionary biology, cellular functions, and roles in human diseases.

    J Mol Biol 2010, 402:491-509. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Forrest AR, Kanamori-Katayama M, Tomaru Y, Lassmann T, Ninomiya N, Takahashi Y, de Hoon MJ, Kubosaki A, Kaiho A, Suzuki M, et al.: Induction of microRNAs, mir-155, mir-222, mir-424 and mir-503, promotes monocytic differentiation through combinatorial regulation.

    Leukemia 2010, 24:460-466. PubMed Abstract | Publisher Full Text OpenURL

  23. Xiao C, Calado DP, Galler G, Thai TH, Patterson HC, Wang J, Rajewsky N, Bender TP, Rajewsky K: MiR-150 controls B cell differentiation by targeting the transcription factor c-Myb.

    Cell 2007, 131:146-159. PubMed Abstract | Publisher Full Text OpenURL

  24. Li QJ, Chau J, Ebert PJ, Sylvester G, Min H, Liu G, Braich R, Manoharan M, Soutschek J, Skare P, et al.: miR-181a is an intrinsic modulator of T cell sensitivity and selection.

    Cell 2007, 129:147-161. PubMed Abstract | Publisher Full Text OpenURL

  25. Perng DW, Yang DM, Hsiao YH, Lo T, Lee OK, Wu MT, Wu YC, Lee YC: miRNA-146a expression positively regulates tumor necrosis factor-alpha-induced interleukin-8 production in mesenchymal stem cells and differentiated lung epithelial-like cells.

    Tissue Eng Part A 2012, 18:2259-2267. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Xie W, Li M, Xu N, Lv Q, Huang N, He J, Zhang Y: miR-181a regulates inflammation responses in monocytes and macrophages.

    PLoS One 2013, 8:e58639. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Fu Y, Yi Z, Wu X, Li J, Xu F: Circulating microRNAs in patients with active pulmonary tuberculosis.

    J Clin Microbiol 2011, 49:4246-4251. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Rock RB, Olin M, Baker CA, Molitor TW, Peterson PK: Central nervous system tuberculosis: pathogenesis and clinical aspects.

    Clin Microbiol Rev 2008, 21:243-261. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Del Brutto OH, Mosquera A: Brainstem tuberculoma mimicking glioma: the role of antituberculous drugs as a diagnostic tool.

    Neurology 1999, 52:210-211. OpenURL

  30. Jacobsen M, Repsilber D, Gutschmidt A, Neher A, Feldmann K, Mollenkopf HJ, Ziegler A, Kaufmann SH: Candidate biomarkers for discrimination between infection and disease caused by Mycobacterium tuberculosis.

    J Mol Med (Berl) 2007, 85:613-621. PubMed Abstract | Publisher Full Text OpenURL