Mortality from idiopathic pulmonary fibrosis: a temporal trend analysis in Brazil, 1979-2014.

OBJECTIVE: To analyze mortality from idiopathic pulmonary fibrosis (IPF) in Brazil over the period 1979-2014.
METHODS: Microdata were extracted from the Brazilian National Ministry of Health Mortality Database. Only deaths for which the underlying cause was coded as International Classification of Diseases version 9 (ICD-9) 515 or 516.3 (until 1995) or as ICD version 10 (ICD-10) J84.1 (from 1996 onward) were included in our analysis. Standardized mortality rates were calculated for the 2010 Brazilian population. The annual trend in mortality rates was analyzed by joinpoint regression. We calculated risk ratios (RRs) by age group, time period of death, and gender, using a person-years denominator.
RESULTS: A total of 32,092 deaths were recorded in the study period. Standardized mortality rates trended upward, rising from 0.24/100,000 population in 1979 to 1.10/100,000 population in 2014. The annual upward trend in mortality rates had two inflection points, in 1992 and 2008, separating three distinct time segments with an annual growth of 2.2%, 6.8%, and 2.4%, respectively. The comparison of RRs for the age groups, using the 50- to 54-year age group as a reference, and for the study period, using 1979-1984 as a reference, were 16.14 (14.44-16.36) and 6.71 (6.34-7.12), respectively. Men compared with women had higher standardized mortality rates (per 100,000 person-years) in all age groups.
CONCLUSION: Brazilian IPF mortality rates are lower than those of other countries, suggesting underdiagnosis or underreporting. The temporal trend is similar to those reported in the literature and is not explained solely by population aging.

INTRODUCTION

Idiopathic pulmonary fibrosis (IPF) is a disease of unknown cause, predominantly occurring in individuals over 60 years of age and usually exhibiting rapid clinical progression and high lethality, with a median survival of less than 5 years. Estimates of prevalence and incidence vary according to the data sources used and to the broadness or narrowness of the criteria for disease definition. - It is known that mortality from IPF has been increasing worldwide, ranging from 4 to 10 deaths per 100,000 population in ten countries in Europe, North America, Asia, and Oceania.

Studying the rates of occurrence of IPF is an aiding tool in developing appropriate health policies and in planning IPF management. Recently, the availability of antifibrotic drugs with the potential to slow down the functional decline associated with IPF has sharpened the need for precision, on the part of specialists, in identifying the disease, given the possibility of a change in the clinical approach. , Because the symptoms and imaging findings of IPF have similarities with those of other interstitial lung diseases, added weight has been given to the differential diagnosis, since other diseases of similar presentation progress differently, require different approaches, and differ in prognosis.

National data reveal that, in the Brazilian state of Rio Grande do Sul, mortality from IPF increased significantly between 1970 and 2000, rising from 0.22/100,000 population to 0.48/100,000 population during that period. In Brazil, another analysis of mortality showed that there was an increase from 0.65/100,000 deaths in 1996 to 1.21/100,000 deaths in 2010.

The objective of the present study was to investigate mortality from IPF as the underlying cause and examine its temporal trend, on the basis of deaths occurring between 1979 and 2014 in Brazil.

METHODS

Study design

This was an ecological study of temporal trends in mortality from IPF as the underlying cause, involving analysis of aggregate data from the Brazilian National Ministry of Health Sistema de Informações de Mortalidade (SIM, Mortality Database), for the period 1979 to 2014. Because of possible inclusion of inaccurate data in the SIM (deaths from other interstitial lung diseases, here referred to as “similars”), cases were defined as deaths for which the underlying cause was attributed to IPF + similars (IPF-S).

Microdata were extracted from the SIM, in an anonymized fashion, for the period January 1, 1979 to December 31, 2014. During that period, two versions of the International Classification of Diseases (ICD) were used: ICD version 9 (ICD-9), from 1979 to 1995; and ICD version 10 (ICD-10), from 1996 onward. Deaths occurring until December 31, 1995 for which the underlying cause was coded as ICD-9 515 (post-inflammatory pulmonary fibrosis) or 516.3 (idiopathic pulmonary fibrosis) in the SIM were included in our analysis. For deaths occurring from January 1, 1996 onward, code J84.1 (other interstitial lung diseases with fibrosis) was used, assuming that the diseases were equivalent to IPF. For each unit of observation, information on year of death, age, and gender was collected.

Data analysis

Annual crude mortality rates were calculated using the population of each year, obtained through census data for years 1980, 1991, 2000, and 2010 or from estimates for intercensal years provided by the Brazilian Institute of Geography and Statistics. The annual crude mortality rates were adjusted to the 2010 census in order to obtain standardized mortality rates per 100,000 population.

The annual percent change (APC) in the rates was calculated using joinpoint regression (Joinpoint Regression Program; Statistical Methodology and Applications Branch and Data Modeling Branch, Surveillance Research Program, National Cancer Institute, Rockville, MD, USA). The inflection points indicate statistically significant changes in the curve.

Mortality was also analyzed by obtaining standardized mortality rates per 100,000 person-years, which were calculated for three variables: a) age-first age group: 0-49 years, followed by groups covering a 5-year age range; b) time-˜5-year time periods of death, with no combination of the period preceding and that following the ICD change; and c) gender. For calculating person-years, the 2010 Brazilian population was used as a reference. For age, the number of person-years was estimated by multiplying the population within each age group by the midpoint of that range. Confidence intervals were calculated assuming a Poisson distribution and using Byar approximation. Risk ratios (RRs) were calculated using the 50- to 54-year age group, the first time period studied (1979-1984), and the male gender as a reference. The corresponding confidence intervals were calculated using Taylor series. These calculations were performed using OpenEpi, version 3.01.

Ethical aspects

Because the present study used anonymized secondary public domain data, available, with unrestricted access, in government digital databases, it did not require approval from an ethics committee.

RESULTS

Between 1979 and 2014, there were 32,092 deaths for which the underlying cause was IPF-S. Of those, 3,169 (9.9%) were deaths of individuals up to 49 years of age. Regarding gender, there were 15,782 and 16,302 deaths of men and women, respectively. In 8 cases, there was no information about gender. Among individuals 75 years of age or older, the absolute number of deaths of women was greater than that of men.

Crude and standardized mortality rates trended upward, with the latter rising from 0.24/100,000 population in 1979 to 1.10/100,000 population in 2014. Figure 1 shows the temporal variation over the study period. There were no differences in the pattern of increase in mortality rates when women and men were analyzed separately.

Figure 1

Crude and standardized mortality rates from idiopathic pulmonary fibrosis + similars, 1979-2014, Brazil. pop.: population.

The temporal trend in standardized mortality rates had two inflection points separating three distinct time segments: 1979-1992; 1992-2008; and 2008-2014. Figure 2 shows the time segments and their corresponding APCs. There was a significant increase in mortality throughout the time period analyzed; however, in the years 1992 to 2008, there was an average annual growth of 6.8%.

Figure 2

Annual trend in standardized mortality rates from idiopathic pulmonary fibrosis + similars, as determined by joinpoint regression, 1979-2014, Brazil. pop.: population; APC: annual percent change. The APC was significantly different from zero at α = 0.05 in all periods.

The annual increase in APCs differed slightly between genders, and the trends for both genders had two inflection points (Table 1). The greatest increase in the rates occurred between 1992 and 2008 in men and between 1992 and 2007 in women.

Table 1

Trends in standardized mortality rates from idiopathic pulmonary fibrosis + similars, by gender, as determined by joinpoint regression, 1979-2014, Brazil.

Gender	Trend 1		Trend 2		Trend 3
	Period	APC (95% CI)	Period	APC (95% CI)	Period	APC (95% CI)
Male	1979-1992	2.0 (1.6-2.4)	1992-2008	6.4 (6.1-6.8)	2008-2014	1.6 (0.5-2.9)
Female	1979-1992	2.1 (1.8-2.4)	1992-2007	7.3 (7.0-7.6)	2007-2014	3.6 (2.7-4.5)

APC: annual percent change.

Analysis of census population changes occurring between 1980 and 2010, during which time there were four demographic censuses, revealed a 68% increase in the population 50 years of age or older, which increased by 61% among men and by 74% among women. In that same period, the increase in standardized mortality rates from IPF-S among individuals 50 years of age or older was 317%.

The RRs, calculated from mortality rates per person-years, are shown in Table 2. There was a strong upward trend in mortality as the age group increased and by time period of death. The rates, adjusted by person-years, showed that mortality increased by a factor of more than 16 when comparing the age groups, using the 50- to 54-year age group as a reference, and by a factor of more than 6, when comparing the time periods of death, using the first time period studied as a reference. There were no gender-related differences in overall standardized mortality rates (95% CI: 0.99-1.04).

Table 2

Mortality risk ratio from idiopathic pulmonary fibrosis + similars, by age group, time period of death, and gender, 1979-2014, Brazil.

Variable	Deaths	Person-years	SMR per 100,000 person-years (95% CI)*	RR (95% CI)**
Age group, years
< 50	3,169	3.793.719.975	0.08 (0.06-0.09)	0.32 (0.29-0.37)
50-54	1,370	527.300.904	0.26 (0.24-0.27)	1.00
55-59	2,056	471.744.483	0.44 (0.41-0.45)	1.68 (1.55-1.79)
60-64	2,880	403.565.378	0.71 (0.68-0.74)	2.75 (2.53-2.89)
65-69	3,644	324.334.270	1.12 (1.09-1.16)	4.32 (4.01-4.55)
70-74	4,688	269.397.864	1.74 (1.69-1.79)	6.70 (6.18-7.00)
75-79	5,142	197.385.496	2.61 (2.54-2.68)	10.03 (9.32-10.54)
80-84	4,356	136.691.704	3.19 (3.09-3.28)	12.27 (11.3-12.81)
85 or older	4,787	114.175.170	4.19 (4.07-4.31)	16.14 (14.44-16.36)
Time period of death
1979-1984	1,308	190.755.799	0.69 (0.64-0.71)	1.00
1985-1990	1,904	190.755.799	1.00 (0.95-1.03)	1.46 (1.36-1.56)
1991-1995	2,147	190.755.799	1.13 (1.09-1.16)	1.64 (1.53-1.75)
1996-2000	3,605	190.755.799	1.89 (1.87-1.91)	2.76 (2.58-2.94)
2001-2005	5,641	190.755.799	2.96 (2.95-3.00)	4.31 (4.07-4.58)
2006-2010	8,704	190.755.799	4.56 (4.48-4.61)	6.65 (6.28-7.05)
2011-2014	8,783	190.755.799	4.60 (4.58-4.63)	6.71 (6.34-7.12)
Gendera
Male	15,782	1.112.062.549	1.42 (1.38-1.45)	1.00
Female	16,302	1.128.472.320	1.44 (1.40-1.46)	1.02 (0.99-1.04)

SMR: standardized mortality rate; and RR: risk ratio. aEight cases lacked information on gender.

*95% CIs were calculated assuming a Poisson distribution and using Byar approximation. &95% CIs were calculated using Taylor series.

Analysis by gender and age group showed that the standardized mortality rates per 100,000 person-years were higher in men, with the differences in RR between men and women being statistically significant in all age groups, except in the 85-year or older age group (Figure 3).

Figure 3

Standardized mortality rates per 100,000 person-years, by gender and age group, 1979-2014, Brazil.

DISCUSSION

The data obtained from the Brazilian National Ministry of Health SIM showed a consistent increase in mortality from IPF-S as the underlying cause over the period 1979-2014. During that period, the standardized mortality rates increased 4.6 times, rising from 0.24 to 1.10 per 100,000 population. The upward trend in mortality from IPF had been described in the Brazilian state of Rio Grande do Sul for the period 1970-2000 and in Brazil for the period 1996-2010. Possibly, the small difference between the study by Rufino et al. and ours is due to the reference population used for standardization of rates. In the present study, we did not include any mention of IPF on other lines of the death certificate because this information was unavailable in the SIM until 1999.

In our study, analysis of the annual trend in mortality revealed a marked annual increase, with the average annual growth being strongest between 1992 and 2008 (6.8%). In comparison, an analysis of mortality in ten countries revealed an average annual growth of 2-3% between 1999 and 2012. The change from ICD-9 to ICD-10 in 1996 was not accompanied by a change in curve inflection, which demonstrates that there was no detectable impact on the reporting of deaths.

The increase in mortality rates, as well as the existence of a period when growth was most marked (1992-2008), is possibly due to a combination of factors. Technical issues that are hard to “measure”, such as the increased availability and use of HRCT after the 1990s, as well as the advances in CT scanner technology, have made it possible to define, with increasing precision, the imaging aspects associated with IPF. , Additionally, the practice of multidisciplinary discussions among medical specialties involved in the investigation of interstitial lung diseases has been adopted, which has facilitated the differential diagnosis of this group of diseases, possibly preventing misclassifications.

Demographic factors have also contributed to explaining the increase in mortality rates. According to census counts, between 1980 and 2010, the Brazilian population over 50 years of age increased by 68%. The standardized mortality rates per person-years ranged from 0.69/100,000 population in the first period (1979-1984) to 4.60/100,000 population in the last period (2011-2014); therefore, there was an almost sevenfold increase between the first and last observation periods. In comparison, in England and Wales, an analysis of mortality from IPF-S as the underlying cause revealed a sixfold increase between the first and last observation periods (total observation period: 1968 to 2008). Although the period analyzed in that study does not exactly correspond to that of the present study, the greater increase observed here could be partially explained by the marked change in Brazil’s population pyramid relative to that of those countries.

The standardized mortality rates per person-years were higher in men in all age groups. However, the male-to-female death ratio was 0.97, in comparison with 1.57 in the study of mortality in England and Wales and with 1.29 in the study of the incidence of IPF in Canada. A likely explanation for this finding is the marked difference in life expectancy between men and women in Brazil compared with industrialized countries: over the period 2005-2010, the difference in life expectancy between men and women was close to 4 years in Great Britain and in Canada, whereas it was 8 years in Brazil. The numerical excess of deaths from IPF-S in women is essentially due to greater female survival and to the greater number of females in the general population. Except in the 85-year or older age group, all RRs were significantly higher in men. It should be noted that, in three referral centers in Brazil, the extent to which the proportion of individuals with IPF was higher among men than women was more pronounced than that found in the present study.

Problems in comparing statistics on IPF lie in the data extraction source and in the criteria for defining the disease. - However, there is a consensus that the incidence of IPF has been gradually rising, approaching that of diseases such as stomach and liver cancer, as well as cervical cancer. , The nationwide study in Canada revealed a higher incidence of IPF in the more industrialized areas of the country. In the literature, there are reports of significantly increased risk of IPF for the following occupational and environmental factors: working in agriculture; raising livestock; working with wood; being exposed to metal dusts; being exposed to dusts containing silica; and smoking. This indicates the need to explore possible occupational and environmental exposures. Recently, it has been demonstrated that deaths from IPF-S and from mesothelioma show a significant linear relationship with historical data on asbestos imports in Great Britain. It should be pointed out that, because of its clinical presentation and, especially, its imaging presentation, asbestosis is one of the main differential diagnoses of IPF. A review of 1,718 cases of lobectomy for lung tumors demonstrated that, among patients with asbestos exposure markers, the proportion of histology results consistent with usual interstitial pneumonia was significantly higher than it was among patients without these markers, it being concluded that there is evidence that asbestos exposure leads to the development of IPF. Given the long-term, historically high consumption of asbestos in Brazil, this merits attention, because the legal and clinical management implications of asbestosis are very different from those of IPF.

Population-based mortality data have the advantages of the strength of large numbers, of not requiring extrapolations, and of making it possible to analyze long historical series. Because of the high lethality of IPF, mortality and incidence rates are close. However, an important limitation associated with the use of mortality data lies in the reliability of the diagnosis. In the present analysis, it was not possible to review the clinical data on the recorded deaths. It is likely that other interstitial lung diseases have been erroneously coded as ICD-9 515, ICD-9 516.3, or as ICD-10 J84.1, which is why we adopted the term IPF-S, similarly to the study by Navaratnam et al. Differences in utilization of ICD codes by physicians who certify deaths and the lack of a specific ICD code for IPF preclude the determination of the actual mortality rate. Another important limitation concerns the completeness of mortality data in the SIM. From 2000 onward, the estimated coverage of the SIM was 100% for the southern region and nearly 100% for the southeastern region. However, for the northern region, in 2000 and 2010, the estimated coverage was 75.3% and 85.4%, respectively. Given the heterogeneity in the completeness of mortality data in Brazil, it is likely that there are regional differences in mortality from IPF-S.

Although IPF is highly lethal, underlying cause-of-death data alone could mask the actual incidence of the disease. In England and Wales, it was estimated that approximately 60% of patients with known IPF had this disease recorded as the underlying cause of death on the death certificate. Despite its limitations, the use of underlying cause-of-death data revealed an upward temporal trend similar to that obtained when analyzing IPF cases derived from outpatient clinical records of chronic disease patients, which are records that have higher diagnostic accuracy.

We conclude that mortality from IPF-S in Brazil shows an upward trend, similar to that reported in other countries, but with low mortality rates, suggesting that the disease is underdiagnosed and/or underreported. The change in Brazil’s population pyramid does not alone explain the observed increase. Other studies with appropriate designs should be conducted to provide a better understanding of the findings.