|
|
||||||||
Papers In Press, published online ahead of print November 1, 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Journal of Lipid Research, Vol. 44, 2109-2119, November 2003
Copyright © 2003 by American Society for Biochemistry and Molecular Biology
Laboratory of Biochemical Genetics and Metabolism, The Rockefeller University, 1230 York Avenue, Box 179, New York, NY 10021
The online version of this article (available at http://www.jlr.org) contains one supplemental table. ![]()
Published, JLR Papers in Press, August 1, 2003. DOI 10.1194/jlr.M300203-JLR200
1 To whom correspondence should be addressed. e-mail: breslow{at}rockefeller.edu
| ABSTRACT |
|---|
|
|
|---|
In summary, these studies identified and confirmed six novel dietary cholesterol-regulated genes, three putative SREBP target genes (calcium/calmodulin-dependent protein kinase 1D, fatty acid binding protein 5, and proprotein convertase subtilisin/kexin 9), and three putative LXR target genes (a disintegrin and metalloprotease domain 11, apoptosis-inhibitory 6, and F-box-only protein 3).
Abbreviations: Abcg5, ATP binding cassette, subfamily G, member 5; Acac, acetyl-CoA carboxylase; Acly, ATP citrate lyase; Adam11, a disintegrin and metalloprotease domain 11; Api6, apoptosis-inhibitory 6; Camk1D, calcium/calmodulin-dependent protein kinase 1D; CYP7A1, cholesterol 7-
-hydroxylase; FABP5, fatty acid binding protein 5; Fbxo3, F-box-only protein 3; HMGCR, HMG-CoA reductase; HMGCS, HMG-CoA synthase; Laptm5, lysosomal-associated protein transmembrane 5; LXR, liver X receptor; Pcsk9, proprotein convertase subtilisin/kexin 9; Rgs16, regulator of G-protein signaling 16; SAA3, serum amyloid A 3; Sqle, squalene epoxidase; SREBP, sterol-regulatory element binding protein; StarD4, StAR-related lipid transfer domain containing 4
Supplementary key words dietary cholesterol oligonucleotide arrays gene expression profiling sterol-regulatory element binding proteins liver X receptor
| INTRODUCTION |
|---|
|
|
|---|
-hydroxylase (CYP7A1) (1, 2). Cholesterol can regulate expression at both transcriptional and posttranscriptional levels, and in the former, two key transcription factors have been implicated, sterol-regulatory element binding protein (SREBP) and the liver X receptor (LXR) (3, 4). The SREBP pathway utilizes three transcription factors, SREBP-1a, -1c, and -2, which are synthesized as transmembrane precursors in the endoplasmic reticulum with an N-terminal basic helix-loop-helix-Zip transcription factor domain facing the cytoplasm. When cellular sterols are low, the SREBP cleavage-activating protein (SCAP)-INSIG-1 complex brings the SREBP to the Golgi, where the site 1 protease cleaves intraluminally and the site 2 protease cleaves intramembraneously to release the N-terminal domain of the protein into the cytoplasm. This N-terminal domain translocates to the nucleus and activates gene transcription (57). Korn et al. (8) showed that the SREBPs are responsive to cellular sterol levels in vivo, as the amounts of nuclear SREBP-1 and SREBP-2 are decreased by a 1% cholesterol/0.5% cholic acid diet. Studies with mice overexpressing truncated nuclear forms of SREBP-1a, -1c, and -2 have revealed that SREBP-2 preferentially activates genes involved in cholesterol biosynthesis and metabolism, such as HMG-CoA synthase (HMGCS) and HMGCR, whereas SREBP-1c and -1a preferentially activate genes involved in fatty acid biosynthesis, such as acetyl-CoA carboxylase (Acac) and ATP citrate lyase (Acly) (912).
LXRs are nuclear hormone receptors that form obligate heterodimers with the retinoid X receptor. When bound to their oxysterol ligands, LXRs activate transcription of target genes (3). The importance of the LXR pathway was demonstrated in part by studies in knockout mice. When both wild-type and LXR
knockout mice were challenged with increased dietary cholesterol, LXR
knockout mice showed a greater increase in liver cholesterol, increased plasma cholesterol, blunted bile acid excretion, and the inability to upregulate CYP7A1, the rate-limiting step in bile acid production (13). It has also been shown that LXRs control genes involved in sterol movement across membranes, including ATP binding cassettes G5 and G8 (Abcg5 and Abcg8) and ABCA1, in lipoprotein metabolism, including cholesterol ester transfer protein, and in fatty acid synthesis, including SREBP-1c, among others (14).
While much information about cholesterol-responsive genes has been obtained from in vitro work and studies with genetically modified mice, it is important to understand the mechanisms an intact organism uses to respond to increased levels of dietary cholesterol. This study aims to identify genes in the mouse liver that are important in this regard by using microarrays. Affymetrix oligonucleotide microarrays were interrogated with liver RNA samples from mice fed 0.0% versus 0.5% cholesterol diets. The combined analysis of three experiments identified 69 genes regulated by dietary cholesterol, 37 downregulated and 32 upregulated. Regulation of a subset of the identified array genes was confirmed using quantitative RT-PCR (Q-PCR) analysis. Further studies are presented regarding the time course of regulation by dietary cholesterol and the role of the SREBP and LXR transcription factors. These studies have allowed the identification of three novel putative SREBP and three novel putative LXR target genes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Plasma and liver cholesterol measurements
Plasma was separated from whole blood by centrifugation. Plasma total cholesterol (mg/dl) was measured enzymatically using the Infinity Cholesterol Reagent (Sigma Diagnostics, 402-100P). Liver total, free, and esterified cholesterol (mg/g liver) were measured by gas chromatography with coprostanol as an internal standard as previously described (16).
Sample preparation for gene expression analysis
Liver tissue in RNA later was homogenized in TRIzol reagent (Invitrogen), and total RNA was isolated according to the manufacturer's instructions. Total RNA was then subjected to RNeasy Cleanup (Qiagen) for microarray and Q-PCR. For microarrays, 20 µg of total RNA was reverse transcribed using Superscript II (Invitrogen) and a poly-dT primer containing the T7 RNA polymerase binding site (Genset Corporation). Second-strand cDNA was then made using Escherichia coli DNA polymerase, E. coli DNA Ligase, and T4 DNA polymerase (Invitrogen). Double-stranded cDNA was purified using phenol-chloroform-isoamyl alcohol on Phase-Lock gel columns (Eppendorf), ethanol precipitated, and resuspended in water. cRNA was synthesized from the cDNA using biotin-labeled ribonucleotides and T7 RNA polymerase (Enzo Bioarray/Affymetrix) and purified on Qiagen RNeasy columns. cRNA was then fragmented for 30 min at 95°C in a solution of 40 mM Tris-acetate pH 8.1, 100 mM potassium acetate, and 30 mM magnesium acetate. For Q-PCR analysis, total RNA was treated with Dnase I (Ambion), and 5 µg was reverse transcribed using Superscript II and a mixture of oligo-dT and random hexamer primers.
Affymetrix oligonucleotide microarrays
Three microarray experiments were performed using cholesterol-fed and control mice according to the protocol above. Six male mice were used for experiment 1 (n = 3 on each diet), and 10 male and 10 female mice (n = 5 on each diet) were used for experiment 2. In these experiments, mice were housed together in cages. Equal amounts of liver RNA samples from mice within a sex and feeding group were pooled together and split into duplicate samples for cDNA and cRNA synthesis and application to the microarrays. In experiment 3, eight male and eight female mice (n = 4 on each diet) were housed in separate metabolic cages to prevent coprophagia and competition for food, and liver RNA from each mouse was processed for an individual microarray. For all three experiments, fragmented cRNA samples were combined with hybridization controls (Affymetrix), herring sperm DNA (Promega), and acetylated BSA (Invitrogen) and hybridized to microarrays for 16 h at 45°C. Each sample was hybridized to the Affymetrix Test3 microarray and the MGU74v2 microarray set (three chips: A, B, and C). This set contains >36,000 probe sets corresponding to >23,000 Unigene clusters (for information, see www.affymetrix.com). Washing, staining, and scanning of the microarrays was performed by the Gene Array Core Facility at Rockefeller University according to standard protocols (Affymetrix).
Data analysis of Affymetrix microarrays
The raw data from the microarrays was analyzed using the Microarray Suite Software version 5.0 (MAS5.0, Affymetrix). Each array was subjected to absolute analysis, and all chips showed proper behavior of spiked controls, signal ratios of 5' to 3' sequences from housekeeping genes approximately equal to 1, and low background and noise values. The signal values for all arrays were scaled using a scaling factor (SF) derived from the formula: SF x mean signal = 250. For experiments 1 and 2, there were duplicate arrays for each diet, so four pair-wise comparisons were performed between feeding groups within a sex. In both experiments, genes with a signal in at least one condition >250 and average fold for the pair-wise comparisons >|1.6| were included in the preliminary analysis. As experiment 3 had four mice per group, this allowed examination of biological variation in gene expression, and statistical analysis was performed to determine significant differences between the groups. From an initial power calculation, utilizing the average variance for all genes on the MGU74v2 chip set in the 0.0% cholesterol male group, it was determined that with four replicates, a fold >|1.6| could be reliably discerned if genes with a signal >250 were used. Therefore, the experiment 3 signals derived from the MAS5.0 absolute analysis were subjected to an ANOVA statistical analysis software package (17). Genes in experiment 3 were determined to be significantly regulated by dietary cholesterol if the signals in at least one group were >250, the fold >|1.6|, and the ANOVA P value < 0.001. Genes were also kept in the final analysis if in experiment 3 their ANOVA P value was between 0.01 and 0.001, and they also met the criteria specified above for the preliminary analysis in experiments 1 or 2. Information regarding each probe set identified in the combined statistical and MAS analysis was downloaded from www.affymetrix.com. This information included the gene names shown in Tables 1 and 2. Other information such as accession number and Unigene cluster can be obtained by searching with the probe set identifier. In cases where the probe set was not assigned to a gene by Affymetrix, the sequence was used to query the mouse genome (www.ensembl.org) and identify a putative gene. This information is in parentheses next to the Affymetrix name in Tables 1 and 2.
|
|
Cloning of calcium/calmodulin-dependent protein kinase 1D and proprotein convertase subtilisin/kexin 9
At the time these studies were performed, two genes of interest had not been cloned. The Affymetrix probe set sequence was used to query the mouse genome, and putative coding sequences were identified and confirmed by RT-PCR from mouse liver. PCR reagents were Advantage cDNA polymerase (Clontech) and Super Taq Plus polymerase (Ambion), primers were from GeneLink (sequences in the supplementary data), and thermal cycling was on a Perkin-Elmer 9700. The 5' ends of both genes were identified by 5' rapid amplification of cDNA ends (RACE) (Ambion 5' RLM-RACE), and the 3' end was deduced by identification of a poly-A signal (ATTAAA) in the genome sequence in agreement with transcript size on Northern blotting. PCR products were TA cloned with pCR-2.1-TOPO (Invitrogen) and the sequence verified. Sequences have been deposited in GenBank [Accession Numbers: proprotein convertase subtilisin/kexin 9 (Pcsk9)-AY273821, and calcium/calmodulin-dependent protein kinase 1D (Camk1D)-AY273822].
| RESULTS |
|---|
|
|
|---|
|
The design of experiment 3 allowed a statistical analysis of the microarray results. Using an ANOVA statistical analysis software package (17), genes were determined to be significantly regulated by the 0.5% cholesterol diet if the signals in at least one group were >250, fold >|1.6|, and with an ANOVA P value < 0.001. Genes were also kept in the analysis if the ANOVA P value in experiment 3 was between 0.01 and 0.001, and if they also met the criteria for MAS5.0 comparison analysis in experiments 1 or 2. This more stringent analysis yielded 48 downregulated probe sets and 32 upregulated probe sets that showed regulation by dietary cholesterol in males and/or females.
Genes regulated by dietary cholesterol
Tables 1 and 2 show the 48 probe sets representing 37 Unigene clusters downregulated by dietary cholesterol and the 32 probe sets representing 32 Unigene clusters upregulated by dietary cholesterol in this study. The tables also show P values for experiment 3 as derived from the ANOVA analysis and the average folds for all pair-wise comparisons from each experiment as derived from the MAS5.0 software. The cholesterol-regulated genes were classified according to the categories specified by Affymetrix.
Among the 37 genes downregulated by dietary cholesterol (Table 1), 19 were classified in cholesterol metabolism, fatty acid metabolism, or the metabolism of other lipids. Eleven of these have previously been shown to be SREBP targets (marked by superscript a). Squalene epoxidase (Sqle), isopentenyl diphosphate delta isomerase, farnesyl diphosphate farnesyl transferase, phosphomevalonate kinase, farnesyl diphosphate synthetase, and HMGCS are genes in the cholesterol biosynthetic pathway (1922). Malic enzyme and acetyl-CoA synthetase are involved in fatty acid synthesis (23, 24). Acyl-CoA binding protein is implicated in acyl-CoA transport and steroidogenesis (25). Delta-5 desaturase is necessary for the conversion of dietary linoleic acid (
3 series) and
-linolenic acid (
6 series) to the biologically active arachidonic acid, eicosopentenoic acid, and docosahexaenoic acid (26). CTP-choline phosphate cytidyltransferase is involved in phosphatidylcholine biosynthesis (27). The other eight genes in the metabolism categories, as well as the remaining 18 genes in the miscellaneous and unclassified categories, are not previously known SREBP targets.
Among the 32 genes upregulated by dietary cholesterol (Table 2), 13 genes were classified in metabolism. Two of these have previously been shown to be LXR targets (marked by superscript a). Abcg5 is involved in sterol transport in the intestine and liver (28). ABCA1 mediates cholesterol efflux, and its mRNA is highly regulated by LXR in the macrophage (29). Although in one previous study, liver ABCA1 mRNA levels were not increased by an LXR agonist (30), another study showed increased liver protein levels in mice fed a high-fat/high-cholesterol diet (31). In the current study, liver ABCA1 mRNA levels were increased by dietary cholesterol only in females. The other 11 genes in the metabolism category and the remaining 19 genes in the other categories (seven in immune/acute phase response, two as transcription factors, and ten unclassified) are not previously known LXR targets.
Confirmation of regulation by RT-PCR analysis for selected array genes
To confirm the results of the microarray study, a selection of genes that appeared to be highly regulated and/or possessed interesting domains suggestive of function were subjected to Q-PCR analysis, along with control SREBP and LXR targets (Table 3). For the downregulated genes, this analysis was performed for two known SREBP-1 target genes (Acac and Acly), three known SREBP-2 target genes (HMGCR, HMGCS, and Sqle), and six other genes. Acac was downregulated 2.3-fold in males only, and Acly was downregulated 2.9-fold in males only. HMGCR was downregulated 7.9-fold in males only, HMGCS was downregulated 4.5- and 2.3-fold in males and females, respectively, and Sqle was downregulated 38.6- and 27.8-fold in males and females, respectively. The six other genes were not known SREBP targets and four were previously cloned. Fatty acid binding protein 5 (FABP5) can bind long-chain fatty acids, but its function is unknown (32). FABP5 was downregulated 6.0- and 3.1-fold in males and females, respectively. Regulator of G protein signaling 16 (Rgs16) is a GTPase-activating protein that attenuates heterotrimeric G-protein signaling (33). Rgs16 was downregulated 1.8-fold in males and <1.6-fold in females, yet this did not reach significance in either case (data not shown). Ubiquitously expressed transcript, which may be a transcription factor, was not significantly regulated in males and females (data not shown). StAR-related lipid transfer domain containing 4 (StarD4), a recently cloned gene putatively involved in intracellular cholesterol transport, has previously been shown by Q-PCR analysis to be downregulated by dietary cholesterol in females (18). This was confirmed in the current study for males, in which there was 2.3-fold downregulation.
|
To confirm the microarray results for the upregulated genes, this analysis was performed for one known LXR target gene (Abcg5) and five other genes. Abcg5 was upregulated 4.0- and 2.2-fold in males and females, respectively. The five other genes were not known LXR targets, and all were previously cloned. F-box-only protein 3 (Fbxo3) was upregulated 2.0-fold in females only (Table 3). Lysosomal-associated transmembrane protein 5 was upregulated 1.8-fold in males only, although this was not significant (data not shown). Serum amyloid A 3 (SAA3), an acute phase reactant, was upregulated 11.4- and 6.0-fold in males and females, respectively. Apoptosis-inhibitory 6 (Api6) is a member of the scavenger receptor family (36) and was upregulated 1.7- and 1.8-fold in males and females, respectively. Finally, a disintegrin and metalloprotease domain 11 (Adam11) family member (37) was the most highly upregulated gene in this study. It was upregulated 11.7- and 10.5-fold in males and females, respectively.
Temporal patterns of gene expression after dietary cholesterol feeding
The effect of dietary cholesterol on gene expression over time was determined by Q-PCR analysis for each of the confirmed genes above except StarD4. The downregulated genes showed two temporal patterns of regulation by dietary cholesterol, as shown in Fig. 2A
. HMGCR, HMGCS, Pcsk9, and Sqle were downregulated 7090% after 1 day of cholesterol feeding and remained down through day 7. In contrast, Acly, Camk1D, and FABP5 were downregulated only 2040% after 1 day of cholesterol feeding, trended toward upregulation on day 2, and then were further downregulated between 5075% through day 7.
|
Expression of downregulated genes in SREBP transgenic mice
To confirm that the selected downregulated genes discussed above were SREBP targets, Q-PCR analysis was performed on liver RNA from transgenic mice overexpressing the transactivation domain of SREBP-1a and SREBP-2 (9, 11). In these mice, SREBP-1a or SREBP-2 activates transcription regardless of cellular cholesterol levels, and target genes should be upregulated relative to levels in littermate controls. In this experiment, the known SREBP-2 target genes, HMGCR, HMGCS, and Sqle, were significantly upregulated in both male and female SREBP-1a and SREBP-2 transgenic mice compared with littermate controls, as expected (Table 4). A similar pattern was observed for Pcsk9, which was significantly upregulated in male and female SREBP-1a transgenic mice by 14.7- and 6.1-fold and in SREBP-2 transgenic mice by 6.4- and 4.0-fold, respectively. The known SREBP-1 target genes, Acac and Acly, were highly upregulated in male and female SREBP-1a mice and marginally or not at all in SREBP-2 mice, as expected. A similar pattern was observed for FABP5, which was significantly upregulated in male and female SREBP-1a transgenic mice by 7.3- and 26.0-fold, respectively, and unchanged in SREBP-2 transgenic mice. Similarly, Camk1D was significantly upregulated in male and female SREBP-1a transgenic mice by 5.6- and 3.1-fold, respectively, and unchanged in SREBP-2 transgenic mice.
|
As shown in Table 2, nine of the 13 metabolism genes upregulated by dietary cholesterol were also upregulated by treatment with the LXR agonist. None of the genes classified as immune/acute phase response genes that were upregulated by dietary cholesterol were upregulated by the LXR agonist. Finally, 6 of the 10 unclassified genes were upregulated by the LXR agonist. Ten other cholesterol upregulated genes were not found as LXR target genes for several reasons: five were on the C chip and not evaluated, three were called absent in the LXR experiment, and two were unchanged. Because only two of the upregulated genes, Abcg5 and ABCA1, are previously demonstrated LXR targets, this experiment suggests 13 additional LXR-regulated genes. Three of these, Adam11, Api6, and Fbxo3, were tested by Q-PCR analysis and confirmed (data not shown).
| DISCUSSION |
|---|
|
|
|---|
8,900 hybridized to liver RNA (24%) (data not shown). Thus, <1% of the probe sets revealed consistent regulation by dietary cholesterol. The expression changes of 15 of the 69 unique genes regulated by dietary cholesterol were examined by Q-PCR analysis and confirmed for 12 genes. The expression of this subset of genes was further examined in a time course of cholesterol feeding and in SREBP transgenic mice by Q-PCR. In addition, the expression of all of the genes was further examined by Affymetrix microarray analysis in mice treated with the LXR agonist. The combined analysis yielded three novel putative SREBP target genes, Camk1D, FABP5, and Pcsk9, and 13 novel putative LXR target genes in mouse liver, three of which were confirmed by Q-PCR, Adam11, Api6, and Fbxo3. This study has two obvious strengths. The first strength was the model we used, which was based on a physiological paradigm, namely cholesterol feeding of an intact animal, whereas previous studies have used cell culture models or genetically manipulated mice, such as SREBP transgenic and LXR knockout mice. The second strength was performing three independent experiments, especially the third experiment, which examined responses of four mice per condition to assess the individual variation and to perform statistical analysis of significant differences between groups. Many microarray analyses rely on a single experiment using one microarray per condition. In the current study, compared with the 48 decreased and 32 increased probe sets listed for the combined analysis, had we taken only the first experiment and used one A chip per condition (0.5% vs. 0.0% cholesterol), the Microarray Suite software would have listed 82 decreased and 150 increased probe sets. Of the 82 listed as decreased, only 18 were on the final list, and of the 150 listed as increased, only 12 were on the final list. It is difficult to assess the sensitivity and specificity for regulated genes identified by microarray analysis for most studies in the literature. In the current study, despite the conservative nature in which genes were selected, three of 15 genes tested by Q-PCR were not confirmed.
This study also has potential disadvantages. First, many studies of regulated gene expression have used cell culture systems, and the current use of an in vivo model clearly is more complex and subject to greater variability. However, the use of an inbred mouse strain and strictly controlled environmental conditions minimizes this concern to the extent possible in a complex living organism. Second, C57BL/6J mice do not respond to cholesterol feeding with an increase in plasma cholesterol levels. However, we show in the current study that dietary cholesterol increases liver total, free, and esterified cholesterol, thus validating our approach for this tissue. Third, a number of well-characterized cholesterol-regulated genes were not identified by our statistical criteria as regulated by 0.5% dietary cholesterol for 1 week. This could be due to the stringency of the qualifying criteria and/or the experimental design. For example, HMGCR and 7-dehydrocholesterol reductase, two known SREBP target genes (20), did not make our list, even though they trended toward being decreased by cholesterol feeding, as they showed extensive array signal variability between mice and did not reach statistical significance. As another example, mevalonate kinase and diphosphomevalonate decarboxylase, also SREBP target genes (20), had very low signal values, presumably due to poor probe set design, and thus were not considered. Finally, the known SREBP target gene LDL receptor (5) and the known LXR target gene CYP7A1 (3) were not regulated by our feeding paradigm as analyzed by microarray or Q-PCR, although they have been shown to be regulated by higher concentrations of dietary cholesterol and by longer feeding times (13, 41). Perhaps if we had used a different concentration of cholesterol or length of feeding, these and other well-characterized cholesterol regulated genes would have been regulated on our microarrays.
In our experimental design, gene expression was analyzed after mice had been fed cholesterol for 1 week. To assess which genes were early and late responders, a subset of 12 responsive genes were examined at earlier time points. Genes involved in cholesterol biosynthesis and regulated by SREBP-2, such as HMGCR, HMGCS, and Sqle, were immediately downregulated at day 1 of feeding. As Pcsk9 also followed this pattern, we predict it may be involved in cholesterol biosynthesis. In contrast, a gene involved in fatty acid metabolism and regulated by SREBP-1, Acly, was downregulated later in the cholesterol feeding time course. Because Camk1D and FABP5 also followed this pattern, we predict they may be involved in fatty acid metabolism. There were several genes that were upregulated after 1 day of cholesterol feeding, such as Abcg5, Adam11, Api6, and Fbxo3. As Abcg5 is involved in reverse cholesterol transport, we speculate that the other three genes may also be important for this process. SAA3 was not induced until a later time point and probably represents a secondary response to the high-cholesterol diet, such as an inflammatory response.
Integration of information derived from the cholesterol feeding study, the LXR agonist study, and the SREBP-1a and SREBP-2 transgenic experiments also allows us to make predictions about the function of some of the novel genes we have identified that were downregulated by dietary cholesterol. For example, cholesterologenic SREBP-2 target genes should have the following profile: 1) downregulated by dietary cholesterol due to the inhibitory effect of sterols on SREBP processing, 2) upregulated in SREBP-2 transgenic mice and SREBP-1a transgenic mice, as SREBP-1a is a potent activator of all SREBP target genes (6), and 3) unchanged by LXR agonist feeding, because LXR transcriptionally activates only SREBP-1c (39). In our study, this has been shown for HMGCR, HMGCS, and Sqle, three known cholesterol biosynthetic genes (see Tables 3 and 4). This is supported by the behavior of known SREBP-2 target genes in our dietary cholesterol and LXR agonist feeding studies (see cholesterol metabolism genes marked by superscript a in Table 1). Therefore, we propose that other genes downregulated by dietary cholesterol and unchanged by the LXR agonist are also cholesterologenic genes regulated by SREBP-2 (Table 1). Some examples include sterol C5 desaturase, sterol C4 methyl oxidase-like, NADP-dependent steroid dehydrogenase-like, StarD4, TM7SF2/lamin B receptor-related and androgen-regulated short-chain dehydrogenase/reductase 1.
In contrast, lipogenic SREBP-1c target genes should have the following profile: 1) downregulated by dietary cholesterol, 2) upregulated in SREBP-1a transgenic mice only, and 3) upregulated by LXR agonist feeding. This has been shown in our study for Acac and Acly, known lipogenic genes (Tables 3 and 4). This is supported by the behavior of known SREBP-1 target genes in our dietary cholesterol and LXR agonist feeding studies (see fatty acid metabolism genes marked by superscript a in Table 1). Therefore, we propose that other genes downregulated by dietary cholesterol and upregulated by the LXR agonist are also lipogenic genes regulated by SREBP-1 (Table 1). An interesting example is FABP5, which we have shown to be regulated by SREBP-1a and not SREBP-2 in this study. FABP5 is a member of the fatty acid binding protein family and is expressed ubiquitously. It has been shown to bind long-chain fatty acids, but its physiological function is unknown (32, 42).
Two of the novel dietary cholesterol downregulated genes cloned here, the kinase Camk1D and the protease Pcsk9, do not fit neatly into this analysis. Camk1D expression was modestly downregulated by dietary cholesterol only in males, upregulated in SREBP-1a but not SREBP-2 transgenic male and female mice, but unchanged by LXR agonist feeding in males. Pcsk9 expression was highly downregulated by dietary cholesterol, highly upregulated in both SREBP-1a and SREBP-2 transgenic mice, but modestly upregulated by LXR agonist feeding. Thus, we cannot apply our algorithm exactly to these genes, but Camk1D expression resembles more a lipogenic gene and Pcsk9 a cholesterologenic gene. Mouse Camk1D has a human homolog, but nothing is known about its function. During review of this manuscript, Abifadel et al. (43) reported that mutations in Pcsk9 are present in all affected members of two families with autosomal dominant hypercholesterolemia (OMIM 603776), but absent in 200 control chromosomes. This strongly implicates mutations in Pcsk9 as causative for some forms of familial hypercholesterolemia. Their findings and our observation of strong regulation by dietary cholesterol imply an important role for Pcsk9 in cholesterol metabolism. This also provides exciting validation of our experimental approach to identifying novel genes important in cholesterol metabolism.
We can also predict the function of some novel genes upregulated by dietary cholesterol. For example, genes involved in reverse cholesterol transport should be upregulated by feeding either cholesterol or the LXR agonist (such genes should not be SREBP targets and were not examined in the SREBP-1a or SREBP-2 transgenic mice). We observe this pattern for ABCA1 and Abcg5, two genes with established roles in reverse cholesterol transport. Therefore, we propose that other genes upregulated by dietary cholesterol and upregulated by the LXR agonist may also be involved in reverse cholesterol transport (Table 2). Interesting examples include the cytochrome P450 3a mouse genes (Cyp3a11 and 3a16), which have been shown to be involved in the formation of lithocholic acid (44), and Api6, a member of the scavenger receptor family (36). The highest upregulated gene in the dietary cholesterol microarray, Adam11, was also very highly induced by the LXR agonist. This gene is a member of the disintegrin and metalloprotease family with an inactive metalloprotease domain (37). ADAMs are membrane bound and some interact with integrins (45). The physiological function of Adam11 remains to be determined.
This study combined the use of microarrays to identify liver genes regulated by dietary cholesterol in mice with investigation of the transcription factor pathways that may be responsible for the effects. In this way, we have identified three putative SREBP targets and three putative LXR targets. Further studies of gene promoters are necessary to confirm that these are indeed SREBP and LXR targets. Additional cell culture and in vivo studies will be necessary to reveal the physiological role these genes play in cholesterol metabolism.
| ACKNOWLEDGMENTS |
|---|
Manuscript received May 15, 2003 and in revised form July 25, 2003.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||