- Research article
- Open Access
- Open Peer Review
Urinary mitochondrial DNA level as a biomarker of tissue injury in non-diabetic chronic kidney diseases
BMC Nephrology volume 19, Article number: 367 (2018)
Urinary mitochondrial DNA (mtDNA) fragment level has been proposed as a biomarker of chronic kidney disease (CKD). In this study, we determine the relation between urinary mtDNA level and rate of renal function deterioration in non-diabetic CKD.
We recruited 102 non-diabetic CKD patients (43 with kidney biopsy that showed non-specific nephrosclerosis). Urinary mtDNA level was measured and compared to baseline clinical and pathological parameters. The patients were followed 48.3 ± 31.8 months for renal events (need of dialysis or over 30% reduction in estimated glomerular filtration rate [eGFR]).
The median urinary mtDNA level was 1519.42 (inter-quartile range 511.81–3073.03) million copy/mmol creatinine. There were significant correlations between urinary mtDNA level and baseline eGFR (r = 0.429, p < 0.001), proteinuria (r = 0.368, p < 0.001), severity of glomerulosclerosis (r = − 0.537, p < 0.001), and tubulointerstitial fibrosis (r = − 0.374, p = 0.014). The overall rate of eGFR decline was − 2.18 ± 5.94 ml/min/1.73m2 per year. There was no significant correlation between the rate of eGFR decline and urinary mtDNA level. By univariate analysis, urinary mtDNA level predicts dialysis-free survival, but the result became insignificant after adjusting for clinical and histological confounding factors.
Urinary mtDNA levels have no significant association with the rate of renal function decline in non-diabetic CKD, although the levels correlate with baseline renal function, proteinuria, and the severity of histological damage. Urinary mtDNA level may be a surrogate marker of permanent renal damage in non-diabetic CKD.
Chronic kidney disease (CKD) is a global public health issue . The prevalence of CKD is estimated to be 8 to 16% worldwide . A considerable proportion of CKD patients eventually progresses to dialysis-dependent end stage renal disease (ESRD), which is an important economic burden to the health care system . Given the increasing prevalence of CKD, the associated burden of CKD is heavy and ever-rising worldwide [3, 4].
CKD is characterized by progressive renal function loss irrespective to the primary cause of kidney damage . However, the underlying mechanism of progressive renal function loss in CKD, it still remains incompletely understood. Emerging evidence shows that mitochondrial dysfunction plays an important role in this regard, especially in the setting of diabetic nephropathy [6,7,8]. Glomerular injury or intra-glomerular hypertension results in proteinuria, which is taken up by renal tubular cells, resulting in the generation of reactive oxygen species and cellular oxidative stress [9, 10]. The end result is mitochondrial damage, which contributes to the tubular cell damage as well as epithelial-mesenchymal transition (EMT) [11, 12]. In experimental models of CKD, aborting mitochondrial dysfunction prevents renal tubular cell EMT  and renal fibrosis .
It has long been known that extracellular cell-free mitochondrial DNA (mtDNA) could easily be detected in various body fluids. Previous studies indicated that mtDNA is released into the systemic circulation after mitochondria is damaging [14, 15]. In addition to blood, dysfunctioned mitochondria in renal tubular cells are also found to release mtDNA into urine, and extracellular mtDNA in urine has recently been explored as the biomarker of various kidney diseases [16,17,18,19,20,21]. Our previous study showed that urinary mtDNA levels are associated with the severity of renal impairment and histological scarring in diabetic kidney disease . In the present study, we determine the relation between urinary mtDNA level and kidney damage in non-diabetic CKD.
This is an observational study approved by the Clinical Research Ethics Committee of the Chinese University of Hong Kong. All study protocols and procedures were in compliance with the Declaration of Helsinki. We recruited patients with non-diabetic CKD. In short, CKD was defined as estimated glomerular filtration rate (eGFR) ≤60 ml/min/1.73m2 or proteinuria ≥1 g/day that persisted for over 3 months. Diabetes was excluded by two fasting plasma glucose tests. A whole-stream early morning urine specimen was collected after written informed consent. Other clinical, biochemical, and histological information were collected by chart review. The eGFR was calculated by the Modification of Diet in Renal Disease (MDRD) formula  and the stage of CKD was defined according to the Kidney Disease Improving Global Outcomes (KDIGO) criteria .
Urine mtDNA preparation
Urine specimens were were processed immediately after collection according to methods described previously . Briefly, protease inhibitors were added, and the samples were centrifuged at 1000-g. The supernatant was then stored at − 80 °C until mtDNA extraction and purification by the QIAamp Blood DNA Midi Kit (Qiagen, Valencia, CA).
Urinary mtDNA level was then measured by digital polymerase chain reaction (PCR) by the chip-based dPCR platform (QuantStudio 3D digital PCR System, Life Technologies, Carlsbad, CA). The method of quantification of urinary mtDNA has been described previously . Briefly, extracted DNA was mixed with digital PCR master mix, assay primers and probe. This PCR mix was then loaded on the digital PCR 20 K chip using the QuantStudio 3D digital PCR Chip Loader. PCR reaction was performed using flat-block ProFlex PCR system according to manufacturer instructions. Each PCR run was performed with a blank control by adding water in place of template. After the thermal cycling, the chips were imaged using the QuantStudio 3D digital PCR Chip Reader and analyzed by the QuantStudio 3D Analysis Suite software according the manufacturers’ instruction.
Morphometric study of kidney biopsy
The degree of renal scarring was assessed by morphometric study on frozen renal biopsy specimen with the Jones’ silver staining. The method of morphometric study has been described previously [27, 28]. In essence, it was a semi-quantitative method using computerized image analysis. The Leica Twin Pro image analysis system (Leica Microsystems, Wetzlar, Germany) was used. Ten glomeruli and 10 randomly selected tubulointerstitial areas were assessed for each patient, and the average percentage of scarred areas were computed.
Recruited patients were followed for at least 24 months. During the follow-up, clinical management of the patients was decided by individual nephrologist and not affected by the study. Serum creatinine, eGFR, and spot urine protein-to-creatinine ratio, were assessed at least every 4 months. The rate of GFR decline was calculated by the least squares regression method. The primary end point was the slope of GFR decline. Secondary outcomes was dialysis-free survival and renal event-free survival. Renal event was defined as over 30% reduction in eGFR as compared to baseline or the need of dialysis.
The software package of SPSS for Windows version 18.0 (SPSS Inc., Chicago, IL) was use to perform the statistical analysis. Continuous data were expressed as mean ± SD or median (inter-quartile range [IQR]) as appropriate. Groups were compared by Chi-square test, independent Student’s t test, Kruskal-Wallis test, or one-way analysis of variance (ANOVA) as appropriate. Spearman’s rank correlation coefficient was use to explore the relation between urinary mtDNA level and clinical or pathological parameters. Urinary mtDNA levels were further sub-grouped to tertiles for further analysis. For renal survival, Kaplan-Meier plot and Cox regression analysis was performed. In addition to urinary mtDNA levels, we used baseline clinical (including eGFR and proteinuria) and histological data (including percentage of glomerulosclerosis and tubulointerstitial fibrosis) as confounding parameters for model construction. These criteria were chosen as covariate because they are well established prognostic markers of CKD. A P value below 0.05 was considered statistically significant. All statistical tests report two-tailed probabilities.
We recruited 102 patients with non-diabetic CKD; 43 had kidney biopsy which showed non-specific nephrosclerosis. The other 59 had CKD of unknown etiology. Table 1 summarizes their demographic, baseline clinical, and biochemical information. All patients received maximal tolerated dose of angiotensin-converting enzyme inhibitor or angiotensin receptor blocker. For the patients with kidney biopsy performed, the median percentage of degree of glomerulosclerosis was 9% (IQR 0 to 55%), and median tubulointerstitial fibrosis 10% (IQR 0 to 30%). There is no statistical significant difference in baseline demographic or clinical data between patients with and without kidney biopsy, except that patients with renal biopsy had a lower serum creatinine and higher urinary mtDNA level than those without biopsy (see Table 1). Mitochondrial DNA was detectable in all urine supernatants, and the median level was 1519.42 million copy per mmol creatinine [IQR 511.81 to 3073.03]. Patients with kidney biopsy had higher urinary mtDNA and lower serum creatinine levels as compared to patients without kidney biopsy (see Table 1).
Relation with baseline parameters
Urinary mtDNA level had a modest but statistically significant correlation with baseline serum creatinine (r = − 0.435, p < 0.001), proteinuria (r = 0.368, p < 0.001), and eGFR (r = 0.429, p < 0.001) (Fig. 1). The correlation was similar between patients with and without kidney biopsy (details not shown). There were also significant correlations between urinary mtDNA level and percentage of glomerulosclerosis (r = − 0.537, p < 0.001) and tubulointerstitial fibrosis (r = − 0.374, p = 0.014).
Relation with renal function decline
The patients were followed for 48.3 ± 31.8 months. The overall rate of eGFR decline was − 2.18 ± 5.94 ml/min/1.73m2 per year. The rate of eGFR decline was similar between patients with and without kidney biopsy (− 1.66 ± 5.01 vs − 2.59 ± 6.61 ml/min/1.73 m2 per year, p = 0.5). There was no significant correlation between urinary mtDNA level and the overall rate of eGFR decline (r = 0.166, p = 0.12). There was no significant difference in the rate of eGFR decline between urinary mtDNA tertiles (Kruskal-Wallis test, p = 0.7) (Fig. 2).
During the study period, 31 patients progress to dialysis dependent renal failure, and another 30 patients had eGFR declined for over 30% from the baseline level. The Kaplan-Meier plot of renal survival is shown in Fig. 3. The dialysis-free survival of urinary mtDNA tertile groups were 56.1, 61.9 and 84.2% at 60 months, respectively (log rank test, p = 0.039). Univariate Cox analysis showed that urinary mtDNA level was a predictor of dialysis-free survival (unadjusted hazard ratio [HR] for each log increase in urinary mtDNA, 0.410, 95%CI 0.233–0.721, p = 0.002). However, after adjusting for clinical and histological confounding factors with multivariate Cox regression analysis, urinary mtDNA level was not an independent predictor of dialysis-free survival. In this model, only baseline eGFR was an independent predictor of dialysis-free survival (adjusted HR 0.881, 95%CI 0.811–0.958, p = 0.003), while urinary mtDNA level had no noticeable effect (adjusted HR 0.804, 95%CI 0.154–4.188, p = 0.8). The result remained similar when renal event-free survival was used for analysis (details not shown).
In this study, we show that mtDNA is easy to detect in the urinary supernatant of patients with non-diabetic CKD. Urinary mtDNA level correlates with the degree of renal dysfunction and histological damage. Although urinary mtDNA levels are associated with renal outcome by univariate analysis, the result was not statistically significant after adjusting for clinical and pathological confounders. The role of urinary mtDNA level as a biomarker of patients with non-diabetic CKD may be limited.
Intrarenal oxidative stress plays a critical role in the initiation and progression of kidney disease . The re-absorbing process in renal tubule is highly energy consuming, so that intra-renal hypoxia is an early manifestation of kidney injury . The end result is mitochondrial dysfunction and damage, which release mtDNA into the extracellular space and systemic circulation. In this regard, Chien et al.  previously showed that diabetic patients had significantly less cellular mtDNA content, and diabetic atherogenesis is associated with decreased amounts of mitochondrial DNA. Mishra et al.  further showed a strong relation between peripheral blood mtDNA damage and diabetic retinopathy, and circulating mtDNA level is a noninvasive biomarker of diabetic retinopathy.
In recent years, urinary cell-free mtDNA level has been explored as a biomarker of kidney diseases. Several studies in the setting of acute kidney injury (AKI) showed that urinary mtDNA level is increased in mice  and human  after renal ischemic injury, and the level of urinary cell-free mtDNA correlates with the duration of ischemia. Published data on urinary cell-free mtDNA level in the setting of CKD are few. Our previous study showed that reduction of intra-renal mtDNA level in diabetic nephropathy correlates with the degree of renal damage , suggesting that intra-renal mitochondrial dysfunction is an important feature of diabetic nephropathy. In another study on biopsy-proved hypertensive nephrosclerosis and IgA nephropathy, we found that urinary mtDNA level correlates with the rate of renal function decline and predicts the risk of doubling of serum creatinine or need of dialysis . In the present study, we extend our observation to non-diabetic kidney diseases with no specific pathological diagnosis.
Our present study showed a significant correlation between urinary mtDNA level and serum creatinine as well as proteinuria. On the first glance, the result may seem counter-intuitive because higher eGFR is expected to be associated with lower proteinuria. However, a lower eGFR probably indicates more permanent renal parenchymal damage and therefore a lower capability to leak mtDNA into the urine. Among the patients with renal biopsy, we observed significant inverse correlations between urinary mtDNA level and the severity of pathological scarring, both in terms of glomerulosclerosis and tubulointerstitial fibrosis. The result is expected because, with increasing fibrosis and declining number of renal parenchymal cells, it is logical to predict that the amount of mtDNA that could be leaked into the urine would reduce. In other words, urinary mtDNA level could be considered as a surrogate marker of the degree of permanent renal parenchymal damage or tissue fibrosis, which also explains its relation with renal survival. Since histological assessment of renal fibrosis is both invasive and subjected to sampling error, our result suggests that urinary mtDNA level may serve as a valuable biomarker of renal damage.
The result of this study is somewhat different from our previous ones [22, 34]. In our previous study on diabetic nephropathy , urinary mtDNA levels inversely correlate with eGFR and positively with tubulointerstitial fibrosis. In another study, urinary mtDNA levels positively correlate with proteinuria and inversely with decline of eGFR as well as renal survival in patients with hypertensive nephrosclerosis and IgA nephropathy . The pattern of correlation of our present study is different from the previous one on diabetic nephropathy . Although the underlying reason is not clear, it is possible that high blood glucose causes metabolic stress and directly triggers the damage to tubular epithelial cells, resulting in the excretion of mitochondrial fragments to the extracellular space and urine. Multiple studies showed that hyperglycemia affects renal tubular cells by several mechanisms: overproduction of reactive oxygen species (ROS), activation of apoptotic pathway, and initiation of autophagy [35,36,37,38,39]. Mitochondria, the powerhouse of all cells, are the center of these events. Urinary mtDNA level is conceivably a surrogate marker of mitochondrial repair and regeneration.
There are other limitations of our study. First, the proportion of patients with biopsy is small, limiting the statistical power of multivariate analysis. There may also be unrecognized selection bias for patients to undergo kidney biopsy. Further studies with a larger sample size from an unselected cohort would be necessary to validate our result. Furthermore, the variability of urine mtDNA level is substantial, which may limit the clinical application of urinary mtDNA for risk stratification of CKD.
In conclusion, our study shows that mtDNA is readily detectable in the urinary supernatant of non-diabetic CKD patients. Urinary mtDNA levels have no significant association with the rate of renal function decline in non-diabetic CKD, although the levels correlate with baseline renal function, proteinuria, and the severity of histological damage. Urinary mtDNA level may be a surrogate marker of permanent renal damage in non-diabetic CKD.
In this study, we find that urinary mtDNA levels have no significant association with the rate of renal function decline in non-diabetic CKD, although the levels correlate with baseline renal function, proteinuria, and the severity of histological damage. Urinary mtDNA level may be a surrogate marker of permanent renal damage in non-diabetic CKD.
Levin A, Tonelli M, Bonventre J, Coresh J, Donner JA, Fogo AB, et al. Global kidney health 2017 and beyond: a roadmap for closing gaps in care, research, and policy. Lancet. 2017;390:1888–917.
Webster AC, Nagler EV, Morton RL, Masson P. Chronic kidney disease. Lancet. 2017;389:1238–52.
GBD 2015 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet. 2016;388:1545–602.
Couser WG, Remuzzi G, Mendis S, Tonelli M. The contribution of chronic kidney disease to the global burden of major noncommunicable diseases. Kidney Int. 2011;80:1258–70.
Bainotti S, Rota E, Bertrero M, Tamburrini O, Balducci A, Formica M. Nephrogenic systemic fibrosis: the first Italian gadolinium-proven case. Clin Nephrol. 2008;70:514–7.
Hallan S, Sharma K. The role of mitochondria in diabetic kidney disease. Curr Diab Rep. 2016;16:61.
Higgins GC, Coughlan MT. Mitochondrial dysfunction and mitophagy: the beginning and end to diabetic nephropathy? Br J Pharmacol. 2014;171:1917–42.
Che R, Yuan Y, Huang S, Zhang A. Mitochondrial dysfunction in the pathophysiology of renal diseases. Am J Physiol Renal Physiol. 2014;306:F367–78.
Yuan YG, Chen Y, Zhang P, Huang SM, Zhu CH, Ding GX, Liu BC, Yang TX, Zhang AH. Mitochondrial dysfunction accounts for aldosterone-induced epithelial-to-mesenchymal transition of renal proximal tubular epithelial cells. Free Radic Biol Med. 2012;53:30–43.
Nangaku M. Mechanisms of tubulointerstitial injury in the kidney: final common pathways to end-stage renal failure. Intern Med. 2004;43:9–17.
Emma F, Montini G, Parikh SM, Salviati L. Mitochondrial dysfunction in inherited renal disease and acute kidney injury. Nat Rev Nephrol. 2016;12:267–80.
Zhu CH, Huang SM, Yuan YG, Ding GX, Chen RH, Liu BC, Yang TX, Zhang AH. Mitochondrial dysfunction mediates aldosterone-induced podocyte damage a Therapeutic Target of PPAR gamma. Am J Pathol. 2011;178:2020–31.
Gong W, Mao S, Yu J, Song J, Jia Z, Huang S, Zhang A. NLRP3 deletion protects against renal fibrosis and attenuates mitochondrial abnormality in mouse with 5/6 nephrectomy. Am J Physiol Renal Physiol. 2016;310:F1081–8.
Zhang Q, Raoof M, Chen Y, Sumi Y, Sursal T, Junger W, Brohi K, Itagaki K, Hauser CJ. Circulating mitochondrial DAMPs cause inflammatory responses to injury. Nature. 2010;464:104–U115.
Oka T, Hikoso S, Yamaguchi O, Taneike M, Takeda T, Tamai T, Oyabu J, Murakawa T, Nakayama H, Nishida K, Akira S, Yamamoto A, Komuro I, Otsu K. Mitochondrial DNA that escapes from autophagy causes inflammation and heart failure. Nature. 2012;485:251–U142.
Nakahira K, Kyung SY, Rogers AJ, Gazourian L, Youn S, Massaro AF, Quintana C, Osorio JC, Wang Z, Zhao Y, Lawler LA, Christie JD, Meyer NJ, Mc Causland FR, Waikar SS, Waxman AB, Chung RT, Bueno R, Rosas IO, Fredenburgh LE, Baron RM, Christiani DC, Hunninghake GM, Choi AM. Circulating mitochondrial DNA in patients in the ICU as a marker of mortality: derivation and validation. PLoS Med. 2013;10:e1001577.
Lee JE, Park H, Ju YS, Kwak M, Kim JI, Oh HY, Seo JS. Higher mitochondrial DNA copy number is associated with lower prevalence of microalbuminuria. Exp Mol Med. 2009;41:253–8.
Zhang Y, Zhao Y, Wen S, Yan R, Yang Q, Chen H. Associations of mitochondrial haplogroups and mitochondrial DNA copy numbers with end-stage renal disease in a Han population. Mitochondrial DNA A DNA Mapp Seq Anal. 2017;28:725–31.
Tin A, Grams ME, Ashar FN, Lane JA, Rosenberg AZ, Grove ML, Boerwinkle E, Selvin E, Coresh J, Pankratz N, Arking DE. Association between mitochondrial DNA copy number in peripheral blood and incident CKD in the atherosclerosis risk in communities study. J Am Soc Nephrol. 2016;27:2467–73.
Wang YC, Lee WC, Liao SC, Lee LC, Su YJ, Lee CT, Chen JB. Mitochondrial DNA copy number correlates with oxidative stress and predicts mortality in nondiabetic hemodialysis patients. J Nephrol. 2011;24:351–8.
Cao H, Ye H, Sun Z, Shen X, Song Z, Wu X, He W, Dai C, Yang J. Circulatory mitochondrial DNA is a pro-inflammatory agent in maintenance hemodialysis patients. PLoS One. 2014;9:e113179.
Wei PZ, Kwan BC, Chow KM, Cheng PM, Luk CC, Li PK, Szeto CC. Urinary mitochondrial DNA level is an indicator of intra-renal mitochondrial depletion and renal scarring in diabetic nephropathy. Nephrol Dial Transplant. 2018;33:784–8.
Levey AS, Bosch JP, Lewis JB, Greene T, Rogers N. Roth D. a more accurate method to estimate glomerular filtration rate from serum creatinine: a new prediction equation. Modification of diet in renal disease study group. Ann Intern Med. 1999;130:461–70.
Kidney Disease Improving Global Outcomes (KDIGO) CKD Work Group. KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney Int Suppl. 2013;3:1–150.
Szeto CC. Urine miRNA in nephrotic syndrome. Clin Chim Acta. 2014;436:308–13.
Whitaker RM, Stallons LJ, Kneff JE, Alge JL, Harmon JL, Rahn JJ, Arthur JM, Beeson CC, Chan SL, Schnellmann RG. Urinary mitochondrial DNA is a biomarker of mitochondrial disruption and renal dysfunction in acute kidney injury. Kidney Int. 2015;88:1336–44.
Bruneval P, Bariety J, Lair MF, Mandet C, Heudes D, Nicoletti A. Mesangial expansion associated with glomerular endothelial cell activation and macrophage recruitment is developing in hyperlipidaemic apoE null mice. Nephrol Dial Transplant. 2002;17:2099–107.
Wang G, Lai FMM, Chow KM, Kwan BCH, Pang WF, Luk CCW, Leung CB, Li PKT, Szeto CC. Urinary mRNA levels of ELR-negative CXC chemokine ligand and extracellular matrix in diabetic nephropathy. Diabetes Metab Res Rev. 2015;31:699–706.
Jha JC, Banal C, Chow BS, Cooper ME, Jandeleit-Dahm K. Diabetes and kidney disease: role of oxidative stress. Antioxid Redox Signal. 2016;25:657–84.
Bhargava P, Schnellmann RG. Mitochondrial energetics in the kidney. Nat Rev Nephrol. 2017;13:629–46.
Chien MC, Huang WT, Wang PW, Liou CW, Lin TK, Hsieh CJ, Weng SW. Role of mitochondrial DNA variants and copy number in diabetic atherogenesis. Genet Mol Res. 2012;11:3339–48.
Mishra M, Lillvis J, Seyoum B, Kowluru RA. Peripheral blood mitochondrial DNA damage as a potential noninvasive biomarker of diabetic retinopathy. Invest Ophthalmol Vis Sci. 2016;57:4035–44.
Hu Q, Ren J, Wu J, Li G, Wu X, Liu S, Wang G, Gu G, Ren H, Hong Z, Li J. Urinary mitochondrial DNA levels identify acute kidney injury in surgical critical illness patients. Shock. 2017;48:11–7.
Wei PZ, Kwan BC, Chow KM, Cheng PM, Luk CC, Lai KB, Li PK, Szeto CC. Urinary mitochondrial DNA level in non-diabetic chronic kidney diseases. Clin Chim Acta. 2018;484:36–9.
Zha D, Cheng H, Li W, Wu Y, Li X, Zhang L, Feng YH, Wu X. High glucose instigates tubulointerstitial injury by stimulating hetero-dimerization of adiponectin and angiotensin II receptors. Biochem Biophys Res Commun. 2017;493:840–6.
Jiao X, Li Y, Zhang T, Liu M, Chi Y. Role of Sirtuin3 in high glucose-induced apoptosis in renal tubular epithelial cells. Biochem Biophys Res Commun. 2016;480:387–93.
Gou R, Chen J, Sheng S, Wang R, Fang Y, Yang Z, Wang L, Tang L. KIM-1 mediates high glucose-induced autophagy and apoptosis in renal tubular epithelial cells. Cell Physiol Biochem. 2016;38:2479–88.
Nie H, Xue X, Liu G, Guan G, Liu H, Sun L, Zhao L, Wang X, Chen Z. Nitro-oleic acid ameliorates oxygen and glucose deprivation/re-oxygenation triggered oxidative stress in renal tubular cells via activation of Nrf2 and suppression of NADPH oxidase. Free Radic Res. 2016;50:1200–13.
Wu Y, Zhang M, Liu R, Zhao C. Oxidative stress-activated NHE1 is involved in high glucose-induced apoptosis in renal tubular epithelial cells. Yonsei Med J. 2016;57:1252–9.
This study was supported in part by Chinese University of Hong Kong research accounts 6901031 and 7101215. The funders of this study do not have any role in study design, data collection, analysis, result interpretation, report writing, or the decision to submit the report for publication. The results presented in this paper have not been published previously in whole or part, except in abstract format.
Availability of data and materials
All study data are available from the corresponding author upon written request.
Ethics approval and consent to participate
This study is approved by the Clinical Research Ethics Committee of the Chinese University of Hong Kong (approval number CRE-2016.480). All patients provided written informed consent to participate in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.