Next Article in Journal
Entomological Surveillance of Aedes Mosquitoes: Comparison of Different Collection Methods in an Endemic Area in RIO de Janeiro, Brazil
Previous Article in Journal
The Natural Alkaloid Tryptanthrin Induces Apoptosis-like Death in Leishmania spp.
Previous Article in Special Issue
Comparison of Immunochromatographic Test (ICT) and Filariasis Test Strip (FTS) for Detecting Lymphatic Filariasis Antigen in American Samoa, 2016
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Network Analysis of Lymphatic Filariasis Serology from Myanmar Shows Benefit of Adding Antibody Testing to Post-MDA Surveillance

1
College of Medicine and Dentistry, James Cook University, Cairns, QLD 4870, Australia
2
College of Public Health, Medical and Veterinary Sciences, James Cook University, Cairns, QLD 4870, Australia
3
School of Public Health, Faculty of Medicine, The University of Queensland, Brisbane, QLD 4072, Australia
4
Department of Medical Research, Ministry of Health and Sports, Yangon 11191, Myanmar
*
Author to whom correspondence should be addressed.
Trop. Med. Infect. Dis. 2022, 7(7), 113; https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7070113
Submission received: 30 April 2022 / Revised: 12 June 2022 / Accepted: 17 June 2022 / Published: 21 June 2022

Abstract

:
The elimination of lymphatic filariasis (LF) is achieved through repeated mass drug administration (MDA) of anti-filarial medications, which interrupts transmission and prevents new infections. Accurate transmission assessments are critical to deciding when to stop MDA. Current methods for evaluating transmission may be insufficiently sensitive, resulting in post-MDA resurgence. We, therefore, evaluated potential diagnostic testing scenarios for post-MDA surveillance. Data were used from two surveys (a household cluster and a cohort) conducted in an area of Mandalay Region, Myanmar, with ongoing transmission following several rounds of MDA. First, age- and sex-adjusted seroprevalence were estimated for the area using the household survey. Next, three Bayesian networks were built from the combined datasets to compare antigens by immunochromatic testing (ICT) and/or Og4C3 enzyme-linked immunosorbent assay (ELISA) and antibody (Ab) detection methods (Wb123 or Bm14 Ab ELISA). The networks were checked for validity and then used to compare diagnostic testing scenarios. The adjusted prevalence from the household survey for antigen, Wb123 Ab and Bm14 Ab were 4.4% (95% CI 2.6–7.3%), 8.7% (5.96–12.5%) and 20.8% (16.0–26.6%), respectively. For the three networks, the True Skill Statistic and Area Under the Receiver Operating Characteristic Curve for antigen, Wb123 and Bm14 Ab were 0.79, 0.68 and 0.55; and 0.97, 0.92 and 0.80, respectively. In the Bayesian network analysis, a positive case was defined as testing positive to one or more infection markers. A missed result was therefore the probability of a positive case having a negative test result to an alternate marker. The probability of a positive case prior to any testing scenario was 17.4%, 16.8% and 26.6% for antigen, Wb123 Ab and Bm14 Ab, respectively. In the antigen-only testing scenario, the probability of a missed positive LF result was 5.2% for Wb123 and 15.6% for Bm14 Ab. The combination of antigen plus Bm14 Ab testing reduced the probability of missing a positive LF case as measured by Wb123 Ab to 0.88%. The combination of antigen plus Wb123 Ab was less successful and yielded an 11.5% probability of a missed positive result by Bm14 Ab testing. Across scenarios, there was a greater discordance between Bm14 and both antigen and Wb123 Ab in the 1–10 age group compared to older ages. These findings suggest that the addition of Bm14 Ab improves the sensitivity of LF testing for current or past infection. The combination of antigen plus Bm14 Ab should therefore be considered for inclusion in post-MDA surveillance to improve the sensitivity of transmission surveys and prevent the premature cessation of MDA.

1. Introduction

Lymphatic filariasis (LF) is a vector-borne helminth infection that is transmitted by a variety of mosquito genera. Adult worms live in the human lymphatic system and produce microfilariae (Mf) which circulate in peripheral blood. Infection with LF damages the lymphatic system resulting in acute adenolymphangitis attacks and chronic swelling of limbs and genitals, which can become irreversible and cause severe morbidity and disability. In 2020, 51.4 million people in 72 endemic countries were estimated to be infected [1].
The Global Programme to Eliminate LF has been delivering annual mass drug administration (MDA) of anti-filarial medication for over 20 years. Population surveys for LF infection are crucial to the programme and are used at three time points: (1) to determine whether areas are endemic and need MDA, (2) to decide when MDA should stop, and (3) to determine whether elimination (interruption of transmission) can be validated. To date, 17 countries have been certified by the World Health Organization (WHO) to have eliminated LF as a public health problem. This requires passing a series of surveillance activities, including Transmission Assessment Surveys (TAS) on 6–7-year-old children [2,3].
Blood tests for surveillance both pre- and post-cessation of MDA and after validation of elimination are usually conducted using tests for circulating filarial antigen. Rapid tests available include the immunochromatographic test (ICT), now superseded by the filariasis test strip (FTS) [4]. An alternative is an enzyme-linked immunosorbent assay (ELISA) such as the Og4C3 test [5]. These tests all detect antigen that is produced by adult worms in the lymphatic system. The rapid tests are done on whole blood, whereas the ELISA antigen tests can be done using serum, plasma or dried blood spots [6]. Since adult worms need to be mature and mated before producing Mf, antigen detection tests are roughly 2–4 times more sensitive than examining blood slides for Mf, although the ratio between the two measures depends on prevalence and recent history of MDA [7].
It is known that antibodies to LF, including Bm14 and Wb123, appear very early in those infected by LF [8], are present at higher prevalence when compared to antigens [9], and persist for many months and even years after treatment [10]. Thus, the presence of antibodies against filarial antigens, especially in young children, may be a more sensitive method of detecting remaining transmission as elimination approaches, thus avoiding premature cessation of MDA.
South-East Asia is highly endemic for LF [11,12], but significant progress in reducing prevalence has been made in recent decades [13]. Three countries in the WHO South-East region (Thailand, Sri Lanka and Maldives) have reached validation of elimination [1].
Myanmar has one of the highest burdens of filariasis in South-East Asia, with 41 million (80% of the population) at risk [11,14]. It is distributed predominantly in the central and western dry zones, where baseline (pre-MDA) prevalence was 20–30% [11,15,16,17]. Wuchereria bancrofti is the sole cause of filariasis in the country, where it is transmitted by the mosquito Culex quinquefasciatus [11,18].
Myanmar commenced an MDA program with diethylcarbamazine and albendazole in 2001 [14,15]. Since then, coverage has been progressively expanded, reaching all endemic districts by 2015 [17]. The program has faced several challenges in medication supply, financing and concerns regarding adverse events, as well as insufficient program reach and participation in some areas [14,15,17,19]. As a result, six non-consecutive rounds of MDA over an 11-year period were completed in the study area at the time of the surveys from 2004 to 2014 inclusive. The program has led to considerable reductions in LF prevalence in many areas of Myanmar [15,17], but areas of persisting transmission remain [15,17,19].
Two recent surveys assessed the cross-sectional prevalence of LF in the Mandalay Region, a remaining transmission hotspot, in 2014–2015 [6,19,20,21,22,23]. Their complete serological data have not yet been reported, nor has any investigation into the utility of each individual infection marker in post-MDA surveillance.
Results from both surveys contribute to the current study, which aims to first generate a representative, age-stratified prevalence estimate for the study area in 2015 by Mf, antigen and antibody (Wb123 and Bm14 Ab). Next, to investigate the predictive value and relationship between different infection markers. Finally, to use Bayesian networks to evaluate the most efficient and sensitive testing methods for post-MDA and post-validation LF surveillance.

2. Methods

2.1. Sample Collection

Samples for serological testing were collected in two surveys in Mandalay Region, Myanmar, in 2014 and 2015. Firstly, from a representative community household cluster survey of those one year and older residing in 24 villages across four townships (Amarapura, Patheingyi, Tada-U, Wundwin) from January to March 2015 [19,20]. Secondly, from the baseline of a cohort study of antigen-positive and antigen-negative young people aged 10–21 years in one of these townships (Amarapura) in October 2014 [21,22].
Fingerprick whole blood samples from participants in both surveys were tested immediately for circulating filarial antigen by Binax Now Filariasis ICT cards (Alere Inc., Waltham, MA, USA) read at 10 min. In survey 1, night blood Mf slides (60 μL total in 3 lines) were taken on ICT positives where possible [19]; no slides were possible in survey 2. Slides were stained with Giemsa and read by two blinded microscopists as described in Dickson et al. [19]. Fingerprick whole blood samples (60 μL) for later antigen and antibody testing from dried blood spots were taken in both surveys from all participants. Survey 2 was a matched sample of antigen-positive and -negative young people screened for the longitudinal study who were invited to a follow-up visit. Venous blood using EDTA anticoagulant was taken at the follow-up visit for plasma. Samples were handled and transported to Australia in July 2015, as previously described [21].
Both surveys collected information on demographics (age and gender), place of residence, and history of taking MDA in prior years. Information on whether MDA was taken at the most recent opportunity (2014 MDA for Survey 1 and 2013 MDA for Survey 2) was available for all participants. Information on ever-taking MDA was available for Survey 1 only.

2.2. Selection of Samples for Testing

Due to limited resources, antibody assays were done on an age-stratified sample (n = 376) of dried blood spots from survey 1. An age-stratified sample was selected to ensure sufficient sample size by age to compare antigen and antibody prevalence in younger age groups, as well as overall population prevalence of serology markers. The samples were classified by age category and randomly selected within these age strata using Excel random number generator, sorting by the resulting number, and selecting of the following. Set A comprises all children 14 years and under; 25% of those aged 15–29 y, 30–44 y, 45–59 y and ≥60 y by group. In addition, all remaining ICT positives in survey 1 (n = 33, i.e., those not randomly selected for SET A) were tested for Bm14 and Wb123 antibodies and are referred to as Set B. Set C comprises all samples from survey 2 that were tested for Bm14 and Wb123 antibody (n = 80 comprising 40 matched pairs (antigen positive and negative) based on gender and within one year of age).

2.3. Serological Testing

Og4C3 ELISA for antigen was performed on eluted fingerprick dried blood spots from both surveys to provide comparable results. The methodology of testing and positive/negative determination was as previously described [6]. Briefly, Og4C3 antigen was detected using the TropBio Og4C3 kits (Cellabs Pty. Ltd., Manly, Sydney, Australia) at 1:16 dilution. The kit supplies positive and negative controls with a cut-off for positivity recommended based on panels of negative control sera. A standard curve was constructed for each ELISA plate using the kit-supplied positive controls assigning an arbitrarily high value of 32,768 units with subsequent 4-fold dilutions. The Softmax Pro v5 software (Molecular Devices, Sunnyvale, CA, USA) was used to convert optical density readings of the test samples to units based on a four-parameter curve. A single cut-off point of >32 units (representing the sixth point of the standard curve from subsequent dilutions, corresponding to the kit recommended cut-off) was used to determine positive readings. The Og4C3 ELISA is less sensitive for dried blood spots (DBS) than serum or plasma, likely due to the recommended dilution for each sample type, which is approximately four-fold higher dilution for the DBS version [6].
Bm14 antibody assays were performed on eluted fingerprick blood samples in both surveys, using Cellabs kits (Cellabs Pty. Ltd., Manly, Sydney, Australia) protocols at 1:100 dilution with methods as previously described [23]. Briefly, a standard curve was constructed on each plate using a highly-positive sample from Papua New Guinea, which was defined by an arbitrarily high value of 1000 antibody units with subsequent 2-fold dilutions. A single cut-off point of >125 units was used to determine positive readings based on previous literature [5,9,24,25].
Wb123 antibody assays were performed using InBios kits (InBios International Inc., Seattle, WA, USA) as described in [26,27]. The kit provides positive and negative reference controls. In previous work, samples were classed as positive if their average optical density (OD) ratio was 9 or more times the negative control (based on the average ratio of the low positive 1:5000 serum control to blank), and results applied (for repeated control sera) to standard curves constructed using two-fold dilutions of a highly positive sample as described for Bm14. The tests were done at 1:50 dilution using eluted fingerprick samples from survey 1 (since there were no venous samples available) and EDTA plasma at same dilution for survey 2 samples.

2.4. Statistical Analysis

Ninety-five percent confidence intervals were estimated using the proportion command with the default logit option. Proportions were compared using the Chi-squared test in Stata 16 (Stata Corporation, College Station, TX, USA). Age and gender standardization were done using the stdize option in the proportion command using the 2014 census data by five-year age and gender groups for the Mandalay region (Myanmar Information Management Unit; http://themimu.info/ (accessed on 1 December 2021)).
Definition of outcomes:
(1)
Antigen positive: Since the Og4C3 ELISA done with dried blood spots is less sensitive than ICT (see above), the definition of “antigen positive” was “positivity for ICT and/or Og4C3”.
(2)
Antibody positive: samples positive for either or both Bm14 and Wb123 antibody was defined as “antibody positive”.
(3)
Any infection marker positive: samples positive for at least one marker (antigen or either antibody) was defined as “any positive”.
(4)
Missed positive: A “missed positive” was defined as a sample that had a negative result but was positive for one or more alternate markers.

2.5. Bayesian Network Analysis

2.5.1. Bayesian Networks

Bayesian networks [28] are a graphical modeling method that represents the relationships between variables in terms of conditional probabilities. For variables (represented as nodes) that are dependent on another variable (the parent node), the conditional probability table (CPT) for each variable gives the probability of that node being in a given state based on the states of the parent nodes, for example, the probability that a sample is positive for antigen given a person’s age and whether they have ever taken MDA (Figure 1). Bayesian networks are increasingly being used in epidemiological studies because they are better suited than regression models for assessing complex systems and outcomes under different scenarios [29,30,31,32]. In this study, they were employed because they can efficiently compare the probability of a missed positive LF result with different diagnostic testing scenarios and assess the influence of various participant variables.

2.5.2. Network Design

Bayesian networks were designed to assess the relationship between different infection markers and their utility in post-MDA surveillance.
Three separate networks were created in Netica (Norsys Software Corporation, Vancouver, BC, Canada): one for each infection marker as the respective outcome node. A naïve network structure (in which the outcome node was the only parent of every other node) was used to assess the effect of covariates (e.g., age) on outcomes. To facilitate scenario analysis, the structure was then modified with additional links: (a) between infection markers to evaluate their relationship, and (b) between ‘Ever taken MDA medication’ and ‘Took MDA medication last year’ to preserve causal logic (e.g., only participants who had ever taken medication could have taken it in the last year). The network structure was validated by researchers with expertise in both filariasis epidemiology and Bayesian networks.

2.5.3. Network Testing

To evaluate the predictive validity of the Bayesian networks [33], each of the three models was tested using cross-validation. A random selection of 80% (n = 391) of the data was used to parameterize the CPTs of the Bayesian networks using the Entropy Maximisation algorithm. Model predictions for the remaining 20% of the records (n = 91) were then evaluated based on sensitivity, specificity and the True Skill Statistic [34], each of which used a classification cut-off of 50% (i.e., the option predicted as most probable was taken as the chosen outcome). A fourth metric, area under the curve (AUC) [35], was used to evaluate the models over a range of cut-offs. The four metrics were evaluated over 100 trials using a different random selection of records for each trial (with replacement between trials). Data from all three sample datasets were included in the analysis.

2.5.4. Scenario Analysis

To examine the relationships between the three infection markers, each network was parameterized using all data from the combined dataset (n = 482). Networks were then used to estimate the probability of a missed positive case (i.e., negative result in a sample positive by other infection markers) during different diagnostic testing scenarios. Additional scenarios were conducted to assess the effect of participant factors and dataset on testing outcomes. Scenarios were undertaken by setting the node states of the relevant nodes and noting the probability of a positive test result on the outcome node (for example, setting the value of Wb123 Ab and Antigen to “negative”) to calculate the probability of a positive Bm14 Ab result.

3. Results

Table 1 describes the characteristics of each dataset and its participants. Dataset A and B originate from Survey 1. Survey 1 was a cross-sectional survey of individuals ≥1 year old across four townships. Dataset A included all individuals ≤14 years and an age-stratified sample (25%) of those ≥15 years. Dataset B incorporates all ICT-positive individuals not already included in Survey 1. Survey 2 was the baseline for a cohort study of age- and sex-matched antigen-positive and -negative young people aged 10–21 years in one of the townships included in Survey 1. The different sampling methods for each dataset resulted in variable median age, infection marker prevalence and MDA participation between datasets.

3.1. Infection Marker Positivity in the Overall Sample

Figure 2 shows a comparison of positivity by ICT and Og4C3 in the total sample of 489 tested from SETS A, B and C. It demonstrates the lower sensitivity of Og4C3 done using dried blood spots at the recommended dilution, which was seen to a similar extent in both surveys. Ag positive was therefore defined as either ICT and/or Og4C3 positive.
Overall, 85 of 489 samples were classed as antigen-positive (17.4%, 95% confidence interval (CI) from 14.3 to 21.0%). After standardization for age and gender, the estimated proportion positive for Ag in the whole sample was 18.1% (95% CI: from 14.6 to 22.2%).
Figure 3 shows the antibody results. In the total sample (SETS A, B and C), the prevalence of antibodies to Bm14 (26.6%, 95% CI from 22.8 to 30.7%) was higher than for Wb123 (16.8%, 95% CI from 13.7 to 20.4%), which was of similar prevalence to antigen (17.4%, as noted above). For any antibody, 28.2% (95% CI from 24.4 to 32.4%) were positive, and 30.9% (95% CI from 26.9 to 35.1%) were positive for any marker of LF infection.
After adjustment for age and gender, the prevalence for Wb123, Bm14, any Ab and any positive were 17.4% (95% CI: from 14.0 to 21.5%), 27.9% (95% CI: from 23.7 to 32.5%), 29.6% (95% CI: from 25.3 to 34.2%) and 32.5% (95% CI: from 28.0 to 37.3%), respectively.
Combining SETS A and B, there were 47 Ag positive persons, of whom 42 had slides available; 17 of these were Mf positive, all for Wuchereria bancrofti. The Mf prevalence among Ag positives was 40.5% (95% CI from 25.6 to 56.7%). Seven of the thirteen antigen-positive, antibody-negative persons had slides taken, of which two (29%) were positive. No Mf were found in those under 15 years. Further details are given elsewhere [19].

3.2. Consumption of MDA Medication in the Overall Sample

Participants in all three datasets reported whether they had taken MDA at the last opportunity. The prevalence of LF markers of infection was much lower in those who had reported having taken MDA medication compared to those who had not (Figure S1). Datasets A and B also had information on whether participants had reported ever taking MDA, which demonstrated a similar association with LF markers.

3.3. Prevalence of Antigen and Antibody in the SET A Samples

The prevalence of Wb123 antibody was generally similar to antigen, while the prevalence for Bm14 or any antibody was higher. In the 376 persons sampled in SET A, overall crude antigen prevalence was 3.7% (95% CI from 2.1 to 6.2%). The corresponding crude prevalence for Wb123 Ab, Bm14 Ab, any Ab and any marker were 8.0% (from 5.6 to 11.2%), 18.9% (from 15.2 to 23.2%), 19.7% (from 16.0 to 24.02%) and 19.9% (from 16.2 to 24.3%), respectively.
After standardization for age and gender (necessary due to differing prevalence by age and gender, see below), the Ag prevalence in SET A was 4.4% (95% CI from 2.6 to 7.3%). The corresponding age and gender standardized prevalence for Wb123 Ab, Bm14 Ab, any Ab and any marker were, respectively, 8.7% (from 5.96 to 12.5%), 20.8% (from 16.0–26.6%), 21.9% (from 17.0 to 27.8%) and 22.5% (from 17.5 to 28.4%).
There were 230 females and 146 males in the SET A sample. The overall prevalence was significantly higher in males for antigen (6.2% versus 2.2%, p = 0.046), Bm14 Ab (24.7% vs. 15.2%, p = 0.023), any Ab (26.0% vs. 15.7%, p = 0.014) and any marker (26.7 vs. 15.7%, p = 0.009). A similar trend was seen for Wb123 Ab (10.3% vs. 6.5%, p = 0.191).
The age- and gender-specific prevalence of antigen in SET A are shown in Figure 4.

3.4. Bayesian Network Analysis for SETS A, B and C Combined

The three Bayesian networks used to evaluate the infection markers are shown in Figure 5 and Figure S2.
The average AUC and True Skill Statistic over 100 trials ranged from 0.80 to 0.97 and from 0.54 to 0.79, respectively (Figure S3). The performance of each network was considered sufficient to warrant their use in the scenario analysis. The Bayesian network for Bm14 Ab was less sensitive than the other two models, indicating that antigen and Wb123 were less able to predict Bm14 Ab.
Using the entire dataset, the probability of a positive test by different diagnostic testing scenarios is shown in Figure 6. Prior to defining any scenario, the probability of a positive result was 17.4%, 16.8% and 26.6% for antigen, Wb123 Ab and Bm14 Ab, respectively. If antigen was negative, there remained a sizeable probability of a positive result by Wb123 Ab (5.2%) and Bm14 Ab (15.6%). If either antibody was tested individually or together, and the results were negative, a considerable proportion was antigen-positive (5.9% for Wb123, 5.0% for Bm14, 3.7% for both). The testing combination of antigen plus Bm14 Ab detected the highest number of cases, with a probability of a missed positive as detected by Wb123 Ab of 0.88%.
Figure S4 displays the results of the different diagnostic testing scenarios by dataset and various participant factors (age, gender, township and history of MDA consumption).
By age, there was a higher probability of a positive result in older groups across all scenarios (Figure S4a). Of note, in the 1–10 age group, there was a greater discordance between Bm14 and both antigen and Wb123 Ab compared to older ages. Males demonstrated a greater probability of a positive result compared to females, but the relationship between markers was similar (Figure S4b).
Among townships, the likelihood of positive result was highest in Amarapura (the site of the cohort study) (Figure S4c). The concordance between Bm14 Ab and both antigen and Wb123 Ab was notably closer in Patheingyi and Wundwin compared to Amarapura and Tada-U.
Across all scenarios, the probability of a positive result was lower if MDA medication had ever been taken or within the last year (Figure S4d,e). Taking MDA medication was associated with a lower antigen and Wb123 Ab, relative to Bm14 Ab.
As expected, the representative sample (survey 1) had a lower likelihood of a positive result compared to the cohort baseline sample (survey 2), where half the individuals were antigen-positive (Figure S4f). All three datasets showed a moderate probability of antibody positivity if antigen-negative and that testing antigen plus Bm14 Ab resulted in the lowest probability of a missed positive result. Dataset C had a greater difference between antigen and Wb123 Ab, but each correlated more closely with Bm14 Ab.

4. Discussion

This study demonstrated a high seroprevalence of lymphatic filariasis in an area of Myanmar with ongoing transmission following several rounds of MDA. A Bayesian network framework was used to evaluate the relative utility of different infection markers for detecting signals of transmission in post-MDA surveillance. It found that antigen-based tests alone missed a considerable proportion of those with positive Wb123 and/or Bm14 Ab, including in the under-10 age group. Similarly, the antibody-based tests Bm14 and Wb123 alone, or in combination, also failed to detect some antigen-positive cases. The combination of antigen plus Bm14 Ab testing yielded the most sensitive marker of current or past infection. These findings suggest that post-MDA surveillance surveys would benefit from the addition of Bm14 Ab to more accurately detect potential signals of transmission and prevent the premature cessation of MDA.
Previous studies have reported on Bm14 antibodies in Myanmar for the SET C samples [23]. This paper is the first to report on Wb123 antibody prevalence and to compare antigen and antibody results for both Bm14 and Wb123 Ab in a population-representative sample (SET A). The high seroprevalence observed is similar to historic prevalence estimates in the region prior to the commencement of MDA [16]. The greater prevalence in males compared to females and with increasing age are consistent with those previously reported [25,36]. The marked antibody prevalence in males aged 41–50 years (60%) compared with a relatively lower prevalence in the 31–40 years group was unexpected and may reflect increased environmental exposure or sampling variation, although neither group was underrepresented.
A key finding of this study was that testing with antigen-based methods alone missed a sizable number of individuals with current or past infection. In post-MDA surveillance, the detection of any prior infection in children born since the commencement of MDA (11 years old in this context but 6–7 years for TAS surveys) is programmatically significant because it signals potential ongoing transmission. In the sample across all ages, the probability that antigen-negatives were positive for Wb123 or Bm14 Ab was 5.2% and 15.6%, respectively. Amongst those aged 10 years and younger, the relative difference between antigen and antibody tests was even greater.
The antibody tests for Bm14 and Wb123 alone, or in combination, also did not detect a considerable proportion of current infections in the Bayesian network. The finding of Mf in 29% of the antigen-positive, antibody-negative participants who had slides examined underscores the potential sensitivity limitations of these antibody assays. Bm14 Ab performed better than Wb123 Ab and detected most Wb123 positives (90%). The lower sensitivity of the Bm14 Ab Bayesian network signifies that antigen and Wb123 Ab were less able to predict Bm14 Ab and suggests that the Bm14 Ab may occur earlier during infection [8] or be more refractory to the effect of MDA. In contrast, Wb123 Ab tended to correspond with antigen tests across different diagnostic scenarios, including MDA consumption. This suggests that there are similarities between the expression of antigen and Wb123 during the filarial lifecycle and/or that they respond similarly to MDA.
The combination of antigen plus Bm14 Ab yielded the most sensitive marker of current or previous infection. This combination performed best both overall and across all subgroup scenarios. This suggests that Bm14 Ab testing should be considered for inclusion in TAS surveys to detect transmission more accurately and therefore prevent the premature cessation of MDA. Since Bm14 Ab is an ELISA-based assay that is not optimal for use in the field, a Bm14 rapid test will be required before it can be used in large-scale surveys. The available Brugia Rapid Test, which is used in non-W. bancrofti areas and detects antibodies to the similar antigen BmR1 [4], may provide a practical alternative. A comparison between Brugia Rapid and the Bm14 ELISA is needed before the rapid test can be recommended for use in TAS or similar surveys.
A notable finding was that MDA consumption was associated with a lower prevalence of all infection markers, with the greatest difference for antigen and Wb123 Ab. Antibody reductions after MDA have been observed previously [37] and probably relate either to a direct consequence of MDA’s filaricidal activity or to a lower probability of LF-exposure in areas where MDA uptake is higher. The greater discordance between Bm14 Ab compared to antigen and Wb123 Ab in Patheingyi and Wundwin townships in the current study likely indicates that transmission has reduced there to a greater degree than in Amarapura and Tada-U, where the markers remain more closely related.
These findings should be interpreted within the context of the study’s limitations. Firstly, resource constraints allowed for only a subset of participants from survey 1 to undergo serological testing, resulting in a conservative overall sample size (489). Secondly, the lack of clear timelines for longitudinal appearance and disappearance of each serological marker in individuals impedes clear conclusions regarding their specificity and sensitivity for detecting current infection. Thirdly, by defining a missed case as a negative result that was positive for one or more alternate markers, without the availability of a gold-standard test, the results do not account for the potential of false positives caused, for example, by testing cut-offs or parasite cross-reactivity. We also recognize that these samples were taken during 2014 and 2015 and may not reflect the current seroprevalence of infection markers. Nonetheless, the results from these robust Bayesian networks, across two separate surveys, demonstrate that either antigen or antibody testing alone may be insufficient for post-MDA surveillance.

5. Conclusions

This study has demonstrated that antigen or antibody testing in isolation is likely insufficient to detect transmission during post-MDA surveillance. The combination of antigen plus Bm14 Ab testing (or an equivalent antibody test) should be considered for inclusion in post-MDA surveys to increase the probability of detecting signals of transmission, particularly in low prevalence settings. Further studies are required to support this evidence and to evaluate comparable antibodies which are appropriate for large-scale field use.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/tropicalmed7070113/s1. Figure S1: Proportion of antigen and antibody positive by age group, shown by history of MDA taking last year; Figure S2: Bayesian networks with Bm14 Ab and Wb123 AB as the outcome nodes; Figure S3: Cross-validation testing of the three Bayesian Networks; Figure S4: Probability of a positive result in the Bayesian network analysis by different testing and subgroup scenarios.

Author Contributions

Conceptualization: B.F.R.D., J.J.R.M. and P.M.G.; sample collection: B.F.R.D., J.D. and L.B.; serological testing: J.J.R.M., K.S.A., K.M.H., M.R., A.A., S.R. and L.B.; data analysis: B.F.R.D., J.J.R.M., H.J.M. and P.M.G.; drafting of manuscript: B.F.R.D., H.J.M., J.J.R.M. and P.M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the generous support of the Australian Institute of International Affairs, Euan Crone, the College of Medicine and Dentistry of James Cook University, The Royal College of Pathologists of Australasia and private donor crowdfunding. James Cook University College of Public Health, Medical and Veterinary Sciences and GlaxoSmithKline contributed to survey costs in Myanmar and laboratory support in Cairns through the WHO Collaborating Centre for Lymphatic Filariasis, Soil Transmitted Helminths and Other Neglected Tropical Diseases (now called WHO Collaborating Centre for Vector-Borne and Neglected Tropical Diseases). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Myanmar Ministry of Health and Sports and The Human Research Ethics Committee at James Cook University (H5849 and H5261).

Informed Consent Statement

Written informed consent was obtained from all participants or their guardians.

Data Availability Statement

Data will be provided by authors upon request.

Acknowledgments

We thank all individuals who offered their time to participate in the study. We also acknowledge the other individuals who assisted in sample and data collection, in particular Koko Win Thu, Thet Wai Nwe, Tint Wai and Yi Yi Win from the Vector Borne Disease Control Unit, San San Win from the World Health Organization, Myint Shwe and the staff of the four township hospitals. We are indebted to Colleen Lau of the University of Queensland for advice and support on Bayesian network analysis. Lastly, we are very grateful to the NIH/NIAID Filariasis Research Reagent Resource Center for providing the ICT positive control.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Global programme to eliminate lymphatic filariasis: Progress report, 2020. Wkly. Epidemiol. Rec. 2021, 41, 497–508. [Google Scholar]
  2. WHO. Global Programme to Eliminate Lymphatic Filariasis. In Monitoring and Epidemiological Assessment of Mass Drug Administration; World Health Organization: Geneva, Switzerland, 2011. [Google Scholar]
  3. WHO. Validation of Elimination of Lymphatic Filariasis as a Public Health Problem; WHO: Geneva, Switzerland, 2017. [Google Scholar]
  4. Pastor, A.F.; Silva, M.R.; Dos Santos, W.J.T.; Rego, T.; Brandao, E.; de-Melo-Neto, O.P.; Rocha, A. Recombinant antigens used as diagnostic tools for lymphatic filariasis. Parasit Vectors 2021, 14, 474. [Google Scholar] [CrossRef] [PubMed]
  5. More, S.J.; Copeman, D.B. A highly specific and sensitive monoclonal antibody-based ELISA for the detection of circulating antigen in bancroftian filariasis. Trop. Med. Parasitol. 1990, 41, 403–406. [Google Scholar] [PubMed]
  6. Masson, J.; Douglass, J.; Roineau, M.; Aye, K.S.; Htwe, K.M.; Warner, J.; Graves, P.M. Relative Performance and Predictive Values of Plasma and Dried Blood Spots with Filter Paper Sampling Techniques and Dilutions of the Lymphatic Filariasis Og4C3 Antigen ELISA for Samples from Myanmar. Trop. Med. Infect. Dis. 2017, 2, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Berg Soto, A.; Xu, Z.; Wood, P.; Sanuku, N.; Robinson, L.J.; King, C.L.; Tisch, D.; Susapu, M.; Graves, P.M. Combining different diagnostic studies of lymphatic filariasis for risk mapping in Papua New Guinea: A predictive model from microfilaraemia and antigenaemia prevalence surveys. Trop. Med. Health 2018, 46, 41. [Google Scholar] [CrossRef] [PubMed]
  8. Hamlin, K.L.; Moss, D.M.; Priest, J.W.; Roberts, J.; Kubofcik, J.; Gass, K.; Streit, T.G.; Nutman, T.B.; Eberhard, M.L.; Lammie, P.J. Longitudinal monitoring of the development of antifilarial antibodies and acquisition of Wuchereria bancrofti in a highly endemic area of Haiti. PLoS Negl. Trop. Dis. 2012, 6, e1941. [Google Scholar] [CrossRef] [Green Version]
  9. Lau, C.L.; Sheel, M.; Gass, K.; Fuimaono, S.; David, M.C.; Won, K.Y.; Sheridan, S.; Graves, P.M. Potential strategies for strengthening surveillance of lymphatic filariasis in American Samoa after mass drug administration: Reducing ‘number needed to test’ by targeting older age groups, hotspots, and household members of infected persons. PLoS Negl. Trop. Dis. 2021, 14, e0008916. [Google Scholar] [CrossRef]
  10. King, C.L.; Weil, G.J.; Kazura, J.W. Single-Dose Triple-Drug Therapy for Wuchereria bancrofti-5-Year Follow-up. N. Engl. J. Med. 2020, 382, 1956–1957. [Google Scholar] [CrossRef]
  11. Dickson, B.F.R.; Graves, P.M.; McBride, W.J. Lymphatic Filariasis in Mainland Southeast Asia: A Systematic Review and Meta-Analysis of Prevalence and Disease Burden. Trop. Med. Infect. Dis. 2017, 2, 32. [Google Scholar] [CrossRef] [Green Version]
  12. Bizhani, N.; Hashemi Hafshejani, S.; Mohammadi, N.; Rezaei, M.; Rokni, M.B. Lymphatic filariasis in Asia: A systematic review and meta-analysis. Parasitol. Res. 2021, 120, 411–422. [Google Scholar] [CrossRef]
  13. Ramaiah Kapa, D.; Mohamed., A.J. Progress and impact of 20 years of a lymphatic filariasis elimination programme in South-East Asia. Int. Health 2021, 13, S17–S21. [Google Scholar] [CrossRef] [PubMed]
  14. World Health Organization, Regional Office for South-East Asia. Towards Eliminating Lymphatic Filariasis: Progress in the South-East Region (2001–2011); World Health Organization, Regional Office for South-East Asia: New Delhi, India, 2013. [Google Scholar]
  15. Aye, N.N.; Lin, Z.; Lon, K.N.; Linn, N.Y.Y.; Nwe, T.W.; Mon, K.M.; Ramaiah, K.; Betts, H.; Kelly-Hope, L.H. Mapping and modelling the impact of mass drug adminstration on filariasis prevalence in Myanmar. Infect. Dis. Poverty 2018, 7, 56. [Google Scholar] [CrossRef] [PubMed]
  16. UNDP; World Bank; WHO Special Programme for Research and Training in Tropical Diseases; WHO-UNICEF Joint Programme for Health Mapping; World Health Organization Division of Tropical Diseases. Research on Rapid Geographical Assessment of Bancroftian Filariasis; World Health Organization: Geneva, Switzerland, 1998. [Google Scholar]
  17. Win, K.M.; Tripathy, J.P.; Maung, T.M.; Oo, T.; Thi, A.; Lon, K.N.; Lin, Z. Rapid progress towards elimination of lymphatic filariasis in endemic regions of Myanmar as a result of 16 years of anti-filarial activities (2001–2016). Trop. Med. Health 2018, 46, 14. [Google Scholar] [CrossRef] [PubMed]
  18. de Meillon, B.; Grab, B.; Sebastian, A. Evaluation of Wuchereria bancrofti infection in Culex pipiens fatigans in Rangoon, Burma. Bull. World Health Organ. 1967, 36, 91–100. [Google Scholar] [PubMed]
  19. Dickson, B.F.R.; Graves, P.M.; Aye, N.N.; Nwe, T.W.; Wai, T.; Win, S.S.; Shwe, M.; Douglass, J.; Bradbury, R.S.; McBride, W.J. The prevalence of lymphatic filariasis infection and disease following six rounds of mass drug administration in Mandalay Region, Myanmar. PLoS Negl. Trop. Dis. 2018, 12, e0006944. [Google Scholar] [CrossRef] [PubMed]
  20. Dickson, B.F.R.; Graves, P.M.; Aye, N.N.; Nwe, T.W.; Wai, T.; Win, S.S.; Shwe, M.; Douglass, J.; Wood, P.; Wangdi, K.; et al. Risk factors for lymphatic filariasis and mass drug administration non-participation in Mandalay Region, Myanmar. Parasit Vectors 2021, 14, 72. [Google Scholar] [CrossRef] [PubMed]
  21. Douglass, J.; Graves, P.; Lindsay, D.; Becker, L.; Roineau, M.; Masson, J.; Aye, N.N.; Win, S.S.; Wai, T.; Win, Y.Y.; et al. Lymphatic Filariasis Increases Tissue Compressibility and Extracellular Fluid in Lower Limbs of Asymptomatic Young People in Central Myanmar. Trop. Med. Infect. Dis. 2017, 2, 50. [Google Scholar] [CrossRef] [Green Version]
  22. Douglass, J.; Dykes, L.; Kelly-Hope, L.; Gordon, S.; Leggat, P.; Aye, N.N.; Win, S.S.; Wai, T.; Win, Y.Y.; Nwe, T.W.; et al. Preventive chemotherapy reverses covert, lymphatic-associated tissue change in young people with lymphatic filariasis in Myanmar. Trop. Med. Int. Health 2019, 24, 463–476. [Google Scholar] [CrossRef]
  23. Masson, J.; Douglass, J.; Roineau, M.; Aye, K.S.; Htwe, K.M.; Warner, J.; Graves, P.M. Concordance between Plasma and Filter Paper Sampling Techniques for the Lymphatic Filariasis Bm14 Antibody ELISA. Trop. Med. Infect. Dis. 2017, 2, 6. [Google Scholar] [CrossRef] [Green Version]
  24. Gass, K.; Beau de Rochars, M.V.; Boakye, D.; Bradley, M.; Fischer, P.U.; Gyapong, J.; Itoh, M.; Ituaso-Conway, N.; Joseph, H.; Kyelem, D. A multicenter evaluation of diagnostic tools to define endpoints for programs to eliminate bancroftian filariasis. PLoS Negl. Trop. Dis. 2012, 6, e1479. [Google Scholar] [CrossRef]
  25. Lau, C.L.; Won, K.Y.; Becker, L.; Soares Magalhaes, R.J.; Fuimaono, S.; Melrose, W.; Lammie, P.J.; Graves, P.M. Seroprevalence and spatial epidemiology of Lymphatic Filariasis in American Samoa after successful mass drug administration. PLoS Negl. Trop. Dis. 2014, 8, e3297. [Google Scholar] [CrossRef] [PubMed]
  26. Lau, C.L.; Sheridan, S.; Ryan, S.; Roineau, M.; Andreosso, A.; Fuimaono, S.; Tufa, J.; Graves, P.M. Detecting and confirming residual hotspots of lymphatic filariasis transmission in American Samoa 8 years after stopping mass drug administration. PLoS Negl. Trop. Dis. 2017, 11, e0005914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Graves, P.M.; Sheridan, S.; Fuimaono, S.; Lau, C.L. Demographic, socioeconomic and disease knowledge factors, but not population mobility, associated with lymphatic filariasis infection in adult workers in American Samoa in 2014. Parasit Vectors 2020, 13, 125. [Google Scholar] [CrossRef] [PubMed]
  28. Fenton, N.; Neil, M. Risk Assessment and Decision Analysis with Bayesian Networks; CRC Press: New York, NY, USA, 2013. [Google Scholar]
  29. Lau, C.L.; Mayfield, H.; Lowry, J.H.; Watson, C.H.; Kama, M.; Nilles, E.J.; Smith, C.H. Unravelling infectious disease eco-epidemiology using Bayesian networks and scenario analysis: A case study of leptospirosis in Fiji. Environ. Model. Softw. 2017, 97, 271–286. [Google Scholar] [CrossRef]
  30. Mayfield, H.J.; Smith, C.S.; Lowry, J.H.; Watson, C.H.; Baker, M.G.; Kama, M.; Nilles, E.J.; Lau, C.L. Predictive risk mapping of an environmentally-driven infectious disease using spatial Bayesian networks: A case study of leptospirosis in Fiji. PLoS Negl. Trop. Dis. 2018, 12, e0006857. [Google Scholar] [CrossRef] [PubMed]
  31. Prodhan, G.; Fenton, N. Extending the range of COVID-19 risk factors in a Bayesian network model for personalised risk assessment. medRxiv 2020. [Google Scholar] [CrossRef]
  32. Lau, C.L.; Mayfield, H.J.; Sinclair, J.E.; Brown, S.J.; Waller, M.; Enjeti, A.K.; Baird, A.; Short, K.R.; Mengersen, K.; Litt, J. Risk-benefit analysis of the AstraZeneca COVID-19 vaccine in Australia using a Bayesian network modelling framework. Vaccine 2021, 39, 7429–7440. [Google Scholar] [CrossRef]
  33. Pitchforth, J.; Mengersen, K. A proposed validation framework for expert elicited Bayesian Networks. Expert Syst. Appl. 2013, 40, 162–167. [Google Scholar] [CrossRef] [Green Version]
  34. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  35. Fielding, A.B.J. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ. Conserv. 1997, 24, 38–49. [Google Scholar] [CrossRef]
  36. Michael, E.; Bundy, D.A.; Grenfell, B.T. Re-assessing the global prevalence and distribution of lymphatic filariasis. Parasitology 1996, 112 Pt 4, 409–428. [Google Scholar] [CrossRef] [PubMed]
  37. Tisch, D.J.; Bockarie, M.J.; Dimber, Z.; Kiniboro, B.; Tarongka, N.; Hazlett, F.E.; Kastens, W.; Alpers, M.P.; Kazura, J.W. Mass drug administration trial to eliminate lymphatic filariasis in Papua New Guinea: Changes in microfilaremia, filarial antigen, and Bm14 antibody after cessation. Am. J. Trop. Med. Hyg. 2008, 78, 289–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. An example Bayesian network showing a sample conditional probability table (CPT) for the child node (Antigen) based on the possible state combinations of the parent nodes (Age and MDA). With the network in its default state, the probabilities shown in each node represent the evidence in the dataset, i.e., 27.1% of the population were aged >40 years, 79.3% took MDA, and 18.6% were antigen-positive.
Figure 1. An example Bayesian network showing a sample conditional probability table (CPT) for the child node (Antigen) based on the possible state combinations of the parent nodes (Age and MDA). With the network in its default state, the probabilities shown in each node represent the evidence in the dataset, i.e., 27.1% of the population were aged >40 years, 79.3% took MDA, and 18.6% were antigen-positive.
Tropicalmed 07 00113 g001
Figure 2. Number of samples positive by ICT and Og4C3 in the combined datasets, Myanmar, from 2014 to 2015.
Figure 2. Number of samples positive by ICT and Og4C3 in the combined datasets, Myanmar, from 2014 to 2015.
Tropicalmed 07 00113 g002
Figure 3. Numbers of samples positive for antigen and antibody in the combined datasets.
Figure 3. Numbers of samples positive for antigen and antibody in the combined datasets.
Tropicalmed 07 00113 g003
Figure 4. Age–specific prevalence for Ag and antibody markers by age group and overall for SET A samples. (a) Both genders; (b) females; (c) males. In (a), the total is standardized for age and gender. In (b,c), the estimates are standardized for age.
Figure 4. Age–specific prevalence for Ag and antibody markers by age group and overall for SET A samples. (a) Both genders; (b) females; (c) males. In (a), the total is standardized for age and gender. In (b,c), the estimates are standardized for age.
Tropicalmed 07 00113 g004aTropicalmed 07 00113 g004b
Figure 5. Bayesian network with antigen as the outcome node.
Figure 5. Bayesian network with antigen as the outcome node.
Tropicalmed 07 00113 g005
Figure 6. Probability of a positive result in the Bayesian network analysis if tested negative under different diagnostic testing scenarios.
Figure 6. Probability of a positive result in the Bayesian network analysis if tested negative under different diagnostic testing scenarios.
Tropicalmed 07 00113 g006
Table 1. Characteristics of the study and participants by dataset.
Table 1. Characteristics of the study and participants by dataset.
Dataset
Survey 1Survey 2
ABC
Data Characteristics
Study dateJanuary–March 2015January–March 2015October 2014
Study typeRepresentative cross-sectionalRepresentative cross-sectionalCohort
Study populationCommunityCommunityCommunity
Participant selectionHousehold members ≥ 1 y.oHousehold members ≥ 1 y.o10–21 y.o age- and sex-matched antigen positive and negative
Sample selection≤ 14 y.o: all
≥15 y.o: age-stratified sample (25%)
All remaining ICT positivesAll
Participant Characteristics
Sample size3763380
Median age (years)
(interquartile range)
20.5 (9, 43)46.5 (35.5, 57)15 (12.5, 18.5)
Proportion female (%)615858
Township (n)Amarapura: 94Amarapura: 25Amarapura: 80
Patheingyi: 50Patheingyi: 1
Tada-U: 112Tada-U: 5
Wundwin: 120Wundwin: 2
Ever taken MDA medication (%)8164a
Took MDA medication last year (%)766145
Antigen positive (ICT and/or Og4C3) (%)3.710047.5
Wb123 Ab positive (%)8.069.736.3
Bm14 Ab positive (%)18.969.745.0
a Data not available.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dickson, B.F.R.; Masson, J.J.R.; Mayfield, H.J.; Aye, K.S.; Htwe, K.M.; Roineau, M.; Andreosso, A.; Ryan, S.; Becker, L.; Douglass, J.; et al. Bayesian Network Analysis of Lymphatic Filariasis Serology from Myanmar Shows Benefit of Adding Antibody Testing to Post-MDA Surveillance. Trop. Med. Infect. Dis. 2022, 7, 113. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7070113

AMA Style

Dickson BFR, Masson JJR, Mayfield HJ, Aye KS, Htwe KM, Roineau M, Andreosso A, Ryan S, Becker L, Douglass J, et al. Bayesian Network Analysis of Lymphatic Filariasis Serology from Myanmar Shows Benefit of Adding Antibody Testing to Post-MDA Surveillance. Tropical Medicine and Infectious Disease. 2022; 7(7):113. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7070113

Chicago/Turabian Style

Dickson, Benjamin F. R., Jesse J. R. Masson, Helen J. Mayfield, Khin Saw Aye, Kyi May Htwe, Maureen Roineau, Athena Andreosso, Stephanie Ryan, Luke Becker, Janet Douglass, and et al. 2022. "Bayesian Network Analysis of Lymphatic Filariasis Serology from Myanmar Shows Benefit of Adding Antibody Testing to Post-MDA Surveillance" Tropical Medicine and Infectious Disease 7, no. 7: 113. https://0-doi-org.brum.beds.ac.uk/10.3390/tropicalmed7070113

Article Metrics

Back to TopTop