|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 11, 8343-8355, March 16, 2007
High Copper Selectively Alters Lipid Metabolism and Cell Cycle Machinery in the Mouse Model of Wilson Disease*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Liver injury is the most common manifestation of WD, although neurological and psychiatric symptoms are also frequently observed (2, 4, 6). WD patients may show progressive hepatic cirrhosis, chronic active hepatitis, or rapidly developing liver failure (7). None of these clinical features is specific to WD, complicating its diagnosis. Remarkably, despite a long history (WD was described in 1912 (8)) and significant progress in characterization of its genetic basis (9-11), the molecular and metabolic changes that accompany initial stages of copper accumulation (and may serve as markers of disease progression) remain poorly characterized.
The effects of accumulated copper on lipid peroxidation, enzyme activity, and DNA stability have been reported in both patients and in animal models of WD (12-14), yet it is still unclear which of these manifestations are the primary effects of accumulated copper and which are the longer term consequences of the disease. Understanding the initial copper-specific changes is particularly important, not only because detecting such changes may serve as a useful diagnostic tool, but also because corrective or supportive treatment can be implemented to overcome observed metabolic alterations and further improve and supplement copper-chelation therapy.
To understand molecular events associated with the early stages of copper overload, we have utilized Atp7b-/- mice, an animal model of WD. We have shown previously that these animals accumulate copper to high levels and have several phenotypic features resembling WD (15). We have also observed that in the 6-week-old Atp7b-/- mice, the hepatic copper is at its highest level, yet the pathology is still minimal (15). At this presymptomatic stage the changes in gene/protein expression as well as associated metabolic alterations in the liver are likely to represent specific responses to copper. In this study, we characterize changes in the hepatic mRNAs, proteins, and metabolites in response to copper accumulation in 6-week-old Atp7b-/- mice. We demonstrate that the effect of copper is selective and involves distinct metabolic pathways, organized in the interconnected network, although oxidative stress is not apparent. The most significant and unexpected are the up-regulation of cell cycle machinery and the down-regulation of cholesterol metabolism; the latter is accompanied by marked decrease of cholesterol in the liver. These novel findings pave the way for detailed mechanistic analysis of WD and other copper-induced pathologies.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
RNA Isolation for MicroarrayImmediately after removal, the liver pieces (
75 mg) were frozen in liquid nitrogen and stored until further use. The total RNA was isolated using TRIzol reagent (Invitrogen) according to the manufacturer's protocol followed by an RNeasy cleanup procedure (Qiagen). The integrity of isolated RNA was electrophoretically verified by ethidium bromide staining and by absorbance ratio (A260 nm/280 nm > 1.8). A total of six liver RNA samples (three biological replicates for WT and three for KO animals) were tested for RNA quality. Each of the six sample targets was hybridized to two different MOE430A Affymetrix GeneChip arrays (technical replicates). The Affymetrix chip contained 22,691 spots corresponding to 13,016 unique Unigene IDs (i.e. unique genes).
Chip Performance AnalysisImage processing and analysis were performed using Affymetrix MAS 5.0 software. The resulting intensities and coordinate information were saved in a CEL file format and then subjected to global scaling with an average target intensity of 350 to allow for direct comparison of hybridization values from different targets. Scaled results for each sample were saved as CHP files, and these data were used to evaluate overall chip performance (supplemental Table 1). The analysis indicated that the parameters describing the quality of RNA, hybridization, and detection were all within acceptable range.
The data in CEL format were imported into GeneSifter.net, a microarray data analysis application and subjected to RMA normalization using published algorithms (17). The RMA-normalized data set was then used to identify changed genes and determine the statistical significance and magnitude of changes (fold change). The significance was established using Wilcoxon t test with Benjamini and Hochberg adjustment for multiple comparisons. The data were annotated with NCBI Unigene IDs and with related Gene Ontology terms by using the Affymetrics NetAffx tool and NIAID-NIH DAVID (apps1.niaid.nih.gov/david/). Of the 310 significantly differentially expressed probe sets, 304 had Unigene IDs and were annotated in mouse Gene Ontology (GO). In a number of cases, when changes in the transcript were detected by more than one probe set, the sequences of these probe sets were manually compared with the sequence of the identified transcript to ensure that the probe sets indeed mapped to the same gene.
The GO analysis was carried out using GenMapp on-line tool (18, 19). For each identified GO term, Gene Map z-statistic was calculated. A positive z-statistic for any particular GO term indicated that there were more genes associated with this GO term than would be expected by random chance. The GO terms with positive z-statistic were then ranked based on the percent of changed genes that belong to the corresponding GO compared with the total number of genes on a chip with the same GO. This analysis identified those GO terms that were most strongly represented among the differentially expressed genes (probe sets). The results were independently confirmed by applying the GO Browser tool in NetAffx to the same set of differentially expressed genes. The gene expression data have been submitted to GEO data base (GEO number GSE5348 [NCBI GEO] ).
The associations between altered genes or pathways were further evaluated using the Ingenuity Pathways Analysis software (Ingenuity® Systems). Affymetrix identifiers of the differentially expressed genes (the fold change of 1.5 or higher) and their corresponding expression values were loaded into the software and mapped to its corresponding gene object (so-called focus genes) in the Ingenuity Pathways Knowledge Base. The significance of the associations between the data set and the canonical pathway and functional annotations was calculated in two ways. First, the number of genes from the data set that map to the pathway was divided by the number of all known genes ascribed to the pathway. Second, the left-tailed Fisher exact test was used to calculate related p values and distinguish those functional/pathway annotations that had more focus genes than expected by chance. The networks of the focus genes were algorithmically generated based on their connectivity.
Real Time PCR (RT-PCR)Total RNA was extracted from livers of 4-, 6-, and 32-week-old WT and KO mice (n = 6 per group) and checked for integrity as described above. One-step RT-PCR was performed with a LightCycler instrument (Roche Applied Science) in a total volume of 20 µl containing 50 ng of total RNA, 0.5 µM each primer, LightCycler RT-PCR reaction mix SYBR Green I (1x), and LightCycler RT-PCR Enzyme Mix (Roche Applied Science). Reverse transcription was performed at 50 °C for 20 min. The denaturation and amplification conditions were 95 °C for 15 s followed by up to 35 cycles of PCR. Each cycle of PCR included denaturation at 95 °C for 15 s and then 10 s of primer annealing at 55 °C and 20 s of extension/synthesis at 72 °C (20 s). The temperature ramp was 20 °C/s, except when heating to 72 °C, when it was 2 °C/s. At the end of the extension, step fluorescence of each sample was measured to allow quantification of the RNA. After amplification a melting curve was obtained by heating at 20 °C/s to 95 °C, cooling at 20 °C/s to 60 °C, and slowly heating at 0.1 °C/s to 90 °C with fluorescence data collection at 0.1 °C intervals.
RT-PCR primers for the internal standard glyceraldehyde-3-phosphate dehydrogenase gene and the target genes were designed using OmigaTM software (Oxford Molecular, Oxford, UK). Gene sequences were taken from the Affymetrix web site based on Affymetrix ID and Unigene number. Primer design followed standard criteria; the complete list of primers is given in the supplemental Table 2. For selected genes QuantiTect® (Qiagen) primer assays were used (supplemental Table 2). Product identity was confirmed by separation on agarose gel and ethidium bromide staining. The abundance of target mRNA was calculated in relation to the glyceraldehyde-3-phosphate dehydrogenase mRNA in the same sample. The amount of target mRNA in Atp7b-/- mice relative to the wild-type mice was quantified using the 2
Ct method (20).
Human liver samples were obtained from eight Wilson disease patients who underwent liver transplantation and eight control livers from other patients who underwent liver resection. The study was approved by the institutional ethics committee of the University of Leipzig (registration number 236-2006) and followed ethical guidelines. RNA from small liver samples was isolated as described above. cDNA was synthesized from 2 µg of RNA by reverse transcription with Super-Script II RNase H- reverse transcriptase (Invitrogen) and random hexamer primers. Gene expression was determined by quantitative fluorogenic RT-PCR (ABI PRISM SDS 7900, Applied Biosystems). Primers and probes were selected to span two exons in order to prevent amplification of genomic DNA (supplemental Table 2). The PCR was prepared in a final volume of 12.5 µl of a reaction mixture containing 2.5 µl of cDNA (diluted 1:15), 5 mM MgCl2, 1.25 µl of 10x AmpliTaq buffer A, 200 µM dNTP (each), 0.3 unit of AmpliTaq Gold (Applied Biosystems), 200 nM 6-carboxyfluorescein-labeled oligonucleotide probe, and 900 nM of each oligonucleotide primer. The cycling conditions were as follows: 95 °C for 10 min and 40 two-step cycles of 95 °C for 15 s and 60 °C for 1 min. Standardization was performed using serial dilutions of linearized plasmid cDNA ranging from 10 to 107 copies. Data were analyzed with the ABI PRISM software. The mRNA expression results are given as the fold change of LDL receptor compared with control and for HMG-CoA reductase as the copy number normalized to 106 copies of
-actin.
Analysis of Soluble Proteins by Two-dimensional Gel ElectrophoresisThe liver tissue (50-100 mg wet weight) from control and Atp7b-/- mice was homogenized in 800 µl of buffer containing 10 mM HEPES, 10 mM NaCl, 1 mM KH2PO4, 5 mM NaHCO3, 5 mM EDTA, 1 mM CaCl2, and 0.5 mM MgCl2 using a glass Dounce homogenizer with tight pestle. Soluble proteins were obtained by centrifuging samples at 100,000 x g for 30 min and collecting the supernatant. The lipids at the top of the supernatant were removed using a 0.22-µm centrifugal filter device (ultrafree-MC, Millipore). The soluble proteins were precipitated by adding 4 volumes of ice-cold acetone, 1 mM HCl. The mixture was vortexed and incubated at -20 °C for 1 h, and proteins were pelleted by centrifugation at 15,000 x g for 10 min, and the pellet was air dried for 5 min. Pellets were then dissolved by vortexing in 8 M urea, and protein content was determined using a BCA assay (Pierce). 0.17 ml of this solution containing 2 mg of protein was mixed with 0.17 ml of 8 M urea, 4% CHAPS, 100 mM DTT, 4% glycerol, and 4% pH 3-10NL IPG buffer (Amersham Biosciences). This solution was used for overnight reswelling of 18-cm pH 3-10NL IPG strips (Amersham Biosciences). The strips were focused using a Protean IEF Cell (Bio-Rad) under the following conditions: 20 °C, 0-500 V for 6 h (rapid), 500-3500 V for 3 h (linear), 3500 V for 8 h, with a 50-mA limit per gel. Following isoelectric focusing, the strips were reduced and alkylated, and the second dimension separation was performed using 12% gels. Gels were stained with Coomassie G-250 as described previously (21).
Protein Identification by Mass SpectrometryThe spots of interest were excised from Coomassie-stained gels, washed twice with doubly deionized water, and then twice with 50% (v/v) 50 mM NH4CO3, 50% acetonitrile. The gel pieces were then incubated in 100% acetonitrile for 2 min after which the liquid was removed, and the pieces were air dried until white in appearance. Rehydration of the gel pieces was performed with digestion buffer containing 0.01 µg/µl trypsin (Princeton Applied Research), 50 mM NH4CO3, and 50 mM CaCl2 followed by the addition of 60 µl of the same buffer without trypsin and incubation at 37 °C overnight. The reaction was stopped by adding 3 µl of 98% formic acid (Aldrich) to
60 µl of digest. A 45-µl aliquot was analyzed by liquid chromatography-electrospray ionization-tandem mass spectroscopy using an ion trap mass spectrometer (LCQ Deca XP Plus, Thermo Electron Corp.). The analysis by liquid chromatography-electrospray ionization-tandem mass spectroscopy was performed using a capillary 180 µm x 12-cm column packed in-house with stable-bond C18 packing material (5 µm, ZORBAX, Agilent Technologies). The samples were applied to the column through a trap column (180 µm x 2.5 cm, packed in-house with the same packing material as above), and the peptides were separated using a linear gradient changing the solvent composition from 2% acetonitrile to 30% acetonitrile over a 30-min period with a constant flow rate of 1.5 µl/min. MS data were acquired in data-dependent mode in which a single survey scan (MS) was followed by up to four sequential data-dependent MS/MS scans on the four most intense peptide ions. A dynamic exclusion feature was used to extend the analysis to less abundant ions. The peptides were identified using SEQUEST (22), from within the Trans Proteomic Pipeline (version 2.71) (23). Scaffold (Proteome Software, Portland, OR) was used to analyze protein identifications derived from MS/MS sequencing results. Scaffold validates peptide identifications using PeptideProphet (24) and derives corresponding protein probabilities using Protein-Prophet (25). Proteins were considered as "identified" when probability scores were
0.99 (protein) and
0.90 (peptide) with a minimum of two peptides identified.
Serum AnalysisBlood was collected by cardiac puncture, and serum was separated by centrifugation after blood coagulation and used immediately for analysis. Serum lipids were quantified after isolation by sequential centrifugation as described previously (26).
Metabolic Profiling of Liver Tissue2 mg of frozen tissue was homogenized for 30 s in an MM301 ball mill (Retsch, Germany). Extraction was carried using 1 ml of a one phasic mixture of isopropyl alcohol:acetonitrile:water (3:3:2, v/v) at -20 °C for 5 min. After centrifugation, the supernatant was completely dried in a SpeedVac concentrator and then derivatized in two steps. First, carbonyl functions were protected by methoximation using 20 µl of a 40 mg/ml solution of methoxyamine hydrochloride in pyridine at 28 °C for 30 min. The samples were then derivatized using 180 µl of N-methyl-N-trimethylsilyltrifluoroacetamide (Macherey-Nagel, Dueren, Germany) at 37 °C for 30 min to increase the volatility. 0.5 µl of this solution was injected into an automatic liner exchange system unit (ALEX, Gerstel, Germany) at a 1:5 split ratio. For every sample, a fresh liner was used to avoid sample carryover. Each liner was deactivated and cleaned using a blank N-methyl-N-trimethylsilyltrifluoroacetamide injection followed by a flash heat ramp. The sample was introduced at 50 °C using an MPS injector (Gerstel, Germany) and heated to 250 °C using a 5 °C/min ramp.
An Agilent 6890 gas chromatography oven (Hewlett-Packard, Atlanta, GA) was coupled to a Pegasus IV time of flight-mass spectrometer from Leco (St. Joseph, MI). An rtx5-SilMS fused silica capillary column (Restek) of 30 m length, 0.25 mm inner diameter, and 0.25-µm film thickness was used for separation with the 0.25-mm inner diameter 10-m IntegraGuard column without film. For the analysis the gas chromatography oven was set to 85 °C with duration of 210 s and a following ramp of 20 °C/min. The target time was 330 °C with duration of 5 min. The transfer line temperature was set to 250 °C. Mass spectra were acquired with a scan range of 83-500 m/z and an acquisition rate of 20 spectra/s. The ionization mode was electron impact at 70 eV. The temperature for the ion source was set to 250 °C. Chromatogram acquisition, data handling, automated peak deconvolution, library search, and retention index calculation was done by the Leco ChromaTOF software (version 2.32) using an in-house custom mass spectral library for compound identification.
Analysis of Mature SREBP-2 in Nuclear ExtractsNuclear extracts were prepared essentially as described by Sheng et al. (27). Mouse livers (
0.7-1.0 g) were homogenized in 3 volumes of buffer 1 (10 mM HEPES, pH 7.6, 25 mM KCl, 0.15 mM spermine, 0.5 mM spermidine, 1.0 mM EDTA, 2.0 M sucrose, 10% glycerol, 1.0 mM DTT and Complete Protease Inhibitor Mixture (Roche Applied Science)) by six passes of a motor-driven Potter-Elvehjem homogenizer. Homogenate was layered over 1.0 ml of buffer 1 and centrifuged for 45 min at 30,000 rpm in a Beckman SW50 rotor. The supernatant was removed, and pelleted nuclei were resuspended in 10 mM HEPES, pH 7.6, 25 mM KCl, 0.15 mM spermine, 0.5 mM spermidine, 1.0 mM EDTA, 1.0 mM EGTA, 20% glycerol, 1.0 mM DTT, Complete Protease Inhibitor Mixture. The suspension was centrifuged for 5 min at 2000 x g at 4 °C. Pelleted nuclei were resuspended in 10 mM HEPES, pH 7.6, 100 mM KCl, 2.0 mM MgCl2, 1.0 mM EDTA, 10% glycerol, 1.0 mM DTT, Complete Protease Inhibitor Mixture. The one-tenth volume of 4.0 M ammonium sulfate, pH 7.9, was added, and the mixture was placed on a rocking shaker at 4 °C for 45 min. The samples were then centrifuged at 20,000 x g at 4 °C for 15 min. The supernatant contained nuclear extract. 10.0 µg of nuclear extract was then separated on a 7.5% Laemmli gel. Proteins were transferred to nitrocellulose in 10 mM CAPS, pH 11.0, 10% methanol. The membrane was blocked overnight in 50% Aquablock (EastCoastBio, North Berwick, ME). Anti-SREBP2 (Affinity Bioreagents, Golden, CO) was used at 1:5,000 dilution, and IRDye800-conjugated goat antirabbit (Rockland Immunochemicals, Gilbertsville, PA) was used at 1:20,000 dilution. Washes after both primary and secondary antibodies were in phosphate-buffered saline supplemented with 0.1% (w/v) SDS, 1% (v/v) Nonidet P-40, and 0.5% sodium deoxycholate as described in Sheng et al. (27). Immunodetection was visualized with an Odyssey infrared scanner (LI-COR, Lincoln, NE).
| RESULTS |
|---|
|
|
|---|
This result provided the first indication that there is a difference in gene expression between the control and Atp7b-/- livers even at the presymptomatic stage of disease. This conclusion was confirmed by quantitative evaluation of the data. The mean values, the range of signal intensities, and the patterns of mRNA signal distribution were compared for all WT and KO samples. These studies demonstrate that the mean and the range of signal amplitudes were very similar for all samples (Fig. 1A). However, the individual patterns of mRNA signal intensities fell into two distinct groups (WT and KO) corresponding to the sample genotype; the similarity between the members within each group was higher than between the groups (Fig. 1B). Altogether, the results in Fig. 1 and supplemental Fig. 1 demonstrate that the difference in the distribution of mRNA intensities between the WT and KO samples could be reliably attributed to the actual change in the individual gene expression rather than be caused by nonspecific variations between samples.
|
|
|
Verification of the Array Data Using Quantitative RT-PCRTo verify the results of oligonucleotide arrays, the RT-PCR analysis was performed for 21 transcripts encoding protein with known function (10 up-regulated and 11 down-regulated; for the list of RT-PCR primers see supplemental Table 2). The transcripts were selected to represent different metabolic pathways (proteins involved in cell cycle, metabolic enzymes, and metalloproteins) and different fold change. Altogether, we observed good correlation with respect to the direction of the change (up or down) (Fig. 2). For example, 9 of 10 tested transcripts, which were predicted to be up-regulated by the arrays, showed significant increase by RT-PCR, with the 10th transcript also showing an increase, which was less significant (Fig. 2A). Down-regulated transcripts also consistently showed a decrease in the RT-PCR analysis, although for 4 of 11 examined transcripts the change was small and/or not statistically significant (Fig. 2B). Weaker agreement between the array and RT-PCR data for the down-regulated transcripts could be due to the metabolic nature of their products and possible nutritional influences (see "Discussion"). In general, the magnitude of changes for all significantly altered transcripts (up or down) was larger in RT-PCR analysis compared with the oligonucleotide array, which is commonly observed in such type of experiments.
|
|
Correlation between the Extent of Changes in the Transcriptome and ProteomeThe relatively small number of transcripts that was changed in response to copper accumulation pointed to specific effects of copper on gene expression/mRNA stability and presumably on liver proteins. To independently evaluate the magnitude of copper-induced changes in the liver proteome, the patterns of soluble proteins from control and Atp7b-/- livers were compared using two-dimensional gel electrophoresis. The overall protein patterns for control and Atp7b-/- samples were very similar (Fig. 3), consistent with the limited effect of copper on liver transcriptome. Several changes, however, were detected. The most obvious change was a marked increase in the amount of small proteins with the apparent molecular mass of 15 kDa in the KO samples (Fig. 3). Mass spectrometry analysis identified these proteins as metallothioneins 1 and 2. This finding was in full agreement with the results of arrays and RT-PCR, both of which showed significant up-regulation of metallothioneins 1 and 2 (Table 1 and Fig. 2).
Copper-induced formation of reactive oxygen species and resultant oxidative stress have been discussed in the literature as major factors in the development of WD pathology. Therefore, one of the unexpected outcomes of our gene array experiments was the apparent lack of significant changes in protein machinery associated with oxidative stress. No change in mRNA level was detected for either catalase or copper zinc-dependent superoxide dismutase (SOD1), the key enzymes involved in peroxide and superoxide detoxification. Also, no up-regulation was detected for the heat shock proteins. A similar lack of changes was observed at the protein level, when SOD1 was examined. SOD1 is an abundant liver protein with known behavior on two-dimensional gels. It was easily detectable in our soluble samples (Fig. 3), and we confirmed the identity of SOD1 by mass spectrometry. No change in the amount of SOD1 in the Atp7b-/- liver was observed, supporting the results of the array experiments.
|
|
Consistently, the genes associated with mitotic cell cycle and lipid metabolism, particularly cholesterol biosynthesis, were among the most changed (Fig. 4). Steroids, glycerol, and alcohol metabolism as well as electron transport were also altered, whereas the transcripts encoding proteins involved in nucleic acid, amino acid, or carbohydrates metabolism were much less affected, if at all. (The complete list of the significantly changed GO terms is shown in supplemental Table 5.) The detailed hierarchical analysis of significantly changed GO terms also demonstrated that the up-regulated pathways were distinct from the down-regulated pathways and largely nonoverlapping.
By far, the largest fraction of up-regulated transcripts is associated with the cell cycle and more specifically with the regulation of mitosis, including chromosome segregation (Fig. 4; supplemental Table 6). These include SMC2 and SMC4 (structure chromosome maintenance), cyclins G1 and G2, and CDC20 and CDC28 protein kinases (for complete list see supplemental Table 6). The examination of hepatic nuclei from two pairs of WT and KO animals using electron microscopy showed that Atp7b-/- nuclei frequently contain electron-dense material that is less apparent in control nuclei (Fig. 5). This result pointed to structural and possibly functional changes in nuclear components of Atp7b-/- hepatocytes.
By contrast, the down-regulated transcripts were mostly associated with lipid metabolism and particularly cholesterol biosynthesis (Fig. 4, Table 2, and supplemental Table 6). These include both the biosynthetic enzymes and regulatory proteins. For example, mRNA for nine proteins involved in cholesterol biosynthesis were all significantly down-regulated; these include cholesterol 7-
-monooxygenase, farnesyl-diphosphate farnesyltransferase, HMG-CoA-reductase, lanosterol synthase, and others. The mRNA levels for main transcription factors controlling lipid and fatty acid synthesis and oxidation, such as sterol regulatory element-binding protein-1c and peroxisome proliferator-activated receptor-
, were decreased, although less strongly. Similarly, the mRNA for SPOT14 homologue (which can be involved in regulating the triglyceride pool and overall efficiency of de novo lipid synthesis) (29) was down-regulated.
|
|
|
We have reasoned that elevated copper may inhibit SREBP-2 maturation or re-localization of activated SREBP-2 into the nuclei, thus blocking cholesterol biosynthesis. However, Western blot analysis of SREPB-2 in nuclear extracts from control and Atp7b-/- livers indicate that mature SREBP-2 is produced and present at higher levels in the nuclei of Atp7b-/- mice (Fig. 7A).
The SREBP-2 Target Genes Are Misregulated in Atp7b-/- Mice and in Human Wilson Disease LiversTo evaluate SREBP-2 activity, we measured the amount of mRNA for low density lipoprotein receptor (LDLR). The transcription of the LDLR gene is directly controlled by SREBP-2 via sterol-binding element 1 and is expected to be up-regulated in response to increased nuclear SREBP-2 (30, 31). Fig. 7B illustrates that despite high levels of mature SREBP-2 in Atp7b-/- nuclei, the amount of LDLR mRNA was significantly lower in KO samples compared with control. These data suggested that elevated copper inhibits SREBP-2 function in Atp7b-/- livers.
|
HMG-CoA reductase, which plays a central role in cholesterol biosynthesis, is another important target of SREBP-2 that was down-regulated in Atp7b-/- mice (Table 2). Analysis of human HMG-CoA reductase mRNA by real time PCR demonstrated that HMG-CoA reductase transcript levels were markedly down-regulated (3.4-fold) in Wilson disease livers compared with control (Fig. 7C).
| DISCUSSION |
|---|
|
|
|---|
We conclude that the dysregulation of lipid metabolism and changes in the cell cycle machinery characterize presymptomatic stages of WD, whereas oxidative stress plays a minor role. Comparison of copper-induced mRNA changes with published mRNA profiles of liver responses to various stresses (32-35) illustrates that the overall pattern of changes in response to copper accumulation (i.e. the combination of most altered pathways) is distinct.
The Up-regulated PathwaysAlthough a large increase in the amount of metallothionein in response to copper accumulation was anticipated, the up-regulation of genes associated with cell cycle and chromosome structure was unexpected and very interesting. The up-regulated SMC2 and SMC4 mRNAs encode key components of the condensin complex, which is required for conversion of interphase chromatin into mitoticlike condensed chromosomes. Significantly, the level of up-regulation for SMC2 and SMC4 is similar, as would be expected for proteins that form heterodimers. The condensins association with chromosome and chromosome condensation is regulated via phosphorylation by CDC2 (36); the mRNA for mouse homologue of CDC2 is elevated in the Atp7b-/- liver (Table 1). Altogether, these transcriptional and presumably protein changes have a subtle but noticeable effect on nuclear structure. The Atp7b-/- nuclei show larger variation in size than control nuclei (our data and see Ref. 15) and frequently show electrondense material, the exact nature of which remains to be established (Fig. 5).
In addition to these structural changes, up-regulation of cyclins G1 and G2, CDC28 kinase, as well as elevated levels of the fetal liver antigen and securin are all suggestive of increased cell proliferation. This conclusion is consistent with our earlier observation of increased nuclear staining with the proliferative marker Ki-67 in the KO livers (15). Up-regulation of cyclin G1 was also observed previously in the presymptomatic LEC rats, another animal model of WD (37), suggesting that the copper-dependent alteration of a cell cycle machinery is a species-independent phenomenon.
Pathways Down-regulated in Atp7b-/- LiverThe most significant effect of accumulated copper is on genes involved in cholesterol metabolism, although fatty acid and bile acid biosynthesis are also affected. The extent of these changes is illustrated in Fig. 8, which demonstrates that the enzymes involved in key steps of cholesterol biosynthesis are all affected (down-regulated) by copper overload. Indeed, direct metabolite measurements revealed a marked decrease of liver cholesterol. The impaired cholesterol biosynthesis is also clearly reflected in the reduction of the VLDL fraction of serum cholesterol, whereas the observed decrease in triglyceride levels can be attributed to increased expression of lipoprotein lipase. It is interesting that recent clinical chemistry measurements in a group of 45 WD patients found that in these patients total cholesterol, LDL cholesterol, and
-tocopherol levels were significantly lower compared with the control group (38), indicating that copper-induced alterations of cholesterols homeostasis are shared by WD patients and Atp7b-/- mice. It would be important to establish whether or not the observed changes in the serum lipid levels can serve as additional biochemical markers facilitating detection of the early stages of WD.
Inhibition of SREBP-2 Function Is Likely to Contribute to the Atp7b-/- PhenotypeThe metabolic control of genes for the cholesterol, fatty acid, and lipid biosynthesis is well described. Low levels of metabolites up-regulate the rate-limiting biosynthetic enzymes through activation of specific transcription factors. For example, a decrease in intracellular cholesterol induces activation of SREBP-2, which in turn up-regulates genes involved in cholesterol metabolism. Our data indicate that this metabolic control is not disrupted by accumulated copper. Lower levels of cholesterol in Atp7b-/- livers are associated with increased levels of activated SREBP-2 in the Atp7b-/- nuclei, illustrating that neither processing of SREBP-2 nor its trafficking to nucleus are inhibited by copper overload. However, the target genes of SREBP-2 do not show expected up-regulation; in fact, they are markedly down-regulated. These results suggest that accumulated copper is likely to inhibit the SREBP-2 function either by directly binding and causing structural changes in SREBP-2 or by modulating SREBP-2 activity via other proteins. The precise mechanism for copper-dependent inhibition remains to be established; it may involve the known property of copper to compete with zinc and disrupt zinc fingers and/or copper-induced oxidation of Cys residues required for structure and function of transcription factors.
|
It is particularly significant that the effects on transcription/RNA stability that we observed in response to copper overload are opposite to the previously reported response to dietary copper deficiency. Specifically, in rats fed a copper-deficient diet, the fatty acid synthesis and assembly into triacylglycerols and phospholipids were up-regulated (42) in contrast to the down-regulation of fatty acid biosynthesis detected in our studies. These observations point to a tight connection between copper and lipid homeostasis, which is likely to be mediated through copper-dependent transcriptional control or modulation of RNA stability.
A Link between Lipid Metabolism and Cell CycleThe effect of copper on a limited number of metabolic pathways suggests that the corresponding pathways could be functionally linked. Indeed, analysis of relationships between the significantly changed genes using the Ingenuity Pathways Analysis software reveals a highly scored network of 36 up- and down-regulated transcripts (Fig. 9). In this network, the connection is observed between lipid metabolism via sterol regulatory element-binding factor 1 and cell cycle via cyclin-dependent kinase inhibitor 1. Search for this link in the literature uncovered recent work in which lipid deprivation was shown to activate sterol regulatory element-binding factor 1 and induce cyclin-dependent kinase inhibitor 1 mRNA and protein (43). In our studies, the decrease in liver lipids also coincided with the elevated amount of cyclin-dependent kinase inhibitor 1 mRNA.
Implications for Understanding the WD PathologyThe relationship between copper and specific metabolic pathways revealed in our study has important implications for understanding the pathology of WD. The toxic role of copper has been commonly associated with the ability of the metal to induce formation of reactive oxygen species. We demonstrate that at the presymptomatic stage of copper accumulation oxidative stress does not appear to play the major role, as no significant up-regulation of transcript and/or proteins involved in oxidative stress response is observed. This observation is reminiscent of findings by Armendariz et al. (44). These investigators performed gene expression profiling of fibroblasts chronically accumulating copper because of ATP7A inactivation, and they found no evidence of oxidative stress-related gene expression despite significant copper accumulation. Thus, copper-induced radical formation, if it occurs, is unlikely to be central to the early stages of WD pathogenesis.
At the same time, some level of oxidative stress cannot be excluded. The expression of Rad51, which is cell cycle-dependent (45), can also be increased as a result of DNA damage, suggesting that accumulated copper may have some effect on DNA structure. Similarly, the up-regulation of glutathione S-tranferase
was noticeable and was also reported in the LEC rat (37). This enzyme metabolizes bilirubin and exhibits glutathione peroxidase activity thereby protecting cells from reactive oxygen species and the products of peroxidation (46). The up-regulated selenoprotein H may also have an antioxidant function (28). Finally, up-regulation of transcripts for proteosome subunit (macropain), ubiquitin-conjugating enzyme E2C, and glyoxalase 1 suggested the need for removal/detoxification of metabolic by-products, which may play a larger pathogenic role at later stages of the disease.
Previously, in the Atp7b-/- animals nursed by homozygous knock-out dams, we observed microsteatosis at 6 weeks, which progressed to focal steatosis and lipid droplets in the nuclei of 12-20-week-old mice (15). In the 6-week-old Atp7b-/- mice nursed by heterozygous dams, the steatosis is less severe and varies from apparent to nondetectable (data not shown). It appears that dysregulation of lipid biosynthesis is accompanied by accumulation of certain (perhaps intermediate) lipid metabolites and that the degree of accumulation may depend on the nutritional status of the animal. If variation in nutrition does play a role in the severity of hepatic response to copper overload, this may explain the larger variation between the RT-PCR and gene array data for metabolic genes compared with the cell cycle genes, measured at 6 weeks when the homeostatic capacity of liver is still significant. A better understanding of the effect of nutrition on pathology development in Atp7b-/- liver may provide important insights into significant variation of WD phenotype in humans.
Proliferation of bile ductules is commonly observed at the intermediate stage of disease in WD patients (47). Similarly, the bile duct abnormalities are very apparent at the later stages of disease in the Atp7b-/- mice (15). Thus, it was interesting to see abnormalities of enzymes involved in bile acid biosynthesis. 7-
-Hydroxylase (CYP7A) involved in the classic neutral pathway was down-regulated, whereas microsomal oxysterol 7-
-hydroxylase (CYP7B1) was up-regulated. Bile acid synthesis defects may cause progressive cholestatic liver disease because of accumulation of potentially toxic intermediates (48). It is tempting to speculate that the dysregulation of bile acid biosynthesis may trigger pathological changes in the bile ducts observed in the Atp7b-/- livers.
In summary, our data provide insight into the initial stages of pathology induced by copper overload in livers lacking functional Atp7b. We demonstrate that the effect of copper is selective and involves specific metabolic pathways. This observation may have important practical implications, because specific cellular responses, particularly those involving metabolites, are more amenable to detection and possibly correction. We have also discovered that oxidative stress is not the major effect of copper at the early stages of the disease. Instead, lipid metabolism and cell cycle machinery are selectively affected by copper overload. The marked effect of copper on the specific metabolic pathways suggests that the extent of pathology may depend on the ability of tissue to correct/sustain these metabolic insults and thus be altered by "modifier" genes and nutritional status, which are different in different individuals.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Tables I-VI and Fig. 1. ![]()
1 Both authors contributed equally to this work. ![]()
2 Recipient of a postdoctoral fellowship from the Deutsche Forschungsgemeinschaft Grant HU 932/1-1. ![]()
3 To whom correspondence should be addressed: Oregon Health & Science University, L224, Sam Jackson Park Rd., Portland, OR 97239. Tel.: 503-494-6953; Fax: 503-494-8393; E-mail: lutsenko{at}ohsu.edu.
4 The abbreviations used are: WD, Wilson disease; VLDL, very low density lipoprotein; LDL, low density lipoprotein; LDLR, low density lipoprotein receptor; WT, wild type; KO, knock-out; DTT, dithiothreitol; CHAPS, 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonic acid; RT, real time; HMG, hydroxymethylglutaryl; MS/MS, tandem mass spectrometry; CAPS, 3-(cyclohexylamino)propanesulfonic acid; GO, Gene Ontology; SREBP-2, sterol regulatory-binding protein 2; SOD1, superoxide dismutase 1. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. S. Mangala, V. Zuzel, R. Schmandt, E. S. Leshane, J. B. Halder, G. N. Armaiz-Pena, W. A. Spannuth, T. Tanaka, M. M.K. Shahzad, Y. G. Lin, et al. Therapeutic Targeting of ATP7B in Ovarian Carcinoma Clin. Cancer Res., June 1, 2009; 15(11): 3770 - 3780. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Leonhardt, R. Gebhardt, J. Mossner, S. Lutsenko, and D. Huster Functional Interactions of Cu-ATPase ATP7B with Cisplatin and the Role of ATP7B in the Resistance of Cells to the Drug J. Biol. Chem., March 20, 2009; 284(12): 7793 - 7802. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. A. Muller and L. W Klomp Novel perspectives in mammalian copper metabolism through the use of genome-wide approaches Am. J. Clinical Nutrition, September 1, 2008; 88(3): 821S - 825S. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Gonzalez, A. Reyes-Jara, M. Suazo, W. J Jo, and C. Vulpe Expression of copper-related genes in response to copper load Am. J. Clinical Nutrition, September 1, 2008; 88(3): 830S - 834S. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. I. Vonk, C. Wijmenga, and B. van de Sluis Relevance of animal models for understanding mammalian copper homeostasis Am. J. Clinical Nutrition, September 1, 2008; 88(3): 840S - 845S. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A Roberts and B. Sarkar Liver as a key organ in the supply, storage, and excretion of copper Am. J. Clinical Nutrition, September 1, 2008; 88(3): 851S - 854S. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |