1 Medicine and
2 Pediatrics, Naomi Berrie Diabetes Center, Columbia University College of Physicians & Surgeons, New York, New York
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Two of these loci, obese (Lepob/Lepob) and diabetes (Leprdb/Leprdb), harbor mutations in the leptin (LEP) and leptin receptor (LEPR) genes, respectively (3,4,5,6). Leptin, a hormone secreted primarily by adipocytes in proportion to adipose mass, serves as a measure of total body energy stores. LEPR is widely expressed in at least five isoforms, but the LEPR-B isoform seems to mediate most, if not all, of leptins energy-regulating effects. Consistent with this hypothesis, the LEPR-B isoform is specifically truncated in the original Leprdb allele (4,5,6) and is most highly expressed within energy-regulating regions of the hypothalamus, including the arcuate, paraventricular, dorsomedial, ventromedial, and lateral nuclei (7,8). In the absence of leptin or a functional LEPR-B, rodents and humans develop a complex metabolic and behavioral syndrome of hyperphagia, hypogonadotropism, decreased energy expenditure, severe insulin resistance, and morbid obesity (9,10). Molecular characterization of leptin and LEPR have allowed identification of an important pathway in the mammalian energy homeostatic system.
Although rare humans and animal strains with mutations in single genes have provided key insights, a more complete understanding of mammalian energy homeostasis and its derangement in diabetes and obesity will require methods and analytical tools that can measure the actions and interactions of thousands of genes simultaneously. Microarrays can be used to monitor the variation and expression of tens of thousands of genes, measuring on a global scale differences between normal and pathologic states (11,12). Early studies using microarrays compared patterns of gene expression among various benign and neoplastic hematopoietic cells, normal and failing hearts (13), and, recently, adipose tissue from lean and genetically obese mice (14). These studies have been useful in identifying patterns of expression that are specific for pathologic states and, in the case of hematopoietic malignancies, identifying clinically important subclasses of lymphomas and leukemias (15,16).
Here we report the production of spotted microarrays and an assessment of their specificity in detecting differences in gene expression and their use to generate expression profiles from livers of Lepob/Lepob and wild-type mice. The liver serves as a primary center for gluconeogenesis, fatty acid metabolism, ketone body production, and lipoprotein trafficking and, therefore, is a key regulator of total body energy homeostasis. In Lepob/Lepob mice, each of these processes has been characterized and found to be perturbed, contributing to the overall leptin-deficient phenotype. These expression data identify sets of genes whose expression is altered in the absence of leptin and corrected by short-term caloric restriction or leptin treatment. These genes represent candidates for mediating specific aspects of the obese phenotype and the leptin-specific component of mammalian energy homeostatic pathways.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Arraying PCR products.
All arraying was performed with a four-pin-head GMS 417 Arrayer (Affymetrix, Santa Clara, CA) according to the manufacturers instructions at a density of 1,300 spots/cm2. The robotics systems are self-contained within a Plexiglas enclosure and maintained at a relative humidity of 4050% and temperature of 2326°C during arraying. The microarrays were stored in a low-humidity environment for >2 weeks before use.
Preprocessing of microarrays.
For ensuring covalent linkage of PCR products to the polylysine surface, microarrays were preprocessed on the day of hybridization. Each microarray was moistened over a bath of heated water (95°C) and exposed to 60 mJ ultraviolet energy. After cross-linking, the microarrays were washed serially in 0.1% SDS, 0.2x SSC, and 100% ethanol and dried by centrifugation in a tabletop centrifuge adapted for this purpose. Dried microarrays were prehybridized immediately.
Generation of fluorescently labeled cDNA from polyA RNA.
Total RNA from 50 mg liver was extracted with the use of an acid-phenol method according to the manufacturers instructions (Trizol; BRL Gibco, Gaithersburg, MD). From each sample, 0.51.0 µg polyA RNA was isolated using a polydeoxythymidine resin (Oligotex; Qiagen). PolyA RNA was transcribed into fluorescently labeled cDNA using an anchored oligo-dT primer (dT18MN), unlabeled deoxynucleotide triphosphates (dNTPs), reverse transcriptase, and deoxyuridine-triphosphate (dUTP) labeled with one of two fluorophores, Cy3 or Cy5 (17). Unincorporated dNTPs were removed, and the sample was concentrated with the use of a minifiltration system (YM-30; Microcon, Bedford, MA).
Hybridization and washing of microarrays.
Microarrays were prehybridized with 15 µl of a formamide-based hybridization solution (0.75 mol/l NaCl, 0.15 mol/l Tris-HCl [pH 7.5], 0.1 mol/l Na2HPO4, 0.1% Na4HP2O7, 2% SDS, 5x Denhardts solution, 100 µg/ml sheared salmon sperm DNA, and 50% formamide) for 1 h at 42°C in a humidified hybridization chamber. For each microarray, Cy5- and Cy3-labeled cDNA from two separate mouse liver RNA extractions were combined with 10 µl hybridization buffer, heated at 95°C for 10 min, and centrifuged for 5 min in a microfuge at maximum speed to remove any particulate matter. The combined sample was placed on a prehybridized microarray and covered with a coverslip. The coverslip was held in place by rubber cement, and the microarray was hybridized for 1214 h at 42°C in a humidified chamber. After hybridization, the coverslip was removed and the microarray was washed at room temperature in serial 3-min washes of 0.1x SSC/0.1% SDS and two washes of 0.1x SSC. The slides were dried by centrifugation (1,000g) in a tabletop centrifuge.
Fluorescent scanning of microarrays.
The microarrays were scanned according to the manufacturers instructions with the use of an Affymetrix 428 Array Scanner. The scan of each microarray generated two 16-bit TIFF image files, one representing the fluorescence emission from each fluorophore.
Extracting data from TIFF images.
Paired image files, generated from scans of microarrays, were analyzed with the use of Imagene software (BioDiscovery, Los Angeles, CA). This software package identifies spots within an array image and calculates mean, total, and median signal intensity values for the pixels that constitute each spot. In addition, the software calculates local background intensities and allows for outlier pixels to be excluded from calculations. In combination with CloneTracker 2.0 (BioDiscovery), the Imagene software correlates each spot with a specific cDNA clone and generates a spreadsheet output file containing mean signal and local background for each spot.
Calculating expression ratios.
After the raw data from a single microarray were extracted from two TIFF images, they were downloaded into a single spreadsheet. Subtracting the mean local background from the mean signal of each spot generated a net signal for each represented cDNA. When the net signal for any spot was both greater than the average background for the entire array and greater than the local background, that spot was used for subsequent calculations. When both of these criteria were not met, a spot and the cDNA that it represents were not likely to give reliable data and, therefore, were not considered in further calculations. The sum of the net signal for each spot on both images was calculated, and the ratio between the sums for each image was defined as the normalization factor. Individual normalized expression ratios were calculated by dividing the net signal of a given spot i from image A by the net signal from the same spot on image B, multiplied by the normalization factor:
Mean normalized ratio for cDNAi = (net signal A)i/(net signal B)i x ( net signal B)/(
net signal A)
Indentifying differentially regulated genes.
Each microarray experiment was performed in triplicate using RNA extracted from a total of six animals. Thus, in comparing wild-type (Lep+/Lep+ C57BL/6J) with obese (Lepob/Lepob C57BL/6J) expression profiles, RNA was extracted separately from the three mice of each genotype to generate labeled cDNA. The labeled cDNA from a single wild-type mouse was combined with cDNA from a single obese mouse and hybridized to a microarray. Normalized mean ratio data were calculated for each cDNA on each microarray, and those sets of genes that showed a consistent 30% or greater increase or decrease in expression between two genotype/treatment groups were identified. The normalized mean ratio data were then used to calculate mean ratios for genes that demonstrated consistent changes.
General animal care.
All mice were females purchased from Jackson Laboratory (Bar Harbor, ME) at 46 weeks of age and acclimatized for at least 1 week before experimental manipulation or killing. The mice were kept individually in Plexiglas, ventilated cages with cob bedding and except where noted had free access to water and a standard pellet diet, which derives 10% of calories from fat (PicoLab Diets). The temperature was kept at 22 ± 1°C, and a 12-h light-dark cycle was maintained. All animals were killed by CO2 asphyxiation 1 h into the light cycle. Columbias Institutional Animal Care and Use Committee guidelines were followed at all times.
Leptin injections.
Recombinant murine leptin was purchased from Peprotech (Princeton, NJ) and resuspended according to the manufacturers instructions in sterile phosphate-buffered saline, pH 7.4 (1 mg/ml), and sterile filtered (0.22-µm cellulose acetate filter). Resuspended leptin was stored at -20°C for no more than 2 weeks. Leptin-treated animals were weighed and received an intraperitoneal injection once daily 1 h before the dark cycle with 0.5 µg/g. Vehicle-treated animals were handled identically to the leptin-treated animals except that they received an injection of phosphate-buffered saline.
Northern blot analysis.
Ten micrograms of total RNA was extracted from the livers of Lepob/Lepob or wild-type C57BL/6J mice, electrophoresed through denaturing formaldehyde agarose gels, blotted to nylon membranes, and probed with 32P-labeled, randomly primed DNA probes using standard protocols (18). Probes were generated from purified PCR products used in arraying and subjected to cycle sequencing with T7 and T3 primers using recommended protocols and an ABI3700 automated sequencing apparatus.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
To assess the specificity of our arrays and protocols in the study of liver expression, we measured the false-positive rate as a function of threshold criteria and experimental repetition. We produced a set of arrays (M1.6K) that represented >1,600 murine genes. The cDNAs were taken randomly from a commercially obtained, sequence-verified, Unigene plasmid library. No effort was made to enrich the pool of cDNAs for liver-specific transcripts or for genes known or suspected to be involved in energy homeostasis. One microgram of polyA RNA was prepared from the liver of an 8-week-old C57BL/6J female mouse. The sample was divided equally; one half was used to generate Cy3-labeled cDNA, the other half was used to generate Cy5-labeled cDNA. The two differentially labeled samples were hybridized simultaneously to a M1.6K array, washed, and scanned with a fluorescent scanning confocal microscope. Using different lasers to excite each of the dyes, we obtained raw fluorescence data for each spot on the array. The pixels (10 µm x 10 µm) that define each spot were determined in a two-step process. Initially, a map that contained the location and size of each spot of DNA was loaded into an analysis program (Imagene, BioDiscovery) and subsequently confirmed visually by us. Thus, for each spot, total, median, mean, and background pixel fluorescent intensity were obtained for the two identical samples that had been labeled with two different fluorophores. From these intensity data, a ratio that represented the relative signal from each dye was calculated for every spot. Each ratio was multiplied by a normalization factor intended to account for systematic bias in the incorporation or scanning of the fluorescent dyes. For these calculations, any ratio <1 was transformed by taking its multiplicative inverse.
If the specificity of the expression arrays were ideal, then the normalized fluorescence ratio calculated for each spot in this experiment would equal 1. Any variation from unity, therefore, represents the technical specificity limits of our arrays and protocols. The distribution of ratios from this experiment was very tight (Fig. 1A). Fewer than 0.5% of spots had a normalized fluorescence ratio of >2. Thus, for a threshold of 2, a value used by many workers, the specificity of our array and protocol is >99%. As would be expected, the specificity of microarray expression analysis is a function of threshold ratioas the threshold is increased, the specificity increases with a concomitant decrease in sensitivity. However, independently repeating the experiment and combining data across array experiments also would be expected to increase the specificity but with less reduction in sensitivity. Repeating this experiment two more times, with polyA RNA prepared from a second and third liver, demonstrated a similar distribution of ratios for each individual array (data not shown). Averaging the expression ratios across all three arrays, while keeping the threshold ratio at 2, increases the specificity to 99.9% (Fig. 1B). If there is no gene-specific bias or systematic error, then another and more conservative method for increasing specificity is to set a minimum ratio threshold across a set of arrays. In our test arrays, for example, there were no spots for which the expression ratio was 1.2 or greater on each of the three microarrays (Fig. 1B). Thus, if we set a minimum expression ratio threshold of 1.3, then we will detect differences with a specificity of >99.9%. Or, stated more directly, a difference in expression of 30% or more across three microarrays will identify genes with true differences in expression, with a false-positive rate of <0.1%.
|
|
|
In human obesity, there is an increased risk of cardiovascular disease and mortality (21,22). However, Lepob/Lepob are relatively resistant to atherogenesis and have elevated concentrations of serum HDL particles (23,24). The basis of these two subphenotypes remains poorly understood. The elevation of HDL and its reduction by low doses of leptin result from altered hepatic metabolism of HDL (25). Among the genes altered in livers Lepob/Lepob mice, we identified four with potential roles in modulating serum cholesterol and thereby atherogenic susceptibility: HDL-binding protein, HMG-CoA synthase 2, HMG-CoA lyase, and apolipoprotein A-IV (Apo A-IV). All show increased expression in livers of obese mice. If the increased expression of any of these four genes mediates the atherogenic-resistant and elevated HDL subphenotype in obese mice, then their expression would be expected to normalize with short-term leptin treatment as HDL levels decrease.
Leptin-regulated gene expression in livers of obese mice.
Treating Lepob/Lepob mice with leptin rapidly reverses many aspects of the obese phenotype (26). Within 1 h of peripheral or intracerebroventricular leptin administration, serum glucose and insulin decline. Within 24 h, food intake decreases more than 50%, and within 3 days, total body weight is significantly lower. To identify genes whose expression is modulated by leptin, we compared the expression profiles of obese mice that were treated for 1 week with leptin with those treated with vehicle. After 1 week of daily intraperitoneal leptin treatment, obese mice lost 17% of their body mass, whereas the vehicle-treated obese mice continued to gain weight (4%). As expected, Lepob/Lepob mice that were pair-fed with leptin-treated mice lost less body mass (13%) than leptin-treated mice (17%). Liver RNA was extracted from three mice in each treatment group and converted individually to fluorescently labeled cDNA. cDNA from each leptin-treated mouse was paired with a vehicle-treated or pair-fed cDNA sample and hybridized to a microarray. Data from the six microarrays identified sets of genes whose expression was altered between obese mice and obese mice that were treated with leptin and between obese mice that were treated with leptin and pair-fed obese mice.
The profiles generated from arrays that compared leptin- and vehicle-treated mice identified 29 genes whose expression was consistently altered 30% or more (Table 2); the comparison of leptin-treated and pair-fed profiles consistently detected such differences in expression of 35 genes (Table 2). Of the 158 genes whose expression differed between livers of obese and wild-type type, 11 showed correction or partial correction with short-term leptin treatment. The number of genes with corrected expression is a small fraction of the total and likely reflects several inadequacies of short-term treatment in mice with long-term leptin deficiency. Specifically, long-term deficiency is likely to have developmental consequences that are not correctable with 1 week of leptin treatment. In addition, once-daily intraperitoneal injections do not mimic the complex responses of ambient leptin to feeding and circadian rhythms (27,28). Finally, the very stringent criteria that we set for identifying differentially expressed genes likely will overlook genes whose expression is only partially corrected by leptin.
|
|
Although hyperphagic and obese, leptin-deficient Lepob/Lepob mice are postulated to exist in a state of perceived starvation (29). Thus, with free access to food, they behave as though they have been calorie-restricted, eating voraciously. However, their hyperphagic response can be enhanced further with fasting (30). This suggests that there is a set of nutrient-sensing pathways that are activated further in the fasted or calorie-restricted obese mice (31). The set of genes (C) that increased in expression in obese and pair-fed obese mice, which are not correctable with short-term leptin treatment, may represent genes that are part of the augmentable hyperphagic response of Lepob/Lepob mice.
The phenotype of calorie-restricted Lepob/Lepob mice has not been explored fully, but in addition to an augmented hyperphagic response, obese mice survive longer than wild-type mice under severe calorie restriction (32). The subset of genes (G)only induced or repressed in calorie-restricted obese mice when compared with leptin-treated miceare candidates for genes in metabolic and signaling pathways that mediate the increased survival phenotype of obese mice. Finally, the direct comparison of expression profiles between mice that gained (obese vehicle) and lost weight (obese leptin) detected 18 genes with consistently altered expression. The expression of these genes (F) is independent of leptin and inversely regulated in negative and positive energy balance states.
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Leptin plays a central role in the regulation of feeding behavior, basal metabolism, glucose utilization, insulin sensitivity, and energy expenditure and thus overall in the regulation of energy homeostasis. In addition, leptin can modulate fertility, lipoprotein metabolism, bone turnover, and immune responses (34). Using microarrays, we identified more than 160 genes whose expression is altered in the livers of leptin-deficient Lepob/Lepob mice. By further comparing expression among leptin-treated and calorie-restricted Lepob/Lepob mice, we correlated the expression of subsets of genes with known phenotypic responses to pharmacologic and dietary interventions. The goal of these studies was to identify candidate genes that mediate specific metabolic subphenotypes of leptin deficiency.
Three genes specifically modulated by leptin (group A) were identified. In obese mice, these genes are upregulated and are corrected by short-term leptin treatment. In comparison with leptin treatment, however, the expression of these genes was not altered by calorie restriction. This result underscores the observation that leptins peripheral effects are not mediated solely by leptins ability to modulate feeding behavior. Indeed, most of leptins pleiotropic effects, including those that regulate insulin resistance, basal energy metabolism, fertility, lipoprotein metabolism, and immune responses, must be mediated, at least in part, by pathways that are independent of food intake and energy balance. Of the genes that we identified that are dependent on leptin but not on nutrient intake, two have not yet been characterized but do share sequence similarity to other genes. The third is an apolipoprotein.
Normally, Apo A-IV is secreted primarily into the portal circulation from the small intestine and rises in response to high-fat meals (35). A lesser amount is produced by the liver. Apo A-IV has been postulated to modulate feeding behavior negatively (36), regulate HDL metabolism (37), and impart atherosclerosis resistance (38,39). Our microarray profiles detected a 10-fold increase in hepatic expression of Apo A-IV in Lepob/Lepob mice. This increase is partially normalized after a short course of daily leptin injections. C57BL/6J mice that overexpress Apo A-IV in the liver have elevated serum HDL and are resistant to developing atherosclerosis, even when fed an atherogenic diet or when made genetically deficient in Apo E (38,39). Conversely, mice homozygous for a null mutation in the Apo A-IV gene have lower serum HDL (37). These data make Apo A-IV an excellent candidate gene for mediating the elevated HDL and atherosclerosis-resistant phenotypes of obese mice. The role of Apo A-IV in mediating these phenotypes can be tested directly. If, in obese mice, Apo A-IV is necessary for the elevated HDL levels and resistance to atherogenesis, then mice that are homozygous for null mutations in leptin and Apo A-IV, <Lepob/Lepob>;<ApoA4-/ApoA4->, will be more susceptible to atherosclerosis than the Lepob/Lepob mice. In addition, mice that carry the lethal yellow agouti allele are resistant to both exogenous leptin and atherogenesis (23), and we would predict that levels of liver expression of Apo A-IV will be elevated in <Ay/a+> animals.
Another of the leptin-regulated transcripts (Unigene cluster Mm.27136) is highly similar (94% nucleotide sequence identity) to a human gene, insulin-induced protein 2 (IIP2). IIP2 is expressed in regenerating livers (40) and is upregulated in hepatocytes treated with insulin (41). These studies show that in hyperinsulinemic obese mice, the IIP2-like transcript is upregulated approximately four-fold. Leptin treatment, which lowers the expression of this gene, also normalizes insulin sensitivity, resulting in decreases in both serum glucose and insulin concentrations. Calorie restriction of obese mice can improve hyperglycemia and modestly reduce insulin concentration, but such restriction does not reverse the underlying insulin resistance and does not normalize IIP2-like expression. Shimomura et al. (42) recently argued that in leptin-deficient states, insulin is unable to downregulate gluconeogenesis and glycogenolysis appropriately, whereas its ability to stimulate lipogenesis remains intact. Specifically, they suggest that within the liver of leptin-deficient mice, insulin signaling diverges along two paths, one that remains insulin-sensitive and stimulates excess fatty acid production and one that is resistant to insulin downregulation and leads to excessive glucose synthesis. If correct, the IIP2-like transcript is a candidate for mediating the insulin-sensitive responses in leptin-deficient mice. An initial test of this hypothesis would measure the transcriptional response of the IIP2-like gene in obese mice treated with insulin. If the IIP2-like gene is part of a pathway that remains insulin-sensitive in the leptin-deficient mice, then its expression should increase further in tandem with the transcription factor SREBP1-c, which remains insulin-responsive in obese mice.
In conclusion, the complex pathways that regulate short- and long-term energy homeostasis remain poorly defined. By defining transcriptional perturbations that occur in the livers of Lepob/Lepob mice, we identified pathways and candidate genes whose hepatic function is altered by leptin deficiency. Using microarrays, we demonstrated upregulation of genes required for lipid metabolism and identified candidate genes that may mediate the metabolic subphenotypes of resistance to atherosclerosis and insulin-responsive fatty acid overproduction. These results suggest that such studies could be expanded fruitfully to include a larger set of genes, to interrogate other energy-regulating tissues (hypothalami, adipose tissue, and skeletal muscle), and to include other murine models of obesity (Ay/a+, Cpefat/Cpefat, diet-induced). Through performance of more comprehensive expression analyses, additional candidate genes that contribute to obesity-related phenotypes and participate in regulating feeding behavior, energy partitioning, and energy expenditure could be identified.
![]() |
ACKNOWLEDGMENTS |
---|
![]() |
FOOTNOTES |
---|
Received for publication 2 April 2001 and accepted in revised form 19 June 2001.
Apo A-IV, apolipoprotein A-IV; dNTP, deoxynucleotide triphosphates; dUTP; deoxyuridine-triphosphate; HMG, hydroxymethylglutaryl; IIP2, insulin-induced protein 2; PCR, polymerase chain reaction; SSC, sodium chloridesodium citrate; SREBP-1, sterol regulatory element-binding protein 1.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|