Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “synonymous codon”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,096 records · Page 61Linked to original sources
PubMed1997

[Molecular evolution of MHC DQA genes. I. The maintenance of interallelic divergence and the influence of GC content on gene structure].

The analyses of the proportion of synonymous and missense nucleotide substitution (PS and PN) in different exons, antigen recognition sites (ARS) and non-ARS of EN2 (NAEN2) of 23 alleles at MHC DQA loci in 7 mammal species gave rise to the following findings. (1) PN was about twice as much as PS in ARS among the alleles at DQA1 of any given species, i.e. 7 alleles at HLA-DQA1 or 8 alleles at IaAa this accords with overdominant selection; (2) PS showed more or less the same as PN in ARS among different loci (DQA1 or DQA2 in different species, or DQA1 and DQA2 in one species) or NAEN2 of all comparative pairs, this conforms the expectation of neutral selection; (3) In exon4 and exon3, not only was the substitution proportion extremely low, but also PS was much higher than PN (the ratio PS over PN is 19.5 in alleles at IaAa of mouse and 4 among alleles at different loci), this coincides obviously with purification selection. The analysis of GC content of MHC DQA showed that its peaks were in the regions corresponding to the middle bulks of some domains, that the highest and constant level was in exon4 and that GC content in the third codon position (GC III content) associates inversely with PS. These results indicate that the specified maintenance mechanisms of interallelic diversity relevant to their functions exist in given exons corresponding to some domains of the same MHC DQA locus and GC III content is an important factor in keeping the structure and function of gene under selection constraint. The method for estimating nucleotide substitution proportion was modified.

Alleles↗
PubMed2005

Using purine skews to predict genes in AT-rich poxviruses.

BACKGROUND: Clusters or runs of purines on the mRNA synonymous strand have been found in many different organisms including orthopoxviruses. The purine bias that is exhibited by these clusters can be observed using a purine skew and in the case of poxviruses, these skews can be used to help determine the coding strand of a particular segment of the genome. Combined with previous findings that minor ORFs have lower than average aspartate and glutamate composition and higher than average serine composition, purine content can be used to predict the likelihood of a poxvirus ORF being a "real gene". RESULTS: Using purine skews and a "quality" measure designed to incorporate previous findings about minor ORFs, we have found that in our training case (vaccinia virus strain Copenhagen), 59 of 65 minor (small and unlikely to be a real genes) ORFs were correctly classified as being minor. Of the 201 major (large and likely to be real genes) vaccinia ORFs, 192 were correctly classified as being major. Performing a similar analysis with the entomopoxvirus amsacta moorei (AMEV), it was found that 4 major ORFs were incorrectly classified as minor and 9 minor ORFs were incorrectly classified as major. The purine abundance observed for major ORFs in vaccinia virus was found to stem primarily from the first codon position with both the second and third codon positions containing roughly equal amounts of purines and pyrimidines. CONCLUSION: Purine skews and a "quality" measure can be used to predict functional ORFs and purine skews in particular can be used to determine which of two overlapping ORFs is most likely to be the real gene if neither of the two ORFs has orthologs in other poxviruses.

Aspartic Acid↗
PubMed2006

Evaluation of six methods for estimating synonymous and nonsynonymous substitution rates.

Methods for estimating synonymous and nonsynonymous substitution rates among protein-coding sequences adopt different mutation (substitution) models with subtle yet significant differences, which lead to different estimates of evolutionary information. Little attention has been devoted to the comparison of methods for obtaining reliable estimates since the amount of sequence variations within targeted datasets is always unpredictable. To our knowledge, there is little information available in literature about evaluation of these different methods. In this study, we compared six widely used methods and provided with evaluation results using simulated sequences. The results indicate that incorporating sequence features (such as transition/transversion bias and nucleotide/codon frequency bias) into methods could yield better performance. We recommend that conclusions related to or derived from Ka and Ks analyses should not be readily drawn only according to results from one method.

Codon↗
PubMed2005

Novel Plasmodium vivax dhfr alleles from the Indonesian Archipelago and Papua New Guinea: association with pyrimethamine resistance determined by a Saccharomyces cerevisiae expression system.

In plasmodia, the dihydrofolate reductase (DHFR) enzyme is the target of the pyrimethamine component of sulfadoxine-pyrimethamine (S/P). Plasmodium vivax infections are not treated intentionally with antifolates. However, outside Africa, coinfections with Plasmodium falciparum and P. vivax are common, and P. vivax infections are often exposed to S/P. Cloning of the P. vivax dhfr gene has allowed molecular comparisons of dhfr alleles from different regions. Examination of the dhfr locus from a few locations has identified a very diverse set of alleles and showed that mutant alleles of the vivax dhfr gene are prevalent in Southeast Asia where S/P has been used extensively. We have surveyed patient isolates from six locations in Indonesia and two locations in Papua New Guinea. We sequenced P. vivax dhfr alleles from 114 patient samples and identified 24 different alleles that differed from the wild type by synonymous and nonsynonymous point mutations, insertions, or deletions. Most importantly, five alleles that carried four or more nonsynonymous mutations were identified. Only one of these highly mutant alleles had been previously observed, and all carried the 57L and 117T mutations. P. vivax cannot be cultured continuously, so we used a yeast assay system to determine in vitro sensitivity to pyrimethamine for a subset of the alleles. Alleles with four nonsynonymous mutations conferred very high levels of resistance to pyrimethamine. This study expands significantly the total number of novel dhfr alleles now identified from P. vivax and provides a foundation for understanding how antifolate resistance arises and spreads in natural P. vivax populations.

Alleles↗
PubMed2002

Different rates of nucleotide substitutions in Wolbachia endosymbionts of arthropods and nematodes: arms race or host shifts?

The genus Wolbachia encompasses intracellular bacteria found in arthropods and in filarial nematodes. In arthropods, Wolbachia is primarily a reproductive parasite and shows relatively frequent horizontal transfer between host species, while in nematodes it appears to be a mutualist and is strictly vertically transmitted. We can expect that different selective pressures are acting on their genomes. Here we present an analysis of three Wolbachia genes, wsp, ftsZ and dnaA. In wsp of arthropod Wolbachia, an excess of non-synonymous substitutions was observed, providing evidence for positive selection. In nematode Wolbachia, no evidence for positive selection was found. Pressure for amino acid variation in wsp of arthropod Wolbachia could derive either from an arms race with the host or from the occurrence of more frequent hosts shifts due to horizontal transmission. In nematode Wolbachia, the lack of positively selected sites could result from the absence of an arms race, or from the homogeneity of the biochemical environment they exist in (ensured by strict vertical transmission). In ftsZ minor differences in substitution patterns were observed between arthropod and nematode Wolbachia, only in the 3'-portion of the gene. dnaA showed comparable patterns of variation in both lineages, with evidence for strong conservation.

Animals↗
PubMed2005

DNA sequencing of CREBBP demonstrates mutations in 56% of patients with Rubinstein-Taybi syndrome (RSTS) and in another patient with incomplete RSTS.

Rubinstein-Taybi syndrome (RSTS) is a distinct dominant disorder characterized by short stature, typical face, broad angulated thumbs and halluces, and mental retardation. The RSTS can be caused by chromosomal microdeletions and molecular mutations in the CREBBP gene; however, relatively few mutations have been reported to date. Here, we aimed to determine the rate of point mutations and other small molecular lesions in true RSTS and possible mild variants, by using genomic DNA sequencing. A consecutive series of patients including 17 patients from our previous study was investigated. We identified 19 causative mutations of CREBBP in a total of 45 patients representing three different diagnostic groups: (a) 17 mutations in 30 patients with unequivocal RSTS (detection rate 56.6%), (b) two mutations in eight patients with features suggestive of RSTS ("moderate or incomplete RSTS", detection rate 25%), and (c) no mutation in seven patients with undiagnosed syndromes and isolated features of RSTS. In general, the mutations were distributed without hot spots and most were unique; however, three recurrent mutations (R370X, R1664H, and N1978S) were identified. Furthermore, we detected 15 different intragenic polymorphisms, including two non-synonymous coding polymorphisms, L551I and Q2208H. We report not only the highest detection rate (56.6%) of CREBBP mutations in patients with RSTS to date, but also the second missense mutation (N1978S) in a patient with moderate or incomplete RSTS. Previous studies have identified cytogenetic deletions in the CREBBP gene in eight to 12% of patients and very recently, Roelfsema et al. reported EP300 gene mutations in three of 92 (3.3%) patients with either true RSTS or different syndromes resembling RSTS. Our 56.6% detection rate of molecular mutations in CREBBP in patients with unequivocal RSTS supports the new concept that RSTS is a genetically heterogeneous disorder and furthermore, indicates that RSTS may be caused by gene/s other than CREBBP in up to 30% of cases.

Adolescent↗
PubMed1982

Evolutionary changes of nucleotide sequences of papova viruses BKV and SV40: they are possibly hybrids.

Complete nucleotide sequences were compared between papova viruses BKV and SV40 and the degrees of sequence divergences were compared between structurally and/or functionally different segments or genes in details. It was shown that the rate of synonymous substitution is not only very high but also approximately uniform among different genes in these viruses as in eukaryotic genes examined to date. While all the non-coding regions including the intron showed marked sequence preservation which is in sharp contrasted with the case of eukaryotic genes where the large bulk of non-coding regions evolve at a rate as rapidly as that of synonymous substitution. It is remarkable that a long continuous stretch of sequence including the putative VPX gene and a 5' half of VP2 gene showed strong homology between BKV and SV40. A close examination of the pattern of base substitutions revealed that this unusual homology was derived by recombination between the two viruses during their evolution. On the basis of the pattern of base substitutions and the bias in code word utilization, we also showed that the putative VPX gene actually could code for a functional polypeptide. In papova viruses, the 3' terminal sequence of VP2/3 gene overlaps with the 5' terminal sequence of VPI gene. The pattern of base substitutions in the overlapping segment was examined in detail in comparison with those in the non-overlapping portions of VP2/3 and VP1 genes. It was shown that the evolutionary mode of the overlapping genes is in good agreement with our previous prediction.

BK Virus↗
PubMed2006

Spatial covariation of mutation and nonsynonymous substitution rates in vertebrate mitochondrial genomes.

Mitochondrial genomes encode fundamental subunits of the basic energy producing machinery of eukaryotic cells that are under strong functional constraint. Paradoxically, these genes evolve rapidly in general, and there is substantial variation in evolutionary rates among genes within genomes. In order to investigate spatial variation in selection intensity, we conducted tests of neutrality using ratios of synonymous to nonsynonymous substitutions (dN/dS = omega) on numerous protein gene segments from fishes and mammals. Values of omega were very low for nearly all genomic regions. However, values of both omega and dN varied in a clinal pattern with increasing distance from the light-strand origin of replication. Spatial heterogeneity of nonsynonymous substitution rates exhibits a significantly positive correlation with variation in mutation rates that are related to the mode of mitochondrial DNA replication. The finding that nonsynonymous substitution rates are proportional to mutation rates is expected if a majority of substitutions are selectively neutral or slightly deleterious. Spatial patterns of among-gene variation in nonsynonymous rates were highly similar between fishes and mammals, suggesting that forces governing mitochondrial gene evolution have remained relatively constant over 450 Myr of vertebrate evolution. Conservation of substitution patterns despite major shifts in thermal habit and metabolic demands among taxa implicates a conserved replication mechanism controlling relative mutation rates as a major determinant of mitochondrial protein evolution.

Animals↗
PubMedDate not supplied

Different levels of variability in subtypes 1b and 4a of hepatitis C viruses.

We performed genetic and phenic analyses to evaluate nucleotide and amino-acid sequences of the amino-terminus of the E1 protein of HCV genotype 1b (extracted from databank) and 4a (characterised in this study). The non-synonymous (ka) mutation analysis demonstrated that the genome of genotype 1b was not saturated by variations, with a rate of transition/transversion (s/v) of 1.5, which is similar to the expected ratio (i.e., 2.0). The s/v ratio in genotype 4a isolates was lower (0.98), indicating saturation due long-term variability. Moreover, the genotype 1b sequences showed a higher number of ka mutations (s+v) (mean of 2.8 per sequence) than genotype 4a (mean of 1.5). The introduction of ka mutations resulted in a higher degree of amino acid variability in genotype 4a. In the genome of genotype 1b, each nucleotide mutation introduced new amino acids, with a Granthan distance of 3.35-42.5, whereas for genotype 4a the distances ranged from 48.8 to 102.1. The phenic analysis also indicated different and complex patterns of amino-acid substitution. Finally, diverse isoelectric points and hydrophobicity were predicted for the two genotypes, with a higher acidity for genotype 4a E1 proteins.

Amino Acid Substitution↗
PubMed2001

Genetic analysis of recent Taiwanese isolates of a variant of coxsackievirus A24.

Epidemics of acute hemorrhagic conjunctivitis (AHC) caused by a variant of coxsackievirus A24 (CA24v) reappeared in Taiwan in 1990 and 1994, following the first two epidemics of 1985--86 and 1988--89. To analyze the genetic diversity of recent CA24v in Taiwan, 7 Taiwanese strains isolated during the 1990--94 period were studied together with one Japanese and two Thai strains isolated in 1993. A fragment of 674 nucleotides between the carboxy terminal 3A and the amino terminal 3D polymerase, including the entire 3C protease (3C(pro)), was amplified by a reverse transcription-polymerase chain reaction (RT-PCR) and the nucleotide sequences were determined. In the 549 nucleotides (183 amino acids) of the entire 3C(pro), we found nucleotide differences at 80 positions between 10 strains and the prototype strain, EH24/70, one of the earliest strains of CA24v. Most of the nucleotide changes were synonymous substitutions and only nine amino acid changes were found. The nucleotide sequence homologies among 71 strains worldwide were 88-100%. These 71 nucleotide sequences were then analyzed by Neighbor-joining method and phylogenetically separated into three distinct genotypes. Genotype I consisted of early strains isolated in 1970--71 from Singapore and Hong Kong. Genotype II included isolates from Singapore and Thailand obtained in 1975. Genotype III comprised strains from the eastern hemisphere isolated in 1985--94 from Japan, Taiwan, China, Hong Kong, Thailand, Singapore, Pakistan and Ghana. They were further divided chronologically into six clusters. The recent isolates from Taiwan obtained in 1985/1986, 1988/1989 and 1990--94 were classified into genotype III Clusters 1, 5, and 6 respectively. The evolutionary rate was re-estimated to be 3 x 10(- 3) 30 years after the emergence of the virus.

Amino Acid Sequence↗
PubMed2002

Phosphoglucose isomerases of hagfish, zebrafish, gray mullet, toad, and snake, with reference to the evolution of the genes in vertebrates.

Phosphoglucose isomerase (PGI) is a protein with multiple functions. To infer its structure changes and evolution in vertebrates, we cloned cDNAs encoding PGI genes from hagfish (Paramyxine yangi), gray mullet (Mugil cephalus), zebrafish (Danio rerio), toad (Bufo melanosticus), and snake (Boiga kraepelini). Only one PGI gene was cloned in each of hagfish, toad, and snake, but two PGI genes were found in zebrafish and gray mullet, respectively. The PGI of hagfish encodes 554 amino acids, in contrast to the PGIs of bonyfishes, toad, and snake which encode 553 amino acids and the PGIs of mammals which encode 558 amino acids. Among 558 aligned amino acid sites, there are 314 sites (56.27%) totally conserved. To see if diversifying selection acts on PGI amino acids of vertebrates, we calculated the pairwise ratio of nonsynonymous versus synonymous substitution per site (Ka/Ks) and the ratio of radical amino acid changes versus conservative amino acid changes per sites (dR/dC) between PGI sequences. The average pairwise ratio between nonsynonymous substitutions per nucleotide (Ka) and synonymous substitutions per nucleotide (Ks) among vertebrate PGI sequences equals 0.047 +/- 0.019. The average pairwise ratio between radical amino acid changes and conservative amino acid changes (dR/dC) among the vertebrate PGIs equal 0.938 +/- 0.158 for charge changes, 0.558 +/- 0.085 for polarity changes, and 0.465 +/- 0.0714 when both polarity and volume are considered. There is no amino acid within the vertebrate PGIs under diversifying selection as analyzed by the method of Yang et al. (2000b). The results suggest that the present vertebrate PGIs are at evolutionary stasis and are being subjected to intense purifying selection. The purifying selection is to maintain polarity and volume of the protein but not the charge groups of amino acids. Phylogenetic analysis reveals that vertebrate PGIs can be classified into three major groups: the mammalian, amphibian-reptilian, and teleostean PGIs. The gene tree suggests that the gene duplication event of PGI in bonyfishes occurred before diversification of Acanthopterygii but after the split of bonyfishes and tetrapods. The evolution of multiple functions of PGI is discussed.

Amino Acid Sequence↗
PubMed1998

Characterizing sequence variation in the VP1 capsid proteins of foot and mouth disease virus (serotype 0) with respect to virion structure.

The VP1 capsid protein of foot and mouth disease virus (FMDV) is highly polymorphic and contains several of the major immunogenic sites important to effective antibody neutralization and subsequent viral clearance by the immune system. Whether this high level of polymorphism is of adaptive value to the virus remains unknown. In this study we examined sequence data from a set of 55 isolates in order to establish the nature of selective pressures acting on this gene. Using the known molecular structure of VP1, the rates and ratios of different types of nonsynonymous and synonymous changes were compared between different parts of the protein. All parts of the protein are subject to purifying selection, but this is greatest amongst those amino acid residues within beta-strands and is significantly reduced at residues exposed on the capsid surface, which include those residues demonstrated by previous mutational analyses to permit the virus to escape from monoclonal antibody binding. The ratios of nonsynonymous substitution resulting in various forms of physicochemically radical and conserved amino acid change were shown to be largely equal throughout these different parts of the protein. There was a consistently higher level of nonsynonymous and charge radical sites in those regions of the gene coding for residues exposed on the outer surface of the capsid and a marked difference in the use of amino acids between surface and nonsurface regions of the protein. However, the analysis is consistent with the hypothesis that the observed sequence variation arises where it is least likely to be disruptive to the higher-order structure of the protein and is not necessarily due to positive Darwinian selection.

Amino Acid Sequence↗
PubMed2005

Size of the protein-coding genome and rate of molecular evolution.

In diploid populations of size N, there will be 2 Nmu mutations per nucleotide (nt) site (or per locus) per generation (mu stands for mutation rate). If either the population or the coding genome double in size, one expects 4 Nmu mutations. What is important is not the population size per se but the number of genes (coding sites), the two being often interconverted. Here we compared the total physical length of protein-coding genomes (n) with the corresponding absolute rates of synonymous substitution (K(S)), an empirical neutral reference. In the classical occupancy problem and in the coupons collector (CC) problem, n was expressed as the mean rate of change (K(CC)). Despite inherently very low power of the approaches involving averaging of rates, the mode of molecular evolution of the total size phenotype of the coding genome could be evidenced through differences between the genomic estimates of K(CC) [K(CC)=1/(ln n + 0.57721) n] and rate of molecular evolution, K(S). We found that (1) the estimates of n and K(S) are reciprocally correlated across taxa (r=0.812; p<< 0.001); (2) the gamete-cell division hypothesis (Chang et al. Proc Natl Acad Sci USA 91:827-831, 1994) can be confirmed independently in terms of K(CC)/K(S) ratios; (3) the time scale of molecular evolution changes with change in mutation rate, as previously shown by Takahata (Proc Natl Acad Sci USA 87:2419-2423, 1990), Takahata et al. (Genetics 130:925-938, 1992), and Vekemans and Slatkin (Genetics 137:1157-1165, 1994); (4) the generation time and population size (Lynch and Conery, Science 302:1401-1404, 2003) effects left their "signatures" at the level of the size phenotype of the protein-coding genome.

Codon↗
PubMed2006

Strand bias in complementary single-nucleotide polymorphisms of transcribed human sequences: evidence for functional effects of synonymous polymorphisms.

BACKGROUND: Complementary single-nucleotide polymorphisms (SNPs) may not be distributed equally between two DNA strands if the strands are functionally distinct, such as in transcribed genes. In introns, an excess of A<-->G over the complementary C<-->T substitutions had previously been found and attributed to transcription-coupled repair (TCR), demonstrating the valuable functional clues that can be obtained by studying such asymmetry. Here we studied asymmetry of human synonymous SNPs (sSNPs) in the fourfold degenerate (FFD) sites as compared to intronic SNPs (iSNPs). RESULTS: The identities of the ancestral bases and the direction of mutations were inferred from human-chimpanzee genomic alignment. After correction for background nucleotide composition, excess of A-->G over the complementary T-->C polymorphisms, which was observed previously and can be explained by TCR, was confirmed in FFD SNPs and iSNPs. However, when SNPs were separately examined according to whether they mapped to a CpG dinucleotide or not, an excess of C-->T over G-->A polymorphisms was found in non-CpG site FFD SNPs but was absent from iSNPs and CpG site FFD SNPs. CONCLUSION: The genome-wide discrepancy of human FFD SNPs provides novel evidence for widespread selective pressure due to functional effects of sSNPs. The similar asymmetry pattern of FFD SNPs and iSNPs that map to a CpG can be explained by transcription-coupled mechanisms, including TCR and transcription-coupled mutation. Because of the hypermutability of CpG sites, more CpG site FFD SNPs are relatively younger and have confronted less selection effect than non-CpG FFD SNPs, which can explain the asymmetric discrepancy of CpG site FFD SNPs vs. non-CpG site FFD SNPs.

Algorithms↗
PubMed2005

Selected DNA repair polymorphisms and gastric cancer in Poland.

Impaired DNA repair capacity may adversely affect cancer risk, particularly in subjects exposed to DNA damaging carcinogens, as found in tobacco smoke, or among subjects deficient for protective factors, as found in fruits and vegetables. We studied tobacco use, fruit and vegetable intake, and common non-synonymous single nucleotide polymorphisms in four DNA repair genes in relation to gastric cancer risk, in a population-based, case-control study of 281 incident gastric cancer cases and 390 controls, in Warsaw, Poland. Multivariate logistic regression analysis was performed to calculate odds ratios (OR) and 95% confidence intervals (CI). Increased risks of gastric cancer were found for smokers (OR=3.1, CI=1.9-5.1 for pack-years>or=40 versus never smokers) and subjects with low fruit intake (OR=2.2, CI=1.3-3.6 for 1st versus 4th quartile); risk associated with vegetable intake was not statistically significant. Allele frequencies among the controls were consistent with those previously reported for the 5 polymorphisms studied: XRCC1-Arg399Gln, XPD-Lys751Gln, MGMT-Ile143Val, Leu84Phe, and XRCC3-Thr241Met. None of the studied polymorphisms were independently associated with gastric cancer risk. Smoking-associated risks, however, were greatest for carriers of the XRCC1-399 ArgArg genotype (Pinteraction=0.004). Risks associated with low intake of fruits or vegetables tended to be modified by selected polymorphisms in XRCC1, XPD and MGMT (Pinteraction=0.1-0.2). Risk modification was not found for the other repair polymorphisms. Selected DNA repair polymorphisms did not have independent effects on gastric cancer risk; however, they may modify smoking- and probably diet-related risks for this disease. These results need replication in larger epidemiological studies of gastric cancer.

Adult↗
PubMed2006

Silent mutations in the gene encoding the p53 protein are preferentially located in conserved amino acid positions and splicing enhancers.

The last release of p53 somatic mutation database contains more than 20,000 of mutation among which 951 are silent (synonymous). This striking amount of silent mutations is much more than what would be expected if synonymous mutations were effectively neutral. The prevalent explanation to reconcile this vast amount of silent mutations with the neutral expectation is that they are just the subproduct of the hypermutability process that affect cancer cells. Some evidences have been presented in this direction, and the explanation has been taken as granted. Assuming that silent mutations are effectively neutral has major implication in the investigation of mutational processes that affect the gene encoding the p53 protein, since on the basis of this assumption they are considered the Null hypothesis, for instance for measuring and comparing among tissues the endogenous mutability. From this it follows that determining whether silent mutations in the p53 gene, and in all disease genes in general, are or not basically mutational noise, is of paramount importance. In this paper we readdress this topic by testing whether there is a relationship between the spatial distribution of silent mutations inside the p53 gene and functional significant features of the gene. For this purpose we divided the population of silent mutations in three groups: those that are found accompanied by other mutations (doublets and multiplest), those that were isolated as singlets, but the same mutation was also isolated as being part of a doublet (or multiplet) in another individual. And the last group is composed by those that were always found as singlets and never as being part of a doublet or a multiplet. This last group was expected to be enriched in functionally significant silent mutations. We found that all silent mutations, but particularly those of the last group, are preferentially located in conserved amino acid positions (i.e. functionally important amino acids) and also tend to be located inside suspected splicing enhancers. Noteworthy, this association remains even after eliminating the possible contribution of mutation hotspots. Besides, we present additional evidence in the direction that these putative splicing enhancers are real functional enhancers.

Amino Acid Sequence↗