Endoparasitoid lifestyle promotes endogenization and domestication of dsDNA viruses.

The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.

From Wolbachia genomics to phenotype: molecular models of cytoplasmic incompatibility must account for the multiplicity of compatibility types.

Wolbachia endosymbionts commonly induce cytoplasmic incompatibility, making infected males' sperm lethal to the embryos unless these are rescued by the same bacterium, inherited from their mother. Causal genes were recently identified but two families of mechanistic models are still opposed. In the toxin-antidote model, interaction between the toxin and the antidote is required for rescuing the embryos. In host modification models, a host factor is misregulated in sperm and rescue occurs through compensation or withdrawal of this modification. While these models have been thoroughly discussed, the multiplicity of compatibility types, that is, the existence of many mutually incompatible strains, as seen in Culex mosquitoes, has not received sufficient attention. To explain such a fact, host modification models must posit that the same embryonic defects can be induced and rescued through a large variety of host targets. Conversely, the toxin-antidote model simply accommodates this pattern in a lock-key fashion, through variations in the toxin-antidote interaction sites.

Natural Selection beyond Life? A Workshop Report.

Natural selection is commonly seen not just as an explanation for adaptive evolution, but as the inevitable consequence of "heritable variation in fitness among individuals". Although it remains embedded in biological concepts, such a formalisation makes it tempting to explore whether this precondition may be met not only in life as we know it, but also in other physical systems. This would imply that these systems are subject to natural selection and may perhaps be investigated in a biological framework, where properties are typically examined in light of their putative functions. Here we relate the major questions that were debated during a three-day workshop devoted to discussing whether natural selection may take place in non-living physical systems. We start this report with a brief overview of research fields dealing with "life-like" or "proto-biotic" systems, where mimicking evolution by natural selection in test tubes stands as a major objective. We contend the challenge may be as much conceptual as technical. Taking the problem from a physical angle, we then discuss the framework of dissipative structures. Although life is viewed in this context as a particular case within a larger ensemble of physical phenomena, this approach does not provide general principles from which natural selection can be derived. Turning back to evolutionary biology, we ask to what extent the most general formulations of the necessary conditions or signatures of natural selection may be applicable beyond biology. In our view, such a cross-disciplinary jump is impeded by reliance on individuality as a central yet implicit and loosely defined concept. Overall, these discussions thus lead us to conjecture that understanding, in physico-chemical terms, how individuality emerges and how it can be recognised, will be essential in the search for instances of evolution by natural selection outside of living systems.

Wolbachia host shifts: routes, mechanisms, constraints and evolutionary consequences.

Wolbachia is one of the most abundant endosymbionts on earth, with a wide distribution especially in arthropods. Effective maternal transmission and the induction of various phenotypes in their hosts are two key features of this bacterium. Here, we review our current understanding of another central aspect of Wolbachia's success: their ability to switch from one host species to another. We build on the proposal that Wolbachia host shifts occur in four main steps: (i) physical transfer to a new species; (ii) proliferation within that host; (iii) successful maternal transmission; and (iv) spread within the host species. Host shift can fail at each of these steps, and the likelihood of ultimate success is influenced by many factors. Some stem from traits of Wolbachia (different strains have different abilities for host switching), others on host features such as genetic resemblance (e.g. host shifting is likely to be easier between closely related species), ecological connections (the donor and recipient host need to interact), or the resident microbiota. Host shifts have enabled Wolbachia to reach its enormous current incidence and global distribution among arthropods in an epidemiological process shaped by loss and acquisition events across host species. The ability of Wolbachia to transfer between species also forms the basis of ongoing endeavours to control pests and disease vectors, following artificial introduction into uninfected hosts such as mosquitoes. Throughout, we emphasise the many knowledge gaps in our understanding of Wolbachia host shifts, and question the effectiveness of current methodology to detect these events. We conclude by discussing an apparent paradox: how can Wolbachia maintain its ability to undergo host shifts given that its biology seems dominated by vertical transmission?

How consistent is RAD-seq divergence with DNA-barcode based clustering in insects?

Promoted by the barcoding approach, mitochondrial DNA is more than ever used as a molecular marker to identify species boundaries. Yet, it has been repeatedly argued that it may be poorly suited for this purpose, especially in insects where mitochondria are often associated with invasive intracellular bacteria that may promote their introgression. Here, we inform this debate by assessing how divergent nuclear genomes can be when mitochondrial barcodes indicate very high proximity. To this end, we obtained RAD-seq data from 92 barcode-based species-like units (operational taxonomic units [OTUs]) spanning four insect orders. In 100% of the cases, the observed median nuclear divergence was lower than 2%, a value that was recently estimated as one below which nuclear gene flow is not uncommon. These results suggest that although mitochondria may occasionally leak between species, this process is rare enough in insects to make DNA barcoding a reliable tool for clustering specimens into species-like units.

(In)exhaustible Suppliers for Evolution? Epistatic Selection Tunes the Adaptive Potential of Nongenetic Inheritance.

Nongenetic inheritance media-from methyl-accepting cytosines to culture-tend to mutate more frequently than DNA sequences. Whether this makes them inexhaustible suppliers for adaptive evolution will depend on the effect of nongenetic mutations (hereafter, epimutations) on fitness-related traits. Here we investigate how these effects might themselves evolve, specifically whether natural selection may set boundaries to the adaptive potential of nongenetic inheritance media because of their higher mutability. In our model, the genetic and epigenetic contributions to a nonneutral phenotype are controlled by an epistatic modifier locus, which evolves under the combined effects of drift and selection. We show that a pure genetic control evolves when the environment is stable-provided that the population is large-such that the phenotype becomes robust to frequent epimutations. When the environment fluctuates, however, selection on the modifier locus also fluctuates and can overall produce a large nongenetic contribution to the phenotype, especially when the epimutation rate matches the rate of environmental variation. We further show that selection on the modifier locus is generally insensitive to recombination, meaning it is mostly direct, that is, not relying on subsequent effects in future generations. These results suggest that unstable inheritance media might significantly contribute to fitness variation of traits subject to highly variable selective pressures but little to traits responding to scarcely variable aspects of the environment. More generally, our study demonstrates that the rate of mutation and the adaptive potential of any inheritance media should not be seen as independent properties.

Categorization of species as native or nonnative using DNA sequence signatures without a complete reference library.

New genetic diagnostic approaches have greatly aided efforts to document global biodiversity and improve biosecurity. This is especially true for organismal groups in which species diversity has been underestimated historically due to difficulties associated with sampling, the lack of clear morphological characteristics, and/or limited availability of taxonomic expertise. Among these methods, DNA sequence barcoding (also known as "DNA barcoding") and by extension, meta-barcoding for biological communities, has emerged as one of the most frequently utilized methods for DNA-based species identifications. Unfortunately, the use of DNA barcoding is limited by the availability of complete reference libraries (i.e., a collection of DNA sequences from morphologically identified species), and by the fact that the vast majority of species do not have sequences present in reference databases. Such conditions are critical especially in tropical locations that are simultaneously biodiversity rich and suffer from a lack of exploration and DNA characterization by trained taxonomic specialists. To facilitate efforts to document biodiversity in regions lacking complete reference libraries, we developed a novel statistical approach that categorizes unidentified species as being either likely native or likely nonnative based solely on measures of nucleotide diversity. We demonstrate the utility of this approach by categorizing a large sample of specimens of terrestrial insects and spiders (collected as part of the Moorea BioCode project) using a generalized linear mixed model (GLMM). Using a training data set of known endemic (n = 45) and known introduced species (n = 102), we then estimated the likely native/nonnative status for 4,663 specimens representing an estimated 1,288 species (412 identified species), including both those specimens that were either unidentified or whose endemic/introduced status was uncertain. Using this approach, we were able to increase the number of categorized specimens by a factor of 4.4 (from 794 to 3,497), and the number of categorized species by a factor of 4.8 from (147 to 707) at a rate much greater than chance (77.6% accuracy). The study identifies phylogenetic signatures of both native and nonnative species and suggests several practical applications for this approach including monitoring biodiversity and facilitating biosecurity.

Global survey of mobile DNA horizontal transfer in arthropods reveals Lepidoptera as a prime hotspot.

More than any other genome components, Transposable Elements (TEs) have the capacity to move across species barriers through Horizontal Transfer (HT), with substantial evolutionary consequences. Previous large-scale surveys, based on full-genomes comparisons, have revealed the transposition mode as an important predictor of HT rates variation across TE superfamilies. However, host biology could represent another major explanatory factor, one that needs to be investigated through extensive taxonomic sampling. Here we test this hypothesis using a field collection of 460 arthropod species from Tahiti and surrounding islands. Through targeted massive parallel sequencing, we uncover patterns of HT in three widely-distributed TE superfamilies with contrasted modes of transposition. In line with earlier findings, the DNA transposons under study (TC1-Mariner) were found to transfer horizontally at the highest frequency, closely followed by the LTR superfamily (Copia), in contrast with the non-LTR superfamily (Jockey), that mostly diversifies through vertical inheritance and persists longer within genomes. Strikingly, across all superfamilies, we observe a marked excess of HTs in Lepidoptera, an insect order that also commonly hosts baculoviruses, known for their ability to transport host TEs. These results turn the spotlight on baculoviruses as major potential vectors of TEs in arthropods, and further emphasize the importance of non-vertical TE inheritance in genome evolution.

The Toxin-Antidote Model of Cytoplasmic Incompatibility: Genetics and Evolutionary Implications.

Wolbachia bacteria inhabit the cells of about half of all arthropod species, an unparalleled success stemming in large part from selfish invasive strategies. Cytoplasmic incompatibility (CI), whereby the symbiont makes itself essential to embryo viability, is the most common of these and constitutes a promising weapon against vector-borne diseases. After decades of theoretical and experimental struggle, major recent advances have been made toward a molecular understanding of this phenomenon. As pieces of the puzzle come together, from yeast and Drosophila fly transgenesis to CI diversity patterns in natural mosquito populations, it becomes clearer than ever that the CI induction and rescue stem from a toxin-antidote (TA) system. Further, the tight association of the CI genes with prophages provides clues to the possible evolutionary origin of this phenomenon and the levels of selection at play.

Comparing 16S rDNA amplicon sequencing and hybridization capture for pea aphid microbiota diversity analysis.

OBJECTIVE: Targeted sequencing of 16S rDNA amplicons is routinely used for microbial community profiling but this method suffers several limitations such as bias affinity of universal primers and short read size. Gene capture by hybridization represents a promising alternative. Here we used a metagenomic extract from the pea aphid Acyrthosiphon pisum to compare the performances of two widely used PCR primer pairs with DNA capture, based on solution hybrid selection. RESULTS: All methods produced an exhaustive description of the 8 bacterial taxa known to be present in this sample. In addition, the methods yielded similar quantitative results, with the number of reads strongly correlating with quantitative PCR controls. Both methods can thus be considered as qualitatively and quantitatively robust on such a sample with low microbial complexity.

How Long Does Wolbachia Remain on Board?

Wolbachia bacteria infect about half of all arthropods, with diverse and extreme consequences ranging from sex-ratio distortion and mating incompatibilities to protection against viruses. These phenotypic effects, combined with efficient vertical transmission from mothers to offspring, satisfactorily explain the invasion dynamics of Wolbachia within species. However, beyond the species level, the lack of congruence between the host and symbiont phylogenetic trees indicates that Wolbachia horizontal transfers and extinctions do happen and underlie its global distribution. But how often do they occur? And has the Wolbachia pandemic reached its equilibrium? Here, we address these questions by inferring recent acquisition/loss events from the distribution of Wolbachia lineages across the mitochondrial DNA tree of 3,600 arthropod specimens, spanning 1,100 species from Tahiti and surrounding islands. We show that most events occurred within the last million years, but are likely attributable to individual level variation (e.g., imperfect maternal transmission) rather than population level variation (e.g., Wolbachia extinction). At the population level, we estimate that mitochondria typically accumulate 4.7% substitutions per site during an infected episode, and 7.1% substitutions per site during the uninfected phase. Using a Bayesian time calibration of the mitochondrial tree, these numbers translate into infected and uninfected phases of approximately 7 and 9 million years. Infected species thus lose Wolbachia slightly more often than uninfected species acquire it, supporting the view that its present incidence, estimated here slightly below 0.5, represents an epidemiological equilibrium.

How and how much does RAD-seq bias genetic diversity estimates?

BACKGROUND: RAD-seq is a powerful tool, increasingly used in population genomics. However, earlier studies have raised red flags regarding possible biases associated with this technique. In particular, polymorphism on restriction sites results in preferential sampling of closely related haplotypes, so that RAD data tends to underestimate genetic diversity. RESULTS: Here we (1) clarify the theoretical basis of this bias, highlighting the potential confounding effects of population structure and selection, (2) confront predictions to real data from in silico digestion of full genomes and (3) provide a proof of concept toward an ABC-based correction of the RAD-seq bias. Under a neutral and panmictic model, we confirm the previously established relationship between the true polymorphism and its RAD-based estimation, showing a more pronounced bias when polymorphism is high. Using more elaborate models, we show that selection, resulting in heterogeneous levels of polymorphism along the genome, exacerbates the bias and leads to a more pronounced underestimation. On the contrary, spatial genetic structure tends to reduce the bias. We confront the neutral and panmictic model to "ideal" empirical data (in silico RAD-sequencing) using full genomes from natural populations of the fruit fly Drosophila melanogaster and the fungus Shizophyllum commune, harbouring respectively moderate and high genetic diversity. In D. melanogaster, predictions fit the model, but the small difference between the true and RAD polymorphism makes this comparison insensitive to deviations from the model. In the highly polymorphic fungus, the model captures a large part of the bias but makes inaccurate predictions. Accordingly, ABC corrections based on this model improve the estimations, albeit with some imprecisions. CONCLUSION: The RAD-seq underestimation of genetic diversity associated with polymorphism in restriction sites becomes more pronounced when polymorphism is high. In practice, this means that in many systems where polymorphism does not exceed 2 %, the bias is of minor importance in the face of other sources of uncertainty, such as heterogeneous bases composition or technical artefacts. The neutral panmictic model provides a practical mean to correct the bias through ABC, albeit with some imprecisions. More elaborate ABC methods might integrate additional parameters, such as population structure and selection, but their opposite effects could hinder accurate corrections.

Optimization of multiplexed RADseq libraries using low-cost adaptors.

Reduced representation genomics approaches, of which RADseq is currently the most popular form, offer the possibility to produce genome wide data from potentially any species, without previous genomic information. The application of RADseq to highly multiplexed libraries (including numerous specimens, and potentially numerous different species) is however limited by technical constraints. First, the cost of synthesis of Illumina adaptors including molecular identifiers (MIDs) becomes excessive when numerous specimens are to be multiplexed. Second, the necessity to empirically adjust the ratio of adaptors to genomic DNA concentration impedes the high throughput application of RADseq to heterogeneous samples, of variable DNA concentration and quality. In an attempt to solve these problems, we propose here some adjustments regarding the adaptor synthesis. First, we show that the common and unique (MID) parts of adaptors can be synthesized separately and subsequently ligated, which drastically reduces the synthesis cost, and thus allows multiplexing hundreds of specimens. Second, we show that self-ligation of adaptors, which makes the adaptor concentration so critical, can be simply prevented by using unphosphorylated adaptors, which significantly improves the ligation and sequencing yield.

The evolution of sex ratio distorter suppression affects a 25 cM genomic region in the butterfly Hypolimnas bolina.

Symbionts that distort their host's sex ratio by favouring the production and survival of females are common in arthropods. Their presence produces intense Fisherian selection to return the sex ratio to parity, typified by the rapid spread of host 'suppressor' loci that restore male survival/development. In this study, we investigated the genomic impact of a selective event of this kind in the butterfly Hypolimnas bolina. Through linkage mapping, we first identified a genomic region that was necessary for males to survive Wolbachia-induced male-killing. We then investigated the genomic impact of the rapid spread of suppression, which converted the Samoan population of this butterfly from a 100:1 female-biased sex ratio in 2001 to a 1:1 sex ratio by 2006. Models of this process revealed the potential for a chromosome-wide effect. To measure the impact of this episode of selection directly, the pattern of genetic variation before and after the spread of suppression was compared. Changes in allele frequencies were observed over a 25 cM region surrounding the suppressor locus, with a reduction in overall diversity observed at loci that co-segregate with the suppressor. These changes exceeded those expected from drift and occurred alongside the generation of linkage disequilibrium. The presence of novel allelic variants in 2006 suggests that the suppressor was likely to have been introduced via immigration rather than through de novo mutation. In addition, further sampling in 2010 indicated that many of the introduced variants were lost or had declined in frequency since 2006. We hypothesize that this loss may have resulted from a period of purifying selection, removing deleterious material that introgressed during the initial sweep. Our observations of the impact of suppression of sex ratio distorting activity reveal a very wide genomic imprint, reflecting its status as one of the strongest selective forces in nature.

Wolbachia divergence and the evolution of cytoplasmic incompatibility in Culex pipiens.

Many insect species harbor Wolbachia bacteria that induce cytoplasmic incompatibility (CI), i.e. embryonic lethality in crosses between infected males and uninfected females, or between males and females carrying incompatible Wolbachia strains. The molecular mechanism of CI remains unknown, but the available data are best interpreted under a modification-rescue model, where a mod function disables the reproductive success of infected males' sperm, unless the eggs are infected and express a compatible resc function. Here we examine the evolution of CI in the mosquito Culex pipiens, harbouring a large number of closely related Wolbachia strains structured in five distinct phylogenetic groups. Specifically, we used a worldwide sample of mosquito lines to assess the hypothesis that genetic divergence should correlate with the divergence of CI properties on a low evolutionary scale. We observed a significant association of Wolbachia genetic divergence with CI patterns. Most Wolbachia strains from the same group were compatible whereas those from different groups were often incompatible. Consistently, we found a strong association between Wolbachia groups and their mod-resc properties. Finally, lines from the same geographical area were rarely incompatible, confirming the conjecture that the spatial distribution of Wolbachia compatibility types should be constrained by selection. This study indicates a clear correlation between Wolbachia genotypes and CI properties, paving the way toward the identification of the molecular basis of CI through comparative genomics.

On the genetic architecture of cytoplasmic incompatibility: inference from phenotypic data.

Numerous insects carry intracellular bacteria that manipulate the insects' reproduction and thus facilitate their own spread. Cytoplasmic incompatibility (CI) is a common form of such manipulation, where a (currently uncharacterized) bacterial modification of male sperm induces the early death of embryos unless the fertilized eggs carry the same bacteria, inherited from the mother. The death of uninfected embryos provides an indirect selective advantage to infected ones, thus enabling the spread of the bacteria. Here we use and expand recently developed algorithms to infer the genetic architecture underlying the complex incompatibility data from the mosquito Culex pipiens. We show that CI requires more genetic determinants than previously believed and that quantitative variation in gene products potentially contributes to the observed CI patterns. In line with population genetic theory of CI, our analysis suggests that toxin factors (those inducing embryo death) are present in fewer copies in the bacterial genomes than antitoxin factors (those ensuring that infected embryos survive). In combination with comparative genomics, our approach will provide helpful guidance to identify the genetic basis of CI and more generally of other toxin/antitoxin systems that can be conceptualized under the same framework.

Is RAD-seq suitable for phylogenetic inference? An in silico assessment and optimization.

INFERRING PHYLOGENETIC RELATIONSHIPS BETWEEN CLOSELY RELATED TAXA CAN BE HINDERED BY THREE FACTORS: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD-seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within-species level will necessarily reduce the number of conserved and non-duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD-seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD-seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD-seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.

Draft genome sequence of the male-killing Wolbachia strain wBol1 reveals recent horizontal gene transfers from diverse sources.

BACKGROUND: The endosymbiont Wolbachia pipientis causes diverse and sometimes dramatic phenotypes in its invertebrate hosts. Four Wolbachia strains sequenced to date indicate that the constitution of the genome is dynamic, but these strains are quite divergent and do not allow resolution of genome diversification over shorter time periods. We have sequenced the genome of the strain wBol1-b, found in the butterfly Hypolimnas bolina, which kills the male offspring of infected hosts during embyronic development and is closely related to the non-male-killing strain wPip from Culex pipiens. RESULTS: The genomes of wBol1-b and wPip are similar in genomic organisation, sequence and gene content, but show substantial differences at some rapidly evolving regions of the genome, primarily associated with prophage and repetitive elements. We identified 44 genes in wBol1-b that do not have homologs in any previously sequenced strains, indicating that Wolbachia's non-core genome diversifies rapidly. These wBol1-b specific genes include a number that have been recently horizontally transferred from phylogenetically distant bacterial taxa. We further report a second possible case of horizontal gene transfer from a eukaryote into Wolbachia. CONCLUSIONS: Our analyses support the developing view that many endosymbiotic genomes are highly dynamic, and are exposed and receptive to exogenous genetic material from a wide range of sources. These data also suggest either that this bacterial species is particularly permissive for eukaryote-to-prokaryote gene transfers, or that these transfers may be more common than previously believed. The wBol1-b-specific genes we have identified provide candidates for further investigations of the genomic bases of phenotypic differences between closely-related Wolbachia strains.

Making (good) use of Wolbachia: what the models say.

Wolbachia, probably the most common animal endosymbiont, infects a wide range of arthropods as well as filarial nematodes. Generally vertically transmitted from mothers to offspring, it has evolved various strategies, ranging from brutal male-killing to mutualism, which facilitate invasion and persistence of the infections within host populations. Current interest in Wolbachia as a potential control agent against harmful nematodes and arthropods makes it important to be able to predict Wolbachia epidemiology and evolutionary trajectory. Here we highlight recent theoretical developments and suggest future modelling and empirical directions for basic and applied research in this domain.