C/EBPδ-induced epigenetic changes control the dynamic gene transcription of S100a8 and S100a9
Abstract
The proinflammatory alarmins S100A8 and S100A9 are among the most abundant proteins in neutrophils and monocytes but are completely silenced after differentiation to macrophages. The molecular mechanisms of the extraordinarily dynamic transcriptional regulation of S100a8 and S100a9 genes, however, are only barely understood. Using an unbiased genome-wide CRISPR/Cas9 knockout (KO)-based screening approach in immortalized murine monocytes, we identified the transcription factor C/EBPδ as a central regulator of S100a8 and S100a9 expression. We showed that S100A8/A9 expression and thereby neutrophil recruitment and cytokine release were decreased in C/EBPδ KO mice in a mouse model of acute lung inflammation. S100a8 and S100a9 expression was further controlled by the C/EBPδ antagonists ATF3 and FBXW7. We confirmed the clinical relevance of this regulatory network in subpopulations of human monocytes in a clinical cohort of cardiovascular patients. Moreover, we identified specific C/EBPδ-binding sites within S100a8 and S100a9 promoter regions, and demonstrated that C/EBPδ-dependent JMJD3-mediated demethylation of H3K27me3 is indispensable for their expression. Overall, our work uncovered C/EBPδ as a novel regulator of S100a8 and S100a9 expression. Therefore, C/EBPδ represents a promising target for modulation of inflammatory conditions that are characterized by S100a8 and S100a9 overexpression.
Editor's evaluation
The study uses an elegant CRISPR/Cas9 screening approach in a myeloid cell line to identify the underlying regulators of the alarmins S100A8 and S100A9, which amplify inflammation. This approach identified the transcription factor C/EBP-δ as a regulator of S100A8 and S100A9 expression in the myeloid cell line and also showed a correlation between the expression levels of C/EBP-δ and the alarmins in patient samples of peripheral blood mononuclear cells. Furthermore, the authors also validate their findings in primary monocytes using mice with genetic C/EBP-δ deletion. This work will be of significant interest to researchers studying the regulation of immune responses and inflammation, and it also highlights how unbiased CRISPR/Cas9 screening can lead to novel mechanistic insights in myeloid cells.
https://doi.org/10.7554/eLife.75594.sa0Introduction
As the first line of immune defence, both monocytes and neutrophils are important for the modulation of the innate immune response. To amplify the immune response at sites of inflammation, the activation of further immune cells is required, mediated by the release of signaling molecules such as chemokines and DAMPs (damage-associated molecular patterns). The two members of the S100 family, S100A8 and S100A9, also termed myeloid-related proteins 8 and 14 (MRP8 and MRP14), respectively, belong to the group of DAMPs or so-called alarmins. Their primary expression is referred to myeloid lineage-derived cells, particularly neutrophils and monocytes, where S100A8 and S100A9 are predominantly present as a heterodimeric complex, also called calprotectin (Austermann et al., 2018).
Intracellularly, S100A8/A9 complexes represent up to 40% of the soluble protein content in neutrophils and about 5% in monocytes (Hessian et al., 1993). However, in mature macrophages, protein and mRNA expression of these factors is completely downregulated. This data indicates that expression of S100A8 and S100A9 is controlled by the most dynamic promoters in myeloid cells. The S100A8/A9 complex interacts with the cytoskeleton in a calcium-dependent manner. Calcium-induced (S100A8/A9)2 tetramer promotes tubulin polymerization and microtubule bundling, thereby affecting transendothelial migration of phagocytes (Leukert et al., 2006). During inflammation or tissue damage, S100A8/A9 is actively secreted by neutrophils and monocytes, and represents the most abundant DAMP/alarmin activating inflammatory processes in infection, cancer, autoimmunity, and cardiovascular diseases. The S100A8/A9 complex is recognized by Toll-like receptor 4 (TLR4), which leads to the production of proinflammatory cytokines and chemokines (Fassl et al., 2015). Accordingly, S100A9 knockout (KO) mice exhibit decreased pathogenic outcomes in several mouse models of disease, such as sepsis (Vogl et al., 2007), autoimmune disease (Loser et al., 2010) or arthritis (van Lent et al., 2008). In addition, S100A8 and S100A9 are highly abundant during infectious diseases and exhibit anti-microbial activities. The S100A8/A9 complex plays a crucial role in host defence against bacterial and fungal pathogens by sequestering manganese and zinc ions which compete with high affinity bacterial transporters to import these essential nutrient metals (Kehl-Fie and Skaar, 2010; Kehl-Fie et al., 2011). In contrast to the proinflammatory role of S100A8/A9, regulatory functions in terms of hyporesponsiveness in phagocytes, resembling a classical endotoxin-induced tolerance, have also been described (Freise et al., 2019).
In humans, S100A8/A9 is the most abundant alarmin in many clinically relevant diseases, and is closely associated with disease activity in rheumatoid arthritis (RA), inflammatory bowel disease, sepsis, cardiovascular diseases, multiple sclerosis, acute lung injury (ALI) and psoriasis (Foell et al., 2004). Altered S100A8/A9 expression has also been found in different cancer types, including gastric, colorectal, breast, lung, prostate and liver cancer (Cross et al., 2005). Despite the high expression in neutrophils and monocytes under inflammatory conditions, and the strong effects of S100A8 and S100A9 on disease activities, transcriptional mechanisms regulating these extreme dynamics of gene expression remain unclear. Identifying the mechanisms regulating S100a8 and S100a9 gene expression may open new insights into the pathological processes involving S100A8/A9 during inflammatory conditions.
So far, several potential transcription factors modulating S100a8 and S100a9 expression have been described (Kuruto-Niwa et al., 1998; Fujiu et al., 2011; Lee et al., 2012; Liu et al., 2016; Yang et al., 2017), but their functional relevance remains unresolved. Many of the stated studies used malignant immortalized cell lines or even cell models whose homologous primary cells do not express these genes at all.
To overcome difficulties of artificial expression and malignant cell lines we used ER (estrogen-regulated) Hoxb8 cells, estrogen dependent transiently immortalized myeloid precursor cells that can be differentiated to primary monocytes and granulocytes upon estrogen-withdrawal (Wang et al., 2006a), and show the physiologically high dynamics of S100A8 and S100A9 mRNA and protein expression during differentiation. In order to detect genes involved in the regulation of S100a8 and S100a9 expression in an unbiased manner, we used a mouse Genome-Scale CRISPR/Cas9 Knockout (GeCKO) library and screened for monocytes with reduced or absent S100A9 expression. We thereby identified the CCAAT/enhancer-binding protein-family member C/EBPδ as a direct transcriptional regulator of S100a8 and S100a9. Furthermore, we found that the epigenetic factor JMJD3 contributes to S100a8 and S100a9 expression in monocytes by erasure of the repressive histone mark H3K27me3 at S100a8 and S100a9 promoter regions. Moreover, we confirmed the biomedical relevance of this network a murine model of ALI and in specific monocyte subpopulations in a clinical cohort of patients with cardiovascular disease.
Results
Genome-wide CRISPR/Cas9 knockout screen reveals C/EBPδ as a regulatory factor of S100A9 expression
To detect genes involved in the regulation of S100a8 and S100a9 expression, we established the mouse GeCKO lentiviral pooled library designed in Cas9-expressing ER-Hoxb8 cells. The used library contained a large mixture of CRISPR sgRNA constructs, where six gRNAs per target gene increase efficiency and enable the analysis of the molecular effects of many thousand genes in one experiment. After infection of Cas9-expressing ER-Hoxb8 precursor cells with CRISPR library lentiviral particles, the cells were differentiated for 3 days in the presence of GM-CSF to induce S100A8 and S100A9 expression. Because we hypothesized that the parallel S100A8 and S100A9 expression is based on a common regulatory mechanism, we assumed that screening of one of the two alarmins was sufficient in the first step. Therefore, cells with no or low S100A9 expression were selected by sorting and considered as hits, whereas the remaining cells functioned as reference cells. To exclude phenotypes that were S100A9low/neg due to general differentiation defects, we pre-gated for CD11b+Ly-6C+ monocytes. DNA of sorted cell pools was purified and analysed by NGS (Figure 1A). Intracellular S100A9-FITC staining of precursor and differentiated Cas9 ER-Hoxb8 control monocytes was used as standard for definition of sorting gates. Differentiated Cas9-library ER-Hoxb8 monocytes showed a wider distribution among the gates, indicating the presence of S100A9low/neg expressing cells due to disruptions of regulatory genes caused by CRISPR/Cas9. A small amount of S100A9neg sorted cells were re-analysed by immunoblotting to validate S100A9 deficiency in this cell population (Figure 1B). Analysis of CRISPR KO library screen using the Model-based Analysis of Genome-wide CRISPR-Cas9 Knockout (MAGeCK) method (Li et al., 2014) resulted in a list of genes for which the respective gRNAs were enriched in the hits sample. The highest number of three gRNAs was found within the top 20 hits targeted Cebpd, a gene encoding for a member of the CCAAT/enhancer-binding protein family, C/EBPδ. Other gene hits, namely Manea, Htr1f, Ifit3, Phf8, Hand1, and Casp4, were targeted by two gRNAs (Figure 1C). From these, those genes known to be involved in gene regulation, such as the transcription factors Cebpd and Hand1, as well as the histone demethylase Phf8, were selected for validation. Another transcription factor in the top 20 hit list, Csrp1, was targeted by only one gRNA but was chosen for validation due to its known functions in gene regulatory mechanisms (Kamar et al., 2017). Targeting Phf8, Hand1, and Csrp1 only slightly influenced S100a8 and S100a9 expression, whereas C/EBPδ deficiency strongly decreased S100a8 and S100a9 levels in differentiated monocytes (Figure 1—figure supplement 1A). Interestingly, none of the transcription factors previously reported to target S100a8 and S100a9 (Kuruto-Niwa et al., 1998; Fujiu et al., 2011; Lee et al., 2012; Liu et al., 2016; Yang et al., 2017) were found within the hit list of our CRISPR KO library screen (Figure 1—source data 1). Nevertheless, to test the published candidate transcription factors ATF3, STAT3, KLF5, IRF7, and C/EBPβ for their effects on S100a8 and S100a9 regulation, we created single KO ER-Hoxb8 cell lines of each individual candidate transcription factor. Deficiency of none of the stated candidate transcription factors affects S100a8 and S1009 expression during monocyte differentiation, whereas C/EBPδ deficiency had a strong attenuating effect on S100a8 and S100a9 expression, as shown on day 2 (Figure 1—figure supplement 1B).
Decreased S100a8 and S100a9 expression in C/EBPδ KO monocytes
We confirmed extraordinarily high dynamics of S100a8 and S100a9 expression during monocyte/macrophage differentiation. ER-Hoxb8-derived monocytes showed about 590-fold increase in S100a8 mRNA expression and about 1800-fold increase of S100a9 mRNA expression on day 2 compared to day 0 of differentiation. At day 5 the S100a8 mRNA expression is already about 70-fold and the S100a9 mRNA expression about 110-fold downregulated compared to day 2 of differentiation (Figure 2A, B).
We confirmed the relevance of C/EBPδ for S100a8 and S100a9 expression by creating independent C/EBPδ-deficient ER-Hoxb8 cells from C/EBPδ KO mice. Not only on differentiation day 3, but already at the very beginning of differentiation, when S100a8 and S100a9 levels start to rise, C/EBPδ-deficient ER-Hoxb8 monocytes showed significantly reduced levels of both S100a8 (Figure 2A) and S100a9 (Figure 2B) mRNAs compared to wildtype (WT) controls. The same effect was detectable in C/EBPδ-deficient ER-Hoxb8 cells that were differentiated into the neutrophilic lineage (Figure 2—figure supplement 1A). Accordingly, Cebpd and S100a8 and S100a9 mRNAs were co-expressed in differentiating WT monocytes and neutrophils, supporting a mechanistic connection (Figure 2—figure supplement 1B, C). WT ER-Hoxb8 monocytes secreted significant S100A8/A9 protein amounts, whereas the supernatant of C/EBPδ KO cells contained up to 80% less S100A8/A9 (Figure 2C). Accordingly, S100a8 and S100a9 expression was decreased upon C/EBPδ deletion in primary bone marrow-derived monocytes (Figure 2D). To confirm that C/EBPδ deficiency affects alarmin expression and disease progression in vivo, we used a mouse model for acute lung inflammation. S100A8/A9 levels are locally increased in acute respiratory distress syndrome patients (Lorenz et al., 2008) and have been shown to play a role in neutrophil recruitment during ALI in mice (Chakraborty et al., 2017). Serum and bronchoalveolar lavage fluid (BALF) from LPS-exposed mice showed highly increased S100A8/A9 levels compared to control mice. In C/EBPδ KO mice, alarmin levels were systemically (serum) and locally (BALF) decreased compared to WT mice at baseline (NaCl) and after LPS exposure (Figure 2E, F), confirming our in vitro data. Accordingly, in BALF neutrophil recruitment (Figure 2G), as well as in IL-1α and IL-6, cytokine production was significantly decreased in LPS-exposed C/EBPδ KO mice in relation to WT littermates (TNF-α shows a similar trend, Figure 2H-J), which highlights the impact of S100A8/A9 on the disease outcome.
Inflammatory capacities and differentiation kinetics are unaffected in C/EBPδ KO monocytes
Although the proinflammatory molecule S100A8/A9 was strongly reduced in the C/EBPδ KO monocytes, these cells exhibited no general alterations of inflammatory functions, indicating a rather specific effect on S100a8 and S100a9 regulation than a general attenuation of inflammatory signaling due to C/EBPδ deficiency. Phagocytosis capacities, examined by using FITC-coupled Latex Beads (Figure 3A, B) and by using Staphylococcus aureus (Figure 3C, D), were even elevated in C/EBPδ KO monocytes. These observations were accompanied by enhanced gene expression of phagocytosis-related PRRs, such as Ptx3, CD209a, and Cd14, in C/EBPδ KO monocytes (Figure 3—figure supplement 1A, B, D). ROS production was not influenced by C/EBPδ deficiency (Figure 3E). Moreover, analysis of differentiation kinetics revealed no exorbitant differences in the quantities of CD11b+Ly-6C+ between WT and C/EBPδ KO cells, neither during differentiation of bone marrow-derived cells (BMCs) nor of ER-Hoxb8 cells (Figure 3F). With regard to monocyte differentiation and polarization capacities, we tested the effect of C/EBPδ deficiency in primary bone marrow-derived mouse monocytes (BMDMs). Whereas expression of M1-monocyte associated markers, such as Tnfa, Il6, Inos, Cd86, and Il1b, was only slightly decreased in LPS- and IFN-γ-treated C/EBPδ KO monocytes in vitro (Figure 3—figure supplement 2A), M2-associated Il10 expression was significantly decreased upon IL-4-stimulation in C/EBPδ KO cells confirming earlier results (Liu et al., 2003). Expression of Cd163 showed no significant effect in C/EBPδ KO cells (Figure 3—figure supplement 2B). Decreased IL-10 expression was also found in sera from LPS-treated and control C/EBPδ KO mice, whereas most other cytokine levels did not differ greatly between WT and C/EBPδ KO mice (Figure 2—figure supplement 2).
Enhanced C/EBPδ expression induces S100a8 and S100a9 expression
To test the impact of C/EBPδ induction on S100 alarmin expression, we infected C/EBPδ-deficient ER-Hoxb8 cells with lentiviral particles carrying a Tet-On system for doxycycline-inducible 3xFlag-C/EBPδ expression (Figure 4A). Doxycycline treatment led to expression of the fusion protein 3xFlag-C/EBPδ, as revealed by western blot analysis (Figure 4B) and by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in comparison to C/EBPδ-deficient cells (Figure 4C). Induction of 3xFlag-C/EBPδ upon doxycycline treatment led to increased S100a8 and S100a9 mRNA levels. Cebpd, S100a8, and S100a9 mRNA levels in doxycycline-treated TRE_3xFlag-C/EBPδ cells were comparable to WT cells at the same differentiation stage (Figure 4D), demonstrating a positive effect of C/EBPδ expression on S100a8 and S100a9 regulation. KO of ATF3, a known regulatory attenuator of Cebpd expression (Litvak et al., 2009), in ER-Hoxb8 monocytes led to S100a8 and S100a9 overexpression (Figure 4E, Figure 1—figure supplement 1B), especially during early stages of differentiation. ATF3 KO cells showed significantly elevated cebpd level, indicating a C/EBPδ-mediated effect on the expression of S100a8 and S100a9 (Figure 4F). In the next step, we created FBXW7-deficient monocytes. FBXW7 is another well-known attenuator of C/EBPδ expression (Balamurugan et al., 2013). Lack of this antagonist resulted in an even higher overexpression of S100a8 and S100a9 (Figure 4G) with huge increases of Cebpd levels (Figure 4H).
To confirm the biomedical relevance of the identified molecular network, we analysed the expression of these genes in peripheral blood mononuclear cells (PBMCs) and monocyte subpopulations of a subset of participants in the BioNRW Study (Witten et al., 2022). Here, we found upregulation of S100A8, S100A9, and CEBPD in PBMCs of stable coronary artery disease/myocardial infarction (sCAD/MI) cases, compared against controls (Figure 5A), together with a positive correlation of S100A8 and S100A9 expression with that of CEBPD in these cells (Figure 5B, C). Moreover, there was also significant upregulation of these three genes specifically in classical monocytes, compared to intermediate and non-classical monocyte subpopulations (Figure 5D, E, for further comparisons, see Figure 5—source data 1). A strong positive correlation between the expression of S100A8 and S100A9 and CEBPD in these monocyte subpopulations was found, suggesting that the expression of these genes is mainly associated with the subset of proinflammatory monocytes. Interestingly, we also found significant albeit milder, negative correlations between the expression of CEBPD and its antagonists FBXW7 and ATF3 in monocytes (Figure 5F, G).
C/EBPδ-binding sites within S100a8 and S100a9 promoter regions
Chromatin immunoprecipitation (ChIP) revealed 3xFlag-C/EBPδ binding on S100a8 and S100a9 promoter regions just before or within the predicted enhancers (Figure 6A). Co-transfection of HEK293T cells with GFP-expressing s100 reporter vectors, together with doxycycline-inducible 3xFlag-C/EBPδ vector (TRE_3xFlag-C/EBPδ) or its backbone lacking the 3xFlag-C/EBPδ construct (TRE_ctrl), was performed to further examine specific C/EBPδ binding (Figure 6B). Doxycycline treatment resulted in 3xFlag-C/EBPδ protein expression after 24 hr post-transfection in 3xFlag-C/EBPδ vector transfected cells (Figure 6C). Transfection of either S100a8 reporter construct (Figure 6D) or S100a9 reporter construct (Figure 6F) together with 3xFlag-C/EBPδ vector led to enhanced GFP expression upon doxycycline treatment in comparison to co-transfection with backbone plasmid (TRE_ctrl). Next, we modified predicted C/EBP-binding sites on S100 promoters by mutagenesis of the corresponding vectors. Again, co-transfection of mutated S100 reporter vectors and doxycycline-dependent 3xFlag-C/EBPδ vector was performed to analyse the relevance of specific C/EBPδ-binding sites. Two sites within the S100a8 promoter region, stated as site 2 and site 3 (Figure 6E), and one within the S100a9 promoter region, stated as site 4 (Figure 6G), caused a reduced or absent GFP expression upon co-transfection when deleted. These binding sites, in turn, were located within the S100a8 and S100a9 promoter regions where C/EBPδ binding was confirmed by ChIP (Figure 6A).
Epigenetic landscape on S100 promoter regions reflects S100a8 and S100a9 expression in monocytes
Regulation of gene expression relies on variable factors; among these are chromatin structure and epigenetic features. To measure changes in chromatin accessibility in monocytic progenitors and in S100a8- and S100a9-expressing monocytes, we performed ATAC-seq of precursor and differentiated WT and C/EBPδ KO ER-Hoxb8 cells. This revealed over 1000 regions with differential peaks in all comparisons (Figure 7A and Figure 7—figure supplement 1). Among the regions with significantly higher ATAC-seq reads in differentiated samples at day 3 of WT cells were the S100a8 and S100a9 promoter and enhancer locations. Interestingly, within these regions chromatin accessibility was significantly decreased in C/EBPδ KO in comparison to WT cells at day 3 of differentiation (Figure 7B, C). Consistent with the changes in chromatin accessibility at S100 promoter regions during differentiation, we also found changes in histone marks by ChIP. H3K27 acetylation (H3K27ac), a marker for active transcription, was increased at differentiation day 3 in monocytes over precursor cells at S100a8 (Figure 7D) and S100a9 loci (Figure 7E) in both, WT and C/EBPδ KO cells. In contrast, tri-methylated H3K27 (H3K27me3), associated with gene silencing, was overrepresented in precursor cells over differentiated cells at the same loci in WT cells, whereas H3K27me3 marks did not decrease over the course of differentiation in C/EBPδ KO cells. Accordingly, tri-methylated H3K27 was increased in C/EBPδ KO monocytes, compared to the WT counterparts (Figure 7F, G).
The histone demethylase JMJD3 drives S100a8 and S100a9 expression in dependency of C/EBPδ
Decreased S100a8 and S100a9 expression in C/EBPδ-deficient day 3 monocytes was mirrored in the epigenetic landscape by only slightly decreased H3K27ac level, but highly increased H3K27me3 level at S100 promoter regions. Erasure of tri-methylation and di-methylation at H3K27 is known to be catalysed by the histone demethylase JMJD3 (JmjC Domain-Containing Protein 3) (Xiang et al., 2007). We found decreased expression of Jmjd3 in differentiated C/EBPδ KO monocytes, compared to WT cells at the same stage (Figure 8A). Further, we used the potent JMJD3 inhibitor GSK-J4 (Kruidenier et al., 2012) to block H3K27 demethylation in differentiating ER-Hobx8 cells, and discovered significantly decreased S100a8 and S100a9 expression in GSK-J4-treated WT cells (Figure 8B). These mRNA quantities were comparable to untreated C/EBPδ-deficient monocytes, whereas the effects on S100a8 and S100a9 expression in GSK-J4-treated C/EBPδ-deficient monocytes, compared to the untreated counterparts, were negligible (Figure 7B). These effects of GSK-J4 on S100a8 and S100a9 expression are in line with increased H3K27me3 marks in GSK-J4-treated WT monocytes, compared to untreated WT cells on both, S100a8 (Figure 8C) and S100a9 (Figure 8D) promoter regions.
Discussion
The ER-Hoxb8 cell system serves as a substitute for murine primary cells of myeloid origin that can be differentiated into phagocytes, such as monocytes and neutrophils. It reflects reliably the differentiation of myeloid progenitor cells to bone marrow-derived monocytes and macrophages. Due to the high expression of S100A8 and S100A9 in the early days of culture around day 3, we describe these stage as monocytes, later days of differentiation with no or low expression of these S100 proteins as macrophages in analogy to the human system. This system allows comparison with in vitro differentiated primary cells (Wang et al., 2006a) and, therefore, provides an experimental cell model for analysis of S100A8 and S100A9 expression. Although the alarmins are regarded as key factors in various inflammatory conditions (Foell and Roth, 2004), cancer types (Cross et al., 2005), and cardiovascular diseases (Frangogiannis, 2019), little is known about their transcriptional regulation. The serum concentrations of alarmins correlate with disease severity and activity and, hence, they are reliable biomarkers for monitoring several inflammatory diseases (Foell et al., 2004; Ehrchen et al., 2009). The expression levels of S100a8 and S100a9 differ extremely during myeloid differentiation and the promoters of their genes represent probably one of the most dynamic regulatory elements in the myeloid lineage. Whereas both proteins are completely absent in myeloid precursor cells, they are highly expressed in monocytes and neutrophils, which suggests that highly dynamic regulatory mechanisms drive S100a8 and S100a9 expression.
The CRISPR/Cas9-mediated KO screening approach based on a lentiviral pooled library has been used so far to investigate various mechanisms, such as immunity-related pathways and cancer-modulating events (Kweon and Kim, 2018). In this study, our unbiased genome-wide screening approach allowed the identification of C/EBPδ as a factor involved in S100A9 regulation during murine monocyte differentiation. We further focused our investigations on Cebpd because this gene was in the top list of the robust rank aggregation (RRA) scores and showed the highest numbers of guide RNAs with efficient effects on S100A9 expression in our screening. This redundancy of independent parameters helped to distinguish true positive from false positive hits. Furthermore, a robust phenotype-of-interest, such as a clear S100A9 protein signal at day 3 of monocyte differentiation, allowed reliable negative selection in Cas9-library monocytes. Selection of remaining cells served as a reference control to distinguish true from false positives. The specificity of our selection procedure was confirmed at the protein level by western blot analysis of sorted cell populations. CRISPR/Cas9-based functional genomic screening has been reported to be highly specific, thereby causing fewer cases of false positives in direct comparison with knockdown analysis by RNA interference (Shalem et al., 2014). We were now able to identify a novel regulator of S100A8 and S100A9 using this unbiased method. By pre-gating on CD11b+Ly-6C+ monocytes, we revealed C/EBPδ as a specific and differentiation-independent regulator of S100A8 and S100A9, excluding pathways linked to general functions or development of phagocytes. Re-analysis of S100A9 expression by immunoblotting of sorted cell populations verified the gates set (Figure 1B). However, validation of four candidate genes from our hit list demonstrates the limitation of the CRISPR screening system. Not every targeted hit, such as Phf8, Hand1, and Csrp1, showed a prominent effect on S100a8 and S100a9 expression. Possibly, the slight differences are rather due to indirect effects on S100 expression than due to deficiency of actual direct regulators. Another interesting gene within our top 20 hit list is the transcription factor creb1 (Figure 1C) that is known to mediate cytokine signaling at least in human neutrophils and could therefore be tested for S100a8 and S100a9 regulation in future experiments (Mayer et al., 2013). Known regulators of S100A8 and S100A9 that are also essential for general differentiation, such as PU.1 (Xu et al., 2021), were very likely excluded due to our monocyte pre-gating strategy. Our data indicate that this screening technique seems not optimal enough to identify all factors of an interrelated regulatory protein complex, but it is sufficient for unbiased detection of single factors for gene regulation, as shown in our study.
We confirmed that the transcription factor C/EBPδ is a direct regulator of S100A8 and S100A9 in murine monocytes using independent approaches. Cebpd, and S100a8 and S100a9 were co-expressed in differentiating monocytes, and induction of C/EBPδ clearly showed that the expression of S100a8 and S100a9 was upregulated by the presence of C/EBPδ. This evidence was further supported by increased S100a8 and S100a9 levels caused by deletion of ATF3 and FBXW7, which are natural inhibitors of C/EBPδ. There is a parallel decrease of Cebpd, and S100a8 and S100a9 expression during later stages of monocyte/macrophage differentiation which may point to a functional relation of C/EBPδ expression and downregulation of S100 expression as well. However, the late expression kinetics of S100a8 and S100a9 in WT and C/EBPδ KO cells clearly point to the presence of additional relevant factors which were not addressed by our screening approach.
The specificity of our approach was further confirmed by the fact that deficiency of several transcription factors, such as STAT3, KLF5, IRF7, and C/EBPβ, described as S100A8 and S100A9 regulators in previous studies (Kuruto-Niwa et al., 1998; Fujiu et al., 2011; Lee et al., 2012; Liu et al., 2016; Yang et al., 2017), did neither affect S100a8 and S100a9 expression in our ER-Hoxb8 monocytes, nor were these listed as gene hits in our CRISPR/Cas9-mediated KO screening approach. Our ChIP data clearly showed that C/EBPδ specifically binds within S100a8 and S100a9 promoter regions. Co-transfection of an inducible C/EBPδ construct and S100a8 and S100a9 reporter constructs not only demonstrated S100 promoter activation due to C/EBPδ expression, but also revealed functional relevance of specific binding sites, via promoter bashing, that are located exactly within the stated promoter regions. The two DNA motifs for specific C/EBPδ responses on the S100a8 promoter regions did not share the core sequence 5’-C/G GCAAT-3’ that we found within the S100a9 promoter region in our study. The latter has been described in three other promoters, the human PPARG2 promoter (Lai et al., 2008), the murine and human CEBPD promoter itself (Wang et al., 2021), and the human COX2 promoter (Wang et al., 2006b). We were able to show that the functionally relevant C/EBPδ-binding sites within the S100 promoters lie within genome regions which switch from closed to open chromatin states during monocyte differentiation, and concomitant induction of S100 expression as examined by ATAC-seq.
Our chromatin accessibility data on S100a8 and S100a9 promoter regions reflected active S100a8 and S100a9 transcription on day 3 cells and accurately mirrored decreased S100 expression of C/EBPδ KO monocytes in relation to WT counterparts. These observations were again reflected and supported by the characterization of the epigenetic landscape using H3K27ac and H3K27me3 marks. The fact that H3K27me3 marks were strongly decreased in WT monocytes, but not in precursors or C/EBPδ-deficient monocytes, showed the indispensability of H3K27 demethylation for S100a8 and S100a9 expression. Moreover, our data demonstrated that the Jumonji C family histone demethylase JMJD3 regulates S100a8 and S100a9 expression by erasure of H3K27me3 in dependency of C/EBPδ, which was confirmed by GSK-J4-mediated inhibition of JMJD3 activities. Neither a link of C/EBPδ nor of S100A8/A9 and JMJD3 has been published yet. It has been shown that histone demethylase activities of recombinant JMJD3 on mono-nucleosome substrates is relatively low in contrast to higher activities on bulk histones (Lan et al., 2007), suggesting that further factors, such as C/EBPδ, are involved in chromatin binding. Several studies highlight JMJD3 as a regulator of innate immune responses, especially via NF-κB-mediated inflammation in macrophages (Na et al., 2016; Na et al., 2017; Davis et al., 2020). Accordingly, knockdown of JMJD3 affected mainly inflammatory response networks in monocytic THP-1 cells (Das et al., 2012) and blocked activation of the NLRP3 inflammasome in bone marrow-derived macrophages (Huang et al., 2020). GSK-J4 treatment of mice attenuated disease progression and inflammatory activities in several mouse models for inflammatory diseases, such as arthritis (Jia et al., 2018), colitis (Huang et al., 2020), and EAE (experimental autoimmune encephalomyelitis) (Doñas et al., 2016). Accordingly, GSK-J4 treatment of our ER-Hoxb8 monocytes reduced expression of the proinflammatory alarmins S100a8 and S100a9, which have been shown to drive the inflammatory process of arthritis (van Lent et al., 2012). With our study, we have taken a step forward to uncover the role of epigenetic features on S100a8 and S100a9 expression and, thereby, on inflammatory conditions in murine monocytes.
We were also able to demonstrate an association of CEBPD, and S100A8 and S100A9 expression in the context of human cardiovascular disease. The expression of these molecules shows a significant positive correlation not only to each other but also to the manifestation of sCAD and MI in patient-derived PBMCs. Moreover, expression of CEBPD, S100A8 and S100A9 showed an even stronger association with classical, proinflammatory monocytes (CD14++CD16-), compared to non-classical (CD14+CD16++) and intermediate (CD14++CD16+) monocytes. The endogenous antagonists of C/EBPδ, ATF3 and especially FBXW7, showed a negative correlation of their expression pattern in these monocyte subpopulations. Interestingly, inflammatory monocytes with phagocytic and proteolytic activities have been reported to show an early peak at infarct sites, which are followed by infiltration of non-classical, anti-inflammatory monocytes (Nahrendorf et al., 2007; Dutta and Nahrendorf, 2015). Genetic deletion of S100A8/A9 was reported to attenuate MI and improve cardiac function in murine models. In contrast, overexpression of S100a9 in mice increased infarct size and mortality, and treatment with recombinant S100 proteins raised influx of immune cells into the infarct area (Li et al., 2019; Sreejit et al., 2020). Moreover, serum concentrations of S100A8/A9 are known to be highly sensitive and prognostic markers for myocardial injury (Aydin et al., 2019). Taken together, these data indicate that the C/EBPδ-S100 alarmin axis drives a clinically relevant pathomechanism in cardiovascular disease and probably other inflammation-driven conditions.
There are several published reports suggesting a biomedical relevance of the link between the C/EBPδ and the S100A8/A9 alarmin under other inflammatory conditions as well. For example, C/EBPδ has been shown to play a role in the pathogenesis of psoriasis (Lan et al., 2020) and in acute inflammatory signaling by regulating COX-2 (Wadleigh et al., 2000), IL-6 (Litvak et al., 2009), and TLR4 (Balamurugan et al., 2013). Analysis of the genome-wide transcription pattern of monocytes revealed IL6 as the top gene induced by S100 alarmin stimulation via interaction with TLR4 (Fassl et al., 2015), and targeted deletion of S100A9 ameliorated inflammation in a murine psoriasis model (Zenz et al., 2008). Additionally, C/EBPδ levels were elevated in mouse models and patients of Alzheimer’s disease (AD) (Li et al., 2004; Ko et al., 2012) and RA (Nishioka et al., 2000; Chang et al., 2012). In mouse models of AD, downregulation (Ha et al., 2010) and deficiency of S100A9 (Kummer et al., 2012) had therapeutic effects on disease activity. Also in human studies, S100A9 was found to be associated with AD pathogenesis (Shepherd et al., 2006). Beyond that, S100A8 and S100A9 are known key players in the pathogenesis of arthritis in murine models (van Lent et al., 2012). Gene expression profiling of blood cells from RA patients receiving anti-TNF-α-based treatment showed that both CEBPD and S100A8 were downregulated by the treatment (Meugnier et al., 2011). Uncontrolled activity of S100A8/A9 alarmins drives TNF-induced arthritis in mice (Vogl et al., 2018). In the context of human RA, the expression and serum concentrations of S100A8/A9 correlate very well with disease activity and are the first predictive marker for disease relapses in juvenile patients, and of the responses to therapy in juvenile and adult patients (Moncrieffe et al., 2013; Choi et al., 2015). However, no direct molecular or functional link between S100A8/A9 and C/EBPδ in arthritis has yet been reported. Using an in vivo model for acute lung inflammation, we were now able to show that C/EBPδ-deficient mice express lower S100A8/A9 level than WT littermates in response to LPS exposure, which is accompanied by a milder disease phenotype. In S100A9 KO mice, neutrophil recruitment in the lung was impeded when using the same mouse model (Chakraborty et al., 2017), suggesting that S100A8/A9 expression mediates effects of C/EBPδ as a key mediator of LPS-induced lung inflammation (Yan et al., 2013). Our data define a new regulatory axis between the transcription factor C/EBPδ and the alarmins S100A8 and S100A9 in myeloid cells which is of relevance in different inflammatory processes and clinical diseases.
Materials and methods
Murine model of acute lung inflammation
Request a detailed protocolWT and C/EBPδ KO (kindly provided by Esta Sterneck, National Cancer Institute, Frederick, MD) (Sterneck et al., 1998) mice were exposed to LPS from Salmonella enterica (0.5 mg/ml, Sigma-Aldrich) in saline or saline (NaCl) only for 45 min inside a cylindrical Pyrex chamber which was connected to a nebulizer. After a 4 hr resting period, mice were anesthetized, blood was collected from heart and BALF was collected for S100A8/A9 measurements by flushing the lungs five times with 0.8 ml NaCl through the trachea. Serum and BALF were saved for later analyses of alarmins and cytokines at –80°C. Migrated polymorphonuclear leukocytes were analysed via Kimura and following Ly-6B.2/Gr-1-antibody (Bio-Rad) staining. Mouse experiments were in accordance with German Animal Welfare Legislation and performed as approved by the North Rhine-Westphalia Office of Nature, Environment and Consumer Protection (LANUV), and the District Government and District Veterinary Office Muenster under the reference number 81-02.04.2019.A445.
Cell culture
Request a detailed protocolER-Hoxb8 cells were generated as described earlier (Wang et al., 2006a) and grown in RPMI medium (Thermo Fisher Scientific) supplemented with 10% FBS (Biowest), 1% penicillin/streptomycin solution (Sigma-Aldrich), 1% glutamine solution (Thermo Fisher Scientific), 40 ng/ml recombinant mouse GM-CSF (ImmunoTools) an 1 µM β-estradiol (Sigma-Aldrich). For differentiation, precursor cells were washed and incubated in estradiol-free medium containing 40 ng/ml GM-CSF for several days. HEK293T were grown in DMEM (Thermo Fisher Scientific) supplemented with 10% FBS (Biowest) and 1% penicillin/streptomycin solution (Sigma-Aldrich), 1% glutamine solution (Thermo Fisher Scientific), and 1% sodium pyruvate (Merck). All cell lines were cultured at 37°C, 5% CO2, and routinely screened and found negative for mycoplasma contamination in a PCR-based assay (PromoCell).
Cell line generation and manipulation
Request a detailed protocolWT, C/EBPδ KO (kindly provided by Esta Sterneck, National Cancer Institute, Frederick, MD) (Sterneck et al., 1998) and Cas9-expressing (Chiou et al., 2015) ER-Hoxb8 cells originated from corresponding mice. PHF8, CSRP1, HAND1, FBXW7, ATF3, STAT3, KLF5, IRF7, and C/EBPβ KO ER-Hoxb8 cells were generated using CRISPR/Cas9 as described earlier (Gran et al., 2018). The oligos for gRNA cloning are listed in Supplementary file 1a. For lentiviral production, the lentiGuide-Puro (for GeCKO screen), lentiCRISPRv2-gRNA (for single KO cell lines), or TRE_3xFlag-C/EBPδ was co-transfected into HEK293T cells, together with the packaging plasmids pCMV-VSV-G (AddGene, #8454) and psPAX2 (AddGene, #12260). For transduction of ER-Hoxb8 cells, cells were incubated with lentiviral particles and 8 µg/ml polybrene (Sigma-Aldrich) for 1 hr upon spinfection and selected for several days using puromycin (InvivoGen). For transfection of HEK293T cells, the cells were seeded 1 day prior to transfection. Then, cells were co-transfected with TRE_3xFlag-C/EBPδ and S100a8 reporter or S100a9 reporter using the Lipofectamine 3000 Transfection Reagent (Thermo Scientific) according to the manufacturer’s manual. For inhibition of JMJD3 activities, cells were treated using 5 µM GSK-J4 HCl (SellekChem) for 3 days. To induce cebpd in TRE_3xFlag-C/EBPδ ER-Hoxb8 cells or transfected HEK293T cells, cells were treated using 2 µg/ml doxycycline (Sigma-Aldrich) for 24 hr.
Isolation of BMDMs
Request a detailed protocolBMDMs were obtained by flushing the femurs from WT and C/EBPδ KO mice. Erythrocytes were depleted by osmotic shock. Cells were washed and separated using a Ficoll gradient (PAN-Biotech). Primary monocytes (M0) were generated by culturing for 3 days in DMEM containing with 10% FBS (Biowest) and 1% penicillin/streptomycin solution (Sigma-Aldrich), 1% glutamine solution (Thermo Fisher Scientific), and 15% of L929 supernatant as a source of macrophage colony-stimulating factor. Stimulation of cells with 50 ng/ml IFN-γ (ImmunoTools) and 10 ng/ml LPS (Sigma-Aldrich) for 24 hr was used for M1-polarization, whereas a 24 hr stimulation with 20 ng/ml IL-4 (Peprotech) led to M2-polarization of BMDMs.
GeCKO-library screening
Request a detailed protocolAmplification of mouse CRISPR Knockout pooled library (GeCKO v2) in lentiGuide-Puro plasmid, purchased from AddGene (#1000000053) (Sanjana et al., 2014), was performed as described (Joung et al., 2017). Cas9-expressing ER-Hoxb8 cells, transduced with library lentiviral particles at an MOI of 0.4, were differentiated to day 3. Intracellular S100A9 was stained with an S100A9-FITC coupled antibody using the Foxp3/Transcription Factor Staining Buffer Set (eBioscience). Cells with no/lower S100A9 expression (hits) and cells with normal S100A9 expression (reference) were sorted using an SH800S Cell Sorter (Sony, Minato, Japan) and DNA was purified by phenol-chloroform extraction. Next-generation sequencing was performed as described earlier (Joung et al., 2017). Briefly, sgRNA library for next-generation sequencing was prepared via PCR using primers amplifying the target region with Illumina adapter sequences (Supplementary file 1b), the purified DNA, and the NEBNext High-Fidelity 2× PCR Master Mix (NEB). PCRs were pooled and purified using the QIAquick PCR Purification Kit (Qiagen). Size and quantity was determined using the Bioanalyzer High Sensitivity DNA Analysis Agilent High Sensitivity DNA Kit (Agilent). Samples were sequenced according to the Illumina user manual with 80 cycles of read 1 (forward) using the NextSeq 500/550 High Output Kit v2.5 (75 Cycles) (Illumina) with the 20% PhiX spike in Illumina PhiX control kit (Illumina).
Cloning and plasmid production
TRE_3xFlag-C/EBPδ and TRE_ctrl
Request a detailed protocolThe pcDNA 3.1 (-) mouse C/EBPδ expression vector (AddGene, #12559) and annealed oligonucleotides (Supplementary file 1c) were digested using XbaI and EcoRI and then ligated. Using primers carrying restriction enzyme recognition sites (Supplementary file 1c), the 3xFlag-C/EBPδ expression cassette was amplified. The resulting amplicon and the pCW57.1 mDux-CA target vector (AddGene, #99284) (Whiddon et al., 2017) were digested using NheI and AgeI and subsequently ligated. TRE_ctrl was produced by digesting TRE_3xFlag-C/EBPδ using NheI and AgeI and by subsequent blunting of ends by 3' overhang removal and fill-in of 3' recessed (5' overhang) ends using DNA Polymerase I, Large (Klenow) Fragment (NEB) prior to ligation.
S100a8 and S100a9 reporter
Request a detailed protocolTo construct S100 reporter vectors, 1500 bp upstream of S100a8 and 1800 bp upstream of S100a9 TSS were amplified from genomic mouse DNA. Using primers carrying restriction enzyme recognition sites (Supplementary file 1d), promoter regions were amplified and cloned into pLenti CMV GFP Blast vector (AddGene, #17445) (Campeau et al., 2009) using XbaI and ClaI. Resulting S100 prom-GFP constructs were cloned into MSCV-PIG-empty vector (AddGene, #105594) (Xu et al., 2018) by digestion with NsiI and ClaI together with the MSCV backbone to exchange IRES-GFP-cassette with S100a8/a9prom-GFP-cassette and subsequent ligation. Proposed C/EBP DNA-binding sites within S100a8 and S100a9 promoter regions were identified using the AliBaba2.1 net-based transcription factor-binding site search tool (Grabe, 2002), and were mutated by deleting 6–7 base pairs using the QuikChange II XL Site-Directed Mutagenesis Kit (Agilent Technologies). The primers used for mutagenesis are listed in Supplementary file 1e. Plasmids were produced in DH5α cells and purified using the PureLinkTM HiPure Plasmid Midiprep Kit (Thermo Scientific).
qRT-PCR
Request a detailed protocolRNA was isolated using a NucleoSpin Extract II Isolation Kit (Macherey Nagel). The mRNA expression of selected genes was measured by qRT-PCR as described earlier (Heming et al., 2018). The primers used are listed in Supplementary file 1f. The relative expression level of each target gene was analysed using the 2−ΔΔCq method and was normalized to GAPDH.
Chromatin immunoprecipitation
Request a detailed protocolFor chromatin preparation, progenitor and differentiated ER-Hoxb8 cells were fixed using 1% formaldehyde for 5 min and reaction was stopped by adding 125 mM glycine. Chromatin was extracted as previously described (Fujita and Fujii, 2013). Approximately 1–5% of chromatin served as the input sample. DNA from input samples was isolated using phenol-chloroform extraction as described earlier (Heming et al., 2018). For immunoprecipitation, 3 µg antibody against Flag (Sigma-Aldrich), H3K27ac (Abcam), H3K27me3 (Cell Signaling Technology), normal Rabbit IgG (Cell Signaling,) or Mouse IgG1, κ Isotype control (BioLegend) was conjugated to 900 µg magnetic Dynabeads-Protein G (Thermo Scientific) at 4°C overnight. Sonicated chromatin was added to AB-conjugated Dynabeads and incubated at 4°C overnight. The Dynabeads were washed as described earlier (Fujita and Fujii, 2013). For elution, Dynabeads were incubated twice with elution buffer (0.05 M NaHCO3, 1% SDS) at 65°C for 15 min. DNA from eluates was isolated using phenol-chloroform extraction as with input samples. Values were taken into account only when the amount of DNA pulled down by using the antibody of interest was more than 5-fold increased over DNA pulled down by using IgG antibodies. The primers used for ChIP-PCR are listed in Supplementary file 1g.
ATAC-seq
Request a detailed protocolPrecursor and day 3 differentiated WT and C/EBPδ KO ER-Hoxb8 cells were harvested, washed, and cryopreserved in 50% FBS/40% growth media/10% DMSO using a freezing container at –80°C overnight. Cells were shipped to Active Motif to perform ATAC-seq as previously described (Buenrostro et al., 2013).
Measurements of S100A8/A9 protein level
Request a detailed protocolThe S100A8/A9 protein concentrations were measured using an in-house S100A8/A9 enzyme-linked immunosorbent assay (ELISA), as previously described (Vogl et al., 2014).
Cytokine measurements by bead-based immunoassay
Request a detailed protocolIL-1α, IFN-γ, GM-CSF, MCP-1, IL-12p70, IL-1β, IL-10, IL-6, IL-27, IL-17A, IFN-β, and TNF-α were analysed using the bead-based immunoassay LEGENDplex mouse inflammation panel according to manufacturer’s instructions (BioLegend). Analytes were measured by flow cytometry using Navios (Beckmann Coulter).
Immunoblotting
Request a detailed protocolCells were lysed in M-PER Mammalian Protein Extraction Reagent (Thermo Scientific) containing a protease inhibitor mixture (Sigma-Aldrich). Protein concentration was determined, and equal amounts (15–30 µg) were run on an SDS-PAGE. After blotting on a nitrocellulose membrane, the membrane was incubated overnight with primary antibodies against: polyclonal rabbit S100A8 and S100A9 antibodies (originated from our own production; Vogl et al., 2014), GAPDH (Cell Signaling Technology), α/β-Tubulin (Cell Signaling Technology), and Flag (Sigma-Aldrich). Membranes were incubated with an HRP-linked secondary antibody (Agilent) for 1 hr. Chemiluminescence signal was detected using ChemiDoc XRS+ (Bio-Rad) together with ImageJ (National Institutes of Health) to quantify the signal intensity.
Phagocytosis
Request a detailed protocolFluoSpheres polystrene microspheres (Thermo Scientific) that were shortly sonicated in a bath sonicator (Latex Beads) or pHrodo Red S. aureus Bioparticles Conjugate (S. aureus Bioparticles) were added to 5 × 105 differentiated cells at a ratio 1:10 for 2 hr at 37°C. The rate of phagocytosis was measured by flow cytometry using Navios (Beckmann Coulter).
Oxidative burst
Request a detailed protocolCells were stimulated using 10 nM PMA (Abcam) for 15 min or left untreated. After incubation, 15 μM DHR123 (Sigma-Aldrich) were added for another 15 min. The fluorescence signal was analysed using flow cytometry (Navios, Beckmann Coulter).
FACS analysis
Request a detailed protocolBMCs and ER-Hoxb8 cells were differentiated using 40 ng/ml rmGM-CSF (ImmunoTools), harvested, stained for CD11b, Ly-6C, and appropriate Isotype control antibodies (all BioLegend) and measured using flow cytometry (Navios, Beckmann Coulter) to determine cell population purity and differentiation kinetics.
RNA-seq
Study population
Request a detailed protocolFor this study, we used bulk mRNA-sequencing (RNA-seq) data of PBMCs and monocytes from two subsets of participants in the German BioNRW Study (Witten et al., 2022). BioNRW actively recruits patients undergoing coronary angiography for the diagnosis and percutaneous coronary intervention of coronary artery disease, as well as age and gender matched healthy control individuals without history of cardiovascular disease, all aged 18–70 years of age. Patients receive standard cardiovascular care and medication (ACE-inhibitor, AT1-receptor blocker, β-blocker, diuretics, statin), according to current guidelines. Here, we included a total of 42 patients with sCAD or acute MI, as well as 39 of the corresponding age and sex matched controls.
The BioNRW Study is conducted in accordance with the guidelines of the Declaration of Helsinki. The research protocol, including the case report forms, was approved by the local ethics committee (#245-12). Written informed consent was obtained from all study participants.
Blood collection and isolation of PBMCs
Request a detailed protocolIn case of MI, blood samples were collected during the first 4 days following the event. EDTA blood was drawn from each subject by venipuncture. Sample processing followed within 2 hr. PBMCs were obtained from 40 ml blood by density gradient centrifugation (Ficoll; Biochrom). Lymphocytes were collected and washed twice with PBS. The pellet was re-suspended in freezing medium Cryo-SFM (PromoCell) and cryopreserved.
Isolation of monocyte subpopulations
Request a detailed protocolAfter washing, PBMCs were stained with anti-human antibodies specific for CD2 (PE, RPA-2.10, T-cell marker), CD14 (APC, M5E2, monocyte subset differentiation), CD15 (PE, HIM1, granulocyte marker), CD16 (PE-Cy7, 3G8, monocyte subset differentiation), CD19 (PE, HIB19, B-cell marker), CD56 (PE, MY31, NK-cell marker), CD335 (PE, 9E2, NK-cell marker), HLA-DR (FITC, TU36, antigen-presenting cells) (all from BD Biosciences), as reported by Cros et al., 2010. Cells were acquired on a FACS LSR II flow cytometer (BD Biosciences) and analysed using FlowJo software version 10 (Treestar Inc). For sorting of monocyte subsets, PBMCs were stained and sorted on a MoFlo Astrios cell sorter (Beckman Coulter). Cells were sorted in 1 ml of Isol-RNA lysis reagent (5-Prime GmbH) and frozen at –80°C. To avoid gender-specific effects, three representative male samples of each BioNRW diagnostic group (sCAD, MI, and control) were selected to be subjected to cell sorting and subsequent RNA isolation.
Differential expression analysis in PBMCs and monocytes
Request a detailed protocolFor mRNA profiling of PBMCs and monocyte subpopulations using RNA-seq, mRNA was enriched using the NEBNext Poly(A) Magnetic Isolation Module (NEB), followed by cDNA NGS library preparation (NEBNext Ultra RNA Library Prep Kit for Illumina, NEB). The size of the resulting libraries was controlled by the use of a Bioanalyzer High Sensitivity DNA Kit (Agilent Technologies) and quantified using the KAPA Library Quantification Kit for Illumina (Roche). Equimolar, appropriately pooled libraries were sequenced in a single read mode (75 cycles) on a NextSeq500 System (Illumina) using v2 chemistry, yielding in an average QScore distribution of 92%≥Q30 score. They were subsequently demultiplexed and converted to FASTQ files using bcl2fastq v2.20 Conversion software (Illumina). Data was quality controlled using FASTQC software and trimmed for adapter sequences using Trimmomatic (Bolger et al., 2014).
General statistics
Request a detailed protocolThe statistical significance of the data was determined using Prism 5.0 software (GraphPad Software, San Diego, CA). Analyses between two groups were performed using an unpaired two-tailed Student’s t test. Comparisons among three or more groups were performed by using one-way ANOVA, followed by Bonferroni’s multiple means tests for comparing all pairs of columns. Differences were considered statistically significant at a probability (p-value) of <0.05.
Bioinformatics analysis
GeCKO-library screening
Request a detailed protocolAnalysis of counting the reads for each gRNA and differential analysis was performed using MaGeCK 0.5.9.3, a computational tool to identify important genes from GeCKO-based screens (Li et al., 2014). A modified RRA method with a redefined ρ value was used. Former RRA computed a significant p-value for genes in the middle of gRNA ranked list and thereby introducing false positives because the assumption of uniformity is not necessarily satisfied in real applications. Thus, top ranked % gRNAs were selected if their negative binomial p-values were smaller than a threshold, such as 0.05. If j of the n gRNAs targeting a gene were selected, then the modified value is defined as ρ=min (p1, p2,…,pj), where j≤n. This modified RRA method could efficiently remove the effect of insignificant gRNAs in the assessment of gene significance. A permutation test where the gRNAs were randomly assigned to genes was performed to compute a p-value based on the ρ values. By default, 100× ng permutations are performed, where ng is the number of genes. We then compute the FDR from the empirical permutation p-values using the Benjamini-Hochberg procedure.
ATAC-seq
Request a detailed protocolSequence analysis was performed by mapping the paired-end 42 bp sequencing reads (PE42) generated by Illumina sequencing (using NextSeq500) to the genome using the BWA algorithm with default settings (‘bwa mem’). Only reads that passed Illumina’s purity filter, aligned with no more than two mismatches, and mapped uniquely to the genome were used in the subsequent analysis. In addition, duplicate reads (‘PCR duplicates’) were removed. BAM files provided by Active Motif were used to perform peak calling with MACS2, using paired-end mode with a bandwidth of 200 and q-value cutoff of 0.01. BedGraph files were converted to BigWig format for visualization with the Integrative Genomics Viewer (IGV) (Thorvaldsdóttir et al., 2013) using bedGraphToBigWig. Differential accessibility analyses for all possible comparisons between WT and/or C/EBPδ KO in day 0 and day 3 groups were carried out with GUAVA v1 (Divate and Cheung, 2018), which implements DESeq2 (Love et al., 2014). Differential peaks were defined as those merged intervals within a window of 5000 bp upstream and 3000 bp downstream of the transcription start site (TSS) showing an adjusted p-value (padj) <0.05 and log2 fold change ≥±1.5.
RNA-seq
Request a detailed protocolThe resulting reads were mapped to the human reference genome builds hg19 (monocytes) or hg38 (PBMCs) using Tophat2 (Kim et al., 2013) or HISAT2 v2.1.0 (Kim et al., 2019), counted by using the R package GenomicAlignments (Lawrence et al., 2013) or HTSeq v0.11.2 (Anders et al., 2015), and followed by differential expression analysis using DEseq2 (Love et al., 2014). The PBMCs dataset used for analysis consisted of 72 individuals, from which 36 were sCAD/MI cases and 36 were controls (21 females and 15 males in each group, mean age: 50.8±12.3 years), while the monocytes dataset contained read counts of classical, intermediate, and non-classical monocyte subpopulations from 9 male individuals (3 MI, 3 sCAD, and 3 controls). One sCAD non-classical monocyte sample had to be excluded from analysis due to low mapping rate; therefore, the monocytes dataset used for analysis contained 26 samples. Genes were considered differentially expressed at adjusted p<0.05 (Benjamini-Hochberg method). R was used to perform Pearson correlation tests and generate plots for the genes of interest from the normalized count data.
Conclusion
Request a detailed protocolWe found that the transcription factor C/EBPδ drives expression of the abundant alarmins S100A8 and S100A9, and demonstrated that C/EBPδ binding to specific sites on S100a8 and S100a9 promoter regions also induced changes in chromatin accessibility via JMJD3-mediated demethylation of H3K27me3 marks, which includes a so far unknown link. Due to the high relevance of S100A8/A9 alarmin expression in many inflammatory diseases, our findings may point to novel molecular targets for innovative anti-inflammatory therapeutic approaches.
Appendix 1
Data availability
GeCKO-screen data are deposited in ID PRJNA754262 at the NCBI’s Database, RNA-seq data are deposited in ID 706411 at the NCBI SRA repository and ATAC-seq data are deposited in ID GSE200730 at the NCBI’s GEO data repository.
-
NCBI BioProjectID PRJNA706411. monocyte subpopulation profiling.
-
NCBI Gene Expression OmnibusID GSE200730. ATAC-seq in precursor and differentiated ER-Hoxb8 cells.
References
-
HTSeq--a Python framework to work with high-throughput sequencing dataBioinformatics (Oxford, England) 31:166–169.https://doi.org/10.1093/bioinformatics/btu638
-
S100 proteins in rheumatic diseasesNature Reviews. Rheumatology 14:528–541.https://doi.org/10.1038/s41584-018-0058-9
-
Biomarkers in acute myocardial infarction: current perspectivesVascular Health and Risk Management 15:1–10.https://doi.org/10.2147/VHRM.S166157
-
Trimmomatic: a flexible trimmer for Illumina sequence dataBioinformatics (Oxford, England) 30:2114–2120.https://doi.org/10.1093/bioinformatics/btu170
-
MRP8/14 serum levels as a strong predictor of response to biological treatments in patients with rheumatoid arthritisAnnals of the Rheumatic Diseases 74:499–505.https://doi.org/10.1136/annrheumdis-2013-203923
-
Gene networking and inflammatory pathway analysis in a JMJD3 knockdown human monocytic cell lineCell Biochemistry and Function 30:224–232.https://doi.org/10.1002/cbf.1839
-
Palmitate-TLR4 signaling regulates the histone demethylase, JMJD3, in macrophages and impairs diabetic wound healingEuropean Journal of Immunology 50:1929–1940.https://doi.org/10.1002/eji.202048651
-
Monocytes in myocardial infarctionArteriosclerosis, Thrombosis, and Vascular Biology 35:1066–1070.https://doi.org/10.1161/ATVBAHA.114.304652
-
Transcriptome assessment reveals a dominant role for TLR4 in the activation of human monocytes by the alarmin MRP8Journal of Immunology (Baltimore, Md 194:575–583.https://doi.org/10.4049/jimmunol.1401085
-
Phagocyte-specific calcium-binding S100 proteins as clinical laboratory markers of inflammationClinica Chimica Acta; International Journal of Clinical Chemistry 344:37–51.https://doi.org/10.1016/j.cccn.2004.02.023
-
Proinflammatory S100 proteins in arthritis and autoimmune diseaseArthritis and Rheumatism 50:3762–3771.https://doi.org/10.1002/art.20631
-
Efficient isolation of specific genomic regions and identification of associated proteins by engineered DNA-binding molecule-mediated chromatin immunoprecipitation (enChIP) using CRISPRBiochemical and Biophysical Research Communications 439:132–136.https://doi.org/10.1016/j.bbrc.2013.08.013
-
Renal collecting duct epithelial cells regulate inflammation in tubulointerstitial damage in miceThe Journal of Clinical Investigation 121:3425–3441.https://doi.org/10.1172/JCI57582
-
AliBaba2: context specific identification of transcription factor binding sitesIn Silico Biology 2:S1–15.
-
MRP-8 and MRP-14, two abundant Ca(2+)-binding proteins of neutrophils and monocytesJournal of Leukocyte Biology 53:197–204.
-
Nutritional immunity beyond iron: a role for manganese and zincCurrent Opinion in Chemical Biology 14:218–224.https://doi.org/10.1016/j.cbpa.2009.11.008
-
Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotypeNature Biotechnology 37:907–915.https://doi.org/10.1038/s41587-019-0201-4
-
High-throughput genetic screens using CRISPR-Cas9 systemArchives of Pharmacal Research 41:875–884.https://doi.org/10.1007/s12272-018-1029-z
-
HDAC1/HDAC3 modulates PPARG2 transcription through the sumoylated CEBPD in hepatic lipogenesisBiochimica et Biophysica Acta 1783:1803–1814.https://doi.org/10.1016/j.bbamcr.2008.06.008
-
Shikonin inhibits CEBPD downregulation in IL‑17‑treated HaCaT cells and in an imiquimod‑induced psoriasis modelMolecular Medicine Reports 22:2263–2272.https://doi.org/10.3892/mmr.2020.11315
-
Software for computing and annotating genomic rangesPLOS Computational Biology 9:118.https://doi.org/10.1371/journal.pcbi.1003118
-
Calcium-dependent tetramer formation of S100A8 and S100A9 is essential for biological activityJournal of Molecular Biology 359:961–972.https://doi.org/10.1016/j.jmb.2006.04.009
-
Different expression ratio of S100A8/A9 and S100A12 in acute and chronic lung diseasesRespiratory Medicine 102:567–573.https://doi.org/10.1016/j.rmed.2007.11.011
-
The p38-MSK1 signaling cascade influences cytokine production through CREB and C/EBP factors in human neutrophilsJournal of Immunology (Baltimore, Md 191:4299–4307.https://doi.org/10.4049/jimmunol.1301117
-
Histone H3K27 Demethylase JMJD3 in Cooperation with NF-κB Regulates Keratinocyte Wound HealingThe Journal of Investigative Dermatology 136:847–858.https://doi.org/10.1016/j.jid.2015.11.029
-
The healing myocardium sequentially mobilizes two monocyte subsets with divergent and complementary functionsThe Journal of Experimental Medicine 204:3037–3047.https://doi.org/10.1084/jem.20070885
-
Improved vectors and genome-wide libraries for CRISPR screeningNature Methods 11:783–784.https://doi.org/10.1038/nmeth.3047
-
Genome-scale CRISPR-Cas9 knockout screening in human cellsScience (New York, N.Y.) 343:84–87.https://doi.org/10.1126/science.1247005
-
Inflammatory S100A9 and S100A12 proteins in Alzheimer’s diseaseNeurobiology of Aging 27:1554–1563.https://doi.org/10.1016/j.neurobiolaging.2005.09.033
-
Integrative Genomics Viewer (IGV): high-performance genomics data visualization and explorationBriefings in Bioinformatics 14:178–192.https://doi.org/10.1093/bib/bbs017
-
Myeloid-related proteins S100A8/S100A9 regulate joint inflammation and cartilage destruction during antigen-induced arthritisAnnals of the Rheumatic Diseases 67:1750–1758.https://doi.org/10.1136/ard.2007.077800
-
Autoinhibitory regulation of S100A8/S100A9 alarmin activity locally restricts sterile inflammationThe Journal of Clinical Investigation 128:1852–1866.https://doi.org/10.1172/JCI89867
-
Transcriptional activation of the cyclooxygenase-2 gene in endotoxin-treated RAW 264.7 macrophagesThe Journal of Biological Chemistry 275:6259–6266.https://doi.org/10.1074/jbc.275.9.6259
-
A hierarchical and collaborative BRD4/CEBPD partnership governs vascular smooth muscle cell inflammationMolecular Therapy. Methods & Clinical Development 21:54–66.https://doi.org/10.1016/j.omtm.2021.02.021
-
Conservation and innovation in the DUX4-family gene networkNature Genetics 49:935–940.https://doi.org/10.1038/ng.3846
-
CCAAT/enhancer-binding protein δ is a critical mediator of lipopolysaccharide-induced acute lung injuryThe American Journal of Pathology 182:420–430.https://doi.org/10.1016/j.ajpath.2012.10.013
-
Activator protein 1 (Fos/Jun) functions in inflammatory bone and skin diseaseArthritis Research & Therapy 10:201.https://doi.org/10.1186/ar2338
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (CRC 1009 B9)
- Johannes Roth
Deutsche Forschungsgemeinschaft (CRC 1009 Z2)
- Johannes Roth
Deutsche Forschungsgemeinschaft (CRC 1009 B8)
- Thomas Vogl
Deutsche Forschungsgemeinschaft (CRU 342 P3)
- Johannes Roth
Deutsche Forschungsgemeinschaft (RO 1190/14-1)
- Johannes Roth
Deutsche Forschungsgemeinschaft (CRU 342 P5)
- Thomas Vogl
Interdisciplinary Center of Clinical Research at the University of Münster (Ro2/023/19)
- Johannes Roth
Interdisciplinary Center of Clinical Research at the University of Münster (Vo2/011/19)
- Thomas Vogl
EU EFRE Bio NRW programme (005-1007-0006)
- Monika Stoll
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank Ursula Nordhues, Eva Nattkemper, Heike Hater, Heike Berheide, Sina Mersmann, Elvira Barg, and Marianne Jansen-Rust for their excellent technical support, and Esta Sterneck (Center for Cancer Research, National Cancer Institute, Frederick, MD) for providing the C/EBPδ KO mice. This work was supported by grants from the Interdisciplinary Center of Clinical Research at the University of Münster (Ro2/023/19, Vo2/011/19), the German Research Foundation CRC 1009 B8, B9, and Z2, CRU 342 P3 and P5 and RO 1190/14–1 (to J Roth and T Vogl) and by the EU EFRE Bio NRW programme (005-1007-0006) to MS. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics
Human subjects: The BioNRW Study is conducted in accordance with the guidelines of the Declaration of Helsinki. The research protocol, including the case report forms, was approved by the local ethics committee (#245-12). Written informed consent was obtained from all study participants.
Mouse experiments were in accordance with German Animal Welfare Legislation and performed as approved by the North Rhine-Westphalia Office of Nature, Environment and Consumer Protection (LANUV) and the District Government and District Veterinary Office Muenster under the reference number 81-02.04.2019.A445.
Copyright
© 2022, Jauch-Speer et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,416
- views
-
- 272
- downloads
-
- 12
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
Membrane proteins are sorted to the plasma membrane via Golgi-dependent trafficking. However, our recent studies challenged the essentiality of Golgi in the biogenesis of specific transporters. Here, we investigate the trafficking mechanisms of membrane proteins by following the localization of the polarized R-SNARE SynA versus the non-polarized transporter UapA, synchronously co-expressed in wild-type or isogenic genetic backgrounds repressible for conventional cargo secretion. In wild-type, the two cargoes dynamically label distinct secretory compartments, highlighted by the finding that, unlike SynA, UapA does not colocalize with the late-Golgi. In line with early partitioning into distinct secretory carriers, the two cargoes collapse in distinct ER-Exit Sites (ERES) in a sec31ts background. Trafficking via distinct cargo-specific carriers is further supported by showing that repression of proteins essential for conventional cargo secretion does not affect UapA trafficking, while blocking SynA secretion. Overall, this work establishes the existence of distinct, cargo-dependent, trafficking mechanisms, initiating at ERES and being differentially dependent on Golgi and SNARE interactions.
-
- Cell Biology
- Microbiology and Infectious Disease
Zika virus (ZIKV) infection causes significant human disease that, with no approved treatment or vaccine, constitutes a major public health concern. Its life cycle entirely relies on the cytoplasmic fate of the viral RNA genome (vRNA) through a fine-tuned equilibrium between vRNA translation, replication, and packaging into new virions, all within virus-induced replication organelles (vROs). In this study, with an RNA interference (RNAi) mini-screening and subsequent functional characterization, we have identified insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) as a new host dependency factor that regulates vRNA synthesis. In infected cells, IGF2BP2 associates with viral NS5 polymerase and redistributes to the perinuclear viral replication compartment. Combined fluorescence in situ hybridization-based confocal imaging, in vitro binding assays, and immunoprecipitation coupled to RT-qPCR showed that IGF2BP2 directly interacts with ZIKV vRNA 3’ nontranslated region. Using ZIKV sub-genomic replicons and a replication-independent vRO induction system, we demonstrated that IGF2BP2 knockdown impairs de novo vRO biogenesis and, consistently, vRNA synthesis. Finally, the analysis of immunopurified IGF2BP2 complex using quantitative mass spectrometry and RT-qPCR revealed that ZIKV infection alters the protein and RNA interactomes of IGF2BP2. Altogether, our data support that ZIKV hijacks and remodels the IGF2BP2 ribonucleoprotein complex to regulate vRO biogenesis and vRNA neosynthesis.