Skip to main content

Analysis of differentially expressed long non-coding RNAs in LPS-induced human HMC3 microglial cells

Abstract

Background

Long non-coding RNAs (lncRNAs) are emerging as key modulators of inflammatory gene expression, but their roles in neuroinflammation are poorly understood. Here, we identified the inflammation-related lncRNAs and correlated mRNAs of the lipopolysaccharide (LPS)-treated human microglial cell line HMC3. We explored their potential roles and interactions using bioinformatics tools such as gene ontology (GO), kyoto encyclopedia of genes and genomes (KEGG), and weighted gene co-expression network analysis (WGCNA).

Results

We identified 5 differentially expressed (DE) lncRNAs, 4 of which (AC083837.1, IRF1-AS1, LINC02605, and MIR3142HG) are novel for microglia. The DElncRNAs with their correlated DEmRNAs (99 total) fell into two network modules that both were enriched with inflammation-related RNAs. However, treatment with the anti-inflammatory agent JQ1, an inhibitor of the bromodomain and extra-terminal (BET) protein BRD4, neutralized the LPS effect in only one module, showing little or even enhancing effect on the other.

Conclusions

These results provide insight into, and a resource for studying, the regulation of microglia-mediated neuroinflammation and its potential therapy by small-molecule BET inhibitors.

Peer Review reports

Background

Microglia, the macrophages of the central nervous system (CNS), play crucial roles in its homeostasis and immune defense [1, 2]. In response to inflammatory stimuli such as lipopolysaccharide (LPS), these cells become polarized to the M1 phenotype and produce pro-inflammatory cytokines and oxidative metabolites such as IL-1β, TNF, IL-6, and nitric oxide [3].

The bromodomain and extra-terminal (BET) family proteins are epigenetic readers that control the inflammatory response by regulating the expression of transcription factors and cytokines in T cells, monocytes, and macrophages, and therefore are potential therapeutic targets [4,5,6]. Small-molecule pan-BET inhibitors such as JQ1 and I-BET151 protected mice against LPS-induced sepsis [7] and attenuated the LPS-stimulated expression of pro-inflammatory genes in murine bone marrow-derived macrophages [8]. Similarly, JQ1 attenuated pro-inflammatory chemokine, cytokine, and interferon response genes in the LPS-treated murine glial cell line BV-2 [9], an observation that we confirmed and expanded in the human microglial cell line HMC3 [10]. Thus, it appears that the BET proteins also regulate neuroinflammation.

Our previous study of LPS/JQ1-treated HMC3 cells [10] determined their transcriptomes by RNA sequencing (RNA-seq), but limited the analysis to differentially expressed mRNAs (DEmRNAs), i.e. to protein-coding genes. However, long non-coding RNAs (lncRNAs), such as lincRNA-Cox2 [11], PACER [12], and THRIL [13], are emerging as additional players in gene regulation [14, 15], including regulation of inflammatory genes [16,17,18]. Accumulating studies have shown that many lncRNAs, including Lethe [19], NEAT1 [20], AS-IL1α [21], and FIRRE [22], play crucial roles in the immune system by regulating excessive or uncontrolled inflammation.

The present study therefore extends our transcriptomic analysis of the LPS/JQ1-treated HMC3 cells into the category of differentially expressed lncRNAs (DElncRNAs). By re-analyzing the published RNA-seq datasets [10], we identified (i) inflammation-related and BET inhibitor-sensitive DElncRNAs, (ii) their correlated DEmRNAs, and (iii) potential functional networks that contain these two classes of transcripts.

Results

Approach

We first identified the LPS-induced DElncRNAs in the dataset GSE155408. Compared to the initial study [10], we added processing methods that are specifically designed to identify DElncRNAs (Fig. 1A); furthermore, we now used a stricter cutoff, which slightly downsized the DEmRNA data set (Supplementary Fig. 1A). We then identified the DEmRNAs that were correlated with the LPS-induced DElncRNAs. In parallel, we assessed whether or how the LPS-induced DE RNAs (including lncRNAs and mRNAs) were affected by JQ1. Next, we performed a network analysis in order to construct interaction modules incorporating the LPS-induced DElncRNAs and their correlated DEmRNAs. This was followed by pathway analysis in order to learn about the potential functional significance of these modules. Finally, we constructed a corresponding protein–protein interaction network.

Fig. 1
figure 1

Workflow and identification of DElncRNAs. A Schematic summarizing project workflow. B Volcano plots visualizing the effects of LPS and/or JQ1 on the lncRNA levels in HMC3 human microglial cells; the treatments were performed in triplicates and the combined results are shown. The gray vertical dot lines indicate log2 FC of ≥ 1.2 and ≤ -1.2. The gray horizontal dot lines indicate –log(padj) ≤ 0.01. C Heat map showing the expression changes of the 5 LPS-induced lncRNAs that qualified as DElncRNAs by the criteria of this study (see Materials and methods), and the heat maps for the same DElncRNAs upon treatment with JQ1 or LPS + JQ1. The color scale represents the log2 FC values. DElncRNAs; differentially expressed lncRNAs, DEmRNAs; differentially expressed mRNAs, log2 FC; log2 fold change, GO; Gene Ontology, KEGG; Kyoto Encyclopedia of Genes and Genomes

Identification of DElncRNAs and DEmRNAs in LPS-treated HMC3 cells

In the control vs. LPS-treated cells, we identified a total of 5 DElncRNAs (all upregulated; Fig. 1B and C) and 99 DEmRNAs (98 upregulated and 1 downregulated; Supplementary Fig. 1A). Conversely, the DElncRNAs and a majority of the DEmRNAs tended to be downregulated or unaffected by JQ1. When LPS and JQ1 were combined, JQ1 further increased the expression of the two DElncRNAs that were most increased by LPS (AC083837.1 and LINC02605), but partially (MIR3142HG) or fully (MIR155HG and IRF1-AS1) neutralized the LPS effect on the other DElncRNAs (Fig. 1C). Likewise, JQ1 counteracted LPS or did not alter its effect for the majority of DEmRNAs, but enhanced the LPS effect in some cases. We note that the above DElncRNAs had not yet been identified in the HMC3 cells; the DEmRNAs are often inflammation- and immunity-related (e.g., CCL20, CSF3, CXCL10, TNF, and CXCL8) (Supplementary Fig. 1B and C). The up- and downregulated DElncRNAs and DEmRNAs are listed in Supplementary Table 1.

Correlations between DElncRNAs and DEmRNAs

The heat map in Fig. 2 visualizes the Pearson correlation coefficient (r) between the DElncRNAs and DEmRNAs; numerical values are listed in Supplementary Table 2. When including all 5 LPS-induced DElncRNAs and all 99 LPS-induced DEmRNAs into the analysis, a total of 211 DElncRNA-DEmRNA pairs (defined as |r|≥ 0.75 and padj ≤ 0.01) were identified. Specifically, MIR155HG, IRF1-AS1, AC083837.1, LINC02605, and MIR3142HG were paired with 73, 57, 36, 21, and 24 DEmRNAs, respectively.

Fig. 2
figure 2

Correlation heat map of the DElncRNAs and DEmRNAs. Each column corresponds to one of the LPS-induced DElnRNAs, and the cells of the rows show for each LPS-induced DEmRNA the Pearson correlation (r value) with that given DElncRNA both in color code (according to the color scale shown on the right) and as the numerical value. A red cells indicate the positive correlation, while blue cells indicate the negative correlation, with higher correlation indicated by dark color intensity as shown by the color scale

Identification of DElncRNA-DEmRNA network modules

In order to find potential interactions, we constructed co-expression networks based on the above correlations between the DElncRNAs and DEmRNAs. The resulting nodes and relations fell into two modules, which we dubbed “large turquoise” and “small cyan” (Fig. 3). The large turquoise module was defined by IRF1-AS1 and MIR155HG (the two DElncRNAs that were least upregulated by LPS); the small cyan module was defined by AC083837.1, LINC02605, and MIR3142HG (the DElncRNAs that were more strongly induced by LPS). Within these modules, MIR155HG has the maximum number of co-expressed genes.

Fig. 3
figure 3

Co-expression networks of DElncRNAs with DEmRNAs. A network was constructed using DElncRNA-DEmRNA pairs with Pearson correlation coefficient (r) value ≥ 0.75. We identified 632 connections between the 5 DElncRNAs and 99 DEmRNAs, which could be organized into a large turquoise (left) and a small cyan (right) module. Purple circles denote DElncRNAs, turquoise and cyan rectangles denote DEmRNAs. The yellow rectangles denote DEmRNAs that have core interactions with DElncRNAs. Solid lines indicate correlations

We determined by chi-square test whether there was a significant difference in how the two modules respond to the JQ1 treatment (Supplementary Table 3). The mRNAs of the large turquoise module were significantly related to JQ1, while the mRNAs of the small cyan module were not related to JQ1. Specifically, out of 65 mRNAs in the large turquoise module, the levels of 64 mRNAs were decreased by JQ1. In contrast, out of 29 mRNAs of the small cyan module, 11 were marginally affected and 18 were increased. In more detail, the mRNA levels of the small cyan module that were highly correlated with MIR3142HG were reduced or marginally affected by JQ1, while the mRNAs that were highly correlated with AC083837.1 and LINC02605 showed increased expression or marginal changes.

To validate the two modules, which were identified on the basis of RNA-seq data, we repeated the LPS/JQ1 treatments of the HMC3 cells and measured the expression levels of the DElncRNAs and selected (i.e., highly correlated, inflammation-related) DEmRNAs. In the large turquoise module, JQ1 mostly neutralized the LPS effect (22 of 23 tested mRNAs, both tested lncRNAs) (Fig. 4 and Supplementary Fig. 2), whereas in the small cyan module, most of the LPS-stimulated RNAs were marginally affected or even further increased by the additional presence of JQ1 (15 of 17 mRNAs, 1 of 2 tested lncRNAs) (Fig. 5 and Supplementary Fig. 2). In this study, primers for qRT-PCR of lncRNAs were designed based on the sequences of IRF1-AS1-203 (ENST00000378953.8), MIR155HG-201 (ENST00000456917.2), and MIR3142HG-201 (ENST00000517927.1) transcripts.

Fig. 4
figure 4

qRT-PCR validation of the large turquoise module. A Co-expression network image highlighting the mRNAs of the large turquoise module. Of the genes in the large turquoise module, only selected mRNAs are shown in bold and in red. B A list of selected (highly correlated, inflammation-related) DElncRNA-DEmRNA pairs in the large turquoise module. C The HMC3 cells were treated with DMSO, LPS, and/or with JQ1, and selected DElncRNAs and DEmRNAs of the large turquoise module assessed by qRT-PCR. The RNA levels are normalized to GAPDH transcript levels. Dark green indicates lncRNA and light green indicates mRNA. The data represent three independent experiments. The values are the mean ± SEM of triplicate experiments (*p < 0.05 and **p < 0.001)

Fig. 5
figure 5

qRT-PCR validation of the small cyan module. A Co-expression network image highlighting the mRNAs of the small cyan module. Of the genes in small cyan module, only selected mRNAs are shown in bold and in red. B A list of selected DElncRNA-(highly correlated, inflammation-related) DEmRNA pairs in the small cyan module. C The HMC3 cells were treated with DMSO, LPS, and/or with JQ1, and selected DElncRNAs and DEmRNAs of the small cyan module assessed by qRT-PCR. The RNA levels are normalized to GAPDH transcript levels. Dark blue indicates lncRNA and light blue indicates mRNA. The data represent three independent experiments. The values are the mean ± SEM of triplicate experiments (*p < 0.05 and **p < 0.001)

To summarize, while the mRNAs of the large turquoise module mostly were decreased by JQ1, the mRNAs of the small cyan module showed a mixed response pattern.

Functional annotations

Having validated the network modules, we performed a functional classification and pathway enrichment analysis of their respective mRNAs. In both the large turquoise (Fig. 6A) and small cyan (Fig. 6C) modules, Gene Ontology (GO) analysis highlighted terms related to inflammatory response such as type I interferon signaling pathway, defense response to virus, and chemokine-mediated signaling pathway. Similarly, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (Fig. 6B and D) found a total of 32 KEGG pathways (FDR ≤ 0.05, DEmRNA counts ≥ 4) predominantly related to inflammation, such as TNF, IL-17, NOD-like receptor, TLR, and NF-kappa B signaling. The main functions enriched in both modules are related to immune and inflammation. However, we found that apoptosis and cell adhesion-related functions were more enriched in the large turquoise module. The GO and KEGG analysis results of the total mRNAs are shown in Supplementary Fig. 1B and C.

Fig. 6
figure 6

GO and KEGG pathway analyses of the DElncRNAs and DEmRNAs. Shown are the GO term (A, C) and KEGG pathway enrichment (B, D) analyses of the large turquoise (A, B) and small cyan (C, D) modules. In the GO term analyses, the numbers of genes and false discovery rate (FDR) values are displayed for the top 5 GO terms in biological process (BP; upper panel), cellular component (CC; middle panel), and molecular function (MF; bottom panel). The blue column is the count value indicating the number of genes enriched in the GO term, and the red line is the -log10 (FDR) value. In the KEGG pathway enrichment analyses, each row represents an enriched function, and the size of the bubble represents the p-value (KOBAS, http://kobas.cbi.pku.edu.cn). The KOBAS algorithm divides the clusters according to the values computed for the enriched pathway, and the color of each bubble represents a different cluster

Protein–protein interaction (PPI) network analysis

The PPI network mathematically calculates the physical interactions between proteins in cells, which allowing a molecular assessment at the molecular and system level. We constructed a PPI network of DEmRNAs, which highly correlated DElncRNAs using the STRING database (Fig. 7). The PPI network contained 63 nodes and 332 edges in the co-expressed mRNAs of the large turquoise module (related to the MIR155HG and IRF1-AS1 lncRNAs) and 29 nodes and 90 edges in the co-expressed mRNA of the small cyan module (related to the AC083837.1, LINC02605, and MIR3142HG lncRNAs). The top 5 mRNAs that had the highest node degrees in the large turquoise module were IL1B (degree = 33), IRF1 (degree = 31), IL6 (degree = 29), CXCL10 (degree = 28), and TNFSF10 (degree = 27). The top 5 mRNAs that had the highest node degrees in the small cyan module were TNF (degree = 21), CXCL8 (degree = 18), NFKBIA (degree = 18), CXCL1 (degree = 12), and CXCL2 (degree = 12).

Fig. 7
figure 7

Protein–protein interaction (PPI) network analysis of DEmRNAs. Shown are the PPI analyses of the large turquoise (A) and small cyan (B) modules. The potential interactions of the proteins encoded by DElncRNA-correlated DEmRNA were determined by STRING online software, using a combined score > 0.7 as cut-off criterion. Lines indicate associations/interactions between genes. The larger the number of connections of a given gene (“node”), the higher the connectivity. A high connectivity indicates the importance of a gene in the PPI network. The color represents the log2 FC values and the size of the bubble represents the p-value

Discussion

In this study, we identified 5 DElncRNAs in LPS-treated HMC3 human microglial cells. 4 of these transcripts (AC083837.1, IRF1-AS1, LINC02605, and MIR3142HG) have not previously been invoked in the microglial inflammatory response. We further found that the BET inhibitor JQ1 can exert both positive and negative effects on the LPS-induced DElncRNAs. Importantly, we were able to correlate the novel DElncRNAs with inflammation-related DEmRNAs, in line with the idea that the DElncRNAs modulate the microglial immune response. More specifically, we found that the LPS-induced DElncRNAs and their correlated mRNAs fall into two modules that we dubbed “large turquoise” (defined by IRF1-AS1 and MIR155HG) and “small cyan” (defined by AC083837.1, LINC02605, and MIR3142HG). Interestingly, JQ1 affected their patterns in opposite ways: Inflammation processes were decreased by JQ1 in the large turquoise module, but marginally affected or even significantly increased in the small cyan module. We performed qRT-PCR to verify the expression of lncRNAs and mRNAs of the divided modules based on the Pearson correlation coefficient (r) of DElncRNAs-DEmRNAs. We found some genes with conflicting DElncRNAs-DEmRNAs correlations as a result of qRT-PCR, which appeared as a probability of false-positive in Pearson's correlation analysis. Pearson correlation analysis is the commonly used analysis method to measure the strength of association in a linear relationship between two datasets based on the mapped read count value. Although the Pearson correlation analysis result is reliable, it can have false positives due to the influence of several variables such as the magnitude of the slope at which the points are clustered and the heteroscedasticity [23], and thus results of RNA-seq analysis and qRT-PCR may be different. Therefore, verification of qRT-PCR of RNA-seq data is necessary.

As a result of co-expression network analysis, PPI network analysis, and functional annotation analysis, we found that genes and biological terms related to immune response, inflammatory responses, cytokine, and chemokine activity were enriched in both modules, and TNFα and NF-κB signaling-related genes were also found in both modules. Interestingly, in terms of inflammatory response, the large turquoise module was enriched with interferon signaling-related functions, and the small cyan module was enriched with the CXC chemokine in particular. Regarding cell proliferation, survival, and metabolism, the large turquoise module was associated with negative regulation of metabolic processes, and positive regulation of apoptotic processes such as TRAIL signaling and necroptosis. On the other hand, the small cyan module showed positive regulation of metabolic processes and more enriched MAPK signaling regulation including activated tak1 mediates p38 mapk activation compared to the large turquoise module.

As summarized above, TNFα-related genes (including also some TNF family members themselves) and NF-κB-related genes (including also NF-κB components) are present in both modules. As such, this is in line with the known roles of the TNFs and NF-κB as global regulators of inflammation, i.e. as regulators of multiple aspects of inflammation. However, the response of the TNFα- or NF-κB-related genes to JQ1 is different between the two modules, such that the TNFα- or NF-κB-related genes are downregulated along with the IFN-related genes in the large turquoise module but upregulated or unchanged along with the CXC chemokine genes in the small cyan module. Hence, different aspects of inflammation regulated by TNFs or NF-κB respond differently to JQ1, hence are under distinct epigenetic controls. Given that the TNFs and NF-κB are each related to more than one signaling pathway [24], it will be of interest to relate the large turquoise vs. small cyan modules to those specific pathways.

A complementary message emerges from the findings that we made with the more specialized inflammatory regulators, namely the IFNs and CXC chemokines. The direct biological effects of IFNs (mediated mainly by JAK/STAT signaling) are related especially, albeit not exclusively, to the modulation of immune system [25]; the direct biological effect of CXC chemokines is especially, albeit not exclusively, to promote neutrophil migration [26].

Here, we found that along with the TNFα/NF-κB genes of the large turquoise module, the IFN (mainly type I and II)-related genes (including IFNB1), which were limited to that same module, were also inhibited by JQ1. This agrees with the literature, since the direct inhibition of IFN response by BET inhibitors has been well documented [27]. Furthermore, we found that along with the TNFα/NF-κB-related genes of the small cyan module, the CXC chemokine genes (encoding CXCL1/2/3/5/8, all known to bind the chemokine receptor CXCR2), which were limited to that same module, were not inhibited by JQ1. We note in this respect that the GO term “activated tak1 mediated p38 mapk signaling activation” was also enriched in the small cyan module, in line with the known role of tak1/p38 signaling in CXC chemokine signaling [28]. Thus, not the type of cytokine per se (TNF vs. IFN vs. CXC chemokine) determines the outcome of JQ1 treatment, but the belonging to the large turquoise vs. small cyan module.

The selectivity of the gene response to I-BET has been reported to be related to the epigenetic status of the responding gene and the mechanism of BET recruitment. Nicodeme E et al. reported that some LPS-induced cytokines and chemokines, such as Tnf, Ccl2-5 and Cxcl1/2, have a highly selective effect on I-BET in bone marrow-derived macrophages. Further study on the mechanism of selectivity for BET inhibitors is needed, and considering the previous reports, selectivity for BET inhibitors may be related to the BET recruitment pathway or histone acetylation level for each gene, and the binding of the BET protein to the gene promoter or super-enhancer.

Of note, the DElncRNAs that we identified in the LPS-stimulated HMC3 cells have been observed in other inflammation-related contexts before. IRF1-AS1 (Lnc-SLC22A5-6) was found to act as a positive modulator of the IFN response in esophageal squamous cell carcinoma (ESCC), functioning as a tumor suppressor by regulating cell proliferation and apoptosis [29]. MIR155HG was increased in human M1 (inflammatory-type) macrophages that were derived from monocytes by treatment with LPS and IFN-gamma [30]. In human B cell lymphoma cells, LPS induced the nuclear translocation of an NF-κB p50/p65 heterodimer that could bind to the MIR155HG promoter, suggesting that MIR155HG was a direct NF-κB target gene [31]. MIR3142HG, the host gene for miR-3142 and miR-146a, was reported to regulate the IL-1β-induced inflammatory response in idiopathic pulmonary fibrosis (IPF) [32] and to be highly expressed in LPS-exposed human pulmonary microvascular endothelial cells (HPMECs) [33]. It is therefore of interest that our study showed that the LPS-induced increase of MIR3142HG was reduced by JQ1. However, the expression levels of MIR3142HG-correlated genes are controversial. In our study, JQ1 was not or only marginally effective on the MIR3142HG-correlated genes. Finally, LINC02605 (IL7-AS; Lnc-ZC2HC1A-1) was reported to regulate the immune response in a cell-type specific manner [34]. IL7-AS is a positive regulator of the IL1β-induced inflammatory response in human A549 epithelial cells, but a negative regulator in LPS-stimulated human THP-1 monocytes and mouse RAW 264.7 macrophages. Additionally, knockdown of IL7-AS increased the IL-6 release in IPF-derived fibroblasts, indicating that it is a negative regulator [32]. Our co-expression network analysis showed that while LINC02605 was induced by LPS, neither LINC02605 nor its correlated genes were significantly affected by JQ1. Interestingly, this was observed for all three DElncRNAs and their correlated genes in the small cyan module. These results are consistent with our previous transcriptome analyses of human and mouse microglial cells and mouse bone marrow-derived macrophages [8,9,10].

Conclusion

Here, we present a comprehensive analysis of inflammation-related DElncRNA and DEmRNA expression profiles and functional networks in the human microglial cell line HMC3. We identified 5 DElncRNAs (including 4 novel ones in microglia) and 99 DEmRNAs. We constructed DElncRNA-DEmRNA co-expression networks, which fell into two separate modules, and investigated their functions and pathways, which – for both modules—turned out as largely known to be inflammation-related. We determined that although considered as an anti-inflammatory agent, the BET inhibitor JQ1 regulates the two modules differently, showing an almost uniform anti-inflammatory effect on one module, but little or even enhancing effect on the other. This interesting result will have to be complemented and validated with methods such as knockdown and overexpression; it is also of interest whether it can be extended to other models of neuroinflammation. Altogether, the RNA expression modules that we identified here provide a resource for further studies of human microglial neuroinflammation through both computational analysis and functional approaches.

Materials and methods

Identification of differentially expressed lncRNAs and mRNAs

For this study, we used the RNA-seq data of our previous paper (GSE155408). To identify DElncRNAs, a comprehensive reference list of known lncRNAs was included in the processing of the RNA-seq data [10, 35, 36]. Briefly, FASTQ data were quality controlled and trimmed with Trimmomatic (version 0.36) [37]. The FASTQ files were aligned using STAR (version 2.7.8) [38] alignment software with the GENCODE Homo sapiens reference sequence GRCh38 (Release 27). DElncRNAs and DEmRNAs were normalized to sequencing depth and RNA composition using the median method with default parameters of DESeq2 [39]. Differential expression analysis of lncRNA and mRNA was conducted using the DESeq2 R package. The DElncRNAs and DEmRNAs were selected with a cutoff of | log2 fold change (log2 FC) | ≥ 1.2, | log2 FC | ≤ -1.2, and adjusted p-value (padj) ≤ 0.05 in LPS-treated HMC3 cells.

Weighted gene co-expression network analysis (WGCNA)

First, the Pearson correlation coefficient (r) values were calculated to assess the similarity of the expression patterns of transcripts. Then, a scale-free network was obtained by weighting the correlation coefficient between transcripts with soft-thresholding power. A module is defined as a cluster of densely interconnected transcripts in terms of co-expression. We considered a |r|≥ 0.75 as a meaningful value. Cytoscape MCODE plug-in (Version 3.4.0, available online: http://www.cytoscape.org/) [40] was applied for visualization of the co-expression networks.

Functional annotation and canonical pathway analysis

Database for Annotation, Visualization, and Integrated Discovery (DAVID, version 6.8) software (http://david.abcc.ncifcrf.gov/home.jsp) was used to analyze the biological functions in the datasets [41]. DAVID uses a modified Fisher's exact p-value to examine gene ontology (GO) enrichment. A false discovery rate (FDR) ≤ 0.05 was used as the criterion for GO term analysis. The KEGG Orthology Based Annotation System (KOBAS, version 3.0) software (http://kobas.cbi.pku.edu.cn/) [42] was used to analyze the enriched KEGG pathways [43] in the datasets. FDR ≤ 0.05 was used as the criterion for KEGG pathway enrichment analysis.

Cell culture and treatment

HMC3 human microglial cells were purchased from the Korean Cell Line Bank (Seoul, Korea). The cells were cultured in minimum essential medium (MEM) supplemented with 10% fetal bovine serum (FBS), 100 IU/ml penicillin, and 10 μg/ml streptomycin and were maintained in a humidified incubator at 37 °C with 95% air/5% CO2. The cells were treated with 100 ng/ml LPS (Sigma-Aldrich, St. Louis, MO, USA) and/or 500 nM JQ1 for 4 h under standard culture conditions. The LPS and JQ1 were dissolved in dimethyl sulfoxide (DMSO; Sigma-Aldrich, St. Louis, MO, USA).

Quantitative RT-PCR

Total RNA extractions and cDNA preparation were performed according to the manufacture’s instruction (Takara, Shiga, Japan). Quantitative Reverse Transcription PCR (qRT-PCR) was performed using an ABI 7500 real-time PCR system (Applied Biosystems Inc., Foster City, CA, USA). The critical threshold (CT) value was normalized by the expression of an internal control, glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Finally, the results were also analyzed using the comparative critical threshold (CT) method. The primers were designed using Primer Bank (http://pga.mgh.harvard.edu/primerbank/index.html) and are listed in Supplementary Table 4.

Protein–protein interaction (PPI) network analysis

The Search Tool for the Retrieval of Interacting Genes (STRING, http://string.embl.de/) [44] was used to construct the PPI network for DEmRNAs (minimum required interaction score > 0.7). The interaction relationships of the proteins encoded by DEmRNAs were searched by STRING online software, and the combined score > 0.7 was used as the cut-off criterion. Cytoscape MCODE plug-in (Version 3.4.0, available online: http://www.cytoscape.org/) was applied for visualization of the protein–protein interaction. Additionally, the network analyzer was used to compute the basic properties of the PPI network, including average clustering co-efficient distribution, closeness centrality, average neighborhood connectivity, node degree distribution, shortest path length distribution, and topological coefficients.

Statistical analysis

All data are expressed as the mean ± standard deviation of the mean (SD). The chi-square test was performed to confirm whether there was a statistically significant relationship between categorical variables by JQ1 treatment. A p-value or padj ≤ 0.05 was considered significant. The statistical analyses were performed using IBM SPSS Statistics ver. 26.0 (IBM Corporation, Armonk, NY, USA). All qRT-PCR data were tested using one-way ANOVA followed by Tukey's honestly significant difference (HSD) post hoc test. Differences for which p ≤ 0.05 were considered significant.

Availability of data and materials

The RNA-seq data used and analyzed during the current study are included in this published article and its supplementary files. The open RNA-seq data (Accession number: GSE155408) is available in the NCBI database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE155408).

Abbreviations

BET:

Bromodomain and extra-terminal

BP:

Biological process

BRD:

Bromodomain-containing protein

CC:

Cellular component

DAVID:

Database for annotation, visualization and integrated discovery

DEmRNA:

Differentially expressed mRNA

DElncRNA:

Differentially expressed lncRNA

FDR:

False discovery rate

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

LncRNA:

Long non-coding RNA

Log2 FC:

Log2 fold change

LPS:

Lipopolysaccharide

MF:

Molecular function

PPI:

Protein–protein interaction

qRT-PCR:

Quantitative reverse transcription-polymerase chain reaction

RNA-seq:

RNA sequencing

STRING:

Search tool for the retrieval of interacting genes

WGCNA:

Weighted gene co-expression network analysis

References

  1. Croese T, Castellani G, Schwartz M. Immune cell compartmentalization for brain surveillance and protection. Nat Immunol. 2021;22(9):1083–92.

    Article  CAS  Google Scholar 

  2. Bennett ML, Bennett FC. The influence of environment and origin on brain resident macrophages and implications for therapy. Nat Neurosci. 2020;23(2):157–66.

    Article  CAS  Google Scholar 

  3. Bernier LP, York EM, MacVicar BA. Immunometabolism in the brain: how metabolism shapes microglial function. Trends Neurosci. 2020;43(11):854–69.

    Article  CAS  Google Scholar 

  4. Chen S, Yang J, Wei Y, Wei X. Epigenetic regulation of macrophages: from homeostasis maintenance to host defense. Cell Mol Immunol. 2020;17(1):36–49.

    Article  CAS  Google Scholar 

  5. Yeh H, Ikezu T. Transcriptional and Epigenetic Regulation of Microglia in Health and Disease. Trends Mol Med. 2019;25(2):96–111.

    Article  CAS  Google Scholar 

  6. Wang N, Wu R, Tang D, Kang R. The BET family in immunity and disease. Signal Transduct Target Ther. 2021;6(1):23.

    Article  Google Scholar 

  7. Nicodeme E, Jeffrey KL, Schaefer U, Beinke S, Dewell S, Chung CW, Chandwani R, Marazzi I, Wilson P, Coste H, et al. Suppression of inflammation by a synthetic histone mimic. Nature. 2010;468(7327):1119–23.

    Article  CAS  Google Scholar 

  8. Das A, Chai JC, Yang CS, Lee YS, Das ND, Jung KH, Chai YG. Dual transcriptome sequencing reveals resistance of TLR4 ligand-activated bone marrow-derived macrophages to inflammation mediated by the BET inhibitor JQ1. Sci Rep. 2015;5:16932.

    Article  CAS  Google Scholar 

  9. Jung KH, Das A, Chai JC, Kim SH, Morya N, Park KS, Lee YS, Chai YG. RNA sequencing reveals distinct mechanisms underlying BET inhibitor JQ1-mediated modulation of the LPS-induced activation of BV-2 microglial cells. J Neuroinflammation. 2015;12:36.

    Article  Google Scholar 

  10. Baek M, Yoo E, Choi HI, An GY, Chai JC, Lee YS, Jung KH, Chai YG. The BET inhibitor attenuates the inflammatory response and cell migration in human microglial HMC3 cell line. Sci Rep. 2021;11(1):8828.

    Article  CAS  Google Scholar 

  11. Carpenter S, Aiello D, Atianand MK, Ricci EP, Gandhi P, Hall LL, Byron M, Monks B, Henry-Bezy M, Lawrence JB, et al. A long noncoding RNA mediates both activation and repression of immune response genes. Science. 2013;341(6147):789–92.

    Article  CAS  Google Scholar 

  12. Krawczyk M, Emerson BM. p50-associated COX-2 extragenic RNA (PACER) activates COX-2 gene expression by occluding repressive NF-kappaB complexes. Elife. 2014;3: e01776.

    Article  Google Scholar 

  13. Li Z, Chao TC, Chang KY, Lin N, Patil VS, Shimizu C, Head SR, Burns JC, Rana TM. The long noncoding RNA THRIL regulates TNFalpha expression through its interaction with hnRNPL. Proc Natl Acad Sci U S A. 2014;111(3):1002–7.

    Article  CAS  Google Scholar 

  14. Palazzo AF, Koonin EV. Functional Long Non-coding RNAs Evolve from Junk Transcripts. Cell. 2020;183(5):1151–61.

    Article  CAS  Google Scholar 

  15. Statello L, Guo CJ, Chen LL, Huarte M. Gene regulation by long non-coding RNAs and its biological functions. Nat Rev Mol Cell Biol. 2021;22(2):96–118.

    Article  CAS  Google Scholar 

  16. Atianand MK, Caffrey DR, Fitzgerald KA. Immunobiology of Long Noncoding RNAs. Annu Rev Immunol. 2017;35:177–98.

    Article  CAS  Google Scholar 

  17. Chen YG, Satpathy AT, Chang HY. Gene regulation in the immune system by long noncoding RNAs. Nat Immunol. 2017;18(9):962–72.

    Article  CAS  Google Scholar 

  18. Chew CL, Conos SA, Unal B, Tergaonkar V. Noncoding RNAs: master regulators of inflammatory signaling. Trends Mol Med. 2018;24(1):66–84.

    Article  CAS  Google Scholar 

  19. Rapicavoli NA, Qu K, Zhang J, Mikhail M, Laberge RM, Chang HY. A mammalian pseudogene lncRNA at the interface of inflammation and anti-inflammatory therapeutics. Elife. 2013;2:e00762.

    Article  Google Scholar 

  20. Hirose T, Virnicchi G, Tanigawa A, Naganuma T, Li R, Kimura H, Yokoi T, Nakagawa S, Benard M, Fox AH, et al. NEAT1 long noncoding RNA regulates transcription via protein sequestration within subnuclear bodies. Mol Biol Cell. 2014;25(1):169–83.

    Article  Google Scholar 

  21. Chan J, Atianand M, Jiang Z, Carpenter S, Aiello D, Elling R, Fitzgerald KA, Caffrey DR. Cutting Edge: A Natural Antisense Transcript, AS-IL1alpha, Controls Inducible Transcription of the Proinflammatory Cytokine IL-1alpha. J Immunol. 2015;195(4):1359–63.

    Article  CAS  Google Scholar 

  22. Lu Y, Liu X, Xie M, Liu M, Ye M, Li M, Chen XM, Li X, Zhou R. The NF-kappaB-Responsive Long Noncoding RNA FIRRE Regulates Posttranscriptional Regulation of Inflammatory Gene Expression through Interacting with hnRNPU. J Immunol. 2017;199(10):3571–82.

    Article  CAS  Google Scholar 

  23. Rousselet GA, Pernet CR. Improving standards in brain-behavior correlation analyses. Front Hum Neurosci. 2012;6:119.

    Article  Google Scholar 

  24. Liu T, Zhang L, Joo D, Sun SC. NF-kappaB signaling in inflammation. Signal Transduct Target Ther. 2017;2:17023.

  25. Kopitar-Jerala N. The role of interferons in inflammation and inflammasome activation. Front Immunol. 2017;8:873.

    Article  Google Scholar 

  26. Han X, Shi H, Sun Y, Shang C, Luan T, Wang D, Ba X, Zeng X. CXCR2 expression on granulocyte and macrophage progenitors under tumor conditions contributes to mo-MDSC generation via SAP18/ERK/STAT3. Cell Death Dis. 2019;10(8):598.

    Article  CAS  Google Scholar 

  27. Gusyatiner O, Bady P, Pham MDT, Lei Y, Park J, Daniel RT, Delorenzi M, Hegi ME. BET inhibitors repress expression of interferon-stimulated genes and synergize with HDAC inhibitors in glioblastoma. Neuro Oncol. 2021;23(10):1680–92.

    Article  CAS  Google Scholar 

  28. Soto-Diaz K, Juda MB, Blackmore S, Walsh C, Steelman AJ. TAK1 inhibition in mouse astrocyte cultures ameliorates cytokine-induced chemokine production and neutrophil migration. J Neurochem. 2020;152(6):697–709.

    Article  CAS  Google Scholar 

  29. Huang J, Li J, Li Y, Lu Z, Che Y, Mao S, Lei Y, Zang R, Zheng S, Liu C, et al. Interferon-inducible lncRNA IRF1-AS represses esophageal squamous cell carcinoma by promoting interferon response. Cancer Lett. 2019;459:86–99.

    Article  CAS  Google Scholar 

  30. Zhang H, Xue C, Wang Y, Shi J, Zhang X, Li W, Nunez S, Foulkes AS, Lin J, Hinkle CC, et al. Deep RNA Sequencing Uncovers a Repertoire of Human Macrophage Long Intergenic Noncoding RNAs Modulated by Macrophage Activation and Associated With Cardiometabolic Diseases. J Am Heart Assoc. 2017;6(11):e007431.

    Article  Google Scholar 

  31. Thompson RC, Vardinogiannis I, Gilmore TD. Identification of an NF-kappaB p50/p65-responsive site in the human MIR155HG promoter. BMC Mol Biol. 2013;14:24.

    Article  CAS  Google Scholar 

  32. Hadjicharalambous MR, Roux BT, Feghali-Bostwick CA, Murray LA, Clarke DL, Lindsay MA. Long Non-coding RNAs Are Central Regulators of the IL-1beta-Induced Inflammatory Response in Normal and Idiopathic Pulmonary Lung Fibroblasts. Front Immunol. 2018;9:2906.

    Article  CAS  Google Scholar 

  33. Wang D, Gu C, Liu M, Liu G, Liu H, Wang Y. Analysis of Long Noncoding RNA Expression Profile in Human Pulmonary Microvascular Endothelial Cells Exposed to Lipopolysaccharide. Cell Physiol Biochem. 2019;52(4):653–67.

    Article  CAS  Google Scholar 

  34. Roux BT, Heward JA, Donnelly LE, Jones SW, Lindsay MA. Catalog of Differentially Expressed Long Non-Coding RNA following Activation of Human and Mouse Innate Immune Response. Front Immunol. 2017;8:1038.

    Article  Google Scholar 

  35. Baek M, Chai JC, Choi HI, Yoo E, Binas B, Lee YS, Jung KH, Chai YG. Comprehensive transcriptome profiling of BET inhibitor-treated HepG2 cells. PLoS One. 2022;17(4):e0266966.

    Article  CAS  Google Scholar 

  36. Choi HI, An GY, Yoo E, Baek M, Binas B, Chai JC, Lee YS, Jung KH, Chai YG. The bromodomain inhibitor JQ1 up-regulates the long non-coding RNA MALAT1 in cultured human hepatic carcinoma cells. Sci Rep. 2022;12(1):7779.

    Article  CAS  Google Scholar 

  37. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  CAS  Google Scholar 

  38. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  Google Scholar 

  39. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  Google Scholar 

  40. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  Google Scholar 

  41. da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  42. Bu D, Luo H, Huo P, Wang Z, Zhang S, He Z, Wu Y, Zhao L, Liu J, Guo J, et al. KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 2021;49(W1):W317–25.

    Article  CAS  Google Scholar 

  43. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  Google Scholar 

  44. Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S, Doncheva NT, Legeay M, Fang T, Bork P, et al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49(D1):D605–12.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Research Foundation of Korea (NRF) Grants 2017M3A9G7073033 and 2020R1A2C1014193 (to Y. G. C.), 2016R1D1A1B04934970 (to K. H. J.), and 2014M3C9A3064693 (to Y. S. L.), 2020R1F1A1063217 (to B. B.) from the Korean government.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: Mina Baek, Kyoung Hwa Jung, Young Gyu Chai. Data curation: Jin Choul Chai, Hae In Choi, Bert Binas, Young Seek Lee, Kyoung Hwa Jung. Funding acquisition: Bert Binas, Young Seek Lee, Kyoung Hwa Jung, Young Gyu Chai. Investigation and validation: Hae In Choi, Eunyoung Yoo. Methodology and Software: Jin Choul Chai. Project administration: Mina Baek. Supervision: Young Gyu Chai. Writing: Mina Baek, Bert Binas. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Young Seek Lee, Kyoung Hwa Jung or Young Gyu Chai.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have declared that no competing interests exist.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Fig. 1.

Identification and functional annotation of DEmRNAs.

Additional file 2: Supplementary Fig. 2.

qRT-PCR validation of the DEmRNAs.

Additional file 3: Supplementary Table 1.

All significantly up- and downregulated DElncRNAs and DEmRNAs.

Additional file 4: Supplementary Table 2.

Pearson correlation coefficient (r) values of DElncRNA-DEmRNA pairs.

Additional file 5: Supplementary Table 3.

Results of the chi-square test for each module.

Additional file 6: Supplementary Table 4.

List of primers used for qRT-PCR.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Baek, M., Chai, J.C., Choi, H.I. et al. Analysis of differentially expressed long non-coding RNAs in LPS-induced human HMC3 microglial cells. BMC Genomics 23, 853 (2022). https://doi.org/10.1186/s12864-022-09083-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-09083-6

Keywords