Descriptive Epidemiology and Survival Analysis of Prolactinomas and Cushing's Disease in Korea.

BACKGROUND: Only a few studies have established the epidemiology of prolactinoma and Cushing's disease in Korea. Furthermore, the incidence of these disease are increasing than before associated with the development of technologies. This study was designed to evaluate the epidemiology of prolactinoma and Cushing's disease and their survival analysis according to treatment.
METHODS: The nationwide, population-based study evaluated incidence and prevalence of prolactinoma and Cushing's disease using de-identified claims data in The Korean Health Insurance Review and Assessment Service database between 2013 and 2017. The survival analysis investigated regarding treatment over a period of 6 years. A log-rank test and Cox proportional hazard regression analysis were used.
RESULTS: The 6,056 patients with newly diagnosed prolactinoma and 584 patients with Cushing's disease were recorded between 2013 and 2017. The annual incidence of prolactinoma was 23.5 cases per million, and its prevalence was 82.5 cases per million, and 2.3 cases per million/year and 9.8 cases per million for Cushing's disease. The survival benefit was insignificant in prolactinoma according to treatment, but treatment of Cushing's disease ameliorated the survival rate significantly.
CONCLUSION: Overall, the incidence of prolactinoma and Cushing's disease was similar with those found previously, but the prevalence of two diseases were inconsistent when compared with the early studies. The present study also proposed necessity of treatment in Cushing's disease for improving the survival rate.

INTRODUCTION

A pituitary adenoma (PA) has been considered to be a rare disease in the past; however, it is now reported more frequently than before owing to the development of technologies associated with radiological and hormonal diagnostic techniques [1-3]. The socioeconomic burden of these diseases is also increasing with the high morbidity and mortality rates [3]. In Korea, previously reported data have revealed a deterioration of the burden of these diseases on health care [4]. Although the epidemiologic study is one of the first tasks for understanding a disease, the studies on PAs in Korea were insufficient, and most of the previous studies analyzed Western populations. In addition, the earlier studies have shown inconsistent results because most of these were based on small-population, cross-sectional design studies, and their data sources were a single surgical center or cancer registry [1-3,5-9].

PAs are classified as macroadenomas (>10 mm in diameter) or microadenomas (≤10 mm in diameter) depending on their diameter [10,11]. PAs are also divided into functional or nonfunctional adenomas depending on their hormonal hyper-secretory capacity [12]. Functional PAs are categorized into several subtypes according to the hormone secreted by each tumor. Each subtype has different morbidity and mortality rates. Unlike a number of epidemiological studies regarding acromegaly [1,6- 9,13-22], there were few studies about morbidity and mortality associated with prolactinomas and Cushing’s disease (CD). A prolactinoma is known to be one of the most frequent subtypes of PA, which accounts for 40% to 50% of PAs. CD is rare, accounting for 4% to 6% of PAs and has a high cardiovascular mortality rate [2,23-26]. Despite their importance in public health care, epidemiologic data on morbidity and mortality of these two diseases is not enough, even with incoherent results. Furthermore, some of non-functional PAs, especially macroadenomas, may be capable of secreting prolactin [27-29]; which makes it difficult to determine the accurate epidemiology of pure prolactinomas. Likewise, a low incidence of CD challenges investigators to differentiate from non-pituitary Cushing’s syndrome (CS); therefore, a well-equipped diagnostic laboratory system and protocol with experienced endocrinologist are mandatory to guarantee diagnostic accuracy leading to the proper analysis of epidemiologic data for CD.

Furthermore, since epidemiologic studies depend on population-specific registries with regional differences in practices, such as diagnostic practices, reporting patterns, and case definitions [30], there are limitations to accepting the results of previous studies based on different populations. Therefore, it was necessary to conduct an analysis using Korean population-based data. The Korean Health Insurance Review and Assessment Service (HIRA), which was introduced in 1989, collects claims data for the purpose of national insurance coverage in Korea. This data includes health care utilization, health screening, and mortality for the entire population of Korea that includes over 50 million people based on the HIRA database, we conducted a nationwide, population-based study evaluated the epidemiology of prolactinoma and CD, along with their recent trends, in Korea,

METHODS

Data collection

The HIRA database includes insurance claims data of the entire Korean population that is covered by the national healthcare insurance system and is open to the public, with the aim of contributing towards establishing knowledge related to health care. As the database is primarily for insurance claims, the establishment of appropriate operational definitions using assessment codes such as laboratory test, procedure, and drug prescription, in addition to the diagnostic code is inevitable for the reliability of our analysis. Therefore, in the present study, we extracted deidentified claims data containing the treatment history of all patients with prolactinomas or CD who visited a healthcare facility. To identify suitable subjects for analysis, we established an operational definition using the International Classification of Disease 10th revision (ICD-10) codes. In order to assure more accurate operational definitions in our study, we cross matched the benefit extension policy (BEP) codes with the HIRA codes. The BEP was established to support medical expenses for patients with rare and incurable diseases such as cancers, rare genetic disorders, and so on [31]. The BEP code is mandatory to exempt certain amount of medical expenses for managing registered diseases. Double comparison with the HIRA dataset and the BEP codes is believed to improve the reliability of data extracted in our study population.

In this study, patients having diagnostic codes (ICD-10) of hyperprolactinemia (E22.1) and benign neoplasm of the pituitary gland (D35.2) with the corresponding BEP codes (V162 for PAs and V112 for hyperprolactinemia) for prolactinomas were identified as patients with prolactinoma. Patients with CD were defined as those who were registered with the diagnostic codes of pituitary-dependent CD (E24.0) and benign neoplasm of the pituitary gland (D35.2) along with the corresponding BEP code for PAs (V162) and CD (V114). In addition, patients registered with the code of other CS or CS unspecified (E24.8) were also included in the study only if they had a code for both CD (E24.0) and benign neoplasm of the pituitary gland (D35.2) at the same time.

The present study included all patients with prolactinoma and CD who were registered in the HIRA database and had a history of receiving outpatient care or hospitalization between January 2009 and December 2017. The annual incidence and prevalence rates were analyzed between January 2013 and December 2017. For the annual incidence rate, the period from 2009 to 2012 was set as the washout period. We used the resident registration demographics established by Statistics Korea as the denominator.

Patients with newly diagnosed prolactinoma or CD between January 2011 and December 2012 were extracted, and the survival rates of these patients were estimated for the period from January 2013 to December 2017. Since one of the previous studies which investigated the epidemiology of PAs in Korea employed the analysis period between 2009 and 2013 [4], the longest period we can manage was 5 years to avoid overlapping the study period. In addition, in order to establish the secure time period for the adequate survival analysis, we were obliged to investigate the group of patients newly diagnosed between 2011 and 2012 with wash out period between 2009 and 2010, different from those used for the analysis of the aforementioned epidemiology. The differences in survival probability according to treatment were compared. The treated group was defined as patients with any record of treatment including the first-generation somatostatin analogues such as octreotide or lanreotide, surgery such as trans-sphenoidal approaches, and radiation therapy including gamma knife surgery after their diagnostic codes were entered, while the untreated group included those without any corresponding treatment record but only with diagnostic codes.

The present study was established in accordance with the guidelines of the Declaration of Helsinki and approved by the Kyung Hee University Hospital Institutional Review Board (KHUH 2018-11-001). Informed consent from the subjects was waived due to the retrospective design of this study.

Statistical analysis

Hazard ratios (HRs) and 95% confidence intervals (CIs) for mortality were evaluated using a log-rank test and Cox proportional hazard regression analysis. Multivariate model analysis was performed after adjusting for age, sex, and income level. The survival period was defined as the time between the first or initial day of observation and the day when a new clinical event occurred or the end of the observation period on December 31, 2017. Patients treated before the initial diagnosis of prolactinomas or CD were excluded from the analysis. The results were considered statistically significant if the P value was <0.05. Statistical analyses were performed using SAS software version 9.4 (SAS Institute, Cary, NC, USA).

RESULTS

Epidemiologic data

Prolactinoma

Between January 2013 and December 2017, 6,056 patients were newly diagnosed with prolactinoma, with approximately 1,000 patients diagnosed annually (Table 1, Fig. 1A). The annual incidence was 23.5 cases per million, and the prevalence was 82.5 cases per million. In 2016, the incidence was 26.0 cases per million, and the prevalence tended to increase during this period; in 2017, 94.4 cases per million were observed (Fig. 1B). Of the newly diagnosed patients, 91.9% were women. The highest incidence was observed in individuals in their 40s, accounting for 36.3% of the total patients, while elderly patients aged 80 years or more demonstrated an incidence rate of <1%, and adolescents accounted for approximately 4.3% of total cases (Table 1).

Table 1.

The Number of Cases with Prolactinoma in Korea According to Sex and Year

	2013	2014	2015	2016	2017
Sex
Male	86 (8)	101 (9.1)	105 (8.8)	89 (6.6)	98 (7.4)
Female	990 (91.9)	1,005 (90.8)	1,092 (91.1)	1,255 (93.1)	1,227 (92.6)
Unknown	1	1	2	4
Age
0–19	44 (4.1)	41 (3.7)	48 (4)	72 (5.3)	64 (4.8)
20–29	213 (19.8)	247 (22.3)	239 (19.9)	276 (20.5)	307 (23.2)
30–39	409 (38)	403 (36.4)	428 (35.7)	492 (36.5)	464 (35)
40–49	277 (25.7)	263 (23.8)	335 (27.9)	335 (24.9)	345 (26)
50–59	89 (8.3)	107 (9.7)	101 (8.4)	115 (8.5)	99 (7.5)
60–69	27 (2.5)	24 (2.2)	30 (2.5)	34 (2.5)	30 (2.3)
70–79	15 (1.4)	18 (1.6)	15 (1.3)	16 (1.2)	12 (0.9)
>79	2	3 (0.3)	1 (0.1)	4 (0.3)	4 (0.3)
Unknown	1	1	2	4
Total	1,077	1,107	1,199	1,348	1,325

Values are expressed as number (%).

Fig. 1.

Incidence rate of prolactinoma in Korea during 2013 to 2017: (A) number of newly diagnosed patients with prolactinoma (B) incidence and prevalence rates of prolactinoma per million.

Cushing’s disease

Between January 2013 and December 2017, a total of 584 newly diagnosed CD patients were recorded, with more than 100 cases each year (Table 2, Fig. 2A). Majority of patients were females, accounting for 88.0% of all cases. The number of patients in each age group were evenly distributed between 20 and 70 years, accounting for approximately 10% to 20% of the total cases (Table 2). The annual incidence of CD was 2.3 cases per million, and its prevalence was 9.3 cases per million (Fig. 2B).

Table 2.

The Number of Cases with Cushing’s Disease in Korea According to Sex and Year

	2013	2014	2015	2016	2017
Sex
Male	14 (12.5)	19 (12.9)	23 (22.8)	22 (20)	25 (21.9)
Female	98 (87.5)	128 (87.1)	78 (77.2)	88 (80)	89 (78.1)
Unknown
Age
0–19	4 (3.6)	8 (5.4)	5 (5)	5 (4.5)	2 (1.8)
20–29	17 (15.2)	14 (9.5)	10 (9.9)	14 (12.7)	14 (12.3)
30–39	30 (26.8)	26 (17.7)	18 (17.8)	21 (19.1)	24 (21.1)
40–49	18 (16.1)	34 (23.1)	23 (22.8)	20 (18.2)	19 (16.7)
50–59	19 (17)	26 (17.7)	23 (22.8)	23 (20.9)	30 (26.3)
60–69	13 (11.6)	24 (16.3)	15 (14.9)	11 (10)	15 (13.2)
70–79	8 (7.1)	13 (8.8)	7 (6.9)	13 (11.8)	9 (7.9)
>79	3 (2.7)	2 (1.4)		3 (2.7)	1 (0.9)
Unknown
Total	112	147	101	110	114

Values are expressed as number (%).

Fig. 2.

Incidence rate of Cushing’s disease in Korea during 2013 to 2017: (A) number of newly diagnosed patients with Cushing’s disease (B) incidence and prevalence rates of Cushing’s disease per million.

Survival analysis

From January 2013 to December 2017, there were a total of 35/1,685 cases with death were noted in patients with prolactinoma, and 20/1,172 cases were shown to undergo either medical or surgical treatment, and 15/511 cases appeared to be untreated. This led to a slightly lower mortality rate in the treated group than those in the untreated group, but this difference was not statistically significant (HR, 1.471; 95% CI, 0.745 to 2.902; P=0.1066) (Fig. 3A). On the other hands, there were a total of 14/117 cases with death were noted in patients with CD, and 9/102 cases were shown to undergo either medical or surgical treatment, and 5/15 cases were found to be untreated. This resulted in a significantly lower mortality rate in the treated group than those in the untreated group (HR, 4.648; 95% CI, 1.391 to 8.183; P=0.0029) (Fig. 3B).

Fig. 3.

Probability of survival in (A) prolactinoma (B) Cushing’s disease between treated and untreated patients (curve difference, P=0.1066 for prolactinoma; P=0.0029 for Cushing’s disease).

DISCUSSION

This study was designed to investigate the epidemiology of prolactinomas and CD in Korea using a nationwide, population-based study based on the HIRA database, which identified 6,056 cases of prolactinomas and 584 cases of CD between 2013 and 2017. The incidence of prolactinomas was 23.5 cases per million, and its prevalence was 82.5 cases per million. The incidence of CD was 2.3 cases per million, and its prevalence was 9.8 cases per million. The survival rate of patients with prolactinoma who were treated was slightly higher than that of untreated patients, but the difference was insignificant. However, the survival rate of treated patients with CD was significantly higher than that of untreated patients.

There are a few population-based studies related to the epidemiology of prolactinomas. Gruppetta et al. [7] and Agustsson et al. [13], reported that the incidence of prolactinoma was 20 cases per million in the population of Malta and Iceland, respectively, which is similar to that in the present study. However, the prevalence in these studies were considerably different from each other, with Gruppetta et al. [7] reporting 350 cases per million and Agustsson et al. [13] reporting 543.7 cases per million. Moreover, a community-based study in Banbury in the UK reported a prevalence of 444 cases per million for prolactinoma [6]. The reason for this inconsistency was possibly because there were differences in the population size targeted by each study, data source, and methods of analysis. In Korea, two studies have been reported which investigated the epidemiology of PAs using the HIRA database. Park et al. [4] reported that the incidence of hyperprolactinemia was 15 cases per million, with the prevalence increasing from 38.5 to 68.6 cases per million between 2009 and 2013, while Oh et al. [25] reported only the incidence, which was 16 cases per million. The disagreement between studies using the same database with respect to the incidence is possibly due to differences in the operational definition of the disease in terms of different combination of diagnostic and BEP codes in each study [4,25]. In addition, previous results in Korea demonstrated that the highest incidence rate was observed in patients in their 30s, and a sharp decline of the incidence rate was observed among older individuals [4]. Similarly, in our study, the patients in their 30s accounted for the largest proportion among the entire patients (36.3%); the high proportion was maintained until the 40s but decreased thereafter. The incidence of females were a 10-fold higher than that of males, but this sex disparity tended to become equalized in patients aged >50 years [1,4,6]. One hypothesis for this phenomenon is that the clinical symptoms and signs related to amenorrhea-galactorrhea syndrome are more definite in women than in men [25].

Since CD is rare and difficult to diagnose, most epidemiologic studies have focused on the entire spectrum of endogenous CS. Gruppetta et al. [7] reported that in Malta the incidence of CD was 0.7 cases per million, and the prevalence was 16.8 cases per million. Agustsson et al. [13] reported that the incidence and prevalence in Iceland were five cases per million for women and 62 cases per million, respectively. The community-based study at Banbury in the UK reported only a prevalence of 12 cases per million [6]. Another community-based study in the residents of Vizcaya in Spain reported that the incidence of CD was 2.4 cases per million, but the prevalence was significantly different at 39.1 cases per million [2]. The incidence and prevalence of CD have been mostly inconsistent according to the studies varied between the studies, possibly due to the characteristically low incidence of CD with different study designs and small size of populations in previous studies. In addition, it should be emphasized for the difficulty in differentiating between CD and CS. In the case of some studies conducted before the 90s or the early 2000s, magnetic resonance imaging was not used to diagnose PAs, and laboratory techniques may not have been sufficiently accurate. Therefore, there are many limitations to the comparison between the results from these studies and the results of our study based on recent data. Female predominance in patients with CD was observed as seen in previous studies [2,9,13,24]. In addition, newly diagnosed patients were evenly distributed between 20 and 60 years of age, indicating the same trend as in previous studies [2].

Despite the low mortality rates observed among patients with prolactinoma who were treated, there was no significant difference in mortality rates when compared with that of untreated patients. Although the most common cause of death in patients with PA is cardiovascular disease (CVD) and its associated conditions [3], the incidence of CVD in patients with prolactinoma has been reported to be lower than those with adrenocorticotropic hormone-, thyroid stimulation hormone-secreting adenomas, or non-functional adenomas [25]. In addition, the mortality rate of patients with prolactinoma was similar to that of the general population [25], which may lead to insignificant difference in mortality regardless of treatment. On the other hand, in patients with CD, significant improvement was observed in the mortality rate of patients in the treated group. CD is known to have a high CVD-associated mortality rate because of impaired glucose metabolism and persistent hypertension caused by neurohormonal effects [2,23-26]. Few studies have demonstrated the long-term mortality in CD, with conflicting results. Colao et al. [23] documented that CVD risk factors such as body mass index, blood pressure, and fasting blood glucose level were increased in patients with CD even after 5 years of successful treatment when compared with that age- and sex-matched controls, which presumably played a role to rather aggravate the CVD risks. However, this was a small-scale study of 15 patients with CD who were cured. A previous population-based study was consistent with the results of our study, which suggested long-term survival after successful transsphenoidal surgery for CD was not markedly different compared to the background population [24], and another case series study was also in close agreement with ours [26].

Our study has certain limitations. First, the HIRA database may not be optimal for estimating the incidence and prevalence of diseases, since the data have been collected for tracking medical service claims and reimbursement [32]. Moreover, the accuracy of the diagnostic code entered by each physician cannot be guaranteed, possibly leading to over- or under-estimation. To improve the reliability and accuracy of the extracted data, an operational definition was established using the simultaneous application of BEP and ICD-10 codes. It has been previously reported that the diagnostic accuracy of operational definitions was improved by combining the HIRA database and BEP codes [25], which showed that re-evaluation of medical records of patients who were extracted by operational definition resulted in the accuracy of diagnostic codes had 93% sensitivity and 94% specificity [25]. Unfortunately, we were not able to apply the same method to evaluate the accuracy of our operational definition as done in the previous study [25], because the raw dataset initially provided by the HIRA was anonymous after de-identification and was allowed to access in the strictly authorized area only, not eligible for comparing with data in our hospital or finding out the detailed individual information of each subjects. Second, our survival analysis was performed after adjusting for age, sex and income levels, but for more accurate analysis it was necessary to perform the survival analysis after adjustment with additional covariates and cause of death, which may possibly undermine the representativeness of the actual clinical course between treated and untreated group. Finally, in patients with CD, a significant improvement in survival was observed in the treated group compared to that in the untreated group; however, a comparative analysis with the general population was not performed. Despite using the study design similar with previous two studies in Korea based on same database, however, it should be of note that our study attempted to investigate the most recent trend of the epidemiology of prolactinoma and CD. In addition, it is meaningful to mention that our study compared the mortality between the treated and untreated group, but not general population to observe the effect of treatment more directly.

In conclusion, the annual incidence of prolactinoma was similar to that reported in previous studies; however, its prevalence was considerably different. In patients with CD, the annual incidence was also similar to that in previous studies although its prevalence was lower than that in other countries. In addition, a significant difference in survival rates was observed between treated and untreated patients with CD, which sheds light on the role of the treatment in improving 281 the survival rate in CD. We are currently planning to conduct a comparative analysis of the mortality rate with the general population with long term follow-up. The epidemiology of other types of PA, including acromegaly, thyrotropin-producing PA, is also being evaluated.