Transcriptional response of peripheral blood mononuclear cells from cattle infected with Mycobacterium bovis

Integrative genomics sheds light on the immunogenetics of tuberculosis in cattle

Mycobacterium bovis causes bovine tuberculosis (bTB), an infectious disease of cattle that represents a zoonotic threat to humans. Research has shown that the peripheral blood (PB) transcriptome is perturbed during bTB disease but the genomic architecture underpinning this transcriptional response remains poorly understood. Here, we analyse PB transcriptomics data from 63 control and 60 confirmed M. bovis-infected animals and detect 2592 differently expressed genes perturbing multiple immune response pathways. Leveraging imputed genome-wide SNP data, we characterise thousands of cis-expression quantitative trait loci (eQTLs) and show that the PB transcriptome is substantially impacted by intrapopulation genomic variation during M. bovis infection. Integrating our cis-eQTL data with bTB susceptibility GWAS summary statistics, we perform a transcriptome-wide association study and identify 115 functionally relevant genes (including RGS10, GBP4, TREML2, and RELT) and provide important new omics data for understanding the host response to mycobacterial infections that cause tuberculosis in mammals.

Assessment of tuberculosis biomarkers in paratuberculosis-infected cattle

Introduction

Mycobacterium bovis and Mycobacterium avium subsp. paratuberculosis, respectively the causative agents of bovine tuberculosis (bTB) and bovine paratuberculosis (PTB), share a high number of antigenic proteins. This characteristics makes the differential diagnosis of the diseases difficult. The interferon gamma (IFN-γ), C-X-C motif chemokine ligand 10 (CXCL10), matrix metallopeptidase 9 (MMP9), interleukin 22 (IL-22) and thrombospondin 1 (THBS1) bovine genes have already been shown to be accurate transcriptional biomarkers of bTB. In order to improve the diagnosis of bTB and PTB, in the present study we evaluated the risk of false positivity of these bTB biomarkers in cattle with PTB.

Material and Methods

The transcription of these genes was studied in 13 PTB-infected cattle, using Mycobacterium avium subsp. paratuberculosis (MAP)-stimulated peripheral blood mononuclear cells (PBMC).

Results

Overall, the levels of IFN-γ, CXCL10, MMP9 and IL-22 transcripts in MAP-stimulated PBMC failed to differentiate animals with PTB from healthy animals. However, as bTB-afflicted cattle do, the MAP-infected group also displayed a lower level of THBS1 transcription than the non-infected animals.

Conclusion

The results of this study add new specificity attributes to the levels of transcription of IFN-γ, CXCL10, MMP9 and IL-22 as biomarkers for bTB.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

Objective

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

Methods

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

Results

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

Conclusion

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

Revisiting tuberculosis screening: An insight to complementary diagnosis and prospective molecular approaches for the recognition of the dormant TB infection in human and cattle hosts

Tuberculosis (TB) is defined as a chronic infection in both human and cattle hosts and many subclinical cases remain undetected. After the pathogen is inhaled by a host, phagocyted bacilli can persist inside macrophages surviving intracellularly. Hosts develop granulomatous lesions in the lungs or lymph nodes, limiting infection. However, bacilli become persister cells. Immunological diagnosis of TB is performed basically by routine tuberculin skin test (TST), and in some cases, by ancillary interferon-gamma release assay (IGRA). The concept of human latent TB infection (LTBI) by M. tuberculosis is recognized in cohorts without symptoms by routine clinical diagnostic tests, and nowadays IGRA tests are used to confirm LTBI with either active or latent specific antigens of M. tuberculosis. On the other hand, dormant infection in cattle by M. bovis has not been described by TST or IGRA testing as complications occur by cross-reactive immune responses to homolog antigens of environmental mycobacteria or a false-negative test by anergic states of a wained bovine immunity, evidencing the need for deciphering more specific biomarkers by new-generation platforms of analysis for detection of M. bovis dormant infection. The study and description of bovine latent TB infection (boLTBI) would permit the recognition of hidden animal infection with an increase in the sensitivity of routine tests for an accurate estimation of infected dairy cattle. Evidence of immunological and experimental analysis of LTBI should be taken into account to improve the study and the description of the still neglected boLTBI.

RNA-Seq Transcriptome Analysis of Peripheral Blood From Cattle Infected With Mycobacterium bovis Across an Experimental Time Course

Bovine tuberculosis, caused by infection with members of the Mycobacterium tuberculosis complex, particularly Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including RNA sequencing, has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analysed the transcriptome of bovine whole peripheral blood samples collected at −1 week pre-infection and +1, +2, +6, +10, and +12 weeks post-infection time points. Differentially expressed genes were catalogued and evaluated at each post-infection time point relative to the −1 week pre-infection time point and used for the identification of putative candidate host transcriptional biomarkers for M. bovis infection. Differentially expressed gene sets were also used for examination of cellular pathways associated with the host response to M. bovis infection, construction of de novo gene interaction networks enriched for host differentially expressed genes, and time-series analyses to identify functionally important groups of genes displaying similar patterns of expression across the infection time course. A notable outcome of these analyses was identification of a 19-gene transcriptional biosignature of infection consisting of genes increased in expression across the time course from +1 week to +12 weeks post-infection.

Assessment of candidate biomarkers to detect resistance to Mycobacterium bovis in Holstein-Friesian cattle

Bovine tuberculosis (bTB) caused by Mycobacterium bovis has a significant economic impact worldwide each year. Control of bTB is based on skin testing and removal of reactors. However, additional strategies are required to control this disorder. Natural disease resistance has been defined as the inherent capacity of an individual to resist disease when exposed to pathogens without previous exposure or immunization. However, little is known about natural disease resistance against Mycobacterium bovis in cattle. In this study, we aimed to identify candidate biomarkers to detect host resistance to M. bovis. We used a microbicidal assay to identify the resistance phenotype. A genomic microarray analysis was carried out on RNA from 2 resistant (R) and 2 susceptible (S) cows. Our results evidenced 69 differentially expressed genes. A subset of six genes that showed differential up (IL1RN), and down-regulation (VNN, GATM, ARHGEF11, NAAA and HSPA2) were selected for further analysis. To further validate the candidate biomarkers, we identified the R phenotype in 31 cattle (9 R and 22 S). Macrophage mRNA was isolated from this group of cattle. Expression of candidate biomarkers was evaluated by qPCR 2^-ΔCt and ROC curves to determine their diagnostic potential. Candidates IL1RN and ARHGEF11 discriminates between R and S cattle. Furthermore, combination of candidates ARHGEF11: VNN: HSPA2 discriminate between R from S with AUC 0.7993 and agreement index of 0.853 (p ≤ 0.01). Our data suggest that candidate biomarkers may support the preliminary screening to identify natural resistance in herds against Mycobacterium bovis in Holstein-Friesian cattle.

Biomarkers of health and welfare: A One Health perspective from the laboratory side

A biomarker is any measurement reflecting an interaction between a biological system and a potential hazard, which may be chemical, physical, or biological. The One World, One Health concept established that human and animal health and the environmental state are highly interconnected, sharing common aspects that can be applied globally in these three components. In this paper, we review how the concept of One Health can be applied to biomarkers of health and welfare, with a special focus on five points that can be applied to any biomarker when it is expected to be used to evaluate the human, animal or environmental health. Three of these points are: (1) the different biomarkers that can be used, (2) the different sample types where the biomarkers can be analysed, and (3) the main methods that can be used for their measurement. In addition, we will evaluate two key points needed for adequate use of a biomarker in any situation: (4) a proper analytical validation in the sample that it is going to be used, and (5) a correct selection of the biomarker. It is expected that this knowledge will help to have a broader idea about the use of biomarkers of health and welfare and also will contribute to a better and more accurate use of these biomarkers having in mind their One Health perspective.

Integrated bioinformatics analysis and validation revealed potential immune-regulatory miR-892b, miR-199b-5p and miR-582-5p as diagnostic biomarkers in active tuberculosis

Tuberculosis (TB) is one of the most prevalent pulmonary diseases caused by Mycobacterium tuberculosis (Mtb). MiRNAs (miRNAs) participate in TB progression by modulating the host-pathogen interaction. Bioinformatics advancements provide basis for exploring novel immunoregulatory miRNAs and their performance as diagnostic biomarkers. Gene and miRNA expression datasets, GSE29190 and GSE54992, were downloaded from Gene Expression Omnibus (GEO) database. Based on fold changes and statistical significance, a total of 7463 differentially expressed mRNAs (DE-mRNAs) and 38 differentially expressed miRNAs (DE-miRNAs) were screened. Function annotation and protein-protein interaction (PPI) network were constructed to reveal underlying mechanisms of TB pathogenesis. Functional annotation identified the MAPK signalling pathway and leukocyte migration as the top enriched processes. The PPI and MGIP networks indicated that chemokine ligands like CXCL1/CXCL2 and receptors, like CCR7 were important down-regulated genes, implying that a protective mechanism against overdue inflammation induced cell death. MiRNA-gene-immune processes (MGIP) network enriched 7 deregulated miRNAs, and their expression was further examined with quantitative real-time PCR (qRT-PCR), in PBMC samples of 20 active TB patients and 20 healthy donors. The diagnostic performance was evaluated with ROC curves. MiR-892b; miR-199b-5p and miR-582-5p were significantly deregulated in TB patients, compared with healthy participants. The best overall performance was from miR-892b, with an area under curve (AUC) of 0.77, 55% sensitivity and 90% specificity. AUC of miR-199b-5p and miR-582-5p were 0.71 and 0.70, respectively. MiR-892b, miR-199b-5p and miR-582-5p could be considered promising novel diagnostic biomarkers for active tuberculosis.

Global Transcriptome Analysis During Adipogenic Differentiation and Involvement of Transthyretin Gene in Adipogenesis in Cattle

Adipose tissue plays central role in determining the gustatory quality of beef, but traditional Chinese beef cattle have low levels of fat content. We applied RNA-seq to study the molecular mechanisms underlying adipocyte differentiation in Qinchuan cattle. A total of 18,283 genes were found to be expressed in preadipocytes and mature adipocytes, respectively. 470 of which were significantly differentially expressed genes (DEGs) [false discovery rate (FDR) values < 0.05 and fold change ≥ 2]. In addition, 4534 alternative splicing (AS) events and 5153 AS events were detected in preadipocytes and adipocytes, respectively. We constructed a protein interaction network, which suggested that collagen plays an important role during bovine adipogenic differentiation. We characterized the function of the most down-regulated DEG (P < 0.001) among genes we have detected by qPCR, namely, the transthyretin (TTR) gene. Overexpression of TTR appears to promote the expression of the peroxisome proliferator activated receptor γ (PPARγ) (P < 0.05) and fatty acid binding Protein 4 (FABP4) (P < 0.05). Hence, TTR appears to be involved in the regulation of bovine adipogenic differentiation. Our study represents the comprehensive approach to explore bovine adipocyte differentiation using transcriptomic data and reports an involvement of TTR during bovine adipogenic differentiation. Our results provide novel insights into the molecular mechanisms underlying bovine adipogenic differentiation.

RNA Sequencing (RNA-Seq) Reveals Extremely Low Levels of Reticulocyte-Derived Globin Gene Transcripts in Peripheral Blood From Horses (Equus caballus) and Cattle (Bos taurus)

RNA-seq has emerged as an important technology for measuring gene expression in peripheral blood samples collected from humans and other vertebrate species. In particular, transcriptomics analyses of whole blood can be used to study immunobiology and develop novel biomarkers of infectious disease. However, an obstacle to these methods in many mammalian species is the presence of reticulocyte-derived globin mRNAs in large quantities, which can complicate RNA-seq library sequencing and impede detection of other mRNA transcripts. A range of supplementary procedures for targeted depletion of globin transcripts have, therefore, been developed to alleviate this problem. Here, we use comparative analyses of RNA-seq data sets generated from human, porcine, equine, and bovine peripheral blood to systematically assess the impact of globin mRNA on routine transcriptome profiling of whole blood in cattle and horses. The results of these analyses demonstrate that total RNA isolated from equine and bovine peripheral blood contains very low levels of globin mRNA transcripts, thereby negating the need for globin depletion and greatly simplifying blood-based transcriptomic studies in these two domestic species.

Assessing genetic architecture and signatures of selection of dual purpose Gir cattle populations using genomic information

Gir is one of the main cattle breeds raised in tropical South American countries. Strong artificial selection through its domestication resulted in increased genetic differentiation among the countries in recent years. Over the years, genomic studies in Gir have become more common. However, studies of population structure and signatures of selection in divergent Gir populations are scarce and need more attention to better understand genetic differentiation, gene flow, and genetic distance. Genotypes of 173 animals selected for growth traits and 273 animals selected for milk production were used in this study. Clear genetic differentiation between beef and dairy populations was observed. Different criteria led to genetic divergence and genetic differences in allele frequencies between the two populations. Gene segregation in each population was forced by artificial selection, promoting isolation, and increasing genetic variation between them. Results showed evidence of selective forces in different regions of the genome. A total of 282 genes were detected under selection in the test population based on the fixation index (Fst), integrated haplotype score (iHS), and cross-population extend haplotype homozygosity (XP-EHH) approaches. The QTL mapping identified 35 genes associated with reproduction, milk composition, growth, meat and carcass, health, or body conformation traits. The investigation of genes and pathways showed that quantitative traits associated to fertility, milk production, beef quality, and growth were involved in the process of differentiation of these populations. These results would support further investigations of population structure and differentiation in the Gir breed.

Differential gene expression in Mycobacterium bovis challenged monocyte-derived macrophages of cattle

A functional genomics approach was used to examine the immune response for transcriptional profiling of PBMC M. bovis infected cattle and healthy control cattle to stimulation with bovine tuberculin (purified protein derivative PPD-b). Total cellular RNA was extracted from non-challenged control and M. bovis challenged MDM for all animals at intervals of 6 h post-challenge, in response to in-vitro challenge with M. bovis (multiplicity of infection 2:1) and prepared for global gene expression analysis using the Agilent Bovine (V2) Gene Expression Microarray, 8 × 60 K. The pattern of expression of these genes in PPD bovine stimulated PBMC provides the first description of an M. bovis specific signature of infection that may provide insights into the molecular basis of the host response to infection. Analysis of these mapped reads showed 2450 genes (1291 up regulated and 1158 down regulated) 462 putative natural antisense transcripts (354 up-regulated and 108 down regulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). The results provided enrichment for genes involved top ten up regulated and down regulated panel of genes, including transcription factors proliferation of T and B lymphocytes. The highest differentially-expressed genes were associated to immune and inflammatory responses, immunity, differentiation, cell growth, apoptosis, cellular trafficking and regulation of lipolysis and thermogenesis. Microarray results were confirmed in infected cattle by RT qPCR to identify potential biomarkers TLR2, CD80, NFKB1, IL8, CXCL6 and ADORA3 of bovine tuberculosis.

In vitro gene expression profile of bovine peripheral blood mononuclear cells in early Mycobacterium bovis infection

The intracellular parasite Mycobacterium bovis (M. bovis) causes tuberculosis in cattle and humans. Understanding the interactions between M. bovis and host cells is essential in developing tools for the prevention, detection, and treatment of M. bovis infection. Gene expression profiles provide a large amount of information regarding the molecular mechanisms underlying these interactions. The present study analyzed changes in gene expression in bovine peripheral blood mononuclear cells (PBMCs) at 0, 4 and 24 h following exposure to M. bovis. Using bovine whole-genome microarrays, a total of 420 genes were identified that exhibited significant alterations in expression (≥2-fold). Significantly enriched genes were identified using the Kyoto Encyclopedia of Genes and Genomes database, of which the highest differentially expressed genes were associated with the immune system, signal transduction, endocytosis, cellular transport, inflammation, and apoptosis. Of the genes associated with the immune system, 84.85% displayed downregulation. These findings support the view that M. bovis inhibits signaling pathways of antimycobacterial host defense in bovine PBMCs. These in vitro data demonstrated that molecular alterations underlying the pathogenesis of tuberculosis begin early, during the initial 24 h following M. bovis infection.

Whole Blood Transcriptome Analysis of Mycoplasma mycoides Subsp. mycoides-Infected Cattle Confirms Immunosuppression but Does Not Reflect Local Inflammation

Contagious bovine pleuropneumonia (CBPP), caused by Mycoplasma mycoides subsp. mycoides (Mmm), is a severe respiratory disease of cattle responsible for major economic losses in sub-Saharan Africa. Disease control relies mainly on the use of empirically attenuated vaccines that provide limited protection. Thus, understanding the virulence mechanisms used by Mmm as well as the role of the host immune system in disease development, persistence, and control is a prerequisite for the development of new, rationally designed control strategies. The aim of this study was to assess the use of whole blood transcriptome analysis to study cattle-Mmm interactions, starting by the characterization of the bovine response to Mmm infection during the acute form of the disease. For that purpose, we compared the transcriptome profile of whole blood from six cattle, before challenge by contact with Mmm-infected animals and at the appearance of first clinical signs, using a bovine microarray. Functional analysis revealed that 680 annotated genes were differentially expressed, with an overwhelming majority of down-regulated genes characterizing an immunosuppression. The main bio-functions affected were "organismal survival", "cellular development, morphology and functions" and "cell-to cell signaling and interactions". These affected functions were consistent with the results of previous in vitro immunological studies. However, microarray and qPCR validation results did not highlight pro-inflammatory molecules (such as TNFα, TLR2, IL-12B and IL-6), whereas inflammation is one of the most characteristic traits of acute CBPP. This global gene expression pattern may be considered as the result, in blood, of the local pulmonary response and the systemic events occurring during acute CBPP. Nevertheless, to understand the immune events occurring during disease, detailed analyses on the different immune cell subpopulations, either in vivo, at the local site, or in vitro, will be required. Whole blood transcriptome analysis remains an interesting approach for the identification of bio-signatures correlating to recovery and protection, which should facilitate the evaluation and validation of novel vaccine formulations.

Transcriptome changes upon in vitro challenge with Mycobacterium bovis in monocyte-derived macrophages from bovine tuberculosis-infected and healthy cows

As innate immune cells, macrophages are expected to respond to mycobacterial infection equally in both Mycobacterium bovis-infected cows and healthy cows. We previously found that monocyte-derived macrophages (MDMs) from M. bovis-infected cows respond differently than MDMs from healthy cows when exposed to in vitro M. bovis challenge. We have now used the Agilent™ Bovine Gene Expression Microarray to examine transcriptional differences between these MDMs. At a high multiplicity of infection (10), in vitro challenge led to changes in several thousands of genes, with dysregulation at multiple orders of magnitude. For example, significant changes were seen for colony stimulating factor 3 (granulocyte) (CSF3), colony stimulating factor 2 (granulocyte-macrophage) (CSF2), and chemokine (C-C motif) ligand 20 (CCL20). Classical macrophage activation was also observed, although to a lesser degree in interleukin 12 (IL12) expression. For macrophages, kallikrein-related peptidase 12 (KLK12) and protease, serine, 2 (trypsin 2) (PRSS2), as well as a secreted protein, acidic, cysteine-rich (osteonectin) (SPARC)-centered matricellular gene network, were differentially expressed in infected animals. Finally, global transcriptome fold-changes caused by in vitro challenge were higher in healthy cows than in tuberculosis-positive cows, suggesting that healthy macrophages responded marginally better to in vitro infection. Macrophages from healthy and already infected animals can both be fully activated during M. bovis infection, yet there are differences between these macrophages: distinct expression pattern in matricellular proteins, and their different responses to in vitro infection.

Role of the Mce1 transporter in the lipid homeostasis of Mycobacterium tuberculosis

Tuberculosis is one of the leading causes of mortality throughout the world. Mycobacterium tuberculosis, the causative agent of human tuberculosis, has developed several strategies involving proteins and other compounds known collectively as virulence factors to subvert human host defences and invade the human host. The Mce proteins are among these virulence-related proteins and are encoded by the mce1, mce2, mce3 and mce4 operons in the genome of M. tuberculosis. It has been proposed that these operons encode ABC-like lipid transporters; however, the nature of their substrates has only been revealed in the case of the Mce4 proteins. Here we found that the knockout of the mce1 operon alters the lipid profile of M. tuberculosis H37Rv and the uptake of palmitic acid. Thin layer chromatography and liquid chromatography-mass spectrometry analysis showed that the mce1 mutant accumulates more mycolic acids than the wild type and complemented strains. Interestingly, this accumulation of mycolic acid is exacerbated when bacteria are cultured in the presence of palmitic acid or arachidonic acid. These results suggest that the mce1 operon may serve as a mycolic acid re-importer.

Study of the in vivo role of Mce2R, the transcriptional regulator of mce2 operon in Mycobacterium tuberculosis

Background

Tuberculosis is one of the leading causes of mortality throughout the world. Mycobacterium tuberculosis, the agent of human tuberculosis, has developed strategies involving proteins and other compounds called virulence factors to subvert human host defences and damage and invade the human host. Among these virulence-related proteins are the Mce proteins, which are encoded in the mce1, mce2, mce3 and mce4 operons of M. tuberculosis. The expression of the mce2 operon is negatively regulated by the Mce2R transcriptional repressor. Here we evaluated the role of Mce2R during the infection of M. tuberculosis in mice and macrophages and defined the genes whose expression is in vitro regulated by this transcriptional repressor.

Results

We used a specialized transduction method for generating a mce2R mutant of M. tuberculosis H37Rv. Although we found equivalent replication of the MtΔmce2R mutant and the wild type strains in mouse lungs, overexpression of Mce2R in the complemented strain (MtΔmce2RComp) significantly impaired its replication. During in vitro infection of macrophages, we observed a significantly increased association of the late endosomal marker LAMP-2 to MtΔmce2RComp-containing phagosomes as compared to MtΔmce2R and the wild type strains. Whole transcriptional analysis showed that Mce2R regulates mainly the expression of the mce2 operon, in the in vitro conditions studied.

Conclusions

The findings of the current study indicate that Mce2R weakly represses the in vivo expression of the mce2 operon in the studied conditions and argue for a role of the proteins encoded in Mce2R regulon in the arrest of phagosome maturation induced by M. tuberculosis.

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro

Background

Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA.

Results

A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology.

Conclusions

This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA.