Gene size reduction in the bacterial aphid endosymbiont, Buchnera

Comparative genomics of Mollicutes-related endobacteria supports a late invasion into Mucoromycota fungi

Diverse members of early-diverging Mucoromycota, including mycorrhizal taxa and soil-associated Mortierellaceae, are known to harbor Mollicutes-related endobacteria (MRE). It has been hypothesized that MRE were acquired by a common ancestor and transmitted vertically. Alternatively, MRE endosymbionts could have invaded after the divergence of Mucoromycota lineages and subsequently spread to new hosts horizontally. To better understand the evolutionary history of MRE symbionts, we generated and analyzed four complete MRE genomes from two Mortierellaceae genera: Linnemannia (MRE-L) and Benniella (MRE-B). These genomes include the smallest known of fungal endosymbionts and showed signals of a tight relationship with hosts including a reduced functional capacity and genes transferred from fungal hosts to MRE. Phylogenetic reconstruction including nine MRE from mycorrhizal fungi revealed that MRE-B genomes are more closely related to MRE from Glomeromycotina than MRE-L from the same host family. We posit that reductions in genome size, GC content, pseudogene content, and repeat content in MRE-L may reflect a longer-term relationship with their fungal hosts. These data indicate Linnemannia and Benniella MRE were likely acquired independently after their fungal hosts diverged from a common ancestor. This work expands upon foundational knowledge on minimal genomes and provides insights into the evolution of bacterial endosymbionts.

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.

Reduced Genome of the Gut Symbiotic Bacterium "Candidatus Benitsuchiphilus tojoi" Provides Insight Into Its Possible Roles in Ecology and Adaptation of the Host Insect

Diverse animals, including insects, harbor microbial symbionts within their gut, body cavity, or cells. The subsocial parastrachiid stinkbug Parastrachia japonensis is well-known for its peculiar ecological and behavioral traits, including its prolonged non-feeding diapause period and maternal care of eggs/nymphs in an underground nest. P. japonensis harbors a specific bacterial symbiont within the gut cavity extracellularly, which is vertically inherited through maternal excretion of symbiont-containing white mucus. Thus far, biological roles of the symbiont in the host lifecycle has been little understood. Here we sequenced the genome of the uncultivable gut symbiont “Candidatus Benitsuchiphilus tojoi.” The symbiont has an 804 kb circular chromosome encoding 606 proteins and a 14.5 kb plasmid encoding 13 proteins. Phylogenetic analysis indicated that the bacterium is closely related to other obligate insect symbionts belonging to the Gammaproteobacteria, including Buchnera of aphids and Blochmannia of ants, and the most closely related to Ishikawaella, an extracellular gut symbiont of plataspid stinkbugs. These data suggested that the symbiont genome has evolved like highly reduced gamma-proteobacterial symbiont genomes reported from a variety of insects. The presence of genes involved in biosynthesis pathways for amino acids, vitamins, and cofactors in the genome implicated the symbiont as a nutritional mutualist, supplementing essential nutrients to the host. Interestingly, the symbiont’s plasmid encoded genes for thiamine and carotenoid synthesis pathways, suggesting the possibility of additional functions of the symbiont for protecting the host against oxidative stress and DNA damage. Finally, possible involvement of the symbiont in uric acid metabolism during diapause is discussed.

Evolutionary Dynamics of Chloroplast Genomes in Low Light: A Case Study of the Endolithic Green Alga Ostreobium quekettii

Some photosynthetic organisms live in extremely low light environments. Light limitation is associated with selective forces as well as reduced exposure to mutagens, and over evolutionary timescales it can leave a footprint on species’ genomes. Here, we present the chloroplast genomes of four green algae (Bryopsidales, Ulvophyceae), including the endolithic (limestone-boring) alga Ostreobium quekettii, which is a low light specialist. We use phylogenetic models and comparative genomic tools to investigate whether the chloroplast genome of Ostreobium corresponds to our expectations of how low light would affect genome evolution. Ostreobium has the smallest and most gene-dense chloroplast genome among Ulvophyceae reported to date, matching our expectation that light limitation would impose resource constraints reflected in the chloroplast genome architecture. Rates of molecular evolution are significantly slower along the phylogenetic branch leading to Ostreobium, in agreement with the expected effects of low light and energy levels on molecular evolution. We expected the ability of Ostreobium to perform photosynthesis in very low light to be associated with positive selection in genes related to the photosynthetic machinery, but instead, we observed that these genes may be under stronger purifying selection. Besides shedding light on the genome dynamics associated with a low light lifestyle, this study helps to resolve the role of environmental factors in shaping the diversity of genome architectures observed in nature.

Obligate insect endosymbionts exhibit increased ortholog length variation and loss of large accessory proteins concurrent with genome shrinkage

Extreme genome reduction has been observed in obligate intracellular insect mutualists and is an assumed consequence of fixed, long-term host isolation. Rapid accumulation of mutations and pseudogenization of genes no longer vital for an intracellular lifestyle, followed by deletion of many genes, are factors that lead to genome reduction. Size reductions in individual genes due to small-scale deletions have also been implicated in contributing to overall genome shrinkage. Conserved protein functional domains are expected to exhibit low tolerance for mutations and therefore remain relatively unchanged throughout protein length reduction while nondomain regions, presumably under less selective pressures, would shorten. This hypothesis was tested using orthologous protein sets from the Flavobacteriaceae (phylum: Bacteroidetes) and Enterobacteriaceae (subphylum: Gammaproteobacteria) families, each of which includes some of the smallest known genomes. Upon examination of protein, functional domain, and nondomain region lengths, we found that proteins were not uniformly shrinking with genome reduction, but instead increased length variability and variability was observed in both the functional domain and nondomain regions. Additionally, as complete gene loss also contributes to overall genome shrinkage, we found that the largest proteins in the proteomes of nonhost-restricted bacteroidetial and gammaproteobacterial species often were inferred to be involved in secondary metabolic processes, extracellular sensing, or of unknown function. These proteins were absent in the proteomes of obligate insect endosymbionts. Therefore, loss of genes encoding large proteins not required for host-restricted lifestyles in obligate endosymbiont proteomes likely contributes to extreme genome reduction to a greater degree than gene shrinkage.

Sequence context of indel mutations and their effect on protein evolution in a bacterial endosymbiont

Indel mutations play key roles in genome and protein evolution, yet we lack a comprehensive understanding of how indels impact evolutionary processes. Genome-wide analyses enabled by next-generation sequencing can clarify the context and effect of indels, thereby integrating a more detailed consideration of indels with our knowledge of nucleotide substitutions. To this end, we sequenced Blochmannia chromaiodes, an obligate bacterial endosymbiont of carpenter ants, and compared it with the close relative, B. pennsylvanicus. The genetic distance between these species is small enough for accurate whole genome alignment but large enough to provide a meaningful spectrum of indel mutations. We found that indels are subjected to purifying selection in coding regions and even intergenic regions, which show a reduced rate of indel base pairs per kilobase compared with nonfunctional pseudogenes. Indels occur almost exclusively in repeat regions composed of homopolymers and multimeric simple sequence repeats, demonstrating the importance of sequence context for indel mutations. Despite purifying selection, some indels occur in protein-coding genes. Most are multiples of three, indicating selective pressure to maintain the reading frame. The deleterious effect of frameshift-inducing indels is minimized by either compensation from a nearby indel to restore reading frame or the indel's location near the 3'-end of the gene. We observed amino acid divergence exceeding nucleotide divergence in regions affected by frameshift-inducing indels, suggesting that these indels may either drive adaptive protein evolution or initiate gene degradation. Our results shed light on how indel mutations impact processes of molecular evolution underlying endosymbiont genome evolution.

Deep sequencing of the transcriptomes of soybean aphid and associated endosymbionts

BACKGROUND

The soybean aphid has significantly impacted soybean production in the U.S. Transcriptomic analyses were conducted for further insight into leads for potential novel management strategies.

METHODOLOGY/PRINCIPAL FINDINGS

Transcriptomic data were generated from whole aphids and from 2,000 aphid guts using an Illumina GAII sequencer. The sequence data were assembled de novo using the Velvet assembler. In addition to providing a general overview, we demonstrate (i) the use of the Multiple-k/Multiple-C method for de novo assembly of short read sequences, followed by BLAST annotation of contigs for increased transcript identification: From 400,000 contigs analyzed, 16,257 non-redundant BLAST hits were identified; (ii) analysis of species distributions of top non-redundant hits: 80% of BLAST hits (minimum e-value of 1.0-E3) were to the pea aphid or other aphid species, representing about half of the pea aphid genes; (iii) comparison of relative depth of sequence coverage to relative transcript abundance for genes with high (membrane alanyl aminopeptidase N) or low transcript abundance; (iv) analysis of the Buchnera transcriptome: Transcripts from 57.6% of the genes from Buchnera aphidicola were identified; (v) identification of Arsenophonus and Wolbachia as potential secondary endosymbionts; (vi) alignment of full length sequences from RNA-seq data for the putative salivary gland protein C002, the silencing of which has potential for aphid management, and the putative Bacillus thuringiensis Cry toxin receptors, aminopeptidase N and alkaline phosphatase.

CONCLUSIONS/SIGNIFICANCE

THIS STUDY PROVIDES THE MOST COMPREHENSIVE DATA SET TO DATE FOR SOYBEAN APHID GENE EXPRESSION: This work also illustrates the utility of short-read transcriptome sequencing and the Multiple-k/Multiple-C method followed by BLAST annotation for rapid identification of target genes for organisms for which reference genome sequences are not available, and extends the utility to include the transcriptomes of endosymbionts.

Structure and dynamics of the operon map of Buchnera aphidicola sp. strain APS

BACKGROUND

Gene expression regulation is still poorly documented in bacteria with highly reduced genomes. Understanding the evolution and mechanisms underlying the regulation of gene transcription in Buchnera aphidicola, the primary endosymbiont of aphids, is expected both to enhance our understanding of this nutritionally based association and to provide an intriguing case-study of the evolution of gene expression regulation in a reduced bacterial genome.

RESULTS

A Bayesian predictor was defined to infer the B. aphidicola transcription units, which were further validated using transcriptomic data and RT-PCR experiments. The characteristics of B. aphidicola predicted transcription units (TUs) were analyzed in order to evaluate the impact of operon map organization on the regulation of gene transcription.On average, B. aphidicola TUs contain more genes than those of E. coli. The global layout of B. aphidicola operon map was mainly shaped by the big reduction and the rearrangements events, which occurred at the early stage of the symbiosis. Our analysis suggests that this operon map may evolve further only by small reorganizations around the frontiers of B. aphidicola TUs, through promoter and/or terminator sequence modifications and/or by pseudogenization events. We also found that the need for specific transcription regulation exerts some pressure on gene conservation, but not on gene assembling in the operon map in Buchnera. Our analysis of the TUs spacing pointed out that a selection pressure is maintained on the length of the intergenic regions between divergent adjacent gene pairs.

CONCLUSIONS

B. aphidicola can seemingly only evolve towards a more polycistronic operon map. This implies that gene transcription regulation is probably subject to weak selection pressure in Buchnera conserving operons composed of genes with unrelated functions.

Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum

Buchnera aphidicola is the primary obligate intracellular symbiont of most aphid species. B. aphidicola and aphids have been evolving in parallel since their association started, about 150 Myr ago. Both partners have lost their autonomy, and aphid diversification has been confined to smaller ecological niches by this co-evolution. B. aphidicola has undergone major genomic and biochemical changes as a result of adapting to intracellular life. Several genomes of B. aphidicola from different aphid species have been sequenced in the last decade, making it possible to carry out analyses and comparative studies using system-level in silico methods. This review attempts to provide a systemic description of the symbiotic function of aphid endosymbionts, particularly of B. aphidicola from the pea aphid Acyrthosiphon pisum, by analyzing their structural genomic properties, as well as their genetic and metabolic networks.

The functional transfer of genes from the mitochondria to the nucleus: the effects of selection, mutation, population size and rate of self-fertilization

The transfer of mitochondrial genes to the nucleus is a recurrent and consistent feature of eukaryotic genome evolution. Although many theories have been proposed to explain such transfers, little relevant data exist. The observation that clonal and self-fertilizing plants transfer more mitochondrial genes to their nuclei than do outcrossing plants contradicts predictions of major theories based on nuclear recombination and leaves a gap in our conceptual understanding how the observed pattern of gene transfer could arise. Here, with a series of deterministic and stochastic simulations, we show how epistatic selection and relative mutation rates of mitochondrial and nuclear genes influence mitochondrial-to-nuclear gene transfer. Specifically, we show that when there is a benefit to having a mitochondrial gene present in the nucleus, but absent in the mitochondria, self-fertilization dramatically increases both the rate and the probability of gene transfer. However, absent such a benefit, when mitochondrial mutation rates exceed those of the nucleus, self-fertilization decreases the rate and probability of transfer. This latter effect, however, is much weaker than the former. Our results are relevant to understanding the probabilities of fixation when loci in different genomes interact.

Evolution of the G+C Content Frontier in the Rat Cytomegalovirus Genome

Within the 230138 bp of the rat cytomegalovirus (RCMV) genome, the G+C content changes abruptly at position 142644, constituting a G+C content frontier. To the left of this point, overall G+C content is 69.2%, and to the right it is only 47.6%. A region of extremely low G+C content (33.8%) is found in the 5 kb immediately to the right of the frontier, in which there are no predicted coding sequences. To the right of position 147501, the G+C content rises and predicted coding sequences reappear. However, these genes are much shorter (average 848 bp, 50% G+C) than those in the left two-thirds of the genome (average 1462 bp, 70% G+C). Whole genome alignment of several viruses indicates that the initial ultra-low G+C region appeared in the common ancestor of the genera Cytomegalovirus and Muromegalovirus, and that the lowering of G+C in the right third has been a subsequent process in the lineage leading to RCMV. The left two-thirds of RCMV has stop codon occurrences at 67.5% of their expected level, based on a modified Markov chain model of stop codon distribution, and the corresponding figure for the right third is 78%. Therefore, despite heavy mutation pressure, selective constraint has operated in the right third of the RCMV genome to maintain a degree of gene length unusual for such low G+C sequences.

The nature and dynamics of bacterial genomes

Though generally small and gene rich, bacterial genomes are constantly subjected to both mutational and population-level processes that operate to increase amounts of functionless DNA. As a result, the coding potential of bacterial genomes can be substantially lower than originally predicted. Whereas only a single pseudogene was included in the original annotation of the bacterium Escherichia coli, we estimate that this genome harbors hundreds of inactivated and otherwise functionless genes. Such regions will never yield a detectable phenotype, but their identification is vital to efforts to elucidate the biological role of all the proteins within the cell.

Recognizing the pseudogenes in bacterial genomes

Pseudogenes are now known to be a regular feature of bacterial genomes and are found in particularly high numbers within the genomes of recently emerged bacterial pathogens. As most pseudogenes are recognized by sequence alignments, we use newly available genomic sequences to identify the pseudogenes in 11 genomes from 4 bacterial genera, each of which contains at least 1 human pathogen. The numbers of pseudogenes range from 27 in Staphylococcus aureus MW2 to 337 in Yersinia pestis CO92 (e.g. 1–8% of the annotated genes in the genome). Most pseudogenes are formed by small frameshifting indels, but because stop codons are A + T-rich, the two low-G + C Gram-positive taxa (Streptococcus and Staphylococcus) have relatively high fractions of pseudogenes generated by nonsense mutations when compared with more G + C-rich genomes. Over half of the pseudogenes are produced from genes whose original functions were annotated as ‘hypothetical’ or ‘unknown’; however, several broadly distributed genes involved in nucleotide processing, repair or replication have become pseudogenes in one of the sequenced Vibrio vulnificus genomes. Although many of our comparisons involved closely related strains with broadly overlapping gene inventories, each genome contains a largely unique set of pseudogenes, suggesting that pseudogenes are formed and eliminated relatively rapidly from most bacterial genomes.

Genome evolution in bacterial endosymbionts of insects

Many insect species rely on intracellular bacterial symbionts for their viability and fecundity. Large-scale DNA-sequence analyses are revealing the forces that shape the evolution of these bacterial associates and the genetic basis of their specialization to an intracellular lifestyle. The full genome sequences of two obligate mutualists, Buchnera aphidicola of aphids and Wigglesworthia glossinidia of tsetse flies, reveal substantial gene loss and an integration of host and symbiont metabolic functions. Further genomic comparisons should reveal the generality of these features among bacterial mutualists and the extent to which they are shared with other intracellular bacteria, including obligate pathogens.

Small genome of Candidatus Blochmannia, the bacterial endosymbiont of Camponotus, implies irreversible specialization to an intracellular lifestyle

Blochmannia (Candidatus Blochmannia gen. nov.) is the primary bacterial endosymbiont of the ant genus CAMPONOTUS: Like other obligate endosymbionts of insects, Blochmannia occurs exclusively within eukaryotic cells and has experienced long-term vertical transmission through host lineages. In this study, PFGE was used to estimate the genome size of Blochmannia as approximately 800 kb, which is significantly smaller than its free-living relatives in the enterobacteria. This small genome implies that Blochmannia has deleted most of the genetic machinery of related free-living bacteria. Due to restricted gene exchange in obligate endosymbionts, the substantial gene loss in Blochmannia and other insect mutualists may reflect irreversible specialization to a host cellular environment.

The genetic properties of the primary endosymbionts of mealybugs differ from those of other endosymbionts of plant sap-sucking insects

Mealybugs (Hemiptera, Coccoidea, Pseudococcidae), like aphids and psyllids, are plant sap-sucking insects that have an obligate association with prokaryotic endosymbionts that are acquired through vertical, maternal transmission. We sequenced two fragments of the genome of Tremblaya princeps, the endosymbiont of mealybugs, which is a member of the beta subdivision of the Proteobacteria. Each of the fragments (35 and 30 kb) contains a copy of 16S-23S-5S rRNA genes. A total of 37 open reading frames were detected, which corresponded to putative rRNA proteins, chaperones, and enzymes of branched-chain amino acid biosynthesis, DNA replication, protein translation, and RNA synthesis. The genome of T. princeps has a number of properties that distinguish it from the genomes of Buchnera aphidicola and Carsonella ruddii, the endosymbionts of aphids and psyllids, respectively. Among these properties are a high G+C content (57.1 mol%), the same G+C content in intergenic spaces and structural genes, and similar G+C contents of the genes encoding highly and poorly conserved proteins. The high G+C content has a substantial effect on protein composition; about one-third of the residues consist of four amino acids with high-G+C-content codons. Sequence analysis of DNA fragments containing the rRNA operon and adjacent regions from endosymbionts of several mealybug species suggested that there was a single duplication of the rRNA operon and the adjacent genes in an ancestor of the present T. princeps. Subsequently, in one mealybug lineage rpS15, one of the duplicated genes, was retained, while in another lineage it decayed. These results extend the diversity of the types of endosymbiotic associations found in plant sap-sucking insects.

Degenerative minimalism in the genome of a psyllid endosymbiont

Psyllids, like aphids, feed on plant phloem sap and are obligately associated with prokaryotic endosymbionts acquired through vertical transmission from an ancestral infection. We have sequenced 37 kb of DNA of the genome of Carsonella ruddii, the endosymbiont of psyllids, and found that it has a number of unusual properties revealing a more extreme case of degeneration than was previously reported from studies of eubacterial genomes, including that of the aphid endosymbiont Buchnera aphidicola. Among the unusual properties are an exceptionally low guanine-plus-cytosine content (19.9%), almost complete absence of intergenic spaces, operon fusion, and lack of the usual promoter sequences upstream of 16S rDNA. These features suggest the synthesis of long mRNAs and translational coupling. The most extreme instances of base compositional bias occur in the genes encoding proteins that have less highly conserved amino acid sequences; the guanine-plus-cytosine content of some protein-coding sequences is as low as 10%. The shift in base composition has a large effect on proteins: in polypeptides of C. ruddii, half of the residues consist of five amino acids with codons low in guanine plus cytosine. Furthermore, the proteins of C. ruddii are reduced in size, with an average of about 9% fewer amino acids than in homologous proteins of related bacteria. These observations suggest that the C. ruddii genome is not subject to constraints that limit the evolution of other known eubacteria.

Psyllid endosymbionts exhibit patterns of co-speciation with hosts and destabilizing substitutions in ribosomal RNA

Eubacterial 16S rDNAs were sequenced from endosymbionts of seven psyllids (Psylloidea) and one whitefly (Aleyrodoidea), to investigate the evolution of endosymbionts and their hosts. Primary endosymbionts from all psyllids formed a highly supported clade, tentatively placed as the sister to whitefly primary endosymbionts, and showing several points of congruence with the host morphological phylogeny. Almost all host taxa yielded an additional eubacterial sequence, related either to known psyllid secondary endosymbionts or to other insect endosymbionts or parasites. The relationships of some secondary endosymbionts also suggested cospeciation with psyllid hosts, or ancient horizontal transfers. All primary endosymbionts, and some secondary endosymbionts, exhibited molecular genetic effects of a long-term, intracellular existence in their biased nucleotide content and decreased stability of rRNA secondary structure.

New insights into symbiotic associations between ants and bacteria

Many ants live in complex mutualistic or parasitic relationships with other insects or plants, some of which are classical examples of the mutual benefit of symbiosis. However, only in the past few years have new insights into the symbiosis of ants and microorganisms been reported. Examples are the symbiosis of ants of the genus Camponotus with intracellular bacteria present in their midgut, and the tripartite relationship of ants of the tribe Attini with an extracellular bacterium essential for the cultivation of their fungus gardens. The analysis of the parasitic and mutualistic interactions of these organisms will allow interesting insights into the evolution of symbiosis and possibly lead to novel strategies of pest control.

Bacterial endosymbionts in animals

Molecular phylogenetic studies reveal that many endosymbioses between bacteria and invertebrate hosts result from ancient infections followed by strict vertical transmission within host lineages. Endosymbionts display a distinctive constellation of genetic properties including AT-biased base composition, accelerated sequence evolution, and, at least sometimes, small genome size; these features suggest increased genetic drift. Molecular genetic characterization also has revealed adaptive, host-beneficial traits such as amplification of genes underlying nutrient provision.