Meta-analysis results
The literature search and selection process was presented in Fig. 1. The comprehensive search generated 4292 potentially relevant studies, of which 4206 articles were excluded based on title and abstract review. We then examined full texts of the remaining 83 studies and excluded 39 of them. To be specific, 17 studies with irrelevant topics, 16 studies in which the study design was not cohort, four studies with incomplete data to calculate the estimates, and two studies being reviews were excluded. Additionally, four articles were included through previous review or reference lists of relevant studies. Finally, 48 eligible articles were included in the meta-analysis.
The characteristics of the included studies were summarized in Supplementary Table S1. 48 articles with a total of 247,575 patients with SLE (more than 178,332 females and 19,900 males) were represented. Among them, 43 studies estimated the cancer incidence in SLE patients with the follow-up period ranging from 1,000 person-years to 157,969 person-years or from 4.7 years to 35.3 years. Seven studies estimated cancer-specific death in SLE patients with the follow-up period ranging from 48 person-years to 91,669 person-years or from 8.1 years to 11.9 years. In addition, a total of 30 human cancers were systematically divided into six systemic groups (digestive cancers, respiratory cancers, reproductive cancers, urinary cancers, hematopoietic cancers, and other cancers).
The relationships between SLE and cancers were shown in Table 1 and Fig. 2. Our results suggested an increased cancer risk in SLE patients (RR = 1.62, 95% CI, 1.47–1.79). Site-specific analysis suggested that SLE were associated with an increased risk of lymphoma, Hodgkin’s lymphoma, non-Hodgkin Lymphoma, leukemia, multiple myeloma, as well as esophagus, colon, anal, hepatobiliary, liver, pancreatic, larynx, lung, cervical, vagina/vulva, renal, bladder, skin (non-melanoma), and thyroid cancer. However, no significant associations were observed between SLE and cancers of stomach, colorectum, rectal, lip, oral cavity and pharynx, breast, ovary, uterus, prostate, melanoma, and brain.
The results of subgroup analysis by region were presented in Table 1 and Fig. 3. SLE patients were associated with an increased risk of overall cancers in Europe (RR = 1.66, 95% CI, 1.35–2.03), America (RR = 1.58, 95% CI, 1.19–2.09), and Asia (RR = 1.57, 95% CI, 1.39–1.77). In region of Europe, we observed an increased risk of lymphoma, Hodgkin’s lymphoma, non-Hodgkin Lymphoma, leukemia, as well as colon, anal, hepatobiliary, liver, pancreatic, larynx, lung, cervical, vagina/vulva, bladder, and skin (non-melanoma) cancer. However, no significant associations were observed between SLE and cancers of esophagus, stomach, colorectum, rectal, lip, oral cavity and pharynx, breast, ovary, uterus, prostate, kidney, brain and thyroid, as well as Multiple myeloma and melanoma. In region of America, we observed an increased risk of liver cancer, lymphoma, Hodgkin’s lymphoma and non-Hodgkin Lymphoma, whereas a deceased risk of prostate cancer in SLE patients. However, no significant associations were observed between SLE and cancers of colon, lung, breast, ovary, cervix, kidney, and melanoma. In region of Asia, SLE patients were associated with an increased risk of esophagus, hepatobiliary, liver, lip, oral cavity and pharynx, lung, cervical, skin (non-melanoma), brain and thyroid cancer, as well as lymphoma, non-Hodgkin Lymphoma, and leukemia. However, no significant associations were observed between SLE and cancers of stomach, colorectum, pancreatic, breast, ovary, uterus, prostate, kidney, and bladder.
We also estimated the cancer-specific death in SLE patients, and the results were presented in Table 1 and Supplementary Figure S1. The cause-specific standard mortality ratio (SMR) in SLE patients was higher for overall cancer (SMR = 1.52, 95% CI, 1.36–1.70), particularly for lung cancer (SMR = 1.52, 95% CI, 2.13–3.10). However, the cause-specific SMR for hematologic cancer was 1.42 (95% CI, 0.97–2.08).
Mendelian randomization results
The instrumental variable of SLE was constructed using 69 SNPs in European ancestry, explaining approximately 34% of the heritability totally, and the F-statistic was 172.36 (F > 100). In addition, 16 SNPs remained as the final genetic variants to construct instrumental variable of SLE in Asian ancestry, explaining over 1% of the variation, and the F-statistic was 10.47 (F > 10). The details of the SNPs enrolled in our study were represented in Supplementary Table S2 and Supplementary Table S3.
The significant associations (15 for SLE and sites-specific cancers in European ancestry, and 12 for SLE and sites-specific cancers in Asian ancestry) identified by the meta-analysis of cohort studies were further estimated by two-sample Mendelian randomization analysis. We searched for the outcome data using MR-Base and only found the data of colon cancer, pancreatic cancer, lung cancer, cervical cancer, bladder cancer, lymphoma, and Non-melanoma skin cancer in European ancestry, liver cancer and lung cancer in Asian ancestry. The details of outcome included in Mendelian randomization analyses were shown in Supplementary Table S4.
Table 2 presented associations of cancer risk in SLE patients as predicted using the instrument variants in Mendelian randomization analysis. Being consistent with the findings in the meta-analysis, the conventional IVW method indicated a causal association between genetically predisposed SLE and lymphoma (OR = 1.0004, 95% CI, 1.0001–1.0007, P = 0.0035), whereas a decreased risk of bladder cancer (OR = 0.9996, 95% CI, 0.9994–0.9998, P = 0.00004) in European ancestry. Furthermore, the MR-Egger and weighted median methods yielded similar results (Supplementary Table S5). However, no association was observed between genetically predisposed SLE and risk of colon, pancreatic, lung, cervical and Non-melanoma skin cancer in European ancestry, liver cancer and lung cancer in Asian ancestry. Sensitivity analyses were performed to evaluate the robustness of the associations, and the results of horizontal pleiotropy and heterogeneity tests were shown in Supplementary Table S5.
Funnel plots, scatter plots and leave-one-out plots concerning relationship between SLE and lymphoma were shown in Supplementary figures S2-4. In addtion, leave-one-out analysis suggested that no single instrument was strongly driving the overall effect of SLE on other specific cancer types, indicating that these results were not sensitive to SNP selection (data not shown).