- 1Chongqing General Hospital, Chongqing, China
- 2Center for Mental Health of Rongchang District, Chongqing, China
- 3Health Committee of Rongchang District, Chongqing, China
Background: The prevalence, patterns, and impact of multimorbidity on health-related quality of life (HRQoL) remain inadequately understood among rural populations in southwest China. This study seeks to fill this knowledge gap by systematically examining these aspects.
Methods: Participants were recruited from the China Multi-Ethnic Cohort (CMEC) study. Incident cases of 13 chronic conditions were documented. Multimorbidity was defined as the presence of two or more chronic conditions in an individual. Principal component factor analysis (PCFA) was performed to identify patterns of multimorbidity. Tobit regression analysis and restricted cubic spline were employed to assess the correlation between multimorbidity patterns and HRQoL.
Results: A total of 2,998 participants were enrolled, with a mean age of 50.65 years (SD = 11.99). The prevalence of multimorbidity was 48.50%. Four multimorbidity patterns were identified by PCFA: circulatory system pattern, digestive system pattern, metabolic syndrome pattern, and hepatobiliary system pattern. All four patterns were negatively correlated with HRQoL, as demonstrated by tobit regression analysis (β = −0.024, β = −0.020, β = −0.007, β = −0.018; all p < 0.001). Restricted cubic spline also demonstrated the negative correlation between different multimorbidity patterns and HRQoL, after adjusting for potential confounding factors. Subgroup analysis in different gender, age, and average yearly family total income also demonstrated these negative correlations.
Conclusion: The prevalence of multimorbidity is relatively high in rural southwest China. Distinct multimorbidity patterns were correlated with poorer HRQoL. These findings enhance the understanding of multimorbidity patterns and may inform the development of tailored primary healthcare services.
Introduction
Multimorbidity is commonly defined as the co-occurrence of multiple long-term health conditions in the same individual, requiring ongoing and diverse treatment (1). The prevalence of chronic diseases is rising partly due to improved health literacy and advancements in diagnostic technologies (2). The co-occurrence of multiple chronic conditions can result in complex interactions that not only significantly reduce the quality of life (3) but also impose a substantial economic burden on both families and society. Compared to individuals with a single chronic condition, patients with multimorbidity face higher risk of mortality (4) and increased healthcare costs, greater demand for medical services, and more frequent polypharmacy (5–7). Research on multimorbidity is crucial for advancing patient self-management, optimizing treatment outcomes, and facilitating the judicious allocation of healthcare facilities and resources, as well as informing policymaking and clinical guidance. Multimorbidity has emerged as a significant public health challenge in China, necessitating urgent attention and resolution.
Recent research shows increasing interest in multimorbidity among the middle-aged and elderly population. In the United States, over half of older adults suffer from three or more chronic conditions (8), while in China, the prevalence of multimorbidity among adults in China is 58.10% (9). Studies consistently link the number of chronic diseases to health-related quality of life (HRQoL) in the older adults (7) with systematic reviews confirming the inverse relationship between multimorbidity conditions and HRQoL (10, 11). While the body of literature on multimorbidity and HRQoL is growing, fewer studies explore multimorbidity patterns affecting HRQoL in the general population and the contributing factors. The main challenge in assessing multimorbidity risk lies in identifying and classifying its subtypes. Numerous studies have linked multimorbidity patterns to functional decline (12), depression (13), dosage frequency (14), and healthcare utilization (15). Common multimorbidity patterns are often identified in older population. Three major multimorbidity patterns were identified from 2015 waves of CHARLS study, including asthma/chronic lungs diseases, asthma, arthritis, or rheumatism/chronic lung diseases, dyslipidemia, hypertension, arthritis, or rheumatism/heart attack (5). Despite limited research on multimorbidity prevalence and patterns in the general population, a review summarized 51 studies and identified 407 profiles, and cardio-metabolic syndrome and mental health were more consistently observed in multimorbidity pattern studies (16).
HRQoL refers to a multidimensional concept that reflects individuals’ experiences across different cultures in relation to their goals, expectations, standards, concerns, and overall life circumstances. It provides a comprehensive evaluation of an individual’s disease burden or health status from physiological, psychological, and social perspectives. It has been extensively applied in assessing health conditions, evaluating the effectiveness of health services, and guiding the development of health service programs, thereby promoting the utilization of medical and healthcare resources. Although HRQoL loss associated with single chronic diseases has been well studied, research focusing on the impact of multimorbidity remains limited (17). Previous studies have identified a negative correlation between multiple chronic conditions and HRQoL (10). However, detailed knowledge about the influence of specific multimorbidity patterns on HRQoL is still lacking.
Previous studies have primarily focused on nationally representative samples of older Chinese adults. Recently, research has shifted from examining single diseases to multimorbidity in relation to health status. However, few studies have explored the correlation between multimorbidity patterns and HRQoL. This study was undertaken to examine the prevalence of multimorbidity and identify its common patterns in the general population of rural southwest China.
Methods
Study design and population
The data used in this study were obtained from the baseline survey of the China Multi-Ethnic Cohort Study (18), which was the largest cohort study conducted by Sichuan University in southwestern China from September 2018 to January 2019. Briefly in Rongchang region (19), participants were recruited via a three-stage stratified random sampling approach. Initially, four streets named Changyuan, Changzhou, Anfu, and Guangshun were randomly chosen from a total of 21 streets or towns. Subsequently, 10 villages were randomly selected from each of these streets. Finally, 50 to 80 individuals were randomly sampled in accordance with the age and sex distribution of the Rongchang population. Ethics approval was granted by the ethics committee of Sichuan University (No. K2016038), and all participants provided written informed consent.
The recruitment of participants in this study was guided by the following inclusion criteria: (1) age between 30 and 79 years at the time of the investigation; (2) having resided in Rongchang for at least 6 months; (3) Han ethnicity; (4) voluntarily participating in the survey, providing written informed consent, agreeing to provide biological samples, and committing to complete follow-up procedures; and (5) having no mental illness or cognitive disorders and possessing intact expressive abilities. Individuals were excluded if they met any of the following criteria: incomplete data regarding general characteristics, questionnaires, physical examinations, or blood biochemical tests. A total of 3,002 individuals were initially recruited during the baseline survey. After excluding those with missing data on general characteristics, the final analysis included 2,998 participants.
Data collection
General characteristics
Age was expressed as mean and standard deviation and grouped into “30–39,” “40–49,” “50–59,” “60–69,” and “70–79” years. Gender was grouped into “males” and “females.” Marital status was grouped into “married/cohabiting” and “separated/divorced/widowed/unmarried.” Education level was grouped into “primary school or below,” “junior middle school,” and “high school or above.” Job status was grouped into “farmers,” “government employee,” “workers,” “sales & service staffs,” and “others.” Average yearly family total income was grouped into “<20,000 yuan,” “20,000–59,999 yuan,” “60,000–99,999 yuan,” and “≥100,000 yuan.”
Definition of healthy life factors
Non-smoking, moderate alcohol intake, healthy waist-to-hip ratio (WHR), adequate physical activity (PA), the dietary approaches to stop hypertension (DASH) score, and adequate sleep were considered six modifiable healthy lifestyle factors based on literature reviews (20, 21). Non-smoking was defined as participants who had never consumed any cigarette products. Moderate alcohol intake was defined as alcohol consumption within 1–14 g/day for females or 1–28 g/day for males (22) according to the self-reported alcohol consumption data. Healthy WHR was defined as <0.90 for males or <0.85 for females (20). WHR is calculated by dividing waist circumference by hip circumference. Physical activity levels, including work, transportation, household chores, and leisure-time activities, were assessed by metabolic equivalent tasks (METs) (23). Adequate PA was defined as having a total METs greater than 3,000 (19). A modified DASH diet score was calculated using a method adapted from prior study (24) with minor adjustments tailored to the CEMC dataset (25). The modified DASH score was constructed based on seven food categories: whole grains, fresh fruits, fresh vegetables, beans, dairy, red meat products, and sodium. Scores ranging from 1 to 5 were allocated according to the quintiles of average consumption for each food group. Specifically, participants in the highest quintile of intake for whole grains, fresh fruits, fresh vegetables, beans, and dairy were assigned a score of 5, while those in the lowest quintile of intake for red meat products and sodium received a score of 1. The overall DASH score for each participant was then calculated by summing the scores from all seven food groups (26). Based on the lowest tertile of the DASH score as the cutoff, participants were categorized into two groups: “<20” and “≥20.” Adequate sleep was defined as a sleep duration of 7–8 h/day (27).
Measurement of chronic diseases
A total of 21 chronic diseases or conditions were collected, including hypertension, diabetes, dyslipidemia, coronary heart disease, stroke, rheumatic heart diseases, pulmonary heart disease, pulmonary tuberculosis, chronic bronchitis/pulmonary emphysema, asthma, chronic hepatitis/cirrhosis, ulcers of the digestive tract, gastroenteritis, gallstones/cholecystitis, bone fracture, rheumatoid or arthritis, intervertebral disc disease, mental and psychological disorders, neurasthenia, brain injury, and cancers. The incidence rate of each chronic diseases or conditions is shown in Supplementary Table 1. To improve the robustness of the principal component factor analysis, chronic diseases with an incidence rate below 1% were excluded from the current study based on previous study (28).
Hypertension (29) was identified as having an average systolic blood pressure (SBP) of ≥140 mmHg or diastolic blood pressure (DBP) of ≥90 mmHg, based on three consecutive blood pressure measurements taken at 5-min intervals while the individual was at rest. It was also defined as having a prior physician diagnosis of hypertension or current use of blood pressure-lowering treatments (e.g., medications and exercise). Diabetes (30) was diagnosed when fasting plasma glucose (FPG) levels were above 7.00 mmol/L, or when there was a prior physician diagnosis of diabetes mellitus or the individual was on treatments to lower FPG levels (e.g., medications and exercise). Dyslipidemia was defined by the presence of any of the following four abnormalities (31): (1) total cholesterol ≥6.2 mmol/L; (2) triacylglycerol≥2.3 mmol/L; (3) low-density lipoprotein cholesterol≥4.1 mmol/L; (4) high-density lipoprotein cholesterol < 1.0 mmol/L.
Hypertension, diabetes, and dyslipidemia were verified through clinical or biometric methods, while other chronic diseases or conditions were self-reported by the participants. Coronary heart disease, rheumatic heart disease, and pulmonary heart disease were grouped into heart diseases. Multimorbidity was defined as the coexistence of multiple chronic conditions greater than or equal to two.
Health-related quality of life
HRQoL was measured using the European Five-Dimensional Five-level Health Scale (EQ-5D-5L) instruments (32), which enhances the ability to differentiate respondents’ health status and enhances objectivity in evaluating their subjective experiences. The scale consists of the EQ-5D health description system and the EQ-5D visual analogue scale. HRQoL serves as a multidimensional metric for assessing physical condition, mental functioning, social capability, and overall personal wellbeing, thereby reflecting the impact of diseases or physical and mental impairments on an individual’s overall health status. The EQ-5D health description system encompasses five dimensions: mobility, self-care, usual activities, pain/discomfort, and anxiety/depression. Each dimension is divided into five levels: no difficulty, slight difficulty, moderate difficulty, severe difficulty, and unable to complete/extreme difficulty. The five response levels (1–5) correspond to each dimension in sequence. Initially, a five-level health profile is generated (e.g., 55,555, indicating extreme difficulty or inability in all five dimensions), allowing the EQ-5D to represent a total of 3,125 (55) distinct health states (33). Subsequently, the EQ-5D converts raw data into health utility values through a utility scoring system developed by Luo (34). The scale ranges from −0.391 to 1, where 0 corresponds to death and 1 denotes perfect health, and a health utility value below 0 indicates a health state deemed worse than death (35). The Chinese version of the EQ-5D-5L has been validated by Chinese researchers and has demonstrated satisfactory reliability and validity. The methodology for deriving EQ-5D-5L values is detailed in the Supplementary Tables 2, 3.
Statistical analysis
Age and HRQoL were reported using mean and standard deviation, despite the absence of a normal distribution. General characteristics and healthy life factors were presented by constituent ratio (%), difference in HRQoL between binary variables was analyzed by Mann–Whitney U test, and differences among multivariate variables were analyzed by Kruskal–Wallis H test. Principal component factor analysis (PCFA) is one of the statistical approaches to identify non-random cluster patterns; we used PCFA with varimax rotation based on the same chronic conditions to z scores (25). After comprehensively considering the Kaiser–Meyer–Olkin measure of sample adequacy, Bartlett’s test of sphericity, eigenvalues, cumulative variance explained, scree plot, and interpretability into considerations, four major multimorbidity patterns were identified. For each multimorbidity pattern, factor scores were calculated for all participants by summing up the standardized chronic conditions weighted by their factor loadings. The multimorbidity patterns were assessed by a split-sample validation to support the robustness of the characteristics. Given the ability of tobit regression to effectively handle samples with limited data ranges and skewed distributions (such as health utility value was skewed and censored at 1) (28, 36). The influencing factors and multimorbidity patterns correlated with HRQoL were conducted by tobit regression.
The statistical significance of the tobit regression model was assessed using the likelihood ratio test. In addition, the model fit was evaluated based on the Akaike information criterion (AIC) and the Bayesian information criterion (BIC), with lower values indicating superior model performance. Restricted cubic spline was performed to examine the correlation between factor scores of different multimorbidity patterns and health utility values, after adjusting for potential confounding factors. To assess multicollinearity among the variables, we employed the variance inflation factor (VIF) and tolerance. A VIF value exceeding 5 was indicative of significant collinearity among the variables. In addition, tolerance was calculated as the reciprocal of the VIF (tolerance = 1/VIF). When the tolerance value was less than 0.2, it suggested the presence of substantial collinearity. Subgroup analyses in different gender, age, and average yearly family total income were also conducted to further assess the robustness of the results by tobit regression and restricted cubic spline. The statistical analysis was performed using R software (version 4.1.1), with a two-tailed p-value of less than 0.05 indicating a statistically significant difference.
Results
General characteristics and HRQoL of the participants
General characteristics and HRQoL of the participants are presented in Table 1. Among the 2,998 participants, the mean age was 50.65 (S.D. = 11.99) years and 50.03% of the participants were females. The mean health utility value was 0.95 (S.D. = 0.11), and 2020 of the participants was in the “perfect health” (health utility value = 1). Statistically significant differences were observed across various demographic variables, including gender, age, marital status, education level, job status, and average yearly family total income (all P<0.05).
The healthy life factors and their difference in HRQoL
The healthy life factors and their difference in HRQoL are presented in Table 2. Among the 2,998 participants, 80.19% of the participants were non-smoking, 43.40% of the participants had a healthy WHR, 69.51% of the participants had adequate physical activities, 71.08% of the participants had a higher DASH score over 20, and 62.74% of the participants had adequate sleep. The HRQoL in non-smoking, healthy WHR, adequate PA, and adequate sleep are statistically significant difference (all P<0.05).
Prevalence of multimorbidity, multimorbidity patterns, and its difference in HRQoL
The cases of each chronic disease or condition, both independently and in the context of multimorbidity, are presented in Table 3. A total of 732 participants (24.42%) were free from chronic disease or condition, while 812 participants (27.08%) were suffering from a single chronic disease or condition. The prevalence of multimorbidity was 48.50%. Participants suffering from rheumatoid arthritis, peptic ulcers, or heart disease exhibited a higher likelihood of experiencing multimorbidity, while those with chronic hepatitis/cirrhosis, dyslipidemia, and hypertension were less likely to have multimorbidity. Spearman’s correlation analysis found a negative correlation between the total number of chronic diseases or conditions and HRQoL (rs = −0.349, P<0.001).

Table 3. Each chronic disease or condition independent of comorbidities (all cases) and by each chronic disease with multimorbidity.
Before performing the principal component factor analysis, we evaluated the suitability of the data by calculating the Kaiser–Meyer–Olkin (KMO) measure of sampling adequacy (KMO = 0.677, cutoff > 0.50) and performing Bartlett’s test of sphericity (χ2 = 1396.26, p < 0.001). Both tests indicated significant correlations among the different chronic diseases or conditions, supporting the appropriateness of applying a factor analytic approach. The findings confirmed that the sample was appropriate for the generating distinct and reliable factors. Principal component factor analysis was conducted on the 13 chronic diseases or conditions, resulting in four factors with eigenvalues greater than 1 (1.904, 1.328, 1.099, and 1.050 for the first through fourth factors, respectively). These factors were extracted following varimax rotation with Kaiser normalization. Collectively, the four factors explained 41.39% of the cumulative variance (Supplementary Table 4). Furthermore, the inflection points observed in the scree plot supported the presence of four factors, further substantiating the appropriateness of the factor extraction (Supplementary Figure 1).
Figure 1A shows the multimorbidity pattern in the analysis, 4 common patterns of multimorbidity were observed, namely, circulatory system pattern (which predominantly consists of chronic bronchitis/pulmonary emphysema, heart disease, rheumatoid arthritis, neurasthenia, and hypertension), digestive system pattern (which predominantly consists of ulcers of the digestive tract, gastroenteritis, neurasthenia, rheumatoid arthritis, and intervertebral disc disease), metabolic syndrome pattern (which predominantly consists of dyslipidemia, diabetes, and hypertension), and hepatobiliary system pattern (which predominantly consists of gallstones/cholecystitis, chronic hepatitis/cirrhosis, bone fracture, and intervertebral disc disease). The loading factor of each multimorbidity pattern is shown in Supplementary Table 4.

Figure 1. Exploratory factor analysis for 13 chronic diseases or conditions [(A) loading factor of four common patterns of multimorbidity; (B) the difference in health utility value from Q1 to Q5 in different multimorbidity pattern; (C) The difference of heath utility value in the same group of different multimorbidity pattern].
Figure 1B shows the difference in health utility value from Q1 to Q5 in different multimorbidity patterns. Q5 always had a relatively lower HRQoL in circulatory system pattern, digestive system pattern, and hepatobiliary system pattern (P<0.05).
Figure 1C shows the difference of health utility value in the same group of different multimorbidity patterns. There were statistically significant difference among Q1, Q2, Q3, and Q5 groups (P<0.05). In Q5 group, circulatory system pattern had the lowest health utility value.
The correlation between multimorbidity patterns and HRQoL by tobit regression and restricted cubic spline
The correlation between multimorbidity patterns and HRQoL by tobit regression is presented in Table 4. Multicollinearity diagnostics revealed the absence of multicollinearity among the variables. Specifically, the VIF values for all variable indicators were below 5, while the corresponding tolerance values exceeded 0.2 (Supplementary Table 5). The model fit of the tobit regression was evaluated using the −2 log likelihood ratio, AIC, and BIC. The corresponding results are presented in Supplementary Table 6. After adjusting for confounding factors, the adjusted model 2 demonstrated that the circulatory system pattern (β = −0.024, P < 0.001), digestive system pattern (β = −0.020, P < 0.001), metabolic syndrome pattern (β = −0.007, P < 0.001), and hepatobiliary system pattern (β = −0.018, P < 0.001) were negatively correlated with HRQoL. Restricted cubic spline also demonstrated the negative correlation between factor score of different multimorbidity patterns and health utility values, after adjusting for potential confounding factors, including basic characteristics and life factors (Supplementary Figure 2). Subgroup analysis in different gender, age, and average yearly family total income by tobit regression (Table 5) and restricted cubic spline analysis in different gender (Supplementary Figure 3), age (Supplementary Figure 4), and average yearly family total income (Supplementary Figure 5) also demonstrated those similar negative correlation.

Table 5. Subgroup tobit regression analysis among the correlation between multimorbidity patterns and HRQoL, [β(95% CI)].
Discussion
To our knowledge, this is the first study to explore the prevalence of multimorbidity and examine its common patterns correlated with HRQoL among the general population in rural southwest China. The mean health utility value among the participants was 0.95, which was higher than the German population (0.90) (37) and South Australian population (0.91) (38). Differences in HRQoL may be attributed to the study populations, assessment methods, and levels of regional economic development. Multimorbidity occurred in 48.50% of the participants aged 30–79 years, which was slightly higher than the older Chinese aged over 50 years (42.40%) (39) and Japanese adult population (29.9%) (14). Although the impact of multiple chronic conditions has long been recognized in European health systems, chronic disease management in most countries still relies on single-disease control programs and fails to develop integrated health management for multiple chronic conditions (40). Therefore, it is crucial to identify the patterns of multiple chronic conditions and develop integrated interventions tailored to various combinations. This strategy aims to address the increasing prevalence of multimorbidity in China and improve overall health status. In the current study, four distinct patterns were identified through principal component factor analysis. These findings aim to meet the growing demand from epidemiologists, clinicians, and policymakers for a deeper understanding of the burden of multimorbidity, particularly its clinical and social implications. Moreover, health economists can apply these results in cost–benefit analyses to inform reimbursement decisions.
The current study demonstrated that individuals with rheumatoid or arthritis, ulcers of the digestive tract, and heart disease had the highest prevalence of multimorbidity, while individuals with chronic hepatitis/cirrhosis, dyslipidemia, and hypertension had the lowest multimorbidity prevalence. Yao (39) et al. reported the 91.60% of multimorbidity prevalence in heart problems, 77.50% in hypertension patients, and 75.90% in arthritis or rheumatism patients. Mino-Leon (41) et al. demonstrated 89.90% of multimorbidity prevalence in cardiac problems and 76.30% in hypertension patients. Marengoni (42) et al. found 97.90% of multimorbidity prevalence in subjects with heart failure, 95.20% in subjects with hip fracture, and 74.60% in subjects with chronic obstructive pulmonary disease. The difference in multimorbidity prevalence may be attributed to variations in disease identification, study populations, or methods of analyzing patterns of multimorbidity. In this study, we found a less reported higher prevalence of multimorbidity (95.97%) in rheumatoid or arthritis than a nationally representative older Chinese (75.90%) (39), which was consistent with a study in middle-aged and elderly populations in Xinjiang (43). Chongqing, located in southwestern China, is renowned for its hot pot, a globally recognized dish characterized by high oil and salt content. The substantial consumption of meat contributes to excessive purine intake. Although the humidity in Chongqing is relatively low, the local population believes that the spicy and numbing nature of hot pot induces sweating, which helps eliminate moisture from the body (19, 24).
The findings suggest that patterns of disease combinations can influence the adverse outcomes of HRQoL. Four common patterns of multimorbidity were identified through PCFA, namely, circulatory system pattern, digestive system pattern, metabolic syndrome pattern, and hepatobiliary system pattern, which were all significantly correlated with reduced HRQoL. Among the identified multimorbidity patterns in the current study, circulatory and metabolic system were replicable patterns in previous studies (16, 44). Differences in multimorbidity patterns may be attributed to variations in study design, sample selection, social demographic characteristics, eligible diseases, and statistical methods. The circulatory system pattern, which primarily includes chronic bronchitis/pulmonary emphysema and heart diseases, observed in current study, is consistent with the findings from previous research studies (39, 42). Previous studies have documented the increasing prevalence of heart diseases in China, accompanied by significant complications and mortality rates (45). Chronic lung disease often causes respiratory distress and persistent coughing, resulting in a reduced quality of life (46). The association between cardiovascular disease and pulmonary disease has been supported by epidemiological (47) and physiologic evidence (48). A Canadian study concluded strong negative correlations between these diseases and the physical components of HRQoL (49). The coexistence of cardiovascular and pulmonary diseases can impair energy levels and social functioning, leading to diminished psychosocial adaptation and resilience (50).
With regard to the metabolic syndrome pattern, including dyslipidemia, diabetes, and hypertension, previous studies have demonstrated that hypertension was commonly combined with diabetes (9). Hypertension combined with diabetes and dyslipidemia was a relatively common pattern of co-morbidity in China, both in the elderly and in young population (51). Evidence from the China Health and Retirement Longitudinal Study demonstrated a higher co-morbidity prevalence of 10.00% in older population aged over 60 years (52). Uncontrolled blood pressure, glucose, and lipid levels heighten the risk of complications including cardiovascular and cerebrovascular diseases, as well as chronic renal insufficiency, which seriously threaten physical state, mental health, and even quality of life (53). Cai (51) et al. found that the 5-year survival rate of patients with hypertension and diabetes combined with dyslipidemia was lower than that of the general population.
In this study, we identified a less commonly reported patterns of multimorbidity termed digestive system pattern, which included ulcers of the digestive tract, gastroenteritis, neurasthenia, rheumatoid arthritis, etc. Ulcers of the digestive tract and gastroenteritis frequently occurred together (54). Ulcers of the digestive tract had the second highest prevalence (94.35%) among all multimorbidities, while gastroenteritis showed a multimorbidity prevalence of 85.44%, which might explain their frequent coexistence. In addition, neurasthenia and rheumatoid or arthritis were also included in the digestive system pattern. Alprazolam, commonly used to treat neurasthenia, can cause gastric cramps, nausea, and vomiting (55). Furthermore, the use of non-steroidal anti-inflammatory drugs (NSAIDs) in the management of rheumatoid arthritis may predispose patients to gastrointestinal ulcers and bleeding (56).
The hepatobiliary system pattern, which primarily consisted of gallstones/cholecystitis, chronic hepatitis/cirrhosis, bone fracture, and intervertebral disc disease, has been rarely reported in previous studies with multimorbidity patterns. Studies have shown that gallstones were associated with sarcopenia (57), while sarcopenia was positively associated with falls or facture (58). Acute cholecystitis was considered as a complication following hip fracture (59). Furthermore, individuals with chronic liver disease face have a heightened risk of developing osteoporosis, with fractures considered potential complications of this condition (60). The interaction of diseases between hepatobiliary system and fracture may partly explain the potential mechanism.
Although the Chinese government made an explicitly prioritization on development of rural healthcare and the launch of initiatives such as the New Rural Cooperative Medical System, rural healthcare continued to face serious challenges. Individuals living rural areas lack access to high-quality healthcare services, including accurate diagnosis, thorough physical examination, and timely treatment (61), which may result in lower self-awareness for health status. Primary healthcare is facilitated through public health and hospital services. Given that HRQoL is strongly influenced by physical wellbeing, psychological functioning, social capability, and overall personal status, these findings suggest that strategies aimed at improving quality of life outcomes are the most promising. The decline in HRQoL represents a critical clinical endpoint in the management of chronic diseases and multimorbidity. The results of the current study can help identify patients with multimorbidity for better resource allocation. Recognizing these patterns can be useful in disease categorization, management, reducing disease burden, and enhance patient health quality.
Therefore, we conclude that multimorbidity is even more critical than initially thought. The identification of multimorbidity patterns suggests that, in chronic disease management, policymakers and clinicians should not focus solely on individual chronic diseases. Instead, they should identify high-risk groups based on multimorbidity patterns and prioritize interventions accordingly. It is recommended to classify individuals with different chronic diseases for early assessment according to their risk of developing multimorbidity patterns and provide individualized management. For example, patients with rheumatoid arthritis should be regularly screened for pre-existing gastric or digestive diseases. In subgroup analyses, we observed that not all metabolic syndrome patterns exhibited statistically significant negative correlations with HRQoL, regardless of gender, age, or average yearly family total income. This finding underscores the complex nature of metabolic comorbidities. In addition, restricted cubic spline analysis among individuals aged 70–79 years showed that the digestive system, metabolic syndrome, and hepatobiliary system patterns each demonstrated non-linear “S” negative correlations with HRQoL. These results further emphasize the complexity of comorbid conditions in older adults.
The findings from the current research remain highly relevant, providing crucial insights for informing public health policies and guiding clinical interventions. This study has several notable strengths. First, it utilizes a large and representative sample drawn from the general population in southwest China. Second, the application of the EQ-5D-5L with the Chinese value set enhances the precision of our assessment, consistent with guidelines recommending country-specific evaluations. However, the potential limitations also should be acknowledged. First, we investigated only the prevalence of 13 common chronic diseases, and some diseases or conditions were underrepresented, potentially leading to an underestimation of the co-morbidity rate. Second, most of chronic conditions were self-reported by participants, which may introduce recall bias and potential misclassification errors as individuals may not accurately remember or describe their medical histories or conditions. Third, prevalent diseases or conditions with high prevalence, such as stroke, asthma, cancer, vision problems, hearing problems, and neurodegenerative diseases, were excluded. This exclusion might have resulted in underestimation attributable to recall bias or inaccuracies. Nonetheless, a systematic review reported a weak positive correlation between the number of chronic diseases included and multimorbidity, although this correlation was not statistically significant (p = 0.11) (62). Finally, detailed information on disease severity or duration time was not included in the current study.
Conclusion
The current study confirmed the prevalence and patterns of multimorbidity among residents in rural areas of southwest China and highlighted the importance of shifting from single disease treatment to multimorbidity prevention and management for rural adults in southwest China. Continued efforts to improve our understanding of multimorbidity patterns will facilitate the development of innovative primary care models and comprehensive home-based care models tailored to the specific needs of individuals with diverse multimorbidity patterns. The goal of improving health related quality of life can only be achieved by fully recognizing the heterogeneity of multimorbidity patterns. In other regions, effective prevention and control of chronic diseases require a comprehensive understanding of these patterns. Specifically, it is crucial to elucidate the distinct characteristics of multimorbidity across demographic groups, including gender and age. Based on this knowledge, targeted strategies for prevention, control, and treatment can be devised. Furthermore, rational allocation of health resources and greater emphasis on primary healthcare are essential to enhance the effectiveness of multimorbidity management.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The studies involving humans were approved by the ethics committee of Sichuan University (No. K2016038). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
NX: Software, Writing – original draft, Investigation, Data curation, Visualization, Validation, Formal analysis, Writing – review & editing. HX: Funding acquisition, Resources, Writing – review & editing, Formal analysis, Investigation, Data curation. XJ: Writing – review & editing, Validation, Supervision, Software, Formal analysis, Resources, Visualization.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the Health Commission of Chongqing and the Science and Technology Bureau of Chongqing (2024MSXM162) and the Chongqing Medical and Pharmaceutical Vocational Education Group (CQZJ202354).
Acknowledgments
It is a tremendous honor to be involved in the China Multi-Ethnic Cohort (CMEC) study. We express our sincere gratitude to all participants and team members for their invaluable assistance in completing this research project. We would like to acknowledge the language polishing service provided by Home for Researchers (www.home-for-researchers.com), which greatly enhanced the quality of our manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1609831/full#supplementary-material
References
1. Xu, X, Mishra, GD, and Jones, M. Evidence on multimorbidity from definition to intervention: an overview of systematic reviews. Ageing Res Rev. (2017) 37:53–68. doi: 10.1016/j.arr.2017.05.003
2. Zemedikun, DT, Gray, LJ, Khunti, K, Davies, MJ, and Dhalwani, NN. Patterns of multimorbidity in middle-aged and older adults: an analysis of the UK biobank data. Mayo Clin Proc. (2018) 93:857–66. doi: 10.1016/j.mayocp.2018.02.012
3. Yao, SS, Meng, X, Cao, GY, Huang, ZT, Chen, ZS, Han, L, et al. Associations between multimorbidity and physical performance in older Chinese adults. Int J Environ Res Public Health. (2020) 17:17. doi: 10.3390/ijerph17124546
4. Singh, K, Patel, SA, Biswas, S, Shivashankar, R, Kondal, D, Ajay, VS, et al. Multimorbidity in south Asian adults: prevalence, risk factors and mortality. J Public Health (Oxf). (2019) 41:80–9. doi: 10.1093/pubmed/fdy017
5. Guo, X, Zhao, B, Chen, T, Hao, B, Yang, T, and Xu, H. Multimorbidity in the elderly in China based on the China health and retirement longitudinal study. PLoS One. (2021) 16:e0255908. doi: 10.1371/journal.pone.0255908
6. Zhang, Q, Han, X, Zhao, X, and Wang, Y. Multimorbidity patterns and associated factors in older Chinese: results from the China health and retirement longitudinal study. BMC Geriatr. (2022) 22:470. doi: 10.1186/s12877-022-03154-9
7. Tran, TN, Lee, S, Oh, CM, and Cho, H. Multimorbidity patterns by health-related quality of life status in older adults: an association rules and network analysis utilizing the Korea National Health and nutrition examination survey. Epidemiol Health. (2022) 44:e2022113. doi: 10.4178/epih.e2022113
8. Zheng, DD, Loewenstein, DA, Christ, SL, Feaster, DJ, Lam, BL, McCollister, KE, et al. Multimorbidity patterns and their relationship to mortality in the US older adult population. PLoS One. (2021) 16:e0245053. doi: 10.1371/journal.pone.0245053
9. Sun, ZJ, Fan, JN, Yu, CQ, Guo, Y, Bian, Z, Pei, P, et al. Prevalence, patterns and long-term changes of multimorbidity in adults from 10 regions of China. Zhonghua Liu Xing Bing Xue Za Zhi. (2021) 42:755–62. doi: 10.3760/cma.j.cn112338-20200305-00259
10. Makovski, TT, Schmitz, S, Zeegers, MP, Stranges, S, and van den Akker, M. Multimorbidity and quality of life: systematic literature review and meta-analysis. Ageing Res Rev. (2019) 53:100903. doi: 10.1016/j.arr.2019.04.005
11. Fortin, M, Lapointe, L, Hudon, C, Vanasse, A, Ntetu, AL, and Maltais, D. Multimorbidity and quality of life in primary care: a systematic review. Health Qual Life Outcomes. (2004) 2:51. doi: 10.1186/1477-7525-2-51
12. Jackson, CA, Jones, M, Tooth, L, Mishra, GD, Byles, J, and Dobson, A. Multimorbidity patterns are differentially associated with functional ability and decline in a longitudinal cohort of older women. Age Ageing. (2015) 44:810–6. doi: 10.1093/ageing/afv095
13. Aoki, T, Yamamoto, Y, Shimizu, S, and Fukuhara, S. Physical multimorbidity patterns and depressive symptoms: a nationwide cross-sectional study in Japan. Fam Med Community Health. (2020) 8:e000234. doi: 10.1136/fmch-2019-000234
14. Aoki, T, Yamamoto, Y, Ikenoue, T, Onishi, Y, and Fukuhara, S. Multimorbidity patterns in relation to polypharmacy and dosage frequency: a nationwide, cross-sectional study in a Japanese population. Sci Rep. (2018) 8:3806. doi: 10.1038/s41598-018-21917-6
15. Gimeno-Miguel, A, Gracia Gutiérrez, A, Poblador-Plou, B, Coscollar-Santaliestra, C, Pérez-Calvo, JI, Divo, MJ, et al. Multimorbidity patterns in patients with heart failure: an observational Spanish study based on electronic health records. BMJ Open. (2019) 9:e033174. doi: 10.1136/bmjopen-2019-033174
16. Busija, L, Lim, K, Szoeke, C, Sanders, KM, and McCabe, MP. Do replicable profiles of multimorbidity exist? Systematic review and synthesis. Eur J Epidemiol. (2019) 34:1025–53. doi: 10.1007/s10654-019-00568-5
17. Jia, H, Lubetkin, EI, Barile, JP, Horner-Johnson, W, DeMichele, K, Stark, DS, et al. Quality-adjusted life years (QALY) for 15 chronic conditions and combinations of conditions among US adults aged 65 and older. Med Care. (2018) 56:740–6. doi: 10.1097/mlr.0000000000000943
18. Zhao, X, Hong, F, Yin, J, Tang, W, Zhang, G, Liang, X, et al. Cohort profile: the China multi-ethnic cohort (CMEC) study. Int J Epidemiol. (2021) 50:721–21l. doi: 10.1093/ije/dyaa185
19. Xiong, H, Zhao, P, Ma, F, Tang, D, and Liu, D. Spicy food intake and overweight/obesity in rural Southwest China: findings from a cross-sectional study. Front Nutr. (2025) 12:1526775. doi: 10.3389/fnut.2025.1526775
20. Zhu, K, Qian, F, Lu, Q, Li, R, Qiu, Z, Li, L, et al. Modifiable lifestyle factors, genetic risk, and incident peripheral artery disease among individuals with type 2 diabetes: a prospective study. Diabetes Care. (2024) 47:435–43. doi: 10.2337/dc23-1503
21. Xiong, Y, Zhang, F, Zhang, Y, Wang, W, Ran, Y, Wu, C, et al. Insights into modifiable risk factors of erectile dysfunction, a wide-angled Mendelian randomization study. J Adv Res. (2024) 58:149–61. doi: 10.1016/j.jare.2023.05.008
22. Zhang, YB, Chen, C, Pan, XF, Guo, J, Li, Y, Franco, OH, et al. Associations of healthy lifestyle and socioeconomic status with mortality and incident cardiovascular disease: two prospective cohort studies. BMJ. (2021) 373:n604. doi: 10.1136/bmj.n604
23. Ainsworth, BE, Haskell, WL, Whitt, MC, Irwin, ML, Swartz, AM, Strath, SJ, et al. Compendium of physical activities: an update of activity codes and MET intensities. Med Sci Sports Exerc. (2000) 32:S498–516. doi: 10.1097/00005768-200009001-00009
24. Chen, L, Tang, W, Wu, X, Zhang, R, Ding, R, Liu, X, et al. Eating spicy food, dietary approaches to stop hypertension (DASH) score, and their interaction on incident stroke in southwestern Chinese aged 30-79: a prospective cohort study. Nutrients. (2023) 15:1222. doi: 10.3390/nu15051222
25. Xiao, X, Qin, Z, Lv, X, Dai, Y, Ciren, Z, Yangla, Y, et al. Dietary patterns and cardiometabolic risks in diverse less-developed ethnic minority regions: results from the China multi-ethnic cohort (CMEC) study. Lancet Reg Health West Pac. (2021) 15:100252. doi: 10.1016/j.lanwpc.2021.100252
26. Fung, TT, Chiuve, SE, McCullough, ML, Rexrode, KM, Logroscino, G, and Hu, FB. Adherence to a DASH-style diet and risk of coronary heart disease and stroke in women. Arch Intern Med. (2008) 168:713–20. doi: 10.1001/archinte.168.7.713
27. Nagayoshi, M, Lutsey, PL, Benkeser, D, Wassel, CL, Folsom, AR, Shahar, E, et al. Association of sleep apnea and sleep duration with peripheral artery disease: the multi-ethnic study of atherosclerosis (MESA). Atherosclerosis. (2016) 251:467–75. doi: 10.1016/j.atherosclerosis.2016.06.040
28. Shi, BW, and Xiong, JY. Study on the impact of multiple chronic conditions on health related quality of life of Chinese elderly. Popul Develop. (2024) 30:120–8.
29. Williams, B, Mancia, G, Spiering, W, Agabiti Rosei, E, Azizi, M, Burnier, M, et al. 2018 ESC/ESH guidelines for the management of arterial hypertension. Eur Heart J. (2018) 39:3021–104. doi: 10.1093/eurheartj/ehy339
30. Xiong, H, Zhao, P, Ma, F, Tang, D, and Liu, D. The effects of interaction between smoking and abdominal obesity on pre-diabetes mellitus. Sci Rep. (2025) 15:1332. doi: 10.1038/s41598-025-85367-7
31. Xiong, H, Liu, D, Tang, D, and Ma, F. Female exposed to the Chinese famine increases the risk of dyslipidemia in later life. Medicine (Baltimore). (2023) 102:e34262. doi: 10.1097/md.0000000000034262
32. Gandhi, M, Ang, M, Teo, K, Wong, CW, Wei, YC, Tan, RL, et al. EQ-5D-5L is more responsive than EQ-5D-3L to treatment benefit of cataract surgery. Patient. (2019) 12:383–92. doi: 10.1007/s40271-018-00354-7
33. Xiong, H, Ma, F, Tang, D, and Liu, D. Correlations among nicotine dependence, health-related quality of life, and depression in current smokers: a cross-sectional study with a mediation model. Front Psych. (2024) 15:1455918. doi: 10.3389/fpsyt.2024.1455918
34. Luo, N, Liu, G, Li, M, Guan, H, Jin, X, and Rand-Hendriksen, K. Estimating an EQ-5D-5L value set for China. Value Health. (2017) 20:662–9. doi: 10.1016/j.jval.2016.11.016
35. Carlton, J, Peasgood, T, Khan, S, Barber, R, Bostock, J, and Keetharuth, AD. An emerging framework for fully incorporating public involvement (PI) into patient-reported outcome measures (PROMs). J Patient Rep Outcomes. (2020) 4:4. doi: 10.1186/s41687-019-0172-8
36. Hu, W, Zhou, L, Chu, J, Sun, N, Chen, X, Liu, S, et al. Estimating population norms for the health-related quality of life of adults in southern Jiangsu Province, China. Sci Rep. (2022) 12:9906. doi: 10.1038/s41598-022-13910-x
37. Grochtdreis, T, Dams, J, König, HH, and Konnopka, A. Health-related quality of life measured with the EQ-5D-5L: estimation of normative index values based on a representative German population sample and value set. Eur J Health Econ. (2019) 20:933–44. doi: 10.1007/s10198-019-01054-1
38. McCaffrey, N, Kaambwa, B, Currow, DC, and Ratcliffe, J. Health-related quality of life measured using the EQ-5D-5L: south Australian population norms. Health Qual Life Outcomes. (2016) 14:133. doi: 10.1186/s12955-016-0537-0
39. Yao, SS, Cao, GY, Han, L, Chen, ZS, Huang, ZT, Gong, P, et al. Prevalence and patterns of multimorbidity in a nationally representative sample of older Chinese: results from the China health and retirement longitudinal study. J Gerontol A Biol Sci Med Sci. (2020) 75:1974–80. doi: 10.1093/gerona/glz185
40. Palladino, R, Pennino, F, Finbarr, M, Millett, C, and Triassi, M. Multimorbidity and health outcomes in older adults in ten European health systems, 2006-15. Health Aff (Millwood). (2019) 38:613–23. doi: 10.1377/hlthaff.2018.05273
41. Mino-León, D, Reyes-Morales, H, Doubova, SV, Pérez-Cuevas, R, Giraldo-Rodríguez, L, and Agudelo-Botero, M. Multimorbidity patterns in older adults: an approach to the complex interrelationships among chronic diseases. Arch Med Res. (2017) 48:121–7. doi: 10.1016/j.arcmed.2017.03.001
42. Marengoni, A, Rizzuto, D, Wang, HX, Winblad, B, and Fratiglioni, L. Patterns of chronic multimorbidity in the elderly population. J Am Geriatr Soc. (2009) 57:225–30. doi: 10.1111/j.1532-5415.2008.02109.x
43. Aiwaiduli, S, Kadeer, G, Taxi, M, Yujuan, L, Nasier, K, and Alimujiang, S. Effects of chronic disease prevalence and comorbidity patterns on SRH status in middle-aged and elderly populations in rural areas. Chin Gen Pract. (2024) 27:3520–8.
44. Aoki, T, Fukuhara, S, Fujinuma, Y, and Yamamoto, Y. Effect of multimorbidity patterns on the decline in health-related quality of life: a nationwide prospective cohort study in Japan. BMJ Open. (2021) 11:e047812. doi: 10.1136/bmjopen-2020-047812
45. Wang, H, Zhang, H, and Zou, Z. Changing profiles of cardiovascular disease and risk factors in China: a secondary analysis for the global burden of disease study 2019. Chin Med J. (2023) 136:2431–41. doi: 10.1097/cm9.0000000000002741
46. Yee, LK, Guiqu, L, Juanjuan, H, Cheng, Z, and Xiangyan, Z. Path analysis of influencing factors of hospitalization expenses for chronic obstructive pulmonary disease patients. Modern Preventive Medicine (2022). 49:605–608+26.
47. Hole, DJ, Watt, GC, Davey-Smith, G, Hart, CL, Gillis, CR, and Hawthorne, VM. Impaired lung function and mortality risk in men and women: findings from the Renfrew and Paisley prospective population study. BMJ (1996). 313:711–715; discussion 15–6. doi: 10.1136/bmj.313.7059.711
48. Van Eeden, S, Leipsic, J, Paul Man, SF, and Sin, DD. The relationship between lung inflammation and cardiovascular disease. Am J Respir Crit Care Med. (2012) 186:11–6. doi: 10.1164/rccm.201203-0455PP
49. Fortin, M, Dubois, MF, Hudon, C, Soubhi, H, and Almirall, J. Multimorbidity and quality of life: a closer look. Health Qual Life Outcomes. (2007) 5:52. doi: 10.1186/1477-7525-5-52
50. Craig, LS, Cunningham-Myrie, CA, Theall, KP, Gustat, J, Hernandez, JH, and Hotchkiss, DR. Multimorbidity patterns and health-related quality of life in Jamaican adults: a cross sectional study exploring potential pathways. Front Med (Lausanne). (2023) 10:1094280. doi: 10.3389/fmed.2023.1094280
51. Guangyun, C, Yanjun, C, Libin, W, Deyi, L, shishi, W, Hong, W, et al. The 5-year survival of elderly hypertension and diabetes patients combined with dyslipidemia. Chin Gen Pract. (2024) 27:91–7.
52. Beijing Hypertension Association, China Association of Gerontology, Geriatrics Beijing Community Health Service, Association, Beijing Community Health Promotion Association, Chinese expert consensus on grassroots prevention and treatment of hypertension combined with type 2 diabetes mellitus and dyslipidemia in adults 2024. Chin Gen Pract. (2024). 27:3453–3475+82. doi: 10.12114/j.issn.1007-9572.2024.0116
53. Wang, LM, Chen, ZH, Zhang, M, Zhao, ZP, Huang, ZJ, Zhang, X, et al. Study of the prevalence and disease burden of chronic disease in the elderly in China. Zhonghua Liu Xing Bing Xue Za Zhi. (2019) 40:277–83. doi: 10.3760/cma.j.issn.0254-6450.2019.03.005
54. Butov, MA, Zhestkova, TV, Esakova, EM, and Efanova, LV. Assessment of ten-year dynamics of cases of hospitalizations of patients with peptic ulcer disease, chronic gastritis and chronic duodenitis. Ter Arkh. (2022) 94:160–5. doi: 10.26442/00403660.2022.02.201365
55. Türkoğlu, S. Paradoxical reactions related to alprazolam. J Child Adolesc Psychopharmacol. (2015) 25:276. doi: 10.1089/cap.2014.0116
56. Henry, D, Lim, LL, Garcia Rodriguez, LA, Perez Gutthann, S, Carson, JL, Griffin, M, et al. Variability in risk of gastrointestinal complications with individual non-steroidal anti-inflammatory drugs: results of a collaborative meta-analysis. BMJ. (1996) 312:1563–6. doi: 10.1136/bmj.312.7046.1563
57. Bhanji, RA, Narayanan, P, Allen, AM, Malhi, H, and Watt, KD. Sarcopenia in hiding: the risk and consequence of underestimating muscle dysfunction in nonalcoholic steatohepatitis. Hepatology. (2017) 66:2055–65. doi: 10.1002/hep.29420
58. Duque, G, Kirk, B, and Arai, H. Strong muscles lead to strong bones: the connection between osteoporosis, sarcopenia, falls and fractures. Bone. (2023) 173:116789. doi: 10.1016/j.bone.2023.116789
59. Yuan, Y, Tian, W, Jin, Z, Wang, L, and Zhu, S. Cholecystitis in older patients following hip fracture: a case series and literature review. BMC Geriatr. (2023) 23:694. doi: 10.1186/s12877-023-04336-9
60. Yang, YJ, and Kim, DJ. An overview of the molecular mechanisms contributing to musculoskeletal disorders in chronic liver disease: osteoporosis, sarcopenia, and osteoporotic sarcopenia. Int J Mol Sci. (2021) 22:22. doi: 10.3390/ijms22052604
61. Zhang, X, Dupre, ME, Qiu, L, Zhou, W, Zhao, Y, and Gu, D. Urban-rural differences in the association between access to healthcare and health outcomes among older adults in China. BMC Geriatr. (2017) 17:151. doi: 10.1186/s12877-017-0538-9
Keywords: prevalence, multimorbidity, multimorbidity pattern, health-related quality of life, correlation
Citation: Xie N, Xiong H and Jiang X (2025) Prevalence, patterns of multimorbidity, and its correlations with health-related quality of life in rural southwest China: a cross-sectional study. Front. Med. 12:1609831. doi: 10.3389/fmed.2025.1609831
Edited by:
Seyed Abolfazl Valizadeh, Shahid Beheshti University, IranReviewed by:
Marcus Cheetham, University Hospital Zurich, SwitzerlandRaheleh Davoodi, Shahid Beheshti University, Iran
Copyright © 2025 Xie, Xiong and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xin Jiang, MjE0MzEwNzU2QHFxLmNvbQ==
†These authors have contributed equally to this work and share first authorship