- Research article
- Open Access
Polymorphisms within the Telomerase Reverse Transcriptase gene (TERT) in four breeds of dogs selected for difference in lifespan and cancer susceptibility
BMC Veterinary Research volume 10, Article number: 20 (2014)
Enzymatic activity of Telomerase Reverse Transcriptase (TERT) is important in maintaining the telomere length and has been implicated in cancer and aging related pathology. Since cancer susceptibility as well as longevity of dogs vary between breeds, this study involved sequencing the entire TERT gene of Canis familiaris from DNA samples obtained from forty dogs, with ten dogs each of four breeds: Shih Tzu, Dachshund, Irish Wolfhound, and Newfoundland, each with different life expectancies and susceptibility to cancer.
We compared the sequences of all forty individuals amongst one another and with the published sequence of canine TERT, and analyzed relationships between members of the same or different breeds. Two separate phylogenetic trees were generated and analyzed from these individuals. Polymorphisms were found most frequently in intronic regions of the gene, although exonic polymorphisms also were observed. In many locations genotypes were observed that were either homozygous for the reference sequence or heterozygous, but the variant homozygous genotype was not observed.
We propose that these homozygous variants are likely to have adverse effects in dogs. It was also found that the polymorphisms did not segregate by breed. Because the four breeds chosen come from geographically and physiologically distinct backgrounds, it can be inferred that the polymorphic diversification of TERT preceded breed derivation.
Telomeres are linked to susceptibility to cancer, the aging process, and the lifespan of individual cells . Also, the dysfunction of a single enzyme involved in telomere maintenance has been shown to be sufficient to cause severe telomere damage, premature tissue aging, and development of neoplastic lesions . Telomere length is controlled by the enzyme telomerase, a component of which is Telomerase Reverse Transcriptase (TERT). Therefore, variations in telomerase activity could result from coding polymorphisms affecting the enzymatic function of TERT, or from non-coding polymorphisms affecting the magnitude or kinetics of TERT expression. Interestingly, it has also been observed that while deficient telomerase activity results in a decreased rate of cancer in mice , over-expression and ubiquitious expression of telomerase result in an increased rate of tumors .
Domestic dogs provide an exceptional experimental model to assess the potential effect of such polymorphisms on organismal lifespan, as it is well recognized that many large breed dogs age faster and have shorter life spans than dogs of small breeds. Furthermore, cancer susceptibility varies widely by breed  making dogs excellent models for studying cancer susceptibility. It has been specifically shown that in dogs, as in humans, cancer cells rely on telomerase to maintain telomeres and extend their lifespan .
For this study, polymorphisms in the TERT sequence were observed in Shih Tzus, Dachshunds, Irish Wolfhounds and Newfoundlands: four breeds that show large differences in size, expected lifespan, cancer susceptibility, and geographic origins. Shih Tzus are small dogs that have an average lifespan of 13 years ; the breed is thought to be “ancient”, and it dates back at least to the year 624 A.D. when it was first depicted in Chinese art . These dogs are thought to be at a decreased risk for cancer on the whole . Dachshunds are small dogs that were first bred in Germany in the early 1600s  and also have an average lifespan of 13 years . Evidence points to them having a decreased risk of cancer compared to the average, but they appear to have a predilection for mammary tumors . Records of Irish Wolfhounds, large dogs first bred in Ireland, date back to 391 A.D. . These dogs only live 7 years on average , and are particularly susceptible to osteosarcoma, even as compared to other large, fast-growing breeds that are typically prone to this cancer . Finally, the large Newfoundland breed was developed in Newfoundland, Canda in the early 1600s  and live 9 years on average . These dogs are thought to be no more or less susceptible to cancer than average .
Although reproductive isolation among breeds was not formally standardized until the mid-19th century with the formation of breed clubs and breed standards , geographic isolation of these breeds in their place of origin makes interbreeding among these breeds unlikely. Furthermore, Parker et al. showed that modern dog breeds are “distinct genetic units” and assigned the dog breeds used in this study to genetically distinct “clusters”.
Our objectives here were to assess conserved polymorphism in TERT, and to establish if they segregated according to the breeds or randomly. In order to do so, we sequenced the TERT gene of 10 dogs each from the aforementioned breeds and compared their sequences and any polymorphism that were found.
Dr. Heidi Parker (NHGRI, NIH, Bethesda, MD) provided de-identified and breed-verified samples from forty dogs including ten Shi Tzu, ten Dachshund, ten Newfoundland, and ten Irish Wolfhounds. All ethical approvals for protocols, including IACUC approval for collection of the samples, were obtained for this study.
We used the genomic sequence of canine TERT (TERT_CANFA ENSCAFT00000017081 CanFam 2.0) [Release 58 May 2010 ] as a reference, including 500 additional nucleotides flanking either end of the gene. The plasmid editing program, ApE (Wayne Davis; Salt Lake City, Utah; version “ApE_OSX_1_1_7.dmg”)  was used to parcel exons, introns, and previously catalogued polymorphisms. As depicted in Additional file 1: Figure S1, the gene was divided into 500 nucleotide-long segments, with each segment overlapping by 250 nucleotides, for a total of 74 sections. The second primer pair could not be generated due to extensive GC pairing. ApE was used to generate forward and reverse primers for each section. We then manually reviewed each list of potential primers per section, with the final primer pairs chosen based on stability, GC content, melting point, and salt content. The primer sequences are listed in Additional file 1: Table S1.
We extracted DNA from lymph node tissue of a healthy dog obtained with IACUC approval, and used it at a concentration of 40 ng/μL to standardize the conditions of PCR amplification. The performance of each pair was verified by conventional PCR and gel electrophoresis. Of 69 primer pairs that successfully amplified DNA, we chose 11 that were interspersed throughout the gene sequence. The amplified DNA was then sequenced in both forward and reverse directions for a total of 22 samples. These partial DNA sequences were examined for length, quality, and how closely they matched the reference canine TERT sequence.
Two sets of twenty dogs were used to derive the sequence: the first set was used to run primer pairs 1 through 40 and 66 through 75, and the second set was used to run primer pairs 41 through 65. A DNA concentration of 20 ng/μL was used for each sample for PCR, with go-Taq polymerase (Promega). PCR was performed at 42-55°C for 40 cycles. Each PCR was verified by carrying out gel electrophoresis of 8 μL of sample in a 2% agarose gel. The remaining 12 μL of successful samples were subjected to capillary sequencing (Biomedical Genomic Center, University of Minnesota).
Sequence analysis and alignment
We compiled the resulting sequences using CodonCode Aligner (CodonCode; Dedham, Massachusetts)  with default parameters (Version 22.214.171.124), aligning each amplicon to the published reference sequence. The consensus sequences of each dog were aligned using clustalx (Version 2.0.11) . We considered three or more identical variants from the original reference sequence to be polymorphisms. At this point each group of 20 dogs was handled separately. Gaps and ambiguous alignment segments were trimmed in Jalview4 (Version 11.0)  prior to input into three complementary programs for phylogenetic tree determination: 1. MrBayes (Version 3.2) , a Bayesian inference algorithm that utilizes Markov Chain Monte Carlo optimization (parameters: default General Time Reversible model of evolution, invariant sites, gamma distributed rates, 2 million generations) 2. Simple DNA parsimony as implemented in dnapars using consense (Majority Rule consensus) from the Phylip package (Version 3.69) . Neighbor Joining algorithm implemented in clustalx. We visualized trees output by the tree-building programs using FigTree .
The DNA samples from individual dogs were prepared and PCR carried out as described in Supplementary Materials. PCR reactions were originally carried out for 74 primer pairs, and it was found that 69 of these worked on the sample dog (Additional file 1: Figure S1). It was verified that 1,245 of the 1,380 PCR samples yielded a clearly defined band of the expected size (Additional file 1: Table S1, Additional file 1: Figure S1).
Clustalx was used to analyze the consensus sequence of all twenty dogs and compared them with the published sequence. As hypothesized, we identified polymorphisms within both exons and introns (Figure 1a and b). Exonic sequences within the TERT gene were more conserved compared to the intronic sequences. Strikingly, for many observed SNPs we found two groups: one predominantly homozyougous for the reference sequence and the other heterozygous. However, there were specific loci where the variant homozygous allele was never seen when a heterozygous allele was present. For example, at bp 19298 of the TERT gene (Figure 1a), located within exon 17, dogs of all breeds surveyed here were either homozygous C or heterozygous C/G (S), but none were homozygous G. It is unlikely this was the result of read contamination, as the sequence trace chromatograms for this section were clean and the peaks were clear (Figure 1b). A further example of this variation pattern can be seen at bp 17273 of the TERT gene (Figure 1c and d), located within an intronic region: dogs of all breeds surveyed here were either homozygous A or heterozygous A/C (M), but none were homozygous C. This variation pattern was observed throughout the gene in at least 115 instances in both sets of 20 samples, with 14 of these instances occurring in exonic regions, and 7 of those exonic variants causing a coding change. The locations of the exonic changes and their effect on codons within the gene are displayed in Table 1. A complete analysis of our sequence and the SNPs obtained is shown in Additional file 1: Figure S2 and S3.
In analyzing the phylogenetic tree obtained from the two separate sets of dogs at the TERT locus, it was found there was no instance that placed all individuals of the same breed in one branch. Figures 2a and b illustrate this: even when as many as three members of the same breed were grouped together there was at least one other dog from a different breed in that branch. This was not an artifact of the phylogenetic inference method used, as nearly all clades within the phylogeny having high confidence (i.e., with posterior probabilities > 50%) were preserved by diverse phylogenetic methods (MrBayes, parsimony and Neighbor-joining) (data not shown).
We found that polymorphisms exist within both exons and introns of the TERT gene. A salient and unexpected observation was that in some locations, individuals were either homozygous to the reference sequence, or had polymorphisms that were heterozygous SNPs. The variant homozygote SNPs at these loci were never observed. This is likely due to the fact that dog breeds as populations are not subject to Hardy-Weinberg Equilibrium: matings are not random, and traits are deliberately selected for or against by human breeders. However, this phenomenon may also be explained by a deleterious effect of the alternative homozygous variant: it could be lethal for a conceptus. This is a definite possibility, given that inbred populations tend to have a higher concentration of deleterious recessive traits than less inbred populations . Alternatively, the small sample size analyzed could be the reason for the lack of variants. In 7 of the 14 exonic occurrences, the variant caused a coding change and thus could prove deleterious. However, it is unknown if the alleles have dominant negative effects.
We expected to observe breed-specific SNPs shared amongst the majority of individuals in each breed. We also hypothesized the specific SNPs shared in each breed would be distinct from other breeds, possibly in correlation to breed size and lifespan. However, it was found that SNPs occurred in individuals without breed specificity, and there was no correlation between dog size and the observed polymorphisms; most strikingly this was seen in the close relation of the Shi Tzu and Irish Wolfhound in Figure 2b (SH2 and IW5, respectively). This lack of distinct breed association in addition to the geographic and genetic isolation of the breeds’ places of origin implies that polymorphic diversification of TERT preceded derivation of the four breeds studied. Despite the founder effects, population bottlenecks, and careful human selection for or against traits involved in breed creation, the TERT gene has remained diverse within and among dog breeds. The gene’s mutation rate and the prevalence of mutations do not seem to have been drastically altered by breed creation.
Access to breeds was limited in this study. Thus, future studies with expanded the sample sizes (individuals and breeds) and where the TERT genes in wolves and indigenous dogs from South East Asia were sequenced would further enhance our understanding of how and if the TERT gene has changed with the derivation of modern dog breeds. This also could help answer whether TERT is resistant to high levels of inbreeding, clarifying the actual distribution of the homozygous variants and their potential effects on viability.
Our results suggest that the TERT gene is not wholly responsible for differences in lifespan and cancer susceptibility seen in dog breeds. Furthermore, the diversity and surprisingly higher number of SNPs in TERT makes us infer that mutations in this sequence are well tolerated. However, given the heterozygous SNPs that were found, we believe that the gene is not immune to the depressive effects of inbreeding. Furthermore, previous work has identified quantitative trait loci in canine chromosome 7 (CFA7) and 15 (CFA15) that were associated with size and longevity . The IGF1 gene in CFA15 was proposed as one candidate having a major influence on these traits. However, the density of SNPs used in this study would not conclusively eliminate TERT as an additional, related or independent modulator.
Mutations in the human TERT gene have been shown to be important in dyskeratosis congenita, aplastic anemia, and idiopathic pulmonary fibrosis, among other disorders . In addition, heterozygous mutations in the TERT gene have been shown to impair telomerase activity through haploinsufficiency . Despite the fact that many of these disorders are heritable, it is not clear as to what extent mutations of TERT contribute to these diseases. Given our findings, we propose that these diseases do not arise solely due to polymorphic variants of the TERT gene, especially since the canine TERT protein is the closest homologue to the human protein known thus far .
Availability of supporting data
The data set(s) supporting the results of this article are included within the article (and its additional files).
Calado RT, Young NS: Telomere diseases. N Engl J Med. 2009, 361: 2353-2365. 10.1056/NEJMra0903373.
Donate LE, Blasco MA: Telomeres in cancer and ageing. Phil Trans R Soc B. 2011, 366: 76-84. 10.1098/rstb.2010.0291.
González-Suárez E, Samper E, Flores JM, Blasco MA: Telomerase-deficient mice with short telomeres are resistant to skin tumorigenesis. Nat Genet. 2000, 26: 114-117. 10.1038/79089.
González-Suárez E, Geserick C, Flores JM, Blasco MA: Antagonist effects of telomerase on cancer and aging in K5-mTert transgenic mice. Oncogene. 2005, 24: 2256-2270. 10.1038/sj.onc.1208413.
Dobson JM: Breed-predispositions to cancer in pedigree dogs. ISRN Vet Sci. 2013, 2013: 941275.
Nasir L: Telomeres and telomerase: biological and clinical importance in dogs. Vet J. 2008, 175 (2): 155-163. 10.1016/j.tvjl.2007.01.024.
Cassidy KM: Dog Longevity. [http://users.pullman.com/lostriver/breeddata.htm]
American kennel club. [http://www.akc.org/index.cfm?nav_area=homepage]
MacVean DW, Monlux AW, Anderson PS, Silberg SL, Roszel JF: Frequency of canine and feline tumors in a defined population. Vet Pathol. 1978, 15: 700-715. 10.1177/030098587801500602.
Urfer SR, Gaillard C, Steiger A: Lifespan and disease predispositions in the Irish wolfhound: a review. Vet Q. 2007, 29 (Suppl 3): 102-111.
Parker HG, Kim LV, Sutter NB, Carlson S, Lorentzen TD, Malek TB, Johnson GS, DeFrance HB, Ostrander EA, Kruglyak L: Genetic structure of the purebred domestic dog. Science. 2004, 304: 1160-1164. 10.1126/science.1097406.
Telomerase reverse transcriptase (Telomerase catalytic subunit). http://aug2010.archive.ensembl.org/Canis_familiaris/Gene/Summary?g=ENSCAFG00000010739;r=34:14294248-14312666;t=ENSCAFT00000017081.
Davis WM: ApE- A plasmid Editor. 2012, http://biologylabs.utah.edu/jorgensen/wayned/ape/.
Sequence Assembly and Alignment Software - CodonCode. CodonCode Corporation;2011. http://www.codoncode.com/.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.
Clamp M, Cuff J, Searle SM, Barton GJ: The jalview java alignment editor. Bioinformatics. 2004, 20: 426-427. 10.1093/bioinformatics/btg430.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Felsenstein J: PHYLIP (Phylogeny Inference Package) Version 3.5c. 1993, http://cmgm.stanford.edu/phylip/.
Molecular evolution, phylogenetics and epidemiology. v 1.4. http://tree.bio.ed.ac.uk/software/figtree/.
Charlesworth D, Willis JH: The genetics of inbreeding depression. Nat Rev Genet. 2009, 10: 783-796. 10.1038/nrg2664.
Chase K, Jones P, Martin A, Ostrander EA, Lark KG: Genetic mapping of fixed phenotypes: disease frequency as a breed characteristic. J Hered. 2009, 100 (Suppl 1): S37-41.
Yamaguchi H, Calado RT, Ly H, Kajigaya S, Baerlocher GM, Chanock SJ, Lansdorp PM, Young NS: Mutations in TERT, the gene for telomerase reverse transcriptase, in aplastic anemia. N Engl J Med. 2005, 352: 1413-1424. 10.1056/NEJMoa042980.
Nasir L, Gault E, Campbell S, Veeramalai M, Glibert D, McFarlane R, Munro A, Argyle D: Isolation and expression of the reverse transcriptase component of the Canis familiaris telomerase ribonucleoprotein (dogTERT). Gene. 2004, 336: 105-113. 10.1016/j.gene.2004.03.032.
We wish to thank the Minnesota Supercomputing Institute for computational infrastructure and support, the Animal Cancer Care and Research Program and Masonic Cancer Center for providing funding, Elaine Ostrander and Heidi Parker at the National Human Genome Research Institute of the NIH for providing samples, Curt’s Media and Mary “Andy” Scott for graphics assistance, Matthew Beckman, PhD for advise and constructive discussions, Milcah Scott for technical support, and Evan Tang for scripting and data organization. CAM was the recipient of an Augsburg College URGO program grant.
The authors declare that they have no competing interests.
CAM, JFM and AB designed the experiments. CAM performed the experiments. KATS performed the bioinformatic analyses. CAM, JFM, KATS and AB analyzed the data and prepared the manuscript. All authors read and approved the final manuscript.