1 Center for Biotechnology and Genomic Medicine, Medical College of Georgia, Augusta, Georgia
2 Department of Pathology, Immunology and Laboratory Medicine, College of Medicine, University of Florida, Gainesville, Florida
3 North Florida/South Georgia Veterans Health System, Gainesville, Florida
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Diabetic nephropathy (DN) is the leading cause of end-stage renal disease in the U.S. and the largest contributor to the total cost of diabetes care. The incidence of end-stage renal disease has risen at an annual rate of 79% for the last decade (1), caused in part by an incomplete understanding of renal disease pathophysiology and limited therapeutic options to prevent disease progression. The development of spontaneous autoimmune diabetes in NOD mice provides a valuable model for human type 1 diabetes. Because these mice are prone to kidney complications, they are also suitable for DN studies (2). The most common diseases, such as diabetes and its complications, have a strong genetic basis as well as an environmental influence. Uncovering the genetic causes and their interaction with the environment will allow the development of novel strategies for the diagnosis, prognosis, prevention, and/or treatment of these devastating diseases. However, conventional approaches of investigation have limited the rate of progress due to the complex nature of these diseases. The recent genomic revolution has created unparalleled opportunities for biomedical research. In this study we assessed the impact of diabetes on the kidney transcriptome by comparing kidneys of a control group of nondiabetic NOD mice, a group of NOD mice at the onset of diabetes, and a group of NOD mice that had been diabetic for 1 month.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
RNA.
Immediately after the mice were killed, kidneys were flash-frozen in liquid nitrogen and stored at -80°C. Total RNA was extracted using a Qiagen mini-kit, and poly(A)+ RNA was purified using an Oligotex mRNA mini-kit. For histology and Northern blotting studies, the two poles from the left kidney were removed, snap-frozen, and stored at -80°C. The center one-third of the left kidney was fixed in 10% buffered formalin for routine histology and immunohistochemistry. The right kidney was first weighed and then prepared for histological evaluations.
Northern blot.
Total RNA was reverse-transcribed into cDNA for PCR amplification of the stearoyl-coenzyme A desaturase 1 (SCD1) probe. The primers for the SCD1 probe are a forward primer, 5'-CAT GCT CCA AGA GAT CTC CA-3', and a reverse primer, 5'-CAG AGC GCT GGT CAT GTA G-3'. After PCR, the probe was purified by electrophoresis, extracted from the gel, and 32P-labeled with the Primeit II Stratagene Random Labeling kit. Total RNA (15 µg) was loaded on a 1.6% agarose gel containing 20% formaldehyde and 1 x MOPS. Electrophoresis was allowed to proceed for 3 h at 150 V, after which the RNA was transferred overnight onto a BrightStar membrane (Ambion) in 20 x sodium chloridesodium citrate. After UV crosslinking and prehybridization, the filter was hybridized overnight at 42°C with the 32P-labeled SCD1 probe in ultrahybe buffer (Ambion).
Histology and immunohistochemistry.
Tissues were imbedded in paraffin, and 3-µm thick sections were stained with hematoxylin and eosin as well as periodic acid Schiff. Frozen tissue sections were fixed in acetone and washed with Tris buffer followed by peroxide (peroxiblock; Zymed, South San Francisco, CA) and serum and avidin blocking. Tissue sections were respectively exposed to different titers of SCD1 primary antibody (Santa Cruz Biotechnology, Santa Cruz, CA). Secondary antibody incubation and detection were performed with the ABC/DAB Elite kit (goat IgG type; Vectorlabs, Burlingame, CA) according to the manufacturer. Optimal SCD1 primary antibody titer was found to be a 1:40 dilution (2.5 µg/ml). Immunohistochemistry with SCD1 was conducted on paraffin sections. The tissue sections were deparaffinated in xylene and washed with 100% ethanol. Endogenous peroxidase activity was blocked with 3% H2O2 in methanol for 10 min followed by hydration in graded ethanol (95 and 70%) and distilled water. In preliminary tests, antigen retrieval in Trilogy (Cell Marque, Hot Springs, AK) at 95°C for 25 min, followed by standard washes, gave optimal results. Staining was comparable to the frozen sections but with better tissue structure. The tissues were then blocked with a serum/avidin solution from a Vectastain ABC/DAB Elite kit (goat IgG type; Vectorlabs). Excess serum was blotted, followed by a 1-h incubation with the SCD1 antibody at a dilution of 1:40 in diluent (Zymed) containing biotin block. The tissues were washed in Tris and incubated with biotinylated secondary anti-goat IgG antibody (Vectastain ABC/DAB Elite kit) according to the manufacturer. After washing with Tris, ABC reagent (Vectastain ABC/DAB Elite kit) followed by diaminobenzidine hydrotetrachloride chromagen (brown color) was subsequently added and washed with distilled water. The tissues were counterstained with hematoxylin, dehydrated, and mounted using Cytoseal-XYL medium (Richard-Allan Scientific).
cDNA library construction.
A subtracted cDNA library was created using the PCR-Select cDNA subtraction system (CloneTech, Palo Alto, CA). Briefly, 2 µg of "diabetic" poly(A)+ mRNA (pooled mRNAs from five diabetic mice at 1719 weeks) was used as the tester and 2 µg of "nondiabetic" poly(A)+ mRNA (pooled mRNA from five nondiabetic mice at 4 weeks) was used as the driver. The amplified fragments were then ligated into the pCR2.1 vector (Invitrogen, Carlsbad, CA), and a cDNA library was prepared in Escherichia coli INV F cells.
Microarrays.
The basic procedure for creating cDNA microarrays contains four major steps: amplification of cDNA clones, array printing, hybridization, and data acquisition and analysis. Details of the procedure can be found online at http://www.genomics.mcg.edu/microarray/.
Amplification of cDNA clones.
A total of 5,760 cDNA clone inserts were directly PCR-amplified from 2 µl bacterial culture. The optimized 100-µl PCR was performed using a primer pair complementary to the library vector on sequences adjacent to the insert cloning sites. Primers were as follows: TAF, 5'-CCG CCA GTG TGA TGG ATA TCT G-3'; and TAR, 5'-TCC ACT AGT AAC GGC CGC CAG-3'. The PCR products were then purified by ethanol precipitation, resuspended in 25 µl of sodium phosphate buffer/SDS (150 mmol/l; pH 8.5/0.01%), and stored at -20°C until needed for printing.
Printing slides.
The glass slides used in the production of the arrays were coated in-house with poly-L-lysine using a standard protocol. Purified PCR products were spotted onto the poly-L-lysineprecoated slides. The post-processing of the microarrays after printing follows the description by the Brown laboratory (available online at http://cmgm.stanford.edu/pbrown/protocols). The microarrays were then stored for up to several months in a dehumidified chamber.
Probe preparation and hybridization.
The reference RNA was a pool that combines equal quantities of total RNA from the kidneys of eight nondiabetic mice at 10 weeks. RNAs were converted to cDNA and labeled with fluorophores. The reference pool cDNA was labeled with Cy5, and each experimental sample of cDNA was labeled with Cy3. The labeling of cDNA was achieved by an amino-allyl coupling strategy (3). Briefly, an amino-allylmodified dUTP [5-(3-aminoallyl)-2'-deoxyuridine 5'-triphosphate (aadUTP)] was incorporated during cDNA synthesis from the RNA template (15 µg) in a reverse transcriptase reaction using Superscript II reverse transcriptase. After the synthesis, a monofunctional NHS-ester Cye dye (either Cy5 or Cy3; Amersham) was coupled to the modified dUTP in a sodium bicarbonate buffer. After coupling the Cye dye to the aadUTP, the unincorporated dyes were removed using a Qia-quick PCR purification kit.
Hybridization.
Each Cy3-labeled probe was combined to an equal quantity of the universal Cy5-labeled reference and hybridized to a microarray. The hybridization was accomplished under a cover slip in a hybridization chamber (Array-It) and allowed to anneal for 16 h at 60°C. Slides were then washed and spun-dried.
Data acquisition.
A GMS 418 Scanner (Affymetrix, Santa Clara, CA) was used to scan the images by scanning the slide twice, the first time at 532 nm and the second time at 635 nm. This process generates two 16-bit tagged image file format (TIFF) image files. Numerical values for each spot were extracted from the images using either Scanalyze by Michael Eisen (4) or Molecularware (Molecularware, Cambridge, MA).
A statistical program was used to identify/flag spots with low-intensity/background ratios (5). This flagging procedure allows us to 1) determine whether the data quality for each spot (gene) was sufficiently good to warrant subsequent analyses and 2) eliminate unreliable elements with expression statistically too close to the background. After "flagging," the data were uploaded to a flat file database, where the gene expression information was linked to the coordinates of the spot on the array.
Statistical and cluster analyses.
Because hybridizations were repeated three times for each RNA sample to improve the accuracy of the measurements, the median of the log (2) of expression ratios from replicates were calculated for each RNA sample using Microsoft Excel 97 SR2. Subsequent statistical analyses were based on the median expression data. We used Students t test (6) to select genes with the highest discrimination power between two groups of samples. Interesting genes were clustered using Cluster software (4) and viewed using TreeView (4).
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
|
Interestingly, when the new-onset and long-term diabetic groups were compared, only three genes reached the statistical threshold of P < 0.001 (Fig. 1B and Table 1). Even after we relaxed the statistical criteria, only three other genes had a P value of 0.0070.001 (Table 1).
Expression of the SCD1 gene and protein in NOD kidneys.
The SCD1 gene (Scd1) encodes a key enzyme in the control of membrane fluidity and lipid metabolism. This gene had the most discriminating power between the diabetic and nondiabetic mice (P < 10-7). Our microarray revealed that Scd1 was downregulated by 3.8-fold in the new-onset diabetic mice and by 4.3-fold in the long-term diabetic mice compared with the control NOD mice. Northern blotting on a separate set of animals that were not used in microarray analysis provided confirmatory evidence for the differential gene expression data from microarray (Fig. 2). The mean expression and SD for the control, new-onset, and long-term diabetic group is 10.2 ± 5.8, 4.1 ± 2.1, and 1.0 ± 0.8, respectively. This represents a 2.5- and 10-fold decrease in the new-onset and long-term diabetic groups compared with controls, consistent with the microarray data.
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
One interesting finding from our study is the relatively small number of genes differentially expressed between the new-onset and long-term diabetic mice, although both diabetic groups are very different from the age-matched controls. Only three genes reached the statistical threshold of P < 0.001. There is also a high probability that these genes are false-positive because of the small number of genes detected in this comparison. We believe that this is less likely to have been caused by our experimental approaches. Rather, it may reflect the paucity of physiological changes between new-onset diabetic mice and those 1 month after diabetes onset. The duration of diabetes in our study may not be long enough to cause significant changes. Therefore, it will be interesting to investigate the diabetic kidneys after longer duration of diabetes.
Consistent with the high metabolic role of the kidney, genes involved in oxidative phosphorylation and ATP synthesis are changed by the presence of diabetes. A number of genes encode enzymes that protect cells against oxidative stress, whereas others were implicated in transcription and translation, lipid and carbohydrate metabolism, and the cell cycle. In addition, we found genes responsible for protein trafficking, constitutive protein processing, and immune function. Finally, a number of channels and transporters were detected by microarray analysis. A number of the genes have been previously studied in the context of type 1 diabetes. We will only discuss a few of the important groups of genes that may relate to DN.
Oxidative phosphorylation.
NADH-dehydrogenase has decreased expression in the diabetic kidney. Interestingly, the activity of NADH-dehydrogenase has been reported to be significantly lower in the kidney of alloxan-induced diabetic rats (9). We found a reduced expression of cytochrome C oxidase subunit VIIb, a regulatory subunit of cytochrome C oxidase in the diabetic NOD mouse. In diabetic rats, the activity of cytochrome C oxidase has been shown to be significantly reduced in the kidneys (9).
Oxidative stress.
Hyperglycemia in diabetes generates free radicals by mechanisms that are thought to involve metal-catalyzed oxidation of glucose, oxidative degeneration, and protein glycation. These observations would imply that the kidney is affected by oxidative stress. It has been suggested that enzymes normally involved in detoxifying free radicals are partially incapacitated by nonenzymatic glycosylation in diabetic individuals. Our study demonstrates that glutathione peroxidase mRNA levels are increased in the diabetic kidney of NOD mice. Consistent with this result, a number of studies have found an increase in the activity of this enzyme. Specifically, alloxan-induced (10) and streptozotocin-induced (1114) diabetic rats have been shown to increase glutathione peroxidase activity in their kidneys when compared with controls. Other studies nevertheless suggest that the activity of glutathione peroxidase may vary during diabetes and depends on the specific location of the enzyme in the kidney. Jachec et al. (15) studied the renal cortex for glutathione peroxidase activity of streptozocin-induced diabetic rats in a time-dependent manner. The activity of glutathione peroxidase was significantly higher in diabetic rats at 5 weeks after onset of the disease. At 10 weeks of diabetes, a sharp decrease in activity was observed compared with 5 weeks, but the enzyme activity was still higher in the diabetic group. At 15 weeks, the decline in the glutathione peroxidase activity in the diabetic rats was such that it had become lower than in the controls.
In contrast to glutathione peroxidase, the expression of Cu/Zn superoxide dismutase 1 (SOD-1) was found to decline significantly in the kidneys of the diabetic NOD mouse. In agreement with this result, Genet et al. (10) noted a decrease in SOD-1 activity by 26% in the kidneys of alloxan-induced diabetic rats. Investigations with streptozotocin-induced diabetic rats, however, seem to indicate that the SOD-1 enzyme activity varies as a function of the kidney substructure and disease time course. A slight decrease in SOD-1 activity in the whole kidney was observed at 5 weeks postdiabetes onset (14), whereas others have demonstrated an increase in SOD-1 activity within 16 weeks after disease onset (14). Jachec et al. (15) followed SOD-1 activity over time in the diabetic rats. The observed SOD-1 activity levels in the renal cortex at 5 weeks postdiabetes onset were comparable between nondiabetic controls and diabetic rats. After 10 weeks postdiabetes onset, a slight increase in the activity of the enzyme occurred in the diabetic rats, whereas at 15 weeks postdisease onset, a decrease was observed in the diabetic rats.
Possible regulatory pathways in DN.
Changes in the transcriptome or proteome during the development of DN may include functional, pathogenetic, or compensatory functions. The alterations are best understood when taken in the context of the stages of the pathological process. The time points we have studied are generally considered early in the course of DN (1621). Nevertheless, striking hypertrophy occurs by 2 days (20,22,23) in the streptozocin model of diabetes. Rapid nephromegaly is also evident in diabetic NOD mice (2). This increase is principally caused by ornithine decarboxylasedependent hypertrophy (19,21,23) and a subsequent increase in the activity of diverse cellular components, from NaK ATPase to transforming growth factor-ß (TGF-ß). In this regard, our results indicate a >30% reduction in ornithine decarboxylase (ODC)-antizyme (Oaz1) message (Table 1). This is consistent with a complex mechanism for increased ODC activity via a double-negative regulatory mechanism involving transcriptome and proteome (24) with decreased transcription of the protein inhibitor of ODC.
In hyperglycemia, there is an increased filtered load of glucose that is reabsorbed by Na-coupled glucose transport (19,20). Secondary to Na-coupled glucose resorption, there is a decreased salt load in the distal nephron, and tubuloglomerular feedback (16) produces increased filtration (25). This response occurs within rapid physiological time (26). Although it is logical that the increased work of transport might contribute to hypertrophy, the counterintuitive events are that hypertrophy occurs before hyperfiltration. Therefore, other mechanisms, such as solute and volume shifts in the very early (<24 h) period (21) or a direct effect of glucose on the tubular epithelium (27), must initiate hypertrophy (22). Other diabetic complications of salt and volume regulation are exemplified by sodium and volume expansion exacerbated by sodium restriction (18) and hyperkalemic-hyperchloremic acidosis (17). It is this period of stable and persistent renal structural and functional abnormality that characterizes our study.
Many of the differentially expressed genes in this time period are regulated by TGF-ß and peroxisome proliferatoractivated receptor (PPAR) isoform signaling pathways (Table 1). TGF-ß isoforms are upregulated in this early phase of diabetes, so this is an unlikely mechanism for these genes, which are downregulated (2,28) by a direct feed forward or other positive regulation (29). In contrast, a variety of metabolic effects of diabetes are linked to decreased PPAR activities. Thiazolidinadiones (TZDs) and fibrates are agonists of PPARs and correct these effects (3032). The PPARs are potential candidates for differential downregulation of many genes in DN. We selected SCD1 for further study because it has pleotropic effects on cell membranes to lipid metabolism and may explain some of the effects seen in DN. In particular, a potential regulatory loop is suggested. The conversion of saturated long-chain fatty acids to unsaturated fatty acid is rate limiting and catalyzed by SCD1 (33). The unsaturated fatty acids are natural agonists of the PPARs (3032), which in turn regulate SCD1. Uncertainty in completing this potential regulatory cycle appears at this point because the TZD group of PPAR agonists downregulate SCD1 in cultured preadipocites, and arachidonic acid has a similar effect in mouse liver, but SCD1 is increased in mouse skeletal muscle (34); the clofibrate type of PPAR agonists upregulate SCD1 in mouse liver (35,36). Directed studies will be important in determining kidney-specific events. Nevertheless, if the composite effect of PPAR activity is positively trophic for SCD1 in the kidney, this could produce a feedback loop in DN.
SCD1 transcription is also positively regulated by TGF-ß in vitro and inhibited by overexpression of Smad7. Our study demonstrates a reduction of the FKBP class of transcripts (Table 1), which encode protein that bind and inhibit TGF-ß. This is a potential mechanism for increased TGF-ß signaling. Again, renal-specific effects need to be established to define the interactions of these various trophic influences. Because there is variable expression of SCD1 throughout the nephron (Fig. 3), cell-specific interactions will need to be addressed.
We discussed only a few of the potential regulatory pathways suggested by the results of this study. Nevertheless, the approach that we have taken to dissecting the kidney transcriptome in diabetes is powerful in that it reduces investigator bias in data collection and permits the concurrent analysis of a large number of genes. Certainly, it has raised more questions and provides a valuable paradigm for moving forward with this complex disease.
![]() |
FOOTNOTES |
---|
Address correspondence and reprint requests to Dr. Jin-Xiong She, Center for Biotechnology and Genomic Medicine, Medical College of Georgia, 1120 15th St., PV6B108, Augusta, GA 30912. E-mail: jshe{at}mail.mcg.edu
Received for publication March 7, 2003 and accepted in revised form April 22, 2003
aadUTP, 5-(3-aminoallyl)-2'-deoxyuridine 5'-triphosphate; DN, diabetic nephropathy; ODC, ornithine decarboxylase; PPAR, peroxisome proliferatoractivated receptor; SCD1, stearoyl-coenzyme A desaturase 1; SOD-1, superoxide dismutase 1; TGF-ß, transforming growth factor-ß; TZD, thiazolidinadione
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|