 Research article
 Open Access
 Published:
Heterogeneities in dengue spatialtemporal transmission in Brazilian cities and its influence on the optimal age of vaccination
BMC Public Health volume 19, Article number: 155 (2019)
Abstract
Background
The development of a safe and effective vaccine is considered crucial for dengue transmission control since vetor control has been failed; some potential candidates are currently in test, and in this context theoretical studies are necessary to evaluate vaccination strategies such as the age groups that should be vaccinated, the percentage of the population at risk, and the target geographic regions to make dengue control feasible and optimal.
Methods
A partial differential model is used to mimics dengue transmission in human population in order to estimate the optimal vaccination age, using data collected from dengue reported cases in ten cities of Brazil from 2001 to 2014. For this purpose, the basic reproduction number of the disease was minimized assuming a singledose vaccination strategy, equal vaccine efficacy for all circulating serotypes, and no vaccine failure. Numerical methods were used to assess the optimal vaccination age and its confidence age range.
Results
The results reveal complex spatialtemporal patterns associated to the disease transmission, highlighting the heterogeneity in defining the target population for dengue vaccination. However, the values obtained for the optimal age of vaccination, as targeting individuals under 13 years old, are compatible with the ones reported in similar studies in Brazil. The results also show that the optimal age for vaccination in general does not match with the age of the highest number of cases.
Conclusions
The variation of the optimal age for vaccination across the country reflects heterogeneities in dengue spatialtemporal transmission in Brazilian cities, and can be used to define the target population and cities to optimize vaccination strategies in a context of high cost and low quantity of available vaccine.
Background
Control methods of infectious diseases must take into account ecological, immunological, and behavioral aspects of the disease agent, as well as its transmission mode and the types of hosts. For this end, analysis of spatialtemporal patterns of infectious diseases reports, together with laboratory and field experiments, are important tools for understanding diseases transmission, and designing control strategies, such as sanitary measures, health policies, vaccination, or vector control in the case of vectorborne diseases.
One of the most efficient control tools has been the vaccine. Several diseases that plagued humanity for centuries, such as poliomyelitis and measles, have been currently controlled in some regions or countries, thanks to vaccination routine or campaigns. The most notorious example is smallpox, which was eradicated in 1979, by a collaborative global vaccination program [1].
For vectorborne diseases, the control is based mainly on the reduction of the vector population. In particular, for dengue, all efforts have been focused on Aedes mosquitoes control, by means of the use of insecticides that eliminate adult stage, larvicides and destruction of oviposition containers to remove immature stage, and by the use of biological techniques [2–4]. Unfortunately, increasing urbanization, waste generation, human displacement, mosquito resistance to insecticides, as well as mosquito fast spread across world, have contributed to vector control failure [2].
The development of safe and effective vaccines is considered crucial for better expectations for dengue control. The principal obstacle is the existence of four different serotypes, and the potential risk of dengue hemorrhagic fever (DHF) associated to secondary infections with heterologous serotype [5]. Due to these denguespecific complexities, vaccine development focuses on the generation of a tetravalent vaccine aimed to provide longterm protection against all virus serotypes. Some vaccine candidates are currently in test, and only one was registered until now in some countries [6]. It is worth to mention that, at the end of 2017, there was an increment in the number of hospitalizations for severe dengue in vaccinated individuals who had never had dengue before [7].
Theoretical studies, such as mathematical modeling, are useful to provide tools for the development of vaccination strategies that aim to promote herd immunity. Using dengue data, these models can be useful to decide which are the age groups that should be vaccinated, the percentage of the population at risk that must be vaccinated to control disease transmission, and the target geographic regions to optimize disease control.
Mathematical models have been proposed to evaluate possible dengue vaccination strategies using different techniques. Billings et al. (2008) [8] used a system of ordinary differential equations (ODE) to evaluate the efficacy of a singlestrain vaccine assuming antibodydependent enhancement. The authors considered several possible scenarios to evaluate the vaccine efficacy in the presence of two serotypes, resulting in a diagram showing effective vaccination rates versus strains persistence and extinction. Amaku et al. (2012) [9] used dengue serological data from Recife City, Brazil, and a system of timedelayed differential equations to estimate the optimal vaccination age; they found that it should be vaccinated children between 3 and 14 years, and 80% of vaccination coverage has to be achieved. On the other hand, an agentbased model was developed by Chao et al. [10] to simulate the epidemiology of dengue transmission in a semirural area of Thailand; they obtained that, for a fixed number of doses, vaccinating children from 2 to 14 years old would reduce dengue infection in the total population more than covering both children and adults (2 to 46 years old).
An agestructured multistrain model was done by the authors of reference [11] to design scenarios for the potential impact of a dengue vaccine on a population. Using data from Southern Vietnam, the authors showed that seasonality and short crossprotection against infection ranging from 6 to 17 months are a keystone to produce the observed disease periodicity. Also they argued that vaccination reduces disease burden by reducing the magnitude and frequency of outbreaks. The same kind of model was used to set up the vaccination strategy that minimizes the incidence of DHF [12]: in Thailand the optimal strategy is to vaccine children from 0.5 to 12 years old while in Brazil, it is better to vaccine adults from 18 to 34 years old. Recently, using a system of partial differential equations for human population, and delay differential equations for vector population, as well as dengue incidence data in Brazil, Maier et al. [13] provided an estimation of the optimal vaccination age when different assumptions for vaccine efficacy and risk of infection have been taking into account; minimizing hospitalization or mortality due to dengue as measures of risk of dengue infection, wide ranges of values of optimal vaccination ages have been obtained, regarding both the serotypes in circulation and the model assumptions.
The objective of this study is to estimate the priority vaccination age range against the four serotypes of dengue for several Brazilian cities that had heterogeneities in spatialtemporal dengue transmission for the period of 2001 to 2014. For this end, we use the partial differential equations (PDE) model proposed by CruzPacheco et al. [14]. The variation of the optimal age for vaccination across the country can be used to define the target population, taking into account dengue epidemics in that period in different cities, to optimize vaccination strategies in a context of high cost and low quantity of available vaccine.
Methods
The method we used to evaluate the optimal vaccination age [14] in ten selected Brazilians cities, as well as in the whole country, is based on continuous compartamental models. These models are applied for large enough populations, since they are based on the assumption that susceptible and infectious are well mixed. Also vectorborne diseases have a higher prevalence in regions with high density of human and vectors. For these reasons, we establish the following criteria to select the cities used in the present study: (i) the city has to have more than 500,000 inhabitants; (ii) it has to present a cumulative incidence of dengue fever during the period of 2001 to 2014 greater than 500 per 100,000 inhabitants.
The number of dengue cases by year (2001 to 2014) and by age (< 1 year, 1 to 4 years, 5 year interval from 5 to 59 years, > 60 years) were obtained from Brazilian Notified Disease Information System (SINAN). The midpoints of these age classes were used in the analyses.
In the epidemiological model, the human population is divided into Susceptible, Infected and Recovered (SIR model), and the vector population is divided into Susceptible and Infected (SI model). The coupling between the two population is done by the transmission rate between infected humans (or susceptible vectors) and susceptible vectors (or infected humans).
Supposing that the increase of human protection under an agevaccination distribution is equivalent to the decrease of the net reproductive number R_{0}(V) constrained by the cost of vaccination, CruzPacheco et al. [14] showed that minimizing R_{0}(V) is equivalent to maximize the function H(a) given by:
where f(a) is the probability that a susceptible human of age a acquires dengue virus from an infected vector and it is obtained from the distribution of dengue cases by age. In the above expression, K is given by
where b is the mosquito biting rate, μ_{v} is the vector mortality rate, α_{v} is the infection rate from vectors to humans, μ_{h} is the human mortality rate, γ_{h} is the human recovery rate, and c is the cost associated to the vaccination schedule. Finally, N_{v} and N_{h} are, respectively, the total number of vectors and humans whose values depend on each city. As in reference [14], we assume that b=0.5 days ^{−1}, α_{v}=0.75, μ_{h}=0.00004 days ^{−1}, c=1, N_{v}/N_{h}=1, μ_{v}=0.06 days ^{−1}, and γ_{h}=0.125 days ^{−1}.
The optimal vaccination age corresponds to the age for which H(a) assumes its maximum value. All parameters in the expression for K are constant (do not depend on age a), therefore their values are irrelevant for the evaluation of the optimal vaccination age. The exception is the μ_{h} parameter that also appears inside of the integral multiplying a; the optimal vaccination age can be sensitive to it. We are assuming that everybody is vaccinated at the same age, vaccine efficacy is constant along the age classes, and the same for all serotypes; also that vaccineinduced immunity is lifelong.
Using data obtained from SINAN [15], the optimal vaccination age for each city in each year was obtained by the following algorithm: (i) evaluation of f(a) from the proportion of dengue reported cases; (ii) approximation of the integral that appears in equation (1) by the trapezoidal rule resulting in a curve of H(a) that is a function of a; (iii) evaluation of the derivative of this curve by the method of finite differences to obtain another curve that is also a function of a, h(a):=dH(a)/da; iv) the use of Newton’s method to find the roots of h(a). As we have a discrete set of data instead of a curve, we used two different approaches: 1) a linear interpolation in the interval of h(a_{1}) and h(a_{2}) where h(a_{1})h(a_{2})<0, to estimate \({\bar a}\in ]a_{1},a_{2}[\), the optimal vaccination age; 2) a cubic spline to interpolate all the curve before searching for h(a_{1})h(a_{2})<0. The results obtained from these two approaches were compared to define the confidence interval of the optimal vaccination range. We decided to report the results obtained from the linear interpolation, except in cases where the linear interpolation failed on getting the optimum value of \({\bar a}\).
In Brazil, the serology for dengue is done in a small fraction of the reported cases during the epidemic and nonepidemic period, when dengue incidence is below than 300 per 100.000 habitants [16]. After that, all diagnoses are clinical, based on the individual symptoms, and area of circulation of the dengue virus. Overall, it is estimated from the literature that the percentage of dengue asymptomatic individuals can reach up to 80% of the cases [17]. Therefore, we assumed that the age distribution of dengue cases is not perfect, and we reported the optimal vaccination age as well as its confidence age range. In order to construct a confidence interval for the optimal vaccination age, we generated 100,000 samples obtained from the dengue reported cases, assuming an error of 15% in each one of them. This value of 15% was chosen after a systematic study that involved the error coming from the use of numerical methods (derivation, integration, root finding, and interpolation) to find the optimal vaccination age. Besides the uncertainty on the reported cases of dengue, variation on the human mortality rate can affect the optimal vaccination age (see Eqs. (1) and (2)). In order to investigate that effect, the same procedure was applied to estimate the optimal vaccination age, as well as its confidence age range, assuming different values for human mortality rate. Since from 2001 to 2014, the life expectancy in Brazil changes from 70.5 to 75 years [18], we used the minimal (Table 1) and the maximal (Table 2) values of life expectancy in the estimations for each city at each year.
Results
Only 38 from 5565 cities, corresponding to 0.7% of the Brazilian cities, meet the minimum population criterion described in the section of Methods. Among them, 10 from the 38 cities are selected since they also follow the criterion of having cumulative incidence above or equal to 500 per 100,000 habitants during the studied period. The chosen cities belong to different states and they are localized in three of the five geographic regions of Brazil (see Fig. 1): Rio de Janeiro (RJ), Ribeirão Preto (RP), Campinas (Ca) and Belo Horizonte (BH) belong to the Southwest Region; Fortaleza (Fo), Natal (Na) and Maceió (Ma) are localized in the Northeast Region; Goiânia (Go), Cuiabá (Cu) and Campo Grande (CG) are from the Midwest Region. These cities are responsible for 22.8% (1,779,350/7,818,624) of dengue cases reported over the 14 years period.
Figure 2 illustrates the age distribution of dengue cases in the whole country and in three other studied cities selected from each region of the country. Distinct patterns can be observed among the cities and for the same city in different years. In particular, we can see that the age profile for the proportion of dengue cases in Goiânia (Fig. 2b) is homogeneous during the 14 years period. Also, 25% of the cases occur in individuals that are under 19 years old. For Fortaleza (Fig. 2c), we can see that between 2005 and 2010 there was a substantial increase of the proportion of dengue cases over small age classes (until 14 years old). Among the studied cities, Rio de Janeiro (Fig. 2d) had the most heterogeneous age profile for the proportions of dengue cases.
In order to illustrate the procedure used to obtain the optimal vaccination age and its confidence interval we show in Fig. 3 the proportion of dengue cases by age, the curve H(a) as a function of age, and the derivative of the curve H(a) with respect to age, h(a). The three curves correspond to the proportion of dengue cases in Maceió reported in 2014. For this data set, we estimated the value of the optimal vaccination age as 8.0 years (Fig. 3b) and from 6.7 to 9.3 its confidence age range (Fig. 3c).
Table 1 summarizes the results obtained for the optimal age and its confidence age range across the cities and years assuming μ_{h}=1/70.5 years ^{−1}. We also present the result for the whole country. The confidence age range varies from less than 1 year (Campo Grande and Cuiabá, 2004, and Fortaleza, from 2007 to 2011) to 20 years in Campinas (2005). The optimal vaccination age varies from 2 years old to 19 years old in Campinas (2004). Concerning the whole country, the minimum and maximal values for the confidence intervals should be from 2 years old (2008) to 10 years old (2001, 2004, 20132014), meanwhile the optimal age is from 3 years old (2008) to 8 years old (2001, 20032004, 20132014).
Table 2 shows the confidence age range of the optimal vaccination age across the cities and years assuming the other extreme value (1/75.0 years ^{−1}) of the human mortality rate μ_{h}. The confidence age range varies from less than 1 year old (Campo Grande, 2004, Fortaleza, from 2007 to 2011, and Campinas, 2001) to 20 years old in Campinas (2005). The optimal vaccination age varies from 2 years old to 18 years old in Campinas (2004). Concerning the whole country, the only difference regards the maximal value of the confidence interval that decrease from 10 (<10.5) years old in Table 1 to 9 (> 9.5) years old in Table 2 (2001, 2004, 20132014).
Figure 4 shows the behaviour of the optimal vaccination age for the ten Brazilian cities during the study period (2001 to 2014) through a box plot representation of the optimal age exhibited in Table 1. Although the results exhibited in Table 2 are slightly different for optimal age and confidence interval, the variation from the minimal to the maximal value of μ_{h} does not produce significant changes in Fig. 4, which compiles, after all, the optimal age and its dispersion and skewness for each city taking into account the studied period. Concerning the different scenarios of the studied cities, Fig. 4 also shows that the median value of the optimal vaccination age is around 9 years old in Belo Horizonte, Campinas, Ribeirão Preto, Campo Grande and Goiânia; moreover their dispersion and skeewess are smaller than for the rest of the studied cities. On the other hand, the median value of the optimal vaccination age is around 6 years old in Fortaleza, Natal, Rio de Janeiro and Cuiabá; the largest dispersion is observed in Natal, meanwhile Maceió presents the largest skeewess.
Discussion
In this work, a partial differential equation model with agestructure was used to evaluate the optimal vaccination age in ten selected Brazilian cities as well as in the whole country. It is worth to note that, currently, an efficient vaccine with no restriction does not exist. The vaccine licensed in some countries, including Brazil, is indicated to aged from 9 to 45 years old, and recommended only to individuals with previous dengue infection, due to the highest risk of hospitalizations in dengue naive individuals [19, 20].
Considering scenarios where everybody is vaccinated at the same age, vaccine efficacy is constant along the age classes, and the same for all serotypes; the priority age range of vaccination against dengue, in the 10 Brazilian cities of this study, varies from less than one year old to 20 years old. The optimal vaccination age indicates that vaccination has to target individuals under 19 years old according to the results presented in Table 1 (that suffered a slightly reduction to 18.7 years old in Table 2). If the outliers are not considered, according to Fig. 4, vaccination has to target individuals under 13 years old; therefore we have assumed age of 13 as the maximum age for dengue vaccination.
As far as we know, this is the first time that a mathematical model, based on PDE, is applied to dengue epidemics using actual age distribution of cases in different cities, promoting a general and comparative scenario of many epidemics in those Brazilian cities during a time interval of more than 10 years (from 2001 to 2014). The variation of the optimal age for vaccination across the country highlights the complex spatialtemporal behaviour of the disease, which is evident in the results illustrated in Fig. 4 as well as in Table 1. For instance, in the case of Fortaleza, comparing Fig. 2c to results presented in Table 1 and Fig. 4, the decreasing of optimal age seems to be related to the increasing of children infected by dengue, which is also observed in other cities of the Northeast of Brazil [21].
The results are compatible with others from the literature. For instance, the authors in reference [9] applied an age dependent seroprevalence model to dengue data of a Northeast Brazil city, founding a similar age range for vaccination (3 to 14 years old). Due to lack of age distribution of dengue cases by serotype for each studied city, our model does not consider the four serotypes. However, our results are also compatible with others from the literature that considers cocirculation of dengue serotypes, since the authors obtain that the target age range is from under 1 year to 9 years old, for the whole country [13]. That comparison was allowed due to the values of optimal age (from 3 to 8 years old) and confidence interval (from 2 to 10 years old) estimated for the whole country in Table 1. Although that result for the whole country should not be applied for all cities with dengue cases, it could be a coarse approximation for other cities (besides the ten studied cities) that have large population and large incidence. It is a kind of ’spatial’ mean field result to provide a general view of the country with the aim of comparing with the results in the literature.
The model, as well as the method of estimating an optimal age, was previously applied to dengue data in Mexico (the whole country) in 2014 using an analytical fitted function for the age distribution [14]. Since the PDE model allows the variables to change continuously with the age, it is possible to calculate the age interval for priority of dengue vaccination for each city. In this work we take advantage of the numerical approach to evaluate the proportion of dengue cases by age instead of founding an analytical function. As CruzPacheco et al. [14] we find out another important point that should be emphasized: the optimal age for vaccination does not match with the age of highest number of cases, showing up that this quantity is not a good criterion to prioritize the population to be firstly vaccinated.
Since the results are based on dengue reported cases in 10 Brazilian cities from 2001 to 2014, from the theoretical point of view, we could say that if a safe and effective vaccine should be applied before 2001 to the target population, it would have been possible to avoid the cases of dengue in children under the age of 13 (taking out the outliers above the maximal optimal age among the 10 cities in Fig. 4), as well as to reduce the incidence of dengue in adults, although this is only a plausible hypothesis. In addition, as the contingent of children under the age of 13 in Brazil is very large (around 42 millions), this would require great operational efforts from the National Immunization Program, as well as great financial resources.
Considering the whole country, these results indicate that the confidence age range of vaccination against dengue is from 2 to 10 years old, which implies that the target population is reduced to 27 millions of individuals which is still a large population for a vaccination campaign, mainly if it is needed more than one dose. Nowadays, it is not so far to dispose of a safe and effective vaccine [5, 6] and when available, possibly, there will be no scale production able to satisfy a demand of such magnitude. Moreover, it should be noted that more than a hundred countries have been affected by dengue epidemics [22], and most of them will also demand for vaccine. For these reasons, it is necessary to adapt vaccination strategies to the cost and quantity of the vaccine eventually available. Due to the limitations explained above, in the particular case of Brazil, a country with continental dimension with more than 200 millions of habitants, an alternative to reduce the target population should be to apply vaccine in areas that had already large epidemics as well as areas with higher risk of severe forms of the disease, because large proportions of the population already have antibodies against one or more serotypes of the dengue virus [7].
The method presented in this work provides another alternative using dengue notification data in different cities that meet the epidemiological diversity of dengue. For instance, using data of dengue epidemics during 2014, the confidence age range for vaccination should be from 2 to 10 years old in Rio de Janeiro, meanwhile in Fortaleza it is from 3 to 7 years old, that would correspond to fewer vaccine doses (see Table 1). Moreover, based on the statistical analysis, it is possible to consider the past history of the disease for each city in order to estimate the optimal age of vaccination and its variation along a time interval (see Fig. 4), setting up the priority age vaccination for each city based on the dispersion of the optimal age for the previous longterm time series of dengue cases. These results reveal heterogeneous scenarios of priority age vaccination for the analysed cities in a large country as Brazil. Finally, it is worth to notice that the decision should be subordinated not only to the age range defined by the model, but also to the age health restrictions of the vaccine.
The main advantage of using mathematical models to define the target population for vaccination is the possibility of addressing and checking different scenarios that may provide valuable information for the Health Public System. In addition, due to patterns variability of the age distribution of cases, the numerical approach is more suitable to be used with epidemiological data frequently updated. Finally it worth to say that it is a powerful method to consider more complex situations such as different vaccination schedules, taking advantage of the variation of parameters and functions with age.
Conclusion
In conclusion, the findings of this study indicate that the model and methods used here, coupling mathematical methods with actual epidemics data, may constitute a valuable tool to define the age range to be prioritized for dengue vaccination, attending to the necessity of the application of minimum doses of vaccine necessary to obtain herd immunity capable of reducing transmission of the dengue virus, and consequently to control epidemics.
Abbreviations
 BH:

Belo Horizonte
 Ca:

Campinas
 CG:

Campo Grande
 Cu:

Cuiabá
 DHF:

Dengue hemorrhagic fever
 Fo:

Fortaleza
 Go:

Goiânia
 IBGE:

Brazilian institute of geography and statistics
 Ma:

Maceió
 Na:

Natal
 ODE:

Ordinary differential equations
 PDE:

Partial differential equations
 RJ:

Rio de Janeiro
 RP:

Ribeirão Preto
 SI:

Susceptible, infected
 SINAN:

Brazilian Notified Disease Information System
 SIR:

Susceptible, infected, recovered
References
 1
Bailey NT. The Mathematical Theory of Infectious Diseases and its Applications, 255. London: Charles Griffin and Company LTD; 1975.
 2
Gubler DJ. Dengue, urbanization and globalization: the unholy trinity of the 21(st) century. Trop Med Health. 2011; 39:3–11.
 3
Esteva L, Rivas G, Yang HM. Modelling parasitism and predation of mosquitoes by water mites. J Math Biol. 2006; 53:540–55.
 4
Ferreira CP, Yang HM, Esteva L. Assesing the suitability of sterile insect technique applied to Aedes aegypti. J Biol Syst. 2008; 16:565–77.
 5
Halstead BS, Heinz FX, Barrett ADT, Roehrigd TJ. Dengue virus: molecula basis of cell entry and pathogenesis. Vaccine. 2003; 23:849–56.
 6
Wichmanna O, Vannicea K, Asturiasc EJ, Lunae EJA, Longinif I, Lopezg AL, et al. Liveattenuated tetravalent dengue vaccines: The needs and challenges of postlicensure evaluation of vaccine safety and effectiveness. Vaccine. 2017; 35:5535–42.
 7
World Health Organisation (WHO). http://www.who.int/vaccinesafety/committee, Accesed 7 Dec 2017.
 8
Billings L, Fiorillo A, Schwarts IB. Vaccinations in disease models with antibodydependent enhancement. Math Biosci. 2008; 211:265–81.
 9
Amaku M, Coudeville L, Massad E. Designing a vaccination strategy against dengue. Rev Inst Med Trop. 2012; 54:S18–S21.
 10
Chao DL, Halstead SB, Halloran ME, Longini Jr. Controlling Dengue with Vaccines in Thailand. PLoS Negl Trop Dis. 2012; 6:e1876.
 11
Coudeville L, Garnett GP. Transmission dynamics of the four dengue serotypes in Southern Vietnam and the potential impact of vaccination. PLoS ONE. 2012; 7:e51244.
 12
Mbah MLN, Durham DP, Medlock J, Galvani AP. Country and agespecific optimall allocation of dengue vaccines. J Theor Biol. 2014; 342:15–22.
 13
Maier SB, Huang X, Massad E, Amaku M, Burattini MN, Greenhalgh D. Analysis of the optimal vaccination age for dengue in Brazil with a tetravalent dengue vaccine. Math Biosci. 2017; 294:15–32.
 14
Cruz G, Esteva L, Vargas C. Vaccination Strategies for SIR VectorTransmitted Diseases. Bull Math Biol. 2014; 76:2073–90.
 15
Ministério da Saúde. SINAN. Brazilian Notified Disease Information System. Brasl̀lia: Ministério da Saúde. (Acessed 2016 Jan). Available at http://www2.datasus.gov.br/DATASUS.
 16
Teixeira MG, Costa MN, Barreto ML, Barreto FR. Epidemiologia da dengue. In: Dengue Teorias e Práticas. Rio de Janiero: Fiocruz: 2015. p. 293–315.
 17
Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013; 496:504–7.
 18
World Development Indicators. Available in: http://databank.worldbank.org/data/reports.aspx?source=worlddevelopmentindicators. Accessed 19 Nov 2018.
 19
The Lancet Infectious Diseases 18 (2018), P123 (Editorial) The dengue vaccine dilemma. https://doi.org/10.1016/S14733099(18)300239.
 20
World Health Organization (WHO). http://www.who.int/immunization/diseases/dengue/q_and_a_dengue_vaccine_Dengvaxia/en/ Accessed 18 Oct 2018.
 21
Cavalcanti LP, Vilar D, SouzaSantos R, Teixeira MG. Change in Age Pattern of Persons with Dengue, Northeastern Brazil. Emerg Infect Dis. 2011; 17:132–3.
 22
World Health Organization (WHO). 2017. http://www.who.int/denguecontrol/epidemiology/en/. Accessed 2 Feb 2018.
 23
Brazilian Institute of Geography and Statistics (IBGE). Available in: https://portaldemapas.ibge.gov.br. Accessed 20 Dec 2018.
Acknowledgments
The authors thank the two anonymous reviewers for their careful reading and comments that helped on the improvement of the manuscript. The authors also acknowledge the support for the research foundations referred in funding section.
Funding
LLC had a scholarship from Brazilian Research Foundation (CNPq) for supporting her PhD thesis about dengue vaccination strategies, CPF had grant 16/237383 from São Paulo Research Foundation (FAPESP), STRP had grant Pronex (PNX 0006/2009) and International Cooperation (INT0002/2016), both from Bahia Research Foundation (FAPESB). MGT, MCNC and STRP were supported by INCTCitecs 57 386/20089 from MCTCNPqBrazil. LE was supported by Project IN113716 PAPIITUNAM. The funding support alows the authors to perform scientific visit in which the design of the study, analysis and discussion of the results and partial writing of the manuscript were done.
Availability of data and materials
The datasets analysed during the current study are available in the Sistema de Informação de Agravo de Notificação (SINAN)  Dados do Sistema Único de Saúde (DATASUS), [http://www2.datasus.gov.br/DATASUS/index.php?area=0203].
Author information
Affiliations
Contributions
STRP, MGT, MCNC proposed the main idea of the work based on a previous work of LE. LLC, MGT and MCNC contributed to deseign and acquisition of data as well as the epidemiological discussion of the results. LLC and CPF manipulated the data, generated the results, made the figures and tables. CPF worked on the program code. LE generated the first version of the manuscript. STRP and CPF coordinated the discussion that lead to the critical revised version of the manuscript, including the cover letter with the answers to reviewer’s comments. STRP, LE and CPF are responsible for the analysis of the results, stablishing the link between the mathematical results and their epidemiological meaning. All the authors contributed on the process of writing and discussing of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The ethics committee of Instituto de Saúde Coletiva da Universidade Federal da Bahia (CEP/ISC) with reference number 1.099.503 aproved the study done. Also, the used data from DATASET (reference 15) are secondary data which are public available (registered data).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional information
Authors’ information
LLC (lucianacardim@yahoo.com.br), STRP (suani@ufba.br), MGT (magloria@ufba.br), MCNC (mcncosta@ufba.br), LE (lesteva@ciencias.unam.mx), CPF (pio@ibb.unesp.br).
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.
About this article
Cite this article
Cardim, L., Pinho, S., Teixeira, M. et al. Heterogeneities in dengue spatialtemporal transmission in Brazilian cities and its influence on the optimal age of vaccination. BMC Public Health 19, 155 (2019). https://doi.org/10.1186/s1288901964269
Received:
Accepted:
Published:
Keywords
 Mathematical model
 Age profile
 Vaccination strategies