Abstract
Background:
Diabetic peripheral neuropathy (DPN), the most common microvascular complication of type-2 diabetes mellitus (T2DM), results in nontraumatic lower-limb amputations. When DPN is not detected early, disease progression is irreversible. Thus, biomarkers for diagnosing DPN are needed.
Methods:
We analyzed three data sets of T2DM DPN: two for mouse models (GSE70852 and GSE34889) and one for a human model (GSE24290). We found common differentially expressed genes (DEGs) in the two mouse data sets and validated them in the human data set. To identify the phenotypic function of the DEGs, we overexpressed them in zebrafish embryos. Clinical information and serum samples of T2DM patients with and without DPN were obtained from the Korea Biobank Network. To assess the plausibility of DEGs as biomarkers of DPN, we performed an enzyme-linked immunosorbent assay.
Results:
Among the DEGs, only
Conclusion:
In the pathogenesis of DPN,
Introduction
By 2030, 1 in 10 adults aged over 20 years globally will be living with diabetes, 1 which is expected to increase health expenditure from USD 760 billion in 2019 to USD 825 billion in 2030. As type-2 diabetes mellitus (T2DM) progresses, metabolic disorders, such as hyperglycemia, insulin resistance, and obesity, cause microvascular disruptions that significantly increase morbidity and mortality in T2DM patients. 2 Diabetic peripheral neuropathy (DPN) is the most common microvascular complication of diabetes. DPN is characterized by a loss of sensation beginning distally in the lower extremities, accompanied by pain and significant morbidity. The incidence and prevalence of DPN vary widely because of differences in definitions used among studies. 3 In epidemiological studies, the prevalence of DPN has been reported as 6100 per 100,000 people, and up to 51% of T2DM patients suffer from DPN. 4
Individuals with DPN complain of unbearable stabbing, burning, and tingling sensations in the limbs, in addition to depression, anxiety, and sleep deprivation. 5 Moreover, sensory deficits and poor wound healing often result in foot ulceration and nontraumatic lower-limb amputations. 6 Amputations in patients with diabetes affect their quality of life and are associated with substantially increased mortality, especially major amputations (5-year mortality rate of 44–68%, which is greater than those of prostate and breast cancer). 7
Early manifestations of DPN often go undiagnosed until the disease is well-developed, at which point they appear to be irreversible. Screening for early symptoms and signs of DPN is critical in clinical practice to improve patient prognosis and quality of life and reduce healthcare costs. The gold standard methods for detecting DPN are nerve conduction studies. 8 However, these studies are often impractical to implement as they are labor-intensive, time-consuming, and costly. Because of the lack of treatments for nerve damage, preventing DPN is the best approach for improving patient prognosis; thus, it is essential to identify a biomarker for early detection of DPN.
Biomarkers provide precise diagnostic and predictive data and objectively identify pathogenic processes. A recent study reported that inflammation, oxidative stress, lipid and carbohydrate metabolism, axonogenesis, mitochondrion, and peroxisome proliferator–activated receptor signaling events are dysregulated in DPN and are important pathomechanisms in DPN development. 9 Molecules associated with these pathways have been suggested as predictive biomarkers. However, traditional biomarker discovery processes have limitations, and few have been adopted in the clinic.
With the accumulation of gene expression data and the development of computational analyses, researchers, particularly in oncology, have begun investigating biomarkers for developing immunotherapeutic and other targeting agents to maximize the efficacy of treatments, while minimizing toxicity. 10 The aim of this study was to identify a biomarker for DPN by gathering gene expression data and performing computational analysis, experimental validation in animals, and final validation using immunoaffinity assays.
Materials and methods
Microarray data
We searched transcriptomics data from the Gene Expression Omnibus data set repository (GEO data set; https://www.ncbi.nlm.nih.gov/gds) using the search query “diabetes mellitus neuropathy.” We excluded data for animals with type-1 diabetes mellitus and samples from the central nervous system. We downloaded data sets using the getGEO function in GEO query
In silico analyses
A linear model for microarray analysis (LIMMA)
Maintenance of adult zebrafish and embryos
Wild-type zebrafish AB and
In vitro transcription for mRNA injection
For overexpression experiments, full length was cloned into the pcs2p + vector. mRNA was synthesized using an mMESSAGE mMACHINE SP6 (Invitrogen, Carlsbad, CA, USA). The synthesized mRNA was injected into single-cell-stage zebrafish embryos. Because secretory leukocyte protease inhibitor (SLPI) is not evolutionally conserved in vertebrates, we did not synthesize
Live imaging of zebrafish embryos
Zebrafish embryos were anesthetized with tricaine for observation at 5 days after fertilization. Anesthetized embryos were observed using a Leica M165C microscope (Wetzlar, Germany).
Quantitative reverse-transcription polymerase chain reaction
Total RNA from zebrafish embryos was extracted using TRIzol reagent (Molecular Research Center, Cincinnati, OH, USA). Embryos were homogenized using a pestle and dissolved using TRIzol reagent. Total RNA was isolated, and cDNA was synthesized as previously described. 16 A quantitative reverse-transcription polymerase chain reaction (qRT-PCR) was performed using PowerUP SYBR Green Master Mix (Thermo Fisher Scientific, Waltham, MA, USA). The expression level of the target gene was normalized to that of β-actin as an endogenous control.
Confocal imaging of transgenic zebrafish
To observe peripheral nervous system development in zebrafish embryos,
Statistical analysis of data from zebrafish experiments
Statistical analyses were performed using Student’s
Patient serum and data acquisition
Serum from T2DM patients diagnosed with and without DPN was obtained from the Korea Biobank Network. Serum from 20 T2DM patients diagnosed with DPN and 9 T2DM patients without DPN was obtained from Kangwon National University Hospital. Serum from five T2DM patients without DPN was obtained from the Gyeongsang National University Hospital. A nerve conduction study confirmed the diagnosis of DPN. Information on baseline characteristics including age, sex, HbA1c, and comorbidities, such as hypertension, coronary artery disease, cerebrovascular infarction, and cancer was also obtained. We examined the levels of NPY and SLPI in the patient serum to identify a serum biomarker for DPN. The Institutional Review Board of the Gachon University Gil Medical Center (approval id. GCIRB2020-044) approved this study.
Enzyme-linked immunosorbent assay
Serum NPY was measured by enzyme-linked immunosorbent assay (ELISA) (Cat. No. NBP2-76680, Novus Biologicals, Littleton, CO, USA). Standards or samples were added to each well, followed immediately by the addition of a biotinylated detection antibody and incubation for 45 min at 37°C.
After washing, the horseradish peroxidase–conjugated NPY antibody was added and incubated for 30 min at 37°C. After washing, the substrate reagent was added to each well for 15 min at 37°C. Stop solution was added to each well, and the absorbance was read at 450 nm using a microplate reader (Spectra MAX M2e, Molecular Devices, San Jose, CA, USA).
We prepared an SLPI-targeted capture antibody for quantitative sandwich enzyme immunoassay (Cat. No. DY1274-05, R&D Systems, Minneapolis, MN, USA) according to the manufacturers’ instructions. Standards or serum samples were added to each well and incubated for 2 h at room temperature. The plates were washed to remove any residual contaminants, and the detection antibody was added to each well. After aspiration/washing, the working solution of streptavidin-horseradish peroxidase was added and incubated for 20 min at room temperature. The plates were washed, and substrate solution was added per well for incubation for 20 min at room temperature.
Finally, stop solution was added to each well. The optical density was immediately determined using a microplate reader (Molecular Devices Spectra MAX M2e) at an absorbance of 450 nm.
Statistical analysis of patient characteristics and ELISA results
Continuous variables were analyzed using Mann–Whitney
Results
Identification of DEGs related to DPN
Using a threshold of|log2FC| > 1.5 and BH-adjusted

Diagram of analysis procedure: data collection, analysis, common DEGs selection, and experimental validation.

DEGs in each GEO data set and expression of selected genes: (a) DEGs between two groups were identified via the LIMMA package and visualized using a volcano plot. Significantly expressed genes are represented as red dots. (b) Expression of
Overexpression of npy induced the development of the peripheral nerve system in zebrafish embryos
Zebrafish are commonly used as an
The zebrafish ortholog of human NPY,

Overexpression of
The effect of
NPY level was higher in the serum of T2DM patients with DPN than in that of those without DPN
There were 20 T2DM patients with DPN and 14 T2DM patients without DPN. The mean ages of T2DM patients with and without DPN were 69.5 and 61.8 years, respectively (
In ELISA, the mean levels of NPY in the serum of T2DM patients with and without DPN were 2953 and 2882 pg/mL, respectively (

Results of ELISA in type-2 diabetes patients: (a) serum level of NPY in type-2 diabetes patients with diabetic PN and without diabetic PN. (*
Discussion
We analyzed three microarray data sets acquired from GEO (two for mouse sciatic nerve with T2DM DPN and one for human sural nerve with DPN). Among the common DEGs from the two mouse data sets, only
Despite the importance of DPN, current pharmacotherapies are limited to resolving symptoms rather than the disease mechanisms, as the precise pathogenic mechanisms of DPN are unknown. 18 The known pathomechanisms of DPN include oxidative stress, inflammation, microvascular alterations, nerve degeneration, and regrowth. Oxidative stress plays a key role in the pathogenesis of DPN. 19 Thus, biomarkers for systemic oxidative stress have been investigated.
A large population-based Cooperative Health Research study was conducted in Augsburg, F4, and included a large proportion of individuals with prediabetes and T2DM. This study showed that myeloperoxidase was positively associated with DPN and that baseline superoxide dismutase levels were related to the incidence of DPN.
20
Methylglyoxal is a highly reactive metabolite formed during pyruvate metabolism. High baseline blood methylglyoxal was significantly associated with DPN risk in a Danish cohort.
21
Cross-sectional studies revealed a positive correlation between proinflammatory cytokines (interleukin-6 and tumor necrosis factor
Biomarkers are a cornerstone of medical diagnostics and are used in various applications, such as disease diagnosis, prognosis, and patient stratification. An array of molecular pathways involved in DPN pathogenesis has been reported; however, there are no disease-targeted treatments or biomarkers for predicting the onset or progression of DPN. In this context, hypothesis-free approaches involving multiomics technologies have been employed to identify biomarkers associated with pathways involved in DPN. Hur and colleagues 13 examined gene expression profiles in sural nerve samples from two groups of patients with progressing or nonprogressing DPN. A literature-derived co-citation network of the DEGs revealed five subnetworks centered on apolipoprotein E, c-Jun, leptin, serpin peptidase inhibitor E type 1, and peroxisome proliferator-activated receptor-gamma. Ridge regression models, including 14 DEGs, correctly identified the progression status of 92% of the patients with DPN.
In the clinic, peripheral nerve biopsies, such as for the sural nerve, are rarely performed; thus, human tissue sources in omics studies are limited. McGregor and colleagues
25
used a systems biology approach and demonstrated a highly conserved pathway across human and murine DPN models. Hinder and colleagues
26
examined gene expression changes in the sciatic nerve and dorsal root ganglia of
We identified
Several studies have investigated the function of NPY and its receptors in regulating food intake, blood pressure, seizure activity, anxiety, depression, and the pathogenesis of neurodegenerative disorders.28,29 A growing body of research has suggested that NPY and its receptors influence brain activity and drive feeding behaviors; research in this area is relatively young. Howell and colleagues 30 found that NPY promoted dentate gyrus neural precursor cell proliferation. They also showed that NPY1 receptor knockout mice exhibited reduced proliferation of precursor cells and the generation of immature neuroblasts in the dentate gyrus compared to wild-type mice.
Furthermore, NPY treatment promoted the differentiation of newly generated stem cells toward hippocampal progenitors. 31 Studies have shown that peripheral nerve injury causes an increase in NPY expression in adult dorsal root ganglion neurons in mice and rats.32,33 It appears to mediate antinociceptive actions via the NPY1 and NPY2 receptors in sciatic nerve injury models of neuropathic pain and plantar incision models of postoperative pain.34,35
However, the roles of NPY in central nervous system neurogenesis and antinociception and its functions in the peripheral nervous system have not been widely examined. We found one study showing that in the peripheral nervous system, NPY promotes angiogenesis through the NPY2 and NPY5 receptors during ischemic revascularization and wound healing. 36
Matyal and colleagues proved that plasma NPY levels were increased in patients with diabetes and nondiabetic subjects undergoing cardiac surgery under cardiopulmonary bypass.
37
Atrial tissue
We compared the serum levels of NPY between patients with T2DM DPN and without T2DM DPN and showed that the serum level of NPY was increased in patients with T2DM DPN. As the AUC of NPY was 0.7179, it may not suffice as a sole biomarker for detecting DPN in T2DM patients. Moreover, early-stage DPN is usually asymptomatic, leading to a delay in diagnosis and poor prognosis. If the serum level of NPY is serially investigated for T2DM patients and is revealed to be useful for detecting DPN in the early stage, it could be used in combination with current diagnostic tools and could impact T2DM patients’ prognosis. For this purpose, further prospective study should be performed on a larger cohort of T2DM patients.
The consistency of mRNA expression and serum levels of NPY in T2DM patients with DPN suggests that NPY is a key biomarker for detecting T2DM DPN. Based on the phenotypic changes observed in zebrafish overexpressing NPY, NPY may affect nerve regeneration in the peripheral nervous system.
This study had some limitations. Our clinical samples for the ELISA were relatively small. When we requested for serum samples from T2DM patients with DPN and without DPN from the Korea Biobank, only samples from 20 T2DM patients with DPN and 14 T2DM patients without DPN were available. Even though our sample size was not adequate (low statistical power) to reach a definite conclusion, baseline characteristics and comorbidities were not different between the two groups. Thus, our results could provide insights into the possibility of using NPY as a biomarker for DPN. In addition, we did not determine the molecular mechanism of NPY in T2DM DPN. However, our study suggests that NPY is mechanistically involved in nerve regeneration and can be used as a biomarker of neuropathy in patients with T2DM DPN.
Conclusion
We overexpressed
Supplemental Material
sj-docx-1-taj-10.1177_20406223211041936 – Supplemental material for Neuropeptide Y: a potential theranostic biomarker for diabetic peripheral neuropathy in patients with type-2 diabetes
Supplemental material, sj-docx-1-taj-10.1177_20406223211041936 for Neuropeptide Y: a potential theranostic biomarker for diabetic peripheral neuropathy in patients with type-2 diabetes by Noo Ree Cho, Yeuni Yu, Chang-Kyu Oh, Dai Sik Ko, Kyungjae Myung, Yoonsung Lee, Hee Sam Na and Yun Hak Kim in Therapeutic Advances in Chronic Disease
Footnotes
Acknowledgements
The biospecimens and data used in this study were provided by the Biobank of Gyeongsang National University Hospital and Kangwon University Hospital, a member of the Korea Biobank Network. N.R.C., Y.Y., and C.-K.O. contributed equally to this study.
Author contributions
N.R.C. and Y.H.K. made substantial contributions to the conception and design. D.S.K., K.M., and Y.L. collected data, and N.R.C. and Y.Y. performed statistical analysis. N.R.C. and H.S.N. analyzed and interpreted the data. C.-K.O. produced data using zebrafish. N.R.C. and C.-K.O. wrote the manuscript. D.S.K., H.S.N., and Y.H.K. made critical revision.
Conflict of interest statement
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by grants from the Basic Science Research Program (grant no. NRF-2020R1C1C1003741), Medical Research Center (grant no. NRF-2018R1A5A2023879) through the National Research Foundation of Korea Grant funded by the Korean government, the Institute for Basic Science (grant no. IBS-R022-D1), and Korean Medical Institute (KMI).
Ethical statement
The Institutional Review Board of the Gachon University Gil Medical Center (approval id. GCIRB2020-044) approved this study. The acquisition of informed consent was waived by the Institutional Review Board of the Gachon University Gil Medical Center as only unidentifying information was collected. All experiments using zebrafish embryos were performed according to the guidelines of the Ulsan National Institute of Science and Technology Institutional Animal Care and Use Committee (approval number: UNISTIACUC-15-14, date: 2016-10-11).
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
