Back to Journals » Clinical Epidemiology » Volume 11

Tumor size and postoperative kidney function following radical nephrectomy

Authors Ellis RJ , White VM , Bolton DM, Coory MD, Davis ID , Francis RS, Giles GG , Gobe GC, Neale RE, Wood ST, Jordan SJ 

Received 12 December 2018

Accepted for publication 2 March 2019

Published 6 May 2019 Volume 2019:11 Pages 333—348

DOI https://doi.org/10.2147/CLEP.S197968

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 3

Editor who approved publication: Professor Irene Petersen



Robert J Ellis,1–4 Victoria M White,5,6 Damien M Bolton,7,8 Michael D Coory,8 Ian D Davis,9,10 Ross S Francis,2–4 Graham G Giles,5,8 Glenda C Gobe,3,4 Rachel E Neale,1 Simon T Wood,3–4,10 Susan J Jordan1,3

1Department of Population Health, QIMR Berghofer Medical Research Institute, Brisbane, QLD, Australia; 2Department of Nephrology, Princess Alexandra Hospital, Brisbane, QLD, Australia; 3Faculty of Medicine, University of Queensland, Brisbane, QLD, Australia; 4Kidney Disease Research Collaborative, Translational Research Institute, Brisbane, QLD, Australia; 5Cancer Council Victoria, Melbourne, VIC, Australia; 6School of Psychology, Deakin University, Melbourne, VIC, Australia; 7Department of Urology, Austin Health, Melbourne, VIC, Australia; 8Faculty of Medicine, Dentistry and Health Sciences, University of Melbourne, Melbourne, VIC, Australia; 9Eastern Health Clinical School, Monash University and Eastern Health, Melbourne, VIC, Australia; 10Department of Urology, Princess Alexandra Hospital, Melbourne, VIC, Australia

Background: Chronic kidney disease (CKD) following nephrectomy for kidney tumors is common, and both patient and tumor characteristics may affect postoperative kidney function. Several studies have reported that surgery for large tumors is associated with a lower likelihood of postoperative CKD, but others have reported CKD to be more common before surgery in patients with large tumors.
Objective: The aim of this study was to clarify inconsistencies in the literature regarding the prognostic significance of tumor size for postoperative kidney function.
Study design and setting: We analyzed data from 944 kidney cancer patients managed with radical nephrectomy between January 2012 and December 2013, and 242 living kidney donors who underwent surgery between January 2011 and December 2014 in the Australian states of Queensland and Victoria. Multivariable logistic regression was used to assess the primary outcome of CKD upstaging. Structural equation modeling was used to evaluate causal models, to delineate the influence of patient and tumor characteristics on postoperative kidney function.
Results: We determined that a significant interaction between age and tumor size (P=0.03) led to the observed inverse association between large tumor size and CKD upstaging, and was accentuated by other forms of selection bias. Subgrouping patients by age and tumor size demonstrated that all patients aged ≥65 years were at increased risk of CKD upstaging, regardless of tumor size. Risk of CKD upstaging was comparable between age-matched living donors and kidney cancer patients.
Conclusion: Larger tumors are unlikely to confer a protective effect with respect to postoperative kidney function. The reason for the previously reported inconsistency is likely a combination of the analytical approach and selection bias.

Keywords: renal cell carcinoma, glomerular filtration rate, selection bias, tumor size, kidney cancer, living kidney donors

Introduction

In Australia, despite increased uptake up nephron-sparing approaches, radical nephrectomy remains the most common management approach for kidney masses suspicious of malignancy.1 Nephron mass reduction is associated with an increased risk of subsequent kidney functional decline and development of chronic kidney disease (CKD). In patients managed with radical nephrectomy, a complex interplay between tumor and patient characteristics determines the pre-to-postoperative change in kidney function. This interplay is not well described.

Patient characteristics associated with postoperative kidney function are those related to the overall health and functional reserve of the kidney, predominantly patient age and comorbidities, which cause damage to the kidneys.2,3 The associated bilateral chronic fibrotic cortical damage, with concomitant nephron dysfunction, reduces the ability of the remaining nephrons to compensate for increased functional demand following resection of the affected kidney.2,3

The role of tumor characteristics in determining postoperative kidney function most likely relate to the influence of the tumor on preoperative kidney function, as a consequence of the physical effects of tumor expansion into kidney parenchyma. An expanding tumor is likely to cause dysfunction and eventual loss of adjacent nephrons. Consequently, the glomerular filtration rate (GFR) of the affected kidney is likely to decline. This nephron loss is suggested to generate functional demand, which may be compensated for by increased single-nephron GFR of the remaining nephrons.4 Tumor size and complexity are likely determinants of this, and some studies have reported that risk of postoperative CKD is greater for patients with small tumors managed with radical nephrectomy, compared with larger tumors.5,6 Conversely, preoperative kidney function has been reported to be worse in patients with larger tumors.5,7 This presents a paradox: patients with larger tumors are at increased risk of CKD before nephrectomy, but apparently decreased the risk of CKD after nephrectomy.

The objective of this study was, therefore, to comprehensively investigate the association between tumor size and kidney function following radical nephrectomy, using data from population-based samples of people undergoing nephrectomy for renal cell carcinoma (RCC) or for kidney donation, in order to resolve the apparent paradox stated above. Living kidney donors were included as a comparison group, as their kidney function is not affected by a preexisting tumor and they are generally free from significant underlying health conditions affecting kidney function. We hypothesized that associations would vary according to patient age because of the strong associations between age and CKD,8 and evidence that compensation for nephron mass reduction is modified by age.9 We also aimed to thoroughly evaluate potential sources of bias in our results.

Patients and methods

Study population

The present analysis included data from an Australian kidney cancer patterns of care study, which captured relevant information on all patients diagnosed with incident RCC between January 2012 and December 2013 who were residents of the Australian states of Queensland and Victoria.10 There were 2,323 patients aged ≥18 years who were notified to either Queensland or Victorian Cancer Registries with newly diagnosed RCC. Patients were excluded if they: were not managed surgically (n=229), underwent partial nephrectomy (n=657), had missing values for preoperative (n=270) or postoperative (n=119) estimated GFR (eGFR), or had a preoperative eGFR <45 mL/min per 1.73 m2 (n=104). After exclusion criteria were applied, data for 944 patients were available for analysis (Figure 1).

Figure 1 Study participants. Flow diagram demonstrating the application of exclusion criteria in study cohorts.Abbreviations: CKD, chronic kidney disease; eGFR, estimated glomerular filtration rate; RCC, renal cell carcinoma.

We included a sample of living kidney donors as a comparison group. All living donors notified to the Australia and New Zealand Dialysis and Transplant (ANZDATA) registry who donated a kidney in Queensland or Victoria between January 2011 and December 2014 were identified (n=461). Of these, 242 patients had data on 12-month follow-up kidney function available. Missing data were presumed to be missing at random.

Ethical considerations

Approval to access medical records was obtained under the Queensland Public Health Act and the human research ethics committees of Cancer Council Victoria, QIMR Berghofer Medical Research Institute, and Metro South Hospital and Health Service. Accessing ANZDATA records was designated low risk and approved by the University of Queensland.

Exposures

The Queensland and Victorian Cancer Registries provided sufficient identifying information to study staff to access medical records of patients with RCC. Data on age, sex, serum creatinine concentration, comorbidities; and tumor size, location, and histology for patients managed with RCC, were abstracted. Tumor size was recorded as the largest diameter in the pathology report, and grouped as ≤70 mm and >70 mm (the upper bound for classifying T1 tumors).11 Age was categorized as <65 and ≥65 years, as this threshold has prognostic significance in patients with CKD.12 We calculated eGFR using the CKD-EPI equation.13 To evaluate comorbidity burden, a Charlson comorbidity index was calculated for RCC patients, excluding parameters relating to the RCC;14 this was categorized as low-medium (0–1) and high (≥2). As albuminuria measurements and tumor complexity data were not available for most RCC patients, they were not considered. Demographic data for living donors were supplied by ANZDATA. Living donors were assumed to fall into the low-medium category for comorbidities, due to extensive medical screening and exclusions applied to potential kidney donors in Australia.15

Outcomes

Postoperative kidney function was recorded at approximately 12 months after surgery for the majority of patients. We considered postoperative eGFR as a continuous outcome and CKD upstaging as a categorical outcome. CKD upstaging was defined as new-onset postoperative eGFR <60 or <45 mL/min per 1.73 m2 (if preoperative eGFR was ≥60 or ≥45 mL/min per 1.73 m2, respectively), corresponding to new-onset stage 3a or 3b (or greater) CKD.16

Statistical analysis

Patient and tumor characteristics, and exclusion criteria, were first compared descriptively between groups. To replicate the approach used by previous studies, we then used univariable and multivariable logistic regression to investigate the association between tumor size and CKD upstaging in RCC patients. Multivariable models were adjusted only for potential confounders (not mediators), which we identified using directed acyclic graphs. The association between other patient and tumor characteristics and these outcomes was also evaluated in a similar fashion, in order to contextualize results. As previous studies have demonstrated that CKD upstaging after nephrectomy is relatively common,10 we also used log-binomial regression to estimate relative risks; however, convergence was not achievable in some multivariable models, so more parsimonious models were used. We were therefore not able to use this as our main analytic strategy.

Previous studies that have evaluated the association between tumor size and postoperative kidney function adjusted for preoperative eGFR,5,6 which we considered to be a potential mediator of this association. As a sensitivity analysis, we developed an additional model which adjusted for preoperative eGFR.

CKD upstaging may not be the best measure to investigate the relationship between tumor size and postoperative CKD, as many patients with larger tumors have already experienced significant preoperative reduction in kidney function due to tumor expansion, and thus nephrectomy may result in only minimal further decline. To address this, we also compared groups (those with large and small tumors) using linear regression analysis, considering postoperative eGFR as the outcome.

As certain patient characteristics (eg, age, sex, comorbidity burden, and preoperative eGFR) may affect the ability of the kidney to compensate for surgical removal of nephrons, we also assessed whether there was an interaction between these parameters and tumor size, by adding a first-order multiplicative interaction term to logistic regression models. A significant interaction between age and tumor size was noted. To investigate this further, patients were subgrouped by clinically relevant thresholds for both age and tumor size. We used logistic and linear regression models to compare the association between subgroup membership and CKD upstaging, or postoperative eGFR, respectively. We then added living kidney donors grouped by the same age thresholds (designating donors aged <65 years the reference group) in order to see how the results were affected by comparisons with patients who did not have a kidney tumor. There were 129 RCC patients aged >75 years that were not included in this analysis, as there were no comparable living donors of this age. Models were initially adjusted only for confounders, and subsequently also adjusted for preoperative eGFR. Two additional sensitivity analyses were conducted evaluating postoperative eGFR using linear regression analysis, considering broadened inclusion criteria. Patients with a preoperative eGFR >30 mL/min per 1.73 m2, and patients with any preoperative eGFR were included in these two subsequent analyses.

As we were particularly interested in the relationship between tumor and patient characteristics, and postoperative eGFR (and whether tumor size directly leads to lower risk of postoperative CKD), we also used structural equation modeling to investigate and visualize this relationship. Because the effect of tumor size on postoperative eGFR at 12 postoperative months would either work indirectly through its influence on preoperative eGFR, or through both direct and indirect pathways, we considered two prespecified models in this analysis. The first assumed there was a direct causal path between tumor size and postoperative eGFR; the second did not consider this pathway, instead assuming that the only causal path between tumor size and postoperative eGFR was indirect, via preoperative eGFR. The relationship between variables was compared across the two models.

All analysis was performed using Stata 13.0 (Stata Corp, College Station, TX, USA).

Results

Clinical characteristics

Characteristics of the study population are presented in Table 1 and Table S1. Exclusion criteria for RCC patients, compared across subgroups, are presented in Table S2. Kidney function was recorded at a median (interquartile range) follow-up time of 12.0 (9.8–13.3) and 12.1 (11.0–13.2) months for patients with RCC and living kidney donors, respectively. Donors were younger on average, with higher preoperative eGFR than RCC patients. The male-to-female ratio of donors was approximately equal, whereas there was a male predominance in patients with RCC. When compared on tumor size, patients with tumors ≥70 mm generally had lower preoperative eGFR and smaller pre-to-postoperative ΔeGFR. When comparing postoperative eGFR between patients subgrouped by age and indication/tumor size, older patients tended to have lower postoperative eGFR values, but there were no major differences in the postoperative eGFR of patients of similar age grouped by tumor size (Figure 2).

Table 1 Clinical characteristics of patients grouped by indication/tumor size, and age

Figure 2 Comparison of (A) pre- and (B) postoperative estimated glomerular filtration rate (eGFR) by patients grouped by age and tumor size/indication. Box and whisker plot of pre- and postoperative eGFR, with patients subgrouped by age and indication. Postoperative kidney function was recorded at approximately 12 months after surgery for the majority of patients.

Associations with CKD upstaging

Tumor size was the only tumor characteristics associated with CKD upstaging; tumors ≥70 mm were inversely associated (adjusted odds ratio [aOR]: 0.6, 95% CI: 0.4–0.8) when adjustment was made only for potential confounders (Table 2). After adjustment was made for preoperative eGFR (mediator), this inverse association became stronger (aOR: 0.4, 95% CI: 0.3–0.6).

Table 2 Associations between patient and tumor characteristics, and chronic kidney disease (CKD) upstaging in 944 patients who underwent radical nephrectomy for kidney tumors, considering the interaction between tumor size and patient characteristics

The patient characteristics most strongly associated with CKD upstaging were age (aOR per 5-year increase: 1.5, 95% CI: 1.4–1.6), male sex (aOR: 1.4, 95% CI: 1.0–1.9), and preoperative eGFR (aOR per 5-unit decrease: 1.2, 95% CI: 1.2–1.3) (Table 2).

Interaction between patient age and tumor size

There was a significant interaction between age and tumor size (P=0.03). Compared with patients aged <65 years with tumors <70 mm, patients ≥65 years with tumors <70 mm were at increased risk of CKD upstaging (aOR: 5.2, 95% CI: 3.5–7.5), as were patients aged ≥65 years with tumors ≥70 mm (aOR: 1.9, 95% CI: 1.0–3.1). There was no statistically significant difference in the risk of CKD upstaging in patients aged <65 years with tumors ≥70 mm (aOR: 0.8, 95% CI: 0.5–1.1). The association between tumor size and postoperative kidney function stratified by age is presented in Table S3.

Comparisons with living kidney donors

To investigate interactions between age and tumor size further, analyses were performed including living kidney donors, considering donors aged <65 years as the reference group, as by definition postoperative kidney function in these patients was not impacted by tumor factors (Figure 3A; Table 3). Compared with younger donors, patients aged <65 years with tumors of any size had similar risk of CKD upstaging. Donors aged ≥65 years and patients in the same age group with smaller tumors were at similarly increased risk of CKD upstaging compared with donors aged <65 years (aOR: 6.2, 95% CI: 2.5–15.0; and aOR: 6.1, 95% CI: 3.7–10.0, respectively). Patients with larger tumors aged ≥65 years were also at increased risk of CKD upstaging, but to a lesser degree (aOR: 1.9, 95% CI: 1.0–3.5). After adjustment for preoperative eGFR, the association between tumor size and CKD upstaging became inverse for patients with tumors ≥70 mm, while associations for patients with smaller tumors remained similar to age-matched living kidney donors (Figure 3B). Log-binomial regression models estimating relative risk rather than odds ratios showed similar patterns to the logistic regression analysis, but with substantially smaller differences in estimates between the groups of patients aged ≥65 years.

Table 3 Comparisons of postoperative kidney function between donor and tumor nephrectomy, grouped by age and tumor size

Figure 3 Forest plot showing odds of chronic kidney disease upstaging in patients grouped by tumor size/indication. (A) Forest plot showing associations between patients grouped by indication/tumor size and age, with adjustment made only for potential confounders (sex and Charlson comorbidity index). (B) The same model as (A), with adjustment also made for preoperative estimated glomerular filtration rate (eGFR). The estimates remain relatively similar for all groups, except large tumors, where the effect size reverses following adjustment for eGFR. Postoperative kidney function was recorded at approximately 12 months after surgery for the majority of patients.

Linear regression analyses evaluating postoperative eGFR demonstrated a similar pattern of results to the logistic regression analysis, where the β coefficient was less negative for patients with larger tumors who were aged ≥65 years, compared with both age-matched patients with smaller tumors and living donors (Table 3). A sensitivity analysis was conducted, broadening inclusion to patients with a preoperative eGFR >30 mL/min per 1.73 m2, and any preoperative eGFR. This resulted in a smaller difference between the β coefficient for older patients with smaller compared with larger tumors (percentage difference narrowed from 110% to 40%) (Table S4).

Path analysis

To further investigate associations between patient age/tumor size, and postoperative eGFR, we used structural equation modeling to see how this relationship varied when considering indirect causal pathways directed through preoperative eGFR. We considered two models. The first assumed a direct causal pathway between tumor size and postoperative eGFR, as well as an indirect path mediated through preoperative eGFR. When considering direct effects, tumor size was negatively associated with preoperative eGFR and positively associated with postoperative eGFR (Figure S1A). When considering total effects, larger tumors were associated with a higher postoperative eGFR (β: 0.5, 95% CI: 0.1, 0.9; Table S5). The second model assumed no causal pathway between tumor size and postoperative eGFR. The direct effect of tumor size on preoperative eGFR was essential of the same magnitude as in the previous model (Figure S1B) but the total effect was reversed (β: −0.4, 95% CI: −0.6, −0.2), such that larger tumors were associated with lower postoperative eGFR. Age was negatively associated with pre- and postoperative eGFR in both models, and the magnitude of this estimate was essentially unchanged. Sensitivity analyses were conducted to determine if the interaction between age and tumor size affected these models; no major deviations were observed.

Discussion

Our goal was to clarify whether large tumors were associated with lower risk of CKD upstaging and better postoperative kidney function following radical nephrectomy for RCC. This was to address an apparent contradiction in the literature: that patients with large tumors are at higher risk of CKD before nephrectomy,7 but apparently decreased the risk of CKD postoperatively.5 We attempted to resolve this through causal modeling, concluding that large tumors do not confer a protective effect, and that it may not be appropriate to consider a direct causal pathway between tumor size and postoperative eGFR. We have proposed a number of interrelated contributing biases, which may have led to the reported negative association between large tumors and postoperative CKD.

In our cohort, we initially replicated the apparent paradox that patients with larger tumors tended to have lower preoperative eGFR values (Table 1), yet larger tumors were inversely associated with CKD upstaging (Table 2), which is consistent with other reports. In a retrospective cohort study of 271 Japanese patients, those with tumors <70 mm were at increased risk of new-onset CKD compared with patients with tumors ≥70 mm.6 This was also observed in a retrospective study of 1,371 patients managed at a single Korean center, where, compared with having a tumor >70 mm, having smaller tumors (<40 mm and 40–70 mm) was associated with a greater likelihood of new-onset CKD (aOR: 2.4, 95% CI: 1.6–3.6; and 2.2, 95% CI: 1.5–3.4, respectively).5,6 These authors also showed that larger tumors were associated with lower preoperative eGFR,5 a finding supported by results of a study of 1,569 RCC patients from the United States, which showed that 52% of patients with a tumor ≥70 mm had CKD before surgery (eGFR <60 mL/min per 1.73 m2 and/or urine albumin >30 mg/dL), compared with 40% and 47% of patients with tumors <40 and 40–70 mm, respectively.7

We then showed that, when evaluating the association between tumor size and postoperative CKD upstaging/eGFR, this association varied according to the age of the patient. Additional analyses confirmed that there was no statistically significant difference in the risk of CKD upstaging in all patients with RCC younger than 65 years, regardless of tumor size; and that this risk was comparable to living kidney donors in the same age group. This potentially represents a practical example of Simpson’s paradox, which occurs when there is an interaction between two exposures in a nonrandomized cohort that can potentially lead to incongruent findings between aggregate and disaggregate results.17,18 Although a total reversal of effect was not seen in our data (Table S3), when considering the comparisons between living kidney donors and patients with kidney cancer, the difference was great enough to suggest that larger tumors were not associated with a lower risk of CKD upstaging in younger patients (Figure 3A), and that aggregate results were likely influenced by this interaction.

We did identify that the positive association between older age and CKD upstaging was not as pronounced in patients with larger tumors (Figure 3A). This does not seem biologically consistent, as older patients have limited capacity to compensate to nephron reduction, and if tumor size was driving contralateral compensation, it would be expected that this effect would be more pronounced in the younger age group.9,19,20 Interestingly, we did note this counterintuitive effect was less obvious when comparing relative risk to odds ratio, which could indicate that the odds ratio was exaggerated discrepantly between subgroups when quantifying effect size (Table 3).21 A substantial part of the reason for this counterintuitive effect is that older patients with larger tumors were more likely to experience declines in kidney function before surgery, and therefore already experienced CKD upstaging before undergoing nephrectomy. This could be a consequence of both the fact that older patients generally tend to have a lower eGFR, and because larger tumors probably caused reductions in kidney function prior to surgery, due to secondary nephron loss. Thus, surgical removal of the affected kidney had very little effect on postoperative kidney function. This hypothesis is supported by the observation that patients in this subgroup had the smallest pre-to-postoperative ΔeGFR compared with any other subgroup (Table 1), and the fact that 27% of older patients with larger tumors had a preoperative eGFR <60 mL/min per 1.73 m2, compared with 16% of older patients with small tumors. Another potential contributor to this effect is selection bias, introduced because a higher proportion of patients were excluded because they had a preoperative eGFR <45 mL/min per 1.73 m2 in the subgroup of older patients with larger tumors. This type of selection bias was first reported by Berkson, and occurs when an exposure has an association with both the outcome of interest and the likelihood of a patient being included in a study/subgroup.22 To partially address this, we evaluated postoperative eGFR using linear regression analysis, expanding inclusion to all patients and those with a preoperative eGFR >30 mL/min per 1.73 m2. We found the estimates for older patients with small and large tumors were closer in value. This supports the assertion that larger tumors in older patients do not reduce the risk of CKD; notwithstanding, this cannot be stated definitively due to limitations of this study, and it merits further investigation. It is also possible that this finding is a consequence of survivor-treatment bias, where patients with multiple comorbidities affecting kidney function were more likely to die before developing or being diagnosed with RCC in the older group, in contrast to the younger group, where these patients would more likely have been included.23 We did not address survivor-treatment bias in this study.

As other studies evaluating the association between tumor size and postoperative CKD adjusted for preoperative eGFR in their analysis, this may have contributed to the strong associations reported in the literature.5,6 We also noted that adjustment for preoperative eGFR in multivariable models caused a stronger inverse association between tumor size and CKD upstaging/postoperative eGFR, and adjustment reverses the direction of the effect in both subgroups of patients with large tumors in the logistic regression analysis. We suggest that previous findings could have been contributed to by collider-stratification bias, a type of selection bias that occurs when conditioning on an exposure (preoperative eGFR) that can have two or more common causes (eg, tumor size, age, and a variety of unmeasured patient factors).24 As a possible explanation, consider that preoperative eGFR is influenced by two reasonably independent groups of variables: patient- and tumor-derived characteristics. Patient characteristics (eg, older age) cause reductions in preoperative eGFR due to chronic, bilaterally symmetric pathological changes to the kidney.3 Conversely, tumor characteristics (eg, tumor size) exert their effect predominantly on preoperative eGFR through ipsilateral nephron reduction, which has a null effect on the function of the contralateral kidney, or (if compensation is present) increases the function of the contralateral kidney before nephrectomy. When conditioning on preoperative eGFR while investigating the effect of tumor size, a spurious backdoor path between tumor size and unmeasured patient characteristics (eg, various unmeasured comorbidities which are associated with CKD) may have been generated, leading to biased estimates (Figure 4).25 Another potential cause for the effect reversal in the logistic regression models is that the mediating effect of preoperative eGFR was exaggerated due to the commonality of the outcome, which explains why this effect was resolved when a log-binomial model was used.26

Figure 4 Potential role of collider-stratification bias. This directed acyclic graph (DAG) depicts the hypothesized causal relationship between preoperative estimated glomerular filtration rate (eGFR) and postoperative eGFR, confounded by both tumor size and other risk factors for chronic kidney disease (CKD). In this model, the relationship between tumor size and postoperative eGFR is shown to be mediated by preoperative eGFR, as larger tumors tend to cause preoperative reductions in kidney function. Physiologically, it would be expected that the direct effect of tumor size on postoperative eGFR is quite small in magnitude, because once a tumor has been excised, it should not continue to influence kidney function. Also depicted in this model are other risk factors for CKD, most of which were unmeasured, which would cause reductions in both pre- and postoperative eGFR. Unlike tumor size, other risk factors for CKD will probably lead to ongoing deterioration in kidney function. Therefore, preoperative eGFR becomes a collider in this DAG. When evaluating the association between tumor size and postoperative eGFR, adjusting for this collider may result in a biased estimate, because a spurious causal pathway is opened (Tumor Size → Preoperative eGFR ← CKD Risk Factors → Postoperative eGFR). This becomes a problem because an artificial comparison is generated. The existence of a low preoperative eGFR can be caused by a large tumor, CKD risk factors, or both; however, if a patient has a low preoperative eGFR caused by a growing tumor, it becomes less likely that the preoperative eGFR is caused by CKD risk factors. This is not taken into account by the model, which assumes a low preoperative eGFR has the same likelihood of causing low postoperative eGFR, regardless of the underlying reason (because the CKD risk factors are largely unmeasured, and not accounted for in the model). Consequently, patients with low preoperative eGFR caused by something that is unlikely to be associated with ongoing functional deterioration (a large tumor) are compared with patients who have a low preoperative eGFR caused by something that is likely to be associated with ongoing functional deterioration (CKD risk factors). This results in larger tumors being inappropriately seen as protective.

Aside from analytical issues, it is unclear if a direct path between tumor size and CKD upstaging is biologically plausible. Theoretically, tumor size could affect CKD risk through hemodynamic/structural adaptation in the contralateral kidney following nephron reduction in the kidney affected by the expanding tumor.27 It is likely that tumor size does not continue to affect eGFR at 12 postoperative months, as the tumor has been removed and the compensatory process is complete. This is supported by data showing that at 12 months following radical nephrectomy there were no differences in either eGFR or functional volume of the remaining kidney between patients treated for small compared with large tumors .28

Our structural models showed that, when a direct path between tumor size and postoperative eGFR was considered, larger tumors were associated with a higher postoperative eGFR; but, when only a direct path between tumor size and postoperative eGFR was considered, larger tumors were associated with a lower postoperative eGFR (Figure S1B). This complete reversal in the total effects between these two models tends to support the hypothesis that the pathway between tumor size and postoperative eGFR is non-causal.29 Notwithstanding, this analysis is only exploratory in nature, and limited by the observational design of our study.

The strengths of our study lie in its large size and population-based sampling strategy. It is limited by missing data on tumor complexity and albuminuria, likely underestimation of comorbidities, use of a single follow-up eGFR value, and a reasonably short follow-up duration. There were also some missing data for both pre- and postoperative eGFR, particularly from the state of Victoria, which led to a number of the patients being excluded from this study. These exclusions could potentially have led to further selection bias, if data were not missing at random. The presumed reason for these data being missing was that the clinical record was not accessible to investigators (due to the fact that there is a larger number of private pathology providers in Victoria, compared with Queensland, which investigators were not able to access), and that missingness was not related to patient characteristics or management, which makes it less likely that the conclusions of this manuscript were significantly affected. Notwithstanding, this presumption is difficult to test, and should be considered as a limitation in the interpretation of our results.

Our analyses correct the erroneous assumption that larger-sized kidney tumors reduce the likelihood of having CKD after nephrectomy; previous findings to the contrary appear to have been an artifact of the analytical approach. Patient characteristics, with age being an important indicator, determine the likelihood of postoperative CKD. In practical terms, it is important to consider patient characteristics rather than tumor size when assessing the risk of postoperative CKD for patients managed with radical nephrectomy.

Acknowledgments

We acknowledge the work of Leah Laurenson and Christine Hill for managing and coordinating data collection, in Victoria and Queensland, respectively. We thank the IMPROVE investigators for their role in ascertainment of Victorian data. We would like to acknowledge the assistance of the hospitals and health services in both states that facilitated data access. Data collection was funded through a translational research grant from the Victorian Cancer Agency (EOI09_E36) and a project grant from Cancer Council Queensland. RJE was supported by an Australian Government Research Training Stipend. IDD, REN, and SJJ were supported by fellowships from the Australian National Health and Medical Research Council (APP1102604, APP1043029, and APP1061341, respectively). The IMPROVE Investigators are: Graham Giles, Victoria White, Michael Coory, Jeremy Millar, Damien Bolton, Ian Davis, Ingrid Winship, Ian Roos, Melissa Southey, Gianlucca Severi, Miles Prince, Simon Harrison, John Seymour, Michael Jefford, David Hill, and Anne Kay. Data from living kidney donors were supplied by the Australia and New Zealand Dialysis and Transplant Registry (ANZDATA). The interpretation and reporting of these data are the responsibility of the authors and in no way should be seen as an official policy or interpretation of ANZDATA.

Author contributions

RJE and SJJ were responsible for the study design and data analysis. VMW, DMB, MDC, IDD, GGG, REN, STW, and SJJ were involved in conception and design, and acquisition of data for Queensland and Victorian kidney cancer patients. RJE, RSF, and GCG were involved in the conception and design, and acquisition of data for living kidney donors. All authors contributed significantly to the interpretation of results, drafting or critically revising the manuscript, gave final approval of the final version, and agree to be accountable for all aspects of the work.

Disclosure

The authors report no conflicts of interest related to this work.

References

1. Ta AD, Bolton DM, Dimech MK, et al. Contemporary management of renal cell carcinoma (RCC) in Victoria: implications for longer term outcomes and costs. BJU Int. 2013;112(Suppl 2):36–43. doi:10.1111/bju.12204

2. Klarenbach S, Moore RB, Chapman DW, Dong J, Braam B. Adverse renal outcomes in subjects undergoing nephrectomy for renal tumors: a population-based analysis. Eur Urol. 2011;59(3):333–339. doi:10.1016/j.eururo.2010.11.013

3. Salvatore SP, Cha EK, Rosoff JS, Seshan SV. Nonneoplastic renal cortical scarring at tumor nephrectomy predicts decline in kidney function. Arch Pathol Lab Med. 2013;137(4):531–540. doi:10.5858/arpa.2012-0070-OA

4. Shirasaki Y, Tsushima T, Saika T, Nasu Y, Kumon H. Kidney function after nephrectomy for renal cell carcinoma. Urology. 2004;64(1):43–47. doi:10.1016/j.urology.2004.02.039

5. Jeon HG, Choo SH, Sung HH, et al. Small tumour size is associated with new-onset chronic kidney disease after radical nephrectomy in patients with renal cell carcinoma. Eur J Cancer. 2014;50(1):64–69. doi:10.1016/j.ejca.2013.08.018

6. Ohno Y, Nakashima J, Ohori M, et al. Impact of tumor size on renal function and prediction of renal insufficiency after radical nephrectomy in patients with renal cell carcinoma. J Urol. 2011;186(4):1242–1246. doi:10.1016/j.juro.2011.05.087

7. Dey S, Hamilton Z, Noyes SL, et al. Chronic kidney disease is more common in locally advanced renal cell carcinoma. Urology. 2017;105:101–107. doi:10.1016/j.urology.2017.03.033

8. Levey AS, Coresh J. Chronic kidney disease. Lancet. 2012;379(9811):165–180. doi:10.1016/S0140-6736(11)60178-5

9. Anderson RG, Bueschen AJ, Lloyd LK, Dubovsky EV, Burns JR. Short-term and long-term changes in renal function after donor nephrectomy. J Urol. 1991;145(1):11–13.

10. Ahn T, Ellis RJ, White VM, et al. Predictors of new-onset chronic kidney disease in patients managed surgically for T1a renal cell carcinoma: an Australian population-based analysis. J Surg Oncol. 2018;117(7):1597–1610. doi:10.1002/jso.25037

11. Amin MB, Edge SB, Greene F, et al. AJCC Cancer Staging Manual. 8thed. New York: Springer; 2017.

12. White SL, Polkinghorne KR, Atkins RC, Chadban SJ. Comparison of the prevalence and mortality risk of CKD in Australia using the CKD epidemiology collaboration (CKD-EPI) and Modification of Diet in Renal Disease (MDRD) study GFR estimating equations: the AusDiab (Australian diabetes, obesity and lifestyle) study. Am J Kidney Dis. 2010;55(4):660–670. doi:10.1053/j.ajkd.2009.12.011

13. Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–612.

14. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–383.

15. Levidiotis V. Live kidney donors - assessment and follow up. Aust Fam Physician. 2009;38(5):316–320.

16. Kidney Disease: Improving Global Outcomes CKD Workgroup. KDIGO 2012 clinical practice guideline for the evaluation and management of chronic kidney disease. Kidney Int Suppl. 2013;3(1):1–150.

17. Simpson EH. The interpretation of interaction in contingency tables. J Royal Statistics Soc B. 1951;13(2):238–241.

18. Baker SG, Kramer BS. Good for women, good for men, bad for people: Simpson‘s paradox and the importance of sex-specific analysis in observational studies. J Womens Health Gend Based Med. 2001;10(9):867–872. doi:10.1089/152460901753285769

19. Denic A, Glassock RJ, Rule AD. Structural and functional changes with the aging kidney. Adv Chronic Kidney Dis. 2016;23(1):19–28. doi:10.1053/j.ackd.2015.08.004

20. Aschinberg LC, Koskimies O, Bernstein J, Nash M, Edelmann CM, Spitzer A. The influence of age on the response to renal parenchymal loss. Yale J Biol Med. 1978;51(3):341–345.

21. Davies HT, Crombie IK, Tavakoli M. When can odds ratios mislead? BMJ. 1998;316(7136):989–991.

22. Berkson J. Limitations of the application of fourfold table analysis to hospital data. Biometrics. 1946;2(3):47–53.

23. Glesby MJ, Hoover DR. Survivor treatment selection bias in observational studies: examples from the AIDS literature. Ann Intern Med. 1996;124(11):999–1005.

24. Cole SR, Platt RW, Schisterman EF, et al. Illustrating bias due to conditioning on a collider. Int J Epidemiol. 2010;39(2):417–420. doi:10.1093/ije/dyp334

25. Pearl J. Lord‘s paradox revisited (Oh Lord! Kumbaya!). J Causal Inference. 2016;4(2):20160021. doi:10.1515/jci-2016-0021

26. VanderWeele TJ. Mediation analysis: a practitioner‘s guide. Annu Rev Public Health. 2016;37:17–32. doi:10.1146/annurev-publhealth-032315-021402

27. Peters G. Introduction: history and problems of compensatory adaptation of renal functions and of compensatory hypertrophy of the kidney. Yale J Biol Med. 1978;51(3):235–245.

28. Park BH, Sung HH, Jeong BC, et al. Tumor size is associated with compensatory hypertrophy in the contralateral kidney after radical nephrectomy in patients with renal cell carcinoma. Int Urol Nephrol. 2016;48(6):977–983. doi:10.1007/s11255-016-1250-y

29. Pearl J. Understanding Simpson’s paradox. Am Stat. 2014;68(1):8–13. doi:10.1080/00031305.2014.876829

Supplementary materials

Table S1 Clinical characteristics of nephrectomy patients grouped by indication

Table S2 Application of exclusion criteria between patients grouped by age and tumor size in 1,981 patients diagnosed with RCCa

Table S3 Association between tumor size and postoperative kidney function, stratified by age

Table S4 Comparisons of postoperative kidney function patients grouped by age and tumor size with narrowed exclusion criteria

Table S5 Direct, indirect and total effects of various exposures on pre- and postoperative eGFR

Figure S1 Structural models with and without a direct causal pathway between tumor size and postoperative estimated glomerular filtration rate (eGFR). (A) Structural equation model evaluating the relative contributions of the most significant patient and tumor characteristics on postoperative eGFR; this model assumes there is a direct causal effect of tumor size on postoperative eGFR. The net total effect of tumor size on postoperative eGFR is positive (β: 0.5, 95% CI: 0.4, 0.5). (B) The same model as A, except there is no direct causal association between tumor size and postoperative eGFR. The net total effect of tumor size on postoperative eGFR reverses direction (β: −0.4, 95% CI: −0.6, −0.2).

Creative Commons License © 2019 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.