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.

Stochasticity, determinism, and contingency shape genome evolution of endosymbiotic bacteria.

Evolution results from the interaction of stochastic and deterministic processes that create a web of historical contingency, shaping gene content and organismal function. To understand the scope of this interaction, we examine the relative contributions of stochasticity, determinism, and contingency in shaping gene inactivation in 34 lineages of endosymbiotic bacteria, Sodalis, found in parasitic lice, Columbicola, that are independently undergoing genome degeneration. Here we show that the process of genome degeneration in this system is largely deterministic: genes involved in amino acid biosynthesis are lost while those involved in providing B-vitamins to the host are retained. In contrast, many genes encoding redundant functions, including components of the respiratory chain and DNA repair pathways, are subject to stochastic loss, yielding historical contingencies that constrain subsequent losses. Thus, while selection results in functional convergence between symbiont lineages, stochastic mutations initiate distinct evolutionary trajectories, generating diverse gene inventories that lack the functional redundancy typically found in free-living relatives.

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 (Ne) 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 nonmodel 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, nonmodel 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.

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.

Phylogenomics reveals the origin of mammal lice out of Afrotheria.

Mammals host a wide diversity of parasites. Lice, comprising more than 5,000 species, are one group of ectoparasites whose major lineages have a somewhat patchwork distribution across the major groups of mammals. Here we explored patterns in the diversification of mammalian lice by reconstructing a higher-level phylogeny of these lice, leveraging whole genome sequence reads to assemble single-copy orthologue genes across the genome. The evolutionary tree of lice indicated that three of the major lineages of placental mammal lice had a single common ancestor. Comparisons of this parasite phylogeny with that for their mammalian hosts indicated that the common ancestor of elephants, elephant shrews and hyraxes (that is, Afrotheria) was the ancestral host of this group of lice. Other groups of placental mammals obtained their lice via host-switching out of these Afrotherian ancestors. In addition, reconstructions of the ancestral host group (bird versus mammal) for all parasitic lice supported an avian ancestral host, indicating that the ancestor of Afrotheria acquired these parasites via host-switching from an ancient avian host. These results shed new light on the long-standing question of why the major groups of parasitic lice are not uniformly distributed across mammals and reveal the origins of mammalian lice.

Independent evolution of highly variable, fragmented mitogenomes of parasitic lice.

The mitochondrial genomes (mitogenomes) of bilaterian animals are highly conserved structures that usually consist of a single circular chromosome. However, several species of parasitic lice (Insecta: Phthiraptera) possess fragmented mitogenomes, where the mitochondrial genes are present on separate, circular chromosomes. Nevertheless, the extent, causes, and consequences of this structural variation remain poorly understood. Here, we combined new and existing data to better understand the evolution of mitogenome fragmentation in major groups of parasitic lice. We found strong evidence that fragmented mitogenomes evolved many times within parasitic lice and that the level of fragmentation is highly variable, including examples of heteroplasmic arrangements. We also found a significant association between mitochondrial fragmentation and signatures of relaxed selection. Mitochondrial fragmentation was also associated with changes to a lower AT%, possibly due to differences in mutation biases. Together, our results provide a significant advance in understanding the process of mitogenome fragmentation and provide an important perspective on mitochondrial evolution in eukaryotes.

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.

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.

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.

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.

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.