Article Text
Abstract
Introduction Diabetes mellitus is known to increase the risk of cancer. Fasting blood glucose (FBG) levels can be changed over time. However, the association between FBG trajectory and cancer risk has been insufficiently studied. This research aims to examine the relationship between FBG trajectories and cancer risk in the Korean population.
Research design and methods We analyzed data from the National Health Insurance Service–National Health Screening Cohort collected between 2002 and 2015. Group-based trajectory modeling was performed on 256,271 Koreans aged 40–79 years who had participated in health examinations at least three times from 2002 to 2007. After excluding patients with cancer history before 2008, we constructed a cancer-free cohort. The Cox proportional hazards model was applied to examine the association between FBG trajectories and cancer incidence by cancer type, after adjustments for covariates. Cancer case was defined as a person who was an outpatient thrice or was hospitalized once or more with a cancer diagnosis code within the first year of the claim.
Results During the follow-up time (2008–2015), 18,991 cancer cases were identified. Four glucose trajectories were found: low-stable (mean of FBG at each wave <100 mg/dL), elevated-stable (113–124 mg/dL), elevated-high (104–166 mg/dL), and high-stable (>177 mg/dL). The high-stable group had a higher risk of multiple myeloma, liver cancer and gastrointestinal cancer than the low-stable group, with HR 4.09 (95% CI 1.40 to 11.95), HR 1.68 (95% CI 1.25 to 2.26) and HR 1.27 (95% CI 1.11 to 1.45), respectively. In elevated-stable trajectory, the risk increased for all cancer (HR 1.08, 95% CI 1.02 to 1.16) and stomach cancer (HR 1.24, 95% CI 1.07 to 1.43). Significant associations were also found in the elevated-high group with oral (HR 2.13, 95% CI 1.01 to 4.47), liver (HR 1.50, 95% CI 1.08 to 2.08) and pancreatic cancer (HR 1.99, 95% CI 1.20 to 3.30).
Conclusions Our study highlights that the uncontrolled high glucose level for many years may increase the risk of cancer.
- Public Health
- Blood Glucose
- Risk Factors
Data availability statement
Data are available upon reasonable request. The HEALS dataset is available on request from the National Health Insurance Sharing Service, https://nhiss.nhis.or.kr/.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
It is well established that diabetes mellitus (DM) elevates the risk of cancer. However, most published data often required reinterpretation because DM is not a single disease, but a group of metabolic disorders characterized by chronic hyperglycemia.
Hyperglycemia may appear at the stage of pre-diabetes. Further studies are needed to investigate the impact of hyperglycemia, including fasting blood glucose (FBG) levels, on cancer incidence and to develop preventive strategies.
There exists a major research gap because most studies measured FBG levels only at baseline, which may not reflect FBG levels over an extended period and may lead to inaccurate estimations.
WHAT THIS STUDY ADDS
Our study highlights that the uncontrolled high glucose level for many years may increase the risk of cancer, particularly in the case of digestive cancers.
FBG trajectories may be a predictive factor for the development of cancer.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
Further studies are needed to consider longer follow-up times, larger sample sizes, and more frequent FBG measurements to build a more comprehensive understanding of this area.
Introduction
Diabetes represents an important problem for socioeconomic and public healthcare systems on a global scale, and its burden is anticipated to grow.1 In South Korea, diabetes is a major contributor to the cause of morbidity and mortality,2 impacting approximately 6 million Korean adults as of 2021.3 The association between cancer and diabetes mellitus (DM) has been extensively studied; however, not all studies have yielded consistent findings regarding increased cancer risk in individuals with DM.4–6 Most published data often required reinterpretation because of the truth that DM is not a single disease, but a group of metabolic disorders characterized by chronic hyperglycemia.4 There is a long-standing controversy regarding the underlying mechanisms connecting hyperglycemia and cancer.7–9 In addition, each type of DM has different mechanisms in affected individuals.10 Furthermore, numerous factors that contribute to cancer incidence, such as obesity, diet, physical exercise, medication usage, have been identified.4 Consequently, it is inappropriate to consider patients with DM as a homogeneous group.4
As an indication of elevated blood glucose, hyperglycemia may appear at the stage of pre-diabetes or impaired fasting glucose (IFG), in which the level of blood glucose during fasting is higher than the normal stage but does not meet the diagnostic criteria for DM (FBG is between the range of 5.6–6.9 mmol/L and 100–125 mg/dL).11 This suggests that pre-diabetes or IFG may potentially elevate the risk of cancer development. The age-adjusted prevalence of IFG was 25.0% in 2013 among Korean adults.2 Additionally, approximately 15.8 million Korean adults are estimated to have pre-diabetes and are at risk of developing DM in 2021,3 indicating a potential threat posed by hyperglycemia in increasing the risk of cancer. Therefore, further studies are needed to investigate the impact of hyperglycemia, including FBG levels, on cancer incidence and to develop preventive strategies.
Many prospective cohort studies have established a connection between FBG and overall cancer risk,12 especially gastrointestinal cancer, including stomach,12 13 colorectal,14 liver,15 and pancreas,16 which can be mediated by chronic hyperglycemia,17 18 hyperinsulinemia and elevated level of bioavailable IFG.19 However, FBG levels in individuals may change during follow-up, and most studies measured FBG only at baseline, which may not reflect FBG levels over an extended period and lead to inaccurate estimations. In particular, patients with DM tend to actively manage their blood sugar levels by using glucose-lowering agents or lifestyle modifications, which could introduce bias without accounting for these factors. Considering those limitations, our study will examine FBG trajectories to capture changes in FBG levels over a follow-up period.
Trajectory analysis has recently been applied in epidemiological research to identify similar longitudinal changes over time.20 Among the different trajectory approaches, group-based trajectory modeling (GBTM) is often selected by researchers owing to its efficiency, fewer errors, and easy of explanation because it is less complex.21 However, few studies have employed GBTM to track FBG trajectories as well as examine their association with cancer risk based on FBG patterns.22 23 Given this evidence, it would be meaningful to use GBTM to investigate the link between FBG trajectory and the risk of cancer.
In this research, we conducted a population-based cohort study to investigate FBG trajectories in the Korean population using the GBTM for the purpose of investigating the relationship with several cancer types in both sexes.
Methods
Data source and study population
Our analysis used data derived from a population-wide cohort at a national level from the National Health Insurance Service–Health Screening Cohort (NHIS-HEALS) database. The NHIS is a single mandatory system covering Korea’s entire population and provides a general health screening program to all beneficiaries every 1 or 2 years. The NHIS-HEALs is a cohort of randomly selected 10% of screening participants during 2002–2003 and followed up until 2015.24 The cohort included 514,866 people aged 40–79 years with essential information, including sociodemographic factors, self-reported health behaviors, clinical laboratory results, and healthcare usage based on insurance claim data.24
From the cohort entry 2002–2003, 264,777 subjects who underwent health examinations in three waves, 2002–2003, 2004–2005, and 2006–2007, were included in further analyses. Additionally, 7448 patients were diagnosed with cancer and 539 patients who died before 2008 were excluded to create a cancer-free cohort study. Furthermore, 519 individuals with missing age, sex, and FBG data were excluded. Consequently, our final study cohort comprised 256,271 individuals (figure 1).
Exposure and covariates
Fasting blood glucose
As defined by WHO, FBG levels are measured by taking a blood sample from participants who had fasted for a minimum of 8 hours,25 measured mmol/L or mg/dL. According to the American Diabetes Association, there are three levels of FBG: normal, pre-diabetes, and diabetes. The pre-diabetes stage includes individuals with FBG concentration of 100–125 mg/dL or 5.6–6.9 mmol/L. Individuals were considered to have diabetes if they had higher FBG levels than the pre-diabetes criterion (FBG level ≥126 mg/dL or 7 mmol/L). On the other hand, they were considered to be in the normal stage (FBG level <100 mg/dL or 5.6 mmol/L).11
Covariates
Data on covariates, including sex, age, income level, body mass index (BMI), smoking behaviors, alcohol consumption, exercise, Charlson Comorbidity Index (CCI), family history of DM, presence of DM, use of antidiabetic medications, and presence of chronic viral hepatitis B and C infection, were collected from the first wave in 2002–2003. CCI was treated as continuous variable, while other covariates were considered categorical including age (40–49, 50–59, 60–69, 70–79 years old), income (four quartiles, using income-based insurance premiums as a proxy for income), smoking behaviors (never, former, and current smoker), drinking behaviors (rarely, 2–3 times/month, 1–2 times/week, 3–4 times/week, and almost every day), BMI (underweight: <18.5, normal: 18.5–22.9, overweight: 23–24.9, and obese: ≥25 kg/m226), exercise (no physical activity, 1–2, 3–4, 5–6 times/week, and every day), presence of DM (yes or no), family history of DM (yes or no), and antidiabetic medications (yes or no). Presence of chronic viral hepatitis B and C infections (yes or no) were adjusted to analyze liver cancer status.
DM diagnosis was determined as meeting one of the two following criteria or both: having a diagnosis code of DM (International Classification of Diseases 10th edition (ICD-10) codes), specifically code E10 for type 1 diabetes, E11 to E14 for type 2 diabetes and the prescription of any diabetic agents, or having FBG level ≥126 mg/dL.27 The use of antidiabetic medication was determined if an individual had at least one claim per year for the relevant prescription.27
Case ascertainment
Incident cancers were defined by using primary diagnosis, validated with high sensitivity (>90%) and accuracy (>80% consistency with cancer registry data).28 An individual who was an inpatient once or outpatient three or more times (ICD code: C00-C97) within the first year of the claim was considered a newly diagnosed patient with cancer. The date of cancer diagnosis was determined as the date on which cancer was first detected.28 We determined the incidence of 17 specific types of cancer for both sex in the NHIS-HEALS database (according to the GLOBOCAN cancer dictionary), all cancer combined and gastrointestinal cancer (ICD code: C15–C25) including esophagus, stomach, colon, rectum, anus, liver, gallbladder, and pancreatic.
Statistical methods
Group-based trajectory modeling
To calculate changes in FBG, we employed GBTM, a statistical method that determines separate trajectories within a population.20 To adhere to best practices in latent modeling, which typically require data from multiple time points—at least three occasions, we monitored FBG changes across three waves spanning from 2002 to 2007. The package ‘lcmm’ in R was used for GBTM analysis.
FBG level was treated as a continuous variable. To determine the best model, we initially checked all the quadratic models from one to six groups, considering that previous studies identified up to five FBG trajectories.22 23 The optimum number of groups was built on the largest Bayesian Information Criterion (BIC).20 However, BIC increases when more groups were added to the model. The best-fit model was then identified based on the percentage of group membership (≥1%), model parsimony, as well as distinctiveness (online supplemental table 1). The model with four groups was finally chosen after accessing all the above criteria.
Supplemental material
After selecting the number of groups, we determined the polynomial of each group by comparing models with different functional forms based on the significant level (p value <0.05), beginning with the highest polynomial of each group. Finally, the model with four quadratic trajectories was chosen as the best-fitted model (figure 2): 87.78% (N=224,957) of individuals maintained a low FBG concentration (mean FBG at each wave was smaller than 100 mg/dL, referred to as ‘low-stable group’); 7.17% (N=18,385) of individuals kept a moderate FBG concentration (mean FBG at each wave ranged from 113 to 124 mg/dL, referred to as ‘elevated-stable group’); 1.82% (N=4667) of participants initially were at pre-diabetes stage and then progressed to the diabetes stage (mean FBG at each wave ranged from 104 to 166 mg/dL, referred to as ‘elevated-high group’); and 3.22% (N=8262) of participants maintained a high FBG concentration (mean FBG at each wave was greater than 177 mg/dL, which is called as ‘high-stable group’).
Cox proportional hazard models
To compute HRs and 95% CIs for the FBG index relative to the reference group, we used Cox proportional hazard models. The low-stable group served as the reference for comparison with the other trajectories. The censored cases comprised individuals who either died or did not experience an event of interest during the period from 2008 to 2015. The time-to-event (measured in years) was defined as the period starting from January 1, 2008, and ending at either the date of cancer diagnosis date, censoring date, or the end of the following period on December 31, 2015, whichever event came first.
For confounding adjustments, we included the following covariates: sex, age, income level, smoking behaviors, alcohol consumption, BMI, exercise, family history of DM, CCI, DM diagnosis, and antidiabetic medications. The presence of chronic viral hepatitis B and C infections (ICD-10: B18.*) were adjusted to assess liver cancer status.
Cox models were stratified according to types of DM and medication at baseline by including information of DM diagnosis as well as using antidiabetic drugs, and then adjusting for covariates. Analyses of the Cox proportional hazard model were performed using SAS Enterprise Guide V.7.1 (SAS Institute, Cary, North Carolina, USA).
Results
General characteristics of the participants and cancer incidence
The mean age of the 256,271 participants included was 51.27 years (SD=8.87). The mean BMI (kg/m2) was 23.97 (SD=2.89), and approximately 62.26% of the participants had excess weight (BMI ≥23 kg/m2). More than half of the participants (52.31%) did not engage in any physical activity. Approximately 62.81% of the participants were non-smokers; around 52.68% reported rarely drinking at the baseline. Nearly 6.26% of participants had a family history of DM (table 1).
Participants in the low-stable group accounted for the highest proportion of the age between 40 and 49 years old. In contrast, a larger percentage of DM, family history of diabetes, and medication users were observed from examinees in the high-stable group than in the other groups (table 1).
By December 31, 2015, 18,991 cancer cases were reported. Stomach cancer (17.20%), colorectal cancer (13.44%), lung cancer (10.82%), and thyroid cancer (12.61%) account for the highest rates of cancers (online supplemental table 2).
Cancer risks according to FBG trajectories
Relative to the low-stable group, HR adjusted for all covariates in the elevated-stable group was 1.08 (95% CI 1.02 to 1.16) for all cancers combined. For each type of cancer, the risk of developing multiple myeloma and malignant plasma cell neoplasms was the strongest in the high-stable group (HR 4.09 (95% CI 1.40 to 11.95)), followed by oral and pancreatic cancer (HR 2.13 (95% CI 1.01 to 4.47) and HR 1.99 (95% CI 1.20 to 3.30), respectively). The risk increased by 50%–68% for liver cancer in elevated-high and high-stable groups. We found significant associations in stomach cancer in elevated-stable (HR 1.24 (95% CI 1.07 to 1.43)) and high-stable group (HR 1.25 (95% CI 1.00 to 1.57)). Significant associations were also observed in gastrointestinal cancer in elevated-stable (HR 1.22 (95% CI 1.12 to 1.34)) and high-stable group (HR 1.27 (95% CI 1.11 to 1.45)) (table 2) (online supplemental table 3).
Regarding the results of subgroup analysis, we divided DM diagnosis into two types: type 1 and type 2. Approximately 85.98% of individuals did not have diagnosis of DM. Most patients with DM had type 2 diabetes, whereas only 2.67% had type 1 diabetes. Of patients with type 1 diabetes, 78.69% were in the low-stable group, indicating that these patients adhered to use medications and were able to control their FBG level. In contrast, approximately 18.61% of patients with type 2 diabetes maintain normal FBG levels. A significant association was observed in patients without DM in the elevated-stable group (HR 1.09 (95% CI 1.01 to 1.18)). There was no significant difference in cancer risk in the other trajectory groups and in the DM group (online supplemental table 4).
Glucose-lowering agents were divided into seven subgroups existing in 2002–2003 NHIS-HEALS data: insulin, sulfonylureas, metformin, meglitinide, thiazolidinedione, acarbose, and combined therapy which included two or more classes of antidiabetic agents. Sulfonylureas and metformin were the most commonly used drugs for treatment (>87%). No significant association was observed among medication users except for participants taking medications from acarbose in the elevated-high group. In this group, the risk of developing cancer was approximately 2.5-fold higher than that of the reference group (HR 2.54 (95% CI 1.04 to 6.19)) (online supplemental table 5).
Discussion
We applied trajectory analysis to examine changes in FBG levels at three different time points between 2002 and 2007, which enabled us to identify four distinct trajectory groups. The most prevalent group was the low-stable group, which constituted 87.78% of the sample population. The elevated-stable, high-stable, and elevated-high groups accounted for 7.17%, 3.22%, and 1.82% of the population, respectively. Notably, previous research exploring FBG trajectories reported varying numbers of trajectories, typically ranging from four to five.22 23 29 Similarly, a prior study conducted in Korea by Jeon et al also investigated FBG trajectories and identified four distinct patterns.23 However, the distribution and proportion of each trajectory group in our study differed from those observed in the aforementioned study. Specifically, four groups included low-stable, moderate-stable, elevated-upward, and high-upward with the percentage of each group were 47.9%, 44.1%, 6.7%, and 14.0%, respectively.23 This difference is due to dissimilarities in the targeted population, sample size, and follow-up duration. Notably, our study used a larger sample size and included participants aged between 40 and 79 years who were at high risk, whereas the study of Jeon et al recruited examinees starting from 20 years of age.23
Our analysis identified positive relations between FBG trajectories and the risks of overall cancers, gastrointestinal cancer, as well as several specific cancers (multiple myeloma and malignant plasma cell neoplasms, oral, stomach, liver, and pancreatic cancer). A previous study in China also explored FBG trajectories and their association with overall cancers and digestive cancers.22 Significant associations with all cancers and digestive cancers were observed in low-stable and elevated-stable groups.22 However, because of the limited number of cancer cases, a detailed analysis of site-specific cancers could not be extensively conducted.22
An elevated risk for overall cancers was observed in the elevated-stable group, which provides a potential suggestion for further investigations into the prediabetic population in the future. While past research consistently indicates a positive association between a history of diabetes and cancer risk, suggesting that the high-stable group may be closely linked to an elevated risk of total cancer, our study produced conflicting results.4 5 This discrepancy may be attributed to the predominance of diabetic individuals in this group, comprising over 90% of the sample. While both adjusted HR for all cancer combined in model 1 and model 2 in table 2 showed significant results, this became non-significant after adjusting for diabetes diagnosis. Regarding increase cancer risk in patients with elevated FBG, a prospective cohort study found a significant association between FBG levels ranging from 110 to 125 mg/dL and cancer risk in males, but not in females.30 Other studies have also documented an elevated risk of cancer in prediabetic individuals, accompanied by higher mortality from certain types of cancer.12 19 31 Several biological mechanisms underlie this association. Hyperglycemia-induced accumulation of advanced glycation end-products and the generation of oxidative stress may promote cancer development.32 Furthermore, hyperinsulinemia and increased level of bioavailable IGF, which are associated with insulin resistance,33 may contribute to cancer cell proliferation.34
Few studies may evaluate the association with specific cancers such as multiple myeloma and malignant plasma cell neoplasms because of insufficient numbers of cancer cases. In this study, we observed a strong positive association with high-stable group for these specific cancers; however, the estimations were based on the small number of cases. To our knowledge, positive associations between FBG and oral cancer have not been described previously, which is required for further investigation.
Gastrointestinal cancers, including pancreatic, stomach, and liver cancer, are widely acknowledged to be closely linked to elevated FBG levels.13 15 16 In line with previous research, our study revealed significant findings in these specific cancer types, although the patterns varied among different groups. The potential results may be contributed by the high number of patients with diabetes in high-stable group. Furthermore, some individuals may have diabetes but may not have sought medical attention, leading to an underdiagnosis despite displaying elevated FBG level in those patients. These may affect to the estimation in final model after adjusting for diabetes diagnosis covariate.
Our subgroup analysis revealed a high overall cancer risk in the elevated-stable group among subjects without a diagnosis of DM. The risk increased in individuals with elevated FBG levels but did not meet the diagnostic criteria for DM, according to the current American Diabetes Association definition of pre-diabetes.11 It suggests that individuals with persistently uncontrolled high glucose levels over an extended period may have an increased risk of cancer. This finding has important public health implications because, in the Korean population, the prevalence of pre-diabetes tended to increase over time from 21.5% in 2006 to 25.0% in 2013.2 Considering the substantial prevalence of pre-diabetes and its association with cancer risk, implementing interventions targeting this large population could have a major impact on public health.
The cancer risk among the diabetic population may be influenced by certain therapies used for diabetic treatment, potentially through their impact on insulin levels.35 In the case of acarbose, treatment did not show significant differences in cancer incidence in previous studies.35 Contrastingly, we did not find a significant association between cancer risk and most types of therapy, except for that with acarbose. We observed that patients with DM in the elevated-high group had a greater risk of cancer than those in the low-stability group, even when using acarbose to manage their blood sugar levels. Acarbose blocks the breakdown of starches in the intestine, thereby slowing the rise in blood glucose levels after a meal.11 However, acarbose is not a primary choice in clinical practice, following the guidelines for using glucose-lowering medication use,36 which may explain the limited number of cases using acarbose in our research. Thus, our evaluation was limited and lacked sufficient statistical power, resulting in divergent findings compared with the findings of previous research.
In our analysis, cancer risk did not increase when treating metformin in all FBG patterns. Numerous studies have suggested that metformin plays a protective role in cancer development. A nested case-control study involving female users of oral antidiabetic drugs conducted in the UK reported that prolonged use of metformin linked with a reduced risk of developing breast cancer, compared with non-use of metformin.37 Additionally, each gram added to the metformin dose has been associated with a 42% decline in cancer mortality.37 The mechanism is far from clear, but it has been observed that adenosine monophosphate-activated protein kinase (AMPK)—activated by metformin— leads to a strong suppression of cell proliferation in malignant as well as non-malignant cells.38 This effect may be mediated through cell-cycle regulation and inhibition of protein synthesis.38
Our study had several notable strengths. First, we used data from the NHIS-HEALS, which covered a large sample size and a long follow-up duration (7.6 years). Second, we employed trajectory analysis, specifically the GBTM, to identify the FBG patterns. Moreover, this study was based on measuring FBG, rather than self-reported information, which helped to mitigate recall bias or misclassification. Furthermore, our study stands out as one of the few to incorporate medication use as a covariate in the analysis.
However, this study had some limitations. The first was the relatively short timeframe over which FBG changes were measured (6 years—three waves). Although no existing publications have captured lifelong FBG measurements, this limitation should be acknowledged. It is important to note that certain cancer types, such as breast, prostate, or testicular cancer, were not available in the NHIS-HEALS dataset owing confidentiality, thus limiting our analysis in these specific cancer categories. However, this is one of the few studies to assess cancer risk related to FBG trajectories. Third, associations with DM diagnosis as well as antidiabetic medications may not have been fully evaluated, as the limited number of cases after subgroup analysis could not provide enough statistical power. Fourth, although some important covariates were adjusted, residual confounders may exist, for example, family history of cancer. Fifth, our dataset lacked certain tests that reflect long-term glucose metabolism, such as HbA1c or oral glucose tolerance tests, which may result in misclassifications within the pre-diabetes and diabetes groups. Lastly, our cohort study’s population may not fully represent the entire Korean population because it included participants aged greater than 40 who underwent multiple check-ups between 2002 and 2007. Thus, the characteristics of patients who participated in screening programs may differ from those of patients who did not participate.
In conclusion, our analysis will contribute to the understanding of the link between cancer risk and changes in FBG levels. Our findings revealed relationships between FBG patterns and the risk of overall cancer, especially gastrointestinal cancer, and several specific types of cancer such as multiple myeloma and malignant plasma cell neoplasms, oral, liver, stomach, and pancreatic cancer. These results suggested that individuals in the pre-diabetes stage may be at an increased risk of developing cancer. Future studies should consider longer follow-up times, larger sample sizes, and more frequent FBG measurements to gain a more comprehensive understanding of this area.
Data availability statement
Data are available upon reasonable request. The HEALS dataset is available on request from the National Health Insurance Sharing Service, https://nhiss.nhis.or.kr/.
Ethics statements
Patient consent for publication
Ethics approval
This study protocol was approved by the Institutional Review Board of the National Cancer Center, Korea (NCC2022-0180) and the Korea NHIS Medical Information Disclosure Committee (NHIS-2022-2-185). As this study used anonymous secondary data, informed consent was waived.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Contributors The contribution of the authors were as follows: Guarantor: J-KO, TMTK. Conceptualization: J-KO, TMTK. Methodology: TMTK, TTB. Formal analysis: TMTK. Investigation: J-KO, TMTK. Writing—original draft preparation: TMTK. Writing—review and editing: J-KO, HYK, EP. All authors approved the final version of the manuscript.
Funding This work was supported by the National Cancer Center (grant number: NCC-2210862). The funder had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
Author note This study was designed and conducted according to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.