|
|
||||||||
Papers In Press, published online ahead of print June 1, 2005
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



* Center for Molecular and Human Genetics, Columbus Children's Research Institute, Ohio State University, Columbus, OH
Division of Human Cancer Genetics, Comprehensive Cancer Center, Ohio State University, Columbus, OH
Department of Pediatrics, Ohio State University, Columbus, OH
The online version of this article (available at http://www.jlr.org) contains an additional three tables. ![]()
Published, JLR Papers in Press, April 1, 2005. DOI 10.1194/jlr.M400462-JLR200
1 To whom correspondence should be addressed. e-mail: hermang{at}pediatrics.ohio-state.edu
| ABSTRACT |
|---|
|
|
|---|
Of the 21 remaining genes, 16 are known genes, some of which have been implicated previously in cholesterol homeostasis or lipid-mediated signaling, and 5 are uncharacterized cDNA clones.
Abbreviations: Bpa, bare patches; CHILD, congenital hemidysplasia with ichthyosiform nevus and limb defects; 7DHC, 7-dehydrocholesterol; E13.5, embryonic day 13.5; ER, endoplasmic reticulum; FACS, fluorescence-activated cell sorting; GFP, green fluorescent protein; HMGR, 3-hydroxy-3-methylglutaryl-CoA reductase; LDS, lipid-depleted serum; MEF, mouse embryonic fibroblast; NSDHL, NADH sterol dehydrogenase-like; SCAP, sterol-regulatory element binding protein cleavage-activating protein; SREBP, sterol-regulatory element binding protein; wt, wild type
Supplementary key words cholesterol biosynthesis congenital hemidysplasia with ichthyosiform nevus and limb defects syndrome microarray bare patches sterol development
| INTRODUCTION |
|---|
|
|
|---|
30 enzymatic reactions that constitute the cholesterol biosynthetic pathway. Primary regulation of the pathway occurs by sterol-mediated feedback of transcription of the genes encoding cholesterogenic enzymes (reviewed in 1). Posttranscriptional mechanisms have also been described, including the regulated degradation in the endoplasmic reticulum (ER) of HMG-CoA reductase (HMGR), which catalyzes the rate-limiting step in the pathway (2). Further control involves regulation of the cellular uptake, efflux, and sites of storage of cholesterol within the cell (3). The precision with which cholesterol homeostasis is achieved and maintained may be attributable to its potentially deleterious effects on the organism when levels are increased and the tremendous energy expenditure involved in its de novo production. Most important, however, cholesterol is an essential constituent of cell membranes. Localized regions of increased cholesterol in the plasma membrane help define lipid rafts and caveolae that function in intercellular signaling and endocytosis, respectively (4). Covalent modification of active hedgehog proteins by cholesterol to provide membrane anchoring is required for their proper signaling functions during embryonic development (5). Furthermore, cholesterol is the precursor for the synthesis of steroid hormones and bile acids.
During the past decade, seven human developmental malformation syndromes have been identified that arise from mutations in genes encoding enzymes that function in the postsqualene steps of cholesterol biosynthesis (6, 7). These disorders are marked by several common features, including growth and mental retardation, major developmental malformations, and skeletal defects. However, they are distinguished by their incidence and severity, following a general trend of decreasing incidence and increasing severity, including prenatal lethality, for genes that function in earlier steps of the pathway. It has been suggested that deficiency of cholesterol or total sterols, toxicity from the excessive accumulation of specific sterol intermediates, distinct effects on hedgehog protein modification and signaling, and abnormal feedback regulation of early steps in the pathway, with concomitant effects on isoprenoids that are generated as side products, all may contribute to the observed phenotypes and disease pathogenesis.
Spontaneous mouse mutants or animals generated by targeted gene disruption exist for six of the seven human disorders, providing potential models to examine disease pathogenesis (7, 8). One of these mouse models is the X-linked, male-lethal bare patches (Bpa) mouse resulting from mutations in the NADPH steroid dehydrogenase-like (Nsdhl) gene (9). Nsdhl encodes a ubiquitously expressed 362 amino acid protein that functions as a 3ß-hydroxysterol dehydrogenase and is involved in the removal of C-4 methyl groups at an intermediate step in the conversion of lanosterol to cholesterol. The gene was first identified by positional cloning in the mouse. Subsequently, mutations in human NSDHL were discovered in patients with congenital hemidysplasia with ichthyosiform nevus and limb defects (CHILD) syndrome, an extremely rare X-linked, male-lethal malformation syndrome associated with unilateral ichthyosiform skin lesions and limb reduction defects (10, 11).
All of the seven known mutant murine Nsdhl alleles produce a characteristic striping of the coat in heterozygous females that follows the lines of X inactivation, and all of the mutations are lethal by embryonic day 13.5 (E13.5) in hemizygous males (9, 12, 13). The Bpa1H allele, resulting from the nonsense mutation K103X, produces the most severe phenotype in surviving females, with asymmetric dwarfing and skeletal dysplasia, early postnatal patchy hyperkeratotic skin eruptions, and occasional microphthalmia and/or cataracts. The majority of hemizygous Bpa1H males die by E9.5. Evidence from sterol analyses of affected male embryos for two less severe Nsdhl alleles suggests that the male lethality is not attributable to simple cholesterol deficiency, because cholesterol and total sterol levels in affected male embryos are comparable to those of wild-type (wt) littermates, probably as a result of maternal transfer across the placenta (13).
To begin to examine the cellular consequences of an enzymatic block in cholesterol synthesis at the level of the sterol C-4 demethylase complex, we have isolated mouse embryonic fibroblasts (MEFs) from Bpa1H embryos. We describe here comparisons of gene expression profiles from mutant and normal MEFs using microarray and real-time PCR analyses.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Generation and culture of primary murine embryonic fibroblasts
To generate primary cultures of embryonic fibroblasts, Bpa1H females were mated to a backcross N5 GFP transgenic male. Pregnant females were sacrificed at E15.5, and yolk sac DNA was purified for PCR genotyping of the embryos. The sex of the embryos was determined using primers for Smcx as described previously (18), and Nsdhl genotype was inferred using the closely linked microsatellite DXMit1 (19).
All tissue culture reagents were purchased from Invitrogen Life Technologies (Carlsbad, CA) unless stated otherwise. After removing the head and internal organs, the embryos were minced and incubated for 30 min at 37°C in 1 mg/ml collagenase A (Roche Applied Science, Indianapolis, IN) dissolved in MEM
medium. The disaggregated cells were pelleted by centrifugation and plated at 5 x 105 cells/75 cm flask in MEM
supplemented with 10% FBS (Hy-Clone, Logan, UT), penicillin, streptomycin, and 2 mM L-glutamine. Cells were grown at 37°C in a 5% CO2 atmosphere and replated at 5 x 105 cells/flask every 3 days.
Lipid-depleted serum (LDS) was prepared by adding 2 mg of silicon dioxide (Cab-O-SilTM M-5; Acros, Morris Plains, NJ) to 100 ml of FBS and stirring overnight at 4°C (20). The serum was passed through a 0.2 µM sterile filter, and 15 mg of insulin-transferrin-sodium selenite media supplement (catalog number I 1884; Sigma, St. Louis, MO) was added per 100 ml of serum. The treated serum was 97% cholesterol-depleted as determined using the InfinityTM cholesterol assay (Thermo Electron Corp., Waltham, MA).
MEFs that expressed GFP were counted using a Coulter EpicsXL flow cytometer equipped with EXPO32 acquisition and analysis software (Beckman Coulter, Fullerton, CA). GFP-negative (GFP) cells were sorted from mixed populations of heterozygous female MEFs at passage 3 using a FACSVantage DiVa instrument (BD Biosciences, San Jose, CA). Sorted cells were plated at 5 x 105/75 cm flask and replated when they reached confluence after 4 days. Subconfluent cultures of MEFs at passage 5 were used for microarray analysis. For growth in LDS, cells were first grown in normal medium for 24 h after replating, then rinsed twice with PBS, and LDS medium was added. Incubation in LDS medium was continued at 37°C for 24, 36, or 48 h, at which time the cells were harvested for RNA isolation.
RT-PCR, sequencing, and microarrays
Total RNA was isolated from MEFs at passage 5 using TRIzol Reagent (Invitrogen, Carlsbad, C) according to the manufacturer's protocol. RNA was further purified using an RNeasy Mini Kit (Qiagen, Valencia, CA). The quality and concentration of RNA samples were assayed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). RT-PCR and DNA sequencing of cDNA to test for the monoallelic expression of Nsdhl in sorted cell populations were performed as described (21) using a forward primer from exon 2 (NsdhlE2F, 5'-CCAGCTGATCATGGTGAATC-3') and a reverse primer from exon 5 (NsdhlE5R, 5'-TGGCACTGCTGGTTAAAATGAGTTT-3'), generating a 400 bp PCR product that spans the site of the Bpa1H mutation in exon 4.
For microarray analysis, cDNA was synthesized from 40 µg of total RNA using amino-allyl dUTP (Sigma), and either Cy3 or Cy5 dye (Amersham Biosciences, Piscataway, NJ) was then conjugated to the cDNA samples (22). Approximately 0.5 µg of purified target cDNAs from mutant and wt cells was combined with control target cDNAs in 500 µl of hybridization buffer (In Situ Hybridization Kit Plus; Agilent Technologies) and hybridized to glass slide microarrays containing 60-mer oligonucleotides representing
22,000 mouse gene sequences (Mouse Development Oligo Microarray; Agilent Technologies). Hybridizations were performed with constant rotation in Agilent hybridization chambers for 16 h at 60°C. Slides were washed in 6x SSC (1x SSC is 150 mM sodium chloride, 15 mM sodium citrate) with 0.005% Triton X-102 for 10 min at room temperature and in ice-cold 0.1x SSC, 0.0005% Triton X-102 for 5 min. Fluorescence signal was quantified using an Affymetrix 428 Array Scanner (Affymetrix, Santa Clara, CA) and GenePix 4.0 software (Axon Instruments, Union City, CA).
A total of 15 microarray experiments were performed, each involving cohybridization of wt and Bpa1H cDNA samples from littermates. Six primary cultures of sorted GFP MEFs from three different litters were used and designated 1W, 1B, 2W, 2B, 3W, and 3B, where the number specifies the litter and W or B signifies MEFs expressing the wt or Bpa1H Nsdhl allele, respectively. Five experiments (four with 1W vs. 1B and one with 3W vs. 3B) were performed with RNA from cells grown in normal medium; four experiments (two using 1W vs. 1B, one using 2W vs. 2B, and one using 3W vs. 3B) included RNA from cells grown in LDS for 24 h; two experiments (one using 1W vs. 1B and one using 3W vs. 3B) included samples from cells grown in LDS for 36 h; and four experiments (three using 1W vs. 1B and one using 3W vs. 3B) contained RNA from cells grown in LDS for 48 h. Independently isolated RNA samples were used in all of the experiments except two, one in normal medium and one after 48 h in LDS, in which the same RNA sample pairs, 1W and 1B with reversed Cy3 and Cy5 labeling, were hybridized to different microarrays.
Statistical analysis of microarray data was performed using the software package S-Plus version 6.2.1 (Insightful Corp., Seattle, WA). Before statistical analysis, the median of the background corrected fluorescence intensities of the two dyes (Cy3 and Cy5), as well as the ratios of the median intensities of the two dyes, were calculated for each spot on the array. Array spots were screened for poor quality or low signal intensity compared with the local background (% > B635 + 2SD < 25 for the Cy5 channel and % > B532 + 2SD < 25 for the Cy3 channel), and those identified were excluded from the analysis and treated as missing values (denoted NA in supplementary Tables I and II Array data sets were then normalized using an intensity-dependent normalization method by fitting a locally weighted least-square fit (LOWESS) to the M versus A plot, where M = log intensity ratio log2 (Cy3/Cy5) and A = mean log intensity log2 (Cy3 x Cy5)1/2 (23). Data for color inversion experiments were combined and averaged as a single experiment. A one-sample t-test was performed to identify upregulated or downregulated genes in the Bpa1H cells compared with the wt cells for replicate experiments from each growth condition time point separately, as well as for the three LDS time points combined, under the null hypothesis that the mean log ratio of test over reference is zero. Cluster analysis was performed using GeneCluster software, and the results were displayed graphically using TREEVIEW version 1.60 software (24). Microarray probes that met the following criteria were included in the cluster analysis: they showed a mean fold difference between the Bpa1H and wt signal greater than 1.5 for at least one of the experimental time points with a P value of less than 0.1, and there were less than 80% missing values in the 13 experiments.
To identify the most consistent and statistically significant differences in gene expression between Bpa1H and wt cells associated with lipid depletion, data from the three LDS time point experiments (24, 36, and 48 h) were combined. Although the yield of RNA from samples from the 48 h time point was lower than from the other time points, there was no difference in the quality of the 48 h RNA samples, as determined by analysis on the Bioanalyzer. P values were assigned to genes that displayed a mean fold difference of greater than 1.5 and had less than 50% missing values for the nine LDS experiments. Lists of all of the upregulated and downregulated microarray probes, ranked by P value, with the calculated fold difference in individual experiments expressed as a log2 value, as well as the mean fold difference and P values calculated for the combined LDS experiments, are provided in supplementary Tables I and II.
Real-time PCR
cDNA was synthesized from 20 µg of total RNA from MEFs using SuperScript II reverse transcriptase with oligo-dT and random hexamer primers (Invitrogen) using standard protocols. After cDNA samples were brought to a final volume of 200 µl in water, a dilution series of 1:2, 1:4, and 1:8 was made for each sample, and aliquots were stored at 20°C. Amplification reactions included 2 µl of cDNA with 0.5 µM forward and reverse primers in SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA) in a total volume of 25 µl. Primers were designed to span introns if possible and to generate PCR products between 80 and 120 bp. Primer sequences are provided in supplementary Table III. Real-time PCR amplification was performed using an ABI Prism 7700 Sequence Detector equipped with Sequence Detector version 1.6.3 software (Perkin Elmer, Wellesley, MA) for data acquisition and analysis. Amplification reactions of the three cDNA dilutions (1:2, 1:4, and 1:8) were performed in parallel, and the results using the cDNA from wt MEFs grown in normal medium were used to generate a standard curve for each primer pair. The relative signal quantity for the remaining samples was calculated using this standard curve. RNA samples were analyzed from two Bpa1H (1B and 3B) and two wt (1W and 3W) cultures of sorted MEFs, each of which had been grown in either normal medium or LDS medium for 24 h. RNA aliquots from MEF cultures 3B and 3W were from the same samples used for microarray analysis, whereas RNAs from 1B and 1W were isolated from independent cultures of these cells.
| RESULTS |
|---|
|
|
|---|
The early lethality of affected Bpa1H male embryos, however, prohibited the expansion of sufficient viable MEFs for our experiments. We circumvented this problem by taking advantage of the phenomenon of random and stable X inactivation in mammalian females as follows: Bpa1H females were mated to males carrying a ubiquitously expressed X-linked GFP transgene [TgN(GFPX)4Nagy] that undergoes normal, random inactivation in the early postimplantation female embryo (14). Resulting female embryos all carried a paternally inherited X chromosome with a GFP transgene and wt Nsdhl allele, but they expressed both in approximately half of their cells as a result of random X inactivation. The remaining cells expressed a maternal Bpa1H or wt Nsdhl allele, depending on which X chromosome was inherited from the Bpa1H dam (Fig. 1).
|
|
Under standard culture conditions, the Bpa1H and wt cells were essentially indistinguishable in terms of morphology and growth rate. Because the Bpa1H cells would be expected to be fully dependent on the uptake of exogenous cholesterol from the culture medium, we compared the growth of parallel cultures of wt and Bpa1H cells in medium supplemented with LDS (Fig. 2E). After 24 h in LDS medium, the appearance of the two cultures was roughly equivalent. However, at 36 h, some of the mutant cells had rounded up and detached from the surface, and by 48 h in LDS medium, approximately half as many Bpa1H cells as wt cells could be seen attached to the surface, presumably as a result of cell death caused by cholesterol deficiency. The wt cells showed little to no cell death in LDS medium and subsequently grew to confluence at a normal rate.
Microarray analysis
RNA samples were isolated from three independent pairs of early-passage, sorted, primary cell cultures of Bpa1H and wt MEFs grown in both normal medium and LDS medium (see Materials and Methods). Because the Bpa1H MEFs appeared to grow normally in standard medium, we reasoned that the effects of the mutation on global gene expression might be more readily detectable under conditions of cholesterol deprivation. We chose three time points of growth in LDS medium (24, 36, and 48 h), because the observed cell death of the mutant MEFs over time would be expected to result in increasing secondary effects on gene expression associated with apoptosis or altered cell cycle regulation. In contrast, altered expression patterns in early time points or in all of the time points might reflect effects that were more directly related to Nsdhl loss of function.
Fluorescently labeled cDNA samples derived from mutant and wt RNA were cohybridized to oligonucleotide arrays (Agilent Mouse Development) representing
22,000 mouse genes from the NIA 15k and 7.4k cDNA clone sets (25). A total of 15 hybridizations representing 13 independent experiments and data sets comparing Bpa1H versus wt RNA were performed: five with cells grown in normal medium, four with cells in LDS medium for 24 h, two with cells in LDS medium for 36 h, and four with cells grown in LDS medium for 48 h (see Materials and Methods). To gain an overview of the microarray results, cluster analysis was performed on
1,525 genes that showed at least 1.5-fold higher or lower expression with P < 0.1 in the Bpa1H MEFs relative to wt cells in at least one of the time points (data not shown). One cluster, which included 25 genes, was characterized by upregulation of expression in the mutant cells in all three LDS time points relative to samples grown in normal medium (Fig. 3A). Not surprisingly, most of the genes in this cluster are known to be directly involved in cholesterol synthesis or in its regulation. In contrast, another cluster was characterized by genes that showed the greatest change in expression only in the samples from 48 h in LDS medium. It included several genes that are involved in cell cycle control, such as Cdc20, Ccne2, and CcnA2, probably reflecting deleterious secondary effects of cholesterol deficiency that accumulate over time (Fig. 3B).
|
, a multifunctional regulator of phosphoproteins (26, 27). The lower level of Nsdhl transcript in the Bpa1H-expressing cells may be a result of accelerated decay of the message attributable to premature termination of translation caused by the nonsense mutation (28). Based on searches of mouse expressed sequence tag databases (University of California Santa Cruz Genome Browser, http://genome.ucsc.edu/, mouse genome May 4, 2004 assembly), the oligo probe for Ywhaz (Agilent feature A_65_P09812) apparently represents sequence from the 3' untranslated region of a rare alternative transcript that extends the message
500 bp and uses a downstream polyadenylation signal. The functional significance of this, or several other variant Ywhaz transcripts with alternative 5' and 3' untranslated regions, is not known. A separate oligo probe from the coding region of Ywhaz (Agilent feature A_65_P13182) showed no significant differences between Bpa1H and wt samples, suggesting that the majority of Ywhaz message was unaffected in the mutant cells. Of those genes identified as upregulated in Bpa1H versus wt samples from the three time points of growth in LDS medium, a majority were involved in lipid metabolism. When the data for all three time points in LDS were combined and analyzed together, 44 genes had P values of <0.001 (Table 1). Four of these, Cyp51, Lss, StarD4, and Ctsd, were represented twice, because of the presence of two independent probes for the same gene on the microarray. Eleven of the 44 genes encode enzymes directly involved in cholesterol synthesis, 3 function as regulators of cholesterol synthesis, 7 are involved in fatty acid synthesis or metabolism, 2 are lipid transporters, 5 represent anonymous cDNAs, and 16 are genes with miscellaneous functions. The fold increase detected by the microarray analysis for most of the genes was generally 1.52.5. This relatively modest difference between the Bpa1H and wt cells was attributable in part to the fact that the wt cells also induced cholesterogenic genes when grown in LDS, as demonstrated by real-time PCR results (see below), thus partially masking the total increase in the expression of these genes in the mutant cells.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
To begin to investigate disease pathogenesis, we screened for changes in gene expression in fibroblasts cultured from Bpa1H embryos, taking advantage of an X-linked GFP transgene and random X inactivation to isolate pure populations of Bpa1H MEFs. The ratio of GFP+ to GFP MEFs in the initial unsorted cultures was not significantly different between wt and mutant embryos, indicating that the subpopulation expressing the mutant Bpa1H null allele proliferated at approximately the same rate as wt fibroblasts in the functionally mosaic female embryos, at least until E15.5. When cultured in normal medium, sorted populations of Bpa1H MEFs also appeared morphologically normal and grew at a rate comparable to wt cells. When cultured in LDS medium, however, the mutant MEFs showed substantial cell death by 48 h, presumably as a result of cholesterol deficiency, whereas wt cells continued to grow normally. In contrast, when unsorted populations of GFP+ and Bpa1H MEFs were cultured in LDS medium, cell death was not observed, and no change in the ratio of GFP+ to GFP cells was seen after 72 h (data not shown), suggesting that cells expressing the Bpa1H Nsdhl allele were rescued by cholesterol synthesized and exported by cells expressing the wt allele. The apparent ability of wt MEFs to provide sufficient cholesterol for the survival of Bpa1H cells would explain the lack of negative selection against Bpa1H fibroblasts seen in the E15.5 embryos. The use of the X-linked GFP transgene as an in vivo marker will allow further investigation into whether other cell types, including those affected in heterozygous embryos and adults, demonstrate a selective advantage for cells expressing the wt Nsdhl allele.
Finally, the death of Nsdhl null fibroblasts after 48 h of culture in LDS medium strongly suggests that the lanosterol derivatives that accumulate in these cells cannot support one or more vital functions of cellular cholesterol and/or that they themselves are toxic. In contrast, fibroblasts from several Smith-Lemli-Opitz patients homozygous for null alleles of DHCR7 survived for at least 1 week when cultured in LDS medium, suggesting that 7-dehydrocholesterol (7DHC) can substitute for at least some key functions of cholesterol, including maintaining the integrity of cell membranes (31).
Based on homology with yeast, it is predicted that the NSDHL protein would function in a complex with a sterol C-4 methyloxidase, encoded by the SC4MOL gene in human (32), and a 3ß-ketosterol reductase. During the course of our studies, Marijanovic et al. (33) demonstrated that the protein encoded by the human HSD17B7 locus functions as a cholesterol biosynthetic 3-ketosterol reductase. In Saccharomyces cerevisiae, the ERG28 gene encodes a regulatory protein for the complex that may act as a scaffold and anchor the C-4 demethylase enzymes to the ER membrane (34, 35). Consistent with this model, our results showed upregulation of Sc4mol, Hsd17b7, and Orf11, the murine ortholog of ERG28 (34), in Bpa1H cells grown in LDS. The increased expression of Orf11 provides the first functional evidence for a conserved role of the protein in mammalian cholesterol metabolism.
An initial question in understanding the pathogenesis associated with inherited defects in cholesterol biosynthesis is determining how the cell responds to both a deficiency of cholesterol and an abnormal accumulation of specific sterol intermediates at the blocked step in the pathway. Analysis of cultured human fibroblasts isolated from Smith-Lemli-Opitz patients that carry mutations in the DHCR7 gene showed that when 7DHC levels are significantly increased upon prolonged culture in delipidated medium, the activity of HMGR is suppressed, despite the presence of significantly lower total sterol levels in the mutant cells (31). Indeed, 7DHC was more effective in suppressing HMGR than cholesterol itself. Similarly, mice that are homozygous for a targeted null allele of Dhcr7, in which tissue cholesterol is decreased 5- to 6-fold and 7DHC is increased 30- to 40-fold, demonstrate decreased levels of HMGR protein and activity, probably as a result of accelerated proteolysis of the enzyme (36). Levels of mRNA for several cholesterol regulatory proteins and biosynthetic enzymes, including HMGR, remain unchanged between Dhcr7/ and wt animals. In contrast, cultured fibroblasts from patients with deficiency of mevalonate kinase, the enzyme subsequent to HMGR in the pathway, accumulate excess mevalonate and demonstrate a 6-fold increase in HMGR activity (20). Thus, mutant cells may show distinct responses to disruption of cholesterol biosynthesis depending on which step in the pathway is affected. Our analysis of Bpa1H MEFs showed that, unlike Dhcr7 mutants, loss of NSDHL function caused higher expression of 11 cholesterogenic enzymes, including HMGR, even after 24 h of growth in LDS medium. There was also a concomitant upregulation of genes for seven enzymes involved in fatty acid metabolism. Upregulation of fatty acid biosynthetic enzymes was reported in transgenic mice overexpressing sterol-regulatory element binding protein 2 (SREBP2) (37, 38), consistent with the lack of absolute specificity of SREBP1 and SREBP2 for fatty acid and cholesterol biosynthetic genes, respectively (38, 39).
In addition to the lipogenic enzymes, we detected modestly increased expression of the two regulatory genes, Srebp2 and Insig1. Srebp1 was not represented in the array, whereas expression of Insig2 was unchanged (data not shown). The upregulation of Srebp2 may account for the increased expression of the cholesterol biosynthetic enzymes, as well as the LDL receptor, because most of these genes have been shown to be direct targets of transcriptional activation by SREBP2 (3741). Although our study did not measure levels of nuclear SREBP2 in the MEFs, it has been demonstrated that SREBP2 activation and nuclear import is stimulated by cholesterol depletion (reviewed in 1, 2). Thus, we hypothesize that the higher level of Srebp2 mRNA in Bpa1H versus wt MEFs in LDS is responsible for or at least contributes to the increased expression of cholesterogenic enzymes observed in these cells. The observed increased expression of Insig1 was somewhat surprising given that the sterol-mediated binding of INSIGs to SREBP cleavage-activating protein (SCAP) results in the retention of inactive SREBPs in the ER (42). However, the expression of Insig1 itself is directly upregulated by SREBPs, whereas Insig2 is insensitive to regulation by SREBPs, except in the liver. Thus, the increased expression of Srebp2 observed in Bpa1H MEFs may account for the modest upregulation of Insig1 without concomitant increased expression of Insig2.
Among the genes showing upregulation in Bpa1H MEFs were two, Map17 and StarD4, that encode proteins whose precise function is unclear but have previously been implicated in cholesterol metabolism. MAP17 (for membrane-associated protein 17) plays a regulatory role in the uptake of HDL cholesterol through its interaction with PDZK1 in the liver (43). Among several proteins known to associate with PDZK1 at the plasma membrane is scavenger receptor class B type I, the receptor that mediates both the uptake of HDL cholesterol in the liver and cholesterol efflux from a variety of peripheral cell types (reviewed in 44). Overexpression of MAP17 in the livers of transgenic mice results in 2-fold increased plasma cholesterol levels and dramatic reductions in the scavenger receptor class B type I and PDZK1 proteins (43). MAP17 transcript levels are increased in livers of transgenic mice that overexpress either SREBP-1a or SREBP-2, but not in Scap knockout mice, suggesting that the MAP17 gene is either a direct or indirect target of SREBPs (38).
STARD4 is a member of the START domain family of proteins based on its homology to STAR (for steroidogenic acute regulatory protein), a protein that binds cholesterol and participates in its translocation from the outer to the inner mitochondrial membrane in steroidogenic cells (reviewed in 45). The gene was initially identified by microarray analysis as downregulated in livers of mice fed a high-cholesterol diet (41, 46). Although the exact function of the protein remains unknown, Stard4 expression in mouse 3T3-L1 cells decreases in the presence of increased levels of exogenous cholesterol and increases in the presence of lovastatin, suggesting a role in cholesterol homeostasis (46). Furthermore, the gene appears to be a direct target of SREBP regulation based on its increased expression in the livers of transgenic mice that overexpress SREBP and in Scap knockout mice (38).
Finally, one can envision how several of the other upregulated genes might be perturbed by a complete block in cholesterol biosynthesis. Cathepsin D is an endosomal aspartate protease that is activated by the sphingolipid ceramide and plays a role in ceramide-induced apoptosis as well as in epidermal differentiation (47, 48). Ceramide is generated in response to a variety of stress conditions and signaling molecules and influences cellular proliferation, differentiation, and apoptosis (reviewed in 49). Ceramide may also affect cell function through its role in membrane structure and microdomains at the Golgi, endosomes, and plasma membrane (reviewed in 50). The predicted amino acid sequence of one of the anonymous cDNAs identified as upregulated in our microarray analysis (NM_026708) showed 42% sequence similarity (E-value of 1E-18) to a domain in the yeast Lag1 gene that functions in ceramide synthesis (5153).
In light of the fully penetrant early lethality of the Bpa1H mutation in hemizygous males, our microarray analysis showed surprisingly few differences in gene expression between mutant and wt cells when the MEFs were cultured in normal serum. This result suggests that the abnormalities seen in vivo arise from cellular defects that are not associated with dramatic changes in mRNA level, that localized cholesterol deficiency may in fact occur in a subset of tissues of the mutant embryo but is not detectable by sterol analysis of the whole embryo, or that embryonic fibroblasts derived from midgestation embryos are less sensitive to the loss of NSDHL function than other cell types of the earlier, gastrulating embryo. Finally, the sorted Bpa1H MEFs used in this study will provide a unique and valuable tool for the in vitro analysis of possible perturbation of lipid raft structure, intracellular sterol trafficking, or signaling pathways.
| ACKNOWLEDGMENTS |
|---|
Manuscript received November 17, 2004 and in revised form February 1, 2005 and in re-revised form March 17, 2005.
| REFERENCES |
|---|
|
|
|---|