Open Access

A survey of green plant tRNA 3'-end processing enzyme tRNase Zs, homologs of the candidate prostate cancer susceptibility protein ELAC2

BMC Evolutionary Biology201111:219

DOI: 10.1186/1471-2148-11-219

Received: 8 May 2011

Accepted: 23 July 2011

Published: 23 July 2011

Abstract

Background

tRNase Z removes the 3'-trailer sequences from precursor tRNAs, which is an essential step preceding the addition of the CCA sequence. tRNase Z exists in the short (tRNase ZS) and long (tRNase ZL) forms. Based on the sequence characteristics, they can be divided into two major types: bacterial-type tRNase ZS and eukaryotic-type tRNase ZL, and one minor type, Thermotoga maritima (TM)-type tRNase ZS. The number of tRNase Zs is highly variable, with the largest number being identified experimentally in the flowering plant Arabidopsis thaliana. It is unknown whether multiple tRNase Zs found in A. thaliana is common to the plant kingdom. Also unknown is the extent of sequence and structural conservation among tRNase Zs from the plant kingdom.

Results

We report the identification and analysis of candidate tRNase Zs in 27 fully sequenced genomes of green plants, the great majority of which are flowering plants. It appears that green plants contain multiple distinct tRNase Zs predicted to reside in different subcellular compartments. Furthermore, while the bacterial-type tRNase ZSs are present only in basal land plants and green algae, the TM-type tRNase ZSs are widespread in green plants. The protein sequences of the TM-type tRNase ZSs identified in green plants are similar to those of the bacterial-type tRNase ZSs but have distinct features, including the TM-type flexible arm, the variant catalytic HEAT and HST motifs, and a lack of the PxKxRN motif involved in CCA anti-determination (inhibition of tRNase Z activity by CCA), which prevents tRNase Z cleavage of mature tRNAs. Examination of flowering plant chloroplast tRNA genes reveals that many of these genes encode partial CCA sequences. Based on our results and previous studies, we predict that the plant TM-type tRNase ZSs may not recognize the CCA sequence as an anti-determinant.

Conclusions

Our findings substantially expand the current repertoire of the TM-type tRNase ZSs and hint at the possibility that these proteins may have been selected for their ability to process chloroplast pre-tRNAs with whole or partial CCA sequences. Our results also support the coevolution of tRNase Zs and tRNA 3'-trailer sequences in plants.

Keywords

tRNA precursor (pre-tRNA) endonuclease tRNase Z ELAC2 post-transcriptional processing green plant phylogenetics

Background

tRNA 3'-end maturation is a process through which the 3'-trailer sequence of precursor tRNAs (pre-tRNAs) is removed, and processed tRNAs acquire the CCA end which is absolutely essential for tRNA aminoacylation and protein synthesis (for reviews, see [13]). In prokaryotes, this process can be either exonucleolytic or endonucleolytic depending on whether the 3'-CCA sequence is genomically encoded. CCA-containing pre-tRNAs are generally processed by the exonucleases that tend to stop removing nucleotides from the 3'-end upon encountering the transcriptionally encoded CCA, whereas CCA-less pre-tRNAs are processed by a 3'-endonuclease termed tRNase Z (also termed RNase Z or 3'-tRNase; for reviews, see [47]) that cleaves immediately after the N73 discriminator nucleotide (the first unpaired base after the acceptor stem) to allow subsequent addition of the CCA sequence.

Unlike prokaryotic pre-tRNAs, eukaryotic nuclear and organellar pre-tRNAs generally lack the 3'-CCA sequence (which is added post-transcriptionally) and their 3'-trailer sequences are removed by tRNase Z. Also unlike prokaryotic pre-tRNAs, eukaryotic nuclear pre-tRNAs contain oligo (U) at their 3'-ends, which are recognized and bound by the La protein (for reviews see [2, 8]). In the budding yeast Saccharomyces cerevisiae and fission yeast Schizosaccharomyces pombe, the endonucleolytic cleavage of nuclear pre-tRNAs requires the presence of the yeast La protein [9, 10]. In the absence of the yeast La protein, the 3'-trailer sequence of nuclear pre-tRNAs is trimmed by 3'-exoribonucleases including Rex1p [11]. However, organellar pre-tRNAs lack terminal oligo (U). Furthermore, unlike nuclear pre-tRNAs which are typically monocistronic, most organellar pre-tRNAs are polycistronic [12, 13].

tRNase Z is present in all kingdoms of life. It exists in two forms: tRNase ZS [300-400 amino acids (aa)] and tRNase ZL (700-800 aa), which are encoded by different genes. It is believed that the tRNase ZL gene has evolved from a tandem duplication of the tRNase ZS gene, followed by divergence of the sequence [14]. In prokaryotes, only tRNase ZS is identified. By contrast, all eukaryotes possess tRNase ZL, and some have both forms.

The species distribution of tRNase Z is complex. The majority of eukaryotic species analyzed to date, including S. cerevisiae, the fruit fly Drosophila melanogaster and the nematode worm Caenorhabditis elegans contain a single tRNase ZL [1517]. In contrast, S. pombe have two tRNase ZLs [18, 19]. Interestingly, two tRNase ZLs and two tRNase ZSs have been experimentally identified in the flowering plant Arabidopsis thaliana [20]. In humans, one tRNase ZS (also termed ELAC1) and one tRNase ZL (also termed ELAC2) are found [14]. Our BLAST searches against public genomic and expressed sequence tag (EST) databases reveal that with few exceptions, vertebrates contain one tRNase ZL and one tRNase ZS (a detailed description of tRNase Z protein distribution in the animal kingdom will be provided elsewhere).

tRNase Z belongs to the metallo-β-lactamase (MBL) superfamily [14, 2124]. The typical MBL domain contains five conserved sequence motifs termed Motifs I-V. Motifs I and IV each harbor an invariant Asp, Motif II (HxHxDH), which is also called the His motif, is the signature motif of the superfamily, whereas Motifs III and V each contain a conserved His residue. Structural studies of tRNase ZSs from E. coli, T. maritima and B. subtilis [2528] and mutation analyses of tRNase Zs from a variety of species [2935] reveal that the His and Asp residues of Motifs II-V form the active site for coordination of two catalytic zinc ions. In particular, the Asp residue of Motif II may participate in both zinc ion coordination and act as a general base to generate a hydroxide ion for nucleophilic attack on the scissile phosphodiester bond at the cleavage site [25, 29]. The Asp residue of Motif I is also catalytically important and appears to stabilize the catalytic site [33].

Besides tRNase Zs, some nucleic acid processing enzymes are also members of the MBL superfamily. Most of these proteins belong to the β-CASP (MBL-associated CISF Artemis SNM1/PSO2) subfamily of the MBL [23]. This subfamily includes the 73-kD subunit of the cleavage and polyadenylation specificity factor (CPSF-73) and its yeast homolog Ysh1p, which are involved in endonucleolytic cleavage of pre-mRNA, the Intergrator complex subunit 11 (Int11) involved in the 3'-end formation of small nuclear RNAs (snRNA) [36], bacterial RNase J, which participates in rRNA 5'-end maturation [37] and RNA decay [38], and the eukaryotic Pso/Snm1/Artemis proteins, which function in DNA repair and V(D)J recombination [23]. However, unlike tRNase Zs, β-CASP proteins contain conserved β-CASP sequence motifs in place of Motif V.

tRNase Z is distinguished from other MBL members by their unique substrate binding domain termed the flexible arm (also termed the exosite). Based on flexible arm type, there are two major types (bacterial- and eukaryotic-types) and one minor type [T. maritima (TM)-type] of tRNase Zs [39]. The bacterial-type tRNase Zs, which are present predominantly in bacteria, possess the bacterial-type flexible arm. The bacterial-type flexible arm is ~55 aa in length and contains the Gly- and Pro-rich GP motif (GxPxGP, sometimes GxPPGP) [39]. The eukaryotic-type tRNase Zs, which are found only in eukaryotes, contain the ELAC2-type flexible arm. This type of flexible arm harbors the GP motif and is ~62 aa long, which is slightly longer than the bacterial-type flexible arm.

The TM-type tRNase Z was believed to be the minor type at the time of discovery since it was found only in T. maritima and A. thaliana [39]. The flexible arm found in TM-type tRNase Zs appears to be shorter (~30 aa) and lacks the GP motif but instead contains one short basic residue-rich region [39]. In addition, both the bacterial- and eukaryotic-type tRNase Zs contain the PxKxRN, HEAT and HST motifs, which form part of loop structures, whereas the TM-type tRNase Z lacks these motifs [33, 40, 41]. The PxKxRN motif has been suggested to function in CCA anti-determination (tRNase Z activity is inhibited by 3'-CCA) [25, 33], whereas the HEAT and HST motifs have been suggested to play a role in facilitating proton transfer at the final stage of reaction [25, 29, 40].

tRNase Z has diverse functions besides its primary role in tRNA 3'-end processing. This is perhaps best exemplified by ELAC2, which serves a multitude of functions within cells. Recent studies have shown that ELAC2 is involved in the generation of MALAT1, a cancer-associated long noncoding RNA which participates in regulation of pre-mRNA splicing [42], tRNA-derived small RNAs [43, 44], and viral microRNAs (miRNAs) [45, 46]. Overexpression of ELAC2 delays cell cycle progression, suggesting that ELAC2 may be involved in cell cycle control either directly or indirectly via its role as tRNA processing enzyme [47]. ELAC2 also potentiates TGF-β(transforming growth factor-β/Smad-induced transcription response, indicating a role for ELAC2 in TGF-β/Smad signaling mediated growth arrest [48]. Interestingly, a recent study has shown that destruction of human mitochondria through depletion of mitochondrial DNA results in down-regulation of ELAC2 and a delay in cell cycle progression [49]. Since ELAC2 may be involved in cell cycle regulation, it is likely that ELAC2 may link mitochondrial function and cell cycle control. It is important to note that ELAC2 is a candidate prostate cancer susceptibility gene as its mutations are associated with prostate cancer [14]. However, the underlying mechanisms are unknown. In S. cerevisiae, either inactivating mutations or overexpression of tRNase ZL causes a petite phenotype, suggesting that the action of tRNase ZL may be related to mitochondrial function [15]. In addition, the S. cerevisiae tRNase ZL has also been suggested to play a role in 35S rRNA processing [50].

The study of tRNase Z evolution has been facilitated by the increasing availability of genome sequences. A previous study showed that only tRNase ZS is found in bacteria and that its presence in bacteria is widespread [6]. We recently reported on a systematic survey of tRNase Zs in fungi [51]. Our analysis reveals that while the majority of fungal species contain one tRNase ZL, all four sequenced Schizosaccharomyces species contain two distinct tRNase ZLs either demonstrated or predicted to be localized to the nucleus and mitochondria, respectively. In addition, the presence of tRNase ZS in fungi is restricted to the phylum Basidiomycota and the basal fungal phyla.

Green plants (Viridiplantae) represent a monophyletic group of land plants and green algae that evolved near the base of the tree of eukaryotic life. Flowering plants (angiosperms), which are typically polyploidy, represent the largest, most diverse and most evolutionary advanced phylum of land plants making up 90% of the plant kingdom. It can be divided into two major groups: dicotyledons (dicots), which accounts for the majority of the angiosperm species, and monocotyledons (moncots). At present, there are at least 27 sequenced and annotated genomes representing the major taxonomic groups within green plants, although the majority of them are those of flowering plants. The public availability of these genome sequences enabled us to identify tRNase Zs in green plants and to study their evolution.

In this study, we undertook a comprehensive survey of candidate tRNase Zs from annotated green plant genomes. To understand the evolutionary relationships among green plant tRNase Zs, we further conducted a phylogenetic analysis of these newly identified candidates. Finally, we presented a detailed sequence analysis of tRNase Zs with the intent of further delineating the distinct features of green plant tRNase Zs.

Results

Identification of candidate green plant tRNase Zs

To extend our previous study of tRNase Z diversity and evolution, we searched public genome databases for putative green plant tRNase Zs with significant matches to known bacterial and eukaryotic tRNase Zs. Since most of candidate sequences identified from the databases are computationally generated without subsequent manual annotation, it is likely that many predictions may contain errors. Therefore, we verified each candidate. We first validated each prediction by reciprocal searches against the GenBank. In back-searches, a candidate was confirmed if reverse BLAST also gave tRNase Z hits in the top matches. Accuracy of prediction was further evaluated by multiple sequence alignment. All discordant candidate sequences were checked manually for possible errors including sequencing errors, intron mispredictions and existence of gaps in the genome sequences. We found that many candidate sequences are apparently incomplete or contain annotation errors. For example, the predicted coding sequence of the monkeyflower (Mimulus guttatus) MguTRZ2 (Phytozome accession no. mgv1a024577 m.g) in the database was incomplete lacking the N-terminal region. We were able to predict this region from the genomic DNA sequence based on sequence similarity. The predicted full-length coding region of MguTRZ2 has 364 aa. As another example, the sequence annotated as the candidate castor bean (Ricinus communis) tRNase ZL (Phytozome accession no. 30146.t000117) appears to be mispredicted due to the presence of sequence gaps. Thus, this sequence was excluded from the list.

Several incorrect predictions are apparently caused by the presence of the non-canonic GC-AG splice site pairs. While the GT donor splice site is a canonical 5'-splice site for introns in eukaryotic genes, the GC donor splice sites account for the majority of the non-canonical donor splice site. Thus, for those that cannot be accurately predicted by the conventional FGENESH program, we carried out gene prediction using FGENESH _GC, which is a new version of the FGENESH program including noncanonical GC dinucleotide in donor splice sites. Indeed, by doing so, we could predict some exons encoding missing conserved motifs. For example, the HEAT motif was originally missing in several annotated candidate tRNase ZLs from flowering plants including papaya (Carica papaya CpaTRZ3), cassava (Manihot esculenta MesTRZ4), Medicago (Medicago truncatula MtrTRZ3) and black cottonwood (Populus trichocarpa PtrTRZ3). After re-evaluation of intron splice sites using FGENESH_GC, we were able to recover their HEAT motifs. Because the intron sequences of some candidate sequences appear to be extremely difficult to predict correctly, we could not conclusively rule out the possibility of errors in certain candidates.

In total, 54 candidate tRNase ZS and 32 candidate tRNase ZL were identified from 27 green plant species including 21 flowering plants, 1 moss, 1 lycophyte and 4 green algae (Additional file 1). The names of these candidates follow the A. thaliana tRNase Z nomenclature [52]. Of these, only tRNase Zs from A. thaliana have been experimentally characterized [20]. While most of these sequenced genomes examined belong to the flowering plants, the availability of the genome sequences from two basal land plants (the bryophyte Physcomitrella patens and the lycophyte Selaginella moellendorffii) and four green algae (Chlamydomonas reinhardtii, Volvox carteri, Micromonas pusilla and Ostreococcus lucimarinus) allows for evaluation of differences between flowering plant tRNase Zs and those from basal land plants and green algae.

The flowering plants examined to date appear to contain multiple tRNase Zs. Foxtail millet (Setaria italica) contains the largest number of tRNase Zs (5) so far identified in a flowering plant. Most flowering plant species have a single tRNase ZL. In contrast, six flowering plants including two Arabidopsis species (A. thaliana and Arabidopsis lyrata) harbor two tRNase ZLs. It should be noted that the presence of two tRNase ZLs are not species-specific since the plant species containing two tRNase ZLs come from diverse taxonomic groups.

Unlike most flowering plants that possess two tRNase ZSs, two members of the Panicoideae subfamily of grasses, sorghum (Sorghum bicolor) and foxtail millet have three tRNase ZSs and four tRNase ZSs, respectively. The presence of multiple tRNase ZSs appears not to be grass-specific, since three other grass plants including rice (Oryza sativa), Brachypodium (Brachypodium distachyon) and maize (Zea mays) contain only two tRNase ZSs. Calculation of the percentage identity and similarity between candidate tRNase Zs from these two grass species shows strong conservation of the proteins at the amino acid level with the most identity (93%) and similarity (95%) between S. bicolor SbiTRZ2 and S. italica SitTRZ2 (Additional files 2 and 3). To our surprise, all candidate tRNase ZSs identified in flowering plants are highly similar and belong to the TM-type tRNase ZS (see below for a detailed discussion).

The number and type of tRNase Zs appear to be highly variable in the two primitive plant species and four green algae (Table 1 and Additional file 1). The largest number of tRNase Zs is found in the moss P. patens, which contains two TM-type and one bacterial-type tRNase ZSs, and one tRNase ZL. In contrast, two green algae C. reinhardtii and V. carteria have the least number of tRNase Zs, comprising one TM-type tRNase ZS and one tRNase ZL. The two other green algae M. pusilla and O. lucimarinus contain one TM-type and one bacterial-type tRNase ZSs and one tRNase ZL. The lycophyte S. moellendorffii, which has the smallest genome size of any land plant reported, contains two tRNase ZLs in addition to one TM-type tRNase ZS.
Table 1

Distribution of candidate tRNase Zs from representative green plants

Species

Protein name

   Type

Accession number

   Database

   No. aa+

DICOTS

     

Aquilegia coerulea

AcoTRZ1

   TM-type tRNase ZS

AcoGoldSmith_v1.019418 m.g

   Phytozome

   283

Aquilegia coerulea

AcoTRZ2

   TM-type tRNase ZS

AcoGoldSmith_v1.007495 m.g

   Phytozome

   356

Aquilegia coerulea

AcoTRZ3

   tRNase ZL

AcoGoldSmith_v1.000749 m.g

   Phytozome

   997*

Aquilegia coerulea

AcoTRZ4

   tRNase ZL

AcoGoldSmith_v1.024294 m.g

   Phytozome

   848*

Arabidopsis thaliana

AthTRZ1

   TM-type tRNase ZS

NP_177608.2

   NCBI

   280

Arabidopsis thaliana

AthTRZ2

   TM-type tRNase ZS

NP_178532.2

   NCBI

   354

Arabidopsis thaliana

AthTRZ3

   tRNase ZL

NP_175628.2

   NCBI

   890

Arabidopsis thaliana

AthTRZ4

   tRNase ZL

NP_188247.2

   NCBI

   942

Carica papaya

CpaTRZ1

   TM-type tRNase ZS

evm.TU.supercontig_17.217

   Phytozome

   306

Carica papaya

CpaTRZ2

   TM-type tRNase ZS

evm.TU.supercontig_132.51

   Phytozome

   359*

Carica papaya

CpaTRZ3

   tRNase ZL

evm.TU.supercontig_34.31

   Phytozome

   943*

Citrus clementina

CclTRZ1

   TM-type tRNase ZS

clementine0.9_017657 m.g

   Phytozome

   288

Citrus clementina

CclTRZ2

   TM-type tRNase ZS

clementine0.9_014182 m.g

   Phytozome

   355

Citrus clementina

CclTRZ3

   tRNase ZL

clementine0.9_001700 m.g

   Phytozome

   938

Cucumis sativus

CsaTRZ1

   TM-type tRNase ZS

Cucsa.369470

   Phytozome

   301*

Cucumis sativus

CsaTRZ2

   TM-type tRNase ZS

Cucsa.185180

   Phytozome

   347

Cucumis sativus

CsaTRZ3

   tRNase ZL

Cucsa.359360

   Phytozome

   977

Eucalyptus grandis

EgrTRZ1

   TM-type tRNase ZS

Egrandis_v1_0.020509 m.g

   Phytozome

   309

Eucalyptus grandis

EgrTRZ2

   TM-type tRNase ZS

Egrandis_v1_0.021229 m.g

   Phytozome

   357*

Eucalyptus grandis

EgrTRZ3

   tRNase ZL

Egrandis_v1_0.001658 m.g

   Phytozome

   971

Glycine max

GmaTRZ1

   TM-type tRNase ZS

Glyma20g01480

   Phytozome

   279

Glycine max

GmaTRZ2

   TM-type tRNase ZS

ACU18735.1

   NCBI

   353

Glycine max

GmaTRZ3

   tRNase ZL

Glyma13g43270

   Phytozome

   923*

Glycine max

GmaTRZ4

   tRNase ZL

Glyma15g02070

   Phytozome

   931*

Mimulus guttatus

MguTRZ1

   TM-type tRNase ZS

mgv1a010642 m.g

   Phytozome

   306

Mimulus guttatus

MguTRZ2

   TM-type tRNase ZS

mgv1a024577 m.g

   Phytozome

   364*

Mimulus guttatus

MguTRZ3

   tRNase ZL

mgv1a000815 m.g

   Phytozome

   976

Vitis vinifera

VviTRZ1

   TM-type tRNase Z S

XP_002273058.1

   NCBI

   290

Vitis vinifera

VviTRZ2

   TM-type tRNase Z S

XP_002277241.1

   NCBI

   354

Vitis vinifera

VviTRZ3

   tRNase ZL

XP_002278956.1

   NCBI

   951

MONOCOTS

     

Oryza sativa japonica

OsaTRZ1

   TM-type tRNase ZS

NP_001046280.1

   NCBI

   302

Oryza sativa japonica

OsaTRZ2

   TM-type tRNase ZS

LOC_Os09g30466

   Phytozome

   365

Oryza sativa japonica

OsaTRZ3

   tRNase ZL

B9EUI7

   UniProt

   964

Zea mays

ZmaTRZ1

   TM-type tRNase ZS

GRMZM2G147727

   Phytozome

   302

Zea mays

ZmaTRZ2

   TM-type tRNase ZS

ACG47847.1

   NCBI

   357

Zea mays

ZmaTRZ3

   tRNase ZL

GRMZM2G379286

   Phytozome

   930

MOSSES

     

Physcomitrella patens

PpaTRZ1

   TM-type tRNase ZS

XP_001785449.1

   NCBI

   295

Physcomitrella patens

PpaTRZ2

   TM-type tRNase ZS

Pp1s17_61V6

   Phytozome

   420

Physcomitrella patens

PpaTRZ3

   Bacterial-type tRNase ZS

Pp1s409_46V6

   Phytozome

   336*

Physcomitrella patens

PpaTRZ4

   tRNase ZL

Pp1s126_39V6

   Phytozome

   984*

Physcomitrella patens

PpaTLP1

   tRNase ZL-like

Pp1s3_568V6

   Phytozome

   923*

LYCOPHYTES

     

Selaginella moellendorffii

SmoTRZ1

   TM-type tRNase ZS

167717

   Phytozome

   299

Selaginella moellendorffii

SmoTRZ2

   tRNase ZL

111409

   Phytozome

   771*

Selaginella moellendorffii

SmoTRZ3

   tRNase ZL

97243

   Phytozome

   819*

Selaginella moellendorffii

SmoTLP1

   tRNase ZS-like

404484

   Phytozome

   303

Selaginella moellendorffii

SmoTLP2

   tRNase-like

XP_002967989

   NCBI

   346

Selaginella moellendorffii

SmoTLP3

   tRNase ZL-like

416540

   Phytozome

   731

GREEN ALGAE

     

Chlamydomonas reinhardtii

CreTRZ1

   TM-type tRNase ZS

Cre10.g431400

   Phytozome

   335

Chlamydomonas reinhardtii

CreTRZ2

   tRNase ZL

Cre01.g039900

   Phytozome

   880*

Chlamydomonas reinhardtii

CreTLP1

   tRNase ZS-like

Cre12.g492250

   Phytozome

   425*

Chlamydomonas reinhardtii

CreTLP2

   tRNase ZS -like

Cre07.g321900

   Phytozome

   378

Chlamydomonas reinhardtii

CreTLP3

   tRNase ZL-like

Cre01.g068350

   Phytozome

   710*

Chlamydomonas reinhardtii

CreTLP4

   tRNase ZL-like

Cre09.g394700

   Phytozome

   1073

Micromonas pusilla

MpuTRZ1

   TM-type tRNase ZS

5692

   JGI

   383*

Micromonas pusilla

MpuTRZ2

   Bacterial-type tRNase ZS

18700

   JGI

   398*

Micromonas pusilla

MpuTRZ3

   tRNase ZL

54013

   JGI

   953*

Micromonas pusilla

MpuTLP1

   tRNase ZL-like

C1N4P0

   UniProt

   866

Ostreococcus lucimarinus

OluTRZ1

   TM-type tRNase ZS

26345

   JGI

   313

Ostreococcus lucimarinus

OluTRZ2

   Bacterial-type tRNase ZS

XP_001416263

   NCBI

   369*

Ostreococcus lucimarinus

OluTRZ3

   tRNase ZL

A4S9E8

   UniProt

   781*

Ostreococcus lucimarinus

OluTLP1

   tRNase ZL-like

A4RR41

   UniProt

   780

Volvox carteri

VcaTRZ1

   TM-type tRNase Z S

80918

   Phytozome

   308*

Volvox carteri

VcaTRZ2

   tRNase ZL

118548

   Phytozome

   881*

Volvox carteri

VcaTLP1

   tRNase ZS-like

103834

   Phytozome

   445*

Volvox carteri

VcaTLP2

   tRNase ZL-like

100314

   Phytozome

   723*

+ The number of amino acids in plant tRNase Z and tRNase Z-like proteins.

*Indicates that mispredicted sequences obtained from the databases have been corrected.

Interestingly, tRNase Z-like proteins (TLP) which apparently lack one or more conserved motifs of tRNase Zs necessary for the enzymatic activity of the protein are widespread in basal land plant and green algal genomes that have been analyzed here (Table 1 and Additional file 1). Many tRNase Z-like proteins appear to lack the flexible arm but contain all other conserved motifs of tRNase Zs (Figure 1 and see below for a more detailed discussion of the motifs). At one extreme, one S. moellendorffii tRNase Z-like protein (SmoTLP2) appears to lack all of the conserved motifs.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig1_HTML.jpg
Figure 1

Alignment of bacterial-type and tRNase Z S -like candidates in green plants. Bacterial-type tRNase ZSs are from P. patens (PpaTRZ3), O. lucimarinus (OluTRZ2), M. pusilla (MpuTRZ2), Synechocystis sp. PCC 6803 (slr0050), E. coli (EcoTrz1) [72], B. subtilis (BsuTrz1) [60] and humans (HsaTrz1) [14]. Plant tRNase ZS-like (TLP) candidates are from S. moellendorffii (SmoTLP1 and SmoTLP2), V. carteri (VcaTLP1), C. reinhardtii (CreTLP1 and CreTLP2). TM-type tRNase ZSs from T. maritima (TmaTrz1) [30] and A. thaliana (AthTRZ1) [20] are included for comparison. The alignment was constructed using Clustal W [68]. Identical residues are on a black ground and conserved residues shared in gray. Also indicated above the alignment are the conserved motifs of tRNase Zs, which are labeled according to references [33, 40, 41]. The numbers in brackets indicate the length of the region in the protein, which are species-specific and could not be correctly aligned. Hyphens represent gaps introduced into sequences for maximum alignment. Amino acid residues predicted to be critical for activity are indicated by a star.

Based on their sizes, tRNase Z-like proteins can be divided into tRNase ZS-like and tRNase ZL-like proteins which are comparable in size to tRNase ZSs and tRNase ZLs, respectively. The number and form of tRNase Z-like proteins vary among the species, being largest in C. reinhardtii (two tRNase ZS-like and two tRNase ZL-like proteins) and next largest in S. moellendorffii (two tRNase ZS-like and one tRNase ZL-like proteins). The basal land plant P. patens and the two green algae M. pusilla and O. lucimarinus contain one tRNase ZL-like protein, whereas the green alga V. carteria contains one tRNase ZS-like and one tRNase ZL-like proteins. In contrast, black cottonwood (P. trichocarpa) appears to be the only species among the flowering plant genomes examined that contains the tRNase ZL-like protein. This species contains one tRNase ZL-like protein in which the second His in the His motif is mutated to Gln (data not shown).

Prediction of subcellular localization of candidate tRNase Zs from flowering plants

To help understand the function of candidate tRNase Zs, we predicted in silico the subcellular localization of each of flowering plant tRNase Zs using different bioinformatic prediction programs. Most flowering plants have two tRNase ZSs (tRNase ZS1 and tRNase ZS2). All tRNase ZS1 proteins apparently lack any predictable signal sequences and therefore predicted to be cytoplasmic proteins, with the exception of O. sativa tRNase ZS1 (OsaTRZ1), which is predicted to be localized in the chloroplasts (Table 2). On the other hand, all tRNase ZS2 except P. trichocarpa tRNase ZS2 (PtrTRZ2), contain a predicted chloroplast targeting signal (Table 2). The predicted chloroplast targeting signal of PtrTRZ2 is unusual in that it is only seven residues long since the majority of the known chloroplast targeting signals consist of N-terminal 20-80 amino acid residues [53]. In the two flowering plants (S. bicolor and S. italica) carrying multiple tRNase ZSs, only one tRNase ZS (tRNase ZS2) from each species is predicted to be chloroplast-localized, whereas the rest are predicted to be cytoplasmic (Table 2).
Table 2

Prediction of the chloroplast targeting signals in candidate TM-type tRNase ZSs from flowering plants

Protein name

Form

Chl

Chloroplast targeting signal

AcoTrz1

AcoTrzS1

-

 

AcoTrz2

AcoTrzS2

C

MKISLGISTSPILLPFSKPTLKIPHNRVCIP (31)

AlyTrz1

AlyTrzS1

-

 

AlyTrz2

AlyTrzS2

C

MQLSSFPISPSKIFPFTKYHAPPVIIHQLAAQIQSHRSNFVSPVKVSGYFSSISRAIEEEEEYRKAR (67)

AthTrz1

AthTrzS1

-

 

AthTrz2

AthTrzS2

C

MQLSSSFPISPPKIFPSTKHHKPPVITHQLAAQIQSNRRHFVSPVKVSGYFSSISRAIEEEEEYRKAR (68)

BdiTrz1

BdiTrzS1

-

 

BdiTrz2

BdiTrzS2

C

MATISLFSLPPLRLTRGLLSPSSGPASRFQTL (32)

CpaTrz1

CpaTrzS1

-

 

CpaTrz2

CpaTrzS2

C

MLTSVPILPSKLQSVFPFGHHHYHSISHKPTPIHQFSFLV (40)

CclTrz1

CclTrzS1

-

 

CclTrz2

CclTrzS2

C

MQLSLPTSPSKLPTIFPFHPSSIPKTPQSPHHLSLQSHVGPLNALKSAGFLSSISRAIDEEEEYRKAR (68)

CsiTrz1

CsiTrz S1

-

 

CsiTrz2

CsiTrz S2

C

MQLSLPTSPSKLPTIFPFHPSSIPKTPQSPHHLSLQSHVGPLNALKSAGFLSSISRAIDEEEEYRKAR (68)

CsaTrz1

CsaTrzS1

-

 

CsaTrz2

CsaTrzS2

C

MQISIPISPLRPPQVFPFHQPLLHTPKPPGVALQSHLNPVNSFRDSGLLSTIGVE (55)

EgrTrz1

EgrTrzS1

-

 

EgrTrz2

EgrTrzS2

C

MRISIPISPSKSPQIFPFHHHHRRRHQIPDLRRRPPALSLPASAVNPLKSSGYLSTIHR (59)

GmaTrz1

GmaTrzS1

-

 

GmaTrz2

GmaTrzS2

C

MQISLSDLAFKTPQLFPIHHPIFPPKPPLNHQVST (35)

MesTrz1

MesTrzS1

-

 

MesTrz2

MesTrzS2

C

MQTSFPVSPSKIPSIFPFNHPILHNPTAPTNHHQLPLQTYLKRPINSTSLKSSGFLSAIGRAIEEEEEYKKAR (73)

MtrTrz1

MtrTrzS1

-

 

MtrTrz2

MtrTrzS2

C

MQIPLTNPTFKPPQIFPFHHTIPSPKQQLSFPT (33)

MguTrz1

MguTrzS1

-

 

MguTrz2

MguTrzS2

C

MITLNPTTSNLLNLHPFHPTPATHKHHLPPQTRRSVIACVSDAVV (45)

OsaTrz1

OsaTrzS1

C

MANSGKSSPAATSTTAPPPGRPKAKAPPLTVEGYPVEGISIGGQETCVIFPT (52)

OsaTrz2

OsaTrzS2

C

MAATSLLSLPSLRLTHRLLVPASSSAPASRSQFQTL (36)

PtrTrz1

PtrTrzS1

-

 

PtrTrz2

PtrTrzS2

C

MQISIPL (7)

PpeTrz1

PpeTrzS1

-

 

PpeTrz2

PpeTrzS2

C

MQIFLTISPCKAPLILPFHPPISKTPKTQRTSLKTLSKPYRLARSQVLRKGV (52)

RcoTrz1

RcoTrzS1

-

 

RcoTrz2

RcoTrzS2

C

MQTSLPISHSKFPSIFPFNHPISHKPTTTR (30)

SitTrz1

SitTrzS1

-

 

SitTrz2

SitTrzS2

C

MATASLFSPPSLRLLSRTTARLSRFQTLAARKPP (34)

SitTrz3

SitTrzS3

-

 

SitTrz4

SitTrzS4

-

 

SbiTrz1

SbiTrzS1

-

 

SbiTrz2

SbiTrzS2

C

MATASLFSLPSLRALSRTSARCSRFQTLAARKPVESSSSTATS (43)

SbiTrz3

SbiTrzS3

-

 

VviTrz1

VviTrzS1

-

 

VviTrz2

VviTrzS2

C

MQISLPFSTSKVPYLSPLPNPTPQPPLTIPKPHPKSYIT (39)

ZmaTrz1

ZmaTrzS1

-

 

ZmaTrz2

ZmaTrzS2

C

MVTASLFSLPSLRVLSRTSAHGPRFQILAARKPVESSTATSGSRRGGSKGAGLLSVLDR (59)

Identification of the putative chloroplast targeting signals was made by using ChloroP http://www.cbs.dtu.dk/services/ChloroP/. Chl stands for the chloroplast. C indicates chloroplast localization. "-", the localization could not be predicted. The numbers refer to amino acid positions starting from the N-terminus.

Besides tRNase ZS, flowering plants also contain one (tRNase ZL1) or two (tRNase ZL1 and tRNase ZL2) tRNase ZLs. Most tRNase ZL1 proteins contain a predicted mitochondrial targeting signal between two putative translational initiation sites at the N-terminus and also have a nuclear localization signal (Table 3). Some tRNase ZL2 proteins have a putative mitochondrial targeting signal, while others are predicted to have both nuclear and mitochondrial targeting signals. The length of the putative mitochondrial targeting signals found in these candidates are within the expected size-range (20-80 aa) [54]. Consistent with the prediction, A. thaliana tRNase ZL1 (AthTRZ3) was found in both the nucleus and the mitochondria [20]. However, A. thaliana tRNase ZL2 (AthTRZ4) predicted to have both nuclear and mitochondrial targeting signals is actually found only in the mitochondria [20].
Table 3

Prediction of the nuclear and mitochondrial targeting sequences in tRNase ZLs from flowering plants

Protein name

From

Nuc

Mito

Predicted mitochondrial targeting sequences

AcoTrz3

AcoTrzL1

N

-

 

AcoTrz4

AcoTrzL2

-

M

MPQISSNLKLFFSKTNHSPLFQFSFKASFCSSFLLSSKTPYKPLSSVTVISSSSSSSSSSRKGPKFPPLRSRS (73)

AlyTrz3

AlyTrzL1

N

M

MINSMPYLHKNLRLLLLLSSKSSPFPLSLRPFSPRSFSLSTLFS (44)

AlyTrz4

AlyTrzL2

N

M

MLTSSMPHRHVPQNLSLFGFSPLKSSSFALFLRPFSLYPPIFASSSPSPSRRPPRTAGYRRS (62)

AthTrz3

AthTrzL1

N

M

MINSMPYLHKNLRLLRLLSSKSSPFPLSLRPFSPRSFSLSTLFS (44)

AthTrz4

AthTrzL2

N

M

MLTSSMPQNLSLFGFSPLKSSSFALILRPFSLYPPIFASSSPAPSRRPPRTAGYRRSGPSPPRRK (65)

BdiTrz3

BdiTrzL1

-

M

MPQVAAPLRLLLPLSQTLAPPAPLLHLSRRLLSFCSPASFRRAASLRALAYRRSRHPEPRRG (62)

CpaTrz3

CpaTrzL1

N

M

MFLIYPNLRLLLNPPLFLFSKPNSTPLSLFTVFASSSHKRHRSVSYRDSPFGLHRRRRNFTTFKERD (67)

CclTrz3

CclTrzL1

N

M

MPFITPNLRLLFSSSSSSLFPLKLSVPLLSTKPTNRHRSLFTILSYSKRQRSTPFPQQNQRRN (63)

CsiTrz3

CsiTrzL1

N

M

MPFITPNLRLLFSSSSSSSLFPLKLSVPLLSTKPTNRHHSLFTILSYSKRQRSTPFPQQNQRRN(64)

CsaTrz3

CsaTrzL1

N

M

MPLPHLSTLRFLFFSPSKLPFSPSLYSPKSHSLFTVLASSPPKRRRSATAPPSLNFKRRN (60)

EgrTrz3

EgrTrzL1

N

M

MPCVYSNLRLLFSSSATAAAATATAAASPFLSPLKLKLRRPSSSSSAFPLLPLPPLSSLRS (61)

GmaTrz3

GmaTrzL1

N

M

MAQVSKFGYFLLHSSLPKPSNIQFRSLLTVLASSSKRHRRK (41)

GmaTrz4

GmaTrzL2

N

M

MAQVSKFGHLLLHSSLPKPSNSNIQFRSLLTLLASSSKRHRSIPPFRRKS (50)

MesTrz3

MesTrzL1

N

M

MPQISNLRFLLSPIKPSLPFPFSKPKPYSLFTVLCSSSSSRRHRTTPNHQSLNFRSRS (58)

MesTrz4

MesTrzL2

-

M

MPNVLNFKLCLSDLLTQCCKICLQHYILTFSSLSHVLVSGKLQSSNPQTYFSINNPLPRSKNSRT (65)

MtrTrz3

MtrTrzL1

N

M

MAQILNFRNFLFLPSYKPTTHFRLRFLSTLVSSSSRRSNINAPPLHLRRRS (51)

MguTrz3

MguTrzL1

N

M

MPQSTNLRLLLSSANCHRRHPFSAASNFFPKHLSFSSSFQFFLKPQFKTREIPLLFATFSSYSKKPYATNNNSNNNNKNSRSFNRN (86)

OsaTrz3

OsaTrzL1

-

M

MPQLPSPLRRLLPLSQTLAAATPAPLLHLSRRLFSSSSSPSPSPSPRAACLRALAYRGGQAGGGGRRGHHNNLLRRG (77)

PtrTrz3

PtrTrzL1

N

M

MSHISNLRLLLSPLNPTLRFPFSSKHRPYSLLTILSSSSPYPKRRHRTTPNHPSLNFRSRS (61)

PpeTrz3

PpeTrzL1

N

M

MPQVTNLRLLFFSPFPRLSLSSLSFKPLKPRTLFTALASSYRKRHRPIPNQSPNTGARN (59)

SitTrz5

SitTrzL1

-

M

MPQVAAPLRLLLPLSQTLAPAPLLHLSRRLFTSSSPSFGRAASLRALAYRRHHHPRRGSSTLRK (64)

SbiTrz4

SbiTrzL1

N

M

MPQVAGPLRRLLPLSQTLASAPAPLLHLSRRLLSSCSPASFGRAASLRALAYRRRRHPEPRRG (63)

VviTrz3

VviTrzL1

N

M

MPHLTSFRLLYCSPLLSPFKSPFLSFSTLSKSKSPLLNPPSFFTVLSSSSGRYPKLRRHPHHLRRRN (67)

ZmaTrz3

ZmaTrzL1

N

M

MPQVAAPLRLLLPLSQTLAAPAPLLHHSRRLLSSCPASLSRAAGLRALAYRRRRHTEPRRG (61)

The putative nuclear localization signals were predicted using PSORT http://psort.hgc.jp/form.html, while the putative mitochondrial targeting signals were predicted using MITOPROT http://ihg2.helmholtz-muenchen.de/ihg/mitoprot.html. Nuc: stand for the nucleus; Mito: mitochondria. N indicates nuclear localization and M denotes mitochondrial localization. "-", the localization could not be predicted. The numbers refer to amino acid positions starting from the N-terminus.

Phylogenetic analysis

To gain insights into the evolutionary relationship among plant tRNase Zs, we reconstructed the phylogeny of 86 candidate sequences using Bayesian phylogenetics. Although most of plant species are flowering plants, they are taxonomically diverse. In addition to tRNase Zs from plants, we also include tRNase ZS from T. maritima as outgroup. Phylogenetic analysis reveals the presences of two well-supported clades: one formed by including all TM-type tRNase ZSs, and the other all tRNase ZLs and bacterial-type tRNase ZSs (Figure 2). The latter clade is further divided into two subclades, of which one contains all tRNase ZLs, and the other bacterial-type tRNase ZSs. Notably, two tRNase ZLs (AlyTRZ4 and AthTRZ4) from the two Arabidopsis species (A. thaliana and A. lyrata) together form a group sister to a group formed by another two tRNase ZLs (AlyTRZ3 and AthTRZ3) from the same two species. Likewise, the two tRNase ZLs found in some species including Aquilegia coerulea (AcoTRZ3 and AcoTRZ4), Glycine max (GmaTRZ3 and GmaTRZ4) and S. moellendorffii (SmoTRZ2 and SmoTRZ3) are sister to each other with a posterior probability value of 1. These observations suggest that tRNase ZL gene duplication took place in certain species.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig2_HTML.jpg
Figure 2

Phylogram showing phylogenetic relationships among candidate green plant tRNase Zs. Posterior probabilities, which are indicated at the nodes, are generated by the Bayesian analysis. The scale bar indicates 0.1 nucleotide substitutions per site. For each protein, the species, the accession number and the database can be found in Additional file 1. Taxonomic designations are indicated on the right side of the tree.

Conservation of candidate green plant tRNase ZSs

To assess the presence and conservation of motifs in candidate green plant tRNase Zs, we performed multiple sequence alignments of identified tRNase Zs sequences. Candidate tRNase ZSs are analyzed first. A list of aligned tRNase ZSs from representative green plants is shown in Figures 1 and 3. For comparison, tRNase ZSs from T. maritima, B. subtilis, E. coli, the cyanobacterium Synechocystis sp. PCC 6803 and humans are included as needed. A full list of all aligned green plant tRNase ZSs is presented in Additional file 4.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig3_HTML.jpg
Figure 3

Amino acid sequence alignment of the representative TM-type green plant and non-green plant tRNase Z S s. Plant TM-type tRNase ZSs are from C. reinhardtii (CreTRZ1), V. carteri (VcaTRZ1), M. pusilla (MpuTRZ1) O. lucimarinus (OluTRZ1), P. patens (PpaTRZ1 and PpaTRZ2), S. moellendorffii (SmoTRZ1), O. sativa japonica (OsaTRZ1 and OsaTRZ2), Z. mays (ZmaTRZ1 and ZmaTRZ2), A. thaliana (AthTRZ1 and AthTRZ2) [20], C. sativus (CsaTRZ1 and CsaTRZ2), G. max (GmaTRZ1 and GmaTRZ2), C. clementina (CclTRZ1 and CclTRZ2), M. guttatus (MguTRZ1 and MguTRZ2), A. coerulea (AcoTRZ1), V. vinifera (VviTRZ1 and VviTRZ2) and E. grandis (EgrTRZ1 and EgrTRZ2). Non-green plant tRNase ZSs are from T. maritima (TmaTRZ1) [30], B. subtilis (BsuTRZ1) [60], humans (HsaTRZ1) [14]. Protein accession numbers are shown in Table 1. Alignment presentations are as described in the legend to Figure 1.

Sequence comparison reveals that except for the N-terminal regions, tRNase ZSs from the flowering plants appear to be more similar to each other than to either the basal land plants or green algae (Figure 3). Furthermore, the sequences of the green plant tRNase ZSs are highly divergent from those of T. maritima, B. subtilis, E. coli and human tRNase ZSs. For example, A. thaliana tRNase ZS1 (AthTRZ1) and tRNase ZS2 (AthTRZ2) exhibit only 17% - 24% and 16% - 21% identity (25% - 35% and 27% - 31% similarity), respectively, with those from T. maritima, B. subtilis, E. coli and humans. Overall, the sequence conservation between the green plant tRNase ZSs and non-plant species is largely confined to the highly conserved motifs of tRNase Zs.

Unexpectedly, careful examination of the sequences of green plant tRNase ZSs reveals that most of the proteins possess several unique features that distinguish them from bacterial-type tRNase ZSs and thus justify their classification as TM-type tRNase ZS. The sequence logos for motifs unique to TM-type tRNase ZSs are presented in Figure 4. First, most green plant tRNase ZSs harbor the TM type flexible arm. The plant TM-type flexible arms show only weak protein sequence homology to the bacterial-type flexible arms, and have distinctive features including the absence of the GP motif and the presence of a consensus sequence KLKxxYxxLxGxxIxxLK, here termed the KL motif (Figures 3 and 4). This Lys- and Leu-rich motif was previously unappreciated in A. thaliana, likely due to the limited number of plant sequences available at the time analysis was performed [39]. As might be expected, there are variations in the consensus sequence (Figures 3 and 4).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig4_HTML.jpg
Figure 4

Sequence logos of unique motifs found in candidate TM-type tRNase Z S s from green plants. The sequence logos of motifs specific to candidate TM-type tRNase ZSs from green plants were derived from the 51 alignments. The sequence logos were created using WebLogo http://weblogo.berkeley.edu. The height of each amino acid indicates the level of conservation at that position. Amino acids are colored as follows: red, basic; blue, acidic; orange; polar; and green, hydrophobic.

Second, most green plant tRNase ZSs lack the PxKxRN motif normally present on the amino side of Motif I. Instead, they share a Gly-rich consensus sequence EGxSxxG in an analogous position to the PxKxRN motif (Figure 3). In some cases, variants in the consensus sequence are found (Figures 3 and 4). Notably, there is very weak sequence similarity between this motif and the corresponding region in T. maritima tRNase ZS.

Third, most green plant tRNase ZSs contain conserved xExT and HxH motifs in place of the HEAT and HST motifs, respectively. The xExT motif is related to the HEAT motif except that the invariant His is replaced by a Thr/Met/Val/Leu residue and the invariant Ala is often replaced by Ser, or occasionally Met, Thr and Cys (Figures 3 and 4). In contrast, the HxH motif is related to the HST motif, except that the conserved Ser is mostly replaced by Thr, and the conserved Thr is substituted by His (Figures 3 and 4). Based on the structural and mutagenesis studies, it has been suggested that the Glu of the HEAT motif and His of the HST motif play a role in facilitating proton transfer at the final stage of reaction [25, 29, 40]. In particular, the His residue is implicated as the proton donor. It is highly likely that the conserved Glu-His pair from the xExT and HxH motifs in the TM-type tRNase ZSs may also participate in the terminal proton transfer reaction.

Multiple sequence alignment of bacterial-type tRNase ZSs from the basal plant P. patens (PpaTRZ3) and green algae reveals that in addition to the well conserved Motifs I-V, these candidates possess a bacterial-type flexible arm containing the GP, HEAT and HST motifs, and the variant PxKxRN motif, in which Lys is replaced with Leu (Figure 1).

Conservation of tRNase ZLs in green plants

Since tRNase ZL can be divided into the N-terminal and C-terminal halves, which are related by weak sequence similarity, we aligned these two halves separately. The alignment of the N-terminal and C-terminal halves of representative candidate tRNase ZLs from diverse green plant species are shown in Figures 5 and 6, respectively. Two non-plant eukaryotic tRNase ZLs from D. melanogaster and humans were included for comparison. A complete list of all aligned green plant tRNase ZLs identified here is given in Additional file 5.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig5_HTML.jpg
Figure 5

Sequence alignment of N-terminal halves of representative green plant and non-green plant tRNase Z L s. Green plant tRNase ZLs are from O. sativa japonica (OsaTRZ3), Z. mays (ZmaTRZ3), A. thaliana (AthTRZ3 and AthTRZ4) [20], M. guttatus (MguTRZ3), G. max (GmaTRZ3), C. sativus (CsaTRZ3), C. clementina (CclTRZ3), C. papaya (CpaTRZ3), E. grandis (EgrTRZ3), V. vinifera (VviTRZ3), A. coerulea (AcoTRZ3 and AcoTRZ4), S. moellendorffii (SmoTRZ2 and SmoTRZ3), P. patens (PpaTRZ4), O. lucimarinus (OluTRZ3), M. pusilla (MpuTRZ3), Chlamydomonas reinhardtii (CreTRZ2) and V. carteri (VcaTRZ2). Non-green plant D. melanogaster (DmeTRZ1) [16] and humans (HsaTRZ2) [14] are included for comparison. The annotation of the alignment is as described in the legend to Figure 1.

https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig6_HTML.jpg
Figure 6

Sequence alignment of C-terminal halves of representative green plant and non-green plant tRNase Z L s. Same legend as in Figure 5.

Except for the N-terminal regions, candidate tRNase ZLs from land plants share a high degree of overall sequence similarity to each other. In contrast, candidate tRNase ZLs from the four green algae show relatively low overall sequence similarity to each other and to the land plant tRNase ZLs. In particular, their N-terminal halves are highly divergent compared to those from the land plants. Furthermore, sequence similarity between land plant tRNase ZLs and their homologs in green algae, D. melanogaster and humans is limited to the conserved motifs identified in tRNase Zs.

The N-terminal halves of all green plant tRNase ZLs contain a conserved eukaryotic-type flexible arm domain containing the GP motif. In addition, the N-terminal halves of the proteins also contain a pseudo-PxKxRN motif and pseudo-Motifs I and II, which differ from their original motifs in many positions including critical residues for tRNase Z functions (Figure 5). These pseudo-motifs were previously identified in candidate fungi tRNase ZLs [51]. In contrast, the C-terminal halves of green plant tRNase ZLs harbor conserved Motifs I-V and the PxKxRN, HEAT and HST motifs (Figure 6). These motifs appear in the same relative order in all the sequences in which they are present.

Many chloroplast tRNA genes encode a portion of the 3'-CCA sequence

It has been suggested that many chloroplast tRNA genes encode partial CCA sequences [55]. However, this conclusion was primarily based on the examination of the 3'-flanking sequences of tRNA genes from three chloroplast genomes. To determine if this conclusion holds true when more chloroplast genome sequences are now available for analysis, we examined the presence of tRNA genes encoding whole or partial CCA sequences in 15 chloroplast genomes. The results obtained were tabulated in Table 4.
Table 4

The distribution of the trinucleotide sequences immediately after the discriminator nucleotide in chloroplast tRNA genes

Species

74CCA76

74CCN76

74CNN76

Total genes examined

%+

Arabidopsis thaliana

0

1

11

37

32

Brachypodium distachyon

1

1

8

38

26

Carica papaya

0

1

8

37

24

Citrus sinensis

2

1

7

37

27

Cucumis sativus

0

1

7

37

22

Eucalyptus grandis

0

1

8

36

25

Glycine max

0

1

9

37

27

Manihot esculenta

0

1

7

37

22

Medicago truncatula

0

0

5

29

17

Oryza sativa japonica

1

2

7

38

26

Populus trichocarpa

0

0

10

37

27

Prunus persica

0

1

8

37

24

Sorghum bicolor

0

2

7

38

24

Vitis vinifera

0

1

8

37

24

Zea mays

1

2

7

37

27

Nucleotide sequences of chloroplast tRNA genes and their 3' flanking regions were obtained from the GenBank genome database http://www.ncbi.nlm.nih.gov/genomes/genlist.cgi?taxid=2759&type=4&name=Eukaryotae%20Organelles, the Chloroplast Genome Database (ChloroplastDB; http://chloroplast.cbio.psu.edu/organism.cgi), and the Organelle Genome Database (GOBASE; http://gobase.bcm.umontreal.ca/searches/gene.php). The discriminator nucleotide was predicted using tRNAscan-SE 1.21 http://lowelab.ucsc.edu/tRNAscan-SE/. +Percentage of chloroplast tRNA genes encoding 74CNN76, 74CCN76 and 74CCA76.

Indeed, we found that many tRNA genes in the chloroplast genomes encode partial CCA sequences, which could serve as part of the CCA sequence (Table 4). For instance, in the A. thaliana chloroplast genome, 31% of tRNA genes encode the whole or partial CCA sequences. Moreover, many chloroplast tRNA genes encode the first base of the CCA sequence. For example, of the 37 A. thaliana chloroplast tRNA genes examined, 11 had C after the discriminator.

Discussion

The presence of multiple tRNase Zs in green plants

Unlike C. elegans, D. melanogaster, humans and most fungal species examined to date, green plants are unique in that they possess multiple tRNase Zs. Of the 27 complete green plant genomes analyzed, the majority of them encode two tRNase ZSs and one or two tRNase ZLs. It is reasonable to expect that the existence of multiple tRNase Zs would be common to green plants. This phenomenon appears to have arisen from genome doubling (polyploidy), which seems to be a driving force in plant evolution and variation [56]. However, there seems to be no correlation between the number of tRNase Zs and genome size (Table 5).
Table 5

The genome size and ploidy level of representative green plants and the number of tRNase Zs encoded by their genomes

Species

Common name

No. tRNase Z#

Ploidy level

Genome size (Mbp)

Refs

Arabidopsis thaliana

Mouse-ear cress

4

Diploid

125

[73]

Brachypodium distachyon

Purple false brome

3

Diploid

272

[74]

Carica papaya

Papaya

3

Diploid

372

[75]

Cucumis sativus

Cucumber

3

Diploid

367

[76]

Glycine max

Soybean

4

Diploid

1100

[77]

Oryza sativa japonica

Rice

3

Diploid

389

[78]

Populus trichocarpa

Black cottonwood or poplar

3

Diploid

485

[79]

Ricinus communis

Castor bean

3

Diploid

350

[80]

Setaria italica

Foxtail millet

5

Diploid

490

[81]

Sorghum bicolor

Sorghum

4

Diploid

730

[82]

Vitis vinifera

Grape

3

Diploid

505

[83]

Zea mays

Maize

3

Diploid

2300

[84]

The genome size and ploidy level were obtained from the Plant DNA C-values database http://www.kew.org/cvalues/.

#The number of tRNase Zs

An unexpected observation in this study is that most green plant tRNase ZSs, which are clearly distinct from the bacterial-type tRNase ZS, represent the TM-type tRNase ZS. A comparison of motifs found in three different types of tRNase Zs from green plants is provided in Figure 7. Our phylogenetic analysis also supports the existence of both the TM- and bacterial-types of tRNase ZSs in green plants (Figure 2). The TM-type tRNase ZS was previously found only in the hyperthermophilic bacterium T. maritima and the flowering plant A. thaliana, and thus, it was originally thought to be a minor type. Our data greatly expand the repertoire of this type of tRNase Z. Although it remains to be determined if the TM-type tRNase ZS is also widespread in other taxonomic groups, the prevalence of the TM-type tRNase ZS in green plants suggests that this type might be plant-specific.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-219/MediaObjects/12862_2011_Article_1825_Fig7_HTML.jpg
Figure 7

Schematic representation of the three different types of tRNase Zs identified in green plants. The N-terminal half (eukaryotic-type-N) and C-terminal half (eukaryotic-type-C) of the eukaryotic-type tRNase ZL are shown separately. Motifs are indicated by colored boxes: red, specific to TM-type; green, specific to both bacterial- and eukaryotic-type; and blue, common to all three types. N and C denote N and C terminus, respectively. Relative positions of these motifs are not drawn to scale. The predicted functions of each motif are also indicated.

In contrast to the TM-type tRNase ZS, which is widespread in the green plants, the bacterial-type tRNase ZS is only present in the basal land plants and green algae. Sequence analysis reveals that the two types of tRNase ZSs are related to but distinct from each other, suggesting that they likely arose from the same ancestral tRNase ZS gene, followed by sequence divergence. Since the flowering plants possess only the TM-type tRNase ZS, the bacterial-type tRNase ZS has apparently been lost during the course of plant evolution whereas the TM-type tRNase ZS has been maintained. Interestingly, the bacterial-type tRNase ZS is also found in humans. Although the precise function of human tRNase ZS is unknown, it may play a role in degradation of a subset of microRNAs (miRNA) in the cytoplasm [57].

It is likely that multiple tRNase Zs found in the flowering plants are localized in different subcellular compartments. In most of the flowering plants examined, the two TM-type tRNase ZSs (tRNase ZS1 and tRNase ZS2) are predicted to be either in the cytosol or in the chloroplast, whereas one of tRNase ZL (tRNase ZL1) is predicted to contain both nuclear and mitochondrial targeting signals. The second tRNase ZL (tRNase ZL2) found in some plants either contains or lacks a predicted mitochondrial targeting signal. These predictions are generally supported by recent subcellular localization studies of A. thaliana tRNase Zs [20]. A. thaliana tRNase ZS1 (AthTRZ1) and tRNase ZS2 (AthTRZ2) are localized in the cytoplasm and chloroplast, respectively, whereas A. thaliana tRNase ZL1 (AthTRZ3) and tRNase ZL2 (AthTRZ4) are targeted to both the nucleus and mitochondria, and the mitochondria, respectively [20].

Based on their predicted subcellular localization, it is most likely that tRNase ZS2 and tRNase ZL1 are involved in the 3'-ends processing of pre-tRNAs in the chloroplasts, and both the nucleus and mitochondria, respectively, whereas tRNase ZL2 either serves as a back-up for mitochondrial pre-tRNA 3'-end processing or plays a specialized, yet to be discovered function. The function of tRNase ZS1 is unknown. It may play a role in repair of incorrectly processed 3'-ends of tRNAs in the cytosol. Alternatively, tRNase ZS1 could also participate in nuclear pre-tRNA 3'-end processing, as suggested by the observation that deletion of both A. thaliana tRNase ZS1 (AthTRZ1) and tRNase ZL1 (AthTRZ3) results in a lethal phenotype [20].

The plant TM-type tRNase ZSs may not recognize the CCA sequence as an anti-determinant

The CCA triplet following the discriminator generally inhibits tRNase Z activity, with the first C at nucleotide position 74 having the greatest effect [25, 33, 58]. This anti-determinant property of the CCA sequence is thought to prevent removal of the 3'-terminal CCA sequence from mature tRNAs, and thus avoids futile cycles of CCA addition and removal [33, 5860]. However, unlike other types of tRNase Zs, the TM-type tRNase Z does not appear to recognize the CCA sequence that is downstream the discriminator and present in the trailer sequence as an anti-determinant. The recombinant T. maritima tRNase ZS cleaves after the CCA sequence which is encoded by all but one tRNA gene, leaving the CCA sequence intact [30]. Similarly, the recombinant A. thaliana tRNase ZS can process pre-tRNAs with the 3'-terminal CCA sequences embedded in the trailer sequence, albeit at a different position [61]. Like the recombinant proteins, a partially purified spinach chloroplast tRNase Z fraction can also cleave pre-tRNAs containing complete or partial CCA sequences after the first C74 regardless of the sequence of the flanking region [55]. Importantly, the remainder of the CCA sequence can be added by chloroplast tRNA nucleotidyltransferease [55].

However, A. thaliana tRNase ZS appears to be unusual in that it can also cleave off the 3'-terminal CCA sequence from mature tRNAs in vitro [61]. Since mature tRNAs must be protected from counterproductive cleavage by tRNase ZS in vivo, it is likely that the plant chloroplast protein may acquire the ability to recognize the 3'-terminal CCA sequence as the mature tRNA 3'-end via cofactors [61].

Sequence determinants potentially involved in cleavage specificity

The mechanisms responsible for cleavage site selection and CCA inhibition are not fully understood, but appear to involve a combination of sequence features. First, the flexible arm has been suggested to contribute to cleavage site selection and the inhibitory effect of the CCA sequence [62]. This tRNase Z-specific element is located on the opposite side of the active site and binds primarily the D and T loops of the pre-tRNA [28, 39, 63]. Notably, the flexible arms of the T. maritima and plant tRNase ZSs lack the GP motif but contain the KL motif. Additionally, they are significantly smaller in length relative to those of other tRNase Zs possessing the anti-determinant function. Structural studies reveal that although TM- and bacterial-types of flexible arms share a similar overall structure which is composed of a compact globular domain and an extended two-stranded stalk and protrudes from the protein core, they have different globular domains [2528, 64]. In the TM-type flexible arm, the globular domain consists of one very short α-helix, one long helix and one 310-helix, whereas in the bacterial-type flexible arm, it is composed of two α-helices, two β-strands and one 310-helix. It would be interesting to know how the differences in the sequence feature, length and topology of the flexible arm may contribute to cleavage specificity by the enzymes.

Another possible motif involved in selection of the cleavage site has been suggested to be Motif I which appears to participate in binding the acceptor stem of pre-tRNA substrates [33, 62]. In vitro studies using the recombinant T. maritima tRNase ZS has suggested that the Ser31 and Thr33 residues of Motif I are involved in the cleavage site selection, with the former residue being more critical [30]. However, only individual, but not simultaneous mutations of these two non-Gln residues to Gln (which are found at the corresponding positions in other tRNase Zs with the CCA anti-determinant) in T. maritima tRNase ZS affects the cleavage site selection [62].

In addition to the flexible arm and Motif I, the PxKxRN motif, in particular, the two basic residues in the motif, has also been suggested to be involved in the cleavage specificity and, by inference, CCA inhibition of tRNase Zs, since this motif is absent from T. maritima and A. thaliana [33]. Strikingly, tRNase ZSs found in all flowering plants examined so far also lack this motif (Figure 3). However, since the archaeon Methanococcus jannaschii tRNase ZS harbors a PxKxRN motif but cannot recognize the CCA sequence as an anti-determinant in vitro, suggesting that other sequence elements may also be involved in the cleavage specificity of tRNase Z [61]. Alternatively, it has been suggested that a high enzyme concentration used in the assay may contribute to the lack of a CCA anti-determinant effect with the M. jannaschii tRNase ZS [33].

Why is the TM-type tRNase ZSdeveloped for chloroplast pre-tRNA 3'-end processing?

The discovery that the TM-type tRNase ZS is widespread in green plants raises a question as to why plants adopt the TM-type tRNase ZS over other types of tRNase Z for the 3'-end processing of chloroplast pre-tRNAs. One possible explanation is that the TM-type tRNase ZS identified in plants has evolved to adapt to chloroplast pre-tRNA 3'-end processing. It has previously been suggested that many plant chloroplast tRNA genes encode C74 based on a limited number of available chloroplast tRNA genes [55]. To see if this conclusion could be extended to more flowering plants, we have examined the 3'-flanking region for tRNA genes in additional chloroplast genomes. Indeed, many chloroplast tRNA genes seem to encode partial CCA sequences (Table 4).

The development of the TM-type tRNase ZS may be particularly important for chloroplasts, which have a limited but sufficient number of self-encoded tRNA species and do not import tRNAs [65, 66]. As all possible codons are used in the chloroplast protein-encoding genes, all chloroplast tRNA species appear to be used in protein synthesis [66]. Thus, efficient 3'-end processing of each pre-tRNA by tRNase Z may be critical in ensuring maximum efficiency in chloroplast protein synthesis.

The presence of candidate tRNase Z-like proteins in green plants

tRNase ZS-like proteins have previously been found in cyanobacteria including Synechocystis sp. PCC 6803 [67] and fungi [51], whereas tRNase ZL-like proteins have not been reported. In Synechocystis, one tRNase ZS-like protein (sll1036) has been identified. This protein does not exhibit any tRNase Z activity in vitro, consistent with the lack of some of the most conserved motifs of tRNase Zs in the protein sequence [67]. Since candidate tRNase Z-like proteins found in basal land plants and green algae either lack all or some of the conserved motifs of tRNase Zs, they most likely possess no tRNase Z activity.

It has been suggested that tRNase ZL has arisen from the fusion of duplicated tRNase ZS genes with further sequence diversification. It is possible that tRNase ZS-like proteins found in the basal land plants and green algae may represent relics of original tRNase ZS that were mutated during diversification of eukaryotic tRNase Z genes. Alternatively, these proteins may play species-specific functions as suggested for cyanobacterial tRNase ZS-like proteins [67]. The predominant presence in the basal land plants and green algae of tRNase Z-like proteins also suggests that duplication of ancestral tRNase ZS genes may occur early in green plant evolution.

Conclusions

This study represents the first large-scale identification and analysis of green plant tRNase Zs. Our survey of current plant genome databases shows that green plants are represented by multiple tRNase Zs, which include one or two tRNase ZLs and two tRNase ZSs. One tRNase ZL is predicted to participate in 3'-end processing of nuclear and mitochondrial pre-tRNA, whereas the other is likely to provide a backup for mitochondrial pre-tRNA processing. It appears that most tRNase ZSs, which is widespread throughout the green plants, belong to a minor but highly distinct type of tRNase ZS (TM-type). In contrast, the typical bacterial-type tRNase ZS is restricted to basal land plants and green algae. The apparent lack of the bacterial-type tRNase ZS in flowering plants suggest that while both types were present in the basal land plants, the bacterial-type tRNase ZS was discarded in favor of TM-type during plant evolution. Based on our results and previous studies, we propose that like T. maritima tRNase ZS, TM-type tRNase ZSs found in green plants seem not to recognize the CCA sequence as an anti-determinant and that the rise of this type of tRNase ZS appears to accommodate the 3'-end processing of chloroplast pre-tRNAs with partial or whole CCA sequences. This unusual property of green plant tRNase ZSs is likely due to multiple sequence determinants including the TM-type-specific flexible arm comprising the KL motif, Motif I and lack of the PxKxRN motif.

However, it should be noted that bioinformatics analysis alone cannot resolve possible differences in cleavage specificity among TM-type tRNase ZSs. The complete understanding of the mechanisms of the cleavage specificities of the TM-type tRNase ZSs awaits the determination of the structures of these enzymes bound to pre-tRNA and the accumulation of more precise biochemical data.

Methods

Plant genome database search and protein sequence analysis

Candidate tRNase Zs were identified by BLAST and PSI-BLAST searches against the genome databases using known tRNase Zs as query sequences. The databases used include Phytozome http://www.phytozome.net/, the NCBI nonredundant protein sequence database http://blast.ncbi.nlm.nih.gov/Blast.cgi, Joint Genome Institute (JGI; http://www.jgi.doe.gov/ and Universal Protein Resource (Uniprot; http://www.uniprot.org/. An E-value cutoff of 0.001 was used in all searches. With this value, no β-CASP protein or other MBL protein was found. The resulting sequences were subject to validation as described [51]. The splicing pattern was verified using the FGENESH and FGENESH_GC programs provided at the Softberry website http://linux1.softberry.com/berry.phtml?topic=fgenesh. Prediction of subcellular localization of proteins was made using web-based prediction programs such as MITOPROT http://ihg2.helmholtz-muenchen.de/ihg/mitoprot.html, PSORT http://psort.hgc.jp/form.html and ChloroP http://www.cbs.dtu.dk/services/ChloroP/. Multiple sequence alignments were done by Clustal W [68].

Phylogenetic analysis

Phylogenetic analysis was performed using the Bayesian approach, with T. maritima tRNase ZS (TmaTrz1) as an outgroup as described [51]. Briefly, full-length amino acid sequences of candidate plant tRNase Zs and TmaTrz1 were aligned by using Clustal W implemented in MEGA 5.0 [69]. After excluding gaps and the ambiguous sites, we used ProtTest 2.4 [70] to choose the most appropriate evolutionary model for our data set. The phylogenies were estimated by Bayesian inference with MrBayes 3.1.2 [71] using a mixture of the fixed amino acid models and I + G distribution. Statistical confidence was assessed by using Markov Chain Monte Carlo (MCMC) sampling approaches. Two separate runs including a total of four independent tree searches were conducted. All searches consisted of one 'cold' and three 'heated' Markov chains estimated for 107 generations, and every 1000 generations were sampled. The burn-in parameter was estimated by plotting -lnL against the generation number using Tracer 1.4.1 http://beast.bio.ed.ac.uk/Tracer, and the retained trees were used to estimate the consensus tree and the Bayesian posterior probabilities.

List of abbreviations

pre-tRNAs: 

tRNA precursors

miRNA: 

microRNAs

tRNase Z: 

tRNA 3' endonuclease

tRNase ZS

the short form of tRNase Z

tRNase ZL

the long form of tRNase Z

β-CASP: 

MBL-associated CISF Artemis SNM1/PSO2

CPSF-73: 

the cleavage and polyadenylation specificity factor

MBL: 

the metallo-β-lactamase superfamily

TGF-β

transforming growth factor-β

TM: 

T. maritima.

Declarations

Acknowledgements

We thank three anonymous reviewers for improving the quality of the manuscript, Jingjing Peng and Wei Zhao for valuable discussion. This work was supported in part by grants from the National Science Foundation of China (31070703), the Priority Academic Program Development of Jiangsu Higher Education Institutions and Nanjing Normal University (2007104XGQ0148).

Authors’ Affiliations

(1)
Laboratory of Yeast Genetics and Molecular Biology, School of Life Sciences, Nanjing Normal University
(2)
Jiangsu Key Laboratory for Microbes and Genomics, School of Life Sciences, Nanjing Normal University
(3)
Jiangsu Engineering and Technology Research Center for Industrialization of Microbial Resources, School of Life Sciences, Nanjing Normal University
(4)
Jiangsu Key Laboratory for Biodiversity and Biotechnology, School of Life Sciences, Nanjing Normal University

References

  1. Phizicky EM, Hopper AK: tRNA biology charges to the front. Genes Dev. 2010, 24: 1832-1860. 10.1101/gad.1956510.View ArticlePubMedPubMed CentralGoogle Scholar
  2. Bayfield MA, Yang R, Maraia RJ: Conserved and divergent features of the structure and function of La and La-related proteins (LARPs). Biochim Biophys Acta. 2010, 1799: 365-378.View ArticlePubMedPubMed CentralGoogle Scholar
  3. Maraia RJ, Lamichhane TN: 3' processing of eukaryotic precursor tRNAs. WIREs RNA. 2011, 2: 362-375. 10.1002/wrna.64.View ArticlePubMedPubMed CentralGoogle Scholar
  4. Ceballos M, Vioque A: tRNase Z. Protein Pept Lett. 2007, 14: 137-145. 10.2174/092986607779816050.View ArticlePubMedGoogle Scholar
  5. Hartmann RK, Gossringer M, Spath B, Fischer S, Marchfelder A: Chapter 8 The Making of tRNAs and More - RNase P and tRNase Z. Prog Nucleic Acid Res Mol Biol. 2009, 85C: 319-368.Google Scholar
  6. Redko Y, Li de Lasierra-Gallay I, Condon C: When all's zed and done: the structure and function of RNase Z in prokaryotes. Nat Rev Microbiol. 2007, 5: 278-286. 10.1038/nrmicro1622.View ArticlePubMedGoogle Scholar
  7. Spath B, Canino G, Marchfelder A: tRNase Z: the end is not in sight. Cell Mol Life Sci. 2007, 64: 2404-2412. 10.1007/s00018-007-7160-5.View ArticlePubMedGoogle Scholar
  8. Wolin SL, Cedervall T: The La protein. Annu Rev Biochem. 2002, 71: 375-403. 10.1146/annurev.biochem.71.090501.150003.View ArticlePubMedGoogle Scholar
  9. Van Horn DJ, Yoo CJ, Xue D, Shi H, Wolin SL: The La protein in Schizosaccharomyces pombe: a conserved yet dispensable phosphoprotein that functions in tRNA maturation. RNA. 1997, 3: 1434-1443.PubMedPubMed CentralGoogle Scholar
  10. Yoo CJ, Wolin SL: The yeast La protein is required for the 3' endonucleolytic cleavage that matures tRNA precursors. Cell. 1997, 89: 393-402. 10.1016/S0092-8674(00)80220-2.View ArticlePubMedGoogle Scholar
  11. Copela LA, Fernandez CF, Sherrer RL, Wolin SL: Competition between the Rex1 exonuclease and the La protein affects both Trf4p-mediated RNA quality control and pre-tRNA maturation. RNA. 2008, 14: 1214-1227. 10.1261/rna.1050408.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Schafer B: RNA maturation in mitochondria of S. cerevisiae and S. pombe. Gene. 2005, 354: 80-85.View ArticlePubMedGoogle Scholar
  13. Sugita M, Sugiura M: Regulation of gene expression in chloroplasts of higher plants. Plant Mol Biol. 1996, 32: 315-326. 10.1007/BF00039388.View ArticlePubMedGoogle Scholar
  14. Tavtigian SV, Simard J, Teng DH, Abtin V, Baumgard M, Beck A, Camp NJ, Carillo AR, Chen Y, Dayananth P, Desrochers M, Dumont M, Farnham JM, Frank D, Frye C, Ghaffari S, Gupte JS, Hu R, Iliev D, Janecki T, Kort EN, Laity KE, Leavitt A, Leblanc G, McArthur-Morrison J, Pederson A, Penn B, Peterson KT, Reid JE, Richards S, et al: A candidate prostate cancer susceptibility gene at chromosome 17p. Nat Genet. 2001, 27: 172-180. 10.1038/84808.View ArticlePubMedGoogle Scholar
  15. Chen Y, Beck A, Davenport C, Chen Y, Shattuck D, Tavtigian SV: Characterization of TRZ1, a yeast homolog of the human candidate prostate cancer susceptibility gene ELAC2 encoding tRNase Z. BMC Mol Biol. 2005, 6: 12-10.1186/1471-2199-6-12.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Dubrovsky EB, Dubrovskaya VA, Levinger L, Schiffer S, Marchfelder A: Drosophila RNase Z processes mitochondrial and nuclear pre-tRNA 3' ends in vivo. Nucleic Acids Res. 2004, 32: 255-262. 10.1093/nar/gkh182.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Smith MM, Levitan DJ: The Caenorhabditis elegans homolog of the putative prostate cancer susceptibility gene ELAC2, hoe-1, plays a role in germline proliferation. Dev Biol. 2004, 266: 151-160. 10.1016/j.ydbio.2003.10.016.View ArticlePubMedGoogle Scholar
  18. Zhao Z, Su W, Yuan S, Huang Y: Functional conservation of tRNase ZL among Saccharomyces cerevisiae, Schizosaccharomyces pombe and humans. Biochem J. 2009, 422: 483-492. 10.1042/BJ20090743.View ArticlePubMedGoogle Scholar
  19. Gan X, Yang J, Li J, Yu H, Dai H, Liu J, Huang Y: The fission yeast Schizosaccharomyces pombe has two distinct tRNase ZLs encoded by two different genes and differentially targeted to the nucleus and mitochondria. Biochem J. 2011, 435: 103-111. 10.1042/BJ20101619.View ArticlePubMedGoogle Scholar
  20. Canino G, Bocian E, Barbezier N, Echeverria M, Forner J, Binder S, Marchfelder A: Arabidopsis encodes four tRNase Z enzymes. Plant Physiol. 2009, 150: 1491-1502.View ArticleGoogle Scholar
  21. Aravind L: An evolutionary classification of the metallo-beta-lactamase fold proteins. In Silico Biol. 1999, 1: 69-91.PubMedGoogle Scholar
  22. Daiyasu H, Osaka K, Ishino Y, Toh H: Expansion of the zinc metallo-hydrolase family of the beta-lactamase fold. FEBS Lett. 2001, 503: 1-6. 10.1016/S0014-5793(01)02686-2.View ArticlePubMedGoogle Scholar
  23. Dominski Z: Nucleases of the metallo-beta-lactamase family and their role in DNA and RNA metabolism. Crit Rev Biochem Mol Biol. 2007, 42: 67-93. 10.1080/10409230701279118.View ArticlePubMedGoogle Scholar
  24. Schiffer S, Rosch S, Marchfelder A: Assigning a function to a conserved group of proteins: the tRNA 3'-processing enzymes. EMBO J. 2002, 21: 2769-2777. 10.1093/emboj/21.11.2769.View ArticlePubMedPubMed CentralGoogle Scholar
  25. de la Sierra-Gallay IL, Pellegrini O, Condon C: Structural basis for substrate binding, cleavage and allostery in the tRNA maturase RNase Z. Nature. 2005, 433: 657-661. 10.1038/nature03284.View ArticleGoogle Scholar
  26. Ishii R, Minagawa A, Takaku H, Takagi M, Nashimoto M, Yokoyama S: Crystal structure of the tRNA 3' processing endoribonuclease tRNase Z from Thermotoga maritima. J Biol Chem. 2005, 280: 14138-14144. 10.1074/jbc.M500355200.View ArticlePubMedGoogle Scholar
  27. Kostelecky B, Pohl E, Vogel A, Schilling O, Meyer-Klaucke W: The crystal structure of the zinc phosphodiesterase from Escherichia coli provides insight into function and cooperativity of tRNase Z-family proteins. J Bacteriol. 2006, 188: 1607-1614. 10.1128/JB.188.4.1607-1614.2006.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Li de la Sierra-Gallay I, Mathy N, Pellegrini O, Condon C: Structure of the ubiquitous 3' processing enzyme RNase Z bound to transfer RNA. Nat Struct Mol Biol. 2006, 13: 376-377. 10.1038/nsmb1066.View ArticlePubMedGoogle Scholar
  29. Minagawa A, Takaku H, Ishii R, Takagi M, Yokoyama S, Nashimoto M: Identification by Mn2+ rescue of two residues essential for the proton transfer of tRNase Z catalysis. Nucleic Acids Res. 2006, 34: 3811-3818. 10.1093/nar/gkl517.View ArticlePubMedPubMed CentralGoogle Scholar
  30. Minagawa A, Takaku H, Takagi M, Nashimoto M: A novel endonucleolytic mechanism to generate the CCA 3' termini of tRNA molecules in Thermotoga maritima. J Biol Chem. 2004, 279: 15688-15697. 10.1074/jbc.M313951200.View ArticlePubMedGoogle Scholar
  31. Schilling O, Vogel A, Kostelecky B, Natal da Luz H, Spemann D, Spath B, Marchfelder A, Troger W, Meyer-Klaucke W: Zinc- and iron-dependent cytosolic metallo-beta-lactamase domain proteins exhibit similar zinc-binding affinities, independent of an atypical glutamate at the metal-binding site. Biochem J. 2005, 385: 145-153. 10.1042/BJ20040773.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Vogel A, Schilling O, Meyer-Klaucke W: Identification of metal binding residues for the binuclear zinc phosphodiesterase reveals identical coordination as glyoxalase II. Biochemistry. 2004, 43: 10379-10386. 10.1021/bi049703+.View ArticlePubMedGoogle Scholar
  33. Zareen N, Hopkinson A, Levinger L: Residues in two homology blocks on the amino side of the tRNase Z His domain contribute unexpectedly to pre-tRNA 3' end processing. RNA. 2006, 12: 1104-1115. 10.1261/rna.4206.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Zareen N, Yan H, Hopkinson A, Levinger L: Residues in the conserved His domain of fruit fly tRNase Z that function in catalysis are not involved in substrate recognition or binding. J Mol Biol. 2005, 350: 189-199. 10.1016/j.jmb.2005.04.073.View ArticlePubMedGoogle Scholar
  35. Spath B, Kirchner S, Vogel A, Schubert S, Meinlschmidt P, Aymanns S, Nezzar J, Marchfelder A: Analysis of the functional modules of the tRNA 3' endonuclease (tRNase Z). J Biol Chem. 2005, 280: 35440-35447. 10.1074/jbc.M506418200.View ArticlePubMedGoogle Scholar
  36. Baillat D, Hakimi MA, Naar AM, Shilatifard A, Cooch N, Shiekhattar R: Integrator, a multiprotein mediator of small nuclear RNA processing, associates with the C-terminal repeat of RNA polymerase II. Cell. 2005, 123: 265-276. 10.1016/j.cell.2005.08.019.View ArticlePubMedGoogle Scholar
  37. Britton RA, Wen T, Schaefer L, Pellegrini O, Uicker WC, Mathy N, Tobin C, Daou R, Szyk J, Condon C: Maturation of the 5' end of Bacillus subtilis 16S rRNA by the essential ribonuclease YkqC/RNase J1. Mol Microbiol. 2007, 63: 127-138. 10.1111/j.1365-2958.2006.05499.x.View ArticlePubMedGoogle Scholar
  38. Mathy N, Benard L, Pellegrini O, Daou R, Wen T, Condon C: 5'-to-3' exoribonuclease activity in bacteria: role of RNase J1 in rRNA maturation and 5' stability of mRNA. Cell. 2007, 129: 681-692. 10.1016/j.cell.2007.02.051.View ArticlePubMedGoogle Scholar
  39. Schilling O, Spath B, Kostelecky B, Marchfelder A, Meyer-Klaucke W, Vogel A: Exosite modules guide substrate recognition in the ZiPD/ElaC protein family. J Biol Chem. 2005, 280: 17857-17862.View ArticlePubMedGoogle Scholar
  40. Karkashon S, Hopkinson A, Levinger L: tRNase Z catalysis and conserved residues on the carboxy side of the His cluster. Biochemistry. 2007, 46: 9380-9387. 10.1021/bi700578v.View ArticlePubMedPubMed CentralGoogle Scholar
  41. Levinger L, Hopkinson A, Desetty R, Wilson C: Effect of changes in the flexible arm on tRNase Z processing kinetics. J Biol Chem. 2009, 284: 15685-15691. 10.1074/jbc.M900745200.View ArticlePubMedPubMed CentralGoogle Scholar
  42. Wilusz JE, Freier SM, Spector DL: 3' end processing of a long nuclear-retained noncoding RNA yields a tRNA-like cytoplasmic RNA. Cell. 2008, 135: 919-932. 10.1016/j.cell.2008.10.012.View ArticlePubMedPubMed CentralGoogle Scholar
  43. Haussecker D, Huang Y, Lau A, Parameswaran P, Fire AZ, Kay MA: Human tRNA-derived small RNAs in the global regulation of RNA silencing. RNA. 2010, 16: 673-695. 10.1261/rna.2000810.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Lee YS, Shibata Y, Malhotra A, Dutta A: A novel class of small RNAs: tRNA-derived RNA fragments (tRFs). Genes Dev. 2009, 23: 2639-2649. 10.1101/gad.1837609.View ArticlePubMedPubMed CentralGoogle Scholar
  45. Bogerd HP, Karnowski HW, Cai X, Shin J, Pohlers M, Cullen BR: A mammalian herpesvirus uses noncanonical expression and processing mechanisms to generate viral MicroRNAs. Mol Cell. 2010, 37: 135-142. 10.1016/j.molcel.2009.12.016.View ArticlePubMedPubMed CentralGoogle Scholar
  46. Diebel KW, Smith AL, van Dyk LF: Mature and functional viral miRNAs transcribed from novel RNA polymerase III promoters. RNA. 2010, 16: 170-185. 10.1261/rna.1873910.View ArticlePubMedPubMed CentralGoogle Scholar
  47. Korver W, Guevara C, Chen Y, Neuteboom S, Bookstein R, Tavtigian S, Lees E: The product of the candidate prostate cancer susceptibility gene ELAC2 interacts with the gamma-tubulin complex. Int J Cancer. 2003, 104: 283-288. 10.1002/ijc.10945.View ArticlePubMedGoogle Scholar
  48. Noda D, Itoh S, Watanabe Y, Inamitsu M, Dennler S, Itoh F, Koike S, Danielpour D, Ten Dijke P, Kato M: ELAC2, a putative prostate cancer susceptibility gene product, potentiates TGF-beta/Smad-induced growth arrest of prostate cells. Oncogene. 2006, 25: 5591-5600. 10.1038/sj.onc.1209571.View ArticlePubMedGoogle Scholar
  49. Mineri R, Pavelka N, Fernandez-Vizarra E, Ricciardi-Castagnoli P, Zeviani M, Tiranti V: How do human cells react to the absence of mitochondrial DNA?. PloS one. 2009, 4: e5713-10.1371/journal.pone.0005713.View ArticlePubMedPubMed CentralGoogle Scholar
  50. Peng WT, Robinson MD, Mnaimneh S, Krogan NJ, Cagney G, Morris Q, Davierwala AP, Grigull J, Yang X, Zhang W, Mitsakakis N, Ryan OW, Datta N, Jojic V, Pal C, Canadien V, Richards D, Beattie B, Wu LF, Altschuler SJ, Roweis S, Frey BJ, Emili A, Greenblatt JF, Hughes TR: A panoramic view of yeast noncoding RNA processing. Cell. 2003, 113: 919-933. 10.1016/S0092-8674(03)00466-5.View ArticlePubMedGoogle Scholar
  51. Zhao W, Yu H, Li S, Huang Y: Identification and analysis of candidate fungal tRNA 3'-end processing endonucleases tRNase Zs, homologs of the putative prostate cancer susceptibility protein ELAC2. BMC Evol Biol. 2010, 10: 272-10.1186/1471-2148-10-272.View ArticlePubMedPubMed CentralGoogle Scholar
  52. Vogel A, Schilling O, Spath B, Marchfelder A: The tRNase Z family of proteins: physiological functions, substrate specificity and structural properties. Biol Chem. 2005, 386: 1253-1264. 10.1515/BC.2005.142.PubMedGoogle Scholar
  53. Jarvis P, Robinson C: Mechanisms of protein import and routing in chloroplasts. Curr Biol. 2004, 14: R1064-1077. 10.1016/j.cub.2004.11.049.View ArticlePubMedGoogle Scholar
  54. Pfanner N, Neupert W: The mitochondrial protein import apparatus. Annu Rev Biochem. 1990, 59: 331-353. 10.1146/annurev.bi.59.070190.001555.View ArticlePubMedGoogle Scholar
  55. Oommen A, Li XQ, Gegenheimer P: Cleavage specificity of chloroplast and nuclear tRNA 3'-processing nucleases. Mol Cell Biol. 1992, 12: 865-875.View ArticlePubMedPubMed CentralGoogle Scholar
  56. Adams KL, Wendel JF: Polyploidy and genome evolution in plants. Curr Opin Plant Biol. 2005, 8: 135-141. 10.1016/j.pbi.2005.01.001.View ArticlePubMedGoogle Scholar
  57. Takahashi M, Takaku H, Nashimoto M: Regulation of the human tRNase ZS gene expression. FEBS Lett. 2008, 582: 2532-2536. 10.1016/j.febslet.2008.06.020.View ArticlePubMedGoogle Scholar
  58. Mohan A, Whyte S, Wang X, Nashimoto M, Levinger L: The 3' end CCA of mature tRNA is an antideterminant for eukaryotic 3'-tRNase. RNA. 1999, 5: 245-256. 10.1017/S1355838299981256.View ArticlePubMedPubMed CentralGoogle Scholar
  59. Nashimoto M: Distribution of both lengths and 5' terminal nucleotides of mammalian pre-tRNA 3' trailers reflects properties of 3' processing endoribonuclease. Nucleic Acids Res. 1997, 25: 1148-1154. 10.1093/nar/25.6.1148.View ArticlePubMedPubMed CentralGoogle Scholar
  60. Pellegrini O, Nezzar J, Marchfelder A, Putzer H, Condon C: Endonucleolytic processing of CCA-less tRNA precursors by RNase Z in Bacillus subtilis. EMBO J. 2003, 22: 4534-4543. 10.1093/emboj/cdg435.View ArticlePubMedPubMed CentralGoogle Scholar
  61. Schiffer S, Rosch S, Marchfelder A: Recombinant RNase Z does not recognize CCA as part of the tRNA and its cleavage efficieny is influenced by acceptor stem length. Biol Chem. 2003, 384: 333-342. 10.1515/BC.2003.039.View ArticlePubMedGoogle Scholar
  62. Minagawa A, Ishii R, Takaku H, Yokoyama S, Nashimoto M: The flexible arm of tRNase Z is not essential for pre-tRNA binding but affects cleavage site selection. J Mol Biol. 2008, 381: 289-299. 10.1016/j.jmb.2008.05.016.View ArticlePubMedGoogle Scholar
  63. Hopkinson A, Levinger L: Effects of conserved D/T loop substitutions in the pre-tRNA substrate on tRNase Z catalysis. RNA Biol. 2008, 5: 104-111. 10.4161/rna.5.2.6086.View ArticlePubMedGoogle Scholar
  64. Ishii R, Minagawa A, Takaku H, Takagi M, Nashimoto M, Yokoyama S: The structure of the flexible arm of Thermotoga maritima tRNase Z differs from those of homologous enzymes. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2007, 63: 637-641. 10.1107/S1744309107033623.View ArticlePubMedPubMed CentralGoogle Scholar
  65. Michaud M, Cognat V, Duchene AM, Marechal-Drouard L: A global picture of tRNA genes in plant genomes. Plant J. 2010, 66: 80-93.View ArticleGoogle Scholar
  66. Pfitzinger H, Guillemaut P, Weil JH, Pillay DT: Adjustment of the tRNA population to the codon usage in chloroplasts. Nucleic Acids Res. 1987, 15: 1377-1386. 10.1093/nar/15.4.1377.View ArticlePubMedPubMed CentralGoogle Scholar
  67. Ceballos-Chavez M, Vioque A: Sequence-dependent Cleavage Site Selection by RNase Z from the Cyanobacterium Synechocystis sp. PCC 6803. J Biol Chem. 2005, 280: 33461-33469. 10.1074/jbc.M504691200.View ArticlePubMedGoogle Scholar
  68. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
  69. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011Google Scholar
  70. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.View ArticlePubMedGoogle Scholar
  71. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
  72. Ezraty B, Dahlgren B, Deutscher MP: The RNase Z homologue encoded by Escherichia coli elaC gene is RNase BN. J Biol Chem. 2005, 280: 16542-16545. 10.1074/jbc.C500098200.View ArticlePubMedGoogle Scholar
  73. Initiative TAG: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.View ArticleGoogle Scholar
  74. Initiative TIB: Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 463: 763-768.
  75. Ming R, Hou S, Feng Y, Yu Q, Dionne-Laporte A, Saw JH, Senin P, Wang W, Ly BV, Lewis KL, Salzberg SL, Feng L, Jones MR, Skelton RL, Murray JE, Chen C, Qian W, Shen J, Du P, Eustice M, Tong E, Tang H, Lyons E, Paull RE, Michael TP, Wall K, Rice DW, Albert H, Wang ML, Zhu YJ, et al: The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008, 452: 991-996. 10.1038/nature06856.View ArticlePubMedPubMed CentralGoogle Scholar
  76. Huang S, Li R, Zhang Z, Li L, Gu X, Fan W, Lucas WJ, Wang X, Xie B, Ni P, et al: The genome of the cucumber, Cucumis sativus L. Nat Genet. 2009, 41: 1275-1281. 10.1038/ng.475.View ArticlePubMedGoogle Scholar
  77. Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, et al: Genome sequence of the Palaeopolyploid soybean. Nature. 463: 178-183.
  78. Project IRGS: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.View ArticleGoogle Scholar
  79. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, et al: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313: 1596-1604. 10.1126/science.1128691.View ArticlePubMedGoogle Scholar
  80. Chan AP, Crabtree J, Zhao Q, Lorenzi H, Orvis J, Puiu D, Melake-Berhan A, Jones KM, Redman J, Chen G, Cahoon EB, Gedil M, Stanke M, Haas BJ, Wortman JR, Fraser-Liggett CM, Ravel J, Rabinowicz PD: Draft genome sequence of the oilseed species Ricinus communis. Nat Biotechnol. 28: 951-956.
  81. Doust AN, Kellogg EA, Devos KM, Bennetzen JL: Foxtail millet: a sequence-driven grass model system. Plant Physiol. 2009, 149: 137-141. 10.1104/pp.108.129627.View ArticlePubMedPubMed CentralGoogle Scholar
  82. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, et al: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457: 551-556. 10.1038/nature07723.View ArticlePubMedGoogle Scholar
  83. Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, Fitzgerald LM, Vezzulli S, Reid J, Malacarne G, Iliev D, Coppola G, Wardell B, Micheletti D, Macalma T, Facci M, Mitchell JT, Perazzolli M, Eldredge G, Gatto P, Oyzerski R, Moretto M, Gutin N, Stefanini M, Chen Y, Segala C, Davenport C, Dematte L, Mraz A, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PloS one. 2007, 2: e1326-10.1371/journal.pone.0001326.View ArticlePubMedPubMed CentralGoogle Scholar
  84. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA, Minx P, Reily AD, Courtney L, Kruchowski SS, Tomlinson C, Strong C, Delehaunty K, Fronick C, Courtney B, Rock SM, Belter E, Du F, Kim K, Abbott RM, Cotton M, Levy A, Marchetto P, Ochoa K, Jackson SM, Gillam B, et al: The B73 maize genome: complexity, diversity, and dynamics. Science. 326: 1112-1115.

Copyright

© Fan et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.