Evolution of symbiotic organs and endosymbionts in lygaeid stinkbugs

We investigated seed bugs of the genus Nysius (Insecta: Hemiptera: Lygaeidae) for their symbiotic bacteria. From all the samples representing 4 species, 18 populations and 281 individuals, specific bacterial 16S rRNA gene sequences were consistently identified, which formed a distinct clade in the Gammaproteobacteria. In situ hybridization showed that the bacterium was endocellularly localized in a pair of large bacteriomes that were amorphous in shape, deep red in color, and in association with gonads. In the ovary of adult females, the endosymbiont was also localized in the ‘infection zone'' in the middle of each germarium and in the ‘symbiont ball'' at the anterior pole of each oocyte, indicating vertical transmission of the endosymbiont through the ovarial passage. Phylogenetic analyses based on bacterial 16S rRNA, groEL and gyrB genes consistently supported a coherent monophyly of the Nysius endosymbionts. The possibility of a sister relationship to ‘Candidatus Kleidoceria schneideri'', the bacteriome-associated endosymbiont of a lygaeid bug Kleidocerys resedae, was statistically rejected, indicating independent evolutionary origins of the endosymbionts in the Lygaeidae. The endosymbiont genes consistently exhibited AT-biased nucleotide compositions and accelerated rates of molecular evolution, and the endosymbiont genome was only 0.6 Mb in size. The endosymbiont phylogeny was congruent with the host insect phylogeny, suggesting strict vertical transmission and host–symbiont co-speciation over evolutionary time. Based on these results, we discuss the evolution of bacteriomes and endosymbionts in the Heteroptera, most members of which are associated with gut symbiotic bacteria. The designation ‘Candidatus Schneideria nysicola'' is proposed for the endosymbiont clade.     

Diverse Phage-Encoded Toxins in a Protective Insect Endosymbiont

The lysogenic bacteriophage APSE infects “Candidatus Hamiltonella defensa,” a facultative endosymbiont of aphids and other sap-feeding insects. This endosymbiont has established a beneficial association with aphids, increasing survivorship following attack by parasitoid wasps. Although APSE and “Ca. Hamiltonella defensa” are effectively maternally transmitted between aphid generations, they can also be horizontally transferred among insect hosts, which results in genetically distinct “Ca. Hamiltonella defensa” strains infecting the same aphid species and sporadic distributions of both APSE and “Ca. Hamiltonella defensa” among hosts. Aphids infected only with “Ca. Hamiltonella defensa” have significantly less protection than those infected with both “Ca. Hamiltonella defensa” and APSE. This protection has been proposed to be connected to eukaryote-targeted toxins previously discovered in the genomes of two characterized APSE strains. In this study, we have sequenced partial genomes from seven additional APSE strains to address the evolution and extent of toxin variation in this phage. The APSE lysis region has been a hot spot for nonhomologous recombination of novel virulence cassettes. We identified four new toxins from three protein families, Shiga-like toxin, cytolethal distending toxin, and YD-repeat toxins. These recombination events have also resulted in reassortment of the downstream lysozyme and holin genes. Analysis of the conserved APSE genes flanking the variable toxin cassettes reveals a close phylogenetic association with phage sequences from two other facultative endosymbionts of insects. Thus, phage may act as a conduit for ongoing gene exchange among heritable endosymbionts.     

Isolation, Pure Culture, and Characterization of “Candidatus Arsenophonus arthropodicus,” an Intracellular Secondary Endosymbiont from the Hippoboscid Louse Fly Pseudolynchia canariensis

Members of the genus Arsenophonus comprise a large group of bacterial endosymbionts that are widely distributed in arthropods of medical, veterinary, and agricultural importance. At present, little is known about the role of these bacteria in arthropods, because few representatives have been isolated and cultured in the laboratory. In the current study, we describe the isolation and pure culture of an Arsenophonus endosymbiont from the hippoboscid louse fly Pseudolynchia canariensis. We propose provisional nomenclature for this bacterium in the genus Arsenophonus as “Candidatus Arsenophonus arthropodicus.” Phylogenetic analyses indicate that “Candidatus Arsenophonus arthropodicus” is closely related to the Arsenophonus endosymbionts found in psyllids, whiteflies, aphids, and mealybugs. The pure culture of this endosymbiont offers new opportunities to examine the role of Arsenophonus in insects. To this end, we describe methods for the culture of “Candidatus Arsenophonus arthropodicus” in an insect cell line and the transformation of this bacterium with a broad-host-range plasmid.     

Complete Genome Sequence of “Candidatus Portiera aleyrodidarum” BT-QVLC,an Obligate Symbiont That Supplies Amino Acids and Carotenoids to Bemisia tabaci

The genome of “Candidatus Portiera aleyrodidarum,” the primary endosymbiont of the whitefly Bemisia tabaci (Mediterranean species), is reported. It presents a reduced genome (357 kb) encoding the capability to synthetize, or participate in the synthesis of, several amino acids and carotenoids, being the first insect endosymbiont capable of supplying carotenoids.     

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis

Background

Insertion sequences (ISs) are approximately 1 kbp long “jumping” genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (“cut-and-paste”) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium ‘Nostoc azollae’ 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of ‘Nostoc azollae’ 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria.

Results

Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (“local hopping”) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of ‘Nostoc azollae’ 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style.

Conclusions

The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium ‘Nostoc azollae’ 0708 and in other obligate prokaryotic symbionts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1386-7) contains supplementary material, which is available to authorized users.     

A Novel Human-Infection-Derived Bacterium Provides Insights into the Evolutionary Origins of Mutualistic Insect–Bacterial Symbioses

Despite extensive study, little is known about the origins of the mutualistic bacterial endosymbionts that inhabit approximately 10% of the world''s insects. In this study, we characterized a novel opportunistic human pathogen, designated “strain HS,” and found that it is a close relative of the insect endosymbiont Sodalis glossinidius. Our results indicate that ancestral relatives of strain HS have served as progenitors for the independent descent of Sodalis-allied endosymbionts found in several insect hosts. Comparative analyses indicate that the gene inventories of the insect endosymbionts were independently derived from a common ancestral template through a combination of irreversible degenerative changes. Our results provide compelling support for the notion that mutualists evolve from pathogenic progenitors. They also elucidate the role of degenerative evolutionary processes in shaping the gene inventories of symbiotic bacteria at a very early stage in these mutualistic associations.     

Arsenophonus and Sodalis Symbionts in Louse Flies: an Analogy to the Wigglesworthia and Sodalis System in Tsetse Flies

Symbiosis between insects and bacteria result in a variety of arrangements, genomic modifications, and metabolic interconnections. Here, we present genomic, phylogenetic, and morphological characteristics of a symbiotic system associated with Melophagus ovinus, a member of the blood-feeding family Hippoboscidae. The system comprises four unrelated bacteria representing different stages in symbiosis evolution, from typical obligate mutualists inhabiting bacteriomes to freely associated commensals and parasites. Interestingly, the whole system provides a remarkable analogy to the association between Glossina and its symbiotic bacteria. In both, the symbiotic systems are composed of an obligate symbiont and two facultative intracellular associates, Sodalis and Wolbachia. In addition, extracellular Bartonella resides in the gut of Melophagus. However, the phylogenetic origins of the two obligate mutualist symbionts differ. In Glossina, the mutualistic Wigglesworthia appears to be a relatively isolated symbiotic lineage, whereas in Melophagus, the obligate symbiont originated within the widely distributed Arsenophonus cluster. Although phylogenetically distant, the two obligate symbionts display several remarkably similar traits (e.g., transmission via the host''s “milk glands” or similar pattern of genome reduction). To obtain better insight into the biology and possible role of the M. ovinus obligate symbiont, “Candidatus Arsenophonus melophagi,” we performed several comparisons of its gene content based on assignments of the Cluster of Orthologous Genes (COG). Using this criterion, we show that within a set of 44 primary and secondary symbionts, “Ca. Arsenophonus melophagi” is most similar to Wigglesworthia. On the other hand, these two bacteria also display interesting differences, such as absence of flagellar genes in Arsenophonus and their presence in Wigglesworthia. This finding implies that a flagellum is not essential for bacterial transmission via milk glands.     

Evolutionary Relationships of Three New Species of Enterobacteriaceae Living as Symbionts of Aphids and Other Insects

Ecological studies on three bacterial lineages symbiotic in aphids have shown that they impose a variety of effects on their hosts, including resistance to parasitoids and tolerance to heat stress. Phylogenetic analyses of partial sequences of gyrB and recA are consistent with previous analyses limited to 16S rRNA gene sequences and yield improved confidence of the evolutionary relationships of these symbionts. All three symbionts are in the Enterobacteriaceae. One of the symbionts, here given the provisional designation “Candidatus Serratia symbiotica,” is a Serratia species that has acquired a symbiotic lifestyle. The other two symbionts, here designated “Candidatus Hamiltonella defensa” and “Candidatus Regiella insecticola,” are sister groups to one another and together show a relationship to species of Photorhabdus.     

Evolutionary Relationships of “Candidatus Riesia spp.,” Endosymbiotic Enterobacteriaceae Living within Hematophagous Primate Lice

The primary endosymbiotic bacteria from three species of parasitic primate lice were characterized molecularly. We have confirmed the characterization of the primary endosymbiont (P-endosymbiont) of the human head/body louse Pediculus humanus and provide new characterizations of the P-endosymbionts from Pediculus schaeffi from chimpanzees and Pthirus pubis, the pubic louse of humans. The endosymbionts show an average percent sequence divergence of 11 to 15% from the most closely related known bacterium “Candidatus Arsenophonus insecticola.” We propose that two additional species be added to the genus “Candidatus Riesia.” The new species proposed within “Candidatus Riesia” have sequence divergences of 3.4% and 10 to 12% based on uncorrected pairwise differences. Our Bayesian analysis shows that the branching pattern for the primary endosymbionts was the same as that for their louse hosts, suggesting a long coevolutionary history between primate lice and their primary endosymbionts. We used a calibration of 5.6 million years to date the divergence between endosymbionts from human and chimpanzee lice and estimated an evolutionary rate of nucleotide substitution of 0.67% per million years, which is 15 to 30 times faster than previous estimates calculated for Buchnera, the primary endosymbiont in aphids. Given the evidence for cospeciation with primate lice and the evidence for fast evolutionary rates, this lineage of endosymbiotic bacteria can be evaluated as a fast-evolving marker of both louse and primate evolutionary histories.     

Comparison of the Genome Sequences of “Candidatus Portiera aleyrodidarum” Primary Endosymbionts of the Whitefly Bemisia tabaci B and Q Biotypes

“Candidatus Portiera aleyrodidarum” is the primary endosymbiont of whiteflies. We report two complete genome sequences of this bacterium from the worldwide invasive B and Q biotypes of the whitefly Bemisia tabaci. Differences in the two genome sequences may add insights into the complex differences in the biology of both biotypes.     

Assessment of Survival and Body Size Variation of Culicoides imicola (Diptera: Ceratopogonidae) as Functions of “Candidatus Cardinium” (Bacteroidetes) Infection Status

“Candidatus Cardinium hertigii” (Bacteroidetes) is a maternally inherited endosymbiont known from several arthropods. Its mechanisms for persistence in host populations are mostly reproductive manipulation, though it has been occasionally reported to improve fitness parameters in several hosts. In Culicoides (Diptera: Ceratopogonidae) biting midges, the prevalence of “Candidatus Cardinium” infection was documented as moderate, with no detectable sex bias. We therefore investigated whether “Candidatus Cardinium” affects important fitness parameters, such as survival and body size, in Culicoides imicola, a dominant vector species. Field-collected midges were trapped and analyzed for survival under different environmental conditions and antibiotic treatment, taking into account “Candidatus Cardinium” infection status and parity status (i.e., parous or nulliparous). Additionally, wing lengths were measured as a proxy parameter for body size and analyzed together with “Candidatus Cardinium” infection data. The findings revealed no difference in survival of Culicoides infected with “Candidatus Cardinium” and that of uninfected midges in both parity states and under all tested conditions: optimal, starvation, heat, and antibiotic treatment. Beyond survival, no wing length difference was found for “Candidatus Cardinium”-infected versus uninfected midges. In aggregate, these findings support our conclusion that “Candidatus Cardinium” does not have an overt effect on the survival and size of adult C. imicola midges. “Candidatus Cardinium” may affect immature stages or may alter adult reproductive performance.     

Endosymbiotic Bacteroidales Bacteria of the Flagellated Protist Pseudotrichonympha grassii in the Gut of the Termite Coptotermes formosanus

A unique lineage of bacteria belonging to the order Bacteroidales was identified as an intracellular endosymbiont of the protist Pseudotrichonympha grassii (Parabasalia, Hypermastigea) in the gut of the termite Coptotermes formosanus. We identified the 16S rRNA, gyrB, elongation factor Tu, and groEL gene sequences in the endosymbiont and detected a very low level of sequence divergence (<0.9% of the nucleotides) in the endosymbiont population within and among protist cells. The Bacteroidales endosymbiont sequence was affiliated with a cluster comprising only sequences from termite gut bacteria and was not closely related to sequences identified for members of the Bacteroidales attached to the cell surfaces of other gut protists. Transmission electron microscopy showed that there were numerous rod-shaped bacteria in the cytoplasm of the host protist, and we detected the endosymbiont by fluorescence in situ hybridization (FISH) with an oligonucleotide probe specific for the 16S rRNA gene identified. Quantification of the abundance of the Bacteroidales endosymbiont by sequence-specific cleavage of rRNA with RNase H and FISH cell counting revealed, surprisingly, that the endosymbiont accounted for 82% of the total bacterial rRNA and 71% of the total bacterial cells in the gut community. The genetically nearly homogeneous endosymbionts of Pseudotrichonympha were very abundant in the gut symbiotic community of the termite.     

Phylogeny of Symbiotic Genes and the Symbiotic Properties of Rhizobia Specific to Astragalus glycyphyllos L.

The phylogeny of symbiotic genes of Astragalus glycyphyllos L. (liquorice milkvetch) nodule isolates was studied by comparative sequence analysis of nodA, nodC, nodH and nifH loci. In all these genes phylograms, liquorice milkvetch rhizobia (closely related to bacteria of three species, i.e. Mesorhizobium amorphae, Mesorhizobium septentrionale and Mesorhizobium ciceri) formed one clearly separate cluster suggesting the horizontal transfer of symbiotic genes from a single ancestor to the bacteria being studied. The high sequence similarity of the symbiotic genes of A. glycyphyllos rhizobia (99–100% in the case of nodAC and nifH genes, and 98–99% in the case of nodH one) points to the relatively recent (in evolutionary scale) lateral transfer of these genes. In the nodACH and nifH phylograms, A. glycyphyllos nodule isolates were grouped together with the genus Mesorhizobium species in one monophyletic clade, close to M. ciceri, Mesorhizobium opportunistum and Mesorhizobium australicum symbiovar biserrulae bacteria, which correlates with the close relationship of these rhizobia host plants. Plant tests revealed the narrow host range of A. glycyphyllos rhizobia. They formed effective symbiotic interactions with their native host (A. glycyphyllos) and Amorpha fruticosa but not with 11 other fabacean species. The nodules induced on A. glycyphyllos roots were indeterminate with apical, persistent meristem, an age gradient of nodule tissues and cortical vascular bundles. To reflect the symbiosis-adaptive phenotype of rhizobia, specific for A. glycyphyllos, we propose for these bacteria the new symbiovar “glycyphyllae”, based on nodA and nodC genes sequences.     

Enhanced Mutation Rate,Relaxed Selection,and the “Domino Effect” are associated with Gene Loss in Blattabacterium,A Cockroach Endosymbiont

Intracellular endosymbionts have reduced genomes that progressively lose genes at a timescale of tens of million years. We previously reported that gene loss rate is linked to mutation rate in Blattabacterium, however, the mechanisms causing gene loss are not yet fully understood. Here, we carried out comparative genomic analyses on the complete genome sequences of a representative set of 67 Blattabacterium strains, with sizes ranging between 511 and 645 kb. We found that 200 of the 566 analyzed protein-coding genes were lost in at least one lineage of Blattabacterium, with the most extreme case being one gene that was lost independently in 24 lineages. We found evidence for three mechanisms influencing gene loss in Blattabacterium. First, gene loss rates were found to increase exponentially with the accumulation of substitutions. Second, genes involved in vitamin and amino acid metabolism experienced relaxed selection in Cryptocercus and Mastotermes, possibly triggered by their vertically inherited gut symbionts. Third, we found evidence of epistatic interactions among genes leading to a “domino effect” of gene loss within pathways. Our results highlight the complexity of the process of genome erosion in an endosymbiont.     

The Pine Bark Adelgid,Pineus strobi,Contains Two Novel Bacteriocyte-Associated Gammaproteobacterial Symbionts

Bacterial endosymbionts of the pine bark adelgid, Pineus strobi (Insecta: Hemiptera: Adelgidae), were investigated using transmission electron microscopy, 16S and 23S rRNA-based phylogeny, and fluorescence in situ hybridization. Two morphologically different symbionts affiliated with the Gammaproteobacteria were present in distinct bacteriocytes. One of them (“Candidatus Annandia pinicola”) is most closely related to an endosymbiont of Adelges tsugae, suggesting that they originate from a lineage already present in ancient adelgids before the hosts diversified into the two major clades, Adelges and Pineus. The other P. strobi symbiont (“Candidatus Hartigia pinicola”) represents a novel symbiont lineage in members of the Adelgidae. Our findings lend further support for a complex evolutionary history of the association of adelgids with a phylogenetically diverse set of bacterial symbionts.     

A Novel Bacteroidetes Symbiont Is Localized in Scaphoideus titanus, the Insect Vector of Flavescence Dorée in Vitis vinifera

Flavescence dorée (FD) is a grapevine disease that afflicts several wine production areas in Europe, from Portugal to Serbia. FD is caused by a bacterium, “Candidatus Phytoplasma vitis,” which is spread throughout the vineyards by a leafhopper, Scaphoideus titanus (Cicadellidae). After collection of S. titanus specimens from FD-contaminated vineyards in three different areas in the Piedmont region of Italy, we performed a survey to characterize the bacterial microflora associated with this insect. Using length heterogeneity PCR with universal primers for bacteria we identified a major peak associated with almost all of the individuals examined (both males and females). Characterization by denaturing gradient gel electrophoresis confirmed the presence of a major band that, after sequencing, showed a 97 to 99% identity with Bacteroidetes symbionts of the “Candidatus Cardinium hertigii” group. In addition, electron microscopy of tissues of S. titanus fed for 3 months on phytoplasma-infected grapevine plants showed bacterial cells with the typical morphology of “Ca. Cardinium hertigii.” This endosymbiont, tentatively designated ST1-C, was found in the cytoplasm of previtellogenic and vitellogenic ovarian cells, in the follicle cells, and in the fat body and salivary glands. In addition, cell morphologies resembling those of “Ca. Phytoplasma vitis” were detected in the midgut, and specific PCR assays indicated the presence of the phytoplasma in the gut, fat body and salivary glands. These results indicate that ST1-C and “Ca. Phytoplasma vitis” have a complex life cycle in the body of S. titanus and are colocalized in different organs and tissues.     

Pleiotropic Impact of Endosymbiont Load and Co-Occurrence in the Maize Weevil Sitophilus zeamais

Individual traits vary among and within populations, and the co-occurrence of different endosymbiont species within a host may take place under varying endosymbiont loads in each individual host. This makes the recognition of the potential impact of such endosymbiont associations in insect species difficult, particularly in insect pest species. The maize weevil, Sitophilus zeamais Motsch. (Coleoptera: Curculionidae), a key pest species of stored cereal grains, exhibits associations with two endosymbiotic bacteria: the obligatory endosymbiont SZPE (“Sitophilus zeamais Primary Endosymbiont”) and the facultative endosymbiont Wolbachia. The impact of the lack of SZPE in maize weevil physiology is the impairment of nutrient acquisition and energy metabolism, while Wolbachia is an important factor in reproductive incompatibility. However, the role of endosymbiont load and co-occurrence in insect behavior, grain consumption, body mass and subsequent reproductive factors has not yet been explored. Here we report on the impacts of co-occurrence and varying endosymbiont loads achieved via thermal treatment and antibiotic provision via ingested water in the maize weevil. SZPE exhibited strong effects on respiration rate, grain consumption and weevil body mass, with observed effects on weevil behavior, particularly flight activity, and potential consequences for the management of this pest species. Wolbachia directly favored weevil fertility and exhibited only mild indirect effects, usually enhancing the SZPE effect. SZPE suppression delayed weevil emergence, which reduced the insect population growth rate, and the thermal inactivation of both symbionts prevented insect reproduction. Such findings are likely important for strain divergences reported in the maize weevil and their control, aspects still deserving future attention.     

“Candidatus Accumulibacter” Population Structure in Enhanced Biological Phosphorus Removal Sludges as Revealed by Polyphosphate Kinase Genes

We investigated the fine-scale population structure of the “Candidatus Accumulibacter” lineage in enhanced biological phosphorus removal (EBPR) systems using the polyphosphate kinase 1 gene (ppk1) as a genetic marker. We retrieved fragments of “Candidatus Accumulibacter” 16S rRNA and ppk1 genes from one laboratory-scale and several full-scale EBPR systems. Phylogenies reconstructed using 16S rRNA genes and ppk1 were largely congruent, with ppk1 granting higher phylogenetic resolution and clearer tree topology and thus serving as a better genetic marker than 16S rRNA for revealing population structure within the “Candidatus Accumulibacter” lineage. Sequences from at least five clades of “Candidatus Accumulibacter” were recovered by ppk1-targeted PCR, and subsequently, specific primer sets were designed to target the ppk1 gene for each clade. Quantitative real-time PCR (qPCR) assays using “Candidatus Accumulibacter”-specific 16S rRNA and “Candidatus Accumulibacter” clade-specific ppk1 primers were developed and conducted on three laboratory-scale and nine full-scale EBPR samples and two full-scale non-EBPR samples to determine the abundance of the total “Candidatus Accumulibacter” lineage and the relative distributions and abundances of the five “Candidatus Accumulibacter” clades. The qPCR-based estimation of the total “Candidatus Accumulibacter” fraction as a proportion of the bacterial community as measured using 16S rRNA genes was not significantly different from the estimation measured using ppk1, demonstrating the power of ppk1 as a genetic marker for detection of all currently defined “Candidatus Accumulibacter” clades. The relative distributions of “Candidatus Accumulibacter” clades varied among different EBPR systems and also temporally within a system. Our results suggest that the “Candidatus Accumulibacter” lineage is more diverse than previously realized and that different clades within the lineage are ecologically distinct.