Mitochondrial genome fragmentation is correlated with increased rates of molecular evolution

While mitochondrial genome content and organization is quite diverse across all Eukaryotes, most bilaterian animal mitochondrial genomes (mitogenomes) exhibit highly conserved gene content and organisation, with genes typically encoded on a single circular chromosome. However, many species of parasitic lice (Insecta: Phthiraptera) are among the notable exceptions, having mitogenomes fragmented into multiple circular chromosomes. To better understand the process of mitogenome fragmentation, we conducted a large-scale genomic study of a major group of lice, Amblycera, with extensive taxon sampling. Analyses of the evolution of mitogenome structure across a phylogenomic tree of 90 samples from 53 genera revealed evidence for multiple independent origins of mitogenome fragmentation, some inferred to have occurred less than five million years ago. We leveraged these many independent origins of fragmentation to compare the rates of DNA substitution and gene rearrangement, specifically contrasting branches with fragmented and non-fragmented mitogenomes. We found that lineages with fragmented mitochondrial genomes had significantly higher rates of mitochondrial sequence evolution. In addition, lineages with fragmented mitochondrial genomes were more likely to have mitogenome gene rearrangements than those with single-chromosome mitochondrial genomes. By combining phylogenomics and mitochondrial genomics we provide a detailed portrait of mitogenome evolution across this group of insects with a remarkably unstable mitogenome structure, identifying processes of molecular evolution that are correlated with mitogenome fragmentation.

Parasite escape mechanisms drive morphological diversification in avian lice

Organisms that have repeatedly evolved similar morphologies owing to the same selective pressures provide excellent cases in which to examine specific morphological changes and their relevance to the ecology and evolution of taxa. Hosts of permanent parasites act as an independent evolutionary experiment, as parasites on these hosts are thought to be undergoing similar selective pressures. Parasitic feather lice have repeatedly diversified into convergent ecomorphs in different microhabitats on their avian hosts. We quantified specific morphological characters to determine (i) which traits are associated with each ecomorph, (ii) the quantitative differences between these ecomorphs, and (iii) if there is evidence of displacement among co-occurring lice as might be expected under louse-louse competition on the host. We used nano-computed tomography scan data of 89 specimens, belonging to four repeatedly evolved ecomorphs, to examine their mandibular muscle volume, limb length and three-dimensional head shape data. Here, we find evidence that lice repeatedly evolve similar morphologies as a mechanism to escape host defences, but also diverge into different ecomorphs related to the way they escape these defences. Lice that co-occur with other genera on a host exhibit greater morphological divergence, indicating a potential role of competition in evolutionary divergence.

Horizontal transmission maintains host specificity and codiversification of symbionts in a brood parasitic host

In host-symbiont systems, interspecific transmissions create opportunities for host switches, potentially leading to cophylogenetic incongruence. In contrast, conspecific transmissions often result in high host specificity and congruent cophylogenies. In most bird-feather mite systems, conspecific transmission is considered dominant, while interspecific transmission is supposedly rare. However, while mites typically maintain high host specificity, incongruent cophylogenies are common. To explain this conundrum, we quantify the magnitude of conspecific vs. interspecific transmission in the brood parasitic shiny cowbird (Molothrus bonariensis). M. bonariensis lacks parental care, allowing the assessment of the role of horizontal transmission alone in maintaining host specificity. We found that despite frequent interspecific interactions via foster parental care, mite species dispersing via conspecific horizontal contacts are three times more likely to colonize M. bonariensis than mites transmitted vertically via foster parents. The results highlight the previously underappreciated rate of transmission via horizontal contacts in maintaining host specificity on a microevolutionary scale. On a macroevolutionary scale, however, host switches were estimated to have occurred as frequently as codivergences. This suggests that macroevolutionary patterns resulting from rare events cannot be easily generalized from short-term evolutionary trends.

Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict

Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.

Host body size, not host population size, predicts genome-wide effective population size of parasites

The effective population size (N_e) of an organism is expected to be generally proportional to the total number of individuals in a population. In parasites, we might expect the effective population size to be proportional to host population size and host body size, because both are expected to increase the number of parasite individuals. However, among other factors, parasite populations are sometimes so extremely subdivided that high levels of inbreeding may distort these predicted relationships. Here, we used whole-genome sequence data from dove parasites (71 feather louse species of the genus Columbicola) and phylogenetic comparative methods to study the relationship between parasite effective population size and host population size and body size. We found that parasite effective population size is largely explained by host body size but not host population size. These results suggest the potential local population size (infrapopulation or deme size) is more predictive of the long-term effective population size of parasites than is the total number of potential parasite infrapopulations (i.e., host individuals).

Draft genome assemblies of the avian louse Brueelia nebulosa and its associates using long-read sequencing from an individual specimen

Sequencing high molecular weight (HMW) DNA with long-read and linked-read technologies has promoted a major increase in more complete genome sequences for non-model organisms. Sequencing approaches that rely on HMW DNA have been limited to larger organisms or pools of multiple individuals, but recent advances have allowed for sequencing from individuals of small-bodied organisms. Here, we use HMW DNA sequencing with PacBio long-reads and TELL-Seq linked-reads to assemble and annotate the genome from a single individual feather louse (Brueelia nebulosa) from a European Starling (Sturnus vulgaris). We assembled a genome with a relatively high scaffold N50 (637 kb) and with BUSCO scores (96.1%) comparable to louse genomes assembled from pooled individuals. We annotated a number of genes (10,938) similar to the human louse (Pediculus humanus) genome. Additionally, calling phased variants revealed that the Brueelia genome is more heterozygous (∼1%) then expected for a highly obligate and dispersal-limited parasite. We also assembled and annotated the mitochondrial genome and primary endosymbiont (Sodalis) genome from the individual louse, which showed evidence for heteroplasmy in the mitogenome and a reduced genome size in the endosymbiont compared to its free-living relative. Our study is a valuable demonstration of the capability to obtain high-quality genomes from individual small, non-model organisms. Applying this approach to other organisms could greatly increase our understanding of the diversity and evolution of individual genomes.

Genomic Approaches to Uncovering the Coevolutionary History of Parasitic Lice

Simple Summary

New sequencing technologies have now made it possible to sequence entire genomes for a diversity of life on earth. Parasites comprise nearly half of all species. Lice are one important group of parasites of birds and mammals, including humans. Genome sequencing approaches have been applied to this group of parasites to uncover patterns of diversification. These patterns can be compared to the patterns of diversification in their hosts. Key findings from these studies have revealed that parasitic lice likely originated on birds and then switched to mammals multiple times. Within groups of birds and mammals, the evolutionary trees of lice match those for mammal hosts more than those for birds. Genomic approaches have also revealed that individual birds and mammals harbor distinct populations of lice. Thus, these new techniques allow for the study of patterns of diversification at a wide variety of scales.

Abstract

Next-generation sequencing technologies are revolutionizing the fields of genomics, phylogenetics, and population genetics. These new genomic approaches have been extensively applied to a major group of parasites, the lice (Insecta: Phthiraptera) of birds and mammals. Two louse genomes have been assembled and annotated to date, and these have opened up new resources for the study of louse biology. Whole genome sequencing has been used to assemble large phylogenomic datasets for lice, incorporating sequences of thousands of genes. These datasets have provided highly supported trees at all taxonomic levels, ranging from relationships among the major groups of lice to those among closely related species. Such approaches have also been applied at the population scale in lice, revealing patterns of population subdivision and inbreeding. Finally, whole genome sequence datasets can also be used for additional study beyond that of the louse nuclear genome, such as in the study of mitochondrial genome fragmentation or endosymbiont function.

High levels of inbreeding with spatial and host-associated structure in lice of an endangered freshwater seal

Host-specialist parasites of endangered large vertebrates are in many cases more endangered than their hosts. In particular, low host population densities and reduced among-host transmission rates are expected to lead to inbreeding within parasite infrapopulations living on single host individuals. Furthermore, spatial population structures of directly-transmitted parasites should be concordant with those of their hosts. Using population genomic approaches, we investigated inbreeding and population structure in a host-specialist seal louse (Echinophthirius horridus) infesting the Saimaa ringed seal (Phoca hispida saimensis), which is endemic to Lake Saimaa in Finland, and is one of the most endangered pinnipeds in the world. We conducted genome resequencing of pairs of lice collected from 18 individual Saimaa ringed seals throughout the Lake Saimaa complex. Our analyses showed high genetic similarity and inbreeding between lice inhabiting the same individual seal host, indicating low among-host transmission rates. Across the lake, genetic differentiation among individual lice was correlated with their geographic distance, and assignment analyses revealed a marked break in the genetic variation of the lice in the middle of the lake, indicating substantial population structure. These findings indicate that movements of Saimaa ringed seals across the main breeding areas of the fragmented Lake Saimaa complex may in fact be more restricted than suggested by previous population-genetic analyses of the seals themselves.

Molecular phylogenetics of the avian feather louse Philopterus-complex (Phthiraptera: Philopteridae)

The avian feather louse Philopterus-complex (Phthiraptera: Ischnocera: Philopteridae) currently contains 12 genera that have been grouped together because of shared morphological characteristics. Although previously lumped into a single genus (Philopterus), more recent morphological treatments have separated the group into several different genera. Here we evaluate the status of these genera using DNA sequence data from 118 ingroup specimens belonging to ten genera in the Philopterus-complex: Australophilopterus Mey, 2004, Cinclosomicola Mey 2004, Clayiella Eichler, 1940, Corcorides Mey, 2004, Mayriphilopterus Mey, 2004, Paraphilopterus Mey 2004, Philopteroides Mey 2004, Philopterus Nitzsch, 1818, Tyranniphilopterus Mey, 2004, and Vinceopterus Gustafsson, Lei, Chu, Zou, and Bush, 2019. Our sampling includes 97 new louse-host association records. Our analyses suggest that the genus Debeauxoecus Conci, 1941, parasitic on pittas (Aves: Pittidae), is outside of the Philopterus-complex, and that there is strong support for the monophyly of a group containing the remaining genera from the complex. Some diverse genera, such as Philopterus (sensu stricto) and Mayriphilopterus are supported as monophyletic, whereas the genera Australophilopterus, Philopteroides, and Tyranniphilopterus are not. The present study is the largest phylogenetic reconstruction of avian lice belonging to the Philopterus-complex to date and suggests that further generic revision is needed in the group to integrate molecular and morphological information.

Long-distance dispersal of pigeons and doves generated new ecological opportunities for host-switching and adaptive radiation by their parasites

Adaptive radiation is an important mechanism of organismal diversification and can be triggered by new ecological opportunities. Although poorly studied in this regard, parasites are an ideal group in which to study adaptive radiations because of their close associations with host species. Both experimental and comparative studies suggest that the ectoparasitic wing lice of pigeons and doves have adaptively radiated, leading to differences in body size and overall coloration. Here, we show that long-distance dispersal by dove hosts was central to parasite diversification because it provided new ecological opportunities for parasites to speciate after host-switching. We further show that among extant parasite lineages host-switching decreased over time, with cospeciation becoming the more dominant mode of parasite speciation. Taken together, our results suggest that host dispersal, followed by host-switching, provided novel ecological opportunities that facilitated adaptive radiation by parasites.

The interplay between host biogeography and phylogeny in structuring diversification of the feather louse genus Penenirmus

Parasite diversification is influenced by many of the same factors that affect speciation of free-living organisms, such as biogeographic barriers. However, the ecology and evolution of the host lineage also has a major impact on parasite speciation. Here we explore the interplay between biogeography and host-association on the pattern of diversification in a group of ectoparasitic lice (Insecta: Phthiraptera: Penenirmus) that feeds on the feathers of woodpeckers, barbets, and honeyguides (Piciformes) and some songbirds (Passeriformes). We use whole genome sequencing of 41 ingroup and 12 outgroup samples to develop a phylogenomic dataset of DNA sequences from a reference set of 2395 single copy ortholog genes, for a total of nearly four million aligned base positions. The phylogenetic trees resulting from both concatenated and gene-tree/species-tree coalescent analyses were nearly identical and highly supported. These trees recovered the genus Penenirmus as monophyletic and identified several major clades, which tended to be associated with one major host group. However, cophylogenetic analysis revealed that host-switching was a prominent process in the diversification of this group. This host-switching generally occurred within single major biogeographic regions. We did, however, find one case in which it appears that a rare dispersal event by a woodpecker lineage from North America to Africa allowed its associated louse to colonize a woodpecker in Africa, even though the woodpecker lineage from North America never became established there.

Genome-Resolved Metagenomic Analyses Reveal the Presence of a Putative Bacterial Endosymbiont in an Avian Nasal Mite (Rhinonyssidae; Mesostigmata)

Rhinonyssidae (Mesostigmata) is a family of nasal mites only found in birds. All species are hematophagous endoparasites, which may damage the nasal cavities of birds, and also could be potential reservoirs or vectors of other infections. However, the role of members of Rhinonyssidae as disease vectors in wild bird populations remains uninvestigated, with studies of the microbiomes of Rhinonyssidae being almost non-existent. In the nasal mite (Tinaminyssus melloi) from rock doves (Columba livia), a previous study found evidence of a highly abundant putatively endosymbiotic bacteria from Class Alphaproteobacteria. Here, we expanded the sample size of this species (two different hosts- ten nasal mites from two independent samples per host), incorporated contamination controls, and increased sequencing depth in shotgun sequencing and genome-resolved metagenomic analyses. Our goal was to increase the information regarding this mite species and its putative endosymbiont. We obtained a metagenome assembled genome (MAG) that was estimated to be 98.1% complete and containing only 0.9% possible contamination. Moreover, the MAG has characteristics typical of endosymbionts (namely, small genome size an AT bias). Overall, our results support the presence of a potential endosymbiont, which is the first described for avian nasal mites to date, and improve the overall understanding of the microbiota inhabiting these mites.

The reduced genome of a heritable symbiont from an ectoparasitic feather feeding louse

Background

Feather feeding lice are abundant and diverse ectoparasites that complete their entire life cycle on an avian host. The principal or sole source of nutrition for these lice is feathers. Feathers appear to lack four amino acids that the lice would require to complete development and reproduce. Several insect groups have acquired heritable and intracellular bacteria that can synthesize metabolites absent in an insect’s diet, allowing insects to feed exclusively on nutrient-poor resources. Multiple species of feather feeding lice have been shown to harbor heritable and intracellular bacteria. We expected that these bacteria augment the louse’s diet with amino acids and facilitated the evolution of these diverse and specialized parasites. Heritable symbionts of insects often have small genomes that contain a minimal set of genes needed to maintain essential cell functions and synthesize metabolites absent in the host insect’s diet. Therefore, we expected the genome of a bacterial endosymbiont in feather lice would be small, but encode pathways for biosynthesis of amino acids.

Results

We sequenced the genome of a bacterial symbiont from a feather feeding louse (Columbicola wolffhuegeli) that parasitizes the Pied Imperial Pigeon (Ducula bicolor) and used its genome to predict metabolism of amino acids based on the presence or absence of genes. We found that this bacterial symbiont has a small genome, similar to the genomes of heritable symbionts described in other insect groups. However, we failed to identify many of the genes that we expected would support metabolism of amino acids in the symbiont genome. We also evaluated other gene pathways and features of the highly reduced genome of this symbiotic bacterium.

Conclusions

Based on the data collected in this study, it does not appear that this bacterial symbiont can synthesize amino acids needed to complement the diet of a feather feeding louse. Our results raise additional questions about the biology of feather chewing lice and the roles of symbiotic bacteria in evolution of diverse avian parasites.

Patterns of Microbiome Variation Among Infrapopulations of Permanent Bloodsucking Parasites

While interspecific variation in microbiome composition can often be readily explained by factors such as host species identity, there is still limited knowledge of how microbiomes vary at scales lower than the species level (e.g., between individuals or populations). Here, we evaluated variation in microbiome composition of individual parasites among infrapopulations (i.e., populations of parasites of the same species living on a single host individual). To address this question, we used genome-resolved and shotgun metagenomic data of 17 infrapopulations (balanced design) of the permanent, bloodsucking seal louse Echinophthirius horridus sampled from individual Saimaa ringed seals Pusa hispida saimensis. Both genome-resolved and read-based metagenomic classification approaches consistently show that parasite infrapopulation identity is a significant factor that explains both qualitative and quantitative patterns of microbiome variation at the intraspecific level. This study contributes to the general understanding of the factors driving patterns of intraspecific variation in microbiome composition, especially of bloodsucking parasites, and has implications for understanding how well-known processes occurring at higher taxonomic levels, such as phylosymbiosis, might arise in these systems.

The assembled and annotated genome of the pigeon louse Columbicola columbae, a model ectoparasite

The pigeon louse Columbicola columbae is a longstanding and important model for studies of ectoparasitism and host-parasite coevolution. However, a deeper understanding of its evolution and capacity for rapid adaptation is limited by a lack of genomic resources. Here, we present a high-quality draft assembly of the C. columbae genome, produced using a combination of Oxford Nanopore, Illumina, and Hi-C technologies. The final assembly is 208 Mb in length, with 12 chromosome-size scaffolds representing 98.1% of the assembly. For gene model prediction, we used a novel clustering method (wavy_choose) for Oxford Nanopore RNA-seq reads to feed into the MAKER annotation pipeline. High recovery of conserved single-copy orthologs (BUSCOs) suggests that our assembly and annotation are both highly complete and highly accurate. Consistent with the results of the only other assembled louse genome, Pediculus humanus, we find that C. columbae has a relatively low density of repetitive elements, the majority of which are DNA transposons. Also similar to P. humanus, we find a reduced number of genes encoding opsins, G protein-coupled receptors, odorant receptors, insulin signaling pathway components, and detoxification proteins in the C. columbae genome, relative to other insects. We propose that such losses might characterize the genomes of obligate, permanent ectoparasites with predictable habitats, limited foraging complexity, and simple dietary regimes. The sequencing and analysis for this genome were relatively low cost, and took advantage of a new clustering technique for Oxford Nanopore RNAseq reads that will be useful to future genome projects.

Combining Nuclear and Mitochondrial Loci Provides Phylogenetic Information in the Philopterus Complex of Lice (Psocodea: Ischnocera: Philopteridae)

The Philopterus Complex includes several lineages of lice that occur on birds. The complex includes the genera Philopterus (Nitzsch, 1818; Psocodea: Philopteridae), Philopteroides (Mey, 2004; Psocodea: Philopteridae), and many other lineages that have sometimes been regarded as separate genera. Only a few studies have investigated the phylogeny of this complex, all of which are based on morphological data. Here we evaluate the utility of nuclear and mitochondrial loci for recovering the phylogeny within this group. We obtained phylogenetic trees from 39 samples of the Philopterus Complex (Psocodea: Philopteridae), using sequences of two nuclear (hyp and TMEDE6) and one mitochondrial (COI) marker. We evaluated trees derived from these genes individually as well as from concatenated sequences. All trees show 20 clearly demarcated taxa (i.e., putative species) divided into five well-supported clades. Percent sequence divergence between putative species (~5-30%) for the COI gene tended to be much higher than those for the nuclear genes (~1-15%), as expected. In cases where species are described, the lineages identified based on molecular divergence correspond to morphologically defined species. In some cases, species that are host generalists exhibit additional underlying genetic variation and such cases need to be explored by further future taxonomic revisions of the Philopterus Complex.

Structure, gene order, and nucleotide composition of mitochondrial genomes in parasitic lice from Amblycera

Parasitic lice have unique mitochondrial (mt) genomes characterized by rearranged gene orders, variable genome structures, and less AT content compared to most other insects. However, relatively little is known about the mt genomes of Amblycera, the suborder sister to all other parasitic lice. Comparing among nine different genera (including representative of all seven families), we show that Amblycera have variable and highly rearranged mt genomes. Some genera have fragmented genomes that vary considerably in length, whereas others have a single mt chromosome. Notably, these genomes are more AT-biased than most other lice. We also recover genus-level phylogenetic relationships among Amblycera that are consistent with those reported from large nuclear datasets, indicating that mt sequences are reliable for reconstructing evolutionary relationships in Amblycera. However, gene order data cannot reliably recover these same relationships. Overall, our results suggest that the mt genomes of lice, already know to be distinctive, are even more variable than previously thought.

Dendritic and parallel processing of visual threats in the retina control defensive responses

Approaching predators cast expanding shadows (i.e., looming) that elicit innate defensive responses in most animals. Where looming is first detected and how critical parameters of predatory approaches are extracted are unclear. In mice, we identify a retinal interneuron (the VG3 amacrine cell) that responds robustly to looming, but not to related forms of motion. Looming-sensitive calcium transients are restricted to a specific layer of the VG3 dendrite arbor, which provides glutamatergic input to two ganglion cells (W3 and OFFα). These projection neurons combine shared excitation with dissimilar inhibition to signal approach onset and speed, respectively. Removal of VG3 amacrine cells reduces the excitation of W3 and OFFα ganglion cells and diminishes defensive responses of mice to looming without affecting other visual behaviors. Thus, the dendrites of a retinal interneuron detect visual threats, divergent circuits downstream extract critical threat parameters, and these retinal computations initiate an innate survival behavior.

Disentangling lousy relationships: Comparative phylogenomics of two sucking louse lineages parasitizing chipmunks

The evolution of obligate parasites is often interpreted in light of their hosts' evolutionary history. An expanded approach is to examine the histories of multiple lineages of parasites that inhabit similar environments on a particular host lineage. Western North American chipmunks (genus Tamias) have a broad distribution, a history of divergence with gene flow, and host two species of sucking lice (Anoplura), Hoplopleura arboricola and Neohaematopinus pacificus. From total genomic sequencing, we obtained sequences of over 1100 loci sampled across the genomes of these lice to compare their evolutionary histories and examine the roles of host association in structuring louse relationships. Within each louse species, clades are largely associated with closely related chipmunk host species. Exceptions to this pattern appear to have a biogeographic component, but differ between the two louse species. Phylogenetic relationships among these major louse clades, in both species, are not congruent with chipmunk relationships. In the context of host associations, each louse lineage has a different evolutionary history, supporting the hypothesis that host-parasite assemblages vary both across the landscape and with the taxa under investigation. In addition, the louse Hoplopleura erratica (parasitizing the eastern Tamias striatus) is embedded within H. arboricola, rendering it paraphyletic. This phylogenetic result, together with comparable divergences within H. arboricola, indicate a need for taxonomic revision. Both host divergence and biogeographic components shape parasite diversification as demonstrated by the distinctive diversification patterns of these two independently evolving lineages that parasitize the same hosts.

Comparing rates of introgression in parasitic feather lice with differing dispersal capabilities

Organisms vary in their dispersal abilities, and these differences can have important biological consequences, such as impacting the likelihood of hybridization events. However, there is still much to learn about the factors influencing hybridization, and specifically how dispersal ability affects the opportunities for hybridization. Here, using the ecological replicate system of dove wing and body lice (Insecta: Phthiraptera), we show that species with higher dispersal abilities exhibited increased genomic signatures of introgression. Specifically, we found a higher proportion of introgressed genomic reads and more reticulated phylogenetic networks in wing lice, the louse group with higher dispersal abilities. Our results are consistent with the hypothesis that differences in dispersal ability might drive the extent of introgression through hybridization.

Jorge Doña, Andrew Sweet and Kevin Johnson find that dove lice species with higher dispersal abilities have stronger genomic signatures of introgression. By using sequence data from multiple species of both wing and body lice from the same species of hosts, the authors are able to control for nearly all factors besides dispersal ability, demonstrating the power of this study system.

Lousy grouse: Comparing evolutionary patterns in Alaska galliform lice to understand host evolution and host-parasite interactions

Understanding both sides of host-parasite relationships can provide more complete insights into host and parasite biology in natural systems. For example, phylogenetic and population genetic comparisons between a group of hosts and their closely associated parasites can reveal patterns of host dispersal, interspecies interactions, and population structure that might not be evident from host data alone. These comparisons are also useful for understanding factors that drive host-parasite coevolutionary patterns (e.g., codivergence or host switching) over different periods of time. However, few studies have compared the evolutionary histories between multiple groups of parasites from the same group of hosts at a regional geographic scale. Here, we used genomic data to compare phylogenomic and population genomic patterns of Alaska ptarmigan and grouse species (Aves: Tetraoninae) and two genera of their associated feather lice: Lagopoecus and Goniodes. We used whole-genome sequencing to obtain hundreds of genes and thousands of single-nucleotide polymorphisms (SNPs) for the lice and double-digest restriction-associated DNA sequences to obtain SNPs from Alaska populations of two species of ptarmigan. We found that both genera of lice have some codivergence with their galliform hosts, but these relationships are primarily characterized by host switching and phylogenetic incongruence. Population structure was also uncorrelated between the hosts and lice. These patterns suggest that grouse, and ptarmigan in particular, share habitats and have likely had historical and ongoing dispersal within Alaska. However, the two genera of lice also have sufficient dissimilarities in the relationships with their hosts to suggest there are other factors, such as differences in louse dispersal ability, that shape the evolutionary patterns with their hosts.

Mitochondrial genomes of Columbicola feather lice are highly fragmented, indicating repeated evolution of minicircle-type genomes in parasitic lice

Most animals have a conserved mitochondrial genome structure composed of a single chromosome. However, some organisms have their mitochondrial genes separated on several smaller circular or linear chromosomes. Highly fragmented circular chromosomes (“minicircles”) are especially prevalent in parasitic lice (Insecta: Phthiraptera), with 16 species known to have between nine and 20 mitochondrial minicircles per genome. All of these species belong to the same clade (mammalian lice), suggesting a single origin of drastic fragmentation. Nevertheless, other work indicates a lesser degree of fragmentation (2–3 chromosomes/genome) is present in some avian feather lice (Ischnocera: Philopteridae). In this study, we tested for minicircles in four species of the feather louse genus Columbicola (Philopteridae). Using whole genome shotgun sequence data, we applied three different bioinformatic approaches for assembling the Columbicola mitochondrial genome. We further confirmed these approaches by assembling the mitochondrial genome of Pediculus humanus from shotgun sequencing reads, a species known to have minicircles. Columbicola spp. genomes are highly fragmented into 15–17 minicircles between ∼1,100 and ∼3,100 bp in length, with 1–4 genes per minicircle. Subsequent annotation of the minicircles indicated that tRNA arrangements of minicircles varied substantially between species. These mitochondrial minicircles for species of Columbicola represent the first feather lice (Philopteridae) for which minicircles have been found in a full mitochondrial genome assembly. Combined with recent phylogenetic studies of parasitic lice, our results provide strong evidence that highly fragmented mitochondrial genomes, which are otherwise rare across the Tree of Life, evolved multiple times within parasitic lice.

Extensive in situ radiation of feather lice on tinamous

Tinamous host the highest generic diversity of lice of any group of birds, as well as hosting representatives of all four avian feather louse ecomorphs. Although the generic diversity of tinamou feather lice is well documented, few attempts have been made to reconstruct the phylogenetic relationships among these lice. To test whether tinamou feather lice form a monophyletic group as a whole, we used whole-genome sequencing to estimate a higher-level phylogeny of tinamou feather lice, together with a broad diversity of other avian feather louse groups. In total, we analysed sequences from over 1000 genes for 48 genera of avian lice using both concatenated and coalescent approaches to estimate the phylogeny of this diverse group of avian feather lice. Although the body louse ecomorph of tinamou feather lice formed a monophyletic group, they did not strictly form a monophyletic group together with the other three ecomorphs of tinamou feather lice. In particular, a clade comprised of several feather louse genera, mainly from South America, is nested phylogenetically within tinamou lice, which also have their main centre of diversity in South America. These results suggest in situ radiation of these parasites in South America.

Putting the genome in insect phylogenomics

Next-generation sequencing technologies provide a substantial increase in the size of molecular phylogenetic datasets that can be obtained for studies of insect systematics. Several new genome reduction approaches are leveraging these technologies to generate large phylogenomic datasets: targeted amplicon sequencing, target capture, and transcriptome sequencing. Although cost effective, these approaches provide limited data for questions outside of phylogenetics. For many groups of insects, sequencing the entire genome at modest coverage is feasible. Using these genomic reads, an automated Target Restricted Assembly Method (aTRAM) can use the results of blast searches to assemble thousands of single copy ortholog genes across a group of interest. These locally assembled genes can then be compiled into very large phylogenomic datasets. These genomic libraries have the advantage in that they also contain reads from the mitochondrial genome and symbiont genomes, as well the entire insect genome, and can be leveraged for additional studies beyond phylogenetics.

Extensive host-switching of avian feather lice following the Cretaceous-Paleogene mass extinction event

Nearly all lineages of birds host parasitic feather lice. Based on recent phylogenomic studies, the three major lineages of modern birds diverged from each other before the Cretaceous-Paleogene (K-Pg) mass extinction event. In contrast, studies of the phylogeny of feather lice on birds, indicate that these parasites diversified largely after this event. However, these studies were unable to reconstruct the ancestral avian host lineage for feather lice. Here we use genome sequences of a broad diversity of lice to reconstruct a phylogeny based on 1,075 genes. By comparing this louse evolutionary tree to the avian host tree, we show that feather lice began diversifying on the common ancestor of waterfowl and landfowl, then radiated onto other avian lineages by extensive host-switching. Dating analyses and cophylogenetic comparisons revealed that two of three lineages of birds that diverged before the K-Pg boundary acquired their feather lice after this event via host-switching.

Persistence of single species of symbionts across multiple closely-related host species

Some symbiont species are highly host-specific, inhabiting only one or a very few host species, and typically have limited dispersal abilities. When they do occur on multiple host species, populations of such symbionts are expected to become genetically structured across these different host species, and this may eventually lead to new symbiont species over evolutionary timescales. However, a low number of dispersal events of symbionts between host species across time might be enough to prevent population structure and species divergence. Overall, processes of evolutionary divergence and the species status of most putative multi-host symbiont systems are yet to be investigated. Here, we used DNA metabarcoding data of 6,023 feather mites (a total of 2,225 OTU representative sequences) from 147 infracommunities (i.e., the assemblage consisting of all mites of different species collected from the same bird host individual) to investigate patterns of population genetic structure and species status of three different putative multi-host feather mite species Proctophyllodes macedo Vitzthum, 1922, Proctophyllodes motacillae Gaud, 1953, and Trouessartia jedliczkai (Zimmerman, 1894), each of which inhabits a variable number of different closely related wagtail host species (genus Motacilla). We show that mite populations from different host species represent a single species. This pattern was found in all the mite species, suggesting that each of these species is a multi-host species in which dispersal of mites among host species prevents species divergence. Also, we found evidence of limited evolutionary divergence manifested by a low but significant level of population genetic structure among symbiont populations inhabiting different host species. Our study agrees with previous studies showing a higher than expected colonization opportunities in host-specific symbionts. Indeed, our results support that these dispersal events would allow the persistence of multi-host species even in symbionts with limited dispersal capabilities, though additional factors such as the geographical structure of some bird populations may also play a role.

Rapid experimental evolution of reproductive isolation from a single natural population

Ecological speciation occurs when local adaptation generates reproductive isolation as a by-product of natural selection. Although ecological speciation is a fundamental source of diversification, the mechanistic link between natural selection and reproductive isolation remains poorly understood, especially in natural populations. Here, we show that experimental evolution of parasite body size over 4 y (approximately 60 generations) leads to reproductive isolation in natural populations of feather lice on birds. When lice are transferred to pigeons of different sizes, they rapidly evolve differences in body size that are correlated with host size. These differences in size trigger mechanical mating isolation between lice that are locally adapted to the different sized hosts. Size differences among lice also influence the outcome of competition between males for access to females. Thus, body size directly mediates reproductive isolation through its influence on both intersexual compatibility and intrasexual competition. Our results confirm that divergent natural selection acting on a single phenotypic trait can cause reproductive isolation to emerge from a single natural population in real time.

Simultaneous radiation of bird and mammal lice following the K-Pg boundary

The diversification of parasite groups often occurs at the same time as the diversification of their hosts. However, most studies demonstrating this concordance only examine single host-parasite groups. Multiple diverse lineages of ectoparasitic lice occur across both birds and mammals. Here, we describe the evolutionary history of lice based on analyses of 1107 single-copy orthologous genes from sequenced genomes of 46 species of lice. We identify three major diverse groups of lice: one exclusively on mammals, one almost exclusively on birds and one on both birds and mammals. Each of these groups radiated just after the Cretaceous-Paleogene (K-Pg) boundary, the time of the mass extinction event of the dinosaurs and rapid diversification of most of the modern lineages of birds and mammals.

Host defense triggers rapid adaptive radiation in experimentally evolving parasites

Adaptive radiation occurs when the members of a single lineage evolve different adaptive forms in response to selection imposed by competitors or predators. Iconic examples include Darwin's finches, Caribbean anoles, and Hawaiian silverswords, all of which live on islands. Although adaptive radiation is thought to be an important generator of biodiversity, most studies concern groups that have already diversified. Here, we take the opposite approach. We experimentally triggered diversification in the descendants of a single population of host-specific parasites confined to different host "islands." We show rapid adaptive divergence of experimentally evolving feather lice in response to preening, which is a bird's main defense against ectoparasites. We demonstrate that host defense exerts strong phenotypic selection for crypsis in lice transferred to different colored rock pigeons (Columba livia). During four years of experimental evolution (∼60 generations), the lice evolved heritable differences in color. Strikingly, the observed color differences spanned the range of phenotypes found among congeneric lice adapted to other species of birds. To our knowledge, this is the first real-time demonstration that microevolution is fast enough to simulate millions of years of macroevolutionary change. Our results further indicate that host-mediated selection triggers rapid divergence in the adaptive radiation of parasites, which are among the most diverse organisms on Earth.

Feather mites play a role in cleaning host feathers: New insights from DNA metabarcoding and microscopy

Parasites and other symbionts are crucial components of ecosystems, regulating host populations and supporting food webs. However, most symbiont systems, especially those involving commensals and mutualists, are relatively poorly understood. In this study, we have investigated the nature of the symbiotic relationship between birds and their most abundant and diverse ectosymbionts: the vane-dwelling feather mites. For this purpose, we studied the diet of feather mites using two complementary methods. First, we used light microscopy to examine the gut contents of 1,300 individual feather mites representing 100 mite genera (18 families) from 190 bird species belonging to 72 families and 19 orders. Second, we used high-throughput sequencing (HTS) and DNA metabarcoding to determine gut contents from 1,833 individual mites of 18 species inhabiting 18 bird species. Results showed fungi and potentially bacteria as the main food resources for feather mites (apart from potential bird uropygial gland oil). Diatoms and plant matter appeared as rare food resources for feather mites. Importantly, we did not find any evidence of feather mites feeding upon bird resources (e.g., blood, skin) other than potentially uropygial gland oil. In addition, we found a high prevalence of both keratinophilic and pathogenic fungal taxa in the feather mite species examined. Altogether, our results shed light on the long-standing question of the nature of the relationship between birds and their vane-dwelling feather mites, supporting previous evidence for a commensalistic-mutualistic role of feather mites, which are revealed as likely fungivore-microbivore-detritivore symbionts of bird feathers.

Phylogenomics and the evolution of hemipteroid insects

Hemipteroid insects (Paraneoptera), with over 10% of all known insect diversity, are a major component of terrestrial and aquatic ecosystems. Previous phylogenetic analyses have not consistently resolved the relationships among major hemipteroid lineages. We provide maximum likelihood-based phylogenomic analyses of a taxonomically comprehensive dataset comprising sequences of 2,395 single-copy, protein-coding genes for 193 samples of hemipteroid insects and outgroups. These analyses yield a well-supported phylogeny for hemipteroid insects. Monophyly of each of the three hemipteroid orders (Psocodea, Thysanoptera, and Hemiptera) is strongly supported, as are most relationships among suborders and families. Thysanoptera (thrips) is strongly supported as sister to Hemiptera. However, as in a recent large-scale analysis sampling all insect orders, trees from our data matrices support Psocodea (bark lice and parasitic lice) as the sister group to the holometabolous insects (those with complete metamorphosis). In contrast, four-cluster likelihood mapping of these data does not support this result. A molecular dating analysis using 23 fossil calibration points suggests hemipteroid insects began diversifying before the Carboniferous, over 365 million years ago. We also explore implications for understanding the timing of diversification, the evolution of morphological traits, and the evolution of mitochondrial genome organization. These results provide a phylogenetic framework for future studies of the group.

Composition and distribution of lice (Insecta: Phthiraptera) on Colombian and Peruvian birds: New data on louse-host association in the Neotropics

The diversity of permanent ectoparasites is likely underestimated due to the difficulty of collecting samples. Lice (Insecta: Phthiraptera) are permanent ectoparasites of birds and mammals; there are approximately 5,000 species described and many more undescribed, particularly in the Neotropics. We document the louse genera collected from birds sampled in Peru (2006-2007) and Colombia (2009-2016), from 22 localities across a variety of ecosystems, ranging from lowland tropical forest and Llanos to high elevation cloud forest. We identified 35 louse genera from a total of 210 bird species belonging to 37 avian families and 13 orders. These genera belong to two suborders and three families of lice: Amblycera, families Menoponidae (present on 131 bird species) and Ricinidae (39 bird species); and Ischnocera, family Philopteridae (119 bird species). We compared our bird-louse associations with data in Price et al. (2003) and recently published Neotropical studies. The majority of bird-louse associations (51.9%) were new, with most of these coming from Passeriformes, the most diverse avian order, with the most poorly known louse fauna. Finally, we found geographical variation in louse infestation and prevalence rates. With this study, we report the first comprehensive documentation of bird-louse associations for Colombia and substantially increase the known associations documented for Peru.

Paraclaravis, a new genus for the Purple-winged and Maroon-chested Ground-doves (Aves: Columbidae)

Previous molecular phylogenetic analyses and new analyses reported here demonstrate that the genus Claravis is not monophyletic and comprises two lineages, one with the species pretiosa Ferrari-Pérez, 1886 (Blue Ground-dove), and the other with two species: geoffroyi Temminck, 1811 (Purple-winged Ground-dove) and mondetoura Bonaparte, 1856 (Maroon-chested Ground-dove). Because the generic name Claravis is typified by C. pretiosa (Ferrari-Pérez, 1886), a new genus, Paraclaravis gen. nov., is described for geoffroyi Temminck, 1811 and mondetoura Bonaparte, 1856.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Phoretic dispersal influences parasite population genetic structure

Dispersal is a fundamental component of the life history of most species. Dispersal influences fitness, population dynamics, gene flow, genetic drift and population genetic structure. Even small differences in dispersal can alter ecological interactions and trigger an evolutionary cascade. Linking such ecological processes with evolutionary patterns is difficult, but can be carried out in the proper comparative context. Here, we investigate how differences in phoretic dispersal influence the population genetic structure of two different parasites of the same host species. We focus on two species of host-specific feather lice (Phthiraptera: Ischnocera) that co-occur on feral rock pigeons (Columba livia). Although these lice are ecologically very similar, "wing lice" (Columbicola columbae) disperse phoretically by "hitchhiking" on pigeon flies (Diptera: Hippoboscidae), while "body lice" (Campanulotes compar) do not. Differences in the phoretic dispersal of these species are thought to underlie observed differences in host specificity, as well as the degree of host-parasite cospeciation. These ecological and macroevolutionary patterns suggest that body lice should exhibit more genetic differentiation than wing lice. We tested this prediction among lice on individual birds and among lice on birds from three pigeon flocks. We found higher levels of genetic differentiation in body lice compared to wing lice at two spatial scales. Our results indicate that differences in phoretic dispersal can explain microevolutionary differences in population genetic structure and are consistent with macroevolutionary differences in the degree of host-parasite cospeciation.

aTRAM 2.0: An Improved, Flexible Locus Assembler for NGS Data

Massive strides have been made in technologies for collecting genome-scale data. However, tools for efficiently and flexibly assembling raw outputs into downstream analytical workflows are still nascent. aTRAM 1.0 was designed to assemble any locus from genome sequencing data but was neither optimized for efficiency nor able to serve as a single toolkit for all assembly needs. We have completely re-implemented aTRAM and redesigned its structure for faster read retrieval while adding a number of key features to improve flexibility and functionality. The software can now (1) assemble single- or paired-end data, (2) utilize both read directions in the database, (3) use an additional de novo assembly module, and (4) leverage new built-in pipelines to automate common workflows in phylogenomics. Owing to reimplementation of databasing strategies, we demonstrate that aTRAM 2.0 is much faster across all applications compared to the previous version.

Host and parasite morphology influence congruence between host and parasite phylogenies

Comparisons of host and parasite phylogenies often show varying degrees of phylogenetic congruence. However, few studies have rigorously explored the factors driving this variation. Multiple factors such as host or parasite morphology may govern the degree of phylogenetic congruence. An ideal analysis for understanding the factors correlated with congruence would focus on a diverse host-parasite system for increased variation and statistical power. In this study, we focused on the Brueelia-complex, a diverse and widespread group of feather lice that primarily parasitise songbirds. We generated a molecular phylogeny of the lice and compared this tree with a phylogeny of their avian hosts. We also tested for the contribution of each host-parasite association to the overall congruence. The two trees overall were significantly congruent, but the contribution of individual associations to this congruence varied. To understand this variation, we developed a novel approach to test whether host, parasite or biogeographic factors were statistically associated with patterns of congruence. Both host plumage dimorphism and parasite ecomorphology were associated with patterns of congruence, whereas host body size, other plumage traits and biogeography were not. Our results lay the framework for future studies to further elucidate how these factors influence the process of host-parasite coevolution.

Brochosomins and other novel proteins from brochosomes of leafhoppers (Insecta, Hemiptera, Cicadellidae)

Brochosomes (BS) are secretory granules resembling buckyballs, produced intracellularly in specialized glandular segments of the Malpighian tubules and forming superhydrophobic coatings on the integuments of leafhoppers (Hemiptera, Cicadellidae). Their composition is poorly known. Using a combination of SDS-PAGE, LC-MS/MS, next-generation sequencing (RNAseq) and bioinformatics we demonstrate that the major structural component of BS of the leafhopper Graphocephala fennahi Young is a novel family of 21-40-kDa secretory proteins, referred to herein as brochosomins (BSM), apparently cross-linked by disulfide bonds. At least 28 paralogous BSM were identified in a transcriptome assembly of this species, most of which were detected in BS. Multiple additional BS-associated proteins (BSAP), possibly loosely attached to the outer and inner surfaces of BS, were also identified; some of these were glycine-, tyrosine- and proline-rich. BSM and BSAP together accounted for half of the 100 most expressed transcripts in the Malpighian tubules of G. fennahi. Except for several minor BSAP possibly related to cyclases, BSM and BSAP had no homologs among known proteins, thus representing taxonomically restricted gene families (orphans). Searching in 50 whole-body transcriptome assemblies of Hemiptera found homologs of BSM in representatives of all five families of the superfamily Membracoidea (Cicadellidae, Myerslopiidae, Aetalionidae, Membracidae, and Melizoderidae), but not in other lineages. Among the identified proteins only BSM were shared in common between all 17 surveyed leafhoppers known to produce BS. Combined CHN elemental and aminoacid analyses estimated the total protein content of BS from the integument of G. fennahi to be 60-70%.

Phylogenomics from Whole Genome Sequences Using aTRAM

Novel sequencing technologies are rapidly expanding the size of data sets that can be applied to phylogenetic studies. Currently the most commonly used phylogenomic approaches involve some form of genome reduction. While these approaches make assembling phylogenomic data sets more economical for organisms with large genomes, they reduce the genomic coverage and thereby the long-term utility of the data. Currently, for organisms with moderate to small genomes ($<$1000 Mbp) it is feasible to sequence the entire genome at modest coverage ($10-30\times$). Computational challenges for handling these large data sets can be alleviated by assembling targeted reads, rather than assembling the entire genome, to produce a phylogenomic data matrix. Here we demonstrate the use of automated Target Restricted Assembly Method (aTRAM) to assemble 1107 single-copy ortholog genes from whole genome sequencing of sucking lice (Anoplura) and out-groups. We developed a pipeline to extract exon sequences from the aTRAM assemblies by annotating them with respect to the original target protein. We aligned these protein sequences with the inferred amino acids and then performed phylogenetic analyses on both the concatenated matrix of genes and on each gene separately in a coalescent analysis. Finally, we tested the limits of successful assembly in aTRAM by assembling 100 genes from close- to distantly related taxa at high to low levels of coverage.Both the concatenated analysis and the coalescent-based analysis produced the same tree topology, which was consistent with previously published results and resolved weakly supported nodes. These results demonstrate that this approach is successful at developing phylogenomic data sets from raw genome sequencing reads. Further, we found that with coverages above $5-10\times$, aTRAM was successful at assembling 80-90% of the contigs for both close and distantly related taxa. As sequencing costs continue to decline, we expect full genome sequencing will become more feasible for a wider array of organisms, and aTRAM will enable mining of these genomic data sets for an extensive variety of applications, including phylogenomics. [aTRAM; gene assembly; genome sequencing; phylogenomics.].

Integrating phylogenomic and population genomic patterns in avian lice provides a more complete picture of parasite evolution

Parasite diversity accounts for most of the biodiversity on earth, and is shaped by many processes (e.g., cospeciation, host switching). To identify the effects of the processes that shape parasite diversity, it is ideal to incorporate both deep (phylogenetic) and shallow (population) perspectives. To this end, we developed a novel workflow to obtain phylogenetic and population genetic data from whole genome sequences of body lice parasitizing New World ground-doves. Phylogenies from these data showed consistent, highly resolved species-level relationships for the lice. By comparing the louse and ground-dove phylogenies, we found that over long-term evolutionary scales their phylogenies were largely congruent. Many louse lineages (both species and populations) also demonstrated high host-specificity, suggesting ground-dove divergence is a primary driver of their parasites' diversity. However, the few louse taxa that are generalists are structured according to biogeography at the population level. This suggests dispersal among sympatric hosts has some effect on body louse diversity, but over deeper time scales the parasites eventually sort according to host species. Overall, our results demonstrate that multiple factors explain the patterns of diversity in this group of parasites, and that the effects of these factors can vary over different evolutionary scales. The integrative approach we employed was crucial for uncovering these patterns, and should be broadly applicable to other studies.

Development and Characterization of 12 Novel Polymorphic Microsatellite Loci for the Mammal Chewing Louse Geomydoecus aurei (Insecta: Phthiraptera) and a Comparison of Next-Generation Sequencing Approaches for Use in Parasitology

Next-generation sequencing methodologies open the door for evolutionary studies of wildlife parasites. We used 2 next-generation sequencing approaches to discover microsatellite loci in the pocket gopher chewing louse Geomydoecus aurei for use in population genetic studies. In one approach, we sequenced a library enriched for microsatellite loci; in the other approach, we mined microsatellites from genomic sequences. Following microsatellite discovery, promising loci were tested for amplification and polymorphism in 390 louse individuals from 13 pocket gopher hosts. In total, 12 loci were selected for analysis (6 from each methodology), none of which exhibited evidence of null alleles or heterozygote deficiencies. These 12 loci showed adequate genetic diversity for population-level analyses, with 3-9 alleles per locus with an average H_E per locus ranging from 0.32 to 0.70. Analysis of Molecular Variance (AMOVA) indicated that genetic variation among infrapopulations accounts for a low, but significant, percentage of the overall genetic variation, and individual louse infrapopulations showed F_ST values that were significantly different from zero in the majority of pairwise infrapopulation comparisons, despite all 13 infrapopulations being taken from the same locality. Therefore, these 12 polymorphic markers will be useful at the infrapopulation and population levels for future studies involving G. aurei. This study shows that next-generation sequencing methodologies can successfully be used to efficiently obtain data for a variety of evolutionary questions.

Phylogenomics using Target-Restricted Assembly Resolves Intrageneric Relationships of Parasitic Lice (Phthiraptera: Columbicola)

Parasitic "wing lice" (Phthiraptera: Columbicola) and their dove and pigeon hosts are a well-recognized model system for coevolutionary studies at the intersection of micro- and macroevolution. Selection on lice in microevolutionary time occurs as pigeons and doves defend themselves against lice by preening. In turn, behavioral and morphological adaptations of the lice improve their ability to evade host defense. Over macroevolutionary time wing lice tend to cospeciate with their hosts; yet, some species of Columbicola have switched to new host species. Understanding the ecological and evolutionary factors that influence coadaptation and codiversification in this system will substantially improve our understanding of coevolution in general. However, further work is hampered by the lack of a robust phylogenetic framework for Columbicola spp. and their hosts. Previous attempts to resolve the phylogeny of Columbicola based on sequences from a few genes provided limited support. Here, we apply a new approach, target restricted assembly, to assemble 977 orthologous gene sequences from whole-genome sequence data generated from very small, ethanol-preserved specimens, representing up to 61 species of wing lice. Both concatenation and coalescent methods were used to estimate the species tree. These two approaches yielded consistent and well-supported trees with 90% of all relationships receiving 100% support, which is a substantial improvement over previous studies. We used this new phylogeny to show that biogeographic ranges are generally conserved within clades of Columbicola wing lice. Limited inconsistencies are probably attributable to intercontinental dispersal of hosts, and host switching by some of the lice. [aTRAM; coalescent; coevolution; concatenation; species tree.].

Primates, Lice and Bacteria: Speciation and Genome Evolution in the Symbionts of Hominid Lice

Insects with restricted diets rely on symbiotic bacteria to provide essential metabolites missing in their diet. The blood-sucking lice are obligate, host-specific parasites of mammals and are themselves host to symbiotic bacteria. In human lice, these bacterial symbionts supply the lice with B-vitamins. Here, we sequenced the genomes of symbiotic and heritable bacterial of human, chimpanzee, gorilla, and monkey lice and used phylogenomics to investigate their evolutionary relationships. We find that these symbionts have a phylogenetic history reflecting the louse phylogeny, a finding contrary to previous reports of symbiont replacement. Examination of the highly reduced symbiont genomes (0.53–0.57 Mb) reveals much of the genomes are dedicated to vitamin synthesis. This is unchanged in the smallest symbiont genome and one that appears to have been reorganized. Specifically, symbionts from human lice, chimpanzee lice, and gorilla lice carry a small plasmid that encodes synthesis of vitamin B5, a vitamin critical to the bacteria-louse symbiosis. This plasmid is absent in an old world monkey louse symbiont, where this pathway is on its primary chromosome. This suggests the unique genomic configuration brought about by the plasmid is not essential for symbiosis, but once obtained, it has persisted for up to 25 My. We also find evidence that human, chimpanzee, and gorilla louse endosymbionts have lost a pathway for synthesis of vitamin B1, whereas the monkey louse symbiont has retained this pathway. It is unclear whether these changes are adaptive, but they may point to evolutionary responses of louse symbionts to shifts in primate biology.

Host sympatry and body size influence parasite straggling rate in a highly connected multihost, multiparasite system

Parasite lineages commonly diverge when host lineages diverge. However, when large clades of hosts and parasites are analyzed, some cases suggest host switching as another major diversification mechanism. The first step in host switching is the appearance of a parasite on an atypical host, or “straggling.” We analyze the conditions associated with straggling events. We use five species of colonially nesting seabirds from the Galapagos Archipelago and two genera of highly specific ectoparasitic lice to examine host switching. We use both genetic and morphological identification of lice, together with measurements of spatial distribution of hosts in mixed breeding colonies, to test: (1) effects of local host community composition on straggling parasite identity; (2) effects of relative host density within a mixed colony on straggling frequency and parasite species identity; and (3) how straggling rates are influenced by the specifics of louse attachment. Finally, we determine whether there is evidence of breeding in cases where straggling adult lice were found, which may indicate a shift from straggling to the initial stages of host switching. We analyzed more than 5,000 parasite individuals and found that only ~1% of lice could be considered stragglers, with ~5% of 436 host individuals having straggling parasites. We found that the presence of the typical host and recipient host in the same locality influenced straggling. Additionally, parasites most likely to be found on alternate hosts are those that are smaller than the typical parasite of that host, implying that the ability of lice to attach to the host might limit host switching. Given that lice generally follow Harrison's rule, with larger parasites on larger hosts, parasites infecting the larger host species are less likely to successfully colonize smaller host species. Moreover, our study supports the general perception that successful colonization of a novel host is extremely rare, as we found only one nymph of a straggling species, which may indicate successful reproduction.