Mutation exposed: a neutral explanation for extreme base composition of an endosymbiont genome

Genome-centric view of the microbiome in a new deep-sea glass sponge species Bathydorus sp

Sponges are widely distributed in the global ocean and harbor diverse symbiotic microbes with mutualistic relationships. However, sponge symbionts in the deep sea remain poorly studied at the genome level. Here, we report a new glass sponge species of the genus Bathydorus and provide a genome-centric view of its microbiome. We obtained 14 high-quality prokaryotic metagenome-assembled genomes (MAGs) affiliated with the phyla Nitrososphaerota, Pseudomonadota, Nitrospirota, Bdellovibrionota, SAR324, Bacteroidota, and Patescibacteria. In total, 13 of these MAGs probably represent new species, suggesting the high novelty of the deep-sea glass sponge microbiome. An ammonia-oxidizing Nitrososphaerota MAG B01, which accounted for up to 70% of the metagenome reads, dominated the sponge microbiomes. The B01 genome had a highly complex CRISPR array, which likely represents an advantageous evolution toward a symbiotic lifestyle and forceful ability to defend against phages. A sulfur-oxidizing Gammaproteobacteria species was the second most dominant symbiont, and a nitrite-oxidizing Nitrospirota species could also be detected, but with lower relative abundance. Bdellovibrio species represented by two MAGs, B11 and B12, were first reported as potential predatory symbionts in deep-sea glass sponges and have undergone dramatic genome reduction. Comprehensive functional analysis indicated that most of the sponge symbionts encoded CRISPR-Cas systems and eukaryotic-like proteins for symbiotic interactions with the host. Metabolic reconstruction further illustrated their essential roles in carbon, nitrogen, and sulfur cycles. In addition, diverse putative phages were identified from the sponge metagenomes. Our study expands the knowledge of microbial diversity, evolutionary adaption, and metabolic complementarity in deep-sea glass sponges.

Adaptation of Translational Machinery in Malaria Parasites to Accommodate Translation of Poly-Adenosine Stretches Throughout Its Life Cycle

Malaria is caused by unicellular apicomplexan parasites of the genus Plasmodium, which includes the major human parasite Plasmodium falciparum. The complex cycle of the malaria parasite in both mosquito and human hosts has been studied extensively. There is tight control of gene expression in each developmental stage, and at every level of gene synthesis: from RNA transcription, to its subsequent translation, and finally post-translational modifications of the resulting protein. Whole-genome sequencing of P. falciparum has laid the foundation for significant biological advances by revealing surprising genomic information. The P. falciparum genome is extremely AT-rich (∼80%), with a substantial portion of genes encoding intragenic polyadenosine (polyA) tracks being expressed throughout the entire parasite life cycle. In most eukaryotes, intragenic polyA runs act as negative regulators of gene expression. Recent studies have shown that translation of mRNAs containing 12 or more consecutive adenosines results in ribosomal stalling and frameshifting; activating mRNA surveillance mechanisms. In contrast, P. falciparum translational machinery can efficiently and accurately translate polyA tracks without activating mRNA surveillance pathways. This unique feature of P. falciparum raises interesting questions: (1) How is P. falciparum able to efficiently and correctly translate polyA track transcripts, and (2) What are the specifics of the translational machinery and mRNA surveillance mechanisms that separate P. falciparum from other organisms? In this review, we analyze possible evolutionary shifts in P. falciparum protein synthesis machinery that allow efficient translation of an AU rich-transcriptome. We focus on physiological and structural differences of P. falciparum stage specific ribosomes, ribosome-associated proteins, and changes in mRNA surveillance mechanisms throughout the complete parasite life cycle, with an emphasis on the mosquito and liver stages.

Selective advantages favour high genomic AT-contents in intracellular elements

Extrachromosomal genetic elements such as bacterial endosymbionts and plasmids generally exhibit AT-contents that are increased relative to their hosts’ DNA. The AT-bias of endosymbiotic genomes is commonly explained by neutral evolutionary processes such as a mutational bias towards increased A+T. Here we show experimentally that an increased AT-content of host-dependent elements can be selectively favoured. Manipulating the nucleotide composition of bacterial cells by introducing A+T-rich or G+C-rich plasmids, we demonstrate that cells containing GC-rich plasmids are less fit than cells containing AT-rich plasmids. Moreover, the cost of GC-rich elements could be compensated by providing precursors of G+C, but not of A+T, thus linking the observed fitness effects to the cytoplasmic availability of nucleotides. Accordingly, introducing AT-rich and GC-rich plasmids into other bacterial species with different genomic GC-contents revealed that the costs of G+C-rich plasmids decreased with an increasing GC-content of their host’s genomic DNA. Taken together, our work identifies selection as a strong evolutionary force that drives the genomes of intracellular genetic elements toward higher A+T contents.

Genomes of endosymbiotic bacteria are commonly more AT-rich than the ones of their free-living relatives. Interestingly, genomes of other intracellular elements like plasmids or bacteriophages also tend to be richer in AT than the genomes of their hosts. The AT-bias of endosymbiotic genomes is commonly explained by neutral evolutionary processes. However, since A+T nucleotides are both more abundant and energetically less expensive than G+C nucleotides, an alternative explanation is that selective advantages drive the nucleotide composition of intracellular elements. Here we provide strong experimental evidence that intracellular elements, whose genome is more AT-rich than the genome of the host, are selectively favoured on the host level. Thus, our results emphasize the importance of selection for shaping the DNA base composition of extrachromosomal genetic elements.

Plasmodium parasites of birds have the most AT-rich genes of eukaryotes

The genomic architecture of organisms, including nucleotide composition, can be highly variable, even among closely-related species. To better understand the causes leading to structural variation in genomes, information on distinct and diverse genomic features is needed. Malaria parasites are known for encompassing a wide range of genomic GC-content and it has long been thought that Plasmodium falciparum, the virulent malaria parasite of humans, has the most AT-biased eukaryotic genome. Here, I perform comparative genomic analyses of the most AT-rich eukaryotes sequenced to date, and show that the avian malaria parasites Plasmodium gallinaceum, P. ashfordi, and P. relictum have the most extreme coding sequences in terms of AT-bias. Their mean GC-content is 21.21, 21.22 and 21.60 %, respectively, which is considerably lower than the transcriptome of P. falciparum (23.79 %) and other eukaryotes. This information enables a better understanding of genome evolution and raises the question of how certain organisms are able to prosper despite severe compositional constraints.

Genome Evolution in the Obligate but Environmentally Active Luminous Symbionts of Flashlight Fish

The luminous bacterial symbionts of anomalopid flashlight fish are thought to be obligately dependent on their hosts for growth and share several aspects of genome evolution with unrelated obligate symbionts, including genome reduction. However, in contrast to most obligate bacteria, anomalopid symbionts have an active environmental phase that may be important for symbiont transmission. Here we investigated patterns of evolution between anomalopid symbionts compared with patterns in free-living relatives and unrelated obligate symbionts to determine if trends common to obligate symbionts are also found in anomalopid symbionts. Two symbionts, "Candidatus Photodesmus katoptron" and "Candidatus Photodesmus blepharus," have genomes that are highly similar in gene content and order, suggesting genome stasis similar to ancient obligate symbionts present in insect lineages. This genome stasis exists in spite of the symbiont's inferred ability to recombine, which is frequently lacking in obligate symbionts with stable genomes. Additionally, we used genome comparisons and tests of selection to infer which genes may be particularly important for the symbiont's ecology compared with relatives. In keeping with obligate dependence, substitution patterns suggest that most symbiont genes are experiencing relaxed purifying selection compared with relatives. However, genes involved in motility and carbon storage, which are likely to be used outside the host, appear to be under increased purifying selection. Two chemoreceptor chemotaxis genes are retained by both species and show high conservation with amino acid sensing genes, suggesting that the bacteria may actively seek out hosts using chemotaxis toward amino acids, which the symbionts are not able to synthesize.

Endosymbiont evolution: predictions from theory and surprises from genomes

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

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

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

Neurospora and the dead-end hypothesis: genomic consequences of selfing in the model genus

It is becoming increasingly evident that adoption of different reproductive strategies, such as sexual selfing and asexuality, greatly impacts genome evolution. In this study, we test theoretical predictions on genomic maladaptation of selfing lineages using empirical data from the model fungus Neurospora. We sequenced the genomes of four species representing distinct transitions to selfing within the history of the genus, as well as the transcriptome of one of these, and compared with available data from three outcrossing species. Our results provide evidence for a relaxation of purifying selection in protein-coding genes and for a reduced efficiency of transposable element silencing by Repeat Induced Point mutation. A reduction in adaptive evolution was also identified in the form of reduced codon usage bias in highly expressed genes of selfing Neurospora, but this result may be confounded by mutational bias. Potentially counteracting these negative effects, the nucleotide substitution rate and the spread of transposons is reduced in selfing species. We suggest that differences in substitution rate relate to the absence, in selfing Neurospora, of the asexual pathway producing conidia. Our results support the dead-end theory and show that Neurospora genomes bear signatures of both sexual and asexual reproductive mode.

Evolution of optimal codon choices in the family Enterobacteriaceae

The Enterobacteriaceae are a large family of Proteobacteria that include many well-known prokaryotic genera, such as Escherichia, Yersinia and Salmonella. The main ideas of synonymous codon usage (CU) evolution and translational selection have been deeply influenced by studies with these bacterial groups. In this work we report the analysis of the CU pattern of completely sequenced bacterial genomes that belong to the Enterobacteriaceae. The effect of selection in translation acting at the levels of speed and accuracy, and phylogenetic trends within this group are described. Preferred (optimal) codons were identified. The evolutionary dynamics of these codons were studied and following a Bayesian approach these preferences were traced back to the common ancestor of the family. We found that there is some level of variation in selection among the analysed micro-organisms that is probably associated with lineage-specific trends. The codon bias was largely conserved across the evolutionary time of the family in highly expressed genes and protein conserved regions, suggesting a major role of negative selection. In this sense, the results support the idea that the extant CU bias is finely tuned over the ancestral well-conserved pool of tRNAs.

Keeping it simple: flowering plants tend to retain, and revert to, simple leaves

• A wide range of factors (developmental, physiological, ecological) with unpredictable interactions control variation in leaf form. Here, we examined the distribution of leaf morphologies (simple and complex forms) across angiosperms in a phylogenetic context to detect patterns in the directions of changes in leaf shape. • Seven datasets (diverse angiosperms and six nested clades, Sapindales, Apiales, Papaveraceae, Fabaceae, Lepidium, Solanum) were analysed using maximum likelihood and parsimony methods to estimate asymmetries in rates of change among character states. • Simple leaves are most frequent among angiosperm lineages today, were inferred to be ancestral in angiosperms and tended to be retained in evolution (stasis). Complex leaves slowly originated ('gains') and quickly reverted to simple leaves ('losses') multiple times, with a significantly greater rate of losses than gains. Lobed leaves may be a labile intermediate step between different forms. The nested clades showed mixed trends; Solanum, like the angiosperms in general, had higher rates of losses than gains, but the other clades had higher rates of gains than losses. • The angiosperm-wide pattern could be taken as a null model to test leaf evolution patterns in particular clades, in which patterns of variation suggest clade-specific processes that have yet to be investigated fully.

A slowly evolving host moves first in symbiotic interactions

Symbiotic relationships, both parasitic and mutualistic, are ubiquitous in nature. Understanding how these symbioses evolve, from bacteria and their phages to humans and our gut microflora, is crucial in understanding how life operates. Often, symbioses consist of a slowly evolving host species with each host only interacting with its own subpopulation of symbionts. The Red Queen hypothesis describes coevolutionary relationships as constant arms races with each species rushing to evolve an advantage over the other, suggesting that faster evolution is favored. Here, we use a simple game theoretic model of host-symbiont coevolution that includes population structure to show that if the symbionts evolve much faster than the host, the equilibrium distribution is the same as it would be if it were a sequential game where the host moves first against its symbionts. For the slowly evolving host, this will prove to be advantageous in mutualisms and a handicap in antagonisms. The result follows from rapid symbiont adaptation to its host and is robust to changes in the parameters, even generalizing to continuous and multiplayer games. Our findings provide insight into a wide range of symbiotic phenomena and help to unify the field of coevolutionary theory.

Massive genomic decay in Serratia symbiotica, a recently evolved symbiont of aphids

All vertically transmitted bacterial symbionts undergo a process of genome reduction over time, resulting in tiny, gene-dense genomes. Comparison of genomes of ancient bacterial symbionts gives only limited information about the early stages in the transition from a free-living to symbiotic lifestyle because many changes become obscured over time. Here, we present the genome sequence for the recently evolved aphid symbiont Serratia symbiotica. The S. symbiotica genome exhibits several of the hallmarks of genome evolution observed in more ancient symbionts, including elevated rates of evolution and reduction in genome size. The genome also shows evidence for massive genomic decay compared with free-living relatives in the same genus of bacteria, including large deletions, many pseudogenes, and a slew of rearrangements, perhaps promoted by mobile DNA. Annotation of pseudogenes allowed examination of the past and current metabolic capabilities of S. symbiotica and revealed a somewhat random process of gene inactivation with respect to function. Analysis of mutational patterns showed that deletions are more common in neutral DNA. The S. symbiotica genome provides a rare opportunity to study genome evolution in a recently derived heritable symbiont.

Evidence of selection upon genomic GC-content in bacteria

The genomic GC-content of bacteria varies dramatically, from less than 20% to more than 70%. This variation is generally ascribed to differences in the pattern of mutation between bacteria. Here we test this hypothesis by examining patterns of synonymous polymorphism using datasets from 149 bacterial species. We find a large excess of synonymous GC→AT mutations over AT→GC mutations segregating in all but the most AT-rich bacteria, across a broad range of phylogenetically diverse species. We show that the excess of GC→AT mutations is inconsistent with mutation bias, since it would imply that most GC-rich bacteria are declining in GC-content; such a pattern would be unsustainable. We also show that the patterns are probably not due to translational selection or biased gene conversion, because optimal codons tend to be AT-rich, and the excess of GC→AT SNPs is observed in datasets with no evidence of recombination. We therefore conclude that there is selection to increase synonymous GC-content in many species. Since synonymous GC-content is highly correlated to genomic GC-content, we further conclude that there is selection on genomic base composition in many bacteria.

Forces that influence the evolution of codon bias

The frequencies of alternative synonymous codons vary both among species and among genes from the same genome. These patterns have been inferred to reflect the action of natural selection. Here we evaluate this in bacteria. While intragenomic variation in many species is consistent with selection favouring translationally optimal codons, much of the variation among species appears to be due to biased patterns of mutation. The strength of selection on codon usage can be estimated by two different approaches. First, the extent of bias in favour of translationally optimal codons in highly expressed genes, compared to that in genes where selection is weak, reveals the long-term effectiveness of selection. Here we show that the strength of selected codon usage bias is highly correlated with bacterial growth rate, suggesting that selection has favoured translational efficiency. Second, the pattern of bias towards optimal codons at polymorphic sites reveals the ongoing action of selection. Using this approach we obtained results that were completely consistent with the first method; importantly, the frequency spectra of optimal codons at polymorphic sites were similar to those predicted under an equilibrium model. Highly expressed genes in Escherichia coli appear to be under continuing strong selection, whereas selection is very weak in genes expressed at low levels.

Evidence for homologous recombination in intracellular chemosynthetic clam symbionts

Homologous recombination is a fundamental mechanism for the genetic diversification of free-living bacteria. However, recombination may be limited in endosymbiotic bacteria, as these taxa are locked into an intracellular niche and may rarely encounter sources of foreign DNA. This study tested the hypothesis that vertically transmitted endosymbionts of deep-sea clams (Bivalvia: Vesicomyidae) show little or no evidence of recombination. Phylogenetic analysis of 13 loci distributed across the genomes of 14 vesicomyid symbionts revealed multiple, well-supported inconsistencies among gene tree topologies, and maximum likelihood-based tests rejected a hypothesis of shared evolutionary history (linkage) among loci. Further, multiple statistical methods confirmed the presence of recombination by detecting intragenic breakpoints in two symbiont loci. Recombination may be confined to a subset of vesicomyid symbionts, as some clades showed high levels of genomic stability, whereas others showed clear patterns of homologous exchange. Notably, a mosaic genome is present in symB, a symbiont lineage shown to have been acquired laterally (i.e., nonvertically) by Vesicomya sp. JdF clams. The majority of loci analyzed here supported a tight sister clustering of symB with the symbiont of a host species from the Mid-Atlantic Ridge, whereas others placed symB in a clade with symA, the dominant phylotype of V. sp. JdF clams. This result raises the hypothesis that lateral symbiont transfer between hosts may facilitate recombination by bringing divergent symbiont lineages into contact. Together, the data show that homologous recombination contributes to the diversification of vesicomyid clam symbionts, despite the intracellular lifestyle of these bacteria.

Global distribution and evolution of a toxinogenic Burkholderia-Rhizopus symbiosis

Toxinogenic endobacteria were isolated from a collection of Rhizopus spp. representing highly diverse geographic origins and ecological niches. All endosymbionts belonged to the Burkholderia rhizoxinica complex according to matrix-assisted laser desorption ionization-time of flight biotyping and multilocus sequence typing, suggesting a common ancestor. Comparison of host and symbiont phylogenies provides insights into possible cospeciation and horizontal-transmission events.

The temporal dynamics of slightly deleterious mutations in Escherichia coli and Shigella spp

The Shigella are recently emerged clones of Escherichia coli, which have independently adopted an intracellular pathogenic lifestyle. We examined the molecular evolutionary consequences of this niche specialization by comparing the normalized, directional frequency profiles of unique polymorphisms within 2,098 orthologues representing the intersection of five E. coli and four Shigella genomes. We note a surfeit of AT-enriching changes (GC-->AT), transversions, and nonsynonymous changes in the Shigella genomes. By examining these differences within a temporal framework, we conclude that our results are consistent with relaxed or inefficient selection in Shigella owing to a reduced effective population size. Alternative interpretations, and the interesting exception of Shigella sonnei, are discussed. Finally, this analysis lends support to the view that nucleotide composition typically does not lie at mutational equilibrium but that selection plays a role in maintaining a higher GC content than would result solely from mutation bias. This argument sheds light on the enrichment of adenine and thymine in the genomes of bacterial endosymbionts where purifying selection is very weak.

The proteomic constraint and its role in molecular evolution

Recently, the concept of a "Proteomic Constraint" was introduced to explain the frequency of genetic code deviations in mitochondrial genomes. The Proteomic Constraint was proposed to be proportional to the size of the mitochondrially encoded proteome, hence small proteomes are expected to experience smaller total numbers of errors resulting from genetic code deviations, leading to less likelihood of causing lethality. The concept is now extended to encompass several other aspects of the genetic information system. When the Proteomic Constraint is small, it is proposed that there is little selective pressure to evolve or maintain error correction mechanisms, as a result of the smaller total number of errors that accumulate. Conversely, a large Proteomic Constraint is proposed to result in a correspondingly large selective pressure to evolve or maintain error correction mechanisms. Differences in the size of the Proteomic Constraint can help to explain differences in replicational, transcriptional, and translational fidelities between genomes. A key piece of evidence is the existence of negative power law relationships between proteome size and error rates; these are demonstrated to be diagnostic of the action of the Proteomic Constraint. The Proteomic Constraint is argued to be a major factor determining mutation rates in a diverse range of DNA genomes, implying that mutation rates are clock like. A small Proteomic Constraint partly explains why RNA viruses possess high mutation rates. A reduced Proteomic Constraint in intracellular pathogenic bacteria predicts a drift upwards in mutation rates. Differences in the Proteomic Constraint also appear to be linked to differences in recombination rates between eukaryotes. In addition, a reduced Proteomic Constraint may explain features of resident genomes, such as loss of DNA repair pathways, increased substitution rates, and AT biases, in addition to the occurrence of genetic code deviations. Thus, it is argued that the Proteomic Constraint is a universal factor that influences a wide range of properties of the genetic information system.

An effect of 16S rRNA intercistronic variability on coevolutionary analysis in symbiotic bacteria: molecular phylogeny of Arsenophonus triatominarum

The genes of ribosomal RNA are the most popular and frequently used markers for bacterial phylogeny and reconstruction of insect-symbiont coevolution. In primary symbionts, such as Buchnera and Wigglesworthia, genome economization leads to the establishment of a single copy of these sequences. In phylogenetic studies, they provide sufficient information and yield phylogenetic trees congruent with host evolution. In contrast, other symbiotic lineages (e.g., the genus Arsenophonus) carry a higher number of rRNA copies in their genomes, which may have serious consequences for phylogenetic inference. In this study, we show that in Arsenophonus triatominarum the degree of heterogeneity can affect reconstruction of phylogenetic relationships and mask possible coevolution between the symbiont and its host. Phylogenetic arrangement of individual rRNA copies was used, together with a calculation of their divergence time, to demonstrate that the incongruent 16S rDNA trees and low nucleotide diversity in the secondary symbiont could be reconciled with the coevolutionary scenario.

Population structure, levels of selection, and the evolution of intracellular symbionts

Many obligately intracellular symbionts exhibit a characteristic set of genetic changes that include an increase in substitution rates, loss of many genes, and apparent destabilization of many proteins and structural RNAs. Authors have suggested that these changes are due to increased mutation rates, or, more commonly, decreased effective population size due to population bottlenecks at the symbiont or, perhaps, host level. I propose that the increase in substitution rates and accumulation of deleterious mutations is a consequence of the population structure imposed on the endosymbionts by strict host association, loss of horizontal transmission and potentially conflicting levels of selection. I analyze a population genetic model of endosymbiont evolution, and demonstrate that substitution rates will increase, and the effect of those substitutions on endosymbiont fitness will become more deleterious as horizontal transmission among hosts decreases. Additionally, I find that there is a critical level of horizontal transmission below which natural selection cannot effectively purge deleterious mutations, leading to an expected loss of fitness over time. This critical level varies across loci with the degree of correlation between host and endosymbiont fitness, and may help explain differential retention and loss of certain genes.

The rise and fall of deleterious mutation

It is well established that selection is less efficient in small populations than in large ones. Here we review the impact of this effect by considering the gradual selective purging of deleterious mutation over time. We outline an approach to explore the dynamics of this process, and highlight its profound implications.

Molecular evidence for polyphyletic origin of the primary symbionts of sucking lice (phthiraptera, anoplura)

Based on 16S rDNA analyses, the primary symbionts of sucking lice were found to form a polyphyletic assemblage of several distant lineages that have arisen several times within Enterobacteriaceae and at least once within Legionellaceae. Another independent lineage of endosymbiotic enterobacteria inhabits a sister group of the sucking lice, Rhynchophthirina. The inspection of 16S rDNA supports the symbiotic nature of the investigated bacteria; they display a typical trait of degenerative processes, an increased AT content (Adenine-Thymine content) in comparison with free-living bacteria. The calculation of divergence time between the closest anopluran and rhynchophthirine symbionts further support their independent origin. The results shown here, together with evidence from other groups, indicate that the significance of primary symbionts for blood-feeding insects should be reconsidered.

Significant positive correlation between the recombination rate and GC content in the human pseudoautosomal region

This paper establishes that recombination drives the evolution of GC content in a significant way. Because the human P-arm pseudoautosomal region (PAR1) has been shown to have a high recombination rate, at least 20-fold more frequent than the genomic average of approximately 1 cM/Mb, this region provides an ideal system to study the role of recombination in the evolution of base composition. Nine non-coding regions of PAR1 are analyzed in this study. We have observed a highly significant positive correlation between the recombination rate and GC content (rho = 0.837, p < or = 0.005). Five regions that lie in the distal part of PAR1 are shown to be significantly higher than genomic average divergence. By comparing the intra- and inter-specific AT->GC -GC->AT ratios, we have detected no fixation bias toward GC alleles except for L254915, which has excessive AT-->GC changes in the human lineage. Thus, we conclude that the high GC content of the PAR1 genes better fits the biased gene conversion (BGC) model.

Lineage-specific variations of congruent evolution among DNA sequences from three genomes, and relaxed selective constraints on rbcL in Cryptomonas (Cryptophyceae)

Background

Plastid-bearing cryptophytes like Cryptomonas contain four genomes in a cell, the nucleus, the nucleomorph, the plastid genome and the mitochondrial genome. Comparative phylogenetic analyses encompassing DNA sequences from three different genomes were performed on nineteen photosynthetic and four colorless Cryptomonas strains. Twenty-three rbcL genes and fourteen nuclear SSU rDNA sequences were newly sequenced to examine the impact of photosynthesis loss on codon usage in the rbcL genes, and to compare the rbcL gene phylogeny in terms of tree topology and evolutionary rates with phylogenies inferred from nuclear ribosomal DNA (concatenated SSU rDNA, ITS2 and partial LSU rDNA), and nucleomorph SSU rDNA.

Results

Largely congruent branching patterns and accelerated evolutionary rates were found in nucleomorph SSU rDNA and rbcL genes in a clade that consisted of photosynthetic and colorless species suggesting a coevolution of the two genomes. The extremely accelerated rates in the rbcL phylogeny correlated with a shift from selection to mutation drift in codon usage of two-fold degenerate NNY codons comprising the amino acids asparagine, aspartate, histidine, phenylalanine, and tyrosine. Cysteine was the sole exception. The shift in codon usage seemed to follow a gradient from early diverging photosynthetic to late diverging photosynthetic or heterotrophic taxa along the branches. In the early branching taxa, codon preferences were changed in one to two amino acids, whereas in the late diverging taxa, including the colorless strains, between four and five amino acids showed changes in codon usage.

Conclusion

Nucleomorph and plastid gene phylogenies indicate that loss of photosynthesis in the colorless Cryptomonas strains examined in this study possibly was the result of accelerated evolutionary rates that started already in photosynthetic ancestors. Shifts in codon usage are usually considered to be caused by changes in functional constraints and in gene expression levels. Thus, the increasing influence of mutation drift on codon usage along the clade may indicate gradually relaxed constraints and reduced expression levels on the rbcL gene, finally correlating with a loss of photosynthesis in the colorless Cryptomonas paramaecium strains.