- Research article
- Open Access
Genome-wide landscape of liver X receptor chromatin binding and gene regulation in human macrophages
BMC Genomics volume 13, Article number: 50 (2012)
The liver X receptors (LXRs) are oxysterol sensing nuclear receptors with multiple effects on metabolism and immune cells. However, the complete genome-wide cistrome of LXR in cells of human origin has not yet been provided.
We performed ChIP-seq in phorbol myristate acetate-differentiated THP-1 cells (macrophage-type) after stimulation with the potent synthetic LXR ligand T0901317 (T09). Microarray gene expression analysis was performed in the same cellular model. We identified 1357 genome-wide LXR locations (FDR < 1%), of which 526 were observed after T09 treatment. De novo analysis of LXR binding sequences identified a DR4-type element as the major motif. On mRNA level T09 up-regulated 1258 genes and repressed 455 genes. Our results show that LXR actions are focused on 112 genomic regions that contain up to 11 T09 target genes per region under the control of highly stringent LXR binding sites with individual constellations for each region. We could confirm that LXR controls lipid metabolism and transport and observed a strong association with apoptosis-related functions.
This first report on genome-wide binding of LXR in a human cell line provides new insights into the transcriptional network of LXR and its target genes with their link to physiological processes, such as apoptosis.
The gene expression microarray and sequence data have been submitted collectively to the NCBI Gene Expression Omnibus http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo under accession number GSE28319.
The nuclear receptors liver X receptor (LXR) α and β (encoded by the genes NR1H3 and NR1H2) are transcription factors that act as sensors for oxidized cholesterol (oxysterols) [1–3]. High expression levels of LXRα in metabolic active tissues fit with the central role of the receptor in lipid metabolism, while LXRβ is more ubiquitously expressed . Interestingly, both LXRs are found in various cells of the immune system such as macrophages, dendritic cells and lymphocytes, which suggests an important function in the innate and adaptive immune response [5–7]. In macrophages the accumulation of excess lipoprotein-derived cholesterol activates LXR and triggers the induction of a transcriptional program for cholesterol efflux, such as ATP-binding cassette transporter (ABC) A1 (ABCA1) and ABCG1, while in parallel the receptor transrepresses inflammatory genes, such as inducible nitric oxide synthase, interleukin 1β and monocyte chemotactic protein-1 [8–11].
Oxysterols and intermediates of the biosynthetic cholesterol pathway have been identified as the natural ligands for LXR, while T0901317 (T09) is a potent synthetic LXR agonist with an EC50 of about 50 nM [12, 13]. LXRs bind to DNA as a heterodimer with the nuclear receptor retinoid X receptor (RXR) on direct repeats (DRs) of (A/G)GGTCA core binding motifs with four intervening nucleotides (DR4) . These DR4-type response elements (REs) have been identified in the regulatory regions of a number of primary LXR target genes [15–17]. Recently, the first genome-wide views of LXR binding were obtained in a murine macrophage cell line  and in mouse liver . In the murine macrophage study overexpressed biotin-tagged LXRβ was used for the chromatin immunoprecipitation (ChIP), followed by massive parallel sequencing (ChIP-Seq). De novo motif analysis identified DR4-type REs as the most highly enriched binding sequence, but only 6.3% of the LXR peaks contained such a DR4-type RE. Also motifs for the transcription factors PU.1 and AP-1 were co-enriched .
The ability of LXRs to integrate lipid metabolism and immune functions places them in an ideal position to tailor macrophage responses. Therefore, we performed a genome-wide analysis of T09-induced LXR binding in macrophage-type phorbol myristate acetate (PMA)-differentiated THP-1 human monocytic leukemia cells by ChIP-Seq using a highly specific anti-LXR antibody. In total we identified 1357 chromosomal LXR binding locations, of which 526 were observed after the T09 treatment. At the mRNA level the ligand induced 1258 genes and repressed 455 genes. We studied the LXR cistrome and identified LXR-enriched genomic regions as well as individual LXR target genes. Binding of LXR is focused on 112 genomic regions with individual constellations of binding sites and target genes. We could confirm that LXR controls functions related to lipid metabolism and observed a strong link to apoptosis. Therefore, this first report on genome-wide binding of LXR in a human cell line provides new insights into LXR target genes and their link to physiological processes, such as apoptosis.
Genome-wide binding of LXR in a human macrophage-type cell line
THP-1 human monocytic leukemia cells were treated for 3 days with PMA to induce a human macrophage-type model and subsequently treated for 60 min with the potent synthetic LXR agonist T09 or vehicle DMSO. ChIP assays were performed using an antibody specific for LXRα and β (and non-specific IgG as negative control). The LXR antibody was successfully used in regular ChIP assays [20, 21] and also applied in a very recent LXR ChIP-seq study in mouse liver . Its specificity is further demonstrated in Western blot analysis (see additional file 1: Figure S1), where the antibody recognizes both LXRα and LXRβ protein in liver from wild type mice. LXRα is not recognized in LXRα-/- mice, LXRβ is not recognized in LXRβ-/- mice and both LXR bands disappear in the LXRαβ-/- mice. Furthermore, the antibody recognizes both human LXRα and LXRβ when they are overexpressed in HeLa cells. This indicates that the antibody is specific to both mouse and human LXRα and LXRβ.
Purified chromatin samples were sequenced using a Solexa GAII platform. In order to detect genomic LXR binding locations, we used Bowtie software  for the read sequence alignment and the MACS program  for detection of statistically significant pileups of fragments when comparing to IgG. The number and overlap of detected LXR binding locations in and between the T09- and vehicle-treated cells are shown in Figure 1A (see additional file 2: Table S1 for all binding locations). Since the use of a single criterion for the selection of a representative peak set could be restrictive for the further analysis, we considered three criteria, each with different stringency.
The first high stringency criterion with false discovery rate (FDR) < 1%, fold enrichment (FE) > 4 and raw P-value < 10-10 showed a higher number of peaks in the ligand-treated sample (151 peaks) than in the vehicle-treated sample (137 peaks), suggesting that the ligand slightly induces LXR binding among the high confidence binding locations. Interestingly, the peaks of both T09- and vehicle-treated samples were enriched to exon and intron regions, but ligand treatment clearly increased the density of LXR binding close to the transcription start site (TSS) (see additional file 3: Figure S2). Due to these very stringent criteria the total number of peaks is not very high (202 peak locations), but they have proportionally higher overlap between T09- and vehicle-treated samples (42.6%) than the two following criteria showing overlaps of only 28.1% and 19.4%.
The second stringency level using a single criterion of a FDR < 1% increased the number of LXR binding sites in the T09-treated sample to 526 peaks and even more (1212 peaks) in the vehicle-treated sample. In total, this represents 1357 LXR binding locations. Loosening the stringency to FDR < 5% further increases the total number of LXR binding locations in the vehicle-treated sample (6903 peaks) compared to the ligand treatment (2817 peaks). In general, ligand treatment seems to increase the number of high confidence peaks and to decrease the number of lower confidence peaks suggesting that liganded LXR is focusing on a fewer number of genomic binding sites of possibly increased importance.
We then applied the de novo motif detection tool MEME  on the sequences within ± 100 bp of the summits of the 1357 LXR peaks with a FDR < 1%. From the top 10 results obtained with the TOMTOM tool  in the T09-treated peaks, only one sequence motif resembled a known transcription factor binding site with similarity score E < 105. This motif has a MEME score of E = 1.4 × 10-12 and the TOMTOM tool indicated its high similarity to two DR4-type LXR consensus matrices reported in the Transfac database (M00766 with E = 4.95 × 10-10 and M00647 with E = 1.60 × 10-8) (Figure 1B). As a major difference, the first two nucleotides of the spacer (positions 7 and 8) seem to form a consensus "CT" in our matrix, whereas in the best matching Transfac matrix these positions are not clearly enriched. The motifs reported by MEME of the vehicle-treated sample included several compositionally biased hits and none of them could be recognized by TOMTOM.
Further analysis of all 1357 peak sequences with FDR < 1% using the matrix screening function of the regulatory sequence analysis tools (RSAT) web server  showed that a site similar to the de novo derived DR4-type RE or its version with spacer positions 7 and 8 set equal to any nucleotide (see Methods for details) can be found in 7.4% (39 of 526 peaks) of the T09-treated peak set and in 7.2% (87 of 1212) of the vehicle-treated peak set with similarity P < 10-4. When the similarity threshold was reduced to P < 10-3 the number of DR4-type REs could be increased to 42.6% and 41.3% in T09- and vehicle-treated samples, respectively (data not shown). Due to relatively low percentage of DR4-type REs, we screened the high stringency LXR peak sequences for the nuclear receptor half-site motif RGKTCA in a DR 0-6, everted repeat (ER) 0-12 and inverted repeat (IR) 0-6 arrangement. This resulted in enrichment for DR1-, DR4- and IR1-type REs in the T09-treated sample and DR4- and ER9-type REs in the vehicle-treated sample (see additional file 4: Figure S3). Direct LXR-RXR heterodimer binding to DR1- and ER9-type REs has not yet been found, but it had been reported on an IR1-type RE .
Screening of the LXR peak sequences with all 459 non-redundant JASPAR database matrices (JASPAR core, version 10/2009), without DR4-type LXR RE among them, resulted in the following five most frequent motifs in the T09-treated sample (with similarity P < 10-4): EWSR1-FLI1 (15.4%), KLF4 (13.7%), RREB1 (13.7%), SP1 (13.1%) and the DR1-type RE PPARG:RXR (12.0%). In the vehicle-treated sample these were EWSR1-FLI1 (17.6%), MYF (11.8%), SP1 (11.8%), RREB1 (11.7%) and SPIB (10.9%). Among the top motifs SPIB represents a purine-rich ETS-type motif, which is also recognized by the transcription factor PU.1, previously found in the LXR ChIP-Seq analysis from the murine macrophage cell line RAW264.7 . The EWSR1-FLI1 matrix also represents a binding site for ETS-type transcription factors with various functions, such as cell cycle regulation and cell migration . RREB1 is known to relate to RAS-mediated cell differentiation  and SP1 and KLF4 both represent binding sites for KLF-family proteins, where KLF4 is known to drive the differentiation of immune cells . As the motif PPARG:RXR was the fifth most frequent in the T09-treated sample and the PPARG gene was also detected as up-regulated with a FC = 1.4 in our microarray at 4 h, we decided to search peaks with FDR < 1% for motifs PPARG::RXRA (DR1-type PPRE based on ChIP-Seq performed in mouse with identifier MA0065.2) and PPARG (Pal3-type PPRE based on Selex screening performed in human with identifier MA0066.1), available in JASPAR. These covered together 16.5% (87 of 526) of the T09-treated and 14.3% (173 of 1212) of the vehicle-treated peak set with a similarity P < 10-4. Interestingly, most of the peaks containing a PPARγ binding motif did not contain any DR4-type RE motif (79 of 87 peaks (90.8%) in the T09-treated and 162 of 173 peaks (93.6%) in the vehicle-treated peak set). These results indicate that some of the LXR peaks observed in ChIP-Seq data could be explained by an indirect DNA binding of LXR via other transcription factors, such as PPARγ, or by cooperative direct DNA binding of LXR together with some of the mentioned transcription factors.
A large part (42.6%) of the LXR binding locations in the high stringency set (Figure 1A) represent the well understood case of LXR being present on its genomic binding site both before and after the ligand treatment. This is illustrated at the genomic region of the well-known LXR target gene ABCA1 showing four peaks from the stringent set and one additional from the FDR < 1% set (Figure 1C). This observation corroborates our previous report of LXR binding sites on the ABCA1 gene by regular ChIP assays . Moreover, of the five sites observed at the ABCA1 gene, three contain DR4-type REs and four correspond to those detected previously in the first mouse LXR ChIP-Seq study . All five peaks were also occupied by LXR in the absence of ligand, but after T09 treatment the intensity of LXR binding increased. Interestingly, the five LXR locations and the ABCA1 TSS are contained within the same block, which is flanked by CTCF binding sites (data from in CD4+ T cells ) that are known as chromatin barrier insulators. This suggests that this genomic region displays the complete set of binding sites needed to understand the LXR regulation of the ABCA1 gene.
In summary, depending on chosen stringency criteria we detected in a human macrophage-type cell line between 202 and 8139 genomic LXR binding sites. The only identified de novo motif in the peak areas was a DR4-type RE, but matrix screening also identified binding sites for other transcription factors, such as PPARγ. The example of the ABCA1 locus indicates that our ChIP-Seq data can fully explain the LXR regulation of a T09 target gene.
Spatial clusters of LXR binding locations and regulated genes
Next we studied the genome-wide clustering of LXR binding sites and the location of target genes within the clustered chromosomal regions with enrichment of binding locations. We stimulated macrophage-type PMA-treated THP-1 cells for 4 h with T09 or vehicle (DMSO), extracted total RNA and analyzed it on Illumina Human HT-12 v3 Expression BeadChip gene expression microarrays. Statistical analysis detected 1713 regulated genes (P < 0.01), out of which 1258 were up and 455 down-regulated (see additional file 5: Table S2). The response already after 4 h suggests that these genes are primary LXR target genes.
The genome-wide organization of the high stringency set of 202 LXR binding locations (Figure 1A) and of the 1713 T09-regulated genes in their physical context is shown in Figure 2 using window density based visualization. This visualization shows the densities of ChIP-Seq peaks and regulated genes across the genome within 1 Mb windows. The window densities are weighted by the FEs of ChIP-Seq peaks and the log2 fold changes (FC) of the differentially expressed (DE) genes, which emphasizes high peaks or extremely differentially expressed genes. Subsequent analysis of the window density data using segmentation (for details, see Methods) to detect the exact borders of the peak-enriched regions resulted in the indicated 112 distinct genomic areas. Regions with ≥ 2 peaks and ≥ 3 T09-regulated genes are highlighted in red and listed in Table 1 (the complete list of all regions with extended information is summarized in additional file 6: Table S3).
Generally, the distribution of LXR binding locations in the genome-wide view (Figure 2) correlates with the density of all genes (r = 0.55) and slightly with the density of differentially expressed genes (r = 0.32), but the number of LXR peaks in the 112 hotspots does not correlate with the proportion of DE genes in these regions nor with the density of DE genes in these regions (see additional file 7: Figure S4). In a more detailed view, we correlated each of the 112 hotspot regions with the number of measured up-regulated (Figure 3A, left panel) and down-regulated (Figure 3A, right panel) genes with the number of expected up- or down-regulated genes. Using a binomial test with threshold P < 0.05 to indicate statistical significance, we identified in both groups several regions that contain an unexpected high number of DE genes. Region R26 in chromosome 5 (Figure 3B) and region R56 on chromosome 11 (Figure 3C) are examples of this type of analysis. Region R26 contains six up- and five down-regulated genes and two high stringency LXR locations present in both T09- and vehicle-treated samples. The high stringency LXR location close to the T09 target genes EGR1 (early growth response protein) and ETF1 (eukaryotic translation termination factor 1) is shown in higher resolution (Figure 3B). Region R56 contains eight down- and two up-regulated genes and five high stringency LXR locations in the T09 sample and four in the DMSO sample, two of which are overlapping. Here the region upstream of the T09 target genes GPR137 (G protein-couple receptor 137), TRMT112 (tRNA methyltransferase 11-2 homolog) and PRDX5 (peroxiredoxin 5) is shown in higher resolution (Figure 3C).
However, there are also several differences in the location of LXR binding sites and T09 target genes. One of the most obvious differences is the differential expression of several genes on the X chromosome (Figure 2). With the exception of region R112 located at the 3' end of the X chromosome no LXR binding location in neither the high stringent nor in the FDR < 1% peak set could be detected. These genes without any LXR binding locations nearby may either represent effects of T09 that are not mediated by LXR, secondary targets of LXR or effects that are not visible in the ChIP-Seq dataset at the investigated time point. It should be noted that T09 has been shown to be also a ligand of other nuclear receptor superfamily members, such as farnesoid X receptor , pregnane X receptor , retinoid acid receptor-related orphan receptor  and androgen receptor . However, based on our microarray analysis the genes encoding for farnesoid X receptor and pregnane X receptor, NR1H4 and NR1I2, are not expressed in PMA-differentiated THP-1 cells used in this study.
Taken together, the genome-wide correlation of the 202 high stringency LXR locations with 1258 up-regulated and 455 down-regulated T09 target genes indicated 112 genomic hotspots for the actions of LXR. Each of these hotspot regions provides a genomic scenario, where up to seven LXR binding sites may explain the regulation of up to eleven T09 target genes.
Identification of direct LXR target genes
In order to discriminate those T09 target genes that are directly regulated by LXR compared to indirect LXR targets or no LXR targets at all, we examined the co-location of LXR ChIP-Seq peaks with the TSS of DE genes within proximal (0 to ± 100 kb) and distal (± 0.1 to ± 1 Mb) regions of peak summits (Figure 4A). For the proximal regions we used the looser criterion with FDR < 5%, whereas for the distal region we used the more stringent criterion FDR < 1%. From the 1713 DE T09 target genes 814 showed a LXR site within their proximal region and 799 a LXR peak within their distal region, 550 of which are overlapping (Figure 4B). This selection filtered out 38% (of which 85% were up-regulated) of the T09-regulated genes resulting in 1063, 706 up-regulated and 357 down-regulated, more probable direct LXR target genes (Figure 4C).
The 1063 genes with a LXR peak within 1 Mb of their TSS are marked in the list of all T09 target genes (see additional file 5: Table S2) and those that show an adjusted P < 0.001 for differential expression and a proximal T09-induced LXR binding site are summarized in Table 2. From this selection we show the locus of the most up-regulated genes, NACA (nascent polypeptide-associated complex alpha subunit), flanked by the PTGES3 (prostaglandin E synthase 3) and PRIM1 (primase, DNA, polypeptide 1), which together form the core of region R61 (Figure 4D). Another exemplified locus is the most down-regulated gene, SLC3A2 (solute carrier family 3A2), located in region R55 (Figure 4E). In both examples there are dominant T09-induced tandem LXR peaks controlling the respective region, i.e. both the up-regulated (PTGES3, NACA and PRIM1) and the down-regulated (SLC3A2) genes appear to follow the same mechanism of a T09-induced LXR binding. In addition, we made a comparison with the LXR locations found previously in mouse macrophages  in a distance of 100 kb from the target gene TSS regions (see additional file 5: Table S2). Among all 1713 differentially expressed human genes only 7.9% (136 of 1713) showed LXR binding in the corresponding mouse region. Furthermore, among the 1063 selected target genes only 9.9% (105 of 1063) contained an LXR peak in the corresponding genomic region in mouse macrophages. Finally, from the human genes with proximal LXR peaks (Table 2) 25.8% (8 of 31) of their mouse homolog also showed LXR binding.
Next we compared altered expression of LXR target genes to the occurrence of DR4-type REs within the proximal LXR peaks (Figure 4A). Although DR4-type REs with very high similarity scores are not very common (Figure 1B), they seem to be enriched to the peaks in the vicinity of DE genes, especially when located nearby up-regulated genes (Figure 5A). This is consistent with the principal activation mechanism of LXR. The six up-regulated genes ABCA1, ABCG1, SMPDL3A (sphingomyelin phosphodiesterase, acid-like 3A), NR1H3, SCD (stearoyl-CoA desaturase) and TATDN2 (TatD DNase domain containing 2) and the three down-regulated genes CNNM4 (cyclin M4), HARS (histidyl-tRNA synthetase) and PUF60 (poly-U binding splicing factor 60 KDa) have LXR binding location with a motif highly similar (P < 10-6) to a DR4-type RE. The ABCA1 gene is one of the best-known LXR target genes and contains even three DR4-type REs within its regulatory region (Figure 1C). Another example, the SMPDL3A gene, carries two DR4-type REs within an LXR twin peak very close to its TSS (Figure 5B). Using quantitative real-time PCR (qPCR) and RNA from independently performed stimulation experiments of PMA-differentiated THP-1 cells with the synthetic LXR ligands T09 and GW3965 (GW), we validated nine representative LXR targets genes: the known targets ABCA1, ABCG1, MYLIP (myosin regulatory light chain interact), NR1H3 and SCD and the novel targets PPARG, SMPDL3A, ADM (adrenomedullin) and ACSL3 (acyl-CoA synthetase long-chain family member 3) (see additional file 8: Figure S5). After 4 h of ligand treatment all nine genes were significantly up-regulated by both LXR ligands.
In summary, 1063 of the 1713 T09 responding genes are probable direct targets of LXR, since 77% of them have at least one LXR binding site within 100 kb of their TSS and the further 23% show at least one LXR peak within 1 Mb distance. Interestingly, highly regulated LXR target genes have a higher probability to contain a DR4-type REs within the LXR peak associated with these genes.
Association of functions to LXR target genes
In order to gain further insight into the possible common functions of the 1063 presumed direct LXR target genes identified, we performed a functional annotation analysis using the DAVID tool  and Gene Ontology (GO) biological process terms. This resulted in 78 GO terms with FDR < 5% mostly related to general themes, such as metabolism, translation and RNA processing (see additional file 9: Table S4). Among these, 71 genes are related to regulation of programmed cell death (rank 42) and 73 genes to regulation of apoptosis (rank 53). Interestingly, the immune-related terms are not enriched significantly, as the best ranking immune-related term "somatic diversification of immune receptors via germline recombination" is only related to five LXR target genes with FDR = 57.3%. We also analyzed separately a more specific gene list, a subset of 1063 genes, with P < 0.001 for differential expression (Table 3). Also here we mostly found more general functions, such as translation and RNA processing, and only two terms related to lipid handling. None of the terms were directly related to immune function but the two terms "cellular response to stress" and "response to virus" were detected. Separate analysis of genes with different binding profiles of LXR in their vicinity (being present either only in T09 samples, DMSO samples or both) did not show any significant enrichment of additional pathways (data not shown). This suggests a homogeneous function of LXR targets independent of possible different regulatory mechanisms.
In order to detect the non-redundant sets of genes and associated biological themes, we performed a clustering analysis using enriched gene categories from the databases GO, Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome and CGAP tissue EST expression for the 1063 true LXR target genes. In order to ease the later inspection of results, we limited the number of genes used in the clustering analysis to the 150 most up-regulated and the 76 most down-regulated, preserving the same ratio as between all 706 and 357 up- and down-regulated LXR target genes. For the annotations, we used all that were enriched with FDR < 50%, corresponding to the raw P ≤ 0.036 in the GO enrichment analysis. We were able to use this loose criterion and obtain adequate data, because the clustering is robust to false positive annotations. The genes and associated annotations were clustered and visualized in parallel using the heatmap.2 function available in the R-package gplots . Agglomerative hierarchical clustering with an asymmetric binary distance measure was used by treating each association between a gene and an annotation as 1 and a lack of association as 0. This resulted in eight visually homogeneous clusters of associated genes and annotations (Figure 6): translation-related genes (cluster 1), oxidation- and diabetes-related genes expressed in liver (cluster 2), mRNA processing-related genes (cluster 3), nitrogen metabolism-related genes (cluster 4), programmed cell death regulation-related genes (cluster 5), ubiquitin system genes with relation to cell cycle and apoptosis (cluster 6), genes related to intracellular transport including cholesterol transport (cluster 7) and general ubiquitin system-related genes (cluster 8).
Taken together, functional annotation analysis of the 1063 true LXR target genes using GO, KEGG, Reactome and CGAP tissue EST expression databases resulted in eight clusters, in which the functions apoptosis and lipid transport are found, but no link to immune functions were observed.
This study provides the first genome-wide view of LXR binding patterns in a human cell line using ChIP-Seq assays. We performed this study in PMA-differentiated THP-1 cells, a macrophage-type cellular system, which increases our understanding not only of the well-known role of LXR in lipid metabolism and transport, but also the receptor's assumed role in innate immunity and other physiological processes.
The highly stringent analysis of the ChIP-Seq peaks obtained from PMA-treated THP-1 cells after a 60 min stimulation with the synthetic LXR ligand T09 or its vehicle DMSO yielded a surprisingly low number of 202 genome-wide binding locations of the receptor. However, from the 151 LXR sites observed in the presence of T09, 57% (86) are also bound in the absence of ligand. This supports the canonical model for nuclear receptor binding being valid for most members of the superfamily . According to this model the receptor binds genomic DNA already in the absence of ligand, probably in a complex with co-repressor and histone deacetylase proteins, and locally represses the chromatin structure. The addition of ligand induces a conformational change in the ligand-binding domain of the receptor, which then leads to dissociation of the co-repressor proteins and the recruitment of co-activators that open chromatin structure. Alternatively, co-activators act as mediators building a bridge to the basal transcription machinery, which leads to the activation of RNA polymerase II and gene transcription. According to the data presented here, this model seems to apply for a number of known LXR binding sites close to target genes, such as ABCA1.
A lower stringency in the detection of LXR ChIP-Seq peaks increased the number of sites substantially, to a total of 1357 (FDR < 1%) or even 8139 (FDR < 5%). However, with this increase of putative LXR binding sites the percentage of overlapping sites reduced to 28.1 and 19.4%, respectively, while the percentage of apparent binding sites in the absence of ligand increased to 61.2 and 65.4%, respectively. We assume that most of the latter sites are not genomic locations from which LXR initiates gene activation, but rather unspecific contact points of the genome with no functional impact. However, the tendency that ligand stimulation reduces the number of sites to less than half (from 1212 to 526 and from 6903 to 2817) suggests that in its ligand-activated state LXR makes a more focused selection of genomic targets.
Surprisingly, only 7.3% of the LXR peak summit sequences contain a DR4-type RE, which, based on in vitro studies, is the only high affinity DNA binding site for LXR-RXR heterodimers. This fits with the observations of the two presently published mouse LXR ChIP-Seq studies, where only 6.3%  or up to 8%  of the LXR peaks contained a DR4-type RE. In a very recent ChIP-on-chip report, 2035 LXRβ-RXR binding sites were identified within promoter regions of normal human epidermal keratinocytes . Although 1666 (82%) of these rather large genomic fragments contain a kind of DR4-type RE, only 142 (7.0%) of them are highly scored, resulting in a very similar percentage as in our human and the mouse ChIP-Seq studies. For comparison, in our ChIP-Seq study on genomic vitamin D receptor (VDR) locations in undifferentiated THP-1 cells , we found only 31% of the VDR peak summits containing a VDR-specific DR3-type RE. We obtained the same result also by a re-analysis of VDR ChIP-Seq data from human lymphoblast cell lines [42, 43]. These observations indicate that either the binding specificity of LXR (and VDR) is in the chromatin context far different than in in vitro assays (maybe due to cooperative binding with other transcription factors, such as PPARγ) or LXR is not directly contacting DNA, but sits "piggyback" on another DNA-binding transcription factor. Interestingly, in our T09 microarray the PPARG gene is 1.4-fold up-regulated, a value that was also confirmed by qPCR with two different LXR ligands. In parallel, we found for 16.5% of all LXR peaks in the T09-treated sample a PPARγ binding site. This suggests that PPARγ may be involved in the regulation of a subset of the LXR target genes. Boergesen et al. reported very recently that in mouse liver many genomic binding sites of LXR are also bound by PPARα, which is the predominant PPAR subtype in liver . They suggest that PPARα helps LXR in recognizing its preferred genomic locations, a concept that may apply also for PPARγ in macrophages. Furthermore, Boergesen et al. also found a higher percentage of DR1-type REs below their LXR ChIP-Seq peaks than DR4-type REs, but in sum this can explain only less than 25% of all LXR locations. Therefore, they assume that both LXR and PPAR are far more promiscuous in the recognition of their binding sites as suggested by in vitro studies.
Nevertheless, we found DR4-type REs as functional LXR binding sites in a number of prominent and highly ligand-responsive LXR target genes, such as ABCA1, ABCG1, NR1H3 and SCD. Moreover, in a reasonable number of gene regulatory scenarios two or more genomic LXR binding sites seem to work together in the regulation of a gene cluster (for example, NACA, GPR137 and their respective neighboring genes) or individual genes (for example, ABCA1, SLC3A2 and SMPDL3A). These multiple RE arrangements confirm observations from single gene analyses with other nuclear receptors, such as the VDR [44–46].
More important than the actual number of individual LXR binding sites may be their spatial clustering along the genome. We identified 112 genomic regions of an average size of 1.7 Mb that contain clusters of highly stringent LXR binding sites and up to 11 T09 target genes. In total we identified 432 T09 target genes within these genomic LXR hotspots. Although each of these regions displays a rather individual arrangement of LXR binding sites in relation to up- and down-regulated genes, they seem to represent the core of the genome-wide activities of LXR. However, 13 of the 112 regions do not have any T09 target gene, suggesting that these stringent LXR locations may have effects on more distant target genes.
In total our microarray analysis of the 4 h T09-stimulated human macrophage-type cells lists 1713 genes, 73% of which are up-regulated, i.e. more than double as many genes are up-regulated than down-regulated. Out of these T09 target genes, 432 (25%) are found within the 112 genomic LXR hotspots, 814 genes (48%) have an LXR binding site within 100 kb distance from the respective gene's TSS and for a further 249 genes (15%) an LXR location within 1 Mb of the TSS could be identified. This rate is similar to other reports comparing ChIP-Seq and microarray results . Nevertheless, this raises the question about the regulation of the 650 T09 target genes without an LXR peak. These genes may be secondary LXR targets. As only 15% of them are down-regulated, a trans-repression mechanism is not very likely.
Our study confirmed a number of known primary LXR target genes, such as ABCA1, ABCG1, MYLIP, NR1H3 and SCD, but we identified also a number of previously unknown, novel LXR targets. The most up-regulated LXR target gene, NACA, encodes for the nascent polypeptide associated complex alpha protein, which is associated to translation and protein folding related processes and when depleted is responsible of triggering endoplasmic reticulum stress-driven apoptosis in hypoxic cells . The PTGES3 gene, which is co-located with the NACA gene, encodes for the co-chaperone protein prostaglandin E synthase 3 and is required together with the primary chaperone proteins for proper folding and functioning of the glucocorticoid receptor and other steroid receptors . The most down-regulated gene, SLC3A2, is related to various processes, for instance in the migration of leukocytes from blood to the central nervous system , but also to serum cholesterol levels . Also the SMPDL3A gene is an interesting new LXR target, but it has not yet been well characterized as there is only one report linking SMPDL3A to bladder tumorigenesis .
LXR actions have been related to a number of autoimmune diseases, such as multiple sclerosis [52, 53], rheumatoid arthritis [54, 55] and experimental autoimmune encephalomyelitis . Interestingly, the association of LXR with autoimmune and metabolic diseases is also one of the major results of the annotation analysis of the 1063 T09 responding genes that we consider as true LXR target genes. The analysis showed the expected relation to lipid metabolism and transport genes, but did not provide any significant link to genes related to innate immunity. However, the LXR regulation of inflammatory cytokines is generally observed in experimental settings, where lipopolysaccharide is used for stimulation and where longer LXR ligand treatment times are allowed . Under these conditions inflammatory cytokines are induced via the transcription factors AP1 and NF-κB and can be inhibited via tethered LXR. In this study, we used PMA for the differentiation of THP-1 cells and T09 treatment of only 4 h. Interestingly, under our experiment conditions, we observed a strong association with apoptosis as more than 70 genes within the 1063 candidates are related to this physiological process. The link of LXRs to apoptosis has is already been reported, not only in macrophages , but also in pancreatic β-cells  and different cancer cells [60, 61].
We present here the first genome-wide view on LXR locations in a human macrophage-type cellular system. The core action of LXR is focused on 112 genomic hotspots that contain 432 target genes. In total 1063 of the 1713 LXR target genes can be explained by a direct action of LXR, many of which have not been reported before. These genes are related to lipid metabolism and transport and to apoptosis, but not directly to immune functions.
THP-1 human monocytic leukemia cells were grown in RPMI 1640 supplemented with 10% fetal calf serum, 2 mM L-glutamine, 0.1 mg/ml streptomycin and 100 U/ml penicillin and the cells were kept at 37°C in a humidified 95% air/5% CO2 incubator. For differentiation into macrophage-type cells the THP-1 cells were incubated for 3 days with 20 nM PMA. Prior to stimulation with 1 μM of the synthetic LXR ligands T09 or GW (both from Sigma-Aldrich) or vehicle (dimethyl sulfoxide (DMSO), final concentration 0.1%) the medium was replaced by RPMI 1640 supplemented with 5% lipid-depleted fetal calf serum, 2 mM L-glutamine, 0.1 mg/ml streptomycin and 100 U/ml penicillin and 20 nM PMA.
Generation of LXR antibodies
Rat LXRα and LXRβ proteins were purified as previously described . Rabbits were immunized using a standard immunization program at Agrisera (Vännäs, Sweden). In brief, four injections with a total of 0.25 mg of both LXRα and LXRβ protein were performed and serum was collected after 15 weeks. Purified LXRα and LXRβ proteins (1.4 mg) in 0.2 M NaHCO3, 0.5 M NaCl, pH 8.3 buffer were coupled on a N-hydroxysuccinimide activated matrix column and the immunized rabbit serum was added to the column, which was washed according to standard procedure. Elution of anti-LXR antibody was performed in 10 cycles and all fractions were pooled.
PMA-differentiated THP-1 cells (per condition nine 175 cm2 flasks, each 2 × 107 cells, density of 1 × 106 cells/ml) were treated for 60 min with 1 μM T09 or vehicle (DMSO). Then nuclear proteins were cross-linked to DNA by adding formaldehyde directly to the medium to a final concentration of 1% following incubation for 10 min at room temperature on a rocking platform. Cross-linking was stopped by adding glycine to a final concentration of 0.15 M and incubating at room temperature for 10 min on a rocking platform. The cells were collected by centrifugation and washed twice with ice cold PBS (140 mM NaCl, 2.7 mM KCl, 1.5 mM KH2PO4, 8.1 mM Na2HPO4.2H2O). Nuclear extraction was performed by adding 500 μl PIPES buffer (5 mM Pipes pH 8.0, 85 mM KCl, 0.5% Nonidet P-40, protease inhibitors), incubating for 10 min on ice and removing cytoplasmic components by centrifugation. Nuclear pellets were dissolved in 500 μl SDS lysis buffer (1% SDS, 10 mM EDTA, protease inhibitors, 50 mM Tris-HCl, pH 8.1) and incubated 10 min on ice. Lysates were sonicated at high power with 22 × 30 s pulses in a Bioruptor (Diagenode, Liège, Belgium) to result in DNA fragments of 100 to 600 bp. Cellular debris were removed by centrifugation. Aliquots of 100 μl of the lysate were diluted 1:10 in ChIP dilution buffer (0.01% SDS, 1.1% Triton X-100, 1.2 mM EDTA, 167 mM NaCl, protease inhibitors, 16.7 mM Tris-HCl, pH 8.1) and 2 μg of anti-LXR antibody  or non-specific anti-IgG rabbit (sc-2027, Santa Cruz Biotechnologies, Santa Cruz, CA, USA) were added and the samples were incubated for overnight at 4°C on a rotating platform. The immunocomplexes were collected using 60 μl of BSA-coated protein A agarose bead slurry (Millipore) for 3 h at 4°C with rotation. The beads were washed sequentially for 4 min in rotating platform with 1 ml of the following buffers: low salt wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 150 mM NaCl, 20 mM Tris-HCl, pH 8.1), high salt wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 500 mM NaCl, 20 mM Tris-HCl, pH 8.1) and LiCl wash buffer (0.25 M LiCl, 1% Nonidet P-40, 1% sodium deoxycholate, 1 mM EDTA, 10 mM Tris-HCl, pH 8.1). Finally, the beads were washed twice with 1 ml TE buffer (1 mM EDTA, 10 mM Tris-HCl, pH 8.0) and the immune complexes were eluted twice using 200 μl elution buffer (1% SDS, 100 mM NaHCO3) for 15 min at room temperature with rotation. The supernatants were combined and the immune complexes were reverse cross-linked overnight at 65°C in the presence of proteinase K (Fermentas) in a final concentration of 0.1 mg/ml. DNA was extracted with the ChIP DNA Clean & Concentrator Kit (Zymo Research Cooperation, HiSS Diagnostics, Freiburg, Germany) according to manufacturer's instructions and eluted in 40 μl nuclease-free H2O. The ChIP templates were sequenced using a Solexa Gene Analyzer II platform (Illumina) at 36 bp read length using standard manufacturer protocols at the Genomics Core Facility in Heidelberg, Germany.
ChIP-seq data analysis
Some of our following in-house bioinformatics tools were already described recently [41, 63]. Alignment of sequence reads produced by T09-treated anti-LXR immunoprecipitated sample, vehicle treated anti-LXR immunoprecipitated sample and the IgG immunoprecipitated negative control sample against the reference genome of version hg19 was done using Bowtie software version 0.12.2 . Command line arguments used with Bowtie were: bowtie -n 1 -m 1 -e 70 -l 28 -k 1 -t -p 8 -q -S --best hg19 input_file_name output_file_name. MACS program version 184.108.40.206  was used for finding statistically significant peaks from the alignment sequences using the following arguments: macs --pvalue = 1e-3 --nomodel --wig -t input_sample_file_name -c input_control_file_name --tsize = 36 --format = BAM --name = analysis_topic --mfold = 13 --shiftsize = 250 --bw = 250 --verbose = 3. Subsequent refinement of MACS peaks was done using the PeakSplitter program  with arguments: PeakSplitter -p peak_folder_name -w aligned_read_wig_folder_name -o output_folder -c 5 -v 0.6 -f. An in-house R script was further used to calculate FEs, P-values and FDR estimates for the found subpeaks using an identical approach to MACS 220.127.116.11. Since splitting the original peaks into several subpeaks may produce large sets of weaker flanking residual peaks, distorting for example the genomic and FE distributions of the peaks, we kept for further analysis only those peaks that either had FDR < 1% or were the best subpeak (by FDR) within their parental MACS peak.
Further data analyses were conducted using an in-house R pipeline containing tools for identifying peaks overlapping in two samples, the analysis of genomic distribution of peaks and the integration of peak and gene expression datasets. In the analysis of overlapping peaks, we require that the narrower of the two overlapping peaks shares at least 30% overlap with the broader peak. Under these conditions also any weaker peak (still fulfilling the criterion of raw P < 10-3 reported by the MACS program) in one sample can have FDR < 5%, FDR < 1% or high stringency peak sets, when the overlapping peak in the other sample fulfills the criterion. Since many genomic positions cannot be uniquely assigned to the selected 10 types of genomic elements used in the analysis of peak distributions therein, in the genomic element analysis a prioritization scheme was employed, where the peaks were uniquely overlapped to the elements in a step-wise exclusive scheme, starting from coding elements (UTRs and coding exons), and then moving to introns and outward from the gene.
For comparison with the ChIP-Seq data from a mouse macrophage cell line  all peak coordinates from that study were mapped to the human hg19 genome version using Batch Coordinate Conversion tool available at the UCSC Genome Browser .
De novo analysis of LXR binding locations was performed using stand alone version of MEME  on sequences within ± 100 bp of the summits of the LXR peaks. Peak sets with FDR < 1% and FDR < 5% with different FE cutoffs FE > 1, FE > 2 etc. from T09- and vehicle-treated samples were analyzed separately in a batch run. The analysis of peak sequences from the T09-treated sample resulted in DR4-type REs in the top 10 of the MEME results with FDR < 1% peaks, when using the cutoff FE > 2 (motif shown in Figure 1B) or higher (tested until FE > 10). DR4-type REs could not be detected, when similar de novo analysis for the peak sequences obtained from the vehicle-treated sample or all peak sequences with FDR < 5% from the T09-treated sample were performed.
Identification of DR4-type REs within LXR peak regions was performed using the RSAT matrix scan tool available at http://rsat.ulb.ac.be/rsat. Two matrices were used as a model for a DR4-type RE: the de novo detected matrix (Figure 1B) and the same matrix modified from the positions 7 and 8 within the spacer to resemble more the DR4-type RE known in the literature. The modification was made by setting at these positions the frequency of any nucleotide equal (3). The analysis of JASPAR matrices (JASPAR core, version 10/2009) was performed in similar manner. For the background model, the input peak sequences were used to take into account the nucleotide content within these regions. Also for the analysis of DR-, ER- and IR-type REs the RSAT dna-pattern tool was used. One mismatch was allowed for an individual nuclear receptor binding site.
Total RNA was extracted using High Pure RNA Isolation kit (Roche). cDNA synthesis was performed using Transcriptor First Strand cDNA Synthesis Kit (Roche), using 1 μg of total RNA as a template and 50 pmol oligo(dT)18 primers. qPCR was performed using a LightCycler® 480 System (Roche). The reactions were performed using 4 pmol of reverse and forward primers, 4 μl cDNA template and FastStart SYBR Green Master (Roche) in a total volume of 10 μl. In the PCR reaction the DNA templates were pre-denaturated at for 10 min at 95°C, followed by amplification steps cycles of 20 s denaturation at 95°C, 15 s annealing at primer specific temperatures (see additional file 10: Table S5), 15 s elongation at 72°C and a final elongation for 10 min at 72°C. PCR product quality was monitored using post-PCR melt curve analysis. Fold inductions were calculated using the formula 2-(ΔΔCt), where ΔΔCt is ΔCt(stimulus)-ΔCt(solvent), ΔCt is Ct(target gene)-Ct(contol gene) and the Ct is the cycle, at which the threshold is crossed. Relative expression levels of the target genes were normalized to the internal control gene RPLP0.
Total RNA was checked for RNA integrity using a Biorad Experion automated electrophoresis system (Biorad, Nazareth, Belgium). None of the RNA samples showed any sign of degradation. Triplicate samples were analyzed with HumanHT-12 v3 Expression BeadChips from Illumina (San Diego, CA, USA) at the Finnish Microarray Centre (Turku, Finland). Raw data are available at GEO under accession GSE28319. Analyses of microarray data were performed using R statistical software version 2.11  with associated libraries from Bioconductor project version 2.6 . Data were normalized using VST transformation and RSN normalization used as standard approach for Illumina arrays. Normalized data were filtered in order to remove probes without detected signal for any of the samples. Probe sets that were not linked to any known or predicted human gene were also filtered out. Linear Models for Microarray Data (limma) package  using linear model fitting for statistical testing with empirical Bayes variance smoothing procedure was applied to detection of differentially expressed genes. Obtained P-values were corrected for multiple testing using the Benjamini & Hochberg FDR procedure . For downstream analysis, GO  biological process terms were tested for enrichment.
Spatial clustering of LXR binding locations
For the analysis of spatial clusters of LXR binding locations we first performed a density analysis of genomic coordinates of LXR peak summit locations in T09- and vehicle-treated samples and the starting coordinates of up- and down-regulated genes. This was made using the standard R function "density" with a 1 Mb size for the sliding window with 0.5 Mb steps over each chromosome using the default Gaussian window kernel function. The density values resulting from each of these windows were weighted using FE values for LXR binding locations or logarithmic FCs between T09- and vehicle-treated samples for the locations of regulated genes. In order to detect the exact borders of the LXR hotspot regions, we performed a subsequent clustering of the density data representing the LXR binding locations by using methods developed originally for the analysis of array CGH data implemented in R package SegClust [71, 72]. First, the density vectors representing LXR binding site distributions in T09- and vehicle-treated samples were combined by taking the maximum value for the two samples in each genomic location. The resulting combined density data were used for the "segmean" function, implemented in SegClust, performing a dynamic programming-based algorithm for finding the optimal breakpoints in terms of changes of mean for a fixed number of breakpoints K. The method was applied using various K starting from 1 ending to the number of analyzed data points. The results were further used for the "segselect" function, which uses an adaptive model selection  between outcomes with different numbers of breakpoints K. A result optimal in terms of this model selection was chosen as representative for each chromosome. The resulting chromosomal regions were selected for containing at least one LXR binding location within the high stringency set of peaks resulting in 112 separate regions. The region borders were widened 100 kb upstream of the 5'-border and 100 kb downstream of the 3'-border of each region to cover putative target genes in these regions.
ATP-binding cassette transporter
acyl-CoA synthetase long-chain family member 3
direct repeat spaced by 4 nucleotides
early growth response protein
eukaryotic translation termination factor 1
Ewing sarcoma breakpoint region 1
false discovery rate
friend leukemia virus integration 1
G protein-coupled receptor 137
krüppel-like factor 4
liver X receptor
myosin regulatory light chain interact
nascent polypeptide-associated complex alpha subunit
phorbol myristate acetate
peroxisome proliferator-activated receptor
peroxiredoxin 5, PRIM1, primase, DNA, polypeptide 1
prostaglandin E synthase 3
poly-U binding splicing factor 60 KDa
position weight matrix
quantitative real-time PCR
ribosomal protein large P0
Ras-responsive element-binding protein 1
regulatory sequence analysis tools
retinoid X receptor
solute carrier family 3A2
sphingomyelin phosphodiesterase: acid-like 3A
Sp1 transcription factor
Spi-B transcription factor
TatD DNase domain containing 2
tRNA methyltransferase 11-2 homolog
transcription start site
vitamin D receptor.
Apfel R, Benbrook D, Lernhardt E, Ortiz MA, Salbert G, Pfahl M: A novel orphan receptor specific for a subset of thyroid hormone-responsive elements and its interaction with the retinoid/thyroid hormone receptor subfamily. Mol Cell Biol. 1994, 14: 7025-7035.
Janowski BA, Willy PJ, Devi TR, Falck JR, Mangelsdorf DJ: An oxysterol signalling pathway mediated by the nuclear receptor LXRα. Nature. 1996, 383: 728-731. 10.1038/383728a0.
Willy PJ, Umesono K, Ong ES, Evans RM, Heyman RA, Mangelsdorf DJ: LXR, a nuclear receptor that defines a distinct retinoid response pathway. Genes Dev. 1995, 9: 1033-1045. 10.1101/gad.9.9.1033.
Peet DJ, Turley SD, Ma W, Janowski BA, Lobaccaro JM, Hammer RE, Mangelsdorf DJ: Cholesterol and bile acid metabolism are impaired in mice lacking the nuclear oxysterol receptor LXRα. Cell. 1998, 93: 693-704. 10.1016/S0092-8674(00)81432-4.
Joseph SB, Bradley MN, Castrillo A, Bruhn KW, Mak PA, Pei L, Hogenesch J, O'Connell R M, Cheng G, Saez E, Miller JF, Tontonoz P: LXR-dependent gene expression is important for macrophage survival and the innate immune response. Cell. 2004, 119: 299-309. 10.1016/j.cell.2004.09.032.
Bensinger SJ, Bradley MN, Joseph SB, Zelcer N, Janssen EM, Hausner MA, Shih R, Parks JS, Edwards PA, Jamieson BD, Tontonoz P: LXR signaling couples sterol metabolism to proliferation in the acquired immune response. Cell. 2008, 134: 97-111. 10.1016/j.cell.2008.04.052.
Töröcsik D, Barath M, Benko S, Szeles L, Dezso B, Poliska S, Hegyi Z, Homolya L, Szatmari I, Lanyi A, Nagy L: Activation of liver X receptor sensitizes human dendritic cells to inflammatory stimuli. J Immunol. 2010, 184: 5456-5465. 10.4049/jimmunol.0902399.
Chawla A, Boisvert WA, Lee C, Laffitte BA, Barak Y, Joseph SB, Liao D, Nagy L, Edwards PA, Curtiss LK, Evans RM, Tontonoz P: A PPARγ-LXR-ABCA1 pathway in macrophages is involved in cholesterol efflux and atherogenesis. Mol Cell. 2001, 7: 161-171. 10.1016/S1097-2765(01)00164-2.
Joseph SB, McKilligin E, Pei L, Watson MA, Collins AR, Laffitte BA, Chen M, Noh G, Goodman J, Hagger GN, Tran J, Tippin TK, Wang X, Lusis AJ, Hsueh WA, Law RE, Collins JL, Willson TM, Tontonoz P: Synthetic LXR ligand inhibits the development of atherosclerosis in mice. Proc Natl Acad Sci USA. 2002, 99: 7604-7609. 10.1073/pnas.112059299.
Castrillo A, Joseph SB, Vaidya SA, Haberland M, Fogelman AM, Cheng G, Tontonoz P: Crosstalk between LXR and toll-like receptor signaling mediates bacterial and viral antagonism of cholesterol metabolism. Mol Cell. 2003, 12: 805-816. 10.1016/S1097-2765(03)00384-8.
Ghisletti S, Huang W, Ogawa S, Pascual G, Lin ME, Willson TM, Rosenfeld MG, Glass CK: Parallel SUMOylation-dependent pathways mediate gene- and signal-specific transrepression by LXRs and PPARγ. Mol Cell. 2007, 25: 57-70. 10.1016/j.molcel.2006.11.022.
Repa JJ, Turley SD, Lobaccaro JA, Medina J, Li L, Lustig K, Shan B, Heyman RA, Dietschy JM, Mangelsdorf DJ: Regulation of absorption and ABC1-mediated efflux of cholesterol by RXR heterodimers. Science. 2000, 289: 1524-1529. 10.1126/science.289.5484.1524.
Schultz JR, Tu H, Luk A, Repa JJ, Medina JC, Li L, Schwendner S, Wang S, Thoolen M, Mangelsdorf DJ, Lustig KD, Shan B: Role of LXRs in control of lipogenesis. Genes Dev. 2000, 14: 2831-2838. 10.1101/gad.850400.
Quack M, Frank C, Carlberg C: Differential nuclear receptor signalling from DR4-type response elements. J Cell Biochem. 2002, 86: 601-612. 10.1002/jcb.10247.
Luo Y, Tall AR: Sterol upregulation of human CETP expression in vitro and in transgenic mice by an LXR element. J Clin Invest. 2000, 105: 513-520. 10.1172/JCI8573.
Kaplan R, Zhang T, Hernandez M, Gan FX, Wright SD, Waters MG, Cai TQ: Regulation of the angiopoietin-like protein 3 gene by LXR. J Lipid Res. 2003, 44: 136-143. 10.1194/jlr.M200367-JLR200.
Uehara Y, Miura S, von Eckardstein A, Abe S, Fujii A, Matsuo Y, Rust S, Lorkowski S, Assmann G, Yamada T, Saku K: Unsaturated fatty acids suppress the expression of the ATP-binding cassette transporter G1 (ABCG1) and ABCA1 genes via an LXR/RXR responsive element. Atherosclerosis. 2007, 191: 11-21. 10.1016/j.atherosclerosis.2006.04.018.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK: Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010, 38: 576-589. 10.1016/j.molcel.2010.05.004.
Boergesen M, Pedersen TA, Gross B, van Heeringen SJ, Hagenbeek D, Bindesbøll C, Caron S, Lalloyer F, Steffensen KR, Nebb H, Gustafsson JA, Stunnenberg HG, Staels B, Mandrup S: Genome-wide profiling of LXR, RXR and PPARα in mouse liver reveals extensive sharing of binding sites. Mol Cell Biol. 2012, 32: 852-867.
Jakobsson T, Venteclef N, Toresson G, Damdimopoulos AE, Ehrlund A, Lou X, Sanyal S, Steffensen KR, Gustafsson JA, Treuter E: GPS2 is required for cholesterol efflux by triggering histone demethylation, LXR recruitment, and coregulator assembly at the ABCG1 locus. Mol Cell. 2009, 34: 510-518. 10.1016/j.molcel.2009.05.006.
Venteclef N, Jakobsson T, Ehrlund A, Damdimopoulos A, Mikkonen L, Ellis E, Nilsson LM, Parini P, Jänne OA, Gustafsson JA, Steffensen KR, Treuter E: GPS2-dependent corepressor/SUMO pathways govern anti-inflammatory actions of LRH-1 and LXRβ in the hepatic acute phase response. Genes Dev. 2010, 24: 381-395. 10.1101/gad.545110.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.
Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.
Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS: Quantifying similarity between motifs. Genome Biol. 2007, 8: R24-10.1186/gb-2007-8-2-r24.
Turatsinze JV, Thomas-Chollier M, Defrance M, van Helden J: Using RSAT to scan genome sequences for transcription factor binding sites and cis-regulatory modules. Nat Protoc. 2008, 3: 1578-1588. 10.1038/nprot.2008.97.
Landrier JF, Grober J, Demydchuk J, Besnard P: FXRE can function as an LXRE in the promoter of human ileal bile acid-binding protein (I-BABP) gene. FEBS Lett. 2003, 553: 299-303. 10.1016/S0014-5793(03)01033-0.
Guillon N, Tirode F, Boeva V, Zynovyev A, Barillot E, Delattre O: The oncogenic EWS-FLI1 protein binds in vivo GGAA microsatellite sequences with potential transcriptional activation function. PLoS One. 2009, 4: e4932-10.1371/journal.pone.0004932.
Thiagalingam A, De Bustros A, Borges M, Jasti R, Compton D, Diamond L, Mabry M, Ball DW, Baylin SB, Nelkin BD: RREB-1, a novel zinc finger protein, is involved in the differentiation response to Ras in human medullary thyroid carcinomas. Mol Cell Biol. 1996, 16: 5335-5345.
Lebson L, Gocke A, Rosenzweig J, Alder J, Civin C, Calabresi PA, Whartenby KA: Cutting edge: The transcription factor Kruppel-like factor 4 regulates the differentiation of Th17 cells independently of RORγt. J Immunol. 185: 7161-7164.
Venkateswaran A, Laffitte BA, Joseph SB, Mak PA, Wilpitz DC, Edwards PA, Tontonoz P: Control of cellular cholesterol efflux by the nuclear oxysterol receptor LXRα. Proc Natl Acad Sci USA. 2000, 97: 12097-12102. 10.1073/pnas.200367697.
Cuddapah S, Jothi R, Schones DE, Roh TY, Cui K, Zhao K: Global analysis of the insulator binding protein CTCF in chromatin barrier regions reveals demarcation of active and repressive domains. Genome Res. 2009, 19: 24-32.
Houck KA, Borchert KM, Hepler CD, Thomas JS, Bramlett KS, Michael LF, Burris TP: T0901317 is a dual LXR/FXR agonist. Mol Genet Metab. 2004, 83: 184-187. 10.1016/j.ymgme.2004.07.007.
Mitro N, Vargas L, Romeo R, Koder A, Saez E: T0901317 is a potent PXR ligand: implications for the biology ascribed to LXR. FEBS Lett. 2007, 581: 1721-1726. 10.1016/j.febslet.2007.03.047.
Kumar N, Solt LA, Conkright JJ, Wang Y, Istrate MA, Busby SA, Garcia-Ordonez RD, Burris TP, Griffin PR: The benzenesulfoamide T0901317 [N-(2,2,2-trifluoroethyl)-N-[4-[2,2,2-trifluoro-1-hydroxy-1-(trifluorometh yl)ethyl]phenyl]-benzenesulfonamide] is a novel retinoic acid receptor-related orphan receptor-alpha/gamma inverse agonist. Mol Pharmacol. 2010, 77: 228-236. 10.1124/mol.109.060905.
Chuu CP, Chen RY, Hiipakka RA, Kokontis JM, Warner KV, Xiang J, Liao S: The liver X receptor agonist T0901317 acts as androgen receptor antagonist in human prostate cancer cells. Biochem Biophys Res Commun. 2007, 357: 341-346. 10.1016/j.bbrc.2007.03.116.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.
Bolker B, Bonebakker L, Gentleman R, Huber W, Liaw A, Lumley T, Maechler M, Magnusson A, Moeller S, Schwartz M, Venables B: gplots: various R programming tools for plotting data. R package version 280. 2010
Moras D, Gronemeyer H: The nuclear receptor ligand-binding domain: structure and function. Curr Opin Cell Biol. 1998, 10: 384-391. 10.1016/S0955-0674(98)80015-X.
Shen Q, Bai Y, Chang KC, Wang Y, Burris TP, Freedman LP, Thompson CC, Nagpal S: Liver X receptor-retinoid X receptor (LXR-RXR) heterodimer cistrome reveals co-ordination of LXR and AP1 signaling in keratinocytes. J Biol Chem. 2011, 286: 14554-14563. 10.1074/jbc.M110.165704.
Heikkinen S, Väisänen S, Pehkonen P, Seuter S, Benes V, Carlberg C: Nuclear hormone 1α,25-dihydroxyvitamin D3 elicts a genome-wide shift in the locations of VDR chromatin accupancy. Nucl Acids Res. 2011, 39: 9181-9193. 10.1093/nar/gkr654.
Ramagopalan SV, Heger A, Berlanga AJ, Maugeri NJ, Lincoln MR, Burrell A, Handunnetthi L, Handel AE, Disanto G, Orton SM, Watson CT, Morahan JM, Giovannoni G, Ponting CP, Ebers GC, Knight JC: A ChIP-seq defined genome-wide map of vitamin D receptor binding: associations with disease and evolution. Genome Res. 2010, 20: 1352-1360. 10.1101/gr.107920.110.
Heikkinen S, Seuter S, Carlberg C: The first genome-wide view on vitamin D receptor locations and their mechanistic implications. Anticancer Res. 2012, 32: 271-282.
Väisänen S, Dunlop TW, Sinkkonen L, Frank C, Carlberg C: Spatio-temporal activation of chromatin on the human CYP24 gene promoter in the presence of 1α,25-dihydroxyvitamin D3. J Mol Biol. 2005, 350: 65-77. 10.1016/j.jmb.2005.04.057.
Saramäki A, Banwell CM, Campbell MJ, Carlberg C: Regulation of the human p21(waf1/cip1) gene promoter via multiple binding sites for p53 and the vitamin D3 receptor. Nucleic Acids Res. 2006, 34: 543-554. 10.1093/nar/gkj460.
Meyer MB, Goetsch PD, Pike JW: A downstream intergenic cluster of regulatory enhancers contributes to the induction of CYP24A1 Expression by 1α,25-dihydroxyvitamin D3. J Biol Chem. 2010, 285: 15599-15610. 10.1074/jbc.M110.119958.
Hotokezaka Y, van Leyen K, Lo EH, Beatrix B, Katayama I, Jin G, Nakamura T: alphaNAC depletion as an initiator of ER stress-induced apoptosis in hypoxia. Cell Death Differ. 2009, 16: 1505-1514. 10.1038/cdd.2009.90.
Pratt WB, Morishima Y, Murphy M, Harrell M: Chaperoning of glucocorticoid receptors. Handb Exp Pharmacol. 2006, 111-138.
Callahan MK, Ransohoff RM: Analysis of leukocyte extravasation across the blood-brain barrier: conceptual and technical aspects. Curr Allergy Asthma Rep. 2004, 4: 65-73. 10.1007/s11882-004-0046-9.
Igl W, Johansson A, Wilson JF, Wild SH, Polasek O, Hayward C, Vitart V, Hastie N, Rudan P, Gnewuch C, Schmitz G, Meitinger T, Pramstaller PP, Hicks AA, Oostra BA, van Duijn CM, Rudan I, Wright A, Campbell H, Gyllensten U: Modeling of environmental effects in genome-wide association studies identifies SLC2A2 and HP as novel loci influencing serum cholesterol levels. PLoS Genet. 2010, 6: e1000798-10.1371/journal.pgen.1000798.
Wright KO, Messing EM, Reeder JE: Increased expression of the acid sphingomyelinase-like protein ASML3a in bladder tumors. J Urol. 2002, 168: 2645-2649. 10.1016/S0022-5347(05)64236-X.
Liu X, Steffensen KR, Sanna A, Arru G, Fois ML, Rosati G, Sotgiu S, Link H, Gustafsson JA, Huang YM: Anti-inflammatory nuclear receptor superfamily in multiple sclerosis patients from Sardinia and Sweden. Neurobiol Dis. 2005, 20: 961-968. 10.1016/j.nbd.2005.06.001.
Giorelli M, Livrea P, Minervini MG, Trojano M: Immunomodulatory properties of increased levels of liver X receptor beta in peripheral blood mononuclear cells from multiple sclerosis patients. Exp Neurol. 2007, 204: 759-766. 10.1016/j.expneurol.2007.01.013.
Chintalacharuvu SR, Sandusky GE, Burris TP, Burmer GC, Nagpal S: Liver X receptor is a therapeutic target in collagen-induced arthritis. Arthritis Rheum. 2007, 56: 1365-1367. 10.1002/art.22528.
Asquith DL, Miller AM, Hueber AJ, McKinnon HJ, Sattar N, Graham GJ, McInnes IB: Liver X receptor agonism promotes articular inflammation in murine collagen-induced arthritis. Arthritis Rheum. 2009, 60: 2655-2665. 10.1002/art.24717.
Hindinger C, Hinton DR, Kirwin SJ, Atkinson RD, Burnett ME, Bergmann CC, Stohlman SA: Liver X receptor activation decreases the severity of experimental autoimmune encephalomyelitis. J Neurosci Res. 2006, 84: 1225-1234. 10.1002/jnr.21038.
Joseph SB, Castrillo A, Laffitte BA, Mangelsdorf DJ, Tontonoz P: Reciprocal regulation of inflammation and lipid metabolism by liver X receptors. Nat Med. 2003, 9: 213-219. 10.1038/nm820.
Seres L, Cserepes J, Elkind NB, Torocsik D, Nagy L, Sarkadi B, Homolya L: Functional ABCG1 expression induces apoptosis in macrophages and other cell types. Biochim Biophys Acta. 2008, 1778: 2378-2387. 10.1016/j.bbamem.2008.06.010.
Choe SS, Choi AH, Lee JW, Kim KH, Chung JJ, Park J, Lee KM, Park KG, Lee IK, Kim JB: Chronic activation of liver X receptor induces beta-cell apoptosis through hyperactivation of lipogenesis: liver X receptor-mediated lipotoxicity in pancreatic beta-cells. Diabetes. 2007, 56: 1534-1543. 10.2337/db06-1059.
Mehrotra A, Kaul D, Joshi K: LXR-alpha selectively reprogrammes cancer cells to enter into apoptosis. Mol Cell Biochem. 2011, 349: 41-55. 10.1007/s11010-010-0659-3.
Pommier AJ, Alves G, Viennois E, Bernard S, Communal Y, Sion B, Marceau G, Damon C, Mouzat K, Caira F, Baron S, Lobaccaro JM: Liver X Receptor activation downregulates AKT survival signaling in lipid rafts and induces apoptosis of prostate cancer cells. Oncogene. 2010, 29: 2712-2723. 10.1038/onc.2010.30.
Toresson G, Schuster GU, Steffensen KR, Bengtsson M, Ljunggren J, Dahlman-Wright K, Gustafsson JA: Purification of functional full-length liver X receptor beta produced in Escherichia coli. Protein Expr Purif. 2004, 35: 190-198. 10.1016/j.pep.2004.01.007.
Pehkonen P, Heikkinen S, Carlberg C: Analysis of genome-wide data on transcription factor binding and their target gene regulation: a case study on liver X receptors. Proc Int Workshop Comp Sys Biol. 2011, 2011: 149-152.
Salmon-Divon M, Dvinge H, Tammoja K, Bertone P: PeakAnalyzer: genome-wide annotation of chromatin binding and modification loci. BMC Bioinformatics. 2010, 11: 415-10.1186/1471-2105-11-415.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12: 996-1006.
R-Development-Core-Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria, ISBN 2009, 3-900051-900007-900050
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article 3
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Picard F, Robin S, Lavielle M, Vaisse C, Daudin JJ: A statistical approach for array CGH data analysis. BMC Bioinformatics. 2005, 6: 27-10.1186/1471-2105-6-27.
Picard F, Robin S, Lebarbier E, Daudin JJ: A segmentation/clustering model for the analysis of array CGH data. Biometrics. 2007, 63: 758-766. 10.1111/j.1541-0420.2006.00729.x.
Alberti S, Schuster G, Parini P, Feltkamp D, Diczfalusy U, Rudling M, Angelin B, Björkhem I, Pettersson S, Gustafsson JA: Hepatic cholesterol metabolism and resistance to dietary cholesterol in LXRbeta-deficient mice. J Clin Invest. 2001, 107: 565-573. 10.1172/JCI9794.
This work was supported by the Academy of Finland, the Juselius Foundation, Télévie (Luxembourg), the Luxembourgish National Research Fund (FNR) and the University of Luxembourg (all to CC). LW-S had a FLARE fellowship (FNR/INTER/07/01) and E.T and K.R.S (522-2008-3745) were supported by the Swedish Research Council.
Kind thanks to the members of the Genomics Core Facility of the EMBL for massive parallel sequencing service and the Finnish Microarray and Sequencing Centre in Turku, Finland for microarray services.
PP, LW-S and CC designed the project, LW-S, JD and JR performed the experiments, PP, AW-B and SH analyzed the data, ET and KRS provided material and participated in the discussion and PP and CC wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1. Validation of LXR antibody. Protein extracts from livers of WT, LXRα-/-, LXRβ-/- and LXRαβ-/- mice  (A) or from HeLa cells, which were transfected with the empty pSG5 expression vector as control and pSG5 expressing human LXRα or LXRβ, respectively, using FuGENE 6 transfection reagent (Roche) (B), were separated on 8% sodium dodecyl sulfate (SDS) polyacrylamide gels. Western blotting using anti-LXR antibody was then performed using standard procedures. (PDF 81 KB)
Additional file 2: Table S1. LXR binding locations in a human macrophage-type cell line. Columns indicate the chromosome, start and end locations of ChIP-Seq peaks, peak length and summit, P-values from Poisson distribution, FE in comparison to IgG, FDR value and information on which sample peak was detected (T09- or vehicle-treated). Information regarding to any peak present in both samples is presented in two separate rows. (XLS 1 MB)
Additional file 3: Figure S2. Distribution of genomic LXR binding sites to genomic elements. A Distribution of LXR binding locations from the high stringent set of peaks to different genomic elements. Unique peaks represent binding locations present only in one of the two samples (disjoint areas of Venn diagram in Figure 1A) and in both samples (joint set of Venn diagram in Figure 1A). B Distribution of peaks around TSSs of closest genes. (PDF 1 MB)
Additional file 4: Figure S3. RE types present within the 202 high stringency LXR peak set. Proportions of high stringency set of peaks containing direct repeats (DRs), everted repeats (ERs) and inverted repeats (IRs). Search has been made with RSAT DNA-pattern tool  using RGKTCA half-site with indicated number of spacings. (PDF 146 KB)
Additional file 5: Table S2. Differentially expressed genes in a human macrophage-type cell line. Columns indicate the chromosome, start of the T09 target gene, log2 expression after T09- and vehicle-treatment, adjusted P-value for gene expression difference, gene name and FC. The detected LXR binding locations within proximal (± 100 kb) and distal (± 1 Mb) area from the respective gene TSS are displayed on the right side of the table. (XLS 378 KB)
Additional file 6: Table S3. Detailed information on all 112 LXR hotspots. Numbers of peaks and genes are indicated similarly as in Table 1. Additional columns indicate the tendency of each region having vast majority (> 2 times) of peaks in T09- (or vehicle-) treated sample versus vehicle- (or T09-) treated sample, the tendency of each region having majority (> 2 times) of up- (or down-) regulated genes versus down- (up-) regulated genes. Moreover, in separate columns up- and down-regulated genes, number of all genes, gene density and enriched GO terms (P < 0.05) are shown. Regions already shown in Table 1 are in bold. (XLS 69 KB)
Additional file 7: Figure S4. Enrichment statistics of LXR binding locations. The number of LXR peaks in one of the 112 hotspot regions (Figure 2) is compared with the number of all genes (A), the number of DE genes with adjusted P < 0.01 (B), the proportion of DE genes (C) and the density of DE genes (D). (PDF 217 KB)
Additional file 8: Figure S5. LXR target gene validation. PMA-differentiated THP-1 cells were treated for 4 h with vehicle (DMSO), 1 μM GW3965 (GW) or 1 μM T0901317 (T09), total RNA was extracted and qPCR was performed with primers specific for selected genes. The data were normalized to the expression of the housekeeping gene RPLP0 and fold inductions were calculated in reference to vehicle control. Columns indicate the means of four independent cell treatments and the bars represent standard deviations. Student's t-test was performed to determine the significance of the stimulation in reference to vehicle-treated control (** P < 0.01; *** P < 0.001). (PDF 41 KB)
Additional file 9: Table S4. Enrichment analysis for all LXR target genes. The analysis was performed using the DAVID tool  for all 1063 LXR target genes with P < 0.01 for differential expression. The columns indicate GO identifier, term name, count of associated target genes and the P-value and FDR for the enrichment. (XLS 56 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Pehkonen, P., Welter-Stahl, L., Diwo, J. et al. Genome-wide landscape of liver X receptor chromatin binding and gene regulation in human macrophages. BMC Genomics 13, 50 (2012) doi:10.1186/1471-2164-13-50
- False Discovery Rate
- Differentially Express Gene
- Regulatory Sequence Analysis Tool