Changes in H3K27ac at Gene Regulatory Regions in Porcine Alveolar Macrophages Following LPS or PolyIC Exposure

Porcine cis-acting lnc-CAST positively regulates CXCL8 expression through histone H3K27ac

The chemokine CXCL8, also known as the neutrophil chemotactic factor, plays a crucial role in mediating inflammatory responses and managing cellular immune reactions during viral infections. Porcine reproductive and respiratory syndrome virus (PRRSV) primarily infects pulmonary alveolar macrophages (PAMs), leading to acute pulmonary infections. In this study, we explored a novel long non-coding RNA (lncRNA), termed lnc-CAST, situated within the Cxcl8 gene locus. This lncRNA was found to be highly expressed in porcine macrophages. We observed that both lnc-CAST and CXCL8 were significantly upregulated in PAMs following PRRSV infection, and after treatments with lipopolysaccharide (LPS) or lipoteichoic acid (LTA). Furthermore, we noticed a concurrent upregulation of lnc-CAST and CXCL8 expression in lungs of PRRSV-infected pigs. We then determined that lnc-CAST positively influenced CXCL8 expression in PAMs. Overexpression of lnc-CAST led to an increase in CXCL8 production, which in turn enhanced the migration of epithelial cells and the recruitment of neutrophils. Conversely, inhibiting lnc-CAST expression resulted in reduced CXCL8 production in PAMs, leading to decreased migration levels of epithelial cells and neutrophils. From a mechanistic perspective, we found that lnc-CAST, localized in the nucleus, facilitated the enrichment of histone H3K27ac in CXCL8 promoter region, thereby stimulating CXCL8 transcription in a cis-regulatory manner. In conclusion, our study underscores the pivotal critical role of lnc-CAST in regulating CXCL8 production, offering valuable insights into chemokine regulation and lung damage during PRRSV infection.

Neutrophils secrete exosome-associated DNA to resolve sterile acute inflammation

Acute inflammation, characterized by a rapid influx of neutrophils, is a protective response that can lead to chronic inflammatory diseases when left unresolved. Secretion of LTB 4 -containing exosomes is required for effective neutrophil infiltration during inflammation. In this study, we show that neutrophils release nuclear DNA in a non-lytic, rapid, and repetitive manner, via a mechanism distinct from suicidal NET release and cell death. The packaging of nuclear DNA occurs in the lumen of nuclear envelope (NE)-derived multivesicular bodies (MVBs) that harbor the LTB 4 synthesizing machinery and is mediated by the lamin B receptor (LBR) and chromatin decondensation. Disruption of secreted exosome-associated DNA (SEAD) in a model of sterile inflammation in mouse skin amplifies and prolongs the presence of neutrophils, impeding the onset of resolution. Together, these findings advance our understanding of neutrophil functions during inflammation and the physiological significance of NETs, with implications for novel treatments for inflammatory disorders.

Genetic variation is a key determinant of chromatin accessibility and drives differences in the regulatory landscape of C57BL/6J and 129S1/SvImJ mice

Most common genetic variants associated with disease are located in non-coding regions of the genome. One mechanism by which they function is through altering transcription factor (TF) binding. In this study, we explore how genetic variation is connected to differences in the regulatory landscape of livers from C57BL/6J and 129S1/SvImJ mice fed either chow or a high-fat diet. To identify sites where regulatory variation affects TF binding and nearby gene expression, we employed an integrative analysis of H3K27ac ChIP-seq (active enhancers), ATAC-seq (chromatin accessibility) and RNA-seq (gene expression). We show that, across all these assays, the genetically driven (i.e. strain-specific) differences in the regulatory landscape are more pronounced than those modified by diet. Most notably, our analysis revealed that differentially accessible regions (DARs, N = 29635, FDR < 0.01 and fold change > 50%) are almost always strain-specific and enriched with genetic variation. Moreover, proximal DARs are highly correlated with differentially expressed genes. We also show that TF binding is affected by genetic variation, which we validate experimentally using ChIP-seq for TCF7L2 and CTCF. This study provides detailed insights into how non-coding genetic variation alters the gene regulatory landscape, and demonstrates how this can be used to study the regulatory variation influencing TF binding.

Sulforaphane impedes mitochondrial reprogramming and histone acetylation in polarizing M1 (LPS) macrophages

M1 (LPS) macrophages are characterized by a high expression of pro-inflammatory mediators, and distinct metabolic features that comprise increased glycolysis, a broken TCA cycle, or impaired OXPHOS with augmented mitochondrial ROS production. This study investigated whether the phytochemical sulforaphane (Sfn) influences mitochondrial reprogramming during M1 polarization, as well as to what extent this can contribute to Sfn-mediated inhibition of M1 marker expression in murine macrophages. The use of extracellular flux-, metabolite-, and immunoblot analyses as well as fluorescent dyes indicative for mitochondrial morphology, membrane potential or superoxide production, demonstrated that M1 (LPS/Sfn) macrophages maintain an unbroken TCA cycle, higher OXPHOS rate, boosted fusion dynamics, lower membrane potential, and less superoxide production in their mitochondria when compared to control M1 (LPS) cells. Sustained OXPHOS and TCA activity but not the concomitantly observed high dependency on fatty acids as fuel appeared necessary for M1 (LPS/Sfn) macrophages to reduce expression of nos2, il1β, il6 and tnfα. M1 (LPS/Sfn) macrophages also displayed lower nucleo/cytosolic acetyl-CoA levels in association with lower global and site-specific histone acetylation at selected pro-inflammatory gene promoters than M1 (LPS), evident in colorimetric coupled enzyme assays, immunoblot and ChIP-qPCR analyses, respectively. Supplementation with acetate or citrate was able to rescue both histone acetylation and mRNA expression of the investigated M1 marker genes in Sfn-treated cells. Overall, Sfn preserves mitochondrial functionality and restricts indispensable nuclear acetyl-CoA for histone acetylation and M1 marker expression in LPS-stimulated macrophages.

Transcriptomic analysis of intestinal organoids, derived from pigs divergent in feed efficiency, and their response to Escherichia coli

BACKGROUND

There is increasing interest in using intestinal organoids to study complex traits like feed efficiency (FE) and host-microbe interactions. The aim of this study was to investigate differences in the molecular phenotype of organoids derived from pigs divergent for FE as well as their responses to challenge with adherent and invasive Escherichia coli (E. coli).

RESULTS

Colon and ileum tissue from low and high FE pigs was used to generate 3D organoids and two dimensional (2D) monolayers of organoid cells for E. coli challenge. Genome-wide gene expression was used to investigate molecular differences between pigs that were phenotypically divergent for FE and to study the difference in gene expression after challenge with E. coli. We showed, (1) minor differences in gene expression of colon organoids from pigs with low and high FE phenotypes, (2) that an E. coli challenge results in a strong innate immune gene response in both colon and ileum organoids, (3) that the immune response seems to be less pronounced in the colon organoids of high FE pigs and (4) a slightly stronger immune response was observed in ileum than in colon organoids.

CONCLUSIONS

These findings demonstrate the potential for using organoids to gain insights into complex biological mechanisms such as FE.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Intranasal administration of Acinetobacter lwoffii in a murine model of asthma induces IL-6-mediated protection associated with cecal microbiota changes

BACKGROUND

Early-life exposure to certain environmental bacteria including Acinetobacter lwoffii (AL) has been implicated in protection from chronic inflammatory diseases including asthma later in life. However, the underlying mechanisms at the immune-microbe interface remain largely unknown.

METHODS

The effects of repeated intranasal AL exposure on local and systemic innate immune responses were investigated in wild-type and Il6-/- , Il10-/- , and Il17-/- mice exposed to ovalbumin-induced allergic airway inflammation. Those investigations were expanded by microbiome analyses. To assess for AL-associated changes in gene expression, the picture arising from animal data was supplemented by in vitro experiments of macrophage and T-cell responses, yielding expression and epigenetic data.

RESULTS

The asthma preventive effect of AL was confirmed in the lung. Repeated intranasal AL administration triggered a proinflammatory immune response particularly characterized by elevated levels of IL-6, and consequently, IL-6 induced IL-10 production in CD4+ T-cells. Both IL-6 and IL-10, but not IL-17, were required for asthma protection. AL had a profound impact on the gene regulatory landscape of CD4+ T-cells which could be largely recapitulated by recombinant IL-6. AL administration also induced marked changes in the gastrointestinal microbiome but not in the lung microbiome. By comparing the effects on the microbiota according to mouse genotype and AL-treatment status, we have identified microbial taxa that were associated with either disease protection or activity.

CONCLUSION

These experiments provide a novel mechanism of Acinetobacter lwoffii-induced asthma protection operating through IL-6-mediated epigenetic activation of IL-10 production and with associated effects on the intestinal microbiome.

Predicting enhancer-promoter interaction based on epigenomic signals

Introduction: The physical interactions between enhancers and promoters are often involved in gene transcriptional regulation. High tissue-specific enhancer-promoter interactions (EPIs) are responsible for the differential expression of genes. Experimental methods are time-consuming and labor-intensive in measuring EPIs. An alternative approach, machine learning, has been widely used to predict EPIs. However, most existing machine learning methods require a large number of functional genomic and epigenomic features as input, which limits the application to different cell lines. Methods: In this paper, we developed a random forest model, HARD (H3K27ac, ATAC-seq, RAD21, and Distance), to predict EPI using only four types of features. Results: Independent tests on a benchmark dataset showed that HARD outperforms other models with the fewest features. Discussion: Our results revealed that chromatin accessibility and the binding of cohesin are important for cell-line-specific EPIs. Furthermore, we trained the HARD model in the GM12878 cell line and performed testing in the HeLa cell line. The cross-cell-lines prediction also performs well, suggesting it has the potential to be applied to other cell lines.

Inhibition of EZH2 Causes Retrotransposon Derepression and Immune Activation in Porcine Lung Alveolar Macrophages

Alveolar macrophages (AMs) form the first defense line against various respiratory pathogens, and their immune response has a profound impact on the outcome of respiratory infection. Enhancer of zeste homolog 2 (EZH2), which catalyzes the trimethylation of H3K27 for epigenetic repression, has gained increasing attention for its immune regulation function, yet its exact function in AMs remains largely obscure. Using porcine 3D4/21 AM cells as a model, we characterized the transcriptomic and epigenomic alterations after the inhibition of EZH2. We found that the inhibition of EZH2 causes transcriptional activation of numerous immune genes and inhibits the subsequent infection by influenza A virus. Interestingly, specific families of transposable elements, particularly endogenous retrovirus elements (ERVs) and LINEs which belong to retrotransposons, also become derepressed. While some of the derepressed ERV families are pig-specific, a few ancestral families are known to be under EZH2-mediated repression in humans. Given that derepression of ERVs can promote innate immune activation through "viral mimicry", we speculate that ERVs may also contribute to the coinciding immune activation in AMs after the inhibition of EZH2. Overall, this study improves the understanding of the EZH2-related immune regulation in AMs and provides novel insights into the epigenetic regulation of retrotransposons in pigs.

Assessment of DNA methylation in porcine immune cells reveals novel regulatory elements associated with cell-specific gene expression and immune capacity traits

Background

Genetics studies in the porcine immune system have enhanced selection practices for disease resistance phenotypes and increased the efficacy of porcine models in biomedical research; however limited functional annotation of the porcine immunome has hindered progress on both fronts. Among epigenetic mechanisms that regulate gene expression, DNA methylation is the most ubiquitous modification made to the DNA molecule and influences transcription factor binding as well as gene and phenotype expression. Human and mouse DNA methylation studies have improved mapping of regulatory elements in these species, but comparable studies in the pig have been limited in scope.

Results

We performed whole-genome bisulfite sequencing to assess DNA methylation patterns in nine pig immune cell populations: CD21⁺ and CD21⁻ B cells, four T cell fractions (CD4⁺, CD8⁺, CD8⁺CD4⁺, and SWC6γδ⁺), natural killer and myeloid cells, and neutrophils. We identified 54,391 cell differentially methylated regions (cDMRs), and clustering by cDMR methylation rate grouped samples by cell lineage. 32,737 cDMRs were classified as cell lowly methylated regions (cLMRs) in at least one cell type, and cLMRs were broadly enriched in genes and regions of intermediate CpG density. We observed strong correlations between differential methylation and expression across immune cell populations, with cell-specific low methylation disproportionately impacting genes exhibiting enriched gene expression in the same cell type. Motif analysis of cLMRs revealed cell type-specific enrichment of transcription factor binding motifs, indicating that cell-specific methylation patterns may influence accessibility by trans-acting factors. Lastly, cDMRs were enriched for immune capacity GWAS SNPs, and many such overlaps occurred within genes known to influence immune cell development and function (CD8B, NDRG1).

Conclusion

Our DNA methylation data improve functional annotation of the porcine genome through characterization of epigenomic regulatory patterns that contribute to immune cell identity and function, and increase the potential for identifying mechanistic links between genotype and phenotype.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08773-5.

PFKM inhibits doxorubicin-induced cardiotoxicity by enhancing oxidative phosphorylation and glycolysis

Heart failure (HF) is a global pandemic which affects about 26 million people. PFKM (Phosphofructokinase, Muscle), catalyzing the phosphorylation of fructose-6-phosphate, plays a very important role in cardiovascular diseases. However, the effect of PFKM in glycolysis and HF remains to be elucidated. H9c2 rat cardiomyocyte cells were treated with doxorubicin (DOX) to establish injury models, and the cell viability, apoptosis and glycolysis were measured. Quantitative reverse transcription-polymerase chain reaction (RT-PCR) and immunoblotting were used for gene expression. DOX treatment significantly inhibited PFKM expression in H9c2 cells. Overexpression of PFKM inhibited DOX-induced cell apoptosis and DOX-decreased glycolysis and oxidative phosphorylation (OXPHOS), while silencing PFKM promoted cell apoptosis and inhibited glycolysis and OXPHOS in H9c2 cells. Moreover, PFKM regulated DOX-mediated cell viability and apoptosis through glycolysis pathway. Mechanism study showed that histone deacetylase 1 (HDAC1) inhibited H3K27ac-induced transcription of PFKM in DOX-treated cells and regulated glycolysis. PFKM could inhibit DOX-induced cardiotoxicity by enhancing OXPHOS and glycolysis, which might benefit us in developing novel therapeutics for prevention or treatment of HF.

Study of microRNA expression in Salmonella Typhimurium-infected porcine ileum reveals miR-194a-5p as an important regulator of the TLR4-mediated inflammatory response

Infection with Salmonella Typhimurium (S. Typhimurium) is a common cause of food-borne zoonosis leading to acute gastroenteritis in humans and pigs, causing economic losses to producers and farmers, and generating a food security risk. In a previous study, we demonstrated that S. Typhimurium infection produces a severe transcriptional activation of inflammatory processes in ileum. However, little is known regarding how microRNAs regulate this response during infection. Here, small RNA sequencing was used to identify 28 miRNAs differentially expressed (DE) in ileum of S. Typhimurium-infected pigs, which potentially regulate 14 target genes involved in immune system processes such as regulation of cytokine production, monocyte chemotaxis, or cellular response to interferon gamma. Using in vitro functional and gain/loss of function (mimics/CRISPR-Cas system) approaches, we show that porcine miR-194a-5p (homologous to human miR-194-5p) regulates TLR4 gene expression, an important molecule involved in pathogen virulence, recognition and activation of innate immunity in Salmonella infection.

Genome-wide association study of disease resilience traits from a natural polymicrobial disease challenge model in pigs identifies the importance of the major histocompatibility complex region

Infectious diseases cause tremendous financial losses in the pork industry, emphasizing the importance of disease resilience, which is the ability of an animal to maintain performance under disease. Previously, a natural polymicrobial disease challenge model was established, in which pigs were challenged in the late nursery phase by multiple pathogens to maximize expression of genetic differences in disease resilience. Genetic analysis found that performance traits in this model, including growth rate, feed and water intake, and carcass traits, as well as clinical disease phenotypes, were heritable and could be selected for to increase disease resilience of pigs. The objectives of the current study were to identify genomic regions that are associated with disease resilience in this model, using genome-wide association studies and fine-mapping methods, and to use gene set enrichment analyses to determine whether genomic regions associated with disease resilience are enriched for previously published quantitative trait loci, functional pathways, and differentially expressed genes subject to physiological states. Multiple quantitative trait loci were detected for all recorded performance and clinical disease traits. The major histocompatibility complex region was found to explain substantial genetic variance for multiple traits, including for growth rate in the late nursery (12.8%) and finisher (2.7%), for several clinical disease traits (up to 2.7%), and for several feeding and drinking traits (up to 4%). Further fine mapping identified 4 quantitative trait loci in the major histocompatibility complex region for growth rate in the late nursery that spanned the subregions for class I, II, and III, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects, explaining 0.8–27.1% of genetic variance for growth rate and for multiple clinical disease traits. This single-nucleotide polymorphism was located in the enhancer of TRIM39 gene, which is involved in innate immune response. The major histocompatibility complex region was pleiotropic for growth rate in the late nursery and finisher, and for treatment and mortality rates. Growth rate in the late nursery showed strong negative genetic correlations in the major histocompatibility complex region with treatment or mortality rates (−0.62 to −0.85) and a strong positive genetic correlation with growth rate in the finisher (0.79). Gene set enrichment analyses found genomic regions associated with resilience phenotypes to be enriched for previously identified disease susceptibility and immune capacity quantitative trait loci, for genes that were differentially expressed following bacterial or virus infection and immune response, and for gene ontology terms related to immune and inflammatory response. In conclusion, the major histocompatibility complex and other quantitative trait loci that harbor immune-related genes were identified to be associated with disease resilience traits in a large-scale natural polymicrobial disease challenge. The major histocompatibility complex region was pleiotropic for growth rate under challenge and for clinical disease traits. Four quantitative trait loci were identified across the class I, II, and III subregions of the major histocompatibility complex for nursery growth rate under challenge, with 1 single-nucleotide polymorphism in the major histocompatibility complex class I subregion capturing the largest effects. The major histocompatibility complex and other quantitative trait loci identified play an important role in host response to infectious diseases and can be incorporated in selection to improve disease resilience, in particular the identified single-nucleotide polymorphism in the major histocompatibility complex class I subregion.

Epigenomic and transcriptomic analyses reveal differences between low-grade inflammation and severe exhaustion in LPS-challenged murine monocytes

Emerging studies suggest that monocytes can be trained by bacterial endotoxin to adopt distinct memory states ranging from low-grade inflammation to immune exhaustion. While low-grade inflammation may contribute to the pathogenesis of chronic diseases, exhausted monocytes with pathogenic and immune-suppressive characteristics may underlie the pathogenesis of polymicrobial sepsis including COVID-19. However, detailed processes by which the dynamic adaption of monocytes occur remain poorly understood. Here we exposed murine bone-marrow derived monocytes to chronic lipopolysaccharide (LPS) stimulation at low-dose or high-dose, as well as a PBS control. The cells were profiled for genome-wide H3K27ac modification and gene expression. The gene expression of TRAM-deficient and IRAK-M-deficient monocytes with LPS exposure was also analyzed. We discover that low-grade inflammation preferentially utilizes the TRAM-dependent pathway of TLR4 signaling, and induces the expression of interferon response genes. In contrast, high dose LPS uniquely upregulates exhaustion signatures with metabolic and proliferative pathways. The extensive differences in the epigenomic landscape between low-dose and high-dose conditions suggest the importance of epigenetic regulations in driving differential responses. Our data provide potential targets for future mechanistic or therapeutic studies.

Reference Transcriptomes of Porcine Peripheral Immune Cells Created Through Bulk and Single-Cell RNA Sequencing

Pigs are a valuable human biomedical model and an important protein source supporting global food security. The transcriptomes of peripheral blood immune cells in pigs were defined at the bulk cell-type and single cell levels. First, eight cell types were isolated in bulk from peripheral blood mononuclear cells (PBMCs) by cell sorting, representing Myeloid, NK cells and specific populations of T and B-cells. Transcriptomes for each bulk population of cells were generated by RNA-seq with 10,974 expressed genes detected. Pairwise comparisons between cell types revealed specific expression, while enrichment analysis identified 1,885 to 3,591 significantly enriched genes across all 8 cell types. Gene Ontology analysis for the top 25% of significantly enriched genes (SEG) showed high enrichment of biological processes related to the nature of each cell type. Comparison of gene expression indicated highly significant correlations between pig cells and corresponding human PBMC bulk RNA-seq data available in Haemopedia. Second, higher resolution of distinct cell populations was obtained by single-cell RNA-sequencing (scRNA-seq) of PBMC. Seven PBMC samples were partitioned and sequenced that produced 28,810 single cell transcriptomes distributed across 36 clusters and classified into 13 general cell types including plasmacytoid dendritic cells (DC), conventional DCs, monocytes, B-cell, conventional CD4 and CD8 αβ T-cells, NK cells, and γδ T-cells. Signature gene sets from the human Haemopedia data were assessed for relative enrichment in genes expressed in pig cells and integration of pig scRNA-seq with a public human scRNA-seq dataset provided further validation for similarity between human and pig data. The sorted porcine bulk RNAseq dataset informed classification of scRNA-seq PBMC populations; specifically, an integration of the datasets showed that the pig bulk RNAseq data helped define the CD4CD8 double-positive T-cell populations in the scRNA-seq data. Overall, the data provides deep and well-validated transcriptomic data from sorted PBMC populations and the first single-cell transcriptomic data for porcine PBMCs. This resource will be invaluable for annotation of pig genes controlling immunogenetic traits as part of the porcine Functional Annotation of Animal Genomes (FAANG) project, as well as further study of, and development of new reagents for, porcine immunology.

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.

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.