Early over expression of messenger RNA for multiple genes, including insulin, in the Pancreatic Lymph Nodes of NOD mice is associated with Islet Autoimmunity
© Regnault et al; licensee BioMed Central Ltd. 2009
Received: 9 May 2009
Accepted: 2 October 2009
Published: 2 October 2009
Autoimmune diabetes (T1D) onset is preceded by a long inflammatory process directed against the insulin-secreting β cells of the pancreas. Deciphering the early autoimmune mechanisms represents a challenge due to the absence of clinical signs at early disease stages. The aim of this study was to identify genes implicated in the early steps of the autoimmune process, prior to inflammation, in T1D. We have previously established that insulin autoantibodies (E-IAA) predict early diabetes onset delineating an early phenotypic check point (window 1) in disease pathogenesis. We used this sub-phenotype and applied differential gene expression analysis in the pancreatic lymph nodes (PLN) of 5 weeks old Non Obese Diabetic (NOD) mice differing solely upon the presence or absence of E-IAA. Analysis of gene expression profiles has the potential to provide a global understanding of the disease and to generate novel hypothesis concerning the initiation of the autoimmune process.
Animals have been screened weekly for the presence of E-IAA between 3 and 5 weeks of age. E-IAA positive or negative NOD mice at least twice were selected and RNAs isolated from the PLN were used for microarray analysis. Comparison of transcriptional profiles between positive and negative animals and functional annotations of the resulting differentially expressed genes, using software together with manual literature data mining, have been performed.
The expression of 165 genes was modulated between E-IAA positive and negative PLN. In particular, genes coding for insulin and for proteins known to be implicated in tissue remodelling and Th1 immunity have been found to be highly differentially expressed. Forty one genes showed over 5 fold differences between the two sets of samples and 30 code for extracellular proteins. This class of proteins represents potential diagnostic markers and drug targets for T1D.
Our data strongly suggest that the immune related mechanisms taking place at this early age in the PLN, correlate with homeostatic changes influencing tissue integrity of the adjacent pancreatic tissue. Functional analysis of the identified genes suggested that similar mechanisms might be operating during pre-inflammatory processes deployed in tissues i) hosting parasitic microorganisms and ii) experiencing unrestricted invasion by tumour cells.
Type 1 diabetes (T1D) is an autoimmune disease characterized by the absence of insulin due to the specific destruction of the insulin-producing β cells of the pancreatic islets. This is a progressive process taking over 20 weeks in the NOD mouse and several years in human patients to be completed . The NOD animal model has been a valuable source of information for several aspects of disease pathogenesis . Genetic studies have contributed to portray the complexity of the disease and have established that multiple loci are carrying genes implicated in T1D in human  and animal models . In human, more than 6 genes contribute to the disease [5, 6] while over 20 loci have been described in the NOD mouse  but only few possible candidate genes have been unequivocally defined  other than the H2g7 Idd1 locus [7, 9].
Despite intensive research, the initial causal events remain elusive since the selection of individual mice at early stages, prior to the overt clinical signs, represents a challenge. Indeed, even though the NOD mice are inbred not all animals develop the disease, with an incidence of 40-90% in females, depending on the colony . Additional hindrances for selecting individual animals that will subsequently develop the disease with certainty are the low penetrance of the implicated genes and the influence of environmental factors. For these reasons, the exact mechanisms taking place prior to the onset of the pancreatic islet-damaging sustained inflammatory processes remain largely unknown.
The aim of our investigation was to evaluate the possibility of the existence of distinct gene expression profiles in order to eventually render possible the study of the early molecular changes, taking place before the onset of inflammation in autoimmune prone mice. One hypothesis is that during the post-weaning period in genetically autoimmune prone individuals (mouse or human), homeostatic changes prompt the immune system not to conform to physiological responses, but to instead trigger pathways that lead to the final autoimmune condition in later life. Exogenous factors including the more or less sustained presence of micro-organisms in otherwise microbe-free tissues may take place during early life stages . This may influence the immune regulatory processes put in place to prevent, attenuate and/or repair the acute or sustained inflammation. It has been reported that early immune stimulation prevents autoimmune diseases, while infections later in life might exacerbate their advent in genetically prone individuals . A single injection of BCG (the Mycobacterium bovis vaccine strain) is highly protective against T1D in the NOD mice  and BB rats  when administered early but it has no beneficial effect given after disease onset; similar to what has been observed in human . These observations lead to the indication that the time "window" of opportunity to prevent disease in human may share similarities with the NOD mice and thus BCG or other immune-system signalling, safe biological preparation, should be given as therapy before the onset of inflammation of the pancreatic islets (insulitis). This early "window" coincides in the NOD mouse model with the early appearance of insulitis, occurring between 4 and 7 weeks  and led to the hypothesis that early initiating mechanisms should be evaluated at this age i.e. over the weaning period. We have previously established that IAA are present early between 3 and 5 weeks of age and represent a diabetes sub-phenotype since their presence in NOD mouse sera, at this early period, correlated with the final disease phenotype . Thus Early-IAA (E-IAA) have been shown to mark the first measurable phenotypic checkpoint in T1D pathogenesis. Animals positive for E-IAA at 3-5 weeks after weaning, develop autoimmune diabetes earlier (between 16 and 20 weeks) than their E-IAA negative littermates . We have used this E-IAA sub-phenotype to select individual mice as diabetes prone and search for changes in gene expression patterns in the PLN, correlating with the presence of E-IAA in the NOD mouse sera.
Overt T1D is considered to require the dominant presence of islet-destructive pro-inflammatory T lymphocytes . The onset of T cell activation is known to deploy within organized secondary lymphoid tissues, draining the tissue that delivers the adequate signals for activation of T and B lymphocytes.
In the NOD mice, any early islet-disruptive remodelling process taking place, whether it implies β cell apoptosis or repair mechanisms , might result in delivery of immune cell-signalling to the PLN. The endpoint of such signalling might include deletion of autoreactive regulatory T cells, and/or expansion of autoreactive pro-inflammatory T lymphocytes as well as T cell-dependent autoreactive B lymphocytes. The later further differentiate as insulin auto-antibody secreting plasma cells . Antigen specific T cells interact with dendritic cells (DC) within intact explanted lymph nodes (LN), and it has been observed that immunological synapse formation and prolonged interactions between these 2 cell types is followed by activation, dissociation and rapid migration of T cells away from the antigenic stimulus area . Even more strikingly, excision of PLNs at 3 weeks of age protects NOD mice against IAA, insulitis and diabetes development, almost completely, but has no effect when performed at 10 weeks of age . Therefore the time frame of 3 to 5 weeks seems to correspond to the initiation of the autoimmune process in the PLN and the related transcriptional profiling can be potentially sort out in this tissue.
We anticipated that at 5 weeks of age, differential gene expression correlated with the presence or absence of E-IAA, in the PLN of NOD mice, could be detected. This timing allowed for prior tests of E-IAA at 3 and 4 weeks. We used microarray analysis to determine gene expression patterns in the PLN. This approach has the potential i) to provide a global profiling of the modulated genes at this first measurable phenotypic checkpoint in T1D pathogenesis, ii) to allow novel hypothesis to be formulated concerning the initiation of the islet-destructive autoimmune process and eventually to further challenge these hypothesis. We identified ectopic high expression of insulin genes in the E-IAA positive PLN. This data confirm the possible cross talk between the pancreas and its adjacent lymph nodes via at least one pancreas specific autoantigen.
Functional annotations of up-regulated transcripts in the PLN of E-IAA positive NOD mice revealed a gene network otherwise found also to operate during inflammatory and tumour processes. Furthermore, we have identified several transcripts coding for extracellular proteins that have to be evaluated as potential disease pathogenesis associated markers and possibly therapeutic targets. This data represent a "freeze frame" of a first identifiable disease check point, associated with the presence of IAA and corresponding to the early steps of the autoimmune process.
Mice were purchased from Taconic farms (NOD/tac). Pregnant females were tested one week before delivery for the presence of IAA. To discard any possibility of trans-placental antibodies affecting our experimental design, based upon the spontaneous appearance of IAA, we have used, when possible, individual mice from progenies of negative females and from litters that they contained both positive and negative individuals.
Animals were tested as previously described , at 3, 4 and 5 weeks of age for the presence of IAA. PLN of IAA positive and negative animals were used for RNA preparation and pancreata have been also isolated for histology. All animals were kept under SPF conditions and all studies were performed under the recommended Laboratory Animal Care committee of the University of Colorado.
IAA were measured in a standard radioimmunoassay as previously described [15, 21] incorporating competition with unlabeled insulin and precipitation with Protein A/G sepharose in a 96 well filtration plate.
Tissues were formalin (10%) fixed and embedded in paraffin. Sections of paraffin-embedded tissues were stained with Hematoxylin & Eosin and with polyclonal guinea pig anti-insulin antibodies (Millipore) followed by incubation with a peroxidase-labeled anti-guinea pig IgG antibody (Kierkegaard & Perry Laboratories Inc.). The presence of insulitis has been evaluated by visualizing more than 10 islets in 3 different sections from the same individual. Similarly the presence of insulin expressed in the PLNs has been evaluated in more than 3 sections of each individual tissue.
Tissues and RNA preparation
PLN have been dissected from six E-IAA negative and three E-IAA positive animals at 5 weeks of age. Tissues were immediately placed in RNA later containing RNases inhibitors (Qiagen) and processed for total RNA preparation on the same day with Qiagen mini-RNA commercial kits (RNeasy) with on-column DNase treatment following the manufacturer's protocol. Quantity, purity and integrity were confirmed initially by spectrophotometry (A260/A280 ratio) and capillary electrophoresis (2100 Bio analyzer; Agilent Technologies, Palo Alto, CA). 4.5 μg of total RNA were used for the target preparations and no additional step of cDNA amplification has been used in our studies.
Target preparations and chip hybridizations have been performed as recommended by the manufacturer MG_U74A_version 2 GeneChips (Affymetrix, Santa Clara, Ca). Each chip contained 12 486 probe sets (Supplemental Materials and Methods). All gene expression data discussed in this report have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession number GSE15582 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE15582.
Chromosomal location enrichment by Map genes to chromosomes
Genes used for hierarchical clustering have been also used to create the Genome view of our transcriptome data by using the dChip software (map of genes tool), as described in the dChip site http://www.dchip.org.
P-values have been calculated for all stretches containing ≤20 selected genes to assess the statistical significance of gene proximity and the significant P-values have been reported in the analysis view. The "tightness" i.e. the rank distance of the genes on the two ends of the stretch was calculated by/n genes against that of n genes randomly put on the chromosome as described http://biosun1.harvard.edu/complab/dchip/map%20chromosome.htm.
Gene Ontology analysis was carried out using DAVID NIAID/NIH online functional annotation tool and ontology terms were retained according to the significance of P-values for any particular GO term. This indicated that there were more genes associated with this GO term than would be expected by random chance . The results were independently confirmed by applying the GO Browser tool in NetAffx to the same set of differentially expressed genes. The associations between altered genes or pathways were further evaluated using the Ingenuity Pathways Analysis software (IngenuityR Systems).
Reverse transcription-coupled real-time PCR
Total RNAs from PLN of E-IAA positive and negative mice were used for Real Time PCR following the ABI software and method. To confirm gene expression levels identified in the arrays by RT-PCR, RNAs were aliquoted from the same samples used for the array hybridizations. Sequences corresponding to the investigated genes have been retrieved from the public database and aligned with Affymetrix target sequences using BLAST software http://blast.ncbi.nlm.nih.gov/Blast.cgi. Primers and TaqMan probes [see Additional File 1] have been designed using the proprietary TaqMan software. We verified that the amplified fragments corresponded to the exact Affymetrix sequences. mRNA from an insulinoma cell line has been used as positive control and 18S rRNA for data normalization.
E-IAA as a marker of autoimmune-related modifications in the PLN of NOD mice
Subphenotypes of PLN used for transcriptome
Gene expression patterns
E-IAA related gene expression patterns in PLN
Our interest in this study was to identify gene expression variations solely due to the presence or absence of IAA, in a set of identical animals of the same age and same genetic background, which however are not synchronized for the evolution of the disease. Albeit the NOD strain is inbred, homeostatic differences inherent to each animal might influence the physiopathological status and in consequence the disease stage. Therefore, RNAs have been prepared from PLN, isolated from individual mice and used without prior amplification. This allows identifying and controlling for inherent biological variability of gene expression, between individual animals that might not be related to the E-IAA sub-phenotype and allowed us avoiding the possibility for any preferential transcript amplification.
Statistically significant differences in gene expression levels between the E-IAA positive and E-IAA negative lymph nodes were determined (adjusted P-values < 0.05 and log2 fold changes >0.42). A total of 177 probe sets were deemed statistically significant by these criteria and corresponded to 165 unique genes, whilst 12 transcripts were identified repeatedly because they were represented on the chip by more than one probe set [see Additional File 3]. Unsupervised hierarchical clustering of the modulated probe sets confirmed clear sample segregation in the two phenotypic groups (Figure 2). Gene expression patterns correlated well with phenotypic evaluation for all individual mice except one: sample A36.4 is classified with the positive group (Figure 2, noted by *) whilst its phenotypic evaluation, as described above (Table 1), has identified it as a negative sample. This discrepancy, as discussed above, correlates with previous findings in human and the NOD mouse, whereas autoantibody negative individuals have still a probability to develop the disease [25, 27].
The biological relevance of these trends in gene expression changes will obviously need to be tested on an individual gene basis.
Validation of microarray data by Real time PCR and immunohistochemistry
Expression of insulin genes in the PLN was not expected since insulin is synthesized, stored and secreted by the β cells of the pancreatic islets in a highly regulated manner in human and mouse. However ectopic expression or so-called illegitimate transcription of self-antigens has been previously reported, mainly in the thymus  and spleen . Cells expressing pancreatic autoantigens are also found in peripheral lymphoid tissues and the same cell can express more than one self antigen [15, 30]. The functional significance of these observations has been questioned and evidence from studies in the mouse strongly suggested that self-antigen-expressing cells, in the thymus, may participate in negative selection and mediate tolerogenic signals, some of which may culminate in the apoptotic death of autoreactive lymphocytes .
Considering this data, we have stained NOD PLN with anti-insulin antibodies by immunohistochemistry (Figure 3F) and used pancreas as positive control (Figure 3E). Sparse distribution of insulin expressing cells was noted in the PLN (Figure 3F), while no staining has been obtained in inguinal lymph nodes used as negative controls (Figure 3G), thus confirming the ectopic insulin gene expression observed in our data. It is interesting to note that insulin genes have been also found expressed in the PLN of NOD mice at 4 weeks in a recent publication  although at lower levels.
Noteworthy that except insulin no other endocrine-specific transcripts (somatostatin, glucagon, somatostatin receptor or Islet amyloidal polypeptide) have been identified in our data set [see Additional File 4]. As expected and similarly to the endocrine specific transcripts, no other transcripts known to be specific either only to the exocrine pancreas or to the developing pancreas were expressed in the PLN at 5 weeks [see Additional File 4]. While expression of insulin genes in the PLN is not a novel finding, the high levels of expression found in the PLN of E-IAA positive NOD mice was not expected and is bewildering. One explanation might be that in the E-IAA positive PLN, APCs expressing insulin might act, locally in the PLN, as insulin-presenting cells. This would contribute to lack of tolerance to this molecule and consequently to the instigation of insulin autoreactive lymphocyte trafficking, leading to infiltration of the islets resulting to T1D. However a role of insulin in local tissue architecture, independent of its autoantigen function, cannot be excluded.
Transcriptional signatures specific to E-IAA sub-phenotype
An important aspect of our experimental design concerns the early selection of NOD mice by the presence of E-IAA, used as a marker of early autoimmune activity, at an age whereas homeostatic changes in the pancreas might take place. It has been demonstrated that increased neonatal beta cell apoptosis is taking place with a peak at 2-3 weeks after birth, in T1D animal models including the BB rat and the NOD mouse . The cross-talk between the pancreas and the adjacent LN during or immediately after this phase might be reflected by discrete gene expression changes in this later tissue.
We undertook a series of exploratory functional data analysis aiming to comprehend the roles of the identified gene expression variations in E-IAA-dependant islet autoimmunity. Forty-one out of the 119 up-regulated genes showed differential expressions over 5 fold [Additional File 3]. This indicates that the autoimmune process as reflected by the presence of IAA in the positive NOD mice is marked by local gene expression changes in the PLN. This might include changes in the immune process, occurring in this lymphoid tissue and/or changes in local tissue architecture and remodelling, not only of the lymph node but also of the pancreas, as mentioned above, assisted or mediated by the genes expressed in the PLN.
i. Cellular compartment localization
Annotations chart report for over-representation of functional categories.
UP & DOWN REGULATED
Gene Count (Ease score P-value)
60 (9E-22) ↑, 2 ↓
11 (0.02) ↑
8↑ 6 (0.02) ↓
35↑, 18 (0.03) ↓
10 (0.04) ↑
65 (1.4E-5) ↑, 9 ↓
41 (8.9E-8) ↑, 3 ↓
27 (2.9E-10) ↑ 1 ↓
Serine-type peptidase activity
25 (8.3E-19) ↑
15 (8.2E-4) ↑ 7 (0.02) ↓
Calcium binding **
10 (NS) ↑
6 (0.0008) ↑
3 (0.05) ↑, 3 ↓
5 (3.6E-4) ↓
Tissue kallikrein activity
8 (8.5E-12) ↑
2 (2.5E-2) ↑
74 (6.4E-4) ↑, 17 (0.02)
68a ↑ 14 (3.6E-6)
27 (1.5E-11) ↑
26 (0.02) ↑
Lipid, Fatty acid & steroid metabolism**
15 (1.31E-4) ↑
mRNA transcription regulation**
9 (0.05) ↑
8 (3.73E-8) ↑
7 (0.01) ↑
B-cell & Antibody-mediated immunity**
9 (1.2E-4) ↓
Response to stress
18 (8.3E-5) ↑ 2.3E-2↓
16 (0.09) ↑
Immunity & Defense**
16 (0.03) ↑ 4 (0.2) ↓
15 (0.002) ↑
Response to external stimulus
14 (4.9E-5) ↑, 6 (0.5) ↓
Response to wounding
10 (5.09E-4) ↑
8 (0.001) ↑
Cytokine & chemokine mediated signaling
7 (0.06) ↑
Other Immune defense**
7 (0.05) ↑
Humoral Immune response
4 (2.2E-3) ↑
3 (1.1E-2) ↓
3 (0.008) ↓
3 (0.008) ↓
Genes coding for secreted proteins include the Spp1 gene coding for osteopontin that has been found to play also a role in autoimmune demyelinating disease (MS) in human . High levels of this pro-inflammatory cytokine have been reported in active relapsing-remitting MS patients . Several other genes coding for secreted proteins with an immune/inflammatory function, include kallikreins genes (Klk1b22, Klk1b9, Klk1b5, KlK1rb9, Klk1b1, Klk16) that have been described as potential biomarkers for certain cancers . Members of the kallikrein large family that belongs to serine-type proteases have been implicated in autoimmune diseases, in particular in a model of Sjögren's syndrome .
In contrast to the high representation of secreted proteins coding genes, nuclear proteins transcripts are more rare in our data set, representing 11% of the total number of the modulated genes (18/165) (Figure 4a). This could indicate that the high levels of mRNA found in the E-IAA positive PLN reflect a natural steady state process of gene transcription, initiated prior to 5 weeks and/or regulation of the mRNAs is taking place at the post transcriptional level.
ii. Functional significance of the PLN gene signatures
Aiming to correlate relevant functions of our transcriptome data with the E-IAA sub-phenotype and the end T1D phenotype, we assessed for functional significance of the identified genes. We applied two different approaches: first the identified differentially expressed genes between E-IAA positive and negative NOD mice have been subjected to Gene Ontology analysis (GO, EASE)  and evaluated for classification according to functional annotations for over-representation, in combination with assessment by PANTHER terms resource . Thus in each category are included not only genes belonging to known functional categories by GO terms but also genes for which annotations are attributed by existing published experimental evidence together with consideration of evolutionary relationships (PANTHER).
In this analysis, terms with the highest levels of significance have been retained (for GO 8.35E-5<p < 0.07 and for PANTHER 1.2E-4<p < 0.02). Several genes have been annotated to Immune related process (Table 2) and several belong to the serine proteases group and include inflammation and immune functions (Table 2) and [see Additional File 5]. This group of proteins is characterized by the presence of a serine in the catalytic domain that is highly conserved between its members. Disorders of the dynamic balance of proteases and their inhibitors have been reported in the blood stream and in lymphocytes of diabetes patients and it has been proposed that increase of proteases concentrations represents an early marker of T1D . Imbalance of the proteolytic system might also change the local homeostasis and interfere with the extracellular matrix of cells in the pancreas that become more permissive to lymphocyte invasion as it has been reported for tumour metastasis . Further experimental evidence is required in order to comprehend the significance of the presence of serine proteases in the PLN of NOD mice.
Putative functional categories of transcriptional signatures in E-IAA positive versus E-IAA negative PLN
Reg2, Reg1, Try4, Pap, Tff2, Klk1, Try10/Try4, Reg3a, Klk1b22, Klk1b9, Klk1b5, KlK1rb9, Klk1b1, Cuzd1, Pla2g1b, Mt1, Nupr1, Mt2, Spp1, Klk1b16, Krt18, Ahsg, Krt8, Vtn, Klk1b5, Plg , Klk1b26, Klk1b4, S100a9, Reg3γ, Hspa1a, Retn, S100a8, Slc23a3, Ghr
Spp1, Rnase 1, Sycn, Dmbt1, Ela1, Ela2A, Egfbp2, Ang, Fga, Crlf1, Crlf2, C1qb, Rnase4, Cfd, Ctrl, Gc , Egf, Ngp
Ighg, Igl-V1, AI324046/IgG1c, Rap1a, Ms4a6b, Crlf3, Cd55, Ctla2b, Cd1d1, Myo5a, Cd74, Cr1/Igk-C, LOC100046793 similar to IgkV28
Cell death and apoptosis
Rnase1, Ins2, Alb, Mt1, Nupr1, Mt2, Spp1, Muc1, Krt18, Hba-a1, Ins1, Hpn, Krt8, Vtn, Plg, Alas2 , Gjb2, Egf, Lrpap1, Gstz1, Ptger3, Fzd3, S100a9, Dsp, Hspa1a, Aass, S100a8, Hbb-b2, Ghr, Zbtb16, Cyp2e1
Pcna, Cd55, Hnrpc, Cd1d1, Cd74, Oaz1, Eef1a1
Prss3, Prss2, Dmbt1, Tff2, Mt1, Nupr1, Mt2, Spp1, Muc1, Ggh/LOC667301, Arhgdig, Tmem97, Itih4, Cyp3a11, Ang, Krt18, Gc, Ahsg, Sel1l, Hpn, Vtn, Actb, Plg , Gjb2, Egf, Lrpap, Lpl, Asns, Fzd3, S100a9, Dsp, Lcat, Acta2, Tmem56, Rhou, Ndrg2 similar, Ces3, S100a8, Ghr, Zbtb16, Tmem30b, Cyp2e1
Hmgn1 , Top2b, Rap1a, Usp1, Pcna, Cd55, Cd1d1, Rps3, Cd74, Oaz1
Prss2, Ins2, Reg3a, Ins1, Alb, Spp1, Ang, Krt18, Gc, Fga, Krt8, Vtn, Plg , Egf, Lrpap1, Lpl, Ptger3, S100a9, Acta2, S100a8, Ghr
Cd74, Cd1d1, Cd55
Prss3, Ins2, Tff2, Ins1, Cckar, Pla2g1b, Spp1, Muc1, Ang, Gc, Fga, Sel1l, Hpn, Vtn, Actb, Plg,, Ctrb1 , Gjb2, Egf, Lrpap1, Aldh1a7, S100a9, Lcat, S100a8, Zbtb16
Nde1 , Top2b, Cd55
Pnliprp1, Amy2-1, Ins2, Cel, Clps, Ins1, Amy1, Pnliprp2, Alb, Fkbp11, Rbp4, Mod1 , Cpn1
Regeneration & Remodelling
Reg2, Reg1, Pap, Reg3a, Spp1, Itih4, Hbb-b2, Hba-a1, Itih2 , Reg3γ, Ptger3, Zbtb16
Mt1, Nupr1, Mt2, Spp1, Krt18 , Hbb-b2, Krt8, Vtn, Plg, Egf, Lrpap1, Gstz1, Lpl, S100a9, Dsp, Acta2, Cyp2e1
Ins2, Ins1, Pla2g1b, Mt1, Mt2, Spp1, Krt18, Krt8, Vtn , Egf, Dsp, Rhou, Ghr
Nde1 , Rap1a, Myo5a
Znrf2 , Hnrpc, Hnrpr, Usp1, Trim59, Ube2e1
Nupr1, Ang, Pcbd1, Hbb-b2, Rnase4 , Foxa3, Zbtb16
Hmgn1 , Spic, Sp4
Nucb2, Olfr93, Car3, Ero1lb, Pah, Gcat, Gv-1 , Copz2, Cdo1, Xist, Pck1
Slx/Xmr , Slain1, Acsl4
AI593999, AE000663, Cckar , 1810009J0RIK, Arhgdig, Gnmt , Pcbd1 , M26005, AV171666, Gamt , AE000663, Gatm , Thrsp, AV250694, AW049643, AV067171
Gt(ROSA)26Sor, Pcid2, AA590345, AI842858, AI848107
In a second approach we used literature mining together with Ingenuity functional annotation tool (Ingenuity® Systems). Combination of the resulting functional annotations for all the genes in the list [see Additional File 6] allowed assignment of more than one biological function to certain genes and gave an overall classification in twelve main functional categories (Table 3). Genes that did not clearly belong to any of the categories were classified together with ESTs, as 'unknown' and comprised less than 12% of the modulated genes (Table 3). In our data set ESTs represented 7% of the identified genes, while in the chips, unknown expressed sequences represent approximately 50% of the 12,486 probe sets spotted (MG_U74A_v2 arrays).
Several transcripts belonged to more than one category (Table 3 and Additional File 6). Given the essential role attributed to the PLN in the activation of islet-specific immune response , it was assumed, at this pre- or peri-inflammatory stage (5 weeks of age) in the NOD mice, that several chemokines or adhesion molecules would be modulated in the autoantibody positive mice. However this has not been the case. While several inflammatory mediators were consistently expressed in all PLN from E-IAA positive and negative mice, few of these molecules were expressed in a significantly regulated manner (data not shown). This concerned genes coding for several chemokines (CCR7, CXCl12, CCR5, CXCl13, CCl21a and CXCl9) that although gave present calls (P) to all samples, did not meet statistical criteria for significant differences between the two groups (data not shown). One hypothesis might be that although all NOD mice possess the necessary requirements (i.e. presence of several chemokines in all samples) for the initiation of the inflammatory process against the pancreatic islets, the autoimmune-driven destructive events are probably dependent upon the levels of pancreatic tissue remodelling taking place. This in turn could be due to either antibody bearing B cells (IAA) or T lymphocyte specific trafficking against pancreatic signals.
Genome view of the transcriptome signatures
Genes on the chromosome 6 are placed in 3 families including the lithostatin family, immunoglobulin light chain and several serine type proteases coding genes. Genes on chromosome 7 concerned mainly the kallikrein coding genes and the insulin-2 gene. Interestingly 4 T1D Idd loci have been described on these two chromosomes by genetic linkage analysis: two loci on mouse chromosome 6, the Idd6  and Idd19  and two Idd loci on mouse chromosome 7: the Idd7  and the Idd27 .
Autoimmune diabetes in human is a heterogeneous disorder whereas several genes influence the apparent clinical phenotype together with environmental factors and homeostatic differences inherent to the organism. This leads to phenotypic variations of the disease observed not only in human but also in the NOD mouse, despite the genetic homogeneity of this inbred strain. Indeed autoantibody responses vary between individual mice by the titers and by the age of appearance. Similarly temporal variation in the severity of the islet-damaging inflammatory process is observed as well as variations in the onset of appearance of the end point disease.
In the present study we used the early presence of IAA in the NOD serum, as the earliest measurable sub-phenotype, aiming first to address for changes occurring in the pancreas-draining lymph nodes, when the autoimmune process has been just initiated. The PLN is a structured secondary lymphoid organ where pancreas derived cells and signals are expected to be delivered. The resulting output of the processing of these signals is either the maintenance or the loss of pancreas integrity. Our second aim was to attempt a stratification of disease check points by studying the first possibly identifiable disease stage, as selected by the presence of E-IAA.
This rational, potentially, allows identification of early genetic components that due to their specific expression in the PLN, at a precise time, may not be identified by genetic studies whereas the final T1D phenotype is usually taken in consideration.
Biological significance of the identified transcriptional signatures and relevance to T1D phenotypes
Phenotypic stratification, as delineated by sub phenotypic attribution for the identification of gene expression variation in autoimmune diabetes, has not been used, to our knowledge, prior to our study. Two recent reports have also used PLN in a similar methodology approach, at early ages in the NOD mice, without however prior phenotypic evaluation [32, 51]. One group studied 4 and 6 weeks old NOD/BDC 2.5 mice , the other group used NOD mice of various ages starting at 10 days and including 4 weeks, to compare with NOD.B10H2b mice .
Tissue regeneration and remodelling
The expression of five genes belonging to the regenerating islet-derived gene family (Pap, Reg1, Reg3a, Reg3g and Reg2), implicated in tissue integrity, maintenance and regeneration were highly up-regulated in our data set [see Additional File 6]. They code for a set of related but distinct proteins . The presence of these proteins in the E-IAA samples is intriguing, especially in the light of recent published reports demonstrating that adult human islets possess a remarkable degree of morphogenetic plasticity . Similarly to Reg genes, Itih-4, a serine-type protease inhibitor with anti-apoptotic and matrix stabilizing functions, is also up-regulated and reported to play a role in liver development and regeneration . These transcriptional signatures, if sorted out as such, might represent a new tool in the understanding of islet neogenesis and tissue preservation by inherent factors and could be valuable against β-cell destruction.
Sustained inflammatory processes and tumour progression
Several of the modulated transcripts belong to complex networks of gene products acting in tissues during inflammatory sustained processes (ie. Spp1, Klks, Mt2), or playing a role in tissue regeneration (ie. Regs) and have been found also to be associated with neoplasia (Table 3), [see Additional File 6]. Certain kallikrein coding genes have been up-regulated in our data set and are considered possible biomarkers for cancer . Similar expression patterns have been found for Muc1 gene that is highly expressed in gastric carcinoma and proposed to be involved in gastric carcinogenesis  as well as for the genes Cckar, Ggh [60, 61] and Ang  that were reported to be implicated in the genesis of various tumours. Ang, in particular, is involved in the Notch signalling pathway related to neurovascular progression of pancreatic cancer and metastasis , as well as in multiple myeloma . It has been reported that circulating levels of angiogenic cytokines can predict tumour progression and prognosis in neuroendocrine carcinomas of the gastro-entero-pancreatic system .
Finally, the gene coding for nucleobindin 2 (Nucb2) implicated in B cell lymphomas  was also up-regulated. It is a DNA binding protein, called also NEFA and its sequence contains a signal peptide, suggesting that it is a secreted or trans-membrane protein .
Infection and inflammation
Several genes known to play a role in various inflammatory conditions were differentially modulated in our data set (Table 3). Such molecules include two elastase genes, Ela1  and the neutrophil elastase Ela2 , known to play a role in peripheral tissues invaded by microbes , the Spp1 gene , and the secretory granule protein syncollin (Sync), expressed by neutrophils with a role in host defence to infections by invading bacteria, fungi, and protozoa . Sync is also expressed in the exocrine pancreas and in the duodenum and colon . Other gene products operating during the inflammatory process include the Egfbp2, part of the kallikrein gene family , the actin binding protein (Actb), the mentioned above Muc1 that is part of signalling proteins regulating cell adhesion/de-adhesion  and two genes coding for metallothionein, Mt1 and Mt2. Increased levels of Muc1 have been found in the small intestine of a cystic fibrosis mouse model . Over expression of metallothioneins has been described in gastric cancer and in intestinal metaplasia and dysplasia . Mice lacking metalothionein are more susceptible to Helicobacter pylori colonization and gastric inflammation, indicating a protective effect of these proteins against H. pylori-induced gastritis .
Interestingly, several lipases coding genes that have been associated with diabetes were up-regulated in the E-IAA positive PLN (Table 3). Pla2g1b gene has been described as a pro inflammatory molecule expressed in lung epithelium hosting Pseudomonas aeruginosa . It is also implicated in glucose uptake by the liver, heart and muscle tissues under high glucose diet . Inhibition of lipase activity was shown to reduce the incidence of type 2 diabetes mellitus and genetic variations in the Clps and Pnlip genes are associated with type 2 diabetes .
Two additional genes coding for pancreatic lipase related proteins Pnliprp1 and Pnliprp2 were highly up-regulated in the E-IAA PLN (Table 3) and [see Additional File 3]. Pnliprp2 protein plays a crucial role in the digestion of dietary fats in milk suckling mammals  while it seems to participate in T cell cytotoxicity. A T cell lipase shares 74% amino acid identity with the pancreatic Pnliprp2 .
Could novel hypothesis be further explored in NOD mice at early post birth period?
Overall the data reported herein suggest that the events leading to islet-damaging autoimmune destruction might be rooted to early processes, some sharing common characteristics with several other inflammatory conditions found also in microbes-hosting tissues (i.e. Pla2g1b, Mt1) and in certain cancers (i.e. Klk genes, Muc1, Reg genes). Noticeably among the cytokine receptor-like factor 1 and 3 (Crlf1, Crlf3) that belong to the same family , Crlf3 has been found to be up-regulated in skin cancer . While Crlf1 expression is up-regulated in the E-IAA positive samples, expression of the Crlf3 gene is down-regulated. Crlf1 mRNA has been reported to be up-regulated by TNF-α, IL-6, and IFN-γ . Polymorphisms in the regulatory region of C1qb gene have been correlated with down-regulation of the murin c1q protein levels and linked to lupus nephritis .
The exact order of events in autoimmune initiation requires additional studies in order to be elucidated. Our hypothesis implies that defects in pancreas remodelling processes, taking place at the weaning period, might be a trigger for immune deregulation leading to autoimmune tissue destruction. This might include tissue regeneration or damage. The developmental stage of the organism as well as agents exogenous to the organism, such as invasive microbes, potentially can contribute to local homeostatic changes. It can be assumed that the identified genes have the potential to re-address the question of the link between environmental triggers, such as sustained microbial signalling processes, with the autoimmune process, at least in T1D.
This hypothesis requires certainly additional integrative studies and experimental evidence to be validated. It is noticeable that none of the common pro-inflammatory molecules such as tumour necrosis factor or interferons have been found in our group of genes. It has been recently reported 8 IFN-α inducible genes to be differentially up-regulated in CD4+ T cells isolated from the PLN of 6 wks, when compared with 2 weeks old NOD/BDC 2.5 mice . These authors concluded that IFN-α initiates T1D in the NOD mice and suggested that pDCs and TLRs play a role in T1D initiation. Though NOD mice differ from strains unable to initiate islet autoimmunity in terms of their response to IFN, our data indicate that such differences are not associated, within the NOD strain, with early expression of IAA, which in its turn is associated with greater risk of diabetes progression, at least at the age studied in this report.
The autoimmune-related transcriptional early landscape of the PLN, in a first glance, seems difficult to be contained due to its complexity. The stepwise nature of the process resulting in T1D, is illustrated in our study by the difficulty to sort out a clear mechanism describing a possible scenario for disease initiation and implicating, preferably, a small number of molecules. Nevertheless, the gene signatures identified by this global analysis of the PLN in the autoantibody positive NOD mice at 5 weeks, together with the functional annotations described herein, represent a valuable tool in designing additional experiments for understanding the irreversible initiation of the autoimmune process in T1D. Experimental investigation of the functions of the identified genes together with further integrative analysis studies have the potential to sear light into the exact autoimmune-related processes that take place at this early age, in autoimmune prone, insulin autoantibody positive animals.
type 1 diabetes
Non Obese Diabetic
pancreatic lymph node
Antigen Presenting Cells
- Idd :
insulin-dependent diabetes locus.
We thank Rebecca Brown for assistance to initial analysis of the data and Geneviève Milon and Eric Prina for critical reading of the manuscript. EM is endeavoured to Genevieve Milon for support and encouraging discussions. This work has been supported by the BDC, University of Colorado and the Institut Pasteur, Paris.
- Leiter EH, Prochazka M, Coleman DL: The non-obese diabetic (NOD) mouse. Am J Pathol. 1987, 128: 380-383.PubMedPubMed CentralGoogle Scholar
- Anderson MS, Bluestone JA: The NOD mouse: A Model of Immune Dysregulation. Annu Rev Immunol. 2005, 23: 447-485. 10.1146/annurev.immunol.23.021704.115643.View ArticlePubMedGoogle Scholar
- Burton PR, Clayton DG, Cardon LR, Craddock N, Deloukas P, Duncanson A, Kwiatkowski DP, McCarthy MI, Ouwehand WH, Samani NJ, et al: Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007, 39: 1329-1337. 10.1038/ng.2007.17.View ArticlePubMedGoogle Scholar
- Melanitou E: The autoimmune contrivance: genetics in the mouse model. Clin Immunol. 2005, 117: 195-206. 10.1016/j.clim.2005.07.007.View ArticlePubMedGoogle Scholar
- Todd JA, Walker NM, Cooper JD, Smyth DJ, Downes K, Plagnol V, Bailey R, Nejentsev S, Field SF, Payne F, et al: Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nat Genet. 2007, 39: 857-864. 10.1038/ng2068.View ArticlePubMedPubMed CentralGoogle Scholar
- Baschal EE, Eisenbarth GS: Extreme genetic risk for type 1A diabetes in the post-genome era. J Autoimmun. 2008, 31: 1-6. 10.1016/j.jaut.2008.03.003.View ArticlePubMedPubMed CentralGoogle Scholar
- Wicker LS, Todd JA, Peterson LB: Genetic control of autoimmune diabetes in the NOD mouse. Annu Rev Immunol. 1995, 13: 179-200. 10.1146/annurev.iy.13.040195.001143.View ArticlePubMedGoogle Scholar
- Wicker LS, Clark J, Fraser HI, Garner VE, Gonzalez-Munoz A, Healy B, Howlett S, Hunter K, Rainbow D, Rosa RL, et al: Type 1 diabetes genes and pathways shared by humans and NOD mice. J Autoimmun. 2005, 25 (Suppl): 29-33.View ArticlePubMedGoogle Scholar
- Tisch R, McDevitt H: Insulin-dependent diabetes mellitus. Cell. 1996, 85: 291-297. 10.1016/S0092-8674(00)81106-X.View ArticlePubMedGoogle Scholar
- Akerblom HK, Knip M: Putative environmental factors in Type 1 diabetes. Diabetes Metab Rev. 1998, 14: 31-67.View ArticlePubMedGoogle Scholar
- Singh B: Stimulation of the developing immune system can prevent autoimmunity. J Autoimmun. 2000, 14: 15-22. 10.1006/jaut.1999.0349.View ArticlePubMedGoogle Scholar
- Silveira PA, Baxter AG: The NOD mouse as a model of SLE. Autoimmunity. 2001, 34: 53-64.View ArticlePubMedGoogle Scholar
- Sadelain MW, Qin HY, Sumoski W, Parfrey N, Singh B, Rabinovitch A: Prevention of diabetes in the BB rat by early immunotherapy using Freund's adjuvant. J Autoimmun. 1990, 3: 671-680. 10.1016/S0896-8411(05)80034-4.View ArticlePubMedGoogle Scholar
- Elliott JF, Marlin KL, Couch RM: Effect of bacille Calmette-Guerin vaccination on C-peptide secretion in children newly diagnosed with IDDM. Diabetes Care. 1998, 21: 1691-1693. 10.2337/diacare.21.10.1691.View ArticlePubMedGoogle Scholar
- Melanitou E, Devendra D, Liu E, Miao D, Eisenbarth GS: Early and quantal (by litter) expression of insulin autoantibodies in the nonobese diabetic mice predict early diabetes onset. J Immunol. 2004, 173: 6603-6610.View ArticlePubMedGoogle Scholar
- Salomon B, Lenschow DJ, Rhee L, Ashourian N, Singh B, Sharpe A, Bluestone JA: B7/CD28 costimulation is essential for the homeostasis of the CD4+CD25+ immunoregulatory T cells that control autoimmune diabetes. Immunity. 2000, 12: 431-440. 10.1016/S1074-7613(00)80195-8.View ArticlePubMedGoogle Scholar
- Hauben E, Roncarolo MG, Nevo U, Schwartz M: Beneficial autoimmunity in Type 1 diabetes mellitus. Trends Immunol. 2005, 26: 248-253. 10.1016/j.it.2005.03.004.View ArticlePubMedGoogle Scholar
- Judkowski V, Krakowski M, Rodriguez E, Mocnick L, Santamaria P, Sarvetnick N: Increased islet antigen presentation leads to type-1 diabetes in mice with autoimmune susceptibility. Eur J Immunol. 2004, 34: 1031-1040. 10.1002/eji.200324563.View ArticlePubMedGoogle Scholar
- Stoll S, Delon J, Brotz TM, Germain RN: Dynamic imaging of T cell-dendritic cell interactions in lymph nodes. Science. 2002, 296: 1873-1876. 10.1126/science.1071065.View ArticlePubMedGoogle Scholar
- Gagnerault MC, Luan JJ, Lotton C, Lepault F: Pancreatic lymph nodes are required for priming of beta cell reactive T cells in NOD mice. J Exp Med. 2002, 196: 369-377. 10.1084/jem.20011353.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu L, Robles DT, Abiru N, Kaur P, Rewers M, Kelemen K, Eisenbarth GS: Early expression of antiinsulin autoantibodies of humans and the NOD mouse: evidence for early determination of subsequent diabetes. Proc Natl Acad Sci USA. 2000, 97: 1701-1706. 10.1073/pnas.040556697.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.View ArticlePubMedPubMed CentralGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- Melanitou E: Functional genomics in early autoimmunity. Ann N Y Acad Sci. 2005, 1050: 64-72. 10.1196/annals.1313.008.View ArticlePubMedGoogle Scholar
- Yu LRD, Abiru N, Rewers M, Kelemen K, Eisenbarth GS: Early expression of anti-insulin autoantibodies of humans and the NOD mouse: Evidence for early determination of subsequent diabetes. Proc Natl Acad Sci USA. 2000, 97: 1701-1706. 10.1073/pnas.040556697.View ArticlePubMedPubMed CentralGoogle Scholar
- Pietropaolo MEG, (Ed): Autoantibodies in human Disease. Molecular Pathology of Type 1 Diabetes mellitus edition. 2001, Basel: Karger
- Wang J, Miao D, Babu S, Yu J, Barker J, Klingensmith G, Rewers M, Eisenbarth GS, Yu L: Prevalence of autoantibody-negative diabetes is not rare at all ages and increases with older age and obesity. J Clin Endocrinol Metab. 2007, 92: 88-92. 10.1210/jc.2006-1494.View ArticlePubMedGoogle Scholar
- Derbinski J, Schulte A, Kyewski B, Klein L: Promiscuous gene expression in medullary thymic epithelial cells mirrors the peripheral self. Nat Immunol. 2001, 2: 1032-1039. 10.1038/ni723.View ArticlePubMedGoogle Scholar
- Zheng X, Yin L, Liu Y, Zheng P: Expression of tissue-specific autoantigens in the hematopoietic cells leads to activation-induced cell death of autoreactive T cells in the secondary lymphoid organs. Eur J Immunol. 2004, 34: 3126-3134. 10.1002/eji.200425177.View ArticlePubMedGoogle Scholar
- Pugliese A, Brown D, Garza D, Murchison D, Zeller M, Redondo MJ, Diez J, Eisenbarth GS, Patel DD, Ricordi C: Self-antigen-presenting cells expressing diabetes-associated autoantigens exist in both thymus and peripheral lymphoid organs. J Clin Invest. 2001, 107: 555-564. 10.1172/JCI10860.View ArticlePubMedPubMed CentralGoogle Scholar
- Hanahan D: Peripheral-antigen-expressing cells in thymic medulla: factors in self-tolerance and autoimmunity. Curr Opin Immunol. 1998, 10: 656-662. 10.1016/S0952-7915(98)80085-X.View ArticlePubMedGoogle Scholar
- Kodama K, Butte AJ, Creusot RJ, Su L, Sheng D, Hartnett M, Iwai H, Soares LR, Fathman CG: Tissue- and age-specific changes in gene expression during disease induction and progression in NOD mice. Clin Immunol. 2008, 129: 195-201. 10.1016/j.clim.2008.07.028.View ArticlePubMedPubMed CentralGoogle Scholar
- Trudeau JD, Dutz JP, Arany E, Hill DJ, Fieldus WE, Finegood DT: Neonatal beta-cell apoptosis: a trigger for autoimmune diabetes?. Diabetes. 2000, 49: 1-7. 10.2337/diabetes.49.1.1.View ArticlePubMedGoogle Scholar
- Chabas D, Baranzini SE, Mitchell D, Bernard CC, Rittling SR, Denhardt DT, Sobel RA, Lock C, Karpuj M, Pedotti R, et al: The influence of the proinflammatory cytokine, osteopontin, on autoimmune demyelinating disease. Science. 2001, 294: 1731-1735. 10.1126/science.1062960.View ArticlePubMedGoogle Scholar
- Vogt MH, Lopatinskaya L, Smits M, Polman CH, Nagelkerken L: Elevated osteopontin levels in active relapsing-remitting multiple sclerosis. Ann Neurol. 2003, 53: 819-822. 10.1002/ana.10606.View ArticlePubMedGoogle Scholar
- Diamandis EP, Yousef GM: Human tissue kallikreins: a family of new cancer biomarkers. Clin Chem. 2002, 48: 1198-1205.PubMedGoogle Scholar
- Takada K, Takiguchi M, Konno A, Inaba M: Autoimmunity against a tissue kallikrein in IQI/Jic Mice: a model for Sjogren's syndrome. J Biol Chem. 2005, 280: 3982-3988. 10.1074/jbc.M410157200.View ArticlePubMedGoogle Scholar
- Thomas PD, Kejariwal A, Campbell MJ, Mi H, Diemer K, Guo N, Ladunga I, Ulitsky-Lazareva B, Muruganujan A, Rabkin S, et al: PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res. 2003, 31: 334-341. 10.1093/nar/gkg115.View ArticlePubMedPubMed CentralGoogle Scholar
- Finotti P: The role played by serine proteases in the development and worsening of vascular complications in type 1 diabetes mellitus. Curr Diabetes Rev. 2006, 2: 295-305. 10.2174/157339906777950624.View ArticlePubMedGoogle Scholar
- Del Rosso M, Fibbi G, Pucci M, D'Alessio S, Del Rosso A, Magnelli L, Chiarugi V: Multiple pathways of cell invasion are regulated by multiple families of serine proteases. Clin Exp Metastasis. 2002, 19: 193-207. 10.1023/A:1015531321445.View ArticlePubMedGoogle Scholar
- Mendiratta SK, Martin WD, Hong S, Boesteanu A, Joyce S, Van Kaer L: CD1d1 mutant mice are deficient in natural T cells that promptly produce IL-4. Immunity. 1997, 6: 469-477. 10.1016/S1074-7613(00)80290-3.View ArticlePubMedGoogle Scholar
- Sonoda KH, Exley M, Snapper S, Balk SP, Stein-Streilein J: CD1-reactive natural killer T cells are required for development of systemic tolerance through an immune-privileged site. J Exp Med. 1999, 190: 1215-1226. 10.1084/jem.190.9.1215.View ArticlePubMedPubMed CentralGoogle Scholar
- Godfrey DI, Kronenberg M: Going both ways: immune regulation via CD1d-dependent NKT cells. J Clin Invest. 2004, 114: 1379-1388.View ArticlePubMedPubMed CentralGoogle Scholar
- Godfrey DI, MacDonald HR, Kronenberg M, Smyth MJ, Van Kaer L: NKT cells: what's in a name?. Nat Rev Immunol. 2004, 4: 231-237. 10.1038/nri1309.View ArticlePubMedGoogle Scholar
- Miyamoto K, Miyake S, Yamamura T: A synthetic glycolipid prevents autoimmune encephalomyelitis by inducing TH2 bias of natural killer T cells. Nature. 2001, 413: 531-534. 10.1038/35097097.View ArticlePubMedGoogle Scholar
- Jaakkola I, Jalkanen S, Hanninen A: Diabetogenic T cells are primed both in pancreatic and gut-associated lymph nodes in NOD mice. Eur J Immunol. 2003, 33: 3255-3264. 10.1002/eji.200324405.View ArticlePubMedGoogle Scholar
- Ghosh S, Palmer SM, Rodrigues NR, Cordell HJ, Hearne CM, Cornall RJ, Prins JB, McShane P, Lathrop GM, Peterson LB, et al: Polygenic control of autoimmune diabetes in nonobese diabetic mice. Nat Genet. 1993, 4: 404-409. 10.1038/ng0893-404.View ArticlePubMedGoogle Scholar
- Melanitou E, Joly F, Lathrop M, Boitard C, Avner P: Evidence for the presence of insulin-dependent diabetes-associated alleles on the distal part of mouse chromosome 6. Genome Res. 1998, 8: 608-620.PubMedPubMed CentralGoogle Scholar
- Serreze DV, Choisy-Rossi CM, Grier AE, Holl TM, Chapman HD, Gahagan JR, Osborne MA, Zhang W, King BL, Brown A, et al: Through regulation of TCR expression levels, an Idd7 region gene(s) interactively contributes to the impaired thymic deletion of autoreactive diabetogenic CD8+ T cells in nonobese diabetic mice. J Immunol. 2008, 180: 3250-3259.View ArticlePubMedGoogle Scholar
- Chen J, Reifsnyder PC, Scheuplein F, Schott WH, Mileikovsky M, Soodeen-Karamath S, Nagy A, Dosch MH, Ellis J, Koch-Nolte F, Leiter EH: "Agouti NOD": identification of a CBA-derived Idd locus on Chromosome 7 and its use for chimera production with NOD embryonic stem cells. Mamm Genome. 2005, 16: 775-783. 10.1007/s00335-005-0007-1.View ArticlePubMedGoogle Scholar
- Li Q, Xu B, Michie SA, Rubins KH, Schreriber RD, McDevitt HO: Interferon-alpha initiates type 1 diabetes in nonobese diabetic mice. Proc Natl Acad Sci USA. 2008, 105: 12439-12444. 10.1073/pnas.0806439105.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim HS, Han MS, Chung KW, Kim S, Kim E, Kim MJ, Jang E, Lee HA, Youn J, Akira S, Lee MS: Toll-like receptor 2 senses beta-cell death and contributes to the initiation of autoimmune diabetes. Immunity. 2007, 27: 321-333. 10.1016/j.immuni.2007.06.010.View ArticlePubMedGoogle Scholar
- Gallucci S, Lolkema M, Matzinger P: Natural adjuvants: endogenous activators of dendritic cells. Nat Med. 1999, 5: 1249-1255. 10.1038/15200.View ArticlePubMedGoogle Scholar
- Huang FP, Platt N, Wykes M, Major JR, Powell TJ, Jenkins CD, MacPherson GG: A discrete subpopulation of dendritic cells transports apoptotic intestinal epithelial cells to T cell areas of mesenteric lymph nodes. J Exp Med. 2000, 191: 435-444. 10.1084/jem.191.3.435.View ArticlePubMedPubMed CentralGoogle Scholar
- Graf R, Schiesser M, Reding T, Appenzeller P, Sun LK, Fortunato F, Perren A, Bimmler D: Exocrine meets endocrine: pancreatic stone protein and regenerating protein--two sides of the same coin. J Surg Res. 2006, 133: 113-120. 10.1016/j.jss.2005.09.030.View ArticlePubMedGoogle Scholar
- Jamal AM, Lipsett M, Sladek R, Laganiere S, Hanley S, Rosenberg L: Morphogenetic plasticity of adult human pancreatic islets of Langerhans. Cell Death Differ. 2005, 12: 702-712. 10.1038/sj.cdd.4401617.View ArticlePubMedGoogle Scholar
- Bhanumathy CD, Tang Y, Monga SP, Katuri V, Cox JA, Mishra B, Mishra L: Itih-4, a serine protease inhibitor regulated in interleukin-6-dependent liver formation: role in liver development and regeneration. Dev Dyn. 2002, 223: 59-69. 10.1002/dvdy.1235.View ArticlePubMedGoogle Scholar
- Borgono CA, Diamandis EP: The emerging roles of human tissue kallikreins in cancer. Nat Rev Cancer. 2004, 4: 876-890. 10.1038/nrc1474.View ArticlePubMedGoogle Scholar
- Silva F, Carvalho F, Peixoto A, Seixas M, Almeida R, Carneiro F, Mesquita P, Figueiredo C, Nogueira C, Swallow DM, et al: MUC1 gene polymorphism in the gastric carcinogenesis pathway. Eur J Hum Genet. 2001, 9: 548-552. 10.1038/sj.ejhg.5200677.View ArticlePubMedGoogle Scholar
- Takata Y, Takiguchi S, Kataoka K, Funakoshi A, Miyasaka K, Kono A: Mouse cholecystokinin type-A receptor gene and its structural analysis. Gene. 1997, 187: 267-271. 10.1016/S0378-1119(96)00765-2.View ArticlePubMedGoogle Scholar
- Cheng ZJ, Harikumar KG, Ding WQ, Holicky EL, Miller LJ: Analysis of the cellular and molecular mechanisms of trophic action of a misspliced form of the type B cholecystokinin receptor present in colon and pancreatic cancer. Cancer Lett. 2005, 222: 95-105. 10.1016/j.canlet.2004.09.008.View ArticlePubMedGoogle Scholar
- Pavel ME, Hassler G, Baum U, Hahn EG, Lohmann T, Schuppan D: Circulating levels of angiogenic cytokines can predict tumour progression and prognosis in neuroendocrine carcinomas. Clin Endocrinol (Oxf). 2005, 62: 434-443. 10.1111/j.1365-2265.2005.02238.x.View ArticleGoogle Scholar
- Buchler P, Gazdhar A, Schubert M, Giese N, Reber HA, Hines OJ, Giese T, Ceyhan GO, Muller M, Buchler MW, Friess H: The Notch signaling pathway is related to neurovascular progression of pancreatic cancer. Ann Surg. 2005, 242: 791-800. 10.1097/01.sla.0000189115.94847.f1. discussion 800-791View ArticlePubMedPubMed CentralGoogle Scholar
- Politou M, Naresh K, Terpos E, Crawley D, Lampert I, Apperley JF, Rahemtulla A: Anti-angiogenic effect of bortezomib in patients with multiple myeloma. Acta H aematol. 2005, 114: 170-173.View ArticleGoogle Scholar
- de Vos S, Hofmann WK, Grogan TM, Krug U, Schrage M, Miller TP, Braun JG, Wachsman W, Koeffler HP, Said JW: Gene expression profile of serial samples of transformed B-cell lymphomas. Lab Invest. 2003, 83: 271-285.View ArticlePubMedGoogle Scholar
- Petersson U, Somogyi E, Reinholt FP, Karlsson T, Sugars RV, Wendel M: Nucleobindin is produced by bone cells and secreted into the osteoid, with a potential role as a modulator of matrix maturation. Bone. 2004, 34: 949-960. 10.1016/j.bone.2004.01.019.View ArticlePubMedGoogle Scholar
- Benarafa C, Cooley J, Zeng W, Bird PI, Remold-O'Donnell E: Characterization of four murine homologs of the human ov-serpin monocyte neutrophil elastase inhibitor MNEI (SERPINB1). J Biol Chem. 2002, 277: 42028-42033. 10.1074/jbc.M207080200.View ArticlePubMedGoogle Scholar
- Adkison AM, Raptis SZ, Kelley DG, Pham CT: Dipeptidyl peptidase I activates neutrophil-derived serine proteases and regulates the development of acute experimental arthritis. J Clin Invest. 2002, 109: 363-371.View ArticlePubMedPubMed CentralGoogle Scholar
- Weinrauch Y, Drujan D, Shapiro SD, Weiss J, Zychlinsky A: Neutrophil elastase targets virulence factors of enterobacteria. Nature. 2002, 417: 91-94. 10.1038/417091a.View ArticlePubMedGoogle Scholar
- Denhardt DT, Guo X: Osteopontin: a protein with diverse functions. FASEB J. 1993, 7: 1475-1482.PubMedGoogle Scholar
- Bach JP, Borta H, Ackermann W, Faust F, Borchers O, Schrader M: The secretory granule protein syncollin localizes to HL-60 cells and neutrophils. J Histochem Cytochem. 2006, 54: 877-888. 10.1369/jhc.5A6792.2006.View ArticlePubMedGoogle Scholar
- Tan S, Hooi SC: Syncollin is differentially expressed in rat proximal small intestine and regulated by feeding behavior. Am J Physiol Gastrointest Liver Physiol. 2000, 278: G308-320.PubMedGoogle Scholar
- Olsson AY, Lundwall A: Organization and evolution of the glandular kallikrein locus in Mus musculus. Biochem Biophys Res Commun. 2002, 299: 305-311. 10.1016/S0006-291X(02)02629-3.View ArticlePubMedGoogle Scholar
- Mukherjee P, Tinder TL, Basu GD, Gendler SJ: MUC1 (CD227) interacts with lck tyrosine kinase in Jurkat lymphoma cells and normal T cells. J Leukoc Biol. 2005, 77: 90-99.PubMedGoogle Scholar
- Malmberg EK, Noaksson KA, Phillipson M, Johansson ME, Hinojosa-Kurtzberg M, Holm L, Gendler SJ, Hansson GC: Increased levels of mucins in the cystic fibrosis mouse small intestine, and modulator effects of the Muc1 mucin expression. Am J Physiol Gastrointest Liver Physiol. 2006, 291: G203-210. 10.1152/ajpgi.00491.2005.View ArticlePubMedGoogle Scholar
- Ebert MP, Gunther T, Hoffmann J, Yu J, Miehlke S, Schulz HU, Roessner A, Korc M, Malfertheiner P: Expression of metallothionein II in intestinal metaplasia, dysplasia, and gastric cancer. Cancer Res. 2000, 60: 1995-2001.PubMedGoogle Scholar
- Tran CD, Huynh H, Berg van den M, Pas van der M, Campbell MA, Philcox JC, Coyle P, Rofe AM, Butler RN: Helicobacter-induced gastritis in mice not expressing metallothionein-I and II. Helicobacter. 2003, 8: 533-541. 10.1046/j.1523-5378.2003.00174.x.View ArticlePubMedGoogle Scholar
- Agassandian M, Miakotina OL, Andrews M, Mathur SN, Mallampalli RK: Pseudomonas aeruginosa and sPLA2 IB stimulate ABCA1-mediated phospholipid efflux via ERK-activation of PPARalpha-RXR. Biochem J. 2007, 403: 409-420. 10.1042/BJ20061364.View ArticlePubMedPubMed CentralGoogle Scholar
- Labonte ED, Kirby RJ, Schildmeyer NM, Cannon AM, Huggins KW, Hui DY: Group 1B phospholipase A2-mediated lysophospholipid absorption directly contributes to postprandial hyperglycemia. Diabetes. 2006, 55: 935-941. 10.2337/diabetes.55.04.06.db05-1286.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindner I, Helwig U, Rubin D, Li Y, Fisher E, Boeing H, Mohlig M, Spranger J, Pfeiffer A, Hampe J, et al: Putative association between a new polymorphism in exon 3 (Arg109Cys) of the pancreatic colipase gene and type 2 diabetes mellitus in two independent Caucasian study populations. Mol Nutr Food Res. 2005, 49: 972-976. 10.1002/mnfr.200500087.View ArticlePubMedGoogle Scholar
- Lowe ME, Kaplan MH, Jackson-Grusby L, D'Agostino D, Grusby MJ: Decreased neonatal dietary fat absorption and T cell cytotoxicity in pancreatic lipase-related protein 2-deficient mice. J Biol Chem. 1998, 273: 31215-31221. 10.1074/jbc.273.47.31215.View ArticlePubMedPubMed CentralGoogle Scholar
- Grusby MJ, Nabavi N, Wong H, Dick RF, Bluestone JA, Schotz MC, Glimcher LH: Cloning of an interleukin-4 inducible gene from cytotoxic T lymphocytes and its identification as a lipase. Cell. 1990, 60: 451-459. 10.1016/0092-8674(90)90596-7.View ArticlePubMedGoogle Scholar
- Elson GC, Graber P, Losberger C, Herren S, Gretener D, Menoud LN, Wells TN, Kosco-Vilbois MH, Gauchat JF: Cytokine-like factor-1, a novel soluble protein, shares homology with members of the cytokine type I receptor family. J Immunol. 1998, 161: 1371-1379.PubMedGoogle Scholar
- Dang C, Gottschling M, Manning K, O'Currain E, Schneider S, Sterry W, Stockfleth E, Nindl I: Identification of dysregulated genes in cutaneous squamous cell carcinoma. Oncol Rep. 2006, 16: 513-519.PubMedGoogle Scholar
- Miura-Shimura Y, Nakamura K, Ohtsuji M, Tomita H, Jiang Y, Abe M, Zhang D, Hamano Y, Tsuda H, Hashimoto H, et al: C1q regulatory region polymorphism down-regulating murine c1q protein levels with linkage to lupus nephritis. J Immunol. 2002, 169: 1334-1339.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/63/prepub