Bioinformatic analysis of the effect of SNPs in the pig TERT gene on the structural and functional characteristics of the enzyme to develop new genetic markers of productivity traits
BMC Genomics volume 24, Article number: 487 (2023)
Telomerase reverse transcriptase (TERT) plays a crucial role in synthesizing telomeric repeats that safeguard chromosomes from damage and fusion, thereby maintaining genome stability. Mutations in the TERT gene can lead to a deviation in gene expression, impaired enzyme activity, and, as a result, abnormal telomere shortening. Genetic markers of productivity traits in livestock can be developed based on the TERT gene polymorphism for use in marker-associated selection (MAS). In this study, a bioinformatic-based approach is proposed to evaluate the effect of missense single-nucleotide polymorphisms (SNPs) in the pig TERT gene on enzyme function and structure, with the prospect of developing genetic markers.
A comparative analysis of the coding and amino acid sequences of the pig TERT was performed with corresponding sequences of other species. The distribution of polymorphisms in the pig TERT gene, with respect to the enzyme’s structural-functional domains, was established. A three-dimensional model of the pig TERT structure was obtained through homological modeling. The potential impact of each of the 23 missense SNPs in the pig TERT gene on telomerase function and stability was assessed using predictive bioinformatic tools utilizing data on the amino acid sequence and structure of pig TERT.
According to bioinformatic analysis of 23 missense SNPs of the pig TERT gene, a predictive effect of rs789641834 (TEN domain), rs706045634 (TEN domain), rs325294961 (TRBD domain) and rs705602819 (RTD domain) on the structural and functional parameters of the enzyme was established. These SNPs hold the potential to serve as genetic markers of productivity traits. Therefore, the possibility of their application in MAS should be further evaluated in associative analysis studies.
Telomeres are short tandem nucleotide repeats with the TTAGGG motif in vertebrates [1, 2] located at the ends of chromosomes. The role of telomeres is to protect chromosomes from the destructive action of DNases and prevent their fusion, which is critical for maintaining the stability of the genome. The synthesis of telomeric repeats is carried out by telomerase, which is a ribonucleoprotein complex consisting of telomerase reverse transcriptase (TERT) and telomerase RNA (TER). TER acts as a template for the synthesis of telomeric tandem repeats which is carried out by TERT. The telomerase complex includes a number of other components that interact with the enzyme and are necessary for its functioning in cells . The efficiency of telomerase is determined by the number of telomeric DNA repeats that it can complete at the ends of telomeres . TERT enzyme is encoded by the corresponding TERT gene which is represented by orthologues in a wide range of biological species accordingly to the Ensembl database .
A decrease in the length of telomeres in mammalian somatic cells is associated with their division and is a consequence of end-replication problem, as well as blocking or decreasing telomerase activity. In contrast, in generative and stem cells, the enzyme activity remains high, which allows for maintaining the size of telomeric regions in a number of cell generations [2, 6]. However, under the influence of adverse environmental factors, a reduction in the number of telomeric repeats is often observed, and the length of telomeres in such cases can be considered a molecular indicator of such impact [7,8,9,10,11]. On the other hand, the shortening of the telomere length can obviously be caused by abnormal insufficient telomerase activity, which either initially does not allow the synthesis of the number of telomeric repeats typical for different cells of the organism or their number is not maintained in the generations of proliferating cells. One of the reasons for the latter may be mutations in the TERT and TER genes, some of which can cause both structural and functional changes in the telomerase complex and deviations in the level of expression of these genes [12,13,14,15].
Interest in studying the organization of telomeres, the structure of telomerase and the activity of this enzyme in cells of different tissues of the organism is associated, first, with the solution of fundamental problems of biology, the study of the mechanisms of cell division and apoptosis, the determination of the life span of the organism, and in the applied aspect, with the study of the causes and processes of the appearance cancerous tumors and a number of other diseases, the search for antitumor agents, the treatment of pathologies caused by abnormalities in the structure of telomeric regions, as well as the assessment of the “genetic age” of cloned animals .
In recent years, several studies have emerged that explore the relationship between telomere size, telomerase activity, and key qualities of livestock. These studies have considered such crucial parameters as animal health, productive lifespan, and resistance to stress in chickens , cattle [9, 17, 18], sheep , pigs , horses and donkeys . Breeding progress in improving these qualities using traditional methods is limited by their low level of heritability and difficulty assessing them. At the same time, one of the approaches to accelerate such progress may be the use of molecular markers, which can be the length of telomeric repeats (the coefficient of their heritability for different species is determined in the range of 0.3 – 0.8) [17, 22,23,24] and genetically determined telomerase activity. In the context of the above, the polymorphisms found in TERT could be the basis for the development of genetic markers capable of showing association with key traits of livestock. This thesis is confirmed by studies where TERT is considered a candidate gene related to a number of physiological and pathological processes [25,26,27,28] and identified in QTLs [29, 30].
Most often, the nature of genetic markers and candidate genes is associated with linkage disequilibrium with causative mutations that are the direct reason of phenotypic changes. The search for such causative mutations is associated with the transition to the study of individual polymorphisms and is a rather difficult task [31,32,33]. According to our suggestion, bioinformatic analysis methods can become a tool that facilitates its solution, which provides a preliminary screening of missense mutations and the selection of those that can potentially be causative. A similar bioinformatics-based approach has previously been successfully used to identify candidate disease genes and predict the potential impact of polymorphic variants on human phenotype [34, 35], and therefore can be applied to the analysis of genetic polymorphisms in animals. The use of this bioinformatic approach can precede experimental testing of polymorphisms, thereby reducing the amount of laboratory work required.
Considering pig TERT in our study as an object for the development of genetic markers, one should pay attention to the fact that many SNPs were found in the gene of this species (674 allelic variants of the gene were known for pig TERT according to the Ensembl database) some of them are the result of missense substitutions. Bioinformatic tools can be used for predictive assessment of their influence and for narrowing the range of SNPs that claim to be genetic markers. Bioinformatic analysis makes it possible to calculate the effects of missense variants on the enzyme structure and function. This will help determine the most promising polymorphisms, which can be further considered probable genetic markers and for which it is advisable to conduct association studies for MAS.
This study proposes an analysis of the structure of the pig TERT gene in comparison with other biological species, as well as a prognostic evaluation of the influence of missense SNPs found in the pig TERT gene on the structural and functional characteristics of the enzyme with the prospect of developing genetic markers for MAS of this animal species.
Comparative analysis of the CDS and AAS of pig TERT with other biological species
TERT gene in pigs is mapped on chromosome 16 and includes 17 exons with a transcript size of 5497 bp. The length of pig TERT AAS is 1130 amino acid residues accordingly to the reference sequence from the Ensembl database . In order to gain a more complete idea about the pig TERT gene, a comparative analysis was carried out with the human TERT gene, which is rather well studied, and corresponding orthologues of some animals that are important model or livestock species. The results of the pairwise alignment of the pig TERT CDS and AAS with the corresponding sequences of other biological species indicate the highest degrees of identity and similarity with those species that are most phylogenetically close: cattle, sheep, goat, horse, donkey (Table 1). High degrees of identity and similarity were also found for the pig TERT CDS and AAS with human TERT. The smallest degrees of identity and similarity of pig TERT were observed with the TERT of those biological species that are phylogenetically distant from the pig, such as rodents (mouse, rat), and with the TERT of other common model organisms (African clawed frog, zebrafish, takifugu).
Distribution of the polymorphisms in the pig TERT gene
Information regarding the number and chromosomal localization of the polymorphisms in the pig TERT gene was obtained from the Ensembl database in which known genetic variants are connected with corresponding rsIDs. The distribution of these polymorphisms was then established in relation to the gene regions responsible for encoding distinct structural-functional domains of the TERT enzyme (Table 2). For this analysis, a comparison was made between the pig TERT gene and its orthologues in cattle and humans. The selection of cattle as a comparative species was based on their phylogenetic proximity to pigs, while the choice of humans was because they are the most extensively studied species, and information on the domain structure of human TERT is available in the UniProt database . This approach was implemented through multiple alignment of pig, human and cattle TERT AASs. Consequently, the expected localization of the four main domains (TEN, TRBD, RTD, CTE) and linker regions  in the pig telomerase reverse transcriptase molecule were determined. Subsequently, the pig polymorphisms were categorized into the respective regions of the pig TERT enzyme.
Each of the regions of the TERT gene, corresponding to distinct structural-functional domains of telomerase reverse transcriptase, is characterized by a certain level of polymorphism. Both synonymous and missense SNPs occur. In order to assess the intraspecific variability of pig TERT, a phylogenetic analysis was carried out based on the CDSs of the TERT gene of 12 pig breeds, for which data on the whole genome sequencing are available. The phylogenetic tree presented in Fig. 1 reflects the results of this analysis.
It can be expected that the polymorphisms found in pig TERT can also affect the structure and telomerase activity of the enzyme. In prospect, polymorphisms with certain impact can be used to develop genetic markers for the biological and productive traits of pigs. In this study, only missense variants (23 pig SNPs with corresponding rsIDs from the Ensembl database) are considered because they change AAS and are able to introduce certain structural changes into the enzyme, so the potential impact of missense variants on TERT functional characteristics can be assessed using bioinformatic analysis methods. Data on pig TERT missense SNPs, their rsIDs, localization on the chromosome, corresponding allelic and amino acid variants are shown in Table 3. In addition, it is indicated which allelic variants according to these missense polymorphisms correspond to each of the analyzed 12 pig breeds.
Prognostic evaluation of the effect of pig TERT gene missense SNPs on telomerase function and stability
The results of evaluating the impact of missense SNPs in the pig TERT gene on the function of the enzyme, obtained with sequence-based methods, are shown in Table 4. Combining a set of sequence-based predictive tools with various evaluation algorithms and synthesizing their results makes it possible to identify SNPs that are highly likely to have an effect on TERT enzyme and may be related to the productive traits of animals. Table 4 shows the numerical values calculated by each of the predictive tools; if a certain threshold value is reached, they are considered to “have an effect” or “be deleterious (damaging)”, depending on the tool developer [37,38,39,40,41,42,43,44]. At the same time, it should be taken into account that the missense SNP with a pronounced effect on a protein molecule can be associated with both positive and negative functional consequences. Therefore, this study considers the division of SNPs by the results of sequence-based prediction into those that have a functional “effect” (or possible functional “effect”) and those that are “neutral” instead of characterizing “effective” SNPs as deleterious (damaging) ones.
Those SNPs for which the presence of an “effect” on TERT protein was predicted by the absolute majority of the tools include rs325294961 (R354G), and rs705602819 (R629W); among the SNPs for which such an “effect” is predicted by the vast majority of tools (or the value of the assessment is close to the threshold), rs789641834 (L158M) rs706045634 (R201P), and rs705219838 (T270I) should be considered.
Another approach to assess the effect of missense SNPs on the functional characteristics of TERT is structural analysis. For this purpose, the three-dimensional structure of the pig TERT was obtained by homology modeling on the basis of the corresponding structure of human TERT, which was made possible by a sufficient degree of identity between human and pig AASs (Table 1). Minimized variants of this structure were deposited and are free available (ModelArchive, accession codes: ma-ydbtw, ma-p89hn). The visualization of pig TERT structure is presented in Fig. 2. The localizations of amino acid substitutions corresponding to the considered missense SNPs in the TERT gene are also presented. It allows drawing a conclusion about the possible effect of each of the substitutions on the secondary and tertiary structures of a particular domain of the TERT protein.
The stability of the tertiary structure of the protein is ensured by its folding mechanism . An important characteristic of protein stability is the level of its folding free energy . In this study, the concept of folding free energy (ΔG) considers the difference between the Gibbs free energy between the unfolded (Gu) and the folded (Gf) states (ΔG = Gu–Gf); more positive values of ΔG indicate greater energy stability of the protein . To determine exactly how each of the studied missense SNPs affects the energy stability of the TERT protein, prediction of changes in folding free energy (ΔΔG) was made as the difference between the folding energy of protein with missense variant (ΔGSNP) and reference one (ΔG): ΔΔG = ΔGSNP – ΔGref. (Table 5). Substitutions that are characterized by negative ΔΔG values, thus, lead to a decrease in TERT protein stability, and substitutions for which the calculated value of folding free energy is positive are stabilizing.
Bioinformatic tools based on different algorithms were used to increase the accuracy of predicting the effect of substitutions on the energy stability of the tertiary structure of TERT. For some of the studied SNPs, a consistent score was obtained across all tools. Thus, it is worth paying attention to rs789641834 (L158M), rs325294961 (R354G), and rs324158660 (I516T), for which the corresponding amino acid substitutions are characterized by consensus negative ΔΔG values, which demonstrates their destabilizing effect on TERT (Table 5). In addition, the G227E, T270A, A997T, P1022L, and E1114K substitutions can be destabilizing. Amino acid substitutions R281W and R606L, which correspond to rs330770291 and rs698738374, stabilize the enzyme due to an increase in folding free energy, and they can probably be considered useful SNPs.
This study compared the degrees of identity and similarity of the pig TERT AAS and CDS with other phylogenetic species. As expected, the CDDs and AASs of the phylogenetically close to pig Cetartiodactyla species (cattle, goat, sheep) had the highest degrees of identity and similarity with corresponding pig TERT sequences. At the same time, rather high degrees of identity and similarity were also found for the pig TERT CDS and AAS with human TERT, which is important because the structural and functional features of the human TERT gene and the protein encoded by it are well-studied and can be used in the analysis of pig TERT. It is notable that the length of human TERT AAS is 1132 amino acid residues, which is very close to pig TERT, but a similar AAS length compared to human is due to a number of indels in the TERT gene and does not indicate a greater similarity of pig TERT to human TERT than to Cetartiodactyla species (cattle, goat, sheep).
It is also worth noting that the reference sequence of pig TERT entered the NCBI  database (NCBI Gene ID: 492280) has some differences from this sequence in Ensembl. Thus, according to NCBI data, the reference pig TERT protein consists of 1131 amino acid residues and, compared to the reference sequence in Ensembl, has the amino acid substitution R354G (the relevant rs325294961 is also considered in this study), as well as a single amino acid insertion (P700_P701insA). This may be the reason for some differences regarding the length of the AAS of pig TERT and the localization of individual SNPs when using the NCBI sequence as a reference.
According to the Ensembl data, there are 4 missense SNPs that correspond to the TEN domain of pig TERT (Table 2), while a total of 6 polymorphisms were detected in this gene region. For comparison, 336 polymorphisms (from which 178 missense) of human TERT and 75 polymorphisms (from which 56 missense) of cattle TERT genes are known in the same region. The significant difference in the number of SNPs found in the TERT genes of pigs, humans, and cattle is attributed to the varying degree of knowledge of these species. It should be noted that the influence of synonymous variants can also take place, for example, on the conformation and function of the protein, affecting post-transcriptional processing, or on the level of its synthesis through the mechanisms of mRNA translation [49, 50]. Functional data confirm that the TEN domain is required for telomerase recruitment to telomeres . In humans, mutations in the DAT (Dissociates Activity of Telomerase) region in the TEN domain render the enzyme unable to function in vivo, but telomerase retains some catalytic activity in vitro . This fact indicates the possibility of the influence of missense SNPs localized in the part of the pig TERT gene that encodes the TEN domain of telomerase on the activity of the enzyme.
The gene region that corresponds to the TRBD domain in the pig TERT contains 12 polymorphisms, out of which 5 are missense SNPs. Due to the low level of research on pig TERT, these data are unlikely to reflect the real level of polymorphism of this area in the pig TERT gene because 330 and 63 polymorphisms have been established accordingly for human and cattle TERT in a similar area.
As for the RTD domain, there are a total of 19 SNPs known, with 5 of them being missense SNPs. According to data obtained for human TERT, this domain contains 7 conserved motifs . Mutations in this domain lead to a decrease in enzyme activity . Moreover, for yeast TERT, it was shown that in addition to mutations that disrupt the function of the enzyme and assembly of the telomerase complex [54, 55], there are mutations that lead to an increase in the length of telomeres . These data give reason to expect that some of the missense SNPs found in the region of the pig TERT gene encoding the RTD domain can significantly affect the structure of the enzyme and its activity.
The gene region that corresponds to the CTE domain in pig telomerase is characterized by 5 missense polymorphisms. The CTE is known to participate in telomerase recruitment, and mutations in this domain do not affect telomerase activity in vitro but do not allow the enzyme to maintain telomere length in vivo .
In this study, a predictive assessment of the expected effect of pig TERT gene missense SNPs on telomerase function and stability was performed, for which sequence-based and structure-based methods were used. The results of evaluation by sequence-based methods showed that of all the considered SNPs, missense ones rs325294961, rs705602819, rs789641834, and rs706045634 can have such an effect.
As for rs325294961, this polymorphism involves an amino acid substitution of arginine to glycine (R354G) in the TRBD domain of telomerase, which may affect the interaction of TERT with telomerase RNA. There is the substitution of the cytosine nucleotide with a guanine nucleotide (CGG/GGG) in the first position of the codon. This agrees with the notion that substitutions of the first nucleotide of the codon have the greatest effect on protein structure. As a result of such a substitution, a charged amino acid is often replaced by an amino acid with the opposite charge . In this case, arginine belongs to the basic amino acids and, thanks to the guanidine group, exhibits strong alkalic properties. This amino acid is able to form multiple hydrogen bonds with phosphate groups of nucleic acids. Glycine is a neutral amino acid, which differs significantly from arginine in terms of its chemical properties . Thus, it is obvious that the substitution of arginine with glycine (R354G) can change the structural characteristics of telomerase and, probably, the affinity of the TRBD domain to telomerase RNA, which, in turn, can affect the catalytic activity of the enzyme. This assumption is confirmed by the evaluation of the effect of the R354G substitution in the TRBD domain, obtained using the bioinformatic tools SIFT, PROVEAN, PolyPhen-2, and SNAP2. According to involved software resources, this mutation is defined as having a significant effect on the enzyme. It is also worth noting that rs200191524 in an adjacent position of human TERT leads to the substitution of arginine for glutamine (R358Q), characterized by association with the phenotype according to ClinVar  data (RCV001219035.5).
The rs705602819 polymorphism corresponds to the amino acid substitution R629W in the RTD domain. In the same way as for the R354G, the substitution of arginine to tryptophan is associated with a substitution of a nucleotide in the first position of the codon (CGG/TGG). Indeed, the specified amino acids differ significantly in their chemical properties. Arginine is an alkalic hydrophilic amino acid; tryptophan is an aromatic amino acid that exhibits hydrophobic properties . Predictive evaluation of all used sequence-based tools indicates a significant effect of the R629W substitution on the function of telomerase. Given that the RTD domain is responsible for the enzymatic reverse transcriptase function of telomerase, it can be expected that rs705602819 has the prospect of being used as a genetic marker associated with the activity of the enzyme. This thesis is also supported by the fact that the analyzed SNP in the pig TERT gene corresponds to rs1194223999 in the human TERT gene, which causes an identical amino acid substitution R631W in the analogous position of the enzyme. In addition, the polymorphisms leading to the substitutions R629W in the pigs and R631W in the humans are in a conserved region with similar AASs between different species, suggesting their similar influence. The rs1194223999 of human TERT gene shows a phenotypic effect consistent with ClinVar data (RCV001172450.2, RCV002411664.1) and scientific publications [61, 62]. In turn, current study bioinformatically predicted the functional effect of rs705602819 in pig TERT gene.
As mentioned above, the predictive scores for rs789641834, rs706045634, and rs705219838 by several sequence-based tools also suggest their possible effect on the functional characteristics of the enzyme (Table 4). The first two of these SNPs are located in the gene region corresponding to the TEN domain of telomerase. The rs789641834 is caused by a substitution of the first nucleotide of the codon (CTG/ATG) and is the cause of the corresponding amino acid substitution L158M in the AAS of the enzyme. Leucine is a typical non-polar aliphatic α-amino acid, and methionine is also a non-polar aliphatic α-amino acid but has a bonded sulfur atom that exhibits hydrophobic properties . When leucine is replaced by methionine, the hydrophobicity of the latter may affect the spatial structure of the protein. Another rs706045634 as assessed by the SIFT, PolyPhen-2, and SNAP2 tools also demonstrates the possibility of influencing the functional properties of telomerase, leading to the amino acid substitution R201P in the same TEN domain. As for the last of the mentioned substitutions with possible effect, the rs705219838 is in the Linker region of the TERT enzyme and leads to the substitution T270I.
Comparison of the results of the predicted impact of the missense SNPs obtained by structure-based methods with the results obtained by sequence-based methods reveals certain regularities in the estimates. Thus, results of the sequence-based methods for rs789641834 (L158M) and rs325294961 (R354G) indicate that there is a potential “effect” of these polymorphisms on TERT functions (Table 4). At the same time, these polymorphisms predictably destabilize the structure of the TERT protein (Table 5). This allows us to make a general tentative conclusion about their significant conditionally deleterious effect on telomerase activity.
According to the results obtained using sequence-based methods, amino acid substitution R201P (rs706045634) refers to those that have a certain effect on the functional properties of telomerase. And according to structure-based methods, it is likely to contribute to increasing the stability of the enzyme. It can be assumed that this substitution affects the functional properties of the enzyme by increasing the stability of the tertiary structure of TERT.
Based on sequence-based methods, the substitutions R281W (rs330770291) and R606L (rs698738374) are estimated as neutral. This prediction is consistent with their assessment by structure-based methods, according to which they are likely stabilizing with respect to the tertiary structure of TERT. An increase in the stability of the TERT structure may contribute to a change in the activity of the enzyme.
Thus, changes in telomerase activity associated with the influence of rs789641834, rs325294961, and rs706045634 (based on both sequence and structure prediction) and rs705602819 (based on strong sequence prediction) can result in changes in processes associated with cell proliferation, which, as can be assumed, affect physical parameters.
Given the important biological role of telomerase and the relationship of its activity, as mentioned above, with the manifestation of traits of livestock animals, the TERT gene is considered in this study as a candidate gene that exhibits biological functional impact. The results of bioinformatic evaluation help to conduct a target selection for associative analysis studies of those SNPs for which the effect on the structure and function of the protein is demonstrated in silico. In this case, out of 23 missense polymorphisms analyzed in the pig TERT gene, 4 SNPs (rs789641834, rs325294961, rs706045634, and rs705602819) are the most promising candidate genetic markers.
Noted SNPs should be tested in associative analysis studies to establish their actual effect on animal performance. If the results of such studies demonstrate their significant association with such traits, they can obviously be used in MAS and can be considered direct genetic markers. The mutations corresponding to them can claim to be classified as causative, since they directly affect the functional and structural characteristics of TERT, which, in turn, changes the activity of the enzyme. Those SNPs for which the corresponding amino acid substitutions, according to the results of bioinformatic analysis, did not reveal an effect on the characteristics of TERT, but in the associative analysis performed for them show an effect on productive traits, obviously, can be considered linkage disequilibrium genetic markers.
This approach seems to be useful for searching genetic markers of candidate genes. The TERT gene in this study served as a model object for its demonstration.
Pig productivity depends on health, productive lifespan, resistance to stress, and environmental factors. The use in breeding programs of genetic markers associated with these key qualities of farm animals will improve the selection and economic aspect of pig breeding. Such markers can be developed based on the polymorphism of the telomerase gene, which is involved in the control of important biological processes.
The organization of the pig TERT gene was analyzed in comparison with TERT in other biological species using bioinformatic methods. Significant similarities between CDS and AAS of the pig TERT with those sequences of other species have been established. The comparison was also made with human TERT in the context of the structural-functional domains of the enzyme and the distribution of SNPs. Due to the high similarity between pig and human TERT, a three-dimensional structure of pig TERT was obtained using the homology modeling method.
According to the results of bioinformatic analysis of 23 missense SNPs of the pig telomerase gene, a predictive effect of rs789641834, rs706045634, rs325294961 and rs705602819 on the structural and functional parameters of the enzyme were established. The TERT gene, given the important role of telomerase enzyme, which it encodes, can be considered a candidate gene involved in the control of biological and productive traits. The bioinformatic assessment of these SNPs enables a targeted selection of polymorphisms that may serve as potential genetic markers, thus the possibility of their application in MAS should be further investigated through associative analysis studies.
Analysis of the primary structure of the TERT gene and protein
The coding sequence (CDS) and the amino acid sequence (AAS) corresponding to the pig TERT gene (Ensembl ID: ENSSSCG00000017118, duroc breed) obtained from Ensembl  were used as references. Canonical CDSs and AASs of the TERT of the following organisms were used for comparative analysis and alignments: human (Ensembl ID: ENSG00000164362); model organisms: mouse (Ensembl ID: ENSMUSG00000021611), rat (Ensembl ID: ENSRNOG00000025327), African clawed frog (NCBI Gene ID: 373635), zebrafish (Ensembl ID: ENSDARG00000042637), takifugu (Ensembl ID: ENSTRUG00000014198); agricultural species: cattle (Ensembl ID: ENSBTAG00000012567), sheep (Ensembl ID: ENSOARG00020010728), goat (Ensembl ID: ENSCHIG00000008229), horse (NCBI Gene ID: 100630695, isoform X1), donkey (NCBI Gene ID: 106835880).
Pairwise alignments of CDSs and AASs were performed using the Needleman-Wunsch algorithm  in the EMBOSS Needle . According to the alignment results, the degrees of identity between pigs and other biological species CDSs were determined, as well as the degrees of identity and the degree of similarity between the pigs and AASs of other biological species (the degree of identity in this study refers to the percentage of identical nucleotides or amino acids in CDSs and AASs respectively; the degree of similarity refers to the percentage of amino acids in AASs with alignment score > 0 accordingly to EBLOSUM62 substitution matrix).
SNPs in the pig TERT gene, along with their corresponding rsIDs, were identified using data from the Ensembl database. To evaluate the distribution of polymorphisms across pig TERT domains, multiple alignment of pig TERT AAS was conducted with corresponding sequences of cattle (a phylogenetically related species) and human (whose TERT domain structure is known, sourced from UniProt ID: O14746) sequences. Alignment was performed using MUSCLE algorithm  in MEGA11 software .
For pig breeds with available whole genome sequencing data, a phylogenetic analysis was conducted using the relevant CDSs. The following breeds were included in the analysis: bamei (Ensembl ID: ENSSSCG00050053174), berkshire (Ensembl ID: ENSSSCG00065076605), duroc (reference pig TERT sequence), hampshire (Ensembl ID: ENSSSCG00035007219), jinhua (Ensembl ID: ENSSSCG00060012419), landrace (Ensembl ID: ENSSSCG00045034770), largewhite (Ensembl ID: ENSSSCG00025035362), meishan (Ensembl ID: ENSSSCG00040053333), pietrain (Ensembl ID: ENSSSCG00055017091), rongchang (Ensembl ID: ENSSSCG00030034123), tibetan (Ensembl ID: ENSSSCG00015056815), USMARC (Ensembl ID: ENSSSCG00070009936). The MUSCLE algorithm in MEGA11 software was used to perform a multiple alignment of the TERT CDSs of these breeds. Subsequently, a phylogenetic tree was constructed using the Maximum Likelihood method and the JTT matrix-based model . Missense polymorphisms in pig TERT gene were included in the further analysis because they can lead to amino acid substitutions in AAS and potentially impact the function of the protein encoded by the TERT gene. For all missense polymorphisms, allelic variants corresponding to each of the studied 12 pig breeds were identified.
Assessment of the SNPs using sequence-based methods
This approach involved evaluating the impact of missense substitutions on TERT protein functions using computational tools that utilize TERT AAS as input data. Sequence homology-based tools such as SIFT  and PROVEAN [38, 39] determine the impact of missense substitutions based on the alignment of the input sequence with homologous sequences and assessment of position conservatism, where there are suitable substitutions. A set of homologous sequences that were selected using BLAST  in the non-redundant UniRef90 database  were used as related sequences in SIFT evaluation. PANTHER  is a predictive tool that evaluates position-specific evolutionary preservation as the length of time (in millions of years) a position in protein has been preserved, i.e. determines the conservatism of the position in which the substitutions take place. PolyPhen-2 [41, 42] and SNAP2 [43, 44] are the tools that utilize various characteristics of substitution to assess its functional impact.
Assessment of the SNPs using structure-based methods
This approach involved evaluating the impact of polymorphisms on the stability of the TERT protein based on the data of the three-dimensional structure of the protein. Since no such experimental structure has been established for the pig TERT, homology modeling was used. Homology modeling of the pig TERT was performed using MODELLER 10.3 software . The cryo-EM map of the human telomerase-DNA-TPP1 complex (PDB ID: 7QXA  was used as a template structure for modeling. The obtained as a result of homology modeling three-dimensional pig structure was minimized in the Amber ff14SB force field  according to the built-in protocol in the UCSF Chimera software . Evaluation of the minimized model was carried out using the ProSA-web service [74, 75]. The three-dimensional structure was visualized in free PyMOL software . To assess the influence of SNPs on the energy stability of the protein, bioinformatic tools with different predictive algorithms were used: I-Mutant 3.0 is the support vector machine (SVM)-based tool ; PoPMuSiC is the server that makes predictions using a linear combination of statistical potentials [78,79,80]: mCSM is the tool based on graph-based structural signatures ; SDM2 is the server that uses conformationally-constrained environment-specific substitution tables for prediction ; DDGun is the tool that makes prediction through a linear combination of scores derived from sequence and structural evolutionary features . Pig TERT protein structure minimized in UCSF Chimera, as mentioned above, was used as input in all services except SDM2; for making a prediction using the SDM2 service, a structure with minimization in the GROMOS96 43B1 force field  performed in the Swiss-PdbViewer 4.1.0  was used to reduce the size of the input file following requirements of the SDM2 server.
Availability of data and materials
Models of the structure of Sus scrofa (pig) telomerase reverse transcriptase (TERT) are available in ModelArchive (modelarchive.org) with the accession codes ma-ydbtw, ma-p89hn. Other datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Amino acid sequence
Reverse transcriptase domain
Telomerase RNA component
Telomerase reverse transcriptase
Telomerase RNA-binding domain
Moyzis RK, Buckingham JM, Cram LS, Dani M, Deaven LL, Jones MD, Meyne J, Ratliff RL, Wu JR. A highly conserved repetitive DNA sequence, (TTAGGG)n, present at the telomeres of human chromosomes. Proc Natl Acad Sci U S A. 1988;85(18):6622–6. https://doi.org/10.1073/pnas.85.18.6622.
Fradiani PA, Ascenzioni F, Lavitrano M, Donini P. Telomeres and telomerase activity in pig tissues. Biochimie. 2004;86(1):7–12. https://doi.org/10.1016/j.biochi.2003.11.009.
Nguyen KTTT, Wong JMY. Telomerase Biogenesis and Activities from the Perspective of Its Direct Interacting Partners. Cancers (Basel). 2020;12(6):1679. https://doi.org/10.3390/cancers12061679.
Zvereva MI, Shcherbakova DM, Dontsova OA. Telomerase: structure, functions, and activity regulation. Biochemistry (Mosc). 2010;75(13):1563–83. https://doi.org/10.1134/s0006297910130055.
Cunningham F, Allen JE, Allen J, Alvarez-Jarreta J, Amode MR, Armean IM, Austine-Orimoloye O, Azov AG, Barnes I, Bennett R, Berry A, Bhai J, Bignell A, Billis K, Boddu S, Brooks L, Charkhchi M, Cummins C, Da Rin Fioretto L, Davidson C, Dodiya K, Donaldson S, El Houdaigui B, El Naboulsi T, Fatima R, Giron CG, Genez T, Martinez JG, Guijarro-Clarke C, Gymer A, Hardy M, Hollis Z, Hourlier T, Hunt T, Juettemann T, Kaikala V, Kay M, Lavidas I, Le T, Lemos D, Marugán JC, Mohanan S, Mushtaq A, Naven M, Ogeh DN, Parker A, Parton A, Perry M, Piližota I, Prosovetskaia I, Sakthivel MP, Salam AIA, Schmitt BM, Schuilenburg H, Sheppard D, Pérez-Silva JG, Stark W, Steed E, Sutinen K, Sukumaran R, Sumathipala D, Suner MM, Szpak M, Thormann A, Tricomi FF, Urbina-Gómez D, Veidenberg A, Walsh TA, Walts B, Willhoft N, Winterbottom A, Wass E, Chakiachvili M, Flint B, Frankish A, Giorgetti S, Haggerty L, Hunt SE, IIsley GR, Loveland JE, Martin FJ, Moore B, Mudge JM, Muffato M, Perry E, Ruffier M, Tate J, Thybert D, Trevanion SJ, Dyer S, Harrison PW, Howe KL, Yates AD, Zerbino DR, Flicek P. Ensembl 2022. Nucleic Acids Res. 2022;50(D1):D988–95. https://doi.org/10.1093/nar/gkab1049.
Yadav PS, Wakil AM. Telomerase Structure and Function, Activity and Its Regulation with Emerging Methods of Measurement in Eukaryotes. In: Morrish TA, editor. Telomerase and non-Telomerase Mechanisms of Telomere Maintenance. London: IntechOpen. 2019. https://doi.org/10.5772/intechopen.89506.
Seluanov A, Chen Z, Hine C, Sasahara TH, Ribeiro AA, Catania KC, Presgraves DC, Gorbunova V. Telomerase activity coevolves with body mass not lifespan. Aging Cell. 2007;6(1):45–52. https://doi.org/10.1111/j.1474-9726.2006.00262.x.
Gorbunova V, Seluanov A. Coevolution of telomerase activity and body mass in mammals: from mice to beavers. Mech Ageing Dev. 2009;130(1–2):3–9. https://doi.org/10.1016/j.mad.2008.02.008.
Brown DE, Dechow CD, Liu WS, Harvatine KJ, Ott TL. Hot topic: association of telomere length with age, herd, and culling in lactating Holsteins. J Dairy Sci. 2012;95(11):6384–7. https://doi.org/10.3168/jds.2012-5593.
Bateson M. Cumulative stress in research animals: Telomere attrition as a biomarker in a welfare context? BioEssays. 2016;38(2):201–12. https://doi.org/10.1002/bies.201500127.
Badmus KA, Idrus Z, Meng GY, Sazili AQ, Mamat-Hamidi K. Telomere Length and Regulatory Genes as Novel Stress Biomarkers and Their Diversities in Broiler Chickens (Gallus gallus domesticus) Subjected to Corticosterone Feeding. Animals (Basel). 2021;11(10):2759. https://doi.org/10.3390/ani11102759.
Mitchell JR, Wood E, Collins K. A telomerase component is defective in the human disease dyskeratosis congenita. Nature. 1999;402(6761):551–5. https://doi.org/10.1038/990141.
Selivanova LS, Volganova KS, Abrosimov AY. Telomerase reverse transcriptase (TERT) promoter mutations in the tumors of human endocrine organs: Biological and prognostic value. Arkh Patol. 2016;78(1):62–9. https://doi.org/10.17116/patol201678162-68. (Russian).
Dai J, Cai H, Zhuang Y, Wu Y, Min H, Li J, Shi Y, Gao Q, Yi L. Telomerase gene mutations and telomere length shortening in patients with idiopathic pulmonary fibrosis in a Chinese population. Respirology. 2015;20(1):122–8. https://doi.org/10.1111/resp.12422.
Grill S, Nandakumar J. Molecular mechanisms of telomere biology disorders. J Biol Chem. 2021;296:100064. https://doi.org/10.1074/jbc.REV120.014017.
Burgstaller JP, Brem G. Aging of Cloned Animals: A Mini-Review. Gerontology. 2017;63(5):417–25. https://doi.org/10.1159/000452444.
Seeker LA, Ilska JJ, Psifidi A, Wilbourn RV, Underwood SL, Fairlie J, Holland R, Froy H, Salvo-Chirnside E, Bagnall A, Whitelaw B, Coffey MP, Nussey DH, Banos G. Bovine telomere dynamics and the association between telomere length and productive lifespan. Sci Rep. 2018;8(1):12748. https://doi.org/10.1038/s41598-018-31185-z.
Iannuzzi A, Albarella S, Parma P, Galdiero G, D’Anza E, Pistucci R, Peretti V, Ciotola F. Characterization of telomere length in Agerolese cattle breed, correlating blood and milk samples. Anim Genet. 2022;53(5):676–9. https://doi.org/10.1111/age.13227.
Froy H, Underwood SL, Dorrens J, Seeker LA, Watt K, Wilbourn RV, Pilkington JG, Harrington L, Pemberton JM, Nussey DH. Heritable variation in telomere length predicts mortality in Soay sheep. Proc Natl Acad Sci USA. 2021;118(15):e2020563118. https://doi.org/10.1073/pnas.2020563118.
Russo V, Berardinelli P, Martelli A, Di Giacinto O, Nardinocchi D, Fantasia D, Barboni B. Expression of telomerase reverse transcriptase subunit (TERT) and telomere sizing in pig ovarian follicles. J Histochem Cytochem. 2006;54(4):443–55. https://doi.org/10.1369/jhc.4A6603.2006.
Argyle D, Ellsmore V, Gault EA, Munro AF, Nasir L. Equine telomeres and telomerase in cellular immortalisation and ageing. Mech Ageing Dev. 2003;124(6):759–64. https://doi.org/10.1016/s0047-6374(03)00104-0.
Broer L, Codd V, Nyholt DR, Deelen J, Mangino M, Willemsen G, Albrecht E, Amin N, Beekman M, de Geus EJ, Henders A, Nelson CP, Steves CJ, Wright MJ, de Craen AJ, Isaacs A, Matthews M, Moayyeri A, Montgomery GW, Oostra BA, Vink JM, Spector TD, Slagboom PE, Martin NG, Samani NJ, van Duijn CM, Boomsma DI. Meta-analysis of telomere length in 19,713 subjects reveals high heritability, stronger maternal inheritance and a paternal age effect. Eur J Hum Genet. 2013;21(10):1163–8. https://doi.org/10.1038/ejhg.2012.303.
Eisenberg DT. Inconsistent inheritance of telomere length (TL): is offspring TL more strongly correlated with maternal or paternal TL? Eur J Hum Genet. 2014;22(1):8–9. https://doi.org/10.1038/ejhg.2013.202.
Faul JD, Mitchell CM, Smith JA, Zhao W. Estimating Telomere Length Heritability in an Unrelated Sample of Adults: Is Heritability of Telomere Length Modified by Life Course Socioeconomic Status? Biodemography Soc Biol. 2016;62(1):73–86. https://doi.org/10.1080/19485565.2015.1120645.
Tong Y, Xiang Y, Li B, Bao S, Zhou Y, Yuan W, Ling Y, Hao D, Zhu H, Sun Z. Association between TERT gene polymorphisms and acute myeloid leukemia susceptibility in a Chinese population: a case-control study. Cancer Cell Int. 2020;20:313. https://doi.org/10.1186/s12935-020-01335-3.
Li Y, Cheang I, Zhang Z, Yao W, Zhou Y, Zhang H, Liu Y, Zuo X, Li X, Cao Q. Prognostic Association of TERC, TERT Gene Polymorphism, and Leukocyte Telomere Length in Acute Heart Failure: A Prospective Study. Front Endocrinol (Lausanne). 2021;12:650922. https://doi.org/10.3389/fendo.2021.650922.
Zhang X, Chen Y, Yan D, Han J, Zhu L. TERT Gene rs2736100 and rs2736098 Polymorphisms are Associated with Increased Cancer Risk: A Meta-Analysis. Biochem Genet. 2022;60(1):241–66. https://doi.org/10.1007/s10528-021-10097-0.
Marinaccio J, Micheli E, Udroiu I, Di Nottia M, Carrozzo R, Baranzini N, Grimaldi A, Leone S, Moreno S, Muzzi M, Sgura A. TERT Extra-Telomeric Roles: Antioxidant Activity and Mitochondrial Protection. Int J Mol Sci. 2023;24(5):4450. https://doi.org/10.3390/ijms24054450.
Durán Aguilar M, Román Ponce SI, Ruiz López FJ, González Padilla E, Vásquez Peláez CG, Bagnato A, Strillacci MG. Genome-wide association study for milk somatic cell score in holstein cattle using copy number variation as markers. J Anim Breed Genet. 2017;134(1):49–59.
Lai E, Danner AL, Famula TR, Oberbauer AM. Genome-Wide Association Studies Reveal Susceptibility Loci for Digital Dermatitis in Holstein Cattle. Animals (Basel). 2020;10(11):2009. https://doi.org/10.3390/ani10112009.
Suchocki T, Egger-Danner C, Schwarzenbacher H, Szyda J. Two-stage genome-wide association study for the identification of causal variants underlying hoof disorders in cattle. J Dairy Sci. 2020;103(5):4483–94. https://doi.org/10.3168/jds.2019-17542.
Johnsson M, Jungnickel MK. Evidence for and localization of proposed causative variants in cattle and pig genomes. Genet Sel Evol. 2021;53(1):67. https://doi.org/10.1186/s12711-021-00662-x.
Häfliger IM, Spengeler M, Seefried FR, Drögemüller C. Four novel candidate causal variants for deficient homozygous haplotypes in Holstein cattle. Sci Rep. 2022;12(1):5435. https://doi.org/10.1038/s41598-022-09403-6.
van Driel MA, Brunner HG. Bioinformatics methods for identifying candidate disease genes. Hum Genomics. 2006;2(6):429–32. https://doi.org/10.1186/1479-7364-2-6-429.
Liu Y, Yeung WSB, Chiu PCN, Cao D. Computational approaches for predicting variant impact: An overview from resources, principles to applications. Front Genet. 2022;13:981005. https://doi.org/10.3389/fgene.2022.981005.
UniProt Consortium. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021;49(D1):D480–9. https://doi.org/10.1093/nar/gkaa1100.
Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11(5):863–74. https://doi.org/10.1101/gr.176601.
Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS One. 2012;7(10):e46688. https://doi.org/10.1371/journal.pone.0046688.
Choi Y, Chan AP. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015;31(16):2745–7. https://doi.org/10.1093/bioinformatics/btv195.
Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, Diemer K, Muruganujan A, Narechania A. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13(9):2129–41. https://doi.org/10.1101/gr.772403.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9. https://doi.org/10.1038/nmeth0410-248.
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;Chapter 7:Unit7.20. https://doi.org/10.1002/0471142905.hg0720s76.
Bromberg Y, Rost B. SNAP: predict effect of non-synonymous polymorphisms on function. Nucleic Acids Res. 2007;35(11):3823–35. https://doi.org/10.1093/nar/gkm238.
Hecht M, Bromberg Y, Rost B. Better prediction of functional effects for sequence variants. BMC Genomics. 2015;16(Suppl 8):1. https://doi.org/10.1186/1471-2164-16-S8-S1.
Jaenicke R. Stability and folding of domain proteins. Prog Biophys Mol Biol. 1999;71(2):155–241. https://doi.org/10.1016/s0079-6107(98)00032-7.
Zhang Z, Wang L, Gao Y, Zhang J, Zhenirovskyy M, Alexov E. Predicting folding free energy changes upon single point mutations. Bioinformatics. 2012;28(5):664–71. https://doi.org/10.1093/bioinformatics/bts005.
Quan L, Lv Q, Zhang Y. STRUM: structure-based prediction of protein stability changes upon single-point mutation. Bioinformatics. 2016;32(19):2936–46. https://doi.org/10.1093/bioinformatics/btw361.
National Center for Biotechnology Information (NCBI). Bethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information. 1988. Available from: https://www.ncbi.nlm.nih.gov/. Cited 10 Jul 2022.
Elf J, Nilsson D, Tenson T, Ehrenberg M. Selective charging of tRNA isoacceptors explains patterns of codon usage. Science. 2003;300(5626):1718–22. https://doi.org/10.1126/science.1083811.
Sauna ZE, Kimchi-Sarfaty C. Understanding the contribution of synonymous mutations to human disease. Nat Rev Genet. 2011;12(10):683–91. https://doi.org/10.1038/nrg3051.
Smith EM, Pendlebury DF, Nandakumar J. Structural biology of telomeres and telomerase. Cell Mol Life Sci. 2020;77(1):61–79. https://doi.org/10.1007/s00018-019-03369-x.
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(22):7775–86. https://doi.org/10.1128/MCB.21.22.7775-7786.2001.
Lue NF, Lin YC, Mian IS. A conserved telomerase motif within the catalytic domain of telomerase reverse transcriptase is specifically required for repeat addition processivity. Mol Cell Biol. 2003;23(23):8440–9. https://doi.org/10.1128/MCB.23.23.8440-8449.2003.
Friedman KL, Heit JJ, Long DM, Cech TR. N-terminal domain of yeast telomerase reverse transcriptase: recruitment of Est3p to the telomerase complex. Mol Biol Cell. 2003;14(1):1–13. https://doi.org/10.1091/mbc.e02-06-0327.
Bosoy D, Peng Y, Mian IS, Lue NF. Conserved N-terminal motifs of telomerase reverse transcriptase required for ribonucleoprotein assembly in vivo. J Biol Chem. 2003;278(6):3882–90. https://doi.org/10.1074/jbc.M210645200.
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(6):1201–11. https://doi.org/10.1016/s1097-2765(01)00268-4.
Banik SS, Guo C, Smith AC, Margolis SS, Richardson DA, Tirado CA, Counter CM. C-terminal regions of the human telomerase catalytic subunit essential for in vivo enzyme activity. Mol Cell Biol. 2002;22(17):6234–46. https://doi.org/10.1128/MCB.22.17.6234-6246.2002.
Fricke M, Gerst R, Ibrahim B, Niepmann M, Marz M. Global importance of RNA secondary structures in protein-coding sequences. Bioinformatics. 2019;35(4):579–83. https://doi.org/10.1093/bioinformatics/bty678.
Nelson D, Cox M. Lehninger Principles of Biochemistry. 4th ed. New York: W.H. Freeman and Company; 2005.
Landrum MJ, Lee JM, Benson M, Brown GR, Chao C, Chitipiralla S, Gu B, Hart J, Hoffman D, Jang W, Karapetyan K, Katz K, Liu C, Maddipatla Z, Malheiro A, McDaniel K, Ovetsky M, Riley G, Zhou G, Holmes JB, Kattman BL, Maglott DR. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018;46(D1):D1062–7. https://doi.org/10.1093/nar/gkx1153.
Basel-Vanagaite L, Dokal I, Tamary H, Avigdor A, Garty BZ, Volkov A, Vulliamy T. Expanding the clinical phenotype of autosomal dominant dyskeratosis congenita caused by TERT mutations. Haematologica. 2008;93(6):943–4. https://doi.org/10.3324/haematol.12317.
Pereboeva L, Hubbard M, Goldman FD, Westin ER. Robust DNA Damage Response and Elevated Reactive Oxygen Species in TINF2-Mutated Dyskeratosis Congenita Cells. PLoS One. 2016;11(2):0148793. https://doi.org/10.1371/journal.pone.0148793.
Needleman SB, Wunsch CD. A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970;48(3):443–53. https://doi.org/10.1016/0022-2836(70)90057-4.
Madeira F, Pearce M, Tivey ARN, Basutkar P, Lee J, Edbali O, Madhusoodanan N, Kolesnikov A, Lopez R. Search and sequence analysis tools services from EMBL-EBI in 2022. Nucleic Acids Res. 2022;50(W1):W276–9. https://doi.org/10.1093/nar/gkac240.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113. https://doi.org/10.1186/1471-2105-5-113.
Tamura K, Stecher G, Kumar S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol Biol Evol. 2021;38(7):3022–7. https://doi.org/10.1093/molbev/msab120.
Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8(3):275–82. https://doi.org/10.1093/bioinformatics/8.3.275.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. https://doi.org/10.1016/S0022-2836(05)80360-2.
Suzek BE, Wang Y, Huang H, McGarvey PB, Wu CH, UniProt Consortium. UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches. Bioinformatics. 2015;31(6):926–32. https://doi.org/10.1093/bioinformatics/btu739.
Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234(3):779–815. https://doi.org/10.1006/jmbi.1993.1626.
Sekne Z, Ghanim GE, van Roon AM, Nguyen THD. Structural basis of human telomerase recruitment by TPP1-POT1. Science. 2022;375(6585):1173–6. https://doi.org/10.1126/science.abn6840.
Case DA, Babin V, Berryman JT, Betz RM, Cai Q, Cerutti DS, Cheatham TE III, Darden TA, Duke RE, Gohlke H, Goetz AW, Gusarov S, Homeyer N, Janowski P, Kaus J, Kolossváry I, Kovalenko A, Lee TS, LeGrand S, Luchko T, Luo R, Madej B, Merz KM, Paesani F, Roe DR, Roitberg A, Sagui C, Salomon-Ferrer R, Seabra G, Simmerling CL, Smith W, Swails J, Walker RC, Wang J, Wolf RM, Wu X, Kollman PA. AMBER 14. San Francisco: University of California; 2014.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera–a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12. https://doi.org/10.1002/jcc.20084.
Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35(Web Server issue):W407-10. https://doi.org/10.1093/nar/gkm290.
Sippl MJ. Recognition of errors in three-dimensional structures of proteins. Proteins. 1993;17(4):355–62. https://doi.org/10.1002/prot.340170404.
The PyMOL Molecular Graphics System. Version 2.0 for Windows. Schrödinger, LLC; 2022. https://pymol.org/2/.
Capriotti E, Fariselli P, Casadio R. I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 2005;33(Web Server issue):W306-10. https://doi.org/10.1093/nar/gki375.
Dehouck Y, Grosfils A, Folch B, Gilis D, Bogaerts P, Rooman M. Fast and accurate predictions of protein stability changes upon mutations using statistical potentials and neural networks: PoPMuSiC-2.0. Bioinformatics. 2009;25(19):2537–43. https://doi.org/10.1093/bioinformatics/btp445.
Dehouck Y, Kwasigroch JM, Gilis D, Rooman M. PoPMuSiC 2.1: a web server for the estimation of protein stability changes upon mutation and sequence optimality. BMC Bioinformatics. 2011;12:151. https://doi.org/10.1186/1471-2105-12-151.
Gonnelli G, Rooman M, Dehouck Y. Structure-based mutant stability predictions on proteins of unknown structure. J Biotechnol. 2012;161(3):287–93. https://doi.org/10.1016/j.jbiotec.2012.06.020.
Pires DE, Ascher DB, Blundell TL. mCSM: predicting the effects of mutations in proteins using graph-based signatures. Bioinformatics. 2014;30(3):335–42. https://doi.org/10.1093/bioinformatics/btt691.
Pandurangan AP, Ochoa-Montaño B, Ascher DB, Blundell TL. SDM: a server for predicting effects of mutations on protein stability. Nucleic Acids Res. 2017;45(W1):W229–35. https://doi.org/10.1093/nar/gkx439.
Montanucci L, Capriotti E, Frank Y, Ben-Tal N, Fariselli P. DDGun: an untrained method for the prediction of protein stability changes upon single and multiple point variations. BMC Bioinformatics. 2019;20(Suppl 14):335. https://doi.org/10.1186/s12859-019-2923-1.
van Gunsteren WF, Billeter SR, Eising AA, Hünenberger PH, Krüger P, Mark AE., Scott WRP, Tironi IG. Biomolecular Simulation: The GROMOS96 Manual and User Guide. 1st ed. Zurich-Groningen: Biomos b. v; 1996. https://www.gromos.net/.
Guex N, Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997;18(15):2714–23. https://doi.org/10.1002/elps.1150181505.
The authors are grateful for the possibility to use for analysis UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311.
The research was funded by the National Academy of Agrarian Sciences of Ukraine (grant registration number: 0121U109836).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Peka, M., Balatsky, V., Saienko, A. et al. Bioinformatic analysis of the effect of SNPs in the pig TERT gene on the structural and functional characteristics of the enzyme to develop new genetic markers of productivity traits. BMC Genomics 24, 487 (2023). https://doi.org/10.1186/s12864-023-09592-y