The evolution of the flagellar assembly pathway in endosymbiotic bacterial genomes

Genome shrinkage is a common feature of most intracellular pathogens and symbionts. Reduction of genome sizes is among the best-characterized evolutionary ways of intracellular organisms to save and avoid maintaining expensive redundant biological processes. Endosymbiotic bacteria of insects are examples of biological economy taken to completion because their genomes are dramatically reduced. These bacteria are nonmotile, and their biochemical processes are intimately related to those of their host. Because of this relationship, many of the processes in these bacteria have been either lost or have suffered massive remodeling to adapt to the intracellular symbiotic lifestyle. An example of such changes is the flagellum structure that is essential for bacterial motility and infectivity. Our analysis indicates that genes responsible for flagellar assembly have been partially or totally lost in most intracellular symbionts of gamma-Proteobacteria. Comparative genomic analyses show that flagellar genes have been differentially lost in endosymbiotic bacteria of insects. Only proteins involved in protein export within the flagella assembly pathway (type III secretion system and the basal body) have been kept in most of the endosymbionts, whereas those involved in building the filament and hook of flagella have only in few instances been kept, indicating a change in the functional purpose of this pathway. In some endosymbionts, genes controlling protein-export switch and hook length have undergone functional divergence as shown through an analysis of their evolutionary dynamics. Based on our results, we suggest that genes of flagellum have diverged functionally as to specialize in the export of proteins from the bacterium to the host.     

The structure and evolution of angiosperm nuclear genomes

Despite several decades of investigation, the organization of angiosperm genomes remained largely unknown until very recently. Data describing the sequence composition of large segments of genomes, covering hundreds of kilobases of contiguous sequence, have only become available in the past two years. Recent results indicate commonalities in the characteristics of many plant genomes, including in the structure of chromosomal components like telomeres and centromeres, and in the order and content of genes. Major differences between angiosperms have been associated mainly with repetitive DNAs, both gene families and mobile elements. Intriguing new studies have begun to characterize the dynamic three-dimensional structures of chromosomes and chromatin, and the relationship between genome structure and co-ordinated gene function.     

Correlated rates of synonymous site evolution across plant genomes

Synonymous substitution rates have been shown to vary among evolutionary lineages of both nuclear and organellar genes across a broad range of taxonomic groups. In animals, rate heterogeneity does not appear to be correlated across nuclear and mitochondrial genes. In this paper, we contrast substitution rates in two plant groups and show that grasses evolve more rapidly than palms at synonymous sites in a mitochondrial, a nuclear, and a plastid gene. Furthermore, we show that the relative rates of synonymous substitution between grasses and palms are similar at the three loci. The correlation in synonymous substitution rates across genes is particularly striking because the three genes evolve at very different absolute rates. In contrast, relative rates of nonsynonymous substitution are not conserved among the three genes.      

Increased rates of sequence evolution in endosymbiotic bacteria and fungi with small effective population sizes

Mutualistic, maternally transmitted endosymbiotic microorganisms undergo severe population bottlenecks at each host generation, resulting in a reduction in effective population size (Ne). Previous studies of Buchnera, the primary endosymbiont of aphids, and of several other species of endosymbiotic bacteria have shown that these species exhibit an increase in the rate of substitution of slightly deleterious mutations, among other predicted effects of increased drift due to small Ne, such as reduced codon bias. However, these studies have been limited in taxonomic scope, and it was therefore not clear whether the increase in rate is a general feature of endosymbiont lineages. Here, we test the prediction that a long-term reduction in Ne causes an increase in substitution rate using DNA sequences of the 16S rRNA gene from 13 phylogenetically independent comparisons between taxonomically diverse endosymbiotic microorganisms and their free-living relatives. Maximum likelihood and distance-based methods both indicate a significant increase in substitution rate in a wide range of bacterial and fungal endosymbionts compared to closely related free-living lineages. We use the same data set to test whether 16S genes from endosymbionts display increased A + T content, another indicator of increased genetic drift, and find that there is no significant difference in base composition between endosymbiont and nonendosymbiont 16S genes. However, analysis of an additional data set of whole bacterial genomes demonstrates that, while host-dependent bacteria have significantly increased genomic A + T content, the base content of the 16S gene tends to vary less than that of the whole genome. It is possible that selection for stability of rRNA is strong enough to overcome the effects of drift toward increased A + T content in endosymbiont 16S genes, despite the reduced effective population sizes of these organisms.     

Conservation of plastid sequences in the plant nuclear genome for millions of years facilitates endosymbiotic evolution

The nuclear genome of eukaryotes contains large amounts of cytoplasmic organelle DNA (nuclear integrants of organelle DNA [norgs]). The recent sequencing of many mitochondrial and chloroplast genomes has enabled investigation of the potential role of norgs in endosymbiotic evolution. In this article, we describe a new polymerase chain reaction-based method that allows the identification and evolutionary study of recent and older norgs in a range of eukaryotes. We tested this method in the genus Nicotiana and obtained sequences from seven nuclear integrants of plastid DNA (nupts) totaling 25 kb in length. These nupts were estimated to have been transferred 0.033 to 5.81 million years ago. The spectrum of mutations present in the potential protein-coding sequences compared with the noncoding sequences of each nupt revealed that nupts evolve in a nuclear-specific manner and are under neutral evolution. Indels were more frequent in noncoding regions than in potential coding sequences of former chloroplastic DNA, most probably due to the presence of a higher number of homopolymeric sequences. Unexpectedly, some potential protein-coding sequences within the nupts still contained intact open reading frames for up to 5.81 million years. These results suggest that chloroplast genes transferred to the nucleus have in some cases several millions of years to acquire nuclear regulatory elements and become functional. The different factors influencing this time frame and the potential role of nupts in endosymbiotic gene transfer are discussed.     

Contrasting patterns of nonneutral evolution in proteins encoded in nuclear and mitochondrial genomes

We report that patterns of nonneutral DNA sequence evolution among published nuclear and mitochondrially encoded protein-coding loci differ significantly in animals. Whereas an apparent excess of amino acid polymorphism is seen in most (25/31) mitochondrial genes, this pattern is seen in fewer than half (15/36) of the nuclear data sets. This differentiation is even greater among data sets with significant departures from neutrality (14/15 vs. 1/6). Using forward simulations, we examined patterns of nonneutral evolution using parameters chosen to mimic the differences between mitochondrial and nuclear genetics (we varied recombination rate, population size, mutation rate, selective dominance, and intensity of germ line bottleneck). Patterns of evolution were correlated only with effective population size and strength of selection, and no single genetic factor explains the empirical contrast in patterns. We further report that in Arabidopsis thaliana, a highly self-fertilizing plant with effectively low recombination, five of six published nuclear data sets also exhibit an excess of amino acid polymorphism. We suggest that the contrast between nuclear and mitochondrial nonneutrality in animals stems from differences in rates of recombination in conjunction with a distribution of selective effects. If the majority of mutations segregating in populations are deleterious, high linkage may hinder the spread of the occasional beneficial mutation.     

Concordant evolution of mitochondrial and nuclear genomes in the wheat pathogen Phaeosphaeria nodorum

We compared patterns of mitochondrial restriction fragment length polymorphism (RFLP) diversity with patterns of nuclear RFLP diversity to investigate the effects of selection, gene flow, and sexual reproduction on the population genetic structure and evolutionary history of the wheat pathogen Phaeosphaeria nodorum. A total of 315 fungal isolates from Texas, Oregon, and Switzerland were analyzed using seven nuclear RFLP probes that hybridized to discrete loci and purified mitochondrial DNA that hybridized to the entire mtDNA genome. Forty-two different mitochondrial haplotypes and 298 different nuclear haplotypes were detected. The two most frequent mtDNA haplotypes were present in every population and represented 32% of all isolates. High levels of gene flow, low levels of population subdivision, no evidence for either host specificity or cyto-nuclear disequilibrium were inferred from the analysis of both genomes. The concordance in estimates of these population genetic parameters from both genomes suggests that the two genomes experienced similar degrees of migration, genetic drift and selection.     

Comparative analyses of DNA methylation and sequence evolution using Nasonia genomes

The functional and evolutionary significance of DNA methylation in insect genomes remains to be resolved. Nasonia is well situated for comparative analyses of DNA methylation and genome evolution, since the genomes of a moderately distant outgroup species as well as closely related sibling species are available. Using direct sequencing of bisulfite-converted DNA, we uncovered a substantial level of DNA methylation in 17 of 18 Nasonia vitripennis genes and a strong correlation between methylation level and CpG depletion. Notably, in the sex-determining locus transformer, the exon that is alternatively spliced between the sexes is heavily methylated in both males and females, whereas other exons are only sparsely methylated. Orthologous genes of the honeybee and Nasonia show highly similar relative levels of CpG depletion, despite ~190 My divergence. Densely and sparsely methylated genes in these species also exhibit similar functional enrichments. We found that the degree of CpG depletion is negatively correlated with substitution rates between closely related Nasonia species for synonymous, nonsynonymous, and intron sites. This suggests that mutation rates increase with decreasing levels of germ line methylation. Thus, DNA methylation is prevalent in the Nasonia genome, may participate in regulatory processes such as sex determination and alternative splicing, and is correlated with several aspects of genome and sequence evolution.     

Endosymbiotic evolution: RNA intermediates in endosymbiotic gene transfer

More than a billion years of endosymbiotic evolution has resulted in extensive gene relocation between the genetic compartments of eukaryotic cells. A new study uses chloroplast genome transformation to shed light on the mechanisms involved.     

Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.     

Faster evolutionary rates in endosymbiotic bacteria than in cospeciating insect hosts

The hypothesis of a universal molecular clock holds that divergent lineages exhibit approximately constant rates of nucleotide substitution over evolutionary time for a particular macromolecule. We compare divergences of ribosomal DNA for aphids (Insecta) and Buchnera, the maternally transmitted, endosymbiotic bacteria that have cospeciated with aphids since initially infecting them over 100 million years ago. Substitution rates average 36 times greater for Buchnera than for their aphid hosts for regions of small-subunit rDNA that are homologous for prokaryotes and eukaryotes. Aphids exhibit 18S rDNA substitution rates that are within the range observed in related insects. In contrast, 16S rDNA evolves about twice as fast in Buchnera as in related free-living bacterial lineages. Nonetheless, the difference between Buchnera and aphids is much greater, suggesting that rates may be generally higher in bacteria. This finding adds to evidence that molecular clocks are only locally rather than universally valid among taxonomic groups. It is consistent with the hypothesis that rates of sequence evolution depend on generation time.     

Relative patterns and rates of evolution in heron nuclear and mitochondrial DNA

Mitochondrial cytochrome b sequence data from 15 species of herons (Aves: Ardeidae), representing 13 genera, were compared with DNA hybridization data of single-copy nuclear DNA (scnDNA) from the same species in a taxonomic congruence assessment of heron phylogeny. The two data sets produced a partially resolved, completely congruent estimate of phylogeny with the following basic structure: (Tigrisoma, Cochlearius, (((Zebrilus, (Ixobrychus, Botaurus)), (((Ardea, Casmerodius), Bubulcus), ((Egretta thula, Egretta caerulea, Egretta tricolor), Syrigma), Butorides, Nycticorax, Nyctanassa)))). Because congruence indicated similar phylogenetic information in the two data sets, we used the relatively unsaturated DNA hybridization distances as surrogates of time to examine graphically the patterns and rates of change in cytochrome b distances. Cytochrome b distances were computed either from whole sequences or from partitioned sequences consisting of transitions, transversions, specific codon site positions, or specific protein-coding regions. These graphical comparisons indicated that unpartitioned cytochrome b has evolved at 5-10 times the rate of scnDNA. Third-position transversions appeared to offer the most useful sequence partition for phylogenetic analysis because of their relatively fast rate of substitution (two times that of scnDNA) and negligible saturation. We also examined lineage-based rates of evolution by comparing branch length patterns between the nuclear and cytochrome b trees. The degree of correlation in corresponding branch lengths between cytochrome b and DNA hybridization trees depended on DNA sequence partitioning. When cytochrome b sequences were not partitioned, branch lengths in the cytochrome b and DNA hybridization trees were not correlated. However, when cytochrome b sequences were reduced to third-position transversions (i.e., unsaturated, relatively fast changing data), branch lengths were correlated. This finding suggests that lineage-based rates of DNA evolution in nuclear and mitochondrial genomes are influenced by common causes.     

Rapid evolution of enormous, multichromosomal genomes in flowering plant mitochondria with exceptionally high mutation rates

Genome size and complexity vary tremendously among eukaryotic species and their organelles. Comparisons across deeply divergent eukaryotic lineages have suggested that variation in mutation rates may explain this diversity, with increased mutational burdens favoring reduced genome size and complexity. The discovery that mitochondrial mutation rates can differ by orders of magnitude among closely related angiosperm species presents a unique opportunity to test this hypothesis. We sequenced the mitochondrial genomes from two species in the angiosperm genus Silene with recent and dramatic accelerations in their mitochondrial mutation rates. Contrary to theoretical predictions, these genomes have experienced a massive proliferation of noncoding content. At 6.7 and 11.3 Mb, they are by far the largest known mitochondrial genomes, larger than most bacterial genomes and even some nuclear genomes. In contrast, two slowly evolving Silene mitochondrial genomes are smaller than average for angiosperms. Consequently, this genus captures approximately 98% of known variation in organelle genome size. The expanded genomes reveal several architectural changes, including the evolution of complex multichromosomal structures (with 59 and 128 circular-mapping chromosomes, ranging in size from 44 to 192 kb). They also exhibit a substantial reduction in recombination and gene conversion activity as measured by the relative frequency of alternative genome conformations and the level of sequence divergence between repeat copies. The evolution of mutation rate, genome size, and chromosome structure can therefore be extremely rapid and interrelated in ways not predicted by current evolutionary theories. Our results raise the hypothesis that changes in recombinational processes, including gene conversion, may be a central force driving the evolution of both mutation rate and genome structure.     

Adaptive evolution in GroEL from distantly related endosymbiotic bacteria of insects

Many symbioses between bacteria and insects resulted from ancient infections followed by strict vertical transmission within host lineages. The strong bottlenecks under which this transmission occurs promote the neutral fixation of slightly deleterious mutations by genetic drift. As predicted by Muller's ratchet, this fixation will drive endosymbiotic bacteria through an irreversible dynamics of fitness loss. The chaperonin GroEL has been proposed as a compensatory mechanism whereby endosymbiotic bacteria of aphids persist. Here, we show that endosymbiotic bacteria of insects from two phylogenetically very distant bacterial phyla have fixed amino acid substitutions by positive selection in functionally important GroEL regions involved in either GroES/peptide binding or in the en bloc movement of the GroEL apical domain. These results, together with the high levels of constitutive expression of GroEL in these endosymbionts, provide valuable insights into the evolution of a molecular mechanism responsible for the maintenance of the symbiotic lifestyle.     

Relative rates of nuclear DNA evolution in human and Old World monkey lineages

Rates of substitution were compared between humans and Old World monkeys for sequences in or adjacent to 19 genes. The comparison of 21,299 sites in noncoding regions indicates that the substitution rate is approximately 43% greater in the Old World monkey lineage. However, 83% of the compared sites are in the region of the beta-globin gene family. Outside this region there is no consistent pattern of rate difference between the two lineages. Comparison of the coding regions of 16 genomically dispersed genes, involving 1,592 synonymous sites and 5,275 nonsynonymous sites, showed a faster rate of substitution in the human lineage at the nonsynonymous sites of the prion gene, but otherwise no evidence of rate difference between the two lineages. It is concluded that rate differences between these two lineages may be specific to certain regions of the genome rather than being a general phenomenon. This conclusion needs to be confirmed by comparison of a larger number of genomically dispersed sequences. It is, however, consistent with the results of DNA-DNA hybridization experiments, which show no difference in evolutionary rate between the two lineages.      

Retrotransposon evolution in diverse plant genomes

Retrotransposon or retrotransposon-like sequences have been reported to be conserved components of cereal centromeres. Here we show that the published sequences are derived from a single conventional Ty3-gypsy family or a nonautonomous derivative. Both autonomous and nonautonomous elements are likely to have colonized Poaceae centromeres at the time of a common ancestor but have been maintained since by active retrotransposition. The retrotransposon family is also present at a lower copy number in the Arabidopsis genome, where it shows less pronounced localization. The history of the family in the two types of genome provides an interesting contrast between "boom and bust" and persistent evolutionary patterns.