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

Variable effects of transient Wolbachia infections on alphaviruses in Aedes aegypti

Wolbachia pipientis (= Wolbachia) has promise as a tool to suppress virus transmission by Aedes aegypti mosquitoes. However, Wolbachia can have variable effects on mosquito-borne viruses. This variation remains poorly characterized, yet the multimodal effects of Wolbachia on diverse pathogens could have important implications for public health. Here, we examine the effects of transient somatic infection with two strains of Wolbachia (wAlbB and wMel) on the alphaviruses Sindbis virus (SINV), O'nyong-nyong virus (ONNV), and Mayaro virus (MAYV) in Ae. aegypti. We found variable effects of Wolbachia including enhancement and suppression of viral infections, with some effects depending on Wolbachia strain. Both wAlbB- and wMel-infected mosquitoes showed enhancement of SINV infection rates one week post-infection, with wAlbB-infected mosquitoes also having higher viral titers than controls. Infection rates with ONNV were low across all treatments and no significant effects of Wolbachia were observed. The effects of Wolbachia on MAYV infections were strikingly strain-specific; wMel strongly blocked MAYV infections and suppressed viral titers, while wAlbB had more modest effects. The variable effects of Wolbachia on vector competence underscore the importance of further research into how this bacterium impacts the virome of wild mosquitoes including the emergent human pathogens they transmit.

Comprehensive identification of GASA genes in sunflower and expression profiling in response to drought

Drought stress poses a critical threat to global crop yields and sustainable agriculture. The GASA genes are recognized for their pivotal role in stress tolerance and plant growth, but little is known about how they function in sunflowers. The investigation aimed to identify and elucidate the role of HaGASA genes in conferring sunflowers with drought tolerance. Twenty-seven different HaGASA gene family members were found in this study that were inconsistently located across eleven sunflower chromosomes. Phylogeny analysis revealed that the sunflower HaGASA genes were divided into five subgroups by comparing GASA genes with those from Arabidopsis, peanut, and soybean, with members within each subgroup displaying similar conserved motifs and gene structures. In-silico evaluation of cis-regulatory elements indicated the existence of specific elements associated with stress-responsiveness being the most abundant, followed by hormone, light, and growth-responsive elements. Transcriptomic data from the NCBI database was utilized to assess the HaGASA genes expression profile in different sunflower varieties under drought conditions. The HaGASA genes expression across ten sunflower genotypes under drought stress, revealed 14 differentially expressed HaGASA genes, implying their active role in the plant's stress response. The expression in different organs revealed that HaGASA2, HaGASA11, HaGASA17, HaGASA19, HaGASA21 and HaGASA26 displayed maximum expression in the stem. Our findings implicate HaGASA genes in mediating sunflower growth maintenance and adaptation to abiotic stress, particularly drought. The findings, taken together, provided a basic understanding of the structure and potential functions of HaGASA genes, setting the framework for further functional investigations into their roles in drought stress mitigation and crop improvement strategies.

Discovery of a novel Wolbachia in Heterodera expands nematode host distribution

Bioinformatics sequence data mining can reveal hidden microbial symbionts that might normally be filtered and removed as contaminants. Data mining can be helpful to detect Wolbachia, a widespread bacterial endosymbiont in insects and filarial nematodes whose distribution in plant-parasitic nematodes (PPNs) remains underexplored. To date, Wolbachia has only been reported a few PPNs, yet nematode-infecting Wolbachia may have been widespread in the evolutionary history of the phylum based on evidence of horizontal gene transfers, suggesting there may be undiscovered Wolbachia infections in PPNs. The goal of this study was to more broadly sample PPN Wolbachia strains in tylenchid nematodes to enable further comparative genomic analyses that may reveal Wolbachia's role and identify targets for biocontrol. Published whole-genome shotgun assemblies and their raw sequence data from 33 Meloidogyne spp. assemblies, seven Globodera spp. assemblies, and seven Heterodera spp. assemblies were analyzed to look for Wolbachia. No Wolbachia was found in Meloidogyne spp. and Globodera spp., but among seven genome assemblies for Heterodera spp., an H. schachtii assembly from the Netherlands was found to have a large Wolbachia-like sequence that, when re-assembled from reads, formed a complete, circular genome. Detailed analyses comparing read coverage, GC content, pseudogenes, and phylogenomic patterns clearly demonstrated that the H. schachtii Wolbachia represented a novel strain (hereafter, denoted wHet). Phylogenomic tree construction with PhyloBayes showed wHet was most closely related to another PPN Wolbachia, wTex, while 16S rRNA gene analysis showed it clustered with other Heterodera Wolbachia assembled from sequence databases. Pseudogenes in wHet suggested relatedness to the PPN clade, as did the lack of significantly enriched GO terms compared to PPN Wolbachia strains. It remains unclear whether the lack of Wolbachia in other published H. schachtii isolates represents the true absence of the endosymbiont from some hosts.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Bacterial community and genome analysis of cytoplasmic incompatibility-inducing Wolbachia in American serpentine leafminer, Liriomyza trifolii

Liriomyza trifolii, an agricultural pest, is occasionally infected by Wolbachia. A Wolbachia strain present in Liriomyza trifolii is associated with cytoplasmic incompatibility (CI) effects, leading to the death of embryos resulting from incompatible crosses between antibiotic-treated or naturally Wolbachia-free strain females and Wolbachia-infected males. In this study, high-throughput sequencing of hypervariable rRNA genes was employed to characterize the bacterial community in Wolbachia-infected L. trifolii without antibiotic treatment. The analysis revealed that Wolbachia dominates the bacterial community in L. trifolii, with minor presence of Acinetobacter, Pseudomonas, and Limnobacter. To elucidate the genetic basis of the CI phenotype, metagenomic sequencing was also conducted to assemble the genome of the Wolbachia strain. The draft-genome of the Wolbachia strain wLtri was 1.35 Mbp with 34% GC content and contained 1,487 predicted genes. Notably, within the wLtri genome, there are three distinct types of cytoplasmic incompatibility factor (cif) genes: Type I, Type III, and Type V cifA;B. These genes are likely responsible for inducing the strong cytoplasmic incompatibility observed in L. trifolii.

Intra-lineage microevolution of Wolbachia leads to the emergence of new cytoplasmic incompatibility patterns

Mosquitoes of the Culex pipiens complex are worldwide vectors of arbovirus, filarial nematodes, and avian malaria agents. In these hosts, the endosymbiotic bacteria Wolbachia induce cytoplasmic incompatibility (CI), i.e., reduced embryo viability in so-called incompatible crosses. Wolbachia infecting Culex pipiens (wPip) cause CI patterns of unparalleled complexity, associated with the amplification and diversification of cidA and cidB genes, with up to 6 different gene copies described in a single wPip genome. In wPip, CI is thought to function as a toxin-antidote (TA) system where compatibility relies on having the right antidotes (CidA) in the female to bind and neutralize the male's toxins (CidB). By repeating crosses between Culex isofemale lines over a 17 years period, we documented the emergence of a new compatibility type in real time and linked it to a change in cid genes genotype. We showed that loss of specific cidA gene copies in some wPip genomes results in a loss of compatibility. More precisely, we found that this lost antidote had an original sequence at its binding interface, corresponding to the original sequence at the toxin's binding interface. We showed that these original cid variants are recombinant, supporting a role for recombination rather than point mutations in rapid CI evolution. These results strongly support the TA model in natura, adding to all previous data acquired with transgenes expression.

Two male-killing Wolbachia from Drosophila birauraia that are closely related but distinct in genome structure

Insects harbour diverse maternally inherited bacteria and viruses, some of which have evolved to kill the male progeny of their hosts (male killing: MK). The fly species Drosophila biauraria carries a maternally transmitted MK-inducing partiti-like virus, but it was unknown if it carries other MK-inducing endosymbionts. Here, we identified two male-killing Wolbachia strains (wBiau1 and wBiau2) from D. biauraria and compared their genomes to elucidate their evolutionary processes. The two strains were genetically closely related but had exceptionally different genome structures with considerable rearrangements compared with combinations of other Wolbachia strains. Despite substantial changes in the genome structure, the two Wolbachia strains did not experience gene losses that would disrupt the male-killing expression or persistence in the host population. The two Wolbachia-infected matrilines carried distinct mitochondrial haplotypes, suggesting that wBiau1 and wBiau2 have invaded D. biauraria independently and undergone considerable genome changes owing to unknown selective pressures in evolutionary history. This study demonstrated the presence of three male-killers from two distinct origins in one fly species and highlighted the diverse and rapid genome evolution of MK Wolbachia in the host.

Two Complete Genomes of Male-Killing Wolbachia Infecting Ostrinia Moth Species Illuminate Their Evolutionary Dynamics and Association with Hosts

Wolbachia is an extremely widespread intracellular symbiont which causes reproductive manipulation on various arthropod hosts. Male progenies are killed in Wolbachia-infected lineages of the Japanese Ostrinia moth population. While the mechanism of male killing and the evolutionary interaction between host and symbiont are significant concerns for this system, the absence of Wolbachia genomic information has limited approaches to these issues. We determined the complete genome sequences of wFur and wSca, the male-killing Wolbachia of Ostrinia furnacalis and Ostrinia scapulalis. The two genomes shared an extremely high degree of homology, with over 95% of the predicted protein sequences being identical. A comparison of these two genomes revealed nearly minimal genome evolution, with a strong emphasis on the frequent genome rearrangements and the rapid evolution of ankyrin repeat-containing proteins. Additionally, we determined the mitochondrial genomes of both species' infected lineages and performed phylogenetic analyses to deduce the evolutionary dynamics of Wolbachia infection in the Ostrinia clade. According to the inferred phylogenetic relationship, two possible scenarios were proposed: (1) Wolbachia infection was established in the Ostrinia clade prior to the speciation of related species such as O. furnacalis and O. scapulalis, or (2) Wolbachia infection in these species was introgressively transferred from a currently unidentified relative. Simultaneously, the relatively high homology of mitochondrial genomes suggested recent Wolbachia introgression between infected Ostrinia species. The findings of this study collectively shed light on the host-symbiont interaction from an evolutionary standpoint.

Variable effects of Wolbachia on alphavirus infection in Aedes aegypti

Wolbachia pipientis (=Wolbachia) has promise as a tool to suppress virus transmission by Aedes aegypti mosquitoes. However, Wolbachia can have variable effects on mosquito-borne viruses. This variation remains poorly characterized, yet the multimodal effects of Wolbachia on diverse pathogens could have important implications for public health. Here, we examine the effects of somatic infection with two strains of Wolbachia (wAlbB and wMel) on the alphaviruses Sindbis virus (SINV), O'nyong-nyong virus (ONNV), and Mayaro virus (MAYV) in Ae. aegypti. We found variable effects of Wolbachia including enhancement and suppression of viral infections, with some effects depending on Wolbachia strain. Both wAlbB- and wMel-infected mosquitoes showed enhancement of SINV infection rates one week post-infection, with wAlbB-infected mosquitoes also having higher viral titers than controls. Infection rates with ONNV were low across all treatments and no significant effects of Wolbachia were observed. The effects of Wolbachia on MAYV infections were strikingly strain-specific; wMel strongly blocked MAYV infections and suppressed viral titers, while wAlbB did not influence MAYV infection. The variable effects of Wolbachia on vector competence underscore the importance of further research into how this bacterium impacts the virome of wild mosquitoes including the emergent human pathogens they transmit.

New insights into the genome and transmission of the microsporidian pathogen Nosema muscidifuracis

Introduction

Nosema is a diverse genus of unicellular microsporidian parasites of insects and other arthropods. Nosema muscidifuracis infects parasitoid wasp species of Muscidifurax zaraptor and M. raptor (Hymenoptera: Pteromalidae), causing ~50% reduction in longevity and ~90% reduction in fecundity.

Methods and Results

Here, we report the first assembly of the N. muscidifuracis genome (14,397,169 bp in 28 contigs) of high continuity (contig N50 544.3 Kb) and completeness (BUSCO score 97.0%). A total of 2,782 protein-coding genes were annotated, with 66.2% of the genes having two copies and 24.0% of genes having three copies. These duplicated genes are highly similar, with a sequence identity of 99.3%. The complex pattern suggests extensive gene duplications and rearrangements across the genome. We annotated 57 rDNA loci, which are highly GC-rich (37%) in a GC-poor genome (25% genome average). Nosema-specific qPCR primer sets were designed based on 18S rDNA annotation as a diagnostic tool to determine its titer in host samples. We discovered high Nosema titers in Nosema-cured M. raptor and M. zaraptor using heat treatment in 2017 and 2019, suggesting that the remedy did not completely eliminate the Nosema infection. Cytogenetic analyses revealed heavy infections of N. muscidifuracis within the ovaries of M. raptor and M. zaraptor, consistent with the titer determined by qPCR and suggesting a heritable component of infection and per ovum vertical transmission.

Discussion

The parasitoids-Nosema system is laboratory tractable and, therefore, can serve as a model to inform future genome manipulations of Nosema-host system for investigations of Nosemosis.

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.

Complete de novo assembly of Wolbachia endosymbiont of Diaphorina citri Kuwayama (Hemiptera: Liviidae) using long-read genome sequencing

Wolbachia, a gram-negative \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{\alpha }$$\end{document}α-proteobacterium, is an endosymbiont found in some arthropods and nematodes. Diaphorina citri Kuwayama, the vector of ‘Candidatus Liberibacter asiaticus’ (CLas), are naturally infected with a strain of Wolbachia (wDi), which has been shown to colocalize with the bacteria pathogens CLas, the pathogen associated with huanglongbing (HLB) disease of citrus. The relationship between wDi and CLas is poorly understood in part because the complete genome of wDi has not been available. Using high-quality long-read PacBio circular consensus sequences, we present the largest complete circular wDi genome among supergroup-B members. The assembled circular chromosome is 1.52 megabases with 95.7% genome completeness with contamination of 1.45%, as assessed by checkM. We identified Insertion Sequences (ISs) and prophage genes scattered throughout the genomes. The proteins were annotated using Pfam, eggNOG, and COG that assigned unique domains and functions. The wDi genome was compared with previously sequenced Wolbachia genomes using pangenome and phylogenetic analyses. The availability of a complete circular chromosome of wDi will facilitate understanding of its role within the insect vector, which may assist in developing tools for disease management. This information also provides a baseline for understanding phylogenetic relationships among Wolbachia of other insect vectors.

Long-Read Assembly and Annotation of the Parasitoid Wasp Muscidifurax raptorellus, a Biological Control Agent for Filth Flies

The parasitoid wasp Muscidifurax raptorellus (Hymenoptera: Pteromalidae) is a gregarious species that has received extensive attention for its potential in biological pest control against house fly, stable fly, and other filth flies. It has a high reproductive capacity and can be reared easily. However, genome assembly is not available for M. raptorellus or any other species in this genus. Previously, we assembled a complete circular mitochondrial genome with a length of 24,717 bp. Here, we assembled and annotated a high-quality nuclear genome of M. raptorellus, using a combination of long-read (104× genome coverage) and short-read (326× genome coverage) sequencing technologies. The assembled genome size is 314 Mbp in 226 contigs, with a 97.9% BUSCO completeness score and a contig N50 of 4.67 Mb, suggesting excellent continuity of this assembly. Our assembly builds the foundation for comparative and evolutionary genomic analysis in the genus of Muscidifurax and possible future biocontrol applications.

Parallel Sequencing of Wolbachia wCer2 from Donor and Novel Hosts Reveals Multiple Incompatibility Factors and Genome Stability after Host Transfers

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 (cif) genes. De novo assembly was successful from Wolbachia cell-enriched singly infected D. simulans embryos, with minimal host and other bacterial genome traces. The low yield 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.

Living in the endosymbiotic world of Wolbachia: A centennial review

The most widespread intracellular bacteria in the animal kingdom are maternally inherited endosymbionts of the genus Wolbachia. Their prevalence in arthropods and nematodes worldwide and stunning arsenal of parasitic and mutualistic adaptations make these bacteria a biological archetype for basic studies of symbiosis and applied outcomes for curbing human and agricultural diseases. Here, we conduct a summative, centennial analysis of living in the Wolbachia world. We synthesize literature on Wolbachia's host range, phylogenetic diversity, genomics, cell biology, and applications to filarial, arboviral, and agricultural diseases. We also review the mobilome of Wolbachia including phage WO and its essentiality to hallmark reproductive phenotypes in arthropods. Finally, the Wolbachia system is an exemplar for discovery-based science education using biodiversity, biotechnology, and bioinformatics lessons. As we approach a century of Wolbachia research, the interdisciplinary science of this symbiosis stands as a model for consolidating and teaching the integrative rules of endosymbiotic life.

Living in the endosymbiotic world of Wolbachia: A centennial review

The most widespread intracellular bacteria in the animal kingdom are maternally inherited endosymbionts of the genus Wolbachia. Their prevalence in arthropods and nematodes worldwide and stunning arsenal of parasitic and mutualistic adaptations make these bacteria a biological archetype for basic studies of symbiosis and applied outcomes for curbing human and agricultural diseases. Here, we conduct a summative, centennial analysis of living in the Wolbachia world. We synthesize literature on Wolbachia's host range, phylogenetic diversity, genomics, cell biology, and applications to filarial, arboviral, and agricultural diseases. We also review the mobilome of Wolbachia including phage WO and its essentiality to hallmark reproductive phenotypes in arthropods. Finally, the Wolbachia system is an exemplar for discovery-based science education using biodiversity, biotechnology, and bioinformatics lessons. As we approach a century of Wolbachia research, the interdisciplinary science of this symbiosis stands as a model for consolidating and teaching the integrative rules of endosymbiotic life.

Comparative Genomics Reveals Factors Associated with Phenotypic Expression of Wolbachia

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.

Shrinking of repeating unit length in leucine-rich repeats from double-stranded DNA viruses

Leucine-rich repeats (LRRs) are present in over 563,000 proteins from viruses to eukaryotes. LRRs repeat in tandem and have been classified into fifteen classes in which the repeat unit lengths range from 20 to 29 residues. Most LRR proteins are involved in protein-protein or ligand interactions. The amount of genome sequence data from viruses is increasing rapidly, and although viral LRR proteins have been identified, a comprehensive sequence analysis has not yet been done, and their structures, functions, and evolution are still unknown. In the present study, we characterized viral LRRs by sequence analysis and identified over 600 LRR proteins from 89 virus species. Most of these proteins were from double-stranded DNA (dsDNA) viruses, including nucleocytoplasmic large dsDNA viruses (NCLDVs). We found that the repeating unit lengths of 11 types are one to five residues shorter than those of the seven known corresponding LRR classes. The repeating units of six types are 19 residues long and are thus the shortest among all LRRs. In addition, two of the LRR types are unique and have not been observed in bacteria, archae or eukaryotes. Conserved strongly hydrophobic residues such as Leu, Val or Ile in the consensus sequences are replaced by Cys with high frequency. Phylogenetic analysis indicated that horizontal gene transfer of some viral LRR genes had occurred between the virus and its host. We suggest that the shortening might contribute to the survival strategy of viruses. The present findings provide a new perspective on the origin and evolution of LRRs.

Mycoavidus sp. Strain B2-EB: Comparative Genomics Reveals Minimal Genomic Features Required by a Cultivable Burkholderiaceae-Related Endofungal Bacterium

Obligate bacterial endosymbionts are critical to the existence of many eukaryotes. Such endobacteria are usually characterized by reduced genomes and metabolic dependence on the host, which may cause difficulty in isolating them in pure cultures. Family Burkholderiaceae-related endofungal bacteria affiliated with the Mycoavidus-Glomeribacter clade can be associated with the fungal subphyla Mortierellomycotina and Glomeromycotina. In this study, a cultivable endosymbiotic bacterium, Mycoavidus sp. strain B2-EB, present in the fungal host Mortierella parvispora was obtained successfully. The B2-EB genome (1.88 Mb) represents the smallest genome among the endofungal bacterium Mycoavidus cysteinexigens (2.64-2.80 Mb) of Mortierella elongata and the uncultured endosymbiont "Candidatus Glomeribacter gigasporarum" (1.37 to 2.36 Mb) of arbuscular mycorrhizal fungi. Despite a reduction in genome size, strain B2-EB displays a high genome completeness, suggesting a nondegenerative reduction in the B2-EB genome. Compared with a large proportion of transposable elements (TEs) in other known Mycoavidus genomes (7.2 to 11.5% of the total genome length), TEs accounted for only 2.4% of the B2-EB genome. This pattern, together with a high proportion of single-copy genes in the B2-EB genome, suggests that the B2-EB genome reached a state of relative evolutionary stability. These results represent the most streamlined structure among the cultivable endofungal bacteria and suggest the minimal genome features required by both an endofungal lifestyle and artificial culture. This study allows us to understand the genome evolution of Burkholderiaceae-related endosymbionts and to elucidate microbiological interactions.IMPORTANCE This study attempted the isolation of a novel endobacterium, Mycoavidus sp. B2-EB (JCM 33615), harbored in the fungal host Mortierella parvispora E1425 (JCM 39028). We report the complete genome sequence of this strain, which possesses a reduced genome size with relatively high genome completeness and a streamlined genome structure. The information indicates the minimal genomic features required by both the endofungal lifestyle and artificial cultivation, which furthers our understanding of genome reduction in fungal endosymbionts and extends the culture resources for biotechnological development on engineering synthetic microbiomes.

Paraburkholderia Symbionts Display Variable Infection Patterns That Are Not Predictive of Amoeba Host Outcomes

Symbiotic interactions exist within a parasitism to mutualism continuum that is influenced, among others, by genes and context. Dynamics of intracellular invasion, replication, and prevalence may underscore both host survivability and symbiont stability. More infectious symbionts might exert higher corresponding costs to hosts, which could ultimately disadvantage both partners. Here, we quantify infection patterns of diverse Paraburkholderia symbiont genotypes in their amoeba host Dictyostelium discoideum and probe the relationship between these patterns and host outcomes. We exposed D. discoideum to thirteen strains of Paraburkholderia each belonging to one of the three symbiont species found to naturally infect D. discoideum: Paraburkholderia agricolaris, Paraburkholderia hayleyella, and Paraburkholderia bonniea. We quantified the infection prevalence and intracellular density of fluorescently labeled symbionts along with the final host population size using flow cytometry and confocal microscopy. We find that infection phenotypes vary across symbiont strains. Symbionts belonging to the same species generally display similar infection patterns but are interestingly distinct when it comes to host outcomes. This results in final infection loads that do not strongly correlate to final host outcomes, suggesting other genetic factors that are not a direct cause or consequence of symbiont abundance impact host fitness.

Transcriptome of Tetranychus urticae embryos reveals insights into Wolbachia-induced cytoplasmic incompatibility

The endosymbiont Wolbachia is known for manipulating host reproduction in selfish ways. However, the molecular mechanisms have not yet been investigated in embryos. Here, we found that Wolbachia had no effect on the number of deposited eggs in Tetranychus urticae Koch (Acari: Tetranychidae) but caused two types of reproductive manipulation: killing uninfected female embryos via cytoplasmic incompatibility (CI) and increasing the hatching ratio of infected female embryos. RNA sequencing analyses showed that 145 genes were differentially expressed between Wolbachia-infected (WI) and Wolbachia-uninfected (WU) embryos. Wolbachia infection down-regulated messenger RNA (mRNA) expression of glutathione S-transferase that could buffer oxidative stress. In addition, 1613 and 294 genes were identified as CI-specific up-/down-regulated genes. Compared to WU and WI embryos, embryos of CI cross strongly expressed genes involved in transcription, translation, tissue morphogenesis, DNA damage and mRNA surveillance. In contrast, most of the genes associated with energy production and metabolism were down-regulated in the CI embryos compared to the WU and WI embryos, which provides some clues as to the cause of death of CI embryos. These results identify several genes that could be candidates for explaining Wolbachia-induced CI. Our data form a basis to help elucidate the molecular consequences of CI in embryos.

Genome Assembly of the A-Group Wolbachia in Nasonia oneida Using Linked-Reads Technology

Wolbachia are obligate intracellular bacteria which commonly infect various nematode and arthropod species. Genome sequences have been generated from arthropod samples following enrichment for the intracellular bacteria, and genomes have also been assembled from arthropod whole-genome sequencing projects. However, these methods remain challenging for infections that occur at low titers in hosts. Here we report the first Wolbachia genome assembled from host sequences using 10× Genomics linked-reads technology. The high read depth attainable by this method allows for recovery of intracellular bacteria that are at low concentrations. Based on the depth differences (714× for the insect and 59× for the bacterium), we assembled the genome of a Wolbachia in the parasitoid jewel wasp species Nasonia oneida. The final draft assembly consists of 1,293, 06 bp in 47 scaffolds with 1,114 coding genes and 97.01% genome completeness assessed by checkM. Comparisons of the five Multi Locus Sequence Typing genes revealed that the sequenced Wolbachia genome is the A1 strain (henceforth wOneA1) previously reported in N. oneida. Pyrosequencing confirms that the wasp strain lacks A2 and B types previously detected in this insect, which were likely lost during laboratory culturing. Assembling bacterial genomes from host genome projects can provide an effective method for sequencing bacterial genomes, even when the infections occur at low density in sampled tissues.

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.

Large Enriched Fragment Targeted Sequencing (LEFT-SEQ) Applied to Capture of Wolbachia Genomes

Symbiosis is a major force of evolutionary change, influencing virtually all aspects of biology, from population ecology and evolution to genomics and molecular/biochemical mechanisms of development and reproduction. A remarkable example is Wolbachia endobacteria, present in some parasitic nematodes and many arthropod species. Acquisition of genomic data from diverse Wolbachia clades will aid in the elucidation of the different symbiotic mechanisms(s). However, challenges of de novo assembly of Wolbachia genomes include the presence in the sample of host DNA: nematode/vertebrate or insect. We designed biotinylated probes to capture large fragments of Wolbachia DNA for sequencing using PacBio technology (LEFT-SEQ: Large Enriched Fragment Targeted Sequencing). LEFT-SEQ was used to capture and sequence four Wolbachia genomes: the filarial nematode Brugia malayi, wBm, (21-fold enrichment), Drosophila mauritiana flies (2 isolates), wMau (11-fold enrichment), and Aedes albopictus mosquitoes, wAlbB (200-fold enrichment). LEFT-SEQ resulted in complete genomes for wBm and for wMau. For wBm, 18 single-nucleotide polymorphisms (SNPs), relative to the wBm reference, were identified and confirmed by PCR. A limit of LEFT-SEQ is illustrated by the wAlbB genome, characterized by a very high level of insertion sequences elements (ISs) and DNA repeats, for which only a 20-contig draft assembly was achieved.

Symbiosis Comes of age at the 10th Biennial Meeting of Wolbachia Researchers

Wolbachia pipientis is an alpha-proteobacterial, obligate intracellular microbe and arguably the most successful infection on our planet, colonizing 40-60% of insect species. Wolbachia are also present in most, but not all, filarial nematodes where they are obligate mutualists and are the targets for anti-filarial drug discovery. Although Wolbachia are related to important human pathogens they do not infect mammals, but instead are well known for their reproductive manipulations of insect populations, inducing the following phenotypes: male-killing, feminization, parthenogenesis induction, or cytoplasmic incompatibility (CI). The most common of these, CI, results in a sperm-egg incompatibility and increases the relative fecundity of infected females in a population. In the last decade, Wolbachia have also been shown to provide a benefit to insects, where the infection can inhibit RNA virus replication within the host. Wolbachia cannot be cultivated outside of host cells and no genetic tools are available in the symbiont, limiting approaches available to its study. This means that many questions fundamental to our understanding of Wolbachia basic biology remained unknown for decades. The tenth biennial international Wolbachia conference, "Wolbachia Evolution, Ecology, Genomics and Cell Biology: A Chronicle of the Most Ubiquitous Symbiont", was held on June 17-22, 2018, Salem, MA USA. In the review below we highlight the new science presented at the meeting, link it to prior efforts to answer these questions across the Wolbachia genus, and the importance to the field of symbiosis. The topics covered in this review are based on the presentations at the conference.

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.

Phylogenetic analysis based on the 16S rDNA, gltA, gatB, and hcpA gene sequences of Wolbachia from the novel host Ceratozetes thienemanni (Acari: Oribatida)

We determined the occurrence of intracellular endosymbionts (Wolbachia, Cardinium, Arsenophonus, Rickettsia, Spiroplasma, Hamiltonella, flavobacteria, and microsporidia) in oribatid mites (Acari: Oribatida) with the use of PCR technique. For the first time we looked for and detected Wolbachia in parthenogenetic oribatid mite Ceratozetes thienemanni Willmann, 1943. The 16S rDNA, gatB, hcpA, and gltA sequences of Wolbachia in C. thienemanni showed the highest similarity (≥ 90%) to the genes of Wolbachia from springtails (Collembola) and oribatid mite Gustavia microcephala. We found the unique sequence 5'-GGGGTAATGGCC-3' in 16S rDNA of Wolbachia from C. thienemanni and collembolan representing group E. The phylogeny of Wolbachia based on the analysis of single genes as well as concatenated alignments of four bacterial loci showed that the bacteria from C. thienemanni belonged to Wolbachia group E, like the endosymbionts from springtail hosts and G. microcephala. Considering coexisting of representatives of Oribatida and Collembola in the same soil habitat and similar food, it is possible that the source of Wolbachia infection was the same. Residues of dead invertebrates could be in organic matter of their soil food, so the scenario of infection transferred by eating of remains of soil cohabitates is also possible. It could explain the similarity and relationship of the Wolbachia in these two arthropod groups. Oribatid mite C. thienemanni is a parthenogenetic mite which is a unique feature in the genus Ceratozetes. Moreover, this species, within the entire genus Ceratozetes, is characterized by the most northerly distribution. It is difficult to determine either it is parthenogenesis or the presence of endosymbionts that are in some way responsible for this kind of evolutionary success. Maybe we are dealing here with a kind of synergy of both factors?

Evolutionary Genetics of Cytoplasmic Incompatibility Genes cifA and cifB in Prophage WO of Wolbachia

The bacterial endosymbiont Wolbachia manipulates arthropod reproduction to facilitate its maternal spread through host populations. The most common manipulation is cytoplasmic incompatibility (CI): Wolbachia-infected males produce modified sperm that cause embryonic mortality, unless rescued by embryos harboring the same Wolbachia. The genes underlying CI, cifA and cifB, were recently identified in the eukaryotic association module of Wolbachia’s prophage WO. Here, we use transcriptomic and genomic approaches to address three important evolutionary facets of the cif genes. First, we assess whether or not cifA and cifB comprise a classic toxin–antitoxin operon in wMel and show that the two genes exhibit striking, transcriptional differences across host development. They can produce a bicistronic message despite a predicted hairpin termination element in their intergenic region. Second, cifA and cifB strongly coevolve across the diversity of phage WO. Third, we provide new domain and functional predictions across homologs within Wolbachia, and show that amino acid sequences vary substantially across the genus. Finally, we investigate conservation of cifA and cifB and find frequent degradation and loss of the genes in strains that no longer induce CI. Taken together, we demonstrate that cifA and cifB exhibit complex transcriptional regulation in wMel, provide functional annotations that broaden the potential mechanisms of CI induction, and report recurrent erosion of cifA and cifB in non-CI strains, thus expanding our understanding of the most widespread form of reproductive parasitism.

Uncovering the hidden players in Lepidoptera biology: the heritable microbial endosymbionts

The Lepidoptera is one of the most widespread and recognisable insect orders. Due to their remarkable diversity, economic and ecological importance, moths and butterflies have been studied extensively over the last 200 years. More recently, the relationship between Lepidoptera and their heritable microbial endosymbionts has received increasing attention. Heritable endosymbionts reside within the host’s body and are often, but not exclusively, inherited through the female line. Advancements in molecular genetics have revealed that host-associated microbes are both extremely prevalent among arthropods and highly diverse. Furthermore, heritable endosymbionts have been repeatedly demonstrated to play an integral role in many aspects of host biology, particularly host reproduction. Here, we review the major findings of research of heritable microbial endosymbionts of butterflies and moths. We promote the Lepidoptera as important models in the study of reproductive manipulations employed by heritable endosymbionts, with the mechanisms underlying male-killing and feminisation currently being elucidated in moths and butterflies. We also reveal that the vast majority of research undertaken of Lepidopteran endosymbionts concerns Wolbachia. While this highly prevalent bacterium is undoubtedly important, studies should move towards investigating the presence of other, and interacting endosymbionts, and we discuss the merits of examining the microbiome of Lepidoptera to this end. We finally consider the importance of understanding the influence of endosymbionts under global environmental change and when planning conservation management of endangered Lepidoptera species.

Transcriptome Sequencing Reveals Novel Candidate Genes for Cardinium hertigii-Caused Cytoplasmic Incompatibility and Host-Cell Interaction

The majority of insects carry maternally inherited intracellular bacteria that are important in their hosts’ biology, ecology, and evolution. Some of these bacterial symbionts cause a reproductive failure known as cytoplasmic incompatibility (CI). In CI, the mating of symbiont-infected males and uninfected females produces few or no daughters. The CI symbiont then spreads and can have a significant impact on the insect host population. Cardinium, a bacterial endosymbiont of the parasitoid wasp Encarsia in the Bacteroidetes, is the only bacterial lineage known to cause CI outside the Alphaproteobacteria, where Wolbachia and another recently discovered CI symbiont reside. Here, we sought insight into the gene expression of a CI-inducing Cardinium strain in its natural host, Encarsia suzannae. Our study provides the first insights into the Cardinium transcriptome and provides support for the hypothesis that Wolbachia and Cardinium target similar host pathways with distinct and largely unrelated sets of genes.

Cytoplasmic incompatibility (CI) is an intriguing, widespread, symbiont-induced reproductive failure that decreases offspring production of arthropods through crossing incompatibility of infected males with uninfected females or with females infected with a distinct symbiont genotype. For years, the molecular mechanism of CI remained unknown. Recent genomic, proteomic, biochemical, and cell biological studies have contributed to understanding of CI in the alphaproteobacterium Wolbachia and implicate genes associated with the WO prophage. Besides a recently discovered additional lineage of alphaproteobacterial symbionts only moderately related to Wolbachia, Cardinium (Bacteroidetes) is the only other symbiont known to cause CI, and genomic evidence suggests that it has very little homology with Wolbachia and evolved this phenotype independently. Here, we present the first transcriptomic study of the CI Cardinium strain cEper1, in its natural host, Encarsia suzannae, to detect important CI candidates and genes involved in the insect-Cardinium symbiosis. Highly expressed transcripts included genes involved in manipulating ubiquitination, apoptosis, and host DNA. Female-biased genes encoding ribosomal proteins suggest an increase in general translational activity of Cardinium in female wasps. The results confirm previous genomic analyses that indicated that Wolbachia and Cardinium utilize different genes to induce CI, and transcriptome patterns further highlight expression of some common pathways that these bacteria use to interact with the host and potentially cause this enigmatic and fundamental manipulation of host reproduction.

IMPORTANCE The majority of insects carry maternally inherited intracellular bacteria that are important in their hosts’ biology, ecology, and evolution. Some of these bacterial symbionts cause a reproductive failure known as cytoplasmic incompatibility (CI). In CI, the mating of symbiont-infected males and uninfected females produces few or no daughters. The CI symbiont then spreads and can have a significant impact on the insect host population. Cardinium, a bacterial endosymbiont of the parasitoid wasp Encarsia in the Bacteroidetes, is the only bacterial lineage known to cause CI outside the Alphaproteobacteria, where Wolbachia and another recently discovered CI symbiont reside. Here, we sought insight into the gene expression of a CI-inducing Cardinium strain in its natural host, Encarsia suzannae. Our study provides the first insights into the Cardinium transcriptome and provides support for the hypothesis that Wolbachia and Cardinium target similar host pathways with distinct and largely unrelated sets of genes.

Comparative genomics provides a timeframe for Wolbachia evolution and exposes a recent biotin synthesis operon transfer

The genus Wolbachia (Alphaproteobacteria) comprises the most abundant inherited intracellular bacteria¹. Despite their relevance as manipulators of human pathogen transmission² and arthropod reproduction³, many aspects of their evolutionary history are not well understood⁴. In arthropods, Wolbachia infections are typically transient on evolutionary timescales^5,6 and co-divergence between hosts and Wolbachia is supposedly rare. Consequently, much of our knowledge of Wolbachia genome evolution derives from very recently diverged strains, and a timescale for Wolbachia is lacking. Here, we investigated the genomes of four Wolbachia strains that have persisted within and co-diverged with their host lineage for ∼2 million years. Although the genomes showed very little evolutionary change on a nucleotide level, we found evidence for a recent lateral transfer of a complete biotin synthesis operon that has the potential to transform Wolbachia-host relationships⁷. Furthermore, this evolutionary snapshot enabled us to calibrate the divergence times of the supergroup A and B Wolbachia lineages using genome-wide data sets and relaxed molecular clock models. We estimated the origin of Wolbachia supergroups A and B to be ∼200 million years ago (Ma), which is considerably older than previously appreciated. This age coincides with the diversification of many insect lineages⁸ that represent most of Wolbachia's host spectrum.