Characterization of two rice genes for nuclear-encoded chloroplast ribosomal protein L12 and phylogenetic analysis of the acquisition of transit peptides and gene duplication

 We have identified two genes coding for chloroplast ribosomal protein L12 encoded in the nuclear genome of rice (Oryza sativa). These genes were designated rpl12-1 and rpl12-2 (rpl12, ribosomal protein L12). Northern analysis with specific probes revealed that both genes are transcribed. The expression of each gene seems to have a different regulatory machinery. It is also suggested that the expression of rpl12-1 is controlled in an organ-specific manner. The deduced amino-acid sequences of the mature peptide parts are more conserved than those of the transit peptide parts in both monocotyledonous and dicotyledonous plants. A phylogenetic tree was constructed using the nucleotide sequences of the transit peptide region of the rpl12s of reported plant. The tree includes estimates of when the transit peptides were acquired, and when the genes were duplicated, in the course of evolution. According to our hypothesis, the nuclear-translocated chloroplast ribosomal protein L12 gene obtained its transit peptide after the divergence of monocots and dicots, then gene duplications occurred independently in monocots and dicots, and subsequently rice and rye branched apart. Received: 4 October 1997 / Accepted: 5 January 1998     

Examination of four HLA class I pseudogenes. Common events in the evolution of HLA genes and pseudogenes.

The HLA class I gene family in lymphoblastoid cell line 721 has been studied in detail and a number of sequences in addition to the classical genes have been identified. The cloning, characterization, and nucleotide sequences of four sequences, all full length HLA class I pseudogenes, are described in this report. These pseudogenes, contained within 5.4-, 5.9-, 7.0-, and 9.2-kb HindIII fragments, each have the class I exon-intron structure as well as class I homology in their 5' and 3' flanking regions. However, all four sequences have one or more substitutions that perturb the coding region, leaving little doubt that they are in fact pseudogenes. Comparisons among these sequences and the HLA class I genes revealed that their homology with the class I genes is patchwork. Thus, although some regions have diverged, other contiguous intron-exon sequences are highly conserved. Comparisons in the 5' regions indicate that the pseudogene promoters more closely resemble the classical HLA promoters than the nonclassical promoters as none of the unique structural features found in the HLA-E, -F, or -G regulatory regions are present in any of the pseudogene promoters. Further comparisons revealed that at least two putative gene conversion events, similar to those hypothesized to have occurred in the evolution of some HLA genes, may have occurred in the evolution of some of the pseudogenes. These and other hypothetical events in the evolution of the class I gene family are discussed.     

Evolutionary dynamics of wheat mitochondrial gene structure with special remarks on the origin and effects of RNA editing in cereals

We investigated the evolutionary dynamics of wheat mitochondrial genes with respect to their structural differentiation during organellar evolution, and to mutations that occurred during cereal evolution. First, we compared the nucleotide sequences of three wheat mitochondrial genes to those of wheat chloroplast, alpha-proteobacterium and cyanobacterium orthologs. As a result, we were able to (1) differentiate the conserved and variable segments of the orthologs, (2) reveal the functional importance of the conserved segments, and (3) provide a corroborative support for the alpha-proteobacterial and cyanobacterial origins of those mitochondrial and chloroplast genes, respectively. Second, we compared the nucleotide sequences of wheat mitochondrial genes to those of rice and maize to determine the types and frequencies of base changes and indels occurred in cereal evolution. Our analyses showed that both the evolutionary speed, in terms of number of base substitutions per site, and the transition/transversion ratio of the cereal mitochondrial genes were less than two-fifths of those of the chloroplast genes. Eight mitochondrial gene groups differed in their evolutionary variability, RNA and Complex I (nad) genes being most stable whereas Complex V (atp) and ribosomal protein genes most variable. C-to-T transition was the most frequent type of base change; C-to-G and G-to-C transversions occurred at lower rates than all other changes. The excess of C-to-T transitions was attributed to C-to-U RNA editing that developed in early stage of vascular plant evolution. On the contrary, the editing of C residues at cereal T-to-C transition sites developed mostly during cereal divergence. Most indels were associated with short direct repeats, suggesting intra- and intermolecular recombination as an important mechanism for their origin. Most of the repeats associated with indels were di- or trinucleotides, although no preference was noticed for their sequences. The maize mt genome was characterized by a high incidence of indels, comparing to the wheat and rice mt genomes.     

A pseudogene cluster in the leader region of the Euglena chloroplast 16S-23S rRNA genes.

The nucleotide sequence of a region (leader region) preceding the 5'-end of 16S-23S rRNA gene region of Euglena gracilis chloroplast DNA was compared with the homologous sequences that code for the 16S-23S rRNA operons of Euglena and E. coli. The leader region shows close homology in sequence to the 16S-23S rRNA gene region of Euglena (Orozco et al. (1980) J. Biol.Chem. 255, 10997-11003) as well as to the rrnD operon of E. coli, suggesting that it was derived from the 16S-23S rRNA gene region by gene duplication. It was shown that the leader region had accumulated nucleotide substitutions at an extremely rapid rate in its entirety, similar to the rate of tRNAIle pseudogene identified in the leader region. In addition, the leader region shows an unique base content which is quite distinct from those of 16S-23S rRNA gene regions of Euglena and E. coli, but again is similar to that of the tRNAIle pseudogene. The above two results strongly suggest that the leader region contains a pseudogene cluster which was derived from a gene cluster coding for the functional 16S-23S rRNA operon possibly by imperfect duplication during evolution of Euglena chloroplast DNA.     

The nucleotide sequence of five ribosomal protein genes from the cyanelles ofCyanophora paradoxa: Implications concerning the phylogenetic relationship between cyanelles and chloroplasts

Summary The nucleotide sequences of the ribosomal protein genesrps18, rps19, rpl2, rpl33, and partial sequence ofrpl22 from cyanelles, the photosynthetic organelles of the protistCyanophora paradoxa, have been determined. These genes form two clusters oriented in opposite and divergent directions. One cluster contains therpl33 andrps18 genes; the other contains therpl2, rps19, andrpl22 genes, in that order. Phylogenetic trees were constructed from both the DNA sequences and the deduced protein sequences of cyanelles,Euglena gracilis and land plant chloroplasts, andEscherichia coli, using parsimony or maximum likelihood methods. In addition, a phylogenetic tree was built from a distance matrix comparing the number of nucleotide substitutions per site. The phylogeny inferred from all these methods suggests that cyanelles fall within the chloroplast line of evolution and that the evolutionary distances between cyanelles and land plant chloroplasts are shorter than betweenE. gracilis chloroplasts and land plant chloroplasts.     

Molecular evolution and phylogeny of the atpB-rbcL spacer of chloroplast DNA in the true mosses.

The nucleotide variation of a noncoding region between the atpB and rbcL genes of the chloroplast genome was used to estimate the phylogeny of 11 species of true mosses (subclass Bryidae). The A+T rich (82.6%) spacer sequence is conserved with 48% of bases showing no variation between the ingroup and outgroup. Rooted at liverworts, Marchantia and Bazzania, the monophyly of true mosses was supported cladistically and statistically. A nonparametric Wilcoxon Signed-Ranks test Ts statistic for testing the taxonomic congruence showed no significant differences between gene trees and organism trees as well as between parsimony trees and neighbor-joining trees. The reconstructed phylogeny based on the atpB-rbcL spacer sequences indicated the validity of the division of acrocarpous and pleurocarpous mosses. The size of the chloroplast spacer in mosses fits into an evolutionary trend of increasing spacer length from liverworts through ferns to seed plants. According to the relative rate tests, the hypothesis of a molecular clock was supported in all species except for Thuidium, which evolved relatively fast. The evolutionary rate of the chloroplast DNA spacer in mosses was estimated to be (1.12 +/- 0.019) x 10(-10) nucleotides per site per year, which is close to the nonsynonymous substitution rates of the rbcL gene in the vascular plants. The constrained molecular evolution (total nucleotide substitutions, K approximately 0.0248) of the chloroplast DNA spacer is consistent with the slow evolution in morphological traits of mosses. Based on the calibrated evolutionary rate, the time of the divergence of true mosses was estimated to have been as early as 220 million years ago.     

Molecular analysis of a 21.1-kb fragment of wheat chloroplast DNA bearing RNA polymerase subunit (rpo) genes.

The entire nucleotide sequence of a 21.1-kb fragment of wheat chloroplast (ct) DNA was determined. This fragment carries 18 intact genes and parts of two additional genes, including the three RNA polymerase genes rpoB, rpoC1, and rpoC2. The gene arrangement of this region is conserved in wheat, rice, and maize, but not in non-grass species. Comparison of these 20 genes in wheat, rice, and maize showed that tRNA genes evolved more slowly than protein-coding genes in the chloroplast genome. Intergenic regions evolved much faster than both types of genes. Although the 19 genes of wheat, except for orf42, showed high identity to those of other plants, there were three novel structural features in the wheat rpoC2 gene; a deletion of 81 bp in the middle region, a variable insertion (408 bp), and a nonsense mutation in the 3' terminal region, resulting in truncation of a sequence of ca. 10 amino acids. An intermolecular recombination between the stretches of CTTAT and CTTTT was suggested as the mechanism of the 81-bp deletion in the wheat rpoC2 gene. Evolutionary distance between the chloroplast genomes of wheat and maize was larger than those between wheat and rice and between rice and maize.     

Noncoding sequences from the slowly evolving chloroplast inverted repeat in addition to rbcL data do not support gnetalean affinities of angiosperms

We developed PCR primers against highly conserved regions of the rRNA operon located within the inverted repeat of the chloroplast genome and used these to amplify the region spanning from the 3' terminus of the 23S rRNA gene to the 5' terminus of the 5S rRNA gene. The sequence of this roughly 500-bp region, which includes the 4.5S rRNA gene and two chloroplast intergenic transcribed spacer regions (cpITS2 and cpITS3), was determined from 20 angiosperms, 7 gymnosperms, and 16 ferns (21,700 bp). Sequences for the large subunit of ribulose bisphosphate carboxylase/oxygenase (rbcL) from the same or confamilial genera were analyzed in both separate and combined data sets. Due to the low substitution rate in the inverted repeat region, noncoding sequences in the cpITS region are not saturated with substitutions, in contrast to synonymous sites in rbcL, which are shown to evolve roughly six times faster than noncoding cpITS sequences. Several length polymorphisms with very clear phylogenetic distributions were detected in the data set. Results of phylogenetic analyses provide very strong bootstrap support for monophyly of both spermatophytes and angiosperms. No support for a sister group relationship between Gnetales and angiosperms in either cpITS or rbcL data was found. Rather, weak bootstrap support for monophyly of gymnosperms studied and for a basal position for the aquatic angiosperm Nymphaea among angiosperms studied was observed. Noncoding sequences from the inverted repeat region of chloroplast DNA appear suitable for study of land plant evolution.      

Recombination diversifies chloroplast trnF pseudogenes in Arabidopsis lyrata

Extensive intraspecific variation in the chloroplast trnL(UAA)-trnF(GAA) spacer of model plant Arabidopsis lyrata is caused by multiple copies of a tandemly repeated trnF pseudogene undergoing parallel independent changes in copy number. Linkage disequilibrium and secondary structure analyses indicate that the diversification of pseudogene copies is driven by complex processes of structurally mediated illegitimate recombination. Disperse repeats sharing similar secondary structures interact, facilitating reciprocal exchange of structural motifs between copies via intramolecular and intermolecular recombinations, forming chimeric sequences and iterative expansion and contraction in pseudogene copy numbers. Widely held assumptions that chloroplast sequence evolution is simple and structural changes are informative are violated. Our findings have important implications for the use of this highly variable region in Brassicaceae studies. The reticulate evolution and nonindependent nucleotide substitution render the pseudogene inappropriate for standard phylogenetic reconstruction, but over short evolutionary timescales they may be useful for assessing gene flow, hybridization and introgression.     

Mitochondrial gene transfer in pieces: fission of the ribosomal protein gene rpl2 and partial or complete gene transfer to the nucleus.

Mitochondrial genes are usually conserved in size in angiosperms. A notable exception is the rpl2 gene, which is considerably shorter in the eudicot Arabidopsis than in the monocot rice. Here, we show that a severely truncated mitochondrial rpl2 gene (termed 5' rpl2) was created by the formation of a premature stop codon early in eudicot evolution. This 5' rpl2 gene was subsequently lost many times from the mitochondrial DNAs of 179 core eudicots surveyed by Southern hybridization. The sequence corresponding to the 3' end of rice rpl2 (termed 3' rpl2) has been lost much more pervasively among the mitochondrial DNAs of core eudicots than has 5' rpl2. Furthermore, where still present in these mitochondrial genomes, 3' rpl2 always appears to be a pseudogene, and there is no evidence that 3' rpl2 was ever a functional mitochondrial gene. An intact and expressed 3' rpl2 gene was discovered in the nucleus of five diverse eudicots (tomato, cotton, Arabidopsis, soybean, and Medicago). In the first three of these species, 5' rpl2 is still present in the mitochondrion, unlike the two legumes, where both parts of rpl2 are present in the nucleus as separate genes. The full-length rpl2 gene has been transferred intact to the nucleus in maize. We propose that the 3' end of rpl2 was functionally transferred to the nucleus early in eudicot evolution, and that this event then permitted the nonsense mutation that gave rise to the mitochondrial 5' rpl2 gene. Once 5' rpl2 was established as a stand-alone mitochondrial gene, it was then lost, and was probably transferred to the nucleus many times. This complex history of gene fission and gene transfer has created four distinct types of rpl2 structures or compartmentalizations in angiosperms: (1) intact rpl2 gene in the mitochondrion, (2) intact gene in the nucleus, (3) split gene, 5' in the mitochondrion and 3' in the nucleus, and (4) split gene, both parts in the nucleus.     

Molecular evolution of chloroplast DNA sequences

Comparative data on the evolution of chloroplast genes are reviewed. The chloroplast genome has maintained a similar structural organization over most plant taxa so far examined. Comparisons of nucleotide sequence divergence among chloroplast genes reveals marked similarity across the plant kingdom and beyond to the cyanobacteria (blue-green algae). Estimates of rates of nucleotide substitution indicate a synonymous rate of 1.1 x 10(-9) substitutions per site per year. Noncoding regions also appear to be constrained in their evolution, although addition/deletion events are common. There have also been evolutionary changes in the distribution of introns in chloroplast encoded genes. Relative to mammalian mitochondrial DNA, the chloroplast genome evolves at a conservative rate.      

Whole chloroplast genome comparison of rice,maize, and wheat: implications for chloroplast gene diversification and phylogeny of cereals

The fully sequenced chloroplast genomes of maize (subfamily Panicoideae), rice (subfamily Bambusoideae), and wheat (subfamily Pooideae) provide the unique opportunity to investigate the evolution of chloroplast genes and genomes in the grass family (Poaceae) by whole-genome comparison. Analyses of nucleotide sequence variations in 106 cereal chloroplast genes with tobacco sequences as the outgroup suggested that (1) most of the genic regions of the chloroplast genomes of maize, rice, and wheat have evolved at similar rates; (2) RNA genes have highly conservative evolutionary rates relative to the other genes; (3) photosynthetic genes have been under strong purifying selection; (4) between the three cereals, 14 genes which account for about 28% of the genic region have evolved with heterogeneous nucleotide substitution rates; and (5) rice genes tend to have evolved more slowly than the others at loci where rate heterogeneity exists. Although the mechanism that underlies chloroplast gene diversification is complex, our analyses identified variation in nonsynonymous substitution rates as a genetic force that generates heterogeneity, which is evidence of selection in chloroplast gene diversification at the intrafamilial level. Phylogenetic trees constructed with the variable nucleotide sites of the chloroplast genes place maize basal to the rice-wheat clade, revealing a close relationship between the Bambusoideae and Pooideae.     

A study of the phylogeny of Brassica rapa,B. nigra,Raphanus sativus,and their related genera using noncoding regions of chloroplast DNA

There are two evolutionary lineages in the genus Brassica: the rapa/oleracea lineage and the nigra lineage. Using nuclear DNA sequences such as the intergenic spacer between 5S rRNA genes and the internal transcribed spacer between 18S and 25S rRNA genes, we and others had previously demonstrated that Raphanus sativus is closely related to the nigra lineage. In the present study, we sequenced the chloroplast noncoding region between trnT and trnF and that between trnD and trnT in seven species and showed that R. sativus is more closely related to the rapa/oleracea lineage than to the nigra lineage. The conflicting results from nuclear DNA and chloroplast DNA support the hypothesis that Raphanus was derived from a hybridization between the rapa/oleracea and the nigra lineages. We estimated the date of this hybridization event to be 60% of the divergence time between the two Brassica lineages. In addition, the pattern and rate of nucleotide substitution were studied. There were more transversions than transitions in these noncoding regions, which have a high AT content. Furthermore, the proportion of transversions among the substitutions at a site increases with increasing A + T content of its two adjacent nucleotides. An influence of immediate 5(') pyrimidine on substitution pattern is also observed when both adjacent bases in the two DNA strands are A or T. The rate of nucleotide substitution in the trnL group I intron is only about one third of the rate in the nearby intergenic spacers in the trnT-trnF fragment. The rate of nucleotide substitution in the rapa/oleracea lineage is at least 1.5 times that in the nigra lineage.     

Fine structural features of the chloroplast genome: comparison of the sequenced chloroplast genomes.

The entire nucleotide sequences of the rice, tobacco and liverwort chloroplast genomes have been determined. We compared all the chloroplast genes, open reading frames and spacer regions in the plastid genomes of these three species in order to elucidate general structural features of the chloroplast genome. Analyses of homology, GC content and codon usage of the genes enabled us to classify them into two groups: photosynthesis genes and genetic system genes. Based on comparisons of homology, GC content and codon usage, unidentified ORFs can also be assigned to each of these groups such that it is possible to speculate about the functions of products which may be produced by these ORFs. The spacer regions and intron sequences were compared and found to have no obvious homology between rice and liverwort or between tobacco and liverwort.     

Molecular evolution of mammalian lactate dehydrogenase-A genes and pseudogenes: association of a mouse processed pseudogene with a B1 repetitive sequence

A mouse genomic clone containing a lactate dehydrogenase-A (LDH-A) processed pseudogene and a B1 repetitive element was isolated, and a nucleotide sequence of approximately 3 kb was determined. The pseudogene and B1 element are flanked by perfect 13-bp repeats, and the B1 sequence starts at 14 nucleotides 3' to the presumptive polyadenylation signal of the pseudogene. The nucleotide sequences of the LDH-A genes and processed pseudogenes from mouse, rat, and human were compared, and a phylogenetic tree was constructed. The rate and pattern of nucleotide substitutions in the LDH-A pseudogenes are similar to previously reported results (Li et al. 1984). The average rate of nucleotide substitutions in the LDH-A pseudogenes is 4.3 X 10(- 9)/site/year. The substitutions of C----T and G----A are most frequent, and A----G substitutions are relatively high. The rate of synonymous substitutions in the LDH-A genes is 5.3 X 10(-9), which is not significantly higher than the average rate of 4.7 X 10(-9) for 35 mammalian genes. The rate of nonsynonymous substitutions in the LDH-A genes is 0.20 X 10(-9), which is considerably lower than the average rate of 0.88 X 10(-9) for 35 mammalian genes. Thus, the mammalian LDH-A gene appears to be highly conserved in evolution.