The protein subunit of telomerase displays patterns of dynamic evolution and conservation across different metazoan taxa
© The Author(s). 2017
Received: 5 January 2017
Accepted: 4 April 2017
Published: 26 April 2017
Most animals employ telomerase, which consists of a catalytic subunit known as the telomerase reverse transcriptase (TERT) and an RNA template, to maintain telomere ends. Given the importance of TERT and telomere biology in core metazoan life history traits, like ageing and the control of somatic cell proliferation, we hypothesised that TERT would have patterns of sequence and regulatory evolution reflecting the diverse life histories across the Animal Kingdom.
We performed a complete investigation of the evolutionary history of TERT across animals. We show that although TERT is almost ubiquitous across Metazoa, it has undergone substantial sequence evolution within canonical motifs. Beyond the known canonical motifs, we also identify and compare regions that are highly variable between lineages, but show conservation within phyla. Recent data have highlighted the importance of alternative splice forms of TERT in non-canonical functions and although animals may share some conserved introns, we find that the selection of exons for alternative splicing appears to be highly variable, and regulation by alternative splicing appears to be a very dynamic feature of TERT evolution. We show that even within a closely related group of triclad flatworms, where alternative splicing of TERT was previously correlated with reproductive strategy, we observe highly diverse splicing patterns.
Our work establishes that the evolutionary history and structural evolution of TERT involves previously unappreciated levels of change and the emergence of lineage specific motifs. The sequence conservation we describe within phyla suggests that these new motifs likely serve essential biological functions of TERT, which along with changes in splicing, underpin diverse functions of TERT important for animal life histories.
KeywordsTelomerase TERT Evolution Alternative splicing Planarian Metazoa
In linear eukaryotic chromosomes, telomeres function as protective nucleoprotein structures that guard chromosome ends against the action of nucleases, DNA damage signalling and nucleolytic degradation [1–3]. Telomeric DNA, consisting entirely of guanine-rich tandem repeats, is maintained by the action of the enzyme telomerase. The activity of telomerase was first eluded with the framing of the “end replication problem” by Watson  and Olovnikov , which noted that the lagging strand nature of DNA replication would lead to progressive shortening of linear eukaryotic chromosomes. In humans, cells in culture have a finite replicative lifespan limiting the number of times a cell can divide, and eventually senesce . This has since been shown to be a key mechanism in protecting against cancer and may also contribute to causing a decline in physiological functions and organismal ageing [7–12]. While together the shortening of chromosome ends and the observation of the “Hayflick limit” led to the notion that telomeres could act as timers of cellular replicative age, it also raised the question as to how chromosomes maintain their ends across generations, through reproduction, development and during the lifespan of long-lived individuals [2, 13–15]. This led to intensive research on both the structure of chromosome ends and a search for the cellular activities that could maintain these structures in light of the “end replication problem”. Eventually, the ribonucleoprotein telomerase was discovered in the ciliated protozoa Tetrahymena thermophila [16–20]. Soon after, the catalytic subunit of this enzyme, telomerase reverse transcriptase (TERT) was also discovered and then studied in human cells [21–27].
However, despite telomeres and telomerase providing the solution for the existence of linear chromosomes, it is clear that the huge variation in life histories, particularly amongst metazoans, would have different demands on an end chromosome maintenance system. For example, short-lived and long-lived metazoans would have very different needs for telomerase and telomere maintenance. This is exemplified by comparing ourselves with single celled protists, where human life history in an evolutionary context means that signalling from chromosome ends must achieve a balance between allowing growth and preventing inevitable somatic mutations from causing cancer and death before reproductive age [28, 29]. Subsequent pleiotropic effects of signalling mechanisms associated with telomeres now appear to be central to the human ageing process and age associated pathologies, including processes that are unconnected to telomere length [30–33]. On this basis we would expect to observe telomere biology adaptations to reflect life history specific requirements and therefore variation between different groups. For instance, telomerase levels in the long-lived rodent, the naked-mole rat show stable expression across different life stages . In long-lived birds, telomerase activity is correlated with proliferative potential and the demands of specific organ types, for example birds with long lifespans have high levels of telomerase in their bone marrow . In other species telomerase does not appear to be related to homeostasis and longevity, for example studies in laboratory mice have demonstrated that ageing is unlikely to be caused by telomere shortening alone as aged mice have long telomeres and loss of telomerase has no immediate effect on ageing [11, 36–39]. Unlike mammals that stop growing upon reaching adulthood, certain animals with indeterminate growth such as some fish [40–42], bivalve molluscs , sea urchins , lobsters  and planarians  express telomerase through the life cycle.
In addition to its canonical role in telomere maintenance, TERT, the protein subunit of the telomerase enzyme has also been shown to have non-canonical functions [47–50]. Not much is known about whether these functions are conserved or lineage specific. Additionally, TERT transcripts in many animal species, including vertebrates, insects, nematodes and planarians are alternatively spliced [46, 51–55] and spliced variants have also been shown to have non-canonical roles, for example in cell proliferation [56, 57]. TERT splicing dysregulation is present in many cancers and it is thought that human TERT (hTERT) splicing could be manipulated for therapeutic purposes [58, 59]. While only a little is known about TERT alternatively spliced (AS) variants in animals, this represents a potential evolutionary mode for lineage specific life history adaptation, and thus broader assessment of TERT alternate splicing is required particularly in lineages spanning short evolutionary time frames.
Given the importance of telomerase function in lineage-specific life history traits, it is surprising that an in depth understanding of the evolutionary history of TERT in the Metazoa is currently lacking from the literature. Here, we make a detailed study of TERT across the Animal Kingdom using extant data to look for evidence (or lack) of patterns of TERT protein and alternative splicing evolution. We identify TERTs across animal phylogeny and perform within and between phyla comparison of the structure of this gene. We are able to confirm the domain content of the ancestral TERT at the base of the Animal Kingdom and then we are able to show dynamic evolution of the protein in different lineages. This includes both the loss of canonical motifs and the invention of lineage specific domains, for example the invention of novel C-terminal extension and N-terminal linker motifs specific to the vertebrate lineage. In some cases we hypothesized that novel motifs may simply compensate for lost canonical motifs, and in others they may relate to lineage specific activities and interactions that remain to be discovered. Next we use available data to assess the levels of alternative splicing, implicated in non-canonical TERT function in humans and described these in different taxa. We find evidence that alternative splicing is likely to be common feature of some lineages, but is not broadly conserved with respect to splicing patterns. Finally, we study the evolution of splicing in one particular related group of animals, the free-living Tricladida planarians [60–62]. We find that it evolves relatively rapidly in this group of highly regenerative animals, suggesting that study of alternative splicing over shorter evolutionary timescales may be required to understand adaptive non-canonical roles of TERT. Taken together our work highlights a previously unappreciated evolutionary and likely functional diversity in this core eukaryotic protein. We conclude that telomere biology, a core synapomorphy of eukaryotes, has undergone significant changes in different lineages and some of this is through the evolution of the TERT protein subunit of telomerase itself.
Results and discussion
TERT distribution across the animal kingdom and reconstruction of the ancestral TERT domain structure
TERT is found in unicellular relatives of metazoans and early-branching metazoans
Metazoans are related to unicellular lineages from the holozoan clade: filastereans, ichthyosporeans and choanoflagellates [63–67]. While we could identify TERT from Capsaspora owczarzaki, we could not find evidence for TERT orthologs in either the choanoflagellates Monosiga brevicollis or Salpingoeca rosetta. Despite this, both M. brevicollis  and S. rosetta  have the 5′-TTAGGG(n)-3′ telomeric repeat. Based on this finding, a canonical TERT is predicted to be present . The observed absence could reflect incomplete genome sequences or, more likely given the absence from two different choanoflagellate genomes, that TERT primary sequence has evolved beyond detection or a related reverse transcriptase activity has instead taken up the role of TERT. From the genome sequences of other unicellular parasitic protozoans, putative TERT sequences have been predicted [70–77] and our phylogenetic analyses confirms that these predictions are all canonical TERT proteins, so it seems likely that loss of a recognisable TERT in choanoflagellates is a derived character rather than the alternative hypothesis of independent evolution of TERT like proteins in different eukaryotic lineages (Fig. 1). TERT sequences from early-branching metazoans reported to date were from the genomes of Hydra magnipapillata and Nematostella vectensis [78, 79]. The genomes of the sponge Amphimedon queenslandica , comb jelly Mnemiopsis leidyi , myxosporean Thelohanellus kitauei  and sea anemone Aiptasia pallida  have been sequenced and from these genomes we identified new TERT sequences for comparative analysis (Fig. 1; Additional file 1: Table S2). While TERT proteins have been identified in many bilaterian phyla [21–25, 34, 46, 52, 55, 84–94], we additionally identified novel TERT sequences that have not been reported elsewhere. These include TERT orthologs from nematodes species of the classes Secernentea, Chromadorea, Enoplea and Rhabditea and multiple orthologs from Insecta, Actinopterygii, and Mammalia (Additional file 1: Table S2). We also identified novel TERT sequences from parasitic platyhelminthes from the Trematoda class, from the newly sequenced tardigrade genomes of Hypsibius dujardini  and Ramazzottius varieornatus  and the crustacean Parhyale hawaiensis  (Additional file 1: Table S2).
TERT is absent from the bdelloid rotifer
Despite the ubiquity of TERT amongst metazoans, its presence is not universal. TERT has been reported to be absent from some animal species; one prominent example being the loss of TERT in Diptera before divergence from Siphonaptera . Like in dipterans, we failed to identify TERT in the bdelloid rotifer Adineta vaga, suggesting the loss of TERT from this group of long-term asexual species [99, 100]. In the absence of TERT and canonical telomere repeats, Drosophila melanogaster uses a retrotransposon-based telomere maintenance in which telomere-specific retrotransposons are reverse transcribed onto chromosome ends . Three non-Long Terminal Repeat (LTR) retrotransposons, HeT-A [102, 103], TART [104, 105] and TAHRE  are specifically used for telomere maintenance in this species. The bdelloid chromosome ends also have specialised telomere-associated retrotransposons, such as the Hebe transposon  and sub-telomeric reverse transcriptase (RT)-related genes . Penelope-like elements, which lack the endonuclease domain and are located in sub-telomeric regions and are found in two bdelloid rotifer species A. vaga and Philodina roseola . Although no evidence yet directly links non-LTR transposons to end chromosome maintenance in A. vaga, it seems plausible that in the absence of TERT, bdelloid rotifers may employ retrotransposon-based mechanisms for telomere maintenance.
The ancestral structure of TERT in metazoans is likely to possess a complete complement of canonical motifs
Taxa-specific divergence of canonical TERT motifs
We next examined evolution of TERT canonical motifs in detail across taxa. Sequence analyses of 11 TERT canonical motifs revealed that the TEN (GQ motif) and TRBD (CP, QFP and T motifs) domains have high sequence conservation across their whole lengths (Additional file 2: Figure S1).
TERT sequence divergence in unicellular relatives of metazoans
Amongst unicellular eukaryotes considered here, only C. owczarzaki TERT has all 11 motifs (Fig. 2). An independent study reported that TERT in Leishmania amazonensis also contained all 11 motifs . However, from our sequence alignment analyses, we were unable to identify clearly defined GQ motifs in Leishmania sp., Trypanosoma sp., Theileria sp., Babesia sp. and Cryptosporidium sp. (Additional file 3: Figure S2). The observation of GQ in L. amazonensis might be an alignment artefact because from their alignment, only eight out of 244 amino acid residues are found to be in common with other TERTs . We also failed to find the CP motif in this group of unicellular eukaryotes, except for C. owczarzaki and Cryptosporidium sp. (Additional file 3: Figure S2). It was reported Giardia lamblia lacks the T motif . However, this could be due to the use of evolutionarily distant TERT sequences from plants, Caenorhabditis elegans and yeast, for the alignments in their earlier study. We retrieved three TERT sequences for the parasitic protozoan G. lamblia, obtained from different animal isolates . Sequence alignment revealed that they share only 72% identity. From an alignment of these G. lamblia TERTs with other protozoans, we were able to identify a putative T motif, suggesting that it hasn’t been lost. The G. lamblia T motif, however, has eight extra residues not present in other unicellular eukaryotes (Additional file 3: Figure S2). Apart from C. owczarzaki and G. lamblia, it appears that the TEN domain has diverged substantially in unicellular eukaryotes with poor conservation to the canonical GQ motif. Within kinetoplastid protozoans from the genera Leishmania and Trypanosoma, we identified a conserved block of residues, (I/V)QQRVxLQF, between the QFP and T motif. Considering that the divergence time of Trypanosoma sp. from other trypanosomatids including Leishmania sp. is well over 300 million years ago (m.y.a) , it is likely this block of residues in TERT is an example of a lineage specific motif.
TERT sequence divergence in early branching metazoans
Amongst early branching metazoans, TERT sequences from the comb jelly M. leidyi and two cnidarians, A. pallida and the myxozoan T. kitauei appear to be incomplete because their sequences do not extend beyond reverse transcriptase canonical motifs D, B and A respectively (Fig. 2). TERT from the sponge A. queenslandica appears to be complete but entirely lacks the E motif. The GQ motif is collectively absent from the comb jelly and cnidarians T. kitauei and A. pallida, but not from H. vulgaris, indicating that multiple independent losses of GQ occurred amongst early branching metazoans (Fig. 2).
TERT sequence divergence in bilaterians
Within Deuterostomia, vertebrates and Strongylocentrotus purpuratus have retained all 11 conserved motifs while Ciona sp. lack the GQ motif (Fig. 2) . The canonical GQ motif has also been lost from protostomes, except from the nematode lineage (Fig. 2, Additional file 4: Figure S3). In addition to the free-living Caenorhabditis sp. , we find that other nematodes have also retained the GQ motif (Additional file 4: Figure S3). As previously reported we find that the early branching nematode Trichinella spiralis TERT has all 11 canonical TERT motifs , so the ancestor of nematode likely possessed all canonical motifs. Studies on human cell lines have shown that within the GQ motif, an additional region known as ‘dissociates activities of telomerase’ (DAT), is essential for TERT activity . Although nematodes have retained the GQ motif, the DAT region within this motif appears to have lost multiple important residues (Additional file 4: Figure S3). The DAT region was proposed to be involved in other aspects of in vivo telomere elongation because mutations in this region have no direct effect on human TERT multimerization, nuclear targeting or template binding . It appears that either nematodes do not have sensu stricto DAT regions, or that the GQ motif in protostomes, including nematodes, have diverged significantly that sequence conservation with the deuterostome GQ can no longer be detected (Additional file 4: Figure S3). The GQ motif in yeast, human cell lines and T. thermophila plays a role in repeat addition processivity [9, 114–117]. It remains to be determined, however, whether the loss or sequence divergence of GQ in protostomes bears any functional implications in species-specific life histories. From the crystal structure of the red flour beetle Tribolium castaneum TERT, it was shown that although GQ is absent in this species, remaining motifs (CP and T) in the N-terminal region of TERT not only share high degrees of structural conservation with T. thermophila TERT [118, 119] but also conferred in vitro activity . Therefore, even in the absence of a canonical GQ sequence across protostomes (Fig. 2), structural and functional conservation of TERT N-terminal region may persist in these phylogenetic groups. The second TERT canonical motif CP has also been also lost in some protostomes: e.g. in molluscs and tardigrades (Fig. 2). Similarly, most arthropods lack the CP motif, but several exceptions exist: the centipede Strigamia maritima, hymenopterans (ants, sawflies, bees but not wasps) and T. castaneum  have retained the CP motif (Fig. 2; Additional file 5: Figure S4). Other TERT motifs are maintained in metazoan taxa.
We next performed functional divergence analyses on monophyletic clusters from the TERT phylogeny using the DIVERGE programme . A coefficient of functional divergence (θ) value that is significantly different from zero, as indicated by the z-score (MFE z-score) for the model-free method (MFE θ) or the likelihood ratio test (LRT θ) for the maximum likelihood estimate (θ ML), suggests evidence for functional divergence between two clusters. MFE θ and θ ML estimates appeared to be statistically significant for five pairwise comparisons: Euglenozoa/Nematoda (MFE θ = 0.56, P value <0.01, θ ML = 0.55, LRT θ = 16.54), Actinopterygii/Euglenozoa (MFE θ = 0.79, P value <0.01, θ ML = 0.74, LRT θ = 22.63), Mammalia/Crustacea (MFE θ = 0.54, P value <0.01, θ ML = 0.66, LRT θ = 33.55), Mammalia/Euglenozoa (MFE θ = 0.76, P value <0.01, θ ML = 0.76, LRT θ = 22.28) and Mammalia/Apicomplexa (MFE θ = 0.61, P value <0.01, θ ML = 0.69, LRT θ = 35.25; Additional file 6: Table S1). Future experimental validations will be required to investigate these results in detail.
Identification of novel lineage-specific motifs in TERT
TERT proteins from most organisms possess two highly variable regions that are not part of the 11 canonical motifs: the C-terminal extension (CTE) and the N-terminal linker region . The CTE region starts after the last canonical E motif, and despite poor sequence homology of CTEs between human, T. castaneum and T. thermophila, it is thought to constitute the RT ‘thumb’ that is made up of three α-helices crucial in template binding and TERT processivity [118, 122–124]. The mechanistic significance of the CTE has remained elusive until recently, when the CTE was shown to be directly involved in differential binding of DNA and a mutated version of the CTE resulted in defective DNA binding and faster DNA dissociation rates .
CTE regions in unicellular relatives of metazoans and early branching metazoans have regions with high conservation
Since the CTE appears to play an important role despite having poor sequence conservation we examined whether this region has evolved to harbour lineage-specific amino acid conservation. We extracted regions after the last canonical E motif from 96 species representing various taxa and investigated these as potential CTEs. Not all species have CTE regions because TERT proteins in some animals do not extend beyond the E motif. We show that unicellular eukaryotes and two early branching metazoan species, H. vulgaris and A. queenslandica possess CTE regions (Additional file 7: Figure S6A). Alignment of unicellular eukaryotes CTE regions revealed that blocks of conservation exist within different protozoan lineages, although not between unicellular species and human TERT (Additional file 7: Figure S6C). Crytosporidium sp. CTE regions are considerably shorter than those of other protozoans. Kinetoplastid protozoans, Leishmania sp. and Trypanosoma sp., have a block of conserved residues made up of 48 amino acids at the end of their TERT proteins (Additional file 7: Figure S6C). In other early-branching metazoans, although the parasitic cnidarian T. kitauei and the comb jelly M. leidyi lack CTE regions, we identified some conserved residues in CTE regions when aligning these regions from A. queenslandica and H.vulgaris (Additional file 7: Figure S6A).
Taxa specific conservation of CTE regions exist in bilaterians
Within bilaterians, it has been reported that nematodes from the Caenorhabditis genus have some of the shortest TERT proteins and appear to have lost their CTE structure [52, 111]. We observe that this is also true for the parasitic roundworm Ancylostoma ceylanicum where the C-terminal region beyond the E motif only has 42 amino acid residues (Additional file 7: Figure S6E). We discovered that most parasitic nematodes from the Chromadorea and Secernentea classes possess intact CTE regions (Additional file 7: Figure S6E). In addition, we show that filarial nematodes B. malayi, O. volvulus and L. loa share a RIAVL(R/K)FLKASLLEKYR(M/V motif (Additional file 7: Figure S6E). The observation that a CTE is present in most parasitic nematodes but not the free-living Caenorhabditis sp. may reflect adaptation of TERT to their very different life histories; i.e. some parasitic nematodes have long generation times and life expectancies and are able to survive for years in their hosts . CTE structures in Arthropoda varied markedly and we performed separate alignments for insects and crustaceans (Additional file 7: Figure S6D and F). We identified two crustacean specific motifs within the CTE region RL(K/Q)x(I/V) and R(L/F)xAL (Additional file 7: Figure S6F). Within insects, ants lack the CTE completely as TERT proteins within this lineage do not extend beyond the E motif (Additional file 7: Figure S6D). Lepidopterans have additional residues at the extreme ends of their TERT proteins conserved within this lineage only; a lepidopteran-specific CTE contains 32 to 36 residues (Additional file 7: Figure S6D). Within Deuterostomia, we annotated highly conserved CTE motifs in tetrapods and fishes (Additional file 7: Figure S6G).
Our findings suggest that the ancestor of metazoans possessed a CTE and that while this feature has been largely retained in most animal lineages, it has been lost multiple times. We have demonstrated that these CTE regions can be incredibly diverse across animal phyla, and we show that novel motifs exist in different lineages. It seems likely that these lineage specific motifs in the CTE region would have significant roles, since functional studies in T. thermophila, yeast and humans have revealed the importance of CTE for TERT catalytic activity related to enzyme multimerization, processivity and DNA binding [124, 125, 127–129].
Conserved residues within the N-terminal linkers of unicellular relatives of metazoans and early-branching metazoans
We next investigated the N-terminal hypermutable linker (hereafter referred to as the linker) of TERT that is defined as the region between the GQ motif and CP motif. The linker region, although evolutionarily divergent, has been shown to be biologically essential for the function of TERT. In vertebrates, the association between TERT and telomerase RNA (TR) is mediated by the VSR motif in the linker region that binds the activation domain of TR [130, 131]. The puffer fish Takifugu rubripes has an additional short alpha-helix upstream of the VSR motif, known as the TFLY motif that is also implicated in binding the template boundary element [119, 132], which is a precise region within the RNA template that is defined by structural elements identified in T. thermophila, Kluyveromyces lactis and mammals [133–136]. Amongst early-branching metazoans, we found pockets of conservation in the linker region; the motif LxxAIF is present in three out of five early-branching metazoan species (Additional file 8: Figure S7A). Unicellular eukaryotes have very divergent N-terminal linkers between phyla but some conservation within each protozoan lineage, the kinetoplastids (Leishmania sp. and Trypanosoma sp.) have the ALxR(T/L)(D/N)V(P/S)RxxL motif and the piroplasmids (Babesia sp. and Theileria sp.) have the LLxGLF motif (Additional file 8: Figure S7C).
Vertebrates have a distinct VL region within their N-terminal linkers
We performed an alignment of the linker regions from representative deuterostome phyla and observed that linker regions have poor conservation and they do not align well. Separately, an alignment of linkers from vertebrates alone revealed that vertebrates share a well-conserved block of 91 amino acid residues directly upstream of the CP motif, which we named the Vertebrata-specific N-terminal linker (VL; Additional file 8: Figure S7F). Within this VL region we noticed that class specific motifs exist. According to the molecular timescale for vertebrate evolution, mammals diverged from birds 310 m.y.a. and marsupials from placental mammals at 173 m.y.a. . Teleosts and tetrapods diverged approximately 450 m.y.a. after the vertebrate genome duplication . Such long evolutionary distances enabled us to distinguish conserved motifs from neutrally evolving sequences. For example, we demonstrated that ray-finned fishes have the FIRTLGFLY and RRxQGxD motif, birds have the NQSL motif and marsupials have the HLFxxKGDPxQQ motif (Additional file 8: Figure S7F). It seems highly likely that these conserved motifs have functional roles that will be important for TERT function in each of these lineages.
Protostomes have truncated N-terminal regions
Extending our analyses to protostomes, we observed that protostomes generally have truncated N-terminal regions. Because arthropods and lophotrochozoans lack the GQ motif altogether, we extracted sequences upstream of the CP motif (for species that have CP motifs) or the QFP motif (for species without the CP motif) as putative linker regions (Additional file 8: Figure S7B, D and E). Insects, including ants, bees and lepidopterans exhibit long stretches of conserved residues (Additional file 8: Figure S7D). Crustaceans from the classes Malacostraca and Branchiopoda share both FPxxHILS and ILxxNxG motifs in their N-terminal linker regions (Additional file 8: Figure S7A). Amongst nematodes, most species with Trichuris suris as the exception (Fig. 2), have truncated N-termini upstream of the CP motif. Nevertheless, some conservation exists within this region, for instance nematodes share a IGxxNxxxx(L/V) motif (Additional file 5: Figure S4A). Overall, although the functional significance of the N-terminal linker is unclear beyond vertebrates, it is apparent that the linker sequence can provide additional surface areas to allow conformational flexibility between the N-terminus and the rest of the TERT protein . It is likely that the conserved linker modifications within different lineages confer some but yet unknown essential biological functions.
TERT exhibits conservation of intron-exon structure and pan-metazoan regulation by alternative splicing
Intron-exon structure of TERT in metazoans
Evidence of alternative splicing of TERT in bilaterians
Introns provide the opportunity for creating new gene products via alternative splicing mechanisms. TERT alternatively spliced products have been implicated in non-canonical functions in human disease processes, particularly cancer and in adaptation to reproductive mode in planarians flatworms . Alternative splicing has been reported in many other taxa [51–54, 89, 92, 93, 141–144]. We sought to determine whether alternative splicing of TERT is common across animals and if AS variants show any conservation. We probed available genomic and expression datasets for the presence of AS variants. We identified 19 animal species with TERT AS variants as evident by computational predictions using RNA-sequencing datasets (Additional file 9: Table S3; Additional file 10). We found AS variants in both protostomes (10 species) and deuterostomes (9 species). From this list, AS variants have been previously described for these five species: Homo sapiens [51, 56, 141, 144], Ornithorhynchus anatinus , Gallus gallus [54, 142], Danio rerio [89, 143] and Solenopsis invicta . We observed extensive patterns of alternative splicing in hymenopteran and hemipteran insects (Additional file 11: Figure S8). Some of the most common splicing events in hymenopterans and hemipterans involved the exclusion of exon 2 (6 out of 9 species) and exon 3 (4 out of 9 species; Additional file 11: Figure S8). All 43 AS variants in insects have retained an open reading frame (ORF; Additional file 12: Figure S9). Of these, 11 have truncated TRBD (motifs CP, QFP and T) domains and nine have truncated RT (motifs 1, 2, A, B′, C, D and E) domains (Additional file 11: Figure S8). The substantial lengths and complete ORFs of many AS variants in insects indicate that they may still be functional at the protein level. In Lophotrochozoa, two AS variants were identified in the mollusc A. californica (Additional file 12: Figure S9). Since the combination of exons 2 and 3 only encode 39 amino acid residues, all AS variants in A. californica appear to retain a complete set of motifs like the full-length TERT (Additional file 12: Figure S9). Within Arthropoda, aphids and lepidopterans lack any introns and hence alternative splicing (Fig. 3). The short-lived nature of these taxa may preclude the need for highly regulated TERT activity.
Amongst vertebrate species, many AS variants have been previously described but to our knowledge, splicing patterns of TERTs in the Amazon molly Poecilia formosa, the green spotted puffer fish Tetraodon nigroviridis, the wild turkey Meleagris gallopavo and the Rhesus macaque Macaca mulatta have not been reported (Additional file 13: Figure S10). AS variants in hTERT involved exons 6 to 8 and we witness the splicing events concerning these exons in five out of nine species (3 mammals, the wild turkey and zebrafish) (Additional file 13: Figure S10), suggesting some conservation of splicing across the vertebrates. The exclusion of the unusually large exon 2 (1.3 kb) is a common 5′ splice variant in primate TERTs  and we observe this splicing event in all vertebrates except the Rhesus macaque (Additional file 13: Figure S10). The lack of observation of the AS variant in Rhesus Macaque may be due to lack of confirmatory data. Overall, many AS variants in mammals and fishes appear to be of substantial lengths, suggesting that they would encode functional proteins given that most of their TRBD and RT motifs are still present (Additional file 9: Table S3). In avians however, many AS variants appear to lack any canonical motifs. The wild turkey has 10 AS variants, eight of these involve the exclusion of multiple exons and/or loss of frame mutations, do not contain any TERT motifs and are therefore likely to be non-coding and be subjected to nonsense-mediated decay (Additional file 13: Figure S10B). This observation of complex splicing patterns generating many short AS variants with premature termination codon has also been reported in chickens [54, 142] (Additional file 13: Figure S10B). It is unknown as to whether these AS variants would participate directly or indirectly in canonical telomerase function or in other non-canonical functions of TERT. But clearly alternative splicing of the TERT protein is a pan-metazoan phenomenon. Given that we do not see broad conservation of splicing patterns (except within the vertebrate lineage) our data suggest, that in combination with the lineage specific sequence evolution we have also described, AS is likely to contribute to lineage specific functions of this core eukaryotic protein.
Free-living planarian species exhibit complex patterns of alternative splicing
While our analysis of alternative splicing on TERT revealed extensive evidence of splicing, data available to assess evolutionary changes of splicing in taxa outside of the vertebrates is relative sparse. A previous study had investigated telomere length, telomerase activity and TERT regulation in the highly regenerative model planarian Schmidtea mediterranea . This study found that the asexually reproducing biotype upregulated active full length splice forms of TERT during regeneration to sufficient levels to maintain telomere length during stem cell division. Another study in the related species Dugesia ryukyuensis demonstrated equivalent telomere length maintenance differences between sexual and asexual strains . In order to use planarians as a study system for understanding how AS could contribute to the evolution of TERT regulation we performed a detailed investigation of TERT gene structure, expression and alternative splicing in Tricladida planarians, a group where most, but not all species, are known to be highly regenerative.
Planarian-specific motifs within the N-terminal linkers and CTE regions of TERT
TERT intron-exon structure in planarians and the prevalence of alternative splicing
The observation of multiple PCR products when using primers that span the entire TERT transcript in planarians was suggestive of TERT alternative splicing (Fig. 4b). We wished to investigate the extent of the regulation by alternative splicing in this group of animals that exist over short evolutionary time frames. It was previously reported that TERT in S. mediterranea is alternatively spliced into four variants . Analysis of TERT intron-exon boundaries using existing genome assembly datasets for S. mediterranea  and Girardia tigrina (unpublished data) revealed that Smed-TERT and Gtig-TERT have 16 and 15 exons respectively. Of these, 11 Smed-TERT and 10 Gtig-TERT introns are conserved with hTERT introns where they share the same positions on the TERT protein (Fig. 4f). As planarian TERTs shared a high degree of structural similarity and due to the lack of genomic data for the nine remaining planarian species, we based our AS variant analysis on an inferred genomic structure generated from the Smed-TERT sequence. From the splicing patterns of cloned AS variants, we could visualise exons that are being excluded and this confirmed that the intron-exon boundaries were accurately inferred based on Smed-TERT sequence (Additional file 17: Figure S13; Additional file 18: Figure S14; Additional file 19). AS variants were identified from all seven Dugesiidae species, two Planariidae species and one Dendrocoelidae species Additional file 17: Figure S13; Additional file 18: Figure S14A and B). We did not find any AS variants from the marine planarian Procerodes plebeja (Procerodidae), an outgroup to all other species (Additional file 18: Figure S14C). The highest number of AS variants was identified in G. tigrina (13 variants) with the lowest being Dendrocoelum lacteum (2 variants; Additional file 17: Figure S13; Additional file 18: Figure S14). Two of the shortest AS variants in G. tigrina did not contain any TERT motifs (Additional file 17: Figure S13). Dugesiidae species in general had more AS variants compared to Dendrocoelidae and Planariidae species. Five AS variants in D. ryukyuensis only have partial QFP motifs and not any other TERT motifs (Additional file 17: Figure S13). Of all 52 AS variants in triclad planarians, the most prevalent pattern is the truncation of the RT domain where 24 out of 52 have partial or missing RT domains (Additional file 17: Figure S13; Additional file 18: Figure S14; Additional file 20: Table S5). We observed 20 out of 52 AS variants containing truncated versions of both TRBD and RT domains and only 7 out of 52 contain a truncated TRBD but intact RT domain (Additional file 20: Table S5; Additional file 17: Figure S13; Additional file 18: Figure S14). Despite sharing conserved positions of introns (as inferred from splicing patterns of each species), there is low degree of conservation in the selection of spliced sites for these species. The most prevalent splicing pattern is the deletion of exons 8 to 10 (represented by orange triangles), which is present in six out of ten species where AS variants were found (Additional file 17: Figure S13; Additional file 18: Figure S14). Deletion of exons 11 and 12 (represented by blue triangles) were identified in five out of ten species. These splicing events would affect the TERT RT domain that is encoded by exons 8 to 13 (Additional file 17: Figure S13; Additional file 18: Figure S14). We note that the least amount of TERT splicing is observed in the least regenerative species, P. plebeja and D. lacteum (Additional file 21: Figure S15; Additional file 22: Figure S16) . The prevalence of TERT splicing and regenerative capacity is broadly correlative at this stage and future functional studies will be required to ascertain whether any direct linkage between splicing and regenerative capacity exists. Overall our data suggest that AS splicing of TERT is evolving rapidly within this group of animals and that both fine grained study of TERT alternative splicing as well as functional study of AS variants are required to understand why TERT has evolved so dynamically during animal evolution.
We performed thorough analyses of TERT canonical and non-canonical motifs, gene structure and alternative splicing across representative metazoan species. We show that although the ancestral TERT protein in Metazoa is likely to possess all 11 canonical motifs, the GQ and CP motifs are prone to lineage specific losses (Fig. 2). Beyond the canonical motifs, we demonstrate that the N-terminal linkers and CTE regions of TERTs are highly divergent across phyla. In these regions, we discovered novel motifs that exhibit high levels of sequence conservation over long evolutionary times indicating that they may serve an unknown but important biological function. For example, the CTE regions are distinct within phyla (Additional file 7: Figure S6) and it seems likely that such differences are reflective of life history diversity as growing evidence points to the importance of the CTE for telomerase processivity in vivo [123, 124, 157]. We also demonstrated that the N-terminal linker regions have poor sequence conservation between metazoan phyla (Additional file 8: Figure S7). The N-terminal linkers of TERT proteins have been implicated to mediate a conserved function in enhancing the translocation of the 3′ end of DNA substrate relative to the RNA template [158, 159]. Deleting residues in this region reduces the translocation of DNA substrate and overall processivity . Future studies will be needed to unravel the biological significance of phyla specific sequences and whether they fine tune telomerase processivity to life history strategy.
Analysis on TERT intron exon structure revealed that the metazoan ancestor is intron rich and share many intron positions with TERTs in the vertebrate lineage (Fig. 3). TERTs in protostomes have undergone significant intron loss; lepidopterans and aphids do not have any introns, perhaps due to short lifespans of these taxa. Alternative splicing of TERT has been extensively studied in many species. We report new observations on splicing events for 14 metazoan species and show that the selection of spliced exons is poorly conserved in these animals (Additional file 11: Figure S8; Additional file 12: Figure S9; Additional file 13: Figure S10). Lastly, we investigated TERT alternative splicing in Tricladida planarians and demonstrated that splicing is dynamic and rapidly evolving within this group of closely related species (Fig. 4, Additional file 17: Figure S13; Additional file 18: Figure S14). It is likely that splicing in planarian TERTs is significant for the highly regenerative and potentially immortal life histories of some of these species and AS variants could possess important non-canonical functions. This study highlights that neither TERT structure nor its regulation is static. A future endeavour will be to unravel mechanistic linkages that connect unique sequence evolution and the dynamic regulation of TERT to actual biological roles.
Data mining of TERT sequences
TERT sequences were retrieved from Ensembl Genome and Ensembl Metazoa , NCBI and Uniprot. Sequences were confirmed using reciprocal blastx against the nr database using the NCBI Blast suite. We retrieved tardigrade TERTs from the Hypsibius dujardini genome browser (http://badger.bio.ed.ac.uk/H_dujardini/blast/index) and the Ramazzottius varieornatus gene models ( http://kumamushi.org/database.html). TERT sequence from the crustacean Parhyale hawaiensis was retrieved from the assembled genome . Additional references for other datasets are included in the main text. Complete species names, accession numbers and database sources are listed in Additional file 1: Table S2.
Multiple sequence alignment of TERT protein sequences (from motif T to motif E) was performed using MAFFT with the following parameters: scoring matrix BLOSUM62, gap penalty 1.53, offset value 0.123 . Phylogenetic tree was constructed using RAxML  employing the WAG matrix and gamma distribution rate model with 1000 bootstrap replicates. A single tree topology was generated based on the best-scoring maximum likelihood tree. The tree Figure was made using FigTree (http://tree.bio.ed.ac.uk/software/figtree/). Branches that contained less than 30% boostrap replicates were collapsed.
Annotation of TERT canonical motifs and non-canonical motifs in C-terminal extension (CTE) and N-terminal linker regions
TERT protein motifs (GQ, CP, QFP, T, 1, 2, A, B′, C, D and E) were annotated based on regional comparisons with hTERT protein annotations in multiple sequence alignments performed using MAFFT with the following parameters: scoring matrix BLOSUM62, gap penalty 1.53, offset value 0.123 . For the alignments of canonical domains, we used published domain annotations from two reference sequences (hTERT and C. elegans TERT). For the alignment of TERT domains in unicellular relatives of metazoans (Additional file 3: Figure S2) and insects (Additional file 5: Figure S4), the hTERT sequence was used as a reference while for nematodes, the C. elegans TERT sequence was used as a reference (Additional file 4: Figure S3). Global alignments of full-length TERT sequences were first performed on the target taxa to allow the annotation of regions that align well. Regions that do not align well were extracted and local alignments were then performed with the reference sequence. For example, we identified the GQ motif in unicellular relatives of metazoans from a local alignment of the N-terminal region (Additional file 3: Figure S2). For CTE and linker regions, local alignments were performed only on selected taxa for the identification of non-canonical motifs. Fasta files for all alignments performed in this study are included (Additional files 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52). Alignment Figures were generated using Geneious (version 7) .
All freshwater flatworm strains were cultured at 20 °C using the planarian water formulation. 1X Montjuic salt solution was prepared using milliQ ddH2O with the following composition: 1.6 mM NaCl, 1 mM CaCl2, 1 mM MgSO4, 0.1 mM MgCl2, 0.1 mM KCl, 1.2 mM NaHCO3 . The marine flatworm Procerodes plebeja was cultured at 14 °C using a salt water formulation made with Tropic Marin sea salt to a salinity of 28-30 ppm. Dendrocoelum lacteum was fed with shrimp while all other flatworms used in this study were fed organic beef liver once a week. All flatworms were starved for 1 week prior to any experimental procedures. The worms were kept in the dark at all times apart from feeding and water changing times.
RNA extraction and reverse transcription in Triclad flatworms
Total RNA was isolated from three animals of each species using the TRIzol reagent (Thermo Fisher Scientific) according to manufacturer’s instructions. cDNA was generated from total RNA extract using the QuantiTect reverse transcription kit (Qiagen) containing the RT primer mix (blend of oligo-dT and random hexamers). A genomic DNA elimination step was also included in this kit performed using the gDNA wipeout buffer.
Degenerate polymerase chain reaction (PCR) and 5′ 3′ RACE for the cloning of TERT genes
TERT transcript sequences for Girardia tigrina , Dugesia japonica , D. lacteum, and P. nigra  were retrieved using tblastn with Smed_TERT as a query from published transcriptome sequences of the respective species. Top hits with the best e-value for each species was used for reciprocal blast against the nr database in NCBI to confirm the identification of TERT. For degenerate PCR, first strand cDNA synthesis was performed as mentioned previously. For 5′ and 3′ RACE, first strand cDNA synthesis was performed on TRIzol extracted total RNA using the SMARTer RACE 5′ 3′ Kit (Clontech).
For the cloning of TERT in all the other Dugesia sp. (Dugesia ryukyuensis, Dugesia tahitientis, Dugesia benazzi) degenerate PCR primers were designed based on the pairwise alignment of Smed-TERT and Djap-TERT: DugesiaFD (5′-CACTGGTGYGAATCACCA-3′) and DugesiaRD (5′- AAAACATCATCMACRTATTG-3′). PCR products were cloned in the pGEMT-easy vector (Promega) followed by colony PCR with M13F (5′- GTAAAACGACGGCCAGT-3′) and M13R (5′- GGAAACAGCTATGACCATG-3′) primers. Positive colonies were selected for Sanger sequencing. Once the middle sequence of TERT transcripts for all three Dugesia sp. were obtained, 5′ and 3′ RACE was performed on cDNA made with the SMARTer RACE 5′ 3′ Kit using gene specific primers listed in Additional file 14: Table S4 according to manufacturer’s instructions.
For the cloning of TERT in Schmidtea lugubris, the same set of degenerate PCR primers (DugesiaFD and DugesiaRD) was used for PCR followed by cloning and sequencing. PCR products were cloned in the pGEMT-easy vector (Promega) followed by colony PCR and sequencing as mentioned previously. The 5′ end of Slug_TERT transcript was cloned using the following primer set: Dugesia5F (5′-AATYGAGMGWMATGGTTT-3′; this primer was designed based on 5′ sequence region in Dugesia sp.) and LugR (5′- CTGAAATTTGTGCCATTG-3′; Slug-TERT gene specific primer). Next, 3′ RACE was performed with the SMARTer RACE 5′ 3′ Kit using gene specific primers listed in Additional file 14: Table S4 according to manufacturer’s instructions.
For the cloning of TERT in Polycelis felina, degenerate PCR primers were designed based on pairwise alignment of Pnig-TERT and Polycelis tenuis TERT (obtained from PlanMine) : PolycelisDF (5′-AATTGGCACMTSTTYCTG-3′) and PolycelisDR (5′-GACTCRTARCAAYTCTT-3′). PCR products were cloned in the pGEMT-easy vector (Promega) followed by colony PCR and sequencing as mentioned previously. The 5′ and 3′ end of Pfel-TERT was obtained using 5′ and 3′ RACE with the SMARTer RACE 5′ 3′ Kit using gene specific primers listed in Additional file 14: Table S4 according to manufacturer’s instructions.
For the cloning of TERT in P. plebeja, degenerate PCR primers were designed based on pairwise alignment of Pnig_TERT and Pfel_TERT: PpleDF (5′- ATGSARTATAAAGGATWTAT-3′) and PpleDR (5′-AAYRTCRTCGACATATTG-3′). PCR products were cloned in the pGEMT-easy vector (Promega) followed by colony PCR and sequencing as mentioned previously. The 5′ and 3′ end of Pple-TERT was obtained using 5′ and 3′ RACE with the SMARTer RACE 5′ 3′ Kit using gene specific primers listed in Additional file 14: Table S4 according to manufacturer’s instructions. All degenerate PCRs were performed using the Advantage 2 Polymerase (Clontech) according to recommended thermal cycling parameters. All RACE products were ran on a gel and gel extraction was performed with the MinElute gel extraction kit (Qiagen) following by cloning and sequencing of the RACE products.
DNA extraction, 18S rDNA PCR and sequencing
High molecular weight genomic DNA (gDNA) was isolated using the phenol-chloroform method followed by ethanol precipitation. PCR was performed according to conditions in Carranza et al.  using the Advantage 2 Polymerase (Clontech). List of 18S primer sequences are provided in Additional file 14: Table S4. PCR products were gel extracted using the MinElute gel extraction kit (Qiagen) and eluted in 10uL of molecular grade water. Gel extracted products were sequenced in both directions using 18S nested primers listed in Additional file 14: Table S4. Sequences were aligned with known 18S sequence data available from NCBI (accession numbers provided in Additional file 14: Table S4). 18S rDNA sequences obtained from this study fell into two groups (18S type I and 18S type II).
Cloning of planarian TERT AS variants
PCR primers were designed to amplify from the start to end of TERT transcript sequences obtained from above. Full list of primers used to amplify TERT AS variants are provided in Additional file 14: Table S4. PCR was performed using Phusion Polymerase (Thermo Fisher Scientific) and all PCR products were ran on a gel. In order to cut out gel pieces containing the AS variants, the gel image was overexposed to enable visualization of low abundance variants. Gel regions containing bands were excised and extracted using the MinElute gel extraction kit (Qiagen). A-tailing was subsequently performed with GoTaq Polymerase (Promega) for TA-cloning purposes since Phusion PCR products were blunt ended. PCR products were cloned in the pGEMT-easy vector (Promega) followed by colony PCR with M13F and M13R primers. Colony PCR products were run on a gel for 2 h using 1.3% agarose in 1X TAE buffer to allow good resolution of bands. Colony PCR gel would display an array of bands with different sizes representing the AS variants. Bands of different sizes were selected for Sanger sequencing. For each species, 25 colonies were selected for sequencing to allow exhaustive identification of TERT AS variants.
Sequence analyses of flatworm TERT AS variants
Eleven planarian species were used in this study. Planarians Schmidtea mediterranea, S. lugubris, G. tigrina, D. japonica, D. benazzi, D. tahitientis, P. felina, P. nigra and D. lacteum were purchased from Sciento (http://www.sciento.co.uk/). D. ryukyuensis was a gift from Dr. Midori Matsumoto (Keio University, Japan) and P. plebeja from Dr. Bernhard Egger (Universität Innsbruck). The gene structures (exon-intron boundaries) for Smed-TERT and Gtig-TERT were identified using blastn against the S. mediterranea  and G. tigrina genome assembly (unpublished data). Both Smed-TERT and Gtig-TERT have almost identical intron positions, with Gtig-TERT having one less exon. Since genome data is not available for all the other flatworm species, we based our analyses on the gene structure of Smed-TERT and Gtig-TERT and refer to them as inferred exon positions. Further analyses on cloned AS variants for species without genomic data confirmed the accurate positioning of intron-exon boundaries based on the excluded exons. The full-length TERT sequences for each species were referred to as ‘isoform 1’. Cloned AS variants were mapped to ‘isoform 1’, which set as a reference sequence for each respective species. From this, we were able to visualize the skipped exons, splice site mutations or retained introns by comparing the AS variant sequences to the reference sequence. All the AS variants including the longest isoform were in silico translated using Geneious (version 7)  to the correct open reading frame. Multiple sequence alignments of all AS protein sequences for each respective species were performed using MAFFT with the following parameters: scoring matrix BLOSUM62, gap penalty 1.53, offset value 0.123 . Annotations of TERT canonical motifs were performed according to positions in hTERT; telomerase RNA-binding domain motifs (QFP and T) and canonical reverse transcriptase motifs (1, 2, A, B′, C, D and E) in all AS protein variants. Further information on the flatworm AS variants were provided in Additional file 20: Table S5.
Annotation of TERT AS variants in vertebrates, insects and mollusc
The annotation of TERT AS variants were performed based on publically available transcriptomic and genomic datasets (Additional file 1: Table S2). To evaluate alternatively spliced exons, we have identified 19 species that contained more than one TERT transcripts in the transcriptomes and these additional TERT transcripts were considered as TERT AS variants. These TERT transcripts were in silico translated according to the correct open reading frame into protein sequences. Subsequent analyses were performed as described in the previous section “Sequence analyses of flatworm TERT AS variants”. The longest transcripts/isoforms were used as reference sequences and all other AS variants were mapped to the reference using Geneious (version 7) . This allowed the visualization of skipped exons, splice site mutations or retained introns by comparing the AS variant sequences to the reference sequence. Detailed descriptions on splicing patterns are available in Additional file 9: Table S3.
Annotation of intron positions on TERT protein sequences
Methods for the identification of conserved intron positions were adapted from two previous reports [139, 165]. The locations of the exon boundaries were first annotated on TERT transcripts based on their cognate genomic sequences (schematic diagram in Additional file 53: Figure S17). Next, the transcripts were in silico translated to protein sequences while retaining information on the location of each exon boundary (also referred here as intron positions). TERT protein sequences were independently aligned with the hTERT sequence. Pairwise alignments with hTERT were performed to allow the generation of a schematic diagram that is scaled according to the hTERT sequence. From the pairwise alignments, intron positions that overlap with hTERT introns were defined as conserved introns, while intron positions that do not overlap with hTERT introns were defined as species-specific introns (Additional file 53: Figure S17).
We thank Dr. Hanae Nodono (Kagoshima University, Japan) for her contribution on TERT identification in planarians. We thank Dr. Midori Matsumoto (Keio University, Japan) for the gift of Dugesia ryukyuensis. We thank Dr. Bernhard Egger (Universität Innsbruck) for the gift of Procerodes plebeja.
This work was supported by the Biotechnology and Biological Sciences Research Council (grant number BBK0075641); Human Frontier Science Program Postdoctoral Fellowship (to A.G.L.); University of Oxford Elizabeth Hannah Jenkinson Grant (to A.G.L. and S.S.) and Clarendon Scholarship (to S.S.).
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article as additional files 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53.
AGL and AAA conceived and designed the experiments. AGL, NP-S, AD, GK and SS performed the experiments. AGL and AAA analysed the data. AGL and AAA wrote the paper with help from all the other authors. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Muller HJ. The production of mutations by X-rays. Proc Natl Acad Sci U S A. 1928;14:714–26.PubMedPubMed CentralView ArticleGoogle Scholar
- Harley CB, Futcher AB, Greider CW. Telomeres shorten during ageing of human fibroblasts. Nature. 1990;345:458–60.PubMedView ArticleGoogle Scholar
- De Lange T. Telomere-related genome instability in cancer. Cold Spring Harb Symp Quant Biol. 2005;70:197–204.PubMedView ArticleGoogle Scholar
- Watson JD. Origin of concatemeric T7 DNA. Nat New Biol. 1972;239:197–201.PubMedView ArticleGoogle Scholar
- Olovnikov AM. A theory of marginotomy. The incomplete copying of template margin in enzymic synthesis of polynucleotides and biological significance of the phenomenon. J Theor Biol. 1973;41:181–90.PubMedView ArticleGoogle Scholar
- Hayflick L, Moorhead PS. The serial cultivation of human diploid cell strains. Exp Cell Res. 1961;25:585–621.PubMedView ArticleGoogle Scholar
- Harley CB, Villeponteau B. Telomeres and telomerase in aging and cancer. Curr Opin Genet Dev. 1995;5:249–55.PubMedView ArticleGoogle Scholar
- Rudolph KL, Chang S, Lee HW, Blasco M, Gottlieb GJ, Greider C, et al. Longevity, stress response, and cancer in aging telomerase-deficient mice. Cell. 1999;96:701–12.PubMedView ArticleGoogle Scholar
- Autexier C, Lue NF. The structure and function of telomerase reverse transcriptase. Annu Rev Biochem. 2006;75:493–517.PubMedView ArticleGoogle Scholar
- Blackburn EH, Greider CW, Szostak JW. Telomeres and telomerase: the path from maize, Tetrahymena and yeast to human cancer and aging. Nat Med. 2006;12:1133–8.PubMedView ArticleGoogle Scholar
- Tomás-Loba A, Flores I, Fernández-Marcos PJ, Cayuela ML, Maraver A, Tejera A, et al. Telomerase reverse transcriptase delays aging in cancer-resistant mice. Cell. 2008;135:609–22.PubMedView ArticleGoogle Scholar
- Bernardes de Jesus B, Vera E, Schneeberger K, Tejera AM, Ayuso E, Bosch F, et al. Telomerase gene therapy in adult and old mice delays aging and increases longevity without increasing cancer. EMBO Mol Med. 2012;4:691–704.PubMedPubMed CentralView ArticleGoogle Scholar
- Counter CM, Avilion AA, LeFeuvre CE, Stewart NG, Greider CW, Harley CB, et al. Telomere shortening associated with chromosome instability is arrested in immortal cells which express telomerase activity. EMBO J. 1992;11:1921–9.PubMedPubMed CentralGoogle Scholar
- Vaziri H, Dragowska W, Allsopp RC, Thomas TE, Harley CB, Lansdorp PM. Evidence for a mitotic clock in human hematopoietic stem cells: loss of telomeric DNA with age. Proc Natl Acad Sci U S A. 1994;91:9857–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Bryan TM, Englezou A, Gupta J, Bacchetti S, Reddel RR. Telomere elongation in immortal human cells without detectable telomerase activity. EMBO J. 1995;14:4240–8.PubMedPubMed CentralGoogle Scholar
- Blackburn EH, Gall JG. A tandemly repeated sequence at the termini of the extrachromosomal ribosomal RNA genes in Tetrahymena. J Mol Biol. 1978;120:33–53.PubMedView ArticleGoogle Scholar
- Szostak JW, Blackburn EH. Cloning yeast telomeres on linear plasmid vectors. Cell. 1982;29:245–55.PubMedView ArticleGoogle Scholar
- Greider CW, Blackburn EH. Identification of a specific telomere terminal transferase activity in Tetrahymena extracts. Cell. 1985;43:405–13.PubMedView ArticleGoogle Scholar
- Greider CW, Blackburn EH. The telomere terminal transferase of Tetrahymena is a ribonucleoprotein enzyme with two kinds of primer specificity. Cell. 1987;51:887–98.PubMedView ArticleGoogle Scholar
- Greider CW, Blackburn EH. A telomeric sequence in the RNA of Tetrahymena telomerase required for telomere repeat synthesis. Nature. 1989;337:331–7.PubMedView ArticleGoogle Scholar
- Kilian A, Bowtell DD, Abud HE, Hime GR, Venter DJ, Keese PK, et al. Isolation of a candidate human telomerase catalytic subunit gene, which reveals complex splicing patterns in different cell types. Hum Mol Genet. 1997;6:2011–9.PubMedView ArticleGoogle Scholar
- Harrington L, Zhou W, McPhail T, Oulton R, Yeung DS, Mar V, et al. Human telomerase contains evolutionarily conserved catalytic and structural subunits. Genes Dev. 1997;11:3109–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Meyerson M, Counter CM, Eaton EN, Ellisen LW, Steiner P, Caddle SD, et al. hEST2, the putative human telomerase catalytic subunit Gene, is up-regulated in tumor cells and during immortalization. Cell. 1997;90:785–95.PubMedView ArticleGoogle Scholar
- Nakamura TM, Morin GB, Chapman KB, Weinrich SL, Andrews WH, Lingner J, et al. Telomerase catalytic subunit homologs from fission yeast and human. Science. 1997;277:955–9.PubMedView ArticleGoogle Scholar
- Lingner J, Hughes TR, Shevchenko A, Mann M, Lundblad V, Cech TR. Reverse transcriptase motifs in the catalytic subunit of telomerase. Science. 1997;276:561–7.PubMedView ArticleGoogle Scholar
- Nugent CI, Lundblad V. The telomerase reverse transcriptase: components and regulation. Genes Dev. 1998;12:1073–85.PubMedView ArticleGoogle Scholar
- Nakayama J, Tahara H, Tahara E, Saito M, Ito K, Nakamura H, et al. Telomerase activation by hTRT in human normal fibroblasts and hepatocellular carcinomas. Nat Genet. 1998;18:65–8.PubMedView ArticleGoogle Scholar
- Rubio MA, Davalos AR, Campisi J. Telomere length mediates the effects of telomerase on the cellular response to genotoxic stress. Exp Cell Res. 2004;298:17–27.PubMedView ArticleGoogle Scholar
- Suram A, Herbig U. The replicometer is broken: telomeres activate cellular senescence in response to genotoxic stresses. Aging Cell. 2014;13:780–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Bischoff C, Petersen HC, Graakjaer J, Andersen-Ranberg K, Vaupel JW, Bohr VA, et al. No association between telomere length and survival among the elderly and oldest old. Epidemiology. 2006;17:190–4.PubMedView ArticleGoogle Scholar
- Kimura M, Cherkas LF, Kato BS, Demissie S, Hjelmborg JB, Brimacombe M, et al. Offspring’s leukocyte telomere length, paternal age, and telomere elongation in sperm. PLoS Genet. 2008;4:e37.PubMedPubMed CentralView ArticleGoogle Scholar
- Steenstrup T. Hjelmborg J v B, Kark JD, Christensen K, Aviv a. The telomere lengthening conundrum--artifact or biology? Nucleic Acids Res. 2013;41:e131.PubMedPubMed CentralView ArticleGoogle Scholar
- Stone RC, Horvath K, Kark JD, Susser E, Tishkoff SA, Aviv A. Telomere length and the cancer–atherosclerosis trade-off. Fisher EMC, editor. PLoS Genet. 2016;12:e1006144.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim EB, Fang X, Fushan AA, Huang Z, Lobanov AV, Han L, et al. Genome sequencing reveals insights into physiology and longevity of the naked mole rat. Nature. 2011;479:223–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Haussmann MF, Winkler DW, Huntington CE, Nisbet ICT, Vleck CM. Telomerase activity is maintained throughout the lifespan of long-lived birds. Exp Gerontol. 2007;42:610–8.PubMedView ArticleGoogle Scholar
- Blasco MA, Lee H-W, Hande MP, Samper E, Lansdorp PM, DePinho RA, et al. Telomere shortening and tumor formation by mouse cells lacking telomerase RNA. Cell. 1997;91:25–34.PubMedView ArticleGoogle Scholar
- Samper E, Flores JM, Blasco MA. Restoration of telomerase activity rescues chromosomal instability and premature aging in Terc−/− mice with short telomeres. EMBO Rep. 2001;2:800–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Blasco MA. Telomeres and human disease: ageing, cancer and beyond. Nat Rev Genet. 2005;6:611–22.PubMedView ArticleGoogle Scholar
- Jaskelioff M, Muller FL, Paik J-H, Thomas E, Jiang S, Adams AC, et al. Telomerase reactivation reverses tissue degeneration in aged telomerase-deficient mice. Nature. 2011;469:102–6.PubMedView ArticleGoogle Scholar
- Charnov EL, Berrigan D. Evolution of life history parameters in animals with indeterminate growth, particularly fish. Evol Ecol. 1991;5:63–8.View ArticleGoogle Scholar
- Klapper W, Heidorn K, Kühne K, Parwaresch R, Krupp G. Telomerase activity in “immortal” fish. FEBS Lett. 1998;434:409–12.PubMedView ArticleGoogle Scholar
- Reznick D, Ghalambor C, Nunney L. The evolution of senescence in fish. Mech Ageing Dev. 2002;123:773–89.PubMedView ArticleGoogle Scholar
- Owen R, Sarkis S, Bodnar A. Developmental pattern of telomerase expression in the sand scallop, Euvola ziczac. Invertebr Biol. 2007;126:40–5.View ArticleGoogle Scholar
- Ebert TA, Russell MP, Gamba G, Bodnar A. Growth, survival, and longevity estimates for the rock-boring sea urchin Echinometra lucunter Lucunter (Echinodermata, Echinoidea) in Bermuda. Bull Mar Sci. 2008;82:381–403.Google Scholar
- Klapper W, Kühne K, Singh KK, Heidorn K, Parwaresch R, Krupp G. Longevity of lobsters is linked to ubiquitous telomerase expression. FEBS Lett. 1998;439:143–6.PubMedView ArticleGoogle Scholar
- Tan TCJ, Rahman R, Jaber-Hijazi F, Felix DA, Chen C, Louis EJ, et al. Telomere maintenance and telomerase activity are differentially regulated in asexual and sexual worms. Proc Natl Acad Sci U S A. 2012;109:4209–14.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith LL, Coller HA, Roberts JM. Telomerase modulates expression of growth-controlling genes and enhances cell proliferation. Nat Cell Biol. 2003;5:474–9.PubMedView ArticleGoogle Scholar
- Parkinson EK, Fitchett C, Cereser B. Dissecting the non-canonical functions of telomerase. Cytogenet Genome Res. 2008;122:273–80.PubMedView ArticleGoogle Scholar
- Choi J, Southworth LK, Sarin KY, Venteicher AS, Ma W, Chang W, et al. TERT promotes epithelial proliferation through transcriptional control of a Myc- and Wnt-related Developmental program. PLoS Genet. 2008;4:e10.PubMedPubMed CentralView ArticleGoogle Scholar
- Low KC, Tergaonkar V. Telomerase: central regulator of all of the hallmarks of cancer. Trends Biochem Sci. 2013;38:426–34.PubMedView ArticleGoogle Scholar
- Ulaner GA, Hu JF, Vu TH, Giudice LC, Hoffman AR. Telomerase activity in human development is regulated by human telomerase reverse transcriptase (hTERT) transcription and by alternate splicing of hTERT transcripts. Cancer Res. 1998;58:4168–72.PubMedGoogle Scholar
- Meier B, Clejan I, Liu Y, Lowden M, Gartner A, Hodgkin J, et al. Trt-1 is the Caenorhabditis elegans catalytic subunit of telomerase. PLoS Genet. 2006;2:e18.PubMedPubMed CentralView ArticleGoogle Scholar
- Sýkorová E, Fajkus J. Structure-function relationships in telomerase genes. Biol Cell. 2009;101:375–92.PubMedView ArticleGoogle Scholar
- Nehyba J, Bose HR Jr, Hrdličková R. The regulation of telomerase by alternative splicing of TERT. INTECH open access Publisher; 2012.Google Scholar
- Withers JB, Ashvetiya T, Beemon KL. Exclusion of Exon 2 is a common mRNA splice variant of primate telomerase reverse Transcriptases. Tian B, editor. PLoS One. 2012;7:e48016.PubMedPubMed CentralView ArticleGoogle Scholar
- Hrdlickova R, Nehyba J, Bose HR. Alternatively spliced telomerase reverse transcriptase variants lacking telomerase activity stimulate cell proliferation. Mol Cell Biol. 2012;32:4283–96.PubMedPubMed CentralView ArticleGoogle Scholar
- Listerman I, Sun J, Gazzaniga FS, Lukas JL, Blackburn EH. The major reverse transcriptase-incompetent splice variant of the human telomerase protein inhibits telomerase activity but protects from apoptosis. Cancer Res. 2013;73:2817–28.PubMedPubMed CentralView ArticleGoogle Scholar
- Wong MS, Chen L, Foster C, Kainthla R, Shay JW, Wright WE. Regulation of telomerase alternative splicing: a target for chemotherapy. Cell Rep. 2013;3:1028–35.PubMedPubMed CentralView ArticleGoogle Scholar
- Wong MS, Wright WE, Shay JW. Alternative splicing regulation of telomerase: a new paradigm? Trends Genet. 2014;30:430–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Riutort M, Field KG, Raff RA, Baguña J. 18S rRNA sequences and phylogeny of platyhelminthes. Biochem Syst Ecol. 1993;21:71–7.View ArticleGoogle Scholar
- Carranza S, Giribet G, Ribera C, Baguñà, Riutort M. Evidence that two types of 18S rDNA coexist in the genome of Dugesia (Schmidtea) mediterranea (Platyhelminthes, Turbellaria, Tricladida). Mol Biol Evol. 1996;13:824–32.PubMedView ArticleGoogle Scholar
- Alvarez-Presas M, Baguñà J, Riutort M. Molecular phylogeny of land and freshwater planarians (Tricladida, Platyhelminthes): from freshwater to land and back. Mol Phylogenet Evol. 2008;47:555–68.PubMedView ArticleGoogle Scholar
- Lang BF, O’Kelly C, Nerad T, Gray MW, Burger G. The closest unicellular relatives of animals. Curr Biol. 2002;12:1773–8.PubMedView ArticleGoogle Scholar
- Ruiz-Trillo I, Roger AJ, Burger G, Gray MW, Lang BF. A Phylogenomic investigation into the origin of Metazoa. Mol Biol Evol. 2008;25:664–72.PubMedView ArticleGoogle Scholar
- Shalchian-Tabrizi K, Minge MA, Espelund M, Orr R, Ruden T, Jakobsen KS, et al. Multigene phylogeny of Choanozoa and the origin of animals. Aramayo R, editor. PLoS One. 2008;3:e2098.PubMedPubMed CentralView ArticleGoogle Scholar
- Torruella G, Derelle R, Paps J, Lang BF, Roger AJ, Shalchian-Tabrizi K, et al. Phylogenetic relationships within the Opisthokonta based on Phylogenomic analyses of conserved single-copy protein domains. Mol Biol Evol. 2012;29:531–44.PubMedView ArticleGoogle Scholar
- Suga H, Chen Z, de Mendoza A, Sebé-Pedrós A, Brown MW, Kramer E, et al. The Capsaspora genome reveals a complex unicellular prehistory of animals. Nat Commun. 2013;4:2325.PubMedPubMed CentralView ArticleGoogle Scholar
- Robertson HM. The choanoflagellate Monosiga Brevicollis karyotype revealed by the genome sequence: telomere-linked helicase genes resemble those of some fungi. Chromosom Res. 2009;17:873–82.View ArticleGoogle Scholar
- Fairclough SR, Chen Z, Kramer E, Zeng Q, Young S, Robertson HM, et al. Premetazoan genome evolution and the regulation of cell differentiation in the choanoflagellate Salpingoeca Rosetta. Genome Biol. 2013;14:R15.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu P, Widmer G, Wang Y, Ozaki LS, Alves JM, Serrano MG, et al. The genome of Cryptosporidium hominis. Nature. 2004;431:1107–12.PubMedView ArticleGoogle Scholar
- Figueiredo LM, Rocha EPC, Mancio-Silva L, Prevost C, Hernandez-Verdun D, Scherf A. The unusually large plasmodium telomerase reverse-transcriptase localizes in a discrete compartment associated with the nucleolus. Nucleic Acids Res. 2005;33:1111–22.PubMedPubMed CentralView ArticleGoogle Scholar
- Hall N, Karras M, Raine JD, Carlton JM, Kooij TWA, Berriman M, et al. A comprehensive survey of the plasmodium life cycle by genomic, transcriptomic, and proteomic analyses. Science. 2005;307:82–6.PubMedView ArticleGoogle Scholar
- Pain A, Renauld H, Berriman M, Murphy L, Yeats CA, Weir W, et al. Genome of the host-cell transforming parasite Theileria annulata compared with T. parva. Science. 2005;309:131–3.PubMedView ArticleGoogle Scholar
- Dreesen O, Li B, Cross GAM. Telomere structure and shortening in telomerase-deficient Trypanosoma brucei. Nucleic Acids Res. 2005;33:4536–43.PubMedPubMed CentralView ArticleGoogle Scholar
- El-Sayed NM, Myler PJ, Blandin G, Berriman M, Crabtree J, Aggarwal G, et al. Comparative genomics of trypanosomatid parasitic protozoa. Science. 2005;309:404–9.PubMedView ArticleGoogle Scholar
- Giardini MA, Lira CBB, Conte FF, Camillo LR, de Siqueira Neto JL, Ramos CHI, et al. The putative telomerase reverse transcriptase component of Leishmania amazonensis: gene cloning and characterization. Parasitol Res. 2006;98:447–54.PubMedView ArticleGoogle Scholar
- Morrison HG, McArthur AG, Gillin FD, Aley SB, Adam RD, Olsen GJ, et al. Genomic minimalism in the early diverging intestinal parasite Giardia lamblia. Science. 2007;317:1921–6.PubMedView ArticleGoogle Scholar
- Chapman JA, Kirkness EF, Simakov O, Hampson SE, Mitros T, Weinmaier T, et al. The dynamic genome of Hydra. Nature. 2010;464:592–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Steele RE, David CN, Technau U. A genomic view of 500 million years of cnidarian evolution. Trends Genet. 2011;27:7–13.PubMedView ArticleGoogle Scholar
- Srivastava M, Simakov O, Chapman J, Fahey B, Gauthier MEA, Mitros T, et al. The Amphimedon Queenslandica genome and the evolution of animal complexity. Nature. 2010;466:720–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Ryan JF, Pang K, Schnitzler CE, Nguyen A-D, Moreland RT, Simmons DK, et al. The genome of the ctenophore Mnemiopsis leidyi and its implications for cell type evolution. Science. 2013;342:1242592.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Y, Xiong J, Zhou Z, Huo F, Miao W, Ran C, et al. The genome of the Myxosporean Thelohanellus kitauei shows adaptations to nutrient acquisition within its fish host. Genome Biol Evol. 2014;6:3182–98.PubMedPubMed CentralView ArticleGoogle Scholar
- Baumgarten S, Simakov O, Esherick LY, Liew YJ, Lehnert EM, Michell CT, et al. The genome of Aiptasia, a sea anemone model for coral symbiosis. Proc Natl Acad Sci. 2015;112:11893–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Bryan TM, Sperger JM, Chapman KB, Cech TR. Telomerase reverse transcriptase genes identified in Tetrahymena Thermophila and Oxytricha trifallax. Proc Natl Acad Sci. 1998;95:8479–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Collins K, Gandhi L. The reverse transcriptase component of the Tetrahymena telomerase ribonucleoprotein complex. Proc Natl Acad Sci. 1998;95:8485–90.PubMedPubMed CentralView ArticleGoogle Scholar
- Greenberg RA, Allsopp RC, Chin L, Morin GB, DePinho RA. Expression of mouse telomerase reverse transcriptase during development, differentiation and proliferation. Oncogene. 1998;16:1723–30.PubMedView ArticleGoogle Scholar
- Martin-Rivera L, Herrera E, Albar JP, Blasco MA. Expression of mouse telomerase catalytic subunit in embryos and adult tissues. Proc Natl Acad Sci. 1998;95:10471–6.PubMedPubMed CentralView ArticleGoogle Scholar
- O’Reilly M, Teichmann SA, Rhodes D. Telomerases. Curr Opin Struct Biol. 1999;9:56–65.PubMedView ArticleGoogle Scholar
- Lau BW-M, Wong AO-L, Tsao GS-W, So K-F, Yip HK-F. Molecular cloning and characterization of the zebrafish (Danio rerio) telomerase catalytic subunit (telomerase reverse transcriptase, TERT). J Mol Neurosci. 2008;34:63–75.PubMedView ArticleGoogle Scholar
- Li Y, Yates JA, Chen JJ-L. Identification and characterization of sea squirt telomerase reverse transcriptase. Gene. 2007;400:16–24.PubMedView ArticleGoogle Scholar
- Cai Y, Ai Y, Zhao Q, Li J, Yang G, Gong P, et al. Cloning and characterization of telomerase reverse transcriptase gene in Trichinella spiralis. Parasitol Res. 2012;110:411–7.PubMedView ArticleGoogle Scholar
- Wurm Y, Wang J, Riba-Grognuz O, Corona M, Nygaard S, Hunt BG, et al. The genome of the fire ant Solenopsis invicta. Proc Natl Acad Sci U S A. 2011;108:5679–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Hrdličková R, Nehyba J, Lim SL, Grützner F, Bose HR. Insights into the evolution of mammalian telomerase: platypus TERT shares similarities with genes of birds and other reptiles and localizes on sex chromosomes. BMC Genomics. 2012;13:216.PubMedPubMed CentralView ArticleGoogle Scholar
- Schumpert C, Nelson J, Kim E, Dudycha JL, Patel RC. Telomerase activity and telomere length in Daphnia. Saretzki G, editor. PLoS One. 2015;10:e0127196.PubMedPubMed CentralView ArticleGoogle Scholar
- Koutsovoulos G, Kumar S, Laetsch DR, Stevens L, Daub J, Conlon C, et al. No evidence for extensive horizontal gene transfer in the genome of the tardigrade Hypsibius dujardini. Proc Natl Acad Sci. 2016;113:5053–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Hashimoto T, Horikawa DD, Saito Y, Kuwahara H, Kozuka-Hata H, Shin-I T, et al. Extremotolerant tardigrade genome and improved radiotolerance of human cultured cells by tardigrade-unique protein. Nat Commun. 2016;7:12808.PubMedPubMed CentralView ArticleGoogle Scholar
- Kao D, Lai AG, Stamataki E, Rosic S, Konstantinides N, Jarvis E, et al. The genome of the crustacean Parhyale hawaiensis, a model for animal development, regeneration, immunity and lignocellulose digestion. elife. 2016;5:1.View ArticleGoogle Scholar
- Mason JM, Randall TA, Capkova Frydrychova R. Telomerase lost? Chromosoma. 2016;125:65–73.PubMedView ArticleGoogle Scholar
- Mark Welch D, Meselson M. Evidence for the evolution of bdelloid rotifers without sexual reproduction or genetic exchange. Science. 2000;288:1211–5.PubMedView ArticleGoogle Scholar
- Gladyshev EA, Arkhipova IR. Genome structure of Bdelloid rotifers: shaped by asexuality or desiccation? J Hered. 2010;101:S85–93.PubMedView ArticleGoogle Scholar
- Pardue M-L, DeBaryshe PG. Retrotransposons provide an evolutionarily robust non-telomerase mechanism to maintain telomeres. Annu Rev Genet. 2003;37:485–511.PubMedView ArticleGoogle Scholar
- Biessmann H, Carter SB, Mason JM. Chromosome ends in Drosophila without telomeric DNA sequences. Proc Natl Acad Sci U S A. 1990;87:1758–61.PubMedPubMed CentralView ArticleGoogle Scholar
- Biessmann H, Champion LE, O’Hair M, Ikenaga K, Kasravi B, Mason JM. Frequent transpositions of Drosophila melanogaster HeT-A transposable elements to receding chromosome ends. EMBO J. 1992;11:4459–69.PubMedPubMed CentralGoogle Scholar
- Casacuberta E, Pardue M-L. Coevolution of the telomeric retrotransposons across Drosophila species. Genetics. 2002;161:1113–24.PubMedPubMed CentralGoogle Scholar
- Casacuberta E, Pardue M-L. Transposon telomeres are widely distributed in the Drosophila genus: TART elements in the virilis group. Proc Natl Acad Sci. 2003;100:3363–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Abad JP, De Pablos B, Osoegawa K, De Jong PJ, Martín-Gallardo A, Villasante A. TAHRE, a novel telomeric retrotransposon from Drosophila melanogaster, reveals the origin of Drosophila telomeres. Mol Biol Evol. 2004;21:1620–4.PubMedView ArticleGoogle Scholar
- Gladyshev EA, Arkhipova IR. A subtelomeric non-LTR retrotransposon Hebe in the bdelloid rotifer Adineta vaga is subject to inactivation by deletions but not 5′ truncations. Mob DNA. 2010;1:12.PubMedPubMed CentralView ArticleGoogle Scholar
- Gladyshev EA, Arkhipova IR. A widespread class of reverse transcriptase-related cellular genes. Proc Natl Acad Sci U S A. 2011;108:20311–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Gladyshev EA, Arkhipova IR. Telomere-associated endonuclease-deficient Penelope-like retroelements in diverse eukaryotes. Proc Natl Acad Sci. 2007;104:9352–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Kelleher C, Teixeira MT, Förstemann K, Lingner J. Telomerase: biochemical considerations for enzyme and substrate. Trends Biochem Sci. 2002;27:572–9.PubMedView ArticleGoogle Scholar
- Malik HS, Burke WD, Eickbush TH. Putative telomerase catalytic subunits from Giardia lamblia and Caenorhabditis elegans. Gene. 2000;251:101–8.PubMedView ArticleGoogle Scholar
- Adam RD, Dahlstrom EW, Martens CA, Bruno DP, Barbian KD, Ricklefs SM, et al. Genome sequencing of Giardia lamblia genotypes A2 and B isolates (DH and GS) and comparative analysis with the genomes of genotypes A1 and E (WB and pig). Genome Biol Evol. 2013;5:2498–511.PubMedPubMed CentralView ArticleGoogle Scholar
- Stevens JR, Noyes HA, Schofield CJ, Gibson W. The molecular evolution of Trypanosomatidae. Adv Parasitol. 2001;48:1–53.PubMedView ArticleGoogle Scholar
- Armbruster BN, Banik SS, Guo C, Smith AC, Counter CM. N-terminal domains of the human telomerase catalytic subunit required for enzyme activity in vivo. Mol Cell Biol. 2001;21:7775–86.PubMedPubMed CentralView ArticleGoogle Scholar
- Moriarty TJ, Marie-Egyptienne DT, Autexier C. Functional Organization of Repeat Addition Processivity and DNA synthesis determinants in the human telomerase Multimer. Mol Cell Biol. 2004;24:3720–33.PubMedPubMed CentralView ArticleGoogle Scholar
- Sekaran VG, Soares J, Jarstfer MB. Structures of telomerase subunits provide functional insights. Biochim Biophys Acta-Proteins Proteomics. 1804;2010:1190–201.Google Scholar
- Wyatt HDM, West SC, Beattie TL. InTERTpreting telomerase structure and function. Nucleic Acids Res. 2010;38:5609–22.PubMedPubMed CentralView ArticleGoogle Scholar
- Gillis AJ, Schuller AP, Skordalakes E. Structure of the Tribolium castaneum telomerase catalytic subunit TERT. Nature. 2008;455:633–7.PubMedView ArticleGoogle Scholar
- Rouda S, Skordalakes E. Structure of the RNA-binding domain of telomerase: implications for RNA recognition and binding. Structure. 2007;15:1403–12.PubMedView ArticleGoogle Scholar
- Mitchell M, Gillis A, Futahashi M, Fujiwara H, Skordalakes E. Structural basis for telomerase catalytic subunit TERT binding to RNA template and telomeric DNA. Nat Struct Mol Biol. 2010;17:513–8.PubMedView ArticleGoogle Scholar
- Gu X. Statistical methods for testing functional divergence after gene duplication. Mol Biol Evol. 1999;16:1664–74.PubMedView ArticleGoogle Scholar
- Huang H, Chopra R, Verdine GL, Harrison SC. Structure of a covalently trapped catalytic complex of HIV-1 reverse transcriptase: implications for drug resistance. Science. 1998;282:1669–75.PubMedView ArticleGoogle Scholar
- Hossain S, Singh S, Lue NF. Functional analysis of the C-terminal extension of telomerase reverse transcriptase. A putative “thumb” domain. J Biol Chem. 2002;277:36174–80.PubMedView ArticleGoogle Scholar
- Huard S, Moriarty TJ, Autexier C. The C terminus of the human telomerase reverse transcriptase is a determinant of enzyme processivity. Nucleic Acids Res. 2003;31:4059–70.PubMedPubMed CentralView ArticleGoogle Scholar
- Tomlinson CG, Holien JK, Mathias JAT, Parker MW, Bryan TM. The C-terminal extension of human telomerase reverse transcriptase is necessary for high affinity binding to telomeric DNA. Biochimie. 2016;128–129:114–21.PubMedView ArticleGoogle Scholar
- Morand S. Life-history traits in parasitic nematodes: a comparative approach for the search of invariants. Funct Ecol. 1996;10:210–8.View ArticleGoogle Scholar
- Friedman KL, Cech TR. Essential functions of amino-terminal domains in the yeast telomerase catalytic subunit revealed by selection for viable mutants. Genes Dev. 1999;13:2863–74.PubMedPubMed CentralView ArticleGoogle Scholar
- Beattie TL, Zhou W, Robinson MO, Harrington L. Functional multimerization of the human telomerase reverse transcriptase. Mol Cell Biol. 2001;21:6151–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Lai CK, Mitchell JR, Collins K. RNA binding domain of telomerase reverse transcriptase. Mol Cell Biol. 2001;21:990–1000.PubMedPubMed CentralView ArticleGoogle Scholar
- Moriarty TJ, Huard S, Dupuis S, Autexier C. Functional multimerization of human telomerase requires an RNA interaction domain in the N terminus of the catalytic subunit. Mol Cell Biol. 2002;22:1253–65.PubMedPubMed CentralView ArticleGoogle Scholar
- Bley CJ, Qi X, Rand DP, Borges CR, Nelson RW, Chen JJ-L. RNA-protein binding interface in the telomerase ribonucleoprotein. Proc Natl Acad Sci U S A. 2011;108:20333–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Harkisheimer M, Mason M, Shuvaeva E, Skordalakes E. A motif in the vertebrate telomerase N-terminal linker of TERT contributes to RNA binding and telomerase activity and processivity. Structure. 2013;21:1870–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Autexier C, Greider CW. Boundary elements of the Tetrahymena telomerase RNA template and alignment domains. Genes Dev. 1995;9:2227–39.PubMedView ArticleGoogle Scholar
- Tzfati Y, Fulton TB, Roy J, Blackburn EH. Template boundary in a yeast telomerase specified by RNA structure. Science. 2000;288:863–7.PubMedView ArticleGoogle Scholar
- Lai CK, Miller MC, Collins K. Template boundary definition in Tetrahymena telomerase. Genes Dev. 2002;16:415–20.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen J-L, Greider CW. Template boundary definition in mammalian telomerase. Genes Dev. 2003;17:2747–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Kumar S, Hedges SB. A molecular timescale for vertebrate evolution. Nature. 1998;392:917–20.PubMedView ArticleGoogle Scholar
- Hedges SB. The origin and evolution of model organisms. Nat Rev Genet. 2002;3:838–49.PubMedView ArticleGoogle Scholar
- Raible F, Tessmar-Raible K, Osoegawa K, Wincker P, Jubin C, Balavoine G, et al. Vertebrate-type intron-rich genes in the marine annelid Platynereis dumerilii. Science. 2005;310:1325–6.PubMedView ArticleGoogle Scholar
- Colbourne JK, Pfrender ME, Gilbert D, Thomas WK, Tucker A, Oakley TH, et al. The ecoresponsive genome of Daphnia pulex. Science. 2011;331:555–61.PubMedPubMed CentralView ArticleGoogle Scholar
- Saebøe-Larssen S, Fossberg E, Gaudernack G. Characterization of novel alternative splicing sites in human telomerase reverse transcriptase (hTERT): analysis of expression and mutual correlation in mRNA isoforms from normal and tumour tissues. BMC Mol Biol. 2006;7:26.PubMedPubMed CentralView ArticleGoogle Scholar
- Chang H, Delany ME. Complicated RNA splicing of chicken telomerase reverse transcriptase revealed by profiling cells both positive and negative for telomerase activity. Gene. 2006;379:33–9.PubMedView ArticleGoogle Scholar
- Imamura S, Uchiyama J, Koshimizu E, Hanai J, Raftopoulou C, Murphey RD, et al. A non-canonical function of Zebrafish telomerase reverse transcriptase is required for Developmental Hematopoiesis. PLoS One. 2008;3:e3364.PubMedPubMed CentralView ArticleGoogle Scholar
- Ulaner GA, Hu J-F, Vu TH, Oruganti H, Giudice LC, Hoffman AR. Regulation of telomerase by alternate splicing of human telomerase reverse transcriptase (hTERT) in normal and neoplastic ovary, endometrium and myometrium. Int J Cancer. 2000;85:330–5.PubMedView ArticleGoogle Scholar
- Tasaka K, Yokoyama N, Nodono H, Hoshi M, Matsumoto M. Innate sexuality determines the mechanisms of telomere maintenance. Int J Dev Biol. 2013;57:69–72.PubMedView ArticleGoogle Scholar
- Robb SMC, Ross E, Sánchez Alvarado A. SmedGD: the Schmidtea mediterranea genome database. Nucleic Acids Res. 2008;36:D599–606.PubMedView ArticleGoogle Scholar
- Nishimura O, Hirao Y, Tarui H, Agata K. Comparative transcriptome analysis between planarian Dugesia Japonica and other platyhelminth species. BMC Genomics. 2012;13:289.PubMedPubMed CentralView ArticleGoogle Scholar
- Wheeler NJ, Agbedanu PN, Kimber MJ, Ribeiro P, Day TA, Zamanian M. Functional analysis of Girardia tigrina transcriptome seeds pipeline for anthelmintic target discovery. Parasit Vectors. 2015;8:34.PubMedPubMed CentralView ArticleGoogle Scholar
- Robb SMC, Gotting K, Ross E, Sánchez AA. SmedGD 2.0: the Schmidtea mediterranea genome database. Genesis. 2015;53:535–46.PubMedPubMed CentralView ArticleGoogle Scholar
- Egger B, Lapraz F, Tomiczek B, Müller S, Dessimoz C, Girstmair J, et al. A Transcriptomic-Phylogenomic analysis of the evolutionary relationships of flatworms. Curr Biol. 2015;25:1347–53.PubMedPubMed CentralView ArticleGoogle Scholar
- Brandl H, Moon H, Vila-Farré M, Liu S-Y, Henry I, Rink JC. PlanMine--a mineable resource of planarian biology and biodiversity. Nucleic Acids Res. 2016;44:D764–73.PubMedView ArticleGoogle Scholar
- Young ND, Nagarajan N, Lin SJ, Korhonen PK, Jex AR, Hall RS, et al. The Opisthorchis viverrini genome provides insights into life in the bile duct. Nat Commun. 2014;5:4378.PubMedPubMed CentralGoogle Scholar
- Howe KL, Bolt BJ, Cain S, Chan J, Chen WJ, Davis P, et al. WormBase 2016: expanding to enable helminth genomic research. Nucleic Acids Res. 2016;44:D774–80.PubMedView ArticleGoogle Scholar
- Tsai IJ, Zarowiecki M, Holroyd N, Garciarrubio A, Sanchez-Flores A, Brooks KL, et al. The genomes of four tapeworm species reveal adaptations to parasitism. Nature. 2013;496:57–63.PubMedPubMed CentralView ArticleGoogle Scholar
- Wasik K, Gurtowski J, Zhou X, Ramos OM, Delás MJ, Battistoni G, et al. Genome and transcriptome of the regeneration-competent flatworm, Macrostomum lignano. Proc Natl Acad Sci. 2015;112:12462–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu S-Y, Selck C, Friedrich B, Lutz R, Vila-Farré M, Dahl A, et al. Reactivating head regrowth in a regeneration-deficient planarian species. Nature. 2013;500:81–4.PubMedView ArticleGoogle Scholar
- Lue NF. Adding to the ends: what makes telomerase processive and how important is it? BioEssays. 2004;26:955–62.PubMedView ArticleGoogle Scholar
- Collins K. Ciliate telomerase biochemistry. Annu Rev Biochem. 1999;68:187–218.PubMedView ArticleGoogle Scholar
- Peng Y, Mian IS, Lue NF. Analysis of telomerase processivity: mechanistic similarity to HIV-1 reverse transcriptase and role in telomere maintenance. Mol Cell. 2001;7:1201–11.PubMedView ArticleGoogle Scholar
- Aken BL, Ayling S, Barrell D, Clarke L, Curwen V, Fairley S, et al. The Ensembl Gene annotation system. Database (Oxford). 2016;2016:baw093.View ArticleGoogle Scholar
- Katoh K, Asimenos G, Toh H. Multiple alignment of DNA sequences with MAFFT. Methods Mol Biol. 2009;537:39–64.PubMedView ArticleGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.PubMedPubMed CentralView ArticleGoogle Scholar
- González-Estévez C, Felix DA, Smith MD, Paps J, Morley SJ, James V, et al. SMG-1 and mTORC1 act antagonistically to regulate response to injury and growth in planarians. PLoS Genet. 2012;8:e1002619.PubMedPubMed CentralView ArticleGoogle Scholar
- Lek M, MacArthur DG, Yang N, North KN. Phylogenetic analysis of gene structure and alternative splicing in alpha-actinins. Mol Biol Evol. 2010;27:773–80.PubMedView ArticleGoogle Scholar