- Open Access
Gene expression profiling of RIP2-knockdown in HD11 macrophages — elucidation of potential pathways (gene network) when challenged with avian pathogenic E.coli (APEC)
BMC Genomics volume 23, Article number: 341 (2022)
Receptor interacting serine/threonine kinase 2 (RIP2), ubiquitous in many tissue/cell types, is the key regulator of immune and inflammatory responses for many diseases, including avian pathogenic E. coli (APEC), which causes a wide variety of localized or systemic infections. However, the molecular mechanisms by which RIP2 drives its transcriptional program to affect immune and inflammatory response upon APEC infection remains poorly understood.
In this study, RNA-seq and bioinformatics analyses were used to detect gene expression and new direct/indirect RIP2 targets in the treatments of wild type HD11 cells (WT), RIP2 knockdown cells (shRIP2), APEC stimulation cells (APEC), and RIP2 knockdown cells combined with APEC infection (shRIP2 + APEC). The results revealed that a total of 4691 and 2605 differentially expressed genes (DEGs) were screened in shRIP2 + APEC vs. APEC and shRIP2 vs. WT, respectively. Functional annotation analysis showed that apoptosis, MAPK, p53, Toll-like receptor, and Nod-like receptor signaling pathways were involved in APEC-induced RIP2 knockdown HD11 cells. By analyzing the enriched pathway and gene networks, we identified that several DEGs, including HSP90AB1, BID, and CASP9 were targeted by RIP2 upon APEC infection.
As a whole, this study can not only provide data support for constructing gene networks of RIP2 knockdown with APEC challenge but also provide new ideas for improving the immune and inflammatory response.
Avian pathogenic E. coli (APEC), a type of extraintestinal pathogenic E. coli, is the causative agent of avian colibacillosis, which can result in significant economic losses due to the mortality and the reduced productivity of affected birds [1,2,3]. Although great progress has been made for the treatment of APEC, the disease is still hard to eradicate due to the diversity of APEC serotypes causing disease and the emergence of antimicrobial resistance among APEC [4, 5]. The pathogenesis of APEC is still not well understood but appears to influence apoptosis of immune cells and tissues injury [6,7,8]. Moreover, it has been demonstrated that APEC shared the similar genomic sequences with uropathogenic E. coli responsible for human urinary tract infections [9, 10], indicating the horizontal gene transfer ability and zoonotic potential of APEC. Therefore, it is critical to study the mechanism of the immune response and effectively prevent the excessive inflammation at genomic level.
Receptor interacting serine/threonine kinase 2 (RIP2) is the key molecule in the regulation of immune responses, inflammation, and cell death against viral and bacterial infections [11,12,13,14]. The N-terminal CARD of RIP2 can bind to the CARD domain of procaspase-1 to activate NFκB signal cascade and initiate the pro-inflammatory pathways [15,16,17]. Moreover, it has been demonstrated that cells deficient of RIP2 could decrease the activation of NFκB, resulting in impaired expression of IL6, TNFα, IP10 and reduced neutrophil infiltration, which alleviate the excessive inflammatory response [18,19,20]. The research of Homer et al.  discovered that RIP2 tyrosine kinase has a dual function in NOD2-dependent autophagy, that is RIP2 both sends a positive autophagy signal through the activation of p38 MAPK and relieves repression of autophagy mediated by the phosphatase PP2A. Currently, RIP2 has been found to play an important role in the response of the avian host to an APEC infection . Although it is considered the master regulator of immune and inflammatory response, it is not clear how the ubiquitous RIP2 can direct the immune or inflammation-specific process.
In this study we attempted to unravel the gene network underlying RIP2 regulation associated with APEC infection via RNA-seq analysis, using the chicken HD11 macrophage cell line. We studied the expression profiles of complete and knockdown of RIP2 HD11 cells with or without 24 h APEC post-infection. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) mechanistic analyses were performed taking into consideration the function of mRNAs, which changed their expression significantly in RNA-seq analysis. Altogether, our results suggest that HSP90AB1, BID, and CASP9 were the new targets downstream of RIP2, whose depletion directly or indirectly caused the suppression of those genes. Thus, HSP90AB1, BID, and CASP9 were involved in the immune and inflammatory response, while precisely RIP2 could down-regulate those genes during APEC stimulation.
RIP2 inhibition level in chicken HD11 macrophages with or without APEC challenge
To identify the inhibition level of RIP2, the protein and mRNA levels were first measured in wild type HD11 cells (WT), HD11 cells with knockdown of RIP2 (shRIP2), APEC-challenged HD11 cells (APEC), and APEC-challenged HD11 cells with knockdown of RIP2 (shRIP2 + APEC). As shown in Fig. 1A and B, a significantly decreased mRNA and protein expression level of RIP2 were observed in the RIP2 knockdown group in comparison to WT (p < 0.05), indicating the expression of RIP2 was indeed inhibited. Also, APEC challenge resulted in the significantly increased mRNA and protein abundance of RIP2 compared with WT (p < 0.05). However, the mRNA and protein level of RIP2 were found to be significantly decreased in the RIP2 knockdown HD11 cells with APEC challenge in comparison to wild type HD11 cells or APEC challenge group (Fig. 1A-C).
Cell viability in knockdown of RIP2 HD11 macrophages with or without APEC challenge
Morphological changes of APEC-challenged HD11 cells with or without knockdown of RIP2 were observed. As displayed in Fig. 1D, there was no effect of RIP2 knockdown on HD11 cells growth and proliferation. After challenge with APEC at 108 cfu/mL for 24 h, cytopathic effects appeared in HD11 cells when compared with the WT, while knockdown of RIP2 gene could significantly alleviate the APEC-challenged cytopathy. Moreover, the viability of chicken HD11 cells toward APEC challenge with or without knockdown of RIP2 was measured using an CCK-8 assay. There was no difference between WT and shRIP2 group (p > 0.05; Fig. 1E), which was consistent with the results of cells morphology. The APEC-challenged chicken HD11 cells had a significantly lower cell survival rate than those in the RIP2 knockdown group at 24 h post-infection (p < 0.05; Fig. 1E), indicating knockdown of RIP2 can effectively inhibit the APEC-challenged cell apoptosis in HD11 cells.
Nitric oxide (NO) production in knockdown of RIP2 HD11 macrophages with or without APEC challenge
Considering the function of chicken HD11 macrophages, NO production in the cell supernatant from the four different groups (WT, shRIP2, APEC, and shRIP2 + APEC) was determined using the Griess reagent kit. As shown in Fig. 1F, there was no significant difference between WT and RIP2 knockdown group for NO production. However, APEC challenge and RIP2 knockdown combined with APEC challenge groups induced significant NO production in HD11 macrophages, with APEC challenged group giving higher levels of NO compared to RIP2 knockdown combined with APEC challenged group (p < 0.05, Fig. 1F). These results suggested that knockdown of RIP2 had the ability to effectively reduce the APEC-induced NO production, although the level was still significantly higher than WT (p < 0.05).
Overview of RNA sequencing data
A total of twelve cDNA libraries were constructed, respectively, from four groups (WT, shRIP2, APEC, and shRIP2 + APEC) with three biological replicates for each. A total of 71,438,338-95,374,638 sequence reads were obtained, and each sample yielded approximately 80,223,040 clean reads (range from 69,393,540 to 92,952,070, Table S3). The percentages of Q30 for the clean reads was more than 99%, and the average of GC content of clean reads was 52.34% (Table S3). Moreover, 92.81–94.10% of the clean reads were found to successfully map to the chicken reference genome, of which 89.16–91.25% were uniquely mapped to genome (Table S4). Additionally, Table 1 shows the information on the fraction of reads mapping to features for each sample, i.e. any expressed parts of the genome. Interestingly, the percentage of reads mapping to coding sequence (CDS) in WT, shRIP2, and shRIP2 + APEC were 59.78, 59.91, and 60.06%, respectively, while that in APEC was 56.95% (Table 1).
Also, sequencing homogeneity was assessed to evaluate the bias of RNA-seq data. Results showed reads were evenly distributed on the genes without 5′ or 3′ bias (Fig. 2A). Then, a multidimensional scaling analysis was performed on the count data, revealing the biological replicates were relatively concentrated in each of the four groups (Fig. 2B). The WT and shRIP2 HD11 macrophages with or without APEC infection were clearly separated (Fig. 2B). Furthermore, the heatmap of sample correlation showed the shRIP2 + APEC group clustered with the WT, and then the shRIP2, finally the APEC group (Fig. 2D), which were consistent with the results of reads distribution to CDS in Table 1. Then, the dynamic range of the expression values was estimated and exhibited as a box plot of logarithmic transformed RPKM values for each sample separately (Fig. 2C), and the RPKM density distribution is presented in Fig. 2E.
Identification of differentially expressed genes (DEGs)
In order to investigate the specific mechanism that RIP2 affected cell immune and inflammatory response upon APEC infection, DEGs were identified from WT and shRIP2 HD11 cells with or without APEC challenge. Pairwise comparisons were performed as follows: shRIP2 (knockdown of RIP2 HD11 cells) vs. WT (complete HD11 macrophages) and shRIP2 + APEC (knockdown of RIP2 HD11 cells combined with avian pathogenic E. coli infection) vs. APEC (avian pathogenic E. coli infection HD11 cells). A total of 2605 DEGs were obtained when adjusted p-value ≤0.05 and |logFC| ≥ 0.58 were set as the cut-off limits (Fig. 3A, B) in shRIP2 vs. WT. Among them, 1133 down-regulated DEGs and 1472 up-regulated DEGs were identified (Fig. 3C). In shRIP2 + APEC vs. APEC, a total of 4691 DEGs were identified (Fig. 3D, E), of which 2340 DEGs were down-regulated and 2351 DEGs were up-regulated (Fig. 3F).
Functional annotation of DEGs in different comparisons
The gene ontology (GO) classification system was used to classify the possible functions of DEGs in different comparisons. A total of 734 genes (28.18%) and 1476 (31.46%) were successfully assigned to at least one GO term annotation in shRIP2 vs. WT (Table S5) and shRIP2 + APEC vs. APEC (Table S6), respectively. The significantly enriched GO terms in shRIP2 vs. WT and shRIP2 + APEC vs. APEC were similar (Fig. 4). For the molecular function category, the top two largest categories were “binding” and “catalytic activity”. According to biological process, the top categories were “cellular process”, “biological regulation”, and “metabolic process”. More remarkable, the immune related GO terms were also identified, including “cell communication”, “immune response”, “signal transduction”, and “response to stimulus”. These results show that knockdown of RIP2 significantly affected immune system processes in chicken HD11 macrophages with or without APEC challenge.
Then, Kyoto Encyclopedia of Genes and Genomes (KEGG) classification system was also performed to identify the possible functions of DEGs. A total of 15 and 55 significantly changed pathways were detected in shRIP2 vs. WT (Table S7) and shRIP2 + APEC vs. APEC (Table S8), respectively, with an adjusted p-value ≤0.05. The significantly enriched pathways were primarily involved in immune system for the two comparisons, which included “Phagosome”, “Lysosome”, “MAPK signaling pathway” (Fig. 5). In addition, the pathways related to cell growth, differentiation, survival, signal transduction were also significantly enriched in the two comparisons, including “Apoptosis”, “ECM-receptor interaction”, “VEGF signaling pathway”, and “Focal adhesion” (Fig. 5). However, the “p53 signaling pathway”, “cell cycle”, “Toll-like receptor signaling pathway”, and “NOD-like receptor signaling pathway” were uniquely identified to be significantly enriched in shRIP2 + APEC vs. APEC (Fig. 5B).
RIP2 is a critical regulator for MAPK signaling and apoptosis pathway in chicken HD11 macrophages with or without APEC challenge
After knockdown of RIP2, the MAPK signaling and apoptosis pathway were significantly activated with or without APEC challenge. A total of 39 and 67 DEGs were involved in “MAPK signaling pathway” in shRIP2 vs. WT and shRIP2 + APEC vs. APEC, respectively, as shown in Tables S9 and S10. Nineteen DEGs were commonly involved in the MAPK signaling pathway in both shRIP2 vs. WT and shRIP2 + APEC vs. APEC. To further corroborate the correlation between RIP2 and MAPK signaling pathway, we used qRT-PCR to detect the expression of the commonly identified DEGs—IL1β, JUN, and CASP3 in the two comparisons, as well as the uniquely expressed DEGs—MAPK9 and FOS in shRIP2 + APEC vs. APEC. The qRT-PCR results were in agreement with their transcript abundance changes determined by RNA-seq (Fig. 6).
Moreover, a total of 25 and 48 DEGs were enriched in “apoptosis pathway” in shRIP2 vs. WT (Table S11) and shRIP2 + APEC vs. APEC (Table S12), respectively. Twelve DEGs were enriched in apoptosis pathway in the two comparisons, 3 (CTSB, CTSC, and PIK3CB) of which were selected and used for qRT-PCR validation. Moreover, the uniquely enriched DEGs— FADD, CFLAR, BID, and CASP9 in the comparison of shRIP2 + APEC vs. APEC were also identified for RNA-seq data reliability. The results showed that knockdown of RIP2 indeed affected the expression of the the selected genes with or without APEC challenge, which presented 100% consistency between RNA-seq and qRT-PCR (Fig. 6). And, more remarkable, three DEGs (JUN, CASP3, and ENSGALG00000031518) were involved in both MAPK signaling and apoptosis pathway in the two comparisons (Fig. S1), suggesting they had important role in APEC challenge and related with RIP2. All these results show that RIP2 is a critical regulator for MAPK signaling and apoptosis pathway with or without APEC infection.
Knockdown of RIP2 impacts the activation of p53 and those genes involved in TLR and NLR signaling pathway upon APEC challenge
Also, we identified the uniquely significantly enriched pathways “p53 signaling pathway”, “Toll-like receptor signaling pathway”, and “NOD-like receptor signaling pathway” in the knockdown of RIP2 HD11 cells toward APEC challenge (shRIP2 + APEC vs. APEC). A total of 26 DEGs were involved in p53 signaling pathway, of which consisted 16 up-regulated DEGs (FC = 1.5 ~ 6.76) and 9 down-regulated DEGs (FC = − 7.18 ~ − 1.5) (Fig. 7A). Additionally, 16 DEGs (10 up-regulated DEGs with FC of 1.5 ~ 4.38 and 6 down-regulated DEGs with FC of − 5.83 ~ − 1.5) were identified to be involved in NLR signaling pathway, which were the downstream genes of RIP2 (Fig. 7B). Moreover, there were 30 DEGs enriched in TLR signaling pathway, including 15 up-regulated DEGs (FC = 1.5 ~ 4.38) and 15 down-regulated DEGs (FC = − 15.95 ~ − 1.5) (Fig. 8). Based on the analyses aforementioned, 8 DEGs were randomly selected for validation by qRT-PCR, which includes CD80, CD86, IL18, CCL5, MyD88, HSP90AB1, TLR7, and IL8. The results showed the significant differences in expression levels of those selected genes determined by qRT-PCR analysis are concordant with the RNA-seq results (Fig. 6), indicating the transcriptome sequencing data are reliable.
New target genes of RIP2 upon APEC challenge
In order to find the new target of RIP2 toward APEC infection, we collected and analyzed the DEGs enriched in the uniquely significantly changed pathways “p53 signaling pathway”, “Toll-like receptor signaling pathway”, and “NOD-like receptor signaling pathway” in the comparison of shRIP2 + APEC vs. APEC (Fig. 9). Compared to APEC challenge, RNA-seq data showed that BID, CASP9, and HSP90AB1 were significantly down-regulated in RIP2 knockdown HD11 cells toward APEC infection. To evaluate the RIP2 role in those genes, we studied the expression of above genes in RIP2 knockdown HD11 cell following APEC challenge. As shown in Fig. 10, the mRNA and protein levels of HSP90AB1 showed remarkable reduction in RIP2 knockdown HD11 cells (p < 0.001) combined APEC challenge compared to the APEC challenge group. Interestingly, compared to the wild type HD11 cells, the mRNA and protein expression level of HSP90AB1 were significantly decreased in the knockdown of RIP2 HD11 macrophages group (Fig. 10C, D, and G), which is consistent with the results of RNA-seq. These results indicate a strong relationship exists between RIP2 and HSP90AB1 with or without APEC challenge. However, interfering RIP2 did not affect the mRNA and protein expression level of BID and CASP9 (Fig. 10A, B, D, E, and F), which is in agreement with the RNA-seq data. APEC infection had the ability to induce a significant difference in the mRNA and protein expression level of BID and CASP9 between RIP2 knockdown HD11 cells and wild type HD11 cells, indicating RIP2 is involved in the regulation of BID and CASP9 in a certain relationship as a result of APEC infection. Taken together, these results showed that the RNA-seq results were considered to be reliable and RIP2 was involved in the regulation of BID, CASP9, and HSP90AB1 in response to APEC infection.
Although RIP2 has been cloned in chicken , the function of RIP2 was only focused on the activation of NFκB signal pathway . This study first demonstrated the essential genes or pathways, as well as the downstream targets regulated by RIP2 with or without APEC challenge through transcriptome analysis. We compared the transcriptomes of APEC-induced chicken HD11 cells with knockdown or complete expression of RIP2. It was found that in RIP2 knockdown HD11 cells, many identified DEGs were grouped into the unique functional categories, p53 (p = 1.67 × 10− 3), TLRs (p = 4.14 × 10− 3), and NLRs (p = 2.54 × 10− 2) signaling, toward APEC infection, as expected in cases of impairment of the key regulator for HD11 cells transcriptional program.
p53 signaling pathway is involved in the regulation of a serials normal life activities of cells, such as apoptosis and senescence [24,25,26]. Normally, p53 has the ability to prevent DNA replication in order to provide time for DNA repair, otherwise it triggers apoptosis, suggesting that p53 plays an important role in monitoring and identifying DNA damage points . Earlier studies have shown that the p53 signaling pathway was significantly suppressed in the comparison of susceptible vs. non-challenged birds at 5 day APEC post-infection in bone marrow and thymus [6, 7], indicating APEC might tend to promote host tumorigenesis. Also, compared to the control group, this pathway was significantly down-regulated in the APEC-induced HD11 cells in the current study, which was consistent with previous studies. However, p53 signaling pathway was significantly up-regulated in the knockdown of RIP2 HD11 cells combined with APEC infection in comparison to APEC infection in the present study (Table S13). Then, we focused on the p53 related DEGs, since the identification of these genes could help to elucidate the molecular mechanism that underlies the function of RIP2 in response to APEC infection. Take SIAH1 as an example, it was up-regulated (Fold changes = 2.3) in shRIP2 + APEC vs. APEC, while down-regulated (Fold changes = − 1.96) in APEC vs. WT. This result was expected, since SIAH1 was a tumor suppressor and played essential role in regulating cell apoptosis . Another two genes of interest that changed in shRIP2 + APEC vs. APEC were PMAIP1 and MDM4. It has been demonstrated that PMAIP1 belonged to the pro-apoptotic subfamily and could determine whether a cell commits to apoptosis . Also, this study demonstrated that MDM4 had the ability to promote cell apoptosis upon genotoxic stress . These two genes were up-regulated in shRIP2 + APEC vs. APEC, indicating knockdown of RIP2 triggered the activation of p53 signaling pathway and further promote the apoptosis upon APEC challenge. These results coincide with the expectations, since RIP2 activity correlated with tumor and metastasis et al. [31,32,33].
Another interesting gene HSP90AB1 is down-regulated in the comparisons of shRIP2 vs. WT and shRIP2 + APEC vs. APEC. We validated these results finding decreased expression levels of both mRNA and protein in RIP2 knockdown HD11 cells in any conditions compared to wild type HD11 cells (Fig. 10C, D, and G), suggesting the existence of a positive and cross-regulation between RIP2 and HSP90AB1 in chicken HD11 cells with or without APEC challenge. Then, bioinformatics analyses showed that HSP90AB1 was involved to the GO group of “response to stimulus” and, most importantly, to the KEGG group of “NLR signaling”, closely related to the immune activity of RIP2. HSP90AB1 is a member of the large family of HSPs and function as molecular chaperones, which is important to signal transduction, protein folding, apoptosis, inflammation, and cell survival [34, 35]. It has been demonstrated that the expression of HSP90AB1 mRNA was often expressed on the tissues of heart, liver, brain, and spleen, regulating the heat stress in chickens . In fact, HSP90AB1 was involved in antiapoptotic and immunostimulatory effects of CpG in both mouse and chicken macrophages [37, 38]. Our transcriptome data showed that the apoptosis and immune related pathways were significantly changed in RIP2 knockdown HD11 cells with or without APEC challenge, which included the gene of interest HSP90AB1. Herein, this study strongly suggests that HSP90AB1 is a downstream target of RIP2 and positively modulated by RIP2.
Moreover, BID and CASP9 were another two important genes that were down-regulated in shRIP2 + APEC vs. APEC, whereas these genes were up-regulated when APEC vs. WT were compared. BID is a member of the proapoptotic BCL2 family, which could stimulate mitochondrial outer membrane permeabilization to induce the release of regulators or other cell-death mediators . CASP9 is the member of cysteine aspartate-specific proteases, serving as an intrinsic initiator of apoptosis [40, 41]. Data obtained by qRT-PCR and western blot showed that the mRNA and protein of BID and CASP9 were both significantly decreased when shRIP2+ APEC vs APEC were compared (Fig. 10A, B, E and F). It has been demonstrated that BID and CASP9 were both involved in the apoptosis pathway that was downstream of RIP2 [42, 43]. Accordingly, our data point to BID and CASP9 as the targets of RIP2, since RIP2-depletion inhibits the expression of BID and CASP9 upon APEC infection. Apoptosis is a safety and important process to remove the damaged DNA, abnormal-proliferated or dedifferentiated cells, which is indispensable for maintaining cellular homeostasis, and normal regulation of the immune system [44, 45]. In the present study, the apoptosis was significantly activated in shRIP2 + APEC vs. APEC, indicating silence of RIP2 could induce apoptosis to inhibit tumorgenesis. These results were reasonable and expected, since knockdown of RIP2 could avoid the excessive tissues/cells injury and inflammatory response [42, 46]. In summary, our data suggest that BID and CASP9 could be the two important factors that determine immune and inflammatory response towards HD11 macrophages in the absence of RIP2 upon APEC challenge, even if the cellular fate in the presence of BID and CASP9 down-regulation is presently unknown. The specific regulation mechanism between RIP2 and BID/CASP9 involved in response to APEC challenge warrants further study.
In summary, this study has provided an analysis of the genetic landscape associated with RIP2 knockdown in HD11 macrophages in the presence or absence of an APEC challenge. In total, 4691 and 2605 DEGs were screened in shRIP2 + APEC vs. APEC and shRIP2 vs. WT, respectively. Functional annotation analysis showed that apoptosis, MAPK, p53, Toll-like receptor, and Nod-like receptor signaling pathways were involved in APEC-induced RIP2 knockdown HD11 cells. By analyzing the enriched pathway and gene networks, we identified that several key DEGs, like HSP90AB1, BID, and CASP9 were targeted by RIP2 upon APEC infection. As a whole, this study can not only provide data support for constructing gene networks of RIP2 knockdown associated with APEC challenge but also provide new ideas for improving the immune and inflammatory response.
Materials and methods
Chicken HD11 macrophage, an immortalized cell line, was selected as experimental materials as it is similar morphology and function to the in vivo macrophages and helpful to understand the host immune response in the early infection during APEC infection. The chicken HD11 cell line was kindly provided by Dr. Xuming Hu (Yangzhou University). The chicken macrophage-like cell line HD11 was maintained in RPMI1640 (Gibco, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (FBS, Gibco, Carlsbad, CA, USA) in a humidified incubator with 5% CO2 at 37 °C, and cells were passaged before 80–90% confluence.
Small hairpin RIP2 (shRIP2) plasmid (Table S1) were synthesized by GenePharma (Shanghai, China). The detailed information for the establishment of shRIP2 HD11 cells can be found in the study of Sun et al. . The Lipofectamine™ 2000 reagent (Invitrogen, Carlsbad, CA, USA) was used for the cell transfection according to the manufacturer’s instructions. After transfection with shRIP2 for 48 h, cells were challenged with or without APEC for 24 h, and collected for further study.
Quantitative real time PCR (qRT-PCR) analysis
Total RNA was isolated from cells of wild type group (WT), RIP2 knockdown group (shRIP2), APEC challenge group (APEC), and RIP2 knockdown with APEC challenge group (shRIP2 + APEC) using Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. Then the RNA was reverse transcribed into cDNA using a Reverse Transcription Kit (Takara, Dalian, China). The One Step SYBR® PrimeScript® PLUS RTRNA PCR Kit (Takara, Dalian, China) was used for cDNA synthesis. qRT-PCR was conducted using a SYBR® Premix Ex Taq™ II kit (Takara, Dalian, China) to evaluate the expression level of GAPDH, RIP2, HSP90AB1, BID, and CASP9. Primer sequences are displayed in Table S2. qRT-PCR thermal cycling conditions were as follows: denaturation for 3 min at 95 °C, 40 cycles of 10 s at 95 °C, 58 °C for 30 s, and then 72 °C for 30 s. Relative expression of above genes were calculated using the 2−ΔΔCt method and GAPDH was utilized as an internal control. The formula of ΔΔCt is (Ct of gene in test group - Ct of GAPDH in test group) - (Ct of gene in control group - Ct of GAPDH in control group).
Cells from each of the four groups were lysed on ice using 200 μL RIPA buffer (Beyotime Biotechnology, Shanghai, China) for 30 min. Next, the lysis mixtures were centrifuged and the supernatants were collected. BCA™ Protein Assay Kit (Pierce, Appleton, WI, USA) was used for quantification of proteins. Then, the isolated proteins were subjected to sodium dodecyl sulfate-polyacrylamide gel (SDS-PAGE) and electrophoretically transferred to PVDF membranes. Afterwards, membranes were blocked in 5% BSA for 2 h at room temperature and then probed with the primary antibodies at 4 °C overnight. The primery antibodies included anti-GAPDH (ab181602, Abcam, MA, USA), anti-RIP2 (70R-10,459, Fitzgerald, MA, USA), anti-HSP90AB1 (CBDH1187, Creative Biolabs, NY, USA), anti-BID (AB10002, MilliporeSigma, ON, Canada), and anti-CASP9 (STJ96979, St John’s Laboratory, London, UK) were used at a dilution of 1:1000. Then the membranes were incubated with secondary antibodies tagged with horseradish peroxidas (Sigma-Aldrich, St. Louis, MI, USA) at a 1:10000 dilution at room temperature for 2 h. Then, immunoblots were visualized by enhanced chemiluminescence (ECL kit, Santa Cruz Biotechnology, Dallas, TX, USA). The blots were visualized by using Image Lab™ Software (Bio-Rad, Hercules, CA, USA).
Nitric oxide (NO) production assay
Chicken HD11 cells from different groups (WT, shRIP2, APEC, and shRIP2 + APEC) were incubated for 24 h, then NO production in the cell supernatant was determined using the Griess reagent kit (Molecular Probes, Carlsbad, CA, USA). The cell supernatant was mixed with Griess reagents and incubated for 30 minutes in the dark, and then measured at 540 nm on a spectrophotometer. The absorbance values were compared to the sodium nitrite standard curve to determine nitrite concentrations (μM).
Cell viability assay
Cell counting kit-8 (CCK-8) was utilized to determine the viability of cells from different groups (WT, shRIP2, APEC, and shRIP2 + APEC). Cells from each of the four groups were placed in three replicates at a density of 1 × 105 cells per well in a 96-well plate with 100 μL of medium and incubated for 48 h. Then, the cells were incubated for 2 h in 10 μL of CCK-8 solution. Absorbance (optical density, OD) was assessed at 450 nm using a spectrophotometer. The experiment was performed in triplicate.
Total RNA was extracted from HD11 cells of different groups (WT, shRIP2, APEC, and shRIP2 + APEC) using an RNA isolation kit (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol. The quality of RNA was analyzed by agarose gel electrophoresis and a Nanodrop™ OneCspectrophotometer (Thermo Fisher Scientific Inc., MA, USA). The RNA Integrity was confirmed by Qseq (Qseq100, Guangding, Taiwan). A total of 2 μg of RNA was used for stranded RNA sequencing library preparation using Ribo-off rRNA depletion kit (Catalog NO. MRZG12324, Illumina, San Diego, CA, USA) following the manufacturer’s instruction. The library products corresponding to 200–500 bps were enriched, quantified and finally sequenced on NovaSeq 6000 sequencer (Illumina, San Diego, CA, USA) with PE150 sequencing platform.
Quality control and differentially expressed genes (DEGs) analysis
For quality control, raw data were first filtered by Trimmomatic (version 0.36) (i.e., removing low-quality reads (Q ≤ 10) and repeated and adaptor sequences (10 bp overlap (AGATCGGAAG)). The obtained clean reads were further treated with in-house scripts to eliminate duplication bias introduced in library preparation and sequencing. Deduplicated reads were mapped to the reference genome of the chicken (Gallus gallus) from Ensembl (https://asia.ensembl.org/info/data/ftp/index.html) using STRA software (version 2.5.3a) with default parameters. Reads mapped to the exon regions of each gene were counted by featureCounts (Subread-1.5.1; Bioconductor) and then RPKM was calculated. Genes differentially expressed between groups were identified using the edgeR package (version 3.12.1). A p-value cutoff of 0.05 and fold-change cutoff of 1.5 were used to judge the statistical significance of gene expression differences.
GO and KEGG pathway enrichment analyses
Statistical analysis was conducted using a one-way ANOVA and Tukey Honestly Significant differences test (HSD; SAS, 2000; Cary, NC) using JMP statistical software (version 15.2.1, SAS Institute). Data are expressed as the mean ± standard error. The statistical significance was defined at p < 0.05. The test results represent the data of four independent experiments.
Availability of data and materials
All data generated or analyzed for this study are included in this article and its supplementary files. The raw sequence reads were deposited into NCBI SRA database under accession no. PRJNA786725 (https://dataview.ncbi.nlm.nih.gov/object/PRJNA786725).
Receptor interacting serine/threonine kinase 2
Avian pathogenic E. coli
Mitogen-activated protein kinase
Heat shock protein 90 alpha family class B member 1
BH3 interacting domain death agonist
Nuclear factor kappa B
Tumor necrosis factor
C-X-C motif chemokine ligand 10
Kyoto encyclopedia of genes and genomes
Reads per kilobase per million
Jun proto-oncogene, AP-1 transcription factor subunit
- MAP K9:
Mitogen-activated protein kinase 9
Fos proto-oncogene, AP-1 transcription factor subunit
Phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta
Fas associated via death domain
CASP8 and FADD like apoptosis regulator
C-C motif chemokine ligand
Myeloid differentiation primary response 88
Toll like receptor 7
Siah E3 ubiquitin protein ligase 1
Phorbol-12-myristate-13-acetate-induced protein 1
Nolan LK, Barnes HJ, Vaillancourt JP, Abdul-Aziz T, Logue CM. Colibacillosis. Dis Poult. 2013:751–805 https://onlinelibrary.wiley.com/doi/abs/10.1002/9781119421481.ch18.
Stokholm NM, Permin A, Bisgaard M, Christensen JP. Causes of mortality in commercial organic layers in Denmark. Avian Dis. 2010;54:1241–50 https://pubmed.ncbi.nlm.nih.gov/21313846/.
Alber A, Morris KM, Bryson KJ, Sutton KM, Monson MS, Chintoan-Uta C, et al. Avian pathogenic Escherichia coli (APEC) strain-dependent immunomodulation of respiratory granulocytes and mononuclear phagocytes in CSF1R-reporter transgenic chickens. Front Immunol. 2020;10:3055 https://www.frontiersin.org/articles/10.3389/fimmu.2019.03055/full.
Nguyen LT, Thuan NK, Tam NT, Huyen-Trang CT, Khanh NP, Bich TN, et al. Prevalence and genetic telationship of predominant escherichia coli serotypes isolated from poultry, wild animals, and environment in the mekong delta, Vietnam. Vet Med Int. 2021;2021 https://pubmed.ncbi.nlm.nih.gov/34804471/.
Azam M, Mohsin M, Johnson TJ, Smith EA, Johnson A, Umair M, et al. Genomic landscape of multi-drug resistant avian pathogenic Escherichia coli recovered from broilers. Vet Microbiol. 2020;247:108766 https://pubmed.ncbi.nlm.nih.gov/32768218/.
Sun H, Liu P, Nolan LK, Lamont SJ. Thymus transcriptome reveals novel pathways in response to avian pathogenic Escherichia coli infection. Poult Sci. 2016;95:2803–14 https://pubmed.ncbi.nlm.nih.gov/27466434/.
Sun H, Liu P, Nolan LK, Lamont SJ. Avian pathogenic Escherichia coli (APEC) infection alters bone marrow transcriptome in chickens. BMC Genomics. 2015;16:1–15 https://pubmed.ncbi.nlm.nih.gov/26369556/.
Sun H, Liu P, Nolan LK, Lamont SJ. Novel pathways revealed in bursa of fabricius transcriptome in response to extraintestinal pathogenic Escherichia coli (ExPEC) infection. PLoS One. 2015;10:e0142570 https://pubmed.ncbi.nlm.nih.gov/26556806/.
Ghorbani AR, Khoshbakht R, Kaboosi H, Shirzad-Aski H, Ghadikolaii FP. Phylogenetic relationship and virulence gene profiles of avian pathogenic and uropathogenic Escherichia coli isolated from avian colibacillosis and human urinary tract infections (UTIs). Iran J Vet Res. 2021;22:203 https://pubmed.ncbi.nlm.nih.gov/34777520/.
da Rocha DT, de Oliveira SF, Borges KA, Furian TQ, do Nascimento VP, de Souza Moraes HL, Salle CTP. Avian pathogenic Escherichia coli (APEC) and uropathogenic Escherichia coli (UPEC): characterization and comparison. J Infect Dev Ctries. 2021;15:962–71 https://pubmed.ncbi.nlm.nih.gov/34343121/.
Hofmann SR, Girschick L, Stein R, Schulze F. Immune modulating effects of receptor interacting protein 2 (RIP2) in autoinflammation and immunity. Clin Immunol. 2021;223:108648 https://pubmed.ncbi.nlm.nih.gov/33310070/.
Humphries F, Yang S, Wang B, Moynagh PN. RIP kinases: key decision makers in cell death and innate immunity. Cell Death Differ. 2015;22:225–36 https://pubmed.ncbi.nlm.nih.gov/25146926/.
Wu XM, Chen WQ, Hu YW, Cao L, Nie P, Chang MX. RIP2 is a critical regulator for NLRs signaling and MHC antigen presentation but not for MAPK and PI3K/Akt pathways. Front Immunol. 2018;9:726 https://pubmed.ncbi.nlm.nih.gov/29692779/.
Park JH, Kim YG, McDonald C, Kanneganti TD, Hasegawa M, Body-Malapel M, et al. RICK/RIP2 mediates innate immune responses induced through Nod1 and Nod2 but not TLRs. J Immunol. 2007;178:2380–6 https://pubmed.ncbi.nlm.nih.gov/17277144/.
Lamkanfi M, Kalai M, Saelens X, Declercq W, Vandenabeele P. Caspase-1 activates nuclear factor of the κ-enhancer in B cells independently of its enzymatic activity. J Biol Chem. 2004;279:24785–93 https://pubmed.ncbi.nlm.nih.gov/15039421/.
Sarkar A, Duncan M, Hart J, Hertlein E, Guttridge DC, Wewers MD. ASC directs NFκB activation by regulating receptor interacting protein-2 (RIP2) caspase-1 interactions. J Immunol. 2006;176:4979–86 https://pubmed.ncbi.nlm.nih.gov/16585594/.
Gong Q, Long Z, Zhong FL, Teo DET, Jin Y, Yin Z, et al. Structural basis of RIP2 activation and signaling. Nat Commun. 2018;9:1–13 https://pubmed.ncbi.nlm.nih.gov/30478312/.
Chin AI, Dempsey PW, Bruhn K, Miller JF, Xu Y, Cheng G. Involvement of receptor-interacting protein 2 in innate and adaptive immune responses. Nature. 2002;416:190–4 https://pubmed.ncbi.nlm.nih.gov/11894097/.
Jeong Y, Kim C, Kim J, Oh S, Lee K, Park J, et al. RIP 2/RICK-dependent cytokine production upon yersinia enterocolitica infection in macrophages with TLR 4 deficiency. Scand J Immunol. 2013;78:401–7 https://pubmed.ncbi.nlm.nih.gov/23952047/.
Reinke S, Linge M, Diebner HH, Luksch H, Glage S, Gocht A, et al. Non-canonical caspase-1 signaling drives RIP2-dependent and TNF-α-mediated inflammation in vivo. Cell Rep. 2020;30:2501–11 https://pubmed.ncbi.nlm.nih.gov/32101731/.
Homer CR, Kabi A, Marina-García N, Sreekumar A, Nesvizhskii AI, Nickerson KP, et al. A dual role for receptor-interacting protein kinase 2 (RIP2) kinase activity in nucleotide-binding oligomerization domain 2 (NOD2)-dependent autophagy. J Biol Chem. 2012;287:25565–76 https://pubmed.ncbi.nlm.nih.gov/22665475/.
Caldwell RB, Kierzek AM, Arakawa H, Bezzubov Y, Zaim J, Fiedler P, et al. Full-length cDNAs from chicken bursal lymphocytes to facilitate gene function analysis. Genome Biol. 2005;6:1–9 https://pubmed.ncbi.nlm.nih.gov/15642098/.
Nachbur U, Stafford CA, Bankovacki A, Zhan Y, Lindqvist LM, Fiil BK, et al. A RIPK2 inhibitor delays NOD signalling events yet prevents inflammatory cytokine production. Nat Commun. 2015;6:1–13 https://pubmed.ncbi.nlm.nih.gov/25778803/.
Zhang K, Wang L, Hong X, Chen H, Shi Y, Liu Y, et al. Pulmonary alveolar stem cell senescence, apoptosis, and differentiation by p53-dependent and independent mechanisms in telomerase-deficient mice. Cells. 2021;10:2892 https://pubmed.ncbi.nlm.nih.gov/34831112/.
Qian Y, Chen X. Senescence regulation by the p53 protein family. Methods Mol Biol. 2013:37–61 https://pubmed.ncbi.nlm.nih.gov/23296650/.
Reinhardt HC, Schumacher B. The p53 network: cellular and systemic DNA damage responses in aging and cancer. Trends Genet. 2012;28:128–36 https://pubmed.ncbi.nlm.nih.gov/22265392/.
Hung JC, Wu JL, Li HC, Chiu HW, Hong JR. The proapoptotic gene bad regulates brain development via p53-mediated stress signals in zebrafish. Cells. 2021;10:2820 https://pubmed.ncbi.nlm.nih.gov/34831043/.
Wen Y, Yang Z, Song M, Li B, Zhu J, Wang E. SIAH1 induced apoptosis by activation of the JNK pathway and inhibited invasion by inactivation of the ERK pathway in breast cancer cells. Cancer Sci. 2010;101:73–9 https://pubmed.ncbi.nlm.nih.gov/19775288/.
Morsi RZ, Hage-Sleiman R, Kobeissy H, Dbaibo G. Noxa: role in cancer pathogenesis and treatment. Curr Cancer Drug Targets. 2018;18:914–28 https://pubmed.ncbi.nlm.nih.gov/29521234/.
Mancini F, Pieroni L, Monteleone V, Lucà R, Fici L, Luca E, et al. MDM4/HIPK2/p53 cytoplasmic assembly uncovers coordinated repression of molecules with anti-apoptotic activity during early DNA damage response. Oncogene. 2016;35:228–40 https://pubmed.ncbi.nlm.nih.gov/25961923/.
Zare A, Petrova A, Agoumi M, Armstrong H, Bigras G, Tonkin K, et al. RIPK2: new elements in modulating inflammatory breast cancer pathogenesis. Cancers (Basel). 2018;10:184 https://pubmed.ncbi.nlm.nih.gov/29874851/.
Zhang H, Chin AI. Role of Rip2 in development of tumor-infiltrating MDSCs and bladder cancer metastasis. PLoS One. 2014;9:e94793 https://pubmed.ncbi.nlm.nih.gov/24733360/.
Hu Y, Jiao B, Wang C, Wu J. Regulation of temozolomide resistance in glioma cells via the RIP2/NF-κB/MGMT pathway. CNS Neurosci Ther. 2021;27:552–63 https://pubmed.ncbi.nlm.nih.gov/33460245/.
Buchner J, Li J. Structure, function and regulation of the hsp90 machinery. Biom J. 2013;36:106 https://pubmed.ncbi.nlm.nih.gov/23806880/.
Radli M, Veprintsev DB, Rüdiger SGD. Production and purification of human Hsp90β in Escherichia coli. PLoS One. 2017;12:e0180047 https://pubmed.ncbi.nlm.nih.gov/28651008/.
Chen ZY, Gan JK, Xiao X, Jiang LY, Zhang XQ, Luo QB. The association of SNPs in Hsp90β gene 5′ flanking region with thermo tolerance traits and tissue mRNA expression in two chicken breeds. Mol Biol Rep. 2013;40:5295–306 https://pubmed.ncbi.nlm.nih.gov/23793829/.
Kuo CC, Liang CM, Lai CY, Liang SM. Involvement of heat shock protein (Hsp) 90β but not Hsp90α in antiapoptotic effect of CpG-B oligodeoxynucleotide. J Immunol. 2007;178:6100–8 https://pubmed.ncbi.nlm.nih.gov/17475835/.
Bhat A, Gomis S, Potter A, Tikoo SK. Role of Hsp90 in CpG ODN mediated immunostimulation in avian macrophages. Mol Immunol. 2010;47:1337–46 https://pubmed.ncbi.nlm.nih.gov/20096933/.
Yip KW, Reed JC. Bcl-2 family proteins and cancer. Oncogene. 2008;27:6398–406 https://pubmed.ncbi.nlm.nih.gov/18955968/.
Kim YR, Kim KM, Yoo NJ, Lee SH. Mutational analysis of CASP1, 2, 3, 4, 5, 6, 7, 8, 9, 10, and 14 genes in gastrointestinal stromal tumors. Hum Pathol. 2009;40:868–71 https://pubmed.ncbi.nlm.nih.gov/19269008/.
Kumar S. Mechanisms mediating caspase activation in cell death. Cell Death Differ. 1999;6:1060–6 https://pubmed.ncbi.nlm.nih.gov/10578174/.
Yang Q, Tian S, Liu Z, Dong W. Knockdown of RIPK2 inhibits proliferation and migration, and induces apoptosis via the NFκB signaling pathway in gastric cancer. Front Genet. 2021;12:84 https://pubmed.ncbi.nlm.nih.gov/33633788/.
Cai X, Yang Y, Xia W, Kong H, Wang M, Fu W, et al. RIP2 promotes glioma cell growth by regulating TRAF3 and activating the NF-κB and p38 signaling pathways. Oncol Rep. 2018;39:2915–23 https://pubmed.ncbi.nlm.nih.gov/29693188/.
Lessene G, Czabotar PE, Colman PM. BCL-2 family antagonists for cancer therapy. Nat Rev Drug Discov. 2008;7:989–1000 https://pubmed.ncbi.nlm.nih.gov/19043450/.
Seo HW, Rengaraj D, Choi JW, Park KJ, Song G, Han JY. The expression profile of apoptosis-related genes in the chicken as a human epithelial ovarian cancer model. Oncol Rep. 2011;25:49–56 https://pubmed.ncbi.nlm.nih.gov/21109956/.
Wu S, Kanda T, Nakamoto S, Imazeki F, Yokosuka O. Knockdown of receptor-interacting serine/threonine protein kinase-2 (RIPK2) affects EMT-associated gene expression in human hepatoma cells. Anticancer Res. 2012;32:3775–83 https://pubmed.ncbi.nlm.nih.gov/22993319/.
Sun H, Li N, Tan J, Li H, Zhang J, Qu L, et al. Transcriptional regulation of RIP2 gene by NFIB is associated with cellular immune and inflammatory response to APEC infection. Int J Mol Sci. 2022;23:3814. https://doi.org/10.3390/ijms23073814.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13 https://pubmed.ncbi.nlm.nih.gov/19033363/.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57 https://pubmed.ncbi.nlm.nih.gov/19131956/.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30 https://pubmed.ncbi.nlm.nih.gov/10592173/.
We would like to thank Xuming Hu who provided the HD11 cells for this experiment.
This research was supported by the National Natural Science Foundation of China (Grant No. 31802053), The Natural Science Foundation of Jiangsu Province (Grant No. BK20180907), the China Postdoctoral Science Foundation (2019 M661950), and Jiangsu Postdoctoral Science Foundation (137070510).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing financial interest in the work described.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
shRNA information for RIP2. Table S2. Primers for candidate genes in qRT-PCR validation experiment. Table S3. Characteristics of RNA sequencing data before or after quality control. Table S4. The reads mapping information for each sample. Table S5. The significantly enriched gene ontology (GO) in the comparison of shRIP2 vs. WT with adjusted p-value ≤0.05. Table S6. The significantly enriched gene ontology (GO) in the comparison of shRIP2 + APEC vs. APEC with adjusted p-value ≤0.05. Table S7. The significantly enriched pathways in the comparison of shRIP2 vs. WT with adjusted. p_value ≤0.05. Table S8. The significantly enriched pathways in the comparison of shRIP2 + APEC vs. APEC with adjusted p-value ≤0.05. Table S9. The differentially expressed genes involved in the MAPK signaling pathway in the comparison of shRIP2 vs. WT. Table S10. The differentially expressed genes involved in the MAPK signaling pathway in the comparison of shRIP2 + APEC vs. APEC. Table S11. The differentially expressed genes involved in the apoptosis pathway in the comparison of shRIP2 vs. WT. Table S12. The differentially expressed genes involved in the apoptosis pathway in the comparison of shRIP2 + APEC vs. APEC.
Veen diagram for the MAPK signaling pathway in the contrast of knockdown of RIP2 HD11 cells (shRIP2) vs. wild type HD11 cells (WT) and knockdown of RIP2 combined with avian pathogenic E. coli infection HD11 cells (shRIP2 + APEC) vs. avian pathogenic E. coli infection HD11 cells (APEC), as well as the apoptosis pathway in the contrast of shRIP2 vs. WT and shRIP2 + APEC vs. APEC. Figure S2. Full-length blot of GAPDH protein expression level for Fig. 1B. Figure S3. Full-length blot of RIP2 protein expression level for Fig. 1B. Figure S4. Full-length blot of BID protein expression level for Fig. 10D. Figure S5. Full-length blot of CASP9 protein expression level for Fig. 10D. Figure S6. Full-length blot of HSP90AB1 protein expression level for Fig. 10D. Figure S7. Full-length blot of GAPDH protein expression level for Fig. 10D.
About this article
Cite this article
Sun, H., Yang, Y., Cao, Y. et al. Gene expression profiling of RIP2-knockdown in HD11 macrophages — elucidation of potential pathways (gene network) when challenged with avian pathogenic E.coli (APEC). BMC Genomics 23, 341 (2022). https://doi.org/10.1186/s12864-022-08595-5
- HD11 macrophages
- Gene network