The diversity and evolution of Wolbachia ankyrin repeat domain genes

Pan-Genome Analysis of Wolbachia, Endosymbiont of Diaphorina citri, Reveals Independent Origin in Asia and North America

Wolbachia, a group of Gram-negative symbiotic bacteria, infects nematodes and a wide range of arthropods. Diaphorina citri Kuwayama, the vector of Candidatus Liberibacter asiaticus (CLas) that causes citrus greening disease, is naturally infected with Wolbachia (wDi). However, the interaction between wDi and D. citri remains poorly understood. In this study, we performed a pan-genome analysis using 65 wDi genomes to gain a comprehensive understanding of wDi. Based on average nucleotide identity (ANI) analysis, we classified the wDi strains into Asia and North America strains. The ANI analysis, principal coordinates analysis (PCoA), and phylogenetic tree analysis supported that the D. citri in Florida did not originate from China. Furthermore, we found that a significant number of core genes were associated with metabolic pathways. Pathways such as thiamine metabolism, type I secretion system, biotin transport, and phospholipid transport were highly conserved across all analyzed wDi genomes. The variation analysis between Asia and North America wDi showed that there were 39,625 single-nucleotide polymorphisms (SNPs), 2153 indels, 10 inversions, 29 translocations, 65 duplications, 10 SV-based insertions, and 4 SV-based deletions. The SV-based insertions and deletions involved genes encoding transposase, phage tail tube protein, ankyrin repeat (ANK) protein, and group II intron-encoded protein. Pan-genome analysis of wDi contributes to our understanding of the geographical population of wDi, the origin of hosts of D. citri, and the interaction between wDi and its host, thus facilitating the development of strategies to control the insects and huanglongbing (HLB).

Wolbachia populations across organs of individual Culex pipiens: highly conserved intra-individual core pangenome with inter-individual polymorphisms

Wolbachia is a maternally inherited intracellular bacterium that infects a wide range of arthropods including mosquitoes. The endosymbiont is widely used in biocontrol strategies due to its capacity to modulate arthropod reproduction and limit pathogen transmission. Wolbachia infections in Culex spp. are generally assumed to be monoclonal but the potential presence of genetically distinct Wolbachia subpopulations within and between individual organs has not been investigated using whole genome sequencing. Here we reconstructed Wolbachia genomes from ovary and midgut metagenomes of single naturally infected Culex pipiens mosquitoes from Southern France to investigate patterns of intra- and inter-individual differences across mosquito organs. Our analyses revealed a remarkable degree of intra-individual conservancy among Wolbachia genomes from distinct organs of the same mosquito both at the level of gene presence-absence signal and single-nucleotide polymorphisms (SNPs). Yet, we identified several synonymous and non-synonymous substitutions between individuals, demonstrating the presence of some level of genomic heterogeneity among Wolbachia that infect the same C. pipiens field population. Overall, the absence of genetic heterogeneity within Wolbachia populations in a single individual confirms the presence of a dominant Wolbachia that is maintained under strong purifying forces of evolution.

'Candidatus Tisiphia' is a widespread Rickettsiaceae symbiont in the mosquito Anopheles plumbeus (Diptera: Culicidae)

Symbiotic bacteria can alter host biology by providing protection from natural enemies, or alter reproduction or vectoral competence. Symbiont-linked control of vector-borne disease in Anopheles has been hampered by a lack of symbioses that can establish stable vertical transmission in the host. Previous screening found the symbiont 'Candidatus Tisiphia' in Anopheles plumbeus, an aggressive biter and potential secondary vector of malaria parasites and West Nile virus. We screened samples collected over 10-years across Germany and used climate databases to assess environmental influence on incidence. We observed a 95% infection rate, and that the frequency of infection did not fluctuate with broad environmental factors. Maternal inheritance is indicated by presence in the ovaries through FISH microscopy. Finally, we assembled a high-quality 1.6 Mbp draft genome of 'Ca. Tisiphia' to explore its phylogeny and potential metabolic competence. The infection is closely related to strains found in Culicoides biting midges and shows similar patterns of metabolism, providing no evidence of the capacity to synthesize B-vitamins. This infection offers avenues for onward research in anopheline mosquito symbioses. Additionally, it provides future opportunity to study the impact of 'Ca. Tisiphia' on natural and transinfected hosts, especially in relation to reproductive fitness and vectorial competence and capacity.

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Male-killer symbiont screening reveals novel associations in Adalia ladybirds

While male-killing bacteria are known to infect across arthropods, ladybird beetles represent a hotspot for these symbioses. In some host species, there are multiple different symbionts that vary in presence and frequency between populations. To further our understanding of spatial and frequency variation, we tested for the presence of three male-killing bacteria: Wolbachia , Rickettsia and Spiroplasma , in two Adalia ladybird species from a previously unexplored UK population. The two-spot ladybird, A. bipunctata, is known to harbour all three male-killers, and we identified Spiroplasma infection in the Merseyside population for the first time. However, in contrast to previous studies on two-spot ladybirds from continental Europe, evidence from egg-hatch rates indicates the Spiroplasma strain present in the Merseyside population does not cause embryonic male-killing. In the related ten-spot ladybird, A. decempunctata, there is only one previous record of a male-killing symbiont, a Rickettsia , which we did not detect in the Merseyside sample. However, PCR assays indicated the presence of a Spiroplasma in a single A. decempunctata specimen. Marker sequence indicated that this Spiroplasma was divergent from that found in sympatric A. bipunctata. Genome sequencing of the Spiroplasma -infected A. decempunctata additionally revealed the presence of cobionts in the form of a Centistes parasitoid wasp and the parasitic fungi Beauveria. Further study of A. decempunctata from this population is needed to resolve whether it is the ladybird or wasp cobiont that harbours Spiroplasma , and to establish the phenotype of this strain. These data indicate first that microbial symbiont phenotype should not be assumed from past studies conducted in different locations, and second that cobiont presence may confound screening studies aimed to detect the frequency of a symbiont in field collected material from a focal host species.

Horizontally transferred genes as RNA interference targets for aphid and whitefly control

RNA interference (RNAi)-based technologies are starting to be commercialized as a new approach for agricultural pest control. Horizontally transferred genes (HTGs), which have been transferred into insect genomes from viruses, bacteria, fungi or plants, are attractive targets for RNAi-mediated pest control. HTGs are often unique to a specific insect family or even genus, making it unlikely that RNAi constructs targeting such genes will have negative effects on ladybugs, lacewings and other beneficial predatory insect species. In this study, we sequenced the genome of a red, tobacco-adapted isolate of Myzus persicae (green peach aphid) and bioinformatically identified 30 HTGs. We then used plant-mediated virus-induced gene silencing (VIGS) to show that several HTGs of bacterial and plant origin are important for aphid growth and/or survival. Silencing the expression of fungal-origin HTGs did not affect aphid survivorship but decreased aphid reproduction. Importantly, although there was uptake of plant-expressed RNA by Coccinella septempunctata (seven-spotted ladybugs) via the aphids that they consumed, we did not observe negative effects on ladybugs from aphid-targeted VIGS constructs. To demonstrate that this approach is more broadly applicable, we also targeted five Bemisia tabaci (whitefly) HTGs using VIGS and demonstrated that knockdown of some of these genes affected whitefly survival. As functional HTGs have been identified in the genomes of numerous pest species, we propose that these HTGs should be explored further as efficient and safe targets for control of insect pests using plant-mediated RNA interference.

One to host them all: genomics of the diverse bacterial endosymbionts of the spider Oedothorax gibbosus

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analysed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, 'Candidatus Tisiphia' (formerly Torix group Rickettsia), Cardinium and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome-resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, 'Ca. Tisiphia' and Cardinium endosymbionts. We provide insights into the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, 'Ca. Tisiphia' and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on Earth and show the usefulness of transposable elements as an evolutionary 'contact-tracing' tool.

Genomic and Phenotypic Comparisons Reveal Distinct Variants of Wolbachia Strain wAlbB

The intracellular bacterium Wolbachia inhibits virus replication and is being harnessed around the world to fight mosquito-borne diseases through releases of mosquitoes carrying the symbiont. Wolbachia strains vary in their ability to invade mosquito populations and suppress viruses in part due to differences in their density within the insect and associated fitness costs. Using whole-genome sequencing, we demonstrate the existence of two variants in wAlbB, a Wolbachia strain being released in natural populations of Aedes aegypti mosquitoes. The two variants display striking differences in genome architecture and gene content. Differences in the presence/absence of 52 genes between variants include genes located in prophage regions and others potentially involved in controlling the symbiont's density. Importantly, we show that these genetic differences correlate with variation in wAlbB density and its tolerance to heat stress, suggesting that different wAlbB variants may be better suited for field deployment depending on local environmental conditions. Finally, we found that the wAlbB genome remained stable following its introduction in a Malaysian mosquito population. Our results highlight the need for further genomic and phenotypic characterization of Wolbachia strains in order to inform ongoing Wolbachia-based programs and improve the selection of optimal strains in future field interventions. IMPORTANCE Dengue is a viral disease transmitted by Aedes mosquitoes that threatens around half of the world population. Recent advances in dengue control involve the introduction of Wolbachia bacterial symbionts with antiviral properties into mosquito populations, which can lead to dramatic decreases in the incidence of the disease. In light of these promising results, there is a crucial need to better understand the factors affecting the success of such strategies, in particular the choice of Wolbachia strain for field releases and the potential for evolutionary changes. Here, we characterized two variants of a Wolbachia strain used for dengue control that differ at the genomic level and in their ability to replicate within the mosquito. We also found no evidence for the evolution of the symbiont within the 2 years following its deployment in Malaysia. Our results have implications for current and future Wolbachia-based health interventions.

Widespread phages of endosymbionts: Phage WO genomics and the proposed taxonomic classification of Symbioviridae

Wolbachia are the most common obligate, intracellular bacteria in animals. They exist worldwide in arthropod and nematode hosts in which they commonly act as reproductive parasites or mutualists, respectively. Bacteriophage WO, the largest of Wolbachia’s mobile elements, includes reproductive parasitism genes, serves as a hotspot for genetic divergence and genomic rearrangement of the bacterial chromosome, and uniquely encodes a Eukaryotic Association Module with eukaryotic-like genes and an ensemble of putative host interaction genes. Despite WO’s relevance to genome evolution, selfish genetics, and symbiotic applications, relatively little is known about its origin, host range, diversification, and taxonomic classification. Here we analyze the most comprehensive set of 150 Wolbachia and phage WO assemblies to provide a framework for discretely organizing and naming integrated phage WO genomes. We demonstrate that WO is principally in arthropod Wolbachia with relatives in diverse endosymbionts and metagenomes, organized into four variants related by gene synteny, often oriented opposite the putative origin of replication in the Wolbachia chromosome, and the large serine recombinase is an ideal typing tool to distinguish the four variants. We identify a novel, putative lytic cassette and WO’s association with a conserved eleven gene island, termed Undecim Cluster, that is enriched with virulence-like genes. Finally, we evaluate WO-like Islands in the Wolbachia genome and discuss a new model in which Octomom, a notable WO-like Island, arose from a split with WO. Together, these findings establish the first comprehensive Linnaean taxonomic classification of endosymbiont phages, including non-Wolbachia phages from aquatic environments, that includes a new family and two new genera to capture the collective relatedness of these viruses.

Temperature effects on cellular host-microbe interactions explain continent-wide endosymbiont prevalence

Endosymbioses influence host physiology, reproduction, and fitness, but these relationships require efficient microbe transmission between host generations to persist. Maternally transmitted Wolbachia are the most common known endosymbionts,¹ but their frequencies vary widely within and among host populations for unknown reasons.^2,3 Here, we integrate genomic, cellular, and phenotypic analyses with mathematical models to provide an unexpectedly simple explanation for global wMel Wolbachia prevalence in Drosophila melanogaster. Cooling temperatures decrease wMel cellular abundance at a key stage of host oogenesis, producing temperature-dependent variation in maternal transmission that plausibly explains latitudinal clines of wMel frequencies on multiple continents. wMel sampled from a temperate climate targets the germline more efficiently in the cold than a recently differentiated tropical variant (∼2,200 years ago), indicative of rapid wMel adaptation to climate. Genomic analyses identify a very narrow list of wMel alleles-most notably, a derived stop codon in the major Wolbachia surface protein WspB-that underlie thermal sensitivity of cellular Wolbachia abundance and covary with temperature globally. Decoupling temperate wMel and host genomes further reduces transmission in the cold, a pattern that is characteristic of host-microbe co-adaptation to a temperate climate. Complex interactions among Wolbachia, hosts, and the environment (GxGxE) mediate wMel cellular abundance and maternal transmission, implicating temperature as a key determinant of Wolbachia spread and equilibrium frequencies, in conjunction with Wolbachia effects on host fitness and reproduction.^4,5 Our results motivate the strategic use of locally selected wMel variants for Wolbachia-based biocontrol efforts, which protect millions of individuals from arboviruses that cause human disease.⁶.

A Combined Analysis of Transcriptome and Proteome Reveals the Inhibitory Mechanism of a Novel Oligosaccharide Ester against Penicillium italicum

Blue mold caused by Penicillium italicum is one of the most serious postharvest diseases of citrus fruit. The aim of this study was to investigate the inhibitory effect of a novel oligosaccharide ester, 6-O-β-L-mannopyranosyl-3-O-(2-methylbutanoyl)-4-O-(8-methyldecanoyl)-2-O-(4-methyl-hexanoyl) trehalose (MTE-1), against P. italicum. Scanning electron microscopy (SEM) and transmission electron microscopy (TEM), along with transcriptome and proteome analysis also, were conducted to illuminate the underlying mechanism. Results showed that MTE-1 significantly inhibited P. italicum growth in vitro in a dose-dependent manner. Moreover, MTE-1 suppressed the disease development of citrus fruit inoculated with P. italicum. Furthermore, ultrastructure observation, as well as transcriptome and proteome analysis, indicated that MTE-1 treatment damaged the cell wall and plasma membrane in spores and mycelia of P. italicum. In addition, MTE-1 regulated genes or proteins involved in primary metabolism, cell-wall metabolism, and pathogenicity. These results demonstrate that MTE-1 inhibited P. italicum by damaging cell walls and membranes and disrupting normal cellular metabolism. These findings contribute to the understanding of the possible molecular action of MTE-1. Finally, MTE-1 also provides a new natural strategy for controlling diseases in postharvest fruit.

Label-free proteomic analysis reveals differentially expressed Wolbachia proteins in Tyrophagus putrescentiae: Mite allergens and markers reflecting population-related proteome differences

Tyrophagus putrescentiae is an astigmatid mite of great economic, medical and veterinary importance. The microbiome, especially intracellular bacteria, may affect allergy/allergen expression. We targeted Wolbachia proteins, allergen comparisons and markers in Wolbachia-mite interactions in three mite populations. A decoy database was constructed by proteogenomics using the T. putrescentiae draft genome, Wolbachia transcriptome assembly and current T. putrescentiae-related sequences in GenBank. Among thousands of mite-derived proteins, 18 Wolbachia proteins were reliably identified. We suggest that peroxiredoxin, bacterioferritin, ankyrin repeat domain-containing protein and DegQ family serine endoprotease indicate a higher-level bacterium-bacterium-host interaction. We produced evidence that the host-Wolbachia interaction is modulated through pattern recognition receptors (PRRs), mannose-binding lectins/mannose receptors, the cholinergic anti-inflammatory pathway with TNF-α, and others. We observed Tyr p 3 suppression in mites with Wolbachia, linking trypsin to PRR modulation. Nine out of the 12 current WHO/IUIS official allergens were reliably identified, but the remaining three allergens, Tyr p 1, 8 and 35, were detected as only trace hits. This study provides numerous markers for further Wolbachia-host interaction research. For accuracy, mite allergens should be considered according to abundance in species, but mite populations/strains, as well as their microbiome structure, may be key factors. SIGNIFICANCE: The astigmatid mites occurring in homes are significant producers of allergens that are highly dangerous to humans and domesticated animals. Mites are tightly associated with microorganisms that affect their biology and consequently allergy signatures. Mite populations were found to be infected with certain intracellular bacteria, but some populations lacked an intracellular bacterium. Our previous research showed that some populations of Tyrophagus putrescentiae are infected with Wolbachia, but some populations host additional bacteria of interest. Thus, there are not only interactions between the mites and Wolbachia but also likely an additional level of interaction that can be found in the interaction between different bacteria in the mites. These "higher-level" signatures and consequences that bacteria affect, including allergen production, are not understood in mites. In this study, we identified Wolbachia-specific proteins in mites for the first time. This study provides Wolbachia- and mite-derived markers that can be clues for describing "higher-level" mite-bacterium-bacterium interactions. Indeed, the microbiome contribution to allergies can potentially be derived directly from bacterial proteins, especially if they are abundant.

Isolation and characterization of a novel bacteriophage WO from Allonemobius socius crickets in Missouri

Wolbachia are endosymbionts of numerous arthropod and some nematode species, are important for their development and if present can cause distinct phenotypes of their hosts. Prophage DNA has been frequently detected in Wolbachia, but particles of Wolbachia bacteriophages (phage WO) have been only occasionally isolated. Here, we report the characterization and isolation of a phage WO of the southern ground cricket, Allonemobius socius, and provided the first whole-genome sequence of phage WO from this arthropod family outside of Asia. We screened A. socius abdomen DNA extracts from a cricket population in eastern Missouri by quantitative PCR for Wolbachia surface protein and phage WO capsid protein and found a prevalence of 55% and 50%, respectively, with many crickets positive for both. Immunohistochemistry using antibodies against Wolbachia surface protein showed many Wolbachia clusters in the reproductive system of female crickets. Whole-genome sequencing using Oxford Nanopore MinION and Illumina technology allowed for the assembly of a high-quality, 55 kb phage genome containing 63 open reading frames (ORF) encoding for phage WO structural proteins and host lysis and transcriptional manipulation. Taxonomically important regions of the assembled phage genome were validated by Sanger sequencing of PCR amplicons. Analysis of the nucleotides sequences of the ORFs encoding the large terminase subunit (ORF2) and minor capsid (ORF7) frequently used for phage WO phylogenetics showed highest homology to phage WOAu of Drosophila simulans (94.46% identity) and WOCin2USA1 of the cherry fruit fly, Rhagoletis cingulata (99.33% identity), respectively. Transmission electron microscopy examination of cricket ovaries showed a high density of phage particles within Wolbachia cells. Isolation of phage WO revealed particles characterized by 40–62 nm diameter heads and up to 190 nm long tails. This study provides the first detailed description and genomic characterization of phage WO from North America that is easily accessible in a widely distributed cricket species.

Transcriptional Response of Wolbachia to Dengue Virus Infection in Cells of the Mosquito Aedes aegypti

Aedes aegypti transmits one of the most significant mosquito-borne viruses, dengue virus (DENV). The absence of effective vaccines and clinical treatments and the emergence of insecticide resistance in A. aegypti necessitate novel vector control strategies. A new approach uses the endosymbiotic bacterium Wolbachia pipientis to reduce the spread of arboviruses. However, the Wolbachia-mediated antiviral mechanism is not well understood. To shed light on this mechanism, we investigated an unexplored aspect of Wolbachia-virus-mosquito interaction. We used RNA sequencing to examine the transcriptional response of Wolbachia to DENV infection in A. aegypti Aag2 cells transinfected with the wAlbB strain of Wolbachia. Our results suggest that genes encoding an endoribonuclease (RNase HI), a regulator of sigma 70-dependent gene transcription (6S RNA), essential cellular, transmembrane, and stress response functions and primary type I and IV secretion systems were upregulated, while a number of transport and binding proteins of Wolbachia, ribosome structure, and elongation factor-associated genes were downregulated due to DENV infection. Furthermore, bacterial retrotransposon, transposable, and phage-related elements were found among the up- and downregulated genes. We show that Wolbachia elicits a transcriptional response to virus infection and identify differentially expressed Wolbachia genes mostly at the early stages of virus infection. These findings highlight Wolbachia's ability to alter its gene expression in response to DENV infection of the host cell. IMPORTANCE Aedes aegypti is a vector of several pathogenic viruses, including dengue, Zika, chikungunya, and yellow fever viruses, which are of importance to human health. Wolbachia is an endosymbiotic bacterium currently used in transinfected mosquitoes to suppress replication and transmission of dengue viruses. However, the mechanism of Wolbachia-mediated virus inhibition is not fully understood. While several studies have shown mosquitoes' transcriptional responses to dengue virus infection, none have investigated these responses in Wolbachia, which may provide clues to the inhibition mechanism. Our results suggest changes in the expression of a number of functionally important Wolbachia genes upon dengue virus infection, including those involved in stress responses, providing insights into the endosymbiont's reaction to virus infection.

Genome Features of Asaia sp. W12 Isolated from the Mosquito Anopheles stephensi Reveal Symbiotic Traits

Asaia bacteria commonly comprise part of the microbiome of many mosquito species in the genera Anopheles and Aedes, including important vectors of infectious agents. Their close association with multiple organs and tissues of their mosquito hosts enhances the potential for paratransgenesis for the delivery of antimalaria or antivirus effectors. The molecular mechanisms involved in the interactions between Asaia and mosquito hosts, as well as Asaia and other bacterial members of the mosquito microbiome, remain underexplored. Here, we determined the genome sequence of Asaia strain W12 isolated from Anopheles stephensi mosquitoes, compared it to other Asaia species associated with plants or insects, and investigated the properties of the bacteria relevant to their symbiosis with mosquitoes. The assembled genome of strain W12 had a size of 3.94 MB, the largest among Asaia spp. studied so far. At least 3585 coding sequences were predicted. Insect-associated Asaia carried more glycoside hydrolase (GH)-encoding genes than those isolated from plants, showing their high plant biomass-degrading capacity in the insect gut. W12 had the most predicted regulatory protein components comparatively among the selected Asaia, indicating its capacity to adapt to frequent environmental changes in the mosquito gut. Two complete operons encoding cytochrome bo₃-type ubiquinol terminal oxidases (cyoABCD-1 and cyoABCD-2) were found in most Asaia genomes, possibly offering alternative terminal oxidases and allowing the flexible transition of respiratory pathways. Genes involved in the production of 2,3-butandiol and inositol have been found in Asaia sp. W12, possibly contributing to biofilm formation and stress tolerance.

Transmission of the wMel Wolbachia strain is modulated by its titre and by immune genes in Drosophila melanogaster (Wolbachia density and transmission)

Wolbachia are common intracellular endosymbionts of arthropods, but the interactions between Wolbachia and arthropods are only partially understood. The fruit fly Drosophila melanogaster is a model insect for understanding Wolbachia-host interactions. Here the native wMel strain of D. melanogaster was isolated and then different initial titres of wMel were artificially transferred back into antibiotics-treated fruit flies. Our purpose was to examine the interactions between the injected wMel in a density gradient and the recipient host during trans-generational transmission. The results showed that the trans-generational transmission rates of wMel and titres of wMel exhibited a fluctuating trend over nine generations, and the titres of wMel displayed a similar fluctuating trans-generational trend. There was a significant positive correlation between the transmission rate and the titre of wMel. Reciprocal crossings between wMel-transinfected and uninfected fruit flies revealed that wMel could induce cytoplasmic incompatibility (CI) at different initial titres, but the intensity of CI was not significantly correlated with the initial titre of wMel. Quantitative PCR analysis showed that the immune genes Drsl5 and Spn38F displayed a significant transcriptional response to wMel transfection, with an obvious negative correlation with the titre of wMel at the 3rd and 4th generations. Furthermore, RNA interference-mediated knockdown of Drsl5 and Spn38F elicited a drastic increase in the titre of wMel. In combination, our study suggests that the trans-generational transmission of wMel is modulated by its density, and the immune genes are involved in the regulation of Wolbachia density.

Diminutive, degraded but dissimilar: Wolbachia genomes from filarial nematodes do not conform to a single paradigm

Wolbachia are alpha-proteobacteria symbionts infecting a large range of arthropod species and two different families of nematodes. Interestingly, these endosymbionts are able to induce diverse phenotypes in their hosts: they are reproductive parasites within many arthropods, nutritional mutualists within some insects and obligate mutualists within their filarial nematode hosts. Defining Wolbachia 'species' is controversial and so they are commonly classified into 17 different phylogenetic lineages, termed supergroups, named A-F, H-Q and S. However, available genomic data remain limited and not representative of the full Wolbachia diversity; indeed, of the 24 complete genomes and 55 draft genomes of Wolbachia available to date, 84 % belong to supergroups A and B, exclusively composed of Wolbachia from arthropods. For the current study, we took advantage of a recently developed DNA-enrichment method to produce four complete genomes and two draft genomes of Wolbachia from filarial nematodes. Two complete genomes, wCtub and wDcau, are the smallest Wolbachia genomes sequenced to date (863 988 bp and 863 427 bp, respectively), as well as the first genomes representing supergroup J. These genomes confirm the validity of this supergroup, a controversial clade due to weaknesses of the multilocus sequence typing approach. We also produced the first draft Wolbachia genome from a supergroup F filarial nematode representative (wMhie), two genomes from supergroup D (wLsig and wLbra) and the complete genome of wDimm from supergroup C. Our new data confirm the paradigm of smaller Wolbachia genomes from filarial nematodes containing low levels of transposable elements and the absence of intact bacteriophage sequences, unlike many Wolbachia from arthropods, where both are more abundant. However, we observe differences among the Wolbachia genomes from filarial nematodes: no global co-evolutionary pattern, strong synteny between supergroup C and supergroup J Wolbachia, and more transposable elements observed in supergroup D Wolbachia compared to the other supergroups. Metabolic pathway analysis indicates several highly conserved pathways (haem and nucleotide biosynthesis, for example) as opposed to more variable pathways, such as vitamin B biosynthesis, which might be specific to certain host-symbiont associations. Overall, there appears to be no single Wolbachia-filarial nematode pattern of co-evolution or symbiotic relationship.

Transcriptomic Analysis of the Early Strobilar Development of Echinococcus granulosus

Echinococcus granulosus has a complex life cycle involving two mammalian hosts. The transition from one host to another is accompanied by changes in gene expression, and the transcriptional events that underlie this transition have not yet been fully characterized. In this study, RNA-seq was used to compare the transcription profiles of samples from E. granulosus protoscoleces induced in vitro to strobilar development at three time points. We identified 818 differentially expressed genes, which were divided into eight expression clusters formed over the entire 24 h period. An enrichment of gene transcripts with molecular functions of signal transduction, enzymes, and protein modifications was observed upon induction and developmental progression. This transcriptomic study provides insights for understanding the complex life cycle of E. granulosus and contributes for searching for the key genes correlating with the strobilar development, which can be used to identify potential candidates for the development of anthelmintic drugs.

Wolbachia Endosymbiont of the Horn Fly (Haematobia irritans irritans): a Supergroup A Strain with Multiple Horizontally Acquired Cytoplasmic Incompatibility Genes

The horn fly, Haematobia irritans irritans, is a hematophagous parasite of livestock distributed throughout Europe, Africa, Asia, and the Americas. Welfare losses on livestock due to horn fly infestation are estimated to cost between $1 billion and $2.5 billion (U.S. dollars) annually in North America and Brazil. The endosymbiotic bacterium Wolbachia pipientis is a maternally inherited manipulator of reproductive biology in arthropods and naturally infects laboratory colonies of horn flies from Kerrville, TX, and Alberta, Canada, but it has also been identified in wild-caught samples from Canada, the United States, Mexico, and Hungary. Reassembly of PacBio long-read and Illumina genomic DNA libraries from the Kerrville H. i. irritans genome project allowed for a complete and circularized 1.3-Mb Wolbachia genome (wIrr). Annotation of wIrr yielded 1,249 coding genes, 34 tRNAs, 3 rRNAs, and 5 prophage regions. Comparative genomics and whole-genome Bayesian evolutionary analysis of wIrr compared to published Wolbachia genomes suggested that wIrr is most closely related to and diverged from Wolbachia supergroup A strains known to infect Drosophila spp. Whole-genome synteny analyses between wIrr and closely related genomes indicated that wIrr has undergone significant genome rearrangements while maintaining high nucleotide identity. Comparative analysis of the cytoplasmic incompatibility (CI) genes of wIrr suggested two phylogenetically distinct CI loci and acquisition of another cifB homolog from phylogenetically distant supergroup A Wolbachia strains, suggesting horizontal acquisition of these loci. The wIrr genome provides a resource for future examination of the impact Wolbachia may have in both biocontrol and potential insecticide resistance of horn flies.IMPORTANCE Horn flies, Haematobia irritans irritans, are obligate hematophagous parasites of cattle having significant effects on production and animal welfare. Control of horn flies mainly relies on the use of insecticides, but issues with resistance have increased interest in development of alternative means of control. Wolbachia pipientis is an endosymbiont bacterium known to have a range of effects on host reproduction, such as induction of cytoplasmic incompatibility, feminization, male killing, and also impacts vector transmission. These characteristics of Wolbachia have been exploited in biological control approaches for a range of insect pests. Here we report the assembly and annotation of the circular genome of the Wolbachia strain of the Kerrville, TX, horn fly (wIrr). Annotation of wIrr suggests its unique features, including the horizontal acquisition of additional transcriptionally active cytoplasmic incompatibility loci. This study provides the foundation for future studies of Wolbachia-induced biological effects for control of horn flies.

Transgenic Testing Does Not Support a Role for Additional Candidate Genes in Wolbachia Male Killing or Cytoplasmic Incompatibility

Wolbachia are widespread bacterial endosymbionts that manipulate the reproduction of diverse arthropods to spread through a population and can substantially shape host evolution. Recently, reports identified three prophage WO genes (wmk, cifA, and cifB) that transgenically recapitulate many aspects of reproductive manipulation in Drosophila melanogaster. Here, we transgenically tested 10 additional gene candidates for CI and/or male killing in flies. The results yield no evidence for the involvement of these gene candidates in reproductive parasitism, bolstering the evidence for identification of the cif and wmk genes as the major factors involved in their phenotypes. In addition, evidence supports new hypotheses for prediction of male-killing phenotypes or lack thereof based on wmk transcript length and copy number. These experiments inform efforts to understand the full basis of reproductive parasitism for basic and applied purposes and lay the foundation for future work on the function of an interesting group of Wolbachia and phage WO genes.

Endosymbiotic bacteria in the genus Wolbachia remarkably infect nearly half of all arthropod species. They spread in part because of manipulations of host sexual reproduction that enhance the maternal transmission of the bacteria, including male killing (death of infected males) and unidirectional cytoplasmic incompatibility (CI; death of offspring from infected fathers and uninfected mothers). Recent discoveries identified several genes in prophage WO of Wolbachia (wmk, cifA, and cifB) that fully or partially recapitulate male killing or CI when transgenically expressed in Drosophila melanogaster. However, it is not yet fully resolved if other gene candidates contribute to these phenotypes. Here, we transgenically tested 10 additional gene candidates for their involvement in male killing and/or CI. The results show that despite sequence and protein architecture similarities or comparative associations with reproductive parasitism, transgenic expression of the candidates does not recapitulate male killing or CI. Sequence analysis across Wmk and its closest relatives reveals amino acids that may be important to its function. In addition, evidence is presented to propose new hypotheses regarding the relationship between wmk transcript length and its ability to kill a given host, as well as copy number of wmk homologs within a bacterial strain, which may be predictive of host resistance. Together, these analyses continue to build the evidence for identification of wmk, cifA, and cifB as the major genes that have thus far been shown to cause reproductive parasitism in Wolbachia, and the transgenic resources provide a basis for further functional study of phage WO genes.

IMPORTANCEWolbachia are widespread bacterial endosymbionts that manipulate the reproduction of diverse arthropods to spread through a population and can substantially shape host evolution. Recently, reports identified three prophage WO genes (wmk, cifA, and cifB) that transgenically recapitulate many aspects of reproductive manipulation in Drosophila melanogaster. Here, we transgenically tested 10 additional gene candidates for CI and/or male killing in flies. The results yield no evidence for the involvement of these gene candidates in reproductive parasitism, bolstering the evidence for identification of the cif and wmk genes as the major factors involved in their phenotypes. In addition, evidence supports new hypotheses for prediction of male-killing phenotypes or lack thereof based on wmk transcript length and copy number. These experiments inform efforts to understand the full basis of reproductive parasitism for basic and applied purposes and lay the foundation for future work on the function of an interesting group of Wolbachia and phage WO genes.

Trends in Symbiont-Induced Host Cellular Differentiation

Bacteria participate in a wide diversity of symbiotic associations with eukaryotic hosts that require precise interactions for bacterial recognition and persistence. Most commonly, host-associated bacteria interfere with host gene expression to modulate the immune response to the infection. However, many of these bacteria also interfere with host cellular differentiation pathways to create a hospitable niche, resulting in the formation of novel cell types, tissues, and organs. In both of these situations, bacterial symbionts must interact with eukaryotic regulatory pathways. Here, we detail what is known about how bacterial symbionts, from pathogens to mutualists, control host cellular differentiation across the central dogma, from epigenetic chromatin modifications, to transcription and mRNA processing, to translation and protein modifications. We identify four main trends from this survey. First, mechanisms for controlling host gene expression appear to evolve from symbionts co-opting cross-talk between host signaling pathways. Second, symbiont regulatory capacity is constrained by the processes that drive reductive genome evolution in host-associated bacteria. Third, the regulatory mechanisms symbionts exhibit correlate with the cost/benefit nature of the association. And, fourth, symbiont mechanisms for interacting with host genetic regulatory elements are not bound by native bacterial capabilities. Using this knowledge, we explore how the ubiquitous intracellular Wolbachia symbiont of arthropods and nematodes may modulate host cellular differentiation to manipulate host reproduction. Our survey of the literature on how infection alters gene expression in Wolbachia and its hosts revealed that, despite their intermediate-sized genomes, different strains appear capable of a wide diversity of regulatory manipulations. Given this and Wolbachia's diversity of phenotypes and eukaryotic-like proteins, we expect that many symbiont-induced host differentiation mechanisms will be discovered in this system.

Molecular host mimicry and manipulation in bacterial symbionts

It is common among intracellular bacterial pathogens to use eukaryotic-like proteins that mimic and manipulate host cellular processes to promote colonization and intracellular survival. Eukaryotic-like proteins are bacterial proteins with domains that are rare in bacteria, and known to function in the context of a eukaryotic cell. Such proteins can originate through horizontal gene transfer from eukaryotes or, in the case of simple repeat proteins, through convergent evolution. Recent studies of microbiomes associated with several eukaryotic hosts suggest that similar molecular strategies are deployed by cooperative bacteria that interact closely with eukaryotic cells. Some mimics, like ankyrin repeats, leucine rich repeats and tetratricopeptide repeats are shared across diverse symbiotic systems ranging from amoebae to plants, and may have originated early, or evolved independently in multiple systems. Others, like plant-mimicking domains in members of the plant microbiome are likely to be more recent innovations resulting from horizontal gene transfer from the host, or from microbial eukaryotes occupying the same host. Host protein mimics have only been described in a limited set of symbiotic systems, but are likely to be more widespread. Systematic searches for eukaryote-like proteins in symbiont genomes could lead to the discovery of novel mechanisms underlying host-symbiont interactions.

Wolbachia Horizontal Transmission Events in Ants: What Do We Know and What Can We Learn?

While strict vertical transmission insures the durability of intracellular symbioses, phylogenetic incongruences between hosts and endosymbionts suggest horizontal transmission must also occur. These horizontal acquisitions can have important implications for the biology of the host. Wolbachia is one of the most ecologically successful prokaryotes in arthropods, infecting an estimated 50–70% of all insect species. Much of this success is likely due to the fact that, in arthropods, Wolbachia is notorious for manipulating host reproduction to favor transmission through the female germline. However, its natural potential for horizontal transmission remains poorly understood. Here we evaluate the fundamental prerequisites for successful horizontal transfer, including necessary environmental conditions, genetic potential of bacterial strains, and means of mediating transfers. Furthermore, we revisit the relatedness of Wolbachia strains infecting the Panamanian leaf-cutting ant, Acromyrmex echinatior, and its inquiline social parasite, Acromyrmex insinuator, and compare our results to a study published more than 15 years ago by Van Borm et al. (2003). The results of this pilot study prompt us to reevaluate previous notions that obligate social parasitism reliably facilitates horizontal transfer and suggest that not all Wolbachia strains associated with ants have the same genetic potential for horizontal transmission.

Genome organisation and comparative genomics of four novel Wolbachia genome assemblies from Indian Drosophila host

Wolbachia has long been known to share an endosymbiotic relationship with its host as an obligate intracellular organism. Wolbachia diversity as different supergroups is found to be host-specific in most cases except a few, where the host species is seen to accommodate multiple strains. Besides, the Wolbachia genome must have undergone several changes in response to the evolving host genome in order to adapt and establish a strong association with its host, thus making a distinctive Wolbachia-host alliance. The present study focusses on four novel genome assembly and genome-wide sequence variations of Indian Wolbachia strains, i.e. wMel and wRi isolated from two different Drosophila hosts. The genome assembly has an average size of ~ 1.1 Mb and contains ~ 1100 genes, which is comparable with the previously sequenced Wolbachia genomes. The comparative genomics analysis of these genomes and sequence-wide comparison of some functionally significant genes, i.e. ankyrin repeats, Wsp and T4SS, highlight their sequence similarities and dissimilarities, further supporting the strain-specific association of Wolbachia to its host. Interestingly, some of the sequence variations are also found to be restricted to only Indian Wolbachia strains. Further analysis of prophage and their flanking regions in the Wolbachia genome reveals the presence of several functional genes which may assist the phage to reside inside the bacterial host, thus providing a trade-off for the endosymbiont-host association. Understanding this endosymbiont genome in different eco-geographical conditions has become imperative for the recent use of Wolbachia in medical entomology as a vector-control agent.

The draft genome of strain cCpun from biting midges confirms insect Cardinium are not a monophyletic group and reveals a novel gene family expansion in a symbiont

Background

It is estimated that 13% of arthropod species carry the heritable symbiont Cardinium hertigii. 16S rRNA and gyrB sequence divides this species into at least four groups (A–D), with the A group infecting a range of arthropods, the B group infecting nematode worms, the C group infecting Culicoides biting midges, and the D group associated with the marine copepod Nitocra spinipes. To date, genome sequence has only been available for strains from groups A and B, impeding general understanding of the evolutionary history of the radiation. We present a draft genome sequence for a C group Cardinium, motivated both by the paucity of genomic information outside of the A and B group, and the importance of Culicoides biting midge hosts as arbovirus vectors.

Methods

We reconstructed the genome of cCpun, a Cardinium strain from group C that naturally infects Culicoides punctatus, through Illumina sequencing of infected host specimens.

Results

The draft genome presented has high completeness, with BUSCO scores comparable to closed group A Cardinium genomes. Phylogenomic analysis based on concatenated single copy core proteins do not support Cardinium from arthropod hosts as a monophyletic group, with nematode Cardinium strains nested within the two groups infecting arthropod hosts. Analysis of the genome of cCpun revealed expansion of a variety of gene families classically considered important in symbiosis (e.g., ankyrin domain containing genes), and one set—characterized by DUF1703 domains—not previously associated with symbiotic lifestyle. This protein group encodes putative secreted nucleases, and the cCpun genome carried at least 25 widely divergent paralogs, 24 of which shared a common ancestor in the C group. The genome revealed no evidence in support of B vitamin provisioning to its haematophagous host, and indeed suggests Cardinium may be a net importer of biotin.

Discussion

These data indicate strains of Cardinium within nematodes cluster within Cardinium strains found in insects. The draft genome of cCpun further produces new hypotheses as to the interaction of the symbiont with the midge host, in particular the biological role of DUF1703 nuclease proteins that are predicted as being secreted by cCpun. In contrast, the coding content of this genome provides no support for a role for the symbiont in provisioning the host with B vitamins.

Genome expansion of an obligate parthenogenesis-associated Wolbachia poses an exception to the symbiont reduction model

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5492-9) contains supplementary material, which is available to authorized users.

Male-killing toxin in a bacterial symbiont of Drosophila

Several lineages of symbiotic bacteria in insects selfishly manipulate host reproduction to spread in a population ¹ , often by distorting host sex ratios. Spiroplasma poulsonii^2,3 is a helical and motile, Gram-positive symbiotic bacterium that resides in a wide range of Drosophila species ⁴ . A notable feature of S. poulsonii is male killing, whereby the sons of infected female hosts are selectively killed during development^1,2. Although male killing caused by S. poulsonii has been studied since the 1950s, its underlying mechanism is unknown. Here we identify an S. poulsonii protein, designated Spaid, whose expression induces male killing. Overexpression of Spaid in D. melanogaster kills males but not females, and induces massive apoptosis and neural defects, recapitulating the pathology observed in S. poulsonii-infected male embryos^5-11. Our data suggest that Spaid targets the dosage compensation machinery on the male X chromosome to mediate its effects. Spaid contains ankyrin repeats and a deubiquitinase domain, which are required for its subcellular localization and activity. Moreover, we found a laboratory mutant strain of S. poulsonii with reduced male-killing ability and a large deletion in the spaid locus. Our study has uncovered a bacterial protein that affects host cellular machinery in a sex-specific way, which is likely to be the long-searched-for factor responsible for S. poulsonii-induced male killing.

In Vitro Culture of the Insect Endosymbiont Spiroplasma poulsonii Highlights Bacterial Genes Involved in Host-Symbiont Interaction

Endosymbiotic bacteria associated with eukaryotic hosts are omnipresent in nature, particularly in insects. Studying the bacterial side of host-symbiont interactions is, however, often limited by the unculturability and genetic intractability of the symbionts. Spiroplasma poulsonii is a maternally transmitted bacterial endosymbiont that is naturally associated with several Drosophila species. S. poulsonii strongly affects its host’s physiology, for example by causing male killing or by protecting it against various parasites. Despite intense work on this model since the 1950s, attempts to cultivate endosymbiotic Spiroplasma in vitro have failed so far. Here, we developed a method to sustain the in vitro culture of S. poulsonii by optimizing a commercially accessible medium. We also provide a complete genome assembly, including the first sequence of a natural plasmid of an endosymbiotic Spiroplasma species. Last, by comparing the transcriptome of the in vitro culture to the transcriptome of bacteria extracted from the host, we identified genes putatively involved in host-symbiont interactions. This work provides new opportunities to study the physiology of endosymbiotic Spiroplasma and paves the way to dissect insect-endosymbiont interactions with two genetically tractable partners.

The discovery of insect bacterial endosymbionts (maternally transmitted bacteria) has revolutionized the study of insects, suggesting novel strategies for their control. Most endosymbionts are strongly dependent on their host to survive, making them uncultivable in artificial systems and genetically intractable. Spiroplasma poulsonii is an endosymbiont of Drosophila that affects host metabolism, reproduction, and defense against parasites. By providing the first reliable culture medium that allows a long-lasting in vitro culture of Spiroplasma and by elucidating its complete genome, this work lays the foundation for the development of genetic engineering tools to dissect endosymbiosis with two partners amenable to molecular study. Furthermore, the optimization method that we describe can be used on other yet uncultivable symbionts, opening new technical opportunities in the field of host-microbes interactions.

Culex pipiens crossing type diversity is governed by an amplified and polymorphic operon of Wolbachia

Culex pipiens mosquitoes are infected with Wolbachia (wPip) that cause an important diversity of cytoplasmic incompatibilities (CIs). Functional transgenic studies have implicated the cidA-cidB operon from wPip and its homolog in wMel in CI between infected Drosophila males and uninfected females. However, the genetic basis of the CI diversity induced by different Wolbachia strains was unknown. We show here that the remarkable diversity of CI in the C. pipiens complex is due to the presence, in all tested wPip genomes, of several copies of the cidA-cidB operon, which undergoes diversification through recombination events. In 183 isofemale lines of C. pipiens collected worldwide, specific variations of the cidA-cidB gene repertoires are found to match crossing types. The diversification of cidA-cidB is consistent with the hypothesis of a toxin–antitoxin system in which the gene cidB co-diversifies with the gene cidA, particularly in putative domains of reciprocal interactions.

Wolbachia causes cytoplasmic incompatibility (CI) between mosquitoes infected with different strains, but the genetic basis of observed CI diversity is unknown. Here, Bonneau et al. sequence Wolbachia from over 100 Culex pipiens lines and show that crossing types match variations of the toxin-antitoxin cidA-cidB genes.

The number of genes encoding repeat domain-containing proteins positively correlates with genome size in amoebal giant viruses

Curiously, in viruses, the virion volume appears to be predominantly driven by genome length rather than the number of proteins it encodes or geometric constraints. With their large genome and giant particle size, amoebal viruses (AVs) are ideally suited to study the relationship between genome and virion size and explore the role of genome plasticity in their evolutionary success. Different genomic regions of AVs exhibit distinct genealogies. Although the vertically transferred core genes and their functions are universally conserved across the nucleocytoplasmic large DNA virus (NCLDV) families and are essential for their replication, the horizontally acquired genes are variable across families and are lineage-specific. When compared with other giant virus families, we observed a near–linear increase in the number of genes encoding repeat domain-containing proteins (RDCPs) with the increase in the genome size of AVs. From what is known about the functions of RDCPs in bacteria and eukaryotes and their prevalence in the AV genomes, we envisage important roles for RDCPs in the life cycle of AVs, their genome expansion, and plasticity. This observation also supports the evolution of AVs from a smaller viral ancestor by the acquisition of diverse gene families from the environment including RDCPs that might have helped in host adaption.

Large-Scale Identification of Wolbachia pipientis Effectors

Wolbachia pipientis is an intracellular symbiont of arthropods well known for the reproductive manipulations induced in the host and, more recently, for the ability of Wolbachia to block virus replication in insect vectors. Since Wolbachia cannot yet be genetically manipulated, and due to the constraints imposed when working with an intracellular symbiont, little is known about mechanisms used by Wolbachia for host interaction. Here we employed a bioinformatics pipeline and identified 163 candidate effectors, potentially secreted by Wolbachia into the host cell. A total of 84 of these candidates were then subjected to a screen of growth defects induced in yeast upon heterologous expression which identified 14 top candidates likely secreted by Wolbachia. These predicted secreted effectors may function in concert as we find that their native expression is correlated and is highly upregulated at specific time points during Drosophila development. In addition, the evolutionary histories of some of these predicted effectors are also correlated, suggesting they may function together, or in the same pathway, during host infection. Similarly, most of these predicted effectors are limited to one or two Wolbachia strains—perhaps reflecting shared evolutionary history and strain specific functions in host manipulation. Identification of these Wolbachia candidate effectors is the first step in dissecting the mechanisms of symbiont–host interaction in this important system.

New paradigm in ankyrin repeats: Beyond protein-protein interaction module

Classically, ankyrin repeat (ANK) proteins are built from tandems of two or more repeats and form curved solenoid structures that are associated with protein-protein interactions. These are short, widespread structural motif of around 33 amino acids repeats in tandem, having a canonical helix-loop-helix fold, found individually or in combination with other domains. The multiplicity of structural pattern enables it to form assemblies of diverse sizes, required for their abilities to confer multiple binding and structural roles of proteins. Three-dimensional structures of these repeats determined to date reveal a degree of structural variability that translates into the considerable functional versatility of this protein superfamily. Recent work on the ANK has proposed novel structural information, especially protein-lipid, protein-sugar and protein-protein interaction. Self-assembly of these repeats was also shown to prevent the associated protein in forming filaments. In this review, we summarize the latest findings and how the new structural information has increased our understanding of the structural determinants of ANK proteins. We discussed latest findings on how these proteins participate in various interactions to diversify the ANK roles in numerous biological processes, and explored the emerging and evolving field of designer ankyrins and its framework for protein engineering emphasizing on biotechnological applications.

Insertion sequence polymorphism and genomic rearrangements uncover hidden Wolbachia diversity in Drosophila suzukii and D. subpulchrella

Ability to distinguish between closely related Wolbachia strains is crucial for understanding the evolution of Wolbachia-host interactions and the diversity of Wolbachia-induced phenotypes. A useful model to tackle these issues is the Drosophila suzukii - Wolbachia association. D. suzukii, a destructive insect pest, harbor a non-CI inducing Wolbachia 'wSuz' closely related to the strong CI-inducing wRi strain. Multi locus sequence typing (MLST) suggests presence of genetic homogeneity across wSuz strains infecting European and American D. suzukii populations, although different Wolbachia infection frequencies and host fecundity levels have been observed in both populations. Currently, it is not clear if these differences are due to cryptic wSuz polymorphism, host background, geographical factors or a combination of all of them. Here, we have identified geographical diversity in wSuz in D. suzukii populations from different continents using a highly diagnostic set of markers based on insertion sequence (IS) site polymorphism and genomic rearrangements (GR). We further identified inter-strain diversity between Wolbachia infecting D. suzukii and its sister species D. subpulchrella (wSpc). Based on our results, we speculate that discernible wSuz variants may associate with different observed host phenotypes, a hypothesis that demands future investigation. More generally, our results demonstrate the utility of IS and GRs in discriminating closely related Wolbachia strains.

Mi Casa es Su Casa: how an intracellular symbiont manipulates host biology

Wolbachia pipientis, the most common intracellular infection on the planet, infects 40% of insects as well as nematodes, isopods and arachnids. Wolbachia are obligately intracellular and challenging to study; there are no genetic tools for manipulating Wolbachia nor can they be cultured outside of host cells. Despite these roadblocks, the research community has defined a set of Wolbachia loci involved in host interaction: Wolbachia effectors. Through the use of Drosophila genetics, surrogate systems and biochemistry, the field has begun to define the toolkit Wolbachia use for host manipulation. Below we review recent findings identifying these Wolbachia effectors and point to potential, as yet uncharacterized, links between known phenotypes induced by Wolbachia infection and predicted effectors.

Population genomics of Wolbachia and mtDNA in Drosophila simulans from California

Wolbachia pipientis is an intracellular endosymbiont infecting many arthropods and filarial nematodes. Little is known about the short-term evolution of Wolbachia or its interaction with its host. Wolbachia is maternally inherited, resulting in co-inheritance of mitochondrial organelles such as mtDNA. Here I explore the evolution of Wolbachia, and the relationship between Wolbachia and mtDNA, using a large inbred panel of Drosophila simulans. I compare this to the only other large population genomic Wolbachia dataset from D. melanogaster. I find reduced diversity relative to expectation in both Wolbachia and mtDNA, but only mtDNA shows evidence of a recent selective sweep or population bottleneck. I estimate Wolbachia and mtDNA titre in each genotype, and I find considerable variation in both phenotypes, despite low genetic diversity in Wolbachia and mtDNA. A phylogeny of Wolbachia and of mtDNA suggest a recent origin of the infection derived from a single origin. Using Wolbachia and mtDNA titre as a phenotype, I perform the first association analysis using this phenotype with the nuclear genome and find several implicated regions, including one which contains four CAAX-box protein processing genes. CAAX-box protein processing can be an important part of host-pathogen interactions in other systems, suggesting interesting directions for future research.

Genomes of Candidatus Wolbachia bourtzisii wDacA and Candidatus Wolbachia pipientis wDacB from the Cochineal Insect Dactylopius coccus (Hemiptera: Dactylopiidae)

Dactylopius species, known as cochineal insects, are the source of the carminic acid dye used worldwide. The presence of two Wolbachia strains in Dactylopius coccus from Mexico was revealed by PCR amplification of wsp and sequencing of 16S rRNA genes. A metagenome analysis recovered the genome sequences of Candidatus Wolbachia bourtzisii wDacA (supergroup A) and Candidatus Wolbachia pipientis wDacB (supergroup B). Genome read coverage, as well as 16S rRNA clone sequencing, revealed that wDacB was more abundant than wDacA. The strains shared similar predicted metabolic capabilities that are common to Wolbachia, including riboflavin, ubiquinone, and heme biosynthesis, but lacked other vitamin and cofactor biosynthesis as well as glycolysis, the oxidative pentose phosphate pathway, and sugar uptake systems. A complete tricarboxylic acid cycle and gluconeogenesis were predicted as well as limited amino acid biosynthesis. Uptake and catabolism of proline were evidenced in Dactylopius Wolbachia strains. Both strains possessed WO-like phage regions and type I and type IV secretion systems. Several efflux systems found suggested the existence of metal toxicity within their host. Besides already described putative virulence factors like ankyrin domain proteins, VlrC homologs, and patatin-like proteins, putative novel virulence factors related to those found in intracellular pathogens like Legionella and Mycobacterium are highlighted for the first time in Wolbachia. Candidate genes identified in other Wolbachia that are likely involved in cytoplasmic incompatibility were found in wDacB but not in wDacA.

Genomic evidence for plant-parasitic nematodes as the earliest Wolbachia hosts

Wolbachia, one of the most widespread endosymbionts, is a target for biological control of mosquito-borne diseases (malaria and dengue virus), and antibiotic elimination of infectious filarial nematodes. We sequenced and analyzed the genome of a new Wolbachia strain (wPpe) in the plant-parasitic nematode Pratylenchus penetrans. Phylogenomic analyses placed wPpe as the earliest diverging Wolbachia, suggesting two evolutionary invasions into nematodes. The next branches comprised strains in sap-feeding insects, suggesting Wolbachia may have first evolved as a nutritional mutualist. Genome size, protein content, %GC, and repetitive DNA allied wPpe with mutualistic Wolbachia, whereas gene repertoire analyses placed it between parasite (A, B) and mutualist (C, D, F) groups. Conservation of iron metabolism genes across Wolbachia suggests iron homeostasis as a potential factor in its success. This study enhances our understanding of this globally pandemic endosymbiont, highlighting genetic patterns associated with host changes. Combined with future work on this strain, these genomic data could help provide potential new targets for plant-parasitic nematode control.

Evolution and function of eukaryotic-like proteins from sponge symbionts

Sponges (Porifera) are ancient metazoans that harbour diverse microorganisms, whose symbiotic interactions are essential for the host's health and function. Although symbiosis between bacteria and sponges are ubiquitous, the molecular mechanisms that control these associations are largely unknown. Recent (meta-) genomic analyses discovered an abundance of genes encoding for eukaryotic-like proteins (ELPs) in bacterial symbionts from different sponge species. ELPs belonging to the ankyrin repeat (AR) class from a bacterial symbiont of the sponge Cymbastela concentrica were subsequently found to modulate amoebal phagocytosis. This might be a molecular mechanism, by which symbionts can control their interaction with the sponge. In this study, we investigated the evolution and function of ELPs from other classes and from symbionts found in other sponges to better understand the importance of ELPs for bacteria-eukaryote interactions. Phylogenetic analyses showed that all of the nine ELPs investigated were most closely related to proteins found either in eukaryotes or in bacteria that can live in association with eukaryotes. ELPs were then recombinantly expressed in Escherichia coli and exposed to the amoeba Acanthamoeba castellanii, which is functionally analogous to phagocytic cells in sponges. Phagocytosis assays with E. coli containing three ELP classes (AR, TPR-SEL1 and NHL) showed a significantly higher percentage of amoeba containing bacteria and average number of intracellular bacteria per amoeba when compared to negative controls. The result that various classes of ELPs found in symbionts of different sponges can modulate phagocytosis indicates that they have a broader function in mediating bacteria-sponge interactions.

Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium dendrobatidis

Amphibian populations are experiencing catastrophic declines driven by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Although horizontal gene transfer (HGT) facilitates the evolution and adaptation in many fungi by conferring novel function genes to the recipient fungi, inter-kingdom HGT in Bd remains largely unexplored. In this study, our investigation detects 19 bacterial genes transferred to Bd, including metallo-beta-lactamase and arsenate reductase that play important roles in the resistance to antibiotics and arsenates. Moreover, three probable HGT gene families in Bd are from plants and one gene family coding the ankyrin repeat-containing protein appears to come from oomycetes. The observed multi-copy gene families associated with HGT are probably due to the independent transfer events or gene duplications. Five HGT genes with extracellular locations may relate to infection, and some other genes may participate in a variety of metabolic pathways, and in doing so add important metabolic traits to the recipient. The evolutionary analysis indicates that all the transferred genes evolved under purifying selection, suggesting that their functions in Bd are similar to those of the donors. Collectively, our results indicate that HGT from diverse donors may be an important evolutionary driver of Bd, and improve its adaptations for infecting and colonizing host amphibians.

Comparative Genomics of a Parthenogenesis-Inducing Wolbachia Symbiont

Wolbachia is an intracellular symbiont of invertebrates responsible for inducing a wide variety of phenotypes in its host. These host-Wolbachia relationships span the continuum from reproductive parasitism to obligate mutualism, and provide a unique system to study genomic changes associated with the evolution of symbiosis. We present the genome sequence from a parthenogenesis-inducing Wolbachia strain (wTpre) infecting the minute parasitoid wasp Trichogramma pretiosum. The wTpre genome is the most complete parthenogenesis-inducing Wolbachia genome available to date. We used comparative genomics across 16 Wolbachia strains, representing five supergroups, to identify a core Wolbachia genome of 496 sets of orthologous genes. Only 14 of these sets are unique to Wolbachia when compared to other bacteria from the Rickettsiales. We show that the B supergroup of Wolbachia, of which wTpre is a member, contains a significantly higher number of ankyrin repeat-containing genes than other supergroups. In the wTpre genome, there is evidence for truncation of the protein coding sequences in 20% of ORFs, mostly as a result of frameshift mutations. The wTpre strain represents a conversion from cytoplasmic incompatibility to a parthenogenesis-inducing lifestyle, and is required for reproduction in the Trichogramma host it infects. We hypothesize that the large number of coding frame truncations has accompanied the change in reproductive mode of the wTpre strain.

Host-Associated Genomic Features of the Novel Uncultured Intracellular Pathogen Ca. Ichthyocystis Revealed by Direct Sequencing of Epitheliocysts

Advances in single-cell and mini-metagenome sequencing have enabled important investigations into uncultured bacteria. In this study, we applied the mini-metagenome sequencing method to assemble genome drafts of the uncultured causative agents of epitheliocystis, an emerging infectious disease in the Mediterranean aquaculture species gilthead seabream. We sequenced multiple cyst samples and constructed 11 genome drafts from a novel beta-proteobacterial lineage, Candidatus Ichthyocystis. The draft genomes demonstrate features typical of pathogenic bacteria with an obligate intracellular lifestyle: a reduced genome of up to 2.6 Mb, reduced G + C content, and reduced metabolic capacity. Reconstruction of metabolic pathways reveals that Ca Ichthyocystis genomes lack all amino acid synthesis pathways, compelling them to scavenge from the fish host. All genomes encode type II, III, and IV secretion systems, a large repertoire of predicted effectors, and a type IV pilus. These are all considered to be virulence factors, required for adherence, invasion, and host manipulation. However, no evidence of lipopolysaccharide synthesis could be found. Beyond the core functions shared within the genus, alignments showed distinction into different species, characterized by alternative large gene families. These comprise up to a third of each genome, appear to have arisen through duplication and diversification, encode many effector proteins, and are seemingly critical for virulence. Thus, Ca Ichthyocystis represents a novel obligatory intracellular pathogenic beta-proteobacterial lineage. The methods used: mini-metagenome analysis and manual annotation, have generated important insights into the lifestyle and evolution of the novel, uncultured pathogens, elucidating many putative virulence factors including an unprecedented array of novel gene families.

Comparative Genomics of Two Closely Related Wolbachia with Different Reproductive Effects on Hosts

Wolbachia pipientis are obligate intracellular bacteria commonly found in many arthropods. They can induce various reproductive alterations in hosts, including cytoplasmic incompatibility, male-killing, feminization, and parthenogenetic development, and can provide host protection against some viruses and other pathogens. Wolbachia differ from many other primary endosymbionts in arthropods because they undergo frequent horizontal transmission between hosts and are well known for an abundance of mobile elements and relatively high recombination rates. Here, we compare the genomes of two closely related Wolbachia (with 0.57% genome-wide synonymous divergence) that differ in their reproductive effects on hosts. wVitA induces a sperm-egg incompatibility (also known as cytoplasmic incompatibility) in the parasitoid insect Nasonia vitripennis, whereas wUni causes parthenogenetic development in a different parasitoid, Muscidifurax uniraptor Although these bacteria are closely related, the genomic comparison reveals rampant rearrangements, protein truncations (particularly in proteins predicted to be secreted), and elevated substitution rates. These changes occur predominantly in the wUni lineage, and may be due in part to adaptations by wUni to a new host environment, or its phenotypic shift to parthenogenesis induction. However, we conclude that the approximately 8-fold elevated synonymous substitution rate in wUni is due to a either an elevated mutation rate or a greater number of generations per year in wUni, which occurs in semitropical host species. We identify a set of genes whose loss or pseudogenization in the wUni lineage implicates them in the phenotypic shift from cytoplasmic incompatibility to parthenogenesis induction. Finally, comparison of these closely related strains allows us to determine the fine-scale mutation patterns in Wolbachia Although Wolbachia are AT rich, mutation probabilities estimated from 4-fold degenerate sites are not AT biased, and predict an equilibrium AT content much less biased than observed (57-50% AT predicted vs. 76% current content at degenerate sites genome wide). The contrast suggests selection for increased AT content within Wolbachia genomes.

Comparative Genomics of Candidate Phylum TM6 Suggests That Parasitism Is Widespread and Ancestral in This Lineage

Candidate phylum TM6 is a major bacterial lineage recognized through culture-independent rRNA surveys to be low abundance members in a wide range of habitats; however, they are poorly characterized due to a lack of pure culture representatives. Two recent genomic studies of TM6 bacteria revealed small genomes and limited gene repertoire, consistent with known or inferred dependence on eukaryotic hosts for their metabolic needs. Here, we obtained additional near-complete genomes of TM6 populations from agricultural soil and upflow anaerobic sludge blanket reactor metagenomes which, together with the two publicly available TM6 genomes, represent seven distinct family level lineages in the TM6 phylum. Genome-based phylogenetic analysis confirms that TM6 is an independent phylum level lineage in the bacterial domain, possibly affiliated with the Patescibacteria superphylum. All seven genomes are small (1.0–1.5 Mb) and lack complete biosynthetic pathways for various essential cellular building blocks including amino acids, lipids, and nucleotides. These and other features identified in the TM6 genomes such as a degenerated cell envelope, ATP/ADP translocases for parasitizing host ATP pools, and protein motifs to facilitate eukaryotic host interactions indicate that parasitism is widespread in this phylum. Phylogenetic analysis of ATP/ADP translocase genes suggests that the ancestral TM6 lineage was also parasitic. We propose the name Dependentiae (phyl. nov.) to reflect dependence of TM6 bacteria on host organisms.

Dynamics of Wolbachia pipientis Gene Expression Across the Drosophila melanogaster Life Cycle

Symbiotic interactions between microbes and their multicellular hosts have manifold biological consequences. To better understand how bacteria maintain symbiotic associations with animal hosts, we analyzed genome-wide gene expression for the endosymbiotic α-proteobacteria Wolbachia pipientis across the entire life cycle of Drosophila melanogaster. We found that the majority of Wolbachia genes are expressed stably across the D. melanogaster life cycle, but that 7.8% of Wolbachia genes exhibit robust stage- or sex-specific expression differences when studied in the whole-organism context. Differentially-expressed Wolbachia genes are typically up-regulated after Drosophila embryogenesis and include many bacterial membrane, secretion system, and ankyrin repeat-containing proteins. Sex-biased genes are often organized as small operons of uncharacterized genes and are mainly up-regulated in adult Drosophila males in an age-dependent manner. We also systematically investigated expression levels of previously-reported candidate genes thought to be involved in host-microbe interaction, including those in the WO-A and WO-B prophages and in the Octomom region, which has been implicated in regulating bacterial titer and pathogenicity. Our work provides comprehensive insight into the developmental dynamics of gene expression for a widespread endosymbiont in its natural host context, and shows that public gene expression data harbor rich resources to probe the functional basis of the Wolbachia-Drosophila symbiosis and annotate the transcriptional outputs of the Wolbachia genome.

The Wolbachia WO bacteriophage proteome in the Aedes albopictus C/wStr1 cell line: evidence for lytic activity?

Wolbachia pipientis (Rickettsiales), an obligate intracellular alphaproteobacterium in insects, manipulates host reproduction to maximize invasion of uninfected insect populations. Modification of host population structure has potential applications for control of pest species, particularly if Wolbachia can be maintained, manipulated, and genetically engineered in vitro. Although Wolbachia maintains an obligate mutualism with genome stability in nematodes, arthropods can be co-infected with distinct Wolbachia strains, and horizontal gene transfer between strains is potentially mediated by WO phages encoded within Wolbachia genomes. Proteomic analysis of a robust, persistent infection of a mosquito cell line with wStr from the planthopper, Laodelphax striatellus, revealed expression of a full array of WO phage genes, as well as nine of ten non-phage genes that occur between two distinct clusters of WOMelB genes in the genome of wMel, which infects Drosophila melanogaster. These non-phage genes encode potential host-adaptive proteins and are expressed in wStr at higher levels than phage structural proteins. A subset of seven of the non-phage genes is flanked by highly conserved non-coding sequences, including a putative promoter element, that are not present in a syntenically arranged array of homologs in plasmids from three tick-associated Rickettsia spp. These studies expand our understanding of wStr in a host cell line derived from the mosquito, Aedes albopictus, and provide a basis for investigating conditions that favor the lytic phase of the WO phage life cycle and recovery of infectious phage particles.

Discovery of putative small non-coding RNAs from the obligate intracellular bacterium Wolbachia pipientis

Wolbachia pipientis is an endosymbiotic bacterium that induces a wide range of effects in its insect hosts, including manipulation of reproduction and protection against pathogens. Little is known of the molecular mechanisms underlying the insect-Wolbachia interaction, though it is likely to be mediated via the secretion of proteins or other factors. There is an increasing amount of evidence that bacteria regulate many cellular processes, including secretion of virulence factors, using small non-coding RNAs (sRNAs), but sRNAs have not previously been described from Wolbachia. We have used two independent approaches, one based on comparative genomics and the other using RNA-Seq data generated for gene expression studies, to identify candidate sRNAs in Wolbachia. We experimentally characterized the expression of one of these candidates in four Wolbachia strains, and showed that it is differentially regulated in different host tissues and sexes. Given the roles played by sRNAs in other host-associated bacteria, the conservation of the candidate sRNAs between different Wolbachia strains, and the sex- and tissue-specific differential regulation we have identified, we hypothesise that sRNAs may play a significant role in the biology of Wolbachia, and in particular in its interactions with its host.

Tandem-repeat protein domains across the tree of life

Tandem-repeat protein domains, composed of repeated units of conserved stretches of 20–40 amino acids, are required for a wide array of biological functions. Despite their diverse and fundamental functions, there has been no comprehensive assessment of their taxonomic distribution, incidence, and associations with organismal lifestyle and phylogeny. In this study, we assess for the first time the abundance of armadillo (ARM) and tetratricopeptide (TPR) repeat domains across all three domains in the tree of life and compare the results to our previous analysis on ankyrin (ANK) repeat domains in this journal. All eukaryotes and a majority of the bacterial and archaeal genomes analyzed have a minimum of one TPR and ARM repeat. In eukaryotes, the fraction of ARM-containing proteins is approximately double that of TPR and ANK-containing proteins, whereas bacteria and archaea are enriched in TPR-containing proteins relative to ARM- and ANK-containing proteins. We show in bacteria that phylogenetic history, rather than lifestyle or pathogenicity, is a predictor of TPR repeat domain abundance, while neither phylogenetic history nor lifestyle predicts ARM repeat domain abundance. Surprisingly, pathogenic bacteria were not enriched in TPR-containing proteins, which have been associated within virulence factors in certain species. Taken together, this comparative analysis provides a newly appreciated view of the prevalence and diversity of multiple types of tandem-repeat protein domains across the tree of life. A central finding of this analysis is that tandem repeat domain-containing proteins are prevalent not just in eukaryotes, but also in bacterial and archaeal species.