Abstract
Islet-level oxidative stress has been proposed as a trigger for type 1 diabetes (T1D), and release of cytokines by infiltrating immune cells further elevates reactive oxygen species (ROS), exacerbating β cell duress. To identify genes/mechanisms involved with diabetogenesis at the β cell level, gene expression profiling and targeted follow-up studies were used to investigate islet activity in the biobreeding (BB) rat. Forty-day-old spontaneously diabetic lymphopenic BB DRlyp/lyp rats (before T cell insulitis) as well as nondiabetic BB DR+/+ rats, nondiabetic but lymphopenic F344lyp/lyp rats, and healthy Fischer (F344) rats were examined. Gene expression profiles of BB rat islets were highly distinct from F344 islets and under-expressed numerous genes involved in ROS metabolism, including glutathione S-transferase (GST) family members (Gstm2, Gstm4, Gstm7, Gstt1, Gstp1, and Gstk1), superoxide dismutases (Sod2 and Sod3), peroxidases, and peroxiredoxins. This pattern of under-expression was not observed in brain, liver, or muscle. Compared with F344 rats, BB rat pancreata exhibited lower GST protein levels, while plasma GST activity was found significantly lower in BB rats. Systemic administration of the antioxidant N-acetyl cysteine to DRlyp/lyp rats altered abundances of peripheral eosinophils, reduced severity of insulitis, and significantly delayed but did not prevent diabetes onset. We find evidence of β cell dysfunction in BB rats independent of T1D progression, which includes lower expression of genes related to antioxidative defense mechanisms during the pre-onset period that may contribute to overall T1D susceptibility.
Introduction
Type 1 diabetes (T1D) is an autoimmune disease characterized by immunocyte infiltration of the pancreatic islets (insulitis) and destruction of the insulin-secreting β cells. Diabetes in the biobreeding (BB) rat exhibits many similarities to the human disease, including onset during puberty, insulitis, and lifelong dependency on exogenous insulin (Scott 1990, Crisa et al. 1992).
Diabetes in BB rats is polygenic (Wallis et al. 2009). The MHC (insulin-dependent diabetes mellitus locus 1, Iddm1) contributes the largest genetic risk as it does in humans and in the NOD mouse. In rats, this is the HLA-DQB1 homolog, RT1-B, specifically the RT1u haplotype (Colle 1990, Awata et al. 1995). The DR+/+ and DRlyp/lyp congenic BB rat lines differ by the Iddm2 region (Jacob et al. 1992), where the T cell lymphopenia of the DRlyp/lyp rat arises from a single nucleotide deletion in the Gimap5 gene (Hornum et al. 2002, MacMurray et al. 2002). T1D develops spontaneously in 100% of DRlyp/lyp rats at ∼60 days of age and is elicited through a deficiency in CD4+CD25+ regulatory T (TREG) cells, as adoptive transfer of this population prevents T1D (Lundsgaard et al. 2005, Poussier et al. 2005). DR+/+ rats possess a wild-type Gimap5 and do not spontaneously develop T1D; however, it is possible to induce diabetes through depletion of TREG cells (Mordes et al. 1996, Zipris et al. 2003). Fischer (F344) rats, which normally possess the RT1lv1/lv1 MHC, do not develop T1D even after introgression of Iddm1 and/or Iddm2 (Jacob et al. 1992, Moralejo et al. 2003), indicating that they lack this additional susceptibility. Thus, predisposition for T1D independent of Iddm1 and Iddm2 exists in the BB rat that is phenotypically manifest upon loss of immune regulation.
Oxidative stress is associated with many pathological states and has been implicated as a trigger for T1D (reviewed in Lenzen (2008)). Compared with other tissues, such as kidney or liver, islets possess lower levels of the antioxidant enzymes superoxide dismutase (SOD) and catalase (Lenzen 2008) and may be more susceptible to redox imbalances arising from overproduction of reactive oxygen species (ROS). Furthermore, the release of pro-inflammatory cytokines by infiltrating immune cells can elevate intra-islet ROS as exposure of rodent islets to interleukin 1β (IL1β) induces nitric oxide synthase (iNOS) expression and production of the free radical nitric oxide (Corbett et al. 1992, Eizirik et al. 1993). While cultured human islets express higher levels of antioxidant enzymes than do their rodent counterparts (Welsh et al. 1995), inherited deficiencies in antioxidant defense mechanisms may be relevant to T1D susceptibility as strain-dependent differences in islet sensitivity to IL1β (Reimers et al. 1996) and lower islet levels of catalase and SOD activity have been reported in BB relative to Wistar rats (Sigfrid et al. 2004).
Consistent with the concept that autoimmunity involves deficiencies in self-tolerance as well as target organ-specific factors, we have discovered that β cells of DRlyp/lyp and DR+/+ rats share an islet-level T1D susceptibility. This is reflected by the expression of the eosinophil and mast cell recruiting chemokine eotaxin, beginning by 40 days of age, providing a means for immunocyte recruitment and activation (Hessner et al. 2004, Geoffrey et al. 2006). Here, we further examine tissue-specific T1D susceptibility by comparing the islet transcriptomes of DR+/+ and DRlyp/lyp rats during the pre-onset period to those of F344 and F344lyp/lyp rats (Iddm2 introgressed onto the F344 background). We find BB rat islets under-express genes involved in metabolism of ROS, in particular numerous isoforms of glutathione S-transferase (GST), an important family of enzymes that participate in the detoxification of reactive electrophilic compounds by catalyzing their conjugation to glutathione. Administration of N-acetyl cysteine (NAC), a thiol compound that acts directly as a free radical scavenger and a precursor in glutathione synthesis, to DRlyp/lyp rats, decreased abundances of peripheral eosinophils, reduced severity of insulitis, and delayed but did not prevent diabetes onset.
Materials and methods
Animals and NAC treatment
BB (Bieg et al. 1998), F344, and F344lyp/lyp (Moralejo et al. 2003) rats were maintained at the Pacific Northwest Diabetes Research Institute, The University of Washington-Seattle, and The Medical College of Wisconsin. BB DR+/+ and DRlyp/lyp animals were propagated by mating DRlyp/+ breeders and genotyped as described (MacMurray et al. 2002), housed under specific pathogen-free conditions with standard light/dark cycles, and were fed a regular diet and water ad libitum. All federal (http://grants1.nih.gov/grants/olaw/references/phspol.htm) guidelines for use and care of laboratory animals were followed, and all protocols were approved by the respective Institutional Animal Care and Use Committees (IACUC).
DRlyp/lyp and DR+/+ rats were treated with NAC daily (200 mg/kg, i.p.; Sigma), dissolved in 0.9% sterile NaCl, and neutralized to pH 7 with 1 M NaOH beginning at weaning. Age-matched DRlyp/lyp and DR+/+ rats were sham-treated with equal volumes of 0.9% NaCl.
Animal phenotyping and tissue collection
Weight and blood glucose were measured from day 40 of age until DRlyp/lyp rats developed T1D (blood glucose levels ≥250 mg/dl measured with an Ascensia Elite XL glucometer; Bayer). Animals underwent weekly oral glucose tolerance tests starting from 6 weeks of age. Following a 6-h fasting, glucose solution (2 g/kg body weight) was administered by gavage, followed by blood glucose measurements at 0, 15, 30, 60, 90, and 120 min. Plasma was separated from tail vein blood and stored at −20 °C until assayed for C-peptide by RIA (Linco, St Charles, MO, USA). At onset, DRlyp/lyp and age-matched DR+/+ littermates were killed and pancreata were processed for histological analysis and measurement of insulin content.
Differential cell counts of peripheral blood were performed by an automated cell counter. Percentages of the different leukocyte populations were determined: 50 μl blood was incubated for 15 min in 10 ml red blood cell lysis solution, leukocytes were harvested by centrifugation, spread onto glass slides, stained with Giemsa–Grunwald solution, and percentages of lymphocytes, polymorphonuclear cells, and monocytes were determined (800–1200 total cells per sample).
Flow cytometric analysis of peripheral blood lymphocytes was accomplished as described (Fuller et al. 2006). TREG were identified by surface staining with anti-CD4 (clone OX-35) and anti-CD25 (clone OX-39) followed by intracellular staining with anti-Foxp3 (clone FJK-16s) using an intracellular staining kit from eBioscience (San Diego, CA, USA).
Islet gene expression profiling
Islets were isolated from normoglycemic 40-day-old rats as described (Sweet et al. 2004) and total RNA was extracted with Trizol (Invitrogen). RNA (∼100 ng) was amplified/labeled (Affymetrix two-cycle cDNA synthesis kit, Affymetrix, Santa Clara, CA, USA) and then hybridized to the Affymetrix RG230 2.0 array as per the manufacturers' protocol. Image data were quantified with Affymetrix Expression Console Software and normalized with Robust Multichip Analysis (www.bioconductor.org) to determine signal log ratios. Data files have been deposited in The National Center for Biotechnology Information Gene Expression Omnibus (accession number GSE20214). The statistical significance of differential gene expression was derived through a Student's t-test and false discovery rates (FDR) were determined with Significance Analysis of Microarrays (SAM) software as described (Tusher et al. 2001). More stringent statistical criteria were not applied, as relevant genes and pathways would be investigated further by quantitative RT-PCR (qRT-PCR) and targeted follow-up studies. Ontological pathway analysis was performed with the Database for Annotation, Visualization, and Integrated Discovery (DAVID; Hosack et al. 2003). Hierarchical clustering was conducted with Genesis (Sturn et al. 2002).
Quantitative RT-PCR
Total RNA was extracted from brain, muscle, liver, or cultured cells using TRIzol reagent, and qRT-PCR for Gstm2, Gstm7, Gstt1, Sod2, Sod3, Ccs, Cd68, and FcεR1 was performed as described (Pfaffl 2001, Wang et al. 2008). Primer designs and reaction performance parameters are provided in Table 1.
Quantitative real-time RT-PCR performance parameters and primer designs
Gene name | UniGene | Slopea | R2b | Ec | Primer sequenced |
---|---|---|---|---|---|
Superoxide dismutase copper chaperone (Ccs) | Rn.12311 | −3.41 | 0.999 | 1.96 | 5′-GTGATGGGCTCACTATCTGG-3′ |
5′-CCCTCTGCAGGAAGATGTT-3′ | |||||
Superoxide dismutase 2 mitochondrial, mRNA (Sod2) | Rn.10488 | −3.45 | 0.9981 | 1.95 | 5′-GGCTTGGCTTCAATAAGG-3′ |
5′-GAAGATAGTAAGCGTGCTCC-3′ | |||||
Superoxide dismutase 3, extracellular (Sod3) | Rn.10358 | −3.20 | 0.9973 | 2.05 | 5′-AATGTGCGCAGATAACTCTC-3′ |
5′-CAAGGGGAAGACAAACGA-3′ | |||||
Glutathione S-transferase, mu type 2 (Yb2) (Gstm2) | Rn.625 | −3.47 | 0.9994 | 1.94 | 5′-GATTCGTGTGGACGTTTTG-3′ |
5′-TCAAAGTCAGGGCTGTAGCA-3′ | |||||
Glutathione S-transferase mu 7 (Gstm7) | Rn.6036 | −3.17 | 0.9983 | 2.07 | 5′-CCCTCGAGCTCATGTGATTG-3′ |
5′-CGGTAGGAGCGGGTTACAG-3′ | |||||
Glutathione S-transferase theta 1 (Gstt1) | Rn.11122 | −3.18 | 0.9979 | 2.06 | 5′-CCAGTCTTTGAAGGGCGTCC-3′ |
5′-GGGCGCACAGTCGTGTAATG-3′ | |||||
Cd68 molecule (Cd68) | Rn.12478 | −3.56 | 0.9922 | 1.91 | 5′-TGCTTCTCTTGCGCCAGT-3′ |
5′-GGGCTGGTAGGTTGATTGTC-3′ | |||||
Fc fragment of IgE, high affinity I, receptor for; γ polypeptide (FcεR1g) | Rn.201810 | −3.38 | 0.9964 | 1.97 | 5′-AGAAACCACCCCAATAGCTT-3′ 5′-GTGGATATGGGGACATTGA-3′ |
Slope of standard curve.
Linear regression of standard curve.
Reaction efficiency, determined as previously described (Pfaffl 2001).
Top sequence is forward primer, bottom is reverse.
Rat insulinoma cell culture and cytokine exposure
RINm5f cells (ATCC, Rockville, MD, USA) were maintained in RPMI 1640 medium (Mediatech, Inc., Manassas, VA, USA) containing 10% fetal bovine serum (Atlanta Biologicals, Lawrenceville, GA, USA) and Pen/Strep (Lonza, Walkersville, MD, USA). Medium was changed every 2–3 days and cells were subcultured as required. Cells were allowed to adhere overnight in medium containing 10% charcoal dextran-treated serum. The medium was changed and cells were stimulated with cytokines for 24 h: IL1β (2 ng/ml), IFNγ (30 ng/ml), and TNFα (10 ng/ml) (R&D Systems, Inc., Minneapolis, MN, USA). Following treatment, medium was aspirated and cells were harvested for qRT-PCR.
Histological studies
Pancreata were processed for histology as described (Hessner et al. 2004, Bogdani et al. 2005, Geoffrey et al. 2006). Insulitis was evaluated in hematoxylin- and eosin-stained sections (five to seven sections per pancreas; 200–400 islets) and independently scored by an investigator unaware of the animals' genotype as described (Fuller et al. 2009): Grade 0, no infiltration; Grade 1, leukocytes around ducts and vessels only; Grade 2, leukocytes around islets; Grade 3: leukocytes inside islets without change in β cell morphology; Grade 4, leukocytes inside islets with distorted β cell morphology or islets devoid of β cells.
Consecutive sections were immunostained using antibodies targeting CD45 (leukocyte common antigen), CD3, CD4, CD8 (T cell subsets), CD45R (B-cells), and CD68 (macrophages/monocytes); >1000 cells/pancreas were counted; and percentages were determined. Pancreatic sections were also stained with anti-glucagon (Sigma) and anti-GST-mu (GSTM, Abcam, Cambridge, MA, USA) antibodies as described (Hessner et al. 2004, Geoffrey et al. 2006). Slides were imaged on a Nikon E600 system (Nikon USA, Melville, NY, USA). MetaMorph version 6.3r3 software (Universal Imaging, Buckinghamshire, UK) was used to trace the islet and to capture the fluorescence intensity within the tracing for a minimum of ten islets per section, one section per rat and n≥4 rats per strain.
GST, glutathione reductase, and glutathione measurements
Assays for measurement of plasma GST activity, glutathione reductase levels, and glutathione levels were conducted as per the manufacturers' instructions (Arbor Assays, Ann Arbor, MI, USA). All samples were analyzed in duplicate and results were analyzed using four-parameter logistic curve fitting.
Statistical analysis
Data are expressed as mean±s.d. of n independent experiments unless otherwise noted. The statistical significance of differences in glycemic levels between experimental and control groups was calculated by the Kruskal–Wallis test. Changes in percentage of different cell populations and insulitis score were tested by the Mann–Whitney U test. P<0.05 was considered significant. Survival was analyzed with the Kaplan–Meier method.
Results
Islet gene expression profiling
Given that β cells of 40-day-old BB rats show evidence of eotaxin expression, before detection of insulitis in DRlyp/lyp rats (Hessner et al. 2004, Geoffrey et al. 2006), we compared the islet transcriptomes of normoglycemic 40-day-old DRlyp/lyp and DR+/+ rats vs F344 and F344lyp/lyp rats with the goal of identifying pathways/mechanisms associated with T1D susceptibility common to BB rats, independent of Gimap5 status. Hierarchical clustering of the significantly regulated probe sets illustrates the distinctiveness between these lineages at 40 days of age (Fig. 1A). Genes over-expressed and under-expressed in BB islets (n=2084 and 2545 probe sets respectively; Supplementary Table 1, see section on supplementary data given at the end of this article) were independently evaluated for biological pathway enrichment using the DAVID annotation tool (Supplementary Table 2, see section on supplementary data given at the end of this article). Selected pathway terms that functionally discriminate the islet expression profiles between the BB and the F344 strains are tabulated in Table 2.
Pathway analysis of probe sets differentially regulated by BB vs F344 islets. Genes over-expressed in the BB strains (n=2084 probe sets) and those under-expressed in the BB strains (n=2545 probe sets) were independently evaluated for biological pathway enrichment (GO molecular function and KEGG) using the DAVID annotation tool. Selected terms that functionally discriminate the islet expression profiles between the BB vs non-BB strains are shown, a complete tabulation is provided in Supplementary Table 2
Annotation category | Term | Count | P value |
---|---|---|---|
Over-represented in BB (n=2084 probe sets) | |||
GO molecular function | Protein serine/threonine kinase activity | 64 | 1.3×10−8 |
GO molecular function | Transcription regulator activity | 125 | 6.6×10−6 |
GO molecular function | Zinc ion binding | 159 | 1.7×10−5 |
GO molecular function | Transcription repressor activity | 41 | 2.3×10−5 |
KEGG pathway | Wnt signaling | 28 | 3.9×10−5 |
KEGG pathway | MAPK signaling | 42 | 6.1×10−5 |
KEGG pathway | TGFβ signaling | 16 | 2.6×10−3 |
GO molecular function | SMAD binding | 11 | 2.6×10−3 |
GO molecular function | Guanyl-nucleotide exchange factor activity | 19 | 4.5×10−3 |
GO molecular function | GTPase regulator activity | 35 | 7.1×10−3 |
GO molecular function | MAP kinase activity | 6 | 1.9×10−2 |
KEGG pathway | mTOR signaling | 9 | 5.5×10−2 |
Under-represented in BB (n=2545 probe sets) | |||
KEGG pathway | Oxidative phosphorylation | 51 | 1.5×10−10 |
KEGG pathway | Fructose and mannose metabolism | 15 | 3.7×10−4 |
GO molecular function | Intramolecular oxidoreductase activity | 15 | 4.1×10−4 |
KEGG pathway | Glutathione metabolism | 17 | 1.4×10−3 |
KEGG pathway | Proteasome | 17 | 1.4×10−3 |
GO molecular function | NADH dehydrogenase activity | 9 | 4.2×10−3 |
KEGG pathway | Pentose phosphate | 10 | 9.3×10−3 |
GO molecular function | Antioxidant activity | 13 | 1.3×10−2 |
KEGG pathway | Pyruvate metabolism | 13 | 2.2×10−2 |
GO molecular function | Cytochrome c oxidase activity | 8 | 3.1×10−2 |
GO molecular function | NADH dehydrogenase (ubiquinone) activity | 7 | 3.2×10−2 |
KEGG pathway | Galactose metabolism | 9 | 3.9×10−2 |
The P value defines the significance of the association of a particular biological process with the gene list analyzed.
Genes related to antioxidant activity and the glutathione pathways are under-represented in islets of BB compared with F344 rats
Glutathione metabolism and antioxidant activity were identified as enriched pathway terms when analyzing genes under-expressed in DRlyp/lyp and DR+/+ relative to F344 and F344lyp/lyp islets (Table 2). Within these terms were 27 unique genes related to oxidative defense mechanisms. Among these were members of the GST superfamily (Gstk1, Gstp1, Gstt1, Mgst3, Gstm2, Gstm4, Gstm5, and Gstm7), which function in the detoxification of electrophilic compounds through their conjugation to reduced glutathione. BB rat islets showed lower expression of γ-glutamyltransferases (Ggt1, Ggt5, and Ggt7) and 5-oxoprolinase (Oplah) important in glutathione synthesis/metabolism. Also under-expressed in BB islets were glutathione peroxidase 4 (Gpx4), a phospholipid hydroperoxidase that protects cells against membrane lipid peroxidation through the reduction of hydrogen peroxide, organic hydroperoxides, and lipid peroxides; peroxiredoxin 4 (Prdx4) and Prdx5, peroxidases that respectively control oxidative stress within the endoplasmic reticulum and mitochondria by reducing hydrogen peroxide; as well as superoxide dismutase 2 (Sod2, mitochondrial MnSOD) and Sod3 (extracellular Cu-ZnSOD). Apolipoprotein E (Apoe), which possesses antioxidant properties and can be downregulated under conditions of oxidative stress (Espiritu & Mazzone 2008), was under-expressed by BB rat islets. Ornithine decarboxylase (Odc1) and spermidine synthase (Srm), both important in the synthesis of the antioxidant polyamine, spermidine, were also under-expressed in BB islets. The overall pattern of reduced expression levels of these and other genes related to ROS defense/metabolism are illustrated in Fig. 1B.
Examination of GST expression in islets and other tissues
qRT-PCR of the independent islet RNA isolations from each rat used to generate the islet RNA pools for each strain in the array studies confirmed under-expression of Gstm2, Gstm7, Gstt1, Sod2, Sod3, and copper chaperone for SOD (Ccs, a metallochaperone required for incorporation of copper into Sod1) in BB rat islets. Smaller expression differences were observed in day 40 liver, brain, and skeletal muscle of DRlyp/lyp and DR+/+ vs F344lyp/lyp and F344 rats (Fig. 2A).
Immunofluorescent staining was used to localize and examine relative GST protein levels in pancreatic sections of 5-week-old DRlyp/lyp, DR+/+, F344, and F344lyp/lyp rats. Consistent with the array and RT-PCR expression measurements, significantly lower fluorescence indicative of reduced protein GST protein levels was observed in islets of DRlyp/lyp and DR+/+ rats (Fig. 2B).
Reduced plasma GST activity in BB rats
Significantly lower GST activity levels were detected in plasma of DRlyp/lyp and DR+/+ rats compared with F344 and F344lyp/lyp rats (Fig. 3A). To exclude the possibility that lower GST activity was due to a reduced ability to convert oxidized glutathione to reduced glutathione or a lower level of substrate, plasma glutathione reductase activity as well as total and free glutathione levels were measured. Neither glutathione reductase activity nor the percentage of free glutathione was significantly different between strains (Fig. 3B and C).
Inflammatory processes as a basis for reduced islet GST expression
As the under-expression of genes related to antioxidant defense was observed only in islets, we hypothesized that innate inflammatory processes may underlie reduced GST expression. DR+/+ rats do not develop insulitis nor T1D rats. By histological examination, DRlyp/lyp islets are free of immune cell infiltrates at day 40 but not at time points after 60 days (Hessner et al. 2004). To exclude the possibility that macrophages and/or mast cells, at abundances difficult to capture in routine histological studies, were present in the day 40 islets used for the gene expression studies, we respectively performed qRT-PCR for CD68 and FcεR1. Islets of 40-day-old DRlyp/lyp, DR+/+, F344lyp/lyp, and F344 islets expressed less CD68 and FcεR1 transcript compared with normoglycemic day 60 DRlyp/lyp islets. These differences reached statistical significance (P<0.05) in comparisons between day 60 DRlyp/lyp islets vs day 40 DRlyp/lyp and DR+/+ islets. While not of the level observed in day 60 DRlyp/lyp islets, significantly higher levels (P<0.05) of CD68 transcript were observed when comparing day 40 F344 islets to day 40 DRlyp/lyp and DR+/+ islets (Fig. 4A). These data are consistent with previous studies showing islets of day 40 DRlyp/lyp and DR+/+ rats free of immune infiltrates (Hessner et al. 2004). Further, these results support that the observed gene expression differences were not a consequence of the activity of infiltrating immune cells. However, it cannot be excluded that the activity of resident islet mast cells and macrophages may differ between the strains.
Our previous studies on BB rats revealed the presence of inflammatory processes in both DR+/+ and DRlyp/lyp rats at day 60 as well as evidence these processes begin as early as day 40 in DRlyp/lyp rats (Kaldunski et al. 2010). To examine the influence of inflammatory cytokines and cytokine-induced oxidative stress on the expression of GSTs, RINm5f rat insulinoma cells were treated with IL1β, TNFα, and IFNγ. Relative to untreated cultures, a 24-h cytokine treatment resulted in significant decreases in Gstm2 (−3.3-fold, P<0.001), Gstm7 (−3.6-fold, P<0.001), and Sod3 (−4.3-fold, P<0.05) expression while Sod2 expression was increased (+13.1-fold, P<0.01, Fig. 4B). Similar results were obtained when using the rat INS-1 cell line (data not shown). Given our previous early detection of inflammation in the BB rat (Kaldunski et al. 2010) and the sensitivity of Gst gene expression to cytokine exposure in RINm5f cells, it cannot be excluded that under-expression of GSTs in BB rat islets is the consequence of early pre-onset inflammatory processes. However, increased expression of Sod2 in response to IL1β, TNFα, and IFNγ suggests that β cell intrinsic factors also contribute to the overall under-expressed antioxidant genes.
Delay of BB rat T1D through treatment with NAC
Given the evidence for deficiencies in islet antioxidant capacity in BB rats during the pre-onset period relative to Fischer rats, we tested the effectiveness of the glutathione precursor NAC in delaying diabetes onset in DRlyp/lyp rats. Treatment with NAC resulted in a modest but significant delay in onset relative to sham-treated controls (P=0.033, Fig. 5A).
No differences in growth rate, weight, or organ weight (kidney, spleen, thymus, adrenal gland, and thyroid; P>0.05) were observed between NAC- and sham-treated DRlyp/lyp rats. NAC- and sham-treated DRlyp/lyp rats showed similar progressive weight gain that abruptly ceased a few days before diabetes onset; both groups exhibited a slight increase in blood glucose levels in the 2 days preceding the disease onset, being between 150 and 200 mg/dl, followed by levels >250 mg/dl (Fig. 5B). NAC treatment did not influence blood glucose levels in DR+/+ rats, which remained within the normal range values as in nontreated DR+/+ rats (data not shown).
As under-expression of antioxidant genes is evident in both BB sub-strains, OGTT, at days 42, 49, 56, and 63, was used to assess the effect of NAC treatment on islet function in DRlyp/lyp and DR+/+ rats. While DR+/+ animals remained glucose tolerant at all time points, NAC- and sham-treated DRlyp/lyp rats gradually became glucose intolerant. These data, aligned to diabetes onset (Fig. 5C), show that all DRlyp/lyp rats were glucose tolerant at −4 to −2 weeks. At −1 week, both NAC- and sham-treated DRlyp/lyp rats exhibited higher blood glucose levels than DR+/+ rats (P<0.05). While blood glucose levels were lower in NAC- vs sham-treated DRlyp/lyp rats at 30 through 90 min after the glucose challenge at −1 week, the difference did not reach statistical significance (P>0.5). OGTT did not reveal any change in glucose tolerance when comparing NAC- vs sham-treated DR+/+ rats.
At onset, circulating C-peptide levels and pancreatic insulin content in DRlyp/lyp rats were <10% of those of age-matched DR+/+ rats (P<0.05), and these parameters were not influenced by NAC treatment (data not shown). NAC treatment did not affect peripheral GST levels in DRlyp/lyp or DR+/+ rats. A trend of increased free glutathione levels (∼20%), while not reaching statistical significance, was observed in NAC-treated vs sham-treated DRlyp/lyp and DR+/+ rats (data not shown).
Effect of NAC on circulating leukocytes and insulitis
Flow cytometry and differential counts found DRlyp/lyp rats characteristically lymphopenic (Crisa et al. 1992, Groen et al. 1996, Greiner et al. 1997), possessing <15% the normal T cell count relative to DR+/+ rats. NAC treatment did not alter the abundances of T cells nor the relative abundances of lymphocytes, neutrophils, basophils, and monocytes (data not shown). It is well documented that DRlyp/lyp rats exhibit eosinophilia and islet eosinophil infiltration before diabetes onset (Kurner et al. 1986, Maruta et al. 1989, Eastman et al. 1991, Hessner et al. 2004). The eosinophilia is one of the several DRlyp/lyp phenotypes that parallels immunodysregulation polyendocrinopathy enteropathy X-linked syndrome (IPEX), a human disorder that arises through mutations to FoxP3, a transcription factor required for TREG cell development, and is characterized by autoimmunities that include diabetes (Geoffrey et al. 2006). An increase was observed in the proportion and the number of circulating eosinophils in both the NAC- and the sham-treated DRlyp/lyp rats in the 4 weeks before diabetes onset; however, it was significantly less in the NAC-treated group (Fig. 6A).
Pancreata of NAC- and sham-treated DRlyp/lyp rats at onset, as well as age-matched NAC- and sham-treated DR+/+ rats, were histologically evaluated. Inflammatory cells were scarce or absent in DR+/+ pancreata but were present in both the exocrine and the islet tissue of NAC- and sham-treated DRlyp/lyp rats. Nearly all islets of NAC- and sham-treated DRlyp/lyp rats exhibited insulitis; however, the severity differed between the two groups. Grade 4 insulitis was observed in 58 and 82% of NAC- and sham-treated DRlyp/lyp islets (Fig. 6B) respectively. More islets (34%) in the NAC-treated DRlyp/lyp rats were of Grade 2 or 3 inflammation compared with only 9% of the islets in the sham-treated DRlyp/lyp rats. The effect of NAC on the circulating eosinophil cell abundance was not associated with a reduction in the proportion of islets possessing eosinophils in the peri-islet inflammation nor with the number of eosinophils present in these infiltrates (data not shown).
At diabetes onset, macrophages (CD68+) made up the majority (>90%) of the inflammatory cells present in the exocrine tissue and penetrating the NAC- and sham-treated DRlyp/lyp islets. T cells (CD3+, CD4+, and CD8+) were present but were very few in number while B-cells (CD45R+) were not detected. The percentage of CD68+ cells in the pancreas was lower in NAC-treated DRlyp/lyp rats (Fig. 6C), being ∼50% that observed in the exocrine and endocrine tissues of sham-treated DRlyp/lyp rats. These data suggest that NAC may attenuate the severity of innate pancreatic inflammation, thereby reducing the rate of disease progression in BB rats.
Discussion
β Cells of lymphopenic and non-lymphopenic BB rats exhibit eotaxin expression by 40 days of age, before development of insulitis in DRlyp/lyp rats (Geoffrey et al. 2006). Eotaxin expression by β cells is a means of recruiting eosinophils, mast cells, and other immune cells that bear the receptor CCR3 to the islet. Given that, generally, healthy tissues do not express immune cell recruiting chemokines that promote inflammation, we have considered eotaxin expression a sign of underlying pathology in BB rats. Furthermore, β cell expression of eotaxin correlates with the ability of the BB lineage to develop T1D (either spontaneously in the DRlyp/lyp rat or through induction protocols in the DR+/+ rat). Therefore, we examined islet transcriptomes at day 40 to capture a time after evidence of islet dysfunction, but before evidence of insulitis and/or robust adaptive immunity, to identify factors in the BB rat that may contribute to T1D susceptibility at the β cell level.
Comparison of the day 40 islet transcriptomes of DRlyp/lyp and DR+/+ to F344lyp/lyp and F344 rats found the two lineages highly distinct. Confirming our previous studies, eotaxin transcript (Ccl11) was found more abundant in BB vs F344 islets (1.7-fold, P=0.033). Other transcripts consistent with immune cell recruitment/activation were differentially regulated in BB rat islets, including the neutrophil chemoattractant Cxcl1 (2.6-fold, P=0.016) and the IL22 antagonist Il22ra2 (−2.74-fold, P=0.002). Extending observations of lower catalase and SOD expression in BB rat islets (Pisanti et al. 1988, Sigfrid et al. 2004), gene expression related to multiple ROS defense mechanisms was significantly under-represented in BB vs F344 islets, in particular members of the GST superfamily. This raises the question whether the lower antioxidant gene expression is a consequence of islet-level pathological/inflammatory processes and/or genetic control.
CD68 and FcεR1 transcript levels were not different between day 40 DRlyp/lyp and DR+/+ islets and found lower in day 40 BB islets compared with day 40 F344 and F344lyp/lyp islets. Thus, greater numbers of infiltrated macrophages and mast cells were not co-isolated with the day 40 BB rat islets. While these and histological studies are consistent with an absence of insulitis at day 40, we cannot exclude that co-isolated resident myeloid cells, if any, possessed different activities between the two lineages. Further, we cannot exclude the possibility that activity of immunocytes at the islet periphery, lost during isolation, may be responsible for under-expression of antioxidant genes in BB islets. We found that cytokine treatment of RINm5f cells reduced expression of Gstm2, Gstm7, and Sod3, but the level of Gstt1 was not altered. Consistent with other reports, we also found that cytokine treatment of RINm5F cells increased Sod2 expression (Borg et al. 1992, Rieneck et al. 2000, Cardozo et al. 2001, Sigfrid et al. 2004). Taken together, these results suggest that both β cell intrinsic and extrinsic factors (such as innate immune cell-mediated inflammation) contribute to reduced expression of antioxidant genes in BB vs F344 islets.
Lower total plasma antioxidants have been reported in first-degree relatives of T1D patients (Rocic et al. 1997), supporting a possible genetic basis for under-expression of antioxidant genes that remains to be elucidated in humans and BB rats. Gstm2, Gstm4, and Gstm7 are targets of NRF2 (NFE212) (nuclear factor E2 p45-related factor 2) (Malhotra et al. 2010), a transcription factor that activates genes in response to oxidative stress (Chen & Kong 2004). Although we did not observe differential expression of Nrf2 or its negative regulator Keap1 in the array studies, Park7, a positive regulator of Nrf2 (Clements et al. 2006), was significantly under-expressed in BB rat islets. It is noteworthy that the four mu-class GSTs (Gstm2, Gstm4, Gstm5, and Gstm7) map to Iddm3 (cytoband 2q34 (Jacob et al. 1992, Klaff et al. 1999)), a region syntenic to Idd18.2 in mouse and the 1p13.3 T1D susceptibility region in human, and in crosses between DRlyp/lyp and F344 rats, a single F344 allele at Iddm3 has been sufficient for T1D protection, suggesting the presence of a recessive diabetes susceptibility trait at this locus.
Consistent with other animal studies that associate augmented antioxidant reserve/defense with amelioration of T1D (Prasad 2000, Piganelli et al. 2002, Szabo et al. 2002, Delmastro & Piganelli 2011), treatment of DRlyp/lyp rats with NAC delayed diabetes onset. Reported dosages of NAC vary widely, with our literature searches revealing dosages as high as 1 g/kg per day being administered to adult Wistar rats (Ozaras et al. 2003). Our optimization studies revealed that a dosage 400 mg/kg per day was not well tolerated by weanling rats, which showed slower weight gain during the first week of the NAC administration. While mild toxicity and/or off target effects cannot be completely excluded, at the 200 mg/kg per day NAC dosage, DRlyp/lyp rats were not different from the controls from 21 days of age, when the treatment was started, for either blood glucose or weight gain through at least 46 days of age. The NAC and saline-treated nonlymphopenic/nondiabetic BB control rats were not different in any of the parameters analyzed at any age (45–230 days). As reviewed (Delmastro & Piganelli 2011), ROS potentially acts at multiple levels of diabetogenesis: an islet-level autoimmune trigger, activating innate immune cells and initiating insulitis, promoting T cell-mediated adaptive immune responses, and impairing FOXP3 expression, thereby suppressing TREG activity in autoimmune disorders (Brahmachari & Pahan 2010). NAC treatment reduced the severity of the insulitic lesion in DRlyp/lyp rats but did not increase the frequency of TREG cells, nor did it alter the level of Foxp3 expression in these immune regulators (data not shown). NAC treatment significantly reduced the eosinophilia that develops in DRlyp/lyp rats during the weeks immediately preceding diabetes onset. NAC inhibits apoptosis of many cell types; however, it exerts a pro-apoptotic effect on human peripheral blood eosinophils exposed to inflammatory cytokines by inhibiting NF-κB activation and blocking production of survival factors (Martinez-Losa et al. 2007). Therefore, it is plausible that eosinophil survival was similarly modulated in the NAC-treated DRlyp/lyp rats. While a direct role for eosinophils remains to be proven in both BB rat and human T1D, recently eosinophils of T1D patients were found to express high levels of myeloid α-defensins and myeloperoxidase, suggesting that eosinophils could contribute to the innate inflammatory state that may underlie the development of diabetes (Neuwirth et al. 2012). In our study, it is possible that the modest benefit of systemic NAC administration observed here arose through action at the level of the β cell or immunomodulation or both.
Notably, in the array studies, trypsin 1 (Prss1) and chymotrypsin (Cela2a and Ctrb1) transcripts exhibited decreased abundance in the DRlyp/lyp and DR+/+ compared with the F344 and F344lyp/lyp islet pools, possibly indicating different exocrine contamination in the preparations. Islets of Langerhans consist of α cells (15–20% of total cells, produce glucagon), β cells (65–80% of total cells, produce insulin and islet amyloid polypeptide), delta cells (3–10% of total cells, producing somatostatin), pancreatic polypeptide cells (3–5% of total cells, producing pancreatic polypeptide), and epsilon cells (<1% of total cells, producing ghrelin) (Elayat et al. 1995). Although interrogated by the array, transcripts for glucagon, insulin 1, insulin 2, islet amyloid polypeptide, somatostatin, pancreatic polypeptide, and ghrelin were not differentially expressed. The concordant results of the follow-up RT-PCR studies on antioxidant enzymes (which utilized the six independent islet RNA isolations for each strain), the GST protein staining and plasma GST studies, and the NAC administration to DRlyp/lyp rats suggest that potentially modest differences in exocrine contamination of the islet preparations did not negatively influence the outcome of the study and each result independently supports that deficiencies in antioxidative defense exist in the BB rat.
While day 40 BB rat islets under-express genes related to ROS defense, they express a broader transcriptional signature consisting of both pro-survival/repair and apoptotic pathways (such as under-expression of Akt1, Pdx1, Akt2, and Perp; and overexpression of Hsf2, Ccar1, and Bclaf) suggestive of underlying dysfunction. In pursuing the hypothesis that β cell dysfunction is a prerequisite for T1D in BB rats, these day 40 analyses, before insulitis in DRlyp/lyp rats, serve as justification for longitudinal studies aimed at defining how and when this state arises and determining its mechanistic basis.
Supplementary data
This is linked to the online version of the paper at http://dx.doi.org/10.1530/JOE-12-0385.
Declaration of interest
The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.
Funding
This work was supported by National Institute of Allergy and Infectious Diseases Grants (grant numbers R01AI078713 (M J H) and P01AI42380 (Å L) and DK17047 (DERC Islet Core, The University of Washington)); Juvenile Diabetes Research Foundation International (grant number 1-2008-1026 (M J H)); and The Children's Hospital of Wisconsin Foundation.
Author contribution statement
M B, A M H, Å L, and M J H conceived, designed, interpreted results, and drafted/edited manuscript.; S K, J M F, R G, S J, M L K, S P, and Y-G C acquired data and conducted analyses and revised manuscript. All authors approve the content of this manuscript.
Acknowledgements
The authors thank Dr Paul Robertson and Dr Lydia Bryan for logistic support.
References
Awata T, Guberski DL & Like AA 1995 Genetics of the BB rat: association of autoimmune disorders (diabetes, insulitis, and thyroiditis) with lymphopenia and major histocompatibility complex class II. Endocrinology 136 5731–5735. (doi:10.1210/en.136.12.5731)
Bieg S, Koike G, Jiang J, Klaff L, Pettersson A, MacMurray AJ, Jacob HJ, Lander ES & Lernmark A 1998 Genetic isolation of iddm 1 on chromosome 4 in the biobreeding (BB) rat. Mammalian Genome 9 324–326. (doi:10.1007/s003359900759)
Bogdani M, Suenens K, Bock T, Pipeleers-Marichal M, In't Veld P & Pipeleers D 2005 Growth and functional maturation of β-cells in implants of endocrine cells purified from prenatal porcine pancreas. Diabetes 54 3387–3394. (doi:10.2337/diabetes.54.12.3387)
Borg LA, Cagliero E, Sandler S, Welsh N & Eizirik DL 1992 Interleukin-1β increases the activity of superoxide dismutase in rat pancreatic islets. Endocrinology 130 2851–2857. (doi:10.1210/en.130.5.2851)
Brahmachari S & Pahan K 2010 Myelin basic protein priming reduces the expression of Foxp3 in T cells via nitric oxide. Journal of Immunology 184 1799–1809. (doi:10.4049/jimmunol.0804394)
Cardozo AK, Heimberg H, Heremans Y, Leeman R, Kutlu B, Kruhoffer M, Orntoft T & Eizirik DL 2001 A comprehensive analysis of cytokine-induced and nuclear factor-kappa B-dependent genes in primary rat pancreatic β-cells. Journal of Biological Chemistry 276 48879–48886. (doi:10.1074/jbc.M108658200)
Chen C & Kong AN 2004 Dietary chemopreventive compounds and ARE/EpRE signaling. Free Radical Biology & Medicine 36 1505–1516. (doi:10.1016/j.freeradbiomed.2004.03.015)
Clements CM, McNally RS, Conti BJ, Mak TW & Ting JP 2006 DJ-1, a cancer- and Parkinson's disease-associated protein, stabilizes the antioxidant transcriptional master regulator Nrf2. PNAS 103 15091–15096. (doi:10.1073/pnas.0607260103)
Colle E 1990 Genetic susceptibility to the development of spontaneous insulin-dependent diabetes mellitus in the rat. Clinical Immunology and Immunopathology 57 1–9. (doi:10.1016/0090-1229(90)90017-K)
Corbett JA, Wang JL, Sweetland MA, Lancaster JR Jr & McDaniel ML 1992 Interleukin 1β induces the formation of nitric oxide by β-cells purified from rodent islets of Langerhans. Evidence for the β-cell as a source and site of action of nitric oxide. Journal of Clinical Investigation 90 2384–2391. (doi:10.1172/JCI116129)
Crisa L, Mordes JP & Rossini AA 1992 Autoimmune diabetes mellitus in the BB rat. Diabetes/Metabolism Reviews 8 4–37. (doi:10.1002/dmr.5610080104)
Delmastro MM & Piganelli JD 2011 Oxidative stress and redox modulation potential in type 1 diabetes. Clinical & Developmental Immunology 2011 593863. (doi:10.1155/2011/593863)
Eastman S, Markholst H, Wilson D & Lernmark A 1991 Leukocytosis at the onset of diabetes in crosses of inbred BB rats. Diabetes Research and Clinical Practice 12 113–123. (doi:10.1016/0168-8227(91)90088-U)
Eizirik DL, Bjorklund A & Welsh N 1993 Interleukin-1-induced expression of nitric oxide synthase in insulin-producing cells is preceded by c-fos induction and depends on gene transcription and protein synthesis. FEBS Letters 317 62–66. (doi:10.1016/0014-5793(93)81492-I)
Elayat AA, el-Naggar MM & Tahir M 1995 An immunocytochemical and morphometric study of the rat pancreatic islets. Journal of Anatomy 186 629–637.
Espiritu DJ & Mazzone T 2008 Oxidative stress regulates adipocyte apolipoprotein e and suppresses its expression in obesity. Diabetes 57 2992–2998. (doi:10.2337/db08-0592)
Fuller JM, Kwitek AE, Hawkins TJ, Moralejo DH, Lu W, Tupling TD, Macmurray AJ, Borchardt G, Hasinoff M & Lernmark A 2006 Introgression of F344 rat genomic DNA on BB rat chromosome 4 generates diabetes-resistant lymphopenic BB rats. Diabetes 55 3351–3357. (doi:10.2337/db06-0715)
Fuller JM, Bogdani M, Tupling TD, Jensen RA, Pefley R, Manavi S, Cort L, Blankenhorn EP, Mordes JP & Lernmark A et al. 2009 Genetic dissection reveals diabetes loci proximal to the gimap5 lymphopenia gene. Physiological Genomics 38 89–97. (doi:10.1152/physiolgenomics.00015.2009)
Geoffrey R, Jia S, Kwitek AE, Woodliff J, Ghosh S, Lernmark A, Wang X & Hessner MJ 2006 Evidence of a functional role for mast cells in the development of type 1 diabetes mellitus in the biobreeding rat. Journal of Immunology 177 7275–7286.
Greiner DL, Malkani S, Kanaitsuka T, Bortell R, Doukas J, Rigby M, Whalen B, Stevens LA, Moss J & Mordes JP et al. 1997 The T cell marker RT6 in a rat model of autoimmune diabetes. Advances in Experimental Medicine and Biology 419 209–216.
Groen H, Klatter FA, Brons NH, Mesander G, Nieuwenhuis P & Kampinga J 1996 Abnormal thymocyte subset distribution and differential reduction of CD4+ and CD8+T cell subsets during peripheral maturation in diabetes-prone biobreeding rats. Journal of Immunology 156 1269–1275.
Hessner MJ, Wang X, Meyer L, Geoffrey R, Jia S, Fuller J, Lernmark A & Ghosh S 2004 Involvement of eotaxin, eosinophils, and pancreatic predisposition in development of type 1 diabetes mellitus in the biobreeding rat. Journal of Immunology 173 6993–7002.
Hornum L, Romer J & Markholst H 2002 The diabetes-prone BB rat carries a frameshift mutation in Ian4, a positional candidate of Iddm1. Diabetes 51 1972–1979. (doi:10.2337/diabetes.51.6.1972)
Hosack DA, Dennis G Jr, Sherman BT, Lane HC & Lempicki RA 2003 Identifying biological themes within lists of genes with EASE. Genome Biology 4 R70. (doi:10.1186/gb-2003-4-10-r70)
Jacob HJ, Pettersson A, Wilson D, Mao Y, Lernmark A & Lander ES 1992 Genetic dissection of autoimmune type I diabetes in the BB rat. Nature Genetics 2 56–60. (doi:10.1038/ng0992-56)
Kaldunski M, Jia S, Geoffrey R, Basken J, Prosser S, Kansra S, Mordes JP, Lernmark A, Wang X & Hessner MJ 2010 Identification of a serum-induced transcriptional signature associated with type 1 diabetes in the biobreeding rat. Diabetes 59 2375–2385. (doi:10.2337/db10-0372)
Klaff LS, Koike G, Jiang J, Wang Y, Bieg S, Pettersson A, Lander E, Jacob H & Lernmark A 1999 BB rat diabetes susceptibility and body weight regulation genes colocalize on chromosome 2. Mammalian Genome 10 883–887. (doi:10.1007/s003359901108)
Kurner T, Burkart V & Kolb H 1986 Large increase of cytotoxic/suppressor T-lymphoblasts and eosinophils around manifestation of diabetes in BB rats. Diabetes Research 3 349–353.
Lenzen S 2008 Oxidative stress: the vulnerable β-cell. Biochemical Society Transactions 36 343–347. (doi:10.1042/BST0360343)
Lundsgaard D, Holm TL, Hornum L & Markholst H 2005 In vivo control of diabetogenic T-cells by regulatory CD4+CD25+ T-cells expressing Foxp3. Diabetes 54 1040–1047. (doi:10.2337/diabetes.54.4.1040)
MacMurray AJ, Moralejo DH, Kwitek AE, Rutledge EA, Van Yserloo B, Gohlke P, Speros SJ, Snyder B, Schaefer J & Bieg S et al. 2002 Lymphopenia in the BB rat model of type 1 diabetes is due to a mutation in a novel immune-associated nucleotide (Ian)-related gene. Genome Research 12 1029–1039. (doi:10.1101/gr.412702)
Malhotra D, Portales-Casamar E, Singh A, Srivastava S, Arenillas D, Happel C, Shyr C, Wakabayashi N, Kensler TW & Wasserman WW et al. 2010 Global mapping of binding sites for Nrf2 identifies novel targets in cell survival response through ChIP-Seq profiling and network analysis. Nucleic Acids Research 38 5718–5734. (doi:10.1093/nar/gkq212)
Martinez-Losa M, Cortijo J, Juan G, Ramon M, Sanz MJ & Morcillo EJ 2007 Modulatory effects of N-acetyl-l-cysteine on human eosinophil apoptosis. European Respiratory Journal 30 436–442. (doi:10.1183/09031936.00073706)
Maruta K, Treichel U, Kolb-Bachofen V & Kolb H 1989 Prospective analysis of eosinophilia in spontaneously diabetic BB rats: correlation with islet inflammation but not with diabetes development. Diabetes Research 11 173–176.
Moralejo DH, Park HA, Speros SJ, MacMurray AJ, Kwitek AE, Jacob HJ, Lander ES & Lernmark A 2003 Genetic dissection of lymphopenia from autoimmunity by introgression of mutated Ian5 gene onto the F344 rat. Journal of Autoimmunity 21 315–324. (doi:10.1016/S0896-8411(03)00138-0)
Mordes JP, Bortell R, Doukas J, Rigby M, Whalen B, Zipris D, Greiner DL & Rossini AA 1996 The BB/Wor rat and the balance hypothesis of autoimmunity. Diabetes/Metabolism Reviews 12 103–109.
Neuwirth A, Dobes J, Oujezdska J, Ballek O, Benesova M, Sumnik Z, Vcelakova J, Kolouskova S, Obermannova B & Kolar M et al. 2012 Eosinophils from patients with type 1 diabetes mellitus express high level of myeloid alpha-defensins and myeloperoxidase. Cellular Immunology 273 158–163. (doi:10.1016/j.cellimm.2011.12.001)
Ozaras R, Tahan V, Aydin S, Uzun H, Kaya S & Senturk H 2003 N-acetylcysteine attenuates alcohol-induced oxidative stress in the rat. World Journal of Gastroenterology 9 125–128.
Pfaffl MW 2001 A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Research 29 e45. (doi:10.1093/nar/29.9.e45)
Piganelli JD, Flores SC, Cruz C, Koepp J, Batinic-Haberle I, Crapo J, Day B, Kachadourian R, Young R & Bradley B et al. 2002 A metalloporphyrin-based superoxide dismutase mimic inhibits adoptive transfer of autoimmune diabetes by a diabetogenic T-cell clone. Diabetes 51 347–355. (doi:10.2337/diabetes.51.2.347)
Pisanti FA, Frascatore S & Papaccio G 1988 Superoxide dismutase activity in the BB rat: a dynamic time-course study. Life Science 43 1625–1632. (doi:10.1016/0024-3205(88)90534-6)
Poussier P, Ning T, Murphy T, Dabrowski D & Ramanathan S 2005 Impaired post-thymic development of regulatory CD4+25+ T cells contributes to diabetes pathogenesis in BB rats. Journal of Immunology 174 4081–4089.
Prasad K 2000 Oxidative stress as a mechanism of diabetes in diabetic BB prone rats: effect of secoisolariciresinol diglucoside (SDG). Molecular and Cellular Biochemistry 209 89–96. (doi:10.1023/A:1007079802459)
Reimers JI, Andersen HU, Mauricio D, Pociot F, Karlsen AE, Petersen JS, Mandrup-Poulsen T & Nerup J 1996 Strain-dependent differences in sensitivity of rat β-cells to interleukin 1β in vitro and in vivo: association with islet nitric oxide synthesis. Diabetes 45 771–778. (doi:10.2337/diabetes.45.6.771)
Rieneck K, Bovin LF, Josefsen K, Buschard K, Svenson M & Bendtzen K 2000 Massive parallel gene expression profiling of RINm5F pancreatic islet β-cells stimulated with interleukin-1β. APMIS: Acta Pathologica, Microbiologica, et Immunologica Scandinavica 108 855–872. (doi:10.1111/j.1600-0463.2000.tb00009.x)
Rocic B, Vucic M, Knezevic-Cuca J, Radica A, Pavlic-Renar I, Profozic V & Metelko Z 1997 Total plasma antioxidants in first-degree relatives of patients with insulin-dependent diabetes. Experimental and Clinical Endocrinology & Diabetes 105 213–217. (doi:10.1055/s-0029-1211754)
Scott J 1990 The spontaneously diabetic BB rat: sites of the defects leading to autoimmunity and diabetes mellitus. A review. Current Topics in Microbiology and Immunology 156 1–14.
Sigfrid LA, Cunningham JM, Beeharry N, Hakan Borg LA, Rosales Hernandez AL, Carlsson C, Bone AJ & Green IC 2004 Antioxidant enzyme activity and mRNA expression in the islets of Langerhans from the BB/S rat model of type 1 diabetes and an insulin-producing cell line. Journal of Molecular Medicine 82 325–335. (doi:10.1007/s00109-004-0533-4)
Sturn A, Quackenbush J & Trajanoski Z 2002 Genesis: cluster analysis of microarray data. Bioinformatics 18 207–208. (doi:10.1093/bioinformatics/18.1.207)
Sweet IR, Cook DL, DeJulio E, Wallen AR, Khalil G, Callis J & Reems J 2004 Regulation of ATP/ADP in pancreatic islets. Diabetes 53 401–409. (doi:10.2337/diabetes.53.2.401)
Szabo C, Mabley JG, Moeller SM, Shimanovich R, Pacher P, Virag L, Soriano FG, Van Duzer JH, Williams W & Salzman AL et al. 2002 Part I: pathogenetic role of peroxynitrite in the development of diabetes and diabetic vascular complications: studies with FP15, a novel potent peroxynitrite decomposition catalyst. Molecular Medicine 8 571–580.
Tanaka Y, Gleason CE, Tran PO, Harmon JS & Robertson RP 1999 Prevention of glucose toxicity in HIT-T15 cells and Zucker diabetic fatty rats by antioxidants. PNAS 96 10857–10862. (doi:10.1073/pnas.96.19.10857)
Tusher VG, Tibshirani R & Chu G 2001 Significance analysis of microarrays applied to the ionizing radiation response. PNAS 98 5116–5121. (doi:10.1073/pnas.091062498)
Wallis RH, Wang K, Marandi L, Hsieh E, Ning T, Chao GY, Sarmiento J, Paterson AD & Poussier P 2009 Type 1 diabetes in the BB rat: a polygenic disease. Diabetes 58 1007–1017. (doi:10.2337/db08-1215)
Wang X, Jia S, Geoffrey R, Alemzadeh R, Ghosh S & Hessner MJ 2008 Identification of a molecular signature in human type 1 diabetes mellitus using serum and functional genomics. Journal of Immunology 180 1929–1937.
Welsh N, Margulis B, Borg LA, Wiklund HJ, Saldeen J, Flodstrom M, Mello MA, Andersson A, Pipeleers DG & Hellerstrom C et al. 1995 Differences in the expression of heat-shock proteins and antioxidant enzymes between human and rodent pancreatic islets: implications for the pathogenesis of insulin-dependent diabetes mellitus. Molecular Medicine 1 806–820.
Zipris D, Hillebrands JL, Welsh RM, Rozing J, Xie JX, Mordes JP, Greiner DL & Rossini AA 2003 Infections that induce autoimmune diabetes in BBDR rats modulate CD4+CD25+T cell populations. Journal of Immunology 170 3592–3602.
(MJ Hessner is now at Division of Endocrinology, Department of Pediatrics, Medical College of Wisconsin, 9000 West Wisconsin Avenue, Milwaukee, Wisconsin 53201, USA)