The global burden and attributable risk factors of chronic lymphocytic leukemia in 204 countries and territories from 1990 to 2019: analysis based on the global burden of disease study 2019

<jats:title>Abstract</jats:title><jats:sec>                 <jats:title>Background</jats:title>                 <jats:p>Chronic lymphocytic leukemia (CLL) is the most prevalent subtype of leukemia in Western countries, causing a substantial health burden on patients and society. Comprehensive evaluation of the epidemiological characteristics of CLL is warranted, especially in the current context of global population aging. The main objective of this study is evaluating the disease burden of CLL at global, regional, and national levels from 1990 to 2019. As secondary objectives, we studied the influence of demographic factors and performed risk factor analysis. We hope this study could provide evidence for the evaluation of the effectiveness of previous prevention strategies and the formulation of future global health policies.</jats:p>               </jats:sec><jats:sec>                 <jats:title>Results</jats:title>                 <jats:p>Based on data of CLL between 1990 to 2019 from the Global Burden of Disease (GBD) study 2019, we depicted the age, gender, and regional structure of the CLL burden population and described the impact of social development on the disease burden of CLL. The distribution and changing trends of attributable risk factors were also investigated. The global burden of CLL has increased dramatically. A high incidence has been achieved in males and elder people. Countries and territories with high social-demographic index (SDI) tended to have higher global burden than low-SDI region. Of risk factors, high body mass index and smoking were the major contributors for CLL-related mortality and disability adjusted life-years (DALYs).</jats:p>               </jats:sec><jats:sec>                 <jats:title>Conclusion</jats:title>                 <jats:p>In summary, the global CLL burden continues to rise over the past 30 years. The relocation of medical resource should be considered on a global scale.</jats:p>               </jats:sec><jats:sec>                 <jats:title>Graphical Abstract</jats:title>                                </jats:sec>

section_name: Graphical Abstract

section_content: Chronic lymphocytic leukemia (CLL) represents a prevalent adult leukemia, which is characterized by abnormal accumulation of immunologically incompetent lymphocytes in blood, bone marrow, lymph nodes, and spleen. CLL accounts for 25-30% of all the leukemia in Western Countries [1], with over 100,000 incidence cases and over 40,000 death cases globally reported in 2019. Epidemiological studies found that the incidence of CLL rises exponentially with age and reaches a peak in elderly populations [2]. The incidence of CLL is approximately 2 times higher in males than that in females [3, 4]. Additionally, markable geographical imbalances were found in CLL-related incidence cases. While CLL is the most prevalent adult leukemia in Western countries, it is relatively rare in Asia, even in Asian immigrants moving to the Western hemisphere [4] [5] [6]. Another specific depiction of CLL epidemiology comes from the Arab world. It is reported that CLL incidence of the Jews in Israel is significantly higher than the Arabs in Israel and the Arabs in the surrounding Middle Eastern countries, which suggested that ethnic factors rather than geographical factors may be a more critical contributor [7]. 

Despite of promising results in emerging targeted medications including BCL-2 inhibitor venetoclax and Bruton tyrosine kinase (BTK) inhibitors represented by ibrutinib and zanubrutinib [8] [9] [10] [11], it cannot be neglected that the high-cost treatment and accompanied severe adverse events contributed to a heavy global burden to CLL patients. By far CLL's incidence and mortality are still increasing both in developing and developed countries, with a heterogeneous survival rate correlated with local medical conditions and economic settings [12]. But up to date, previous studies on CLL burden presented several limitations. Some studies only described the epidemiological profile of CLL. The lack of adequate indicators to reflect CLL's disease burden has weakened the intuitive understanding of the hazards of disease burden [13]. Besides, most studies focused on the disease burden of CLL in a localized region, within most underdeveloped areas ignored [14, 15]. In addition, significant differences of disease burden patterns were existed in heterogeneous types of leukemias. Previously, most studies using the Global Burden of Disease (GBD) platform focused on the entire leukemia category, leading to a neglect of the specific internal genetic differences existing across leukemia [16, 17]. Thus, there is no guarantee that the evidence offered by these studies could provide accurate guidance on targeting CLL burden management. In addition to the existence of the above-mentioned scientific gaps, we also considered that apart from genetic factors, social development and medical advance have also caused significant changes in the disease burden patterns of CLL. Therefore, an updated and comprehensive picture of global CLL-contributed disease burden is warranted for assessing the status of public health and formulating future policies of medical resource allocation. 

The GBD study 2019 assessed epidemiologic data about 369 diseases across 204 countries and territories and provided an unprecedented opportunity to understand the trends in the global burden of CLL [18, 19]. In this study, we collected disability adjusted life-years (DALYs) from the GBD study 2019 to intuitively reflect the burden of disease and then calculated estimated annual percentage changes (EAPCs) to evaluate the disease burden changing trends. Based on the depict of overall trend of the disease burden caused by CLL from 1990 to 2019, we investigated the impact of age, gender, region and socio-demographic development level on the CLL burden. Further we performed a risk factor analysis on CLL-contributed disease burden to clarify the distribution and changing trend of risk factors. 

To our best knowledge, this study is the first study to provide a comprehensive description of the epidemiology and global burden of CLL worldwide. Moreover, this study also extended the following hypotheses: (1) The higher the proportion of men in the population structure, the greater the disease burden of CLL; (2) The more serious the degree of aging, the greater the disease burden of CLL; (3) The higher the degree of social development, which could be measured by Social-demographic index (SDI), the greater the disease burden of CLL. Up to date, since the health care policies, epidemiology and disease burden trends of CLL have changed over time, understanding the updated CLLcontributed disease burden is critical to assess the effectiveness of previous prevention strategies. Besides, this study can also provide certain instructive evidence for the design of the follow-up health policies. This study described the age, gender, and regional structure of CLL-burdened population. Policymakers should consider the demographic structure of local CLL-burdened populations and social development level to adjust local medical care policies. In addition, the lifestyle and occupational environment of patients also have impact on their survival. Finally, knowledge of the attributable risk factors of CLL provides theoretical basis and specific directions for CLL prevention and management.
section_name: Background

section_name: Results

section_name: CLL-related incidence, death, DALY burden and corresponding change trends

section_content: Globally, during the last 30 years, CLL-related incidence cases increased significantly from 40,537 in 1990 to 103,467 in 2019, with age-standardized incidence rate (ASIR) rising from 0. 76/100,000 persons in 1990 to 1. 34/100,000 persons in 2019 (EAPC: 1. 86, 95% CI: 1. 79-1. 92) (Table 1 ). Based on SDI-stratified regional analysis, the number of incidence cases and respective ASIR increased in all SDI categories between 1990 and 2019, with high-SDI quintiles exhibiting the highest ASIR in 2019. Of note, the most rapid increase was observed in middle-SDI quintiles. 

In the geographical region levels, Western Europe, High-income North America, and Central Europe displayed the highest ASIR in 2019, while East Asia, Central Europe, and Andean Latin America showed rapidest growth. In the country or territory level, of 204 countries and territories, the USA, China, and India were the 3 countries with the highest incidence cases of CLL in 2019 (Additional file 1: Table S1 ) (Fig. 1a ). Croatia, Monaco, and Slovenia displayed the highest ASIR in 2019 (Additional file 1: Table S2 ) (Fig. 2a ).
section_name: CLL-related incidence and its change trends

section_content: Global deaths cases of CLL had a prompt growth from 21,548 in 1990 to 44,613 in 2019, with age-standardized death rate (ASDR) rising from 0. 40/100,000 persons in 1990 to 0. 58/100,000 persons in 2019 (EAPC: 1. 17, 95% CI: 1. 07-1. 27) (Table 2 ). 

Based on SDI-stratified regional analysis, the number of death cases and corresponding ASDR increased in all SDI categories between 1990 and 2019, with high-SDI quintiles exhibiting the highest ASDR in 2019. Of note, the promptest death was observed in middle-SDI quintiles, consistent with incidence trends. 

In the geographical region levels, the highest ASDR was found in Central Europe, Western Europe, and High-income North America in 2019. East Asia, Central Europe, and Andean Latin America showed the rapidest growth. In the country or territory level, India, the USA, and China were the 3 countries with the highest death cases of CLL in 2019 (Additional file 1: Table S3 ) (Fig. 1b ). Croatia, Latvia, and Lithuania displayed the highest ASDR in 2019 (Additional file 1: Table S4 ) (Fig. 2b ).
section_name: CLL-related death and its change trends

section_content: Global DALY cases of CLL increased rapidly from 492,075 in 1990 to 948,464 in 2019, with age-standardized DALY rate rising from 9. 20/100,000 persons in 1990 to 12. 26/100,000 persons in 2019 (EAPC:0. 92, 95% CI: 0. 90-0. 94) (Table 3 ). 

CLL, chronic lymphocytic leukemia; ASIR, age-standardized incidence rate; EAPC, estimated annual percentage changes; SDI, social-demographic index According to SDI-stratified regional analysis, the number of DALY cases and corresponding age-standardized DALY rate increased in all SDI stratifications from 1990 to 2019. Among all SDI stratifications, high-SDI quintiles had the highest age-standardized DALY rate. Of note, the promptest DALY was observed in middle-SDI quintiles, consistent with incidence, also consistent with incidence trends. 

In the geographical region levels, Central Europe, Western Europe, and Eastern Europe exhibited the highest age-standardized DALY rate in 2019, while East Asia, Central Europe, and Andean Latin America showed the rapidest growth. In the country or territory level, China, India, and the USA were the 3 countries burdened with the highest number of CLL DALY cases in 2019 (Additional file 1: Table S5 ) (Fig. 1c ). Latvia, Croatia, and Poland displayed the highest age-standardized DALY rate in 2019 (Additional file 1: Table S6 ) (Fig. 2c ).
section_name: CLL-related DALY burden and its change trends

section_content: The incidence, death, and DALY burden differ by sex. Males were more likely to suffer from CLL than females (male: female in ASIR = 1. 13:1 in 1990, and 1. 14:1 in 2019) (Fig. 3a ). Whereas, the death gap between males and females continuously shrank and even reversed in the past 30 years (male: female in ASIR = 1. 05:1 in 1990, and 0. 98:1 in 2019) (Fig. 3b ). DALY burden of CLL was heavier in males than females (male: female in age-standardized DALY rate = 1. 12:1 in 1990, and 1. 04:1 in 2019), with an ongoing narrowing gap (Fig. 3c ). 

We then evaluated CLL incidence and ASIR in 3 different age groups: 15-49 years, 50-69 years, and above 70 years in the globe and different regions based on SDI levels. The results revealed that most incidences occurred in the population with 50 years of age or older. Furthermore, in the high-SDI region, the incidence cases of patients above 70 years occupied the most proportion. While in low-SDI region, the incidence cases of patients aged 50-69 years accounted for the highest percentage (Fig. 4a ). In all age groups, patients aged above 70 years displayed the highest incidence rate, especially in the high-SDI region (Fig. 4b ).
section_name: Sex and age patterns of CLL

section_content: We evaluated the relationship between ASIR of CLL in 1990 and corresponding EAPC and found that that the EAPC of ASIR was negatively correlated with ASIR (correlation coefficient = -0. 19, P = 0. 0058), indicating that the CLL incidence of countries and territories with low ASIR could be substantially underestimated (Fig. 5 ). Then, we investigated the correlation between SDI and EAPC values of ASIR, ASDR, and agestandardized DALY rate in 21 geographical regions across the globe. All age standardized ratios (ASRs) values displayed an apparent positive correlation with SDI (correlation coefficient of ASIR = 0. 70, of ASDR = 0. 68, of age-standardized DALY rate = 0. 67, all P values < 0. 0001) (Fig. 6a-c ), indicated that a heavier disease burden was more likely to be found in higher SDI regions.
section_name: The correlation between SDI and CLL's disease burden

section_content: Based on GBD study 2019, four potential CLL-related mortality and DALY attributable risk factors including high body mass index, occupational exposure to benzene, occupational exposure to formaldehyde, and smoking were identified. Among these risk factors, smoking was the strongest risk factor to CLL-mediated death and DALY from 1990 to 2019 at a global scale (Fig. 7a-d ). Of note, compared with high-SDI areas, the proportion of CLL's disease burden attributable to high body mass index in low-SDI areas has a significant upward trend. In addition, although the percent of CLL deaths and DALYs attributed to occupational carcinogen-exposure only accounted for a very small proportion, a significantly higher risk of carcinogen exposure was found in low-SDI regions compared to high-SDI regions. 

CLL, chronic lymphocytic leukemia; DALY, disability adjusted life-year; EAPC, estimated annual percentage changes; SDI, social-demographic index
section_name: CLL burden attributable risk factors

section_content: Our research focused on the current status and trends of the global burden of CLL based on the latest GBD study 2019 database. In this study, we collected the global incidence, mortality, DALY data attributable to CLL and evaluated the epidemiological trends from 1990 to 2019. The results revealed that global burden disease of CLL presented a constant growing trend during the past 30 years. Previous epidemiologic evidence has demonstrated that CLL predominately occurred in the elderly, with median diagnostic age above 70 years-old [2]. In our study, an age distribution of over 50 was found occupied the vast majority of CLL population, which is consistent with the previous reports. Aging affects the hematopoietic system. The aging of hematopoietic stem cells is accompanied by a series of biological changes including the accumulation of genomic damage, epigenetic changes, telomere shortening, and oxidative stress, which are closely associated to the occurrence of abnormal clonal hematopoiesis [20]. According to reports, in CLL, the ability to produce cloned B cells may have been acquired at the stage of hematopoietic stem cells [21]. Therefore, hematopoietic stem cells that have accumulated multiple key mutations accompanying aging are considered to play a central role in the occurrence of CLL. Recurrent Fig. 4 The incidence cases (A) and corresponding asir (B) of CLL across different age groups from 1990 to 2017 in the globe and various regions stratified by SDI. CLL, chronic lymphocytic leukemia; ASIR, age-standardized incidence rate; SDI, socio-demographic index mutated genes such as NOTCH1, MYD88, TP53, etc. have been identified in the occurrence of CLL [22]. Furthermore, based on the SDI-stratified regional analysis, we revealed that the majority of incidence cases were between 50 and 69 in regions with low SDI, while more than half of the incidence cases were over 70 in regions with high SDI. A generally expanding aging of the population in regions with high SDI may account for the differences between two regions. Besides, advanced ages imply a worse prognosis, accompanied with an increased disease burden due to bad health status and poor tolerance to chemotherapy toxicity [23]. Therefore, preferential attention should be paid to the rapid increase of CLL considering the current context of global population aging. 

In the global context, the incidence and mortality of CLL in both genders displayed an increasing trend during the past 30 years, with males presenting a relatively larger proportion. Moreover, we found significantly divergent gender distribution of disease burden among different SDI regions. In regions with low SDI, we found that females account for the majority of incidence and mortality, which is in contrast to the global trends. Although several epidemiological studies in underdeveloped areas have reported similar gender ratio as our investigations, these studies are primarily limited to local place and contained insufficient sample size [24]. Even considering the poor health of local females due to unequal social status and the stronger tendency for females to seek medical assistance compared to males [25], the available evidence is not sufficient to fully explain the contrast incidence. Given the rapid increase in the incidence of CLL in these areas, it is necessary to conduct in-depth multi-center, large-sample CLL epidemiological investigations. Additionally, policies and strategies to improve the health status of women in these regions should be a priority. As for geographical variation factors, the results revealed that the incidence and mortality of CLL were higher in North America, Europe, and Australia, and lower in Asia, which is approximately consistent with previous studies recognized Caucasian race as risk factors [4]. Interestingly, when focused on the Middle East, we found that although CLL generally presented a low incidence in the Arab world, it is clearly manifested as a high burden disease pattern in Israel. Such situation could be explained by the main demographic composition of Israel which is dominated by Jews, and we speculated that other social environmental factors may also contribute to this difference. 

To our knowledge, variations were observed in CLL burden across different SDI quintiles, with apparently heavier burden in regions with higher SDI. Notably, middle-SDI regions presented rapidly increasing incidence and mortality trends compared with high-or low-SDI regions. The possible explanation might be the huge imbalances of local healthcare environment and settings existing worldwide. Wide coverage of cancer screening and demographic characteristics of aging population are prevalent in high-SDI regions, contributing to a relatively stable disease burden growth trend although with high incident rates. In low-SDI regions, a long-standing lack of screening conditions, possible missed diagnosis, and incomplete case reports caused certain detection biases, leading to constant underestimation of the incidence and mortality. Besides, a catch-up development in annual increases of SDI was recorded in some developing countries [19]. In recent years, given that the improvement of basic medical conditions and the emphasis on early prevention in some middle-SDI regions, lost morbidity and mortality due to missed diagnosis and case underreporting are reducing gradually, thus exhibiting a trend of rapidly increasing disease burden [19]. Overall, these findings prompted us to rationally mobilize existing resources to accurately evaluate and reduce the CLL burden in less-developed areas, which required enhanced disease detection and early treatment management. 

We also investigated the risk factors that affected CLL-related mortality and DALY. Smoking is the major contributor of the 4 risk factors across the world during the last 30 years. More than 60 compounds were identified as known carcinogens in tobacco smoke [26]. Although there is lack of definite association between tobacco use and CLL incidence up to date [27], various cohort studies have reported the association between smoking and the occurrence of myeloproliferative tumors [28, 29]. A striking variation in smoking rates was found among countries in a study investigating global trends for tobacco use from 1990 to 2010, with tobacco use in high-income countries effectively controlled, which may be partly attributed to the increase in awareness of smoking cessation and implementation of national tobacco control policies. Unfortunately, smoking prevalence is on the rise in low-income countries, suggesting that some interventions such as tobacco excise taxes and smoking cessation propaganda should be considered [30, 31]. 

High body mass index also behaves as an important risk factor for CLL, with a rapidly increasing trend especially in low-SDI regions. The difference in lifestyle between developed and underdeveloped regions may explain this difference in trends. In the past few decades, dependence on processed foods brought a high-sugar and fat diet across the world. Besides, mechanized production has replaced the original manual labor in most areas. These changes have accelerated a global obesity pandemic, which is even severer among low-and middle-income populations [32]. A meta-analysis on the correlation with obesity during adulthood and risk of lympho-hematopoietic cancers revealed that general adiposity in adulthood and early adulthood may increase the risk of CLL [33]. 

Later, several studies analyzing the impact of obesity on CLL patients indicated that a poorer baseline response to induction treatment, a lower complete remission rate, and a reduced progression-free survival time were observed in obese patients [34, 35]. According to these evidences, additional attention should be paid to advocate healthy diet and reasonable exercise in the public. 

Besides, the exposure risk of benzene and formaldehyde is significantly higher in low-SDI regions than in high-SDI regions. Inhalation is the predominant way for occupational exposure to benzene and formaldehyde. A study evaluating the genetic effects of long-term occupational exposure to formaldehyde showed that long-term exposure to formaldehyde caused higher frequencies of micronuclei in nasal cells and higher frequency of sister chromatid exchanges of peripheral lymphocytes [36]. A meta-analysis accessing the correlation between benzene exposure and leukemia showed that exposure to benzene at work increased the risk of AML and CLL in a dose-response pattern [37]. People at high risk of exposure to benzene and formaldehyde include workers in paint factory, shoe factory, furniture factory, and decorator. With globalization, a large number of manufacturing factories have moved to underdeveloped regions [38]. Meanwhile, strict control of carcinogenic occupations prevalent in developed countries has not been fully implemented in underdeveloped regions. Therefore, special attention should be paid to the risk of occupational exposure to carcinogens in these underdeveloped areas. 

Since there are limited available epidemiological studies on CLL in a global perspective, our research based on GBD study 2019 provides the latest global epidemiological distribution and trends on CLL for future research. However, several limitations are unavoidable in the study. As mentioned above, although GBD covers the data of disease burden in most countries and regions in the world, morbidity and mortality in underdeveloped regions may be underestimated due to missed diagnosis and lack of reliable disease information systems. This detection error is difficult to be completely corrected by the subsequent re-allocation algorithm of GBD. Secondly, GBD includes classifications of disease populations in geographical areas, but lacks ethnicity data to facilitate the analysis of genetic susceptibility. Thirdly, limited potential risk factors of CLL are included in the GBD, hampered further research on the distribution and trends of CLL risk factors.
section_name: Discussion

section_content: In summary, the global burden of CLL has maintained a gradually increasing trend from 1990 to 2019. The disease tended to occur in males, the elderly populations, and people living in high-SDI regions. What cannot be ignored is the rapid growth of the disease burden in middle-SDI regions, which potentially indicated an underestimated incidence and mortality in underdeveloped countries. In addition, of attributable risk factors, smoking presented as the most contributed across the globe, with potential risk of carcinogen exposure containing a prominent issue in low-SDI regions which needs further investigation. Based on the evaluation of the increasing CLL global burden trends and the highly heterogeneous distribution pattern, policy-makers could assess the effectiveness of previous prevention strategies and rationally adjust the follow-up health policies to alleviate the growing burden.
section_name: Conclusion

section_name: Methods

section_content: GBD study 2019 data resources were available online from the Global Health Data Exchange (GHDx) query tool (http:// ghdx. healt hdata. org/ gbd-resul ts-tool), containing a global collection of epidemiological data evaluating burden of disease worldwide with 369 diseases across 204 countries and territories [39]. The data for this study were collected from GBD study 2019, including the following epidemiological data:
section_name: Data source

section_content: CLL-related incidences, deaths, DALYs, corresponding ASRs and the SDI are directly provided by the GBD study 2019. ASRs (including ASIRs, ASDRs, age-standardized DALY rates) could eliminate the interference of age structure and population size. DALY, refers to the total healthy life years lost due to disease from morbidity to death, is obtained by adding up the years lived with disability (YLDs) and the years of life lost (YLLs). SDI is an intuitive indicator which reflects the degree of social development [40]. The value of SDI which ranged from 0 to 1 were calculated through the comprehensive evaluation of fertility, education, and income to reflect the degree of social development. Countries and territories were then stratified into five levels (high SDI, high-middle SDI, middle SDI, low-middle SDI, and low SDI) according to SDI values obtained from GHDx. The attributable risk factors for CLL from the GBD study 2019 include: high body mass index, smoking, exposure to benzene, and exposure to formaldehyde. This study collected the above risk factors and presented the results through data visualization. 

To assess trends in the CLL's disease burden, EAPC values (including EAPC based on ASIR, ASDR, and age-standardized DALY rate per 100,000 persons) were employed in this study. Taking the calendar year as the independent variable x while taking the natural logarithmic transformation of the ASR as the dependent variable y. A regression equation was fitted to the natural logarithm of the rates, i. e., y = α + βx, where x = year and y = ln (ASR). Further EAPC values were calculated according to the formula EAPC = 100* (e β -1). For the 95% confidence intervals obtained from the regression, if there is an overlap within 0 in the 95% confidence intervals, then the corresponding ASR is regarded stable. If the 95% confidence intervals fall entirely below 0, then the corresponding ASR is regarded declined, otherwise considered to be elevated [40]. 

To assess whether the ASIR in low-incidence areas has been underestimated, we calculated the Pearson correlation coefficient between ASIR-based EAPC and ASIR in 204 countries and territories in 1990. Then, to analyze the correlation between the ASR trends and SDI, we calculated the Pearson correlation coefficient between the ASRs and SDI values of 21 geographical regions from 1990 to 2019. The Pearson correlation coefficient, ranged from -1 to 1, was used as a measure of linear correlation between two variables. The closer the absolute value of the correlation coefficient is to 1, the stronger the relevance [39]. 

All data analyses in this study were performed with R software (R software, version 4. 0. 0). The R packages employed include: dplyr, stringr, Rcan, maps, and ggplot2. Dplyr and stringr were used for data collation. Rcan was used to calculate EAPC. Maps and ggplot2 were used for visualization and further correlation coefficient calculation. The map was used to present the current status of CLL burden worldwide. The histogram was used to show the changing trend of CLL-mediated disease burden. The scatter plot and regression curve were used for correlation analysis between ASRs and SDI. All tests were two-tailed, and a P value of less than 0. 05 was considered as statistically significant.
section_name: Statistical analysis and data visualization

section_name: Acknowledgements

section_content: Not applicable.
section_name: Acknowledgements

section_name: Funding

section_content: This work was supported in part by National Natural Science Foundation of China (No. 81800146 ); Key research and development program of Zhejiang Province, China (No. 2021C03123 ); Key international cooperation projects of the National Natural Science Foundation of China (No. 81820108004 ); Youth Natural Science Foundation of Zhejiang Province, China ( LQ18H080001 ).
section_name: Funding

section_name: Availability of data and materials

section_content: The datasets generated and/or analyzed during the current study are available from the Global Health Data Exchange query tool (http:// ghdx. healt hdata. org/ gbd-resul ts-tool).
section_name: Availability of data and materials

section_name: Abbreviations

section_content: CLL: Chronic lymphocytic leukemia; GBD: Global Burden of Disease; SDI: Social-demographic index; DALY: Disability adjusted life-year; BTK: Bruton tyrosine kinase; EAPC: Estimated annual percentage change; ASIR: Age-standardized incidence rate; ASDR: Age-standardized death rate; ASR: Age-standardized rate; GHDx: Global Health Data Exchange; YLD: Year lived with disability; YLL: Year of life lost.
section_name: Abbreviations

section_content: The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12938-021-00973-6.
section_name: Supplementary Information

section_content: Authors' contributions YY, XL, FL contributed to the conception and drafting of the manuscript and figures. JJ and HW reviewed the manuscript and gave final approval of the version to be published. All authors read and approved the final manuscript.
section_name: Additional file 1. Supplementary materials.

section_content: Ethics approval and consent to participate Not applicable.
section_name: Declarations

section_content: Not applicable.
section_name: Consent for publication

section_content: The authors declare that they have no competing interests. 

• fast, convenient online submission • thorough peer review by experienced researchers in your field
section_name: Competing interests

section_content: • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year
section_name: • rapid publication on acceptance

section_content: At BMC, research is always in progress.
section_name: •

section_content: Ready to submit your research Ready to submit your research ? Choose BMC and benefit from:

? Choose BMC and benefit from:
section_name: Learn more biomedcentral.com/submissions

section_content: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
section_name: Publisher's Note