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.

Genomic Diversity in the Endosymbiotic Bacteria of Human Head Lice.

Insects have repeatedly forged symbioses with heritable microbes, gaining novel traits. For the microbe, the transition to symbioses can lead to the degeneration of the symbiont's genome through transmission bottlenecks, isolation, and the loss of DNA repair enzymes. However, some insect-microbial symbioses have persisted for millions of years, suggesting that natural selection slows genetic drift and maintains functional consistency between symbiont populations. By sampling in multiple countries, we examine genomic diversity within a symbiont species, a heritable symbiotic bacterium found only in human head lice. We find that human head louse symbionts contain genetic diversity that appears to have arisen contemporaneously with the appearance of anatomically modern humans within Africa and/or during the colonization of Eurasia by humans. We predict that the observed genetic diversity underlies functional differences in extant symbiont lineages, through the inactivation of genes involved in symbiont membrane construction. Furthermore, we find evidence of additional gene losses prior to the appearance of modern humans, also impacting the symbiont membrane. From this, we conclude that symbiont genome degeneration is proceeding, via gene inactivation and subsequent loss, in human head louse symbionts, while genomic diversity is maintained. Collectively, our results provide a look into the genomic diversity within a single symbiont species and highlight the shared evolutionary history of humans, lice, and bacteria.

The mosquito Aedes aegypti requires a gut microbiota for normal fecundity, longevity and vector competence.

Mosquitoes shift from detritus-feeding larvae to blood-feeding adults that can vector pathogens to humans and other vertebrates. The sugar and blood meals adults consume are rich in carbohydrates and protein but are deficient in other nutrients including B vitamins. Facultatively hematophagous insects like mosquitoes have been hypothesized to avoid B vitamin deficiencies by carryover of resources from the larval stage. However, prior experimental studies have also used adults with a gut microbiota that could provision B vitamins. Here, we used Aedes aegypti, which is the primary vector of dengue virus (DENV), to ask if carryover effects enable normal function in adults with no microbiota. We show that adults with no gut microbiota produce fewer eggs, live longer with lower metabolic rates, and exhibit reduced DENV vector competence but are rescued by provisioning B vitamins or recolonizing the gut with B vitamin autotrophs. We conclude carryover effects do not enable normal function.

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.

Evolutionary genomics of APSE: a tailed phage that lysogenically converts the bacterium Hamiltonella defensa into a heritable protective symbiont of aphids.

BACKGROUND: Most phages infect free-living bacteria but a few have been identified that infect heritable symbionts of insects or other eukaryotes. Heritable symbionts are usually specialized and isolated from other bacteria with little known about the origins of associated phages. Hamiltonella defensa is a heritable bacterial symbiont of aphids that is usually infected by a tailed, double-stranded DNA phage named APSE. METHODS: We conducted comparative genomic and phylogenetic studies to determine how APSE is related to other phages and prophages. RESULTS: Each APSE genome was organized into four modules and two predicted functional units. Gene content and order were near-fully conserved in modules 1 and 2, which encode predicted DNA metabolism genes, and module 4, which encodes predicted virion assembly genes. Gene content of module 3, which contains predicted toxin, holin and lysozyme genes differed among haplotypes. Comparisons to other sequenced phages suggested APSE genomes are mosaics with modules 1 and 2 sharing similarities with Bordetella-Bcep-Xylostella fastidiosa-like podoviruses, module 4 sharing similarities with P22-like podoviruses, and module 3 sharing no similarities with known phages. Comparisons to other sequenced bacterial genomes identified APSE-like elements in other heritable insect symbionts (Arsenophonus spp.) and enteric bacteria in the family Morganellaceae. CONCLUSIONS: APSEs are most closely related to phage elements in the genus Arsenophonus and other bacteria in the Morganellaceae.

Global shifts in species richness have shaped carpet shark evolution.

BACKGROUND: The evolutionary processes that shape patterns of species richness in marine ecosystems are complex and may differ between organismal groups. There has been considerable interest in understanding the evolutionary processes that led to marine species richness being concentrated in specific geographical locations. In this study we focus on the evolutionary history of a group of small-to-medium sized sharks known as carpet sharks. While a few carpet shark species are widespread, the majority of carpet shark species richness is contained within a biodiversity hotspot at the boundary of the Indian and Pacific oceans. We address the significance of this biodiversity hotspot in carpet shark evolution and speciation by leveraging a rich fossil record and molecular phylogenetics to examine the prehistoric distribution of carpet sharks. RESULTS: We find that carpet sharks species richness was greatest in shallow seas connected to the Atlantic Ocean during the Late Cretaceous, but that there was a subsequent loss of biodiversity in Atlantic waters. Fossil evidence from sites in close geographic proximity to the current center of carpet shark diversity are generally restricted to younger geologic strata. CONCLUSIONS: From this data we conclude that (1) center of carpet shark biodiversity has shifted during the last 100 million years, (2) carpet sharks have repeatedly dispersed to nascent habitat (including to their current center of diversity), and (3) the current center of carpet shark biodiversity conserves lineages that have been extirpated from this prehistoric range and is a source of new carpet shark species. Our findings provide insights into the roles of marine biodiversity hotspots for higher-tropic level predators and the methods applied here can be used for additional studies of shark evolution.

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.

The Life of an Insect Endosymbiont from the Cradle to the Grave.

Host-beneficial endosymbioses, which are formed when a microorganism takes up residence inside another cell and provides a fitness advantage to the host, have had a dramatic influence on the evolution of life. These intimate relationships have yielded the mitochondrion and the plastid (chloroplast) - the ancient organelles that in part define eukaryotic life - along with many more recent associations involving a wide variety of hosts and microbial partners. These relationships are often envisioned as stable associations that appear cooperative and persist for extremely long periods of time. But recent evidence suggests that this stable state is often born from turbulent and conflicting origins, and that the apparent stability of many beneficial endosymbiotic relationships - although certainly real in many cases - is not an inevitable outcome of these associations. Here we review how stable endosymbioses form, how they are maintained, and how they sometimes break down and are reborn. We focus on relationships formed by insects and their resident microorganisms because these symbioses have been the focus of significant empirical work over the last two decades. We review these relationships over five life stages: origin, birth, middle age, old age, and death.

Fully automated sequence alignment methods are comparable to, and much faster than, traditional methods in large data sets: an example with hepatitis B virus.

Aligning sequences for phylogenetic analysis (multiple sequence alignment; MSA) is an important, but increasingly computationally expensive step with the recent surge in DNA sequence data. Much of this sequence data is publicly available, but can be extremely fragmentary (i.e., a combination of full genomes and genomic fragments), which can compound the computational issues related to MSA. Traditionally, alignments are produced with automated algorithms and then checked and/or corrected "by eye" prior to phylogenetic inference. However, this manual curation is inefficient at the data scales required of modern phylogenetics and results in alignments that are not reproducible. Recently, methods have been developed for fully automating alignments of large data sets, but it is unclear if these methods produce alignments that result in compatible phylogenies when compared to more traditional alignment approaches that combined automated and manual methods. Here we use approximately 33,000 publicly available sequences from the hepatitis B virus (HBV), a globally distributed and rapidly evolving virus, to compare different alignment approaches. Using one data set comprised exclusively of whole genomes and a second that also included sequence fragments, we compared three MSA methods: (1) a purely automated approach using traditional software, (2) an automated approach including by eye manual editing, and (3) more recent fully automated approaches. To understand how these methods affect phylogenetic results, we compared resulting tree topologies based on these different alignment methods using multiple metrics. We further determined if the monophyly of existing HBV genotypes was supported in phylogenies estimated from each alignment type and under different statistical support thresholds. Traditional and fully automated alignments produced similar HBV phylogenies. Although there was variability between branch support thresholds, allowing lower support thresholds tended to result in more differences among trees. Therefore, differences between the trees could be best explained by phylogenetic uncertainty unrelated to the MSA method used. Nevertheless, automated alignment approaches did not require human intervention and were therefore considerably less time-intensive than traditional approaches. Because of this, we conclude that fully automated algorithms for MSA are fully compatible with older methods even in extremely difficult to align data sets. Additionally, we found that most HBV diagnostic genotypes did not correspond to evolutionarily-sound groups, regardless of alignment type and support threshold. This suggests there may be errors in genotype classification in the database or that HBV genotypes may need a revision.

Both living bacteria and eukaryotes in the mosquito gut promote growth of larvae.

We recently reported that larval stage Aedes aegypti and several other species of mosquitoes grow when living bacteria are present in the gut but do not grow when living bacteria are absent. We further reported that living bacteria induce a hypoxia signal in the gut, which activates hypoxia-induced transcription factors and other processes larvae require for growth. In this study we assessed whether other types of organisms induce mosquito larvae to grow and asked if the density of non-living microbes or diet larvae are fed obviate the requirement for living organisms prior results indicated are required for growth. Using culture conditions identical to our own prior studies, we determined that inoculation density of living Escherichia coli positively affected growth rates of Ae. aegypti larvae, whereas non-living E. coli had no effect on growth across the same range of inoculation densities. A living yeast, alga, and insect cell line induced axenic Ae. aegypti first instars to grow, and stimulated similar levels of midgut hypoxia, HIF-α stabilization, and neutral lipid accumulation in the fat body as E. coli. However, the same organisms had no effect on larval growth if heat-killed. In addition, no axenic larvae molted when fed two other diets, when fed diets supplemented with heat-killed microbes or lysed and heat-killed microbes. Experiments conducted with An. gambiae yielded similar findings. Taken together, our results indicate that organisms from different prokaryotic and eukaryotic groups induce mosquito larvae to grow, whereas no conditions were identified that stimulated larvae to grow in the absence of living organisms.

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.

Culture-Facilitated Comparative Genomics of the Facultative Symbiont Hamiltonella defensa.

Many insects host facultative, bacterial symbionts that confer conditional fitness benefits to their hosts. Hamiltonella defensa is a common facultative symbiont of aphids that provides protection against parasitoid wasps. Protection levels vary among strains of H. defensa that are also differentially infected by bacteriophages named APSEs. However, little is known about trait variation among strains because only one isolate has been fully sequenced. Generating complete genomes for facultative symbionts is hindered by relatively large genome sizes but low abundances in hosts like aphids that are very small. Here, we took advantage of methods for culturing H. defensa outside of aphids to generate complete genomes and transcriptome data for four strains of H. defensa from the pea aphid Acyrthosiphon pisum. Chosen strains also spanned the breadth of the H. defensa phylogeny and differed in strength of protection conferred against parasitoids. Results indicated that strains shared most genes with roles in nutrient acquisition, metabolism, and essential housekeeping functions. In contrast, the inventory of mobile genetic elements varied substantially, which generated strain specific differences in gene content and genome architecture. In some cases, specific traits correlated with differences in protection against parasitoids, but in others high variation between strains obscured identification of traits with likely roles in defense. Transcriptome data generated continuous distributions to genome assemblies with some genes that were highly expressed and others that were not. Single molecule real-time sequencing further identified differences in DNA methylation patterns and restriction modification systems that provide defense against phage infection.

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.

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.

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.].

Two Bacterial Genera, Sodalis and Rickettsia, Associated with the Seal Louse Proechinophthirus fluctus (Phthiraptera: Anoplura).

UNLABELLED: Roughly 10% to 15% of insect species host heritable symbiotic bacteria known as endosymbionts. The lice parasitizing mammals rely on endosymbionts to provide essential vitamins absent in their blood meals. Here, we describe two bacterial associates from a louse, Proechinophthirus fluctus, which is an obligate ectoparasite of a marine mammal. One of these is a heritable endosymbiont that is not closely related to endosymbionts of other mammalian lice. Rather, it is more closely related to endosymbionts of the genus Sodalis associated with spittlebugs and feather-chewing bird lice. Localization and vertical transmission of this endosymbiont are also more similar to those of bird lice than to those of other mammalian lice. The endosymbiont genome appears to be degrading in symbiosis; however, it is considerably larger than the genomes of other mammalian louse endosymbionts. These patterns suggest the possibility that this Sodalis endosymbiont might be recently acquired, replacing a now-extinct, ancient endosymbiont. From the same lice, we also identified an abundant bacterium belonging to the genus Rickettsia that is closely related to Rickettsia ricketsii, a human pathogen vectored by ticks. No obvious masses of the Rickettsia bacterium were observed in louse tissues, nor did we find any evidence of vertical transmission, so the nature of its association remains unclear. IMPORTANCE: Many insects are host to heritable symbiotic bacteria. These heritable bacteria have been identified from numerous species of parasitic lice. It appears that novel symbioses have formed between lice and bacteria many times, with new bacterial symbionts potentially replacing existing ones. However, little was known about the symbionts of lice parasitizing marine mammals. Here, we identified a heritable bacterial symbiont in lice parasitizing northern fur seals. This bacterial symbiont appears to have been recently acquired by the lice. The findings reported here provide insights into how new symbioses form and how this lifestyle is shaping the symbiont genome.

Genome sequence of Candidatus Riesia pediculischaeffi, endosymbiont of chimpanzee lice, and genomic comparison of recently acquired endosymbionts from human and chimpanzee lice.

The obligate-heritable endosymbionts of insects possess some of the smallest known bacterial genomes. This is likely due to loss of genomic material during symbiosis. The mode and rate of this erosion may change over evolutionary time: faster in newly formed associations and slower in long-established ones. The endosymbionts of human and anthropoid primate lice present a unique opportunity to study genome erosion in newly established (or young) symbionts. This is because we have a detailed phylogenetic history of these endosymbionts with divergence dates for closely related species. This allows for genome evolution to be studied in detail and rates of change to be estimated in a phylogenetic framework. Here, we sequenced the genome of the chimpanzee louse endosymbiont (Candidatus Riesia pediculischaeffi) and compared it with the closely related genome of the human body louse endosymbiont. From this comparison, we found evidence for recent genome erosion leading to gene loss in these endosymbionts. Although gene loss was detected, it was not significantly greater than in older endosymbionts from aphids and ants. Additionally, we searched for genes associated with B-vitamin synthesis in the two louse endosymbiont genomes because these endosymbionts are believed to synthesize essential B vitamins absent in the louse's diet. All of the expected genes were present, except those involved in thiamin synthesis. We failed to find genes encoding for proteins involved in the biosynthesis of thiamin or any complete exogenous means of salvaging thiamin, suggesting there is an undescribed mechanism for the salvage of thiamin. Finally, genes encoding for the pantothenate de novo biosynthesis pathway were located on a plasmid in both taxa along with a heat shock protein. Movement of these genes onto a plasmid may be functionally and evolutionarily significant, potentially increasing production and guarding against the deleterious effects of mutation. These data add to a growing resource of obligate endosymbiont genomes and to our understanding of the rate and mode of genome erosion in obligate animal-associated bacteria. Ultimately sequencing additional louse p-endosymbiont genomes will provide a model system for studying genome evolution in obligate host associated bacteria.

What the X has to do with it: differences in regulatory variability between the sexes in Drosophila simulans.

The mechanistic basis of regulatory variation and the prevailing evolutionary forces shaping that variation are known to differ between sexes and between chromosomes. Regulatory variation of gene expression can be due to functional changes within a gene itself (cis) or in other genes elsewhere in the genome (trans). The evolutionary properties of cis mutations are expected to differ from mutations affecting gene expression in trans. We analyze allele-specific expression across a set of X substitution lines in intact adult Drosophila simulans to evaluate whether regulatory variation differs for cis and trans, for males and females, and for X-linked and autosomal genes. Regulatory variation is common (56% of genes), and patterns of variation within D. simulans are consistent with previous observations in Drosophila that there is more cis than trans variation within species (47% vs. 25%, respectively). The relationship between sex-bias and sex-limited variation is remarkably consistent across sexes. However, there are differences between cis and trans effects: cis variants show evidence of purifying selection in the sex toward which expression is biased, while trans variants do not. For female-biased genes, the X is depleted for trans variation in a manner consistent with a female-dominated selection regime on the X. Surprisingly, there is no evidence for depletion of trans variation for male-biased genes on X. This is evidence for regulatory feminization of the X, trans-acting factors controlling male-biased genes are more likely to be found on the autosomes than those controlling female-biased genes.