Your privacy, your choice

We use essential cookies to make sure the site can function. We also use optional cookies for advertising, personalisation of content, usage analysis, and social media.

By accepting optional cookies, you consent to the processing of your personal data - including transfers to third parties. Some third parties are outside of the European Economic Area, with varying standards of data protection.

See our privacy policy for more information on the use of your personal data.

for further information and to change your choices.

Skip to main content

Risk factors of immune-related endocrine toxicities in non-small cell lung cancer patients treated with pembrolizumab and its impact on patient outcomes: a multicenter retrospective study

Abstract

Background

Pembrolizumab has been approved as a first-line treatment for non-small cell lung cancer (NSCLC) patients. However, a percentage of patients discontinue immunotherapy due to immune-related adverse events (irAE). Among these events, immune-related endocrine toxicities (E-irAE) represents the most common form, though their etiology, risk factors, and impact on patient outcomes remain poorly understood.

Materials and methods

This retrospective cohort study was conducted across 5 multiple centers to investigate the outcomes of NSCLC patients who received pembrolizumab treatment between October 1, 2019, and September 30, 2023. E-irAE can occur on the thyroid, pituitary, adrenal glands, pancreas, and parathyroid. So thyroid function, adrenocorticotropic hormone, cortisol, sex hormone and glycaemia were measured at baseline and at regular intervals after the initiation of pembrolizumab treatment.

Results

Our study included a total of 380 NSCLC patients treated with pembrolizumab, 114 patients (30.00%) developed E-irAE. Among them, 107 patients (93.86%) developed immune-related thyroid dysfunction (irTD) (5 cases of combined immune-related hypophysitis (IH)), 4 patients (3.51%) only developed IH, and 3 patients (2.63%) developed type 1 diabetes mellitus. IrTD was found to be independently associated only with monocyte-to-lymphocyte ratio (MLR) (odds ratios (OR) = 0.060, 95% CI 0.000-0.375; p = 0.015) and anti-thyroglobulin antibody (TGAb) (OR = 31.898, 95% CI 1.516-671.367; p = 0.026). Kaplan-Meier Survival Analysis showed that the progression-free survival (PFS) was significant longer in stage IV NSCLC patients with irTD than in those who did not (44.72 weeks vs. 27.79 weeks; hazard ratio (HR) = 0.645, 95% CI 0.440–0.946; p = 0.025), particularly in the subgroup of subclinical hypothyroidism (HR = 0.567, 95% CI 0.324–0.994; p = 0.047). We also found that sex (HR = 0.493, 95% CI 0.291–0.834; p = 0.008) was identified as an independent factor associated with better PFS.

Conclusion

E-irAE are recognized as prevalent common irAE with various phenotypic manifestations. Low MLR and positive TGAb at baseline have been identified as risk factors that increase the likelihood of developing irTD. Sex and the occurrence of irTD were independently associated with improved PFS.

Peer Review reports

Introduction

According to Global Cancer Statistics 2020 [1], lung cancer has the second highest prevalence and the highest mortality rate among malignant tumors. In China, lung cancer has the highest incidence and mortality rates among cancers [2]. The two major histologic subtypes are non-small cell lung cancer (NSCLC) as well as small-cell lung cancer, accounting for approximately 80–85% and 15–20% of cases, respectively. With an enhanced understanding of tumor mechanisms, immunotherapy has become a promising treatment for malignant cancers.

In contemporary clinical practice, various immune checkpoint inhibitors (ICIs) are used, including programmed death-1/programmed death ligand-1 (PD-1/PD-L1) inhibitors and cytotoxic T lymphocyte-associated antigen-4 inhibitors. Pembrolizumab, a humanized monoclonal antibody and one of the most common PD-1 inhibitors [3], has been approved by the Food and Drug Administration and the National Medical Products Administration for advanced NSCLC. Clinical trials such as KEYNOTE-010 [4], KEYNOTE-024 [5], and KEYNOTE-042 [6] have demonstrated that pembrolizumab has superior outcomes in terms of overall survival (OS) and progression-free survival (PFS) compared to chemotherapy. Furthermore, other clinical trials, including KEYNOTE-021 [7], KEYNOTE-189 [8], and KEYNOTE-407 [9] have shown that the combination of pembrolizumab and chemotherapy also results in prolonged OS and PFS in NSCLC patients.

With increased utilization and indications, immune-related adverse event (irAE) has attracted considerable attention due to their high incidence and lethality in severe cases. IrAE can occur in any organ or tissue, immune-related endocrine toxicities (E-irAE) were the most common form of irAE. Importantly, immune-related thyroid dysfunction (irTD) is the most common E-irAE, while immune-related hypophysitis (IH), immune-related diabetes mellitus (DM), immune-related primary adrenal insufficiency (PAI) and hypoparathyroidism occur less frequently. Several investigations have examined potential risk factors associated with the development of irTD, including gender, age, body mass index (BMI), and antibodies. However, there is still controversy regarding the factors that contribute to the development of irTD, as different studies have reached disparate conclusions. Additionally, the association between the occurrence of irTD and the prognosis is still controversial. Therefore, we conducted further research to explore the incidence of irTD, potential risk factors, and the correlation between irTD and PFS among NSCLC patients treated with pembrolizumab in real-world settings.

Materials and methods

Patients and study design

This retrospective study was approved by the local medical ethics committees in 5 centers. The study included patients who received pembrolizumab treatment between October 24, 2019, and September 30, 2023. Relevant clinicopathological characteristics were collected, including age, sex, Eastern Cooperative Oncology Group performance status (ECOG PS), BMI, neoplasm stage, histological type, PD-L1 expression by immunohistochemistry, prior treatment, and baseline peripheral blood parameters.

Since E-irAE mainly include irTD, IH, ICI-DM, ICI-PAI and immune-related hypoparathyroidism, and the prevalence of these disorders is very high, we monitored thyroid function tests (TFTs), thyroid autoantibodies, ACTH, cortisol, and sex hormone, including prolactine (PRL), follicle-stimulating hormone (FSH), luteinizing hormone (LH), estradiol (E2), progesterone (P), testosterone (T) at baseline and every two cycles during treatment. Glycaemia and glycosylated hemoglobin were performed at baseline and every treatment cycle, if a patient develops hyperglycemia, we would evaluate pancreatic insulin reserve determined by C-peptide levels and diabetes related antibody (Ab) (glutamic acid decarboxylase anti-GAD-, insulin Ab-Anti-IA2- and islet cell Ab-ICA-). The diagnosis of IH was supported radiographically with an MRI of the pituitary gland.

The inclusion criteria included: (1) patients with pathologically and/or cytologically confirmed NSCLC; (2) patients with the administration of at least two doses of pembrolizumab, and a minimum of two measurements of tests following the initiation of treatment; (3) patients without endocrine-related disease before starting ICI therapy; Patients were excluded: (1) patients with endocrine malignancies or metastases; (2) patients had previous treatment of any other ICIs; (3) patients with endocrine disease requiring medication at baseline; (4) patients had incomplete thyroid function assessment at baseline.

Definitions

IrTD refers to the identification of new-onset thyroid dysfunction subsequent to immunotherapy, excluding other potential causes. It was categorized into four groups based on the biochemical pattern of thyroid dysfunction. Overt hypothyroidism (Oho) referred to a Thyroid-stimulating hormone (TSH) level above the upper limit of the normal reference range and falling between 10.0 µIU/mL, with FT4 level below the reference range. Subclinical hypothyroidism (Scho) was defined as a elevated TSH level but FT4 level within the reference range. Overt hyperthyroidism (Ohe) referred to a TSH level below the lower reference range with FT4 level above the upper reference interval. Subclinical hyperthyroidism (Sche) was defined as a TSH level below the lower reference range with a normal FT4 level.

IH is a relatively rare irAE, characterized by a lack of specificity in clinical manifestations. Its diagnosis is presumptive, relying on clinical symptoms (e.g., headache and asthenia) and/or hyponatremia and/or at least one pituitary deficit and/or abnormal imaging [10]. DM typically occurs as rapid-onset hyperglycemia, and has a high risk of progressing to diabetic ketoacidosis in absence of prompt diagnosis and treatment [11]. The clinical manifestations of PAI often lack of specificity, it’s typically characterized by low or relatively low plasma cortisol levels with elevated adrenocorticotrophic hormone (ACTH) [12]. Immune-related hypoparathyroidism can induce hypocalcemia, accompanied by nausea, vomiting, paresthesia, abdominal pain, and fatigue [12], it can be diagnosed with the level of parathyroid hormone.

The laboratory cut-offs utilized in hospitals were as follows: reference ranges of 0.35–4.94 mIU/L for TSH, 9.01–19.05 pmol/L for free thyroxine (FT4), 2.63–15.70 pmol/L for free triiodothyronine (FT3) levels, 0-4.11 IU/mL for thyroid peroxidase antibody (TPOAb), and 0-5.61 IU/mL for anti-thyroglobulin antibody (TGAb). Patients who initially experienced transient Ohe/Sche but subsequently developed hypothyroidism during the follow-up period were classified as the Oho/Scho group.

The time to onset of E-irAE was defined as the duration between the administration of the initial dose of ICI therapy and the first documented abnormality in the monitored hormones. PFS was determined as the time interval between the initiation of pembrolizumab and the first occurrence of disease progression or death.

Statistical analysis

Continuous quantitative variables with skewed distribution were compared using the Wilcoxon test, and summarized using the interquartile range (IQR). Categorical variables were compared using Fisher’s exact test or the Chi-square test. Univariate logistic regression analysis was used to evaluate the baseline characteristics and biochemical associations with irTD. Factors with a p-value < 0.05 were included in a multivariate logistic regression model, and odds ratios (ORs) with their 95% confidence intervals (CIs) were presented. PFS for different groups was calculated using the Kaplan-Meier method, and the log-rank test was used for comparison between groups. Univariate and multivariate Cox proportional hazard regression models were used to estimate hazard ratios (HRs) and their 95% CIs. Prespecified subgroup analyses were conducted for patients with Oho, Scho, Ohe, and Sche, compared to patients without irTD. All statistical analyses were performed using IBM SPSS (version 25.0), and a p-value of < 0.05 was considered statistically significant. Figures were created using GraphPad Prism (version 9.5.1).

Results

Patient characteristics

A total of 479 patients with NSCLC who received pembrolizumab treatment in 5 medical centers were included. We excluded 99 patients from the analysis for various reasons, including 5 patients with a history of previous treatment with other ICIs, 41 patients with pre-existing endocrine disease prior to pembrolizumab therapy, 28 patients who lacked TFTs results, and 25 patients who received only a single dose of pembrolizumab. Therefore, a total of 380 patients were enrolled in our study (Fig. 1). The median age of the enrolled patients was 66.00 (59.00-70.25) years, with a majority of male participants (92.1%). Most patients had an ECOG PS of 0–1 (85.0%), and nearly half had a standardized BMI. The majority of NSCLC cases (62.1%) were classified as squamous cell carcinoma, with most patients diagnosed at stage IIIB + IIIC and stage IV (32.1% and 55.3%, respectively). PD-L1 expression data were available for 142 patients, with high expression (≥ 50%) observed in 13.4% of cases. The majority of patients (87.4%) received pembrolizumab as first-line therapy (Table 1).

Fig. 1
figure 1

Flowchart of patients selected

Table 1 Clinical characteristics in NSCLC patients

Incidence, onset and outcomes of IrTD

Within the follow-up, irTD was observed in 107 (28.16%) of all patients. Among these cases, the most prevalent form of irTD was Sche (44.86%), followed by Scho (34.58%), Oho (16.82%), and Ohe (3.74%). The occurrence of irTD was typically observed within weeks to months after the initiation of pembrolizumab (Fig. 2). Interestingly, the onset of hyperthyroidism was earlier than the onset of hypothyroidism (9.43 [6.86–19.71] weeks vs. 33.79 [17.68–69.47] weeks; p < 0.001). The time to Oho onset ranged from 7.43 to 75.57 weeks, for Scho it ranged from 3.29 to 203.29 weeks, for Ohe it ranged from 5.86 to 71.14 weeks, and for Sche it ranged from 2.14 to 86.86 weeks. The median time to Oho onset was 22.64 (15.54–27.57) weeks, while for Scho it was 34.86 (16.00–67.00) weeks, for Ohe it was 8.21 (6.39–25.18) weeks, and for Sche it was 9.79 (6.71–15.11) weeks. Among the 55 patients who presented with hypothyroidism, 25.45% experienced transient hyperthyroidism before the onset of hypothyroidism, with a median time of 8.29 (5.79–13.22) weeks. Notably, patients with Ohe (6/10) demonstrated a higher risk of developing subsequent hypothyroidism compared to those with Sche (5/53).

Fig. 2
figure 2

Relationship between time of irTD onsets and incidence of irTD

Overall, the outcomes of patients with irTD were favorable, with TSH levels returning to normal levels in 5 cases of Oho; while 19 cases of Scho, 2 cases of Ohe, and 34 cases of Sche also showed normalization of TSH levels. Regarding FT4 levels, 11 cases of Oho and 3 cases of Ohe achieved normalization. However, it should be noted that the majority of cases still exhibited abnormalities in autoantibodies. Permanent hypothyroidism requiring levothyroxine replacement therapy was more prevalent among patients with Oho (38.9%) compared to those with Scho (5.4%). Importantly, none of the patients discontinued their medication due to severe irTD.

Incidence, onset and outcomes of IH

Among the 380 patients, IH was the second common type of E-irAE (9 patients; 7.89%), the median time to IH onset was 42.57 (25.14–62.29) weeks. Most patients were squamous cell carcinoma (7 patients; 77.78%), with a mean age of 66.12 years and all of them had ACTH deficiency. 5/9 patients had concurrent irTD, but only case 9 required replacement treatment with levothyroxine. All patients presented with fatigue, followed by headache (77.78%) as the second most common syndrome. Interestingly, case 3 exhibited increased urine output, while case 7 had low blood pressure. Notably, none of them had abnormal findings on pituitary MRI (Table 2).

Table 2 Patients with immune-related hypophysitis during follow up

All patients diagnosed with IH needed replacement treatment with glucocorticoids, and their outcomes were good, none of them required permanent discontinuation of medication due to severe IH. Specifically, all patients achieved normalization of ACTH levels, while only 3 cases demonstrated restoration of cortisol levels to within the normal range (Table 2).

Incidence, onset and outcomes of DM

There were 3 cases (2.63%) of newly developed DM during pembrolizumab treatment in our study, with a median onset time of 30.43 (20.64–34.93) weeks. Furthermore, their fasting blood glucose (FBG) level were found to be higher than 13.9 mmol/L, and C-peptide levels were consistently below the normal range. Interestingly, their HbA1c levels were only slightly elevated above the normal range. All 3 cases exhibited fulminant type 1 diabetes, with case 3 experiencing the additional complication of diabetic ketoacidosis (DKA). Due to two cases being graded as Grade 3 and one case as Grade 4 (according to CTCAE 5.0), all three patients temporarily discontinued pembrolizumab treatment and received symptomatic management with insulin and oral hypoglycemic agents. However, pembrolizumab treatment was restored in all cases once they achieved stable blood glucose control (Table 3).

Table 3 Patients with immune-related diabetes mellitus during follow up

Clinical features and biomarkers

No statistically significant differences were observed between the irTD group and euthyroid group in terms of age, sex, ECOG PS, BMI, histological features, neoplasm stage, PD-L1 TPS, treatment line, and prior treatment. Analysis of baseline peripheral blood results revealed several correlations with the development of irTD. Specifically, monocyte-to-lymphocyte ratio (MLR) (IQR = 0.36[0.26–0.51]; Z = 2.255, p = 0.024), FT3 levels (IQR = 4.12[3.75–4.53]; Z = 1.975, p = 0.048), positive expressions of TPOAb at baseline (positive: 11/76 vs. 8/173; x2 = 7.286, p = 0.007) and TGAb at baseline (positive: 5/19 vs. 2/46; x2 = 4.660, p = 0.031) were found to be statistically associated with the incidence of irTD (Table 4).

Table 4 Baseline peripheral blood results in NSCLC patients

Previous studies have suggested that the occurrence of irTD may be associated with chemotherapy history, baseline TSH levels [13], gender [14], age and BMI [15]. In this study, the Wilcoxon test and Chi-square test analyses revealed statistically significant differences in baseline MLR, FT3, TPOAb, and TGAb at baseline between groups (p < 0.05), and therefore the presence or absence of a history of chemotherapy, baseline MLR, TSH, FT3, FT4, gender, age, BMI, TPOAb, and TGAb were finally included in the binary logistic regression analysis. The analysis demonstrated that was only MLR (OR = 0.060, 95% CI 0.000-0.375; p = 0.015) and TGAb at baseline (OR = 31.898, 95% CI 1.516-671.367; p = 0.026) were associated with the occurrence of irTD (Fig. 3).

Fig. 3
figure 3

Binary logistic regression analysis for independent risk factors of irTD onset

Although there was no statistical significance when comparing the irTD and euthyroid groups, it is worth noting that TSH levels exhibited a significant correlation with the onset of irTD when comparing the different clinical phenotypes separately. Specifically, in the Oho group, TSH levels showed a positive correlation (OR = 1.874, 95% CI 1.094–3.208, p = 0.022), while in the Scho group, a similar positive correlation was observed (OR = 2.214, 95% CI 1.522–3.219, p < 0.001). In contrast, the Sche group demonstrated a negative correlation between TSH levels and irTD (OR = 0.186, 95% CI 0.090–0.386, p < 0.001). Moreover, there was no statistically significant difference in FT4 levels and treatment lines between the groups (Table 5).

Table 5 Clinical and biochemical associations with patients with pembrolizumab associated different subtypes of IrTD during follow-up

Outcomes survival

Kaplan-Meier Survival Analysis found that stage IV NSCLC patients who developed irTD during immunotherapy had a longer PFS compared to those who remained euthyroid (44.72 weeks vs. 27.79 weeks; HR = 0.645, 95% CI 0.440–0.946; p = 0.025)Fig. 4A). However, the difference in PFS between stage IIIB/IIIC NSCLC patients who developed irTD during immunotherapy was not statistically significant compared to those who did not develop irTD (HR = 0.776, 95% CI 0.436–1.381; p = 0.387)Fig. 4C). To further explore the relationship between irTD and prognosis, we subdivided irTD into subclinical and overt categories. In subgroup comparisons, PFS was improved in the Scho group compared to the euthyroid group with statistical significance in stage IV NSCLC patients (HR = 0.567, 95% CI 0.324–0.994; p = 0.047), while no significant difference was observed between the other subtype groups and euthyroid group (Fig. 4B, 4D).

Fig. 4
figure 4

Kaplan-Meier estimates of progression-free survival (PFS) in stage IV NSCLC patients with irTD after pembrolizumab therapy. (A) PFS in stage IV NSCLC patients with irTD relative the patients who remained euthyroid. (B) PFS in stage IV NSCLC patients with Oho, Scho, Ohe, and Sche compared to the patients who remained euthyroid. (C) PFS in stage IIIB/IIIC NSCLC patients with irTD was not significantly difference from the patients who remained euthyroid. (D) PFS in stage IIIB/IIIC NSCLC patients with Oho, Scho, Ohe, and Sche compared to the patients who remained euthyroid

The results of the meta-analysis of Dall’Olio [16] et al. showed that immortal time bias had an real effect on PFS in patients. Therefore, we used 6 weeks as a time point, and patients who experienced death or progression within 6 weeks were included in the euthyroid group. As the result of univariate Cox regression analysis, only sex (HR = 1.995, 95% CI 1.192–3.340; p = 0.009) was found to be associated with PFS in patients with stage IV NSCLC in our study, but no correlation was found between age, ECOG, BMI, histological features, PD-L1 TPS, number of lines of treatment, and immortal time bias and PFS (p>0.05) (Table 6).

Table 6 Univariate Cox regression analysis of independent prognostic factors of PFS for stage IV NSCLC patients

However, previous studies have shown that age, number of lines of treatment, and treatment history may all be associated with PFS in cancer patients [17], therefore we included all the above variables as covariates in the multivariate Cox regression analysis along with the occurrence of irTD, and found that sex (HR = 0.493, 95% CI 0.291–0.834; p = 0.008) and the occurrence of irTD (HR = 0.603, 95% CI 0.397–0.916, p = 0.018) were independent predictors of a favorable prognosis (Fig. 5).

Fig. 5
figure 5

Multivariate Cox regression analysis of independent prognostic factors of PFS for stage IV NSCLC patients

Discussion

E-irAE is a common adverse event induced by PD-1 inhibitors, and severe cases can pose life-threatening risks. Pembrolizumab is one of the most commonly used PD-1 inhibitors, but there is limited research on the E-irAE relevant with it in the real world, data comes mainly from the clinical trials. Our study aims to investigate the incidence, onset and outcomes of E-irAE, identify risk factors associated with the development of irTD, and examine the potential association between irTD and better PFS in patients treated with pembrolizumab in a real-world clinical setting.

Several previous meta-analyses have shown that irTD and IH are the most common E-irAE, while DM, PAI and hypoparathyroidism are rare [18, 19]. Because of the small number of events, statistical inferences were made only for irTD and descriptive analyses were made for other types of E-irAE. Interestingly, in our cohort, although we only observed 2.36% of patients developed IH, the incidence was still higher than it in previous meta-analysis [19]. Additionally, 0.79% of patients developed DM and did not observe PAI or hypoparathyroidism.

The reported incidence of irTD ranges from 10.9 to 19.5%, hypothyroidism ranges from 7.9 to 12.1% and hyperthyroidism from 4.0–7.8% [4,5,6,7,8,9, 20]. In our study, 28.2% of NSCLC patients presented with irTD, with 14.5% manifesting hypothyroidism and 13.7% manifesting hyperthyroidism. The observed frequency of irTD was higher compared to that reported in prior clinical trials but similar to findings from earlier real-world investigations [21, 22]. The above findings indicate that the incidence of irTD in the real world may be underestimated. Moreover, irTD was categorized into four subtypes, with Sche being the most prevalent, followed by Scho, Oho, and Ohe. However, this result differs from previous studies [23, 24], which could be attributed to the inclusion of solely NSCLC patients in our study, whereas prior investigations encompassed diverse tumor types.

Hyperthyroidism typically manifests within a few weeks to 6 months after the initiation of ICIs and has a median onset time that precedes that of hypothyroidism. Consistent with this finding, our study also revealed a disparity in the median onset times of hyperthyroidism and hypothyroidism ( p < 0.001). Further analysis using the Wilcoxon test showed a difference in the median time to irTD onset across different subgroups with statistical significance (p = 0.001), particularly between Sche group and either Oho group or Scho group (p = 0.023; p< 0.001, respectively). However, it is noteworthy that the median time to onset of irTD in our study exceeded that previously reported [25,26,27,28,29,30,31], possibly due to inadequate regular monitoring of thyroid function.

In our study, all patients with hyperthyroidism and the majority of those with hypothyroidism experienced self-limiting symptoms, indicating that the grading of irTD was generally mild, which was consistent with previous reports [5, 6, 8, 20, 32]. We also found that a higher proportion of patients in the Oho subgroup, compared to the Scho subgroup, required levothyroxine replacement therapy for permanent hypothyroidism (7/18 vs. 2/37).

Several articles have investigated clinical factors related to irTD, but none have demonstrated a definitive correlation. In our cohort, age, gender, and BMI did not emerge as risk factors for irTD, possibly due to differences in cancer types, potential data bias, and limited statistical power [14, 15, 23, 33], but similar to other reports on irTD in lung cancer patients treated with nivolumab [21, 22]. Although previous studies have indicated that a low neutrophil-to-lymphocyte ratio (NLR) and platelet-to-lymphocyte ratio (PLR) at baseline may be associated with the occurrence of irAEs [34, 35], we did not observe such correlations in our study. Additionally, we found a low MLR at baseline acted as an independent risk factor for irTD, specifically for Ohe. Our conclusions are similar to those of previous articles suggesting that elevated TSH levels at baseline could serve as a predictor of irTD [23, 36, 37]. However, it is worth noting that we did not find an association between Ohe and TSH levels in the present study, which may be attributed to the insufficient sample size.

The relationship between baseline positivity of TPOAb, TGAb and the occurrence of irTD remains unclear. Previous researches have yielded conflicting findings regarding the potential impact of baseline TGAb and TPOAb positivity on the development of irTD [29, 38, 39], while others have found no correlation [40, 41]. Some studies have suggested that it is the presence of baseline TGAb, rather than TPOAb, that increases the risk of irTD [13, 36], while others have proposed that positive TGAb and TPOAb after ICIs therapy may serve as predictive factors [42, 43]. However, our results demonstrate a possible predictive role of baseline TGAb positivity in relation to irTD. It is hypothesized that irTD may be caused by an antibody-independent pathway or by an unidentified antibody [41].

Although several studies have shown a significant improvement in PFS among patients with irTD [23, 44, 45], there have been inconsistent findings reported in some studies [29]. This inconsistency could potentially be attributed to limited sample sizes. In our study, we confirmed these previous observations through Cox multivariate analysis, which revealed that irTD independently contributed to a prolonged PFS (Fig. 5). However, the underlying mechanisms remain unclear and require further investigation through large-scale studies. Moreover, we found that male patients with NSCLC was another independent protective factor associated with a longer PFS.

This study had some limitations that warrant discussion. Firstly, it is important to acknowledge that our research design was a retrospective study, which may introduce inherent biases, and our sample size may not have been sufficient for subgroup analysis in NSCLC patients with E-irAE. It is noteworthy that 350 of the 380 patients included in this study were male, accounting for 92.1%. The reason for this may be caused by the following: (1) Because gene mutations such as EGFR/ALK/ROS1/RET are very common in female NSCLC in China, many female patients receive targeted therapy instead of immunotherapy; (2) The patients we included were predominantly squamous cancer (62.1%), there is a higher proportion of male than female in squamous lung carcinoma; (3) It may be related to the national situation in China, where the use of pembrolizumab is self-funded, and in the overall gender tendency will be more in favor of male. Secondly, not all patients underwent regular TFTs, sex hormone, cortisol, ACTH tests and blood glucose monitoring before each treatment with pembrolizumab. Thirdly, some patients received pembrolizumab for a short period, resulting in limited availability of post-treatment data. This restricts our ability to capture the full extent of E-irAE in these cases. Fourthly, there were instances where the assessment of antibodies, a possible predictor of irTD, was missing. This limitation may be attributed to insufficient recognition and monitoring of these predictors during the study period. Therefore, to address these limitations and obtain a more comprehensive understanding of E-irAE, further studies with larger sample sizes and standardized monitoring protocols are warranted.

Conclusion

Our retrospective study provided insights into the clinical / biochemical features of E-irAE in Chinese NSCLC patients. Given the increasing use of ICIs, irAEs also become an important part of oncology practice. In our cohort, we observed that TGAb positivity and a low MLR level at baseline were associated with the development of irTD. TSH level at baseline may also predict the onset of irTD, except for Ohe. Additionally, the presence of irTD and male were identified as independent protective factors associated with a longer PFS. Nevertheless, this study still had some limitations. Prospective large-scale studies would be necessary to further evaluate these findings and address the remaining uncertainties surrounding E-irAE in patients with NSCLC.

Data availability

The datasets used during the current study are available from the corresponding author upon reasonable request.

Abbreviations

NSCLC:

Non-small cell lung cancer

irAEs:

Immune-related adverse events

E-irAE:

Immune-related endocrine toxicities

irTD:

Immune-related thyroid dysfunction

ICI:

Immune checkpoint inhibitors

PFS:

Progression free survival

MLR:

Monocyte-to-lymphocyte ratio

PLR:

Platelet-to-lymphocyte ratio

NLR:

Neutrophil-to-lymphocyte ratio

TSH:

Thyroid-stimulating hormone

FT3:

Free Triiodothyronine

FT4:

Free thyroxine

TPOAb:

Thyroid peroxidase antibody

TGAb:

Anti-thyroglobulin antibody

OR:

Odds ratio

CI:

Confidence interval

HR:

Hazard ratio

PD-1/PD-L1:

Programmed death-1/Programmed death ligand-1

OS:

Overall survival

BMI:

Body mass index

ECOG-PS:

Eastern Cooperative Oncology Group-performance status

TFTs:

Thyroid function tests

IQR:

Interquartile range

T1DM:

Type 1 diabetes mellitus

GAD:

Glutamic acid decarboxylase

IA2:

Insulin Antibody

ICA:

Islet cell antibody

ACTH:

Adrenocorticotrophic hormone

FSH:

Follicle-stimulating hormone

LH:

Luteinizing hormon

PRL:

Prolactine

CRP:

C-reactive protein

LDH:

Lactate dehydrogenase

IH:

Hypophysitis

DM:

Immune-related diabetes mellitus

PAI:

Immune-related primary adrenal insufficiency

Oho:

Overt hypothyroidism

Scho:

Subclinical hypothyroidism

Ohe:

Overt hyperthyroidism

Sche:

Subclinical hyperthyroidism

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Zhou C, Wang J, Wang B, Cheng Y, Wang Z, Han B, Lu Y, Wu G, Zhang L, Song Y, et al. [Chinese experts consensus on immune checkpoint inhibitors for Non-small cell lung Cancer (2020 Version)]. Zhongguo Fei Ai Za Zhi. 2021;24(4):217–35.

    PubMed  Google Scholar 

  3. Kwok G, Yau TC, Chiu JW, Tse E, Kwong YL. Pembrolizumab (Keytruda). Hum Vaccin Immunother. 2016;12(11):2777–89.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Herbst RS, Baas P, Kim DW, Felip E, Pérez-Gracia JL, Han JY, Molina J, Kim JH, Arvis CD, Ahn MJ, et al. Pembrolizumab versus docetaxel for previously treated, PD-L1-positive, advanced non-small-cell lung cancer (KEYNOTE-010): a randomised controlled trial. Lancet. 2016;387(10027):1540–50.

    Article  CAS  PubMed  Google Scholar 

  5. Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csőszi T, Fülöp A, Gottfried M, Peled N, Tafreshi A, Cuffe S, et al. Pembrolizumab versus chemotherapy for PD-L1-Positive Non-Small-Cell lung Cancer. N Engl J Med. 2016;375(19):1823–33.

    Article  CAS  PubMed  Google Scholar 

  6. Mok TSK, Wu YL, Kudaba I, Kowalski DM, Cho BC, Turna HZ, Castro G Jr., Srimuninnimit V, Laktionov KK, Bondarenko I, et al. Pembrolizumab versus chemotherapy for previously untreated, PD-L1-expressing, locally advanced or metastatic non-small-cell lung cancer (KEYNOTE-042): a randomised, open-label, controlled, phase 3 trial. Lancet. 2019;393(10183):1819–30.

    Article  CAS  PubMed  Google Scholar 

  7. Langer CJ, Gadgeel SM, Borghaei H, Papadimitrakopoulou VA, Patnaik A, Powell SF, Gentzler RD, Martins RG, Stevenson JP, Jalal SI, et al. Carboplatin and pemetrexed with or without pembrolizumab for advanced, non-squamous non-small-cell lung cancer: a randomised, phase 2 cohort of the open-label KEYNOTE-021 study. Lancet Oncol. 2016;17(11):1497–508.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Gandhi L, Rodríguez-Abreu D, Gadgeel S, Esteban E, Felip E, De Angelis F, Domine M, Clingan P, Hochmair MJ, Powell SF, et al. Pembrolizumab plus chemotherapy in metastatic Non-Small-Cell lung Cancer. N Engl J Med. 2018;378(22):2078–92.

    Article  CAS  PubMed  Google Scholar 

  9. Paz-Ares L, Luft A, Vicente D, Tafreshi A, Gümüş M, Mazières J, Hermes B, Çay Şenler F, Csőszi T, Fülöp A, et al. Pembrolizumab plus chemotherapy for squamous Non-Small-Cell lung Cancer. N Engl J Med. 2018;379(21):2040–51.

    Article  CAS  PubMed  Google Scholar 

  10. Briet C, Albarel F, Kuhn E, Merlen E, Chanson P, Cortet C. Expert opinion on pituitary complications in immunotherapy. Ann Endocrinol (Paris). 2018;79(5):562–8.

    Article  PubMed  Google Scholar 

  11. Stamatouli AM, Quandt Z, Perdigoto AL, Clark PL, Kluger H, Weiss SA, Gettinger S, Sznol M, Young A, Rushakoff R, et al. Collateral damage: Insulin-Dependent diabetes induced with checkpoint inhibitors. Diabetes. 2018;67(8):1471–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Chiloiro S, Bianchi A, Giampietro A, Milardi D, De Marinis L, Pontecorvi A. The changing clinical spectrum of endocrine adverse events in cancer immunotherapy. Trends Endocrinol Metab. 2022;33(2):87–104.

    Article  CAS  PubMed  Google Scholar 

  13. Kimbara S, Fujiwara Y, Iwama S, Ohashi K, Kuchiba A, Arima H, Yamazaki N, Kitano S, Yamamoto N, Ohe Y. Association of antithyroglobulin antibodies with the development of thyroid dysfunction induced by nivolumab. Cancer Sci. 2018;109(11):3583–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Campredon P, Mouly C, Lusque A, Bigay-Game L, Bousquet E, Mazières J, Caron P. Incidence of thyroid dysfunctions during treatment with nivolumab for non-small cell lung cancer: retrospective study of 105 patients. Presse Med. 2019;48(4):e199–207.

    Article  PubMed  Google Scholar 

  15. Pollack R, Ashash A, Cahn A, Rottenberg Y, Stern H, Dresner-Pollak R. Immune checkpoint Inhibitor-induced thyroid dysfunction is associated with higher body mass index. J Clin Endocrinol Metab 2020, 105(10).

  16. Dall’Olio FG, Rizzo A, Mollica V, Massucci M, Maggio I, Massari F. Immortal time bias in the association between toxicity and response for immune checkpoint inhibitors: a meta-analysis. Immunotherapy. 2021;13(3):257–70.

    Article  PubMed  Google Scholar 

  17. Dafni U, Tsourti Z, Vervita K, Peters S. Immune checkpoint inhibitors, alone or in combination with chemotherapy, as first-line treatment for advanced non-small cell lung cancer. A systematic review and network meta-analysis. Lung Cancer. 2019;134:127–40.

    Article  PubMed  Google Scholar 

  18. de Filette J, Andreescu CE, Cools F, Bravenboer B, Velkeniers B. A systematic review and Meta-Analysis of Endocrine-Related adverse events associated with immune checkpoint inhibitors. Horm Metab Res. 2019;51(3):145–56.

    Article  PubMed  Google Scholar 

  19. Barroso-Sousa R, Barry WT, Garrido-Castro AC, Hodi FS, Min L, Krop IE, Tolaney SM. Incidence of endocrine dysfunction following the use of different immune checkpoint inhibitor regimens: A systematic review and Meta-analysis. JAMA Oncol. 2018;4(2):173–82.

    Article  PubMed  Google Scholar 

  20. Boyer M, Şendur MAN, Rodríguez-Abreu D, Park K, Lee DH, Çiçin I, Yumuk PF, Orlandi FJ, Leal TA, Molinier O, et al. Pembrolizumab plus ipilimumab or placebo for metastatic Non-Small-Cell lung Cancer with PD-L1 tumor proportion score ≥ 50%: randomized, Double-Blind phase III KEYNOTE-598 study. J Clin Oncol. 2021;39(21):2327–38.

    Article  CAS  PubMed  Google Scholar 

  21. Baek HS, Jeong C, Shin K, Lee J, Suh H, Lim DJ, Kang MI, Ha J. Association between the type of thyroid dysfunction induced by immune checkpoint inhibitors and prognosis in cancer patients. BMC Endocr Disord. 2022;22(1):89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Yamauchi I, Yasoda A, Matsumoto S, Sakamori Y, Kim YH, Nomura M, Otsuka A, Yamasaki T, Saito R, Kitamura M, et al. Incidence, features, and prognosis of immune-related adverse events involving the thyroid gland induced by nivolumab. PLoS ONE. 2019;14(5):e0216954.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Wu L, Xu Y, Wang X, Cheng X, Zhang Y, Wang Y, Fan X, Zhao H, Liu H, Chai X, et al. Thyroid dysfunction after immune checkpoint inhibitor treatment in a single-center Chinese cohort: a retrospective study. Endocrine. 2023;81(1):123–33.

    Article  CAS  PubMed  Google Scholar 

  24. Muir CA, Menzies AM, Clifton-Bligh R, Tsang VHM. Thyroid toxicity following immune checkpoint inhibitor treatment in advanced Cancer. Thyroid. 2020;30(10):1458–69.

    Article  CAS  PubMed  Google Scholar 

  25. Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, Chow LQ, Vokes EE, Felip E, Holgado E, et al. Nivolumab versus docetaxel in advanced nonsquamous Non-Small-Cell lung Cancer. N Engl J Med. 2015;373(17):1627–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Brahmer J, Reckamp KL, Baas P, Crinò L, Eberhardt WE, Poddubskaya E, Antonia S, Pluzanski A, Vokes EE, Holgado E, et al. Nivolumab versus docetaxel in advanced Squamous-Cell Non-Small-Cell lung Cancer. N Engl J Med. 2015;373(2):123–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Iyer PC, Cabanillas ME, Waguespack SG, Hu MI, Thosani S, Lavis VR, Busaidy NL, Subudhi SK, Diab A, Dadu R. Immune-Related thyroiditis with immune checkpoint inhibitors. Thyroid. 2018;28(10):1243–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lee H, Hodi FS, Giobbie-Hurder A, Ott PA, Buchbinder EI, Haq R, Tolaney S, Barroso-Sousa R, Zhang K, Donahue H, et al. Characterization of thyroid disorders in patients receiving immune checkpoint Inhibition therapy. Cancer Immunol Res. 2017;5(12):1133–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Osorio JC, Ni A, Chaft JE, Pollina R, Kasler MK, Stephens D, Rodriguez C, Cambridge L, Rizvi H, Wolchok JD, et al. Antibody-mediated thyroid dysfunction during T-cell checkpoint Blockade in patients with non-small-cell lung cancer. Ann Oncol. 2017;28(3):583–9.

    Article  CAS  PubMed  Google Scholar 

  30. Yamazaki H, Iwasaki H, Yamashita T, Yoshida T, Suganuma N, Yamanaka T, Masudo K, Nakayama H, Kohagura K, Rino Y, Masuda M. Potential risk factors for Nivolumab-induced thyroid dysfunction. Vivo. 2017;31(6):1225–8.

    CAS  Google Scholar 

  31. Zhou N, Velez MA, Bachrach B, Gukasyan J, Fares CM, Cummings AL, Lind-Lebuffe JP, Akingbemi WO, Li DY, Brodrick PM, et al. Immune checkpoint inhibitor induced thyroid dysfunction is a frequent event post-treatment in NSCLC. Lung Cancer. 2021;161:34–41.

    Article  CAS  PubMed  Google Scholar 

  32. Wu YL, Zhang L, Fan Y, Zhou J, Zhang L, Zhou Q, Li W, Hu C, Chen G, Zhang X, et al. Randomized clinical trial of pembrolizumab vs chemotherapy for previously untreated Chinese patients with PD-L1-positive locally advanced or metastatic non-small-cell lung cancer: KEYNOTE-042 China study. Int J Cancer. 2021;148(9):2313–20.

    Article  CAS  PubMed  Google Scholar 

  33. Muir CA, Clifton-Bligh RJ, Long GV, Scolyer RA, Lo SN, Carlino MS, Tsang VHM, Menzies AM. Thyroid immune-related adverse events following immune checkpoint inhibitor treatment. J Clin Endocrinol Metab. 2021;106(9):e3704–13.

    Article  PubMed  Google Scholar 

  34. Pavan A, Calvetti L, Dal Maso A, Attili I, Del Bianco P, Pasello G, Guarneri V, Aprile G, Conte P, Bonanno L. Peripheral blood markers identify risk of Immune-Related toxicity in advanced Non-Small cell lung Cancer treated with Immune-Checkpoint inhibitors. Oncologist. 2019;24(8):1128–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Peng L, Wang Y, Liu F, Qiu X, Zhang X, Fang C, Qian X, Li Y. Peripheral blood markers predictive of outcome and immune-related adverse events in advanced non-small cell lung cancer treated with PD-1 inhibitors. Cancer Immunol Immunother. 2020;69(9):1813–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Kobayashi T, Iwama S, Yamagami A, Yasuda Y, Okuji T, Ito M, Zhou X, Ando M, Onoue T, Miyata T, et al. Elevated TSH level, TgAb, and prior use of ramucirumab or TKIs as risk factors for thyroid dysfunction in PD-L1 Blockade. J Clin Endocrinol Metab. 2022;107(10):e4115–23.

    Article  PubMed  Google Scholar 

  37. Pollack RM, Kagan M, Lotem M, Dresner-Pollak R. BASELINE TSH LEVEL IS ASSOCIATED WITH RISK OF ANTI-PD-1-INDUCED THYROID DYSFUNCTION. Endocr Pract. 2019;25(8):824–9.

    Article  PubMed  Google Scholar 

  38. Zhou X, Iwama S, Kobayashi T, Ando M, Arima H. Risk of thyroid dysfunction in PD-1 Blockade is stratified by the pattern of TgAb and TPOAb positivity at baseline. J Clin Endocrinol Metab. 2023;108(10):e1056–62.

    Article  PubMed  Google Scholar 

  39. Torimoto K, Okada Y, Nakayamada S, Kubo S, Tanaka Y. Anti-PD-1 antibody therapy induces Hashimoto’s disease with an increase in peripheral blood follicular helper T cells. Thyroid. 2017;27(10):1335–6.

    Article  CAS  PubMed  Google Scholar 

  40. Mazarico I, Capel I, Giménez-Palop O, Albert L, Berges I, Luchtenberg F, García Y, Fernández-Morales LA, De Pedro VJ, Caixàs A, Rigla M. Low frequency of positive antithyroid antibodies is observed in patients with thyroid dysfunction related to immune check point inhibitors. J Endocrinol Invest. 2019;42(12):1443–50.

    Article  CAS  PubMed  Google Scholar 

  41. Yamauchi I, Sakane Y, Fukuda Y, Fujii T, Taura D, Hirata M, Hirota K, Ueda Y, Kanai Y, Yamashita Y, et al. Clinical features of Nivolumab-Induced thyroiditis: A case series study. Thyroid. 2017;27(7):894–901.

    Article  CAS  PubMed  Google Scholar 

  42. Yoon JH, Hong AR, Kim HK, Kang HC. Characteristics of Immune-Related thyroid adverse events in patients treated with PD-1/PD-L1 inhibitors. Endocrinol Metab (Seoul). 2021;36(2):413–23.

    Article  CAS  PubMed  Google Scholar 

  43. Gay S, Rossi G, Corica G, Graziani G, Genova C, Rijavec E, Tagliamento M, Grossi F, Giusti M. RETRACTED ARTICLE: can baseline endocrinological examination and thyroid ultrasound predict the development of thyroid disease in immunotherapy-treated patients? Results from a prospective, single-center, open-label study. Endocrine. 2020;69(1):233.

    Article  CAS  PubMed  Google Scholar 

  44. Kim HI, Kim M, Lee SH, Park SY, Kim YN, Kim H, Jeon MJ, Kim TY, Kim SW, Kim WB, et al. Development of thyroid dysfunction is associated with clinical response to PD-1 Blockade treatment in patients with advanced non-small cell lung cancer. Oncoimmunology. 2017;7(1):e1375642.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Thuillier P, Joly C, Alavi Z, Crouzeix G, Descourt R, Quere G, Kerlan V, Roudaut N. Thyroid dysfunction induced by immune checkpoint inhibitors is associated with a better progression-free survival and overall survival in non-small cell lung cancer: an original cohort study. Cancer Immunol Immunother. 2021;70(7):2023–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all patients, nurses, and physicians who participated in the study.

Funding

We thank all participating study sites for collaboration on this project. This research was supported by Zhejiang Provincial Clinical Research Center for Respiratory Disease (2022E50005 to Jianya Zhou) and Key R&D Program of Zhejiang Respiratory Disease (No.2023C03069).

Author information

Authors and Affiliations

Authors

Contributions

H.Z., J.Z., C.R., C.Y., X.W., X.L. and Y.L. provided the raw data. H.Z. and J.Z. conducted the statistical analysis. The first draft of the manuscript was written by H.Z., J.Z. and H.Z. carefully edited and guided writing of the manuscript. J.Z. and J.Z. are co-corresponding authors. All authors contributed to the article and approved the final manuscript.

Corresponding authors

Correspondence to Jianya Zhou or Jianying Zhou.

Ethics declarations

Ethical approval

All experiments were performed in accordance with relevant guidelines and regulations (such as the Declaration of Helsinki). This study was reviewed and approved by the Ethical Committee of the First Affiliated Hospital, College of Medicine, Zhejiang University on October 24th (Number: 2019 − 1433) and by the Ethical Committee of the Fourth Affiliated Hospital, College of Medicine, Zhejiang University on November 13th (Number: 20231108). Given the retrospective nature of this study, informed consent exemption was obtained for this study from Ethics Committee of the First Affiliated Hospital of Zhejiang University, the Fourth Affiliated Hospital of Zhejiang University, Guangfu Hospital, Affiliated Hospital of Jiaxing University and the First Affiliated Hospital Wenzhou Medical College.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, H., Zheng, J., Ren, C. et al. Risk factors of immune-related endocrine toxicities in non-small cell lung cancer patients treated with pembrolizumab and its impact on patient outcomes: a multicenter retrospective study. BMC Pulm Med 25, 111 (2025). https://doi.org/10.1186/s12890-025-03570-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-025-03570-8

Keywords