ORIGINAL RESEARCH article

Front. Genet., 26 September 2022

Sec. Genetics of Common and Rare Diseases

Volume 13 - 2022 | https://doi.org/10.3389/fgene.2022.965816

Identification and validation of P4HB as a novel autophagy-related biomarker in diabetic nephropathy

  • 1. Department of Nephrology, Qilu Hospital of Shandong University, Jinan, Shandong, China

  • 2. Laboratory of Basic Medical Sciences, Qilu Hospital of Shandong University, Jinan, Shandong, China

  • 3. Department of General Practice, Qilu Hospital of Shandong University, Jinan, Shandong, China

Abstract

Diabetic nephropathy (DN), a frequent microvascular complication of diabetes, has been recognized as a primary cause of chronic kidney disease (CKD) and end-stage renal disease (ESRD). Previous studies found that autophagy of renal tubular epithelial cells plays an important role in DN pathogenesis. Our research aimed to investigate the differentially expressed autophagy-related genes (DEARGs) between DN and healthy renal tubule samples and identify a novel autophagy-related biomarker associated with tubulointerstitial injury in DN. In this study, gene expression profiles of renal tubules from 10 DN patients and 24 healthy controls in the GSE30122 dataset were analyzed, and 43 DEARGs were identified by bioinformatics analysis. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and correlation analysis were performed on DEARGs, and the hub gene prolyl 4-hydroxylase subunit beta (P4HB) was screened by protein–protein interaction and verified by utilizing other datasets and stimulating HK-2 cells under high glucose concentration. We found that the expression of P4HB in renal tubules was correlated with renal function. In summary, our research provided novel insights for comprehension of DN molecular mechanisms and identified P4HB as a novel autophagy-related biomarker of DN.

Introduction

Diabetic nephropathy (DN), a common microvascular complication of diabetes, has been recognized as a primary cause of chronic kidney disease (CKD) and end-stage renal disease (ESRD) in many developed and developing countries (; ). According to recent reports, DN has accounted for 20 %–40 % of patients requiring kidney replacement therapy worldwide, contributing to more than 950,000 deaths worldwide each year (). DN has become a serious global healthcare problem.

Although glomerular injury is a major pathological manifestation of DN, there is growing evidence suggesting that tubulointerstitial pathological changes, such as tubular atrophy, interstitial fibrosis, autophagy, and apoptosis of tubular epithelial cells, play a crucial role in DN development (; C. ; ). Apart from all other reasons, non-enzymatic glycation is considered one of the major reasons behind diabetes, including DN. It was found that DNA-AGEs and auto-antibodies against glycated DNA are related to diabetic microvascular complications, such as DN (). At the same time, it was also reported that some commonly used drugs like ezetimibe and rosuvastatin can show strong antidiabetic and renal protective effects by targeting AGE/RAGE-associated signaling (; ). In addition, some common components in plants such as ellagic acid and lycopene play a protective role in diabetes and DN by reducing the formation of AGEs (; ).

Autophagy is a cellular process in which damaged organelles, protein aggregates, and other macromolecules are degraded in the cytoplasm (). Autophagy dysfunction is associated with pathogenesis of various diseases (; C. ). Several pathways affect the biological function of DN by influencing autophagy. reported that activation of the epidermal growth factor receptor (EGFR) signaling pathway can exacerbate kidney damage by inhibiting autophagy. However, autophagy-related genes (ARGs) in DN are still largely unknown and require further exploration. Exploration and illumination of differentially expressed autophagy-related genes (DEARGs) in DN will provide us with novel biomarkers for treatment of DN. Bioinformatics is a method for efficiently and accurately processing large quantities of data, providing valuable information to patients. Nevertheless, studies on the investigation of the expression of diabetic renal tubular interstitial genes and autophagy through bioinformatics are still lacking.

Woroniecka et al. completed the DN-related dataset GSE30122 for analysis of gene expression differences between DN patients and healthy controls (). In this study, we re-analyzed their dataset from other perspectives to explore DEARGs in DN vs. normal human renal tubular interstitial genes. Subsequently, we performed Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and correlation analysis on DEARGs. Protein–protein interaction (PPI) showed that prolyl 4-hydroxylase subunit beta (P4HB) was identified as an autophagy-related hub gene for DN. Finally, we further validated the upregulated expression of P4HB in DN tubulointerstitium by exploiting another database (GSE104954) () and the European Renal cDNA Bank (ERCB) cohort, as well as establishing an in vitro model, and the expression of P4HB in renal tubules was correlated with renal function. The experiments’ schematic workflow is displayed in Figure 1. Our study suggested that P4HB is a potential key biomarker in the pathogenesis of renal tubular injury in DN.

FIGURE 1

Materials and methods

Data download and preprocessing

The Human Autophagy Database (HADb; http://www.autophagy.lu/index.html) was used to obtain genes involved in autophagy. Genomic and transcriptomic datasets of DN and healthy renal tubule samples were obtained using Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/). Data from GSE30122 (GPL571 platform, Affymetrix Human Genome U133A 2.0 Array) included the data of 10 patients with DN (GSM757014–GSM757023) and 24 control human kidney tubules (GSM757024–GSM757035, GSM758498–GSM758509). The basic characteristics of 10 DN patients are summarized in Supplementary Table S1, and the histological evaluation of the 10 DN patients showed that renal tubular atrophy, interstitial fibrosis, vascular sclerosis, and mesangial matrix dilatation increased significantly (). Data from GSE104954 (GPL22, 945 platform, Affymetrix Human Genome U133 Plus 2.0 Array, and GPL 24,120 platform, Affymetrix Human Genome U133A Array) were used for validation. The “sva” () and “limma” () packages in R software were used to normalize raw data in batches.

Differentially expressed autophagy-related gene analysis

The reproducibility of the GSE30122 data was examined using the principal component analysis (PCA) method. The “limma” package was utilized to investigate the differential expression of genes related to autophagy. Genes were identified as differentially expressed genes based on an adjusted p-value of <0.05 and an absolute fold-change value of >1.5. Heatmaps were created using the “pheatmap” package in R, and the volcano plots and box plots were performed using the Sangerbox tools, a free online platform for data analysis (http://vip.sangerbox.com/home.html).

Gene functions and correlation analysis

Enrichment analyses were conducted for DEARGs. GO and KEGG pathway enrichment analyses were performed using DAVID version 6.8 (https://david.ncifcrf.gov/conversion.jsp), a commonly used tool for detailed analysis and classification of genes and protein functions in bioinformatics research. The GO analysis included cellular composition (CC), biological process (BP), and molecular function (MF). Enrichment results and correlation analysis were drawn by https://www.bioinformatics.com.cn, a free online platform for data analysis and visualization.

Protein–protein interaction network construction

The PPI network was constructed on the basis of the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) online analysis (https://string-db.org/). PPI network visualization and analyses were performed with Cytoscape (Version 3.9.1). Thirty hub genes related to autophagy were identified by the Density of Maximum Neighborhood Component (DMNC), and P4HB was screened as a hub gene. An online platform (https://www.bioinformatics.com.cn) was also used to show the receiver operating curve (ROC) to evaluate the ability of P4HB to discriminate between DN patients from healthy controls.

P4HB validation and correlation analysis with renal function

The differential gene expression of P4HB between DN and healthy renal tubule samples was verified in the GSE104954 dataset and ERCB cohort (31 healthy controls and 17 DN patients) (). Then, the association analysis of P4HB expression and clinical characteristics was validated in ERCB using Pearson’s correlation analysis by using the Nephroseq v5 online database (http://v5.nephroseq.org).

Cell culture and treatments

The human proximal tubular cell line HK-2 was purchased from the National Collection of Authenticated Cell Cultures and cultured in Dulbecco’s modified Eagle’s medium (c11885500BT, Gibco) by adding 10 % fetal bovine serum (10270-106, Gibco) and 1 % penicillin–streptomycin (p1400, Solarbio) at 37°C and 5 % CO2 in a humidified environment. High glucose concentration (30 mM, G8644, Sigma) for 24 h () was used to cause HK-2 cell damage, and the addition of 5.5 mmol/L glucose served as the control.

RNA extraction and RT-qPCR analysis

Total RNA was extracted from cultured HK-2 cells using the RNAfast200 Kit (Fastagen), and cDNA was synthesized using the SureScriptTM First-Stand cDNA Synthesis Kit (QP056T, GeneCopoeia). Real-time quantitative polymerase chain reaction (RT-qPCR) was conducted using SYBR Green reagent (1725201, Bio-Rad) on a Bio-Rad CFX PCR System (Bio-Rad). The procedure was repeated three times for each sample. The primers are given in Table 1. Gene expression analysis was performed using the 2−ΔΔCt method, and expression levels were normalized to those of GAPDH.

TABLE 1

GenePrimer sequence (from 5′ to 3′)
P4HBForward: 5′-CTG​CGG​AAA​AGC​AAC​TTC​GC-3′
Reverse: 5′-CCA​CAC​CAA​GGG​GCA​TAG​AA-3′
NGALForward: 5′-AGC​ACC​AAC​TAC​AAC​CAG​CAT-3′
Reverse: 5′-TTG​GGA​CAG​GGA​AGA​CGA​TG-3′
GAPDHForward: 5′-GCA​CCG​TCA​AGG​CTG​AGA​AC-3′
Reverse: 5′-TGG​TGA​AGA​CGC​CAG​TGG​A-3′

Quantitative PCR primers used in the study.

Protein extraction and Western blotting analysis

Total proteins were extracted by the incubation of cultured HK-2 cells with radioimmunoprecipitation assay (RIPA) buffer (P0013D, Beyotime) and adding 1 % phenylmethylsulfonyl fluoride (PMSF) (329-98-6, Solarbio) and quantified by Bio-RAD assays. The same amount of protein was separated by 10 % sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis (SDS-PAGE), and proteins were transferred to polyvinylidene difluoride (PVDF) membranes. The membranes were then blocked with 5 % skim milk and subsequently incubated with primary antibodies against P4HB (1:1,000 dilution, Cat.137,110, Abcam) and GAPDH (1:6,000 dilution, Cat. 60004-1-Ig, ProteinTech) overnight at 4°C. After washing, the proteins were incubated with horseradish peroxidase (HRP)-conjugated secondary antibody (1:6,000 dilution, Cat. SA00001-2, ProteinTech) for 1 h at room temperature. The bands were visualized using an enhanced chemiluminescence (ECL) system, and the signal intensity was quantitatively processed by ImageJ software.

Immunohistochemistry and immunofluorescent staining

The paraffin-embedded kidney sections of the kidney tissues of three DN patients diagnosed by renal biopsy from the pathology department of Qilu Hospital affiliated with Shandong University and the paraffin-embedded kidney sections of healthy adjacent kidney tissues of three individuals who underwent tumor nephrectomy (no diabetes or other kidney diseases) were immunohistochemically stained by P4HB, as approved by the Research Ethics Committee Qilu Hospital of Shandong University (NO. KYLL-2020(KS)-030). These sections were incubated with an anti-P4HB antibody (1:1,000 dilution, Cat.137110, Abcam) at 4°C overnight. The general two-step method was used for detection and incubation, and the 3, 3′-diaminobenzidine (DAB) chromogenic kit was used for immunohistochemical staining. Images were collected and analyzed using the NIS Element software and Nikon microscope imaging system, and ImageJ software was used for quantitative analysis.

For immunofluorescent staining of cells, the cells were fixed with 4 % paraformaldehyde, incubated with TritonX-100 (1 %) for 20 min, and then incubated with the P4HB primary antibody (1:200 dilution, Cat.137110, Abcam) overnight at 4°C, Then, the cells were incubated with goat anti-rabbit immunoglobulin (Ig)G DyLight 594 (1:500, Cat. A23420, Abbkine Scientific Company) coupled with a fluorescent probe in the dark at room temperature for 1 h, and the nuclei were observed by DAPI (AR1176, Boster Bio) staining. Images were captured by a fluorescence microscope (Olympus).

Statistical analysis

GraphPad Prism 7 and R software version 4.1.3 were used for statistical analysis. A two-tailed unpaired t test was used for comparisons between the two groups. Differences were considered statistically significant at p <0.05. (*p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001).

Results

Identification of differentially expressed autophagy-related genes between the DN and healthy tubule samples

The renal tubular transcriptome data GSE10322 were used for further analysis to investigate the role of ARGs in DN pathogenesis. The distribution of DN and normal samples was shown by the PCA results (Figure 2A). The expression of 222 ARGs in 10 DN patients and 24 normal samples was subsequently analyzed, and 43 ARGs showed differential expression in DN with an adjusted p value <0.05 and an absolute fold-change value of >1.5, including 38 upregulated and 5 downregulated genes (Table 2; Figures 2B–D).

FIGURE 2

TABLE 2

Gene symbollogFCChangep-valueAdj.p-valueChromosome
CASP12.815552955Up1.56195E-163.09266E-1411q23
CCR22.58694031Up1.36828E-101.35459E-083p21.31
PEA151.786883586Up1.39952E-050.000230921q21.1
CASP41.615423503Up1.0762E-074.26175E-0611q22.2-q22.3
RAB11A1.5878937Up0.0002490190.00197222915q22.31
CCL21.499910231Up1.94226E-076.40947E-0617q11.2-q12
FKBP1B1.303245Up1.08E-036.40E-032p23.3
CDKN1B1.297208099Up0.0001919670.00158372712p13.1-p12
ATG121.258743997Up0.0004886660.0032251955q21-q22
FAS1.258208019Up3.27484E-067.20464E-0510q24.1
HSPB81.211358555Up0.0021462350.00965673212q24.23
CAPN21.191207762Up0.0037022170.0145116481q41-q42
HSPA51.16304934Up0.0002737480.0020846949q33.3
RB11.155716399Up5.24217E-050.00061055913q14.2
CASP31.138499883Up2.22484E-081.46839E-064q34
HSP90AB11.122569062Up0.0050940940.018338746p12
IKBKB1.117062221Up8.29611E-084.10657E-068p11.2
TNFSF101.110782311Up6.15161E-071.74003E-053q26
MYC1.04536829Up2.70E-053.89E-048q24.21
APOL10.982670415Up3.40551E-050.00044952822q13.1
ITGA30.926426396Up0.0019076680.00921263917q21.33
SH3GLB10.871043634Up0.0035915730.0145116481p22
ARNT0.835171835Up3.97754E-050.000492221q21
LAMP20.806156698Up6.16193E-050.000658743Xq24
ATIC0.800030145Up2.75244E-050.0003892742q35
ITGB10.788846065Up7.99052E-050.00075339210p11.2
WIPI10.744366027Up0.0041948940.01597286517q24.2
GNAI30.74200647Up0.0011439940.0064717391p13
FKBP1A0.7292738Up0.0092016320.02760489520p13
CALCOCO20.719431492Up0.0003479260.00255145717q21.32
P4HB0.71230372Up0.001338950.00697663617q25
CXCR40.692229997Up7.17148E-050.0007099772q21
TBK10.641395881Up0.0019941290.00940089312q14.1
CFLAR0.618087278Up6.32127E-050.0006587432q33-q34
TM9SF10.614734304Up0.0017425860.00862580114q11.2
CHMP4B0.60073086Up0.0059996510.02073649314q12
GABARAPL20.598015027Up0.0075543470.02492934516q22.1
SPNS10.592296801Up0.0010988470.00639916916p11.2
STK11−0.607962624Down8.98484E-060.000177919p13.3
HDAC6−0.628923381Down1.13142E-050.000203656Xp11.23
SIRT2−0.737375167Down1.05727E-062.61673E-0519q13
GABARAPL1−0.760812687Down8.81349E-050.00079321415q26.1
GRID1−0.919198606Down0.0087458390.02748692210q23.3

43 differentially expressed autophagy-related genes in the tubulointerstitium between DN and normal renal tissue.

Functional enrichment analysis of the differentially expressed autophagy-related genes

GO and KEGG pathway enrichment analyses were performed to further explore the biological functions of these DEARGs. The most involved processes or components in the GO included autophagy and apoptotic process (biological process), autophagosome membrane and membrane raft (cell component), and cysteine-type endopeptidase activity involved in the apoptotic signaling pathway and ubiquitin protein ligase binding (molecular function) (Figures 3A, B; Supplementary Table S2). The DEARGs mainly involved in the autophagy-associated process were shown by the KEGG pathway enrichment analysis (Figures 3C, D; Supplementary Table S3).

FIGURE 3

Hub gene identification and validation

A correlation analysis was conducted to investigate the expression relevance of these DEARGs. The relationship of 43 DEARGs in the GSE30122 dataset is shown in Figure 4A. The human protein interaction database (String) was chosen to identify the interactions among DEARGs, and 30 hub genes were calculated by the DMNC algorithm. P4HB was screened as a hub DEARG in DN (Figure 4B). We utilized another independent dataset GSE104954 and an ERCB cohort to further verify the change in P4HB expression, showing that P4HB expression was significantly upregulated in DN tubule samples (Figures 4C–F). Additionally, the ROC analysis demonstrated that the expression of P4HB showed an excellent diagnostic value for DN patients and healthy controls (Figures 4C–F).

FIGURE 4

Validation of P4HB expression under high glucose stimulation

After determining from other datasets that the expression of P4HB was significantly elevated in DN patients, the transcriptional and protein levels of P4HB were explored after high glucose stimulation and normal glucose treatment of HK-2 cells. RT-qPCR (Figure 5A) and Western blot analysis (Figures 5B, C) further validated that P4HB expression was elevated in HK-2 cells under high glucose (30 mM) stimulation, accompanied by increased expression of NGAL, a marker for renal tubular injury, indicating that there was significant damage in HK-2 cells. Simultaneously, cell immunofluorescence showed that the expression of P4HB in HK-2 cells increased significantly after high glucose stimulation (Figure 5D).

FIGURE 5

Clinical relevance of P4HB expression and immunohistochemical validation of P4HB expression

The correlation analysis of the ERCB cohort was performed to validate the relationship between the expression of P4HB in renal tubular and renal function. P4HB expression was positively correlated with serum creatinine levels (r = 0.351, p = 0.028) and negatively correlated with glomerular filtration rate (GFR) (r = −0.472, p = 0.002) (Figures 6A,B). Additionally, we performed immunohistochemical staining of P4HB in kidney tissues of healthy controls and DN patients, and found that the expression of P4HB in the renal tubules of DN patients was significantly higher than that in normal kidneys (Figures 6C,D).

FIGURE 6

Discussion

DN is one of the most serious microvascular complications of diabetes, accounting for approximately 30 %–40 % of end-stage kidney disease patients worldwide (; L. ). Poor prognosis and low quality of life are distinctive characteristics of patients with DN (). Therefore, novel and satisfactory strategies are urgently needed to treat DN. There is growing evidence that multiple biological functions are engaged in the pathogenesis of DN, such as immunity, inflammation, and autophagy (; ; ; ). In our research, starting from autophagy, for the first time, we identified DEARGs in DN vs. normal human renal tubular interstitial cells by bioinformatics analysis, as well as determined a new autophagy-related biomarker for DN, thus providing novel insights into the tubulointerstitial pathogenesis of DN and contributing to the identification of novel therapeutic targets.

P4HB, a member of the protein disulfide isomerase (PDI) family, is a multifunctional protein capable of catalyzing the generation and reorganization of disulfide bonds (). P4HB, as an autophagy-related gene, can be detected in various diseases that involve inflammation and apoptosis, including cancer, endocrine diseases, and skin diseases. Elevated expression of P4HB has been reported in several solid tumors, such as ovarian cancer (), bladder cancer (; ), and prostate cancer (; ). Additionally, reported that targeting P4HB can reduce inflammation and melanogenesis of the skin. Furthermore, previous studies found that PDIA1 contributes to oxidative maturation of proinsulin in the endoplasmic reticulum to support insulin production and ß-cell health in diet-induced obesity (), indicating that P4HB could indirectly influence insulin production and ß-cell health.

One study has determined that P4HB was substantially elevated as an ARG in kidney renal clear cell carcinoma (KIRC) and showed high diagnostic and prognostic ability (). There are also reports demonstrating that P4HB overexpression was associated with poor prognosis in human KIRC (; ), indicating that overexpression of P4HB is an adverse prognostic factor. Additionally, reported that the self-antigen P4HB located on the cell membrane of kidney cells could be cross-recognized by anti-HU1 (a conserved peptide derived from DNABII proteins) and induce lupus nephritis (LN). Hence, evidence suggests that P4HB might be a key biomarker and therapeutic target for human kidney diseases.

In this study, transcriptomic variations of ARGs in the gene expression profiles of 10 DN and 24 healthy renal tubule samples were analyzed, and 43 DEARGs were identified in DN samples compared with healthy samples. The GO and KEGG enrichment analyses of DEARGs were subsequently performed. These genes were mostly enriched in autophagy-related biological processes, such as autophagy, mitophagy, autophagosome assembly, and macroautophagy. These processes may be associated with various infections and diseases that were inferred from the enrichment analysis of pathways, including influenza A, lipid dysfunction and atherosclerosis, and cancer. Previous studies have also shown that DN was associated with various autophagy-related biological functions (; ). Then, correlation analysis and the PPI network were constructed to further explore the correlation between the DEARG expression. According to the PPI network, DMNC scores from cytoHubba confirmed that P4HB was a key hub gene in DN tubular injury.

The elevation of renal tubular P4HB expression in DN was verified in other datasets and in vitro experiments. We first verified the elevated gene expression by other databases, and the result was consistent with that of our previous study. The ROC analysis indicated that P4HB expression showed an excellent diagnostic value for DN patients and healthy controls. Through in vitro experiments, we found that the expression of P4HB was significantly elevated after high glucose stimulation of HK-2, accompanied by increased expression of NGAL, a marker for kidney tubular injury, indicating that there was significant damage in HK-2 cells (). In the meantime, we found that P4HB expression was positively correlated with serum creatinine levels and negatively correlated with GFR, and immunohistochemistry staining showed that the expression of P4HB in the renal tubules of DN patients was significantly higher than that of normal kidneys. Therefore, P4HB might serve as an autophagy-related biomarker for DN.

This study can provide novel insights and potential targets for further studies on the connection between DN and autophagy. However, our analysis was limited by the number of samples included, as transcriptomic data on DN tubuleinterstitium are restricted and most data were tested with different platforms. More combined samples and clinical information are required to clarify the potential mechanisms of P4HB in DN.

Conclusion

We identified DEARGs in DN vs. normal human renal tubular interstitial cells by bioinformatics analysis for the first time, and P4HB was found and confirmed as an autophagy-related biomarker for DN, thereby providing new insights and potential targets for further studies on the correlation between DN and autophagy.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below:, GSE30122; further inquiries can be directed to the corresponding author.

Ethics statement

Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

FB carried out the assays, made relevant analysis, and wrote the manuscript. FB, KY, YY, YZ, LD, and XA conducted this research. FF, NS, and JF provided materials and skills and assisted in building cell models. LL, HY, and XY designed the experiment and reviewed and approved the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Grant No. 82070746) and the National Natural Science Fund for Outstanding Youth Fund (Grant No. 82000692).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.965816/full#supplementary-material

References

  • 1

    AhamdS.AlouffiS.KhanS.KhanM.AkashaR.AshrafJ. M.et al (2022). Physicochemical characterization of in vitro LDL glycation and its inhibition by ellagic acid (EA): An in vivo approach to inhibit diabetes in experimental animals. Biomed. Res. Int.2022, 5583298. 10.1155/2022/5583298

  • 2

    AhamdS.MmoinuddinShahabU.HabisbS.KhanM. S.AlamK.AliA. (2014). Glycoxidative damage to human DNA: Neo-antigenic epitopes on DNA molecule could be a possible reason for autoimmune response in type 1 diabetes. Glycobiology24 (3), 281291. 10.1093/glycob/cwt109

  • 3

    BikbovB.PurcellC.LeveyA. S.SmithM.AbdoliA.AbebeM.et al (2020). Global, regional, and national burden of chronic kidney disease, 1990-2017: A systematic analysis for the global burden of disease study 2017. Lancet395 (10225), 709733. 10.1016/s0140-6736(20)30045-3

  • 4

    BonomeT.LevineD. A.ShihJ.RandonovichM.Pise-MasisonC. A.BogomolniyF.et al (2008). A gene signature predicting for survival in suboptimally debulked patients with ovarian cancer. Cancer Res.68 (13), 54785486. 10.1158/0008-5472.Can-07-6595

  • 5

    BrezniceanuM.-L.LauC. J.GodinN.ChenierI.DuclosA.EthierJ.et al (2010). Reactive oxygen species promote caspase-12 expression and tubular apoptosis in diabetic nephropathy. J. Am. Soc. Nephrol.21 (6), 943954. 10.1681/asn.2009030242

  • 6

    DikicI.ElazarZ. (2018). Mechanism and medical implications of mammalian autophagy. Nat. Rev. Mol. Cell Biol.19 (6), 349364. 10.1038/s41580-018-0003-4

  • 7

    DingX.-J.ZhangZ.-Y.JinJ.HanJ.-X.WangY.YangK.et al (2020). Salidroside can target both P4HB-mediated inflammation and melanogenesis of the skin. Theranostics10 (24), 1111011126. 10.7150/thno.47413

  • 8

    FuW.LiuY.LiuF.LiuC.LiJ.NiuJ.et al (2021). A novel autoantibody induced by bacterial biofilm conserved components aggravates Lupus nephritis. Front. Immunol.12, 656090. 10.3389/fimmu.2021.656090

  • 9

    GalluzziL.BaehreckeE. H.BallabioA.BoyaP.PedroJ. M. B.-S.CecconiF.et al (2017). Molecular definitions of autophagy and related processes. Embo J.36 (13), 18111836. 10.15252/embj.201796697

  • 10

    GraysonP. C.EddyS.TaroniJ. N.LightfootY. L.MarianiL.ParikhH.et al (2018). Metabolic pathways and immunometabolism in rare kidney diseases. Ann. Rheum. Dis.77 (8), 12261233. 10.1136/annrheumdis-2017-212935

  • 11

    JangI.PottekatA.PoothongJ.YongJ.Lagunas-AcostaJ.CharbonoA.et al (2019). PDIA1/P4HB is required for efficient proinsulin maturation and ß cell health in response to diet induced obesity.Elife8, e44528. 10.7554/eLife.44528

  • 12

    JohansenK. L.ChertowG. M.FoleyR. N.GilbertsonD. T.HerzogC. A.IshaniA.et al (2021). US renal data system 2020 annual data report: Epidemiology of kidney disease in the United States. Am. J. Kidney Dis.77, A7A8. 10.1053/j.ajkd.2021.01.002

  • 13

    JuW.NairV.SmithS.ZhuL.SheddenK.SongP. X. K.et al (2015). Tissue transcriptome-driven identification of epidermal growth factor as a chronic kidney disease biomarker. Sci. Transl. Med.7 (316), 316ra193. 10.1126/scitranslmed.aac7071

  • 14

    LiA.YiB.HanH.YangS.HuZ.ZhengL.et al (2022). Vitamin D-VDR (vitamin D receptor) regulates defective autophagy in renal tubular epithelial cell in streptozotocin-induced diabetic mice via the AMPK pathway. Autophagy18 (4), 877890. 10.1080/15548627.2021.1962681

  • 15

    LiY.PanY.CaouasasakiS. K.WangY.NiuA.FanXet al (2021). Podocyte EGFR inhibits autophagy through upregulation of rubicon in type 2 diabetic nephropathy. Diabetes70 (2), 562576. 10.2337/db20-0660

  • 16

    LiuL.BaiF.SongH.XiaoR.WangY. Z.YangH. M.et al (2022). Upregulation of TIPE1 in tubular epithelial cell aggravates diabetic nephropathy by disrupting PHB2 mediated mitophagy. Redox Biol.50, 102260. 10.1016/j.redox.2022.102260

  • 17

    LyuL.XiangW.ZhengF.HuangT.FengY.YuanJ.et al (2020). Significant prognostic value of the autophagy-related gene P4HB in bladder urothelial carcinoma. Front. Oncol.10, 1613. 10.3389/fonc.2020.01613

  • 18

    MaezawaY.TakemotoM.YokoteK. (2015). Cell biology of diabetic nephropathy: Roles of endothelial cells, tubulointerstitial cells and podocytes. J. Diabetes Investig.6 (1), 315. 10.1111/jdi.12255

  • 19

    NabiR.AlviS. S.AlouffiS.KhanS.AhmadA.KhanM.et al (2021a). Amelioration of neuropilin-1 and RAGE/matrix metalloproteinase-2 pathway-induced renal injury in diabetic rats by rosuvastatin. Arch. Biol. Sci.73 (2), 265278. 10.2298/abs210316021n

  • 20

    NabiR.AlviS. S.ShahA.ChaturvediC. P.FaisalM.AlatA. A.et al (2021b). Ezetimibe attenuates experimental diabetes and renal pathologies via targeting the advanced glycation, oxidative stress and AGE-RAGE signalling in rats. Arch. Physiol. Biochem.2021, 116. 10.1080/13813455.2021.1874996

  • 21

    NoivaR. (1999). Protein disulfide isomerase: The multifunctional redox chaperone of the endoplasmic reticulum. Semin. Cell Dev. Biol.10 (5), 481493. 10.1006/scdb.1999.0319

  • 22

    ParkerH. S.LeekJ. T.FavorovA. V.ConsidineM.XiaX.ChavanS.et al (2014). Preserving biological heterogeneity with a permuted surrogate variable analysis for genomics batch correction. Bioinformatics30 (19), 27572763. 10.1093/bioinformatics/btu375

  • 23

    RitcheieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res.43 (7), e47. 10.1093/nar/gkv007

  • 24

    SatirapojB. (2018). Tubulointerstitial biomarkers for diabetic nephropathy. J. Diabetes Res.2018, 2852398. 10.1155/2018/2852398

  • 25

    SinghD.FebboP. G.RossK.JacksonD. G.ManolaJ.LaddC.et al (2002). Gene expression correlates of clinical prostate cancer behavior. Cancer Cell1 (2), 203209. 10.1016/s1535-6108(02)00030-2

  • 26

    SliekerR. C.Van Der HeijdenA. A. W. A.SiddiquiM. K.Langendoen-GortM.NijpelsG.HeringsR.et al (2021). Performance of prediction models for nephropathy in people with type 2 diabetes: Systematic review and external validation study. Bmj-British Med. J.374, n2134. 10.1136/bmj.n2134

  • 27

    TabrezS.Al-ShaliK. Z.AhmadS. (2015). Lycopene powers the inhibition of glycation-induced diabetic nephropathy: A novel approach to halt the AGE-RAGE axis menace. Biofactors41 (5), 372381. 10.1002/biof.1238

  • 28

    TangS. C. W.YiuW. H. (2020). Innate immunity in diabetic kidney disease. Nat. Rev. Nephrol.16 (4), 206222. 10.1038/s41581-019-0234-4

  • 29

    WangX.BaiY.ZhangF.YangY.FengD.LiA.et al (2020). Targeted inhibition of P4HB promotes cell sensitivity to gemcitabine in urothelial carcinoma of the bladder. Onco. Targets. Ther.13, 95439558. 10.2147/ott.S267734

  • 30

    WelshJ. B.SapinosoL. M.SuA. I.KernS. G.Wang-rodriguezJ.MoskalukC. A.et al (2001). Analysis of gene expression identifies candidate markers and pharmacological targets in prostate cancer. Cancer Res.61 (16), 59745978.

  • 31

    WoronieckaK. I.ParkA. S. D.MohtatD.ThomasD. B.PullmanJ. M.SusztakK. (2011). Transcriptome analysis of human diabetic kidney disease. Diabetes60 (9), 23542369. 10.2337/db10-1181

  • 32

    WuP.XiangT.WangJ.LvR.MaS.YuanL.et al (2021). Identification of immunization-related new prognostic biomarkers for papillary renal cell carcinoma by integrated bioinformatics analysis. BMC Med. Genomics14 (1), 241. 10.1186/s12920-021-01092-w

  • 33

    XieL.LiH.ZhangL.MaX.DangY.GuoJ.et al (2020). Autophagy-related gene P4HB: A novel diagnosis and prognosis marker for kidney renal clear cell carcinoma. Aging12 (2), 18281842. 10.18632/aging.102715

  • 34

    YangC.ChenX.-C.LiZ.-H.WuH.-L.JingK.-P.HuangX.-R.et al (2021a). SMAD3 promotes autophagy dysregulation by triggering lysosome depletion in tubular epithelial cells in diabetic nephropathy. Autophagy17 (9), 23252344. 10.1080/15548627.2020.1824694

  • 35

    YangH.XieT.LiD.DuX.WangT.LiC.et al (2019). Tim-3 aggravates podocyte injury in diabetic nephropathy by promoting macrophage activation via the NF-kappa B/TNF-alpha pathway. Mol. Metab.23, 2436. 10.1016/j.molmet.2019.02.007

  • 36

    YangM.LuoS.JiangN.WangX.HanY.ZhaoH.et al (2021b). DsbA-L ameliorates renal injury through the AMPK/NLRP3 inflammasome signaling pathway in diabetic nephropathy. Front. Physiol.12, 659751. 10.3389/fphys.2021.659751

  • 37

    YuK.LiD.XuF.GuoH.FengF.DingY.et al (2021). Ido1 as a new immune biomarker for diabetic nephropathy and its correlation with immune cell infiltration. Int. Immunopharmacol.94, 107446. 10.1016/j.intimp.2021.107446

  • 38

    ZhanM.UsmanI. M.SunL.KanwarY. S. (2015). Disruption of renal tubular mitochondrial quality control by myo-inositol oxygenase in diabetic kidney disease. J. Am. Soc. Nephrol.26 (6), 13041321. 10.1681/asn.2014050457

  • 39

    ZhangC.JiangJ.WangL.ZhengL.XuJ.QiX.et al (2020). Identification of autophagy-associated biomarkers and corresponding regulatory factors in the progression of colorectal cancer. Front. Genet.11, 245. 10.3389/fgene.2020.00245

  • 40

    ZhangL.LongJ.JiangW.ShiY.HeX.ZhouZ.et al (2016). Trends in chronic kidney disease in China. N. Engl. J. Med.375 (9), 905906. 10.1056/NEJMc1602469

  • 41

    ZhuZ.HeA.LvT.XuC.LinL.LinJ. (2019). Overexpression of P4HB is correlated with poor prognosis in human clear cell renal cell carcinoma. Cancer Biomark.26 (4), 431439. 10.3233/cbm-190450

Summary

Keywords

autophagy, P4HB, bioinformatics analysis, diabetic nephropathy, biomarker

Citation

Bai F, Yu K, Yang Y, Zhang Y, Ding L, An X, Feng F, Sun N, Fan J, Liu L, Yang H and Yang X (2022) Identification and validation of P4HB as a novel autophagy-related biomarker in diabetic nephropathy. Front. Genet. 13:965816. doi: 10.3389/fgene.2022.965816

Received

10 June 2022

Accepted

25 August 2022

Published

26 September 2022

Volume

13 - 2022

Edited by

Saheem Ahmad, University of Hail, Saudi Arabia

Reviewed by

Xiaoyan Bai, Guangdong Provincial People’s Hospital, China

Hao Zhang, Central South University, China

Updates

Copyright

*Correspondence: Xiangdong Yang,

This article was submitted to Genetics of Common and Rare Diseases, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics