|
|
||||||||
1 From the Department of Ophthalmology and the 3 Columbia Genome Center, Columbia University, New York, New York.
| Abstract |
|---|
|
|
|---|
METHODS. The immortal mouse lens epithelial cell line
TN4-1 was conditioned to survive 125 µM H2O2 (H cells) or a combination of both 100 µM tertiary butyl hydroperoxide (TBHP) and 125 µM H2O2 (HT cells), by a methodology previously described. The total RNA was isolated from the different cell lines and analyzed with oligonucleotide mouse expression microarrays. Four microarrays were used for each cell line. Microarray results were confirmed by real-time RT-PCR.
RESULTS. A new cell line resistant to both 125 µM H2O2 and 100 µM TBHP was developed, because cells resistant to H2O2 were killed by TBHP. Analysis of classic antioxidative enzyme activities showed little change between cells that survive H2O2 (H) and those that survive H2O2 and TBHP (HT). Therefore, the global change in gene expression in these cell lines was determined with gene expression microarrays. The fluorescent signal changes of the genes within the three cell lines, H, HT, and control (C), were analyzed by statistical methods including Tukey analysis. It was found that from the 12,422 gene fragments and expressed sequence tags (ESTs) analyzedbased on a one-way ANOVA with a stringent cutoff of one false positive per 1000 genes and correcting for microarray background and noiseapproximately 950 (7.6%) genes had a significant change in expression in comparing the C, H, and HT groups. A small group of antioxidative defense genes were found in this population, including catalase, members of the glutathione (GSH)-S-transferase family, NAD(P)H menadione oxidoreductase 1, and the ferritin light chain. The remaining genes are involved in a broad spectrum of other biological systems. In the HT versus H comparison, only a few genes were found that had increased expression in the HT line compared with expression in the H line, including GSH-S-transferase alpha 3 and hephaestin. Many genes that are frequently considered antioxidative defense genes, including most of the GSH peroxidases, unexpectedly showed little change.
CONCLUSIONS. An unusual and generally unexpected small group of antioxidative defense genes appear to have increased expression in response to H2O2 stress. Cell lines resistant to H2O2 do not appear to survive challenge with another type of peroxide, TBHP, a lipid peroxide prototype. However, acquisition of TBHP resistance by H cells was found to be accompanied by significantly amplified expression of only a few additional antioxidative defense genes. Many of the amplified genes do not appear to be involved with antioxidative systems, reflecting the complexity of the cells response to oxidative stress.
| Introduction |
|---|
|
|
|---|
It has been shown that H2O2 is elevated in the aqueous, vitreous, and lenses in eyes of patients with maturity-onset cataract,19 20 21 and this increase correlates with extensive oxidation of lens components.11 22 23 24 Furthermore, H2O2 can cause opacification of the lens in organ culture at concentrations found in patients with cataract.17 25
To investigate the response to peroxide stress is complex. The lens contains only a single layer of epithelial cells, where most of the metabolic activity resides and a few layers of newly formed fiber cells that still have the biological machinery to make new protein directed by DNA. It is this thin layer of cells that is primarily responsible for the organs defense against oxidative stress. Even this situation is complex, however, in that the cells in the vicinity of the visual axis do not divide; but moving toward the equator, there is a zone in which cell division occurs, and finally, in the bow region, terminal differentiation is initiated, generating the formation of the fiber cells and the gradual disintegration of the protein synthesizing apparatus. Each cell type would be expected to respond in a unique manner to oxidative stress. Furthermore, the cellular response is not simple, involving shutting down some metabolic systems, increasing energy production, stimulating DNA repair systems, activating programmed apoptosis in some cells and stimulating the mobilization of antioxidative defenses. To sort out from this complex and time-dependent variation of gene expression and other mechanisms for modifying cell biology, those genes that specifically are capable of protecting the cell so that it can survive the stress is difficult at best.
Another approach is to investigate a homogeneous immortal cell type, such as the immortal murine lens epithelial cell line
TN4-1. Unlike primary lens cell cultures, there are relatively few changes from one generation to another. Gene expression is relatively constant. Therefore, this cell line has been conditioned to survive peroxide stress.
However, the relationship to the lens of the simian virus (SV)40transformed cells and their conditioned progeny should be considered. Even though the conditioned cells retain morphologic and cell growth features similar to those of unconditioned cells and their appearance is not significantly different from lens cell cultures, their cell biology is clearly different.26 This is also true of primary lens cell cultures, which lose many of their lens characteristics with passage. In this investigation, genes were found that were not expressed in unconditioned cell lines but were activated in the conditioned cell preparations. It would be expected that conditioned cells would have modified gene expression precisely to survive their new environment. The immortal cells also have been modified to take on characteristics required for eternal life. Thus, can information obtained from such modified cells be applied to the lens systems? The definitive answer cannot be forthcoming until the antioxidative defenses of the conditioned cell lines have been defined and the lens, enriched with such defenses, challenged with oxidative stress. If the lens remains transparent, then the exercise has been successful. However, it should be noted that work with other immortal cell systems such as MDCK,27 28 COS,29 30 CHO,31 32 HeLa,33 34 and 3T335 36 are widely used based on the assumption that biological systems are basically similar and that observations made in widely divergent cell types are applicable to each other. Furthermore, in these systems, this assumption has been confirmed repeatedly.
Analysis of the H2O2-conditioned cells indicated that not only were H2O2 detoxification systems amplified but repair systems appeared more effective.26 The unexpected finding that some GSH-S-transferases were markedly increased in activity and that catalase and not GSH peroxidase appeared to be the major H2O2-degrading enzyme indicated that our understanding of how the cell defends itself against H2O2 stress is incomplete. Therefore, to find the key genes involved with antioxidative defense, a global approach is required. The biology of the cell is controlled by gene expression. Thus, a reasonable rationale to assess modification in cell biology resulting from H2O2 conditioning is to determine changes in gene expression. Modification of gene expression in response to H2O2 stress has previously been shown for selective genes,37 38 39 but the entire genome has not been analyzed.
In the present work, murine expression high-density microarrays (Affymetrix, Santa Clara, CA), containing probe sets for approximately 12,400 gene fragments and expressed sequence tags (ESTs), were used to determine differential modification of gene expression in H2O2-conditioned cells (H line) and cells conditioned with both H2O2 and a lipid peroxide prototype, tertiary butyl hydroperoxide (TBHP), the HT line. We first determined that the array analysis is reliable. Then, using a rigorous statistical level of stringency (one false-positive selection in 1000 genes), we compared mRNA from four independent preparations of control cells (C cells) with four independent preparations from both the H and HT lines. An unusual group of antioxidative genes were found to have significant changes in expression. A few other interesting antioxidative defense genes with higher probabilities of error were also examined.
| Materials and Methods |
|---|
|
|
|---|
TN4-1, an immortal murine lens epithelial cell line, generously provided by Paul Russell, National Eye Institute (Bethesda, MD), was used in all work. The H2O2-conditioned cells (H) were prepared as previously described.26
The cells conditioned to withstand both an H2O2 and a TBHP stress (HT cells) were prepared in a similar manner. Briefly, approximately 200,000 cells were subcultured in a 35-mm dish in 2 mL minimum essential medium (MEM, cat. no. 41500-034; Gibco BRL, Grand Island, NY), supplemented with NaHCO3, 2.2 g/L medium (pH 7.2); 100 U penicillin and 100 µg streptomycin/mL medium (cat. no. 1540-122; Gibco BRL); Fungizone, 2.5 mg/L medium (cat. no. 15295-017); and fetal bovine serum (cat. no. SH30070.02; Hyclone, Logan, UT) to give a final concentration of 10%. This is the standard medium. After overnight incubation at 37°C in 5% CO2, the medium was replaced with 4 mL standard medium containing 125 µM H2O2 and various concentrations of TBHP. Every 24 hours, additional aliquots of H2O2 and TBHP were added to give the desired concentration. Culture medium was changed every 3 days. Cells were subcultured when they were approximately 80% confluent and after approximately 16 hours were again exposed to the peroxides. During a period of a few months, the TBHP concentration was gradually increased from a starting concentration of 25 µM to a final concentration of 100 µM, whereas the H2O2 concentration was maintained at 125 µM. H cells and C cells were maintained in an identical manner with 125 µM H2O2 or in standard medium, respectively. Cell viability was determined by cell counting and trypan blue staining, as previously described.26
Enzyme Assays
All enzyme assays were conducted by methodologies previously described.26
Isolation of Total RNA
Total RNA was isolated from 1 x 107 cells removed from 80% confluent cultures using an RNA isolation kit (cat. no. 74104; RNeasy Mini Kit, Qiagen, Valencia, CA).
Preparation of cRNA
Approximately 30 µg total RNA was used to prepare double-stranded cDNA using a commercial system (SuperScript Choice; Gibco BRL) and a T7-(dT)24 primer: 5'-GGCCAGTGAATTGTAATACGACTCACTATA-GGGAGGCGG-(dT)24-3' for first-strand cDNA synthesis. After synthesis of the second-strand cDNA, the double-stranded cDNA was then purified with Phase Lock Gels (Brinkmann Instruments, Inc., Boulder, CO), phenol-chloroform extraction and ethanol precipitation. An RNA transcript labeling kit (product code 42655; Enzo Diagnostic, Inc., Farmingdale, NY) was then used to prepare a biotin-labeled cRNA. The cRNA was purified using a spin column (RNeasy; Qiagen) followed by ethanol precipitation. The cRNA was quantified and examined by gel electrophoresis. The cRNA was then fragmented at 94°C for 35 minutes.
Gene Microarray Analysis
The use of oligonucleotide expression microarrays (Affymetrix) is now a frequently used methodology for determining gene expression. The methodology has been described in numerous publications.40
41
42
Briefly, the fragmented biotin-labeled cRNA was hybridized to a mouse expression microarray (MG-U74A Ver. 2; Affymetrix), at the Genome Center of the College of Physicians and Surgeons, Columbia University, and stained with streptavidin phycoerythrin. The hybridized components were detected by fluorescence laser scanning and confocal localization and then analyzed on computer (Microarray Suite, ver. 5.0, Affymetrix; and Excel; Microsoft, Redmond, WA). Each microarray contained two groups of approximately 16 twenty-fiveresidue oligonucleotides for each of the 12,422 probe sets and ESTs as well as internal standards. Each member of one group was a perfect complimentary match for a particular segment of a given gene fragment or EST. The second group contained a single base mismatch in the middle of each of the oligos. Thus, one probe set contains approximately 32 oligonucleotides for each gene fragment or EST, giving a total of over 400,000 oligonucleotides plus standards on the microarray. The oligonucleotides for each gene fragment are distributed on the microarray and collated by computer for further analysis.
It should be noted that in some cases, there is more than one probe set for the same gene. The probe sets recognize different sections of the gene, including up and down stream and splicing regions. In such cases, there may be considerable differences in hybridization and, thus, the recognized expression of the gene. In most cases, only probe sets showing the strongest fluorescent signal intensities were used. The actual fluorescent signal levels and present, absent, and marginal calls determined by the software program for each of the 12,422 probe sets in each of four control (C), four H2O2-conditioned (H) and four H2O2+TBHPconditioned (HT) microarrays are shown in electronic folder 1, available at http://www.iovs.org/cgi/content/full/43/10/3251/DC1. ANOVA values reflecting the probability of a significant change in expression in the peroxide-conditioned cell lines in comparison to the control group are given in electronic folder 2, available at http://www.iovs.org/cgi/content/full/43/10/3251/DC1. The genes are ranked on the basis of probability. Average signal ratios and Tukey evaluations for the first thousand genes are also presented. The probe sets are identified by Affymetrix numbers as well as by GenBank and Unigene accession numbers (GenBank and Unigene are provided in the public domain by the National Center for Biotechnology Information, Bethesda, MD, and are available at http://www.ncbi.nlm.nih.gov).
Analyses were performed on total RNA isolated from the C, H, or HT lines. For each line, four independently isolated RNA samples were used, each containing approximately 30 µg. Data from each microarray were scaled to give an overall equivalent fluorescence intensity and minimize variation. Average intensity was set at 250. For samesame comparisons, the array with the higher scaling factor was always used as the baseline. The gene array software generates a difference call based on the fluorescence intensities of the perfect match oligos versus their mismatch partners and also minimizes high and low fluorescence intensities, using the Tukey biweight analysis, thus, correcting for fluorescence intensities above the linear range or where the signal is initially below detection. If no artifacts were introduced in the preparation of the microarray, in the hybridization, or in the scaling, then there should be no difference in expression in samesame probe set comparisons.
Figure 1 shows typical scatter graphs in which the signal intensities of all genes on two arrays are plotted against each other. In Figure 1A , a comparison between two control arrays is shown. In theory, there should be no difference in signal intensity of a given probe set in the two arrays, and all points should fall on a line delineated by a ratio of 1. Indeed, the red dots delineating the probe sets present in both arrays were clustered around the no-difference line (a ratio of 1), the clustering becoming tighter as the signal intensity increases. The yellow spots denote probe sets with which the expression was marginal or have an absent call in both arrays. These probe sets were generally in the lower-signal-intensity range. Finally, the blue spots represent probe sets for which the detection call indicated the presence of the probe in one of the two arrays. The other lines indicate the ratio of the probe set intensities. The outermost lines, both above and below the ratio = 1 line, represent a ratio of four, the next lines 3 and the innermost lines 2. Most of the red dots fall within the twofold zone.
|
A major basis for error in measurement of global gene expression is technical. We addressed the technical problems by using four arrays of each cell type and using statistical analysis of the data to select valid results and validate conclusions with RT-PCR. The problem of noise has been examined in the following manner.43 Table 1 shows the results from four independent samesame comparisons. Two from the C groups and two from the H group. In a single comparison, a large number of false positives (either increase or decrease) were found. However, if two comparisons in each group are considered (i.e., the same gene must have a false call either increased or decreased in both sets to be counted), then the number of false-positive results becomes very small: twenty-one false positives in the C group and 22 in the H group. Furthermore, if all four comparisons are considered, there are no false-positive genes in common. Thus, the selection is random and disappears with additional comparisons.
|
Statistical Analysis
A more popular method of selecting and evaluating genes that have had changes in gene expression is to use ANOVA. In this method, the probability that a gene has had a significant change in expression can be determined. Therefore, one-way ANOVA was applied to the gene population with a stringent cutoff of one false positive per 1000 genes. All three populations (12 microarrays) were examined. Thus, a significant change between any two groups would be detected, and the analytical power is considerably increased. The average difference data for the fluorescent signal for each gene was analyzed with a standard linear one-way ANOVA. A Bonferroni correction was applied to the obtained 12,422 probabilities. Genes with corrected P < 0.001 were considered to have been significantly affected by the treatment. To differentiate between effects of H and HT compared with C, as well as to identify differences between H and HT, these genes were subjected to a Tukey test for the contrasts H versus C, HT versus C, and HT versus H. Statistical tests were run using the "R" statistical language (see http://www.r-project.org/). Electronic folder 2 contains the probabilities of the genes. Antioxidative defense genes were selected from the 950 genes with ANOVA P < 0.001.
Real-Time RT-PCR
The statistical analysis of the expression microarray data revealed a small population of antioxidative genes with significant increase in expression in the conditioned cell lines. Confirmation of some of these data were undertaken with real-time RT-PCR.44
45
46
This technique is currently the most accurate method for mRNA quantitation.46
Total RNA was used to synthesize first-strand cDNA. One microgram of total RNA from a given sample was reverse transcribed using oligo-p(dT)15 and avian myeloblastosis virus (AMV) reverse transcriptase (first-strand cDNA synthesis kit for RT-PCR; Roche Molecular Biochemicals, Indianapolis, IN).
Real-time quantitative PCR was then used to determine cDNA and the change in gene transcription. A fluorescein PCR detection system (LightCycler; Roche Molecular Biochemicals, Indianapolis, IN) was used for this purpose.
The reactions were performed in a volume of 20 µL of mixture containing 10 pmol of each oligonucleotide primer and 2 µL green fluorescent dye (LightCycler DNA Master SYBR Green I; Roche Molecular Biochemicals) containing Taq DNA polymerase, reaction buffer, dNTP mixture, and SYBR green I. To determine the amount of the mRNA in a given sample, the equivalent cDNA quantity is analyzed by real-time PCR. This analysis is based on the relationship of the number of PCR cycles necessary to obtain a fluorescent signal to the amount of cDNA. A typical set of data for a glyceraldehyde-3-phosphate dehydrogenase (GAPDH) DNA standard at a number of dilutions is shown in Figure 2A . The fluorescent signals are plotted against cycle number for samples with 10-fold serial dilutions. The solution with highest concentration of DNA is detected with the smallest number of cycles. A standard curve is then constructed by the software program by plotting the log of the concentration versus the cycle number at which the signal is detected (Fig. 2B) .
|
| Results and Discussion |
|---|
|
|
|---|
TN4-1 can be conditioned to survive an H2O2 stress capable of killing nonconditioned cells and causing cataract in vitro,26
these results raise certain fundamental questions, such as the ability of H2O2-conditioned cells (H), to withstand a lipid peroxide stress and how the cell biology has been altered to survive this hostile environment. As shown in Figure 3
, unconditioned control (C) cells exposed to 125 µM TBHP, a lipid peroxide prototype, died only slightly slower than cells exposed to 125 H2O2 when all other conditions were the same. Cells conditioned to survive H2O2, H cells, also die when exposed to 125 µM TBHP, although at a considerably slower rate than C cells. The presence of both TBHP and H2O2 does not significantly increase the rate at which the H cells die. To gain insight into these unexpected results, H cells were exposed to gradually increasing concentrations of TBHP while continuing to add 125 µM H2O2 to the medium. After approximately 6 to 7 months, the cells were able to survive both 125 µM H2O2 and 100 µM TBHP (HT cells). These cells grew almost as rapidly as the H cells under stressing conditions. The experiments indicate that cells that have been conditioned to survive H2O2 are vulnerable to other types of peroxide stress but can be further conditioned to resist both hydrogen peroxide and lipid peroxide challenge.
|
|
The results obtained in comparing the enzyme activities of HT, H, and C cell preparations suggest that the few enzymes examined do not reflect the overall cellular response to oxidative stress and in the case of the HT cells, give an incomplete indication of how the cell has adjusted its biology to survive the additional oxidative stress.
Analysis of mRNA from C and H Cells by Using Oligonucleotide Expression Arrays
To obtain information on the changes in these conditioned cell lines that may indicate how the cells survive particular oxidative stress, analysis of gene expression was initiated. For this purpose, oligonucleotide expression arrays (Affymetrix) were used. With all the cell lines, it was found that 43% to 51% of the genes and ESTs are expressed.
To detect antioxidative defense genes that have significant change in expression, a one-way ANOVA was performed on the total 12,422 genes, with the use of four C, H, and HT microarrays. The ANOVA probability for all genes is shown in electronic folder 2. The analysis is based on the actual fluorescent signals. For selection of the antioxidative defense genes, P
0.001 was usedthat is, the probability of selection of a false positive was equal or less than 1 in 1000. The genes were chosen from the list shown in electronic folder 2. In some cases in which more than one probe set for a given gene was present, it is included, even though one of the probe sets had a high ANOVA probability. Also, in some cases, a gene was selected because it either had a relatively high signal intensity in comparison to a reference group, or it was a member of a well-known antioxidative defense family. The selected 28 genes are shown in Table 3
and are ranked by probability. A brief description of the function of the gene is given as well.
|
A few members of the large GSH-S-transferase family are listed. These enzymes detoxify products of oxidation and xenobiotics, with GSH as a cofactor. Because the different GSH-S-transferase classes have different specificities and cellular locations, it is not surprising to find the expression of more than one affected by peroxide stress. Whereas GSH alone can detoxify the products of reactive oxygen species (ROS), its effectiveness is directly related to its concentration. It is well established that with H2O2 stress, GSH concentrations quickly plummet in lens epithelial cell preparations and in the epithelial cell layer of the lens itself.24 47 Thus, GSH cannot effectively protect the cell through nonenzymatic reactions for more than a very short period. However, by acting as a cofactor for the GSH-S-transferases, where the Km for GSH is approximately 0.1 mM,48 100-fold lower than the normal lens epithelial cell concentration,49 it continues to effectively degrade the toxic products generated by H2O2. Previous work has shown that GSH-S-transferase alpha 1 and alpha 2 have been markedly increased at the protein level in the H cell lines.26 Recently, it has been found that transfection of human lens epithelial cells with an alpha class GSH-S-transferase protects the cells against H2O2-induced lipid peroxidation.50 It is interesting that not all detected GSH-S-transferase transcripts had increased expression. As shown in Table 4 , the theta-1 enzyme mRNA was found to have decreased activity in both the H and HT line in relation to the C cells.
|
Heme oxygenase (HO)-1 is induced by oxidative stress, with induction increasing with diminishing concentrations of GSH.64 65 It degrades heme and produces antioxidant bile pigments. In the degradation of heme, iron is released. It has been suggested that HO participates in a coupled reaction with ferritin, sequestering and oxidizing the released iron and thus protecting the cell from iron catalyzed oxidation.58 The transcriptional regulation of the gene by oxidants has led to the viewpoint that HO provides a cellular defense mechanism to oxidative stress. HO has been shown to be induced in the retina with oxidative stress and is localized exclusively in the Müller cells.66 However, HO was suppressed in the H and HT lines.
Platelet-activating factor (PAF)-acetylhydrolases are members of a structurally distinct family of enzymes related to the phospholipases and specifically hydrolyze phospholipids, including oxidatively damaged components with a short chain at the sn-2 position. The hydrolases may detoxify oxidized phospholipid in conjunction with superoxide dismutase and GSH-Px-1.67 68 The ferredoxins are low-molecular-weight proteins with Fe-S clusters that are ubiquitous in both animals and plants. They are powerful reductants and are involved in a broad range of reactions, including carrying electrons from ferredoxin-reductase to membrane-bound cytochrome P450, which is important in the detoxification of xenobiotics.69 70
NADPH-dependent leukotriene B4 12-hydroxydehydrogenase is involved in the detoxification of
,ß unsaturated aldehydes and ketones. This is accomplished by the reduction of the
,ß carbon-carbon double bond, using NADPH as a cofactor. The enzyme is also known as alkenal/one oxidoreductase. It is effective in reducing cytotoxic lipid peroxidation products.71
It is striking that potent cancer chemopreventive agents upregulate a number of the genes noted on this list, including GSH-S-transferases, epoxide hydrolase, NADPH quinone reductase, and NADPH-dependent leukotriene B4 12-hydroxydehydrogenase.71
Aldehyde dehydrogenase II is a member of the large family of enzymes, including retinaldehyde dehydrogenase,72 73 that degrade xenobiotic aldehydes.74 75 These enzymes are involved in the metabolism of lipid peroxidation products, such as malondialdehyde and hydroxyalkenals as well as N-heterocyclic compounds, N-oxides, and nitrosamines.73 76 Although its overall physiological function remains obscure, it is clear that it can act as an antioxidative defense component. The aldehyde dehydrogenase II gene has been recruited as a crystallin (omega-crystallin) in cephalopods.77
Aldo-keto reductase 1 C13 is a member of the superfamily of monomeric proteins of approximately 320 residues that metabolize a broad range of aldehydes, including xenobiotics using NAD(P)H as a cofactor. Aldose reductase, an enzyme involved in the development of sugar-induced cataract is a member of this group.78 79 NAD(P)H menadione oxidoreductase, detoxifies quinones and quinone-imines and is capable of producing antioxidant forms of ubiquinone and vitamin E after oxidative stress.80 It utilizes both NADH and NADPH as cofactors and uses a two-electron-reaction mechanism.81 It has been suggested that the enzyme regulates the intracellular redox state by controlling the NAD(P)H-to-NAD(P) ratio and is induced by H2O2.82
Epoxide hydrolase 2 is the soluble form of the enzyme that is the major defense of the cell against xenobiotic-derived epoxides,83 84 85 as well as endogenous components such as epoxides of leukotrienes.86 The enzyme forms vicinal 1 to 2 diols that are generally less reactive than the epoxides.
Peroxiredoxins are an emerging family of antioxidant enzymes that regulate the concentrations of reactive oxygen species and protect tissues against oxidative attack.87 Peroxiredoxin 5 is a thioredoxin peroxidase found in mitochondria, peroxisomes, and the cytosol. This peroxiredoxin is unique in forming an intramolecular disulfide intermediate that is reduced by thioredoxin.88
Alcohol dehydrogenase (ADH) 3 appears to be identical with formaldehyde dehydrogenase.89 It acts on the product formed by formaldehyde and GSH, hydroxymethylglutathione. In the presence of nicotinamide adenine dinucleotide (NAD), it forms S-formylglutathione. It can also oxidize long-chain aliphatic alcohols and aldehydes but is not involved in retinaldehyde metabolism.90 It, thus, contributes to the removal of toxic aldehydes possibly formed by oxidative stress.
The copper chaperone for superoxide dismutase (SOD) interacts with the Cu/Zn enzyme (SOD1) delivering copper to the enzyme. Its presence is required for SOD1 activity.91
Glutaredoxin, also known as thiol transferase, is a member of the thiol-disulfide oxidoreductases and has been reported as a dehydroascorbate reductase.92 93 It appears to be important in maintaining ascorbate in its reduced form and in removing GSH and other thiols from thiolated proteins. It has been shown to protect human lens epithelial cells from H2O2 stress.94
Sulfide quinone reductase is an interesting protein that has been proposed to function in the detoxification and utilization of endogenous sulfides and appears to be present in eukaryotes. It is related to a mitochondrial enzyme HMT-2.95
GSH peroxidase (GPX) 4, the selenophospholipid hydroperoxide peroxidase (PHGPX), is an essential enzyme in protecting the cell against lethal peroxidative damage.96 97 It differs from GPX-1 in being a monomeric protein with very different substrate specificity. GPX-4 catalyzes the direct reduction of phospholipid hydroperoxides, whereas GPX-1 is unreactive with such hydroperoxides requiring sn-1 acyl bond hydrolysis before reducing the liberated fatty acid hydroperoxide.98 99
Secretory leukoprotease inhibitor (SLPI) is a low-molecular-weight serine protease inhibitor100
that also has been reported to elevate GSH levels in the lung101
and thus may contribute directly to both antioxidative defense as well as control of deleterious protease reactions. It inhibits activation of nuclear factor (NF)-
B.102
103
Lysosomal thiol reductase is induced by gamma interferon, IFN-
, and is frequently denoted as GILT. It is synthesized as a 224-amino-acid precursor that is transported to endocytic compartments. It has the ability to cleave disulfide bonds by a mechanism similar to that used by members of the thioredoxin family. It contains a CXXC motif reminiscent of the thioredoxin family of proteins. The enzyme requires a reducing agent for activation. Although L-cysteine acts in this capacity, there is no activity with glutathione.104
105
106
Aldehyde oxidase is a molybdenum dimeric enzyme that contains flavin adenine dinucleotide (FAD) and 2 Fe-S clusters and uses O2 as an electron acceptor.107 108 The enzyme has a major role in the oxidation of aldehydes, as well as the detoxification of a wide range of compounds. Recent work suggests that aldehyde oxidase is identical with retinal oxidase, converting retinaldehyde to retinoic acid.109
Macrosialin, a murine counterpart of human CD68, is a glycosylated transmembrane protein that scavenges low-density oxidized lipoprotein. Its overall role as an antioxidative agent is not clear at the present time.110 111
Additional information on these antioxidative defense genes is given in Table 4 . On the left side of the table, the relative fluorescence intensity rescaled from -2 to +2 for each probe set is shown. The weaker the fluorescent signal, the darker the color. Thus, the strongest signal is denoted by white. A standard color bar is presented above the gene data. With the color depiction of the signal intensities, it is possible to assess quickly the variation in the fluorescent signal intensity of a given probe set across all 12 microarrays. The average fold change for the H/C, HT/C and HT/H comparison is given. Its reliability has been evaluated by Tukey analysis in which the probability of a real difference has been set at 0.01 or less and is indicated by asterisks. Finally, the actual corrected signal intensities are given.
Some probe sets such as hephaestin, PAF-acetylhydrolase, and epoxide hydrolase 2 have a very high probability that their gene expression has been amplified, based on both ANOVA and Tukey analysis. Yet, the signal intensities are very low, raising questions concerning their significance.
Four of these genes, aldehyde dehydrogenase II, GSH-S-transferase-theta, heme oxygenase 1, and ceruloplasmin have decreased expression in the conditioned cell preparations. With the latter two genes, expression is down in both the H and HT lines. Only three genes (hephaestin, GSH-S-transferase alpha 3, and alcohol dehydrogenase 3) clearly had amplified expression in the HT but not in the H line. With all of these genes, except alcohol dehydrogenase 3, the signal intensities are low and with the GSH-S-transferase, the ANOVA probability is relatively high. There are also two probe sets, macrosialin and secretory leukocyte protease inhibitor with which the expression is increased significantly, based on Tukey, only with H2O2 conditioning and one, GSH-S-transferase theta 1, where it is decreased. A few genes have amplified expression in the H line and a further increase in the HT line. They include PAF-acetylhydrolase, menadione oxidoreductase I, alcohol dehydrogenase 3, and possibly aldehyde oxidase. Only the first two meet the Tukey test for significance. There are some genes that are essentially not expressed in the control cell preparation but are expressed with stress. Such genes include secretory leukocyte protease inhibitor, hephaestin, and microsomal GSH-S-transferase.
Real-Time RT-PCR
To obtain independent verification of the microarray data, real-time RT-PCR has been initiated with two independent internal standards: glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and mitochondrial ribosomal protein S24 (MRPS24). These housekeeping genes are frequently used as internal standards because of their stability. Similar results were obtained with both internal standards. Electronic folder 3 (available at http://www.iovs.org/cgi/content/full/43/10/3251/DC1) contains the forward and reverse primers used for the real-time RT-PCR analyses. The size of the PCR products are also given. They varied from approximately 250 to 383 base pairs, with one exception of 221 base pairs. Table 5
presents Tukey, microarray, and RT-PCR data for most of the selected antioxidative defense genes. The real-time RT-PCR results, in most cases, confirm the microarray data and the Tukey evaluation. However, the actual values can be significantly different. For example, with PAF-acetylhydrolase, the microarray gave an H/C value of 11, but the RT-PCR result was 5.6 (Table 5)
. Both values indicate a considerable amplification in gene expression. Indeed, the point of the exercise is to select genes that have their expression significantly modified. The RT-PCR results verify that this has occurred with PAF-acetylhydrolase and this conclusion is confirmed by the Tukey analysis.
|
|
Carper et al.,38 used differential display to examine the H and C cell lines. Increased expression was observed in a number of genes reported in their study. Of the genes they investigated, catalase had the greatest increase in expression, up 14-fold, reticulocalbin up 6-fold, and GSH-peroxidase and ferritin (both chains) up 2-fold. Metallothionein 2 was essentially unchanged. These results are in reasonable agreement with the present work, with the exception of the ferritin heavy chain, which appeared unchanged and GSH peroxidases, which may have a slight increase in gene expression.
Expression of Classic Antioxidative Defense Genes
The surprising absence of many of the common antioxidative defenses among the genes with amplified expression led us to survey some of the members of this group (Table 6)
. It was found that none of the GSH peroxidases except phospholipid GSH peroxidase was significantly amplified (ANOVA, P = 0.10.7) lending further support to the viewpoint that GSH-Pxs have limited roles as antioxidant defense enzymes and that the absence of GSH-Px-1 does not compromise animal viability.112
The thioredoxins and peroxiredoxins (except for peroxiredoxin 5) also showed little change. Nor was there a significant amplified expression of the phospholipase A2 family, which is believed to be involved in degradation of membrane lipid peroxides. Although the ferritin light chain appeared to have an increased expression, two different probe sets for the heavy chain showed no significant change in expression. On the basis of the present data, it appears that many of the well-defined antioxidative genes do not have amplified expression in the conditioned lens epithelial cell lines.
|
Based on this work, a relatively small number of antioxidative genes have their expression significantly amplified in response to peroxide stress and an even smaller number are found when already conditioned H cells are subjected to a different peroxide, TBHP. (It is possible that the presence of both H2O2 and TBHP obscures the cells response to just one stressing agent. Therefore, a cell line conditioned to resist only TBHP has been developed and is now being analyzed.) This suggests that there are a few key antioxidative components that determine whether a cell will survive a particular stress and that the same genes may not be involved or may not be sufficient to respond to other closely allied stresses. Some of the identified genes have not been previously considered in examination of the lens antioxidative defenses. Although, in most cases, the results obtained by microarray expression analysis have been verified by real-time RT-PCR, further work is needed at the protein level to confirm the conclusions drawn from this study.
The analysis presented in this article considers approximately 12,400 genes and ESTs. Although there is no certainty about the total number of genes in the mouse genome, it is probable, based on the human genome work, that there are only approximately 30,000 genes in the genome.114 115 116 On this basis, our work has examined approximately one third of the mouse genome. It will be interesting to find whether other important antioxidative systems will be revealed when the remainder of the genome is examined and annotation of genes with unknown function is completed.
| Acknowledgements |
|---|
| Footnotes |
|---|
Supported by grants from the National Eye Institute, Research to Prevent Blindness, Research to Cure Cataract Foundation, and the Department of Ophthalmology, Columbia University.
Submitted for publication August 27, 2001; revised April 30, 2002; accepted June 5, 2002.
Commercial relationships policy: N.
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Abraham Spector, Department of Ophthalmology, Columbia University, 630 West 168th Street, New York, NY 10032; as42{at}columbia.edu.
| References |
|---|
|
|
|---|