Global gene expression and systems biology analysis of bovine monocyte-derived macrophages in response to in vitro challenge with Mycobacterium bovis

A bovine pulmosphere model and multiomics reveal early host response signature in tuberculosis

Early interactions between tubercle bacilli and lung cells are critical in tuberculosis (TB) pathogenesis. Conventional two-dimensional cell cultures fail to replicate the multicellular complexity of lungs. We introduce a three-dimensional pulmosphere model for Mycobacterium tuberculosis infection in bovine systems, demonstrating through comprehensive transcriptome and proteome analyses that these multicellular spheroids closely mimic lung cell diversity, interactions, and extracellular matrix (ECM) composition. Cell viability, hypoxia, and reactive oxygen species assessments over three weeks confirm the model's suitability. To establish infection, we employed M. bovis BCG-an attenuated vaccine strain, and M. tuberculosis H37Rv-a laboratory adapted human clinical strain that is attenuated for cattle infection compared to M. bovis. Both infection upregulated key host pathways; however, M. tuberculosis induced distinct responses, including enhanced ECM receptors expression, neutrophil chemotaxis, interferon signaling, and RIG-1 signaling. A six genes/protein signature- IRF1, CCL5, CXCL8, CXCL10, SERPINE1, and CFB -emerges as an early host response marker to M. tuberculosis infection. Infection with virulent M. bovis and M. orygis revealed a shared upregulated gene signature across Mycobacterium tuberculosis complex species, but with pathogen-specific variations. This study presents a robust ex vivo bovine pulmosphere TB model with implications in biomarkers discovery, high-throughput drug screening, and TB control strategies.

Screening of mRNA markers in early bovine tuberculosis blood samples

Bovine tuberculosis (bTB) is a chronic zoonotic disease caused by Mycobacterium bovis. A large number of cattle are infected with bTB every year, resulting in huge economic losses. How to control bTB is an important issue in the current global livestock economy. In this study, the original transcriptome sequences related to this study were obtained from the dataset GSE192537 by searching the Gene Expression Omnibus (GEO) database. Our differential gene analysis showed that there were obvious biological activities related to immune activation and immune regulation in the early stage of bTB. Immune-related biological processes were more active in the early stage of bTB than in the late. There were obvious immune activation and immune cell recruitment in the early stage of bTB. Regulations in immune receptors are associated with pathophysiological processes of the early stage of bTB. A gene module consisting of 236 genes significantly related to the early stage of bTB was obtained by weighted gene co-expression network analysis, and 18 hub genes were further identified as potential biomarkers or therapeutic targets. Finally, by random forest algorithm and logistic regression modeling, FCRL1 was identified as a representative mRNA marker in early bTB blood. FCRL1 has the potential to be a diagnostic biomarker in early 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.

Mycobacterium bovis induced human tuberculosis in India: Current status, challenges & opportunities

Tuberculosis (TB) caused by Mycobacterium tuberculosis is a leading cause of human deaths due to any infectious disease worldwide. However, infection of Mycobacterium bovis, primarily an animal pathogen, also leads to the development of 'human tuberculosis'. Infected animals have been considered the major source of M. bovis infection and humans get exposed to M. bovis through close contact with infected animals or consumption of contaminated milk, unpasteurized dairy products and improperly cooked contaminated meat. The information on the global distribution of bovine TB (bTB) is limited, but the disease has been reported from all the livestock-producing middle- and low-income countries of the world. In recent years, there is a renewed interest for the control of bTB to minimize human infection worldwide. In India, while the sporadic presence of M. bovis has been reported in domestic animals, animal-derived food products and human beings from different geographical regions of the country, the information on the national prevalence of bTB and transmission dynamics of zoonotic TB is, however, not available. The present article reviewed published information on the status of M. bovis-induced zoonotic TB to highlight the key challenges and opportunities for intervention to minimize the risk of M. bovis infection in humans and secure optimum animal productivity in India.

Transcriptional Profiling of Early and Late Phases of Bovine Tuberculosis

Bovine tuberculosis, caused by Mycobacterium tuberculosis var. bovis (M. bovis), is an important enzootic disease affecting mainly cattle, worldwide. Despite the implementation of national campaigns to eliminate the disease, bovine tuberculosis remains recalcitrant to eradication in several countries. Characterizing the host response to M. bovis infection is crucial for understanding the immunopathogenesis of the disease and for developing better control strategies. To profile the host responses to M. bovis infection, we analyzed the transcriptome of whole blood cells collected from experimentally infected calves with a virulent strain of M. bovis using RNA transcriptome sequencing (RNAseq). Comparative analysis of calf transcriptomes at early (8 weeks) vs. late (20 weeks) aerosol infection with M. bovis revealed divergent and unique profile for each stage of infection. Notably, at the early time point, transcriptional upregulation was observed among several of the top-ranking canonical pathways involved in T-cell chemotaxis. At the late time point, enrichment in the cell mediated cytotoxicity (e.g. Granzyme B) was the predominant host response. These results showed significant change in bovine transcriptional profiles and identified networks of chemokine receptors and monocyte chemoattractant protein (CCL) co-regulated genes that underline the host-mycobacterial interactions during progression of bovine tuberculosis in cattle. Further analysis of the transcriptomic profiles identified potential biomarker targets for early and late phases of tuberculosis in cattle. Overall, the identified profiles better characterized identified novel immunomodulatory mechanisms and provided a list of targets for further development of potential diagnostics for tuberculosis in cattle.

Identifying Bacterial and Host Factors Involved in the Interaction of Mycobacterium bovis with the Bovine Innate Immune Cells

Bovine tuberculosis is an important animal and zoonotic disease caused by Mycobacterium bovis. The innate immune response is the first line of defense against pathogens and is also crucial for the development of an efficient adaptive immune response. In this study we used an in vitro co-culture model of antigen presenting cells (APC) and autologous lymphocytes derived from peripheral blood mononuclear cells to identify the cell populations and immune mediators that participate in the development of an efficient innate response capable of controlling the intracellular replication of M. bovis. After M. bovis infection, bovine immune cell cultures displayed upregulated levels of iNOS, IL-22 and IFN-γ and the induction of the innate immune response was dependent on the presence of differentiated APC. Among the analyzed M. bovis isolates, only a live virulent M. bovis isolate induced an efficient innate immune response, which was increased upon stimulation of cell co-cultures with the M. bovis culture supernatant. Moreover, we demonstrated that an allelic variation of the early secreted protein ESAT-6 (ESAT6 T63A) expressed in the virulent strain is involved in this increased innate immune response. These results highlight the relevance of the compounds secreted by live M. bovis as well as the variability among the assessed M. bovis strains to induce an efficient innate immune response.

Crosstalk between Neospora caninum and the bovine host at the maternal-foetal interface determines the outcome of infection

Neospora caninum is an apicomplexan cyst-forming parasite that is considered one of the main causes of abortion. The pathogenic mechanisms associated with parasite virulence at the maternal-foetal interface that are responsible for the outcome of infection are largely unknown. Here, utilizing placentomes from cattle experimentally infected with high-virulence (Nc-Spain7) and low-virulence (Nc-Spain1H) isolates, we studied key elements of the innate and adaptive immune responses, as well as components of the extracellular matrix (ECM), at 10 and 20 days post-infection (dpi). The low-virulence isolate elicited a robust immune response characterized by upregulation of genes involved in pathogen recognition, chemokines and pro-inflammatory cytokines, crucial for its adequate control. In addition, Nc-Spain1H triggered the expression of anti-inflammatory cytokines and other mechanisms implicated in the maintenance of ECM integrity to ensure foetal survival. In contrast, local immune responses were initially (10 dpi) impaired by Nc-Spain7, allowing parasite multiplication. Subsequently (20 dpi), a predominantly pro-inflammatory Th1-based response and an increase in leucocyte infiltration were observed. Moreover, Nc-Spain7-infected placentomes from animals carrying non-viable foetuses exhibited higher expression of the IL-8, TNF-α, iNOS and SERP-1 genes and lower expression of the metalloproteases and their inhibitors than Nc-Spain7-infected placentomes from animals carrying viable foetuses. In addition, profound placental damage characterized by an alteration in the ECM organization in necrotic foci, which could contribute to foetal death, was found. Two different host-parasite interaction patterns were observed at the bovine placenta as representative examples of different evolutionary strategies used by this parasite for transmission to offspring.

Transcriptomic Profiles of Monocyte-Derived Macrophages in Response to Escherichia coli is Associated with the Host Genetics

Reactive Nitrogen Species (RNS) are a group of bactericidal molecules produced by macrophages in response to pathogens in a process called oxidative burst. Nitric oxide (NO^-) is a member of RNS produced from arginine by inducible Nitric Oxide Synthase (iNOS) enzyme. The activity of iNOS and production of NO^- by macrophages following stimulation is one of the indicators of macrophage polarization towards M1/proinflammatory. Production of NO^- by bovine monocyte-derived macrophage (MDM) and mouse peritoneal macrophages has been shown to be strongly associated with host genetic with the heritability of 0.776 in bovine MDM and 0.8 in mouse peritoneal macrophages. However, the mechanism of genetic regulation of macrophage response has remained less explored. In the current study, the transcriptome of bovine MDMs was compared between two extreme phenotypes that had been classified as high and low responder based on NO^- production. The results showed that 179 and 392 genes were differentially expressed (DE) between high and low responder groups at 3 and 18 hours after exposure to Escherichia coli, respectively. A set of 11 Transcription Factors (TFs) (STAT1, IRF7, SPI1, STAT4, IRF1, HIF1A, FOXO3, REL, NFAT5, HIC1, and IRF4) at 3 hours and a set of 13 TFs (STAT1, IRF1, HIF1A, STAT4, ATF4, TP63, EGR1, CDKN2A, RBL1, E2F1, PRDM1, GATA3, and IRF4) at 18 hours after exposure to E. coli were identified to be differentially regulated between the high and low responder phenotypes. These TFs were found to be divided into two clusters of inflammatory- and hypoxia-related TFs. Functional analysis revealed that some key canonical pathways such as phagocytosis, chemotaxis, antigen presentation, and cell-to-cell signalling are enriched among the over-expressed genes by high responder phenotype. Based on the results of this study, it was inferred that the functional characteristics of bovine MDMs are associated with NO-based classification. Since NO^- production is strongly associated with host genetics, this study for the first time shows the distinct proinflammatory profiles of macrophages are controlled by the natural genetic polymorphism in an outbred population. In addition, the results suggest that genetics can be considered as a new dimension in the current model of macrophage polarization which is currently described by the combination of stimulants, only.

Bovine Derived in vitro Cultures Generate Heterogeneous Populations of Antigen Presenting Cells

Antigen presenting cells (APC) of the mononuclear phagocytic system include dendritic cells (DCs) and macrophages (Macs) which are essential mediators of innate and adaptive immune responses. Many of the biological functions attributed to these cell subsets have been elucidated using models that utilize in vitro-matured cells derived from common progenitors. However, it has recently been shown that monocyte culture systems generate heterogeneous populations of cells, DCs, and Macs. In light of these findings, we analyzed the most commonly used bovine in vitro-derived APC models and compared them to bona fide DCs. Here, we show that bovine monocyte-derived DCs and Macs can be differentiated on the basis of CD11c and MHC class II (MHCII) expression and that in vitro conditions generate a heterologous group of both DCs and Macs with defined and specific biological activities. In addition, skin-migrating macrophages present in the bovine afferent lymph were identified and phenotyped for the first time. RNA sequencing analyses showed that these monophagocytic cells have distinct transcriptomic profiles similar to those described in other species. These results have important implications for the interpretation of data obtained using in vitro systems.

Stimulation of bovine monocyte-derived macrophages with lipopolysaccharide, interferon-ɣ, Interleukin-4 or Interleukin-13 does not induce detectable changes in nitric oxide or arginase activity

BACKGROUND

Bacterial lipopolysaccharide and interferon-γ stimulation of rodent macrophages in vitro induces up-regulation of inducible nitric oxide synthase, whereas interleukin-4 stimulation results in increased activity of arginase-1. Thus different stimulants result in differing macrophage phenotypes, appropriate for responses to a range of pathogens. The current study was conducted in order to determine whether bovine macrophages derived from monocytes and spleen respond similarly.

RESULTS

Lipopolysaccharide and interferon-γ did not induce detectable increases in nitric oxide production by bovine monocyte-derived or splenic macrophages in vitro. Similarly, interleukin-4 and interleukin-13 did not affect arginase activity. However, changes in transcription of genes coding for these products were detected.

CONCLUSION

Differences between macrophage activation patterns exist between cattle and other species and these differences may occur during the post-transcription phase.

Bacterial and viral pathogen-associated molecular patterns induce divergent early transcriptomic landscapes in a bovine macrophage cell line

Background

Pathogens stimulate immune functions of macrophages. Macrophages are a key sentinel cell regulating the response to pathogenic ligands and orchestrating the direction of the immune response. Our study aimed at investigating the early transcriptomic changes of bovine macrophages (Bomacs) in response to stimulation with CpG DNA or polyI:C, representing bacterial and viral ligands respectively, and performed transcriptomics by RNA sequencing (RNASeq). KEGG, GO and IPA analytical tools were used to reconstruct pathways, networks and to map out molecular and cellular functions of differentially expressed genes (DE) in stimulated cells.

Results

A one-way ANOVA analysis of RNASeq data revealed significant differences between the CpG DNA and polyI:C-stimulated Bomac. Of the 13,740 genes mapped to the bovine genome, 2245 had p-value ≤0.05, deemed as DE. At 6 h post stimulation of Bomac, poly(I:C) induced a very different transcriptomic profile from that induced by CpG DNA. Whereas, 347 genes were upregulated and 210 downregulated in response to CpG DNA, poly(I:C) upregulated 761 genes and downregulated 414 genes. The topmost DE genes in poly(I:C)-stimulated cells had thousand-fold changes with highly significant p-values, whereas in CpG DNA stimulated cells had 2–5-fold changes with less stringent p-values. The highest DE genes in both stimulations belonged to the TNF superfamily, TNFSF18 (CpG) and TNFSF10 (poly(I:C)) and in both cases the lowest downregulated gene was CYP1A1. CpG DNA highly induced canonical pathways that are unrelated to immune response in Bomac. CpG DNA influenced expression of genes involved in molecular and cellular functions in free radical scavenging. By contrast, poly(I:C) highly induced exclusively canonical pathways directly related to antiviral immune functions mediated by interferon signalling genes. The transcriptomic profile after poly(I:C)-stimulation was consistent with induction of TLR3 signalling.

Conclusion

CpG DNA and poly(I:C) induce different early transcriptional landscapes in Bomac, but each is suited to a specific function of macrophages during interaction with pathogens. Poly(I:C) influenced antiviral response genes, whereas CpG DNA influenced genes important for phagocytic processes. Poly(I:C) was more potent in setting the inflammatory landscape desirable for an efficient immune response against virus infection.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5411-5) contains supplementary material, which is available to authorized users.

Identification of genes associated with susceptibility to Mycobacterium avium ssp. paratuberculosis (Map) tissue infection in Holstein cattle using gene set enrichment analysis-SNP

Multiple genome-wide association analyses have investigated susceptibility to bovine paratuberculosis, but few loci have been identified across independent cattle populations. A SNP-based gene set enrichment analysis (GSEA-SNP) allows expanded identification of genes with moderate effects on a trait through the enrichment of gene sets instead of identifying only few loci with large effects. Therefore, the objective of this study was to identify genes that were moderately associated with Mycobacterium avium ssp. paratuberculosis (Map) tissue infection using GSEA-SNP in Holstein cattle from the Pacific Northwest (PNW; n = 205) and from the PNW and Northeast (PNW+NE; n = 245) which were previously genotyped with the Illumina BovineSNP50 BeadChip. The GSEA-SNP utilized 4389 gene sets from five databases. For each annotated gene in the UMD3.1 assembly (n = 19,723), the most significant SNP within each gene and its surrounding region (10 kb up- and downstream) was selected as a proxy for that gene. Any gene set with a normalized enrichment score > 2.5 was considered enriched. Thirteen gene sets (8 PNW GSEA-SNP; 5 PNW+NE) were enriched in these analyses and all have functions that relate to nuclear factor kappa beta. Nuclear factor kappa beta is critical to gut immune responses, implicated in host immune responses to other mycobacterial diseases, and has established roles in inflammation as well as cancer. Gene sets and genes moderately associated with Map infection could be used in genomic selection to allow producers to select for less susceptible cattle, lower the prevalence of the disease, and reduce economic losses.

The immunoregulatory effects of co-infection with Fasciola hepatica: From bovine tuberculosis to Johne's disease

Fasciola hepatica (liver fluke) is a parasite prevalent in much of the world that causes the economically-important disease of fasciolosis in livestock. The threat that this disease poses extends beyond its direct effects due to the parasite's immunomodulatory effects. Research at this laboratory is focusing on whether this immunoregulation can, in animals infected with liver fluke, exert a bystander effect on concurrent infections in the host. It has already been established that F. hepatica infection reduces cell mediated immune responses to Mycobacterium bovis in cattle, and that the interaction between the two pathogens can be detected on an epidemiological scale. This review explores the immunological consequences of co-infection between F. hepatica and other bacterial infections. Arguments are presented suggesting that immunity of cattle to Mycobacterium avium subsp. paratuberculosis is also likely to be affected.

The CD4(+) T cell methylome contributes to a distinct CD4(+) T cell transcriptional signature in Mycobacterium bovis-infected cattle

We hypothesised that epigenetic regulation of CD4⁺ T lymphocytes contributes to a shift toward a dysfunctional T cell phenotype which may impact on their ability to clear mycobacterial infection. Combined RNA-seq transcriptomic profiling and Reduced Representation Bisulfite Sequencing identified 193 significantly differentially expressed genes and 760 differentially methylated regions (DMRs), between CD4⁺ T cells from M. bovis infected and healthy cattle. 196 DMRs were located within 10 kb of annotated genes, including GATA3 and RORC, both of which encode transcription factors that promote T_H2 and T_H17 T helper cell subsets respectively. Gene-specific DNA methylation and gene expression levels for the TNFRSF4 and Interferon-γ genes were significantly negatively correlated suggesting a regulatory relationship. Pathway analysis of DMRs identified enrichment of genes involved in the anti-proliferative TGF-β signaling pathway and TGFB1 expression was significantly increased in peripheral blood leukocytes from TB-infected cattle. This first analysis of the bovine CD4⁺ T cell methylome suggests that DNA methylation directly contributes to a distinct gene expression signature in CD4⁺ T cells from cattle infected with M. bovis. Specific methylation changes proximal to key inflammatory gene loci may be critical to the emergence of a non-protective CD4⁺ T cell response during mycobacterial infection in cattle.

Differential Transcriptome Networks between IDO1-Knockout and Wild-Type Mice in Brain Microglia and Macrophages

Microglia in the brain and macrophages in peripheral organs are cell types responsible for immune response to challenges. Indoleamine 2,3-dioxygenase 1 (IDO1) is an immunomodulatory enzyme of the tryptophan pathway that is expressed in the brain. The higher activity of IDO1 in response to immune challenge has been implicated in behavioral disorders. The impact of IDO1 depletion on the microglia transcriptome has not been studied. An investigation of the transcript networks in the brain microglia from IDO1-knockout (IDO1-KO) mice was undertaken, relative to peripheral macrophages and to wild-type (WT) mice under unchallenged conditions. Over 105 transcript isoforms were differentially expressed between WT and IDO1-KO within cell type. Within microglia, Saa3 and Irg1 were over-expressed in IDO1-KO relative to WT. Within macrophages, Csf3 and Sele were over-expressed in IDO1-KO relative to WT. Among the genes differentially expressed between strains, enriched biological processes included ion homeostasis and ensheathment of neurons within microglia, and cytokine and chemokine expression within macrophages. Over 11,110 transcript isoforms were differentially expressed between microglia and macrophages and of these, over 10,800 transcripts overlapped between strains. Enriched biological processes among the genes over- and under-expressed in microglia relative to macrophages included cell adhesion and apoptosis, respectively. Detected only in microglia or macrophages were 421 and 43 transcript isoforms, respectively. Alternative splicing between cell types based on differential transcript isoform abundance was detected in 210 genes including Phf11d, H2afy, and Abr. Across strains, networks depicted a predominance of genes under-expressed in microglia relative to macrophages that may be a precursor for the different response of both cell types to challenges. The detected transcriptome differences enhance the understanding of the role of IDO1 in the microglia transcriptome under unchallenged conditions.

Dual Transcriptome Profiling of Leishmania-Infected Human Macrophages Reveals Distinct Reprogramming Signatures

Macrophages are mononuclear phagocytes that constitute a first line of defense against pathogens. While lethal to many microbes, they are the primary host cells of Leishmania spp. parasites, the obligate intracellular pathogens that cause leishmaniasis. We conducted transcriptomic profiling of two Leishmania species and the human macrophage over the course of intracellular infection by using high-throughput RNA sequencing to characterize the global gene expression changes and reprogramming events that underlie the interactions between the pathogen and its host. A systematic exclusion of the generic effects of large-particle phagocytosis revealed a vigorous, parasite-specific response of the human macrophage early in the infection that was greatly tempered at later time points. An analogous temporal expression pattern was observed with the parasite, suggesting that much of the reprogramming that occurs as parasites transform into intracellular forms generally stabilizes shortly after entry. Following that, the parasite establishes an intracellular niche within macrophages, with minimal communication between the parasite and the host cell later during the infection. No significant difference was observed between parasite species transcriptomes or in the transcriptional response of macrophages infected with each species. Our comparative analysis of gene expression changes that occur as mouse and human macrophages are infected by Leishmania spp. points toward a general signature of the Leishmania-macrophage infectome.

Little is known about the transcriptional changes that occur within mammalian cells harboring intracellular pathogens. This study characterizes the gene expression signatures of Leishmania spp. parasites and the coordinated response of infected human macrophages as the pathogen enters and persists within them. After accounting for the generic effects of large-particle phagocytosis, we observed a parasite-specific response of the human macrophages early in infection that was reduced at later time points. A similar expression pattern was observed in the parasites. Our analyses provide specific insights into the interplay between human macrophages and Leishmania parasites and constitute an important general resource for the study of how pathogens evade host defenses and modulate the functions of the cell to survive intracellularly.

Neonatal lipopolysaccharide exposure does not diminish the innate immune response to a subsequent lipopolysaccharide challenge in Holstein bull calves

The innate immune response following experimental mastitis is quite variable between individual dairy cattle. An inflammatory response that minimizes collateral damage to the mammary gland while still effectively resolving the infection following pathogen exposure is beneficial to dairy producers. The ability of a lipopolysaccharide (LPS) exposure in early life to generate a low-responding phenotype and thus reduce the inflammatory response to a later-life LPS challenge was investigated in neonatal bull calves. Ten Holstein bull calves were randomly assigned to either an early life LPS (ELL) group (n=5) or an early life saline (ELS) group (n=5). At 7d of age, calves received either LPS or saline, and at 32d of age, all calves were challenged with an intravenous dose of LPS to determine the effect of the early life treatment (LPS or saline) on the immune response generated toward a subsequent LPS challenge. Dermal fibroblast and monocyte-derived macrophage cultures from each calf were established at age 20 and 27d, respectively, to model sustained effects from the early life LPS exposure on gene expression and protein production of components within the LPS response pathway. The ELL calves had greater levels of plasma IL-6 and tumor necrosis factor-α than the ELS calves following the early life LPS or saline treatments. However, levels of these 2 immune markers were similar between ELL and ELS calves when both groups were subsequently challenged with LPS. A comparison of the in vitro LPS responses of the ELL and ELS calves revealed similar patterns of protein production and gene expression following an LPS challenge of both dermal fibroblast and monocyte-derived macrophage cultures established from the treatment groups. Whereas an early life exposure to LPS did not result in a dampened inflammatory response toward a later LPS challenge in these neonatal bull calves, the potential that exposure to inflammation or stress in early life or in utero can create an offspring with a low-responding phenotype as an adult is intriguing and has been documented in rodents. Further work is needed to determine if an inflammatory exposure in utero in a dairy animal would result in a low-responding innate immune phenotype.

GOexpress: an R/Bioconductor package for the identification and visualisation of robust gene ontology signatures through supervised learning of gene expression data

Background

Identification of gene expression profiles that differentiate experimental groups is critical for discovery and analysis of key molecular pathways and also for selection of robust diagnostic or prognostic biomarkers. While integration of differential expression statistics has been used to refine gene set enrichment analyses, such approaches are typically limited to single gene lists resulting from simple two-group comparisons or time-series analyses. In contrast, functional class scoring and machine learning approaches provide powerful alternative methods to leverage molecular measurements for pathway analyses, and to compare continuous and multi-level categorical factors.

Results

We introduce GOexpress, a software package for scoring and summarising the capacity of gene ontology features to simultaneously classify samples from multiple experimental groups. GOexpress integrates normalised gene expression data (e.g., from microarray and RNA-seq experiments) and phenotypic information of individual samples with gene ontology annotations to derive a ranking of genes and gene ontology terms using a supervised learning approach. The default random forest algorithm allows interactions between all experimental factors, and competitive scoring of expressed genes to evaluate their relative importance in classifying predefined groups of samples.

Conclusions

GOexpress enables rapid identification and visualisation of ontology-related gene panels that robustly classify groups of samples and supports both categorical (e.g., infection status, treatment) and continuous (e.g., time-series, drug concentrations) experimental factors. The use of standard Bioconductor extension packages and publicly available gene ontology annotations facilitates straightforward integration of GOexpress within existing computational biology pipelines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0971-3) contains supplementary material, which is available to authorized users.

Relationship of bovine TNF-α gene polymorphisms with the risk of bovine tuberculosis in Holstein cattle

Many studies suggest significant genetic variation in the resistance of cattle and humans to infection with Mycobacterium bovis (M. bovis), the causative agent of zoonotic tuberculosis. TNF-α promotes inflammation and induces apoptosis in response to mycobacterial infection. The aim of the present study was to investigate the influence of single nucleotide polymorphisms of the TNF-α gene on bovine tuberculosis (bTB) susceptibility. We genotyped the TNF-α gene in 74 bTB-infected Holstein cows and 90 healthy control animals. The influence in the exon 3 region of TNF-α polymorphisms on bTB susceptibility was subsequently investigated by association analysis. Our finding demonstrated that the g.27534932A>C polymorphism of the TNF-α is associated with bTB in Holstein cattle. The susceptibility of cattle with the g.27534932A>C genotype compared with the CC genotype was 4.11-fold (95% CI, 1.27–13.36; P=0.02) higher. The g.27534932A>C polymorphism located in exon 3 of the TNF-α gene, and the functional consequence was missense. The deduced amino acid sequence for the protein product revealed an arginine to serine conversion at position 159, which may affect initiation of protein synthesis and disrupt normal TNF-α function that protects animals against mycobacterial infection. A significant association was observed with the A allele as a risk factor for bTB susceptibility (OR, 3.84; 95% CI, 1.21–12.17; P=0.02). In conclusion, this is the first report showing that the g.27534932A>C polymorphism may contribute to TNF-α-mediated bTB susceptibility.

Genomic scan for identifying candidate genes for paratuberculosis resistance in sheep

Breeding objectives relating to health, functional traits and welfare need to receive priority in the research programs and selection schemes, but very few reports are available on natural resistant genotypes in livestock, where some important diseases cause severe economic losses and pose serious zoonotic threats. In this study, diagnosis of paratuberculosis was performed on 759 adult sheep, from a single flock, with the serum antibody enzyme-linked immunosorbent assay; 100 sheep were selected among the extreme divergent animals for the S/P ratio obtained from the test, and were genotyped on the Illumina Ovine SNP50K BeadChip. A genome-wide scan was then performed on the individual marker genotypes, in the attempt to identify genomic regions associated with disease resistance in sheep. For each marker, the allelic substitution effect was calculated by regressing the S/P value on the number of copies of the reference allele. The position on the OARv3.1 Genome Assembly was searched for 32 markers, which showed a statistically significant allelic substitution effect (Raw P < 0.0006 and FDR P < 0.09). All markers were located within, or close to, annotated genes. Five of these genes, SEMA3, CD109, PCP4, PRDM2 and ITFG2 are referred in literature to play a role in either disease resistance or cell-mediated immune response.

High-throughput transcriptomics reveals common and strain-specific responses of human macrophages to infection with Mycobacterium abscessus Smooth and Rough variants

Background

Mycobacterium abscessus (MAB) is an emerging pathogen causing pulmonary infections in those with inflammatory lung disorders, such as Cystic Fibrosis (CF), and is associated with the highest fatality rate among rapidly growing mycobacteria (RGM). Phenotypically, MAB manifests as either a Smooth (MAB-S) or a Rough (MAB-R) morphotype, which differ in their levels of cell wall glycopeptidolipids (GPLs) and in their pathogenicity in vivo. As one of the primary immune cells encountered by MAB, we sought to examine the early transcriptional events within macrophages, following infection with both MAB-S or MAB-R.

Results

We sampled the transcriptomes (mRNA and miRNA) of THP-1-derived macrophages infected with both MAB-R and MAB-S at 1, 4 and 24 h post-infection (hpi) using RNA-seq. A core set of 606 genes showed consistent expression profiles in response to both morphotypes, corresponding to the early transcriptional response to MAB. The core response is type I Interferon (IFN)-driven, involving the NF-κB and MAPK signaling pathways with concomitant pro-inflammatory cytokine production, and network analysis identified STAT1, EGR1, and SRC as key hub and bottleneck genes. MAB-S elicited a more robust transcriptional response at both the mRNA and miRNA levels, which was reflected in higher cytokine levels in culture supernatants. The transcriptional profiles of macrophages infected with both morphotypes were highly correlated, however, and a direct comparison identified few genes to distinguish them. Most of the induced miRNAs have previously been associated with mycobacterial infection and overall miRNA expression patterns were similarly highly correlated between the morphotypes.

Conclusions

The report here details the first whole transcriptome analysis of the early macrophage response to MAB infection. The overall picture at the early stages of macrophage infection is similar to that of other mycobacteria, reflected in a core type I IFN and pro-inflammatory cytokine response. Large-scale transcriptional differences in the host response to the different MAB morphotypes are not evident in the early stages of infection, however the subset of genes with distinct expression profiles suggest potentially interesting differences in internal trafficking of MAB within macrophages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2246-1) contains supplementary material, which is available to authorized users.

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.

RNA sequencing provides exquisite insight into the manipulation of the alveolar macrophage by tubercle bacilli

Mycobacterium bovis, the agent of bovine tuberculosis, causes an estimated $3 billion annual losses to global agriculture due, in part, to the limitations of current diagnostics. Development of next-generation diagnostics requires a greater understanding of the interaction between the pathogen and the bovine host. Therefore, to explore the early response of the alveolar macrophage to infection, we report the first application of RNA-sequencing to define, in exquisite detail, the transcriptomes of M. bovis-infected and non-infected alveolar macrophages from ten calves at 2, 6, 24 and 48 hours post-infection. Differentially expressed sense genes were detected at these time points that revealed enrichment of innate immune signalling functions, and transcriptional suppression of host defence mechanisms (e.g., lysosome maturation). We also detected differentially expressed natural antisense transcripts, which may play a role in subverting innate immune mechanisms following infection. Furthermore, we report differential expression of novel bovine genes, some of which have immune-related functions based on orthology with human proteins. This is the first in-depth transcriptomics investigation of the alveolar macrophage response to the early stages of M. bovis infection and reveals complex patterns of gene expression and regulation that underlie the immunomodulatory mechanisms used by M. bovis to evade host defence mechanisms.

MicroRNA profiling of the bovine alveolar macrophage response to Mycobacterium bovis infection suggests pathogen survival is enhanced by microRNA regulation of endocytosis and lysosome trafficking

Mycobacterium bovis, the causative agent of bovine tuberculosis, a major problem for global agriculture, spreads via an airborne route and is taken up by alveolar macrophages (AM) in the lung. Here, we describe the first next-generation sequencing (RNA-seq) approach to temporally profile miRNA expression in primary bovine AMs post-infection with M. bovis. One, six, and forty miRNAs were identified as significantly differentially expressed at 2, 24 and 48 h post-infection, respectively. The differential expression of three miRNAs (bta-miR-142-5p, bta-miR-146a, and bta-miR-423-3p) was confirmed by RT-qPCR. Pathway analysis of the predicted mRNA targets of differentially expressed miRNAs suggests that these miRNAs preferentially target several pathways that are functionally relevant for mycobacterial pathogenesis, including endocytosis and lysosome trafficking, IL-1 signalling and the TGF-β pathway. Over-expression studies using a bovine macrophage cell-line (Bomac) reveal the targeting of two key genes in the innate immune response to M. bovis, IL-1 receptor-associated kinase 1 (IRAK1) and TGF-β receptor 2 (TGFBR2), by miR-146. Taken together, our study suggests that miRNAs play a key role in tuning the complex interplay between M. bovis survival strategies and the host immune response.

Host gene expression for Mycobacterium avium subsp. paratuberculosis infection in human THP-1 macrophages

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease, which causes considerable economic loss in the dairy industry and has a possible relationship to Crohn's disease (CD) in humans. As MAP has been detected in retail pasteurized milk samples, its transmission via milk is of concern. Despite its possible role in the etiology of CD, there have been few studies examining the interactions between MAP and human cells. In the current study, we applied Ingenuity Pathway Analysis to the transcription profiles generated from a murine model with MAP infection as part of a previously conducted study. Twenty-one genes were selected as potential host immune responses, compared with the transcriptional profiles in naturally MAP-infected cattle, and validated in MAP-infected human monocyte-derived macrophage THP-1 cells. Of these, the potential host responses included up-regulation of genes related to immune response (CD14, S100A8, S100A9, LTF, HP and CHCIL3), up-regulation of Th1-polarizing factor (CCL4, CCL5, CXCL9 and CXCL10), down-regulation of genes related to metabolism (ELANE, IGF1, TCF7L2 and MPO) and no significant response of other genes (GADD45a, GPNMB, HMOX1, IFNG and NQO1) in THP-1 cells infected with MAP.

Analytical workflow profiling gene expression in murine macrophages

Comprehensive and simultaneous analysis of all genes in a biological sample is a capability of RNA-Seq technology. Analysis of the entire transcriptome benefits from summarization of genes at the functional level. As a cellular response of interest not previously explored with RNA-Seq, peritoneal macrophages from mice under two conditions (control and immunologically challenged) were analyzed for gene expression differences. Quantification of individual transcripts modeled RNA-Seq read distribution and uncertainty (using a Beta Negative Binomial distribution), then tested for differential transcript expression (False Discovery Rate-adjusted p-value < 0.05). Enrichment of functional categories utilized the list of differentially expressed genes. A total of 2079 differentially expressed transcripts representing 1884 genes were detected. Enrichment of 92 categories from Gene Ontology Biological Processes and Molecular Functions, and KEGG pathways were grouped into 6 clusters. Clusters included defense and inflammatory response (Enrichment Score = 11.24) and ribosomal activity (Enrichment Score = 17.89). Our work provides a context to the fine detail of individual gene expression differences in murine peritoneal macrophages during immunological challenge with high throughput RNA-Seq.

Analysis of the Bovine Monocyte-Derived Macrophage Response to Mycobacterium avium Subspecies Paratuberculosis Infection Using RNA-seq

Johne's disease, caused by infection with Mycobacterium avium subsp. paratuberculosis, (MAP), is a chronic intestinal disease of ruminants with serious economic consequences for cattle production in the United States and elsewhere. During infection, MAP bacilli are phagocytosed and subvert host macrophage processes, resulting in subclinical infections that can lead to immunopathology and dissemination of disease. Analysis of the host macrophage transcriptome during infection can therefore shed light on the molecular mechanisms and host-pathogen interplay associated with Johne's disease. Here, we describe results of an in vitro study of the bovine monocyte-derived macrophage (MDM) transcriptome response during MAP infection using RNA-seq. MDM were obtained from seven age- and sex-matched Holstein-Friesian cattle and were infected with MAP across a 6-h infection time course with non-infected controls. We observed 245 and 574 differentially expressed (DE) genes in MAP-infected versus non-infected control samples (adjusted P value ≤0.05) at 2 and 6 h post-infection, respectively. Functional analyses of these DE genes, including biological pathway enrichment, highlighted potential functional roles for genes that have not been previously described in the host response to infection with MAP bacilli. In addition, differential expression of pro- and anti-inflammatory cytokine genes, such as those associated with the IL-10 signaling pathway, and other immune-related genes that encode proteins involved in the bovine macrophage response to MAP infection emphasize the balance between protective host immunity and bacilli survival and proliferation. Systematic comparisons of RNA-seq gene expression results with Affymetrix(®) microarray data generated from the same experimental samples also demonstrated that RNA-seq represents a superior technology for studying host transcriptional responses to intracellular infection.

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.

Comparative functional genomics and the bovine macrophage response to strains of the mycobacterium genus

Mycobacterial infections are major causes of morbidity and mortality in cattle and are also potential zoonotic agents with implications for human health. Despite the implementation of comprehensive animal surveillance programs, many mycobacterial diseases have remained recalcitrant to eradication in several industrialized countries. Two major mycobacterial pathogens of cattle are Mycobacterium bovis and Mycobacterium avium subspecies paratuberculosis (MAP), the causative agents of bovine tuberculosis (BTB) and Johne's disease (JD), respectively. BTB is a chronic, granulomatous disease of the respiratory tract that is spread via aerosol transmission, while JD is a chronic granulomatous disease of the intestines that is transmitted via the fecal-oral route. Although these diseases exhibit differential tissue tropism and distinct complex etiologies, both M. bovis and MAP infect, reside, and replicate in host macrophages - the key host innate immune cell that encounters mycobacterial pathogens after initial exposure and mediates the subsequent immune response. The persistence of M. bovis and MAP in macrophages relies on a diverse series of immunomodulatory mechanisms, including the inhibition of phagosome maturation and apoptosis, generation of cytokine-induced necrosis enabling dissemination of infection through the host, local pathology, and ultimately shedding of the pathogen. Here, we review the bovine macrophage response to infection with M. bovis and MAP. In particular, we describe how recent advances in functional genomics are shedding light on the host macrophage-pathogen interactions that underlie different mycobacterial diseases. To illustrate this, we present new analyses of previously published bovine macrophage transcriptomics data following in vitro infection with virulent M. bovis, the attenuated vaccine strain M. bovis BCG, and MAP, and discuss our findings with respect to the differing etiologies of BTB and JD.

Key Hub and Bottleneck Genes Differentiate the Macrophage Response to Virulent and Attenuated Mycobacterium bovis

Mycobacterium bovis is an intracellular pathogen that causes tuberculosis in cattle. Following infection, the pathogen resides and persists inside host macrophages by subverting host immune responses via a diverse range of mechanisms. Here, a high-density bovine microarray platform was used to examine the bovine monocyte-derived macrophage transcriptome response to M. bovis infection relative to infection with the attenuated vaccine strain, M. bovis Bacille Calmette-Guérin. Differentially expressed genes were identified (adjusted P-value ≤0.01) and interaction networks generated across an infection time course of 2, 6, and 24 h. The largest number of biological interactions was observed in the 24-h network, which exhibited scale-free network properties. The 24-h network featured a small number of key hub and bottleneck gene nodes, including IKBKE, MYC, NFKB1, and EGR1 that differentiated the macrophage response to virulent and attenuated M. bovis strains, possibly via the modulation of host cell death mechanisms. These hub and bottleneck genes represent possible targets for immuno-modulation of host macrophages by virulent mycobacterial species that enable their survival within a hostile environment.

RNA-seq Transcriptional Profiling of Peripheral Blood Leukocytes from Cattle Infected with Mycobacterium bovis

Bovine tuberculosis, caused by infection with 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 gene expression microarrays and RNA-sequencing (RNA-seq), 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 analyzed the peripheral blood leukocyte (PBL) transcriptome of eight natural M. bovis-infected and eight age- and sex-matched non-infected control Holstein-Friesian animals using RNA-seq. In addition, we compared gene expression profiles generated using RNA-seq with those previously generated using the high-density Affymetrix^® GeneChip^® Bovine Genome Array platform from the same PBL-extracted RNA. A total of 3,250 differentially expressed (DE) annotated genes were detected in the M. bovis-infected samples relative to the controls (adjusted P-value ≤0.05), with the number of genes displaying decreased relative expression (1,671) exceeding those with increased relative expression (1,579). Ingenuity^® Systems Pathway Analysis (IPA) of all DE genes revealed enrichment for genes with immune function. Notably, transcriptional suppression was observed among several of the top-ranking canonical pathways including Leukocyte Extravasation Signaling. Comparative platform analysis demonstrated that RNA-seq detected a larger number of annotated DE genes (3,250) relative to the microarray (1,398), of which 917 genes were common to both technologies and displayed the same direction of expression. Finally, we show that RNA-seq had an increased dynamic range compared to the microarray for estimating differential gene expression.

The quality of methods reporting in parasitology experiments

There is a growing concern both inside and outside the scientific community over the lack of reproducibility of experiments. The depth and detail of reported methods are critical to the reproducibility of findings, but also for making it possible to compare and integrate data from different studies. In this study, we evaluated in detail the methods reporting in a comprehensive set of trypanosomiasis experiments that should enable valid reproduction, integration and comparison of research findings. We evaluated a subset of other parasitic (Leishmania, Toxoplasma, Plasmodium, Trichuris and Schistosoma) and non-parasitic (Mycobacterium) experimental infections in order to compare the quality of method reporting more generally. A systematic review using PubMed (2000-2012) of all publications describing gene expression in cells and animals infected with Trypanosoma spp was undertaken based on PRISMA guidelines; 23 papers were identified and included. We defined a checklist of essential parameters that should be reported and have scored the number of those parameters that are reported for each publication. Bibliometric parameters (impact factor, citations and h-index) were used to look for association between Journal and Author status and the quality of method reporting. Trichuriasis experiments achieved the highest scores and included the only paper to score 100% in all criteria. The mean of scores achieved by Trypanosoma articles through the checklist was 65.5% (range 32-90%). Bibliometric parameters were not correlated with the quality of method reporting (Spearman's rank correlation coefficient <-0.5; p>0.05). Our results indicate that the quality of methods reporting in experimental parasitology is a cause for concern and it has not improved over time, despite there being evidence that most of the assessed parameters do influence the results. We propose that our set of parameters be used as guidelines to improve the quality of the reporting of experimental infection models as a pre-requisite for integrating and comparing sets of data.

Phenotypic and functional heterogeneity of bovine blood monocytes

Murine and human peripheral blood monocytes are heterogeneous in size, granularity, nuclear morphology, phenotype and function. Whether and how bovine blood monocytes follow this pattern was analyzed in this study. Flow cytometrically, classical monocytes (cM) CD14⁺ CD16⁻, intermediate monocytes (intM) CD14⁺ CD16⁺ and nonclassical monocytes (ncM) CD14⁺ CD16⁺ were identified, with cM being the predominant subset (89%). cM showed a significant lower expression of CD172a, intM expressed the highest level of MHC class II molecules, and ncM were low positive for CD163. Compared to cM and intM, ncM showed a significantly reduced phagocytosis capacity, a significantly reduced generation of reactive oxygen species, and reduced mRNA expression of CXCL8, CXCL1 and IL-1β after LPS stimulation. Based on IL-1β secretion after LPS/ATP stimulation, the inflammasome could be activated in cM and intM, but not in ncM. IFNγ increased the expression of CD16 selectively on cM and induced a shift from cM into intM in vitro. In summary, bovine CD172a-positive mononuclear cells define three monocyte subsets with distinct phenotypic and functional differences. Bovine cM and intM share homologies with their human counterparts, whereas bovine ncM are not inflammatory monocytes.

Receptor-mediated recognition of mycobacterial pathogens

Mycobacteria are a genus of bacteria that range from the non-pathogenic Mycobacterium smegmatis to Mycobacterium tuberculosis, the causative agent of tuberculosis in humans. Mycobacteria primarily infect host tissues through inhalation or ingestion. They are phagocytosed by host macrophages and dendritic cells. Here, conserved pathogen-associated molecular patterns (PAMPs) on the surface of mycobacteria are recognized by phagocytic pattern recognition receptors (PRRs). Several families of PRRs have been shown to non-opsonically recognize mycobacterial PAMPs, including membrane-bound C-type lectin receptors, membrane-bound and cytosolic Toll-like receptors and cytosolic NOD-like receptors. Recently, a possible role for intracellular cytosolic PRRs in the recognition of mycobacterial pathogens has been proposed. Here, we discuss currentideas on receptor-mediated recognition of mycobacterial pathogens by macrophages and dendritic cells.

Transcriptome Analysis Reveals Novel Entry Mechanisms and a Central Role of SRC in Host Defense during High Multiplicity Mycobacterial Infection

Mycobacterium tuberculosis (MTB) infects an estimated one-third of the global population and is one of the main causes of mortality from an infectious agent. The characteristics of macrophages challenged by MTB with a high multiplicity of infection (MOI), which mimics both clinical disseminated infection and granuloma formation, are distinct from macrophages challenged with a low MOI. To better understand the cross talk between macrophage host cells and mycobacteria, we compared the transcription patterns of mouse macrophages infected with bacille Calmette-Guérin, H37Ra and M. smegmatis. Attention was focused on the changes in the abundance of transcripts related to immune system function. From the results of a transcriptome profiling study with a high mycobacterial MOI, we defined a pathogen-specific host gene expression pattern. The present study suggests that two integrins, ITGA5 and ITGAV, are novel cell surface receptors mediating mycobacterium entry into macrophages challenged with high MOI. Our results indicate that SRC likely plays a central role in regulating multiple unique signaling pathways activated by MTB infection. The integrated results increase our understanding of the molecular networks behind the host innate immune response and identify important targets that might be useful for the development of tuberculosis therapy.

Global gene expression profiling of monocyte-derived macrophages from red deer (Cervus elaphus) genotypically resistant or susceptible to Mycobacterium avium subspecies paratuberculosis infection

Mycobacterium avium subspecies paratuberculosis (MAP) can cause a chronic inflammatory bowel disease, Johne's disease (JD), in ruminant animals. This study has explored the molecular basis of resistance and susceptibility to this disease in red deer breeds previously confirmed to express polarised phenotypes by experimental infection trials and following natural infection. Monocyte-derived macrophage cultures were obtained from uninfected red deer selected for either a resistant or susceptible phenotype. Cells were infected with MAP in vitro and gene expression analysed by RNA-Seq. Transcriptome analysis revealed a more disrupted gene expression profile in macrophages from susceptible animals compared with cells from resistant animals in terms of the number of genes up- or downregulated. Highly upregulated genes were related to chemotaxis (CXCL10, CSF3, and CCL8) and type 1 interferon signalling (RSAD2, IFIT1, IFIT2, ISG12, ISG15, USP18, and HERC6). Upregulation of these genes was observed to be greater in macrophages from susceptible animals compared to cells from resistant animals in response to in vitro MAP infection. These data support the use of transcriptomic approaches to enable the identification of markers associated particularly with susceptibility to MAP infection.

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.

Transcriptional response of bovine monocyte-derived macrophages after the infection with different Argentinean Mycobacterium bovis isolates

Infection of bovines with Mycobacterium bovis causes important financial hardship in many countries presenting also a risk for humans. M. bovis is known to be adapted to survive and thrive within the intramacrophage environment. In spite of its relevance, at present the information about macrophage expression patterns is scarce, particularly regarding the bovine host. In this study, transcriptomic analysis was used to detect genes differentially expressed in macrophages derived from peripheral blood mononuclear cells at early stages of infection with two Argentinean strains of M. bovis, a virulent and an attenuated strains. The results showed that the number of differentially expressed genes in the cells infected with the virulent strain (5) was significantly lower than those in the cells infected with the attenuated strain (172). Several genes were more strongly expressed in infected macrophages. Among them, we detected encoding transcription factors, anthrax toxin receptor, cell division and apoptosis regulator, ankyrin proteins, cytoskeleton proteins, protein of cell differentiation, and regulators of endocytic traffic of membrane. Quantitative real-time PCR of a selected group of differentially expressed genes confirmed the microarrays results. Altogether, the present results contribute to understanding the mechanisms involved in the early interaction of M. bovis with the bovine macrophage.

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

Mycobacterium bovis is the causative agent of most cases of bovine tuberculosis. The identification of bTB biomarkers in specific stages of the disease will contribute to a better understanding of the immunopathology associated with tuberculosis and will enable their use in disease diagnosis and prognosis. The aim of this study was to evaluate the gene expression profile induced after specific stimulation of bovine peripheral blood mononuclear cells from cattle infected with M. bovis using the Affymetrix® GeneChip® Bovine Genome Array. A total of 172 genes showed differential expression profile that was statistically significant with log2-fold change >2.5 and <−2.5. Twenty-four out of these genes were upregulated and 148 were downregulated in bovine peripheral blood mononuclear cells of M. bovis-infected cattle. The highest differentially-expressed genes were related to immune and inflammatory responses, apoptosis, endocytosis, cellular trafficking and genes encoding proteins involved in cellular matrix degradation. Microarray results were confirmed in another group of infected cattle by RT-qPCR for the CD14, IL-1R, THBS1, MMP9 and FYVE genes. This study confirms previous findings that have shown that M. bovis infection in cattle results in the downregulation of immune response-related genes. Moreover, it validates the use of microarray platforms in combination with RT-qPCR to identify biomarkers of bovine tuberculosis. In addition, we propose CD14, IL-1R, THBS1, MMP9 and FYVE as potential biomarkers of bovine tuberculosis.