Skip to main content

Bioinformatics-based discovery of the urinary BBOX1 mRNA as a potential biomarker of diabetic kidney disease

Abstract

Background

Diabetic kidney disease (DKD) is the leading cause of end-stage kidney disease (ESKD) in the world. Emerging evidence has shown that urinary mRNAs may serve as early diagnostic and prognostic biomarkers of DKD. In this article, we aimed to first establish a novel bioinformatics-based methodology for analyzing the “urinary kidney-specific mRNAs” and verify their potential clinical utility in DKD.

Methods

To select candidate mRNAs, a total of 127 Affymetrix microarray datasets of diabetic kidney tissues and other tissues from humans were compiled and analyzed using an integrative bioinformatics approach. Then, the urinary expression of candidate mRNAs in stage 1 study (n = 82) was verified, and the one with best performance moved on to stage 2 study (n = 80) for validation. To avoid potential detection bias, a one-step Taqman PCR assay was developed for quantification of the interested mRNA in stage 2 study. Lastly, the in situ expression of the selected mRNA was further confirmed using fluorescent in situ hybridization (FISH) assay and bioinformatics analysis.

Results

Our bioinformatics analysis identified sixteen mRNAs as candidates, of which urinary BBOX1 (uBBOX1) levels were significantly upregulated in the urine of patients with DKD. The expression of uBBOX1 was also increased in normoalbuminuric diabetes subjects, while remained unchanged in patients with urinary tract infection or bladder cancer. Besides, uBBOX1 levels correlated with glycemic control, albuminuria and urinary tubular injury marker levels. Similar results were obtained in stage 2 study. FISH assay further demonstrated that BBOX1 mRNA was predominantly located in renal tubular epithelial cells, while its expression in podocytes and urothelium was weak. Further bioinformatics analysis also suggested that tubular BBOX1 mRNA expression was quite stable in various types of kidney diseases.

Conclusions

Our study provided a novel methodology to identify and analyze urinary kidney-specific mRNAs. uBBOX1 might serve as a promising biomarker of DKD. The performance of the selected urinary mRNAs in monitoring disease progression needs further validation.

Background

Diabetic kidney disease (DKD), one of the most common complications of diabetes mellitus (DM), is the leading cause of end-stage kidney disease (ESKD) [1]. As the prognosis for patients diagnosed with DKD is not optimistic, early detection and risk classification of DKD are of paramount importance for providing appropriate and timely therapies to slow the progression to ESKD.

Currently, the microalbuminuria (MA) test is recognized as the most important diagnostic measure of DKD. However, emerging evidence has shown that MA is not an ideal diagnostic tool. A recent autopsy study demonstrated that patients with diabetes mellitus (DM) may have advanced renal pathological changes before the onset of MA [2]. Over 20% of individuals with type 1 diabetes without MA or macroalbuminuria reached stage 3 chronic kidney disease (CKD) [3]. For these reasons, novel markers which can identify patients with DKD, as well as those with a high risk of developing DKD are urgently needed.

Recently, analysis of urinary mRNAs has emerged as a novel non-invasive method for monitoring a variety of kidney diseases, including CKD [4], renal-allograft rejection [5], renal fibrosis [6] and DKD [7]. In previous studies, the candidate mRNAs were selected from molecules reported to have a role in disease development. This literature-based screening approach, however, is often biased and low-throughput. Moreover, it does not account for the interference from other urinary cells (e.g. bladder cells, inflammatory cells and cancer cells), which imposes considerable limitations when interpreting the origin of urinary mRNAs.

To address the above challenges, a novel bioinformatics-based methodology for selecting urinary mRNAs that more accurately reflect the in situ molecular changes in kidney was established in this study. Then, based on this novel approach, a two-stage study to discover potential biomarkers for DKD was further performed.

Methods

Overall design for the bioinformatics analysis

The overall study design is shown in Fig. 1. Microarray datasets of the following tissues were collected: (1) diabetic kidney tissues (including glomerular and tubular compartments), (2) normal kidney tissues (including glomerular and tubular compartments), (3) normal bladder tissues, (4) bladder cancer (BC) tissues, and (5) inflammatory cells from patients with urinary tract infection (UTI). Next, the transcriptome profiles of these tissues were compared and two types of differentially expressed genes (DEGs) were identified as candidate biomarkers: (1) upregulated genes with (a) highest fold changes (FCs) in diabetic kidney tissues compared with normal kidney tissues and (b) low expression in normal bladder, BC tissues and UTI-derived inflammatory cells; (2) upregulated genes with (a) highest FCs in diabetic kidney tissues compared with normal bladder tissues and (b) low expression in BC tissues and UTI-derived inflammatory cells. Next, using a targeted polymerase chain reaction (PCR) array, the urinary expression of the candidate mRNAs was determined in stage 1 study, and the one with best performance moved on to stage 2 study for validation. To avoid potential detection bias, a one-step PCR assay based on Taqman probes was developed for quantification of the interested mRNA in stage 2 study. Lastly, the in situ expression of the interested mRNA was determined by fluorescent in situ hybridization (FISH) assay and bioinformatics analysis.

Fig. 1
figure 1

The overall study design

Datasets collection

The gene expression datasets used in this study were compiled using the gene expression omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/). For datasets of normal human glomeruli, tubules and bladder tissues, the primary search criteria were set to “glomerulus OR glomeruli”, “tubulus OR tubules” and “bladder” respectively. For DKD, the primary search criteria were set to “diabetic nephropathy OR diabetic kidney disease OR DKD”. For UTI-derived inflammatory cells and BC tissues, the primary search criteria were set to “urinary tract infection OR UTI” and “bladder AND (cancer OR tumor)”, respectively. To minimize platform variations, the platform filtering criterion was set to “Affymetrix U133″, which included three types of microarrays (Human Genome U133 Plus 2.0 Array, Human Genome U133A 2.0 Array, and Human Genome U133A&B). The searching criteria for other types of CKD can be found in our previous study [8].

Processing microarray datasets

Data processing was performed using Bioconductor (version 3.4) in R software (version × 64 3.3.2) according to the following steps [9]:

  1. 1.

    Before included in the analysis, the quality of each dataset was examined using the general quality control (QC) stats in the simpleaffy package and RNA degradation analysis in the affy package. Those unqualified datasets were discarded.

  2. 2.

    For analysis involving different types of microarrays, the 22,215 Affymetrix identifiers shared by both the Human Genome U133 Plus 2.0 Array and Human Genome U133A Array were extracted. Datasets of the same category were pooled to expand the sample size. Then the relative log expression graph was used to evaluate the consistency among datasets, and those with significant bias were discarded.

  3. 3.

    The robust multi-array average method in the affy package was used to preprocess the original data. Then the principal component analysis was used to detect potential batch effects of pooled datasets. Lastly, a linear model was fit to the normalized data to obtain an expression measure for each probe set on qualified datasets.

  4. 4.

    The empirical Bayes method was used to select DEGs for each disease. Statistically significant DEGs were defined as those with p-values < 0.05 after adjustment by the Benjamini–Hochberg method and FCs > 2.

Study population

This 2-stage cross-sectional study was approved by the Ethical Committee of Zhong Da Hospital of Southeast University. All participants provided written informed consent. In the stage 1 study, four categories of participants (n = 82) were enrolled: (1) normoalbuminuric patients with DM (NA group: albumin-creatine ratio [ACR] < 30 mg/g); (2) patients with DKD; (3) patients with UTI (UTI group) or BC (BC group); and (4) healthy controls (HCs). Patients with DKD were further divided into three subgroups: (1) those with MA (MA group: ACR 30–300 mg/g); (2) those with overt albuminuria (OA group: ACR > 300 mg/g); and (3) those with ESKD (ESKD group). In the stage-2 study, another 80 participants were enrolled (NA = 20; MA = 20; OA = 20; HC = 20).

DKD was diagnosed according to the Kidney Disease Outcomes Quality Initiative (KDOQI) guidelines of 2007 (10). The inclusion criteria were set as follows: at least 5 years from diagnosis of type 2 diabetes, the presence of diabetic retinopathy, and an elevated ACR. The exclusion criteria were as follows: infection, signs or symptoms of other systemic diseases, or suspected non-diabetic kidney disease. ESKD was defined as the onset of dialysis or estimated glomerular filtration rate (eGFR) < 15 ml/min/1.73 m2. The eGFR was calculated according to modified modification of diet in renal disease (MDRD) equations [10].

The HCs were enrolled from the Zhong Da Hospital Health Care Center, all of whom met the following criteria: (1) no record of abnormal renal function (eGFR < 90 mL/min/1.73 m2); (2) normal routine urinalysis, ACR, and 24 h urinary protein test results; (3) no record of hypertension, diabetes, hyperlipidaemia, or hyperuricaemia; and (4) no family history of kidney diseases. Moreover, all patients with UTI and BC have no evidence of kidney diseases or diabetes.

PCR assay

In the stage-1 study, a targeted PCR array was fabricated to detect differentially expressed mRNAs with highest FCs. The primer sets were designed using PRIMER 5 software according to optimized experimental conditions. Two housekeeping genes (B2M and RPL27) were used to normalize data. Positive PCR control (PPC) and reverse transcription control (RTC) were set into the array to monitor PCR reaction performance. Cycling conditions were set as follows: 95 °C for 10 min, followed by 40 cycles of 15 s at 95 °C and 60 °C for 1 min.

In the stage-2 study, a Taqman one-step PCR assay was developed for quantification of the BBOX1 mRNA. B2M was used as the housekeeping gene. Cycling conditions were set as follows: 95 °C for 10 min, followed by 45 cycles of 20 s at 95 °C and 60 °C for 45 s. All primer sequences can be found in the additional files (Additional file 1: Table S1). All PCR assays were performed using an ABI PRISM7700 system (Applied Biosystems).

In order to further test the reproducibility and sensitivity of the Taqman PCR assay, the cycle threshold (Ct) values of BBOX1 and B2M in different amounts of total urinary RNA (500 ng, 50 ng, 5 ng and 0.5 ng) were measured for three times. The reproducibility was measured by the coefficient of variation (CV) according to the following formula \(c_{v} = \frac{\sigma }{\mu }\), where σ and μ stand for the standard deviation(SD) and mean of repeated measurements, respectively.

Urinary mRNA samples collection and measurement

First morning urine samples were collected and centrifuged at 3000g for 30 min at 4 °C within 2 h of collection to obtain the urinary sediments. The sediments were then resuspended in 1.5 ml DEPC-treated PBS and centrifuged at 12,000g for 5 min at 4 °C. RNAiso Plus (Takara) was added to preserve total RNA, and the samples were stored at − 80 °C until use. Total RNA was extracted according to the manufacturer’s protocol (Invitrogen). Then, RNA concentrations were measured using a NanoDrop 2000 (Thermo) based on the relative absorbance ratio at 260/280 nm. The qualified RNA samples were then reverse transcribed to cDNA according to the manufacturer’s protocol (Takara), which were stored at − 20 °C until use.

Confirmation of in situ mRNA expression

The Ethical Committee of Zhong Da Hospital of Southeast University approved the use of human samples for the experiments outlined in this study. FISH assay was performed on 2-μm-thick sections of diabetic kidney tissues and normal urothelium to determine the in situ mRNA expression levels. Additionally, kidney tissues were co-stained with podocalyxin antibody to detect the expression of BBOX1 mRNA in podocytes.

Briefly, sections were first deparaffinized and dehydrated in dimethylbenzene and ethanol, followed by rinsing twice in distilled water for 5 min each. After pre-treatment with pepsin and permeabilization, the sections were treated with a BBOX1 gene probe mix (Exiqon, sequence:5′-AGTAA TCCAC TCCAA TGTCT GT-3′) overnight at room temperature. To stain podocytes, the slides were additionally incubated with labeled anti-human podocalyxin monoclonal antibodies(Abcam) at a dilution of 1:100 overnight. Nuclei were counterstained with 4′,6-diamidino-2-phenylindole (DAPI) and coverslips were fixed with nail polish. Analysis of fluorescence signals was performed using a Nikon Eclipse C1 epifluorescence microscope with interference filters (AHF Analysentechnik AG).

Statistical analysis

R software (version ×64 3.3.2) and Graphpad Prism 7.0 were used for all other statistical analysis and figure construction, respectively. The Shapiro–Wilk test was used to determine the normality of the data. Numeric results with a normal distribution were presented as the mean ± SD. Non-normal numeric results were presented with the interquartile range (IQR). Analysis of variance (ANOVA) was applied to compare the means of normalized data. For skewed data, the Kruskal–Wallis test was applied. The frequencies were compared using the Chi squared test. The correlation between gene expression levels and clinical parameters were analyzed using Spearman’s rank-order test. The discriminative power of the biomarker was evaluated by generating receiver operating characteristic (ROC) curves. The area under the curve (AUC) was used to assess the overall discriminatory power. An AUC of 0.6–0.7 was considered poor, 0.7–0.8 was considered moderate, 0.8–0.9 was considered good, and > 0.9 was considered excellent. Optimal cut-offs were determined by selecting the data points that maximized the sum of specificity and sensitivity on the ROC curve. Two-tailed p < 0.05 was considered statistically significant.

Results

Characteristics of microarray datasets

After screening and quality control tests, a total of 127 qualified datasets (DKD glomeruli: n = 7, normal glomeruli: n = 40, DKD tubules: n = 11, normal tubules: n = 22, normal bladder: n = 23, UTI: n = 5, BC: n = 19) were included for candidate screening (Table 1). These datasets were obtained from GSE37463 [11], GSE24152(unpublished study), GSE47185 [12], GSE35489 [13], GSE7476 [14], GSE11783 [15], GSE18810(unpublished study), GSE21785 [16] and GSE20602 [17]. All datasets were generated by either the Human Genome U133 Plus 2.0 Array or Human Genome U133A Array. Of note, all subjects with DKD had biopsy-proven diabetic nephropathy. Normal renal tissues were obtained from either pre-transplant living donors or patients after tumor nephrectomy.

Table 1 Included microarray datasets after quality control

Identification of candidate mRNAs

As shown in Fig. 2, after comparing the gene expression patterns in glomeruli from diabetic kidney tissues and others, a total of 9203 upregulated DEGs were identified (DKD glomeruli vs. normal bladder, n = 6025; DKD glomeruli vs. normal glomeruli, n = 278; DKD glomeruli vs. BC, n = 5083; DKD glomeruli vs. UTI-derived polymorphonuclear leukocytes, n = 8036), of which seventy-eight were co-differentially expressed.

Fig. 2
figure 2

The identification of candidate mRNAs. a The Venn diagram of DEGs identified in glomerular compartments of diabetic kidney tissues; b The Venn diagram of DEGs identified in tubular compartments of diabetic kidney tissues; c A display of selected candidate genes in the PCR array. G, glomerular compartments of diabetic kidney tissues; T, tubular compartments of diabetic kidney tissues; g, glomerular compartments of normal kidney tissues; t, tubular compartments of normal kidney tissues; B, normal bladder tissues; BC, bladder cancer tissues; L, polymorphonuclear leukocytes in urinary tract infection

For tubules from diabetic kidney tissues, a total of 4151 upregulated DEGs were identified (DKD tubules vs. normal bladder, n = 388; DKD tubules vs. normal tubules, n = 28; DKD tubules vs. BC, n = 541; DKD tubules vs. UTI-derived polymorphonuclear leukocytes, n = 4079), of which two were co-differentially expressed.

According to the rules described in the Methods section (Overall design for the bioinformatics analysis), sixteen mRNAs were finally selected as candidate biomarkers. Detailed lists of differentially expressed genes are shown in the additional files (Additional file 2: Table S2, Additional file 3: Table S3, Additional file 4: Table S4, Additional file 5: Table S5, Additional file 6: Table S6, Additional file 7: Table S7, Additional file 8: Table S8, Additional file 9: Table S9).

Characteristics of the study population

A total of 82 participants were enrolled in stage 1 study (HC: n = 14; NA: n = 16; MA: n = 11; OA; n = 12, ESKD: n = 13; UTI: n = 8; BC: n = 8). Table 2 shows the demographic and clinical characteristics of included participants. The mean ages of participants were similar in the HC, NA and UTI group, but those in the MA, NA, ESKD and BC group were older. In addition, the renal function of patients in the MA, OA and ESKD group was significantly poorer than those of HCs. Patients in the UTI and BC group all had normal albuminuria, blood glucose levels and eGFR.

Table 2 Basic characteristics of participants in stage-1 study

In the stage-2 study, another 80 participants were enrolled (NA: n = 20; MA: n = 20; OA: n = 20; HC: n = 20). As shown in Table 3, there was no significant difference in age among different groups. Patients in the OA group, but not the NA or MA group, had lower eGFR than HCs.

Table 3 Basic characteristics of participants in stage-2 study

Urinary BBOX1 mRNA (uBBOX1) expression was upregulated in patients with DKD

To verify the differential expression of candidate mRNAs in the urine, the urinary mRNA profile of each participant of stage 1 study was measured using the self-assembly PCR array. As a result, four urinary mRNAs (BBOX1, CCL18, NPHS2, and SLC3A1) were found to be upregulated in the patients with DKD (MA and OA) compared to HCs (Table 4). After adjusting for false positives using the Benjamini–Hochberg method, uBBOX1 remained significantly different (adjusted p = 0.036). Next, its expression among different groups was further examined in stage 1 study. The median relative expression of uBBOX1 was 0.0072 in the HC group (IQR 0.0072 [0.0027–0.018]). Its expression was increased threefold in the NA group (IQR 0.022 [0.0093–0.059], adjusted p = 0.036) and 3.1-fold in the non-ESKD DKD groups (IQR 0.027 [0.010–0.093], adjusted p = 0.0083), respectively (Fig. 3). Interestingly, its expression was not significantly changed in patients with ESKD compared to the HCs (IQR 0.013 [0.0044–0.024]). However, patients in the NA, MA, and OA group all had similar uBBOX1 expression. In addition, uBBOX1 levels of UTI or BC group were not elevated in comparison with those of HCs, and were significantly lower than those of the NA group (UTI vs. NA: adjusted p = 0.012; BC vs. NA: adjusted p = 0.025) and the non-ESKD DKD groups (UTI vs. MA and OA: adjusted p = 0.0047; BC vs. MA and OA: adjusted p = 0.0086).

Table 4 A list of differentially expressed urinary mRNAs
Fig. 3
figure 3

The differential expression of uBBOX1 among different populations in stage 1 study. a uBBOX1 levels were significantly elevated in the NA, MA and OA group, but not the ESKD group; b uBBOX1 levels of the UTI and BC group were significantly lower than those of the NA, MA and OA group. *p < 0.05; **p < 0.01

Correlation between uBBOX1 expression and clinical parameters

The correlation between uBBOX1 expression and clinical parameters was further analyzed with patients with ESKD excluded. As shown in Fig. 4, the levels of uBBOX1 were positively correlated with albuminuria (Spearman’s r = 0.371, p = 0.0017), urinary N-acetyl-beta-d-glucosaminidase (NAG) levels (Spearman’s r = 0.407, p = 0.0005), blood glucose levels (Spearman’s r = 0.323, p = 0.0069) and HbA1c levels (Spearman’s r = 0.292, p = 0.015). However, significant correlation between uBBOX1 and blood pressure or eGFR was not found.

Fig. 4
figure 4

Correlation between uBBOX1 expression and clinical parameters (patients with ESKD not included). uBBOX1 positively correlated with the levels of urinary ACR (a), urinary NAG (b), blood glucose (c) and HbA1c (d)

The discriminative performance of uBBOX1

Next, ROC analysis was performed to further test the overall accuracy of uBBOX1 in discriminating different groups of participants enrolled in the stage 1 study. As shown in Fig. 5, with HCs as the control group and patients with diabetes and non-ESKD DKD as the test group, uBBOX1 had an AUC of 0.762 (p = 0.0084). At its optimal cut-off value of 0.0082, uBBOX1 yielded a specificity of 63.6% and a sensitivity of 82.0%. When patients with UTI and BC were added to the control group, the AUC for uBBOX1 changed to 0.805 (p < 0.0001). Interestingly, the optimal cut-off value remained at 0.0082, where uBBOX1 yielded a specificity of 66.7% and a sensitivity of 82.1%.

Fig. 5
figure 5

The discriminative performance of uBBOX1. a uBBOX1 yielded an AUC of 0.762 in discriminating the test group (NA, OA and MA) from HC group; b uBBOX1 yielded an AUC of 0.805 in discriminating the test group (NA, OA and MA) from non-diabetes groups (HC, UTI and BC). AUC, area under the curve; 95% CI, 95% confidence interval

Performance validation of uBBOX1 in stage 2 study

To avoid potential detection bias, the reproducibility and sensitivity of the one-step Taqman PCR assay were tested. As shown in Table 5, both BBOX1 and B2M mRNAs can be detected in the reaction systems containing as low as 0.5 ng total urinary RNA/10ul with low variance (average CVs for Ct values of BBOX1: 0.26%; average CVs for Ct values of B2M: 0.14%).

Table 5 Reproducibility and sensitivity of the one-step Taqman PCR assay

Then using this one-step Taqman PCR assay, the uBBOX1 levels of all participants in study 2 were measured. As shown in Fig. 6, the expression of uBBOX1 was significantly upregulated in patients of NA (IQR 0.0064 [0.0018–0.017], adjusted p = 0.0441), MA (IQR 0.0086 [0.0018–0.019], adjusted p = 0.0205) and OA group (IQR 0.027 [0.011–0.058], adjusted p = 0.0013) compared with HCs (IQR 0.0030 [0.0014–0.0074]). However, significant difference of uBBOX1 levels within the NA, MA and OA group was not found. In correlation analysis (Fig. 7), uBBOX1 also positively correlated with urinary ACR levels (Spearman’s r = 0.471, p < 0.0001), urinary NAG levels (Spearman’s r = 0.488, p < 0.0001), blood glucose levels (Spearman’s r = 0.293, p = 0.0084), HbA1c levels (Spearman’s r = 0.263, p = 0.019). A weak negative correlation between eGFR and uBBOX1 was also found (Spearman’s r = 0.263, p = 0.019).

Fig. 6
figure 6

The differential expression of uBBOX1 among different populations in stage 2 study. uBBOX1 levels were significantly elevated in NA, MA and OA group, *p < 0.05; **p < 0.01

Fig. 7
figure 7

Correlation between uBBOX1 expression and clinical parameters (stage 2 study). uBBOX1 positively correlated with the levels of urinary ACR (a), urinary NAG (b), blood glucose (c) and HbA1c (d), while negatively correlated with the levels of eGFR (e)

BBOX1 mRNA in situ expression in diabetic kidney tissues

Next, the in situ expression of BBOX1 mRNA in human tissues was examined. As shown in Fig. 8, BBOX1 mRNA was highly expressed in the tubular compartments of diabetic kidney tissues, while its expression in glomerular compartments including podocytes was weak. Moreover, BBOX1 mRNA was poorly expressed in normal bladder tissues, indicating that detached tubular epithelial cells (TECs) primarily contributed to the elevated uBBOX1 levels. Next, we asked whether the tubular expression of BBOX1 mRNA would change in the context of DKD. As shown in Table 6, bioinformatics analysis involving 33 microarray datasets showed that tubular expression of BBOX1 mRNA was not significantly changed in diabetic kidney tissues (FC: 0.745, adjusted p = 0.136).

Fig. 8
figure 8

In situ expression of BBOX1 mRNA in diabetic kidney tissue and normal bladder tissue. BBOX1 mRNA expression in the glomerular compartment of diabetic kidney tissue (a), the tubular compartment of diabetic kidney tissue (b) and the normal bladder tissue (c). DAPI, 4′,6-diamidino-2-phenylindole; PODXL, podocalyxin

Table 6 Relative expression of BBOX1 mRNA in various types of chronic kidney disease compared with normal kidney tissues

The differential expression of BBOX1 mRNA in other types of CKD

To expand the potential use of uBBOX1, the differential expression of BBOX1 mRNA in other types of CKD was further examined using an integrative bioinformatics approach. The basic characteristics of the included datasets have been described in our previous study [8]. As shown in Table 6, compared with normal kidney tissues, tubular BBOX1 mRNA expression was either not significantly changed or slightly upregulated in other types of CKD including hypertensive nephropathy (relative expression: 0.929, adjusted p = 0.438), focal segmental glomerulosclerosis (relative expression: 1.17, adjusted p = 0.134), IgA nephropathy (relative expression: 0.965, adjusted p = 0.662) and membranous nephropathy (relative expression: 1.28, adjusted p = 0.0045).

Discussion

Although MA is the most common diagnostic measure for DKD, recent studies have suggested that it is suboptimal in many cases. Therefore, there is an urgent need to identify novel biomarkers to increase the accuracy of diagnostic tools.

Previous studies have applied a literature-based approach to screen candidate urinary mRNA markers for kidney diseases. This approach is inherently biased and lacks efficiency. In recent years, bioinformatics has emerged as a powerful tool for the high-throughput identification of potential biomarkers [18]. However, owing to the tremendous amount of data generated by high-throughput technology, bioinformatics analysis based on a few samples may be at a high risk of false positive results. In the present study, we developed a novel bioinformatics workflow based on over one hundred microarray datasets and found that uBBOX1, a TEC-specific mRNA, could be used as a potential biomarker of DKD.

The BBOX1 gene (also known as BBH) is located on chromosome 11 and encodes gamma-butyrobetaine hydroxylase 1, which catalyzes the formation of L-carnitine from gamma-butyrobetaine [19]. L-carnitine deficiency is associated with skeletal myopathies, poorer renal and cardiac function, and anemia in the context of CKD [20, 21]. It also exerts antioxidant and anti-inflammatory effects on TECs in various kidney injury models [22, 23]. The kidney, the liver and the brain are three major organs where carnitine biosynthesis is highly activated [19]. Using quantitative PCR, Rigault et al. found that BBOX1 mRNA is highly expressed in the kidney, with an increase in expression of over threefold and ninefold compared to the liver and brain, respectively [19]. Our result further demonstrated that BBOX1 mRNA is predominantly located in TECs, while its expression in podocytes and urothelium is weak. Therefore, detached TECs and fragments are the major origin of uBBOX1. In addition, BBOX1 mRNA expression, as revealed by our bioinformatics analysis, was not significantly changed in TECs. Hence, detecting uBBOX levels provides a novel and fast method for measuring the extent of TEC loss in the urine.

Previous studies have shown that urinary podocyte mRNAs, which reflect the extent of podocyturia, can be useful as a non-invasive marker for DKD [24, 25]. In a cohort of 1143 patients, urinary podocyte mRNAs were markedly increased in patients with glomerular diseases including DKD and even higher in progressors [26]. Recently, tubulopathy is becoming increasingly recognized as a key culprit of DKD—not only during late stages, but also at the onset of disease [27]. As a consequence of tubular injury induced by DM-related pathophysiological disturbances including albuminuria and hyperglycemia, elevated TEC excretion can be found in the urine of patients with DM [28,29,30]. Hyperglycemia can trigger oxidative stress and lead to aggravated injury of TECs even in non-diabetic rats [31]. Consistent with these studies, we observed that uBBOX1 levels were augmented in normoalbuminuric patients with diabetes and patients with DKD, and positively correlated with albuminuria, hyperglycemia and urinary NAG excretion, indicating its potential diagnostic value for DKD. Its correlation with eGFR, however, was relatively weak. Besides, the differences of uBBOX1 levels among the NA, MA and OA group did not reach statistical significance. Aside from the fact that the relative small sample size limits the statistical power, the sharp increase of urinary mRNA upon kidney damage may also contribute to the insignificance [32]. Another interesting phenomenon is that patients with ESKD had a decreased level of uBBOX1 compared with those in earlier stages, which may have been owing to the lack of kidney intrinsic cells in fibrotic kidneys. Taken together, these data suggest that uBBOX1 may act as an early biomarker for DKD.

To expand the potential use of uBBOX1, we also analyzed its in situ expression in other types of CKD by an integrative bioinformatics method. We found that tubular BBOX1 mRNA expression is quite stable in hypertensive nephropathy and various types of glomerulonephritis, suggesting that elevated uBBOX1 may also act as an indicator of tubular injury in these kidney diseases. Further studies are still required to verify this hypothesis.

Other strengths of this study include: (1) we developed a convenient and effective detection method based on Taqman probes for large-scale clinical verification and transformation; (2) we pre-planned to minimize the potential interference from urinary cells in the context of UTI and BC. Unsurprisingly, uBBOX1 expression was not substantially increased in patients with UTI or BC, supporting the efficacy of this top-down study design. Notably, uBBOX1 was found to be highly effective in discriminating patients with diabetes from HCs and those with confounding diseases.

Some limitations of our study should be considered. Firstly, the study design was cross-sectional with a relatively small sample size, which reduced the significance of the study. Additionally, the diagnosis of DKD in our study was based on laboratory parameters rather than pathological parameters. Large-scale prospective studies are still needed to further validate the prognostic role of uBBOX1 for DKD and other CKDs.

Conclusions

In conclusion, our study provided a novel methodology to identify and analyze urinary kidney-specific mRNAs. Urinary BBOX1 mRNA might serve as a promising biomarker of early detection of DKD. The performance of the selected urinary mRNAs in monitoring DKD progression needs further validation.

Abbreviations

DKD:

diabetic kidney disease

ESKD:

end-stage kidney disease

MA:

microalbuminuria

DM:

diabetes mellitus

CKD:

chronic kidney disease

UTI:

urinary tract infection

BC:

bladder cancer

LUT:

lower urinary tract

DEG:

differentially expressed gene

FC:

fold change

qPCR:

quantitative PCR

PPC:

positive PCR control

RTC:

reverse transcription control

OA:

overt albuminuria

ACR:

albumin-creatine ratio

eGFR:

estimated glomerular filtration rate

HC:

healthy control

FISH:

fluorescent in situ hybridization

DAPI:

4′,6-diamidino-2-phenylindole

ROC:

receiver operating characteristic

AUC:

area under the curve

uBBOX1:

urinary BBOX1

NAG:

N-acetyl-beta-d-glucosaminidase

References

  1. Molitch ME, Adler AI, Flyvbjerg A, Nelson RG, So WY, Wanner C, Kasiske BL, Wheeler DC, de Zeeuw D, Mogensen CE. Diabetic kidney disease: a clinical update from Kidney Disease: improving global outcomes. Kidney Int. 2015;87:20–30.

    Article  Google Scholar 

  2. Klessens CQF, Woutman TD, Veraar KAM, Zandbergen M, Valk EJJ, Rotmans JI, Wolterbeek R, Bruijn JA, Bajema IM. An autopsy study suggests that diabetic nephropathy is underdiagnosed. Kidney Int. 2016;90:149–56.

    Article  Google Scholar 

  3. Molitch ME, Steffes M, Sun WJ, Rutledge B, Cleary P, de Boer IH, Zinman B, Lachin J, Grp ES. Development and progression of renal insufficiency with and without albuminuria in adults with type 1 diabetes in the diabetes control and complications trial and the epidemiology of diabetes interventions and complications study. Diabetes Care. 2010;33:1536–43.

    Article  CAS  Google Scholar 

  4. Szeto CC, Chan RW, Lai KB, Szeto CY, Chow KM, Li PK, Lai FM. Messenger RNA expression of target genes in the urinary sediment of patients with chronic kidney diseases. Nephrol Dial Transplant. 2005;20:105–13.

    Article  CAS  Google Scholar 

  5. Li B, Hartono C, Ding R, Sharma VK, Ramaswamy R, Qian B, Serur D, Mouradian J, Schwartz JE, Suthanthiran M. Noninvasive diagnosis of renal-allograft rejection by measurement of messenger RNA for perforin and granzyme B in urine. N Engl J Med. 2001;344:947–54.

    Article  CAS  Google Scholar 

  6. Zhou LT, Cao YH, Lv LL, Ma KL, Chen PS, Ni HF, Lei XD, Liu BC. Feature selection and classification of urinary mRNA microarray data by iterative random forest to diagnose renal fibrosis: a two-stage study. Sci Rep. 2017;7:39832.

    Article  CAS  Google Scholar 

  7. Wang G, Lai FMM, Lai KB, Chow KM, Kwan CHB, Li KTP, Szeto CC. Urinary mRNA expression of ACE and ACE2 in human type 2 diabetic nephropathy. Diabetologia. 2008;51:1062–7.

    Article  CAS  Google Scholar 

  8. Zhou LT, Qiu S, Lv LL, Li ZL, Liu H, Tang RN, Ma KL, Liu BC. integrative bioinformatics analysis provides insight into the molecular mechanisms of chronic kidney disease. Kidney Blood Press Res. 2018;43:568–81.

    Article  CAS  Google Scholar 

  9. Gao S, Xu S, Fang YP, Fang JW. Prediction of core cancer genes using multi-task classification framework. J Theor Biol. 2013;317:62–70.

    Article  Google Scholar 

  10. Ma YC, Zuo L, Chen JH, Luo Q, Yu XQ, Li Y, Xu JS, Huang SM, Wang LN, Huang W, et al. Modified glomerular filtration rate estimating equation for Chinese patients with chronic kidney disease. J Am Soc Nephrol. 2006;17:2937–44.

    Article  Google Scholar 

  11. Berthier CC, Bethunaickan R, Gonzalez-Rivera T, Nair V, Ramanujam M, Zhang WJ, Bottinger EP, Segerer S, Lindenmeyer M, Cohen CD, et al. Cross-species transcriptional network analysis defines shared inflammatory responses in murine and human lupus nephritis. J Immunol. 2012;189:988–1001.

    Article  CAS  Google Scholar 

  12. Ju WJ, Greene CS, Eichinger F, Nair V, Hodgin JB, Bitzer M, Lee YS, Zhu Q, Kehata M, Li M, et al. Defining cell-type specificity at the transcriptional level in human disease. Genome Res. 2013;23:1862–73.

    Article  CAS  Google Scholar 

  13. Reich HN, Tritchler D, Cattran DC, Herzenberg AM, Eichinger F, Boucherot A, Henger A, Berthier CC, Nair V, Cohen CD, et al. A molecular signature of proteinuria in glomerulonephritis. PLoS ONE. 2010;5:e13451.

    Article  Google Scholar 

  14. Mengual L, Burset M, Ars E, Lozano JJ, Villavicencio H, Ribal MJ, Alcaraz A. DNA microarray expression profiling of bladder cancer allows identification of noninvasive diagnostic markers. J Urol. 2009;182:741–8.

    Article  CAS  Google Scholar 

  15. Gamper M, Viereck V, Geissbuhler V, Eberhard J, Binder J, Moll C, Rehrauer H, Moser R. Gene expression profile of bladder tissue of patients with ulcerative interstitial cystitis. BMC Genomics. 2009;10:199.

    Article  Google Scholar 

  16. Lindenmeyer MT, Eichinger F, Sen K, Anders HJ, Edenhofer I, Mattinzoli D, Kretzler M, Rastaldi MP, Cohen CD. Systematic analysis of a novel human renal glomerulus-enriched gene expression dataset. PLoS ONE. 2010;5:e11545.

    Article  Google Scholar 

  17. Neusser MA, Lindenmeyer MT, Moll AG, Segerer S, Edenhofer I, Sen K, Stiehl DP, Kretzler M, Grone HJ, Schlondorff D, Cohen CD. Human nephrosclerosis triggers a hypoxia-related glomerulopathy. Am J Pathol. 2010;176:594–607.

    Article  CAS  Google Scholar 

  18. Li YP, Li YH, Bai ZJ, Pan J, Wang J, Fang F. Identification of potential transcriptomic markers in developing pediatric sepsis: a weighted gene co-expression network analysis and a case-control validation study. J Transl Med. 2017;15:254.

    Article  Google Scholar 

  19. Rigault C, Le Borgne F, Demarquoy J. Genomic structure, alternative maturation and tissue expression of the human BBOX1 gene. Biochim Biophys Acta. 2006;1761:1469–81.

    Article  CAS  Google Scholar 

  20. Abu Ahmad N, Armaly Z, Berman S, Jabour A, Aga-Mizrachi S, Mosenego-Ornan E, Avital A. l-Carnitine improves cognitive and renal functions in a rat model of chronic kidney disease. Physiol Behav. 2016;164:182–8.

    Article  CAS  Google Scholar 

  21. Guarnieri G. Carnitine in maintenance hemodialysis patients. J Ren Nutr. 2015;25:169–75.

    Article  CAS  Google Scholar 

  22. Li S, Wu W, Wu W, Duan X, Kong Z, Zeng G. L-carnitine protects renal tubular cells against calcium oxalate monohydrate crystals adhesion through preventing cells from dedifferentiation. Kidney Blood Press Res. 2016;41:582–92.

    Article  CAS  Google Scholar 

  23. Alabi QK, Akomolafe RO, Olukiran OS, Nafiu AO, Adefisayo MA, Owotomo OI, Omole JG, Olamilosoye KP. Combined administration of L-carnitine and ascorbic acid ameliorates cisplatin-induced nephrotoxicity in rats. J Am Coll Nutr. 2018;37:387–98.

    Article  CAS  Google Scholar 

  24. Zheng M, Lv LL, Ni J, Ni HF, Li Q, Ma KL, Liu BC. Urinary podocyte-associated mRNA profile in various stages of diabetic nephropathy. PLoS ONE. 2011;6:e20431.

    Article  CAS  Google Scholar 

  25. Zheng M, Lv LL, Cao YH, Liu H, Ni J, Dai HY, Liu D, Lei XD, Liu BC. A pilot trial assessing urinary gene expression profiling with an mRNA array for diabetic nephropathy. PLoS ONE. 2012;7:e34824.

    Article  CAS  Google Scholar 

  26. Wickman L, Afshinnia F, Wang SQ, Yang Y, Wang F, Chowdhury M, Graham D, Hawkins J, Nishizono R, Tanzer M, et al. Urine podocyte mRNAs, proteinuria, and progression in human glomerular diseases. J Am Soc Nephrol. 2013;24:2081–95.

    Article  CAS  Google Scholar 

  27. Gilbert RE. Proximal tubulopathy: prime mover and key therapeutic target in diabetic kidney disease. Diabetes. 2017;66:791–800.

    Article  CAS  Google Scholar 

  28. Morcos M, Sayed AAR, Bierhaus A, Yard B, Waldherr R, Merz W, Kloeting I, Schleicher E, Mentz S, el Baki RF, et al. Activation of tubular epithelial cells in diabetic nephropathy. Diabetes. 2002;51:3532–44.

    Article  CAS  Google Scholar 

  29. Grupp C, John H, Hemprich U, Singer A, Munzel U, Muller GA. Identification of nucleated cells in urine using lectin staining. Am J Kidney Dis. 2001;37:84–93.

    Article  CAS  Google Scholar 

  30. Detrisac C, Mayfield R, Colwell J, Garvin A, Sens D. In vitro culture of cells exfoliated in the urine by patients with diabetes mellitus. J Clin Invest. 1983;71:170.

    Article  CAS  Google Scholar 

  31. Wang JY, Yang JH, Xu J, Jia JY, Zhang XR, Yue XD, Chen LM, Shan CY, Zheng MY, Han F, et al. Renal tubular damage may contribute more to acute hyperglycemia induced kidney injury in non-diabetic conscious rats. J Diabetes Complications. 2015;29:621–8.

    Article  Google Scholar 

  32. Fukuda A, Wickman LT, Venkatareddy MP, Wang SQ, Chowdhury MA, Wiggins JE, Shedden KA, Wiggins RC. Urine podocin:nephrin mRNA ratio (PNR) as a podocyte stress biomarker. Nephrol Dial Transplant. 2012;27:4079–87.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

Author contributions: L-TZ and B-CL designed the research; L-TZ, SQ, T-TT, BW, YF and QY performed experiments; L-TZ and L-LL analyzed the data; L-TZ, L-LL, and K-LM interpreted results of experiments; L-TZ, Z-LL and L-HN prepared figures; L-TZ drafted manuscript; L-LL and B-CL edited and revised manuscript; B-CL approved final version of manuscript. All authors read and approved the final manuscript.

Acknowledgements

None.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Consent for publication

The study was approved by the Ethical Committee of Zhong Da Hospital of Southeast University. Written informed consent was obtained from all the participants and the study was conducted in accordance with the Declaration of Helsinki.

Ethics approval and consent to participate

The study was approved by the Ethical Committee of Zhong Da Hospital of Southeast University. Written informed consent was obtained from all the participants and the study was conducted in accordance with the Declaration of Helsinki.

Funding

This study was supported by National Key R&D Program of China (2018YFC1314002), National Natural Scientific Foundation Grants 81670696, Clinic Research Center of Jiangsu Province Grant BL2014080, and Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX17_0173).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bi-Cheng Liu.

Additional files

Additional file 1: Table S1.

Primer sequences of candidate mRNAs.

Additional file 2: Table S2.

Identified DEGs between diabetic glomeruli and normal bladder.

Additional file 3: Table S3.

Identified DEGs between diabetic glomeruli and normal glomeruli.

Additional file 4: Table S4.

Identified DEGs between diabetic glomeruli and UTI-derived inflammatory cells.

Additional file 5: Table S5.

Identified DEGs between diabetic glomeruli and bladder cancer.

Additional file 6: Table S6.

Identified DEGs between diabetic tubules and normal bladder.

Additional file 7: Table S7.

Identified DEGs between diabetic tubules and normal tubules.

Additional file 8: Table S8.

Identified DEGs between diabetic tubules and UTI-derived inflammatory cells.

Additional file 9: Table S9.

Identified DEGs between diabetic tubules and bladder cancer.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, LT., Lv, LL., Qiu, S. et al. Bioinformatics-based discovery of the urinary BBOX1 mRNA as a potential biomarker of diabetic kidney disease. J Transl Med 17, 59 (2019). https://doi.org/10.1186/s12967-019-1818-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-1818-2

Keywords