Introduction

In 1973 the International Agency for Research on Cancer (IARC) declared asbestos in all its commercial forms as established human carcinogen [1].

However, worldwide asbestos exposure has not yet stopped. The Global Burden of Disease (GBD) project estimates that 125 million people are exposed to asbestos globally each year [2], with chrysotile accounting for more than 95% of all the asbestos used globally [3].

Asbestos-induced carcinogenesis is a complex event resulting from different factors including physicochemical characteristics of the fibers (dimension, surface reactivity and chemical composition), time and dose of exposure and, finally, host genetic determinants [4].

Occupational and non-occupational asbestos exposures play a well-known role in determining of asbestosis, pulmonary function decline, lung cancer and mesothelioma [5,6,7]; in addition, several studies investigated relationship between workplace asbestos exposure and others type of cancer.

Although asbestos damages primarily the airways (oral cavity, pharynx, larynx, and lung), some studies suggested a correlation between asbestos occupational exposure and kidney cancer [3].

Kidney cancer is the 14th most common cancer worldwide and ninth most common cancer in Europe. The incidence varies globally, with the highest rates in more developed regions such as North America and Europe and the lowest rates in Asia and Africa [8].

Apart from genetic factors and family history of kidney cancer, established risk factors such as smoking, obesity and hypertension are still prevalent. In addition to such lifestyle factors, environmental agents (such as solvents [9, 10], including in particular trichloroethylene (TCE) [11], pesticides, dusts, diesel and polycyclic aromatic hydrocarbons (PAH) [12]. and medications are thought to be important determinants as well. Physical inactivity, excessive alcohol consumption, unhealthy body weight and poor dietary habits could account for more than 20% of cancer cases [13].

Association between asbestos exposure and kidney cancer remains controversial; to contribute to clarify this issue we conduct an update of previous systematic review and meta-analysis in order to search more recent data about the role of occupational exposure to asbestos fibers in determining kidney cancer.

Methods

We conducted a systematic review and meta-analysis of cohort studies of workers employed in industries entailing asbestos exposure, according to the PRISMA guidelines [14]. A PRISMA checklist is included as Appendix 1. We retained the studies included in the meta-analysis by Sali and Boffetta [15], and performed a systematic search of literature published from 2001 to May 2020. Three different databases were searched: MEDLINE (PubMed), Scopus and Ovid (Embase). We used a string included all types of cancer and exposure to asbestos (without limitation of types of employment or cancer); the string is reported in Appendix 2.

The study protocol was uploaded on the PROSPERO database and is available from the authors.

Selection of studies

Searching through databases, 1443 potentially relevant articles were identified. First, two authors (PB and CZ) independently screened papers the titles and abstracts of the articles and selected those matching the inclusion and exclusion criteria. After that, the same two authors read full texts of papers that were considered potentially eligible and selected those that reported relevant results. In case of disagreements, consensus was reached.

We also searched for additional articles through the lists of recent reviews and metanalyses and those of the articles retained for the meta-analysis.

Inclusion criteria of studies for the meta-analysis were (1) Design: cohort and nested case–control study; (2) Asbestos exposure: workers with exclusive or predominant occupational exposure to asbestos (3) Outcome: incidence or mortality due to kidney cancer (C64-ICD10).

Exclusion criteria consisted in (1) Design: community-based studies, cross-sectional and case–control studies (2) type of exposure to asbestos: we excluded studies of workers who might have occasional, low-level exposure to asbestos, as well as exposure to other known carcinogens, such as seamen, bitumen workers, and mechanics.

In Fig. 1 is shown the flowchart for selection of the studies in meta-analysis. Details on the studies retained in each step of the process are available from the authors.

Fig. 1
figure 1

Flowchart for the identification of articles for the meta-analyses

Extraction of data

For each mortality study we extracted the standardized mortality ratio (SMR) of kidney cancer and its 95% confidence interval (CI); when these measures were not directly available from publications, we calculated them. Similarly, we abstracted or calculated (if not specified in the text) the standardized incidence ratio (SIR) corresponding CI for cancer incidence studies. Incidence studies obtained data from cancer registries, while mortality studies were based on information from death certificates.

In addition, we extracted data about (1) type of industry or occupation of employment (for example shipyard, asbestos textile or asbestos mining); (2) asbestos type, classified as pure chrysotile (Ch), predominantly chrysotile (PCh), amphiboles (Am), and mixed or unspecified (Mix); (3) geographic region, categorized as Europe, North America, and others; (4) source of diagnosis: cancer registry, death certificate; (5) composition of the cohort by sex: female when cohort was formed up to 80% by women, male when there were up to 80% of men, and mixed in case of both male and women comprised less than 80% of the entire cohort; (6) years of employment; (7) proportion of cohort members lost to follow-up; (8) period of first exposure: we distinguished cohorts first exposed to asbestos before 1945, from 1946 to 1960 and after 1961; (9) duration of exposure and time since first exposure (latency).

We also performed quality assessment (QA) according to NIH score for quality assessment tool for Observational Cohort and Cross-Sectional Studies. We modified predefined criteria basing on our object of study [16]; we retained thirteen criteria and classified studies into three categories based on the score obtained: (1) high quality (20.25–15 points); (2) medium quality (15–10 points); (3) low quality (0–10 points).

Statistical analysis

We conducted random-effects [17]. meta-analyses of the RR both combining results for cancer incidence and mortality (when both were available for the same cohort, we used the SIR), and separating them. We also conducted stratified meta-analyses according to the characteristics described above, to explore potential sources of heterogeneity, and a cumulative meta-analysis, according to year of publication of subsequent studies.

To evaluated results stability, we performed sensitivity analyses by quality score and repeated the meta-analysis after excluding one study at a time. Furthermore, we considered the funnel plot and performed the Egger’s regression asymmetry test to assess publication bias [18].

Analyses were performed by STATA16 program [19], using specific commands metan, metabias, and metafunnel.

Results

In total, 927 publications were identified from search conducted among three databases from 2001 to 2020. We found out 243 results from PubMed (MEDLINE), 491 from Scopus and 193 from Embase, of which 414 were removed because they were duplicates.

We screened the titles and abstracts of 513 articles and rejected 382 of them because not relevant (Fig. 1), and retained 131 articles as potentially eligible. After reading full-texts of 122 of them (full texts were not available for 9), we excluded 109 articles that did not meet inclusion criteria.

We included in the review 13 studies retrieved from the search since 2001, in addition to 29 non-overlapping studies out of 33 included in previous meta-analysis [15]. Out of the 42 studies, we abstracted 54 results from 51 separate cohorts (Table 1).

Table 1 Selected characteristics of cohorts included in the meta-analysis

Among the 51 cohorts, 34 comprised predominantly men, six predominantly women, and 11 both men and women. Furthermore, 34 cohorts were from Europe, nine from North America and eight from other countries. The main type of asbestos fiber was chrysotile in 21 cohorts and amphiboles in seven, while in the remaining 23 cohorts it was not possible to distinguish between the two types. Results on kidney cancer mortality were available for 43 cohorts, and those on cancer incidence for 11 cohorts (for three cohorts both types of results were available).

Figure 2 presents the pooled results based on random-effected model. The pooled RR was 1.14 (95% CI 1.04–1.29); eight results were excluded from the meta-analysis because they were based on zero observed events There was limited evidence of heterogeneity between studies (p = 0.3; I2 = 11%). The summary RR of results on cancer incidence was 0.98 (95% CI 0.79–1.22) with no evidence of heterogeneity between studies (p = 0.8, I2 = 0%); that of results on mortality 1.15 (95% CI 0.99–1.34; p = 0.2; I2 = 19%;).

Fig. 2
figure 2

Forest plot of meta-analysis of results on risk of kidney cancer among occupational cohorts exposed to asbestos

Table 2 shows the results of stratified analyses, that revealed no difference by asbestos type, region, sex, or period of first exposure.

Table 2 Results of meta-analyses stratified according to selected characteristics

We performed quality assessment by NIH score for studies included in our analysis published after 2000. Only two studies were considered good quality because of good exposure assessment (several measurements), low proportion of lost to follow up and analysis adjusted for others risk factors of kidney cancer. Eight studies were of moderate quality, only few of them included a satisfactory description of sample size justification or prior measurement of the outcome.

Some authors considered also others circumstances of asbestos exposure. Kjaerheim et al. (2005) studied the effect of exposure to asbestos in drinking water; the SIR for kidney cancer was 1.4 (95% CI 0.6–2.7). Van den Borre et al. (2015) distinguished between manual workers and non-manual workers exposed to asbestos, and found no differences between the two groups, with SMR equal to 1.88 (95% CI 0.23–6.80) and 3.24 (95% CI 0.08–18.06), respectively.

Only two studies reported data about latency or duration of exposure [33, 34].

As shown in Fig. 3, we found no evidence of publication bias in the main results of the cohorts included in the meta-analysis (p-value of Egger’s test, 0.4 for mortality studies and 0.9 for incidence studies).

Fig. 3
figure 3

Funnel plot of Egger’s test to assess publication bias among studies included in the meta-analysis

Discussion

Kidney cancer is the 14th common cancer in the world and, while 50% of kidney cancer are not attributable to known risk factors for that malignancy [8], only few studies considered its association with asbestos exposure. Furthermore, is difficult to observed cohorts of workers exposed only to asbestos fibers, as many of them, entail exposure to other potential kidney cancerogenic risk factors.

Our results show a lack of association between occupational asbestos exposure and kidney cancer; indeed, pooled analysis shows relative risk 1.1 (95% CI 1.0–1.3). This result was confirmed in subgroup analyses stratified by type of asbestos fibers, country, sex or period of first exposure. These results are also consistent with those of a previous meta-analysis of studies published up to 2000 [15].

Our study suffers from several limitations. The most important problem is that seven studies we selected after systematic review of literature showed no cases of kidney cancers among cohorts exposed to asbestos, and were excluded from the meta-analysis, that was based on log-transformed measures of associations. In order to assess the magnitude of the bias away from the null, we performed also a meta-analysis without logarithm transformation both including and excluding studies with no observed cases or deaths, and the results were not materially changed by the exclusion of these latter studies.

Our analysis included a small number of studies based on cancer incidence, and results of studies on cancer mortality might be affected by determinants of survival, such as access to effective treatment. Furthermore, many studies evidenced low accuracy of COD, especially among elderly people who more frequently have a multiple comorbidity before death (62).

The incidence of kidney cancer is often underestimated as it is a tumor that is diagnosed late, sometimes when metastases are already present in other sites and it is not always assumed that it is possible to trace the primary tumor (especially as regards data relating to tumors that have arisen several years ago, when the analysis of the tumor histotype was not so frequent and therefore it was even more difficult to trace the starting site of the tumor).

Furthermore, only three studies [51, 54, 63] considered tobacco smoking as potential confounding factor, and none included obesity, alcohol consumption, diabetes, hypertension, analgesic therapy in the analysis [64].

In addition only two papers (mentioned before) reported data about latency and duration of exposure, for that reason we decided to not perform a dose–response meta-analysis.

Finally, exposure assessment remains an issue in occupation cohort studies based on job titles or industries of employment, without direct exposure measurements. Such approach does not distinguish between circumstances at different levels of exposure, leading to possible underestimate of the effects.

We decided to review the data about risk for kidney cancer among workers exposed to asbestos due to well-known role of asbestos fibers in determining carcinogenic process in other organs and some supporting biological data. Boor et al. showed an increase of renal fibrosis and production of inflammatory mediators locally in rats exposed to amiosite fibers [65]. and that kind of data support the fact that asbestos fibers could be transported by blood to renal cells and there could trigger inflammatory process and, subsequently, genetic alteration.

More evidence of capacity for asbestos fibers to reach kidney has been provided by necropsy studies of humans; moreover, asbestos fibers have been founded in human urine [66].

A strength of our study is the large number of cohort studies included in the meta-analysis. Furthermore, the limited heterogeneity between study results, the consistency of results across different subgroups of studies, and the lack of evidence of publication bias support the main conclusion of our meta-analysis.

In conclusion our updated meta-analysis provided additional evidence to the hypothesis of a lack of association between occupational exposure to asbestos and risk of kidney cancer.