Transcriptome sequencing of microglial cells stimulated with TLR3 and TLR4 ligands
© Das et al. 2015
Received: 1 December 2014
Accepted: 26 June 2015
Published: 10 July 2015
Resident macrophages in the CNS microglia become activated and produce proinflammatory molecules upon encountering bacteria or viruses. TLRs are a phylogenetically conserved diverse family of sensors that drive innate immune responses following interactions with PAMPs. TLR3 and TLR4 recognize viral dsRNA Poly (I:C) and bacterial endotoxin LPS, respectively. Importantly, these receptors differ in their downstream adaptor molecules. Thus far, only a few studies have investigated the effects of TLR3 and TLR4 in macrophages. However, a genome-wide search for the effects of these TLRs has not been performed in microglia using RNA-seq. Gene expression patterns were determined for the BV-2 microglial cell line when stimulated with viral dsRNA Poly (I:C) or bacterial endotoxin LPS to identify novel transcribed genes, as well as investigate how differences in downstream signaling could influence gene expression in innate immunity.
Sequencing assessment and quality evaluation revealed that common and unique patterns of proinflammatory genes were significantly up-regulated in response to TLR3 and TLR4 stimulation. However, the IFN/viral response gene showed a stronger response to TLR3 stimulation than to TLR4 stimulation. Unexpectedly, TLR3 and TLR4 stimulation did not activate IFN-ß and IRF3 in BV-2 microglia. Most importantly, we observed that previously unidentified transcription factors (TFs) (i.e., IRF1, IRF7, and IRF9) and the epigenetic regulators KDM4A and DNMT3L were significantly up-regulated in both TLR3- and TLR4-stimulated microglia. We also identified 29 previously unidentified genes that are important in immune regulation. In addition, we confirmed the expressions of key inflammatory genes as well as pro-inflammatory mediators in the supernatants were significantly induced in TLR3-and TLR4-stimulated primary microglial cells. Moreover, transcriptional start sites (TSSs) and isoforms, as well as differential promoter usage, revealed a complex pattern of transcriptional and post-transcriptional gene regulation upon infection with TLR3 and TLR4. Furthermore, TF motif analysis (-950 to +50 bp of the 5′ upstream promoters) revealed that the DNA sequences for NF-κB, IRF1, and STAT1 were significantly enriched in TLR3- and TLR4-stimulated microglia.
These unprecedented findings not only permit a comparison of TLR3-and TLR4-stimulated genes but also identify new genes that have not been previously implicated in innate immunity.
The immune system includes cells that can recognize Toll-like receptors (TLRs) and respond to TLR ligands by expressing several genes that are implicated in innate immunity, which is the first line of mammalian defense against invading pathogens. TLRs distinguish conserved patterns that are derived from microbial pathogens, which are recognized as pathogen-associated molecular patterns (PAMPs), including TLR3 (viral double-stranded RNA (dsRNA)), TLR7/8 (single-stranded RNA (ssRNA)) and TLR9 (DNA) . The interaction of TLRs with a PAMP, such as bacterial endotoxin lipopolysaccharide (LPS) and viral dsRNA, triggers several signaling cascades that lead to the induction of numerous target genes involved in inflammation . Polyinosinic-polycytidylic acid Poly (I:C) is a synthetic analog of viral dsRNA, which is a PAMP that is generated during the replication of RNA and DNA viruses and is recognized by TLR3, while TLR4 is the receptor that recognizes LPS . TLR3 and TLR4 are unique in that these receptors can transmit signals through the myeloid differentiation factor 88 (MyD88)-independent, TIR-domain–containing adaptor inducing interferon ß (TRIF)-dependent pathway . Next, MyD88 and TRIF induce the key transcription factor (TF) nuclear factor κB (NF-κB) signalosome and mitogen-activated protein kinases (MAPKs); both pathways facilitate the transcription and stabilization of mRNAs for proinflammatory mediators, including cyclooxygenase-2 (COX-2), interleukin-6 (IL-6), IL-1ß, tumor necrosis factor-α (TNF-α), and IL-12 . Although these signaling molecules and their cascades have been explored in macrophages, the dynamic outcome of genome-wide approaches and the effects of TLR3 vs TLR4 in murine microglia remain poorly understood.
Microglia, which are a type of glial cell, are resident macrophages of the brain and spinal cord that act as primary effector cells and these cells play an important role in the brain’s innate immunity, neuronal homeostasis, and neuroinflammatory pathologies . Microglial cells become rapidly activated in response to infection, inflammation or injury, and their activation is associated with the production and secretion of reactive oxygen species (ROS), nitric oxide (NO) and prostaglandin E2 (PGE2), and a variety of proinflammatory cytokines, including IL-1ß, IL-6 and TNF-α . Although microglial activation is considered a protective mechanism involved in the clearance of pathogen infection and in the regulation of tissue repair and recovery, excessive or persistent activation as an uncontrolled immune response stimulates and increases the production of neurotoxic pro-inflammatory mediators, causing neuroinflammation and neurological disorders including Alzheimer’s disease (AD), Parkinson’s disease (PD), brain ischemia and multiple sclerosis (MS) . Currently, the heterogeneity of microglia in terms of their activation phenotypes and their role in neurodegenerative disease progression and repair is an area of active investigation. Microglia express all known TLRs to detect and respond robustly to both Poly (I:C) and LPS . Previous studies have demonstrated that LPS stimulation induces the gene expression of IL-1ß, IL-6, TNF-α, inducible nitric oxide synthase (iNOS) and prostaglandin-endoperoxide synthase 2 (PTGS-2), as well as the production of NO and PGE2 in microglial cell cultures . In contrast, Poly (I:C) activates antiviral and inflammatory genes independent of the adaptor protein MyD88, which is required for all other TLRs . However, the full effect on gene up-regulation and the basis of this synergy is unknown.
Although a few studies have compared the effects of Poly (I:C) vs LPS on murine microglia or macrophages [8, 9], thus far, a genome-wide search for the effects of Poly (I:C) vs LPS using the RNA-seq specifically has not been performed in murine microglia. Therefore, we performed gene array studies and comparative gene expression profiling analyses of BV-2 cells treated with Poly (I:C) and LPS using RNA-seq. In this study we examined BV-2 cell lines as a model of inflammation studies. This is one of the major uses of microglia. Previously, others reports demonstrated that BV-2 cell lines have close resemblance to primary brain microglia [10–12]. Since BV-2 cells are easy to culture, they are an important tool to study not only inflammatory processes , but also phagocytosis . Increasing yield and homogeneity of cells in vitro culture allows for enlarged data output compared to most in vivo experiments. To the best of our knowledge, this study is the first to apply these approaches to assess the differences in response induced by different TLR ligands and their consequences with respect to global gene expression in BV-2 microglial cells.
Cell culture and stimulation
Mouse microglial BV-2 cells were grown in high-glucose Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10 % fetal bovine serum (FBS) (catalog # 26140), 100 IU/ml penicillin, and 10 μg/ml streptomycin (catalog # 15140; Invitrogen, USA). The cells were maintained in a humidified incubator with 95 % air and a 5 % CO2 atmosphere at 37 °C. The cells were incubated with LPS (10 ng/ml, Sigma-Aldrich) and Poly (I:C) (5 μg/ml, Sigma-Aldrich) for the specified times under normal culture conditions. The medium, which contained the appropriate agents, was replaced every other day. The NF-κB pathway Inhibitor BAY11-7082 was purchased from Calbiochem (San Diego, CA, USA). Unless otherwise indicated, BAY11-7082 was freshly dissolved before use in the experiments. In this study, we used 10 μM of BAY11-7082 [14–16]. Primary microglial cells were isolated from 3-day-old ICR mice as previously described [17, 18]. All experimental protocols were conducted in accordance with Institutional Animal Care and Use Committee (IACUC) guidelines and were approved by the IACUC committee at Hanyang University (HY-IACUC-2014-0164A). Briefly, whole brains of neonatal mice were taken; blood vessel and meninges were carefully removed. Then, the whole brains of 12 mice were pooled together, finely minced, and digested with Neural Tissue Dissociation Kit-Postnatal Neurons (Miltenyi Biotec-130-094-802). Next, digested cells pass through 70-μm nylon cell strainer (BD Biosciene) and were seeded in poly-L-lysine-coated T-75 flask in DMEM/nutrient mixture F-12 (DMEM/F12, 1:1) containing 20 % FBS (catalog # 26140), 100 IU/ml penicillin and 10 μg/ml streptomycin (catalog # 15140) from Invitrogen (CA, USA). The cells were maintained in a humidified incubator with a 95 % air/5 % CO2 atmosphere at 37 °C. The medium was changed every 2–3 days. After two weeks in culture, mixed glial cell cultures are shaken at 150 rpm at 37 °C for 45 min, and the glial cell suspension was collected from each flask and seeded on poly-L-lysine coated cell culture plate. Microglial cells were sub-plated and used for further experiments.
Total RNA extraction
Total RNA (~8 μg) was extracted using TRIzol (Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s instructions. Briefly, 200 μl of chloroform was added, and the tubes with the lysis mixture were gently inverted for 5 min. The mixture was centrifuged at 12,000 x g for 15 min at 4 °C, and the clear upper solution was placed into a new tube, to which 500 μl of isopropanol was added. The tubes were inverted before incubation on ice for 1 h. The lysis mixture was centrifuged at 12,000 x g for 10 min at 4 °C, and the isopropanol was decanted. Ice-cold 70 % ethanol was added to the RNA pellet for gentle washing. The ethanol was removed after centrifugation as indicated above for 10 min. The RNA pellets were dried at room temperature for 5–10 min before reconstitution in 20 μl of RNase-free water, and the RNA was treated with RNase-free DNase I (Promega, WI, USA). The RNA quality was assessed using an Agilent 2100 Bioanalyzer with an RNA 6000 Nano Chip (Agilent Technologies, Waldbronn, Germany), and the quantity was determined using a spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA).
Quantitative real-time RT-PCR (qRT-PCR)
List of primers used in qRT-PCR studies
Forward sequence (5′ - >3′)
Reverse sequence (5′ - >3′)
CAG GCG GTG CCT ATG TCT C
CGA TCA CCC CGA AGT TCA GTA G
GAA ATG CCA CCT TTT GAC AGT G
CTG GAT GCT CTC ATC AGG ACA
TGC TGG GTC TGA GTG GGA CT
CCC TAT GGC CCT CAT TCT CAC
ATG CCA ATC ACT CGA ATG CG
TTG TAT CGG CCT GTG TGA ATG
AGG CTT CTG GGC CTT ATG TG
TGC TTC TCT CGC CAG GAA TAC
GAC CAC ACT CTG CCC ACA C
TCC TGG GGT ATT TCC AGA CA
GAG AGC CGA ACG AGG TTC AG
CTT CCA GGT TGA CAC GTC CG
AGT GGA CAT TTT GCT CTC TCG
GTT AAT GAG GTC ACT GCG TTG A
cDNA library preparation for RNA-seq
Total RNA was extracted from each group of BV-2 cells i.e., control 2 h (2 samples), control 4 h (2 samples), LPS 2 h (2 samples), LPS 4 h (2 samples), Poly (I:C) 2 h (2 samples) and Poly (I:C) 4 h (2 samples), using TRIzol® (Life Technologies) according to the manufacturer’s instructions. For RNA-Seq, RNA libraries were created from each group of BV-2 cells using a NEBNext® Ultra™ Directional RNA Library Preparation Kit from Illumina®. The first step in the workflow involved the removal of ribosomal RNA using an RNAMius™ Transcriptome Isolation Kit (Life Technologies). Following purification, total RNA was fragmented into small pieces using divalent cations at elevated temperatures. The cleaved RNA fragments were subjected to first-strand cDNA synthesis using reverse transcriptase and random primers, followed by second-strand cDNA synthesis using DNA polymerase I and RNase H. The cDNA fragments were subsequently processed using an end-repair reaction after the addition of a single ‘A’ base, followed by adapter ligation. The products of these reactions were purified and enriched through PCR to generate the final cDNA library. The cDNA fragments were sequenced using an Illumina HiSeq2500 System (101 cycles, paired-end (PE) lane) (National Instrumentation Center for Environmental Management in Seoul National University). Biological duplicate RNA sequencing was performed on twelve independent RNA samples from the treated BV-2 cells.
Differential gene expression analysis
Raw sequence files were subjected to quality control analysis using FastQC software (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). To ensure high-quality data, we clipped the adapters and trimmed the reads using the FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/). The quality-checked reads for each condition were processed using the TopHat version 2.0.10 software package (Bowtie 2 version 2.2.1 software) as FASTQ files to analyze differentially expressed genes . The reads were mapped to the reference genome [Mus musculus UCSC mm10 sequence], and the alignment files were generated as BAM files. These files were used as the input for Cufflinks software, which is a complementary method used to generate assembled transcripts for each condition; the abundance was evaluated using read data. The “fragments per kilobase per million map reads” (FPKM) values were calculated for each gene to normalize the data . These assemblies are used with the Cuffquant and Cuffdiff tools from the Cufflinks 2.2.1 package to calculate the differential expression levels and to evaluate the statistical significance of the detected alterations . RNA-seq experiments were normalized and visualized using HOMER software (http://homer.salk.edu/homer/) after preparing custom tracks for the UCSC Genome browser (http://genome.ucsc.edu/).
DAVID (Database for Annotation, Visualization and Integrated Discovery) version 6.7 software (http://david.abcc.ncifcrf.gov/home.jsp) was used to determine the most functional annotation of significant genes in the datasets as previously described . The DAVID program calculates a modified Fisher’s exact P value to demonstrate gene ontology (GO) or molecular pathway enrichment. P values less than 0.05 were considered strongly enriched in the annotation category.
Canonical pathway analysis of datasets
Ingenuity Pathway Analysis (IPA) (Ingenuity Systems, http://www.ingenuity.com, Mountain View, CA, USA) was conducted to analyze the most significant canonical pathways in datasets as previously described [19, 24]. The genes from datasets associated with canonical pathways in the Ingenuity Pathways Knowledge Base (IPAKB) were considered for literary analysis. The significance of the associations between datasets and canonical pathways was measured in the following manner: (1) the ratio of the number of genes from the dataset that mapped to a canonical pathway divided by the total number of genes that mapped to the same canonical pathway; and (2) Fisher’s exact test for a P value indicating the probability that the association could be explained by chance. After uploading the datasets, gene identifiers were mapped to corresponding gene objects, and the genes were overlaid onto a global molecular network in the IPAKB. Gene networks were algorithmically generated based on connectivity.
Graphical representation of networks and pathways
The molecules from the normalized filtered RNA-seq dataset were each mapped to corresponding objects in the IPAKB for network generation. A fold change cutoff of up-regulated genes (≥1.5 log2-fold) was set to identify significantly and differentially regulated genes in BV-2 microglial cells at 2 and 4 h after Poly (I:C) and LPS stimulation. The graphical representation of molecular relationships between genes and gene products is based on genes or gene products, which are represented as nodes, and the biological relationship between two nodes is represented as an edge (line). All edges were supported by at least one reference from the literature, textbook or canonical information in the IPAKB. The node color intensity indicates the degree of up-regulation (red). The nodes are displayed using shapes to represent functional classes of gene products.
Transcription factor binding motif enrichment analysis
NCBI reference sequence mRNA accession numbers were subjected to transcription factor binding motif analysis using the web-based software Pscan (http://184.108.40.206/pscan/) . The JASPAR database of transcription binding factor sequences was analyzed using enriched groups of -950 base pair (bp) sequences to +50 bp of the 5′ upstream promoters . The range of -950 to +50 was selected from the range options in Pscan software to obtain the best coverage for a -1000 to +50 bp range.
Enzyme-linked immunosorbent assay (ELISA)
Primary microglial cells were cultured in the same condition as above. Primary microglial cells were treated with Poly (I:C) (5 μg/ml) and (LPS 10 ng/mL), for 2 h and 4 h. After treatment, the concentration of the pro-inflammatory mediators DNMT3L, TNF-α, IL-1ß and CCL4 were determined in cell culture supernatants using the mouse ELISA kit according to the manufacturer’s protocol. The mouse ELISA kit DNMT3L from Elabscience (Wuhan, China) and TNF-α, IL-1ß and CCL4 from Komabiotec (Seoul, Korea).
The data were analyzed using Origin Pro 8 software (Origin Lab Corporation, Northampton, MA, USA). Each value is expressed as the mean ± standard error of the mean (SEM). The statistical analysis was performed using SPSS 17.0 software (SPSS Inc., IL, USA). The data were tested using a one-way ANOVA, followed by Tukey’s HSD post-hoc test. * P <0.05 and ** P <0.001 were considered significant.
RNA-Seq transcriptional profiles of murine microglia in response to TLR3 vs TLR4
Gene network analysis and canonical pathways modulated through TLR3- and TLR4-stimulated BV-2 microglial cells
Differential induction of proinflammatory cytokines, chemokines and interferon response genes by TLR3 and TLR4 stimulation
Different TFs in multiple families of RNA-Seq data are activated by TLR3 and TLR4 stimulation in BV-2 microglial cells
Leads to activation of inflammatory genes by identified TFs in response to TLR3 stimulation
STAT1 predicted to be activated (45 genes)p = 4.76 e-43)
IRF7 predicted to be activated (45 genes) (p = 1.33 e-43)
NFKBIA predicted to be activated (35 genes) (p = 5.31 e-30)
IRF1 predicted to be activated (24 genes) (p = 7.19 e-22)
JUNB predicted to be activated (7 genes) (p = 1.24 e-07)
USP18, TRAFD1, TNF, TAP1, STAT2, SP110, SOCS3, SLFN5, SLFN2, SLFN13, SLFN12L, RSAD2, PTGS2, OASL, NEURL3, MX1/MX2, JAK2, ISG15, IRGM1, IRG1, IRF9, IRF7, IRF1, IL15, IGTP, IFIT2, IFIT1B, IFI35, IFI16, ICAM1, HERC6, EIF2AK2, CXCL10, CMPK2, CDKN1A, CD40, CD274, CD14, CCRL2, CCL2, CASP4, C3, APOBEC3B, CASP2, APOC2
USP18, UBA7, TRIM5, TRIM30A, TRIM21, TREX1, TAP1, STAT2, STAT1, SAMD9L, RTP4, RSAD2, PLSCR1, PHF11, PELI1, PARP14, PARP12, OASL2, OASL, OAS3, OAS2, OAS1, MX1/MX2, JAK2, ISG15, IRGM1, IRGM, IRF9, IRF1, IL15, IGTP, IFIT2, IFIT1B, IFIH1, IFI35, IFI16, HELZ2, DHX58, DAXX, CXCL10, CMPK2, CD40, CCRL2, CASP4, ADAR
TNFAIP3, TNFAIP2, SOCS3, SLFN2, SAA3, RCAN1, RASA3, PTGS2, PLAU, NFKBIE, NFKB2, JUNB, JUN, ISG15, IRG1, IRF1, IL1RN, IL1A, IL15, IFIT1B, IFI16, ICAM1, HMOX1, HLA-A, GCH1, CLU, CEBPD, CCL9, CCL7, CCL3L3, CCL2, CASP4, AMPD3, MXD4, CXCR4
TRIM21, TNF, TLR3, TAP1, STAT2, STAT1, RSAD2, PTGS2, PML, OAS1, JAK2, ISG15, IRF9, IRF7, IL1B, IL18, IL15, IFIT2, IFIH1, IFI35, EIF2AK2, CXCL10, CDKN1A, CD40, CASP2
TNFRSF8, SERPINE1, PLAUR, PLAU, HMOX1, BCL3, ATF3
Leads to activation of inflammatory genes by identified TFs in response to TLR4 stimulation
IRF7 predicted to be activated (43 genes) (p = 3.25 e-38)
STAT1 predicted to be activated (36 genes) (p = 3.02 e-32)
NFKBIA predicted to be activated (35 genes) (p = 5.88 e-25)
IRF1 predicted to be activated (20 genes) (p = 1.10 e-14)
CEBPD predicted to be activated (11 genes) (p = 1.22 e-10)
USP18, XAF1, UBA7, TRIM5, TRIM30A, TRIM21, TREX1, TAP2, TAP1, STAT2, STAT1, SAP30, RSAD2, PHF11, PELI1, PARP14, PARP12, OASL2, OASL, OAS3. OAS2, OAS1, MX1/MX2, ISG15, IRGM1, IRGM, IRF9, IRF1, IRF7, IGTP, IFIT2, IFIT1B, IFIH1, IFI16, HELZ2, DHX58, DAXX, CXCL10, CMPK2, CCRL2, CASP4, ADAR, TLR8
USP18, TRAFD1, TNF, TAP1, STAT2, STAT1, SOCS3, SLFN5, SLFN2, SLFN13, SLFN12L, RSAD2, PTGS2, OASL, NEURL3, MX1/MX2, ISG15, IRGM1, IRG1, IRF9, IRF1, IRF7, IGTP, IFIT2, IFIT1B, IFI16, ICAM1, HERC6, EIF2AK2, CXCL10, CMPK2, CDKN1A, CD274, CD14, CCRL2, CASP4, C3, APOBEC3B,
TNFRSF1B, TNFAIP3, TNFAIP2, SOCS3, SLFN2, SAA3, RASA3, PTGS2, NFKBIE, NFKB2, NFKB1, JUNB, ISG15, IRG1, IRF1, IL1RN, IL1A, IFIT1B, IFI16, ICAM1, HMOX1, HLA-A, GCH1, CP, CFLAR, CEBPD, CEBPB, CCL9, CCL7, CCL2, CCL3L3, CASP4, BTG2, AMPD3, CXCR4
TRIM21, TNF, TLR3, TAP1, TAP2, STAT2, STAT1, RSAD2, PTGS2, OAS1, ISG15, IRF9, IRF7, IL1B, IL18, IFIT2, IFIH1, EIF2AK2, CXCL10, CDKN1A,
TNF, SAA3, PTGS2, PAX3, IKBKE, HP, CEBPB, CD14, CCL3, CCL2, C3
Comparison of the transcriptional and post-transcriptional regulation by TLR3 and TLR4 stimulation in BV-2 microglial cells
Different epigenetic regulators in multiple families of RNA-Seq data are activated by TLR3 and TLR4 stimulation in BV-2 microglial cells
Identification of novel genes synergistically up-regulated by TLR3 and TLR4 stimulation in BV-2 microglial cells
Many studies have investigated the effects of TLR ligands, particularly TLR4, on cytokines, chemokines, and interferon response genes, as well as of different TF families on macrophages, which are important in inflammatory diseases. However, none of these studies has focused on novel genes that are also potentially involved in inflammatory diseases. To gain greater insight into mRNA expression studies, our RNA-seq analysis unexpectedly revealed approximately 29 previously unidentified genes that were particularly up-regulated by TLR ligands. Surprisingly, these 29 novel genes are frequently overexpressed not only in Poly (I:C)-induced microglial responses but also in LPS-induced microglial responses (Fig. 7c and d). Notably, the up-regulated genes that exhibited the significant links with inflammation included CLEC4E, RSAD2, PDE4B, GPR84, IER3, SLC15A3, CASP4, FCGR4, PELI1, PIK2, RNF114 and TREX, which are known to play important roles in inflammatory disorders [33, 40–42]. We observed that the RSAD2, FCGR4, and RNF114 genes displayed higher up-regulation upon Poly (I:C) stimulation than upon LPS stimulation. The reverse situation was true for PDE4B, IER3, SLC15A3, PIK2 and TREX.
Confirmation of differentially expressed genes through qRT-PCR
Comparison of fold change values from RNA-seq data and qRT-PCR in 2 h LPS and Poly I:C treated BV-2 microglia cells
RNA-seq fold change
qRT-PCR fold change
Gene accession ID
Comparison of fold change values from RNA-seq data and qRT-PCR in 4 h LPS and Poly (I:C) treated BV-2 microglia cells
RNA-seq fold change
qRT-PCR fold change
Gene accession ID
Microglial cells become rapidly activated through interactions with pathogens, and their persistent activation is associated with the production and secretion of various pro-inflammatory genes, cytokines and chemokines, which may initiate or amplify neurodegenerative diseases . Previous studies have reported that LPS stimulation induces the gene expression of TNF-α, IL-1β, IL-6, iNOS and PTGS-2 and the production of NO and PGE2 in primary and BV-2 microglial cell cultures . In contrast, Poly (I:C) activates antiviral and inflammatory genes independent of the adaptor protein MyD88, which is required for all other TLRs . However, none of these studies addressed the effects of TLR ligands Poly (I:C) and LPS; thus, their responses were compared using RNA-seq analysis in microglial cells in the present study. This study provides the most comprehensive analysis thus far because, compared to microarrays, this technique can provide unbiased profiles, can be extremely accurate, and can identify novel transcribed regions if a sufficient level of coverage is obtained [43, 44]. Furthermore, these technologies facilitate the determination of the difference between the expression of alternative mature mRNAs from the same precursor and the identification of the differential expression of mRNA isoforms [43–45]. Validation techniques such as qRT-PCR  have corroborated the accuracy of RNA-seq.
The RNA-seq analysis revealed that cytokines/chemokines, antiviral genes, and IRGs associated with inflammation were significantly up-regulated in response to TLR3- and TLR4-induced microglial activation. Both the extent of the fold change and the number of genes were significantly modulated. The activation of TNF-α, IL-1β, IL-1A, PTGS-2, CXCL10 (IP-10), CCL7, CCL12, CCL4, and CCL2, all occurred similarly, with minor variations (Fig. 4a and b). IL-1 is the most widely studied pro-inflammatory gene; the extensively characterized forms of IL-1 are IL-1A and IL-1β. IL-1A and IL-1β play a crucial role in the development of AD and PD, the pathogenic hallmark of which is CNS inflammation. Following CNS damage, IL-1 is rapidly released from activated microglia, and an elevated level of the IL-1 cytokine is an important hallmark of neuroinflammation . PTGS-2 is the key enzyme responsible for brain inflammation, and increased PTGS-2 expression contributes to neurodegeneration . Chemokines are key regulators of inflammation, and the excessive production of these molecules has been associated with disease progression and severe inflammation pathologies, including MS. Conductier et. al.(2010), reported that CCL2 plays a crucial role in neuroinflammatory diseases and is considered a target in the treatment of neuroinflammatory disorders . CCL2 and CCL7 are highly expressed in microglia, astrocytes and other inflammatory cells during MS . CCL12 also plays an inflammatory role because the levels of this chemokine are up-regulated in both microglia and astrocytes when stimulated with the proinflammatory cytokine IL-17 . The expression of CXCL10 is observed during infectious and inflammatory diseases, playing a crucial role in T-cell-mediated inflammation in the CNS. In addition, CXCL10 plays a role in inflammatory demyelinating diseases, such as MS, through the destruction of the myelin sheath or neurons by facilitating leukocyte trafficking in the brain . These patterns are consistent with our studies showing the consequences of TLR3 and TLR4 stimulation with respect to the induction of cytokines/chemokines (Fig. 4a and b).
Various signaling pathways, including phosphoinositide 3-kinase, Jun N-terminal kinase (JNK), p38, NF-κB, extracellular signal-related kinase, and IRF3, have been activated by TLR3/TLR4, leading to the induction of numerous target genes involved in antiviral immunity, including IFN-ß and IP-10 production . Our RNA-seq data revealed that the induction of several genes involved in antiviral immunity and signaling was more potent in TLR3-stimulated BV-2 microglial cells than in TLR4-stimulated BV-2 microglial cells and included IFIT1, ISG15, IFIT2, IFI203, IFI204, IRF1, IRF7, IRF9, STAT1, and STAT2 (Fig. 4a). These results were rather anticipated because TLR3 is unique, recruiting only TRIF and not MyD88. Unexpectedly, we were unable to identify IRF3 and its target gene IFN-ß, but not IP-10, in TLR3/TLR4-stimulated BV-2 microglial cells. Our data show consequences similar to that in a previously published report regarding the induction of IRF3 in DCs and MØs . However, because no IRF3 activation could be observed in BV-2 microglia cells, the mechanism by which the production of IP-10 and IRGs is regulated in these cell types remains unclear. This mechanism is the subject of ongoing investigations. Furthermore, to evaluate the expressions of IRF3, IFN-ß, and IP-10 (CXCL10) on TLR3- and TLR4-induced primary microglial cells we measured the expressions of IRF3, IFN-ß, and IP-10 for 2 h and 4 h time points. Interestingly, we found that IFN-ß, and IP-10 were significantly up-regulated in primary microglial cells (Additional file 3 A and B and Fig. 9a and b). Although, primary microglial cells with these factors did not induce the expression of IRF-3 (Additional file 3 A and B). Nevertheless further studies are warranted to assess genome wide transcriptomic analysis on TLR3-and TLR4-stimulated primary microglial cells. Importantly, we first identified IRF1, IRF7, and IRF9 as significantly up-regulated in response to TLR3- and TLR4-stimulated microglia (Fig. 5a). A previous study reported that IRF7 could be a master regulator of type-I IFN-mediated responses . In addition, IRF1 has been suggested to be a master transcription factor contributing to IRGs . Thus, exploring whether interferon response gene induction mechanisms not relying on IRF3 may exist in cells other than microglia or macrophages would be interesting. Therefore, IRF7, IRF1, and IRF9 most likely do not act in place of IRF3 in BV-2 microglia cells in response to TLR3 and TLR4 stimulation. This possibility is an exciting area that we are keenly pursuing further.
In addition to IRF TFs, we identified NF-κB and STAT, as well as additional TFs (PU.1, KLF7, JUNB, ATF3, and FOXP4) whose roles in microglia activation have not yet been well established (Fig. 5a and b). We observed the significant up-regulation of NF-κB transcription factor, which was induced by Poly (I:C) and LPS stimulation in microglia, with LPS being more potent than Poly (I:C). Similar to our findings, Reimer et. al. (2008) reported that Poly (I:C) and LPS induce the activation of NF-κB, followed by the release of TNF-α . JUNB also plays an important role in controlling activity-dependent gene programs that are critical for nervous system function . However, IRF2, IRF4, IRF6, IRF8, STAT6, KLF1, KLF2, KLF4, and KLF5 were unaffected by Poly (I:C) and LPS stimulation, suggesting the highly selective induction of TFs through TLR3 and TLR4 stimulation in BV-2 microglial cells. To delineate conserved transcription factor-binding motifs further, we performed TF motif analysis on TLR3- and TLR4-stimulated genes in BV-2 microglial cells. The core promoters of co-expressed genes (typically, regulatory regions within -1000 to +50 bp relative to the TSS) can be evaluated for overrepresented cis-regulatory elements after partitioning into suitable modules . Of the two ranges available in Pscan software that are closest to this region of interest (-950 to 50 and -1000 to 0), the -950 to +50 bp range was selected for the analyses. The promoters of differentially expressed genes revealed the enrichment of DNA sequences not only for NF-κB transcription factors but also for IRF1 and STAT1, as shown in Fig. 5c. These analyses converged the first insights into 3 TF binding motifs that may be involved in regulating subset-specific genes in Poly (I:C)- and LPS-stimulated BV-2 microglial cells. Next, we used IPA software to identify the target genes that were directly or indirectly activated by the identified TFs in response to TLR3 and TLR4 stimulation. Importantly, we found that the majority of the cytokines, chemokines whose expression is directly regulated by the identified TFs, including IRF1, IRF7, STAT1, JUNB, CREPD, and NF-κBIA, where JUNB and CREPD target genes were unique in response to TLR3 and TLR4, respectively (Fig. 5d and e and Tables 2 and 3).
We performed the first functional analysis of the transcripts, isoforms and TSSs to gain greater insight into the functional categories for Poly (I:C) and LPS-stimulated inflammatory genes in BV-2 microglial cells. Most of the differentially expressed genes in microglial following Poly (I:C) and LPS stimulation were expressed as several isoforms subjected to transcriptional/post-transcriptional regulation and/or to differential promoter usage. We classified these genes into three primary groups (genes with one isoform and one TSS; genes with more than one isoform and one TSS; and genes with more than one isoform and more than one TSS) (Fig. 6a and b). The first two groups included genes crucial for the innate immune response, which might be under stronger selection to prevent the emergence of new isoforms and/or post-transcriptional regulation. In addition, some genes belonged to the third group , suggesting that these genes could be subjected to positive selection by transcriptional and post-transcriptional regulation in Poly (I:C)- and LPS-stimulated BV-2 microglial cells. Interestingly, we found that spliced and both transcriptional and post-transcriptional regulated genes were greater in Poly (I:C)-stimulated BV-2 microglia than in LPS-stimulated BV-2 microglia after unspliced and transcriptionally regulated genes. The reverse situation was true for LPS-stimulated genes. However, further targeted studies are required to validate this regulation and to establish the potential effects of these genes during TLR3 and TLR4 stimulation in microglia.
Another interesting finding is the increased expression of epigenetic regulators, which involves chemical modifications of DNA cytosine residues and DNA-bound histone proteins without alterations in the DNA sequence; epigenetic regulation is promising as one of the major factors regulating gene expression in response to environmental stimuli . Recent studies have demonstrated that histone demethylase (KDM6B) and histone deacetylases (HDAC1, HDAC2, HDAC3, and HDAC7) potentially regulate proinflammatory gene expression in macrophages [39, 57, 58]. Recently, we showed that the histone demethylase KDM4A was significantly expressed in neuroectodermal stem cells and might play a role in tumorigenic development . Interestingly, among different epigenetic regulators, only histone demethylase KDM4A and DNA methyltransferase DNMT3L were strikingly differentially expressed in Poly (I:C)- and LPS-stimulated BV-2 microglial cells, with LPS being more potent than Poly (I:C), as determined by the RNA-Seq data in the present study (Fig. 7a and b). The mechanism by which KDM4A and DNMT3L become activated following TLR3 or TLR4 receptor stimulation remains unknown; determining how these epigenetic regulators, along with modified TFs, can regulate distinct set of inflammatory genes in microglial cells would be intriguing. We posit that the role of these genes in inflammatory diseases should now be assessed in animal models involving TLR-specific gene deletion or overexpression. However, histone demethylase (KDM6B) and histone deacetylases (HDAC1, HDAC2, HDAC3, and HDAC7), were not expressed in either Poly (I:C)- or LPS-stimulated BV-2 microglial cells.
One of the most striking features is that our RNA-seq analysis was the first to identify several notable differences in the pattern of previously unidentified genes activation induced by these distinct TLR ligands in BV-2 microglial cells. In particular, this technology allowed us to identify over 25 direct TLR3- or TLR4-stimulated genes in these cells. The direct TLR3- or TLR4-stimulated genes known to be important for activity-regulated processes in microglia include the macrophage-inducible C-type lectin (CLEC4E), which regulates immune responses to pathogens ; the secreted endo-lysosomal peptide transporter SLC15A3, which is highly expressed by bone marrow-derived dendritic cells after LPS stimulation and which plays a key role in regulating innate immune responses ; orphan G protein-coupled receptor 84 (GPR84), which is highly expressed in leukocytes, monocytes, and macrophages upon activation by LPS and which plays a critical role in immunological regulation ; and caspase-4, which plays an important role in inflammatory diseases ; (Fig. 7c and d). This finding indicates that these TLR3 or TLR4 targets also mediate the effects of the inflammatory response in mouse microglia. Furthermore, we confirmed the expression of key inflammation- and immunity-related genes as well as cytokines/chemokines in the supernatants were significantly induced in Poly (I:C) and LPS treated primary microglial cells including IL-1ß, IL1A, TNF-α, PTGS2, CCL3, CCL4, CCL7, CXCL10, IRF1, IRF7, JUNB, NF-κBIA, CLEC4E, GPR84, SLC15A3 and KDM4A (Fig. 9a, b, c and d). Future studies determining how these TLR3 or TLR4 target genes are expressed during the immune response in mouse microglia will be useful. This mechanism is the subject of ongoing investigations.
Overall, our RNA-seq data provide novel insights into the transcriptional landscape TLR3 and TLR4 ligands mimic certain aspects of viral infection by triggering key molecular events such as a rapid innate immune response during the early course of the inflammatory response. Regardless of certain boundaries in exactitude, TLR3- and TLR4-stimulated inflammatory gene expression profiles, TF and epigenetic regulator clusters, and microglial cis-regulatory elements predictions offer valuable information for future studies, such as potential gene targets for chromatin immunoprecipitation (ChIP)-seq assays. In the future, this model can be extended to include data from other high-dimensional surveys, such as microRNA, ChIP-seq, and proteomics, providing further insight into global gene regulation in TLR3- and TLR4-stimulated microglial cells.
In summary, using RNA-seq, our study is the first to identify family-wide DEGs of most known and unknown immune genes that characterize the microglial response to two TLR stimuli, Poly (I:C) and LPS and to demonstrate both common and unique patterns. In addition, our data argue that, while both TLR3 and TLR4 have been evolutionarily selected to induce antiviral gene expression, TLR3 seems to be even more specialized than TLR4 to initiate antiviral responses and that the reverse situation was true for the master TF NF-κB. Furthermore, we examined epigenetic regulators, different TFs, and TF-regulated genes, as well as the transcriptional and post-transcriptional regulation of genes based on their isoforms, TSSs and differential promoters in two TLR-stimulated BV-2 microglial cells. These unprecedented results together indicate that both known and unknown DEGs, TFs, and epigenetic regulators identified in the present study may provide new insight regarding the innate immunity to TLR3 and TLR4 in BV-2 microglial cells and may enhance our knowledge of neuroinflammation.
Availability of supporting data
The acquired data were deposited in the NCBI sequence read archive (http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=run_browser) under dataset accession no. SRR1554368, SRR1598823, SRR1554452, SRR1554456, SRR1554369, SRR1598824, SRR1554453, SRR1554457, SRR1692678, SRR1692736, SRR1692771, and SRR1692789.
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP) (2013R1A1A3011026 to K.H.J & 2011-0030049 to Y.G.C).
- O’Neill LA. When signaling pathways collide. positive and negative regulation of toll-like receptor signal transduction. Immunity. 2008;29(1):12–20.PubMedView ArticleGoogle Scholar
- Krishnan J, Selvarajoo K, Tsuchiya M, Lee G, Choi S. Toll-like receptor signal transduction. Exp Mol Med. 2007;39(4):421–38.PubMedView ArticleGoogle Scholar
- Takeda K, Akira S. Toll-like receptors in innate immunity. Int Immunol. 2005;17(1):1–14.PubMedView ArticleGoogle Scholar
- Yamamoto M, Yamazaki S, Uematsu S, Sato S, Hemmi H, Hoshino K, et al. Regulation of Toll/IL-1-receptor-mediated gene expression by the inducible nuclear protein IkappaBzeta. Nature. 2004;430(6996):218–22.PubMedView ArticleGoogle Scholar
- Graeber MB, Streit WJ. Microglia: biology and pathology. Acta Neuropathol. 2010;119(1):89–105.PubMedView ArticleGoogle Scholar
- Merson TD, Binder MD, Kilpatrick TJ. Role of cytokines as mediators and regulators of microglial activity in inflammatory demyelination of the CNS. Neuromolecular Med. 2010;12(2):99–132.PubMedView ArticleGoogle Scholar
- Carpentier PA, Duncan DS, Miller SD. Glial toll-like receptor signaling in central nervous system infection and autoimmunity. Brain Behav Immun. 2008;22(2):140–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee HJ, Kong PJ, Lee SH, Kwon OY, Chun WJ, Kim SS. Differences between lipopolysaccharide and double-stranded RNA in innate immune responses of BV2 microglial cells. Int J Neurosci. 2007;117(7):885–94.PubMedView ArticleGoogle Scholar
- Reimer T, Brcic M, Schweizer M, Jungi TW. poly (I:C) and LPS induce distinct IRF3 and NF-kappaB signaling during type-I IFN and TNF responses in human macrophages. J Leukoc Biol. 2008;83(5):1249–57.PubMedView ArticleGoogle Scholar
- Blasi E, Barluzzi R, Bocchini V, Mazzolla R, Bistoni F. Immortalization of murine microglial cells by a v-raf/v-myc carrying retrovirus. J Neuroimmunol. 1990;27(2-3):229–37.PubMedView ArticleGoogle Scholar
- Horvath RJ, Nutile-McMenemy N, Alkaitis MS, Deleo JA. Differential migration, LPS-induced cytokine, chemokine, and NO expression in immortalized BV-2 and HAPI cell lines and primary microglial cultures. J Neurochem. 2008;107(2):557–69.PubMed CentralPubMedView ArticleGoogle Scholar
- Henn A, Lund S, Hedtjarn M, Schrattenholz A, Porzgen P, Leist M. The suitability of BV2 cells as alternative model system for primary microglia cultures or for animal experiments examining brain inflammation. ALTEX. 2009;26(2):83–94.PubMedGoogle Scholar
- Hirt UA, Leist M. Rapid, noninflammatory and PS-dependent phagocytic clearance of necrotic cells. Cell Death Differ. 2003;10(10):1156–64.PubMedView ArticleGoogle Scholar
- Pierce JW, Schoenleber R, Jesmok G, Best J, Moore SA, Collins T, et al. Novel inhibitors of cytokine-induced IkappaBalpha phosphorylation and endothelial cell adhesion molecule expression show anti-inflammatory effects in vivo. J Biol Chem. 1997;272(34):21096–103.PubMedView ArticleGoogle Scholar
- Mendes Sdos S, Candi A, Vansteenbrugge M, Pignon MR, Bult H, Boudjeltia KZ, et al. Microarray analyses of the effects of NF-κB or PI3K pathway inhibitors on the LPS-induced gene expression profile in RAW264.7 cells: Synergistic effects of rapamycin on LPSinduced MMP9-overexpression. Cell Signal. 2009;21(7):1109–22.PubMedView ArticleGoogle Scholar
- Juliana C, Fernandes-Alnemri T, Wu J, Datta P, Solorzano L, Yu JW, et al. Anti-inflammatory compounds parthenolide and Bay 11-7082 are direct inhibitors of the inflammasome. J Biol Chem. 2010;285(13):9792–802.PubMed CentralPubMedView ArticleGoogle Scholar
- Witting A, Moller T. Microglia cell culture: a primer for the novice. Methods Mol Biol. 2011;758:49–66.PubMedView ArticleGoogle Scholar
- Jung KH, Das A, Chai JC, Kim SH, Morya N, Park KS, et al. 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Das A, Chai JC, Jung KH, Das ND, Kang SC, Lee YS, et al. JMJD2A attenuation affects cell cycle and tumourigenic inflammatory gene regulation in lipopolysaccharide stimulated neuroectodermal stem cells. Exp Cell Res. 2014;S0014-4827(14):00367-X.Google Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.PubMed CentralPubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.PubMedView ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.PubMed CentralPubMedView ArticleGoogle Scholar
- Rudra D, deRoos P, Chaudhry A, Niec RE, Arvey A, Samstein RM, et al. Transcription factor Foxp3 and its protein partners form a complex regulatory network. Nat Immunol. 2012;13(10):1010–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Rosowski EE, Lu D, Julien L, Rodda L, Gaiser RA, Jensen KD, et al. Strain-specific activation of the NF-kappaB pathway by GRA15, a novel Toxoplasma gondii dense granule protein. J Exp Med. 2011;208(1):195–212.PubMed CentralPubMedView ArticleGoogle Scholar
- Zambelli F, Pesole G, Pavesi G. Pscan: finding over-represented transcription factor binding site motifs in sequences from co-regulated or co-expressed genes. Nucleic Acids Res. 2009;37:W247–52.PubMed CentralPubMedView ArticleGoogle Scholar
- Portales-Casamar E, Thongjuea S, Kwon AT, Arenillas D, Zhao X, Valen E, et al. JASPAR 2010: the greatly expanded open-access database of transcription factor binding profiles. Nucleic Acids Res. 2010;38:D105–10.PubMed CentralPubMedView ArticleGoogle Scholar
- Lund S, Christensen KV, Hedtjarn M, Mortensen AL, Hagberg H, Falsig J, et al. The dynamics of the LPS triggered inflammatory response of murine microglia under different culture and in vivo conditions. J Neuroimmunol. 2006;180(1-2):71–87.PubMedView ArticleGoogle Scholar
- Thomas DM, Francescutti-Verbeem DM, Kuhn DM. Gene expression profile of activated microglia under conditions associated with dopamine neuronal damage. FASEB J. 2006;20(3):515–7.PubMedGoogle Scholar
- Kaushik DK, Mukhopadhyay R, Kumawat KL, Gupta M, Basu A. Therapeutic targeting of Kruppel-like factor 4 abrogates microglial activation. J Neuroinflammation. 2012;9:57.PubMed CentralPubMedView ArticleGoogle Scholar
- Fitzgerald DC, O’Brien K, Young A, Fonseca-Kelly Z, Rostami A, Gran B. Interferon regulatory factor (IRF) 3 is critical for the development of experimental autoimmune encephalomyelitis. J Neuroinflammation. 2014;11:130.PubMed CentralPubMedView ArticleGoogle Scholar
- Tanaka S, Ishii A, Ohtaki H, Shioda S, Yoshida T, Numazawa S. Activation of microglia induces symptoms of Parkinson’s disease in wild-type, but not in IL-1 knockout mice. J Neuroinflammation. 2013;10:143.PubMed CentralPubMedView ArticleGoogle Scholar
- Shoemaker ML, Smirnova NP, Bielefeldt-Ohmann H, Austin KJ, van Olphen A, Clapper JA, et al. Differential expression of the type I interferon pathway during persistent and transient bovine viral diarrhea virus infection. J Interferon Cytokine Res. 2009;29(1):23–35.PubMedView ArticleGoogle Scholar
- Sollberger G, Strittmatter GE, Kistowska M, French LE, Beer HD. Caspase-4 is required for activation of inflammasomes. J Immunol. 2012;188(4):1992–2000.PubMedView ArticleGoogle Scholar
- Re F, Strominger JL. IL-10 released by concomitant TLR2 stimulation blocks the induction of a subset of Th1 cytokines that are specifically induced by TLR4 or TLR3 in human dendritic cells. J Immunol. 2004;173(12):7548–55.PubMedView ArticleGoogle Scholar
- Dijkstra IM, de Haas AH, Brouwer N, Boddeke HW, Biber K. Challenge with innate and protein antigens induces CCR7 expression by microglia in vitro and in vivo. Glia. 2006;54(8):861–72.PubMedView ArticleGoogle Scholar
- De Jesus AF, Chen G, Li Z, Grayfer L, Robert J. Susceptibility of Xenopus laevis tadpoles to infection by the ranavirus Frog-Virus 3 correlates with a reduced and delayed innate immune response in comparison with adult frogs. Virology. 2012;432(2):435–43.View ArticleGoogle Scholar
- Lundberg AM, Drexler SK, Monaco C, Williams LM, Sacre SM, Feldmann M, et al. Key differences in TLR3/poly I:C signaling and cytokine induction by human primary cells: a phenomenon absent from murine cell systems. Blood. 2007;110(9):3245–52.PubMedView ArticleGoogle Scholar
- Andreoli F, Barbosa AJ, Parenti MD, Del Rio A. Modulation of epigenetic targets for anticancer therapy: clinicopathological relevance, structural data and drug discovery perspectives. Curr Pharm Des. 2013;19(4):578–613.PubMedView ArticleGoogle Scholar
- Jeong Y, Du R, Zhu X, Yin S, Wang J, Cui H, et al. Histone deacetylase isoforms regulate innate immune responses by deacetylating mitogen-activated protein kinase phosphatase-1. J Leukoc Biol. 2014;95(4):651–9.PubMedView ArticleGoogle Scholar
- Yamasaki S, Ishikawa E, Sakuma M, Hara H, Ogata K, Saito T. Mincle is an ITAM-coupled activating receptor that senses damaged cells. Nat Immunol. 2008;9(10):1179–88.PubMedView ArticleGoogle Scholar
- Wang J, Wu X, Simonavicius N, Tian H, Ling L. Medium-chain fatty acids as ligands for orphan G protein-coupled receptor GPR84. J Biol Chem. 2006;281(45):34457–64.PubMedView ArticleGoogle Scholar
- Nakamura N, Lill JR, Phung Q, Jiang Z, Bakalarski C, de Maziere A, et al. Endosomes are specialized platforms for bacterial sensing and NOD2 signalling. Nature. 2014;509(7499):240–4.PubMedView ArticleGoogle Scholar
- Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, et al. A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008;321(5891):956–60.PubMedView ArticleGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, et al. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008;320(5881):1344–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Richard H, Schulz MH, Sultan M, Nurnberger A, Schrinner S, Balzereit D, et al. Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments. Nucleic Acids Res. 2010;38(10), e112.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilhelm BT, Marguerat S, Watt S, Schubert F, Wood V, Goodhead I, et al. Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution. Nature. 2008;453(199):1239–43.PubMedView ArticleGoogle Scholar
- Kawano T, Anrather J, Zhou P, Park L, Wang G, Frys KA, et al. Prostaglandin E2 EP1 receptors: downstream effectors of COX-2 neurotoxicity. Nat Med. 2006;12(2):225–9.PubMedView ArticleGoogle Scholar
- Conductier G, Blondeau N, Guyon A, Nahon JL, Rovere C. The role of monocyte chemoattractant protein MCP1/CCL2 in neuroinflammatory diseases. J Neuroimmunol. 2010;224(1-2):93–100.PubMedView ArticleGoogle Scholar
- Banisor I, Leist TP, Kalman B. Involvement of beta-chemokines in the development of inflammatory demyelination. J Neuroinflammation. 2005;2(1):7.PubMed CentralPubMedView ArticleGoogle Scholar
- Das Sarma J, Ciric B, Marek R, Sadhukhan S, Caruso ML, Shafagh J, et al. Functional interleukin-17 receptor A is expressed in central nervous system glia and upregulated in experimental autoimmune encephalomyelitis. J Neuroinflammation. 2009;6:14.PubMedView ArticleGoogle Scholar
- Shen Q, Zhang R, Bhat NR. MAP kinase regulation of IP10/CXCL10 chemokine gene expression in microglial cells. Brain Res. 2006;1086(1):9–16.PubMedView ArticleGoogle Scholar
- Honda K, Yanai H, Negishi H, Asagiri M, Sato M, Mizutani T, et al. IRF-7 is the master regulator of type-I interferon-dependent immune responses. Nature. 2005;434(7034):772–7.PubMedView ArticleGoogle Scholar
- Forero A, Giacobbi NS, McCormick KD, Gjoerup OV, Bakkenist CJ, Pipas JM, et al. Simian virus 40 large T antigen induces IFN-stimulated genes through ATR kinase. J Immunol. 2014;192(12):5933–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Malik AN, Vierbuchen T, Hemberg M, Rubin AA, Ling E, Couch CH, et al. Genome-wide identification and characterization of functional neuronal activity-dependent enhancers. Nat Neurosci. 2014;17(10):1330–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang MQ. Computational analyses of eukaryotic promoters. BMC Bioinformatics. 2007;8 Suppl 6:S3.PubMed CentralPubMedView ArticleGoogle Scholar
- Delputte PL, Nauwynck HJ. Porcine arterivirus infection of alveolar macrophages is mediated by sialic acid on the virus. J Virol. 2004;78(15):8094–101.PubMed CentralPubMedView ArticleGoogle Scholar
- De Santa F, Totaro MG, Prosperini E, Notarbartolo S, Testa G, Natoli G. The histone H3 lysine-27 demethylase Jmjd3 links inflammation to inhibition of polycomb-mediated gene silencing. Cell. 2007;130(6):1083–94.PubMedView ArticleGoogle Scholar
- Shakespear MR, Hohenhaus DM, Kelly GM, Kamal NA, Gupta P, Labzin LI, et al. Histone deacetylase 7 promotes Toll-like receptor 4-dependent proinflammatory gene expression in macrophages. J Biol Chem. 2013;288(35):25362–74.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.