Development of the first oligonucleotide microarray for global gene expression profiling in guinea pigs: defining the transcription signature of infectious diseases
© Jain et al.; licensee BioMed Central Ltd. 2012
Received: 22 February 2012
Accepted: 24 September 2012
Published: 2 October 2012
The Guinea pig (Cavia porcellus) is one of the most extensively used animal models to study infectious diseases. However, despite its tremendous contribution towards understanding the establishment, progression and control of a number of diseases in general and tuberculosis in particular, the lack of fully annotated guinea pig genome sequence as well as appropriate molecular reagents has severely hampered detailed genetic and immunological analysis in this animal model.
By employing the cross-species hybridization technique, we have developed an oligonucleotide microarray with 44,000 features assembled from different mammalian species, which to the best of our knowledge is the first attempt to employ microarray to study the global gene expression profile in guinea pigs. To validate and demonstrate the merit of this microarray, we have studied, as an example, the expression profile of guinea pig lungs during the advanced phase of M. tuberculosis infection. A significant upregulation of 1344 genes and a marked down regulation of 1856 genes in the lungs identified a disease signature of pulmonary tuberculosis infection.
We report the development of first comprehensive microarray for studying the global gene expression profile in guinea pigs and validation of its usefulness with tuberculosis as a case study. An important gap in the area of infectious diseases has been addressed and a valuable molecular tool is provided to optimally harness the potential of guinea pig model to develop better vaccines and therapies against human diseases.
KeywordsMicroarray Guinea pig Transcriptional profile Infectious diseases Tuberculosis
Pioneering research of Robert Koch in the guinea pigs (Cavia porcellus) laid the foundation of microbiology and popularized this animal model for the study of various infectious diseases. Since then, guinea pig has proven to be an indispensable tool in studying a range of infectious diseases, such as tuberculosis (TB) (Mycobacterium tuberculosis), anthrax (Bacillus anthracis), diphtheria (Corynebacterium diphtheriae), legionellosis (Legionella pneumophila), sexually transmitted diseases such as chlamydia infection (Chlamydia trachomatis), syphilis (Treponema pallidum) etc.. These studies showed that guinea pigs bear substantial similarity to humans with respect to thymic and bone marrow physiology, innate immune responses, complement system, lung physiology, corticosteroid response, requirement of exogenous source of vitamin C and delayed-type hypersensitivity reaction to infections[1, 2]. Due to these features, guinea pigs emerged as one of the best models for testing bio-defense agents and vaccines and for toxicological studies. However, despite the tremendous contribution of this animal to medical research, the paucity of immunological reagents and molecular tools and the lack of fully annotated guinea pig genome sequence has severely hampered a holistic analysis of host responses in this model. Unlike the mouse model, gene deletion technology (for example, gene knockout and knock-in) and trans-gene expression is also unavailable in case of guinea pigs. Thus, considering the biological relevance of guinea pig model in a large number of infectious diseases, the importance of developing a comprehensive guinea pig microarray cannot be over emphasized. Recently, Tree and colleagues developed an oligonucleotide microarray comprising of 86 genes to study the host response to BCG vaccination in guinea pigs, however, despite its usefulness in immune response studies, it has limited applications due to a small number of genes in the array[3, 4]. In this study, we have overcome this limitation by extending the cross-species hybridization technique to a number of mammalian species such as Human, Mouse, Rat, Macaque, Horse, Cat, Sheep, Pig, Chinchilla, Chimpanzee, Gray tailed opussum and Cattle in order to expand the number of features and have developed a 44 K guinea pig oligonucleotide microarray (GPOM). Further, to validate the array, the transcriptome of guinea pig lungs was analysed post - M. tuberculosis infection. Global gene expression profiling in this study not only facilitated the analysis of immunologically relevant genes but also helped in describing the transcriptional signature of pulmonary granulomas in guinea pigs represented by several key genes and commensurate pathways that are modulated in response to M. tuberculosis infection.
Microarray design and annotation
Probe distribution in 44 K GPOM: The table depicts the number of features derived from various mammalian species that have been used for designing the 44 K GPOM
Sequence data source
Number of Probes in 44 K Array
Human (Homo sapiens)
Agilent Catalogue Arrays
Mouse (Mus musculus)
Agilent Catalogue Arrays
Rat (Rattus norvegicus)
Agilent Catalogue Arrays
Rhesus Monkey (Macaca mulatta)
Agilent Catalogue Arrays
Dog (Canis familiaris)
Agilent Catalogue Arrays
Horse (Equus caballus)
NCBI, mRNA sequences
Cat (Felis catus)
NCBI, mRNA sequences
Sheep (Ovis aries)
NCBI, mRNA sequences
Pig (Sus scrofa)
NCBI, mRNA sequences
Guinea pig (Cavea porcellus)
NCBI, mRNA sequences
NCBI, mRNA sequences
Chimpanzee (Pan troglodyte)
NCBI, mRNA sequences
Gray tailed opussum (Monodelphis domestica)
NCBI, mRNA sequences
Cattle (Bos taurus)
NCBI, mRNA sequences
Total Number of probes
Bacillary load and granulomatous pathology in the lungs of M. tuberculosis infected guinea pigs
Gene expression profiling of pulmonary tuberculosis in guinea pigs by employing GPOM
We next compared the pulmonary gene expression profile of infected and uninfected guinea pigs by employing the GPOM developed in this study. As depicted in Additional file2, clustered heat maps were obtained for all the genes on the 44 K GPOM in case of infected guinea pigs compared to uninfected control. Since, a small perturbation in the gene expression may considerably influence the biological response, a small difference in the fold change in gene expression are also relevant. Hence, genes exhibiting ≥1.5 fold difference in gene expression with a statistical significance of p < 0.05 were considered as differentially regulated. The rationale behind the selection of the cutoff for considerable fold difference while comparing the gene expression is based on standard practice in microarray data analysis[5–10].
Regulation of key metabolic pathways in guinea pig lungs infected with M. tuberculosis
NDUFA12, NDUFV1, UQCRC2, SDHB, NDUFAB1, ATP5C1, PPA1, ATP6V0E1, ATP5J, NDUFB8, ATP5A1, NDUFA7, COX5A, LHPP, UQCRC1, ATP5G3, NDUFA11, ATP5J2, NDUFV3, NDUFV2, NDUFS6, NDUFS8, NDUFS3, ATP6V1D, COX5B, COX6A2, COX4I1, ATP6V1F, ATP6V0D1, COX7A2L, ATP6V0C, ATP6V1E1, ATP6V0B, ATP6V1B2, NDUFB4, NDUFB2, NDUFA3, COX17
ATP5O, ATP5F1, ATP12A, NDUFB5, NDUFB7, NDUFA4
LDHA, ENO1, HK1, ADH5, TPI1, GAPDH, ADH1A, AKR1A1, ALDOA, PGK2, ENO3, GPI
ALDOC, DLAT, PFKL
FH, SDHB, PC, MDH2, IDH1, ACO2, SUCLG1
Lipid, Glycerolipid and Glycerophospholipid metabolism
ABHD2, ABHD8, ACSL5, ADH1A, ADH5, AGPAT2, AKR1A1, CPT2, CEPT1, CYP27A1, CYP2B18A, CYP26B1, DGKA, FDPS, GCDH, HADHA, HADHB, LPL, PLA2G2A, PGS1
AGPAT1, AGPAT3, ACSL1, CPT1B, CYP2F1, CYP2A13, CYP1A2, CYP11A1, CYP2D40, PLD2, PPAP2A, PLA2G2C, PLA2G5, PAFAH1B1, PNPLA3, PTDSS2, LPGAT1
Purine and Pyrimidine metabolism
ADSS, NT5C2, NPR2, AK2, XDH, ADCY4, POLR2E, HPRT1, PRP, DPYD, TXNRD2, NT5C2, CAD
AK1, DCTD, POLA1, POLR2F, POLR2K, PRPS1L1, GUK1, PDE6G, PDE6A, POLR3GL, ENPP1, ZNRD1, PRUNE, TYMS, TK1
Amino acid metabolism
GLUL, GCDH, HADHA, HADHB, PLOD3, TMLHE, SETD1A, SETD8, HSD17B10, IVD, LDHA, MIF, PRDX6
ALDH4A1, ALDH6A1, ASS1, AMT, AGXT, AGMAT, ABAT, AHCYL1, CYP1A2, CKB, DNMT1, DBT, DDC, GLS, GLUD1, GLDC, GNMT, HMGCS1, NOS1, MAT2A, ODC1, WHSC1L1, SUV420H2, SUV420H1
Regulation of key signaling pathways in guinea pig lungs infected with M. tuberculosis
MRAS, MAPK7, RAC1, TGFB1, MAPKAPK5, MAPKAPK2, FGFR3, IL1B, PLA2G2A, RAC2, RAP1B, SOS2, MAPK8IP3, RRAS, HSPB1, CACNA1E, ARRB2, MAP2K2, FGF17, STK4, SRF, PAK2
MEF2C, CACNB2, CACNB4, TAOK3, CACNG4, DUSP9, FGFR4, PPP5C, FGFR2, FGFR1, RAP1A, CACNA1C, FOS, TNF, RPS6KA1, TGFB3, PPM1A, FGF19, ATF4, AKT3, CRKL, MAPK11, FGF8, MAPK1, FGF2, TP53, PPP3CB, MAPK8IP2, FGF12, FGF10, MAP3K4, PLA2G2C, PLA2G5, PRKACG, PRKACA, CHUK, PTPRR, RASA1, CACNA1S, BRAF, MAP2K6, MAP2K1, DUSP2, DUSP4, RASGRP2
RAC1, CUL1, CSNK2A2, PPP2R1A, RAC2, CSNK1A1, PPP2CA, RHOA, PSEN1
PRICKLE1, CSNK1E, LRP5, PLCB3, EP300, TBL1XR1, CTBP1, WNT3, CAMK2G, CAMK2A, CCND1, NKD2, SMAD4, SMAD2, PPP2R5D, TP53, PPP3CB, CSNK2A1, CTNNB1, PRKACG, PRKACA, CAMK2D, LEF1, WNT10A, FBXW11, BTRC, FZD9, FZD3, TCF7L2
VDAC1, ERBB2, ATP2A1, GNAQ, TNNC1, CACNA1E, P2RX2, ADCY4, SLC25A6
PLCB3, ADRA1B, PHKG2, CACNA1C, ADORA2A, CAMK2G, CAMK2A, GNAS, PLN, HTR6, NOS1, RYR1, VDAC2, PPP3CB, SLC8A1, CYSLTR1, P2RX1, CALM2, ATP2B4, ATP2A3, PRKACG, PRKACA, CAMK2D, LTB4R2, CACNA1S, ADRB2, PLCG1
Regulation of key immune response related genes in guinea pig lungs infected with M. tuberculosis
Chemokine signaling pathway
GNB4, RAC1, IL8, GNG5, HCK, RAC2, FOXO3, RAP1B, SOS2, ADRBK1, NCF1, GNAI2, CCL27, VAV3, RHOA, ADCY4, ARRB2, NFKBIB
GNB1, PLCB3, CCL5, CXCL9, RAP1A, CXCR3, GNG8, GNG2, CCL21, AKT3, CRKL, MAPK1, PRKACG, PRKACA, CHUK, CCL11, BRAF, MAP2K1, RASGRP2
Cell adhesion molecules (CAMs)
HLA-DQB1, ITGB7, HLA-A, HLA-B, HLA-B, HLA-DMA, NRXN2, HLA-DPA1, HLA-DRB1, HLA-DRB1, ALCAM, CLDN18, CLDN14
MPZL1, HLA-E, HLA-F, CADM3, CADM1, CD34, CD8A, CD99, CDH3, CLDN4, CLDN6, NCAM1, GLG1, NLGN2, ITGB2, ITGB8, ITGA9, CNTNAP1, SIGLEC1
Cytokine-cytokine receptor interaction
IL1, RAP, IL8, TGFB1, IL1B, CSF1R, IL23A, BMP8A, CCL27, BMP4, BMPR2
GHR, CCL5, CD70, CXCL9, KIT, CD27, CXCR3, TNFα, IL2RG, CNTFR, BMPR1A, IL23R, CCL21, TGFB3, FLT4, BMP5, IFNGR2, LTA, IL9, CCL11, TNFSF10, ACVR2A
Leukocyte trans-endothelial migration
RAC1, CYBA, CYBB, RAC2, RAP1B, ACTG1, NCF1, GNAI2, EZR, ACTN4, VAV3, RHOA, PTPN11, CLDN18, CLDN14, ACTB
RAP1A, CD99, RHOH, THY1, MAPK11, CTNNB1, CLDN4, CLDN6, ITGB2, MLLT4, PLCG1
Natural killer cell mediated cytotoxicity
RAC1, ARAF, HLA-A, HLA-B, RAC2, SOS2, VAV3, LAT, FCER1G, PTPN11, MAP2K2, PTPN6
TYROBP, HLA-E, HLA-F, TNF, IFNGR2, MAPK1, PPP3CB, GZMB, ITGB2, TNFSF10, BRAF, MAP2K1, PLCG1
T cell receptor signaling pathway
SOS2, GRAP2, VAV3, LAT, RHOA, NFKBIB, MAP2K2, PTPN6, PAK2
CD8A, FOS, TNF, NCK1, AKT3, MAPK11, MAPK1, PPP3CB, CHUK, MAP2K1, PLCG1
B cell receptor signaling pathway
RAC1, RAC2, SOS2, CD81, VAV3, LILRB3, NFKBIB, MAP2K2, PTPN6, CR2
FOS, AKT3, MAPK1, PPP3CB, CHUK, MAP2K1
Antigen processing and presentation
HLA-DQB1, HLA-A, HLA-B, HLA-B, CD74, HLA-DMA, HLA-DPA1, HLA-DRB1, HLA-DRB1, HSP90AB1, CTSB, CTSS, HSP90AA1
HLA-E, HLA-F, CD8A, TNF, TAPBP, HSPA5, PSME2
Complement and coagulation cascades
C4A, C3, C1S, SERPING1, CR2
CD55, MASP2, KNG1, SERPINC1, F8, C9, C6, C3AR1, C5AR1
Pathway analysis revealed that more than 196 genes that were differentially regulated in infected samples as compared to the uninfected control samples were associated with general metabolism (Additional file4). In particular, several key genes (Table2) involved in oxidative phosphorylation (44 genes), purine and pyrimidine metabolism (28 genes, respectively), amino acid metabolism (37 genes), glycolysis/gluconeogenesis (15 genes), glycerophospholipid and glycerolipid metabolism (13 and 8 genes, respectively), fatty acid (9 genes) and arachidonic acid metabolism (9 genes), pentose phosphate pathway (9 genes), TCA cycle (8 genes) and ether metabolism (7 genes) were significantly regulated in response to infection (Additional file4).
Among the cell signaling pathways, the most profound modulation was observed in the MAPK signaling pathway with perturbation of 67 genes (Table3). In addition, the expression of several key genes related to Wnt (38 genes), calcium (36 genes), neurotrophin (34 genes), insulin (28 genes), TGFβ (25 genes), ErbB (18 genes), JAK-STAT (17 genes), PPAR (15 genes), VEGF (15 genes), hedgehog (14 genes), p53 (11 genes), notch (8 genes) and phosphatidylinositol signaling systems (4 genes) were also affected (Additional file4). Characteristically, majority of the genes related to MAPK, Wnt and Calcium signaling pathways exhibited a marked down regulation (Table3).
Among the immune response related genes, a significant alteration was observed in the expression of several key genes associated with chemokine signaling pathway (37 genes), cell adhesion molecules (CAMS) (32 genes), cytokine-cytokine receptor interaction (33 genes), leukocyte trans-endothelial migration (27 genes), natural killer cell mediated cytotoxicity (25 genes), T cell receptor signaling (20 genes) and B cell receptor signaling (16 genes) (Table4). In addition, profound modulation was observed in genes associated with the antigen processing and presentation (20 genes), hematopoietic cell lineage (16 genes) and complement and coagulation cascades (14 genes) (Table4 and Additional file4). A marked regulation was also observed in the signaling pathways related to TGF-β (25 genes), Toll-like receptor (19 genes), Fc epsilon RI (18 genes), adipocytokine (13 genes), NOD-like receptor (13 genes), RIG-I-like receptor (12 genes), mTOR (11 genes) and intestinal immune network for IgA production (8 genes) (Additional file4).
Validation of microarray results by real time RT-PCR
For the validation of 44 K GPOM, 5 genes were selected based on the following criteria: (i) differential regulation, (ii) immunological relevance and (iii) availability of cDNA sequence in the NCBI database. Expression of these genes was analyzed on the same RNA samples, which were used for the microarray study by semi- quantitative real time RT-PCR by employing SYBR green PCR Master Mix (Applied Biosystems) with 18S as the internal control. The primer sequences of C3AR1, CAMP, CCL5, IFNγ, C4BPA and 18S rRNA genes were designed using the guinea pig gene specific cDNA sequences available in the public database (NCBI) based on the recommended guidelines for designing real time PCR primers (Primer express software, Applied Biosystems]. Sequences of the primers are described in Supporting Table S2 in Additional file1. The primers were experimentally validated for two standard quality control criteria, (i) single amplicon specificity by analysis of dissociation curves and (ii) consistently high amplification efficiency by analysis of calibration curve. For analysis of real time PCR data ΔΔCt method was employed. Consistent with our microarray results, wherein several probes corresponding to C3AR1, CAMP, CCL5 and IFNγ exhibited a down regulation along with up regulation of C4BPA, real time PCR results also matched the microarray expression profile (Supporting Table S2 in Additional file1).
Guinea pig model has made tremendous contribution towards the understanding of several infectious and non-infectious human diseases[1, 2]. Moreover, in case of TB, guinea pig has emerged as the most preferred and biologically relevant model to investigate TB pathogenesis and therapy. The success of M. tuberculosis as an intracellular pathogen is primarily attributable to its ability to reside in human lungs inside hypoxic granulomas in a dormant stage for years or even decades. However, the conventional mouse model for TB does not form hypoxic granulomas, which serve as the primary host defense mechanism for the containment of infection and is the central feature of TB pathogenesis in humans. Thus, due to its ability to produce hypoxic granulomas guinea pig as a surrogate model fulfills an important niche in the field of TB. However, the lack of a fully annotated genome sequence remains a major impediment towards the realization of full potential of guinea pig as an animal model. Moreover, till date not even a single study has investigated the genome-wide transcriptional response associated with M. tuberculosis infection in guinea pigs. Thus, development of a microarray platform to understand the host responses in guinea pigs is a critical step towards harnessing the true potential of this model.
In the present study, by employing cross-species hybridization technique, we have developed an oligonucleotide microarray with 44,000 features assembled from different mammalian species, which to the best of our knowledge is the first attempt to employ microarray to study global gene expression profile in guinea pigs. In order to demonstrate the utility of this microarray and to gain insight into the host responses, we have carried out the expression profile of guinea pig lungs at 10 weeks post-infection with M. tuberculosis. At this stage of infection, guinea pigs exhibit advanced stage of TB disease with multiple coalescing granulomas along with caseation and liquefaction necrosis in the lungs. Hence, the gene expression profile observed in this study represents a transcriptional signature of advanced progressive TB disease in guinea pigs.
Induction of catabolic processes with consequential ATP accumulation has recently been shown to provide an interface between metabolism and host defense to infection[18–20]. ATP molecules generated in response to injury to airway epithelial cells have been reported to be a critical determinant of cell migration and repair following the injury and have been shown to be associated with the activation of down-stream signaling cascades and induction of IL-1β through the interaction of ATP with purinergic receptors. A few in vitro studies have also indicated the role of ATP mediated macrophage apoptosis in killing of M. tuberculosis. A concurrent up regulation in the expression of oxidative phosphorylation related genes (expected to result in increased ATP levels), purinergic receptors and IL-1β in this study to the best of our knowledge, provides the first in vivo evidence for the involvement of these pathways in TB. Further, the lungs of the infected guinea pigs also exhibited a marked perturbation in the expression of several key genes associated with chemokine signaling (CCL27, CCL5, CXCL9, CXCR3, CCL21 and CCL11), cell adhesion molecules (CAMs) (HLA, ALCAM, MPZL1, CADM3, CADM1, CD34, CD8A, CD99, CDH3, CLDN4, CLDN6, NCAM1, ITGB2, ITGB8 and ITGA9) and cytokine and cytokine receptors (IL1β, IL1RAP, IL2RG, IL8, IL9, IL23A, IL23R, TGFB1, TGFB3, IFNGR2, TNFα, TNFSF10, CSF1R, BMP4, BMP8A, BMPR1A, BMPR2, LTA and ACVR2A), which are known to contribute to leukocyte trans-endothelial migration, inflammation and granulomatous pathology.
Perturbation in the cellular signaling pathways is another typical theme that emerged from our study. The most prominent observation relates to the repression of numerous genes related to MAPK, Wnt and calcium signaling pathways. These observations are consistent with previous studies, which have suggested that modulation of MAPK signaling pathway along with the reduction in the levels of intracellular calcium are some of the important means by which, M. tuberculosis represses phagosome - lysosome maturation and pro-inflammatory responses at the site of infection[22, 23]. MAPK signaling is known to be crucial for the anti-bacterial response of the host and it also represents a strategic target for bacterial subversion tactics. Thus, dampening of the MAPK signaling has emerged as a key to achieve alteration in the antibacterial phenotype of macrophages. Recently, Wnt signaling pathway has been implicated in the generation of long-lived multi-potent memory T cells and in the modulation of inflammatory response of macrophages to M. tuberculosis infection[25, 26], thus repression of Wnt signaling pathway observed in this study suggests a possible mechanism by which, M. tuberculosis inhibits effective T cell memory response.
Another key observation from this study relates to the modulation of several key genes involved in the re-modeling of extracellular matrix such as, COL6A2, COL14A1, COL12A1, MMP24, TIMP1, SERPING1, SERPINB1, ADAMTS1, ADAMTS7, MMP1, PITRM1, SERPINA3N, SERPINB6, SERPINE2, SERPINH1 and CNDP2. These observations are in agreement with the previous studies, which have implicated these genes in dual events associated with tissue remodeling as well as tissue damaging[27–30]. The presence of extensive necrosis along with thick bands of collagen observed in the guinea pig pulmonary granulomas indicates that the balance is heavily tipped in the favor of tissue damaging events. Increased expression of the complement receptor CR2 and numerous key genes involved in phagocytosis and antigen presentation as observed in this study further substantiates that the pathogen exploits the normally effective defense system to its advantage by subverting or co-opting these pathways as has been also reported earlier[17, 22].
This study reports for the first time the development of a 44 K oligonucleotide microarray for guinea pigs and provides an important tool to capture the genome wide transcriptional changes in this model. The transcriptional profiling of M. tuberculosis infected guinea pig lungs not only revealed modulation of key immunologically relevant genes but also demonstrated involvement of novel metabolic and signaling pathways in TB pathogenesis. Moreover, in silico analysis revealed a higher resemblance of guinea pigs to humans in terms of transcriptional response to M. tuberculosis infection when compared to mouse and non-human primates. Development of the 44 K GPOM is thus, a critical step towards characterization of the guinea pig model, which will greatly aid in improving our understanding of host responses to a number of infectious diseases. We believe that optimal use of guinea pig model and further research on its biology would generate tremendous opportunity to understand host-pathogen interaction and thus, help in the development of new therapeutic intervention strategies.
Guinea pig experiments were reviewed and approved by the Institutional Animal Ethics Committee of University of Delhi South Campus, New Delhi, India (Permit number: 159/1999/CPCSEA). All procedures with the infected animals/tissues were performed in a Biosafety level three (BSL-3) containment facility at the University of Delhi South Campus, according to the approved protocols. All animals were routinely cared for according to the guidelines of CPCSEA (Committee for the Purpose of Control and Supervision on Experiments on Animals), India with a 12 hr light/dark cycle (0600–1800) maintained at 25 degrees Celsius with a relative humidity of 50%.
Experimental animals, infection and study design
Pathogen free 6–8 weeks old (200-300 g) female outbred guinea pigs (Dunkin Hartley strain) used for the microarray studies were procured from Disease Free Small Animal House Facility, CCS Haryana Agricultural University, Hisar, India. The infected and uninfected animals were housed separately in individually ventilated cages (2 animals/cage) and were provided with ad libitum food and water in the BSLIII facility at University of Delhi South Campus, India. Guinea pigs were infected by using the method as described previously. Briefly, M. tuberculosis (H37Rv strain, ATCC no. 25618 procured from AIIMS, New Delhi, India) was grown to mid-log phase in Middle Brook 7 H9 media (0.05% Tween80 and 0.5% Glycerol) and stocks were prepared as described. The CFU of stocks was enumerated by plating 10 fold serial dilutions on 7 H11 agar (1XADC and 0.5% glycerol). By using pre-calibrated infection parameters, guinea pigs were infected in the inhalation exposure system, (Glas-col Inc.), which resulted in ~500 bacilli in the lungs of guinea pigs at day 1 post-infection. For enumeration of day 1 CFU, the whole lung homogenates were plated onto 7 H11 agar plates and colonies were counted after 3 weeks of incubation at 37°C.
Necropsy procedure and histopathological evaluation
Guinea pigs were euthanized by carbon dioxide asphyxiation. After aseptically dissecting the animals, for histopathological evaluation, three lung lobes (right caudal, middle and cranial) were removed and fixed in 10% neutral buffered formalin. Left caudal lung lobe was aseptically removed for the measurement of bacillary load. A portion of left cranial lung lobe was stored in RNA later® (Ambion) at −20°C for isolation of RNA to be used for microarray and real time RT-PCR studies.
For histopathological examination, as described previously sections of 5 μm thickness from formalin fixed and paraffin embedded tissues were cut on to glass slides and stained with haematoxylin and eosin. The percent granuloma in lung, type and extent of necrosis, organization of granuloma along with the type of infiltrating cells were assessed. In order to determine the extent of collagen deposition and fibrosis, the lung sections were also stained with Van Gieson stain.
Specific portions of lungs were weighed and homogenized separately in 5 ml saline in a Polytron homogenizer. Appropriate dilutions of the homogenates were inoculated on to MB7H11 agar plates in duplicates and incubated at 37°C in a CO2 incubator for three to four weeks. The number of colonies were counted and expressed as log10 CFU/g of tissue. The detection limit in case of both lungs CFU was 1.0 log10 CFU/g.
Labeling of RNA samples and quality control for 44 K microarray hybridization
To evaluate the gene expression profile associated with pulmonary tuberculosis, total RNA was isolated from lung tissues of M. tuberculosis infected and naive control guinea pigs (n = 3) by using Qiagen’s RNeasy mini kit as per the manufacturer’s recommendations followed by assessment of quality (specific activity) and quantity (yield) by using Bioanalyzer (Agilent Technologies). RNA samples then were labelled with Cy3 by using Agilent Quick-Amp labeling Kit as per the manufacturer’s recommendations. Briefly, 500 ng each of the control and test RNA samples were incubated with reverse transcription mix at 40°C and converted to double stranded cDNA primed by oligo dT with a T7 polymerase promoter. Synthesized double stranded cDNA was then used as template for cRNA generation. cRNA was generated by in vitro transcription along with the incorporation of Cy3 CTP during this step (Agilent Technologies). By using Qiagen’s RNeasy mini kit, fluorescently labeled cRNA was purified followed by the assessment of quality and quantity by using Bioanalyzer (Agilent Technologies).
Hybridization and scanning
Linear amplified Cy3 labeled cRNA were hybridized to guinea pig 44 K microarray. Briefly, 1.65 μg of Cy3 labeled cRNAs were fragmented and hybridized to the array. Fragmentation of labeled cRNA and hybridization were carried out by using the Agilent Gene Expression Hybridization kit. Hybridization was carried out in Agilent’s Surehyb Chambers at 65°C for 16 hours. Following hybridization, the slides were washed by using Agilent Gene Expression wash buffers and scanned at 3μm resolution by using the Agilent Microarray Scanner G2505C under conditions to limit saturation to < 80% and were saved as TIFF images. The features then were extracted with the Feature Extraction Software (Agilent technologies, v10.7). Details of design, development and annotation of microarray are provided in Supporting materials and methods in Additional file1.
Microarray data analysis
The data was analyzed by using GeneSpring GX v11.5 software (Agilent Technologies). Normalization of the data was carried out in GeneSpring GX by using the 75th percentile shift. It subtracts this value from the expression value of each entity and normalizes to specific control samples and fold change was calculated for the infected group relative to the baseline control (cRNA derived from uninfected animals). The fold intensity data then were filtered for significantly regulated (up and down regulated) genes in the treatment group in comparison to the control group based on the following stringent criteria: Expression fold values are provided in terms of log base 2. For up regulated genes, the cutoff used is fold change > 0.6 along with geometric mean fold change > 1 and Flags "detected" in infected samples. For filtering the down regulated genes, the cutoff used is fold change < − 0.6 along with geometric mean fold change < −1 and Flags "detected" in normal uninfected control sample that means a change in expression by ≥ 1.5 fold (up or down). Based on these criteria, only those genes, which showed a significant regulation in triplicates, were further subjected to hierarchical clustering based on Pearson coefficient correlation algorithm to identify significant gene expression patterns. Genes were classified on the basis of functional categories and pathways by using GeneSpringGX software and Genotypic Biointerpreter Biological Analysis software (Genotypic Technology Pvt. Ltd.). The microarray data reported here have been submitted at NCBI’s Gene Expression Omnibus [GEO Accession number: GSE32447].
Real time RT-PCR
For the validation of 44 K GPOM, expression of a few selected genes was analysed by real time RT-PCR by employing SYBR green PCR Master Mix (Applied Biosystems) on the same RNA samples, which were used for the microarray study. Sequences of the primers employed for real time RT-PCR of C3AR1, CAMP, CCL5, IFNγ, C4BPA and 18S rRNA genes and comparison of gene expression pattern with respect to microarray are described in Supporting Table S2 in Additional file1. For the analysis of real time PCR data ΔΔCt method was employed. First, ΔCt value was calculated for each sample as the difference between the Ct values for the gene of interest and the housekeeping gene (18S) in each sample. Then, ΔΔCt value was calculated as the difference between the ΔCt values of an experimental sample and the control sample. Fold change in gene expression were calculated as 2 -ΔΔCt.
Mean differences for percent fold induction in mRNA expression levels were analyzed by student’s t test. Differences were considered statistically significant, when p < 0.05. Based on Pearson coefficient correlation algorithm, the samples were clustered by using hierarchical clustering to identify similar conditions.
RJ, BD and AKT conceived, designed and performed the experiments, analysed the data and wrote the manuscript. All authors read and approved the final manuscript.
Guinea pig oligonucleotide microarray
Reverse Transcription – Polymerase Chain Reaction.
We acknowledge Genotypic Technology Pvt. Ltd Bangalore for the microarray processing and data analysis reported in this publication. Priyanka Chauhan is acknowledged for critical reading of the manuscript. This work was supported by a research grant from the Department of Biotechnology, Ministry of Science and Technology, Government of India, New Delhi. However, the funding agency had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.
- Padilla-Carlin DJ, McMurray DN, Hickey AJ: The guinea pig as a model of infectious diseases. Comp Med. 2008, 58 (4): 324-340.PubMed CentralPubMedGoogle Scholar
- Bloom BE: Guinea pig model of tuberculosis. 1994, Washington DC: ASM PressGoogle Scholar
- Tree JA, Patel J, Thom RE, Elmore MJ, Schafer H, Williams A, Marsh PD: Temporal changes in the gene signatures of BCG-vaccinated guinea pigs in response to different mycobacterial antigens. Vaccine. 2010, 28 (50): 7979-7986. 10.1016/j.vaccine.2010.09.061.View ArticlePubMedGoogle Scholar
- Tree JA, Elmore MJ, Javed S, Williams A, Marsh PD: Development of a guinea pig immune response-related microarray and its use to define the host response following Mycobacterium bovis BCG vaccination. Infect Immun. 2006, 74 (2): 1436-1441. 10.1128/IAI.74.2.1436-1441.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Donigan KA, Tuck D, Schulz V, Sweasy JB: DNA Polymerase beta Variant Ile260Met Generates Global Gene Expression Changes Related to Cellular Transformation. Mutagenesis. 2012, 10.1093/mutage/ges034Google Scholar
- Kaur S, Archer KJ, Devi MG, Kriplani A, Strauss JF, Singh R: Differential Gene Expression in Granulosa Cells from Polycystic Ovary Syndrome Patients with and without Insulin Resistance: Identification of Susceptibility Gene Sets through Network Analysis. J Clin Endocrinol Metab. 2012, 10.1210/jc.2011-3441Google Scholar
- Guo NL, Wan YW, Denvir J, Porter DW, Pacurari M, Wolfarth MG, Castranova V, Qian Y: Multiwalled carbon nanotube-induced gene signatures in the mouse lung: potential predictive value for human lung cancer risk and prognosis. J Toxicol Environ Health A. 2012, 75 (18): 1129-1153. 10.1080/15287394.2012.699852.PubMed CentralView ArticlePubMedGoogle Scholar
- Milan M, Coppe A, Reinhardt R, Cancela LM, Leite RB, Saavedra C, Ciofi C, Chelazzi G, Patarnello T, Bortoluzzi S, et al: Transcriptome sequencing and microarray development for the Manila clam. Ruditapes philippinarum: genomic tools for environmental monitoring. BMC Genomics. 2011, 12: 234-PubMedGoogle Scholar
- Fisher MA, Plikaytis BB, Shinnick TM: Microarray analysis of the Mycobacterium tuberculosis transcriptional response to the acidic conditions found in phagosomes. J Bacteriol. 2002, 184 (14): 4025-4032. 10.1128/JB.184.14.4025-4032.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Lorenz WW, Alba R, Yu YS, Bordeaux JM, Simoes M, Dean JF: icroarray analysis and scale-free gene networks identify candidate regulators in drought-stressed roots of loblolly pine (P. taeda L.). BMC Genomics. 2011, 12: 264-10.1186/1471-2164-12-264.PubMed CentralView ArticlePubMedGoogle Scholar
- Via LE, Lin PL, Ray SM, Carrillo J, Allen SS, Eum SY, Taylor K, Klein E, Manjunatha U, Gonzales J, et al: Tuberculous granulomas are hypoxic in guinea pigs, rabbits, and nonhuman primates. Infect Immun. 2008, 76 (6): 2333-2340. 10.1128/IAI.01515-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Khovidhunkit W, Kim MS, Memon RA, Shigenaga JK, Moser AH, Feingold KR, Grunfeld C: Effects of infection and inflammation on lipid and lipoprotein metabolism: mechanisms and consequences to the host. J Lipid Res. 2004, 45 (7): 1169-1196. 10.1194/jlr.R300019-JLR200.View ArticlePubMedGoogle Scholar
- Bordbar A, Lewis NE, Schellenberger J, Palsson BO, Jamshidi N: Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions. Mol Syst Biol. 2010, 6: 422-PubMed CentralView ArticlePubMedGoogle Scholar
- Oscar B, Bay ML, Besedovsky H, Rey DA: The Immune-Endocrine-Metabolic Unit During Human Tuberculosis. Current Immunology Reviews. 2010, 6 (9): 314-322.Google Scholar
- Kim MJ, Wainwright HC, Locketz M, Bekker LG, Walther GB, Dittrich C, Visser A, Wang W, Hsu FF, Wiehart U, et al: Caseation of human tuberculosis granulomas correlates with elevated host lipid metabolism. EMBO Mol Med. 2010, 2 (7): 258-274. 10.1002/emmm.201000079.PubMed CentralView ArticlePubMedGoogle Scholar
- Day TA, Koch M, Nouailles G, Jacobsen M, Kosmiadi GA, Miekley D, Kuhlmann S, Jorg S, Gamradt P, Mollenkopf HJ, et al: Secondary lymphoid organs are dispensable for the development of T-cell-mediated immunity during tuberculosis. Eur J Immunol. 2010, 40 (6): 1663-1673. 10.1002/eji.201040299.View ArticlePubMedGoogle Scholar
- Mehra S, Pahar B, Dutta NK, Conerly CN, Philippi-Falkenstein K, Alvarez X, Kaushal D: Transcriptional reprogramming in nonhuman primate (rhesus macaque) tuberculosis granulomas. PLoS One. 2010, 5 (8): e12266-10.1371/journal.pone.0012266.PubMed CentralView ArticlePubMedGoogle Scholar
- Wesley UV, Bove PF, Hristova M, McCarthy S, van der Vliet A: Airway epithelial cell migration and wound repair by ATP-mediated activation of dual oxidase 1. J Biol Chem. 2007, 282 (5): 3213-3220.View ArticlePubMedGoogle Scholar
- Placido R, Auricchio G, Falzoni S, Battistini L, Colizzi V, Brunetti E, Di Virgilio F, Mancino G: P2X(7) purinergic receptors and extracellular ATP mediate apoptosis of human monocytes/macrophages infected with Mycobacterium tuberculosis reducing the intracellular bacterial viability. Cell Immunol. 2006, 244 (1): 10-18. 10.1016/j.cellimm.2007.02.001.View ArticlePubMedGoogle Scholar
- Kusner DJ, Barton JA: ATP stimulates human macrophages to kill intracellular virulent Mycobacterium tuberculosis via calcium-dependent phagosome-lysosome fusion. J Immunol. 2001, 167 (6): 3308-3315.View ArticlePubMedGoogle Scholar
- Netea MG, Simon A, van de Veerdonk F, Kullberg BJ, Van der Meer JW, Joosten LA: IL-1beta processing in host defense: beyond the inflammasomes. PLoS Pathog. 2010, 6 (2): e1000661-10.1371/journal.ppat.1000661.PubMed CentralView ArticlePubMedGoogle Scholar
- Malik ZA, Denning GM, Kusner DJ: Inhibition of Ca(2+) signaling by Mycobacterium tuberculosis is associated with reduced phagosome-lysosome fusion and increased survival within human macrophages. J Exp Med. 2000, 191 (2): 287-302. 10.1084/jem.191.2.287.PubMed CentralView ArticlePubMedGoogle Scholar
- Salam N, Gupta S, Sharma S, Pahujani S, Sinha A, Saxena RK, Natarajan K: Protective immunity to Mycobacterium tuberculosis infection by chemokine and cytokine conditioned CFP-10 differentiated dendritic cells. PLoS One. 2008, 3 (8): e2869-10.1371/journal.pone.0002869.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenberger CM, Finlay BB: Phagocyte sabotage: disruption of macrophage signalling by bacterial pathogens. Nat Rev Mol Cell Biol. 2003, 4 (5): 385-396. 10.1038/nrm1104.View ArticlePubMedGoogle Scholar
- Gattinoni L, Zhong XS, Palmer DC, Ji Y, Hinrichs CS, Yu Z, Wrzesinski C, Boni A, Cassard L, Garvin LM, et al: Wnt signaling arrests effector T cell differentiation and generates CD8+ memory stem cells. Nat Med. 2009, 15 (7): 808-813. 10.1038/nm.1982.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaale K, Neumann J, Schneider D, Ehlers S, Reiling N: Wnt signaling in macrophages: augmenting and inhibiting mycobacteria-induced inflammatory responses. Eur J Cell Biol. 2011, 90 (6-7): 553-559. 10.1016/j.ejcb.2010.11.004.View ArticlePubMedGoogle Scholar
- Taylor JL, Hattle JM, Dreitz SA, Troudt JM, Izzo LS, Basaraba RJ, Orme IM, Matrisian LM, Izzo AA: Role for matrix metalloproteinase 9 in granuloma formation during pulmonary Mycobacterium tuberculosis infection. Infect Immun. 2006, 74 (11): 6135-6144. 10.1128/IAI.02048-05.PubMed CentralView ArticlePubMedGoogle Scholar
- Elkington PT, Nuttall RK, Boyle JJ, O'Kane CM, Horncastle DE, Edwards DR, Friedland JS: Mycobacterium tuberculosis, but not vaccine BCG, specifically upregulates matrix metalloproteinase-1. Am J Respir Crit Care Med. 2005, 172 (12): 1596-1604. 10.1164/rccm.200505-753OC.View ArticlePubMedGoogle Scholar
- Elkington PT, O'Kane CM, Friedland JS: The paradox of matrix metalloproteinases in infectious disease. Clin Exp Immunol. 2005, 142 (1): 12-20. 10.1111/j.1365-2249.2005.02840.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Izzo AA, Izzo LS, Kasimos J, Majka S: A matrix metalloproteinase inhibitor promotes granuloma formation during the early phase of Mycobacterium tuberculosis pulmonary infection. Tuberculosis (Edinb). 2004, 84 (6): 387-396. 10.1016/j.tube.2004.07.001.View ArticleGoogle Scholar
- Jain R, Dey B, Dhar N, Rao V, Singh R, Gupta UD, Katoch VM, Ramanathan VD, Tyagi AK: Enhanced and enduring protection against tuberculosis by recombinant BCG-Ag85C and its association with modulation of cytokine profile in lung. PLoS One. 2008, 3 (12): e3869-10.1371/journal.pone.0003869.PubMed CentralView ArticlePubMedGoogle Scholar
- Singh R, Rao V, Shakila H, Gupta R, Khera A, Dhar N, Singh A, Koul A, Singh Y, Naseema M, et al: Disruption of mptpB impairs the ability of Mycobacterium tuberculosis to survive in guinea pigs. Mol Microbiol. 2003, 50 (3): 751-762. 10.1046/j.1365-2958.2003.03712.x.View ArticlePubMedGoogle Scholar