50 million years of genomic stasis in endosymbiotic bacteria

Twenty-nine newly sequenced genomes and a comprehensive genome dataset for the insect endosymbiont Buchnera

Most phloem-feeding insects face nutritional deficiency and rely on their intracellular symbionts to provide nutrients, and most of endosymbiont genomes have undergone reduction. However, the study of genome reduction processes of endosymbionts has been constrained by the limited availability of genome data from different insect lineages. The obligate relationship between aphids and Buchnera aphidicola (hereafter Buchnera) makes them a classic model for studying insect-endosymbiont interaction. Here, we report 29 newly sequenced Buchnera genomes from 11 aphid subfamilies, and a comprehensive dataset based on 90 Buchnera genomes from 14 aphid subfamilies. The dataset shows a significant genomic difference of Buchnera among different aphid lineages. The dataset exhibits a more balanced distribution of Buchnera (from 14 aphid subfamilies) genome sizes, ranging from 400 kb to 600 kb, which can illustrate the genome reduction process of Buchnera. The new genome data provide valuable insights into the microevolutionary processes leading to genomic reduction of insect endosymbionts.

Comparative genomics of the primary endosymbiont Buchnera aphidicola in aphid hosts and their coevolutionary relationships

BACKGROUND

Coevolution between modern aphids and their primary obligate, bacterial endosymbiont, Buchnera aphidicola, has been previously reported at different classification levels based on molecular phylogenetic analyses. However, the Buchnera genome remains poorly understood within the Rhus gall aphids.

RESULTS

We assembled the complete genome of the endosymbiont Buchnera in 16 aphid samples, representing 13 species in all six genera of Rhus gall aphids by shotgun genome skimming method. We compared the newly assembled genomes with those from GenBank to comprehensively investigate patterns of coevolution between the bacteria Buchnera and their aphid hosts. Buchnera genomes were mostly collinear, and the pan-genome contained 684 genes, in which the core genome contained 256 genes with some lineages having large numbers of tandem gene duplications. There has been substantial gene-loss in each Buchnera lineage. We also reconstructed the phylogeny for Buchnera and their host aphids, respectively, using 72 complete genomes of Buchnera, along with the complete mitochondrial genomes and three nuclear genes of 31 corresponding host aphid accessions. The cophylogenetic test demonstrated significant coevolution between these two partner groups at individual, species, generic, and tribal levels.

CONCLUSIONS

Buchnera exhibits very high levels of genomic sequence divergence but relative stability in gene order. The relationship between the symbionts Buchnera and its aphid hosts shows a significant coevolutionary pattern and supports complexity of the obligate symbiotic relationship.

Host-specific co-evolution likely driven by diet in Buchnera aphidicola

BACKGROUND

Russian wheat aphid (Diuraphis noxia Kurd.) is a severe pest to wheat, and even though resistance varieties are available to curb this pest, they are becoming obsolete with the development of new virulent aphid populations. Unlike many other aphids, D noxia only harbours a single endosymbiont, Buchnera aphidicola. Considering the importance of Buchnera, this study aimed to elucidate commonalities and dissimilarities between various hosts, to better understand its distinctiveness within its symbiotic relationship with D. noxia. To do so, the genome of the D. noxia's Buchnera was assembled and compared to those of other aphid species that feed on diverse host species.

RESULTS

The overall importance of several features such as gene length and percentage GC content was found to be critical for the maintenance of Buchnera genes when compared to their closest free-living relative, Escherichia coli. Buchnera protein coding genes were found to have percentage GC contents that tended towards a mean of ~ 26% which had strong correlation to their identity to their E. coli homologs. Several SNPs were identified between different aphid populations and multiple isolates of Buchnera were confirmed in single aphids.

CONCLUSIONS

Establishing the strong correlation of percentage GC content of protein coding genes and gene identity will allow for identifying which genes will be lost in the continually shrinking Buchnera genome. This is also the first report of a parthenogenically reproducing aphid that hosts multiple Buchnera strains in a single aphid, raising questions regarding the benefits of maintaining multiple strains. We also found preliminary evidence for post-transcriptional regulation of Buchnera genes in the form of polyadenylation.

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.

Highly Reduced Complementary Genomes of Dual Bacterial Symbionts in the Mulberry Psyllid Anomoneura mori

The genomes of obligately host-restricted bacteria suffer from accumulating mildly deleterious mutations, resulting in marked size reductions. Psyllids (Hemiptera) are phloem sap-sucking insects with a specialized organ called the bacteriome, which typically harbors two vertically transmitted bacterial symbionts: the primary symbiont "Candidatus Carsonella ruddii" (Gammaproteobacteria) and a secondary symbiont that is phylogenetically diverse among psyllid lineages. The genomes of several Carsonella lineages were revealed to be markedly reduced (158-174‍ ‍kb), AT-rich (14.0-17.9% GC), and structurally conserved with similar gene inventories devoted to synthesizing essential amino acids that are scarce in the phloem sap. However, limited genomic information is currently available on secondary symbionts. Therefore, the present study investigated the genomes of the bacteriome-associated dual symbionts, Secondary_AM (Gammaproteobacteria) and Carsonella_AM, in the mulberry psyllid Anomoneura mori (Psyllidae). The results obtained revealed that the Secondary_AM genome is as small and AT-rich (229,822 bp, 17.3% GC) as those of Carsonella lineages, including Carsonella_AM (169,120 bp, 16.2% GC), implying that Secondary_AM is an evolutionarily ancient obligate mutualist, as is Carsonella. Phylogenomic ana-lyses showed that Secondary_AM is sister to "Candidatus Psyllophila symbiotica" of Cacopsylla spp. (Psyllidae), the genomes of which were recently reported (221-237‍ ‍kb, 17.3-18.6% GC). The Secondary_AM and Psyllophila genomes showed highly conserved synteny, sharing all genes for complementing the incomplete tryptophan biosynthetic pathway of Carsonella and those for synthesizing B vitamins. However, sulfur assimilation and carotenoid-synthesizing genes were only retained in Secondary_AM and Psyllophila, respectively, indicating ongoing gene silencing. Average nucleotide identity, gene ortholog similarity, genome-wide synteny, and substitution rates suggest that the Secondary_AM/Psyllophila genomes are more labile than Carsonella genomes.

Nutrient supplementation by genome-eroded Burkholderia symbionts of scale insects

Hemipterans are known as hosts to bacterial or fungal symbionts that supplement their unbalanced diet with essential nutrients. Among them, scale insects (Coccomorpha) are characterized by a particularly large diversity of symbiotic systems. Here, using microscopic and genomic approaches, we functionally characterized the symbionts of two scale insects belonging to the Eriococcidae family, Acanthococcus aceris and Gossyparia spuria. These species host Burkholderia bacteria that are localized in the cytoplasm of the fat body cells. Metagenome sequencing revealed very similar and highly reduced genomes (<900KBp) with a low GC content (~38%), making them the smallest and most AT-biased Burkholderia genomes yet sequenced. In their eroded genomes, both symbionts retain biosynthetic pathways for the essential amino acids leucine, isoleucine, valine, threonine, lysine, arginine, histidine, phenylalanine, and precursors for the semi-essential amino acid tyrosine, as well as the cobalamin-dependent methionine synthase MetH. A tryptophan biosynthesis pathway is conserved in the symbiont of G. spuria, but appeared pseudogenized in A. aceris, suggesting differential availability of tryptophan in the two host species' diets. In addition to the pathways for essential amino acid biosynthesis, both symbionts maintain biosynthetic pathways for multiple cofactors, including riboflavin, cobalamin, thiamine, and folate. The localization of Burkholderia symbionts and their genome traits indicate that the symbiosis between Burkholderia and eriococcids is younger than other hemipteran symbioses, but is functionally convergent. Our results add to the emerging picture of dynamic symbiont replacements in sap-sucking Hemiptera and highlight Burkholderia as widespread and versatile intra- and extracellular symbionts of animals, plants, and fungi.

The pan-genome of Mycobacterium avium subsp. paratuberculosis (Map) confirms ancestral lineage and reveals gene rearrangements within Map Type S

BACKGROUND

To date genomic studies on Map have concentrated on Type C strains with only a few Type S strains included for comparison. In this study the entire pan-genome of 261 Map genomes (205 Type C, 52 Type S and 4 Type B) and 7 Mycobacterium avium complex (Mac) genomes were analysed to identify genomic similarities and differences between the strains and provide more insight into the evolutionary relationship within this Mycobacterial species.

RESULTS

Our analysis of the core genome of all the Map isolates identified two distinct lineages, Type S and Type C Map that is consistent with previous phylogenetic studies of Map. Pan-genome analysis revealed that Map has a larger accessory genome than Mycobacterium avium subsp. avium (Maa) and Type C Map has a larger accessory genome than Type S Map. In addition, we found large rearrangements within Type S strains of Map and little to none in Type C and Type B strains. There were 50 core genes identified that were unique to Type S Map and there were no unique core genes identified between Type B and Type C Map strains. In Type C Map we identified an additional CE10 CAZyme class which was identified as an alpha/beta hydrolase and an additional polyketide and non-ribosomal peptide synthetase cluster. Consistent with previous analysis no plasmids and only incomplete prophages were identified in the genomes of Map. There were 45 hypothetical CRISPR elements identified with no associated cas genes.

CONCLUSION

This is the most comprehensive comparison of the genomic content of Map isolates to date and included the closing of eight Map genomes. The analysis revealed that there is greater variation in gene synteny within Type S strains when compared to Type C indicating that the Type C Map strain emerged after Type S. Further analysis of Type C and Type B genomes revealed that they are structurally similar with little to no genetic variation and that Type B Map may be a distinct clade within Type C Map and not a different strain type of Map. The evolutionary lineage of Maa and Map was confirmed as emerging after M. hominissuis.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Physiological and evolutionary contexts of a new symbiotic species from the nitrogen-recycling gut community of turtle ants

While genome sequencing has expanded our knowledge of symbiosis, role assignment within multi-species microbiomes remains challenging due to genomic redundancy and the uncertainties of in vivo impacts. We address such questions, here, for a specialized nitrogen (N) recycling microbiome of turtle ants, describing a new genus and species of gut symbiont-Ischyrobacter davidsoniae (Betaproteobacteria: Burkholderiales: Alcaligenaceae)-and its in vivo physiological context. A re-analysis of amplicon sequencing data, with precisely assigned Ischyrobacter reads, revealed a seemingly ubiquitous distribution across the turtle ant genus Cephalotes, suggesting ≥50 million years since domestication. Through new genome sequencing, we also show that divergent I. davidsoniae lineages are conserved in their uricolytic and urea-generating capacities. With phylogenetically refined definitions of Ischyrobacter and separately domesticated Burkholderiales symbionts, our FISH microscopy revealed a distinct niche for I. davidsoniae, with dense populations at the anterior ileum. Being positioned at the site of host N-waste delivery, in vivo metatranscriptomics and metabolomics further implicate I. davidsoniae within a symbiont-autonomous N-recycling pathway. While encoding much of this pathway, I. davidsoniae expressed only a subset of the requisite steps in mature adult workers, including the penultimate step deriving urea from allantoate. The remaining steps were expressed by other specialized gut symbionts. Collectively, this assemblage converts inosine, made from midgut symbionts, into urea and ammonia in the hindgut. With urea supporting host amino acid budgets and cuticle synthesis, and with the ancient nature of other active N-recyclers discovered here, I. davidsoniae emerges as a central player in a conserved and impactful, multipartite symbiosis.

Evolution of an Alternative Genetic Code in the Providencia Symbiont of the Hematophagous Leech Haementeria acuecueyetzin

Strict blood-feeding animals are confronted with a strong B-vitamin deficiency. Blood-feeding leeches from the Glossiphoniidae family, similarly to hematophagous insects, have evolved specialized organs called bacteriomes to harbor symbiotic bacteria. Leeches of the Haementeria genus have two pairs of globular bacteriomes attached to the esophagus which house intracellular "Candidatus Providencia siddallii" bacteria. Previous work analyzing a draft genome of the Providencia symbiont of the Mexican leech Haementeria officinalis showed that, in this species, the bacteria hold a reduced genome capable of synthesizing B vitamins. In this work, we aimed to expand our knowledge on the diversity and evolution of Providencia symbionts of Haementeria. For this purpose, we sequenced the symbiont genomes of three selected leech species. We found that all genomes are highly syntenic and have kept a stable genetic repertoire, mirroring ancient insect endosymbionts. Additionally, we found B-vitamin pathways to be conserved among these symbionts, pointing to a conserved symbiotic role. Lastly and most notably, we found that the symbiont of H. acuecueyetzin has evolved an alternative genetic code, affecting a portion of its proteome and showing evidence of a lineage-specific and likely intermediate stage of genetic code reassignment.

No Transcriptional Compensation for Extreme Gene Dosage Imbalance in Fragmented Bacterial Endosymbionts of Cicadas

Bacteria that form long-term intracellular associations with host cells lose many genes, a process that often results in tiny, gene-dense, and stable genomes. Paradoxically, the some of the same evolutionary processes that drive genome reduction and simplification may also cause genome expansion and complexification. A bacterial endosymbiont of cicadas, Hodgkinia cicadicola, exemplifies this paradox. In many cicada species, a single Hodgkinia lineage with a tiny, gene-dense genome has split into several interdependent cell and genome lineages. Each new Hodgkinia lineage encodes a unique subset of the ancestral unsplit genome in a complementary way, such that the collective gene contents of all lineages match the total found in the ancestral single genome. This splitting creates genetically distinct Hodgkinia cells that must function together to carry out basic cellular processes. It also creates a gene dosage problem where some genes are encoded by only a small fraction of cells while others are much more abundant. Here, by sequencing DNA and RNA of Hodgkinia from different cicada species with different amounts of splitting-along with its structurally stable, unsplit partner endosymbiont Sulcia muelleri-we show that Hodgkinia does not transcriptionally compensate to rescue the wildly unbalanced gene and genome ratios that result from lineage splitting. We also find that Hodgkinia has a reduced capacity for basic transcriptional control independent of the splitting process. Our findings reveal another layer of degeneration further pushing the limits of canonical molecular and cell biology in Hodgkinia and may partially explain its propensity to go extinct through symbiont replacement.

Role of gut symbionts of insect pests: A novel target for insect-pest control

Insects possess beneficial and nuisance values in the context of the agricultural sector and human life around them. An ensemble of gut symbionts assists insects to adapt to diverse and extreme environments and to occupy every available niche on earth. Microbial symbiosis helps host insects by supplementing necessary diet elements, providing protection from predators and parasitoids through camouflage, modulation of signaling pathway to attain homeostasis and to trigger immunity against pathogens, hijacking plant pathways to circumvent plant defence, acquiring the capability to degrade chemical pesticides, and degradation of harmful pesticides. Therefore, a microbial protection strategy can lead to overpopulation of insect pests, which can drastically reduce crop yield. Some studies have demonstrated increased insect mortality via the destruction of insect gut symbionts; through the use of antibiotics. The review summarizes various roles played by the gut microbiota of insect pests and some studies that have been conducted on pest control by targeting the symbionts. Manipulation or exploitation of the gut symbionts alters the growth and population of the host insects and is consequently a potential target for the development of better pest control strategies. Methods such as modulation of gut symbionts via CRISPR/Cas9, RNAi and the combining of IIT and SIT to increase the insect mortality are further discussed. In the ongoing insect pest management scenario, gut symbionts are proving to be the reliable, eco-friendly and novel approach in the integrated pest management.

Complex host/symbiont integration of a multi-partner symbiotic system in the eusocial aphid Ceratovacuna japonica

Some hemipteran insects rely on multiple endosymbionts for essential nutrients. However, the evolution of multi-partner symbiotic systems is not well-established. Here, we report a co-obligate symbiosis in the eusocial aphid, Ceratovacuna japonica. 16S rRNA amplicon sequencing unveiled co-infection with a novel Arsenophonus sp. symbiont and Buchnera aphidicola, a common obligate endosymbiont in aphids. Both symbionts were housed within distinct bacteriocytes and were maternally transmitted. The Buchnera and Arsenophonus symbionts had streamlined genomes of 432,286 bp and 853,149 bp, respectively, and exhibited metabolic complementarity in riboflavin and peptidoglycan synthesis pathways. These anatomical and genomic properties were similar to those of independently evolved multi-partner symbiotic systems, such as Buchnera-Serratia in Lachninae and Periphyllus aphids, representing remarkable parallelism. Furthermore, symbiont populations and bacteriome morphology differed between reproductive and soldier castes. Our study provides the first example of co-obligate symbiosis in Hormaphidinae and gives insight into the evolutionary genetics of this complex system.

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.

Fungal Strains with Identical Genomes Were Found at a Distance of 2000 Kilometers after 40 Years

Heredity and variation are inherent characteristics of species and are mainly reflected in the stability and variation of the genome; the former is relative, while the latter is continuous. However, whether life has both stable genomes and extremely diverse genomes at the same time is unknown. In this study, we isolated Sclerotinia sclerotiorum strains from sclerotium samples in Quincy, Washington State, USA, and found that four single-sclerotium-isolation strains (PB4, PB273, PB615, and PB623) had almost identical genomes to the reference strain 1980 isolated in the west of Nebraska 40 years ago. The genome of strain PB4 sequenced by the next-generation sequencing (NGS) and Pacific Biosciences (PacBio) sequencing carried only 135 single nucleotide polymorphisms (SNPs) and 18 structural variations (SVs) compared with the genome of strain 1980 and 48 SNPs were distributed on Contig_20. Based on data generated by NGS, three other strains, PB273, PB615, and PB623, had 256, 275, and 262 SNPs, respectively, against strain 1980, which were much less than in strain PB4 (532 SNPs) and none of them occurred on Contig_20, suggesting much closer genomes to strain 1980 than to strain PB4. All other strains from America and China are rich in SNPs with a range of 34,391-77,618 when compared with strain 1980. We also found that there were 39-79 SNPs between strain PB4 and its sexual offspring, 53.1% of which also occurred on Contig_20. Our discoveries show that there are two types of genomes in S. sclerotiorum, one is very stable and the other tends to change constantly. Investigating the mechanism of such genome stability will enhance our understanding of heredity and variation.

Convergent evolution of a labile nutritional symbiosis in ants

Ants are among the most successful organisms on Earth. It has been suggested that forming symbioses with nutrient-supplementing microbes may have contributed to their success, by allowing ants to invade otherwise inaccessible niches. However, it is unclear whether ants have evolved symbioses repeatedly to overcome the same nutrient limitations. Here, we address this question by comparing the independently evolved symbioses in Camponotus, Plagiolepis, Formica and Cardiocondyla ants. Our analysis reveals the only metabolic function consistently retained in all of the symbiont genomes is the capacity to synthesise tyrosine. We also show that in certain multi-queen lineages that have co-diversified with their symbiont for millions of years, only a fraction of queens carry the symbiont, suggesting ants differ in their colony-level reliance on symbiont-derived resources. Our results imply that symbioses can arise to solve common problems, but hosts may differ in their dependence on symbionts, highlighting the evolutionary forces influencing the persistence of long-term endosymbiotic mutualisms.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

The Di-Symbiotic Systems in the Aphids Sipha maydis and Periphyllus lyropictus Provide a Contrasting Picture of Recent Co-Obligate Nutritional Endosymbiosis in Aphids

Dependence on multiple nutritional bacterial symbionts forming a metabolic unit has repeatedly evolved in many insect species that feed on nutritionally unbalanced diets such as plant sap. This is the case for aphids of the subfamilies Lachninae and Chaitophorinae, which have evolved di-symbiotic systems in which the ancient obligate nutritional symbiont Buchnera aphidicola is metabolically complemented by an additional nutritional symbiont acquired more recently. Deciphering how different symbionts integrate both metabolically and anatomically in such systems is crucial to understanding how complex nutritional symbiotic systems function and evolve. In this study, we sequenced and analyzed the genomes of the symbionts B. aphidicola and Serratia symbiotica associated with the Chaitophorinae aphids Sipha maydis and Periphyllus lyropictus. Our results show that, in these two species, B. aphidicola and S. symbiotica complement each other metabolically (and their hosts) for the biosynthesis of essential amino acids and vitamins, but with distinct metabolic reactions supported by each symbiont depending on the host species. Furthermore, the S. symbiotica symbiont associated with S. maydis appears to be strictly compartmentalized into the specialized host cells housing symbionts in aphids, the bacteriocytes, whereas the S. symbiotica symbiont associated with P. lyropictus exhibits a highly invasive phenotype, presumably because it is capable of expressing a larger set of virulence factors, including a complete flagellum for bacterial motility. Such contrasting levels of metabolic and anatomical integration for two S. symbiotica symbionts that were recently acquired as nutritional co-obligate partners reflect distinct coevolutionary processes specific to each association.

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.

Aphid hologenomics: current status and future challenges

Aphids are important model organisms in ecological, developmental, and evolutionary studies of, for example, symbiosis, insect-plant interactions, pest management, and developmental polyphenism. Here, we review the recent progress made in the genomics of aphids and their symbionts: hologenomics. The reference genome of Acyrthosiphon pisum has been greatly improved, and chromosome-level assembly is now available. The genomes of over 20 aphid species have been sequenced, and comparative genomic analyses have revealed pervasive gene duplication and dynamic chromosomal rearrangements. Over 120 symbiont genomes (both obligate and facultative) have been sequenced, and modern deep-sequencing technologies have identified novel symbionts. The advances in hologenomics have helped to elucidate the dynamic evolution of facultative and co-obligate symbionts with the ancient obligate symbiont Buchnera.

Divergent paths in the evolutionary history of maternally transmitted clam symbionts

Vertical transmission of bacterial endosymbionts is accompanied by virtually irreversible gene loss that results in a progressive reduction in genome size. While the evolutionary processes of genome reduction have been well described in some terrestrial symbioses, they are less understood in marine systems where vertical transmission is rarely observed. The association between deep-sea vesicomyid clams and chemosynthetic Gammaproteobacteria is one example of maternally inherited symbioses in the ocean. Here, we assessed the contributions of drift, recombination and selection to genome evolution in two extant vesicomyid symbiont clades by comparing 15 representative symbiont genomes (1.017-1.586 Mb) to those of closely related bacteria and the hosts' mitochondria. Our analyses suggest that drift is a significant force driving genome evolution in vesicomyid symbionts, though selection and interspecific recombination appear to be critical for maintaining symbiont functional integrity and creating divergent patterns of gene conservation. Notably, the two symbiont clades possess putative functional differences in sulfide physiology, anaerobic respiration and dependency on environmental vitamin B12, which probably reflect adaptations to different ecological habitats available to each symbiont group. Overall, these results contribute to our understanding of the eco-evolutionary processes shaping reductive genome evolution in vertically transmitted symbioses.

Transitional genomes and nutritional role reversals identified for dual symbionts of adelgids (Aphidoidea: Adelgidae)

Many plant-sap-feeding insects have maintained a single, obligate, nutritional symbiont over the long history of their lineage. This senior symbiont may be joined by one or more junior symbionts that compensate for gaps in function incurred through genome-degradative forces. Adelgids are sap-sucking insects that feed solely on conifer trees and follow complex life cycles in which the diet fluctuates in nutrient levels. Adelgids are unusual in that both senior and junior symbionts appear to have been replaced repeatedly over their evolutionary history. Genomes can provide clues to understanding symbiont replacements, but only the dual symbionts of hemlock adelgids have been examined thus far. Here, we sequence and compare genomes of four additional dual-symbiont pairs in adelgids. We show that these symbionts are nutritional partners originating from diverse bacterial lineages and exhibiting wide variation in general genome characteristics. Although dual symbionts cooperate to produce nutrients, the balance of contributions varies widely across pairs, and total genome contents reflect a range of ages and degrees of degradation. Most symbionts appear to be in transitional states of genome reduction. Our findings support a hypothesis of periodic symbiont turnover driven by fluctuating selection for nutritional provisioning related to gains and losses of complex life cycles in their hosts.

Genomic Analyses of the Fungus Paraconiothyrium sp. Isolated from the Chinese White Wax Scale Insect Reveals Its Symbiotic Character

The Chinese white wax scale, Ericerus pela, is an insect native to China. It harbors a variety of microbes. The Paraconiothyrium fungus was isolated from E. pela and genome sequenced in this study. A fungal cytotoxicity assay was performed on the Aedes albopictus cell line C6/36. The assembled Paraconiothyrium sp. genome was 39.55 Mb and consisted of 14,174 genes. The coding sequences accounted for 50.75% of the entire genome. Functional pathway analyses showed that Paraconiothyrium sp. possesses complete pathways for the biosynthesis of 20 amino acids, 10 of which E. pela lacks. It also had complementary genes in the vitamin B groups synthesis pathways. Secondary metabolism prediction showed many gene clusters that produce polyketide. Additionally, a large number of genes associated with ‘reduced virulence’ in the genome were annotated with the Pathogen–Host Interaction database. A total of 651 genes encoding carbohydrate-active enzymes were predicted to be mostly involved in plant polysaccharide degradation. Pan-specific genomic analyses showed that genes unique to Paraconiothyrium sp. were enriched in the pathways related to amino acid metabolism and secondary metabolism. GO annotation analysis yielded similar results. The top COG categories were ‘carbohydrate transport and metabolism’, ‘lipid transport and metabolism’, and ‘secondary metabolite biosynthesis, transport and catabolism’. Phylogenetic analyses based on gene family and pan genes showed that Paraconiothyrium sp is clustered together with species from the Didymosphaeriaceae family. A multi-locus sequence analysis showed that it converged with the same branch as P. brasiliense and they formed one group with fungi from the Paraconiothyrium genus. To validate the in vitro toxicity of Paraconiothyrium sp., a cytotoxicity assay was performed. The results showed that medium-cultured Paraconiothyrium sp. had no harmful effect on cell viability. No toxins were secreted by the fungus during growth. Our results imply that Paraconiothyrium sp. may establish a symbiotic relationship with the host to supply complementary nutrition to E. pela.

Adaptation to genome decay in the structure of the smallest eukaryotic ribosome

The evolution of microbial parasites involves the counterplay between natural selection forcing parasites to improve and genetic drifts forcing parasites to lose genes and accumulate deleterious mutations. Here, to understand how this counterplay occurs at the scale of individual macromolecules, we describe cryo-EM structure of ribosomes from Encephalitozoon cuniculi, a eukaryote with one of the smallest genomes in nature. The extreme rRNA reduction in E. cuniculi ribosomes is accompanied with unparalleled structural changes, such as the evolution of previously unknown molten rRNA linkers and bulgeless rRNA. Furthermore, E. cuniculi ribosomes withstand the loss of rRNA and protein segments by evolving an ability to use small molecules as structural mimics of degenerated rRNA and protein segments. Overall, we show that the molecular structures long viewed as reduced, degenerated, and suffering from debilitating mutations possess an array of compensatory mechanisms that allow them to remain active despite the extreme molecular reduction.

Many parasitic organisms contain molecular structures that are drastically smaller than analogous structures in non-parasitic organisms. Here the authors describe a cryo-EM structure of the ribosome from E. cuniculi that reveals that it compensated rRNA truncations by evolving the ability to use small molecules as ribosomal building blocks.

Prokaryote pangenomes are dynamic entities

Prokaryote pangenomes are influenced heavily by environmental factors and the opportunity for gene gain and loss events. As the field of pangenome analysis has expanded, so has the need to fully understand the complexity of how eco-evolutionary dynamics shape pangenomes. Here, we describe current models of pangenome evolution and discuss their suitability and accuracy. We suggest that pangenomes are dynamic entities under constant flux, highlighting the influence of two-way interactions between pangenome and environment. New classifications of core and accessory genes are also considered, underscoring the need for continuous evaluation of nomenclature in a fast-moving field. We conclude that future models of pangenome evolution should incorporate eco-evolutionary dynamics to fully encompass their dynamic, changeable nature.

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic strategies.

Microbe Profile: Buchnera aphidicola: ancient aphid accomplice and endosymbiont exemplar

Buchnera aphidicola is an obligate endosymbiont of aphids that cannot be cultured outside of hosts. It exists as diverse strains in different aphid species, and phylogenetic reconstructions show that it has been maternally transmitted in aphids for >100 million years. B. aphidicola genomes are highly reduced and show conserved gene order and no gene acquisition, but encoded proteins undergo rapid evolution. Aphids depend on B. aphidicola for biosynthesis of essential amino acids and as an integral part of embryonic development. How B. aphidicola populations are regulated within hosts remains little known.

Potato leafroll virus reduces Buchnera aphidocola titer and alters vector transcriptome responses

Viruses in the Luteoviridae family, such as Potato leafroll virus (PLRV), are transmitted by aphids in a circulative and nonpropagative mode. This means the virions enter the aphid body through the gut when they feed from infected plants and then the virions circulate through the hemolymph to enter the salivary glands before being released into the saliva. Although these viruses do not replicate in their insect vectors, previous studies have demonstrated viruliferous aphid behavior is altered and the obligate symbiont of aphids, Buchnera aphidocola, may be involved in transmission. Here we provide the transcriptome of green peach aphids (Myzus persicae) carrying PLRV and virus-free control aphids using Illumina sequencing. Over 150 million paired-end reads were obtained through Illumina sequencing, with an average of 19 million reads per library. The comparative analysis identified 134 differentially expressed genes (DEGs) between the M. persicae transcriptomes, including 64 and 70 genes that were up- and down-regulated in aphids carrying PLRV, respectively. Using functional classification in the GO databases, 80 of the DEGs were assigned to 391 functional subcategories at category level 2. The most highly up-regulated genes in aphids carrying PLRV were cytochrome p450s, genes related to cuticle production, and genes related to development, while genes related to heat shock proteins, histones, and histone modification were the most down-regulated. PLRV aphids had reduced Buchnera titer and lower abundance of several Buchnera transcripts related to stress responses and metabolism. These results suggest carrying PLRV may reduce both aphid and Buchnera genes in response to stress. This work provides valuable basis for further investigation into the complicated mechanisms of circulative and nonpropagative transmission.

Comparative Genomics Underlines Multiple Roles of Profftella, an Obligate Symbiont of Psyllids: Providing Toxins, Vitamins, and Carotenoids

The Asian citrus psyllid Diaphorina citri (Insecta: Hemiptera: Psylloidea), a serious pest of citrus species worldwide, harbors vertically transmitted intracellular mutualists, Candidatus Profftella armatura (Profftella_DC, Gammaproteobacteria: Burkholderiales) and Candidatus Carsonella ruddii (Carsonella_DC, Gammaproteobacteria: Oceanospirillales). Whereas Carsonella_DC is a typical nutritional symbiont, Profftella_DC is a unique defensive symbiont with organelle-like features, including intracellular localization within the host, perfect infection in host populations, vertical transmission over evolutionary time, and drastic genome reduction down to much less than 1 Mb. Large parts of the 460-kb genome of Profftella_DC are devoted to genes for synthesizing a polyketide toxin; diaphorin. To better understand the evolution of this unusual symbiont, the present study analyzed the genome of Profftella_Dco, a sister lineage to Profftella_DC, using Diaphorina cf. continua, a host psyllid congeneric with D. citri. The genome of coresiding Carsonella (Carsonella_Dco) was also analyzed. The analysis revealed nearly perfect synteny conservation in these genomes with their counterparts from D. citri. The substitution rate analysis further demonstrated genomic stability of Profftella which is comparable to that of Carsonella. Profftella_Dco and Profftella_DC shared all genes for the biosynthesis of diaphorin, hemolysin, riboflavin, biotin, and carotenoids, underlining multiple roles of Profftella, which may contribute to stabilizing symbiotic relationships with the host. However, acyl carrier proteins were extensively amplified in polyketide synthases DipP and DipT for diaphorin synthesis in Profftella_Dco. This level of acyl carrier protein augmentation, unprecedented in modular polyketide synthases of any known organism, is not thought to influence the polyketide structure but may improve the synthesis efficiency.

Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids

Numerous animal lineages have maternally inherited symbionts that are required for host reproduction and growth. Endosymbionts also pose a risk to their hosts because of the mutational decay of their genomes through genetic drift or to selfish mutations that favor symbiont fitness over host fitness. One model for heritable endosymbiosis is the association of aphids with their obligate bacterial symbiont, Buchnera We experimentally established heteroplasmic pea aphid matrilines containing pairs of closely related Buchnera haplotypes and used deep sequencing of diagnostic markers to measure haplotype frequencies in successive host generations. These frequencies were used to estimate the effective population size of Buchnera within hosts (i.e., the transmission bottleneck size) and the extent of within-host selection. The within-host effective population size was in the range of 10 to 20, indicating a strong potential for genetic drift and fixation of deleterious mutations. Remarkably, closely related haplotypes were subject to strong within-host selection, with selection coefficients as high as 0.5 per aphid generation. In one case, the direction of selection depended on the thermal environment and went in the same direction as between-host selection. In another, a new mutant haplotype had a strong within-host advantage under both environments but had no discernible effect on host-level fitness under laboratory conditions. Thus, within-host selection can be strong, resulting in a rapid fixation of mutations with little impact on host-level fitness. Together, these results show that within-host selection can drive evolution of an obligate symbiont, accelerating sequence evolution.

Adaptive Variation of Buchnera Endosymbiont Density in Aphid Host Acyrthosiphon pisum Controlled by Environmental Conditions

The scarcity of transcriptional regulatory genes in Buchnera aphidicola, an obligate endosymbiont in aphids, suggests the stability of expressed gene patterns and metabolic pathways. This observation argues in favor of the hypothesis that this endosymbiont bacteria might contribute little to the host adaptation when aphid hosts are facing challenging fluctuating environment. Finding evidence for the increased expression or silenced genes involved in metabolic pathways under the pressure of stress conditions and/or a given environment has been challenging for experimenters with this bacterial symbiotic model. Transcriptomic data have shown that Buchnera gene expression changes are confined to a narrow range when the aphids face brutal environmental variations. In this report, we demonstrate that instead of manipulating individual genes, the conditions may act on the relative mass of endosymbiont corresponding to the needs of the host. The control of the fluctuating number of endosymbiont cells per individual host appears to be an unexpected regulatory modality that contributes to the adaptation of aphids to their environment. This feature may account for the success of the symbiotic advantages in overcoming the drastic changes in temperature and food supplies during evolution.

Comparative Genomics Provides Insights into the Genetic Diversity and Evolution of the DPANN Superphylum

DPANN is known as highly diverse, globally widespread, and mostly ectosymbiotic archaeal superphylum. However, this group of archaea was overlooked for a long time, and there were limited in-depth studies reported. In this investigation, 41 metagenome-assembled genomes (MAGs) belonging to the DPANN superphylum were recovered (18 MAGs had average nucleotide identity [ANI] values of <95% and a percentage of conserved proteins [POCP] of >50%, while 14 MAGs showed a POCP of <50%), which were analyzed comparatively with 515 other published DPANN genomes. Mismatches to known 16S rRNA gene primers were identified among 16S rRNA genes of DPANN archaea. Numbers of gene families lost (mostly related to energy and amino acid metabolism) were over three times greater than those gained in the evolution of DPANN archaea. Lateral gene transfer (LGT; ∼45.5% was cross-domain) had facilitated niche adaption of the DPANN archaea, ensuring a delicate equilibrium of streamlined genomes with efficient niche-adaptive strategies. For instance, LGT-derived cytochrome bd ubiquinol oxidase and arginine deiminase in the genomes of “Candidatus Micrarchaeota” could help them better adapt to aerobic acidic mine drainage habitats. In addition, most DPANN archaea acquired enzymes for biosynthesis of extracellular polymeric substances (EPS) and transketolase/transaldolase for the pentose phosphate pathway from Bacteria.

IMPORTANCE The domain Archaea is a key research model for gaining insights into the origin and evolution of life, as well as the relevant biogeochemical processes. The discovery of nanosized DPANN archaea has overthrown many aspects of microbiology. However, the DPANN superphylum still contains a vast genetic novelty and diversity that need to be explored. Comprehensively comparative genomic analysis on the DPANN superphylum was performed in this study, with an attempt to illuminate its metabolic potential, ecological distribution and evolutionary history. Many interphylum differences within the DPANN superphylum were found. For example, Altiarchaeota had the biggest genome among DPANN phyla, possessing many pathways missing in other phyla, such as formaldehyde assimilation and the Wood-Ljungdahl pathway. In addition, LGT acted as an important force to provide DPANN archaeal genetic flexibility that permitted the occupation of diverse niches. This study has advanced our understanding of the diversity and genome evolution of archaea.

Isolation of the Buchnera aphidicola flagellum basal body complexes from the Buchnera membrane

Buchnera aphidicola is an intracellular bacterial symbiont of aphids and maintains a small genome of only 600 kbps. Buchnera is thought to maintain only genes relevant to the symbiosis with its aphid host. Curiously, the Buchnera genome contains gene clusters coding for flagellum basal body structural proteins and for flagellum type III export machinery. These structures have been shown to be highly expressed and present in large numbers on Buchnera cells. No recognizable pathogenicity factors or secreted proteins have been identified in the Buchnera genome, and the relevance of this protein complex to the symbiosis is unknown. Here, we show isolation of Buchnera flagellum basal body proteins from the cellular membrane of Buchnera, confirming the enrichment of flagellum basal body proteins relative to other proteins in the Buchnera proteome. This will facilitate studies of the structure and function of the Buchnera flagellum structure, and its role in this model symbiosis.

The Origin of Niches and Species in the Bacterial World

Niches are spaces for the biological units of selection, from cells to complex communities. In a broad sense, "species" are biological units of individuation. Niches do not exist without individual organisms, and every organism has a niche. We use "niche" in the Hutchinsonian sense as an abstraction of a multidimensional environmental space characterized by a variety of conditions, both biotic and abiotic, whose quantitative ranges determine the positive or negative growth rates of the microbial individual, typically a species, but also parts of the communities of species contained in this space. Microbial organisms ("species") constantly diversify, and such diversification (radiation) depends on the possibility of opening up unexploited or insufficiently exploited niches. Niche exploitation frequently implies "niche construction," as the colonized niche evolves with time, giving rise to new potential subniches, thereby influencing the selection of a series of new variants in the progeny. The evolution of niches and organisms is the result of reciprocal interacting processes that form a single unified process. Centrifugal microbial diversification expands the limits of the species' niches while a centripetal or cohesive process occurs simultaneously, mediated by horizontal gene transfers and recombinatorial events, condensing all of the information recovered during the diversifying specialization into "novel organisms" (possible future species), thereby creating a more complex niche, where the selfishness of the new organism(s) establishes a "homeostatic power" limiting the niche's variation. Once the niche's full carrying capacity has been reached, reproductive isolation occurs, as no foreign organisms can outcompete the established population/community, thereby facilitating speciation. In the case of individualization-speciation of the microbiota, its contribution to the animal' gut structure is a type of "niche construction," the result of crosstalk between the niche (host) and microorganism(s). Lastly, there is a parallelism between the hierarchy of niches and that of microbial individuals. The increasing anthropogenic effects on the biosphere (such as globalization) might reduce the diversity of niches and bacterial individuals, with the potential emergence of highly transmissible multispecialists (which are eventually deleterious) resulting from the homogenization of the microbiosphere, a possibility that should be explored and prevented.

The Metabolic Building Blocks of a Minimal Cell

Defining the essential gene components for a system to be considered alive is a crucial step toward the synthesis of artificial life. Fifteen years ago, Gil and coworkers proposed the core of a putative minimal bacterial genome, which would provide the capability to achieve metabolic homeostasis, reproduce, and evolve to a bacterium in an ideally controlled environment. They also proposed a simplified metabolic chart capable of providing energy and basic components for a minimal living cell. For this work, we have identified the components of the minimal metabolic network based on the aforementioned studies, associated them to the KEGG database and, by applying the MetaDAG methodology, determined its Metabolic Building Blocks (MBB) and reconstructed its metabolic Directed Acyclic Graph (m-DAG). The reaction graph of this metabolic network consists of 80 compounds and 98 reactions, while its m-DAG has 36 MBBs. Additionally, we identified 12 essential reactions in the m-DAG that are critical for maintaining the connectivity of this network. In a similar manner, we reconstructed the m-DAG of JCVI-syn3.0, which is an artificially designed and manufactured viable cell whose genome arose by minimizing the one from Mycoplasma mycoides JCVI-syn1.0, and of "Candidatus Nasuia deltocephalinicola", the bacteria with the smallest natural genome known to date. The comparison of the m-DAGs derived from a theoretical, an artificial, and a natural genome denote slightly different lifestyles, with a consistent core metabolism. The MetaDAG methodology we employ uses homogeneous descriptors and identifiers from the KEGG database, so that comparisons between bacterial strains are not only easy but also suitable for many research fields. The modeling of m-DAGs based on minimal metabolisms can be the first step for the synthesis and manipulation of minimal cells.

The Life of an Insect Endosymbiont from the Cradle to the Grave

Host-beneficial endosymbioses, which are formed when a microorganism takes up residence inside another cell and provides a fitness advantage to the host, have had a dramatic influence on the evolution of life. These intimate relationships have yielded the mitochondrion and the plastid (chloroplast) - the ancient organelles that in part define eukaryotic life - along with many more recent associations involving a wide variety of hosts and microbial partners. These relationships are often envisioned as stable associations that appear cooperative and persist for extremely long periods of time. But recent evidence suggests that this stable state is often born from turbulent and conflicting origins, and that the apparent stability of many beneficial endosymbiotic relationships - although certainly real in many cases - is not an inevitable outcome of these associations. Here we review how stable endosymbioses form, how they are maintained, and how they sometimes break down and are reborn. We focus on relationships formed by insects and their resident microorganisms because these symbioses have been the focus of significant empirical work over the last two decades. We review these relationships over five life stages: origin, birth, middle age, old age, and death.

Bacterial symbionts support larval sap feeding and adult folivory in (semi-)aquatic reed beetles

Symbiotic microbes can enable their host to access untapped nutritional resources but may also constrain niche space by promoting specialization. Here, we reconstruct functional changes in the evolutionary history of the symbiosis between a group of (semi-)aquatic herbivorous insects and mutualistic bacteria. Sequencing the symbiont genomes across 26 species of reed beetles (Chrysomelidae, Donaciinae) spanning four genera indicates that the genome-eroded mutualists provide life stage-specific benefits to larvae and adults, respectively. In the plant sap-feeding larvae, the symbionts are inferred to synthesize most of the essential amino acids as well as the B vitamin riboflavin. The adult reed beetles’ folivory is likely supported by symbiont-encoded pectinases that complement the host-encoded set of cellulases, as revealed by transcriptome sequencing. However, mapping the occurrence of the symbionts’ pectinase genes and the hosts’ food plant preferences onto the beetles’ phylogeny reveals multiple independent losses of pectinase genes in lineages that switched to feeding on pectin-poor plants, presumably constraining their hosts’ subsequent adaptive potential.

Symbiotic microbes in insects can enable their hosts to access untapped nutritional resources. Here, the authors show that symbiotic bacteria in reed beetles can provide essential amino acids to sap-feeding larvae and help leaf-feeding adults to degrade pectin, respectively.

Identification of gene fusion events in Mycobacterium tuberculosis that encode chimeric proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen responsible for causing tuberculosis. The harsh environment in which M. tuberculosis survives requires this pathogen to continuously adapt in order to maintain an evolutionary advantage. However, the apparent absence of horizontal gene transfer in M. tuberculosis imposes restrictions in the ways by which evolution can occur. Large-scale changes in the genome can be introduced through genome reduction, recombination events and structural variation. Here, we identify a functional chimeric protein in the ppe38-71 locus, the absence of which is known to have an impact on protein secretion and virulence. To examine whether this approach was used more often by this pathogen, we further develop software that detects potential gene fusion events from multigene deletions using whole genome sequencing data. With this software we could identify a number of other putative gene fusion events within the genomes of M. tuberculosis isolates. We were able to demonstrate the expression of one of these gene fusions at the protein level using mass spectrometry. Therefore, gene fusions may provide an additional means of evolution for M. tuberculosis in its natural environment whereby novel chimeric proteins and functions can arise.

Host Plants Influence the Symbiont Diversity of Eriosomatinae (Hemiptera: Aphididae)

Eriosomatinae is a particular aphid group with typically heteroecious holocyclic life cycle, exhibiting strong primary host plant specialization and inducing galls on primary host plants. Aphids are frequently associated with bacterial symbionts, which can play fundamental roles in the ecology and evolution of their host aphids. However, the bacterial communities in Eriosomatinae are poorly known. In the present study, using high-throughput sequencing of the bacterial 16S ribosomal RNA gene, we surveyed the bacterial flora of eriosomatines and explored the associations between symbiont diversity and aphid relatedness, aphid host plant and geographical distribution. The microbiota of Eriosomatinae is dominated by the heritable primary endosymbiont Buchnera and several facultative symbionts. The primary endosymbiont Buchnera is expectedly the most abundant symbiont across all species. Six facultative symbionts were identified. Regiella was the most commonly identified facultative symbiont, and multiple infections of facultative symbionts were detected in the majority of the samples. Ordination analyses and statistical tests show that the symbiont community of aphids feeding on plants from the family Ulmaceae were distinguishable from aphids feeding on other host plants. Species in Eriosomatinae feeding on different plants are likely to carry different symbiont compositions. The symbiont distributions seem to be not related to taxonomic distance and geographical distance. Our findings suggest that host plants can affect symbiont maintenance, and will improve our understanding of the interactions between aphids, their symbionts and ecological conditions.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

What Goes Up Might Come Down: the Spectacular Spread of an Endosymbiont Is Followed by Its Decline a Decade Later

Facultative, intracellular bacterial symbionts of arthropods may dramatically affect host biology and reproduction. The length of these symbiont-host associations may be thousands to millions of years, and while symbiont loss is predicted, there have been very few observations of a decline of symbiont infection rates. In a population of the sweet potato whitefly species (Bemisia tabaci MEAM1) in Arizona, USA, we documented the frequency decline of a strain of Rickettsia in the Rickettsia bellii clade from near-fixation in 2011 to 36% of whiteflies infected in 2017. In previous studies, Rickettsia had been shown to increase from 1 to 97% from 2000 to 2006 and remained at high frequency for at least five years. At that time, Rickettsia infection was associated with both fitness benefits and female bias. In the current study, we established matrilines of whiteflies from the field (2016, Rickettsia infection frequency = 58%) and studied (a) Rickettsia vertical transmission, (b) fitness and sex ratios associated with Rickettsia infection, (c) symbiont titer, and (d) bacterial communities within whiteflies. The vertical transmission rate was high, approximately 98%. Rickettsia infection in the matrilines was not associated with fitness benefits or sex ratio bias and appeared to be slightly costly, as more Rickettsia-infected individuals produced non-hatching eggs. Overall, the titer of Rickettsia in the matrilines was lower in 2016 than in the whiteflies collected in 2011, but the titer distribution appeared bimodal, with high- and low-titer lines, and constancy of the average titer within lines over three generations. We found neither association between Rickettsia titer and fitness benefits or sex ratio bias nor evidence that Rickettsia was replaced by another secondary symbiont. The change in the interaction between symbiont and host in 2016 whiteflies may explain the drop in symbiont frequency we observed.

Genome Evolution of the Obligate Endosymbiont Buchnera aphidicola

An evolutionary consequence of uniparentally transmitted symbiosis is degradation of symbiont genomes. We use the system of aphids and their maternally inherited obligate endosymbiont, Buchnera aphidicola, to explore the evolutionary process of genome degradation. We compared complete genome sequences for 39 Buchnera strains, including 23 newly sequenced symbiont genomes from diverse aphid hosts. We reconstructed the genome of the most recent shared Buchnera ancestor, which contained 616 protein-coding genes, and 39 RNA genes. The extent of subsequent gene loss varied across lineages, resulting in modern genomes ranging from 412 to 646 kb and containing 354-587 protein-coding genes. Loss events were highly nonrandom across loci. Genes involved in replication, transcription, translation, and amino acid biosynthesis are largely retained, whereas genes underlying ornithine biosynthesis, stress responses, and transcriptional regulation were lost repeatedly. Aside from losses, gene order is almost completely stable. The main exceptions involve movement between plasmid and chromosome locations of genes underlying tryptophan and leucine biosynthesis and supporting nutrition of aphid hosts. This set of complete genomes enabled tests for signatures of positive diversifying selection. Of 371 Buchnera genes tested, 29 genes show strong support for ongoing positive selection. These include genes encoding outer membrane porins that are expected to be involved in direct interactions with hosts. Collectively, these results indicate that extensive genome reduction occurred in the ancestral Buchnera prior to aphid diversification and that reduction has continued since, with losses greater in some lineages and for some loci.

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.

Diverse deep-sea anglerfishes share a genetically reduced luminous symbiont that is acquired from the environment

Deep-sea anglerfishes are relatively abundant and diverse, but their luminescent bacterial symbionts remain enigmatic. The genomes of two symbiont species have qualities common to vertically transmitted, host-dependent bacteria. However, a number of traits suggest that these symbionts may be environmentally acquired. To determine how anglerfish symbionts are transmitted, we analyzed bacteria-host codivergence across six diverse anglerfish genera. Most of the anglerfish species surveyed shared a common species of symbiont. Only one other symbiont species was found, which had a specific relationship with one anglerfish species, Cryptopsaras couesii. Host and symbiont phylogenies lacked congruence, and there was no statistical support for codivergence broadly. We also recovered symbiont-specific gene sequences from water collected near hosts, suggesting environmental persistence of symbionts. Based on these results we conclude that diverse anglerfishes share symbionts that are acquired from the environment, and that these bacteria have undergone extreme genome reduction although they are not vertically transmitted.

The deep sea is home to many different species of anglerfish, a group of animals in which females often display a dangling lure on the top of their heads. This organ shelters bacteria that make light, a partnership (known as symbiosis) that benefits both parties. The bacteria get a safe environment in which to grow, while the animal may use the light to confuse predators as well as attract prey and mates.

The genetic information of these bacteria has changed since they became associated with their host. Their genomes have become smaller and more specialized, limiting their ability to survive outside of the fish. This phenomenon is also observed in other symbiotic bacteria, but mostly in microorganisms that are directly transmitted from parent to offspring, never having to live on their own. Yet, some evidence suggests that the bacteria in the lure of anglerfish may be spending time in the water until they find a new host, crossing thousands of meters of ocean in the process.

To explore this paradox, Baker et al. looked into the type of bacteria carried by different groups of anglerfish. If each type of fish has its own kind of bacteria, this would suggest that the microorganisms are passed from one generation to the next, and are evolving with their hosts. On the other hand, if the same sort of bacteria can be found in different anglerfish species, this would imply that the bacteria pass from host to host and evolve independently from the fish.

Genetic data analysis showed that amongst six groups of anglerfishes, one species of bacteria is shared across five groups while another is specific to one type of fish. The analyses also revealed that anglerfish and their bacteria are most likely not evolving together. This means that the bacteria must make the difficult journey from host to host by persisting in the deep sea, which was confirmed by finding the genetic information of these bacteria in the water near the fish.

Anglerfish and the bacteria that light up their lure are hard to study, as they live so deep in the ocean. In fact, many symbiotic relationships are equally difficult to investigate. Examining genetic information can help to give an insight into how hosts and bacteria interact across the tree of life.

Conservation of transcriptional elements in the obligate symbiont of the whitefly Bemisia tabaci

Background

Bacterial symbiosis is widespread in arthropods, especially in insects. Some of the symbionts undergo a long-term co-evolution with the host, resulting in massive genome decay. One particular consequence of genome decay is thought to be the elimination of transcriptional elements within both the coding region and intergenic sequences. In the whitefly Bemisia tabaci species complex, the obligate symbiont Candidatus Portiera aleyrodidarum is of vital importance in nutrient provision, and yet little is known about the regulatory capacities of it.

Methods

Portiera genomes of two whitefly species in China were sequenced and assembled. Gene content of these two Portiera genomes was predicted, and then subjected to Kyoto Encyclopedia of Genes and Genomes pathway analysis. Together with two other Portiera genomes from whitefly species available previously, four Portiera genomes were utilized to investigate regulatory capacities of Portiera, focusing on transcriptional elements, including genes related with transcription and functional elements within the intergenic spacers.

Results

Comparative analyses of the four Portiera genomes of whitefly B. tabaci indicate that the obligate symbionts Portiera is similar in different species of whiteflies, in terms of general genome features and possible functions in the biosynthesis of essential amino acids. The screening of transcriptional factors suggests compromised ability of Portiera to regulate the essential amino acid biosynthesis pathways. Meanwhile, thermal tolerance ability of Portiera is indicated with the detection of a σ³² factor, as well as two predicted σ³² binding sites. Within intergenic spacers, functional elements are predicted, including 37 Shine-Dalgarno sequences and 34 putative small RNAs.

Cladogenesis and Genomic Streamlining in Extracellular Endosymbionts of Tropical Stink Bugs

Phytophagous stink bugs are globally distributed and many harbor vertically inherited bacterial symbionts that are extracellular, yet little is known about how the symbiont’s genomes have evolved under this transmission strategy. Genome reduction is common in insect intracellular symbionts but limited genome sampling of the extracellular symbionts of distantly related stink bugs has precluded inferring patterns of extracellular symbiont genome evolution. To address this knowledge gap, we completely sequenced the genomes of the uncultivable bacterial symbionts of four neotropical stink bugs of the Edessa genus. Phylogenetic and comparative analyses indicated that the symbionts form a clade within the Pantoea genus and their genomes are highly reduced (∼0.8 Mb). Furthermore, genome synteny analysis and a jackknife approach for phylogenetic reconstruction, which corrected for long branch attraction artifacts, indicated that the Edessa symbionts were the result of a single symbiotic event that was distinct from the symbiosis event giving rise to Candidatus “Pantoea carbekii,” the extracellular symbiont of the invasive pentatomid stink bug, Halyomorpha halys. Metabolic functions inferred from the Edessa symbiont genomes suggests a shift in genomic composition characteristic of its lifestyle in that they retained many host-supportive functions while undergoing dramatic gene loss and establishing a stable relationship with their host insects. Given the undersampled nature of extracellular insect symbionts, this study is the first comparative analysis of these symbiont genomes from four distinct Edessa stink bug species. Finally, we propose the candidate name “Candidatus Pantoea edessiphila” for the species of these symbionts with strain designations according to their host species.

Phylogeny, Divergent Evolution, and Speciation of Sulfur-Oxidizing Acidithiobacillus Populations

Background

Habitats colonized by acidophiles as an ideal physical barrier may induce genetic exchange of microbial members within the common communities, but little is known about how species in extremely acidic environments diverge and evolve.

Results

Using the acidophilic sulfur-oxidizer Acidithiobacillus as a case study, taxonomic reclassifications of many isolates provides novel insights into their phylogenetic lineage. Whole-genome-based comparisons were attempted to investigate the intra- and inter-species divergence. Recent studies clarified that functional and structural specificities of bacterial strains might provide opportunities for adaptive evolution responding to local environmental conditions. Acidophilic microorganisms play a key role in the acidification of natural waters and thus the formation of extremely acidic environments, and the feedbacks of the latter might confer the distinct evolutionary patterns of Acidithiobacillus spp. Varied horizontal gene transfer events occurred in different bacterial strains, probably resulting in the expansion of Acidithiobacillus genomes. Gene loss as another evolutionary force might cause the adaptive phenotypic diversity. A conceptual model for potential community-dependent evolutionary adaptation was thus proposed to illustrate the observed genome differentiation.

Conclusions

Collectively, the findings shed light on the phylogeny and divergent evolution of Acidithiobacillus strains, and provided a useful reference for evolutionary studies of other extremophiles.

Electronic supplementary material

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