Identification and analysis of candidate fungal tRNA 3'-end processing endonucleases tRNase Zs, homologs of the putative prostate cancer susceptibility protein ELAC2
© Zhao et al; licensee BioMed Central Ltd. 2010
Received: 26 April 2010
Accepted: 6 September 2010
Published: 6 September 2010
tRNase Z is the endonuclease that is responsible for the 3'-end processing of tRNA precursors, a process essential for tRNA 3'-CCA addition and subsequent tRNA aminoacylation. Based on their sizes, tRNase Zs can be divided into the long (tRNase ZL) and short (tRNase ZS) forms. tRNase ZL is thought to have arisen from a tandem gene duplication of tRNase ZS with further sequence divergence. The species distribution of tRNase Z is complex. Fungi represent an evolutionarily diverse group of eukaryotes. The recent proliferation of fungal genome sequences provides an opportunity to explore the structural and functional diversity of eukaryotic tRNase Zs.
We report a survey and analysis of candidate tRNase Zs in 84 completed fungal genomes, spanning a broad diversity of fungi. We find that tRNase ZL is present in all fungi we have examined, whereas tRNase ZS exists only in the fungal phyla Basidiomycota, Chytridiomycota and Zygomycota. Furthermore, we find that unlike the Pezizomycotina and Saccharomycotina, which contain a single tRNase ZL, Schizosaccharomyces fission yeasts (Taphrinomycotina) contain two tRNase ZLs encoded by two different tRNase ZL genes. These two tRNase ZLs are most likely localized to the nucleus and mitochondria, respectively, suggesting partitioning of tRNase Z function between two different tRNase ZLs in fission yeasts. The fungal tRNase Z phylogeny suggests that tRNase ZSs are ancestral to tRNase ZLs. Additionally, the evolutionary relationship of fungal tRNase ZLs is generally consistent with known phylogenetic relationships among the fungal species and supports tRNase ZL gene duplication in certain fungal taxa, including Schizosaccharomyces fission yeasts. Analysis of tRNase Z protein sequences reveals putative atypical substrate binding domains in most fungal tRNase ZSs and in a subset of fungal tRNase ZLs. Finally, we demonstrate the presence of pseudo-substrate recognition and catalytic motifs at the N-terminal halves of tRNase ZLs.
This study describes the first comprehensive identification and sequence analysis of candidate fungal tRNase Zs. Our results support the proposal that tRNase ZL has evolved as a result of duplication and diversification of the tRNase ZS gene.
The endonuclease tRNase Z (also called RNase Z or 3'-tRNase) participates in maturation of tRNA 3'-end by removing the 3'-trailer sequence from tRNA precursors (pre-tRNAs, for reviews, see [1–4]). It belongs to the metallo-β-lactamase (MBL) superfamily, the members of which have diverse functions from hydrolysis and inactivation of β-lactam antibiotics to processing of RNA precursors [5–9]. Other nucleases in the MBL superfamily that act on nuclei acids include members of the β-CASP (MBL-associated CISF Artemis SNM1/PSO2) family: the cleavage and polyadenylation specificity factor 73 kDa subunit (CPSF-73)  and the Integrator complex subunit 11 (Int11) , which are involved in eukaryotic mRNA and small nuclear RNA (snRNA) 3'-end formation, respectively; RNase J, which functions in rRNA maturation and mRNA stability in bacteria  and the eukaryotic Pso2/Snm1/Artemis proteins, which play a role in DNA repair . Although displaying distinct substrate specificity defined by their specific domains, these proteins appear to have a similar catalytic mechanism since their active sites are composed of highly conserved motifs including the histidine motif (HxHxDH, where x is any hydrophobic residues).
There are two forms of tRNase Z. The long form (tRNase ZL) with 800-900 aa (amino acids) is about twice the size of the short form (tRNase ZS) with 300-400 aa . Sequence analysis suggests that tRNase ZL has evolved by gene duplication from tRNase ZS followed by sequence divergence . The species distribution of tRNase Z is not homogenous. tRNase ZS exists in all three domains of life (i.e. Bacteria, Archaea, and Eukarya) whereas tRNase ZL has been found only in eukaryotes so far. The number of tRNase Zs varies among different organisms. The largest number of tRNase Zs was detected in the plant Arabidopsis thaliana (two tRNase ZSs and two tRNase ZLs) . The fission yeast Schizosaccharomyces pombe contains two tRNase ZLs. Unexpectedly, the human genome encodes one tRNase ZS and one tRNase ZL. Human tRNase ZL gene (also termed ELAC2) was originally identified as the first prostate cancer susceptibility gene by positional cloning . However, the mechanism by which specific mutations in human tRNase ZL lead to an increased risk of prostate cancer remains unknown. In contrast, the budding yeast Saccharomyces cerevisiae, the fruit fly Drosophila melanogaster and the nematode worm Caenorhabditis elegans have just one tRNase ZL.
An intriguing question is why species have evolved to have more than one tRNase Z. One explanation is that additional tRNase Zs are targeted to organelles such as mitochondria and chloroplasts in which organelle-encoded pre-tRNAs must also be processed. Indeed, one of two S. pombe tRNase ZLs is targeted to the mitochondria, and has been suggested to play a role in mitochondrial-encoded pre-tRNA processing . In A. thaliana, three of four tRNase ZLs are targeted to organelles . However, it appears that the majority of tRNase ZLs identified thus far are imported both into the nucleus and mitochondria. Another explanation is that additional tRNase Zs may provide a back-up mechanism for nuclear tRNA 3'-end processing. The third explanation is that additional tRNase Zs may have different functions.
Recently, tRNase ZL itself has been either demonstrated or suggested to have additional functions other than tRNA 3'-end processing. For example, human tRNase ZL has been shown to play a role in generation of non-tRNA noncoding RNAs [15, 16] and viral microRNAs (miRNAs) . Moreover, human tRNase ZL has been proposed to cleave a subset of miRNAs in the cytoplasm . In S. cerevisiae, tRNase ZL has been suggested to have additional functions including rRNA biogenesis, mRNA splicing and mitochondrial maintenance . Similarly, our previous study also suggested that the nuclear-localized tRNase ZL in S. pombe may play a role beyond tRNA 3'-end processing .
Our current understanding of tRNase Z evolution is limited since there has been only one comprehensive survey on tRNase ZSs from prokaryotes . Of eukaryotes, the Fungi is a large and diverse kingdom encompassing roughly 1.5 million species and spanning one billion years of evolution . Sequence-based phylogenies show that the Chytridiomycota is the most basal phylum (group) among the Fungi, followed by the Zygomycota, with the Ascomycota and Basidiomycota as two largest phyla that together comprise the subkingdom Dikarya (also referred to as the "Higher Fungi") [21–24]. The Ascomycota (also known as sac fungi, yeasts or ascomycetes) is the largest and most diverse phylum in the Fungi, accounting for approximately 75% of all known fungi. Many popular model organisms such as S. cerevisiae, S. pombe, Neurospora crassa, Aspergillus nidulans and Candida albicans are classified in this phylum. This phylum is further divided into three major monophyletic subphyla (subgroups): Pezizomycotina, Saccharomycotina and Taphrinomycotina . The Pezizomycotina (also known as euascomyces) is the largest subphyla and contains over 90% of total Ascomycota species. They are multicellular filamentous fungi and grow by hyphal extension and branching. In contrast, the Saccharomycotina (also known as true yeasts) comprises the majority of unicellular species. The Taphrinomycotina is thought to be the earliest diverging group sister to the Saccharomycotina and Pezizomycotina. It constitutes a diverse group of organisms including unicellular yeast (for example, Schizosaccharomyces), multicellular filamentous fungi, and dimorphic fungi that can switch between yeast and hyphal growth forms. Like the Pezizomycotina, the Basidiomycota consists of primarily filamentous fungi.
Currently, most of eukaryotic species with sequenced genomes belong to the kingdom Fungi. The public fungal genome databases cover a broad range of fungal taxonomic groups with the majority coming from the Ascomycota and Basidiomycota. The availability of a large number of fungal genome sequences, together with the vast diversity of fungal morphology and lifestyle, provides an opportunity to identify tRNase Zs in the kingdom Fungi and to study evolution of eukaryotic tRNase Z.
In the present study, we performed a comprehensive survey of candidate tRNase Zs from 84 publicly available fungal genomes. To explore the evolutionary relationship among fungal tRNase Zs, we conducted a phylogenetic analysis of predicted fungal tRNase Zs. Finally, we examined their domain architectures. Our results support the view that tRNase ZL comes from tRNase ZS.
Identification of putative fungal tRNase Zs
As part of our efforts to better understand functional and structural diversity of tRNase Zs, we conducted extensive BLAST and PSI-BLAST homology searches against the publicly available fungal genome databases. Currently, the majority of sequenced species belong to the Dikarya, with a much higher proportion of Ascomycota species. Since other fungal phyla are poorly represented in public databases (three Zygomycota, three Microspordia, and three Chytridiomycota species), it is difficult to assess the true diversity of tRNase Z in these basal groups of fungi.
The initial candidates from the BLAST and PSI-BLAST were verified by multiple sequence alignment and reciprocal searches against the GenBank. Protein sequence alignment revealed a number of incorrectly predicted candidates, most likely due to misprediction of exon/intron boundaries or existence of gaps in the genome sequence. For example, the sequence (Broad accession no. CC1G_14814.2) annotated as the candidate Basidiomycete Coprinopsis cinerea tRNase Z in the fungal genome database at the Broad Institute has mispredicted exon/intron junctions. This 946-aa-long candidate is devoid of a histidine motif, which is a signature motif for the MBL superfamily, indicating that the exon encoding the histidine motif was likely mispredicted. After re-evaluating intron splicing pattern of the gene sequence, we were able to predict the exon encoding the histidine motif. The correctly predicted protein has 967 aa, and has the histidine motif. The sequence annotated as the candidate Pezizomycotina Botrytis cinerea (also named Botryotinia fuckeliana) tRNase Z (Broad accession no. BC1G_03733.1) is an example of misprediction due to the presence of sequence gaps in the genome. This sequence has 444 aa. However, examination of the genomic sequence revealed that its 5'-coding sequence contains gaps. Thus, this sequence was excluded. In back-searches, no candidate that shows homology to metallo-β-lactamase was found, but a limited number of candidates were found to show homology to the yeast homolog of CPSF-73 (Ysh1). Such candidates were also excluded from our final list. However, we cannot rule out the possibility that certain candidates may not be correctly predicted despite our efforts devoted to verification of these candidates.
We identified a total of 90 candidate tRNase ZLs and 19 candidate tRNase ZSs proteins from 84 fungal species including 67 Ascomycota, 14 Basidiomycota and 3 Chytridiomycota (Additional file 1). Candidate tRNase Zs from two taxonomic groups, the Zygomycota and Microspordia, were not listed since their full-length protein sequences could not be correctly predicted. Of the proteins identified here, only tRNase ZLs from S. cerevisiae and S. pombe have been experimentally characterized [14, 19, 25–27].
All species of the Ascomycota we have examined lack tRNase ZS. The Pezizomycotina and Saccharomycotina species have a single tRNase ZL. Surprisingly, in contrast to the Pezizomycotina and Saccharomycotina species, all four sequenced Schizosaccharomyces species (S. pombe, Schizosaccharomyces octosporus, Schizosaccharomyces japonicus and a recently described Schizosaccharomyces cryophobus) in the Taphrinomycotina have two tRNase ZLs, which we term tRNase ZL1 and tRNase ZL2, respectively. tRNase ZL1s and tRNase ZL2s have been either shown or predicted to localize to the nucleus and mitochondria, respectively (Additional file 2 and data not shown) . Since in the current databases, all sequenced Taphrinomycotina species come from only Schizosaccharomyces, it would be interesting to see whether species in other genera also contain two tRNase ZLs.
Like Ascomycota species, all sequenced Basidiomycota species (except for Agaricus bisporus) have a single tRNase ZL. However, unlike the situation in the Ascomycota, tRNase ZS was found in all sequenced Basidiomycota species. While the majority of Basidiomycota species have a single tRNase ZS, four Basidiomycota species (A. bisporus, C. cinerea, Laccaria bicolor and Postia placenta) have two tRNase ZSs. Among the Basidiomycota species examined, A. bisporus has the largest number of tRNase Zs (two tRNase ZLs and two tRNase ZSs). The number of tRNase Z seems to be variable in the three sequenced chytrid species. Allomyces macrogynus and Spizellomyces punctatus have two tRNase ZLs whereas Batrachochytrium dendrobatidis appears to have one tRNase ZL. Moreover, tRNase ZS was only identified in S. punctatus. Although we could not correctly predict the full-length tRNase Zs in three sequenced Zygomycota species (Rhizopus oryzae, Mucor circinellodes and Phycomyces blakesleeanus) and in three sequenced Microspordian species (Encephalitozoon cuniculi, Enterocytozoon bieneusi and Nosema ceranae), it is important to note that tRNase ZS appears to exist in all sequenced Zygomycota fungi but not in sequenced Microspordian fungi known for extreme genome reduction and compaction .
Analysis of fungal tRNase ZLs
The sizes of predicted tRNase ZLs vary considerably among fungal species, ranging from 648 to 1140 aa with an average size of ~924 aa. The variation in protein size is due to a high degree of length and sequence variation of N-terminal and C-terminal extensions and many insertions and/or deletions. Remarkably, tRNase ZLs in Sordaria macrospora and three Neurospora species (N. crassa, Neurospora discreta and Neurospora discreta) have a very long N-terminal extension (~200 residues). This feature appears to be family-specific since all these species belong to the family Sordariaceae.
A number of fungal tRNase ZLs have a variable length N-terminal extension predicted to contain a canonical MTS (Additional file 2). In addition, tRNase ZL2s from four Schizosaccharomyces species we have examined also contain a putative MTS in their N-terminal extensions. The N-terminal extensions found in fungal tRNase ZLs may serve as transit sequences for directing the proteins to the mitochondria. It is interesting to note that tRNase ZLs from D. melanogaster and humans also contain a canonical MTS.
Representatives of candidate fungal tRNase Zs used in sequence alignment
Ashbya gossypii (Ago)
Aspergillus nidulans (Ani)
Candida albicans (Cal)
Coccidioides immitis (Cim)
Fusarium graminearum (Fgr)
Neurospora crassa (Ncr)
Pyrenophora tritici-repentis (Ptr)
Saccharomyces cerevisiae (Sce)
Schizosaccharomyces pombe (Spo)
Schizosaccharomyces pombe (Spo)
Sclerotinia sclerotiorum (Ssc)
Yarrowia lipolytica (Yli)
Cryptococcus neoformans (Cne)
Malassezia globosa (Mgl)
Puccinia graminis (Pgr)
Spizellomyces punctatus (Spu)
Agaricus bisporus (Abi)
Agaricus bisporus (Abi)
Coprinopsis cinerea (Cci)
Coprinopsis cinerea (Cci)
Laccaria bicolor (Lbi)
Laccaria bicolor (Lbi)
Malassezia globosa (Mgl)
Melampsora laricis-populina (Mla)
Postia placenta (Ppl)
Postia placenta (Ppl)
Puccinia graminis (Pgr)
Spizellomyces punctatus (Spu)
Although tRNase ZLs from closely related species share the high degree of sequence similarity, the sequence similarity among fungal tRNase ZLs is low. Overall, sequence conservation of fungal tRNase ZLs is largely confined to highly conserved Motifs I-V (Motif II is also called the histidine motif) and the PxKxRN, HEAT and HST loop motifs at the C-terminal halves of the proteins (Figure 2). Except for the PxKxRN loop and Motif I, which were found to play a role in pre-tRNA acceptor stem binding and CCA anti-determination , all other motifs are involved in zinc binding and catalysis [31–33]. Motifs I-V contain invariant histidine and/or aspartate residues essential for the tRNase Z activity. In particular, the histidine and aspartate residues in the histidine motif, together with the histidine residues in Motifs III and V and the aspartate residue in Motif IV are involved in the coordination of the two zinc ions at the catalytic center [1, 3, 4].
All characterized tRNase Zs contain a characteristic domain of 30~50 aa residues, termed a flexible arm (or an exosite), which is important for substrate binding [34, 35]. Based on sequence comparison, three types of flexible arms, termed the zinc-dependent phosphodiesterase (ZiPD)-, ELAC2- and Thermotoga maritima (TM)-type flexible arms, have been found in tRNase Zs. The ZiPD- and ELAC2-type flexible arms were typical for bacterial tRNase ZSs and eukaryotic tRNase ZLs, respectively, whereas the TM-type flexible arm is an atypical one found only in tRNase ZSs from T. maritima and A. thaliana. Interestingly, tRNase ZS in T. maritima itself is an atypical enzyme as it cleaves CCA-containing pre-tRNAs after CCA. Although having sequence and length variations, both ZiPD- and ELAC2-type flexible arms comprise a GP motif rich in glycine and proline residues , followed by a Walker A-like motif . However, unlike the ZiPD- and ELAC2-type flexible arms, the TM-type flexible arm is short and lacks the GP motif. Instead it contains a short stretch of mainly basic amino acids .
Analysis of fungal tRNase ZSs
The phylogenetic distribution of tRNase Z in fungi
tRNase ZL is widespread in fungi and the majority of fungal species appear to have a single tRNase ZL. This latter finding is somewhat unexpected given striking differences in the genome size, life cycle and morphology of species of fungi. The Pezizomycotina and Saccharomycotina belong to later diverging fungi. Their genomes vary considerably in size due to gene gain and loss events including tandem gene duplication, whole-genome duplication and extensive gene loss [37, 38]. The Saccharomycotina genome sizes vary from ~9 (Pichia pastoris) to ~24 Mb (Candida parapsilosis), whereas the genome sizes of the Pezizomycotina fungi range from 23 (Microsporum canis) to 43 Mb (N. crassa). Furthermore, Pezizomycotina and Saccharomycotina fungi range in complexity from unicellular yeasts to filamentous molds. However, despite their remarkable differences in genome size, life cycle and morphology, the Pezizomycotina and Saccharomycotina fungi tend to contain only one tRNase ZL. These results indicate that the diversity of tRNase Z in fungi is not directly proportional to either the difference in genome size or the complexity of the life cycle and the morphology.
In contrast to tRNase ZL, tRNase ZS has a limited phylogenetic distribution. The apparent lack of the tRNase ZS gene in the genomes of Ascomycota species suggests that it has been deleted from the genomes of Ascomycota fungi. It is possible that tRNase ZS existed before the divergence of the Ascomycota from the Basidiomycota, and was subsequently lost after the appearance of a novel structure (tRNase ZL). This is supported by the finding that tRNase ZS is retained in the genomes of Basidiomycota, Chytridiomycota and Zygomycota fungi.
Most of Basidiomycota species examined contain one tRNase ZS. The fungal tRNase ZSs appear to be unique among all known tRNase ZSs in either lacking the flexible arm or having an atypical flexible arm (see discussion below). One explanation of existence of tRNase ZS genes in Basidiomycota species is that these genes may represent pseudo-tRNase ZS genes. Another explanation, which we favor, is that fungal tRNase ZS, at least some of them, may play a back-up or different role. Support for this hypothesis comes from recent studies of tRNase ZS in A. thaliana and humans. In A. thaliana, one tRNase ZS may represent a back-up for the nuclear tRNA 3'-end processing in case of dysfunction of nuclear-localized tRNase ZL, whereas the other plays a role in chloroplasts . In human cells, tRNase ZS is located in the cytosol and likely have substrates other than pre-tRNA .
The flexible arms of candidate fungi tRNase Zs display remarkable diversity
An unexpected and striking result of this analysis is the diversity of the flexible arms within candidate fungal tRNase Zs, particularly tRNase ZSs. The most characteristic features of the typical flexible arm found in tRNase Zs are the GP- and Walker-A like motifs. A subset of fungal tRNase ZLs and all fungal tRNase ZSs appear to lack the GP-motif, and some fungal tRNase ZSs do not seem to have the Walker A-like motif. In the most extreme case, the flexible arm is missing in fungal tRNase ZSs. It is not yet understood why the flexible arms of fungal tRNase Zs display diversity in the primary sequence.
It is interesting to note that four Basidiomycota species contain two candidate tRNase ZSs, one of which lacks the flexible arm. It is likely that these two tRNase ZSs form heterodimers that would look like tRNase ZL, where only the N-terminal half has a flexible arm.
The apparent lack of the GP-motif in some fungal tRNase ZLs and all fungal tRNase ZSs that we have examined raises the question of whether this motif is absolutely required for substrate binding. Structural and biochemical evidence suggests that the GP-motif may not be essential for pre-tRNA binding. To date, the three-dimensional structures of tRNase ZSs from B. subtilis, E. coli and T. maritima have been solved by X-ray crystallography [34, 40–43]. Remarkably, the flexible arms of T. maritima tRNase ZS lacking the GP motif and the other two tRNase ZSs harboring the GP motif have very similar structures, composed of a compact globular domain and an extended two-stranded stalk, which extrude from the tRNase ZS core. However, they have different lengths and globular domains. The globular domains at the end of the flexible arms of B. subtilis and E. coli tRNase ZSs are composed of two α-helices, two β-strands and one 310-helix, whereas the counterpart in T. maritima tRNase ZS consists of one very short α-helix, one long helix and one 310-helix. The conserved GP-motif, particularly the proline residues, appears to add rigidity to two flexible arm helices since it is localized between them . It would be interesting to know how the flexible arm lacking the GP-motif participates in substrate binding.
Recent biochemical studies have also suggested that the GP-motif may not be essential for substrate binding. Single alanine substitutions across the GP motif in D. melanogaster tRNase ZL only moderately affect substrate binding. In contrast, substitution of a conserved leucine residue at the boundary of the globular domain and stalk with alanine almost completely abolishes substrate binding as the globular domain deletion . Similarly, deletion of the GP motif in B. subtilis tRNase ZS does not eliminate pre-tRNA binding but alters the cleavage specificity of the enzyme . These results suggest that the GP motif may be important but not essential for substrate binding.
The evolutionary relationship between tRNase ZS and tRNase ZL
In eukaryotes, tRNase ZL appears to take over tRNase ZS in endonucleolytic 3'-end processing of pre-tRNAs, which raises the question of how it evolves. The protein sequence of tRNase ZS is much more similar to the C-terminal half of tRNase ZL than to the N-terminal half of tRNase ZL. Furthermore, the C-terminal half of tRNase ZL retains all conserved motifs for proper catalytic function but has lost the flexible arm involved in substrate binding, whereas the N-terminal half has lost all active motifs but contains the flexible arm. These observations led to the proposal that tRNase ZL has evolved from tRNase ZS by gene duplication and subsequent sequence divergence . To assess whether phylogenetic evidence exists that is consistent with this notion, we estimated phylogenetic relationships among fungal tRNase Zs by using a Bayesian phylogenetic method. The clustering of all fungal tRNase ZSs with representative bacterial tRNase ZSs support the notion that tRNase ZL comes to eukaryotes through duplication of tRNase ZS gene. Further evidence that the N-terminal half of tRNase ZL is derived from an ancient tRNase ZS comes from our findings that the N-terminal half of fungal tRNase ZL contains candidate pseudo-motifs and that these pseudo-motifs are present in the same relative order as their original motifs appeared in tRNase ZS. These pseudo-motifs likely represent relics of original tRNase ZS motifs that were inactivated during diversification of the eukaryotic tRNase ZL gene.
The reason for the adoption of tRNase ZL over tRNase ZS in eukaryotes is unknown. One possibility is that eukaryotic cells may require more efficient tRNase Z enzymes. Support for this proposal comes from biochemical characterization of human tRNase Zs. In vitro characterization of recombinant human tRNase Zs have shown that tRNase ZL cleaves pre-tRNA significantly more efficiently compared to tRNase ZS . Although strong structural evidence to support that tRNase ZL evolved into a more efficient enzyme than tRNase ZS is still lacking, it is interesting to note that tRNase ZS and tRNase ZL may have different processing center numbers which would make much difference in the efficiency of pre-tRNA 3'-end processing. Three-dimensional structures of three bacterial tRNase ZSs have revealed that the proteins form homodimers [34, 40–43]. In particular, the crystal structure of B. subtilis tRNase ZS in complex with tRNA shows that the dimer has two identical processing centers with two substrate binding and catalytic sites. In contrast, a molecular modeling study has suggested that both the N-terminal and C-terminal halves of human tRNase ZL can fold into two distinct MBL domains with one domain containing a fully functionally catalytic site and the other containing a candidate substrate binding domain .
Schizosaccharomyces fission yeasts have two tRNase ZLs, most likely targeted to the nucleus and mitochondria, respectively
Schizosaccharomyces fission yeasts including S. pombe appear to be unique among the Ascomycota in having two tRNase ZLs (tRNase ZL1 and tRNase ZL2) that appear to be localized to the nucleus and mitochondria, respectively, as suggested in our previous study of tRNase ZL in S. pombe . Our fungal tRNase Z phylogeny shows that the two tRNase ZLs in fission yeasts may have arisen through gene duplication (Figure 1). Although the whole genome duplication found in the Saccharomycotina yeasts does not seem to occur in fission yeasts such as S. pombe , the tRNase ZL gene could be duplicated by other mechanisms such as tandem and segmental gene duplication. The two tRNase ZLs in each Schizosaccharomyces species all have very limited homology with each other (20-23% identity and 31-33% similarity, see Additional file 5), indicating that these two proteins have diverged considerably from each other since their duplication. It is interesting to note that the tRNase ZL gene could also be duplicated in non-fungal eukaryotic species A. thaliana. However, the two plant tRNase ZLs are highly related to each other (69% identity and 72% similarity).
Why do Schizosaccharomyces fission yeasts have two tRNase ZLs? In our previous study, it was found that the nuclear-targeted tRNase ZL1 (SpoTrz1) is involved in nuclear pre-tRNA 3'-end processing in S. pombe . Furthermore, its function can be compensated by either S. cerevisiae or human tRNase ZL. Although the role of mitochondrial-targeted tRNase ZL2 (SpoTrz2) remains to be determined, it is likely that this protein plays an essential role in mitochondrial RNA processing . Based on these results, it is possibly that the presence of two tRNase ZLs reflects that the nuclear and mitochondrial tRNA processing activities are associated with two different tRNase ZLs in Schizosaccharomyces fission yeasts. This may also hold true for wheat and potato since in these plants, enzymes involved in nuclear and mitochondrial tRNA 3'-end processing appear to be different . However, it is important to note that the nuclear and organelle tRNase Z activities in the majority of organisms described to date seem to reside in the same enzyme.
A survey of fungal databases shows that tRNase ZL appears to be universally present in fungi, whereas the presence of tRNase ZS is restricted to certain fungal phyla, indicative of the fundamental role of tRNase ZL in eukaryotic tRNA biogenesis. The apparent lack of tRNase ZS in the Ascomycota suggests that tRNase ZS may have lost before divergence of the Ascomycota and the Basidiomycota. A striking aspect of the tRNase ZL distribution is that there are two different tRNase ZLs in Schizosaccharomyces fission yeasts. These two tRNase ZLs are likely present in different cellular compartments, suggesting functional partitioning between these two proteins. Phylogenetic analysis suggests that tRNase ZS is ancestral to tRNase ZL and that tRNase ZL gene duplications may have occurred in certain fungal taxa, including Schizosaccharomyces fission yeasts. Sequence analysis reveals that the domain architecture of tRNase ZLs is highly conserved among fungi and metazoa. A surprising result of sequence analysis is the sequence diversity in the putative flexible arm of candidate fungal tRNase Zs. Our analysis also reveals pseudo-motifs at the N-terminal halves of tRNase ZLs. These findings support the view that tRNase ZL evolved through duplication and divergence of the tRNase ZS gene.
Fungal genome database search and protein sequence analysis
To identify candidate tRNase Zs, we conducted BLAST and PSI-BLAST searches using the known tRNase Z protein sequences as queries against fungal genomes databases including the National Center for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov/sutils/genom_table.cgi?organism=fungi), the Broad Institute http://www.broadinstitute.org/science/data, the Joint Genome Institute http://genome.jgi-psf.org/pages/fungi/home.jsf, the Genome Center at Washington University http://genome.wustl.edu/ and the Universal Protein Resource http://www.uniprot.org. All candidate sequences were obtained by using the cut-off E-value of 0.01. All candidate proteins were subjected to validation, which was carried out by using a variety tests that evaluate the likelihood of annotation errors and the amino acid sequence conservation within and among taxonomic groups. First, confirmation of true candidate tRNase Zs was done by back-searching individual candidate protein sequence against the GenBank database. Second, the gene sequences for the predicted tRNase Zs were manually checked for possible sequence gaps. Third, multiple protein sequence alignment was used to identify candidate proteins that were discordant due to possible genomic sequencing errors and/or intron misprediction. In most cases, we changed the splicing pattern of candidate tRNase Z either using gene prediction programs Fgenesh http://linux1.softberry.com/berry.phtml?topic=fgenesh%26group=programs%26subgroup=gfind and Geneid http://genome.crg.es/geneid.html or manually to restore the high degree of sequence conservation. Multiple sequence alignment was performed by using Clustal W , and the resulting alignment was further manually examined and adjusted to improve the detection of conserved regions. The putative subcellular localization signals of tRNase Zs were predicted by using the programs MitoProt http://ihg2.helmholtz-muenchen.de/ihg/mitoprot.html and PSORT II http://psort.hgc.jp/.
Phylogenetic analysis of fungal tRNase Zs
Full-length amino acid sequences of tRNase Zs from fungi and two bacteria, B. subtilis and E. coli were aligned by using Clustal W implemented in Mega 4.0 . Conserved regions were selected and ambiguous aligned regions were removed by using the program Gblocks 0.91b . tRNase ZSs from B. subtilis and E. coli were chosen as reference. The phylogenies were estimated by Bayesian inference with MrBayes 3.1.2  using a mixture of the fixed amino acid models and the gamma distribution. Statistical confidence was assessed by using Markov Chain Monte Carlo (MCMC) sampling approaches. Four simultaneous Markov chains were run for one million generations sampling every 1,000 generation in two replicate runs. The first 250 trees were discarded as burn-in and the convergence of the chains was evaluated using AWTY implemented in MrBayes 3.1.2 .
List of abbreviations
- tRNase Z:
tRNA 3' endonuclease
- tRNase ZS:
the short form of tRNase Z
- tRNase ZL:
the long form of tRNase Z
nuclear localization signal
mitochondrial targeting signal
We thank Dr. Jie Yan and two anonymous reviewers for improving the quality of the manuscript. We are grateful to Dr. Guang Yang for the use of his laboratory facilities to construct the phylogenetic tree. This work was supported by grants from the National Science Foundation of China (30771178) and Nanjing Normal University (2007104XGQ0148).
- Ceballos M, Vioque A: tRNase Z. Protein Pept Lett. 2007, 14: 137-145. 10.2174/092986607779816050.View ArticlePubMedGoogle Scholar
- 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
- 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
- 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
- Aravind L: An evolutionary classification of the metallo-beta-lactamase fold proteins. In Silico Biol. 1999, 1: 69-91.PubMedGoogle Scholar
- 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
- 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
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- Sullivan KD, Steiniger M, Marzluff WF: A core complex of CPSF73, CPSF100, and Symplekin may form two different cleavage factors for processing of poly(A) and histone mRNAs. Mol Cell. 2009, 34: 322-332. 10.1016/j.molcel.2009.04.024.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- 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
- 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. 10.1104/pp.109.137950.View ArticleGoogle Scholar
- 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
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (1): 135-142. 10.1016/j.molcel.2009.12.016.PubMed CentralView ArticlePubMedGoogle Scholar
- Elbarbary RA, Takaku H, Uchiumi N, Tamiya H, Abe M, Nishida H, Nashimoto M: Human cytosolic tRNase ZL can downregulate gene expression through miRNA. FEBS Lett. 2009, 583: 3241-3246. 10.1016/j.febslet.2009.09.015.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawksworth DL: Fungal diversity and its implication for genetic resource collections. Stud Mycol. 2004, 50: 9-18.Google Scholar
- Fitzpatrick DA, Logue ME, Stajich JE, Butler G: A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis. BMC Evol Biol. 2006, 6: 99-10.1186/1471-2148-6-99.PubMed CentralView ArticlePubMedGoogle Scholar
- James TY, Kauff F, Schoch CL, Matheny PB, Hofstetter V, Cox CJ, Celio G, Gueidan C, Fraker E, Miadlikowska J, Lumbsch HT, Rauhut A, Reeb V, Arnold AE, Amtoft A, Stajich JE, Hosaka K, Sung GH, Johnson D, O'Rourke B, Crockett M, Binder M, Curtis JM, Slot JC, Wang Z, Wilson AW, Schüssler A, Longcore JE, O'Donnell K, Mozley-Standridge S, et al: Reconstructing the early evolution of Fungi using a six-gene phylogeny. Nature. 2006, 443: 818-822. 10.1038/nature05110.View ArticlePubMedGoogle Scholar
- Liu Y, Leigh JW, Brinkmann H, Cushion MT, Rodriguez-Ezpeleta N, Philippe H, Lang BF: Phylogenomic analyses support the monophyly of Taphrinomycotina, including Schizosaccharomyces fission yeasts. Mol Biol Evol. 2009, 26: 27-34. 10.1093/molbev/msn221.View ArticlePubMedGoogle Scholar
- Robbertse B, Reeves JB, Schoch CL, Spatafora JW: A phylogenomic analysis of the Ascomycota. Fungal Genet Biol. 2006, 43: 715-725. 10.1016/j.fgb.2006.05.001.View ArticlePubMedGoogle Scholar
- Hibbett DS, Binder M, Bischoff JF, Blackwell M, Cannon PF, Eriksson OE, Huhndorf S, James T, Kirk PM, Lücking R, Thorsten Lumbsch H, Lutzoni F, Matheny PB, McLaughlin DJ, Powell MJ, Redhead S, Schoch CL, Spatafora JW, Stalpers JA, Vilgalys R, Aime MC, Aptroot A, Bauer R, Begerow D, Benny GL, Castlebury LA, Crous PW, Dai YC, Gams W, Geiser DM, et al: A higher-level phylogenetic classification of the Fungi. Mycol Res. 2007, 111: 509-547. 10.1016/j.mycres.2007.03.004.View ArticlePubMedGoogle Scholar
- Takaku H, Minagawa A, Takagi M, Nashimoto M: A candidate prostate cancer susceptibility gene encodes tRNA 3' processing endoribonuclease. Nucleic Acids Res. 2003, 31: 2272-2278. 10.1093/nar/gkg337.PubMed CentralView ArticlePubMedGoogle Scholar
- Takaku H, Minagawa A, Takagi M, Nashimoto M: The N-terminal half-domain of the long form of tRNase Z is required for the RNase 65 activity. Nucleic Acids Res. 2004, 32: 4429-4438. 10.1093/nar/gkh774.PubMed CentralView ArticlePubMedGoogle Scholar
- Matsuyama A, Arai R, Yashiroda Y, Shirai A, Kamata A, Sekido S, Kobayashi Y, Hashimoto A, Hamamoto M, Hiraoka Y, Horinouchi S, Yoshida M: ORFeome cloning and global analysis of protein localization in the fission yeast Schizosaccharomyces pombe. Nat Biotechnol. 2006, 24: 841-847. 10.1038/nbt1222.View ArticlePubMedGoogle Scholar
- Keeling PJ, Slamovits CH: Simplicity and complexity of microsporidian genomes. Eukaryot Cell. 2004, 3: 1363-1369. 10.1128/EC.3.6.1363-1369.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- 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
- 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
- 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. 10.1074/jbc.M500591200.View ArticlePubMedGoogle Scholar
- 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
- Dujon B, Sherman D, Fischer G, Durrens P, Casaregola S, Lafontaine I, De Montigny J, Marck C, Neuvéglise C, Talla E, Goffard N, Frangeul L, Aigle M, Anthouard V, Babour A, Barbe V, Barnay S, Blanchin S, Beckerich JM, Beyne E, Bleykasten C, Boisramé A, Boyer J, Cattolico L, Confanioleri F, De Daruvar A, Despons L, Fabre E, Fairhead C, Ferry-Dumazet H, et al: Genome evolution in yeasts. Nature. 2004, 430: 35-44. 10.1038/nature02579.View ArticlePubMedGoogle Scholar
- Kellis M, Birren BW, Lander ES: Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae. Nature. 2004, 428: 617-624. 10.1038/nature02424.View ArticlePubMedGoogle Scholar
- 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
- 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
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- Yan H, Zareen N, Levinger L: Naturally occurring mutations in human mitochondrial pre-tRNASer(UCN) can affect the transfer ribonuclease Z cleavage site, processing kinetics, and substrate secondary structure. J Biol Chem. 2006, 281: 3926-3935. 10.1074/jbc.M509822200.View ArticlePubMedGoogle Scholar
- Wood V, Gwilliam R, Rajandream MA, Lyne M, Lyne R, Stewart A, Sgouros J, Peat N, Hayles J, Baker S, Basham D, Bowman S, Brooks K, Brown D, Brown S, Chillingworth T, Churcher C, Collins M, Connor R, Cronin A, Davis P, Feltwell T, Fraser A, Gentles S, Goble A, Hamlin N, Harris D, Hidalgo J, Hodgson G, Holroyd S, et al: The genome sequence of Schizosaccharomyces pombe. Nature. 2002, 415: 871-880. 10.1038/nature724.View ArticlePubMedGoogle Scholar
- Mayer M, Schiffer S, Marchfelder A: tRNA 3' processing in plants: nuclear and mitochondrial activities differ. Biochemistry. 2000, 39 (8): 2096-2105. 10.1021/bi992253e.View ArticlePubMedGoogle Scholar
- 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
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-552.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Nylander JA, Wilgenbusch JC, Warren DL, Swofford DL: AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics. Bioinformatics. 2008, 24: 581-583. 10.1093/bioinformatics/btm388.View ArticlePubMedGoogle Scholar
- 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 (1): 255-262. 10.1093/nar/gkh182.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (17): 4534-4543. 10.1093/emboj/cdg435.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (1): 151-160. 10.1016/j.ydbio.2003.10.016.View ArticlePubMedGoogle Scholar
- Ezraty B, Dahlgren B, Deutscher MP: The RNase Z homologue encoded by Escherichia coli elaC gene is RNase BN. J Biol Chem. 2005, 280 (17): 16542-16545. 10.1074/jbc.C500098200.View ArticlePubMedGoogle Scholar
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.