|Home | About | Journals | Submit | Contact Us | Français|
A recently developed high-throughput technique that allows multi-locus microsatellite analysis of individual miracidia of Schistosoma mansoni was used to assess the levels of genetic diversity and population structure in 12 infrapopulations of the parasite, each infrapopulation derived from an infected school child from the Mwea area, central Kenya. The mean number of alleles per locus was in the range 8.22 – 10.22, expected heterozygosity in Hardy-Weinberg equilibrium was 0.68 – 0.70, and pair wise FST values ranged from 0.16 −3.98% for the 12 infrapopulations. Although the genetic diversity within each infrapopulation of S. mansoni in this area was generally high, low levels of genetic structure were observed, suggestive of high levels of gene flow among infrapopulations. Private alleles were found in 8 of the 12 infrapopulation, the highest number of private alleles recorded per infrapopulation was 3. Our data suggest that the level of gene flow among infrapopulations of S. mansoni in Mwea is extremely high thus, providing opportunity for spread of rare alleles, including those that may confer character traits such as drug resistance and virulence.
Schistosoma mansoni causes human intestinal schistosomiasis in Africa, west Asia, parts of South America, and the Caribbean region (Savioli et al., 1997; Morgan et al., 2005). Human intestinal schistosomiasis is a debilitating disease characterized by a variety of symptoms including abdominal pain, diarrhoea and in severe cases hepatosplenomegaly and oesophageal varices (van der Werf et al., 2003). Worldwide, over 200 million people, mostly children, harbour schistosome parasites (Bundy, 1997; Chitsulo et al., 2000; Steinmann et al., 2006).
Schistosomes undergo sexual reproduction in the definitive hosts (humans, other primates, rodents) and this type of reproduction allows for reassortment and the perpetuation of parasite genotypic diversity. The eggs laid by the adult female worm pass in the host’s faeces and each will hatch in water to release a miracidium. These free-swimming larvae must find and penetrate an appropriate freshwater snail, in the case of S. mansoni, a snail of the genus Biomphalaria. Once it has penetrated the snail, the miracidium transforms into a mother sporocyst which produces multiple daughter sporocysts through asexual reproduction. These in turn produce cercariae which are released into water and are infective to human hosts, thus completing the life cycle.
Although S. mansoni is considered a uniform species morphologically, differences between parasite strains or populations have been observed in a number of biological characteristics such as infectivity (Thiong’o et al., 1997; Sire et al., 1999), virulence (Imbert-Establet and Combes, 1986), response to treatment (Dias et al., 1982; Gryseels et al., 2006) or fecundity (Imbert- Establet and Combes, 1986). By investigating mitochondrial sequence diversity, Morgan et al. (2005) provided a global overview of the genetic diversity inherent in S. mansoni, identifying five distinct lineages mostly separated by geographic distance, and that revealed considerable genetic variation particularly among East African specimens. Understanding the genetic structure of S. mansoni will lead to a better understanding of the variation in natural populations and the transmission dynamics of schistosomes between hosts and across geography. However, the intravascular location of adult schistosomes makes them routinely unavailable for study thus hindering population genetic studies of S. mansoni in the human definitive host. Most of the population genetic studies of S. mansoni undertaken so far have relied upon parasite materials obtained by indirect sampling procedures that involve passage of the parasite through laboratory animals (Minchella et al., 1995; Barral et al., 1996; Sire et al., 1999; Prugnolle et al., 2002; Morgan et al., 2003; Wang et al., 2006; Agola et al., 2006). Such procedures present ethical and practical problems, and may also result in sampling biases such as bottlenecking, and/or host-induced selection pressure (Gower et al., 2007).
Multi-locus microsatellite analysis of individual schistosome miracidia obtained from infected human hosts has recently been used successfully to characterize S. mansoni populations (Gower et al., 2007; Steinauer et al., 2008; Beltran et al., 2008), circumventing the problems associated with indirect sampling procedures and also providing a new overall approach to study the biology of S. mansoni in the human host.
S. mansoni is endemic in Kenya (Ministry of Health, Kenya, unpublished data), and recent studies of the parasites using adult worms derived from passages through laboratory animals have shown significant genetic diversity in the Kenyan S. mansoni populations, and that there is evidence of restricted gene flow among geographic locations (Morgan et al., 2003, 2005; Agola et al., 2006). Using the sampling procedure for schistosome miracidia recently described by Steinauer et al. (2008), genetic diversity and population structure of S. mansoni was examined using 9 microsatellite loci in infected school children living in the Mwea area, central Kenya, where intense rice farming is undertaken under irrigation, and S. mansoni prevalence is 47% (Kihara et al., 2007). Children are of particular relevance to sample in this respect because they play a major role in transmission of schistosomiasis, and in general, they are more at risk to morbidity and mortality than adults (Waiyaki, 1987; Gillespie, 2001).
The study population comprised school children aged 6–15 years old living in the rice farming irrigation scheme of Mwea in central Kenya, located approximately 100km north-west of Nairobi (Fig 1). A model schistosomiasis control project based on use of praziquantel has been established in Mwea as part of the school health program in the area. Only children diagnosed to have S. mansoni infection based on a Kato-Katz procedure (Katz et al., 1972), and whose parents had given their consent for participation in the study were sampled. A total of twelve children from four schools were recruited into the study. School was chosen as a sampling unit based on the fact that children presumably get infected near school grounds and also because schools reflect the distinct geographical area from which children originate.
A fecal sample of approximately 150 gm was collected in a plastic container from each participating child. To isolate S. mansoni eggs, the fecal sample was placed in a screw capped tube and normal saline added to completely cover the sample. The tube was shaken thoroughly to prepare a uniform suspension. The suspension was then filtered through a series of nested sieves of sizes 710, 425, 212, 45 µm. The eggs were then washed off the last sieve using distilled water into a 1-litre conical flask and distilled water was added to fill up the flask. Parasite egg hatching was induced by exposure to bright light. To collect the hatched miracidia, the flask was covered with a black cloth for 30 min to allow the miracidia (which are phototropic) to swim to the top of the flask. Individual miracidia were then collected directly into each wells of a 96 well plate.
DNA was extracted from schistosomes (miracidia) using a modified method described by Truett et al. (2000). In the modified procedure, 96 well PCR plates were used for DNA extraction instead of PCR tubes. Briefly, 5µl of a HotSHOT lysis reagent (25mM NaOH, 0.2mM disodium ethylenediaminetetraacetic acid (EDTA), pH 12) was added into each well of a 96 well PCR plate, sealed with pierceable aluminum tape, and chilled during specimen collection. An individual S. mansoni miracidium was transferred in a 3 µl volume of water using a 10 µl micropipette into a PCR plate well containing the lysis reagent by piercing the foil with the pipette tip. Pipetting through the foil eliminated the possibility of adding miracidia into the same well twice. Once the miracidia were in the wells, the plates were centrifuged for 20 seconds at 16000 × g, sealed with a silicone sealing mat, and heated in a PCR machine with a heated lid at 95°C for 30 minutes. The plates were periodically removed from the PCR machine and centrifuged to ensure that the liquid containing the miracidia remained at the bottom of the wells. After heating, the plates were cooled, centrifuged, and 5µl of HotSHOT neutralizing reagent (40mM Tris–HCl at pH 5) added into each well. The plates were then sealed and stored below −20°C until DNA amplification was performed. Using this procedure many samples could be processed very rapidly.
A total of 9 previously published microsatellite loci were investigated: AF325695 (Curtis et al., 2001), AF202965, AF202966, AF202968, L46951, M85305, R95529 (Durand et al., 2000), AI395184 and BF936409 (Silva et al., 2006).
DNA amplification was carried out in a multiplex PCR format using two panels, P17 and P22, as described by Steinauer et al. (2008). Panel P17 used primers for AF325695, AF202965, AF202966, AF202968 and L46951 while P22 used primers for AI395184, BF936409, M85305 and R95529. The QIAGEN Multiplex PCR Kit (Qiagen) was used for PCR amplifications in 5.8 µl reactions according to the manufacturer’s directions. The thermal cycling profile included an initial denaturation step at 94 °C for 5 min, followed by 35 cycles of 45 s at 94 °C, annealing temperature specific for each multiplex panel (52 °C for P 17, 48 °C for P 22), 30 s at 72 °C and a final step of 7 min at 72 °C using an Eppendorf Mastercycler Systems thermocycler. Primers were labeled with 6-FAM, HEX and NED dyes (Applied Biosystems). Multiplexed PCR products were diluted in N, N’-dimethyl formamide with Genescan® – 500 [ROX 500] as an internal size standard and genotyped using an ABI3100 automated sequencer (Applied Biosystems) and scored with GeneMapper® v. 4.0 (Applied Biosystems) software. All genotype calls were verified manually.
The mean number of alleles per locus (MNA) detected in each infrapopulation and the expected heterozygosities were used as indicators of the genetic polymorphism within the populations under study. Differences both for the mean number of allele and the expected heterozygosity values were tested for samples of each population using Freiedman non parametric tests. The MNA and allele frequencies were calculated using the program MICROSATELLITE TOOLKIT (http:// animalgenomics.ucd.ie/sdepark/ms-toolkit/). Both the expected and observed heterozygosities were also calculated using MICROSATELLITE TOOLKIT and their statistical significance tested using the chi-square test and an α̍ of 0.05. GENEPOP 3.3 (Raymond and Rousset, 1997) was used to test for deviations from Hardy-Weinberg equilibrium (HW) using exact tests, testing the hypothesis that observed diploid genotypes are the product of random union of gametes. An exact test for linkage disequilibrium between pairs of loci was performed using the Markov chain method in GENEPOP web version 3.4 (Raymond and Rousset, 1997). A sequential Bonferroni correction was also performed to control type 1 error. MICRO-CHECKER software (Van Oosterhout et al., 2004) was used to identify loci that potentially are affected by presence of null alleles, long allele dropout, or scoring errors due to stuttering. Frequency of null alleles was estimated from the heterozygosity deficiency as described by Brookfield (1996) using the formula (HE - HO) / (1 + HE.). Mean estimates of FIS for each population and pairwise Fst between all population pairs were also calculated following the method of Weir and Cockerham (1984). Deviation of FIsand Fst values from zero was tested using a permutation test. All F statistics were carried out using FSTAT 126.96.36.199 (Goudet, 2000). The recoding method of Meirmans (2006) was used to standardize all measures of Fst and only standardized measures of Fst are reported. For PCA (Principle Component Analysis), genetic distance was calculated from the allele data and the genetic distance was plotted as PCA using GenAlEx (Peakall and Smouse, 2006). Phylogenetic trees were constructed using the neighbor-joining (NJ) method (Saitou and Nei, 1987) from matrices of Nei’s genetic, DA, distances (Nei et al., 1983). Bootstrap tests (Felsenstein, 1985) were performed with 1000 replicate resampling. These analyses were done using MEGA (Kumar et al., 2004). Correlation between genetic distance and geographic distance was implemented using GENEPOP 3.3 and the significance of the correlation coefficient was tested using the Mantel test (Mantel, 1967). In order to determine the usefulness of the different markers to distinguish between the twelve infrapopulations the statistical significance of the FST for each of the 9 loci and for each infrapopulation pair was calculated. An unbiased estimate of the P-value of a log-likelihood (G-based) exact test was performed using GENEPOP. Individual population assignment was also undertaken using GENECLASS 1.0.02 (Cornuet et al., 1999).
Schistosomes were defined at two hierarchical levels: 1) all miracidia from an individual patient were designated as an “infrapopulation”, 2) miracidia from all the patients within a school were treated as a “population” in the analysis. Twelve different infrapopulations were analysed, totaling 904 individual S. mansoni miracidia, each of which was examined at 9 loci. Total number of alleles per locus ranged from 2 for AF202966 to 25 for L46951. The largest mean number of alleles was found in the NYG infrapopulation (10.22), while the lowest mean number of alleles occurred in the NGC and MU infrapopulations (8.22). Private alleles were found in eight of the twelve infrapopulations, and all occurred at frequencies above 5%. In all infrapopulations, the multilocus probability tests showed departure from Hardy-Weinberg equilibrium (P < 0.001). Expected heterozygosity ranged from 68% in NGC to 70% in both the NYA and NYD infrapopulations (Table 1) and there were no significant differences between the expected and observed heterozygosities (x2 = 0.1849, DF = 1, P = 0.6671). Friedman non parametric test for both the mean number of alleles and the expected heterozygosity values were nonsignificant. Following bonferroni correction, all locus pair tests were found not to exhibit linkage disequilibrium.
A total of four populations (schools) were analysed. To remove any bias, the mean number of alleles per locus was corrected for sample size by random sampling of 55 individual miracidia with 99 replicate replacements. The largest mean number of alleles was found in the NY (Nyamindi school) population (13.67), while the lowest mean number of alleles occurred in the MU (Mukou school) population (8.22). Private alleles were found in all populations except the KR (Kirogo school) population, and most of them occurred in the NY population (Table 1). In all populations, the multilocus probability tests showed departure from Hardy-Weinberg equilibrium (P < 0.001). Expected heterozygosity ranged from 69% in KR to 70% in the NY, NG (Ngurubani school) and MU populations (Table 1) and there were no significant differences between the expected and observed heterozygosities (x2 = 0.085, DF = 1, P = 0.7706).
For infrapopulations, genetic distances (Da) were highest (0.1178) between NYD and NGB. The lowest Da was observed between NGC and the MU (0.051). A neighbor-joining tree based on Da genetic distances for the twelve infrapopulations of S. mansoni (Fig 2) showed that the NGB and NYC infrapopulations are separated from the other infrapopulations with a bootstrap support of 70%. For populations, genetic distances (DA) were highest (0.0994) between the MU population (Mukou primary school) and KR (Kirogo primary school). The lowest DA (0.0278) was observed between NG and NY populations.
For infrapopulations, all but one (NYF) had positive FIS values (Local fixation index which corresponds to the excess of homozygosity of individuals in a sub-population) that ranged from 0.01 in NYA to 0.15 in NGB (Table 1). Global testing showed that all but one infrapopulation (NYF) had positive FIS values that deviated significantly from Hardy-Weinberg equilibrium. The presence of null alleles was tested and the frequencies ranged from 0.05 to 0.056 which is insignificant (i.e < 0.10) (Brookfield, 1996). Pairwise FST values (measure of allele frequencies among subpopulations) are presented in Table 2 and the highest values were recorded between S. mansoni infrapopulations NYF and NYD (FST = 0.1241) while the lowest were between NYB and NYG (FST = 0.0023). Although FST values were low, pairwise FST between infrapopulations were significantly different (P < 0.001). We also tested the usefulness of different markers in infrapopulation differentiation. All markers were efficient at distinguishing between all the infrapopulations, and the level of differentiation, at all infrapopulation pairs, across the nine loci was significant (P < 0.001). To understand the clustering pattern of S. mansoni infrapopulations, we also carried out genetic distance based Principal Component Analysis (PCA) using GenAlEx. The PCA lacked clear geographical patterns suggesting the absence of strong substructure within the S. mansoni population within Mwea. Assignment of individual to populations using GENECLASS software also indicated the same pattern of random assignment of individuals to population further indicating a lack of genetic structure.
FST values between populations (defined by schools) were generally low. The highest pairwise FST values were recorded between S. mansoni from Mukou school (MU) and Ngurubani school (NG) (FST = 0.0156) while the lowest were between Nyamindi school (NY) and Kirogo school (KR)(FST = 0.0061). All population pairs were significantly differentiated (P < 0.001) but the levels of differentiation were generally low (Table 2).
We applied the Mantel test to assess the correlation between genetic differentiation (FST/ 1- FST) and the logarithm of linear geographical distances (Ln distance) for all S. mansoni collection localities (schools). The analysis revealed a non-significant trend between both variables (r = 0.698; P = 0.383).
S. mansoni infrapopulations (each infrapopulation derived from a single host) from Mwea in central Kenya have a high level of genetic diversity as seen by the large number of alleles detected in each infrapopulation, in agreement with observations of S. mansoni populations from different parts of Kenya (Morgan et al., 2005; Agola et al., 2006). Even though the present study was carried out on a small geographic scale, high levels of genetic diversity of the parasite were still observable. The level of genetic diversity in schistosome populations is likely to be influenced by a variety of factors. On the one hand, repeated exposure of children in a selected number of transmission sites can result in the recirculation of related populations of schistosomes. If infection of humans and snails occurs at isolated transmission sites, then this would lead to development of distinct schistosome lineages as mutation and drift affects isolated populations within hosts and across geography. Also, acquired immunity and /or incomplete host susceptibility may also serve to limit the variety of schistosome seen in an infrapopulation (Brouwer et al., 2001). The rapid turnover of the infected snails induces a temporal succession and renewal of the different larval genotypes within a site and spatial aggregation of infected snails concentrate cercariae with different genotypes within a limited standing water area favouring multigenotype infection of definitive hosts (Theron et al., 2004). On the other hand, opportunities for genetic exchange resulting from snail and especially human movements may lead to coexistence of a large number of genotypes in an area (Brouwer et al., 2001).
Schistosome populations were significantly structured at two hierarchical scales, across geography and among human hosts. To assess the levels of differentiation between the infrapopulations analysed in this study, FST was calculated for all possible infrapopulations pairs. FST was considered appropriate because genetic drift is assumed to be the main factor in genetic differentiation among closely related infrapopulations or for short term evolution (Reynolds et al., 1983; Weir, 1990). Generally, FST calculated between all infrapopulation pairs was low but significant (P < 0.001). One possible explanation for the differences among patients is the presence of family substructure. We did not see negative FIS values as expected if there was strong family structure in the samples. Hence there were very low levels of genetic differentiation within S. mansoni infrapopulations in the Mwea area, suggesting high levels of gene flow between all the populations analysed. However, some infrapopulations from school children from different schools had slightly higher FST value as compared to the other pair wise FST values indicating slightly restricted gene flow (0.0392; between the NYD and NGC populations). A major contributor to gene flow is most likely to be the free movement of the school children in the study area. These children potentially sample a number of transmission sites, thereby becoming ‘genetic mixing bowls’ for the parasite (Curtis and Minchella, 2000). The parasites also live long enough to integrate infections across multiple seasons. The low levels of genetic differentiation can also be attributed to the nature of human population sampled in this study. Cross-sectional surveys of communities that will include other age groups left out in the present study are more likely to reveal population structure than, for example, a population of schoolchildren (Curtis et al., 2002). Efforts were also made to determine whether there was a correlation between genetic differentiation and geographical distance. The calculated interpopulation similarity indices showed that the populations do not follow a pattern of isolation by distance. Similar studies in Brazil (Thiele et al., 2008) also observed in their study of S. mansoni population a lack of significant isolation by distance.
The PCA analysis with GenAlEx failed to find significant structuring within the data, in contrary to the other analyses. This apparent discordance suggests that although there is significant structure, it is not necessarily strong. This is further evidenced by the relatively low genetic distances between all the infrapopulations (0.049- 0.1178) compared to those detected in our previous studies (0.203- 0.609) of S. mansoni populations analysed across the country (Agola et al., 2006), with the highest being recorded between NGB and NYE; distances were also low between populations such as KR and the MU, suggesting that these populations probably share a common gene pool. We hypothesize that although there is significant genetic structure of schistosomes among geographic regions and hosts, extensive human host movement in the area, leads to gene flow among infrapopulations. Hence the human host in this area serves both as a means of dispersal and reservoir of genotypes.
The high levels of genetic diversity reported here are comparable with those that have been detected in other mammalian hosts such as in the rat population in Guadeloupe and also those reported in an earlier study of the genetic diversity of S. mansoni within Kenya (Barral et al., 1996; Agola et al., 2006). Despite the homogeneity of the genotypes observed in this area, with most alleles being shared among the infrapopulations, private alleles (alleles occurring exclusively in one population) were observed in most of the infrapopulations analysed except four infrapopulations (NYD, NYF, NGA and KR infrapopulations). Interestingly, these private alleles occurred at frequencies higher than 5%. We hypothesize that the presence of private alleles in these populations could be as a result of exposure to unique schistosome genotypes found exclusively at one infection foci. An alternative hypothesis could be that since we are samplings offspring and not the adult population, private alleles could arise as a result of few highly fecund worm pairs with unique alleles that are producing more offspring thus creating the high frequencies of private alleles. In all the infrapopulations analysed in the present study, the observed heterozygosity was much lower than expected, if the populations were in Hardy-Weinberg equilibrium. These results are comparable to those of Rodrigues et al. (2002) and Stohler et al. (2004) who detected heterozygote deficiency in field populations of S. mansoni. We hypothesize here that the Wahlund effect (sampling different sites that are combined to represent a single population) is the most probable cause of heterozygote deficiency in the populations analysed in this study. This again supports the idea that there is a lot of movement and sampling of the parasite by school children from several infection foci within the larger Mwea area. Alternatively, the possible cause of heterozygote deficiency could have been due to the presence of null alleles and/or scoring errors. The frequencies of null alleles in all the populations analysed, however, were found to be low (0.5–5.6%) and therefore unlikely to be the cause of the observed heterozygote deficiency in the S. mansoni infrapopulations in Mwea. Long allele drop out is a problem in scoring microsatellites particularly when the quantity of DNA is low, which is an inherent problem with genotyping miracidia but this type of error appeared to be minimal in this study because the analysis with MICROCHECKER did not reveal any problematic loci.
In terms of the potential for the evolution of drug resistance, high genetic diversity suggests that there may be a greater potential for such populations to harbor alleles that may confer drug resistance, and also that the effective population size is large, which would increase the probability of novel mutations that confer resistance. Large population sizes tend to prevent the fixation of such alleles unless they have a strong selective advantage, or population sizes become dramatically reduced, as predicted with large scale drug treatment. Because there is gene flow occurring among subpopulations (infrapopulations) in Mwea, there is an opportunity for the spread of rare alleles, some of which may confer traits such as drug resistance or virulence. It is also important to note that in Mwea there is an ongoing schistosomiasis control program involving mass treatment of school children with praziquantel. Finally, although we do not know if drug resistance alleles occur in this population, one of the important factors that determine how rapidly such alleles would increase in frequency is the size of the treated proportion of the population relative to the untreated proportion, termed “refugium” ( Sissay et al., 2006; Leathwick et al., 2008; Waghorn et al., 2008). Our findings suggest this refugium is relatively large because of the high rates of gene flow and lack of population substructure throughout Mwea and only a subset of school children are targeted for treatment.
We thank the children and teachers of Nyamindi, Kirogo, Ngurubani and Mukou primary schools in Mwea for their participation in the study and cooperation, the staff of the Division of Vector Borne Disease (DVBD), Mwea in particular, Mr. Njoroge for facilitating access to their lab facility during field sample collection, and Dr. Charles Mwandawiro, Jimmy Kihara and the staff of the Japan International Co-operation Agency (JICA)-supported Eastern and Southern Africa Centre of International Parasite Control (ECACIPAC) for allowing us access to their field sites and the children participating in their schistosomiasis control program. We are equally grateful to Mr. George Rosenberg and the staff of the Molecular Biology Facility in the Department of Biology, University of New Mexico (UNM) for help with genotyping. The technical assistance of Joseph Kinuthia, Geoffrey Maina and Stephen Kamau of the Centre for Biotechnology Research and Development of KEMRI is highly appreciated. This research was supported through NIH grant R01-A144913. Approvals for this study were obtained through KEMRI’s scientific and ethical review committees. This paper is being published with the approval of the Director, KEMRI.
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
L. E. Agola, Centre for Biotechnology Research and Development, Kenya Medical Research Institute (KEMRI), P.O. Box 54840, Nairobi, Kenya.
M. L. Steinauer, Department of Biology, University of New Mexico, Albuquerque, NM 87131, USA.
D. N. Mburu, University of Manitoba/ Nairobi project, P.O. Box 19676, Nairobi, Kenya.
B.N. Mungai, Centre for Biotechnology Research and Development, Kenya Medical Research Institute (KEMRI), P.O. Box 54840, Nairobi, Kenya.
I.N. Mwangi, Centre for Biotechnology Research and Development, Kenya Medical Research Institute (KEMRI), P.O. Box 54840, Nairobi, Kenya.
G.N. Magoma, Department of Biochemistry, Jomo Kenyatta University of Agriculture and Technology (JKUAT), P.O. Box 43844, Nairobi, Kenya.
E. S. Loker, Department of Biology, University of New Mexico, Albuquerque, NM 87131, USA.
G. M. Mkoji, Centre for Biotechnology Research and Development, Kenya Medical Research Institute (KEMRI), P.O. Box 54840, Nairobi, Kenya.