Transcriptional Profile of Mycobacterium tuberculosis in an in vitro Model of Intraocular Tuberculosis

Introduction

Intraocular tuberculosis (IOTB) is a unique extrapulmonary manifestation (Gupta et al., 2016) of the disease in TB endemic as well as non-endemic countries (Patel et al., 2013; Lee et al., 2016). No prospective data are available on the population-based prevalence of IOTB either from the low-endemic or the high-endemic TB regions in the world (Gupta et al., 2016). However, among the uveitis clinic population, the prevalence of TB as an etiology of uveitis has varied from 0.5% in the USA, 4% in the People's Republic of China, 6.31% in Italy, 6.9% in Japan, 9.86% in north India, 10.5% in Saudi Arabia, and 11.4% in Iraq where TB is endemic (Shakarchi, 2015). Altogether, the prevalence ranges from 1–4% (non-endemic) to 10–26% (endemic) in different areas of the world (Lee et al., 2016). A recent report from India showed that 5.2% of patients with uveitis also had sputum positive pulmonary or extrapulmonary TB (Gogia et al., 2018). Another retrospective multinational cohort study including 801 patients from 25 ophthalmology referral centers demonstrated good treatment response in TB uveitis patients (Agrawal et al., 2017). More importantly, IOTB is becoming a major health problem, similar to other types of extrapulmonary TB (Dalvin and Smith, 2017), not only due to the increasing prevalence but also due to the identification of drug-resistant cases (Sharma et al., 2018).

Intraocular tuberculosis is a significant cause of inflammation and visual morbidity, particularly in TB endemic countries (Gupta et al., 2015; Shakarchi, 2015). It poses clinical, diagnostic, and therapeutic challenges due to the protean clinical presentations of IOTB within a single organ (Gupta et al., 2015). Unavailability of adequate human samples from the eye, an immune-privileged site, and failure in detecting the presence of live bacteria further makes it difficult to understand the pathogenesis of IOTB (Shechter et al., 2013; Gupta et al., 2015; Lee et al., 2016). In view of the increasing number of cases of IOTB along with the diagnostic dilemma associated with this disease, it is very important to study the molecular pathogenesis of IOTB using the relevant experimental models. Studies using experimental TB models, such as animal and cell lines, have already been proven useful for understanding the various aspects of pulmonary as well as extrapulmonary forms of TB (Talaat et al., 2004; Jain et al., 2006; Be et al., 2008; Doycheva et al., 2010; Krishnan et al., 2010; Scordo et al., 2016; Fonseca et al., 2017; Zhan et al., 2017). Mycobacterium tuberculosis primarily localizes in the lung and is taken up by the alveolar macrophages which are also involved in the transport of bacilli by the hematogenous route (Henderson et al., 1963; Balasubramanian et al., 1994; Danelishvili et al., 2003) to various other organs where it remains dormant until it gets activated and produces extrapulmonary TB disease (Tufariello et al., 2003; Barrios-Payán et al., 2012). So far it is not known how and where on reaching the eye, M. tuberculosis is localized and activates sight-threatening inflammation/uveitis. Although recent clinical reports highlight that M. tuberculosis can affect any tissue of the eye, primarily the posterior part of the eye is involved due to high oxygen tension (Dalvin and Smith, 2017; Moharana et al., 2018). The late-stage IOTB has been found to occur in retina as retinitis and retinal vasculitis (Doycheva et al., 2010; Gupta et al., 2015), and in a clinical sample representing granulomatous uveitis, acid-fast bacilli (AFB) have been shown to be present in the retinal pigment epithelium (RPE) cells (Rao et al., 2006). Thus, the RPE cells—the non-professional phagocytic cells in the eye—have been considered as a probable host for the survival and replication of M. tuberculosis (Gupta et al., 2007), and reactivation of these sequestered organisms may lead to the recurrence of IOTB (Patel et al., 2013). Studies on the intracellular M. tuberculosis in both alveolar macrophages (professional) and alveolar epithelial (non-professional) cells have indicated that (Danelishvili et al., 2003) soon after invasion, M. tuberculosis gets localized in a cytoplasmic compartment known as phagosomes, and acquires the fusion with late endosomal/lysosomal markers but inhibits the biogenesis of phagolysosomes for its intracellular survival (Hasan et al., 1997; Huynh et al., 2007).

Recently, a human fetal retinal pigment epithelial cell line model has been reported to phagocytize an avirulent form of mycobacteria, H37Ra, with the involvement of host toll-like receptors (Nazari et al., 2014). However, a cell line model for IOTB using virulent M. tuberculosis (H37Rv) for elucidating the key determinants (genes/proteins) of disease pathogenesis is still lacking. Identification of transcriptional signatures of the bacteria is key to understanding how the bacteria colonize, invade, or replicate in host microenvironment (Talaat et al., 2004) and has been proven to be a valid approach in understanding the pathogenesis of TB (Ward et al., 2010). Transcriptional changes in pathogen (M. tuberculosis) have been studied in well-defined in vitro culture systems (Lin et al., 2016) as they efficiently mimic the host environment and this has improved our knowledge on the host-pathogen interactions (Nandy et al., 2018). These studies have led to the identification of transcripts and pathways specific to growth conditions involved in the survival of the pathogen in the unfavorable lysosomal microenvironment (Lin et al., 2016). Thus, transcriptomic knowledge is possibly valuable for the design of new drugs, vaccines, and new approaches for controlling TB (Aguilar-Ayala et al., 2017). At present, various critical key gaps exist in the field of IOTB as it is still not known how the mycobacteria alter its transcriptional profile to support its survival and replication in the ocular environment.

Thus, the current study is an attempt to investigate some of these aspects. We have developed a well-defined in vitro RPE cell line model of IOTB and utilized it to examine the intracellular fate of M. tuberculosis using late endosomal/lysosomal cell surface markers and M. tuberculosis transcriptome to identify the key determinants involved in the pathogenesis of IOTB.

Materials and Methods

Cell Line and M. tuberculosis Strain

All the experiments were carried out in accordance with the rules and regulations set forth by the Institute's Ethics Committee (IEC) and Institute's Biosafety Committee (IBC), Post Graduate Institute of Medical Education and Research (PGIMER), Chandigarh. Human samples were obtained from the ocular samples repository for ophthalmology research (#NK/3323/study/186).

M. tuberculosis H37Rv was grown in Sauton's media (0.5 g of potassium dihydrogen phosphate, 0.5 g of magnesium sulfate, 2 g of citric acid, 0.05 g of ferric ammonium citrate, and 4 g of L-asparagine in 900 mL of double distilled water. Once dissolved, 60 mL of glycerol and 0.05% Tween-80 were added and pH of the medium was adjusted to 7.2–7.4 with sodium hydroxide solution).

Adult retinal pigment epithelium-19 (ARPE-19) (ATCC^® CRL-2302™) cell line was maintained in Dulbecco's Modified Eagle Medium Nutrient Mixture F-12 (DMEM/F-12) media (Thermo Scientific #12400024) supplemented with sodium bicarbonate (1.2 g/L), penicillin (100 U/ml), streptomycin (100 μg/ml), 10% FBS, and 5% carbon dioxide (CO₂) at 37°C. The ARPE-19 cells were resuspended in fresh growth media without antibiotics prior to M. tuberculosis infection.

In vitro Model of Intraocular Tuberculosis

H37RV Inoculum

Single cell suspension of H37Rv was prepared using log phase H37Rv culture grown in Sauton's media with minor modifications (Napier et al., 2011). The H37Rv culture was centrifuged at 3,000 rpm for 5 min and the supernatant was removed. The pellet was resuspended and washed twice using phosphate buffer saline supplemented with 0.05% Tween-80 buffer (PBS-T; pH-7.2). The washed pellet thus obtained was suspended in 20 ml of PBS-T buffer with 3 mm glass beads (Merck), vortexed for 15 s to disperse clumps, and the suspension was kept standing for 10 min at room temperature. Ten milliliters of this suspension was collected in fresh 50 ml falcon tubes, mixed with fresh 3 mm of glass beads, vortexed again for 15 s, followed by standing for 10 min at room temperature. This step was repeated three times and the single cell suspension of H37Rv thus obtained was matched with MacFarland standards, and the appropriate dilutions were prepared in fresh antibiotic-free DMEM/F-12 medium for infecting the ARPE-19 cells.

Adherence, Invasion, and Multiplication of M. tuberculosis H37Rv in ARPE-19 Cells

The standard methodology for in vitro infection of epithelial cells with H37Rv was followed with minor modifications (Mehta et al., 1996; Napier et al., 2011). The ARPE-19 monolayer was washed (three times) with antibiotic-free DMEM/F-12 medium followed by addition of fresh antibiotic-free DMEM/F-12 medium supplemented with 2% FBS. The cells were infected with a single cell suspension of H37Rv in the same medium at a multiplicity of infection (MOI) of 10 bacteria per cell (10:1) for 2-h (adherence time-point). The adherence of M. tuberculosis to ARPE-19 cell line was observed through scanning electron microscopy (SEM) and quantified by colony forming units (CFU) assay (Mehta et al., 1996; Ramsugit and Pillay, 2016). Two-hour post infection, the monolayers of the infected and uninfected cells were washed three times with warm antibiotic-free DMEM/F-12 medium and replaced by fresh antibiotic-free DMEM/F-12 medium supplemented with 2% FBS and amikacin (20 μg/mL) to kill the extracellular bacilli at 37°C under 5% CO₂ for 2-h (0-day post infection; dpi). At this time-point, the invasion of bacilli in the ARPE-19 cells was monitored by transmission electron microscopy (TEM) and quantitatively confirmed by CFU enumeration. Following 2 h of amikacin treatment, all the cultured cells (infected and uninfected) were maintained in fresh antibiotic-free DMEM/F-12 medium supplemented with 1% FBS and amikacin (2 μg/mL) (to preclude any extracellular growth of M. tuberculosis) to monitor the intracellular replication at 1−, 3−, and 5-dpi by the CFU assay. Both the infected and uninfected ARPE-19 cells were incubated at 37°C under 5% CO₂ at all the time-points.

Scanning Electron Microscopic Examination

The ARPE-19 cell line was cultured at a cell density of 2 × 10⁵ cells/well on 8 mm coverslips in 24-well tissue culture plates. At 15- and 30-min of infection, both the infected and uninfected ARPE-19 cells on the coverslips were fixed using a fixative solution [glutaraldehyde 2.5% (v/v), paraformaldehyde fixative 2% (w/v), phosphate buffer 0.1 M] for 24 h. After fixation, the coverslips were washed with PBS, dehydrated using ethanol, air dried, platinum coated, and observed under JSM-IT300LV scanning electron microscope (JEOL, Tokyo Japan) at Central Sophisticated Instruments Cell (CSIC), PGIMER, Chandigarh, India. The H37Rv inoculum used for infection was also fixed on coverslips and processed for SEM imaging as a positive control.

Transmission Electron Microscopic Examination

To visualize the internalized bacilli within the RPE cells, TEM studies were conducted. The ARPE-19 cells were seeded at a density of 1 × 10⁶ cells in T75 tissue culture flasks. At 0-, 1-, and 3-dpi, the infected and uninfected ARPE-19 cells were washed, de-adhered, and the cell pellets were collected by centrifugation at 2,000 rpm for 10 min. The cell pellets were then fixed in fixative solution in 0.2 M Sorensen's phosphate buffer (sodium dihydrogen phosphate 0.2 M, sodium hydrogen phosphate 0.2 M, sucrose 7.5%, and glutaraldehyde 4%) for 2–4 h followed by a mild spin and fixation in 1% osmium tetroxide and dehydration in ethanol series. The cells were then infiltrated with EPON mixture and propylene oxide (1:1) for 2 h at room temperature and then embedded in EPON resin, followed by polymerization at 60°C for 24 h. Sections of 0.5 μm thickness were cut using an ultramicrotome (Reichert-Jung, Leica, Germany) and stained with 0.5% toluidine blue to confirm the presence of cells. Further, 60 nm ultra-thin sections were cut and mounted on nickel grids (300 mesh). The sections were double-stained with uranyl acetate and lead citrate and then examined by TEM at the Department of Histopathology, PGIMER, Chandigarh, India and photographed. The H37Rv inoculum used for infection during the study was also fixed and processed similarly for TEM analysis.

Colony Forming Units Enumeration

The ARPE-19 cell line was cultured at a cell density of 2 × 10⁵ cells/well in 24-well tissue culture plates. The infected and uninfected ARPE-19 cells at the indicated time points were lysed by treatment with 0.5 ml of 0.1% Triton X-100 for 5 min, and serial dilutions of the infected and non-infected cell lysates (1:10, 1:100, and 1:1,000) prepared in 7H9 broth were plated on 7H11 agar supplemented with 10% OADC. H37Rv used during infection (inoculum) was also plated at similar dilutions. Colonies were counted after 3–4 weeks of incubation at 37°C/5% CO₂. The colonies of the respective time points at each dilution were counted and averaged (from two independent experiments with two independent wells at each time-point) and the percent adherence and invasion were calculated with respect to the initial inoculums as explained in the legend of Table 2. Intracellular multiplication was calculated in terms of fold CFU (mean CFU at each time-point divided by the day 0 CFU) (Scharn et al., 2016).

Cell Viability Assay

The effect of M. tuberculosis infection (MOI 10:1) on the viability of ARPE-19 cells at different time-points was studied by MTT (3-[4,5-methylthiazol-2-yl]-2,5-diphenyl-tetrazolium bromide) assay. The RPE cells (10⁴ cells/100 μL/well) were plated in a 96-well tissue culture plate. The uninfected and infected ARPE-19 cells (at 0-, 1-, 3-, and 5-dpi) were washed and incubated with 100 μL of fresh DMEM/F-12 medium containing 10 μL of MTT solution (5 mg/mL) at 37°C for 4 h. After incubation, 100 μL of a freshly prepared solution containing 0.01 M HCl and 10% SDS was added to the wells to dissolve formazan crystals. Following overnight incubation, the absorbance was measured on an ELISA plate reader at 570 nm, and the viability of the uninfected cells at the respective time-points was considered as 100%. At each time-point, the wells containing DMEM/F-12 medium, MTT and SDS/HCl solution without ARPE-19 cells were used as a blank. The data presented are the mean ± SEM of O.D. at A570 (uninfected and infected wells), and percent cytolysis was calculated by dividing the mean A570 of the test wells by mean A570 of the control cells × 100.

Intracellular Fate of M. tuberculosis in the Retinal Pigment Epithelial Cells

To understand the intracellular fate of M. tuberculosis inside the RPE cells, colocalization of M. tuberculosis and late endosomal/lysosomal markers (LAMP-1, lysosomal-associated membrane protein-1 and LAMP-2, lysosomal-associated membrane protein-2) were studied during the course of infection in the ARPE-19 cells (Hasan et al., 1997; Huynh et al., 2007).

Labeling of M. tuberculosis

The H37Rv was labeled before infection with PKH26 dye using MINI 26 Sigma kit which stably incorporates a yellow-orange fluorescent dye (PKH26) into the lipid regions of the cell membrane. A single cell suspension containing 2 × 10⁷ bacilli was washed and suspended in a 2X cell suspension solution (1 mL of diluent C and 1 mL 2X dye solution). The mix was incubated for 1–5 min with periodic mixing at ambient temperature (20–25°C) in the dark. After the completion of labeling, the bacterial pellets were washed twice with PBS and resuspended in antibiotic-free DMEM/F-12 medium for infection of ARPE-19 cells.

Confocal Microscopy

The ARPE-19 cells were cultured on 14 mm glass coverslips in a 24-well tissue culture plate at a cell density of 2 × 10⁵ cells and were infected with PKH26-labeled H37Rv at MOI (10:1) as explained above. Post 2 h of amikacin (20 μg/ml) treatment, the infected and uninfected monolayer cells were washed three times with warm antibiotic-free DMEM/F-12 medium and the coverslips were collected to study the intracellular bacilli at different time points (0-, 1-, and 3-dpi). The ARPE-19 monolayer cells on the coverslips were fixed with 4% paraformaldehyde for 15 min, permeabilized with 0.05% saponin in PBS for 5 min, and blocked with 5% serum in PBS for 2 h at room temperature or overnight at 4°C. The cells were further incubated with primary antibodies for endosomal/lysosomal markers, i.e., anti-human -LAMP-1 and -LAMP-2 at a concentration of 1 and 2 μg/ml, respectively for, 2 h followed by washing thrice with PBS and incubation with FITC labeled anti-rabbit IgG secondary antibody. The ARPE-19 cells were washed again 3 times and 4′6′-diamidino-2-phenylindole dihydrochloride (DAPI) was added to the cells at a concentration of 2 ng/ml and incubated for 1 min to stain the nuclei. The coverslips were mounted in glycerol: saline (v/v, 1:1) on glass slides and the images were acquired with 100X oil immersion objective and standard filters using a confocal microscope (FLUOVIEW FV1000). Further images were captured and processed using FV10-ASW 3.0 software (Olympus) followed by colocalization analysis using Coloc2-ImageJ plugin. All the steps were performed at ambient temperature (20–25°C). Three coverslips at each time-point were analyzed each for infected and uninfected cells.

Transcriptome Profiling of M. tuberculosis in in vitro Model of Intraocular Tuberculosis

The ARPE-19 monolayers were infected with M. tuberculosis at MOI of 10:1 as explained above and 3-dpi RNA was isolated from intracellular bacilli. A pool of five T225 flasks, each containing 2 × 10⁷ monolayer cells, was infected with H37Rv, and the total RNA was isolated as one sample experiment to study the whole M. tuberculosis genome through microarray. The RNA isolated from M. tuberculosis H37Rv inoculum used for infection was taken as the reference for comparison and calculation of the change in mycobacterial gene expression in ARPE cells. The data were obtained and analyzed from three independent experiments.

M. tuberculosis Ribonucleic Acid Isolation

For the isolation of M. tuberculosis RNA from the infected ARPE-19 cells, the eukaryotic RNA was removed by treating the infected monolayer cells with GTC solution (4 M guanidine thiocyanate, 0.5% sodium-N-lauryl sarcosine, 25 mM tri-sodium citrate, 0.1 M β-2-mercaptoethanol, and 0.5% Triton-80) for 2–5 min. The lysate was centrifuged at 5,000 g for 20 min at 4°C to collect the pellet to which 1–2 mL of the TRIzol solution was added followed by addition of 0.1 mm zirconia beads for bead-beating. The TRIzol solution was collected by centrifugation at 16,000 g for 1 min at 4°C to which 1/5th volume of chloroform was added followed by incubation at RT for 10 min. After centrifugation, 250 μL of chilled isopropanol and 250 μL 0.1 M sodium acetate were added to the aqueous phase and kept at −20°C for 2 h. The samples were thereafter centrifuged at 12,000 g for 15 min at 4°C to pellet down the precipitates of RNA. Washing of RNA pellet was done using 1 mL of 75% chilled ethanol followed by centrifugation at 7,500 g for 5 min at 4°C. The pellet was air dried for 3–5 min followed by addition of nuclease-free water to the RNA pellet. RNA was quantitated and treated with Turbo DNase (2 U/μl; Ambion) as per manufacturer's instructions. The DNase-treated RNA of all the samples was purified using the Mini Elute Clean-up kit (Qiagen) and RNA amount, purity, and quality were determined by nanodrop and Agilent 2100 bioanalyzer.

Microarray Sample Preparation and Array Processing

Input RNA (100 ng) from experimental and control samples were amplified using T7 RNA polymerase and simultaneously labeled with Cyanine 3-CTP using One-Color Microarray-Based Low Input Quick Amp WT Labeling Protocol (Agilent Technologies). Using the RNeasy Mini kit, the labeled/amplified cRNA samples were purified and quantitated through NanoDrop ND-1000 UV-VIS Spectrophotometer. Cyanine 3 dye concentration (pmol/μL), RNA absorbance ratio (260/280 nm), and cRNA concentration (ng/μL) were recorded to determine the yield and specific activity of each reaction. The yield (μg of cRNA) was calculated using the formula, concentration of cRNA × 30 μL (elution volume)/1,000, and the specific activity (pmol Cy3 per μg cRNA) was calculated as the concentration of Cy3/concentration of cRNA × 1,000. The samples with a specific activity of more than 15 were used for hybridization to custom design the microarray slides to study the M. tuberculosis gene expression profile. One microgram of labeled cRNA was used from each experimental and control sample for hybridization to Agilent whole M. tuberculosis Genome Oligo Microarrays slides (G2509F, GE 8X15K, printed with Agilent SurePrint^® technology, USA) using Agilent Microarray Hybridization Chamber Kit. After hybridization, the array slides were washed and processed for scanning on the Agilent Sure Scan Microarray Scanner using one-color scan setting for 8X15K array slides. The scanned image (.tiff) file was analyzed with Agilent Feature Extraction (FE) software (version 12.0.3.1) using default parameters. The raw data for each hybridization (experimental and control) as well as the final processed data (normalized data) have been submitted as per Minimum Information about a Microarray Experiment (MIAME) guidelines and is available at the Gene Expression Omnibus website (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE115292) with accession number GSE115292.

Microarray Data Analysis

The.tiff image file was utilized to generate a quality control (QC) report for each sample (experimental and control) as per the Agilent FE protocol. For each sample, from FE software, QC metric sets, thresholds, and charting tools were determined, and the resulting text files (.txt) of the control and experimental samples were analyzed for data normalization, statistical significance, and gene expression analysis using Agilent Technologies GeneSpring^® GX software (version 14.9). Significantly modulated genes were selected based on Benjamini–Hochberg procedure (Benjamini and Hochberg, 1995, 2000; Keselman et al., 2002; Pawitan et al., 2005) with a fold change >2 and p < 0.05, based on previous microarray studies with M. tuberculosis (Sharma et al., 2017).

Validation of M. tuberculosis

Transcripts M. tuberculosis gene-specific primers (Supplementary Table 1) were designed using NCBI Primer-BLAST tool. The microarray expression data was validated through quantitative RT-PCR on Qiagen Rotor-Gene Q real-time PCR machine using SYBR Green dye. Unamplified RNA from experimental (two independent experiments) and control (H37Rv) samples were used for cDNA synthesis (Thermo cDNA synthesis kit, #AB1453). Each sample cDNA along with the specific no enzyme control (NEC) was run in duplicate using specific cycling conditions for each gene. M. tuberculosis 16s gene transcript level was used as endogenous control. After completion, real-time software (Rotor-Gene Q Series Software) was used to determine the amplification plot (C_T, cycle threshold) and to analyze the melt curve for the product specificity. Relative quantitation of target gene expression was carried out using the comparative C_T method (ΔΔC_T) (Livak and Schmittgen, 2001). The C_T of the transcripts which showed the difference of less than 5C_T in comparison to its respective NEC was corrected using the formula (Laurell et al., 2012)

$\begin{array}{l}\left\{-\text{Lo}{\text{g}}_2(2^{\wedge }-\text{RT}+)-(2^{\wedge }-\text{RT}-)\right\},\end{array}$

where,

RT+ is C_T of transcripts observed in cDNA sample,

RT– is C_T of transcripts observed in no enzyme control (NEC) sample.

Further, Graph-Pad Prism software was used for unpaired two-tailed t-test for statistical analysis and p < 0.05 was considered statistically significant.

Mycobacterium tuberculosis Transcripts in Human Intraocular Tuberculosis Samples

We were further interested to study if the M. tuberculosis transcriptional signatures identified in the in vitro model of IOTB are also present under in vivo conditions to evaluate their utility as potential molecular biomarkers in human IOTB. For this, qualitative real-time PCR was used as a tool which is considered to be advantageous over conventional PCR (Peng et al., 2016) and can also be manipulated as a diagnostic approach. Thus, three M. tuberculosis transcripts (Rv1230c, Rv1971, and Rv3872) observed to be upregulated in ARPE-19 cells were also examined in the human vitreous samples from two different categories (possible and confirmed IOTB). Two hundred microliters of each vitreous sample (n = 10; possible = 3, and confirmed = 7) and H37Rv (positive control) was utilized for total RNA isolation as explained above. The cDNA was synthesized (Thermo cDNA synthesis kit, #AB1453) from isolated RNA along with no enzyme control (NEC) reactions. The two reactions differed from each other in terms of presence (reverse transcriptase, RT+) or absence (no enzyme control NEC, RT–) of the reverse transcriptase enzyme. Additionally, no template controls were also run. All the samples along with the positive and negative controls were then analyzed for specific melt curves and C_T values. A C_T value < 35 was taken as positive amplification. The amplification was considered due to DNA (D+) if the C_T value in RT– sample was lesser than or equal to RT+ sample and due to RNA (R+) if the C_T in RT+ sample was lesser than or equal to RT– sample.

Results

In the present study, an in vitro RPE cell line model of IOTB has been established to monitor the adherence, invasion, intracellular replication, and the fate of M. tuberculosis. Using this in vitro model of IOTB, the whole genome microarray was carried out to determine the M. tuberculosis transcriptional signatures to understand the key determinants involved in the survival of mycobacteria in ocular environment and pathogenesis of IOTB.

In vitro Model of Intraocular Tuberculosis

M. tuberculosis Adheres, Invades, and Multiplies in the ARPE-19 Cell Line

The adherence of M. tuberculosis to RPE cells was initially studied by SEM analysis at 15- and 30- min post-infection. The scanning micrographs (Figure 1) of the uninfected ARPE-19 cells represented an intact monolayer (Figure 1A), while the micrographs of infected ARPE-19 cells showed M. tuberculosis adhering to monolayer cells within 15-min of infection (Figure 1B). At 30-min post-infection, cell surface changes and disruption particularly at the site of M. tuberculosis attachment were observed (Figure 1C).

FIGURE 1

Figure 1. Scanning electron microscopy images of uninfected and M. tuberculosis infected ARPE-19 cells. (A) Electron micrograph of uninfected ARPE-19 monolayer cells. (B) The white arrow indicates adherent M. tuberculosis bacilli at 15 min post infection and the inset shows SEM micrographs of in vitro grown M. tuberculosis used for infection. (C) White arrows indicate M. tuberculosis adhering to ARPE-19 cells at 30 min and cell surface changes at the site of attachment of bacilli. The experiment was done in three independent wells and the representative figures are shown.

Based on the surface interaction of H37Rv with the ARPE-19 cells, further studies were planned to monitor the invasion of M. tuberculosis in ARPE-19 cells using TEM analysis. The mycobacterial cells were visualized as electron dense and electron transparent layer with thick capsular outer layer beyond the cell wall, typical of mycobacterial species as reported earlier (Takade et al., 2003). The uninfected ARPE-19 cell line (0-day) showed a large number of cytoplasmic extensions, high nucleus to cytoplasm ratio, and scanty organelles (Figure 2A). Post amikacin treatment (0-dpi, invasion), multiple bacilli were found in a single vacuolar compartment, indicating internalization/phagocytosis of M. tuberculosis by the ARPE-19 cells (Figure 2B). At 1-dpi, bacilli were similarly located inside vacuoles in the cytoplasm (Figure 2C) whereas, at 3-dpi, micrographs showed an increase in the number of vacuolar compartments with multiple bacilli (Figure 2D). Additionally, at 1- and 3-dpi disrupted cellular morphology with degenerated cytoplasm was evident in the infected ARPE cells (Figures 2C,D) while no such changes were observed in the morphology of uninfected ARPE-19 cells at 1- and 3-day (data not shown).

FIGURE 2

Figure 2. Transmission electron microscopy images of uninfected and M. tuberculosis infected ARPE-19 cells. (A) Electron micrograph showing a cross-section of uninfected ARPE-19 cells. (B–D) Infected ARPE-19 cells after 0- (B), 1- (C), and 3-dpi (D). The figure (B) in inset represents M. tuberculosis H37Rv cells. Black arrows indicate the vacuolar membrane surrounding the multiple bacilli in ARPE-19 cells. The experiment was done twice and the cells were collected from independent flasks, and the representative figures are shown. dpi, day post infection.

Additionally, the adherence and invasion of M. tuberculosis to ARPE-19 cells were quantitatively confirmed by CFU enumeration. The adherence of M. tuberculosis to ARPE-19 cells was found to be 22.7% of the initial inoculum [86.2 ± 33.8 (× 10³) CFU] within 2 h of incubation with total CFU of 19.7 ± 6.3 (×10³). The percent invasion of M. tuberculosis in the RPE cell line as studied after the completion of amikacin treatment (0-dpi) indicated that ARPE-19 cells had taken up 8.4% of the initial inoculum (Table 1).

TABLE 1

Table 1. Mycobacterium tuberculosis adherence, invasion, and replication in in vitro model of intraocular tuberculosis.

CFU enumeration at later time points of infection demonstrated the intracellular multiplication of M. tuberculosis in ARPE-19 cells. The infected ARPE-19 cells showed 5.7 ± 0.8 (×10³) CFU at 1-dpi, 30.2 ± 10.4 (×10³) CFU at 3-dpi, and 47.8 ± 12.1 (×10³) CFU at 5-dpi (Table 1 and Figure 3A). The increase in intracellular M. tuberculosis growth over time was calculated by an increase in fold CFU at 1-, 3-, and 5-dpi in comparison to day 0.The bacterial replication was found to be increased by four-fold at 3-dpi while it was increased up to 6.5-fold at 5-dpi, and this increase in bacterial growth was significant (p < 0.05) in comparison to 0-dpi (Figure 3B).These results indicate that M. tuberculosis H37Rv adheres, invades, and utilizes RPE cells for intracellular replication.

FIGURE 3

Figure 3. Colony forming units (CFU) and fold CFU of M. tuberculosis in ARPE-19 cells at different time-points after infection. ARPE-19 monolayer cells were infected with H37Rv at MOI 10:1. (A) CFU counts are mean ± standard error mean (SEM). (B) Fold CFU was calculated by dividing the mean CFU at each time-point (1-, 3-, and 5-dpi) by the day 0 CFU. **p < 0.01 and *p < 0.05 in comparison to 0-dpi. Two independent experiments were performed and cells were infected in duplicate for each time-point. dpi, day post-infection.

Cytotoxicity of M. tuberculosis to Retinal Pigment Epithelial Cell Line

To further confirm the observations obtained from electron microscopy and to determine the effect of intracellularly multiplying bacilli on the viability of ARPE-19 cells, an MTT assay was performed. The mean optical density (O.D.) obtained at 570 nm (A570) for uninfected and infected ARPE-19 cells at different time-points are plotted in Figure 4A. A significant decrease in A570 at the start of intracellular replication (1 dpi; p < 0.05) with further decrease at 3 dpi, as compared to the uninfected cells (p < 0.05), indicated a loss of cell viability in the infected cell line. Under similar conditions, A570 was maintained in case of uninfected cells till 5 dpi indicating no significant change in cell viability. In terms of percent cytolysis, M. tuberculosis caused 3.5, 8.8, 28.8, and 36.9% loss in ARPE-19 cell viability at 0−, 1−, 3−, and 5-dpi, respectively, with significant loss of cells occurring at 1-dpi (p < 0.05) and 3-dpi (p < 0.05) in comparison to their respective controls (Figure 4B).

FIGURE 4

Figure 4. Percent cytotoxicity caused by M. tuberculosis to ARPE-19 cells. (A) Mean ± SEM O.D. values at 570 nm (A570) of uninfected and infected ARPE-19 cells at different time-points (0-,−1-, 3-, 5-dpi). (B) Cytolysis of cells was quantified using MTT dye at different time-points. *p < 0.05 in comparison to 0- dpi. The percent (%) cytotoxicity was calculated as = (Mean O.D. of control cells - Mean O.D. of infected cells)/Mean O.D. of control cells × 100. dpi, day post infection; SEM, standard error of the mean.

Intracellular Fate of M. tuberculosis in Retinal Pigment Epithelial Cells

As the ability of mycobacteria to survive in host cells is attributed to colocalization with endosomal/lysosomal compartments (Hasan et al., 1997; Huynh et al., 2007), after confirming the adherence, invasion, and intracellular multiplication of M. tuberculosis in RPE cells, we were further interested to determine the intracellular fate of M. tuberculosis. For this purpose, the ARPE-19 cells were infected with PKH26-labeled M. tuberculosis bacilli followed by confocal microscopic analysis at 0-, 1-, and 3-dpi.

Mycobacterium tuberculosis Colocalizes With LAMP-1and LAMP-2 in ARPE-19 Cell Line

To examine the localization of M. tuberculosis to late endosomes/lysosomes, we used antibodies against LAMP-1/LAMP-2 proteins, the most abundant proteins in lysosomal membranes. Post amikacin treatment (0-dpi), the PKH26-labeled bacilli showed partial colocalization with both LAMP-1 (Figure 5) and LAMP-2 (Figure 6) marker proteins (green: LAMP-1/LAMP-2 and red: M. tuberculosis) as revealed by an analysis performed using Coloc2-ImageJ plugin that provided Pearson's coefficient (R) values of 0.897 and 0.797, respectively (where zero is no colocalization, 1 means perfect colocalization). With LAMP-1, the R values corresponded to 0.730 and 0.503 at 1- and 3-dpi, respectively (Figure 5), whereas with LAMP-2, the R-values corresponded to 0.536 and 0.858 at 1- and 3-dpi, respectively (Figure 6). A number of intracellular bacilli were also found in free form, without any colocalization with LAMP-1/LAMP-2 at 3-dpi (Figures 5, 6; row D, column iv), suggesting the survival and high fold replication of M. tuberculosis in the cytoplasm of ARPE-19 cells. Overall, these results suggest that following invasion (0-dpi), M. tuberculosis gets recruited to the endocytic pathway and partially colocalizes with the late endosome/lysosomes markers during its replication at 1- and 3-dpi within the RPE cells. A similar type of intracellular fate of M. tuberculosis has also been reported in other professional and non-professional phagocytic cells (Hasan et al., 1997; Huynh et al., 2007).

FIGURE 5

Figure 5. LAMP-1 staining for lysosomal localization of M. tuberculosis in RPE (ARPE-19) cells. Confocal micrographs of uninfected (UI) (A) and infected (I) ARPE-19 cells (MOI 10:1) at 0-dpi (B), 1-dpi (C), and 3-dpi (D) showing staining for (i) M. tuberculosis labeled with PKH26 dye (red channel); (ii) LAMP-1 antibody using secondary FITC-IgG (green channel); (iii) DAPI for nuclei (blue channel). In infected ARPE-19 cells (rows B–D; column iv), white arrowheads show merging of two channels (red and green) and the Pearson's R coefficient for the colocalization of M. tuberculosis with LAMP-1. Zero is no colocalization, and 1 means perfect colocalization. Dpi, day post-infection; UI, uninfected; I, infected.

FIGURE 6

Figure 6. LAMP-2 staining for lysosomal localization of M. tuberculosis in RPE (ARPE-19) cells. Confocal micrographs of uninfected (UI) (A) and infected (I) ARPE-19 cells (MOI 10:1) at 0-dpi (B), 1-dpi (C), and 3-dpi (D) showing staining for (i) M. tuberculosis labeled with PKH26 dye (red channel); (ii) LAMP-2 antibody using secondary FITC-IgG (green channel); (iii) DAPI for nuclei (blue channel). In infected ARPE-19 cells (rows B–D; column iv), white arrowheads show merging of two channels (red and green) and the Pearson's R coefficient for the colocalization of M. tuberculosis with LAMP-2. Zero is no colocalization, and 1 means perfect colocalization. Dpi, day post-infection; UI, uninfected; I, infected.

Transcriptional Profile of Intracellular M. tuberculosis in in vitro Model of Intraocular Tuberculosis

The decrease in viable CFU count at 1-dpi in comparison to 0-dpi (Figure 3A) followed by a steady increase in intracellular CFU load till 5-dpi (Figure 3B) indicates the adaptation of few intracellular bacilli after initial infection of ARPE-19 cells. This survival and adaptation of M. tuberculosis in an in vitro RPE model of infection could be due to altered bacterial physiological behavior that could be inferred through global transcriptome analysis (Rohde et al., 2012). Additionally, such analysis could also be significant for identifying the metabolic checkpoints which are involved in successively averting the host hostile environment and may be exploited for therapeutics for avoiding the replication of M. tuberculosis (Du et al., 2016).

Transcriptional Signatures of M. tuberculosis Inside RPE Cells

We analyzed M. tuberculosis transcriptional signatures by whole genome microarray of intracellular bacilli at 3-dpi in the ARPE-19 cells. The heat map shows the signal intensities in H37Rv control (rows 1, 2) and infected (rows 3–8) experimental sets (Figure 7A). Hierarchical clustering of microarray data demonstrated that the infected and control samples were overlapping to their respective technical sets, and the principal components analysis further confirmed the hierarchical clustering as the correlation coefficient value of each data set were similar and were thus treated as experimental replicates (Supplementary Figure 1). Among all the genes encoded by M. tuberculosis genome, the transcripts corresponding to 1,040 genes were found to be differentially regulated (cut off 2.0) in comparison to H37Rv control. Out of these differentially regulated transcripts, 469 transcripts were found to be upregulated whereas 571 transcripts were downregulated. The volcano graph (Figure 7B) shows the upregulated (blue box), downregulated (red box), and unchanged (gray boxes) transcripts. The upregulation and downregulation of differentially regulated transcripts were in the range of +1.15 to +10.1 Log₂ fold change and −4.2 to −1.15 Log₂ fold change, respectively (Supplementary Table 2; GEO accession number # GSE115292). Comparison of differentially regulated transcripts with the earlier reported data (Jain et al., 2006; Rohde et al., 2012; Ryndak et al., 2015) from other in vitro models of M. tuberculosis infection indicated that many of these transcripts were uniquely regulated in in vitro model of IOTB (Supplementary Table 2).

FIGURE 7

Figure 7. Transcriptome analysis of intracellular M. tuberculosis in in vitro model of intraocular tuberculosis. (A) The heat map shows the 2 control samples (in vitro M. tuberculosis as inoculum used during infection; rows 1, 2) and 6 test samples (three biological replicates, each with two technical replicates. Rows 3-5 is biological and row 6-8 is technical set) from 3 independent experiments. (B) Volcano plot shows the upregulated genes (blue boxes), the downregulated genes (red boxes), and the unchanged genes (gray boxes). The x-axis denotes log2 fold change whereas the Y-axis denotes the log₁₀ of p-values. (C) The pie chart shows the percentage of differentially regulated genes in each functional category as per the TubercuList. (D) The table shows the number of genes upregulated and downregulated in each functional category. The transcripts at 2.0-fold change cut-off were considered to be differentially regulated at p < 0.05. To avoid the false positive selection of transcripts, multiple testing correction was performed using Benjamini–Hochberg procedure in Gene Spring software. To check the statistical number of genes differentially regulated in each functional category, Hypergeometric probability (value in parentheses) was applied and bold values indicate the statistically significant at p < 0.05. (E) Validation of microarray data on qRT-PCR by analyzing the relative expression of 9 genes (5 upregulated and 4 downregulated) in two independent experiments (3-dpi), with technical replicates as compared to in vitro grown M. tb H37Rv. 16s rRNA was used as reference gene for normalization. Y-axis values (Log₂ fold change) of ≥1 indicate upregulation and values ≤ -1 indicate down-regulation. Each bar represents mean ± SD values for each of the genes with two biological and two technical replicates. *p < 0.05; **p < 0.01 by student's t-test.

The differentially expressed M. tuberculosis transcripts in RPE cells were functionally categorized and organized as per the TubercuList database (Lew et al., 2011). The highest percentage of transcripts that were differentially regulated encode for proteins corresponding to the categories of conserved hypothetical (25%), cell wall and cell process (22%), and intermediary metabolism and respiration (22%). Around 4–7% of the differentially regulated transcripts belonged to the PE and PPE families, regulatory proteins, virulence and detoxification, information pathways, or lipid metabolism. Remaining 1–3% transcripts belong to insertion sequences, phages, or stable RNAs functional categories (Figure 7C). Each functional group contained a different number of upregulated and downregulated transcripts (Figure 7D). The hypergeometric test (p < 0.05) was applied to determine the functional categories that were enriched in differentially expressed transcripts (Figure 7D) (Leong and Kipling, 2009; Qureshi and Sacan, 2013; Sharma et al., 2017). The annotated list of identified M. tuberculosis transcripts with a potential role in invasion, transcriptional mechanism, environmental adaptation, and survival inside RPE cells at 3-dpi has been provided in Table 2 according to their functional categories as sorted by the TubercuList database (http://tuberculist.epfl.ch).

TABLE 2

Table 2. Predicted/known functions of M. tuberculosis gene products whose transcripts were upregulated at time of replication of bacteria inside ARPE-19 cells, an in vitro model of intraocular tuberculosis.

Correlating Physiological State and the Transcriptional Signatures of M. tuberculosis in in vitro Model of IOTB

To further study the successful adaptation, survival, and replication of M. tuberculosis in the host RPE cells, we analyzed the involvement of altered transcripts in M. tuberculosis metabolic pathways. Identifying the genetic network of metabolism is important in understanding the physiology and pathogenesis of M. tuberculosis in a hostile host environment (Rohde et al., 2012; Warner, 2014; Du et al., 2016). Using BioCyc database (Karp et al., 2017) we identified the changes in metabolic pathways and defined the transcriptional signatures involved in the modulation of cellular functions such as biosynthesis, energy, central dogma, regulation, cellular process, cell exterior, and response to a stimulus. The complete SmartTable (BioCyc dashboard display) of cellular functions along with the defined transcripts with Log₂ fold change can be visualized using the following link (https://biocyc.org/dashboard/dashboard.html?st=biocyc13-29738-3741434995). Overall, it was observed that the intracellular M. tuberculosis in RPE cells upregulates the transcripts from the cellular process and cell exterior classes as their protein products are involved in host interaction, adhesion, cell cycle, and biogenesis/organization of cell wall, thus suggesting the adaptation and symbiosis of mycobacteria with the host RPE cells. The upregulation of transcripts encoding for the proteins involved in energy metabolism suggests the survival of bacteria inside the RPE cells. The upregulation of transcripts whose protein products are involved in biosynthesis of amino acid, fatty acid and lipids, carbohydrates, DNA, RNA, and proteins demonstrated the growth and replication (Supplementary Table 3). Thus, the transcriptional profile of M. tuberculosis during the replication stage after the adaption in host RPE cell demonstrates the alteration in transcripts controlling the physiological behavior of M. tuberculosis inside the RPE cells for its survival and intracellular replication, thus leading to productive infection.

Validation of M. tuberculosis Transcripts Identified in in vitro Model of Intraocular Tuberculosis

To validate the microarray data, we analyzed the expression of the top-most upregulated (Rv1230c) and downregulated (Rv2623) transcripts along with few other transcripts from the upregulated (Rv1971, Rv3872, and Rv2674) and downregulated (Rv3876, Rv1965, and Rv3875) categories by qRT-PCR. Quantitative RT-PCR showed similar upregulation and downregulation profile of the transcripts as detected on a microarray platform. Rv1230c, Rv1971, Rv3872, and Rv2674 were 2.5, 2.8, 7.4, and 6.5 Log₂ fold upregulated, respectively, while Rv2623, Rv3876, Rv1965, and Rv3875 were −6.2, −0.4, −4.2, and −7.6 Log₂ fold downregulated, respectively (Figure 7E). Altogether, the qRT-PCR results were in concordance with the microarray data.

Mycobacterium tuberculosis Transcripts in Vitreous Samples From Intraocular Tuberculosis Patients

Three M. tuberculosis transcripts showing upregulated expression in infected ARPE-19 cells were also tested qualitatively by RT-PCR for in vivo expression in human vitreous samples. Based on the C_T value (<35) and specific melt curves as explained in methodology, Rv1230c was identified in 6/7 confirmed human IOTB samples and in all the three possible IOTB cases; whereas, Rv1971 was detected in 6/7 confirmed and in 1/3 possible IOTB cases (Table 3). However, Rv3872 was not detected in any of the human IOTB samples. For all the three genes, H37Rv RNA showed C_T value <35 and specific product melt curves, whereas in case of all the negative control no C_T and transcript-specific melt curve was obtained.

TABLE 3

Table 3. Mycobacterium tuberculosis transcripts in vitreous samples from intraocular tuberculosis patients.

Discussion

An in vitro model of IOTB was established using the RPE cell line ARPE-19. Retinal pigment epithelial cells are the primary barrier cells in the eye to counter any invading pathogen and have been chosen as an experimental model in other diseases also, such as endophthalmitis caused by Bacillus cereus, to understand the host-pathogen interactions (Callegan et al., 2007). Supporting the earlier reported clinical and histopathological evidence for localization of AFB in RPE cells (Rao et al., 2006), the experimental studies with the ARPE-19 cell line in the present study also demonstrated that M. tuberculosis can utilize these cells for its survival. The scanning electron micrographs of the ARPE-19 cell line infected with M. tuberculosis showed the adherence of bacilli to these cells along with morphological cell surface changes (Figure 1). Similar morphological changes have been reported in the airway and alveolar epithelial cells after the adherence and uptake of M. tuberculosis (García-Pérez et al., 2003, 2008, 2012; Hall-Stoodley et al., 2006). These surface changes in non-professional phagocytic cells usually triggered by bacterial proteins (Hall-Stoodley et al., 2006; Ham et al., 2011; Vir et al., 2014) are considered as the mechanism evolved for the uptake and invasion of the pathogens such as M. tuberculosis, Salmonella typhimurium, and Listeria monocytogenes (Garcia-del Portillo and Finlay, 1994; Menon et al., 2003). Following the invasion, the intracellular localization of M. tuberculosis in the ARPE-19 cell line was monitored by TEM analysis (Figure 2) which showed M. tuberculosis bacilli surrounded by vacuolar membranes in the cytoplasm of ARPE-19 cells (Figure 2B). The presence of bacteria within the vacuoles in the RPE cell lines was in concordance with a similar study conducted in A549 alveolar epithelial cell line where multiple bacilli within the vacuoles were seen (Mehta et al., 1996). At 3-dpi, the RPE cells showed numerous vacuoles containing multiple bacilli, degenerated cytoplasm with disrupted morphology (Figure 2D) in comparison to 0-dpi (Figure 2B), suggesting that the increased mycobacterial load is cytotoxic to RPE cells (Figure 4B). The adherence and invasion along with multiplication of M. tuberculosis in ARPE-19 cells were also quantified using CFU enumeration which showed a significant percentage of M. tuberculosis adhering to and invading the RPE cells (Figure 3A and Table 1). Further, an increase in fold CFU of intracellular bacilli observed till 5-dpi (Figure 3B) confirmed the intracellular replication similar to that reported earlier in alveolar non-professional (Mehta et al., 1996; Dobos et al., 2000) and professional phagocytic cells (Sharma et al., 2012; Scharn et al., 2016). Thus, it is evident that the RPE cells, the non-professional phagocytic cells in the eye, provide a permissive habitat to virulent M. tuberculosis like other professional phagocytic cells (Russell et al., 2009). With the significant increase in intracellular growth, a loss in cellular viability of infected host cells was seen at 3-dpi and at 5-dpi (Figure 4B) in comparison to uninfected cells. Earlier, in A549 cells also, it has been reported that cytotoxicity at 3- and 5-dpi was related to the cytotoxic phenotype of increased mycobacterial growth (Mehta et al., 1996; Dobos et al., 2000). The loss in cellular viability of infected cells during M. tuberculosis replication in professional phagocytic cells is considered as a mechanism so that the dead cells are utilized as a bait for new neighboring cells for the dissemination of infection (Cambier et al., 2014; Mahamed et al., 2017). However, no such mechanism is known yet in the case of RPE cells.

The canonical pathway (endosomal-lysosomal trafficking) of intracellular fate of M. tuberculosis has already been elucidated via various in vitro models of professional (macrophages, neutrophils, monocytes, and dendritic cells) (López de Armentia et al., 2016) and non-professional phagocytic cells (epithelial, endothelial cells) (Baltierra-Uribe et al., 2014; López de Armentia et al., 2016). Mycobacterium tuberculosis in both types of phagocytic cells is known to block the final step of the maturation process involving fusion of phagosomes with lysosomes, thus inhibiting the formation of phagolysosome and resulting in its intracellular multiplication and survival (Seto et al., 2012). It has also been established that M. tuberculosis inhibits the fusion of lysosomes with phagosomes through the selective exclusion of the GTPase Rab7 and lysosomal-associated membrane protein 1 (LAMP-1) or LAMP-2 coupled with the retention of Rab5 on the phagosome (Via et al., 1997; López de Armentia et al., 2016). If these mechanisms are employed by tubercle bacilli for the intracellular survival in RPE cells also is still not known. The co-localization of PKH26-labeled M. tuberculosis with LAMP-1 (Figure 5) and LAMP-2 (Figure 6) during intracellular replication at 1- and 3-dpi was to lesser extent as compared to 0-dpi in ARPE-19 cell line. Decrease in colocalization with late endosomal markers during the multiplication stage of M. tuberculosis in the RPE cells suggests the deceleration of phagosomal maturation. Similarly, in human monocyte-derived dendritic cells and macrophages, tubercle bacilli were found translocating from phagolysosomes into the cytosol 2 days after infection (Van der Wel et al., 2007). This could be due to an attempt at multiplying M. tuberculosis for its cytosolic localization (Van der Wel et al., 2007), as bacteria in the cytosol are considered to have a superior capacity for overpowering the host cellular autophagy, a lysosomal degradation pathway for cytoplasmic materials (Jamwal et al., 2016). Intracellular fate of M. tuberculosis in non-professional phagocytic cells is dependent on more than one endocytic routes (Baltierra-Uribe et al., 2014). A detailed labeling protocol for early and late endosomes with various Rab family proteins will further help to elucidate the trafficking events of M. tuberculosis in RPE cells more appropriately.

Delineating the pathogen and/or host genetic factors is not only important for understanding the co-evolution of host-pathogen but also to have an insight in the capabilities of the pathogens in successively modulating the hostile host environment into a favorable niche (Reiling et al., 2018). Studying the transcriptional signatures of pathogens in an appropriate host is therefore known to aid in studying the disease pathogenesis (Ward et al., 2010). Therefore, we analyzed the mycobacterial transcriptional signatures in in vitro ARPE cell line model of IOTB in an attempt to understand the M. tuberculosis adaptation and survival in host RPE cells. The microarray technology was utilized as it is a well-established tool to study the whole-genome transcriptome leading to the availability of extensive information in the field of functional genomics (Nookaew et al., 2012). The identified transcripts were functionally categorized using the TubercuList database (Supplementary Table 2). The role of protein products of these transcripts in the physiological status of intracellularly replicating M. tuberculosis inside host environment was further described by studying the modulation in various metabolic pathways (through BioCyc database; Supplementary Table 3). Out of 1,040 differentially expressed transcripts by intracellularly replicating M. tuberculosis in ARPE-19 cells, we have primarily focused on those transcripts whose products (proteins) have known or hypothetical functions in the invasion, adaptation, survival, and replication in the host environment. The functional annotation of these transcripts is based on the studies carried out in other experimental conditions (in vivo or in vitro) of M. tuberculosis infection and such commonly expressed transcripts have been considered as hallmark transcriptional response of M. tuberculosis (Talaat et al., 2004; Rohde et al., 2012), thus supporting the validity of transcripts identified in the present study. In addition, many transcripts were uniquely altered in the in vitro model of IOTB (Supplementary Table 2) as these transcripts have not been reported earlier in other in vitro M. tuberculosis infection models (Jain et al., 2006; Rohde et al., 2012; Ryndak et al., 2015). It is important to understand the role of these uniquely identified transcripts in the pathogenesis of IOTB and could be exploited for new diagnostics/therapeutics for IOTB.

Bacterial pathogens are known to subvert the hostile host environment by regulating its transcription regulatory network which primarily includes the activation of two-component system and/or alternative sigma factors leading to modulation in the bacterial metabolism and physiology (Leistikow et al., 2010; Flentie et al., 2016). We found that Rv2028c transcript (from Rv2028-Rv2031 operon) was upregulated and its protein is known to regulate the different metabolic pathways and sub-cellular processes like the two-component system (DevS-DevR) under stress condition (Mushtaq et al., 2015). Another upregulated transcript of Rv0902c is also a part of two-component regulatory system (Rv0903c-Rv0902c) and is known to play a role in environmental adaptation and intracellular multiplication inside murine (Chandolia et al., 2014) and human (Graham and Clark-Curtiss, 1999; Haydel and Clark-Curtiss, 2004) macrophages and also controls mycobacterial metabolism and viability in professional phagocytic cells (Haydel et al., 2012). Additionally, the upregulation of 2 sigma factors transcripts of sigJ (Rv3328c) and sigM (Rv3911) and another possible sigma factor Rv1364c transcripts were observed during M. tuberculosis intracellular replication in ARPE-19 cells. Upregulation of ribosomal (Rv2056c, Rv3459c, Rv3456c, Rv1630, Rv0720, Rv0979A, Rv1642, and Rv0700), transcriptional (Rv0823c, Rv1994c, Rv3765c, Rv3160c, Rv1963c, Rv3416, Rv0232, Rv0602c, Rv1129c, Rv2324, Rv3183, and Rv3833), and translation initiation (Rv3462c) transcripts whose protein products have a role in RNA metabolism of mycobacteria (Supplementary Table 3) also suggested the active replicative state of M. tuberculosis inside the RPE cells. All these transcripts of sigma factors, ribosomal and translation regulators belonging to information and regulatory protein categories (Supplementary Table 2) were found to be uniquely regulated in the in vitro model of IOTB. However, among the transcriptional regulators, Rv0823c, Rv3416, and Rv1129c transcripts have earlier been defined in the primary macrophage model also (Rohde et al., 2012) and the Rv3833 transcript is known to be upregulated in another immune privileged in vitro central nervous system-TB model (Jain et al., 2006). Rv1129c has been identified as a transcription factor that is directly involved in the regulation of enzymes of methylcitrate cycle for cholesterol consumption during intracellular growth (Griffin et al., 2012) and Rv3833 is a member of AraC transcriptional regulator family which regulates the genes involved in carbon metabolism, stress responses, and virulence (Martin and Rosner, 2001; Frota et al., 2004).

Further, transcripts encoding proteins with a role in the adaptive process (Rv1026) or in stress (Rv3667) condition (Lew et al., 2011) were also upregulated in M. tuberculosis infected RPE cells. Rv1026 has been defined with a key role in maintaining the M. tuberculosis intracellular growth and persister-like features inside macrophages, thus leading to antibiotic tolerance (Chuang et al., 2015). Rv3667 is known to aid in reducing the metabolic stress experienced by bacterium inside the macrophages (Lee et al., 2013). To further maintain the balance between several stress and physiological conditions, a multifaceted integrated network known as the MprAB system is involved (He et al., 2006; Pang et al., 2007). This system regulates the diverse subnetwork of gene panel encoding for NADH dehydrogenase complex (Du et al., 2016). Five transcripts of probable NADH dehydrogenase complex (Rv3153: nuoI, Rv3147: nuoC, Rv3149: nuoE, Rv3154: nuoJ, and Rv3155: nuoK) involved in aerobic respiration were upregulated (Supplementary Table 3) during intracellular M. tuberculosis replication in ARPE-19 cells. Rv3147 and Rv3149 were uniquely regulated while Rv3153 and Rv3155 transcripts are also known to be upregulated in infected A549 cells (Ryndak et al., 2015), and Rv3154 is upregulated in both A549 (Ryndak et al., 2015) and macrophage models (Rohde et al., 2012). Besides, the ATPase family has a critical role in the response of mycobacteria against the toxic substances in the phagosomal environment during active and latent infection (Novoa-Aponte and Soto Ospina, 2014). Two transcripts from the ATPase family (CtpF: Rv1997 and CtpG:1992c) were found to be upregulated in the in vitro model of IOTB infection and Rv1992c transcript is reported to be upregulated in the macrophage model also (Rohde et al., 2012). Around 310 transcripts whose protein products constitute for cell wall component proteins, plasma membrane proteins, and cell wall biogenesis/organization were upregulated during the replication stage marking the growth of M. tuberculosis inside the host RPE cells (Supplementary Table 3). The MprAB system is also known to be involved in the activation of cellular processes associated with cell regrowth and in the regulation of cell wall component genes (He et al., 2006; Pang et al., 2007). However, the transcripts of MprAB was not found to be modulated in our study, which signifies that the change in these cell wall component protein is independent of the MprAB system and demands further analysis for regulation of these transcripts (Du et al., 2016). Certain transcripts whose proteins are involved in interaction with the host and performing symbiosis function during cellular processes were also found to be expressed (Supplementary Table 3). These results suggest that M. tuberculosis during its intracellular adaptation and growth is able to sustain the host-derived pressures by evolving its metabolic processes (Rohde et al., 2012).

Most genes upregulated from the cell wall and cell process and conserved hypothetical category are non-essential genes of M. tuberculosis during in vitro growth. It has been postulated that non-essential M. tuberculosis genes may have a potential role in in vivo survival and/or virulence (Jain et al., 2006). Rv1230c was among the top upregulated transcripts from the cell wall and cell processes functional category and is known to be involved in interacting with bacterial histone-like proteins which regulate many other genes involved in stress response and virulence factors (Gupta et al., 2014) and is also important for mycobacterial survival (Katsube et al., 2007). Presence of Rv1230c transcript (RNA) in 3 confirmed and 1 possible human IOTB cases (Table 3) suggests the presence of viable mycobacteria in vitreous samples of human IOTB cases. Interestingly, the protein encoded by Rv1230c has been reported as a potential drug target (Anand and Chandra, 2014). Two resuscitation-promoting factors (RpfC: Rv1884c and RpfE: Rv2450c) upregulated from the cell wall and cell processes functional category have earlier been reported to regulate the bacterial replication under stress conditions or modulate the host immune responses and contribute to mycobacterial pathogenesis (Rosser et al., 2017). Another upregulated gene Rv2235 from the same category is predicted to be involved in the survival of mycobacteria inside the human macrophages (Miller and Shinnick, 2001). Among differentially regulated transcripts from PE and PPE families, the majority of the proteins encoded by the upregulated transcripts have an unknown function. However, it is reported that genes of this category occurs only in mycobacteria and may hold a key to command the immune pathogenesis (Brennan, 2017). Additionally, due to the presence of these proteins (antigens) on the mycobacterial cell surface, they elicit good antibodies and T-cell response (Copin et al., 2014; Brennan, 2017). PPE10 (Rv0442c) transcript whose protein product is involved in multiple cellular functions, such as in interaction with host and symbiosis functions (Supplementary Table 3), was upregulated in RPE cells. It is also found to be upregulated in infected A549 cells (Ryndak et al., 2015) and is involved in the maintenance of mycobacterial capsule layer with a significant role in virulence and immune modulation (Ates et al., 2016).

During mycobacterial growth, many biosynthesis pathways are active which are known to be involved in the synthesis of amino acids, fatty acids, and lipids as well as folate biosynthesis. Thus, the upregulation of transcripts whose protein products are involved in these metabolic activities (Supplementary Table 3) highlights the active metabolic physiological state of bacteria inside the host environment during replication. Rv2672 encodes for a mycobacterial secretory protein with lipase and protease activities and hydrolyzes the host lipid for its nutrient requirements and helps in survival of M. tuberculosis inside lipid-rich foamy macrophages (Singh et al., 2017). Genes (Rv0462: lpdC and Rv1005c: pabB) of folate biosynthesis are predicted to be essential for mycobacterial survival (Griffin et al., 2011) and enzymes of folate metabolism are targets of potent antitubercular agents (Minato et al., 2015). Rv0462 transcript is also known to be upregulated in M. tuberculosis infected A549 cells (Ryndak et al., 2015) and is involved in biosynthesis, energy, cellular process, and cell exterior pathways (Supplementary Table 3). Rv0462 contributes to induction of host innate and adaptive immunity (Heo et al., 2011) and plays a direct role in restricting the phagosome maturation (Philips, 2008). In parallel, the upregulation of many transcripts (Rv2238c, Rv2874, Rv0462, Rv1676, and Rv3673c) whose protein products have a role in oxidant detoxifications suggests the removal of toxic byproducts from fatty acid metabolism or another metabolic pathway (Supplementary Table 2) generated during active mycobacterial growth. It was observed that Rv2238c, Rv1676, and Rv3673c were uniquely upregulated in the in vitro model of IOTB while Rv2874 is also known to be upregulated in the macrophage model (Rohde et al., 2012). Overall, 25 M. tuberculosis genes were upregulated from the functional category of virulence and detoxification (Figure 7D). Among these genes, four genes (Rv1971, Rv0169, Rv0592, and Rv1967) belong to mammalian cell entry (mce) family protein that plays a crucial role in invasion and virulence of M. tuberculosis (Zhang and Xie, 2011). Recently, Rv1971 from this family was also reported to be upregulated in sputum samples of pulmonary TB patients (Sharma et al., 2017) and bioinformatically observed to be a highly antigenic protein involved in immune responses (Zhang et al., 2009). Interestingly, this gene was also detected in human IOTB samples (Table 3). Another mce gene Rv0169 has been shown to be involved in the entry and survival of M. tuberculosis inside the nonphagocytic HeLa cells (Arruda et al., 1993; Saini et al., 2008), in cytoskeletal rearrangements (Chitale et al., 2001), and is also predicted to play a role in nutrient deficient environment inside the host cells (Pandey and Sassetti, 2008; Saini et al., 2008). Besides, the protein encoded by Rv0169 is located on mycobacterial cell wall with its cell entry epitope flank on the surface (Das et al., 2003; Mitra et al., 2005) suggesting its role in host-pathogen interactions (Shimono et al., 2003). Rv1755c upregulated during intracellular replication of mycobacteria has predicted protein product with virulence function and is involved in the pathogenesis of M. tuberculosis for intracellular survival by altering the cell signaling events or directly causing cytotoxicity to host cells (Castro-Garza et al., 2016). Altogether, the upregulation of these transcripts controlling the various metabolic states of mycobacteria reflects the adaptation of M. tuberculosis to host immunity and environmental stresses in different conditions.

Besides, most of the significantly downregulated (through hypergeometric statistical evaluation with p < 0.05, Figure 7D) genes were from the functional categories of conserved hypothetical (135 genes), cell wall and cell process (133 genes), intermediary metabolism and respiration (133 genes), and lipid metabolism (51 genes), an indicator of decrease in metabolic state of M. tuberculosis soon after infection and internalization in host cells, which has been attributed to precursor of latency (Jain et al., 2006). Thus, analysis of the transcriptional profile of M. tuberculosis in RPE cells demonstrates the modulation of several transcripts whose protein products have earlier been reported to have a role in intracellular adaptation, survival, and multiplication in different in vitro experimental models of professional and non-professional phagocytic cells. These results not only favor the role of RPE cells as potential host cells for the survival of mycobacteria in the ocular environment but also suggest the role of the identified transcripts as potential targets for anti-tuberculous drugs and could also be exploited in IOTB detection. These findings indicate that M. tuberculosis growth is a programmed regulation of various metabolic pathways genes which helps in the survival and replication of bacteria inside RPE cells. Although the microarray technology used in the present study lead to the identification of a large number of differentially regulated M. tuberculosis transcripts in in vitro model of IOTB, it may not have detected all the differentially regulated transcripts as pre-designed complement probes against each gene (of M. tuberculosis genome) were used. The advancement in transcriptomic studies with technologies like RNA sequencing may add on to the current study.

In conclusion, this study demonstrates the invasion and replication of virulent strain of M. tuberculosis (H37Rv) in RPE cells, and the mycobacterial transcripts identified in the intracellular environment of RPE cells correspond to the pathways which have a crucial role in invasion, host adaptation, replication, and survival. As many of these transcript products are regulators of important metabolic pathways during M. tuberculosis growth and survival inside RPE cells, the results of this study may be important for unraveling the many unknown facets of IOTB pathogenesis.

Author Contributions

SA, IV, and AG conceptualized and designed the study. SA performed the experiments. SA and IV carried out data analysis, manuscript writing and editing. VG, RB and NS helped with studying human samples. US supported the TEM studies. SL arranged for resources and the training of SA for this study under the NIH funded Fogarty International Centre, USA (#D43 TW009588). AG and IV arranged for funding and resources.

Funding

This work was supported by Department of Biotechnology, Ministry of Science & Technology, Government of India [Grant #BT/01/CEIB/11/I/02 to IV].

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Dr. C. S. Rayyat for taking TEM images, Department of Histopathology, PGIMER, Chandigarh, India.

Mr. Virender Singh Mahta, Senior technologist, for taking TEM images, Department of CSIC, PGIMER, Chandigarh, India.

Mr. Gurpreet Singh, Senior technologist, for taking Confocal images, Department of Translational and Regenerative Medicine, PGIMER, Chandigarh, India.

Dr. Michelle Ryndak, Department of Pathology, New York University Langone Medical Center, New York, NY-10010, United States of America for training in microarray technology funded by National Institute of Health, Fogarty International Centre, USA (Grant # D43 TW009588 to SL).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2018.00330/full#supplementary-material

Supplementary Table 1. Mycobacterium tuberculosis genes specific primers used for a real time-polymerase chain reaction.

Supplementary Table 2. Transcriptional signatures of M. tuberculosis in an in vitro model of IOTB. M. tuberculosis transcripts differentially regulated in each functional category in an in vitro model of IOTB at 3-day post infection.

Supplementary Table 3. Mycobacterium tuberculosis proteins of biosynthesis, central dogma, energy, cellular process, and cell exterior pathways whose transcripts were upregulated (>±1.15 Log2 FC) during intracellular multiplication of bacteria at 3-days post-infection in in vitro model of intraocular tuberculosis.

Supplementary Figure 1. Hierarchical clustering and principal components analysis of microarray data. (A) The hierarchical clustering of all the samples (infected and control) shows the replicability of the data obtained from each set. Six infected sets (3 biological and 3 technical replicate), isolated from infected ARPE-19 cells at 3-dpi and 2 control sets (1 biological and 1 technical) were evaluated for intracellular M. tuberculosis transcriptional signatures. (B) Principal components analysis of each sample (infected and control) was analyzed for correlation coefficient. The correlation coefficient data demonstrated the correlation between each respective set. The control sample (H37Rv RNA) was similar to its technical set, while all the infected sets (intracellular H37Rv RNA) showed similar correlation and thus were treated as replicates. (C) Values of co-relation coefficient obtained for each identifier (control and infected sample set). Dpi, days post infection.

SG14484452_252632310045_1_1(Experiment-1),

SG14484452_252632310045_1_2(Experiment-2),

SG14484452_252632310045_1_3(Experiment-3),

SG14484452_252632310045_1_4(Control-1),

SG14484452_252632310045_2_1(technical repeat-1),

SG14484452_252632310045_2_2(technical repeat-2),

SG14484452_252632310045_2_3(technical repeat-3), and

SG14484452_252632310045_2_4 (control technical repeat-1).

References

Agrawal, R., Gunasekeran, D.V., Grant, R., Agarwal, A., Kon, O.M., Nguyen, Q.D., et al. (2017). Clinical features and outcomes of patients with tubercular uveitis treated with antitubercular therapy in the Collaborative Ocular Tuberculosis Study (COTS)−1. JAMA Ophthalmol. 135:1318. doi: 10.1001/jamaophthalmol.2017.4485

PubMed Abstract | CrossRef Full Text | Google Scholar

Aguilar-Ayala, D. A., Tilleman, L., Van Nieuwerburgh, F., Deforce, D., Palomino, J. C., Vandamme, P., et al. (2017). The transcriptome of Mycobacterium tuberculosis in a lipid-rich dormancy model through RNAseq analysis. Sci. Rep. 7:17665. doi: 10.1038/s41598-017-17751-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Anand, P., and Chandra, N. (2014). Characterizing the pocketome of Mycobacterium tuberculosis and application in rationalizing polypharmacological target selection. Sci. Rep. 3:6356. doi: 10.1038/srep06356

CrossRef Full Text | Google Scholar

Arruda, S., Bomfim, G., Knights, R., Huima-Byron, T., and Riley, L. W. (1993). Cloning of an M. tuberculosis DNA fragment associated with entry and survival inside cells. Science 261, 1454–1457. doi: 10.1126/science.8367727

PubMed Abstract | CrossRef Full Text | Google Scholar

Ates, L. S., van der Woude, A. D., Bestebroer, J., van Stempvoort, G., Musters, R. J. P., Garcia-Vallejo, J. J., et al. (2016). The ESX-5 system of pathogenic mycobacteria is involved in capsule integrity and virulence through its substrate PPE10. PLOS Pathog. 12:e1005696. doi: 10.1371/journal.ppat.1005696

PubMed Abstract | CrossRef Full Text | Google Scholar

Balasubramanian, V., Wiegeshaus, E. H., Taylor, B. T., and Smith, D. W. (1994). Pathogenesis of tuberculosis: pathway to apical localization. Tuber. Lung Dis. 75, 168–178. doi: 10.1016/0962-8479(94)90002-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Baltierra-Uribe, S. L., García-Vásquez, M. J., Castrejón-Jiménez, N. S., Estrella-Piñón, M. P., Luna-Herrera, J., and García-Pérez, B. E. (2014). Mycobacteria entry and trafficking into endothelial cells. Can. J. Microbiol. 60, 569–577. doi: 10.1139/cjm-2014-0087

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrios-Payán, J., Saqui-Salces, M., Jeyanathan, M., Alcántara-Vazquez, A., Castañon-Arreola, M., Rook, G., et al. (2012). Extrapulmonary locations of Mycobacterium tuberculosis DNA during latent infection. J. Infect. Dis. 206, 1194–1205. doi: 10.1093/infdis/jis381

PubMed Abstract | CrossRef Full Text | Google Scholar

Be, N. A., Lamichhane, G., Grosset, J., Tyagi, S., Cheng, Q., Kim, K. S., et al. (2008). Murine model to study the invasion and survival of Mycobacterium tuberculosis in the central nervous system. J. Infect. Dis. 198, 1520–1528. doi: 10.1086/592447

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300.

Google Scholar

Benjamini, Y., and Hochberg, Y. (2000). On the adaptive control of the false discovery rate in multiple testing with independent statistics. J. Educ. Behav. Stat. 25, 60–83. doi: 10.3102/10769986025001060

CrossRef Full Text | Google Scholar

Brennan, M. J. (2017). The enigmatic PE/PPE multigene family of mycobacteria and tuberculosis vaccination. Infect. Immun. 85:e00969-16. doi: 10.1128/IAI.00969-16.

PubMed Abstract | CrossRef Full Text | Google Scholar

Callegan, M. C., Gilmore, M. S., Gregory, M., Ramadan, R. T., Wiskur, B. J., Moyer, A. L., et al. (2007). Bacterial endophthalmitis: Therapeutic challenges and host-pathogen interactions. Prog. Retin. Eye Res. 26, 189–203. doi: 10.1016/j.preteyeres.2006.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Cambier, C. J., Falkow, S., and Ramakrishnan, L. (2014). Host evasion and exploitation schemes of Mycobacterium tuberculosis. Cell 159, 1497–1509. doi: 10.1016/j.cell.2014.11.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Castro-Garza, J., González-Salazar, F., Quinn, F. D., Karls, R. K., De La Garza-Salinas, L. H., Guzmán-de la Garza, F. J., et al. (2016). An acidic sphingomyelinase Type C activity from Mycobacterium tuberculosis. Rev. Argent. Microbiol. 48, 21–26. doi: 10.1016/j.ram.2016.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandolia, A., Rathor, N., Sharma, M., Saini, N. K., Sinha, R., Malhotra, P., et al. (2014). Functional analysis of mce4A gene of Mycobacterium tuberculosis H37Rv using antisense approach. Microbiol. Res. 169, 780–787. doi: 10.1016/j.micres.2013.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitale, S., Ehrt, S., Kawamura, I., Fujimura, T., Shimono, N., Anand, N., et al. (2001). Recombinant Mycobacterium tuberculosis protein associated with mammalian cell entry. Cell. Microbiol. 3, 247–254. doi: 10.1046/j.1462-5822.2001.00110.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chuang, Y. M., Bandyopadhyay, N., Rifat, D., Rubin, H., Bader, J. S., and Karakousis, P. C. (2015). Deficiency of the novel exopolyphosphatase Rv1026/PPX2 leads to metabolic downshift and altered cell wall permeability in Mycobacterium tuberculosis. MBio 6:e02428-14. doi: 10.1128/mBio.02428-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Copin, R., Coscollá, M., Seiffert, S. N., Bothamley, G., Sutherland, J., Mbayo, G., et al. (2014). Sequence diversity in the pe_pgrs genes of Mycobacterium tuberculosis is independent of human T cell recognition. MBio 5:e00960-13. doi: 10.1128/mBio.00960-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Dalvin, L. A., and Smith, W. M. (2017). Intraocular manifestations of Mycobacterium tuberculosis: a review of the literature. J. Clin. Tuberc. Other Mycobact. Dis. 7, 13–21. doi: 10.1016/j.jctube.2017.01.003

CrossRef Full Text | Google Scholar

Danelishvili, L., McGarvey, J., Li, Y.-J., and Bermudez, L. E. (2003). Mycobacterium tuberculosis infection causes different levels of apoptosis and necrosis in human macrophages and alveolar epithelial cells. Cell. Microbiol. 5, 649–660. doi: 10.1046/j.1462-5822.2003.00312.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, A. K., Mitra, D., Harboe, M., Nandi, B., Harkness, R. E., Das, D., et al. (2003). Predicted molecular structure of the mammalian cell entry protein Mce1A of Mycobacterium tuberculosis. Biochem. Biophys. Res. Commun. 302, 442–447. doi: 10.1016/S0006-291X(03)00116-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobos, K. M., Spotts, E. A., Quinn, F. D., and King, C. H. (2000). Necrosis of lung epithelial cells during infection with Mycobacterium tuberculosis is preceded by cell permeation. Infect. Immun. 68, 6300–6310. doi: 10.1128/IAI.68.11.6300-6310.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Doycheva, D., Pfannenberg, C., Hetzel, J., Deuter, C. M. E., Pavesio, C., Kempf, V. A. J., et al. (2010). Presumed tuberculosis-induced retinal vasculitis, diagnosed with positron emission tomography (18F-FDG-PET/CT), aspir- ation biopsy, and culture. Ocul. Immunol. Inflamm. 18, 194–199. doi: 10.3109/09273948.2010.483318

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, P., Sohaskey, C. D., and Shi, L. (2016). Transcriptional and physiological changes during Mycobacterium tuberculosis reactivation from non-replicating persistence. Front. Microbiol. 7:1346. doi: 10.3389/fmicb.2016.01346

PubMed Abstract | CrossRef Full Text | Google Scholar

Flentie, K., Garner, A. L., and Stallings, C. L. (2016). Mycobacterium tuberculosis transcription machinery: ready to respond to host attacks. J. Bacteriol. 198, 1360–1373. doi: 10.1128/JB.00935-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Fonseca, K. L., Rodrigues, P. N. S., Olsson, I. A. S., and Saraiva, M. (2017). Experimental study of tuberculosis: From animal models to complex cell systems and organoids. PLoS Pathog. 13:e1006421. doi: 10.1371/journal.ppat.1006421

PubMed Abstract | CrossRef Full Text | Google Scholar

Frota, C. C., Papavinasasundaram, K. G., Davis, E. O., and Colston, M. J. (2004). The AraC family transcriptional regulator Rv1931c plays a role in the virulence of Mycobacterium tuberculosis. Infect. Immun. 72, 5483–5486. doi: 10.1128/IAI.72.9.5483-5486.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia-del Portillo, F., and Finlay, B. B. (1994). Salmonella invasion of nonphagocytic cells induces formation of macropinosomes in the host cell. Infect. Immun. 62, 4641–4645.

PubMed Abstract | Google Scholar

García-Pérez, B. E., De la Cruz-López, J. J., Castañeda-Sánchez, J. I., Muñóz-Duarte, A. R., Hernández-Pérez, A. D., Villegas-Castrejón, H., et al. (2012). Macropinocytosis is responsible for the uptake of pathogenic and non-pathogenic mycobacteria by B lymphocytes (Raji cells). BMC Microbiol. 12:246. doi: 10.1186/1471-2180-12-246

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Pérez, B. E., Hernández-González, J. C., García-Nieto, S., and Luna-Herrera, J. (2008). Internalization of a non-pathogenic mycobacteria by macropinocytosis in human alveolar epithelial A549 cells. Microb. Pathog. 45, 1–6. doi: 10.1016/j.micpath.2008.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Pérez, B. E., Mondragón-Flores, R., and Luna-Herrera, J. (2003). Internalization of Mycobacterium tuberculosis by macropinocytosis in non-phagocytic cells. Microb. Pathog. 35, 49–55. doi: 10.1016/S0882-4010(03)00089-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Gogia, V., Venkatesh, P., Garg, S. P., Takkar, B., and Sheemar, A. (2018). Patterns of uveitis in patients with proven systemic (pulmonary and extrapulmonary) tuberculosis. Int. Ophthalmol. doi: 10.1007/s10792-018-0989-9. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, J. E., and Clark-Curtiss, J. E. (1999). Identification of Mycobacterium tuberculosis RNAs synthesized in response to phagocytosis by human macrophages by selective capture of transcribed sequences (SCOTS). Proc. Natl. Acad. Sci. U.S.A. 96, 11554–11559. doi: 10.1073/pnas.96.20.11554

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffin, J. E., Gawronski, J. D., DeJesus, M. A., Ioerger, T. R., Akerley, B. J., and Sassetti, C. M. (2011). High-resolution phenotypic profiling defines genes essential for mycobacterial growth and cholesterol catabolism. PLoS Pathog. 7:e1002251. doi: 10.1371/journal.ppat.1002251

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffin, J. E., Pandey, A. K., Gilmore, S. A., Mizrahi, V., Mckinney, J. D., Bertozzi, C. R., et al. (2012). Cholesterol catabolism by Mycobacterium tuberculosis requires transcriptional and metabolic adaptations. Chem. Biol. 19, 218–227. doi: 10.1016/j.chembiol.2011.12.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, A., Sharma, A., Bansal, R., and Sharma, K. (2015). Classification of intraocular tuberculosis. Ocul. Immunol. Inflamm. 23, 7–13. doi: 10.3109/09273948.2014.967358

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, B., Agrawal, R., Swampillai, A. J., Lim, R. H., Kee, A., Gunasekaran, D., et al. (2016). Ocular manifestations of tuberculosis: an update. Expert Rev. Ophthalmol. 11, 145–154. doi: 10.1586/17469899.2016.1152887

CrossRef Full Text | Google Scholar

Gupta, M., Sajid, A., Sharma, K., Ghosh, S., Arora, G., Singh, R., et al. (2014). HupB, a nucleoid-associated protein of Mycobacterium tuberculosis, is modified by serine/threonine protein kinases in vivo. J. Bacteriol. 196, 2646–2657. doi: 10.1128/JB.01625-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, V., Gupta, A., and Rao, N. A. (2007). Intraocular tuberculosis-an update. Surv. Ophthalmol. 52, 561–587. doi: 10.1016/j.survophthal.2007.08.015

CrossRef Full Text | Google Scholar

Hall-Stoodley, L., Watts, G., Crowther, J. E., Balagopal, A., Torrelles, J. B., Robison-Cox, J., et al. (2006). Mycobactenum tuberculosis binding to human surfactant proteins A and D, fibronectin, and small airway epithelial cells under shear conditions. Infect. Immun. 74, 3587–3596. doi: 10.1128/IAI.01644-05

PubMed Abstract | CrossRef Full Text | Google Scholar

Ham, H., Sreelatha, A., and Orth, K. (2011). Manipulation of host membranes by bacterial effectors. Nat. Rev. Microbiol. 9, 635–646. doi: 10.1038/nrmicro2602

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasan, Z., Schlax, C., Kuhn, L., Lefkovits, I., Young, D., Thole, J., et al. (1997). Isolation and characterization of the mycobacterial phagosome: segregation from the endosomal/lysosomal pathway. Mol. Microbiol. 24, 545–553.

PubMed Abstract | Google Scholar

Haydel, S. E., and Clark-Curtiss, J. E. (2004). Global expression analysis of two-component system regulator genes during Mycobacterium tuberculosis growth in human macrophages. FEMS Microbiol. Lett. 236, 341–347. doi: 10.1111/j.1574-6968.2004.tb09667.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haydel, S. E., Malhotra, V., Cornelison, G. L., and Clark-Curtiss, J. E. (2012). The prrAB two-component system is essential for Mycobacterium tuberculosis viability and is induced under nitrogen-limiting conditions. J. Bacteriol. 194, 354–361. doi: 10.1128/JB.06258-11

PubMed Abstract | CrossRef Full Text | Google Scholar

He, H., Hovey, R., Kane, J., Singh, V., and Zahrt, T. C. (2006). MprAB is a stress-responsive two-component system that directly regulates expression of sigma factors SigB and SigE in Mycobacterium tuberculosis. J. Bacteriol. 188, 2134–2143. doi: 10.1128/JB.188.6.2134-2143.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Henderson, H. J., Dannenberg, A. M., and Lurie, M. B. (1963). Phagocytosis of tubercle bacilli by rabbit pulmonary alveolar macrophages and its relation to native resistance to tuberculosis. J. Immunol. 91, 553–556.

PubMed Abstract | Google Scholar

Heo, D. R., Shin, S. J., Kim, W. S., Noh, K. T., Park, J. W., Son, K. H., et al. (2011). Mycobacterium tuberculosis lpdC, Rv0462, induces dendritic cell maturation and Th1 polarization. Biochem. Biophys. Res. Commun. 411, 642–647. doi: 10.1016/j.bbrc.2011.07.013

CrossRef Full Text | Google Scholar

Huynh, K. K., Eskelinen, E.-L., Scott, C. C., Malevanets, A., Saftig, P., and Grinstein, S. (2007). LAMP proteins are required for fusion of lysosomes with phagosomes. EMBO J. 26, 313–324. doi: 10.1038/sj.emboj.7601511

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, S. K., Paul-Satyaseela, M., Lamichhane, G., Kim, K. S., and Bishai, W. R. (2006). Mycobacterium tuberculosis invasion and traversal across an in vitro human blood-brain barrier as a pathogenic mechanism for central nervous system tuberculosis. J. Infect. Dis. 193, 1287–1295. doi: 10.1086/502631

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamwal, S. V., Mehrotra, P., Singh, A., Siddiqui, Z., Basu, A., and Rao, K. V. S. (2016). Mycobacterial escape from macrophage phagosomes to the cytoplasm represents an alternate adaptation mechanism. Sci. Rep. 6:23089. doi: 10.1038/srep23089

PubMed Abstract | CrossRef Full Text | Google Scholar

Karp, P. D., Billington, R., Caspi, R., Fulcher, C. A., Latendresse, M., Kothari, A., et al. (2017). The BioCyc collection of microbial genomes and metabolic pathways. Brief. Bioinform. doi: 10.1093/bib/bbx085. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Katsube, T., Matsumoto, S., Takatsuka, M., Okuyama, M., Ozeki, Y., Naito, M., et al. (2007). Control of cell wall assembly by a histone-like protein in mycobacteria. J. Bacteriol. 189, 8241–8249. doi: 10.1128/JB.00550-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Keselman, H. J., Cribbie, R., and Holland, B. (2002). Controlling the rate of Type I error over a large set of statistical tests. Br. J. Math. Stat. Psychol. 55, 27–39. doi: 10.1348/000711002159680

PubMed Abstract | CrossRef Full Text | Google Scholar

Krishnan, N., Robertson, B. D., and Thwaites, G. (2010). The mechanisms and consequences of the extra-pulmonary dissemination of Mycobacterium tuberculosis. Tuberculosis 90, 361–366. doi: 10.1016/j.tube.2010.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Laurell, H., Iacovoni, J. S., Abot, A., Svec, D., Maoret, J. J., Arnal, J. F., et al. (2012). Correction of RT-qPCR data for genomic DNA-derived signals with ValidPrime. Nucleic Acids Res. 40:e51. doi: 10.1093/nar/gkr1259

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C., Agrawal, R., and Pavesio, C. (2016). Ocular tuberculosis - a clinical conundrum. Ocul. Immunol. Inflamm. 24, 237–242. doi: 10.3109/09273948.2014.985387

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, W., VanderVen, B. C., Fahey, R. J., and Russell, D. G. (2013). Intracellular Mycobacterium tuberculosis exploits host-derived fatty acids to limit metabolic stress. J. Biol. Chem. 288, 6788–6800. doi: 10.1074/jbc.M112.445056

PubMed Abstract | CrossRef Full Text | Google Scholar

Leistikow, R. L., Morton, R. A., Bartek, I. L., Frimpong, I., Wagner, K., and Voskuil, M. I. (2010). The Mycobacterium tuberculosis DosR regulon assists in metabolic homeostasis and enables rapid recovery from nonrespiring dormancy. J. Bacteriol. 192, 1662–1670. doi: 10.1128/JB.00926-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Leong, H. S., and Kipling, D. (2009). Text-based over-representation analysis of microarray gene lists with annotation bias. Nucleic Acids Res. 37:e79. doi: 10.1093/nar/gkp310

PubMed Abstract | CrossRef Full Text | Google Scholar

Lew, J. M., Kapopoulou, A., Jones, L. M., and Cole, S. T. (2011). TubercuList−10 years after. Tuberculosis (Edinb) 91, 1–7. doi: 10.1016/j.tube.2010.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W., De Sessions, P. F., Teoh, G. H. K., Mohamed, A. N. N., Zhu, Y. O., Koh, V. H. Q., et al. (2016). Transcriptional profiling of Mycobacterium tuberculosis exposed to in vitro lysosomal stress. Infect. Immun. 84, 2505–2523. doi: 10.1128/IAI.00072-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

López de Armentia, M., Amaya, C., and Colombo, M. (2016). Rab GTPases and the autophagy pathway: bacterial targets for a suitable biogenesis and trafficking of their own vacuoles. Cells 5:11. doi: 10.3390/cells5010011

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahamed, D., Boulle, M., Ganga, Y., McArthur, C., Skroch, S., Oom, L., et al. (2017). Intracellular growth of Mycobacterium tuberculosis after macrophage cell death leads to serial killing of host cells. Elife 6:e22028. doi: 10.7554/eLife.22028.

CrossRef Full Text | Google Scholar

Martin, R. G., and Rosner, J. L. (2001). The AraC transcriptional activators. Curr. Opin. Microbiol. 4, 132–137. doi: 10.1016/S1369-5274(00)00178-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Mehta, P. K., King, C. H., White, E. H., Murtagh, J. J., and Quinn, F. D. (1996). Comparison of in vitro models for the study of Mycobacterium tuberculosis invasion and intracellular replication. Infect. Immun. 64, 2673–2679.

Google Scholar

Menon, A., Shroyer, M. L., Wampler, J. L., Chawan, C. B., and Bhunia, A. K. (2003). In vitro study of Listeria monocytogenes infection to murine primary and human transformed B cells. Comp. Immunol. Microbiol. Infect. Dis. 26, 157–174. doi: 10.1016/S0147-9571(02)00039-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, B. H., and Shinnick, T. M. (2001). Identification of two Mycobacterium tuberculosis H37Rv ORFs involved in resistance to killing by human macrophages. BMC Microbiol. 1:26. doi: 10.1186/1471-2180-1-26

PubMed Abstract | CrossRef Full Text | Google Scholar

Minato, Y., Thiede, J. M., Kordus, S. L., McKlveen, E. J., Turman, B. J., and Baughn, A. D. (2015). Mycobacterium tuberculosis folate metabolism and the mechanistic basis for para-aminosalicylic acid susceptibility and resistance. Antimicrob. Agents Chemother. 59, 5097–5106. doi: 10.1128/AAC.00647-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitra, D., Saha, B., Das, D., Wiker, H. G., and Das, A. K. (2005). Correlating sequential homology of Mce1A, Mce2A, Mce3A and Mce4A with their possible functions in mammalian cell entry of Mycobacterium tuberculosis performing homology modeling. Tuberculosis 85, 337–345. doi: 10.1016/j.tube.2005.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Moharana, B., Bansal, R., Singh, R., Sharma, A., Gupta, V., and Gupta, A. (2018). Enhanced depth imaging by high-resolution spectral domain optical coherence tomography in tubercular multifocal serpiginoid choroiditis. Ocul. Immunol. Inflamm. doi: 10.1080/09273948.2018.1465101. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Mushtaq, K., Sheikh, J. A., Amir, M., Khan, N., Singh, B., and Agrewala, J. N. (2015). Rv2031c of Mycobacterium tuberculosis: a master regulator of Rv2028-Rv2031 (HspX) operon. Front. Microbiol. 6:351. doi: 10.3389/fmicb.2015.00351

PubMed Abstract | CrossRef Full Text | Google Scholar

Nandy, A., Mondal, A. K., Pandey, R., Arumugam, P., Dawa, S., Jaisinghani, N., et al. (2018). Adipocyte model of Mycobacterium tuberculosis infection reveals differential availability of iron to bacilli in the lipid-rich caseous environment. Infect. Immun. 86, e00041–e00018. doi: 10.1128/IAI.00041-18

CrossRef Full Text | Google Scholar

Napier, R. J., Rafi, W., Cheruvu, M., Powell, K. R., Zaunbrecher, M. A., Bornmann, W., et al. (2011). Imatinib-Sensitive tyrosine kinases regulate mycobacterial pathogenesis and represent therapeutic targets against tuberculosis. Cell Host Microbe 10, 475–485. doi: 10.1016/j.chom.2011.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazari, H., Karakousis, P. C., and Rao, N. a (2014). Replication of Mycobacterium tuberculosis in retinal pigment epithelium. JAMA Ophthalmol. 132, 724–729. doi: 10.1001/jamaophthalmol.2014.270

PubMed Abstract | CrossRef Full Text | Google Scholar

Nookaew, I., Papini, M., Pornputtapong, N., Scalcinati, G., Fagerberg, L., Uhlén, M., et al. (2012). A comprehensive comparison of RNA-Seq-based transcriptome analysis from reads to differential gene expression and cross-comparison with microarrays: a case study in Saccharomyces cerevisiae. Nucleic Acids Res. 40, 10084–10097. doi: 10.1093/nar/gks804

PubMed Abstract | CrossRef Full Text | Google Scholar

Novoa-Aponte, L., and Soto Ospina, C. Y. (2014). Mycobacterium tuberculosis P-type ATPases: possible targets for drug or vaccine development. Biomed Res. Int. 2014:296986. doi: 10.1155/2014/296986

CrossRef Full Text | Google Scholar

Pandey, A. K., and Sassetti, C. M. (2008). Mycobacterial persistence requires the utilization of host cholesterol. Proc. Natl. Acad. Sci. U.S.A. 105, 4376–4380. doi: 10.1073/pnas.0711159105

CrossRef Full Text | Google Scholar

Pang, X., Vu, P., Byrd, T. F., Ghanny, S., Soteropoulos, P., Mukamolova, G. V., et al. (2007). Evidence for complex interactions of stress-associated regulons in an mprAB deletion mutant of Mycobacterium tuberculosis. Microbiology 153, 1229–1242. doi: 10.1099/mic.0.29281-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Patel, S. S., Saraiya, N. V., Tessler, H. H., and Goldstein, D. A. (2013). Mycobacterial ocular inflammation: delay in diagnosis and other factors impacting morbidity. JAMA Ophthalmol. 131, 752–758. doi: 10.1001/jamaophthalmol.2013.71

PubMed Abstract | CrossRef Full Text | Google Scholar

Pawitan, Y., Michiels, S., Koscielny, S., Gusnanto, A., and Ploner, A. (2005). False discovery rate, sensitivity and sample size for microarray studies. Bioinformatics 21, 3017–3024. doi: 10.1093/bioinformatics/bti448

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, C., Wang, P., Xu, X., Wang, X., Wei, W., Chen, X., et al. (2016). Development of a qualitative real-time PCR method to detect 19 targets for identification of genetically modified organisms. Springerplus 5:889. doi: 10.1186/s40064-016-2395-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Philips, J. A. (2008). Mycobacterial manipulation of vacuolar sorting. Cell. Microbiol. 10, 2408–2415. doi: 10.1111/j.1462-5822.2008.01239.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Qureshi, R., and Sacan, A. (2013). Weighted set enrichment of gene expression data. BMC Syst. Biol. 7:S10. doi: 10.1186/1752-0509-7-S4-S10

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramsugit, S., and Pillay, M. (2016). Identification of Mycobacterium tuberculosis adherence-mediating components: a review of key methods to confirm adhesin function. Iran. J. Basic Med. Sci. 19, 579–584.

PubMed Abstract | Google Scholar

Rao, N. A., Saraswathy, S., and Smith, R. E. (2006). Tuberculous uveitis: Distribution of Mycobacterium tuberculosis in the retinal pigment epithelium. Arch. Ophthalmol. 124, 1777–1779. doi: 10.1001/archopht.124.12.1777

PubMed Abstract | CrossRef Full Text | Google Scholar

Reiling, N., Homolka, S., Kohl, T. A., Steinhäuser, C., Kolbe, K., Schütze, S., et al. (2018). Shaping the niche in macrophages: Genetic diversity of the M. tuberculosis complex and its consequences for the infected host. Int. J. Med. Microbiol. 308, 118–128. doi: 10.1016/j.ijmm.2017.09.009

CrossRef Full Text | Google Scholar

Rohde, K. H., Veiga, D. F. T., Caldwell, S., Balázsi, G., and Russell, D. G. (2012). Linking the transcriptional profiles and the physiological states of Mycobacterium tuberculosis during an extended intracellular infection. PLoS Pathog. 8:e1002769. doi: 10.1371/journal.ppat.1002769

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosser, A., Stover, C., Pareek, M., and Mukamolova, G. V. (2017). Resuscitation-promoting factors are important determinants of the pathophysiology in Mycobacterium tuberculosis infection. Crit. Rev. Microbiol. 43, 621–630. doi: 10.1080/1040841X.2017.1283485

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, D. G., VanderVen, B. C., Glennie, S., Mwandumba, H., and Heyderman, R. S. (2009). The macrophage marches on its phagosome: dynamic assays of phagosome function. Nat. Rev. Immunol. 9, 594–600. doi: 10.1038/nri2591

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryndak, M. B., Singh, K. K., Peng, Z., and Laal, S. (2015). Transcriptional profile of Mycobacterium tuberculosis replicating in type II alveolar epithelial cells. PLoS ONE 10:e0123745. doi: 10.1371/journal.pone.0123745

PubMed Abstract | CrossRef Full Text | Google Scholar

Saini, N. K., Sharma, M., Chandolia, A., Pasricha, R., Brahmachari, V., and Bose, M. (2008). Characterization of Mce4A protein of Mycobacterium tuberculosis: role in invasion and survival. BMC Microbiol. 8:200. doi: 10.1186/1471-2180-8-200

PubMed Abstract | CrossRef Full Text | Google Scholar

Scharn, C. R., Collins, A. C., Nair, V. R., Stamm, C. E., Marciano, D. K., Graviss, E. A., and Shiloh, M.U. (2016). Heme oxygenase-1 regulates inflammation and mycobacterial survival in human macrophages during Mycobacterium tuberculosis infection. J. Immunol. 196, 4641–4649. doi: 10.4049/jimmunol.1500434

CrossRef Full Text | Google Scholar

Scordo, J. M., Knoell, D. L., and Torrelles, J. B. (2016). Alveolar epithelial cells in Mycobacterium tuberculosis infection: active players or innocent bystanders? J. Innate Immun. 8, 3–14. doi: 10.1159/000439275

CrossRef Full Text | Google Scholar

Seto, S., Tsujimura, K., and Koide, Y. (2012). Coronin-1a inhibits autophagosome formation around Mycobacterium tuberculosis-containing phagosomes and assists mycobacterial survival in macrophages. Cell. Microbiol. 14, 710–727. doi: 10.1111/j.1462-5822.2012.01754.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shakarchi, F. I. (2015). Ocular tuberculosis: current perspectives. Clin. Ophthalmol. 9, 2223–2227. doi: 10.2147/OPTH.S65254

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, K., Gupta, A., Sharma, M., Sharma, A., Singh, R., Aggarwal, K., et al. (2018). MTBDRplus for the rapid diagnosis of ocular tuberculosis and screening of drug resistance. Eye 32, 451–456. doi: 10.1038/eye.2017.214

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, M., Bose, M., Abhimanyu Sharma, L., Diwakar, A., Kumar, S., and Banavalikar, J. N. (2012). Intracellular survival of Mycobacterium tuberculosis in macrophages is modulated by phenotype of the pathogen and immune status of the host. Int. J. Mycobacteriol. 1, 65–74. doi: 10.1016/j.ijmyco.2012.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, S., Ryndak, M. B., Aggarwal, A. N., Yadav, R., Sethi, S., Masih, S., et al. (2017). Transcriptome analysis of mycobacteria in sputum samples of pulmonary tuberculosis patients. PLoS ONE 12:e0173508. doi: 10.1371/journal.pone.0173508

PubMed Abstract | CrossRef Full Text | Google Scholar

Shechter, R., London, A., and Schwartz, M. (2013). Orchestrated leukocyte recruitment to immune-privileged sites: absolute barriers versus educational gates. Nat. Rev. Immunol. 13, 206–218. doi: 10.1038/nri3391

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimono, N., Morici, L., Casali, N., Cantrell, S., Sidders, B., Ehrt, S., et al. (2003). Hypervirulent mutant of Mycobacterium tuberculosis resulting from disruption of the mce1 operon. Proc. Natl. Acad. Sci. U.S.A. 100, 15918–15923. doi: 10.1073/pnas.2433882100

CrossRef Full Text | Google Scholar

Singh, K. H., Jha, B., Dwivedy, A., Choudhary, E., Arpitha, G. N., Ashraf, A., et al. (2017). Characterization of a secretory hydrolase from Mycobacterium tuberculosis sheds critical insight into host lipid utilization by M. tuberculosis. J. Biol. Chem. 292, 11326–11335. doi: 10.1074/jbc.M117.794297

PubMed Abstract | CrossRef Full Text | Google Scholar

Takade, A., Umeda, A., Matsuoka, M., Yoshida, S., Nakamura, M., and Amako, K. (2003). Comparative studies of the cell structures of Mycobacterium leprae and M. tuberculosis using the electron microscopy freeze-substitution technique. Microbiol. Immunol. 47, 265–270. doi: 10.1111/j.1348-0421.2003.tb03394.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Talaat, A. M., Lyons, R., Howard, S. T., and Johnston, S. A. (2004). The temporal expression profile of Mycobacterium tuberculosis infection in mice. Proc. Natl. Acad. Sci. U.S.A. 101, 4602–4607. doi: 10.1073/pnas.0306023101

CrossRef Full Text | Google Scholar

Tufariello, J. M., Chan, J., and Flynn, J. L. (2003). Latent tuberculosis: mechanisms of host and bacillus that contribute to persistent infection. Lancet. Infect. Dis. 3, 578–590. doi: 10.1016/S1473-3099(03)00741-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Van der Wel, N., Hava, D., Houben, D., Fluitsma, D., van Zon, M., Pierson, J., et al. (2007). M. tuberculosis and M. leprae translocate from the phagolysosome to the cytosol in myeloid cells. Cell 129, 1287–1298. doi: 10.1016/j.cell.2007.05.059

PubMed Abstract | CrossRef Full Text

Via, L. E., Deretic, D., Ulmer, R. J., Hibler, N. S., Huber, L. A., and Deretic, V. (1997). Arrest of mycobacterial phagosome maturation is caused by a block in vesicle fusion between stages controlled by rab5 and rab7. J. Biol. Chem. 272, 13326–13331. doi: 10.1074/jbc.272.20.13326

PubMed Abstract | CrossRef Full Text | Google Scholar

Vir, P., Gupta, D., Agarwal, R., and Verma, I. (2014). Interaction of alveolar epithelial cells with CFP21, a mycobacterial cutinase-like enzyme. Mol. Cell. Biochem. 396, 187–199. doi: 10.1007/s11010-014-2154-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, S. K., Abomoelak, B., Marcus, S. A., and Talaat, A. M. (2010). Transcriptional profiling of Mycobacterium tuberculosis during infection: Lessons learned. Front. Microbiol. 1:121. doi: 10.3389/fmicb.2010.00121

PubMed Abstract | CrossRef Full Text | Google Scholar

Warner, D. F. (2014). Mycobacterium tuberculosis metabolism. Cold Spring Harb. Perspect. Med. 5:a021121. doi: 10.1101/cshperspect.a021121

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhan, L., Tang, J., Sun, M., and Qin, C. (2017). Animal models for tuberculosis in translational and precision medicine. Front. Microbiol. 8:717. doi: 10.3389/fmicb.2017.00717

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, F., and Xie, J.-P. (2011). Mammalian cell entry gene family of Mycobacterium tuberculosis. Mol. Cell. Biochem. 352, 1–10. doi: 10.1007/s11010-011-0733-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G., Zhang, L., Zhang, M., Pan, L., Wang, F., Huang, J., et al. (2009). Screening and assessing 11 Mycobacterium tuberculosis proteins as potential serodiagnostical markers for discriminating TB patients from BCG vaccinees. Genomics Proteomics Bioinform. 7, 107–115. doi: 10.1016/S1672-0229(08)60039-X

PubMed Abstract | CrossRef Full Text | Google Scholar