Email updates

Keep up to date with the latest news and content from BMC Medical Genomics and BioMed Central.

Open Access Research article

Literature-based discovery of diabetes- and ROS-related targets

Junguk Hur12, Kelli A Sullivan2, Adam D Schuyler4, Yu Hong2, Manjusha Pande23, David J States5, H V Jagadish13 and Eva L Feldman123*

Author Affiliations

1 Bioinformatics Program, University of Michigan, Ann Arbor, MI 48109, USA

2 Department of Neurology, University of Michigan, Ann Arbor, MI 48109, USA

3 National Center for Integrative Biomedical Informatics, University of Michigan, Ann Arbor, MI 48109, USA

4 Department of Molecular, Microbial and Structural Biology, University of Connecticut Health Center, Farmington, CT 06030, USA

5 School of Health Information Science, University of Texas, Houston, TX 77030, USA

For all author emails, please log on.

BMC Medical Genomics 2010, 3:49  doi:10.1186/1755-8794-3-49

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1755-8794/3/49


Received:6 May 2010
Accepted:27 October 2010
Published:27 October 2010

© 2010 Hur et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Reactive oxygen species (ROS) are known mediators of cellular damage in multiple diseases including diabetic complications. Despite its importance, no comprehensive database is currently available for the genes associated with ROS.

Methods

We present ROS- and diabetes-related targets (genes/proteins) collected from the biomedical literature through a text mining technology. A web-based literature mining tool, SciMiner, was applied to 1,154 biomedical papers indexed with diabetes and ROS by PubMed to identify relevant targets. Over-represented targets in the ROS-diabetes literature were obtained through comparisons against randomly selected literature. The expression levels of nine genes, selected from the top ranked ROS-diabetes set, were measured in the dorsal root ganglia (DRG) of diabetic and non-diabetic DBA/2J mice in order to evaluate the biological relevance of literature-derived targets in the pathogenesis of diabetic neuropathy.

Results

SciMiner identified 1,026 ROS- and diabetes-related targets from the 1,154 biomedical papers (http://jdrf.neurology.med.umich.edu/ROSDiabetes/ webcite). Fifty-three targets were significantly over-represented in the ROS-diabetes literature compared to randomly selected literature. These over-represented targets included well-known members of the oxidative stress response including catalase, the NADPH oxidase family, and the superoxide dismutase family of proteins. Eight of the nine selected genes exhibited significant differential expression between diabetic and non-diabetic mice. For six genes, the direction of expression change in diabetes paralleled enhanced oxidative stress in the DRG.

Conclusions

Literature mining compiled ROS-diabetes related targets from the biomedical literature and led us to evaluate the biological relevance of selected targets in the pathogenesis of diabetic neuropathy.

Background

Diabetes is a metabolic disease in which the body does not produce or properly respond to insulin, a hormone required to convert carbohydrates into energy for daily life. According to the American Diabetes Association, 23.6 million children and adults, approximately 7.8% of the population in the United States, have diabetes [1]. The cost of diabetes in 2007 was estimated to be $174 billion [1]. The micro- and macro-vascular complications of diabetes are the most common causes of renal failure, blindness and amputations leading to significant mortality, morbidity and poor quality of life; however, incomplete understanding of the causes of diabetic complications hinders the development of mechanism-based therapies.

In vivo and in vitro experiments implicate a number of enzymatic and non-enzymatic metabolic pathways in the initiation and progression of diabetic complications [2] including: (1) increased polyol pathway activity leading to sorbitol and fructose accumulation, NAD(P)-redox imbalances and changes in signal transduction; (2) non-enzymatic glycation of proteins yielding "advanced glycation end-products" (AGEs); (3) activation of protein kinase C (PKC), initiating a cascade of intracellular stress responses; and (4) increased hexosamine pathway flux [2,3]. Only recently has a link among these pathways been established that provides a unified mechanism of tissue damage. Each of these pathways directly and indirectly leads to overproduction of reactive oxygen species (ROS) [2,3].

ROS are highly reactive ions or small molecules including oxygen ions, free radicals and peroxides, formed as natural byproducts of cellular energy metabolism. ROS are implicated in multiple cellular pathways such as mitogen-activated protein kinase (MAPK) signaling, c-Jun amino-terminal kinase (JNK), cell proliferation and apoptosis [4-6]. Due to the highly reactive properties of ROS, excessive ROS may cause significant damage to proteins, DNA, RNA and lipids. All cells express enzymes capable of neutralizing ROS. In addition to the maintenance of antioxidant systems such as glutathione and thioredoxins, primary sensory neurons express two main detoxifying enzymes: superoxide dismutase (SOD) [7] and catalase [8]. SOD converts superoxide (O2-) to H2O2, which is reduced to H2O by glutathione and catalase [8]. SOD1 is the main form of SOD in the cytoplasm; SOD2 is located within the mitochondria. In neurons, SOD1 activity represents approximately 90% of total SOD activity and SOD2 approximately 10% [9]. Under diabetic conditions, this protective mechanism is overwhelmed due to the substantial increase in ROS, leading to cellular damage and dysfunction [10].

The idea that increased ROS and oxidative stress contribute to the pathogenesis of diabetic complications has led scientists to investigate different oxidative stress pathways [7,11]. Inhibition of ROS or maintenance of euglycemia restores metabolic and vascular imbalances and blocks both the initiation and progression of complications [12,13]. Despite the significant implications and extensive research into the role of ROS in diabetes, no comprehensive database regarding ROS-related genes or proteins is currently available.

In the present study, a comprehensive list of ROS- and diabetes-related targets (genes/proteins) was compiled from the biomedical literature through text mining technology. SciMiner, a web-based literature mining tool [14], was used to retrieve and process documents and identify targets from the text. SciMiner provides a convenient web-based platform for target-identification within the biomedical literature, similar to other tools including EBIMed [15], ALI BABA [16], and PolySearch [17]; however, SciMiner is unique in that it searches full text documents, supports free-text PubMed query style, and allows the comparison of target lists from multiple queries.

The ROS-diabetes targets collected by SciMiner were further tested against randomly selected non-ROS-diabetes literature to identify targets that are significantly over-represented in the ROS-diabetes literature. Functional enrichment analyses were performed on these targets to identify significantly over-represented biological functions in terms of Gene Ontology (GO) terms and pathways.

In order to confirm the biological relevance of the over-represented ROS-diabetes targets, the gene expression levels of nine selected targets were measured in dorsal root ganglia (DRG) from mice with and without diabetes. DRG contain primary sensory neurons that relay information from the periphery to the central nervous system (CNS) [7,10,18]. Unlike the CNS, DRG are not protected by a blood-nerve barrier, and are consequently vulnerable to metabolic and toxic injury [19]. We hypothesize that differential expression of identified targets in DRG would confirm their involvement in the pathogenesis of diabetic neuropathy.

Methods

Defining ROS-diabetes literature

To retrieve the list of biomedical literature associated with ROS and diabetes, PubMed was queried using ("Reactive Oxygen Species"[MeSH] AND "Diabetes Mellitus"[MeSH]). This query yielded 1,154 articles as of April 27, 2009. SciMiner, a web-based literature mining tool [14], was used to retrieve and process the abstracts and available full text documents to identify targets (full text documents were available for approximately 40% of the 1,154 articles). SciMiner-identified targets, reported in the form of HGNC [HUGO (Human Genome Organization) Gene Nomenclature Committee] genes, were confirmed by manual review of the text.

Comparison with human curated data (NCBI Gene2PubMed)

The NCBI Gene database provides links between Gene and PubMed. The links are the result of (1) manual curation within the NCBI via literature analysis as part of generating a Gene record, (2) integration of information from other public databases, and (3) GeneRIF (Gene Reference Into Function) in which human experts provide a brief summary of gene functions and make the connections between citation (PubMed) and Gene databases. For the 1,154 ROS-diabetes articles, gene-paper associations were retrieved from the NCBI Gene database. Non-human genes were mapped to homologous human genes through the NCBI HomoloGene database. The retrieved genes were compared against the SciMiner derived targets. Any genes missed by SciMiner were added to the ROS-diabetes target set.

Protein-protein interactions among ROS-diabetes targets

To indirectly examine the association of literature derived targets (by SciMiner and NCBI Gene2PubMed) with ROS and diabetes, protein-protein interactions (PPIs) among the targets were surveyed. This was based on an assumption that targets are more likely to have PPIs with each other if they are truly associated within the same biological functions/pathways. A PPI network of the ROS-diabetes targets was generated using the Michigan Molecular Interactions (MiMI, http://mimi.ncibi.org/ webcite) database [20] and compared against 100 PPI networks of randomly drawn sets (the same number of the ROS-diabetes target set) from HUGO. A standard Z-test and one sample T-test were used to calculate the statistical significance of the ROS-diabetes PPI network with respect to the random PPI networks.

Functional enrichment analysis

Literature derived ROS-diabetes targets (by SciMiner and NCBI Gene2PubMed) were subject to functional enrichment analyses to identify significantly over-represented biological functions in terms of Gene Ontology [21], pathways (Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/ webcite) [22] and Reactome http://www.reactome.org/ webcite[23]). Fisher's exact test [24] was used to calculate the statistical significance of these biological functions with Benjamini-Hochberg (BH) adjusted p-value < 0.05 [25] as the cut-off.

Over-represented ROS-diabetes targets

Defining background corpora

To identify a subset of targets that are highly over-represented within the ROS-diabetes targets, the frequency of each target (defined as the number of documents in which the target was identified divided by the number of total documents in the query) was compared against the frequencies in randomly selected background corpora. Depending on how the background set is defined, over-represented targets may vary widely; therefore, to maintain the background corpora close to the ROS and diabetes context, documents were selected from the same journal, volume, and issue of the 1,154 ROS-diabetes documents, but were NOT indexed with "Reactive Oxygen Species"[MeSH] nor "Diabetes Mellitus"[MeSH]. For example, one of the ROS-diabetes articles (PMID: 18227068), was published in the Journal of Biological Chemistry, Volume 283, Issue 16. This issue contained 85 papers, 78 of which were not indexed with either "Reactive Oxygen Species"[MeSH] or "Diabetes Mellitus"[MeSH] indexed. One of these 78 papers was randomly selected as a background document. Three sets of 1,154 documents were selected using this approach and processed using SciMiner. Identified targets were confirmed by manual review for accuracy.

Identifying significantly over-represented targets

ROS-diabetes targets were tested for over-representation against targets identified from the three background sets. Fisher's exact test was used to determine if the frequency of each target in the ROS-diabetes target set was significantly different from that of the background sets. Any targets with a BH adjusted p-value < 0.05 in at least two of the three comparisons were deemed to be an over-represented ROS-diabetes target. Functional enrichment analyses were performed on these over-represented ROS-diabetes targets as described above.

Selecting targets for real-time RT-PCR

A subset of targets were selected for RT-PCR from the top 10 over-represented ROS-diabetes targets excluding insulin and NADPH oxidase 5 (NOX5), which does not have a mouse ortholog. Nitric oxide synthase 1 (NOS1), the main generator of nitric oxide, ranked at the 15th position and was additionally selected for inclusion in the test set.

Differential gene expression by real-time RT-PCR

Mice

DBA/2J mice were purchased from the Jackson Laboratory (Bar Harbor, ME). Mice were housed in a pathogen-free environment and cared for following the University of Michigan Committee on the Care and Use of Animals guidelines. Mice were fed AIN76A chow (Research Diets, New Brunswick, NJ). Male mice were used for this study.

Induction of diabetes

Two treatment groups were defined: control (n = 4) and diabetic (n = 4). Diabetes was induced at 13 weeks of age by low-dose streptozotocin (STZ) injections, 50 mg/kg/day for five consecutive days. All diabetic mice received LinBit sustained release insulin implants (LinShin, Toronto, Canada) at 8 weeks post-STZ treatment. Insulin implants were replaced every 4 weeks, at 12 and 16 weeks post-STZ treatment. At 20 weeks post-STZ treatment, mice were euthanized by sodium pentobarbital overdose and DRG were harvested as previously described [26].

Real-time RT-PCR

The gene expression of the selected nine literature-derived ROS-diabetes targets in DRG was measured using real-time RT-PCR in duplicate. The amount of mRNA isolated from each DRG was normalized to an endogenous reference [Tbp: TATA box binding protein; Δ cycle threshold (CT)].

Results

Identification of ROS-diabetes targets

A total of 1,021 unique targets were identified by SciMiner from the 1,154 ROS-diabetes papers defined by the query of ("Reactive Oxygen Species"[MeSH] AND "Diabetes Mellitus"[MeSH]) and confirmed by manual review. Table 1 contains the top 10 most frequently mentioned targets in the ROS-diabetes papers. Insulin was the most frequently mentioned target, followed by superoxide dismutase 1 and catalase.

Table 1. Top 10 most frequent ROS-diabetes targets

The NCBI Gene2PubMed database, containing expert-curated associations between the NCBI Gene and PubMed databases, revealed 90 unique genes associated with the 1,154 ROS-diabetes papers (Additional File 1). SciMiner identified 85 out of these 90 targets, indicating a 94% recall rate. Five targets missed by SciMiner were added to the initial ROS-diabetes target set to result in 1,026 unique targets (Additional File 2).

Additional file 1. The list of 90 genes from the NCBI Gene2PubMed database for the ROS-Diabetes literature (1,154 papers).

Format: XLS Size: 35KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 2. The list of 1,026 ROS-Diabetes targets.

Format: XLS Size: 229KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

PPI network of the ROS-diabetes targets

The PPI network among the ROS-diabetes targets was evaluated using MiMI interaction data. This was based on the assumption that targets commonly related to a certain topic are more likely to have frequent interactions with each other. One hundred PPI networks were generated for comparison using the same number of genes (1,026) randomly selected from the complete HUGO gene set (25,254). The PPI network of the ROS-diabetes targets was significantly different from the randomly generated networks indicating their strong association with the topic "ROS and Diabetes". Table 2 demonstrates that the mean number of targets with any PPI interaction in the randomly generated target sets was 528.9 (approximately 52% of 1,026 targets), while the number of targets with any PPI interaction in the ROS-diabetes target was 983 (96%). The number of targets interacting with each other was also significantly different between the random networks (mean = 155.4) and the ROS-diabetes network (mean = 879). Figure 1 illustrates the distributions of these measurements from the 100 random networks with the ROS-diabetes set depicted as a red vertical line. It is obvious that the PPI network of the ROS-diabetes targets is significantly different from the random networks.

Table 2. Summary of 100 randomly generated PPI networks

thumbnailFigure 1. Histograms of randomly generated PPI networks. The histograms illustrate the distributions of 100 randomly generated networks, while the red line indicates the ROS-diabetes targets. The network of the ROS-diabetes targets is significantly different from the 100 randomly generated networks, indicating the overlap of ROS-diabetes targets with respect to the topic "Reactive Oxygen Species and Diabetes".

Functional enrichment analyses of the ROS-diabetes targets

Functional enrichment analyses of the 1,026 ROS-diabetes targets were performed to identify over-represented biological functions of the ROS-diabetes targets. After Benjamini-Hochberg correction, a total of 189 molecular functions, 450 biological processes, 73 cellular components and 341 pathways were significantly enriched in the ROS-diabetes targets when compared against all the HUGO genes (see Additional Files 3, 4, 5 and 6 for the full lists). Table 3 lists the top 3 most over-represented GO terms and pathways ranked by p-values of Fisher's exact test: e.g., apoptosis, oxidoreductase activity and insulin signaling pathway.

Additional file 3. The enriched Molecular Functions Gene Ontology Terms in the 1,026 ROS-Diabetes targets.

Format: XLS Size: 91KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 4. The enriched Biological Processes Gene Ontology Terms in the 1,026 ROS-Diabetes targets.

Format: XLS Size: 95KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 5. The enriched Cellular Components Gene Ontology Terms in the 1,026 ROS-Diabetes targets.

Format: XLS Size: 61KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 6. The enriched pathways in the 1,026 ROS-Diabetes targets.

Format: XLS Size: 104KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 3. Enriched functions of 1,026 ROS-diabetes targets

Identification of over-represented ROS-diabetes targets

To identify the ROS-diabetes targets highly over-represented in ROS-diabetes literature, three sets of background corpora of the same size (n = 1,154 documents) were generated using the same journal, volume and issue approach. The overlap among the three background sets in terms of documents and identified targets are illustrated in Figure 2. Approximately 90% of the selected background documents were unique to the individual set, while 50% of the identified targets were identified in at least one of the three background document sets. The frequencies of the identified targets were compared among the background sets for significant differences. None of the targets had a BH adjusted p-value < 0.05, indicating no significant difference among the targets from the three different background sets (See Additional File 7).

thumbnailFigure 2. Venn diagrams of document compositions and identified targets of the randomly generated background sets. Approximately 90% of the selected background documents were unique to individual set (A), while 50% of the identified targets were identified in at least one of the three background document sets (B).

Additional file 7. Comparisons of target frequencies among three background sets.

Format: XLS Size: 22KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Comparisons of the ROS-diabetes targets against these background sets revealed 53 highly over-represented ROS-diabetes targets as listed in Table 4. These 53 targets were significant (p-value < 0.05) against all three background sets and significant following Benjamini-Hochberg multiple testing correction (BH adjusted p-value < 0.05) against at least two of the three background sets. SOD1 was the most over-represented in the ROS-diabetes targets.

Table 4. 53 targets over-represented in ROS-diabetes literature

Functional enrichment analyses of the over-represented ROS-diabetes targets

Functional enrichment analyses of the 53 ROS-diabetes targets were performed to identify over-represented biological functions. Following Benjamini-Hochberg correction, a total of 65 molecular functions, 209 biological processes, 26 cellular components and 108 pathways were significantly over-represented when compared against all the HUGO genes (see Additional Files 8, 9, 10 and 11 for the full lists). Table 5 shows the top 3 most significantly over-represented GO terms and pathways ranked by p-values of Fisher's exact test. GO terms related to oxidative stress such as "superoxide metabolic process", "superoxide release", "electron carrier activity" and "mitochondrion" were highly over-represented in the 53 ROS-diabetes targets.

Additional file 8. The enriched Molecular Functions Gene Ontology Terms in the Over-represented 53 ROS-Diabetes targets.

Format: XLS Size: 46KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 9. The enriched Biological Processes Gene Ontology Terms in the Over-represented 53 ROS-Diabetes targets.

Format: XLS Size: 95KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 10. The enriched Cellular Components Gene Ontology Terms in the Over-represented 53 ROS-Diabetes targets.

Format: XLS Size: 66KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 11. The enriched pathways in the Over-represented 53 ROS-Diabetes targets.

Format: XLS Size: 75KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 5. Enriched functions of the 53 over-represented targets in diabetes

Gene expression change in diabetes

Two groups of DBA/2J mice exhibited significantly different levels of glycosylated hemoglobin (%GHb). The mean ± SEM were 6.2 ± 0.3 for the non-diabetic control group and for 14.0 ± 0.8 for the diabetic group (p-value < 0.001), indicative of prolonged hyperglycemia in the diabetic group [26]. DRG were harvested from these animals for gene expression assays. Nine genes were selected from the top ranked ROS-diabetes targets: superoxide dismutase 1 (Sod1), catalase (Cat), xanthine dehydrogenase (Xdh), protein kinase C alpha (Prkca), neutrophil cytosolic factor 1 (Ncf1), nitric oxide synthase 3 (Nos3), superoxide dismutase 2 (Sod2), cytochrome b-245 alpha (Cyba), and nitric oxide synthase 1 (Nos1). Eight genes exhibited differential expression between diabetic and non-diabetic mice (p-value < 0.05) as shown in Figure 3. Cat, Sod1, Sod2, Prkca, and Nos1 expression levels were decreased, while Ncf1, Xdh, and Cyba expression levels were increased in diabetes.

thumbnailFigure 3. Gene expression levels of selected ROS-diabetes targets in DRG examined by real-time RT-PCR. Expression levels are relative to Tbp, an internal control (error bar = SEM) (*, p < 0.05; **, p < 0.01; ***, p < 0.001). Eight (Cat, Sod1, Ncf1, Xdh, Sod2, Cyba, Prkca, and Nos1) out of the nine selected ROS-diabetes genes were significantly regulated by diabetes.

Discussion

Reactive oxygen species (ROS) are products of normal energy metabolism and play important roles in many other biological processes such as the immune response and signaling cascades [4-6]. As mediators of cellular damage, ROS are implicated in pathogenesis of multiple diseases including diabetic complications [27-30]. With the aid of literature mining technology, we collected 1,026 possible ROS-related targets from a set of biomedical literature indexed with both ROS and diabetes.

Fifty-three targets were significantly over-represented in the ROS-diabetes papers when compared against three background sets. Depending on how the background set is defined, the over-represented targets may vary widely. An ideal background set would be the entire PubMed set; however, this is not possible due to limited access to full texts and intense data processing. An alternative method would be to use only abstracts in PubMed, but this may not fully represent the literature. Using only the abstracts, our target identification method resulted in 21 (39%) of the 53 key ROS-diabetes targets (Additional File 12), suggesting the benefit of rich information in full text documents. In the present study, background documents were randomly selected from the same journal, volume, and issue of the 1,154 ROS-diabetes documents, which were not indexed with "Reactive Oxygen Species"[MeSH] nor "Diabetes Mellitus"[MeSH]. This approach maintained the background corpora not far from the ROS and diabetes context.

Additional file 12. The Key 53 ROS-Diabetes Targets Identifiable Using Only the Abstracts.

Format: XLS Size: 23KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The gene expression levels of nine targets selected from the 53 over-represented ROS-diabetes targets were measured in diabetic and non-diabetic DRG. Our laboratory is particularly interested in deciphering the underlying mechanisms of diabetic neuropathy, a major complication of diabetes. Data published by our laboratory both in vitro and in vivo confirm the negative impact of oxidative stress in complication-prone neuron tissues like DRG [7,10,18,31]. In an effort to obtain diabetic neuropathy specific targets, SciMiner was employed to further analyze a subset of the ROS-diabetes papers (data not shown). Nerve growth factor (NGF) was identified as the most over-represented target in this subset when compared to the full ROS-diabetes set; however, NGF did not have statistical significance (BH adjusted p-value = 0.06). The relatively small numbers of papers and associated targets may have contributed to this non-significance. Therefore, the candidate targets for gene expression validation were selected from among the 53 over-represented ROS-diabetes targets derived from the full ROS-diabetes corpus.

Among the tested genes, the expression levels of Cat, Sod1, Sod2, Prkca, and Nos1 were decreased, while the expression levels of Ncf1, Xdh, and Cyba were increased under diabetic conditions. Cat, Sod1, and Sod2 are responsible for protecting cells from oxidative stress by destroying superoxides and hydrogen peroxides [8-11]. Decreased expression of these genes may result in oxidative stress [32]. Increased expression of Cyba and Ncf1, subunits of superoxide-generating nicotinamide adenine dinucleotide phosphate (NADPH) oxidase complex [30], also supports enhanced oxidative stress. Xdh and its inter-convertible form, Xanthine oxidase (Xod), showed increased activity in various rat tissues under oxidative stress conditions with diabetes [33], and also showed increased expression in diabetic DRG in the current study.

Unlike the above concordant genes, protein kinase C and nitric oxide synthases did not exhibit predicted expression changes in diabetes. Protein kinase C activates NADPH oxidase, further promoting oxidative stress in the cell [34,35]. Decreased expression of Prkca in our diabetic DRG is not parallel with expression levels of other enzymes expected to increase oxidative stress. Between the two nitric oxide synthases tested in the present study, Nos1 (neuronal) expression was significantly decreased (p-value < 0.001) in diabetes, while Nos3 (endothelial) expression was not significant (p-value = 0.06). The neuronal Nos1 is expected to play a major role in producing nitric oxide, another type of highly reactive free radical. Thus, with some exceptions, the majority of the differentially expressed genes in DRG show parallel results to the known activities of these targets in diabetes, suggesting enhanced oxidative stress in the diabetic DRG.

Assessment of antioxidant enzyme expression in diabetes has yielded a variety of results [36-40] depending upon the duration of diabetes, the tissue studied and other factors. In diabetic mice and rats, it is commonly reported that superoxide dismutases are down-regulated [37-40], where data regarding catalase are variable [36,40]. PKC is activated in diabetes, but most papers that examined mRNA demonstrated that its expression is largely unchanged [41].

Among the 53 over-represented ROS-diabetes targets, SOD1 was the most over-represented and was differentially expressed under diabetic and non-diabetic conditions. To the best of our knowledge, no published study has investigated the role of SOD1 in the onset and/or progression of diabetic neuropathy. Mutations of SOD1 have long been associated with the inherited form of amyotrophic lateral sclerosis (ALS) [42] and the theory of oxidative stress-based aging [43]. Early reports indicate that knockout of the SOD1 gene does not affect nervous system development [44], although recovery following injury is slow and incomplete [45,46]. With respect to diabetes, SOD1 KO accelerates the development of diabetic nephropathy [47] and cataract formation [48]. Thus, examining the SOD1 KO mouse as a model of diabetic neuropathy would be a reasonable follow-up study.

One limitation of the current approach using literature mining technology is incorrect or missed identification of the mentioned targets within the literature. Based on a performance evaluation using a standard text set BioCreAtIvE (Critical Assessment of Information Extraction systems in Biology) version 2 [49], SciMiner achieved 87.1% recall (percentage identification of targets in the given text), 71.3% precision (percentage accuracy of identified target) and 75.8% F-measure (harmonious average of recall and precision = (2 × recall × precision)/(recall + precision)) before manual revision [14]. In order to improve the accuracy of SciMiner's results, each target was manually reviewed and corrected by checking the sentences in which each target was identified. Approximately, 120 targets (~10% of the initially identified targets from the ROS-diabetes papers) were removed during the manual review process. The overall accuracy is expected to improve through the review process; however, the review process did not address targets missed by SciMiner, since we did not thoroughly review individual papers. Instead, 5 missed targets, whose associations with ROS-diabetes literature were available in the NCBI Gene2PubMed database, were added to the final ROS-diabetes target list (Additional File 2).

Conclusions

The present approach enabled us to collect a comprehensive list of ROS and diabetes related targets and led us to confirm the biological relevance to diabetic neuropathy of the selected ROS-diabetes targets. Using SciMiner to identify significantly enriched targets is applicable to other disease topics of interest by providing a more focused subset of literature for review and by highlighting targets common to multiple manuscripts.

List of abbreviations

NAD(P): nicotinamide adenine dinucleotide (phosphate); AGE: advanced glycoation end-products; PKC: protein kinase C; ROS: reactive oxygen species; MAPK: mitogen-activated protein kinase; JNK: c-Jun amino-terminal kinase; SOD1: superoxide dismutase 1; SOD2: superoxide dismutase 2; CAT: catalase; XDH: xanthine dehydrogenase; NCF1: neutrophil cytosolic factor 1; NOS3: nitric oxide synthase 3; CYBA: cytochrome b-245 alpha; NOS1: nitric oxide synthase 1; ALS: amyotrophic lateral sclerosis; BioCreAtIvE: Critical Assessment of Information Extraction systems in Biology; MiMI: Michigan Molecular Interactions; KEGG: Kyoto Encyclopedia of Genes and Genomes; STZ: streptozotocin; GO: gene ontology; DRG: dorsal root ganglia; CNS: central nervous system; HGNC: HUGO (human genome organization) nomenclature committee; PPI: protein-protein interaction; SEM: standard error mean.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JH participated in the study design, performed the literature mining and functional enrichment analyses, and drafted the manuscript. KAS participated in the study design and drafted the manuscript. ADS participated in the statistical analysis. YH carried out the quantitative RT-PCR assay. MP participated in the manuscript revision. DJS participated in the study design and manuscript revision. HVJ participated in the study design and manuscript revision. ELF participated in the study design and manuscript revision. All authors read and approved the final manuscript.

Acknowledgements

The authors acknowledge the expertise of Mr. John Hayes, Dr. Sang Su Oh and Ms. Carey Backus with regard to animal models of diabetes and neuropathy phenotyping methods.

Funding: This work was supported by National Institutes of Health (R01-LM008106 to DJS, U54-DA021519 (NCIBI) to DJS, HVJ, and ELF); Animal Models of Diabetic Complications Consortium (AMDCC U01-DK076160) to ELF; Program for Neurology Research and Discovery, and a Rackham Graduate School Predoctoral Fellowship (to JH).

References

  1. American-Diabetes-Association: Economic costs of diabetes in the U.S. In 2007.

    Diabetes Care 2008, 31(3):596-615. PubMed Abstract | Publisher Full Text OpenURL

  2. Edwards JL, Vincent AM, Cheng HT, Feldman EL: Diabetic neuropathy: mechanisms to management.

    Pharmacol Ther 2008, 120(1):1-34. PubMed Abstract | Publisher Full Text OpenURL

  3. Folli F, Guzzi V, Perego L, Coletta DK, Finzi G, Placidi C, La Rosa S, Capella C, Socci C, Lauro D, Tripathy D, Jenkinson C, Paroni R, Orsenigo E, Cighetti G, Gregorini L, Staudacher C, Secchi A, Bachi A, Brownlee M, Fiorina P: Proteomics reveals novel oxidative and glycolytic mechanisms in type 1 diabetic patients' skin which are normalized by kidney-pancreas transplantation.

    PLoS One 2010, 5(3):e9923. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Erol A: Insulin resistance is an evolutionarily conserved physiological mechanism at the cellular level for protection against increased oxidative stress.

    Bioessays 2007, 29(8):811-818. PubMed Abstract | Publisher Full Text OpenURL

  5. Pan JS, Hong MZ, Ren JL: Reactive oxygen species: a double-edged sword in oncogenesis.

    World J Gastroenterol 2009, 15(14):1702-1707. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Sarsour EH, Kumar MG, Chaudhuri L, Kalen AL, Goswami PC: Redox control of the cell cycle in health and disease.

    Antioxid Redox Signal 2009, 11(12):2985-3011. PubMed Abstract | Publisher Full Text OpenURL

  7. Vincent AM, Russell JW, Low P, Feldman EL: Oxidative stress in the pathogenesis of diabetic neuropathy.

    Endocr Rev 2004, 25(4):612-628. PubMed Abstract | Publisher Full Text OpenURL

  8. Schwedhelm E, Maas R, Troost R, Boger RH: Clinical pharmacokinetics of antioxidants and their impact on systemic oxidative stress.

    Clin Pharmacokinet 2003, 42(5):437-459. PubMed Abstract | Publisher Full Text OpenURL

  9. Shaw PJ: Molecular and cellular pathways of neurodegeneration in motor neurone disease.

    J Neurol Neurosurg Psychiatry 2005, 76(8):1046-1057. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Vincent AM, Kato K, McLean LL, Soules ME, Feldman EL: Sensory neurons and schwann cells respond to oxidative stress by increasing antioxidant defense mechanisms.

    Antioxid Redox Signal 2009, 11(3):425-438. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Vincent AM, Feldman EL: New insights into the mechanisms of diabetic neuropathy.

    Rev Endocr Metab Disord 2004, 5(3):227-236. PubMed Abstract | Publisher Full Text OpenURL

  12. Brownlee M: Biochemistry and molecular cell biology of diabetic complications.

    Nature 2001, 414(6865):813-820. PubMed Abstract | Publisher Full Text OpenURL

  13. Feldman EL: Oxidative stress and diabetic neuropathy: a new understanding of an old problem.

    J Clin Invest 2003, 111(4):431-433. PubMed Abstract | PubMed Central Full Text OpenURL

  14. Hur J, Schuyler AD, States DJ, Feldman EL: SciMiner: web-based literature mining tool for target identification and functional enrichment analysis.

    Bioinformatics 2009, 25(6):838-840. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Rebholz-Schuhmann D, Kirsch H, Arregui M, Gaudan S, Riethoven M, Stoehr P: EBIMed--text crunching to gather facts for proteins from Medline.

    Bioinformatics 2007, 23(2):e237-244. PubMed Abstract | Publisher Full Text OpenURL

  16. Plake C, Schiemann T, Pankalla M, Hakenberg J, Leser U: AliBaba: PubMed as a graph.

    Bioinformatics 2006, 22(19):2444-2445. PubMed Abstract | Publisher Full Text OpenURL

  17. Cheng D, Knox C, Young N, Stothard P, Damaraju S, Wishart DS: PolySearch: a web-based text mining system for extracting relationships between human diseases, genes, mutations, drugs and metabolites.

    Nucleic Acids Res 2008, (36 Web Server):W399-405. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Russell JW, Berent-Spillson A, Vincent AM, Freimann CL, Sullivan KA, Feldman EL: Oxidative injury and neuropathy in diabetes and impaired glucose tolerance.

    Neurobiol Dis 2008, 30(3):420-429. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Jimenez-Andrade JM, Herrera MB, Ghilardi JR, Vardanyan M, Melemedjian OK, Mantyh PW: Vascularization of the dorsal root ganglia and peripheral nerve of the mouse: implications for chemical-induced peripheral sensory neuropathies.

    Mol Pain 2008, 4:10. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Tarcea VG, Weymouth T, Ade A, Bookvich A, Gao J, Mahavisno V, Wright Z, Chapman A, Jayapandian M, Ozgur A, Tian Y, Cavalcoli J, Mirel B, Patel J, Radev D, Athey B, States D, Jagadish HV: Michigan molecular interactions r2: from interacting proteins to pathways.

    Nucleic Acids Res 2009, (37 Database):D642-646. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.

    Nat Genet 2000, 25(1):25-29. PubMed Abstract | Publisher Full Text OpenURL

  22. Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes.

    Nucleic Acids Res 2000, 28(1):27-30. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Vastrik I, D'Eustachio P, Schmidt E, Joshi-Tope G, Gopinath G, Croft D, de Bono B, Gillespie M, Jassal B, Lewis S, Matthews L, Wu G, Birney E, Stein L: Reactome: a knowledge base of biologic pathways and processes.

    Genome Biol 2007, 8(3):R39. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Fisher RA: On the interpretation of x2 from contingency tables, and the calculation of P.

    Journal of the Royal Statistical Society 1922, 85(1):87-94. Publisher Full Text OpenURL

  25. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing.

    Journal of the Royal Statistical Society Series B (Methodological) 1995, 57(1):289-300. OpenURL

  26. Sullivan KA, Hayes JM, Wiggin TD, Backus C, Su Oh S, Lentz SI, Brosius F, Feldman EL: Mouse models of diabetic neuropathy.

    Neurobiol Dis 2007, 28(3):276-285. PubMed Abstract | Publisher Full Text OpenURL

  27. Figueroa-Romero C, Sadidi M, Feldman EL: Mechanisms of disease: the oxidative stress theory of diabetic neuropathy.

    Rev Endocr Metab Disord 2008, 9(4):301-314. PubMed Abstract | Publisher Full Text OpenURL

  28. Hakim FA, Pflueger A: Role of oxidative stress in diabetic kidney disease.

    Med Sci Monit 2010, 16(2):RA37-48. PubMed Abstract OpenURL

  29. Ilieva H, Polymenidou M, Cleveland DW: Non-cell autonomous toxicity in neurodegenerative disorders: ALS and beyond.

    J Cell Biol 2009, 187(6):761-772. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Lee HB, Yu MR, Yang Y, Jiang Z, Ha H: Reactive oxygen species-regulated signaling pathways in diabetic nephropathy.

    J Am Soc Nephrol 2003, 14(8 Suppl 3):S241-245. PubMed Abstract | Publisher Full Text OpenURL

  31. Vincent AM, Hayes JM, McLean LL, Vivekanandan-Giri A, Pennathur S, Feldman EL: Dyslipidemia-induced neuropathy in mice: the role of oxLDL/LOX-1.

    Diabetes 2009, 58(10):2376-2385. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Cetinkalp S, Delen Y, Karadeniz M, Yuce G, Yilmaz C: The effect of 1alpha,25(OH)2D3 vitamin over oxidative stress and biochemical parameters in rats where Type 1 diabetes is formed by streptozotocin.

    J Diabetes Complications 2009, 23(6):401-408. PubMed Abstract | Publisher Full Text OpenURL

  33. Aliciguzel Y, Ozen I, Aslan M, Karayalcin U: Activities of xanthine oxidoreductase and antioxidant enzymes in different tissues of diabetic rats.

    J Lab Clin Med 2003, 142(3):172-177. PubMed Abstract | Publisher Full Text OpenURL

  34. Iaccio A, Collinet C, Gesualdi NM, Ammendola R: Protein kinase C-alpha and -delta are required for NADPH oxidase activation in WKYMVm-stimulated IMR90 human fibroblasts.

    Arch Biochem Biophys 2007, 459(2):288-294. PubMed Abstract | Publisher Full Text OpenURL

  35. Fontayne A, Dang PM, Gougerot-Pocidalo MA, El-Benna J: Phosphorylation of p47phox sites by PKC alpha, beta II, delta, and zeta: effect on binding to p22phox and on NADPH oxidase activation.

    Biochemistry 2002, 41(24):7743-7750. PubMed Abstract | Publisher Full Text OpenURL

  36. Bagi Z, Koller A, Kaley G: PPARgamma activation, by reducing oxidative stress, increases NO bioavailability in coronary arterioles of mice with Type 2 diabetes.

    Am J Physiol Heart Circ Physiol 2004, 286(2):H742-748. PubMed Abstract | Publisher Full Text OpenURL

  37. Fukuda M, Nakamura T, Kataoka K, Nako H, Tokutomi Y, Dong YF, Ogawa H, Kim-Mitsuyama S: Potentiation by candesartan of protective effects of pioglitazone against type 2 diabetic cardiovascular and renal complications in obese mice.

    J Hypertens 2010, 28(2):340-352. PubMed Abstract | Publisher Full Text OpenURL

  38. Fukuda M, Nakamura T, Kataoka K, Nako H, Tokutomi Y, Dong YF, Yasuda O, Ogawa H, Kim-Mitsuyama S: Ezetimibe ameliorates cardiovascular complications and hepatic steatosis in obese and type 2 diabetic db/db mice.

    J Pharmacol Exp Ther 2010. PubMed Abstract | Publisher Full Text OpenURL

  39. Ghosh S, Khazaei M, Moien-Afshari F, Ang LS, Granville DJ, Verchere CB, Dunn SR, McCue P, Mizisin A, Sharma K, Laher I: Moderate exercise attenuates caspase-3 activity, oxidative stress, and inhibits progression of diabetic renal disease in db/db mice.

    Am J Physiol Renal Physiol 2009, 296(4):F700-708. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Matsunami T, Sato Y, Sato T, Ariga S, Shimomura T, Yukawa M: Oxidative stress and gene expression of antioxidant enzymes in the streptozotocin-induced diabetic rats under hyperbaric oxygen exposure.

    Int J Clin Exp Pathol 2009, 3(2):177-188. PubMed Abstract | PubMed Central Full Text OpenURL

  41. Kameyama N, Arisawa S, Ueyama J, Kagota S, Shinozuka K, Hattori A, Tatsumi Y, Hayashi H, Takagi K, Wakusawa S: Increase in P-glycoprotein accompanied by activation of protein kinase Calpha and NF-kappaB p65 in the livers of rats with streptozotocin-induced diabetes.

    Biochim Biophys Acta 2008, 1782(5):355-360. PubMed Abstract | Publisher Full Text OpenURL

  42. Turner BJ, Talbot K: Transgenics, toxicity and therapeutics in rodent models of mutant SOD1-mediated familial ALS.

    Prog Neurobiol 2008, 85(1):94-134. PubMed Abstract | Publisher Full Text OpenURL

  43. Perez VI, Bokov A, Van Remmen H, Mele J, Ran Q, Ikeno Y, Richardson A: Is the oxidative stress theory of aging dead?

    Biochim Biophys Acta 2009, 1790(10):1005-1014. PubMed Abstract | PubMed Central Full Text OpenURL

  44. Ohlemiller KK, McFadden SL, Ding DL, Flood DG, Reaume AG, Hoffman EK, Scott RW, Wright JS, Putcha GV, Salvi RJ: Targeted deletion of the cytosolic Cu/Zn-superoxide dismutase gene (Sod1) increases susceptibility to noise-induced hearing loss.

    Audiol Neurootol 1999, 4(5):237-246. PubMed Abstract | Publisher Full Text OpenURL

  45. Kawase M, Murakami K, Fujimura M, Morita-Fujimura Y, Gasche Y, Kondo T, Scott RW, Chan PH: Exacerbation of delayed cell injury after transient global ischemia in mutant mice with CuZn superoxide dismutase deficiency.

    Stroke 1999, 30(9):1962-1968. PubMed Abstract | Publisher Full Text OpenURL

  46. Reaume AG, Elliott JL, Hoffman EK, Kowall NW, Ferrante RJ, Siwek DF, Wilcox HM, Flood DG, Beal MF, Brown RH, Scott RW Jr, Snider WD: Motor neurons in Cu/Zn superoxide dismutase-deficient mice develop normally but exhibit enhanced cell death after axonal injury.

    Nat Genet 1996, 13(1):43-47. PubMed Abstract | Publisher Full Text OpenURL

  47. DeRubertis FR, Craven PA, Melhem MF: Acceleration of diabetic renal injury in the superoxide dismutase knockout mouse: effects of tempol.

    Metabolism 2007, 56(9):1256-1264. PubMed Abstract | Publisher Full Text OpenURL

  48. Olofsson EM, Marklund SL, Behndig A: Enhanced diabetes-induced cataract in copper-zinc superoxide dismutase-null mice.

    Invest Ophthalmol Vis Sci 2009, 50(6):2913-2918. PubMed Abstract | Publisher Full Text OpenURL

  49. Hirschman L, Yeh A, Blaschke C, Valencia A: Overview of BioCreAtIvE: critical assessment of information extraction for biology.

    BMC Bioinformatics 2005, 6(Suppl 1):S1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1755-8794/3/49/prepub