A method for estimating the numbers of synonymous and nonsynonymous substitutions per site

A method for estimating the numbers of synonymous (Ks) and nonsynonymous (Ka) substitutions per site is proposed. The method is based on the Li's (J Mol. Evol. 36:96–99, 1993) and Pamilo and Bianchi's (Mol. Biol. Evol. 10:271–281, 1993) method, but a putative source of bias is solved. It is proposed that the number of synonymous substitutions that are actually transitions or transversions should be computed by separating the twofold degenerate sites into two types of sites, 2S-fold and 2V-fold, where only transitional and transversional substitutions are synonymous, respectively. Kimura's (J. Mol. Evol. 16:111–120, 1980) two-parameter correcting method for multiple substitutions at a site is then applied using the overall observed synonymous transversion frequency to estimate both the numbers of synonymous transversional (Bs) and transitional (As) substitutions per site. This approach, therefore, also minimizes stochastic errors. Computer simulations indicate that the method presented gives more accurate Ks and Ka estimates than the aforementioned methods. Furthermore, the obtention of confidence intervals for divergence estimates by computer simulation is proposed.     

A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

How Common Are Intragene Windows with <Emphasis Type="Italic">K</Emphasis><Subscript>A</Subscript> > <Emphasis Type="Italic">K</Emphasis><Subscript>S</Subscript> Owing to Purifying Selection on Synonymous Mutations?

One method for diagnosing the mode of sequence evolution considers the ratio of nonsynonymous substitutions per nonsynonymous site (K _A) to the corresponding figure for synonymous substitutions (K _S). A ratio (K _A/K _S) greater than unity is taken as evidence for positive selection. This, however, need not necessarily be the case. Notably, there is one instance of a high intragenic K _A/K _S peak, revealed by sliding window analysis and observed in two pairwise comparisons, better accounted for by localised purifying selection on synonymous mutations that affect splicing. Is this example exceptional? To address this we isolate intragenic domains with K _A/K _S > 1 from more than 1000 long mouse-rat orthologues. Approximately one K _A/K _S > 1 peak is found per 12–15 kb of coding sequence. Surprisingly, low synonymous substitution rates underpin more incidences than do high nonsynonymous rates. Several reasons, however, prevent us from supposing that the low synonymous rates reflect purifying selection on synonymous mutations. First, for many peaks, the null that the peak is no higher than expected given the underlying rates of evolution, cannot be rejected. Second, of 18 statistically significant incidences with unusually low K _S values, only 3 are repeatable across independent comparisons. At least two of these are within alternatively spliced exons. We conclude that repeatable statistically significant intragenic domains of low intragenic K _S are rare. As so few K _A/K _S peaks reflect increased rates of protein evolution and so few hold statistical support, we additionally conclude that sliding window analysis to infer domains of positive selection is highly error-prone.     

Mitochondrial DNA polymorphisms as risk factors for Parkinson’s disease and Parkinson’s disease dementia

The activity of complex I of the mitochondrial respiratory chain has been found to be decreased in patients with Parkinsons disease (PD), but no mutations have been identified in genes encoding complex I subunits. Recent studies have suggested that polymorphisms in mitochondrial DNA (mtDNA)-encoded complex I genes (MTND) modify susceptibility to PD. We hypothesize that the risk of PD is conveyed by the total number of nonsynonymous substitutions in the MTND genes in various mtDNA lineages rather than by single mutations. To test this possibility, we determined the number of nonsynonymous substitutions of the seven MTND genes from 183 Finns. The differences in the total number of nonsynonymous substitutions and the nonsynonymous to synonymous substitution rate ratio (K_a/K_s) of MTND genes between the European mtDNA haplogroup clusters (HV, JT, KU, IWX) were analysed by using a statistical approach. Patients with PD (n=238) underwent clinical examination together with mtDNA haplogroup analysis and the clinical features between patient groups defined by the number of nonsynonymous substitutions were compared. Our analysis revealed that the haplogroup clusters HV and KU had a lower average number of amino acid replacements and a lower K_a/K_s ratio in the MTND genes than clusters JT and IWX. Supercluster JTIWX with the highest number of amino acid replacements was more frequent among PD patients and even more frequent among patients with PD who developed dementia. Our results suggest that a relative excess of nonsynonymous mutations in MTND genes in supercluster JTWIX is associated with an increased risk of PD and the disease progression to dementia.     

A genetic component of extinction risk in mammals

Genetic factors may play an important role in species extinction but their actual effect remains poorly understood, particularly because of a strong and potentially masking effect expected from ecological traits. We investigated the role of genetics in mammal extinction taking both ecological and genetic factors into account. As a proxy for the role of genetics we used the ratio of the rates of nonsynonymous (amino acid changing) to synonymous (leaving the amino acid unchanged) nucleotide substitutions, K_a / K_s. Because most nonsynonymous substitutions are likely to be slightly deleterious and thus selected against, this ratio is a measure of the inefficiency of selection: if large (but less than 1), it implies a low efficiency of selection against nonsynonymous mutations. As a result, nonsynonymous mutations may accumulate and thus contribute to extinction. As a proxy for the role of ecology we used body mass W, with which most extinction‐related ecological traits strongly correlate. As a measure of extinction risk we used species’ affiliation with the five levels of extinction threat according to the IUCN Red List of Threatened Species. We calculated K_a / K_s for mitochondrial protein‐coding genes of 211 mammalian species, each of which was characterized by body mass and the level of threat. Using logistic regression analysis, we then constructed a set of logistic regression models of extinction risk on ln(K_a / K_s) and lnW. We found that K_a / K_s and body mass are responsible for a 38% and a 62% increase in extinction risk, respectively. Given that the standard error of these values is 13%, the contribution of genetic factors to extinction risk in mammals is estimated to be one‐quarter to one‐half of the total of ecological and genetic effects. We conclude that the effect of genetics on extinction is significant, though it is almost certainly smaller than the effect of ecological traits. Synthesis Mutation provides the material for evolution. However, most mutations that play a role in evolution are slightly deleterious and thus may contribute to extinction. We assess the role of mitochondrial DNA mutations in mammalian extinction risk and find it to be one‐quarter to one‐half of the total of mutation and body mass effects, where body mass represents an integral measure of extinction‐related ecological traits. Genetic factors may be all the more important, because ecological traits associated with large body mass would both promote and protect from extinction, while mutation accumulation caused by low effective population size seems to have no counterbalance.     

Sequence variation in the outer-surface-protein genes of Borrelia burgdorferi

Borrelia burgdorferi is a spirochete pathogen transmitted among warm- blooded hosts by ixodid ticks. Frequency-dependent selection for variant outer-surface proteins might be expected to arise in this species, since rare variants are more likely to avoid immune surveillance in previously infected hosts. We sequenced the OspA and OspB genes of nine North American strains and compared them with nine strains previously described. For each gene, the mean number of synonymous substitutions per synonymous site and the mean number of nonsynonymous substitutions per nonsynonymous site show only a twofold excess of silent mutations. Synonymous rates vary widely along the OspB protein. Some regions show a significant excess of silent substitutions, while divergence in other regions is constrained by biased base composition or selection. The presence, in antigenically important regions of the protein, of significant variation among strains, as well as evidence for recombination among strains, should be considered in attempts to develop vaccines against this disease.      

Variable Substitution Rates of the 18 Domain Sequences in Artemia Hemoglobin

The Artemia hemoglobin is a dimer comprising two nine-domain covalent polymers in quaternary association. Each polymer is encoded by a gene representing nine successive globin domains which have different sequences and are presumed to have been copied originally from a single-domain gene. Two different polymers exist as the result of a complete duplication of the nine-domain gene, allowing the formation of either homodimers or the heterodimer. The total population size of 18 domains comprising nine corresponding pairs, coupled with the probability that they reflect several hundred million years of evolution in the same lineage, provides a unique model in which the process of gene multiplication can be analyzed. The outcome has important implications for the reliability of local molecular clocks. The two polymers differ from each other at 11.7% of amino acid sites; however when corresponding individual domains are compared between polymers, amino acid substitution fluctuates by a factor of 2.7-fold from lowest to highest. This variation is not obvious at the DNA level: Domain pair identity values fluctuate by 1.3-fold. Identity values are, however, uncorrected for multiple substitutions, and both silent and nonsilent changes are pooled. Therefore, to determine the variability in relative substitution rates at the DNA level, we have used the method of Li (1993, J Mol Evol 36:96–99) to determine estimates of nonsynonymous (K _A ) and synonymous (K _S ) substitutions per site for the nine pairs of domains. As expected, the overall level of silent substitutions (K _S of 56.9%) far exceeded nonsilent substitutions (K _A of 6.7%); however, for corresponding domain pairs, K _A fluctuates by 2.3-fold and K _S by 1.7-fold. The large discrepancies reflected in the expressed protein have accrued within a single lineage and the implication is that divergence dates of different genera based on amino acid sequences, even with well-studied proteins of reasonable size, can be wrong by a factor well in excess of 2. Received: 4 June 1997 / Accepted: 17 December 1997     

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (μ _s P _s)/(μ _v P _v), where μ _s and μ _v denote mutation rate of transitions and transversions, respectively, andP _s andP _v denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions,P _s/P_v ≈ 1 at fourfold degenerate sites, so that thes/v ratio at fourfold degenerate sites is approximately μ _s /μ _v , which is a measure of mutational contribution to transition bias. Similarly, thes/v ratio at nondegenerate sites is also an estimate of μ _s /μ _v if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochromeb (cyt-b) in pocket gophers, thes/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, thes/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the μ _s /μ _v ratio is about 2, theP _s/P_v ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content ons/v ratio of protein-coding genes. Correspondence to: X. Xia     

Evolution of Hypervariable Region 1 of Hepatitis C Virus in Primary Infection

The hypervariable region 1 (HVR-1) of the putative envelope encoding E2 region of hepatitis C virus (HCV) RNA was analyzed in sequential samples from three patients with acute type C hepatitis infected from different sources to address (i) the dynamics of intrahost HCV variability during the primary infection and (ii) the role of host selective pressure in driving viral genetic evolution. HVR-1 sequences from 20 clones per each point in time were analyzed after amplification, cloning, and purification of plasmid DNA from single colonies of transformed cells. The intrasample evolutionary analysis (nonsynonymous mutations per nonsynonymous site [K_a], synonymous mutations per synonymous site [K_s], K_a/K_s ratio, and genetic distances [gd]) documented low gd in early samples (ranging from 2.11 to 7.79%) and a further decrease after seroconversion (from 0 to 4.80%), suggesting that primary HCV infection is an oligoclonal event, and found different levels and dynamics of host pressure in the three cases. The intersample analysis (pairwise comparisons of intrapatient sequences; rK_a, rK_s, rK_a/rK_s ratio, and gd) confirmed the individual features of HCV genetic evolution in the three subjects and pointed to the relative contribution of either neutral evolution or selective forces in driving viral variability, documenting that adaptation of HCV for persistence in vivo follows different routes, probably representing the molecular counterpart of the viral fitness for individual environments.