The process of genome shrinkage in the obligate symbiont Buchnera aphidicola

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.

Coordination of host and endosymbiont gene expression governs endosymbiont growth and elimination in the cereal weevil Sitophilus spp

BACKGROUND

Insects living in nutritionally poor environments often establish long-term relationships with intracellular bacteria that supplement their diets and improve their adaptive and invasive powers. Even though these symbiotic associations have been extensively studied on physiological, ecological, and evolutionary levels, few studies have focused on the molecular dialogue between host and endosymbionts to identify genes and pathways involved in endosymbiosis control and dynamics throughout host development.

RESULTS

We simultaneously analyzed host and endosymbiont gene expression during the life cycle of the cereal weevil Sitophilus oryzae, from larval stages to adults, with a particular emphasis on emerging adults where the endosymbiont Sodalis pierantonius experiences a contrasted growth-climax-elimination dynamics. We unraveled a constant arms race in which different biological functions are intertwined and coregulated across both partners. These include immunity, metabolism, metal control, apoptosis, and bacterial stress response.

CONCLUSIONS

The study of these tightly regulated functions, which are at the center of symbiotic regulations, provides evidence on how hosts and bacteria finely tune their gene expression and respond to different physiological challenges constrained by insect development in a nutritionally limited ecological niche. Video Abstract.

Frenemies on the reef? Resolving the coral-Endozoicomonas association

Stony corals are poster child holobionts due to their intimate association with diverse microorganisms from all domains of life. We are only beginning to understand the diverse functions of most of these microbial associates, including potential main contributors to holobiont health and resilience. Among these, bacteria of the elusive genus Endozoicomonas are widely perceived as beneficial symbionts based on their genomic potential and their high prevalence and ubiquitous presence in coral tissues. Simultaneously, evidence of pathogenic and parasitic Endozoicomonas lineages in other marine animals is emerging. Synthesizing the current knowledge on the association of Endozoicomonas with marine holobionts, we challenge the perception of a purely mutualistic coral-Endozoicomonas relationship and propose directions to elucidate its role along the symbiotic spectrum.

Selection for oligotrophy among bacteria inhabiting host microbiomes

IMPORTANCE

Understanding how natural selection has historically shaped the traits of microbial populations comprising host microbiomes would help predict how the functions of these microbes may continue to evolve over space and time. Numerous host-associated microbes have been found to adapt to their host, sometimes becoming obligate symbionts, whereas free-living microbes are best known to adapt to their surrounding environment. Our study assessed the selective pressures of both the host environment and the surrounding external environment in shaping the functional potential of host-associated bacteria. Despite residing within the resource-rich microbiome of their hosts, we demonstrate that host-associated heterotrophic bacteria show evidence of trait selection that matches the nutrient availability of their broader surrounding environment. These findings illustrate the complex mix of selective pressures that likely shape the present-day function of bacteria found inhabiting host microbiomes. Our study lends insight into the shifts in function that may occur as environments fluctuate over time.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

Linking prokaryotic genome size variation to metabolic potential and environment

While theories and models have appeared to explain genome size as a result of evolutionary processes, little work has shown that genome sizes carry ecological signatures. Our work delves into the ecological implications of microbial genome size variation in benthic and pelagic habitats across environmental gradients of the brackish Baltic Sea. While depth is significantly associated with genome size in benthic and pelagic brackish metagenomes, salinity is only correlated to genome size in benthic metagenomes. Overall, we confirm that prokaryotic genome sizes in Baltic sediments (3.47 Mbp) are significantly bigger than in the water column (2.96 Mbp). While benthic genomes have a higher number of functions than pelagic genomes, the smallest genomes coded for a higher number of module steps per Mbp for most of the functions irrespective of their environment. Some examples of this functions are amino acid metabolism and central carbohydrate metabolism. However, we observed that nitrogen metabolism was almost absent in pelagic genomes and was mostly present in benthic genomes. Finally, we also show that Bacteria inhabiting Baltic sediments and water column not only differ in taxonomy, but also in their metabolic potential, such as the Wood-Ljungdahl pathway or the presence of different hydrogenases. Our work shows how microbial genome size is linked to abiotic factors in the environment, metabolic potential and taxonomic identity of Bacteria and Archaea within aquatic ecosystems.

Weevil Carbohydrate Intake Triggers Endosymbiont Proliferation: A Trade-Off between Host Benefit and Endosymbiont Burden

Nutritional symbioses between insects and intracellular bacteria (endosymbionts) are a major force of adaptation, allowing animals to colonize nutrient-poor ecological niches. Many beetles feeding on tyrosine-poor substrates rely on a surplus of aromatic amino acids produced by bacterial endosymbionts. This surplus of aromatic amino acids is crucial for the biosynthesis of a thick exoskeleton, the cuticle, which is made of a matrix of chitin with proteins and pigments built from tyrosine-derived molecules, providing an important defensive barrier against biotic and abiotic stress. Other endosymbiont-related advantages for beetles include faster development and improved fecundity. The association between Sitophilus oryzae and the Sodalis pierantonius endosymbiont represents a unique case study among beetles: endosymbionts undergo an exponential proliferation in young adults concomitant with the cuticle tanning, and then they are fully eliminated. While endosymbiont clearance, as well as total endosymbiont titer, are host-controlled processes, the mechanism triggering endosymbiont exponential proliferation remains poorly understood. Here, we show that endosymbiont exponential proliferation relies on host carbohydrate intake, unlike the total endosymbiont titer or the endosymbiont clearance, which are under host genetic control. Remarkably, insect fecundity was preserved, and the cuticle tanning was achieved, even when endosymbiont exponential proliferation was experimentally blocked, except in the context of a severely unbalanced diet. Moreover, a high endosymbiont titer coupled with nutrient shortage dramatically impacted host survival, revealing possible environment-dependent disadvantages for the host, likely due to the high energy cost of exponentially proliferating endosymbionts. IMPORTANCE Beetles thriving on tyrosine-poor diet sources often develop mutualistic associations with endosymbionts able to synthesize aromatic amino acids. This surplus of aromatic amino acids is used to reinforce the insect's protective cuticle. An exceptional feature of the Sitophilus oryzae/Sodalis pierantonius interaction is the exponential increase in endosymbiotic titer observed in young adult insects, in concomitance with cuticle biosynthesis. Here, we show that host carbohydrate intake triggers endosymbiont exponential proliferation, even in conditions that lead to the detriment of the host survival. In addition, when hosts thrive on a balanced diet, endosymbiont proliferation is dispensable for several host fitness traits. The endosymbiont exponential proliferation is therefore dependent on the nutritional status of the host, and its consequences on host cuticle biosynthesis and survival depend on food quality and availability.

Metabolic relay gene of aphid and primary symbiont as RNAi target loci for aphid control

INTRODUCTION

Aphids form a stable and mutually beneficial relationship with their primary symbiont Buchnera aphidicola, which play an important role in providing the missing nutrients to the host aphid. Based on the genome sequence of wheat aphid Siotobion miscanthi and its primary symbiont Buchnera that we obtained in our previously study, we identified a metabolic relay gene, ilvA, involved in the isoleucine synthesis pathway between aphids and Buchnera.

METHOD

In this study, we identified the location and sequence structure of ilvA gene in aphid genome, the expression level in different instars and tissues of aphids, and the effect of reducing ilvA expression on the growth and development of aphids by bioinformatics analysis, quantitative PCR, RNAi and bioassay experiments.

RESULT

Our study showed that ilvA was expressed at the highest level in the 2^nd instar of the aphid, while the expression of this gene was significantly higher in the aphid bacteriocytes than in other tissues. Notably, this gene is localized on the aphid sex chromosome and remains highly conserved and collinearity across different aphid genomes. Knocking down the expression of ilvA reduced the aphid body weight and production. However, the indices of mortality decreased slightly, but were not significantly different, compared to the control.

DISCUSSION

The results show that the relay genes between aphids and their symbionts in the metabolism of essential nutrients have potential roles in the growth and development of aphids, meanwhile, providing target loci and new ideas for RNAi-based aphid green control strategies.

Facultative symbionts are potential agents of symbiont-mediated RNAi in aphids

Aphids are major crop pests, and they can be controlled through the application of the promising RNA interference (RNAi) techniques. However, chemical synthesis yield of dsRNA for RNAi is low and costly. Another sustainable aphid pest control strategy takes advantage of symbiont-mediated RNAi (SMR), which can generate dsRNA by engineered microbes. Aphid host the obligate endosymbiont Buchnera aphidicola and various facultative symbionts that not only have a wide host range but are also vertically and horizontally transmitted. Thus, we described the potential of facultative symbionts in aphid pest control by SMR. We summarized the community and host range of these facultative symbionts, and then reviewed their probable horizontal transmitted routes and ecological functions. Moreover, recent advances in the cultivation and genetic engineering of aphid facultative symbionts were discussed. In addition, current legislation of dsRNA-based pest control strategies and their safety assessments were reviewed.

Nutrition or nature: using elementary flux modes to disentangle the complex forces shaping prokaryote pan-genomes

Background

Microbial pan-genomes are shaped by a complex combination of stochastic and deterministic forces. Even closely related genomes exhibit extensive variation in their gene content. Understanding what drives this variation requires exploring the interactions of gene products with each other and with the organism’s external environment. However, to date, conceptual models of pan-genome dynamics often represent genes as independent units and provide limited information about their mechanistic interactions.

Results

We simulated the stochastic process of gene-loss using the pooled genome-scale metabolic reaction networks of 46 taxonomically diverse bacterial and archaeal families as proxies for their pan-genomes. The frequency by which reactions are retained in functional networks when stochastic gene loss is simulated in diverse environments allowed us to disentangle the metabolic reactions whose presence depends on the metabolite composition of the external environment (constrained by “nutrition”) from those that are independent of the environment (constrained by “nature”). By comparing the frequency of reactions from the first group with their observed frequencies in bacterial and archaeal families, we predicted the metabolic niches that shaped the genomic composition of these lineages. Moreover, we found that the lineages that were shaped by a more diverse metabolic niche also occur in more diverse biomes as assessed by global environmental sequencing datasets.

Conclusion

We introduce a computational framework for analyzing and interpreting pan-reactomes that provides novel insights into the ecological and evolutionary drivers of pan-genome dynamics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02052-3.

The C-Terminal Repeat Units of SpaA Mediate Adhesion of Erysipelothrix rhusiopathiae to Host Cells and Regulate Its Virulence

Simple Summary

Erysipelothrix rhusiopathiae is an important zoonotic pathogen, which poses a serious harm to the pig industry. We aimed to evaluate the genomic differences between virulent and avirulent strains to study the pathogenic mechanism of Erysipelothrix rhusiopathiae. The results showed that the spaA gene of avirulent strain lacked 120bp, encoding repeat units at the C-terminal of SpaA, the virulence of the virulent strain with this 120 bp deletion was attenuated, and the mutant strain decreased adhesion to porcine iliac artery endothelial cells.

Abstract

Erysipelothrix rhusiopathiae is a causative agent of erysipelas in animals and erysipeloid in humans. However, current information regarding E. rhusiopathiae pathogenesis remains limited. Previously, we identified two E. rhusiopathiae strains, SE38 and G4T10, which were virulent and avirulent in pigs, respectively. Here, to further study the pathogenic mechanism of E. rhusiopathiae, we sequenced and assembled the genomes of strains SE38 and G4T10, and performed a comparative genomic analysis to identify differences or mutations in virulence-associated genes. Next, we comparatively analyzed 25 E. rhusiopathiae virulence-associated genes in SE38 and G4T10. Compared with that of SE38, the spaA gene of the G4T10 strain lacked 120 bp, encoding repeat units at the C-terminal of SpaA. To examine whether these deletions or splits influence E. rhusiopathiae virulence, these 120 bp were successfully deleted from the spaA gene in strain SE38 by homologous recombination. The mutant strain ΔspaA displayed attenuated virulence in mice and decreased adhesion to porcine iliac artery endothelial cells, which was also observed using the corresponding mutant protein SpaA’. Our results demonstrate that SpaA-mediated adhesion between E. rhusiopathiae and host cells is dependent on its C-terminal repeat units.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel’s lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

The evolutionary forces driving genome size in bacteria and archaea have been subject to debate during the last decades. Typically, independent comparative analyses have suggested that unique variables, such as the strength of selection, environmental complexity, and mutation rate, are the main drivers of this trait, without considering for potential biases derived from shared ancestry. Here, we applied a phylogeny-based statistical approach to assess how tightly genome size in bacteria and archaea is linked to evolutionary history. Moreover, we also evaluated the predictability of genome size from the strength of purifying selection and ecological strategy on a broad diversity of bacteria and archaea genomes under a phylogenetic comparative framework. Our approach indicates that despite the ability of bacteria and archaea to rapidly exchange genes, a strong phylogenetic signal to genome size distributions can be recovered at broad phylogenetic scales.

Tightly Constrained Genome Reduction and Relaxation of Purifying Selection during Secondary Plastid Endosymbiosis

Endosymbiosis, the establishment of a former free-living prokaryotic or eukaryotic cell as an organelle inside a host cell, can dramatically alter the genomic architecture of the endosymbiont. Plastids or chloroplasts, the light-harvesting organelle of photosynthetic eukaryotes, are excellent models to study this phenomenon because plastid origin has occurred multiple times in evolution. Here, we investigate the genomic signature of molecular processes acting through secondary plastid endosymbiosis-the origination of a new plastid from a free-living eukaryotic alga. We used phylogenetic comparative methods to study gene loss and changes in selective regimes on plastid genomes, focusing on green algae that have given rise to three independent lineages with secondary plastids (euglenophytes, chlorarachniophytes, and Lepidodinium). Our results show an overall increase in gene loss associated with secondary endosymbiosis, but this loss is tightly constrained by the retention of genes essential for plastid function. The data show that secondary plastids have experienced temporary relaxation of purifying selection during secondary endosymbiosis. However, this process is tightly constrained, with selection relaxed only relative to the background in primary plastids. Purifying selection remains strong in absolute terms even during the endosymbiosis events. Selection intensity rebounds to pre-endosymbiosis levels following endosymbiosis events, demonstrating the changes in selection efficiency during different origin phases of secondary plastids. Independent endosymbiosis events in the euglenophytes, chlorarachniophytes, and Lepidodinium differ in their degree of relaxation of selection, highlighting the different evolutionary contexts of these events. This study reveals the selection-drift interplay during secondary endosymbiosis and evolutionary parallels during organellogenesis.

Evolutionarily recent dual obligatory symbiosis among adelgids indicates a transition between fungus- and insect-associated lifestyles

Adelgids (Insecta: Hemiptera: Adelgidae) form a small group of insects but harbor a surprisingly diverse set of bacteriocyte-associated endosymbionts, which suggest multiple replacement and acquisition of symbionts over evolutionary time. Specific pairs of symbionts have been associated with adelgid lineages specialized on different secondary host conifers. Using a metagenomic approach, we investigated the symbiosis of the Adelges laricis/Adelges tardus species complex containing betaproteobacterial ("Candidatus Vallotia tarda") and gammaproteobacterial ("Candidatus Profftia tarda") symbionts. Genomic characteristics and metabolic pathway reconstructions revealed that Vallotia and Profftia are evolutionary young endosymbionts, which complement each other's role in essential amino acid production. Phylogenomic analyses and a high level of genomic synteny indicate an origin of the betaproteobacterial symbiont from endosymbionts of Rhizopus fungi. This evolutionary transition was accompanied with substantial loss of functions related to transcription regulation, secondary metabolite production, bacterial defense mechanisms, host infection, and manipulation. The transition from fungus to insect endosymbionts extends our current framework about evolutionary trajectories of host-associated microbes.

Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

An assessment of genome annotation coverage across the bacterial tree of life

Although gene-finding in bacterial genomes is relatively straightforward, the automated assignment of gene function is still challenging, resulting in a vast quantity of hypothetical sequences of unknown function. But how prevalent are hypothetical sequences across bacteria, what proportion of genes in different bacterial genomes remain unannotated, and what factors affect annotation completeness? To address these questions, we surveyed over 27 000 bacterial genomes from the Genome Taxonomy Database, and measured genome annotation completeness as a function of annotation method, taxonomy, genome size, 'research bias' and publication date. Our analysis revealed that 52 and 79 % of the average bacterial proteome could be functionally annotated based on protein and domain-based homology searches, respectively. Annotation coverage using protein homology search varied significantly from as low as 14 % in some species to as high as 98 % in others. We found that taxonomy is a major factor influencing annotation completeness, with distinct trends observed across the microbial tree (e.g. the lowest level of completeness was found in the Patescibacteria lineage). Most lineages showed a significant association between genome size and annotation incompleteness, likely reflecting a greater degree of uncharacterized sequences in 'accessory' proteomes than in 'core' proteomes. Finally, research bias, as measured by publication volume, was also an important factor influencing genome annotation completeness, with early model organisms showing high completeness levels relative to other genomes in their own taxonomic lineages. Our work highlights the disparity in annotation coverage across the bacterial tree of life and emphasizes a need for more experimental characterization of accessory proteomes as well as understudied lineages.

Growing Ungrowable Bacteria: Overview and Perspectives on Insect Symbiont Culturability

Insects are often involved in endosymbiosis, that is, the housing of symbiotic microbes within their tissues or within their cells. Endosymbionts are a major driving force in insects' evolution, because they dramatically affect their host physiology and allow them to adapt to new niches, for example, by complementing their diet or by protecting them against pathogens. Endosymbiotic bacteria are, however, fastidious and therefore difficult to manipulate outside of their hosts, especially intracellular species. The coevolution between hosts and endosymbionts leads to alterations in the genomes of endosymbionts, limiting their ability to cope with changing environments. Consequently, few insect endosymbionts are culturable in vitro and genetically tractable, making functional genetics studies impracticable on most endosymbiotic bacteria. However, recently, major progress has been made in manipulating several intracellular endosymbiont species in vitro, leading to astonishing discoveries on their physiology and the way they interact with their host. This review establishes a comprehensive picture of the in vitro tractability of insect endosymbiotic bacteria and addresses the reason why most species are not culturable. By compiling and discussing the latest developments in the design of custom media and genetic manipulation protocols, it aims at providing new leads to expand the range of tractable endosymbionts and foster genetic research on these models.

Characterization of a novel Pantoea symbiont allows inference of a pattern of convergent genome reduction in bacteria associated with Pentatomidae

Phytophagous stink bugs typically harbor nutritional symbiotic bacteria in their midgut, to integrate their unbalanced diet. In the Pentatomidae, most symbionts are affiliated to the genus Pantoea, and are polyphyletic. This suggests a scenario of an ancestral establishment of symbiosis, followed by multiple symbiont replacement events by akin environmental bacteria in different host lineages. In this study, a novel Pantoeaspecies ('CandidatusPantoea persica') was characterized from the gut of the pentatomid Acrosternum arabicum, and shown to be highly abundant in a specific portion of the gut and necessary for the host development. The genome of the symbiont (2.9 Mb), while presenting putative host-supportive metabolic pathways, including those for amino acids and vitamin synthesis, showed a high level of pseudogenization, indicating ongoing genome reduction. Comparative analyses with other free-living and symbiotic Pantoea highlighted a convergent pattern of genome reduction in symbionts of pentatomids, putatively following the typical phases modelized in obligate nutritional symbionts of insects. Additionally, this system has distinctive traits, as hosts are closely related, and symbionts originated multiple independent times from closely related free-living bacteria, displaying convergent and independent conspicuous genome reduction. Due to such peculiarities, this may become an ideal model to study genome evolutionary processes in insect symbionts.

The study of host-microbiome (co)evolution across levels of selection

Microorganismal diversity can be explained in large part by selection imposed from both the abiotic and biotic environments, including-in the case of host-associated microbiomes-interactions with eukaryotes. As such, the diversity of host-associated microbiomes can be usefully studied across a variety of scales: within a single host over time, among host genotypes within a population, between populations and among host species. A plethora of recent studies across these scales and across diverse systems are: (i) exemplifying the importance of the host genetics in shaping microbiome composition; (ii) uncovering the role of the microbiome in shaping key host phenotypes; and (iii) highlighting the dynamic nature of the microbiome. They have also raised a critical question: do these complex associations fit within our existing understanding of evolution and coevolution, or do these often intimate and seemingly cross-generational interactions follow novel evolutionary rules from those previously identified? Herein, we describe the known importance of (co)evolution in host-microbiome systems, placing the existing data within extant frameworks that have been developed over decades of study, and ask whether there are unique properties of host-microbiome systems that require a paradigm shift. By examining when and how selection can act on the host and its microbiome as a unit (termed, the holobiont), we find that the existing conceptual framework, which focuses on individuals, as well as interactions among individuals and groups, is generally well suited for understanding (co)evolutionary change in these intimate assemblages. This article is part of the theme issue 'The role of the microbiome in host evolution'.

Quorum sensing sets the stage for the establishment and vertical transmission of Sodalis praecaptivus in tsetse flies

Bacterial virulence factors facilitate host colonization and set the stage for the evolution of parasitic and mutualistic interactions. The Sodalis-allied clade of bacteria exhibit striking diversity in the range of both plant and animal feeding insects they inhabit, suggesting the appropriation of universal molecular mechanisms that facilitate establishment. Here, we report on the infection of the tsetse fly by free-living Sodalis praecaptivus, a close relative of many Sodalis-allied symbionts. Key genes involved in quorum sensing, including the homoserine lactone synthase (ypeI) and response regulators (yenR and ypeR) are integral for the benign colonization of S. praecaptivus. Mutants lacking ypeI, yenR and ypeR compromised tsetse survival as a consequence of their inability to repress virulence. Genes under quorum sensing, including homologs of the binary insecticidal toxin PirAB and a putative symbiosis-promoting factor CpmAJ, demonstrated negative and positive impacts, respectively, on tsetse survival. Taken together with results obtained from experiments involving weevils, this work shows that quorum sensing virulence suppression plays an integral role in facilitating the establishment of Sodalis-allied symbionts in diverse insect hosts. This knowledge contributes to the understanding of the early evolutionary steps involved in the formation of insect-bacterial symbiosis. Further, despite having no established history of interaction with tsetse, S. praecaptivus can infect reproductive tissues, enabling vertical transmission through adenotrophic viviparity within a single host generation. This creates an option for the use of S. praecaptivus in the biocontrol of insect disease vectors via paratransgenesis.

RAD51 Gene Family Structure and Function

Accurate DNA repair and replication are critical for genomic stability and cancer prevention. RAD51 and its gene family are key regulators of DNA fidelity through diverse roles in double-strand break repair, replication stress, and meiosis. RAD51 is an ATPase that forms a nucleoprotein filament on single-stranded DNA. RAD51 has the function of finding and invading homologous DNA sequences to enable accurate and timely DNA repair. Its paralogs, which arose from ancient gene duplications of RAD51, have evolved to regulate and promote RAD51 function. Underscoring its importance, misregulation of RAD51, and its paralogs, is associated with diseases such as cancer and Fanconi anemia. In this review, we focus on the mammalian RAD51 structure and function and highlight the use of model systems to enable mechanistic understanding of RAD51 cellular roles. We also discuss how misregulation of the RAD51 gene family members contributes to disease and consider new approaches to pharmacologically inhibit RAD51.

Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

Evolutionary dynamics of the chromatophore genome in three photosynthetic Paulinella species

The thecate amoeba Paulinella is a valuable model for understanding plastid organellogenesis because this lineage has independently gained plastids (termed chromatophores) of alpha-cyanobacterial provenance. Plastid primary endosymbiosis in Paulinella occurred relatively recently (90–140 million years ago, Mya), whereas the origin of the canonical Archaeplastida plastid occurred >1,500 Mya. Therefore, these two events provide independent perspectives on plastid formation on vastly different timescales. Here we generated the complete chromatophore genome sequence from P. longichromatophora (979,356 bp, GC-content = 38.8%, 915 predicted genes) and P. micropora NZ27 (977,190 bp, GC-content = 39.9%, 911 predicted genes) and compared these data to that from existing chromatophore genomes. Our analysis suggests that when a basal split occurred among photosynthetic Paulinella species ca. 60 Mya, only 35% of the ancestral orthologous gene families from the cyanobacterial endosymbiont remained in chromatophore DNA. Following major gene losses during the early stages of endosymbiosis, this process slowed down significantly, resulting in a conserved gene content across extant taxa. Chromatophore genes faced relaxed selection when compared to homologs in free-living alpha-cyanobacteria, likely reflecting the homogeneous intracellular environment of the Paulinella host. Comparison of nucleotide substitution and insertion/deletion events among different P. micropora strains demonstrates that increases in AT-content and genome reduction are ongoing and dynamic processes in chromatophore evolution.

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.

Functions of the Microbiota for the Physiology of Animal Metaorganisms

Animals are usually regarded as independent entities within their respective environments. However, within an organism, eukaryotes and prokaryotes interact dynamically to form the so-called metaorganism or holobiont, where each partner fulfils its versatile and crucial role. This review focuses on the interplay between microorganisms and multicellular eukaryotes in the context of host physiology, in particular aging and mucus-associated crosstalk. In addition to the interactions between bacteria and the host, we highlight the importance of viruses and nonmodel organisms. Moreover, we discuss current culturing and computational methodologies that allow a deeper understanding of underlying mechanisms controlling the physiology of metaorganisms.

The Domestication of a Large DNA Virus by the Wasp Venturia canescens Involves Targeted Genome Reduction through Pseudogenization

Polydnaviruses (PDVs) are compelling examples of viral domestication, in which wasps express a large set of genes originating from a chromosomally integrated virus to produce particles necessary for their reproductive success. Parasitoid wasps generally use PDVs as a virulence gene delivery system allowing the protection of their progeny in the body of parasitized host. However, in the wasp Venturia canescens an independent viral domestication process led to an alternative strategy as the wasp incorporates virulence proteins in viral liposomes named virus-like particles (VLPs), instead of DNA molecules. Proteomic analysis of purified VLPs and transcriptome sequencing revealed the loss of some viral functions. In particular, the genes coding for capsid components are no longer expressed, which explains why VLPs do not incorporate DNA. Here a thorough examination of V. canescens genome revealed the presence of the pseudogenes corresponding to most of the genes involved in lost functions. This strongly suggests that an accumulation of mutations that leads to gene specific pseudogenization precedes the loss of viral genes observed during virus domestication. No evidence was found for block loss of collinear genes, although extensive gene order reshuffling of the viral genome was identified from comparisons between endogenous and exogenous viruses. These results provide the first insights on the early stages of large DNA virus domestication implicating massive genome reduction through gene-specific pseudogenization, a process which differs from the large deletions described for bacterial endosymbionts.

Random Genetic Drift and Selective Pressures Shaping the Blattabacterium Genome

Estimates suggest that at least half of all extant insect genera harbor obligate bacterial mutualists. Whereas an endosymbiotic relationship imparts many benefits upon host and symbiont alike, the intracellular lifestyle has profound effects on the bacterial genome. The obligate endosymbiont genome is a product of opposing forces: genes important to host survival are maintained through physiological constraint, contrasted by the fixation of deleterious mutations and genome erosion through random genetic drift. The obligate cockroach endosymbiont, Blattabacterium - providing nutritional augmentation to its host in the form of amino acid synthesis - displays radical genome alterations when compared to its most recent free-living relative Flavobacterium. To date, eight Blattabacterium genomes have been published, affording an unparalleled opportunity to examine the direction and magnitude of selective forces acting upon this group of symbionts. Here, we find that the Blattabacterium genome is experiencing a 10-fold increase in selection rate compared to Flavobacteria. Additionally, the proportion of selection events is largely negative in direction, with only a handful of loci exhibiting signatures of positive selection. These findings suggest that the Blattabacterium genome will continue to erode, potentially resulting in an endosymbiont with an even further reduced genome, as seen in other insect groups such as Hemiptera.

mTOR Complex 1 Implicated in Aphid/Buchnera Host/Symbiont Integration

Obligate nutritional endosymbioses are arguably the most intimate of all interspecific associations. While many insect nutritional endosymbioses are well studied, a full picture of how two disparate organisms, a bacterial endosymbiont and a eukaryotic host, are integrated is still lacking. The mTOR pathway is known to integrate nutritional conditions with cell growth and survival in eukaryotes. Characterization and localization of amino acid transporters in aphids suggest the mTOR pathway as a point of integration between an aphid host and its amino acid-provisioning endosymbiont Buchnera aphidicola. The mTOR pathway is unannotated in aphids and unstudied in any nutritional endosymbiosis. We annotated mTOR pathway genes in two aphid species, Acyrthosiphon pisum and Myzus persicae, using both BLASTp searches and Hidden Markov Models. Using previously collected RNAseq data we constructed new reference transcriptomes for bacteriocyte, gut, and whole insect tissue for three lines of M. persicae. Annotation of the mTOR pathway identified homologs of all known invertebrate mTOR genes in both aphid species with some duplications. Differential expression analysis showed that genes specific to the amino acid-sensitive mTOR Complex 1 were more highly expressed in bacteriocytes than genes specific to the amino acid-insensitive mTOR Complex 2. Almost all mTOR genes involved in sensing amino acids showed higher expression in bacteriocytes than in whole insect tissue. When compared to gut, the putative glutamine/arginine sensing transporter ACYPI000333, an ortholog of SLC38A9, showed 6.5 times higher expression in bacteriocytes. Our results suggest that the mTOR pathway may be functionally important in mediating integration of Buchnera into aphid growth and reproduction.

Biology of Fungi and Their Bacterial Endosymbionts

Heritable symbioses, in which endosymbiotic bacteria (EB) are transmitted vertically between host generations, are an important source of evolutionary novelties. A primary example of such symbioses is the eukaryotic cell with its EB-derived organelles. Recent discoveries suggest that endosymbiosis-related innovations can be also found in associations formed by early divergent fungi in the phylum Mucoromycota with heritable EB from two classes, Betaproteobacteria and Mollicutes. These symbioses exemplify novel types of host-symbiont interactions. Studies of these partnerships fuel theoretical models describing mechanisms that stabilize heritable symbioses, control the rate of molecular evolution, and enable the establishment of mutualisms. Lastly, by altering host phenotypes and metabolism, these associations represent an important instrument for probing the basic biology of the Mucoromycota hosts, which remain one of the least explored filamentous fungi.

Coexistence of novel gammaproteobacterial and Arsenophonus symbionts in the scale insect Greenisca brachypodii (Hemiptera, Coccomorpha: Eriococcidae)

Scale insects are commonly associated with obligate, intracellular microorganisms which play important roles in complementing their hosts with essential nutrients. Here we characterized the symbiotic system of Greenisca brachypodii, a member of the family Eriococcidae. Histological and ultrastructural analyses have indicated that G. brachypodii is stably associated with coccoid and rod-shaped bacteria. Phylogenetic analyses have revealed that the coccoid bacteria represent a sister group to the secondary symbiont of the mealybug Melanococcus albizziae, whereas the rod-shaped symbionts are close relatives of Arsenophonus symbionts in insects - to our knowledge, this is the first report of the presence of Arsenophonus bacterium in scale insects. As a comparison of 16S and 23S rRNA genes sequences of the G. brachypodii coccoid symbiont with other gammaprotebacterial sequences showed only low similarity (∼90%), we propose the name 'Candidatus Kotejella greeniscae' for its tentative classification. Both symbionts are transovarially transmitted from one generation to the next. The infection takes place in the neck region of the ovariole. The bacteria migrate between follicular cells, as well as through the cytoplasm of those cells to the perivitelline space, where they form a characteristic 'symbiont ball'. Our findings provide evidence for a polyphyletic origin of symbionts of Eriococcidae.

Non-Random Inversion Landscapes in Prokaryotic Genomes Are Shaped by Heterogeneous Selection Pressures

Inversions are a major contributor to structural genome evolution in prokaryotes. Here, using a novel alignment-based method, we systematically compare 1,651 bacterial and 98 archaeal genomes to show that inversion landscapes are frequently biased toward (symmetric) inversions around the origin–terminus axis. However, symmetric inversion bias is not a universal feature of prokaryotic genome evolution but varies considerably across clades. At the extremes, inversion landscapes in Bacillus–Clostridium and Actinobacteria are dominated by symmetric inversions, while there is little or no systematic bias favoring symmetric rearrangements in archaea with a single origin of replication. Within clades, we find strong but clade-specific relationships between symmetric inversion bias and different features of adaptive genome architecture, including the distance of essential genes to the origin of replication and the preferential localization of genes on the leading strand. We suggest that heterogeneous selection pressures have converged to produce similar patterns of structural genome evolution across prokaryotes.

A sRNA in a reduced mutualistic symbiont genome regulates its own gene expression

Similar to other nutritional endosymbionts that are obligate for host survival, the mutualistic aphid endosymbiont, Buchnera, has a highly reduced genome with few regulatory elements. Until recently, it was thought that aphid hosts were primarily responsible for regulating their symbiotic relationship. However, we recently revealed that Buchnera displays differential protein regulation, but not mRNA expression. We also identified a number of conserved small RNAs (sRNAs) that are expressed among Buchnera taxa. In this study, we investigate whether differential protein regulation in Buchnera is the result of post-transcriptional gene regulation via sRNAs. We characterize the sRNA profile of two Buchnera life stages: (i) when Buchnera is transitioning from an extracellular proliferating state in aphid embryos and (ii) when Buchnera is in an intracellular nonproliferating state in aphid bacteriocytes (specialized symbiont cells). Overall, we identified 90 differentially expressed sRNAs, 97% of which were upregulated in aphid embryos. Of these sRNAs, the majority were predicted to be involved in the regulation of various metabolic processes, including arginine biosynthesis. Using a heterologous dual expression vector, we reveal for the first time that a Buchnera antisense sRNA can post-transcriptionally interact with its cognate Buchnera coding sequence, carB, a gene involved in arginine biosynthesis. These results corroborate our in vivo RNAseq and proteomic data, where the candidate antisense sRNA carB and the protein CarB are significantly upregulated in aphid embryos. Overall, we demonstrate that Buchnera may regulate gene expression independently from its host by utilizing sRNAs.

Repeated replacement of an intrabacterial symbiont in the tripartite nested mealybug symbiosis

Stable endosymbiosis of a bacterium into a host cell promotes cellular and genomic complexity. The mealybug Planococcus citri has two bacterial endosymbionts with an unusual nested arrangement: the γ-proteobacterium Moranella endobia lives in the cytoplasm of the β-proteobacterium Tremblaya princeps These two bacteria, along with genes horizontally transferred from other bacteria to the P. citri genome, encode gene sets that form an interdependent metabolic patchwork. Here, we test the stability of this three-way symbiosis by sequencing host and symbiont genomes for five diverse mealybug species and find marked fluidity over evolutionary time. Although Tremblaya is the result of a single infection in the ancestor of mealybugs, the γ-proteobacterial symbionts result from multiple replacements of inferred different ages from related but distinct bacterial lineages. Our data show that symbiont replacement can happen even in the most intricate symbiotic arrangements and that preexisting horizontally transferred genes can remain stable on genomes in the face of extensive symbiont turnover.

Losing Complexity: The Role of Simplification in Macroevolution

Macroevolutionary patterns can be produced by combinations of diverse and even oppositional dynamics. A growing body of data indicates that secondary simplifications of molecular and cellular structures are common. Some major diversifications in eukaryotes have occurred because of loss and minimalisation; numerous episodes in prokaryote evolution have likewise been driven by the reduction of structure. After examining a range of examples of secondary simplification and its consequences across the tree of life, we address how macroevolutionary explanations might incorporate simplification as well as complexification, and adaptive as well as nonadaptive dynamics.

Loss of Conserved Noncoding RNAs in Genomes of Bacterial Endosymbionts

The genomes of intracellular symbiotic or pathogenic bacteria, such as of Buchnera, Mycoplasma, and Rickettsia, are typically smaller compared with their free-living counterparts. Here we showed that noncoding RNA (ncRNA) families, which are conserved in free-living bacteria, frequently could not be detected by computational methods in the small genomes. Statistical tests demonstrated that their absence is not an artifact of low GC content or small deletions in these small genomes, and thus it was indicative of an independent loss of ncRNAs in different endosymbiotic lineages. By analyzing the synteny (conservation of gene order) between the reduced and nonreduced genomes, we revealed instances of protein-coding genes that were preserved in the reduced genomes but lost cis-regulatory elements. We found that the loss of cis-regulatory ncRNA sequences, which regulate the expression of cognate protein-coding genes, is characterized by the reduction of secondary structure formation propensity, GC content, and length of the corresponding genomic regions.

Shifts in Buchnera aphidicola density in soybean aphids (Aphis glycines) feeding on virus-infected soybean

Vertically transmitted bacterial symbionts are common in arthropods. Aphids undergo an obligate symbiosis with Buchnera aphidicola, which provides essential amino acids to its host and contributes directly to nymph growth and reproduction. We previously found that newly adult Aphis glycines feeding on soybean infected with the beetle-transmitted Bean pod mottle virus (BPMV) had significantly reduced fecundity. We hypothesized that the reduced fecundity was attributable to detrimental impacts of the virus on the aphid microbiome, namely Buchnera. To test this, mRNA sequencing and quantitative real-time PCR were used to assay Buchnera transcript abundance and titre in A. glycines feeding on Soybean mosaic virus-infected, BPMV-infected, and healthy soybean for up to 14 days. Our results indicated that Buchnera density was lower and ultimately suppressed in aphids feeding on virus-infected soybean. While the decreased Buchnera titre may be associated with reduced aphid fecundity, additional mechanisms are probably involved. The present report begins to describe how interactions among insects, plants, and plant pathogens influence endosymbiont population dynamics.

Endosymbiont evolution: predictions from theory and surprises from genomes

Genome data have created new opportunities to untangle evolutionary processes shaping microbial variation. Among bacteria, long-term mutualists of insects represent the smallest and (typically) most AT-rich genomes. Evolutionary theory provides a context to predict how an endosymbiotic lifestyle may alter fundamental evolutionary processes--mutation, selection, genetic drift, and recombination--and thus contribute to extreme genomic outcomes. These predictions can then be explored by comparing evolutionary rates, genome size and stability, and base compositional biases across endosymbiotic and free-living bacteria. Recent surprises from such comparisons include genome reduction among uncultured, free-living species. Some studies suggest that selection generally drives this streamlining, while drift drives genome reduction in endosymbionts; however, this remains an hypothesis requiring additional data. Unexpected evidence of selection acting on endosymbiont GC content hints that even weak selection may be effective in some long-term mutualists. Moving forward, intraspecific analysis offers a promising approach to distinguish underlying mechanisms, by testing the null hypothesis of neutrality and by quantifying mutational spectra. Such analyses may clarify whether endosymbionts and free-living bacteria occupy distinct evolutionary trajectories or, alternatively, represent varied outcomes of similar underlying forces.

A molecular survey of Australian and North American termite genera indicates that vertical inheritance is the primary force shaping termite gut microbiomes

BACKGROUND

Termites and their microbial gut symbionts are major recyclers of lignocellulosic biomass. This important symbiosis is obligate but relatively open and more complex in comparison to other well-known insect symbioses such as the strict vertical transmission of Buchnera in aphids. The relative roles of vertical inheritance and environmental factors such as diet in shaping the termite gut microbiome are not well understood.

RESULTS

The gut microbiomes of 66 specimens representing seven higher and nine lower termite genera collected in Australia and North America were profiled by small subunit (SSU) rRNA amplicon pyrosequencing. These represent the first reported culture-independent gut microbiome data for three higher termite genera: Tenuirostritermes, Drepanotermes, and Gnathamitermes; and two lower termite genera: Marginitermes and Porotermes. Consistent with previous studies, bacteria comprise the largest fraction of termite gut symbionts, of which 11 phylotypes (6 Treponema, 1 Desulfarculus-like, 1 Desulfovibrio, 1 Anaerovorax-like, 1 Sporobacter-like, and 1 Pirellula-like) were widespread occurring in ≥50% of collected specimens. Archaea are generally considered to comprise only a minority of the termite gut microbiota (<3%); however, archaeal relative abundance was substantially higher and variable in a number of specimens including Macrognathotermes, Coptotermes, Schedorhinotermes, Porotermes, and Mastotermes (representing up to 54% of amplicon reads). A ciliate related to Clevelandella was detected in low abundance in Gnathamitermes indicating that protists were either reacquired after protists loss in higher termites or persisted in low numbers across this transition. Phylogenetic analyses of the bacterial communities indicate that vertical inheritance is the primary force shaping termite gut microbiota. The effect of diet is secondary and appears to influence the relative abundance, but not membership, of the gut communities.

CONCLUSIONS

Vertical inheritance is the primary force shaping the termite gut microbiome indicating that species are successfully and faithfully passed from one generation to the next via trophallaxis or coprophagy. Changes in relative abundance can occur on shorter time scales and appear to be an adaptive mechanism for dietary fluctuations.

Taking the pseudo out of pseudogenes

Pseudogenes are defined as fragments of once-functional genes that have been silenced by one or more nonsense, frameshift or missense mutations. Despite continuing increases in the speed of sequencing and annotating bacterial genomes, the identification and categorisation of pseudogenes remains problematic. Even when identified, pseudogenes are considered to be rare and tend to be ignored. On the contrary, pseudogenes are surprisingly prevalent and can persist for long evolutionary time periods, representing a record of once-functional genetic characteristics. Most importantly, pseudogenes provide an insight into prokaryotic evolutionary history as a record of phenotypic traits that have been lost. Focusing on the intracellular and symbiotic bacteria in which pseudogenes predominate, this review discusses the importance of identifying pseudogenes to fully understand the abilities of bacteria, and to understand prokaryotes within their evolutionary context.

Widespread expression of conserved small RNAs in small symbiont genomes

Genome architecture of a microbe markedly changes when it transitions from a free-living lifestyle to an obligate symbiotic association within eukaryotic cells. These symbiont genomes experience numerous rearrangements and massive gene loss, which is expected to radically alter gene regulatory networks compared with those of free-living relatives. As such, it remains unclear whether and how these small symbiont genomes regulate gene expression. Here, using a label-free mass-spec quantification approach we found that differential protein regulation occurs in Buchnera, a model symbiont with a reduced genome, when it transitions between two distinct life stages. However, differential mRNA expression could not be detected between Buchnera life stages, despite the presence of a small number of putative transcriptional regulators. Instead a comparative analysis of small RNA expression profiles among five divergent Buchnera lineages, spanning a variety of Buchnera life stages, reveals 140 novel intergenic and antisense small RNAs and 517 untranslated regions that were significantly expressed, some of which have been conserved for ∼65 million years. In addition, the majority of these small RNAs exhibit both sequence covariation and thermodynamic stability, indicators of a potential structural RNA role. Together, these data suggest that gene regulation at the post-transcriptional level may be important in Buchnera. This is the first study to empirically identify Buchnera small RNAs, and we propose that these novel small RNAs may facilitate post-transcriptional regulation through translational inhibition/activation, and/or transcript stability. Ultimately, post-transcriptional regulation may shape metabolic complementation between Buchnera and its aphid host, thus impacting the animal's ecology and evolution.

Natural selection for operons depends on genome size

In prokaryotes, genome size is associated with metabolic versatility, regulatory complexity, effective population size, and horizontal transfer rates. We therefore analyzed the covariation of genome size and operon conservation to assess the evolutionary models of operon formation and maintenance. In agreement with previous results, intraoperonic pairs of essential and of highly expressed genes are more conserved. Interestingly, intraoperonic pairs of genes are also more conserved when they encode proteins at similar cell concentrations, suggesting a role of cotranscription in diminishing the cost of waste and shortfall in gene expression. Larger genomes have fewer and smaller operons that are also less conserved. Importantly, lower conservation in larger genomes was observed for all classes of operons in terms of gene expression, essentiality, and balanced protein concentration. We reached very similar conclusions in independent analyses of three major bacterial clades (α- and β-Proteobacteria and Firmicutes). Operon conservation is inversely correlated to the abundance of transcription factors in the genome when controlled for genome size. This suggests a negative association between the complexity of genetic networks and operon conservation. These results show that genome size and/or its proxies are key determinants of the intensity of natural selection for operon organization. Our data fit better the evolutionary models based on the advantage of coregulation than those based on genetic linkage or stochastic gene expression. We suggest that larger genomes with highly complex genetic networks and many transcription factors endure weaker selection for operons than smaller genomes with fewer alternative tools for genetic regulation.

Signs of neutralization in a redundant gene involved in homologous recombination in Wolbachia endosymbionts

Genomic reduction in bacterial endosymbionts occurs through large genomic deletions and long-term accumulation of mutations. The latter process involves successive steps including gene neutralization, pseudogenization, and gradual erosion until complete loss. Although many examples of pseudogenes at various levels of degradation have been reported, neutralization cases are scarce because of the transient nature of the process. Gene neutralization may occur due to relaxation of selection in nonessential genes, for example, those involved in redundant functions. Here, we report an example of gene neutralization in the homologous recombination (HR) pathway of Wolbachia, a bacterial endosymbiont of arthropods and nematodes. The HR pathway is often depleted in endosymbiont genomes, but it is apparently intact in some Wolbachia strains. Analysis of 12 major HR genes showed that they have been globally under strong purifying selection during the evolution of Wolbachia strains hosted by arthropods, supporting the evolutionary importance of the HR pathway for these Wolbachia genomes. However, we detected signs of recent neutralization of the ruvA gene in a subset of Wolbachia strains, which might be related to an ancestral, clade-specific amino acid change that impaired DNA-binding activity. Strikingly, RuvA is part of the RuvAB complex involved in branch migration, whose function overlaps with the RecG helicase. Although ruvA is experiencing neutralization, recG is under strong purifying selection. Thus, our high phylogenetic resolution suggests that we identified a rare example of targeted neutralization of a gene involved in a redundant function in an endosymbiont genome.

Comparative genomics reveals surprising divergence of two closely related strains of uncultivated UCYN-A cyanobacteria

Marine planktonic cyanobacteria capable of fixing molecular nitrogen (termed 'diazotrophs') are key in biogeochemical cycling, and the nitrogen fixed is one of the major external sources of nitrogen to the open ocean. Candidatus Atelocyanobacterium thalassa (UCYN-A) is a diazotrophic cyanobacterium known for its widespread geographic distribution in tropical and subtropical oligotrophic oceans, unusually reduced genome and symbiosis with a single-celled prymnesiophyte alga. Recently a novel strain of this organism was also detected in coastal waters sampled from the Scripps Institute of Oceanography pier. We analyzed the metagenome of this UCYN-A2 population by concentrating cells by flow cytometry. Phylogenomic analysis provided strong bootstrap support for the monophyly of UCYN-A (here called UCYN-A1) and UCYN-A2 within the marine Crocosphaera sp. and Cyanothece sp. clade. UCYN-A2 shares 1159 of the 1200 UCYN-A1 protein-coding genes (96.6%) with high synteny, yet the average amino-acid sequence identity between these orthologs is only 86%. UCYN-A2 lacks the same major pathways and proteins that are absent in UCYN-A1, suggesting that both strains can be grouped at the same functional and ecological level. Our results suggest that UCYN-A1 and UCYN-A2 had a common ancestor and diverged after genome reduction. These two variants may reflect adaptation of the host to different niches, which could be coastal and open ocean habitats.