Skip to main content

Regional years of life lost, years lived with disability, and disability-adjusted life-years for severe mental disorders in Guangdong Province, China: a real-world longitudinal study

Abstract

Background

To understand the magnitude and spatial–temporal distribution of the regional burden attributable to severe mental disorders is of great essential and high policy relevance. The study aimed to address the burden of severe mental disorders by evaluating the years of life lost, years lived with disability, and disability-adjusted life-years (DALYs) in Guangdong, China.

Methods

We undertook a longitudinal study based on a multicenter database established by the Health Commission of Guangdong, involving a total of 21 prefectures and four economic regions in the Guangdong province. A total of 520,731 medical records from patients with severe mental disorders were collected for 2010–2020. Data were analyzed via an integrated evaluation framework by synthesizing prevalence estimates, epidemiological adjustment as well as comorbidity assessment to develop internally consistent estimates of DALY. DALY changes during 2010–2020 were decomposed by population growth and aging and further grouped by Socio-demographic Index (SDI). DALYs were projected to 2030 by the weighted median annualized rate of change in 2010–2020.

Results

In 2010–2020, the average DALYs for severe mental disorders reached 798,474 (95% uncertainty interval [UI]: 536,280–1,270,465) person-years (52.2% for males, and 47.8% for females). Severe mental disorders led to a great amount of disease burden, especially in Guangzhou, Shenzhen, and Foshan cities. Schizophrenia and mental retardation with mental disorders were the two leading sources of the burden ascribed to severe mental disorders. Population growth and aging could be accountable for the increasing burden of severe mental disorders. Economic regions with higher SDI carried a greater burden but had lower annualized rates of change in DALYs. The overall burden of severe mental disorders is projected to rise modestly over the next decade.

Conclusions

The findings urge prioritization of initiatives focused on public mental health, prevention strategies, health resources reallocation, and active involvement of authorities to effectively address the anticipated needs.

Background

Over the past three decades, China’s rapid economic growth was accompanied by an epidemiological transition from communicable to non-communicable diseases [1, 2]. Among the non-communicable diseases, mental disorders in China have historically been a low priority compared with cancers and cardiovascular diseases. Specifically, resource allocation, both financial and human, for mental disorders has been much scarce due to the inadequate understanding of public health policies and insufficient design of the healthcare delivery system [3, 4]. However, mental disorders are one of the primary sources of disease burden for population health [5, 6], leading to ongoing direct and indirect treatment costs [7, 8]. Previous research demonstrated that the lifetime prevalence of mental disorders in China may have reached 15.8–16.6% [9, 10], whereas the average treatment-seeking rate is only less than 6% [11, 12]. The severe gaps between detection and treatment of mental disorders may partly be ascribed to the incomplete knowledge about the magnitude of losses of healthy life-years in patients before resource allocation for diagnosis, treatment as well as management of these disorders [7]. Therefore, a better understanding of the magnitude of the disease burden attributed to mental disorders is crucial to prioritize the management of these conditions.

Findings from the Global Burden of Disease (GBD) study have resulted in a paradigm shift from determining the health loss due to mental disorders by prevalence rates, to conceptualizing such losses based on disability-adjusted life-years (DALYs) [7, 13]. The GBD 2017 reported the years of life lost (YLLs), years lived with disability (YLDs), and DALYs ascribed to mental health disorders in China’s 34 provinces as well as special administrative regions [2]. Mental health disorders ranked second place for non-fatal health outcomes measured by YLDs, and severe mental disorders such as schizophrenia are the leading causes of DALYs. Among the 34 provinces, Guangdong ranked first place for DALYs attributable to mental health disorders, with an estimated DALY of 1.84 (95% uncertainty interval [UI]: 1.38–2.37) million person-years by 2017. Given limited health resources as well as current management and intervention programs, the Chinese healthcare delivery system focused predominantly on severe mental disorders such as schizophrenia, schizoaffective disorder, mental retardation with mental disorders, and so forth [14,15,16]. Due to the rapid economic development, aging population as well as social-economic disparities in Guangdong [17, 18], there is a pressing need to quantify the regional spatial–temporal trends of DALYs for severe mental disorders to support evidence-based decision-making in public mental health. However, a systematic quantification of the magnitude as well as spatial–temporal distribution for regional burden attributable to severe mental disorders in Guangdong is not readily available.

Since 2010, Guangdong Health Commission has established the Guangdong Mental Health System (GDMHS), a population-based surveillance network for collecting the real-world medical records of such patients [18]. Data were gathered through efforts of consistent methodology including coherent diagnostic nomenclature, fully structured diagnostic interviews, and sophisticated household surveys [19, 20]. In this study, using the well-established GDMHS, we aim to estimate the regional disease burden associated with severe mental disorders in Guangdong and to present a detailed account of changes in the burden due to mental disorders, from 2010 to 2020. We also aim to project the estimated burden up to 2030, to inform prevention strategies as well as relevant policy-making to reduce the burden attributable to mental disorders.

Methods

Study design and population

This longitudinal study was conducted to evaluate the regional disease burden associated with severe mental disorders using mental health surveillance records from GDMHS. Within the GDMHS platform, the longitudinal data concerning information on diagnosis, medication, risk of violent behavior, etc., were collected and maintained in response to the National Health Commission of China [21], and it has extensive coverage of all of the 21 prefecture cities in Guangdong province. The province is located in southern China and adjacent to the Hong Kong and Macao Special Administrative Region. It remains the largest economic body in China and has a high population density of 115.21 million permanent residents [22]. According to the gross domestic product by prefecture, Guangzhou, Shenzhen, and Foshan are recognized as the most economically advantaged cities among the 21 prefectures. These prefectures were further categorized into four economic regions, namely, the Greater Bay Area, the North region, the West region, and the East region. The study complies with the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) statement.

Covariates

Variables available on the GDMHS platform included: (a) de-identified individual-level demographic information such as the residential address, age, biological gender, education level, marital and employment status; (b) a follow-up medical record of diagnosis for mental disorders and relevant comorbidities (shown in Additional file 1: Table S1); (c) family and genetic history of diagnosed mental disorders; and d) survival status (coded in the 10th revision of the International Classification of Diseases [23], ICD-10). Information on extracted variables was summarized in Table 1.

Table 1 Descriptive statistics of the collected study records (N = 520,731)

Data collection and processing

We collected valid registered medical records from GDMHS, covering January 1, 2010, to December 24, 2020. From the surveillance-based system, de-identified medical information was extracted on mental health disorders diagnosed by ICD-10 codes (F00-F99). Each record from GDMHS was validated and confirmed by the Health Commission of Guangdong upon data entry [2, 9]. The registered medical records consisted of 520,731 individuals that were clinically diagnosed with severe mental disorders including schizophrenia, schizoaffective disorder, paranoid psychosis, bipolar disorder, psychotic disorder due to epilepsy, and mental retardation with mental disorders [3, 24]. These disorders were anchored according to the lawful definition of severe mental disorders in the Chinese context [14, 25], and more details concerning information from GDMHS can be found in Additional file 1: Figure S1 and elsewhere [18]. Our collected dataset did not contain personally identifiable information and individuals could not be contacted, and thereby individual informed consent was not required. For data pre-processing, repeated records within the same year were treated as one single record, given the need for annual estimates of the burden metrics [26]. We collected city-specific population data and relevant information from the Guangdong Statistical Yearbook (2010–2020) [22] and the Global Health Data Exchange site of the Institute for Health Metrics and Evaluation. [27]

Data analysis

Burden of disease estimation

DALYs are calculated as the sum of YLL due to premature mortality and YLD due to non-fatal disability, thus incorporating both mortality and morbidity information. Using the DALY metric, the burden of disease across different cities with varying population structures can be directly compared [28,29,30]. We also presented age-standardized DALYs, which allowed comparisons over time and across different areas. The age-standardized rates were based on the Chinese population reported by the GBD 2019. [31]

As indicated by previous research [28, 32, 33], YLDs are the main contributors to DALYs attributed to mental disorders. In the study, underlying causes of death from the database did not register mental disorders as direct causes of death (Table 1). Thus, we only estimated prevalence from the collected records. Next, YLDs were estimated by multiplying the prevalence cases by disability weight (DW) [2, 33]. DW quantifies the magnitude of health loss associated with mental disorder on a scale from 0 (perfect health) to 1 (equivalent to death) and was extracted from the latest GBD (Additional file 1: Table S1) [29]. YLDs attributed to each mental disorder were then computed by applying the age-sex-city-year-specific cause distribution to the estimates [34]. For YLL, we used the indirect method by applying the estimated YLDs to the YLL/YLD ratios. This method was often used in the burden of disease study where information related to burden estimates was unavailable [35, 36]. The YLL/YLD ratios were summarized in Additional file 1: Table S2.

Prevalence estimates and epidemiological adjustment

To represent the continuous survival status of patients, and to better reflect the years lived with non-fatal mental disorders, incident counts of newly-diagnosed cases from the follow-up records were accumulated over years to generate crude prevalence counts by age, sex, and prefecture city. Deceased patients were not included in the estimation of crude prevalence. Age-sex-city-specific crude prevalence was calculated as the ratio between the annual accumulated number of diagnosed patients and the population in each city.

Next, to generate internally consistent epidemiological estimates, the age-sex-city-specific prevalence was adjusted using an epidemiological system dynamic model similar to the GBD method [29, 37]. A two-compartment framework with susceptible (S) and condition (C) compartments was implemented to account for the process of incidence, remission, and mortality. Based on the estimated number of patients with the condition and the number of susceptible populations, the adjusted age-sex-city-specific prevalence of mental disorders was calculated. Further information describing the epidemiological adjustment can be found in the appendix (Additional file 1: Table S12).

Comorbidity assessment

In line with the GBD study, we applied comorbidity assessment to the disability weight of mental disorders in the calculation of YLDs [29, 38]. Under the analysis framework, we performed comorbidity simulations separately for age-sex-city-specific estimates to account for the comorbidities amongst individuals. We synthesized a population of 20,000 individuals (\(n\)) and exposed the synthetic population to the binomial probability of having all the different comorbidities recorded in the GDMHS database. Using the synthesized population, a total of 11 comorbidities extracted from the follow-up records (Additional file 1: Table S1) were considered in the comorbidity assessment. Comorbidity-adjusted DW, as well as comorbidity-adjusted YLD rates attributed to severe mental disorders, were calculated, yielding the age-sex-city-year-specific YLD estimates. Next, we reported 95% uncertainty intervals (UI) by conducting a bootstrap analysis with 1000 repetitions, after which 95% uncertainty intervals were generated as the 2.5th and 97.5th values of the posterior distribution.

Decomposition and projection of DALY over time

To explore the underlying reasons for changes in disease burden over time, we implemented the time decomposition method proposed by Charlson [7] for DALY ascribed to mental health disorders. A total of two scenarios were set, with one to account for population growth and the other to explain the influence of population aging. Further details were summarized in the appendix. In the population growth scenario, the percentage changes between DALYs and scenario expected values reflect the change in absolute numbers due to population growth. The percentage change for the second scenario represents the change in absolute numbers due to changing age structure.

As a final step, we projected DALYs ascribed to mental disorders for each region of Guangdong to 2030, while implementing methods of the weighted median annualized rate of change [39, 40]. The projection process was summarized in the appendix methods and Additional file 1: Table S13. To be consistent with the GBD estimates, weighted median annualized rates of change in DALYs were categorized by the Socio-demographic Index (SDI), a composite indicator of development status ranging from 0 to 1 [31]. Details on SDI can be found elsewhere [41]. Five SDI quintiles (high, high-middle, middle, low-middle, and low SDI state groups) were empirically calculated. We conducted an exploratory analysis of the associations between the SDI and disease burden indicator, for logarithm and annualized rate of change of DALYs, based on Pearson correlation. Analyses were performed in R (version 3.6) and Python (version 3.8). A p-value of less than 0.05 was considered to be statistically significant.

Results

Table 1 shows the summary descriptive statistics in the study. From January 1, 2010, to December 24, 2020, records for a total of 520,731 individuals that were clinically diagnosed with severe mental disorders in Guangdong were extracted from GDMHS. The mean (SD) age at study enrollment was 39.7 (15.7) years, and most patients were aged 18 ~ 44 (54.8%). A total number of 229,448 (44.1%) individuals were females, while 291,283 (55.9%) were males. Sex-specific and residence-specific summary descriptive statistics were further depicted in Additional file 1: Table S3.

DALY of severe mental disorders

For overall DALYs of severe mental disorders, the average number of DALYs in 2010–2020 was estimated to be 798,474 (95% UI: 536,280, 1,270,465), of which 52.2% were estimated in males. We also estimated DALYs for mental disorders without comorbidity adjustment, which would shrink down to a number of 274,789 (95%UI: 191,823, 359,284) person-years. Table 2 demonstrates that the top three prefectures whose comorbidity-adjusted DALYs ranked the most were Guangzhou (183,090 with 95% UI: 122,811, 292,087), Shenzhen (150,599 with 95% UI: 100,871, 240,217), and Foshan (122,803 with 95% UI: 82,277, 195,698). Additional file 1: Table S4 presented the DALY estimates by urban and rural residents. Though the ranking of DALY burden by residence remains largely the same as the overall ranking, the number of DALY of severe mental disorders was significantly higher in rural residents across a few cities in the non-Greater Bay Area (namely, Zhanjiang, Qingyuan, Meizhou, Jieyang, Shaoguan, Maoming, Heyuan, and Yunfu), with non-overlapping 95% uncertainty intervals. Decomposition of DALYs into YLL and YLD by city, sex, and residence was shown in Additional file 1: Tables S5–S6.

Table 2 Average disability-adjusted life-years (DALYs) for patients with severe mental disorders in Guangdong, China, during 2010–2020, by sex and prefectural city

Trends and patterns from 2010 to 2020

Table 3 shows the changes in DALY absolute numbers and age-standardized rates between 2010 and 2020 for four economic regions of Guangdong. For the overall changes, these regions have seen a dramatic increase in the DALY burden from 2010 to 2020. However, increases in DALYs may be accounted for by different scenarios. For instance, even though Greater Bay Area remained the top region with the highest DALYs during 2010–2020, the proportions of change in absolute numbers that were due to population growth (39.96%) and population aging (39.87%) were relatively smaller than those in the non-Greater Bay Area (ranging from 65.75% ~ 76.32%). Besides, standardization of DALY rates across four regions reflected the deteriorating burden of disease over 2010–2020, regardless of age-sex-specific structures, particularly for the non-Greater Bay Area. Sex-specific and residence-specific scenario decomposition were shown in Additional file 1: Tables S7–S10.

Table 3 Change in absolute DALYa numbers and age-standardizedb rates (per 100,000) for patients with severe mental disorders in Guangdong, China, between 2010 and 2020, by economic region

Changing pattern of DALY burden was further illustrated in Fig. 1 (Fig. 1a for annualized rates of change by SDI quintile groups, and 1B for change of logarithm of DALYs). We found a modest negative correlation between SDI and annualized rates of change in DALYs ascribed to severe mental disorders across 21 prefectures and 4 economic regions of Guangdong (r = -0.49, p = 0.015). The East and West regions were similar in terms of SDI state groups. The logarithm of DALYs changed drastically for the non-Greater Bay Area from 2010 to 2020. Most prefectures in Greater Bay Area had higher SDIs, whereas prefecture cities in non-Greater Bay Areas clustered at low- to middle-SDI state groups. We also found a significant positive association of the estimated SDIs with the logarithm of DALYs across all prefectures (For 2010, r = 0.52, p = 0.009; for 2020, r = 0.62, p = 0.002). Regional and prefectural SDIs were presented in Additional file 1: Table S11. Similar patterns for the association of SDI with DALYs were observed for schizophrenia and mental retardation with mental disorders (Additional file 1: Figures S4–S5).

Fig. 1
figure 1

Weighted median annualized rates of change (a) and change of logarithm of values (b) in disability-adjusted life-years (DALYs) for patients with severe mental disorders in Guangdong, China, during 2010–2020, by SDI group. Notes: SDI = Socio-demographic Index. SDI groups were generated by calculating quintiles at the set of percentiles p = {0, 0.2, 0.4, 0.6, 0.8, 1.0}, in which five categories of SDI were yielded (namely, low SDI < 0.11, Low-middle SDI < 0.14, Middle SDI < 0.24, high-middle SDI < 0.52, and High SDI ≥ 0.56). In (a), the median of SDI during 2010–2020 was plotted, whereas the blue line reflects the linear relationship between estimated SDI and weighted median annualized rates of change in DALYs ascribed to severe mental disorders. In (b), the logarithm of DALYs across 21 prefecture cities was further aggregated by the median for the four economic regions. Arrows indicate the direction of change in the logarithm of DALY, where comorbidity-adjusted values were reported

Cause-specific averaged DALYs during 2010–2020

The averaged DALYs attributable to specific severe mental disorders were demonstrated in Additional file 1: Figure S3. The burden of severe mental disorders in the younger age groups under 20 years was primarily due to mental retardation with mental disorders, which accounted for around 17.0% (95% UI: 12.3, 24.4) of severe mental disorders for patients aged under 20 years. In contrast, the mental burden in patients in working-age groups (aged 20–69 years) was accounted for by schizophrenia, which took up around 12.4% (95% UI: 11.5, 13.7) for those aged 20–69 years. For all age groups combined, the absolute number and age-standardized rates of DALYs for schizophrenia were generally higher than those for mental retardation with mental disorders. The burden of schizophrenia cost an average number of 400,811 DALYs (95% UI: 271,132, 625,005) and 409.83 person-years per 100,000 population (95% UI: 272.40, 643.91), respectively, followed by mental retardation with mental disorders (with absolute number as well as age-standardized rates being 238,078 DALYs [95% UI: 157,498, 393844] and 211.39 person-years [95% UI: 139.21, 350.89] per 100,000, respectively).

Regional patterns through 2030

For the course of future burden in the next decade, the number of DALY attributed to severe mental disorders for the Greater Bay Area and North regions is expected to remain at a high level (Fig. 2), albeit the smaller annualized rates of change in DALY (0.035% and 0.031%, respectively for Greater Bay Area and North region). In the East and West regions, rising patterns were also observed in 2010–2020 but with slightly higher annualized rates of change (0.05% and 0.059%, respectively), leading to the increase of DALYs.

Fig. 2
figure 2

Temporal trends (a–c) for age-standardized DALYa rates and forecasts (d) for absolute DALYs ascribed to severe mental disorders in Guangdong, China, during 2010–2030, by economic region. aDALY: disability-adjusted life-years, for which comorbidity-adjusted estimates were reported. Age-standardized rates were calculated based on the Chinese population reported by the GBD 2019 [31]. Notes: In (a–c), region-specific patterns of age-standardized DALYs ascribed to severe mental disorders during 2010–2030 were depicted. The Greater Bay Area refers to Guangzhou(1), Foshan(2), Zhaoqing(5), Shenzhen(3), Dongguan(4), Huizhou(6), Zhuhai(7), Zhongshan(8), Jiangmen(9). The North Guangdong region refers to Shaoguan(10), Qingyuan(11), Meizhou(12), and Heyuan(13). The West Guangdong region refers to Zhanjiang(14), Yangjiang(15), Maoming(16) and Yunfu(17). Finally, the East Guangdong region refers to Chaozhou(18), Shantou(19), Jieyang(20), and Shanwei(21). In panel (D), temporal trends between 2010 and 2030 were shown, for four economic regions of Guangdong. Uncertainty intervals of forecasts for DALYs were presented in shaded colors

Discussion

Severe mental disorders have become a major public health concern in China. This study reports for the first time the current size and future burden of severe mental disorders for the Guangdong province, the most populous province and the key economic contributor to China. Schizophrenia and mental retardation with mental disorders accounted for a great proportion of severe mental disorders, which consume most health resources and bring great economic loss to patients and society [18]. Overall DALYs attributable to severe mental disorders in Guangdong may increase in the coming decades, where the DALY burden for the Greater Bay Area remains the heaviest. While the shift to a much greater number of DALYs currently occurred in the working-age groups, with an oversize population, aging structure, and extension of life expectancy in the next decades [2], policy-makers should remain aware of how aging shapes future health needs and understand the importance of prioritizing resources for intervention, treatment, rehabilitation, and prevention of severe mental disorders. Indeed, the number of DALYs represents the burden of disease that the health systems in a region must face. To the best of our knowledge, this is the first regional burden of disease study implementing the standard DALY metric while also incorporating the integrated information from a large-scale population-based medical system in China.

Our findings reiterate the complexity of management and intervention programs for patients with severe mental disorders because these patients may have mental disorder comorbidity [42]. As shown in Table 1, mental disorders not belonging to the six types of severe mental disorders were also registered in the GDMHS database, including behavioral syndromes, stress-related disorders, disorders of psychological development, and so forth. Comorbidities of these mental disorders nevertheless contributed to the overall disease burden for patients and the health services system as a whole [43]. Although mechanisms about how each mental disorder interacts with one another are beyond the scope of the current paper, implications drawn from our study are well applicable to the medical treatment and rehabilitation of patients with severe mental disorders. Due to the dynamics of synergistically related social-economic determinants, biological conditions, as well as complex conditions of mental health patients [44], health professionals may need to keep track of patients’ mental disorder comorbidities that are manifested either in sequence or simultaneously, while tailoring treatment strategies for precise medicine therapies. Moreover, the development of innovative medical treatment and the constant need to evaluate comorbidity dynamics over time in patients with severe mental disorders may be highly relevant to establishing clinical training programs for mental health professionals. It is worth noting that, there remains a pressing need for the cultivation of a greater number of mental health professionals, as supported by our foresight of a steady increase in the overall burden of severe mental disorders in the next decade. The aggregated burden of severe mental disorders over years will greatly increase the needs and demands for mental health services in the population [45]. However, as suggested by previous research [7, 46, 47], the shortage of human resources in public mental health for the Chinese population would not see improvement if training programs for mental health professionals and policy initiatives for mental health care remain unreformed. Active involvement of authorities in public mental health may also help remove barriers to accessing equitable and high-quality mental health services for the public [46]. Thereby, expanding the mental health workforce with qualified training programs, and higher involvement of authorities are the key priorities for public mental health in the coming decades in China.

DALYs ascribed to severe mental disorders that manifest predominantly during adulthood were generally higher in the high-SDI quintile group than in the low- and middle-SDI quintiles, which is consistent with GBD studies [13]. We found that the estimated SDI states were positively correlated with the DALY burden, while negatively correlated with the annualized rates of change in DALY. These relations reflect the association between education, income, as well as fertility level of a specific area, and the DALY indicator. As previously reported [18, 48, 49], lack of education and poverty are well-known determinants for mental disorders. Within the Greater Bay Area, Guangzhou, Shenzhen, and Foshan are more economically-empowered [22], the medical conditions and health services are therefore more readily available. As stated by Tan and colleagues [18], severe mental disorders in economically-developed cities are better managed and treated. We suspect that this may account for the lower annualized rates of change among some economically-empowered cities. However, due to the higher population density and an oversize population (illustrated in Additional file 1: Figure S1) [10, 50], the absolute numbers of DALYs in the high-SDI state group are notable, indicating the historically high burden in the Greater Bay Area.

As for the less developed regions, such as the East and West regions, even though the population density is lower, their smaller SDI quintiles indicate the poorer performance and progress in educational attainment and economic status. Thereby, we infer that the poorer SDI states might shed light on the unequal distribution of social and health resources for mental disorders. In agreement with a few studies [51, 52], unequal economic development and health services in poorer regions might account for the deteriorating burden of diseases and explain the more prominent annualized rates of change in regions with lower SDIs. The finding that rural residents in some non-Greater Bay cities during 2010–2020 carried a higher disease burden than their urban counterparts may further account for the pronounced overall change of severe mental burden in these less developed regions. Though the Chinese government is pursuing a rural industry re-vitalization strategy to promote economic development in rural areas and to close the urban–rural gap in the healthcare systems [53], previous publications [54, 55] found that the level of urban–rural coordination development in the non-Greater Bay Areas of Guangdong remained relatively low while there being much room to promote the performance of public health services in rural areas. Nevertheless, our study would highly appeal to the proper redistribution of limited resources among various regions, such that balanced allocation of necessary resources for the urban and rural areas can be coordinated; and also appeal to the more active involvement of political, economic, social, and health professional stakeholders at the local communities, especially in the relatively less developed regions and rural areas. [18, 56, 57]

Population growth and aging may have pivotal roles in the dynamic trends of DALY. As hinted by several GBD studies [7, 58], these factors would fundamentally alter the distribution of disease burden associated with mental disorders. In Guangdong, though the total population has been increasing over the last few decades, most prefecture cities also rapidly entered the aging societies [59]. Our study shows that the population growth and aging explained 39.96% and 39.87%, respectively, of the increase of severe mental burden in the Greater Bay Area. For the non-Greater Bay Areas and especially the rural areas, the contribution of population and aging seems even greater, with a proportion exceeding 60%. To explain the aging effects from the spatial distribution of patients, we found that the numbers of the elderly in the rural areas were higher than those in the urban areas (shown in Additional file 1: Table S3). Since the Greater Bay Area is well recognized as an economically-developed region [17, 18], we infer that it may attract more migrant workers while the more elderly people and those with mental health disorders are more likely to be left behind in the rural communities in the non-Greater Bay Areas. It may in part be held accountable for the higher population aging effect among the non-Greater Bay Areas, much similar to previous findings on the societal aging in rural China [60]. These findings constantly remind the authorities and relevant stakeholders of the importance of considering regional population structure and aging when planning for future health needs. Specifically, we anticipate the rising trends for overall DALYs in Guangdong to persist into the future, where the DALY burden for the Greater Bay Area remains the heaviest. The projected rising burden shall warrant urgent prioritization of mental health programs focused on population health, prevention strategies, and continuous medical support in the basic and specialized public health services programs tailored for patients with severe mental disorders. However, we noted that the estimates on the scenario decomposition by population growth and by changing age structure were much alike. The similarities may in part reflect the crucial role of the aging effect in exacerbating the mental health burden, while also suggesting the less prominent role of population growth due to birth and death. Given the more pressing issues of entering aging societies in China [59] and the rising future burden, it is high time to plan for mental public health services with a healthy aging agenda and to promote health at the later stage of life for patients with severe mental health disorders.

In comparison with the GBD study [2], our estimated numbers of DALYs for patients with severe mental disorders took up nearly half of the GBD estimates (43.4%), the latter of which reported DALY estimates ascribed to both severe and non-severe mental disorders. These findings not only confirm the importance of severe mental disorders but also suggest that schizophrenia may be one of the leading mental disorders and that mental burden is primarily due to the non-fatal burden of YLD [2, 28, 29, 32]. In contrast to the GBD 2017 [2], we observed slightly higher estimates for the age-standardized DALY rates of schizophrenia in Guangdong, and also found overlapping uncertainty intervals between our estimates and those reported by the GBD 2017 (307.33 person-years [95%UI: 231.98, 376.55] per 100,000). This is probably due to our comprehensive evaluation of the disease patterns and comorbidities presented by the real-world GDMHS database (shown in Additional file 1: Table S1) and it put out a crucial signal for drawing the attention of relevant stakeholders to care for the coexistent conditions and multi-system diseases in severe mental patients. The overlapping uncertainty intervals nonetheless reflect the consistency and robustness of these age-standardized results. Compared with other provinces and municipalities at a similar level of development (such as Beijing, Shanghai, and Tianjin), Guangdong has slightly higher age-standardized DALY rates of schizophrenia, albeit their overlapping uncertainty intervals. The similarity in the scales of these uncertainty intervals may add credibility to the findings of the present study.

The findings in our study strengthen the previous research in several ways. Firstly, a key strength of this study is the very large number of follow-up data collected over the last decade. The follow-up data is a population-based database and offers a great amount of validated information on the clinical as well as epidemiological features of patients with severe mental disorders [18]. With its extensive coverage of over 4900 mental health-related institutes, DALYs estimated from the database would reflect the urgency of the burden of severe mental disorders in a more thorough manner. Secondly, to our knowledge, this is the first and latest regional burden of disease study detailing the DALY burden for sub-regions of Guangdong. Next, the study enriched the burden of disease study by combining the GBD comorbidity assessment method with locally observed data. We found that the difference between comorbidity-adjusted DALYs and non-adjusted DALYs was almost as twice much as the non-adjusted burden estimates. This nonetheless highlights the importance of consideration of local disease patterns and comorbidities [61] when estimating the DALY burden.

However, some limitations should be noted in our study. Firstly, YLL could not be estimated directly from the database, since the underlying causes of death were primarily due to diseases other than mental disorders (Table 1). Secondly, we did not conduct a complete risk analysis, since it is beyond the scope of the current paper. However, we have utilized the novel SDI indicator [40] to interpret the changing patterns in DALYs. Annualized rates of change in DALYs categorized by SDI quintiles were similar to those reported in GBD. Finally, we could only estimate DALYs based on country-level disability weight from the GBD. Since no regional disability weights were publicly available, we argue future research in evaluating the regional disability weights may be needed. Notwithstanding, the proposed comprehensive evaluation framework and modeling schemes developed as part of the research can yet be applied to the other regional burden of disease study throughout different areas of the world.

Conclusions

This large-scale real-world follow-up study confirms the necessity of continuous support for prioritizing mental health services among patients with severe mental disorders. The study highlights the urgent need to redistribute health resources and prioritize mental health programs focused on population health and prevention strategies in places with low SDI states. With a drastic increase in population size, a transition towards aging societies, and an oversize population with severe mental disorders, the overall mental burden may become an enormous challenge for the authorities as well as relevant stakeholders. Our regional burden of disease study can be used by policy-makers and relevant stakeholders to plan for and improve the mental health strategies and services able to address the current as well as the future burden of mental disorders.

Availability of data and materials

The data that support the findings of this study are available from the corresponding authors, [HL and XL], upon reasonable request.

Abbreviations

GBD:

Global Burden of Disease

GATHER statement:

Guidelines for Accurate and Transparent Health Estimates Reporting statement

GDMHS:

Guangdong Mental Health System

YLL:

Year of life lost

YLD:

Year lived with disability

DALY:

Disability-adjusted life-year

SDI:

Socio-demographic Index

UI:

Uncertainty interval

ICD-10:

The 10th revision of the International Classification of Diseases

References

  1. Yang G, Wang Y, Zeng Y, et al. Rapid health transition in China, 1990–2010: findings from the Global Burden of Disease Study 2010. The Lancet. 2013;381(9882):1987–2015.

    Article  Google Scholar 

  2. Zhou M, Wang H, Zeng X, et al. Mortality, morbidity, and risk factors in China and its provinces, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. The Lancet. 2019.

  3. Liang D, Mays VM, Hwang W. Integrated mental health services in China: challenges and planning for the future. Health Policy Plann. 2018;33(1):107–22.

    Article  Google Scholar 

  4. Liu J, Ma H, He Y, et al. Mental health system in China: history, recent service reform and future challenges. World Psychiatry. 2011;10(3):210–6.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Qin X, Wang S, Hsieh C. The prevalence of depression and depressive symptoms among adults in China: estimation based on a National Household Survey. China Econ Rev. 2018;51:271–82.

    Article  Google Scholar 

  6. Naber D, Hansen K, Forray C, et al. Qualify: a randomized head-to-head study of aripiprazole once-monthly and paliperidone palmitate in the treatment of schizophrenia. Schizophr Res. 2015;168(1–2):498–504.

    Article  PubMed  Google Scholar 

  7. Charlson FJ, Baxter AJ, Cheng HG, Shidhaye R, Whiteford HA. The burden of mental, neurological, and substance use disorders in China and India: a systematic analysis of community representative epidemiological studies. The Lancet. 2016;388(10042):376–89.

    Article  Google Scholar 

  8. Yu Y, Peng M, Bai X, et al. Schizophrenia, social support, caregiving burden and household poverty in rural China. Soc Psych Psych Epid. 2020;55(12):1571–80.

    Article  Google Scholar 

  9. Huang Y, Wang Y, Wang H, et al. Prevalence of mental disorders in China: a cross-sectional epidemiological study. The Lancet Psychiatry. 2019;6(3):211–24.

    Article  PubMed  Google Scholar 

  10. Zhao Z, Huang Y, Li J, et al. An epidemiological survey of mental disorders in Guangzhou area. Chin J Nervous Ment Dis. 2009;35(9):530–4.

    CAS  Google Scholar 

  11. Demyttenaere K, Bruffaerts R, Posada-Villa J, et al. Prevalence, severity, and unmet need for treatment of mental disorders in the World Health Organization World Mental Health Surveys. J Am Med Assoc. 2004;291(21):2581–90.

    Article  Google Scholar 

  12. Phillips MR, Zhang J, Shi Q, et al. Prevalence, treatment, and associated disability of mental disorders in four provinces in China during 2001–05: an epidemiological survey. The Lancet. 2009;373(9680):2041–53.

    Article  Google Scholar 

  13. Sagar R, Dandona R, Gururaj G, et al. The burden of mental disorders across the states of India: the Global Burden of Disease Study 1990–2017. The Lancet Psychiatry. 2020;7(2):148–61.

    Article  Google Scholar 

  14. NADPC (National administration of disease prevention and control). An integrated group of joint prevention and control mechanisms for the outbreak of new coronavirus pneumonia of the state council of China. Notice on Strengthening the Treatment and Management of Patients with Severe Mental Disorders during the Outbreak of the New Coronary Pneumonia (in Chinese). http://www.nhc.gov.cn/xcs/zhengcwj/202002/f315a6bb2955474c8ca0b33b0c356a32.shtml. Accessed 12 March 2022.

  15. Guan L, Liu J, Wu XM, et al. Unlocking patients with mental disorders who were in restraints at home: a national follow-up study of China’s new public mental health initiatives. PLoS ONE. 2015;10(4):e121425.

    Google Scholar 

  16. Good BJ, Good MD. Significance of the 686 Program for China and for global mental health. Shanghai Arch Psychiatry. 2012;24(3):175–7.

    PubMed  PubMed Central  Google Scholar 

  17. Ho CY, Li D. Spatial dependence and divergence across Chinese cities. Rev Dev Econ. 2010;14(2):386–403.

    Article  Google Scholar 

  18. Tan W, Lin H, Lei B, et al. The psychosis analysis in real-world on a cohort of large-scale patients with schizophrenia. BMC Med Inform Decis. 2020;20(3):132.

    Article  CAS  Google Scholar 

  19. Liu Z, Huang Y, Lv P, et al. The China mental health survey: II. Design and field procedures. Soc Psych Psych Epid. 2016;51(11):1547–57.

    Article  Google Scholar 

  20. Lu L. Review and foresight of the development of mental health science in China. J Peking Univ (Health Sci). 2019;51(3):379–83.

    Google Scholar 

  21. NHC (The National Health Commission of the People's Republic of China). Notice of the Health Commission on Printing and Distributing the Practices of Treatment and Management of Severe Mental illness (2018 Edition). http://www.gov.cn/gongbao/content/2018/content_5338247.htm. Accessed 12 April 2022 (in Chinese).

  22. GDSTATS (Guangdong Provincial Bureau of Statistics). Guangdong Statistical Yearbook (2010–2020). http://stats.gd.gov.cn/gdtjnj/. Accessed 12 March 2022.

  23. WHO (World Health Organization). International Classification of Diseases, Tenth Revision (ICD-10). 2015.

  24. Liu Y, Liu X, Wen H, et al. Risk behavior in patients with severe mental disorders: a prospective study of 121,830 patients managed in rural households of western China. BMC Psychiatry. 2018;18(1):134.

    Article  PubMed  PubMed Central  Google Scholar 

  25. SCPRC (State Council of the People's Republic of China). Release of Management for Reporting the Incidence of Severe Mental Disorders (Trial). http://www.gov.cn/jrzg/2013-09/02/content_2479713.htm. Accessed 12 March 2022.

  26. Lin X, Du Z, Liu Y, Hao Y. The short-term association of ambient fine particulate air pollution with hypertension clinic visits: A multi-community study in Guangzhou, China. Sci Total Environ. 2021;774:145707.

    Article  CAS  PubMed  Google Scholar 

  27. IHME (Institute for Health Metrics and Evaluation). Global Health Data Exchange (GHDx). http://ghdx.healthdata.org/. Accessed 12 March 2022.

  28. van Hoeken D, Hoek HW. Review of the burden of eating disorders: mortality, disability, costs, quality of life, and family burden. Curr Opin Psychiatry. 2020;33(6):521–7.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Vos T, Lim SS, Abbafati C, et al. Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. The Lancet. 2020;396(10258):1204–22.

    Article  Google Scholar 

  30. Bemme D, Souza NA. Global mental health and its discontents: An inquiry into the making of global and local scale. Transcult Psychiatry. 2014;51(6):850–74.

    Article  PubMed  Google Scholar 

  31. Murray CJL, Aravkin AY, Zheng P, et al. Global burden of 87 risk factors in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. The Lancet. 2020;396(10258):1223–49.

    Article  Google Scholar 

  32. Bonadiman CSC, Malta DC, de Azeredo Passos VM, Naghavi M, Melo APS. Depressive disorders in Brazil: results from the Global Burden of Disease Study 2017. Popul Health Metr. 2020;18(Suppl 1):6.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Phanthunane P, Vos T, Whiteford H, Bertram M, Udomratn P. Schizophrenia in Thailand: prevalence and burden of disease. Popul Health Metr. 2010;8(1):24.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Devleesschauwer B, Haagsma JA, Angulo FJ, et al. Methodological framework for World Health Organization estimates of the global burden of foodborne disease. PLoS ONE. 2015;10(12):e142498.

    Article  CAS  Google Scholar 

  35. Fernández de Larrea-Baz N, Álvarez-Martín E, Morant-Ginestar C, et al. Burden of disease due to cancer in Spain. BMC Public Health. 2009;9(1):1–11.

    Article  Google Scholar 

  36. Prayitno L, Sugiharto M, Pratiwi NL. Ischemic stroke caused by metabolite factors in Indonesia (GBD study 2017). Indian J Forensic Med Toxicol. 2020;14(4).

  37. Flaxman AD, Vos DT, Murray CJ. An integrative meta-regression framework for descriptive epidemiology. Washington: University of Washington Press; 2015.

    Google Scholar 

  38. Vos T, Flaxman AD, Naghavi M, et al. Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. The Lancet. 2012;380(9859):2163–96.

    Article  Google Scholar 

  39. Fullman N, Barber RM, Abajobir AA, et al. Measuring progress and projecting attainment on the basis of past trends of the health-related Sustainable Development Goals in 188 countries: an analysis from the Global Burden of Disease Study 2016. The Lancet. 2017;390(10100):1423–59.

    Article  Google Scholar 

  40. Lim SS, Fullman N, Murray CJ, Mason-Jones AJ. Measuring the health-related Sustainable Development Goals in 188 countries: a baseline analysis from the Global Burden of Disease Study. The Lancet. 2015;2016:1–38.

    Google Scholar 

  41. Vos T, Abajobir AA, Abate KH, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. The Lancet. 2017;390(10100):1211–59.

    Article  Google Scholar 

  42. Scott KM, McGee MA, Oakley Browne MA, Wells JE. Mental disorder comorbidity in Te Rau Hinengaro: The New Zealand Mental Health Survey. Aust N Z J Psychiatry. 2006;40(10):875–81.

    Article  PubMed  Google Scholar 

  43. Andrews G, Henderson S, Hall W. Prevalence, comorbidity, disability and service utilisation: overview of the Australian National Mental Health Survey. Br J Psychiatry. 2001;178(2):145–53.

    Article  CAS  PubMed  Google Scholar 

  44. Tsai AC. Syndemics: A theory in search of data or data in search of a theory? Soc Sci Med. 2018;206:117–22.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Fang M, Hu SX, Hall BJ. A mental health workforce crisis in China: A pre-existing treatment gap coping with the COVID-19 pandemic challenges. Asian J Psychiatr. 2020;54:102265.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Xiang Y, Ng CH, Yu X, Wang G. Rethinking progress and challenges of mental health care in China. World Psychiatry. 2018;17(2):231–2.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Hu X, Rohrbaugh R, Deng Q, He Q, Munger KF, Liu Z. Expanding the mental health workforce in China: narrowing the mental health service gap. Psychiat Serv. 2017;68(10):987–9.

    Article  Google Scholar 

  48. Herman H, Jané-Llopis E. Mental health promotion in public health. Promot Educ. 2005;12(2_suppl):42–7.

    Article  Google Scholar 

  49. Davis J, Eyre H, Jacka FN, et al. A review of vulnerability and risks for schizophrenia: beyond the two hit hypothesis. Neurosci Biobehav Rev. 2016;65:185–94.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Long J, Huang G, Liang W, et al. The prevalence of schizophrenia in mainland China: evidence from epidemiological surveys. Acta Psych Scand. 2014;130(4):244–56.

    Article  CAS  Google Scholar 

  51. Wong DFK, Zhuang XY, Pan JY, He XS. A critical review of mental health and mental health-related policies in China: more actions required. Int J Soc Welf. 2014;23(2):195–204.

    Article  Google Scholar 

  52. Fan VS, Mahadevan R, Leung J. Effect of income inequality, community infrastructure and individual stressors on adult depression. Health Promot Int. 2021;36(1):46–57.

    Article  PubMed  Google Scholar 

  53. Liu Y, Zang Y, Yang Y. China’s rural revitalization and development: theory, technology and management. J Geogr Sci. 2020;30(12):1923–42.

    Article  Google Scholar 

  54. Tang D, Li B, Qiu Y, Zhao L. Research on urban and rural coordination development and its driving force based on the space-time evolvement taking Guangdong Province as an example. Land. 2020;9(8):253.

    Article  Google Scholar 

  55. Tian M, Feng D, Chen X, et al. China’s rural public health system performance: a cross-sectional study. PLoS ONE. 2013;8(12):e83822.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Gong W, Zhang C, Xu DR, et al. The association between a free medicine program and functioning in people with schizophrenia: a cross-sectional study in Liuyang, China. PeerJ. 2020;8:e8929.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Cheng HG, Shidhaye R, Charlson F, et al. Social correlates of mental, neurological, and substance use disorders in China and India: a review. The Lancet Psychiatry. 2016;3(9):882–99.

    Article  PubMed  Google Scholar 

  58. Murray CJ, Lopez AD. The global burden of disease: a comprehensive assessment of mortality and disability from diseases, injuries, and risk factors in 1990 and projected to 2020: summary. World Health Organization, 1996.

  59. Feng J, Hong G, Qian W, Hu R, Shi G. Aging in China: an international and domestic comparative study. Sustainability-Basel. 2020;12(12):5086.

    Article  Google Scholar 

  60. Fang Y, Wang H, Ouyang Z. An overview report on the situation of elderly people in urban and rural China: 2000–2015. China Popul Dev Stud. 2019;2(3):323–45.

    Article  Google Scholar 

  61. Zhou Z, Wang C, Yang H, Wang X, Zheng C, Wang J. Health-related quality of life and preferred health-seeking institutions among rural elderly individuals with and without chronic conditions: a population-based study in Guangdong Province, China. BioMed Res Int. 2014.

Download references

Acknowledgements

We thank the individuals and workers at the multiple CHSCs across different cities in Guangdong. We thank the editor and reviewers for their professional suggestions which greatly improve the manuscript. Professor Yuantao Hao gratefully acknowledges the support of K.C.Wong Education Foundation.

Funding

This work was supported in part by the Guangdong Basic and Applied Basic Research Foundation [Grant Numbers 2020A1515011294, 2020A1515110230, 2021A1515011765, 2021A1515011591], and by the China Postdoctoral Science Foundation [Grant Number 2021M693594]. Professor Yuantao Hao gratefully acknowledges the support of K.C.Wong Education Foundation. The sponsor of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

XL, LC, and WT conceived and designed the research. WT and LC contributed to the acquisition of data. XL and LC contributed to data curation and its analysis. XL, LC, and WT wrote the paper. YZ, JX, YH, FJ, BH, JG, and SW contributed to the visualization and investigation of the manuscript. XL and HL is also the guarantor of the manuscript. All co-authors contributed to revisions of the manuscript, and have reviewed and approved the final version of the manuscript for publication. The corresponding author had full access to all of the data in the study and the final responsibility for publication. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Haicheng Lin or Xiao Lin.

Ethics declarations

Ethics approval and consent to participate

The burden of disease study was approved by the Human Ethics Committee at Sun Yat-sen University (SYSU-SPH-2020126).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Supplementary Information

Additional file 1

. Section 1. Methods in details. Section 2. Table S1. Comorbidity list recorded by the GDMHS with disability weights extracted from the Global Burden of Disease (GBD) studies. Table S2. Ratios of years of life lost (YLL) to years lived with disability (YLD) for mental disorders, extracted from the Global Burden of Disease (GBD) studies, for 2010-2020. Table S3. Descriptive statistics of the collected study records (N = 520731), by sex and residence. Table S4. Average disability-adjusted life-years (DALYs) for patients with severe mental disorders in Guangdong, China, during 2010-2020, by residence and prefectural city. Table S5. Average years of life lost (YLLs) and years lived with disability (YLDs) for severe mental disorders in Guangdong province during 2010 ~ 2020, by sex and prefectural city. Table S6. Average years of life lost (YLLs) and years lived with disability (YLDs) for severe mental disorders in Guangdong province during 2010 ~ 2020, by residence and prefectural city. Table S7. Change in absolute DALY numbers and age-standardized rates (per 100000) for severe mental disorders in Guangdong province, China, between 2010 and 2020, by economic region, for females. Table S8. Change in absolute DALY numbers and age-standardized rates (per 100000) for severe mental disorders in Guangdong province, China, between 2010 and 2020, by economic region, for males. Table S9. Change in absolute DALY numbers and age-standardized rates (per 100000) for severe mental disorders in Guangdong province, China, between 2010 and 2020, by economic region, for the urban residents. Table S10. Change in absolute DALY numbers and age-standardized rates (per 100000) for severe mental disorders in Guangdong province, China, between 2010 and 2020, by economic region, for the rural residents. Table S11. Socio-demographic Index (SDI) for each prefecture and economic region in Guangdong province, China, between 2010 and 2020. Table S12. Evaluation metrics for different generalized additive model (GAM) based age smoothing functions, implemented in the epidemiological adjustment for prevalence estimates. Table S13. A range of ω used in the weighted median annualized rate of change formula, and the corresponding RMSE metrics for each city. Figure S1. Location of each local community health service center (CHSC) and psychiatric hospital within the Guangdong Mental Health center network medical System (GDMHS), in Guangdong, China, 2020. Figure S2. The estimated age-specific duration for mental disorders in Guangdong province, by sex, and by city. Figure S3. The average number of disorder-specific DALYs (person-years) for each economic region in Guangdong province, China, during 2010 ~ 2020, by age. Figure S4. Weighted median annualized rates of change (A) and change of logarithm of values (B) in disabilityadjusted life-years (DALYs) for patients with schizophrenia in Guangdong province, China, from 2010 to 2020, by SDI group. Figure S5. Weighted median annualized rates of change (A) and change of logarithm of values (B) in disability-adjusted life-years (DALYs) for patients with severe mental retardation in Guangdong province, China, from 2010 to 2020, by SDI group.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tan, W., Chen, L., Zhang, Y. et al. Regional years of life lost, years lived with disability, and disability-adjusted life-years for severe mental disorders in Guangdong Province, China: a real-world longitudinal study. glob health res policy 7, 17 (2022). https://doi.org/10.1186/s41256-022-00253-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41256-022-00253-3

Keywords