Alveolar Macrophage Chromatin Is Modified to Orchestrate Host Response to Mycobacterium bovis Infection

Extensive differential DNA methylation between tuberculosis skin test positive and skin test negative cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), represents a significant problem for the agriculture industry as well as posing a risk for human health. Current diagnostic tests for bTB target the cell-mediated immune (CMI) response to infection with M. bovis, primarily through screening of animals with the tuberculin skin test. Epigenetic modifications have been shown to alter the course of the immune response and differentially methylated regions (DMRs) might also influence the outcome of the skin test in cattle. Whole Genome Bisulphite Sequencing (WGBS) was used to profile DNA methylation levels from peripheral blood of a group of cattle identified as test positive for M. bovis (positive for the single intradermal comparative tuberculin test (SICTT) and/or the interferon-γ release assay compared to a test negative control group [n = 8/group, total of 16 WGBS libraries]. Although global methylation profiles were similar for both groups across the genome, 223 DMRs and 159 Differentially Promoter Methylated Genes (DPMGs) were identified between groups with an excess of hypermethylated sites in SICTT positive cattle (threshold > 15% differential methylation). Genes located within these DMRs included the Interleukin 1 receptor (IL1R1) and MHC related genes (BOLA and BOLA-DQB). KEGG pathway analysis identified enrichment of genes involved in Calcium and MAPK signalling, as well as metabolism pathways. Analysis of DMRs in a subset of SICTT negative cattle that were IFN-γ positive showed differential methylation of genes including Interleukin 10 Receptor, alpha (IL10RA), Interleukin 17 F (IL17F) and host defence peptides (DEFB and BDEF109). This study has identified a number of immune gene loci at which differential methylation is associated with SICTT test results and the degree of methylation could influence effective host immune responses.

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

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.

Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity

Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.

Surveying the Epigenetic Landscape of Tuberculosis in Alveolar Macrophages

Tuberculosis (TB) remains the leading cause of bacterial disease-related death and is among the top 10 overall causes of death worldwide. The complex nature of this infectious lung disease has proven difficult to treat, and significant research efforts are now evaluating the feasibility of host-directed, adjunctive therapies. An attractive approach in host-directed therapy targets host epigenetics, or gene regulation, to redirect the immune response in a host-beneficial manner. Substantial evidence exists demonstrating that host epigenetics are dysregulated during TB and that epigenetic-based therapies may be highly effective to treat TB. However, the caveat is that much of the knowledge that exists on the modulation of the host epigenome during TB has been gained using in vitro, small-animal, or blood-derived cell models, which do not accurately reflect the pulmonary nature of the disease. In humans, the first and major target cells of Mycobacterium tuberculosis are alveolar macrophages (AM). As such, their response to infection and treatment is clinically relevant and ultimately drives the outcome of disease. In this review, we compare the fundamental differences between AM and circulating monocyte-derived macrophages in the context of TB and summarize the recent advances in elucidating the epigenomes of these cells, including changes to the transcriptome, DNA methylome, and chromatin architecture. We will also discuss trained immunity in AM as a new and emerging field in TB research and provide some perspectives for the translational potential of targeting host epigenetics as an alternative TB therapy.

Stem cell-derived porcine macrophages as a new platform for studying host-pathogen interactions

BACKGROUND

Infectious diseases of farmed and wild animals pose a recurrent threat to food security and human health. The macrophage, a key component of the innate immune system, is the first line of defence against many infectious agents and plays a major role in shaping the adaptive immune response. However, this phagocyte is a target and host for many pathogens. Understanding the molecular basis of interactions between macrophages and pathogens is therefore crucial for the development of effective strategies to combat important infectious diseases.

RESULTS

We explored how porcine pluripotent stem cells (PSCs) can provide a limitless in vitro supply of genetically and experimentally tractable macrophages. Porcine PSC-derived macrophages (PSCdMs) exhibited molecular and functional characteristics of ex vivo primary macrophages and were productively infected by pig pathogens, including porcine reproductive and respiratory syndrome virus (PRRSV) and African swine fever virus (ASFV), two of the most economically important and devastating viruses in pig farming. Moreover, porcine PSCdMs were readily amenable to genetic modification by CRISPR/Cas9 gene editing applied either in parental stem cells or directly in the macrophages by lentiviral vector transduction.

CONCLUSIONS

We show that porcine PSCdMs exhibit key macrophage characteristics, including infection by a range of commercially relevant pig pathogens. In addition, genetic engineering of PSCs and PSCdMs affords new opportunities for functional analysis of macrophage biology in an important livestock species. PSCs and differentiated derivatives should therefore represent a useful and ethical experimental platform to investigate the genetic and molecular basis of host-pathogen interactions in pigs, and also have wider applications in livestock.

Effect of co-positivity for brucellosis and tuberculosis on milk yield and fertility of Holstein cows

This study aimed to determine whether cows detected as tuberculosis (bTB) reactors and seropositive to brucellosis (bBR), as well as co-positive to bBR and bTB (bBR-bTB) and with a complete lactation before slaughter, were associated with reduced milk production and fertility. A total of 8068 productive and reproductive records of high-yielding Holstein cows from a single large dairy herd with a high prevalence of bTB and bBR were collected from 2012 to 2015. Lactation derived either from calving (n = 6019) or hormonally induced lactation (n = 2049), and all cows received growth hormone throughout lactation. For cows not induced into lactation, pregnancy rate to first service for healthy cows (C; 26.6%) was higher (P < 0.01) than bBR (15.2%), bTB (15.8%), and bBR-bTB (1.3%) cows. For induced cows, pregnancy rate to first service did not differ significantly among C, bBR, and bTB (14.5-17.3%) cows, but the percentage success of first service was extremely low (1.3%; P < 0.01) in bBR-bTB cows. Services per pregnancy (only pregnant cows) were lowest for C (3.3 ± 2.9; P < 0.01) and highest (6.4 ± 3.4) for bBR-bTB non-induced cows. This variable was lowest for C (2.9 ± 2.5; P < 0.01) and highest for bBR-bTB non-induced cows (6.3 ± 3.1). Pregnancy rate to all services did not differed for C (79.5%), bBR (76.7%), and bTB (75.9%) but was lower (58.9%; P < 0.01) for bBR-bTB non-induced cows. For induced cows this variable was highest for bBR (53.3%) and lowest for bBR-bTB (34.1%; P < 0.01) non-induced cows. 305-d milk production was increased by 4%, and total milk yield by 7% in TB-positive cows compared to that of the negative cows non-induced hormonally into lactation. This study showed the negative impact of the co-positivity for bTB and bBR on the reproductive efficiency of Holstein cows, although positive bTB and bBR tests enhanced milk yield.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

Background

Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection.

Results

The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively.

Conclusions

The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07643-w.

Impacts of Epigenetic Processes on the Health and Productivity of Livestock

The dynamic changes in the epigenome resulting from the intricate interactions of genetic and environmental factors play crucial roles in individual growth and development. Numerous studies in plants, rodents, and humans have provided evidence of the regulatory roles of epigenetic processes in health and disease. There is increasing pressure to increase livestock production in light of increasing food needs of an expanding human population and environment challenges, but there is limited related epigenetic data on livestock to complement genomic information and support advances in improvement breeding and health management. This review examines the recent discoveries on epigenetic processes due to DNA methylation, histone modification, and chromatin remodeling and their impacts on health and production traits in farm animals, including bovine, swine, sheep, goat, and poultry species. Most of the reports focused on epigenome profiling at the genome-wide or specific genic regions in response to developmental processes, environmental stressors, nutrition, and disease pathogens. The bulk of available data mainly characterized the epigenetic markers in tissues/organs or in relation to traits and detection of epigenetic regulatory mechanisms underlying livestock phenotype diversity. However, available data is inadequate to support gainful exploitation of epigenetic processes for improved animal health and productivity management. Increased research effort, which is vital to elucidate how epigenetic mechanisms affect the health and productivity of livestock, is currently limited due to several factors including lack of adequate analytical tools. In this review, we (1) summarize available evidence of the impacts of epigenetic processes on livestock production and health traits, (2) discuss the application of epigenetics data in livestock production, and (3) present gaps in livestock epigenetics research. Knowledge of the epigenetic factors influencing livestock health and productivity is vital for the management and improvement of livestock productivity.

Recent advances in non-specific immune memory against bovine tuberculosis

Bovine tuberculosis is an important worldwide disease mainly related to cattle, although it also affects other mammals, including humans. In recent years, there have been considerable advances in the knowledge of the immune response mechanisms underlying the interaction of Mycobacterium bovis, the main agent of bovine tuberculosis, with its hosts. In this review we describe the most recent findings on the cattle immune response to M. bovis, particularly regarding trained innate immune responses and γδ T cells, that could support the development of vaccines and diagnostic tools to control this disease.

Genome-Wide Histone Modifications and CTCF Enrichment Predict Gene Expression in Sheep Macrophages

Alveolar macrophages function in innate and adaptive immunity, wound healing, and homeostasis in the lungs dependent on tissue-specific gene expression under epigenetic regulation. The functional diversity of tissue resident macrophages, despite their common myeloid lineage, highlights the need to study tissue-specific regulatory elements that control gene expression. Increasing evidence supports the hypothesis that subtle genetic changes alter sheep macrophage response to important production pathogens and zoonoses, for example, viruses like small ruminant lentiviruses and bacteria like Coxiella burnetii. Annotation of transcriptional regulatory elements will aid researchers in identifying genetic mutations of immunological consequence. Here we report the first genome-wide survey of regulatory elements in any sheep immune cell, utilizing alveolar macrophages. We assayed histone modifications and CTCF enrichment by chromatin immunoprecipitation with deep sequencing (ChIP-seq) in two sheep to determine cis-regulatory DNA elements and chromatin domain boundaries that control immunity-related gene expression. Histone modifications included H3K4me3 (denoting active promoters), H3K27ac (active enhancers), H3K4me1 (primed and distal enhancers), and H3K27me3 (broad silencers). In total, we identified 248,674 reproducible regulatory elements, which allowed assignment of putative biological function in macrophages to 12% of the sheep genome. Data exceeded the FAANG and ENCODE standards of 20 million and 45 million useable fragments for narrow and broad marks, respectively. Active elements showed consensus with RNA-seq data and were predictive of gene expression in alveolar macrophages from the publicly available Sheep Gene Expression Atlas. Silencer elements were not enriched for expressed genes, but rather for repressed developmental genes. CTCF enrichment enabled identification of 11,000 chromatin domains with mean size of 258 kb. To our knowledge, this is the first report to use immunoprecipitated CTCF to determine putative topological domains in sheep immune cells. Furthermore, these data will empower phenotype-associated mutation discovery since most causal variants are within regulatory elements.