Csf2 and Ptgs2 Epigenetic Dysregulation in Diabetes-prone Bicongenic B6.NODC11bxC1tb Mice

Mouse strains

Ethics statement: All procedures were conducted according to Institutional Animal Care and Use Committee (IACUC) approved protocols at the University of Florida (UF), Kennedy Space Center Space Life Sciences Laboratory, and Sanford-Burnham Medical Research Institute (SBMRI) at Lake Nona, Florida. This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All efforts were made to minimize suffering and provide humane treatment for the animals used in this study.

Three- to 20-week-old NOD, C57BL/6, and C57L female and male mice (The Jackson Laboratory, Bar Harbor, ME, USA) were used in these studies. At least two mice of each strain were used for tissues per analysis, and each analysis was run at minimum in triplicate. These strains were maintained as breeding stock in the University of Florida College of Medicine Pathology Department Specific Pathogen Free (SPF) colony, in microisolator cages with food and water ad libitum. Figure 1 gives a schematic representation of the regions on chromosomes contributed by NOD, C57BL/6, and C57L parental strains in the NOD.LC11 ¹⁷ and B6.NODC11^16,18 congenic strains used in this study. Subcongenic strains of the B6.NODC11 and NOD.LC11 mice (labeled with letters after C11 in their names, Fig. 1) were bred by backcross to C57BL/6 then intercross as in speed congenic protocols previously described.^18,20

OPEN IN VIEWER

Figure 1.

Summary of congenic mouse chromosomal maps for chromosomes 11, 1, and 17. Color indicates the genetic origin of the sequences on Chr 11 (A), Chr 1 (B), and Chr 17 (C) and in the multi-congenic strains (D): black = C57BL/6, grey = C57L, and white = NOD. The location of Idd and other regions/genes of interest are listed on the left side of the maps and highlighted across the maps by the shaded horizontal areas. Strain names and shorthand identifiers (in parentheses) are listed at the top of the maps. (E) Summary of the current phenotypic findings for each strain including data from the current study as well as our previously published work ¹³ , ¹⁵ , ¹⁶ , ¹⁹ and that of others. ¹⁷ , ¹⁸ Histological phenotype scoring by three blinded reviewers are indicated as Peri = infiltration around pancreatic inlets; invasive = destructive islet infiltration;‘–’ = no infiltrate seen.

NOD background subcongenic strain NOD.LC11 (NOD.L-(D11Mit314-D11Mit42)/McdfJ; now available from Jackson Laboratory, stock #008053) and its subcongenic recombinant derivatives NOD.LC11b (NOD.L-(D11Mit314-Csf2)/McdfJ; Jackson Lab stock #008055) and NOD.LC11e (NOD.L-(Csf2-D11Mit339)/McdfJ; Jackson Lab stock #008056) were developed at the University of Virginia. ¹⁷ NOD. LC11e mice carry a C57L interval covering the Csf2 gene for GM-CSF and Idd4 interval, but excluding the promoter region (Fig. 2). NOD.LC11b contains a C57L region, which includes the Idd4.3 region down through the Csf2 gene but excludes all the regions telomeric and downstream of the Csf2 gene.

OPEN IN VIEWER

Figure 2.

Sequence analysis of Csf2 promoter region and definition of the STAT5 binding site polymorphisms involved in NOD myeloid cell phenotypes and chromosome 11 diabetes susceptibility. Genomic DNA from each congenic and parental strain was prepared from liver and amplified in PCR using primers designed to amplify the –969 bp to –3 bp sequence upstream of the Csf2 gene. ⁶ , ⁷ Boxed sequences indicate potential STAT protein consensus binding sequences (CTTNNN(N)GAA) and microsatellite polymorphism seen in NOD compared to the GT repeat seen in control strains C57BL/6 and C57L. Dashed areas indicate gaps in sequence of some strains (ie, base pairs not found in some strains) indicated by CLUSTAL alignment. Vertical lines indicate sites of transition between NOD and control strain DNA sequences in congenic strains. N = NOD genetic sequence; B = control strain (C57BL/6 for B6.NODC11b strain; C57L for NOD.LC11E and NOD.LC11B strains). Asterisks (*) indicate sequence homology across all strains.

The B6.NODC11b and B6.NODC1tb subcongenic mouse strains were derived for this study from backcrossing the B6.NODC11 and B6.NODC1 mice in the University of Florida Department of Pathology SPF Mouse colony, which were originally derived by Yui and Wakeland, ¹ with C57BL/6J mice from Jackson Laboratory for three to five generations then intercrossing the offspring for three generations to homozygosity. Each subcongenic line was intercrossed for at least three more generations to establish the strains depicted in Figure 1. These subcongenic mice do not develop diabetes but have periductal insulitis (data not shown).

The B6.NOD congenic mice, along with NOD and C57BL/6 mice purchased from Jackson Laboratories, were originally bred and housed in the University of Florida SPF colony. The strains used in these studies were maintained at the University of Florida, Kennedy Space Center, and later at the new SBMRI at Lake Nona SPF vivarium facilities.

Diabetes incidence in the NOD.LC11 and B6.NODC11 and B6.NODC1 strains was measured by blood and urine glycemia and confirmed by pathological review of immune infiltration and beta cell destruction in pancreatic tissues as previously described.^16–18 Genotype analysis on tail or ear biopsy tissue was routinely performed at weaning of subcongenic litters to confirm maintenance of NOD genetic intervals and C57BL/6 background genetics at all other Idd loci were done using microsatellite polymorphism polymerase chain reaction (PCR) analysis as described by Yui et al ¹⁸ and McDufe. ¹⁷

Development of multicongenic strains B6.NODC11 bxC1tb and B6.NODC11bxC1tbxC17

Subcongenic line mice, B6.NODC11b, and B6.NODC1tb mice were cross-bred for two to three generations and the offspring that were homozygotic NOD at each allele were bred to their siblings for five more generations to establish the bicongenic line, B6.NODC11bxC1tb (sperm available through Jackson Laboratory). The line was then bred for at least five more generations to stabilize the strain and has been maintained by homozygotic offspring breeding for over 15 generations, with breeder changes at least every 6 months. B6.NODC11bxC1tb mice were genotyped at all Idd intervals on chromosomes 1–19 originally described for their single congenic ancestors ¹⁸ and found to be only carrying NOD genotype alleles at the intervals depicted in maps of Figure 1 and sequence in Figure 2, including two possible STAT CCTNNN(N)GAA consensus binding sites (one found in NOD, the other in C57BL/6). This strain bred well and was exceptionally good at nurturing of litters, including being excellent at fostering other strain pups. There was a low incidence of stillbirth (3 litters over 15 generations), congenial cranial bone deformity in later generation (3 mice over 15 generations), sterile hermaphrodites (2/15 generations), and even lower incidence of cataract development (1 mouse/15th generation). Invasive, destructive insulitis scored by three blinded pathologist readers for this strain was 3/38 females (8%) and 5/30 males (16%), starting at 9 weeks of age through to 30 weeks of age. Hyperglycemia was seen in 2/18 females (11%) and 7/23 males (30%).

B6.NODC17 mice maintained for over 20 generations in the University of Florida Department of Pathology SPF Mouse colony, which were originally derived by Yui and Wakeland, ¹⁸ were bred to B6.C11bxC1tb mice for three to five generations and Chr11/Chr1/Chr17 multi-congenic line mice were established out of heterozygotic breeding of sibling pairs for two to three generations and then maintained by breeding homozygotic offspring in sibling mating for an additional five generations. The line has been maintained for over 12 generations to create the triple congenic strain, B6.NODC11bxC1tbxC17.

This strain bred poorly and was very poor at nurturing of litters, especially prone to male parent infanticide and female failure to lactate. There was a high incidence of stillbirth and miscarriages (18 litters over 12 generations). Invasive, destructive insulitis scored by three blinded pathologist readers for this strain was 4/18 females (22%) and 6/24 males (25%), starting at 9 weeks of age through to 30 weeks of age. Hyperglycemia was seen in 5/21 females (24%) and 18/34 males (53%).

Tissue collection, cell culture, and sample preparation

Subcongenic and bicongenic strain mice were analyzed at weaning (3–4 weeks) through to 30 weeks of age when retired from breeding. Blood monocytes were collected from 3- to 4-week-old pups at the time of tail or ear tissue biopsy, from tail or cheek vein collection on live adult mice or from cardiac puncture postmortem after euthanasia. Mice were euthanized by overanesthetizing and cervical dislocation. After blood collection, the peritoneal cavity was filled with ice-cold RPMI medium supplemented with 10% fetal calf serum and 1% antibiotic/antimycotic mix (Cellgro-Mediatech, Herndon, VA, USA) using a 20-gauge needle and syringe. The lavage fluid and cells were withdrawn and washed with cold media by centrifugation. Liver tissue was collected after lavage for use in genomic DNA preparation, along with salivary gland and pancreatic tissues for histology. After lavage, long bones from the hind limbs were then collected and bone marrow flushed out of the bones using a 30-gauge needle and syringe filled with ice-cold RPMI medium supplemented with 10% fetal calf serum and 1% antibiotic/antimycotic mix (Cellgro-Mediatech).

The bone marrow and lavage cells were washed with cold media and then the red blood cells in samples were lysed by incubation in sterile cold 0.84% NH₄Cl buffer. The remaining cells were then plated on tissue culture dishes and fed with fresh sterile medium for culture at 37°C/5% CO2. After 1 hour of incubation, lavage cell cultures were washed clear of nonadherent cells and the remaining adherent macrophages were refed with medium. The bone marrow and macrophage cultures were maintained in medium alone or with 1000 U/mL of GM-CSF (Biosource/Invitrogen, Carlsbad, CA, USA) or with 2 mg/mL anti-GM-CSF–blocking antibodies (Endogen/ Pierce, Rockville, IL, USA) for 15 minutes to 24 hours at 37°C/5% CO2. After the media harvest, the cells were fixed in situ with 1% (v/v C_f) formaldehyde (methanol-free, Sigma-Aldrich, St Louis, MO, USA) added in the remaining media for 10 minutes at 37°C, then washed with 1 × phosphate-buffered saline (PBS) and sonicated in sodium dodecyl sulfate (SDS) lysis buffer (Upstate/Millipore, Chicago, IL, USA) + protease inhibitors (Roche, Indianapolis, IN, USA) for later analysis in chromatin immunoprecipitation (ChIP).^19,21

Pancreatic histology was confirmed in tissues collected at necropsy and processed through hematoxylin and eosin staining by either the University of Florida Pathology Core or the SBMRI Histology Core. Slides were examined by three blinded readers at the University of Florida and compared with the coded identification by the authors kept at SBMRI.

Phosphotyrosine STAT5 flow cytometric analysis

Aliquots of macrophage and peripheral blood cells were taken to confirm phenotypic identification and phosphotyrosine STAT5 analysis by flow cytometry as previously described.^16,19

GM-CSF and PGE2 production analysis

After incubation, half of media volume from bone marrow and peritoneal macrophage cultures was collected and frozen at – 80°C to measure GM-CSF concentration by Luminex (Upstate Biotech Beadlyte kits, Upstate/Millipore) and enzyme-linked immuosorbent assay (ELISA; BD Biosciences OptiEIA kits, San Diego, CA, USA) as well as PGE2 by ELISA (GE/ Amersham, Piscataway, NJ, USA) as previously described.^13,15

PGS2/COX2 western blot analysis

Peritoneal macrophage cultures were also extracted for protein analysis to detect PGS2/COX2. Cultures of 5 × 10⁶ lavage cells were incubated for 1 hour at 37°C/5% CO2 to allow for macrophage adhesion. Nonadherent cells were washed off with sterile 1 × PBS and the adherent cells extracted by sonication in SDS lysis buffer (Upstate/Millipore) + protease inhibitors (Roche). Extracts were frozen at – 80°C until analysis. Thawed extracts were resuspended with 5 × Leammli sample buffer (Tris-HCL [pH 6.5], 50% glycerol, 1% SDS, 1% 2-mercaptoethanol, 0.1% bromophenol blue, made with reagents from Fisher Scientific, Pittsburgh, PA, USA) and heated to 100°C for 3 minutes. Fifty microliters of each sample was run on 4%–20% acrylamide gels and transferred to Odyssey polyvinylidene floride (PVDF) membranes for western blot analysis using anti-COX2 polyclonal antibodies from Cayman Chemicals (Ann Arbor, MI, USA) and detection using the Li-Cor anti-Rabbit IR700 dye conjugates fluorescence at 700 nm on a Li-Cor Odyssey IR Imaging system (Li-Cor Biosciences, Lincoln, NE, USA).

Sequence analysis of Csf2 gene promoter and Ptgs2 enhancer region

Genomic DNA from each strain was prepared from liver of each strain and amplified in PCR using Roche Biosciences Master Mix reagents and primers (3′ CTA AAG CAT GTT TCT TGG CTA; 5′ AAA TAA GGT CCA GCC CAA TG; Integrated DNA Technologies, Coralville, IA, USA) designed to amplify the –3 to –969 bp sequence upstream of the Csf2 gene.^6,7 The amplified DNA was gel purified using Qiagen gel extraction reagents and phenol/chloroform extracted (Qiagen, Valencia, CA, USA). The amplified fragment was then used as template in a Big Dye PCR amplification reaction (Applied Biosystems, Foster City, CA, USA) and sequenced using an AB capillary sequence analyzer (Applied Biosystems). ChromasLite and ClustalW freeware were used for the sequence analysis and alignment. Overlapping sequence region is depicted in Figure 2. Sequence analysis of the region between Il10 and Ptgs2 genes on mouse Chr 1 did not reveal any sequence polymorphisms unique to the NOD mouse or its congenics as compared with C57BL/6 and C57L mice other than the already reported length polymorphisms between these strains detectable by Massachusetts Institute of Technology (MIT) marker analysis.^17,18

ChIP and ChIP multiplex format analysis

Four to five million cells from bone marrow cultures or ex vivo peritoneal macrophages were fixed in situ with 1% formaldehyde in 1 × PBS (methanol-free, Sigma-Aldrich) for 10 minutes at 37°C prior to lysis with SDS lysis buffer (Upstate/Millipore). The fixed lysates were sonicated to disrupt membranes and shear chromatin to approximately 1,000-bp fragments then frozen until analysis. Once thawed, the samples were divided into aliquots for each run of the analysis. For single-ChIP analyses, aliquots used for IP were precleared with salmon sperm DNA Protein A or G agarose beads (Upstate/Millipore), then incubated overnight at 4°C with antityrosine phosphorylated STAT5 (STAT5Ptyr) antibodies (Upstate/Millipore).

For multiplex format analyses, macrophage extracts were thawed and precleared of nonspecific DNA and Protein G binding proteins by incubation in situ with salmon sperm-Protein G–coated sepharose beads (UpState/Millipore) and bead-free supernatants aliquoted equally over eight PCR microtubes in a strip/capped format (Eppendorf, Westbury NY, USA). The aliquots were diluted with ChIP dilution buffer (16.7 mM Tris–HCl [pH 8.1], 0.01% SDS, 1.1% Triton X-100, 1.2 mM ethylene diamine tetraacetate (EDTA), 167 mM NaCl; Sigma-Aldrich and Fisher Scientific reagents) to a total volume of 200 μL each. Salmon sperm-Protein G sepharose beads (Upstate/ Millipore) were preincubated with 1 μg/30 μL beads suspension in a 96-tube Multiplex format using parallel ChIP isolations with one tube of each of the following Upstate/Millipore ChIP–certified antibodies: anti-STAT5Ptyr, anti-RNA polymerase II, antiacetylated histone H3, anti-P300/anti-CBP acetylase, and anti-SMRTe deacetylase. Thirty microliters of the antibody-bead mixes was added to six of the seven tubes of sample, with the final tube being left unprecipitated to act as a total DNA control.

After incubation, the antibody-bound chromatin complexes were precipitated using salmon sperm DNA Protein G agarose beads and washed extensively with a series of increasing stringency buffers (low salt, high salt, LiCl, Tris-EDTA (TE); ChIP reagent kit, Upstate/Millipore). A nonspecific antibody control (mouse IgG, Upstate/Millipore) and a no extract sham IP were run as negative controls in both single and multiplex format ChIP analyses. Total cell and ChIP extract aliquots were dissociated from the beads in 1% SDS, 0.1 M bicarbonate buffer (Fisher Scientific). NaCl was then added to a final concentration of 500 mM and the samples incubated for 4 hours at 65°C to reverse the formaldehyde crosslinks. DNA was purified from these aliquots for PCR and real-time PCR amplification of DNA sequences from Csf2 promoters, which have been shown to be epigenetic regulatory sites for inducible Csf2 expression,^6,7 and from the Ptgs2 enhancer STAT5 binding site. ²

Real-time PCR was run using Sybr Green Master Mixes (Applied Biosystems or BioRad, Hercules, CA, USA) with the same hot-start protocol designed to remove secondary structure in the DNA template. A cycle profile of 98°C 5 minutes, 94°C 30 sec, 55°C–60°C (dependent on the primer set used) 30 sec, 72°C 30 sec, for 45 cycles, was used.

Real-time amplification quantitation was compared on the basis of the R-value calculated as R = 2^{((Nonspecific Ig ChIP ΔCt)–-(anti–STAT5Ptyr ChIP ΔCt))}. ²¹ Statistical analyses of data were performed using Prism 5/6 software (GraphPad, San Diego, CA, USA).

Sequence analysis of Csf2 promoter region defines STAT5 binding site polymorphisms

Sequence analysis of the Csf2 promoter region was done in all congenic strains and compared in a ClustalW sequence alignment analysis (Figs. 1 and 2). B6.NODC11b and the NOD.LC11e mice share an overlapping NOD genetic region in the Csf2 promoter that does not extend into the Csf2 coding region (Fig. 2).

NOD STAT5Ptyr and GM-CSF phenotypes segregate with the NOD Csf2 promoter

NOD.LC11 and NOD. LC11b mice develop a strong peri-islet insulitis, which remains noninvasive despite having autoreactive T cells present in their T cell population; however, the NOD.LC11e mice developed invasive and destructive insulitis that resulted in overt diabetes. ¹⁷ In contrast, pancreatic tissue analyzed from B6.NODC11b and B6.NODC1tb mice shows very little infiltration, mainly limited to periductal areas of the pancreas near islets (data not shown).

Both strains exhibit increased STAT5 phosphorylation (Fig. 3A and C) and GM-CSF expression (Fig. 3B and D), indicating that the overlap region delineates the genetic locus most likely responsible for the GM-CSF and STAT5 phenotypes seen in the NOD. This shared B6.NODC11b and NOD.LC11e region contains at least one NOD unique potential STAT5 GAS binding (TTCNNNGAA) and two altered STAT6 GAS binding sites (TTCNNNNGAA) ²² and includes the regions previously defined as epigenetic regulatory sites.^6,7 Of note, the levels of STAT5Ptyr in Chr11 congenic strains on both NOD and nonautoimmune genetic backgrounds were significantly higher than in nonautoimmune strains but did not reach the levels seen in the NOD (Fig. 3A and C); whereas, the GM-CSF production of these strains were comparable to that seen in NOD myeloid cells (Fig. 3B and D). These data suggest that other NOD genetic components may be influencing the activation and/or persistence of STAT5Ptyr.

OPEN IN VIEWER

Figure 3.

STAT5Ptyr GM-CSF and PGE2 phenotypes for B6.NOD and NOD. L congenic mouse bone marrow cells and peritoneal macrophages. Bone marrow cells (A and C) and peritoneal macrophages (B and D) were isolated from female mice of each strain immediately processed for STAT5Ptyr flow cytometry (A and B) or cultured for 24 hours at 37°C/5% CO₂ without stimulation. Cell-free supernatants from these unstimulated cultures were then analyzed by Luminex or ELISA for GM-CSF (B and D). Both bone marrow and macrophages from NOD and strains carrying the Csf2 promoter sequences from the NOD ancestor (NOD, NOD.LC11e, and B6.NODC11b) exhibited elevated levels of STAT5Ptyr and GM-CSF in their CD11b+ myeloid cells compared with nonautoimmune strains or strains that had this sequence from non-NOD ancestors, C57BL/6 or C57L, and NOD.LC11b mice. Means were compared pairwise with NOD and with C57BL/6 using Mann–Whitney U-tests or with all strains via one-way analysis of variance. STAT5Ptyr data are representative of multiple runs of the analysis with a minimum of 2 and a maximum of 66 samples per strain tested. GM-CSF Luminex/ELISA data are representative of multiple runs of the analysis with a minimum of 2 and a maximum of 25 samples per strain tested.

STAT5 chromatin binding at the Csf2 promoter in myeloid cells was also defined by the presence of NOD sequence in this region (Fig. 4). Untreated ex vivo macrophages of NOD mice had enhanced STAT5Ptyr chromatin binding (>56,000-fold above the C57BL/6 control strain) at the site as did macrophages of the B6.NODC11b strain (30-fold higher than the C57BL/6 strain) (Fig. 4A). In contrast, NOD.LC11e congenic mouse macrophages did not have STAT5Ptyr binding at the Csf2 site. These data indicate that NOD genetic sequence in the promoter of Csf2 not its coding sequence influence the ability of STAT5Ptyr to bind to this region and potential influence GM-CSF expression and activation of STAT5.

OPEN IN VIEWER

Figure 4.

B6.NOD and NOD. L congenic mice STAT5 binding at the Csf2 promoter and Ptgs2 enhancer regions. Peritoneal macrophages were isolated from female mice of each strain and cultured for 24 hours at 37°C/5% CO₂ without stimulation. Cultured cells were then fixed and extracted for ChIP analysis using anti-STAT5Ptyr antibodies to isolate STAT5-associated chromatin amplified in real-time PCRs using primers specific for (A) the Chr 11 Csf2 promoter region^6,7,19 and (B) the Il10-Ptgs2 Chr1 region.^2,18 Bars depict mean log R-values with standard error of the mean error bars. Data are representative of at least three runs of the analysis. (C) Cell-free supernatants were collected for PGE2 analysis by ELISA. Data are representative of at least two runs of the analysis. Bars depict mean with standard error of the mean error bars. Due to unequal variance in the subgroups, no accurate statistical testing was possible. Fold increases for each parameter and each strain as compared to the C57BL/6 nonautoimmune control strain are indicated on the graphs.

NOD Csf2 promoter polymorphisms affects STAT5 binding at the Ptgs2 enhancer

The high GM-CSF production of NOD and T1D human myeloid cells can enhance the PGS2/COX2 production of these cells, even without additional activation stimuli, such as by lipopolysaccharide (endotoxin) (LPS).^13,15 ChIP analysis of macrophage STAT5Ptyr binding at an enhancer site upstream of the Ptgs2 gene indicated that NOD sequence at the Chr11 Csf2 promoter affected STAT5Ptyr binding at the Chr1 interval between Il10-Ptgs2 even in B6.NODC11b mice, which do not have NOD sequence on Chr1 (>38,000-fold higher than seen in the C57BL/6 strain, Fig. 4B). However, mice that had NOD sequence at the Chr 1 Il10-Ptgs2 region without NOD sequence at the Csf2 promoter (ie, NOD.LC11b, NOD.LC11e, and B6.NODCt1b) also showed STAT5Ptyr binding above the level seen in C57BL/6J and C57LJ control macrophages (Fig. 4B, data not shown). These findings give evidence that both sites can influence the ability of STAT5Ptyr to bind at the Chr1 interval.

PGE2 production was measured as an indicator of PGS/ COX expression. The production of this prostaglandin in unstimulated peritoneal macrophages mirrored the STAT5 binding at the Ptgs2 enhancer (Fig. 4C). PGE2 production was increased in macrophages of mice that had NOD genetic sequence at the Ptgs2 site on Chr1. Thus, GM-CSF–activated STAT5 binding at the Csf2 promoter may play a role in proinflammatory dysregulation of the Ptgs2 gene promoting aberrant expression of PGS2/COX2 in NOD macrophages, other factors besides GM-CSF are influencing PGS2/COX2 expression.

Bicongenic B6.NODC11bxC1tb has STAT5-mediated interaction between Chr 11 Ud4.3 and Chr 1 Idd5 with aberrant PGS2/COX2 expression and autoimmune diabetes pathology

To directly test the epigenetic interaction of the Idd4.3 Chr 11 locus with the Chr 1 locus in the activation of aberrant PGS2/COX2 expression, we cross-bred the recombinant B6.NODC11b strain mice with B6.NODC1tb mice to create a bicongenic strain, B6.NODC11bxC1tb. This strain has two recombination-shortened NOD Idd loci (130.8–149.7 Mb, of Idd5 on Chr 1 and at 53.85 Mb of Idd4.3 on Chr11) on C57BL/6 genetic background. These bicongenic mice have NOD defined regulatory regions that are dysregulated through cytokine-induced gene activation. Our analyses of these mice confirmed that the STAT5-mediated interaction between these two loci promoted persistent STAT5Ptyr- (Fig. 5A) and STAT5Ptyr-mediated binding epigenetic phenotypes at these loci (Fig. 6), as seen in the NOD ancestor. In addition, the B6.NODC11bxC1tb mice with only these two NOD-derived loci developed an invasive myeloid cell infiltration, hyperglycemia, and apparent immunopathological islet beta cell destruction in approximately 12% to 22% overall (8%–12% [females] and 16%–30% [males]) of the animals by age 30 weeks (Fig. 5C–G).

OPEN IN VIEWER

Figure 5.

Phenotypic analyses of multicongenic B6.NODC11bxC1tb and B6.NODC11bxC1tbxC17 mice. (A) STAT5Ptyr activation in macrophages, (B) GM-CSF production by myeloid bone marrow cells (left) and macrophages (right), and (C) blood glucose levels of NOD, C56BL/6, and congenic mice. Survival curves of diabetes incidence in NOD, B6.NODC11bxC1tb, and B6.NODC11bxC1tbxC17 mice based on (D) blood glucose. (E) Invasive beta cell destructive insulitis. (F) Invasive insulitis in a B6.NODC11bxC1tb male mouse shown at 20 × magnification. (G) Macrophage marker F4/80 immunohistochemical staining in islets of B6.NODC11bxC1tb mice shown at 20 × (upper) and 100 × (lower) magnification. (H) Western blot analysis of PGS2/ COX2 Expression in 5 × 10⁶ unactivated NOD, C57BL/6, B6.NODC11bxC1tb, and B6.NODC11bxC1tbxC17 mouse peritoneal macrophages. (I) ELISA analysis of macrophage PGE2 production of NOD, C57BL/6, and congenic mice. Note that multiple isoforms of PGS2/COX2 appear to be present with different isoforms predominant in the different strains (eg, NOD and C57BL/6 have a longer isoform predominant compared to the congenic strains). STAT5Ptyr, GM-CSF, and PGE2 data are representative of at least three runs of each assay and survival/diabetic phenotype testing are the results from continuous monitoring of 20–30 animals per stain over 6 months in a longitudinal study. Bars (A, B, and H) and horizontal line (C) indicate means with standard error of the mean error bars. P values for pairwise and one-way analyses of variance are depicted above the data in graphs (A, B, C, and I).

OPEN IN VIEWER

Figure 6.

Multicongenic mice multiplex ChIP analysis of STAT5 binding at the Csf2 promoter and Ptgs2 enhancer regions. (A) Multiplex ChIP analyses for protein association with the Chr11 Csf2 promoter in NOD, C57BL/6, B6.NODC11bxC1tb, and B6.NODC11bxC1tbxC17 mouse unstimulated peritoneal macrophages cultured for 2 hours at 37°C/5% CO₂. (B) Multiplex ChIP analyses done as in (A) for epigenetic factor binding/action at Il10-Ptgs2 Chr1 locus containing enhancer elements for Ptgs2. Parallel ChIP analyses in the multiplex were precipitated with antibodies against STAT5Ptyr (STAT5P), RNA polymerase II (POLII), acetylated histone H3 (ACETYL H3), CBP/P300 histone acetylase (CBP/P300), and SMRT/NCoR histone deacetylase (SMRT). Bar pattern key indicates protein associated with the chromatin; bars indicate mean log R-value of real-time PCR-specific amplification with standard error of the mean error bars. Mann-Whitney U-tests were performed between all pairs of data, indicating that no statistically significant differences were found in any of the data except in CBP/P300 and SMRT bindings indicated: *P = 0.0519 between B6.NODC11bxC1tb and C57BL/6 CBP/ P300 and SMRT binding; **P = 0.0293 between B6.NODC11bxC1tbxC17 and C57BL/6 SMRT binding.

Unlike their B6.NODC11b congenic progenitors, the GM-CSF production in the B6.NODC11bxC1tb mice was not as high as in the NOD (Fig. 5B). Unstimulated macrophage PGS2/COX2 expression was seen in NOD, bicongenic, and tri-congenic strains, but not in the nondiabetic control strains (Fig. 5H). PGE2 production by B6.NODC11bxC1tb macrophages were higher than that seen in the NOD (Fig. 5I). These data suggest that a specific level of GM-CSF expression and/or other genes (eg, Irf-1, Il4, Il3) also regulated by the regional epigenetic control in the Chr 11 Idd4.3 or in the Chr 1 Idd5 regions (eg, Il10) may be affected by STAT5Ptyr binding to influence Ptgs2 expression.

In NOD and B6.NODC11bxC1tb macrophages, STAT5Ptyr binding at both the Csf2 promoter and Ptgs2 enhancer was found in association with CBP/P300 histone acetylase, acetylated histone H3, and RNA polymerase II at these same sites (Fig. 6). In B6.NODC11bxC1tb macrophage, deacetylase SMRT/NCoR was detected in association with STAT5Ptyr binding at the Csf2 promoter but not at the Ptgs2 enhancer. In contrast, C57BL/6 macrophages had markedly lower binding of all but SMRT/NCoR deacetylase at the Ptgs2 site. Overall binding differences were not statistically significant with the exception of SMRT/NCoR binding being higher in B6.NODC11bxC1tbxC17 cells than in the C57BL/6 cells (P = 0.0293, Mann–Whitney U-test, 95% confidence level). Both CBP/P300 and SMRT/NCoR binding in the B6.NODC11bxC1tb appeared modestly higher than seen in the C57BL/6 (P = 0.0519, Mann–Whitney U-test, 90% confidence level).

Congenic addition of Idd1 to B6.NODC11xC1tb background enhances penetrance but is not required for autoimmune diabetes pathology

The unique MHC of the NOD is not required for the immunopathology seen in the B6.NODC11bxC1tb mice. Congenic breeding addition of the NOD MHC locus (Idd1 Chr17; B6.NODC11bxC1tbxC17 mice, Fig. 1) did increase the incidence of diabetes to overall 23% to 42% (25 to 53% of males and 22 to 24% of females), as measured by hyperglycemia, invasive immune cell infiltration, and beta cell loss in the pancreas (Figs. 1, 5, and 6). The addition of the Idd1 locus to the B6.NODC11bxClt1b background did enhance the STAT5Ptyr activation in monocytes (Fig. 5A) but decreased GM-CSF production (Fig. 5B) and STAT5Ptyr binding at the Csf2 promoter and the Ptgs2 enhancer compared to B6.C11bxC1tb macrophages (Fig. 6). B6.NODC11bxC1tbXC17 had lower STAT5Ptyr binding at both Csf2 and Ptgs2 sites than either NOD or B6.NODC11bxC1tb progenitors, but the presence of CBP/ P300, acetylated histone H3, and RNA polymerase II in B6.NODC11bxC1tbXC17 macrophages at both Csf2 and Ptgs2 regulatory sites suggest these cells still have the potential for epigenetic and transcriptional activity in these regions. The binding of SMRT/NCoR deacetylase binding at these sites is enhanced in the NOD, bicongenic, and tri-congenic macrophages as well, suggesting that deacetylation regulation is also affected by the genetic changes seen in these diabetic-prone strains.