<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD with MathML3 v1.3 20210610//EN" "https://jats.nlm.nih.gov/archiving/1.3/JATS-archivearticle1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" dtd-version="1.3" article-type="research-article" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">advances-in-tropical-medicine</journal-id>
      <journal-title-group>
        <journal-title>Advances in Tropical Medicine</journal-title>
      </journal-title-group>
      <issn publication-format="electronic">3068-4048</issn>
      <publisher>
        <publisher-name>Directive Publications</publisher-name>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.52338/aitm.2026.5949</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Research</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>Age-Dependent Distributions Of Aedes-Transmitted Infections</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <name>
            <surname>Coutinho</surname>
            <given-names>Francisco Antônio Bezerra</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author" corresp="yes">
          <name>
            <surname>Amaku</surname>
            <given-names>Marcos</given-names>
          </name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <name>
            <surname>Kallas</surname>
            <given-names>Esper Georges</given-names>
          </name>
          <xref ref-type="aff" rid="aff3">3</xref>
        </contrib>
        <contrib contrib-type="author">
          <name>
            <surname>Massad</surname>
            <given-names>Eduardo</given-names>
          </name>
          <xref ref-type="aff" rid="aff4">4</xref>
        </contrib>
        <aff id="aff1"><label>1</label>Faculdade de Medicina, Universidade de São Paulo, Brasil; Instituto de Física, Universidade de São Paulo, Brasil</aff>
        <aff id="aff2"><label>2</label>Faculdade de Medicina, Universidade de São Paulo, Brasil; Faculdade de Medicina Veterinária e Zootecnia, Universidade de São Paulo, São Paulo, Brasil</aff>
        <aff id="aff3"><label>3</label>Faculdade de Medicina, Universidade de São Paulo, Brasil; Instituto Butantan, São Paulo, Brasil</aff>
        <aff id="aff4"><label>4</label>Instituto Butantan, São Paulo, Brasil; Faculdade de Saúde Pública, Universidade de São Paulo, São Paulo, Brasil UNASP- Centro Univesitarion Adventista de Sao Paulo</aff>
      </contrib-group>
      <pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-08-19">
        <day>19</day>
        <month>08</month>
        <year>2026</year>
      </pub-date>
      <volume>18</volume>
      <issue>1</issue>
      <history>
        <date date-type="received" iso-8601-date="2026-06-17">
          <day>17</day>
          <month>06</month>
          <year>2026</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>© 2026 The Author(s). Published by Directive Publications.</copyright-statement>
        <copyright-year>2026</copyright-year>
        <copyright-holder>The Author(s)</copyright-holder>
        <license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by/4.0/">
          <ali:license_ref>https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
          <license-p>This is an open access article distributed under the terms of the CC BY 4.0 licence (https://creativecommons.org/licenses/by/4.0/).</license-p>
        </license>
      </permissions>
      <self-uri content-type="html" xlink:href="https://directivepublications.org/advances-in-tropical-medicine/age-dependent-distributions-of-aedestransmitted-infections"/>
      <self-uri content-type="pdf" xlink:href="https://directivepublications.org/pdf/1805/age-dependent-distributions-of-aedes-transmitted-infections.pdf"/>
      <abstract>
        <p>We propose a new approach to analyze the age-distribution of reported cases for vector-transmitted infections. Using officially reported cases of dengue, Zika, chikungunya, malaria and leishmaniasis for distinct geographical areas, in different periods. Data were treated in special but well- known procedure, transforming the raw data into a density age-dependent distribution and fitting a special continuous function to it. We found that the proportion of age-dependent cases with respect to the total number of cases in a given year (or any transmission season) is probably determined by the ecological interactions between vectors and hosts. The age-distribution of the proportion of cases for the three Aedes-related infections are essentially the same independently of the magnitude of the outbreak and the geographical region considered. On the other hand, for the infections transmitted by other vectors, the age-distributions of the proportion of cases are entirely different. During specific outbreaks, the ratio between the age distribution of the proportion of officially reported cases and the total number of cases for Aedes transmitted infections such as dengue, chikungunya and zika is independent of the size of the outbreak, the size of the studied population, the period when the outbreak occurs; and the geographical region considered. Our results also suggest that the age-distribution of cases is mainly due to the interaction between vectors and their hosts. Since the interaction between the vectors and male and female hosts may be different, we investigate this and, in fact, the age-distribution is slightly different.</p>
      </abstract>
      <kwd-group kwd-group-type="author">
        <kwd>Age-dependent distribution</kwd>
        <kwd>Chikungunya</kwd>
        <kwd>Dengue</kwd>
        <kwd>Leishmaniasis</kwd>
        <kwd>Malaria</kwd>
        <kwd>Zika</kwd>
      </kwd-group>
      <funding-group>
        <funding-statement>This study was partially supported by Fundacao Nutantan and CNPq.</funding-statement>
      </funding-group>
    </article-meta>
  </front>
  <body>
    <sec id="introduction">
      <label>1</label>
      <title>Introduction</title>
      <p>Vector-borne infections frequently display strong age- and time-dependent patterns in incidence, prevalence, and disease severity<sup>[<xref ref-type="bibr" rid="ref-1">1</xref>]</sup>. Understanding the age distribution of cases is central to epidemiological inference, estimation of transmission parameters, such as the average age of infection (also referred to in the literature as the average age of prime infection<sup>[<xref ref-type="bibr" rid="ref-2">2</xref>]</sup>), and more importantly, in the design of targeted control strategies (work in preparation).</p>
      <p>In this paper, we consider the following diseases: one transmitted by <italic>Anopheles</italic> mosquitoes (malaria), one transmitted by phlebotomine sandfly (visceral leishmaniasis) and three transmitted by <italic>Aedes</italic> mosquitoes in urban areas (dengue, zika and chikungunya). We focus on the age distribution of infection in the human hosts during a given outbreak. Our main hypothesis is that the age distribution of cases is primarily dependent by the contact patterns between vectors and hosts.</p>
      <p>Because <italic>Aedes</italic> mosquitoes transmit dengue, Zika and chikungunya, if the above hypothesis is correct, one would expect that the age distribution of cases for these three infection to be similar. We demonstrate that this pattern is indeed observed empirically (see <bold>Figure (<xref ref-type="fig" rid="figure-10">10</xref>)</bold> below). In contrast, infections transmitted by different vectors exhibit markedly different age distributions, which is also supported by empirical observations (see <bold>Figure (<xref ref-type="fig" rid="figure-11">11</xref>)</bold> below).</p>
      <p>Daily activity biting habits of vectors, occupational behavior of hosts, and housing conditions<sup>[<xref ref-type="bibr" rid="ref-3">3</xref>,<xref ref-type="bibr" rid="ref-4">4</xref>]</sup>, influence human-vector contact. Differences in vector behavior—such as indoor versus outdoor biting—may further shape the age distribution of cases5. In this paper, however, we analyze the age distribution of Aedes-transmitted infections in large and medium-size urban centers and for Brazil as a whole. In Brazilian urban centers the transmission is mainly peri-domestic<sup>[<xref ref-type="bibr" rid="ref-6">6</xref>,<xref ref-type="bibr" rid="ref-7">7</xref>]</sup>. Furthermore, in these settings, the infection ‘travels’ throughout the city carried by infected hosts<sup>[<xref ref-type="bibr" rid="ref-8">8</xref>]</sup>. This is clearly illustrated by dengue transmission in São Paulo city<sup>[<xref ref-type="bibr" rid="ref-8">8</xref>]</sup> (approximately 12 million inhabitants distributed among 1,700 neighborhoods).</p>
      <p>For infections transmitted by insect vectors in hyperendemic areas, cases tend to cluster in children and adolescents because of early exposure and the accumulation of serotype-specific immunity<sup>[<xref ref-type="bibr" rid="ref-9">9</xref>]</sup>. However, in settings where transmission is relatively low, as in the locations considered in this study, the burden of symptomatic disease may occur primarily among adults because a large fraction of the adult population remains immunologically naïve<sup>[<xref ref-type="bibr" rid="ref-10">10</xref>]</sup>. Thus, age-dependent exposure to vectors is a major determinant of the age distribution of cases.</p>
      <p>In this paper, we hypothesize that, in Brazil at least, during specific outbreaks, the ratio between the age distribution of the proportion of officially reported cases, the total number of cases for Aedes-transmitted infections such as dengue, chikungunya, and zika is <bold>independent</bold> of:</p>
      <list list-type="bullet">
        <list-item>
          <p>The size of the outbreak;</p>
        </list-item>
        <list-item>
          <p>The size of the studied population;</p>
        </list-item>
        <list-item>
          <p>The period when the outbreak occurs; and</p>
        </list-item>
        <list-item>
          <p>The geographical region considered.</p>
        </list-item>
      </list>
      <p>The above hypotheses are supported by a proper analysis of empirical data as demonstrated in this paper.</p>
      <p>In summary, the age-dependent distribution of vector-borne infections emerges from a combination of transmission intensity, immunity acquisition, behavioral exposure, and demographic structure. Analyzing age patterns provides valuable insights into epidemiological processes and supports the design of targeted interventions, such as vaccination programs or vector control strategies focused on the most affected age groups. These issues will be addressed in a future work.</p>
      <sec id="testing-the-hypothesis">
        <label>1.1</label>
        <title>Testing the hypothesis</title>
        <sec id="the-available-data">
          <title>The available data</title>
          <p>To test our hypothesis empirical evidence is required. We used age-dependent incidence of dengue for different geographical regions and years to illustrate our method. The data were obtained from the National System of Compulsory Notification of Infections of the Brazilian Ministry of Health (SINAN)<sup>[<xref ref-type="bibr" rid="ref-11">11</xref>]</sup>, which provides the dengue incidences stratified by ages for each year from 2001 to 2026.</p>
          <p>As demonstrated by <bold>Figures (<xref ref-type="fig" rid="figure-5">5</xref>), (<xref ref-type="fig" rid="figure-7">7</xref>)</bold> and <bold>(9)</bold>, which show the outbreaks of dengue in Brazil as a whole, and in cities of Rio de Janeiro and Sao Jose do Rio Preto, respectively, dengue is far from a steady state. The number of cases stratified by age is provided in SINAN<sup>[<xref ref-type="bibr" rid="ref-11">11</xref>]</sup> and is illustrated in <bold><xref ref-type="table" rid="table-1">Table 1</xref></bold> for the city Sao Jose do Rio Preto, State of Sao Paulo, South Eastern Brazil, in 2019.</p>
          <p>To calculate the age density of cases, we introduced the binning procedure for density estimation shown in Figure <xref ref-type="fig" rid="figure-1">1</xref> (see<sup>[<xref ref-type="bibr" rid="ref-12">12</xref>]</sup> and Appendix I for details). For this purpose, we considered the midpoints of each age interval<sup>[<xref ref-type="bibr" rid="ref-13">13</xref>]</sup>. The number of cases in each interval was divided by the width of the age interval (in years).</p>
          <p>For example, in the age interval 20-39 years, the number of reported cases shown in <xref ref-type="table" rid="table-1">table 1</xref>, was 12,213. This corresponds to a density of 12,213 divided by 20, that is, 610.65. As mentioned above, reader unfamiliar with these procedures, are referred to Appendix I for detailed explanation.</p>
          <table-wrap id="table-1">
            <label>Table 1</label>
            <caption>
              <p>Age Dependence of the Number of Cases in a Dengue Outbreak in Sao Jose do Rio Preto, 2019.</p>
            </caption>
            <table>
              <thead>
                <tr>
                  <th scope="col"><bold>Age</bold></th>
                  <th scope="col"><bold>Number of Cases</bold></th>
                </tr>
              </thead>
              <tbody>
                <tr>
                  <th scope="row">&lt;1</th>
                  <td>186</td>
                </tr>
                <tr>
                  <th scope="row">1-4</th>
                  <td>568</td>
                </tr>
                <tr>
                  <th scope="row">5-9</th>
                  <td>1573</td>
                </tr>
                <tr>
                  <th scope="row">10-14</th>
                  <td>2321</td>
                </tr>
                <tr>
                  <th scope="row">15-19</th>
                  <td>2701</td>
                </tr>
                <tr>
                  <th scope="row">20-39</th>
                  <td>12213</td>
                </tr>
                <tr>
                  <th scope="row">40-59</th>
                  <td>9359</td>
                </tr>
                <tr>
                  <th scope="row">60-64</th>
                  <td>1622</td>
                </tr>
                <tr>
                  <th scope="row">65-69</th>
                  <td>1106</td>
                </tr>
                <tr>
                  <th scope="row">70-79</th>
                  <td>1163</td>
                </tr>
                <tr>
                  <th scope="row">&gt;80</th>
                  <td>422</td>
                </tr>
                <tr>
                  <th scope="row">Total</th>
                  <td>33237</td>
                </tr>
              </tbody>
            </table>
          </table-wrap>
          <p>We fitted a density function, a Gaussian Mesa Function<sup>[<xref ref-type="bibr" rid="ref-14">14</xref>]</sup>, to the data so that we obtain the yearly density function as shown in Figure <xref ref-type="fig" rid="figure-2">2</xref> (see Appendix II for details).</p>
          <fig id="figure-1">
            <label>Figure 1</label>
            <caption>
              <p>Age density of dengue cases for the 2019 outbreak in Sao Jose do Rio Preto.</p>
            </caption>
            <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-1.png" mimetype="image" mime-subtype="png"/>
          </fig>
          <p>We used the same procedure for all years outbreaks from 2004-2024 for the three geographical regions analyzed, that is, Brazil as a whole (population 203,080,756), Rio de Janeiro (population 6,211,223) and Sao Jose do Rio Preto (population 480,393) (Figure <xref ref-type="fig" rid="figure-3">3</xref>).</p>
          <fig id="figure-2">
            <label>Figure 2</label>
            <caption>
              <p>Fitted age-dependent density function to the data as in Figure <xref ref-type="fig" rid="figure-1">1</xref>.</p>
            </caption>
            <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-2.png" mimetype="image" mime-subtype="png"/>
          </fig>
          <p>Next, we calculated the <italic>age density proportion of cases</italic>, dividing each point of the above curve by the total number of cases. The results are shown below in Figures (<xref ref-type="fig" rid="figure-4">4</xref>), (<xref ref-type="fig" rid="figure-6">6</xref>) and (<xref ref-type="fig" rid="figure-8">8</xref>), for each site for all years.</p>
          <p>We also applied the above procedures to Chikungunya Fever, Zika Fever, Malaria and Visceral Leishmaniosis for Brazil as a whole in the year 2017 (see Figure <xref ref-type="fig" rid="figure-11">11</xref>). The same method was applied to Chikungunya and Zika for other cities and other periods and the results are the same.</p>
          <fig id="figure-3">
            <label>Figure 3</label>
            <caption>
              <p>Map of the locations studied.</p>
            </caption>
            <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-3.png" mimetype="image" mime-subtype="png"/>
          </fig>
        </sec>
        <sec id="preliminary-mathematical-considerations">
          <title>Preliminary mathematical considerations</title>
          <p>Let us denote <italic>X(a,Y)</italic> the double density of new cases of individuals with ages between <italic>a</italic> and <italic>a+da</italic> occurring between the times <italic>Y</italic> and <italic>Y+dY</italic>. The total number of cases at a specific year <italic>Y</italic>, summed over all ages <italic>Λ(Y)</italic>, is given by:</p>
          <disp-formula id="eq-1"><label>(1)</label><graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/eq-1.png" mimetype="image" mime-subtype="png"><alt-text>Λ(Y)=∫0∞ X(a,Y)da</alt-text></graphic></disp-formula>
          <p>The right-hand side of equation <xref ref-type="disp-formula" rid="eq-1">(1)</xref> is the sum of incident cases occurring in a specific year Y across all ages.</p>
          <p>Rearranging equation <xref ref-type="disp-formula" rid="eq-1">(1)</xref>, we obtain</p>
          <disp-formula id="eq-2"><label>(2)</label><graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/eq-2.png" mimetype="image" mime-subtype="png"/></disp-formula>
          <p>The fraction <inline-formula><inline-graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/ieq-51-0.png" mimetype="image" mime-subtype="png"/></inline-formula> in equation <xref ref-type="disp-formula" rid="eq-2">(2)</xref> represents the proportion of cases at time Y, occurring in individuals aged between a and a+da. We observed the previously overlooked fact, that this quantity is approximately a function of age a only, independently of both location and the year (see <bold>Figures (<xref ref-type="fig" rid="figure-4">4</xref>), (<xref ref-type="fig" rid="figure-6">6</xref>)</bold> and <bold>(8)</bold>):</p>
          <disp-formula id="eq-3"><label>(3)</label><graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/eq-3.png" mimetype="image" mime-subtype="png"/></disp-formula>
          <p>It is remarkable that, for the sites and the years analyzed, this function is independent of both the magnitude of the outbreak and its geographical location, whether within a municipality or across distant regions of the country (see <bold>Figure (<xref ref-type="fig" rid="figure-4">4</xref>), (<xref ref-type="fig" rid="figure-6">6</xref>)</bold> and <bold>(8)</bold> below). <bold>Figure <xref ref-type="fig" rid="figure-10">10</xref></bold> shows the average distribution of the proportion cases across the three regions analyzed, demonstrating their similarity. This observation supports our hypothesis that the distribution depends primarily on the age-dependent exposure of humans to mosquitoes’ bites. Furthermore, as shown below, the age distribution of proportion of cases is similar for infections transmitted by aedes mosquitoes and is very different for other vectors (see <bold>Figure <xref ref-type="fig" rid="figure-11">11</xref></bold>). For dengue, at least, we show below that there is a small difference between the distribution in the proportion of cases between men and women, with the latter showing a slightly larger average age of acquiring the infection. This is to be expected because the behavior of each sex with respect to the exposure to the mosquitoes’ bites is probably different (<bold>Figure (<xref ref-type="fig" rid="figure-12">12</xref>)</bold> to <bold>(14)</bold>).</p>
          <p>Finally, the average age of prime infection, <inline-formula><inline-graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/ieq-55-0.png" mimetype="image" mime-subtype="png"/></inline-formula> is given by:</p>
          <disp-formula id="eq-4"><label>(4)</label><graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/eq-4.png" mimetype="image" mime-subtype="png"/></disp-formula>
          <p>where L is the average lifespan of the population (sometimes referred to as ‘life expectancy’).</p>
        </sec>
      </sec>
    </sec>
    <sec id="results">
      <label>2</label>
      <title>Results</title>
      <p>As mentioned above, <bold>Figure <xref ref-type="fig" rid="figure-4">4</xref></bold> shows the age distribution of the proportion of dengue cases in each year in Brazil for the outbreaks occurring from 2014 to 2024.</p>
      <p>Note that there is a small variation in the distribution of cases with age, although the outbreaks varied considerably from year to year, as shown in Figure <xref ref-type="fig" rid="figure-5">5</xref>. The average age (see below) at which individuals acquire dengue in Brazil, for the period analyzed, follows from the average of the curves in Figure <xref ref-type="fig" rid="figure-4">4</xref> and results in 34 years.</p>
      <fig id="figure-4">
        <label>Figure 4</label>
        <caption>
          <p>Age distribution of the proportion of dengue cases for the 2014-2024 outbreaks in Brazil as a whole. Red line represents the average.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-4.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>Figure <xref ref-type="fig" rid="figure-6">6</xref> shows the age distribution of the proportion of dengue cases in each year in Rio de Janeiro city (RJ) for the outbreaks occurring from 2014 to 2024.</p>
      <fig id="figure-5">
        <label>Figure 5</label>
        <caption>
          <p>Dengue outbreaks of dengue in Brazil as a whole for the period 2014-2014. Small arrows indicate the peak of each yearly outbreak.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-5.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>Note again, that there is a small variation in the age distribution of cases, independently of the vastly different sizes of the outbreaks in that period, as shown in Figure <xref ref-type="fig" rid="figure-7">7</xref>.</p>
      <fig id="figure-6">
        <label>Figure 6</label>
        <caption>
          <p>Age distribution of the proportion of dengue cases for the 2014-2024 outbreaks in Rio de Janeiro city. Red line represents the average.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-6.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>Figure <xref ref-type="fig" rid="figure-8">8</xref> shows the age distribution of the proportion of dengue cases in each year in Sao Jose do Rio Preto city (SJRP) for the outbreaks occurring from 2014 to 2024.</p>
      <fig id="figure-7">
        <label>Figure 7</label>
        <caption>
          <p>Dengue outbreaks of dengue in Rio de Janeiro city for the period 2014-2024. Small arrows indicate the peak of each yearly outbreak.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-7.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>There is a slightly larger variation in the age distribution of proportion of cases (probably due to stochastic variations due to the smaller size of this site) in spite of the great variability in the size of the outbreaks, as shown in Figure <xref ref-type="fig" rid="figure-9">9</xref>.</p>
      <fig id="figure-8">
        <label>Figure 8</label>
        <caption>
          <p>Age distribution of the proportion of dengue cases for the 2014-2024 outbreaks in Sao Jose do Rio Preto city. Red line represents the average.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-8.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>The outbreaks in this city are slightly out of phase with respect to the Figures (<xref ref-type="fig" rid="figure-5">5</xref>) and (<xref ref-type="fig" rid="figure-7">7</xref>). This is probably due to different ecological interactions between vectors and hosts in this city.</p>
      <p>Figure <xref ref-type="fig" rid="figure-10">10</xref> shows the average distribution over the period analyzed of proportion of cases for Brazil (Br), Rio de Janeiro (RJ) and Sao Jose do Rio Preto (SJRP).</p>
      <fig id="figure-9">
        <label>Figure 9</label>
        <caption>
          <p>Dengue outbreaks of dengue in Sao Jose do Rio Preto city for the period 2014-2024. Small arrows indicate the peak of each yearly outbreak.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-9.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>It is remarkable that the age distribution of cases are similar for the three geographical areas studied.</p>
      <p>In Figure <xref ref-type="fig" rid="figure-11">11</xref> we show the age distribution of the proportion of cases for Dengue, Chikungunya Fever, Zika Fever, Malaria and Visceral Leishmaniosis for Brazil as a whole in the year 2017.</p>
      <fig id="figure-10">
        <label>Figure 10</label>
        <caption>
          <p>Average age distribution of the proportion of dengue cases in Brazil as a whole, Rio de Janeiro city and Sao Jose do Rio Preto city.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-10.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>Note the similarity between the three <italic>Aedes aegypti</italic> transmitted infections, which differ greatly from malaria and visceral leishmaniosis because the two latter infections are transmitted by different vectors<sup>[<xref ref-type="bibr" rid="ref-15">15</xref>,<xref ref-type="bibr" rid="ref-16">16</xref>]</sup> with different human age exposition to each of them. We are aware that urban yellow fever is also transmitted by <italic>Aedes aegypti</italic> in Africa<sup>[<xref ref-type="bibr" rid="ref-17">17</xref>]</sup> but in Brazil, it is only transmitted by sylvatic vectors<sup>[<xref ref-type="bibr" rid="ref-18">18</xref>]</sup> and therefore will not be analyzed in this paper.</p>
      <p>Next, we stratified the proportion of cases of dengue in 2024 by sex for the three geographical regions considered, and used the 2024 outbreak as an illustration.</p>
      <p>Figure <xref ref-type="fig" rid="figure-12">12</xref> shows the sex-related results for Brazil as a whole. The number (size of the outbreak) of reported cases in 2024 were 2,912,187 and 3,504,989 dengue cases for men and women, respectively.</p>
      <fig id="figure-11">
        <label>Figure 11</label>
        <caption>
          <p>Distribution of age-dependent proportion of dengue, chikungunya, zika fever, malaria and visceral leishmaniosis cases in Brazil as a whole.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-11.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>The average age of infection was 33 years and 37 years for men and women, respectively.</p>
      <p>Figure <xref ref-type="fig" rid="figure-13">13</xref> shows the results for the city of Rio de Janeiro. The number of reported cases of dengue (size of the outbreak) were 48,116 and 61,822 cases for men and women, respectively.</p>
      <fig id="figure-12">
        <label>Figure 12</label>
        <caption>
          <p>Distribution of age-dependent proportion of dengue cases for men and women in Brazil as a whole.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-12.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>The average age of infection was 29 years and 35 years for men and women, respectively.</p>
      <p>Figure <xref ref-type="fig" rid="figure-14">14</xref> shows the results for the city of Sao Jose do Rio Preto. The number of reported dengue cases were (size of the outbreak) 16,612 and 19,703 cases for men and women, respectively.</p>
      <fig id="figure-13">
        <label>Figure 13</label>
        <caption>
          <p>Distribution of age-dependent proportion of dengue cases for men and women in Rio de Janeiro city.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-13.png" mimetype="image" mime-subtype="png"/>
      </fig>
      <p>The average age of infection was 35 years and 38 years for men and women, respectively.</p>
      <p>Note that the differences between sexes are similar in the three populations analyzed.</p>
      <p>We think that the small differences observed between sexes do not justify separate control strategies for each sex to be considered in a future work.</p>
      <fig id="figure-14">
        <label>Figure 14</label>
        <caption>
          <p>Distribution of age-dependent proportion of dengue cases for men and women in Sao Jose do Rio Preto city.</p>
        </caption>
        <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-14.png" mimetype="image" mime-subtype="png"/>
      </fig>
    </sec>
    <sec id="discussion">
      <label>3</label>
      <title>Discussion</title>
      <p>In this paper, we propose a new hypothesis related to the empirical observation that the age distribution of the proportion of reported cases for infections transmitted by Aedes mosquitoes are essentially independent of transmission intensity and geographical area. In fact, the age distribution of the proportion of cases is remarkably similar for the three main infections transmitted by Aedes aegypti - dengue, chikungunya and Zika, in Brazilian urban centers.</p>
      <p>Because the distribution depends mainly on the exposure of hosts to mosquito’s bites, it differs for infections transmitted by other vectors (Figure <xref ref-type="fig" rid="figure-11">11</xref>) and may differ in other countries where this exposure patterns may vary.</p>
      <p>As mentioned above, using the officially reported number of cases for three geographical regions, we found that the age distribution of cases does not differ significantly across sites or across years. This finding may be used in designing age-target interventions and optimizing disease-control strategies, which will be explored in a future work.</p>
      <p>As shown in the literature<sup>[<xref ref-type="bibr" rid="ref-8">8</xref>]</sup>, outbreaks in cities occur in different neighborhoods and travels along to other areas of the city. Therefore, our results for a given city assumes that each outbreak occurs in a region where there is a considerable fraction of susceptible individuals. Note also that dengue has at least four different serotypes, differences that were not considered in our analysis.</p>
      <p>Some important limitations of our approach are worth mentioning. The first is that by basing the calculations on the officially reported number of cases, the model considers neither underreporting nor recovering from infection. However, underreporting is important only if it modifies the function ∅ (a) (see equation <xref ref-type="disp-formula" rid="eq-3">(3)</xref>), for instance if underreporting occurs predominantly in a certain age interval. This study is valid only for reported cases. However, to use our results to design an intervention, underreporting is important only because it is not necessary to vaccinate everybody in a given age interval since some of them may be already immune by previous infections. The result is that the intervention may be overestimated. Recovering from infection does not modify the model’s performance since it would not change the function ∅ (a).</p>
      <p>Also important to mention is that our results do not allow the proposal of interventions intended to limit the size of the vector population<sup>[<xref ref-type="bibr" rid="ref-19">19</xref>,<xref ref-type="bibr" rid="ref-20">20</xref>]</sup>.</p>
      <p>Another limitation of our approach is that, for dengue, there are 4 serotypes. In an ideal scenario, all the above calculations should be done stratifying the cases by serotype. However, this limitation, although important, is valid for dengue only and the same theoretical approach is valid for other Aedes-mediated vector-borne infections. We think that the similarity of the age distribution of cases between the three infections justify not distinguishing serotypes for the case of dengue.</p>
      <p>The most important finding of this paper is that the age dependent distribution of the proportion of cases does not vary, neither with the intensity of the transmission nor with the geographical area in Brazil. This function may be completely different in other countries, depending on the age distribution for which exposition to the mosquitoes’ bites occurs, because this may modify the form of ∅ (a) and may even modify the shape and size of the outbreak (see Equation <xref ref-type="disp-formula" rid="eq-3">(3)</xref>). For instance, contrasting with the average age of dengue acquisition of 34 years in Brazil, in Thailand, this average is lower, about 24 years<sup>[<xref ref-type="bibr" rid="ref-21">21</xref>]</sup>.</p>
      <p>Finally, our approach completely ignores stochasticity. One consequence of its deterministic nature is that ∅ (a) remains essentially unmodified, even if the number of cases is so small that it does not allow infection in all ages. For instance, suppose that in a certain outbreak we have less than eighty cases in total. Then, some ages intervals will have no cases but our model predicts fractions of cases for all ages.</p>
    </sec>
    <sec id="conclusions">
      <label>4</label>
      <title>Conclusions</title>
      <p>Our results allow us to conclude that our main hypothesis, namely that the ratio between the age distribution of the proportion of officially reported cases and the total number of cases for Aedes-transmitted infections is independent of: the size of the outbreak; the size of the studied population; the period of the outbreak occurs; and the geographical region considered, was confirmed.</p>
      <p>Our results also suggest that the age-distribution of cases is mainly due to the interaction between vectors and their hosts. Since the interaction between the vectors and male and female hosts may be different, we investigate this and, in fact, the age-distribution is slightly different.</p>
      <p>This is actually misleading and quite meaningless. The graph seems to indicate that, for instance, at age 20 we have about 5,000 cases, at age 25 we have about 7,000 cases, at age 30 we have about 12,000 cases and at age 35 we have about 11,000. Only considering 4 ages, we have about 34,000 cases, which exceeds the total number of observed cases.</p>
      <p>However, if the ordinate (the y-axis) is chosen as density of cases rather than number of cases, this inconsistency disappears. By density, we mean the number of cases in each age interval divided by the length of the interval. The age-dependent density of cases is shown in <xref ref-type="table" rid="table-a1">Table A1</xref>.</p>
    </sec>
  </body>
  <back>
    <app-group>
      <app id="appendix-a">
        <label>A</label>
        <title>Appendix A</title>
        <p>In this appendix, we explain that treatment of the original data from SINAN<sup>[<xref ref-type="bibr" rid="ref-11">11</xref>]</sup>. The binning of SINAN<sup>[<xref ref-type="bibr" rid="ref-11">11</xref>]</sup> is reproduced in <xref ref-type="table" rid="table-1">Table 1</xref> of the main text, which has unequal age intervals.</p>
        <p>Our problem is how fit a continuous curve to these data. A first, and quite natural but wrong way, of doing it is to plot the raw data from <xref ref-type="table" rid="table-1">Table 1</xref>, with age as abscissa (the x- axis) and the number of cases of each age interval in the middle of the interval, as shown in Figure <xref ref-type="fig" rid="figure-a1">A1</xref>.</p>
        <fig id="figure-a1">
          <label>Figure A1</label>
          <caption>
            <p>Plot of the original number of cases against age.</p>
          </caption>
          <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-a1.png" mimetype="image" mime-subtype="png"/>
        </fig>
        <p>This is actually misleading and quite meaningless. The graph seems to indicate that, for instance, at age 20 we have about 5,000 cases, at age 25 we have about 7,000 cases, at age 30 we have about 12,000 cases and at age 35 we have about 11,000. Only considering 4 ages, we have about 34,000 cases, which exceeds the total number of observed cases.</p>
        <p>However, if the ordinate (the y-axis) is chosen as density of cases rather than number of cases, this inconsistency disappears. By density, we mean the number of cases in each age interval divided by the length of the interval. The age-dependent density of cases is shown in <xref ref-type="table" rid="table-a1">Table A1</xref>.</p>
        <p>Figure <xref ref-type="fig" rid="figure-a2">A2</xref> (Figure <xref ref-type="fig" rid="figure-1">1</xref> of the main text) shows the density plot. The ordinate now is the derivative of the number of cases with respect to ages.</p>
        <table-wrap id="table-a1">
          <label>Table A1</label>
          <caption>
            <p>Age Dependence of the Density of Cases in a Dengue Outbreak in Sao Jose do Rio Preto, 2019.</p>
          </caption>
          <table>
            <thead>
              <tr>
                <th scope="col"><bold>Age</bold></th>
                <th scope="col"><bold>Density of Cases</bold></th>
              </tr>
            </thead>
            <tbody>
              <tr>
                <th scope="row">&lt;1</th>
                <td>93</td>
              </tr>
              <tr>
                <th scope="row">1-4</th>
                <td>142</td>
              </tr>
              <tr>
                <th scope="row">5-9</th>
                <td>314.6</td>
              </tr>
              <tr>
                <th scope="row">10-14</th>
                <td>464.2</td>
              </tr>
              <tr>
                <th scope="row">15-19</th>
                <td>540.2</td>
              </tr>
              <tr>
                <th scope="row">20-39</th>
                <td>610.65</td>
              </tr>
              <tr>
                <th scope="row">40-59</th>
                <td>467.95</td>
              </tr>
              <tr>
                <th scope="row">60-64</th>
                <td>324.4</td>
              </tr>
              <tr>
                <th scope="row">65-69</th>
                <td>221.2</td>
              </tr>
              <tr>
                <th scope="row">70-79</th>
                <td>116.3</td>
              </tr>
              <tr>
                <th scope="row">&gt;80</th>
                <td>42.2</td>
              </tr>
            </tbody>
          </table>
        </table-wrap>
        <fig id="figure-a2">
          <label>Figure A2</label>
          <caption>
            <p>Plot of density of cases with respect to ages.</p>
          </caption>
          <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-a2.png" mimetype="image" mime-subtype="png"/>
        </fig>
        <p>Each point in this graph, say at age a, represents the number of cases in a vanishingly small interval between a and a+da, so that the total number of cases is obtained by taking the integral of the curve, that is, the area below the curve. This can be checked by integrating the Gaussian Mesa Function<sup>[<xref ref-type="bibr" rid="ref-14">14</xref>]</sup> that fits the density (see Appendix II below), as in Figure <xref ref-type="fig" rid="figure-2">2</xref> of the main text.</p>
      </app>
      <app id="appendix-b">
        <label>B</label>
        <title>Appendix B</title>
        <p>In the main text, and in the Appendix A, we mentioned that we used a Gaussian Mesa Function<sup>[<xref ref-type="bibr" rid="ref-14">14</xref>]</sup> to fit the density of cases and the distribution of the proportion of cases with respect to age.</p>
        <p>In this appendix, we briefly explain what is this function.</p>
        <p>The Gaussian mesa function (also called a flat-top Gaussian or mesa beam profile) is a modification of the classical Gaussian function designed to produce a central plateau with smooth Gaussian edges. Such functions are widely used in optics, laser physics, and beam-shaping applications where a nearly uniform intensity distribution is desired while preserving smooth, differentiable boundaries.</p>
        <p>The shape and the analytic form of this function is given below (Figure <xref ref-type="fig" rid="figure-b1">B1</xref>).</p>
        <fig id="figure-b1">
          <label>Figure B1</label>
          <caption>
            <p>Shape and analytic form of the Gaussian mesa function.</p>
          </caption>
          <graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/fig-figure-b1.png" mimetype="image" mime-subtype="png"/>
        </fig>
        <p>As shown in the figure, parameter A represents the maximum density, which occurs in the interval</p>
        <p>μ- <inline-formula><inline-graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/ieq-132-0.png" mimetype="image" mime-subtype="png"/></inline-formula> ≤t≤μ+ <inline-formula><inline-graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/ieq-132-3.png" mimetype="image" mime-subtype="png"/></inline-formula> , σ<sub>1</sub> is half the width at <inline-formula><inline-graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/ieq-132-1.png" mimetype="image" mime-subtype="png"/></inline-formula> of Gaussian on the left, σ<sub>2</sub> is half the width at <inline-formula><inline-graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/ieq-132-2.png" mimetype="image" mime-subtype="png"/></inline-formula> of Gaussian on the right and μ is the abscissa at half the plateau.</p>
        <p>As an example, the parameters that best fit the data in Figure <xref ref-type="fig" rid="figure-2">2</xref> of the main text are: A=630 cases/year; σ_<sub>L</sub>=0.042 years; σ_<sub>1</sub>=14.9 years; σ_<sub>2</sub>=28.6 years and μ=26 years.</p>
        <p>The reader can check that</p>
        <disp-formula id="eq-u1"><graphic xlink:href="https://directivepublications.org/storage/ftv2/1805-b9b54401376839735a3d/eq-u1.png" mimetype="image" mime-subtype="png"/></disp-formula>
        <p>where, as mentioned in the main text (Equation <xref ref-type="disp-formula" rid="eq-4">(4)</xref>), L is the average lifespan of the population (sometimes referred to as ‘life expectancy’).</p>
      </app>
    </app-group>
    <sec id="decl-ethics" sec-type="ethics-statement">
      <title>Ethics statement</title>
      <p>This manuscript is a hypothesis-based work and does not include any experiments involving human participants, animals, or identifiable human data. Hence, ethical approval and informed consent were not applicable.</p>
    </sec>
    <sec id="decl-funding" sec-type="funding-information">
      <title>Funding</title>
      <p>This study was partially supported by Fundacao Nutantan and CNPq.</p>
    </sec>
    <sec id="decl-contributions" sec-type="author-contributions">
      <title>Author contributions</title>
      <p>All the authors contributed equally for this work.</p>
    </sec>
    <fn-group>
      <title>Competing interests</title>
      <fn id="decl-coi" fn-type="coi-statement">
        <p>The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.</p>
      </fn>
    </fn-group>
    <ref-list>
      <title>References</title>
      <ref id="ref-1">
        <label>1</label>
        <mixed-citation>V. Botosso, A.R. Precioso, A. Wilder-Smith, D. Oliveira, F. Oliveira, C. Oliveira, C. Soares, L. Oliveira, R. Santos, C. Utescher, F. Coutinho, E. Massad, Seroprevalence of Zika in Brazil stratified by age and geographic distribution, Epidemiol. Infect. 151 (2023) 1–16. https://doi.org/<pub-id pub-id-type="doi">10.1017/S0950268823001814</pub-id>.</mixed-citation>
      </ref>
      <ref id="ref-2">
        <label>2</label>
        <mixed-citation>E. Massad, M.N. Burattini, R.S. Azevedo Neto, H.M. Yang, F.A.B. Coutinho, D.M.T. Zanetta, A model-based design of a vaccination strategy against rubella in a non-immunized community of São Paulo, Brazil, Epidemiol. Infect. 112 (1994) 579–594.</mixed-citation>
      </ref>
      <ref id="ref-3">
        <label>3</label>
        <mixed-citation>C. Braga, C.F. Luna, C.M.T. Martelli, W. de Souza, M.T. Cordeiro, N. Alexander, et al., Seroprevalence and risk factors for dengue infection in socio-economically distinct areas of Recife, Brazil, Acta Trop. 113 (3) (2010) 234–240.</mixed-citation>
      </ref>
      <ref id="ref-4">
        <label>4</label>
        <mixed-citation>M. Amaku, L. Coudeville, E. Massad, Designing a vaccination strategy against dengue, Rev. Inst. Med. Trop. Sao Paulo 54 (Suppl 18) (2012) S18–21.</mixed-citation>
      </ref>
      <ref id="ref-5">
        <label>5</label>
        <mixed-citation>R.C. Reiner, T.A. Perkins, C.M. Barker, et al., A systematic review of mathematical models of mosquito-borne pathogen transmission, J. R. Soc. Interface 10 (81) (2013) 20120921.</mixed-citation>
      </ref>
      <ref id="ref-6">
        <label>6</label>
        <mixed-citation>M. Guzman, G. Kouri, Dengue: an update, Lancet Infect. Dis. 2 (11) (2002) 33–42.</mixed-citation>
      </ref>
      <ref id="ref-7">
        <label>7</label>
        <mixed-citation>A.B.I. Cecílio, E.S. Campanelli, K.P.R. Souza, I.B. Figueiredo, M.C. Resende, Natural vertical transmission by Stegomyia albopicta as dengue vector in Brazil, Braz. J. Biol. 69 (1) (2009). https://doi.org/<pub-id pub-id-type="doi">10.1590/S1519-69842009000100015</pub-id>.</mixed-citation>
      </ref>
      <ref id="ref-8">
        <label>8</label>
        <mixed-citation>M. Amaku, F. Azevedo, M.N. Burattini, G.E. Coelho, F.A.B. Coutinho, D. Greenhalgh, L.F. Lopez, R.S. Motitsuki, A. Wilder-Smith, E. Massad, Magnitude and frequency variations of vector-borne infection outbreaks using the Ross-Macdonald model: explaining and predicting outbreaks of dengue fever, Epidemiol. Infect. 144 (16) (2016) 3435–3450.</mixed-citation>
      </ref>
      <ref id="ref-9">
        <label>9</label>
        <mixed-citation>M.G. Guzman, S.B. Halstead, H. Artsob, et al., Dengue: a continuing global threat, Nat. Rev. Microbiol. 8 (12) (2010) S7–S16.</mixed-citation>
      </ref>
      <ref id="ref-10">
        <label>10</label>
        <mixed-citation>S.B. Halstead, Dengue, Lancet 370 (9599) (2007) 1644–1652.</mixed-citation>
      </ref>
      <ref id="ref-11">
        <label>11</label>
        <mixed-citation>SINAN-MOH, Sistema Nacional de Agravos de Notificação, Ministério da Saúde do Brasil, 2025. <ext-link ext-link-type="uri" xlink:href="https://datasus.saude.gov.br/acesso-a-informacao/doencas-e-agravos-de-notificacao-de-2007-em-diante-sinan/">https://datasus.saude.gov.br/acesso-a-informacao/doencas-e-agravos-de-notificacao-de-2007-em-diante-sinan/</ext-link>.</mixed-citation>
      </ref>
      <ref id="ref-12">
        <label>12</label>
        <mixed-citation>J. Millman, S. Seely, Electronics, McGraw-Hill Book Company, New York, 1951.</mixed-citation>
      </ref>
      <ref id="ref-13">
        <label>13</label>
        <mixed-citation>Cole Notes, Statistics &amp; Data, Coles Publishing, Toronto, 2009.</mixed-citation>
      </ref>
      <ref id="ref-14">
        <label>14</label>
        <mixed-citation>F.A.B. Coutinho, M. Amaku, F.C. Boulos, J.A. de Sousa Moreira, E.N.C. de Barros, E.G. Kallas, E. Massad, Estimating the Size of the Aedes Mosquitoes' Population Involved in Outbreaks of Dengue and Chikungunya Using a Mathematical Model, Bull. Math. Biol. 87 (9) (2025) 124.</mixed-citation>
      </ref>
      <ref id="ref-15">
        <label>15</label>
        <mixed-citation>T.M.P. Oliveira, G.Z. Laporta, E.S. Bergo, L.S.M. Chaves, J.L.F. Antunes, A.S. Bickersmith, J.E. Conn, E. Massad, M.A.M. Sallum, Vector role and human biting activity of Anophelinae mosquitoes in different landscapes in the Brazilian Amazon, Parasit. Vectors 14 (1) (2021) 236. https://doi.org/<pub-id pub-id-type="doi">10.1186/s13071-021-04725-2</pub-id>.</mixed-citation>
      </ref>
      <ref id="ref-16">
        <label>16</label>
        <mixed-citation>H.J. Shimozako, J. Wu, E. Massad, Mathematical modelling for Zoonotic Visceral Leishmaniasis dynamics: A new analysis considering updated parameters and notified human Brazilian data, Infect. Dis. Model. 2 (2) (2017) 143–160.</mixed-citation>
      </ref>
      <ref id="ref-17">
        <label>17</label>
        <mixed-citation>A. Wilder-Smith, E. Massad, Estimating the number of unvaccinated Chinese workers against yellow fever in Angola, BMC Infect. Dis. 18 (1) (2018) 185.</mixed-citation>
      </ref>
      <ref id="ref-18">
        <label>18</label>
        <mixed-citation>M.E. Cano, G.A. Marti, A. Balsalobre, E. Muttis, E.A. Bruno, G. Rossi, M.V. Micieli, Database of Sabethes and Haemagogus (Diptera: Culicidae) in Argentina: Sylvatic Vectors of the Yellow Fever Virus, J. Med. Entomol. 58 (4) (2021) 1762–1770. https://doi.org/<pub-id pub-id-type="doi">10.1093/jme/tjab059</pub-id>.</mixed-citation>
      </ref>
      <ref id="ref-19">
        <label>19</label>
        <mixed-citation>E. Massad, F.A. Coutinho, L.F. Lopez, D.R. da Silva, Modeling the impact of global warming on vector-borne infections, Phys. Life Rev. 8 (2) (2011) 169–199.</mixed-citation>
      </ref>
      <ref id="ref-20">
        <label>20</label>
        <mixed-citation>C.P. Ferreira, S.P. Lyra, F. Azevedo, D. Greenhalgh, E. Massad, Modelling the impact of the long-term use of insecticide-treated bed nets on Anopheles mosquito biting time, Malar. J. 16 (1) (2017) 373.</mixed-citation>
      </ref>
      <ref id="ref-21">
        <label>21</label>
        <mixed-citation>A.T. Huang, S. Takahashi, H. Salje, D.A.T. Cummings, Assessing the role of multiple mechanisms increasing the age of dengue cases in Thailand, Proc. Natl. Acad. Sci. USA 119 (20) (2022) e2115790119.</mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>
