|
|
||||||||
1 Department of Biochemistry, University of Cambridge, Cambridge
2 Mammalian Genetics Unit, Medical Research Council (MRC) Harwell, Oxfordshire
3 The Mary Lyon Centre, MRC Harwell, Harwell, Oxfordshire
4 Safety Assessment, GlaxoSmithKline, Ware, Herts, United Kingdom
| ABSTRACT |
|---|
|
|
|---|
metabonomics; metabolic syndrome; biofluids; Zucker rat; db/db mouse; nuclear magnetic resonance spectroscopy; leptin resistance
| INTRODUCTION |
|---|
|
|
|---|
To improve the understanding of the early stages of T2DM development, a number of animal models have been developed. We have considered two such models: the C57BL/KsJ db/db mouse and the obese Zucker (fa/fa) rat. The db/db mouse has a single-gene autosomal recessive defect in the leptin receptor gene. Leptin is a hormone produced primarily by white fat cells and is involved in regulating body weight and energy homoeostasis (31). In the db/db mouse, the concentration of insulin steadily increases with age until
810 wk, after which the concentration declines to below control nondiabetic levels (35). The db/db mouse produces clinical signs of leptin resistance, hyperphagia, obesity, and subsequent insulin resistance. The relevance of this model to human T2DM has been debated, as mutations in this gene in humans are rare (40), but the subsequent T2DM is characteristic of the metabolic changes in humans with severe disease. The obese Zucker (fa/fa) rat (9, 47) also has a single-gene autosomal recessive mutation in the leptin receptor (43), causing clinical signs of leptin resistance, obesity, hyperlipidemia, hyperinsulinemia, and, post-6 wk of age, fasting hyperglycemia and T2DM (32). The heterozygous (fa/+) lean genotype rats remain normoglycemic. Because the relevance of the db/db mouse and Zucker rat models of T2DM has been questioned, this study aims to determine any metabolic similarities or differences among the two rodent models and human subjects with T2DM, and in particular to identify any important perturbations in metabolism common to all three species.
Here we describe the application of 1H-NMR spectroscopy-based metabolomics, combined with multivariate and univariate statistics, to investigate the urinary metabolic profiles in two animal models of T2DM. We have compared these metabolic changes with perturbations observed in a human population of unmedicated diabetic patients who have good daily dietary control over their blood glucose concentrations by following the guidelines on diet issued by the American Diabetes Association. Diabetes can lead to pathological concentrations of several metabolites in plasma that are ultimately detected in urine (27, 33), an effect exacerbated by glycosuria. High blood glucose concentrations cause impaired solute reabsorption from the tubular lumen and reduced efficiency of the epithelium (7, 27). Multivariate analysis has been applied to delineate the effects of hyperglycemia and general metabolic stress responses from other potentially more discriminatory metabolic pathways.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Zucker rat urine collection.
Obese Zucker (fa/fa) rats and lean (fa/+) rats were maintained in stable colonies at GlaxoSmithKline (Ware, UK; temperature 21 ± 2°C, relative humidity 55 ± 10%, and fluorescent lighting 06.0018.00 GMT). Rats were fed on standard laboratory chow (2014 Teklad Global 14% Protein Rodent Maintenance Diet, Harlan) ad libitum. Urine samples were collected from male rats (n = 8 for lean and fatty Zucker rats) housed in individual urine collection cages after acclimatization and collected over solid CO2 at two intervals (08 and 824 h) at 1, 4, 8, and 12 wk. Urine samples (128) were collected and stored at 80°C. Body mass and food and water consumption were monitored daily. All fa/fa rats exhibited a marked diabetic phenotype confirmed by blood biochemistry analysis (data not shown).
Human urine collection.
For the clinical study, participants gave informed consent that met all criteria required for entry into the study (see supplementary information for details). The study was conducted in accordance with "good clinical practice" and all applicable regulatory requirements, including the 1996 version of the Declaration of Helsinki. GlaxoSmithKline provided the investigators with any relevant document(s)/data that were needed for IEC/IRB review and approval of the study. For the human studies, midstream urine (
15 ml) samples were collected and frozen from each volunteer. In total, 84 samples were collected from 12 healthy volunteers (7 time points, 8 males and 4 females) and 50 samples from 30 T2DM patients (13 time points, 17 males and 13 females) with well-controlled blood glucose maintained at normal concentrations by diet, following the guidelines issued by the American Diabetes Association, rather than medication. The healthy subjects were aged 1855 yr, had a body mass index (BMI)
19 and
30 kg/m2 and a body mass
50 kg and
113 kg, and were free from any major disease or pregnancy. The T2DM patients were aged 3065 yr (mean 56 ± 9 yr), had a BMI >25 and <40 kg/m2, weighed between 65 and 140 kg (mean 95 ± 19 kg), and were taking at most one oral anti-diabetic drug. T2DM patients agreed to stop treatment with oral anti-diabetic agents during the study. Subjects went through a washout period of 4 wk before sample collection and abstained from alcohol during the study; diet was controlled throughout the study. Clinical biochemical measurements for the T2DM patients during the study were taken, with most analytes within normal range; cholesterol, glycosylated hemoglobin, high-density lipoprotein (HDL), and triglycerides were slightly elevated (Supplemental Table 3). Stringent selection criteria were also placed on blood pressure and renal function to exclude diabetic patients with complications associated with extreme obesity, high blood pressure, and renal dysfunction (Supplemental Materials). Two outlying diabetic patients were identified as a result of over-the-counter medicine (paracetamol) in the samples. These data were excluded from further analysis.
NMR acquisition and processing methods.
When placed in a strong magnetic field, each chemically distinct proton in a solution will exhibit a unique chemical shift that depends on the exact stereochemical environment surrounding that proton. The intensity of this signal depends on the concentration of the proton, and hence of the metabolite, in the solution. By assigning each chemical shift to a metabolite and analyzing the relative changes in signal intensity between the disease and control samples, we could monitor the changes in metabolite concentration for a wide range of metabolites simultaneously.
Aliquots of either 200-µl (mice) or 400-µl (rat and human) urine samples were made up to 600 µl with phosphate buffer (0.2 M, pH 7.4) and any precipitate removed by centrifugation. In total, 500 µl of supernatant were transferred to 5-mm NMR tubes with 100 µl of sodium 3-trimethylsilyl-(2,2,3,3-2H4)-1-propionate (TSP)/D2O/sodium azide solution (0.05% wt/vol TSP in D2O and 1% wt/vol sodium azide).
NMR spectra of the db/db mice urine samples were recorded on a Varian INOVA spectrometer (Varian, Palo Alto, CA) at a proton frequency of 400.1 MHz, rat urine spectra on a Bruker DRX600 spectrometer (Bruker BioSpin, Rheinstetten, Germany) at a proton frequency of 600.1 MHz, and human urine spectra on a Bruker DRX700 NMR spectrometer at a proton frequency of 700.1 MHz (see Supplemental Materials). The one-dimensional (1D) NOESY pulse sequence with water presaturation was used throughout. NMR spectra were assigned with reference to the literature (15, 30) or confirmed by 2D spectroscopy including homonuclear 1H-1H correlation spectroscopy (COSY), 1H-13C heteronuclear signal quantum correlation (HSQC), and 1H-13C heteronuclear multiple bond correlation (HMBC) spectroscopy.
Spectra were processed using ACD/1D NMR Manager 8.0 with Intelligent Bucketing Integration (Advanced Chemistry Development, Toronto, ON, Canada). Spectra were integrated 0.209.30 ppm excluding water (4.245.04 ppm), glucose (3.193.99 ppm, 5.215.27 ppm), and urea (5.046.00 ppm). Intelligent bucketing ensures that bucket edges do not coincide with peak maxima, preventing resonances from being split across separate integral regions; a 0.04-ppm bucket width and a 50% looseness factor were used. All spectra were normalized to total area excluding the water, urea, and glucose regions.
Confirmation of nucleotide identification.
Because of the relatively low concentration of nucleotides precluding definitive assignment by 2D spectra, LC-MS was used to confirm their presence. Urine samples were analyzed for N-methyl-2-pyridone-5-carboxamide (2PY) and N-methyl-4-pyridone-5-carboxamide (4PY) using an API-3000 HPLC/MS/MS (Applied Biosystems/MDS SCIEX, Foster City, CA) with Agilent 1100 binary solvent delivery system HPLC and 100 x 2.1 mm inner diameter, 4-m Synergi Polar column (Phenomenex, Torrance, CA). The lower limit of detection for all analytes was 1.0 µg/ml. Data were acquired and quantified using Analyst Version 1.4 (Applied Biosystems/MDS SCIEX, Foster City, CA) and SMS2000 Version 1.4 (GlaxoSmithKline, in-house).
Chemometric analysis of the data.
Multivariate data analysis was carried out using SIMCA-P+ 10.0 (Umetrics, Umea, Sweden). Data were mean centered and Pareto scaled (weighted by 1/
of the mean-centered variable) before analysis. Pareto scaling is a compromise between mean centering, which may fail to pick out small changes in metabolite concentrations, and scaling to unit variance, which gives equal weight to baseline imperfections, noise, and defined signals in the NMR spectrum.
Principal components analysis (PCA) was used to examine inherent clustering and correlations within the data. A principal component (PC) is a weighted linear combination of each of the original NMR variables so that the original data matrix is compressed into a smaller number of "latent variables," typically three to four PCs for NMR data. The weight given to each variable within a PC describes how influential that variable is in relation to the other variables. The information that is not captured in the first PC forms the residuals through which the second PC is calculated; all PCs are mutually orthogonal.
Partial least squares (PLS) techniques were used to assess correlation between the observed NMR data and other factors such as age. A PLS model is expressed as a set of X-scores (NMR spectral regions) and Y-score vectors (e.g., age) with corresponding X- and Y-weight vectors for a set of PLS model dimensions. Each dimension expresses a linear relation between an X-score and Y-score vector with the weight vectors describing how the X- and Y-variables are combined to give the X- and Y-score vectors. The model corresponds to fitting the lower-dimensional line, plane, or hyperplane simultaneously to the X- and Y-data as points in multidimensional space that best approximate the original data. The PLS model can be used to estimate the Y-variables corresponding to a given set of X-variables, (e.g., estimate the age of an organism given an NMR spectrum). In PLS-discriminant analysis (PLS-DA), dummy variables representing the class (e.g., disease) of each sample form the Y-matrix. PLS-DA was used for classification of samples where the PCA models were dominated by effects such as species or gender differences.
To assess which metabolite regions were responsible for a given classification, a twofold procedure was used combining both multivariate data analysis (MVDA) and univariate data analysis (UVDA) approaches. For each PLS-DA component responsible for class separation, cross-validation jack-knifing was used to calculate the magnitude and standard error for each variable coefficient, and Student's t-test was used to determine which coefficients differed significantly from zero. For the MVDA component of the combined analysis, a coefficient scored 0.15 if it had P < 0.20, 0.50 if P < 0.10, and 1.00 if P < 0.05.
Univariate statistical tests were carried out on bucketed NMR spectra. Each bucket was treated as an independent variable with two statistical samples: control and disease. Evidence for a statistically significant difference between the distributions of the two samples was based on the parametric Student's t-test (assumes normal distribution) and F-test (comparison of variances), and the nonparametric Kruskal-Wallis (generalized analysis of variance with no assumption of distribution form) and Kolmogorov-Smirnov tests (comparison of distribution shape). A threshold P value of 0.05 was used for each test. In the combined analysis, for each test passed the variable scored 0.25, giving a maximum score of 1.00. For the combined MVDA/UVDA analysis, metabolite regions scoring above 1.0 out of a maximum of 2.0 were considered significant.
For all models, disease/control classification was verified using a Cooman's plot (13), and the model was validated by random permutation of the Y-matrix (see Supplemental Materials) (14).
| RESULTS |
|---|
|
|
|---|
|
Similar metabolites were shown to be important for disease/control discrimination for both males and females when each gender was modeled separately (data not shown). This finding confirms earlier conclusions from the PCA and PLS-DA of the full data set where disease/control separation was the dominant source of variation (PC1) with gender separation orthogonal (PC2). This absence of disease/gender correlation was confirmed by the lack of any disease/gender interaction observed in the Kruskal-Wallis test results.
PLS models were built to describe the age trend for the wildtype/heterozygous samples and for the db/db samples for both male and female samples together (Supplemental Fig. 1, D and E) and separately (not shown). For each PLS model, the mouse data for either db/db or the wildtype/heterozygous group from all age groups combined were used to build the models. To examine whether these age trends were the same for the wildtype/heterozygous and homozygous groups, each model was used to predict the ages of the group not used to build the PLS model. The models built using the wildtype/heterozygous data sets were able to predict the ages of the wildtype/heterozygous samples but were unable to predict the ages of the db/db mice (data not shown) and vice versa. This indicated that the change in urinary metabolite profile with age is different for the db/db mice and wildtype/heterozygous mice. The urinary excretion of ß-hydroxybutyrate and acetone increased with age in db/db mice, whereas the TCA cycle intermediates citrate, 2-oxoglutarate, and fumarate as well as allantoin, creatine, N-methylnicotinamide (NMN amide), hippurate, meta-hydroxyphenyl-propionic acid (mHPPA), and indoxyl sulfate decreased with age.
Analysis of the Zucker rat data set.
The urinary profiles of the Zucker fatty and lean rats were examined by MVDA and UVDA of 1H-NMR spectra excluding the glucose resonances (Fig. 2A). For male Zucker rats, PCA (Q2 = 0.60, R2 = 0.7, 4 components) and PLS-DA (R2 = 0.95, Q2 = 0.94, 2 components) models separated the lean (fa/+) and fatty (fa/fa) rats (Fig. 2B and Supplemental Fig. 2C). While the disease separation dominated PC1, PC2 characterized differences in day/night metabolism. The fa/fa rats were distinguished clearly from fa/+ rats by increased concentrations of a range of metabolites including
-hydroxy-n-butyrate, lactate, fumarate, taurine, betaine, citrate, free fatty acids, DMA, N,N-dimethylglycine (DMG), cis-aconitate, valine, lysine, glutamine/glutamate, succinate/malate, and formate, while 2-oxoglutarate, phenylacetyl-glycine (PAG), hippurate, allantoin, unassigned N-acetyl groups (including glycoprotein), proline, ornithine, and creatinine were decreased (Supplemental Table 5).
|
Human diabetes patients.
PLS-DA models (R2 = 0.67, Q2 = 0.51, and 3 components; males and females excluding outliers; Fig. 3B and Supplemental Fig. 3C) of the NMR urine profile from T2DM patients compared with that of healthy subjects (Fig. 3A) showed a large number of metabolites contributed to the separation of the models. These metabolites included amino acids, tryptophan/tryptamine, 2-oxoisovalerate, alanine, ornithine, leucine, isoleucine, valine, and histidine. In addition, increases in relative concentration were observed for acetoacetate, acetate, n-butyrate,
-hydroxy-n-butyrate, DMA, DMG, NMN amide, and NAA, with decreases for creatinine, N-acetyl groups (including glycoproteins), N-methylnicotinate (NMN acid), aminohippurate, hippurate, PAG, allantoin, fumarate, and succinate. MVDA and UVDA showed that
60 chemical shift regions were responsible for robust classification of diabetic patients (Supplemental Table 6).
|
| DISCUSSION |
|---|
|
|
|---|
|
Kidney function in T2DM.
The creatinine depletion observed in urine from the animal models and T2DM patients may have arisen because of a variety of reasons such as changes in muscle mass, creatinine reabsorption, cell leakage, and changes in caloric intake. However, a decline in the renal glomerular filtration rate (GFR) associated with diabetic nephropathy (8, 29, 39) would not only explain this observation but also a variety of other metabolite changes. For example, reduced urinary allantoin excretion is also indicative of reduced GFR as allantoin is not reabsorbed across the proximal tubule, and thus its urinary concentration is thought to accurately reflect glomerular filtration (5, 19). Reduced urinary creatinine has also been observed in older Zucker rats (>13 wk) exhibiting albuminuria and glomerulosclerosis (41) with altered renal tubular function and morphology (dilation with proteinaceous casts, loss of functional microvilli on epithelium) (6). Hypertension and acquired dopamine D1 receptor dysfunction have also been reported in Zucker rats (3). Reduced GFR and substantial glomerular pathology, including mesangial matrix expansion and albuminuria, were reported for the db/db mouse along with nephropathy (28) and a decline in creatinine clearance after 5 mo of age (1).
Thus the observed changes in urinary metabolite concentration in the rodent models suggest early stage renal effects. However, the selection criteria employed for human subjects ensure no significant renal dysfunction for this data set. Because the changes in urinary metabolite concentration in the rodent models mirror those in the human samples, it is unlikely that the changes observed in the mouse and rat urine result solely from renal dysfunction, and most likely only represent early impairment of renal dysfunction.
Proteinuria in db/db mouse.
High levels of protein were present in some of the mouse samples, especially in the db/db samples. The time scale of this effect and its presence also in C57/B16 wildtype mice suggested that it was independent of altered glomerular function (see above). Because of the broad resonances associated with these proteins, it was not possible to either identify individual proteins or reliably quantify these resonances, although their contribution was higher in male diabetic mice compared with all the other groups.
TCA cycle and fatty acid ß-oxidation intermediates.
Increases in relative concentrations of the TCA cycle intermediates citrate, malate, fumarate, and cis-aconitate were observed in the fa/fa rats and db/db mice but not in humans. Urinary excretion of these metabolites has previously been correlated inversely and nonspecifically with renal and hepatic toxicity as well as general stress (e.g., calorie restriction or neurological disease) (4, 12, 20, 24). The db/db mice and fa/fa rats have glycosuria, and it is possible that the correlated increases in urinary concentrations of TCA cycle intermediates reflect either systemic stress produced by hyperglycemia or local effects on kidney tubular transport.
Conversely, in human T2DM patients, relative decreases in malate, fumarate, and succinate were observed. Most human subjects exhibited more modest hyperglycemia compared with the two rodent models, suggesting that they had good dietary control of diabetes and consequently reduced manifestations of accompanying pathologies (hepatic gluconeogenesis, insulin insensitivity, and renal dysfunction).
In diabetic patients and the animal models of T2DM, ketone bodies and fatty acids were increased in concentration. In humans this included acetate, acetoacetate, n-butyrate,
-hydroxy-n-butyrate, and ß-hydroxybutyrate, which result from increased ß-oxidation (45). In fa/fa rats, there was increased urinary excretion of short/medium chain fatty acids, although acetate and acetoacetate were decreased. These effects indicate an impairment of adipose tissue storage of circulating fatty acids and inhibition of hepatic fatty acid esterification. Both arise from insulin insensitivity causing increased concentrations of nonesterified fatty acids in blood plasma (34) and ultimately increased partial ß-oxidation of fatty acids in the liver and skeletal muscle producing short chain fatty acids and ketone bodies (16).
Amino acid metabolism.
End-stage diabetes, especially type 1 diabetes, is typified by the conversion of protein into glucose via gluconeogenesis. In the db/db mouse, observed metabolic perturbations indicated decreased excretion of amino acids, possibly reflecting increased clearance of these metabolites by the liver for gluconeogenesis, exacerbating the control of hyperglycemia. While there were comparable decreases in leucine, isoleucine, valine, histidine, and tryptophan in diabetic patients, an increase in glutamine and ornithine excretion was detected, suggesting that the perturbations in amino acid metabolism were more complex than simply an increase in gluconeogenesis in the liver.
Because of resonance overlap with glucose, the changes in taurine excretion could only be inferred using visual inspection. However, with this approach, taurine was observed to be increased in both animal models and the human patients. Although urinary taurine concentration is primarily regulated by renal reabsorption, the reabsorption is in turn regulated by taurine availability. Thus the increased excretion of taurine may arise from altered renal reabsorption of taurine as a result of reduced GFR or possibly as a general stress response, particularly following damage to the liver (22). Hypertaurinuria has also been detected following xenobiotic-induced perturbations in hepatic protein synthesis and other hepatotoxicities (4, 42).
Nucleotide metabolism.
In T2DM patients, a relative increase in NMN amide and 2PY with a decrease in NMN acid was observed with similar changes in the db/db mouse and fa/fa rat. NMN amide is involved in the tryptophan-NAD+ pathway, which supplies pyridine nucleotides to the liver (44). NMN amide can be further metabolized to 2PY and 4PY. 2PY formation appears to predominate over 4PY in humans and vice versa in rodents. NMN amide and 4PY have been suggested as urinary and plasma biomarkers of peroxisome proliferation in rats (10, 36). Plasma 2PY also increases with age in humans and in those with renal problems, possibly as a consequence of decreased renal excretion (38). Allantoin, a degradation product of nucleotide metabolism, is decreased in all three species, suggestive of a profound change in nucleotide metabolism during T2DM. It is interesting to note the commonalities between the perturbations reported in nucleotide metabolism in peroxisome proliferation and those observed here for T2DM. Peroxisome proliferation can only be confirmed by electron microscopy of the liver tissue, which was not carried out on either of the rodent models and of course is precluded from the routine analysis of human patients.
Analysis of the correlation of urinary nucleotide concentrations with clinical measures of T2DM revealed that changes in nucleotide metabolism reflect the direct metabolic consequences of the disease and are not the result of secondary complications (see Supplemental Materials for a detailed discussion of the correlation analysis). This included a negative correlation between the ratio of 2PY/NMN acid with HDL (R2 = 0.38 ± 0.12 for men and R2 = 0.49 ± 0.08 for women), although this correlation showed a distinct gender bias. Interestingly, the correlation between adiponectin and 2PY showed a positive correlation in males but a negative correlation in females, underlying the distinct interaction between gender and disease for T2DM.
Age progression and gender effects on metabolism in diabetes.
By use of PLS to model metabolic changes with aging, the mouse and rat diabetic animals demonstrated different aging trends compared with control animals. The older diabetic animals were characterized by increased ß-hydroxybutyrate and acetone and decreased citrate, 2-oxoglutarate, and fumarate relative to younger diabetic animals, with these trends not detected in the control group, suggesting that, in the older animals, increased insulin insensitivity has led to mitochondrial metabolism of acetyl-CoA increasingly via the 3-hydroxy-3-methyl-glutaryl-CoA cycle compared with the TCA cycle. Older db/db mice were characterized by decreased NMN amide, hippurate, mHPPA, and indoxyl sulfate relative to younger mice. Other effects observed include diurnal metabolite changes, consistent with the hormonal and sleep/wake cycles in Zucker rats.
A combined model of T2DM.
A single PLS-DA model of T2DM was built using the NMR data from the mouse and rat models and from the human T2DM patients (Fig. 4B). In the scores plot, there is distinct species separation between the mouse, rat, and human samples with the separation being primarily due to relative changes in concentration of 2-oxoglutarate/creatine, citrate, DMA, and creatinine, with rats exhibiting a relative decrease in 2-oxoglutarate/creatine, citrate, and DMA and humans an increase in creatinine. However, all three groups demonstrate clear disease/wildtype separation along a common axis, indicating a similar disease process in each of the three species, with the differences among species being perpendicular to the differences associated with disease in our analysis.
The pathways highlighted as being affected by T2DM according to our metabolomic analysis (Fig. 4A) tally closely with gene transcript changes detected in the liver of db/db mice during treatment with metformin, a widely used hypoglycemic agent used in the treatment of T2DM acting to reduce hepatic gluconeogenesis and increase peripheral sensitivity to insulin (23). Heishi et al. (23) report changes in glycolysis/gluconeogenesis, TCA cycle, pyruvate metabolism, bile acid biosynthesis, lysine degradation, arginine and proline metabolism, and tryptophan metabolism; all of these pathways were similarly perturbed by T2DM in this study. Along with the expected changes in hepatic glycolysis/gluconeogenesis, TCA cycle, and pyruvate and fatty acid metabolism, significant changes in hepatic amino acid metabolism were observed including tryptophan metabolism. Tryptophan is a precursor for nicotinate from which NMN amide, 2PY, and 4PY are produced. A summary of the metabolite changes observed in all three models of T2DM is given in Table 1.
|
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
* R. M. Salek and M. L. Maguire contributed equally to this work. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. L. Griffin and A. Vidal-Puig Current challenges in metabolomics for diabetes research: a vital functional genomic tool or just a ploy for gaining funding? Physiol Genomics, June 1, 2008; 34(1): 1 - 5. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |