Transfer RNAs and tRNA genes of Vicia faba chloroplasts

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

A mutation hotspot in the chloroplast genome of a conifer (Douglas-fir: Pseudotsuga) is caused by variability in the number of direct repeats derived from a partially duplicated tRNA gene

We determined the DNA sequence of a 2.7-kb cpDNA XbaI fragment from douglas-fir [Pseudotsuga menziesii (Mirb.) Franco]. RFLPs revealed by the 2.7-kb XbaI clone were observed to vary up to 1 kb among species within the genus Pseudotsuga and up to 200 bp among trees of P. menziesii. The polymerase chain reaction (PCR) allowed the locus of polymorphism to be identified, and the variable region was then sequenced in a second Douglas-fir tree, a single tree of a related species, Japanese Douglas-fir (P. japonica), and in a species lacking a mutation hotspot in the region, Pinus radiata (Monterey pine). The locus of polymorphism is characterized by hundreds of base pairs of imperfect, tandem direct repeats flanked by a partially duplicated and an intact trn Y-GUA gene. The duplication is direct in orientation and consists of 43 bp of the 3' end of trnY and 25 bp of its 3' flanking sequence. Tandem repeats show high sequence similarity to a 27-bp region of the trnY gene that overlaps one end of the duplication. The two trees of Douglas-fir sequenced differed by a single tandem repeat unit, whereas these trees differed from the Japanese Douglas-fir sequenced by approximately 34 repeat units. Repetitive DNA in the Pseudotsuga cpDNA hotspot was most likely generated at the time of the partial trnY gene duplication and these sequences expanded by slipped-strand mispairing and unequal crossing-over.

Codon recognition mechanisms in plant chloroplasts

In chloroplasts, all 61 sense codons are found in chloroplast (cp) DNA sequences coding for proteins. However among the sequenced cp tRNAs or tRNA genes, tRNAs with anticodons complementary to codons CUU/C (Leu), CCU/C (Pro), GCU/C (Ala) and CGC/A/G (Arg) [or CGC/A (Arg) in Marchantia] have not been found. In this paper we show that cp tRNA(Ala)(U*GC) cp tRNA(Pro)(U*GG) and cp tRNA(Arg)(ICG) are able to decode the corresponding four-codon family. In the case of leucine codons CUU/C, we show that 'U:U and U:C wobble' mechanisms can operate to allow the reading of these codons by cp tRNA(Leu)(UAm7G).

Rice mitochondrial genome contains a rearranged chloroplast gene cluster

We have previously reported the isolation and partial sequence analysis of a rice mitochondrial DNA fragment (6.9 kb) which contains a transferred copy of a chloroplast gene cluster coding for the large subunit of ribulose-1,5-bisphosphate carboxylase (rbcL), beta and epsilon subunits of ATPase (atpB and atpE), methionine tRNA (trnM) and valine tRNA (trnV). We have now completely sequenced this 6.9 kb fragment and found it to also contain a sequence homologous to the chloroplast gene coding for the ribosomal protein L2 (rpl2), beginning at a site 430 bp downstream from the termination codon of rbcL. In the chloroplast genome, two copies of rpl2 are located at distances of 20 kb and 40 kb, respectively, from rbcL. We have sequenced these two copies of rice chloroplast rpl2 and found their sequences to be identical. In addition, a 151 bp sequence located upstream of the chloroplast rpl2 coding region is also found in the 3' noncoding region of chloroplast rbcL and other as yet undefined locations in the rice chloroplast genome. Hybridization analysis revealed that this 151 bp repeat sequence identified in rice is also present in several copies in 11 other plant species we have examined. Findings from these studies suggest that the translocation of rpl2 to the rbcL gene cluster found in the rice mitochondrial genome might have occurred through homologous recombination between the 151 bp repeat sequence present in both rpl2 and rbcL.

A model for the evolution of the Vicia faba chloroplast genome

The Vicia faba chloroplast genome lacks inverted repeat sequences and contains only one set of ribosomal RNA genes. The genetic organization has been altered by inversions, relative to the typical arrangement of most higher plant chloroplast genomes. The Vicia faba plastid genome thus represents one of the more interesting results of chloroplast genomic evolution. The present study employs small DNA probes and Southern blot hybridizations to investigate the steps involved in the evolution of the Vicia faba chloroplast genome. The data from heterologous hybridizations between chloroplast DNA of Brassica napus (a conserved genome) and of Vicia faba led to three observations: 1) The inverted repeat segment closest to the psbA gene was deleted prior to the rearrangements. 2) A quarter of the ancestral small single copy region was lost during the deletion. 3) The genetic organization observed in Vicia faba resulted from three inversions after the deletion event. Our findings, combined with previous observations, helped devise a stepwise model for the evolution of the Vicia faba chloroplast genome. The area of the small single copy region absent from the Vicia faba chloroplast chromosome lacks in vivo transcription activity in Brassica napus.

Sequence organization of the chloroplast ribosomal spacer of Chlamydomonas reinhardii: uninterrupted tRNAile and tRNAala genes and extensive secondary structure

The 1805 bp spacer between the chloroplast ribosomal 16S and 7S RNA genes of Chlamydomonas reinhardii has been sequenced. It contains the genes of tRNA ala and tRNA ile which are both uninterrupted. The spacer includes several short direct and inverted repeats and a large palindromic structure which maps in the region where DNA rearrangements have occurred in other Chlamydomonas species.

Divergence of chloroplast gene organization in three legumes: Pisum sativum, Vicia faba and Phaseolus vulgaris

Isolated chloroplasts from Pisum sativum were found to contain at least 32 tRNA species. Hybridization of in vitro labeled, identified, chloroplast tRNAs to Pisum chloroplast DNA fragments revealed the locations of the tRNA genes on the circular chloroplast genome. Comparison of this gene map to the maps of Vicia faba and Phaseolus vulgaris showed that the chloroplast genomes of Pisum and Phaseolus are otherwise more closely related than either genome is to the chloroplast genome of Vicia. Furthermore, the results suggest how possible recombination events could be involved in the evolution of these three closely related, but divergent, chloroplast genomes.

Chloroplast transfer RNAs and tRNA genes of wheat

Fractionation (by two-dimensional polyacrylamide gel electrophoresis) of total tRNA from wheat chloroplasts yields about 33 RNA spots. Of these, 30 have been identified by aminoacylation as containing tRNAs specific for 17 amino acids.Hybridization of labeled individual tRNAs to cloned chloroplast DNA fragments has revealed the location of at least nine pairs of tRNA genes in the segments of the inverted repeat, at least twelve tRNA genes in the large single copy region and one tRNA gene in the small single copy region.A comparison of this wheat chloroplast tRNA gene map to that of maize and of other higher plants suggests that gene rearrangements have occurred during evolution, even within cereal chloroplast DNA. These rearrangements have taken place within the inverted repeat, within the large single copy region and between the inverted repeat and the large single copy region.

The intergenic region between the Vicia faba chloroplast tRNA(CAALeu) and tRNA(UAALeu) genes contains a partial copy of the split tRNA(UAALeu) gene

A cluster of three tRNA genes located on fragment Bam6a from Vicia faba chloroplast DNA has been sequenced: it contains the genes for tRNA(CAALeu), tRNA(UAALeu) and tRNA(Phe). The two tRNA(Leu) genes are separated by 443 bp and are transcribed divergently from different DNA strands. The intergenic region contains a series of short repeats and a partial copy of the split tRNA(UAALeu) gene which includes 100 bp of the 5' flanking region, 35 bp of the 5'exon and the first 42 bp of the intron. It is possible that some of these duplications occurred upon the rearrangement of the two tRNA(Leu) genes in broad bean (and in pea) or upon the deletion of one copy of the inverted repeat, since in all other higher plant chloroplast genomes studied so far these two tRNA(Leu) genes are located far apart on the genome, one being in the inverted repeat region, the other one in the large single copy region. The tRNA(Phe) and tRNA(UAALeu) are encoded by the same DNA strand, and separated by 110 bp.