Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.

Maize as a model for the evolution of plant nuclear genomes

The maize genome is replete with chromosomal duplications and repetitive DNA. The duplications resulted from an ancient polyploid event that occurred over 11 million years ago. Based on DNA sequence data, the polyploid event occurred after the divergence between sorghum and maize, and hence the polyploid event explains some of the difference in DNA content between these two species. Genomic rearrangement and diploidization followed the polyploid event. Most of the repetitive DNA in the maize genome is retrotransposable elements, and they comprise 50% of the genome. Retrotransposon multiplication has been relatively recent—within the last 5–6 million years—suggesting that the proliferation of retrotransposons has also contributed to differences in DNA content between sorghum and maize. There are still unanswered questions about repetitive DNA, including the distribution of repetitive DNA throughout the genome, the relative impacts of retrotransposons and chromosomal duplication in plant genome evolution, and the hypothesized correlation of duplication events with transposition. Population genetic processes also affect the evolution of genomes. We discuss how centromeric genes should, in theory, contain less genetic diversity than noncentromeric genes. In addition, studies of diversity in the wild relatives of maize indicate that different genes have different histories and also show that domestication and intensive breeding have had heterogeneous effects on genetic diversity across genes.

Genetic fidelity under harsh conditions: Analysis of spontaneous mutation in the thermoacidophilic archaeon Sulfolobus acidocaldarius

Microbes whose genomes are encoded by DNA and for which adequate information is available display similar genomic mutation rates (average 0.0034 mutations per chromosome replication, range 0.0025 to 0.0046). However, this value currently is based on only a few well characterized microbes reproducing within a narrow range of environmental conditions. In particular, no genomic mutation rate has been determined either for a microbe whose natural growth conditions may extensively damage DNA or for any member of the archaea, a prokaryotic lineage deeply diverged from both bacteria and eukaryotes. Both of these conditions are met by the extreme thermoacidophile Sulfolobus acidocaldarius. We determined the genomic mutation rate for this species when growing at pH 3.5 and 75°C based on the rate of forward mutation at the pyrE gene and the nucleotide changes identified in 101 independent mutants. The observed value of about 0.0018 extends the range of DNA-based microbes with rates close to the standard rate simultaneously to an archaeon and to an extremophile whose cytoplasmic pH and normal growth temperature greatly accelerate the spontaneous decomposition of DNA. The mutations include base pair substitutions (BPSs) and additions and deletions of various sizes, but the S. acidocaldarius spectrum differs from those of other DNA-based organisms in being relatively poor in BPSs. The paucity of BPSs cannot yet be explained by known properties of DNA replication or repair enzymes of Sulfolobus spp. It suggests, however, that molecular evolution per genome replication may proceed more slowly in S. acidocaldarius than in other DNA-based organisms examined to date.

Genes linked by fusion events are generally of the same functional category: A systematic analysis of 30 microbial genomes

Recent work in computational genomics has shown that a functional association between two genes can be derived from the existence of a fusion of the two as one continuous sequence in another genome. For each of 30 completely sequenced microbial genomes, we established all such fusion links among its genes and determined the distribution of links within and among 15 broad functional categories. We found that 72% of all fusion links related genes of the same functional category. A comparison of the distribution of links to simulations on the basis of a random model further confirmed the significance of intracategory fusion links. Where a gene of annotated function is linked to an unclassified gene, the fusion link suggests that the two genes belong to the same functional category. The predictions based on fusion links are shown here for Methanobacterium thermoautotrophicum, and another 661 predictions are available at http://fusion.bu.edu.

A minimal gene set for cellular life derived by comparison of complete bacterial genomes.

The recently sequenced genome of the parasitic bacterium Mycoplasma genitalium contains only 468 identified protein-coding genes that have been dubbed a minimal gene complement [Fraser, C.M., Gocayne, J.D., White, O., Adams, M.D., Clayton, R.A., et al. (1995) Science 270, 397-403]. Although the M. genitalium gene complement is indeed the smallest among known cellular life forms, there is no evidence that it is the minimal self-sufficient gene set. To derive such a set, we compared the 468 predicted M. genitalium protein sequences with the 1703 protein sequences encoded by the other completely sequenced small bacterial genome, that of Haemophilus influenzae. M. genitalium and H. influenzae belong to two ancient bacterial lineages, i.e., Gram-positive and Gram-negative bacteria, respectively. Therefore, the genes that are conserved in these two bacteria are almost certainly essential for cellular function. It is this category of genes that is most likely to approximate the minimal gene set. We found that 240 M. genitalium genes have orthologs among the genes of H. influenzae. This collection of genes falls short of comprising the minimal set as some enzymes responsible for intermediate steps in essential pathways are missing. The apparent reason for this is the phenomenon that we call nonorthologous gene displacement when the same function is fulfilled by nonorthologous proteins in two organisms. We identified 22 nonorthologous displacements and supplemented the set of orthologs with the respective M. genitalium genes. After examining the resulting list of 262 genes for possible functional redundancy and for the presence of apparently parasite-specific genes, 6 genes were removed. We suggest that the remaining 256 genes are close to the minimal gene set that is necessary and sufficient to sustain the existence of a modern-type cell. Most of the proteins encoded by the genes from the minimal set have eukaryotic or archaeal homologs but seven key proteins of DNA replication do not. We speculate that the last common ancestor of the three primary kingdoms had an RNA genome. Possibilities are explored to further reduce the minimal set to model a primitive cell that might have existed at a very early stage of life evolution.

Similarities and dissimilarities of phage genomes.

Genomic similarities and contrasts are investigated in a collection of 23 bacteriophages, including phages with temperate, lytic, and parasitic life histories, with varied sequence organizations and with different hosts and with different morphologies. Comparisons use relative abundances of di-, tri-, and tetranucleotides from entire genomes. We highlight several specific findings. (i) As previously shown for cellular genomes, each viral genome has a distinctive signature of short oligonucleotide abundances that pervade the entire genome and distinguish it from other genomes. (ii) The enteric temperate double-stranded (ds) phages, like enterobacteria, exhibit significantly high relative abundances of GpC = GC and significantly low values of TA, but no such extremes exist in ds lytic phages. (iii) The tetranucleotide CTAG is of statistically low relative abundance in most phages. (iv) The DAM methylase site GATC is of statistically low relative abundance in most phages, but not in P1. This difference may relate to controls on replication (e.g., actions of the host SeqA gene product) and to MutH cleavage potential of the Escherichia coli DAM mismatch repair system. (v) The enteric temperate dsDNA phages form a coherent group: they are relatively close to each other and to their bacteria] hosts in average differences of dinucleotide relative abundance values. By contrast, the lytic dsDNA phages do not form a coherent group. This difference may come about because the temperate phages acquire more sequence characteristics of the host because they use the host replication and repair machinery, whereas the analyzed lytic phages are replicated by their own machinery. (vi) The nonenteric temperate phages with mycoplasmal and mycobacterial hosts are relatively close to their respective hosts and relatively distant from any of the enteric hosts and from the other phages. (vii) The single-stranded RNA phages have dinucleotide relative abundance values closest to those for random sequences, presumably attributable to the mutation rates of RNA phages being much greater than those of DNA phages.

Conserved motifs in prokaryotic and eukaryotic polypeptide release factors: tRNA-protein mimicry hypothesis.

Translation termination requires two codon-specific polypeptide release factors in prokaryotes and one omnipotent factor in eukaryotes. Sequences of 17 different polypeptide release factors from prokaryotes and eukaryotes were compared. The prokaryotic release factors share residues split into seven motifs. Conservation of many discrete, perhaps critical, amino acids is observed in eukaryotic release factors, as well as in the C-terminal portion of elongation factor (EF) G. Given that the C-terminal domains of EF-G interacts with ribosomes by mimicry of a tRNA structure, the pattern of conservation of residues in release factors may reflect requirements for a tRNA-mimicry for binding to the A site of the ribosome. This mimicry would explain why release factors recognize stop codons and suggests that all prokaryotic and eukaryotic release factors evolved from the progenitor of EF-G.

The genomic and physical organization of Ty1-copia-like sequences as a component of large genomes in Pinus elliottii var. elliottii and other gymnosperms.

A DNA sequence, TPE1, representing the internal domain of a Ty1-copia retroelement, was isolated from genomic DNA of Pinus elliottii Engelm. var. elliottii (slash pine). Genomic Southern analysis showed that this sequence, carrying partial reverse transcriptase and integrase gene sequences, is highly amplified within the genome of slash pine and part of a dispersed element >4.8 kbp. Fluorescent in situ hybridization to metaphase chromosomes shows that the element is relatively uniformly dispersed over all 12 chromosome pairs and is highly abundant in the genome. It is largely excluded from centromeric regions and intercalary chromosomal sites representing the 18S-5.8S-25S rRNA genes. Southern hybridization with specific DNA probes for the reverse transcriptase gene shows that TPE1 represents a large subgroup of heterogeneous Ty1-copia retrotransposons in Pinus species. Because no TPE1 transcription could be detected, it is most likely an inactive element--at least in needle tissue. Further evidence for inactivity was found in recombinant reverse transcriptase and integrase sequences. The distribution of TPE1 within different gymnosperms that contain Ty1-copia group retrotransposons, as shown by a PCR assay, was investigated by Southern hybridization. The TPE1 family is highly amplified and conserved in all Pinus species analyzed, showing a similar genomic organization in the three- and five-needle pine species investigated. It is also present in spruce, bald cypress (swamp cypress), and in gingko but in fewer copies and a different genomic organization.

In vitro synthesis of oncogenic human papillomaviruses requires episomal genomes for differentiation-dependent late expression.

Human papillomavirus (HPV) types 16, 18, 31, and 51 are the etiologic agents of many anogenital cancers including those of the cervix. These "high risk" HPVs specifically target genital squamous epithelia, and their lytic life cycle is closely linked to epithelial differentiation. We have developed a genetic assay for HPV functions during pathogenesis using recircularized cloned HPV 31 genomes that were transfected together with a drug resistance marker into monolayer cultures of normal human foreskin keratinocytes, the natural host cell. After drug selection, cell lines were isolated that stably maintained HPV 31 DNA as episomes and underwent terminal differentiation when grown in organotypic raft cultures. In differentiated rafts, the expression of late viral genes, amplification of viral DNA, and production of viral particles were detected in suprabasal cells. This demonstrated the ability to synthesize HPV 31 virions from transfected DNA templates and allowed an examination of HPV functions during the vegetative viral life cycle. We then used this system to investigate whether an episomal genome was required for the induction of late viral gene expression. When an HPV 31 genome (31E1*) containing a missense mutation in the E1 open reading frame was transfected into normal human keratinocytes, the mutant viral sequences were found to integrate into the host cell chromosomal DNA with both early and late regions intact. While high levels of early viral gene transcription were observed, no late gene expression was detected in rafts of cell lines containing the mutant viral genome despite evidence of terminal differentiation. Therefore, the induction of late viral gene expression required that the viral genomes be maintained as extrachromosomal elements, and terminal differentiation alone was not sufficient. These studies provide the basis for a detailed examination of HPV functions during viral pathogenesis.

The amino-terminal domain of the prokaryotic enhancer-binding protein XylR is a specific intramolecular repressor.

The mechanism under which the signal-reception amino-terminal portion (A domain) of the prokaryotic enhancer-binding protein XylR controls the activity of the regulator has been investigated through complementation tests in vivo, in which the various protein segments were produced as independent polypeptides. Separate expression of the A domain repressed the otherwise constitutive activity of a truncated derivative of XylR deleted of its A domain (XylR delta A). Such inhibition was not released by m-xylene, the natural inducer of the system. Repression caused by the A domain was specific for XylR because it did not affect activation of the sigma 54 promoter PnifH by a derivative of its cognate regulator, NifA, deleted of its own A domain. The A domain was also unable to repress the activity of a NifA-XylR hybrid protein resulting from fusing two-thirds of the central domain of NifA to the carboxyl-terminal third of XylR, which includes its DNA-binding domain. The inhibitory effect caused by the A domain of XylR on XylR delta A seems, therefore, to result from specific interactions in trans between the two truncated proteins and not from mere hindering of an activating surface.

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines.

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.

Different cellular backgrounds confer a marked advantage to either mutant or wild-type mitochondrial genomes.

After the introduction of mitochondria with a mixture of mutant and wild-type mitochondrial DNA (mtDNA) into a human rho degree cell line (143B.206), Yoneda et al. [Yoneda, M., Chomyn, A., Martinuzzi, A., Hurko, O. & Attardi, G. (1992) Proc. Natl. Acad. Sci. USA 89, 11164-11168] observed a shift in the proportion of the two mitochondrial genotypes in a number of cybrid clones. In every case where a shift was observed, there was an increase in the proportion of mutant mtDNA. By using the same cell line (143B.206 rho degree), we also generated cybrids that were either stable in their mitochondrial genotype or showed an increase in the proportion of mutant mtDNA. However, temporal analysis of the same mutant mtDNA type in another rho degree cell line revealed a quite distinct outcome. Those clones that showed a change shifted toward higher levels of wild-type rather than mutant mtDNA. These results indicate that the nuclear genetic background of the recipient (rho degree) cell can influence the segregation of mutant and wild-type mitochondrial genomes in cell cybrids.

A prokaryotic origin for light-dependent chlorophyll biosynthesis of plants.

Flowering plants require light for chlorophyll synthesis. Early studies indicated that the dependence on light for greening stemmed in part from the light-dependent reduction of the chlorophyll intermediate protochlorophyllide to the product chlorophyllide. Light-dependent reduction of protochlorophyllide by flowering plants is contrasted by the ability of nonflowering plants, algae, and photosynthetic bacteria to reduce protochlorophyllide and, hence, synthesize (bacterio) chlorophyll in the dark. In this report, we functionally complemented a light-independent protochlorophyllide reductase mutant of the eubacterium Rhodobacter capsulatus with an expression library composed of genomic DNA from the cyanobacterium Synechocystis sp. PCC 6803. The complemented R. capsulatus strain is capable of synthesizing bacteriochlorophyll in the light, thereby indicating that a chlorophyll biosynthesis enzyme can function in the bacteriochlorophyll biosynthetic pathway. However, under dark growth conditions the complemented R. capsulatus strain fails to synthesize bacteriochlorophyll and instead accumulates protochlorophyllide. Sequence analysis demonstrates that the complementing Synechocystis genomic DNA fragment exhibits a high degree of sequence identity (53-56%) with light-dependent protochlorophyllide reductase enzymes found in plants. The observation that a plant-type, light-dependent protochlorophyllide reductase enzyme exists in a cyanobacterium indicates that light-dependent protochlorophyllide reductase evolved before the advent of eukaryotic photosynthesis. As such, this enzyme did not arise to fulfill a function necessitated either by the endosymbiotic evolution of the chloroplast or by multicellularity; rather, it evolved to fulfill a fundamentally cell-autonomous role.

Identification of two flavivirus-like genomes in the GB hepatitis agent.

A subtractive PCR methodology known as representational difference analysis was used to clone specific nucleotide sequences present in the infectious plasma from a tamarin infected with the GB hepatitis agent. Eleven unique clones were identified, seven of which were examined extensively. All seven clones appeared to be derived from sequences exogenous to the genomes of humans, tamarins, Saccharomyces cerevisiae, and Escherichia coli. In addition, sequences from these clones were not detected in plasma or liver tissue of tamarins prior to their inoculation with the GB agent. These sequences were detected by reverse transcription-PCR in acute-phase plasma of tamarins inoculated with the GB agent. Probes derived from two of the seven clones detected an RNA species of > or = 8.3 kb in the liver of a GB-agent-infected tamarin by Northern blot hybridization. Sequence analysis indicated that five of the seven clones encode polypeptides that possess limited amino acid identity with the nonstructural proteins of hepatitis C virus. Extension of the sequences found in the seven clones revealed that plasma from an infected tamarin contained two RNA molecules > 9 kb long. Limited sequence identity with various isolates of hepatitis C virus and the relative positions of putative RNA helicases and RNA-dependent RNA polymerases in the predicted protein products of these molecules suggested that the GB agent contains two unique flavivirus-like genomes.