Global, Regional, and National Trends of Chagas Disease from 1990 to 2019: Comprehensive Analysis of the Global Burden of Disease Study

Background: Chagas disease (CD) is a neglected tropical disease, endemic in Latin America, but due to migration and environmental changes it has become a global public health issue. Objectives: To assess the global prevalence and disability-adjusted life years due to CD using findings from the Global Burden of Disease Study 2019. Methods: The Global Burden of Disease data was obtained from the Global Burden of Disease Collaborative Network; results were provided by the Institute for Health Metrics and Evaluation. The prevalence and disability-adjusted life-years (DALYs) were described at a global, regional, and national level, including data from 1990 to 2019. Results: Globally, CD prevalence decreased by 11.3% during the study period, from 7,292,889 cases estimated in 1990 to 6,469,283 in 2019. Moreover, the global DALY rate of CD decreased by 23.7% during the evaluated period, from 360,872 in 1990 to 275,377 in 2019. In addition, significant differences in the burden by sex, being men the most affected, age, with the elderly having the highest burden of the disease, and sociodemographic index (SDI), with countries with the lowest SDI values having the highest prevalence of the disease, were observed. Finally, the prevalence trends have followed different patterns according to the region, with a sustained decrease in Latin America, compared to an increasing trend in North America and Europe until 2010. Conclusion: The global burden of CD has changed in recent decades, with a sustained decline in the number of cases. Although the majority of cases remain concentrated in Latin America, the increase observed in countries in North America and Europe highlights the importance of screening at-risk populations and raising awareness of this neglected tropical disease.


INTRODUCTION
Chagas Disease (CD), also known as American trypanosomiasis, is a neglected tropical disease caused by the protozoa Trypanosoma cruzi, which is endemic in 21 countries of continental Latin America. CD's geographic distribution was mostly driven by the distribution of the vector species in the Americas regions below 1500m elevation. However, due to migration, it has become a global health threat [1]. According to the World Health Organization, six to eight million people are estimated to be infected with T. cruzi, and another 70 million are estimated to be at risk of infection [2]. This infectious disease has two phases: acute and chronic; the latter is most commonly associated with myocardial involvement, known as Chronic Chagas Cardiomyopathy (CCM) [3]. The pathophysiology of the progressive myocardial damage remains poorly understood; however, recent evidence suggests an essential role of T. cruzi persistence (both in blood and tissues) and a pro-inflammatory immunological profile in the host, which is reflected in sustained oxidative stress, and natural killer/CD8+ T-cell cytotoxicity [4,5].
The conjunction of heart failure and systemic embolisms in Chronic Chagas Cardiomyopathy (CCM), the most severe form of organ involvement in CD, has been associated with a significant burden in the affected patients [6]. CCM has been associated with higher morbidity and mortality rates compared to other cardiomyopathies, mainly due to severe myocardial injury, which not only leads to heart failure, but also favors conduction disorders and new-onset arrhythmias such as atrial fibrillation, which, together with the presence of myocardial aneurysms, leads to a significantly increased risk of embolic events [7]. On the other hand, it is estimated that CD leads to a global expenditure of USD 627.5 million per year, mainly due to health care costs. Precisely, it is estimated that the cost per patient rises from $200 in the early stages of the disease to $6,000 in the chronic forms [8]. As a result, CD has been ranked as the parasitic disease with the highest attributable burden worldwide [9].
Although some mathematical models and systematic reviews have aimed to estimate CD's global burden, a detailed quantitative analysis of this disease burden is unfortunately not available [10]. We aimed to characterize in detail the global state of CD epidemiology by using the Global Burden of Disease (GBD) data from 1990 to 2019, highlighting the differences by sex and the most affected ages across regions and nations.

DATA SOURCES
The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) is an international initiative aimed to provide a systematic scientific evaluation of published data regarding disease and injury incidence, prevalence, and mortality for a wide list of diseases and injuries. We used the GBD 2019 results to assess the global burden of Chagas Disease (CD) by extracting data on CD prevalence, death due to CD, disability-adjusted life-years (DALYs), and their uncertainty intervals from GBD data sources. In summary, the GBD study performed an epidemiologic assessment of 333 diseases by age and sex globally for the years 1990 to 2019. The method by which the GBD estimates Chagas Disease burden in endemic countries is based on seroprevalence data with an integrative Bayesian meta-regression technique, which assesses a generalized negative binomial model for all epidemiological data through DisMod-MR 2.1. On the other hand, in non-endemic countries, the estimates are calculated by estimating the importation of prevalent cases via immigration. Data from the Pew Research Trust's Center (https://www. pewresearch.org/hispanic/2014/04/29/2012-statistical-information-on-immigrants-in-unitedstates/) and the International Migrant Stock of the United Nations are used to estimate the number of foreign-born individuals by age, sex, and country of origin for each state. Finally, the state-specific estimates of the number of people born in CD-endemic countries and the estimated prevalence from each country are used to generate state-specific CD prevalence estimates.
We obtained the GBD data from the Global Burden of Disease Collaborative Network; results were provided by the Institute for Health Metrics and Evaluation (Seattle, WA). In this context, the GBD cause list is organized hierarchically into four levels, Chagas disease is in level 2 group of "Neglected tropical diseases and malaria" and level 1 group of "Communicable, maternal, neonatal, and nutritional diseases." Furthermore, Chagas disease's case definition was based

MEASURES OF BURDEN
The measures of burden assessed in the present study included prevalence, death due to CD, and DALYs due to CD. At first, DALYs were extracted directly from the GBD. This measure is calculated by adding years lived with any short-or long-term disability and years of life lost (YLL) due to premature death. The years lived with disability (YLD) is a burden measure that represents the number of years lived with CD considering the related disability weights to reflect CD's underlying severity. Both YLL and YLD were estimated specifically by age, sex, and country. All measures were reported as number counts and age-standardized rates per 100,000 population. This age standardization was based on the World Health Organization world population standard age structure.

STATISTICAL ANALYSIS
Global, regional and national data from the countries with a prevalence rate per 100.000 population higher than 0.1 was presented. Results of prevalence, deaths, DALYs, YLLs, and YLDs were presented in absolute numbers and age-standardized rates per 100,000 population, along with their respective 95% uncertainty intervals (UIs). We also reported the percentage of change in the values between 1990 and 2019 to show the variations in this period. Finally, we performed a literature review in Pubmed and LILACS to assess new primary data concerning the prevalence of CD in endemic countries to generate an updated estimate and compare it with the GBD data and other estimates published in the literature. Furthermore, we analyzed the CD prevalence rate per 100,000 population between countries according to their sociodemographic index classification (SDI), a composite indicator of development status developed by GBD researchers (https://ghdx.healthdata.org/record/ihme-data/gbd-2019-sociodemographic-index-sdi-1950-2019). The Kruskal-Wallis one-way analysis of variance test was used to compare the prevalence rate among the different SDI categories. A p-value < 0.05 was considered significant. All tests were two tailed. All statistical analyses were performed using R Statistical Software (version 4.2.1; R Foundation for Statistical Computing, Vienna, Austria). difference represented a 11.3% decrease in prevalence over the last three decades. Similarly, the age-standardized prevalence rate per 100,000 population in 2019 decreased to 83.61 (95% UI = 73.13-95.31) representing a 42% reduction of this value with respect to 1990 (Supplementary Figure 1). Bolivia remained the country with the highest prevalence rate per 100,000 population (4993.53; 95% UI = 4540.19 -5483.11), followed by Venezuela (1654,73; 95% UI = 1439.50-1897.81) and Argentina (1524.07; 95% UI = 1343. 16 -1736.00) (Supplementary Figure 2). On the other hand, Brazil remained the country with the higher number of estimated CD cases, with around 2,164,570 cases (33.5% of the total). Furthermore, Argentina was the endemic country with the largest reduction in the prevalence number (-54.68%), followed Chile (-50.95%) and Uruguay (-49.95%), while Mexico (39.58%), Brazil (12.80%) and Honduras (14.12%) were the endemic countries with the largest increase in estimated number of cases in the period between 1990 and 2019. Conversely, none of the endemic countries had an increase in the prevalence rate per 100,000 population. On the other hand, the estimates suggest that Spain was the country with the largest absolute increase in the number of estimated CD cases (25,081 cases), followed by the United States of America (10,048 cases) and Italy (5542 cases). Finally, the U.S.A remained the non-endemic country with the largest estimated number of cases by 2019 (63,553.17; 95% UI = 55,365.83-72,997.92), being California the state with the highest prevalence number (13,600; 95% UI = 11,761.33-15730.08). Table 1 and Figure 1 present a summary of the prevalence data by country and U.S.A states.

BURDEN BY REGIONS
As expected, the Latin American region has suffered the highest burden of the disease, with a total of 6,354,220 individuals with CD diagnosis by 2019, accounting for a prevalence rate per 100,000 population of 933.76 (95% UI 817.66-1064.29). Nevertheless, this prevalence rate has dropped in the last 30 years, as shown in Supplementary Figure 6. A different scenario occurred in North America's region, which showed a progressive increase in the prevalence rate until 2010, followed by a sustained drop (Supplementary Figure 7). Finally, CD prevalence in the European region showed a similar trend, with a substantial increase from 1990 until 2010, followed again by a drop in the cases rate after this point (Supplementary Figure 8). Figure 2 shows a summary of the prevalence trends in the most affected countries/U.S. states by region. Table 1 and Supplementary Table 1 summarize the prevalence rates and numbers and their change between 1990 and 2019 by country and U.S. state.
On the other hand, DALY rates per 100,000 population showed a similar trend in the Latin American region during this period, with a sustained drop in its value from 116.72 (95% UI 48.85-145.08) in 1990 to 40.54 (95% UI 27.04-67.75) in 2019. Regarding the North American region, a steady increase in the DALY rate was observed until 2015 (0.39 DALYs per 100,000 population; 95% UI 0.28-0.5), followed by a steep decline thereafter. Finally, a relatively stable trend was observed in the European region, with an initial drop in the DALY rate after 1994 followed by a period without significant changes in this value.  Finally, we evaluated a potential association between the sociodemographic index (SDI) classification of the assessed countries and the CD prevalence rate. Countries classified as low-middle SDI had the highest CD prevalence rate (median value: 691 cases per 100,000 population), followed by those in the middle SDI group (median value: 345 cases per 100,000 population). On the other hand, the countries in the high-middle and high SDI groups had significantly lower prevalence rates (7.88 cases and 2.15 cases per 100,000 population, respectively). Figure 3 summarizes the CD prevalence rate of the evaluated countries according to the respective SDI category, highlighting a significantly larger rate in the low-middle SDI group compared to the high-middle (p-value: 0.02) and high SDI (p-value < 0.001) groups, along with a significantly higher prevalence in the middle SDI group compared to the high SDI one (p-value: 0.012). Finally, Figure 3 also highlights potential outliers among the classification groups, including Argentina, Chile, and Uruguay in the high-middle SDI group and Bolivia in the Low-middle SDI group.

BURDEN BY AGE AND SEX
Supplementary Figure 9 summarizes the estimated prevalence rate per 100,000 population by age and sex in 1990 and 2019. At first, a male-to-female ratio of around 0.9:1 was observed by 1990, having women a slightly higher prevalence mainly across ages below 60 years with a peak in the age group of 45 to 49 years (220.48 per 100,000 population [95% UI: 192.32-251.63]) and a progressive reduction in the prevalence rate after this point. On the other hand, males showed a peak in the age group of 95 plus years (302. 31 Figure 4 shows a similar trend of DALYs by sex across age groups, highlighting a higher burden in men in almost all categories and a male-tofemale ratio that increases progressively as age increases. It is relevant to highlight the differences between regions regarding the gender gaps observed. For example, we observed a large DALYs rate male-to-female ratio in Latin America, which remained nearly unchanged during the study period (1.68:1 in 1990 and 1.69:1 in 2019), evidencing that men have substantially higher DALY rates than women in this region (Supplementary Figure 10). On the other hand, gender differences in the North America's region were less stable, with a DALY rate male-to-female ratio of 1.22:1 in 1990 that was followed by a peak by 2015 (2.43:1) and a subsequent reduction until 2019 (1.5:1) (Supplementary Figure 11). In parallel with the observed reduction in the rate of DALYs in the European region, a reduction in the gender gap was observed in this region, from 1.64:1 in 1990 to 1.31:1 in 2019 (Supplementary Figure 12).
Finally, a sustained decrease in the rate of deaths attributed to CD over time was observed in both sexes. Nevertheless, males had a higher mortality rate per 100,000 population across the almost three decades of analyzed estimates, with a male-to-female ratio of 1.6:1 in 1990 and 1.5:1 in 2019 (Supplementary Figure 13). Moreover, the countries with the largest maleto-female ratio for mortality in 2019 were El Salvador (2.4:1), Ecuador (2.4:1), and Argentina (2.2:1).

DISCUSSION
According to the GBD data, there has been a sustained decrease in Chagas Disease burden, measured in prevalence, deaths, and DALYs due to CD worldwide from 1990 to 2019. Nevertheless, this trend has not been the same for every affected region, with a more abrupt reduction of the burden in Latin America, while data from North America and Europe reveal a more stable trend, with increasing prevalences until 2010. Moreover, we highlighted significant differences in the burden by sex, being men the most affected, age, with the elderly having the highest burden of the disease, and sociodemographic index (SDI), with countries with the lowest SDI values having the highest prevalence of the disease. These findings have relevant public health implications; however, they need to be assessed in light of the current published evidence estimating CD's burden in each region. The main finding of the present study is the substantial reduction in the prevalence of the disease and the burden attributed to it worldwide. Although interregional differences were observed, the reduction in the number of cases and deaths secondary to CD may reflect the impact of decades of public policies aimed at its mitigation. Initially, the elimination of vectorborne transmission became the goal of governments of endemic countries, which led vector eradication initiatives since the 1970s. These programs included the restoration of homes to prevent triatomine infestation, as well as the use of synthetic insecticides for vector control, among others [11]. The impact of these programs was initially measured in the number of municipalities with triatomine infestation, and in the case of Brazil, in which a significant decrease in this number was observed, especially in those infested by Triatoma infestans. This decrease was also associated with a reduction in the number of cases of acute Chagas disease and in the prevalence of infection in children under five years of age [12]. Furthermore, climate change has been postulated as one of the potential determinants of CD transmission in the region and the world. Despite this, studies that have analyzed the impact of temperature change on parameters such as the life cycle, fertility and fecundity of vectors have yielded divergent results [13][14][15][16][17]. Therefore, there is a clear need for studies that comprehensively evaluate the relevant ecological and climatic factors in the parasite-vector interaction.
The GBD Study estimations summarized in the present study differ vastly from what has been calculated by other studies published in the literature. At first, the prevalence reported in the U.S. (in 1990 and 2019) was substantially lower than that estimated by the study of Bern et al., which estimated a prevalence number of 300,167 cases of individuals with T. cruzi infection living in the U.S. by 2009, with 30,000-45,000 cases of cardiomyopathy [6]. Similar results were reported in the study of Manne-Goehler et al., with a national estimate of 238,091 cases as of 2012 [9]. The differences in these estimations could be related to the varied underlying estimates of immigration and seroprevalence by region used in each study. For example, the study of Bern et al. estimated the number of Latin American migrants from CD endemic countries using information from the Pew Hispanic Center and the reports of unauthorized immigrant populations published by the Department of Homeland Security [6]. On the other hand, Manne-Goehler et al. used information from the American Community Survey to estimate the number of residents born in CD endemic countries [9]. Both studies relied on data from the World Health Organization to estimate the prevalence of T. cruzi infection in 2006 and 2015 [6,9].
The results of this study also differ from previous estimations of CD prevalence in the European Region. For example, the study of Gascon et al. estimated that a total of 47,743 foreign-born individuals infected by T. cruzi could be living in Spain by 2009 [18]. This study also used data from the PAHO/WHO to estimate CD's prevalence and burden in each endemic country. On the other hand, the systematic review of Navarro et al. estimated that there were around 34,202 Bolivian immigrants with T. cruzi infection living in Spain by 2010, highlighting a much higher total prevalence when considering immigrants from other nationalities [19]. Nevertheless, this review included studies conducted in hospital-based settings; therefore, the results can potentially overestimate the actual population prevalence due to selection bias. A way to prevent this overestimation is performing separate estimates from antenatal care screenings, primary health care/community-based studies, and even blood bank data. However, it is important to highlight that studies based on blood bank data tend to underestimate CD prevalence [20].
Considering that the PAHO/WHO estimates may not be the only sources of information regarding the prevalence of CD in endemic regions, we reviewed the studies that have provided national estimates of this prevalence in endemic countries, evaluating three studies assessing the nationwide prevalence of CD in endemic countries. Our results suggest that even the previously mentioned studies, which relied on PAHO/WHO prevalence estimates, could underestimate the actual prevalence of CD in non-endemic countries, as the estimates from the studies included have suggested higher prevalence values compared to those reported by the PAHO/WHO in 2006 and 2015 [21,22]. At first, the meta-analysis of Arnal  was also reported by Olivera MJ et al. in Colombia [25]. Considering only the estimated migrants from these three countries (Colombia, Brazil, and Mexico) in the U.S. (according to the 2019 American Community Survey) and the updated prevalence data, an estimate of 406,748 individuals from these endemic countries that are currently living in the U.S. could be infected with the parasite. Adding this value to the estimated number of CD cases in immigrants from other endemic countries using the PAHO/WHO data on CD prevalence would yield an estimate of 472,397 CD cases in the U.S. alone. This updated data suggest that the burden of the disease could be even higher than what has been estimated before, highlighting the importance of having new up-to-date national estimates on CD prevalence in endemic regions.
Furthermore, important differences in the burden by sex were observed, highlighting a higher burden of CD in males. At first, prevalence differences can be explained by biological, sociocultural, and behavioral factors, such as the exposure to the vector as a result of occupational activities, as there is evidence that suggests a potentially relevant role of vector transmission in working areas [26,27]. Moreover, male sex has been identified as an independent risk factor for progression from the indeterminate form of the disease to chronic Chagas cardiomyopathy, which may promote the diagnosis due to the onset of HF symptoms [28]. Regarding the burden of the disease, several studies have observed a higher risk of morbimortality in males due to cardiac complications of CD [29]. Recently, significant differences in myocardial damage by sex have been observed, suggesting potentially different pathophysiological mechanisms involved in the development and progression of cardiomyopathy. In the study of Assunção Jr et al., in which gender differences in myocardial damage assessed by cardiovascular magnetic resonance were evaluated, a significant association between gender and myocardial dysfunction was observed, with men showing greater myocardial fibrosis and lower LV ejection fraction values. Moreover, the transmural pattern of involvement of the myocardial tissue along with the myocardial gray zone were more prevalent among males [30]. In addition, access to health services may represent a relevant factor differentiating outcomes by sex, highlighting the hesitancy of males in searching for health services until the appearance of severe incapacitating symptoms or complications, such as chronic Chagas cardiomyopathy [31][32][33]. This leads to a late diagnosis of the disease, potentially explaining the higher prevalence in older ages in males. On the other hand, the higher DALY rate in the elderly may be related to the chronic nature of the disease and the increasing prevalence of chronic comorbidities in older ages [34].
Finally, an aspect that has not been thoroughly assessed is the DALYs associated with the disease. The studies of Martins-Melo et al. used the GBD study data to estimate the burden of CD and other neglected tropical diseases in Brazil, highlighting an important reduction in the burden from 1990 to 2016, as reported in our study [35,36]. On the other hand, the burden associated with CD in terms of DALYs has not been studied in non-endemic regions, with the only information available coming from the computational simulation model of Lee

LIMITATIONS
The present study suffers from several limitations. At first, it was based on the GBD study data to generate the estimates provided. Although the methodology of the GBD is considered reliable, we must acknowledge that it can be limited by the quality of the available data, which can vary across countries and time, relying on mathematical models to assess the burden in countries with limited information. In this context, the GBD used the data of CD prevalence by country from the Pan American Health Organization/World Health Organization to estimate the number of T. cruzi infection cases; however, these reports may provide inaccurate estimations, as some were based on extrapolations from mathematical models using vector distribution or from blood donor screening programs. Specifically, the high prevalence of asymptomatic individuals with CD, even in those with cardiomyopathy, may directly impact the prevalence estimates.
Moreover, the information reported was based on multiple assumptions, for example, that the state-wise distribution of foreign-born persons is similar for all countries within a region and that the prevalence of CD in the immigrant population is similar to that of the country of origin. Furthermore, the global trends observed in this study may not correlate with micro-level trends for some regions, such as the United States of America, with may show an important subnational variability in the burden of CD. Therefore, the CD trends observed in this study must be interpreted with caution, especially for those estimates with larger uncertainty intervals.

CONCLUSIONS
The global burden of CD has been dramatically changing in the last decades, with a sustained decrease in the number of cases and the disease's burden. Despite the vast majority of the cases still concentrate in Latin America, a significant prevalence is reported in developed countries of North America and the European Region. Although it provides valuable detailed information for health policymakers, data from the GBD study can potentially underestimate the prevalence of the disease in non-endemic countries. There is a relevant need to carry out representative seroprevalence studies and screen migrant populations to provide more precise estimates of the disease prevalence in non-endemic countries. Moreover, raising the awareness of this neglected tropical disease in these regions is critical for mitigating its impact and allowing earlier diagnoses, therefore, potentially improving the clinical outcomes of CDinfected patients.

DATA ACCESSIBILITY STATEMENT
To download the data used in these analyses, please visit the Global Health Data Exchange GBD 2019 website (http://ghdx.healthdata.org/gbd-results-tool).

ADDITIONAL FILE
The additional file for this article can be found as follows: