SciELO - Scientific Electronic Library Online

 
vol.61Adverse events following Quadrivalent HPV vaccination reported in Sao Paulo State, Brazil, in the first three years after introducing the vaccine for routine immunization (March 2014 to December 2016)Usefulness of 3’- 5’ IS6110-RFLP genotyping and spoligotyping of Mycobacterium tuberculosis isolated in a tertiary hospital: a retrospective study detecting unsuspected epidemiological events índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

Compartilhar


Revista do Instituto de Medicina Tropical de São Paulo

versão On-line ISSN 1678-9946

Rev. Inst. Med. trop. S. Paulo vol.61  São Paulo  2019  Epub 12-Set-2019

http://dx.doi.org/10.1590/s1678-9946201961047 

ORIGINAL ARTICLE

Triatomine dispersion rates and their association with socioeconomic and environmental conditions in Northeastern Brazil, from 2009 to 2013

Thiago Bernardo-Pedro1 
http://orcid.org/0000-0002-7768-0193

Danielle Misael de Sousa2 

Simone Patrícia Carneiro de Freitas2 

Assilon Lindoval Carneiro de Freitas3 

Jacenir Reis dos Santos-Mallet2 

Wagner de Souza Tassinari1  4 

1Universidade Federal Rural do Rio de Janeiro, Instituto de Veterinária, Departamento de Parasitologia Animal, Rio de Janeiro, Rio de Janeiro, Brazil

2Fundação Oswaldo Cruz, Instituto Oswaldo Cruz, Laboratório Interdisciplinar de Vigilância Entomológica em Diptera e Hemiptera, Rio de Janeiro, Rio de Janeiro, Brazil

3Secretaria da Saúde do Estado do Ceará, 20ª Coordenadoria Regional de Saúde, Crato, Ceará, Brazil

4Universidade Federal Rural do Rio de Janeiro, Instituto de Ciências Exatas, Departamento de Matemática, Seropédica, Rio de Janeiro, Brazil


ABSTRACT

Triatomines are vectors of Trypanosoma cruzi , the etiological agent of Chagas disease, which affects between 6 and 12 million people in Latin America, with an incidence rate of 12 thousand cases per year. In the Ceara State, the predominance of the caatinga biome, coupled with a large rural area with precarious human habitations, provides several s shelters for these insects. In this study, we determined the spatiotemporal distribution of triatomine dispersion rates in the Cariri region, Southern Ceara and ascertained the possible association between these rates with socioeconomic and environmental factors. Dispersion rates (number of positive localities/number of searched localities × 100) were analyzed regarding 13 municipalities from the Ceara State, from 2009 to 2013. Socioeconomic and environmental variables collected from national research institutes were associated with the dispersion rates and their local empirical Bayesian estimates. All the municipalities recorded dispersion rates over 10% in all years, and 11 municipalities had average rates over 40% for the period of study. Significant differences were observed among the municipality means. The highest rates were observed in Antonina do Norte and Potengi. According to the correlation analysis, the proportion between the occupied population and the total population showed a significant negative correlation, as well as the percentage of the population who lives under adequate sanitary conditions. Both, the percentage of revenues from external sources and the percentage of urban households in reforested blocks had a significant positive correlation. Our results show that socioeconomic and environmental variables can be factors that contribute to both, the maintenance and the reduction of the elevated dispersion rates observed in the study area. Similar researches that encompass more municipalities from that region may reinforce Chagas disease surveillance and control in the Northeast of Brazil.

Key words: Triatominae; Socioeconomic factors; Public health; Spatial analysis; Bayesian analysis

INTRODUCTION

Chagas disease, also known as American Trypanosomiasis, is caused by Trypanosoma cruzi , a protozoan first described by Carlos Chagas in the Brazilian state of Minas Gerais, in 1909. This disease is a zoonosis and is endemic in 21 countries in the Americas, where it affects approximately 6 million people. There are 70 million individuals exposed to the risk of infection in the Americas, and the incidence of Chagas disease is around 30 thousand cases per year. Although mortality has significantly declined, the disease can cause irreversible and chronic consequences on the heart, digestive system and nervous system1 .

The insects that transmit the Trypanosoma cruzi protozoan belong to the subfamily Triatominae (Hemiptera, Reduviidae), composed by exclusively hematophagous members, in all nymph stages as well as in adulthood2 .

In the Ceara State, Northeastern Brazil, the first human cases of Chagas disease were diagnosed by xenodiagnosis, in 19423 . In the 1980s, the estimated prevalence of human Chagas disease in Ceara was 0.84%, and there were seropositive people in 93 out of 141 municipalities4 . From 1990 to 1997, 35 individuals died of the disease in the state5 . The last national survey on the seroprevalence of Chagas disease that occurred from 2001 to 2008, found six children infected with T. cruzi in Ceara. In two of these cases, the conditions indicated a probable vector transmission6 . Estimates of the World Health Organization (WHO) based on data from 2010 indicated a total of 46 cases of Chagas disease resulting from vector transmission in Brazil, each year7 , although it is estimated that only 10-20% of Chagas disease cases are actually notified, corroborating the Brazilian history of underreporting of Chagas and other diseases of compulsory notification8 .

There are 152 species allocated to the subfamily Triatominae, grouped into 18 genera containing five tribes9 . The Ceara State presents land predominantly covered by the biome caatinga, with a large rural area and precarious human habitations. This ecosystem provides shelters to some triatomines of public health importance, such as Triatoma brasiliensis Neiva 1911, Triatoma pseudomaculata Corrêa & Espínola 1964, Panstrongylus lutzi (Neiva & Pinto 1923), Panstrongylus megistus Burmeister 1835 and Rhodnius nasutus Stål 185910 .

The aim of this study was to analyze the triatomine dispersal rates from the capture records of the control program in the municipalities of the Cariri region, Ceara State, from 2009 to 2013, and to verify whether these data are associated with socioeconomic and environmental conditions of the municipalities. The development of studies in the Northeast region is justified by the high indexes of triatomine infestation5 , 11 , 12 . In addition, secondary species for transmission of Chagas disease such as P. lutzi and R. nasutus seem to have started a peridomicile and intradomicile cycle5 , 13 so that the entomological surveillance is fundamental to manage actions for Chagas disease control.

MATERIALS AND METHODS

Study area

The Cariri region is situated in the Southern Ceara State, Northeast of Brazil, between latitudes 6º41’40’’N and 7º39’00’’S, and longitudes 39º37’20’’E and 40º30’00’’W, in the semi-arid region of the state. This region comprises 13 municipalities, belonging to the 20th Regional Health Coordination [ Coordenadoria Regional de Saúde (CRES)] in Ceara. They cover an area of 8,885,668km2 , with a population of 328,410 people, according to the Brazilian Institute of Geography and Statistics [ Instituto Brasileiro de Geografia e Estatística (IBGE)]14 . The climate is characterized by scarce and irregular rainfall, little cloud cover, strong solar radiation, high evaporation rates and mean temperature around 27 ºC. The vegetation is predominantly the caatinga. The municipalities are: Altaneira, Antonina do Norte, Araripe, Assare, Campos Sales, Crato (main city), Farias Brito, Nova Olinda, Potengi, Salitre, Santana do Cariri, Tarrafas and Varzea Alegre.

Entomological surveillance performed from 1998 to 2008 in the Cariri region showed high indices of triatomine infestation in rural localities. In that period, T. brasiliensis and T. pseudomaculata were captured in the domestic and peridomestic environments of all the municipalities and at all stages of development11 , thus the area stands out from the other areas of the state.

Data source

Vector control-surveillance staff from the municipalities collected triatomines in the peridomestic and domestic environments, and sent the specimens to the Regional Health Coordination (CRES). The CRES, in turn, elaborated a database concerning the dispersion rate ( I d = number of positive localities/number of researched localities × 100) for the five years of study (2009-2013) and for each municipality.

The socioeconomic and environmental variables were retrieved from national researches, at the municipal level, made by national research institutes [(Brazilian Institute of Geography and Statistics (IBGE), the National Institute for Educational Studies and Research (INEP), the National Treasury (STN) and the United Nations Development Programme (PNUD)]. The following variables were selected: average monthly salary of formal workers15 , percentage of occupied population15 , percentage of population with nominal monthly income of at least half the minimum salary14 , 6-14 schooling years rate in regular schools14 , IDEB score in the early years of the elementary school16 , IDEB score in the final years of the elementary school16 , PIB per capita17 , percentage of revenues from external sources18 , Municipal Human Development Index19 , percentage of population with adequate sanitary sewage at home14 , percentage of urban households located in reforested blocks14 , percentage of urban households located in blocks with paving, curb and sidewalk14 . These were secondary data, and approval from the research ethics committee for studies involving human participants was not required, according to the National Health Council, Resolution 510/201620 .

Statistical analysis

First, an exploratory data analysis was performed. Then, we prepared choropleth maps related to the altitude and the triatomine dispersion rates ( I d ) recorded in the municipalities.

For each year, the municipalities were divided into three groups: low, midrange or high dispersion rate. The group definition took into consideration all the values presented by the 13 municipalities throughout the period 2009-2013, which were divided into two quantiles. These quantiles were used to form the groups, as they analyzed the triatomine dispersion rates in that region, at the study moment.

As the geographical analysis units (in this case, the municipalities) present varied sizes, the estimation of crude rates is not appropriate, since the small number of observations in some units leads to unrepresentative estimates with high variability. Bayesian smoothing may be the solution for this problem, including the geographical localization of the municipality as an additional information, promoting better estimates. These spatial statistical methodologies are called bayesian empirical, considering that auxiliary data are obtained from the sample itself21 .

Due to the high variability of the dispersion rates among municipalities, we calculated the dispersion rates based on global and local Bayes estimators22 . The Shapiro-Wilk normality test was used to verify the normality of dispersion rates and their respective empirical Bayesian estimates, as well as the other variables. To compare the dispersion rate with its respective empirical Bayesian estimates, the Pearson correlation test was applied. The Analysis of Variance test (ANOVA) was applied to verify whether there was variation of the dispersion throughout the years and among the municipalities studied23 , 24 .

Afterwards we tested the association between socioeconomic and environmental variables and the dispersion rates adjusted by the Bayes estimator.

Choropleth maps of socioeconomic and environmental variables that showed significant Pearson’s linear correlation were plotted. For all statistical tests, the significance level was set at 5% (α).

The freely available R statistical package version 3.3.225 was used to run the statistical analysis. The Quantum GIS version 2.18 (Las Palmas)26 was used to elaborate the maps and the GeoDa software version 1.1227 was used to calculate the Bayesian estimates and run the spatial statistical analysis.

RESULTS

Among the 13 municipalities studied, Potengi was the only one to present dispersion rates over 50% in all the years of the research. The municipalities of Assare and Varzea Alegre showed rates above half in four of the five years studied. Salitre presented rates of at least 50% for three years. Santana do Cariri and Tarrafas showed rates above half in two years. Municipalities that have not shown dispersion rates over 50% in any of the studied years were Crato and Nova Olinda. The municipality with the highest mean dispersion rate for the five years was Antonina do Norte ( I d = 92.31), followed by Potengi ( I d = 78.38), while Nova Olinda ( I d = 20.62) showed the lowest mean, followed by Crato ( I d = 28.39) ( Figure 1 ).

Figure 1 Triatomine dispersion rates for the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State, from 2009 to 2013. The boxes represent the variability of dispersion rates for each municipality, and the points represent outlier rates. 

There was a slight variation of the mean dispersion rate throughout the studied years ( Figure 2 ). The municipalities of Antonina do Norte, Campos Sales, Crato and Tarrafas did not send sent complete collection data to the CRES.

Figure 2 Violin plots of triatomine dispersion rates for the five studied years in the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State. 

According to the chart presented in Figure 2 and the result of the Analysis of Variance (ANOVA), there is not a significant difference ( p-value = 0.92) among the dispersion rates, throughout the study. However, for the same period, there is heterogeneity of the dispersion rates among municipalities ( p-value < 0.01). According to the Tukeytest, we observed that the municipalities of Antonina do Norte and Potengi presented the highest dispersion rates for the period, they do not differ from each other, but they differ from almost all the other municipalities. The same occurs to Nova Olinda and Crato, as they showed the lowest rates and do not differ from each other, but they differ from all the other municipalities .

The municipalities of Potengi, Assare, Antonina do Norte and Varzea Alegre showed the highest dispersion rates in the region. The first three are bordered and are in the North middle portion of the study area. Potengi and Assare are located West of the Araripe Plateau and belong to the Chapada do Araripe micro-region. Both municipalities are situated at over 450 m altitude ( Figure 3 ). The municipalities of Antonina do Norte and Varzea Alegre belong to the Varzea Alegre micro-region and are at over 300 m altitude. Varzea Alegre is the second most populous municipality of the research area, with 38.434 inhabitants.

Figure 3 Map of the mean altitude of the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State. 

Araripe, Campos Sales and Salitre represent the highest municipalities from the study area, once they are between 570 and 680 m in altitude and showed midrange dispersion rates for the region. The municipalities of Crato and Nova Olinda showed the lowest rates, always under 40%. Both municipalities belong to the Cariri Metropolitan Region (CMR) and are in the most economically active portion of the study area, being Crato the most populous among the 13 municipalities studied, with 121,428 inhabitants.

The municipality of Potengi was the only one included in the high dispersion rate group, every year. Assare and Varzea Alegre showed high rates in four of the five years. The municipalities of Crato and Nova Olinda presented low dispersion ratesin all the studied years. In 2009, six municipalities formed the high dispersion rate group, against only three in the last year of the study ( Figure 4 ).

Figure 4 Maps of triatomine dispersion rates from 2009 to 2013 and the mean for the period, in the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State. 

As a high variability was found for the dispersion rates, we performed global and local empirical Bayesian estimates for all the studied years ( Figure 5 ). These estimates were compared to each other and also to the crude values, to check whether there were correlations. There was a positive significant (α = 5%) correlation between the global Bayesian estimates and the crude rates (ρ rho = 0.95 ), between the local Bayesian estimates and the crude rates (ρ rho = 0.92 ), and between global and local estimates (ρ rho = 0.97 ).

Figure 5 Triatomine dispersion rates and global and local Bayesian estimates for the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State, from 2009 to 2013; “gebayes” represents the global Bayesian estimates, “lebayes” represents the local Bayesian estimates, and “raw data” represents the original dispersion rates. 

The Bayesian estimates promoted a subtle change regarding the themed maps. The municipality of Farias Brito changed from midrange to low dispersion rate. In addition, regarding the local empirical Bayes, the municipality of Varzea Alegre changed from high to midrange rate, what probably has occurred due to spatial effects of local Bayes ( Figure 6 ).

Figure 6 Maps of global and local Bayesian estimates of the triatomine dispersion rates for the average period 2009-2013, in the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State. 

The correlation analysis between the socioeconomic and environmental variables and the local Bayesian estimates of the dispersion rates showed a negative significant (α = 5%) correlation for a 6-14 schooling years rate in 2009 and 2013. On the other hand, the percentage of revenues from external sources showed a positive significant correlation for all the studied years, as well as for the mean of the period. Regarding territory and environment, the percentage of population with adequate sanitary sewage at home showed a negative significant correlation in 2010 and 2012, as well as for the mean. On the contrary, the percentage of urban households located in reforestation blocks showed a positive significant correlation in 2010, 2013 and also for the mean of the period ( Table 1 ).

Table 1 Pearson’s linear correlation coefficients and the respective 95% confidence intervals (CI) to represent the association between socioeconomic and environmental variables with local Bayesian estimates of triatomine dispersion rates, in the 13 municipalities forming the 20th Regional Health Coordination (CRES) of the Ceara State, from 2009 to 2013. 

Pearson’s Correlation [95% CI]

2009 2010 2011 2012 2013 Mean
Salary -0.27 [-0.71; 0.32] 0.19 [-0.46; 0.70] 0.12 [-0.51; 0.67] -0.05 [-0.66; 0.59] 0.03 [-0.61; 0.64] -0.09 [-0.61; 0.47]
Occupied population -0.61* [-0.87; -0.09] -0.63* [-0.89; -0.06] -0.60* [-0.88; -0.01] -0.25 [-0.76; 0.44] -0.35 [-0.80; 0.35] -0.59* [-0.86; -0.06]
Income of ½ salary 0.55* [0.05; 0.84] 0.21 [-0.44; 0.72] 0.31 [-0.34; 0.77] -0.28 [-0.77; 0.42] 0.16 [-0.52; 0.71] 0.30 [-0.29; 0.73]
Schooling -0.62* [-0.87; -0.12] -0.48 [-0.84; 0.15] -0.49 [-0.84; 0.14] -0.14 [-0.71; 0.53] -0.69* [-0.92; -0.11] -0.53 [-0.83; 0.02]
IDEB early years -0.21 [-0.68; 0.38] -0.49 [-0.84; 0.14] -0.41 [-0.81; 0.24] -0.89* [-0.97; -0.61] -0.63* [-0.90; -0.01] -0.48 [-0.81; 0.09]
IDEB final years -0.16 [-0.65; 0.42] -0.44 [-0.82; 0.21] -0.44 [-0.82; 0.20] -0.62 [-0.90; 0.08] -0.67* [-0.91; -0.08] -0.43 [-0.79; 0.14]
PIB -0.59* [-0.86; -0.06] -0.25 [-0.74; 0.40] -0.40 [-0.80; 0.25] -0.23 [-0.75; 0.46] -0.38 [-0.81; 0.32] -0.45 [-0.80; 0.13]
Revenues from external sources 0.57* [0.03; 0.85] 0.77* [0.33; 0.93] 0.79* [0.37; 0.94] 0.67* [0.07; 0.91] 0.68* [0.09; 0.91] 0.80* [0.44; 0.93]
MHDI -0.81* [-0.94; -0.47] -0.41 [-0.81; 0.24] -0.52 [-0.85; 0.10] -0.24 [-0.75; 0.45] -0.51 [-0.86; 0.16] -0.56* [-0.85; -0.02]
Sanitary sewage -0.49 [-0.82; 0.07] -0.74* [-0.93; -0.26] -0.71* [-0.92; -0.20] -0.64* [-0.90; -0.01] -0.58 [-0.88; 0.07] -0.57* [-0.85; -0.03]
Deforestation 0.54 [-0.01; 0.84] 0.67* [0.12; 0.90] 0.46 [-0.18; 0.83] 0.56 [-0.09; 0.88] 0.66* [0.05; 0.91] 0.64* [0.14; 0.88]
Urbanization -0.03 [-0.57; 0.52] -0.07 [-0.64; 0.55] -0.09 [-0.65; 0.53] -0.17 [-0.72; 0.51] -0.15 [-0.71; 0.52] -0.15 [-0.64; 0.43]

*significance level set at 5%; IDEB (Basic Education Development Index); PIB (Gross Domestic Product); MHDI (Municipal Human Development Index)

Potengi was the municipality with the lowest 6 to 14 schooling years rate (92.7%). On the contrary, Campos Sales and Farias Brito had the highest rates (98.3% and 99.3%, respectively). Regarding the revenues from external sources, Antonina do Norte leads with the highest percentage (97.7%), while the municipalities of Crato, Araripe and Nova Olinda showed the lowest values in the studied area (82.4%, 84.9%and 85.6%, respectively). The municipality with lowest sewage cover was Tarrafas (2.2%), followed by Potengi (7.8%). However, Potengi had the highest reforestation rate (96.3%), followed by Tarrafas (94.4%) ( Figure 7 ).

Figure 7 Maps of occupied population; IDEB score in the early years; revenues from external sources; MHDI; sanitary sewage and afforestation, according to IBGE, for the 13 municipalities that form the 20th Regional Health Coordination (CRES) of the Ceara State. 

DISCUSSION

In a study conducted between 1998 and 2008, regarding the same municipalities of the present research, Gonçalves et al. 11 have also found the highest dispersion rates in Antonina do Norte, Assare, Potengi and Varzea Alegre. These authors believed that these municipalities presented these results due to the presence of specific factors such as environment conditions, type of vegetation, territorial proximity and the use of tree trunks to build shelters for domestic animals. All of these factors may influence the dispersion of triatomine bugs to areas that are not covered by the surveillance system, contributing to the maintenance of triatomines in the region13 .

The municipalities of Araripe, Campos Sales and Salitre are among the four with the highest altitudes in the study area, and they presented dispersion rates between 40% and 52% for the average period of 2009-2013. Most of the geographical area of these municipalities is on the Araripe Plateau, situated within the caatinga biome, among the states of Ceara, Pernambuco and Piaui. This plateau offers a wide variety of phytophysiognomies and an environmental dynamic that is distinct from the other lowland areas within the caatinga. The area is under heavy anthropogenic pressure due to agricultural expansion, disorderly occupation and hunting28 . Lima et al. 29 , when investigating dwellings and outbuildings in rural areas in the Sergipe State, calculated the triatomine dispersion rate in 66.7%, with triatomines found in both, domestic and peridomestic environments. In the latter, the authors frequently found corrals, waste deposits and chicken coops, and many of which consisted of makeshift structures. This dispersion rate matches the rates presented by most of the municipalities in our study.

The municipality of Farias Brito, whose dispersion rate was around 40% in the period 2009-2013, presents a deforestation level over 80%, according to IBGE14 . Freitas et al. 10 reported, during the research in the municipality, that the species R. nasutus was found associated with bird nests in palm trees, as macauba ( Acrocomia aculeata ), a quite common palm tree in the region, have become the species natural habitat. Still according to these authors, the municipality of Farias Brito is predominantly covered by hyper xerophilic caatinga, where shrubs like marmeleiro ( Croton sp.) and jurema-preta ( Mimosa tenuiflora ) predominate. These species of plants are used by the rural population as a domestic energy source since their wood is used to produce fire to cook and to build shelters for animals10 . The increased deforestation and anthropogenic invasion of natural ecotopes of Chagas disease can lead to the urbanization of insects that find shelters in deforested areas of municipalities and rural peridomiciles.

In a study conducted by Dias et al. 30 in the Chapada do Araripe, 521 specimens of R. nasutus were collected in five palm species: babassu ( Attalea speciosa ), buriti ( Mauritia flexuosa ), carnauba ( Copernicia prunifera ), catole ( Syagrus oleracea ) and macauba-barriguda ( Acrocomia intumescens ). From the 50 palms that were cut down and dissected, 86% harbored at least one triatomine specimen or presented with exuvia or eggs. Four hundred and fifteen triatomines (79.7%) were found in babassu and catole specimens.

When analyzing natural ecotopes of T. brasiliensis populations in the Jaguaribe Valley, in Ceara, Valença-Barbosa et al. 31 collected 157 triatomines in 21 specimens of xiquexique ( Pilosocereus gounellei ), a shrubby cactus. In total, 44 xiquexique shrubs were investigated corresponding to an infection rate of 47.7%. The authors suggested that P. gounellei is an important ecotope of wild T. brasiliensis populations in the caatinga, and that bugs from cacti might be involved in house infestation, in lowland areas where no rocky outcrops occur.

Jurema-preta stands out among the plant species found in the Cariri region, due to its high density and high geographical coverage. Boles and boughs of M. tenuiflora are widely used by rural populations, either for protection and isolation of domestic animals or as energy source. However, triatomines in the young form can be found under the bark of the plant. In addition, its storage in the peridomestic environment attracts synanthropic animals, as possums and rodents, easing the triatomine dispersion from the natural to the peridomestic environment13 . The common use of endemic palm trees and shrubs as energy sources, in the construction of shelters for animals, or yet as an income source through the commercialization of their byproducts, places the rural populations from the Cariri region in a vulnerable position regarding the vector transmission of Chagas disease. The introduction and storage of the above-mentioned plants in the peridomestic and/or domestic environments potentiates the risk of infection by T. cruzi in these communities, also contributing to the domiciliation process of species that are still considered as wild.

In a study conducted in the Pernambuco State, bordering Ceara, Silva et al. 32 reported the occurrence of 3,323 triatomines in both, the peridomestic and domestic environments of 185 municipalities. Six species were identified, and T. brasiliensis and T. pseudomaculata represented 73.4% of the specimens collected.

In the municipalities of Nova Olinda and Crato, the mean dispersion rate between 2009-2013 was under 30%, representing the lowest mean for the 13 municipalities during the study period. Both municipalities belong to CMR, where the urbanization level and population growth are far more marked than the other municipalities of the 20th CRES. According to IBGE33 , the participation of agriculture in the CMR economic structure is of only 2.84%. Besides that, the population has grown 14% from 2000 to 2010. The marked urbanization process in that metropolitan region has deviated the population from typically rural habits to the cities, decreasing the risk of vector transmission of Chagas disease.

The sampling technique used in the study is not the most appropriate considering that, according to Valença-Barbosa et al. 12 , the active searches made by trained and motivated agents have commonly detected around 40% to 60% of the infestation foci, while non-detected foci are not eliminated. Abad-Franch et al. 34 estimated the sensibility of the active research conducted by control agents by about 20% for localities with low intensity of infestation and around 40% for high intensity localities. Although the vector control-surveillance staff of the municipalities are trained by the Regional Health Coordination, the dispersion rates presented in the study may be underestimated due to the low sensibility of the sampling technique.

From 2007 to 2011, the SVS35 registered the capture of almost 800 thousand triatomines in the peridomestic and domestic environments. It is considered that this record encompasses less than 10% of all triatomines inside or around dwellings. According to the SVS36 , the entomological surveillance for Chagas disease must be implemented all over the country, primarily supported by two pillars: passive surveillance, when the population participates in the triatomine notification process; and active surveillance, made by municipal control teams in partnership with regional health coordination offices, without the need of prior notification.

From this context, the use of spatial statistical analysis techniques to study the occurrence of vectors of Chagas disease has been widely disclosed, both for the recognition of the importance of these methods and for accelerating softwares development12 , 34 .

Thereby, identifying more vulnerable areas for the occurrence of synanthropic triatomines has been a valuable tool to reorient actions related to prevention, control and epidemiological surveillance of Chagas disease37 . According to our results, the municipalities of Antonina do Norte, Assare, Potengi and Varzea Alegre showed high dispersion rates associated with high socioeconomic vulnerability. These municipalities must receive special attention from the Regional Health Coordination.

To fight the coexistence of areas with various degrees of risk, vector control and surveillance actions should be adjusted to the established risk. Aside from taking into consideration the operational capacity of the municipalities, the stratification of a traditionally endemic area is based on a set of variables that could potentially influence the infestation (or reinfestation) process and, consequently, the vector transmission of Chagas disease in the domestic environment35 .

Prevention and control interventions targeted at elevated risk areas represent effective strategies to limit possible outbreaks, not just for reducing the operational costs but also for supporting the control of neglected tropical diseases38 .

CONCLUSION

The results of this study suggest that socioeconomic and environmental variables may partially explain the high triatomine dispersion rates observed in some municipalities from the Cariri region. Children schooling rate, as well as the municipality dependency on external resources, proved important factors related to the maintenance of high triatomine dispersion rates. The deforestation level of the municipalities also seemed to play an important role in triatomine dispersion, since many triatomine species are commonly found associated twith palm trees and some shrub species.

The municipalities belonging to the Cariri Metropolitan Mesoregion showed the lowest dispersion rates, suggesting that the urbanization process has apparently reduced the risk of the vector transmission of Chagas disease to that population, despite the social problems generated by this type of process. Meanwhile, the use of some plant species in the rural areas of Ceara State as domestic energy sources, for protection and isolation of domestic animals, and for commercialization of their byproducts, allows the maintenance of triatomines in the peridomestic and domestic environments.

Finally, we emphasize that studies addressing the relationship between triatomine dispersion, socioeconomic conditions and urbanization processes are the key for an efficient control of the transmission of Chagas disease in the Brazilian Northeast.

ACKNOWLEDGMENTS

To the Coordenadoria Regional de Saúde , Crato, CE (CRES) for its technical support during the field collections, to Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for grants and to IOC/FIOCRUZ for its scientific support.

REFERENCES

1. Pan American Health Organization. Chagas disease: fact sheets. [cited 2019 Jul 29]. Available from: http://www.paho.org/hq/index.php ?option=com_topics&view=rdmore&cid=5965iem=chagasdisease&cat=communication&type=fact-sheets5965&Itemid=40743&lang=enLinks ]

2. Lent H, Wygodzisnky P. Revision of the Triatominae (Hemiptera, Reduviidae) and their significance as vectors of Chagas disease. Bull Am Mus Nat Hist. 1979; 163:123-520. [ Links ]

3. Alencar JE. História natural da doença de Chagas no estado do Ceará. Fortaleza: Universidade Federal do Ceará; 1987. [ Links ]

4. Camargo ME, Silva GR, Castilho EA, Silveira AC. Inquérito sorológico da prevalência de infecção chagásica no Brasil, 1975/1980. Rev Inst Med Trop São Paulo. 1984;26:192-204. [ Links ]

5. Dias JC, Machado EM, Fernandes AL, Vinhaes MC. Esboço geral e perspectivas da doença de Chagas no Nordeste do Brasil. Cad Saude Publica. 2000;16 Suppl 2:13-34. [ Links ]

6. Ostermayer AL, Passos AD, Silveira AC, Ferreira AW, Macedo V, Prata AR. O inquérito nacional de soroprevalência de avaliação do controle da doença de Chagas no Brasil (2001-2008). Rev Soc Bras Med Trop. 2011;44 Suppl 2:108-21. [ Links ]

7. World Health Organization. Chagas disease in Latin America: an epidemiological update based on 2010 estimates. Wkly Epidemiol Rec. 2015;90:33-43. [ Links ]

8. Dias JC, Ramos Jr AN, Gontijo ED, Luquetti A, Shikanai-Yasuda MA, Coura JR, et al. II Consenso Brasileiro em Doença de Chagas, 2015. Epidemiol Serv Saude. 2016; 25 Num Esp:7-86. [ Links ]

9. Jurberg J, Rodrigues JM, Moreira FF, Dale C, Cordeiro IR, Lamas Júnior VD, et al. Atlas iconográfico dos triatomíneos do Brasil (vetores da doença de Chagas). Rio de Janeiro: Instituto Oswaldo Cruz; 2014. [ Links ]

10. Freitas AL, Freitas SP, Gonçalves TC, Lima Neto AS. Vigilância entomológica dos vetores da doença de Chagas no município de Farias Brito, estado do Ceará – Brasil. Cad Saude Colet. 2007;15:231-40. [ Links ]

11. Gonçalves TC, Freitas AL, Freitas SP. Surveillance of Chagas disease vectors in municipalities of the state of Ceará, Brazil. Mem Inst Oswaldo Cruz. 2009; 104:1159-64. [ Links ]

12. Valença-Barbosa C, Lima MM, Sarquis O, Bezerra CM, Abad-Franch F. Modeling disease vector occurrence when detection is imperfect II: drivers of site-occupancy by synanthropic Triatoma brasiliensis in the Brazilian northeast. PLoS Negl Trop Dis. 2014;8:e2861. [ Links ]

13. Freitas SP, Freitas AL, Prazeres SM, Gonçalves TC. Influência de hábitos antrópicos na dispersão de Triatoma pseudomaculata Corrêa & Espínola, 1964, através de Mimosa tenuiflora (Willdenow) (Mimosaceae) no estado do Ceará, Brasil. Cad Saude Publica. 2004;20:333-6. [ Links ]

14. Instituto Brasileiro de Geografia e Estatística. Censo 2010: resultados. [cited 2019 Jul 29]. Available from: https://censo2010.ibge.gov.br/resultados.htmlLinks ]

15. Instituto Brasileiro de Geografia e Estatística. Cadastro Central de Empresas - CEMPRE. [cited 2019 Jul 29]. Available from: https://sidra.ibge.gov.br/pesquisa/cempre/tabelas/brasil/2015Links ]

16. Brasil. Ministério da Educação. Instituto Nacional de Estudos e Pesquisas Educacionais Anísio Teixeira. IDEB: resultados e metas. [cited 2019 Jul 29]. Available from: http://ideb.inep.gov.br/Links ]

17. Instituto Brasileiro de Geografia e Estatística. Diretoria de Pesquisas. Coordenação de Contas Nacionais. Produto Interno Bruto dos municípios 2010. Rio de Janeiro: IBGE; 2012. [cited 2019 Jul 29]. Available from: https://biblioteca.ibge.gov.br/visualizacao/livros/liv62930.pdfLinks ]

18. Brasil. Ministério da Fazenda. Secretaria do Tesouro Nacional. Balanço do Setor Público Nacional (BSPN). [cited 2019 Jul 29]. Available from: http://www.tesouro.fazenda.gov.br/pt_PT/balanco-do-setor-publico-nacional-bspn-Links ]

19. Programa das Nações Unidas para o Desenvolvimento. Atlas do desenvolvimento humano no Brasil. [cited 2019 Jul 29]. Available from: http://www.atlasbrasil.org.brLinks ]

20. Brasil. Ministério da Saúde. Conselho Nacional de Saúde. Resolução nº 510, de 7 de abril de 2016. Diário Oficial da União, Brasília, 24 maio 2016. Seção 1:44-6. [cited 2019 Jul 29]. Available from: http://conselho.saude.gov.br/resolucoes/2016/Reso510.pdfLinks ]

21. Marshall RJ. Mapping disease and mortality rates using empirical bayes estimators. J R Stat Soc Ser C Appl Stat. 1991;40:283-94. [ Links ]

22. Druck S, Carvalho MS, Câmara G, Monteiro AV, editores. Análise espacial de dados geográficos. Brasília: EMBRAPA; 2004. [ Links ]

23. Vieira S. Análise de variância: Anova. São Paulo: Atlas; 2006. [ Links ]

24. Bussab WO, Morettin PA. Estatística básica. 9aed. São Paulo: Saraiva; 2017. [ Links ]

25. The R Foundation. The R Project for Statistical Computing. [cited 2019 Jul 29]. Available from: https://www.R-project.org/Links ]

26. Athan T, Dassau O, Ghisla A, Homann M, Macho W, Engel CA, et al. Quantum GIS user guide: version 1.7 “Wroclaw”. [cited 2019 Aug 8]. Available from: http://download.osgeo.org/qgis/doc/manual/qgis-1.7.0_user_guide_en.pdfLinks ]

27. Anselin L, Syabri I, Kho Y. GeoDa: an introduction to spatial data analysis. Geogr Anal. 2006;38:5-22. [ Links ]

28. Novaes RL, Felix S, Souza RF. Save Caatinga from drought disaster. Nature. 2013;498:170. [ Links ]

29. Lima AF, Jeraldo VL, Silveira MS, Madi RR, Santana TB, Melo CM. Triatomines in dwellings and outbuildings in an endemic area of Chagas disease in northeastern Brazil. Rev Soc Bras Med Trop. 2012;45:701-6. [ Links ]

30. Dias FB, Bezerra CM, Machado EM, Casanova C, Diotaiuti L. Ecological aspects of Rhodnius nasutus Stål, 1859 (Hemiptera: Reduviidae: Triatominae) in palms of the Chapada do Araripe in Ceará, Brazil. Mem Inst Oswaldo Cruz. 2008;103:824-30. [ Links ]

31. Valença-Barbosa C, Lima MM, Sarquis O, Bezerra CM, Abad-Franch F. A common Caatinga cactus, Pilosocereus gounellei, is an important ecotope of wild Triatoma brasiliensis populations in the Jaguaribe Valley of Northeastern Brazil. Am J Trop Med Hyg. 2014;90:1059-62. [ Links ]

32. Silva MB, Menezes KR, Siqueira AM, Balbino VQ, Lorosa ES, Farias MC, et al. Importância da distribuição geográfica dos vetores da doença de Chagas em Pernambuco, Brasil, em 2012. Rev Patol Trop. 2015;44:195-206. [ Links ]

33. Instituto Brasileiro de Geografia e Estatística. Sinopse do Censo Demográfico 2010: Ceará. [cited 2019 Jul 29]. Available from: https://censo2010.ibge.gov.br/sinopse/index.php?dados=29&uf=23Links ]

34. Abad-Franch F, Valença-Barbosa C, Sarquis O, Lima MM. All that glisters is not gold: sampling-process uncertainty in disease-vector surveys with false-negative and false-positive detections. PLoS Negl Trop Dis. 2014;8:e3187. [ Links ]

35. Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Doença de Chagas aguda no Brasil: série histórica de 2000 a 2013. Bol Epidemiol. 2015;46(21):1-9. [ Links ]

36. Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Guia de vigilância em saúde. Brasília: Ministério da Saúde; 2014. [ Links ]

37. Vinhaes MC, Oliveira SV, Reis PO, Lacerda Sousa AC, Silva RA, Obara MT, et al. Assessing the vulnerability of Brazilian municipalities to the vectorial transmission of Trypanosoma cruzi using multi-criteria decision analysis. Acta Trop. 2014;137:105-10. [ Links ]

38. Martins-Melo FR, Ramos Jr AN, Alencar CH, Heukelbach J. Trends and spatial patterns of mortality related to neglected tropical diseases in Brazil. Parasite Epidemiol Control. 2016;1:56-65. [ Links ]

Recebido: 20 de Maio de 2019; Aceito: 26 de Julho de 2019

Correspondence to: Thiago Bernardo-Pedro Universidade Federal Rural do Rio de Janeiro, Instituto de Veterinária, Departamento de Parasitologia Animal, BR-465, Km 7, CEP 23890-000, Seropédica, RJ, Brazil E-mail:thiago.b.pedro@gmail.com

AUTHORS’ CONTRIBUTION

WST and TBP designed the study and wrote the protocol; DMS and ALCF collected data; TBP and WST conducted the analysis; TBP drafted the manuscript; WST and DMS contributed to the manuscript; SPCF and JRSM were responsible for reviewing the manuscript.

Creative Commons License  This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.