Long-term glucocorticoid (GC) exposure causes immunosuppression; therefore, the risk of cancer may be increased in long-term GC users. We investigated whether long-term GC use is associated with a higher risk of cancer in the population without cancer. A population-based cohort study using data from the National Health Insurance Service was conducted among the South Korean adult population in 2010. Long-term GC users were defined as those who were prescribed a continuous supply of oral GC for ≥30 days. The primary endpoint was a new cancer diagnosis from January 1, 2011, to December 31, 2015. Among 770,880 individuals included in the analysis, 1,602 (0.2%) were long-term GC users and 36,157 (4.7%) were newly diagnosed with cancer from January 2011 to December 2015. In the multivariable Cox regression analysis, the risk of cancer among long-term GC users was 1.23-fold higher than that of the unexposed individuals [95% confidence interval (CI), 1.06–1.43; P = 0.007]. In the competing risk analyses, the risks of liver cancer and lung cancer were 1.46-fold (95% CI, 1.03–2.07; P = 0.034) and 1.52-fold (95% CI, 1.04–2.21; P = 0.029) higher in the long-term GC users than that of the unexposed individuals, respectively. We found that long-term GC exposure might be associated with a higher risk of overall cancer, and this association was more evident for lung and liver cancer risk. However, because there might be unmeasured and potential confounders in this study, the results should be interpreted carefully, and future studies should be performed to confirm these findings.

Impact:

Long-term glucocorticoid therapy might be associated with a higher cancer risk. This association was more evident for lung and liver cancer risk. Our findings suggest that long-term prescriptions of glucocorticoids should be administered carefully considering the risk of cancer.

Cancer is one of the major causes of death worldwide (1). From 2006 to 2016, there were 17.2 million patients with cancer worldwide, 8.9 million deaths, and the global incidence of cancer increased by 28% (2). In South Korea, the cancer prevalence has increased very fast, and 229,180 Koreans were newly diagnosed and 78,194 died from cancer in 2016 (3). Furthermore, the incidence of cancer is expected to increase; therefore, cancer prevention is emphasized to achieve a significant reduction in the global burden of disease (4).

Glucocorticoids (GC) are a class of pleiotropic steroid hormones that are commonly prescribed for 1.2% and 1.0% of patients with chronic conditions in the United States and the United Kingdom, respectively (5, 6). GCs have both immunosuppressive and potent anti-inflammatory effects (7), and they inhibit the immune response and production of prostaglandins and leukotrienes, which are the two main products of inflammation (8). Based on its immunosuppressive and anti-inflammatory effects, GCs have been prescribed to patients with chronic obstructive pulmonary disease (COPD) or asthma (9), autoimmune disease (10), rheumatic disease (11), and allergic diseases (12). GCs have also been used in the management of patients with cancer to relieve symptoms such as pain, dyspnea, neurological manifestations, and discomfort from inflammatory reactions and lymphedema (13). GCs have also been widely used in clinical oncology during chemotherapy and radiotherapy to alleviate side effects in patients with nonhematological cancer (14). However, long-term GC therapy is known to cause serious adverse effects such as immunosuppression, and it could impair immunotherapy outcomes (15, 16). Therefore, it is possible that long-term GC use might increase the incidence of cancer; however, there is not enough evidence regarding this issue.

This study aimed to investigate whether long-term GC use is associated with a higher risk of cancer in the population without cancer. In addition, we assessed whether the cancer risk associated with oral GC use is dependent on the daily dose.

Data source and ethical statement

The sample cohort database of the National Health Insurance Service (NHIS) was developed to provide data for research, including health surveys and medical surveys of the Korean population. The validated database comprised a stratified random sample of approximately one million people who were registered with the NHIS in 2002, and it was designed to be representative of the national population in terms of demographic and socioeconomic variables. The cohort was dynamic and followed-up until the end of 2015. It was supplemented with additional cohort data, including data regarding infants, to allow for attrition due to death and loss to follow-up. Each year, using stratified extraction methods to ensure that the cohort was representative of the national population, people were added to the cohort to replace those who had died or emigrated during the previous year. The sample cohort database included information regarding individuals' demographics, socioeconomic status, healthcare use, medical history, and cause of death until the end of 2015 (17). The study protocol was approved by the Institutional Review Board of Seoul National University Bundang Hospital (X-1903-531-901) and the Health Insurance Review and Assessment Service (NHIS-2019-2-140).

Study population

We assessed a sample cohort of 1,000,000 individuals included in the 2010 cohort database of the NHIS. Next, we excluded individuals who died in 2010 and those who had emigrated between 2011 and 2015 because they might have been diagnosed with cancer after emigrating. In addition, individuals with a history of cancer in 2010 were excluded, as our study focused on new cancer diagnoses between 2011 and 2015.

Long-term glucocorticoid use in 2010 as an exposure variable

Individuals prescribed oral GC (prednisolone, methylprednisolone, or dexamethasone) regularly and continuously over ≥30 days were defined as long-term GC users, based on previous report (18). The individuals who were not prescribed any GC or prescribed GC for <30 days were classified as unexposed individuals. Because the longest prescription period of oral GC was from 3 to 7 days in a year for short-term GC users in South Korea, we considered them as unexposed individuals, too. Three types of oral GCs were considered, prednisolone, methylprednisolone, and dexamethasone, and the daily dosages of GCs were divided into two groups, ≥5 mg/day (high-dosage) and <5 mg/day (low-dosage) of prednisolone. Using this dosage calculation, 4 mg/day of methylprednisolone and 0.75 mg/day of dexamethasone were considered equivalent to 5 mg/day of prednisolone, given the potency of GCs (19). Our rationale for dividing the groups using a cut-off value of 5 mg/day was based on previous studies that have reported no significant association with dosages <5 mg/day of prednisolone and all-cause mortality for patients with rheumatoid arthritis (20, 21). As a result, long-term GC users receiving <5 and ≥5 mg/day of prednisolone were classified as the low-dosage group and high-dosage group, respectively.

Development of cancer as the dependent variable

According to the International Classification of Diseases 10th edition (ICD-10) codes, newly registered diagnoses of any malignancy (C00-C96) from January 1, 2011 to December 31, 2015, were defined as the development of cancer in this study. In details, the cancers were categorized as follows: gastric cancer (C16), esophageal cancer (C15), colorectal cancer (C18C20), gall bladder and biliary tract cancer (C23–C24), head and neck cancer (C00–C14), brain cancer (C71), liver cancer (C22), pancreatic cancer (C25), lung cancer (C34), bone and articular cartilage cancer (C40–C41), neoplasms of the breast and genital organs (C50–C63), urinary tract cancer (C64–C68), thyroid cancer (C73), and lymphoma or leukemia (C81–C96). The time to cancer diagnosis was calculated starting from January 1, 2011, to the date of diagnosis of cancer, as registered officially in the ICD-10 system. In South Korea, all patients diagnosed with any C-code cancer should be registered in the NHIS database to receive financial coverage of 95% of the total charges for cancer treatment from the NHIS. Therefore, all patients diagnosed with cancer were registered in NHIS database.

Confounding variables

The following data were considered confounders in this study: demographic information (age and sex); socioeconomic information [income level and place of residence in 2010 (Seoul, metropolitan cities, and other)]; the 2010 Charlson comorbidity index, which was calculated using registered ICD-10 diagnostic codes in the NHIS database between 2009 and 2010 (Supplementary Table S1); other comorbidities [hypertension (I10), COPD (J44), asthma (J45), allergic disease (T78.4), and autoimmune disease (Supplementary Table S2)]; and total number of hospital and outpatient visit days in 2010. The entire cohort was divided into four groups according to age (18–30, 31–50, 51–70, and ≥71 years old). The cohort in 2010 was also divided into five groups according to the income level (0%–20%, 20%–40%, 40%–60%, 60%–80%, and 80%–100%). Regarding the classification of place of residence, Seoul, the capital city, was assigned a separate category, and the cities of Incheon, Kwangju, Busan, Ulsan, Daegu, and Daejeon were classified as metropolitan cities. The number of hospital days included the number of hospital outpatient clinic visits and days spent admitted to the hospital. For example, if an individual visited a hospital outpatient clinic five times and was admitted to the hospital for 3 days, then the individual was considered to have 8 hospital visit days. In the analysis, the number of hospital and outpatient clinic visit days was categorized into five groups (0 days, 1–7 days, 8–29 days, 30–90 days, and >90 days) and four groups (0–7 days, 8–29 days, 30–90 days, and >90 days), respectively.

Study endpoint

The primary endpoint of this study was the new development of cancer from January 1, 2011, to December 31, 2015, among patients registered in the 2010 NHIS sample cohort of South Korea. All individuals without any cancer diagnosis were followed up until December 31, 2015, if they did not die. Individuals who died during follow-up were censored as noncases at the date of death.

Statistical analysis

We performed multivariable Cox regression analyses to determine whether long-term GC exposure in 2010 was associated with the development of cancer from 2011 to 2015, and compared the results with those of unexposed individuals. All covariates except the Charlson comorbidity index were included in the multivariable model for adjustment to avoid multicollinearity with other comorbidities that were used to calculate the Charlson comorbidity index (Supplementary Table S1). In addition, subgroup analysis according to the daily dosage of GC was performed to investigate whether the daily dosage of GC might affect the main results of this study. In this analysis, the individuals were classified according to the daily dosage of GCs, as follows: the low daily dosage group (<5 mg/day of prednisolone) and high daily dosage group (≥5 mg/day of prednisolone). Next, in the competing risk analyses, we constructed in detail 14 multivariable Cox regression models for the development of cancer during 2011–2015. The 14 types of cancer were each set as endpoints in the multivariable Cox models, and the duration to cancer diagnosis from January 1, 2011, was used for the time-to-event competing risk analyses. These competing risk analyses were performed to investigate what type of cancer risk was most closely related to long-term GC exposure. In addition, subgroup analyses according to sex were performed to analyze the hazard of overall cancer and all 14 cancer types. We also fitted a multivariable Cox regression model for the development of cancer during 2012–2015 (not for 2011–2015) in the entire cohort to assess whether reverse causation bias could have affected the result because there was a short latency time between long-term GC exposure and development of cancer in 2011. Finally, subgroup analyses according to Charlson comorbidity index (0–2 and >2), COPD, asthma, autoimmune disease, rheumatic disease other than autoimmune disease, allergic disease, and liver disease in 2010 were performed to assess whether the comorbid status of individuals or indication bias of long-term GC exposure affected the results of this study. The results of the Cox regression models are presented as hazard ratios (HR) with 95% confidence intervals (CI), and it was confirmed that there was no multicollinearity in all multivariable models of the entire cohort (variance inflation factor of <2.0). It was also confirmed that the central assumption of the Cox proportional assumption was met for all multivariable models using log-minus-log plots.

To enhance the robustness of this study, we performed propensity score (PS) matching, which is an efficient method of reducing the effects of confounders in observational studies (22), using the nearest neighbor method with a 1:5 ratio, without replacement, and a caliper width of 0.2. All covariates were included in the PS model, and a logistic regression analysis was performed to calculate the PS as a logistic model. The absolute value of the standardized mean difference (ASD) was used to evaluate the balance between long-term GC users and the control group before and after PS matching. Through PS matching, we aimed to minimize the ASD of all confounders between the two groups to less than 0.1. Since some important information such as body mass index (BMI), alcohol consumption and smoking status in 2010 were available in only 20% of all individuals (17), they were not included in the PS modeling. They were compared after PS matching among individuals whose data was were available, using Chi-square test. After confirming a good degree of balance between the two groups in the PS-matched cohort, we performed conditional Cox regression via a stratified Cox regression analysis of the development of cancer from 2011 to 2015. All statistical analyses were performed using R version 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria); P < 0.05 was considered statistically significant.

The 2010 sample cohort comprised 1,000,000 individuals. We excluded 173,091 individuals younger than 18 years, 4,487 individuals who died in 2010, 208 individuals who emigrated from 2011 to 2015, and 51,334 individuals who had a history of cancer in 2010. Therefore, 770,880 individuals were included in the analysis (Fig. 1). Among them, 1,602 (0.2%) were long-term GC users and 769,278 (99.8%) were in the unexposed individuals. In the 2010 cohort, 36,157 (4.7%) individuals were newly diagnosed with cancer from January 2011 to December 2015. The median duration from January 1, 2011, to the date of the cancer diagnosis was 2.4 years (interquartile range: 1.1–3.7 years).

Figure 1.

Flow chart showing the selection of individuals in the 2010 cohort and reasons for exclusion.

Figure 1.

Flow chart showing the selection of individuals in the 2010 cohort and reasons for exclusion.

Close modal

Cancer development in the entire cohort

Table 1 shows the results of the multivariable Cox regression analyses for cancer development during 2011 to 2015 for the entire sample cohort in 2010 (n = 770,880). In the multivariable model, the risk for cancer development among long-term GC users was 1.23-fold higher than that of the unexposed individuals (HR = 1.23; 95 % CI, 1.06–1.43; P = 0.004; model 1). In addition, the risk of cancer development among high daily dosage long-term GC users was 1.29-fold higher (HR = 1.29; 95% CI, 1.03–1.62; P = 0.030; model 2) than that of the unexposed individuals, whereas that of the low daily dosage long-term GC users was not significantly different (P = 0.096). Supplementary Table S3 shows the results of multivariable Cox regression analyses of the development of cancer during 2012 to 2015 (not for 2011–2015). In the multivariable model, the risk for cancer development among long-term GC users was 1.20-fold higher than that of the unexposed individuals (HR = 1.20; 95% CI, 1.01–1.43; P = 0.042; model 1). In addition, the risk of cancer development among high daily dosage long-term GC users was 1.28-fold higher (HR = 1.28; 95% CI, 1.02–1.61; P = 0.031; model 2) than that of the unexposed individuals, whereas that of the low daily dosage long-term GC users was not significantly different (P = 0.075).

Table 1.

Multivariable Cox regression analysis of development of cancer during 2011 to 2015 in the 2010 sample cohort of South Korea (n = 770,880).

Multivariable model
VariableHazard ratio (95% CI)P value
Age, year 
 18–30  
 31–50 2.53 (2.42–2.65) <0.001 
 51–70 5.34 (5.10–5.59) <0.001 
 ≥ 71 7.82 (7.43–8.24) <0.001 
Sex: Male (vs. female) 1.47 (1.44–1.50) <0.001 
Income level 
 0%–20% (Lowest income level)  
 20%–40% 1.01 (0.97–1.05) 0.643 
 40%–60% 1.04 (1.00–1.07) 0.062 
 60%–80% 0.99 (0.95–1.02) 0.471 
 80%–10% (Highest income level) 1.06 (1.02–1.10) 0.001 
Residence 
 Capital city (Seoul)  
 Metropolitan city 1.02 (0.99–1.05) 0.259 
 Other area 1.00 (0.97–1.02) 0.737 
Comorbidities 
 Hypertension 0.99 (0.96–1.02) 0.366 
 DM with chronic complication 1.14 (1.09–1.18) <0.001 
 DM without chronic complication 1.16 (1.13–1.20) <0.001 
 Cerebrovascular disease 1.06 (1.02–1.10) 0.002 
 Myocardial infarction 1.12 (1.03–1.20) 0.005 
 Congestive heart failure 0.99 (0.93–1.04) 0.579 
 Chronic obstructive pulmonary disease 1.32 (1.24–1.40) <0.001 
 Asthma 1.09 (1.06–1.12) <0.001 
 Other chronic pulmonary disease 1.11 (1.07–1.14) <0.001 
 Allergic disease 0.96 (0.87–1.06) 0.437 
 Autoimmune disease 1.07 (1.02–1.13) 0.010 
 Dementia 1.03 (0.90–1.18) 0.710 
 Hemi- or paraplegia 0.97 (0.88–1.07) 0.543 
 Renal disease 1.17 (1.08–1.28) <0.001 
 Mild liver disease 1.30 (1.27–1.34) <0.001 
 Severe liver disease 1.53 (1.40–1.67) <0.001 
 Peptic ulcer disease 1.12 (1.09–1.14) <0.001 
 Rheumatic disease (not autoimmune disease) 1.05 (1.00–1.11) 0.059 
 Peripheral vascular disease 1.00 (0.97–1.04) 0.821 
Hospital visit in 2010, day 
 0  
 1–7 1.08 (1.04–1.11) <0.001 
 8–30 1.16 (1.11–1.21) <0.001 
 30–90 1.19 (1.10–1.29) <0.001 
 >90 1.01 (0.91–1.13) 0.839 
Outpatient clinic visit in 2010, day 
 0–7  
 8–30 1.27 (1.24–1.31) <0.001 
 30–90 1.28 (1.23–1.34) <0.001 
 >90 1.04 (0.81–1.34) 0.756 
Long-term GC users (model 1) 
 Unexposed individuals  
 Long-term GC users 1.23 (1.06–1.43) 0.007 
Long-term GC users (model 2) 
 Unexposed individuals  
 Low daily dosage of chronic GC user 1.18 (0.97–1.44) 0.096 
 High daily dosage of chronic GC user 1.29 (1.03–1.62) 0.030 
A model adjusted only for age (model 3) 
 Unexposed individuals  
 Long-term GC users 1.52 (1.31–1.76) <0.001 
Multivariable model
VariableHazard ratio (95% CI)P value
Age, year 
 18–30  
 31–50 2.53 (2.42–2.65) <0.001 
 51–70 5.34 (5.10–5.59) <0.001 
 ≥ 71 7.82 (7.43–8.24) <0.001 
Sex: Male (vs. female) 1.47 (1.44–1.50) <0.001 
Income level 
 0%–20% (Lowest income level)  
 20%–40% 1.01 (0.97–1.05) 0.643 
 40%–60% 1.04 (1.00–1.07) 0.062 
 60%–80% 0.99 (0.95–1.02) 0.471 
 80%–10% (Highest income level) 1.06 (1.02–1.10) 0.001 
Residence 
 Capital city (Seoul)  
 Metropolitan city 1.02 (0.99–1.05) 0.259 
 Other area 1.00 (0.97–1.02) 0.737 
Comorbidities 
 Hypertension 0.99 (0.96–1.02) 0.366 
 DM with chronic complication 1.14 (1.09–1.18) <0.001 
 DM without chronic complication 1.16 (1.13–1.20) <0.001 
 Cerebrovascular disease 1.06 (1.02–1.10) 0.002 
 Myocardial infarction 1.12 (1.03–1.20) 0.005 
 Congestive heart failure 0.99 (0.93–1.04) 0.579 
 Chronic obstructive pulmonary disease 1.32 (1.24–1.40) <0.001 
 Asthma 1.09 (1.06–1.12) <0.001 
 Other chronic pulmonary disease 1.11 (1.07–1.14) <0.001 
 Allergic disease 0.96 (0.87–1.06) 0.437 
 Autoimmune disease 1.07 (1.02–1.13) 0.010 
 Dementia 1.03 (0.90–1.18) 0.710 
 Hemi- or paraplegia 0.97 (0.88–1.07) 0.543 
 Renal disease 1.17 (1.08–1.28) <0.001 
 Mild liver disease 1.30 (1.27–1.34) <0.001 
 Severe liver disease 1.53 (1.40–1.67) <0.001 
 Peptic ulcer disease 1.12 (1.09–1.14) <0.001 
 Rheumatic disease (not autoimmune disease) 1.05 (1.00–1.11) 0.059 
 Peripheral vascular disease 1.00 (0.97–1.04) 0.821 
Hospital visit in 2010, day 
 0  
 1–7 1.08 (1.04–1.11) <0.001 
 8–30 1.16 (1.11–1.21) <0.001 
 30–90 1.19 (1.10–1.29) <0.001 
 >90 1.01 (0.91–1.13) 0.839 
Outpatient clinic visit in 2010, day 
 0–7  
 8–30 1.27 (1.24–1.31) <0.001 
 30–90 1.28 (1.23–1.34) <0.001 
 >90 1.04 (0.81–1.34) 0.756 
Long-term GC users (model 1) 
 Unexposed individuals  
 Long-term GC users 1.23 (1.06–1.43) 0.007 
Long-term GC users (model 2) 
 Unexposed individuals  
 Low daily dosage of chronic GC user 1.18 (0.97–1.44) 0.096 
 High daily dosage of chronic GC user 1.29 (1.03–1.62) 0.030 
A model adjusted only for age (model 3) 
 Unexposed individuals  
 Long-term GC users 1.52 (1.31–1.76) <0.001 

Abbreviation: DM, diabetes mellitus.

Competing risk analyses

Table 2 shows the results of the competing risk analyses according to the primary cancer sites, which were obtained from the multivariable Cox regression analyses for the development of cancer during 2011–2015. First, the risk of overall cancer in male long-term GC users was 1.28-fold higher (HR: 1.28, 95% CI: 1.06–1.25; P = 0.001) than that in the unexposed individuals, while that was not significantly different (P = 0.144) in female GC users. Additionally, the risks of liver cancer and lung cancer were 1.46-fold (HR:1.46; 95% CI: 1.03–2.07; P = 0.034) and 1.52-fold (HR:1.52; 95% CI: 1.04–2.21; P = 0.029) higher in the long-term GC users than that in the unexposed individuals, respectively.

Table 2.

Multivariable Cox regression analysis for development of cancer according to primary site during 2011 to 2015.

Multivariable model
Type of cancer (long-term GC users vs. unexposed individuals; reference)HR (95% CI)P value
Overall cancer (n = 36,157, long-term GC user: 176) 1.23 (1.06–1.43) 0.007 
 Male 1.28 (1.06–1.25) 0.001 
 Female 1.18 (0.95–1.47) 0.144 
Gastric cancer (C16) (n = 2,909–long-term GC user: 11) 1.04 (0.57–1.88) 0.903 
 Male 1.10 (0.49–2.46) 0.822 
 Female 1.00 (0.41–2.41) 0.994 
Esophageal cancer (C15) (n = 199, long-term GC user: 2) 2.44 (0.60–10.00) 0.215 
 Male 3.10 (0.75–12.77) 0.118 
 Female 0.00 (0.00-) 0.969 
Colorectal cancer (C18–C20) (n = 4,181, long-term GC user: 18) 1.04 (0.66–1.66) 0.859 
 Male 1.14 (0.59–2.20) 0.707 
 Female 0.98 (0.51–1.89) 0.947 
GB and biliary tract cancer (C23–C24) (n = 671, long-term GC user: 1) 0.34 (0.05–2.41) 0.278 
 Male 0.00 (0.00-) 0.912 
 Female 0.57 (0.08–4.05) 0.569 
Head and Neck cancer (C00–C14) (n = 218, long-term GC user: 2) 2.17 (0.53–8.89) 0.283 
 Male 0.00 (0.00-) 0.943 
 Female 5.99 (1.40–20.71) 0.016 
Brain cancer (C71) (n = 189, long-term GC user: 1) 1.41 (0.19–10.23) 0.737 
 Male 5.04 (0.68–37.66) 0.115 
 Female 0.00 (0.00-) 0.975 
Liver cancer (C22) (n = 6,320, long-term GC user: 32) 1.46 (1.03–2.07) 0.034 
 Male 1.61 (1.00–2.61) 0.050 
 Female 1.32 (0.79–2.20) 0.288 
Pancreatic cancer (C25) (n = 1,845, long-term GC user: 12) 1.51 (0.85–2.67) 0.161 
 Male 1.63 (0.67–3.96) 0.280 
 Female 1.46 (0.69–3.08) 0.327 
Lung cancer (C34) (n = 3,136, long-term GC user: 28) 1.52 (1.04–2.21) 0.029 
 Male 1.40 (0.82–2.37) 0.219 
 Female 1.67 (0.98–2.85) 0.058 
Bone, articular cartilage cancer (C40–C41) (n = 87, long-term GC user: 0) 0.00 (0.00-) 0.985 
 Male 0.00 (0.00-) 0.994 
 Female 0.00 (0.00-) 0.985 
Neoplasms of breast and genital organs (C50–C63) (n = 9,584, long-term GC user: 40) 1.17 (0.85–1.60) 0.333 
 Male 1.32 (0.93–1.89) 0.124 
 Female 0.84 (0.44–1.63) 0.612 
Urinary tract cancer (C64–C68) (n = 1,542, long-term GC user: 6) 0.90 (0.40–2.02) 0.798 
 Male 0.96 (0.31–3.02) 0.949 
 Female 0.83 (0.27–2.61) 0.756 
Thyroid cancer (C73) (n = 2,258, long-term GC user: 11) 1.62 (0.89–2.94) 0.117 
 Male 2.02 (0.50–8.23) 0.326 
 Female 1.55 (0.80–3.01) 0.196 
Lymphoma or Leukemia (C81–C96) (n = 960, long-term GC user: 1) 0.22 (0.03–1.59) 0.134 
 Male 0.00 (0.00-) 0.896 
 Female 0.37 (0.05–2.65) 0.322 
Other site cancer (n = 2,058, long-term GC user: 11) 1.21 (0.67–2.20) 0.528 
 Male 0.98 (0.36–2.63) 0.965 
 Female 1.40 (0.66–2.96) 0.381 
Multivariable model
Type of cancer (long-term GC users vs. unexposed individuals; reference)HR (95% CI)P value
Overall cancer (n = 36,157, long-term GC user: 176) 1.23 (1.06–1.43) 0.007 
 Male 1.28 (1.06–1.25) 0.001 
 Female 1.18 (0.95–1.47) 0.144 
Gastric cancer (C16) (n = 2,909–long-term GC user: 11) 1.04 (0.57–1.88) 0.903 
 Male 1.10 (0.49–2.46) 0.822 
 Female 1.00 (0.41–2.41) 0.994 
Esophageal cancer (C15) (n = 199, long-term GC user: 2) 2.44 (0.60–10.00) 0.215 
 Male 3.10 (0.75–12.77) 0.118 
 Female 0.00 (0.00-) 0.969 
Colorectal cancer (C18–C20) (n = 4,181, long-term GC user: 18) 1.04 (0.66–1.66) 0.859 
 Male 1.14 (0.59–2.20) 0.707 
 Female 0.98 (0.51–1.89) 0.947 
GB and biliary tract cancer (C23–C24) (n = 671, long-term GC user: 1) 0.34 (0.05–2.41) 0.278 
 Male 0.00 (0.00-) 0.912 
 Female 0.57 (0.08–4.05) 0.569 
Head and Neck cancer (C00–C14) (n = 218, long-term GC user: 2) 2.17 (0.53–8.89) 0.283 
 Male 0.00 (0.00-) 0.943 
 Female 5.99 (1.40–20.71) 0.016 
Brain cancer (C71) (n = 189, long-term GC user: 1) 1.41 (0.19–10.23) 0.737 
 Male 5.04 (0.68–37.66) 0.115 
 Female 0.00 (0.00-) 0.975 
Liver cancer (C22) (n = 6,320, long-term GC user: 32) 1.46 (1.03–2.07) 0.034 
 Male 1.61 (1.00–2.61) 0.050 
 Female 1.32 (0.79–2.20) 0.288 
Pancreatic cancer (C25) (n = 1,845, long-term GC user: 12) 1.51 (0.85–2.67) 0.161 
 Male 1.63 (0.67–3.96) 0.280 
 Female 1.46 (0.69–3.08) 0.327 
Lung cancer (C34) (n = 3,136, long-term GC user: 28) 1.52 (1.04–2.21) 0.029 
 Male 1.40 (0.82–2.37) 0.219 
 Female 1.67 (0.98–2.85) 0.058 
Bone, articular cartilage cancer (C40–C41) (n = 87, long-term GC user: 0) 0.00 (0.00-) 0.985 
 Male 0.00 (0.00-) 0.994 
 Female 0.00 (0.00-) 0.985 
Neoplasms of breast and genital organs (C50–C63) (n = 9,584, long-term GC user: 40) 1.17 (0.85–1.60) 0.333 
 Male 1.32 (0.93–1.89) 0.124 
 Female 0.84 (0.44–1.63) 0.612 
Urinary tract cancer (C64–C68) (n = 1,542, long-term GC user: 6) 0.90 (0.40–2.02) 0.798 
 Male 0.96 (0.31–3.02) 0.949 
 Female 0.83 (0.27–2.61) 0.756 
Thyroid cancer (C73) (n = 2,258, long-term GC user: 11) 1.62 (0.89–2.94) 0.117 
 Male 2.02 (0.50–8.23) 0.326 
 Female 1.55 (0.80–3.01) 0.196 
Lymphoma or Leukemia (C81–C96) (n = 960, long-term GC user: 1) 0.22 (0.03–1.59) 0.134 
 Male 0.00 (0.00-) 0.896 
 Female 0.37 (0.05–2.65) 0.322 
Other site cancer (n = 2,058, long-term GC user: 11) 1.21 (0.67–2.20) 0.528 
 Male 0.98 (0.36–2.63) 0.965 
 Female 1.40 (0.66–2.96) 0.381 

Subgroup analyses

Table 3 shows the results of multivariable Cox regression analyses of the development of cancer in the subgroups. In individuals whose Charlson comorbidity index was 0 to 2 in 2010, the risk of cancer was 1.34-fold (HR = 1.34; 95% CI, 1.05–1.71; P = 0.020) higher among long-term GC users than among unexposed individuals. In individuals whose Charlson comorbidity index was >2 in 2010, the risk of cancer was 1.23-fold (HR = 1.23; 95% CI, 1.02–1.48; P = 0.031) higher in long-term GC users than in unexposed individuals. In the other subgroup analyses according to comorbidities related to long-term GC exposure, the risk of cancer among long-term GC users who were patients with asthma was 1.31-fold higher than that of unexposed individuals who were asthma patients (HR = 1.31; 95% CI, 1.03–1.66; P = 0.027). All other subgroups of long-term GC users (COPD, autoimmune disease, rheumatic disease other than autoimmune disease, allergic disease, and liver disease) did not show any significant associated between the comorbidity and the risk of overall cancer, compared with unexposed individuals (all P > 0.05). Supplementary Table S4 shows the risk of liver cancer among the long-term GC users with liver disease, and the risk of lung cancer among the long-term GC users with COPD, and asthma. Among the patients with asthma, the risk of lung cancer in the long-term GC users was 2.13-fold higher (HR = 2.13; 95% CI, 1.32–3.46; P = 0.002) than that in the unexposed individuals. However, the risks of lung cancer and liver cancer in the long-term GC users among the patients with COPD (P = 0.091) and those with liver disease (P = 0.177) were not significantly higher than that in unexposed individuals, respectively.

Table 3.

Subgroup analyses.

Multivariable model
VariableHazard ratio (95% CI)P value
Charlson comorbidity index 0–2 (n = 672,608) 
 Unexposed individuals  
 Long-term GC users 1.34 (1.05–1.71) 0.020 
Charlson comorbidity index >2 (n = 98,272) 
 Unexposed individuals  
 Long-term GC users 1.23 (1.02–1.48) 0.031 
Chronic obstructive pulmonary disease (n = 8,399) 
 Unexposed individuals  
 Long-term GC users 1.12 (0.77–1.62) 0.556 
Asthma (n = 87,577) 
 Unexposed individuals  
 Long-term GC users 1.31 (1.03–1.66) 0.027 
Autoimmune disease (n = 17,746) 
 Unexposed individuals  
 Long-term GC users 1.37 (0.94–1.98) 0.099 
Rheumatic disease other than autoimmune disease (n = 19,745) 
 Unexposed individuals  
 Long-term GC users 1.32 (1.00–1.73) 0.051 
Allergic disease (n = 6,659) 
 Unexposed individuals  
 Long-term GC users 0.31 (0.04–2.24) 0.246 
Liver disease (n = 120,334) 
 Unexposed individuals  
 Long-term GC users 1.26 (1.00–1.59) 0.052 
Multivariable model
VariableHazard ratio (95% CI)P value
Charlson comorbidity index 0–2 (n = 672,608) 
 Unexposed individuals  
 Long-term GC users 1.34 (1.05–1.71) 0.020 
Charlson comorbidity index >2 (n = 98,272) 
 Unexposed individuals  
 Long-term GC users 1.23 (1.02–1.48) 0.031 
Chronic obstructive pulmonary disease (n = 8,399) 
 Unexposed individuals  
 Long-term GC users 1.12 (0.77–1.62) 0.556 
Asthma (n = 87,577) 
 Unexposed individuals  
 Long-term GC users 1.31 (1.03–1.66) 0.027 
Autoimmune disease (n = 17,746) 
 Unexposed individuals  
 Long-term GC users 1.37 (0.94–1.98) 0.099 
Rheumatic disease other than autoimmune disease (n = 19,745) 
 Unexposed individuals  
 Long-term GC users 1.32 (1.00–1.73) 0.051 
Allergic disease (n = 6,659) 
 Unexposed individuals  
 Long-term GC users 0.31 (0.04–2.24) 0.246 
Liver disease (n = 120,334) 
 Unexposed individuals  
 Long-term GC users 1.26 (1.00–1.59) 0.052 

Development of cancer after propensity score adjustment

After PS matching, a total of 9,592 individuals (1,599 long-term GC users and 7,991 controls) were included in the final analysis. The demographic and clinical characteristics of the two groups are shown in Table 4, and all ASDs between the two groups after PS matching were <0.1, demonstrating a good degree of balance. However, as shown in Supplementary Table S5, the distribution BMI and alcohol consumption in long-term GC users and unexposed individuals were statistically different in the PS-matched cohort (P < 0.05), while smoking status in the two groups was not statistically different (P = 0.626). Table 5 shows the results of cancer newly diagnosed from 2011 to 2015 in the PS-matched cohort of 2010. In the stratified Cox regression analysis, the risk for cancer development among long-term GC users was 1.25-fold higher than that of the unexposed individuals (HR = 1.25; 95% CI, 1.06–1.48; P = 0.008).

Table 4.

Results of comparison of characteristics between long-term GC users and unexposed individuals before and after PS matching.

Entire cohort (n = 770,880)PS-matched cohort (n = 9,592)
Long-term GC usersUnexposed individualsASDLong-term GC usersUnexposed individualsASD
Variablesn = 1,602n = 769,278n = 1,599n = 7,991
Age, year 
 18–30 109 (6.8) 185,282 (24.1)  109 (6.8) 625 (7.8)  
 31–50 438 (27.3) 335,178 (43.6) 0.36 438 (27.4) 2,216 (27.7) <0.01 
 51–70 660 (41.2) 192,971 (25.1) 0.33 658 (41.2) 3,296 (41.2) <0.01 
 ≥ 71 395 (24.7) 55,847 (7.3) 0.40 394 (24.6) 1,854 (23.2) 0.03 
 Sex: male 609 (38.0) 383,518 (49.9) 0.24 609 (38.1) 3,085 (38.6) 0.01 
Income level 
 0%–20% (Lowest income level) 215 (13.4) 107,196 (13.9)  214 (13.4) 1,080 (13.5)  
 20%–40% 209 (13.0) 126,083 (16.4) 0.10 209 (13.1) 1,093 (13.7) 0.02 
 40%–60% 447 (27.9) 172,054 (22.4) 0.12 445 (27.8) 2,143 (26.8) 0.02 
 60%–80% 315 (19.7) 166,538 (21.6) 0.05 315 (19.7) 1,597 (20.0) <0.01 
 80%–10% (Highest income level) 416 (26.0) 197,407 (25.7) <0.01 416 (26.0) 2,078 (26.0) <0.01 
Residence 
 Capital city (Seoul) 273 (17.0) 163,416 (21.2)  273 (17.1) 1,414 (17.7)  
 Metropolitan city 359 (22.4) 196,593 (25.6) 0.08 357 (22.3) 1,872 (23.4) 0.03 
 Other area 970 (60.5) 409,269 (53.2) 0.15 969 (60.6) 4,705 (58.9) 0.03 
 Charlson comorbidity index 2.9 (2.3) 1.0 (1.5) 0.84 2.9 (2.2) 2.8 (2.3) 0.07 
 Hypertension 712 (44.4) 129,971 (16.9) 0.55 709 (44.3) 3,436 (43.0) 0.03 
 DM with chronic complication 180 (11.2) 28,781 (3.7) 0.24 178 (11.1) 847 (10.6) 0.02 
 DM without chronic complication 428 (26.7) 72,437 (9.4) 0.39 425 (26.6) 2,027 (25.4) 0.03 
 Cerebrovascular disease 238 (14.9) 38,170 (5.0) 0.28 236 (14.8) 1,183 (14.8) <0.01 
 Myocardial infarction 55 (3.4) 5,839 (0.8) 0.15 54 (3.4) 251 (3.1) 0.01 
 Congestive heart failure 142 (8.9) 14,012 (1.8) 0.25 139 (8.7) 666 (8.3) <0.01 
 Chronic obstructive pulmonary disease 161 (10.0) 8,238 (1.1) 0.30 158 (9.9) 755 (9.4) 0.01 
 Asthma 466 (29.1) 87,111 (11.3) 0.39 463 (29.0) 2,157 (27.0) 0.04 
 Other chronic pulmonary disease 289 (18.0) 91,099 (11.8) 0.16 289 (18.1) 1,395 (17.5) 0.02 
 Allergic disease 26 (1.6) 6,633 (0.9) 0.06 26 (1.6) 141 (1.8) 0.01 
 Autoimmune disease 220 (13.7) 17,526 (2.3) 0.33 220 (13.8) 1,333 (16.7) 0.08 
 Dementia 11 (0.7) 1,637 (0.2) 0.06 11 (0.7) 52 (0.7) <0.01 
 Hemi- or paraplegia 32 (2.0) 4,564 (0.6) 0.10 32 (2.0) 152 (1.9) <0.01 
 Renal disease 76 (4.7) 4,720 (0.6) 0.19 74 (4.6) 294 (3.7) 0.04 
 Mild liver disease 508 (31.7) 117,060 (15.2) 0.35 505 (31.6) 2,436 (30.5) 0.02 
 Severe liver disease 22 (1.4) 4,738 (0.6) 0.07 21 (1.3) 111 (1.4) <0.01 
 Peptic ulcer disease 791 (49.4) 171,815 (22.3) 0.54 788 (49.3) 3,745 (46.9) 0.05 
 Rheumatic disease other than autoimmune disease 519 (32.4) 19,226 (2.5) 0.64 516 (32.3) 2,332 (29.2) 0.06 
 Peripheral vascular disease 369 (23.0) 56,626 (7.4) 0.37 367 (23.0) 1,751 (21.9) 0.02 
Hospital visit in 2010, day 
 0 1,091 (68.1) 673,644 (87.6)  1,091 (68.2) 5,669 (70.9)  
 1–7 219 (13.7) 60,734 (7.9) 0.17 219 (13.7) 1,061 (13.3) 0.01 
 8–30 206 (12.9) 25,595 (3.3) 0.28 205 (12.8) 894 (11.2) 0.05 
 30–90 71 (4.4) 5,611 (0.7) 0.18 69 (4.3) 293 (3.7) 0.03 
 >90 15 (0.9) 3,694 (0.5) 0.05 15 (0.9) 74 (0.9) <0.01 
Outpatient clinic visit in 2010, day 
 0–7 258 (16.1) 487,842 (63.4)  258 (16.1) 1,516 (19.0)  
 8–30 1,056 (65.9) 251,705 (32.7) 0.70 1,056 (66.0) 5,165 (64.6) 0.03 
 30–90 280 (17.5) 29,176 (3.8) 0.36 279 (17.4) 1,275 (16.0) 0.04 
 >90 8 (0.5) 555 (0.1) 0.06 6 (0.4) 35 (0.4) <0.01 
Entire cohort (n = 770,880)PS-matched cohort (n = 9,592)
Long-term GC usersUnexposed individualsASDLong-term GC usersUnexposed individualsASD
Variablesn = 1,602n = 769,278n = 1,599n = 7,991
Age, year 
 18–30 109 (6.8) 185,282 (24.1)  109 (6.8) 625 (7.8)  
 31–50 438 (27.3) 335,178 (43.6) 0.36 438 (27.4) 2,216 (27.7) <0.01 
 51–70 660 (41.2) 192,971 (25.1) 0.33 658 (41.2) 3,296 (41.2) <0.01 
 ≥ 71 395 (24.7) 55,847 (7.3) 0.40 394 (24.6) 1,854 (23.2) 0.03 
 Sex: male 609 (38.0) 383,518 (49.9) 0.24 609 (38.1) 3,085 (38.6) 0.01 
Income level 
 0%–20% (Lowest income level) 215 (13.4) 107,196 (13.9)  214 (13.4) 1,080 (13.5)  
 20%–40% 209 (13.0) 126,083 (16.4) 0.10 209 (13.1) 1,093 (13.7) 0.02 
 40%–60% 447 (27.9) 172,054 (22.4) 0.12 445 (27.8) 2,143 (26.8) 0.02 
 60%–80% 315 (19.7) 166,538 (21.6) 0.05 315 (19.7) 1,597 (20.0) <0.01 
 80%–10% (Highest income level) 416 (26.0) 197,407 (25.7) <0.01 416 (26.0) 2,078 (26.0) <0.01 
Residence 
 Capital city (Seoul) 273 (17.0) 163,416 (21.2)  273 (17.1) 1,414 (17.7)  
 Metropolitan city 359 (22.4) 196,593 (25.6) 0.08 357 (22.3) 1,872 (23.4) 0.03 
 Other area 970 (60.5) 409,269 (53.2) 0.15 969 (60.6) 4,705 (58.9) 0.03 
 Charlson comorbidity index 2.9 (2.3) 1.0 (1.5) 0.84 2.9 (2.2) 2.8 (2.3) 0.07 
 Hypertension 712 (44.4) 129,971 (16.9) 0.55 709 (44.3) 3,436 (43.0) 0.03 
 DM with chronic complication 180 (11.2) 28,781 (3.7) 0.24 178 (11.1) 847 (10.6) 0.02 
 DM without chronic complication 428 (26.7) 72,437 (9.4) 0.39 425 (26.6) 2,027 (25.4) 0.03 
 Cerebrovascular disease 238 (14.9) 38,170 (5.0) 0.28 236 (14.8) 1,183 (14.8) <0.01 
 Myocardial infarction 55 (3.4) 5,839 (0.8) 0.15 54 (3.4) 251 (3.1) 0.01 
 Congestive heart failure 142 (8.9) 14,012 (1.8) 0.25 139 (8.7) 666 (8.3) <0.01 
 Chronic obstructive pulmonary disease 161 (10.0) 8,238 (1.1) 0.30 158 (9.9) 755 (9.4) 0.01 
 Asthma 466 (29.1) 87,111 (11.3) 0.39 463 (29.0) 2,157 (27.0) 0.04 
 Other chronic pulmonary disease 289 (18.0) 91,099 (11.8) 0.16 289 (18.1) 1,395 (17.5) 0.02 
 Allergic disease 26 (1.6) 6,633 (0.9) 0.06 26 (1.6) 141 (1.8) 0.01 
 Autoimmune disease 220 (13.7) 17,526 (2.3) 0.33 220 (13.8) 1,333 (16.7) 0.08 
 Dementia 11 (0.7) 1,637 (0.2) 0.06 11 (0.7) 52 (0.7) <0.01 
 Hemi- or paraplegia 32 (2.0) 4,564 (0.6) 0.10 32 (2.0) 152 (1.9) <0.01 
 Renal disease 76 (4.7) 4,720 (0.6) 0.19 74 (4.6) 294 (3.7) 0.04 
 Mild liver disease 508 (31.7) 117,060 (15.2) 0.35 505 (31.6) 2,436 (30.5) 0.02 
 Severe liver disease 22 (1.4) 4,738 (0.6) 0.07 21 (1.3) 111 (1.4) <0.01 
 Peptic ulcer disease 791 (49.4) 171,815 (22.3) 0.54 788 (49.3) 3,745 (46.9) 0.05 
 Rheumatic disease other than autoimmune disease 519 (32.4) 19,226 (2.5) 0.64 516 (32.3) 2,332 (29.2) 0.06 
 Peripheral vascular disease 369 (23.0) 56,626 (7.4) 0.37 367 (23.0) 1,751 (21.9) 0.02 
Hospital visit in 2010, day 
 0 1,091 (68.1) 673,644 (87.6)  1,091 (68.2) 5,669 (70.9)  
 1–7 219 (13.7) 60,734 (7.9) 0.17 219 (13.7) 1,061 (13.3) 0.01 
 8–30 206 (12.9) 25,595 (3.3) 0.28 205 (12.8) 894 (11.2) 0.05 
 30–90 71 (4.4) 5,611 (0.7) 0.18 69 (4.3) 293 (3.7) 0.03 
 >90 15 (0.9) 3,694 (0.5) 0.05 15 (0.9) 74 (0.9) <0.01 
Outpatient clinic visit in 2010, day 
 0–7 258 (16.1) 487,842 (63.4)  258 (16.1) 1,516 (19.0)  
 8–30 1,056 (65.9) 251,705 (32.7) 0.70 1,056 (66.0) 5,165 (64.6) 0.03 
 30–90 280 (17.5) 29,176 (3.8) 0.36 279 (17.4) 1,275 (16.0) 0.04 
 >90 8 (0.5) 555 (0.1) 0.06 6 (0.4) 35 (0.4) <0.01 

Note: Data are presented as number (percentage) or mean (SD).

Abbreviation: DM, diabetes mellitus.

Table 5.

Development of cancer from 2011 to 2015 before and after PS adjustment.

VariablesDevelopment of cancer (%)HR (95% CI)P value
Before PS adjustment 
Unexposed individuals 35,981 of 769,278 (4.7)  
Long-term GC users 176 of 1,602 (11.0) 2.44 (2.11–2.83) <0.001 
After PS adjustment 
Unexposed individuals 762 of 7,991 (9.5)  
Long-term GC users 175 of 1,599 (10.9) 1.25 (1.06–1.48) 0.008 
VariablesDevelopment of cancer (%)HR (95% CI)P value
Before PS adjustment 
Unexposed individuals 35,981 of 769,278 (4.7)  
Long-term GC users 176 of 1,602 (11.0) 2.44 (2.11–2.83) <0.001 
After PS adjustment 
Unexposed individuals 762 of 7,991 (9.5)  
Long-term GC users 175 of 1,599 (10.9) 1.25 (1.06–1.48) 0.008 

This population-based cohort study showed that long-term GC use was associated with a higher incidence of cancer compared with that of the unexposed individuals. This association was more evident for lung and liver cancer risk. These findings suggested that long-term GC therapy might be a risk factor for cancer, and that long-term GC therapy should be performed carefully because of the cancer risk.

Regarding humoral immunity, long-term GC therapy was reported to be related to a reversible decrease in B-cell counts and antibody responses, even in low dosages (23). Furthermore, long-term GC therapy is known to affect innate immunity by inhibiting macrophage differentiation and suppressing the production of IL1, IL6, TNF, and the pro-inflammatory prostaglandins and leukotrienes by macrophages (24). Because both suppressed innate immunity with enhanced humoral immunity were known as related factors with increased overall cancer risk (25), the influence of immunosuppression may be present in long-term GC users, which might have affected the results of this study. The long-term GC use also inhibits key inflammatory markers such as reactive protein/TNFα, which is known to play a role in carcinogenesis and surpass the antitumor immune responses (26). Moreover, long-term GC therapy is known to suppress the tumoricidal activities of activated macrophages (7); therefore, the risk of cancer might be increased in long-term GC users.

The results of the subgroup analyses of the daily dosage of GC showed that the HR for cancer was higher and significant among the high daily dosage group of long-term GC users, which suggested that the risk of cancer with long-term GC use might be dosage-dependent. In general, the adverse effects of long-term GC therapy are known to be dosage-dependent (27), and our previous retrospective cohort study reported that high daily dosage GC use was significantly associated with a higher risk for 30-day mortality for critically ill patients (28). However, our results should be interpreted carefully, because the 95% CIs in HRs of low daily dosage (0.99–1.50) and high daily dosage GC users (1.14–1.84) were overlapped, suggesting there was no significant difference in HRs for overall cancer development between low daily dosage and high daily dosage GC users. Therefore, the effect of daily dosage in long-term GC users on cancer risk should be evaluated in the future study.

In competing risk analyses, the associations between long-term GC exposure and elevated cancer risk were significant for lung and liver cancer in this study. Recent population-based cohort studies reported that long-term exposure to inhaled corticosteroids in patients with COPD was associated with a lower risk of lung cancer (29, 30). There were differences between these previous studies and ours because we focused on long-term oral GC therapy, whereas previous studies focused on inhaled GC therapy (29, 30). COPD is a well-known risk factor for lung cancer, and inhaled GCs comprise a useful treatment option for preventing inflammatory reactions, which is an important pathophysiologic mechanism underlying the development of lung cancer (31). Furthermore, inhaled GCs are known to have fewer severe adverse effects than orally administered GC (32). In a nested case–control study in Taiwan, the higher risk of squamous lung cancer was significantly associated with recent dose increase of GC among patients with asthma or COPD (33). The study included both the inhaled GC use and oral GC use, and it suggested that long-term oral GC exposure might have a potential risk factor for development of lung cancer. However, the relationship between long-term oral GC exposure and lung cancer risk in population without cancer or patients with COPD should be evaluated in the future with more studies

For liver cancer, there was a very recent study that reported that long-term exposure of oral GC was associated with higher risk of hepatocellular carcinoma in patients who had any one of comorbidities including alcohol-related disease, chronic liver disease, and diabetes mellitus (34). Although the previous study included on the patients with higher risk of liver cancer, our study focused on the individuals without cancer, suggesting that the association of long-term GC exposure might be a potential risk factor for development of liver cancer in healthy individuals without cancer. More study is needed to confirm these associations of long-term GC exposure with liver and lung cancer risk in future.

This population-based cohort study had several limitations. First, some important data, such as BMI, alcoholic consumption, and smoking status, were not included in the analysis for adjustment because they were not present in the entire cohort dataset. In fact, there were statistical differences in distribution of BMI and alcohol consumption in PS-matched cohort as shown in Supplementary Table S4, and it might have effect on the result in this study. Second, we defined the underlying diseases using the ICD-10 codes registered in the NHIS database. The diseases specified by the ICD-10 codes may have differed from the actual underlying diseases. Third, because we assessed the population without cancer in sample cohort, which consisted of many healthy adults, the proportion of long-term GC users was relatively small (0.2%). Fourth, there might have been residual confounders in this study that were not included in the PS model or multivariable adjustment, which could have affected the results of this study. Finally, our findings may have been influenced by the indication bias of GC. For example, although we included some comorbidities such as rheumatic disease and COPD for both PS modeling and multivariable adjustment, the severity of the diseases may differ, leading to differences in GC use. This may have affected the incidence of cancer in long-term GC users.

In conclusion, this population-based cohort study conducted in South Korea showed that long-term GC exposure might be associated with a higher risk of overall cancer. This association was more evident for lung and liver cancer risk. However, there might be unmeasured and potential confounders in this study, therefore, the results should be interpreted carefully, and future studies should be performed to confirm these findings.

No potential conflicts of interest were disclosed.

T.K. Oh: Conceptualization, data curation, formal analysis and writing-original draft. I.-A. Song: Conceptualization, supervision, methodology, writing-review and editing.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Bray
F
,
Ferlay
J
,
Soerjomataram
I
,
Siegel
RL
,
Torre
LA
,
Jemal
A
. 
Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries
.
CA Cancer J Clin
2018
;
68
:
394
424
.
2.
Global Burden of Disease Cancer Collaboration
,
Fitzmaurice
C
,
Akinyemiju
TF
,
Al Lami
FH
,
Alam
T
,
Alizadeh-Navaei
R
, et al
Global, regional, and National Cancer Incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2016: a systematic analysis for the global burden of disease study
.
JAMA Oncol
2018
;
4
:
1553
68
.
3.
Jung
KW
,
Won
YJ
,
Kong
HJ
,
Lee
ES
. 
Cancer statistics in korea: incidence, mortality, survival, and prevalence in 2016
.
Cancer Res Treat
2019
;
51
:
417
30
.
4.
Global Burden of Disease Cancer Collaboration
,
Fitzmaurice
C
,
Allen
C
,
Barber
RM
,
Barregard
L
,
Bhutta
ZA
, et al
Global, regional, and National Cancer Incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015: a systematic analysis for the global burden of disease study
.
JAMA Oncol
2017
;
3
:
524
48
.
5.
Overman
RA
,
Yeh
JY
,
Deal
CL
. 
Prevalence of oral glucocorticoid usage in the United States: a general population perspective
.
Arthritis Care Res
2013
;
65
:
294
8
.
6.
van Staa
TP
,
Leufkens
HG
,
Abenhaim
L
,
Begaud
B
,
Zhang
B
,
Cooper
C
. 
Use of oral corticosteroids in the United Kingdom
.
QJM
2000
;
93
:
105
11
.
7.
Boumpas
DT
,
Chrousos
GP
,
Wilder
RL
,
Cupps
TR
,
Balow
JE
. 
Glucocorticoid therapy for immune-mediated diseases: basic and clinical correlates
.
Ann Intern Med
1993
;
119
:
1198
208
.
8.
Goppelt-Struebe
M
,
Wolter
D
,
Resch
K
. 
Glucocorticoids inhibit prostaglandin synthesis not only at the level of phospholipase A2 but also at the level of cyclo-oxygenase/PGE isomerase
.
Br J Pharmacol
1989
;
98
:
1287
95
.
9.
Keatings
VM
,
Jatakanon
A
,
Worsdell
YM
,
Barnes
PJ
. 
Effects of inhaled and oral glucocorticoids on inflammatory indices in asthma and COPD
.
Am J Respir Crit Care Med
1997
;
155
:
542
8
.
10.
Tait
AS
,
Butts
CL
,
Sternberg
EM
. 
The role of glucocorticoids and progestins in inflammatory, autoimmune, and infectious disease
.
J Leukoc Biol
2008
;
84
:
924
31
.
11.
Berthelot
JM
. 
Comments about the article by Mouterde et al. entitled "Indications of glucocorticoids in early arthritis and rheumatoid arthritis: recommendations for clinical practice based on data from the literature and expert opinions"
.
Joint Bone Spine
2010
;
77
:
597
603
. 
Low-dose prednisone and biologics: allies rather than competitors?
Joint Bone Spine
2012
;
79
:
103
4
;
author reply 104–6
.
12.
Dimeloe
S
,
Nanzer
A
,
Ryanna
K
,
Hawrylowicz
C
. 
Regulatory T cells, inflammation and the allergic response-The role of glucocorticoids and Vitamin D
.
J Steroid Biochem Mol Biol
2010
;
120
:
86
95
.
13.
Herr
I
,
Pfitzenmaier
J
. 
Glucocorticoid use in prostate cancer and other solid tumours: implications for effectiveness of cytotoxic treatment and metastases
.
Lancet Oncol
2006
;
7
:
425
30
.
14.
Gennari
A
,
Salvadori
B
,
Tognoni
A
,
Conte
PF
. 
Rapid intravenous premedication with dexamethasone prevents hypersensitivity reactions to paclitaxel
.
Ann Oncol
1996
;
7
:
978
9
.
15.
Wolfe
F
,
Caplan
L
,
Michaud
K
. 
Treatment for rheumatoid arthritis and the risk of hospitalization for pneumonia: associations with prednisone, disease-modifying antirheumatic drugs, and anti-tumor necrosis factor therapy
.
Arthritis Rheum
2006
;
54
:
628
34
.
16.
Flint
TR
,
Janowitz
T
,
Connell
CM
,
Roberts
EW
,
Denton
AE
,
Coll
AP
, et al
Tumor-induced IL-6 reprograms host metabolism to suppress anti-tumor immunity
.
Cell Metab
2016
;
24
:
672
84
.
17.
Lee
J
,
Lee
JS
,
Park
SH
,
Shin
SA
,
Kim
K
. 
Cohort profile: the National Health Insurance Service-National Sample Cohort (NHIS-NSC), South Korea
.
Int J Epidemiol
2017
;
46
:
e15
.
18.
Longui
CA
. 
Glucocorticoid therapy: minimizing side effects
.
J Pediatr
2007
;
83
:
S163
77
.
19.
Meikle
AW
,
Tyler
FH
. 
Potency and duration of action of glucocorticoids. Effects of hydrocortisone, prednisone and dexamethasone on human pituitary-adrenal function
.
Am J Med
1977
;
63
:
200
7
.
20.
del Rincon
I
,
Battafarano
DF
,
Restrepo
JF
,
Erikson
JM
,
Escalante
A
. 
Glucocorticoid dose thresholds associated with all-cause and cardiovascular mortality in rheumatoid arthritis
.
Arthritis Rheumatol
2014
;
66
:
264
72
.
21.
Movahedi
M
,
Costello
R
,
Lunt
M
,
Pye
SR
,
Sergeant
JC
,
Dixon
WG
. 
Oral glucocorticoid therapy and all-cause and cause-specific mortality in patients with rheumatoid arthritis: a retrospective cohort study
.
Eur J Epidemiol
2016
;
31
:
1045
55
.
22.
Rosenbaum
PR
,
Rubin
DB
. 
Reducing bias in observational studies using subclassification on the propensity score
.
J Am Statist Assoc
1984
;
79
:
516
24
.
23.
Fedor
ME
,
Rubinstein
A
. 
Effects of long-term low-dose corticosteroid therapy on humoral immunity
.
Ann Allergy Asthma Immunol
2006
;
97
:
113
6
.
24.
McEwen
BS
,
Biron
CA
,
Brunson
KW
,
Bulloch
K
,
Chambers
WH
,
Dhabhar
FS
, et al
The role of adrenocorticoids as modulators of immune function in health and disease: neural, endocrine and immune interactions
.
Brain Res Brain Res Rev
1997
;
23
:
79
133
.
25.
Tan
T-T
,
Coussens
LM
. 
Humoral immunity, inflammation and cancer
.
Curr Opin Immunol
2007
;
19
:
209
16
.
26.
Ham
B
,
Fernandez
MC
,
D'Costa
Z
,
Brodt
P
. 
The diverse roles of the TNF axis in cancer progression and metastasis
.
Trends Cancer Res
2016
;
11
:
1
27
.
27.
Walsh
LJ
,
Wong
CA
,
Oborne
J
,
Cooper
S
,
Lewis
SA
,
Pringle
M
, et al
Adverse effects of oral corticosteroids in relation to dose in patients with lung disease
.
Thorax
2001
;
56
:
279
84
.
28.
Oh
TK
,
Song
IA
,
Lee
JH
,
Lim
C
,
Jeon
YT
,
Bae
HJ
, et al
Effect of preadmission glucocorticoid therapy on 30-day mortality in critically ill patients: a retrospective study of a mixed ICU population in a tertiary hospital
.
Ann Intensive Care
2019
;
9
:
8
.
29.
Raymakers
AJN
,
Sadatsafavi
M
,
Sin
DD
,
FitzGerald
JM
,
Marra
CA
,
Lynd
LD
. 
Inhaled corticosteroids and the risk of lung cancer in COPD: a population-based cohort study
.
Eur Respir J
2019
;
53
:
1801257
.
30.
Liu
SF
,
Kuo
HC
,
Lin
MC
,
Ho
SC
,
Tu
ML
,
Chen
YM
, et al
Inhaled corticosteroids have a protective effect against lung cancer in female patients with chronic obstructive pulmonary disease: a nationwide population-based cohort study
.
Oncotarget
2017
;
8
:
29711
21
.
31.
Durham
AL
,
Adcock
IM
. 
The relationship between COPD and lung cancer
.
Lung Cancer
2015
;
90
:
121
7
.
32.
Donahue
JG
,
Weiss
ST
,
Livingston
JM
,
Goetsch
MA
,
Greineder
DK
,
Platt
R
. 
Inhaled steroids and the risk of hospitalization for asthma
.
JAMA
1997
;
277
:
887
91
.
33.
Jian
ZH
,
Huang
JY
,
Lin
FC
,
Nfor
ON
,
Jhang
KM
,
Ku
WY
, et al
The use of corticosteroids in patients with COPD or asthma does not decrease lung squamous cell carcinoma
.
BMC Pulm Med
2015
;
15
:
154
.
34.
Lai
SW
,
Lin
CL
,
Liao
KF
. 
Synergistic effect of oral corticosteroids use on risk of hepatocellular carcinoma in high risk populations
.
Eur J Intern Med
2018
;
52
:
73
7
.