High proportions of deleterious polymorphisms in constrained human genes

Previous studies on human mitochondrial genomes showed that the ratio of intra-specific diversities at nonsynonymous-to-synonymous positions was two to ten times higher than the ratio of interspecific divergences at these positions, suggesting an excess of slightly deleterious nonsynonymous polymorphisms. However, such an overabundance of nonsynonymous single nucleotide polymorphisms (SNPs) was not found in human nuclear genomes. Here, genome-wide estimates using >14,000 human-chimp nuclear genes and 1 million SNPs from four human genomes showed a significant proportion of deleterious nonsynonymous SNPs (～ 15%). Importantly, this study reveals a negative correlation between the magnitude of selection pressure and the proportion of deleterious SNPs on human genes. The proportion of deleterious amino acid replacement polymorphisms is 3.5 times higher in genes under high purifying selection compared with that in less constrained genes (28% vs. 8%). These results are explained by differences in the extent of contribution of mildly deleterious mutations to diversity and substitution.     

The genomic rate of adaptive amino acid substitution in Drosophila

The proportion of amino acid substitutions driven by adaptive evolution can potentially be estimated from polymorphism and divergence data by an extension of the McDonald-Kreitman test. We have developed a maximum-likelihood method to do this and have applied our method to several data sets from three Drosophila species: D. melanogaster, D. simulans, and D. yakuba. The estimated number of adaptive substitutions per codon is not uniformly distributed among genes, but follows a leptokurtic distribution. However, the proportion of amino acid substitutions fixed by adaptive evolution seems to be remarkably constant across the genome (i.e., the proportion of amino acid substitutions that are adaptive appears to be the same in fast-evolving and slow-evolving genes; fast-evolving genes have higher numbers of both adaptive and neutral substitutions). Our estimates do not seem to be significantly biased by selection on synonymous codon use or by the assumption of independence among sites. Nevertheless, an accurate estimate is hampered by the existence of slightly deleterious mutations and variations in effective population size. The analysis of several Drosophila data sets suggests that approximately 25% +/- 20% of amino acid substitutions were driven by positive selection in the divergence between D. simulans and D. yakuba.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Systematic Mapping of Protein Mutational Space by Prolonged Drift Reveals the Deleterious Effects of Seemingly Neutral Mutations

Systematic mappings of the effects of protein mutations are becoming increasingly popular. Unexpectedly, these experiments often find that proteins are tolerant to most amino acid substitutions, including substitutions in positions that are highly conserved in nature. To obtain a more realistic distribution of the effects of protein mutations, we applied a laboratory drift comprising 17 rounds of random mutagenesis and selection of M.HaeIII, a DNA methyltransferase. During this drift, multiple mutations gradually accumulated. Deep sequencing of the drifted gene ensembles allowed determination of the relative effects of all possible single nucleotide mutations. Despite being averaged across many different genetic backgrounds, about 67% of all nonsynonymous, missense mutations were evidently deleterious, and an additional 16% were likely to be deleterious. In the early generations, the frequency of most deleterious mutations remained high. However, by the 17th generation, their frequency was consistently reduced, and those remaining were accepted alongside compensatory mutations. The tolerance to mutations measured in this laboratory drift correlated with sequence exchanges seen in M.HaeIII’s natural orthologs. The biophysical constraints dictating purging in nature and in this laboratory drift also seemed to overlap. Our experiment therefore provides an improved method for measuring the effects of protein mutations that more closely replicates the natural evolutionary forces, and thereby a more realistic view of the mutational space of proteins.     

Quantifying the slightly deleterious mutation model of molecular evolution

We have attempted to quantify the frequency and effects of slightly deleterious mutations (SDMs), those that have selective effects close to the reciprocal of the effective population size of a species, by comparing the level of selective constraint in protein-coding genes of related species that have different present-day effective population sizes. In our two comparisons, the species with the smaller effective population size showed lower constraint, implying that SDMs had become fixed. The fixation of SDMs was supported by the observation of a higher fraction of radical to conservative amino acid substitutions in species with smaller effective population sizes. The fraction of strongly deleterious mutations (which rarely become fixed) is >70% in most species. Only approximately 10% or fewer of mutations seem to behave as SDMs, but SDMs could comprise a substantial fraction of mutations in protein-coding genes that have a chance of becoming fixed between species.     

Human-specific amino acid changes found in 103 protein-coding genes

We humans have many characteristics that are different from those of the great apes. These human-specific characters must have arisen through mutations accumulated in the genome of our direct ancestor after the divergence of the last common ancestor with chimpanzee. Gene trees of human and great apes are necessary for extracting these human-specific genetic changes. We conducted a systematic analysis of 103 protein-coding genes for human, chimpanzee, gorilla, and orangutan. Nucleotide sequences for 18 genes were newly determined for this study, and those for the remaining genes were retrieved from the DDBJ/EMBL/GenBank database. The total number of amino acid changes in the human lineage was 147 for 26,199 codons (0.56%). The total number of amino acid changes in the human genome was, thus, estimated to be about 80,000. We applied the acceleration index test and Fisher's synonymous/nonsynonymous exact test for each gene tree to detect any human-specific enhancement of amino acid changes compared with ape branches. Six and two genes were shown to have significantly higher nonsynonymous changes at the human lineage from the acceleration index and exact tests, respectively. We also compared the distribution of the differences of the nonsynonymous substitutions on the human lineage and those on the great ape lineage. Two genes were more conserved in the ape lineage, whereas one gene was more conserved in the human lineage. These results suggest that a small proportion of protein-coding genes started to evolve differently in the human lineage after it diverged from the ape lineage.     

Estimating the distribution of selection coefficients from phylogenetic data using sitewise mutation-selection models

Estimation of the distribution of selection coefficients of mutations is a long-standing issue in molecular evolution. In addition to population-based methods, the distribution can be estimated from DNA sequence data by phylogenetic-based models. Previous models have generally found unimodal distributions where the probability mass is concentrated between mildly deleterious and nearly neutral mutations. Here we use a sitewise mutation-selection phylogenetic model to estimate the distribution of selection coefficients among novel and fixed mutations (substitutions) in a data set of 244 mammalian mitochondrial genomes and a set of 401 PB2 proteins from influenza. We find a bimodal distribution of selection coefficients for novel mutations in both the mitochondrial data set and for the influenza protein evolving in its natural reservoir, birds. Most of the mutations are strongly deleterious with the rest of the probability mass concentrated around mildly deleterious to neutral mutations. The distribution of the coefficients among substitutions is unimodal and symmetrical around nearly neutral substitutions for both data sets at adaptive equilibrium. About 0.5% of the nonsynonymous mutations and 14% of the nonsynonymous substitutions in the mitochondrial proteins are advantageous, with 0.5% and 24% observed for the influenza protein. Following a host shift of influenza from birds to humans, however, we find among novel mutations in PB2 a trimodal distribution with a small mode of advantageous mutations.     

Evolutionary divergence and salinity-mediated selection in halophilic archaea.

Halophilic (literally salt-loving) archaea are a highly evolved group of organisms that are uniquely able to survive in and exploit hypersaline environments. In this review, we examine the potential interplay between fluctuations in environmental salinity and the primary sequence and tertiary structure of halophilic proteins. The proteins of halophilic archaea are highly adapted and magnificently engineered to function in an intracellular milieu that is in ionic balance with an external environment containing between 2 and 5 M inorganic salt. To understand the nature of halophilic adaptation and to visualize this interplay, the sequences of genes encoding the L11, L1, L10, and L12 proteins of the large ribosome subunit and Mn/Fe superoxide dismutase proteins from three genera of halophilic archaea have been aligned and analyzed for the presence of synonymous and nonsynonymous nucleotide substitutions. Compared to homologous eubacterial genes, these halophilic genes exhibit an inordinately high proportion of nonsynonymous nucleotide substitutions that result in amino acid replacement in the encoded proteins. More than one-third of the replacements involve acidic amino acid residues. We suggest that fluctuations in environmental salinity provide the driving force for fixation of the excessive number of nonsynonymous substitutions. Tinkering with the number, location, and arrangement of acidic and other amino acid residues influences the fitness (i.e., hydrophobicity, surface hydration, and structural stability) of the halophilic protein. Tinkering is also evident at halophilic protein positions monomorphic or polymorphic for serine; more than one-third of these positions use both the TCN and the AGY serine codons, indicating that there have been multiple nonsynonymous substitutions at these positions. Our model suggests that fluctuating environmental salinity prevents optimization of fitness for many halophilic proteins and helps to explain the unusual evolutionary divergence of their encoding genes.     

The amino-acid mutational spectrum of human genetic disease

Background  

Nonsynonymous mutations in the coding regions of human genes are responsible for phenotypic differences between humans and for susceptibility to genetic disease. Computational methods were recently used to predict deleterious effects of nonsynonymous human mutations and polymorphisms. Here we focus on understanding the amino-acid mutation spectrum of human genetic disease. We compare the disease spectrum to the spectra of mutual amino-acid mutation frequencies, non-disease polymorphisms in human genes, and substitutions fixed between species.     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.     

Reevaluation of amino acid variability of the human immunodeficiency virus type 1 gp120 envelope glycoprotein and prediction of new discontinuous epitopes

To elucidate the evolutionary mechanisms of the human immunodeficiency virus type 1 gp120 envelope glycoprotein at the single-site level, the degree of amino acid variation and the numbers of synonymous and nonsynonymous substitutions were examined in 186 nucleotide sequences for gp120 (subtype B). Analyses of amino acid variabilities showed that the level of variability was very different from site to site in both conserved (C1 to C5) and variable (V1 to V5) regions previously assigned. To examine the relative importance of positive and negative selection for each amino acid position, the numbers of synonymous and nonsynonymous substitutions that occurred at each codon position were estimated by taking phylogenetic relationships into account. Among the 414 codon positions examined, we identified 33 positions where nonsynonymous substitutions were significantly predominant. These positions where positive selection may be operating, which we call putative positive selection (PS) sites, were found not only in the variable loops but also in the conserved regions (C1 to C4). In particular, we found seven PS sites at the surface positions of the alpha-helix (positions 335 to 347 in the C3 region) in the opposite face for CD4 binding. Furthermore, two PS sites in the C2 region and four PS sites in the C4 region were detected in the same face of the protein. The PS sites found in the C2, C3, and C4 regions were separated in the amino acid sequence but close together in the three-dimensional structure. This observation suggests the existence of discontinuous epitopes in the protein's surface including this alpha-helix, although the antigenicity of this area has not been reported yet.     

Insertional polymorphisms of full-length endogenous retroviruses in humans

Human endogenous retrovirus K (HERV-K) is distinctive among the retroviruses in the human genome in that many HERV-K proviruses were inserted into the human germline after the human and chimpanzee lineages evolutionarily diverged [1, 2]. However, all full-length endogenous retroviruses described to date in humans are sufficiently old that all humans examined were homozygous for their presence [1]. Moreover, none are intact; all have lethal mutations [1, 3, 4]. Here, we describe the first endogenous retroviruses in humans for which both the full-length provirus and the preintegration site alleles are shown to be present in the human population today. One provirus, called HERV-K113, was present in about 30% of tested individuals, while a second, called HERV-K115, was found in about 15%. HERV-K113 has full-length open reading frames (ORFs) for all viral proteins and lacks any nonsynonymous substitutions in amino acid motifs that are well conserved among retroviruses. This is the first such endogenous retrovirus identified in humans. These findings indicate that HERV-K remained capable of reinfecting humans through very recent evolutionary times and that HERV-K113 is an excellent candidate for an endogenous retrovirus that is capable of reinfecting humans today.     

Positive and negative selection on the human genome.

The distinction between deleterious, neutral, and adaptive mutations is a fundamental problem in the study of molecular evolution. Two significant quantities are the fraction of DNA variation in natural populations that is deleterious and destined to be eliminated and the fraction of fixed differences between species driven by positive Darwinian selection. We estimate these quantities using the large number of human genes for which there are polymorphism and divergence data. The fraction of amino acid mutations that is neutral is estimated to be 0.20 from the ratio of common amino acid (A) to synonymous (S) single nucleotide polymorphisms (SNPs) at frequencies of > or =15%. Among the 80% of amino acid mutations that are deleterious at least 20% of them are only slightly deleterious and often attain frequencies of 1-10%. We estimate that these slightly deleterious mutations comprise at least 3% of amino acid SNPs in the average individual or at least 300 per diploid genome. This estimate is not sensitive to human population history. The A/S ratio of fixed differences is greater than that of common SNPs and suggests that a large fraction of protein divergence is adaptive and driven by positive Darwinian selection.     

Inference of expression-dependent negative selection based on polymorphism and divergence in the human genome

There is a mounting evidence for the correlation between the gene expression pattern and sequence divergence. However, little is known about the relationship between the gene expression pattern and polymorphism. We compiled the gene expression, polymorphism, and divergence data from the public databases of the human genome. The ratios of nonsynonymous (A) to synonymous (S) substitutions in polymorphism and divergence in the human genome were strongly influenced by the expression pattern and breadth of genes and showed strong correlations. Among the tissues we analyzed, the brain-expressed genes have the smallest and the liver-expressed genes have the largest proportion of amino acid changes both in polymorphism and divergence. The analysis implies that negative selection is the primary factor affecting expression-dependent gene evolution and the prevalent but nonuniform distribution of slightly deleterious mutations in the genome. Although the genes under relaxed negative selection evolved faster than the other genes, these genes are even more liable to slightly deleterious mutations in the population. On the other hand, nonneutral mutations in the highly conservative genes, such as brain-expressed and housekeeping genes, are largely deleterious and eliminated before they enter the population.     

Mutations that cause amino acid substitutions at the invariant positions in homeodomain of OSH3 KNOX protein suggest artificial selection during rice domestication

KNOX homeodomain (HD) proteins encoded by KNOTTED1-like homeobox genes (KNOX genes) are considered to work as important regulators for plant developmental and morphogenetic events. We found that OSH3, one of the KNOX genes isolated from a cultivar of Oryza sativa (Nipponbare), encodes a novel HD, which has two amino acid substitutions at invariant positions. Sequence analysis of OSH3 from various domesticated and wild species of rice has revealed that these substitutions are distributed only in Japonica and Javanica type of O. sativa, two groups of domesticated rice in Asia. Surprisingly, nucleotide sequences in the first intron are almost conserved in the rice strains that have the substitutions at the invariant amino acids. Overexpression studies revealed that these invariant amino acids are critical for the function of OSH3 in vivo. The facts that these substitutions occurred specifically at the functionally important amino acids and the sequences are conserved in intron where neutral mutations accumulate suggest the substitutions at the invariant positions of OSH3 have been fixed by artificial selections during domestication. Based on these observations, we hypothesize that OSH3 is responsible for one of the traits that are selectively introduced during the domestication of most of Japonica and a part of Javanica type of rice.     

Genetic polymorphism and sequence evolution of an alternatively spliced exon of the glial fibrillary acidic protein gene,GFAP

Isoform GFAPepsilon of the human cytoskeletal protein GFAP carries, as the result of alternative splicing of exon 7a of GFAP, a novel 42-amino-acid-long C-terminal region with binding capacity for the presenilin proteins. Here we show that exon 7a is present in a variety of mammals but absent from GFAP of chicken and fish. Comparison of the mouse and human GFAP exons showed an increased rate of nonsynonymous nucleotide substitutions in exon 7a compared to the other exons. This resulted in 10 nonconservative and 2 conservative amino acid substitutions and suggests that exon 7a has evolved under different functional constraints. Exons 7a of humans and higher primates are 100% identical apart from alanine codon 426, which is conserved in only 9% of the human alleles, while 21 and 70% of the alleles, respectively, have a valine or a threonine codon at that position. Threonine represents a potential phosphorylation site, and positive selection of that effect could explain the high allele frequency.     

Structure, function, and evolution of the family of superoxide dismutase proteins from halophilic archaebacteria.

The protein sequences of seven members of the superoxide dismutase (SOD) family from halophilic archaebacteria have been aligned and compared with each other and with the homologous Mn and Fe SOD sequences from eubacteria and the methanogenic archaebacterium Methanobacterium thermoautotrophicum. Of 199 common residues in the SOD proteins from halophilic archaebacteria, 125 are conserved in all seven sequences, and 64 of these are encoded by single unique triplets. The 74 remaining positions exhibit a high degree of variability, and for almost half of these, the encoding triplets are connected by at least two nonsynonymous nucleotide substitutions. The majority of nucleotide substitutions within the seven genes are nonsynonymous and result in amino acid replacement in the respective protein; silent third-codon-position (synonymous) substitutions are unexpectedly rare. Halophilic SODs contain 30 specific residues that are not found at the corresponding positions of the methanogenic or eubacterial SOD proteins. Seven of these are replacements of highly conserved amino acids in eubacterial SODs that are believed to play an important role in the three-dimensional structure of the protein. Residues implicated in formation of the active site, catalysis, and metal ion binding are conserved in all Mn and Fe SODs. Molecular phylogenies based on parsimony and neighbor-joining methods coherently group the halophile sequences but surprisingly fail to distinguish between the Mn SOD of Escherichia coli and the Fe SOD of M. thermoautotrophicum as the outgroup. These comparisons indicate that as a group, the SODs of halophilic archaebacteria have many unique and characteristic features. At the same time, the patterns of nucleotide substitution and amino acid replacement indicate that these genes and the proteins that they encode continue to be subject to strong and changing selection. This selection may be related to the presence of oxygen radicals and the inter- and intracellular composition and concentration of metal cations.     

Patterns of intra- and interhost nonsynonymous variation reveal strong purifying selection in dengue virus

Considerable uncertainty surrounds the evolutionary rates of and selection pressures acting on arthropod-borne RNA viruses (arboviruses). In particular, it is unclear why arboviruses such as dengue virus show substantial genetic variation within individual humans and mosquitoes yet low long-term rates of amino acid substitution. To address this question, I compared patterns of nonsynonymous variation in populations of dengue virus sampled at different levels of evolutionary divergence. Although nonsynonymous variation was abundant in viral populations within individual humans, there was a marked reduction in the frequency of nonsynonymous mutations in interhost comparisons. Moreover, intrahost genetic variation corresponded to a random pattern of mutation, and most of the sites that exhibited nonsynonymous variation within hosts were invariant at deeper phylogenetic levels. This loss of long-term nonsynonymous variation is the signature of extensive purifying selection such that more than 90% of all nonsynonymous mutations are deleterious. Consequently, although arboviruses are able to successfully adapt to diverse cell types, they are characterized by a high rate of deleterious mutation.     

Evolution of the major histocompatibility complex: independent origin of nonclassical class I genes in different groups of mammals

The class I major histocompatibility complex genes are composed of classical and nonclassical genes, the latter being largely nonfunctional. To understand the evolutionary relationships of the two groups of class I genes, a phylogenetic analysis of DNA sequences was conducted using 45 genes from six mammalian and one avian species. The results indicate that nonclassical genes in one species are more closely related to classical genes from the same species than to nonclassical genes from a species belonging to a different order or family. This indicates that the differentiation of classical and nonclassical genes occurs rather rapidly in the genome. Classical genes are apparently duplicated with a high frequency in the evolutionary process, and many of the duplicated genes seem to degenerate into nonclassical genes as a result of deleterious mutation. The nonclassical Qa genes in the mouse have sequences homologous to regulatory sequences involved in the universal expression of classical class I genes, but they have accumulated numerous nucleotide substitutions in these sequences. The pattern of nucleotide substitution in nonclassical genes is different from that in classical genes. In nonclassical genes, the rate of nonsynonymous substitution is higher in the antigen recognition site than in other gene regions, as is true of classical genes. However, unlike the case of classical genes, the nonsynonymous rate does not always exceed the synonymous rate in the antigen recognition site. Nonclassical proteins further differ from classical proteins in having amino acid replacements in conserved antigen recognition site positions. These observations are consistent with the hypothesis that nonclassical genes have originated from classical genes but have lost classical class I function because of deleterious mutation.