Multimorbidity patterns in the working age population with the top 10% medical cost from exhaustive insurance claims data of Japan Health Insurance Association

Abstract

Although the economic burden of multimorbidity is a growing global challenge, the contribution of multimorbidity in patients with high medical expenses remains unclear. We aimed to clarify multimorbidity patterns that have a large impact on medical costs in the Japanese population. We conducted a cross-sectional study using health insurance claims data provided by the Japan Health Insurance Association. Latent class analysis (LCA) was used to identify multimorbidity patterns in 1,698,902 patients who had the top 10% of total medical costs in 2015. The present parameters of the LCA model included 68 disease labels that were frequent among this population. Moreover, subgroup analysis was performed using a generalized linear model (GLM) to assess the factors influencing annual medical cost and 5-year mortality. As a result of obtaining 30 latent classes, the kidney disease class required the most expensive cost per capita, while the highest portion (28.6%) of the total medical cost was spent on metabolic syndrome (MetS) classes, which were characterized by hypertension, dyslipidemia, and type 2 diabetes. GLM applied to patients with MetS classes showed that cardiovascular diseases or complex conditions, including malignancies, were powerful determinants of medical cost and mortality. MetS was classified into 7 classes based on real-world data and accounts for a large portion of the total medical costs. MetS classes with cardiovascular diseases or complex conditions, including malignancies, have a significant impact on medical costs and mortality.

Introduction

Japan implements a universal medical care insurance system in which all citizens subscribe to healthcare insurance systems. People in the working-age population (under the age of 70) receive treatment at 30% of the total direct medical costs. The remaining 70% is covered by the employee’s insurance premium, of which the employee and employer each pay 50% [1]. Furthermore, medical expenses incurred during treatment are set at a uniform price nationwide, and patients can freely access any health care provider. This greatly contributes to the stabilization of people’s lives and maintenance of their health. However, it is difficult to maintain this system because of the imbalance between the increase in medical expenditures and the low economic growth caused by an aged society with fewer children. Prolonging this situation will force insurance fees to increase, placing a financial burden on subscribers and companies. In the worst-case scenario, the medical insurance system might go bankrupt because of an increase in dismissals and unemployment. Therefore, to sustain the medical insurance system, the growth of high medical expenses must be curtailed.

Recently, many developed countries have entered an aging society and have attempted to control health care costs for elderly people. Japan is one of the countries with the highest aging population, and there are several studies on health care costs focusing on the elderly population [2, 3]. Total medical expenses among the working-age population (younger than 65 years of age) account for approximately 40%, and these costs are expected to increase due to an increase in noncommunicable diseases [4]; thus, health care cost containment among the working-age population is also an important challenge to sustain the medical insurance system.

A small proportion of patients in the United States account for the majority of total medical expenditures [5]. Similarly, the top 10% of patients account for 60% of total annual health care costs in the working-age population in Japan [6]. If patients with mild symptoms are prevented from moving to the high-cost population, it might be possible to curb future increases in medical costs. The high-cost population is expected to include not only patients with a rare disease requiring an expensive drug or surgery due to an unexpected accident, but also those who have multiple chronic conditions, in other words, “multimorbidity” [7].

Multimorbidity is defined as two or more coexisting chronic conditions in an individual [8], and it is a growing global challenge with substantial effects on individuals, caregivers, and society. This condition affects health outcomes (e.g., increased mortality or deterioration of physical function and quality of life), patient burden (e.g., increase in the number of consultations or risk of polypharmacy), and health economic burden (e.g., increased medical expenses due to emergency hospitalization and increased consultation opportunities) [9, 10]. Although a single disease-oriented model is not suitable for the treatment of multimorbidity [11], patients in Japan are often prone to fragmented treatment. To develop a new treatment model for multimorbidity, it is important to clarify multimorbidity patterns in Japanese individuals. A recent study emphasized the need for primary prevention in the younger population because the risk of mortality was higher in younger and middle-aged adult populations with multimorbidity than in older populations [12]. Therefore, clarifying multimorbidity patterns among the working-age population with high medical costs would help health managers and policymakers to curb the increase in future medical expenses and patients with critical diseases.

The present study aimed to clarify the contribution of multimorbidity among the working-age population with the highest medical cost using latent class analysis (LCA) and the patterns that have a large impact on medical expenses. We also analyzed exhaustive insurance claims data from the Japan Health Insurance Association and provided a greater understanding of multimorbidity in high-cost patients.

Methods

Ethics approval

Anonymized data were used in this study. This study (S1 Checklist) was approved by the Ethics Committees of Keio University (approval number [no.]: 2020–05) and Tokyo Medical and Dental University (approval no.: M2020-385) and was exempt from the need to obtain informed consent from participants.

Study design and data source

We conducted a cross-sectional study using exhaustive health insurance claims data provided by the Japan Health Insurance Association, which comprises the largest working-age population in Japan. This database consists of approximately 40 million people younger than 75 years of age who are formally or informally employed at small to medium enterprises and their family members. The claim data include information on total direct medical costs including patients’ co-payment comprised from inpatient services, outpatient services, and outpatient prescription drugs. However, socioeconomic status (SES), such as education level and income, are not included. Currently, only academic researchers selected by open recruitment can use data for the 7 years from 2015 to 2021.

Diagnosed disease prevalence was represented by the International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD10), which was documented between April 2015 and March 2016. The annual medical cost for each patient was calculated by summing the outpatient, hospital, drug, and dental costs during the same period (JPY was translated to USD at the 2015 rate of 1 USD = 120 JPY).

Study population

The study population comprised insurance affiliates aged 18–64 years as of April 1, 2015 (n = 20,803,665). The exclusion criteria were as follows: insurance affiliates younger than 18 years of age and older than 65 years of age as of April 1, 2015 (n = 1,454,866) and those who took out insurance or withdrew in the middle of the 2015 fiscal year (n = 2,359,770).

Analyses

LCA was used to identify multimorbidity patterns in the population with high medical costs. LCA is a model-based approach to classification that permits a statistical evaluation of how well the proposed LCA model represents the data [13]. The present parameters of the LCA model were disease labels based on the ICD10 diagnosis code. ICD10 diagnosis codes were indicated with the first 3 characters, and the following codes were excluded: (1) codes, the frequency of which in the present population was less than 1%, which were excluded to obtain consistent and clinically interpretable patterns of association; (2) transient disease codes such as common cold; (3) symptom codes (R00-R99); and (4) non-disease codes (S00-U99). As a result, 252 diagnostic codes were extracted (S1 Table). Furthermore, these diagnosis codes were classified into 68 disease labels based on pathological similarity by an expert panel (Table 1), and the labels were used in LCA as a binary variable. Labels assigned to the identified latent disease classes were determined with reference to the highest item-response probability in the class, higher probability than the other classes, treatment content, and sex proportion. The item-response probabilities of the disease labels are shown in a heatmap for each class for visual clarity. The prevalence rates of some disease labels (e.g., oral disease, eye disease, chronic inflammation of the upper and lower airways, and skin and subcutaneous tissue) were high in most classes because the present study focused on patients with high medical costs, and the prevalence rates of these diseases in the present study were higher than those in the general population, including healthy individuals. Therefore, if the probability of the class was below the prevalence rate in all individuals, the box in the heatmap was represented as colorless.

[Figure omitted. See PDF.]

Multimorbidity is formally defined as 2 or more chronic conditions [14]; thus, we selected 46 chronic conditions from 68 disease labels in Table 2 and estimated the multimorbidity (S2 Table). In chronic diseases, malignancy was grouped into one category, and the following diseases not important in terms of multimorbidity were excluded: oral disease, eye disease, and skin and subcutaneous tissue disease.

[Figure omitted. See PDF.]

In addition, subgroup analysis was performed and focused on metabolic syndrome (MetS) classes, which were characterized by the MetS triad (type 2 diabetes, dyslipidemia, and hypertensive disease) because these medical costs accounted for the largest share of the total. The generalized linear model (GLM) with a log-link function and gamma distribution was generated to assess the factors influencing the total cost and marginal effects of the total medical cost [15]. Multiple logistic regression analysis was used to estimate the odds ratios for the associations with the 5-year mortality rate. A 95% confidence level was used to determine statistically significant variables for the total cost and mortality.

Representative values are shown as medians [quartile Q1, Q3]. LCA was performed using JMP Pro 16 (SAS Institute Japan Co., Ltd., Cary, NC, USA). R 4.3.0 statistical software was used for the other analyses, and average marginal effect of cost was calculated using the margins package.

Results

Demographic characteristics of participants

After excluding the inapplicable persons, we extracted the top 10% of patients with high medical costs and finally analyzed the data of 1,698,902 patients (671,450 women and 1,027,452 men). The distribution of annual medical costs in all insured populations is shown in S1 Fig. The demographic characteristics of the study population are shown in Table 2. The median number of co-existing chronic conditions in an individual, out of the selected list of 46 chronic conditions, was 7 [5, 10]. Overall, 95.6% and 91.1% of all patients had 2 or more conditions and 3 or more conditions, respectively; thus, most of them were in a multimorbidity state. The total medical cost was 1,397,311,088,550 JPY (11,644,259,071 USD) and that per capita was 822,480 JPY (6,854 USD). Of the total costs, the proportions of inpatient cost and drug cost were 39.4% and 21.0%, respectively. Men incurred higher medical cost per capita [879,930 JPY (7,333 USD)] than women [734,570 JPY (6,121 USD)], whereas the number of outpatient visits for women seemed to be higher than that for men. The median length of hospital stay was 0 [0, 5], and 31.6% of patients were hospitalized for more than 3 days. Five-year retirement and mortality rates were 46.1% and 1.7%, respectively. The comparisons of demographic characteristics between the top 10% of patients with high medical costs and patients with low to medium costs (i.e., lower 90%) or 0 costs are shown in S3 and S4 Tables.

Determining the number of latent classes

In the preliminary analysis, we started with a 2-class model and successively increased the number of latent models to 40 to select the optimal number of latent classes. Model fitness was examined using the Akaike information criterion (AIC) and Bayesian information criterion (BIC). However, the AIC and BIC did not reach a minimum until 40 classes and seemed to continue to decline thereafter (S2 Fig). Too many classes would make the interpretation difficult; thus, we examined 5 models of 20, 25, 30, 35, and 40 classes, and assigned labels to the identified latent disease classes in each model (S3 Fig). Common classes that appeared in all models were as follows: single chronic condition, chronic inflammation of airways, mental disease, digestive disease, motor disorder, diseases specific to women or men, MetS, connective tissue disease, malignancy, and kidney disease. Classes of urologic disease did not appear in the 20-class model, whereas some classes were difficult to interpret in the 35- and 40- class models. Moreover, malignancy classes were more clearly identified in the 30-class model than in the 25-class model. Consequently, a 30-class model was selected as the final model.

Characteristics of the 30 latent class models

Table 3 summarizes the label names and the features of each class, and Table 4 shows characteristics of the 30 latent class models as follows: single chronic condition (6.9%); otorhinolaryngology diseases (5.8%); 3 types of mental disease (6.7%); liver disease (3.3%); 2 types of digestive diseases (7.9%); 2 types of motor disorders (7.8%); urologic disease (2.0%); 4 types of diseases specific to women (11.9%); diseases specific to men (2.6%); 7 types of MetSs (31.8%); connective tissue disease (2.2%); 5 types of malignancies (10.0%); and kidney disease (1.2%). A heatmap of all classes is shown in S4 Fig. Class 1 was named as the “single chronic condition.” Patients in this class frequently underwent surgery for any fractures and had only one chronic condition (1 [1, 2]); therefore, they were considered as patients with non-multimorbidity. Mental disease classes included patients with a high prevalence of mood disorder, neurotic disorder, and sleep disorder simultaneously and suggested a higher 5-year retirement rate (51.4–54.7%) than that of other classes. If more than 95% of patients within the class were women or men, that class was named as a sex-specific disease class. One of these sex-specific classes included younger women (33 [29, 37] years) who received perinatal care, and the other classes consisted of patients with multimorbidity. There were 7 MetS classes that showed high frequencies of type 2 diabetes, dyslipidemia, and hypertension simultaneously. Malignancy classes were classified into breast cancer, lung cancer, cirrhosis/hepatoma, gastrointestinal cancer including stomach cancer or colorectal cancer, and cancer except for genital cancers. Although the breast cancer class consisted mostly of females (98.3%), the prevalence of female genital tract cancer was not substantial and even males can develop breast cancer infrequently; thus, the breast cancer class was separated from the sex-specific disease class. Patients in these malignancy classes showed a higher 5-year mortality rate than those in the other classes. The kidney disease class had the smallest number of patients in all classes and showed that the prevalence rates of chronic kidney disease (CKD) and hypertension were 90% and 95%, respectively.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Latent class ratios by age and sex

Fig 1 shows the percentage of each patient per class as a 100% stacked vertical bar chart by age group and sex. Among men, the single chronic condition class included a large proportion of patients aged 18–29 years, and the percentage of patients in the MetS classes increased with age and reached more than 50% after age 50–59 years. Among women, patients with diseases specific to their sex accounted for nearly half for those aged 18–29 and 30–39 years. The percentages of breast cancer, MetS, and motor disorder classes increased from age 40s to 60s. As a feature common to both sexes, mental disease classes were larger in the young population than in the old population.

[Figure omitted. See PDF.]

The total medical cost and cost per capita by latent classes

The total cost and cost per person for each class are shown in Fig 2. Although total cost tended to be higher in the class with a large number of patients, some classes of MetS, i.e., the malignancy class and kidney disease class, showed higher total cost due to a high cost per capita. In particular, the kidney disease class showed the highest total cost and cost per person among all classes, despite the small number of patients.

[Figure omitted. See PDF.]

Vertical bars represent the total cost, and the yellow line represents the cost per capita. 1 USD = 120 JPY.

Subgroup analysis of the seven MetS classes

Herein, 31.8% of patients were classified into some type of MetS class, and their combined medical costs accounted for 28.6% of the total. To understand the characteristics of the MetS class in more detail, a heatmap focusing on the 7 MetS classes is shown in Fig 3, in which only highly frequent diseases in the MetS class are shown. The notable features of each MetS class were as follows: MetS class 1, only MetS triad (type 2 diabetes, dyslipidemia, and hypertensive disease); MetS class 2, triad plus infection, chronic inflammation of the upper and lower airways, oral disease, gastrointestinal tract disease, and other locomotive diseases; MetS class 3, triad plus motor disorder; MetS class 4, triad plus liver disease and kidney disease; MetS class 5, triad plus cardiovascular diseases; MetS class 6, triad plus other metabolic diseases, cardiovascular diseases, gastrointestinal disease, and other locomotive diseases; and MetS class 7, triad plus complex conditions including malignancies.

[Figure omitted. See PDF.]

Thick red represents a high probability.

Table 5 shows the results of the GLM used to determine the variables affecting the annual total medical cost and 5-year mortality in the MetS classes. Total medical costs in MetS classes 5 and 6 were more than 20% higher than that in MetS class 1. The highest category in the total number of diagnosis codes out of the list of 68 disease labels (more than 30) increased the total medical cost by 1.66 times compared with the lowest category (less than 10), and its average marginal effect (466,361 yen) was larger than that of any MetS class. The 5-year mortality rate in MetS class 7 was the highest of all classes (odds ratio = 2.47).

[Figure omitted. See PDF.]

Discussion

This study is the first attempt to clarify multimorbidity patterns among working-age patients with high medical costs using LCA. The strength of our study is that the present results of multimorbidity were obtained based on exhaustive insurance claims data and various types of diseases collected from health insurance claims data, which provide reliable results from a representative group of Japanese people.

One reason for the different prevalence of multimorbidity is that there is no consensus on the number of comorbidities used to define multimorbidity. Lee et al. demonstrated that the prevalence of multimorbidity defined as two or more comorbidities was estimated to be higher than that defined as three or more comorbidities [16]. Therefore, when comparing the prevalence rates of previous studies, we need to pay attention to the difference in the definition of multimorbidity. Since the present study focused on patients with high medical costs or severe diseases, the prevalence rate of multimorbidity, which was defined as both 2 or more and 3 or more chronic conditions in this study, was higher than that reported in previous studies that included mild patients [17–20]. Anderson et al. showed that 46.3% of the top 10% high-cost patients had 3 or more chronic conditions [21], and Zulman et al. reported that 77% of the top 5% high-cost patients had multimorbidity defined as 3 or more chronic conditions [7]; the proportion of their patients with multimorbidity was smaller than ours defined as 3 or more chronic conditions (91.1%). The prevalence rates of multimorbidity could not be compared between studies because they depend on the method used for evaluating conditions (e.g., self-report, questionnaire, or health insurance claims data). In particular, the prevalence of multimorbidity tends to be estimated higher if the number of diseases used to define the multimorbidity increases [22]. We included multiple chronic conditions (46 diseases), and there was a possibility that the number of conditions tended to be higher in our study than in other studies. Juul-Larsen et al. recommended including at least the 29 most prevalent chronic conditions when using the chronic condition measurement guide to study multimorbidity [23]. Furthermore, Holzer et al. showed that studies using classifications with under 25 or over 75 chronic conditions tended to yield lower prevalence estimates and proposed choosing a list of chronic conditions that contained 25–75 single conditions [24]. Consequently, it is possible that the prevalence of multimorbidity was underestimated in previous studies that considered a small number of diseases (less than 25) in the analysis.

The present results showed that MetS is particularly important in controlling high medical costs. The prevalence of the MetS class size was 20% among patients in their 30s, but this gradually increased and eventually reached 50% or more after their 50s. The Japanese Ministry of Health, Labour and Welfare introduced a specific health checkup focusing on people older than 40 years of age with MetS in 2008. However, our results suggested that young men in their 30s should also be provided with health education to prevent the aggravation of metabolic disorders. The increase in MetS among women aged approximately 50 years was related to changes in the sex hormone balance associated with menopause [25]. The decrease in estrogen levels in menopausal women is associated with an increased appetite [26], loss of subcutaneous fat, and an increase in intra-abdominal fat, which causes metabolic abnormalities in menopausal women [27]. Although epidemiological studies have suggested an association between osteoarthritis and natural menopause, the causal relationship between sex hormones and osteoarthritis remains unclear [28, 29]. The present results support that obesity [30] and MetS might mediate postmenopausal osteoarthritis. These detailed mechanisms need to be further examined using longitudinal data.

The concept of MetS has been proposed on the basis of pathophysiologies leading to co-occurrences. No studies have examined the actual patterns of MetS complications using real-world data. We identified 7 MetS patterns in Japanese patients with high medical costs. First, patients in MetS class 1 were younger and had fewer comorbidities and lower inpatient costs and 5-year mortality than those in the other MetS classes, indicating that MetS class 1 consisted of patients with mild conditions. MetS class 2 showed almost the same medical cost per person and mortality as MetS class 1 but included some complications characteristic of obesity [31]. In MetS class 3, most patients had motor disorders, and one of the reasons for this might be that muscle dysfunction and mobility impairment are common in patients with long-term complications of diabetes [32]. Conversely, osteoarthritis also causes low physical activity, consequently leading to MetS. It is important for patients with MetS to monitor osteoarthritis early and carefully [33]. MetS classes 4 and 7 included complications in vital organs including the liver and kidney. Evidence linking non-alcoholic fatty liver disease (NAFLD) and CKD has attracted considerable scientific interest in the last decade [34] and our results might support this relationship. However, Akahane et al. reported that NAFLD was not independently associated with CKD in Japanese people [35]. The causality between NAFLD and CKD should be examined in future studies. Finally, MetS classes 5, 6, and 7 showed a higher medical cost per person, higher mortality, and higher rate of inpatients cost due to diseases requiring surgery or many other comorbidities. It is important to control future high medical costs to prevent progression from mild MetS to these critical MetS classes. Herein, we defined the MetS class only if the MetS triad was simultaneously represented, but there were some classes with MetS dyad. Hence, the medical costs associated with metabolic disorders may be higher.

The proportions of mental disease classes were larger among young men and women, and there were 3 types of mental disease classes: mental disease 1, characterized only by mental illness and associated diseases (e.g., schizophrenia, mood disorder, neurotic disorder, and sleep disorder); mental disease 2, and concurrent metabolic diseases such as diabetes and dyslipidemia in addition to mood disorder and sleep disorder; and mental disease 3, a high proportion of women in addition to mood disorder and sleep disorder. Mental disorders are associated with diabetes [36], with a higher prevalence among women than men [37]. A recent study showed that insulin resistance positively predicted incident major depressive disorder [38], and depression was reported to be a risk factor for diabetes in a meta-analysis [39]. Furthermore, Tong et al. showed that the number of multimorbidities increased with the levels of depression, especially in female, young, and middle-aged individuals [40]. In contrast to mental disease class 1, classes 2 and 3 might have developed secondary to other diseases or sex-specific factors. Total costs in the mental disease classes were less, but retirement rates for 5 years were higher than those in other classes. In Japan, mandatory retirement has been adopted since the age of 60. The median age of each mental disease class was in the 40s, and thus, the high retirement rates in these classes have little to do with retirement age. Mental health intervention among the young to middle-aged populations is also important to maintain labor productivity.

Our results showed that the cost per person in the kidney class was by far the most expensive. One of the reasons for this is the high prevalence of dialysis, which is the second highest in the world [41]. Therefore, prevention of dialysis among young or middle-aged people will be important in future health policies. Furthermore, CKD is associated with multimorbidity [42]. Patients in the kidney disease class had the highest number of comorbidities, and half of them were diagnosed with cardiovascular diseases, which further explains the high medical costs.

To address the issues raised by multimorbidity, it is essential to adopt the concept of integrated care that places patients at the center of their healthcare journey. However, the development of integrated care in Japan is limited compared to other countries as there is an absence of a formal gatekeeping system in the country [43]. Previously, most physicians in Japan underwent postgraduate training in specialized clinical departments, leading to an emphasis on specialization and a shortage of primary care physicians. In recent years, the Ministry of Health, Labour and Welfare has reviewed the postgraduate training system to increase primary care physicians, and the number of primary care physicians in Japan might gradually increase in the future. By increasing the number of primary care physicians and promoting a more interconnected healthcare system, Japanese medical policy can transition towards a comprehensive and patient-focused approach to healthcare delivery.

This study has several limitations. First, because of the cross-sectional design, we could not consider time since diagnosis. Second, severity of illness was not obtained from the health insurance claims database. Third, we could not examine the effect of socioeconomic status on multimorbidity patterns because the present database did not include information on income and educational background. Even in Japan where individuals can access health care at 30% of the total cost, a low-income individual has poorer access to outpatient care and more serious health conditions than their higher income counterparts [44], which might cause higher risk of multimorbidity. Future studies need to examine the other insurance associations that cover individuals who are self-employed as well as other daily wage earners. A low education level has been reported to be associated with an increased risk of multimorbidity [45]; in particular, this trend seems to be stronger among younger generations (less than 55 years) [46]. Further studies including SES or other factors such as lifestyle, clinical variables, and functional variables are necessary to completely understand multimorbidity. Fourth, although visceral obesity is a prerequisite for MetS, the obesity code (E66) was not considered for identification of the MetS class. Fifth, the prevalence rate of E66 in this population was 1.7%, which was considerably low compared with the prevalence of obesity (defined in Japan as a body mass index of 25.0 kg/m² or more; male individuals, 28.4%; female individuals, 18.7%) [47]. This is related to the fact that few doctors provide the disease name of obesity when claiming insurance, and there should be more obese patients even if there is no E66. Sixth, we did not include older adults older than 65 years of age for the following reasons. The present database mostly comprises a working age population and includes only 5% of adults older than 65 years of age because the Japanese law stipulates that the employment obligations by enterprises is up to the age of 65. In addition, the co-payment amount is 30% of the total medical cost for those younger than 70 years of age, whereas that is 20% for those older than 70 years of age if their income is not the same as that of the working age population. Kato et al. [48] examined the effects of different co-payments on the utilization of health care by a regression discontinuity design and demonstrated that a lower co-payment increased outpatient expenditure. Because of this, an increase in healthcare service utilization might be induced by not only multimorbidity but also the co-payment model in adults older than 70 years of age in the present database.

Conclusion

Most working-age patients with high medical costs suffered from multimorbidity in Japan. Multimorbidity patterns leading to high medical costs per person were MetS, malignancy, and kidney disease classes. Considering the large number of patients and high total medical expenses, dealing with MetS is particularly important to control the high medical costs. In particular, it is necessary to prevent the aggravation of MetS comorbidities, such as cardiovascular diseases or malignancies, in Japan’s future medical policy.

Supporting information

S1 Checklist. Strengthening the reporting of observational studies in epidemiology statement—A checklist of items that should be included in reports of observational studies.

https://doi.org/10.1371/journal.pone.0291554.s001

(DOC)

S1 Fig. Distribution of annual medical costs in the insured population who were subscribers in 2015 (n = 16,989,029).

The top 1%, 5%, 10%, 20%, 30% of patients accounted for 26.1%, 46.9%, 59.0%, 73.9% and 83.3% of the total annual medical costs. 1 USD = 120 JPY.

https://doi.org/10.1371/journal.pone.0291554.s002

(DOCX)

S2 Fig. Model fit for the latent class analysis (BIC and AIC).

BIC = Bayesian Information Criterion; AIC = Akaike Information Criterion.

https://doi.org/10.1371/journal.pone.0291554.s003

(DOCX)

S3 Fig. A 100% stacked vertical bar of the patient count in each class for 20, 25, 30, 35, and 40 latent class models.

https://doi.org/10.1371/journal.pone.0291554.s004

(DOCX)

S4 Fig. Heatmap for item-response probabilities of disease labels in the 30 latent class models.

Thick red represents a high probability.

https://doi.org/10.1371/journal.pone.0291554.s005

(DOCX)

S1 Table. Prevalences of each ICD10 code among patients with high medical costs.

ICD10 = International Statistical Classification of Diseases and Related Health Problems 10th Revision.

https://doi.org/10.1371/journal.pone.0291554.s006

(XLSX)

S2 Table. Forty-six chronic disease labels classified by the ICD10 codes.

https://doi.org/10.1371/journal.pone.0291554.s007

(XLSX)

S3 Table. Prevalences of disease labels in patients with the top 10% cost and low to medium cost.

https://doi.org/10.1371/journal.pone.0291554.s008

(XLSX)

S4 Table. Demographic characteristics of the patients with the top 10% cost and low to medium cost.

https://doi.org/10.1371/journal.pone.0291554.s009

(XLSX)

Acknowledgments

The authors would like to thank Takehiko Baba and Yosuke Ihara for their technical assistance and knowledge of insurance societies. We are also grateful to Nobuki Ando, the chairman of the Japan Health Insurance Association, and all the relevant staff for providing the opportunity to conduct this study.

Citation: Nishida Y, Anzai T, Takahashi K, Kozuma T, Kanda E, Yamauchi K, et al. (2023) Multimorbidity patterns in the working age population with the top 10% medical cost from exhaustive insurance claims data of Japan Health Insurance Association. PLoS ONE 18(9): e0291554. https://doi.org/10.1371/journal.pone.0291554

About the Authors:

Yuki Nishida

Roles: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft

Affiliations: Department of Biostatistics, M&D Data Science Center, Tokyo Medical and Dental University, Tokyo, Japan, Graduate School of Health Management, Keio University, Yokohama, Kanagawa, Japan, Sports Medicine Research Center, Keio University, Yokohama, Kanagawa, Japan

ORICD: https://orcid.org/0000-0002-8056-1515

Tatsuhiko Anzai

Roles: Conceptualization, Data curation, Formal analysis, Methodology, Writing – review & editing

Affiliation: Department of Biostatistics, M&D Data Science Center, Tokyo Medical and Dental University, Tokyo, Japan

ORICD: https://orcid.org/0000-0003-1073-2996

Kunihiko Takahashi

Roles: Conceptualization, Funding acquisition, Supervision, Writing – review & editing

Affiliation: Department of Biostatistics, M&D Data Science Center, Tokyo Medical and Dental University, Tokyo, Japan

ORICD: https://orcid.org/0000-0003-2387-7772

Takahide Kozuma

Roles: Conceptualization, Writing – review & editing

Affiliation: Department of Internal Medicine, School of Medicine, Keio University, Tokyo, Japan

Eiichiro Kanda

Roles: Conceptualization, Funding acquisition, Writing – review & editing

Affiliation: Medical Science, Kawasaki Medical School, Okayama, Japan

Keita Yamauchi

Roles: Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing

Affiliation: Graduate School of Health Management, Keio University, Yokohama, Kanagawa, Japan

Fuminori Katsukawa

Roles: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing

E-mail: [email protected]

Affiliation: Sports Medicine Research Center, Keio University, Yokohama, Kanagawa, Japan

ORICD: https://orcid.org/0000-0002-6122-8696

References

1. Ministry of Health, Labour and Welfare (Japan). https://www.mhlw.go.jp/english/org/policy/p34-35.html.

2. Kojima G. Increased healthcare costs associated with frailty among community-dwelling older people: A systematic review and meta-analysis. Arch Gerontol Geriatr. 2019;84: 103898. pmid:31228673

3. Kasajima M, Eggleston K, Kusaka S, Matsui H, Tanaka T, Son BK, et al. Projecting prevalence of frailty and dementia and the economic cost of care in Japan from 2016 to 2043: a microsimulation modelling study. Lancet Public Health. 2022;7: e458–e468. pmid:35487231

4. Bloom DE, Cafiero ET, Jané-Llopis E, Abrahams-Gessel S, Bloom LR, Fathima S, et al. The Global Economic Burden of Noncommunicable Diseases. Geneva: World Economic Forum; 2011.

5. Cohen SB, Uberoi N. Differentials in the Concentration in the Level of Health Expenditures across Population Subgroups in the U.S., 2010. Statistical Brief (Medical Expenditure Panel Survey (US)). Rockville (MD): Agency for Healthcare Research and Quality (US); 2001.

6. Osawa I, Goto T, Yamamoto Y, Tsugawa Y. Machine-learning-based prediction models for high-need high-cost patients using nationwide clinical and claims data. NPJ Digit Med. 2020;3: 148. pmid:33299137

7. Zulman DM, Pal Chee C, Wagner TH, Yoon J, Cohen DM, Holmes TH, et al. Multimorbidity and healthcare utilisation among high-cost patients in the US Veterans Affairs Health Care System. BMJ Open. 2015;5: e007771. pmid:25882486

8. van den Akker M, Buntinx F, Roos S, Knottnerus JA. Problems in determining occurrence rates of multimorbidity. J Clin Epidemiol. 2001;54: 675–679. pmid:11438407

9. Wang L, Si L, Cocker F, Palmer AJ, Sanderson K. A Systematic Review of Cost-of-Illness Studies of Multimorbidity. Appl Health Econ Health Policy. 2018;16:15–29. pmid:28856585

10. Tran PB, Kazibwe J, Nikolaidis GF, Linnosmaa I, Rijken M, van Olmen J. Costs of multimorbidity: a systematic review and meta-analyses. BMC Med. 2022;20:234. pmid:35850686

11. Albreht T, Dyakova M, Schellevis FG, Van den Broucke S. Many diseases, one model of care? J Comorb. 2016;6: 12–20. pmid:29090167

12. Zou S, Wang Z, Bhura M, Tang K. Association of multimorbidity of non-communicable diseases with mortality: a 10-year prospective study of 0.5 million Chinese adults. Public Health. 2022;205: 63–71. pmid:35240390

13. Nylund-Gibson K, Garber AC, Singh J, Witkow MR, Nishina A, Bellmore A. The Utility of Latent Class Analysis to Understand Heterogeneity in Youth Coping Strategies: A Methodological Introduction. Behavioral Disorders. 2022;0: 01987429211067214.

14. World Health Organization. Multimorbidity: Technical series on safer primary care 2016. https://apps.who.int/iris/bitstream/handle/10665/252275/9789241511650-eng.pdf.

15. Malehi AS, Pourmotahari F, Angali KA. Statistical models for the analysis of skewed healthcare cost data: a simulation study. Health Econ Rev. 2015;5: 11. pmid:26029491

16. Lee ES, Lee PSS, Xie Y, Ryan BL, Fortin M, Stewart M. The prevalence of multimorbidity in primary care: a comparison of two definitions of multimorbidity with two different lists of chronic conditions in Singapore. BMC Public Health. 2021;21: 1409. pmid:34271890

17. Hajat C, Siegal Y, Adler-Waxman A. Clustering and Healthcare Costs With Multiple Chronic Conditions in a US Study. Front Public Health. 2020; 8:607528. pmid:33553094

18. Nedzinskienė L, Jurevičienė E, Visockienė Ž, Ulytė A, Puronaitė R, Kasiulevičius V, et al. Structure and Distribution of Health Care Costs across Age Groups of Patients with Multimorbidity in Lithuania. Int J Environ Res Public Health. 2021;18. pmid:33803298

19. Heins M, Korevaar J, Schellevis F, Rijken M. Identifying multimorbid patients with high care need—A study based on electronic medical record data. Eur J Gen Pract. 2020;26: 189–195. pmid:33337928

20. Violán C, Foguet-Boreu Q, Roso-Llorach A, Rodriguez-Blanco T, Pons-Vigués M, Pujol-Ribera E, et al. Burden of multimorbidity, socioeconomic status and use of health services across stages of life in urban areas: a cross-sectional study. BMC Public Health. 2014;14: 530. pmid:24885174

21. Anderson M, Revie CW, Stryhn H, Neudorf C, Rosehart Y, Li W, et al. Defining ’actionable’ high- costhealth care use: results using the Canadian Institute for Health Information population grouping methodology. Int J Equity Health. 2019;18: 171. pmid:31707981

22. Simard M, Rahme E, Calfat AC, Sirois C. Multimorbidity measures from health administrative data using ICD system codes: A systematic review. Pharmacoepidemiol Drug Saf. 2022;31: 1–12. pmid:34623723

23. Juul-Larsen HG, Christensen LD, Andersen O, Bandholm T, Kaae S, Petersen J. Development of the "chronic condition measurement guide": a new tool to measure chronic conditions in older people based on ICD-10 and ATC-codes. Eur Geriatr Med. 2019;10: 431–444. pmid:34652799

24. Holzer BM, Siebenhuener K, Bopp M, Minder CE. Evidence-based design recommendations for prevalence studies on multimorbidity: improving comparability of estimates. Popul Health Metr. 2017;15: 9. pmid:28270157

25. Hallajzadeh J, Khoramdad M, Izadi N, Karamzad N, Almasi-Hashiani A, Ayubi E, et al. Metabolic syndrome and its components in premenopausal and postmenopausal women: a comprehensive systematic review and meta-analysis on observational studies. Menopause. 2018;25(10): 1155–1164. pmid:29787477

26. Hirschberg AL. Sex hormones, appetite and eating behaviour in women. Maturitas. 2012;71: 248–256. pmid:22281161

27. Lizcano F, Guzmán G. Estrogen Deficiency and the Origin of Obesity during Menopause. Biomed Res Int. 2014;2014: 757461. pmid:24734243

28. de Klerk BM, Schiphof D, Groeneveld FP, Koes BW, van Osch GJ, van Meurs JB, et al. No clear association between female hormonal aspects and osteoarthritis of the hand, hip and knee: a systematic review. Rheumatology (Oxford). 2009;48: 1160–1165. pmid:19608726

29. Dennison EM. Osteoarthritis: The importance of hormonal status in midlife women. Maturitas. 2022;165: 8–11. pmid:35841775

30. Guh DP, Zhang W, Bansback N, Amarsi Z, Birmingham CL, Anis AH. The incidence of co-morbidities related to obesity and overweight: a systematic review and meta-analysis. BMC Public Health. 2009;9: 88. pmid:19320986

31. Kivimäki M, Strandberg T, Pentti J, Nyberg ST, Frank P, Jokela M, et al. Body-mass index and risk of obesity-related complex multimorbidity: an observational multicohort study. Lancet Diabetes Endocrinol. 2022;10: 253–263. pmid:35248171

32. Bianchi L, Volpato S. Muscle dysfunction in type 2 diabetes: a major threat to patient’s mobility and independence. Acta Diabetol. 2016;53: 879–889. pmid:27393005

33. Liu SY, Zhu WT, Chen BW, Chen YH, Ni GX. Bidirectional association between metabolic syndrome and osteoarthritis: a meta-analysis of observational studies. Diabetol Metab Syndr. 2020;12: 38. pmid:32399062

34. Marcuccilli M, Chonchol M. NAFLD and Chronic Kidney Disease. Int J Mol Sci. 2016;17: 562. pmid:27089331

35. Akahane T, Akahane M, Namisaki T, Kaji K, Moriya K, Kawaratani H, et al. Association between Non-Alcoholic Fatty Liver Disease and Chronic Kidney Disease: A Cross-Sectional Study. J Clin Med. 2020;9. pmid:32481684

36. Roy T, Lloyd CE. Epidemiology of depression and diabetes: a systematic review. J Affect Disord. 2012;142 Suppl: S8–S21. pmid:23062861

37. Salk RH, Hyde JS, Abramson LY. Gender differences in depression in representative national samples: Meta-analyses of diagnoses and symptoms. Psychol Bull. 2017;143: 783–822. pmid:28447828

38. Watson KT, Simard JF, Henderson VW, Nutkiewicz L, Lamers F, Nasca C, et al. Incident Major Depressive Disorder Predicted by Three Measures of Insulin Resistance: A Dutch Cohort Study. Am J Psychiatry. 2021;178: 914–920. pmid:34551583

39. Yu M, Zhang X, Lu F, Fang L. Depression and Risk for Diabetes: A Meta-Analysis. Can J Diabetes. 2015;39: 266–272. pmid:25773933

40. Tong L, Pu L, Guo X, Sun M, Guo F, Zhao S, et al. Multimorbidity study with different levels of depression status. J Affect Disord. 2021;292: 30–35. pmid:34091380

41. USRD System. 2020 USRDS Annual Data Report: Epidemiology of kidney disease in the United States. Bethesda: National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases; 2020.

42. Gurgel do Amaral MS, Reijneveld SA, Meems LMG, Almansa J, Navis GJ, de Winter AF. Multimorbidity prevalence and patterns and their associations with health literacy among chronic kidney disease patients. J Nephrol. 2022;35: 1709–1719. pmid:34985613

43. Kaneko M, Shinoda S, Shimizu S, Kuroki M, Nakagami S, Chiba T, et al. Fragmentation of ambulatory care among older adults: an exhaustive database study in an ageing city in Japan. BMJ Open. 2022;12: e061921. pmid:35953252

44. Fujita M, Sato Y, Nagashima K, Takahashi S, Hata A. Income Related Inequality of Health Care Access in Japan: A Retrospective Cohort Study. PLoS One. 2016;11: e0151690. pmid:26978270

45. Pathirana TI, Jackson CA. Socioeconomic status and multimorbidity: a systematic review and meta-analysis. Aust N Z J Public Health. 2018;42(2):186–94. Epub 2018/02/15. pmid:29442409.

46. Afshar S, Roderick PJ, Kowal P, Dimitrov BD, Hill AG. Multimorbidity and the inequalities of global ageing: a cross-sectional study of 28 countries using the World Health Surveys. BMC Public Health. 2015;15: 776. pmid:26268536

47. Ministry of Health, Labour and Welfare (Japan): The National Health and Nutrition Survey in Japan, 2015. 2017.

48. Kato H, Goto R, Tsuji T, Kondo K. The effects of patient cost-sharing on health expenditure and health among older people: Heterogeneity across income groups. Eur J Health Econ. 2022;23: 847–861. pmid:34779932