Click
here to close Hello! We notice that
you are using Internet Explorer, which is not supported by Echinobase
and may cause the site to display incorrectly. We suggest using a
current version of Chrome,
FireFox,
or Safari.
Mol Biol Evol
2021 Jan 04;381:215-228. doi: 10.1093/molbev/msaa203.
Show Gene links
Show Anatomy links
Monophyletic Origin and Divergent Evolution of Animal Telomerase RNA.
Logeswaran D
,
Li Y
,
Podlevsky JD
,
Chen JJ
.
Abstract
Telomerase RNA (TR) is a noncoding RNA essential for the function of telomerase ribonucleoprotein. TRs from vertebrates, fungi, ciliates, and plants exhibit extreme diversity in size, sequence, secondary structure, and biogenesis pathway. However, the evolutionary pathways leading to such unusual diversity among eukaryotic kingdoms remain elusive. Within the metazoan kingdom, the study of TR has been limited to vertebrates and echinoderms. To understand the origin and evolution of TR across the animal kingdom, we employed a phylogeny-guided, structure-based bioinformatics approach to identify 82 novel TRs from eight previously unexplored metazoan phyla, including the basal-branching sponges. Synthetic TRs from two representative species, a hemichordate and a mollusk, reconstitute active telomerase in vitro with their corresponding telomerase reverse transcriptase components, confirming that they are authentic TRs. Comparative analysis shows that three functional domains, template-pseudoknot (T-PK), CR4/5, and box H/ACA, are conserved between vertebrate and the basal metazoan lineages, indicating a monophyletic origin of the animal TRs with a snoRNA-related biogenesis mechanism. Nonetheless, TRs along separate animal lineages evolved with divergent structural elements in the T-PK and CR4/5 domains. For example, TRs from echinoderms and protostomes lack the canonical CR4/5 and have independently evolved functionally equivalent domains with different secondary structures. In the T-PK domain, a P1.1 stem common in most metazoan clades defines the template boundary, which is replaced by a P1-defined boundary in vertebrates. This study provides unprecedented insight into the divergent evolution of detailed TR secondary structures across broad metazoan lineages, revealing ancestral and later-diversified elements.
Fig. 1. A phylogeny-assisted approach for TR identification. (A) The workflow of phylogeny-assisted reiterative homology search for novel TRs. An initial PWM and secondary structure pattern derived from known TRs was used to search for novel TR candidates in the genomes or transcriptomes of closely related species. The newly identified TR sequences were used to improve the PWM and structural pattern for searching TR candidates in distantly related species. (B) Comparison of vertebrate and echinoderm TR secondary structures. Two structural domains, T-PK and box H/ACA domains, are conserved in both vertebrate and echinoderm TRs, whereas the vertebrate CR4/5 domain that contains P6.1 stem-loop is replaced with a functionally equivalent eCR4/5 domain in echinoderms. The size ranges of vertebrate and echinoderm TRs are shown. (C) PWMs of TR sequences of the pseudoknot and H/ACA domains. PWMs are derived from the multiple-sequence alignment of 42 vertebrate TRs (Chen et al. 2000) and 13 echinoderm TRs (Li et al. 2013; Podlevsky et al. 2016b). Conserved sequence motifs and base-paired regions within each domain are indicated underneath the matrices
Fig. 2. Secondary structures of deuterostome TRs. (A) Evolutionary relationship of deuterostome phyla and classes with TR identified. The numbers of TRs identified in this and previous studies for each class are indicated. An asterisk denotes the classes for which the TR secondary structure of a representative species is shown in BâD. Fungi represents the outgroup in the phylogenetic tree. (BâD) Representative TR secondary structures determined by phylogenetic comparative sequence analyses are shown for Asymmetron lucayanum (lancelet) from phylum chordata (B), Saccoglossus kowalevskii (acorn worm) from phylum hemichordata (C), and Acanthaster planci (sea star) from phylum Echinodermata (D). The conserved TR structural domains are shaded in blue. Universal covariations (thick lines), invariant residues, and residues with >80% conservation are based on the sequence alignment of 55 previously identified animal TRs and 82 novel metazoan TRs identified in this study. Group-specific covariations (thick lines) are indicated and based on the sequence alignment of TRs from individual groups including 50 chordate TRs (42 previously identified and eight novel), 23 echinoderm (13 previously identified and ten novel), and four acorn worm TRs.
Fig. 3. Functional validation and characterization of acorn worm TR. (A) Direct telomerase activity assay of acorn worm telomerase reconstituted in vitro. (top) Template sequence of acorn worm TR (open box) with base-pairing of six permuted telomeric DNA primers aâf. Sequence and number of expected nucleotides added are depicted for each primer. (bottom) Direct primer-extension assay of acorn worm telomerase. Acorn worm telomerase was in vitro reconstituted from T7 transcribed SkoTR (436ânt) and SkoTERT synthesized in RRL. The reconstituted acorn worm telomerase was analyzed with six permuted telomeric DNA primers (lanes 2â7). A reaction omitting SkoTR was included as a negative control (lane 1). A 32P end-labeled 18-mer oligonucleotide was added to each reaction as recovery control (r.c.) prior to ethanol precipitation of DNA products. Numbers to the right of the gel indicate the number of repeats or nucleotides added to the primer. (B) Two essential fragments of SkoTR. The T-PK and CR4/5 fragments of SkoTR were synthesized separately and assembled with SkoTERT in RRL, followed by telomerase activity assay. The schematic secondary structures of the SkoTR T-PK (top) and CR4/5 (bottom) fragments. Nucleotide numbers denote the 5â²- and 3â²-ends of the T7 transcribed RNA fragments, T-PK and CR4/5. Positions of two highly conserved U residues (U289 and U290) in the P6.1 loop of CR4/5 domain are shown. (C) Minimal requirement of TR domains for telomerase activity. T7 transcribed SkoTR fragments, T-PK (nt 1â184) and CR4/5 (nt 237â307), were assembled with in vitro synthesized SkoTERT and analyzed for activity. The CR4/5 fragments with a P6.1 substitution (U289C or U290C) were assembled with T-PK fragment and SkoTERT and assayed for activity with the primer (GGGTTA)3. A 32P end-labeled 18-mer oligonucleotide was added to each reaction prior to ethanol precipitation of DNA products as r.c. The numbers of repeats added to the primer are shown to the right of the gel.
Fig. 4. Secondary structures of protostome TRs. (A) Evolutionary relationship of protostome phyla and classes with TR identified. The numbers of TRs identified in this study for each class are indicated. An asterisk denotes the classes for which the TR secondary structure of a representative species is shown in BâD. Fungi represents the outgroup in the phylogenetic tree. (BâD) Representative TR secondary structures determined by phylogenetic comparative sequence analyses are shown; Pomacea diffusa (apple snail) from phylum mollusca (B), Eisenia fetida (earth worm) from phylum annelida (C), and Lingula anatina (tongue shell) from phylum brachipoda (D). The conserved TR structural domains are shaded in blue. The eCR4/5 domain is shaded in red. Universal covariations (thick lines), invariant residues, and residues with >80% conservation are based on the sequence alignment of 55 previously identified animal TRs and 82 novel metazoan TRs identified in this study. Group-specific covariations (thick lines) are indicated and based on the sequence alignment of TRs from individual groups including 21 mollusca, 7Â annelida, 1 brachiopda, and 1 phoronida TRs, respectively.
Fig. 5. Functional validation and characterization of apple snail TR. (A) Direct telomerase activity assay of acorn worm telomerase reconstituted in vitro. (top) Template sequence of acorn worm TR (open box) with base-pairing of six permuted telomeric DNA primers aâf. Sequence and number of expected nucleotides added are depicted for each primer. (bottom) Direct primer-extension assay of apple snail telomerase. Apple snail telomerase was in vitro reconstituted from T7 transcribed PdiTR (408ânt) and PdiTERT synthesized in RRL. The reconstituted acorn worm telomerase was analyzed with six permuted telomeric DNA primers (lanes 1â6). A 32P end-labeled 18-mer oligonucleotide was added to each reaction as recovery control (r.c.) prior to ethanol precipitation of DNA products. Numbers to the right of the gel indicate the number of repeats or nucleotides added to the primer. (B) Two essential fragments of PdiTR. The T-PK and CR4/5 fragments of PdiTR were synthesized separately and assembled with PdiTERT in RRL, followed by telomerase activity assay. The schematic secondary structures of the PdiTR T-PK (top) and eCR4/5 (bottom) fragments. The eCR4/5 domain consists of three stems, P4âP6. Nucleotide numbers denote the 5â²- and 3â²-ends of the T7 transcribed TR fragments, T-PK, P4/5/6, P5, and P6. (C) Minimal requirement of TR domains for telomerase activity. T7 transcribed PdiTR fragments, T-PK (nt 1â146), P4/5/6 (nt 162â335), P5 (nt 191â244), or P6 (nt 245â309), were assembled with in vitro synthesized PdiTERT and analyzed for activity using the primer (GGGTTA)3. The PdiTR fragments included in each reaction are indicated above the gel. A 32P end-labeled 18-mer oligonucleotide was added to each reaction prior to ethanol precipitation of DNA products as r.c. The number of nucleotides added to the primer is shown to the right of the gel. (D) The effect of P1.1 position on template boundary definition. Two PdiTR mutants, L1 and L2, with a single adenosine residue inserted immediately after positions 11 and 32, respectively, were assembled in vitro with PdiTERT and assayed for telomerase activity using primer (GGGTTA)3.
Fig. 6. Secondary structures of basal metazoan TRs. (A) Evolutionary relationship of basal metazoan phyla and classes with TR identified. The numbers of TRs identified in this study for each class are indicated. An asterisk denotes the classes for which the TR secondary structure of a representative species is shown in BâD. Fungi represents the outgroup in the phylogenetic tree. (BâD) Representative TR secondary structures determined by phylogenetic comparative sequence analyses are shown for Acropora digitifera (digitate coral) from phylum cnidaria (B), Trichoplax adhaerens (flat animal) from phylum placozoa (C), and Amphimedon queenslandica (common sponge) from phylum porifera (D). The conserved TR structural domains are shaded. Universal covariations (thick lines), invariant residues, and residues with >80% conservation are based on the sequence alignment of 55 previously identified animal TRs and 82 novel metazoan TRs identified in this study. Group-specific covariations are indicated and based on the sequence alignment of TRs from individual groups including 26 cnidaria, 1 placozoa, and 3 porifera species.
Fig. 7. Monophyletic origin and divergent evolution of TR structural elements and biogenesis pathways. Simplified phylogenetic trees of major metazoan (A) and eukaryotic (B) lineages shown. Branch lengths do not represent evolutionary distance. (A) Evolution of TR structural elements across metazoan lineages. Loss of RNA secondary structural elements P1.1, P2.1, or P6.1 is indicated in oval. Gain of P2a or P2a.1 is shown in rectangle. Specific events of loss and gain of TRÂ structural elements are marked along the respective metazoan lineages. (B) Divergence of TR biogenesis pathway and transcription machinery across eukaryotes. The TR transcription machinery in each eukaryotic kingdom is indicated as Pol II or Pol III. The size range of TRs from each group is indicated. Biogenesis pathway of respective eukaryotic clades are shown to the right.
Altschul,
Basic local alignment search tool.
1990, Pubmed
Altschul,
Basic local alignment search tool.
1990,
Pubmed
Barquist,
Studying RNA Homology and Conservation with Infernal: From Single Sequences to RNA Families.
2016,
Pubmed
Blackburn,
Telomerase: an RNP enzyme synthesizes DNA.
2011,
Pubmed
Bley,
RNA-protein binding interface in the telomerase ribonucleoprotein.
2011,
Pubmed
Casacuberta,
Drosophila: Retrotransposons Making up Telomeres.
2017,
Pubmed
Chen,
Functional analysis of the pseudoknot structure in human telomerase RNA.
2005,
Pubmed
Chen,
Secondary structure of vertebrate telomerase RNA.
2000,
Pubmed
Chen,
Determinants in mammalian telomerase RNA that mediate enzyme processivity and cross-species incompatibility.
2003,
Pubmed
Chen,
Template boundary definition in mammalian telomerase.
2003,
Pubmed
Chen,
A critical stem-loop structure in the CR4-CR5 domain of mammalian telomerase RNA.
2002,
Pubmed
Dandjinou,
A phylogenetically based secondary structure for the yeast telomerase RNA.
2004,
Pubmed
Egan,
Specificity and stoichiometry of subunit interactions in the human telomerase holoenzyme assembled in vivo.
2010,
Pubmed
Greider,
A telomeric sequence in the RNA of Tetrahymena telomerase required for telomere repeat synthesis.
1989,
Pubmed
Gunisova,
Identification and comparative analysis of telomerase RNAs from Candida species reveal conservation of functional elements.
2009,
Pubmed
Gupta,
The Trypanosoma brucei telomerase RNA (TER) homologue binds core proteins of the C/D snoRNA family.
2013,
Pubmed
Hemmrich,
Compagen, a comparative genomics platform for early branching metazoan animals, reveals early origins of genes regulating stem-cell differentiation.
2008,
Pubmed
Hoover,
DNAWorks: an automated method for designing oligonucleotides for PCR-based gene synthesis.
2002,
Pubmed
Huang,
Structural basis for protein-RNA recognition in telomerase.
2014,
Pubmed
Jády,
Human telomerase RNA and box H/ACA scaRNAs share a common Cajal body-specific localization signal.
2004,
Pubmed
Jansson,
Structural basis of template-boundary definition in Tetrahymena telomerase.
2015,
Pubmed
Jiang,
Structure of Tetrahymena telomerase reveals previously unknown subunits, functions, and interactions.
2015,
Pubmed
Jiang,
Structure of Telomerase with Telomeric DNA.
2018,
Pubmed
Lai,
Template boundary definition in Tetrahymena telomerase.
2002,
Pubmed
Leonardi,
TER1, the RNA subunit of fission yeast telomerase.
2008,
Pubmed
Li,
Identification of purple sea urchin telomerase RNA using a next-generation sequencing based approach.
2013,
Pubmed
,
Echinobase
Mason,
Stem-loop IV of tetrahymena telomerase RNA stimulates processivity in trans.
2003,
Pubmed
Mitchell,
A box H/ACA small nucleolar RNA-like domain at the human telomerase RNA 3' end.
1999,
Pubmed
Mitchell,
Human telomerase activation requires two independent interactions between telomerase RNA and telomerase reverse transcriptase.
2000,
Pubmed
Musgrove,
New perspectives on telomerase RNA structure and function.
2018,
Pubmed
Nawrocki,
Infernal 1.1: 100-fold faster RNA homology searches.
2013,
Pubmed
Nguyen,
Cryo-EM structure of substrate-bound human telomerase holoenzyme.
2018,
Pubmed
Podlevsky,
The functional requirement of two structural domains within telomerase RNA emerged early in eukaryotes.
2016,
Pubmed
Podlevsky,
It all comes together at the ends: telomerase structure, function, and biogenesis.
2012,
Pubmed
Podlevsky,
Evolutionary perspectives of telomerase RNA structure and function.
2016,
Pubmed
Podlevsky,
Structure and function of echinoderm telomerase RNA.
2016,
Pubmed
,
Echinobase
Qi,
The common ancestral core of vertebrate and fungal telomerase RNAs.
2013,
Pubmed
Reichow,
The structure and function of small nucleolar ribonucleoproteins.
2007,
Pubmed
Seto,
Saccharomyces cerevisiae telomerase is an Sm small nuclear ribonucleoprotein particle.
1999,
Pubmed
Shay,
Telomeres and telomerase: three decades of progress.
2019,
Pubmed
Song,
The conserved structure of plant telomerase RNA provides the missing link for an evolutionary pathway from ciliates to humans.
2019,
Pubmed
Stemmer,
DNA shuffling by random fragmentation and reassembly: in vitro recombination for molecular evolution.
1994,
Pubmed
Tesmer,
Two inactive fragments of the integral RNA cooperate to assemble active telomerase with the human protein catalytic subunit (hTERT) in vitro.
1999,
Pubmed
Theimer,
Structure of the human telomerase RNA pseudoknot reveals conserved tertiary interactions essential for function.
2005,
Pubmed
Theimer,
Structural and functional characterization of human telomerase RNA processing and cajal body localization signals.
2007,
Pubmed
Tseng,
The H/ACA complex disrupts triplex in hTR precursor to permit processing by RRP6 and PARN.
2018,
Pubmed
Tzfati,
Template boundary in a yeast telomerase specified by RNA structure.
2000,
Pubmed
Venteicher,
A human telomerase holoenzyme protein required for Cajal body localization and telomere synthesis.
2009,
Pubmed
Vítková,
The evolutionary origin of insect telomeric repeats, (TTAGG)n.
2005,
Pubmed
Wang,
Architecture and assembly of mammalian H/ACA small nucleolar and telomerase ribonucleoproteins.
2004,
Pubmed
Witkin,
Holoenzyme proteins required for the physiological assembly and activity of telomerase.
2004,
Pubmed
Xie,
Structure and function of the smallest vertebrate telomerase RNA from teleost fish.
2008,
Pubmed