The use of mass spectrometry in a proteome-centered multiomics study of human pituitary adenomas.

A pituitary adenoma (PA) is a common intracranial neoplasm, and is a complex, chronic, and whole-body disease with multicausing factors, multiprocesses, and multiconsequences. It is very difficult to clarify molecular mechanism and treat PAs from the single-factor strategy model. The rapid development of multiomics and systems biology changed the paradigms from a traditional single-factor strategy to a multiparameter systematic strategy for effective management of PAs. A series of molecular alterations at the genome, transcriptome, proteome, peptidome, metabolome, and radiome levels are involved in pituitary tumorigenesis, and mutually associate into a complex molecular network system. Also, the center of multiomics is moving from structural genomics to phenomics, including proteomics and metabolomics in the medical sciences. Mass spectrometry (MS) has been extensively used in phenomics studies of human PAs to clarify molecular mechanisms, and to discover biomarkers and therapeutic targets/drugs. MS-based proteomics and proteoform studies play central roles in the multiomics strategy of PAs. This article reviews the status of multiomics, multiomics-based molecular pathway networks, molecular pathway network-based pattern biomarkers and therapeutic targets/drugs, and future perspectives for personalized, predeictive, and preventive (3P) medicine in PAs.

Human growth hormone proteoform pattern changes in pituitary adenomas: Potential biomarkers for 3P medical approaches.

RELEVANCE: Human growth hormone (hGH) is synthesized, stored, and secreted by somatotroph cells in the pituitary gland, and promotes human growth and metabolism. Compared to a normal pituitary, a GH-secreting pituitary adenoma can secrete excessive GH to cause pathological changes in body tissues. GH proteoform changes would be associated with GH-related disease pathogenesis. PURPOSE: This study aimed to elucidate changes in GH proteoforms between GH-secreting pituitary adenomas and control pituitaries for the predictive diagnostics, targeted prevention, and personalization of medical services. METHODS: The isoelectric point (pI) and relative molecular mass (Mr) are two basic features of a proteoform that can be used to effectively array and detect proteoforms with two-dimensional gel electrophoresis (2DGE) and 2DGE-based western blot. GH proteoforms were characterized with liquid chromatography (LC) and mass spectrometry (MS). Phosphoproteomics, ubiquitinomics, acetylomics, and bioinformatics were used to analyze post-translational modifications (PTMs) of GH proteoforms in GH-secreting pituitary adenoma tissues and control pituitaries. RESULTS: Sixty-six 2D gel spots were found to contain hGH, including 46 spots (46 GH proteoforms) in GH-secreting pituitary adenomas and 35 spots (35 GH proteoforms) in control pituitaries. Further, 35 GH proteoforms in control pituitary tissues were matched with 35 of 46 GH proteoforms in GH-secreting pituitary adenoma tissues; and 11 GH proteoforms were presented in only GH-secreting pituitary adenoma tissues but not in control pituitary tissues. The matched 35 GH proteoforms showed quantitative changes in GH-secreting pituitary adenomas compared to the controls. The quantitative levels of those 46 GH proteoforms in GH-secreting pituitary adenomas were significantly different from those 35 GH proteoforms in control pituitaries. Meanwhile, different types of PTMs were identified among those GH proteoforms. Phosphoproteomics identified phosphorylation at residues Ser77, Ser132, Ser134, Thr174, and Ser176 in hGH. Ubiquitinomics identified ubiquitination at residue Lys96 in hGH. Acetylomics identified acetylation at reside Lys171 in hGH. Deamination was identified at residue Asn178 in hGH. CONCLUSION: These findings provide the first hGH proteoform pattern changes in GH-secreting pituitary adenoma tissues compared to control pituitary tissues, and the status of partial PTMs in hGH proteoforms. Those data provide in-depth insights into biological roles of hGH in GH-related diseases, and identify hGH proteoform pattern biomarkers for treatment of a GH-secreting pituitary adenoma in the context of 3P medicine -predictive diagnostics, targeted prevention, and personalization of medical services. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13167-021-00232-7.

SILAC quantitative proteomics analysis of ivermectin-related proteomic profiling and molecular network alterations in human ovarian cancer cells.

The antiparasitic agent ivermectin offers more promises to treat a diverse range of diseases. However, a comprehensive proteomic analysis of ivermectin-treated ovarian cancer (OC) cells has not been performed. This study sought to identify ivermectin-related proteomic profiling and molecular network alterations in human OC cells. Stable isotope labeling with amino acids in cell culture (SILAC)-based quantitative proteomics was used to study the human OC TOV-21G cells. After TOV-21G cells underwent 10 passages in SILAC-labeled growth media, TOV-21G cells were treated with 10 ml of 20 µmol/L ivermectin in cell growing medium for 24 h. The SILAC-labeled proteins were digested with trypsin; tryptic peptides were identified with mass spectrometry (MS). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was used to mine signaling pathway alterations with ivermectin-related proteins in TOV-21G cells. Gene ontology (GO) analysis was used to explore biological functions of ivermectin-related proteins, including biological processes (BPs), cellular components (CCs), and molecular functions (MFs). The protein-protein interaction network was analyzed with molecular complex detection (MCODE) to identify hub modules. In total, 4,447 proteins were identified in ivermectin-treated TOV-21G cells. KEGG analysis revealed 89 statistically significant signaling pathways. Interestingly, the clustering analysis of these pathways showed that ivermectin was involved in various cancer pathogenesis processes, including modulation of replication, RNA metabolism, and translational machinery. GO analysis revealed 69 statistically significant CCs, 87 MFs, and 62 BPs. Furthermore, MCODE analysis identified five hub modules, including 147 hub molecules. Those hub modules involved ribosomal proteins, RNA-binding proteins, cell-cycle progression-related proteins, proteasome subunits, and minichromosome maintenance proteins. These findings demonstrate that SILAC quantitative proteomics is an effective method to analyze ivermectin-treated cells, provide the first ivermectin-related proteomic profiling and molecular network alterations in human OC cells, and provide deeper insights into molecular mechanisms and functions of ivermectin to inhibit OC cells.

Quantitative Analysis of Ubiquitinated Proteins in Human Pituitary and Pituitary Adenoma Tissues.

Protein ubiquitination is an important post-translational modification that is associated with multiple diseases, including pituitary adenomas (PAs). Protein ubiquitination profiling in human pituitary and PAs remains unknown. Here, we performed the first ubiquitination analysis with an anti-ubiquitin antibody (specific to K-ε-GG)-based label-free quantitative proteomics method and bioinformatics to investigate protein ubiquitination profiling between PA and control tissues. A total of 158 ubiquitinated sites and 142 ubiquitinated peptides in 108 proteins were identified, and five ubiquitination motifs were found. KEGG pathway network analysis of 108 ubiquitinated proteins identified four statistically significant signaling pathways, including PI3K-AKT signaling pathway, hippo signaling pathway, ribosome, and nucleotide excision repair. R software Gene Ontology (GO) analysis of 108 ubiquitinated proteins revealed that protein ubiquitination was involved in multiple biological processes, cellular components, and molecule functions. The randomly selected ubiquitinated 14-3-3 zeta/delta protein was further analyzed with Western blot, and it was found that upregulated 14-3-3 zeta/delta protein in nonfunctional PAs might be derived from the significantly decreased level of its ubiquitination compared to control pituitaries, which indicated a contribution of 14-3-3 zeta/delta protein to pituitary tumorigenesis. These findings provided the first ubiquitinated proteomic profiling and ubiquitination-involved signaling pathway networks in human PAs. This study offers new scientific evidence and basic data to elucidate the biological functions of ubiquitination in PAs, insights into its novel molecular mechanisms of pituitary tumorigenesis, and discovery of novel biomarkers and therapeutic targets for effective treatment of PAs.

Prolactin Variants in Human Pituitaries and Pituitary Adenomas Identified With Two-Dimensional Gel Electrophoresis and Mass Spectrometry.

Human prolactin (hPRL) plays multiple roles in growth, metabolism, development, reproduction, and immunoregulation, which is an important protein synthesized in a pituitary. Two-dimensional gel electrophoresis (2DE) is an effective method in identity of protein variants for in-depth insight into functions of that protein. 2DE, 2DE-based PRL-immunoblot, mass spectrometry, and bioinformatics were used to analyze hPRL variants in human normal (control; n = 8) pituitaries and in five subtypes of pituitary adenomas [NF- (n = 3)-, FSH+ (n = 3)-, LH+ (n = 3)-, FSH+/LH+ (n = 3)-, and PRL+ (n = 3)-adenomas]. Six hPRL variants were identified with different isoelectric point (pI)-relative molecular mass (Mr ) distribution on a 2DE pattern, including variants V1 (pI 6.1; 26.0 kDa), V2 (pI 6.3; 26.4 kDa), V3 (pI 6.3; 27.9 kDa), V4 (pI 6.5; 26.1 kDa), V5 (pI 6.8; 25.9 kDa), and V6 (pI 6.7; 25.9 kDa). Compared to controls, except for variants V2-V6 in PRL-adenomas, V2 in FSH+-adenomas, and V3 in NF--adenomas, the other PRL variants were significantly downregulated in each subtype of pituitary adenomas. Moreover, the pattern of those six PRL variants was significantly different among five subtypes of pituitary adenomas relative to control pituitaries. Different hPRL variants might be involved in different types of PRL receptor-signaling pathways in a given condition. Those findings clearly revealed the existence of six hPRL variants in human pituitaries, and the pattern changes of six hPRL variants among different subtypes of pituitary adenomas, which provide novel clues to further study the functions, and mechanisms of action, of hPRL in human pituitary and in PRL-related diseases, and the potential clinical value in pituitary adenomas.

Nitroproteins in Human Astrocytomas Discovered by Gel Electrophoresis and Tandem Mass Spectrometry.

Protein tyrosine nitration is involved in the pathogenesis of highly fatal astrocytomas, a type of brain cancer. To understand the molecular mechanisms of astrocytomas and to discover new biomarkers/therapeutic targets, we sought to identify nitroproteins in human astrocytoma tissue. Anti-nitrotyrosine immunoreaction-positive proteins from a high-grade astrocytoma tissue were detected with two-dimensional gel electrophoresis (2DGE)-based nitrotyrosine immunoblots, and identified with liquid chromatography-tandem mass spectrometry (LC-MS/MS). Fifty-seven nitrotyrosine immunopositive protein spots were detected. A total of 870 proteins (nitrated and non-nitrated) in nitrotyrosine-immunopositive 2D gel spots were identified, and 18 nitroproteins and their 20 nitrotyrosine sites were identified with MS/MS analysis. These nitroproteins participate in multiple processes, including drug-resistance, signal transduction, cytoskeleton, transcription and translation, cell proliferation and apoptosis, immune response, phenotypic dedifferentiation, cell migration, and metastasis. Among those nitroproteins that might play a role in astrocytomas was nitro-sorcin, which is involved in drug resistance and metastasis and might play a role in the spread and treatment of an astrocytoma. Semiquantitative immune-based measurements of different sorcin expressions were found among different grades of astrocytomas relative to controls, and a semiquantitative increased nitration level in high-grade astrocytoma relative to control. Nitro-ß-tubulin functions in cytoskeleton and cell migration. Semiquantitative immunoreactivity of ß-tubulin showed increased expression among different grades of astrocytomas relative to controls and semiquantitatively increased nitration level in high-grade astrocytoma relative to control. Each nitroprotein was rationalized and related to the corresponding functional system to provide new insights into tyrosine nitration and its potential role in the pathogenesis of astrocytoma formation. Graphical Abstract ᅟ.

EPMA position paper in cancer: current overview and future perspectives.

At present, a radical shift in cancer treatment is occurring in terms of predictive, preventive, and personalized medicine (PPPM). Individual patients will participate in more aspects of their healthcare. During the development of PPPM, many rapid, specific, and sensitive new methods for earlier detection of cancer will result in more efficient management of the patient and hence a better quality of life. Coordination of the various activities among different healthcare professionals in primary, secondary, and tertiary care requires well-defined competencies, implementation of training and educational programs, sharing of data, and harmonized guidelines. In this position paper, the current knowledge to understand cancer predisposition and risk factors, the cellular biology of cancer, predictive markers and treatment outcome, the improvement in technologies in screening and diagnosis, and provision of better drug development solutions are discussed in the context of a better implementation of personalized medicine. Recognition of the major risk factors for cancer initiation is the key for preventive strategies (EPMA J. 4(1):6, 2013). Of interest, cancer predisposing syndromes in particular the monogenic subtypes that lead to cancer progression are well defined and one should focus on implementation strategies to identify individuals at risk to allow preventive measures and early screening/diagnosis. Implementation of such measures is disturbed by improper use of the data, with breach of data protection as one of the risks to be heavily controlled. Population screening requires in depth cost-benefit analysis to justify healthcare costs, and the parameters screened should provide information that allow an actionable and deliverable solution, for better healthcare provision.

EVpedia: a community web portal for extracellular vesicles research.

MOTIVATION: Extracellular vesicles (EVs) are spherical bilayered proteolipids, harboring various bioactive molecules. Due to the complexity of the vesicular nomenclatures and components, online searches for EV-related publications and vesicular components are currently challenging. RESULTS: We present an improved version of EVpedia, a public database for EVs research. This community web portal contains a database of publications and vesicular components, identification of orthologous vesicular components, bioinformatic tools and a personalized function. EVpedia includes 6879 publications, 172 080 vesicular components from 263 high-throughput datasets, and has been accessed more than 65 000 times from more than 750 cities. In addition, about 350 members from 73 international research groups have participated in developing EVpedia. This free web-based database might serve as a useful resource to stimulate the emerging field of EV research. AVAILABILITY AND IMPLEMENTATION: The web site was implemented in PHP, Java, MySQL and Apache, and is freely available at http://evpedia.info.

Mass spectrometry analysis of nitrotyrosine-containing proteins.

Oxidative stress plays important roles in a wide range of diseases such as cancer, inflammatory disease, neurodegenerative disorders, etc. Tyrosine nitration in a protein is a chemically stable oxidative modification, and a marker of oxidative injuries. Mass spectrometry (MS) is a key technique to identify nitrotyrosine-containing proteins and nitrotyrosine sites in endogenous and synthetic nitroproteins and nitropeptides. However, in vivo nitrotyrosine-containing proteins occur with extreme low-abundance to severely challenge the use of MS to identify in vivo nitroproteins and nitrotyrosine sites. A preferential enrichment of nitroproteins and/or nitropeptides is necessary before MS analysis. Current enrichment methods include immuno-affinity techniques, chemical derivation of the nitro group plus target isolations, followed with tandem mass spectrometry analysis. This article reviews the MS techniques and pertinent before-MS enrichment techniques for the identification of nitrotyrosine-containing proteins. This article reviews future trends in the field of nitroproteomics, including quantitative nitroproteomics, systems biological networks of nitroproteins, and structural biology study of tyrosine nitration to completely clarify the biological functions of tyrosine nitration.

Heterogeneity analysis of the proteomes in clinically nonfunctional pituitary adenomas.

BACKGROUND: Clinically nonfunctional pituitary adenomas (NFPAs) without any clinical elevation of hormone and with a difficulty in its early-stage diagnosis are highly heterogeneous with different hormone expressions in NFPA tissues, including luteinizing hormone (LH)-positive, follicle-stimulating hormone (FSH)-positive, LH/FSH-positive, and negative (NF). Elucidation of molecular mechanisms and discovery of biomarkers common and specific to those different subtypes of NFPAs will benefit NFPA patients in early-stage diagnosis and individualized treatment. METHODS: Two-dimensional gel electrophoresis (2DGE) and PDQuest image analyses were used to compare proteomes of different NFPA subtypes (NF-, LH-, FSH-, and LH/FSH-positive) relative to control pituitaries (Con). Differentially expressed proteins (DEPs) were characterized with mass spectrometry (MS). Each set of DEPs in four NFPA subtypes was evaluated with overlap analysis and signaling pathway network analysis with comparison to determine any DEP and pathway network that are common and specific to each NFPA subtype. RESULTS: A total of 93 differential protein-spots were determined with comparison of each NFPA type (NF-, LH-, FSH-, and LH/FSH-positive) versus control pituitaries. A total of 76 protein-spots were MS-identified (59 DEPs in NF vs. Con; 65 DEPs in LH vs. Con; 63 DEPs in FSH vs. Con; and 55 DEPs in LH/FSH vs. Con). A set of DEPs and pathway network data were common and specific to each NFPA subtype. Four important common pathway systems included MAPK-signaling abnormality, oxidative stress, mitochondrial dysfunction, and cell-cycle dysregulation. However, these pathway systems were, in fact, different among four NFPA subtypes with different protein-expression levels of most of nodes, different protein profiles, and different pathway network profiles. CONCLUSIONS: These result data demonstrate that common and specific DEPs and pathway networks exist in four NFPA subtypes, and clarify proteome heterogeneity of four NFPA subtypes. Those findings will help to elucidate molecular mechanisms of NFPAs, and discover protein biomarkers to effectively manage NFPA patients towards personalized medicine.

Identification of the proteomic variations of invasive relative to non-invasive non-functional pituitary adenomas.

The incomplete surgery section of invasive non-functional pituitary adenomas (NFPAs) carries the increased risks of complications and requires adjuvant radiotherapy and medications. It is necessary to clarify the molecular mechanisms and markers of invasiveness to guide the management of NFPA patients. The study aimed to proteomic variations of invasive and non-invasive NFPAs and sought the protein markers for invasive NFPAs. Invasive (n = 4) and non-invasive (n = 4) NFPA tissues were analyzed (n = 3-5/each tissue) with 2DE and PDQuest software. Twenty-four high-resolution 2DE gels were quantitatively compared to determine differentially expressed proteins (DEPs) between invasive and non-invasive NFPAs. Approximately 1200 protein spots were detected in each 2DE map, and 103 differential spots (64 upregulated and 39 downregulated) were identified. Among those 103 differential spots, 57 DEPs (30 upregulated and 27 downregulated) were characterized with peptide mass fingerprint and MS/MS. Gene-ontology (GO) and ingenuity pathway analyses of those DEPs revealed pathway networks including mitochondrial dysfunction, oxidative stress, mitogen-activated protein kinase signaling abnormality, TR/RXR activation, proteolysis abnormality, ketogenesis and ketolysis, cyclin-dependent kinase C signaling abnormality, and amyloid processing that were significantly associated with invasive characteristics of invasive NFPA. Those data demonstrate that proteomic variations exist between invasive and non-invasive NFPAs. 2DE-based comparative proteomics is an effective approach to identify proteomic variations and pathway network variations. Those findings will serve as a basis to understand the molecular mechanisms of invasive NFPAs and to discover protein markers to effectively manage patients with invasive NFPAs.

Deimination of linker histones links neutrophil extracellular trap release with autoantibodies in systemic autoimmunity.

Autoantibodies to nuclear antigens arise in human autoimmune diseases, but a unifying pathogenetic mechanism remains elusive. Recently we reported that exposure of neutrophils to inflammatory conditions induces the citrullination of core histones by peptidylarginine deiminase 4 (PAD4) and that patients with autoimmune disorders produce autoantibodies that recognize such citrullinated histones. Here we identify histone H1 as an additional substrate of PAD4, localize H1 within neutrophil extracellular traps, and detect autoantibodies to citrullinated H1 in 6% of sera from patients with systemic lupus erythematosus and Sjögren's syndrome. No preference for deiminated H1 was observed in healthy control sera and sera from patients with scleroderma or rheumatoid arthritis. We map binding to the winged helix of H1 and determine that citrulline 53 represents a key determinant of the autoantibody epitope. In addition, we quantitate RNA for H1 histone subtypes in mature human neutrophils and identify citrulline residues by liquid chromatography and tandem mass spectrometry. Our results indicate that deimination of linker histones generates new autoantibody epitopes with enhanced potential for stimulating autoreactive human B cells.-Dwivedi, N., Neeli, I., Schall, N., Wan, H., Desiderio, D. M., Csernok, E., Thompson, P. R., Dali, H., Briand, J.-P., Muller, S., Radic, M. Deimination of linker histones links neutrophil extracellular trap release with autoantibodies in systemic autoimmunity.

EVpedia: an integrated database of high-throughput data for systemic analyses of extracellular vesicles.

Secretion of extracellular vesicles is a general cellular activity that spans the range from simple unicellular organisms (e.g. archaea; Gram-positive and Gram-negative bacteria) to complex multicellular ones, suggesting that this extracellular vesicle-mediated communication is evolutionarily conserved. Extracellular vesicles are spherical bilayered proteolipids with a mean diameter of 20-1,000 nm, which are known to contain various bioactive molecules including proteins, lipids, and nucleic acids. Here, we present EVpedia, which is an integrated database of high-throughput datasets from prokaryotic and eukaryotic extracellular vesicles. EVpedia provides high-throughput datasets of vesicular components (proteins, mRNAs, miRNAs, and lipids) present on prokaryotic, non-mammalian eukaryotic, and mammalian extracellular vesicles. In addition, EVpedia also provides an array of tools, such as the search and browse of vesicular components, Gene Ontology enrichment analysis, network analysis of vesicular proteins and mRNAs, and a comparison of vesicular datasets by ortholog identification. Moreover, publications on extracellular vesicle studies are listed in the database. This free web-based database of EVpedia (http://evpedia.info) might serve as a fundamental repository to stimulate the advancement of extracellular vesicle studies and to elucidate the novel functions of these complex extracellular organelles.

Pituitary adenoma nitroproteomics: current status and perspectives.

Oxidative stress is extensively associated with tumorigenesis. A series of studies on stable tyrosine nitration as a marker of oxidative damage were performed in human pituitary and adenoma. This paper reviews published research on the mass spectrometry characteristics of nitropeptides and nitroproteomics of pituitary controls and adenomas. The methodology used for nitroproteomics, the current status of human pituitary nitroproteomics studies, and the future perspectives are reviewed. Enrichment of those low-abundance endogenous nitroproteins from human tissues or body fluid samples is the first important step for nitroproteomics studies. Mass spectrometry is the essential approach to determine the amino acid sequence and locate the nitrotyrosine sites. Bioinformatics analyses, including protein domain and motif analyses, are needed to locate the nitrotyrosine site within the corresponding protein domains/motifs. Systems biology techniques, including pathway analysis, are necessary to discover signaling pathway networks involving nitroproteins from the systematically global point of view. Future quantitative nitroproteomics will discover pituitary adenoma-specific nitroprotein(s). Structural biology techniques such as X-ray crystallography analysis will solidly clarify the effects of tyrosine nitration on structure and functions of a protein. Those studies will eventually address the mechanisms and biological functions of tyrosine nitration in pituitary tumorigenesis and will discover nitroprotein biomarkers for pituitary adenomas and targets for drug design for pituitary adenoma therapy.

LC-MS/MS identification of the one-carbon cycle metabolites in human plasma.

The one-carbon cycle is composed of four major biologically important molecules: methionine (L-Met), S-adenosylmethionine (SAM), S-adenosylhomocysteine (SAH), and homocysteine (Hcy). In addition to these key metabolites, there are multiple enzymes, vitamins, and cofactors that play essential roles in the cascade of the biochemical reactions that convert one metabolite into another in the cycle. Simultaneous quantitative measurement of four major metabolites can be used to detect possible aberrations in this vital cycle. Abnormalities in the one-carbon cycle might lead to hyper- or hypomethylation, homocystinemia, liver dysfunction, and accumulation of white-matter hyperintensities in the human brain. Previously published methods describe evaluation of several components of the one-carbon cycle, but none to our knowledge demonstrated simultaneous measurement of all four key molecules (L-Met, SAM, SAH, and Hcy). We describe a novel analytical method suitable for simultaneous identification and quantification of L-Met, SAM, SAH, and Hcy with LC-MS/MS. Moreover, we tested this method to identify these metabolites in human plasma collected from patients with multiple sclerosis and healthy individuals. In a pilot feasibility study, our results indicate that patients with multiple sclerosis showed abnormalities in the one-carbon cycle.

Molecular phenotyping of immune cells from young NOD mice reveals abnormal metabolic pathways in the early induction phase of autoimmune diabetes.

Islet leukocytic infiltration (insulitis) is first obvious at around 4 weeks of age in the NOD mouse--a model for human type 1 diabetes (T1D). The molecular events that lead to insulitis and initiate autoimmune diabetes are poorly understood. Since TID is caused by numerous genes, we hypothesized that multiple molecular pathways are altered and interact to initiate this disease. We evaluated the molecular phenotype (mRNA and protein expression) and molecular networks of ex vivo unfractionated spleen leukocytes from 2 and 4 week-old NOD mice in comparison to two control strains. Analysis of the global gene expression profiles and hierarchical clustering revealed that the majority (~90%) of the differentially expressed genes in NOD mice were repressed. Furthermore, analysis using a modern suite of multiple bioinformatics approaches identified abnormal molecular pathways that can be divided broadly into 2 categories: metabolic pathways, which were predominant at 2 weeks, and immune response pathways, which were predominant at 4 weeks. Network analysis by Ingenuity pathway analysis identified key genes/molecules that may play a role in regulating these pathways. These included five that were common to both ages (TNF, HNF4A, IL15, Progesterone, and YWHAZ), and others that were unique to 2 weeks (e.g. MYC/MYCN, TGFB1, and IL2) and to 4 weeks (e.g. IFNG, beta-estradiol, p53, NFKB, AKT, PRKCA, IL12, and HLA-C). Based on the literature, genes that may play a role in regulating metabolic pathways at 2 weeks include Myc and HNF4A, and at 4 weeks, beta-estradiol, p53, Akt, HNF4A and AR. Our data suggest that abnormalities in regulation of metabolic pathways in the immune cells of young NOD mice lead to abnormalities in the immune response pathways and as such may play a role in the initiation of autoimmune diabetes. Thus, targeting metabolism may provide novel approaches to preventing and/or treating autoimmune diabetes.