A New Method to Predict the Effect of an Intervention in the Host Population to Reduce the Magnitude of an Outbreak of a Vector-Borne Infection

Technische Samenvatting: Een nieuwe methode om het effect van een interventie in de gastheerpopulatie te voorspellen om de omvang van een uitbraak van een door vectoren overdraagbare infectie te verminderen

Probleemstelling
De schatting van de effectiviteit van beheersmaatregelen voor door vectoren overdraagbare infecties, zoals dengue, wordt bemoeilijkt door het niet-stationaire karakter van deze ziekten. In tegen tegenstelling tot direct overdraagbare infecties, die vaak onder aannames van een stationaire toestand (steady-state) worden gemodelleerd, vertonen door vectoren overdraagbare infecties jaarlijkse uitbraken die worden gedreven door seizoensgebonden variaties in vectorpopulaties en menselijke mobiliteit. Standaard wiskundige methoden kampen met aanzienlijke uitdagingen wanneer de aanname van een stationaire toestand wordt losgelaten, omdat ze het oplossen van partiële differentiaalvergelijkingen vereisen met onbekende tijdsafhankelijke parameters (bijv. fluctuerende omvang van de vectorpopulatie en menselijke mobiliteitspatronen). Bovormendien is het voorspellen van de omvang van een specifieke jaarlijkse uitbraak moeilijk, maar openbare gezondheidsplanners hebben methoden nodig om de impact van interventies (zoals vaccinatie) te schatten, ongeacht of een uitbraak groot of klein is.

Methodologie
De auteurs stellen een deterministisch model voor dat de impact van een interventie op de omvang van een door vectoren overdraagbare uitbraak schat zonder te vertrouwen op aannames van een stationaire toestand. De methode rust op een specifieke empirische observatie: voor door Aedes overgedragen infecties in Braziliaanse stedelijke centra is de leeftijdsverdeling van de proportie gerapporteerde gevallen onafhankelijk van zowel de intensiteit van de transmissie als de geografische locatie.

De kerncomponenten van de methodologie zijn:

1. Data-input: Het model maakt gebruik van officieel gerapporteerde, op leeftijd gestratificeerde incidentiegegevens ($Inc_H(a, t)$) van het Nationaal Systeem voor Verplichte Melding van Infecties (SINAN).
2. Leeftijdsverdeling-invariantie: Het model gaat ervan uit dat de proportie gevallen per leeftijd, aangeduid als $\phi(a)$, een functie is van leeftijd alleen, onafhankelijk van het uitbraakjaar of de locatie. Dit stelt het model in staat om de totale omvang van de uitbraak te ontkoppelen van de leeftijdsstructuur van de gevallen.
3. Simulatie van Interventie: De methode simuleert een interventie (bijv. vaccinatie) die een fractie $p$ van de gevallen binnen een specifiek leeftijdsinterval $[a_1, a_2]$ "verwijdert".
   • Vector-impact: Het berekent de reductie in het totaal aantal geïnfecteerde muggen ($\mathcal{I}_{MT}$) na de interventie, aangeduid met een factor $\alpha(a_1, a_2)$. Deze factor wordt afgeleid van de reductie in menselijke gevallen, gebaseerd op de premisse dat het aantal geïnfecteerde vectoren proportioneel is aan het aantal menselijke gevallen.
   • Gastheer-impact: Het berekent de reductie in het totaal aantal vatbare mensen ($\mathcal{S}_{HT}$) na de interventie, aangeduid met een factor $\beta(a_1, a_2, T)$. Dit vereist kennis van de pre-interventie dichtheid van vatbare individuen, bij voorkeur verkregen via seroprevalentie-onderzoeken.
4. Berekening van Effectiviteit: De totale effectiviteit ($Eff_T$) van de interventie wordt berekend als:
   $$Eff_T = 1 - \alpha(a_1, a_2) \beta(a_1, a_2, T)$$
   Dit vertegenwoordigt de reductie in het totaal aantal gevallen ($\Lambda_T$) vergeleken met een scenario zonder interventie.
5. Groepsimmuniteitsdrempel: Het model schat ook het basisreproductiegetal ($R_0$) door een exponentiële functie aan te passen aan de initiële groeifase van een uitbraak. De kritieke proportie om te immuniseren ($p_c$) wordt vervolgens afgeleid met de standaardrelatie $p_c = 1 - 1/R_0$, gecorrigeerd voor de effectiviteit van het vaccin.

Casestudy en Resultaten
De methode werd geïllustreerd met dengue-gegevens uit São José do Rio Preto (SJRP), Brazilië, die de uitbraak van 2019 beslaan.

• Interventie-scenario: Een hypothetische vaccinatiecampagne gericht op 100% van de populatie in de leeftijd van 10–20 jaar werd gesimuleerd.
• Parameters: Gebruikmakend van seroprevalentiegegevens van Chiaravalloti-Neto et al. (2019) om de vatbare populatie te schatten, leverde de simulatie een interventie-effectiviteit ($Eff_T$) op van 28%. Dit resultaat was consistent, of het nu werd berekend via de afgeleide formule of door het gesimuleerde post-interventie aantal gevallen te vergelijken met het pre-interventie aantal.
• Groepsimmunitet: Voor de uitbraak van 2019 in SJRP werd de geschatte $R_0$ 5,76. Bijgevolg werd de groepsimmuniteitsdrempel berekend waarbij ongeveer 83% van de vatbare populatie geïmmuniseerd moet worden om de infectie te elimineren.
• Vaccin-simulatie: De auteurs simuleerden de impact van het Butantan-DV vaccin (74,7% effectiviteit) over verschillende leeftijdsintervallen (12–59 jaar). Ze toonden aan hoe verschillende dekkingsgraden en leeftijdsdoelstellingen zich vertalen in effectieve immuunbescherming, waarbij de trade-offs in vaccinatiestrategieën werden benadrukt.

Belangrijkste Bijdragen en Claims
Het artikel beweert een nieuwe formalisering te bieden voor het vergelijken van de omvang van uitbraken met en zonder interventie, onafhankelijk van de temporele variaties in transmissie-intensiteit. De primaire bijdragen zijn:

1. Onafhankelijkheid van Stationaire Toestand: De methode werkt zonder een stationaire toestand aan te nemen, waardoor deze toepasbaar is op de fluctuerende dynamiek van seizoensgebonden door vectoren overdraagbare uitbraken.
2. Robuustheid tegenover Uitbraakgrootte: Omdat de methode steunt op de invariante leeftijdsverdeling van gevallen, kan het de effectiviteit van interventies voor uitbraken van elke omvang (groot of klein) binnen de bestudeerde regio's schatten.
3. Strategische Optimalisatie: Het kader maakt het mogelijk om optimale vaccinatiestrategieën te bepalen door verschillende leeftijdsintervallen en dekkingsgraden te simuleren om de effectiviteit te maximaliseren.

Beperkingen en Betekenis
De auteurs benoemen expliciet enkele beperkingen en matigende claims met betrekking tot de toepasbaarheid van de methode:

• Data-afhankelijkheid: De methode leunt zwaar op de kwaliteit van officiële meldingssystemen. Bias in rapportage kan de conclusies compromitteren.
• Schatting van Vatbaarheid: Een kritieke beperking is de noodzaak van nauwkeurige gegevens over het aantal of de proportie vatbare individuen. In de praktijk vereist dit meestal seroprevalentie-onderzoeken, die logistiek moeilijk en duur zijn. Zonder de ware proportie van de vatbare groep te kennen, kan de parameter $\beta$ niet precies worden bepaald, wat de mogelijkheid beperkt om de exacte effectiviteit in real-world "blinde" interventies te berekenen.
• Vectorspecifiteit: De bevinding van de invariantie van de leeftijdsverdeling wordt beweerd geldig te zijn voor door Aedes overgedragen infecties in stedelijke centra, maar kan mogelijk niet gelden voor infecties die worden overgedragen door andere vectoren (bijv. malaria).
• Stadsgrootte: Er wordt opgemerkt dat de analyse geldig is voor middelgrote tot kleine steden; in grote steden kan de ruimtelijke heterogeniteit van uitbraken en menselijke beweging de uniformiteit van de leeftijdsverdeling compliceren.

Concluderend presenteert het artikel een deterministisch instrument om de effectiviteit van interventies tegen door vectoren overdraagbare uitbraken te schatten, mits betrouwbare leeftijdsgestratificeerde casusgegevens en seroprevalentiegegevens (om de vatbare pool te definiëren) beschikbaar zijn. Het claimt niet de omvang van uitbraken te voorspellen, maar eerder de relatieve reductie in de omvang van de uitbraak te kwantificeren die bereikt kan worden door specifieke interventiestrategieën.