Gene conversion generates hypervariability at the variable regions of kallikreins and their inhibitors.

The two mechanisms for generating hypervariability at the reactive center of serine proteases and their inhibitors are gene conversion followed by natural selection and natural selection for point mutation. One way to clarify the effects of these two mechanisms is to calculate separately the number of nonsynonymous substitutions and that of synonymous substitutions at the variable regions and at the conserved regions. Our data analysis shows that not only the number of nonsynonymous substitutions but also the number of synonymous substitutions at the variable regions exceed the corresponding numbers at the conserved regions. Thus gene conversion has provided needed variability at the variable regions of serine proteases and their inhibitors. Natural selection has helped perpetuate such variability.     

Molecular evolution of cadherin-related neuronal receptor/protocadherin(alpha) (CNR/Pcdh(alpha)) gene cluster in Mus musculus subspecies

The mouse cadherin-related neuronal receptor/protocadherin (CNR/Pcdh) gene clusters are located on chromosome 18. We sequenced single-nucleotide polymorphisms (SNPs) of the CNR/Pcdh(alpha)-coding region among 12 wild-derived and four laboratory strains; these included the four major subspecies groups of Mus musculus: domesticus, musculus, castaneus, and bactrianus. We detected 883 coding SNPs (cSNPs) in the CNR/Pcdh(alpha) variable exons and three in the constant exons. Among all the cSNPs, 586 synonymous (silent) and 297 nonsynonymous (amino acid exchanged) substitutions were found; therefore, the K(a)/K(s) ratio (nonsynonymous substitutions per synonymous substitution) was 0.51. The synonymous cSNPs were relatively concentrated in the first and fifth extracellular cadherin domain-encoding regions (ECs) of CNR/Pcdh(alpha). These regions have high nucleotide homology among the CNR/Pcdh(alpha) paralogs, suggesting that gene conversion events in synonymous and homologous regions of the CNR/Pcdh(alpha) cluster are related to the generation of cSNPs. A phylogenetic analysis revealed gene conversion events in the EC1 and EC5 regions. Assuming that the common sequences between rat and mouse are ancestral, the GC content of the third codon position has increased in the EC1 and EC5 regions, although biased substitutions from GC to AT were detected in all the codon positions. In addition, nonsynonymous substitutions were extremely high (11 of 13, K(a)/K(s) ratio 5.5) in the laboratory mouse strains. The artificial environment of laboratory mice may allow positive selection for nonsynonymous amino acid variations in CNR/Pcdh(alpha) during inbreeding. In this study, we analyzed the direction of cSNP generation, and concluded that subspecies-specific nucleotide substitutions and region-restricted gene conversion events may have contributed to the generation of genetic variations in the CNR/Pcdh genes within and between species.     

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.     

Analysis of VH gene diversity in rainbow trout (Oncorhynchus mykiss): both nonsynonymous and synonymous nucleotide changes are more frequent in CDRs than in FRs

 Forty-six immunoglobulin VH gene sequences of rainbow trout were compiled to analyze the extent of variations and the frequency of nucleotide changes in CDRs and FRs. The results show that the frequency of nonsynonymous (amino acid replacing) changes (Ka) are on average 4.9 times higher in complementarity determining regions (CDRs) than in FRs, thus contributing more diversity in CDRs. Unexpectedly, however, the frequency of synonymous (silent) changes (Ks) show the same tendency: it was 5.3 times higher in CDRs than in framework regions (FRs). The distribution of Ks/Ka ratios of each comparison shows no segregation between CDRs and FRs. The same analysis applied to five germline VH genes of Heterodontus francisci shows the same result as was found with the rainbow trout. In contrast, the results from mouse data show that, while the CDR/FR ratio for Ka is much higher (7.4), the CDR/FR ratio for Ks is only slightly higher (1.8). The distribution of Ks/Ka ratios in mouse indicates clear segregation between CDRs and FRs. This suggests that CDR germline diversity is largely generated by gene conversion in VHs of lower vertebrates such as rainbow trout or shark. This mechanism might be advantageous to lower vertebrates in generating V gene diversity faster than other mechanisms such as point mutation and selection. Received: 10 July 1996 / Revised: 9 September 1996     

Analysis of polymorphisms in the major expressed class I locus (B-FIV) of the chicken

 We analyzed the polymorphic nature of eleven alleles expressed by the major class I locus (B-FIV) in chickens. Similar to mammalian class I loci, the nucleotide substitutions with high variability occur in exons 2 and 3 encoding the α1 and α2 domains. However, the nonsynonymous to synonymous ratio of nucleotide substitutions in exon 3 encoding the α helix and β sheets is reversed compared with HLA. The region of exon 3 encoding the α2 helix demonstrates a much lower nonsynonymous to synonymous ratio, suggesting evolutionary selection of a more conserved α2 helix in B-FIV compared with HLA. Amino acid residues with high Wu-Kabat variability are typically located in positions predicted to impact antigen presentation. B-FIV amino acid residues predicted to interact with the CDR1α region of the T-cell receptor (Tcr) demonstrate less variability than in mouse and human class I alleles. The combination of a reduced nonsynonymous to synonymous ratio in exon 3 encoding the α2 helix and the limited variability in CDR1α contact residues is discussed with regard to concerted evolution between a minimal major histocompatibility complex and compaction of Tcr variable gene segments in the chicken. Received: 18 Juli 1997 /  Revised: 21 November 1997     

Evolution of Rh Blood Group Genes Have Experienced Gene Conversions and Positive Selection

There are two tightly linked loci (D and CE) for the human Rh blood group. Their gene products are membrane proteins having 12 transmembrane domains and form a complex with Rh50 glycoprotein on erythrocytes. We constructed phylogenetic networks of human and nonhuman primate Rh genes, and the network patterns suggested the occurrences of gene conversions. We therefore used a modified site-by-site reconstruction method by using two assumed gene trees and detected 9 or 11 converted regions. After eliminating the effect of gene conversions, we estimated numbers of nonsynonymous and synonymous substitutions for each branch of both trees. Whichever gene tree we selected the branch connecting hominoids and Old World monkeys showed significantly higher nonsynonymous than synonymous substitutions, an indication of positive selection. Many other branches also showed higher nonsynonymous than synonymous substitutions; this suggests that the Rh genes have experienced some kind of positive selection. Received: 16 March 1999 / Accepted: 17 June 1999     

Synonymous Substitutions in the Xdh Gene of Drosophila: Heterogeneous Distribution along the Coding Region

The Xdh (rosy) region of Drosophila subobscura has been sequenced and compared to the homologous region of D. pseudoobscura and D. melanogaster. Estimates of the numbers of synonymous substitutions per site (Ks) confirm that Xdh has a high synonymous substitution rate. The distributions of both nonsynonymous and synonymous substitutions along the coding region were found to be heterogeneous. Also, no relationship has been detected between Ks estimates and codon usage bias along the gene, in contrast with the generally observed relationship among genes. This heterogeneous distribution of synonymous substitutions along the Xdh gene, which is expression-level independent, could be explained by a differential selection pressure on synonymous sites along the coding region acting on mRNA secondary structure. The synonymous rate in the Xdh coding region is lower in the D. subobscura than in the D. pseudoobscura lineage, whereas the reverse is true for the Adh gene.     

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.      

On the Pattern of Polymorphisms at Major Histocompatibility Complex Loci

The pattern of polymorphisms at major histocompatibility complex loci was studied by computer simulations and by DNA sequence analysis. Two types of selection, overdominance plus short-term selection and maternal–fetal incompatibility, were simulated for a gene family with intra- and interlocus gene conversion. Both types of selection were found to be consistent with the observed patterns of polymorphisms. It was also found that the more interlocus conversion occurs, the higher the divergence becomes at both nonsynonymous and synonymous sites. The ratio of nonsynonymous-to-synonymous divergence among alleles decreases as the interlocus conversion rate increases. These results agree with the interpretation that the rate of interlocus conversion is lower in human genes than in genes of other nonprimate mammals. This is because, in the latter, synonymous divergence at the ARS (antigen recognition site) is often higher than that at the non-ARS, whereas in the former, this is not so. Also, the ratio of nonsynonymous to synonymous substitutions at the ARS tends to be higher in human genes than in other mammalian genes. The main difference between overdominance plus short-term selection and maternal–fetal interaction is that the number of alleles and heterozygosity per locus are higher in the latter than in the former under the presumed selection intensities. However, the average divergence among alleles tends to be lower in the latter than in the former under similar conditions. Received: 30 September 1997 / Accepted: 15 December 1997     

Potential evolutionary influences on overlapping reading frames in the bovine leukemia virus pXBL region

Bovine leukemia virus contains a pXBL region encoding the 3' parts of four regulatory proteins (Tax, Rex, G4, R3) in overlapping reading frames. Here we report the pXBL polymorphisms of 30 isolates from four countries. Rates of overall and synonymous substitutions were consistently lower, and nucleotide/amino acid composition bias and codon bias higher, in more-overlapped than in less-overlapped regions. Ratios of nonsynonymous/synonymous substitutions were lowest in the tax gene and its subregions. The 5' parts of the four genes showed selection patterns corresponding to their genomic context outside of the pXBL region. Longer G4 variants due to a natural stop codon mutation had additional triple overlap with reduced sequence variability. These data support the concept that a higher level of overlapping in coding regions correlates with greater evolutionary constraint. Tax, the most conserved among the four regulatory proteins, showed purifying selection consistent with its importance in the viral life cycle.     

A method for detecting positive selection at single amino acid sites

A method was developed for detecting the selective force at single amino acid sites given a multiple alignment of protein-coding sequences. The phylogenetic tree was reconstructed using the number of synonymous substitutions. Then, the neutrality was tested for each codon site using the numbers of synonymous and nonsynonymous changes throughout the phylogenetic tree. Computer simulation showed that this method accurately estimated the numbers of synonymous and nonsynonymous substitutions per site, as long as the substitution number on each branch was relatively small. The false-positive rate for detecting the selective force was generally low. On the other hand, the true-positive rate for detecting the selective force depended on the parameter values. Within the range of parameter values used in the simulation, the true-positive rate increased as the strength of the selective force and the total branch length (namely the total number of synonymous substitutions per site) in the phylogenetic tree increased. In particular, with the relative rate of nonsynonymous substitutions to synonymous substitutions being 5.0, most of the positively selected codon sites were correctly detected when the total branch length in the phylogenetic tree was > or = 2.5. When this method was applied to the human leukocyte antigen (HLA) gene, which included antigen recognition sites (ARSs), positive selection was detected mainly on ARSs. This finding confirmed the effectiveness of the present method with actual data. Moreover, two amino acid sites were newly identified as positively selected in non-ARSs. The three-dimensional structure of the HLA molecule indicated that these sites might be involved in antigen recognition. Positively selected amino acid sites were also identified in the envelope protein of human immunodeficiency virus and the influenza virus hemagglutinin protein. This method may be helpful for predicting functions of amino acid sites in proteins, especially in the present situation, in which sequence data are accumulating at an enormous speed.     

Genome-wide patterns of nucleotide substitution reveal stringent functional constraints on the protein sequences of thermophiles

To test the hypothesis that the proteins of thermophilic prokaryotes are subject to unusually stringent functional constraints, we estimated the numbers of synonymous and nonsynonymous nucleotide substitutions per site between 17,957 pairs of orthologous genes from 22 pairs of closely related species of Archaea and Bacteria. The average ratio of nonsynonymous to synonymous substitutions was significantly lower in thermophiles than in nonthermophiles, and this effect was observed in both Archaea and Bacteria. There was no evidence that this difference could be explained by factors such as nucleotide content bias. Rather, the results support the hypothesis that proteins of thermophiles are subject to unusually strong purifying selection, leading to a reduced overall level of amino acid evolution per mutational event. The results show that genome-wide patterns of sequence evolution can be influenced by natural selection exerted by a species' environment and shed light on a previous observation that relatively few of the mutations arising in a thermophilic archaeon were nucleotide substitutions in contrast to indels.     

Positive Darwinian selection operating on the immunoglobulin heavy chain of Antarctic fishes

The cooling of the Southern Ocean to the freezing point of seawater (-1.9 degrees C) over the past 25 million years played a dominant selective role in the evolution of the Antarctic fish fauna. During this period, the perciform suborder Notothenioidei, which is largely endemic to the Antarctic, diversified and developed numerous cold-adapted characters. In this report, we provide compelling evidence that the immunoglobulin heavy chain (IgH) of the notothenioid fishes has undergone adaptive selection. Two and four IgH clones were isolated, respectively, from spleen cDNA libraries prepared from the Antarctic icefish Chaenocephalus aceratus and the yellowbelly rockcod Notothenia coriiceps. The transmembrane region of the membrane form of the rockcod IgM heavy chain was located at the end of the second constant (C(H)) domain, in contrast to other teleost IgMs in which the transmembrane region is located at the end of the third constant domain. Phylogenetic analyses of C(H) regions revealed that rates of nonsynonymous nucleotide substitution were higher than rates of synonymous nucleotide substitution. Many of the nonsynonymous substitutions introduced charge changes, consistent with positive Darwinian selection acting to adapt the structure of the notothenioid immunoglobulins. The rates of nonsynonymous nucleotide substitutions were higher than the rates of synonymous nucleotide substitutions in complementarity determining regions of variable regions, suggesting that diversity at antigen binding sites is enhanced by genomic and/or somatic selection. Results of Southern blot hybridization experiments were consistent with a translocon type of IgH gene organization reminiscent of bony fishes and tetrapods.     

New methods for estimating the numbers of synonymous and nonsynonymous substitutions

New methods for estimating the numbers of synonymous and nonsynonymous substitutions per site were developed. The methods are unweighted pathway methods based on Kimura's two-parameter model. Computer simulations were conducted to evaluate the accuracies of the new methods, Nei and Gojobori's (NG) method, Miyata and Yasunaga's (MY) method, Li, Wu, and Luo's (LWL) method, and Pamilo, Bianchi, and Li's (PBL) method. The following results were obtained: (1) The NG, MY, and LWL methods give overestimates of the number of synonymous substitutions and underestimates of the number of nonsynonymous substitutions. The major cause for the biased estimation is that these three methods underestimate the number of synonymous sites and overestimate the number of nonsynonymous sites. (2) The PBL method gives better estimates of the numbers of synonymous and nonsynonymous substitutions than those obtained by the NG, MY, and LWL methods. (3) The new methods also give better estimates of the numbers of synonymous and nonsynonymous substitutions than those obtained by the NG, MY, and LWL methods. In addition, estimates of the numbers of synonymous and nonsynonymous sites obtained by the new methods are reasonably accurate. (4) In some cases, the new methods and the PBL method give biased estimates of substitution numbers. However, from the number of nucleotide substitutions at the third position of codons, we can examine whether estimates obtained by the new methods are good or not, whereas we cannot make an examination of estimates obtained by the PBL method. (5) When there are strong transition/transversion and nucleotide-frequency biases like mitochondrial genes, all of the above methods give biased estimates of substitution numbers. In such cases, Kondo et al.'s method is recommended to be used for estimating the number of synonymous substitutions, although their method cannot estimate the number of nonsynonymous substitutions and is time-consuming. These results, particularly result (1), call for reexaminations of some genes. This is because evolutionary pictures of genes have often been discussed on the basis of results obtained by the NG, MY, and LWL methods, which are favorable for the neutral theory of molecular evolution.     

Synonymous and nonsynonymous substitutions in mammalian genes and the nearly neutral theory

The nearly neutral theory of molecular evolution predicts larger generation-time effects for synonymous than for nonsynonymous substitutions. This prediction is tested using the sequences of 49 single-copy genes by calculating the average and variance of synonymous and nonsynonymous substitutions in mammalian star phylogenies (rodentia, artiodactyla, and primates). The average pattern of the 49 genes supports the prediction of the nearly neutral theory, with some notable exceptions.The nearly neutral theory also predicts that the variance of the evolutionary rate is larger than the value predicted by the completely neutral theory. This prediction is tested by examining the dispersion index (ratio of the variance to the mean), which is positively correlated with the average substitution number. After weighting by the lineage effects, this correlation almost disappears for nonsynonymous substitutions, but not quite so for synonymous substitutions. After weighting, the dispersion indices of both synonymous and nonsynonymous substitutions still exceed values expected under the simple Poisson process. The results indicate that both the systematic bias in evolutionary rate among the lineages and the episodic type of rate variation are contributing to the large variance. The former is more significant to synonymous substitutions than to nonsynonymous substitutions. Isochore evolution may be similar to synonymous substitutions. The rate and pattern found here are consistent with the nearly neutral theory, such that the relative contributions of drift and selection differ between the two types of substitutions. The results are also consistent with Gillespie's episodic selection theory.     

Pattern of Nucleotide Substitutions in Growth Hormone-Prolactin Gene Family: A Paradigm for Evolution by Gene Duplication

The growth hormone-prolactin gene family in mammals is an interesting example of evolution by gene duplication. Divergence among members of duplicated gene families and among species was examined by using reported gene sequences of growth hormone, prolactin and their receptors. Sequence divergence among species was found to show a general tendency in which a generation-time effect is pronounced for synonymous substitutions but not so for nonsynonymous substitutions. Divergence among duplicated genes is characterized by the relatively high rate of nonsynonymous substitutions, i.e., the rate is close to that of synonymous ones. In view of the stage- and tissue-specific expression of duplicated genes, some of the amino acid substitutions among duplicated genes is likely to be caused by positive Darwinian selection.     

Nucleotide variability at the acetyl coenzyme A carboxylase gene and the signature of herbicide selection in the grass weed Alopecurus myosuroides (Huds.)

Acetyl coenzyme A carboxylase (ACCase) is the target of highly effective herbicides. We investigated the nucleotide variability of the ACCase gene in a sample of 18 black-grass (Alopecurus myosuroides [Huds.]) populations to search for the signature of herbicide selection. Sequencing 3,396 bp encompassing ACCase herbicide-binding domain in 86 individuals revealed 92 polymorphisms, which formed 72 haplotypes. The ratio of nonsynonymous versus synonymous substitutions was very low, in agreement with ACCase being a vital metabolic enzyme. Within black grass, most nonsynonymous substitutions were related to resistance to ACCase-inhibiting herbicides. Differentiation between populations was strong, in contrast to expectations for an allogamous, annual plant. Significant H tests revealed recent hitchhiking events within populations. These results were consistent with recent and local positive selection. We propose that, although they have only been used since at most 15 black-grass generations, ACCase-inhibiting herbicides have exerted a positive selection targeting resistant haplotypes that has been strong enough to have a marked effect upon ACCase nucleotide diversity. A minimum-spanning network of nonrecombinant haplotypes revealed multiple, independent apparitions of resistance-associated mutations. This study provides the first evidence for the signature of ongoing, recent, pesticide selection upon variation at the gene encoding the targeted enzyme in natural plant populations.