Search

Hybridization and the consequent introgression of genomic elements is an important source of genetic diversity for biological lineages. This is particularly evident in young clades in which hybrid incompatibilities are still incomplete and mixing between species is more likely to occur. Drosophila paulistorum, a representative of the Neotropical Drosophila willistoni subgroup, is a classic model of incipient speciation. The species is divided into six semispecies that show varying degrees of pre-and post-mating incompatibility with each other. In the present study, we investigate the mitochondrial evolutionary history of D. paulistorum and the willistoni subgroup. For that, we perform phylogenetic and comparative analyses of the complete mitochondrial genomes and draft nuclear assemblies of 25 Drosophila lines of the willistoni and saltans species groups. Our results show that the mitochondria of D. paulistorum are polyphyletic and form two non-sister clades that we name alpha and beta. Identi-fication and analyses of nuclear mitochondrial insertions further reveal that the willistoni subgroup has an alpha-like mitochondrial ancestor and strongly suggest that both the alpha and beta mitochondria of D. paulistorum were acquired through introgression from unknown fly lineages of the willistoni subgroup. We also uncover multiple mito-chondrial introgressions across D. paulistorum semispecies and generate novel insight into the evolution of the species.

Wolbachia is a maternally inherited obligate endosymbiont that can induce a wide spectrum of effects in its host, ranging from mutualism to reproductive parasitism. At the genomic level, recombination within and between strains, transposable elements, and horizontal transfer of strains between host species make Wolbachia an evolutionarily dynamic bacterial system. The invasive cherry fruit fly Rhagoletis cingulata arrived in Europe from North America similar to 40 years ago, where it now co-occurs with the native cherry pest R. cerasi. This shared distribution has been proposed to have led to the horizontal transfer of different Wolbachia strains between the two species. To better understand transmission dynamics, we performed a comparative genome study of the strain wCin2 in its native United States and invasive European populations of R. cingulata with wCer2 in European R. cerasi. Previous multilocus sequence genotyping (MLST) of six genes implied that the source of wCer2 in R. cerasi was wCin2 from R. cingulata. However, we report genomic evidence discounting the recent horizontal transfer hypothesis for the origin of wCer2. Despite near identical sequences for the MLST markers, substantial sequence differences for other loci were found between wCer2 and wCin2, as well as structural rearrangements, and differences in prophage, repetitive element, gene content, and cytoplasmic incompatibility inducing genes. Our study highlights the need for whole-genome sequencing rather than relying on MLST markers for resolving Wolbachia strains and assessing their evolutionary dynamics.

Wolbachia is a widespread, vertically transmitted bacterial endosymbiont known for manipulating arthropod reproduction. Its most common form of reproductive manipulation is cytoplasmic incompatibility (CI), observed when a modification in the male sperm leads to embryonic lethality unless a compatible rescue factor is present in the female egg. CI attracts scientific attention due to its implications for host speciation and in the use of Wolbachia for controlling vector-borne diseases. However, our understanding of CI is complicated by the complexity of the phenotype, whose expression depends on both symbiont and host factors. In the present study, we perform a comparative analysis of nine complete Wolbachia genomes with known CI properties in the same genetic host background, Drosophila simulans STC. We describe genetic differences between closely related strains and uncover evidence that phages and other mobile elements contribute to the rapid evolution of both genomes and phenotypes of Wolbachia. Additionally, we identify both known and novel genes associated with the modification and rescue functions of CI. We combine our observations with published phenotypic information and discuss how variability in cif genes, novel CI-associated genes, and Wolbachia titer might contribute to poorly understood aspects of CI such as strength and bidirectional incompatibility. We speculate that high titer CI strains could be better at invading new hosts already infected with a CI Wolbachia, due to a higher rescue potential, and suggest that titer might thus be a relevant parameter to consider for future strategies using CI Wolbachia in biological control.

Cytoplasmic incompatibility is a selfish reproductive manipulation induced by the endosymbiont Wolbachia in arthropods. In males Wolbachia modifies sperm, leading to embryonic mortality in crosses with Wolbachia-free females. In females, Wolbachia rescues the cross and allows development to proceed normally. This provides a reproductive advantage to infected females, allowing the maternally transmitted symbiont to spread rapidly through host populations. We identified homologs of the genes underlying this phenotype, cifA and cifB, in 52 of 71 new and published Wolbachia genome sequences. They are strongly associated with cytoplasmic incompatibility. There are up to seven copies of the genes in each genome, and phylogenetic analysis shows that Wolbachia frequently acquires new copies due to pervasive horizontal transfer between strains. In many cases, the genes have subsequently acquired loss-of-function mutations to become pseudogenes. As predicted by theory, this tends to occur first in cifB, whose sole function is to modify sperm, and then in cifA, which is required to rescue the cross in females. Although cif genes recombine, recombination is largely restricted to closely related homologs. This is predicted under a model of coevolution between sperm modification and embryonic rescue, where recombination between distantly related pairs of genes would create a self-incompatible strain. Together, these patterns of gene gain, loss, and recombination support evolutionary models of cytoplasmic incompatibility.

The application of Wolbachia in insect pest and vector control requires the establishment of genotypically stable host associations. The cytoplasmic incompatibility (CI) inducing Wolbachia strain wCer2 naturally occurs in the cherry fruit fly Rhagoletis cerasi as co-infection with other strains and was transferred to other fruit fly species by embryonic microinjections. We obtained wCer2 genome data from its native and three novel hosts, Drosophila simulans, Drosophila melanogaster, and Ceratitis capitata and assessed its genome stability, characteristics, and CI factor (cit) genes. De novo assembly was successful from Wolbachia cell-enriched singly infected D. simulans embryos, with minimal host and other bacterial genome traces. The lowyield of Wolbachia sequence reads from total genomic extracts of one multiply infected R. cerasi pupa and one singly infected C. capitata adult limited de novo assemblies but was sufficient for comparative analyses. Across hosts wCer2 was stable in genome synteny and content. Polymorphic nucleotide sites were found in wCer2 of each host; however, only one nucleotide was different between R. cerasi and C capitata, and none between replicated D. simulans lines. The wCer2 genome is highly similar to wAu (D. simulans), wMel (D. melanogaster), and wRec (Drosophila recens). In contrast to wMel and wRec (each with one cif gene pair) and wAu (without any cif genes), wCer2 has three pairs of Type I cif genes, and one Type V cifB gene without a cifA complement. This may explain previously reported CI patterns of wCer2, including incomplete rescue of its own CI modification in three novel host species.

Background: The Neotropical fruit fly Drosophila paulistorum (Diptera: Drosophilidae) is a species complex in statu nascendi comprising six reproductively isolated semispecies, each harboring mutualistic Wolbachia strains. Although wild type flies of each semispecies are isolated from the others by both pre- and postmating incompatibilities, mating between semispecies and successful offspring development can be achieved once flies are treated with antibiotics to reduce Wolbachia titer. Here we use RNA-seq to study the impact of Wolbachia on D. paulistorum and investigate the hypothesis that the symbiont may play a role in host speciation. For that goal, we analyze samples of heads and abdomens of both sexes of the Amazonian, Centro American and Orinocan semispecies of D. paulistorum. Results: We identify between 175 and 1192 differentially expressed genes associated with a variety of biological processes that respond either globally or according to tissue, sex or condition in the three semispecies. Some of the functions associated with differentially expressed genes are known to be affected by Wolbachia in other species, such as metabolism and immunity, whereas others represent putative novel phenotypes involving muscular functions, pheromone signaling, and visual perception. Conclusions: Our results show that Wolbachia affect a large number of biological functions in D. paulistorum, particularly when present in high titer. We suggest that the significant metabolic impact of the infection on the host may cause several of the other putative and observed phenotypes. We also speculate that the observed differential expression of genes associated with chemical communication and reproduction may be associated with the emergence of pre- and postmating barriers between semispecies, which supports a role for Wolbachia in the speciation of D. paulistorum.

Background: Theory predicts that dependency within host-endosymbiont interactions results in endosymbiont genome size reduction. Unexpectedly, the largest Wolbachia genome was found in the obligate, parthenogenesis-associated wFol. In this study, we investigate possible processes underlying this genome expansion by comparing a re-annotated wFol genome to other Wolbachia genomes. In addition, we also search for candidate genes related to parthenogenesis induction (PI). Results: Within wFol, we found five phage WO regions representing 25.4% of the complete genome, few pseudogenized genes, and an expansion of DNA-repair genes in comparison to other Wolbachia. These signs of genome conservation were mirrored in the wFol host, the springtail F. candida, which also had an expanded DNA-repair gene family and many horizontally transferred genes. Across all Wolbachia genomes, there was a strong correlation between gene numbers of Wolbachia strains and their hosts. In order to identify genes with a potential link to PI, we assembled the genome of an additional PI strain, wLcla. Comparisons between four PI Wolbachia, including wFol and wLcla, and fourteen non-PI Wolbachia yielded a small set of potential candidate genes for further investigation. Conclusions: The strong similarities in genome content of wFol and its host, as well as the correlation between host and Wolbachia gene numbers suggest that there may be some form of convergent evolution between endosymbiont and host genomes. If such convergent evolution would be strong enough to overcome the evolutionary forces causing genome reduction, it would enable expanded genomes within long-term obligate endosymbionts.

In "Retrotransposons Are the Major Contributors to the Expansion of the Drosophila ananassae Muller F Element," Leung et al. (2017) improved contigs attributed to the Muller F element from the original CAF1 assembly, and used them to conclude that most of the sequence expansion of the fourth chromosome of D. ananassae is due to a higher transposon load than previously thought, but is not due to Wolbachia DNA integrations. While we do not disagree with the first conclusion, the authors base their second conclusion on the lack of homology detected between their improved CAF1 genome assembly attributed to D. ananassae and reference Wolbachia genomes. While the consensus CAF1 genome assembly lacks any sequence similarity to the reference genome of the Wolbachia endosymbiont of Drosophila melanogaster (wMel), numerous studies from multiple laboratories provide experimental support for a large lateral/horizontal gene transfer (LGT) of a Wolbachia genome into this D. ananassae line. As such, we strongly suspect that the original whole genome assembly was either constructed after the removal of all Wolbachia reads, or that Wolbachia sequences were directly removed from the contigs in the CAF1 assembly. Hence, Leung et al. (2017) could not have identified the Wolbachia LGT using the CAF1 assembly. This manuscript by Leung et al. (2017) highlights that an assembly of the Wolbachia sequence reads and their mate pairs was erroneously attributed solely to the Wolbachia endosymbiont, albeit before we understood the extent of LGT in D. ananassae. As such, we recommend that the sequences deposited at the National Center for Biotechnology Information (NCBI) under PRJNA13365 should not be attributed to Wolbachia endosymbiont of D. ananassae, but should have their taxonomy reclassified by NCBI as "Unclassified sequences." As our knowledge about genome biology improves, we need to reconsider and reanalyze earlier genomes removing the prejudice introduced from now defunct paradigms.

Genome reduction is a hallmark of symbiotic genomes, and the rate and patterns of gene loss associated with this process have been investigated in several different symbiotic systems. However, in long-term host-associated coevolving symbiont clades, the genome size differences between strains are normally quite small and hence patterns of large-scale genome reduction can only be inferred from distant relatives. Here we present the complete genome of a Coxiella-like symbiont from Rhipicephalus turanicus ticks (CRt), and compare it with other genomes from the genus Coxiella in order to investigate the process of genome reduction in a genus consisting of intracellular host-associated bacteria with variable genome sizes. The 1.7-Mb CRt genome is larger than the genomes of most obligate mutualists but has a very low protein-coding content (48.5%) and an extremely high number of identifiable pseudogenes, indicating that it is currently undergoing genome reduction. Analysis of encoded functions suggests that CRt is an obligate tick mutualist, as indicated by the possible provisioning of the tick with biotin (B7), riboflavin (B2) and other cofactors, and by the loss of most genes involved in host cell interactions, such as secretion systems. Comparative analyses between CRt and the 2.5 times smaller genome of Coxiella from the lonestar tick Amblyomma americanum (CLEAA) show that many of the same gene functions are lost and suggest that the large size difference might be due to a higher rate of genome evolution in CLEAA generated by the loss of the mismatch repair genes mutSL. Finally, sequence polymorphisms in the CRt population sampled from field collected ticks reveal up to one distinct strain variant per tick, and analyses of mutational patterns within the population suggest that selection might be acting on synonymous sites. The CRt genome is an extreme example of a symbiont genome caught in the act of genome reduction, and the comparison between CLEAA and CRt indicates

Background: Lateral gene transfer (LGT) from bacterial Wolbachia endosymbionts has been detected in similar to 20% of arthropod and nematode genome sequencing projects. Many of these transfers are large and contain a substantial part of the Wolbachia genome. Results: Here, we re-sequenced three D. ananassae genomes from Asia and the Pacific that contain large LGTs from Wolbachia. We find that multiple copies of the Wolbachia genome are transferred to the Drosophila nuclear genome in all three lines. In the D. ananassae line from Indonesia, the copies of Wolbachia DNA in the nuclear genome are nearly identical in size and sequence yielding an even coverage of mapped reads over the Wolbachia genome. In contrast, the D. ananassae lines from Hawaii and India show an uneven coverage of mapped reads over the Wolbachia genome suggesting that different parts of these LGTs are present in different copy numbers. In the Hawaii line, we find that this LGT is underrepresented in third instar larvae indicative of being heterochromatic. Fluorescence in situ hybridization of mitotic chromosomes confirms that the LGT in the Hawaii line is heterochromatic and represents similar to 20% of the sequence on chromosome 4 (dot chromosome, Muller element F). Conclusions: This collection of related lines contain large lateral gene transfers composed of multiple Wolbachia genomes that constitute >2% of the D. ananassae genome (similar to 5 Mbp) and partially explain the abnormally large size of chromosome 4 in D. ananassae.

The importance of host-specialization to speciation processes in obligate host-associated bacteria is well known, as is also the ability of recombination to generate cohesion in bacterial populations. However, whether divergent strains of highly recombining intracellular bacteria, such as Wolbachia, can maintain their genetic distinctness when infecting the same host is not known. We first developed a protocol for the genome sequencing of uncultivable endosymbionts. Using this method, we have sequenced the complete genomes of the Wolbachia strains wHa and wNo, which occur as natural double infections in Drosophila simulans populations on the Seychelles and in New Caledonia. Taxonomically, wHa belong to supergroup A and wNo to supergroup B. A comparative genomics study including additional strains supported the supergroup classification scheme and revealed 24 and 33 group-specific genes, putatively involved in host-adaptation processes. Recombination frequencies were high for strains of the same supergroup despite different host-preference patterns, leading to genomic cohesion. The inferred recombination fragments for strains of different supergroups were of short sizes, and the genomes of the co-infecting Wolbachia strains wHa and wNo were not more similar to each other and did not share more genes than other A- and B-group strains that infect different hosts. We conclude that Wolbachia strains of supergroup A and B represent genetically distinct clades, and that strains of different supergroups can co-exist in the same arthropod host without converging into the same species. This suggests that the supergroups are irreversibly separated and that barriers other than host-specialization are able to maintain distinct clades in recombining endosymbiont populations. Acquiring a good knowledge of the barriers to genetic exchange in Wolbachia will advance our understanding of how endosymbiont communities are constructed from vertically and horizontally transmitted genes.

Ankyrin repeat domain-encoding genes are common in the eukaryotic and viral domains of life, but they are rare in bacteria, the exception being a few obligate or facultative intracellular Proteobacteria species. Despite having a reduced genome, the arthropod strains of the alphaproteobacterium Wolbachia contain an unusually high number of ankyrin repeat domain-encoding genes ranging from 23 in wMel to 60 in wPip strain. This group of genes has attracted considerable attention for their astonishing large number as well as for the fact that ankyrin proteins are known to participate in protein-protein interactions, suggesting that they play a critical role in the molecular mechanism that determines host-Wolbachia symbiotic interactions. We present a comparative evolutionary analysis of the wMel-related ankyrin repeat domain-encoding genes present in different Drosophila-Wolbachia associations. Our results show that the ankyrin repeat domain-encoding genes change in size by expansion and contraction mediated by short directly repeated sequences. We provide examples of intragenic recombination events and show that these genes are likely to be horizontally transferred between strains with the aid of bacteriophages. These results confirm previous findings that the Wolbachia genomes are evolutionary mosaics and illustrate the potential that these bacteria have to generate diversity in proteins potentially involved in the symbiotic interactions.

The multiple displacement amplification method has revolutionized genomic studies of uncultured bacteria, where the extraction of pure DNA in sufficient quantity for next-generation sequencing is challenging. However, the method is problematic in that it amplifies the target DNA unevenly, induces the formation of chimeric reads and also amplifies contaminating DNA. Here, we have tested the reproducibility of the multiple displacement amplification method using serial dilutions of extracted genomic DNA and intact cells from the cultured endosymbiont Bartonella australis. The amplified DNA was sequenced with the Illumina sequencing technology, and the results were compared to sequence data obtained from unamplified DNA in this study as well as from a previously published genome project. We show that artifacts such as the extent of the amplification bias, the percentage of chimeric reads and the relative fraction of contaminating DNA increase dramatically for the smallest amounts of template DNA. The pattern of read coverage was reproducibly obtained for samples with higher amounts of template DNA, suggesting that the bias is non-random and genome-specific. A re-analysis of previously published sequence data obtained after amplification from clonal endosymbiont populations confirmed these predictions. We conclude that many of the artifacts associated with the use of the multiple displacement amplification method can be alleviated or much reduced by using multiple cells as the template for the amplification. These findings should be particularly useful for researchers studying the genomes of endosymbionts and other uncultured bacteria, for which a small clonal population of cells can be isolated.

BACKGROUND: The evolutionary importance of horizontal gene transfer (HGT) from Wolbachia endosymbiotic bacteria to their eukaryotic hosts is a topic of considerable interest and debate. Recent transfers of genome fragments from Wolbachia into insect chromosomes have been reported, but it has been argued that these fragments may be on an evolutionary trajectory to degradation and loss. RESULTS: We have discovered a case of HGT, involving two adjacent genes, between the genomes of Wolbachia and the currently Wolbachia-uninfected mosquito Aedes aegypti, an important human disease vector. The lower level of sequence identity between Wolbachia and insect, the transcription of all the genes involved, and the fact that we have identified homologs of the two genes in another Aedes species (Ae. mascarensis), suggest that these genes are being expressed after an extended evolutionary period since horizontal transfer, and therefore that the transfer has functional significance. The association of these genes with Wolbachia prophage regions also provides a mechanism for the transfer. CONCLUSION: The data support the argument that HGT between Wolbachia endosymbiotic bacteria and their hosts has produced evolutionary innovation.

BACKGROUND: Wolbachia are obligate endosymbiotic bacteria maternally transmitted through the egg cytoplasm that are responsible for several reproductive disorders in their insect hosts, such as cytoplasmic incompatibility (CI) in infected mosquitoes. Species in the Culex pipiens complex display an unusually high number of Wolbachia-induced crossing types, and based on present data, only the wPip strain is present. RESULTS: The sequencing of the wPip strain of Wolbachia revealed the presence of 60 ankyrin repeat domain (ANK) encoding genes and expression studies of these genes were carried out in adult mosquitoes. One of these ANK genes, pk2, is shown to be part of an operon of three prophage-associated genes with sex-specific expression, and is present in two identical copies in the genome. Another homolog of pk2 is also present that is differentially expressed in different Cx. pipiens group strains. A further two ANK genes showed sex-specific regulation in wPip-infected Cx. pipiens group adults. CONCLUSION: The high number, variability and differential expression of ANK genes in wPip suggest an important role in Wolbachia biology, and the gene family provides both markers and promising candidates for the study of reproductive manipulation.

We are interested in quantifying the contribution of gene acquisition, loss, expansion and rearrangements to the evolution of microbial genomes. Here, we discuss factors influencing microbial genome divergence based on pair-wise genome comparisons of closely related strains and species with different lifestyles. A particular focus is on intracellular pathogens and symbionts of the genera Rickettsia, Bartonella and BUCHNERA: Extensive gene loss and restricted access to phage and plasmid pools may provide an explanation for why single host pathogens are normally less successful than multihost pathogens. We note that species-specific genes tend to be shorter than orthologous genes, suggesting that a fraction of these may represent fossil-orfs, as also supported by multiple sequence alignments among species. The results of our genome comparisons are placed in the context of phylogenomic analyses of alpha and gamma proteobacteria. We highlight artefacts caused by different rates and patterns of mutations, suggesting that atypical phylogenetic placements can not a priori be taken as evidence for horizontal gene transfer events. The flexibility in genome structure among free-living microbes contrasts with the extreme stability observed for the small genomes of aphid endosymbionts, in which no rearrangements or inflow of genetic material have occurred during the past 50 millions years (1). Taken together, the results suggest that genomic stability correlate with the content of repeated sequences and mobile genetic elements, and thereby indirectly with bacterial lifestyles.

Comparative genome analyses of close relatives have yielded exciting insight into the sources of microbial genome variability with respect to gene content, gene order and evolution of genes with unknown functions. The genomes of free-living bacteria often carry phages and repetitive sequences that mediate genomic rearrangements in contrast to the small genomes of obligate host-associated bacteria. This suggests that genomic stability correlates with the genomic content of repeated sequences and movable genetic elements, and thereby with bacterial lifestyle. Genes with unknown functions present in a single species tend to be shorter than conserved, functional genes, indicating that the fraction of unique genes in microbial genomes has been overestimated.