Structure, function, and evolution of the tRNA endonucleases of Archaea: an example of subfunctionalization

Complexes of tRNA and maturation enzymes: shaping up for translation

Several significant structures of transfer ribonucleic acid (tRNA) maturation enzymes complexed with precursor tRNA or fragments thereof have been published recently, providing detailed knowledge of enzyme-tRNA recognition and catalytic strategies. In addition to reinforcing the general principles of RNA-protein interaction, the new structures highlight both the features of composite RNA recognition by multiple enzyme subunits and the pronounced RNA structural flexibility in or near the active site in all cases. These structural principles provide plausible explanations for the exquisite specificity and catalytic power of these enzymes and, in the case of evolutionary adaptation, for the ability of some enzymes to develop novel specificities.

In silico screening of archaeal tRNA-encoding genes having multiple introns with bulge-helix-bulge splicing motifs

In archaeal species, several transfer RNA genes have been reported to contain endogenous introns. Although most of the introns are located at anticodon loop regions between nucleotide positions 37 and 38, a number of introns at noncanonical sites and six cases of tRNA genes containing two introns have also been documented. However, these tRNA genes are often missed by tRNAscan-SE, the software most widely used for the annotation of tRNA genes. We previously developed SPLITS, a computational tool to identify tRNA genes containing one intron at a noncanonical position on the basis of its discriminative splicing motif, but the software was limited in the detection of tRNA genes with multiple introns at noncanonical sites. In this study, we initially updated the system as SPLITSX in order to correctly predict known tRNA genes as well as novel ones with multiple introns. By a comprehensive search for tRNA genes in 29 archaeal genomes using SPLITSX, we listed 43 novel candidates that contain introns at noncanonical sites. As a result, 15 contained two introns and three contained three introns within the respective putative tRNA genes. Moreover, the candidates completely complemented all the codons of two archaeal species of uncultured methanogenic archaeon, RC-I and Thermofilum pendens Hrk 5, with novel candidates that were not detectable by tRNAscan-SE alone.

Genomic analysis of the uncultivated marine crenarchaeote Cenarchaeum symbiosum

Crenarchaeota are ubiquitous and abundant microbial constituents of soils, sediments, lakes, and ocean waters. To further describe the cosmopolitan nonthermophilic Crenarchaeota, we analyzed the genome sequence of one representative, the uncultivated sponge symbiont Cenarchaeum symbiosum. C. symbiosum genotypes coinhabiting the same host partitioned into two dominant populations, corresponding to previously described a- and b-type ribosomal RNA variants. Although they were syntenic, overlapping a- and b-type ribotype genomes harbored significant variability. A single tiling path comprising the dominant a-type genotype was assembled and used to explore the genomic properties of C. symbiosum and its planktonic relatives. Of 2,066 ORFs, 55.6% matched genes with predicted function from previously sequenced genomes. The remaining genes partitioned between functional RNAs (2.4%) and hypotheticals (42%) with limited homology to known functional genes. The latter category included some genes likely involved in the archaeal-sponge symbiotic association. Conversely, 525 C. symbiosum ORFs were most highly similar to sequences from marine environmental genomic surveys, and they apparently represent orthologous genes from free-living planktonic Crenarchaeota. In total, the C. symbiosum genome was remarkably distinct from those of other known Archaea and shared many core metabolic features in common with its free-living planktonic relatives.

Three-dimensional structure determined for a subunit of human tRNA splicing endonuclease (Sen15) reveals a novel dimeric fold

Splicing of eukaryal intron-containing tRNAs requires the action of the heterotetrameric splicing endonuclease, which is composed of two catalytic subunits, Sen34 and Sen2, and two structural subunits, Sen15 and Sen54. Here we report the solution structure of the human tRNA splicing endonuclease subunit HsSen15. To facilitate the structure determination, we removed the disordered 35 N-terminal and 14 C-terminal residues of the full-length protein to produce HsSen15(36-157). The structure of HsSen15(36-157), the first for a subunit of a eukaryal splicing endonuclease, revealed that the protein possesses a novel homodimeric fold. Each monomer consists of three alpha-helices and a mixed antiparallel/parallel beta-sheet, arranged in a topology similar to that of the C-terminal domain of Methanocaldococcus jannaschii endonuclease. The dimeric interface is dominated by a beta-barrel structure, formed by face-to-face packing of two, three-stranded beta-sheets. Each of the beta-sheets results from reciprocal parallel pairing of one beta-strand from one subunit with two other beta-strands from the symmetric subunit. The structural model provides insights into the functional assembly of the human tRNA splicing endonuclease.

Divergence of the Dof gene families in poplar, Arabidopsis, and rice suggests multiple modes of gene evolution after duplication

It is widely accepted that gene duplication is a primary source of genetic novelty. However, the evolutionary fate of duplicated genes remains largely unresolved. The classical Ohno's Duplication-Retention-Non/Neofunctionalization theory, and the recently proposed alternatives such as subfunctionalization or duplication-degeneration-complementation, and subneofunctionalization, each can explain one or more aspects of gene fate after duplication. Duplicated genes are also affected by epigenetic changes. We constructed a phylogenetic tree using Dof (DNA binding with one finger) protein sequences from poplar (Populus trichocarpa) Torr. & Gray ex Brayshaw, Arabidopsis (Arabidopsis thaliana), and rice (Oryza sativa). From the phylogenetic tree, we identified 27 pairs of paralogous Dof genes in the terminal nodes. Analysis of protein motif structure of the Dof paralogs and their ancestors revealed six different gene fates after gene duplication. Differential protein methylation was revealed between a pair of duplicated poplar Dof genes, which have identical motif structure and similar expression pattern, indicating that epigenetics is involved in evolution. Analysis of reverse transcription-PCR, massively parallel signature sequencing, and microarray data revealed that the paralogs differ in expression pattern. Furthermore, analysis of nonsynonymous and synonymous substitution rates indicated that divergence of the duplicated genes was driven by positive selection. About one-half of the motifs in Dof proteins were shared by non-Dof proteins in the three plants species, indicating that motif co-option may be one of the forces driving gene diversification. We provided evidence that the Ohno's Duplication-Retention-Non/Neofunctionalization, subfunctionalization/duplication-degeneration-complementation, and subneofunctionalization hypotheses are complementary with, not alternative to, each other.

Archaeal pre-mRNA splicing: a connection to hetero-oligomeric splicing endonuclease

Eukaryotic Cbf5 is a protein subunit of the small nucleolar RNA-protein complex. Previously, we identified, in archaeal homologs of cbf5 of the crenarchaea, Aeropyrum pernix, Sulfolobus solfataricus, and Sulfolobus tokodaii, the first examples of introns of archaeal protein-coding genes. Here, we report the immunological detection of Cbf5 protein of S. tokodaii, the product of the spliced cbf5 mRNA. The hetero-oligomeric splicing endonuclease activity from recombinant S. tokodaii subunits cleaved at the exon-intron boundaries of cbf5 pre-mRNA fragments,suggesting that synthesis of full-length Cbf5 protein requires this activity. Database searches and PCR screens identified additional cbf5 introns in some, but not all sequenced crenarchaeal genomes. The predicted secondary structures of exon-intron boundaries of many of the newly identified intron-containing cbf5 pre-mRNAs contained relaxed forms of the bulge-helix-bulge motif similar to that of S. tokodaii. These observations are consistent with previous reports indicating that subunit composition of the splicing endonuclease contributes to substrate specificity.

RNA recognition and cleavage by a splicing endonuclease

The RNA splicing endonuclease cleaves two phosphodiester bonds within folded precursor RNAs during intron removal, producing the functional RNAs required for protein synthesis. Here we describe at a resolution of 2.85 angstroms the structure of a splicing endonuclease from Archaeglobus fulgidus bound with a bulge-helix-bulge RNA containing a noncleaved and a cleaved splice site. The endonuclease dimer cooperatively recognized a flipped-out bulge base and stabilizes sharply bent bulge backbones that are poised for an in-line RNA cleavage reaction. Cooperativity arises because an arginine pair from one catalytic domain sandwiches a nucleobase within the bulge cleaved by the other catalytic domain.

The RNA polymerase III-dependent family of genes in hemiascomycetes: comparative RNomics, decoding strategies, transcription and evolutionary implications

We present the first comprehensive analysis of RNA polymerase III (Pol III) transcribed genes in ten yeast genomes. This set includes all tRNA genes (tDNA) and genes coding for SNR6 (U6), SNR52, SCR1 and RPR1 RNA in the nine hemiascomycetes Saccharomyces cerevisiae, Saccharomyces castellii, Candida glabrata, Kluyveromyces waltii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and the archiascomycete Schizosaccharomyces pombe. We systematically analysed sequence specificities of tRNA genes, polymorphism, variability of introns, gene redundancy and gene clustering. Analysis of decoding strategies showed that yeasts close to S.cerevisiae use bacterial decoding rules to read the Leu CUN and Arg CGN codons, in contrast to all other known Eukaryotes. In D.hansenii and C.albicans, we identified a novel tDNA-Leu (AAG), reading the Leu CUU/CUC/CUA codons with an unusual G at position 32. A systematic 'p-distance tree' using the 60 variable positions of the tRNA molecule revealed that most tDNAs cluster into amino acid-specific sub-trees, suggesting that, within hemiascomycetes, orthologous tDNAs are more closely related than paralogs. We finally determined the bipartite A- and B-box sequences recognized by TFIIIC. These minimal sequences are nearly conserved throughout hemiascomycetes and were satisfactorily retrieved at appropriate locations in other Pol III genes.

Leafing through the genomes of our major crop plants: strategies for capturing unique information

Crop plants not only have economic significance, but also comprise important botanical models for evolution and development. This is reflected by the recent increase in the percentage of publicly available sequence data that are derived from angiosperms. Further genome sequencing of the major crop plants will offer new learning opportunities, but their large, repetitive, and often polyploid genomes present challenges. Reduced-representation approaches - such as EST sequencing, methyl filtration and Cot-based cloning and sequencing - provide increased efficiency in extracting key information from crop genomes without full-genome sequencing. Combining these methods with phylogenetically stratified sampling to allow comparative genomic approaches has the potential to further accelerate progress in angiosperm genomics.

Buffering of crucial functions by paleologous duplicated genes may contribute cyclicality to angiosperm genome duplication

Genome duplication followed by massive gene loss has permanently shaped the genomes of many higher eukaryotes, particularly angiosperms. It has long been believed that a primary advantage of genome duplication is the opportunity for the evolution of genes with new functions by modification of duplicated genes. If so, then patterns of genetic diversity among strains within taxa might reveal footprints of selection that are consistent with this advantage. Contrary to classical predictions that duplicated genes may be relatively free to acquire unique functionality, we find among both Arabidopsis ecotypes and Oryza subspecies that SNPs encode less radical amino acid changes in genes for which there exists a duplicated copy at a "paleologous" locus than in "singleton" genes. Preferential retention of duplicated genes encoding long complex proteins and their unexpectedly slow divergence (perhaps because of homogenization) suggest that a primary advantage of retaining duplicated paleologs may be the buffering of crucial functions. Functional buffering and functional divergence may represent extremes in the spectrum of duplicated gene fates. Functional buffering may be especially important during "genomic turmoil" immediately after genome duplication but continues to act approximately 60 million years later, and its gradual deterioration may contribute cyclicality to genome duplication in some lineages.

The heteromeric Nanoarchaeum equitans splicing endonuclease cleaves noncanonical bulge-helix-bulge motifs of joined tRNA halves

Among the tRNA population of the archaeal parasite Nanoarchaeum equitans are five species assembled from separate 5' and 3' tRNA halves and four species derived from tRNA precursors containing introns. In both groups an intervening sequence element must be removed during tRNA maturation. A bulge-helix-bulge (BHB) motif is the hallmark structure required by the archaeal splicing endonuclease for recognition and excision of all introns. BHB motifs are recognizable at the joining sites of all five noncontinuous tRNA species, although deviations from the canonical BHB motif are clearly present in at least two of them. Here, we show that the N. equitans splicing endonuclease cleaves tRNA precursors containing normal introns, as well as all five noncontinuous precursor tRNAs, at the predicted splice sites, indicating the enzyme's dual role in the removal of tRNA introns and processing of tRNA halves to be joined in trans. The cleavage activity on a set of synthetic canonical and noncanonical BHB constructs showed that the N. equitans splicing endonuclease accepts a broader range of substrates than the homodimeric Archaeoglobus fulgidus enzyme. In contrast to the A. fulgidus endonuclease, the N. equitans splicing enzyme possesses two different subunits. This heteromeric endonuclease type, found in N. equitans, in all Crenarchaeota, and in Methanopyrus kandleri, is able to act on the noncanonical tRNA introns present only in these organisms, which suggests coevolution of enzyme and substrate.

Coevolution of tRNA intron motifs and tRNA endonuclease architecture in Archaea

Members of the three kingdoms of life contain tRNA genes with introns. The introns in pre-tRNAs of Bacteria are self-splicing, whereas introns in archaeal and eukaryal pre-tRNAs are removed by splicing endonucleases. We have studied the structures of the endonucleases of Archaea and the architecture of the sites recognized in their pre-tRNA substrates. Three endonuclease structures are known in the Archaea: a homotetramer in some Euryarchaea, a homodimer in other Euryarchaea, and a heterotetramer in the Crenarchaeota. The homotetramer cleaves only the canonical bulge-helix-bulge structure in its substrates. Variants of the substrate structure, termed bulge-helix-loops, appear in the pre-tRNAs of the Crenarcheota and Nanoarcheota. These variant structures can be cleaved only by the homodimer or heterotetramer forms of the endonucleases. Thus, the structures of the endonucleases and their substrates appear to have evolved together.

Structural characterization of the catalytic subunit of a novel RNA splicing endonuclease

The RNA splicing endonuclease is responsible for recognition and excision of nuclear tRNA and all archaeal introns. Despite the conserved RNA cleavage chemistry and a similar enzyme assembly, currently known splicing endonuclease families have limited RNA specificity. Different from previously characterized splicing endonucleases in Archaea, the splicing endonuclease from archaeum Sulfolobus solfataricus was found to contain two different subunits and accept a broader range of substrates. Here, we report a crystal structure of the catalytic subunit of the S.solfataricus endonuclease at 3.1 angstroms resolution. The structure, together with analytical ultracentrifugation analysis, identifies the catalytic subunit as an inactive but stable homodimer, thus suggesting the possibility of two modes of functional assembly for the active enzyme.

Functional reconstitution of a crenarchaeal splicing endonuclease in vitro

Sulfolobus tokodaii strain 7 is one of Crenarchaea whose entire genome has been sequenced. The genome sequence revealed that it possesses two open reading frames (ORFs) that are homologous to EndA, a protein responsible for splicing endonuclease activity in Archaea. Interestingly, one of the two ORFs lacks a putative catalytic amino acid residue for the nuclease activity. To investigate their functions, the two ORF products were individually expressed in Escherichia coli, partially purified, and tested for their nuclease activities in vitro. Using in vitro transcribed tRNA precursor as a substrate, we found that the two ORF products are concurrently required to cleave exon-intron junctions. Our finding implies that the splicing endonuclease for the organism is a multi-subunit complex composed of the two endA gene products.