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

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.

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.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages

Background

Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne’s disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (−); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip.

Results

The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/−) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10^− 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids.

Conclusion

In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07487-4.

Weighted Gene Co-Expression Network Analysis Identifies Key Modules and Hub Genes Associated with Mycobacterial Infection of Human Macrophages

Tuberculosis (TB) is still a leading cause of death worldwide. Treatments remain unsatisfactory due to an incomplete understanding of the underlying host–pathogen interactions during infection. In the present study, weighted gene co-expression network analysis (WGCNA) was conducted to identify key macrophage modules and hub genes associated with mycobacterial infection. WGCNA was performed combining our own transcriptomic results using Mycobacterium aurum-infected human monocytic macrophages (THP1) with publicly accessible datasets obtained from three types of macrophages infected with seven different mycobacterial strains in various one-to-one combinations. A hierarchical clustering tree of 11,533 genes was built from 198 samples, and 47 distinct modules were revealed. We identified a module, consisting of 226 genes, which represented the common response of host macrophages to different mycobacterial infections that showed significant enrichment in innate immune stimulation, bacterial pattern recognition, and leukocyte chemotaxis. Moreover, by network analysis applied to the 74 genes with the best correlation with mycobacteria infection, we identified the top 10 hub-connecting genes: NAMPT, IRAK2, SOCS3, PTGS2, CCL20, IL1B, ZC3H12A, ABTB2, GFPT2, and ELOVL7. Interestingly, apart from the well-known Toll-like receptor and inflammation-associated genes, other genes may serve as novel TB diagnosis markers and potential therapeutic targets.

Transcription Profiling of Monocyte-Derived Macrophages Infected In Vitro With Two Strains of Streptococcus agalactiae Reveals Candidate Pathways Affecting Subclinical Mastitis in Cattle

Macrophages are key cells of innate immune response and serve as the first line of defense against bacteria. Transcription profiling of bacteria-infected macrophages could provide important insights on the pathogenicity and host defense mechanisms during infection. We have examined transcription profiles of bovine monocyte-derived macrophages (bMDMs) isolated from the blood of 12 animals and infected in vitro with two strains of Streptococcus agalactiae. Illumina sequencing of RNA from 36 bMDMs cultures exposed in vitro to either one of two sequence types of S. agalactiae (ST103 or ST12) for 6 h and unchallenged controls was performed. Analyses of over 1,656 million high-quality paired-end sequence reads revealed 5,936 and 6,443 differentially expressed genes (p < 0.05) in bMDMs infected with ST103 and ST12, respectively, versus unchallenged controls. Moreover, 588 genes differentially expressed between bMDMs infected with ST103 versus ST12 were identified. Ingenuity pathway analysis of the differentially up-regulated genes in the bMDMs infected with ST103 revealed significant enrichment for granulocyte adhesion and diapedesis, while significant enrichment for the phagosome formation pathway was found among down-regulated genes. Moreover, Ingenuity pathway analysis of the differentially up-regulated genes in the bMDMs infected with ST12 showed significant enrichment for type 1/type 2 T helper cell activation, while the complement activation pathway was overrepresented in the down-regulated genes. Our study identified pathogen-induced regulation of key genes and pathways involved in the immune response of macrophages against infection but also likely involved in bacterial evasion of the host immune system. These results may contribute to better understanding of the mechanisms underlying subclinical infection such as bovine streptococcal mastitis.

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.

Establishment and characterization of a telomerase-immortalized porcine bronchial epithelial cell line

Primary porcine bronchial epithelial cells (PBECs) are an ideal model to study the molecular and pathogenic mechanisms of various porcine respiratory pathogens. However, the short lifespan of primary PBECs greatly limit their application. Here, we isolated and cultured primary PBECs and established immortalized PBECs by transfecting primary PBECs with the pEGFP-hTERT recombinant plasmid containing human telomerase reverse transcriptase (hTERT). Immortalized PBECs (hTERT-PBECs) retained the morphological and functional features of primary PBECs as indicated by cytokeratin 18 expression, telomerase activity assay, proliferation assays, karyotype analysis, and quantitative reverse-transcriptase polymerase chain reaction. Compared to primary PBECs, hTERT-PBECs had higher telomerase activity, extended replicative lifespan, and displayed enhanced proliferative activity. Moreover, this cell line is not transformed in vitro and does not exhibit a malignant phenotype in vivo, suggesting that it can be safely used in further studies. Besides, hTERT-PBECs were susceptible to swine influenza virus of H3N2 subtype and porcine circovirus type 2. In conclusion, the immortalized hTERT-PBECs represent a valuable in vitro model, which can be widely used in the study of porcine respiratory pathogenic infections.

A Markov chain representation of the multiple testing problem

The problem of multiple hypothesis testing can be represented as a Markov process where a new alternative hypothesis is accepted in accordance with its relative evidence to the currently accepted one. This virtual and not formally observed process provides the most probable set of non null hypotheses given the data; it plays the same role as Markov Chain Monte Carlo in approximating a posterior distribution. To apply this representation and obtain the posterior probabilities over all alternative hypotheses, it is enough to have, for each test, barely defined Bayes Factors, e.g. Bayes Factors obtained up to an unknown constant. Such Bayes Factors may either arise from using default and improper priors or from calibrating p-values with respect to their corresponding Bayes Factor lower bound. Both sources of evidence are used to form a Markov transition kernel on the space of hypotheses. The approach leads to easy interpretable results and involves very simple formulas suitable to analyze large datasets as those arising from gene expression data (microarray or RNA-seq experiments).

RNA-Seq revealed ATF3-regulated inflammation induced by silica

BACKGROUND

Millions of workers are exposed to dust containing silica. Chronic and over-exposure to silica will lead to silicosis, which is an irreversible and sometimes fatal lung disease. The disordered physiological processes of silicosis consist of accumulation of silica particles in the alveoli of the lung. Then, the ingestion of the silica particles by macrophages was followed by an inflammatory response. Up till now, the chest radiographs remain the key tool in diagnosing and assessing the extent of silicosis. However, concerns exist regarding the sensitivity and specificity of the technique. Therefore, there is still a need to develop a biomarker for silicosis for early detection of silicosis.

METHOD

In this study, RNA-Seq was applied to detect the gene expression changes when silica was exposed to macrophages at different time intervals. RNA-Seq provides a broader dynamic range, increased specificity and sensitivity, and easier detection of rare and low-abundance transcripts. Bioinformatics tools such as the Database for Annotation, Visualization and Integrated Discovery (DAVID) and Gene Functional Classification Tool and Search Tool for the Retrieval of Interacting Genes (STRING) were applied for data analysis. Quantitative PCR was used to validate the results.

RESULTS

Our results showed that regulation of transcription factors was the dominant activated pathway in early exposure of silica to macrophages, followed by inflammatory responses which were the main mechanisms in silicosis. One of the findings was the upregulation of activating transcription factor 3 (ATF3) during silica exposure. When ATF3 expression was inhibited by siRNA, the production of cytokines IL-1β, IL-6 and TNF was further increased.

CONCLUSION

This indicated that ATF3 may be a potential early diagnostic biomarker for silicosis and ATF3 acts as a repressor in inflammatory responses induced by silica.

P-value calibration in multiple hypotheses testing

As p-values are the most common measures of evidence against a hypothesis, their calibration with respect to null hypothesis conditional probability is important in order to match frequentist unconditional inference with the Bayesian ones. The Selke, Bayarri and Berger calibration is one of the most popular attempts to obtain such a calibration. This relies on the theoretical sampling null distribution of p-values, which is the well-known Uniform(0,1), but arising only for specific sampling models. We generalize this calibration by considering a sampling null distribution estimated from the data. It is possible to obtain such an empirical null distribution, for instance, in the context of multiple testing in which many p-values come from the null model. Such a context is purely instrumental for the purposes of p-value calibration, and multiple testing still needs to be considered with appropriate techniques. The new calibration proposed here still remains a simple analytic formula like the original one under the Uniform(0,1) and basically provides a stronger interpretation framework for the widely used p-value. Copyright © 2017 John Wiley & Sons, Ltd.

Profiling of the silica-induced molecular events in lung epithelial cells using the RNA-Seq approach

Silicosis is a prolonged, irreversible and incurable occupational disease, and there is a significant number of newly diagnosed cases every year in Hong Kong. Due to the long latency of the disease, the diagnosis can be missed until detailed clinical examination at a later stage. For a better control of this deadly disease, detailing the pro-inflammatory and fibrotic events in the macrophage would be instrumental in understanding the pathogenesis of the disease and essential for the significant biomarkers discovery. In this in vitro study, human cell line model A549 lung epithelial cells were used. The immediate molecular events underneath the activation of quartz silica polymorphs were followed in a time course of 0, 0.5, 2, 8, 16 and 24 h. The transcriptome library was prepared and subjected to RNA-Seq analysis. Data analysis was performed by pathway analysis tools and verified by real-time PCR. The results showed that triggered genes were mainly found in the immune response and inflammatory pathways. An interesting finding was the association of the DNA-binding protein inhibitor (ID) family in the silica exposure to lung cells. The linkage of ID1, ID2 and ID3 to cancer may rationalize themselves to be the markers indicating an early response of silicosis. However, further studies are required to consolidate the roles of these genes in silicosis. Copyright © 2017 John Wiley & Sons, Ltd.

An Optimized Method for Extracting Bacterial RNA from Mouse Skin Tissue Colonized by Mycobacterium ulcerans

Bacterial transcriptome analyses during host colonization are essential to decipher the complexity of the relationship between the bacterium and its host. RNA sequencing (RNA-seq) is a promising approach providing valuable information about bacterial adaptation, the host response and, in some cases, mutual tolerance underlying crosstalk, as recently observed in the context of Mycobacterium ulcerans infection. Buruli ulcer is caused by M. ulcerans. This neglected disease is the third most common mycobacterial disease worldwide. Without treatment, M. ulcerans provokes massive skin ulcers. A healing process may be observed in 5% of Buruli ulcer patients several months after the initiation of disease. This spontaneous healing process suggests that some hosts can counteract the development of the lesions caused by M. ulcerans. Deciphering the mechanisms involved in this process should open up new treatment possibilities. To this end, we recently developed the first mouse model for studies of the spontaneous healing process. We have shown that the healing process is based on mutual tolerance between the bacterium and its host. In this context, RNA-seq seems to be the most appropriate method for deciphering bacterial adaptation. However, due to the low bacterial load in host tissues, the isolation of mycobacterial RNA from skin tissue for RNA-seq analysis remains challenging. We developed a method for extracting and purifying mycobacterial RNA whilst minimizing the amount of host RNA in the sample. This approach was based on the extraction of bacterial RNA by a differential lysis method. The challenge in the development of this method was the choice of a lysis system favoring the removal of host RNA without damage to the bacterial cells. We made use of the thick, resistant cell wall of M. ulcerans to achieve this end.

Circulating microRNAs as Potential Biomarkers of Infectious Disease

microRNAs (miRNAs) are a class of small non-coding endogenous RNA molecules that regulate a wide range of biological processes by post-transcriptionally regulating gene expression. Thousands of these molecules have been discovered to date, and multiple miRNAs have been shown to coordinately fine-tune cellular processes key to organismal development, homeostasis, neurobiology, immunobiology, and control of infection. The fundamental regulatory role of miRNAs in a variety of biological processes suggests that differential expression of these transcripts may be exploited as a novel source of molecular biomarkers for many different disease pathologies or abnormalities. This has been emphasized by the recent discovery of remarkably stable miRNAs in mammalian biofluids, which may originate from intracellular processes elsewhere in the body. The potential of circulating miRNAs as biomarkers of disease has mainly been demonstrated for various types of cancer. More recently, however, attention has focused on the use of circulating miRNAs as diagnostic/prognostic biomarkers of infectious disease; for example, human tuberculosis caused by infection with Mycobacterium tuberculosis, sepsis caused by multiple infectious agents, and viral hepatitis. Here, we review these developments and discuss prospects and challenges for translating circulating miRNA into novel diagnostics for infectious disease.

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.

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.

Transcriptional Profiling of Ileocecal Valve of Holstein Dairy Cows Infected with Mycobacterium avium subsp. Paratuberculosis

Johne’s disease is a chronic infection of the small intestine caused by Mycobacterium avium subspecies paratuberculosis (MAP), an intracellular bacterium. The events of pathogen survival within the host cell(s), chronic inflammation and the progression from asymptomatic subclinical stage to an advanced clinical stage of infection, are poorly understood. This study examines gene expression in the ileocecal valve (ICV) of Holstein dairy cows at different stages of MAP infection. The ICV is known to be a primary site of MAP colonization and provides an ideal location to identify genes that are relevant to the progression of this disease. RNA was prepared from ICV tissues and RNA-Seq was used to compare gene transcription between clinical, subclinical, and uninfected control animals. Interpretation of the gene expression data was performed using pathway analysis and gene ontology categories containing multiple differentially expressed genes. Results demonstrated that many of the pathways that had strong differential gene expression between uninfected control and clinical cows were related to the immune system, such as the T- and B-cell receptor signaling, apoptosis, NOD-like receptor signaling, and leukocyte transendothelial migration pathways. In contrast, the comparison of gene transcription between control and subclinical cows identified pathways that were primarily involved in metabolism. The results from the comparison between clinical and subclinical animals indicate recruitment of neutrophils, up regulation of lysosomal peptidases, increase in immune cell transendothelial migration, and modifications of the extracelluar matrix. This study provides important insight into how cattle respond to a natural MAP infection at the gene transcription level within a key target tissue for infection.

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.

The Transcriptional Foundations of Sp110-mediated Macrophage (RAW264.7) Resistance to Mycobacterium tuberculosis H37Ra

Human tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a leading global health problem, causing 1.3 million deaths each year. The nuclear body protein, Sp110, has been linked to TB resistance and previous work showed that it enhances macrophage apoptosis upon Mtb infection. Here, we report on the role of Sp110 in transcriptional regulation of macrophage responses to Mtb through integrated transcriptome and mechanistic studies. Transcriptome analysis revealed that Sp110 regulates genes involved in immune responses, apoptosis, defence responses, and inflammatory responses. Detailed investigation revealed that, in addition to apoptosis-related genes, Sp110 regulates cytokines, chemokines and genes that regulate intracellular survival of Mtb. Moreover, Sp110 regulates miRNA expression in macrophages, with immune and apoptosis-related miRNAs such as miR-125a, miR-146a, miR-155, miR-21a and miR-99b under Sp110 regulation. Additionally, our results showed that Sp110 upregulates BCL2 modifying factor (Bmf) by inhibiting miR-125a, and forced expression of Bmf induces macrophage apoptosis. These findings not only reveal the transcriptional basis of Sp110-mediated macrophage resistance to Mtb, but also suggest potential regulatory roles for Sp110 related to inflammatory responses, miRNA profiles, and the intracellular growth of Mtb.

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.

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.

Different Blood Cell-Derived Transcriptome Signatures in Cows Exposed to Vaccination Pre- or Postpartum

Periparturient cows have been found to reveal immunosuppression, frequently associated with increased susceptibility to uterine and mammary infections. To improve understanding of the causes and molecular regulatory mechanisms accounting for this phenomenon around calving, we examined the effect of an antigen challenge on gene expression modulation on cows prior to (BC) or after calving (AC) using whole transcriptome sequencing (RNAseq). The transcriptome analysis of the cows’ blood identified a substantially higher number of loci affected in BC cows (2,235) in response to vaccination compared to AC cows (208) and revealed a divergent transcriptional profile specific for each group. In BC cows, a variety of loci involved in immune defense and cellular signaling processes were transcriptionally activated, whereas protein biosynthesis and posttranslational processes were tremendously impaired in response to vaccination. Furthermore, energy metabolism in the blood cells of BC cows was shifted from oxidative phosphorylation to the glycolytic system. In AC cows, the number and variety of regulated pathways involved in immunomodulation and maintenance of immnunocompetence are considerably lower after vaccination, and upregulation of arginine degradation was suggested as an immunosuppressive mechanism. Elevated transcript levels of erythrocyte-specific genes involved in gas exchange processes were a specific transcriptional signature in AC cows pointing to hematopoiesis activation. The divergent and substantially lower magnitude of transcriptional modulation in response to vaccination in AC cows provides evidence for a suppressed immune capacity of early lactating cows on the molecular level and demonstrates that an efficient immune response of cows is related to their physiological and metabolic status.

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.

Gene expression network analysis reveals new transcriptional regulators as novel factors in human ischemic cardiomyopathy

BACKGROUND

Ischemic cardiomyopathy (ICM) is characterized by transcriptomic changes that alter cellular processes leading to decreased cardiac output. Because the molecular network of ICM is largely unknown, the aim of this study was to characterize the role of new transcriptional regulators in the molecular mechanisms underlying the responses to ischemia.

METHODS

Myocardial tissue explants from ICM patients and control (CNT) subjects were analyzed by RNA-Sequencing (RNA-Seq) and quantitative Real-Time PCR.

RESULTS

Enrichment analysis of the ICM transcriptomic profile allowed the characterization of novel master regulators. We found that the expression of the transcriptional regulators SP100 (-1.5-fold, p < 0.05), CITED2 (-3.8-fold, p < 0.05), CEBPD (-4.9-fold, p < 0.05) and BCL3 (-3.3-fold, p < 0.05) were lower in ICM than in CNT. To gain insights into the molecular network defined by the transcription factors, we identified CEBPD, BCL3, and HIF1A target genes in the RNA-Seq datasets. We further characterized the biological processes of the target genes by gene ontology annotation. Our results suggest that CEBPD-inducible genes with roles in the inhibition of apoptosis are downregulated and that BCL3-repressible genes are involved in the regulation of cellular metabolism in ICM. Moreover, our results suggest that CITED2 downregulation causes increased expression of HIF1A target genes. Functional analysis of HIF1A target genes revealed that hypoxic and stress response genes are activated in ICM. Finally, we found a significant correlation between the mRNA levels of BCL3 and the mRNA levels of both CEBPD (r = 0.73, p < 0.001) and CITED2 (r = 0.56, p < 0.05). Interestingly, CITED2 mRNA levels are directly related to ejection fraction (EF) (r = 0.54, p < 0.05).

CONCLUSIONS

Our data indicate that changes in the expression of SP100, CITED2, CEBPD, and BCL3 affect their transcription regulatory networks, which subsequently alter a number of biological processes in ICM patients. The relationship between CITED2 mRNA levels and EF emphasizes the importance of this transcription factor in ICM. Moreover, our findings identify new mechanisms used to interpret gene expression changes in ICM and provide valuable resources for further investigation of the molecular basis of human cardiac ischemic response.

Dual RNA sequencing reveals the expression of unique transcriptomic signatures in lipopolysaccharide-induced BV-2 microglial cells

Microglial cells become rapidly activated through interactions with pathogens, and the persistent activation of these cells is associated with various neurodegenerative diseases. Previous studies have investigated the transcriptomic signatures in microglia or macrophages using microarray technologies. However, this method has numerous restrictions, such as spatial biases, uneven probe properties, low sensitivity, and dependency on the probes spotted. To overcome this limitation and identify novel transcribed genes in response to LPS, we used RNA Sequencing (RNA-Seq) to determine the novel transcriptomic signatures in BV-2 microglial cells. Sequencing assessment and quality evaluation showed that approximately 263 and 319 genes (≥ 1.5 log₂-fold), such as cytokines and chemokines, were strongly induced after 2 and 4 h, respectively, and the induction of several genes with unknown immunological functions was also observed. Importantly, we observed that previously unidentified transcription factors (TFs) (irf1, irf7, and irf9), histone demethylases (kdm4a) and DNA methyltransferases (dnmt3l) were significantly and selectively expressed in BV-2 microglial cells. The gene expression levels, transcription start sites (TSS), isoforms, and differential promoter usage revealed a complex pattern of transcriptional and post-transcriptional gene regulation upon infection with LPS. In addition, gene ontology, molecular networks and pathway analyses identified the top significantly regulated functional classification, canonical pathways and network functions at each activation status. Moreover, we further analyzed differentially expressed genes to identify transcription factor (TF) motifs (−950 to +50 bp of the 5’ upstream promoters) and epigenetic mechanisms. Furthermore, we confirmed that the expressions of key inflammatory genes as well as pro-inflammatory mediators in the supernatants were significantly induced in LPS treated primary microglial cells. This transcriptomic analysis is the first to show a comparison of the family-wide differential expression of most known immune genes and also reveal transcription evidence of multiple gene families in BV-2 microglial cells. Collectively, these findings reveal unique transcriptomic signatures in BV-2 microglial cells required for homeostasis and effective immune responses.

Transcriptome characterization of three wild Chinese Vitis uncovers a large number of distinct disease related genes

BACKGROUND

Grape is one of the most valuable fruit crops and can serve for both fresh consumption and wine production. Grape cultivars have been selected and evolved to produce high-quality fruits during their domestication over thousands of years. However, current widely planted grape cultivars suffer extensive loss to many diseases while most wild species show resistance to various pathogens. Therefore, a comprehensive evaluation of wild grapes would contribute to the improvement of disease resistance in grape breeding programs.

RESULTS

We performed deep transcriptome sequencing of three Chinese wild grapes using the Illumina strand-specific RNA-Seq technology. High quality transcriptomes were assembled de novo and more than 93% transcripts were shared with the reference PN40024 genome. Over 1,600 distinct transcripts, which were absent or highly divergent from sequences in the reference PN40024 genome, were identified in each of the three wild grapes, among which more than 1,000 were potential protein-coding genes. Gene Ontology (GO) and pathway annotations of these distinct genes showed those involved in defense responses and plant secondary metabolisms were highly enriched. More than 87,000 single nucleotide polymorphisms (SNPs) and 2,000 small insertions or deletions (indels) were identified between each genotype and PN40024, and approximately 20% of the SNPs caused nonsynonymous mutations. Finally, we discovered 100 to 200 highly confident cis-natural antisense transcript (cis-NAT) pairs in each genotype. These transcripts were significantly enriched with genes involved in secondary metabolisms and plant responses to abiotic stresses.

CONCLUSION

The three de novo assembled transcriptomes provide a comprehensive sequence resource for molecular genetic research in grape. The newly discovered genes from wild Vitis, as well as SNPs and small indels we identified, may facilitate future studies on the molecular mechanisms related to valuable traits possessed by these wild Vitis and contribute to the grape breeding programs. Furthermore, we identified hundreds of cis-NAT pairs which showed their potential regulatory roles in secondary metabolism and abiotic stress responses.

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.

The Sweden Cancerome Analysis Network - Breast (SCAN-B) Initiative: a large-scale multicenter infrastructure towards implementation of breast cancer genomic analyses in the clinical routine

BACKGROUND

Breast cancer exhibits significant molecular, pathological, and clinical heterogeneity. Current clinicopathological evaluation is imperfect for predicting outcome, which results in overtreatment for many patients, and for others, leads to death from recurrent disease. Therefore, additional criteria are needed to better personalize care and maximize treatment effectiveness and survival.

METHODS

To address these challenges, the Sweden Cancerome Analysis Network - Breast (SCAN-B) consortium was initiated in 2010 as a multicenter prospective study with longsighted aims to analyze breast cancers with next-generation genomic technologies for translational research in a population-based manner and integrated with healthcare; decipher fundamental tumor biology from these analyses; utilize genomic data to develop and validate new clinically-actionable biomarker assays; and establish real-time clinical implementation of molecular diagnostic, prognostic, and predictive tests. In the first phase, we focus on molecular profiling by next-generation RNA-sequencing on the Illumina platform.

RESULTS

In the first 3 years from 30 August 2010 through 31 August 2013, we have consented and enrolled 3,979 patients with primary breast cancer at the seven hospital sites in South Sweden, representing approximately 85% of eligible patients in the catchment area. Preoperative blood samples have been collected for 3,942 (99%) patients and primary tumor specimens collected for 2,929 (74%) patients. Herein we describe the study infrastructure and protocols and present initial proof of concept results from prospective RNA sequencing including tumor molecular subtyping and detection of driver gene mutations. Prospective patient enrollment is ongoing.

CONCLUSIONS

We demonstrate that large-scale population-based collection and RNA-sequencing analysis of breast cancer is feasible. The SCAN-B Initiative should significantly reduce the time to discovery, validation, and clinical implementation of novel molecular diagnostic and predictive tests. We welcome the participation of additional comprehensive cancer treatment centers.

TRIAL REGISTRATION

ClinicalTrials.gov identifier NCT02306096.

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.