Skip to main content

Effects of anthropogenic landscape changes on the abundance and acrodendrophily of Anopheles (Kerteszia) cruzii, the main vector of malaria parasites in the Atlantic Forest in Brazil

Abstract

Background

The mosquito Anopheles (Kerteszia) cruzii is the main vector of human and simian malaria in the Atlantic Forest. This species is usually abundant in the forests where it occurs, preferring to live and feed on canopies, behaviour known as acrodendrophily. However, in several studies and locations this species has been observed in high density near the ground in the forest. In this study, it was hypothesized that factors associated with anthropogenic landscape changes may be responsible for the variation in abundance and acrodendrophily observed in An. cruzii.

Methods

The study was conducted in a conservation unit in the city of São Paulo, Brazil. Monthly entomological collections were performed from March 2015 to April 2017, and the resulting data were used with data from another study conducted in the same area between May 2009 and June 2010. Mosquitoes were collected from five sites using CDC and Shannon traps. Landscape composition and configuration metrics were measured, and generalized linear mixed-effect models were used to investigate the relationship between these metrics and variations in the abundance and acrodendrophily of An. cruzii.

Results

The model that showed the best fit for the relationship between landscape metrics and An. cruzii abundance suggests that an increase in the proportion of forest cover leads to an increase in the abundance of this mosquito, while the model that best explained variations in An. cruzii acrodendrophily suggests that an increase in total forest-edge length leads to greater activity by this species at ground level.

Conclusion

While the data indicate that changes in landscape due to human activities lead to a reduction in An. cruzii abundance, such changes may increase the contact rate between this species and humans living on the edges of forest fragments where An. cruzii is found. Future studies should, therefore, seek to elucidate the effect of these landscape changes on the dynamics of Plasmodium transmission in the Atlantic Forest, which according to some studies includes the participation of simian hosts.

Background

Autochthonous malaria in the Atlantic Forest occurs mainly in the Southern and Southeastern regions of Brazil, where it is also called “bromeliad malaria” because the immature forms of the main mosquito vectors develop in the water that collects in bromeliads (Bromeliaceae), an abundant group of plants in this biome [1,2,3,4]. In recent decades there has been a low incidence of outbreaks of autochthonous malaria in the Atlantic Forest, and most cases have been asymptomatic with low circulating parasite loads [5, 6]. The main parasite involved in the transmission of malaria in the Atlantic Forest is Plasmodium vivax and, less frequently, Plasmodium malariae and Plasmodium falciparum [4].

The species known to be vectors of bromeliad malaria belong to the subgenus Kerteszia of the genus Anopheles, and are distributed along the Atlantic coast of Brazil although some can also be found in the Amazon region [7]. Anopheles (Kerteszia) cruzii is considered the species most associated with the transmission of human and simian Plasmodium in the Atlantic Forest [8,9,10,11,12]. Originally it had a widespread distribution in Brazil, extending from the Southern to Northeastern regions of the country [7]. However, with increasing deforestation in the Atlantic Forest its distribution has decreased [2, 4, 13].

Immature forms of An. cruzii usually develop in epiphytic or terrestrial bromeliads in shady places under the forest canopy and are very abundant in humid forests on coastal slopes [14, 15]. The species tends to be more abundant in hot, rainy periods but can be observed frequently in the forest even in drier, cooler seasons [16,17,18]. Anopheles cruzii females exhibit increased biting activity in the first few hours after evening twilight but can be found actively biting throughout the night or even in daylight hours in the forest, in peridomestic areas or even inside houses [16].

The term “acrodendrophily” refers to the tendency of certain wild species to live and preferentially feed in the canopy. The term was proposed by Garnham et al. [19], when they observed that some mosquitoes were collected more frequently in the upper stratum of the Kaimosi Forest in Kenya. Many studies have demonstrated that An. cruzii has acrodendrophilic behaviour although in several studies and locations this species has been observed in high density at ground level in the forest [20,21,22,23,24]. In locations where An. cruzii is more acrodendrophilic, human cases of malaria do not seem to occur or are rarer, even in the presence of monkeys infected with simian Plasmodium species. During surveys in Southern and Southeastern Brazil, Deane et al. [25, 26] noted that in the Serra da Cantareira Forest, in the state of São Paulo, more than 99% of An. cruzii individuals were collected in the canopy. Although 62% of the howler monkeys (Alouatta guariba clamitans) tested in the study area were infected by Plasmodium, only one natural, accidental human infection due to Plasmodium simium was observed. The same occurred in the municipality of Campo Alegre, in the State of Santa Catarina, where approximately 90% of the An. cruzii collected were collected in the canopies and 43% of the howler monkeys tested were infected with Plasmodium spp. but no human cases were observed [26]. A similar pattern was observed in studies in the municipalities of Guaíba, in the state of Rio Grande do Sul, and Santa Leopoldina, in the state of Espírito Santo [27, 28]. However, near the municipality of Joinville, Santa Catarina, where there have been cases of simian and human malaria, 42% of the An. cruzii collected were collected at ground level [26].

Deane et al. [21] tested the hypothesis of reproductive isolation among populations of An. cruzii collected in traps placed in the canopy and at ground level in a forest near the municipality of São Francisco do Sul, Santa Catarina. Using mark-recapture techniques with fluorescent dust, the authors noticed that there was a vertical movement of An. cruzii between the canopy and ground level in the forest as mosquitoes previously collected and marked in the canopy were later collected in ground traps and vice versa. This observation reinforced the hypothesis that under certain circumstances An. cruzii may carry Plasmodium species between humans and non-human primates, establishing a cycle of zoonotic malaria transmission in these regions. Recently, it has been shown that human infection by P. simium or other variant forms of P. vivax that circulate among monkeys in the Atlantic Forest may be more common than previously thought [29], reinforcing the need for further investigations of ecological aspects of zoonotic transmission of bromeliad malaria, especially in areas where An. cruzii occurs in high abundance.

Variations in the acrodendrophilic behaviour of An. cruzii could be due to the increased human presence in the Atlantic Forest and consequent changes in the environment. Although these changes have led to a degradation in the habitats that favour the development of this species, they can lead to increased contact between humans and females of the local populations of this vector as the increased human presence represents a greater supply of blood meals at ground level in the forest and, potentially, a reduced supply of other vertebrate hosts in the canopy. Landscape changes due to increased human presence can, therefore, be expected to lead to a reduction in An. cruzii abundance, but to favour an increase in the relative proportion of this mosquito searching for blood sources at ground level. To test this hypothesis, entomological surveys were performed and landscape metrics were measured in a remnant of Atlantic Forest on the outskirts of the megacity of São Paulo where several cases of human malaria have been reported in the last decades [11, 30, 31].

Methods

Study area

The Capivari-Monos EPA (Environmental Protection Area) is approximately 40 km from the center of the city of São Paulo, Brazil, and extends over the first hills near the crest of the Serra do Mar mountain range at altitudes varying from 740 to 850 m above sea level. The climate is a super-humid, oceanic, tropical climate with average annual temperatures of around 19 °C and rainfall of between 1600 and 2200 mm [32]. The original vegetation is dense montane ombrophilous forest found in remnants in varying degrees of conservation, from areas in regeneration dating from the 1950s, when logging was stopped, to areas recently degraded as a result of urban and rural expansion. The district of Engenheiro Marsilac lies within the EPA and has a population of around 10 thousand inhabitants and a population density of approximately 41 inhabitants per km2. The population, most of whom have low incomes, live mainly in rural settlements [32]. Cases of autochthonous malaria have been recorded in the district in the last decades as well as in neighbouring municipalities such as Embu-Guaçu, São Bernardo do Campo, Itanhaém and Juquitiba and the Serra do Mar State Park [11, 30, 31].

To select the sampling sites, landscape variations in the region were considered and were included areas with different degrees of anthropogenic influence (Fig. 1). The areas where specimens were collected were (1) Embura—a village surrounded by small farms and the EPA forest (23° 53.036′ S/46° 44.485′ W); (2) Marsilac—a village surrounded by the EPA forest and near a railway line (23° 54.395′ S/46° 42.486′ W); (3) transition zone—private property near Marsilac village constituting a transitional area between a rural environment and the EPA forest (23° 54.556′ S/46° 42.167′ W); (4) Cachoeira do Marsilac—private property in the EPA forest next to a waterfall with a visitation area (23° 56.378′ S/46° 41.659′ W); and (5) Evangelista de Souza—well-preserved EPA forest near a railway station (23° 56.140′ S/46° 38.090′ W).

Fig. 1
figure 1

Study sites in the Capivari-Monos EPA, São Paulo, Brazil: (1) Embura village, (2) Marsilac village, (3) Transition zone, (4) Cachoeira do Marsilac and (5) Evangelista de Souza. The areas were classified according to the map of Atlantic Forest biome remnants in the municipality of São Paulo (available at http://geosampa.prefeitura.sp.gov.br/PaginasPublicas/_SBC.aspx). Green represents dense ombrophilous forest; blue represents heterogeneous forest; pink, natural fields; and white, areas where there is human activity (roads, rural properties or villages). Crosses inside the circles indicate collection points. The inner circle represents a 500 m buffer and the dashed circle a 1 km buffer around the collection points. The map was created using QGIS v2.18.9 (http://www.qgis.org)

Field collections

Culicidae collections were carried out monthly from March 2015 to April 2017. To determine the frequency of An. cruzii in the canopy and at ground level, automatic CDC traps with CO2 bait were used at each collection point, one trap in the tree canopy in heights between 10 and 13 m, depending on the selected branch of the tree (measured according to the length of the ropes used to hoist the traps), and one trap about 1 m from the ground. Traps were installed early in the afternoon and removed the following morning after approximately 18 h of exposure (Fig. 2a, b). Considering that between the years of 2013 to 2015 in the region of study the sunset ranged from 5:28 p.m. (May–June) to 6:59 p.m. (January) and the sunrise ranged from 5:12 a.m. (November–December) to 6:50 a.m. (June–July) [33], the data obtained by CDC traps were not corrected according to seasonality changes in day-length [34], since the 18-h exposure period of the CDC traps ensured that all traps were operating during the hours of higher biting activity of An. cruzii [16, 17].

Fig. 2
figure 2

Collection techniques used in the study: a CDC trap installed near the ground; b CDC trap installed in the canopy of a tree; and c Shannon trap. Images: Laboratory for Entomology in Public Health, University of São Paulo (LESP/FSP/USP)

Shannon traps were set up near the CDC traps to collect mosquitoes during the first 2 h after evening twilight (when many mosquito species, including An. cruzii, exhibit increased biting activity) [16, 17]. In each Shannon trap, two individuals with manual battery-powered aspirators collected mosquitoes on the outer and inner surface of the tent (Fig. 2c).

For the Evangelista de Souza site, data on An. cruzii abundance and frequency at ground level and in the canopy were obtained from Ribeiro et al. [35] and Duarte et al. [11], who performed field collections from May 2009 to June 2010. These authors used the same collection methods and traps as those used in the present study and very kindly provided all the information on their collections in exhaustive detail to the present authors. For Embura village, the information from the studies cited above and the collections in the present study were included.

Specimens were identified morphologically in the Laboratory for Entomology in Public Health, School of Public Health, University of São Paulo (LESP/FSP/USP), with taxonomic keys for Culicidae [2, 15, 36].

Landscape analysis

To investigate whether landscape changes influence variations in An. cruzii abundance and acrodendrophilic behaviour, landscape composition and configuration metrics were calculated and used as explanatory variables. The locations of the Shannon traps were georeferenced in the field and plotted with QGIS 2.18 (http://www.qgis.org) on the map of Atlantic Forest biome remnants in the municipality of São Paulo (a 1:5000-scale orthophoto mosaic) available at http://geosampa.prefeitura.sp.gov.br/PaginasPublicas/_SBC.aspx. In order to investigate whether variations in the abundance and acrodendrophily of An. cruzii may be better predicted by a more local or wider landscape scale, buffers extending 500 m and 1 km were created around each collection point to define the surrounding landscape. These scales were chosen based on observed flight radius of An. cruzii, which accordingly to Ferreira et al. [37] was about 1000 m. Five classes of vegetation or land use were observed: dense ombrophilous forest, heterogeneous forest, natural fields, anthropogenic areas and villages. The classes “dense ombrophilous forest” and “heterogeneous forest”, environments that favour the development of An. cruzii, were grouped into a single class called “forest cover”. Similarly, the classes “anthropogenic areas” and “villages” were grouped into a class called “anthropogenic use”. Landscape composition was measured by considering the relative abundance of each landscape class within each buffer (except for the class “natural fields”, which was considered of minor importance). Landscape configuration, considered here as the degree of fragmentation of the class “forest cover” in each landscape matrix, was measured as the number of “forest cover” fragments and their total edge length in kilometres (Table 2). Landscape configuration metrics were determined with Fragstats v4.2 [38].

Data analysis

The abundance and frequency of An. cruzii collected at ground level were considered response variables. Mosquitoes collected in Shannon traps were used to calculate abundance because this technique allows the number of mosquitoes per human host/hour during peak blood-feeding periods (the first 2 h after evening twilight) to be measured. Therefore, the response variable “abundance” used in the statistical models refers to the average number of An. cruzii specimens per human/hour in the Shannon trap calculated for each sample and site. To investigate variations in An. cruzii acrodendrophily, the data from the CDC traps were used. Mosquito frequency at ground level was calculated by dividing the number of An. cruzii in the trap at ground level by the total number of individuals of this species collected at the site in the same sample (number of specimens collected in the CDC trap at ground level plus the number of specimens collected in the CDC traps in the canopy).

Since the samples collected at each site represent pseudo-replicates (repeated measures over time at the same site), generalized linear mixed-effect models were used [39]. Thus, the fixed effect was represented by the landscape composition and configuration metrics, and the random effect by the different years and months when collections were carried out. The variables year and month were considered as different random factors, since monthly variations reflect the effect of seasonality, while annual variations are less predictable, and may reflect atypical climatic or environmental conditions [40]. Because of the low number of sampling sites, it was decided to test models with only one predictive variable in the fixed effect. For the abundance models it was opted for Poisson errors (log link), and for the ground/canopy frequency binomial errors (logit link) was used. An information-theoretical approach based on the Akaike Information Criterion corrected for small samples (AICc) was applied to select the most plausible statistical models [40]. The models with the smallest AICc were considered the best fit, and ∆AICc ≤ 2 was adopted as the cutoff to select models with more empirical support. The strength of evidence in favour of each model was evaluated using Akaike weights [41]. The selected models were checked for independence of residuals, over dispersion and presence of zero-inflated data [42]. In all cases the models were an adequate fit for the expected behaviour. All analyses were performed with the R computational environment [43] and the lme4 [44], bbmle [45], DHARMa [46] and ggplot2 [47] packages.

Results

A total of 15,764 mosquitoes belonging to 80 species/taxa in 15 genera were collected between March 2015 and April 2017 (Additional file 1). Among these, 6823 specimens of An. cruzii were identified, of which 781 individuals were collected in CDC traps (11.4% of the total) and the remainder in Shannon traps. Based on these data and the data from collections in Embura village and Evangelista de Souza between 2009 and 2010, An. cruzii represents 48% of all the mosquitoes collected in the Capivari-Monos EPA in these studies. However, the relative abundance of this species varied between the sites, ranging from approximately 5% in Embura village to more than 74% of all the mosquitoes collected in Evangelista de Souza. Average abundance of An. cruzii per human/hour in the Shannon traps varied from 0.5 in Embura to 57.5 in Evangelista de Souza. Average frequency of this species at ground level (based on CDC traps) varied from 0.18 (18%) in Evangelista de Souza to 0.58 (58%) in Embura village (Table 1).

Table 1 Results of mosquito collections at each site: total number of specimens, total number of Anopheles cruzii, relative abundance of An. cruzii, mean abundance of An. cruzii per human/hour in Shannon traps and mean proportion of An. cruzii collected at ground level in CDC traps

Turning to the landscape metrics, Evangelista de Souza and Cachoeira do Marsilac had higher values for forest cover and lower values for areas that had undergone anthropogenic changes, forest fragments and total edge length than the other sites. Table 2 shows the landscape metrics for the areas in the vicinity of each collection point.

Table 2 Landscape metrics for each of the five sites in the Capivari-Monos EPA

Among the mixed-effect models proposed to evaluate the relationship between landscape and An. cruzii abundance, the best fit was observed for the model that included the variable “forest cover 1 km” (∆AICc = 0, weight = 0.961), in which a positive predictive relationship was observed and an increase in forest cover led to an increase in An. cruzii abundance (Fig. 3). The second-best fit was for the model with the variable “total edge 500 m” (ΔAICc = 6.6, weight = 0.036), which had a negative predictive relationship with An. cruzii abundance. The null model showed the worst fit of all the proposed models (ΔAICc = 663.7, weight = < 0.001) (Table 3).

Fig. 3
figure 3

Observed values (points) and predicted values (lines) for abundance of Anopheles cruzii per human/hour as a function of the proportion of forest cover within a 1 km buffer of each site. The black line represents the average prediction of the model for the fixed effect. The colored lines represent the variations in the intercept due to random effects (collection month and year)

Table 3 Models proposed to predict the effect of landscape variations on Anopheles cruzii abundance (mean number of mosquitoes collected per human/hour) and acrodendrophily (measured as the proportion of mosquitoes collected at ground level)

In terms of the relationship between landscape and An. cruzii acrodendrophily, the model with the best fit was the one with the variable “total edge 500 m” (ΔAICc = 0, weight = 0.886), in which there was a positive predictive relationship between edge length and An. cruzii frequency at ground level (Fig. 4). The variable “forest cover 1 km” also showed a good fit (ΔAICc = 5.9, weight = 0.047) and a negative predictive relationship with An. cruzii frequency at ground level. This model had the second-best fit of those proposed. Again, the null model had one of the worst fits (ΔAICc = 31.5, weight = < 0.001); however, the difference in AICc between this model and the best model was not as high as the corresponding difference for the abundance models (Table 3).

Fig. 4
figure 4

Observed values (points) and predicted values (lines) for proportion of Anopheles cruzii in CDC traps at ground level as a function of the total edge length in kilometers within a 500 m buffer around each site. The black line represents the average prediction of the model for the fixed effect. The colored lines represent the variations in the intercept due to random effects (collection month and year)

Discussion

It was hypothesized in this study that factors associated with anthropogenic landscape changes may be responsible for the observed variation in An. cruzii abundance and acrodendrophily. The models with the best fits suggest that a loss of natural vegetation leads to a reduction in An. cruzii abundance, while an increase in the edge effect due to fragmentation and suppression of forest areas may favour greater activity by this species near ground level, which in turn may increase the contact rate between these mosquitoes and humans living on the edges of forest fragments where this species occurs.

As observed in many studies, An. cruzii tends to be the dominant, or at least one of the most abundant, species of Culicidae in the humid forests of the coastal slopes of Southeastern and Southern Brazil, and its population density is directly related to the abundance of bromeliads [16, 20, 23, 48, 49]. Recently, Chaves et al. [50] studied the influence of landscape on the distribution and abundance of An. cruzii and Anopheles bellator. They found that landscape configuration and composition seem to play a significant role in the abundance of both species, the former being more abundant in dense forest areas and the latter more prevalent in restinga (sandy, salty soil close to the sea covered with characteristic herbaceous plants). They also found that neither of these species occurs in rural environments, where bromeliads tend to be sparse or absent. Dorvillé [51] analysed several ecological studies on Culicidae in Southeastern Brazil and concluded that mosquitoes of the subgenus Kerteszia can be used as bioindicator species as they are highly susceptible to environmental degradation and their presence reflects a high degree of preservation. In studies between the years 2010 and 2013 in parks in built-up areas of the city of São Paulo, none of the specimens collected were from the Kerteszia subgenus although there were a significant number of bromeliads in isolated, preserved forest fragments in the study area [52, 53], another indicator of the susceptibility of this group to environmental changes. Previously, Ribeiro et al. [35] had also observed this difference in An. cruzii abundance between preserved and degraded areas of the Capivari-Monos EPA (Evangelista de Souza and Embura village, respectively).

There is a consensus among the various studies regarding the preference of An. cruzii for blood feeding in the tree canopy [2, 8, 14]. Nevertheless, variations in the biting activity of this mosquito in the canopy and at ground level have been observed in different locations, and few studies have sought to understand the factors that lead to these variations [20, 22,23,24, 26, 54]. Investigating the vertical stratification of mosquitoes in forest areas in the northeast of the state of São Paulo, Forattini et al. [20] observed differences in An. cruzii acrodendrophily between the study sites and an increase in acrodendrophily in nighttime collections. The authors suggested that the proximity of one of the sites to areas inhabited by humans may have led to a small reduction in acrodendrophily and that the increase in acrodendrophily during the night may have been related to the greater number of animals, particularly birds, at rest in the tree canopy.

In a study conducted between 1981 and 1982 in Serra dos Órgãos National Park in the state of Rio de Janeiro, Guimarães et al. [22] also observed an increase in acrodendrophily in this species at night and a more equitable distribution of individuals between the canopy and ground level during daylight hours. In the same study, the authors observed that throughout the year the average relative humidity of the upper stratum of the forest was lower than in the lower stratum, which could partly explain the preference of An. cruzii for flying and feeding in the treetops as the condensation of water vapor on the body of these insects may lead to a reduction in flight activity. In contrast, Deane [8] observed that in forest fragments in plateau regions of the Atlantic Forest farther inland, such as Cantareira State Park in São Paulo, An. cruzii showed marked acrodendrophilic behaviour, while in coastal mountains where the relative humidity is high, populations of this species seemed to feed more frequently near the ground.

These observations suggest that local variations in An. cruzii acrodendrophily are influenced by environmental and microclimatic factors, genetic variations between populations in different regions and the possibility that this mosquito actually represents a cryptic species complex [55,56,57]. The observations of the present study suggest that the acrodendrophilic behaviour of An. cruzii varies with changes in the composition and configuration of the surrounding landscape and that there is an increase in the blood-feeding activity of this mosquito in the lower stratum of the forest in response to an increase in edge effect in areas with greater human activity following deforestation and fragmentation of the native forest cover in these places. The results of the present study should be viewed with some caution as the number of sampling points was small and the models for the frequency of this mosquito at ground level do not have a much better fit than the null model (unlike the models used to explain the abundance). However, two factors lend support to our findings: the fact that the increased presence of human hosts appears to increase mosquito activity near ground level, and the fact that the microclimatic variations caused by the edge effect may also lead to greater number of mosquitoes near the ground.

The increase in human activity in areas on the edges of forest fragments or even inside the forest in itself represents an increase in the number of hosts available to the vector and to other mosquitoes in search of blood meal sources. It would, therefore, not be surprising to find that An. cruzii moves from the canopy to the ground more often in situations where humans and their domestic animals represent the nearest and most abundant sources of blood meals. Guimarães et al. [58] studied the feeding habits of mosquitoes collected in the canopy and at ground level in the Serra dos Órgãos National Park and found that most species, including An. cruzii, exhibited opportunistic behaviour, feeding on humans and other animals used as bait. Chaves et al. [59] reviewed several studies on mosquito feeding habits and used a null model of species co-occurrence to test whether mosquito feeding behaviour tended to be random (no preference for a particular species), segregated (certain mosquitoes feed only on a particular host) or aggregate (mosquitoes tend to feed primarily on a particular host). They found that in studies where data were collected at multiple sites there was a trend toward aggregate behaviour and concluded that contact between mosquitoes and hosts depends more on the availability of a given host than on an innate preference of the mosquito species.

As for the effect of microclimatic variations, factors such as luminosity, wind speed, temperature and humidity are known to vary from the edge to the interior of forest fragments [60,61,62]. In a study in the Atlantic Forest, Magnago et al. [63] showed that low humidity, high winds and higher temperatures are more common in edge environments than inside the forest. Based on this, a second hypothesis is that the reduction in humidity due to the edge effect in more fragmented areas may favour increased An. cruzii activity in the lower stratum, which agrees with the observations of Veloso et al. [64] and Guimarães et al. [65] regarding the effect of air humidity on the flight activity of Kerteszia mosquitoes.

The main limitations of the present study include the small number of sites investigated and the fact that these are all in the same area, the Capivari-Monos EPA. In addition, it was not possible to observe the number of mosquitoes collected per human/hour in the tree canopy because Shannon traps could not be set up on suspended platforms with human collectors for logistical reasons. Clearly, it is not known whether the results that would have been obtained had it been possible to set Shannon traps up with collectors would differ from those actually obtained; however, there can be no doubt that the presence of human collectors at ground level and in the canopy would be more appropriate for the present study as CDC traps do not fully mimic the presence of humans or simian hosts.

Although competing models were tested for two different landscape scales (500 m and 1 km), there was no prior evidence to believe that a larger landscape scale could better predict the variations in the abundance and acrodendrophily of An. cruzii than smaller, more local scale, or vice versa. Interestingly, the results suggest that a larger forest cover scale (1 km buffer) better predicts variations in the mosquito population abundance; while a smaller landscape scale (500 m buffer) better predicts variations in the acrodendrophily behaviour. However, it was not investigated how these variables would respond to larger or smaller landscape scales besides those measured. It has recently been shown that environmental variables measured at more local scales can have a great influence on spatial heterogeneity of the abundance of mosquitoes in forested urban environments [66]. In a more local scale the presence and quantity of bromeliads is a factor directly related to the abundance of An. cruzii [64]. Although this type of plant is quite common throughout the study area, the distribution and abundance of bromeliads at the sampling sites was not investigated, which may be considered one of the limitations of the present study. Nevertheless, such variations may be indirectly reflected in the abundance of adult mosquitoes that was lower in the more modified areas where bromeliads tend to be less abundant.

In a study about the infectivity of Anopheles mosquitoes at the Capivari-Monos EPA, Duarte et al. [11] found An. cruzii specimens naturally infected by both P. vivax and P. malariae. One of the areas where these were found was Embura village, where An. cruzii abundance was lower. Interestingly, species found in low abundance, such as Anopheles triannulatus, Anopheles strodei and Anopheles lutzi, were also found naturally infected by P. vivax and P. malariae in sites that had been subjected to anthropogenic change. As the role of these species of mosquitoes in malaria transmission in the Atlantic Forest is not yet known, other studies should be conducted in these areas to understand the vector–host transmission dynamics. Information on the prevalence of Plasmodium species in the An. cruzii populations collected in Cachoeira do Marsilac, the transition zone and Marsilac village was obtained from studies carried out in parallel with the present study. However, some of the laboratory work for these studies is still being carried out at the time of writing and data are therefore not yet available to compile with the data on abundance and acrodendrophily from the present study. This analysis should be done as part of a future study.

In the past, bromeliad malaria was an endemic disease of great epidemiological importance in the Southeastern and Southern regions of Brazil. Although it is now under control as a result of the substantial effort made, it has not been totally eliminated as there are still transmission foci in various places. Nevertheless, many of these foci probably do not come to the attention of the health authorities because of the high proportion of asymptomatic and oligosymptomatic cases. Furthermore, the symptoms can be confused with those of other diseases [29, 67,68,69]. Anthropogenic changes in the landscape and the consequent reduction in biodiversity are important factors in the emergence of malaria outbreaks in the Amazon region, where the phenomenon known as frontier malaria is now acknowledged to exist [70,71,72]. However, little is known about the effect of such changes on the dynamics of malaria transmission in the Atlantic Forest.

Conclusion

The data obtained in this study indicate that anthropogenic changes in the landscape lead to a reduction in the abundance of An. cruzii but can increase the contact rate between these mosquitoes and humans living at the edges of forest fragments where this species is found. Future studies should, therefore, seek to elucidate the effects of these landscape changes on the dynamics of Plasmodium transmission in the Atlantic Forest, which according to some studies includes the participation of simian hosts. The development of predictive models that seek to improve the understanding of how malaria vectors respond to changes in landscape composition and configuration can provide important information to assist planning and targeting of prevention and control actions.

Abbreviations

EPA:

Environmental Protection Area

CDC:

Centers for Disease Control and Prevention

AICc:

Akaike Information Criterion corrected for small samples

km:

kilometre

m:

metre

mm:

millimetre

°C:

degrees celsius

CO2 :

carbon dioxide

References

  1. Reitz R. Bromeliaceas: e a malaria-bromelia endemica. In: Herbário Barbosa Rodrigues, Flora Ilus Catarinense. Itajaí (SC), Brazil; 1983. p. 1–58.

  2. Forattini OP. Culicidologia Médica, v.2: Identificação, Biologia e Epidemiologia. edusp, Culicidologia Médica. São Paulo; 2002.

  3. Laporta GZ, Sallum MAM. Coexistence mechanisms at multiple scales in mosquito assemblages. BMC Ecol. 2014;14:30.

    Article  Google Scholar 

  4. Pina-Costa A, Brasil P, di Santi SM, de Araujo MP, Suárez-Mutis MC, Santelli ACF, et al. Malaria in Brazil: what happens outside the Amazonian endemic region. Mem Inst Oswaldo Cruz. 2014;109:618–33.

    Article  Google Scholar 

  5. Curado I, dos Santos Malafronte R, de Castro Duarte AMR, Kirchgatter K, Branquinho MS, Bianchi Galati EA. Malaria epidemiology in low-endemicity areas of the Atlantic Forest in the Vale do Ribeira, São Paulo, Brazil. Acta Trop. 2006;100:54–62.

    Article  Google Scholar 

  6. Cerutti C, Boulos M, Coutinho AF, Hatab MDCLD, Falqueto A, Rezende HR, et al. Epidemiologic aspects of the malaria transmission cycle in an area of very low incidence in Brazil. Malar J. 2007;6:33.

    Article  Google Scholar 

  7. Zavortink TJ. Mosquito studies (Diptera, Culicidae) XXIX. A review of the subgenus Kerteszia of Anopheles. Contrib Am Entomol Inst. 1973;9:1–54.

    Google Scholar 

  8. Deane LM. Simian malaria in Brazil. Mem Inst Oswaldo Cruz. 1992;87:1–20.

    Article  Google Scholar 

  9. Curado I, Duarte AMRC, Lal AA, Oliveira SG, Kloetzel JK. Antibodies anti bloodstream and circumsporozoite antigens (Plasmodium vivax and Plasmodium malariae/P. brasilianum) in areas of very low malaria endemicity in Brazil. Mem Inst Oswaldo Cruz. 1997;92:235–43.

    Article  CAS  Google Scholar 

  10. Rezende HR, Soares RM, Cerutti C Jr, Alves IC, Natal D, Urbinatti PR, et al. Entomological characterization and natural infection of anophelines in an area of the Atlantic Forest with autochthonous malaria cases in mountainous region of Espírito Santo State, Brazil. Neotrop Entomol. 2009;38:272–80.

    Article  Google Scholar 

  11. Duarte AMRC, Pereira DM, de Paula MB, Fernandes A, Urbinatti PR, Ribeiro AF, et al. Natural infection in anopheline species and its implications for autochthonous malaria in the Atlantic Forest in Brazil. Parasit Vectors. 2013;6:58.

    Article  Google Scholar 

  12. Neves A, Urbinatti PR, dos Santos Malafronte R, Fernandes A, da Silva Paganini W, Natal D. Malaria outside the Amazon region: natural Plasmodium infection in anophelines collected near an indigenous village in the Vale do Rio Branco, Itanhaem, SP, Brazil. Acta Trop. 2013;125:102–6.

    Article  Google Scholar 

  13. Laporta GZ, Ramos DG, Ribeiro MC, Sallum MAM. Habitat suitability of Anopheles vector species and association with human malaria in the Atlantic Forest in south-eastern Brazil. Mem Inst Oswaldo Cruz. 2011;106:239–45.

    Article  Google Scholar 

  14. Marrelli MT, Malafronte RS, Sallum MAM, Natal D. Kerteszia subgenus of Anopheles associated with the Brazilian Atlantic rainforest: current knowledge and future challenges. Malar J. 2007;6:127.

    Article  Google Scholar 

  15. Consoli RAGB, Lourenço-de-Oliveira R. Principais mosquitos de importância sanitária no Brasil. Rio de Janeiro: FIOCRUZ; 1994.

    Book  Google Scholar 

  16. Forattini OP, Alves AC, Natal D, Santos JL. Observações sobre atividade de mosquitos Culicidae em mata primitiva da encosta no Vale do Ribeira, São Paulo, Brasil. Rev Saude Publica. 1986;20:1–20.

    Article  CAS  Google Scholar 

  17. Tubaki RM, Carreri-Bruno GC, Glasser CM, Ciaravolo RMC. Biting activity of Anopheles (Kerteszia) cruzii (Diptera, Culicidae) in domiciliary habitats in the southern Atlantic Forest, Peruibe, State of São Paulo, Brazil. Rev Bras Entomol. 1993;37:569–75.

    Google Scholar 

  18. Guimarães AE, Gentile C, Lopes CM, Sant’Anna A. Ecologia de mosquitos em áreas do Parque Nacional da Serra da Bocaina. II—Freqüencia mensal e fatores climáticos. Rev Saude Publica. 2001;35:392–9.

    Article  Google Scholar 

  19. Garnham PCC, Harper JO, Highton RB. The mosquitos of the Kaimosi Forest, Kenya colony, with special reference to yellow fever. Bull Entomol Res. 1946;36:473–96.

    Article  CAS  Google Scholar 

  20. Forattini OP, Lopes OS, Rabello EX. Investigações sobre o comportamento de formas adultas de mosquitos silvestres no Estado de São Paulo, Brasil. Rev Saude Publica. 1968;2:111–73.

    Article  CAS  Google Scholar 

  21. Deane LM, Ferreira Neto JA, Lima MM. The vertical dispersion of Anopheles (Kerteszia) cruzi in a forest in southern Brazil suggests that human cases of malaria of simian origin might be expected. Mem Inst Oswaldo Cruz. 1984;79:461–3.

    Article  CAS  Google Scholar 

  22. Guimarães AÉ, Arlé M, Machado RNM. Mosquitos no Parque Nacional da Serra dos Órgãos, Estado do Rio de Janeiro, Brasil. II. Distribuição vertical. Mem Inst Oswaldo Cruz. 1985;80:171–85.

    Article  Google Scholar 

  23. Ueno HM, Forattini OP, Kakitani I. Distribuição vertical e sazonal de Anopheles (Kerteszia) em Ilha Comprida, SP. Rev Saude Publica. 2007;41:269–75.

    Article  Google Scholar 

  24. Buery JC, Rezende HR, Natal L, Da Silva LS, De Menezes RMT, Fux B, et al. Ecological characterisation and infection of anophelines (Diptera: Culicidae) of the Atlantic Forest in the southeast of Brazil over a 10 year period: has the behaviour of the autochthonous malaria vector changed? Mem Inst Oswaldo Cruz. 2018;113:11–118.

    Article  Google Scholar 

  25. Deane LM, Deane MP, Ferreira Neto J. Studies on transmission of simian malaria and on a natural infection of man with Plasmodium simium in Brazil. Bull World Health Organ. 1966;35:805.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Deane LM, Deane MP, Ferreira Neto J, Almeida FB. On the transmission of simian malaria in Brazil. Rev Inst Med Trop Sao Paulo. 1971;13:311–9.

    CAS  PubMed  Google Scholar 

  27. Deane LM. Estudo sobre a malária simiana no estado do Espirito Santo. Ver Bras Biol. 1968;28:531–6.

    CAS  Google Scholar 

  28. Deane LM, Ferreira Neto JA. Malaria em macacos do estado do Rio Grande do Sul. Observações preliminares. Rev Inst Med Trop Sao Paulo. 1969;11:299–305.

    CAS  PubMed  Google Scholar 

  29. Brasil P, Zalis MG, de Pina-Costa A, Siqueira AM, Júnior CB, Silva S, et al. Outbreak of human malaria caused by Plasmodium simium in the Atlantic Forest in Rio de Janeiro: a molecular epidemiological investigation. Lancet Glob Heal. 2017;15:e1038–46.

    Article  Google Scholar 

  30. Couto RD, Latorre MDRDDO, Di Santi SM, Natal D. Autochthonous malaria notified in the State of São Paulo: clinical and epidemiological characteristics from 1980 to 2007. Rev Soc Bras Med Trop. 1980;2010(43):52–8.

    Google Scholar 

  31. CVE—Centro de Vigilância Epidemiológica. Casos Autóctones de Malária no Estado de São Paulo. 2018. http://www.saude.sp.gov.br/cve-centro-de-vigilancia-epidemiologica-prof.-alexandre-vranjac/zoonoses/malaria/dados-estatisticos. Accessed 25 Nov 2018.

  32. São Paulo (municipality). Secretaria do Verde e Meio Ambiente. Plano de Manejo—Área de Proteção Ambiental (APA) Capivari-Monos. São Paulo; 2011. http://www.prefeitura.sp.gov.br/cidade/secretarias/meio_ambiente. Accessed 20 Sept 2018.

  33. Instituto de Astronomia, Geofísica e Ciências Atmosféricas da Universidade de São Paulo. Nascer e ocaso do Sol (2011–2025). 2019. http://www.iag.usp.br/astronomia/nascer-e-ocaso-do-sol. Accessed 4 Mar 2019.

  34. Chaves LF. Globally invasive, withdrawing at home: Aedes albopictus and Aedes japonicus facing the rise of Aedes flavopictus. Int J Biometeorol. 2016;60:1727–38.

    Article  Google Scholar 

  35. Ribeiro AF, Urbinatti PR, De Castro Duarte AMR, De Paula MB, Pereira DM, Mucci LF, et al. Mosquitoes in degraded and preserved areas of the Atlantic Forest and potential for vector-borne disease risk in the municipality of São Paulo, Brazil. J Vector Ecol. 2012;37:316–24.

    Article  Google Scholar 

  36. Lane J. Neotropical Culicidae, vol. 2. São Paulo: University of São Paulo; 1953.

    Google Scholar 

  37. Ferreira E, Corrêa R, Tomich A, Sa F. Study of the range of flight of Anopheles (Kerteszia) cruzii and of Anopheles (Kerteszia) bellator in Guaratuba, on the coast of the State of Parana, Brazil. Rev Bras Malariol Doencas Trop. 1969;21:819–22.

    CAS  PubMed  Google Scholar 

  38. McGarigal K, Cushman SA, Neel MC, Ene E. FRAGSTATS v4: Spatial Pattern Analysis Program for Categorical and Continuous Maps. Univ Massachusettes, Amherst, MA; 2012. http://www.umass.edu/landeco/research/fragstats/fragstats.html. Accessed 22 Sept 2018.

  39. Chaves LF. An entomologist guide to demystify pseudoreplication: data analysis of field studies with design constraints. J Med Entomol. 2010;47:291–8.

    Article  Google Scholar 

  40. Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.

    Google Scholar 

  41. Johnson JB, Omland KS. Model selection in ecology and evolution. Trends Ecol Evol. 2004;19:101–8.

    Article  Google Scholar 

  42. Dupuy J-F. Statistical methods for overdispersed count data. 1st ed. Amsterdam: Elsevier; 2018.

    Google Scholar 

  43. R Development Core Team. R—a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018. http://www.R-project.org. Accessed 22 Sept 2018.

  44. Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, et al. Package lme4: Linear mixed-effects models using Eigen and S4. R Found Stat Comput. 2016;1:1–113.

    Google Scholar 

  45. Bolker B, Team RDC. bbmle: tools for general maximum likelihood estimation. R Packag. version. 2018.

  46. Hartig F. DHARMa: residual diagnostics for hierarchical (multi-level/mixed) regression models. R Packag. version 0.2.0. 2018.

  47. Wickham H. ggplot2: elegant graphics for data analysis. Berlin: Springer; 2016.

    Book  Google Scholar 

  48. Bona A, Navarro-Silva M. Diversidade de Culicidae durante os períodos crepusculares em bioma de Floresta Atlântica e paridade de Anopheles cruzii (Diptera: Culicidae). Rev Bras Zool. 2008;25:40–8.

    Article  Google Scholar 

  49. Guedes MLP, Navarro-Silva MA. Mosquito community composition in dynamic landscapes from the Atlantic Forest biome (Diptera, Culicidae). Rev Bras Entomol. 2014;58:88–94.

    Article  Google Scholar 

  50. Chaves LSM, de Sá ILR, Bergamaschi DP, Sallum MAM. Kerteszia Theobald (Diptera: Culicidae) mosquitoes and bromeliads: a landscape ecology approach regarding two species in the Atlantic rainforest. Acta Trop. 2016;164:303–13.

    Article  Google Scholar 

  51. Dorville LFM. Mosquitoes as bioindicators of forest degradation in Southeastern Brazil, a statistical evaluation of published data in the literature. Stud Neotrop Fauna Environ. 1996;31:68–78.

    Article  Google Scholar 

  52. Medeiros-Sousa AR, Ceretti W, Urbinatti PR, de Carvalho GC, de Paula MB, Fernandes A, et al. Mosquito fauna in municipal parks of São Paulo City, Brazil: a preliminary survey. J Am Mosq Control Assoc. 2013;29:275–9.

    Article  Google Scholar 

  53. Ceretti-Junior W, de Oliveira Christe R, Rizzo M, Strobel RC, de Matos Junior MO, de Mello MHSH, et al. Species composition and ecological aspects of immature mosquitoes (Diptera: Culicidae) in bromeliads in urban parks in the city of São Paulo, Brazil. J Arthropod Borne Dis. 2016;10:102–12.

    PubMed  Google Scholar 

  54. dos Santos-Neto LG, Lozovei AL. Aspectos ecológicos de Anopheles cruzii e Culex ribeirensis (Diptera, Culicidae) da Mata Atlântica de Morretes, Paraná, Brasil. Rev Bras Entomol. 2008;52:105–11.

    Article  Google Scholar 

  55. Rona LD, Carvalho-Pinto CJ, Gentile C, Grisard EC, Peixoto AA. Assessing the molecular divergence between Anopheles (Kerteszia) cruzii populations from Brazil using the timeless gene: further evidence of a species complex. Malar J. 2009;8:60.

    Article  Google Scholar 

  56. Dias GR, Fujii TTS, Fogel BF, Lourenço-De-Oliveira R, Silva-Do-Nascimento TF, Pitaluga AN, et al. Cryptic diversity in an Atlantic Forest malaria vector from the mountains of South-East Brazil. Parasit Vectors. 2018;11:36.

    Article  Google Scholar 

  57. Multini LC, Wilke ABB, Marrelli MT. Urbanization as a driver for temporal wing-shape variation in Anopheles cruzii (Diptera: Culicidae). Acta Trop. 2019;190:30–6.

    Article  Google Scholar 

  58. Guimarães AE, Arlé M, Machado RNM. Mosquitos no Parque Nacional da Serra dos Orgãos, Estado de Rio do Janeiro, Brasil. IV. Preferencia alimentar. Mem Inst Oswaldo Cruz. 1987;82:277–85.

    Article  Google Scholar 

  59. Chaves LF, Harrington LC, Keogh CL, Nguyen AM, Kitron UD. Blood feeding patterns of mosquitoes: random or structured? Front Zool. 2010;7:3.

    Article  Google Scholar 

  60. Chen J, Franklin JF, Spies TA. Growing-season microclimatic gradients from clearcut edges into old-growth Douglas-fir forests. Ecol Appl. 1995;5:74–86.

    Article  Google Scholar 

  61. Didham RK, Lawton JH. Edge structure determines the magnitude of changes in microclimate and vegetation structure in tropical forest fragments. Biotropica. 1999;31:17–30.

    Google Scholar 

  62. Ewers RM, Banks-Leite C. Fragmentation impairs the microclimate buffering effect of tropical forests. PLoS ONE. 2013;8:e58093.

    Article  CAS  Google Scholar 

  63. Magnago LFS, Rocha MF, Meyer L, Martins SV, Meira-Neto JAA. Microclimatic conditions at forest edges have significant impacts on vegetation structure in large Atlantic forest fragments. Biodivers Conserv. 2015;24:2305–18.

    Article  Google Scholar 

  64. Veloso HP, Fontana Junior P, Klein RM, Siqueira Jaccoud RJ. Os anofelinos do sub-gênero Kerteszia em relação à distribuição das bromeliáceas em comunidades florestais do município de Brusque, Estado de Santa Catarina. Mem Inst Oswaldo Cruz. 1956;54:1–86.

    Article  Google Scholar 

  65. Guimarães AÉ, Victório VMN. Mosquitos no Parque Nacional da Serra dos Orgãos, Estado do Rio de Janeiro, Brasil. III. Preferência horária para hematofagia. Mem Inst Oswaldo Cruz. 1986;81:93–103.

    Article  Google Scholar 

  66. Chaves LF, Imanishi N, Hoshi T. Population dynamics of Armigeres subalbatus (Diptera: Culicidae) across a temperate altitudinal gradient. Bull Entomol Res. 2015;105:589–97.

    Article  CAS  Google Scholar 

  67. Brasil P, Costa AP, Longo CL, Da Silva S, Ferreira-Da-Cruz MF, Daniel-Ribeiro CT. Malaria, a difficult diagnosis in a febrile patient with sub-microscopic parasitaemia and polyclonal lymphocyte activation outside the endemic region, in Brazil. Malar J. 2013;12:402.

    Article  Google Scholar 

  68. Maselli LM, Levy D, Laporta GZ, Monteiro AM, Fukuya LA, Ferreira-Da-Cruz MF, et al. Detection of Plasmodium falciparum and Plasmodium vivax subclinical infection in non-endemic region: implications for blood transfusion and malaria epidemiology. Malar J. 2014;13:224.

    Article  Google Scholar 

  69. Alencar FEC, Malafronte RDS, Cerutti Junior C, Natal Fernandes L, Buery JC, Fux B, et al. Assessment of asymptomatic Plasmodium spp. infection by detection of parasite DNA in residents of an extra-Amazonian region of Brazil. Malar J. 2018;17:113.

    Article  Google Scholar 

  70. Sawyer DR. Frontier malaria in the Amazon Region of Brazil: types of malaria situations and some implications for control. Brasília: PHO/WHO/TDR; 1988.

    Google Scholar 

  71. Castro MC, Monte-Mor RL, Sawyer DO, Singer BH. Malaria risk on the Amazon frontier. Proc Natl Acad Sci USA. 2006;103:2452–7.

    Article  Google Scholar 

  72. Chaves LSM, Conn JE, López RVM, Sallum MAM. Abundance of impacted forest patches less than 5 km2 is a key driver of the incidence of malaria in Amazonian Brazil. Sci Rep. 2018;8:7077.

    Article  Google Scholar 

Download references

Authors’ contributions

Study design: ARM-S, MTM, AMRCD, LFM. Field collections: ARM-S, WC-J, ROC. Data analysis and scientific writing: ARM-S. Scientific review: MTM, AMRCD, LFM, WC-J, ROC. All authors read and approved the final manuscript.

Acknowledgements

We would like to express our gratitude to the field and laboratory teams at the Superintendency for the Control of Endemic Diseases, São Paulo Zoonosis Control Center, and the School of Public Health, São Paulo University: João Carlos do Nascimento, Paulo Frugoli dos Santos, Luis Milton Bonafé, Antônio Waldomiro de Oliveira, Laércio Molinari, Gabriel Marcelino Neto, Luiz Sposito Jr, Renildo Souza Teixeira, Aristides Fernandez, Daniel Pagotto Vendrami, Gabriela Cristina de Carvalho, Ramon Wilk da Silva, Eduardo Evangelista de Souza and Amanda Alves Camargo.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

This work was funded by the São Paulo Research Foundation (FAPESP) Grants Nos. 2014/50444-5 and 2014/10919-4. ARM-S was supported by the São Paulo Research Foundation (FAPESP) Ref. No. 2015/18630-6.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Antônio Ralph Medeiros-Sousa.

Additional file

Additional file 1.

Species and number of individuals collected in the Capivari-Monos Environmental Protection Area by collection. Collections made from March 2015 to April 2017.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Medeiros-Sousa, A.R., de Oliveira Christe, R., de Castro Duarte, A.M.R. et al. Effects of anthropogenic landscape changes on the abundance and acrodendrophily of Anopheles (Kerteszia) cruzii, the main vector of malaria parasites in the Atlantic Forest in Brazil. Malar J 18, 110 (2019). https://doi.org/10.1186/s12936-019-2744-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-019-2744-8

Keywords