 Research
 Open Access
 Published:
Estimating helminth burdens using sibship reconstruction
Parasites & Vectors volume 12, Article number: 441 (2019)
Abstract
Background
Sibship reconstruction is a form of parentage analysis that can be used to identify the number of helminth parental genotypes infecting individual hosts using genetic data on only their offspring. This has the potential to be used for estimating individual worm burdens when adult parasites are otherwise inaccessible, the case for many of the most globally important human helminthiases and neglected tropical diseases. Yet methods of inferring worm burdens from sibship reconstruction data on numbers of unique parental genotypes are lacking, limiting the method’s scope of application.
Results
We developed a novel statistical method for estimating female worm burdens from data on the number of unique female parental genotypes derived from sibship reconstruction. We illustrate the approach using genotypic data on Schistosoma mansoni (miracidial) offspring collected from schoolchildren in Tanzania. We show how the bias and precision of worm burden estimates critically depends on the number of sampled offspring and we discuss strategies for obtaining sufficient sample sizes and for incorporating judiciously formulated prior information to improve the accuracy of estimates.
Conclusions
This work provides a novel approach for estimating individuallevel worm burdens using genetic data on helminth offspring. This represents a step towards a wider scope of application of parentage analysis techniques. We discuss how the method could be used to assist in the interpretation of monitoring and evaluation data collected during mass drug administration programmes targeting human helminthiases and to help resolve outstanding questions on key population biological processes that govern the transmission dynamics of these neglected tropical diseases.
Background
Human helminthiases are caused by persistent parasitic infections associated with chronic poor health, morbidity and mortality. Many of the most globally prevalent and disabling helminth infections [1] are included in the World Health Organization (WHO) list of neglected tropical diseases (NTDs). These include the soiltransmitted helminthiases, ascariasis, hookworm and trichuriasis; the filariases, lymphatic filariasis and onchocerciasis, and the trematodiases, including the foodborne trematodiases and schistosomiasis [2,3,4]. They affect at least one billion people worldwide, predominantly the poorest of the poor [5]. The severity of helminthiases and the contribution to transmission of an infected individual is, at least in part, related to the number of parasites within an individual host at a particular time, the socalled worm burden [6, 7].
In many dioecious helminths, adult parasites live within their definitive host during their entire lifespan and fertile female parasites produce eggs or larvae (transmission stages) that are either released directly into the external environment (e.g. soiltransmitted helminth and trematode infections) or migrate to tissues where they can infect insect vectors (e.g. filarial infections). In humans, it is generally impossible to measure worm burdens directly in natural endemic communities (except for chemoexpulsion techniques for soiltransmitted helminths, see for examples [8,9,10,11] and more recently [12]) because adult parasites reside in inaccessible locations. Therefore, it is generally only feasible to perform indirect inference, by counting eggs or larvae. Yet this approach is hampered by high variability in counts made from a single individual (both from repeated samples taken at a single time and among samples taken over a short period of time, see for examples [9, 13,14,15,16]) and by potential densitydependent constraints on parasite fecundity that can render the wormegg/larva relationship highly nonlinear [6, 7] and in some cases geographically variable [17].
Parentage analysis using molecular data is used widely in molecular ecology [18,19,20]. Sibship reconstruction is a category of parentage analysis which can be used to estimate the number of parents when genetic data are available on offspring only [21,22,23,24]. Essentially, data on neutral genetic markers are used to divide offspring into groups of full siblings (monogamous mating) or groups of full siblings and half siblings (polygamous mating) to reconstruct and identify unique (male and/or female) parental genotypes. Hence the technique has potential to be used as a method of estimating worm burdens ([25], and see for examples [26,27,28]). We note that estimates of unique parental genotypes (worm burdens) can be further used in combination with estimates of the frequency of full and half siblings (also from sibship reconstruction) to estimate the effective population size (N_{e}) [29,30,31] and the effective number of breeders (N_{b}) [24] (using a random sample of individuals from a single cohort in a population with overlapping generations [30]). The (uncertain) relationship between N_{e} and census population size is discussed in detail in Palstra & Fraser [32].
Sibship reconstruction techniques depend fundamentally on the rules of Mendelian inheritance and frequently on the assumption that parental genotypes comprising a host’s worm burden exist at HardyWeinberg equilibrium [18]. Inbreeding, population bottlenecks and linkage disequilibrium of neutral markers proximate to genes under natural selection often violate these assumptions and are well discussed in the literature [33, 34], including parentage methods that relax the assumption of HardyWeinberg equilibrium and account for inbreeding [34,35,36]. But statistical questions that arise when using sibship reconstruction to infer the number of fecund adults have received little attention. While it is intuitive that in dioecious species, the number of unique female parental genotypes identified by sibship reconstruction will be limited (to a maximum) by the number of sampled offspring (eggs/larvae, i.e. each offspring can have only a single mother)—and consequently that the number of full sibling families is a minimum estimate of the female worm burden [24]—the statistical properties of how the number of unique parental genotypes in individual hosts relates to the underlying true number of fecund female adults (female worm burden) has not been explored. A clearer exposition of this statistical relationship is therefore essential to interpret accurately how parental genotypes inferred by sibship reconstruction relate to worm burdens.
Here, we develop a statistical approach to estimate the number and associated uncertainty of fecund female worms from data on the number of unique female parental genotypes identified by sibship reconstruction. We illustrate the technique using data from a recent study [28] where sibship reconstruction was used to identify unique parental genotypes of female Schistosoma mansoni infecting schoolchildren in Tanzania. We discuss potential applications in the context of monitoring and evaluation of mass drug administration programmes and resolving outstanding questions on the fundamental population biology of human helminthiases.
Methods
We define N as the number of fecund female schistosomes in a host. We also refer to N as the (female) worm burden. We assume that the N worms in a host are genetically unrelated and therefore each has a unique genotype of neutral genetic markers (e.g. microsatellites [28]). We define m as the number of miracidal offspring (hatched from eggs) sampled from a host and n as the number of unique female parental genotypes identified by sibship reconstruction [20, 21]. We further assume that the pool of eggs (hatched to miracidia) to sample from is large compared to N and that female worms within a host are equally fecund (i.e. each worm makes an equal contribution to the population of offspring). Under these assumptions, the probability of identifying n female parental genotypes from a sample of m miracidia is described by the unique items distribution [37],
where (N)_{n} is the falling factorial,
and \(\left\{ {\begin{array}{*{20}c} m \\ n \\ \end{array} } \right\}\) is a Stirling number of the second kind.
The expected value of n is given by [37],
such that the bias in E(n) as an underestimate of N expressed as a percentage is [E(n)−N]/N.
Inference on N from n identified parental genotypes is derived from the posterior probability f(Nn,m) using Bayes’ theorem,
where m is omitted for brevity and f(N) denotes the prior probability of a host having a worm burden N (i.e. the prior probability of N fecund female worms). We sampled from the posterior distribution f(Nn) using the following sampling importance resampling algorithm [38, 39]:

i.
draw a random sample of N from an integer uniform distribution g(Nn) for N ∈ [n, N_{max}]

ii.
calculate weights associated with each value of N, given by w(N) = f(nN)f(N)

iii.
resample N with replacement using the weights calculated in step 2
We set N_{max} (the maximum conceivable number of fecund female worms) to 350 in accordance with the autopsy observations made by Cheever [40]. This required the factorial in Equation 1 to be evaluated using Stirling’s approximation for N > 170,
We illustrate our statistical approach using genotypic data on schistosomes from schoolchildren in Tanzania collected in 2005, 2006 and 2010 [28]. In this case, n corresponds to the inferred number of unique female S. mansoni genotypes within each individual child (estimated by sibship reconstruction using multiplexed microsatellite genotypic data), and m corresponds to the number of sampled miracidia per child (ranging from 1 to 20). The number of fecund female worms N was estimated using a weakly informative (uniform) and informative negative binomial priors (W = 45 and k = 0.5 or k = 1). The analysis was performed in R [41] version 3.5.1.
Results
Relationship between number of unique parental genotypes and worm burden
Intuitively, the observed number n of unique female parental genotypes is a biased underestimate of the true number of female parental genotypes N present within a host, here defined as the fecund female worm burden. Clearly, n cannot be greater than the number of sampled (miracidial) offspring m [24] (i.e. when no full or half sibling pairs are identified from a sample of m offspring each of the m offspring has a different mother). The degree of bias can be quantified using the properties of the unique items distribution [37] (Fig. 1) under the assumptions that the number of offspring is large compared to N and that female worms within a host are equally fecund (i.e. there is an equal probability of sampling offspring from any female, see Methods). The relationship between the expected number (mean) of unique female parental genotypes E(n), the number of sampled offspring m and the fecund female worm burden N is shown in Fig. 2. The bias in n as an estimate of N is strongly dependent on the ratio m/N, such that to achieve less than a 5% underestimate of N one typically requires a sample of at least three times as many offspring as fecund female worms (Fig. 2b).
Estimating worm burdens
The purpose of sibship reconstruction in this context is to estimate a host’s female worm burden when it is unknown (i.e. N is unknown). We make inference on N and associated uncertainty from its posterior distribution, given an observed number of unique parental genotypes n and sampled (miracidial) offspring m (see Methods for details). The expected value (mean) E(N) and 95% confidence intervals are depicted in Fig. 3 for increasing numbers of identified unique female parental genotypes n and sampled (miracidial) offspring m. This relationship is based on the mild (weakly informative) a priori belief that hosts cannot harbour more than 350 female worms, but are otherwise equally likely to harbour any number between 1 and 350 worms. This maximum was chosen based on an autopsy study [40] that counted adult female S. mansoni directly from 103 people (Fig. 4). Without this constraint, the upper confidence intervals in Fig. 3 as n → m would be unbounded, tending to infinity.
The autopsy data shown in Fig. 4 follow a pattern of overdispersion that is typical of helminths and are well described by a negative binomial distribution. Hence, when estimating N in a population, it may be reasonable to assume that a randomly selected individual is more likely to have a low than a high worm burden. We incorporated this a priori belief using an informative negative binomial prior distribution, parameterised with a mean W and overdispersion parameter k, exploring the influence of W and k on estimates of N (Fig. 5). The values W = 45 and W = 20 chosen for this analysis were informed by the autopsy data; the mean W = 45 was estimated by fitting a negative binomial distribution to the data and the mean W = 20 was chosen because the autopsies were undertaken on the bodies of patients who had been terminally ill, many of whom had suffered severe schistosomiasis. Thus, we considered that the sample (with a mean W = 45) was likely biased towards heavier schistosome infections. The values of k = 0.5 and k = 1 were chosen as plausible considering both the autopsy data (Fig. 4) and values typically observed for other helminthiases [6, 42]. It is clear from Fig. 5 that higher values of W result in higher estimated values of N and associated degrees of uncertainty (compare Fig. 5c, d with Fig. 5a, b). Higher assumed (a priori) aggregation (lower k) results similarly in higher estimated values of N and associated degrees of uncertainty (compare Fig. 5a, c with Fig. 5b, d).
Illustration using genotypic data on schistosomes from schoolchildren in Tanzania
Gower et al. [28] used multiplexed microsatellite genotypic data of miracidia (hatched from eggs) sampled from 151 schoolchildren in Tanzania to identify by sibship reconstruction unique female S. mansoni genotypes within each individual child. Fullpedigree likelihood methods were used to infer sibship using the COLONY software package [43] and the number of miracidia sampled per child ranged from 1 to 20. Samples were collected in 2005, 6 months before the start of mass drug administration with praziquantel, in 2006 and in 2010.
The estimated number of fecund female worms N from each child in 2005, 2006 and 2010 are shown in Fig. 6, using the weakly informative (uniform) and informative negative binomial priors (W = 45 and k = 0.5 or k = 1). The results show that the number of unique female genotypes n for each child was substantially lower than the corresponding estimated posterior number of fecund female worms, N. For identical values of n, the posterior mean E(N) and associated uncertainty decreases as the number of miracidia m increases. As n becomes close to m, the estimated posterior of N becomes dominated by the prior distribution.
The mean and degree of uncertainty of the posterior estimates of N were considerably greater using a weakly informative prior (Fig. 6a, i.e. using a uniform prior that individuals are equally likely to harbour any number of worms between 1 and 350), compared to those using informative negative binomial prior distributions. The particular parameterisation of the negative binomial prior results in different posterior distributions of N as illustrated using W = 45 and k = 0.5 (Fig. 6b) or W = 45 and k = 1 (Fig. 6c).
Discussion
Parentage analysis by sibship reconstruction can be used to identify unique parental genotypes from genetic information on helminth offspring, which are typically more accessible than adult parasites, particularly for human helminthiases. We have developed a novel statistical approach to estimate—with associated measures of uncertainty—the number of fecund adult helminths from the number of unique parental genotypes identified by sibship reconstruction. We have illustrated the approach using genotypic data on S. mansoni miracidia collected from schoolchildren in Tanzania, highlighting the critical importance of (miracidial) offspring sample sizes to the precision of female worm burden estimates. This work provides a statistical exposition of using parentage analysis to estimate helminth worm burdens and thus a step towards a more robust application of this technique.
The potential of sibship reconstruction approaches to estimate the number of fecund adult helminths infecting a host when direct counts are impossible has been recognised for at least 15 years [25], but has only been applied relatively recently [24, 26,27,28] (see also [44]). Hitherto, the majority of studies have only used sibship reconstruction to identify the number of unique (fecund) parental genotypes [26,27,28] or analogously the number of fullsibling families (indicating the minimum number of worm pairs present within a host) [24], but have rarely attempted to relate these results to the underlying number of fecund female worms. The approach outlined in this paper provides a method to quantify the uncertainty (with credible intervals) of these estimates and has a number of potential population biological and epidemiological applications.
Many of the human helminthiases on the WHO list of NTDs are targeted for control or elimination, principally by mass drug administration (MDA) of anthelmintics to affected communities [45]. Epidemiological studies and routine monitoring and evaluation activities to assess the impact of MDA typically rely on egg counts or other indirect measures of infection intensity (such as detection of eggs/larvae by more sensitive molecular methods, e.g. [12, 46]). But the interpretation of such data can be complicated by uncertain and nonlinear [6, 7] or geographically variable [17] relationships with worm burden. Estimates of worm burden could therefore provide a useful approach to complement and aid the interpretation of such data. For example, in populations where egg or larva counts are higher than expected after years of MDA (in schistosomiasis these are referred to as ‘hotspots’, see for example [47, 48]), worm burden estimates could be used to distinguish programmatic deficiencies (e.g. poor coverage or missed MDA rounds) associated with higher than expected worm burdens from population biological processes associated with asexpected low worm burdens but higher than expected egg/larva counts (as occurs from the relaxation of densitydependent fecundity as worm burdens are reduced [49]).
The importance of densitydependent fecundity in interpreting egg count data is particularly pertinent to schistosomiasis because it remains unclear whether this fundamental population biological process operates in either of the two most globally important species, S. mansoni and S. haematobium [40, 50,51,52,53]. The use of sibship reconstruction to estimate worm burdens could help resolve this longrunning debate by exploring the association between egg counts and inferred female worm burdens, albeit using robust statistical methodologies to account for the likely substantial degree of measurement error (uncertainty) associated with the covariate (worm burden) estimate ([54] and see for example [55]). Moreover, the approach could also be used to revisit the relationship between worm burden and the relatively new antigenbased diagnostics for schistosomiasis, inference on which has been hitherto restricted to comparison with egg counts (see [56] and references therein).
The sample size of (miracidial) offspring is, however, a key limitation to the precision with which worm burdens can be estimated. Defining adequate sample sizes to determine populationlevel genetic diversity of human helminthiases has received attention, and in particular for schistosomiasis [57, 58]. But the sampling strategy required to make individuallevel parentage inference is a different proposition. Clearly, the more parasite offspring that can be collected and genotyped, the more robust the inference on worm burdens and thus ideally one would collect and genotype as many offspring (eggs/larva) as possible. This is feasible for soiltransmitted helminthiases and trematodiases, including schistosomiasis, by whole stool sampling [59, 60] or by using much greater quantities of stool than are used for routine diagnosis (e.g. [61]). This is important both in lightly infected individuals, where routine methods for counting eggs may frequently find no infection due to poor sensitivity [62,63,64,65], and in heavily infected individuals where sampling effort may otherwise be compromised by the greater ease of collecting parasite material from smaller quantities of stool.
One potential sampling approach would be to sample offspring proportionally to parasitological or other molecular indicators of the intensity of infection per individual. Although indirect measures of quantification are unreliable indicators of worm burden [9, 12, 65], which indeed is a key premise of using parentage analysis to estimate worm burdens, they provide some, albeit noisy, information at an individual level with which to motivate desired sample sizes. For example, one could chose a minimum baseline offspring sample size of 10 (which would be sufficient to estimate with reasonable precision female worm burdens of 3–4) and increase this value in proportion to the percentile of observed intensity indicators (such as eggs per gram of faeces per individual) obtained from a group or population of individuals.
The choice of prior distribution is a further important consideration when estimating worm burdens from the results of sibship reconstruction. The most cautious approach is to assume that any number of worms is equally as likely, and indeed that there is no upper limit on how many worms a host may harbour (here 350 female worms was considered a maximum for S. mansoni based on the human autopsy data [40], Fig. 4, but also to bound the posterior distribution of female worm burden as n → m). This will inevitably lead to the greatest uncertainty in worm burden estimates (Fig. 6). Using a negative binomial prior is well justified based on the wealth of empirical evidence from a variety of human and animal helminth infections [6, 42]. The difficulty arises with parameterising this prior distribution; both the mean and the degree of overdispersion will undoubtedly vary considerably among settings and particular contexts. The best approach is likely to assemble estimates made under different plausible prior assumptions. But it is also important to note that the influence of the prior is most pronounced in cases where the number of identified parental genotypes approaches the number of sampled offspring. Therefore, in the majority of hosts (with low worm burdens), the choice of prior may have limited substantive impact if relatively large offspring sample sizes are achieved.
Aside from the statistical properties of the relationship between parental genotypes, the number of sampled offspring and the underlying (female) worm burden, there exist a variety of population biological and genetic assumptions not considered here explicitly that may affect the accuracy of sibship reconstruction [20]. Methods for identifying unique parental genotypes from offspring genetic data are probabilistic [21] and thus will have some inherent uncertainty. Accurate inference particularly depends on the assumption of either a monogamous or polygamous mating system. Parental genotypes were identified from the Tanzanian data under the common assumption that schistosomes are strictly monogamous [6, 66] (as opposed to the ubiquitous assumption of polygamy among other human helminthiases [6]) although in model systems mate changes and polygamy do occur [67]. Studies have also shown that mating competition occurs between different species (and even different genotypes), and increased polygamy exists in high selective pressure environments [68, 69]. The assumption of a strictly monogamous mating system implies that the number of identified female genotypes is equal to the number of male genotypes, and therefore, only fullsibs sharing the same mother and father can be present in the sampled offspring. Under the assumption of a polygamous mating system, halfsib pairs sharing the same mother but different fathers can exist. Therefore, if strict monogamy does not hold, halfsib pairs may be erroneously missed, potentially overestimating the number of unique parental genotypes because of a failure to infer sibling relationships among the sampled offspring. Notwithstanding, a recent empirical comparison of the results of sibship reconstruction under the assumption of either a monogamous or a polygamous schistosome mating system found the assumption to have relatively little impact [24].
An important statistical limitation of this work is the assumption of equally fecund female worms within a host, i.e. that each worm makes an equal contribution to the population of offspring. A previous study used sibship reconstruction to quantify family structure in miracidial offspring (i.e. numbers of miracidia per fullsib family) and indicated that the reproductive success of breeding schistosomes was likely skewed, resulting in differential representation of each family in the offspring pool [24]. The probability of identifying exactly n female parental genotypes from a sample of m miracidia will be affected by this unequal contribution as in turn will the estimated number of N parental genotypes. Changing this assumption to account for variable reproductive output will be explored in future studies.
The assumptions of Mendelian inheritance and HardyWeinberg equilibrium of parental genotypes [18] are also limitations of parentage analysis techniques (but see [34,35,36]) for methods that relax the assumption of HardyWeinberg equilibrium). These assumptions may be violated by many factors, particularly during MDA programmes [70,71,72]. Even before MDA, assumptions of random mixing may not be upheld if genetically related parasites are transmitted together in socalled ‘clumps’ or ‘packets’, as described for directlytransmitted helminth infections [73,74,75]. This process, alongside other host and genotypedependent immunity factors [76], would promote inbreeding [77] and departure from HardyWeinberg equilibrium. For schistosomiasis, although asexual reproduction within intermediate snail hosts may also seemingly enhance inbreeding within the definitive (human) hosts, the mixing in water bodies of freeliving cercariae released from snails will likely act as a counteractive balance.
In practice, simulation studies have shown that inbreeding or relatedness among male and female parents has only a small effect on the accuracy of sibship reconstruction [36]. Therefore, for dioecious species (like schistosomes) it has been recommended that sibship reconstruction is performed under the assumptions of HardyWeinberg equilibrium (no inbreeding), except when there is strong evidence to the contrary and the level of inbreeding is high [78]. Moreover, the statistical relationship between the estimated number of fecund female parasites and the number of unique parental genotypes identified from a finite sample of (miracidial) offspring will be unaffected by the specific assumptions used for sibship reconstruction. Notwithstanding, if inbreeding is suspected as being high, sibship reconstruction can be implemented with or without assumptions of HardyWeinberg equilibrium to see how the identified number of unique parental genotypes may differ. The statistical approach presented here could be applied to different numbers of identified genotypes to estimate the corresponding number (and associated uncertainty) of fecund female worms as a form of sensitivity analysis.
Ultimately, validation of the approach outlined in this work would require that predictions be tested against directly observed adult (female) worm burdens. This is possible for human soiltransmitted helminths because adult worms can be expelled and counted following anthelmintic treatment [8,9,10,11,12]. For schistosomiasis and many other helminthiases, validation in related animal parasites that can be counted by dissection (e.g. S. bovis infections in cattle) or laboratory rodent studies [79,80,81] is the most amenable option. A drawback of this approach is that the population processes that shape the population genetics, and therefore the validity of the core assumptions inherent to sibship reconstruction, may be different in animal compared to human populations, especially populations under the influence of MDA programmes or frequent treatment.
Conclusions
An emerging use of sibship reconstruction is to identify the number of unique parental genotypes of human helminth infections from genetic information on their offspring. This is an important application of a welldeveloped parentage analysis technique because in many human helminthiases it is not possible to access adult worms, while it is relatively easy to sample their offspring. We have developed a statistical method to relate the number of identified parental genotypes to the underlying fecund female worm burden, highlighting the critical importance of offspring sample size on the bias and precision of worm burden estimates. The technique has potential applications in aiding the interpretation of routine monitoring and evaluation data collected during helminth control programmes and could contribute to resolving some outstanding population biological questions, particularly on the operation of densitydependent processes.
Availability of data and materials
Data supporting the conclusions are included within the article. The code is fully operational under R version 3.5.1 [41] and is freely available for download from https://github.com/minesneves
References
 1.
Kuper H. Neglected tropical diseases and disabilitywhat is the link? Trans R Soc Trop Med Hyg. 2019. https://doi.org/10.1093/trstmh/trz001.
 2.
Hotez P, Molyneux DH, Fenwick A, Kumaresan J, Sachs SE, Sachs JD, Savioli L. Control of neglected tropical diseases. New Engl J Med. 2007;6:1018–27.
 3.
Hotez PJ, Brindley PJ, Bethony JM, King CH, Pearce EJ, Jacobson J. Helminth infections: the great neglected tropical diseases. J Clin Invest. 2008;118:1311–21.
 4.
Neglected Tropical Diseases. https://www.who.int/neglected_diseases/diseases/en/.
 5.
Hotez PJ, Fenwick A, Savioli L. Rescuing the bottom billion through control of neglected tropical diseases. Lancet. 2009;373:1570–5.
 6.
Anderson RM, May RM. Infectious diseases of humans: dynamics and control. Oxford: Oxford University Press; 1991.
 7.
Anderson RM, May RM. Helminth infections of humans: mathematical models, population dynamics, and control. Adv Parasitol. 1985;24:1–101.
 8.
Forrester JE, Scott ME. Measurement of Ascaris lumbricoides infection intensity and the dynamics of expulsion following treatment with mebendazole. Parasitology. 1990;100:303–8.
 9.
Anderson RM, Schad GA. Hookworm burdens and faecal egg counts: an analysis of the biological basis of variation. Trans R Soc Trop Med Hyg. 1985;79:812–25.
 10.
Elkins DB, HaswellElkins M, Anderson RM. The epidemiology and control of intestinal helminths in the Pulicat Lake region of southern India. i. Study design and pre and posttreatment observations on Ascaris lumbricoides infection. Trans R Soc Trop Med Hyg. 1986;80:774–92.
 11.
Bundy DA, Thompson DE, Cooper ES. Rate of expulsion of Trichuris trichiura with multiple and single dose regimens of albendazole infection. Trans R Soc Trop Med Hyg. 1985;79:641–4.
 12.
Easton AV, Oliveira RG, O’Connell EM, Kepha S, Mwandawiro CS, Njenga SM, et al. Multiparallel qPCR provides increased sensitivity and diagnostic breadth for gastrointestinal parasites of humans: fieldbased inferences on the impact of mass deworming. Parasit Vectors. 2016;9:38.
 13.
Picq JJ, Jardel JP. A method of evaluating microfilaria densities of Onchocerca volvulus Leuckart, 1893, in onchoceriasis patients assessment of microfilarial densities by site and levels of prevelance in skin biopsies; variations of microfilarial densities over a 24 hour period). Bull World Health Organ. 1974;51:145–53.
 14.
Hall A. Quantitative variability of nematode egg counts in faeces: a study among rural kenyans. Trans R Soc Trop Med Hyg. 1975;75:682–7.
 15.
Utzinger J, Booth M, N’Goran EK, Müller I, Tanner M, Lengeler C. Relative contribution of daytoday and intraspecimen variation in faecal egg counts of Schistosoma mansoni before and after treatment with praziquantel. Parasitology. 2001;122:537–44.
 16.
Engels D, Sinzinkayo E, Gryseels B. Daytoday egg count fluctuation in Schistosoma mansoni infection and its operational implications. Am J Trop Med Hyg. 1996;54:319–24.
 17.
Hall A, Holland C. Geographical variation in Ascaris lumbricoides fecundity and its implications for helminth control. Parasitol Today. 2000;16:540–4.
 18.
Ashley MV, BergerWolf TY, Caballero IC, Chaovalitwongse W, DasGupta B, Sheikh S. Full sibling reconstruction in wild populations from microsatellite genetic markers. In: Computational Biology: New Research. Hauppauge, New York: Nova Publishers; 2009. p. 231–58.
 19.
Jones AG, Ardren WR. Method of parentage analysis in natural populations. Mol Ecol. 2003;12:2511–23.
 20.
Jones AG, Small CM, Paczolt KA, Ratterman NL. A practical guide to methods of parentage analysis. Mol Ecol Resour. 2010;10:6–30.
 21.
Wang JL. Sibship reconstruction from genetic data with typing errors. Genetics. 2004;166:1963–79.
 22.
Smith BR, Herbinger CM, Merry HR. Accurate partition of individuals into fullsib families from genetic data without parental information. Genetics. 2011;158:1329–38.
 23.
Thomas SC, Hill WG. Sibship reconstruction in hierarchical population structures using Markov chain Monte Carlo techniques. Genet Res. 2002;79:227–34.
 24.
Steinauer M, Christie M, Blouin M, Agloa L, Mwangi I, Maina G, et al. Noninvasive sampling of schistosomes from humans requires correcting for family structure. PLoS Negl Trop Dis. 2013;7:2456.
 25.
Criscione CD, Poulin R, Blouin MS. Molecular ecology of parasites: elucidating ecological and microevolutionary processes. Mol Ecol. 2005;14:2247–57.
 26.
Aemero M, Boissier J, Climent D, Moné H, Mouahid G, Berhe N, et al. Genetic diversity, multiplicity of infection and population structure of Schistosoma mansoni isolates from human hosts in ethiopia. BMC Genetics. 2015;16:137.
 27.
Gower CM, Gabrielli A, Sacko M, Dembelé R, Golan R, Emery A, et al. Population genetics of Schistosoma haematobium: development of novel microsatellite markers and their application to schistosomiasis control in Mali. Parasitology. 2011;138:978–94.
 28.
Gower CM, Gehre F, Marques SR, Lamberton PHL, Lwambo N, Webster J. Phenotypic and genotypic monitoring of Schistosoma mansoni in Tanzanian schoolchildren five years into a preventative chemotherapy national control programme. Parasit Vectors. 2017;10:593.
 29.
Ackerman MW, Hand BK, Waples RK, Luikart G, Waples RS, Steele CA, et al. Effective number of breeders from sibship reconstruction: empirical evaluations using hatchery steelhead. Evol Appl. 2017;10:146–60.
 30.
Wang J. A new method for estimating effective population sizes from a single sample of multilocus genotypes. Mol Ecol. 2009;18:2148–64.
 31.
Criscione CD, Blouin MS. Effective sizes of macroparasite populations: a conceptual model. Trends Parasitol. 2005;21:212–7.
 32.
Palstra FP, Fraser DJ. Effective/census population size ratio estimation: acompendium and appraisal. Ecol Evol. 2012;2:2357–65.
 33.
Wang J. Unbiased relatedness estimation in structured populations. Genetics. 2011;187:887–901.
 34.
Huang K, Mi R, Dunn DW, Wang T, Li B. Performing parentage analysis in the presence of inbreeding and null alleles. Genetics. 2018;210:1467–81.
 35.
Wang J. A simulation module in the computer program colony for sibship and parentage analysis. Mol Ecol Resour. 2013;13:734–9.
 36.
Wang J, Santure AW. Parentage and sibship inference from multilocus genotype data under polygamy. Genetics. 2009;181:1579–94.
 37.
Mendelson AF, Zuluaga MA, Hutton BF, Ourselin S. What is the distribution of the number of unique original items in a bootstrap sample? https://arxiv.org/abs/1602.05822.
 38.
Rubin DB. Using the sir algorithm to simulate posterior distributions. In: Bernardo JM, DeGroot MH, Lindlep DV, Smith AFM, editors. Bayesian Statistics 3. Oxford: Oxford University Press; 1988. p. 395–402.
 39.
Smith AFM, Gelfand AE. Bayesian statistics without tears: a samplingresampling perspective. Am Stat. 1992;46:84–8.
 40.
Cheever AW. A quantitative postmortem study of schistosomiasis mansoni in man. Am J Trop Med Hyg. 1968;18:38–64.
 41.
R Core Team: R: A Language and Environment for Statistical Computing. Vienna, Austria. R Foundation for Statistical Computing. 2018. https://www.Rproject.org/.
 42.
Wilson K, Bjørnstad ON, Dobson AP, Merler S, Poglayen G, Randolph SE, et al. Heterogeneities in macroparasite infections—patterns and processes. In: Hudson PJ, Rizzoli A, Grenfell BT, Heesterbeek H, Dobson AP, editors. The Ecology of Wildlife Diseases. Oxford: Oxford University Press; 2002. p. 6–44.
 43.
Jones OR, Wang J. Colony: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour. 2010;10:551–5.
 44.
Hildebrandt J, Eisenbarth A, Renz A, Streit A. Reproductive biology of Onchocerca ochengi, a nodule forming filarial nematode in zebu cattle. Vet Parasitol. 2014;205:318–29.
 45.
WHO. Accelerating work to overcome the global impact of neglected tropical diseases: a roadmap for implementation. Geneva: World Health Organization; 2012. https://www.who.int/neglected_diseases/NTD_RoadMap_2012_Fullversion.pdf.
 46.
Thiele EA, Cama VA, Lakwo T, Mekasha S, Abanyie F, Sleshi M, et al. Detection of Onchocerca volvulus in skin snips by microscopy and realtime polymerase chain reaction: implications for monitoring and evaluation activities. Am J Trop Med Hyg. 2016;94:906–11.
 47.
Wiegand RE, Mwinzi PNM, Montgomery SP, Chan YL, Andiego K, Omedo M, et al. Persistent hotspot of Schistosoma mansoni infection in a fiveyear randomized trial of praziquantel preventative chemotherapy strategies. J Infect Dis. 2017;216:1425–33.
 48.
Kittur N, Binder S, Campbell CH, King CH, Kinung’hi S, Olsen A, et al. Defining persistent hotspots: areas that fail to decrease meaningfully in prevalence after multiple years of mass drug administration with praziquantel for control of schistosomiasis. Am J Trop Med Hyg. 2017;97:1810–7.
 49.
Churcher TS, Filipe JA, Basànēz MG. Density dependence and the control of helminth parasites. J Anim Ecol. 2006;75:1313–20.
 50.
Medley G, Anderson RM. Densitydependent fecundity in Schistosoma mansoni infections in man. Trans R Soc Trop Med Hyg. 1985;79:532–4.
 51.
Wertheimer SP, Vermund SH, Lumey LH, Singer B. Lack of demonstrable densitydependent fecundity of schistosomiasis mansoni: analyses of Egyptian quantitative human autopsies. Am J Trop Med Hyg. 1987;37:79–84.
 52.
Gryseels B, De Vlas SJ. Worm burdens in schistosome infections. Parasitol Today. 1996;12:115–9.
 53.
Polman K, De Vlas SJ, Van Lieshout L, Deelder AM, Gryseels B. Evaluation of densitydependent fecundity in human Schistosoma mansoni infections by relating egg counts to circulating antigens through Deming regression. Parasitology. 2001;122:161–7.
 54.
Carroll RJ, Ruppert D, Stefanski LA, Crainiceanu CM. Measurement error in nonlinear models. London: Chapman and Hall; 2006.
 55.
Walker M, Little M, Wagner KS, SoumbeyAlley E, Boatin BA, Basànēz M. Densitydependent mortality of the human host in onchocerciasis: relationships between microfilarial load and excess mortality. PLoS Negl Trop Dis. 2012;6:1578.
 56.
Prada JM, Touloupou P, Adriko M, Tukahebwa EM, Lamberton PHL, Hollingsworth TD. Understanding the relationship between egg and antigenbased diagnostics of Schistosoma mansoni infection pre and posttreatment in uganda. Parasit Vectors. 2018;11:21.
 57.
Jarne P, Théron A. Genetic structure in natural populations of flukes and snails: a practical approach and review. Parasitology. 2001;123:27–40.
 58.
French CTS, Basáñez MG, Norton AJ, Lwambo NJS. Reductions in genetic diversity of Schistosoma mansoni populations under chemotherapeutic pressure: the effect of sampling approach and parasite population definition. Acta Trop. 2013;128:196–205.
 59.
Barbosa L, Barros B, de Souza Rodrigues M, Silva LK, Reis MG, Blanton RE. The effect of sample size on estimates of genetic differentiation and effective population size for Schistosoma mansoni populations. Int J Parasitol. 2018;48:1149–54.
 60.
Krauth SJ, Coulibaly JT, Knopp S, Traoré M, N’Goran EK, Utzinger J. An indepth analysis of a piece of shit: distribution of Schistosoma mansoni and hookworm eggs in human stool. PLoS Negl Trop Dis. 2012;12:1969.
 61.
Wang T, Shrivastava J, Johansen MV, Zhang SQ, Wang FF, Webster JP. Does multiple hosts mean multiple parasites? population genetic structure of Schistosoma japonicum between definitive host species. Int J Parasitol. 2006;36:1317–25.
 62.
Lamberton PHL, Kabatereine NB, Oguttu DW, Fenwick A, Webster JP. Sensitivity and specificity of multiple KatoKatz thick smears and a circulating cathodic antigen test for Schistosoma mansoni diagnosis pre and postrepeatedpraziquantel treatment. PLoS Negl Trop Dis. 2014;8:3139.
 63.
da Frota SM, Carneiro TR, Queiroz JA, Alencar LM, Heukelbach J. Combination of KatoKatz faecal examinations and ELISA to improve accuracy of diagnosis of intestinal schistosomiasis in a lowendemic setting in Brazil. Acta Trop. 2011;120:138–41.
 64.
Colley DG, Binder S, Campbell C, King CH, Tchuem Tchuenté LA, N’Goran EK, et al. A fivecountry evaluation of a pointofcare circulating cathodic antigen urine assay for the prevalence of Schistosoma mansoni. Am J Trop Med Hyg. 2013;88:426–32.
 65.
Walker M, Hall A, Basáñez MG. Ascaris lumbricoides: new epidemiological insights and mathematical models. In: Holland CV, editor. Ascaris: the neglected parasite. Amsterdam: Elsevier; 2013. p. 155–201.
 66.
May RM. Togetherness among schistosomes: its effects on the dynamics of the infection. Math Biosci. 1977;35:301–43.
 67.
Steinauer ML. The sex lives of parasites: investigating the mating system and mechanisms of sexual selection of the human pathogen Schistosoma mansoni populations. Int J Parasitol. 2009;39:1157–63.
 68.
Webster BL, Southgate VR, Tchuenté LAT. Mating interactions between Schistosoma haematobium and S. mansoni. J Helminthol. 1999;73:351–6.
 69.
Webster BL, Southgate VR. Mating interactions of Schistosoma haematobium and S. intercalatum with their hybrid offspring. Parasitology. 2003;126:327–38.
 70.
Webster JP, Gower CM, Norton AJ. Evolutionary concepts in predicting and evaluating the impact of mass chemotherapy schistosomiasis control programmes on parasites and their hosts. Evol Appl. 2008;1:66–83.
 71.
Churcher TS, Schwab AE, Prichard RK, Basáñez M. An analysis of genetic diversity and inbreeding in Wuchereria bancrofti: implications for the spread and detection of drug resistance. PLoS Negl Trop Dis. 2008;2:211.
 72.
Diawara A, Halpenny CM, Churcher TS, Mwandawiro C, Kihara J, Kaplan RM, et al. Association between response to albendazole treatment and betatubulin genotype frequencies in soiltransmitted helminths. PLoS Negl Trop Dis. 2013;7:2247.
 73.
Cornell SJ, Isham VS, Grenfell BT. Stochastic and spatial dynamics of nematode parasites in farmed ruminants. Proc R Soc Lond B. 2004;271:1243–50.
 74.
Walker M, Hall A, Basáñez MG. Trickle or clumped infection process? An analysis of aggregation in the weights of the parasitic roundworm of humans, Ascaris lumbricoides. Int J Parasitol. 2010;40:1373–80.
 75.
Walker M, Hall A, Basáñez MG. Trickle or clumped infection process? a stochastic model for the infection process of the parasitic roundworm of humans, Ascaris lumbricoides. Int J Parasitol. 2010;40:1381–8.
 76.
Van den Broeck F, Meurs L, Raeymaekers JA, Boon N, Dieye TN, Volckaert FA, et al. Inbreeding within human Schistosoma mansoni: do hostspecific factors shape the genetic composition of parasite populations? Heredity (Edinb). 2014;113:32–41.
 77.
Cornell SJ, Isham VS, Smith G, Grenfell BT. Spatial parasite transmission, drug resistance, and the spread of rare genes. Proc Natl Acad Sci USA. 2003;100:7401–5.
 78.
Wang J. User’s Guide for Software COLONY Version 2.0.6.5. London: Zoological Society of London; 2018.
 79.
Bøgh HO, Zhu XQ, Qian BZ, Gasser RB. Scanning for nucleotide variations in mitochondrial DNA fragments of Schistosoma japonicum by singlestrand conformation polymorphism. Parasitology. 1999;118:73–82.
 80.
Merenlender AM, Woodruff DS, Upatham ES, Viyanant V, Yuan HC. Large genetic distance between Chinese and Philippine Schistosoma japonicum. J Parasitol. 1987;73:861–3.
 81.
He YX, Hu YQ, Yu QF, Ni CH, Xue HC, Qiu LS, Xie M. Strain complex of Schistosoma japonicum in the mainland of China. Southeast Asian J Trop Med Public Health. 1994;25:232–42.
Acknowledgements
We are extremely grateful to Dr Charlotte Gower for access to data originally published in Parasites & Vectors [28].
Funding
MIN acknowledges funding support from a Royal Veterinary College, University of London PhD Studentship.
Author information
Affiliations
Contributions
MIN, JPW and MW conceived and conducted the analysis. MIN and MW wrote the draft and MIN, JPW and MW wrote and approved the final manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Neves, M.I., Webster, J.P. & Walker, M. Estimating helminth burdens using sibship reconstruction. Parasites Vectors 12, 441 (2019). https://doi.org/10.1186/s1307101936871
Received:
Accepted:
Published:
Keywords
 Parentage analysis
 Sibship reconstruction
 Worm burden
 Schistosomiasis
 Neglected tropical diseases