Understanding the evolution of parasites is important to both basic and applied evolutionary biology. Knowledge of the genetic structure of parasite populations is critical for our ability to predict how an infection can spread through a host population and for the design of effective control methods. However, very little is known about the genetic structure of most human parasites, including the human louse (Pediculus humanus). This species is composed of two ecotypes: the head louse (Pediculus humanus capitis De Geer), and the clothing (body) louse (Pediculus humanus humanus Linnaeus). Hundreds of millions of head louse infestations affect children every year, and this number is on the rise, in part because of increased resistance to insecticides. Clothing lice affect mostly homeless and refugee-camp populations and although they are less prevalent than head lice, the medical consequences are more severe because they vector deadly bacterial pathogens. In this study we present the first assessment of the genetic structure of human louse populations by analyzing the nuclear genetic variation at 15 newly developed microsatellite loci in 93 human lice from 11 sites in four world regions. Both ecotypes showed heterozygote deficits relative to Hardy–Weinberg equilibrium and high inbreeding values, an expected pattern given their parasitic life history. Bayesian clustering analyses assigned lice to four distinct genetic clusters that were geographically structured. The low levels of gene flow among louse populations suggested that the evolution of insecticide resistance in lice would most likely be affected by local selection pressures, underscoring the importance of tailoring control strategies to population-specific genetic makeup and evolutionary history. Our panel of microsatellite markers provides powerful data to investigate not only ecological and evolutionary processes in lice, but also those in their human hosts because of the long-term coevolutionary association between lice and humans.
Citation: Ascunce MS, Toups MA, Kassu G, Fane J, Scholl K, Reed DL (2013) Nuclear Genetic Diversity in Human Lice (Pediculus humanus) Reveals Continental Differences and High Inbreeding among Worldwide Populations. PLoS ONE 8(2): e57619. https://doi.org/10.1371/journal.pone.0057619
Editor: Colin Dale, University of Utah, United States of America
Received: September 13, 2012; Accepted: January 24, 2013; Published: February 27, 2013
Copyright: © 2013 Ascunce et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported in part by grants to DLR from the National Science Foundation (www.nsf.gov: DEB 0717165 and DEB 0845392) and the University of Florida (www.research.ufl.edu: SEED Fund Program). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
The study of genetic diversity within and among parasite populations can provide knowledge of parasite evolutionary history, identify genes under selection, and elucidate the origin and spread of disease. Because parasites often show long-term associations with their hosts, parasites can also be used to infer host evolutionary history. In fact, several human parasites have provided insight into aspects of human evolution that were unclear from studies of direct human evidence, such as fossil or molecular data –. Among them, the human louse (Pediculus humanus) is thought to be an ancient parasite based on archeological remains, the worldwide co-distribution with humans, and coevolutionary studies –. This single species of louse is composed of two distinct morphological, behavioral, and ecological types: the head louse (Pediculus humanus capitis De Geer), and the clothing (body) louse (Pediculus humanus humanus Linnaeus) (see  for a review). Human lice are blood-sucking, wingless, host-specific ectoparasites of humans that are both obligate (cannot live off the host) and permanent (complete their life cycle on a single host species). Thus, these parasites are inextricably tied to their host in ecological and evolutionary time. Moreover, the coevolution of lice and their primate hosts over the last 25-million-year (MY) is well documented , , . Although we have learned about the long-term evolutionary history of human lice, we know less about the genetic structure of living populations. The few studies that have looked at this were based on mitochondrial or a small number of nuclear DNA sequences, and those studies provided limited knowledge of the genetic structure of human louse populations worldwide. Understanding the geographic distribution of louse diversity worldwide is important for a variety of reasons, and is the focus of this study.
There are hundreds of millions of head louse infestations every year affecting mostly children of 3 to14 years of age . During recent years head louse infestations have increased globally, in part because of increased resistance to insecticidal shampoos . Most prominent is resistance to synthetic pyrethroids such as phenothrin and permethrin, which are the most widely used insecticides for louse control. Pyrethroid resistance in human lice is due to three knockdown (kdr)-type point mutations (M815I, T917I, and L920) in the voltage-sensitive sodium channel alpha-subunit gene , . The kdr-type resistance has been reported in many countries among head louse populations , , and recently in clothing lice from France . In order to better understand how these resistance alleles are spread, we must also have a firm understanding of neutral genetic variation and patterns of gene flow among populations of lice. In that regards, the mechanisms of how these resistance alleles spread will be improved by a better understanding of global genetic structure, which will be best revealed using highly polymorphic genetic markers like microsatellites.
Clothing lice affect predominantly homeless and refugee-camp populations – and are less prevalent than head lice but far more serious because they vector at least three deadly bacterial pathogens, those responsible for epidemic typhus (Rickettsia prowazekii), trench fever (Bartonella quintana), and relapsing fever (Borrelia recurrentis). In the last 30 years, several outbreaks of louse-borne diseases have occurred, such as epidemic typhus in Burundi, which infected over 45,000 people . It once was believed that only clothing lice vectored these bacterial pathogens, however head lice have been found to carry Bartonella quintana in Nepal , the United States , and Ethiopia , . Moreover, experimental infections showed that head lice can vector louse-borne diseases , . Therefore, understanding the genetic structure of both head lice and clothing lice worldwide is of critical importance to our understanding of the risk of epidemic disease.
The genetic diversity of human lice has been widely studied using the mitochondrial (mt) cytochrome c oxidase 1 gene (COX1) revealing the presence of three deeply divergent mtDNA clades or haplogroups, named A, B, and C , ,  (Figure 1). Haplogroup A is the most common and has a worldwide distribution, whereas B and C are geographically restricted , . Haplogroup B is found in the New World, Europe and Australia, whereas haplogroup C has only been found in Nepal and Ethiopia. Multi-Spacer-Typing (MST) markers, which included a set of four intergenic spacers, showed that lice from Clade A formed two geographic clusters, one containing all the Clade A lice outside Africa, and the other including African lice , . Although the MST technique proved useful for differentiating African and non-African haplotype A lice, further refinement of global louse genetic structure requires markers with greater resolution. To that end we have mined the Pediculus humanus genome to identify new target DNA sequences to develop as microsatellite markers. We genotyped 93 human lice from 11 geographic sites distributed throughout the globe using 15 newly developed microsatellite loci. This panel of microsatellites was successful in uncovering strong signals of genetic structure that corresponded to geography and ecotype.
Materials and Methods
The Institutional Review Board of the University of Florida exempted the study from review (Exemption of Protocol #2009-U-0422) and waived the need for written informed consent of the participants. This exemption is issued based on the United States Department of Health and Human Services (HHS) regulations. Specifically, HHS regulation 45 CFR part 46 applies to research activities involving human subjects. Because louse removal was voluntary and no information was recorded that would allow patients to be identified directly or through identifiers linked to them, the University of Florida waived the need for written informed consent from the participants.
Sampling details are provided in Methods S1. Briefly, 75 human head lice were collected from different individuals at 10 localities throughout the world (Table 1). Clothing lice came from two sites: Canada and Nepal. Canadian clothing lice (N = 16) were collected from a single homeless person. The two clothing lice from Nepal (N = 2) were collected from two persons.
Data harvesting, screen for repeat motifs and primer design
An overview of our primer development strategy and multiplex optimization is presented in Figure 2 and we included a detailed description for microsatellite development in Methods S1. Briefly, assembled genome sequence data for Pediculus humanus USDA strain (PhumU1, 2007) were obtained from VectorBase (http://www.vectorbase.org/) ,  and screened for tandem repeat motifs. Candidate sequences were subject to primer design using Primer3Plus .
Multiple co-amplification: Multiplexes
Primer pairs that produced amplicons of the expected size and demonstrated allelic variation were selected for further testing and multiplex optimization (Table S1). We used fluorescently labeled universal primers: M13 (5’-CACGACGTTGTAAAACGAC-3’) and CAG (5’-CAGTCGGGCGTCATCA-3’), each labeled with a unique fluorescent tag (e.g. FAM, VIC, NED, PET) to co-amplify multiple loci. Locus-specific primers were modified by adding the matching 5′ universal primer sequence tails.
Microsatellite diversity analysis
The final optimized four multiplexes were tested on the DNA extracted from 93 human lice from 11 sites (Table 1; Methods S1). Evidence for large allelic drop out and null alleles was determined using Microchecker version 2.2.3 . We used the software Arlequin version 22.214.171.124  to determine number of alleles per locus, observed and expected heterozygosity (HO and HE), and mean FIS estimates (an index of the inbreeding of individuals resulting from the non-random union of gametes within a subpopulation) over loci per populations. Confidence intervals for FIS estimates were calculated by bootstrapping over loci using 1,023 random permutations. Genotypic disequilibrium among loci was estimated using FSTAT version 126.96.36.199 , .
Population structure was inferred with a Bayesian clustering approach implemented in the STRUCTURE software  (http://pritch.bsd.uchicago.edu/structure.html). This method uses the individual multi-locus genotypic data to evaluate models assuming different numbers of genetic clusters (K) based on the posterior probabilities given the data and model. All simulations used 50,000 Markov chain Monte Carlo (MCMC) generations in the burn-in phase and 100,000 generations in the data collection phase. Ten independent runs, using default parameters for each K, to ensure equilibration during burn-in and consistency in estimation of the posterior probabilities. Selection of the number of distinct clusters was based on the evaluation of the ΔK statistic . The ten STRUCTURE runs at each K produced nearly identical individual membership coefficients. The run with the highest likelihood of the data given the parameter values for the predominant clustering pattern (i.e. the mode) at each K was used for plotting with DISTRUCT  (http://www.stanford.edu/group/rosenberglab/distruct.html). A series of STRUCTURE analyses were conducted: S1) the worldwide dataset for all 15 loci, S2) the worldwide dataset for 14 loci (we excluded locus M2-13 due to some missing data), S3) the worldwide dataset for 14 loci plus the mitochodrial haplogroup coded as another marker, and S4) local datasets for sites with 10 or more lice (Canada, New York, Honduras and Cambodia). We ran up to K = 15 for the global (S1, S2 and S3), and K = 5 for the local (S4) datasets. The mitochondrial haplogroup information was obtained from Ascunce et al. (in prep.).
Pairwise estimates of population structure and gene flow
Principal coordinate analysis (PCA) was conducted using pairwise genetic difference between individuals calculated in GenAlEx ,  (http://biology.anu.edu.au/GenAlEx) to validate and further define genetic clusters for these lice. For sample sizes exceeding 10 lice, we also estimated population pairwise values of the Weir and Cockerham (1984)  analogue of FST, ΘWC, in FSTAT 188.8.131.52 ,  and gene flow (Nem) based on the private alleles method , as implemented in the online program Genepop version 4.0.10 ,  (http://genepop.curtin.edu.au/). These two estimates: ΘWC and Nem were also evaluated considering the genetic clusters inferred from STRUCTURE simulations.
To detect outlier loci under selection we used the program BayeScan  (http://cmpg.unibe.ch/software/bayescan/). BayeScan is a hierarchical Bayesian method that assumes that allele frequencies within populations follow a multinomial-Dirichlet distribution –. It estimates population-specific FST coefficients, therefore allowing for different demographic histories and different amounts of genetic drift between populations. BayeScan incorporates the uncertainty on allele frequencies due to small sample sizes. The estimation of model parameters was automatically tuned on the basis of 20 short pilot runs of 5,000 iterations. The sample size was set to 5,000 and the thinning interval to 10, resulting in a total chain length of 100,000 iterations. Four independent runs were performed for each of the two datasets to account for the consistency of the detected outliers. The loci were ranked according to their estimated posterior probability and all loci with a value over 0.76 were retained as outliers. This corresponds to a Bayes Factor >3, which provides substantial support for acceptance of the model .
The average relatedness (r) among lice was determined using the software Relatedness version 5.0  (http://www.gsoftnet.us/GSoft.html) for the following groups of lice: R1) clothing lice from a single homeless person in Canada; R2) head lice from people in New York; and a single head louse from each child in one of two orphanages in: R3) Cambodia, and R4) Honduras.
We found a total of 295,733 di-, tri- and tetranucleotide tandem repeat motifs sequences. Tri-nucleotide motifs are the most abundant microsatellites, making up 62% of perfect microsatellites (Figure S1). Comparative genomic studies have revealed a great heterogeneity in microsatellite abundance and composition across taxa. Particularly, some other arthropods such as Aedes aegypti  have also shown skews toward tri-nucleotide repeats.
From the approximately 150 primer pairs tested, 33% gave clear amplifications and showed allelic diversity. A final set of 15 microsatellite loci were thoroughly optimized and validated over 93 human lice from around the world (Tables S1, S2, and S3). All microsatellite loci were highly polymorphic, with an average of 10 alleles per locus and an average HO and HE of 0.2748 and 0.7136, respectively (Table 2). In populations with more than 10 lice (Canada, New York, Honduras and Cambodia), Micro-Checker analysis found no evidence for scoring errors due to stuttering or large-allele dropout for each locus in each population. However, Micro-Checker identified null alleles in multiple loci at each geographic location. In these same populations, genotype proportions deviated significantly from Hardy-Weinberg expectations due to heterozygote deficiencies (Table S4). Values of FIS were significantly positive for a large number of loci, indicating also heterozygote deficiencies within these populations (Figure 3). Overall all four sites showed high values of FIS from 0.232 in Canada to 0.767 in New York (Table 2). Moreover, HO and FIS were significantly correlated within New York and Cambodia populations (Spearman rank correlation tests; P<0.001). Out of 519 tests for linkage disequilibrium, none of the tests were found to be significant after Bonferroni correction (adjusted at the 5% level; P<0.000043). We excluded M2-13 loci from further analyses because of its relatively high percentage of missing data (26%).
For monomorphic loci estimates were not calculated and they are shown in the graphs as gaps. Stars indicate significant departure from Hardy-Weinberg equilibrium (P<0.05).
Nuclear genetic clusters at the worldwide scale
Population structure was inferred with a Bayesian clustering approach implemented in the STRUCTURE software . In all the three STRUCTURE analysis, all worldwide human lice were assigned to four genetic clusters (K = 4), one defined by clothing lice from Canada, the other head lice from North America and Europe, a third cluster was composed of head lice from Honduras, and the fourth cluster included Asian lice (both head and clothing lice) (Figure 4A) (STRUCTURE with 15 microsatellite loci is not shown). In another STRUCTURE analysis, we incorporated the mitochondrial haplogroup data for each louse (Figure 4B). This STRUCTURE analysis showed that there is no correlation between mitochondrial haplogroups (A and B) and nuclear genetic clusters, at least among the current samples. We also employed the multivariate technique Principal Coordinate Analysis (PCA) implemented in GenAlEx ,  that allows the visualization of the spatial distribution of the genetic differences among all samples. We found similar results as the STRUCTURE analyses, where one cluster included head lice from North America and Europe, while all Asian and Central American lice comprised a second cluster with the exception of the clothing louse from Nepal that showed an intermediate position between this group and clothing lice from Canada (Figure 4).
A) the worldwide dataset for 14 loci (we excluded M2-13 due to missing data), and B) the worldwide dataset for all 14 loci plus the mitochondrial haplogroup coded as an additional locus. In the bar plot, each individual is represented by a single vertical line and the length of each color segment represents the proportion of membership (Q) to the four clusters. In pannel B, for each louse sample we added the mitochondrial haplogroup as filled black squares for Clade A and the open diamonds representing Clade B. Distributions of points in the first two dimensions resulting from principal coordination analyses (PCA) conducted using pairwise genetic distance comparisons of the same dataset used for the STRUCTURE analyses are below each STRUCTURE plot.
Nuclear genetic clusters at the local/cluster scale
For Canada, New York, Honduras, and Cambodia populations we further analyzed their genetic substructure by analyzing each population individually using STRUCTURE. These results revealed an increase in the number of regional genetic clusters in New York (K = 3), and in Cambodia (K = 3) (Figure 5). Although Evanno's method cannot evaluate K = 1 as the most likely number of clusters, we found that populations from Canada and Honduras showed admixture for all individuals when K = 2. This can be interpreted as evidence supporting Canada and Honduras as a single genetic cluster, respectively, at least with the current number of microsatellite markers analyzed.
Pairwise estimates of population structure and gene flow
In order to estimate gene flow, we calculated FST and Nem across human louse populations with more than 10 lice (Table S5). Both FST and Nem produced results concordant with the STRUCTURE analysis. There were high values of genetic differentiation (FST, above 0.4087) between the clothing louse population (Canada) and any of the remaining head louse populations. However, these results should be taken with caution and more clothing louse populations should be analyzed to corroborate these results. Effective rates of nuclear gene flow (Nem) using the private allele method yielded values less than 1.0. The highest value was found between Central America and Asia (0.7961), which may represent a more recent ancestry.
The microsatellite data set was used in a global analysis for outlier detection using BayeScan . A single locus (T2_7) showed evidence of balancing selection with a posterior probability of 0.9962, log10 (PO) of 2.4185. This corresponded to a Bayes Factor larger than 100, which provides ‘decisive’ evidence for selection .
For the fine scale genetic structure analyses we examined the patterns of relatedness (r) among clothing lice from a single homeless person, and head lice from each population. These analyses can determine whether a few genotypes are responsible for the production of the majority of offspring. Although average population relatedness values (r) were negative suggesting outbreeding and unrelated populations  (Table 3), a detailed analysis of pairwise comparison between individuals within each site showed a variable range of relatedness (Figure S2). Roughly, 50% of lice in each site had r values lower than 0, thus they are very different from each other. The second largest group of lice had r values up to 0.4, showing some relatedness. All sites showed the presence of some full-siblings (r = 0.5) and the highest level of inbreeding with identical genotypic profiles was found in New York (r = 1). These high relatedness values found in the New York samples parallel the evidence of heterozygote deficits and high inbreeding, reflecting a biological phenomenon.
Life history traits and genetic diversity
Understanding the processes that shape the genetic structure of parasite populations is critical in predicting how an infection can spread through a host population and for the design of effective control methods. Here, we provide the first assessment of genetic structure in human louse populations from around the world using microsatellite loci. It is known that parasite life history strategies, such as dispersal and mode of transmission, directly affect their genetic structure. In our study, a large number of loci showed significant heterozygote deficits relative to Hardy–Weinberg equilibrium plus high FIS-values (Table 2, Table S4). This genome-wide pattern would suggest that inbreeding is common in human lice, a pattern expected given its parasitic life history. In an earlier study using five microsatellite loci analyzing head and clothing lice from doubly-infested individuals , the authors also found that certain loci were consistently out of Hardy-Weinberg equilibrium reflecting a population-specific phenomenon rather than microsatellite locus-specific issue. Alternatively, deficiency of heterozygotes could also been explained by null alleles, however we successfully mapped all PCR primers used in this study to assembled genomes from clothing lice, mtDNA Clade A head lice, and mtDNA Clade B head lice, and found no nucleotide mismatches at our priming sites from the clothing louse genome and only one mismatch in each of three primers from the head louse genomes. The Wahlund effect (population substructure) could also cause significant heterozygote deficits relative to Hardy–Weinberg equilibrium. Indeed, in New York three genetic clusters were detected through STRUCTURE analysis.
Interestingly, we found a weaker genome-wide effect in terms of heterozygote deficits relative to Hardy–Weinberg equilibrium and null alleles in clothing lice than among head lice (Table 2, Table S4). Some life history traits of clothing lice may lead to higher effective population sizes (Ne) than in head louse populations. For example, clothing louse females lay up to 300 eggs along the seams or hems of clothes compared with 150 eggs laid (and glued) on human hair shafts by head louse females.
Another factor that can contribute to the difference in heterozygosity is selection for insecticide resistance. Among head louse populations, intense selection by insecticides may have resulted in periodic population bottlenecks. These demographic events are expected to have genome-wide effects by reducing genetic polymorphisms in the louse genome. Whereas insecticide resistance is well known for head lice, it is less common in clothing lice and only recently has a study performed kdr-allele typing for clothing lice from France . Further studies including sympatric clothing and head louse populations are needed to discern the role of these factors in shaping the genetic structure and diversity of human lice.
Gene flow and the evolution of insecticide resistance
Gene flow can have broad implications in the evolution of drug resistance. For example, the widespread resistance to the most common antimalarial drugs (chloroquine and sulfadoxine-pyrimethamine) is thought to be the result of few drug-resistance alleles coupled with high gene flow in vectors (mosquitoes) and human hosts –. Alternatively, convergent insecticide resistance phenotypes can evolve through identical or different mutations in independent populations. For example, the same mutations that confer insecticide resistance in the E3 esterases of blowflies and houseflies can arise de novo and be selected for in separate populations . In human lice, because of the strong genetic structure with low gene flow among populations, it seems more likely that human lice will follow a model of parallel adaptive evolution, where new resistance alleles evolve in different populations. However, the study of the evolution of resistance in human lice is in its infancy and many basic questions still need to be addressed. Pyrethroid resistance occurs worldwide with variable resistant kdr-like haplotype frequencies in each geographic area , . Moreover, a recent study highlighted that some kdr-alleles were not correlated with treatment failure, thus other factors may be involved . In other insects, functional genomic analysis has shown that insecticide resistance could constitute a multigenic trait involving large parts of the insects' genome –. To test the model of resistance evolution in human lice, regional comprehensive studies are needed combining phenotypic and genotypic analysis of resistance coupled with neutral markers, including microsatellite loci linked to the voltage-sensitive sodium channel alpha-subunit gene as well as other genes involved in resistance. Interestingly, with our microsatellite markers panel, we detected one outlier locus that could be mapped at a unique location in the reference genome and was in the vicinity of a putative gene coding for a carbonic anhydrase enzyme. These enzymes seem to be involved in pH regulation (alkalization mechanisms) of a lepidopteran caterpillar and larval mosquito guts –. Thus, we suggest that this gene is a promising candidate for future functional analysis.
Clothing and head louse differentiation
The species status of P. humanus (whether head and clothing lice represent one species or two) has been a topic for debate for over a century (see  for a review). However, the difficulty in determining whether a group of organisms constitutes an independently evolving lineage is particularly compounded in parasites. In our study, the clothing lice (Canada and Nepal) grouped more closely with the Central America-Asia cluster probably because of close ancestry (Figure 4). These results are consistent with the idea that clothing lice evolved from head louse ancestors, invading the body region only recently with the advent of clothing use in modern humans , , . Further, studies have shown that clothing lice emerged from only one of the three mitochondrial haplogroups (Clade A) roughly 83,000 years ago . Although clothing lice belong to a single mtDNA clade, they appear to have evolved locally (in situ) throughout the world from head louse populations . In contrast, Leo et al. (2005) used five microsatellite loci to address the issue of species status within P. humanus by analyzing head and clothing lice from doubly-infested individuals . The authors found that head and clothing lice from each human host formed separate nuclear clusters, and each host represented a unique genetic cluster . Their conclusion should be taken with caution due to the small sample size, restricted geographic distribution of the samples and few loci used in their study. In our study, we are also puzzled by the distant genetic relationship between Canada (clothing) and New York (head) lice. Could the strong pyrethroid-based pediculicide treatment of head lice in developed countries  account for the large differentiation found in this study between clothing lice from Canada and head lice from New York? Could that be the case for the louse population from China and Nepal studied by Leo et al. (2005) ? Understanding the presence or absence of gene flow between head and clothing lice is of paramount importance, because in addition to clothing lice, researchers are now finding an increasing number of cases of louse-borne bacteria in head lice.
Human louse genetic diversity and human evolution
One of the more exciting aspects of our understanding of nuclear diversity in lice is its application to human evolution. Bayesian clustering analyses assigned lice to four distinct genetic clusters: the first cluster consisting of clothing lice from Canada, the other cluster included head lice from North America and Europe, a third cluster was composed of head lice from Honduras, while the fourth cluster included Asian lice. How this geographic structure reflects human migrations requires greater sampling. Although preliminary, our study suggests that the Central America-Asian cluster is mirroring the (human host) colonization of the New World if Central American lice were of Native American origin and Asia was the source population for the first people of the Americas as has been suggested (Figure 6) , . The USA head louse population might be of European decent, explaining its clustering with lice from Europe. Within the New World, the major difference between USA and Honduras may reflect the history of the two major human settlements of the New World: the first peopling of America and the European colonization after Columbus (Figure 6).
Colored circles on map indicate collecting sites, with the color of each circle corresponding to the majority nuclear genetic cluster to which sampled individuals were assigned. Large colored circles are sites with 16 or more lice, small colored circles represent sites with one to three lice. Thick grey arrows indicate proposed migrations of anatomically modern humans out of Africa into Europe, Asia and the Americas, as well as the more recent European colonization of the New World. Colored arrows represent hypothetical human louse co-migrations. The bottom panel is the plot from STRUCTURE corresponding to the assignment of 93 lice from 11 geographical sites (from Figure 4A).
This study provides preliminary evidence that microsatellites are effective in determining population genetic structure among human louse populations at a scale that could help us better understand patterns of human migration worldwide and might also provide insight into interactions between archaic hominids and anatomically modern humans, Homo sapiens. For example, previous studies  suggested that louse mtDNA haplogroup A has had a long history associated with the host lineage that led to anatomically modern humans, Homo sapiens (Figure 1). Studies of modern human expansion out of Africa show the footprint of serial founder effects on the genetic diversity of human populations as revealed by the human pattern of increased genetic distance and decreased diversity with distance from Africa . The microsatellite loci developed in this study are ideal markers to measure louse genetic diversity and how it parallels to human diversity. Louse mitochondrial haplogroup B is found in the New World, Europe and Australia but not in Africa. Reed et al.  suggested that its evolutionary origins might lie with archaic hominids from Eurasia (i.e., Homo neanderthalensis) and that they became associated with modern humans via a host switch during periods of overlap. If true, then examining nuclear markers for lice with haplotypes A (host: H. sapiens) and B (past host: H. neanderthalensis; current host: H. sapiens) would permit a test for admixture in lice and could provide a time frame for when louse admixture occurred (i.e., host switching). This host switching would give us a time period during which H. sapiens and H. neanderthalensis co-occurred. Whether modern humans interbred with archaic hominids is a much-debated question –. However, genomic evidence using ancient DNA from Neanderthals and modern DNA from living humans suggests that interbreeding might have occurred between 47,000 and 65,000 years ago . A less clear question is whether anatomically modern humans overlapped or interbred with other archaic hominids in Asia or Africa. Recent genetic studies have also supported both Asian and African archaic admixtures. If the older haplogroup C lice evolved on archaic hominins in Asia or Africa, then the study of haplogroup C lice could also provide compelling evidence of close-proximity interactions of modern and archaic hominins in Asia or Africa.
Understanding the processes that shape the genetic structure of human parasite populations is important to both basic and applied evolutionary biology. In human lice, knowing the processes and mechanisms that maintain and generate human louse genetic diversity within and among host populations is critical for our ability to design effective control methods and to predict how louse-borne diseases can spread through human populations. In this study, we showed that human louse populations are genetically structured based on geography (Figure 6). The high degree of genetic structure may lead to the evolution of different resistance alleles among different populations, thus suggesting the need of regional epidemiological studies to control lice. Furthermore, this work has shown that the study of the genetic diversity of human lice could help us better understand patterns of human migrations worldwide at different temporal scales, and could also be used to test hypotheses about human evolution such as ecological interactions between modern and archaic hominins.
Microsatellite abundance (counts): Data are shown for di- (diagonal lines), tri- (grey), and tetra- (solid black) per number of repeat motifs.
Average pairwise relatedness (r) of lice within sites. Relatedness values can range from 1 (individuals are identical for all alleles assessed) to –1 (individuals have no alleles in common).
This document includes: Supporting Information Methods and Supporting Information References.
Characteristics of 27 microsatellite loci developed using the genome data of clothing louse (PhumU1, 2007).
Master mixes of primers for each multiplex. All primers working solutions are 10 µM.
Total PCR mix, quantities are indicated for individual tubes and for whole plates. Labeled tail primer concentration is 10 µM.
We greatly appreciate the contribution of Katie Shepherd (Founder & CEO of The Shepherd Institute for Lice Solutions), Kazunori Yoshizawa (Hokkaido University, Japan), and Douglas D. Colwell (Lethbridge Research Centre, Canada) and collectors worldwide for providing lice used in this study. We would like to thank Jim Austin (University of Florida) and Drew Kitchen (University of Iowa) for comments on earlier version of this manuscript.
Conceived and designed the experiments: MSA MAT DLR. Performed the experiments: MSA MAT GK JF. Analyzed the data: MSA. Contributed reagents/materials/analysis tools: KS DLR. Wrote the paper: MSA DLR.
- 1. Reed DL, Toups MA, Light JE, Allen JM, Flanagin S (2009) Lice and other parasites as markers of primate evolutionary history. In: Huffman MA, Chapman CA, editors. Primate Parasite Ecology: The Dynamics and Study of Host-Parasite Relationships. Cambridge: Cambridge University Press. p. 531.
- 2. Nozais J-P (2003) The origin and dispersion of human parasitic diseases in the old world (Africa, Europe and Madagascar). Mem Inst Oswaldo Cruz 98 Suppl 113–19 Available: http://www.ncbi.nlm.nih.gov/pubmed/12687757.
- 3. Dittmar K, Araújo A, Reinhard KJ (2012) The Study of Parasites Through Time: Archaeoparasitology and Paleoparasitology. In: Grauer AL, editor. A companion to paleopathology. Chichester, West Sussex. Malden, MA: Wiley-Blackwell.
- 4. Araújo A, Ferreira LF, Guidon N, Maues da Serra Freire N, Reinhard KJ, et al. (2000) Ten thousand years of head lice infection. Parasitol Today 16: 269 Available: http://www.sciencedirect.com/science/article/pii/S016947580001694X.
- 5. Burgess IF (1995) Human lice and their management. Adv Parasitol 36: 271–342 Available: http://www.ncbi.nlm.nih.gov/pubmed/7484466.
- 6. Burgess IF (2004) Human lice and their control. Annu Rev Entomol 49: 457–481 Available: http://www.ncbi.nlm.nih.gov/pubmed/14651472.
- 7. Reed DL, Smith VS, Hammond SL, Rogers AR, Clayton DH (2004) Genetic analysis of lice supports direct contact between modern and archaic humans. PLoS Biol 2: e340 Available: http://www.ncbi.nlm.nih.gov/pubmed/15502871.
- 8. Reed DL, Light JE, Allen JM, Kirchman JJ (2007) Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biology 5: 7 Available: http://www.biomedcentral.com/1741-7007/5/7.
- 9. Light JE, Toups MA, Reed DL (2008) What's in a name: The taxonomic status of human head and body lice. Mol Phylogenet Evol 47: 1203–1216 Available: http://www.sciencedirect.com/science/article/pii/S1055790308001115.
- 10. Raoult D, Reed DL, Dittmar K, Kirchman JJ, Rolain J, et al. (2008) Molecular Identification of Lice from Pre-Columbian Mummies. J Infect Dis 197: 535–543 Available: http://jid.oxfordjournals.org/content/197/4/535.full.
- 11. Gratz NG (1997) Human lice: their prevalence, control and resistance to insecticides, A review 1985-1997. Geneva, Switzerland: WHOPES, CTD, WHO.
- 12. Lee SH, Yoon K-S, Williamson MS, Goodson SJ, Takano-Lee M, et al. (2000) Molecular analysis of kdr-like resistance in permethrin-resistant strains of head lice, Pediculus capitis. Pestic Biochem Physiol 66: 130–143 Available: http://www.sciencedirect.com/science/article/pii/S0048357599924604.
- 13. Lee SH, Clark JM, Ahn YJ, Lee W-J, Yoon KS, et al. (2010) Molecular mechanisms and monitoring of permethrin resistance in human head lice. Pestic Biochem Physiol 97: 109–114 Available: http://www.sciencedirect.com/science/article/pii/S0048357509000571.
- 14. Clark JM (2009) Determination, mechanism and monitoring of knockdown resistance in permethrin-resistant human head lice, Pediculus humanus capitis. J Asia Pac Entomol 12: 1–7 Available: http://www.ncbi.nlm.nih.gov/pubmed/20161186.
- 15. Hodgdon HE, Yoon KS, Previte DJ, Kim HJ, Aboelghar GE, et al. (2010) Determination of knockdown resistance allele frequencies in global human head louse populations using the serial invasive signal amplification reaction. Pest Manag Sci 66: 1031–1040 Available: http://www.ncbi.nlm.nih.gov/pubmed/20564731.
- 16. Drali R, Benkouiten S, Badiaga S, Bitam I, Rolain JM, et al. (2012) Detection of a knockdown resistance mutation associated with permethrin resistance in the body louse Pediculus humanus corporis by use of melting curve analysis genotyping. J Clin Microbiol 50: 2229–2233 Available: http://www.ncbi.nlm.nih.gov/pubmed/22573588.
- 17. Rydkina EB, Roux V, Gagua EM, Predtechenski AB, Tarasevich IV, et al. (1999) Bartonella quintana in body lice collected from homeless persons in Russia. Emerging Infect Dis 5: 176–178 Available: http://www.ncbi.nlm.nih.gov/pubmed/10081691.
- 18. Seki N, Sasaki T, Sawabe K, Sasaki T, Matsuoka M, et al. (2006) Epidemiological studies on Bartonella quintana infections among homeless people in Tokyo, Japan. Jpn J Infect Dis 59: 31–35 Available: http://www.ncbi.nlm.nih.gov/pubmed/16495631.
- 19. Foucault C, Brouqui P, Raoult D (2006) Bartonella quintana characteristics and clinical management. Emerg Infect Dis 12: 217–223.
- 20. Raoult D, Roux V, Ndihokubwayo JB, Bise G, Baudon D, et al. (1997) Jail fever (epidemic typhus) outbreak in Burundi. Emerging Infect Dis 3: 357–360 Available: http://www.ncbi.nlm.nih.gov/pubmed/9284381.
- 21. Sasaki T, Poudel SKS, Isawa H, Hayashi T, Seki N, et al. (2006) First molecular evidence of Bartonella quintana in Pediculus humanus capitis (Phthiraptera: Pediculidae), collected from Nepalese children. J Med Entomol 43: 110–112.
- 22. Bonilla DL, Kabeya H, Henn J, Kramer VL, Kosoy MY (2009) Bartonella quintana in body lice and head lice from homeless persons, San Francisco, California, USA. Emerging Infect Dis 15: 912–915 Available: http://www.ncbi.nlm.nih.gov/pubmed/19523290.
- 23. Angelakis E, Diatta G, Abdissa A, Trape J, Mediannikov O, et al. (2011) Altitude-dependent Bartonella quintana genotype C in head lice, Ethiopia. Emerging Infect Dis 17: 2357–2359 Available: http://www.ncbi.nlm.nih.gov/pubmed/22172306.
- 24. Cutler S, Abdissa A, Adamu H, Tolosa T, Gashaw A (2012) Bartonella quintana in Ethiopian lice. Comp Immunol Microbiol Infect Dis 35: 17–21 Available: http://www.ncbi.nlm.nih.gov/pubmed/22019400.
- 25. Goldberger J, Anderson JF (1912) The transmission of Typhus fever, with especial reference to transmission by the head louse (Pediculus capitis). Public Health Reports (1896-1970) 27: 297–307 .
- 26. Murray ES, Torrey SB (1975) Virulence of Rickettsia prowazeki for head lice*. Ann N Y Acad Sci 266: 25–34 .
- 27. Kittler R, Kayser M, Stoneking M (2003) Molecular evolution of Pediculus humanus and the origin of clothing. Curr Biol 13: 1414–1417 Available: http://www.ncbi.nlm.nih.gov/pubmed/12932325.
- 28. Kittler R, Kayser M, Stoneking M (2004) Molecular evolution of Pediculus humanus and the origin of clothing. Curr Biol 14: 2309 Available: http://www.sciencedirect.com/science/article/pii/S0960982204009856.
- 29. Light JE, Allen JM, Long LM, Carter TE, Barrow L, et al. (2008) Geographic distributions and origins of human head lice (Pediculus humanus capitis) based on mitochondrial data. J Parasitol 94: 1275–1281 Available: http://www.bioone.org/doi/abs/10.1645/GE-1618.1.
- 30. Veracx A, Boutellis A, Merhej V, Diatta G, Raoult D (2012) Evidence for an African cluster of human head and body lice with variable colors and interbreeding of lice between continents. PLoS ONE 7: e37804 Available: http://www.ncbi.nlm.nih.gov/pubmed/22662229.
- 31. Li W, Ortiz G, Fournier P, Gimenez G, Reed DL, et al. (2010) Genotyping of human lice suggests multiple emergencies of body lice from local head louse populations. PLoS Negl Trop Dis 4: e641 Available: http://www.ncbi.nlm.nih.gov/pubmed/20351779.
- 32. Lawson D, Arensburger P, Atkinson P, Besansky NJ, Bruggner RV, et al. (2009) VectorBase: a data resource for invertebrate vector genomics. Nucleic Acids Res 37: D583–587 Available: http://www.ncbi.nlm.nih.gov/pubmed/19028744.
- 33. Kirkness EF, Haas BJ, Sun W, Braig HR, Perotti MA, et al. (2010) Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle. Proc Natl Acad Sci USA 107 (27) : 12168–12173 Available: http://www.pnas.org/content/107/27/12168.long.
- 34. Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, et al. (2007) Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res 35: W71–74 Available: http://www.ncbi.nlm.nih.gov/pubmed/17485472.
- 35. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4: 535–538 Available: http://doi.wiley.com/10.1111/j.1471-8286.2004.00684.x.
- 36. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567 Available: http://www.ncbi.nlm.nih.gov/pubmed/21565059.
- 37. Goudet J (1995) FSTAT (Version 1.2): A computer program to calculate F-Statistics. J Hered 86: 485–486 Available: http://jhered.oxfordjournals.org/content/86/6/485.short.
- 38. Goudet J (2002) FSTAT, a program to estimate and test gene diversities and fixation indices Version 184.108.40.206. Available: http://www2.unil.ch/popgen/softwares/fstat.htm. (Updated from Goudet, 1995).
- 39. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959 Available: http://www.ncbi.nlm.nih.gov/pubmed/10835412.
- 40. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620 Available: http://www.ncbi.nlm.nih.gov/pubmed/15969739.
- 41. Rosenberg NA (2004) distruct: a program for the graphical display of population structure. Mol Ecol Notes 4: 137–138 Available: http://doi.wiley.com/10.1046/j.1471-8286.2003.00566.x.
- 42. Peakall R, Smouse PE (2006) GenAlEx 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6: 288–295 Available: http://www.blackwell-synergy.com/doi/abs/10.1111/j.1471-8286.2005.01155.x.
- 43. Peakall R, Smouse P (2012) GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research - an update. Bioinformatics. Available: http://www.ncbi.nlm.nih.gov/pubmed/22820204.
- 44. Weir BS, Cockerman CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38: 1358–1370.
- 45. Barton NH, Slatkin M (1986) A quasi-equilibrium theory of the distribution of rare alleles in a subdivided population. Heredity (Edinb) 56(3): 409–415 Available: http://www.ncbi.nlm.nih.gov/pubmed/3733460.
- 46. Raymond M, Rousset F (1995) GENEPOP (Version 1.2): Population Genetics Software for Exact Tests and Ecumenicism. J Hered 86: 248–249 Available: http://jhered.oxfordjournals.org/content/86/3/248.short.
- 47. Rousset F (2008) GENEPOP’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour 8: 103–106 Available: http://www.ncbi.nlm.nih.gov/pubmed/21585727.
- 48. Foll M, Gaggiotti O (2008) A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics 180: 977–993 Available: http://www.ncbi.nlm.nih.gov/pubmed/18780740.
- 49. Balding DJ, Nichols RA (1995) A method for quantifying differentiation between populations at multi-allelic loci and its implications for investigating identity and paternity. Genetica 96: 3–12 Available: http://www.ncbi.nlm.nih.gov/pubmed/7607457.
- 50. Rannala B, Hartigan JA (1996) Estimating gene flow in island populations. Genet Res 67: 147–158 Available: http://www.ncbi.nlm.nih.gov/pubmed/8801187.
- 51. Balding DJ (2003) Likelihood-based inference for genetic correlation coefficients. Theor Popul Biol 63: 221–230 Available: http://www.ncbi.nlm.nih.gov/pubmed/12689793.
- 52. Jeffreys H (1961) Theory of probability. Oxford: Oxford University Press.
- 53. Queller DC, Goodnight KF (1989) Estimating relatedness using genetic markers. Evolution 43: 258–275.
- 54. Pannebakker BA, Niehuis O, Hedley A, Gadau J, Shuker DM (2010) The distribution of microsatellites in the Nasonia parasitoid wasp genome. Insect Mol Biol 19: 91–98 Available: http://dx.doi.org/10.1111/j.1365-2583.2009.00915.x.
- 55. Leo NP, Hughes JM, Yang X, Poudel SKS, Brogdon WG, et al. (2005) The head and body lice of humans are genetically distinct (Insecta: Phthiraptera, Pediculidae): evidence from double infestations. Heredity (Edinb) 95: 34–40 Available: http://www.ncbi.nlm.nih.gov/pubmed/15931254.
- 56. Nair S, Williams JT, Brockman A, Paiphun L, Mayxay M, et al. (2003) A selective sweep driven by pyrimethamine treatment in southeast asian malaria parasites. Mol Biol Evol 20: 1526–1536 Available: http://www.ncbi.nlm.nih.gov/pubmed/12832643.
- 57. Roper C, Pearce R, Bredenkamp B, Gumede J, Drakeley C, et al. (2003) Antifolate antimalarial resistance in southeast Africa: a population-based analysis. Lancet 361: 1174–1181 Available: http://www.ncbi.nlm.nih.gov/pubmed/12686039.
- 58. Mita T, Tanabe K, Kita K (2009) Spread and evolution of Plasmodium falciparum drug resistance. Parasitol Int 58: 201–209 Available: http://www.ncbi.nlm.nih.gov/pubmed/19393762.
- 59. Hartley CJ, Newcomb RD, Russell RJ, Yong CG, Stevens JR, et al. (2006) Amplification of DNA from preserved specimens shows blowflies were preadapted for the rapid evolution of insecticide resistance. Proc Natl Acad Sci USA 103: 8757–8762 Available: http://www.ncbi.nlm.nih.gov/pubmed/16723400.
- 60. Bialek R, Zelck UE, Fölster-Holst R (2011) Permethrin treatment of head lice with knockdown resistance-like gene. N Engl J Med 364: 386–387 Available: http://www.ncbi.nlm.nih.gov/pubmed/21268748.
- 61. Oakeshott JG, Home I, Sutherland TD, Russell RJ (2003) The genomics of insecticide resistance. Genome Biol 4: 202 Available: http://www.ncbi.nlm.nih.gov/pubmed/12540295.
- 62. Pedra JHF, McIntyre LM, Scharf ME, Pittendrigh BR (2004) Genome-wide transcription profile of field- and laboratory-selected dichlorodiphenyltrichloroethane (DDT)-resistant Drosophila. Proc Natl Acad Sci USA 101: 7034–7039 Available: http://www.ncbi.nlm.nih.gov/pubmed/15118106.
- 63. Figueroa CC, Prunier-Leterme N, Rispe C, Sepulveda F, Fuentes-Contreras E, et al. (2007) Annotated expressed sequence tags and xenobiotic detoxification in the aphid Myzus persicae (Sulzer). Insect Sci 14: 29–45 Available: http://dx.doi.org/10.1111/j.1744-7917.2007.00123.x.
- 64. Turbeck BO, Foder B (1970) Studies on a carbonic anhydrase from the midgut epithelium of larvae of lepidoptera. Biochim Biophys Acta 212: 139–149 Available: http://www.ncbi.nlm.nih.gov/pubmed/4993464.
- 65. Corena MP, Seron TJ, Lehman HK, Ochrietor JD, Kohn A, et al. (2002) Carbonic anhydrase in the midgut of larval Aedes aegypti: cloning, localization and inhibition. J Exp Biol 205: 591–602 Available: http://www.ncbi.nlm.nih.gov/pubmed/11907049.
- 66. Linser PJ, Smith KE, Seron TJ, Neira Oviedo M (2009) Carbonic anhydrases and anion transport in mosquito midgut pH regulation. J Exp Biol 212: 1662–1671 Available: http://www.ncbi.nlm.nih.gov/pubmed/19448076.
- 67. Toups MA, Kitchen A, Light JE, Reed DL (2011) Origin of clothing lice indicates early clothing use by anatomically modern humans in Africa. Mol Biol Evol 28: 29–32.
- 68. Torroni A, Sukernik RI, Schurr TG, Starikorskaya YB, Cabell MF, et al. (1993) mtDNA variation of aboriginal Siberians reveals distinct genetic affinities with Native Americans. Am J Hum Genet 53: 591–608 Available: http://www.ncbi.nlm.nih.gov/pubmed/7688933.
- 69. Kolman CJ, Sambuughin N, Bermingham E (1996) Mitochondrial DNA analysis of Mongolian populations and implications for the origin of New World founders. Genetics 142: 1321–1334 Available: http://www.ncbi.nlm.nih.gov/pubmed/8846908.
- 70. Ramachandran S, Deshpande O, Roseman CC, Rosenberg NA, Feldman MW, et al. (2005) Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc Natl Acad Sci USA 102: 15942–15947 Available: http://www.ncbi.nlm.nih.gov/pubmed/16243969.
- 71. Abi-Rached L, Jobin MJ, Kulkarni S, McWhinnie A, Dalva K, et al. (2011) The Shaping of Modern Human Immune Systems by Multiregional Admixture with Archaic Humans. Science 334: 89–94 Available: http://www.sciencemag.org/cgi/doi/10.1126/science.1209202.
- 72. Hammer MF, Woerner AE, Mendez FL, Watkins JC, Wall JD (2011) Genetic evidence for archaic admixture in Africa. Proc Natl Acad Sci USA 108: 15123–15128 Available: http://www.pnas.org/cgi/doi/10.1073/pnas.1109300108.
- 73. Lalueza-Fox C, Gilbert MTP (2011) Paleogenomics of Archaic Hominins. Curr Biol 21: R1002–R1009 Available: http://linkinghub.elsevier.com/retrieve/pii/S096098221101270X.
- 74. Alves I, Srámková Hanulová A, Foll M, Excoffier L (2012) Genomic data reveal a complex making of humans. PLoS Genet 8: e1002837 Available: http://www.ncbi.nlm.nih.gov/pubmed/22829785.
- 75. Yang MA, Malaspinas A-S, Durand EY, Slatkin M (2012) Ancient structure in Africa unlikely to explain Neanderthal and Non-African genetic similarity. Mol Biol Evol 29(10): 2987–2995 Available: http://mbe.oxfordjournals.org/cgi/doi/10.1093/molbev/mss117.
- 76. Sankararaman S, Patterson N, Li H, Pääbo S, Reich D (2012) The date of interbreeding between Neandertals and modern humans. PLoS Genetics 8(10): e1002947 Available: http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.1002947.