- Research article
- Open Access
Evolutionary history of mammalian sucking lice (Phthiraptera: Anoplura)
© Light et al; licensee BioMed Central Ltd. 2010
Received: 17 May 2010
Accepted: 22 September 2010
Published: 22 September 2010
Sucking lice (Phthiraptera: Anoplura) are obligate, permanent ectoparasites of eutherian mammals, parasitizing members of 12 of the 29 recognized mammalian orders and approximately 20% of all mammalian species. These host specific, blood-sucking insects are morphologically adapted for life on mammals: they are wingless, dorso-ventrally flattened, possess tibio-tarsal claws for clinging to host hair, and have piercing mouthparts for feeding. Although there are more than 540 described species of Anoplura and despite the potential economical and medical implications of sucking louse infestations, this study represents the first attempt to examine higher-level anopluran relationships using molecular data. In this study, we use molecular data to reconstruct the evolutionary history of 65 sucking louse taxa with phylogenetic analyses and compare the results to findings based on morphological data. We also estimate divergence times among anopluran taxa and compare our results to host (mammal) relationships.
This study represents the first phylogenetic hypothesis of sucking louse relationships using molecular data and we find significant conflict between phylogenies constructed using molecular and morphological data. We also find that multiple families and genera of sucking lice are not monophyletic and that extensive taxonomic revision will be necessary for this group. Based on our divergence dating analyses, sucking lice diversified in the late Cretaceous, approximately 77 Ma, and soon after the Cretaceous-Paleogene boundary (ca. 65 Ma) these lice proliferated rapidly to parasitize multiple mammalian orders and families.
The diversification time of sucking lice approximately 77 Ma is in agreement with mammalian evolutionary history: all modern mammal orders are hypothesized to have diverged by 75 Ma thus providing suitable habitat for the colonization and radiation of sucking lice. Despite the concordant timing of diversification events early in the association between anoplurans and mammals, there is substantial conflict between the host and parasite phylogenies. This conflict is likely the result of a complex history of host switching and extinction events that occurred throughout the evolutionary association between sucking lice and their mammalian hosts. It is unlikely that there are any ectoparasite groups (including lice) that tracked the early and rapid radiation of eutherian mammals.
Lice (Insecta: Phthiraptera) are obligate, permanent ectoparasites of birds and mammals, entirely dependent upon their vertebrate hosts for survival. Four phthirapteran suborders are recognized: the chewing louse suborders Amblycera, Ischnocera, and Rhynchophthirina, and the sucking louse suborder Anoplura . As a group, chewing lice parasitize birds and mammals, and all have chewing mouthparts that they use to feed upon the skin (feathers, fur, and dander) and sometimes the blood of their hosts . Sucking lice, in contrast, parasitize only eutherian mammals and they are morphologically adapted for life on their mammal hosts: they are wingless, dorso-ventrally flattened, possess adaptive tibio-tarsal claws that are used to cling to host hair, and have modified piercing mouthparts for feeding. These ectoparasitic insects are one of only a handful of haematophagous arthropod groups that use their highly derived mouthparts to feed directly from host blood vessels .
The blood-feeding habits of sucking lice and the close association they have with their mammalian hosts are hypothesized to have evolved via a particular sequence of events [4–6]. Early in their evolutionary history, sucking louse ancestors had simple chewing mouthparts and were free-living associates of the nests and burrows of vertebrates. These nests and burrows served as protective habitats as well as a source of unlimited food supplies such as fungi, dung, and organic debris, specifically sloughed skin, fur, and feathers [5–7]. Over time, some of these nest associates became more directly dependent on their hosts and transitioned from opportunistic associates to obligate parasites. These parasitic species fed directly from their hosts, ingesting more nutritious and easier to digest blood (compared to organic debris such as sloughed skin and feathers) and subsequently developed specialized mouthparts modified to obtain blood meals [2, 7]. This succession of events from free-living nest associates to obligate parasites can be seen by examining the closest living relatives of the Anoplura, members of other phthirapteran suborders, and the bark and book lice (order Psocoptera). The Psocoptera are closely related to the Phthiraptera (together they form the superorder Psocodea), and these non-parasitic insects often interact with vertebrate taxa, living in the nests, burrows, or among the fur and feathers of mammals and birds and use their chewing mouthparts to feed on fungi or organic matter [7–10]. Within the Phthiraptera, phylogenetic studies have shown sucking lice to be monophyletic and derived, nested within the chewing lice and sister to the Rhynchophthirina, a small suborder of chewing lice (3 known species) parasitic on warthogs, bush pigs, and elephants [1, 11–13]. Rhynchophthirina species have modified chewing mouthparts attached to the end of a long proboscis that are used to break through the skin of their hosts allowing pools of blood to form. These chewing lice then use their mouthparts to feed on the blood collected in these pools. Thus, it is likely that sucking lice evolved from a blood-feeding Rhynchophthirina-like ancestor with the highly modified anopluran mouthparts derived from the ancestral chewing mouthparts found in all other lice [14–16].
Generalized Anoplura host associations*
(Number of Species)
Common Names of Major Host Groups
Seals, Walrus (Carnivora)
Earless Seals (Carnivora)
Earless Seals (Carnivora)
Sea Lions (Carnivora)
Giant Flying Squirrels (Rodentia)
Flying Squirrels (Rodentia)
Giant Flying Squirrels (Rodentia)
Artiodactyls, Equids (Artiodactyla, Perissodactyla)
Rodents, Pikas (Rodentia, Lagomorpha)
Asian Rock Squirrel (Rodentia)
Spiny Rats (Rodentia)
Artiodactyls (Artiodactyla), Canids (Carnivora)
Bovids, Cervids (Artiodactyla)
Camels, Lamas (Artiodactyla)
Elephant Shrews (Macroscelidea)
New World Primates (Primates)
Old World Primates (Primates)
Chinchilla Rats (Rodentia)
Spiny Rats (Rodentia)
Tree Shrews (Scandentia)
Heteromyid Rodents (Rodentia)
Rabbits and Hares (Lagomorpha)
Dwarf Lemurs (Primates)
Bush Babies (Primates)
Pygmy Dormice (Rodentia)
Squirrels, Murids (Rodentia)
Rodents, Shrews (Rodentia, Soricomorpha)
Pouched Rats (Rodentia)
Tree Shrews (Scandentia)
Cane, Dassie Rats (Rodentia)
Pygmy Dormice (Rodentia)
Old World Primates (Primates)
To date, there have been only two studies that have attempted to reconstruct phylogenetic relationships among anopluran families [21, 22]. Kim and Ludwig  studied 15 taxa based on 22 morphological characters whereas Kim  examined 47 taxa with 39 morphological characters. Both studies were based on small morphological datasets with few characters, the phylogenetic utility of which has been questioned . While there have been multiple studies examining anopluran relationships within genera and among apparently closely related genera and families [26–29], a higher-level phylogeny of sucking lice is lacking. This study is the first to use molecular data and estimates of divergence times to elucidate the evolutionary history of this unique haematophagous group in relation to their mammalian hosts.
Taxon Sampling, Data Collection, and Phylogenetic Analyses
Lice were obtained from 8 of the 15 sucking louse families (Additional File 1). Unfortunately, louse data from the remaining 7 families could not be obtained due to specimen rarity (some anopluran families have extremely narrow host ranges and some are monotypic, making collection almost impossible) or failure to amplify specimens in the laboratory. Because of PCR failure as well as availability of data from GenBank, some molecular data could not be collected. The genes 18S, EF-1α, and COI were not collected from 11, 7, and 2 specimens, respectively, and, except for one sample of Pedicinus pictus (Pedicinus pictus 2), none of the specimens analyzed were missing data from more than one molecular marker (Additional File 1).
For each gene examined, phylogenetic analyses (MP, ML, and Bayesian) yielded similar topologies, although nodal support and placement of outgroup taxa varied depending on the gene (See Additional Files 2, 3, and 4). Previous studies have noted that the third codon positions of the mitochondrial COI gene tend to be saturated and homoplasious . Saturation plots supported slight saturation of third positions in the COI data analyzed herein; however, phylogenetic analyses including and excluding third codon positions did not result in significant differences to tree topologies or branch lengths with the exception of placement of some of the outgroup taxa (data available upon request). The three genes used in this study were phylogenetically informative at different areas of the phylogeny, similar to previous findings . With 58, 126, and 220 parsimony informative sites for 18S, EF-1α, and COI, respectively, the nuclear genes (18S and EF-1α) provided slightly more resolution basally whereas the mitochondrial marker, COI, provided more resolution at terminal nodes (See Additional Files 2, 3, and 4). Although basal resolution was generally lacking and there were overall differences in resolution for each gene, topologies resulting from phylogenetic analyses of individual genes were not in strong conflict.
All phylogenetic trees reconstructed in our analyses support a monophyletic Anoplura, sister to the chewing louse suborder Rhychophthirina (Figure 2). Several anopluran families (Hoplopleuridae and Polyplacidae), genera (Hoplopleura and Pterophthirus), and species (H. ferrisi) were not monophyletic. Phylogenetic constraints forcing the families and genera to be monophyletic were significantly worse than the best tree (ML Shimodaira-Hasegawa tests and examination of Bayesian suboptimal trees; P < 0.05). The remaining louse families were all monophyletic, although it is important to note that taxon sampling for many of these groups was low (Figure 2 and Additional Files 2, 3, and 4). Notably, primate lice (excluding the polyplacid louse Lemurpediculus verruculosus) belonging to the anopluran families Pedicinidae, Pediculidae, and Pthiridae formed a highly supported monophyletic group (Figure 2).
For the most part, lice did not form monophyletic groups according to host associations. There are two clades of lice parasitizing artiodactyl mammals (Linognathus and Haematopinus, belonging to the families Linognathidae and Haematopinidae, respectively), probably resulting from two independent colonization events of the hosts (Figure 2). Lice parasitizing rodents were far from monophyletic, and instead were scattered throughout the phylogeny in five distinct clades (Figure 2). One rodent-louse clade containing the polyplacid genus Neohaematopinus also included the lice Haemodipsus and Sathrax. These two louse genera do not parasitize rodents and instead are associated with rabbits and hares, and tree shrews, respectively. Within Rodentia, lice parasitizing the families Cricetidae, Muridae, and Sciuridae also were not monophyletic (Figure 2). Sciurid lice (Hoplopleura, Linognathoides, and Neohaematopinus) are distributed among three clades and cricetid (Hoplopleura and Neohaematopinus) and murid (Polyplax and Hoplopleura) lice are each distributed across two clades. Some host lineages, however, were parasitized by monophyletic lineages of lice. These host groups include heteromyid rodents (parasitized by the louse genus Fahrenholzia), carnivores (although the sample size of lice parasitizing carnivores was extremely small), and primates (except Cheirogaleidae; Figure 2).
Estimates of Divergence Times
This study represents the first phylogenetic hypothesis of sucking louse relationships using molecular data. The genes selected for use in this study (18S, EF-1α, and COI) have proven phylogenetically informative for lice in other studies [1, 11–13, 27, 28, 32]. However, even when used in combination, these genes fail to completely resolve relationships at a higher level across Anoplura, which is not unexpected given the relatively small number of species sampled in this study. There is good support for many of the louse lineages, but it is unclear how these major lineages are related because some branches are short and lack support (Figure 2). Future studies will need to increase taxon sampling as well as include additional molecular markers to better resolve phylogenetic relationships among sucking louse species.
The mitochondrial gene COI has been a staple in louse phylogenetic work (see references above); however, rapid rates of evolution and data saturation in this marker are a concern especially for higher-level studies [13, 26]. Sequence divergence within a louse morpho-species can be rather high, upwards to 15% uncorrected p-distances [27, 32, 33], and we found similar results in this study. For example, the two Ancistroplax crocidurae specimens were 13.5% divergent, the Lemurpediculus verruculosus specimens were 14.7% divergent, and the Linognathoides marmotae specimens were 12.5% divergent (all uncorrected p-distances). While these divergences may be indicative of cryptic species, it is more likely that these high numbers for the COI gene are typical for louse lineages because sequences for the nuclear markers were identical. While COI saturation was not an obvious problem in the current study, the phylogenetic placement of Polyplax serrata appears to be the result of elevated rates of evolution in this mitochondrial marker (See Figure 2 and Additional Files 3 and 4). Future studies may encounter similar problems; therefore, inclusion of other data would be helpful when attempting to resolve higher-level relationships among lice. These additional data could be molecular or morphological (see below); however, finding appropriate molecular markers for use within Phthiraptera has often been difficult  and it is likely that novel approaches will be necessary to resolve phylogenetic relationships in this insect group.
Unfortunately, not all anopluran families could be included in this study, resulting in an incomplete picture of phylogenetic relationships in this louse lineage. To better elucidate the evolutionary history of sucking lice, it is imperative that future studies increase the diversity of taxon sampling (in addition to utilizing additional data; see above). Even with incomplete taxonomic sampling, a few aspects of anopluran systematics are apparent, specifically the lack of monophyly of the families Hoplopleuridae and Polyplacidae (Figure 2). This lack of monophyly is not surprising; Hoplopleuridae and Polyplacidae are the two largest sucking louse families, with 162 and 193 described species, respectively (Table 1). Representing such a large number of species, it is likely that these two louse families are taxonomic hodgepodges in need of additional examination and substantial revision. Revision of these two families (as well as genera within both families) will not be possible without more complete taxon sampling. At a lower taxonomic level, the results presented in the current study are in agreement with previous research that focused on smaller anopluran groups [27–29, 34].
Currently, the morphology-based classification of Anoplura detailed by Kim and Ludwig  and with modifications by Durden and Musser  is followed by most researchers studying sucking lice. The molecular phylogeny reported here (Figure 2) agrees with aspects of this morphology-based classification such as the distinct familial lineages of Anoplura associated with pinnipeds (anopluran family Echinophthiriidae), bovids (Linognathidae), bovids and suids (Haematopinidae), hominids (Pediculidae and Pthiridae) and cercopithecids (Pedicinidae). The differences between the two phylogenies (compare Figures 1 and 2) are intriguing; however, it is possible that morphological features may support some of the molecular-based relationships proposed here. For example, the molecular phylogeny places the hoplopleurid genus Pterophthirus within the genus Hoplopleura (Figure 2). Morphologically, the only difference between these two genera is the extension, to varying degrees, of the second pair of paratergal plates on the abdomen. Perhaps the varying extensions of the second pair or paratergal plates evolved more than once within Hoplopleura and, as such, it may not warrant the recognition of Pterophthirus as a distinct genus . The wide separation between hoplopleurid genera Hoplopleura (including Pterophthirus) and Ancistroplax in the molecular phylogeny (Figure 2) also has morphological ramifications. Kim and Ludwig  recognized two subfamilies within the Hoplopleuridae, the only anopluran family for which they recognized subfamilies. Members of the subfamily Hoplopleurinae (genera Hoplopleura, Pterophthirus, and Paradoxophthirus) have a large continuous sternite on abdomonal segment 2 that physically connects with the corresponding paratergal plates, situated laterally. However, in members of the subfamily Haematopinoidinae (genera Ancistroplax, Haematopinoides, and Schizophthirus), the abdominal segment two sternite is clearly divided medially resulting in two separate plates. It is feasible that this morphological difference actually defines two distinct families rather than subfamilies as supported by the large separation between the two clades in Figure 2. As such, it would be beneficial for future researchers to include other relevant hoplopleurid genera in their molecular phylogenetic reconstructions of Anoplura evolutionary history to determine if the two Hoplopleura subfamilies remain genetically distinct.
The most obvious differences between the morphological (Figure 1) and molecular (Figure 2) anopluran phylogenies involve the family Polyplacidae, which is monophyletic based on morphological data but paraphyletic based on molecular data. Interestingly, the morphological definition of the Polyplacidae is quite variable. Notwithstanding the features they share with all other anopluran families, the only morphological characters that are common to all members of the Polyplacidae, as currently recognized, are the presence of 5 antennal segments, 6 pairs of spiracles on the abdomen, small forelegs, and the absence of a notal pit on the thorax . However, none of these characters are synapomorphies for Polyplacidae. Statements reflecting the morphological variability of Polyplacidae in current descriptions include: "antennae...usually sexually dimorphic," "thorax with mesothoracic phragma usually present," "abdomen with paratergites usually highly developed...and occasionally represented by small sclerites or completely lacking," "tergal and sternal plates usually highly developed and at times reduced or lacking," "male...with variously shaped basal apodeme, parameres and pseudopenis," and "female with...spermatheca usually indistinct" . This extreme morphological variability within the Polyplacidae may actually encompass more than one family as suggested by the separate polyplacid lineages shown in Figure 2. The molecular data suggest that rigorous taxonomic reassessment of what is currently treated as Polyplacidae is warranted. In fact, all five of the separate polyplacid lineages shown in Figure 2 correspond with distinct morphological characters that could be used to define separate families and other suprageneric taxa if future taxonomic reevaluation supports such action.
Without complete taxonomic sampling at the family level, it is difficult to compare Kim's  morphological hypothesis of anopluran relationships (Figure 1) to the molecular phylogeny (Figure 2). A few additional differences, however, are apparent. For one, morphological data support a sister relationships between Pedicinidae and Pthiridae, and these two louse families are closely related to Pecaroecidae, Haematopinidae, and Hybophthiridae, all to the exclusion of Pediculidae . Molecular data, however, support monophyly of all non-polyplacid primate lice (families Pedicinidae, Pediculidae, and Pthiridae) with a relatively distant relationship to the Haematopinidae (Figure 2). Monophyly of these three primate louse families has been found in previous molecular studies [28, 29]; however, molecular data from the families Hybophthiridae, and Pecaroecidae will be needed for a more rigorous comparison to the morphological study of Kim . Interestingly, Kim  noted a relatively close relationship among the families Hoplopleuridae, Linognathidae, and Polyplacidae (Figure 1). The molecular data presented herein also support a close relationship among these three families (or at least specific clades within Hoplopleuridae and Polyplacidae; Figure 2) and it will be interesting to see if these relationships hold with additional investigations. Although a more comprehensive morphological study is currently underway (Smith and Light, unpubl. data), additional molecular data and better taxon sampling will be necessary to properly compare morphology and molecules and it is likely that these data will both agree that substantial taxonomic revision of Anoplura will be necessary.
Host Associations and the Origin of Anoplura
Using a Bayesian approach implemented in the program BEAST v1.5.3 [31, 36], we estimated a late Cretaceous diversification of sucking lice (approximately 77 Ma; Figure 3 and Additional File 5). In this analysis, we utilized the 3-gene data set and three calibration points, and we allowed substitution and clock models to be unlinked. Divergence time estimates varied, however, if clock models or clock and substitution models were linked across the three data partitions, and when calibration points were not used concurrently or used as hard bounds. When clock models or clock and substitution models were linked across the three data partitions, estimates of divergence times tended to be much more recent, with an Eocene or Paleocene origin of the Anoplura, approximately 55 Ma. However, linking substitution or clock models across partitioned data sets is the equivalent of performing analyses on concatenated data sets (i.e., not applying different models of evolution to each partition). Since Bayes Factors comparisons performed herein support partitioned analyses, we concentrate our discussion below on the estimates obtained with unlinked model parameters.
Of our three calibration points, one was located basally whereas the other two were located more terminally on the louse phylogeny (Figure 3). Previous studies have found that the use of single calibration points, especially when placed either basally or terminally, can result in erroneous estimations of divergence times [28, 37–39]. Similarly, our analyses using only the basal calibration resulted in overestimates of ages at terminal nodes and analyses using only the terminal primate-louse calibrations resulted in underestimates of the ages at basal nodes. In fact, these over- and underestimates were extreme, 63 Ma for the split between Old World and hominoid-associated lice, and 33 Ma for the age of Anoplura, respectively. Other studies have found that use of hard bounds is often ill-advised (especially when there is a lack of confidence of the exact ages of the calibration) and that fossils provide poor hard maximum bounds [38, 40, 41]. Analyses using only upper bounds also resulted in underestimates for the age of Anoplura, thus providing additional support for the simultaneous use of the three calibration points. Further analyses incorporating increased taxon sampling and additional calibrations (if available) will be necessary to test the hypotheses presented here.
Because the fossil record for lice is so poor [10, 17, 18, 42], our calibrations were based on well-documented cospeciation events between these parasites and their mammalian hosts. The primate-louse calibrations have been used in several previous studies [28, 29] and have generally proven useful to better understand louse evolutionary history. For a basal calibration representing the split between Rhynchophthirina and Anoplura, we chose the time of basal diversification in placental mammals from Bininda-Emonds et al.  because we believed it reasonable that sucking lice could not have diversified until they had appropriate hosts (i.e., placental mammals) to colonize [5, 20]. Additionally, the recent discovery of two louse fossils at 44 Ma  and 100 Ma  adds further weight to support an ancient origin of lice, as does results of a recent study examining diversification times across all suborders of lice utilizing these fossils. This study found that sucking lice diversified approximately 75 Ma [Smith VS, Ford T, Johnson KP, Johnson PCD, Yoshizawa K, Light JE: Multiple lineages of lice pass through the K-Pg boundary, Submitted]. This result is similar to our own findings, further supporting use of a basal calibration point of 94-109 Ma. It is important to note that because some of the outgroup taxa used in this study (e.g., Bovicola, Felicola, and Neotrichodectes) also parasitize eutherian mammals, the Bininda-Emonds et al.  calibration could have been placed at the root of the phylogeny. However, these chewing lice are recently derived lineages within Ischnocera and likely radiated as a result of a recent host switch to eutherian mammals .
Anoplurans are one of only a handful of haematophagous arthropod groups, and as a result of their blood-feeding habits these ectoparastic insects can have severe effects on their hosts. Heavy infestations of sucking lice can cause host anemia, weight loss, damage to hides and fur (due to itching, scratching, and stains from louse feces), and general irritability, costing some livestock industries millions of dollars each year [7, 49–51]. Sucking lice can also serve as vectors of pathogenic organisms, transmitting the causative agents for wildlife and livestock diseases such as swine pox, anaplasmosis, dermatomycosis, Lebombo virus, and seal heartworm . Some anopluran species play important roles in their associations with humans. Head, body, and pubic lice (genera Pediculus and Pthirus, respectively) parasitize millions of people each year and Pediculus species are important in human epidemiology, serving as vectors of the causative agents of epidemic diseases such as trench fever (Bartonella quintana), relapsing fever (Borrelia recurrentis), and louse-borne typhus (Rickettsia prowazekii; ). Despite the potential economical and medical implications of sucking louse infestations, this study represents the first attempt to examine higher-level anopluran relationships using molecular data. Although this study produces novel findings regarding anopluran relationships, future studies with more extensive taxon sampling will be necessary to provide a better view of the evolutionary history of sucking lice.
Our analyses indicate that sucking lice diversified in the late Cretaceous, approximately 77 Ma, and soon after the Cretaceous-Paleogene boundary approximately 65 Ma, these lice proliferated rapidly to parasitize multiple mammalian orders and families (Figure 3). These dates are in agreement with mammalian evolutionary history: all modern mammal orders are hypothesized to have diverged by 75 Ma  thus providing suitable habitat for the colonization and radiation of sucking lice. Despite the concordant timing of diversification events early in the association between anoplurans and mammals, there is substantial conflict between the host and parasite phylogenies (Figure 4). This conflict may be because the free-living psocopteran ancestors of lice colonized and began feeding on the blood of their vertebrate hosts several times as host orders were diverging before the K-Pg boundary [20, 52]. Additionally, host and parasite topological differences may be the result of a complex history of host switching and extinction events that occurred throughout the evolutionary association between sucking lice and their mammalian hosts. For example, multiple mammal lineages that radiated during the Early Tertiary (60-65 Ma) are now extinct as are their ectoparasitic lice (if present; although some of these lice may have successfully switched to an extant host). Furthermore, many of the modern mammal families did not diversify until recently, during the late Eocene to the Miocene (approximately 10-40 Ma), thus providing ample new habitats for lice to colonize.
It is unlikely that there are any ectoparasite groups (including lice) that tracked the early and rapid radiation of eutherian mammals (red and yellow boxes in Figure 4). This is in some ways unfortunate because ectoparasites could potentially be used as independent data points to infer host evolutionary history, something that could have been valuable to help elucidate remaining questions regarding the mammal phylogenetic tree [53, 54]. Given the multiple host radiations and host extinction events , a combination of historical events such as extinction, host switching, and cospeciation likely dominated the evolutionary association between ectoparasites and their mammalian hosts. The sucking louse data presented herein support this scenario, and we expect that additional data from this group as well as other parasite lineages to further elucidate the complicated evolutionary history shared between parasites and their mammalian hosts.
Taxon Sampling and Data Collection
Data from 65 anopluran taxa (including 22 representatives from GenBank) representing 8 families and 48 species were included in the molecular analysis (Additional File 1). Seven outgroup taxa, representing the chewing louse suborders Rhynchophthirina and Ischnocera, also were included in the analysis. DNA was isolated from louse specimens using the DNAeasy Tissue Kit (QIAGEN Inc., Valencia, California) using louse specific protocols [11, 55]. After DNA extraction, lice were mounted on slides and retained as vouchers. Voucher specimens are deposited in the Texas A&M University Insect Collection, voucher number 684.
Due to the large number of taxa obtained from previous studies, we focused our laboratory work on three genes that are well represented in GenBank. Portions of the nuclear 18S rRNA (18S; 460 or 508 base pairs [bp] depending on alignment methodology, see below) and elongation factor 1 alpha (EF-1α; 345 bp) genes and the mitochondrial cytochrome c oxidase subunit 1 (COI; 381 bp) gene were amplified and sequenced using primers NS1 and NS2a , EF1For3 and Cho10 , and L6625 and H7005 , respectively. Double-stranded PCR amplifications, PCR purification, and sequencing of these genes were undertaken following protocols detailed in Light and Reed . Sequences were edited using Sequencher v. 4.2.2 (Gene Codes Corporation, Ann Arbor, Michigan) and primer sequences were removed and sequences trimmed in reference to the translated protein sequence using Se-AL v2.01a11  and MacClade 4.0 . The protein coding genes EF-1α and COI were aligned by eye using Se-Al v2.0a11  and louse 18S rRNA sequences were aligned using CLUSTAL W  and MUSCLE . Phylogenetic analyses of CLUSTAL W and MUSCLE alignments of 18S yielded similar topologies and branch lengths. Results presented herein rely solely on the MUSCLE 18S alignment (results based on CLUSTAL alignments are available upon request). Longer sequences obtained from GenBank were pruned for maximum overlap with sequences generated herein. All sequences are available in GenBank (Additional File 1) and alignments are available on TreeBase (http://purl.org/phylo/treebase/phylows/study/TB2:S10679; Submission ID 10669).
Phylogenetic analyses of individual and combined genes were performed using maximum parsimony (MP), maximum likelihood (ML) and Bayesian approaches. Equally weighted MP searches were performed with 10 random addition replicates and tree bisection-reconnection branch swapping using PAUP*4.0b10 . To assess nodal support, nonparametric bootstrap analyses were performed [200 pseudoreplicates and 10 random sequence additions; ]. To generate the best ML and Bayesian trees, Modeltest  and MrModelTest  were used to examine models of nucleotide substitution (56 and 24, respectively) and to choose a best-fit model of sequence evolution . Models of evolution providing the best approximation of the data using the fewest parameters were chosen for subsequent analyses according to the Akaike Information Criterion [67, 68]. The general time reversible (GTR) model, including among-site rate variation (G) and invariable sites [I; [69, 70]] was chosen as the best model of evolution in both ModelTest and MrModelTest for the 18S rRNA, EF-1α (three codon positions combined and the third codon position), COI (three codon positions combined and the first and second codon positions), and the combined 3-gene data set. The GTR+G, SYM, and HKY+G models were chosen as the best model of evolution for the first EF-1α codon position, second EF-1α codon position, and third COI codon position, respectively. Full heuristic ML searches were conducted using the best-fit model in PAUP* 4.0b10  and GARLI , and full heuristic ML bootstrap (100 pseudoreplicates) searches were conducted using the best-fit model in GARLI .
Bayesian phylogenetic analyses were performed in MrBayes 3.12 . Model parameters were treated as unknown variables with uniform priors and were estimated as part of the analysis. Bayesian analyses were initiated from random starting trees, run for 10 million generations with 4 incrementally heated chains , and sampled at intervals of 1000 generations. Two independent Bayesian analyses were run to avoid entrapment on local optima, and log-likelihood scores were compared for convergence so that burn-in generations (the first 3000 trees) could be discarded. Tracer v1.4  was used to evaluate stability of all parameter estimates following removal of burn-in generations.
The 3-gene data set also was examined with partitioned Bayesian phylogenetic analyses. Individual genes and codon positions were defined as partitions a priori, and each partition was assigned its own substitution model according to MrModelTest (see above). Partitioning schemes included non-partitioned, partitioned by gene, and partitioned by gene and by codon position of the protein coding genes COI and EF-1α. Partitioned Bayesian analyses were performed as described above. Bayesian partitioning schemes were compared using Bayes factors , which were computed using the harmonic means of the likelihoods calculated from the sump command within MrBayes. A difference of 2ln Bayes factor > 10 was used as the minimum value to discriminate between analysis schemes [30, 75].
Bayesian Estimation of Species Trees [BEST; [76, 77]] was used to simultaneously estimate gene and species trees while allowing for independent evolutionary processes for each locus. This type of approach is useful when multiple molecular markers are being used to infer a species tree and, if there are short branch lengths in the species tree, gene trees may not match the species tree resulting in incorrect inferences of species relationships [76, 78]. Because missing data in individual gene trees can result in analytical complications , all taxa that were missing data from an entire gene region were removed prior to BEST analyses. To aid with computational time and to focus on the ingroup taxa, all outgroup taxa except for Haematomyzus elephantis were removed from the data set resulting in a total of 47 individuals analyzed. BEST analyses were initiated from random starting trees, run for 60 million generations and sampled at intervals of 1000 generations. The gene mutation prior and the prior distribution for the effective population size parameter were set at (0.5, 1.5) and 0.05, respectively. The posterior distribution of species trees (post burn-in) was summarized with a 50% majority-rule consensus tree to obtain posterior probability values for species relationships.
Alternative phylogenetic hypotheses were compared statistically using the Shimodaira-Hasegawa tests as implemented in PAUP*4.0b10 [80–82]. Additionally, suboptimal trees from the Bayesian non-partitioned and partitioned analyses were examined to assess alternative phylogenetic hypotheses. The frequency of the Markov chain Monte Carlo trees in agreement with an alternative hypothesis equals the probability of that alternative hypothesis being correct . The probability of trees agreeing with alternative subfamily hypotheses was calculated by applying constraint-based filter trees implemented in PAUP*4.0b10 .
Estimates of Divergence Times
Divergence times of the sucking louse 3-gene data set were estimated using the Bayesian approach implemented in BEAST v1.5.3 [31, 36]. BEAST uses a Bayesian relaxed molecular clock while incorporating tree uncertainty in the MCMC process to infer divergence times. Before estimating divergence times, a likelihood ratio test was performed on the louse 3-gene data set using PAUP* 4.0b10  to determine if the sequence data departed significantly from clocklike behavior. In BEAST, a Yule process speciation prior and an uncorrelated lognormal model of rate variation were implemented in each analysis . Posterior probability distributions of node ages were obtained for the 3-gene alignment with analyses performed in a concatenated and partitioned framework (all model parameters were unlinked across partitions). Best-fit models of nucleotide substitution for each data set were the same as those identified above as part of the phylogenetic analyses using MrModelTest . Two separate MCMC analyses were run for 30,000,000 generations (burnin 10%) with parameters sampled every 1000 steps. Independent runs were combined using LogCombiner v.1.5.3 . Tracer v1.5  was used to measure the effective sample size of each parameter (all resulting effective sample sizes exceeded 100) and calculate the mean and upper and lower bounds of the 95% highest posterior density interval (95% HPD) for divergence times. Tree topologies were assessed using TreeAnnotator v.1.5.3  and FigTree v.1.3.1 . Bayes factors of non-partitioned and partitioned 3-gene data sets were assessed using Tracer v1.5 .
Unfortunately, fossil calibrations for lice are lacking [10, 17, 18]. But, given their parasitic nature and the potential ancient origin of lice considering the recently described 100 Ma book louse fossil [Liposcelididae; ], it is reasonable to assume that sucking lice originated on their mammalian hosts by the end of the Cretaceous, 85-100 Ma [5, 20]. Bininda-Emonds et al.  determined that time of basal diversification in placental mammals was 94-109 Ma; therefore, a calibration with a mean of 101 (and standard deviation of 3.5) was used to represent the basal split between the Rhychophthirina and Anoplura. Additionally, some studies have documented cospeciation (i.e., roughly contemporaneous speciation events) between sucking lice and their hosts, specifically primates and primate lice . The sucking louse tree can therefore be calibrated by placing host fossil dates on the corresponding node of the louse tree. Corresponding to the split between Old World monkeys and apes , a calibration of 20-25 Ma was used to represent the split between Old World monkey lice (Pedicinus spp.) and hominoid-associated lice . Furthermore, a calibration of 5-7 Ma (corresponding to the split between chimpanzees and humans) was used to represent the split between human Pediculus lice (Pediculus humanus) and the chimpanzee louse (Pediculus schaeffi; ). These three calibration points were used in combination, as well as individually, to cross-check the other calibration points, in the divergence dating analyses. Node constraints were assigned a normal prior distribution, with the standard deviations encompassing the minimum and maximum age of each calibration. Utilizing a normal distribution allows for uncertainty in the calibration estimates , which is important for our data because calibrations were taken from the host fossil record rather than from louse fossils.
We thank A. Abramov, L. Barrow, M. Branham, S. Bush, J. Cook, J. Demboski, T. Manning, P. Myers, R. Ream, J. Skinner, N. Whiteman, and S. Zohdy for providing samples of lice to use in this study. Many of the samples used in this study were collected as part of JEL's dissertation research (supported by an NSF Dissertation Improvement Grant; DEB 0308878). M. Branham and J. Zaspel provided microscopical assistance. S. Jansa offered helpful feedback, and members of the Light and G. Voelker labs at Texas A&M University made valuable comments on earlier drafts of the manuscript. We thank M. Gitzendanner and the Florida Museum of Natural History Phyloinformatics Cluster for High Performance Computing in the Life Sciences and the University of Florida Genetics Institute's Fisher Bioinformatics Cluster for analytical support. Four anonymous reviewers offered helpful comments and suggestions that greatly improved the manuscript. This is publication number 184 of the Center for Biosystematics and Biodiversity, Texas A&M University. This work was supported by a grant from the National Science Foundation (DEB 0717165) to JEL, VSS, LAD, and DLR.
- Johnson KP, Yoshizawa K, Smith VS: Multiple origins of parasitism in lice. Proc R Soc Lond Ser B-Biol Sci. 2004, 271: 1771-1776. 10.1098/rspb.2004.2798.View ArticleGoogle Scholar
- Marshall AG: The Ecology of Ectoparasitic Insects. 1981, London: Academic PressGoogle Scholar
- Snodgrass RE: The feeding apperatus of biting and sucking insects affecting man and animals. Smithsonian Miscellaneous Collections. 1944, 104: 1-113.Google Scholar
- Hopkins GHE: The host-associations of the lice of mammals. P Zool Soc Lond. 1949, 119: 387-581. 10.1111/j.1096-3642.1949.tb00888.x.View ArticleGoogle Scholar
- Lyal CHC: Phylogeny and classification of the Psocodea, wtih particular reference to the lice (Psocodea, Phthiraptera). Syst Entomol. 1985, 10: 145-165. 10.1111/j.1365-3113.1985.tb00525.x.View ArticleGoogle Scholar
- Kim KC: Evolution and host associations of Anoplura. Coevolution of Parasitic Arthropods and Mammals. Edited by: Kim KC. 1985, New York: John Wiley and Sons, 197-231.Google Scholar
- Lehane MJ: The Biology of the Blood Sucking Insects. 2005, New York: Cambridge University Press, 2View ArticleGoogle Scholar
- Pearlman JV: Some African Psocoptera found on rats. Entomologist. 1960, 92: 246-250.Google Scholar
- Mockford EL: Some Psocoptera from plumage of birds. Proc Entomol Soc Wash. 1967, 69: 307-309.Google Scholar
- Grimaldi D, Engel MS: The Evolution of the Insects. 2005, New York: Cambridge University PressGoogle Scholar
- Cruickshank RH, Johnson KP, Smith VS, Adams RJ, Clayton DH, Page RDM: Phylogenetic analysis of partial sequences of elongation factor 1 alpha identifies major groups of lice (Insecta : Phthiraptera). Mol Phylogenet Evol. 2001, 19: 202-215. 10.1006/mpev.2001.0928.View ArticlePubMedGoogle Scholar
- Barker SC, Whiting M, Johnson KP, Murrell A: Phylogeny of the lice (Insecta, Phthiraptera) inferred from small subunit rRNA. Zool Scr. 2003, 32: 407-414. 10.1046/j.1463-6409.2003.00120.x.View ArticleGoogle Scholar
- Yoshizawa K, Johnson KP: How stable is the "polyphyly of lice" hypothesis (Insecta: Psocodea)?: A comparison of phylogenetic signal in multiple genes. Mol Phylogenet Evol. 2010, 55: 939-951. 10.1016/j.ympev.2010.02.026.View ArticlePubMedGoogle Scholar
- Cummings BF: Note on the mouthparts of a species of Polyplax (Anoplura) and on the relationship between Anoplura and Mallophaga. Ann Mag Nat Hist (8). 1915, 15: 256-259.Google Scholar
- Stojanovich CJ: The head and mouthparts of the sucking lice (Insecta: Anoplura). Microentomol. 1945, 10: 1-46.Google Scholar
- Symmons S: Comparative anatomy of the Mallophagan head. Trans Zool Soc Lond. 1952, 27: 349-436. 10.1111/j.1096-3642.1952.tb00232.x.View ArticleGoogle Scholar
- Wappler T, Smith VS, Dalgleish RC: Scratching an ancient itch: an Eocene bird louse fossil. Proc R Soc Lond Ser B-Biol Sci. 2004, 271: S255-S258. 10.1098/rsbl.2003.0158.View ArticleGoogle Scholar
- Grimaldi D, Engel MS: Fossil Liposcelididae and the lice ages (Insecta : Psocodea). Proc R Soc Lond Ser B-Biol Sci. 2006, 273: 625-633. 10.1098/rspb.2005.3337.View ArticleGoogle Scholar
- Kim KC, Ludwig HW: Parallel evolution, cladistics, and classification of parasitic Psocodea. Ann Entomol Soc Am. 1982, 75: 537-548.View ArticleGoogle Scholar
- Bininda-Emonds ORP, Cardillo M, Jones KE, MacPhee RDE, Beck RMD, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A: The delayed rise of present-day mammals. Nature. 2007, 446: 507-512. 10.1038/nature05634.View ArticlePubMedGoogle Scholar
- Kim KC, Ludwig HW: The family classification of the Anoplura. Syst Entomol. 1978, 3: 249-284. 10.1111/j.1365-3113.1978.tb00120.x.View ArticleGoogle Scholar
- Kim KC: Evolutionary parallelism in Anoplura and eutherian mammals. Biosystematics of Haematophagous Insects. Edited by: Service MW. 1988, Oxford: Claredon Press, 37: 91-114.Google Scholar
- Durden LA, Musser GG: The sucking lice (Insecta, Anoplura) of the world - a taxonomic checklist with records of mammalian hosts and geographical distributions. Bull Amer Mus Nat Hist. 1994, 218: 1-90.Google Scholar
- Durden LA, Musser GG: The mammalian hosts of the sucking lice (Anoplura) of the world: a host-parasite list. Bull Soc Vector Ecol. 1994, 19: 130-168.Google Scholar
- Kim KC: Blood-sucking lice (Anoplura) of small mammals: True parasites. Micromammals and macroparasites: from evolutionary ecology to management. Edited by: Morand S, Krasnov BR, Poulin R. 2006, Tokyo: Springer, 141-160. full_text.View ArticleGoogle Scholar
- Light JE, Hafner MS: Codivergence in heteromyid rodents (Rodentia : Heteromyidae) and their sucking lice of the genus Fahrenholzia (Phthiraptera : Anoplura). Syst Biol. 2008, 57: 449-465. 10.1080/10635150802169610.View ArticlePubMedGoogle Scholar
- Light JE, Hafner MS: Phylogenetics and host associations of Fahrenholzia sucking lice (Phthiraptera : Anoplura). Syst Entomol. 2007, 32: 359-370. 10.1111/j.1365-3113.2006.00367.x.View ArticleGoogle Scholar
- Light JE, Reed DL: Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura). Mol Phylogenet Evol. 2009, 50: 376-390. 10.1016/j.ympev.2008.10.023.View ArticlePubMedGoogle Scholar
- Reed DL, Light JE, Allen JM, Kirchman JJ: Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biol. 2007, 7 (5): 7-10.1186/1741-7007-5-7.View ArticleGoogle Scholar
- Brandley MC, Schmitz A, Reeder TW: Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scincid lizards. Syst Biol. 2005, 54: 373-390. 10.1080/10635150590946808.View ArticlePubMedGoogle Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: 699-710. 10.1371/journal.pbio.0040088.View ArticleGoogle Scholar
- Johnson KP, Reed DL, Parker SLH, Kim D, Clayton DH: Phylogenetic analysis of nuclear and mitochondrial genes supports species groups for Columbicola (Insecta: Phthiraptera). Mol Phylogenet Evol. 2007, 45: 506-518. 10.1016/j.ympev.2007.07.005.View ArticlePubMedGoogle Scholar
- Malenke JR, Johnson KP, Clayton DH: Host specialization differentiates cryptic species of feather-feeding lice. Evolution. 2009, 63: 1427-1438. 10.1111/j.1558-5646.2009.00642.x.View ArticlePubMedGoogle Scholar
- Smith VS, Light JE, Durden LA: Rodent louse diversity, phylogeny, and cospeciation in the Manu Biosphere Reserve, Peru. Biol J Linnean Soc. 2008, 95: 598-610. 10.1111/j.1095-8312.2008.01069.x.View ArticleGoogle Scholar
- Johnson PT: Sucking lice of Venezuelan rodents, with remarks on related species (Anoplura). Brigham Young Univ Sci Bull, Biol Ser. 1972, 17: 1-62.Google Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMedGoogle Scholar
- Rutschmann F, Eriksson T, Salim KA, Conti E: Assessing calibration uncertainty in molecular dating: the assignment of fossils to alternative calibration points. Syst Biol. 2007, 56: 591-608. 10.1080/10635150701491156.View ArticlePubMedGoogle Scholar
- Yang ZH, Rannala B: Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006, 23: 212-226. 10.1093/molbev/msj024.View ArticlePubMedGoogle Scholar
- Yang ZH, Yoder AD: Comparison of likelihood and Bayesian methods for estimating divergence times using multiple loci and calibration points, with application to a radiation of cute-looking mouse lemur species. Syst Biol. 2003, 52: 705-716. 10.1080/10635150390235557.View ArticlePubMedGoogle Scholar
- Rannala B, Yang ZH: Inferring speciation times under an episodic molecular clock. Syst Biol. 2007, 56: 453-466. 10.1080/10635150701420643.View ArticlePubMedGoogle Scholar
- Ho SYW, Phillips MJ: Accounting for calibration uncertainty in phylogenetic estimation of evolutionary divergence times. Syst Biol. 2009, 58: 367-380. 10.1093/sysbio/syp035.View ArticlePubMedGoogle Scholar
- Dalgleish RC, Palma RL, Price RD, Smith VS: Fossil lice (Insecta: Phthiraptera) reconsidered. Syst Entomol. 2006, 31: 648-651. 10.1111/j.1365-3113.2006.00342.x.View ArticleGoogle Scholar
- Taylor J, Purvis A: Have mammals and their chewing lice diversified in parallel?. Tangled Trees: phylogeny, cospeciation, and coevolution. Edited by: Page RDM. 2002, Chicago: University of Chicago Press, 240-261.Google Scholar
- Musser GG, Durden LA, Holden ME, Light JE: Systematic review of endemic Sulawesi squirrels (Rodentia, Sciuridae), with descriptions of new species of associated sucking lice (Insecta, Anoplura), and phylogenetic and zoogeographic assessments of sciurid lice. Bull Amer Mus Nat Hist. 2010, 2010: 1-260. 10.1206/695.1.View ArticleGoogle Scholar
- Lyal CHC: Coevolution of trichodectid lice (Insecta, Phthiraptera) and their mammalian hosts. J Nat Hist. 1987, 21: 1-28. 10.1080/00222938700770011.View ArticleGoogle Scholar
- Hafner MS, Nadler SA: Phylogenetic trees support the coevolution of parasites and their hosts. Nature. 1988, 322: 258-259. 10.1038/332258a0.View ArticleGoogle Scholar
- Klompen JSH, Keirans JE, Durden LA: Three new species of ticks (Ixodida: Argasidae: Carios) from fruit bats (Chiroptera: Pteropodidae) in the Australasian region with notes on host associations. Acarologia. 1994, 36: 25-40.Google Scholar
- Hafner JC, Light JE, Hafner DJ, Hafner MS, Reddington E, Rogers DS, Riddle BR: Basal clades and molecular systematics of heteromyid rodents. J Mamm. 2007, 88: 1129-1145. 10.1644/06-MAMM-A-413R1.1.View ArticleGoogle Scholar
- Buxton PA: The Louse. An account of the lice which infest man, their medical importance and control. 1947, London: Edward-Arnold and Co, 2Google Scholar
- Clifford CM, Yunker CE, Corwin MD: Control of louse Polyplax serrata with systemic insecticides administered in silastic rubber implants. J Econ Entomol. 1967, 60: 1210-1212.View ArticlePubMedGoogle Scholar
- Sosna CB, Medleau L: External parasites - life-cycles, transmission, and the pathogenesis of disease. Vet Med. 1992, 87: 538-Google Scholar
- Brown JW, Payne RB, Mindell DP: Nuclear DNA does not reconcile 'rocks' and 'clocks' in Neoaves: a comment on Ericson et al. Biol Lett. 2007, 3: 257-260. 10.1098/rsbl.2006.0611.PubMed CentralView ArticlePubMedGoogle Scholar
- Whiteman NK, Parker PG: Using parasites to infer host popuation history: a new rationale for parasite conservation. Anim Conserv. 2005, 8: 175-181. 10.1017/S1367943005001915.View ArticleGoogle Scholar
- Nieberding C, Morand S, Libois R, Michaux JR: A parasite reveals cryptic phylogeographic history of its host. Proc R Soc Lond Ser B-Biol Sci. 2004, 271: 2559-2568. 10.1098/rspb.2004.2930.View ArticleGoogle Scholar
- Johnson KP, Clayton DH: Coevolutionary history of ecological replicates: comparing phylogenies of wing and body lice to Columbiform hosts. Tangled Trees: phylogeny, cospeciation and coevolution. Edited by: Page RDM. 2002, Chicago: University of Chicago Press, 262-286.Google Scholar
- Danforth BN, Ji S: Elongation factor-1α occurs as two copies in bees: implications for phylogenetic analysis of EF-1α sequences in insects. Mol Biol Evol. 1998, 15: 225-235.View ArticlePubMedGoogle Scholar
- Hafner MS, Sudman PD, Villablanca FX, Spradling TA, Demastes JW, Nadler SA: Disparate rates of molecular evolution in cospeciating hosts and parasites. Science. 1994, 265: 1087-1090. 10.1126/science.8066445.View ArticlePubMedGoogle Scholar
- Rambaut A: Se-AL: Sequence Alignment Editor. [http://tree.bio.ed.ac.uk/software/seal/.]
- Maddison DR, Maddison WP: MacClade. 2003, Sunderland, MA, USA: Sinauer Associates, 4.06.OSXGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: Clustal-W-improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). 2002, Sunderland, Massachusetts: Sinauer Associates, 4Google Scholar
- Felsenstein J: Confidence-limits on phylogenies - an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Nylander JAA: MrModeltest v2. Program distributed by the author. 2004, Evolutionary Biology Centre, Uppsala University, 2Google Scholar
- Posada D: jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Rannala B: Phylogenetic methods come of age: Testing hypotheses in an evolutionary context. Science. 1997, 276: 227-232. 10.1126/science.276.5310.227.View ArticlePubMedGoogle Scholar
- Posada D, Buckley TR: Model selection and model averaging in phylogenetics: Advantages of akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol. 2004, 53: 793-808. 10.1080/10635150490522304.View ArticlePubMedGoogle Scholar
- Gu X, Fu YX, Li WH: Maximum-likelihood-estimation of the heterogeneity of substitution rate among nucleotide sites. Mol Biol Evol. 1995, 12: 546-557.PubMedGoogle Scholar
- Yang ZH: Maximum-likelihood phylogenetic estimation from DNA-sequences with variable rates over sites-approximate methods. J Mol Evol. 1994, 39: 306-314. 10.1007/BF00160154.View ArticlePubMedGoogle Scholar
- Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence data sets under the maximum likelihood criterion. Ph.D. dissertation. 2006, The University of Texas at AustinGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Rambaut A, Drummond AJ: TRACER. 2004, Oxford: University of OxfordGoogle Scholar
- Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 53: 47-67. 10.1080/10635150490264699.View ArticlePubMedGoogle Scholar
- Brown JM, Lemmon AR: The importance of data partitioning and the utility of bayes factors in Bayesian phylogenetics. Syst Biol. 2007, 56: 643-655. 10.1080/10635150701546249.View ArticlePubMedGoogle Scholar
- Edwards SV, Liu L, Pearl DK: High-resolution species trees without concatenation. Proc Natl Acad Sci USA. 2007, 104: 5936-5941. 10.1073/pnas.0607004104.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu L, Pearl DK: Species trees from gene trees: Reconstructing Bayesian posterior distributions of a species phylogeny using estimated gene tree distributions. Syst Biol. 2007, 56: 504-514. 10.1080/10635150701429982.View ArticlePubMedGoogle Scholar
- Degnan JH, Rosenberg NA: Discordance of species trees with their most likely gene trees. PLoS Genet. 2006, 2: 762-768. 10.1371/journal.pgen.0020068.View ArticleGoogle Scholar
- Leaché AD: Species tree discordance traces to phylogeographic clade boundaries in North American fence lizards (Sceloporus). Syst Biol. 2009, 58: 547-559. 10.1093/sysbio/syp057.View ArticlePubMedGoogle Scholar
- Kishino H, Hasegawa M: Evaluation of the maximum-likelihood estimate of the evolutionary tree topologies from DNA-sequence data, and the branching order in Hominoidea. J Mol Evol. 1989, 29: 170-179. 10.1007/BF02100115.View ArticlePubMedGoogle Scholar
- Goldman N, Anderson JP, Rodrigo AG: Likelihood-based tests of topologies in phylogenetics. Syst Biol. 2000, 49: 652-670. 10.1080/106351500750049752.View ArticlePubMedGoogle Scholar
- Shimodaira H: An approximately unbiased test for phylogenetic tree selection. Syst Biol. 2002, 51: 492-508. 10.1080/10635150290069913.View ArticlePubMedGoogle Scholar
- Ihlen PG, Ekman S: Outline of phylogeny and character evolution in Rhizocarpon (Rhizocarpaceae, lichenized Ascomycoa) based on nuclear ITS and mitochondrial SSU ribosomal DNA sequences. Biol J Linnean Soc. 2002, 77: 535-546. 10.1046/j.1095-8312.2002.00127.x.View ArticleGoogle Scholar
- Rambaut A: FigTree. [http://tree.bio.ed.ac.uk/software/figtree/]
- Kumar S, Filipski A, Swarna V, Walker A, Hedges SB: Placing confidence limits on the molecular age of the human-chimpanzee divergence. Proc Natl Acad Sci USA. 2005, 102: 18842-18847. 10.1073/pnas.0509585102.PubMed CentralView ArticlePubMedGoogle Scholar
- Stauffer RL, Walker A, Ryder OA, Lyons-Weiler M, Hedges SB: Human and ape molecular clocks and constraints on paleontological hypotheses. J Hered. 2001, 92: 469-474. 10.1093/jhered/92.6.469.View ArticlePubMedGoogle Scholar
- Ho SYW: Calibrating molecular estimates of substitution rates and divergence times in birds. J Avian Biol. 2007, 38: 409-414.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.