Acessibilidade / Reportar erro

Impact of redistributing deaths by ill-defined causes in oral and oropharyngeal cancer mortality in Brazil

Abstract

Less-than-optimal reliability of mortality information systems regarding the underlying cause of death can mask the reality of oral (OC) and oropharyngeal cancer (OPC) mortality. This study aimed to assess the impact on the magnitude and temporal trends of OC and OPC mortality in Brazil of two statistical approaches to redistribute deaths with ill-defined underlying causes. We analyzed deaths with ill-defined causes in Brazil by macro-region, between 1996-2018. The Mortality Information System provided official information on deaths. Two correction methods were applied: the EF method, which proportionally reallocates deaths classified as R00-R99 in the ICD-10 to the remaining specific causes of death according to the proportion of deaths with certified causes; and the GBD method, which considers the concept of garbage codes, redistributing deaths from several ICD-10 chapters according to previously established coefficients. For the trend analysis of mortality (certified and redistributed), the Prais-Winsten method was carried out. The OC and OPC death rates had an evident increase after the redistribution by the two techniques in all regions of the country; the increase was higher using the GBD method. In the Northeast and North regions, this method more than doubled the certified death rates. The redistribution methods also changed time series trends. In epidemiological studies of mortality from OC and OPC, it is necessary to redistribute deaths from ill-defined causes when analyzing data from less-than-optimal information systems. The choice of the correction method is critical; epidemiological studies must manage it as a methodological decision that has significant impacts on results.

Mouth Neoplasms; Oropharyngeal Neoplasms; Mortality; Health Information Systems

Introduction

Oral (OC) and oropharyngeal (OPC) cancers are a relevant cause of mortality worldwide. In 2017, they caused about 4% of all cancer deaths in the world.11. Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, Abdel-Rahman O, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-Adjusted life-years for 29 cancer groups, 1990 to 2017: A systematic analysis for the global burden of disease study. JAMA Oncol. 2019 Dec;5(12):1749-68. https://doi.org/10.1001/jamaoncol.2019.2996
https://doi.org/10.1001/jamaoncol.2019.2...
This percentage accounted for more than 220,000 deaths in 2018, according to the International Agency for Research on Cancer (IARC). Mortality from these malignancies varies widely by region, and OC is the second most common cancer leading to death in regions like South Central Asia.22. Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Piñeros M, et al. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer. 2019 Apr;144(8):1941-53. https://doi.org/10.1002/ijc.31937
https://doi.org/10.1002/ijc.31937...
The mortality rate form OC and OPC estimated by IARC for Brazil in 2020 was 4.9 deaths per 100,000 inhabitants in men in 2020, which was much higher than in Latin America and the Caribbean: 2.6/100,000.33. International Agency for Research on Cancer. Cancer today: data visualization tool for exploring the global cancer burden in 2020 [Internet]. 2021 [cited 2021 Feb 25]. Available from: https://gco.iarc.fr/today/home
https://gco.iarc.fr/today/home...

Tracking oral cancer is one of the priorities of the World Health Organization’s non-communicable diseases and universal health coverage agenda.44. World Health Organization. Oral health: achieving better oral health as part of the universal health coverage and noncommunicable disease agendas towards 2030. Geneva: World Health Organization; 2020 [cited 2021 Feb 25]. Available from: https://apps.who.int/gb/ebwha/pdf_files/EB148/B148_8-en.pdf
https://apps.who.int/gb/ebwha/pdf_files/...
Monitoring the magnitude and time trends of OC and OPC mortality is a tool for understanding their dynamics across time and defining public health policy and programming. This measure becomes more critical in the absence of reliable records of cancer incidence, which is the reality in several countries, such as Brazil.55. Maruthappu M, Watkins J, Noor AM, Williams C, Ali R, Sullivan R, et al. Economic downturns, universal health coverage, and cancer mortality in high-income and middle-income countries, 1990-2010: a longitudinal analysis. Lancet. 2016 Aug;388(10045):684-95. https://doi.org/10.1016/S0140-6736(16)00577-8
https://doi.org/10.1016/S0140-6736(16)00...
Besides, monitoring mortality trends is critical to understand recent changes in the etiology of OC and OPC (i.e., the increasing involvement of HPV virus as a risk factor for OPC).

In Brazil, the Mortality Information System (SIM) gathers the data from death records across the country since 1975/76.66. Pinto LF, Freitas MP, Figueiredo AW. National Information and Population Survey Systems: selected contributions from the Ministry of Health and the IBGE for analysis of Brazilian state capitals over the past 30 years. Cien Saúde Colet. 2018 Jun;23(6):1859-70. https://doi.org/10.1590/1413-81232018236.05072018
https://doi.org/10.1590/1413-81232018236...
In the last decades, researchers have analyzed the quality of the SIM, observing the two leading indicators of reliability of an information system: the coverage and completeness of the cause of death registration. The results show remarkable improvements in the quality of data in the last 20 years.77. Lima EE, Queiroz BL. Evolution of the deaths registry system in Brazil: associations with changes in the mortality profile, under-registration of death counts, and ill-defined causes of death. Cad Saúde Publica. 2014 Aug;30(8):1721-30. https://doi.org/10.1590/0102-311X00131113
https://doi.org/10.1590/0102-311X0013111...
,88. Mikkelsen L, Phillips DE, AbouZahr C, Setel PW, Savigny D, Lozano R, et al. A global assessment of civil registration and vital statistics systems: monitoring data quality and progress. Lancet. 2015 Oct;386(10001):1395-406. https://doi.org/10.1016/S0140-6736(15)60171-4
https://doi.org/10.1016/S0140-6736(15)60...
However, SIM reliability varies between regions of the country, and some consistency issues persist,99. Mathers CD, Fat DM, Inoue M, Rao C, Lopez AD. Counting the dead and what they died from: an assessment of the global status of cause of death data. Bull World Health Organ. 2005 Mar;83(3):171-7.,1010. Szwarcwald CL. Strategies for improving the monitoring of vital events in Brazil. Int J Epidemiol. 2008 Aug;37(4):738-44. https://doi.org/10.1093/ije/dyn130
https://doi.org/10.1093/ije/dyn130...
especially in the most impoverished areas.1111. Teixeira RA, Naghavi M, Guimarães MD, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol. 2019 Nov;22(22 Suppl 3):e19002. https://doi.org/10.1590/1980-549720190002.supl.3
https://doi.org/10.1590/1980-54972019000...

The development of strategies to deal with the reliability of mortality information systems is an international concern, and this is no different in Brazil.1212. Naghavi M, Makela S, Foreman K, O’Brien J, Pourmalek F, Lozano R. Algorithms for enhancing public health utility of national causes-of-death data. Popul Health Metr. 2010 May;8(1):9. https://doi.org/10.1186/1478-7954-8-9
https://doi.org/10.1186/1478-7954-8-9...
,1313. Szwarcwald CL, Morais-Neto OL, Frias PG, Souza Júnior PR, Escalante JJ, Lima RB, et al. Busca ativa de óbitos e nascimentos no Nordeste e na Amazônia Legal: estimação das coberturas do SIM e do Sinasc nos municípios brasileiros. In. Ministério da Saúde (BR). Saúde Brasil 2010: uma análise da situação de saúde e de evidências selecionadas de impacto de ações de vigilância em saúde. Brasília, DF: Ministério da Saúde; 2011. p. 79-98. In 2016, Brazil joined an international initiative to qualify information on causes of death, the Data for Health Initiative, along with 19 other countries. This collaboration triggered a series of SIM evaluation studies and interventions to qualify the records, mainly aimed at correcting “garbage codes” – according to Naghavi et al.,1212. Naghavi M, Makela S, Foreman K, O’Brien J, Pourmalek F, Lozano R. Algorithms for enhancing public health utility of national causes-of-death data. Popul Health Metr. 2010 May;8(1):9. https://doi.org/10.1186/1478-7954-8-9
https://doi.org/10.1186/1478-7954-8-9...
these are “all deaths assigned to codes that should be redistributed to enhance the validity of public health analysis”. However, these initiatives were primarily aimed at directly reclassifying ill-defined deaths, and although it produces highly reliable data, it is not an easily implemented alternative.1414. França EB. Garbage codes assigned as cause-of-death in health statistics. Rev Bras Epidemiol. 2019 Nov;22(3 Suppl 3):e19001. https://doi.org/10.1590/1980-549720190001.supl.3
https://doi.org/10.1590/1980-54972019000...

Various statistical approaches to dealing with this challenge by redistributing ill-defined deaths are currently proposed in the literature, and they differ substantially from one another.1515. França E, Teixeira R, Ishitani L, Duncan BB, Cortez-Escalante JJ, Morais Neto OL, et al. Ill-defined causes of death in Brazil: a redistribution method based on the investigation of such causes. Rev Saude Publica. 2014 Aug;48(4):671-81. https://doi.org/10.1590/S0034-8910.2014048005146
https://doi.org/10.1590/S0034-8910.20140...

16. Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012 Dec;380(9859):2095-128. https://doi.org/10.1016/S0140-6736(12)61728-0
https://doi.org/10.1016/S0140-6736(12)61...

17. World Health Organization. WHO methods and data sources for country-level causes of death 2000-2015. Geneva: World Health Organization; 2017. (Glob Heal Estim Tech Pap WHO/HIS/IER/GHE/20163).

18. Voss T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, 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. Lancet. 2020 Oct;396(10258):1204-22. https://doi.org/10.1016/S0140-6736(20)30925-9
https://doi.org/10.1016/S0140-6736(20)30...
-1919. Teixeira RA, Ishitani LH, Marinho F, Pinto Junior EP, Katikireddi SV, Malta DC. Proposta metodológica para redistribuição de óbitos por causas garbage nas estimativas de mortalidade para doenças crônicas não transmissíveis. Rev Bras Epidemiol. 2021;24:E210004.SUPL.1. https://doi.org/10.1590/1980-549720210004.supl.1
https://doi.org/10.1590/1980-54972021000...
Some methodologies assume that only deaths with ill-defined codes (Chapter XVIII of the ICD-10) are redistributed, while others extend redistribution to codes from other chapters – i.e., adopting the concept of “garbage codes”. The strategies for redistribution also vary in methodological complexity.2020. Bigoni A, Cunha AR, Antunes JL. Redistributing deaths by ill-defined and unspecified causes on cancer mortality in Brazil. Rev Saude Publica. 2021 Dec;55:106. https://doi.org/10.11606/s1518-8787.2021055003319
https://doi.org/10.11606/s1518-8787.2021...
These redistribution methods are usually used for epidemiological studies and are not directly applied to data in the information systems. However, influential epidemiological studies that feed publicly available databases, like the Global Burden of Disease Study, provide global estimates already corrected for garbage codes.1818. Voss T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, 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. Lancet. 2020 Oct;396(10258):1204-22. https://doi.org/10.1016/S0140-6736(20)30925-9
https://doi.org/10.1016/S0140-6736(20)30...
In Brazil, the SIM provides data with the registered cause of death, and the Brazilian Ministry of Health makes data available (by ICD-10 chapter) through an official website after applying a method for sub-record correction; however, data corrected for garbage codes are not available in official databases or websites.

Despite the interest in redistribution strategies, there is a lack of knowledge about how each correction strategy impacts results. In studies on OC and OPC mortality, it is not rare to find methodologies that entirely neglect the correction of ill-defined deaths, and, in most of them, correction procedures are applied without justifying their choice or discussing their impact. The objective of this study was to assess the impact on the magnitude and temporal trends of OC and OPC mortality of two statistical indirect approaches of redistribution of deaths with ill-defined underlying causes. The two techniques are presented in detail in the literature, which allows their replicability, and have been applied in cancer mortality studies worldwide and in Brazil.

Methodology

This study analyzed deaths with ill-defined causes that occurred in Brazil between 1996 and 2018. The SIM of Brazil, maintained by the Ministry of Health, provided official information on deaths, broken down by place of residence. The macro-region division is an official territorial organization of Brazil that distributes the twenty-six states and the Brazilian Federal District in five major regions: North, Northeast, Southeast, South, and Midwest.

The lower limit of the analyzed period is the year in which the SIM started using the International Classification of Diseases - tenth revision (ICD-10) for coding the underlying causes of death. We collected all deaths whose primary cause was oral cancer and oropharyngeal cancer, based on the classification proposed by Chaturvedi et al.2121. Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. 2008 Feb;26(4):612-9. https://doi.org/10.1200/JCO.2007.14.1713
https://doi.org/10.1200/JCO.2007.14.1713...
– the codes are detailed in Table 1. Deaths that did not have information on age, sex, or town of residence were allocated based on the known distribution of these characteristics by region.

Table 1
International Classification of Diseases, tenth revision (ICD-10) codes included in the categories of “oral cancer” and “oropharyngeal cancer”.

We obtained data on resident population from the Demographic Census carried out by the Brazilian Institute of Geography and Statistics in 1991, 2000, and 2010. Based on these data, we calculated within-census estimates using linear trends. Mortality rates due to oral or oropharyngeal cancer were calculated per 100,000 inhabitants for each macro-region. These rates were adjusted by sex and age (in 5-year intervals) using the direct method and the standard distribution of population defined by the World Health Organization.2222. Ahmad OB, Boschi-pinto C, Lopez AD. Age standardization of rates: a new WHO standard. GPE Discuss Pap Ser. 2001 [cited 2021 Feb 25];(31):1-14. Available from: http://www.who.int/healthinfo/paper31.pdf
http://www.who.int/healthinfo/paper31.pd...

Two different methods were applied to distribute deaths with ill-defined causes. The first one was based on the method proposed by França et al.,1515. França E, Teixeira R, Ishitani L, Duncan BB, Cortez-Escalante JJ, Morais Neto OL, et al. Ill-defined causes of death in Brazil: a redistribution method based on the investigation of such causes. Rev Saude Publica. 2014 Aug;48(4):671-81. https://doi.org/10.1590/S0034-8910.2014048005146
https://doi.org/10.1590/S0034-8910.20140...
the “EF method”. This technique proportionally reallocates deaths classified as R00-R99 of the ICD-10 – chapter XVIII – to the remaining specific causes of death according to the proportion of deaths with certified causes. However, it considers that deaths from cancer have a lower chance of being registered with “R” codes and proposes the redistribution of only half of the proportion that would be assigned to this disease. Following the description of Bigoni, Cunha, Antunes2020. Bigoni A, Cunha AR, Antunes JL. Redistributing deaths by ill-defined and unspecified causes on cancer mortality in Brazil. Rev Saude Publica. 2021 Dec;55:106. https://doi.org/10.11606/s1518-8787.2021055003319
https://doi.org/10.11606/s1518-8787.2021...
(2021), this technique can be summarized by the following formula:

0.5 N D r ( N D c A D N D r N E )

Where NDr is the number of deaths from the ICD-10 chapter XVIII; NDc is the certified number of deaths by oral or oropharyngeal cancer; AD is the total number of deaths; NE is the number of deaths by external causes – which would have been correctly recorded because autopsy is mandatory in Brazil for these deaths.

The second method for redistributing ill-defined deaths was based on the methodology of the 2010 Global Burden of Disease Study (GBD). In this technique (the “GBD method”), the concept of poorly registered deaths is expanded, and the terminology “garbage codes” is adopted.1616. Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012 Dec;380(9859):2095-128. https://doi.org/10.1016/S0140-6736(12)61728-0
https://doi.org/10.1016/S0140-6736(12)61...
This concept covers codes from several chapters of the ICD-10, and the GBD 2010 indicates which groups of garbage codes should be redistributed (and in what proportion) for each underlying cause of death studied. Considering only the application for cancer, the method can be summarized by the following formula, as proposed by Bigoni, Cunha and Antunes:2020. Bigoni A, Cunha AR, Antunes JL. Redistributing deaths by ill-defined and unspecified causes on cancer mortality in Brazil. Rev Saude Publica. 2021 Dec;55:106. https://doi.org/10.11606/s1518-8787.2021055003319
https://doi.org/10.11606/s1518-8787.2021...

N D c + N D g i C i

Where NDc is the number of deaths by oral cancer or oropharyngeal cancer (certified); i is the garbage code group; NDgi is the number of deaths by a garbage code group; and C is the coefficient, which indicates the proportion of the group that must be redistributed to cancer. The groups of garbage codes and the coefficients proposed to correct deaths from oral cancer and cancer of other parts of the pharynx and oropharynx are presented in the GBD 2010 publication.1616. Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012 Dec;380(9859):2095-128. https://doi.org/10.1016/S0140-6736(12)61728-0
https://doi.org/10.1016/S0140-6736(12)61...
For the correction of oral and oropharyngeal cancer cases only, we adapted this method by deducting the proportion of salivary gland and hypopharyngeal neoplasms from the coefficient.

For the trend analysis of mortality, a generalized linear regression using the Prais-Winsten method was carried out – the death rate logarithm was the dependent variable. This analysis was made for rates based on certified deaths and for rates corrected by the two redistribution methods, by macro-region. The Prais-Winsten method performs first-order autocorrelation correction: a common problem in the analysis of temporal series. It also quantifies the yearly rate change, allowing the calculation of the Annual Percent Change (APC) and its 95% confidence interval (CI95%):2323. Antunes JL, Waldman EA. Trends and spatial distribution of deaths of children aged 12-60 months in São Paulo, Brazil, 1980-98. Bull World Health Organ. 2002;80(5):391-8. https://doi.org/10.1590/S0042-96862002000500010
https://doi.org/10.1590/S0042-9686200200...

A P C = ( 1 + 10 b 1 ) 100 %
C I 95 % lower = ( 1 + 10 b1lower ) 100 %
C I 95 % upper = ( 1 + 10 b1upper ) 100 %

Where b1 represents the regression coefficient and b1lower/b1upper are its CI95% limits. Mortality was classified as increasing when APC and its CI95% were positive, decreasing, when they were negative, and stationary, when they included the zero.

The initial and final magnitudes of the rates were determined: the former corresponded to the average rate for 1996 to 1998 and the latter to the average rate for 2016 to 2018. Finally, we obtained the Rate Ratio (RR) between the certified and corrected rates (corrected rates divided by rates of certified data) and calculated the APC of RR.

Results

Death rates of oral and oropharyngeal cancer in Brazil from 1996 to 2018 had an evident impact after the redistribution of ill-defined deaths by the two methods; the increase in rates was higher using the GBD technique. The discrepancy between the certified and corrected rates was higher during the first half of the series and in the North and Northeast regions (Figures 1 and 2).

Figure 1
Oral cancer mortality, Brazil, 1996-2018. Certified and redistributed death rates by macro-region, both sexes, per 100,000 inhabitants.

Figure 2
Oropharyngeal cancer mortality, Brazil, 1996-2018. Certified and redistributed death rates by macro-region, both sexes, per 100,000 inhabitants.

Regarding oral cancer, the Southeast region had the highest initial (certified: 1.62/100,000 inhabitants; EF: 1.71; GBD: 2.34) and final (certified: 1.44; EF: 1.47; GBD: 1.90) death rates. Mortality underwent a decreasing trend in the Southeast in all scenarios and this was the only region where the trend in the three scenarios was the same. In general, trends assessed for corrected estimates had lower figures of APC. When the APC was positive, correction methods generated a less pronounced increase in rates than for non-corrected rates. When the APC was negative, correction methods generated a more pronounced decrease (Table 2).

Table 2
Certified versus redistributed mortality of oral and oropharyngeal cancer by the EF and GBD methods by macroregion, for 100,000 inhabitants: initial and final magnitude, annual percent change (APC), and trends. Brazil, 1996-2018.

The highest final (certified: 1.56; EF: 1.63; GBD: 2.20) and initial (certified: 1.47; EF: 1.50; GBD: 1.88) death rates of oropharyngeal cancer were identified in the Southeast region. In the South and Midwest regions, the trend identified by the GBD method was different from the others. Certified deaths in the Northeast had the highest APC (4.34%; increasing trend); mortality redistributed by the GBD method in the Southeast had the lowest APC (-0.84; decreasing trend). With few exceptions, the initial and final death rates of oral cancer were higher than those of oropharyngeal cancer (Table 2).

The correction by the GBD method generated the highest death rates. In the Northeast and North, for both types of cancer, this method more than doubled certified death rates (initial magnitude). The correction by the EF method had a lower impact on the magnitude of rates, only discretely increasing them. For the two types of cancer, both at the beginning and end of the study period, redistributing ill-defined deaths mainly impacted the North and Northeast regions by increasing the resulting rates. The rate ratio – RR comparing certified and redistributed death rates had a statistically significant decreasing trend for both methods and both types of cancer, in all regions (Table 3).

Table 3
Certified versus redistributed mortality for oral and oropharyngeal cancer by the EF and GBD methods by macroregion, for 100,000 inhabitants: Rate Ratio (RR) of the final and initial magnitudes (redistributed/certified), annual percent change (APC), and trend of RR. Brazil, 1996-2018.

Discussion

This study assessed the application of two indirect methods of correcting ill-defined deaths from OC and OPC in Brazil from 1996 to 2018. In comparison with certified deaths, both methods – EF and GBD – resulted in higher rates and, in some regions, they even changed the rate’s temporal trend. The impact of the correction method was different, indicating that this is not a mere methodological formality. In countries like Brazil, where there’s no nationwide information system on cancer incidence, it is even more important to recognize the difference between confirmed and corrected rates; it is mainly the mortality data that allow us to understand the dynamics of this disease. As far as we know, this is the first study that analyzes techniques for correcting ill-defined deaths for OC and OPC specifically.

The two methods we applied have essential methodological differences. The EF method allows the redistribution of ill-defined deaths, which are those registered with codes from chapter XVIII of ICD-10 (the “R” codes), while the GBD technique redistributes a range of other causes defined as “garbage codes”, in addition to the R codes. This category covers deaths registered with codes that should not be listed as the primary cause of death because a) they are immediate or intermediate causes of the disease chain that led to death or b) they are causes that do not lead to death. This is the probable reason why the EF technique’s results were more conservative than those of the GBD: the concept of ill-defined deaths is less comprehensive. Considering these differences, it is essential to highlight that in Brazil, 5.4% of deaths registered in the SIM in 2015 had an R code as the underlying cause, while 26.34% were registered with non-R garbage codes.1111. Teixeira RA, Naghavi M, Guimarães MD, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol. 2019 Nov;22(22 Suppl 3):e19002. https://doi.org/10.1590/1980-549720190002.supl.3
https://doi.org/10.1590/1980-54972019000...

Considering the EF technique, the percentage that we redistributed from ill-defined deaths to deaths by OC and OPC corresponds to an estimate for all cancers because unlike the GBD method, this method does not propose coefficients by cancer type. Almost half of the OC and OPC cases are diagnosed in advanced stages – stages III or IV2424. McGurk M, Chan C, Jones J, O’regan E, Sherriff M. Delay in diagnosis and its effect on outcome in head and neck cancer. Br J Oral Maxillofac Surg. 2005 Aug;43(4):281-4. https://doi.org/10.1016/j.bjoms.2004.01.016
https://doi.org/10.1016/j.bjoms.2004.01....
. Large tumors may spread from their primary site, thus preventing their identification. Considering this particularity, perhaps a different proportion of poorly registered deaths should be redistributed to OC and OPC; this proportion would be larger than that defined for all types of cancer by the EF technique. This assumption is anchored in the evidence that “other and unspecified parts of the mouth” – which is the most nonspecific code related to oral topographies – is the category with the highest proportion of mortality among oral cancers in Brazil.2525. Cunha AR, Prass TS, Hugo FN. Mortality from oral and oropharyngeal cancer in Brazil: impact of the National Oral Health Policy. Cad Saude Publica. 2019 Nov;35(12):e00014319. https://doi.org/10.1590/0102-311x00014319
https://doi.org/10.1590/0102-311x0001431...

The impact of redistribution methods was also observed for the APC, which is a significant finding because monitoring incidence and mortality trends is a strategic tool for OC and OPC’s surveillance related to risk factors.2626. Menezes FD, Fernandes GA, Antunes JL, Villa LL, Toporcov TN. Global incidence trends in head and neck cancer for HPV-related and -unrelated subsites: A systematic review of population-based studies. Oral Oncol. 2021 Apr;115:105177. https://doi.org/10.1016/j.oraloncology.2020.105177
https://doi.org/10.1016/j.oraloncology.2...
With the absence of population-based cancer records that inform the tumor’s etiology and the increased involvement of the HPV virus as a risk factor for OPC, comparing the populational behavior of OC (which is site not related to HPV) with that of OPC (HPV-related) provides clues about their etiology.2626. Menezes FD, Fernandes GA, Antunes JL, Villa LL, Toporcov TN. Global incidence trends in head and neck cancer for HPV-related and -unrelated subsites: A systematic review of population-based studies. Oral Oncol. 2021 Apr;115:105177. https://doi.org/10.1016/j.oraloncology.2020.105177
https://doi.org/10.1016/j.oraloncology.2...

27. Chaturvedi AK, Engels EA, Pfeiffer RM, Hernandez BY, Xiao W, Kim E, et al. Human papillomavirus and rising oropharyngeal cancer incidence in the United States. J Clin Oncol. 2011 Nov;29(32):4294-301. https://doi.org/10.1200/JCO.2011.36.4596
https://doi.org/10.1200/JCO.2011.36.4596...
-2828. Conway DI, Purkayastha M, Chestnutt IG. The changing epidemiology of oral cancer: definitions, trends, and risk factors. Br Dent J. 2018 Nov;225(9):867-73. https://doi.org/10.1038/sj.bdj.2018.922
https://doi.org/10.1038/sj.bdj.2018.922...
This comparison serves as an alternative assessment of the virus’s potential carcinogenic effect.2121. Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. 2008 Feb;26(4):612-9. https://doi.org/10.1200/JCO.2007.14.1713
https://doi.org/10.1200/JCO.2007.14.1713...
Our findings indicate that analyses based only on death certificate data can mask trends and lead to mistaken etiological hypotheses.

Deaths with ill-defined primary cause are more frequent among older adults than in the other age groups.77. Lima EE, Queiroz BL. Evolution of the deaths registry system in Brazil: associations with changes in the mortality profile, under-registration of death counts, and ill-defined causes of death. Cad Saúde Publica. 2014 Aug;30(8):1721-30. https://doi.org/10.1590/0102-311X00131113
https://doi.org/10.1590/0102-311X0013111...
,1111. Teixeira RA, Naghavi M, Guimarães MD, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol. 2019 Nov;22(22 Suppl 3):e19002. https://doi.org/10.1590/1980-549720190002.supl.3
https://doi.org/10.1590/1980-54972019000...
Considering that the correction methods used general coefficients to redistribute ill-defined deaths (i.e., the same for different age groups), they may be underestimating the number of deaths that should be redistributed to older adults. Since OC and OPC deaths also mostly affect older adults, this inaccuracy may be relevant in the context of our results, meaning that mortality might be even higher than what has been obtained with the correction of estimates by the application of the EF and GBD methods.

In the initial part of the study period, the GBD method more than doubled mortality from oral cancer in the North and Northeast, while the correction’s impact was less pronounced in the other regions, which indicated a higher proportion of garbage codes in those two regions. A recent study that corrected deaths from noncommunicable diseases in Brazilian municipalities also identified clusters with higher proportions of garbage codes in the North and Northeast1919. Teixeira RA, Ishitani LH, Marinho F, Pinto Junior EP, Katikireddi SV, Malta DC. Proposta metodológica para redistribuição de óbitos por causas garbage nas estimativas de mortalidade para doenças crônicas não transmissíveis. Rev Bras Epidemiol. 2021;24:E210004.SUPL.1. https://doi.org/10.1590/1980-549720210004.supl.1
https://doi.org/10.1590/1980-54972021000...
. However, in the final period of the study, the difference between corrected and registered deaths was more homogeneous between regions, indicating that these two regions partially compensated for their disadvantage in the quality of death registration data. For oral cancer, correction methods changed the time trend of rates in the North, Northeast, South, and Midwest regions. This means that the qualification of death registries in this historical series has been masking the real mortality trend from oral cancer in these four regions. The results from the Southeast, which were the exception in this analysis, may be reflecting two phenomena: a) lower impact of redistribution methods on temporal trends due to the lower proportion of garbage codes since the start of the series – the results expressed in Table 2 indicated that the Southeast region presented the smallest difference between the certified and corrected rates at the beginning of the series; and b) a stronger real temporal trend (decrease) than in the other regions, which not even the qualification of the registry was able to mask.

The results of the rate ratio that compares certified and redistributed death rates indicate a relevant reduction in the number of ill-defined deaths in the country between 1996 and 2018, regardless of the correction technique. This result is in line with other evidence of qualification of death records in Brazil in recent decades,1111. Teixeira RA, Naghavi M, Guimarães MD, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol. 2019 Nov;22(22 Suppl 3):e19002. https://doi.org/10.1590/1980-549720190002.supl.3
https://doi.org/10.1590/1980-54972019000...
and it does not appear to be due to chance. Several initiatives aimed at qualifying the death registry worldwide and in Brazil are presented in the literature, ranging from advanced methods for assessing the problem to direct training of registrars, supported by international collaboration.1414. França EB. Garbage codes assigned as cause-of-death in health statistics. Rev Bras Epidemiol. 2019 Nov;22(3 Suppl 3):e19001. https://doi.org/10.1590/1980-549720190001.supl.3
https://doi.org/10.1590/1980-54972019000...
,2929. Phillips DE, Lozano R, Naghavi M, Atkinson C, Gonzalez-Medina D, Mikkelsen L, et al. A composite metric for assessing data on mortality and causes of death: the vital statistics performance index. Popul Health Metr. 2014 May;12(1):14. https://doi.org/10.1186/1478-7954-12-14
https://doi.org/10.1186/1478-7954-12-14...
,3030. Mikkelsen L, Moesgaard K, Hegnauer M, Lopez AD. ANACONDA: a new tool to improve mortality and cause of death data. BMC Med. 2020 Mar;18(1):61. https://doi.org/10.1186/s12916-020-01521-0
https://doi.org/10.1186/s12916-020-01521...
The reduction in the number of poorly registered deaths was more significant in the country’s poorest regions, in which death records are historically more precarious (North and Northeast), suggesting that these initiatives are reaching the weakest points of the system.

We adopted the Chaturvedi classification to define which ICD-10 codes should comprise the OC and OPC2121. Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. 2008 Feb;26(4):612-9. https://doi.org/10.1200/JCO.2007.14.1713
https://doi.org/10.1200/JCO.2007.14.1713...
categories, which is a limitation of this study because it is not the only classification internationally accepted2828. Conway DI, Purkayastha M, Chestnutt IG. The changing epidemiology of oral cancer: definitions, trends, and risk factors. Br Dent J. 2018 Nov;225(9):867-73. https://doi.org/10.1038/sj.bdj.2018.922
https://doi.org/10.1038/sj.bdj.2018.922...
. Also, in the GBD method, we employed the correction technique used in the GBD 2010 Study. Since then, the Institute of Health Metrics and Evaluation has developed more sophisticated versions of the statistical corrections for garbage codes.1818. Voss T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, 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. Lancet. 2020 Oct;396(10258):1204-22. https://doi.org/10.1016/S0140-6736(20)30925-9
https://doi.org/10.1016/S0140-6736(20)30...
Our rationale for choosing this less recent technique is its applicability: the new versions require advanced technical expertise and make reproducibility difficult for researchers external to the project.

This study does not allow us to define which of the two techniques is more accurate. For this, they would have to be compared to a gold standard method, in which each ill-defined death is investigated through analysis of health service records or verbal autopsy. Brazilian researchers developed a study with a methodology similar to that used for infant mortality, which produced a correction method that the Brazilian Ministry of Health adopted as a correction strategy.1313. Szwarcwald CL, Morais-Neto OL, Frias PG, Souza Júnior PR, Escalante JJ, Lima RB, et al. Busca ativa de óbitos e nascimentos no Nordeste e na Amazônia Legal: estimação das coberturas do SIM e do Sinasc nos municípios brasileiros. In. Ministério da Saúde (BR). Saúde Brasil 2010: uma análise da situação de saúde e de evidências selecionadas de impacto de ações de vigilância em saúde. Brasília, DF: Ministério da Saúde; 2011. p. 79-98. The results of applying this method are available in the country’s official databases; however, they do not provide corrections by cancer subtype, but only by ICD-10 chapter, making the comparison with the data of the present study not feasible. We emphasize that indirect correction methods are also imperfect and, although more feasible than direct reclassifications, are based on generalized concepts that can result in inaccuracies. The available redistribution methods cannot replace registry qualification efforts.

Conclusions

The present study suggests that mortality by OC and OPC is higher than informed by Brazil’s Mortality Information System. Even considering the rates of the most recent period (2016 to 2018), a significant portion of deaths due to OC and OPC is hidden by low quality death certification. This problem is probably not restricted to Brazil; it is the likely scenario in all countries with a non-negligible proportion of poorly registered deaths. This study concludes that, in epidemiological studies of OC and OPC mortality, it is necessary to redistribute deaths from ill-defined causes when analyzing data from less-than-optimal information systems. The choice of correction method is critical and must be managed in epidemiological studies as a methodological decision that may have a significant impact on results. Researchers should justify, discuss, and, if possible, even compare correction methods before presenting mortality data for OC and OPC to avoid misrepresenting the burden of this disease and the challenges revealed by data analysis.

Acknowledgements

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, Brazil) – Finance Code 001.

References

  • 1
    Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, Abdel-Rahman O, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-Adjusted life-years for 29 cancer groups, 1990 to 2017: A systematic analysis for the global burden of disease study. JAMA Oncol. 2019 Dec;5(12):1749-68. https://doi.org/10.1001/jamaoncol.2019.2996
    » https://doi.org/10.1001/jamaoncol.2019.2996
  • 2
    Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Piñeros M, et al. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer. 2019 Apr;144(8):1941-53. https://doi.org/10.1002/ijc.31937
    » https://doi.org/10.1002/ijc.31937
  • 3
    International Agency for Research on Cancer. Cancer today: data visualization tool for exploring the global cancer burden in 2020 [Internet]. 2021 [cited 2021 Feb 25]. Available from: https://gco.iarc.fr/today/home
    » https://gco.iarc.fr/today/home
  • 4
    World Health Organization. Oral health: achieving better oral health as part of the universal health coverage and noncommunicable disease agendas towards 2030. Geneva: World Health Organization; 2020 [cited 2021 Feb 25]. Available from: https://apps.who.int/gb/ebwha/pdf_files/EB148/B148_8-en.pdf
    » https://apps.who.int/gb/ebwha/pdf_files/EB148/B148_8-en.pdf
  • 5
    Maruthappu M, Watkins J, Noor AM, Williams C, Ali R, Sullivan R, et al. Economic downturns, universal health coverage, and cancer mortality in high-income and middle-income countries, 1990-2010: a longitudinal analysis. Lancet. 2016 Aug;388(10045):684-95. https://doi.org/10.1016/S0140-6736(16)00577-8
    » https://doi.org/10.1016/S0140-6736(16)00577-8
  • 6
    Pinto LF, Freitas MP, Figueiredo AW. National Information and Population Survey Systems: selected contributions from the Ministry of Health and the IBGE for analysis of Brazilian state capitals over the past 30 years. Cien Saúde Colet. 2018 Jun;23(6):1859-70. https://doi.org/10.1590/1413-81232018236.05072018
    » https://doi.org/10.1590/1413-81232018236.05072018
  • 7
    Lima EE, Queiroz BL. Evolution of the deaths registry system in Brazil: associations with changes in the mortality profile, under-registration of death counts, and ill-defined causes of death. Cad Saúde Publica. 2014 Aug;30(8):1721-30. https://doi.org/10.1590/0102-311X00131113
    » https://doi.org/10.1590/0102-311X00131113
  • 8
    Mikkelsen L, Phillips DE, AbouZahr C, Setel PW, Savigny D, Lozano R, et al. A global assessment of civil registration and vital statistics systems: monitoring data quality and progress. Lancet. 2015 Oct;386(10001):1395-406. https://doi.org/10.1016/S0140-6736(15)60171-4
    » https://doi.org/10.1016/S0140-6736(15)60171-4
  • 9
    Mathers CD, Fat DM, Inoue M, Rao C, Lopez AD. Counting the dead and what they died from: an assessment of the global status of cause of death data. Bull World Health Organ. 2005 Mar;83(3):171-7.
  • 10
    Szwarcwald CL. Strategies for improving the monitoring of vital events in Brazil. Int J Epidemiol. 2008 Aug;37(4):738-44. https://doi.org/10.1093/ije/dyn130
    » https://doi.org/10.1093/ije/dyn130
  • 11
    Teixeira RA, Naghavi M, Guimarães MD, Ishitani LH, França EB. Quality of cause-of-death data in Brazil: garbage codes among registered deaths in 2000 and 2015. Rev Bras Epidemiol. 2019 Nov;22(22 Suppl 3):e19002. https://doi.org/10.1590/1980-549720190002.supl.3
    » https://doi.org/10.1590/1980-549720190002.supl.3
  • 12
    Naghavi M, Makela S, Foreman K, O’Brien J, Pourmalek F, Lozano R. Algorithms for enhancing public health utility of national causes-of-death data. Popul Health Metr. 2010 May;8(1):9. https://doi.org/10.1186/1478-7954-8-9
    » https://doi.org/10.1186/1478-7954-8-9
  • 13
    Szwarcwald CL, Morais-Neto OL, Frias PG, Souza Júnior PR, Escalante JJ, Lima RB, et al. Busca ativa de óbitos e nascimentos no Nordeste e na Amazônia Legal: estimação das coberturas do SIM e do Sinasc nos municípios brasileiros. In. Ministério da Saúde (BR). Saúde Brasil 2010: uma análise da situação de saúde e de evidências selecionadas de impacto de ações de vigilância em saúde. Brasília, DF: Ministério da Saúde; 2011. p. 79-98.
  • 14
    França EB. Garbage codes assigned as cause-of-death in health statistics. Rev Bras Epidemiol. 2019 Nov;22(3 Suppl 3):e19001. https://doi.org/10.1590/1980-549720190001.supl.3
    » https://doi.org/10.1590/1980-549720190001.supl.3
  • 15
    França E, Teixeira R, Ishitani L, Duncan BB, Cortez-Escalante JJ, Morais Neto OL, et al. Ill-defined causes of death in Brazil: a redistribution method based on the investigation of such causes. Rev Saude Publica. 2014 Aug;48(4):671-81. https://doi.org/10.1590/S0034-8910.2014048005146
    » https://doi.org/10.1590/S0034-8910.2014048005146
  • 16
    Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012 Dec;380(9859):2095-128. https://doi.org/10.1016/S0140-6736(12)61728-0
    » https://doi.org/10.1016/S0140-6736(12)61728-0
  • 17
    World Health Organization. WHO methods and data sources for country-level causes of death 2000-2015. Geneva: World Health Organization; 2017. (Glob Heal Estim Tech Pap WHO/HIS/IER/GHE/20163).
  • 18
    Voss T, Lim SS, Abbafati C, Abbas KM, Abbasi M, Abbasifard M, 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. Lancet. 2020 Oct;396(10258):1204-22. https://doi.org/10.1016/S0140-6736(20)30925-9
    » https://doi.org/10.1016/S0140-6736(20)30925-9
  • 19
    Teixeira RA, Ishitani LH, Marinho F, Pinto Junior EP, Katikireddi SV, Malta DC. Proposta metodológica para redistribuição de óbitos por causas garbage nas estimativas de mortalidade para doenças crônicas não transmissíveis. Rev Bras Epidemiol. 2021;24:E210004.SUPL.1. https://doi.org/10.1590/1980-549720210004.supl.1
    » https://doi.org/10.1590/1980-549720210004.supl.1
  • 20
    Bigoni A, Cunha AR, Antunes JL. Redistributing deaths by ill-defined and unspecified causes on cancer mortality in Brazil. Rev Saude Publica. 2021 Dec;55:106. https://doi.org/10.11606/s1518-8787.2021055003319
    » https://doi.org/10.11606/s1518-8787.2021055003319
  • 21
    Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. 2008 Feb;26(4):612-9. https://doi.org/10.1200/JCO.2007.14.1713
    » https://doi.org/10.1200/JCO.2007.14.1713
  • 22
    Ahmad OB, Boschi-pinto C, Lopez AD. Age standardization of rates: a new WHO standard. GPE Discuss Pap Ser. 2001 [cited 2021 Feb 25];(31):1-14. Available from: http://www.who.int/healthinfo/paper31.pdf
    » http://www.who.int/healthinfo/paper31.pdf
  • 23
    Antunes JL, Waldman EA. Trends and spatial distribution of deaths of children aged 12-60 months in São Paulo, Brazil, 1980-98. Bull World Health Organ. 2002;80(5):391-8. https://doi.org/10.1590/S0042-96862002000500010
    » https://doi.org/10.1590/S0042-96862002000500010
  • 24
    McGurk M, Chan C, Jones J, O’regan E, Sherriff M. Delay in diagnosis and its effect on outcome in head and neck cancer. Br J Oral Maxillofac Surg. 2005 Aug;43(4):281-4. https://doi.org/10.1016/j.bjoms.2004.01.016
    » https://doi.org/10.1016/j.bjoms.2004.01.016
  • 25
    Cunha AR, Prass TS, Hugo FN. Mortality from oral and oropharyngeal cancer in Brazil: impact of the National Oral Health Policy. Cad Saude Publica. 2019 Nov;35(12):e00014319. https://doi.org/10.1590/0102-311x00014319
    » https://doi.org/10.1590/0102-311x00014319
  • 26
    Menezes FD, Fernandes GA, Antunes JL, Villa LL, Toporcov TN. Global incidence trends in head and neck cancer for HPV-related and -unrelated subsites: A systematic review of population-based studies. Oral Oncol. 2021 Apr;115:105177. https://doi.org/10.1016/j.oraloncology.2020.105177
    » https://doi.org/10.1016/j.oraloncology.2020.105177
  • 27
    Chaturvedi AK, Engels EA, Pfeiffer RM, Hernandez BY, Xiao W, Kim E, et al. Human papillomavirus and rising oropharyngeal cancer incidence in the United States. J Clin Oncol. 2011 Nov;29(32):4294-301. https://doi.org/10.1200/JCO.2011.36.4596
    » https://doi.org/10.1200/JCO.2011.36.4596
  • 28
    Conway DI, Purkayastha M, Chestnutt IG. The changing epidemiology of oral cancer: definitions, trends, and risk factors. Br Dent J. 2018 Nov;225(9):867-73. https://doi.org/10.1038/sj.bdj.2018.922
    » https://doi.org/10.1038/sj.bdj.2018.922
  • 29
    Phillips DE, Lozano R, Naghavi M, Atkinson C, Gonzalez-Medina D, Mikkelsen L, et al. A composite metric for assessing data on mortality and causes of death: the vital statistics performance index. Popul Health Metr. 2014 May;12(1):14. https://doi.org/10.1186/1478-7954-12-14
    » https://doi.org/10.1186/1478-7954-12-14
  • 30
    Mikkelsen L, Moesgaard K, Hegnauer M, Lopez AD. ANACONDA: a new tool to improve mortality and cause of death data. BMC Med. 2020 Mar;18(1):61. https://doi.org/10.1186/s12916-020-01521-0
    » https://doi.org/10.1186/s12916-020-01521-0

Publication Dates

  • Publication in this collection
    10 Oct 2022
  • Date of issue
    2022

History

  • Received
    3 Sept 2021
  • Accepted
    2 May 2022
  • Reviewed
    1 June 2022
Sociedade Brasileira de Pesquisa Odontológica - SBPqO Av. Prof. Lineu Prestes, 2227, 05508-000 São Paulo SP - Brazil, Tel. (55 11) 3044-2393/(55 11) 9-7557-1244 - São Paulo - SP - Brazil
E-mail: office.bor@ingroup.srv.br